Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reservoir Computing”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

UT-GOM2-2 Drilling Fluid (Technical Note)

It is proposed that a salt-saturated, water-based mud might improve borehole stability for UT-GOM2-2 relative to a fresh-water-based mud. However, the primary objective for UT-GOM2-2 is to safely and successfully acquire uncompromised hydrate-bearing sediment cores. This can only be done if conditions stay within the hydrate stability zone, and well away from the hydrate stability boundary. The addition of salt to the drilling mud, shifts the hydrate stability boundary closer to estimated conditions. This shift, by some estimates of pressure and temperature, where conditions now fall outside of the hydrate stability zone, may result in borehole enlargement and the release of free gas into the borehole. This shift by all estimates, even when coring conditions stay inside the hydrate stability zone, shrinks the window between the estimated conditions and the hydrate stability boundary. This shift is likely to compromise the hydrate-bearing cores while they are being recovered from the bottom-hole to the rig floor. In this report we present two examples: 1) Assuming temperature and pressure from in-situ estimates: Hydrate within the target reservoirs (Orange and Blue) at Walker Ridge Block 313 (WR 313) will be outside the hydrate stability zone in the presence of a drilling mud with 9.5 wt.% salinity (the salinity of the proposed 10.5 ppg salt-based mud) at in situ pressure and temperature. Thus, the hydrate will be unstable. A 10.5 ppg salt-based mud may result in dissociation of the hydrate into its components: water and gas. 10.5 ppg salt-based mud may enlarge the borehole, release free gas into the borehole, and compromise the cores. 2) Assuming temperatures equal to the measured LWD borehole temperatures at this location: Hydrate within the target reservoirs will be stable with a salinity of 9.5 wt.% (the salinity of the proposed 10.5 ppg salt-based mud), but possibly unstable (just at the methane hydrate stability boundary) for a salinity of 17.2 wt.% (the salinity of the proposed 13.5 ppg salt-based kill mud). if the borehole has the same temperature as recorded during previous LWD drilling at this location. A 13.5 ppg salt-based kill mud may destabilize the borehole. More importantly, a 10.5 ppg salt-based mud does not provide a large enough window between the estimated conditions and the hydrate stability boundary. Core temperatures and pressures fluctuate from the bottom-hole conditions as the core is brought up from the bottom-hole to the rig floor. These fluctuations are likely to cause the hydrate in the cores to touch or cross the stability boundary. Therefore, cores captured using a 10.5 ppg salt-based mud are likely to be compromised during core recovery to the rig floor.

03 NATURAL GAS↗
DOE OSTI2026

A robust framework for frictional fault contact in geological formations using a stabilized augmented Lagrangian approach

Numerical simulations are essential to evaluate the performance and safety of engineered subsurface systems such as geological carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately modeling the frictional contact behavior along fault surfaces. This problem involves inequality constraints that arise from the physics of frictional slip, requiring specialized numerical methods to handle the resulting highly nonlinear and path-dependent behavior. Here, in this work, we address this challenge using an Augmented Lagrangian Method (ALM) implemented via the Uzawa algorithm. The formulation employs mixed finite element spaces, combining low-order piecewise linear displacements within the 3D domain cells with piecewise constant tractions defined on the fault surfaces. Furthermore, to ensure stability and satisfy the inf-sup condition, the discrete displacement space is enriched with face bubble functions on both sides of the contact interfaces. This approach offers several advantages over other stabilization techniques that rely on additional terms, and it integrates naturally in the Uzawa framework.

58 GEOSCIENCES↗
DOE OSTI2022

A robust deep learning workflow to predict multiphase flow behavior during geological C O 2 sequestration injection and Post-Injection periods

Simulation of multiphase flow in porous media is essential to manage the geologic CO 2 sequestration (GCS) process, and physics-based simulation approaches usually take prohibitively high computational cost due to the nonlinearity of the coupled physics. This paper contributes to the development and evaluation of a deep learning workflow that accurately and efficiently predicts the temporal-spatial evolution of pressure and CO 2 plumes during injection and post-injection periods of GCS operations. Based on a Fourier Neural Operator, the deep learning workflow takes input variables or features including rock properties, well operational controls and time steps, and predicts the state variables of pressure and CO 2 saturation. To further improve the predictive fidelity, separate deep learning models are trained for CO 2 injection and post-injection periods due to the difference in primary driving force of fluid flow and transport during these two phases. We also explore different combinations of features to predict the state variables. We use a realistic example of CO 2 injection and storage in a 3D heterogeneous saline aquifer, and apply the deep learning workflow that is trained from physics-based simulation data and emulate the physics process. Through this numerical experiment, we demonstrate that using two separate deep learning models to distinguish post-injection from injection period generates the most accurate prediction of pressure, and a single deep learning model of the whole GCS process including the cumulative injection volume of CO 2 as a deep learning feature, leads to the most accurate prediction of CO 2 saturation. For the post-injection period, it is key to use cumulative CO 2 injection volume to inform the deep learning models about the total carbon storage when predicting either pressure or saturation. The deep learning workflow not only provides high predictive fidelity across temporal and spatial scales, but also offers a speedup of 250 times compared to full physics reservoir simulation, and thus will be a significant predictive tool for engineers to manage the long-term process of GCS.

58 GEOSCIENCES↗
DOE OSTI2021

A comprehensive analysis of transient pressure and rate data from CO 2 storage projects in a depleted pinnacle reef oil field complex, Michigan, USA

