Currently, volume fracturing of horizontal wells is the main technology for shale oil development. A large amount of fracturing fluid is injected into the formation, but the flowback efficiency is very low. Besides, the impact of fluid retention on productivity is not fully clear. There is still a debate about fast-back or slow-back after fracturing, and the formulation of a reasonable cleanup scheme is lacking a theoretical basis. To illustrate the injected-fluid recovery and production performance of shale oil wells, an integrated workflow involving a complex fracture model and oil-water production simulation was presented, enabling a confident history match of flowback data. Then, the impacts of pumping rate, slick water ratio, cluster spacing, stage spacing and flowback rate were quantitatively analyzed. The results show that the pumping rate is negatively correlated with injected-fluid recovery, but positively correlated with oil production. A high ratio of slick water would induce a quite complex fracture configuration, resulting in a rather low flowback efficiency. Meanwhile, the overall conductivity of the fracture networks would also be reduced, as well as the productivity, which indicates that there is an optimal ratio for hybrid fracturing fluid. Due to the fracture interference, the design of stage or cluster spacing is not the smaller the better, and needs to be combined with the actual reservoir conditions. In addition, the short-term flowback efficiency and oil production increase with the flowback rate. However, considering the damage of pressure sensitivity to long-term production, a slow-back mode should be adopted for shale oil wells. The study results may provide support for the design of a fracturing scheme and the optimization of the flowback schedule for shale oil reservoirs.

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.

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.

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 effectiveness of hydraulic fractures is one of the critical factors determining the production performance of shale oil wells. To investigate the influence of fluid–shale interaction during formation stimulation on the effectiveness of hydraulic fractures and shale oil production, first, a series of conductivity tests were performed on shale cores after soaking in different fluids, revealing the evolution rules of conductivity of propped fracture under different fluid-shale interaction mechanisms. Furthermore, by considering the fluid–shale interaction in the reservoir numerical simulation model, the impact of hydraulic fracture effectiveness evolution on shale oil well productivity and its differences in reservoirs with different permeability anisotropy were analyzed. The research results indicate that the softening of fracture surfaces caused by fluid-shale interaction causes fracture conductivity damage, and the longer the soaking time, the more severe the damage. Generally, the degree of damage on the fracture conductivity caused by water–shale interaction is higher than that caused by supercritical CO2 shale interaction under the same soaking time. When considering the damage of hydraulic fracture effectiveness caused by fluid–shale interaction, the production of shale oil wells is significantly reduced. For shale reservoirs with strong permeability anisotropy, the cumulative oil production loss caused by fluid–shale interaction is more severe than that with weak permeability anisotropy.

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.