SEARCH · Engineering Papers

Results for “Bottom hole”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1985

SSME seal test program: Test results for hole-pattern damper seals

The results consisting of direct and transverse force coefficients are presented for thirteen, hole-pattern, damper-seal configurations. The designation damper seal refers to a seal which uses a deliberately roughened stator nd smooth rotor, to increase the net damping force developed by a seal. The designation hole-pattern refers to a stator roughness pattern which is developed by a pattern of round holes while are milled into the stator. All seals tested use the same smooth rotor and have the same constant minimum clearance. The seal tests examined the following major design options: (1) hole-area density, i.e., the proportion of stator surface area consumed by holes; and (2) hole depth, particularly the ratio of hole depth to minimum clearance. In addition, limited data were taken to examine the influence of in-line versus staggered hole patterns and flat-bottomed versus spherical-bottomed holes.

Childs, D. W.↗

DOE OSTI2023

Low-Temperature Geothermal Geospatial Datasets: An Example from Alaska

This project is a component of a broader effort focused on geothermal heating and cooling (GHC) with the aim of illustrating the numerous benefits of incorporating GHC and geothermal heat exchange (GHX) into community energy planning and national decarbonization strategies. To better assist private sector investment, it is currently necessary to define and assess the potential of low-temperature geothermal resources. For shallow GHC/GHX fields, there is no formal compilation of subsurface characteristics shared among industry practitioners that can improve system design and operations. Alaska is specifically noted in this work, because heretofore, it has not received a similar focus in geothermal potential evaluations as the contiguous United States. The methodology consists of leveraging relevant data to generate a baseline geospatial dataset of low-temperature resources (less than 150 degrees C) to compare and analyze information accessible to anyone trying to understand the potential of GHC/GHX and small-scale low-temperature geothermal power in Alaska (e.g., energy modelers, communities, planners, and policymakers). Importantly, this project identifies data related to (1) the evaluation of GHC/GHX in the shallow subsurface, and (2) the evaluation of low-temperature geothermal resource availability. Additionally, data is being compiled to assess repurposing of oil and gas wells to contribute co-produced fluids toward the geothermal direct use and heating and cooling resource potential. In this work we identified new data from three different datasets of isolated geothermal systems in Alaska and bottom-hole temperature data from oil and gas wells that can be leveraged for evaluation of low-temperature geothermal resource potential. The goal of this project is to facilitate future deployment of GHC/GHX analysis and community-led programs and update the low-temperature geothermal resources assessment of Alaska. A better understanding of shallow potential for GHX will improve design and operations of highly efficient GHC systems. The deployment and impact that can be achieved for low-temperature geothermal resources will contribute to decarbonization goals and facilitate widespread electrification by shaving and shifting grid loads.

15 GEOTHERMAL ENERGY↗

DOE OSTI2020

Analysis of Geothermal Resources in Three Texas Counties

This project updates the geothermal resources beneath our oil and gas fields, as part of the research for the Texas GEO project. This report "Analysis of Geothermal Resources in Three Texas Counties" (October 2020) improves on previous mapping of the Texas resources for the counties of Crockett (West Texas), Jackson (central Gulf Coast) and Webb (South Texas). Through additional bottom-hole temperatures (BHT), the number of well sites increased from 532 to 5,410 in total for these counties. The improved methodology to calculate formation temperatures from 3.5 km (11,500 ft) to 10 km (32,800 ft) includes thermal conductivity values more closely related to the actual county geological formations, incorporated radiogenic heat production of formations, and the related mapped depth to basement. The results show deep temperatures as hotter than previously calculated, with temperatures of 150 degrees Celcius possible for Webb County between depths of 2.6 - 5.1 km, Jackson County between depths 3.0 - 5.4 km, and Crockett County between depths of 2.7 - 8.0 km.

10 km↗

DOE OSTI2018

SECARB Pressure Monitoring Data

