Search NASA⌕ Search

SEARCH · Search NASA

Results for “local electric field mapping”

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

Uncovering microstructure and composition susceptibility of high permittivity ceramic granules to plasma-induced arc breakdown

Effective diversion of surge currents is vital to prevent unwanted damage to sensitive electronics. Among the most successful and efficient strategies relies on a dielectric stimulated arc breakdown mechanism with high permittivity ceramic granules in a spark-gap geometry. Although generally regarded as a self-healing process, substantial energy deposition may occur that, over time, diminishes the ability to withstand repeated electrical assaults. We investigate the susceptibility of lead–magnesium–niobate–lead titanate (PMN–PT) granule microstructure and composition changes following many exposures to high voltage impulses resulting in arc breakdown. Scanning electron microscopy and energy-dispersive spectroscopy mapping reveal a broad range of thermal and mechanical defects entailing thermal reduction of constituent PMN–PT metal ions and recasting due to rapid eruption of volatile species. Additionally, evidence of local melting and microcracking are apparent that can have deleterious impact on the proper function of the granules, namely, the ability to concentrate electric fields across air gaps to establish and sustain discharge pathways. We propose that the localized nature of damage and stochasticity associated with the dielectric stimulated breakdown mechanism may allow granules to maintain functionality provided no permanent conduction paths are established.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Development of a Multi-Robot System for Autonomous Inspection of Nuclear Waste Tank Pits

This paper introduces the overall design plan, development timeline, and preliminary progress of the Autonomous Pit Exploration System project. This project aims to develop an advanced multi-robot system for the efficient inspection of nuclear waste-storage tank pits. The project is structured into three phases: Phase 1 involves data collection and interface definition in collaboration with Hanford Site experts and university partners, focusing on tank riser geometry and hardware solutions. Phase 2 includes the selection of sensors and robot components, detailed mechanical design, and prototyping. Phase 3 integrates all components into a cohesive system managed by a master control package which also incorporates digital twin and surrogate models, and culminates in comprehensive testing and validation at a simulated tank pit at the Idaho National Laboratory. Additionally, the system’s communication design ensures coordinated operation through shared data, power, and control signals. For transportation and deployment, an electric vehicle (EV) is chosen to support the system for a full 10 h shift with better regulatory compliance for field deployment. A telescopic arm design is selected for its simple configuration and superior reach capability and controllability. Preliminary testing utilizes an educational robot to demonstrate the feasibility of splitting computational tasks between edge and cloud computers. Successful simultaneous localization and mapping (SLAM) tasks validate our distributed computing approach. More design considerations are also discussed, including radiation hardness assurance, SLAM performance, software transferability, and digital twinning strategies.

Nuclear waste management↗
DOE OSTI2024

Mapping twist-tuned multiband topology in bilayer WSe 2

Semiconductor moiré superlattices have been shown to host a wide array of interaction-driven ground states. However, twisted homobilayers have been difficult to study in the limit of large moiré wavelengths, where interactions are most dominant. In this study, we conducted local electronic compressibility measurements of twisted bilayer WSe 2 (tWSe 2 ) at small twist angles. We demonstrated multiple topological bands that host a series of Chern insulators at zero magnetic field near a “magic angle” around 1.23°. Using a locally applied electric field, we induced a topological quantum-phase transition at one hole per moiré unit cell. Finally, our work establishes the topological phase diagram of a generalized Kane-Mele-Hubbard model in tWSe 2 , demonstrating a tunable platform for strongly correlated topological phases.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Grassmannian diffusion maps based surrogate modeling via geometric harmonics

