Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reservoir simulation”

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

Proxy-Based Bayesian Inversion Of Poroelastic Simulations To Interpret Strain Tensor Data Measured During Well Testing

The long runtimes of 3D poroelastic numerical simulators makes it impractical to interpret deformation datasets using many inversion schemes. Recent advances in instrumentation have made it possible to measure the strain tensor during well testing, but the lack of robust inversion methods is limiting the ability to interpret these data. We have developed an inversion workflow that reduces the number of computations required to complete a Bayesian inversion using DREAMzs. The workflow trains a KNN model using output from the poroelastic simulator, and then uses the KNN model as a proxy for the simulator during inversion. The workflow also includes a strategy for ensuring the results from the proxy model converge to the results from the simulator, ensuring the accuracy of the final results. An idealized example configured to represent a well test in a deep aquifer is used to verify that the workflow correctly identifies parameters and characterizes noise. Field data measured using strainmeters during an injection test at an oil reservoir in Oklahoma are used to evaluate performance with a real dataset. The workflow identified 265 history matching solutions out of 1240 total simulation runs (21% acceptance ratio), and the results are used to characterize posterior parameter distribution and evaluate the prediction uncertainty. This approach makes it feasible to invert strain data measured during well testing and this has the potential to improve the characterization of aquifers and reservoirs.

Roudini, Soheil↗
DOE OSTI2022

Quantifying Northern High Latitude Gross Primary Productivity (GPP) Using Carbonyl Sulfide (OCS)

The northern high latitude (NHL, 40°N to 90°N) is where the second peak region of gross primary productivity (GPP) other than the tropics. The summer NHL GPP is about 80% of the tropical peak, but both regions are still highly uncertain. Carbonyl sulfide (OCS) provides an important proxy for photosynthetic carbon uptake. Here we optimize the OCS plant uptake fluxes across the NHL by fitting atmospheric concentration simulation with the GEOS-CHEM global transport model to the aircraft profiles acquired over Alaska during NASA's Carbon in Arctic Reservoirs Vulnerability Experiment (2012–2015). We use the empirical biome-specific linear relationship between OCS plant uptake flux and GPP to derive the six plant uptake OCS fluxes from different GPP data. Such GPP-based fluxes are used to drive the concentration simulations. We evaluate the simulations against the independent observations at two ground sites of Alaska. The optimized OCS fluxes suggest the NHL plant uptake OCS flux of -247 Gg S year -1 , about 25% stronger than the ensemble mean of the six GPP-based OCS fluxes. GPP-based OCS fluxes systematically underestimate the peak growing season across the NHL, while a subset of models predict early start of season in Alaska, consistent with previous studies of net ecosystem exchange. The OCS optimized GPP of 34 PgC yr -1 for NHL is also about 25% more than the ensembles mean from six GPP data. Further work is needed to fully understand the environmental and biotic drivers and quantify their rate of photosynthetic carbon uptake in Arctic ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Above-cloud concentrations of cloud condensation nuclei help to sustain some Arctic low-level clouds

Abstract. Previous studies have found that low-level Arctic clouds often persist for long periods even in the face of very low surface cloud condensation nuclei (CCN) concentrations. Here, we investigate whether these conditions could occur due to continuous entrainment of aerosol particles from the free troposphere (FT). We use an idealized large eddy simulation (LES) modeling framework, where aerosol concentrations are low in the boundary layer (BL) but increased up to 50× in the free troposphere. We find that the tests with higher tropospheric aerosol concentrations simulated clouds, which persisted for longer and maintained higher liquid water paths (LWPs). This is due to direct entrainment of the tropospheric aerosol into the cloud layer, which results in a precipitation suppression from the increase in cloud droplet number and in stronger cloud-top radiative cooling, which causes stronger circulations maintaining the cloud in the absence of surface forcing. Together, these two responses result in a more well-mixed boundary layer with a top that remains in contact with the tropospheric aerosol reservoir and can maintain entrainment of those aerosol particles. The surface aerosol concentrations, however, remained low in all simulations. The free-tropospheric aerosol concentration necessary to maintain the clouds is consistent with concentrations that are frequently seen in observations.

