Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reservoir modeling”

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 343 records · Page 19

Io's heat flow from infrared radiometry: 1983-1993

We report the following results from a decade of infrared radiometry of Io: (1) The average global heat flow is more than approx. 2.5 W/sq.m, (2) large warm (less than or equal to 200 K) volcanic regions dominate the global heat flow, (3) smal high-temperature (greater than or = 300 K) 'hotspots' contribute little to the average heat flow, (4) thermal anomalies on the leading hemisphere contribute about half of the heat flow, (5) a substantial amount of heat is radiated during Io's night, (6) high-temperature (greater than or = 600 K) 'outbursts' occurred during approx. 4% of the nights we observed, (7) 'Loki' is the brightest, persistent, infrared emission feature, and (8) some excess emission is always present at the longitude of Loki, but its intensity and other characteristics change between apparitions. Observations of Io at M(4.8 micrometer), 8.7 micrometer, N(10 micrometer), and Q(20 micrometer) with the Infrared Telescope Facility presented here were collected during nine apparitions between 1983 and 1993. These measurements provide full longitudinal coveraged as well as an eclipse observation and the detection of two outbursts. Reflected sunlight, passive thermal emission, and radiation from thermal anomalies all contribute to the observed flux densities. We find that a new thermophysical model is required to match all the data. Two key elements of this model are (1) a 'thermal reservoir' unit which lowers daytime temperatures, and (2) the 'thermal pedestal effect' which shifts to shorter wavelengths the spectral emission due to the reradiation of solar energy absorbed by the thermal anomalies. The thermal anomalies are modeled with a total of 10 source components at five locations. Io's heat flow is the sum of the power from these components.

Veeder, Glenn J.↗

Quantifying Water Storage Change and Land Subsidence Induced by Reservoir Impoundment Using GRACE, Landsat, and GPS Data

The construction of hydropower dams is a common strategy to support a country's increasing need for electricity and river water management for industry and agriculture. Although the hydrological and geophysical impacts of water relocation are usually assessed prior to impoundment, their accuracy is generally limited due to the lack of in situ observations, especially in a remote area. This study presents a workflow to quantify the terrestrial water storage change (TWS) and land subsidence induced by a reservoir's water impoundment using multiple satellite observations (GRACE, Landsat), land surface models (CABLE, GLDAS, NCEP, ECMWF), and GPS data. The study site is the Bakun Dam, located in Sarawak, Malaysia, which is the largest hydropower dam in Southeast Asia. Commencing operation in late 2010, the dam induced a change of water mass and lake surface area that was clearly observed by GRACE and Landsat observations, respectively. During the 17-month impounding period (from August 2010 to December 2011), GRACE observed a dramatic increase of approximately 200 mm equivalent water height, while Landsat detected an increased lake extent of around 600 km2. In this paper, a forward model is developed to determine the increased water surface level corresponding to GRACE observations, estimated to be about 120 m. In contrast to GRACE, the TWS derived from land surface models cannot capture the increased TWS, due to the lack of reservoir routing algorithms in the models. In addition, the land subsidence was calculated using the disk load model constructed based on the GRACE-derived lake level and Landsat-derived lake extent; the result is validated with the GPS data from BIN1 station, located at the western coast of Borneo. The commencement stage of the Bakun Dam induces the large-scale land subsidence, which causes the GPS-BIN1 station to subside by ~9 mm, and move toward the Bakun Lake by ~4 mm. Computation of the surface displacements directly from GRACE spherical harmonic coefficient data fails to capture the subsidence feature, mainly due to the truncation error. Overall, this study demonstrates that evaluating GRACE in conjunction with Landsat, LSMs, and GPS data allows the exploitation of the gravity signal at a much smaller spatial scale than its intrinsic resolution. Benefiting from global coverage, the newly developed satellite-based algorithm is a valuable tool for assessing the impacts of reservoir operation on hydrological and geophysical changes from local to regional scales.

Natthachet Tangdamrongsub↗

Assessment of two neglected effects in the analysis of an oil pumping ring seal

Many factors have been found to affect the performance of a pumping ring seal. In this paper, two effects (elevated reservoir temperature and rod deformation), both heretofore neglected, are assessed through the use of a thermoelastohydrodynamic numerical model of the pumping ring. Elevated reservoir temperatures are found to result in an increase in the amount of lubricant pumped while deformation of the translating rod is shown to cause a reduction in the lubricant pumped.

