Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reservoir flow model”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Influence of sequential stimulation practices on geochemical alteration of shale

Water-based hydraulic fracturing fluids (HFFs) can chemically interact with formation shale, resulting in altered porosity and permeability of the host rock. Experimental investigations of spatial and temporal shale-HFF interactions are helpful in interpreting chemical compositions of the injectate, as well as predicting alteration of hydraulic properties in the reservoir due to mineral dissolution and precipitation. Most bench-top experiments designed to study shale-HFF chemical interactions, either using batch reactors or flow-through setups, are carried out assuming that the acid spearhead has already become mixed with neutral HFFs. During operations, however, HFFs are typically injected according to a sequenced pumping schedule, starting with a concentrated acid spearhead, followed by multiple additions of near-neutral pH HFFs containing chemical amendments and proppant. In this study, we use geochemical modeling to consider whether this pre-mixed experimental protocol provides results directly comparable to a sequential discrete fluid-shale interaction protocol. Our results show that for the batch system, the transient evolution in major ion concentrations is faster with the sequential procedure. After 2 h of reaction time, the two protocols converge to the same aqueous concentrations. In a flow-through geometry, the pre-mixed model predicts extensive chemical alteration close to the injection point but negligible alteration downstream. In contrast, the sequential model predicts mineral reactions over hundreds of meters along the flow path. The extent of shale alteration in the sequential model at a given location depends on shale mineralogy and where the acid spearhead resides during the shut-in period. The predictive model developed in this study can help experimentalists to design bench-top tests and operators to better translate the results of laboratory experiments into practical applications.

Li, Qingyun↗
DOE OSTI2020

Simulation of coupled multiphase flow and geomechanics in porous media with embedded discrete fractures

In fractured natural formations, the equations governing fluid flow and geomechanics are strongly coupled. Hydrodynamical properties depend on the mechanical configuration, and they are therefore difficult to accurately resolve using uncoupled methods. In recent years, significant research has focused on discretization strategies for these coupled systems, particularly in the presence of complicated fracture network geometries. In this work, we explore a finite-volume discretization for the multiphase flow equations coupled with a finite-element scheme for the mechanical equations. Fractures are treated as lower dimensional surfaces embedded in a background grid. Interactions are captured using the embedded discrete fracture model (EDFM) and the embedded finite element method (EFEM) for the flow and the mechanics, respectively. This nonconforming approach significantly alleviates meshing challenges. EDFM considers fractures as lower dimension finite volumes that exchange fluxes with the rock matrix cells. The EFEM method provides, instead, a local enrichment of the finite-element space inside each matrix cell cut by a fracture element. Both the use of piecewise constant and piecewise linear enrichments are investigated. They are also compared to an extended finite element approach. One key advantage of EFEM is the element-based nature of the enrichment, which reduces the geometric complexity of the implementation and leads to linear systems with advantageous properties. Synthetic numerical tests are presented to study the convergence and accuracy of the proposed method. It is also applied to a realistic scenario, involving a heterogeneous reservoir with a complex fracture distribution, to demonstrate its relevance for field applications.

58 GEOSCIENCES↗
DOE OSTI2024

An integrated experimental–modeling approach to identify key processes for carbon mineralization in fractured mafic and ultramafic rocks

Abstract Controlling atmospheric warming requires immediate reduction of carbon dioxide (CO2) emissions, as well as the active removal and sequestration of CO2 from current point sources. One promising proposed strategy to reduce atmospheric CO2 levels is geologic carbon sequestration (GCS), where CO2 is injected into the subsurface and reacts with the formation to precipitate carbonate minerals. Rapid mineralization has recently been reported for field tests in mafic and ultramafic rocks. However, unlike saline aquifers and depleted oil and gas reservoirs historically considered for GCS, these formations can have extremely low porosities and permeabilities, limiting storage volumes and reactive mineral surfaces to the preexisting fracture network. As a result, coupling between geochemical interactions and the fracture network evolution is a critical component of long-term, sustainable carbon storage. In this paper, we summarize recent advances in integrating experimental and modeling approaches to determine the first-order processes for carbon mineralization in a fractured mafic/ultramafic rock system. We observe the critical role of fracture aperture, flow, and surface characteristics in controlling the quantity, identity, and morphology of secondary precipitates and present where the influence of these factors can be reflected in newly developed thermo-hydro-mechanical–chemical models. Our findings provide a roadmap for future work on carbon mineralization, as we present the most important system components and key challenges that we are overcoming to enable GCS in mafic and ultramafic rocks.