Environmental Sciences & Ecology↗
NASA NTRS1991

Numerical analysis of rarefied slit flows. II - Navier-Stokes simulations

The model problem of pressure-driven flow of a rarefied monatomic gas through a two-dimensional slit is analyzed via full Navier-Stokes numerical simulation. Parametric solutions are generated for slit-height based Knudsen number ranging from continuum to transitional flow and for reservoir pressure ratios leading to subsonic and supersonic flow. The change in the structure of the flowfield near the slit as a function of pressure ratio and Knudsen number are quantified from a purely continuum standpoint. The choice of numerical domain size, boundary conditions and treatment of the slit are also discussed. As expected, comparison with a Direct Simulation Monte Carlo solution for a highly rarefied case shows large differences in the predicted mass flow. The cause of these differences can be quantified through detailed comparison of the local flowfield properties. For the larger pressure ratio cases qualitative trends with increasing rarefaction are discussed, including the change in the sonic line shape in the slit and the in total mass flow.

Wadsworth, D. C.↗
DOE OSTI2022

Carbon dioxide-enhanced metal release from kerogen

Heavy metals released from kerogen to produced water during oil/gas extraction have caused major environmental concerns. To curtail water usage and production in an operation and to use the same process for carbon sequestration, supercritical CO 2 (scCO 2 ) has been suggested as a fracking fluid or an oil/gas recovery agent. It has been shown previously that injection of scCO 2 into a reservoir may cause several chemical and physical changes to the reservoir properties including pore surface wettability, gas sorption capacity, and transport properties. Using molecular dynamics simulations, we here demonstrate that injection of scCO 2 might lead to desorption of physically adsorbed metals from kerogen structures. This process on one hand may impact the quality of produced water. On the other hand, it may enhance metal recovery if this process is used for in-situ extraction of critical metals from shale or other organic carbon-rich formations such as coal.

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

Utilizing GEOPHIRES-X Beyond Electricity

The GEOPHIRES tool is a techno-economic simulator for evaluating the thermal performance and cost-competitiveness of geothermal plants for electricity, heating, and/or cooling. The tool combines reservoir, wellbore, and surface plant cost and performance models to estimate overall techno-economic metrics such as net present value or levelized cost of electricity, heating, or cooling. We recently upgraded the tool to an object-oriented Python framework, presented in an accompanying paper. As part of the upgrade, we enhanced the capability to simulate the performance of geothermal plants for heating and cooling, which is the topic of this paper. Specifically, we (1) integrated absorption chillers to investigate the performance of utilizing geothermal heat for cooling, (2) integrated a heat pump module to boost the geothermal temperature and thermal output, (3) integrated a district heating module to estimate heating demand for a district based on local weather data, and simulated heat supply with geothermal energy and peaking boilers, and (4) integrated GEOPHIRES as an engine in the dGeo simulator to perform a geospatial analysis of geothermal district heating feasibility across a large region (e.g., a state or the entire United States) utilizing resource and thermal demand maps. This paper presents background information and case studies for several of these heating and cooling end-use options in GEOPHIRES.

absorption chiller↗
DOE OSTI2021

Simulation of Hurricane Harvey flood event through coupled hydrologic-hydraulic models: Challenges and next steps

Using the 2017 Hurricane Harvey flood event as a test case, this study set up a series of sensitivity analyses to highlight three challenges associated with large-scale flood inundation modeling, including (a) model parameterization, (b) errors in digital elevation models, and (c) effects of reservoir retention. Driven by radar-based hourly rainfall data, a series of hydrologic-hydraulic models including the VIC hydrologic model, RAPID routing model, and Flood2D-GPU hydrodynamic model are set up over Harris County, Texas, to simulate flood inundation and hazards. The results demonstrate the importance of hydrologic parameters in improving flood modeling. For a large flood event such as Hurricane Harvey, the effect of the initial water depths is insignificant. The Manning's n values may increase the peak water depth by ~1%, the flood extents by 65km 2 , and the high danger zone by ~6%. On the contrary, the bathymetry correction factors may reduce the flood extent by ~1.4% and the high-danger zone by ~4%. Reducing the reservoir storage capacity to 1% may increase the flood extent by ~4% and the high-danger zone by ~17%. This study may provide supporting information to guide and prioritize the development of future high-performance computing hydrodynamic large-scale flood simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Potential Formation Damage: An Integrated Reservoir Characterization Study of the Naturally Fractured Carbonate Middle Duperow Formation at the Kevin Dome, Montana