Bottom-hole, above zone monitoring interval, and injection zone pressure data collected during the SECARB project in Cranfield, Mississippi to assess the relationship between pressure field and multiphase field. Submission includes 10-second interval data from Detailed Area of Study wells: CFU31-F1 (injector), CFU31-F2 (observation), CFU31-F3 (observation) and Ella G Lees no. 7 (observation) well located west of the DAS. Associated Publications: Joy, C. A., 2011, The effects of pressure variation and chemical reactions on the elasticity of the lower Tuscaloosa sandstone of the Cranfield Field Mississippi, The University of Texas at Austin, Master’s thesis, 97 p. Kim, S., and Hosseini, S. A., 2013, Above-zone pressure monitoring and geomechanical analysis of a field scale CO2 injection, Cranfield Mississippi, Greenhouse Gases: Science and Technology, doi:10.1002/ghg.1388. Kim, S., and Hosseini, S. A., 2017, Study on the ratio of pore-pressure/stress changes during fluid injection and its implications for CO2 geologic storage: Journal of Petroleum Science and Engineering, v. 149, p. 138-150, doi:10.1016/j.petrol.2016.10.037. Mathias, S. A., Gluyas, J. G., Gonzalez Martinez de Miguel, G. J., and Hosseini, S. A., 2011, Role of partial miscibility on pressure buildup due to constant rate injection of CO2 into closed and open brine aquifers: Water Resources Research, v. 47, W12525, 11 p., doi:10.1029/2011WR011051. Meckel, T. A., Zeidouni, M., Hovorka, S. D., and Hosseini, S.A., 2013, Assessing sensitivity to well leakage from three years of continuous reservoir pressure monitoring during CO2 injection at Cranfield, MS, USA: International Journal of Greenhouse Gas Control, [insert volume no., page numbers], doi:10.1016/j.ijggc.2013.01.019. Nicot, J.-P., Oldenburg, C. M., Bryant, S. L., and Hovorka, S. D., 2009, Pressure perturbations from geologic carbon sequestration: area-of-review boundaries and borehole leakage driving forces, in Energy Procedia (v. 1, no.1), Proceedings of 9th International Conference on Greenhouse Gas Control Technologies, GHGT9, 16–20 November, Washington DC, p. 47–54. Tao, Q., Bryant, S. L., and Meckel, T. A., 2013, Modeling above-zone measurements of pressure and temperature for monitoring CCS sites: International Journal of Greenhouse Gas Control, v. 18, p. 523–530, doi:10.1016/j.ijggc.2012.08.011.

above zone monitoring interval↗

DOE OSTI2014

Apparatus for advancing a wellbore using high power laser energy

Delivering high power laser energy to form a borehole deep into the earth using laser energy. Down hole laser tools, laser systems and laser delivery techniques for advancement, workover and completion activities. A laser bottom hole assembly (LBHA) for the delivery of high power laser energy to the surfaces of a borehole, which assembly may have laser optics, a fluid path for debris removal and a mechanical means to remove earth.

47 OTHER INSTRUMENTATION↗

NASA NTRS1986

Ultrasonic level sensors for liquids under high pressure

An ultrasonic level sensor of novel design continuously measures the level of a liquid subjected to a high pressure (up to about 40 MPa), as is sometimes required for the effective transfer of the liquid. The sensor operates as a composite resonator fabricated from a standard high-pressure plug. A flat-bottom hole is machined into the plug along its center line. An ultrasonic transducer is bonded rigidly to the interior surface of the bottom wall, while the exterior surface is in contact with the liquid. Although the bottom wall is designed to satisfy the pressure code, it is still sufficiently thin to permit ready excitation of the axisymmetric plate modes of vibration. The liquid level is measured by a conventional pulse-echo technique. A prototype sensor was tested successfully in a 2300-l water vessel at pressures up to about 37 MPa. A spectral analysis of the transmitted pulse reveals that the flexural, extensional, thickness-shear, and radial plate modes are excited into vibration, but none of these appears to be significantly affected by the pressurization of the liquid.

Zuckerwar, A. J.↗

DOE OSTI2018

Discovery and analysis of a blind geothermal system in southeastern Gabbs Valley, Western Nevada

