Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “deep well”

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 109 records · Page 6

In Situ X-ray Measurements to Follow the Crystallization of BaTiO 3 Thin Films during RF-Magnetron Sputter Deposition

Here, we report on adding an important dimension to the fundamental understanding of the evolution of the thin film micro structure evolution. Thin films have gained broad attention in their applications for electro-optical devices, solar-cell technology, as well storage devices. Deep insights into fundamental functionalities can be realized via studying crystallization microstructure and formation processes of polycrystalline or epitaxial thin films. Besides the fundamental aspects, it is industrially important to minimize cost which intrinsically requires lower energy consumption at increasing performance which requires new approaches to thin film growth in general. Here, we present a state of the art sputtering technique that allows for time-resolved in situ studies of such thin film growth with a special focus on the crystallization via small angle scattering and X-ray diffraction. Focusing on the crystallization of the example material of BaTiO 3 , we demonstrate how a prototypical thin film forms and how detailed all phases of the structural evolution can be identified. The technique is shaped to enable a versatile approach for understanding and ultimately controlling a broad variety of growth processes, and more over it demonstrate how to in situ investigate the influence of single high temperature sputtering parameters on the film quality. It is shown that the whole evolution from nucleation, diffusion adsorption and grain growth to the crystallization can be observed during all stages of thin film growth as well as quantitatively as qualitatively. This can be used to optimize thin-film quality, efficiency and performance.

36 MATERIALS SCIENCE↗

Delineating Variabilities of Groundwater Level Prediction Across the Agriculturally Intensive Transboundary Aquifers of South Asia

We report the groundwater depletion in South Asia’s Himalayan, transboundary Indus-Ganges-Brahmaputra-Meghna (IGBM) rivers basin is among the highest globally. Given the high irrigation demand and population, groundwater sustainability requires an improved understanding of groundwater systems for the accurate prediction of groundwater levels (GWLs). However, the prediction of groundwater system behaviors is a significant challenge since it is dominated by spatiotemporal and subsurface depth-dependent drivers. Earlier studies that address the challenges are mainly based on the short spatial and temporal extent and/or do not separate the renewable (i.e., shallow) vs nonrenewable (i.e., deeper) groundwater signals. Here, we first identified the variable importance of spatial and depth-dependent drivers on GWL in the IGBM basin. Our results indicate a greater influence of anthropogenic factors (i.e., widespread pumping and increased population) in most parts of the IGBM basin, except in the precipitation-dominated basin of the Brahmaputra. Our next purpose was to delineate a multifactorial approach for GWL prediction using the two most used machine learning models (i.e., support vector machine and feed-forward neural network) in the literature. In general, the machine learning model outputs show a good match in comparison to the GWL from the observation wells (n = 2303 distributed across India and Bangladesh) with some limitations in areas with increased groundwater irrigation. We separately compared the results from shallow (<35 m) and deep (>35 m) observation wells, emphasizing the significance of deep groundwater pumping. Our approach highlights the importance of spatiotemporal to multidepth factors in GWL prediction and can be adopted in other parts of the globe to predict GWLs.

54 ENVIRONMENTAL SCIENCES↗

Effective climate sensitivity distributions from a 1D model of global ocean and land temperature trends, 1970–2021

