Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shell effects”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Ab initio single-neutron spectroscopic overlaps in lithium isotopes

We calculate single-neutron spectroscopic overlaps for lithium isotopes in the framework of the ab initio symmetry-adapted no-core shell model. We report the associated neutron-nucleus asymptotic normalization coefficients (ANCs) and spectroscopic factors (SFs) that are important ingredients in many reaction cross-section calculations. While spectroscopic factors have been traditionally extracted from experimental cross sections, their sensitivity on the type of reactions, the energy, and the underlying models point to the need for determining SFs from first-principle structure considerations. As illustrative examples, we present 6 Li + n, 7 Li ⁢+ n, and 8 Li ⁢+ n, and we show that the results are in a good agreement with those of other ab initio methods, where available, including the quantum Monte Carlo approach. Here, we compare ANCs and SFs to available experimentally deduced values, with a view toward expanding this study to heavier nuclei and to extracting intercluster effective interactions for input into analyses of existing and future experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Degradation of Platinum-Cobalt Alloy PEMFC Cathode Catalysts in Catalyst-Ionomer Inks

In this paper we report on studies of the effects of catalyst-ionomer ink composition: ionomer equivalent weight (EW), ink solvent, and ink mixing on a commercial PtCo alloy PEMFC cathode catalyst. X-ray absorption spectroscopy (XAS) and wide-angle X-ray scattering (WAXS) are utilized to determine catalyst atomic structure, catalyst crystallite composition, and extent of loss of Co into the ionomer-solvent phase. Three different n-propanol to water weight ratios (7:3, 5:5, and 3:7) and two different ionomers (3 M 800 EW and Nafion 1000 EW) were studied. Cobalt loss from the PtCo particles was found to increase with increasing water content in the inks and to be more extensive for the 800 EW ionomer inks, causing loss of contraction of the PtCo lattice and formation of a Pt shell-PtCo alloy core intraparticle structure.

25 ENERGY STORAGE↗
DOE OSTI2023

Electron collisional excitation cross-section measurements and modeling for select Ni-like to Ge-like gold transitions

We have experimentally determined the electron collisional excitation cross-sections for several 3d→4f and 3d→5f excitations in Ni- to Ge-like Au at energies of ~ 0.4, 1, 2, and 3 keV above threshold energy, E T , for the 3d→4f excitations ( E T ~ 2.5 keV) and ~ 0.2, 1, and 2 keV above threshold energy for the 3d→5f excitations ( E T ~ 3.3 keV). The cross-section measurements are possible by using the GSFC micro-calorimeter to record emission spectra from beam plasmas created in the Livermore EBIT-I electron beam ion trap. The cross-sections are experimentally determined from the ratio of the measured intensities of the collisionally excited lines to the intensities of the radiative recombination lines in monoenergetic electron distribution EBIT-I plasmas. The effects of polarization and Auger processes in the beam plasmas are accounted for in the cross-section determination. Experimentally determined cross-sections are compared with those from HULLAC, DWS, and FAC calculations. Finally, the measurements exhibit significant differences with the calculations of these excitation cross-sections.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Superconductivity in a strange metal

