Search NASA⌕ Search

SEARCH · Search NASA

Results for “active source”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

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

High Spatial Melt Rate Variability Near the Totten Glacier Grounding Zone Explained by New Bathymetry Inversion

Totten Glacier is a fast-moving East Antarctic outlet with the potential for significant future sea-level contributions. We deployed four autonomous phase-sensitive radars on its ice shelf to monitor ice-ocean interactions near its grounding zone and made active source seismic observations to constrain gravity-derived bathymetry models. We observe an asymmetry in basal melting with mean melt rates along the grounding zone differing by up to 20 m/a. Our new bathymetry model reveals that this melt rate asymmetry coincides with an asymmetry in water column thickness and that the low-melting ice-shelf portion is shielded from the main cavity circulation. A 2-year record yields year-to-year melt rate variability of 7–9 m/a with no seasonal cycle. Our results highlight the key role of bathymetry near grounding lines for accurate modeling of ice-shelf melt, and the importance of sustained multi-year monitoring, especially at ice-shelf cavities where the dominant melt rate drivers vary primarily inter-annually.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Moisture-induced autonomous surface potential oscillations for energy harvesting

A variety of autonomous oscillations in nature such as heartbeats and some biochemical reactions have been widely studied and utilized for applications in the fields of bioscience and engineering. Here, we report a unique phenomenon of moisture-induced electrical potential oscillations on polymers, poly([2-(methacryloyloxy)ethyl] dimethyl-(3-sulfopropyl) ammonium hydroxide-co-acrylic acid), during the diffusion of water molecules. Chemical reactions are modeled by kinetic simulations while system dynamic equations and the stability matrix are analyzed to show the chaotic nature of the system which oscillates with hidden attractors to induce the autonomous surface potential oscillation. Using moisture in the ambient environment as the activation source, this self-excited chemoelectrical reaction could have broad influences and usages in surface-reaction based devices and systems. As a proof-of-concept demonstration, an energy harvester is constructed and achieved the continuous energy production for more than 15,000 seconds with an energy density of 16.8 mJ/cm 2 . A 2-Volts output voltage has been produced to power a liquid crystal display toward practical applications with five energy harvesters connected in series.

42 ENGINEERING↗
DOE OSTI2022

Direct neutrino-mass measurement with sub-electronvolt sensitivity

Since the discovery of neutrino oscillations, we know that neutrinos have non-zero mass. However, the absolute neutrino-mass scale remains unknown. Here we report the upper limits on effective electron anti-neutrino mass, m ν , from the second physics run of the Karlsruhe Tritium Neutrino experiment. In this experiment, m ν is probed via a high-precision measurement of the tritium β-decay spectrum close to its endpoint. This method is independent of any cosmological model and does not rely on assumptions whether the neutrino is a Dirac or Majorana particle. By increasing the source activity and reducing the background with respect to the first physics campaign, we reached a sensitivity on m ν –2 at a 90% confidence level (CL). The best fit to the spectral data yields m$^2_ν$ = (0.26 ± 0.34) eV 2 c –4 , resulting in an upper limit of m ν < 0.9 eV c –2 at 90% CL. By combining this result with the first neutrino-mass campaign, we find an upper limit of m ν < 0.8 eV c –2 at 90% CL.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Physics informed neural network can retrieve rate and state friction parameters from acoustic monitoring of laboratory stick-slip experiments

Various machine learning (ML) and deep learning (DL) techniques have been recently applied to the forecasting of laboratory earthquakes from friction experiments. The magnitude and timing of shear failures in stick-slip cycles are predicted using features extracted from the recorded ultrasonic or acoustic emission (AE) signals. In addition, the Rate and State Friction (RSF) constitutive laws are extensively used to model the frictional behavior of faults. In this work, we use data from shear experiments coupled with passive acoustic (variance, kurtosis, and AE rate) interleaved with active source ultrasonic monitoring (transmitted wave amplitude) to develop physics-informed neural network (PINN) models incorporating the RSF law and AE rate generation equation with wave amplitude serving as a proxy for friction state variable. This PINN framework allows learning RSF parameters from stick-slip experiments rather than measuring them through a series of velocity step experiments. We observe that when the stick-slip cycles are irregular, the PINN models outperform the data-driven DL models. Transfer learning (TL) PINN models are also developed by pre-training on data collected at one normal stress level followed by forecasting shear failures and retrieving RSF parameters at other stress levels (i.e., with different recurrence intervals) after retraining on a limited amount of new data. Our findings suggest that TL models perform better compared to standalone models. Both standalone and TL PINN-estimated RSF parameters and their ground truth values show excellent agreements thus demonstrating that RSF parameters can be retrieved from laboratory stick-slip experiments using the corresponding acoustic data and that the transmitted wave amplitude provides a good representation of the evolving frictional state during stick-slips.