In this study, we integrate geologic and engineering data of a naturally fractured carbonate reservoir at the Kevin Dome, Montana. Well test data are correlated with core description, geochemical and lithology study to determine the flow behavior and communication within the injection test interval and to the surrounding area. Based on a dual-continuum geologic model, numerical brine injection simulations are carried out to validate the interpretation results from our well test analytical models and forecast the probability of CO 2 injection success using current reservoir properties. As a result, our well test analytical models as well as lithology/core description suggest that fluid flow may be mainly restricted to the injection interval and the assumption of radial (horizontal) flow may be appropriate. The well test models also indicate that there is potentially formation damage with a positive skin factor although prior to brine injection well tests, well stimulation through acid treatment was performed. Our numerical simulation results appear to confirm this formation damage by showing additional pressure buildup in the injection data during later test periods. To explain this, acid may have dissolved dolomite then dolomite or calcite may have been formed again further into the matrix/fracture system. Another possible explanation is mechanical clogging of the fractures due to acid dissolving dolomite and dislodging fine grains. Our work also predicts that if no additional well stimulation is performed, the project will have a lower probability of successfully injecting 1 million tons of CO 2 into the Middle Duperow formation over 4 years.

42 ENGINEERING↗
DOE OSTI2020

Self-Assembly of Silica Nanoparticles at Water–Hydrocarbon Interfaces: Insights from In Operando Small-Angle X-ray Scattering Measurements and Molecular Dynamics Simulations

To achieve tunable controls on the interactions of siliceous materials in subsurface environments bearing water and hydrocarbons, it is essential to determine the influence of hydrocarbon–water interfaces on the self-assembly of the silica particles. The hydrophilic silica particle have a tendency to aggregate on the water front of the interface. Self-assembly of silica nanoparticles proceeds via the migration of these nanoparticles to the water–hydrocarbon interface followed by aggregation at the interface. Fractal-like morphologies of assembled silica nanoparticles at water–hydrocarbon interfaces are observed. Rapid assembly of the hydrophilic silica nanoparticles at water–hydrocarbon interfaces corresponds to an overall reduction in the surface tension of water–toluene and water–heptane systems. These studies demonstrate the silica aggregation that is the precursor of silica polymerization, or nucleation and growth is influenced by the presence of hydrocarbons in subsurface geologic environments. These insights were derived from in operando ultrasmall and small-angle X-ray scattering (USAXS/SAXS) measurements, cryo-scanning electron microscopy (Cryo-SEM) imaging, and classical molecular dynamics simulations. As a result, these studies are intended to inform current and future efforts aimed at tuning silica reactivity in subsurface geologic environments to enhance permeability, the design and use of silica-based proppants to enhance fractures, and the development of effective strategies to control silica-based scaling behavior in subsurface reservoirs.

42 ENGINEERING↗
DOE OSTI2020

On the use of adjoints in the inversion of observed quasi-static deformation

An adjoint-based conjugate gradient algorithm provides an efficient means for imaging sources of deformation within the Earth, such as volume stresses associated with fluid flow in aquifers and reservoirs. For time intervals over which the overburden deforms elastically, one can calculate the gradient elements for a single model update using just two numerical simulations. The first is a forward run that is used to compute the residuals associated with the given iteration. The second simulation is to evaluate the application of the adjoint operator to the residuals. In this adjoint calculation, the residual displacements are applied as sources at the measurement locations, driving the deformation in the simulation. The volume stress on the source grid blocks, in response to the residual displacements, provide the gradient components. We apply this technique to satellite-based interferometric synthetic aperture radar (InSAR) line-of-sight displacements that were observed over an oil reservoir in California’s Central Valley. Furthermore, we find that the adjoint-based gradient estimates, requiring 18 CPU seconds, agree with conventional numerical calculations that take over 3700 CPU seconds to compute. Conjugate gradient algorithms utilizing the conventional approach and adjoint-based gradient computations give roughly the same reductions in misfit and similar final estimates of reservoir volume change.

