• Energy Research

The global energy transition is a widespread phenomenon that requires international exchange of experiences and mutual learning. Germany’s success in its first phase of energy transition can be attributed to its adoption of smart energy technology and implementation of electricity futures and spot marketization, which enabled the achievement of multiple energy spatial–temporal complementarities and overall grid balance through energy conversion and reconversion technologies. While China can draw from Germany’s experience to inform its own energy transition efforts, its 11-fold higher annual electricity consumption requires a distinct approach. We recommend a clean energy system based on smart sector coupling (ENSYSCO) as a suitable pathway for achieving sustainable energy in China, given that renewable energy is expected to guarantee 85% of China’s energy production by 2060, requiring significant future electricity storage capacity. Nonetheless, renewable energy storage remains a significant challenge. We propose four large-scale underground energy storage methods based on ENSYSCO to address this challenge, while considering China’s national conditions. These proposals have culminated in pilot projects for large-scale underground energy storage in China, which we believe is a necessary choice for achieving carbon neutrality in China and enabling efficient and safe grid integration of renewable energy within the framework of ENSYSCO.

Conventionally, high-pressure water-based fluids have been injected for hydraulic stimulation of unconventional petroleum resources such as tight gas reservoirs. Apart from improving productivity, water-based frac-fluids have caused environmental and technical issues. As a result, much of the interest has shifted towards alternative frac-fluids. In this regard, n-heptane, as an alternative frac-fluid, is proposed. It necessitates the development of a multi-phase and multi-component (MM) numerical simulator for hydraulic fracturing. Therefore fracture, MM fluid flow, and proppant transport models are implemented in a thermo-hydro-mechanical (THM) coupled FLAC3D-TMVOCMP framework. After verification, the model is applied to a real field case study for optimization of wellbore x in a tight gas reservoir using n-heptane as the frac-fluid. Sensitivity analysis is carried out to investigate the effect of important parameters, such as fluid viscosity, injection rate, reservoir permeability etc., on fracture geometry with the proposed fluid. The quicker fracture closure and flowback of n-heptane compared to water-based fluid is advantageous for better proppant placement, especially in the upper half of the fracture and the early start of natural gas production in tight reservoirs. Finally, fracture designs with a minimum dimensionless conductivity of 30 are proposed.

Although hydraulic fracturing has been massively studied and applied as a key technique to enhance the gas production from tight formations, some problems and uncertainties exist to accurately predict and analyze the fracture behavior in complex reservoirs, especially in the naturally fractured reservoirs like shale reservoirs. This paper presents a full 3D numerical model (FLAC3D) to study hydraulic fracturing behavior under the impact of preexisting orthogonal natural fractures. In this numerical model, the hydraulic fracture propagation direction is assumed perpendicular to the minimum principal stress and activated only by tensile failure, whereas the preexisting natural fractures can be activated by tensile or shear failure or a combination of them, and only tensile failure can open the natural fracture as well. The newly developed model was used to study the impact of preexisting orthogonal natural fractures on hydraulic fracturing behavior, based on a multistage hydraulic fracturing operation in a naturally fractured reservoir from the Barnett Shale formation, northwest of Texas in USA. In this multistage operation, two more representative stages, i.e., stage 1 with a relatively large horizontal stress anisotropy of 3.3 MPa and stage 4 with a comparatively small one of 1.3 MPa, were selected to conduct the simulation. Based on the numerical results, one can observe that the interaction between hydraulic and natural fracture is driven mainly by induced stress around fracture tip. Besides, the horizontal stress anisotropy plays a key role in opening the natural fracture. Thus, no significant opened fracture is activated on natural fracture in stage 1, while in stage 4 an opened fracture invades to about 90 m into the first natural fracture. Conversely, the hydraulic fracture length in stage 1 is much longer than in stage 4, as some fluid volume is stored in the opened natural fracture in stage 4. In this work, the shear failure on natural fractures is treated as the main factor for inducing the seismic events. And the simulated seismic events, i.e., shear failure on natural fractures, are very comparable with the measured seismic events.