In a gas of charged particles with the density n, mass m, and charge e, the electrical conductivity σ is given by the Drude formula: σ = ne 2 /mΓ, where Γ is the scattering rate. Standard metals are well described by Landau's Fermi Liquid (FL) theory, in which the electric current is carried by quasi-particles, low-energy excitations of the FL that resemble electrons with some effective mass m*. The scattering rate Γ= 1/τ= v F /l, where t is the (momentum) relaxation) time, l is the mean-free path, and v F is the Fermi velocity, can be expressed ("Matthiessen's Rule") as a sum of contributions from various scattering channels: Γ = Γ 0 + Γ el-el + Γ el-ph + …, where Γ 0 describes scattering on lattice imperfections, Γ el-el the electron-electron scattering, Γ el-ph the electron-phonon scattering, etc. Of these, Γ 0 = v F /l 0 , where l 0 is the average distance between the defects, is temperature-independent. Γ el-el should scale as T 2 because of Fermi statistics; for two electrons to scatter on one another, both must come from the "Debye shell" of the width k B T/E F , where k B is the Boltzmann constant and E F is the Fermi energy. Γ el-ph typically grows as T 5 , so we expect this to overwhelm the other channels at a high enough T. However, since l cannot be shorter than the distance between the atoms, the total Γ saturates at Mott-Ioffe-Regel (MIR) limit, roughly v F /a 0 , where v F is the Fermi velocity and a 0 is the lattice constant. The resistivity should also saturate at low T, at ρ 0 = m*v F /ne 2 l 0 T. FL theory also describes other electronic properties; e.g., it predicts that in the magnetic field B, the resistivity of the metal should increase with B 2 , because σ(B) = σ(B=0)/(1 + (ω c /Γ) 2 ), where ω c = eB/m*.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Aerosol Optical Tweezers Elucidate the Chemistry, Acidity, Phase Separations, and Morphology of Atmospheric Microdroplets

Aerosol particles represent unique chemical environments because of their high surface area-to-volume ratio that promotes the effects of interfacial chemistry in confined environments. Properties such as viscosity, diffusivity, water content, pH, and morphology—following liquid–liquid phase separation—can strongly alter how a particle interacts with condensable vapors and reactive trace gases, thus modifying its continual evolution and environmental effects. Our understanding of this chemical evolution of atmospheric particulate matter and its environmental impacts is largely limited by our ability to directly observe how these critical particle properties respond to the addition or reactive uptake of new chemical components. Aerosol optical tweezers (AOT) stably trap particles in focused laser beams, providing positional control and the retrieval of many of these critical properties required to understand and predict the chemistry of aerosolized microdroplets. The analytical power of the AOT stems from the retrieval of the cavity-enhanced Raman spectrum induced by the trapping laser. Analysis of the whispering gallery modes (WGMs) that resonate as a standing wave around the droplet’s interface, provide high accuracy measurements of the droplet’s size, refractive index (and thus a measurement of composition), and can distinguish between core–shell, partially engulfed, and homogeneous morphologies. In this study, we have advanced the ability to determine the properties of the core and shell phases in biphasic droplets, including obtaining high-accuracy pH measurements. These capabilities were applied to perform AOT physical chemistry experiments on authentic secondary organic aerosol (SOA) produced directly in the AOT chamber by ozonolysis of terpene vapors. The propensity of the SOA to phase separate as a shell from a wide range of nonpolar to polar core phases was observed, along with the discovery of a stable emulsified state of SOA particles in an aqueous salt droplet. Micron-thick SOA shells did not impede the gain or loss of water or squalane from the core to the surrounding air, indicating no significant diffusional limitations to condensational growth or partitioning even under dry conditions. These experiments formed the foundation of a new framework that predicts how the phase-separated morphology of complex aerosols containing organic carbon evolves during continual atmospheric oxidation processes. Increases in oxidation state will quickly drive conversion from a partially engulfed to core–shell morphology that has dramatically different chemical reactivity since the core phase is completely concealed by the shell. The recent advances in the experimental capabilities of the AOT technique such as presented here enable novel experimental methodologies that provide insights into the chemistry and multidimensional properties of aerosol microdroplets, and how these coevolve and respond to continual chemical reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Acceleration of Oxygen Ions in Dipolarization Events: 1. CPS Distributions

Abstract Using an MHD simulation of near tail reconnection associated with a flow burst and the collapse (dipolarization) of the inner tail in combination with test particle tracing we study the acceleration and flux increases of energetic oxygen ions (O + ). The characteristic orbits, distributions, and acceleration mechanisms are governed by the dimensionless parameter σ = ω ci t n , where ω ci is the ion gyro frequency and t n a characteristic Alfvén time of the MHD simulation. For σ < 1, central plasma sheet (CPS) populations after the passage of the dipolarization front are found to resemble half‐shells in velocity space oriented toward dusk. They originate from within the CPS and are energized typically by a single encounter of the region of enhanced cross‐tail electric field associated with the flow burst. For larger σ values ( σ > 1) the O + distributions resemble more closely those of protons, consisting of two counter‐streaming field‐aligned beams and an, albeit more tenuous and irregular, ring population perpendicular to the magnetic field. The existence of the beams, however, depends on suitable earthward moving source populations in the plasma sheet boundary layer or the adjacent lobes. The acceleration to higher energies is found to indicate a charge dependence, consistent with a dominance of more highly charged ions at energies of a few hundred keV. As in earlier simulations, the simulated fluxes show large anisotropies and nongyrotropic effects, phase bunching, and spatially and temporally localized beams.

