Search NASA⌕ Search

SEARCH · Search NASA

Results for “geophysical log”

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.

98 records · Page 6

Pore characterization of the Marcellus Shale by nitrogen adsorption and prediction of its gas storage capacity

In a shale gas reservoir, pore characterization is an important factor used to determine gas storage capacity. However, the nanometer (nm)-scale pore system in shale is difficult to explore by traditional optical, scanning electron microscopy, or even nuclear magnetic resonance well logging. We have investigated the pore structure and storage capacity of the Marcellus Shale through integration of petrophysical analysis from laboratory and well-logging data and nitrogen adsorption. The isotherm of Marcellus Shale is a composite isotherm, which has features of type I, type II, and type IV isotherms with type H4 of the hysteresis loop, suggesting slit-like pores developed in the Marcellus Shale. Quantitative analysis of pore volumes from the nitrogen adsorption indicates that density porosity may be more properly used to approximate the shale porosity and estimate the shale gas volume. In addition, the specific surface area, micropore, and mesopore volumes have a positive relationship with the kerogen volume and total organic content (TOC). By using the Langmuir and Brunauer-Emmet-Teller models, the simulated result indicates that the higher adsorbed quantity of the Marcellus Shale could be the result of the increase of micropore volume contributed, by the increase of kerogen or TOC content. The proposed equations rapidly compute TOC, a key parameter to predict gas storage capacity in overmature shale such as the Marcellus Shale.

Geochemistry & Geophysics↗

UT-GOM2-2 Preliminary Report: Terrebonne Basin Northern Gulf of Mexico, 30 July-28 September 2023

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.

03 NATURAL GAS↗

Characterization of elastic mechanical properties of Tuscaloosa Marine Shale from well logs using the vertical transversely isotropic model

To avoid steep declines in the Tuscaloosa Marine Shale (TMS) production, wells are fracture-stimulated to release the hydrocarbons trapped in the matrix of the formation. An accurate estimation of Young’s modulus and Poisson’s ratio is essential for hydraulic fracture propagation. In addition, ignoring the highly heterogeneous and anisotropic character of TMS can lead to erroneous stress values, which subsequently affect hydraulic fracture width estimates and the overall hydraulic fracturing process. We have developed an empirical 1D geomechanical model that takes into account VTI anisotropy, and it is used to characterize the elastic mechanical properties of TMS in two wells. In the analyzed formation, the vertical Poisson’s ratio is less than the horizontal Poisson’s ratio, which suggests the necessity of an alternative to the ANNIE equations. The stiffness coefficients [Formula: see text] and [Formula: see text] were estimated using the relationships developed from the ultrasonic core data available for the two TMS. Further, correlations between the static and dynamic properties from laboratory tests were used to improve the minimum horizontal stress calculation. We compare VTI Young’s moduli, Poisson’s ratios, and minimum horizontal stress with the isotropic solution. VTI modeling improves the estimation of the elastic mechanical properties. The isotropic solution underestimates the minimum horizontal stress in the formation. Moreover, it was shown that the 20 ft shale interval below the TMS base is characterized by a low Young’s modulus (the vertical Young’s modulus is equal to 20 GPa, whereas the horizontal Young’s modulus is equal to 40 GPa) and may be a frac barrier.

Geochemistry & Geophysics↗

Deep learning to estimate permeability using geophysical data

