Shale gas reservoirs, as representative reservoirs in the Sichuan Basin, have attracted widespread attention regarding development. Using gas reservoir numerical simulation to assist development has greatly improved the work efficiency of workers. However, traditional gas reservoir numerical simulation is widely criticized for its inability to effectively integrate with geological and engineering factors. In this study, we proposed a geological engineering integration method that considers pre-fracturing parameters. We further applied it to a typical well (N03) in a certain block of the Sichuan Basin. The reliability of the method was determined through historical fitting. Based on the N03 geological model, the optimization range of fracturing construction parameters in adjacent areas was determined with the goal of maximizing EUR. Recommended values for widely distributed construction parameter combinations of Class II reservoirs were provided through orthogonal analysis. The influence order of fracturing construction parameters is (1) sand addition strength, (2) cluster spacing, (3) construction displacement, (4) fracture fluid strength, and (5) horizontal segment length. Finally, we compared the simulated data with the actual case. The results showed that an integrated numerical simulation method including geological and engineering factors can comprehensively and accurately assist in reservoir development.

In addition to main fractures, a large number of secondary fractures are formed after the volumetric fracturing of shale gas wells. The secondary fracture properties are so complex, that it is difficult to identify and diagnose by direct monitoring methods. In this study, a new approach to model and diagnose secondary fracture properties is presented. First, a new pressure decline model, which is composed of four interconnected domains, i.e., wellbore, main fractures, secondary fractures, and reservoir matrix pores, is built. Then, the fracturing fluid pumping and post-fracturing soaking processes are simulated. The simulated pressure derivatives reflect five fracture-dominated flow regimes, which correspond to multiple alternating positive and negative slopes of the pressure decline derivative. The results of sensitivity simulation show that the density, permeability, and width of secondary fractures are the main controlling factors affecting the size ratio. Finally, based on the simulated pressure decline characteristics, a diagnostic method for the identification and analysis of secondary fracture properties is formed. This method is then applied to three platform wells in the Changning shale gas field in China. This study builds the correlation between the secondary fracture properties and the shut-in pressure decline characteristics, and also provides a theoretical method for comprehensive post-fracturing evaluation of shale gas horizontal wells.

The dynamic characteristics of shale gas wells are complexly affected by the gas–water two-phase flow. Based on the special flow mechanism of gas–water two-phase flow in shale gas reservoir, this paper establishes a mathematical model for gas–water two-phase flow in shale gas multi-stage fractured horizontal wells, introduces the eigenvalue method and orthogonal transformation, and obtains the analytical solution of the two-phase flow model. The gas–water two-phase flow rules and main influence factors of shale gas wells were identified, further combined with the flowback test characteristics and data of the shale gas wells in southern Sichuan, the characteristic parameters for the evaluation of the gas well flowback effect were determined, and an index system was established for the evaluation of shale gas well flowback. The analysis result shows that the shale gas well flowback effect has a good relationship with its production capacity, which is mainly reflected in the flowback characteristic parameters such as gas breakthrough time, gas breakthrough flowback rate, 30 d flowback rate, and maximum production flowback rate. The shale gas wells with lower flowback factors have a better production capacity than those with higher flowback factors. The flowback evaluation index system can accurately forecast the shale gas well production capacity in its initial stage, and furthermore offer guidance to promptly ascertaining the block development potential and formulating the development schemes.

The irregular outer boundaries of reservoirs and the associated heterogeneous high-permeability zones formed by complex geological environment pose significant challenges in reservoir characterization and performance prediction. These irregular geometries, which are commonly encountered in field applications but often oversimplified in conventional models, can substantially influence fluid flow dynamics and transient pressure behavior. To solve this critical issue, this paper presents a semi-analytical model for studying the transient pressure behavior of irregular heterogeneous reservoirs, focusing on the dynamic interactions between hydraulic fractures and the surrounding matrix. The model integrates Green’s function solutions for matrix flow with finite difference methods to simulate fluid flow within complex fracture networks, capturing the heterogeneity of the reservoir and the irregularity of its boundaries. Specifically, the reservoir is divided into locally homogeneous blocks, and the flow within each block is solved using bounded Green’s functions, while the fracture networks are discretized and solved using finite difference methods. This proposed model significantly reduces computational complexity compared to traditional numerical simulations, while maintaining high accuracy. Subsequently, we conducted comprehensive parameter sensitivity analyses. The calculational results show that a multi-fractured horizontal well in an irregular heterogeneous reservoir can observe the following flow regimes: bilinear flow, elliptical flow, and boundary-dominated flow. Longer fractures and higher conductivity enhance fracture flux by increasing the contact area and reducing flow resistance, respectively. However, these positive impacts are constrained by drainage area limitations as production progresses.