Pressure and rate data are commonly recorded as part of a basic monitoring program in CCS projects. This paper discusses the application of multiple analytical techniques to interpret pressure and rate transient data from CO 2 injection and storage operations. The techniques of interest, i.e., injection-falloff analysis, injectivity/productivity index analysis and pressure pulse arrival time analysis, are commonly used in the oil and gas industry to assess reservoir properties, but not well known in the CCS literature (especially the last two). Injection-falloff analysis involves log-log pressure derivative plotting for the falloff data and history-matching of the entire injection-falloff sequence to determine permeability. In the injectivity/productivity index analysis, rate-normalized pressure buildup is plotted against material balance time or ratio of cumulative injection to injection rate to determine the injectivity index (ratio of injection rate to stabilized pressure buildup) which can be related to the permeability-thickness product. The arrival time analysis identifies the arrival of a pressure disturbance (~0.1 psi change from ambient) to determine the hydraulic diffusivity from which permeability can be estimated. The applicability of these techniques is demonstrated via illustrative examples from multiple wells in different pinnacle carbonate reefs undergoing CO 2 -EOR in Northern Michigan. The paper ends with a discussion of the relative merits of each interpretive technique, as well as recommendations that could be useful for other field projects.

42 ENGINEERING↗
DOE OSTI2018

Quantitative Characterization of Impacts of Coupled Geomechanics and Flow on Safe and Permanent Geological Storage of CO 2 in Fractured Aquifers

