Search NASA⌕ Search

SEARCH · Search NASA

Results for “subsurface flow”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Functional Response of a Near-Surface Soil Microbial Community to a Simulated Underground CO 2 Storage Leak

Understanding the impacts of leaks from geologic carbon sequestration, also known as carbon capture and storage, is key to developing effective strategies for carbon dioxide (CO 2 ) emissions management and mitigation of potential negative effects. Here, we provide the first report on the potential effects of leaks from carbon capture and storage sites on microbial functional groups in surface and near-surface soils. Using a simulated subsurface CO 2 storage leak scenario, we demonstrate how CO 2 flow upward through the soil column altered both the abundance (DNA) and activity (mRNA) of microbial functional groups mediating carbon and nitrogen transformations. These microbial responses were found to be seasonally dependent and correlated to shifts in atmospheric conditions. While both DNA and mRNA levels were affected by elevated CO 2 , they did not react equally, suggesting two separate mechanisms for soil microbial community response to high CO 2 levels. The results did not always agree with previous studies on elevated atmospheric (rather than subsurface) CO 2 using FACE (Free-Air CO 2 Enrichment) systems, suggesting that microbial community response to CO 2 seepage from the subsurface might differ from its response to atmospheric CO 2 increases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A multiphysics coupling framework for exascale simulation of fracture evolution in subsurface energy applications

Predicting the evolution of fractured media is challenging due to coupled thermal, hydrological, chemical and mechanical processes that occur over a broad range of spatial scales, from the microscopic pore scale to field scale. We present a software framework and scientific workflow that couples the pore scale flow and reactive transport simulator Chombo-Crunch with the field scale geomechanics solver in GEOS to simulate fracture evolution in subsurface fluid-rock systems. This new multiphysics coupling capability comprises several novel features. An HDF5 data schema for coupling fracture positions between the two codes is employed and leverages the coarse resolution of the GEOS mechanics solver which limits the size of data coupled, and is, thus, not taxed by data resulting from the high resolution pore scale Chombo-Crunch solver. The coupling framework requires tracking of both before and after coarse nodal positions in GEOS as well as the resolved embedded boundary in Chombo-Crunch. We accomplished this by developing an approach to geometry generation that tracks the fracture interface between the two different methodologies. The GEOS quadrilateral mesh is converted to triangles which are organized into bins and an accessible tree structure; the nodes are then mapped to the Chombo representation using a continuous signed distance function that determines locations inside, on and outside of the fracture boundary. The GEOS positions are retained in memory on the Chombo-Crunch side of the coupling. The time stepping cadence for coupled multiphysics processes of flow, transport, reactions and mechanics is stable and demonstrates temporal reach to experimental time scales. The approach is validated by demonstration of 9 days of simulated time of a core flood experiment with fracture aperture evolution due to invasion of carbonated brine in wellbore-cement and sandstone. We also demonstrate usage of exascale computing resources by simulating a high resolution version of the validation problem on OLCF Frontier.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Efficacy of acetate-amended biostimulation for uranium sequestration: Combined analysis of sediment/groundwater geochemistry and bacterial community structure

Systematic flow-through column experiments were conducted using sediments and ground water collected from different subsurface localities at the U.S. Department of Energy's Integrated Field Research Challenge site in Rifle, Colorado. The principal purpose of this study is to gain a better understanding of the interactive effects of groundwater geochemistry, sediment mineralogy, and indigenous bacterial community structures on the efficacy of uranium removal from the groundwater with/without acetate amendment. Overall, we find that the subtle variations in the sediments' mineralogy, particle size, redox conditions, as well as contents of metal(loid) co-contaminants showed a pronounced effect on the associated bacterial population and composition, which mainly determines the system's performance with respect to uranium removal. Positive relationship was identified between the abundance of dissimilatory sulfate-reduction genes (i.e., drsA), markers of sulfatereducing bacteria, and the sediments' propensity to sequester aqueous uranium. In contrast, no obvious connections were observed between the abundance of common iron-reducing bacteria, e.g., Geobacter spp., and the sediments' ability to sequester uranium. In the sediments with low bacterial biomass and the absence of sulfate-reducing conditions, abiotic adsorption onto mineral surfaces such as phyllosilicates likely played a relatively major role in the attenuation of aqueous uranium; however, in these scenarios, acetate amendment induced detectable rebounds in the effluent uranium concentrations. The results of this study suggest that reductive immobilization of uranium can be achieved under predominantly sulfate-reducing conditions, and provide insight into the integrated roles of various biogeochemical components in long-term uranium sequestration.

