Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geologic 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 1,225 records · Page 68

Effects of anisotropic shale creep on the stress and permeability evolution of a geological nuclear waste repository

To ensure long-term safety and performance, geological nuclear waste repositories require low-permeability barriers such as bentonite buffers and/or shale host rock. Shale is not only known for its low permeability but also for its trend to undergo time-dependent deformation (i.e., creep), which could heal damage, but the effects of shale creep on the long-term performance of nuclear waste repositories have not been clearly understood. In particular, the anisotropic nature of shale (i.e., bedding) could have a significant effect on its creep behavior, and consequently, on the long-term performance of nuclear waste repositories. In this research, numerical simulations were carried out with the objective of showing the effects of anisotropic shale creep on the stress and permeability evolution of a generic geological nuclear waste repository in shale. The TOUGH-FLAC simulator, a thermo-hydromechanically (THM) coupled numerical code, was used for the simulations. To achieve the objective, comparisons were performed between the results of anisotropic shale creep simulations and those of different simulation cases, namely, no creep (i.e., elastic), isotropic creep, and long-term creep shale cases. Results of the comparisons show that the elastic and isotropic creep shale cases respectively led to the overestimation and underestimation of stress and permeability in the repository, whereas the long-term creep shale case, which accumulated greater creep in later periods than in earlier periods, helped to keep large shear and tensile stresses from developing while maintaining compressive spherical stress, resulting in consistently low permeability levels. These results indicate that performance assessments with elastic and isotropic creep formation models will provide the upper and lower bound estimates of stress and permeability, while more reasonable estimates will be given by an anisotropic creep formation model, and that shale with long-term creep characteristics will be beneficial in many aspects of the safety and performance of nuclear waste repositories.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Divide and conquer: separating the two probabilities in seismic phase picking

There are two fundamental probabilities in the seismic phase picking process—the probability of the existence of a seismic phase (detection probability) and the probability associated with the phase arrival time estimation (timing probability). The nearly ubiquitous approach in developing deep learning phase picking models is to use a kernel, such as a truncated Gaussian, to mask the labelled phase arrival time and train a segmentation model. Once a model is trained, the times of the peaks in the output are taken as phase arrival times (picks), and the height of the peaks are taken as ‘probability’ of the picks. Here, we show that this ‘probability’ represents neither the detection nor the timing probability because this approach forces the output to follow the shape of the kernel. We introduce an approach using two models to estimate these two distinct probabilities. We use a binary classifier with a calibrated confidence to address the detection probability and a multiclass classifier to obtain a probability mass function to address the timing probability. This new approach can make the deep learning-based phase picking process more interpretable and provide options to logically control seismic monitoring workflows.

58 GEOSCIENCES↗

Response of Global SSTs and ENSO to the Atlantic and Pacific Meridional Overturning Circulations

Consequences from a slowdown or collapse of the Atlantic meridional overturning circulation (AMOC) could include modulations to El Niño–Southern Oscillation (ENSO) and development of the Pacific meridional overturning circulation (PMOC). Despite potential ramifications to the global climate, our understanding of the influence of various AMOC and PMOC states on ENSO and global sea surface temperatures (SSTs) remains limited. Five multicentennial, fully coupled model simulations created with the Community Earth System Model were used to explore the influence of AMOC and PMOC on global SSTs and ENSO. We found that the amplitude of annual cycle SSTs across the tropical Pacific decreases and ENSO amplitude increases as a result of an AMOC shutdown, irrespective of PMOC development. However, active deep overturning circulations in both the Atlantic and Pacific basins reduce ENSO amplitude and variance of monthly SSTs globally. Furthermore, the underlying physical reasons for changes to global SSTs and ENSO are also discussed, with the atmospheric and oceanic mechanisms that drive changes to ENSO amplitude differing based on PMOC state. These results suggest that if climate simulations projecting AMOC weakening are realized, compounding climate impacts could occur given the far-reaching ENSO teleconnections to extreme weather and climate events. More broadly, these results provide us with insight into past geologic era climate states, when PMOC was active.

