Search NASA⌕ Search

SEARCH · Search NASA

Results for “river inversion”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

River Plume Liftoff Dynamics and Surface Expressions

The water surface expression of liftoff and its dependence on discharge are examined using numerical simulations with the Regional Ocean Modeling System (ROMS). Liftoff is the process by which buoyant river water separates from the bed and flows over denser saltwater. During low-discharge conditions liftoff occurs in the river and is accompanied by a change in the surface slope. During high-discharge conditions liftoff occurs outside the mouth and generates a ridge on the water surface. The location and height of the ridge can be described by analytical equations in terms of discharge, shelf slope, and river mouth aspect ratio. The offshore distance and height of the ridge are proportional to the river discharge and vary inversely with river mouth aspect ratio. For steep shelf slopes liftoff occurs close to the river mouth and generates a large ridge. The ridge is modified, but not eliminated, by the presence of tides. The water surface slope change at the ridge peak is large enough to be detected by the upcoming Surface Water and Ocean Topography (SWOT) altimeter and can be used to identify the liftoff location during high discharge. However, during low discharge the water surface slope change at the liftoff location is too small to be detected by SWOT. These results indicate that remote measurements of the presence or absence of the ridge may be useful to distinguish between low and high flows, and remote measurements of the ridge location or height could be used to estimate freshwater discharge.

42 ENGINEERING↗
DOE OSTI2021

Organic sulfur fluxes and geomorphic control of sulfur isotope ratios in rivers

Pyrite oxidation plays a critical role in the relationship between weathering and climate, and its impact on the global carbon cycle has previously been constrained through inversion models utilizing observations of river sulfate ($SO^{2–}_{4}$) and its 34 S/ 32 S isotope ratio (δ 34 S SO4 ). However, measurements from some rivers have suggested that SSO4 can be substantially impacted by processes such as microbial sulfate reduction and/or sulfur assimilation and cycling, rather than simply reflecting a weighted mixture of lithologic sulfur sources. To study the prevalence and controls on $SO^{2–}_{4}$ transformations, in this study we measured dissolved major element concentrations and δ 34 S SO4 in river water samples from throughout western Iceland. Our analyses focused on samples from a small catchment hosting the Efri Haukadalsá river, a system with relatively uniform and isotopically constrained basaltic bedrock. We also measured sediment δ 34 S and sulfur speciation using sulfur K-edge X-ray absorption spectroscopy on sediment and vegetation samples from this catchment. Values of dissolved δ 34 S SO4 in the Efri Haukadalsá ranged from 2.5‰ to 23.7‰ and had a linear relationship with Cl – /$SO^{2–}_{4}$ ratios, indicating that $SO^{2–}_{4}$ predominantly derived from basalt weathering and meteoric precipitation. The lower δ 34 S SO4 values were found in fluvial valleys with V-shaped cross sections, while higher values of δ 34 S SO4 occurred in U-shaped, glacially eroded valleys with thick alluvial fills blanketing the valley floor. Spectroscopic observations identified organic sulfur phases in suspended river sediment, floodplain deposits, and vegetation. Mass balance calculations quantified the organic sulfur flux as less than 10% of $SO^{2–}_{4}$ export, and sediment δ 34 S values were comparable to river δ 34 S SO4 . We interpreted these isotopic and chemical patterns as reflecting differences in the availability of unweathered bedrock across the Efri Haukadalsá catchment, with V-shaped valleys having greater access to fresh sulfide-bearing minerals than alluviated U-shaped valleys; this interpretation is in contrast to one in which the elevated δ 34 S SO4 values reflect fractionation during sulfur transformations along alluvial reaches. These results validated the application of river inversion models for constraining weathering fluxes and affirmed that pyrite oxidation globally, even in the presence of river sulfur cycling, modulates the abundance of atmospheric carbon dioxide.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Boundary Layer Height Variation over the Savannah River Site during a Heat Wave

