In this study, we take into account the (2 + 1)-dimensional Boussinesq equation, a nonlinear evolution partial differential equation that describes how gravity waves move across the surface of the ocean. The symmetry reductions and group invariant precise solutions are systematically determined using the Lie symmetry analysis. We derive the precise multiple wave solutions using the multiple exp-function method, and then, using the multiplier method, we give the conservation laws. The dynamics of complicated waves and their interplay are faithfully recreated by the findings.