Time-dependent problems modeled by hyperbolic partial differential equations can be reformulated in terms of boundary integral equations and solved via the boundary element method. In this context, the analysis of damping phenomena that occur in many physics and engineering problems is a novelty. Starting from a recently developed energetic space-time weak formulation for 1D damped wave propagation problems rewritten in terms of boundary integral equations, we develop here an extension of the so-called energetic boundary element method for the 2D case. Several numerical benchmarks, whose numerical results confirm accuracy and stability of the proposed technique, already proved for the numerical treatment of undamped wave propagation problems in several dimensions and for the 1D damped case, are illustrated and discussed.