Time-lapse electrical resistivity tomography (ERT) is a popular geophysical method to estimate three-dimensional (3D) permeability fields from electrical potential difference measurements. Traditional inversion and data assimilation methods are used to ingest this ERT data into hydrogeophysical models to estimate permeability. Due to ill-posedness and the curse of dimensionality, existing inversion strategies provide poor estimates and low resolution of the 3D permeability field. Recent advances in deep learning provide us with powerful algorithms to overcome this challenge. This paper presents a deep learning (DL) framework to estimate the 3D subsurface permeability from time-lapse ERT data. To test the feasibility of the proposed framework, we train DL-enabled inverse models on simulation data. Each measurement in both synthetic and field data is standardized by removing the mean and scaling the time-series to unit variance. This pre-processing step is necessary to bring simulation data closer to field observations. Subsurface process models based on hydrogeophysics are used to generate this synthetic data. Training performed on limited simulation data resulted in the DL model over-fitting. An advanced data augmentation based on mixup is implemented to generate additional training samples to overcome this issue. This mixup technique creates weakly labeled (low-fidelity) samples from strongly labeled (high-fidelity) data. The weakly labeled training data is then used to develop DL-enabled inverse models and reduce over-fitting. As both time-lapse ERT (1133048 features/realization) and 3D permeability (585453 features/realization) data samples are from a high-dimensional space, principal component analysis (PCA) is employed to reduce dimensionality. Encoded ERT and encoded permeability are generated using the trained PCA estimators. A deep neural network is then trained to map the encoded ERT to encoded permeability. This mixup training and unsupervised learning allowed us to build a fast and reasonably accurate DL-based inverse model under limited simulation data. Results show that proposed weak supervised learning can capture salient spatial features in the 3D permeability field. Quantitatively, the average mean squared error (in terms of the natural log) on the strongly labeled training, validation, and test datasets is less than 0.5. The R 2 -score (global metric) is greater than 0.75, and the percent error in each cell (local metric) is less than 10%. Finally, an added benefit in terms of computational cost is that the proposed DL-based inverse model is at least O(10 4 ) times faster than running a forward model once it is trained. Data generation, DL model training, and hyperparameter tuning to identify optimal neural network architectures utilized high-performance computing resources while the DL inference is performed on a standard laptop. Approximately, O(10 5 ) processor hours are used for generating data and DL tuning and training. We acknowledge that the data generation and DL model development are expensive. But once a DL model is trained, it can be re-used for inversion rapidly for the given system, with set physics and domain. Note that traditional inversion may require multiple forward model simulations (e.g., in the order of 10 to 1000), which are very expensive. This computational savings ≈ O(10 5 ) – O(10 7 )) makes the proposed DL-based inverse model attractive for subsurface imaging and real-time ERT monitoring applications due to fast and yet reasonably accurate estimations of permeability field.

58 GEOSCIENCES↗

Gas hydrate quantification in Walker Ridge block 313, Gulf of Mexico, from full-waveform inversion of ocean-bottom seismic data

The Gulf of Mexico (GOM) Joint Industry Project Leg 2 logging-while-drilling data in Walker Ridge lease block 313 (WR313) in the GOM detected gas hydrate in coarse- and fine-grained sediments at sites WR313-G and WR313-H. The coarse-grained units are thin ([Formula: see text]) and highly saturated, whereas the fine-grained unit is thick (approximately 200 m) with low saturation and fracture-filling gas hydrate. Unlike its coarse-grained counterpart, the seismic character of the fine-grained unit does not clearly indicate the presence of gas hydrate, which would likely have remained undiscovered in the absence of drilling. In this paper, through frequency-domain acoustic full-waveform inversion (FWI) of ocean-bottom seismometer data along a 2D multichannel seismic transect near sites WR313-G and WR313-H, we detect and quantify gas hydrate in the fine-grained unit. Key results are as follows: First, the base of the gas hydrate stability zone, which is not obvious in the reflection profile, can be discerned in the FWI results. Second, the gas hydrate in the fine-grained unit is mainly confined to the area between two sets of opposite-dipping normal faults implying that the fault architecture may be partially responsible for this gas hydrate accumulation and distribution.

Geochemistry & Geophysics↗

DOE Award No.: DE-FE0023919 Phase 4 Scientific/Technical Report