54 ENVIRONMENTAL SCIENCES↗

NNSA-IAEC Science Area V Environmental ISR: Waste Management and Subsurface Science (Final Report FY18-FY20)

In August 2017, the U.S. National Nuclear Security Administration (NNSA) and the Israel Atomic Energy Commission (IAEC) launched a joint waste management program to evaluate the feasibility of siting a radioactive waste repository at intermediate depths in Israel. The bilateral collaboration involves three National Laboratories (LLNL, LANL, and Sandia), and the Nuclear Research Center, Negev, Geological Survey of Israel, and Ben Gurion University in Israel. This research is designed to assist with the evaluation of locations for an intermediate nuclear waste repository in Israel. This final report is a high-level summary of the results from the past three years by the team that worked on the WM-2: Radionuclide colloid-facilitated transport in fractured carbonate rock project for Science Area V, Environmental ISR: Subsurface Science and Waste Management. We have completed all of our tasks and met all of our milestones. MW-2 had three integrated tasks. Task 1 was laboratory experiments investigating colloid transport of radionuclides in natural fractured carbonate rocks from the Avdat formation, Israel. Task 2 was field experiments injecting radionuclide analogues pre-sorbed to clay colloids into fractured carbonate rock, and task 3 was numerical modeling of laboratory and field experimental results to assess the importance of colloids in the migration of relevant radionuclides. This task was incorporated into tasks 1 and 2. A 4 th task was initiated in FY20, to investigate the role of organics in facilitating radionuclide transport under the same conditions explored in task 1 due to the high organic content of the rocks being investigated for the immediate borehole locations. Preliminary results will be summarized here and the work will continue in FY21-23. The overall objective of this research is to evaluate the role of colloids (naturally occurring < 1 micron particles) in facilitating the transport of long-lived radionuclides from a nuclear waste repository situated in fractured carbonate rocks. Currently little data exists on radionuclide transport in carbonate rocks, but this is the main rock type available for siting a repository in Israel. Laboratory experiments were carried out in both Israel and the U.S., field experiments have taken place in Israel, and reactive transport modeling involved LLNL, LANL and Israel partner institutions. LLNL hosted a graduate student, Emily Tran, from Ben Gurion University in FY18 and FY19 and much of the work presented here was part of her PhD research.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Interpretation of hydrothermal conditions, production-injection induced effects, and evidence for enhanced geothermal system-type heat exchange in response to >30 years of production at Roosevelt Hot Springs, Utah, USA

The Roosevelt Hot Springs hydrothermal system is located at the base of the Mineral Mountains in southwestern Utah on the eastern side of the Basin and Range. Hydrothermal activity is related to relatively recent bimodal magmatism, and the system is hosted in coarsely crystalline rock made of Oligocene–Miocene granitoids and Precambrian gneiss. The hydrothermal plume covers ~5 km 2 , with a maximum temperature of 268 °C at ~750 m depth, and a vertically extensive fault-fracture mesh east of the Opal Mound fault controls the upflow of hydrothermal fluids. Power generation (currently 38 MWe gross) began in 1984, and up through 2016, four wells were used for fluid production, and three wells were used for edge-field injection. Chemical analyses of produced fluids show that modern reservoir fluid compositions are similar to but more concentrated than those at the start of production, having near-neutral pH, total dissolved solids of 7000–10,000 mg/kg, and ionic ratios of Cl/HCO 3 ~50–100, Cl/SO4 ~50–100, and Na/K ~4–5. Chemical geothermometers indicate equilibration temperatures that mainly range between 240° and 300 °C. Early production induced a steep drop in pressure (~3.0–3.5 MPa), which was accompanied by a 250–300 m lowering of piezometric levels in wells and development of a shallow steam zone across the system. Hydrothermal fluid compositions evolved continuously in response to production-related steam-loss and injection breakthrough, which is reflected by gradual increases in chloride of up to 35% and stable isotope ratios of up to ~2‰δ 18 O and ~10‰δD. Simple mixing model calculations suggest that there has been a significant amount, ~10–20 MW th , of sustained multi-decadal heat mining and enhanced geothermal system (EGS)–type heat transfer by the injectate as it returns to the production zone. Overall, the two factors that have sustained long-term power production (currently 38 MWe gross) are the increased upflow of deep chloride water and, to a lesser extent, the mining of heat at <1 km depth.