Abstract Current theoretically based Earth system models (ESMs) produce Effective Climate Sensitivities (EffCS) that range over a factor of three, with 80% of those models producing stronger global warming trends for 1970–2021 than do observations. To make a more observationally based estimate of EffCS, a 1D time-dependent forcing-feedback model of temperature departures from energy equilibrium is used to match measured ranges of global-average surface and sub-surface land and ocean temperature trends during 1970–2021. In response to two different radiative forcing scenarios, a full range of three model free parameters are evaluated to produce fits to a range of observed surface temperature trends (± 2σ) from four different land datasets and three ocean datasets, as well as deep-ocean temperature trends and borehole-based trend retrievals over land. Land-derived EffCS are larger than over the ocean, and EffCS is lower using the newer Shared Socioeconomic Pathways (SSP245, 1.86 °C global EffCS, ± 34% range 1.48–2.15 °C) than the older Representative Concentration Pathway forcing (RCP6, 2.49 °C global average EffCS, ± 34% range 2.04–2.87 °C). The strongest dependence of the EffCS results is on the assumed radiative forcing dataset, underscoring the role of radiative forcing uncertainty in determining the sensitivity of the climate system to increasing greenhouse gas concentrations from observations alone. The results are consistent with previous observation-based studies that concluded EffCS during the observational period is on the low end of the range produced by current ESMs.

54 ENVIRONMENTAL SCIENCES↗

Solving high-dimensional partial integral differential equations: The finite expression method

Partial integro-differential equations (PIDEs) have broad applications in the sciences, from electro-magnetism to options pricing. Here, in this paper, we introduce a new finite expression method (FEX) to solve PIDEs. This approach builds upon the original FEX and its inherent advantages with new advances: 1) A novel method of parameter grouping is proposed to reduce the number of coefficients in high-dimensional function approximation; 2) A Taylor series approximation method is implemented to significantly improve the computational efficiency and accuracy of the evaluation of the integral terms of PIDEs. The new FEX based method, denoted FEX-PG to indicate the addition of the parameter grouping (PG) step to the algorithm, provides both high accuracy and interpretable numerical solutions, with the outcome being an explicit equation that facilitates intuitive understanding of the underlying solution structures. These features are often absent in traditional methods, such as finite element methods (FEM) and finite difference methods, as well as in deep learning-based approaches. To benchmark our method against recent advances, we apply the new FEX-PG to solve benchmark PIDEs in the literature. In high-dimensional settings, FEX-PG exhibits strong and robust performance, achieving relative errors on the order of single precision machine epsilon, significantly outperforming existing approaches based on neural networks.

Combinatorial optimization↗

Prandtl number effects on extreme mixing events in forced stratified turbulence

Relatively strongly stratified turbulent flows tend to self-organise into a ‘layered anisotropic stratified turbulence’ (LAST) regime, characterised by relatively deep and well-mixed density ‘layers’ separated by relatively thin ‘interfaces’ of enhanced density gradient. Understanding the associated mixing dynamics is a central problem in geophysical fluid dynamics. It is challenging to study LAST mixing, as it is associated with Reynolds numbers $Re := UL/\nu \gg 1$ and Froude numbers $Fr :=(2{\rm \pi} U)/(L N) \ll 1$ ( $U$ and $L$ being characteristic velocity and length scales, $\nu$ the kinematic viscosity and $N$ the buoyancy frequency). Since a sufficiently large dynamic range (largely) unaffected by stratification and viscosity is required, it is also necessary for the buoyancy Reynolds number $Re_{b} := \epsilon /(\nu N^{2}) \gg 1$ , where $\epsilon$ is the (appropriately volume-averaged) turbulent kinetic energy dissipation rate. This requirement is exacerbated for oceanically relevant flows, as the Prandtl number $Pr := \nu /\kappa = {O}(10)$ in thermally stratified water (where $\kappa$ is the thermal diffusivity), thus leading (potentially) to even finer density field structures. We report here on four forced fully resolved direct numerical simulations of stratified turbulence at various Froude ( $Fr=0.5, 2$ ) and Prandtl ( $Pr=1, 7$ ) numbers forced so that $Re_{b}=50$ , with resolutions up to $30\,240 \times 30\,240 \times 3780$ . We find that, as $Pr$ increases, emergent ‘interfaces’ become finer and their contribution to bulk mixing characteristics decreases at the expense of the small-scale density structures populating the well-mixed ‘layers’. However, extreme mixing events (as quantified by significantly elevated local destruction rates of buoyancy variance $\chi _0$ ) are always preferentially found in the (statically stable) interfaces, irrespective of the value of $Pr$ .

