Search NASA⌕ Search

SEARCH · Search NASA

Results for “leaky wellbore”

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.

10 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Graphene nanoplatelets reinforced cement as a solution to leaky wellbores reinforcing weak points in hydrated Portland cement with graphene nanoparticles improves mechanical and chemical durability of wellbore cements

In this work, to improve the performance of wellbore cement under subsurface environments, graphene nanoplatelets (GNPs) were added in various percentages to Class-H cement slurry. Microstructural characterization of cement slurries cured at 90 °C and 95% RH indicates that GNP modifies the microstructure of hydrated cement by reinforcing pore spaces. As a result, the mechanical properties, such as Young's modulus and axial peak stress, obtained from high temperature triaxial compression tests are significantly improved based on different percentages of added GNPs. Furthermore, the hydrated 1 × 2 inch GNP- Portland cement cores tested under simulated deep wellbore conditions of high-temperature and high pressure (HTHP), appear to have a ductile-like behavior, when compared to a typical brittle nature of Portland cement pastes. From our observations and published data on graphene resistance to fracture, GNP addition to cement is evidently enhancing the flexibility of cement. These improvements would reduce the risks associated with wellbore cement deterioration and a consequent leakage in fossil fuel production, geothermal energy production, CO 2 storage as well as long-term sealing materials in plugging and abandonment of all wellbores at the end of their service life.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Risk-based area of review estimation in overpressured reservoirs to support injection well storage facility permit requirements for CO 2 storage projects

This paper by the Energy & Environmental Research Center presents a workflow and modeling approach for delineating a risk-based area of review (AOR) to support a U.S. Environmental Protection Agency (EPA) Class VI permit for a carbon dioxide (CO 2 ) storage project. The approach combines semianalytical solutions for estimating formation fluid leakage through a hypothetical leaky wellbore with the results of numerical reservoir simulations to define the AOR. The modeling utilizes 1) semianalytical solutions from the peer-reviewed literature for formation fluid leakage through abandoned wellbores by Raven (1990) and Avci (1994), 2) a FORTRAN model compiled and described in Cihan et al. (2011, 2012) called ASLMA (Analytical Solution for Leakage in Multilayered Aquifers), and 3) a computational framework for estimating a risk-based AOR first proposed by Oldenburg et al. (2014, 2016). Therefore, the approach builds upon well-established research and underlying hydrogeological principles that have been upheld for nearly three decades. Moreover, the ASLMA model has been broadly applied to an array of storage projects. The work presented herein extends these earlier works using a custom wrapper written in the software environment, R (R Core Team, 2020), which was developed to perform multiple runs of the ASLMA model using given ranges for one or more input parameters. In addition, the current work simulates the pressure buildup within the storage reservoir in response to CO 2 injection using a compositional simulator to better accommodate the temporospatial evolution of pressure buildup within the storage reservoir that is more accurately modeled using a heterogeneous geologic model and a compositional simulator that accounts for the multiphase interactions. The workflow is demonstrated using a case study for a 180,000-metric-ton-per-year storage project located in the PCOR (Plains CO 2 Reduction) Partnership region. For the storage project evaluated here, under the scenario where the leaky wellbore is open to a saline aquifer (thief zone) between the overlying seal (cap rock) and the underground sources of drinking water (USDW), the risk-based AOR essentially collapses to the areal extent of the CO 2 plume in the storage reservoir because the pressure buildup in the storage reservoir beyond the CO 2 plume is insufficient to drive formation fluids up a hypothetical leaky wellbore into the USDW. However, even under the conservative assumption that the leaky wellbore is not open to a thief zone, beyond the areal extent of the CO 2 plume, the incremental leakage is less than 400 m 3 over 20 years, which represents ~0.0001% or less of the total volume of water contained within the USDW rock volume. As discussed in the text, the threshold criterion for defining the risk-based AOR is site-specific and should be informed by the results of the sensitivity analysis and available site characterization data. The approach outlined in this paper is designed to be protective of USDWs and, therefore, comply with the Safe Drinking Water Act requirements and provisions for the U.S. EPA Class VI Underground Injection Control (UIC) Program (Class VI Rule) and North Dakota Administrative Code Chapter 43-05-01.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Relative permeabilities for two-phase flow through wellbore cement fractures

