Influences of clean fracturing fluid viscosity and horizontal in-situ stress difference on hydraulic fracture propagation and morphology in coal seam

The viscosity of fracturing fluid and in-situ stress difference are the two important factors that affect the hydraulic fracturing pressure and propagation morphology. In this study, raw coal was used to prepare coal samples for experiments, and clean fracturing fluid samples were prepared using CTAB surfactant. A series of hydraulic fracturing tests were conducted with an in-house developed triaxial hydraulic fracturing simulator and the fracturing process was monitored with an acoustic emission instrument to analyze the influences of fracturing fluid viscosity and horizontal in-situ stress difference on coal fracture propagation. The results show that the number of branched fractures decreased, the fracture pattern became simpler, the fractures width increased obviously, and the distribution of AE event points was concentrated with the increase of the fracturing fluid viscosity or the horizontal in-situ stress difference. The acoustic emission energy decreases with the increase of fracturing fluid viscosity and increases with the increase of horizontal in situ stress difference. The low viscosity clean fracturing fluid has strong elasticity and is easy to be compressed into the tip of fractures, resulting in complex fractures. The high viscosity clean fracturing fluids are the opposite. Our experimental results provide a reference and scientific basis for the design and optimization of field hydraulic fracturing parameters.

Implications for fault reactivation and seismicity induced by hydraulic fracturing

Evaluating the physical mechanisms that link hydraulic fracturing(HF) operations to induced earthquakes and the anticipated form of the resulting events is significant in informing subsurface fluid injection operations. Current understanding supports the overriding role of the effective stress magnitude in triggering earthquakes, while the impact of change rate of effective stress has not been systematically addressed. In this work, a modified critical stiffness was brought up to investigate the likelihood, impact,and mitigation of induced seismicity during and after hydraulic fracturing by developing a poroelastic model based on rate-and-state fraction law and linear stability analysis. In the new criterion, the change rate of effective stress was considered a key variable to explore the evolution of this criterion and hence the likelihood of instability slip of fault. A coupled fluid flow-deformation model was used to represent the entire hydraulic fracturing process in COMSOL Multiphysics. The possibility of triggering an earthquake throughout the entire hydraulic fracturing process, from fracturing to cessation, was investigated considering different fault locations, orientations, and positions along the fault. The competition between the effects of the magnitude and change rate of effective stress was notable at each fracturing stage. The effective stress magnitude is a significant controlling factor during fracturing events, with the change rate dominating when fracturing is suddenly started or stopped. Instability dominates when the magnitude of the effective stress increases(constant injection at each fracturing stage) and the change rate of effective stress decreases(the injection process is suddenly stopped). Fracturing with a high injection rate, a fault adjacent to the hydraulic fracturing location and the position of the junction between the reservoir and fault are important to reduce the Coulomb failure stress(CFS) and enhance the critical stiffness as the significant disturbance of stresses at these positions in the coupled process. Therefore,notable attention should be given to the injection rate during fracturing, fault position, and position along faults as important considerations to help reduce the potential for induced seismicity. Our model was verified and confirmed using the case of the Longmaxi Formation in the Sichuan Basin, China, in which the reported microseismic data were correlated with high critical stiffness values. This work supplies new thoughts of the seismic risk associated with HF engineering.

Effect of hydraulic fracturing induced stress field on weak surface activation during unconventional reservoir development

Unconventional reservoirs usually contain many weak surfaces such as faults,laminae and natural fractures,and effective activation and utilization of these weak surfaces in reservoirs can significantly improve the extraction effect.In hydraulic fracturing,when the artificial fracture approaches the natural fracture,the natural fracture would be influenced by both the original in-situ stress field and the hydraulic fracturing-induced stress field.In this paper,the hydraulic fracturing-induced stress field is calculated based on the relative position of hydraulic fracture and natural fracture,the original in-situ stress,the net pressure inside the hydraulic fracture and the pore pressure of the formation.Furthermore,the stability model of the natural fracture is established by combining the Mohr-Coulomb rupture criterion,and extensive parametric studies are conducted to explore the impact of each parameter on the stability of the natural fracture.The validity of the proposed model is verified by comparing with the reservoir characteristics and fracturing process of the X-well 150e155 formation in the Songliao Basin.It is found that the stress field induced by the hydraulic fracture inhibits the activation of the natural fracture after the artificial fracture crossed the natural fracture.Therefore,for similar reservoirs as X-well 150e155,it is suggested to connect natural fractures with hydraulic fractures first and then activate natural fractures which can effectively utilize the natural fractures and form a complex fracture network.

