Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “well casing”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

NASA Trapezoidal Wing Simulation Using Stress-w and One- and Two-Equation Turbulence Models

The Wilcox 2006 stress-omega model (also referred to as WilcoxRSM-w2006) has been implemented in the NASA Langley code CFL3D and used to study a variety of 2-D and 3-D configurations. It predicted a variety of basic cases reasonably well, including secondary flow in a supersonic rectangular duct. One- and two-equation turbulence models that employ the Boussinesq constitutive relation were unable to predict this secondary flow accurately because it is driven by normal turbulent stress differences. For the NASA trapezoidal wing at high angles of attack, the WilcoxRSM-w2006 model predicted lower maximum lift than experiment, similar to results of a two-equation model.

Rodio, J. J.↗
NASA NTRS2017

The First Fermi-LAT SNR Catalog and Cosmic Ray Implications

While supernova remnants (SNRs) are widely thought to be powerful cosmic-ray accelerators, indirect evidence comes from a small number of well-studied cases. Here we systematically determine the gamma-ray emission detected by the Fermi Large Area Telescope (LAT) from all known Galactic SNRs, disentangling them from the sea of cosmic-ray generated photons in the Galactic plane. Using LAT data we have characterized the 1-100 GeV emission in 279 regions containing SNRs, accounting for systematic uncertainties caused by source misattribution and instrumental response. We classified 30 sources as SNRs, using spatial overlap with the radio emission position. For all the remaining regions we evaluated upper limits on SNRs' emission. In the First Fermi-LAT SNR Catalog there is a study of the common characteristics of these SNRs, such as comparisons between GeV, radio and TeV quantities. We show that previously satisfactory models of SNRs' GeV emission no longer adequately describe the data. To address the question of cosmic ray (CR) origins, we also examine the SNRs' maximal CR contribution assuming the GeV emission arises solely from proton interactions. Improved breadth and quality of multiwavelength (MW) data, including distances and local densities, and more, higher resolution gamma-ray data with correspondingly improved Galactic diffuse models will strengthen this constraint.

de Palma, F.↗
NASA NTRS2016

The First Fermi-LAT SNR Catalog and Cosmic Ray Implications

While supernova remnants (SNRs) are widely thought to be powerful cosmicray accelerators, indirect evidence comes from a small number of well-studied cases. Here we systematically determine the gamma-ray emission detected by the ermi Large Area Telescope (LAT) from all known Galactic SNRs, disentangling them from the sea of cosmic-ray generated photons in the Galactic plane. Using LAT data we have characterized the 1-100 GeV emission in 279 regions containing SNRs, accounting for systematic uncertainties caused by source misattribution and instrumental response. We classified 30 sources as SNRs, using spatial overlap with the radio emission position. For all the remaining regions we evaluated upper limits on SNRs' emission. In the First Fermi-LAT SNR Catalog there is a study of the common characteristics of these SNRs, such as comparisons between GeV, radio and TeV quantities. We show that previously satisfactory models of SNRs' GeV emission no longer adequately describe the data. To address the question of cosmic ray (CR) origins, we also examine the SNRs' maximal CR contribution assuming the GeV emission arises solely from proton interactions. Improved breadth and quality of multiwavelength (MW) data, including distances and local densities, and more, higher resolution gamma-ray data with correspondingly improved Galactic diffuse models will strengthen this constraint.

de Palma, F.↗
NASA NTRS2018

Generation of Performance Model for the Aeolian Wind Tunnel (AWT) Rotor at Reduced Pressure