This study assessed the potential for high-temperature (≥130°C) blind geothermal systems in southeastern Gabbs Valley by employing a geothermal play fairway analysis approach involving integration of geologic, geophysical, and geochemical datasets. Gabbs Valley is a complex, tectonically active basin within the Great Basin on the boundary between the transtensional central Walker Lane domain and extensional Basin and Range province. The termination of the Petrified Springs fault, a major dextral fault of the central Walker Lane, in an array of normal faults indicates that southeastern Gabbs Valley occupies a displacement transfer zone, which is a favorable structural setting for geothermal activity. The displacement transfer zone is a structurally complicated area with many faults of varying geometries, kinematics, and ages. A substantial northwest-trending gravity high within the south-central part of southeastern Gabbs Valley is produced from offsets along concealed northwest-striking dextral-normal faults that intersect strands of north-northeast-to east-northeast-striking normal faults. This area of tightly spaced faults and fault intersections is a particularly favorable structural setting within the broader-scale displacement transfer zone. Multiple lines of direct and indirect evidence suggest the presence of a relatively high-temperature (≥130°C) blind geothermal system in this area, including collocated intersecting gravity gradients, magnetic-low, low-resistivity, and 2-m temperature anomaly. Potentially related, water samples from agricultural wells ~7 km northwest of the 2-m anomaly yield geothermometers indicating subsurface fluid temperatures of 130-140°C. Six temperature-gradient holes were drilled to target the extent of the shallow-temperature and geophysical anomalies to define the geothermal system in southeastern Gabbs Valley. Two wells contained high temperatures exceeding boiling with bottom-hole temperatures ii of 114.5°C and 124.9°C, respectively, and the remaining wells displayed elevated to background temperatures ranging from 79.2°C to 28.7°C. The observed temperature gradients for the two hottest drill holes necessitate drill-hole intercepts of convecting hydrothermal fluids. The 2-m survey is spatially and thermally consistent with the temperature-gradient holes and constrains the core of upflow at a shallow depth interval (<250m) to an area <0.5-1 km west-southwest of the hottest drill hole. The zone of inferred upwelling occupies the area corresponding to the center of the magnetic-low. This project has several significant implications. The newly discovered blind geothermal system in southeastern Gabbs Valley may be capable of supporting a binary-type geothermal power plant. This detailed study provides an initial validation of the play fairway approach to geothermal exploration. It also demonstrates both the broad applicability of the play fairway exploration strategy and the large untapped potential for commercial-grade blind geothermal systems in many of the basins in the Great Basin region.

15 GEOTHERMAL ENERGY↗

DOE OSTI2014

Methods and apparatus for removal and control of material in laser drilling of a borehole

The removal of material from the path of a high power laser beam during down hole laser operations including drilling of a borehole and removal of displaced laser effected borehole material from the borehole during laser operations. In particular, paths, dynamics and parameters of fluid flows for use in conjunction with a laser bottom hole assembly.

47 OTHER INSTRUMENTATION↗

DOE OSTI2016

Methods and apparatus for removal and control of material in laser drilling of a borehole

15 GEOTHERMAL ENERGY↗

DOE OSTI2016

West Flank Coso FORGE: 3D Temperature Model

x,y,z data of the 3D temperature model for the West Flank Coso FORGE site. Model grid spacing is 250m. The temperature model for the Coso geothermal field used over 100 geothermal production sized wells and intermediate-depth temperature holes. At the near surface of this model, two boundary temperatures were assumed: (1) areas with surface manifestations, including fumaroles along the northeast striking normal faults and northwest striking dextral faults with the hydrothermal field, a temperature of ~104 deg C was applied to datum at +1066 meters above sea level elevation, and (2) a near-surface temperature at about 10 meters depth, of 20 deg C was applied below the diurnal and annual conductive temperature perturbations. These assumptions were based on heat flow studies conducted at the CVF and for the Mojave Desert. On the edges of the hydrothermal system, a 73 deg C/km temperature gradient contour was established using conductive gradient data from shallow and intermediate-depth temperature holes. This contour was continued to all elevation datums between the 20 deg C surface and -1520 meters below mean sea level. Because the West Flank is outside of the geothermal field footprint, during Phase 1, the three wells inside the FORGE site were incorporated into the preexisting temperature model. To ensure a complete model was built based on all the available data sets, measured bottom-hole temperature gradients in certain wells were downward extrapolated to the next deepest elevation datum (or a maximum of about 25% of the well depth where conductive gradients are evident in the lower portions of the wells). After assuring that the margins of the geothermal field were going to be adequately modelled, the data was contoured using the Kriging method algorithm. Although the extrapolated temperatures and boundary conditions are not rigorous, the calculated temperatures are anticipated to be within ~6 deg C (20 deg F), or one contour interval, of the observed data within the Coso geothermal field. Based on a lack of temperature data west of 74-2TCH, the edges of this model still seem to have an effect on West Flank modeled temperatures.