Abstract A novel surrogate model based on the Grassmannian diffusion maps (GDMaps) and utilizing geometric harmonics (GH) is developed for predicting the response of complex physical phenomena. The method utilizes GDMaps to obtain a low‐dimensional representation of the underlying behavior of physical/mathematical systems with respect to uncertain input parameters. Using this representation, GH, an out‐of‐sample extension technique, is employed to create a global map from the input parameter space to a Grassmannian diffusion manifold. GH is further employed to locally map points on the diffusion manifold onto the tangent space of a Grassmann manifold. The exponential map is then used to project the points in the tangent space onto the Grassmann manifold, where reconstruction of the full solution is performed. The performance of the proposed surrogate model is verified with three examples. The first problem is a toy example used to illustrate the technique. In the second example, errors associated with the various mappings are assessed by studying response predictions of the electric potential of a dielectric cylinder in a homogeneous electric field. The last example applies the method for uncertainty prediction in the strain field evolution in a model amorphous material using the shear transformation zone theory of plasticity.

42 ENGINEERING↗
DOE OSTI2022

Relevance of E × B drifts for particle and heat transport in divertors

Radial electric fields up to ~4 kV m –1 are observed in the boundary between the private flux region (PFR) and the scrape-off layer (SOL) driving E × B drifts between the inner and outer targets at speeds up to 2.8 km s –1 in the Tokamak à configuration variable divertor. The resulting E × B fluxes, located in a narrow region ($\Delta {\rho _\Psi } < 0.012$ in normalized radius or $\Delta $ R – R sep <4 mm mapped to the outer midplane) are equivalent to around 20% of the total heat and particle flux to the divertor targets (inner + outer). At the peak E r , the E × B poloidal transport is equivalent to parallel flows with M ∥ ~ 3. In the snowflake divertor with a second X-point in the outer SOL, the drifts in the PFR-SOL boundary were equivalent to around 30% of the total heat and particle flux to the divertor targets and cover a region ~50% wider than in the single null ($\Delta {\rho _\Psi }$ ~ 0.018, $\Delta $R – R sep ~ 6 mm). Furthermore, the location of the PFR-SOL boundary drift shifts radially in the E ∥ × B direction when reversing the toroidal field direction. Peaks in density and electron pressure have been identified near the primary X-point along with large gradients in density, temperature, and potential, the latter resulting in a local electric field ~2.7 kV m –1 which drives a drift (1.9 km s –1 ) upwards towards the closed flux surfaces. Floating potential (V f ) magnitudes up to 75 V (~2 kT e ) were measured, indicating that V f and parallel currents should not be neglected when estimating plasma potential.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Assessing Low-Temperature Geothermal Play Types: Relevant Data and Play Fairway Analysis Methods

The U.S. Department of Energy (DOE) Geothermal Technologies Office (GTO) is supporting the Geothermal Heating and Cooling Geospatial Datasets and Analysis project conducted by the National Renewable Energy Laboratory (NREL) as part of a broader effort to demonstrate the multi-faceted value of integrating geothermal power and geothermal heating and cooling (GHC) technologies into national decarbonization plans and community energy plans. Currently, there is a need to establish baseline low-temperature geothermal resource datasets and evaluate methods for deploying these technologies to provide the basis for supporting private sector investment. This project is focused on collecting baseline datasets, updating conceptual models, and creating Play Fairway Analysis (PFA) workflows for low-temperature (<150 degrees Celsius) geothermal resources of different geothermal play types (i.e., sedimentary basin, orogenic belts, and radiogenic geothermal play types) that could be used for geothermal heating and cooling (GHC), combined heat and power (CHP), and other geothermal direct uses (GDU) applications. Low-temperature geothermal resources are defined as reservoirs - natural or engineered - with temperatures <150 degrees Celsius. While the focus in the NREL effort is on GHC, resources at the upper end of this temperature range can also be used for small-scale power generation. This project does not include Ground Source Heat Pumps (GSHPs) technologies because they can be effectively developed almost anywhere. Low-temperature geothermal resources have not been studied as extensively as higher- to medium-temperature geothermal resources, but there is recent interest in improving understanding of these types of resources with an uptick of interest in geothermal technologies for decarbonizing heating and cooling systems. In addition, Enhanced Geothermal Systems (EGS) and other emerging technologies for exploiting petrothermal resources have opened the possibility of utilizing deep sedimentary basin systems, where porous media provide permeability and high temperatures can be reached at great depths. This project takes the approach of classifying low- temperature geothermal resources by geothermal play type (GPT). We defined and characterized three major classes of low-temperature GPT: sedimentary basins, orogenic systems, and radiogenic systems. We develop methodologies for evaluating and analyzing the potential for these resources building off the PFA approach to de-risking geothermal exploration and characterization. The proposed PFA approach for low-temperature geothermal resources includes: 1) identifying relevant data (e.g., datasets such bottom-hole temperatures from oil and gas wells, heat flow data, Quaternary faults and stress field data, geophysical data, etc.); 2) grouping and weighting of relevant datasets into PFA criteria (e.g., geological, risk, and economic criteria); 3) uncertainty quantification; 4) developing favorability or common risk maps for low-temperature geothermal resources to identify potential locations for more focused data collection; and 5) estimating electric power generation and heating potential at those locations using the GeoRePORT Resource Size Assessment Tool (RSAT). This project should facilitate future deployment of GHC, CHP, and GDU by providing data, tools, and a workflow applicable to low-temperature geothermal resources. Increased deployment of GHC and GDU will help achieve national and local decarbonization goals.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Cross-correlation spectra in interacting quantum dot systems