Mechanics↗

Photoluminescence Dynamics Defined by Exciton Trapping Potential of Coupled Defect States in DNA-Functionalized Carbon Nanotubes

Chemical reactions between semiconducting single-wall carbon nanotubes (SWCNTs) and single-stranded DNA (ssDNA) achieve spatially patterned covalent functionalization sites and create coupled fluorescent quantum defects on the nanotube surface, tailoring SWCNT photophysics for applications such as single-photon emitters in quantum information technologies. The evaluation of relaxation dynamics of photoluminescence (PL) from those coupled quantum defects is essential for understanding the nanotube electronic structure and beneficial to the design of quantum light emitters. Here, we measured the PL decay for ssDNA-functionalized SWCNTs as a function of the guanine content of the ssDNA oligo that dictates the red-shifting of their PL emission peaks relative to the band-edge exciton. We then correlate the observed dependence of PL decay dynamics on energy red-shifts to the exciton potential energy landscape, which is modeled using first-principles approaches based upon the morphology of ssDNA-altered SWCNTs obtained by atomic force microscopy (AFM) imaging. In this work, our simulations illustrate that the multiple guanine defects introduced within a single ssDNA strand strongly interact to create a deep exciton trapping well, acting as a single hybrid trap. The emission decay from the distinctive trapping potential landscape is found to be biexponential for ssDNA-modified SWCNTs. We attributed the fast time component of the biexponential PL decay to the redistribution of exciton population among the lowest energy bright states and a manifold of dark states emerging from the coupling of multiple guanine defects. The long lifetime component in the biexponential decay, on the other hand, is attributed to the redistribution of exciton population among different exciton trapping sites that arise from the binding of multiple ssDNA strands along the nanotube axis. AFM measurements indicate that those trapping sites are separated on average by ~8 nm along the nanotube axis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum reactive scattering calculations for the cold and ultracold Li + LiNa → Li 2 + Na reaction

A first-principles based quantum dynamics study of the Li + LiNa(ν = 0, j = 0) → Li 2 (ν', j') + Na reaction is reported for collision energies spanning the ultracold (1 nK) to cold (1 K) regimes. A full-dimensional ab initio potential energy surface for the ground electronic state of Li 2 Na is utilized that includes an accurate treatment of the long-range interactions. The Li + LiNa reaction is barrierless and exoergic and exhibits a deep attractive potential well that supports complex formation. Thus, significant reactivity occurs even for collision temperatures approaching absolute zero. The reactive scattering calculations are based on a numerically exact time-independent quantum dynamics methodology in hyperspherical coordinates. Total and rotationally resolved rate coefficients are reported at 56 collision energies and include all contributing partial waves. Several shape resonances are observed in many of the rotationally resolved rate coefficients and a small resonance feature is also reported in the total rate coefficient near 50 mK. Furthermore, of particular interest, the angular distributions or differential cross sections are reported as a function of both the collision energy and scattering angle. Unique quantum fingerprints (bumps, channels, and ripples) are observed in the angular distributions for each product rotational state due to quantum interference and shape resonance contributions. The Li + LiNa reaction is under active experimental investigation so that these intriguing features could be verified experimentally when sufficient product state resolution becomes feasible for collision energies below 1 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An obscured quasar census with the 4MOST IR AGN survey: design, predicted properties, and scientific goals

