Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical reservoir models”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Numerical study of the effect of propped surface area and fracture conductivity on shale gas production: Application for multi-size proppant pumping schedule design

Hydraulic fracturing is a technique extensively used in the oil and gas industry, where water, proppant (sand) and additives are injected into unconventional reservoirs to enhance the recovery of shale hydrocarbon. Although some previous studies have developed pumping schedules that maximize gas production for a single-size proppant, there are very few studies that consider the effect of varying proppant diameters across pumping stages on shale gas production. Motivated by this, we carried out an extensive sensitivity analysis to determine the effect of different proppant diameters on the average fracture conductivity (FC), average propped surface area (PSA) and cumulative shale gas production volume. We found out that the cumulative shale gas production volume depends on both the average PSA and average FC. We also found out that small-diameter proppant resulted in higher average PSA and lower average FC, whereas large-diameter proppant resulted in lower average PSA and higher average FC. Hence, we designed a multi-size proppant pumping schedule considering both of these parameters into account for simultaneously propagating multiple fractures to maximize shale gas production from unconventional reservoirs. Since the size of injected proppant particles determines the average PSA and average FC for the propped hydraulic fractures, we developed a novel framework called Sequentially Interlinked Modeling Structure (SIMS) to predict the average PSA, average FC and cumulative shale gas production volume at the end of 10 years for a given pumping schedule. Then, we used this SIMS framework to obtain a multi-size proppant pumping schedule that maximizes shale gas production. In closing, we demonstrated that the obtained pumping schedule gives a cumulative shale gas production volume greater than the values obtained from the existing pumping schedules.

03 NATURAL GAS↗
DOE OSTI2023

Field Results from New Tensor Borehole Optical Fiber Strainmeter Installations in Oklahoma and Utah

The time evolving strain field contains a wealth of information that can be used to interpret subsurface behavior. For example, injecting or removing fluids from reservoirs or aquifers causes deformation that can be used as a diagnostic signal in some cases, while it can interfere with geodetic interpretations in other cases. We've previously demonstrated the feasibility of measuring the strain tensor at a depth of 30m caused by injection into a reservoir at 530m. The observed strain signals were interpreted using four independent analytic and numerical methods that resulted in estimates of the poroelastic properties and geometry of the reservoir that was consistent with data from well logs. However, studies like these are only possible if these deformations can be reliably measured. Years of lab and field work has culminated in the development of a novel borehole strainmeter capable of resolving multiple components of strain using embedded optical fibers configured as Michelson interferometers. It features four horizontal gauges separated by 45° to resolve the horizontal strain tensor as well as a vertical strain gauge and a sixth null component for state-of-health monitoring. The downhole sensing package also includes an open pipe through its center for grout circulation during single-trip deployments and a fully welded stainless steel exterior for robustness and longevity. These instruments have a resolution of 2x10-13 strain that can easily measure the solid earth tides. Preliminary data are available from four strainmeters in shale at our Oklahoma site and four in compacted sand and gravel in Utah. These are deployed from 40-60m, except one of the strainmeters in Oklahoma is deployed at 500m. The data include strains from the initial grout curing, comparisons to predicted earth tide models and in-situ calibration results, barometric pressure admittances and spectral analyses as well as signals from underground injections and surface waves from teleseismic events. Preliminary analyses indicate behavior consistent with other strainmeter deployments, and comparison to data from a Gladwin strainmeter at the Oklahoma site validate the performance of the new design. Analyses from a suite of six well tests at the Oklahoma site show for the first time how the strain tensor field varies with location during well testing.

DeWolf, Scott↗
DOE OSTI2020

Gas Production from Hot Water Circulation through Hydraulic Fractures in Methane Hydrate-Bearing Sediments: THC-Coupled Simulation of Production Mechanisms

Methane hydrates, widely found in permafrost and deep marine sediments, have great potential as a future energy 6 source. Conventional production schemes perform poorly for challenging hydrate reservoirs with low permeability. We propose an 7 efficient production scheme by combining hydraulic fracturing from horizontal wells and hot water circulation through fractures. A 8 fully coupled thermo-hydro-chemical (THC) model is developed to simulate the key physical processes during gas production from 9 a hydrate reservoir representative of typical geological settings in Shenhu, South China Sea. We found that the gas production 10 process has two distinct stages divided by thermal breakthrough: a relatively short prebreakthrough stage and a postbreakthrough 11 stage yielding stable gas production. Heat advection along and near the hydraulic fracture dominates the prebreakthrough stage, 12 whereas conduction-driven thermal recovery in the volume around fractures dominates the postbreakthrough stage. We identified 13 that the steady-state injection temperature has a strong effect on the performance of the proposed scheme while the fluid mass 14 circulation rate has a moderate impact beyond a threshold. The proposed scheme proves to be efficient and robust over a range of 15 reservoir conditions with respect to initial hydrate saturation and intrinsic permeability, including their spatial heterogeneities, 16 thereby offering a promising solution for challenging reservoir conditions.

gas hydrate, fracture flow, gas production, numeri↗
DOE OSTI2023

The α-WC(0001) Surface as a Hydrogen Sponge: A First Principle Study of H 2 Dissociation and Formation of Low and High Coverages