58 GEOSCIENCES↗
DOE OSTI2026

A path to intelligent watersheds: coordinating the data to decision pipeline

Operations of multi-reservoir systems are challenged in-part by the interplay of complex physical processes functioning within the watershed. The employment of intelligent systems can be of aid by linking environmental sensing, information technology, data analytics, simulation and decision support to achieve a data-to-decision flow of information. A further challenge is that watershed resources are managed for multiple purposes requiring some level of coordination among numerous resource managers, asset operators and users. System intelligence in this context relies on shared community platforms (data portals, community models), and coordinated communication between decision makers. Opportunities to enrich watershed intelligence has been the subject of a roadmapping exercise for the Department of Energy’s Water Power Technologies Office which has relied on broad stakeholder engagement. Initial phases of engagement involved personal interviews and a series of virtual group meetings, which focused on identifying opportunities to improve the intelligence of the physical infrastructure within our watersheds—examples of feedback include improved sensing of snowpack and runoff, data standards for facilitated data sharing, and better forecasting tools. The latter phase of engagement involved the conduct of a case study in the Upper Colorado River basin where key stakeholders were interviewed to map how their decisions are informed by intelligence from other basin stakeholders. Our presentation will highlight the interdisciplinary flow of information in complex watershed systems and identify physical and institutional opportunities toward the strategic operation of water infrastructure.

Colorado River↗
DOE OSTI2020

Multi‐scale simulation of rock compaction through breakage models with microstructure evolution

Abstract. Regional subsidence due to fluid depletion includes the interaction among multiple physical processes. Specifically, rock compaction is governed by coupled hydro-mechanical feedbacks involving fluid flow, effective stress change and pore collapse. Although poroelastic models are often used to explain the delay between depletion and subsidence, recent evidence indicates that inelastic effects could alter the rock microstructure, thus exacerbating coupling effects. Here, a constitutive law built within the framework of Breakage Mechanics is proposed to account for the inherent connection between rock microstructure, hydraulic conductivity, and pore compaction. Furthermore, it is embedded into a 1-D hydromechanical coupled finite element analysis (FEA) to explore the effects of micro-structure rearrangement on the development of reservoir compaction. Numerical examples with the proposed model are compared with simulations under constant hydraulic conductivity to illustrate the model capability to capture the non-linear processes of reservoir compaction induced by fluid depletion.

Buscarnera, Giuseppe↗
DOE OSTI2020

Using Machine Learning to Predict Future Temperature Outputs in Geothermal Systems

Optimizing the power output, and economic value, of geothermal power plants over decades of operation is a major challenge in renewable energy. Optimizing the output requires the ability to predict the mass flow rates and the output temperatures of production wells based on the inputs of injection wells, as well as the time history of the system. Machine Learning (ML) that incorporates the known physics of geothermal systems is one possible solution to this challenge. In this work, we explore the ability of ML algorithms to predict future temperature outputs based on historical data. Considering the challenges with obtaining an empirical dataset from field data that is large enough to enable reliable ML, we propose an alternate approach: developing a high-fidelity reservoir model and using computational resources to build a dataset that enables ML. As a first step towards achieving this goal, we present preliminary results from applying ML to predict the temperature timeseries of simple modeled geothermal systems. We describe the application of relevant state-of-the-art ML approaches, such as the Long Short-Term Memory (LSTM) networks and Convolutional Neural Networks (CNN), to extract temporal structures in the model data. We assess the accuracy of the forecasts we obtain, compare the selected approaches, and share the lessons learned that would inform the process of training and utilizing ML algorithms for larger and more complex geothermal systems.