China has abundant shale gas resources with good exploration value and development potential, making it a recent hotspot for exploration and development. It is widely agreed that large-scale hydraulic fracturing is essential for reservoir enhancement in shale formations. However, the evolution of fractures during hydraulic fracturing is highly complex, necessitating research on the influence of various factors on the vertical propagation of hydraulic fractures. Based on geological and engineering parameters from the Luzhou block in southern Sichuan, this study employed the finite element method (FEM) and the cohesive element method to establish a coupled fluid-solid model for the vertical propagation of hydraulic fractures. Numerical simulations were conducted to investigate the interaction between hydraulic fractures and natural weak planes, clarifying the mechanisms involved. This study elucidates how different rock and natural weak plane properties affect the vertical propagation of hydraulic fractures and draws diagrams illustrating these interactions. The research indicated three fracture distribution patterns after the intersection of hydraulic fractures with natural weak planes: passive fractures, ‘I’-shaped fractures, and crossing fractures. The main fractures in these patterns exhibit initial damage and damage evolution characterized by tensile failure. Specifically, in passive fractures, the initial damage and damage evolution of natural weak planes manifest as shear failure. In ‘I’-shaped fractures, the initial damage in natural weak planes is characterized by shear failure, with damage evolution showing tensile failure. Crossing fractures show minimal damage in the weak planes. Under conditions of high natural weak plane cohesive strength, high Young’s modulus, low interlayer rock cohesive strength, high vertical stress difference, low interlayer stress difference, and high intersection angles, crossing fractures tend to form. Conversely, conditions of low natural weak plane cohesive strength, low Young’s modulus, high interlayer rock cohesive strength, low vertical stress difference, high interlayer stress difference, and low intersection angles favor the formation of ‘I’-shaped fractures. Passive fractures form under conditions of low natural weak plane cohesive strength and high vertical stress difference. This study found that Poisson’s ratio has a minimal effect on the vertical expansion of hydraulic fractures under the studied conditions, with natural weak plane strength being the primary control factor for fracture patterns. These findings enhance the theoretical foundation for the vertical propagation of hydraulic fractures in deep shale formations, facilitating the development and implementation of strategies for enhancing production in shale reservoirs with natural weak planes and better optimizing production in different types of shale reservoirs.

Oxidizing stimulation of organic-rich shale reservoirs, as a supplement of hydraulic fracturing, was proposed to enhance shale gas recovery. Previous publications revealed that the interaction between organic-rich shale and oxidative fluids causes the components’ dissolution, which induces lots of pores and microfractures, resulting in rock microfracturing without confined pressure and associated increments of the matrix permeability, and improving unpropped fracture conductivity. However, the enhancement of shale gas recovery with oxidative fluids still lacks an implementation clue targeted for specific engineering problems. In recent years, water–rock interaction inducing microfractures indicates a positive effect of retained fracturing fluid on the stimulation after the fracturing operation, which sheds light in the enhancement of shale gas production by shut-in. The objectives of this study are to provide a new perspective whereby the shut-in performance to enhance shale gas recovery could be increased by the injection of oxidative fluids into the formation during the fracturing operation. Firstly, the mechanisms of shut-in performance increased by oxidative dissolution, which illustrate the increment of the density of fracture networks, the improvement of fracture network conductivity, and the promotion of gas desorption and diffusivity, are demonstrated. Then, the feasibility of using oxidative fluids to increase shut-in performance, which follows the geological and engineering characteristics of organic-rich shale reservoirs, is evaluated. Finally, according to the analysis of production performance for two typical types of shale gas wells, in which one is a low gas production and a high fracturing fluid recovery (LGP-HFR) and the other is a high gas production and a low fracturing fluid recovery (HGP-LFR), a shut-in strategy with oxidative fluids to enhance shale gas recovery is developed. This indicates that the injection of oxidative fluids during the fracturing operation may become a promising and cost-effective approach to enhance shale gas recovery.