This is the Phase 4 Report for the ‘Deepwater Methane Hydrate Characterization and Scientific Assessment or Genesis of Methane Hydrates in the Gulf of Mexico (GOM2)’ research project (DOE Award No. DE-FE0023919). The report summarizes activities from October 1, 2019 to September 30, 2020. The project is led by the University of Texas at Austin (UT). The project objective is to gain insight into the nature, formation, occurrence and physical properties of methane hydrate-bearing sediments for the purpose of methane hydrate resource appraisal through the planning and execution of drilling, coring, logging, testing and analytical activities that assess the geologic occurrence, regional context, and characteristics of marine methane hydrate deposits in the Gulf of Mexico outer continental shelf (OCS). We published a dedicated American Association of Petroleum Geologists Bulletin volume describing initial results from the UT-GOM2-1 expedition in Sept. 2020. This is part 1 of a multi-volume commitment by AAPG to this project. We further confirmed that the natural gas in hydrate at GC-955 was formed by primary microbial processes (>76.1 %). The in-situ effective permeability hydrate-bearing sandy silts at the GC-955 reservoir ranges from 0.1 md (1.0×10-16 m2) to 2.4 md (2.4×10-15m2) in cores with 83% to 93% hydrate saturation. The intrinsic permeability (the single phase permeability) is estimated from reconstituted samples to be ~12 md (1.2×10-14 m2) to ~41 md (4.1×10-14 m2). We used observation and models to interpret that the core degradation that is found in pressure cores is due to dissociation of the methane hydrate in the outer circumference of the core and dissolution of that methane into the fresh pore water that the core is stored with. We are designing approaches to minimize this core loss in the future. We spent an enormous amount of effort to further improve the ability of the pressure coring tool (the PCTB) to pressure seal correctly. We completed upgrading the upper section of the PCTB to address poor pressure. We successful tested the modifications at Geotek’s test facility in Salt Lake City (Bench Test II). We completed a Land Test of the PCTB at the Schlumberger Cameron Test and Training Facility (CTTF). The tool did not seal in 6 out of 7 tests and we clearly demonstrated that cuttings were wedging in the ball valve assembly, keeping the ball valve from sealing. We reproduced the failure mechanism observed during the land test at Salt Lake City and confirmed the sensitivity of the ball valve assembly to grit. Geotek designed and tested 9 modifications to address this issue and the PCTB is now 100% successfully sealing in the presence of grit. Our science expedition is scheduled for spring 2022 and we are fully focused on preparing for this. UT and Ohio State completed a Shallow Hazard Assessment report for each proposed UT-GOM2-2 drilling location, pursuant to 30 CFR 250.214(f) and 250.244 (f). The Shallow Hazard Reports will accompany the UT-GOM2-2 Exploration Plan that is submitted to BOEM, and completes the geological and geophysical analysis for UT-GOM2-2 permitting efforts. We updated the UT-GOM2-2 Operations Plan (Version 1). We completed the UT-GOM2-2 Science and Sample Distribution Plan (Version 1). We evaluated the scope, budget, and schedule that would result from using a commercial vessel. We developed detailed drilling schedule, mud volume, and resource estimates. We developed a vessel specification document and a well plan, and sent these documents to prospective vessel contractors.

03 NATURAL GAS↗

Experimental study of the partitioning of some platinum group elements (Pd and Ir) between orthopyroxene and silicate melt

Past experiments and observations on natural samples have largely focused on the roles of olivine and chromite in controlling the behaviour of the platinum-group elements (PGE) during melting and solidification, whereas other phases, such as pyroxene, have gone largely uncharacterized. Here, to address this, experiments have been done to measure the partitioning of Pd (with a subset of results for Ir), between orthopyroxene and silicate melt at 1340 °C, 0.1 MPa and log fO 2 of FMQ - 1 to FMQ + 6 (FMQ = Fayalite-Magnetite-Quartz). The X-ray Absorption Near-Edge Structure (XANES) was measured in a subset of experiment glasses. Glass concentrations of Pd (corrected to unit Pd activity) increase from ~6 to ~650 ug/g over the fO 2 range of experiments. The slope of the solubility-fO 2 relation is consistent with Pd 1+ as the dominant oxidation state, with evidence for Pd 0 and Pd 2+ at the lowest and highest experiment fO 2 , respectively. Consistent with this result, the XANES reveal spectral features similar to Pd 0 and Pd 2+ spectral reference materials (specRM) at the most reduced and oxidized synthesis conditions, respectively. Other lines of evidence require the presence of a third melt species, here interpreted to be Pd 1+ . Values of orthopyroxene/melt partition coefficients for Pd (D Pd Opx/melt ) are 0.0051 (+/-0.006) at log fO 2 < ΔFMQ + 3, increasing with fO 2 to a maximum of 0.013 at ~FMQ + 6. Sodium partition coefficients, expected to be similar to Pd, range from 0.0061 (+/-0.00061) at FMQ + 3, increase to 0.007–0.009 at higher fO 2 , but with no clear systematic trend. A value for D Ir Opx-melt of ~0.6 was measured at ~FMQ + 4, indicating significantly more compatible behaviour for Ir relative to Pd. Partitioning results are interpreted in the context of the Blundy-Wood elastic strain model in which the variation in partitioning is related to ionic radius mismatch to an optimal crystallographic site size. Based on the trend in ionic radius with oxidation state, the estimated ionic radius of Pd 1+ in octahedral coordination is similar to Na 1+ , and comparison to previous orthopyroxene-melt partitioning experiments suggests D Pd1+ opx/melt and DNaopx/melt should be nearly identical, consistent with the results of this study. The ionic radius of VI-fold Pd 2+ is close to the optimal M2 site size, so an increased proportion of this species with fO 2 accounts for the larger values of D Pd opx/melt at the highest fO 2 investigated. The much larger partition coefficient for Ir is consistent with the presence of Ir 2+ , whose estimated ionic radius is close to Fe 2+ and Mg 2+ , as well as predictions for the optimal M1 site size. With the assumption that D Pd opx/melt = D Na opx/melt , combined with a revised value for the Pd content of the primitive mantle, a melting model is presented that better reproduces the Pd concentration of high degree melts from sulfide-free mantle sources.

