Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed state estimation”

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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The High Fraction of Thin Disk Galaxies Continues to Challenge ΛCDM Cosmology

Any viable cosmological framework has to match the observed proportion of early- and late-type galaxies. In this contribution, we focus on the distribution of galaxy morphological types in the standard model of cosmology (Lambda cold dark matter, ΛCDM). Using the latest state-of-the-art cosmological ΛCDM simulations known as Illustris, IllustrisTNG, and EAGLE, we calculate the intrinsic and sky-projected aspect ratio distribution of the stars in subhalos with stellar mass M * > 10 10 M ⊙ at redshift z = 0. There is a significant deficit of intrinsically thin disk galaxies, which however comprise most of the locally observed galaxy population. Consequently, the sky-projected aspect ratio distribution produced by these ΛCDM simulations disagrees with the Galaxy And Mass Assembly (GAMA) survey and Sloan Digital Sky Survey at ≥12.52σ (TNG50-1) and ≥14.82σ (EAGLE50) confidence. The deficit of intrinsically thin galaxies could be due to a much less hierarchical merger-driven build-up of observed galaxies than is given by the ΛCDM framework. It might also arise from the implemented sub-grid models, or from the limited resolution of the above-mentioned hydrodynamical simulations. We estimate that an 8 5 times better mass resolution realization than TNG50-1 would reduce the tension with GAMA to the 5.58σ level. Finally, we show that galaxies with fewer major mergers have a somewhat thinner aspect ratio distribution. Given also the high expected frequency of minor mergers in ΛCDM, the problem may be due to minor mergers. In this case, the angular momentum problem could be alleviated in Milgromian dynamics because of a reduced merger frequency arising from the absence of dynamical friction between extended dark matter halos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Elliptically-Contoured Tensor-variate Distributions with Application to Image Learning

Statistical analysis of tensor-valued data has largely used the tensor-variate normal (TVN) distribution that may be inadequate for data arising from distributions with heavier or lighter tails. We study a general family of elliptically contoured (EC) TV distributions and derive its characterizations, moments, marginal, and conditional distributions. We describe procedures for maximum likelihood estimation from data that are (1) uncorrelated draws from an EC distribution, (2) from a scale mixture of the TVN distribution, and (3) from an underlying but unknown EC distribution, for which we extend Tyler’s robust estimator. A detailed simulation study highlights the benefits of choosing an EC distribution over the TVN for heavier-tailed data. We develop TV classification rules using discriminant analysis and EC errors and show that they better predict cats and dogs from images in the Animal Faces-HQ dataset than the TVN-based rules. A novel tensor-on-tensor regression and TV analysis of variance (TANOVA) framework under EC errors is also demonstrated to better characterize gender, age, and ethnic origin than the usual TVN-based TANOVA in the celebrated labeled faces of the wild dataset.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2019

Assessing Dual-Polarization Radar Estimates of Extreme Rainfall During Hurricane Harvey

Hurricane Harvey hit the Texas Gulf Coast as a major hurricane on 25 August 2017 before exiting the state as a tropical storm on 29 August 2017. Left in its wake was historic flooding, with some locations measuring more than 60 in. (150 cm) of rain over a 5-day period. The WSR-88D radar (KHGX) maintained operations for the entirety of the event. Rain gauge data from the Harris County Flood Warning System (HCFWS) was used for validation with the full radar dataset to retrieve daily and event-total precipitation estimates for the period 25–29 August2017. The KHGX precipitation estimates were then compared with the HCFWS gauges. Three different hybrid polarimetric rainfall retrievals were used, along with attenuation-based retrieval that employs the radar-observed differential propagation. An advantage of using a attenuation-based retrieval is its immunity to partial beam blockage and calibration errors in reflectivity and differential reflectivity. All of the retrievals are susceptible to changes in the observed drop size distribution (DSD). No in situ DSD data were available over the study area, so changes in the DSD were interpreted by examining the observed radar data. We examined the parameter space of two key values in the attenuation retrieval to test the sensitivity of the rain retrieval. Selecting a value ofa50.015 andb50.600 provided the best overall results, relative to the gauges, but more work needs to be done to develop an automated technique to account for changes in the ambient DSD.