58 GEOSCIENCES↗

Reconstructing atmospheric H 2 over the past century from bi-polar firn air records

Abstract. Historical atmospheric H2 levels were reconstructed using firn air measurements from two sites in Greenland (NEEM and Summit) and two sites in Antarctica (South Pole and Megadunes). A joint reconstruction based on the two Antarctic sites yields H2 levels monotonically increasing from about 330 ppb in 1900 to 550 ppb in the late 1990s, leveling off thereafter. These results are similar to individual reconstructions published previously (Patterson et al., 2020, 2021). Interpretation of the Greenland firn air measurements is complicated by challenges in modeling enrichment induced by pore close-off at these sites. We used observations of neon enrichment at NEEM and Summit to tune the parameterization of enrichment induced by pore close-off in our firn air model. The joint reconstruction from the Greenland data shows H2 levels rising 30 % between 1950 and the late 1980s, reaching a maximum of 530 ppb. After 1990, reconstructed atmospheric H2 levels over Greenland are roughly constant, with a small decline of 3 % over the next 25 years. The reconstruction shows good agreement with the available flask measurements of H2 at high northern latitudes.

Geology↗

Discovering Hidden Geothermal Signatures using Unsupervised Machine Learning

Discovering hidden geothermal resources is a very challenging task. It requires the mining of large datasets, including various diverse data attributes representing subsurface hydrogeological and geothermal conditions. The commonly used Play Fairway Analysis (PFA) typically relies on subject-matter expertise to analyze site or regional data to estimate geothermal conditions and prospectivity. Here, we demonstrate an alternative approach based on machine learning (ML) to process a geothermal dataset of Southwest New Mexico (SWNM). The study region includes low- and medium-temperature hydrothermal systems. However, most of these systems are poorly characterized because of insufficient existing data and limited past explorative studies. This study aims to discover hidden patterns and relationships in the SWNM geothermal dataset to better understand regional hydrothermal conditions. This is achieved by applying an unsupervised machine learning algorithm based on non-negative matrix factorization coupled with customized k-means clustering (NMFk). NMFk can automatically identify (1) hidden (latent) signatures characterizing datasets, (2) the optimal number of these signatures, (3) dominant data attributes associated with each signature, and (4) spatial distribution of the extracted signatures. Here, NMFk is applied to analyze 18 geological, geophysical, hydrogeological, geothermal attributes at 44 locations in SWNM. NMFk successfully finds data patterns and identifies the spatial associations of hydrothermal signatures with the four physiographic provinces in SWNM (Colorado Plateau, Volcanic Field, Basin and Range, and the Rio Grande rift). The algorithm identified up to 5 hydrothermal signatures in the SWNM datasets that differentiate between low- and medium-temperature hydrothermal systems in different provinces. Also, the algorithm identifies two medium-temperature hydrothermal systems in SWNM that require further exploration for geothermal resource development. Based on our analyses, 12 of the attributes are important to identify medium-temperature hydrothermal systems, and the remaining six attributes are critical to characterize low-temperature hydrothermal systems. Based on the obtained results, we identify potential physiographic provinces for further exploration to characterize them as geothermal resources. The resulting NMFk model can be applied to predict geothermal conditions and their uncertainties at new SWNM locations based on limited data from unexplored areas.

58 GEOSCIENCES↗

