Search NASA⌕ Search

SEARCH · Search NASA

Results for “WELL LOGGING”

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.

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Assessment of CO 2 Storage Resources in Depleted Oil and Gas Fields in the Ship Shoal Area, Gulf of Mexico

The Gulf of Mexico is one of the most important regions in the United States for energy resources and infrastructure. Gulf of Mexico federal offshore oil and gas production accounts for 17% of total U.S. crude oil production and 5% of total U.S. gas production (EIA Gulf of Mexico Fact Sheet). This region presents an excellent combination of high need and significant opportunity for large scale geologic storage of CO 2 . The Ship Shoal Area is located about 20 miles offshore Louisiana within the Gulf Coast federal waters. Miocene and Pliocene sediments in the Ship Shoal Area are proven to provide excellent and secure traps for oil and gas. The Ship Shoal Area contains a large number of depleted oil and gas fields either currently abandoned, or planned for abandonment by 2025, which may provide very significant potential CO 2 storage capacity. GeoMechanics Technologies conducted a comprehensive research project to better characterize Neogene sediments in the Ship Shoal Area for high volume CO 2 storage. The data generated from literature, well data, well log, and formation evaluation were input into Rockwork 16 geologic software to create a detailed geologic model for the Ship Shoal study area. The geologic grids were then fed into TOUGH2, the gas migration model software and FLAC 3D, the geomechanical model software. The fluid flow modeling for both the Miocene and Pliocene simulation from both Block 107 and Block 84 indicated that Pliocene and Miocene are a good reservoir for the CO 2 sequestration. Thirty years of CO 2 injection (at a rate of 1 million metric tons CO 2 per year) and 30 years of observation simulations were run with the fluid flow modeling. The 3D geomechanical modeling for Ship Shoal Block 107 and Block 84 were developed to evaluate induced surface deformation and potential fault reactivation after 30 years CO 2 injection at base of Pliocene and at upper Miocene injection locations. Evaluation results indicate low to no risks for fault slips or fault reactivation. Twelve well bores and 76 well bores for Ship Shoal Blocks 84 and 107 respectively were reviewed for their cement history. Most wells (57 out of 76 and 7 out of 12 wells in Ship Shoal Block 107 and 84 respectively) have good integrity. Nineteen and 5 wells (in Ship Shoal Block 107 and 84 respectively) with no top plug, incomplete cement or Plug and Abandonment information are given yellow cautionary indicators. These cautionary wells may provide leakage paths of CO 2 through the well bores to the USDWs. Using our Quantitative Risk & Decision Analysis Tool (QRDAT) for caprock integrity evaluation, we compared Ship Shoal’s risk to that of In Salah, Sleipner, Kevin Dome, Loudon, Illinois Industrial CCS and Wilmington Graben. We found the risk at the Ship Shoal Blocks 84 and 107 fields are similar to the known CO 2 active sequestration sites, but lower than the Wilmington Graben turbidities offshore California studied site. The risk of natural seismicity in the Gulf of Mexico is relatively low. Geomechanics Technologies has documented the top 25 CO 2 emission sources within the close proximity of the Ship Shoal Block 84 and 107 fields. All the offshore and onshore pipelines have been digitized and can be viewed in an interactive website (http://www.geomechanicstech.com/shipshoal.html). We also performed a feasibility study on the potential for converting existing oil or gas pipelines for CO 2 transport. There are abandoned, idled and retired onshore and offshore pipelines. However, there is no standard specification for maximum pipeline pressures needed for CO 2 transport. It is a function of design, materials, and testing provided; thus it will be the responsibility of the pipeline operator to correctly determine, maintain and operate within the limits of the pipelines. There are a few transit corridors extending from onshore Louisiana (Cailou Bay) to offshore trunk-lines. The cost for constructing a pipeline has increased about 46% from 2015 to an average cost of $5,064,046 per mile in 2016. We use the NETL CO 2 storage resource mass estimate to calculate the potential resources for Ship Shoal Block 84 and 107 fields using the depleted oil and gas reservoir volume and sand volume generated from our geologic modeling. The estimated storage resource results are greater when using the sand volumes derived through geologic modeling than the oil and gas reservoir data. The difference is due to the depleted oil and gas reservoir information not accounting for the water-flooded sand located below the oil/gas-water contact or unproductive sand units. Resource calculation using the sand volume obtained through geologic modeling overestimate the storage capacity as the model accounts for all sand within the formation, not just the interconnected sand.

02 PETROLEUM↗
DOE OSTI2020

Wellbore injectivity response to step-rate CO 2 injection: Coupled thermo-poro-elastic analysis in a vertically heterogeneous formation

Safe and permanent carbon geological storages require an elaborate analysis of geomechanical stability. Subsurface injection of CO 2 changes the local temperature and pore pressure and, further, alters the stress due to thermo-poro-elastic responses. A permanent CO 2 storage testing was conducted in the water leg of Cranfield reservoir in Mississippi, USA, where hosted CO 2 enhanced-oil recovery activities. During CO 2 injection, the injection rate was ramped up twice but the bottom-hole pressure did not increase with the imposed injection rates as expected. This unexpected field observation suggests the possibility of an open-mode fracture development at the injector. However, the injector response and potential fracture development have not been rigorously interpreted with detailed near-wellbore temperature and pore pressure change upon the CO 2 injection. In this study, we performed history matching of CO 2 injection using coupled thermo-poro-elastic reservoir simulation to examine the possibility of fracturing. We built a reservoir model including vertical heterogeneity in both petrophysical and geomechanical properties estimated from well-logging analysis and laboratory experiments. The simulation results show that CO 2 injection changes stresses and support the hypothesis of development of an open-mode fracture at the injector during the Cranfield test. The near-injector region exhibits a large temperature reduction up to 55 °C with ensuing effective horizontal stress reduction up to 9.1 MPa. However, a caprock integrity issue is unlikely because of the horizontal stress contrast within layers and low hydraulic communication with the injection zone. This work indicates that the injectant temperature should be considered in the design of high-rate CO 2 injector.

58 GEOSCIENCES↗
DOE OSTI2016

Progress in Monitoring Strategies for Risk Reduction in Geologic CO2 Storage

Abstract: This geophysical monitoring review is intended to lay the ground work for an approach to integrating monitoring strategies with risk assessment and risk reduction and to summarize work conducted by U.S. Department of Energy's National Risk Assessment Partnership (NRAP) participants that addresses key questions related to monitoring strategies. Commercial, demonstration and pilot geologic CO2 storage (GCS) projects have provided a wealth of practical experience in monitoring technologies. In this manuscript we review a few key projects, specific geophysical methods, well-based methods, geophysical workflows and some suggested best practices. These practices are often driven and governed by regulatory requirements (e.g. U.S. Environmental Protection Agency [US-EPA] Class VI underground injection regulation) and the need for effective reservoir management of pressures and volumes as well as the site-specific risks to containment, plume migration, and underground sources of drinking water. Reducing uncertainty in monitoring is likely to lead to reducing uncertainty in risk. One of our conclusions is the need for field experiment testing that is risk-driven. It is notable that GCS pilot projects have been designed for safe storage, and to date they have been successful. Technical challenges at all geologic CO2 storage sites include assuring injectivity and capacity predictions are reliable, the fate and subsurface topology of the CO2 plume is acceptable, regulatory permitting and compliance are accomplished in a timely manner, needed business assurances are supported by quantitative technical data, and the storage operation is acceptable to stakeholders. There is risk associated with all of these challenges necessitating a monitoring program utilizing both direct and indirect monitoring methods. Because of the direct association of risk and monitoring in GCS, geophysical monitoring is viewed as a key component in approaching risk assessment. To address monitoring in the context of risk assessment the U.S. Department of Energy's National Risk Assessment Partnership (NRAP) has undertaken a review of work related to monitoring storage sites. In this review we provide a framework for monitoring and then present some specific monitoring tools that may be useful for regulatory compliance, as well as technologies that are applicable to leakage detection.

Harbert, William↗
NASA NTRS2012

Characterization of the Vectron PX-570 Crystal Oscillator for Use in Harsh Environments

Computing hardware, data-acquisition systems, communications systems, and many electronic control systems require well-controlled timing signals for proper and accurate operation. These signals are, in most cases, provided by circuits that employ crystal oscillators due to availability, cost, ease of operation, and accuracy. In some cases, the electronic systems are expected to survive and operate under harsh conditions that include exposure to extreme temperatures. These applications exist in terrestrial systems as well as in aerospace products. Well-logging, geothermal systems, and industrial process control are examples of ground-based applications, while distributed jet engine control in aircraft, space-based observatories (such as the James Webb Space Telescope), satellites, and lunar and planetary landers are typical environments where electronics are exposed to harsh operating conditions. To ensure these devices produce reliable results, the digital heartbeat from the oscillator must deliver a stable signal that is not affected by external temperature or other conditions. One such solution is a recently introduced commercial-off-the-shelf (COTS) oscillator, the PX-570 series from Vectron International. The oscillator was designed for high-temperature applications and as proof, the crystal oscillator was subjected to a wide suite of tests to determine its ruggedness for operation in harsh environments. The tests performed by Vectron included electrical characterization under wide range of temperature, accelerated life test/aging, shock and vibration, internal moisture analysis, ESD threshold, and latch-up testing. The parametric evaluation was performed on the oscillator's frequency, output signal rise and fall times, duty cycle, and supply current over the temperature range of -125 C to +230 C. The evaluations also determined the effects of thermal cycling and the oscillator's re-start capability at extreme hot and cold temperatures. These thermal cycling and restart tests were performed at the NASA Glenn Research Center. Overall, the crystal oscillator performed well and demonstrated very good frequency stability. This paper will discuss the test procedures and present details of the performance results.

Li, Jacob↗
DOE OSTI2012

Gas Hydrate Characterization in the GoM using Marine EM Methods

In spite of the importance of gas hydrate as a low-carbon fuel, a possible contributor to rapid climate change, and a significant natural hazard, our current understanding about the amount and distribution of submarine gas hydrate is somewhat poor; estimates of total volume vary by at least an order of magnitude, and commercially useful concentrations of hydrate have remained an elusive target. This is largely because conventional geophysical tools have intrinsic limitations in their ability to quantitatively image hydrate. It has long been known from well logs that gas hydrate is resistive compared to the host sediments, and electrical and electromagnetic methods have been proposed and occasionally used to image hydrates. This project seeks to expand our capabilities to use electromagnetic methods to explore for gas hydrate in the marine environment. An important basic science aspect of our work was to quantify the resistivity of pure gas hydrate as a function of temperature at seafloor pressures. We designed, constructed, and tested a highpressure cell in which hydrate could be synthesized and then subjected to electrical conductivity measurements. Impedance spectroscopy at frequencies between 20 Hz and 2 MHz was used to separate the effect of the blocking electrodes from the intrinsic conductivity of the hydrate. We obtained very reproducible results that showed that pure methane hydrate was several times more resistive than the water ice that seeded the synthesis, 20,000 {Ohm}m at 0{degrees} C, and that the activation energy is 30.6 kJ/mol over the temperature range of -15 to 15{degrees} C. Adding silica sand to the hydrate, however, showed that the addition of the extra phase caused the conductivity of the assemblage to increase in a counterintuitive way. The fact that the increased conductivity collapsed after a percolation threshold was reached, and that the addition of glass beads does not produce a similar increase in conductivity, together suggest that while the surface of the gas hydrate grains are not intrinsically conductive, the presence of sand does increase their conductivity. In the field component of this project, we carried out an 18day cruise on the R.V. Roger Revelle in the Gulf of Mexico from 7th-26th October 2008 to collect controlled-source electromagnetic (CSEM) data over four hydrate prospects; blocks AC 818, WR 313, GC 955, and MC 118. During these surveys we deployed 30 ocean bottom electromagnetic (OBEM) recorders a total of 94 times at four survey areas and towed the Scripps Undersea Electromagnetic Source Instrument (SUESI) a total of 103 hours. SUESI transmission was 200 A on a 50 m dipole antenna at heights of 70-100 m above the seafloor. We also towed a neutrally buoyant 3-axis electric field recorder behind the SUESI antenna at a constant offset of 300 m. The use of a towed receiver that is "flown" above the seafloor allowed us to operate in areas where seafloor infrastructure such as wellheads, pipelines, and installed scientific equipment existed. We reduced the data to apparent resistivity psuedosections. The most compelling results come from the hydrate observatory at MC 118, where a localized resistivity anomaly is clearly identified under the southeast crater in an otherwise uniform 1 {Ohm}m background. The data from MC 118 also show that authigenic carbonate does not necessarily express itself as a confounding resistor, as was feared at the start of this project. While the results from the other prospects are much more complicated, the data are well correlated with known geology, and line to line agreement is good. Although these data are not amenable to 1D inversion as was initially hoped, we expect to use a newly developed 2D CSEM inversion code to continue to get useful information from this rich data set.

03 NATURAL GAS↗
DOE OSTI2017

Impact of Three-Phase Relative Permeability and Hysteresis Models on Forecasts of Storage Associated With CO 2 -EOR

Geological CO 2 sequestration in conjunction with enhanced oil recovery (CO 2 -EOR) includes complex multiphase flow processes compared to CO 2 storage in deep saline aquifers. Two of the most critical factors affecting multiphase flow in CO 2 -EOR are three-phase relative permeability and associated hysteresis, both of which are difficult to measure and are usually represented by numerical interpolation models. The purpose of this research is to improve understanding of (1) the relative impacts of different three-phase relative permeability models and hysteresis models on CO 2 trapping mechanisms, and (2) uncertainty associated with these two factors. Four different three-phase relative permeability models and three hysteresis models were applied to simulations of an active CO 2 -EOR site, the SACROC unit located in western Texas. To eliminate potential bias of deterministic parameters, we utilized a sequential Gaussian simulation technique to generate 50 realizations to describe heterogeneity of porosity and permeability, based on data obtained from well logs and seismic survey. Simulation results of forecasted CO 2 storage suggested that (1) the choice of three-phase relative permeability model and hysteresis model led to noticeable impacts on forecasted CO 2 sequestration capacity; (2) impacts of three-phase relative permeability models and hysteresis models on CO 2 trapping are small during the CO 2 -EOR injection period, and increase during the post-EOR CO 2 injection period; (3) the specific choice of hysteresis model is more important relative to the choice of three-phase relative permeability model; and (4) using the recommended three-phase WAG (Water-Alternating-Gas) hysteresis model may increase the impact of three-phase relative permeability models and uncertainty due to heterogeneity.

42 ENGINEERING↗
DOE OSTI2017

High-Dimensional Intrinsic Interpolation Using Gaussian Process Regression and Diffusion Maps

This article considers the challenging task of estimating geologic properties of interest using a suite of proxy measurements. The current work recast this task as a manifold learning problem. In this process, this article introduces a novel regression procedure for intrinsic variables constrained onto a manifold embedded in an ambient space. The procedure is meant to sharpen high-dimensional interpolation by inferring non-linear correlations from the data being interpolated. The proposed approach augments manifold learning procedures with a Gaussian process regression. It first identifies, using diffusion maps, a low-dimensional manifold embedded in an ambient high-dimensional space associated with the data. It relies on the diffusion distance associated with this construction to define a distance function with which the data model is equipped. This distance metric function is then used to compute the correlation structure of a Gaussian process that describes the statistical dependence of quantities of interest in the high-dimensional ambient space. The proposed method is applicable to arbitrarily high-dimensional data sets. Here, it is applied to subsurface characterization using a suite of well log measurements. The predictions obtained in original, principal component, and diffusion space are compared using both qualitative and quantitative metrics. Considerable improvement in the prediction of the geological structural properties is observed with the proposed method.

58 GEOSCIENCES↗
DOE OSTI2017

Probabilistic Risk Assessment of CO2 Trapping Mechanisms in a Sandstone CO2-EOR Field in Northern Texas, USA

This study employs a polynomial chaos expansion (PCE) technique to develop the reduced order models (ROMs) integrated with the Monte Carlo simulations to quantify uncertainties of CO 2 storage associated with enhanced oil recovery (EOR) in a commercial-scale CO 2 -EOR field. The Morrow B sandstone at the Farnsworth Unit (FWU) EOR field in northern Texas was selected as a case study. A 3-D reservoir model, including simulation of a five-spot injection/production well pattern, was constructed based on the data from well logs and cores, and an updated seismic interpretation in the FWU. Porosity and permeability were considered as sources of uncertainty. The 1000-realization heterogeneous random fields of porosity and permeability were generated based on the variogram analysis on the core measurements of the FWU. Three model outputs of interest include the amount of CO 2 trapped by three trapping mechanisms, i.e. hydrodynamic trapping, oil dissolution trapping, and aqueous dissolution trapping. The total of 25 runs were conducted with the simulations of a 20-year CO 2 -EOR period, and a 20-year monitoring period to develop the ROMs between uncertain model input parameters (i.e., porosity and permeability in this study) and model outputs at the end of simulation for each grid cell using PCE approach. The high coefficient of determination (R 2 ) and small normalized root mean square error (NRMSE) between the simulated (using reservoir modeling) and predicted (using ROMs) CO 2 storage indicate that the ROMs are acceptable and reliable to be used for the predictions. Given the 1000 Monte Carlo samplings of the model input parameters, cumulative distribution functions (CDFs) and uncertainty bounds (10 th and 90 th percentiles) of model outputs were estimated based on ROMs. At the end of simulation, the most injected CO 2 at the FWU was stored by hydrodynamic trapping (between 121,400 and 166,860 tonnes), followed by oil dissolution trapping (between 55,303 and 76,816 tonnes), and aqueous dissolution trapping (between 1979 and 2751 tonnes).

58 GEOSCIENCES↗
DOE OSTI2014

A cell-local finite difference discretization of the low-order quasidiffusion equations for neutral particle transport on unstructured quadrilateral meshes

We present a quasidiffusion (QD) method for solving neutral particle transport problems in Cartesian XY geometry on unstructured quadrilateral meshes, including local refinement capability. Neutral particle transport problems are central to many applications including nuclear reactor design, radiation safety, astrophysics, medical imaging, radiotherapy, nuclear fuel transport/storage, shielding design, and oil well-logging. The primary development is a new discretization of the low-order QD (LOQD) equations based on cell-local finite differences. The accuracy of the LOQD equations depends on proper calculation of special non-linear QD (Eddington) factors from a transport solution. In order to completely define the new QD method, a proper discretization of the transport problem is also presented. The transport equation is discretized by a conservative method of short characteristics with a novel linear approximation of the scattering source term and monotonic, parabolic representation of the angular flux on incoming faces. Analytic and numerical tests are used to test the accuracy and spatial convergence of the non-linear method. All tests exhibit O(h{sup 2}) convergence of the scalar flux on orthogonal, random, and multi-level meshes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Coupling fracture facies with in-situ permeability measurements to generate stochastic simulations of tight carbonate aquifer properties: Example from the Lower Cretaceous aquifer, Northern Provence, SE France

The relationships between fracture facies and permeability distribution have been explored using a well-constrained 3D implicit structural model of hemipelagic low porosity/low permeability carbonate from the Northern Provence region, SE France. Fracture and permeability facies were determined using an extensive dataset of 99 hydrogeological wells. Data processing was undertaken using a step-by-step approach, involving: (i) identification of the fracture facies based on well logs detailing structure and stratigraphy; (ii) determination of permeability facies from the a priori correlation between the dimension of the hydraulic radius of influence (deduced from slug test analyses) and the type of reservoir heterogeneity (fissure, fracture, fault zone, etc.); (iii) three dimensional plot of fracture and permeability facies in the geological model using a variographic analysis of data. Thirty-three sequential indicator simulations (SIS) based on geostatistic analyses were realised on both fracture and permeability facies. Finally, a connectivity algorithm was developed to compute the probability of connection between selected infiltration areas and the major aquifer springs via moderate- to high-permeability geological bodies. Key results are: (i) fault zones have a greater role in controlling permeability facies distribution than on fracture facies repartition; (ii) there is little correlation between permeability and fracture facies distributions; (iii) connectivity results highlight the compartmentalisation of aquifers in the Cadarache area, the extensions of permeable geological bodies being limited by the N130 faults.

58 GEOSCIENCES↗
DOE OSTI2023

Aperture size distribution, length, and preferential location of bed-parallel veins in shale

Bed-parallel, calcite-filled veins (BPVs) are common in shale formations, and although they have been widely described in other studies, little is known about their population aperture size distribution. To address this knowledge gap, we analyzed BPV sizes in outcrops and cores from the Vaca Muerta Formation, Neuquén Basin, Argentina; in two cores from the Marcellus Formation, Appalachian Basin, northeast Pennsylvania; and in one core from the Wolfcamp Shale, Delaware Basin, West Texas. Nine out of ten aperture size populations follow a negative exponential distribution, with one following a weak power law. Bed-parallel vein size distribution and intensity vary among formations and within the same shale. We define three groups of distributions: (1) Vaca Muerta outcrops, with the highest BPV intensity and the largest BPVs (cumulative frequency of 4.9 BPVs per meter [BPVs/m] for apertures 0.265 mm to 8.7 cm); (2) Vaca Muerta cores with a similar BPV intensity overall but with no apertures wider than 1.2 cm; and (3) Vaca Muerta, Wolfcamp, and Marcellus cores with the fewest BPVs (cumulative frequency up to 0.63 BPVs/m) and very few wider than 1 cm. Aperture and length in two outcrop data sets are weakly positively correlated and follow power laws with exponents of 0.44 and 0.49. Mechanical interfaces at boundaries between different lithologies exert a strong control on BPV location, with 65–75% of observed interfaces having BPVs along them. Only 25–30% of the BPVs occur at observed material interfaces, however, and unless subtle, unobserved mechanical layering is present, other factors must also control location. BPV intensity and organic richness (TOC) from Vaca Muerta well logs are correlated in some instances but not in others, indicating TOC is not always a good proxy for BPV location or intensity. Furthermore, these findings provide useful information for modeling of hydraulic fracture treatments where BPVs may influence development of the stimulated fracture network, for example by limiting height growth.

58 GEOSCIENCES↗
DOE OSTI2019

High frequency (4th order) sequence stratigraphy of Early Miocene deltaic shorelines, offshore Texas and Louisiana

Well logs and 3D seismic reflection data are integrated to image Lower Miocene depocenters of the upper Texas and westernmost Louisiana coastal and offshore areas. Although previously interpreted at a large scale (3rd order cycles), the detailed stratigraphy and depositional history of the early Miocene succession has not been fully developed. The upper lower Miocene interval from Robulus L (MFS10) to Amphistegina B (MFS9) has been further subdivided into five 4th order cycles (each 200 to 300 m thick) to provide finer scale stratigraphic detail and interpretation of depositional environments. Results of the finer scale correlations are presented in a series of sandstone percent maps. Overall, the maps display a strike-elongate (sub-parallel to the present day coastline) orientation for the sandstone bodies which thin to the southeast (basinward). Net sandstone thickness is relatively low in older cycles and increases in younger cycles suggesting an overall progradation of the deltaic system into the mapped area through time. Continued marine regression and southward movement of the shoreline are also confirmed by seismic amplitude 2 maps. After a major transgression associated with Robulus L deltaic progradation occurred under rising sea level conditions. Deltaic sediments prograded southward and sandstone brought to the upper Texas coast was distributed laterally by longshore currents to form strandplain coast lines. Maximum regression occurred during the interval MFS9_2 to MFS9_3 when a deltaic depocenter formed offshore Texas in east High Island area. However, sandstone thickness progressively decreased laterally to the east (offshore Louisiana) where more shale was present and marine processes reworked deltaic derived sandstones into shore parallel bars. There was a clear retreat of the shoreline during the youngest cycle marking the beginning of transgression associated with Amphistegina B . Detailed stratigraphic interpretation (at a 4th order scale) shows more variability in the dominant shoreline processes during early Miocene than previously thought. Recognition of sedimentary systems variability at a 4th order scale is critical for improved hydrocarbon exploration and in understanding potential future CO 2 storage in the area.

58 GEOSCIENCES↗
DOE OSTI2020

Theoretical power-law relationship between permeability and formation factor

Permeability k and formation factor F are important petrophysical properties that have broad applications to reservoir engineering. Although well-log analysis has advanced reservoir characterizations, direct measurement of k is still not available. Accordingly, resistivity logs and more particularly formation factor have been frequently used to estimate k. In this study, we present a theoretical framework to relate k to F. We use power-law scaling from continuum percolation theory to develop a theoretic relationship between k and F, which is similar in form to a published power-law equation with empirical coefficients. Our result gives a physical interpretation to this empirical relationship. We also compare power-law scaling laws from percolation theory with data derived from experimental measurements and numerical simulations in the literature. Finally, using the Curve Fitting toolbox of MATLAB, the parameters of the power-law equations are determined, and the significance of the non-universal power laws are discussed.

42 ENGINEERING↗
DOE OSTI2018

The Volume of Gas Hydrate-Bound Gas in the Northern Gulf of Mexico

The northern Gulf of Mexico is known to host gas hydrate in submarine sediments. Estimating the amount of gas hydrate for both carbon cycle and resource interest has been ongoing for more than three decades. A large range of estimates (from 0.2 to 680 trillion cubic meters (TCM)) for hydrate-bound natural gas at standard temperature and pressure exists. We bring a new perspective to resource estimates by using ~800 publicly available petroleum industry well logs assessed for natural gas hydrate. Our resulting probabilistic range of the gas hydrate-bound natural gas in the northern Gulf of Mexico ranges from 37 TCM (10th percentile estimate) to 78 TCM (90th percentile estimate) of gas hydrate-bound natural gas, with a mean estimate of 56 TCM. This suggests that the gas hydrate resource density of the northern Gulf of Mexico is 4 times lower than that previously estimated by the Bureau of Ocean Energy Management. Finally, our results also indicate that over half of the gas hydrate in the area is strategically hosted in sand reservoirs, which is significantly higher than that previously estimated.

58 GEOSCIENCES↗
DOE OSTI2019

Hydrate-filled Fracture Formation at Keathley Canyon 151, Gulf of Mexico, and Implications for Non-vent Sites

Near-vertical hydrate-filled fractures are found in subseafloor marine muds, at both advective methane vent sites and at sites without obvious methane and fluid advection (non-vent sites). At non-vent sites, the mechanisms that transport methane to the fractures and control how hydrate-filled fractures form are not well understood. However, these mechanisms are important to establish, as most of Earth's natural gas hydrate is likely bound in marine mud systems. In this work, we focus on understanding the origin of hydrate and how fracture form at non-vent sites by examining previously hydrate-bearing fractures in conventional cores taken from Keathley Canyon 151, U.S. northern Gulf of Mexico, drilled by the Gas Hydrate Joint Industry Project in 2005. We combine information from well logs, sediment cores, and science party results and add new X-ray computed tomography of archival sections and scanning electron microcopy of core samples to develop a conceptual model. We propose that locally generated microbial methane is transported via diffusion from small pores in marine mud into biomineralized burrows with larger pore size in a process called short-range migration. Hydrate forms in burrows once the methane diffuses into them and the dissolved methane concentration exceeds the solubility threshold. When hydrate fills a burrow, heave from additional hydrate growth places stress on the burrow edges, expands the fracture, and creates additional void space in which methane can diffuse and continue forming hydrate. Fractures slowly propagate in the direction of maximum principal stress.

58 GEOSCIENCES↗
DOE OSTI2020

Laboratory Electrical Conductivity of Marine Gas Hydrate

Methane hydrate was synthesized from pure water ice and flash frozen seawater, with varying amounts of sand or silt added. Electrical conductivity was determined by impedance spectroscopy, using equivalent circuit modeling to separate the effects of electrodes and to gain insight into conduction mechanisms. Silt and sand increase the conductivity of pure hydrate; we infer by contaminant NaCl contributing to conduction in hydrate, to values in agreement with resistivities observed in well logs through hydrate–saturated sediment. The addition of silt and sand lowers the conductivity of hydrate synthesized from seawater by an amount consistent with Archie's law. Here, all samples were characterized using cryogenic scanning electron microscopy and energy dispersive spectroscopy, which show good connectivity of salt and brine phases. Electrical conductivity measurements of pure hydrate and hydrate mixed with silt during pressure–induced dissociation support previous conclusions that sediment increases dissociation rate.

58 GEOSCIENCES↗
DOE OSTI2023

Assessing the potential of composite confining systems for secure and long-term CO 2 retention in geosequestration

A potential geologic target for CO 2 storage should ensure secure containment of injected CO 2 . Traditionally, this objective has been achieved by targeting reservoirs with overlying seals-regionally extensive, low permeability units that have been proven capable of retaining buoyant fluid accumulations over geologic time. However, considering that the amount of CO 2 is limited by a decadal injection period, vertical migration of CO 2 can be effectively halted by a composite system of discontinuous shale/silt/mudstone barriers in bedded sedimentary rocks. Here, we studied the impact of depositional architectures in a composite confining system on CO 2 migration and confinement at reservoir scale. We stochastically generated lithologically heterogeneous reservoir models containing discontinuous barriers consistent with statistical distributions of net-sand-to-gross-shale ratio (NTG) and horizontal correlation lengths derived from well log data and observations of producing hydrocarbon fields in Southern Louisiana. We then performed an extensive suite of reservoir simulations of CO 2 injection and post-injection to evaluate the sensitivity of CO 2 plume migration and pressure response of the composite system to a series of geologic and fluid parameters including the lateral continuity of barriers, NTG, permeability anisotropy within the sand body, and capillary pressure contrast between the sand and shale facies. The results indicate that lateral continuity of barriers and NTG are the dominant controls on CO 2 plume geometry and pressure build-up in the reservoir, while the impact of NTG is particularly pronounced. The significance of intraformational barriers becomes apparent as they facilitate the local capillary trapping of CO 2 . Those barriers improve the pore space occupancy by promoting a more dispersed shape of the plume and ultimately retard the buoyancy-driven upward migration of the plume post injection.

58 GEOSCIENCES↗