Two-color spin-noise spectroscopy of interacting electron spins in singly charged semiconductor quantum dots provides information on the interquantum dot interactions. We investigate the spin cross-correlation function in a quantum dot ensemble employing a modified semiclassical approach. Spin-correlation functions are calculated using a Hamilton quaternion approach that maintains local quantum mechanical properties of the spins. This method takes into account the effects of the nuclear-electric quadrupolar interactions, the randomness of the coupling constants, and the variation of the electron g factor on the spin-noise power spectra. We demonstrate that the quantum dot ensemble can be mapped on an effective two-quantum dot problem and discuss how the characteristic length scale of the interdot interaction modifies the low-frequency cross-correlation spectrum. Here, we argue that details on the interaction strength distribution can be extracted from the cross-correlation spectrum when applying a longitudinal or a transversal external magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Atomic step disorder on polycrystalline surfaces leads to spatially inhomogeneous work functions

Structural disorder causes materials’ surface electronic properties, e.g., work function (ϕ), to vary spatially, yet it is challenging to prove exact causal relationships to underlying ensemble disorder, e.g., roughness or granularity. For polycrystalline Pt, nanoscale resolution photoemission threshold mapping reveals a spatially varying ϕ=5.70±0.03 eV over a distribution of (111) vicinal grain surfaces prepared by sputter deposition and annealing. With regard to field emission and related phenomena, e.g., vacuum arc initiation, a salient feature of the ϕ distribution is that it is skewed with a long tail to values down to 5.4 eV, i.e., far below the mean, which is exponentially impactful to field emission via the Fowler–Nordheim relation. We show that the ϕ spatial variation and distribution can be explained by ensemble variations of granular tilts and surface slopes via a Smoluchowski smoothing model wherein local ϕ variations result from spatially varying densities of electric dipole moments, intrinsic to atomic steps, that locally modify ϕ. Atomic step-terrace structure is confirmed with scanning tunneling microscopy (STM) at several locations on our surfaces, and prior works showed STM evidence for atomic step dipoles at various metal surfaces. From our model, we find an atomic step edge dipole μ=0.12 D/edge atom, which is comparable to values reported in studies that utilized other methods and materials. Our results elucidate a connection between macroscopic ϕ and the nanostructure that may contribute to the spread of reported ϕ for Pt and other surfaces and may be useful toward more complete descriptions of polycrystalline metals in the models of field emission and other related vacuum electronics phenomena, e.g., arc initiation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Photophysics and Electronic Structure of Lateral Graphene/MoS 2 and Metal/MoS 2 Junctions