This is the final scientific report for the award DE-FE0023305, entitled “Quantitative Characterization of Impacts of Coupled Geomechanics and Flow on Safe and Permanent Geological Storage of CO 2 in Fractured Aquifers.” The work has been divided into six tasks. In Task 2, we characterized rock properties, which are important when developing a quantitative approach for understanding and predicting geomechanical effects on large-scale CO 2 injection and long-term storage in the subsurface. Rock properties of interest for this characterization include porosity, permeability, elastic constants, strength, and heat capacity. We measured rock properties for three different rock types: concrete, sandstone and shale. These properties were obtained from the acoustic measurement, permeability and porosity measurement, Brazilian test, the uniaxial compressive test, and heat capacity measurement. Then, using the Brazilian technique to fracture samples, we measured their permeability under brine injection for varying effective stresses, followed by similar measurements of permeability under sc-CO 2 injection. Permeability to brine and sc-CO 2 decreased as the effective stress increased. The apparent permeability to sc-CO 2 was an order of magnitude less than that for brine, a consequence of CO 2 being the non-wetting fluid. In Task 3, we developed understandings and correlations for CO 2 injection pressure induced fracturing. We initially fractured four eight inch cubes of concrete with brine to establish a standard test procedure. Then, additional concrete samples were fractured using brine. Confining stresses were exerted on the samples and brine was injected at a constant rate into a borehole that was previously drilled into the concrete cube. The borehole pressure was measured and fracture initiation was identified as a peak in the pressure profile. That was followed by fracture propagation, which was identified as a plateau in the pressure profile, and finally the fracture reached the sample surface. The sample surfaces were photographed before and after fracturing, which enabled determination of where the fracture reached the sample surface. Afterwards, aqueous food dye solution was injected into the fractured concrete samples to color the fracture planes. The borehole was filled with food dye solution in advance, a gas pulse was transmitted to the borehole, and the borehole was then shut in to allow the food dye to be driven into existing fractures. Then, the concrete samples were broken down by high-pressure gas injection to reveal the geometry and morphology of the hydraulic fractures. In addition, acoustic measurements were conducted on multiple locations on each face of the concrete samples before and after injection for comparison. This comparison was used to confirm changes in the internal structure of the concrete sample, namely the formation of fractures. There were six concrete samples fractured with brine. Afterwards, concrete samples were fracture using sc-CO 2 . The main change in experimental procedure was the inclusion of temperature control, since in the field, we mostly encounter temperature and pressure conditions above the supercritical point of CO 2 . The concrete samples were pre-heated in an oven to elevate their temperature, and thermal tape was used to heat the CO 2 to a temperature above its supercritical temperature. Two concrete samples were fractured using liquid CO 2 , four concrete samples were fractured using gaseous CO 2 , and twenty two concrete samples were fractured using sc-CO 2 . For the concrete samples fractured using sc-CO 2 , three had pre-existing hydraulic fractures on the surface near the wellbore, and two were water saturated (as opposed to containing only air). Two concrete samples were composite samples, with a high permeability ball in the center, in order to represent CO 2 injection into a high permeability zone surrounded by a low permeability sealing formation. Later, we fractured five shale samples, obtained from the Niobrara shale outcrop, with sc-CO 2 . The injection rate was constant, confining stresses were applied, the shale samples were pre-heated in an oven to the desired temperature, and the sample was then fractured. An analysis of the failure type was done for many of the concrete samples, based on the predicted break down pressure for tensile and shear failure, and it was found that most samples fractured due to shear failure. Also, the presence of induced fractures has little to no effect on the fracture initiation pressure, but significantly changes the fracture morphology. The CO 2 fracturing in water saturated samples behaved much similar to brine induced fracturing. There was a significant difference of breakdown pressure for injection of CO 2 and brine, where the breakdown pressure for CO 2 induced fracturing is generally around the minimum horizontal stress and that for brine is much higher. Finally, the fracture orientation for most samples with large confining stress differential was perpendicular to the minimum horizontal stress. When the stress differential, especially the difference between the two directions with the smallest stresses was small, the fracture orientation was determined not only by the minimum stress direction, but also the magnitude of the confining stress, the tensile strength of the rock, and the breakdown pressure. In Task 4, we modified our coupled flow-geomechanical models to model fracture growth and propagation in storage formations and caprocks. These flow-geomechanical models are TOUGH2-CSM and TOUGH2-FLAC. The TOUGH2-CSM fluid and heat flow formulation is based on the TOUGH2 formulation of mass and energy conservation equations that govern fluid and heat flow in general multiphase, multicomponent, multi-porosity systems. The TOUGH2-CSM simplified geomechanical formulation is based on the linear theory of elasticity applied to multi-porosity non-isothermal (thermo-multi-poroelastic) media. We previously derived, from the fundamental relations of the linear theory of elasticity, an equation relating mean stress, pore pressures, temperatures, and body force (the Mean Stress equation) that we added to the TOUGH2-CSM fluid and heat flow equations. We extended the TOUGH2-CSM simplified geomechanical formulation to calculate the entire stress tensor by deriving equations for stress tensor components from derivatives of the Cartesian thermo-multi-poroelastic Navier equation components. In addition, we derived equations for stress tensor components in rz-coordinates. We solved these geomechanical equations using the integral finite difference method. This method utilizes momentum fluxes obtained by “factoring” a divergence operator out of each geomechanical equation. In Task 5, we did a literature survey to determine suitable correlations for fracture initiation, growth, and propagation. After surveying numerous papers, we modified TOUGH2-CSM to model stress dependent fracture initiation and growth using the Mohr-Coulomb criterion for shear failure in faults and caprock and the condition of negative minimum effective stress for tensile failure. Tensile failure results in a fracture whose width depends on the difference between fracture pressure and minimum stress and whose propagation is based on the difference between the stress intensity factor at the fracture tip and the rock toughness. The TOUGH-FLAC model is based on the FLAC3D ubiquitous joint model, extended for the modeling of permeability changes induced by changes in effective normal stress across fractures as well permeability changes caused by shearing along fractures. In Task 6, we validated our coupled flow-geomechanical models using analytical solutions and problems from the literature. We validated the TOUGH2-CSM modifications using analytical solutions to the displacement from a uniform load on a semi-infinite elastic medium and the two-dimensional Mandel-Cryer effect. Those were followed by a sample problem to test our stress tensor calculations, which was for injection into a single-phase reservoir with constant properties. Comparisons of our simulator to published results were done on for the depletion of a single-phase reservoir with stress dependent porosity, two examples from the In Salah gas project, an axisymmetric baserock–reservoir–caprock system with a normal faulting stress regime that had thermally induced shear stresses, simulations of CO 2 leakage through caprock in a two-dimensional reservoir in Cartesian coordinates, and predictions of caprock failure from CO 2 injection into an axisymmetric reservoir. Finally, we simulated CO 2 injection pressure-induced fracturing from one of our laboratory studies, Sample 39, and obtained a fracture profile. We also validated the TOUGH2-FLAC model for fluid driven fracture growth against solutions based on the KGD model and a case with an inclined fracture that was loaded from the boundaries to achieve wing-crack propagation. The strain softening tensile behavior and softening of modulus considering a simple damage approach was verified by a simulation representing hydraulic fracturing stress measurement around a vertical well. Finally, the model was validated against deep fracture zone opening and surface uplift at In Salah with inverses analysis used to improve the match of simulation with field data. In Task 7, we developed a scheme based on inverse modeling that can be used to identify caprock leakage. Injection of fluid into a reservoir results in a time varying pressure profile that depends on the properties of the rock and the fluid. The presence of significant caprock leakage would effect this pressure profile and we used inverse modeling, namely the Levenberg-Marquardt method, to determine the leakage location from the pressure profile. A two-dimensional and a three-dimensional example based on published data were presented. The published simulation was run and its results were considered as “field” data. In the two-dimensional example, a simulation of CO 2 injection into a two-dimensional aquifer-caprock system, the caprock had a vertical fault through which leakage could occur. The location of this leakage was obtained from our inverse modeling scheme. In the three-dimensional example, the simulation domain had four geological layers, one of which was the injection zone that contained a horizontal injection well. We simulated fractures as high permeability gridblocks, introduced a fracture that spanned the caprock depth, and ran a simulation to obtain “field” data. Then, using our inverse modeling scheme, we obtained the location of that fracture.

58 GEOSCIENCES↗
DOE OSTI2018

Geocellular Model of Mt. Simon Sandstone for University of Illinois at Urbana-Champaign DDU feasibility study