Numerical investigation on the effect of depletion-induced stress reorientation on infill well hydraulic fracture propagation

Depletion-induced stress change causes the redistribution of stress field in reservoirs,which can lead to the reorientation of principal stresses.Stress reorientation has a direct impact on fracture propagation of infill wells.To understand the effect of stress reorientation on the propagation of infill well’s fractures,an integrated simulation workflow that combines the reservoir flow calculation and the infill well hydraulic fracturing modeling is adopted.The reservoir simulation is computed to examine the relationship between the extent of stress reversal region and reservoir properties.Then,the hydraulic fracturing model considering the altered stress field for production is built to characterize the stress evolution of secondary fracturing.Numerical simulations show that stress reorientation may occur due to the decreasing of the horizontal stresses in an elliptical region around the parent well.Also,the initial stress difference is the driving factor for stress reorientation.However,the bottom hole pressure,permeability and other properties connected with fluid flow control timing of the stress reorientation.The decrease of the horizontal stresses around the parent well lead to asymmetrical propagation of a hydraulic fracture of the infill well.The study provides insights on understanding the influence of stress reorientation to the infill well fracturing treatment and interference between parent and infill wells.

Three-dimensional distinct element modeling of fault reactivation and induced seismicity due to hydraulic fracturing injection and backflow

This paper presents a three-dimensional fully hydro-mechanical coupled distinct element study on fault reactivation and induced seismicity due to hydraulic fracturing injection and subsequent backflow process,based on the geological data in Horn River Basin,Northeast British Columbia,Canada.The modeling results indicate that the maximum magnitude of seismic events appears at the fracturing stage.The increment of fluid volume in the fault determines the cumulative moment and maximum fault slippage,both of which are essentially proportional to the fluid volume.After backflow starts,the fluid near the joint intersection keeps flowing into the critically stressed fault,rather than backflows to the wellbore.Although fault slippage is affected by the changes of both pore pressure and ambient rock stress,their contributions are different at fracturing and backflow stages.At fracturing stage,pore pressure change shows a dominant effect on induced fault slippage.While at backflow stage,because the fault plane is under a critical stress state,any minor disturbance would trigger a fault slippage.The energy analysis indicates that aseismic deformation takes up a majority of the total deformation energy during hydraulic fracturing.A common regularity is found in both fracturing-and backflow-induced seismicity that the cumulative moment and maximum fault slippage are nearly proportional to the injected fluid volume.This study shows some novel insights into interpreting fracturing-and backflowinduced seismicity,and provides useful information for controlling and mitigating seismic hazards due to hydraulic fracturing.

Mutual impact of true triaxial stress, borehole orientation and bedding inclination on laboratory hydraulic fracturing of Lushan shale

Unconventional resources like shale gas has been the focus of intense research and development for two decades. Apart from intrinsic geologic factors that control the gas shale productivity (e.g. organic matter content, bedding planes, natural fractures, porosity and stress regime among others), external factors like wellbore orientation and stimulation design play a role. In this study, we present a series of true triaxial hydraulic fracturing experiments conducted on Lushan shale to investigate the interplay of internal factors (bedding, natural fractures and in situ stress) and external factors (wellbore orientation) on the growth process of fracture networks in cubic specimens of 200 mm in length. We observe relatively low breakdown pressure and fracture propagation pressure as the wellbore orientation and/or the maximum in situ stress is subparallel to the shale bedding plane. The wellbore orientation has a more prominent effect on the breakdown pressure, but its effect is tapered with increasing angle of bedding inclination. The shale breakdown is followed by an abrupt response in sample displacement, which reflects the stimulated fracture volume. Based on fluid tracer analysis, the morphology of hydraulic fractures (HF) is divided into four categories. Among the categories, activation of bedding planes (bedding failure, BF) and natural fractures (NF) significantly increase bifurcation and fractured areas. Under the same stress regime, a horizontal wellbore is more favorable to enhance the complexity of hydraulic fracture networks. This is attributed to the relatively large surface area in contact with the bedding plane for the horizontal borehole compared to the case with a vertical wellbore. These findings provide important references for hydraulic fracturing design in shale reservoirs.