58 GEOSCIENCES↗
DOE OSTI2012

INITIATION AND DEVELOPMENT OF THE WHITE-LIGHT AND RADIO CORONAL MASS EJECTION ON 2001 APRIL 15

The 2001 April 15 event was one of the largest of the last solar cycle. A former study established that this event was associated with a coronal mass ejection (CME) observed both at white light and radio frequencies. This radio CME is illuminated by synchrotron emission from relativistic electrons. In this paper, we investigate the relation of the radio CME to its extreme-ultraviolet (EUV) and white-light counterpart and reach four main conclusions. (1) The radio CME corresponds to the white-light flux rope cavity. (2) The presence of a reconnecting current sheet behind the erupting flux rope is framed, both from below and above, by bursty radio sources. This reconnection is the source of relativistic radiating electrons which are injected down along the reconnected coronal arches and up along the flux rope border forming the radio CME. (3) Radio imaging reveals an important lateral overexpansion in the low corona; this overexpansion is at the origin of compression regions where type II and III bursts are imaged. (4) Already in the initiation phase, radio images reveal large-scale interactions of the source active region (AR) with its surroundings, including another AR and open magnetic fields. Thus, these complementary radio, EUV, and white-light data validate the flux rope eruption model of CMEs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

AN EXTREME-ULTRAVIOLET WAVE ASSOCIATED WITH A SURGE

Taking advantage of the high temporal and spatial resolution observations from the Solar Dynamics Observatory, we present an extreme-ultraviolet (EUV) wave associated with a surge on 2010 November 13. Due to the magnetic flux cancelation, some surges formed in the source active region (AR). The strongest surge produced our studied event. The surge was deflected by the nearby loops that connected to another AR, and disrupted the overlying loops that slowly expanded and eventually evolved into a weak coronal mass ejection (CME). The surge was likely associated with the core of the CME. The EUV wave happened after the surge deflected. The wave departed far from the flare center and showed a close location relative to the deflected surge. The wave propagated in a narrow angular extent, mainly in the ejection direction of the surge. The close timing and location relations between the EUV wave and the surge indicate that the wave was closely associated with the CME. The wave had a velocity of 310-350 km s{sup -1}, while the speeds of the surge and the expanding loops were about 130 and 150 km s{sup -1}, respectively. All of the results suggest that the EUV wave was a fast-mode wave and was most likely triggered by the weak CME.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

PRISM3D: a 3-D reference seismic model for Iberia and adjacent areas

SUMMARY We present PRISM3D, a 3-D reference seismic model of P- and S-wave velocities for Iberia and adjacent areas. PRISM3D results from the combination of the most up-to-date earth models available for the region. It extends horizontally from 15°W to 5°E in longitude, 34°N to 46°N in latitude and vertically from 3.5 km above to 200 km below sea level, and is modelled on a regular grid with 10 and 0.5 km of grid node spacing in the horizontal and vertical directions, respectively. It was designed using models inferred from local and teleseismic body-wave tomography, earthquake and ambient noise surface wave tomography, receiver function analysis and active source experiments. It includes two interfaces, namely the topography/bathymetry and the Mohorovičić (Moho) discontinuity. The Moho was modelled from previously published receiver function analysis and deep seismic sounding results. To that end we used a probabilistic surface reconstruction algorithm that allowed to extract the mean of the Moho depth surface along with its associated standard deviation, which provides a depth uncertainty estimate. The Moho depth model is in good agreement with previously published models, although it presents slightly sharper gardients in orogenic areas such as the Pyrenees or the Betic-Rif system. Crustal and mantle P- and S-wave wave speed grids were built separately on each side of the Moho depth surface by weighted average of existing models, thus allowing to realistically render the speed gradients across that interface. The associated weighted standard deviation was also calculated, which provides an uncertainty estimation on the average wave speed values at any point of the grid. At shallow depths (<10 km), low P and S wave speeds and high VP/VS are observed in offshore basins, while the Iberian Massif, which covers a large part of western Iberia, appears characterized by a rather flat Moho, higher than average VP and VS and low VP/VS. Conversely, the Betic-Rif system seems to be associated with low VP and VS, combined with high VP/VS in comparison to the rest of the study area. The most prominent feature of the mantle is the well known high wave speed anomaly related to the Alboran slab imaged in various mantle tomography studies. The consistency of PRISM3D with previous work is verified by comparing it with two recent studies, with which it shows a good general agreement.The impact of the new 3-D model is illustrated through a simple synthetic experiment, which shows that the lateral variations of the wave speed can produce traveltime differences ranging from –1.5 and 1.5 s for P waves and from –2.5 and 2.5 s for S waves at local to regional distances. Such values are far larger than phase picking uncertainties and would likely affect earthquake hypocentral parameter estimations. The new 3-D model thus provides a basis for regional studies including earthquake source studies, Earth structure investigations and geodynamic modelling of Iberia and its surroundings.