The NASA Jet Propulsion Laboratory (JPL) designed the Mars Helicopter (MH) in collaboration with AeroVironment Inc., NASA Ames Research Center, and NASA Langley Research Center to explore the possibility of a vertical takeoff and landing (VTOL) Unmanned Aerial Vehicle (UAV) for flight on Mars. A 40-inch-diameter Aeolian Wind Tunnel (AWT) rotor, roughly approximating the proposed MH design by JPL, was tested in forward flight at Mars atmospheric pressure at the NASA Ames Planetary Aeolian Laboratory (PAL) in support of MH research efforts. This report describes the generation of the rotor model used to correlate with that experimental effort as reported by Ament and Koning. The 40-inch-diameter rotor was 3D-scanned and transformed into an airfoil deck. The scanned rotor airfoil sections are analyzed using C81 Generator (C81Gen) to generate the sectional aerodynamic coefficients for comprehensive analyses. A mid-fidelity computational fluid dynamics (CFD) simulation using Rotorcraft CFD (RotCFD) is pursued to efficiently estimate rotor hover and forward flight performance. Simulations at two pressures, 7 mbar (approximate Martian atmospheric pressure) and 1018 mbar (1 atmosphere), are performed to gain an understanding of the performance differences and Reynolds number effects observed. Experimental 1-atmosphere thrust for single- and dual-rotor isolated hover cases correlate well with the modeled rotor. Performance results at reduced pressure (7 mbar) show a drastic decrease in lift for equivalent RPMs tested at 1 atmosphere. Although this is primarily due to pressure reduction, Reynolds number effects also contribute to this decrease, as airfoil lift and drag coefficients are affected when compared with 1-atmosphere results. Further, simulated rotor power coefficient shows drastic increases at reduced pressures, attributed to laminar boundary layer separation, as described in Koning et al. for the MH rotor analysis. PAL experimental Martian Surface Wind Tunnel (MARSWIT) results are presented in the paper by Ament and Koning. The very low Reynolds number range is currently not well understood and presents various challenges for both experimentation and simulation.

Koning, Witold J. F.↗
NASA NTRS2004

MISR Calibration and implications for low-light-level aerosol retrieval over dark water.

The authors push the limits of the vicarious calibration method’s accuracy, aiming to assess absolute, camera-to-camera, and band-toband radiometry. Patterns repeated over many well-constrained cases lend confidence to the results, at a few percent accuracy, as do additional vicarious calibration tests performed with multiplatform observations taken during the Chesapeake Lighthouse and Aircraft Measurements for Satellites (CLAMS) campaign. Conclusions are strongest in the red and green bands, but are too uncertain to accept for the near-infrared. MISR nadir-view and MODIS low-light-level absolute reflectances differ by about 4% in the blue and green bands, with MISR reporting higher values.

Clark, Dennis↗
NASA NTRS2022

Supporting Hazard Analysis for Wildfire Response Using fmdtools and MIKA

The System Wide Safety (SWS) Safety Demonstrator (SD) Series drives development of an increasingly capable In-Time Aviation Safety Management System (IASMS) focusing on humanitarian applications, starting with wildfire response (SD-1). The goals of this report are to (1) provide an early hazard analysis and mitigation evaluation of wildfire response to support these efforts and (2) provide a demonstration of capabilities of the Fault Model Design Tools (fmdtools) and Manager for Intelligent Knowledge Access (MIKA) tools. fmdtools provides a modeling, simulation, and resiliency analysis framework in which a wildfire response model, the System Modeling and Analysis of Resiliency in Scalable Traffic Management for Emergency Response Operations (SMARt-STEReO), is built. MIKA is an intelligent knowledge manager with several capabilities, including assisting in hazard analysis by extracting and analyzing hazards from historical incident reports. The following topics are covered in the report: Understanding Wildfire Hazard Dynamics. We provide a description and simulated examples of how hazards occur in the SMARt-STEReO model of wildfire response and their effect on its outcome. This provides a common mental model and focuses the analysis presented in the remainder of the report. Wildfire Hazard Identification. MIKA identifies wildfire hazards from three relevant datasets: the ICS-209-PLUS, SAFECOM, and SAFENET. Hazards are manually organized into a taxonomy and MIKA analyzes each hazard’s effects, likelihood, severity, and risk. Evaluating Mitigation Strategies. The SMARt-STEReO wildfire response model built in fmdtools evaluates a subset of identified hazards. Specifically, we simulate the effect of communications faults and equipment faults on operator safety, the effect of changing winds and flammability, and a scenario with multiple ignition points and heavy smoke. Tool Limitations and Usage Considerations. We provide a discussion of appropriate tool use cases as well as limitations and considerations for usage. The tool findings are used to synthesize recommendations for wildfire response operations, which can be captured as part of an IASMS. Key recommendations are as follows: Hazards are identified from a broad spectrum of sources including aircraft subsystems, operational sources, and ground crew operations. Highest risk operational environment hazards identified are Evacuations. The highest risk manned aerial operations hazard categorized is Jumper Operations Mishap. Ground crew hazards that are highest risk are Burns, Cargo Operations Overhead, Dehydration, Entrapment, Falling Objects, Heart Attacks, Heat Exhaustion, Inadequate Training or Certification, Vehicle Breakdown, and Vehicle Collision. Modelled containment failures arise from a mismatch between the difficulty of the firefighting scenario and the capacity (e.g., speed, effectiveness, awareness) of the response. In firefighting scenarios where containment is possible (e.g., because the fire does not spread too quickly), these mismatches can occur because of a change in environmental conditions (e.g., wind, flammability, etc) or because of planning, equipment, or communications faults. Improvements to communications increase the capacity of the firefighting response by reducing the time needed to respond to the fire. While surveillance does not increase this capacity by itself, it increases operator safety by increasing state awareness, enabling firefighters to evade approaching fires. Increasing both has a synergistic effect. In general, these performance and resilience increases generalize over fault scenarios as well as unforeseen changes to circumstances (i.e., wind, aridity, etc.). However, these improvements need to be designed so as not to make the system prone to persistent large-scale communications outages, which can reduce performance.

