In fluid dynamics, the generation and bursting of surface bubble liquid films in surfactant-laden environments involve complex phenomena, one of which is marginal pinching at the film's foot—a crucial yet inadequately understood aspect due to experimental limitations. Considering its profound impact on liquid film drainage and lifetime, we utilize high-resolution numerical simulations incorporating a weakly compressible scheme and adaptive mesh refinement to dissect the marginal pinching dynamics with unprecedented detail. Our approach elucidates the pinching dynamics and tracks the evolution of film thickness during the critical late-stage drainage process. By leveraging detailed geometric data from pinched regions, we significantly refine existing drainage models, enhancing their predictive accuracy regarding rupture thickness and film lifetime across various viscosities, surfactant concentrations, and bubble sizes. This refined model demonstrates robust alignment with our simulation results. Furthermore, we establish a quantifiable relationship between the prefactor governing reverse flow induced by the Marangoni effect and surfactant concentration. The methodologies and findings of this study provide foundational knowledge that paves the way for optimized industrial processes and an enhanced understanding of natural phenomena.