Search NASA⌕ Search

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.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

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↗
DOE OSTI2018

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↗
NASA NTRS2010

Exit Presentation: Infrared Thermography on Graphite/Epoxy

This slide presentation reports on the internship project that was accomplished during the summer of 2010. The objectives of the project were to: (1) Simulate Flash Thermography on Graphite/Epoxy Flat Bottom hole Specimen and thin void specimens, (2) Obtain Flash Thermography data on Graphite/Epoxy flat bottom hole specimens, (3) Compare experimental results with simulation results, Compare Flat Bottom Hole Simulation with Thin Void Simulation to create a graph to determine size of IR Thermography detected defects

Comeaux, Kayla↗
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↗
NASA NTRS2005

Thermographic Characterization of Impact Damage in SiC/SiC Composite Materials

SiC/SiC composite materials targeted as turbine components for next-generation aero-engines are being investigated at NASA Glenn Research Center. In order to examine damage mechanisms in these materials, SiC/SiC coupons were impacted with 1.59 mm diameter steel spheres at increasing velocities from 115 m/s to 400 m/s. Pulsed thermography, a nondestructive evaluation technique that monitors the thermal response of a sample over time, was utilized to characterize the impact damage. A thermal standard of similar material was fabricated to aid in the interpretation of the thermographic data and to provide information regarding thermography system detection capabilities in 2.4 mm thick SiC/SiC composite materials. Flat bottom holes at various depths with aspect ratios greater than 2.5 were detectable in the thermal images. In addition, the edges of holes at depths of 1.93 mm into the sample were not as resolvable as flat bottom holes closer to the surface. Finally, cooling behavior was characterized in SiC/SiC materials and used to determine impact damage depth within an 8.5% error of a known depth.

Fox, Dennis S.↗
NASA NTRS2014

Flash Infrared Thermography Contrast Data Analysis Technique

This paper provides information on an IR Contrast technique that involves extracting normalized contrast versus time evolutions from the flash thermography inspection infrared video data. The analysis calculates thermal measurement features from the contrast evolution. In addition, simulation of the contrast evolution is achieved through calibration on measured contrast evolutions from many flat-bottom holes in the subject material. The measurement features and the contrast simulation are used to evaluate flash thermography data in order to characterize delamination-like anomalies. The thermal measurement features relate to the anomaly characteristics. The contrast evolution simulation is matched to the measured contrast evolution over an anomaly to provide an assessment of the anomaly depth and width which correspond to the depth and diameter of the equivalent flat-bottom hole (EFBH) similar to that used as input to the simulation. A similar analysis, in terms of diameter and depth of an equivalent uniform gap (EUG) providing a best match with the measured contrast evolution, is also provided. An edge detection technique called the half-max is used to measure width and length of the anomaly. Results of the half-max width and the EFBH/EUG diameter are compared to evaluate the anomaly. The information provided here is geared towards explaining the IR Contrast technique. Results from a limited amount of validation data on reinforced carbon-carbon (RCC) hardware are included in this paper.

Koshti, Ajay↗
NASA NTRS1995

Thermographic imaging for high-temperature composite materials: A defect detection study

The ability of a thermographic imaging technique for detecting flat-bottom hole defects of various diameters and depths was evaluated in four composite systems (two types of ceramic matrix composites, one metal matrix composite, and one polymer matrix composite) of interest as high-temperature structural materials. The holes ranged from 1 to 13 mm in diameter and 0.1 to 2.5 mm in depth in samples approximately 2-3 mm thick. The thermographic imaging system utilized a scanning mirror optical system and infrared (IR) focusing lens in conjunction with a mercury cadmium telluride infrared detector element to obtain high resolution infrared images. High intensity flash lamps located on the same side as the infrared camera were used to heat the samples. After heating, up to 30 images were sequentially acquired at 70-150 msec intervals. Limits of detectability based on depth and diameter of the flat-bottom holes were defined for each composite material. Ultrasonic and radiographic images of the samples were obtained and compared with the thermographic images.

Roth, Don J.↗
NASA NTRS2014

Infrared Contrast Analysis Technique for Flash Thermography Nondestructive Evaluation

The paper deals with the infrared flash thermography inspection to detect and analyze delamination-like anomalies in nonmetallic materials. It provides information on an IR Contrast technique that involves extracting normalized contrast verses time evolutions from the flash thermography infrared video data. The paper provides the analytical model used in the simulation of infrared image contrast. The contrast evolution simulation is achieved through calibration on measured contrast evolutions from many flat bottom holes in the subject material. The paper also provides formulas to calculate values of the thermal measurement features from the measured contrast evolution curve. Many thermal measurement features of the contrast evolution that relate to the anomaly characteristics are calculated. The measurement features and the contrast simulation are used to evaluate flash thermography inspection data in order to characterize the delamination-like anomalies. In addition, the contrast evolution prediction is matched to the measured anomaly contrast evolution to provide an assessment of the anomaly depth and width in terms of depth and diameter of the corresponding equivalent flat-bottom hole (EFBH) or equivalent uniform gap (EUG). The paper provides anomaly edge detection technique called the half-max technique which is also used to estimate width of an indication. The EFBH/EUG and half-max width estimations are used to assess anomaly size. The paper also provides some information on the "IR Contrast" software application, half-max technique and IR Contrast feature imaging application, which are based on models provided in this paper.

Koshti, Ajay↗
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↗
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 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 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↗