David Billings Wolff↗
DOE OSTI2024

Tracking the Sun: Pricing and Design Trends for Distributed Photovoltaic Systems in the United States, 2024 Edition [Slides]

Berkeley Lab’s annual Tracking the Sun report describes trends among grid-connected, distributed solar photovoltaic (PV) and paired PV+storage systems in the United States. For the purpose of this report, distributed solar includes residential systems, roof-mounted non-residential systems, and ground-mounted systems up to 5 MW-AC. Ground-mounted systems larger than 5 MW-AC are covered in Berkeley Lab’s companion annual report, Utility-Scale Solar. The latest edition of the report is based on 3.7 million systems installed through year-end 2023, representing close to 80% of systems installed to date. The report describes and discusses key trends related to: -Project characteristics, including system size, module efficiencies, roof-coverage ratios, prevalence of paired PV with storage, use of module-level power electronics, third-party ownership, mounting configurations, panel orientation, and customer segmentation -Median installed-price trends, both nationally and by state -Variability in pricing according to system size, state, installer, equipment type, and other factors, relying on both descriptive analysis and a multi-variate regression to estimate the effects of key pricing drivers for residential systems installed in 2023.

14 SOLAR ENERGY↗
NASA NTRS2011

An Estimate of the North Atlantic Basin Tropical Cyclone Activity for the 2011 Hurricane Season

Estimates are presented for the expected level of tropical cyclone activity for the 2011 North Atlantic Basin hurricane season. It is anticipated that the frequency of tropical cyclones for the North Atlantic Basin during the 2011 hurricane season will be near to above the post-1995 means. Based on the Poisson distribution of tropical cyclone frequencies for the current more active interval 1995-2010, one computes P(r) = 63.7% for the expected frequency of the number of tropical cyclones during the 2011 hurricane season to be 14 plus or minus 3; P(r) = 62.4% for the expected frequency of the number of hurricanes to be 8 plus or minus 2; P(r) = 79.3% for the expected frequency of the number of major hurricanes to be 3 plus or minus 2; and P(r) = 72.5% for the expected frequency of the number of strikes by a hurricane along the coastline of the United States to be 1 plus or minus 1. Because El Nino is not expected to recur during the 2011 hurricane season, clearly, the possibility exists that these seasonal frequencies could easily be exceeded. Also examined are the effects of the El Nino-Southern Oscillation phase and climatic change (global warming) on tropical cyclone seasonal frequencies, the variation of the seasonal centroid (latitude and longitude) location of tropical cyclone onsets, and the variation of the seasonal peak wind speed and lowest pressure for tropical cyclones.

Wilson, Robert M.↗
NASA NTRS2019

Assessing Dual-Polarization Radar Estimates of Extreme Rainfall During Hurricane Harvey

Hurricane Harvey hit the Texas Gulf Coast as a major hurricane on August 25, 2017 before exiting the state as a tropical storm on August 29, 2017. Left in its wake was historic flooding, with some locations measuring more than 60 inches of rain over a five-day period. The WSR-88D radar (KHGX) maintained operations for the entirety of the event. Rain gauge data from the Harris County Flood Warning System (HCFWS) was used for validation with the full radar data set to retrieve daily and event-total precipitation estimates for the period August 25-29, 2017. The KHGX precipitation estimates were then compared to the HCFWS gauges. Three different hybrid polarimetric rainfall retrievals were used, along with attenuation-based retrieval that employs the radar-observed differential propagation. An advantage of using a attenuation-based retrieval is its immunity to partial beam blockage and calibration errors in reflectivity and differential reflectivity. All of the retrievals are susceptible to changes in the observed Drop Size Distribution (DSD). No in situ DSD data were available over the study area, so changes in the DSD were interpreted by examining the observed radar data. We examined the parameter space of two key values in the attenuation retrieval to test the sensitivity of the rain retrieval. Selecting a value of α=0.015 and β=0.600 and β=0.625 provided the best overall results, relative to the gauges, but more work needs to be done to develop an automated technique to account for changes in the ambient DSD.

