The natural gas hydrates detected in Shenhu Area, South China Sea are considered to have great potential for exploitation and to become a research hotspot. However, details of the dynamic ...accumulation and evolution processes of gas hydrate over a long time in this area remain unclear. In this study, we developed a numerical model to predict the accumulation of gas hydrates in the marine sediments by introducing the biogenic methanogenesis module to the framework of the existing simulator TOUGH+HYDRATE. A one-dimensional dynamic model was applied to reproduce the formation of gas hydrate at site SH2 in Shenhu Area. The burial of sediments and associated phenomena (e.g., evolution of temperature, sediment compaction and consequent reduction in sediment porosity and permeability, fluid expulsion) are taken into account. Modeling results indicate that with the conventional conversion coefficient (CC) of organic carbon to methane of 0.15, the amount of in-situ biogenic methane is not sufficient to form the highly saturated hydrate detected at this location and the remaining methane (91.4%) is supplied by the upward flow of fluids. The tectonic movement that took place about 1.5 Ma created pathways for the fluids to migrate vertically. The pore water and methane gas fluxes in the upward fluids were determined as 2 × 10−10 kg/s·m2 and 1 × 10−11 kg/s·m2, respectively, using the proved hydrate saturation profile and measured pore water chlorinity from recovered core samples. The distribution of hydrate in the sediments is significantly influenced by pore water salinity. The inhibitory effect of salt on hydrate formation, which shifts the pressure of three-phase (i.e., solid hydrate, methane gas and pore water) equilibrium to actual pore water pressure, could be the major reason for the appearance of a three-phase co-existence zone in the hydrate concentrated zone. Additionally, the water flux and diffusivity of salt have a profound effect on the formation of gas hydrate by affecting the salinity of pore water in the gas hydrate concentrated zone. The results of this study are of great significance for guiding marine gas hydrate exploration and resource evaluation.
•Biogenic methane generation module was embedded in the TOUGH+Hydrate code.•The fluxes of pore water and methane in the upward fluids flow were determined using the measured data.•The formation mechanism for the three-phase (water, methane gas and hydrate) co-existence zone was ascertained.•The influence of pore water salinity on hydrate formation and distribution was highlighted.
Display omitted
•A dimensionality reduction model (DRM) is established on the observable tracer breakthrough curve.•DRM does not rely on prior fracture geometry measurements, which is more feasible ...in practice.•DRM allows for integrated fracture network characterization and transport prediction.•DRM concept is raised from the visualization of flow and transport in 3D-printing fracture model.
Dimensionality reduction models (DRM) are widely utilized for well-posed inversion of highly heterogeneous aquifer properties. However, these models predominantly rely on predefined aquifer properties, which are often unknown in realistic aquifers; moreover, DRM for fracture network characterization is still lacking. In this study, we develop a novel DRM based on the tracer breakthrough curve (TBC) obtained from cross-hole tracer test. This model enables integrated characterization of fracture networks and prediction of flow and transport in fractured aquifers. Firstly, the relationship between TBC shape and flow and transport in fracture networks is elucidated through 3D-printing physical modeling and numerical simulations. It is observed that in complex fracture networks, flow and transport dominantly occur through a limited number of preferential flow paths under pumping scenarios involving simultaneous injection and extraction. The preferential flow paths remain stable under different pumping rates, with their number aligning with the number of inflection points, including concentration peaks in TBC and evident slowness point during the decline of tracer concentration. Building upon this finding, DRM composed of parallel fractures with the number of fractures equal to the number of inflection points in TBC is constructed. Parallel fractures in DRM represent preferential flow paths for flow and transport within the original fracture network. Parameters including sizes, apertures, and separation of parallel fractures are then inferred by TBC inversion. It is tested by the 3D-printing physical model that under pumping rates differential to that implemented in the tracer test for DRM parameters estimation, the established DRM is capable of predicting the tracer concentration and outflow temperature at the extraction well with a relative error lower than 10%. Moreover, the implementation of proposed DRM concept with synthetic large fracture network demonstrates its genericity for the characterization of fracture networks at field scale. This new DRM concept provides a feasible and effective way for fracture network characterization and flow and transport prediction.
Thermo‐hydro‐mechanical‐chemical (THMC) interactions are prevalent during CO2 geological sequestration (CGS). In this study, a sequential coupling THMC numerical simulation program was constructed, ...which can be used to explore the following issues of CGS: fluid and heat flow, solute transport; stresses, displacements and rock failures related to geo‐mechanical effects; equilibrium and kinetic chemical reactions; chemical damage to mechanical properties of the rock. Then, the coupling program was applied to the Ordos CGS Project to study the formation response under the multi‐field interaction caused by CO2 injection. The simulation results show that the mechanical process dominates the short CO2 injection period. Specifically, the formation's permeability near the injection well increases by 43%, due to the reduction of effective stress, which significantly promotes the lateral migration of CO2. When the injection rate exceeds 0.15 million tons per year, the cohesion of the reservoir rock is not enough to resist the shear force inside the rock and rock failure may occur. During the subsequent long‐term sequestration period (200 years), the influence of mineral reactions gradually increases. Due to calcite dissolution, the shear modulus of caprock is predicted to decrease by 7.6%, which will to some extent increase the risk of rock failure.
As an important source of low-carbon, clean fossil energy, natural gas hydrate plays an important role in improving the global energy consumption structure. Developing the hydrate industry in the ...South China Sea is important to achieving ‘carbon peak and carbon neutrality’ goals as soon as possible. Deep-water areas subjected to the action of long-term stress and tectonic movement have developed complex and volatile terrains, and as such, the morphologies of hydrate-bearing sediments (HBSs) fluctuate correspondingly. The key to numerically simulating HBS morphologies is the establishment of the conceptual model, which represents the objective and real description of the actual geological body. However, current numerical simulation models have characterized HBSs into horizontal strata without considering the fluctuation characteristics. Simply representing the HBS as a horizontal element reduces simulation accuracy. Therefore, the commonly used horizontal HBS model and a model considering the HBS’s fluctuation characteristics with the data of the SH2 site in the Shenhu Sea area were first constructed in this paper. Then, their production behaviors were compared, and the huge impact of the fluctuation characteristics on HBS production was determined. On this basis, the key parameters affecting the depressurization production of the fluctuating HBSs were studied and optimized. The research results show that the fluctuation characteristics have an obvious influence on the hydrate production of HBSs by affecting their temperatures and pressure distributions, as well as the transmission of the pressure drop and methane gas discharge. Furthermore, the results show that the gas productivity of fluctuating HBSs was about 5% less than that of horizontal HBSs. By optimizing the depressurization amplitude, well length, and layout location of vertical wells, the productivity of fluctuating HBSs increased by about 56.6%.
Widely employed in hydrate exploitation, the single well method is utilized to broaden the scope of hydrate decomposition. Optimizing the well structure and production strategy is necessary to ...enhance gas recovery efficiency. Complex wells represented by the multilateral wells have great application potential in hydrate mining. This study focused on the impact of multilateral well production methods on productivity, taking the Nankai Trough in Japan as the study area. The spatial distribution of physical parameters such as porosity, permeability, and hydrate saturation in the Nankai Trough has significant heterogeneity. For model accuracy, the Sklearn machine learning and Kriging interpolation methods were used to construct a three-dimensional heterogeneous geological model to describe the structure and physical property parameters in the study area of the hydrate reservoir. The numerical simulation model was solved using the TOUGH + Hydrate program and fitted with the measured data of the trial production project to verify its reliability. Finally, we set the multilateral wells for hydrate high saturation area to predict the gas and water production of hydrate reservoir with different exploitation schemes. The main conclusions are as follows: ① The Sklearn machine learning and Kriging interpolation methods can be used to construct a three-dimensional heterogeneous geological model for limited site data, and the fitting effect of the heterogeneous numerical simulation model is better than that of the homogeneous numerical simulation model. ② The multilateral well method can effectively increase the gas production rate from the hydrate reservoir compared with the traditional single well method by approximately 8000 m3/day on average (approximately 51.8%). ③ In the high saturation area, the number of branches of the multilateral well were set to 2, 3, and 4, and the gas production rate was increased by approximately 51.8%, 52.5%, and 53.5%. Considering economic consumption, the number of branching wells should be set at 2–3 in the same layer.
Most hydrate-bearing sediments worldwide exhibit mixed pore- and fracture-filling hydrates. Due to the high exploitation value, pore-filling hydrate production is the focus of current hydrate ...production research, and there is a lack of systematic research on the decomposition of fracture-filling hydrates and their effects on the evolution of temperature and pressure in hydrate-bearing sediments. If only the decomposition characteristics of pore-filling hydrates are studied while the fracture-filling hydrates decomposition and its effects on the hydrate-bearing sediments production process are ignored, the obtained research results would be inconsistent with the actual situation. Therefore, in this study, the effects of fracture-filling hydrates with different dipping angles on the hydrate production process were studied, and the necessity of considering the phenomenon of mixed pore- and fracture-filling hydrates in hydrate-bearing sediments was illustrated. On this basis, the simulation of a typical site (GMGS2-16) with mixed pore- and fracture-filling hydrates was constructed, and the production process was researched and optimized. The results indicated that: (a) fracture-filling hydrates formed in shallow fine-grained sediments and gradually approached the area of pore-filling hydrates, before a stable mixed zone was formed; (b) the occurrence of fracture-filling hydrates was conducive to the hydrate-bearing sediment depressurization production, and the promoting effect of the fracture-filling hydrate with smaller dipping angles was stronger; and (c) depressurization combined with heat injection could effectively compensate for the local low temperature and secondary hydrate caused by the mass decomposition of fracture-filled hydrates.
Gas hydrates are likely to become an important strategic resource with commercial development prospects. It is therefore of great significance to realize the long-term and efficient production of ...methane hydrate reservoirs. Previous studies have shown that the lithological characteristics of hydrate reservoirs have a significant impact on reservoir productivity by influencing the evolution of seepage parameters in the process of hydrate production. The porosity (Φ) and initial hydrate saturation (SH) affect the amount of hydrate decomposition and pressure transfer, and also indirectly affect the reservoir temperature field. The permeability (k) directly affects the rate of pressure-drop transmission and methane gas discharge. Due to the differences in seepage parameters caused by different reservoir lithology, a sandy hydrate reservoir (SHR) in Japan and a clayey silt hydrate reservoir (CHR) in China were found to have different gas production rates and the spatial evolution characteristics of the temperature and pressure fields varied in gas hydrate production tests. Therefore, to ensure the long-term and efficient production of the CHR in China, two models were established for a comparative analysis based on a numerical simulation. The two models were depressurizing models of the CHR of the W11 drilling site in the Shenhu Sea area of the South China Sea and the SHR of the AT1 drilling site in the Eastern Nankai Trough of Japan. Both models considered the heterogeneity of seepage parameters, and the TOUGH+HYDARATE (T+H) code was used in subsequent calculations. Four key results were obtained: (a) The order of the significance levels of the lithological parameters on productivity was k > SH > Φ in the CHR and SH > k > Φ in the SHR. (b) The heat conduction and heat convection in the CHR were weaker than in the SHR, which made it difficult to recover the low-temperature area caused by hydrate decomposition. (c) The exploitation of a high k hydrate reservoir should be given priority when the other initial conditions were the same in both the CHR and SHR. (d) The exploitation of both the CHR and SHR should not only rely on the hydrate content or seepage capacity to determine the reservoir exploitation potential, but the combined effect of the two parameters should be fully considered.
To further understand the interactions of CO2-brine-rock at geological time scales, in this study, a 1D reactive transport model of CO2 intrusion into sandstone of the Longtan Formation (P2l) in the ...Huangqiao area, China, was constructed based on site-specific data. The simulation time is consistent with the retention time of CO2 in the Longtan sandstone Formation and is set to 20 Ma. The reactive transport model is calibrated and revised using the measured data for sandstone samples from Well X3 (i.e., the natural analogue). By comparing the simulation results with measured data for the natural analogue, the long-term geochemical reactions are investigated. The simulation results indicate that the brine-rock interactions induced by CO2 can be roughly divided into two stages. First, susceptible minerals (e.g., chlorite, ankerite, calcite, and feldspar minerals) dissolve rapidly under acidic conditions formed by the dissolution of CO2. The precipitation of siderite is facilitated by the dissolution of ankerite and chlorite. Smectite-Ca and dawsonite precipitate due to the dissolution of anorthite and albite, respectively. Dawsonite begins to convert into smectite-Na when albite is completely dissolved. As the reactions continue, intermediate products (i.e., illite, smectite-Na, and smectite-Ca) generated in the first stage become the reactants and subsequently react with CO2 and brine. These three clay minerals are not stable under acidic conditions and transform into kaolinite and paragenetic quartz in the later stage of reaction. Comparing the simulation results of the Base Case with the measured data for the natural analogue and inspired by previous studies, the scour of kaolinite is supposed to have occurred in this region and is considered in the revised model by introducing a coefficient of the scour of kaolinite (i.e., Case 2). The simulation results of Case 2 fit well with the measured data on mineral assemblage, and the trend of the sandstone porosity growth caused by the CO2-brine-rock reaction is captured by our simulation results. The combination of numerical simulation and natural analogue study indicates that the joint effects of long-term CO2-brine-rock reactions and scour of kaolinite increase the pore space of the host rock and result in an increase in quartz content in the sandstone.
Geothermal energy is clean and independent to the weather and seasonal changes. In China, the huge demanding of clean energy requires the geothermal energy exploitation in the reservoir with depth ...larger than 1000 m. Before the exploitation, it is necessary to estimate the potential geothermal energy production from deep reservoirs by numerical modeling, which provides an efficient tool for testing alternative scenarios of exploitation. We here numerically assess the energy production in a liquid-dominated middle-temperature geothermal reservoir in the city of Tianjin, China, where the heat and fluid transport in the heterogeneous reservoir and deep wellbores are calculated. It is concluded that the optimal injection/production rate of the typical geothermal doublet well system is 450 m3/h, with the distance between geothermal doublet wells of 850 m. The outflow temperature and heat extraction rate can reach 112°C and 43.5 MW, respectively. Through decreasing injection/production rate lower than 450 m3/h and optimizing layout of the injection well and production well (avoiding the high permeability zone at the interwell sector), the risk of heat breakthrough can be reduced. If the low permeability zone in the reservoir is around injection well, it usually leads to abnormal high wellhead pressure, which may be solved by stimulation technique to realize stable operation. The methodology employed in this paper can be a reference for a double-well exploitation project with similar conditions.
CO2 is now considered as a novel heat transmission fluid to extract geothermal energy. It can achieve both the energy exploitation and CO2 geological sequestration. The migration pathway and the ...process of fluid flow within the reservoirs affect significantly a CO2 plume geothermal (CPG) system. In this study, we built three-dimensional wellbore-reservoir coupled models using geological and geothermal conditions of Qingshankou Formation in Songliao Basin, China. The performance of the CPG system is evaluated in terms of the temperature, CO2 plume distribution, flow rate of production fluid, heat extraction rate, and storage of CO2. For obtaining a deeper understanding of CO2-geothermal system under realistic conditions, heterogeneity of reservoir’s hydrological properties (in terms of permeability and porosity) is taken into account. Due to the fortissimo mobility of CO2, as long as a highly permeable zone exists between the two wells, it is more likely to flow through the highly permeable zone to reach the production well, even though the flow path is longer. The preferential flow shortens circulation time and reduces heat-exchange area, probably leading to early thermal breakthrough, which makes the production fluid temperature decrease rapidly. The analyses of flow dynamics of CO2-water fluid and heat may be useful for future design of a CO2-based geothermal development system.