Tungsten carbide (WC) displays a Pt-like behavior in catalysis, applied in hydrogenation processes. Numerous theoretical studies have modeled the behavior and use of adsorbed hydrogen without obtaining a general picture, missing basic links between H 2 dissociation and generation of high surface coverage (Θ H >0.5 ML). Here, in this study, the capability of C- and W-terminations of the α-WC(0001) surface is analyzed to dissociate several H 2 molecules to produce coverages, Θ H , ranging from low to very high values (0.13<Θ H <2.00 ML). Density functional theory and an ab initio atomistic thermodynamic were used to achieve the conditions for H 2 dissociation. The WC–C surface has higher capacity to dissociate H 2 molecules than WC–W. However, both surfaces can reach full surface coverage, Θ H =1 ML, at mild ambient conditions, T=300 K and P=1 atm, and even up to 500 K at low and high pressures. The H-adatoms on WC–W are more labile than on WC–C. The binding of adsorbates is hindered at high Θ H , implying a need to modulate Θ H according to the application. The results give the basis to understand the capabilities of WC-based catalysts in hydrogenation-related reactions, with the advantage of WC being a hydrogen reservoir at mild practical catalytic conditions.

03 NATURAL GAS↗
DOE OSTI2024

Reactive transport modeling of the Aquifer Thermal Energy Storage (ATES) system at Stockton University, New Jersey during seasonal operations

Hydrogeochemical processes associated with Aquifer Thermal Energy Storage (ATES) operations can often impact the system performance owing to mineral precipitation either at the wellbore or in the aquifer owing to changes in temperature and fluid disequilibria. Although failure of ATES systems due to mineral precipitation ("fouling") is common, predictive reactive-transport models have rarely been applied to plan their design and operation. Here, the objective of this study is to develop a reactive-transport model by coupling thermal, hydrological, and chemical (THC) processes to evaluate effects of introduced atmospheric oxygen on water chemistry, mineral precipitation/dissolution, porosity, and permeability changes associated with an ATES system at Stockton University (New Jersey, USA). The THC model builds on a Thermal-Hydrological-Mechanical (THM) model of the site that evaluated system failure owing to possible fracturing in the caprock or around the wellbore. The causes of the system failure are not known – potential causes include hydraulic fracturing owing to elevated pump pressures that took place, a flow pathway created by one of the boreholes, or a pre-existing natural hydrologic connection between the upper unconfined aquifer and the ATES aquifer, any of which could have led to oxygenated water entering the reservoir and causing the observed Fe-oxide fouling on well screens. The THC model is used to evaluate some of the hypotheses and observations regarding system failure owing to geochemical processes. The reactive-transport code TOUGHREACT V4 was used to model the THC processes during seasonal heating and cooling operations at the Stockton ATES site over 6 years of operation. In the THC simulations, the primary effects on geochemistry were observed when the injection water is saturated with atmospheric oxygen. Simulations show greater precipitation of goethite near the cold wells as compared to the warm wells. Although volume fractions of Fe-hydroxides were relatively small, the model was aimed at processes in the aquifer at the scale of meters and larger rather than at the scale of mm or cm (i.e., a well screen). Kaolinite is the dominant precipitating phase, also around the cold wells. Illite dissolves near the cold wells and precipitates near the warm wells. There is a net decrease in the porosity near the cold wells and increase near the warm wells, although a slight amount of thermal contraction near the cold wells and expansion near the warm wells is responsible for a significant proportion of the porosity change. Owing to the coarse discretization of the numerical grid near the wells (compared to the screen thickness) the magnitude of permeability changes at the wellbore are likely underestimated. The reactive transport model in this study can be used for characterization of aquifers, optimizing the operational parameters (temperature, pressure, pH etc.), and planning of mitigation strategies for ATES systems.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Efficient and generalizable nested Fourier-DeepONet for three-dimensional geological carbon sequestration

Geological carbon sequestration (GCS) involves injecting CO2 into subsurface geological formationsfor permanent storage. Numerical simulations could guide decisions in GCS projects by predictingCO 2 migration pathways and the pressure distribution in storage formation. However, these simula-tions are often computationally expensive due to highly coupled physics and large spatial-temporalsimulation domains. Surrogate modelling with data-driven machine learning has become a promis-ing alternative to accelerate physics-based simulations. Among these, the Fourier neural operator(FNO) has been applied to three-dimensional synthetic subsurface models. Despite its good accuracyin simulating CO 2 plume migration, it requires large computational resources in training and alsolacks generalizability. Here, to further improve performance, we have developed a nested Fourier-DeepONet by combining the expressiveness of the FNO with the modularity of a deep operatornetwork (DeepONet). This new framework is twice as efficient as a nested FNO for training and has atleast 80% lower GPU memory requirement due to its flexibility to treat temporal coordinates sepa-rately. These performance improvements are achieved without compromising prediction accuracy.In addition, the generalization and extrapolation ability of nested Fourier-DeepONet beyond thetraining range has been thoroughly evaluated. Nested Fourier-DeepONet outperformed the nestedFNO for extrapolation in time with more than 50% reduced error. It also exhibited good extrapolationaccuracy beyond the training range in terms of reservoir properties, number of wells, and injectionrate.