Stress tensor determination by modified hydraulic tests on pre-existing fractures:Method and stress constraints

The hydraulic testing of pre-existing fractures(HTPF)is one of the most promising in situ stress measurement methods,particularly for three-dimensional stress tensor determination.However,the stress tensor determination based on the HTPF method requires at least six tests or a minimum of 14-15 tests(under different conditions)for reliable results.In this study,we modified the HTPF method by considering the shear stress on each pre-existing fracture,which increased the number of equations for the stress tensor determination and decreased the number of tests required.Different shear stresses were attributed to different fractures by random sampling;therefore,the stress tensors were obtained by searching for the optimal solution using the least squares criterion based on the Monte Carlo method.Thereafter,we constrained the stress tensor based on the tensile strength criterion,compressive strength criterion,and vertical stress constraints.The inverted stress tensors were presented and analyzed based on the tensorial nature of the stress using the Euclidean mean stress tensor.Two stress-measurement campaigns in Weifang(Shandong Province,China)and Mercantour road tunnel(France)were implemented to highlight the validity and efficiency of the modified HTPF(M-HTPF)method.The results showed that the M-HTPF method can be applied for stress tensor inversion using only three to four tests on pre-existing fractures,neglecting the stress gradient.The inversion results were confined to relatively small distribution dispersions and were significantly reliable and stable due to the shear stresses on the fractures and the stress constraints employed.The M-HTPF method is highly feasible and efficient for complete stress tensor determination in a single borehole.

New development of hydraulic fracturing technique for in-situ stress measurement at great depth of mines

In-situ stress measurement using the hydraulic fracturing technique was made at Wanfu Coal Mine in Shandong Province, China. To solve problems caused by great measuring depth and extra thick overburden soil layers in the mine, a series of improved techniques were developed for the traditional hydraulic fracturing technique and equipment to increase their pressure-enduring ability and to ensure safe and flexible removal of the sealing packers with other experimental apparatus. Successful in-situ stress measurement at 37 points within 7 boreholes, which were mostly over 1000 m deep, was completed. Through the measurement, detailed information of in-situ stress state has been provided for mining design of the mine. The improved hydraulic fracturing technique and equipment also provide reliable tools for in-situ stress measurement at great depth of other mines.

Effect of ground stress on hydraulic fracturing of methane well

Most of the coal reservoirs in China are of low-permeability, so hydraulic fracturing is widely used to improve the per- meability in the extraction of gas by ground drilling. The ground stress around the well was analyzed by using theory of elasticity. The pressure when the well fractured is formulated and the effect of ground stress on pressure is discussed. The effect of ground-stress-differences on hydraulic fracturing was analyzed by using the numerical software RFPA2D-Flow in reference to the tectonic stress in Jincheng coal area. The results show that: 1) the position where initial fracture appears is random and fracture branches emerge when the fractures expand if ground stresses in any two directions within a horizontal plane are equal; 2) other- wise, the fractures expand in general along the direction of maximum ground stress and the critical pressure decreases with in- creasing ground-stress-differences and 3) the preferred well-disposition pattern is diamond shaped. The preferred well spacing is 250 m×300 m. This study can provide a reference for the design of wells.

Multistage hydraulic fracturing of a horizontal well for hard roof related coal burst control:Insights from numerical modelling to field application

Multistage hydraulic fracturing of horizontal wells(MFHW)is a promising technology for controlling coal burst caused by thick and hard roofs in China.However,challenges remain regarding the MFHW control mechanism of coal burst and assessment of the associated fracturing effects.In this study,these challenges were investigated through numerical modelling and field applications,based on the actual operating parameters of MFHW for hard roofs in a Chinese coal mine.A damage parameter(D)is proposed to assess the degree of hydraulic fracturing in the roof.The mechanisms and effects of MFHW for controlling coal burst are analyzed using microseismic(MS)data and front-abutment stress distribution.Results show that the degree of fracturing can be categorized into lightly-fractured(D≤0.3),moderately fractured(0.3<D≤0.6),well-fractured(0.6<D≤0.9),and over-fractured(0.9<D≤0.95).A response stage in the fracturing process,characterized by a slowdown in crack development,indicates the transition to a wellfractured condition.After MFHW,the zone range and peak value of the front-abutment stress decrease.Additionally,MS events shift from near the coal seam to the fractured roof layers,with the number of MS events increases while the average MS energy decreases.The MFHW control mechanisms of coal bursts involve mitigating mining-induced stress and reducing seismic activity during longwall retreat,ensuring stresses remain below the ultimate stress level.These findings provide a reference for evaluating MFHW fracturing effects and controlling coal burst disasters in engineering.