GEOTHERMAL ENERGY↗
DOE OSTI2021

Predicting Cement-Sheath Integrity with Consideration of Initial State of Stress and Thermoporoelastic Effects

In conventional wellbore-integrity analysis, the cement sheath’s initial state of stress and transient thermoporoelastic effects are often neglected. However, the initial state of stress is prerequisite information for accurately predicting the safe operating conditions that prevent a cemented well from being damaged. In addition, transient thermoporoelastic effects can have a profound effect on when damage will occur. Here in this paper, we propose a model that includes these effects to predict the safe operating pressures and temperatures that will prevent cement-sheath failure. For the initial state of stress, we proposed an empirical model using measurements. Subsequent stress changes are evaluated by a fully coupled transient thermoporoelastic model to analyze the mechanical behavior of the cement sheath. We predict the safe operating envelope (SOE) for shear, tensile, and debonding cement-sheath failures caused by pressure and temperature perturbations after the cement sets. Our model predicts that pore pressure is a key factor for cement failure, especially for rapid temperature changes. If the formation is low permeability, the transient pore pressures are amplified, causing the risk of damage to increase. Compared with conventional thermoelastic models, the thermoporoelastic model predicts a smaller SOE when heating the internal casing fluid and a larger envelope when cooling the internal casing fluid. Finally, the heating rate was considered with respect to field applications. The heating rate was also considered, and slower heating/cooling rates can prevent damage to the cement sheath. Finally, the thermoporoelastic model was applied to explain several laboratory and field experiments and achieved good matches.

02 PETROLEUM↗
DOE OSTI2022

Thermal and solubility effects on fault leakage during geologic carbon storage

Geologic carbon storage (GCS) is a promising method for reducing anthropogenic CO 2 emissions to the atmosphere. To safely deploy GCS in the field, it is necessary to assess risks and the effect of uncertainty on safe storage. The effect of uncertainty can be quantified using batches of simulations, but the high computational costs of high-resolution simulations necessitate use of reduced-order models (ROMs). Previous work involves ROMs for quantifying the risk of different potential leakage paths from storage reservoirs to shallow formations. However, previous studies on development of fault-leakage ROMs have limited numbers of uncertain parameters and do not explicitly examine impacts of CO 2 solubility and thermal stresses on fault reactivation, which can generate high-permeability pathways and compromise CO 2 storage. Here, we analyze an ensemble of simulations considering CO 2 leakage from a storage reservoir to a shallow aquifer through a fault while varying a number of uncertain parameters related to thermo-hydro-mechanical properties and CO 2 injection. We show the effects of solubility on: free-phase CO 2 -leakage rates, brine-leakage rates, and poroelastic fault destabilization. We find that CO 2 solubility is more important for estimating free-phase CO 2 -leakage rates compared to brine-leakage rates or poroelastic fault destabilization. We also find that thermal stresses and overpressures have different spatial distributions within the fault, indicating that the spatial variability of overpressures due to variation in flow parameters does not necessarily make the spatial variability of thermal stresses negligible. We suggest the use of the CO 2 phase-change path as a variable in future fault-leakage ROMs.

03 NATURAL GAS↗
DOE OSTI2022

The Source Permeator System and Tritium Transport in the TEX PbLi Loop

