Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 253 records · Page 14
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↗
NASA NTRS2013

An Elevated Reservoir of Air Pollutants over the Mid-Atlantic States During the 2011 DISCOVER-AQ Campaign: Airborne Measurements and Numerical Simulations

During a classic heat wave with record high temperatures and poor air quality from July 18 to 23, 2011, an elevated reservoir of air pollutants was observed over and downwind of Baltimore, MD, with relatively clean conditions near the surface. Aircraft and ozonesonde measurements detected approximately 120 parts per billion by volume ozone at 800 meters altitude, but approximately 80 parts per billion by volume ozone near the surface. High concentrations of other pollutants were also observed around the ozone peak: approximately 300 parts per billion by volume CO at 1200 meters, approximately 2 parts per billion by volume NO2 at 800 meters, approximately 5 parts per billion by volume SO2 at 600 meters, and strong aerosol optical scattering (2 x 10 (sup 4) per meter) at 600 meters. These results suggest that the elevated reservoir is a mixture of automobile exhaust (high concentrations of O3, CO, and NO2) and power plant emissions (high SO2 and aerosols). Back trajectory calculations show a local stagnation event before the formation of this elevated reservoir. Forward trajectories suggest an influence on downwind air quality, supported by surface ozone observations on the next day over the downwind PA, NJ and NY area. Meteorological observations from aircraft and ozonesondes show a dramatic veering of wind direction from south to north within the lowest 5000 meters, implying that the development of the elevated reservoir was caused in part by the Chesapeake Bay breeze. Based on in situ observations, Community Air Quality Multi-scale Model (CMAQ) forecast simulations with 12 kilometers resolution overestimated surface ozone concentrations and failed to predict this elevated reservoir; however, CMAQ research simulations with 4 kilometers and 1.33 kilometers resolution more successfully reproduced this event. These results show that high resolution is essential for resolving coastal effects and predicting air quality for cities near major bodies of water such as Baltimore on the Chesapeake Bay and downwind areas in the Northeast.

ozone↗
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↗
NASA NTRS2019

Prebiotic Chemistry of Pluto

We present the case for the presence of complex organic molecules, such as amino acids and nucleobases, formed by abiotic processes on the surface and in near-subsurface regions of Pluto. Pluto's surface is tinted with a range of non-ice substances with colors ranging from light yellow to red to dark brown; the colors match those of laboratory organic residues called tholins. Tholins are broadly characterized as complex, macromolecular organic solids consisting of a network of aromatic structures connected by aliphatic bridging units (e.g., Imanaka et al.,2004; Materese et al.,2014, 2015). The synthesis of tholins in planetary atmospheres and in surface ices has been explored in numerous laboratory experiments, and both gas- and solid-phase varieties are found on Pluto. A third variety of tholins, exposed at a site of tectonic surface fracturing called Virgil Fossae, appears to have come from a reservoir in the subsurface. Eruptions of tholin-laden liquid H2O from a subsurface aqueous repository appear to have covered portions of Virgil Fossae and its surroundings with a uniquely colored deposit (D.P. Cruikshank, personal communication) that is geographically correlated with an exposure of H2O ice that includes spectroscopically detected NH3 (C.M. Dalle Ore, personal communication). The subsurface organic material could have been derived from presolar or solar nebula processes, or might have formed in situ. Photolysis and radiolysis of a mixture of ices relevant to Pluto's surface composition (N2, CH4, CO) have produced strongly colored, complex organics with a significant aromatic content having a high degree of nitrogen substitution similar to the aromatic heterocycles pyrimidine and purine (Materese et al.,2014, 2015; Cruikshank et al.,2016). Experiments with pyrimidines and purines frozen in H2O-NH3 ice resulted in the formation of numerous nucleobases, including the biologically relevant guanine, cytosine, adenine, uracil, and thymine (Materese et al.,2017). The red material associated with the H2O ice may contain nucleobases resulting from energetic processing on Pluto's surface or in the interior. Some other Kuiper Belt objects also exhibit red colors similar to those found on Pluto and may therefore carry similar inventories of complex organic materials. The widespread and ubiquitous nature of similarly complex organic materials observed in a variety of astronomical settings drives the need for additional laboratory and modeling efforts to explain the origin and evolution of organic molecules. Pluto observations reveal complex organics on a small body that remains close to its place of origin in the outermost regions of the Solar System.

Cruikshank, D. P.↗
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↗
NASA NTRS2017

Imaging a Free-Geometry Magma Reservoir Using FEM-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluidfilled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗
NASA NTRS2017

Imaging the Complex Geometry of a Magma Reservoir Using Fem-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluid-filled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗
NASA NTRS1997

A blast of gas in the latest Paleocene: simulating first-order effects of massive dissociation of oceanic methane hydrate

Carbonate and organic matter deposited during the latest Paleocene thermal maximum is characterized by a remarkable -2.5% excursion in delta 13C that occurred over approximately 10(4) yr and returned to near initial values in an exponential pattern over approximately 2 x 10(5) yr. It has been hypothesized that this excursion signifies transfer of 1.4 to 2.8 x 10(18) g of CH4 from oceanic hydrates to the combined ocean-atmosphere inorganic carbon reservoir. A scenario with 1.12 x 10(18) g of CH4 is numerically simulated here within the framework of the present-day global carbon cycle to test the plausibility of the hypothesis. We find that (1) the delta 13C of the deep ocean, shallow ocean, and atmosphere decreases by -2.3% over 10(4) yr and returns to initial values in an exponential pattern over approximately 2 x 10(5) yr; (2) the depth of the lysocline shoals by up to 400 m over 10(4) yr, and this rise is most pronounced in one ocean region; and (3) global surface temperature increases by approximately 2 degrees C over 10(4) yr and returns to initial values over approximately 2 x 10(6) yr. The first effect is quantitatively consistent with the geologic record; the latter two effects are qualitatively consistent with observations. Thus, significant CH4 release from oceanic hydrates is a plausible explanation for observed carbon cycle perturbations during the thermal maximum. This conclusion is of broad interest because the flux of CH4 invoked during the maximum is of similar magnitude to that released to the atmosphere from present-day anthropogenic CH4 sources.

Non-NASA Center↗
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↗