Integration of semiconducting transition metal dichalcogenides (TMDs) into functional optoelectronic circuitries requires an understanding of the charge transfer across the interface between the TMD and the contacting material. In this work, we use spatially resolved photocurrent microscopy to demonstrate electronic uniformity at the epitaxial graphene/molybdenum disulfide (EG/MoS 2 ) interface. A 10× larger photocurrent is extracted at the EG/MoS 2 interface when compared to the metal (Ti/Au)/MoS 2 interface. This is supported by semi-local density functional theory (DFT), which predicts the Schottky barrier at the EG/MoS 2 interface to be ~2× lower than that at Ti/MoS 2 . We provide a direct visualization of a 2D material Schottky barrier through combination of angle-resolved photoemission spectroscopy with spatial resolution selected to be ~300 nm (nano-ARPES) and DFT calculations. A bending of ~500 meV over a length scale of ~2-3 μm in the valence band maximum of MoS 2 is observed via nano-ARPES. We explicate a correlation between experimental demonstration and theoretical predictions of barriers at graphene/TMD interfaces. Spatially resolved photocurrent mapping allows for directly visualizing the uniformity of built-in electric fields at heterostructure interfaces, providing a guide for microscopic engineering of charge transport across heterointerfaces. This simple probe-based technique also speaks directly to the 2D synthesis community to elucidate electronic uniformity at domain boundaries alongside morphological uniformity over large areas.

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

Continuum gyrokinetic simulations of edge plasmas in single-null geometries

The first continuum gyrokinetic calculations of electrostatic ion scale turbulence are presented for the case of a diverted tokamak geometry. The simulation model solves the long-wavelength limit of the full-F gyrokinetic equation for ion species coupled to the quasi-neutrality equation for electrostatic potential variations, where a fluid model is used for an electron response. Furthermore, in order to facilitate simulations of highly-anisotropic microturbulence in the presence of strong magnetic shear and a magnetic X-point, a numerical algorithm utilizing a locally field-aligned multiblock coordinate system has been developed and implemented in the 5D finite-volume code COGENT. In this approach, the toroidal direction is divided into blocks, such that within each block, the cells are field-aligned and a non-matching grid interface is allowed at block boundaries. The toroidal angle corresponds to the “coarse” field-aligned coordinate, whereas the poloidal cross section, comprised of the radial and poloidal directions, is finely gridded to resolve short-scale perpendicular turbulence structures and to support accurate re-mapping (interpolation) at block boundaries. The 5D simulations explore cross-separatrix ion scale turbulence in the presence of a self-consistent radial electric field and address the effects of magnetic-shear stabilization in the X-point region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Imaging Charged Exciton Localization in van der Waals WSe 2 /MoSe 2 Heterobilayers

Exciton localization in transition-metal dichalcogenide monolayers is behind a variety of interesting phenomena and applications, including broad-spectrum solar cells and singlephoton emissions. Strain fields at the periphery of topographically distinct features such as nanoscopic bubbles were recently associated with localized charge-neutral excitons. Here, we use tip-enhanced photoluminescence (PL) to visualize excitons in WSe 2 /MoSe 2 heterobilayers (HBL). We find strong optical emission from charged excitons, particularly positively charged trions, in HBL supported by interlayer charge transfer. Our results reveal strong trion confinement, with a localization length scale comparable to the trion size, at the apex region inside individual nanoscopic bubbles. Nano-PL mapping also shows sub-10 nm spatial variations in the localized trion emission spectra, which stem from atomic-scale potential energy fluctuations. These findings demonstrate the possibility of confining charged exciton complexes that are electrically tunable, opening up further opportunities to probe many-body exciton physics, and to explore additional possible sites for strong exciton localization that can lead to quantum emissions.

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

Exploring Local Vibrational Structure in Protein-Bound Chlorophyll a: Isotope-Enrichment Experiments and Electrostatic Analysis