Arroucau, P.↗
DOE OSTI2021

Phase-weighted slant stacking for surface wave dispersion measurement

SUMMARY Surface wave retrieval from ambient noise records using seismic interferometry techniques has been widely used for multiscale shear wave velocity (Vs) imaging. One key step during Vs imaging is the generation of dispersion spectra and the extraction of a reliable dispersion curve from the retrieved surface waves. However, the sparse array geometry usually affects the ability for high-frequency (>1 Hz) seismic signals’ acquisition. Dispersion measurements are degraded by array response due to sparse sampling and often present smeared dispersion spectra with sidelobe artefacts. Previous studies usually focus on interferograms’ domain (e.g. cross-correlation function) and attempt to enhance coherent signals before dispersion measurement. We propose an alternative technique to explicitly deblur dispersion spectra through use of a phase-weighted slant-stacking algorithm. Numerical examples demonstrate the strength of the proposed technique to attenuate array responses as well as incoherent noise. Three different field examples prove the flexibility and superiority of the proposed technique: the first data set consists of ambient noise records acquired using a nodal seismometer array; the second data set utilizes distributed acoustic sensing (DAS) and a marine fibre-optic cable to acquire a similar ambient noise data set; the last data set is a vibrator-based active-source surface wave data. The enhanced dispersion measurements provide cleaner and higher-resolution spectra without distortions which will assist both human interpreters as well as ML algorithms in efficiently picking curves for subsequent Vs inversion.

58 GEOSCIENCES↗
DOE OSTI2018

Ground motion response to an ML 4.3 earthquake using co-located distributed acoustic sensing and seismometer arrays

The PoroTomo research team deployed two arrays of seismic sensors in a natural laboratory at Brady Hot Springs, Nevada in March 2016. The 1500 m (length) by 500 m (width) by 400 m (depth) volume of the laboratory overlies a geothermal reservoir. The surface Distributed Acoustic Sensing (DAS) array consisted of 8700 m of fiber-optic cable in a shallow trench, including 340 m in a well. The conventional seismometer array consisted of 238 three- component geophones. The DAS cable was laid out in three parallel zig-zag lines with line segments approximately 100 meters in length and geophones were spaced at approximately 60- meter intervals. Both DAS and conventional geophones recorded continuously over 15 days during which a moderate-sized earthquake with a local magnitude of 4.3 was recorded on March 21, 2016. Its epicenter was approximately 150-km south-southeast of the laboratory. Several DAS line segments with co-located geophone stations were used to compare signal-to-noise (SNR) ratios in both time and frequency domains and to test relationships between DAS and geophone data. The ratios were typically within a factor of five of each other with DAS SNR often greater for P-wave but smaller for S-wave relative to geophone SNR. The SNRs measured for an earthquake can be better than for active sources, because the earthquake signal contains more low frequency energy and the noise level is also lower at those lower frequencies. Amplitudes of the sum of several DAS strain-rate waveforms matched the finite difference of two geophone waveforms reasonably well, as did the amplitudes of DAS strain waveforms with particle-velocity waveforms recorded by geophones. Similar agreement was found between DAS and geophone observations and synthetic strain seismograms. In conclusion, the combination of good SNR in the seismic frequency band, high-spatial density, large N, and highly accurate time control among individual sensors suggests that DAS arrays have potential to assume a role in earthquake seismology.

60 signal-to-noise ratio↗
DOE OSTI2016

Using Atmospheric Dispersion Theory to Inform the Design of a Short-lived Radioactive Particle Release Experiment