David B Wolff↗
NASA NTRSDate not supplied

Flux Comparison of Master-8 and Ordem 3.1 Modelled Space Debris Population

With ESA’s Meteoroid And Space debris Terrestrial Environment Reference (MASTER-8) model and NASA’s Orbital Debris Engineering Model(ORDEM) 3.1, the two premier orbital debris engineering models have been officially released. The two models come with significant enhancements and now represent state-of-the-art orbital debris modelling for their respective agencies. Both models provide the community with estimates of the space debris environment from low Earth orbit (LEO) up to at least geostationary altitude. The MASTER population is an event-based simulation of all known events that generate debris and objects that are part of the U.S. Space Surveillance Network (SSN) catalog, which provides coverage of objects with diameters down to approximately 10cm in LEO and 1 m in geosynchronous Earth orbit (GEO). Different models are used to simulate the artificial objects and their orbit evolution over time. These models are called “sources” since they assign an origin to each individual object and consist of fragments, solid rocket motor (SRM) remainders, sodium-potassium(NaK) droplets, paintflakes, ejecta, and multi-layer insulation (MLI) fragments. The objects from each source are characterized by having individual release mechanisms, as well as orbital distributions, material composition, size, and mass distributions. Dedicated radar and telescope observation data is used to calibrate the model for objects larger than 1 cm in LEO and larger than 10 cm in GEO. For calibrating the small-sized objects, below 1 cm, impact data from returned surfaces are analyzed. Because the >1cm object population is dominated by fragments, the fragmentation event database was updated to include new events, as well as re-evaluate past events. Special attention was drawn to re-evaluating theFengyun-1C anti-satellite test from 2007 and Cosmos-Iridium collision event from 2009 since these events shape the fragment population because of their severity. After 2009, the two largest fragmentations in terms of number of tracked debris are the Briz-M explosion in 2012 and the NOAA-16explosion in 2015. In total, there are 261 confirmed fragmentations in the database up to November 2016.The baseline population for ORDEM 3.1 is based on the U.S. SSN catalog, and observational datasets from radar, in situ, and optical sources provide a foundation from which the model populations are statistically extrapolated to smaller size regions. These regions are not well-covered by the SSN catalog yet may pose the greatest threat to operational spacecraft. The NASA Standard Satellite Breakup Model is used to generate fragments greater than 1mm from collisions and explosions, and these fragment populations are scaled using ground-based radar data. Specific major debris-producing events, including the Fengyun-1C, Iridium 33, and Cosmos2251 debris clouds, and unique populations, such as NaK droplets, were re-examined, modelled, and added to the ORDEM environment separately. Optical measurement data is used to model the GEO population down to 10 cm. The debris environment is propagated using NASA’s LEO-to-GEO Environment Debris model, and future explosions of intact objects and collisions involving objects greater than 10 cm are assessed statistically. The environment from a few millimetres down to 10 𝜇m is modelled using a special degradation model where small particles are generated from intact spacecraft and rocket bodies, then the populations are scaled to fit in situ cratering data from Space Shuttle returned surfaces. Fragments smaller than 10 cm are differentiated based on material density categories, i.e., high-, medium-, and low-density, to better characterize the potential debris risk posed to upper stages and spacecraft. This paper will discuss the MASTER and ORDEM approaches for modelling populations and compare fluxes for specific orbits, including sun-synchronous, ISS-altitude, geosynchronous transfer, and GEO. In the end, a conclusion is drawn towards the importance of having multiple fundamentally different, yet validated, models to estimate the space debris population.

Andre Horstmann↗
DOE OSTI2024