Physical Insights From the Multidecadal Prediction of North Atlantic Sea Surface Temperature Variability Using Explainable Neural Networks

Abstract North Atlantic sea surface temperatures (NASST), particularly in the subpolar region, are among the most predictable in the world's oceans. However, the relative importance of atmospheric and oceanic controls on their variability at multidecadal timescales remain uncertain. Neural networks (NNs) are trained to examine the relative importance of oceanic and atmospheric predictors in predicting the NASST state in the Community Earth System Model 1 (CESM1). In the presence of external forcings, oceanic predictors outperform atmospheric predictors, persistence, and random chance baselines out to 25‐year leadtimes. Layer‐wise relevance propagation is used to unveil the sources of predictability, and reveal that NNs consistently rely upon the Gulf Stream‐North Atlantic Current region for accurate predictions. Additionally, CESM1‐trained NNs successfully predict the phasing of multidecadal variability in an observational data set, suggesting consistency in physical processes driving NASST variability between CESM1 and observations.

Geology↗

Sensitive Response of Atmospheric Oxidative Capacity to the Uncertainty in the Emissions of Nitric Oxide (NO) From Soils in Amazonia

Abstract Soils are a major source of nitrogen oxides, which in the atmosphere help govern its oxidative capacity. Thus the response of soil nitric oxide (NO) emissions to forcings such as warming or forest loss has a meaningful impact on global atmospheric chemistry. We find that the soil emission rate of NO in Amazonia from a common inventory is biased low by at least an order of magnitude in comparison to tower‐based observations. Accounting for this regional bias decreases the modeled global methane lifetime by 1.4%–2.6%. In comparison, a fully deforested Amazonia, representing a 37% decrease in global emissions of isoprene, decreases methane lifetime by at most 4.6%, highlighting the sensitive response of oxidation rates to changes in emissions of NO compared to those of terpenes. Our results demonstrate that improving our understanding of soil NO emissions will yield a more accurate representation of atmospheric oxidative capacity.

Geology↗

Foundations of Molecular 'Isotomics'

