Developing piezoelectric-based plate-like metamaterials necessitates an effective modeling method to elucidate the omnidirectional wave properties of piezoelectric coupled inclusions on a thin plate. The commonly used methods, such as the transfer-matrix method and finite element method, are inadequate for analyzing the transmission and reflection of omnidirectional waves in a two-dimensional elastic medium. Unlike the existing methods, the multiple scattering method employs Bessel functions as the displacement-basis methods which can accurately describe the propagation and scattering characteristics of omnidirectional waves. This paper develops a novel T-matrix formulation to support the multiple scattering method, representing the input-output relationship between incident and reflected waves from a piezoelectric shunt inclusion on a host thin plate. The piezoelectric shunt inclusion comprises a varying-thickness substrate bonded with piezoelectric shunting patches. The T-matrix of the piezoelectric shunt inclusion is formulated by integrating the wave-based method with Rayleigh-Ritz method. The derived T-matrix is then used to semi-analytically analyze the far-field scattering and reflection properties of a single inclusion. Numerical results capture the scattering properties resulting from trapped mode resonances of the piezoelectric shunt inclusion. Additionally, the capability of the piezoelectric shunt damping to design and tune multiple critical coupling conditions for axisymmetric modes of thin plates is parametrically investigated by varying the values of inductors and resistors.