58 GEOSCIENCES↗
DOE OSTI2014

Monitoring CO 2 sequestration into deep saline aquifer and associated salt intrusion using coupled multiphase flow modeling and time lapse electrical resistivity tomography

Successful geological storage and sequestration of carbon dioxide (CO2) require efficient monitoring of the migration of CO2 plume during and after large-scale injection in order to verify the containment of the injected CO2 within the target formation and to evaluate potential leakage risk. Field studies have shown that surface and cross-borehole electrical resistivity tomography (ERT) can be a useful tool in imaging and characterizing solute transport in heterogeneous subsurface. In this synthetic study, we have coupled a 3-D multiphase flow model with a parallel 3-D time-lapse ERT inversion code to explore the feasibility of using time-lapse ERT for simultaneously monitoring the migration of CO2 plume in deep saline formation and potential brine intrusion into shallow fresh water aquifer. Direct comparisons of the inverted CO2 plumes resulting from ERT with multiphase flow simulation results indicate the ERT could be used to delineate the migration of CO2 plume. Detailed comparisons on the locations, sizes and shapes of CO2 plume and intruded brine plumes suggest that ERT inversion tends to underestimate the area review of the CO2 plume, but overestimate the thickness and total volume of the CO2 plume. The total volume of intruded brine plumes is overestimated as well. However, all discrepancies remain within reasonable ranges. Our study suggests that time-lapse ERT is a useful monitoring tool in characterizing the movement of injected CO2 into deep saline aquifer and detecting potential brine intrusion under large-scale field injection conditions.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2019

Reactive Transport Simulation of Fracture Channelization and Transmissivity Evolution

Underground fractures serve as flow conduits, and they may produce unwanted migration of water and other fluids in the subsurface. An example is the migration and leakage of greenhouse gases in the context of geologic carbon sequestration. This study has generated new understanding about how acids erode carbonate fracture surfaces and the positive feedback between reaction and flow. A two-dimensional reactive transport model was developed and used to investigate the extent to which geochemical factors influence fracture permeability and transmissivity evolution in carbonate rocks. The only mineral modeled as reactive is calcite, a fast-reacting mineral that is abundant in subsurface formations. The X-ray computed tomography dataset from a previous experimental study of fractured cores exposed to carbonic acid served as a testbed to benchmark the model simulation results. The model was able to capture not only erosion of fracture surfaces but also the specific phenomenon of channelization, which produces accelerating transmissivity increase. Results corroborated experimental findings that higher reactivity of the influent solution leads to strong channelization without substantial mineral dissolution. Simulations using mineral maps of calcite in a specimen of Amherstburg limestone demonstrated that mineral heterogeneity can either facilitate or suppress the development of flow channels depending on the spatial patterns of reactive mineral. In these cases, fracture transmissivity may increase rapidly, increase slowly, or stay constant, and for all these possibilities, the calcite mineral continues to dissolve. Collectively, these results illustrate that fluid chemistry and mineral spatial patterns need to be considered in predictions of reaction-induced fracture alteration and risks of fluid migration.

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

The EGS Collab Project: Stimulation and Simulation