15 GEOTHERMAL ENERGY↗

DOE OSTI2023

Downhole Sensing and Event-Driven Sensor Fusion for Depth-of-Cut Based Autonomous Fault Response and Drilling Optimization

Achieving robust and efficient drilling is a critical part of reducing the cost of geothermal energy exploration and extraction. Drilling performance is often evaluated using one or more of three key metrics: depth of cut (DOC), rate of penetration (ROP), and mechanical specific energy (MSE). All three of these quantities are related to each other. DOC refers to the depth a bit penetrates into rock during drilling. This is an important quantity for estimating bit behavior. ROP is the simply the DOC multiplied by the rotational rate, and represents how quickly the drill bit is advancing through the ground. ROP is often the parameter used for drilling control and optimization. Finally, MSE provides insight into drilling efficiency and rock type. MSE calculations rely on ROP, drilling force, and drilling torque. Surface-based sensors at the top of the drill are often used to measure all these quantities. However, top-hole measurements can deviate substantially from the behavior at the bit due to lag, vibrations, and friction. Therefore, relying only on top-hole information can lead to suboptimal drilling control. In this work, we describe recent progress towards estimating ROP, DOC, and MSE using down-hole sensing. We assume down-hole measurements of torque, weight-on-bit (WOB). Our hypothesis is that these measurements can provide more rapid and accurate measures of drilling performance. We show how a multi-layer perceptron (MLP) machine learning algorithm can provide rapid and accurate performance when evaluated on experimental data taken from Sandia’s Hard Rock Drilling Facility. In addition, we implement our algorithms on an embedded system intended to emulate a bottom-hole-assembly for sensing and estimation. Our experimental results show that DOC can be estimated accurately and in real-time. These estimates when combined with measurements for rotary speed, torque, and force can provide improved estimates for ROP and MSE. These results have the potential to enable better drilling assessment, improved control, and extended component lifetimes.

15 GEOTHERMAL ENERGY↗

DOE OSTI2021

Borehole breakout modeling in arkose and granite rocks

The existence of a deep borehole in the Earth’s crust disturbs the local stresses and creates a stress concentration that may result in breakout and damage to the borehole. Maintaining wellbore integrity mitigates environmental impacts such as groundwater contamination, gas leakage to the atmosphere, and fluid spills and seepage at the surface. In this paper, the stability of deep boreholes (5 km) is examined by laboratory experiments and numerical models in the context of nuclear waste disposal in Israel. Two rock types in southern Israel are considered: the crystalline basement (granite) and the Zenifim Formation (arkose). A series of room-temperature triaxial rock deformation experiments were conducted at different confining pressures. This mechanical characterization was then used to parameterize the elastic properties and damage behavior of the rocks. This facilitated modeling the stability of the deep boreholes by two different formulations of damage rheology: a dynamic-oriented formulation used to model deformation immediately after the creation of the open hole and a quasi-static formulation used to model longer stress corrosion regime. The calibrated modeling results indicate greater stability with Zenifim arkose than the crystalline granite for deep borehole conditions despite the granite having a greater triaxial compressive strength. Dissipation associated with dilation and porous compaction in the arkose during deformation plays a significant stabilizing role in the borehole compared to crystalline rocks. These results suggest that common strength-based borehole stability assessment may lead to inaccurate predictions. Three-dimensional modeling of bottom-hole stress conditions and the effects of transient borehole geometry show conventional two-dimensional analysis may not be conservative when predicting borehole damage.

58 GEOSCIENCES↗

DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗

DOE OSTI2016

Co-Optimization of CO 2 -EOR and Storage Processes in Mature Oil Reservoirs