The Hydraulic Fracturing Optimization for Stacked Tight Gas Reservoirs Using Multilayers and Multiwells Fracturing Strategies

Based on a geology-engineering sweet spot evaluation,the high-quality reservoir zones and horizontal well landing points were determined.Subsequently,fracture propagation and production were simulated with a multilayer fracturing scenario.The optimal hydraulic fracturing strategy for themultilayer fracturing networkwas determined by introducing a vertical asymmetry factor.This strategy aimed to minimize stress shadowing effects in the vertical direction while maximizing the stimulated reservoir volume(SRV).The study found that the small vertical layer spacing of high-quality reservoirs and the presence of stress-masking layers(with a stress difference of approximately 3∼8 MPa)indicate that interlayer stress interference from multilayers and multiwells fracturing between neighboring developed formations could affect the longitudinal propagation of the reservoirs.In addition,this study investigated well spacing optimization by comparing uniformly spaced wells(100–300 m)with asymmetric spaced wells(200 m upper layer,250 m lower layer).Numerical results indicated that asymmetric spaced well placement yielded the largest stimulated reservoir volume(SRV)of 73,082 m^(3),representing a 65.42%increase compared to 100 m spaced wells.Furthermore,four different hydraulic fracturing sequences(interlayer,up-down,down-up,and center-edge)were compared for multilayer and multiwell networks.The center-edge sequence exhibited the lowest vertical asymmetry factor(0.71)and the least stress shadowing effects compared to the other sequences(0.78 for interlayer,0.75 for up-down,and 0.76 for down-up).This sequence also achieved the largest SRV(486,194m^(3)),representing an 11.87%increase compared to the down-up sequence.Therefore,the center-edge fracturing sequence is recommended formultilayer development in this block.These results offer valuable insights for optimizing well placement and fracturing sequence design in multilayer well networks,ultimately advancing the development of multilayer fracturing technology in the region.

Numerical study of hydraulic fracturing in the sectorial well-factory considering well interference and stress shadowing

In the Changqing Oilfield in northwest China, when traditional petroleum exploitation encounters forestry reserves or water source protection areas, sectorial well-factory design is proposed. The most distinct feature of a sectorial well-factory is the deviation of the well from the minimum horizontal principal stress, resulting in hydraulic fracture deflection after the initiation, along with possible well interference (i.e., fracture hit) and fracture coalescence in the oblique wells. Four indexes describing well deflection are then proposed according to fracture morphology. Several fracturing designs, including stage arrangement, fracturing sequences, and fracturing techniques are applied to study the feasibility of the sectorial well-factory design. The results show that the “gradual” or “sparse” stage arrangement, large injection rate, and simultaneous multifracture treatment can help to optimize the fracture morphology and stimulation design. However, the subsequent stress shadowing effect usually adversely affects the fracturing of adjacent wells. With a small initial horizontal stress difference, large injection rate and staggered stage arrangement can achieve ideal stimulation performance. Our results can provide a guidance for optimizing stimulation design in unconventional well-factory while taking into account environmental protection.

A synthetical geoengineering approach to evaluate the largest hydraulic fracturing-induced earthquake in the East Shale Basin, Alberta

On 2019-03-04,the largest induced earthquake(ML4.18)occurred in the East Shale Basin,Alberta,and the underlying physical mechanisms have not been fully understood.This paper proposes a synthetical geoengineering methodology to comprehensively characterize this earthquake caused by hydraulic fracturing.Based on 3D structural,petrophysical,and geomechanical models,an unconventional fracture model is constructed by considering the stress shadow between adjacent hydraulic fractures and the interactions between hydraulic and natural fractures.Coupled poroelastic simulations are conducted to reveal the triggering mechanisms of induced seismicity.It is found that four vertical basement-rooted faults were identified via focal mechanisms analysis.The brittleness index(BI)along two horizontal wells has a high magnitude(BI>0.5),indicating the potential susceptibility of rock brittleness.Due to the presence of overpressure,pre-existing faults in the Duvernay Formation are highly susceptible to fault reactivation.The occurrence of the earthquake clusters has been attributed to the fracturing fluid injection during the west 38^(th)-39^(th) stage and east 38^(th) stage completions.Rock brittleness,formation overpressure,and large fracturing job size account for the nucleation of earthquake clusters,and unconventional natural-hydraulic fracture networks provide fluid flow pathways to cause fault reactivation.This workflow can be used to mitigate potential seismic risks in unconventional reservoirs in other fields.

