It is well known that the collapse of heterogeneous multi-cavity near the wall will induce the fluctuation of the load field. To address this problem, the Lattice Boltzmann Method (LBM) is applied to model the three-phase coupling between gas-liquid-solid. The objective is to investigate the evolution of heterogeneous double bubbles and the spatial-temporal distribution characteristics of wall loads induced near the wall. In this study, the pseudopotential Multi-Relaxation-Time Lattice Boltzmann Model (MRT-LBM) and the Carnahan-Starling Equation of State (C-S-EOS) with an extended format for the external force term are used. The effects of the distance of the bubble to the wall, the pressure differences between the inside and outside of the bubble, and the relative size of the bubble on the dynamic evolution and the load distribution characteristics of heterogeneous multi-bubbles near the wall are investigated in order to determine the influence of these factors. Under a two-dimensional pressure field, the collapse process of double cavitation bubbles is visualized. Through the flow field, the morphological changes of the cavitation bubble collapse near the wall are also described. Various parameters are found to have an influence on the evolution of double cavitation bubbles near the wall and the resulting load field. The study employs the Lattice Boltzmann Method and the Potential Model for the analysis of the heterogeneous bubble collapses in the near wall region.