Development of Enhanced Geothermal Systems (EGS) requires accurate predictions of flow rates and temperatures at production wells over time. The concept of EGS is simple; however, complications from heterogeneity and complexity of fracture pathways may lead to channeling, short-circuiting, and premature thermal breakthrough. In the EGS Collab project, we are establishing a suite of highly monitored and well characterized intermediate-scale (~10-20 m) field test beds along with fracture stimulation and interwell flow tests to better understand processes that control formation of effective subsurface heat exchangers. EGS Collab experiments will provide a means of testing tools and concepts that could later be employed under geothermal reservoir conditions at FORGE or elsewhere. Key to the project is using numerical simulations in the experiment design and interpretation of results. Our tests will be well-controlled, in situ experiments focused on rock fracture behavior and permeability enhancement. Our pre- and post-test modeling results of each field experiment compared with detailed measurements will allow for improved model prediction and validation. We will use comprehensive instrumentation to collect high-quality and high-resolution geophysical and other fracture characterization and fluid flow data. We will analyze these data and compare them with models and field observations to further elucidate the basic relationships between stress, induced seismicity, and permeability enhancement. We will also observe and quantify other key governing parameters that impact permeability, and will attempt to understand how these parameters might change throughout the development and operation of an EGS project with the ultimate goal of enabling commercial viability of EGS. Our first set of experiments will be performed at the Sanford Underground Research Facility (SURF) in South Dakota. Our team is well underway with building the first field experiment test bed planned for this project, which is supported by the US Department of Energy’s Geothermal Technologies Office.

EGS Collab↗
DOE OSTI2018

Feasibility of a Deep Direct-Use Geothermal System at the University of Illinois Urbana-Champaign

This study assesses the feasibility of using deep direct-use (DDU) geothermal energy in agricultural research facilities on the University of Illinois at Urbana-Champaign campus to exploit low-temperature sedimentary basins, such as the Illinois Basin. Subsurface components of the system include extraction and injection wells and downhole pumps. Surface equipment includes heat pumps/exchangers, and fluid transport and monitoring systems. Two geologic formations in the region exhibit a potential as sources for geothermal energy, based on pre initial temperatures and flow rates of fluids. The St. Peter and Mt. Simon Sandstones lie at depths of 634 and 1,280 m, respectively. Geocellular modeling is used to characterize the reservoirs. A St. Peter Sandstone model was made for an area south of the campus. Petrophysical and geothermal properties used are based on data from the closest wells penetrating the formations. Characterization of the Mt. Simon Sandstone is in progress and is not discussed here. Extraction and injection flows simulated with different wellbore configurations provide estimates of fluid flow out of and into the reservoir. The models are used to optimize flow rates, bottomhole pressure, and temperature of the produced fluid. Individual wellbore models simulate subsurface heat loss and gain, providing guidance on the optimal type and amount of insulation in the wellbore. Design of the surface facilities will address aspects of fluid delivery, heat exchange, capital operating costs, heat loss, and corrosion. Heat capacity and flow rates are assessed to estimate life-cycle costs and benefits, including the environmental benefits of reducing greenhouse gases and water use and increased energy efficiency. A preliminary analysis of surface configurations for the DDU system (including cascading applications) based on building heat loads is being conducted to identify multiple system designs that will maximize performance, energy efficiency, and cost recovery.

15 GEOTHERMAL ENERGY↗
DOE OSTI2012

Mineral Dissolution and Secondary Precipitation on Quartz Sand in Simulated Hanford Tank Solutions Affecting Subsurface Porosity