The naturally occurring rare isotopes are versions of common elements, such as hydrogen, carbon and oxygen, that contain a larger than usual number of neutrons in their atomic nuclei and therefore are higher in mass than the common atoms of that element. Isotopes exist for most elements and are found in most natural and synthetic materials, but are uneven in their distribution because chemical and physical processes are isotope-selective (e.g., a chemical reaction may proceed more rapidly for one isotope than for another). For this reason, abundances of isotopes in a material of interest can provide a record, or ‘signature’ of various features of that material’s origin and history. These signatures have been used in the geo, life, chemical and physical sciences in a wide variety of ways over close to 8 decades. However, many such applications struggle to reach unique interpretations of isotopic data because multiple factors combine to control a given sample’s overall isotopic content. That is, the factors controlling isotopic content are too numerous and complex to fully constrain from a simple measurement of a material’s isotope abundances. However, the distribution of isotopes within materials, at molecular scales potentially provides a vastly larger number and diversity of constraints on the chemical and physical processes that comprise a material’s history. The rare isotopes may be concentrated into one atomic position in a molecule relative to another, some proportion of molecules in a sample may contain two or more rare isotopes, and those multiply-isotope-substituted forms of molecules may also have uneven distributions of those isotopes across individual atomic sites. For these reasons, even small, seemingly simple molecules, such as sugars, amino acids or drug compounds, actually exist in a vast number of isotopically unique forms (often millions or more), and each one of those forms is in some sense an independent ‘vote’ on that sample’s history. This project has focused on opening this rich archive of information by enabling the creation of routinely and widely applicable ways of measuring and interpreting isotopic structures of molecules. This work has included the development of core technologies and analytical methods, advancing fundamental understanding of the physical and chemical properties of isotopic versions of molecules, and conducting proof of concept studies of illustrative geochemical, cosmochemical and forensic problems in order to show how these technologies, methods and principles come together to solve problems in new ways. A key to the success of this project was the adaptation of ‘Fourier transform mass spectrometry’ (FTMS) to the task of precisely measuring proportions of the rare, naturally occurring isotopic forms of molecules. FTMS is a highly specialized form of mass spectrometry that traps ions within magnetic or electrostatic cavities and, effectively, ‘listens’ (through registering of subtle electrical signals) to the harmonic signals they make while rapidly orbiting within those cavities. These signals have periods that are a function of their mass and strength (or ‘loudness’) that is proportional to their abundances. Thus, these signals constrain relative amounts of molecules that differ in their mass due to various isotopic substitutions. This technology has been essential to the identification of organic molecules in the life, chemical and environmental sciences for over 4 decades, but generally has lacked the control, stability and precision to meaningfully measure rare isotope forms of molecules. This project’s most fundamental contribution has been to modify FTMS, both in terms of hardware and methods, to enable such measurements. The raw data of molecular isotopic structure is tremendously voluminous and complex, so another important activity of this project has been developing the theoretical and data-science tools needed to interpret the data generated by this new form of isotopic measurement. A particularly challenging part of this task has been predicting molecular isotopic structure, as only through the comparison of measurements with predictions can we make progress on hypothesis driven research questions. We have attacked this this prediction task through a combination of first-principles chemical-physics models of the effects of isotope substitution on molecule properties and data-science models that permit us to generalize that chemical physics to cases that have not yet been studied by detailed chemical physics theory. The proof of concept applications we have pursued over the course of this study include biological reactions of amino acids and other biomolecules, non-biological synthesis of organic molecules in extra-terrestrial settings such as meteorites, petroleum geoscience questions concerning the origin and evolution of natural gas, oil and kerogen compounds, and forensic questions such as the sourcing of chemical weapons. The successes of these applications have laid the groundwork for the next phase of this field’s development, which will include larger scale and more ambitious studies of molecular isotopic structure as a means of diagnosing human diseases, such as cancer, and reconstructing detailed interpretations of the origin and evolution of organic molecules in modern and geological environments.

Cesar, Jaime↗

Stable isotope equilibria in the dihydrogen-water-methane-ethane-propane system. Part 1: Path-integral calculations with CCSD(T) quality potentials

Isotopic compositions of alkanes are typically assumed to be kinetically controlled, but recently is has been proposed that alkanes can isotopically equilibrate for both C and H isotopes during natural gas generation. Evaluation of this requires knowledge of the isotopic equilibrium between alkanes and other common hydrogen and carbon bearing species. Here, in this study, we calculate isotopic equilibria within and between gaseous dihydrogen (H 2 ), water (H 2 O), methane (CH 4 ), ethane (C 2 H 6 ) and propane (C 3 H 8 ), including isotope fractionation among molecules, clumped isotope effects, as well as among sites of propane (i.e., the site-specific isotope effects) from 0°C to 500°C using a path-integral method paired with high-level descriptions of molecular potentials and the diagonal correction to the Born Oppenheimer approximation. While path-integral calculations with high- level CCSD(T) potentials are available for the isotopic equilibria involving methane, the path-integral calculations for ethane and propane have only been performed based on lower-level descriptions of the molecular potentials. We analyze the relative importance of various approximations that are commonly employed when isotopic equilibria are evaluated. We find that clumped isotope effects can be calculated to the same accuracy using computationally inexpensive combination of the Bigeleisen-Mayer-Urey model with the molecular potential from density functional theory. In contrast, fractionation and site preferences of both deuterium and carbon-13 benefit from the use of the higher level CCSD(T) potentials and accounting for anharmonic effects. Additionally, for fractionation and site preference of deuterium corrections to Born-Oppenheimer approximation can also be important.

03 NATURAL GAS↗

Thermal loading effects on chalk hydromechanical behavior for nuclear waste disposal

