Prediction of the maximum allowable bottom hole pressure in CO2 injection wells
Not Available
Guo, Boyun↗
SEARCH · Search NASA
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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Interpretation of bottom-hole temperature data from CO 2 injection projects
Although bottom-hole temperature (BHT) data is routinely collected with bottom-hole pressure (BHP) data during CO 2 injection in geological storage projects, it is seldom analyzed in a quantitative manner. However, the BHP and BHT signals do exhibit correlated behavior. This work proposes a practical approach for using BHT data from warmback (no injection) periods, in conjunction with the amplitude of BHP and BHT changes during the preceding injection event, to estimate formation permeability. The methodology is demonstrated using data from the AEP Mountaineer project in USA where CO 2 injection was carried out into a saline aquifer. Estimated permeability values from the proposed approach are in good agreement with those obtained from transient pressure analysis of BHP data from the falloff period.
42 ENGINEERING↗
SECARB Bottom Hole Gravity Data
Borehole gravity measurements obtained during the SECARB project at the Cranfield oil site in Mississippi from CFU31-F2 and CFU31-F3 wells. Data was used to calculate density changes within the Cranfield reservoir and to test borehole gravity performance compared to a variety of other methods for monitoring the injected CO2 plume. Associated Publications: Dodds, K., Krahenbuhl, R., Reitz, A., Li, Y., Hovorka, S. D., 2013, Evaluation of time lapse borehole gravity for CO2 plume detection SECARB Cranfield: International Journal of Greenhouse Gas Control.
Borehole↗
Directional antennas for electromagnetic mapping in a borehole
A bottom hole assembly used for a field operation is disclosed herein. The bottom hole assembly can include at least one directional antenna disposed on an outer surface of a first tubing pipe of a tubing string, where the at least one directional antenna receives a first electric current from at least one power source, where the first electric current generates a first magnetic field that radiates from the at least one directional antenna into a formation. The bottom hole assembly can also include at least one receiver disposed on a second tubing pipe of the tubing string, where the at least one receiver receives the first magnetic field returning from the formation.
42 ENGINEERING↗
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↗
Low Temperature Geothermal Play Fairway Analysis For The Appalachian Basin: Phase 1 Revised Report November 18, 2016
Geothermal energy is an attractive sustainable energy source. Yet project developers need confirmation of the resource base to warrant their time and financial resources. The Geothermal Play Fairway Analysis of the Appalachian Basin evaluated risk metrics that communicate the favorability of potential low-temperature geothermal energy resources in reservoirs more than 1000 m below the surface. This analysis is focused on the direct use of the heat, rather than on electricity production. Four risk factors of concern for direct-use geothermal plays in the Appalachian Basin portions of New York, Pennsylvania, and West Virginia are examined individually, and then in combination: 1) thermal resource quality, 2) natural reservoir quality, 3) induced seismicity, and 4) utilization opportunities. Uncertainty in the risk estimation is quantified. Based on these metrics, geothermal plays in the Appalachian Basin were identified as potentially viable for a variety of direct-use-heat applications. The methodologies developed in this project may be applied in other sedimentary basins as a foundation for low temperature (50-150 °C), direct use geothermal resource, risk, and uncertainty assessment. Three methods with which to combine the four risk factors were used. Among these, the averaging of the individual risk factors indicates the most favorable counties within the study area are the West Virginia counties of Monongalia, Harrison, Lewis (dubbed the Morgantown–Clarksburg play fairway), Putnam, and Kanawha (Charleston play fairway), the New York counties of Chemung and Steuben plus adjacent Bradford county in Pennsylvania (Corning–Ithaca play fairway), and the Pennsylvania counties of Mercer, Crawford, Erie, and Warren, and adjacent Chautauqua county in New York (together, the Meadville–Jamestown play fairway). These higher priority regions are surrounded by broader medium priority zones. Also worthy of additional exploration is a broad region near Pittsburgh Pennsylvania, for which the available geological data are insufficient to fully analyze the geological risks but yet the population is high. First, to assess the spatial variation in the depth to which one would need to drill to obtain geothermal temperatures that are useful to a future project, the project used bottom-hole temperature data from Appalachian Basin oil and gas exploration. These bottom hole temperature data are abundant but of low quality. Second, the project examined the potential for sufficient water flow rates through rocks to harvest heat from a geothermal well field, considering only natural reservoirs. This analysis provides a very incomplete picture of spatial variability of natural reservoirs because the oil and gas reservoir data lack key properties and are spatially biased toward those locations with profitable amounts of hydrocarbons in the rock pore spaces. Third, in light of the fact that earthquake activity has been induced in several states by subsurface work related to the oil and gas industry, this project examined the potential for similar activity in the Appalachian Basin. Acknowledging that data for such a task are insufficient, we utilized what was available: records of seismic activity, regional estimates of the orientations of stress in the rocks, and locations and orientations of zones of lateral change in rock properties at depths down to several kilometers below Earth’s surface. With these data, we created a first approximation of spatially variable risks for induced earthquakes. Because no data existed with which to test the reliability of these methods, the results have a high degree of uncertainty. Fourth, we examined the spatial variability of the above-the-ground factors that contribute to the economical viability of projects to tap low-temperature geothermal resources for direct-use. We worked principally with population density as a regionally known variable that would impact the cost of district heating. The resulting maps omit the costs of producing the hot water from the ground, because the below-ground costs are directly coupled to the thermal resource risk factor and natural reservoir risk factor – later analyses of those costs will be needed. The result of the district heating analysis is highly skewed: few census locations yielded a low estimated surface cost. The team also identified more than 165 prospects for high value direct-use geothermal energy opportunities such as industrial sites, university campuses, and federal facilities, among others. At the closure of this regional analysis, the most significant technical uncertainties are 1) reservoir distribution and capacities; 2) validity of thermal resource maps, and 3) the holistic estimation of Levelized Cost of Heat for favorable geological situations.
15 GEOTHERMAL ENERGY↗
Deep Borehole Field Test Research Activities at LBNL
The goal of the U.S. Department of Energy Used Fuel Disposition’s (UFD) Deep Borehole Field Test is to drill two 5 km large-diameter boreholes: a characterization borehole with a bottom-hole diameter of 8.5 inches and a field test borehole with a bottom-hole diameter of 17 inches. These boreholes will be used to demonstrate the ability to drill such holes in crystalline rocks, effectively characterize the bedrock repository system using geophysical, geochemical, and hydrological techniques, and emplace and retrieve test waste packages. These studies will be used to test the deep borehole disposal concept, which requires a hydrologically isolated environment characterized by low permeability, stable fluid density, reducing fluid chemistry conditions, and an effective borehole seal. During FY16, Lawrence Berkeley National Laboratory scientists conducted a number of research studies to support the UFD Deep Borehole Field Test effort. This work included providing supporting data for the Los Alamos National Laboratory geologic framework model for the proposed deep borehole site, conducting an analog study using an extensive suite of geoscience data and samples from a deep (2.5 km) research borehole in Sweden, conducting laboratory experiments and coupled process modeling related to borehole seals, and developing a suite of potential techniques that could be applied to the characterization and monitoring of the deep borehole environment. The results of these studies are presented in this report.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
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↗
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↗
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↗
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↗
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.
15 GEOTHERMAL ENERGY↗
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↗
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↗
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↗
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↗
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↗
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↗