ABSTRACT We present the 4MOST (4-metre Multi-Object Spectroscopic Telescope) infrared (IR) AGN survey, the first large-scale optical spectroscopic survey characterizing mid-infrared (MIR) selected obscured active galactic nuclei (AGNs). The survey targets $\approx 212\,000$ obscured IR AGN candidates over $\approx 10\,000 \rm \: deg^2$ down to a magnitude limit of $r_{\rm AB}=22.8 \, \rm mag$ and will be $\approx 100 \times$ larger than any existing obscured IR AGN spectroscopic sample. We select the targets using an MIR colour criterion applied to the unWISE catalogue from the WISE (Wide-field Infrared Survey Explorer) all-sky survey, and then apply a $r-W2\ge 5.9 \rm \: mag$ cut; we demonstrate that this selection will mostly identify sources obscured by $N_{\rm H}>10^{22} \rm \: cm^{-2}$. The survey complements the 4MOST X-ray survey, which will follow up $\sim 1\,\rm M$ eROSITA (extended ROentgen Survey with an Imaging Telescope Array)-selected (typically unobscured) AGN. We perform simulations to predict the quality of the spectra that we will obtain and validate our MIR–optical colour-selection method using X-ray spectral constraints and UV-to-far-IR spectral energy distribution (SED) modelling in four well-observed deep-sky fields. We find that: (1) $\approx 80-87{{\ \rm per\ cent}}$ of the WISE-selected targets are AGN down to $r_{\rm AB}=22.1-22.8 \: \rm mag$ of which $\approx 70{{\ \rm per\ cent}}$ are obscured by $N_{\rm H}>10^{22} \: \rm cm^{-2}$, and (2) $\approx 80{{\ \rm per\ cent}}$ of the 4MOST IR AGN sample will remain undetected by the deepest eROSITA observations due to extreme absorption. Our SED-fitting results show that the 4MOST IR AGN survey will primarily identify obscured AGN and quasars ($\approx 55{{\ \rm per\ cent}}$ of the sample is expected to have $L_{\rm AGN,IR}>10^{45} \rm \: erg \: s^{-1}$) residing in massive galaxies ($M_{\star }\approx 10^{10}-10^{12} \rm \: M_{\odot }$) at $z\approx 0.5-3.5$ with $\approx 33{{\ \rm per\ cent}}$ expected to be hosted by starburst galaxies.

Andonie, Carolina (ORCID:0000000255804298)↗

Comparing Mapper Graphs of Artificial Neuron Activations

The mapper graph is a popular tool from topological data analysis that provides a graphical summary of point cloud data. It has been used to study data from cancer research, sports analytics, neurosciences, and machine learning. In particular, mapper graphs have been used recently to visualize the topology of high-dimensional artificial neural activations from convolutional neural networks and large language models. However, a key question that arises from using mapper graphs across applications is how to compare mapper graphs to study their structural differences. In this paper, we introduce a distance between mapper graphs using tools from optimal transport. We demonstrate the utility of such a distance by studying the topological changes of neural activations across convolutional layers in deep learning, as well as by capturing the loss of structural information for multiscale mapper.

mapper graphs, computational topology, machine lea↗

LGAD-in-CMOS Design Studies

We describe a process flow for a CMOS sensor with intrinsic gain. The design is based on a base wafer which is prepared with closely spaced buried n and p doping layers that form an amplifying high field gain region when the device is biased. A deep n-well forms the collecting electrode. The buried n layer serves to moderate the internal field so that the CMOS transistors within the n-well remain within their operating region.

43 PARTICLE ACCELERATORS↗

Electric Drive Technologies Consortium (EDTC)/ Cost competitive, high-Performance, highly Reliable (CPR) Power Devices on 4H-SiC (Final Report)