58 GEOSCIENCES↗
DOE OSTI2023

Probabilistic Evaluation of Geomechanical Risks in CO2 Storage: An Exploration of Caprock Integrity Metrics Using a Multilaminate Model

The probabilistic uncertainty assessment of geomechanical risk—specifically, caprock failure—attributable to CO2 injection, as presented in a simplified hypothetical geological model, was the focus of this study. Our approach amalgamates the implementation of a multilaminate model, the creation of a response surface model in conjunction with the Box–Behnken sampling design, the execution of associated numerical modeling experiments, and the utilization of Monte Carlo simulations. Probability distributions to encapsulate the inherent variability (elastic and mechanical properties of the caprock and reservoir) and uncertainty in prediction estimates (vertical displacement, total strain, and F value) were employed. Our findings reveal that the Young modulus of the caprock is a key factor controlling equivalent total strain but is insufficient as a stand-alone indicator of caprock integrity. It is confirmed that the caprock can accommodate significant deformation without failure, if it possesses a low Young’s modulus and high mechanical strength properties, such as the friction angle and uniaxial compressive strength. Similarly, vertical displacement was found to be an unreliable indicator for caprock integrity, as caprock failure can occur across a broad spectrum of vertical displacements, particularly when both the Young modulus and mechanical strength properties have wide ranges. This study introduces the F value as the most dependable indicator for caprock failure, although it is a theoretical attribute (the shortest distance between the Mohr circle and the nearest failure envelope used to measure the sensitivity to failure) and not physically measurable in the field. Deviatoric stress levels were found to vary based on stress regimes, with the maximum levels observed under extensive and compressive stress regimes. In conjunction with the use of the response surface method, this study demonstrates the efficacy of the multilaminate framework and the Mohr–Coulomb constitutive model in providing a simplified, yet effective, probabilistic model of the mechanical behavior of caprock failure, reducing mathematical and computational complexities.

Energy & Fuels↗
NASA NTRS2007

Initial Testing of the Stainless Steel NaK-Cooled Circuit (SNaKC)

An actively pumped alkali metal flow circuit, designed and fabricated at the NASA Marshall Space Flight Center, is currently undergoing testing in the Early Flight Fission Test Facility (EFF-TF). Sodium potassium (NaK) was selected as the primary coolant. Basic circuit components include: simulated reactor core, NaK to gas heat exchanger, electromagnetic liquid metal pump, liquid metal flowmeter, load/drain reservoir, expansion reservoir, test section, and instrumentation. Operation of the circuit is based around the 37-pin partial-array core (pin and flow path dimensions are the same as those in a full core), designed to operate at 33 kWt. This presentation addresses the construction, fill and initial testing of the Stainless Steel NaK-Cooled Circuit (SNaKC).

Garber, Anne↗
DOE OSTI2024

A Novel Workflow for Coupled Simulation of Hydraulic Stimulation with Simultaneous Injection of Proppant