Highly alkaline nuclear waste solutions have been released from underground nuclear waste storage tanks and pipelines into the vadose zone at the U.S. Department of Energy’s Hanford Site in Washington, causing mineral dissolution and re-precipitation upon contact with subsurface sediments. High pH caustic NaNO3 solutions with and without dissolved Al were reacted with quartz sand through flow-through columns stepwise at 45, 51, and 89°C to simulate possible reactions between leaked nuclear waste solution and primary subsurface mineral. Upon reaction, Si was released from the dissolution of quartz sand, and nitrate-cancrinite [Na8Si6Al6O24(NO3)2] precipitated on the quartz surface as a secondary mineral phase. Both steady-state dissolution and precipitation kinetics were quantified, and quartz dissolution apparent activation energy was determined. Mineral alteration through dissolution and precipitation processes results in pore volume and structure changes in the subsurface porous media. In this study, the column porosity increased up to 40.3% in the pure dissolution column when no dissolved Al was present in the leachate, whereas up to a 26.5% porosity decrease was found in columns where both dissolution and precipitation were observed because of the presence of Al in the input solution. The porosity change was also confirmed by calculation using the dissolution and precipitation rates and mineral volume changes.

alkaline nuclear waste solution↗
DOE OSTI2025

Factors controlling injection-induced rupture of intersecting faults during geological sequestration of CO 2

This study addresses coupled multiphase fluid flow and geomechanics effects on potential fault activation associated with subsurface CO 2 injection around intersecting faults. An enhanced fault-representation model is used to capture geomechanical responses of two intersecting faults with finite length during CO 2 injection. The faults are embedded in a strike-slip stress regime of a caprock-reservoir-basement system with the faults represented by zero-thickness interfaces with adjacent finite-thickness damage zones. A sensitivity analysis is conducted to study the effect of fault permeability, slip-weakening behavior, well location relative to the orientation of faults, and well placement (the number and location of injection wells). Five metrics (pressure, CO 2 plume, shear state on the fault, as well as shear displacement and stress path at selected fault monitoring points) are selected to assess CO 2 migration and reactivation of intersecting faults. The results show that induced ruptures are favored by low permeability faults due to high pressure buildup and by slip-weakening behavior resulting from fault strength reduction. The location of one injection well relative to fault orientation determines the magnitude of changes in effective normal stress and shear stress, affecting the location of induced ruptures. Well placement (two injection wells used in the paper) dominates pressure diffusion around the intersection and tips of faults. This redistributes changes in effective normal stress caused by each injection well, influencing the spatial distribution of ruptures along faults. A larger injection volume induces far-field ruptures that are controlled by stress transfer within the injection layer. The findings presented here can provide valuable insights into engineering operations for a long-term, safe, and reliable geologic CO 2 storage.

Fault permeability↗
DOE OSTI2025

Joint state-parameter estimation for the reduced fracture model via the united filter

Here, in this paper, we introduce an effective United Filter method for jointly estimating the solution state and physical parameters in flow and transport problems within fractured porous media. Fluid flow and transport in fractured porous media are critical in subsurface hydrology, geophysics, and reservoir geomechanics. Reduced fracture models, which represent fractures as lower-dimensional interfaces, enable efficient multi-scale simulations. However, reduced fracture models also face accuracy challenges due to modeling errors and uncertainties in physical parameters such as permeability and fracture geometry. To address these challenges, we propose a United Filter method, which integrates the Ensemble Score Filter (EnSF) for state estimation with the Direct Filter for parameter estimation. EnSF, based on a score-based diffusion model framework, produces ensemble representations of the state distribution without deep learning. Meanwhile, the Direct Filter, a recursive Bayesian inference method, estimates parameters directly from state observations. The United Filter combines these methods iteratively: EnSF estimates are used to refine parameter values, which are then fed back to improve state estimation. Numerical experiments demonstrate that the United Filter method surpasses the state-of-the-art Augmented Ensemble Kalman Filter, delivering more accurate state and parameter estimation for reduced fracture models. This framework also provides a robust and efficient solution for PDE-constrained inverse problems with uncertainties and sparse observations.

Bayesian inference↗
DOE OSTI2016

Calibration of an agricultural-hydrological model (RZWQM2) using surrogate global optimization