Keith, T. G.↗

Williston Basin Resource Study for Commercial-Scale Subsurface Hydrogen Storage

This closeout presentation summarizes a Department of Energy–funded study that evaluated whether large amounts of hydrogen can be safely stored underground in the North Dakota portion of the Williston Basin. The project combined lab testing, computer simulations, and basin‑wide analysis to assess saline formations, depleted oil and gas reservoirs, and salt formations for hydrogen storage capacity, recovery, and risks. Results show that underground hydrogen storage is technically feasible across multiple formation types, with depleted oil and gas reservoirs offering higher recovery and saline formations providing large long‑term storage potential. The study also identifies key challenges—such as wellbore material durability and gas purity management—and recommends pilot projects and further site‑specific studies to support future commercialization.

08 HYDROGEN↗

Thermal Management Tools for Propulsion System Trade Studies and Analysis

Energy-related subsystems in modern aircraft are more tightly coupled with less design margin. These subsystems include thermal management subsystems, vehicle electric power generation and distribution, aircraft engines, and flight control. Tighter coupling, lower design margins, and higher system complexity all make preliminary trade studies difficult. A suite of thermal management analysis tools has been developed to facilitate trade studies during preliminary design of air-vehicle propulsion systems. Simulink blocksets (from MathWorks) for developing quasi-steady-state and transient system models of aircraft thermal management systems and related energy systems have been developed. These blocksets extend the Simulink modeling environment in the thermal sciences and aircraft systems disciplines. The blocksets include blocks for modeling aircraft system heat loads, heat exchangers, pumps, reservoirs, fuel tanks, and other components at varying levels of model fidelity. The blocksets have been applied in a first-principles, physics-based modeling and simulation architecture for rapid prototyping of aircraft thermal management and related systems. They have been applied in representative modern aircraft thermal management system studies. The modeling and simulation architecture has also been used to conduct trade studies in a vehicle level model that incorporates coupling effects among the aircraft mission, engine cycle, fuel, and multi-phase heat-transfer materials.

McCarthy, Kevin↗

Tracking the Impacts of Precipitation Phase Changes Through the Hydrologic Cycle in Snowy Regions: From Precipitation to Reservoir Storage

Cool season precipitation plays a critical role in regional water resource management in the western United States. Throughout the twenty-first century, regional precipitation will be impacted by rising temperatures and changing circulation patterns. Changes to precipitation magnitude remain challenging to project; however, precipitation phase is largely dependent on temperature, and temperature predictions from global climate models are generally in agreement. To understand the implications of this dependence, we investigate projected patterns in changing precipitation phase for mountain areas of the western United States over the twenty-first century and how shifts from snow to rain may impact runoff. We downscale two bias-corrected global climate models for historical and end-century decades with the Weather Research and Forecasting (WRF) regional climate model to estimate precipitation phase and spatial patterns at high spatial resolution (9 km). For future decades, we use the RCP 8.5 scenario, which may be considered a very high baseline emissions scenario to quantify snow season differences over major mountain chains in the western U.S. Under this scenario, the average annual snowfall fraction over the Sierra Nevada decreases by >45% by the end of the century. In contrast, for the colder Rocky Mountains, the snowfall fraction decreases by 29%. Streamflow peaks in basins draining the Sierra Nevada are projected to arrive nearly a month earlier by the end of the century. By coupling WRF with a water resources model, we estimate that California reservoirs will shift towards earlier maximum storage by 1–2 months, suggesting that water management strategies will need to adapt to changes in streamflow magnitude and timing.

Melissa L. Wrzesien↗

Facies Analysis of the Prairie Du Chien Group in the Illinois Basin and Analogous Rocks in Missouri and Kentucky