The geocellular model of the Mt. Simon Sandstone was constructed for the University of Illinois at Urbana-Champaign DDU feasibility study. Starting with the initial area of review (18.0 km by 18.1 km [11.2 miles by 11.3 miles]) the boundaries of the model were trimmed down to 9.7 km by 9.7 km (6 miles by 6 miles) to ensure that the model enclosed a large enough volume so that the cones of depression of both the production and injection wells would not interact with each other, while at the same time minimizing the number of cells to model to reduce computational time. The grid-cell size was set to 61.0 m by 61.0 m (200 feet by 200 feet) for 160 nodes in the X and Y directions. Within the model, 67 layers are represented that are parameterized with their sediment/rock properties and petrophysical data. The top surface of the Mt. Simon Sandstone was provided by geologists working on the project, and the average thickness of the formation was taken from the geologic prospectus they provided. An average thickness of 762 m (2500 feet) was used for the Mt. Simon Sandstone, resulting in 60 layers for the model. Petrophysical data was taken from available rotary sidewall core data (Morrow et al., 2017). As geothermal properties (thermal conductivity, specific heat capacity) are closely related to mineralogy, specifically the percentage of quartz, available mineralogical data was assembled and used with published data of geothermal values to determine these properties (Waples and Waples, 2004; Robertson, 1988). The Mt. Simon Sandstone was divided into three separate units (lower, middle, upper) according to similar geothermal and petrophysical properties, and distributed according to available geophysical log data and prevailing interpretations of the depositional/diagenetic history (Freiburg et al. 2016). Petrophysical and geothermal properties were distributed through geostatistical means according to the associated distributions for each lithofacies. The formation temperature was calculated, based on data from continuous temperature geophysical log from a deep well drilled into the Precambrian basement at the nearby Illinois Basin Decatur Project (IBDP) where CO2 is currently being sequestered (Schlumberger, 2012). Salinity values used in the model were taken from regional studies of brine chemistry in the Mt. Simon Sandstone, including for the IBDP (e.g., Panno et al. 2018). After being reviewed by the project's geologists, the model was then passed onto the geological engineers to begin simulations of the geothermal reservoir and wellbores.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Science-informed Machine Learning to Increase Recovery Efficiency in Unconventional Reservoirs

Production of hydrocarbons from fractured, unconventional reservoirs is inherently inefficient. But machine learning offers a pathway both to increasing recovery efficiency at a site and to improving forecasts of production, thereby improving the economics of operations in unconventional reservoirs. Los Alamos—in partnership with DOE, NETL, and WVU—has been developing a science-informed workflow and platform for optimizing pressure-drawdown at a site, which will allow an operator to make reservoir-management decisions that optimize recovery in consideration of future production. This work relies on a hybridization of physics-based prediction and machine learning, whereby accurate synthetic data (in combination with available site data) can enable the application of machine learning methods for rapid forecasting and optimization. The physics-based prediction is built upon experimental and theoretical work to determine transport characteristics in shale at various scales, with an emphasis on materials from MSEEL-I; this fundamental shale R&D was conducted in partnership with DOE, NETL, and several other national labs. This work has resulted from a coordinated leveraging of developments across several projects within DOE FE30, along with internal investments from Los Alamos via LDRD. The development has utilized data from the MSEEL–I site for calibration and demonstration; however, the workflow and platform are readily extendable to operations at other sites, plays, and basins. This machine-learning method can aid operators to improve both recovery efficiency and competitiveness; to this end, future work would quantify processes for other plays/basins and integrate production details with economics.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2018

Determinants of the efficacy of HIV latency-reversing agents and implications for drug and treatment design

HIV eradication studies have focused on developing latency-reversing agents (LRAs). However, it is not understood how the rate of latent reservoir reduction is affected by different steps in the process of latency reversal. Furthermore, as current LRAs are host-directed, LRA treatment is likely to be intermittent to avoid host toxicities. Few careful studies of the serial effects of pulsatile LRA treatment have yet been done. This lack of clarity makes it difficult to evaluate the efficacy of candidate LRAs or predict long-term treatment outcomes. Here, we constructed a mathematical model that describes the dynamics of latently infected cells under LRA treatment. Model analysis showed that, in addition to increasing the immune recognition and clearance of infected cells, the duration of HIV antigen expression (i.e., the period of vulnerability) plays an important role in determining the efficacy of LRAs, especially if effective clearance is achieved. Patients may benefit from pulsatile LRA exposures compared with continuous LRA exposures if the period of vulnerability is long and the clearance rate is high, both in the presence and absence of an LRA. Overall, the model framework serves as a useful tool to evaluate the efficacy and the rational design of LRAs and combination strategies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

MAL33 drives natural variation in maltose metabolism in Saccharomyces eubayanus

Maltose is one of the most abundant sugars in brewer’s wort, and its efficient utilization is critical for successful fermentation. However, maltose consumption varies naturally among Saccharomyces eubayanus strains isolated from different host trees, such as Quercus and Nothofagus. To identify the genetic determinants underlying these phenotypic differences, we performed bulk segregant analysis (BSA) and quantitative trait loci (QTL) mapping using an F 2 offspring derived from QC18 (Quercus-associated) and CL467.1 (Nothofagus-associated) strains. QTL mapping identified two significant genomic regions on subtelomeric loci of chromosomes V-R and XVI-L, each containing complete MAL loci composed of MAL32 (encoding maltase), MAL31 (transporter), and MAL33 (transcriptional activator) genes. Comparative polymorphism analyses identified mutations in MAL32 and MAL33 of QC18, including frameshift mutations resulting in premature stop codons. Functional validation demonstrated that the heterologous expression of MAL33 ChrV from CL467.1 fully restored maltose utilization in QC18, indicating the functional presence of MAL33 cis-regulatory sequences and MAL32 and MAL31 genes in QC18. While structural protein predictions identified truncation and impaired functionality in the maltose-responsive activation domain of Mal33p from QC18, overexpression of QC18’s own MAL33 ChrV allele also improved maltose metabolism, suggesting dosage-dependent transcriptional limitations rather than complete functional loss. These results indicate that allelic variations in the maltose-responsive activation domain of Mal33p result in differences in maltose consumption between strains. Here, we hypothesized that reduced maltose metabolism in QC18 is an adaptive response to the distinct sugar composition in Quercus robur bark, contrasting with the starch-rich environment of Nothofagus pumilio. These findings highlight subtelomeric MAL gene diversity as a reservoir of genetic variation, representing a key evolutionary mechanism that influences maltose adaptation among natural Saccharomyces isolates.

