Shale oil resources, noted for their broad distribution and significant reserves, are increasingly recognized as vital supplements to traditional oil resources. In response to the high fracturing costs and swift decline in productivity associated with shale oil horizontal wells, this research introduces a novel approach utilizing CO₂ for enhanced shale oil recovery in radial boreholes. A compositional numerical simulation method is built accounted for component diffusion, adsorption, and non-Darcy flow, to explore the viability of this technique. The study examines how different factors—such as initial reservoir pressure, permeability, numbers of radial boreholes, and their branching patterns—influence oil production and CO₂ storage. Our principal conclusions indicate that with a constant CO₂ injection rate, lower initial reservoir pressures predominantly lead to immiscible oil displacement, hastening the occurrence of CO₂ gas channeling. Therefore, maintaining higher initial or injection pressures is critical for effective miscible displacement in CO₂-enhanced recovery using radial boreholes. Notably, the adsorption of CO₂ in shale oil results in the displacement of lighter hydrocarbons, an effect amplified by competitive adsorption. While CO₂ diffusion tends to prompt earlier gas channeling, its migration towards areas of lower concentration within the reservoir reduces the extent of channeling CO₂. Nonetheless, when reservoir permeability falls below 0.01 mD, the yield from CO₂-enhanced recovery using radial boreholes is markedly low. Hence, selecting high-permeability “sweet spot” regions within shale oil reservoirs for the deployment of this method is advisable. To boost oil production, utilizing longer and broader radial boreholes, increasing the number of boreholes, or setting the phase angle to 0° are effective strategies. Finally, by comparing the production of shale oil enhanced by CO₂ with that of a dual horizontal well fracturing system enhanced by CO₂, it was found that although the former’s oil production is only 50.6% of the latter, its cost is merely 11.1%, thereby proving its economic viability. These findings present a new perspective for the economically efficient extraction of shale oil, offering potential guidance for industrial practices.

Multistage hydraulic fracturing is one of the most prevalent approaches for shale reservoir development. Due to the complexity of constructing reservoir environments for experiments, numerical simulation is a vital method to study flow behavior under reservoir conditions. In this paper, we propose a numerical model that considers a multistage fractured horizontal well with a hilly-terrain trajectory in a shale reservoir with the presence of natural fractures. The model was constructed based on the MATLAB Reservoir Simulation Toolbox and used the Embedded Discrete Fractured Model (EDFM) to describe the interrelationship between the matrix, fractures, and wellbore. The model was then applied to an actual condensate gas well producing from a shale reservoir, and the effects of reservoir parameters on the simulation data were studied based on this well case. The simulation results were highly consistent with the actual production data, which validates the accuracy of this model and proves its potential for predicting future production trends. We extended the discussion to two examples with extreme well trajectories by reviewing the inflow contribution of each hydraulic fracture with respect to fracture pressure, and the changes in static pressure with time. In conclusion, the proposed model is capable of providing simulation results close to reality and thus guiding field design and operation in the fracturing and drilling process.

The use of Horizontal wells is a common method of underground natural gas storage (UGS), but there is still a need to discuss whether they are more suitable for slotted liner or perforation completions. To address this issue, a numerical model is developed to predict the gas injection rate of horizontal wells while considering the skin factor. Here, a novel uncoupled iteration method is employed to determine the skin factor deriving from turbulence in each time step when the bottom hole pressure is fixed. The uncoupled method begins with an estimate of the initial gas injection rate, which is then used to calculate a turbulent skin factor. This turbulent skin factor is then used to update the gas injection rate, iterating continuously until convergence is achieved. The effects of slotted liner and perforation design parameters, formation damage, and injection pressure on the skin factor are analyzed. The main findings suggest that the error in the gas injection rate calculated by the non-coupled model compared with the coupled model is only 0.6%, yet it can reduce the number of sub-iterations to 1/10 of that required by the coupled model. Moreover, the uncoupled model can provide results within four steps, even when the convergence condition is 10⁻¹⁴. The open area and perforation density play a significant role in determining the connection degree between the horizontal well and the reservoir, with a larger perforation density resulting in a negative skin factor. Perforations are more suitable than slotted liners for reservoirs with severe formation damage, and the difference in skin factor between the two can reach a value of 40.87 when the ratio of the damage zone’s permeability to that of the normal reservoir zone is 0.05. It is easier to reduce turbulence damage in slotted liner completions than perforation completions, with the turbulence damage of the slotted liner being only 15.9% of that of the perforation. However, to avoid damage it is crucial to prevent the screen tube from being plugged in, as it might otherwise rise to three to ten times the original level. This study provides a theoretical basis and practical reference for the application of slotted liner and perforation method in UGS horizontal wells.