Safe disposal of heat-generating nuclear waste depends on host rock stability under thermal, hydrological, and mechanical stresses. This study investigates the effect of thermal loading on mechanical behavior of the shallowly buried Ghareb formation chalk through triaxial and hydrostatic constant strain rate and creep tests at temperatures up to 100 ˚C and effective pressures up to 20.7 MPa. Experimental results show that thermal loading reduces the elastic moduli of chalk by 50–75%, and a transition occurs above 60 ˚C where creep rates increase rapidly. Water saturation nearly doubles the thermally induced strain compared to dry conditions and strongly decreases material rigidity. Thermal loading also leads to significant pore pressure increases under undrained conditions and leads to reductions in the apparent permeability during drained conditions. Laboratory experimental data were used to parameterize and develop a preliminary constitutive model for predicting future deformation during repository operations in the Ghareb. The strongly coupled effects – mechanical weakening, fluid pressure fluctuations, and permeability modification – demonstrate that elevated repository temperatures will have a pronounced effect on the near field Ghareb behavior during waste disposal operations. The findings indicate that the coupled interactions must be considered in predictive models and repository design to ensure long-term nuclear waste isolation and safety.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Subsurface microbial communities as a tool for characterizing regional-scale groundwater flow

Subsurface microbial community distribution patterns are influenced by biogeochemical and groundwater fluxes and may inform hydraulic connections along groundwater-flow paths. This study examined the regional-scale microbial community of the Death Valley Regional Flow System and evaluated whether subsurface communities can be used to identify groundwater-flow paths between recharge and discharge areas. Samples were collected from 36 sites in three groundwater basins: Pahute Mesa–Oasis Valley (PMOV), Ash Meadows (AM), and Alkali Flat–Furnace Creek Ranch (AFFCR). Microbial diversity within and between communities varied by location, and communities were separated into two overall groups that affiliated with the AM and PMOV/AFFCR basins. Network analysis revealed patterns between clusters of common microbes that represented groundwaters with similar geochemical conditions and largely corroborated hydraulic connections between recharge and discharge areas. Null model analyses identified deterministic and stochastic ecological processes contributing to microbial community assemblages. Most communities were more different than expected and governed by dispersal limitation, geochemical differences, or undominating processes. However, certain communities from sites located within or near the Nevada National Security Site were more similar than expected and dominated by homogeneous dispersal or selection. Overall, the (dis)similarities between the microbial communities of DVRFS recharge and discharge areas supported previously documented hydraulic connections between: (1) Spring Mountains and Ash Meadows; (2) Frenchman and Yucca Flat and Amargosa Desert; and (3) Amargosa Desert and Death Valley. However, only a portion of the flow path between Pahute Mesa and Oasis Valley could be supported by microbial community analyses, likely due to well-associated artifacts in samples from the two Oasis Valley sites. This study demonstrates the utility of combining microbial data with hydrologic, geologic, and water-chemistry information to comprehensively characterize groundwater systems, highlighting both strengths and limitations of this approach.

54 ENVIRONMENTAL SCIENCES↗

Hyperspectral Reflectance for Measuring Canopy‐Level Nutrients and Photosynthesis in a Salt Marsh