58 GEOSCIENCES↗

Sulfur speciation in dacitic melts using X-ray absorption near-edge structure spectroscopy of the S K -edge (S-XANES): Consideration of radiation-induced changes and the implications for sulfur in natural arc systems

Abstract The synchrotron technique of micro X-ray absorption near-edge structure spectroscopy at the sulfur K-edge (S-XANES) provides a unique opportunity to measure the proportion of different oxidation states of sulfur (S) in silicate glasses. Although applied extensively in the analysis of basaltic silicate glasses, few S-XANES studies have investigated variations in S oxidation states with fO2 in felsic silicate glasses. In addition, no study has systematically compared the S-XANES results obtained from the same samples at different photon flux densities to quantify the relationship between exposure time and changes in S speciation in silicate glass, as has been done for Fe and V. This study evaluates observed differences in S speciation measured in experimentally produced H2O-saturated dacitic glasses over a range of reducing to oxidizing conditions (from log fO2 = ΔFMQ-0.7 to ΔFMQ+3.3; FMQ is the fayalite-magnetite-quartz mineral redox buffer) and equilibrated at 1000 °C and 300 MPa. S-XANES spectra were collected at three different photon flux densities using three microspectroscopy beamlines. As is observed in S-XANES analyses of basaltic silicate glasses, beam-induced changes to the S6+/ΣS are observed as a function of photon flux density and beam exposure time. Our results demonstrate that silicate glasses of dacitic composition undergo beam-induced photo-reduction in samples equilibrated at ΔFMQ > +1.75 and photo-oxidation if equilibrated at ΔFMQ < +1. The time required to observe beam-induced changes in the spectra varies as a function of flux density, and our study establishes an upper photon density limit at ~1.0 × 1012 photons/µm2. The S6+/ΣS calculated from spectra collected below this absorbed photon limit at intermediate flux densities (~1–4 × 109 photons/s per µm2) are affected by beam damage, as no conditions were found to be completely free of beam-induced changes. However, the S6+/ΣS ratios calculated below the limit at intermediate flux densities are consistent with thermodynamic constraints, demonstrating that S6+/ΣS ratios calculated from S-XANES spectra can be considered reliable for estimating the oxygen fugacity. Our results carry important implications for the S budget of felsic magmas and dissolution mechanisms in evolved melts. While our results from all three flux densities show the presence of S4+ dissolved in relatively oxidized (ΔFMQ > +1.75) dacitic glass, even in the spectra exposed to the lowest photon densities, we are unable to rule out the possibility that the S4+ signal is the result of instantaneous X-ray irradiation induced beam damage using S-XANES alone. When our spectra are compared to S-XANES spectra from basaltic silicate glasses, important differences exist in the solubility of S2– and S6+ between dacitic silicate melts, pointing to differences in solubility mechanisms as melt composition changes. This study highlights the need for further investigation into beam damage systematics, presence of S4+, and the solubility mechanisms of different oxidation states of S as silicate melt composition changes.

Geochemistry & Geophysics↗