Birn, J.↗
DOE OSTI2020

Nonideal mixing effects in warm dense matter studied with first-principles computer simulations

Here, we study nonideal mixing effects in the regime of warm dense matter (WDM) by computing the shock Hugoniot curves of BN, MgO, and MgSiO 3 . First, we derive these curves from the equations of state (EOS) of the fully interacting systems, which were obtained using a combination of path integral Monte Carlo calculations at high temperature and density functional molecular dynamics simulations at lower temperatures. We then use the ideal mixing approximation at constant pressure and temperature to rederive these Hugoniot curves from the EOS tables of the individual elements. We find that the linear mixing approximation works remarkably well at temperatures above ~2 × 10 5 K, where the shock compression ratio exceeds ~3.2. The shape of the Hugoniot curve of each compound is well reproduced. Regions of increased shock compression, which emerge because of the ionization of L and K shell electrons, are well represented, and the maximum compression ratio of the Hugoniot curves is reproduced with high precision. Some deviations are seen near the onset of the L shell ionization regime, where ionization equilibrium in the fully interacting system cannot be well reproduced by the ideal mixing approximation. This approximation also breaks down at lower temperatures, where chemical bonds play an increasingly important role. However, the results imply that the equilibrium properties of binary and ternary mixtures in the regime of WDM can be derived from the EOS tables of the individual elements. This significantly simplifies the characterization of binary and ternary mixtures in the WDM and plasma phases, which otherwise requires large numbers of more computationally expensive first-principles computer simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Multi-channel, multi-template event reconstruction for SuperCDMS data using machine learning

SuperCDMS SNOLAB uses kilogram-scale germanium and silicon detectors to search for dark matter. Each detector has Transition Edge Sensors (TESs) patterned on the top and bottom faces of a large crystal substrate, with the TESs electrically grouped into six phonon readout channels per face. Noise correlations are expected among a detector's readout channels, in part because the channels and their readout electronics are located in close proximity to one another. Moreover, owing to the large size of the detectors, energy deposits can produce vastly different phonon propagation patterns depending on their location in the substrate, resulting in a strong position dependence in the readout-channel pulse shapes. Both of these effects can degrade the energy resolution and consequently diminish the dark matter search sensitivity of the experiment if not accounted for properly. We present a new algorithm for pulse reconstruction, mathematically formulated to take into account correlated noise and pulse shape variations. This new algorithm fits N readout channels with a superposition of M pulse templates simultaneously - hence termed the N$\times$M filter. We describe a method to derive the pulse templates using principal component analysis (PCA) and to extract energy and position information using a gradient boosted decision tree (GBDT). We show that these new N$\times$M and GBDT analysis tools can reduce the impact from correlated noise sources while improving the reconstructed energy resolution for simulated mono-energetic events by more than a factor of three and for the 71Ge K-shell electron-capture peak recoils measured in a previous version of SuperCDMS called CDMSlite to $<$ 50 eV from the previously published value of $\sim$100 eV. These results lay the groundwork for position reconstruction in SuperCDMS with the N$\times$M outputs.

Albakry, M. F. [British Columbia U.; TRIUMF]↗
DOE OSTI2023

Atomic Layers of B2 CuPd on Cu Nanocubes as Catalysts for Selective Hydrogenation