Environmental associations of Ophidiomyces ophidiicola , the causative agent of ophidiomycosis in snakes

Emerging pathogenic fungi have become a topic of conservation concern due to declines observed in several host taxa. One emerging fungal pathogen, Ophidiomyces ophidiicola, is well documented as the causative agent of ophidiomycosis, otherwise known as snake fungal disease (SFD). O. ophidiicola has been found to cause disease in a variety of snake species across the United States, including the eastern massasauga (Sistrurus catenatus), a federally threatened rattlesnake species. Most work to date has involved detecting O. ophidiicola for diagnosis of infection through direct sampling of snakes, and attempts to detect O. ophidiicola in the abiotic environment to better understand its distribution, seasonality, and habitat associations are lacking. We collected topsoil and groundwater samples from four macrohabitat types across multiple seasons in northern Michigan at a site where Ophidiomyces infection has been confirmed in eastern massasauga. Using a quantitative PCR (qPCR) assay developed for O. ophidiicola, we detected Ophidiomyces DNA in topsoil but observed minimal to no detection in groundwater samples. Detection frequency did not differ between habitats, but samples grouped seasonally showed higher detection during mid-summer. We found no relationships of detection with hypothesized environmental correlates such as soil pH, temperature, or moisture content. Furthermore, the distribution of Ophidiomyces positive samples across the site was not linked to estimated space use of massasaugas. Our data suggests that season has some effect on the presence of Ophidiomyces. Differences in presence between habitats may exist but are likely more dependent on the time of sampling and currently uninvestigated soil or biotic parameters. These findings build on our understanding of Ophidiomyces ecology and epidemiology to help inform where and when snakes may be exposed to the fungus in the environment.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2013

The GLAS Algorithm Theoretical Basis Document for Precision Orbit Determination (POD)

The Geoscience Laser Altimeter System (GLAS) was the sole instrument for NASA's Ice, Cloud and land Elevation Satellite (ICESat) laser altimetry mission. The primary purpose of the ICESat mission was to make ice sheet elevation measurements of the polar regions. Additional goals were to measure the global distribution of clouds and aerosols and to map sea ice, land topography and vegetation. ICESat was the benchmark Earth Observing System (EOS) mission to be used to determine the mass balance of the ice sheets, as well as for providing cloud property information, especially for stratospheric clouds common over polar areas. The GLAS instrument operated from 2003 to 2009 and provided multi-year elevation data needed to determine changes in sea ice freeboard, land topography and vegetation around the globe, in addition to elevation changes of the Greenland and Antarctic ice sheets. This document describes the Precision Orbit Determination (POD) algorithm for the ICESat mission. The problem of determining an accurate ephemeris for an orbiting satellite involves estimating the position and velocity of the satellite from a sequence of observations. The ICESatGLAS elevation measurements must be very accurately geolocated, combining precise orbit information with precision pointing information. The ICESat mission POD requirement states that the position of the instrument should be determined with an accuracy of 5 and 20 cm (1-s) in radial and horizontal components, respectively, to meet the science requirements for determining elevation change.

Rim, Hyung Jin↗
DOE OSTI2021

Repeal of the Clean Power Plan: Social Cost and Distributional Implications

The Clean Power Plan (CPP) was repealed due to concerns about the “unnecessary, costly burdens” it may impose on electric utilities, thereby delaying efforts to reduce carbon dioxide emissions (CO 2 ) from the electricity sector. This paper examines the greenhouse gas and welfare implications of this repeal while incorporating the presence of the state renewable portfolio standards (RPS) in the US as the status quo. We assess the carbon abatement and welfare costs with the CPP relative to two alternative baselines: a no–policy baseline and a pre–existing policy baseline with the RPS. The CPP is implemented as a regional mass–based standard, a regional rate–based standard, or as a national mass–based standard with trading of emissions across regions over the 2022–2030 period. We find that the incremental discounted welfare costs per metric ton of CO 2 that would have been abated by the CPP relative to the RPS would be substantially lower than the global social cost of CO 2 . However, the overall costs of carbon abatement with the CPP added to the RPS would have become higher than the social cost of carbon when estimated relative to a no–policy baseline, except with a national mass–based CPP. Across all policy combinations and choice of baselines, the aggregate welfare costs were lowest under a national mass–based standard and highest under the regional rate–based standard. Here, we also find that the CPP would have imposed large welfare costs on consumers and fossil fuel producers while benefiting the renewable fuel producers.