With the number of heat waves expect to increase in both in their severity and frequency in the coming years, an insight into their impact on planetary boundary layer (PBL) heights is essential for one development of future emergency and safety plans. Using High-Resolution Rapid Refresh (HRRR) atmospheric model output, an analysis is done on the effects a heat wave has on PBL heights and temperature inversions across the Savannah River Site (SRS). This study showed that due to its adjacency to anthropogenic land use, SRS experiences a non-negligible gradient in its PBL height, with higher heights to the northwest. The singular location PBL heigh measurement and forecasts for SRS do not account for this gradient, which has potential implications on pollutant and smoke dispersion onsite. During a heat wave this gradient is present but PBL heights are less impacted than surrounding developed areas, possibly attributable to the woody wetlands that surround the Savannah River's banks. Additionally, the heat wave appears to strongly enhance the nocturnal temperature inversions that the site and surrounding areas experience. This is possibly due to increased stability from higher pressure conditions, which are often associated with heat waves. This influence on nocturnal inversions, however, is limited to their intensity, leaving their timing relatively unaltered. The impact a heat wave has on PBL heights and temperature inversions across SRS is seen to be more dependent on adjacent land use and features rather than the unique land use types onsite.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Self‐Potential Tomography Preconditioned by Particle Swarm Optimization—Application to Monitoring Hyporheic Exchange in a Bedrock River

Abstract A self‐potential (SP) data‐inversion algorithm was developed and tested on an analytical model of electrical‐potential profile data attributed to single and multiple polarized electrical sources. The developed algorithm was then validated by an application to SP‐monitoring field data measured on the floodplain of East Fork Poplar Creek, Oak Ridge, Tennessee, to image electrical sources in areas conducive to preferential flow into the flood plain from the bedrock‐lined riverbed. The algorithm combined stochastic source‐localization by particle‐swarm‐optimization (PSO) of electrical sources characterized by simplified geometries with source tomography by regularized weighted least‐squares minimization of a quadratic objective function. Prior information was incorporated by preconditioning the tomography algorithm by PSO results. Variable percentages of random noise were added to analytical‐model data to evaluate the algorithm performance. Results indicated that true parameters of single‐source models were inverted and approximated with small residual error, whereas inversion of analytical‐model data representing multiple electrical sources accurately approximated the locations of the sources but miscalculated some parameters because of the non‐uniqueness of the inverse‐model solution. Source tomography applied to analytical model data during testing produced a spatially continuous parameter field that identified the locations of point‐scale synthetic dipole sources of electrical current flow with varying degrees of accuracy depending on the prior information incorporated into the tomography. When applied to SP‐monitoring field data, the algorithm imaged electrical sources within a known fault that intersects the bedrock riverbed and flood plain of East Fork Poplar Creek and depicted dynamic electrical conditions attributed to hyporheic exchange.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

Evolution of the Large Scale Circulation, Cloud Structure and Regional Water Cycle Associated with the South China Sea Monsoon During May-June, 1998

In this paper, changes in the large-scale circulation, cloud structures and regional water cycle associated with the evolution of the South China Sea (SCS) monsoon in May-June 1998 were investigated using data from the Tropical Rainfall Measuring Mission (TRMM) and field data from the South China Sea Monsoon Experiment (SCSMEX). Results showed that both tropical and extratropical processes strongly influenced the onset and evolution of the SCS monsoon. Prior to the onset of the SCS monsoon, enhanced convective activities associated with the Madden and Julian Oscillation were detected over the Indian Ocean, and the SCS was under the influence of the West Pacific Anticyclone (WPA) with prevailing low level easterlies and suppressed convection. Establishment of low-level westerlies across Indo-China, following the development of a Bay of Bengal depression played an important role in building up convective available potential energy over the SCS. The onset of SCS monsoon appeared to be triggered by the equatorward penetration of extratropical frontal system, which was established over the coastal region of southern China and Taiwan in early May. Convective activities over the SCS were found to vary inversely with those over the Yangtze River Valley (YRV). Analysis of TRMM microwave and precipitation radar data revealed that during the onset phase, convection over the northern SCS consisted of squall-type rain cell embedded in meso-scale complexes similar to extratropical systems. The radar Z-factor intensity indicated that SCS clouds possessed a bimodal distribution, with a pronounced signal (less than 30dBz) at a height of 2-3 km, and another one (less than 25 dBz) at the 8-10 km level, separated by a well-defined melting level indicated by a bright band at around 5-km level. The stratiform-to-convective cloud ratio was approximately 1:1 in the pre-onset phase, but increased to 5:1 in the active phase. Regional water budget calculations indicated that during the active phase, the SCS was a strong sink (E-P much less than 0) of atmospheric moisture, with the primary source of moisture coming from regions further west over Indo-China and the eastern Indian Ocean. Before onset and during the break, the SCS was a moisture source (E-P greater than ) to the overlying atmosphere. In particular, the SCS provided the bulk of moisture to the torrential rain over the YRV in mid-June 1998.