Lee, Jonathan E. [Department of Chemical and Envir↗
DOE OSTI2022

Sodium diffusion in heterogeneous porous media: Connecting laboratory experiments and simulations

Sodium has been suggested as a tracer for brine in reservoir formations where a significant amount of sodium ion is found mainly in the aqueous phase. Fortunately, the most abundant sodium isotope, 23 Na, is Nuclear Magnetic Resonance (NMR) active, making it possible to study the structure and dynamical properties of both bulk and pore confined fluid systems. We investigated the diffusion of key dissolved sodium species in bulk solution and porous media as a function of composition, ionic strength, and permeability of the carbonate samples using NMR relaxometry and pulsed-field gradient (PFG) NMR techniques. We use Indiana limestone as an example of natural porous media and water as a freely diffusible tracer and carrier. We demonstrate apparent diffusion measurement of sodium based on changes in spin–spin relaxation time (T 2 ) signal in pore confinement. The diffusion rate of sodium decreases with increasing counter-ion size. This effect is greater at higher ionic strengths and lower chemical potentials in porous media. The reactive transport code, CrunchFlow, was used to complement the NMR experiments to simulate diffusion behavior in porous media. Combining these two methods provides a powerful approach to estimating effective diffusion coefficients in heterogeneous matrices. The modeling considers the influence of physical properties (porosity and tortuosity) and chemical properties (geochemical composition and chemical gradients). The results presented in this work highlight the advantage of measuring apparent diffusivity using NMR T 2 relaxometry in conjunction with numerical simulation to derive effective diffusivity and the corresponding matrix properties (i.e., tortuosity) of the system.

58 GEOSCIENCES↗
DOE OSTI2026

Pore-Scale Study on the Positive Feedback Between Stress and Porosity Caused by Pressure Solution in Porous Media

Pressure solution is an important process in the evolution of sedimentary rocks, which provide storage space for most of our petroleum resources. It directly influences the generation, migration, and storage of petroleum fluids in subsurface sedimentary rocks. Here, in this paper, we develop a pore-scale, mechanochemical model to demonstrate a possible positive feedback between the local porosity and pore surface stress, in which a higher local porosity causes a higher local pore surface stress, thus enhancing pressure solution and consequently further increasing the local porosity. Pore surface stress represents stress on a solid grain adjacent to a pore. Specifically, the pore-scale, mechanochemical model directly simulates the stress distribution over solid and pore surfaces using a finite element model. The dissolution of solids at the solid-pore interfaces under a far-from-equilibrium condition is simulated using a first-order kinetics model that accounts for the local stress distribution. The updated pore geometry, caused by pore surface dissolution, is then used in the stress simulation in the next numerical iteration. Two types of porous media, the Oriskany sandstone and an artificial porous medium with spherical pores, were tested in the mechanochemical simulation. The positive stress-porosity feedback during pressure solution was observed in both samples. In addition, the model quantitatively illustrated the distribution of local mineral dissolution rates on all pore surfaces, as well as its relation to the effective mineral dissolution rate of the entire sample. Based on the comparison between the two porous media, the local mineral dissolution was regulated by pore space distribution, geometry, and coalescence during pressure solution. This work is the first that uses direct, pore-scale numerical simulation to demonstrate the positive stress-porosity feedback during pressure solution, which has the potential to advance the understanding of the mechanical-chemical (MC) coupling in many geological processes that are relevant to subsurface energy systems, such as the recovery of petroleum hydrocarbons and geothermal energy.

CT scanning↗
DOE OSTI2024

Quarterly Research Performance Progress Report (Q8)

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. Once satisfied with the results we have started by targeting the modeling of the DFITs at conducted at well 58-32. We have a identified a specific test (cycle 4 in zone 2) as the most interesting to be model with GEOS hydraulic fracturing module. Thus, we have first produced results with an isothermal model and adjusted model parameters to get a satisfying match with field pressure data. The, we have added thermal effects and compared the modeling results with and without thermal effects to estimate how thermal effects may influence test interpretation. Models seem to suggest that, for small volumes of fluid, thermal effects are moderate. In Task 2, we have adapted GEOS phase-field formulation to be able to simulate near-wellbore hydraulic fracture nucleation and propagation. We have devised a novel formulation that, compared to other existing ones, incorporates rock strengths. We have submitted a journal publication about our work. We are currently employing this phase-field formulation to model the experiments taking place at U Pitt and help us understand the effect of various parameters. In Task 3, we have built a model of the region surround well 16A and have started modeling stage 3 stimulation because of its simpler planar geometry. After calibrating simulation parameters using known analytical solutions, we have simulated the stage 3 stimulation using our isothermal hydraulic fracturing module, varying the permeability field, the stress conditions including different physics to get a better understanding of the numerical challenges and of the effects of varying these parameters on the simulation results. In Task 4 laboratory experimentation, a set of specialized drilling and injection tools has been customized and constructed to accommodate an inclined well with an orientation of up to 30 degrees relative to material anisotropy or principal stress axes. These inclined samples have also undergone thermal stress and hydraulic fracturing at a temperature of 190 degrees Celsius. Furthermore, both vertical and deviated sampling testing setups enable an extended analysis of post-peak pressure behaviors, facilitating post-test pressure analyses such as the G-function, step rate, and fracture reopening measurements. Thus, the key components of in-situ stress estimation can be extracted and validated through our experiment, providing a solid foundation for validating existing in-situ stress estimation theories or proposing new ones. Simultaneously, we are integrating computer vision techniques with traditional experimental fracture observation methods such as multi-overcore/slicing and water-penetration fracture observation. This combination will prove beneficial in populating the hydraulic fracture patterns database, generated under challenging EGS conditions. This approach aims to deepen our understanding of the complexities in EGS reservoirs and pave the way for future data-driven investigations. Additionally, PITT has also equipped the ELE International compression machine, which is now prepared for conducting indirect tensile and fracture toughness tests. These tests will aid in characterizing how rock fabrics influence the resulting fracture patterns. Additionally, we have completed the required personnel training and gained access to Scanning Electron Microscopy (SEM) and Energy Dispersive Spectroscopy (EDS) for conducting more detailed characterization and analysis of rock fabrics, as well as the examination of thermal and hydraulically induced fracture patterns. Thus, the PITT team has effectively demonstrated the capabilities of our experimental apparatuses in exploring the thermal effects, well deviation angles, material anisotropy, and operational choices (such as circulation rate, injection fluid viscosity, and injection rate) and their impact on pressure responses and fracture trajectories under the Utah FORGE conditions.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Completion Effects on Diagnosing Multistage Facture Treatments with Distributed Temperature Sensing

Summary Distributed temperature sensing (DTS) is a valuable tool to diagnose multistage hydraulic fracture treatments. When a stage interval is shut in, the clusters that take more fluid during pumping warm up more slowly. Therefore, the fluid volume injected into each cluster can be quantitatively interpreted by numerical inversion of the warm-back temperature behavior. This general concept assumes that the different warm-back behavior is controlled by only the injected fluid volume; however, recent observations of DTS data indicate that completion configurations significantly influence the warm-back behavior. This paper investigates the completion effects on the DTS interpretation. In ideal conditions, when a stage is fractured, the upstream stage intervals should show an almost uniform temperature that is close to the injected fluid temperature. This is due to the high fluid velocity of injected fluid in the wellbore, and the upstream intervals have not been perforated (noncommunicating intervals), so the only heat transfer is heat conduction between the wellbore fluid and the surrounding reservoir. But the field DTS data show considerably irregular variations in temperature along the upstream stage intervals. These variations are caused by the completion effects. The nonuniform temperature profile is caused by different heat transfer behavior induced by completion hardware along the production casing string, such as joints, clamps, and blast protectors, and by the sensing cable location in the cement, as well as the cement quality. Because the varying heat transfer behavior impacts the warm-back behavior as well as the temperature profile, the completion effects need to be considered in DTS interpretation. A method of DTS interpretation considering the completion effects to diagnose multistage fracture treatments was developed. Because the heat transfer between a wellbore and a reservoir depends on the overall heat transfer coefficient describing heat conduction through the completion in a forward model, this parameter needs to be tuned all along the wellbore. To calibrate the completion effect, the temperature inversion is conducted using the temperature measured at a stage interval that is upstream of a stage interval currently being treated. Because the interpreted stage interval is not perforated at that time, the thermal behavior at the noncommunicating interval is governed by only the heat conduction through the completion environment. Once the effective values of the overall heat transfer coefficient are estimated along the interpreted stage interval, they can be assumed to be constant physical parameters. Then, the fluid volume distribution is interpreted by using the effective overall heat transfer coefficient profile along each interval. This study provides a field application of the developed interpretation method. The new interpretation method provides more accurate diagnosis of fracture treatments by DTS interpretation.

Engineering↗
DOE OSTI2023

Transient Pressure Interference During CO2 Injection in Saline Aquifers

Abstract CO2 injection in subsurface geological formations (e.g. deep saline aquifers) causes pressure perturbations over a large area surrounding the injection well. Observation wells are widely considered in geologic CO2 storage (GCS) projects where the pressure perturbation induced by CO2 injection is measured. In this work, we use analytical and numerical modeling tools along with field data to examine the pressure behavior in GCS projects before and after CO2 arrival at an observation well. Prior to CO2 arrival, a baseline pressure trend is established which corresponds to single-phase brine flow across the observation well (approximated by Theis solution). Therefore, analysis of early-time pressure data is straightforward, provides the single-phase flow characteristics (mobility and storativity), and helps establishing a baseline pressure change that can be extended beyond the single-phase flow period at the observation well. Upon CO2 arrival, a departure from this baseline trend is expected. For the pressure to detect the CO2 arrival at an observation well, the departure from baseline pressure behavior must be significant and well above the background noise levels. We use existing analytical models to determine the strength of the expected pressure departure signal from the baseline trend upon CO2 arrival. The strength of the expected pressure departure is found to be directly proportional to the mobility ratio. Accordingly, we establish a criterion to determine whether the pressure at an observation well can detect the CO2 arrival. We present an analysis approach through application to synthetic and field data and show the characteristic pressure behavior before and after CO2 arrival. We show that while generally the pressure can be either above or below the expected baseline pressure trend, it would be likely above the baseline upon CO2 arrival. This is because the mobility ratio becomes less than unity after CO2 arrival. We show that depending on the reservoir characteristics, changes in the pressure trend may or may not be sufficient to detect the CO2 arrival.

Engineering↗
DOE OSTI2024

SQMS Quantum R&D in Machine Learning, Optimization and Sensing beyond Fundamental Physics Applications

This newly formed team at SQMS under the Ecosystem Thrust is looking to develop capabilities impacting societal advances outside the core domain of HEP and condensed matter physics. We explicitly leverage the experimental and algorithmic innovations developed across all groups as well as connect to broad-scope external projects of the diverse team of PIs. As the inaugural set of projects, we are studying numerically quantum machine learning models inspired by efficiently trainable echo-state and orthogonal neural networks and developing designs for related experiments to be performed on quantum processors based on SQMS SRF cQED technology and Rigetti s transmon arrays. Investigated models exploit ideas and lessons learned from multiple prior work by SQMS team members in a variety of internal and external activities [R1]. Target initial applications include noisy signal processing, potentially captured by quantum sensors or noisy QPUs, as well as simulation and classification of healthcare data. For instance, image reconstruction of the brain s electrical properties by solving the inverse Maxwell equation problem with uncertainty [R2] through a hybrid quantum-classical physics-informed architecture for time-dependent processes [R3]. The group is also investigating the application and development of novel quantum sensors based on magnetic levitation of a superconducting sphere coupled to a superconducting qubit. This coupling enables high-precision measurements of the position of the sphere, which can be used for sensitive detection of forces, enabling practical applications such as gravimetry for geophysics analysis, or accelerometry for GPS-denied navigation [R4] [R1] Rieffel, Eleanor G., Ata Akbari Asanjan, M. Sohaib Alam, Namit Anand, David E. Bernal Neira, Sophie Block, Lucas T. Brady et al. "Assessing and advancing the potential of quantum computing: A NASA case study." Future Generation Computer Systems (2024). [R2] Yu, X., Serrall s, J.E., Giannakopoulos, I.I., Liu, Z., Daniel, L., Lattanzi, R. and Zhang, Z., 2023. Pifon-ept: Mr-based electrical property tomography using physics-informed fourier networks. IEEE Journal on Multiscale and Multiphysics Computational Techniques. [R3] Wudarski, Filip, Daniel OConnor, Shaun Geaney, Ata Akbari Asanjan, Max Wilson, Elena Strbac, P. Aaron Lott, and Davide Venturelli. "Hybrid quantum-classical reservoir computing for simulating chaotic systems." arXiv preprint arXiv:2311.14105 (2023). [R4] Higgins, Gerard, Saarik Kalia, and Zhen Liu. "Maglev for dark matter: Dark-photon and axion dark matter sensing with levitated superconductors." Physical Review D 109.5 (2024): 055024.

Venturelli, Davide↗
DOE OSTI2024

Dynamic modeling studies of basin-scale pressure interference and CO 2 plume evolution in multi-well geologic CO 2 storage

This research employs numerical modeling to analyze how CO 2 plumes and pressure fronts evolve when CO 2 is injected into a single storage formation from multiple projects located in proximity. Here, this analysis also evaluates the extent to which injection well spacing alleviates pressure buildup in the absence of active pressure management tactics. The simulation approach was based on a single, homogenous saline aquifer in which CO 2 injection occurs under a one-injector baseline case and several multi-well cases where well spacing varies. Analysis results show that the extent of pressure buildup is in the range of tens or a few hundreds of kilometers and contingent upon the defining pressure buildup demarcating the front edge. For the geological setting evaluated in this paper, our analysis suggests that without active basin pressure management strategies, commercial-scale projects would likely need to be sited far apart to avoid pressure interference from one another. Analysis results show the radius of CO 2 plume varies approximately from 2 to 3 km from injection wells (each injecting 1 Mt/year for 30 years) depending on cases and modeling parameters assumed. Given the pressure interference, this paper thus draws attention to the importance of greater coordination among storage operators and regulatory stakeholders. Because this analysis assumes a very specific geologic setting, this exploratory analysis bears further investigations across other geologic scenarios.

42 ENGINEERING↗
DOE OSTI2025

Modeling gas, hydrates, and slope stability on the U.S. Atlantic margin during Pleistocene glacial cycles

Changes in temperature and sea level can cause dissociation of methane hydrates in shallow marine sediments, leading to seafloor destabilization. Along the U.S. Atlantic margin, there exists a well-documented history of slope failure and numerous recorded occurrences of gas seeps. Several studies have linked slope failure in the region to gas seepage and hydrate dissociation driven by glacial-interglacial transitions, but this linkage has not been quantitatively demonstrated. Along the shelf edge, in an area where shallow methane gas seeps have been identified, we modeled methane gas and hydrate formation using a one-dimensional fluid flow model. Methane gas formation was modeled over the last 120,000 years to simulate a glacial-interglacial cycle. Here, we ran this model at 16,044 individual locations in the region between 29° N – 45° N and 82° W – 66° W at a resolution of 1 x 1 arcminutes, focusing specifically on water depths between 200 and 1000 m that bracket the seafloor outcrop of the base of the hydrate stability zone. Using historic temperature and pressure records from the last 120,000 years, sediment properties in the area, and factor of safety calculations, we found that hydrate dissociation alone is unlikely to cause sediment failure in the region, implying that an additional driving force would be necessary for failure to occur.

Hydrate dissociation↗
DOE OSTI2021

Radial structure and formation of the Milky Way disc

The formation of the Galactic disc is an enthusiastically debated issue. Numerous studies and models seek to identify the dominant physical process(es) that shaped its observed properties; for example, satellite accretion, starburst, quenching, gas infall, and stellar radial migration. Taking advantage of the improved coverage of the inner Milky Way provided by the SDSS DR16 APOGEE catalogue and of the ages published in the APOGEE-AstroNN Value Added Catalogue (VAC), we examined the radial evolution of the chemical and age properties of the Galactic stellar disc with the aim of better constraining its formation. Using a sample of 199 307 giant stars with precise APOGEE abundances and APOGEE-AstroNN ages, selected in a ±2 kpc layer around the galactic plane, we assessed the dependency as a function of guiding radius of (i) the median metallicity, (ii) the ridge lines of the [Fe/H] – [Mg/Fe] and age–[Mg/Fe] distributions, and (iii) the age distribution function (ADF). The giant star sample allows us to probe the radial behaviour of the Galactic disc from R g = 0 to 14–16 kpc. The thick disc [Fe/H] – [Mg/Fe] ridge lines follow closely grouped parallel paths, supporting the idea that the thick disc did form from a well-mixed medium. However, the ridge lines present a small drift in [Mg/Fe], which decreases with increasing guiding radius. At sub-solar metallicity, the intermediate and outer thin disc [Fe/H] – [Mg/Fe] ridge lines follow parallel sequences shifted to lower metallicity as the guiding radius increases. We interpret this pattern as the signature of a dilution of the interstellar medium from R g ~ 6 kpc to the outskirts of the disc, which occurred before the onset of the thin disc formation. The APOGEE-AstroNN VAC provides stellar ages for statistically significant samples of thin disc stars from the Galactic centre up to R g ~ 14 kpc. An important result provided by this dataset is that the thin disc presents evidence of an inside-out formation up to R g ~ 10 – 12 kpc. Moreover, about ~7 Gyr ago, the [Mg/Fe] ratio in the outer thin disc (R g > 10 kpc) was higher by about ~0.03–0.05 dex than in the more internal regions of the thin disc. This could be the fossil record of a pollution of the outer disc gas reservoir by the thick disc during its starburst phase.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Robust Carbon Dioxide Plume Imaging Using Joint Tomographic Inversion of Seismic Onset Time and Distributed Pressure and Temperature Measurements (Final Report)

We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.

02 PETROLEUM↗
DOE OSTI2024

HarDWR - Harmonized Water Rights Records

A dataset within the Harmonized Database of Western U.S. Water Rights (HarDWR). For a detailed description of the database, please see the meta-record v2.0. Changelog v2.0 - Recalculated based on data sourced from WestDAAT - Changed using a Site ID column to identify unique records to using aa combination of Site ID and Allocation ID - Removed the Water Management Area (WMA) column from the harmonized records. The replacement is a separate file which stores the relationship between allocations and WMAs. This allows for allocations to contribute to water right amounts to multiple WMAs during the subsequent cumulative process. - Added a column describing a water rights legal status - Added "Unspecified" was a water source category - Added an acre-foot (AF) column - Added a column for the classification of the right's owner v1.02 - Added a .RData file to the dataset as a convenience for anyone exploring our code. This is an internal file, and the one referenced in analysis scripts as the data objects are already in R data objects. v1.01 - Updated the names of each file with an ID number less than 3 digits to include leading 0s v1.0 - Initial public release Description Here we present an updated database of Western U.S. water right records. This database provides consistent unique identifiers for each water right record, and a consistent categorization scheme that puts each water right record into one of seven broad use categories. These data were instrumental in conducting a study of the multi-sector dynamics of inter-sectoral water allocation changes though water markets (Grogan et al., *in review*). Specifically, the data were formatted for use as input to a process-based hydrologic model, Water Balance Model (WBM), with a water rights module (Grogan et al., *in review*). While this specific study motivated the development of the database presented here, water management in the U.S. West is a rich area of study (e.g., Anderson and Woosly, 2005; Tidwell, 2014; Null and Prudencio, 2016; Carney et al., 2021) so releasing this database publicly with documentation and usage notes will enable other researchers to do further work on water management in the U.S. West. We produced the water rights database presented here in four main steps: (1) data collection, (2) data quality control, (3) data harmonization, and (4) generation of cumulative water rights curves. Each of steps (1)-(3) had to be completed in order to produce (4), the final product that was used in the modeling exercise in Grogan et al. (*in review*). All data in each step is associated with a spatial unit called a Water Management Area (WMA), which is the unit of water right administration utilized by the state in which the right came from. Steps (2) and (3) required use to make assumptions and interpretation, and to remove records from the raw data collection. We describe each of these assumptions and interpretations below so that other researchers can choose to implement alternative assumptions an interpretation as fits their research aims. Motivation for Changing Data Sources The most significant change has been a switch from collecting the raw water rights directly from each state to using the water rights records presented in WestDAAT, a product of the Water Data Exchange (WaDE) Program under the Western States Water Council (WSWC). One of the main reasons for this is that each state of interest is a member of the WSWC, meaning that WaDE is partially funded by these states, as well as many universities. As WestDAAT is also a database with consistent categorization, it has allowed us to spend less time on data collection and quality control and more time on answering research questions. This has included records from water right sources we had previously not known about when creating v1.0 of this database. The only major downside to utilizing the WestDAAT records as our raw data is that further updates are tied to when WestDAAT is updated, as some states update their public water right records daily. However, as our focus is on cumulative water amounts at the regional scale, it is unlikely most records updates would have a significant effect on our results. The structure of WestDAAT led to several important changes to how HarWR is formatted. The most significant change is that WaDE has calculated a field known as `SiteUUID`, which is a unique identifier for the Point of Diversion (POD), or where the water is drawn from. This separate from `AllocationNativeID`, which is the identifier for the allocation of water, or the amount of water associated with the water right. It should be noted that it is possible for a single site to have multiple allocations associated with it and for an allocation to be able to be extracted from multiple sites. The site-allocation structure has allowed us to adapt a more consistent, and hopefully more realistic, approach in organizing the water right records than we had with HarDWR v1.0. This was incredibly helpful as the raw data from many states had multiple water uses within a single field within a single row of their raw data, and it was not always clear if the first water use was the most important, or simply first alphabetically. WestDAAT has already addressed this data quality issue. Furthermore, with v1.0, when there were multiple records with the same water right ID, we selected the largest volume or flow amount and disregarded the rest. As WestDAAT was already a common structure for disparate data formats, we were better able to identify sites with multiple allocations and, perhaps more importantly, allocations with multiple sites. This is particularly helpful when an allocation has sites which cross WMA boundaries, instead of just assigning the full water amount to a single WMA we are now able to divide the amount of water between the number of relevant WMAs. As it is now possible to identify allocations with water used in multiple WMAs, it is no longer practical to store this information within a single column. Instead the stAllocationToWMATab.csv file was created, which is an allocation by WMA matrix containing the percent Place of Use area overlap with each WMA. We then use this percentage to divide the allocation's flow amount between the given WMAs during the cumulation process to hopefully provide more realistic totals of water use in each area. However, not every state provides areas of water use, so like HarDWR v1.0, a hierarchical decision tree was used to assign each allocation to a WMA. First, if a WMA could be identified based on the allocation ID, then that WMA was used; typically, when available, this applied to the entire state and no further steps were needed. Second was the spatial analysis of Place of Use to WMAs. Third was a spatial analysis of the POD locations to WMAs, with the assumption that allocation's POD is within the WMA it should belong to; if an allocation still had multiple WMAs based on its POD locations, then the allocation's flow amount would be divided equally between all WMAs. The fourth, and final, process was to include water allocations which spatially fell outside of the state WMA boundaries. This could be due to several reasons, such as coordinate errors / imprecision in the POD location, imprecision in the WMA boundaries, or rights attached with features, such as a reservoir, which crosses state boundaries. To include these records, we decided for any POD which was within one kilometer of the state's edge would be assigned to the nearest WMA. Other Changes WestDAAT has Allowed In addition to a more nuanced and consistent method of assigning water right's data to WMAs, there are other benefits gained from using the WestDAAT dataset. Among those is a consistent categorization of a water right's legal status. In HarDWR v1.0, legal status was effectively ignored, which led to many valid concerns about the quality of the database related to the amounts of water the rights allowed to be claimed. The main issue was that rights with legal status' such as "application withdrawn", "non-active", or "cancelled" were included within HarDWR v1.0. These, and other water rights status' which were deemed to not be in use have been removed from this version of the database. Another major change has been the addition of the "unspecified water source category. This is water that can come from either surface water or groundwater, or the source of which is unknown. The addition of this source category brings the total number of categories to three. Due to reviewer feedback, we decided to add the acre-foot (AF) column so that the data may be more applicable to a wider audience. We added the ownerClassification column so that the data may be more applicable to a wider audience. File Descriptions The dataset is a series of various files organized by state sub-directories. In addition, each file begins with the state's name, in case the file is separate from its sub-directory for some reason. After the state name is the text which describes the contents of the file. Here is each file described in detail. Note that st is a placeholder for the state's name. stFullRecords_HarmonizedRights.csv: A file of the complete water records for each state. The column headers for each of this type of file are: state - The name of the state to which the allocations belong to. FIPS - The two digit numeric state ID code. siteID - The site location ID for POD locations. A site may have multiple allocations, which are the actual amount of water which can be drawn. In a simplified hypothetical, a farm stead may have an allocation for "irrigation" and an allocation for "domestic" water use, but the water is drawn from the same pumping equipment. It should be noted that many of the site ID appear to have been added by WaDE, and therefore may not be recognized by a given state's water rights database. allocationID - The allocation ID for the water right. For most states this is the water right ID, and what is recommended to use should a right be looked up on a given state's water rights database. The water amounts associated with these IDs tend to be finer scaled than those associated with siteID. It should be noted that some allocations may be extracted from multiple sites, particularly for larger Places of Use. ownerClassification - A classification of the types of owners for water rights. The most common is `Private` which incorporates a wide range of entities. Several classifications would be grouped into a government category, most of which are for the U.S. Federal Government. These allocations could be listed as "Federal", "United States of America", or as the names of any number of federal agencies. The last major grouping of entities is for "Native American"s. priorityDate - The date we use as the water right priority date for our modeling analysis. This is the legal priority date when it is available. However, for some rights, specifically from California and New Mexico, we used a pseudo priority date (e.g. well completion date or start of well drilling date) when a legal priority date was not available. The most questionable dates come from New Mexico, where the only date associated with certain water right records was the date the allocation was recorded in the database. As the allocation record creation tended to be within a few months of the filing of the application of the water right, from manually double checking the water rights, and our analysis focuses on aggregating water rights on the timescale of years, we determined it was acceptable to use such dates to include as many records as possible. primaryBeneficialUse - From the numerous state water use categories, WaDE categorized them into 21 categories WestDAAT. This column is the original WaDE category for the primary water use at the PoD site. allocationBeneficialUse - From the numerous state water use categories, WaDE categorized them into 21 categories for WestDAAT. This column is the original WaDE category

Economics↗
DOE OSTI2024

Idaho National Laboratory Integrated Multisite SSHAC Level 3: Probabilistic Volcanic Hazards Assessment

The Idaho National Laboratory (INL) resides on the eastern Snake River Plain (ESRP), part of the Snake River Plain (SRP) with a complex origin and geologic history of volcanism. Much of the Quaternary (last 2.58 million years) volcanism within 400 km of INL is genetically associated with a major thermal anomaly referred to as the Yellowstone hotspot, which is currently located more than 180 km northeast of INL. Near INL, local volcanic sources include silicic domes near its southern border, and numerous dike-fed basaltic vents of the ESRP, some of which are near or within the INL boundaries. Quantitative probabilistic assessments of screened-in volcanic hazardous phenomena for nine different facility complexes at the INL are presented for a Senior Seismic Hazard Analysis Committee (SSHAC) Level 3 (SL3) study. The comprehensive, integrated multisite SSHAC study consists of a single regional Probabilistic Volcanic Hazards Assessment (PVHA) that pertains to all nine INL facility complexes, with site-specific information developed at each respective area of interest (AOI), referred to as the "INL facility AOI" (Figure ES-1). Hazard products are generated for specified INL facility AOIs for use by multiple stakeholders from different agencies in risk-informed decision-making regarding site selection, operations, and design of nuclear facilities at INL, consistent with U.S. Department of Energy (DOE) and U.S. Nuclear Regulatory Commission (NRC) regulatory guidance. The study also serves as the basis for future periodic safety assessments required for existing DOE facilities, such as 10-year evaluations of natural-phenomena hazards. Elements of the INL PVHA including initial characterization, screening, quantitative assessments of eruption and hazard potential, and consideration of facility needs are conducted using the three phases of the SSHAC process: evaluation, integration, and documentation. As per regulatory guidance, the SSHAC framework provided the necessary processes and procedures for the PVHA Technical Integration (TI) team, at three workshops, five formal working meetings and many TI team remote meetings, to conduct initial characterization, screen the volcanic hazards, create PVHA model inputs, exercise those models in the PVHA, consider facility-specific volcanic hazard needs, and perform final hazard calculations. The Participatory Peer Review Panel (PPRP) provided independent oversight and performed process and technical reviews of the PVHA throughout its duration. The study included an extensive New Data Collection and Analyses (NDCA) program developed by the PVHA TI team to reduce uncertainties in hazard-significant elements in the PVHA model. NDCA activities generated 20 reports providing important contributory datasets and results to the project database for characterizing the ESRP. For example, a report compiling the dimensions of ESRP shield volcanoes and lava fields was used to construct volcanic footprints (areas of impact), was compared with data from INL subsurface cores, and was used to validate the results of lava-flow inundation modeling on the contemporary terrain. Another example is the acquisition of aeromagnetic data over INL and its surrounding area, with maps of buried magmatic features (e.g., subsurface volcanoes and swarms of feeder dikes) that informed the PVHA conceptual model of volcanism. The SSHAC evaluation process was used to conduct all elements of the PVHA including initial characterization and screening. Existing data and NDCA activities provided the foundation to develop the tectonomagmatic conceptual model of volcanism for the region of geographical interest in the SRP and Yellowstone hotspot volcanic system, and for considering volcanoes in the western US. Considering Quaternary volcanic sources active during this period, the screening approach identified and evaluated magma compositions, types of eruptive and intrusive phenomena, types of hazardous phenomena, and proximity of sources to INL. The PVHA TI team evaluated 18 types of magmatic sources in terms of 20 potentially hazardous phenomena, resulting in 360 screening decisions. The screening process resulted in 140 screened-in hazardous volcanic phenomena for Quaternary volcanic sources 1) proximal to INL facility complexes in the ESRP (64 basaltic and 54 silicic), 2) regional sources associated with the Yellowstone caldera system and Blackfoot Reservoir volcanic field (7), and 3) more distal sources from thirteen Cascade volcanoes and two volcanoes at Long Valley caldera (CA).

58 GEOSCIENCES↗