Local chlorophyll (Chl) vibrations play an essential role in biological photosynthesis by facilitating rapid energy transfer between pigments. In research studies, they also provide a useful spectroscopic probe of the local protein environment that surrounds each pigment. However, measuring the complete vibrational absorption spectrum of a protein-bound Chl molecule is much more difficult than, for example, Chl in neat solvent due to overlap with protein vibrations that typically drown out Chl vibrational features. Resonance Raman and fluorescence spectroscopies provide a way around this problem for Franck–Condon active vibrations, but these often rely on cryogenic measurement conditions and fail to capture vibrational signatures from, e.g., ester-group vibrations that lack coupling to a convenient electronic transition. In the present contribution, we use 13 C-enrichment of the protein backbone to shift protein background signals to lower frequency, providing a largely clean spectral window in which to study local Chl a C═O stretch modes. A room-temperature absorption spectrum for Chl a in the water-soluble chlorophyll protein (WSCP) of Lepidium virginicum is thus extracted as a difference between protein-plus-pigment and protein-only vibrational spectra. Excellent agreement in the molecular fingerprint region with the vibrational spectrum of Chl a in organic solvents confirms that the resulting spectrum represents the response of the protein-bound Chl a molecule. Furthermore, the ester group resonance is observed to shift in response to the S53P mutation that eliminates a 17 3 ester group hydrogen bond from the protein environment. Finally, we analyze these experimental results using MD-based electrostatic analysis, finding that electric-field mapping at the C atom of the ester group provides a satisfactory explanation for the observed frequency shifts between organic solvent and protein environment. MD analysis further suggests that a red-shifted ester peak observed experimentally for the S53P mutant results from solvation of the 17 3 Chl ester group due increase water penetration into the Chl-binding pocket relative to wild-type WSCP. Here, we anticipate that these results will prove useful both for benchmarking future simulation work and as a reference for interpreting Chl vibrational spectra as a probe of pigment–protein interactions.

14 SOLAR ENERGY↗
DOE OSTI2025

Exploring Domain-Wall Pinning in Ferroelectrics via Automated High-Throughput Atomic Force Microscopy

Domain-wall dynamics in ferroelectric materials are strongly position-dependent, since each polar interface is locked into a unique local microstructure. This necessitates spatially resolved studies of wall pinning using scanning-probe microscopy techniques. The pinning centers and pre-existing domain walls are usually sparse within the image plane, precluding the use of dense hyperspectral imaging modes and requiring time-consuming human experimentation. Here, a large-area epitaxial PbTiO 3 film on cubic KTaO 3 was investigated to quantify the electric-field-driven dynamics of the polar–strain domain structures using ML-controlled automated piezoresponse force microscopy. Analysis of 1500 switching events reveals that domain-wall displacement depends not only on field parameters but also on the local ferroelectric–ferroelastic configuration. For example, twin boundaries in polydomains regions, like a 1 – /c+ ∥ a 2 – /c – , stay pinned up to a certain level of bias magnitude and change only marginally as the bias increases from 20 to 30 V, whereas single-variant boundaries, like the a 2 + /c + ∥ a 2 – /c – stack, are already activated at 20 V. These statistics on the possible ferroelectric and ferroelastic wall orientations, together with the automated high-throughput AFM workflow, can be distilled into a predictive map that links domain configurations to pulse parameters. Here, this microstructure-specific rule set forms the foundation for the design of ferroelectric memories.

automated scanning probe microscopy↗
DOE OSTI2021

Magnetic field sensors for detection of trapped flux in superconducting radio frequency cavities

Superconducting radio-frequency (SRF) cavities are fundamental building block of modern particle accelerators. They operate at liquid helium temperatures 2-4 K to achieve very high quality factors 1e10-1e11. Trapping of magnetic flux within the superconductor is a significant contribution to the residual RF losses which limit the achievable quality factor. Suitable diagnostic tools are in high demand to understand the mechanisms of flux trapping in technical superconductors and the fundamental components of such diagnostic tool are magnetic field sensors. We have studied the performance of commercially available Hall probes, anisotropic magnetoresistive sensors and flux-gate magnetometers with respect to their sensitivity and capability to detect localized, low magnetic flux amplitudes, of the order of few tens of magnetic flux quantum at liquid helium temperature. Although Hall probes have the lowest magnetic field sensitivity (~96$ nV$/microT at 2 K), their physical dimensions are such that they have the ability to detect the lowest number of trapped vortices among the three type of sensors. Hall probes and anisotropic magnetoresistive sensors have been selected to be used in a setup to map regions of trapped flux on the surface of a single-cell SRF cavity.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