Hazard analysis↗
DOE OSTI2020

Moving Zero Energy Office Buildings to the Mainstream: Establishing Design Guidelines and Energy Targets: Preprint

Office buildings represent a significant portion of the commercial sector and are ideal candidates to lead the market towards achieving zero energy or zero energy ready buildings. The variety of sizes and uses for office buildings make them challenging to determine set design guidance that can be deployed across the market. A project committee for ASHRAE’s Zero Energy - Advanced Energy Design Guide (ZE-AEDG) for small and medium office buildings , comprised of members from five major national associations with market-pull in commercial buildings, created design guidance for achieving zero energy office buildings for new construction and major retrofit. This design guidance focuses on the process, providing a pathway for both office owners and design teams to achieve zero energy status. Energy use intensity targets range from 16 kBtu/ft²·yr to 36 kBtu/ft²·yr for the different climate zones were established such that office buildings can be “zero-ready,” that is, having low energy consumption that can easily be offset with renewables, preferably on-site. The design guidance was rooted in extensive energy modeling, as well as case studies of actual office buildings that have achieved high levels of energy performance and met zero energy or zero energy ready status. A unique aspect of the zero-energy guide is that it is not compared against a baseline but focuses on achieving an absolute target. The paper presents process of creating the targets, the selection of strategies to achieve the targets, and sample case studies that show how the cost-sensitive solutions are achievable for the design community. It also presents outreach strategies to deploy market achievable zero energy measures.

buildings↗
DOE OSTI2024

Status of the Ce+BAF project

Jefferson Lab is proposing to add positron beams to the 12 GeV Continuous Electron Beam Accelerator Facility (CEBAF). A team of accelerator and nuclear scientists was formed in 2018 to develop the physics case as well as a concept for the generation, production and delivery of Continuous (CW) polarized positron beams to the experimental halls, up to the full 12 GeV. A layout of the proposed concept will be shown. We will report on the ongoing efforts in the positron generation and capture, target design, beam transport and expected properties of the e+ beam on the experimental targets at 12GeV. This project is supported by the U.S. Department of Energy, Office of Science, Office of Nuclear Physics under contract DE-AC05-06OR23177

Roblin, Yves↗
DOE OSTI2024

TCF - Sustainable Well Cement for Geothermal, Thermal Recovery, and Carbon Storage wells