This article presents an optimization methodology for CO 2 enhanced oil recovery in partially depleted reservoirs. A field-scale compositional reservoir flow model was developed for assessing the performance history of an active CO 2 flood and for optimizing both oil production and CO 2 storage in the Farnsworth Unit (FWU), Ochiltree County, Texas. A geological framework model constructed from geophysical, geological, and engineering data acquired from the FWU was the basis for all reservoir simulations and the optimization method. An equation of state was calibrated with laboratory fluid analyses and subsequently used to predict the thermodynamic minimum miscible pressure (MMP). Initial history calibrations of primary, secondary and tertiary recovery were conducted as the basis for the study. After a good match was achieved, an optimization approach consisting of a proxy or surrogate model was constructed with a polynomial response surface method (PRSM). The PRSM utilized an objective function that maximized both oil recovery and CO 2 storage. Experimental design was used to link uncertain parameters to the objective function. Control variables considered in this study included: water alternating gas cycle and ratio, production rates and bottom-hole pressure of injectors and producers. Other key parameters considered in the modeling process were CO 2 purchase, gas recycle and addition of infill wells and/or patterns. The PRSM proxy model was ‘trained’ or calibrated with a series of training simulations. This involved an iterative process until the surrogate model reached a specific validation criterion. A sensitivity analysis was first conducted to ascertain which of these control variables to retain in the surrogate model. A genetic algorithm with a mixed-integer capability optimization approach was employed to determine the optimum developmental strategy to maximize both oil recovery and CO 2 storage. The proxy model reduced the computational cost significantly. The validation criteria of the reduced order model ensured accuracy in the dynamic modeling results. The prediction outcome suggested robustness and reliability of the genetic algorithm for optimizing both oil recovery and CO 2 storage. The reservoir modeling approach used in this study illustrates an improved approach to optimizing oil production and CO 2 storage within partially depleted oil reservoirs such as FWU. Lastly, this study may serve as a benchmark for potential CO 2 –EOR projects in the Anadarko basin and/or geologically similar basins throughout the world.

02 PETROLEUM↗

DOE OSTI2019

On producing CO 2 from subsurface reservoirs: simulations of liquid-gas phase change caused by decompression

Carbon dioxide (CO 2 ) extraction from deep reservoirs is currently important in CO 2 enhanced oil recovery (EOR) and may become more important in the future if interim CO 2 storage becomes common. In late 2014, we were involved in a production test of liquid CO 2 from the Middle Duperow dolostone at Kevin Dome, Montana. The test resulted in lowering the temperature at the well bottom to ~2 °C, and showed that the well and reservoir had very low CO 2 productivity. We have used the CO 2 modeling capabilities of the TOUGH codes to simulate the test and to show that liquid CO 2 in the reservoir changes to gas phase as the pressure is lowered in the well during production testing. The associated phase change and decompression combine to drastically lower the bottom-hole temperature, creating the potential for water ice or CO 2 hydrate to form. By hypothesizing a relatively high-permeability damage zone near the well surrounded by lower permeability reservoir rock, we can match the observed pressure, temperature, and production rate. Moving from the Kevin Dome test to the question of CO 2 extraction from deep reservoirs in general, we carried out a parametric study to investigate the effects of reservoir depth and transmissivity on CO 2 production rate for a prototypical reservoir. Simulations show that depth and high transmissivity favor productivity. Complex phase changes within the ranges of P-T encountered in typical CO 2 production wells affect production rates. The results of our parametric study may be useful for the preliminary feasibility assessment of CCO 2 extraction from deep reservoirs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2017

Error modeling for surrogates of dynamical systems using machine learning: Machine-learning-based error model for surrogates of dynamical systems

A machine learning–based framework for modeling the error introduced by surrogate models of parameterized dynamical systems is proposed. The framework entails the use of high-dimensional regression techniques (eg, random forests, and LASSO) to map a large set of inexpensively computed “error indicators” (ie, features) produced by the surrogate model at a given time instance to a prediction of the surrogate-model error in a quantity of interest (QoI). This eliminates the need for the user to hand-select a small number of informative features. The methodology requires a training set of parameter instances at which the time-dependent surrogate-model error is computed by simulating both the high-fidelity and surrogate models. Using these training data, the method first determines regression-model locality (via classification or clustering) and subsequently constructs a “local” regression model to predict the time-instantaneous error within each identified region of feature space. We consider 2 uses for the resulting error model: (1) as a correction to the surrogate-model QoI prediction at each time instance and (2) as a way to statistically model arbitrary functions of the time-dependent surrogate-model error (eg, time-integrated errors). We then apply the proposed framework to model errors in reduced-order models of nonlinear oil-water subsurface flow simulations, with time-varying well-control (bottom-hole pressure) parameters. The reduced-order models used in this work entail application of trajectory piecewise linearization in conjunction with proper orthogonal decomposition. Moreover, when the first use of the method is considered, numerical experiments demonstrate consistent improvement in accuracy in the time-instantaneous QoI prediction relative to the original surrogate model, across a large number of test cases. When the second use is considered, results show that the proposed method provides accurate statistical predictions of the time- and well-averaged errors.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Co-optimization of CO 2 -EOR and Storage Processes under Geological Uncertainty