A Synoptic VLBI Technique for Localizing Nonrepeating Fast Radio Bursts with CHIME/FRB

We demonstrate the blind interferometric detection and localization of two fast radio bursts (FRBs) with subarcminute precision on the 400 m baseline between the Canadian Hydrogen Intensity Mapping Experiment (CHIME) and the CHIME Pathfinder. In the same spirit as Very Long Baseline Interferometry (VLBI), the telescopes were synchronized to separate clocks, and the channelized voltage (herein referred to as baseband) data were saved to a disk with correlation performed offline. The simultaneous wide field of view and high sensitivity required for blind FRB searches implies a high data rate—6.5 terabits per second (Tb/s) for CHIME and 0.8 Tb s{sup −1} for the Pathfinder. Since such high data rates cannot be continuously saved, we buffer data from both telescopes locally in memory for ≈40 s, and write to the disk upon receipt of a low-latency trigger from the CHIME Fast Radio Burst Instrument (CHIME/FRB). The ≈200 deg{sup 2} field of view of the two telescopes allows us to use in-field calibrators to synchronize the two telescopes without needing either separate calibrator observations or an atomic timing standard. In addition to our FRB observations, we analyze bright single pulses from the pulsars B0329+54 and B0355+54 to characterize systematic localization errors. Our results demonstrate the successful implementation of key software, triggering, and calibration challenges for CHIME/FRB Outriggers: cylindrical VLBI outrigger telescopes which, along with the CHIME telescope, will localize thousands of single FRB events with sufficient precision to unambiguously associate a host galaxy with each burst.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Geostatistical Mapping of Salinity Conditioned on Borehole Logs, Montebello Oil Field, California

We present a geostatistics-based stochastic salinity estimation framework for the Montebello Oil Field that capitalizes on available total dissolved solids (TDS) data from groundwater samples as well as electrical resistivity (ER) data from borehole logging. Data from TDS samples (n = 4924) was coded into an indicator framework based on falling below four selected thresholds (500, 1000, 3000, and 10,000 mg/L). Collocated TDS-ER data from the surrounding groundwater basin were then employed to produce a kernel density estimator to establish conditional probabilities for ER data (n = 8 boreholes) falling below the selected TDS thresholds within the Montebello Oil Field area. Directional variograms were estimated from these indicator coded data, and 500 TDS realizations from conditional indicator simulation were generated for the subsurface region above the Montebello Oil Field reservoir. Simulations were summarized as 3D maps of median TDS, most likely salinity class, and probability for exceeding each of the specified TDS thresholds. Results suggested TDS was below 500 mg/L in most of the study area, with a trend toward higher values (500 to 1000 mg/L) to the southwest; consistent with the average regional groundwater flow direction. Discrete localized zones of TDS greater than 1000 mg/L were observed, with one of these zones in the greater than 10,000 mg/L range; however, these areas were not prevalent. The probabilistic approach used here is adaptable and is readily modified to include additional data and types and can be employed in time-lapse salinity modeling through Bayesian updating.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Nano-photocurrent Mapping of Local Electronic Structure in Twisted Bilayer Graphene

In this work, we report a combined nano-photocurrent and infrared nanoscopy study of twisted bilayer graphene (TBG) enabling access to the local electronic phenomena at length scales as short as 20 nm. We show that the photocurrent changes sign at carrier densities tracking the local superlattice density of states of TBG. We use this property to identify domains of varying local twist angle by local photothermoelectric effect. Consistent with the photocurrent study, infrared nanoimaging experiments reveal optical conductivity features dominated by twist-angle-dependent interband transitions. Our results provide a fast and robust method for mapping the electronic structure of TBG and suggest that similar methods can be broadly applied to probe electronic inhomogeneities of Moiré superlattices in other van der Waals heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Deep learning to estimate permeability using geophysical data

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

58 GEOSCIENCES↗