Multiple fluids are likely to exist in fractures and flow paths associated with leaky wellbores, including liquids (e.g., crude oil) and gases (e.g., gas exsolved from liquid). These fluids occupy and move through different portions of the pore spaces within the fractures depending on many factors, including fluid properties, fracture size, and the amount of the different fluids. Upward leakage of any phase, through the fracture, can contaminate water-bearing formations, create hazardous surface conditions, and compromise the functionality of the wellbore. Early signs of wellbore leaks may be expressed by anomalous pressure behavior at surface monitoring points on cavern storage wells. These pressure anomalies are difficult to interpret, necessitating knowledge of the factors that affect the multiphase flow in fractures and porous media. These parameters are critical to modeling multiphase flow in fractures. This insight can guide further diagnosis and maximize leak remediation. Here, our study focuses on the relationship of the liquid–gas relative permeabilities for representative variable-aperture wellbore cement fracture. To obtain the relative permeability of each phase, two-phase flow tests were conducted where both fluids were flowing simultaneously through a fractured wellbore cement specimen under a range of factors, namely (1) aperture size, (2) capillary numbers, and (3) viscosity ratio. The flow experiments were conducted under a range of confining stresses and flow velocities, using nitrogen gas and silicone oils (of different viscosities) in a specially designed pressure vessel. The sum of gas and oil relative permeabilities were found to be less than one under all conditions, which indicates that the presence of one phase affects the permeability of the other phase, and vice versa. Since the gas phase flow conditions include a significant inertial flow component in addition to viscous flow, the inertial flow coefficients at different saturation states are presented. The factors affecting the relationship between the relative permeabilities are discussed in detail. A new mathematical model for estimating the relative permeability of wellbore cement fracture is presented and experimentally validated.

58 GEOSCIENCES↗
DOE OSTI2023

Enabling site-specific well leakage risk estimation during geologic carbon sequestration using a modular deep-learning-based wellbore leakage model

Geologic carbon sequestration (GCS) is a promising technology for mitigating net carbon emissions and growing climate concern by storing CO 2 in reservoirs. Oil and gas brownfields are an attractive option for CO 2 storage, but these sites have many historical wellbores from petroleum production and can be a potential leakage pathway for CO 2 or formation brine. Therefore, risk management of GCS operations requires an assessment of potential well leakage. Due to the high uncertainty of the system, stochastic approaches are ideal for quantifying the range of risk behaviors, but they must be computationally efficient in the face of complex physics. Here, we develop a new physics-centric deep learning wellbore model to predict the leakage of CO 2 and brine through leaky wellbores. Multi-physics numerical simulations were used to generate data sets, and physics-informed features were introduced. Neural networks were optimized with an automated searching algorithm. Feature analysis quantifies the impact of each feature on model prediction and confirms the role of physics-inspired parameters. The model shows high predictive performance across a wide range of geologic and injection conditions and well attributes. In conclusion, a case study illustrates how the model is applied to assess well leakage in GCS operations.

58 GEOSCIENCES↗
DOE OSTI2020

Visco-inertial gas flow through wellbore cement fractures

Understanding the nature of fluid flow through fractured wellbore cement is fundamental for evaluating the leakage potential and risk assessments of leaky wellbores. In this study, the conditions that require considering visco-inertial flow for describing the gas flow through wellbore cement fractures were investigated. Nitrogen gas flow tests were conducted on fractured cement samples under varying pressure conditions and flow rates, covering both viscous and visco-inertial flow regimes. The data substantially deviated from Darcy's law at higher flowrates and were well-fit to Forchheimer's equation for visco-inertial flow. The inertial coefficient and critical Reynolds number were expressed as a function of the hydraulic aperture. Finally, the empirical function obtained from the experiments was used as an input to numerical simulations which showed the significant role of visco-inertial flow in wellhead pressure build-up and leakage rates, and demonstrated the importance of visco-inertial flow when modeling gas flow through wellbore cement fractures.

03 NATURAL GAS↗
DOE OSTI2022

Wellbore cement fracture permeability as a function of confining stress and pore pressure

Cemented annulus fractures are a major leakage path in a wellbore system, and their permeability plays an important role in the behavior of fluid flow through a leaky wellbore. The permeability of these fractures is affected by changing conditions including the external stresses acting on the fracture and the fluid pressure within the fracture. Laboratory gas flow experiments were conducted in a triaxial cell to evaluate the permeability of a wellbore cement fracture under a wide range of confining stress and pore pressure conditions. For the first time, an effective stress law that considers the simultaneous effect of confining stress and pore pressure was defined for the wellbore cement fracture permeability. Here the results showed that the effective stress coefficient (λ) for permeability increased linearly with the Terzaghi effective stress ( -p) with an average of λ = 1 in the range of applied pressures. The relationship between the effective stress and fracture permeability was examined using two physical-based models widely used for rock fractures. The results from the experimental work were incorporated into numerical simulations to estimate the impact of effective stress on the interpreted hydraulic aperture and leakage behavior through a fractured annular cement. Accounting for effective stress-dependent permeability through the wellbore length significantly increased the leakage rate at the wellhead compared with the assumption of a constant cemented annulus permeability.

58 GEOSCIENCES↗
DOE OSTI2025

Deep-learning-enhanced assessment of wellbore barrier effectiveness in geologic storage systems with intermediate aquifers