In this work, permeation is investigated for the introduction of hydrogen isotopes into PbLi for the Tritium Extraction eXperiment (TEX). TEX is a forced-convection PbLi loop under construction at Idaho National Laboratory (INL) that will test the vacuum permeator (VP) method of tritium extraction from PbLi. The source permeator (SP) delivers atomic hydrogen (H, D, and T) from a gas-phase reservoir into the PbLi via a permeable dense metal membrane. A modular and a fixed source permeator system are presented. In the modular design, PbLi flows through the inside of a tubular membrane and gas-phase hydrogen is introduced on the outside of the membrane. Atomic hydrogen permeates radially inward through the membrane into the PbLi. In the fixed design, PbLi flows into an expansion chamber with closed-ended tubular membranes inserted. Gas-phase hydrogen is introduced on the inside of the closed-ended membranes and atomic hydrogen permeates radially outward into the flowing PbLi. Hydrogen-transport models based on steady-state mass transport through PbLi and permeation through the metal membrane were developed to assess the operation of the source permeator relative to experimental variables and to allow understanding of uncertain parameter effects, such as PbLi hydrogen-transport properties and the effective hydrogen permeability of the vacuum permeator. This modeling effort considers iron as the source permeator material and vanadium as the vacuum permeator material.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Warming of the Willamette River, 1850–present: the effects of climate change and river system alterations

Abstract. Using archival research methods, we recovered and combined data from multiple sources to produce a unique, 140-year record of daily water temperature (Tw) in the lower Willamette River, Oregon (1881–1890, 1941–present). Additional daily weather and river flow records from the 1850s onwards are used to develop and validate a statistical regression model of Tw for 1850–2020. The model simulates the time-lagged response of Tw to air temperature and river flow and is calibrated for three distinct time periods: the late 19th, mid-20th, and early 21st centuries. Results show that Tw has trended upwards at 1.1 ∘C per century since the mid-19th century, with the largest shift in January and February (1.3 ∘C per century) and the smallest in May and June (∼ 0.8 ∘C per century). The duration that the river exceeds the ecologically important threshold of 20 ∘C has increased by about 20 d since the 1800s, to about 60 d yr−1. Moreover, cold-water days below 2 ∘C have virtually disappeared, and the river no longer freezes. Since 1900, changes are primarily correlated with increases in air temperature (Tw increase of 0.81 ± 0.25 ∘C) but also occur due to alterations in the river system such as depth increases from reservoirs (0.34 ± 0.12 ∘C). Managed release of water affects Tw seasonally, with an average reduction of up to 0.56 ∘C estimated for September. River system changes have decreased variability (σ) in daily minimum Tw by 0.44 ∘C, increased thermal memory, reduced interannual variability, and reduced the response to short-term meteorological forcing (e.g., heat waves). These changes fundamentally alter the response of Tw to climate change, posing additional stressors on fauna.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Temperature-responsive smart tracers for field-measurement of inter-well thermal evolution: Heterogeneous kinetics and field demonstration

Temperature-responsive smart tracers enable advanced warning of a “premature thermal breakthrough” and may therefore improve reservoir management and reduce financial risk. A successful calculation of inter-well fluid temperatures requires that temperature-dependent kinetics are known, which typically result from homogeneous batch reactor experiments. However, recent meso-scale field experiments at the Altona Field Laboratory involving fluid flow in a discretely-fractured reservoir suggest that silica-fluid interfaces may accelerate hydrolysis kinetics. Here, the Arrhenius parameters of phenyl acetate hydrolysis are investigated under heterogeneous reaction conditions in packed-bed column experiments. The breakthrough curve of the reaction product is compared to an inert reference tracer (carbon-cored nanoparticles). Temperatures experienced during field testing ranging from 10 to 40 °C were studied using a phosphate buffer solution to maintain a near-neutral pH, as found in the reservoir. The empirically-determined pre-exponential factor and activation energy values are subsequently used in a first-order kinetic model to improve calculations of effective reservoir temperatures for these meso-scale field tests. Here, the results suggest that hydrolysis reactions can directly monitor inter-well reservoir temperatures, but only if the influence of solid-fluid interactions are carefully accounted for.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Fracture Sustainability in Enhanced Geothermal Systems: Experimental and Modeling Constraints