, Q42↗
DOE OSTI2022

Laser diagnostics to characterize the in-flame growth of platinum nanoparticles manufactured by the reactive spray deposition technology

Reactive Spray Deposition Technology (RSDT) is an atmospheric pressure, flame-based advanced manufacturing method used to fabricate Membrane Electrode Assemblies (MEAs) for proton exchange fuel cells and water electrolyzers. RSDT combines the flame synthesis of catalyst nanoparticles and their deposition onto the membrane in one step. The properties of the synthesized nanoparticles, such as their Size Distribution Function (SDF), determine the manufactured electrodes’ performance, which can be evaluated when the deposition is complete via ex-situ characterization of the products. Efforts to improve RSDT for manufacturing state-of-the-art MEAs can be significantly enhanced with integrated laser diagnostics that enable in situ measurement of the synthesized catalyst nanoparticles. This paper reports the implementation of laser diagnostics in an RSDT facility and evaluates their potential to assist the manufacturing process. Laser Light Scattering (LS) and Laser-Induced Incandescence (LII) measurements are performed in the oxygen-rich zone of the flame at various distances from the flame fuel jet nozzle downstream of its luminous region. The measurements quantitatively track the evolution in size and volume fraction of platinum-based nanoparticles in two flames that yield different catalytic properties in the manufactured electrodes. The profiles of the measured nanoparticles’ volume fraction along the flame axis can be estimated a priori so that the average nanoparticle sizes can be measured in quasi-real time via LS. Nanoparticles experience an extremely slow growth rate while being convected at distances from 150 nm to 300 nm from the fuel nozzle. Complementary characterizations of the synthesized nanoparticles are performed ex-situ via High-Angle Annular Darkfield Scanning Transmission Electron Microscopy (HAADF-STEM) image analyses of samples collected on grids at a fixed distance from the fuel jet nozzle. Comparing the results from laser and microscopy techniques not only cross-validates the findings but also provides the parameters to infer the bimodal SDF in-situ and yields evidence that the coagulation efficiency of the synthesized nanoparticles has extremely low values in the investigated zone of the flames.

Rayleigh Light Scattering (LS)↗
DOE OSTI2026

Evaluation of chromatographic resins for high efficacy separation of radium from chemical analogues

The separation of Radium (Ra) from interfering chemical analogues is of major importance for both the accurate measurement of Ra-226 in environmental samples and the preparation of radium targets destined for medical isotope production. The adsorption behavior of Ra, Barium (Ba), and Strontium (Sr) was investigated on multiple crown-ether-based proprietary extraction chromatographic resins developed by TrisKem Int. Among these, the TK101 resin, characterized by a novel room-temperature ionic liquid (RTIL) dual extractant–solvent system, demonstrated enhanced affinity and selectivity, particularly at low acid concentrations, achieving a D w of 3477 ± 73 for Ra in 0.04 M HNO₃, whereas the maximum Ra D w obtained with the traditional Sr resin in 8 M HNO 3 was only 12.33 ± 0.4. In 0.04 M HNO 3 , a Ra distribution ratio (D w ) of 776 was achieved within 1 min, indicating rapid uptake kinetics. The adsorption capacities were estimated at 19.6 ± 0.1 mg/g for Ba and 12.0 ± 0.2 mg/g for Sr. However, competing alkali metals significantly reduced alkaline earth metal uptake, indicating that the TK101 resin is best suited for radium purification following a preconcentration step to minimize matrix effects in high salinity samples. Dynamic column studies confirmed the strong chromatographic performance of the TK101 resin, with no significant Ba or Sr breakthrough in the Ra fraction and a Ra recovery of 87 ± 2.8%. These results demonstrate the potential of the TK101 resin for selective radium separation in analytical and radiochemical applications.