Primary cementing is the most important operation performed on a subterranean well. Cement, placed in the annulus between the casing and the formations, serves as a hydraulic seal preventing fluids and gas migrations, protecting steel casing from corrosion, and supporting the well structure. Poor cementing jobs can be directly responsible for wells not reaching their full capacity, casing corrosion, compromised well integrity and, in the worst-case scenarios, well collapse. Geothermal, thermal recovery, and carbon storage wells offer environments that are especially difficult for cements to survive while the required lifespan of these wells can be years. In these wells currently used cements cannot provide durable well integrity, and new solutions are necessary. The necessity of stabilizing the electric grid, increasing its flexibility, and providing energy on demand will further expand the market of durable cement solutions for applications in high-temperature underground reservoir thermal energy storage wells (HT RTES). The design challenges of special cement systems for such wells originate from chemistry limitations of currently used well cements, aggressive environments, temperature and repeated shock conditions associated with them, and very weak formations that are not uncommon in these wells. During the last 70 years the most common systems for thermal-well construction have been Ordinary Portland Cement (OPC) (absolute majority of the wells), silica-lime system, and high aluminum cement-containing formulations. The major issues of calcium-silicate hydrates, that form during the hydration of OPC resulting in hardened material, are their poor chemical resistance to acids due to the calcium interactions with acid anions followed by its eventual dissolution, even in mild acids, such as carbonic acid from CO 2 dissolution, and inadequate bonding to the steel causing serious casing corrosion problems (Sugama & Pyatina, 2019). Performance of OPC-based cement may be significantly improved with organic additives; however, those are limited by their temperature stability. An alternative to calcium-silicate cement was developed by BNL and commercialized by Halliburton as ThermaLock TM cement. Calcium-aluminate-phosphate-based chemistry of the material allowed overcoming acid-resistance problems of OPC, especially in CO 2 -rich geothermal environments. In 2012-2015, building upon this experience, BNL developed Thermal Shock Resistant Cement (TSRC) with the support of Geothermal Technology Office (GTO) of DOE. TSRC has high-temperature stable chemistry, superior properties of cement-metal casing bond and corrosion protection (4-times better than currently used OPC-based formulation), significantly outperforms common well-cements in thermal-shock tests, has self-healing properties, and like ThermaLock TM is CO 2 resistant incorporating carbonate ions into stable hydrates (Gill et al., 2012).

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Machine learning based rate optimization under geologic uncertainty

We propose a novel approach for rate optimization during a waterflood under geologic uncertainty in reservoir properties such as permeability and porosity. The traditional approach typically involves several runs of the forward simulator. This may not scale well when the optimization is to be performed at the full field-level and over multiple geologic realizations. A machine-learning (ML) based approach which is quick and scalable for rate optimization over multiple geologic realizations is proposed instead. The training data for the model is generated by running the forward simulator with randomly assigned well rates using multiple geologic realizations. A reduced order representation of the permeability heterogeneity in each of the realizations is derived using a grid connectivity transformation (GCT). This step involves finding basis functions corresponding to the different modal frequencies of the grid connectivity represented by the grid Laplacian. The projection of the heterogeneous property field along these basis functions gives the basis coefficients that form the reduced order representation. Subsequently, for each training datapoint, streamlines are traced and the minimum time of flight (TOF) representing the tracer breakthrough time at each producer is recorded. The basis coefficients and well rates are fed to a machine learning model as input and the minimum TOF at the producers forms the output of the model. This trained model can then be used along with an optimizer for computing the optimal injection rates to maximize the injection sweep efficiency. This corresponds to minimizing the variance in the minimum TOF within each well group. Different architectures of neural network are tested using 5-fold cross validation to decide the best ML model to compute the streamline time of flight. The trained model is used to perform well rate optimization over multiple realizations of geology by using a risk tolerance penalty. The optimal well rates thus obtained are compared with two cases: a) equal well rates assigned to all injectors and producers and b) well rates obtained by optimizing over a single realization without considering the uncertainty in geology. The optimal well rates are seen to offer better oil recovery and sweep efficiency than both cases.

02 PETROLEUM↗
NASA NTRSDate not supplied

Summary of Propagation Cases of the Third AIAA Sonic Boom Prediction Workshop