Robust calibration of an agricultural-hydrological model is critical for simulating crop yield and water quality and making reasonable agricultural management. However, calibration of the agricultural-hydrological system models is challenging because of model complexity, the existence of strong parameter correlation, and significant computational requirements. Therefore, only a limited number of simulations can be allowed in any attempt to find a near-optimal solution within an affordable time, which greatly restricts the successful application of the model. The goal of this study is to locate the optimal solution of the Root Zone Water Quality Model (RZWQM2) given a limited simulation time, so as to improve the model simulation and help make rational and effective agricultural-hydrological decisions. To this end, we propose a computationally efficient global optimization procedure using sparse-grid based surrogates. We first used advanced sparse grid (SG) interpolation to construct a surrogate system of the actual RZWQM2, and then we calibrate the surrogate model using the global optimization algorithm, Quantum-behaved Particle Swarm Optimization (QPSO). As the surrogate model is a polynomial with fast evaluation, it can be efficiently evaluated with a sufficiently large number of times during the optimization, which facilitates the global search. We calibrate seven model parameters against five years of yield, drain flow, and NO 3 -N loss data from a subsurface-drained corn-soybean field in Iowa. Results indicate that an accurate surrogate model can be created for the RZWQM2 with a relatively small number of SG points (i.e., RZWQM2 runs). Compared to the conventional QPSO algorithm, our surrogate-based optimization method can achieve a smaller objective function value and better calibration performance using a fewer number of expensive RZWQM2 executions, which greatly improves computational efficiency.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Permeability Evolution and Frictional Stability of Fabricated Fractures With Specified Roughness

Roughness is widely observed on natural fractures, and its impact on the potential for induced seismicity and associated fluid migration in the subsurface remains unclear. Here we perform fracture shearing and fluid flow experiments on artificially fabricated fractures with specified roughness to investigate the role of fracture roughness on frictional properties and permeability evolution. Given the experimental conditions, we observe that rough fractures show high roughness ratio S q /L w and return higher frictional strength due to the presence of cohesive interlocking asperities. Rough fracture surfaces show velocity strengthening behavior in the initial shearing stage, which may evolve to velocity neutral and velocity weakening at greater displacements—suggesting a dynamic weakening that rough fractures become less stable with shearing. The surface roughness exerts a dominant control on permeability evolution over the entire shearing history. Permeability declines monotonically by about 2 orders of magnitude for smooth fractures. For high roughness fractures, the permeabilities evolve episodically due to cycled compaction and dilation during shearing. With a slip distance of 6 to 8 mm, permeability of the rough surface may enhance up to an order of magnitude, but significant permeability reduction may also occur for rough samples when asperities are highly worn with gouge clogging flow paths. However, there is no obvious correlation between permeability evolution and frictional behavior for rough fracture samples when fractures are subject to sudden sliding velocity changes.

58 GEOSCIENCES↗
DOE OSTI2015

Self-organized iron-oxide cementation geometry as an indicator of paleo-flows

Widespread iron oxide precipitation from groundwater in fine-grained red beds displays various patterns, including nodulation, banding and scallops and fingers. Hematite nodules have been reported also from the Meridiani Planum site on Mars and interpreted as evidence for the ancient presence of water on the red planet. Here we show that such patterns can autonomously emerge from a previously unrecognized Ostwald ripening mechanism and they capture rich information regarding ancient chemical and hydrologic environments. A linear instability analysis of the reaction-transport equations suggests that a pattern transition from nodules to bands may result from a symmetry breaking of mineral dissolution and precipitation triggered by groundwater advection. Round nodules tend to develop under nearly stagnant hydrologic conditions, while repetitive bands form in the presence of persistent water flows. Since water circulation is a prerequisite for a sustainable subsurface life, a Martian site with iron oxide precipitation bands, if one were found, may offer a better chance for detecting extraterrestrial biosignatures on Mars than would sites with nodules.

58 GEOSCIENCES↗
DOE OSTI2026

Modeling supercritical CO 2 flow and mineralization in reactive host rocks with PFLOTRAN v7.0