Enhanced Geothermal Systems (EGS) are located in impermeable host rock. Therefore, to achieve promising circulation, hydraulic stimulation is a key factor. Stimulation involves injecting a combination of fluids and proppant particles into existing fractures under high pressure. The effectiveness of stimulation depends on a variety of different factors, including the pressurization scheme and the optimal choice of carrying fluid and proppant particles to be used. For example, high-density particles can settle out of the injected slurry quickly, whereas lighter particles as well as a dense carrying fluid allow for a long transport into the reservoir. A numerical tool is therefore essential to take all these factors into account and to make robust predictions about the performance of the hydraulic stimulation. We present a novel computational application for coupled numerical simulation of hydraulic stimulation with concurrent injection of proppant in a fractured geothermal reservoir. This application is based on the open-source MOOSE framework that has been enriched with a set of equations to account for the circulation of proppant-laden slurries. It accounts for both the fluid and proppant components of the mixture, as well as particle-driven processes that include gravity settling, particle-particle interactions, and strong density and viscosity contrasts. A particular challenge is the mutual coupling with fracture mechanical processes in a propagating fracture using the FE code Defmod. Aperture and fluid pressure on one side and slurry density and viscosity changes on the other side interact and affect each other. Their basic understanding is therefore essential for the success of the stimulation. We compare our application to common numerical benchmark problems for both slurry flow and hydraulic stimulation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Simulation insights into wetting properties of hydrogen-brine-clay for hydrogen geo-storage

Hydrogen geo-storage is attracting substantial interdisciplinary interest as a cost-effective and sustainable option for medium- and long-term storage. Hydrogen can be stored underground in diverse formations, including aquifers, salt caverns, and depleted oil and gas reservoirs. The wetting dynamics of the hydrogen-brine-rock system are critical for assessing both structural and residual storage capacities, and ensuring containment safety. Through molecular dynamics simulations, we explore how varying concentrations of cushion gases (CO 2 or CH 4 ) influence the wetting properties of hydrogen-brine-clay systems under geological conditions (15 MPa and 333 K). We employed models of talc and the hydroxylated basal face of kaolinite (kaoOH) as clay substrates. Our findings reveal that the effect of cushion gases on hydrogen-brine-clay wettability is strongly dependent on the clay-brine interactions. Notably, CO 2 and CH 4 reduce the water wettability of talc in hydrogen-brine-talc systems, while exerting no influence on the wettability of hydrogen-brine-kaoOH systems. Detailed analysis of free energy of cavity formation near clay surfaces, clay-brine interfacial tensions, and the Willard-Chandler surface for gas-brine interfaces elucidate the molecular mechanisms underlying wettability changes. Our simulations identify empirical correlations between wetting properties and the average free energy required to perturb a flat interface when clay-brine interactions are less dominant. Here, our thorough thermodynamic analysis of rock-fluid and fluid-fluid interactions, aligning with key experimental observations, underscores the utility of simulated interfacial properties in refining contact angle measurements and predicting experimentally relevant properties. These insights significantly enhance the assessment of gas geo-storage potential. Prospectively, the approaches and findings obtained from this study could form a basis for more advanced multiscale simulations that consider a range of geological and operational variables, potentially guiding the development and improvement of geo-storage systems in general, with a particular focus on hydrogen storage.

25 ENERGY STORAGE↗
DOE OSTI2024

Techno-Economic Analysis and Market Potential of Geological Thermal Energy Storage (GeoTES) Charged With Solar Thermal and Heat Pumps