Abstract Salt marsh ecosystems are underrepresented in process‐based models due to their unique location across the terrestrial–aquatic interface. Particularly, the role of leaf nutrients on canopy photosynthesis (F A ) remains unclear, despite their relevance for regulating vegetation growth. We combined multiyear information of canopy‐level nutrients and eddy covariance measurements with canopy surface hyperspectral remote sensing (CSHRS) to quantify the spatial and temporal variability of F A in a temperate salt marsh. We found that F A showed a positive relationship with canopy‐level N at the ecosystem scale and for areas dominated by Spartina cynosuroides , but not for areas dominated by short S. alterniflora . F A showed a positive relationship with canopy‐level P, K, and Na, but a negative relationship with Fe, for areas associated with S. cynosuroides , S. alterniflora , and at the ecosystem scale. We used partial least squares regression (PLSR) with CSHRS and found statistically significant data–model agreements to predict canopy‐level nutrients and F A . The red‐edge electromagnetic region and ∼770 nm showed the highest contribution of variance in PLSR models for canopy‐level nutrients and F A , but we propose that underlying sediment biogeochemistry can complicate interpretation of reflectance measurements. Our findings highlight the relevance of spatial variability in salt marshes vegetation and the promising application of CSHRS for linking information of canopy‐level nutrients with F A . We call for further development of canopy surface hyperspectral methods and analyses across salt marshes to improve our understanding of how these ecosystems will respond to global environmental change.

Environmental Sciences & Ecology↗

Interpretation of Hydraulic Fractures Based on Microseismic Response in the Marcellus Shale, Monongalia County, West Virginia, USA: Implications for Shale Gas Production

Summary Hydraulic fracturing is critical for extracting shale gas in the subsurface. The treatment technique of multistage hydraulic fracturing is widely used to maximize production. In multistage hydraulic fracturing, not all pumping stages make the same contribution for production, although the designed stimulation process is almost same in every pumping stage. In this study, we characterize the microseismic responses of hydraulic fracturing. Significant variations of microseismic characteristics are observed among the different pumping stages. Pre-existing natural fractures are examined along the horizontal well in selected stages. Combined with the Mohr’s circle analysis, results of the fracture study show that induced hydraulic fractures can be captured by pre-existing natural fractures. Induced hydraulic fractures are simulated by the unconventional fracture model (UFM), and the result reveals that the stress-shadow effect diverts the direction of hydraulic fractures. The diverted hydraulic fracture affects the development of the hydraulic fracture network, which has influence on production. The treatment of two-step pumping is investigated by application of hydraulic diffusivity. Hydraulic fracturing performance and production could be optimized by the treatment of two-step pumping in a single stage. The second pumping creates new fractures and fills the fractures with additional proppants to maintain production for a long duration.

Energy & Fuels↗

Optimization of Water-Alternating-CO2 Injection Field Operations Using a Machine-Learning-Assisted Workflow

Summary This paper will present a robust workflow to address multiobjective optimization (MOO) of carbon dioxide (CO2)-enhanced oil recovery (EOR)-sequestration projects with a large number of operational control parameters. Farnsworth unit (FWU) field, a mature oil reservoir undergoing CO2 alternating water injection (CO2-WAG) EOR, will be used as a field case to validate the proposed optimization protocol. The expected outcome of this work would be a repository of Pareto-optimal solutions of multiple objective functions, including oil recovery, carbon storage volume, and project economics. FWU’s numerical model is used to demonstrate the proposed optimization workflow. Because using MOO requires computationally intensive procedures, machine-learning-based proxies are introduced to substitute for the high-fidelity model, thus reducing the total computation overhead. The vector machine regression combined with the Gaussian kernel (Gaussian-SVR) is used to construct proxies. An iterative self-adjusting process prepares the training knowledge base to develop robust proxies and minimizes computational time. The proxies’ hyperparameters will be optimally designed using Bayesian optimization to achieve better generalization performance. Trained proxies will be coupled with multiobjective particle swarm Optimization (MOPSO) protocol to construct the Pareto-front solution repository. The outcomes of this workflow will be a repository containing Pareto-optimal solutions of multiple objectives considered in the CO2-WAG project. The proposed optimization workflow will be compared with another established methodology using a multilayer neural network (MLNN) to validate its feasibility in handling MOO with a large number of parameters to control. Optimization parameters used include operational variables that might be used to control the CO2-WAG process, such as the duration of the water/gas injection period, producer bottomhole pressure (BHP) control, and water injection rate of each well included in the numerical model. It is proved that the workflow coupling Gaussian-SVR proxies and the iterative self-adjusting protocol is more computationally efficient. The MOO process is made more rapid by squeezing the size of the required training knowledge base while maintaining the high accuracy of the optimized results. The outcomes of the optimization study show promising results in successfully establishing the solution repository considering multiple objective functions. Results are also verified by validating the Pareto fronts with simulation results using obtained optimized control parameters. The outcome from this work could provide field operators an opportunity to design a CO2-WAG project using as many inputs as possible from the reservoir models. The proposed work introduces a novel concept that couples Gaussian-SVR proxies with a self-adjusting protocol to increase the computational efficiency of the proposed workflow and to guarantee the high accuracy of the obtained optimized results. More importantly, the workflow can optimize a large number of control parameters used in a complex CO2-WAG process, which greatly extends its utility in solving large-scale MOO problems in various projects with similar desired outcomes.

