Geological Carbon Sequestration in the Washita-Fredricksburg Formation
Geological Carbon Sequestration in the Washita-Fredricksburg Formation
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Geological Carbon Sequestration in the Washita-Fredricksburg Formation
ERTS-1 images of West Aswan area were interpreted in terms of geology, drainage, and structure. Twenty-two geological units were distinguished on ERTS-1 images in West Aswan area covering geological formations and erosional levels within some formations ranging from the Precambrian to the Quaternary. Apart from the distinction of Aswan monumental granite the investigated area shows very interesting exposures of sedimentary rocks ranging from the Cretaceous to the Quaternary. Of special interest is the delineation of the iron-ore member of the Nubian Sandstone and the phosphate-bearing formation. The tracing of the geological formations from south to north and the distinction of the varied geological units within the Pliocene and Quaternary, and the discussion on the origin of tufa are of particular significance. Also, the tracing on these images of major fractures and faults intercepting Aswan Dam Reservoir and their significance on the seepage and possible future development of diversion channels from reservoir is emphasized.
The author has identified the following significant results. During preliminary analysis of the ERTS-1 RBV band-2 images it was possible to detect many geologic features and to delineate several geologic formation boundaries. Also, the spatial frequency characteristics of some different geologic formations are distinguishable from each other. At the present time there is not enough data to provide any statistical evaluation of these differences, but it is hoped that more quantitative results will be produced within the next few months. Multispectral analysis will begin when the retrospective order for the other spectral bands of the first images is received. Since many geologic formation boundaries are visible on the single band that has been analyzed so far, even greater visibility is expected with the spectral information available in the additional bands.
Grid-scale energy storage has been identified by the U.S. Department of Energy’s (DOE) Energy Storage Grand Challenge as a necessary technology to support the continued build-out of intermittent renewable energy resources required to attain a carbon-free energy future. To meet this goal, the 2018 Department of Energy Research and Innovation Act mandated the creation of a comprehensive program to accelerate the development and commercialization of next-generation energy storage technologies. One of numerous energy storage technology options is the storage of excess energy as heated geothermal brine in suitable geologic formations. This concept, known as reservoir thermal energy storage (RTES), geologic thermal energy storage (GeoTES), aquifer thermal energy storage (ATES), etc., relies on the storage of thermal energy in geologic formations for recovery and use in large-scale direct use geothermal (e.g., district heating, industrial processes, etc.) and electrical power generation applications. This thermal energy is derived from excess or waste heat from any high-temperature heat source, such as concentrated solar or from conventional thermal/nuclear generation. As such, RTES can potentially play a significant role in meeting the energy storage shortfall in the coming decades. RTES can provide energy arbitrage through both the storage and production of thermal energy stored in geologic formations for direct use applications and can serve as a source of hot fluids that can be used to generate electricity to support peak demand ramping, thus easing stress on transmission and distribution. This energy storage option has geographic benefits in that energy can be stored locally or regionally depending on the various needs/loads. RTES can also be located across an enormous geographic area, without the need for a traditional hydrothermal resource but where thermal gradients and hydrogeology allow economic exploitation of subsurface heat. The work conducted for this project includes (1) a review of lessons learned from past high-temperature RTES international projects; (2) geochemical experimental investigation and numerical simulations of potential domestic sedimentary reservoirs and (3) development of a thermo-hydrological-mechanical (THM) numerical simulation tool for optimizing formation properties and design parameters to maximize thermal energy storage performance.
Sedimentary geological formations are known to be great candidates for geological carbon sequestration. Published studies suggest the southeast of the United States contains many formations suitable for carbon storage. The Cassville 1 Stratigraphic Borehole well could act as a potential carbon reservoir for nearby energy resource facilities in Georgia, United States. Although studies have shown that porous formations are adequate for geological carbon sequestration, it is important to understand possible geochemical reactions between CO 2 and the targeted geological formation before injecting any fluids. In this study, a sandstone sample from the Cassville 1 well is being considered for geological carbon sequestration in the Conasauga Group in Northwest Georgia. Here, the collected sandstone sample, consisting of quartz, K-feldspar, micas, kaolinite, and carbonate minerals such as calcite and dolomite, has a 6% porosity. Leveraging the formation composition and porosity, a one-dimensional continuum reactive transport model was built using CrunchFlow to assess possible geochemical reactions between injected CO 2 and the geological formation. Simulation results show that the carbonate minerals, calcite and dolomite, dissolve during the injection period of 10,000 days, increasing formation porosity from 6% to as much as 30%. The rate and extent of carbonate mineral dissolution and resulting porosity increase are highly sensitive to mineral reactive surface area values. No evidence of mineral precipitation was observed, suggesting that dissolution reactions will control porosity evolution during the CO 2 injection period.
Abstract Multiscale heterogeneity and insufficient characterization data for a specific subsurface formation of interest render predictions of multi‐phase fluid flow in geologic formations highly uncertain. Quantification of the uncertainty propagation from the geomodel to the fluid‐flow response is typically done within a probabilistic framework. This task is computationally demanding due to, for example, the slow convergence of Monte Carlo simulations (MCS), especially when computing the tails of a distribution that are necessary for risk assessment and decision‐making under uncertainty. The frozen streamlines method (FROST) accelerates probabilistic predictions of immiscible two‐phase fluid flow problems; however, FROST relies on MCS to compute the travel‐time distribution, which is then used to perform the transport (phase saturation) computations. To alleviate this computational bottleneck, we replace MCS with a deterministic equation for the cumulative distribution function (CDF) of travel time. The resulting CDF‐FROST approach yields the CDF of the saturation field without resorting to sampling‐based strategies. Our numerical experiments demonstrate the high accuracy of CDF‐FROST in computing the CDFs of both saturation and travel time. For the same accuracy, it is about 5 and 10 times faster than FROST and MCS, respectively.
Underground reservoirs offer a promising option for medium- and long-term hydrogen storage that could potentially be safer and more cost-effective than surface alternatives. Currently, there are three underground hydrogen storage (UHS) facilities in the United States, which operate solely in salt caverns and store hydrogen primarily for industrial purposes. Salt formations favorable for geologic storage are geographically limited in the U.S. If a hydrogen economy is to develop, expansion of UHS to other geologic formations such as depleted gas fields and saline aquifers will be necessary to increase the availability of hydrogen storage. While research into the feasibility of storing hydrogen in these formations is ongoing, tools and methodologies are needed for characterizing the hydrogen storage potential in these reservoirs. We introduce the Hydrogen Estimator for Logistical Planning (HELP) tool from the Subsurface Hydrogen Assessment, Storage, and Technology Acceleration (SHASTA) project (https://edx.netl.doe.gov/shasta). SHASTA-HELP is a suite of web-based tools that have numerous capabilities for UHS project planning.
Monitoring microseismic activities during CO 2 injection into geologic formations is important for ensuring the safety of the storage operations. The resulting data provide insight into the response of the storage formation to CO 2 injection and can be used to infer the underlying rock flow and mechanical properties. In this paper, assimilation of microseismic data is performed for dynamic characterization of the storage formation by using a stochastic simulation model to forecast the microseismic response of a geologic formation during CO 2 injection. Two modeling approaches are adopted to predict the space-time distribution of the injection-induced microseismicity. The first model is based on pore pressure relaxation assumption, while the second model uses coupled flow and geomechanics simulation to establish the complex physical relation between the storage formation properties and the corresponding microseismic responses during CO 2 injection. The stochastic predictive models in each case are used in ensemble data assimilation frameworks to estimate rock properties from the observed microseismic data. Two data assimilation methods are considered: (i) a new ensemble-based stochastic point process filter (EnPPF) that can directly integrate discrete microseismic events, and (ii) a variant of ensemble smoother, known as the ensemble smoother with multiple data assimilation (ES-MDA), which requires continuous representation of microseismic events for assimilation. The two methods are successfully applied to a geologically realistic model of the Farnsworth Field in Texas, with complex geologic flow units and interacting fault systems.
This user’s manual guides the use of the National Energy Technology Laboratory’s (NETL) CO₂ Storage prospeCtive Resource Estimation Excel aNalysis (CO₂-SCREEN) tool, which was developed to aid users screening geologic formations for prospective CO₂ storage resources. This manual is specific to the CO₂-SCREEN 4.0 version which is based in Python. CO₂-SCREEN applies U.S. Department of Energy (DOE) methods and equations for estimating prospective CO₂ storage resources for saline formations, shale formations, and residual oil zones (ROZ). CO₂-SCREEN was developed to be substantive and user-friendly and provide a consistent method for calculating prospective CO₂ storage resources. CO₂-SCREEN uses a Java based graphical user interface for data inputs and uses Python to calculate prospective CO₂ storage resources.
Near-surface disposal of radioactive waste in shales is a promising option to safeguard the population and environment. However, natural faults intersecting these geological formations can potentially affect the long-term isolation of the repositories. This paper characterizes the physical properties and mineralogy of the internal fault core structure intersecting the Opalinus Clay formation, a host rock under investigation for nuclear waste storage at the Mont Terri Laboratory (Switzerland). We have performed porosity, density, microstructural and mineralogical measurements in different sections of the fault, including intact clays, scaly clays and fault gouge. Mercury intrusion porosimetry analysis reveal a gouge that has a pore network dominated by nanopores of less than 10 nm, yet a high-porosity (21%) and low grain density (2.62 g/cm 3 ) when compared to the intact rock (14.2%, and 2.69 g/cm 3 ). Thus, a more permeable internal fault core structure with respect to the surrounding rock is deduced. Further, we describe the OPA fault gouge as a discrete fault structure having the potential to act as a preferential, yet narrow, and localized channel for fluid-flow if compared to the surrounding rock. Since the fault gouge is limited to a millimetres-thick structure, we expect the barrier property of the geological formation is almost not affected.
The Geoanalytical Economic Evaluation of Saline Storage (GEESS) geodatabase was developed to support the United States Department of Energy (DOE) and National Energy Technology Laboratory (NETL) in their geologic carbon storage efforts by characterizing saline geologic formations present in the FECM/NETL CO2 Saline Storage Cost Model (CO2_S_COM) [1]. Using publicly available literature and data, the GEESS geodatabase characterizes 57 geologic formations across the lower-48 U.S. states in what are called Fully Integrated Geodatabases (FIGs). The FIG is a vector polygon feature containing thousands or tens of thousands of individual polygons, which each contain discrete geologic parameter values. A list of the critical geologic parameters that are characterized in the GEESS geodatabase are described in the “Processing Steps and Workflow” part of the ReadMe file, as well as the Data Catalog accompanying the GEESS geodatabase. The FIG is the basis of the GEESS geodatabase and is the direct representation of the collected geologic data. In addition to the FIG, the GEESS system contains grid files. Due to the complexity of the FIGs, grids are used to sample the geologic data so they can be exercised within CO2_S_COM. The grid files contain the geologic data sampled from the FIG, as well as estimates of “Plume Uncertainty Diameter” and “First-year Break-even Price of CO2” derived from CO2_S_COM based on the GEESS grid data.
The author has identified the following significant results. Comparisons of imagery of three seasons, late summer-fall, winter, and spring indicate that for this region fall imagery is the best for overall geologic analysis. Winter scenes with light to moderate snow cover provide excellent topographic detail owing to snow enhancement, lower sun angle, and clarity of the atmosphere. Spring imagery has considerable reduction of tonal contrast owing to the low reflecting heavy green grass cover which subdues lithologic effects; heavy snow cover in the uplands masks topography. Mapping of geologic formations is impractical in most cases. Separation into tonal units can provide some general clues on structure. A given tonal unit can include parts of several geologic formations and different stratigraphic units can have the same tonal signature. Drainage patterns and anomalies provide the most consistent clues for detecting folds, monoclines, and homoclines. Vegetation only locally reflects lithology and structure. False color infrared 9 x 9 transparencies are the most valuable single imagery. Where these can be supplemented by U-2 color infrared for more detailed work, a tremendous amount of information is available. Adequately field checking such a large area just in one scene is the major logistic problem even in a fairly well known region.
Born Scattering Inversion (BSI) systems and methods are disclosed. A BSI system may be incorporated in a well system for accessing natural gas, oil and geothermal reserves in a geologic formation beneath the surface of the Earth. The BSI system may be used to generate a three-dimensional image of a proppant-filled hydraulically-induced fracture in the geologic formation. The BSI system may include computing equipment and sensors for measuring electromagnetic fields in the vicinity of the fracture before and after the fracture is generated, adjusting the parameters of a first Born approximation model of a scattered component of the surface electromagnetic fields using the measured electromagnetic fields, and generating the image of the proppant-filled fracture using the adjusted parameters.
Born Scattering Inversion (BSI) systems and methods are disclosed. A BSI system may be incorporated in a well system for accessing natural gas, oil and geothermal reserves in a geologic formation beneath the surface of the Earth. The BSI system may be used to generate a three-dimensional image of a proppant-filled hydraulically-induced fracture in the geologic formation. The BSI system may include computing equipment and sensors for measuring electromagnetic fields in the vicinity of the fracture before and after the fracture is generated, adjusting the parameters of a first Born approximation model of a scattered component of the surface electromagnetic fields using the measured electromagnetic fields, and generating the image of the proppant-filled fracture using the adjusted parameters.
Distributed acoustic sensing (DAS) was originally intended to measure oscillatory strain at frequencies of 1 Hz or more on a fiber optic cable. Recently, measurements at much lower frequencies have opened the possibility of using DAS as a dynamic strain sensor in boreholes. A fiber optic cable mechanically coupled to a geologic formation will strain in response to hydraulic stresses in pores and fractures. A DAS interrogator can measure dynamic strain in the borehole, which can be related to fluid pressure through the mechanical compliance properties of the formation. Because DAS makes distributed measurements, it is capable of both locating hydraulically active features and quantifying the fluid pressure in the formation. We present field experiments in which a fiber optic cable was mechanically coupled to two crystalline rock boreholes. The formation was stressed hydraulically at another well using alternating injection and pumping. The DAS instrument measured oscillating strain at the location of a fracture zone known to be hydraulically active. Rock displacements of less than 1 nm were measured. Laboratory experiments confirm that displacement is measured correctly. These results suggest that fiber optic cable embedded in geologic formations may be used to map hydraulic connections in three-dimensional fracture networks. A great advantage of this approach is that strain, an indirect measure of hydraulic stress, can be measured without beforehand knowledge of flowing fractures that intersect boreholes. The technology has obvious applications in water resources, geothermal energy, CO 2 sequestration, and remediation of groundwater in fractured bedrock.
System and methods are disclosed for transmitting and receiving electromagnetic pulses through a geological formation. A preferably programmable transmitter having an all-digital portion in a preferred embodiment may be operated at frequencies below 1 MHz without loss of target resolution by transmitting and over sampling received long PN codes. A gated and stored portion of the received signal may be correlated with the PN code to determine distances of interfaces within the geological formation, such as the distance of a water interfaces from a wellbore. The received signal is oversampled preferably at rates such as five to fifty times as high as a carrier frequency. In one method of the invention, an oil well with multiple production zones may be kept in production by detecting an approaching water front in one of the production zones and shutting down that particular production zone thereby permitting the remaining production zones to continue operating.
The Karnali highway is a vital transport link and the only primary roadway that connects the remote Karnali region to the lowlands in Mid-Western Nepal. Every year there are reports of landslides blocking the road, making this area largely inaccessible. However, little effort has focused on systematically identifying landslides and landslide-prone areas along this highway. In this study, landslides were mapped with an object-based approach from very high-resolution optical satellite imagery obtained by the DigitalGlobe constellation in 2012 and PlanetScope in 2018. Landslides ranging from 10 to 30,496 sq.m were detected within a 3 km buffer along the highway. Most of the landslides were located at lower elevations (between 500–1500 m) and on steep south-facing slopes. Landslides tended to cluster closer to the highway, near drainage channels and away from faults. Landslides were also most prevalent within the Kuncha Formation geologic class, and the forested and agricultural land cover classes. A susceptibility map was then created using a logistic regression methodology to highlight patterns in landslide activity. The landslide susceptibility map showed a good prediction rate with an area under the curve (AUC) of 0.90. A total of 33% of the study arealies in high/very high susceptibility zones. The map highlighted the lower elevated areas between Bangesimal and Manma towns with the Kuncha Formation geologic class as being the most hazardous. The banks of the Karnali River, its tributaries and areas near the highway were also highly susceptible to landslides. The results highlight the potential of very high-resolution optical imagery for documenting detailed spatial information on landslide occurrence, which enables susceptibility assessment in remote and data scarce regions such as the Karnali highway.
California's total annual greenhouse gas (GHG) emissions (425.3 MtCO 2 e) in 2018 were about 6.4% of the US total (6,677 MtCO 2 e) and around 1% of global emissions. About 39% of 2018 GHG emissions in California were from the industrial and electrical sectors. Many of these emissions were from large stationary point sources and were suitable for carbon capture retrofit with subsequent storage of the captured carbon dioxide (CO 2 ) in geological formations. Previous studies of California found suitable geology and CO 2 storage resource. This study refines and furthers prior work using a three-stage screening process of oil fields, gas fields, and underground natural gas storage (UGS) sites by combining criteria from previous studies while excluding sites that pose technical risk or are located in regions with surface restrictions including sensitive habitats and dense populations. In the first stage, 129 CO 2 storage sites in California were identified using qualification criteria based upon formation properties including geological conditions and pore pressure. The second stage identified sensitive sites by applying conservative screens including seismic activity, faulting, population density, restricted lands, and sensitive habitats. During the third stage, 61 CO 2 potential storage sites were identified by subtraction of stage 2 areas from stage 1. The potential storage volume in the third stage ranged from 1.0 to 2.0 GtCO 2 . Finally, we applied a scoring system with seven parameters to rank the 61 potential sites based on subsurface technical criteria. The scored sites are classified as high priority, medium priority, and sites for future study. Prospective CO 2 storage sites with high and moderate priority were selected and linked to CO 2 sources. There are 14 prospective sites (above 20 MtCO 2 storage resource per site) with a total storage resource of 1024 MtCO 2 distributed in Northern and Southern California. Of these sites, there are 9 potential CO 2 -EOR sites and 1 depleted oil field with a total estimated CO 2 storage volume of ~800 MtCO 2 in the Southern San Joaquin and Ventura Basin. These 10 prospective sites with a storage resource greater than 20 MtCO 2 could potentially deliver more than 20 years of storage with an average injection rate of 40 MtCO 2 /year. The remaining 4 highly prospective sites are in Northern California. Additionally, study results suggest that saline formations should be re-evaluated in concert with storage in oil, gas, and natural gas storage reservoirs.