Application of machine learning to characterize gas hydrate reservoirs in permafrost settings
AGU 2021 Fall Meeting, Virtual, December 13-17, 2021.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
AGU 2021 Fall Meeting, Virtual, December 13-17, 2021.
Abstract not provided.
AGU Fall Meeting 2022, Chicago, IL, December 12-16, 2022
Knowledge of petrophysical and geomechanical properties of gas hydrate-bearing sediments are essential for predicting reservoir response due to gas production from gas hydrate reservoirs. That information will be critical parameters for designing production well completion such as specification of depressurization pump, water storage tank, and mesh size of sand screen. In December 2018, Stratigraphic Test Well Hydrate-01 was drilled in the western part of the Prudhoe Bay Unit, Alaska North Slope as part of the technical planning effort for a future long-term production test by collaborative team of DOE/NETL, USGS, and MH21-S (Boswell et al., 2020, Collett et al., 2020, Okinaka et al., 2020). Data set of logging-while-drilling (LWD) were acquired (Haines et al., 2020, Suzuki et al., 2019) and core sampling depth was determined on-site.Side-wall pressure coring was conducted to recover gas hydrate-bearing sediments from two reservoir sections named Unit-B and Unit-D. A total of 34 cores were successfully recovered by 5 runs of a wire-line deployed pressure corer (CoreVault® System - Halliburton). Core analysis scheme of this project are shown in Figure 1. All cores were quenched in liquid nitrogen while at high pressure in the laboratory of Stratum Reservoir, LLC. at Anchorage (Figure 1, a)). And the cores were removed from the pressure corer autoclave with temperature support of dry ice and stored under liquid nitrogen at atmospheric pressure. 19 damaged cores were processed to index property measurements directly such as grain size, grain density. 4 of another 17 cores were depressurized and trimmed for making a plug to analyze petrophysical properties of host sediments. Unsteady-state permeability test was conducted to obtain relative water permeability to gas and core scale NMR T2 distribution measurement was performed for evaluating pore size distribution at Houston (Figure 1, b)). Remained high quality 13 cores were preserved with gas hydrate for advanced laboratory analysis. National Institute of Advanced Industrial Science and Technology, as a part of the Japanese National Hydrate Research Program (MH21-S, funded by Ministry of Economy, Trade and Industry), received the samples at Sapporo, Japan for advanced core analysis. High-resolution X-ray CT was used to analyze the quality of the samples, which showed undisturbed lithological layers. Cores were lathed into cylindrical shape and distributed for multi property measurements (Figure 1, c)).<p>As a result, sediment from Unit-D is categorized as silty sand at ~37% porosity with ~80% gas hydrate saturation. An average hydration number n = 6.16 was measured by Raman spectroscopy. An average intrinsic permeability of ~400 mD and in situ effective permeability (with hydrate) on the order of ~10 mD. The Unit B recovered cores consisted of well sorted sand at ~40% porosity with ~95% gas hydrate saturation. An average intrinsic permeability of ~1 Darcy and in situ effective permeability on the order of ~30 mD was measured for the Unit B cores. There was a small permeability reduction due to porosity loss with increasing effective stress that simulated consolidation behavior along with depressurization in the highly permeable sandy sediment. The apparent minimum change in porosity and permeability may be caused by the low compressibility of quartz sand grains in the recovered cores. XRD and thermal conductivity analysis also suggested high quartz content. Triaxial compression tests established internal friction angles based on the Mohr-Coulomb's failure criterion, which are 40° for hydrate-bearing sediment and 29.8° for hydrate free sediment.</p>
Abstract Gas hydrate is an ice‐like form of water and low molecular weight gas stable at temperatures of roughly −10°C to 25°C and pressures of ~3 to 30 MPa in geologic systems. Natural gas hydrates sequester an estimated one sixth of Earth's methane and are found primarily in deepwater marine sediments on continental margins, but also in permafrost areas and under continental ice sheets. When gas hydrate is removed from its stability field, its breakdown has implications for the global carbon cycle, ocean chemistry, marine geohazards, and interactions between the geosphere and the ocean‐atmosphere system. Gas hydrate breakdown can also be artificially driven as a component of studies assessing the resource potential of these deposits. Furthermore, geologic processes and perturbations to the ocean‐atmosphere system (e.g., warming temperatures) can cause not only dissociation, but also more widespread dissolution of hydrate or even formation of new hydrate in reservoirs. Linkages between gas hydrate and disparate aspects of Earth's near‐surface physical, chemical, and biological systems render an assessment of the rates and processes affecting the persistence of gas hydrate an appropriate Centennial Grand Challenge. This paper reviews the thermodynamic controls on methane hydrate stability and then describes the relative importance of kinetic, mass transfer, and heat transfer processes in the formation and breakdown (dissociation and dissolution) of gas hydrate. Results from numerical modeling, laboratory, and some field studies are used to summarize the rates of hydrate formation and breakdown, followed by an extensive treatment of hydrate dynamics in marine and cryospheric gas hydrate systems.
Here we use 3D seismic data to show that three rows of seafloor gas mounds can be traced downward to leak points that lie at the hydrate-gas contact within three individual dipping coarse-grained sand bodies in the Terrebonne Basin, Walker Ridge block 313, northern Gulf of Mexico. We predict the overpressure within the sand bodies by assuming that the gas pressure at the vent points equal the least principal stress. We interpret that free gas accumulates at the base of the hydrate stability zone, trapped by the overlying sand which has a high methane hydrate saturation. The free gas accumulates until the gas pressure at the base of the hydrate stability zone reaches the least principal stress in the overlying mudrock, whereupon hydraulic fractures form and fluids are vented to the surface. The warm rising fluids and perhaps localized exothermic formation of hydrate raise the local salinity and temperature. This process progressively shifts the base of the hydrate stability zone to shallower depths and dissociates the hydrate seal within the sand, which creates new leak points and results in the observed migration of the seafloor vents. Within the southwest corner of the Terrebonne Basin, this process has repeated multiple times within the Blue, Orange, and Green sands. Here this study shows how the hydrate stability zone can influence the location of fluid expulsion and in turn be affected by the warm, rising, saline fluids sealed by hydrate.
The Ulleung Basin Gas Hydrate field expeditions in 2007 (UBGH1) and 2010 (UBGH2) sought to assess the Basin's gas hydrate resource potential. Coring operations in both expeditions recovered evidence of gas hydrate, primarily as fracture-filling (or vein type) morphologies in mainly silt-sized, fine-grained sediment, but also as pore-occupying hydrate in the coarser-grained layers of interbedded sand and fine-grained systems. A commonality across many of these occurrences is the presence of diatoms in the fine-grained sediment. Here we tested fine-grained sediment (median grain size <12.5 μm) associated with hydrate occurrences at four UBGH2 sites (UBGH2-2-2, UBGH2-3, UBGH2-6 and UBGH2-11) to investigate potential impacts of diatoms on efforts to extract methane from hydrate, or to tap hydrocarbon reservoirs beneath hydrate-bearing sediment. Two key considerations are: the extent to which diatoms control sediment mechanical properties, and the extent to which pore-water freshening, which occurs as gas hydrate breaks down during resource extraction, alters the diatom control on sediment mechanical properties. We conducted experiments to measure sediment index properties, sedimentation behavior and compressibility to address these considerations. We relied on scanning electron microscope (SEM) imagery and X-ray powder diffraction (XRD) to characterize the sediment mineralogy. Our high-level findings are that at the ~20–45% (by volume) diatom concentrations observed at these UBGH2 sites, sediment compressibility increases with diatom content, but diatoms only appear to increase porosity and permeability at the highest diatom concentration (~45%). Our measurements suggest in situ compression indices of 0.35–0.55 and permeabilities on the order of 0.01milliDarcies (1 × 10 -17 m 2 ) can be anticipated at these sites. Importantly, these properties are not expected to vary significantly upon pore water freshening that accompanies gas hydrate dissociation during production.
Accurate modeling of gas hydrate reservoir productivity and geomechanical risks associated with subsurface dissociation of natural gas hydrates (NGH) requires the determination of model parameters through physical testing on natural hydrate-bearing sediments (HBS). This involves investigating the hydro-mechanical behavior of undisturbed hydrate samples from nature under in situ conditions using pressure core characterization and analysis, which provides a unique opportunity for research. By employing state-of-the-art micro computed tomography imagery on cryogenically preserved, hydrate-bearing sediment samples, we can determine hydrate saturation as well as permeability with and without the presence of hydrates in the sediment. Furthermore, utilizing a machine learning based image segmentation technique, it is possible to extract pore space and grain information. Subsections of the entire image volume were used to determine anisotropic permeabilities using a finite-difference method Stokes solver (FDMSS). Additionally, permeability measurements on whole pressure and temperature preserved hydrate-bearing core were analyzed by utilizing the National Energy Technology Laboratory’s (NETL) Pressure Core Characterization and X-ray CT Visualization Tool (PCXT) to manipulate, cut, and analyze pressure preserved sediment. Permeabilities were measured under a broad range of vertical stress states to simulate expected pressure changes during production scenarios, and the results show that permeabilities derived from images are in agreement with those from traditional core derived experiments. The collected stress-dependent permeability, permeability anisotropy, and corresponding gas hydrate saturations provide valuable input into numerical simulations of reservoir productivity. These properties have been proven to be key parameters determining a long-term reservoir response under depressurization.
Reliable estimations of the relative permeability of gas and water in hydrate-bearing sediments (HBS) and the dependency of the relative permeability on hydrate saturation are critical to predict the productivity of a hydrate reservoir. Yet, this remains poorly estimated owing to lack of experimental data associated with difficulties in conducting multiphase flow experiments in HBS. Recognizing the experimental challenges, this study intends to develop and validate a new experimental system and procedure of unsteady-state relative permeability test that can generate reliable and reproducible flow measurements in HBS. Furthermore, gas hydrate is considered as a part of solid matrix in the sediment, so one of the challenges is to maintain a constant hydrate saturation, which is achieved in this experimental study using tight pressure-temperature (P-T) control near the hydrate stability boundary. The measured differential pressure across the specimen, methane injection flow rate, and volume of displaced brine are used to calculate the relative permeability by adopting a conventional Buckley-Leverett theory-based interpretation method. Residual brine saturation calculated for the hydrate-bearing specimen is higher than that of hydrate-free specimen, presumably due to decrease in pore size, increase in heterogeneity of solid matrix, and increase in size distribution of solid matrix and pore in the presence of hydrates. Further studies are necessary to represent the results of the unsteady-state flow experiment in HBS with a gas hydrate-dependent relative permeability model.
Gas molecule clustering within nanopores holds significance in the fields of nanofluidics, biology, gas adsorption/desorption, and geological gas storage. However, the intricate roles of nanoconfinement and surface chemistry that govern the formation of gas clusters remain inadequately explored. In this study, through free energy calculation in molecular simulations, we systematically compared the tendencies of H 2 and CO 2 molecules to aggregate within hydrated hydrophobic pyrophyllite and hydrophilic gibbsite nanopores. The results indicate that nanoconfinement enhances gas dimer formation in the nanopores, irrespective of surface chemistry. However, surface hydrophilicity prohibits the formation of gas clusters larger than dimers, while large gas clusters form easily in hydrophobic nanopores. Despite H 2 and CO 2 both being non-polar, the larger quadrupole moment of CO 2 leads to a stronger preference for dimer/cluster formation compared to H 2 . Here, our results also indicate that gases prefer to enter the nanopores as individual molecules, but exit the nanopores as dimers/clusters. This investigation provides a mechanistic understanding of gas cluster formation within nanopores, which is relevant to various applications, including geological gas storage.
Artificial Neural Networks (ANN) were used to assess methane hydrate occurrence and saturation in marine sediments offshore India. The ANN analysis classifies the gas hydrate occurrence into three types: methane hydrate in pore space, methane hydrate in fractures, or no methane hydrate. Further, predicted saturation characterizes the volume of gas hydrate with respect to the available void volume. Log data collected at six wells, which were drilled during the India National Gas Hydrate Program Expedition 02 (NGHP-02), provided a combination of well log measurements that were used as input for machine learning (ML) models. Well log measurements included density, porosity, electrical resistivity, natural gamma radiation, and acoustic wave velocity. Combinations of well logs used in the ML models provide good overall balanced accuracy (0.79 to 0.86) for the prediction of the gas hydrate occurrence and good accuracy (0.68 to 0.92) for methane hydrate saturation prediction in the marine accumulations against reference data. The accuracy scores indicate that the ML models can successfully predict reservoir characteristics for marine methane hydrate deposits. In conclusion, the results indicate that the ML models can either augment physics-driven methods for assessing the occurrence and saturation of methane hydrate deposits or serve as an independent predictive tool for those characteristics.
In the summer and fall of 2023, the Gulf of Mexico Deepwater Hydrate Coring Expedition (UT-GOM2-2) drilled, cored, made downhole measurements, and analyzed samples from the seafloor to the base of the gas hydrate stability zone in one location (Site H, WR313) in the Terrebonne basin, deepwater Gulf of Mexico. Analyses of data and samples from the expedition will inform biological, geochemical, and geomechanical models to constrain the role of gas hydrates in the carbon cycle and the potential for gas hydrates as an energy resource. Pressure and conventional cores were collected continuously to a depth of 155.1 meters below the seafloor (mbsf). At deeper depths, cores were taken periodically from hydrate-bearing sands and their bounding muds to a total depth of 861.3 mbsf. 162.6 m of conventional core and 54.8 m of pressure core were obtained. Twelve temperature measurements were made between 27.1 and 144.5 mbsf to determine the geothermal gradient. At the seafloor, more than 4 m of sandy silt of unknown origin was encountered. Beneath this sand, to a depth of ~200 mbsf, the section was composed of interbedded mud and biogenic carbonate ooze. The biogenic ooze correlated to low density and high porosity intervals observed in the previously acquired logging while drilling (LWD) data and as measured. Calcareous nannofossil biostratigraphy constrains the entire record to the Pleistocene (< 0.91 million years) with a pronounced increase in sedimentation rate with depth. Beneath 200 mbsf, the section was predominantly composed of mud with two thicker, hydrate-bearing coarse-grained intervals, which are commonly known as the Blue and Orange sands. The dissolved gas concentration was quantified from pressure cores. In the shallow section, dissolved methane concentration increased below the sulfate-methane transition zone (SMTZ) and reaches saturation (the limit of solubility for methane) at 147 mbsf. Gas expansion was very common in conventional and depressurized pressure (conventionalized) cores below the SMTZ. At deeper depths, the methane concentration within muds bounding the Blue and Orange reservoirs was generally found to be less than saturation. The dissolved and hydrate gas composition is consistent with a microbial source, containing greater than 99.99% methane and only trace concentrations of ethane, propane, and butane. The methane to ethane ratio (C 1 /C 2 ) and the methane to ethane plus propane (C 1 /(C 2 +C 3 )) decrease with depth down to at least 678 mbsf, mainly driven by the increase in ethane with depth. It is unclear if this trend continues through the Orange sand interval. The δ 13 C isotopic signature of methane ranges between -69.9 and -78.5 ‰ Vienna Pee Dee Belemnite (VPDB). Pressure core recovery of all sandy intervals was poor. However, pressure core logs of the Orange sand show intervals of low density and high velocity, which are indicative of high hydrate saturation. One core from within the Orange sand was composed of interbedded graded sandy silt and mud. The sandy silts from this core are composed of mainly quartz and feldspar with some lithics. Most of the recovered pressure core samples are maintained at near in-situ pressure and temperature (within the hydrate stability field) at the University of Texas Pressure Core Center awaiting analysis. In the shallow section, samples will be used to determine the flux of organic carbon through the basin system, find the rate at which that carbon was consumed, and understand the microbial population responsible for these processes. In the deeper section, samples from in and around the hydrate reservoirs will be used to determine the petrophysical properties of the reservoir and bounding seals in these systems.
Abstract The formation/precipitation and deposition of pipeline solids, such as gas hydrates, asphaltenes, and waxes have long plagued production fields. Given the vast differences in chemistries of these solids, any current prevention or mitigation strategy, particularly for cases where multiple issues are a concern, is likely to involve an extensive assortment of undefined chemical additives that are both costly and add complexity to the system. Surface treatments (coatings), on the other hand, present a relatively new viable option for management strategies. A chemically and physically robust surface treatment with the ability to address deposition issues for multiple pipeline solids could not only decrease the operating expenditures for a field through material cost savings and obviation of downstream separation, but could also simplify produced fluids by eliminating additional chemicals from the mixture. The purpose of this study is to explore the feasibility of a particular surface treatment as part of a solids management strategy. This work utilizes an omniphobic surface treatment to probe its effects on gas hydrate, asphaltene, and wax deposition. Specifically, high pressure rocking cells are employed to study gas hydrate nucleation and deposition. A bench-scale flowloop filled with crude oil and heptane is used to quantify the deposition of crude oil and asphaltenes after a set time period. Lastly, a mechanical shear device measures the adhesion force of wax deposits on untreated/treated surfaces. The gas hydrate rocking cell tests demonstrate an increase in induction time and occasional elimination of hydrate nucleation with the surface treatment. Moreover, the same apparatus indicates that the critical shear to avoid hydrate deposition may be lowered in the presence of the surface treatment compared to untreated pristine carbon steel coupons. A custom-built asphaltene flowloop then establishes that this surface treatment is effective in dramatically reducing the deposition of the aliphatic components of a crude oil, with a moderate reduction in the asphaltene fraction. Finally, mechanical adhesion force measurements for solidified paraffin wax display trends in agreement with the asphaltene results. The surface treatment on both pristine and corroded surfaces is able to reduce the adhesion of wax deposits to values below a pristine, untreated surface.
Fast 3D X-ray imaging has proved to provide crucial insights into multi-phase dynamic processes in various geomaterials. In particular, many in situ imaging experiments have been made to study gas-hydrate formation in porous sandy samples. Such imaging is challenging for the methane gas-hydrate formation in coal samples because of the coal micro-porosity structure and lower X-ray contrast. Here we present results of the first dynamic in situ micro-computed tomography experiment of methane-hydrate formation in coal samples. Synchrotron phase-contrast tomography techniques allowed to achieve necessary contrast levels to separate all the materials of interest (gas, water, coal, and gas hydrate) in reconstructed images with high spatial and temporal resolution. The imaging results are compared to the ones from a similar tomographic experiment with sand samples. Methane-hydrate formation is accompanied by the water movements caused by cryogenic water suction that happens in sequences of short fast movements with longer equilibrium states in between. Only one type of the hydrate formation was observed in coal (growth as shells on grain boundaries) as opposed to three types in sand (shells, growth into gas pockets, and inside water volumes). In particular, this leads to a slower hydrateformation speed in coal. For the coal sample, we also observed water extraction from grains, and interpreted it as competitive sorption of methane. Lastly, we visualized the dynamic behavior of this water extraction via microchannels inside the coal grains, and performed nano-tomography imaging of these channels for a better understanding of this phenomenon.
In the summer and fall of 2023, the University of Texas (UT) Deepwater Hydrate Coring Expedition (UT-GOM2-2) drilled, cored, made downhole measurements, and analyzed samples from the seafloor to the base of the gas hydrate stability zone at Site H, in the Walker Ridge Protracted Area Block 313 (Site H, WR313), in the Terrebonne Basin, deepwater Gulf of America (Gulf of Mexico). Analyses of data and samples from the expedition will inform biological, geochemical, and geomechanical models to constrain the role of gas hydrates in the carbon cycle and the potential for gas hydrates as an energy resource.
Gas hydrates form at high pressure and low temperatures in marine sediments and permafrost regions of the earth. Despite forming in nanoporous structures, gas hydrates have been extensively studied only in bulk. Understanding nucleation and growth of gas hydrates in nonporous confinement can help create ways for storage and utilization as a future energy source. Herein, we introduce a new method for studying crystal orientation/tilt during tetrahydrofuran (THF) hydrate crystallization under the influence of nano-confinement using polarized Raman spectroscopy. Uniform cylindrical nanometer size pores of anodic aluminum oxide (AAO) are used as a model nano-confinement, and hydrate experiments are performed in a glass microsystem for control of the flash hydrate nucleation kinetics and analysis via in situ polarized Raman spectroscopy. The average THF hydrate crystal tilt of 56 ± 1° and 30.5 ± 0.5° were observed for the 20 nm and 40 nm diameter pores, respectively. Crystal tilt observed in 20 and 40-nanometer-size pores was proportional to the pore diameter, resulting in lower tilt relative to the axis of the confinement at larger diameter pores. Finally, the results indicate that the hydrates nucleation and growth mechanism can depend on the nanoconfinement size. A 1.6 ± 0.01 °C to 1.8 ± 0.01 °C depression in melting point compared to the bulk is predicted using the Gibbs–Thomson equation as a direct effect of nucleation in confinement on the hydrate properties.
Abstract. The pore habits of gas hydrate in natural sediment matrices provide essential clues for understanding physical (mechanical, thermal, hydraulic, and electrical) properties of hydrate-bearing sediments, yet there are no tools that can directly visualize the pore habits of natural gas hydrate other than indirect interpretation based on core-scale or field-scale observations. A significant challenge is to obtain a mini-core from pressure cores retrieved from natural reservoirs for high-resolution micro-CT (computed tomography) scans while maintaining pressure and temperature conditions required for stability of gas hydrate during all operational steps including manipulation, cutting, transferring, sub-coring and CT scanning. We present a new set of tools for pore-scale micro-CT imaging of natural hydrate-bearing sediments while maintaining pressure and temperature control. The tests with laboratory-prepared cores and pressure cores successfully demonstrate the capability of this set of tools to subsample a mini-core from pressure cores, transfer the mini-core to an X-ray transparent core holder, and conduct micro-CT scans. Successfully obtained CT images prove the functionality of this set of tools.