The meaningful utilization of artificially created multiple fractures in tight formations is associated with the performance behavior of such flow channels, especially in the case of thermal energy extraction from sedimentary geothermal system. In this study, an innovative idea is presented to develop a numerical model for geothermal energy production based on concrete physical performance of an artificially created tensile multi-fracture system in a simplified manner. The state-of-the-art software FLAC3Dplus-TOUGH2MP-TMVOC are integrated to develop a coupled thermo-hydro-mechanical (THM) fictive model for constructing a multi-fracture scheme and estimating heat extraction performance. By incorporating the actual fracture width of newly created subsequent fracture under the effect of stress shadow, cubic law is implemented for fluid flow and geothermal energy production. The results depict that fracture spacing plays a vital role in the energy contribution through multiple fractures. Afterwards, a field case study to design huge multiple hydraulic fractures was performed in the geothermal well GB X1 in North Germany. The attenuation of fracture propagation becomes more significant when massive multiple fracturing operation is performed especially in the case of lower fracture spacing. The fictive model results will be extended to study the geothermal utilization of the North German basin through massive multiple fractures in our future work.

The utilization of geothermal energy can reduce CO2 emissions into the atmosphere. The reinjection of cooled return water from a geothermal field by a closed loop system is an important strategy for maintaining the reservoir pressure and prolonging the depletion of the geothermal reservoir by avoiding problems, e.g., water level drawdown, ground subsidence, and thermal pollution. However, the drawdown of water injectivity affected by physical and chemical clogging may occur in sandstone aquifers, and the reservoir temperature may be strongly affected by the reinjection of large amounts of cooled geothermal water, thus resulting in early thermal breakthrough at production wells and a decrease in production efficiency. In addition to the injection of cooled geothermal water, the injection of CO2 can be used to maintain the reservoir pressure and increase the injectivity of the reservoir by enhancing water–rock interactions. However, the thermal breakthrough and cooling effect of the geothermal reservoir may become complex when both CO2 and cooled geothermal water are injected into aquifers. In this paper, a simplified small-scale multilayered geological model is established based on a low-medium geothermal reservoir in Binhai district, Tianjin. The ECO2N module of the TOUGH2MP simulator is used to numerically simulate temperature and pressure responses in the geothermal reservoir while considering different treatment strategies (e.g., injection rates, temperatures, well locations, etc.). The simulation results show that a high injection pressure of CO2 greatly shortens the CO2 and thermal breakthrough at the production well. A much lower CO2 injection pressure is helpful for prolonging hot water production by maintaining the reservoir pressure and eliminating the cooling effect surrounding the production wells. Both pilot-scale and commercial-scale cooled water reinjection rates are considered. When the water production rate is low (2 kg/s), the temperature decrease at the production well is negligible at a distance of 500 m between two wells. However, when both the production and reinjection rates of cooled return water are increased to 100 m3/h, the temperature decrease in the production well exceeds 10 °C after 50 years of operation.

Carbon dioxide (CO2)-enhanced oil recovery (EOR) has great potential and opportunity for further development, and it is one of the vital carbon capture, utilization, and storage (CCUS) technologies. However, strong heterogeneity is one of the several challenges in developing reservoirs, especially for China’s continental tight oil reserves. This study investigates the effects of heterogeneous porosity and permeability on CO2 flooding evolution in low-permeable tight formation. We simulated CO2-EOR using a numerical model developed on the platform of TOUGH2MP-TMVOC to evaluate the effect of different levels of heterogeneity on oil production, gas storage, and flow behaviors in a tight reservoir, controlled by standard deviation and correlation length. A comparison of nine cases reveals that porosity heterogeneity commonly intensifies flow channeling, and there is an oil production decline with higher standard deviation and longer correlation length of porosity field. In addition, the porosity correlation length has a negligible effect on reservoir performance when the standard deviation is relatively low. Furthermore, strong heterogeneity also has a negative impact on the storage capacity of CO2 and oil production. Notably, as the standard deviation was raised to 0.1, a small sweep region arose with the early CO2 breakthrough, which led to a worse flooding effect. Finally, this study exemplifies that a higher injection/production rate and CO2 alternating N2 injection strategies can improve oil recovery in highly heterogeneous reservoirs.