Horizontal wells within the roof are an effective method to develop gas in broken and soft coal seams, and layer-penetrating fracturing is a key engineering method for the stimulating of horizontal wells within the roof of a coal seam. To understand the propagation law of fracture in the composite roof of coal seams, this study conducted research using numerical simulation and physical similarity simulation methods. Furthermore, engineering experiments were carried out at the Panxie coal mine in the Huainan Mining Area and the Luling coal mine in Huaibei Mining Area, to further validate this technology. The numerical simulation results indicated that fracture within the coal seam roof can propagate from the roof to the target coal seam, effectively fracturing the coal seam. Due to the coal seam’s plasticity being greater than that of the roof mudstone, the coal seam forms a broader fracture than the roof. With the increase in pseudo roof mudstone thickness and being under constant fracturing displacement, the energy consumed by the pseudo roof mudstone during fracturing causes a decrease in pore pressure when fracture extends to the coal seam, resulting in a reduction in fracture width. Therefore, the pseudo roof mudstone is an adverse factor for the expansion of hydraulic fracturing. Physical similarity simulation results demonstrated that when horizontal boreholes were arranged within the siltstone of the coal seam roof, were under reasonable vertical distance and high flow rate fracturing via fluid injection conditions, and if the coal seam had a thin pseudo roof mudstone, the fracture could propagate through the direct roof-pseudo roof interface and the pseudo roof-coal seam interface, extending to the lower coal seam. The fracture form was curved and had irregular vertical fractures, indicating that hydraulic fracturing can achieve production enhancement and the transformation of soft and hard coal seams. However, when the coal seam had a thick pseudo roof mudstone, the mudstone posed strong resistance to hydraulic fracturing, making it difficult for the fracture to propagate to the lower coal seam. Therefore, the pseudo roof mudstone plays a detrimental role in hydraulic fracturing and the production enhancement of coal seams. The engineering verification conducted at Panxie coal mine and Luling coal mine showed that by utilizing a construction drainage rate of 7.5 cubic meters per minute at Panxie coal mine, the maximum fracture length reached 218.3 m, with a maximum fracture height of 36.8 m. The maximum daily gas production of a single well reached 1450 cubic meters per day, with a total gas extraction volume of 43.62 × 10⁴ cubic meters across 671 days. At Luling coal mine, utilizing a construction drainage rate of 10 cubic meters per minute, the maximum fracture length reached 169.1 m, with a maximum fracture height of 20.5 m. The maximum daily gas production of a single well reached 10,775 cubic meters per day, with a total gas extraction volume of 590 × 10⁴ cubic meters for 1090 days. This indicated that the fracture within the roof of coal seams can penetrate the composite roof of coal seams and extend to the interior of the coal seams, achieving the purpose of transforming fractured and low-permeability coal seams and providing an effective mode of gas extraction.

Identifying hazardous karst collapse pillars (KCPs) is critical for ensuring safe coal mining operations. While previous studies have focused primarily on physical detection, the spatial clustering characteristics of KCPs have often been overlooked. This study proposes a spatial hotspot identification method based on Moran’s index and applies it to the Wangpo Coal Mine in Shanxi, China. The method integrates morphological feature analysis of KCPs with a combined weighting scheme using the Analytic Hierarchy Process (AHP) and Entropy Weight Method (EWM). A spatial distribution index (SDI) was constructed through geographic information system (GIS) overlay analysis and spatial coordinate calibration. Global Moran’s I (0.1110, p < 0.05) indicates a statistically significant positive spatial autocorrelation of KCP distribution. Local Moran’s I further reveals 11 spatially significant KCPs, including 5 high-high clusters. Geological interpretation shows that these high-risk KCPs are predominantly located near the intersections of faults and folds, highlighting the structural control on KCP formation. The proposed method provides a quantitative and spatially interpretable approach for KCP risk identification and has potential for application to other geohazards exhibiting spatial aggregation patterns.

Studying the identification of hazardous karst collapse pillars(KCPs) is vital for ensuring the safe mining of coal resources. However, current study on identifying KCPs primarily emphasizes physical detection, overlooking the spatial aggregation patterns. In this study, we proposed a hotspot identification method for KCPs using Moran's index. and carry out experiments in Wangpo Coal Mine of Shanxi, China. The method involves evaluating the KCPs by considering their morphological characteristics and using a combination of the Analytic Hierarchy Process (AHP) and Entropy Weighting Method (EWM) for quantitative assessment. Then the spatial distribution index of the KCPs(SDI) is determined through Geographic Information System(GIS) overlay analysis and coordinate calibration. The hotspots analysis resulted in a global Moran's index value of 0.1110, indicating a positive spatial correlation of the SDI in the study area. Local Moran's index is further used to identify hotspots of KCPs. A total of 11 special KCPs were identified within the study area, including 5 high-high cluster KCPs. Finally, we analysed the geological impact of fault and fold distributions on KCPs development in these high-risk areas. The results reveal that the development patterns of KCPs at these specific points are analyzed in relation to the distribution of geological structures, especially at the intersections of faults and folds, where high-risk KCPs are most likely to develop.