evolutionary plasticity↗
DOE OSTI2017

Numerical Analysis of Experiments on Thermally Induced Dissociation of Methane Hydrates in Porous Media

Numerical simulation is essential for the prediction and evaluation of hydrocarbon reservoir performance. Numerical simulators developed for the description of the behavior of hydrates under production and the corresponding flow of fluids and heat accounting for all known processes are powerful, but they need validation through comparison to field or experimental data in order to instill confidence in their predictions. In this study, we analyze by means of numerical simulation the results of an experiment of methane hydrate dissociation by thermal stimulation in unconsolidated porous media heated through the vessel walls. The physics captured by the model include multicomponent heat and mass transfer, multiphase flow through porous media, and the phase behavior of the CH 4 + H 2 O system involved in methane hydrate formation and dissociation. The set of governing equations consists of the mass and energy conservation equations coupled with constitutive relationships, i.e., the dissolution of gas in H 2 O, relative permeability and capillary pressure models, composite thermal conductivity models, and methane hydrate phase equilibria. The model geometry describes accurately the hydrate reactor used in a recent experimental study investigating methane hydrate dissociation behavior [Chong et. al. Appl. Energy 2016, 177, 409-421]. The cumulative gas production is estimated and validated against three tests of experimental data involving different boundary temperatures, showing a good agreement between observations and numerical predictions. The predicted evolution of the spatial distributions of different phases over time shows that hydrate dissociation progresses inward from the reactor boundary to the center, methane gas accumulates to the top of the reactor because of buoyancy, and water migrates down to the bottom of the reactor because of gravity. A sharp hydrate dissociation front is predicted, and the estimated location of hydrate dissociation front suggests a linear relationship with the square root of time. A sensitivity analysis on the thermal conductivity of sand under fully saturated conditions is conducted to elucidate its effect on the gas production behavior. In addition, the energy efficiency ratio computed from the simulation of this boundary-wall heating technique varies from 14.0 to 16.2. Deviations between observations and predictions of the evolution of the temperature profile are attributed to initial heterogeneous distribution of the hydrate phase in the hydrate reactor.

42 ENGINEERING↗
DOE OSTI2025

Pseudo-equilibrium theory for extrinsic doping control of the topological semimetal Cd 3 As 2

The standard approach for predicting defect equilibria from first principles assumes that the solid-state system is initially in a thermodynamic equilibrium with the external atomic reservoirs. This “growth step” is then often followed by a temperature quench in a “pseudo-equilibrium” in which some or all defect concentrations are frozen in until only the Fermi level E F remains to be equilibrated. However, this protocol does not account for the possibility of site exchanges which can create important defect redistributions as long as short-range defect migration is kinetically permissible. To model this redistribution, we developed an approach to solve for the non-equilibrium chemical potentials as a function of temperature while maintaining the overall defect stoichiometry. We then apply this approach to the Dirac semimetal Cd 3 As 2 to model extrinsic doping with group 1/11 and 14 elements. Undoped Cd 3 As 2 exhibits an undesirable mismatch between E F and the Dirac point. This unintentional electron doping originates from intrinsic defects and is difficult to overcome through adjustment of synthesis conditions alone. Employing our pseudo-equilibrium modeling, we identify extrinsic doping strategies for realizing doping-balanced Cd 3 As 2 at the relatively low temperatures accessible in thin-film growth of this material.

36 MATERIALS SCIENCE↗
NASA NTRS2015

Inconsistent Strategies to Spin up Models in CMIP5: Implications for Ocean Biogeochemical Model Performance Assessment

During the fifth phase of the Coupled Model Intercomparison Project (CMIP5) substantial efforts were made to systematically assess the skill of Earth system models. One goal was to check how realistically representative marine biogeochemical tracer distributions could be reproduced by models. In routine assessments model historical hindcasts were compared with available modern biogeochemical observations. However, these assessments considered neither how close modeled biogeochemical reservoirs were to equilibrium nor the sensitivity of model performance to initial conditions or to the spin-up protocols. Here, we explore how the large diversity in spin-up protocols used for marine biogeochemistry in CMIP5 Earth system models (ESMs) contributes to model-to-model differences in the simulated fields. We take advantage of a 500-year spin-up simulation of IPSL-CM5A-LR to quantify the influence of the spin-up protocol on model ability to reproduce relevant data fields. Amplification of biases in selected biogeochemical fields (O2, NO3, Alk-DIC) is assessed as a function of spin-up duration. We demonstrate that a relationship between spin-up duration and assessment metrics emerges from our model results and holds when confronted with a larger ensemble of CMIP5 models. This shows that drift has implications for performance assessment in addition to possibly aliasing estimates of climate change impact. Our study suggests that differences in spin-up protocols could explain a substantial part of model disparities, constituting a source of model-to- model uncertainty. This requires more attention in future model intercomparison exercises in order to provide quantitatively more correct ESM results on marine biogeochemistry and carbon cycle feedbacks.

histories↗
DOE OSTI2015