4H-Silicon carbide (4H-SiC) is a wide bandgap semiconductor that offers superior material properties over silicon, including higher critical electric field, thermal conductivity, and electron saturation velocity. These advantages make 4H-SiC highly attractive for high-voltage, high-efficiency power electronics. However, realizing the full potential of SiC requires device technologies that are not only high-performing but also manufacturable and reliable under real-world operating conditions. This report summarizes the outcomes of a five-year R&D effort funded by the U.S. Department of Energy (DOE) under the Electric Drive Technologies Consortium (EDTC), focused on developing cost-competitive, high-performance, and highly reliable (CPR) power devices on 4H-SiC substrates. The program targeted scalable and manufacturable 1.2 kV-class SiC MOSFETs optimized for next-generation electric vehicles, renewable energy systems, and industrial power conversion. The project delivered transformative advancements in SiC power device performance and ruggedness. Particularly, Specific on-resistance (R on,sp ) was reduced by up to 37%, from ~4.0 m$\Omega \cdot$cm 2 in earlier designs to an industry-leading 2.40 m$\Omega \cdot$cm 2 , driven by optimized doping, refined JFET widths, and layout engineering. Breakdown voltages (BV) exceeded 1600 V, marking improvement over legacy baselines, and demonstrating the robustness of newly implemented junction profiles and edge terminations. Short-circuit withstand time (SCWT) saw a remarkable 4$\times$ increase, from ~2 $\mu$s to over 8 $\mu$s, achieved through the successful deployment of deep P-well structures (~1.8–2.0 $\mu$m) via channeling implantation. This innovative process breakthrough enabled precise junction formation without MeV-class implantation tools, reduced leakage under high field stress, and allowed even the shortest-channel devices (down to 0.3 $\mu$m) to achieve both high BV and excellent ruggedness—breaking the traditional trade-off between conduction efficiency and blocking capability. Several novel architectures pushed the performance envelope further. JBSFETs—featuring embedded Schottky portions—eliminated bipolar degradation and drastically reduced third-quadrant leakage, while Ladder MOSFETs introduced a clever orthogonal conduction path that achieved a 15.4% reduction in R on,sp over standard linear designs. Switching performance reached new benchmarks: short-channel devices showed a 31% reduction in total switching energy compared to 0.5 $\mu$m counterparts, while maintaining manageable gate drive requirements. Layout-optimized structures not only improved transconductance but also accelerated switching transitions, pointing to real-world benefits in converter-level efficiency. The devices also passed rigorous reliability validation. Stress-tested across TDDB, HTGB, HTRB, HVP, and burn-in, the devices screened under 30 V/10 hr and 43 V/1 s protocols consistently exhibited tighter lifetime distributions and long-term oxide robustness. These screening techniques proved effective in identifying latent defects and ensuring deployment-grade reliability. Meanwhile, advanced 3D TCAD simulations revealed and resolved electric field hotspots—particularly in HEXFET corners—where fields exceeding 4.8 MV/cm were mitigated through geometry-aware layout corrections. Overall, the results of this project demonstrate a manufacturable and scalable SiC power device platform that addresses key DOE performance targets for efficient, robust, and reliable 1.2kV 4H-SiC Power Devices. The developed technologies represent a meaningful step forward in the commercial readiness of high-voltage SiC solutions and provide a strong foundation for continued advancement in wide bandgap power electronics.

42 ENGINEERING↗

Optimizing Alabama’s CO 2 Storage in Shelby County (Project OASIS) - Deliverable Task 7: Risk Assessment

The “Optimizing Alabama’s CO 2 Storage in Shelby County: Project OASIS” CarbonSAFE Phase II Project seeks to build on regional data sets that demonstrate that the subsurface within Shelby County, Alabama has the potential to store commercial volumes of CO 2 safely, permanently, and economically. The primary target CO 2 Storage Complex is the deep Ketona Dolomite located within a 140 square mile area of the Valley and Ridge Region of Alabama. This deep saline reservoir is beneath a confining system encompassing at least 6,500 ft of shales and other low permeability sediments. Project OASIS drilled a deep stratigraphic test well to confirm the geological properties of the confining system and storage reservoir(s) within the Storage Complex. The geological data was incorporated into numerical models to establish the areal extent of the CO 2 plume and help design the storage site and its monitoring system. The project is managed by the Southern States Energy Board (SSEB), an interstate compact organization consisting of governors and state legislative leaders from sixteen southern states, Puerto Rico, and the U.S. Virgin Islands, as well as an appointee by the President of the United States. The organizational compact provides it with access to state government organizations and legislatures. The Board also maintains an Associate Members program comprised of energy resource companies, utilities, trade groups, academic R&D science and technology experts and energy consultants. Further, the SSEB staff is experienced in managing and coordinating complex energy and environmental programs, from research programs to full-scale design and demonstrations of new and innovative technologies. Project OASIS is a public-private partnership of six entities with multiple principal investigators (PIs). SSEB’s Lead PI and Co-PI are responsible for all aspects of project performance in accordance with the DOE-NETL Cooperative Agreement. SSEB has issued subgrants to Advanced Resources International, Inc., Alabama A&M University, Auburn University, Crescent Resource Innovation, and Oklahoma State University. Advanced Resources International, Inc., issued subgrants to Baker Hughes and Loudon Technical Services for field services

