The anode porous transport layer (PTL) plays a crucial role in improving the performance of proton exchange membrane (PEM) electrolyzer. However, oxygen accumulation may cover the interface between the anode catalytic layer and PTL, decreasing mass transfer and electrolysis efficiency. In this study, we employed a comprehensive simulation approach, combining lattice Boltzmann simulation with the Immersed Boundary Method (IBM) and Phase Field Model (PFM), to investigate oxygen removal from the anode PTL. The simulation can accurately capture bubble morphology and gas saturation, demonstrating the four stages of gas removal: invasion, exhaust, breakup and retraction. Notably, the invasion stage contributes to approximately 3/4 of the total gas removal time, whereas the invasion time is significantly influenced by the critical throat in PTL. Removing the critical throat may lead to a 20% reduction in the invasion time. To determine the critical throat, we further simulate the breakthrough of a single pore by a bubble. At lower current densities (1.4 A/cm2), reducing interfacial tension facilitates bubble penetration through smaller pores but weakens the longitudinal permeability, leading to a longer invasion stage. However, at higher current densities (14 A/cm2), the breakthrough is predominantly governed by the critical throat. This suggests that optimizing the PTL structure is more effective in mitigating bubble coverage than adjusting interfacial tension under high current density conditions. These findings may provide some guidelines for enhancing the performance of PEM electrolyzer under high current density conditions.