Recombination enhances HIV-1 envelope diversity by facilitating the survival of latent genomic fragments in the plasma virus population

HIV-1 is subject to immune pressure exerted by the host, giving variants that escape the immune response an advantage. Virus released from activated latent cells competes against variants that have continually evolved and adapted to host immune pressure. Nevertheless, there is increasing evidence that virus displaying a signal of latency survives in patient plasma despite having reduced fitness due to long-term immune memory. We investigated the survival of virus with latent envelope genomic fragments by simulating within-host HIV-1 sequence evolution and the cycling of viral lineages in and out of the latent reservoir. Our model incorporates a detailed mutation process including nucleotide substitution, recombination, latent reservoir dynamics, diversifying selection pressure driven by the immune response, and purifying selection pressure asserted by deleterious mutations. We evaluated the ability of our model to capture sequence evolution in vivo by comparing our simulated sequences to HIV-1 envelope sequence data from 16 HIV-infected untreated patients. Empirical sequence divergence and diversity measures were qualitatively and quantitatively similar to those of our simulated HIV-1 populations, suggesting that our model invokes realistic trends of HIV-1 genetic evolution. Moreover, reconstructed phylogenies of simulated and patient HIV-1 populations showed similar topological structures. Our simulation results suggest that recombination is a key mechanism facilitating the persistence of virus with latent envelope genomic fragments in the productively infected cell population. Recombination increased the survival probability of latent virus forms approximately 13-fold. Prevalence of virus with latent fragments in productively infected cells was observed in only 2% of simulations when we ignored recombination, while the proportion increased to 27% of simulations when we allowed recombination. We also found that the selection pressures exerted by different fitness landscapes influenced the shape of phylogenies, diversity trends, and survival of virus with latent genomic fragments. Furthermore, our model predicts that the persistence of latent genomic fragments from multiple different ancestral origins increases sequence diversity in plasma for reasonable fitness landscapes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Two-Stage Fracturing Wastewater Management in Shale Gas Development

Here, management of shale gas wastewater treatment, disposal, and reuse has become a significant environmental challenge, driven by an ongoing boom in development of U.S. shale gas reservoirs. Systems-analysis based decision support is helpful for effective management of wastewater, and provision of cost-effective decision alternatives from a whole-system perspective. Uncertainties are inherent in many modeling parameters, affecting the generated decisions. In order to effectively deal with the recourse issue in decision making, in this work a two-stage stochastic fracturing wastewater management model, named TSWM, is developed to provide decision support for wastewater management planning in shale plays. Using the TSWM model, probabilistic and nonprobabilistic uncertainties are effectively handled. The TSWM model provides flexibility in generating shale gas wastewater management strategies, in which the first-stage decision predefined by decision makers before uncertainties are unfolded is corrected in the second stage to achieve the whole-system’s optimality. Application of the TSWM model to a comprehensive synthetic example demonstrates its practical applicability and feasibility. Optimal results are generated for allowable wastewater quantities, excess wastewater, and capacity expansions of hazardous wastewater treatment plants to achieve the minimized total system cost. The obtained interval solutions encompass both optimistic and conservative decisions. Trade-offs between economic and environmental objectives are made depending on decision makers’ knowledge and judgment, as well as site-specific information. In conclusion, the proposed model is helpful in forming informed decisions for wastewater management associated with shale gas development.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2021

Analysis and Integration of the Hydraulic Fracturing Test Site-2 (HFTS-2) Comprehensive Dataset

Hydraulic Fracturing Test Site-2 (HFTS-2) is a field-based research experiment performed in the Permian (Delaware) Basin. The unique aspect of this program was the acquisition of a unique, comprehensive, diagnostic dataset. Additionally, shorter parent wells drilled three years before the child wells offered clear distinction between the stages influenced by parent-child effects and the stages without any effects. The goal of this study was to analyze and integrate this comprehensive diagnostic dataset to understand the areal and vertical extent of hydraulic fractures (HF). The paper also provides insights on the effects of parent wells’ depletion on child well HF geometry based on various monitoring methods and subsurface models. Areal and vertical coverage for all HFTS-2 wells during stimulation and depletion was estimated based on analysis and interpretation of diagnostics and advanced modeling results. HFTS-2 diagnostics included microseismic (MS), pre- and post-stimulation logs and cores, bottomhole gauges, and fiber optic (FO) data. The diagnostics results (MS, FO, image logs) were integrated and used to calibrate subsurface models. Additional field tests were designed and implemented for depletion monitoring. The tailored program for monitoring depletion included vertical and slant well pressures, interference testing, and a vertical strain depletion trial. Areal Coverage: Conventional MS (and FO MS) were used to compute HF dimensions, which were compared with diagnostics (FO strain, gauge, image logs) observations and calibrated subsurface models. A post-production interference test did not show offset well communication. Vertical Coverage: Vertical coverage during stimulation was monitored using a vertical monitoring well. The stronger mechanical strain signals showed good correlation with MS event intensities, geomechanical properties, and gauge inferences. Vertical depletion was estimated based on vertical/slant well gauges and strain depletion tests. Parent-Child Effects: Diagnostics and calibrated subsurface models show asymmetry in child well fracture geometries for stages that overlap parent wells. Child well image logs serve as a good indicator for parent Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/2-21URTC/D021S031R004/2477423/urtec-2021-5241-ms.pdf/1 by Carol Worster on 28 February 2022 URTeC 5241 well HF tracking. Child well MS events had an eastward bias, in line with pre-stimulation image logs, and was confirmed by parent well frac hits. Novel/Additive Information: The dataset presents a unique, over-constrained problem space to compare independent techniques to arrive at HF metrics (i.e., stimulation height and/or half-length), unlike a single source dataset, in which calibration is done using available data to guide predictions. Here, the asymmetry in HF geometry seen in the stages influenced by parent-child effects offers unique insights into well spacing and landing, which are key capital decisions the unconventional resources industry is seeking to optimize.

