We investigate the statistical properties of kinetic ϵu and thermal ϵθ energy dissipation rates in two-dimensional (2D) thermal vibrational convection (TVC). Direct numerical simulations were conducted in a unit aspect ratio box across a dimensionless angular frequency range of 103≤ω≤107 for amplitudes 0.001≤a≤0.1, with a fixed Prandtl number of 4.38. Our findings indicate ϵu is primarily associated with the characteristics of the vibration force, while ϵθ is more related to the large-scale columnar structures. Both energy dissipation rates exhibit a power-law relationship with the oscillational Reynolds number Reos. ϵu exhibits a scaling relation as ⟨ϵu⟩V,t∼a−1Reos0.93±0.01, while ϵθ exhibits two distinct scaling behaviors, i.e., ⟨ϵθ⟩V,t∼a−1Reos1.97±0.04 for Reos<Reos,cr and ⟨ϵθ⟩V,t∼a−1Reos1.31±0.02 for Reos>Reos,cr, where the fitted critical oscillational Reynolds number is approximately Reos,cr≈80. The different scaling of ⟨ϵθ⟩V,t is determined by the competition between the thermal boundary layer and the oscillating boundary layer. Moreover, the probability density functions (PDFs) of both dissipation rates deviate significantly from the lognormal distribution and exhibit a bimodal shape. By partitioning the contributions from the boundary layer and bulk regions, it is shown that the bulk contributes to the small and moderate dissipation rates, whereas the high dissipation rates are predominantly contributed by the boundary layer. As Reos increases, the heavy tail of the PDFs becomes more pronounced, revealing an enhanced level of small-scale intermittency. This small-scale intermittency is mainly caused by the influence of BL due to vibration. Our study provides insight into the small-scale characteristics of 2D TVC, highlighting the non-trivial scaling laws and intermittent behavior of energy dissipation rates with respect to vibration intensity.