Locating Orphaned Wells Using Georeferenced Historical Aerial Imagery
For presentation at the EnerGIS conference, Canonsburg, PA, September 21, 2023.
SEARCH · Engineering Papers
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.
For presentation at the EnerGIS conference, Canonsburg, PA, September 21, 2023.
The Radiation Protection Center (RPC) of the Iraqi Ministry of Environment continues to evaluate the potential health impacts associated with the Adaya Burial Site, which is located 33 kilometers (20.5 miles) southwest of Mosul. This report documents the radiological analyses of 16 groundwater samples collected from wells located in the vicinity of the Adaya Burial Site and at other sites in northern Iraq. The Adaya Burial Site is a high-risk dump site because a large volume of radioactive material and contaminated soil is located on an unsecure hillside above the village of Tall ar Ragrag. The uranium activities for the 16 water samples in northern Iraq are considered to be naturally occurring and do not indicate artificial (man-made) contamination. With one exception, the alpha spectrometry results for the 16 wells that were sampled in 2019 indicate that the water quality concerning the three uranium isotopes (Uranium-233/234, Uranium-235/236, and Uranium-238) was acceptable for potable purposes (drinking and cooking). However, Well 7 in Mosul had a Uranium-233/234 activity concentration that slightly exceeded the World Health Organization guidance level. Eight of the 16 wells are located in the villages of Tall ar Ragrag and Adaya and had naturally occurring uranium concentrations. Wells in the villages of Tall ar Ragrag and Adaya are located near the Adaya Burial Site and should be sampled on an annual schedule. The list of groundwater analytes should include metals, total uranium, isotopic uranium, gross alpha/beta, gamma spectroscopy, organic compounds, and standard water quality parameters. Our current understanding of the hydrogeologic setting in the vicinity of the Adaya Burial Site is solely based on villager's domestic wells, topographic maps, and satellite imagery. To better understand the hydrogeologic setting, a Groundwater Monitoring Program needs to be developed and should include the installation of twelve groundwater monitoring wells in the vicinity of Tall ar Ragrag and the Adaya Burial Site. Characterization of the limestone aquifer and overlying alluvium is needed. RPC should continue to support health assessments for the villagers in Tall ar Ragrag and Adaya. Collecting samples for surface water (storm water), airborne dust, vegetation, and washway sediment should be conducted on a routine basis. Human access to the Adaya Burial Site needs to be strictly limited. Livestock access on or near the burial site needs to be eliminated. The surface-water exposure pathway is likely a greater threat than the groundwater exposure pathway. Installation of a surface-water diversion or collection system is recommended in order to reduce the potential for humans and livestock to come in contact with contaminated water and sediment. To reduce exposure to villagers, groundwater treatment should be considered if elevated uranium or other contaminants are detected in drinking water. Installing water-treatment systems would likely be quicker to accomplish than remediation and excavation of the Adaya Burial Site. The known potential for human exposure to uranium and metals (such as arsenic, chromium, selenium, and strontium) at the Adaya Burial Site is serious. Additional characterization , mitigation, and remediation efforts should be given a high priority.
This dataset has results and the model associated with the publication Ciulla et al., (2024). It contains a U-Net semantic segmentation model (unet_model.h5) and associated code implemented in tensorflow 2.0 for the model training and identification of oil and gas well symbols in USGS historical topographic maps (HTMC). Given a quadrangle map (7.5 minutes), downloadable at this url: https://ngmdb.usgs.gov/topoview/, and a list of coordinates of the documented wells present in the area, the model returns the coordinates of oil and gas symbols in the HTMC maps. For reproducibility of our workflow, we provide a sample map in California and the documented well locations for the entire State of California (CalGEM_AllWells_20231128.csv) downloaded from https://www.conservation.ca.gov/calgem/maps/Pages/GISMapping2.aspx. Additionally, the locations of 1,301 potential undocumented orphaned wells identified using our deep learning framework or the counties of Los Angeles and Kern in California, and Osage and Oklahoma in Oklahoma are provided in the file found_potential_UOWs.zip. The results of the visual inspection of satellite imagery in Osage County is in the file visible_potential_UOWs.zip. The dataset also includes a custom tool to validate the detected symbols in the HTMC maps (vetting_tool.py). More details about the methodology can be found in the associated paper: Ciulla, F., Santos, A., Jordan, P., Kneafsey, T., Biraud, S.C., and Varadharajan, C. (2024) A Deep Learning Based Framework to Identify Undocumented Orphaned Oil and Gas Wells from Historical Maps: a Case Study for California and Oklahoma. Accepted for publication in Environmental Science and Technology. The geographical coordinates provided correspond to the locations of potential undocumented orphaned oil and gas wells (UOWs) extracted from historical maps. The actual presence of wells need to be confirmed with on-the-ground investigations. For your safety, do not attempt to visit or investigate these sites without appropriate safety training, proper equipment, and authorization from local authorities. Approaching these well sites without proper personal protective equipment (PPE) may pose significant health and safety risks. Oil and gas wells can emit hazardous gasses including methane, which is flammable, odorless and colorless, as well as hydrogen sulfide, which can be fatal even at low concentrations. Additionally, there may be unstable ground near the wellhead that may collapse around the wellbore. This dataset was prepared as an account of work sponsored by the United States Government. While this document is believed to contain correct information, neither the United States Government nor any agency thereof, nor the Regents of the University of California, nor any of their employees, makes any warranty, express or implied, or assumes any legal 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 its trade name, trademark, manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or the Regents of the University of California. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof or the Regents of the University of California.
The subsurface confinement of anthropogenic carbon dioxide (CO2) demands robust risk assessment methodologies to identify potential leakage pathways. Legacy wells within the Area of Review (AoR) represent one potential leakage pathway. Robust methodologies require enormous amounts of data, which are not available for many old legacy wells. This study strategically categorizes 4386 legacy wells within the AoR of a potential CO2 storage site in the Illinois basin and identifies the high-risk wells by leveraging publicly available data—reports and well logs submitted to state regulatory agencies. Wells were categorized based on their proximity to the injection well location, depth, the mechanical integrity of well barriers, and the accessibility to these wells throughout the project lifecycle. Wells posing immediate risks were identified, guiding prioritized corrective actions and monitoring plans. Out of 4386 wells, 54 have high priority for corrective action, 10 have medium priority, and the remainder are of low priority. Case study results from the Illinois basin demonstrate the effectiveness and applicability of this approach, to assess the risk associated with legacy wells within the AoR of potential CO2 storage site, strategically categorizing over 4000 such wells despite data limitations.
This report documents the results of "A Play-Based Exploration of CO₂ Storage in the Illinois Basin" (DE-FE0032366), a project funded by the U.S. Department of Energy Office of Fossil Energy and Carbon Management and conducted by the Illinois State Geological Survey (ISGS) at the University of Illinois Urbana-Champaign in partnership with Visage Energy. The project adapted play-based exploration (PBE), a systematic basin-scale evaluation methodology from the petroleum industry, to screen areas of Illinois for commercial geologic carbon storage (GCS) in Cambro-Ordovician strata. The traditional play concept was expanded to encompass three play element groups, subsurface geologic factors, surface features, and societal factors, yielding 24 play elements with defensible suitability criteria applied through a five-tier classification scheme. An integrated geospatial database was assembled from ISGS, MGSC, MRCI, NATCARB, and public data sources, supported by significant data-improvement work including correction of legacy well locations, digitization of more than 1,300 well construction records using the DOE CATALOG team's OGRRE tool, compilation of a statewide 2D seismic database, and production of a refined fault and fold geodatabase.
Both known and unmapped plugged and abandoned wells are potential leakage pathways for CO 2 from geologic carbon sequestration (GCS) sites. Although many abandoned wells have cement bridge plugs installed to prevent leakage, the seal between the cement and the inner casing wall is subject to failure. In this study, we carried out detailed T2Well simulations of cases of sudden non-Darcy flow of CO 2 and brine leakage up the gap between a cement plug and the inner steel casing wall that becomes a fully connected flow path during the post-injection period. The goal of our study was two-fold: (1) to understand the dynamics, rates, and the characteristic temporal signals associated with the onset of leakage through various gap-aperture sizes, and (2) to suggest potential monitoring strategies based on the findings. Simulation results show that the leakage of CO 2 and brine upward is transient with interesting phase interference behavior. Time-dependent oscillatory flows with varying pressure, temperature, and flow rates of CO 2 and brine show strong dependence on gap aperture. Phase-change and decompression lead to very low temperatures at the top of the well for gap apertures larger than 4 mm suggesting that remote thermal monitoring at the ground surface may be an effective way of monitoring even if well locations are not known a priori. Pressure in the well is also indicative of CO 2 leakage. The temporal patterns of changing temperature and pressure may be useful diagnostic signals for leakage detection. Finally, these transient leakage signals may provide information on the cause of leakage and/or characteristics of the flow path that could inform effective remediation design and execution approaches.
Our goal is to develop, apply and validate a holistic thermal, hydrologic, mechanical, and chemical (THMC) workflow that also includes evaluation of induced seismic slip in EGS reservoirs. We will integrate experimental and modelling approaches to reduce parameteruncertainty and better predict and mitigate seismic hazard at Utah FORGE and future EGS sites.We propose a novel approach that incorporates 3D physics-based Earthquake simulations in THMC models, herein referred to as “THMC+E” models. This capability will enable improvedengineering decisions at Utah FORGE and move EGS operations toward repeatable, robust, economically viable, and socially accepted development. Utah FORGE management and future EGS operators may employ results of THMC+E models for decision making purposes throughout the lifetime of the field operation. For example, before production well installation, our THMC+E models will predict circulation scenarios and related seismic hazard for a suite ofpossible well locations and flow rates, thus enabling evaluation of optimal production well placement. Such efforts will be conducted throughout the lifetime of the project, whereby additional laboratory experiments will constrain key model parameters and machine learning (ML) will reduce the size of the parameter space and the associated uncertainty. THMC+E simulations will enable exploration various circumstances that may hinder EGS success and develop mitigation strategies.
In this work, we investigate by means of numerical simulation a planned year-long field test of depressurization-induced production from a permafrost-associated hydrate reservoir on the Alaska North Slope at the site of the recently drilled Hydrate-01 Stratigraphic Test Well. The main objective of this study is to assess quantitatively the impact of temporary interruptions (well shut-ins) on the expected fluid production performance from the B1 Sand of the stratigraphic Unit B during controlled depressurization over different time scales, as well as on other relevant aspects of the system response that have the potential to significantly affect the design of the field test. We consider eight different cases of depressurization, including (a) rapid depressurization over a 60-day period to a terminal bottomhole pressure P W of 2.8 MPa and (b) a slower depressurization rate to a final P W of 0.6 MPa at the end of the year-long production test, in addition to (c) a multi-step depressurization regime and (d) a quasi-linear continuous depressurization strategy. The results of the study indicate that shut-ins obviously reduce gas release and production during and immediately after their occurrence, but their longer-term effects are strongly dependent on the depressurization regime and on the time of observation, covering the entire range of potential outcomes. Shut-ins (a) have a universally strong negative effect when quasi-linear depressurization is involved regardless of the length of the production period, and (b) have a strong positive effect in multi-step depressurization schemes that becomes apparent earlier for large initial pressure drops, but (c) can also appear to have practically no effect for slow stepwise depressurization at the end of the year-long production test. Shut-ins lead to a rapid reformation of hydrates, even to the point of disappearance of a free gas phase in the reservoir. Rapid depressurization regimes lead to early maximum rates of hydrate dissociation and gas production, while the maximum rates occur at the end of the production test for the cases of slower depressurization. Shut-ins do not appear to have a significant impact on water production, as the cessation of production is followed by higher rates production when depressurization resumes. Similarly, (a) the fraction of produced CH 4 originating from exsolution from the water, (b) the water-to-gas ratio, and (c) the rate of replenishment of produced water by boundary inflows do not appear significantly affected by shut-ins, the effects of which seem to be temporary in the majority of the cases. The study confirmed the superiority of multi-step depressurization methods as the most effective strategies for hydrate dissociation and gas production and showed that two observation wells (located at distances of 30 and 50 m from the production well) are appropriately positioned and both able to capture the P, T, and S G behavior during the fluid production and shut-ins in any of the eight cases we investigated.
This data package contains pumping data (.txt), parameter matrices, and R code (.R, .RData) to perform bootstrapping for parameter selection for the bioclogging model development. The pumping data were collected from the Russian River Riverbank Filtration site located in Sonoma County, California from 2010-2017 from three riverbank collection wells located alongside the study site. The pumping data is directly correlated with water table oscillations, so the code performs these correlations and simulates stochastic versions of water table oscillations. See Metadata Description.pdf for full details on dataset production. This dataset must be used with the R programming language. This dataset and R code is associated with the publication "Influence of Hydrological Perturbations and Riverbed Sediment Characteristics on Hyporheic Zone Respiration of CO2 and N-2"This research was supported by the Jane Lewis Fellowship from the University of California, Berkeley, the Sonoma County Water Agency (SCWA), the Roy G. Post Foundation Scholarship, the U.S. Department of Energy, Office of Science Graduate Student Research (SCGSR) Program, U.S. Department of Energy, Office of Science, Office of Biological and Environmental Research under award DE-AC02-05CH11231, and the UFZ-Helmholtz Centre for Environmental Research, Leipzig, Germany.
A new web-based toolset is being developed to support ongoing remediation optimization efforts and implementation of an adaptive site management strategy for the 200 West Area Pump-and-Treat (P&T) system at the Hanford Site. This toolset, comprising the well performance index tool and the well optimization pre-screening tool, will offer a user-friendly interface to predict and optimize the P&T well network’s performance at a preliminary level. Efforts in fiscal year (FY) 2023 focused on three main components: updating the existing deep learning model for predicting P&T performance, designing and developing a prototype of a web-based performance index tool, and initiating the conceptual design of the well optimization pre-screening tool. The well performance index tool is based on a pre-trained deep learning model that allows users to select a target contaminant and well screen length, then visualize the predicted performance of potential new wells across the site. The well optimization pre-screening tool includes two separate modules: the pre-computed scenario viewer, which organizes and visualizes offline optimization simulation results, and the quick analysis module, which provides real-time model prediction using user-specified well locations. In FY24, the plan is to add web-based applications to SOCRATES for both the well performance prediction tool and the optimization prescreening tool, with accompanying user and theory guides. These tools are intended to enable an accessible, easily applied, and transparent approach to remedy planning and decision-making.
This paper develops a formalism for optimizing nozzle location/configuration with respect to combustion stability of high-frequency transverse modes in a can combustor. The stability of these acoustically noncompact flames was assessed using the Rayleigh integral (RI). Several key control parameters influence RI—flame angle, swirling strength, nozzle location, as well as nozzle location with respect to the acoustic mode shape. In this study, we consider a N-around-1 configuration such as typically used in a multi-nozzle can system and study the overall stability of this system for different natural transverse modes. Typically, such nozzles are distributed in a uniformly circular manner for which we study the overall RI, and for cases where RI > 0, we optimize the nozzle distribution that can reduce and minimize RI. For a fixed geometry such as a circular configuration, the analysis shows how the flame's parameters must vary across the different nozzles, to result in a relatively stable system. Additionally, for a fixed set of flame parameters, the analysis also indicates the noncircular distribution of the N nozzles that minimizes RI. Overall, the analysis aims to provide insights on designing nozzle locations around the center nozzle for minimal amplification of a given transverse mode.
This Storage Field Development plan presents the Storage Complex characterization results, construction, monitoring, and operational plans, and costs associated with the proposed One Earth Sequestration Carbon Capture and Storage (OES-CCS) site in McLean County, Illinois, near Gibson City. The proposed storage complex, known as the Mt. Simon Storage Complex, comprises the Cambrian Mt. Simon Sandstone reservoir and the primary seal, the Cambrian Eau Claire Formation. The lowermost Underground Source of Drinking Water (USDW) identified for the site is the Ordovician St. Peter Sandstone. Geologic characterization of the Mt. Simon Storage Complex at the OES-CCS site was performed by the Illinois Storage Corridor CarbonSAFE Phase III project, which also prepared and submitted three UIC Class VI applications to construct three injection wells; the permit applications were submitted and are in the federal EPA review process. A characterization well, OEE #1, was drilled to collect site-specific data. These data were analyzed and used to develop the UIC Class VI applications. The OEE #1 well will be converted to an in-zone monitoring (IZM) well for the injection phase. The proposed buildout for the OES-CCS site includes (1) three injection wells (OES #1, OES #2, and OES #3), (2) two IZM wells, (3) two above confining zone (ACZ) monitoring wells, one of which will be used to monitor the lowermost USDW, (4) capture and compression facilities, and (5) transportation facilities, i.e., pipelines. A pre-operational testing program was proposed in the Class VI permit application and will be employed at the site pending approval. Additional pre-injection (baseline), syn-injection, and post-injection monitoring and site care procedures will be followed by OES to ensure that injection activities are protective of human health and the environment. Injection is scheduled to begin in 2025, distributed across the three injection wells in accordance with the permit operating conditions. One Earth Sequestration intends to inject up to 90 million tonnes of CO 2 over a period of approximately 20 years. Injection will begin at approximately 0.5 million tonnes of CO 2 annually and ramp up to a maximum of 4.5 million tonnes annually. Daily injection rates are expected to range from 1,400 to 1,500 tonnes per day initially and reach a maximum of approximately 4,225 tonnes per day, depending on site geology and injectivity at each injection well location, and CO 2 availability. The costs associated with the OES-CCS project include pre-operational costs (e. g. additional seismic data acquisition and well drilling), capture and transportation facility and equipment costs, predicted field operating expenditures (OpEx), and decommissioning and post-injection site care (PISC) costs. The risks associated with project activities, such as site construction, injection operations, and verification of secure storage were evaluated, and mitigation strategies proposed to alleviate those risks.
In cases of complex but only partially known geology and a lack of spatial control in observation well locations, water table elevation estimation is very challenging. In some cases, auxiliary information, such as observations of the movement of tracers, operation of injection and extraction wells, and calibration of groundwater models against the historical elevation data, can be combined with expert judgement to estimate flow directions in areas of sparse data and to aid in the production of more reliable contour maps (and associated rasters) than could be produced by relying on sparse well elevation data alone. Given a historical sequence of these raster maps, the question arises how to automate, to the extent possible, the process of producing new raster maps to reflect data from previous times, the current data and the operation of expert judgement. One solution is to adopt a Bayesian point of view and to regard the historical well elevation data, auxiliary information and historical raster maps as prior information. The well elevations for water table wells, as well as those for injection/extraction wells and the data associated with other relevant variables, can be viewed as predictors for the raster surface. From this prior information, we can, conditional on the values of these predictors for a new time period, compute an expected value map and a standard deviation map for the new raster. These then can be taken to specify a prior predictive distribution for the pixels in the new raster map. Then we condition the pixels, corresponding to water level observation wells within the raster, on the observed values in those wells (which in general will differ from the regression estimate) for the new time period. Given the smoothness of the water table surface, we then smooth the surface of deviations from the mean surface, based on the variograms of the historical rasters, and add this smoothed surface to the regression mean surface. The error structure of the produced raster map is defined by the regression error structure and the error due to smoothing based on the estimated variograms. This methodology has been developed and is being further refined for groundwater monitoring and remediation at LANL. It is a very flexible method that can also be applied with a variety of other predictors applied to model the water level wells in the area of interest over the historical record. The smoothness of the spatial process and its possible evolution over time can then be estimated from the residuals from this regression. This can be augmented by expert hydrogeological opinion based on site topography and hydrogeology. (authors)
Catalysis from single active sites is analyzed using methods developed from single molecule kinetics. Using a stochastic Markov state description, the observable properties of general catalytic networks of reactions are expressed using an eigenvalue decomposition of the transition matrix for the Markov process. By the use of a sensitivity analysis, the necessary eigenvalues and eigenvectors are related to the energies of controlling barriers and wells located along the reaction routes. A generalization of the energetic span theory allows the eigenvalues to be computed from several activation energies corresponding to distinct barrier-well pairings. The formalism is demonstrated accurately for model problems for a physically realistic mechanism for an alkene hydrogenation reaction on a single atom catalyst.
Well integrity is critically important to maintain to minimize the environmental impacts of oil and gas development and other subsurface energy operations. The Wattenberg Field of Colorado—a top producing field with >40,000 wells—has one of the most robust publicly reported well integrity programs in the country. Here, in this study, we analyzed annular pressure and annular-fluid geochemical test results collected from Wattenberg wells through the end of 2019 to characterize the frequency and spatial variability of integrity issues in the field and understand their drivers. Estimated frequencies of integrity issues among tested wells were 8.2-17.1% between 1955 and 2019 and 6.1-11.4% in 2019 alone. The frequency of integrity issues was nearly four times greater in wells located above the Longmont Wrench Fault Zone. Potential drivers of integrity issues were identified using ensemble decision tree models trained with a broad set of relevant information. Models show that well integrity issues are spatially clustered on regional and sub-regional scales and suggest the relatively high frequency of integrity issues observed is likely attributed to geologic factors. These findings are valuable for regulatory agencies and operators seeking to inform well integrity monitoring, plugging, and emissions reduction efforts and design future subsurface energy projects.
In this document, disclosure data from 1,228 well stimulations occurring May 2015 to October 2019 in California are evaluated. This evaluation updates a previous study that was based on 618 well stimulations occurring May 2015 to June 2016 (Stringfellow, W.T., Camarillo, M.K., and Jordan, P. 2017, Status of Well Stimulation in California Since Implementation of SB-4 Regulations, Berkeley National Laboratory, Berkeley, CA). While the goal of the previous study was to provide a summary of information obtained since passage of California Senate Bill No. 4—Oil and Gas, the goal of the current evaluation is to provide a basis for studying chemical indicators that could potentially be used to detect aquifer contamination. We evaluate chemical indicators is a separate document using the data contained herein as well as data from monitoring wells located within oil fields. Compared with the data used by Stringfellow et al. 2017, the current evaluation is based on a more diverse data set. While the previous study contained data for only two oil fields with more than 30 well stimulations, the current data set contains four oil fields with more than 30 well stimulations, making these data more ideal for comparisons using statistical tests. More producers and a new oil field, Buena Vista Nose, are represented in the current data set. Some well stimulation practices remain relatively unchanged since the previous study. Well stimulation is mostly occurring in Kern County with the exception of a single well stimulation in Orange County. Almost all well stimulations are hydraulic fracturing. The current data set contains two acid fracturing treatments and one matrix acidizing treatment. The median number of chemicals added per well stimulation—excluding water and proppant—was previously 21; the number in the expanded data set is 20. Median water added per well stimulation was previously 89,000 gallons and is now 98,000 gallons. While this water use represents a 10% increase, water use for hydraulic fracturing in California is still lower than water use in other oil and gas fields. Chemical formulations appear to have shifted since the previous study. This shift may be related to a change in the predominant service company. We identified 26 chemicals for which the frequency of use has changed by more than 30%. The total number of chemicals used has also expanded. Previously, 178 unique chemicals were identified as being added to well stimulation fluids. The total number of unique chemicals is now 205. Analytical data available for base and recovered fluids has expanded. The number of base fluid samples has increased from 12 to 35 and the number of recovered fluid samples has increased from 1,078 to 2,166. The expanded data set is more ideal for an evaluation of indicator chemicals and comparisons across different formations. The analytical data for recovered fluid samples indicates that many ions, radioactive constituents, and organics are consistently observed, making these chemicals ideal for consideration as potential indicators of aquifer contamination.
This study aims to develop a methodology for calibrating subsurface stress changes through time-lapse vertical seismic profiling (VSP) integration. The selected study site is a region around the injector well located within Farnsworth field unit (FWU), where there is an ongoing CO 2 -enhanced oil recovery (EOR) operation. In our study, a site-specific rock physics model was created from extensive geological, geophysical, and geomechanical characterization through 3D seismic data, well logs, and core assessed as part of the 1D MEM conducted on the characterization well within the study area. Here, the Biot-Gassmann workflow was utilized to combine the rock physics and reservoir simulation outputs to determine the seismic velocity change due to fluid substitution. Modeled seismic velocities attributed to mean effective stress were determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from ultrasonic seismic velocity measurements. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer (PSO) was utilized to minimize a penalty function created between the modeled seismic velocities and the observed time-lapse VSP dataset. The successful execution of this workflow has affirmed the suitability of acoustic time-lapse measurements for 4D-VSP geomechanical stress calibration pending measurable stress sensitivities within the anticipated effective stress changes and the availability of suitable and reliable datasets for petroelastic modeling. © 2023 Society of Chemical Industry and John Wiley & Sons, Ltd.
Carbon capture, utilization, and storage (CCUS) describes a set of technically viable processes to separate carbon dioxide (CO 2 ) from industrial byproduct streams and inject it into deep geologic formations for long-term storage. Legacy wells located within the spatial domain of new injection and production activities represent potential pathways for fluids (i.e., CO 2 and aqueous phase) to leak through compromised components (e.g., through fractures or micro-annulus pathways). The finite element (FE) method is a well-established numerical approach to simulate the coupling between multi-phase fluid flow and solid phase deformation interactions that occur in a compromised well system. We assumed the spatial domain consists of a three-phases system: a solid, liquid, and gas phase. For flow in the two fluids phases, we considered two sets of primary variables: the first considering capillary pressure and gas pressure (PP) scheme, and the second considering liquid pressure and gas saturation (PS) scheme. Fluid phases were coupled with the solid phase using the full coupling (i.e., monolithic coupling) and iterative coupling (i.e., sequential coupling) approaches. The challenge of achieving numerical stability in the coupled formulation in heterogeneous media was addressed using the mass lumping and the upwinding techniques. Numerical results were compared with three benchmark problems to assess the performance of coupled FE solutions: 1D Terzaghi’s consolidation, Liakopoulos experiments, and the Kueper and Frind experiments. We found good agreement between our results and the three benchmark problems. For the Kueper and Frind test, the PP scheme successfully captured the observed experimental response of the non-aqueous phase infiltration, in contrast to the PS scheme. These exercises demonstrate the importance of fluid phase primary variable selection for heterogeneous porous media. We then applied the developed model to the hypothetical case of leakage along a compromised well representing a heterogeneous media. Considering the mass lumping and the upwinding techniques, both the monotonic and the sequential coupling provided identical results, but mass lumping was needed to avoid numerical instabilities in the sequential coupling. Additionally, in the monolithic coupling, the magnitude of primary variables in the coupled solution without mass lumping and the upwinding is higher, which is essential for the risk-based analyses.