The present research focuses on the mechanical properties and stress evolution of gas-bearing soft coal seams during drilling, which are affected by a multitude of complex factors such as high ground stress, gas pressure, and pre-existing fractures. In this study, a combination of PFC2D (Particle Flow Code in 2 Dimensions) numerical simulation and theoretical analysis is employed to investigate the borehole mechanics and fracture evolution characteristics under diverse complex conditions and to determine the factors influencing different forms of borehole failure in soft coal seams. The principal outcomes are as follows: (1) At a horizontal displacement of 0.1 m from the borehole orifice of the soft coal seam, a stress peak value of 13.9 MPa is attained; the peak value of the coal body contact force is 15.8 MPa; the peak value of the displacement is 0.008 m; and the porosity of the coal body around the borehole ranges from 0.14 to 0.35. (2) With an increase in the number of pre-existing fractures, the inclination progressively aligns with that of the pre-existing fractures. Maximum values of contact force (5.13–51.9 MPa), stress (3.19–37.2 MPa), shape dimension, and fracture angle (140–150°) are achieved under the highest lateral pressure coefficient and gas pressure (1.5 MPa). (3) The borehole energy is directly proportional to the number of pre-existing fractures, the lateral pressure coefficient, and gas pressure. The number of pre-existing fractures has the most significant impact on the damage degree, followed by the lateral pressure coefficient and then the gas pressure. (4) Two types of failure are identified: fracture-dominated failure, which is controlled by the geometric distribution of pre-existing fractures, and stress-dominated failure, wherein the failure zone gradually extends both upward and downward with an increasing lateral pressure coefficient.

This study addresses the abnormal emission of pressure-relieved methane under high-intensity mining conditions by integrating geostatistical inversion, FLAC^3D-COMSOL coupled numerical simulations, and stable carbon–hydrogen isotopic tracing. Focusing on the 12023 working face at Wangxingzhuang Coal Mine, we established a heterogeneous methane reservoir model to analyze the mechanical responses of surrounding rock, permeability evolution, and gas migration patterns under mining intensities of 2–6 m/d. Key findings include the following: (1) When the working face advanced 180 m, vertical stress in concentration zones increased significantly with mining intensity, peaking at 12.89% higher under 6 m/d compared to 2 m/d. (2) Higher mining intensities exacerbated plastic failure in floor strata, with a maximum depth of 47.9 m at 6 m/d, enhancing permeability to 223 times the original coal seam. (3) Isotopic fingerprinting and multi-method validation identified adjacent seams as the dominant gas source, contributing 77.88% of total emissions. (4) Implementing targeted long directional drainage boreholes in floor strata achieved pressure-relief gas extraction efficiencies of 34.80–40.95%, reducing ventilation air methane by ≥61.79% and maintaining return airflow methane concentration below 0.45%. This research provides theoretical and technical foundations for adaptive gas control in rapidly advancing faces through stress–permeability coupling optimization, enabling the efficient interception and resource utilization of pressure-relieved methane. The outcomes support safe, sustainable coal mining practices and advance China’s Carbon Peak and Neutrality goals.