In this project, we developed a techno-economic analysis (TEA) model that can be used to evaluate the viability of a proposed Geological Thermal Energy Storage (GeoTES) design. This MATLAB-based model integrates distinct subsystem models for the reservoir, wells, power cycle, and solar field to capture their distinct characteristics. It applies this approach in simulating GeoTES storage and dispatch operations for durations ranging from hourly to seasonal. Using cases studies based on GeoTES designs provided by industry partners - Premier Resource Management (PRM) and EarthBridge Energy - we validated the TEA model estimations of system performance and costs (such as thermal and electrical power/energy inflow and outflow, capital costs, and levelized costs of energy and storage) for both concentrating solar thermal (CST) and Carnot Battery (CB) pairings with GeoTES (CST-GeoTES and CB-GeoTES). For the CST-GeoTES case, the model was validated against the proposed system designed by PRM. It showed good agreement with PRM's estimations when well and pump costs derived from PRM's estimations were used. When GETEM-based costs were used, there was a slight overprediction due to GETEM's project/site agnostic assumption of these costs. From a sensitivity analysis perspective, the levelized cost of electricity (LCOE) of the CST-GeoTES case was most sensitive to well flow rate and the charging temperature. An optimal design scenario resulted in an LCOE of 0.11 $\$$/kWhe. CST-GeoTES can also provide a source of heat to meet seasonal demands. With 12-hour and 24-hour levelized cost of heat (LCOH) of 0.018 $\$$/kWhth and 0.022 $\$$/kWhth, respectively, CST-GeoTES could be competitive in the California market with an average industrial price of natural gas in California between 0.041-0.047 $\$$/kWhth. The levelized cost of storage (LCOS) for CST-GeoTES depends on the energy storage duration. Although the LCOS is relatively higher for shorter durations (e.g., ~0.50 $\$$/kWhe for 1 hour of storage), it is an order of magnitude lower (0.06 $\$$/kWhe) for longer storage durations and competitive with lithium-ion batteries (beyond 12 hours of storage) and molten-salt thermal energy storage (beyond 32 hours). Energy. Three options were explored and applied to the EarthBridge case study: (1) A Carnot Battery design using R125 working fluid with both hot and cold storage; (2) A Carnot Battery design using R125 working fluid with only hot storage; (3) A Carnot Battery using a commercially available heat pump with carbon dioxide (CO2) working fluid and hot storage only. The CB-GeoTES with cold storage only had a slight (round-trip) efficiency advantage over the system without (43.4% vs. 42.8%). This is because the cold storage is limited by the freezing point of water, so the cold storage is not much colder than the environment. The system using commercially available technologies was the least efficient - partly because different cycles were used in the heat pump (CO2) and heat engine (binary cycle) which leads to some inefficiencies. Using the commercially available design, the levelized cost of energy (LCOS) from the model (0.10 $\$$/kWhe) was higher than that estimated by EarthBridge (0.068 $\$$/kWhe). This is because of the low round-trip (38.7%) efficiency of the commercially available design. Sensitivity analysis reveals that the model is most sensitive to electricity price. Including electricity price in the TEA for CB-GeoTES leads to an increase in LCOS from the base value to 0.25 $\$$/kWhe. To determine storage sites suitable for GeoTES, we gathered and analyzed geological, petrophysical, and geophysical data of oil and gas reservoir and aquifers in California and Texas. We down-selected possible sites based on cut-off values for site characteristics (e.g., reservoir temperature, formation thickness, permeability, porosity, depth, and brine salinity) and preliminary costs. Using this approach, the Carrizo-Wilcox, Yegua-Jackson, and Dockum brackish aquifers in Texas were identified as having the highest suitability. Similarly, in the central California region, the White Wolf, Belridge South Tulare, and Belridge South Reef Ridge were the most suitable. Going further, we assessed the storage potential in the selected sites. To do this we developed distributions of reservoir characteristic data and applied a Monte Carlo-based analysis to account for intrinsic uncertainty in the acquired data. The analysis revealed that the Carrizo-Wilcox aquifer had the highest storage potential with a mean capacity of 554 TWhth (i.e., 63 TWhe). The estimated capacity serves as an upper limit of storage potential given that not all fields in the basin will be developed. We participated in multiple outreach activities including conference presentations, panel session discussions, and the facilitation of a GeoTES workshop at the NREL Golden campus.

15 GEOTHERMAL ENERGY↗
NASA NTRS1998

Issues in the reconstruction of environmental doses on the basis of thermoluminescence measurements in the Techa riverside

The potential of thermoluminescence measurements of bricks from the contaminated area of the Techa river valley, Southern Urals, Russia, for reconstructing external exposures of affected population groups has been studied. Thermoluminescence dating of background samples was used to evaluate the age of old buildings available on the river banks. The anthropogenic gamma dose accrued in exposed samples is determined by subtracting the natural radiation background dose for the corresponding age from the accumulated dose measured by thermoluminescence. For a site in the upper Techa river region, where the levels of external exposures were extremely high, the depth-dose distribution in bricks and the dependence of accidental dose on the height of the sampling position were determined. For the same site, Monte Carlo simulations of radiation transport were performed for different source configurations corresponding to the situation before and after the construction of a reservoir on the river and evacuation of the population in 1956. A comparison of the results provides an understanding of the features of the measured depth-dose distributions and height dependencies in terms of the source configurations and shows that bricks from the higher sampling positions are likely to have accrued a larger fraction of anthropogenic dose from the time before the construction of the reservoir. The applicability of the thermoluminescent dosimetry method to environmental dose reconstruction in the middle Techa region, where the external exposure was relatively low, was also investigated.