Lau, William K.-M.↗
DOE OSTI2023

Arctic Permafrost Thawing Enhances Sulfide Oxidation

Permafrost degradation is altering biogeochemical processes throughout the Arctic. Thaw-induced changes in organic matter transformations and mineral weathering reactions are impacting fluxes of inorganic carbon (IC) and alkalinity (ALK) in Arctic rivers. However, the net impact of these changing fluxes on the concentration of carbon dioxide in the atmosphere (pCO 2 ) is relatively unconstrained. Resolving this uncertainty is important as thaw-driven changes in the fluxes of IC and ALK could produce feedbacks in the global carbon cycle. Enhanced production of sulfuric acid through sulfide oxidation is particularly poorly quantified despite its potential to remove ALK from the ocean-atmosphere system and increase pCO 2 , producing a positive feedback leading to more warming and permafrost degradation. In this work, we quantified weathering in the Koyukuk River, a major tributary of the Yukon River draining discontinuous permafrost in central Alaska, based on water and sediment samples collected near the village of Huslia in summer 2018. Using measurements of major ion abundances and sulfate (SO 4 2- ) sulfur ( 34 S/ 32 S) and oxygen ( 18 O/ 16 O) isotope ratios, we employed the MEANDIR inversion model to quantify the relative importance of a suite of weathering processes and their net impact on pCO 2 . Calculations found that approximately 80% of SO 4 2- in mainstem samples derived from sulfide oxidation with the remainder from evaporite dissolution. Moreover, 34 S/ 32 S ratios, 13 C/ 12 C ratios of dissolved IC, and sulfur X-ray absorption spectra of mainstem, secondary channel, and floodplain pore fluid and sediment samples revealed modest degrees of microbial sulfate reduction within the floodplain. Weathering fluxes of ALK and IC result in lower values of pCO 2 over timescales shorter than carbonate compensation (~10 4 yr) and, for mainstem samples, higher values of pCO 2 over timescales longer than carbonate compensation but shorter than the residence time of marine (~10 7 yr). Furthermore, the absolute concentrations of and Mg 2+ in the Koyukuk River, as well as the ratios of SO 4 2- and Mg 2+ to other dissolved weathering products, have increased over the past 50 years. Through analogy to similar trends in the Yukon River, we interpret these changes as reflecting enhanced sulfide oxidation due to ongoing exposure of previously frozen sediment and changes in the contributions of shallow and deep flow paths to the active channel. Overall, these findings confirm that sulfide oxidation is a substantial outcome of permafrost degradation and that the sulfur cycle responds to permafrost thaw with a timescale-dependent feedback on warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Variational encoder geostatistical analysis (VEGAS) with an application to large scale riverine bathymetry