The search for highly active and selective catalysts with high precious metal atom utilization efficiency has attracted increasing interest in both the fundamental synthesis of materials and important industrial reactions. Here, in this work, we report the synthesis of Pd–Cu nanocubes with a Cu core and an ordered B2 intermetallic CuPd shell with controllable atomic layers on the surface (denoted as Cu/B2 CuPd), which can efficiently and robustly catalyze the selective hydrogenation of acetylene (C 2 H 2 ) to ethylene (C 2 H 4 ) under mild conditions. The optimized Cu/B2 CuPd with a Pd loading of 9.5 at. % exhibited outstanding performance in the C 2 H 2 semi-hydrogenation with 100% C 2 H 2 conversion and 95.2% C 2 H 4 selectivity at 90 °C. We attributed this outstanding performance to the core/shell structure with a high surface density of active Pd sites isolated by Cu in the B2 intermetallic matrix, representing a structural motif of single-atom alloys (SAAs) on the surface. The combined experimental and computational studies further revealed that the electronic states of Pd and Cu are modulated by SAAs from the synergistic effect between Pd and Cu, leading to enhanced performance compared with pristine Pd and Cu catalysts. This study provides a new synthetic methodology for making single-atom catalysts with high precious metal atom utilization efficiency, enabling simultaneous tuning of both geometric and electronic structures of Pd active sites for enhanced catalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electrical Control of the Nuclear Spin States of Rare-Earth Adatoms

Rare-earth adatoms on surfaces have been studied for potential atomic-scale magnetic storage, quantum sensing, and quantum computing applications. Despite accumulating experimental efforts, a comprehensive description of the electronic configurations of the adatoms remains elusive. Here, we investigate two charge states and several electronic configurations, including 5d and 6s valence shells, for a Sm adatom on a MgO substrate using multiconfigurational ab initio methods, for the possibility of using the Sm nuclear spin levels as qubits. For the configurations in a neutral charge state, we find that the electronic ground state is a singlet, and thus the hyperfine interaction associated with the 147 Sm nucleus is absent, which may greatly enhance nuclear spin coherence time. The degeneracy of the nuclear levels is lifted by the nuclear quadrupole interaction. We show that the splitting of the nuclear levels can be controlled by a static electric field, and that Rabi oscillations between the nuclear levels can be induced by a time-dependent electric field. For the configurations in a singly charged state, electronic Kramers doublets are formed. The electronic configurations including an unpaired 6s orbital exhibit a strong hyperfine Stark effect due to a large Fermi contact contribution to the hyperfine interaction. In these configurations, electric-field-induced Rabi oscillations between the electronic-nuclear levels can occur at frequencies up to 3 orders of magnitude higher than those for the neutral charge state. The proposed system may be experimentally observed within scanning tunneling microscopy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Model output from simulations of manganese-carbon interactions in temperate forest soil profiles

This archive contains model output, code, and scripts for simulations of coupled manganese-carbon cycling in temperate forest soil profiles. These model results were generated as part of a study investigating how manganese availability influences soil organic carbon stocks and demonstrating a new model framework for coupling carbon and manganese cycling. The simulations were in support of a manuscript: "Modeling interactive effects of manganese bioavailability, nitrogen deposition, and warming on soil carbon storage." The study addresses the research questions: How does Mn bioavailability, as driven by subsurface mineral properties, pH, and redox status, affect temperate forest soil organic carbon and litter carbon stocks?How is the relationship between Mn bioavailability and carbon cycling affected by changes in temperature and nitrogen deposition?"Model simulations were conducted in a reactive transport modeling framework using PFLOTRAN coupled to python. Multiple model simulations testing different Mn-bearing mineral solubilities, hydrological patterns, nitrogen deposition rates, and temperatures are included. Soil properties, including total and exchangeable Mn concentrations, are based on values reported for the Susquehanna Shale Hills Critical Zone Observatory (SSHCZO), a temperate forested watershed in central Pennsylvania, U.S.A where Mn cycling through vegetation has been documented.File formats include netCDF (.nc), python script (.py), shell script (.sh), plain text PFLOTRAN input file (.in), and plain text PFLOTRAN database file (.dat), and gzipped tar archive (tar.gz).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Carrier Dynamics of Polar, Semipolar, and Nonpolar InGaN/GaN LEDs Measured by Small-Signal Electroluminescence

