In this paper, a three-dimensional numerical simulation of the motion behavior of bubbles in complex porous medium channels in a large density ratio gas-liquid system is conducted based on the lattice Boltzmann method. The Eötvös number ( Eo), contact angle ( θ) and Reynolds number ( Re) are systematically discussed with emphasis on the law of their coupling effect affecting bubble velocity, morphological evolution and stagnation phenomenon. The results show that the increase of contact angle will reduce the bubble velocity but intensify the velocity fluctuations, making the bubbles tend flat, while the increase of Eonumber significantly suppresses the influence of the contact angle, stabilizes the bubble velocity, and makes its shape close to a bullet head shape. When the contact angle is large ( θ> 90°) and the Eonumber is small ( Eo< 10), the adhesion force is significantly enhanced and the bubbles will stagnate inside the porous medium. Renumber and contact angle compete in the generation of resistance, and have mutually reinforcing effects on the average velocity of bubbles and interface evolution. The larger contact angle makes the deformation of the bubble tail intensify and becomes unstable, and as the Renumber further increases, the tail tentacles are more likely to break, forming residual bubbles. It is also found in this work that the coupling between Eonumber and Renumber significantly affects bubble behavior in motion and morphological evolution. Under the conditions of high Eonumber ( Eo≥ 25) and high Renumber ( Re≥14), the bubble velocity increases with the Eonumber rising, and the trend becomes more significant as the Renumber increases; while under the conditions of low Eonumber ( Eo< 25) and low Renumber ( Re< 14), the speed change pattern is completely opposite. This phenomenon is due to the high instability of bubble morphology under the conditions of high Eonumber and high Renumber, which affects the buoyancy and speed performance. The research results provide important guidance for optimizing the flow behavior of bubbles in porous medium.