For geologic systems where carbon dioxide (CO 2 ) is injected underground, existing wells represent potential pathways for fluid migration. Here, this study introduces a novel deep learning model to quantify the likelihood and potential magnitude of fluid migration through wellbores at sites with intermediate aquifers or thief zones between the injection units and underground drinking water sources. Synthetic datasets, generated using reservoir simulations, captured a wide range of subsurface conditions, well attributes, operational parameters, and fluid migration scenarios. Among the regression models developed to predict brine and CO 2 leakage rates and CO 2 saturations along leaky wellbores, convolutional neural network (CNN) outperformed both Light Gradient Boosting Machine and deep neural network. Additionally, a CNN-based classification model was created to predict whether brine and CO 2 would leak along a wellbore, further improving performance over regression alone. The best models were integrated into the National Risk Assessment Partnership Open-source Integrated Assessment Model for rapid, stochastic assessment of storage system containment and leakage risks. A case study demonstrated the model’s ability to simulate fluid migration through existing wells with multiple intermediate aquifers. This computationally efficient wellbore model offers value in support of site performance evaluation and risk-informed decision making by stakeholders.

CO2 leakage↗
DOE OSTI2022

Dynamic risk assessment for geologic CO 2 sequestration

At a geologic CO 2 sequestration (GCS) site, geologic uncertainty usually leads to large uncertainty in the predictions of properties that influence metrics for leakage risk assessment, such as CO 2 saturations and pressures in potentially leaky wellbores, CO 2 /brine leakage rates, and leakage consequences such as changes in drinking water quality in groundwater aquifers. The large uncertainty in these risk-related system properties and risk metrics can lead to over-conservative risk management decisions to ensure safe operations of GCS sites. The objective of this work is to develop a novel approach based on dynamic risk assessment to effectively reduce the uncertainty in the predicted risk-related system properties and risk metrics. We demonstrate our framework for dynamic risk assessment on two case studies: a 3D synthetic example and a synthetic field example based on the Rock Springs Uplift (RSU) storage site in Wyoming, USA. Results show that the U.S. National Risk Assessment Partnership’s Open Source Integrated Assessment Model (NRAP-Open-IAM) coupled with a conformance evaluation can be used to effectively quantify and reduce the uncertainty in the predictions of risk-related system properties and risk metrics in GCS.

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 OSTI2021

Evaluation of Electrical Resistivity Tomography to Monitor the Transport of Past Releases Beneath Tank Farms

Underground storage tanks at the Hanford Site, in southeastern Washington State, hold radioactive waste generated from four decades of plutonium production. The 149 single-shell tanks and the 28 double-shell tanks have all exceeded their initial design life of approximately 25 years. At least 67 tanks are assumed to have leaked in the past, resulting in radioactive releases into the vadose zone. Gamma ray logging within dry monitoring wells is currently the primary method for tracking the migration of leaked tank waste through the vadose zone. While this approach provides an accurate assessment of radioactive contamination, that information is only provided near (within ~1m) the borehole, leaving most of the vadose zone unmonitored, particularly the important region directly beneath the tank. This report describes a numerical study that investigates the feasibility and performance of time-lapse 3D electrical resistivity tomography (ERT) for long-term monitoring of a hypothetical tank waste location and migration through the vadose zone. ERT is a method of remotely imaging the bulk electrical conductivity (EC) of the subsurface, which is significantly impacted by the presence of conductive solid and liquid tank waste. The release of liquid tank wastes increases subsurface fluid conductivity and saturation over time, creating a target to use time-lapse ERT for long-term monitoring. Although the presence of metallic infrastructure can cause ERT interference, recent advancements in ERT data processing enable the deleterious effects of buried metallic infrastructure (e.g. pipes, wellbore casings, tanks) to be removed to better determine the liquid tank waste migration over time. Three hypothetical realistic scenarios were simulated in the ERT evaluation. The first two scenarios assume the same leak amount and rate (i.e., between 1/1/1951 and 12/31/1951 at the rate of 347 m 3 per year) but different leaky tanks. Scenario 1 assumes leaks under tank B-102, which is located on the edge of the B-tank farm and surrounded by a few metallic infrastructure including cased pipes/wells/tanks. Scenario 2 assumes leaks under tank B-108, which is located near the center of the B-tank farm and surrounded by larger amount of metallic infrastructure than B-102. Scenario 3 assumes the same metallic infrastructure as B-102, with a more recent contaminant leak that was simulated to have occurred between 1/1/2018 and 12/31/2023 at a rate of 1.89 m 3 per year. The leak time in Scenarios 1 and 2 corresponds to a historical overfill event in 1951 and Scenario 3 corresponds to a recent found tank leak in 2019. In each scenario, a “true” bulk EC model vs. time reflecting contaminant migration was generated. ERT data was simulated from these “true” bulk EC models and a time-lapse ERT inversion produced “imaged” bulk EC vs. time. Three electrode configurations in two, four and eight boreholes surrounding the leak tank were used in the ERT simulations in each scenario. These borehole configurations were considered logistically feasible and cost-effective for monitoring. The hypothetical ERT boreholes are assumed to have non-metallic casing. By comparing the “imaged” bulk EC with the “true” bulk EC, it was demonstrated that the three configurations of wells used (two, four, and eight wells) were able to successfully monitor the migration of tank leaks through the vadose zone, with bulk EC resolution increasing with the number of down borehole ERT arrays for all scenarios. Therefore, the use of eight boreholes to perform ERT monitoring beneath the tanks provided the best spatiotemporal information.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