The carrier dynamics in InGaN/GaN light-emitting diodes (LEDs) are directly tied to their efficiency and maximum modulation speed, which are important metrics for solid-state lighting, displays, and optical communication. In this work, we measure the carrier dynamics of a variety of InGaN/GaN LEDs using small-signal electroluminescence methods [1]. A rate equation approach and associated small-signal circuit are used to model carrier injection, recombination in the active region, recombination in the cladding regions, and carrier escape. The model is fit to the measured optical frequency response (S21) and input impedance (S11) of the LEDs to extract the various carrier lifetimes, the carrier density, and the radiative and non-radiative recombination rates. We specifically study planar nonpolar and semipolar LEDs, which show record-high modulation speeds for III-nitride LEDs and present the modulation characteristics of core-shell nanowire-based LEDs. The planar nonpolar m-plane ($101\bar{0}$) micro-LEDs achieve a record-high -3dB modulation bandwidth for a III-nitride LED of 1.5 GHz [2]. The -3dB response of an electrically injected nanowire-based micro-LED with nonpolar facets is also reported, showing a -3dB bandwidth of 1.2 GHz [3]. The high speed is attributed to the shorter carrier lifetime associated with the nonpolar orientation. We also study the carrier dynamics in semipolar ($20\bar{2}\bar{1}$) LEDs for various temperatures [4]. Finally, we present carrier dynamics measurements on commercial-grade c-plane epitaxy for various active region designs, including a wavelength series and a growth quality series. The wavelength series offers insight into the contributions of the quantum confined Stark effect (QCSE) and InGaN material quality on the green gap [5]. The growth quality series investigates the role of non-radiative centers on the LED performance. Extraction of the carrier dynamics using small-signal electroluminescence offers insight into the factors limiting the efficiency and high-speed performance of III-nitride emitters and can be leveraged to ultimately improve the devices.

LEDs, InGaN, electroluminescence, bandwidth, modul↗
DOE OSTI2022

Critical role of magnetic moments in heavy-fermion materials: Revisiting mysteries of SmB 6

Heavy-fermion family exhibits fascinating and often puzzling properties due to the presence of open-shell f ions and the complexity of the associated charge, orbital, and spin degrees of freedom. SmB 6 is a prototypical heavy-fermion compound that is electrically insulating but yet it displays quantum oscillations, which are a telltale signature of the metallic state. Adding to the enigma is the possibility that SmB 6 is a topological Kondo insulator. Here, by treating the spin degree of freedom on an equal footing with other degrees of freedom using the parameter-free strongly-constrained and appropriately-normed (SCAN) density functional, we explore the ground-state electronic structure of SmB6. A number of competing magnetic phases lying very closely in energy are found, indicating the key role of spin fluctuations in the material. The computed band structure, crystal-field splittings in the f-electron complex, the heavy effective electron mass at the Fermi energy, and the large specific heat are all in good agreement with the corresponding experimental results. In particular, our predicted FS explains the experimentally observed bulk quantum oscillations as well as the low electrical conductivity of SmB 6 . The topological Kondo state of SmB6 is shown to be robust regardless of its magnetic configuration. The excellent performance of SCAN in heavy-fermion systems is explained in terms of its ability to treat self-interaction errors and symmetry breaking within the framework of the density functional theory. Our study provides a new approach for modeling heavy-fermion materials

36 MATERIALS SCIENCE↗
DOE OSTI2021

Comparative performance of bio-based coatings formulated with cellulose, chitin, and chitosan nanomaterials suitable for fruit preservation