58 GEOSCIENCES↗
DOE OSTI2025

Modeling hydraulic fracture opening and closure with proppant transport and settlement

Hydraulic fracturing is a widely used reservoir stimulation technique for improving fluid circulation in rock formations with extremely low permeability, particularly in enhanced geothermal systems (EGS). To better understand the complex processes involved and improve hydraulic stimulation performance, we have developed ELK (ELectrical fracKing), a MOOSE-based 3D finite element application designed to model the behavior of proppant-fluid mixtures in propagating fractures. ELK integrates both the fluid and proppant components, incorporating particle-driven processes such as gravity settling, particle-particle interactions, and strong density and viscosity contracts, in addition to conventional fluid-driven fracture propagation. In this contribution, we extend ELK to model propped fracture closure, which occurs after the injection phase due to a dramatic drop in the effective stress on the fracture plane. During the shut-in, flowback, and production periods, the fracture width decreases, with the closure behavior depending on proppant concentration.. At low concentrations, closure follows a nonlinear joint law linked to the stiffness of asperities in the fracture walls. While at high concentrations, it is controlled by the properties of packed proppant bed. The extended ELK application is validated against several benchmark examples, including the propagation of an inclined frictional crack, fracture opening and sliding in response to fluid injection, and flowback analysis. We believe that ELK’s enhanced capabilities can serve as a valuable tool for the design and optimization of EGS deployment.

15 - GEOTHERMAL ENERGY↗
DOE OSTI2020

North Dakota Integrated Carbon Storage Complex Feasibility Study. Final report

In spring 2017, the Energy & Environmental Research Center (EERC) initiated an effort to determine the feasibility of developing a commercial-scale CO 2 geologic storage complex able to store 50+ million tonnes (Mt) of CO 2 in central North Dakota safely, permanently, and economically. The objective was to fulfill the goals of the U.S. Department of Energy (DOE) Carbon Storage Assurance Facility Enterprise (CarbonSAFE) Initiative and address technical and nontechnical challenges specific to commercial-scale deployment of a CO2 storage project. The findings clearly show that the concept of capturing CO 2 from a lignite-fired electrical generation facility in central North Dakota and safely and permanently storing the CO 2 in the deep subsurface is indeed technically, economically, and socially feasible. This project evaluated two study areas and their respective geologic storage complexes located adjacent to separate coal-fired facilities in North Dakota: the Basin Electric Power Cooperative (BEPC)-owned Great Plains Synfuels Plant (GPSP) and the Minnkota Power Cooperative (Minnkota)-owned Milton R. Young Station (MRYS). These locations, one with CO 2 capture in place and an existing CO 2 pipeline, are bolstered by progressive North Dakota pore space ownership and long-term liability laws. These elements and a motivated team created an ideal synergistic scenario for ensuring success of the CarbonSAFE Initiative and promoting North Dakota’s statewide vision for carbon management. The project included drilling two new geologic characterization wells, integrating an existing 3-D seismic survey, creating a geologic model subsequently used for injection simulation, a risk assessment, public outreach, and generating a site development plan based on results. In addition, the performance of select National Risk Assessment Partnership tools was evaluated. The geologic characterization wells were drilled ~5600 feet deep to the Broom Creek Formation; ~350 feet of core was retrieved from each well. The core included the Broom Creek (targeted injection zone) and a portion of the overlying Opeche Shale (seal). The Flemmer-1 well, west of Beulah, North Dakota, yielded 169 feet of sandstone. The BNI-1 well located south of Center, North Dakota, yielded 124 feet of sandstone. In each case, laboratory analysis of the sandstone showed permeability in the 300–1000-mD range, with porosity of 20%–30%. The Flemmer-1 well was sited within the boundaries of an existing 3-D seismic survey. Colocating the well with the seismic survey maximized the relationship between new and legacy data and developed a first-of-its-kind interpretation of the geologic fabric of the Broom Creek. Geologic characterization data were integrated into a 5544-mi 2 geocellular model that encompassed both new wells and stratigraphy from the surface to the Amsden Formation (underlying the Broom Creek). The model was later expanded vertically to include the deeper Black Island–Deadwood interval to examine its potential viability as a storage target. The geocellular model provided the foundation for dynamic simulation of CO 2 into the Broom Creek. Results of the simulation suggest that the Broom Creek could accept the DOE target rate of 2 Mt/yr of CO 2 into as few as two wells. To bracket the expected capture from MRYS, simulations were also investigated for a 4-Mt/yr rate near MRYS. Although more wells are needed (two additional), the Broom Creek still has the storage resource to accept the CO 2 at the increased rate. A risk assessment exercise was conducted to identify and assess technical and nontechnical risks that could prevent potential candidate storage complexes within the study area from serving as commercial storage sites. The assessment identified and evaluated six technical risk categories: 1) CO 2 injectivity, 2) storage capacity, 3) lateral migration of CO 2 , 4) lateral pressure propagation, 5) vertical migration of CO 2 or formation brine, and 6) induced seismicity. Following two rounds of analysis and scoring, no risks were determined to preclude continued efforts to develop carbon capture, utilization, and storage (CCUS) in central North Dakota. The risk assessment results will be used to guide future site characterization, modeling and simulation, and monitoring activities. A specific nontechnical strategic risk based on challenges that may be realized in amalgamation of pore space resulted in vertically expanding the geologic model to incorporate the potential for stacked storage in multiple saline reservoirs. By including the Black Island–Deadwood interval (the basal sedimentary reservoir in this region), the amalgamated areal extent could be reduced by as much as 45%. Working with a smaller geographic area reduces risks and costs associated with monitoring and pore space leasing. An economic evaluation incorporating capture; transport (<5 mi); Class VI wells; permitting; and monitoring, verification, and accounting suggests implementing commercial-scale CCUS is economically attractive if the federal tax benefits of 45Q are included. This is validated by Minnkota’s continued pursuit of CO 2 capture and geologic storage at MRYS through its Project Tundra initiative, indicating that there is a business case for CCUS in central North Dakota. Currently, North Dakota is the only state with underground injection control (UIC) Class VI primacy. Built into the North Dakota Century Code is a series of regulatory requirements that guide the process to obtain a Class VI CO 2 storage facility permit. As part of this project, a site development plan was compiled to assure compliance with North Dakota’s requirements to permit a commercial-scale CO 2 storage operation and includes a prospective time line encompassing a general breakdown of activities. In total, an estimated 30 months would be needed to execute the necessary steps to attain a North Dakota CO 2 storage facility permit. Outreach was an integral part of the project and encompassed any project-related activity that had contact or exposure beyond the project team. The goals of outreach were to foster an environment from which stakeholders could make informed decisions about the project and gauge community receptiveness to a CCUS project. A consistent set of messages and outreach products were developed in conjunction with an outreach advisory board that integrated project partners and team members. To gauge public acceptability of geologically storing CO 2 , 5611 households in the project area were invited to participate in an online survey. The survey results indicate that the public attitude regarding CCUS is neutral to positive, with strong sentiment that CO 2 capture and storage may be an approach to maintain the economic vitality of the region. To achieve project objectives, critical support in the form of financial backing, engineering evaluations, site access, outreach collaboration, operations data, risk assessment/evaluation, and software access was provided by BEPC, the North Dakota Industrial Commission Lignite Research Council, ALLETE Clean Energy, BNI Energy, North American Coal Corporation, Minnkota, Prairie Public Broadcasting, Computer Modelling Group Ltd., and Schlumberger.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects by Incorporating Machine Learning Techniques