Experiment on Hydraulic Fracturing in Rock and Induced Earthquake

Experiment on rock hydraulic fracturing strength under different confining pressures was conducted on a series of test specimens with various pre-cracks prepared from 7 types of rock. Combining the data of an actual reservoir-induced earthquake with the experimental results of the contemporary tectonic stress field according to the theory of rock strength and the principle and method of rock fracture mechanics, the authors tentatively investigated the earthquakes induced by pore-water pressure in rock and obtained the initial results as follows: （1） One type of induced earthquake may occur in the case of larger tectonic stress on such weak planes that strike in similar orientation of principle tectonic compressional stress in the shallows of the rock mass; the pore-water pressure σp may generate tensile fracture on them and induce small earthquakes; （2） Two types of induced earthquake may occur in the case of larger tectonic stress, i.e.,① on such weakness planes that strike in similar orientation of principle tectonic compressional stress, σ1, in the shallows of the rockmass, the pore-water pressure, σp, may generate tensile fracture on them and induce small earthquakes; ② When the tectonic stress approximates the shear strength of the fracture, the pore-water pressure σp may reduce the normal stress, σn, on the fracture face causing failure of the originally stable fracture, producing gliding fracture and thus inducing an earthquake. σp may also increase the fracture depth, leading to an induced earthquake with the magnitude larger than the previous potential magnitude; （3） There is a depth limit for each type of rock mass, and no induced earthquake will occur beyond this limit.

Maximum Magnitude of Seismicity Induced by a Hydraulic Fracturing Stage in a Shale Reservoir: Insights from Numerical Simulations

A key unknown limiting assessment of risk posed by inducing anomalous seismicity during hydraulic fracturing is the potential maximum magnitude of an event. To provide insights into the variation in maximum magnitude that can be induced by a hydraulic fracturing stage, worst-case scenarios were simulated in 2D using coupled hydro-geomechanical models. The sensitivity of the magnitude to the hydro-geomechanical properties of the fault and matrix rock were quantitatively compared through parametric analysis. Our base model predicts a maximum event with moment magnitude (Mw) 4.31 and Mw values range from 3.97 to 4.56 for the series of simulations. The highest magnitude is predicted for the model with a longer fault and the lowest magnitude for the model with a smaller Young’s modulus. For our models, the magnitude is most sensitive to changes in the Young’s modulus and length of the fault and least sensitive to changes in the initial reservoir pressure (i.e. pore pressure) and the Poisson’s ratio.

Hydraulic Fracturing Stress Measurements and Their Implication for the Design of the First Concrete Lined Pressure Shaft in Sri Lanka

The Urea Oya Multipurpose Development Project （UMDP） is a water transfer, hydropower and irrigation in the south-eastern part of the central highland region of Sri Lanka. During geotechnical site investigation program 42 hydraulic fracture tests and 42 impression packer tests were carried out in 3 boreholes to about 840 m depth. Based on the stress measurements the minimum and maximum horizontal stress ratios were calculated. In situ stress computations at all the tests were based on the assumption that the principal stress components were vertical （σv） and horizontal （σH and σh, the maximum and minimum, respectively）. The results of the measurements had a direct impact on the design of the major openings bearing a high overburden--the underground powerhouse and the transformer cavern--and revealed a significant optimization potential concerning the selection of the lining system of the pressure shaft.

Weakening effects of hydraulic fracture in hard roof under the influence of stress arch

For the problem of hydraulic fracture propagation when weakening the hard roof in fully mechanized top-coal caving stope of ultra-thick coal seam, based on the stress arch theory and the fracture mechanics, a two-dimensional model for hydraulic fracture of the roof in the stope was established to investigate the propagation laws of hydraulic fracture. The result shows that, after mining, the principal stress direction of overlaying rock deflects to form the stress arch, whose arrow height and arch thickness increase with the increase of the mining width and the side pressure coefficient. Within the influence range of stress arch, the hydraulic fracture in hard roof deflects towards the stope direction in the course of propagation and forms the ‘‘arch" fracture, which cuts off the roof below the fracture in a laminated way. The deflection angle of hydraulic fracture increases with the increase of the mining width, but decreases with the increase of the side pressure coefficient and the fractured horizon. This research can provide theoretical basis for the application of hydraulic fracturing method in the stope roof weakening.