Understanding the flow and reactivity of CO 2 injected into geological reservoirs is important for many subsurface applications including secure geologic carbon storage (GCS), critical mineral extraction, enhanced geothermal systems (EGS), and enhanced oil recovery (EOR). Traditionally, subsurface CO 2 injection for GCS applications has focused on geologic formations with favorable subsurface configurations for CO 2 migration and trapping through non-reactive mechanisms such as structural, solubility, and petrophysical trapping. Recently, CO 2 -reactive rocks such as mafic and ultramafic basalts have been investigated for their potential to react with injected CO 2 in situ to simultaneously dissolve host rock minerals and mineralize CO 2 as carbonates. Engineering rapid CO 2 mineralization in the subsurface is attractive because of the increased density of stored CO 2 , the additional safety factors associated with solidification, and the potential to extract valuable critical minerals. Here we present recent developments in the parallel flow and reactive transport simulator PFLOTRAN to model coupled CO 2 -brine flow and reactive transport for a wide range of injection and production applications involving reactive CO 2 -brine systems. These developments are based on the well established and trusted CO 2 flow capabilities in the STOMP-CO 2 simulator. New capabilities added to PFLOTRAN include new CO 2 -brine equations of state with optional thermal coupling, several new constitutive relationships like capillary pressure smoothing and scanning path hysteresis, a fully implicit well model, and native linkage with PFLOTRAN's well-established reactive transport libraries. A series of benchmarks between PFLOTRAN and STOMP-CO 2 verify the newly developed CO 2 -brine flow capabilities. Demonstrations of coupled CO 2 -brine flow modeling and reactive transport show how CO 2 mineralization can be engineered in reactive host rocks. Finally, an example use case involving copper leaching by CO 2 and critical mineral extraction is presented to showcase the strengths of this new implementation. Several limitations still remain, including limited availability of field data to parameterize models. Future work should constrain the evolution of mineral surface area during mineralization and the temperature and/or pH dependence of geochemical reactions for specific systems of interest.

Critical Minerals↗
DOE OSTI2015

Radiochemically-supported microbial communities. A potential mechanism for biocolloid production of importance to actinide transport

The work described here revealed the presence of diverse microbial communities located across 19 subsurface sites at the NNSS/NTTR and nearby locations. Overall, the diversity of microorganisms was high for subsurface habitats and variable between sites. As of this writing, preparations are being made to combine the Illumina sequences and 16S rRNA clone libraries with other non-NNSS/NTTR well sites of Southern Nevada Regional Flow System for a publication manuscript describing our very broad landscape scale survey of subsurface microbial diversity. Isolates DRI-13 and DRI-14 remain to be fully characterized and named in accordance with the conventions established by Bergey's Manual of Systematic Bacteriology. In preparation to be published, these microorganisms will be submitted to the American Type Culture Collection (ATCC) and the Deutsche Sammlung von Mikroorganismen und Zellkulturen GmbH (DSMZ).It is anticipated that the data resulting from this study in combination with other data sets that will allow us to produce a number of publications that will be impactful to the subsurface microbiology community.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Predictive Inverse Model for Advective Heat Transfer in a Short–Circuited Fracture: Dimensional Analysis, Machine Learning, and Field Demonstration

Identifying fluid flow maldistribution in planar geometries is a well–established problem in subsurface science/engineering. Of particular importance to the thermal performance of enhanced (or “engineered”) geothermal systems is identifying the existence of nonuniform (i.e., heterogeneous) permeability and subsequently predicting advective heat transfer. Here, machine learning via a genetic algorithm (GA) identifies the spatial distribution of an unknown permeability field in a two–dimensional Hele–Shaw geometry (i.e., parallel plates). The inverse problem is solved by minimizing the L2 norm between simulated residence time distribution (RTD) and measurements of an inert tracer breakthrough curve (BTC) (C–Dot nanoparticle). Principal component analysis (PCA) of spatially correlated permeability fields enabled reduction of the parameter space by more than a factor of 10 and restricted the inverse search to reservoir–scale permeability variations. Thermal experiments and tracer tests conducted at the mesoscale Altona Field Laboratory (AFL) demonstrate that the method accurately predicts the effects of extreme flow channeling on heat transfer in a single bedding–plane rock fracture. However, this is only true when the permeability distributions provide adequate matches to both tracer RTD and frictional pressure loss. Without good agreement to frictional pressure loss, it is still possible to match a simulated RTD to measurements, but subsequent predictions of heat transfer are grossly inaccurate. Here, the results of this study suggest that it is possible to anticipate the thermal effects of flow maldistribution, but only if both simulated RTDs and frictional pressure loss between fluid inlets and outlets are in good agreement with measurements.