The overall objective of this project is to develop a framework for reliable characterization and prediction of the state of stress in the overburden and underburden (including the basement) in CO 2 storage reservoirs using machine learning and integrated geomechanics and geophysical methods. Specifically, we propose to develop workflow encompassing of technologies and/or methods to predict stress and pressure changes due to CO 2 injection in an active tertiary recovery site and their impacts on subtle fault activation, fractures and occurrence of microseismic events and compare responses to field observations. In this project, we anticipate using dataset from the Farnsworth field Unit (FWU) which is operated by Purdure Petroleum. A novel elastic-waveform VSP inversion technique will be used to estimate high-resolution spatial and temporal changes of elastic moduli in CO 2 storage reservoirs, which will be combined with velocity-stress relationship derived from laboratory tests to obtain subsurface pressure and stress. Clustered microseismic data will be jointly inverted for improved focal mechanisms. Least-squares reverse-time migration of microseismic waveform data will be performed to directly image fracture/fault zones. Additionally, a deep neural network machine learning technique with convolutional and recurrent layers will be used for learning the spectro-temporal structures in microseismic waveforms. The results of this geotechnical data analysis will be integrated to develop a high-resolution 3D mechanical earth model extending from the overburden sealing formations to the underburden including the basement. Mechanical properties will be derived through integration of mechanical logs, tests, available results from chemo-mechanical laboratory tests, and elastic inversion of seismic data using a combination of Bayesian and stochastic methods as well as machine learning technique. Failure features (faults/fractures) will be represented and/or modeled based on seismic and core data analysis. A transient hydrodynamic-geomechanical model will be developed through coupling with the calibrated FWU reservoir simulation model. The full physics coupled model will be used to train a reduced order proxy model using machine learning algorithm for estimating stress which will then be used with appropriate constitutive relationships and forward seismological models to simulate pressure changes and induced microseismicity. An advanced optimization framework will be developed to perform a history match to minimize error between field observations and simulated. The history matched proxy model will be verified against the full-physics equivalent. The field observations that will be used in the coupled model calibration process include pressure/stress inverted from VSP, moment magnitude from microseismic analysis, real time downhole pressure measurements, production and injection data. Parameter sensitivity and uncertainty analysis will be performed to characterize the impact of model parameter uncertainty on stress estimates. The proposed project will have significant impact on future field implementation of the proposed technology. Because the project field site is an ongoing CO 2 EOR development, the value of the new technology will be demonstrated in an operational context and evaluated as a viable risk mitigation strategy. Cost/benefit will be evaluated together with the various commercial incentives for CO 2 sequestration available to oil and gas operators. The extensive available dataset and ongoing data acquisition under the SWP Phase III work plan provides flexibility for investigation of multiple approaches and reduces technical risk.

58 GEOSCIENCES↗