Atmospheric dispersion theory can be used to predict ground deposition of particulates downwind of a radionuclide release. This paper utilizes standard formulations found in Gaussian plume models to inform the design of an experimental release of short-lived radioactive particles into the atmosphere. Specifically, a source depletion algorithm is used to determine the optimum particle size and release height that maximizes the near-field deposition while minimizing the both the required source activity and the fraction of activity lost to long-distance transport. The purpose of the release is to provide a realistic deposition pattern that might be observed downwind of a small-scale vent from an underground nuclear explosion. The deposition field will be used, in part, to investigate several techniques of gamma radiation survey and spectrometry that could be utilized by an On-Site Inspection team under the verification regime of the Comprehensive Nuclear-Test-Ban Treaty.

Atmospheric dispersion↗
DOE OSTI2018

Precise measurement of α K and α T for the 39.8-keV E 3 transition in Rh 103 : Test of internal-conversion theory

Neutron-activated sources of 103 Ru and 103 Pd both share the isomeric first-excited state in 103 Rh as a daughter product. From independent measurements of both decays, we have measured the K shell and total internal-conversion coefficients, α K and α T , for the 39.8-keV E3 transition, which de-excites that state in 103 Rh, to be 141.1(23) and 1428(13), respectively. When compared with Dirac-Fock calculations, our new results disagree with the version of the theory that ignores the K-shell atomic vacancy, which is consistent with our conclusion drawn from a series of measurements on high multipolarity transitions in nuclei with higher Z. Calculations that include the atomic vacancy indicate that the transition actually has a small M4 component with mixing ratio δ = 0.023(5).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Size-Tuned Plastic Flow Localization in Irradiated Materials at the Submicron Scale

Three-dimensional discrete dislocation dynamics (3D-DDD) simulations reveal that, with reduction of sample size in the submicron regime, the mechanism of plastic flow localization in irradiated materials transitions from irradiation-controlled to an intrinsic dislocation source controlled. Furthermore, the spatial correlation of plastic deformation decreases due to weaker dislocation interactions and less frequent cross slip as the system size decreases, thus manifesting itself in thinner dislocation channels. A simple model of discrete dislocation source activation coupled with cross slip channel widening is developed to reproduce and physically explain this transition. In order to quantify the phenomenon of plastic flow localization, we introduce a “deformation localization index,” with implications to the design of radiation-resistant materials.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

Temperature dependence of APD-based PET scanners

Purpose: Solid state detectors such as avalanche photodiodes (APDs) are increasingly being used in PET detectors. One of the disadvantages of APDs is the strong decrease of their gain factor with increasing ambient temperature. The light yield of most scintillation crystals also decreases when ambient temperature is increased. Both effects lead to considerable temperature dependence of the performance of APD-based PET scanners. In this paper, the authors propose a model for this dependence and the performance of the LabPET8 APD-based small animal PET scanner is evaluated at different temperatures.Methods: The model proposes that the effect of increasing temperature on the energy histogram of an APD-based PET scanner is a compression of the histogram along the energy axis. The energy histogram of the LabPET system was acquired at 21 °C and 25 °C to verify the validity of this model. Using the proposed model, the effect of temperature on system sensitivity was simulated for different detector temperature coefficients and temperatures. Subsequently, the effect of short term and long term temperature changes on the peak sensitivity of the LabPET system was measured. The axial sensitivity profile was measured at 21 °C and 24 °C following the NEMA NU 4-2008 standard. System spatial resolution was also evaluated. Furthermore, scatter fraction, count losses and random coincidences were evaluated at different temperatures. Image quality was also investigated.Results: As predicted by the model, the photopeak energy at 25 °C is lower than at 21 °C with a shift of approximately 6% per °C. Simulations showed that this results in an approximately linear decrease of sensitivity when temperature is increased from 21 °C to 24 °C and energy thresholds are constant. Experimental evaluation of the peak sensitivity at different temperatures showed a strong linear correlation for short term (2.32 kcps/MBq/°C = 12%/°C, R = −0.95) and long term (1.92 kcps/MBq/°C = 10%/°C , R = −0.96) temperature changes. Count rate evaluation showed that although the total count rate is consistently higher at 21 °C than at 24 °C for different source activity concentrations, this is mainly due to an increase in scattered and random coincidences. The peak total count rate is 400 kcps at both temperatures but is reached at lower activity at 21 °C. The peak true count rate is 138 kcps (at 100 MBq) at 21 °C and 180 kcps (at 125 MBq) at 24 °C. The peak noise equivalent count rate is also lower at 21 °C (70 kcps at 70 MBq) than at 24 °C (100 kcps at 100 MBq). At realistic activity levels, the scatter fraction is lower at higher temperatures, but at the cost of a strong decrease in true count rate.Conclusions: A model was proposed for the temperature dependence of APD-based PET scanners and evaluated using the LabPET small animal PET scanner. System sensitivity and count rate performance are strongly dependent on ambient temperature while system resolution is not. The authors’ results indicate that it is important to assure stable ambient temperature to obtain reproducible results in imaging studies with APD-based PET scanners.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

