SEARCH · Search NASA
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.
Not Available
This is a presentation featuring an overview of the Geoanalytical Economic Evaluation of Saline Storage (GEESS) project and latest results. The GEESS project has worked towards characterizing 57 geologic saline formations targeted for geologic carbon sequestration (GCS) using publicly available datasets. The GEESS system consists of high spatial resolution datasets (up to a 5 km grid spacing) that characterize critical geologic parameters such as depth, thickness, porosity, permeability, fracture pressure, and more. Further, GEESS geologic data were exercised using the FECM/NETL Saline Storage Cost Model (CO2_S_COM) to estimate CO2 plume sizes and the first-year break-even price of CO2 at the grid point level. GEESS is now available on NETL’s Energy Data Exchange (EDX).
Despite the long history of hydraulic fracturing, the optimal number of fracturing stages during multi-stage fracture stimulation in horizontal wells can be a complex subject. In addition to the increased expense of multistage fracturing in horizontal wells, increasing the number of fracturing stages does not always correlate with a commensurate increase in production. The problem is the application of a uniform fracture stimulation design to all stages with no consideration for geological variations along the wellbore. The result is an inefficient use and costly waste of energy and water. Optimization of the fracturing process requires an understanding of the cause-and-effect relationship between fracturing parameters and local geological properties at a given location along the wellbore. Realizing that the generalized rock mechanics theories and hypotheses are not truly applicable to fractured and laminated shales, quantifiable impacts of a shale’s geomechanical and depositional features are a prerequisite for design and implementation of optimized hydraulic fractures. The overarching goal of this project is to understand and define the relationships of shale geology and fracture dynamics using detailed field data that includes coring of the fracture domain. Finally, resource recovery from shale formations is estimated to be less than fifteen percent. Research proposed in this project will establish the foundations for investigating enhanced recovery techniques for increased resource recovery in existing fracture treated wells. Natural gas as an EOR fluid is in the initial stages of being used broadly in the oil and gas industry and as such, many aspects of the process need to be researched and addressed prior to widespread acceptance. In addition, re-fracturing techniques, particularly using cemented liners also need to be further investigated. This project explores several of these EOR techniques in greater detail, not just at the theoretical level, but through actual field test sites. This report delves into the research program, identified tasks, several studies, observations, results, and conclusions under the Hydraulic Fracturing Test Site (HFTS) project. This test project is subdivided into three distinct phases and each phase is discussed in detail under separate sections. The original test site for Phase 1 as well as the Phase 2 site are in Midland Basin in Texas and operated by Laredo Petroleum. In the first two phases, the program targeted the Wolfcamp formation in Midland Basin. Phase 3 of the study targeted the EagleFord play in South Texas with Devon Energy as the operator. Significant conclusions from this study include: 1) lateral and vertical fracture growth mapping as well as proppant transport behavior, 2) Use of field gas and its applicability for huff n puff EOR injection in the Midland Basin, 3) Development and validation of proppant log, both at core scale as well as in-fill well implementation by analyzing drilling mud return samples, 4) Applicability of liner refracturing including best practices for enhanced unit level productivity. This report will detail all these studies as well as technical results and conclusions from said studies. These include general subsurface characterization work, completion designs as implemented, studies looking at fracture geometry, high resolution microseismic study, advanced crosswell strain mapping studies including drainage characterization, core characterization results, proppant analysis (core and observation wells), geochemistry, etc. to name a few.
The ability to accurately model the subsurface transport of radionuclides is fundamental to the remote detection and characterization of underground nuclear explosion (UNE) events. Developing more sophisticated transport models presents a significant opportunity to enhance monitoring capabilities, particularly in the reliable prediction of signature migration. Experimentally determined characterization of geologic materials associated with transport properties is the pertinent base information for such robust model development and calibration. Here, we report results from an unprecedented study demonstrating changes to the pore and fracture network structures in geological materials in response to UNEs over nanometer to micrometer length scales. Volcanic tuffs of five different lithological formations from pre- and post-UNE environments were collected from the Nevada National Security Site. Combined ultra-small and small-angle neutron scattering techniques were used to characterize the tuff pore structure. The results demonstrate measurable differences in the specific surface area and porosity of samples pre- and post-shot from texturally similar lithological formations, indicating that pore properties can serve as a direct physical signature of a UNE. The results also provide experimentally determined transport parameters in support of advanced model development through the integration of gas migration, hydrodynamic simulations, and geologic framework models.
Geologic disposal is a promising solution for a safe permanent isolation of accumulated high-level nuclear waste from nuclear powerplants. The behavior of host rock is highly coupled thermally, hydromechanically, and chemically. Numerical simulations of such coupled phenomena for the extremely long term (>100,000 years) and large length scale (>1 km) of geologic disposal remain to be computationally challenging. In this study, a methodology has been developed to approximate stress and stress-induced permeability change in host rock using only thermo-hydrological (TH) variables. A coupled thermo-hydromechanical (THM) simulation is carried out using TOUGH-FLAC simulator to model THM behaviors of a generic nuclear waste repository, in order to evaluate the performance of the developed methodology, which is implemented in a coupled TH simulation. Results show that stress and permeability change estimated by the developed methodology in the TH-coupled simulation match those calculated in the THM-coupled simulation over the simulated timeframe of over 10,000 years. Details about errors in stress and permeability estimates accrued by the developed methodology are discussed here. The developed methodology will help incorporate stress-induced permeability change into existing TH simulators for the long-term radionuclide transport in geologic disposal.
This paper will present the characterization, scaling, and design of an intermediate-scale field test of a fracture thermal energy storage system (FTES). Seasonal storage of thermal energy has the potential to both significantly reduce the total energy requirements for heating and cooling of buildings, but also allow for the flexibility to store thermal energy from intermittent sources. With approximately half of global energy consumption currently being used for heating and cooling, this represents an important path to reducing greenhouse gas (GHG) emissions. The concept of FTES is to drill two or more wells into a low permeability formation, generally at a depth of less than a few hundred meters, and then generate hydraulic fractures to create flow paths for water to circulate between the wells. Hot or cold thermal energy can then be stored in the surrounding rock mass by circulating hot or cold water through the fractures, which will heat or cool the rock mass. To recover the stored energy, ambient temperature water can be then circulated through the fractures, which will then be heated or cooled by the rock mass. Fractures inherently have a very large ratio of surface area to volume. This allows for very high heat fluxes to and from the rock mass to be achieved despite the relatively low thermal conductivity of most geologic formations. Because large fractures can be made with low-cost equipment and with only inexpensive and environmentally safe materials such as water and sand, the cost to construct even large FTES systems is expected to be quite low. This paper will present what the performance requirements, size, and operating conditions of a full-scale system to operate a commercial building. The paper will describe how these full-scale system characteristics will be used as a design basis for an intermediate-scale FTES test to be conducted at the Sanford Underground Research Facility (SURF) in Lead, SD.
The inclusion of special interface elements between standard finite elements is becoming a common technique to model discontinuities in porous media using the finite element method (FEM). Recently, it was proved that high aspect ratio (HAR) elements can reproduce very satisfactorily the main features of behavior observed in a single hydraulically-induced fracture. In this paper, the mesh fragmentation technique (MFT) is proposed to extend the application of HAR elements to simulate more general problems involving the formation of hydraulic fractures in rocks by introducing this type of element in-between the standard FE of a typical mesh. The proposed framework makes use of standard FE techniques. Therefore, it is relatively easy to upgrade an existing FE program for continuous porous media to deal with evolving discontinuities. The proposed approach is verified and validated against available analytical and numerical solutions. The crack patterns predicted by the MFT are also compared with the ones observed in the laboratory. The formation of multiple fractures is also analyzed. In all cases, the proposed technique was able to properly capture the main features of rock behavior subjected to hydraulic fracturing.
Fractures are conduits for fluid flow in low-permeability geological formations. Multiphase flow properties of fractures are important in natural processes and in engineering applications such as the evaluation of the sealing capacity of caprocks and productivity of hydrocarbon-bearing tight rocks. Investigations of flow and transport through fractures typically focus on the effects of fracture geometric and mechanical factors such as aperture, roughness, and compressibility. The wettability of the fracture surfaces and its influence on microscale interfacial phenomena and macroscale effective transport properties are seldom studied. In this study, we investigated the effect of heterogeneous wetting properties on the displacement of water by supercritical CO2 through a series of lattice Boltzmann method simulations. The results show the evolution of the CO2 plume within a fracture is controlled by both the roughness of the aperture field and the wetting distribution. We combined these factors into a capillary pressure map that can be related to the macroscopic flow behavior of the fracture. We observed that heterogeneous wetting distributions promote the residual trapping of water where lower capillary pressures allowed for isolated water pockets in higher capillary pressure zones. Analysis of fracture unsteady relative permeability shows the effect of wetting on permeability evolution and provides support for the viscous-coupling relative permeability model. Finally, analysis of the steady-state relative permeability and saturation demonstrates a strong correlation between permeability and the standard deviation of the capillary pressure field. Thus, characterizing the distribution of wetting properties of fractures is crucial to understanding multiphase fracture flow and transport properties.
The efficacy of geological carbon sequestration is reliant on the integrity of the caprock and its resistance to physical and chemical alteration. Caprocks with high abundance of reactive carbonates like calcite are susceptible to acid-promoted dissolution and can result in structural weakening. This work investigates the effect of acidified brine flow through an artificially fractured, high-carbonate (30 % by XRD) shale under differential compressive stress. Cylindrical samples were cut in half vertically and milled to create an artificial fracture with interlocking asperities and open channels. Samples were sheared with a single applied stress in a custom flow cell housed within an industrial CT scanner. Further, either acidic (pH 4) or reservoir-simulated (pH 9.5) brine was flowed through the artificial fracture for 7–8 days under reservoir pressure and room temperature. Model simulations indicate flow mainly occurred in open channels, with limited flow between overlapping asperities. Analysis of fracture surfaces by optical and scanning electron microscopy show increased surface alteration and roughness after exposure to pH 4 versus pH 9.5 brine indicating mineral dissolution/loss, and this effect is greater in areas that receive the highest brine flows. Similarly, surface analysis by scratch testing shows fracture toughness decreases more after exposure to acidic versus reservoir-simulated brine, with the greatest alteration in areas of highest acidic brine flows. Despite weakening, no shear slip occurred. Overall, the results indicate that acidified brine can result in significant physical and geomechanical alteration of irregular fracture surfaces in shale caprock, with greatest effects in preferential flow regions.
The Hydraulic Fracturing Test Site 1 (HFTS-1) was a field study performed in the Wolfcamp Formation in the West Texas Permian (Midland) Basin, USA, with a focus on improving the efficiency of hydraulic fracturing. Investigating site-specific rock-fluid geochemical interactions during hydraulic fracturing is an important step to understanding the impact on formation shale porosity, permeability, and long-term shale gas production. During field operations in this region, hydraulic fracturing fluid (HFF) injection usually starts with a concentrated acid spearhead for rapid rock dissolution, followed by the injection of near-neutral pH slickwater containing chemicals and proppants. A multistep sequential injection approach was used to investigate different stages of rock-fluid interactions. The carbonate content in the host rock is important when acid spearhead is considered, as carbonate mineral dissolution is rapid and can result in porosity and permeability changes in the shale matrix. Here, in this study, we designed flow-through experiments using fractured carbonate-rich and clay-rich Wolfcamp shale cores with (1) a short-time acid soaking step and (2) a long-term slickwater flow-through step to simulate the injection method used at HFTS-1. The fluid chemistry was analyzed. A thorough mineralogical progression [e.g., Calcium (Ca) dissolution and iron (Fe) redox progression] in the cores during HFF injection was also characterized and imaged by synchrotron microprobe. Reactive transport modeling was performed based on the experimental setup. The results showed that the acid spearhead is a crucial step in creating a reaction front by mineral dissolution, especially in carbonate-rich shales. A slight layer of ferrihydrite precipitated during the slickwater flow-through period. This study provides insights into potential geochemical impact due to hydraulic fracturing operations in the Permian Basin.
Abstract Fluids with different densities often coexist in subsurface fractures and lead to variable‐density flows that control subsurface processes such as seawater intrusion, contaminant transport, and geologic carbon sequestration. In nature, fractures have dip angles relative to gravity, and density effects are maximized in vertical fractures. However, most studies on flow and transport through fractures are often limited to horizontal fractures. Here, we study the mixing and transport of variable‐density fluids in vertical fractures by combining three‐dimensional (3D) pore‐scale numerical simulations and visual laboratory experiments. Two miscible fluids with different densities are injected through two inlets at the bottom of a fracture and exit from an outlet at the top of the fracture. Laboratory experiments show the emergence of an unstable focused flow path, which we term a “runlet.” We successfully reproduce the unstable runlet using 3D numerical simulations and elucidate the underlying mechanisms triggering the runlet. Dimensionless number analysis shows that the runlet instability arises due to the Rayleigh‐Taylor instability (RTI), and flow topology analysis is applied to identify 3D vortices that are caused by the RTI. Even under laminar flow regimes, fluid inertia is shown to control the runlet instability by affecting the size and movement of vortices. Finally, we confirm the emergence of a runlet in rough‐walled fractures. Since a runlet dramatically affects fluid distribution, residence time, and mixing, the findings in this study have direct implications for the management of groundwater resources and subsurface applications.
This dataset contains black-and-white photographic images of individual core segments from well Preston 119 and photomicrographs taken of petrographic thin sections from well Clay 513. Minipermeameter measurement results of 2,260 samples for fracture and matrix permeability on core segments for well Preston 119. These results are summarized in figures, spreadsheets, and PowerPoints. Microsoft Excel will be necessary to examine the permeability data and software capable of reading a .PDF file will be needed to view the file of annotated core CT scans. Microsoft Powerpoint will be necessary to view the photomicrograph presentation.
This data set consists of Computed Tomographic (CT) data for five sheared Eau Claire Formation core samples with complex and heterogeneous lithology. Each sample was prefractured, housed in a specially adapted core holder, and sheared in incremental fashion. CT scans were taken before and after each shearing event.
Geological carbon sequestration (GCS) will play a critical role in decarbonization and in facilitating the transition to clean energy systems. Because CO 2 is highly mobile, ensuring its safe and permanent injection into subsurface geological formations involves monitoring over larger spatial domains and longer time periods than is typical for hydrocarbon reservoirs. This can benefit from simulation tools capable of modeling key CO 2 trapping mechanisms, particularly those optimized for speed and scalability on high-performance computing systems. Using isothermal versions of the SPE11B and SPE11C benchmark cases, we conduct a mesh refinement study simulating CO 2 injection into kilometer-scale rock formations at centimeter resolution with the GEOS open-source simulation framework. We focus on how mesh refinement improves the accuracy of convective mixing in both 2D and 3D simulations. The computational costs associated with achieving a converged solution highlight the need for predictive upscaling techniques. A systematic performance scaling analysis—including both central processing unit (CPU) and graphics processing unit (GPU) architectures—complements the “Results” section.
Site characterization for underground injection and storage of gigatonne-scale CO₂ requires reliable and cost-effective methods to detect and characterize faults and fractures and to assess their stress state and fault activation potential. This is critical, as wastewater injection and disposal have been shown to activate faults and induce earthquakes, and CO₂ leakage remains a key concern for long-term storage. In this project, we developed seismic methods to detect and characterize large-scale sedimentary and crystalline basement faults and associated small-scale fractures below conventional seismic imaging resolution using multicomponent (9C) surface seismic data. Machine learning was used to automatically interpret large-scale faults, providing key information for estimating the maximum magnitude of potential induced earthquakes. High-fidelity imaging was achieved by exploiting redundancy across multiple elastic wave modes, where independent images from different modes and frequencies cross-validate each other. We also used our nonlinear signal comparison (NLSC) method for ground roll removal, improving data quality in complex near-surface conditions. The methods were validated using field data acquired in central Montana. Results show that basement faults extend into the sedimentary section and that small-scale fractures are widespread above the basement. The inferred stress orientation is consistent with regional stress data, and the estimated maximum induced earthquake magnitude is small (Mw ~2.3). The developed workflow provides a practical approach for fault and fracture characterization and for assessing induced seismicity and leakage risk. It is directly applicable to CO₂ storage site selection and to other subsurface systems.
In this study, we integrate geologic and engineering data of a naturally fractured carbonate reservoir at the Kevin Dome, Montana. Well test data are correlated with core description, geochemical and lithology study to determine the flow behavior and communication within the injection test interval and to the surrounding area. Based on a dual-continuum geologic model, numerical brine injection simulations are carried out to validate the interpretation results from our well test analytical models and forecast the probability of CO 2 injection success using current reservoir properties. As a result, our well test analytical models as well as lithology/core description suggest that fluid flow may be mainly restricted to the injection interval and the assumption of radial (horizontal) flow may be appropriate. The well test models also indicate that there is potentially formation damage with a positive skin factor although prior to brine injection well tests, well stimulation through acid treatment was performed. Our numerical simulation results appear to confirm this formation damage by showing additional pressure buildup in the injection data during later test periods. To explain this, acid may have dissolved dolomite then dolomite or calcite may have been formed again further into the matrix/fracture system. Another possible explanation is mechanical clogging of the fractures due to acid dissolving dolomite and dislodging fine grains. Our work also predicts that if no additional well stimulation is performed, the project will have a lower probability of successfully injecting 1 million tons of CO 2 into the Middle Duperow formation over 4 years.
A quantitative assessment of how faults control the migration of geofluids is critical in many areas of geosciences. Here we integrated geological fieldwork, quantitative analysis of the fractures distribution and numerical modeling to build a geometrical representation of a fault zone and to characterize its hydraulic properties. Our target is a fault located in the Majella Mountain (Italy). We collected 21 scan lines across the fault profile in order to characterize its architecture. The numerical modeling of the fracture network of the damage zones and their hydraulic parameters was performed using both commercial (Move ® ) and open source software (dfnWorks and PFLOTRAN). Move ® was used to build a representative model of the fault zone using fracture spacing as a proxy, and to model the hydraulic parameters of the different fault domains. dfnWorks and PFLOTRAN were employed to infer the hydraulic parameters of the damage zones of the fault and then upscale these properties to an equivalent continuum domain, suitable for fluid flow simulations through the whole fault zone. Our findings show how even in a relatively small area it is possible to describe changes in terms of hydraulic properties of a fault zone and to build models capable to represent these variations.
In most geologic formations, the vertical gradient of the minimum horizontal stress (𝜎 h ) exceeds the hydrostatic gradient, often driving hydraulic fractures to propagate upward. However, observations indicate that the upward growth of real hydraulic fractures is less pronounced than theoretical predictions based on smoothly varying 𝜎 h fields. In fact, the layered structure of sedimentary rocks hinders fracture propagation across layers, partially explaining the limited height growth observed in practice (Zoback et al., 2022). While crystalline rocks lack the pervasive layered fabric of sedimentary formations, they still possess structural fabric formed over their geologic history, resulting in inherently "rough" in situ stress fields. Recent studies have identified another key factor: the "roughness" of in situ stress, characterized by temporal, relatively short-wavelength fluctuations superimposed on the overall stress gradient. This stress roughness leads to apparent toughness anisotropy, where the vertical toughness appears significantly larger than the horizontal toughness (P. Fu et al., 2019; Dontsov & Suarez-Rivera, 2021). Neglecting the effects of rock fabric and stress roughness can yield inaccurate predictions of hydraulic fracture geometry and growth rates at larger length scales.