Funded in 2023 by the U.S. Department of Energy’s Phase II Carbon Storage Assurance Facility Enterprise (CarbonSAFE) initiative, a Heidelberg Materials cement plant in Mitchell, Indiana, is currently being evaluated as a potential Carbon Capture and Storage (CCS) subsurface injection site. The Heidelberg CCS project targets the middle to upper Prairie du Chien Group (Early Ordovician) in southwestern Indiana. Assessment of reservoir feasibility requires collection of field data, seismic surveys, well-log correlation, geologic modeling, characterization well drilling, well testing, and reservoir simulation. However, the proposed Heidelberg CCS site is in a data-limited region, lacking both outcrop analogs and deep wells penetrating the target interval, which makes geologic modelling difficult prior to drilling a characterization well. To directly address this problem, the present study was undertaken to understand the sedimentologic composition and stratigraphic architecture of the Prairie du Chien Group from analogous outcrops and cores in the Illinois Basin and adjacent regions.

Ali, Shah Bilawal [Univ. of Illinois at Urbana-Cha↗

Evaluating proxies for the drivers of natural gas productivity using machine-learning models

We report the extensive development of unconventional reservoirs using horizontal drilling and multistage hydraulic fracturing has generated large volumes of reservoir characterization and production data. The analysis of this abundant data using statistical methods and advanced machine-learning (ML) techniques can provide data-driven insights into well performance. Most predictive modeling studies have focused on the impact that different well completion and stimulation strategies have on well production but have not fully exploited the available in situ rock property data to determine its role in reservoir productivity. We have used machine-learning techniques to rank rock mechanical properties, microseismic attributes, and stimulation parameters in the order of their significance for predicting natural gas production from an unconventional reservoir. The data for this study came from a hydraulically fractured well in the Marcellus Shale in Monongalia County, West Virginia. The data classes included measurements aggregated by well completion stage that included (1) gas production, (2) well-log-derived measurements including bulk density, elastic moduli, shear impedance, compressional impedance, brittleness, and gamma measurements, (3) microseismic attributes, (4) long-period long-duration (LPLD) event counts, (5) fracture counts, and (6) stimulation parameters that included the fluid injection volume and average pumping pressure. To identify observable proxies for the drivers of gas production, we evaluated five commonly used ML approaches including multivariate adaptive regression spline, Gaussian mixture model, random forest, gradient boosting, and neural network. We selected five variables including LPLD event count, seismogenic b-value, hydraulic diffusivity, cumulative moment, and fluid volume as the features most likely to impact gas productivity at the stage level in the study area. The data-driven selection of these parameters for their importance in determining gas production can help reservoir engineers design more effective hydraulic-fracture treatments in the Marcellus Shale and other similar unconventional reservoirs. Plain language summary: We use machine-learning methods and data-driven selection of reservoir parameters to rank and better understand their importance in determining gas production, which can help reservoir engineers design more effective hydraulic-fracture treatments in the Marcellus Shale and other similar unconventional reservoirs.

58 GEOSCIENCES↗

Zircon Lu-Hf systematics: Evidence for the episodic development of Archaean greenstone belts

A combined U-Th-Pb and Lu-Hf isotopic study of zircons was undertaken in order to determine the provenance and age of an Archean granite-greenstone terrain and to test the detailed application of the Lu-Hf system in various Archean zircons. The eastern Wawa subprovince of the Superior province consists of the low grade Michipicoten and Gamitagama greenstone belts and the granitic terrain. The Hf isotopic data indicate that the typical lithological features of a greenstone belt cycle could be accommodated in a crustal growth model that involved decreasing depth of melting in three isotopically distinct reservoirs: mantle, lower crust and upper crust. The model age of the sources given by the intersection of the lower crustal curve with the bulk earth evolution curve is about 2900 My, in good agreement with the zircon U-Pb basement age. This linear array also has a similar intersection age to that of Proterozoic carbonatite complexes. The general convergence of the other reservoir vectors around this age suggests that mantle depletion, crustal extraction and intracrustal differentiation were all part of the same episodic event. It is also apparent that recycling of older basement was important in the formation of many of the later greenstone belt rocks.

Smith, P. E.↗

SNC Meteorites and Martian Reservoirs

Jones first suggested that the inverse covariation of initial epsilon (Nd-143) and Sr-87/Sr-86 of the shergottites could be explained by interaction between mantle-derived magmas with another isotopic reservoir(s) (i.e., assimilation or contamination). In that model, magmas were generated in a source region that was isotopically similar to the Nakhla source and the second reservoir(s) was presumed to be crust. The text also permitted the second reservoir to be another type of mantle, but I can confirm that a second mantle reservoir was never seriously considered by that author. Other features of this model were that (i) it occurred at a particular time, 180 m.y. ago, and (ii) the interacting reservoirs had been separated at approximately 4.5 b.y. In a later paper Jones explored this mixing model more quantitatively and concluded that magmas from a Nakhla-like source region at 180 m.y. would fall on or near an isotopic Nd-Sr-Pb hyperplane defined by the shergottites. This criterion was a necessary prerequisite for the parent magma(s) of the shergottites to have initially been Nakhla-like isotopically. At this juncture, it is perhaps worthwhile to note that this mixing model was not presented to explain geochemical variations but as a justification for a 180 m.y. crystallization age for the shergottites and a 1.3 b.y. crystallization age for the nakhlites. In the mid-1980's crystallization ages estimated for Nakhla ranged from approximately 1.3 b.y to 4.5 b.y. Similarly, preferred crystallization ages for the shergottites ranged from 360 m.y., to 1.3 b.y., to 4.5 b.y. In all these models, the 180 m.y. event seen in the shergottites was deemed to be metamorphic. The fit between the Nakhla-like source region and the shergottite hyperplane was a validation both of the 1.3 b.y. igneous age of Nakhla and the 180 m.y. igneous age of the shergottites.

Jones, J. H.↗

EGS Collab Experiment 1: 3D Seismic Velocity Model and Updated Microseismic Catalog Using Transfer-Learning Aided Double-Difference Tomography

This package contains a 3D Seismic velocity model and an updated microseismic catalog associated with a proceedings paper (Chai et al., 2020) published in the 45th Workshop on Geothermal Reservoir Engineering. The 3D_seismic_velocity_model text file contains x (m), y(m), z(m), P-wave velocity (km/s), P-wave velocity quality indicator (1 for well-constrained; 0 for poorly constrained), S-wave velocity (km/s), and S-wave velocity quality indicator (1 for well-constrained; 0 for poorly constrained). The Updated_MEQ_catalog text file contains event origin time, x(m), y(m), z(m), error in x (m), error in y (m), error in z (m), and RMS misfit (millisecond). The 3D_seismic_P-wave_velocity_model animation file shows slices of the 3D P-wave velocity model. The 3D_seismic_S-wave_velocity_model animation file shows slices of the 3D S-wave velocity model. The Interactive_MEQ_locations API file is an interactive visualization of the updated microseismic event locations. The visualization allows users to view the event locations by dragging, rotating, and zooming in. References: Chai, C., Maceira, M., Santos-Villalobos, H. J., Venkatakrishnan, S. V., Schoenball, M., and EGS Collab Team, 2020, Automatic Seismic Phase Picking Using Deep Learning for the EGS Collab Project, in PROCEEDINGS, 45th Workshop on Geothermal Reservoir Engineering, edited, Stanford University, Stanford, California, 45, 1266-1276.

15 GEOTHERMAL ENERGY↗

Optimal carbon storage reservoir management through deep reinforcement learning

Model-based optimization plays a central role in energy system design and management. The complexity and high-dimensionality of many process-level models, especially those used for geosystem energy exploration and utilization, often lead to formidable computational costs when the dimension of decision space is also large. This work adopts elements of recently advanced deep learning techniques to solve a sequential decision-making problem in applied geosystem management. Specifically, a deep reinforcement learning framework was formed for optimal multiperiod planning, in which a deep Q-learning network (DQN) agent was trained to maximize rewards by learning from high-dimensional inputs and from exploitation of its past experiences. To expedite computation, deep multitask learning was used to approximate high-dimensional, multistate transition functions. Both DQN and deep multitask learning are pattern based. As a demonstration, the framework was applied to optimal carbon sequestration reservoir planning using two different types of management strategies: monitoring only and brine extraction. Both strategies are designed to mitigate potential risks due to pressure buildup. Results show that the DQN agent can identify the optimal policies to maximize the reward for given risk and cost constraints. Finally, experiments also show that knowledge the agent gained from interacting with one environment is largely preserved when deploying the same agent in other similar environments.

15 GEOTHERMAL ENERGY↗

Impacts of irregularly-distributed acidified brine flow on geo-chemo-mechanical alteration in an artificial shale fracture under differential stress

The efficacy of geological carbon sequestration is reliant on the integrity of the caprock and its resistance to physical and chemical alteration. Caprocks with high abundance of reactive carbonates like calcite are susceptible to acid-promoted dissolution and can result in structural weakening. This work investigates the effect of acidified brine flow through an artificially fractured, high-carbonate (30 % by XRD) shale under differential compressive stress. Cylindrical samples were cut in half vertically and milled to create an artificial fracture with interlocking asperities and open channels. Samples were sheared with a single applied stress in a custom flow cell housed within an industrial CT scanner. Further, either acidic (pH 4) or reservoir-simulated (pH 9.5) brine was flowed through the artificial fracture for 7–8 days under reservoir pressure and room temperature. Model simulations indicate flow mainly occurred in open channels, with limited flow between overlapping asperities. Analysis of fracture surfaces by optical and scanning electron microscopy show increased surface alteration and roughness after exposure to pH 4 versus pH 9.5 brine indicating mineral dissolution/loss, and this effect is greater in areas that receive the highest brine flows. Similarly, surface analysis by scratch testing shows fracture toughness decreases more after exposure to acidic versus reservoir-simulated brine, with the greatest alteration in areas of highest acidic brine flows. Despite weakening, no shear slip occurred. Overall, the results indicate that acidified brine can result in significant physical and geomechanical alteration of irregular fracture surfaces in shale caprock, with greatest effects in preferential flow regions.

58 GEOSCIENCES↗

Recurrent neural networks for short-term and long-term prediction of geothermal reservoirs

Accurate prediction of geothermal reservoir responses to alternative energy production scenarios is critical for optimizing the development of the underlying resources. While the conventional physics-based models offer a comprehensive prediction tool, data-driven models provide an efficient alternative to build fit-for-purpose predictive models by extracting and using the statistical patterns in the collected data to make predictions. The recurrent neural network (RNN) is a data-driven model that is commonly applied to predict time series sequences. This paper presents a variant of RNN that also utilizes the efficiency of convolutional neural networks (CNN) for the prediction of energy production from geothermal reservoirs. Specifically, a CNN–RNN architecture is developed that takes historical well controls as input (features) and their corresponding production response data as output (labels) to learn an input-output mapping that can predict the future well production responses/performance for any given future well control inputs. The model is paired with a labeling scheme to handle real field disturbances that create data gaps. In addition to the model structure, we introduce a thorough workflow for applying the model, which includes data pre-processing, feature selection, as well as different training strategies for short-term and long-term prediction. Finally, the performance and accuracy of the model are evaluated by applying it to multiple datasets, including a field reservoir model.

15 GEOTHERMAL ENERGY↗

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↗

GEOPHIRES v2.0

GEOPHIRES v2.0 is a geothermal techno-economic simulation tool upgraded from GEOPHIRES v1.2. The name stands for "GEOthermal energy for Production of Heat and electricity (“IR”) Economically Simulated". GEOPHIRES combines reservoir, wellbore, and surface plant technical models with cost correlations and levelized cost models to estimate the capital and operation and maintenance costs, instantaneous and lifetime energy production, and overall levelized cost of energy of a geothermal plant. In addition to electricity generation, direct-use heat applications and combined heat and power or cogeneration can be modeled. Various upgrades have been implemented to develop GEOPHIRES v2.0, including: updating the built-in cost correlations, coupling to the external reservoir simulator TOUGH2, enhancing the built-in wellbore simulator, converting the programming language to Python, and making the code open-source. Users are encouraged to build upon the tool, e.g. by implementing their own correlations or coupling to other simulators.

Beckers, Koenraad↗

ResORR: A Globally Scalable and Satellite Data-Driven Algorithm for River Flow Regulation Due to Reservoir Operations

We propose a globally scalable algorithm, ResORR (Reservoir Operations driven River Regulation), to predict regulated river flow and tested it over the heavily regulated basin of the Cumberland River in the US. ResORR was found able to model regulated river flow due to upstream reservoir operations of the Cumberland River. Over a mountainous basin dominated by high rainfall, ResORR was effective in capturing extreme flooding modified by upstream hydropower dam operations. On average, ResORR improved regulated river flow simulation by more than 50% across all performance metrics when compared to a hydrologic model without a regulation module. ResORR is a timely software algorithm for understanding human regulation of surface water as satellite-estimated reservoir state is expected to improve globally with the recently launched Surface Water and Ocean Topography (SWOT) mission.

River Regulation↗