A summary is provided for the propagation portion of the Third AIAA Sonic Boom Prediction Workshop held January 5, 2020 in conjunction with the AIAA SciTech 2020 conference. Near-field pressure waveforms for two cases were supplied at multiple azimuthal angles. Ground signatures across the breadth of the sonic boom primary carpet and their corresponding loudness metrics were requested. Each case had some required, as well as some optional, runs. The required cases included atmospheric profiles with measured data including winds. Optional cases included a level-acceleration sonic boom focusing case as well as prediction under standard atmospheric conditions. Participants provided ground signatures along with several loudness metrics using their best sonic boom prediction practices. All the participants' submissions are compared and discussed for each case. Noise or loudness measures are calculated and detailed comparisons and statistical analyses are performed and presented. Compared to previous workshops, the selected cases in the current workshop have much lower boom loudness levels. As was observed in the past workshops, the variation in the loudness measures and spread between participants' submissions increased as the computation proceeded from under-track locations towards the lateral cut-off. Lessons learned during this workshop are discussed and recommendations are made for potential improvements and possible subsequent workshops as we collectively attempt to refine our analysis methods.

SBPW3↗
DOE OSTI2020

Shear strength and permeability of the cement-casing interface

Here, the shear strength and hydraulic permeability of the interface between well cement and casing was investigated using a triaxial direct shear apparatus. For the first time, these experiments provide measurements under controlled stress conditions with fluid flow measurements along the interface. The low cohesion (1.1 ± 1.1 MPa) and the high friction angle (43.4 ± 2.0°) indicates that the shear strength of the interface is provided by friction. This implies that the state of stress of the cement is critical to well integrity. The hydraulic aperture of the undamaged cement-steel samples was 6.8 ± 1.0 microns. Shear damage to the interface caused a decrease (-20 %) in hydraulic aperture for samples aged up to 1 month, and an increase (+300 %) for samples cured for two years. We performed numerical simulations to estimate the leakage potential from a carbon storage operation. This model predicts negligible leakage amounts (47 tonnes) in a shear-damaged well for the modeled injection of ~1.26 million tonnes of CO 2 . Thus, our measurements indicate that the cement-casing interface is not a significant leakage pathway in its intact or damaged state, and that shear-driven failure scenarios for this interface are not a significant risk to CO 2 storage security.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1988

Numerical Simulation of Subsonic and Transonic Propeller Flow

The numerical simulation of 3-D transonic flow about a system of propeller blades is investigated. In particular, it is shown that the use of helical coordinates significantly simplifies the form of the governing equation when the propeller system is assumed to be surrounded by an irrotational flow field of an inviscid fluid. The unsteady small disturbance equation, valid for lightly loaded blades and expressed in helical coordinates, is derived from the general blade-fixed potential equation, given for an arbitrary coordinate system. The use of a coordinate system which inherently adapts to the mean flow results in a disturbance equation requiring relatively few terms to accurately model the physics of the flow. Furthermore, the helical coordinate system presented here is novel in that it is periodic in the circumferential direction while, simultaneously, maintaining orthogonal properties at the mean blade locations. The periodic characteristic allows a complete cascade of blades to be treated, and the orthogonality property affords straightforward treatment of blade boundary conditions. An ADI numerical scheme is used to compute the solution of the steady flow as an asymptotic limit of an unsteady flow. As an example of the method, solutions are presented for subsonic and transonic flow about a 5 percent thick bicircular arc blade of an 8-bladed cascade. Both high and low advance ratio cases are computed and include a lifting as well as nonlifting cases. The nonlifting solutions obtained are compared to solutions from a Euler code.

Snyder, Aaron↗
NASA NTRS2013

Comparative Study on High-Order Positivity-preserving WENO Schemes

The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more di cult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium ow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

Kotov, Dmitry V.↗
NASA NTRS2013

Comparative Study on High-Order Positivity-preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order WENO schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more difficult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium flow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

WENO↗
DOE OSTI2024

Evaluating offshore legacy wells for geologic carbon storage: A case study from the Galveston and Brazos areas in the Gulf of Mexico

