Fei Wang

This paper describes an innovatively designed experimental method for fracturing fluid energy storage to explore the energy storage mechanism during the well shut-in process of fractured shale reservoirs. By improving the existing core clamp and adding fracturing fluid cavities and large volume intermediate containers to simulate artificial fractures and remote shale reservoirs, the pressure changes in the core during the well shut-in process were monitored under the conditions of a real oil–water ratio and real pressure distribution to explore the energy storage law of the shut-in fluid in fractured shale reservoirs. Compared to the 0.62 MPa energy storage obtained from traditional energy storage experiments (without artificial fractures or remote shale reservoirs), the experimental scheme proposed in this paper achieved a 2.45 MPa energy storage, consistent with the field’s monitoring results. The energy storage effects of four fracturing fluids were compared, namely pure CO₂, CO₂ pre-fracturing fluid, slickwater pre-fracturing fluid, and pure slickwater fracturing fluid. Due to the characteristics of a high expansion coefficient and low interfacial tension of pure CO₂, the energy storage effect was the best, and the pressure equilibrium time was the shortest. Considering factors such as comprehensive economy and energy storage efficiency, the optimal range for CO₂ pre-injection is between 20% and 30%. Based on the optimization criterion of energy storage pressure balance, it is recommended that the optimal CO₂ shut-in time be 5 h and the slickwater be 12.8 h. Considering the economic, sand carrying, and energy storage effects, and other factors, CO₂ pre-storage has the best imbibition effect, and the optimal CO₂ pre-storage range is 20~30%. The research results provide theoretical support for energy storage fracturing construction in other shale oil reservoirs of the same type.

Spontaneous imbibition (SI) effectively enhances oil recovery in shale reservoirs, significantly changing well shut-in and flowback design. This study conducted SI and displacement experiments to simulate the well shut-in and flowback stages so that the mechanism of imbibition and displacement between crude oil and fracture fluid can be discussed. In addition, the relative contribution to oil recovery of different types of pores in various stages and the effect of wettability were determined with low-field nuclear magnetic resonance (LF-NMR) via each sample’s T₂ transverse relaxation time at each time. The experimental results show that shale has multiscale pore structure characteristics combined with micropores, small mesopores, and mesopores. During the SI process, crude oil is displaced from micropores by fracture fluid at first, and then a large amount of oil production comes from small mesopores. Oil recovery of water-wet core samples is approximately 40.7%. Oil recovery of oil-wet core samples is about 26%. The wettability significantly affects the imbibition and displacement oil recovery of samples. For the process of SI, oil recovered from small mesopores takes the lead in the complete sample recovery. For the displacement process, oil recovered from small mesopores and mesopores take the lead in the complete sample recovery. After displacement, only 12% of fracture fluid flooded from the samples. This research, demonstrating the imbibition and displacement characteristics of terrestrial shale and several relevant affecting factors, contributes to understanding the fracturing fluid retention mechanism in shale reservoirs and provides crucial theoretical foundations for the development of shale oil reservoirs.

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 existing pump-stopping pressure drop models for the hydraulic fracturing operation of tight oil reservoirs only consider the main hydraulic fracture and the single-phase flow of fracturing fluid. In this paper, a new pump-stopping pressure drop model for fracturing operation based on coupling calculation of the secondary fracture and oil-water two-phase flow is proposed. The physical model includes the horizontal wellbore, the fracture network and the tight oil reservoir. Through the numerical simulation and calculation, the wellbore afterflow performance, the crossflow performance between the main hydraulic fracture and the secondary fracture, the fracturing fluid leakoff and the oil-water replacement after termination of pumping are obtained. The pressure drop characteristic curve is drawn out by the bottom-hole flow pressure calculated through the numerical simulation, and a series of analyses are carried out on the calculated pressure drop curve, which is helpful to diagnose the -oil-water two-phase flow state and the fracture closure performance under the control of the fracture network after hydraulic fracturing pumping. Finally, taking a multi-stage fractured horizontal well in a tight oil reservoir in the Junggar basin, China as an example, the pump-stopping pressure drop data of each stage after hydraulic fracturing are analyzed. Through the history fitting of the pressure drop characteristic curve, the key parameters such as fracture network parameters, which include the half-length of main hydraulic fracture, the conductivity of main hydraulic fracture and the density of secondary fracture, the fracture closure pressure are obtained by inversion, thus, the hydraulic fracturing effect of fractured horizontal well in tight oil reservoirs is further quantified.

CO₂-based fracturing is widely introduced to stimulate shale oil reservoirs for its multiple advantages. However, the range of CO₂ entering the matrix around fractures and CO₂-oil replacement capacity between matrix and fractures cannot be fully explained. To address this issue, a radial constant volume diffusion experiment on shale cores was designed in this study, and the pressure drop curve history was matched through numerical model to determine the composition effective diffusion coefficient. A field-scale numerical model was established, in which a series of certain grids were used to explicitly characterize fracture and quantify the prosess of CO₂ mass transfer and oil replacement. Based on the field-scale numerical model, the process of shut-in, flow back, and oil production was simulated. The distribution of CO₂ in fractured shale oil formation and its impact on crude oil during shut-in stage and flow back stage were investigated. This study concludes that CO₂ gradually exchanges the oil in matrix into fractures and improve the fluidity of oil in matrix until the component concentrations of the whole reservoir reaches equilibrium during the shut-in process. Finally, about 30∼35 mole % of CO₂ in fractures exchanges for oil in matrix. The range of CO₂ entering the matrix around fractures is only 1.5 m, and oil in matrix beyond this distance will not be affected by CO₂. During the process of flow back and production, the CO₂ in fracture flows back quickly, but the CO₂ in matrix is keeping dissolved in oil and will not be quickly produced. It is conclued that the longest possible shut-in time is conducive to making full use of the CO₂-EOR mechanism in fractured shale oil reservoirs. However, due to the pursuit of economic value, a shut-in time of 10 days is the more suitable choice. This work can provide a better understanding of CO₂ mass transfer mechanism in fractured shale oil reservoirs. It also provides a reference for the evaluation of the shut-in time and production management after CO₂ fracturing.