Sustainable nanomaterials (SNMs) from wood, sugarcane and crab shell were prepared and used to coat selected fruits. The properties of SNMs and selected fruits were characterized and strawberry was used as an example to test antifungal activity and freshness preservation of the SNMs. The SNMs with their nano-structured morphology form strong shear-thinning dispersions for easy spraying on fruit surfaces. The fruit surface free energy was influenced by its surface morphology, predominant surface wax components, and cutin monomers. The antifungal activity of SNMs was influenced by their surface functional groups and particle size (crystals vs fibers). The coblend of wood nanocrystals (WCNCs) and chitosan nanofiber (CSNFs) exhibited the best antifungal property, which was comparable with the performance of the fungicide thiabendazole (80 mg L -1 ). The weight loss and color change of the WCNC/CSNF coated strawberries decreased by nearly half compared with the control samples, showing coating effectiveness on preserving fruit freshness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Quantifying Uncertainties In Neutron-Alpha Scattering With Chiral Nucleon-Nucleon And Three-Nucleon Forces

We report that modern ab initio theory combined with high-quality nucleon-nucleon (NN) and three-nucleon (3N) interactions from chiral effective field theory (EFT) can provide a predictive description of low-energy light-nuclei reactions relevant for astrophysics and fusion-energy applications. However, the high cost of computations has so far impeded a complete analysis of the uncertainty budget of such calculations. Starting from NN potentials up to fifth order (N 4 LO) combined with leading-order 3N forces, we study how the order-by-order convergence of the chiral expansion and confidence intervals for the 3N contact and contact-plus-one-pion-exchange low-energy constants (c E and c D ) contribute to the overall uncertainty budget of many-body calculations of neutron- 4 He (n-α) elastic scattering. We compute structure and reaction observables for three-, four-, and five-nucleon systems within the ab initio frameworks of the no-core shell model and no-core shell model with continuum. Using a small set of design runs, we construct a Gaussian process model (GPM) that acts as a statistical emulator for the theory. With this, we gain insight into how uncertainties in the 3N low-energy constants propagate throughout the calculation and determine the Bayesian posterior distribution of these parameters with Markov-Chain Monte Carlo. We find rapidly converging n-α phase shifts with respect to the chiral order. With the adopted leading-order 3N force, calculations based on the NN interaction at N 4 LO of Entem, Machleidt, and Nosyk are unable to reproduce the experimental phase shifts in the 3/2 - channel within the estimated chiral truncation errors. Closer agreement with empirical data is found when using an older parametrization of the NN interaction at order N 3 LO, and the position and width of the P-wave resonances can be used to reduce the uncertainty of the 3N low-energy constants. The present results point to a lack of spin-orbit strength when the newer parametrization of the chiral NN force up to fifth order is combined with the leading-order 3N force. The inclusion of higher-order 3N-force terms may be required to recover the missing strength. GPMs can act as fast and accurate emulators of ab initio many-body calculations of low-energy scattering and reactions of light nuclei, opening the way to a robust quantification of theoretical uncertainties grounded in the description of the underlying chiral Hamiltonian.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Charge-lattice coupling and the dynamic structure of the U–O distribution in UO 2+x

The different structures and behaviors of UO 2+x observed in crystallographic and local structure measurements were examined by extended X-ray absorption fine structure (EXAFS) measurements of pristine UO 2.0 , p + and He 2+ irradiated UO 2.0 , and, at multiple temperatures, bulk U 4 O 9 and U 3 O 7 and thin film U 4 O 9-δ on an epitaxial substrate. The disorder caused by irradiation is mostly limited to increased widths of the existing U–O/U pair distributions, with any new neighbor shells being minor. As has been previously reported, the disorder caused by oxidative addition to U 4 O 9 and U 3 O 7 is much more extensive, resulting in multisite U–O distributions and greater reduction of the U–U amplitude with different distributions in bulk and thin-film U 4 O 9 . This includes the significant spectral feature near R = 1.2 Å for all U 4 O 9 and U 3 O 7 samples fit with a U-oxo type moiety with a U–O distance around 1.7 Å. In addition to indicating that these anomalies only occur in mixed valence materials, this work confirms the continuous rearrangement of the U–O distributions from 10 to 250 K. Although these variations of the structure are not observed in crystallography, their prominence in the EXAFS indicates that the dynamic structure underlying these effects is an essential factor of these materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
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↗