Energy & Fuels↗

Hybrid Natural Gas Geothermal Combined Cycle Power Plant Analysis

Low temperature geothermal resources, including those associated with oil and gas production, are an underutilized source of low carbon energy. The present work investigates coupling of low-temperature geothermal resources with concentrating solar and/or natural gas energy sources to increase the number of locations at which power generation from low temperature geothermal resources would be technically and economically viable. Stand-alone and hybrid geothermal power cycles are simulated using SimTech IPSEpro process modeling software. Design point strategies for a hybrid power cycle that may operate with either a single heat source or two simultaneous heat input sources are considered. Additionally, off-design power plant operation is investigated to consider the impacts of the heat source availability and ambient temperature variations. The off-design analysis utilizes a modeling tool that predicts power plant performance at each time step as a function of topping cycle heat input (from gas or solar), bottoming cycle heat input (from geothermal), and ambient temperature. Addition of a steam topping cycle to an organic Rankine cycle geothermal power plant provides opportunities to increase the efficiency and power output relative to a stand-alone geothermal power plant. Additionally, use of the waste heat from gas turbine power generation in a geothermal bottoming cycle provides opportunities to increase the amount of power generation associated with each unit of carbon dioxide emitted. This paper will describe the hybrid plant configuration evaluated, discuss the predicted power cycle performance, and compare with stand-alone natural gas and geothermal power generation cases. The power cycle investigated is expected to be applicable for use with conventional hydrothermal resources as well as with geologic thermal energy storage applications and/or enhanced geothermal systems. The steam topping cycle could use a concentrated solar heat source for fully renewable hybrid plant configuration. A plant initially constructed with a natural gas topping cycle heat source could be converted to a solar heat source part way though the power cycle operational life to achieve life cycle carbon emission reductions.

concentrating solar power↗

Ambient field seismology in critical zone hydrological sciences

Passive ambient noise monitoring is an emerging tool in environmental seismology, leveraging the ambient seismic field to assess temporal variations in shallow subsurface properties. This review focuses on the potential and challenges of using scattered coda waves from noise correlation functions to monitor critical zone dynamics. The sensitivity of seismic velocities to various environmental factors, including precipitation, snowmelt, atmospheric pressure, and groundwater fluctuations, underscores the method’s versatility. While coda waves excel in detecting subtle changes due to their scattered nature, ballistic waves provide higher spatial resolution, albeit with challenges in source stability. Advances in seismic sensing, including distributed acoustic sensing and low-cost geophone networks, have enabled high-resolution monitoring of hydrological processes, subsurface deformation, and seismic hazards. Integrating seismic data with hydrological models provides insights into water storage, pore pressure changes, and soil moisture dynamics. However, limitations in spatial resolution, calibration with ground truth data, and coupled effects between environmental factors remain key challenges. This review emphasizes the importance of interdisciplinary approaches in refining methodologies, enhancing sensor deployments, and addressing data gaps. Passive seismic monitoring offers opportunities to understand critical zone processes and their broader impacts on seismic hazards and environmental sustainability.

58 GEOSCIENCES↗