Enhanced geothermal systems (EGS) offer the potential for a much larger energy source than conventional hydrothermal systems. Hot, low-permeability rocks are prevalent at depth around the world, but the challenge of extracting thermal energy depends on the ability to create and sustain open fracture networks. Laboratory experiments were conducted using a suite of selected rock cores (granite, metasediment, rhyolite ash-flow tuff, and silicified rhyolitic tuff) at relevant pressures (uniaxial loading up to 20.7 MPa and fluid pressures up to 10.3 MPa) and temperatures (150–250 °C) to evaluate the potential impacts of circulating fluids through fractured rock by monitoring changes in fracture aperture, mineralogy, permeability, and fluid chemistry. Because a fluid in disequilibrium with the rocks (deionized water) was used for these experiments, there was net dissolution of the rock sample: this increased with increasing temperature and experiment duration. Thermal-hydrological-mechanical-chemical (THMC) modeling simulations were performed for the rhyolite ash-flow tuff experiment to test the ability to predict the observed changes. These simulations were performed in two steps: a thermal-hydrological-mechanical (THM) simulation to evaluate the effects of compression of the fracture, and a thermal-hydrological-chemical (THC) simulation to evaluate the effects of hydrothermal reactions on the fracture mineralogy, porosity, and permeability. Furthermore, these experiments and simulations point out how differences in rock mineralogy, fluid chemistry, and geomechanical properties influence how long asperity-propped fracture apertures may be sustained. Such core-scale experiments and simulations can be used to predict EGS reservoir behavior on the field scale.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Pore-Scale Dynamics of Liquid CO2–Water Displacement in 2D Axisymmetric Porous Micromodels Under Strong Drainage and Weak Imbibition Conditions: High-Speed μPIV Measurements

Resolving pore-scale transient flow dynamics is crucial to understanding the physics underlying multiphase flow in porous media and informing large-scale predictive models. Surface properties of the porous matrix play an important role in controlling such physics, yet interfacial mechanisms remain poorly understood, in part due to a lack of direct observations. This study reports on an experimental investigation of the pore-scale flow dynamics of liquid CO 2 and water in two-dimensional (2D) circular porous micromodels with different surface characteristics employing high-speed microscopic particle image velocimetry (μPIV). The design of the micromodel minimized side boundary effects due to the limited size of the domain. The high-speed μPIV technique resolved the spatial and temporal dynamics of multiphase flow of CO 2 and water under reservoir-relevant conditions, for both drainage and imbibition scenarios. When CO 2 displaced water in a hydrophilic micromodel (i.e., drainage), unstable capillary fingering occurred and the pore flow was dominated by successive pore-scale burst events (i.e., Haines jumps). When the same experiment was repeated in a nearly neutral wetting micromodel (i.e., weak imbibition), flow instability and fluctuations were virtually eliminated, leading to a more compact displacement pattern. Energy balance analysis indicates that the conversion efficiency between surface energy and external work is less than 30%, and that kinetic energy is a disproportionately smaller contributor to the energy budget. This is true even during a Haines jump event, which induces velocities typically two orders of magnitude higher than the bulk velocity. These novel measurements further enabled direct observations of the meniscus displacement, revealing a significant alteration of the pore filling mechanisms during drainage and imbibition. While the former typically featured burst events, which often occur only at one of the several throats connecting a pore, the latter is typically dominated by a cooperative filling mechanism involving simultaneous invasion of a pore from multiple throats. This cooperative filling mechanism leads to merging of two interfaces and releases surface energy, causing instantaneous high-speed events that are similar, yet fundamentally different from, burst events. Finally, pore-scale velocity fields were statistically analyzed to provide a quantitative measure of the role of capillary effects in these pore flows.

Li, Yaofa↗
DOE OSTI2021

Experimental study of three-dimensional CO 2 -water drainage and fracture-matrix interactions in fractured porous media