In this article, federal offshore waters in the Gulf of Mexico are of interest for large-scale geologic carbon storage (GCS). However, more than 80,000 offshore oil and gas wells exist in the region, which could impact the integrity of sealing intervals. In this study, we propose a screening methodology for ranking offshore legacy wells based on the challenge they may present to GCS. The methodology relies on the review of well regulatory records to 1) identify leakage pathways and assess the potential hazards that wells pose to planned GCS operations, 2) evaluate well features that impact the accessibility of wells to determine the feasibility of potential corrective actions, and 3) rank wells based on the overall challenge they may pose for GCS. We demonstrate our framework by evaluating the construction and abandonment of 156 wells across eight areas of interest (AOIs) in shallow federal waters along the Texas Gulf Coast. The majority (99.3 %) of wells considered were constructed and plugged in a manner that did not isolate prospective GCS targets in the Upper and Lower Miocene formations and may potentially require a challenging or uncertain corrective action prior to GCS. Dataset trends suggest that the observed well construction and plugging designs may be common in shallow offshore federal waters along the Texas Gulf Coast. Consequently, operators pursuing offshore GCS projects in the region may consider selecting areas that avoid challenging wells or performing robust evaluations of legacy well leakage risks to plan corrective action prior to CO 2 injection.

58 GEOSCIENCES↗
DOE OSTI2024

3D modeling of deep borehole electromagnetic measurements with energized casing source for fracture mapping at the Utah Frontier Observatory for Research in Geothermal Energy

Here, we present a 3D numerical modelling analysis evaluating the deployment of a borehole electromagnetic measurement tool to detect and image a stimulated zone at the Utah Frontier Observatory for Research in Geothermal Energy geothermal site. As the depth to the geothermal reservoir is several kilometres and the size of the stimulated zone is limited to several 100 m, surface-based controlled-source electromagnetic measurements lack the sensitivity for detecting changes in electrical resistivity caused by the stimulation. To overcome the limitation, the study evaluates the feasibility of using a three-component borehole magnetic receiver system at the Frontier Observatory for Research in Geothermal Energy site. To provide sufficient currents inside and around the enhanced geothermal reservoir, we use an injection well as an energized casing source. To efficiently simulate energizing the injection well in a realistic 3D resistivity model, we introduce a novel modelling workflow that leverages the strengths of both 3D cylindrical-mesh-based electromagnetic modelling code and 3D tetrahedral-mesh-based electromagnetic modelling code. The former is particularly well-suited for modelling hollow cylindrical objects like casings, whereas the latter excels at representing more complex 3D geological structures. In this workflow, our initial step involves computing current densities along a vertical steel-cased well using a 3D cylindrical electromagnetic modelling code. Subsequently, we distribute a series of equivalent current sources along the well's trajectory within a complex 3D resistivity model. We then discretize this model using a tetrahedral mesh and simulate the borehole electromagnetic responses excited by the casing source using a 3D finite-element electromagnetic code. This multi-step approach enables us to simulate 3D casing source electromagnetic responses within a complex 3D resistivity model, without the need for explicit discretization of the well using an excessive number of fine cells. We discuss the applicability and limitations of this proposed workflow within an electromagnetic modelling scenario where an energized well is deviated, such as at the Frontier Observatory for Research in Geothermal Energy site. Using the workflow, we demonstrate that the combined use of the energized casing source and the borehole electromagnetic receiver system offer measurable magnetic field amplitudes and sensitivity to the deep localized stimulated zone. The measurements can also distinguish between parallel-fracture anisotropic reservoirs and isotropic cases, providing valuable insights into the fracture system of the stimulated zone. Besides the magnetic field measurements, vertical electric field measurements in the open well sections are also highly sensitive to the stimulated zone and can be used as additional data for detecting and imaging the target. We can also acquire additional multiple-source data by grounding the surface electrode at various locations and repeating borehole electromagnetic measurements. This approach can increase the number of monitoring data by several factors, providing a more comprehensive dataset for analysing the deep-localized stimulated zone. The numerical analysis indicates that it is feasible to use the combination of the energized casing and downhole electromagnetic measurements in monitoring localized stimulated zone at large depths.

58 GEOSCIENCES↗