NASA Discipline Radiation Health↗
DOE OSTI2021

Machine learning based rate optimization under geologic uncertainty

We propose a novel approach for rate optimization during a waterflood under geologic uncertainty in reservoir properties such as permeability and porosity. The traditional approach typically involves several runs of the forward simulator. This may not scale well when the optimization is to be performed at the full field-level and over multiple geologic realizations. A machine-learning (ML) based approach which is quick and scalable for rate optimization over multiple geologic realizations is proposed instead. The training data for the model is generated by running the forward simulator with randomly assigned well rates using multiple geologic realizations. A reduced order representation of the permeability heterogeneity in each of the realizations is derived using a grid connectivity transformation (GCT). This step involves finding basis functions corresponding to the different modal frequencies of the grid connectivity represented by the grid Laplacian. The projection of the heterogeneous property field along these basis functions gives the basis coefficients that form the reduced order representation. Subsequently, for each training datapoint, streamlines are traced and the minimum time of flight (TOF) representing the tracer breakthrough time at each producer is recorded. The basis coefficients and well rates are fed to a machine learning model as input and the minimum TOF at the producers forms the output of the model. This trained model can then be used along with an optimizer for computing the optimal injection rates to maximize the injection sweep efficiency. This corresponds to minimizing the variance in the minimum TOF within each well group. Different architectures of neural network are tested using 5-fold cross validation to decide the best ML model to compute the streamline time of flight. The trained model is used to perform well rate optimization over multiple realizations of geology by using a risk tolerance penalty. The optimal well rates thus obtained are compared with two cases: a) equal well rates assigned to all injectors and producers and b) well rates obtained by optimizing over a single realization without considering the uncertainty in geology. The optimal well rates are seen to offer better oil recovery and sweep efficiency than both cases.

02 PETROLEUM↗
NASA NTRSDate not supplied

Flight Test Exploration of Integrated Wildfire Response Operations with Crewed and Uncrewed Air Assets

Since 2016, Japan Aerospace Exploration Agency (JAXA) and NASA have partnered to investigate the safe and efficient integration of uncrewed aircraft systems (UAS) in disaster relief operations. Current disaster response operations would benefit from the use of UAS but are limited by the need for safe airspace integration technologies and enabling regulations. The work presented in this paper focuses on the safe integration of UAS and helicopters in active wildfire disaster response operations. The authors propose the use of mission planning and traffic management technologies, including onboard pilot situation awareness support technologies, to coordinate airspace operations between crewed and uncrewed aircraft (e.g., UAS). A representative wildfire scenario including a helicopter and UAS was developed. The helicopter mission included a flight path that mimicked the siphoning of water from a natural water reservoir (lake) and a subsequent water drop over vegetation along the progressing edge of a fire perimeter. Simulated UAS were assumed to conduct aerial ignition flights in the vicinity as part of a controlled burn of the vegetation ahead of the fire front to create a natural barrier between the fire and a nearby populated area. A situation awareness tool was developed and used to support the helicopter pilots’ task to stay clear of the volume of airspace where the UAS was conducting its operation while efficiently conducting the fire suppression mission. A flight test using JAXA’s experimental helicopter was conducted in August 2022 to validate the developed technology. The advantages of the developed technology were demonstrated, and pilot feedback confirmed that real-time awareness of the UAS operation airspace helped the flight planning and execution of a safe mission that minimized the flight time loss due to the potential interference with the UAS. Post-flight interviews with the pilots indicated no sufficient increase in the pilots’ workload and provided ideas for future technology development. Future research will provide recommendations to standards development organizations to support considerations for crewed-uncrewed interactions in disaster-response operations.

disaster response↗