Fracturing fluid imbibition and retention are treated as a main mechanism for oil production from shale reservoirs. However, the oil–water exchange phenomenon during post-fracturing soaking periods has not been thoroughly studied. In this study, a water–oil flow model is built to investigate the water imbibition and oil drainage phenomenon in hydraulically fractured shale. With the developed numerical simulator, the main characteristics of post-fracturing soaking, that is, pressure diffusion, water imbibition, and especially, the oil–water exchange behavior are simulated. Three key time points, that is, oil–water exchange equilibrium, steady exchange efficiency, and oil breakthrough in fracture are found. The oil–water exchange efficiency and exchange volume are also calculated. Moreover, the proposed model is validated by field wellhead pressure dynamics, indicating a relevance of time between the oil–water exchange efficiency and the wellhead pressure falloff derivatives. Finally, the influences of shale permeability, wettability, fracture complexity, and oil viscosity on the oil–water exchange behavior are investigated. Results indicate that the matrix of oil-wet shale almost does not suck water and discharge oil, and only the oil in natural fractures exchanges with the water in hydraulic fractures. The water-wet shale with low permeability, high oil viscosity, and few natural fractures needs extra soaking time to achieve good oil–water exchange performance. The suitable soaking period for the water-wet base case in this study is from 17.25 to 169 days, among which 64 days is the optimal soaking time.

Water‐rock interaction cannot be ignored for shale reservoirs with high‐salinity formation brine and complex rock composition, and stimulated by massive slick‐water fracturing treatment. However, there have been few studies on the flowback model fully coupled with different effects of water‐rock interaction. This paper presents the development of a coupled hydro‐chemical‐mechanical model for modeling water‐rock interaction in fractured shale during the post‐fracturing flowback period. The model considers distinguishing water‐rock interaction phenomena, that is, mineral dissolution, clay swelling and chemical osmosis, and accounts for multi‐phase flow in a fractured shale reservoir. The coupling and solution method and a numerical simulator were developed. Numerical simulation indicates that the swelling volume of clay minerals occupies the pores and leads to a decline in matrix porosity, while mineral dissolution increases both the matrix porosity and the solute concentration in the aqueous phase in matrix pores. Clay swelling mainly affects the shape of the porosity ratio profiles. The effect of mineral dissolution becomes increasingly stronger as flowback progresses. Mineral dissolution mainly affects the relative positions of the porosity ratio curves with the progress of flowback. The water‐rock interaction coupled flowback modeling and the numerical simulation results in this study quantify the effects of chemical osmosis, clay swelling and mineral dissolution. Results from this study provide new insights into the mechanisms of fracturing‐fluid flowback and value to flowback transient analysis.

Production prediction plays an important role in decision making, development planning, and economic evaluation during the exploration and development period. However, applying traditional methods for production forecasting of newly developed wells in the conglomerate reservoir is restricted by limited historical data, complex fracture propagation, and frequent operational changes. This study proposed a Gated Recurrent Unit (GRU) neural network-based model to achieve batch production forecasting in M conglomerate reservoir of China, which tackles the limitations of traditional decline curve analysis and conventional time-series prediction methods. The model is trained by four features of production rate, tubing pressure (TP), choke size (CS), and shut-in period (SI) from 70 multistage hydraulic fractured horizontal wells. Firstly, a comprehensive data preprocessing is implemented, including excluding unfit wells, data screening, feature selection, partitioning data set, z-score normalization, and format conversion. Then, the four-feature model is compared with the model considering production only, and it is found that with frequent oilfield operations changes, the four-feature model could accurately capture the complex variance pattern of production rate. Further, Random Forest (RF) is employed to optimize the prediction results of GRU. For a fair evaluation, the performance of the proposed model is compared with that of simple Recurrent Neural Network (RNN) and Long Short-Term Memory (LSTM) neural network. The results show that the proposed approach outperforms the others in prediction accuracy and generalization ability. It is worth mentioning that under the guidance of continuous learning, the GRU model can be updated as soon as more wells become available.

Post-fracturing well shut-in is traditionally due to the elastic closure of hydraulic fractures and proppant compaction. However, for shale gas wells, the extension of shut-in time may improve the post-fracturing gas production due to formation energy supplements by fracturing-fluid imbibition. This paper presents a methodology using numerical simulation to simulate the hydrodynamic equilibrium phenomenon of a hydraulically fractured shale gas reservoir, including matrix imbibition and fracture network crossflow, and further optimize the post-fracturing shut-in time. A mathematical model, which can describe the fracturing-fluid hydrodynamic transport during the shut-in process, and consider the distinguishing imbibition characteristics of a hydraulically fractured shale reservoir, i.e., hydraulic pressure, capillarity and chemical osmosis, is developed. The key concept, i.e., hydrodynamic equilibrium time, for optimizing the post-fracturing shut-in schedule, is proposed. The fracturing-fluid crossflow and imbibition profiles are simulated, which indicate the water discharging and sucking equilibrium process in the coupled fracture–matrix system. Based on the simulation, the hydrodynamic equilibrium time is calculated. The influences of hydraulic pressure difference, capillarity and chemical osmosis on imbibition volume, and hydrodynamic equilibrium time are also investigated. Finally, the optimal shut-in time is determined if the gas production rate is pursued and the fracturing-fluid loss is allowable. The proposed simulation method for determining the optimal shut-in time is meaningful to the post-fracturing shut-in schedule.