In this work, we develop a mesoscopic lattice Boltzmann Bhatnagar-Gross-Krook (BGK) model to solve (2 + 1)-dimensional wave equation with the nonlinear damping and source terms. Through the Chapman-Enskog multiscale expansion, the macroscopic governing evolution equation can be obtained accurately by choosing appropriate local equilibrium distribution functions. We validate the present mesoscopic model by some related issues where the exact solution is known. It turned out that the numerical solution is in very good agreement with exact one, which shows that the present mesoscopic model is pretty valid, and can be used to solve more similar nonlinear wave equations with nonlinear damping and source terms, and predict and enrich the internal mechanism of nonlinearity and complexity in nonlinear dynamic phenomenon.