and nuclear chemistry↗
DOE OSTI2023

Granger Causality for prediction in Dynamic Mode Decomposition: Application to power systems

Here, the dynamic mode decomposition (DMD) technique extracts the dominant modes characterizing the innate dynamical behavior of the system within the measurement data. For appropriate identification of dominant modes from the measurement data, the DMD algorithm necessitates ensuring the quality of the input measurement data sequences. On that account, for validating the usability of the dataset for the DMD algorithm, the paper proposed two conditions: Persistence of excitation (PE) and the Granger Causality Test (GCT). The virtual data sequences are designed with the hankel matrix representation such that the dimensions of the subspace spanning the essential system modes are increased with the addition of new state variables. The PE condition provides the lower bound for the trajectory length, and the GCT provides the order of the model. Satisfying the PE condition enables estimating an approximate linear model, but the predictability with the identified model is only assured with the temporal causation among data searched with GCT. The proposed methodology is validated with the application for coherency identification (CI) in a multi-machine power system (MMPS), an essential phenomenon in transient stability analysis. The significance of PE condition and GCT is demonstrated through various case studies implemented on 22 bus six generator system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1991

Models of recurrent strike-slip earthquake cycles and the state of crustal stress

Numerical models of the strike-slip earthquake cycle, assuming a viscoelastic asthenosphere coupling model, are examined. The time-dependent simulations incorporate a stress-driven fault, which leads to tectonic stress fields and earthquake recurrence histories that are mutually consistent. Single-fault simulations with constant far-field plate motion lead to a nearly periodic earthquake cycle and a distinctive spatial distribution of crustal shear stress. The predicted stress distribution includes a local minimum in stress at depths less than typical seismogenic depths. The width of this stress 'trough' depends on the magnitude of crustal stress relative to asthenospheric drag stresses. The models further predict a local near-fault stress maximum at greater depths, sustained by the cyclic transfer of strain from the elastic crust to the ductile asthenosphere. Models incorporating both low-stress and high-stress fault strength assumptions are examined, under Newtonian and non-Newtonian rheology assumptions. Model results suggest a preference for low-stress (a shear stress level of about 10 MPa) fault models, in agreement with previous estimates based on heat flow measurements and other stress indicators.

Lyzenga, Gregory A.↗
DOE OSTI2025

Evaluating switch lifetime in soft-switched single-stage differential-mode SST

The reliability of semiconductor switches in single-stage differential-mode solid-state transformers (DM-SSTs) has not been systematically evaluated under soft-switching operation and realistic grid conditions. This paper presents a switch-level reliability analysis for soft-switched and hard-switched DM-SST configurations by integrating converter-specific power loss modeling with empirical lifetime prediction. Analytical derivation of device current profiles specific to the DM-SST is used to characterize electrothermal stress, which is then mapped to lifetime using degradation models obtained from power cycling tests (PCTs). Applied to realistic SST load profiles and grid voltage variations, this approach provides a probabilistic prediction of switch lifetime for the DM-SST. Lifetime estimates for both SiC MOSFETs and Si IGBTs are presented, offering insight into device degradation under converter operating conditions. The results quantify the reliability benefits of soft switching in single-stage SSTs, highlighting how switching dynamics influence long-term switch degradation.

14 SOLAR ENERGY↗
NASA NTRS2018

A Singular Value Decomposition Framework for Retrievals with Vertical Distribution Information from Greenhouse Gas Column Absorption Spectroscopy Measurements