42 ENGINEERING↗
DOE OSTI2020

Predicting Water Cycle Characteristics from Percolation Theory and Observational Data

The fate of water and water-soluble toxic wastes in the subsurface is of high importance for many scientific and practical applications. Although solute transport is proportional to water flow rates, theoretical and experimental studies show that heavy-tailed (power-law) solute transport distribution can cause chemical transport retardation, prolonging clean-up time-scales greatly. However, no consensus exists as to the physical basis of such transport laws. In percolation theory, the scaling behavior of such transport rarely relates to specific medium characteristics, but strongly to the dimensionality of the connectivity of the flow paths (for example, two- or three-dimensional, as in fractured-porous media or heterogeneous sediments), as well as to the saturation characteristics (i.e., wetting, drying, and entrapped air). In accordance with the proposed relevance of percolation models of solute transport to environmental clean-up, these predictions also prove relevant to transport-limited chemical weathering and soil formation, where the heavy-tailed distributions slow chemical weathering over time. The predictions of percolation theory have been tested in laboratory and field experiments on reactive solute transport, chemical weathering, and soil formation and found accurate. Recently, this theoretical framework has also been applied to the water partitioning at the Earth’s surface between evapotranspiration, ET, and run-off, Q, known as the water balance. A well-known phenomenological model by Budyko addressed the relationship between the ratio of the actual evapotranspiration (ET) and precipitation, ET/P, versus the aridity index, ET 0 /P, with P being the precipitation and ET 0 being the potential evapotranspiration. Existing work was able to predict the global fractions of P represented by Q and ET through an optimization of plant productivity, in which downward water fluxes affect soil depth, and upward fluxes plant growth. In the present work, based likewise on the concepts of percolation theory, we extend Budyko’s model, and address the partitioning of run-off Q into its surface and subsurface components, as well as the contribution of interception to ET. Using various published data sources on the magnitudes of interception and information regarding the partitioning of Q, we address the variability in ET resulting from these processes. The global success of this prediction demonstrated here provides additional support for the universal applicability of percolation theory for solute transport as well as guidance in predicting the component of subsurface run-off, important for predicting natural flow rates through contaminated aquifers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Imaging of a fluid injection process using geophysical data — A didactic example

In many subsurface industrial applications, fluids are injected into or withdrawn from a geologic formation. It is of practical interest to quantify precisely where, when, and by how much the injected fluid alters the state of the subsurface. Routine geophysical monitoring of such processes attempts to image the way that geophysical properties, such as seismic velocities or electrical conductivity, change through time and space and to then make qualitative inferences as to where the injected fluid has migrated. The more rigorous formulation of the time-lapse geophysical inverse problem forecasts how the subsurface evolves during the course of a fluid-injection application. Using time-lapse geophysical signals as the data to be matched, the model unknowns to be estimated are the multiphysics forward-modeling parameters controlling the fluid-injection process. Properly reproducing the geophysical signature of the flow process, subsequent simulations can predict the fluid migration and alteration in the subsurface. The dynamic nature of fluid-injection processes renders imaging problems more complex than conventional geophysical imaging for static targets. This work intents to clarify the related hydrogeophysical parameter estimation concepts.

58 GEOSCIENCES↗