20 FOSSIL-FUELED POWER PLANTS↗

Modeling impacts of ice-nucleating particles from marine aerosols on mixed-phase orographic clouds during 2015 ACAPEX field campaign

Abstract. A large fraction of annual precipitation over the western United States comes from wintertime orographic clouds associated with atmospheric rivers (ARs). Transported African and Asian dust and marine aerosols from the Pacific Ocean may act as ice-nucleating particles (INPs) to affect cloud and precipitation properties over the region. Here we explored the effects of INPs from marine aerosols on orographic mixed-phase clouds and precipitation at different AR stages for an AR event observed during the 2015 ACAPEX field campaign under low dust (<0.02 cm−3) conditions. Simulations were conducted using the chemistry version of the Weather Research and Forecasting Model coupled with the spectral-bin microphysics at 1 km grid spacing, with ice nucleation connected with dust and marine aerosols. By comparing against airborne and ground-based observations, accounting for marine INP effects improves the simulation of AR-precipitation. The marine INPs enhance the formation of ice and snow, leading to less shallow warm clouds but more mixed-phase and deep clouds, as well as to a large spillover effect of precipitation after AR landfall. The responses of cloud and precipitation to marine INPs vary with the AR stages, with more significant effects before AR landfall and post-AR than after AR landfall, mainly because the moisture and temperature conditions change with the AR evolution. This work suggests weather and climate models need to consider the impacts of marine INPs since their contribution is notable under low dust conditions despite the much lower relative ice nucleation efficiency of marine INPs.

54 ENVIRONMENTAL SCIENCES↗

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University: Preprint

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ("CUBO") was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ('CUBO') was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗

Results from an Aeromagnetic Survey to Detect Steel-Cased Wells at a Marcellus Shale Well Site in Washington County, Pennsylvania

