This study introduces the distributed-order fractional version of the nonlinear two-dimensional Sobolev equation. The orthonormal Chebyshev cardinal polynomials are used to construct a numerical method for this equation. To this end, some derivative matrices related to these polynomials are obtained. The proposed approach turns to solve this equation into solving a nonlinear system of algebraic equations by approximating the unknown solution using the expressed polynomials and employing their derivative matrices. The applicability and validity of this method are examined by solving three examples.