The evaluation of a 2D diode array in “magic phantom” for use in high dose rate brachytherapy pretreatment quality assurance

Purpose: High dose rate (HDR) brachytherapy is a treatment method that is used increasingly worldwide. The development of a sound quality assurance program for the verification of treatment deliveries can be challenging due to the high source activity utilized and the need for precise measurements of dwell positions and times. This paper describes the application of a novel phantom, based on a 2D 11 × 11 diode array detection system, named “magic phantom” (MPh), to accurately measure plan dwell positions and times, compare them directly to the treatment plan, determine errors in treatment delivery, and calculate absorbed dose. Methods: The magic phantom system was CT scanned and a 20 catheter plan was generated to simulate a nonspecific treatment scenario. This plan was delivered to the MPh and, using a custom developed software suite, the dwell positions and times were measured and compared to the plan. The original plan was also modified, with changes not disclosed to the primary authors, and measured again using the device and software to determine the modifications. A new metric, the “position–time gamma index,” was developed to quantify the quality of a treatment delivery when compared to the treatment plan. The MPh was evaluated to determine the minimum measurable dwell time and step size. The incorporation of the TG-43U1 formalism directly into the software allows for dose calculations to be made based on the measured plan. The estimated dose distributions calculated by the software were compared to the treatment plan and to calibrated EBT3 film, using the 2D gamma analysis method. Results: For the original plan, the magic phantom system was capable of measuring all dwell points and dwell times and the majority were found to be within 0.93 mm and 0.25 s, respectively, from the plan. By measuring the altered plan and comparing it to the unmodified treatment plan, the use of the position–time gamma index showed that all modifications made could be readily detected. The MPh was able to measure dwell times down to 0.067 ± 0.001 s and planned dwell positions separated by 1 mm. The dose calculation carried out by the MPh software was found to be in agreement with values calculated by the treatment planning system within 0.75%. Using the 2D gamma index, the dose map of the MPh plane and measured EBT3 were found to have a pass rate of over 95% when compared to the original plan. Conclusions: The application of this magic phantom quality assurance system to HDR brachytherapy has demonstrated promising ability to perform the verification of treatment plans, based upon the measured dwell positions and times. The introduction of the quantitative position–time gamma index allows for direct comparison of measured parameters against the plan and could be used prior to patient treatment to ensure accurate delivery.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-635: Process Mapping of Eye Plaque Brachytherapy

Purpose: To apply a risk-based assessment and analysis technique (AAPM TG 100) to eye plaque brachytherapy treatment of ocular melanoma. Methods: The role and responsibility of personnel involved in the eye plaque brachytherapy is defined for retinal specialist, radiation oncologist, nurse and medical physicist. The entire procedure was examined carefully. First, major processes were identified and then details for each major process were followed. Results: Seventy-one total potential modes were identified. Eight major processes (corresponding detailed number of modes) are patient consultation (2 modes), pretreatment tumor localization (11), treatment planning (13), seed ordering and calibration (10), eye plaque assembly (10), implantation (11), removal (11), and deconstruction (3), respectively. Half of the total modes (36 modes) are related to physicist while physicist is not involved in processes such as during the actual procedure of suturing and removing the plaque. Conclusion: Not only can failure modes arise from physicist-related procedures such as treatment planning and source activity calibration, but it can also exist in more clinical procedures by other medical staff. The improvement of the accurate communication for non-physicist-related clinical procedures could potentially be an approach to prevent human errors. More rigorous physics double check would reduce the error for physicist-related procedures. Eventually, based on this detailed process map, failure mode and effect analysis (FMEA) will identify top tiers of modes by ranking all possible modes with risk priority number (RPN). For those high risk modes, fault tree analysis (FTA) will provide possible preventive action plans.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

A Monte Carlo study on {sup 223}Ra imaging for unsealed radionuclide therapy