Estimation of CO 2 storage capacity in fractured porous reservoirs requires a better understanding of the CO 2 -water flow fundamentals in fracture-matrix systems. There are few core-flood experiments on three-dimensional CO 2 -water drainage in a fracture-matrix system, yet none have examined the impacts of the capillary continuity of fractures and fracture-matrix interactions. In this study, twelve drainage experiments were conducted in four fracture-matrix columns, each comprising a vertical stack of a cylindrical rock core, a filter paper serving as an analogue to a horizontal fracture, and a ceramic plate. Here, the core sample was surrounded by open space at the top and circumferential sides to model fractures, allowing for 3-D CO 2 -water drainage in the rock core and displaced water draining across the horizontal fracture. X-ray computed tomography was conducted to visualize the dynamic invasion/drainage processes in four rock core samples showing contrasts in anisotropy, permeability, and heterogeneity. Experimental results show (1) the equilibrium CO 2 saturations in the rock matrix vary from 0.10 to 0.60 at controlled capillary pressures up to 200 kPa, (2) the CO 2 saturations in the matrix increase with water saturation in the horizontal fracture resulting in better capillary continuity for water to drain across the fracture, and (3) the fracture water saturation can be enhanced by the non-uniform fracture capillary pressure and countercurrent flow of CO 2 and water across the fracture-matrix interface. In the case of high fracture water saturation, matrix CO 2 saturation largely depends on matrix anisotropy and heterogeneity. The core-scale experimental results contribute to understand the fracture-matrix interactions and CO 2 storage efficiency in fractured porous media.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Isotopic fractionation as in-situ sensor of subsurface reactive flow and precursor for rock failure

Greater utilization of subsurface reservoirs perturbs in-situ chemical-mechanical conditions with wide ranging consequences from decreased performance to project failure. Understanding the chemical precursors to rock deformation is critical to reducing the risks of these activities. To address this need, we investigated the coupled flow-dissolution- precipitation-adsorption reactions involving calcite and environmentally-relevant solid phases. Experimentally, we quantified (1) stable isotope fractionation processes for strontium during calcite nucleation and growth, and during reactive fluid flow; (2) consolidation behavior of calcite assemblages in the common brines. Numerically, we quantified water weakening of calcite using molecular dynamics simulations; and quantified the impact of calcite dissolution rate on macroscopic fracturing using finite element models. With microfluidic experiments and modeling, we show the effect of local flow fields on the dissolution kinetics of calcite. Taken together across a wide range of scales and methods, our studies allow us to separate the effects of reaction, flow, and transport, on calcite fracturing and the evolution of strontium isotopic signatures in the reactive fluids.

42 ENGINEERING↗
DOE OSTI2021

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

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

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Performance of Closed Loop Venturi Aspirated Aeration System: Experimental Study and Numerical Analysis with Discrete Bubble Model

In wastewater treatment plants, aeration plays a significant role as it encourages aerobic respiration of microbes, which are necessary to break down carbonaceous matter in the waste stream. This process can account for the majority of energy use in wastewater treatment plants. The aeration process is also necessary in odor control in lagoons and in the aquaculture industry. Generally, the aeration process is accomplished with compressors or blowers which can be of low efficiency due to ideal gas laws. This study introduces the idea of increasing aeration efficiency by looping water from a reservoir through a piping network which includes a venturi aspirator at its inlet. For this purpose, an experimental study has been conducted in a laboratory setup with a pump which pulls water from a tank, passes it through a Venturi aspirator and a helical piping network intended to increase bubble residence time, before depositing it back into the bulk fluid tank. This same process is modeled computationally, using a discrete bubble method (DBM), with good agreement with experiments. The overall purpose here is to determine the optimal configurations for standard aeration efficiency (SAE) and the standard oxygen transfer rate (SOTR). A parametric study has been implemented using the DBM based on different hydraulic and flow parameters. The model is also used to predict the SAE of a hypothetical aeration system. Results indicate that it is possible to achieve SAE values in the range of surface aerators or submersed jet aerators using the proposed aeration system with less complex components, thereby decreasing overall costs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
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↗