This paper presents an integrated numerical framework to co-optimize EOR and CO 2 storage performance in the Farnsworth field unit (FWU), Ochiltree County, Texas. The framework includes a field-scale compositional reservoir flow model, an uncertainty quantification model and a neural network optimization process. The reservoir flow model has been constructed based on the field geophysical, geological, and engineering data. A laboratory fluid analysis was tuned to an equation of state and subsequently used to predict the thermodynamic minimum miscible pressure (MMP). A history match of primary and secondary recovery processes was conducted to estimate the reservoir and multiphase flow parameters as the baseline case for analyzing the effect of recycling produced gas, infill drilling and water alternating gas (WAG) cycles on oil recovery and CO 2 storage. A multi-objective optimization model was defined for maximizing both oil recovery and CO 2 storage. The uncertainty quantification model comprising the Latin Hypercube sampling, Monte Carlo simulation, and sensitivity analysis, was used to study the effects of uncertain variables on the defined objective functions. Uncertain variables such as bottom hole injection pressure, WAG cycle, injection and production group rates, and gas-oil ratio among others were selected. The most significant variables were selected as control variables to be used for the optimization process. A neural network optimization algorithm was utilized to optimize the objective function both with and without geological uncertainty. The vertical permeability anisotropy (Kv/Kh) was selected as one of the uncertain parameters in the optimization process. The simulation results were compared to a scenario baseline case that predicted CO 2 storage of 74%. The results showed an improved approach for optimizing oil recovery and CO 2 storage in the FWU. The optimization process predicted more than 94% of CO 2 storage and most importantly about 28% of incremental oil recovery. The sensitivity analysis reduced the number of control variables to decrease computational time. A risk aversion factor was used to represent results at various confidence levels to assist management in the decision-making process. The defined objective functions were proved to be a robust approach to co-optimize oil recovery and CO 2 storage. The Farnsworth CO 2 project will serve as a benchmark for future CO 2 –EOR or CCUS projects in the Anadarko basin or geologically similar basins throughout the world.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Investigation of enhanced CO 2 storage in deep saline aquifers by WAG and brine extraction in the Minnelusa sandstone, Wyoming

Geological CO 2 sequestration in deep saline aquifers has been extendedly investigated to reach the goal of carbon neutral as large amount of CO 2 can be reduced in a short time. Here, this study investigates the feasibility of water-alternating-gas (WAG) injection and brine extraction on enhancing CO 2 storage in a target deep saline aquifer, Minnelusa Sandstone in the Powder River Basin (PRB) of Wyoming. An integrated numerical model is developed to account for four scenarios: (1) CO 2 continuous injection (CI) through one injection well, (2) WAG injection through one injection well, (3) CO 2 CI through one injection well and brine extraction through one producer, (4) WAG injection through one injection well and brine extraction through one producer. Simulation results suggest that WAG injection and brine extraction corporately make it possible to enhance CO 2 injectivity while securing CO 2 storage safety. For instance, WAG injection considerably reduces structural trapping contribution while enhances dissolved and residual trapping contributions. As for brine extraction, it can decrease the maximum averaged reservoir pressure by 37% and 22% for CI and WAG injection, respectively. Besides, sensitive analyses of the operational parameters for the fourth scenarios are performed. Results reveal that, when the amount of total CO 2 injection has been predetermined, CO 2 injection time and rate within each period of WAG should be small whereas it is preferable to have large water injection time and rate. This is to secure desirable dissolved and residual trapping contributions (store CO 2 with safety). Sensitive degree results confirm that CO 2 injection rate, water injection time, and CO 2 injection time of WAG injection exhibit significant effects on CO 2 storage whereas the impacts of producer bottom-hole pressure (BHP), water injection rate, and CO 2 -water injection time are relatively weaker. This study not only sheds light on achieving double-win goal: enhance CO 2 injectivity and store CO 2 with safety, but also can be a critical reference for other CO 2 storage in deep saline aquifers.

58 GEOSCIENCES↗