A numerical study of fracture initiation under different loads during hydraulic fracturing

The fracture initiation behavior for hydraulic fracturing treatments highlighted the necessity of proposing fracture criteria that precisely predict the fracture initiation type and location during the hydraulic fracturing process.In the present study,a Mohr-Coulomb criterion with a tensile cut-off is incorporated into the finite element code to determine the fracture initiation type and location during the hydraulic fracturing process.This fracture criterion considers the effect of fracture inclination angle,the internal friction angle and the loading conditions on the distribution of stress field around the fracture tip.The results indicate that the internal friction angle resists the shear fracture initiation.Moreover,as the internal friction angle increases,greater external loads are required to maintain the hydraulic fracture extension.Due to the increased pressure of the injected water,the tensile fracture ultimately determines the fracture initiation type.However,the shear fracture preferentially occurs as the stress anisotropy coefficient increases.Both the maximum tensile stress and equivalent maximum shear stress decrease as the stress anisotropy coefficient increases,which indicates that the greater the stress anisotropy coefficient,the higher the external loading required to propagate a new fracture.The numerical results obtained in this paper provide theoretical supports for establishing basis on investigating of the hydraulic fracturing characteristics under different conditions.

Numerical simulation of fault activity owing to hydraulic fracturing

We built a three-dimensional model to simulate the disturbance of the stress field near the reverse fault in Zhaziao, Leyi Township owing to hydraulic fracturing. The pore pressure, and shear and normal stresses during fracturing are analyzed in detail. Input rock mechanics parameters are taken from laboratory test data of shale samples from the study area. The simulation results suggest that after 16 hours of fluid injection, the pore-pressure variation can activate the reverse fault, i.e., we observe reverse slip, and the shear stress and displacement on the fault plane increase with time. The biggest stress–strain change occurs after one hour of fluid injection and the yield point appears about 0.5 h after injection. To observe the stress evolution in each section, the normal displacement on the boundary is constrained and the fault plane is set as nonpermeable. Thus, the sliding is limited and the shear displacement is only in the scale of millimeters, and the calculated magnitude of the induced earthquakes is between Mw-3.5 and Mw-0.2. The simulation results suggest that fluid water injection results in inhomogeneous fracturing. The main ruptured areas are around the injection positions, whereas the extent of rupturing and cracks in other areas are relatively small. Nevertheless, nonnegligible fault activation is recorded. Sensitivity analysis of the key parameters suggests that the pore pressure is most sensitive to the maximum unbalanced force and the internal friction angle strongly affects the fault slip. Finally, the comparison between the effective normal stress and the maximum and minimum principal stresses on the fault plane explains the fault instability, i.e., the Mohr circle moves towards the left with decreasing radius reduces and intersects the critical slip envelope, and causes the fault to slip.

Fixed-length roof cutting with vertical hydraulic fracture based on the stress shadow effect:A case study

Pre-driven longwall retracement roadway(PLRR)is commonly used in large mine shaft.The support crushing disasters occur frequently during the retracement,and roof management is necessary.Taking the 31107 panel as research background,the roof breaking structure of PLRR is analyzed.It is concluded that the roof cutting with vertical hydraulic fracture(HF)at a specified position,that is,fixed-length roof cutting,can reduce support load and keep immediate roof intact.The extended finite element method(XFEM)is applied to simulate hydraulic fracturing.The results show that both the axial and transverse hydraulic fracturing cannot effectively create vertical HFs.Therefore,a novel construction method of vertical HF based on the stress shadow effect(SSE)is proposed.The stress reversal region and HF orientation caused by the prefabricated hydraulic fracture(PF)are verified in simulation.The sub-vertical HFs are obtained between two PFs,the vertical extension range of which is much larger than that of directional hydraulic fracturing.The new construction method was used to determine the field plan for fixed-length roof cutting.The roof formed a stable suspended structure and deformation of the main PLRR was improved after hydraulic fracturing.