Pennsylvania has a 150-year history of oil and gas production—the longest of any state—and this enduring activity has resulted in the drilling of more than 300,000 recorded wells. However, unknown wells likely exist because innumerable wells were drilled during Pennsylvania’s intense early oil and gas history when incomplete records were kept of well locations. There is concern that early wells are likely to be ineffectively sealed because there were no laws that required plugging when the wells were abandoned. Today, many undocumented and unplugged wells are thought to be in areas of emerging shale gas and shale oil development where open wellbores can provide a pathway for undesired upward migration of fluids and gas from hydraulically fractured reservoirs. Due to this concern, Pennsylvania regulators have asked operators to locate orphaned and abandoned wells within a 1,000-ft buffer of proposed new wells. The objective of this report is to demonstrate that high-resolution aeromagnetic surveys, historic air photos, and Light Detection and Ranging (LiDAR) imagery can be rapid and effective methods to reconnoiter large, forested areas of moderate terrain for the presence of abandoned wells. These well-finding methods were evaluated at a proposed Marcellus Shale gas drilling site in Washington County, Pennsylvania, where the methods collectively located 18 confirmed wells: 15 wells were identified from aeromagnetic surveys, two wells were identified from inspection of historical air photos, and one well was identified by evaluation of state-wide LiDAR imagery. Only six wells were previously known, and their locations, as recorded in Pennsylvania’s statewide oil and gas wells database (PA/IRIS/WIS), were often too inaccurate for the wells to be found in the dense underbrush. Twelve wells identified in this study were abandoned, unmarked, and undocumented. Aeromagnetic surveys locate wells by detecting the unique magnetic signature of vertical, steel well casing, which is depicted on magnetic maps as a “bull’s eye” type anomaly that is centered directly over the well. However, when wells were drilled and found to be sub-economic, their casing was sometimes pulled and salvaged for reuse. Such wellbores provide no magnetic response and go undetected if all casing was removed. Oftentimes attempts to retrieve well casing were not 100% successful. For example, historical records for one well in the study area indicate that the well was completed in 1902 as a dry hole and that, to the extent possible, the casing was pulled for reuse. However, a section of 10-in. diameter steel casing was not recovered and remains at an unknown depth in the wellbore. This well was easily detected by the aeromagnetic survey although only deep casing remained in the well. To mitigate for the likelihood that wellbores exist where most or all casing has been removed, this study augmented aeromagnetic data with historic air photos and digital terrain models generated from LiDAR datasets—both databases are publicly available at no cost for areas within Pennsylvania. These complementary methods located three wells where the aeromagnetic anomaly, although present, was subtle and overlooked. Together, these methods determined accurate locations for six known wells within the study area and located 12 previously unknown wells. Although it is not certain that these methods successfully located all wells in the study area, the application of these methods does represent a significant improvement over relying on existing databases for well locations. For the Appendix to the report, see: https://www.netl.doe.gov/energy-analysis/details?id=b46c417a-7c9e-4d25-b810-e6248b0217f4</p>

04 OIL SHALES AND TAR SANDS↗

Geothermal Heat Recovery Complex: Large-Scale, Deep Direct-Use System in a Low-Temperature Sedimentary Basin (Final Report)

A feasibility study of using deep direct-use (DDU) geothermal energy to heat agricultural research facilities (ARFs) was conducted at the University of Illinois at Urbana-Champaign (U of IL) and its similar application to military facilities in the Illinois Basin (ILB). The geothermal energy system (GES) investigated utilizes low-temperature (30–90°C; 90–190°F) geothermal fluid (i.e., brine) from an extraction well that is part of a deep, two-well (doublet) system that extends to the bottom of the ILB. The geothermal reservoir modeled, the Mt. Simon Sandstone (MSS), is about 1,280 m (4,200 feet) deep and 457 m thick (1,500 feet) beneath the U of IL. The DDU GES surface infrastructure includes heat exchangers connected in-parallel to pipelines carrying the geothermal fluid and fresh cold and hot water. Analysis of the GES indicated that the MSS can provide a baseload of 2 MMBtu/hr to heat the ARFs by extracting 954 m³/d (6,000 barrels/day [bbl/d]) of geothermal fluid that has a temperature of 44–46 °C (111–115 °F). In addition to analyzing the levelized cost of heat (LCOH) and life cycle costs, the environmental effects of the DDU GES were evaluated, including reduced greenhouse gas (GHG) emissions and water consumption. Multiple system designs were evaluated and then ranked based on their maximum heating performance, energy efficiency, and cost recovery. This feasibility study identified the key components of the fully-integrated DDU technology that can be implemented, both technically and economically. The results and information from this study provides end-users and policy makers with guidance for additional research on the specific components of DDU technology such that its widespread use can provide an uninterruptible energy source, increase resilience from extreme weather conditions, reduce U.S. dependency on fossil fuels, and reduce greenhouse gas (GHG) emissions. The site-specific part of this study gives U of IL administrators a realistic and pragmatic assessment of the financial resources necessary to add a DDU GES in the MSS to the campus’ energy portfolio.

15 GEOTHERMAL ENERGY↗