We review the singular value decomposition (SVD) framework and use it for quantifying and discerning vertical information in greenhouse gas retrievals from column integrated absorption measurements. While the commonly used traditional Bayesian optimal estimation (OE) assumes a prior distribution in order to regularize the inversion problem, the SVD approach identifies principal components that can be retrieved from the measurement without explicitly specifying a prior mean and prior covariance matrix. We review the SVD method, explicitly recognize the use of an uninformative prior and show it to incur no bias from the choice of the prior. We also make the connection between the SVD method and the pseudo-inverse, which makes it more intuitive and easy to understand. We illustrate the use of the SVD method on an integrated path differential absorption CO2 lidar measurement model and verify our derivations and bias-free properties versus optimal estimation using numerical simulations. In contrast, traditional OE retrievals exhibit bias when the prior mean used in the retrieval differs from the true mean. Hence, the SVD method is particularly useful for situations in which knowledge of the prior mean and prior covariance of the true state (e.g., greenhouse gas profiles) is inadequate.

Ramanathan, Anand K.↗
DOE OSTI2026

Hourly natural gas usage surrounding NIST NEB and NWB GHG monitoring stations in 2023

This dataset includes hourly natural gas usage surrounding two National Institute of Standards and Technology (NIST) greenhouse gas (GHG) monitoring stations, Northwest Baltimore (NWB; 39.3445°N, 76.6851°W) and Northeast Baltimore (NEB; 39.3154°N, 76.5830°W), for the year 2023 in units of therms. The gas usage was provided by the local gas distribution company and includes hourly data averaged across groups of 15 or more addresses to maintain anonymity. The hourly data was averaged spatially within a radius of 1 km from each monitoring station, including all groups containing data from any address within 1km of the monitoring station. This dataset only includes gas usage from service points equipped with advanced metering infrastructure (AMI). Between the two monitoring sites, this dataset includes usage from a total of 3649 service points, which we estimate is at least 44% of the service points in the domain.

Kenion, Helen C. R. [School of Environment and Sus↗
DOE OSTI2023

Real-Time Monitoring of Fracture Dynamics with a Contrast Agent-Assisted Electromagnetic Method

In collaboration with the Advanced Energy Consortium, our team has previously demonstrated that the placement of electrically active proppants (EAPs) in a hydraulic fracture surveyed by electromagnetic (EM) methods can enhance the imaging of the stimulated reservoir volumes during hydraulic fracturing. That work culminated in constructing a well-characterized EAP-filled fracture anomaly at the Devine field pilot site (DFPS). In subsequent laboratory studies, we observed that the electrical conductivity of our EAP correlates with changes in pressure, salinity, and flow. Thus, we postulated that the EAP could be used as an in-situ sensor for the remote monitoring of these changes in previously EAP-filled fractures. This paper presents our latest field data from the DFPS to demonstrate such correlations at an intermediate pilot scale. We conducted surface-based EM surveys during freshwater (200 ppm) and saltwater (2,500 ppm) slug injections while running surfaced-based EM surveys. Simultaneously, we measured the following: 1) bottomhole pressure and salinity in five monitoring wells; 2) injection rate using high-precision data loggers; 3) distributed acoustic sensors in four monitoring wells; and 4) tiltmeter data on the survey area. We demonstrated that injections into an EAP-filled fracture could be successfully coupled with real-time electric field measurements on the surface, leading to remote monitoring of dynamic changes within the EAP-filled fracture. Furthermore, by comparing the electrical field traces with the bottomhole pressure, flow rate, and salinity, we concluded that the observed electric field in our study is influenced by fracture dilation and flow rate. Salinity effect was observed when saltwater was injected. EM simulations solely based on assumptions of fracture conductivity changes during injection did not reproduce all of the measured electric field magnitudes. Preliminary estimates showed that including streaming potential in our geophysical model may be needed to reduce the simulation mismatch. The methods developed and demonstrated during this study will lead to a better understanding of the extent of fracture networks, formation stress states, fluid leakoff and invasion, characterizations of engineered fracture systems, and other applications where monitoring subsurface flow tracking is deemed important.

02 PETROLEUM↗