Estimation of riverbed profiles, also known as bathymetry, plays a vital role in many applications, such as safe and efficient inland navigation, prediction of bank erosion, land subsidence, and flood risk management. The high cost and complex logistics of direct bathymetry surveys, i.e, depth imaging, have encouraged the use of indirect measurements such as surface flow velocities. However, estimating high-resolution bathymetry from indirect measurements is an inverse problem that can be computationally challenging. Here, we propose a reduced-order model (ROM) based approach that utilizes a variational autoencoder (VAE), a type of deep neural network with a narrow layer in the middle, to compress bathymetry and flow velocity information and accelerate bathymetry inverse problems from flow velocity measurements. In our application, the shallow-water equations (SWE) with appropriate boundary conditions (BCs), e.g., the discharge and/or the free surface elevation, constitute the forward problem, to predict flow velocity. Then, ROMs of the SWEs are constructed on a nonlinear manifold of low dimensionality through a variational encoder and the bathymetry inversion problem is derived on the low-dimensional latent space in a Hierarchical Bayesian setting. Further, the reformulation allows variational inference with a small number (e.g., $\mathscr{O}$ (100) of ROM runs and efficient uncertainty quantification. We have tested our inversion approach on a one-mile reach of the Savannah River, GA, USA. Once the neural network is trained (offline stage), the proposed technique can perform the inversion operation orders of magnitude faster than traditional inversion methods that are commonly based on linear projections, such as principal component analysis (PCA), or the principal component geostatistical approach (PCGA). Furthermore, tests show that the algorithm can estimate the bathymetry with good accuracy even with sparse flow velocity measurements.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1973

Power law time dependence of river flood decay and its relationship to long term discharge frequency distribution

Investigations have continued into the possibility that significant information on stream flow rates can be obtained from aerial and satellite imagery of river meander patterns by seeking a correlation between the meander and discharge spectra of rivers. Such a correlation could provide the basis for a simple and inexpensive technique for remote sensing of the water resources of large geographical areas, eliminating the need for much hydrologic recording. The investigation of the nature of the meander and discharge spectra and their interrelationship can also contribute to a more fundamental understanding of the processes of both river meander formation and drainage of large basins. It has been found that floods decay with an inverse power law dependence on time. The exponent of this dependence varies from river to river and even from station to station along the same river. This power law time dependence makes possible the forecasting of river discharge with an uncertainty of about 5% for as long as a month following the flood peak.

Schubert, G.↗
DOE OSTI2021

Optical Classification of Lower Amazon Waters Based on In Situ Data and Sentinel-3 Ocean and Land Color Instrument Imagery

Optical water types (OWTs) were identified from an in situ dataset of concomitant biogeochemical and optical parameters acquired in the Amazon River and its tributaries, in the Lower Amazon region, at different hydrological conditions from 2014 to 2017. A seasonal bio-optical characterization was performed. The k-means classification was applied to the in situ normalized reflectance spectra (rn(λ)), allowing the identification of four OWTs. An optical index method was also applied to the rn(λ) defining the thresholds of the OWTs. Next, level-3 Sentinel-3 Ocean and Land Color Instrument images representative of the seasonal discharge conditions were classified using the identified in situ OWTs as reference. The differences between Amazon River and clearwater tributary OWTs were dependent on the hydrological dynamics of the Amazon River, also showing a strong seasonal variability. Each OWT was associated with a specific bio-optical and biogeochemical environment assessed from the corresponding absorption coefficient values of colored dissolved organic matter (aCDOM) and particulate matter (ap), chlorophyll-a and suspended particulate matter (SPM) concentrations, and aCDOM/ap ratio. The rising water season presented a unique OWT with high SPM concentration and high relative contribution of ap to total absorption compared to the other OWTs. This bio-optical characterization of Lower Amazon River waters represents a first step for developing remote sensing inversion models adjusted to the optical complexity of this region.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1973

Remote sensing of stream flow rates - Correlation of meander and discharge spectra

Results of a study of river meander patterns and discharges, in which attempt was made to correlate the discharge spectrum of a river with the river meander power spectrum determined from aerial and satellite imagery. Some significant characteristics of both the discharge and the meander spectra have been discovered. Discharge frequency spectra based on long-term records of daily streamflow are found to have an inverse power-law dependence on discharge. This is shown to reflect the short-term decay of individual floods which are found to have an inverse power-law dependence on time. Meander power spectra for a number of river reaches, digitized from aerial photography, also show significant structure, the power spectral density having an inverse power-law dependence on wave number over one or more portions of the spectrum with breaks in the spectra at characteristic wave numbers. A number of examples of typical discharge and meander spectra are shown.

Lingenfelter, R. E.↗
NASA NTRS2002

Evolution of Topography in Glaciated Mountain Ranges

This thesis examines the response of alpine landscapes to the onset of glaciation. The basic approach is to compare fluvial and glacial laudscapes, since it is the change from the former to the latter that accompanies climatic cooling. This allows a detailed evaluation of hypotheses relating climate change to tectonic processes in glaciated mountain belts. Fieldwork was carried out in the eastern Sierra Nevada, California, and the Sangre de Cristo Range, Colorado, alongside digital elevation model analyses in the western US, the Southern Alps of New Zealand, and the Himalaya of northwestern Pakistan. hypothesis is overstated in its appeal to glacial erosion as a major source of relief production and subsequent peak uplift. Glaciers in the eastern Sierra Nevada and the western Sangre de Cristos have redistributed relief, but have produced only modest relief by enlarging drainage basins at the expense of low-relief topography. Glaciers have lowered valley floors and ridgelines by similar amounts, limiting the amount of "missing mass' that can be generated, and causing a decrease in drainage basin relief. The principal response of glaciated landscapes to rapid rock uplift is the development of towering cirque headwalls. This represents considerable relief production, but is not caused by glacial erosion alone. Large valley glaciers can maintain their low gradient regardless of uplift rate, which supports the "glacial buzzsaw" hypothesis. However, the inability of glaciers to erode steep hillslopes as rapidly can cause mean elevations to rise. Cosmogenic isotope dating is used to show that (i) where plucking is active, the last major glaciation removed sufficient material to reset the cosmogenic clock; and (ii) former glacial valley floors now stranded near the crest of the Sierra Nevada are at varying stages of abandonment, suggesting a cycle of drainage reorganiszation and relief inversion due to glacial erosion similar to that observed in river networks. Glaciated landscapes are quite distinct from their fluvial counterparts in both landforms and processes. Given the scarcity of purely fluvial, active mountain ranges, it is essential that glacial erosion be considered amongst the processes sculpting active orogenic belts.

Brocklehurst, Simon H.↗
DOE OSTI2020

Using onset times from frequent seismic surveys to understand fluid flow at the Peace River Field, Canada

With our limited knowledge of the relationship between changes in the state of an aquifer or reservoir and the corresponding changes in the elastic moduli, that is the rock physics model, hampers the effective use of time-lapse seismic observations for estimating flow properties within the Earth. A central problem is the complicated dependence of the magnitude of time-lapse changes on the saturation, pressure, and temperature changes within an aquifer or reservoir. We describe an inversion methodology for reservoir characterization that uses onset times, the calendar time of the change in seismic attributes, rather than the magnitude of the changes. We find that onset times are much less sensitive than magnitudes to the rock physics model used to relate time-lapse observations to changes in saturation, temperature and fluid pressure. We apply the inversion scheme to observations from daily monitoring of enhanced oil recovery at the Peace River field in Canada. An array of 1492 buried hydrophones record seismic signals from 49 buried sources. Time-shifts for elastic waves traversing the reservoir are extracted from the daily time-lapse cubes. In our analysis 175 images of time-shifts are transformed into a single map of onset times, leading to a substantial reduction in the volume of data. These observations are used in conjunction with bottom hole pressure data to infer the initial conditions prior to the injection, and to update the reservoir permeability model. The combination of a global and local inversion scheme produces a collection of reservoir models that are best described by three clusters. The updated model leads to a nearly 70 percent reduction in seismic data misfit. The final set of solutions successfully predict the observed normalized pressure history during the soak and flow-back into the wells between 82 and 175 days into the cyclic steaming operation.

58 GEOSCIENCES↗
DOE OSTI2026

3D seismic imaging of a fracture damage zone controlling reservoir compartmentalization at the Raft River EGS using multi-azimuth walkaway VSP

Accurate imaging of steeply dipping fracture zones in crystalline enhanced geothermal systems (EGS) is critical for constraining permeability architecture and guiding stimulation design. However, such structures remain poorly resolved by conventional surface seismic methods. We present a fully three-dimensional (3D) elastic-waveform inversion-migration workflow applied to multi-azimuth walkaway vertical seismic profiling (VSP) data acquired in a deviated borehole at the Raft River EGS. The workflow integrates first-arrival traveltime tomography, multi-scale elastic waveform inversion (EWI), and elastic least-squares reverse-time migration (ELSRTM) to recover high-resolution compressional-and shear-wave velocity models and to image structural discontinuities in the crystalline basement. The results reveal a laterally continuous low-velocity anomaly, with shear-wave velocity reductions of 25-30%, consistent with fractureinduced mechanical weakening. Two steeply dipping discontinuities bound a 50-80 m wide fracture damage zone. Independent constraints from microseismic clustering and geochemical compartmentalization corroborate the geometry and structural significance of this feature. Synthetic modeling further confirms that structures of this scale are resolvable (~30 m). These findings indicate that the Narrows structure is a distributed fracture damage zone rather than a discrete fault plane. By resolving fracture-zone geometry at the tens-of-meters scale using a single borehole, this workflow provides a practical and transferable approach for improving structural characterization, reducing uncertainty in permeability architecture, and supporting reservoir modeling and stimulation design in fractured crystalline EGS reservoirs.

58 GEOSCIENCES↗
DOE OSTI2026

Determining the Solubility Behavior of Kogarkoite in Simulated Nuclear Waste

Kogarkoite (Na 3 FSO 4 ) is a sparingly soluble fluoride–sulfate double salt that has been identified in high level nuclear waste sludge at the Hanford Site and, more recently, in sludge batch compilation samples at the Savannah River Site (SRS). Due to its complex dissolution behavior, which exhibits an inverse dependence on sodium ion activity, the presence of this mineral poses significant challenges to waste retrieval and processing. Incomplete dissolution during sludge washing can lead to the retention of fluoride and sulfate in the high-level waste feed, potentially causing the formation of corrosive, immiscible molten salt layers, known as "glass gall,” in vitrification melters. Current efforts to optimize flowsheet parameters and wash-water volumes are hindered by the absence of a commercially available, certified reference material, which prevents the accurate calibration of analytical methods and the verification of dissolution kinetics. To address this critical gap, this research focuses on the laboratory synthesis of pure Kogarkoite to serve as a standard for comprehensive solubility and washing performance testing. A coupled synthesis and simulant campaign was executed using an evaporative crystallization protocol designed to replicate the dynamic concentration effects observed in tank farm operations. Thirteen simulant matrices were prepared by dissolving systematically varied ratios of sodium fluoride (NaF) and sodium sulfate (Na 2 SO 4 ) in deionized water under three distinct caustic regimes: 0.0 g (control), 4.0 g (~1 M), and 12.0 g (~3 M) sodium hydroxide (NaOH). While thermodynamic equilibrium models suggest that high-caustic environments should favor the stability of the double salt7, results from this evaporative study at 25 0 C revealed a distinct kinetic divergence. Simulants with high hydroxide loading predominantly yielded large, blocky crystals of sodium sulfate decahydrate (Na 2 SO 4 .10H 2 O). Successful synthesis of pure Kogarkoite was achieved exclusively in specific NaOH-free compositional windows, where the precipitate manifested as fine, opaque granular aggregates. Ion chromatography (IC) analysis confirmed phase purity through the simultaneous stoichiometric depletion of both fluoride and sulfate from the supernatant. This successful synthesis establishes a reproducible route to generate bulk Kogarkoite, enabling the subsequent phase of quantitative dissolution testing using inhibited water to optimize sludge-batch assembly.

Sarker, Md Sharif [Florida International Univ. (FI↗
DOE OSTI2023

Inverse aqueous transport modeling for emergency response

ALGE is a three-dimensional, finite-difference aqueous transport model that simulates pollutant fate and transport in lakes, rivers, bays, and estuaries by solving the prognostic equations of mass, momentum, and energy. Its current modeling capabilities include transport of dissolved tracer for a series of predefined basins across the continental United States. Recently, an inverse method (also known as backtracking) has been added to ALGE to provide a possible source of a pollutant should one be detected by a sensor in a body of water and a source is not known. This inverse method is a three step process that uses an algorithm to inverse the flow. We demonstrate the new model’s capabilities through simulating the 2021 Piney Point spill in Tampa Bay, Florida (USA). This involves moving tracer backwards from its detection points, encompassing a potential source area, and applying Bayes’ Theorem and $\frac{𝜒}{𝑄}$ to reduce the area within which the true source could be located.

hydrological modeling↗
DOE OSTI2023

Soil architecture and subsurface stratigraphic heterogeneities across land-lake interfaces along Lake Erie.

Uploaded raw datasets from electrical resistivity tomography (ERT), ground penetrating radar (GPR), electromagnetic induction (EMI) and insitu sensors were acquired across land-lake interfaces along Lake Erie. These datasets were useful to delineate the soil architecture and subsurface stratigraphic heterogeneities. We combined multiple geophysical methods with lithostratigraphic logs from borehole cores and data from soil and borehole sensors to investigate spatial heterogeneities across land-lake interfaces. The apparent electrical conductivity (Eca) maps matched closely with the soil maps from public database and also detected additional soil units that were missing in the traditional soil maps. This implies that electromagnetic induction (EMI) could be relied upon for non invasive characterization of soils in sampling restricted sites where only non-invasive measurements are feasible. Results from electrical resistivity tomography (ERT) and Ground penetrating radar (GPR) are consistent with the surficial geology of the study area, and revealed that the vertical variabilities at the sites are mainly driven by stratigraphy.The ERT data (.stg) were acquired from three different sites portage river (PTR), crane creek (CRC) and old woman creek (OWC). The ".stg" files can be opened with ERT inversion software such as EarthImager, ResIpy or Res2DInv. The GPR data (.GPZ) were also acquired from the three sites, and the ".GPZ" data files can be opened with EKKO_project 5 software from Sensors & Software Inc., or with any other standard GPR software. The EM data (.csv) were acquired from the three sites while the soil sensor data (.csv) here were acquired only from crane creek site. The ".csv) files can be opened with MS Excel or any other spread sheet software.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Robust Carbon Dioxide Plume Imaging Using Joint Tomographic Inversion of Seismic Onset Time and Distributed Pressure and Temperature Measurements (Final Report)

We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.

02 PETROLEUM↗