Purpose: Radium-223 ({sup 223}Ra), an α-emitting radionuclide, is used in unsealed radionuclide therapy for metastatic bone tumors. The demand for qualitative {sup 223}Ra imaging is growing to optimize dosimetry. The authors simulated {sup 223}Ra imaging using an in-house Monte Carlo simulation code and investigated the feasibility and utility of {sup 223}Ra imaging. Methods: The Monte Carlo code comprises two modules, HEXAGON and NAI. The HEXAGON code simulates the photon and electron interactions in the tissues and collimator, and the NAI code simulates the response of the NaI detector system. A 3D numeric phantom created using computed tomography images of a chest phantom was installed in the HEXAGON code. {sup 223}Ra accumulated in a part of the spine, and three x-rays and 19 γ rays between 80 and 450 keV were selected as the emitted photons. To evaluate the quality of the {sup 223}Ra imaging, the authors also simulated technetium-99m ({sup 99m}Tc) imaging under the same conditions and compared the results. Results: The sensitivities of the three photopeaks were 147 counts per unit of source activity (cps MBq{sup −1}; photopeak: 84 keV, full width of energy window: 20%), 166 cps MBq{sup −1} (154 keV, 15%), and 158 cps MBq{sup −1} (270 keV, 10%) for a low-energy general-purpose (LEGP) collimator, and those for the medium-energy general-purpose (MEGP) collimator were 33, 13, and 8.0 cps MBq{sup −1}, respectively. In the case of {sup 99m}Tc, the sensitivity was 55 cps MBq{sup −1} (141 keV, 20%) for LEGP and 52 cps MBq{sup −1} for MEGP. The fractions of unscattered photons of the total photons reflecting the image quality were 0.09 (84 keV), 0.03 (154 keV), and 0.02 (270 keV) for the LEGP collimator and 0.41, 0.25, and 0.50 for the MEGP collimator, respectively. Conversely, this fraction was approximately 0.65 for the simulated {sup 99m}Tc imaging. The sensitivity with the LEGP collimator appeared very high. However, almost all of the counts were because of photons that penetrated or were scattered in the collimator; therefore, the proportions of unscattered photons were small. Conclusions: Their simulation study revealed that the most promising scheme for {sup 223}Ra imaging is an 84-keV window using an MEGP collimator. The sensitivity of the photopeaks above 100 keV is too low for {sup 223}Ra imaging. A comparison of the fractions of unscattered photons reveals that the sensitivity and image quality are approximately two-thirds of those for {sup 99m}Tc imaging.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Solving multiphysics-based inverse problems with learned surrogates and constraints

Abstract Solving multiphysics-based inverse problems for geological carbon storage monitoring can be challenging when multimodal time-lapse data are expensive to collect and costly to simulate numerically. We overcome these challenges by combining computationally cheap learned surrogates with learned constraints. Not only does this combination lead to vastly improved inversions for the important fluid-flow property, permeability, it also provides a natural platform for inverting multimodal data including well measurements and active-source time-lapse seismic data. By adding a learned constraint, we arrive at a computationally feasible inversion approach that remains accurate. This is accomplished by including a trained deep neural network, known as a normalizing flow, which forces the model iterates to remain in-distribution, thereby safeguarding the accuracy of trained Fourier neural operators that act as surrogates for the computationally expensive multiphase flow simulations involving partial differential equation solves. By means of carefully selected experiments, centered around the problem of geological carbon storage, we demonstrate the efficacy of the proposed constrained optimization method on two different data modalities, namely time-lapse well and time-lapse seismic data. While permeability inversions from both these two modalities have their pluses and minuses, their joint inversion benefits from either, yielding valuable superior permeability inversions and CO 2 plume predictions near, and far away, from the monitoring wells.

Yin, Ziyi (ORCID:0000000250248771)↗
DOE OSTI2021

The feasibility of MT tipper data to monitor CO2 storage sites

Monitoring carbon storage sites using geophysical techniques is a critical component to the success and safety of storage programs. Currently, the primary methods of monitoring such sites are seismic and, to a much lesser extent, electromagnetics using active sources. The cost of such methods, especially seismic, can be prohibitively expensive. Natural source lectromagnetics, or magnetotellurics (MT), represents a low-cost, and underutilized method with the potential to aid CO2 monitoring efforts. Specifically, the tipper of MT data gives insight to the dimensionality of the subsurface and is able to detect the expansion front of a CO2 plume in a saline reservoir. We analyze the feasibility of using the tipper to monitor two shallow CO2 plumes and conclude that the tipper may be a suitable method for long-term monitoring.<br>

Kohnke, Colton↗