Results for “Geologic Map” — Page 10

DOE OSTI2020

Height Above Nearest Drainage (HAND) and Hydraulic Property Table for CONUS

The continental flood inundation mapping (CFIM) framework is a high-performance computing (HPC)-based computational framework for the Height Above Nearest Drainage (HAND)-based inundation mapping methodology. Using the 10m Digital Elevation Model (DEM) data produced by U.S. Geological Survey (USGS) 3DEP (the 3-D Elevation Program) and the NHDPlus hydrography dataset produced by USGS and the U.S. Environmental Protection Agency (EPA), a hydrological terrain raster called HAND is computed for HUC6 units in the conterminous U.S. (CONUS). The value of each raster cell in HAND is an approximation of the relative elevation between the cell and its nearest water stream. Derived from HAND, a hydraulic property table is established to calculate river geometry properties for each of the 2.7 million river reaches covered by NHDPlus (5.5 million kilometers in total length). This table is a lookup table for water depth given an input stream flow value. Such lookup is available between water depth 0m and 25m at 1-foot interval. The flood inundation map is then computed by using HAND and this lookup table based on the near real-time water forecast from the National Water Model (NWM) at the National Oceanic and Atmospheric Administration (NOAA). HAND and the Hydraulic Property Table version 0.2. is created to correspond to data updates in the USGS 1/3 arcsec DEM, the USGS National Hydrography Dataset, including its Water Boundary Dataset, and the NHDPlus medium resolution dataset. This dataset comprises 331 HUC6 units for CONUS (excluding the five great lakes units), each is a downloadable zip file. Version 0.1 was computed at the National Center for Supercomputing Applications (NCSA) at the University of Illinois at Urban-Champaign in 2016 and is currently hosted at the Texas Advanced Computing Center (TACC). Please cite this data DOI and the following publications and data DOIs when you use HAND and the hydraulic property table: Liu, Yan Y., David R. Maidment, David G. Tarboton, Xing Zheng, and Shaowen Wang. 'A CyberGIS integration and computation framework for high-resolution continental-scale flood inundation mapping.' JAWRA Journal of the American Water Resources Association 54, no. 4 (2018): 770-784. DOI: 10.1111/1752-1688.12660 Zheng, Xing, David G. Tarboton, David R. Maidment, Yan Y. Liu, and Paola Passalacqua. 'River channel geometry and rating curve estimation using height above the nearest drainage.' JAWRA Journal of the American Water Resources Association 54, no. 4 (2018): 785-806. DOI: 10.1111/1752-1688.12661 Liu, Y. (2018). Height Above Nearest Drainage (HAND) for CONUS - v0.1, HydroShare, https://doi.org/10.4211/hs.69f7d237675c4c73938481904358c789 LICENSE FOR USE -- MAPS AND DATA DISCLAIMER This resource is shared under the Creative Commons Attribution CC BY, http://creativecommons.org/licenses/by/4.0/ MAPS AND DATA DISCLAIMER The Oak Ridge National Laboratory (ORNL) shall not be held liable for improper or incorrect use of the data described or information contained on this map or associated series of maps. The data and related map graphics are not legal, land survey or engineering documents and are not intended to be used as such. ORNL gives no warranty, express or implied, as to the accuracy, reliability, utility or completeness of this information. The user of these maps and data assumes all responsibility and risk for the use of the maps and data. ORNL disclaims all warranties, representations or endorsements either express or implied, with regard to the information contained in this map product, including, but not limited to, all implied warranties of merchantability, fitness for a particular purpose or non-infringement. This preliminary map product is for research and review purposes only. It is not intended to be used for emergency management operational or life safety decisions at the local or regional governmental level or by the general public. Users requiring information regarding hazardous conditions or meteorological conditions for specific geographic areas should consult directly with their city or county emergency management office.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Height Above Nearest Drainage (HAND) and Hydraulic Property Table for CONUS - Version 0.21 (20200601)

The continental flood inundation mapping (CFIM) framework is a high-performance computing (HPC)-based computational framework for the Height Above Nearest Drainage (HAND)-based inundation mapping methodology. Using the 10m Digital Elevation Model (DEM) data produced by U.S. Geological Survey (USGS) 3DEP (the 3-D Elevation Program) and the NHDPlus hydrography dataset produced by USGS and the U.S. Environmental Protection Agency (EPA), a hydrological terrain raster called HAND is computed for HUC6 units in the conterminous U.S. (CONUS). The value of each raster cell in HAND is an approximation of the relative elevation between the cell and its nearest water stream. Derived from HAND, a hydraulic property table is established to calculate river geometry properties for each of the 2.7 million river reaches covered by NHDPlus (5.5 million kilometers in total length). This table is a lookup table for water depth given an input stream flow value. Such lookup is available between water depth 0m and 25m at 1-foot interval. The flood inundation map is then computed by using HAND and this lookup table based on the near real-time water forecast from the National Water Model (NWM) at the National Oceanic and Atmospheric Administration (NOAA). HAND and the Hydraulic Property Table version v0.21 is an upgrade to v0.2. Version 0.21 applies NHD HR etching to DEMs for all HUC6 units before running the HAND workflow. This new version also fixed the issue of void areas along HUC boundary. This dataset comprises 331 HUC6 units for CONUS (excluding the five great lakes units), each is a downloadable zip file. Version 0.1 was computed at the National Center for Supercomputing Applications (NCSA) at the University of Illinois at Urban-Champaign in 2016 and is currently hosted at the Texas Advanced Computing Center (TACC). Please cite this data DOI and the following publications and data DOIs when you use HAND and the hydraulic property table: Liu, Yan Y. Height Above Nearest Drainage (HAND) and Hydraulic Property Table for CONUS - Version 0.2. (20200301). Oak Ridge Leadership Computing Facility. (2020). DOI: 10.13139/ORNLNCCS/1608331 Liu, Yan Y., David R. Maidment, David G. Tarboton, Xing Zheng, and Shaowen Wang. 'A CyberGIS integration and computation framework for high-resolution continental-scale flood inundation mapping.' JAWRA Journal of the American Water Resources Association 54, no. 4 (2018): 770-784. DOI: 10.1111/1752-1688.12660 Zheng, Xing, David G. Tarboton, David R. Maidment, Yan Y. Liu, and Paola Passalacqua. 'River channel geometry and rating curve estimation using height above the nearest drainage.' JAWRA Journal of the American Water Resources Association 54, no. 4 (2018): 785-806. DOI: 10.1111/1752-1688.12661 Liu, Y. (2018). Height Above Nearest Drainage (HAND) for CONUS, HydroShare, https://doi.org/10.4211/hs.69f7d237675c4c73938481904358c789 LICENSE FOR USE -- MAPS AND DATA DISCLAIMER This resource is shared under the Creative Commons Attribution CC BY, http://creativecommons.org/licenses/by/4.0/ MAPS AND DATA DISCLAIMER The Oak Ridge National Laboratory (ORNL) shall not be held liable for improper or incorrect use of the data described or information contained on this map or associated series of maps. The data and related map graphics are not legal, land survey or engineering documents and are not intended to be used as such. ORNL gives no warranty, express or implied, as to the accuracy, reliability, utility or completeness of this information. The user of these maps and data assumes all responsibility and risk for the use of the maps and data. ORNL disclaims all warranties, representations or endorsements either express or implied, with regard to the information contained in this map product, including, but not limited to, all implied warranties of merchantability, fitness for a particular purpose or non-infringement. This preliminary map product is for research and review purposes only. It is not intended to be used for emergency management operational or life safety decisions at the local or regional governmental level or by the general public. Users requiring information regarding hazardous conditions or meteorological conditions for specific geographic areas should consult directly with their city or county emergency management office.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

Landslide hazard mapping with selected dominant factors: A study case of Penang Island, Malaysia

Landslide is one of the destructive natural geohazards in Malaysia. In addition to rainfall as triggering factos for landslide in Malaysia, topographical and geological factors play important role in the landslide susceptibility analysis. Conventional topographic factors such as elevation, slope angle, slope aspect, plan curvature and profile curvature have been considered as landslide causative factors in many research works. However, other topographic factors such as diagonal length, surface area, surface roughness and rugosity have not been considered, especially for the research work in landslide hazard analysis in Malaysia. This paper presents landslide hazard mapping using Frequency Ratio (FR) and the study area is Penang Island of Malaysia. Frequency ratio approach is a variant of probabilistic method that is based on the observed relationships between the distribution of landslides and each landslide-causative factor. Landslide hazard map of Penang Island is produced by considering twenty-two (22) landslide causative factors. Among these twenty-two (22) factors, fourteen (14) factors are topographic factors. They are elevation, slope gradient, slope aspect, plan curvature, profile curvature, general curvature, tangential curvature, longitudinal curvature, cross section curvature, total curvature, diagonal length, surface area, surface roughness and rugosity. These topographic factors are extracted from the digital elevation model of Penang Island. The other eight (8) non-topographic factors considered are land cover, vegetation cover, distance from road, distance from stream, distance from fault line, geology, soil texture and rainfall precipitation. After considering all twenty-two factors for landslide hazard mapping, the analysis is repeated with fourteen dominant factors which are selected from the twenty-two factors. Landslide hazard map was segregated into four categories of risks, i.e. Highly hazardous area, Hazardous area, Moderately hazardous area and Not hazardous area. The maps was assessed using ROC (Rate of Curve) based on the area under the curve method (AUC). The result indicates an increase of accuracy from 77.76% (with all 22 factors) to 79.00% (with 14 dominant factors) in the prediction of landslide occurrence.

58 GEOSCIENCES↗

DOE OSTI2024

Geothermal Play Fairway Analysis of Low-Temperature Resources for Sedimentary Basin Geothermal Play Types: An Example in the Denver Basin

This project is part of a nationwide effort to highlight the advantages of incorporating low-temperature geothermal resource evaluation into the implementation of combined heat and power (CHP), and geothermal direct use (GDU) technologies (e.g., space heating and/or cooling). The initiative aims to hasten the nation's decarbonization process by exploring the potential for using low-temperature geothermal resources (< 150 Degrees Celsius) in selected sedimentary basins that have several population centers. The Play Fairway Analysis (PFA) techniques were modified from earlier studies of sedimentary basin geothermal play types (SBGPTs) that assessed the viability of low-temperature resources. The decision-making process for leveraging low-temperature geothermal resources for GDU and CHP applications is complex and considers a variety of factors, including geological, economic, and risk criteria. This study covers workflows, relevant datasets, python code, and both common and composite maps used to create low-temperature geothermal resource favorability maps for the Denver Basin, which extends across Colorado, Nebraska, and Wyoming. The replication of these methodologies in other SBGPTs can evaluate potential for low-temperature resources. The proposed geothermal PFA approach for low-temperature geothermal resources includes: (1) identifying available relevant data and grouping data sets into PFA criteria (e.g., geological, economic, and risk criteria); (2) analyzing data gaps enable future focalized exploration; (3) performing uncertainty quantification; (4) weighting relevant data; (5) developing favorability and common risk maps for low-temperature geothermal resources to identify potential locations for more focused data collection. This project will facilitate future deployment of CHP and GDU by providing data, tools, and a workflow applicable to low-temperature geothermal resources in sedimentary basins.

15 GEOTHERMAL ENERGY↗

DOE OSTI2016

Low Temperature Geothermal Play Fairway Analysis For The Appalachian Basin: Phase 1 Revised Report November 18, 2016

Geothermal energy is an attractive sustainable energy source. Yet project developers need confirmation of the resource base to warrant their time and financial resources. The Geothermal Play Fairway Analysis of the Appalachian Basin evaluated risk metrics that communicate the favorability of potential low-temperature geothermal energy resources in reservoirs more than 1000 m below the surface. This analysis is focused on the direct use of the heat, rather than on electricity production. Four risk factors of concern for direct-use geothermal plays in the Appalachian Basin portions of New York, Pennsylvania, and West Virginia are examined individually, and then in combination: 1) thermal resource quality, 2) natural reservoir quality, 3) induced seismicity, and 4) utilization opportunities. Uncertainty in the risk estimation is quantified. Based on these metrics, geothermal plays in the Appalachian Basin were identified as potentially viable for a variety of direct-use-heat applications. The methodologies developed in this project may be applied in other sedimentary basins as a foundation for low temperature (50-150 °C), direct use geothermal resource, risk, and uncertainty assessment. Three methods with which to combine the four risk factors were used. Among these, the averaging of the individual risk factors indicates the most favorable counties within the study area are the West Virginia counties of Monongalia, Harrison, Lewis (dubbed the Morgantown–Clarksburg play fairway), Putnam, and Kanawha (Charleston play fairway), the New York counties of Chemung and Steuben plus adjacent Bradford county in Pennsylvania (Corning–Ithaca play fairway), and the Pennsylvania counties of Mercer, Crawford, Erie, and Warren, and adjacent Chautauqua county in New York (together, the Meadville–Jamestown play fairway). These higher priority regions are surrounded by broader medium priority zones. Also worthy of additional exploration is a broad region near Pittsburgh Pennsylvania, for which the available geological data are insufficient to fully analyze the geological risks but yet the population is high. First, to assess the spatial variation in the depth to which one would need to drill to obtain geothermal temperatures that are useful to a future project, the project used bottom-hole temperature data from Appalachian Basin oil and gas exploration. These bottom hole temperature data are abundant but of low quality. Second, the project examined the potential for sufficient water flow rates through rocks to harvest heat from a geothermal well field, considering only natural reservoirs. This analysis provides a very incomplete picture of spatial variability of natural reservoirs because the oil and gas reservoir data lack key properties and are spatially biased toward those locations with profitable amounts of hydrocarbons in the rock pore spaces. Third, in light of the fact that earthquake activity has been induced in several states by subsurface work related to the oil and gas industry, this project examined the potential for similar activity in the Appalachian Basin. Acknowledging that data for such a task are insufficient, we utilized what was available: records of seismic activity, regional estimates of the orientations of stress in the rocks, and locations and orientations of zones of lateral change in rock properties at depths down to several kilometers below Earth’s surface. With these data, we created a first approximation of spatially variable risks for induced earthquakes. Because no data existed with which to test the reliability of these methods, the results have a high degree of uncertainty. Fourth, we examined the spatial variability of the above-the-ground factors that contribute to the economical viability of projects to tap low-temperature geothermal resources for direct-use. We worked principally with population density as a regionally known variable that would impact the cost of district heating. The resulting maps omit the costs of producing the hot water from the ground, because the below-ground costs are directly coupled to the thermal resource risk factor and natural reservoir risk factor – later analyses of those costs will be needed. The result of the district heating analysis is highly skewed: few census locations yielded a low estimated surface cost. The team also identified more than 165 prospects for high value direct-use geothermal energy opportunities such as industrial sites, university campuses, and federal facilities, among others. At the closure of this regional analysis, the most significant technical uncertainties are 1) reservoir distribution and capacities; 2) validity of thermal resource maps, and 3) the holistic estimation of Levelized Cost of Heat for favorable geological situations.

15 GEOTHERMAL ENERGY↗

DOE OSTI2023

A phased workflow to define permit‐ready locations for large volume CO 2 injection and storage

Abstract To‐date, only two UIC Class VI permits have been issued by the US Environmental Protection Agency. We illustrate a four‐phase workflow to first identify regional storage resources and then down‐select sites to yield permit‐ready locations that can accept and store large volumes of CO 2 . Specific permit requirements should guide objectives and define deliverables of respective workflow phases. In the first phase we used available regional data and screened structure and injection zones to locate resources that match CO 2 volumes planned to be captured. Available data were also used to assess presence and depth of usable groundwater, the key resource being protected via permitting. We then used advanced, closed‐form, analytical solutions (EASiTool) to estimate CO 2 injectivity into each hydrologically connected injection compartment. In the second phase we acquired and conditioned additional wireline logs and leased available seismic datasets. We interpreted the depositional systems from wireline well‐log character and mapped sandbody geometry to interpolate injection and confining‐zone distribution. Using available data, we mapped faults and locations of freshwater and overpressure (or other capacity‐limiting geologic parameters) in more detail. In the third phase, we used the augmented geologic data to develop a static model for the selected area, extracted the areas of highest interest, and generated and ran dynamic (flow) models. In a fourth phase, we reduced major uncertainties identified in earlier phases. Our case study indicates that to complete preparation of a permit application requires (1) improved lithologic characterization information (thicknesses and horizontal and vertical connectivity) and (2) better definition of poorly defined local faults. © 2023 The Authors. Greenhouse Gases: Science and Technology published by Society of Chemical Industry and John Wiley & Sons Ltd.

58 GEOSCIENCES↗

DOE OSTI2024

Characterization of Most Promising Sequestration Formations in the Rocky Mountain Region

The project Characterization of Most Promising Sequestration Formations in the Rocky Mountain Region is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The Rocky Mountain Carbon Capture and Storage (RMCCS) project investigated multiple geologic formations and characterized a local site on the Colorado Plateau for future CCS opportunities. The RMCCS project focused on the Cretaceous Dakota, Jurassic Entrada, and Pennsylvanian Weber Sandstones, the three largest regional formations. All formations in this project are potential CO2 storage resources for future power plants, natural gas processing plants, cement plants, and oil shale development projects. The area adjacent to Craig, Colorado, (Sand Wash Basin) was the area selected for detailed geologic characterization on the RMCCS project. The basin was selected in part because the geology can be extrapolated to other sites on the Colorado Plateau. Field mapping and seismic surveys were conducted to identify and evaluate the basin's structural configuration. A 9,745-foot deep characterization well was drilled to collect 131 feet of core and a suite of geophysical well log data. Petrophysical tests on samples of core were used to calibrate geophysical log data, which can be used to obtain storage resource estimates and evaluate associated uncertainty as well as simulate the hydrologic behavior of injected CO2. A detailed analysis of the primary formations (Dakota, Entrada and Weber sandstones) yielded a more accurate CO2 storage resource assessment for these formations within the Colorado Plateau; RMCCS estimates indicate a total CO2 storage resource of more than 38,000 million metric tons. The characterization of the Sand Wash Basin (2-D seismic surveys, multiple well logs and lithological, petrophysical and geochemical analyses) allowed for a detailed 3-D model to be constructed. The model served as the framework for analyses ranging from CO2 storage resource, injectivity, and subsurface flow to uncertainty estimates to evaluation of risk.

2-D seismic↗

DOE OSTI2014

Modeling CO 2 Sequestration in Saline Aquifer and Depleted Oil Reservoir To Evaluate Regional CO 2 Sequestration Potential of Ozark Plateau Aquifer System, South-Central Kansas

1. Drilled, cored, and logged three wells to the basement and collecting more than 2,700 ft of conventional core; obtained 20 m i2 of multicomponent 3D seismic imaging and merged and reprocessed more than 125 mi 2 of existing 3D seismic data for use in modeling CO 2 - EOR oil recovery and CO 2 storage in five oil fields in southern Kansas. 2. Determined the technical feasibility of injecting and sequestering CO 2 in a set of four depleted oil reservoirs in the Cutter, Pleasant Prairie South, Eubank, and Shuck fields in southwest Kansas; of concurrently recovering oil from those fields; and of quantifying the volumes of CO2 sequestered and oil recovered during the process. 3. Formed a consortium of six oil operating companies, five of which own and operate the four fields. The consortium became part of the Southwest Kansas CO 2 -EOR Initiative for the purpose of sharing data, knowledge, and interest in understanding the potential for CO 2 -EOR in Kansas. 4. Built a regional well database covering 30,000 mi 2 and containing stratigraphic tops from ~90,000 wells; correlated 30 major stratigraphic horizons; digitized key wells, including wireline logs and sample logs; and analyzed more than 3,000 drill stem tests to establish that fluid levels in deep aquifers below the Permian evaporites are not connected to the surface and therefore pressures are not hydrostatic. Connectivity with the surface aquifers is lacking because shale aquitards and impermeable evaporite layers consist of both halite and anhydrite. 5. Developed extensive web applications and an interactive mapping system that do the following: a. Facilitate access to a wide array of data obtained in the study, including core descriptions and analyses, sample logs, digital (LAS) well logs, seismic data, gravity and magnetics maps, structural and stratigraphic maps, inferred fault traces, earthquakes, Class I and II disposal wells, and surface lineaments. b. Provide real-time analysis of the project dataset, including automated integration and viewing of well logs, core, core analyses, brine chemistry, and stratigraphy using the Java Profile app. A cross-section app allows for the display of log data for up to four wells at a time. 6. Integrated interpretations from the project’s interactive web-based mapping system to gain insights to aid in assessing the efficacy of geologic CO 2 storage in Kansas and insights toward understanding recent seismicity to aid in evaluating induced vs. naturally occurring earthquakes. 7. Developed a digital type-log system, including web-based software to modify and refine stratigraphic nomenclature to provide stakeholders a common means for communication about the subsurface. 8. Contracted use of a nuclear magnetic resonance (NMR) log and ran it slowly to capture response and characterize larger pores common for carbonate reservoirs. Used NMR to extend core analyses to apply permeability, relative permeability to CO 2 , and capillary pressure to the major rock types, each uniquely expressed as a reservoir quality index (RQI), present in the Mississippian and Arbuckle rocks. 9. Characterized and evaluated the possible role of microbes in dense brines. Used microbes to compliment H/O stable isotopes to fingerprint brine systems. Used perforation/swabbing to obtain samples from multiple hydrostratigraphic units and confirmed equivalent results using less expensive drill stem tests (DST). 10. Used an integrated approach from whole core, logs, tests, and seismic to verify and quantify properties of vuggy, brecciated, and fractured carbonate intervals. 11. Used complex geocellular static and dynamic models to evaluate regional storage capacity using large parallel processing. 12. Carbonates are complex reservoirs and CO 2 -EOR needs to move to the next generation to increase effectiveness of CO 2 and efficiency and safety of the injection.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Modeling CO 2 Sequestration in Saline Aquifer and Depleted Oil Reservoir To Evaluate Regional CO 2 Sequestration Potential of Ozark Plateau Aquifer System, South-Central Kansas

1. Drilled, cored, and logged three wells to the basement and collecting more than 2,700 ft of conventional core; obtained 20 m i2 of multicomponent 3D seismic imaging and merged and reprocessed more than 125 mi 2 of existing 3D seismic data for use in modeling CO 2 - EOR oil recovery and CO 2 storage in five oil fields in southern Kansas. 2. Determined the technical feasibility of injecting and sequestering CO 2 in a set of four depleted oil reservoirs in the Cutter, Pleasant Prairie South, Eubank, and Shuck fields in southwest Kansas; of concurrently recovering oil from those fields; and of quantifying the volumes of CO2 sequestered and oil recovered during the process. 3. Formed a consortium of six oil operating companies, five of which own and operate the four fields. The consortium became part of the Southwest Kansas CO 2 -EOR Initiative for the purpose of sharing data, knowledge, and interest in understanding the potential for CO 2 -EOR in Kansas. 4. Built a regional well database covering 30,000 mi 2 and containing stratigraphic tops from ~90,000 wells; correlated 30 major stratigraphic horizons; digitized key wells, including wireline logs and sample logs; and analyzed more than 3,000 drill stem tests to establish that fluid levels in deep aquifers below the Permian evaporites are not connected to the surface and therefore pressures are not hydrostatic. Connectivity with the surface aquifers is lacking because shale aquitards and impermeable evaporite layers consist of both halite and anhydrite. 5. Developed extensive web applications and an interactive mapping system that do the following: a. Facilitate access to a wide array of data obtained in the study, including core descriptions and analyses, sample logs, digital (LAS) well logs, seismic data, gravity and magnetics maps, structural and stratigraphic maps, inferred fault traces, earthquakes, Class I and II disposal wells, and surface lineaments. b. Provide real-time analysis of the project dataset, including automated integration and viewing of well logs, core, core analyses, brine chemistry, and stratigraphy using the Java Profile app. A cross-section app allows for the display of log data for up to four wells at a time. 6. Integrated interpretations from the project’s interactive web-based mapping system to gain insights to aid in assessing the efficacy of geologic CO 2 storage in Kansas and insights toward understanding recent seismicity to aid in evaluating induced vs. naturally occurring earthquakes. 7. Developed a digital type-log system, including web-based software to modify and refine stratigraphic nomenclature to provide stakeholders a common means for communication about the subsurface. 8. Contracted use of a nuclear magnetic resonance (NMR) log and ran it slowly to capture response and characterize larger pores common for carbonate reservoirs. Used NMR to extend core analyses to apply permeability, relative permeability to CO 2 , and capillary pressure to the major rock types, each uniquely expressed as a reservoir quality index (RQI), present in the Mississippian and Arbuckle rocks. 9. Characterized and evaluated the possible role of microbes in dense brines. Used microbes to compliment H/O stable isotopes to fingerprint brine systems. Used perforation/swabbing to obtain samples from multiple hydrostratigraphic units and confirmed equivalent results using less expensive drill stem tests (DST). 10. Used an integrated approach from whole core, logs, tests, and seismic to verify and quantify properties of vuggy, brecciated, and fractured carbonate intervals. 11. Used complex geocellular static and dynamic models to evaluate regional storage capacity using large parallel processing. 12. Carbonates are complex reservoirs and CO 2 -EOR needs to move to the next generation to increase effectiveness of CO 2 and efficiency and safety of the injection.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Modeling CO 2 Sequestration in Saline Aquifer and Depleted Oil Reservoir To Evaluate Regional CO 2 Sequestration Potential of Ozark Plateau Aquifer System, South-Central Kansas

1. Drilled, cored, and logged three wells to the basement and collecting more than 2,700 ft of conventional core; obtained 20 m i2 of multicomponent 3D seismic imaging and merged and reprocessed more than 125 mi 2 of existing 3D seismic data for use in modeling CO 2 - EOR oil recovery and CO 2 storage in five oil fields in southern Kansas. 2. Determined the technical feasibility of injecting and sequestering CO 2 in a set of four depleted oil reservoirs in the Cutter, Pleasant Prairie South, Eubank, and Shuck fields in southwest Kansas; of concurrently recovering oil from those fields; and of quantifying the volumes of CO2 sequestered and oil recovered during the process. 3. Formed a consortium of six oil operating companies, five of which own and operate the four fields. The consortium became part of the Southwest Kansas CO 2 -EOR Initiative for the purpose of sharing data, knowledge, and interest in understanding the potential for CO 2 -EOR in Kansas. 4. Built a regional well database covering 30,000 mi 2 and containing stratigraphic tops from ~90,000 wells; correlated 30 major stratigraphic horizons; digitized key wells, including wireline logs and sample logs; and analyzed more than 3,000 drill stem tests to establish that fluid levels in deep aquifers below the Permian evaporites are not connected to the surface and therefore pressures are not hydrostatic. Connectivity with the surface aquifers is lacking because shale aquitards and impermeable evaporite layers consist of both halite and anhydrite. 5. Developed extensive web applications and an interactive mapping system that do the following: a. Facilitate access to a wide array of data obtained in the study, including core descriptions and analyses, sample logs, digital (LAS) well logs, seismic data, gravity and magnetics maps, structural and stratigraphic maps, inferred fault traces, earthquakes, Class I and II disposal wells, and surface lineaments. b. Provide real-time analysis of the project dataset, including automated integration and viewing of well logs, core, core analyses, brine chemistry, and stratigraphy using the Java Profile app. A cross-section app allows for the display of log data for up to four wells at a time. 6. Integrated interpretations from the project’s interactive web-based mapping system to gain insights to aid in assessing the efficacy of geologic CO 2 storage in Kansas and insights toward understanding recent seismicity to aid in evaluating induced vs. naturally occurring earthquakes. 7. Developed a digital type-log system, including web-based software to modify and refine stratigraphic nomenclature to provide stakeholders a common means for communication about the subsurface. 8. Contracted use of a nuclear magnetic resonance (NMR) log and ran it slowly to capture response and characterize larger pores common for carbonate reservoirs. Used NMR to extend core analyses to apply permeability, relative permeability to CO 2 , and capillary pressure to the major rock types, each uniquely expressed as a reservoir quality index (RQI), present in the Mississippian and Arbuckle rocks. 9. Characterized and evaluated the possible role of microbes in dense brines. Used microbes to compliment H/O stable isotopes to fingerprint brine systems. Used perforation/swabbing to obtain samples from multiple hydrostratigraphic units and confirmed equivalent results using less expensive drill stem tests (DST). 10. Used an integrated approach from whole core, logs, tests, and seismic to verify and quantify properties of vuggy, brecciated, and fractured carbonate intervals. 11. Used complex geocellular static and dynamic models to evaluate regional storage capacity using large parallel processing. 12. Carbonates are complex reservoirs and CO 2 -EOR needs to move to the next generation to increase effectiveness of CO 2 and efficiency and safety of the injection.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Modeling CO 2 Sequestration in Saline Aquifer and Depleted Oil Reservoir To Evaluate Regional CO 2 Sequestration Potential of Ozark Plateau Aquifer System, South-Central Kansas

1. Drilled, cored, and logged three wells to the basement and collecting more than 2,700 ft of conventional core; obtained 20 m i2 of multicomponent 3D seismic imaging and merged and reprocessed more than 125 mi 2 of existing 3D seismic data for use in modeling CO 2 - EOR oil recovery and CO 2 storage in five oil fields in southern Kansas. 2. Determined the technical feasibility of injecting and sequestering CO 2 in a set of four depleted oil reservoirs in the Cutter, Pleasant Prairie South, Eubank, and Shuck fields in southwest Kansas; of concurrently recovering oil from those fields; and of quantifying the volumes of CO2 sequestered and oil recovered during the process. 3. Formed a consortium of six oil operating companies, five of which own and operate the four fields. The consortium became part of the Southwest Kansas CO 2 -EOR Initiative for the purpose of sharing data, knowledge, and interest in understanding the potential for CO 2 -EOR in Kansas. 4. Built a regional well database covering 30,000 mi 2 and containing stratigraphic tops from ~90,000 wells; correlated 30 major stratigraphic horizons; digitized key wells, including wireline logs and sample logs; and analyzed more than 3,000 drill stem tests to establish that fluid levels in deep aquifers below the Permian evaporites are not connected to the surface and therefore pressures are not hydrostatic. Connectivity with the surface aquifers is lacking because shale aquitards and impermeable evaporite layers consist of both halite and anhydrite. 5. Developed extensive web applications and an interactive mapping system that do the following: a. Facilitate access to a wide array of data obtained in the study, including core descriptions and analyses, sample logs, digital (LAS) well logs, seismic data, gravity and magnetics maps, structural and stratigraphic maps, inferred fault traces, earthquakes, Class I and II disposal wells, and surface lineaments. b. Provide real-time analysis of the project dataset, including automated integration and viewing of well logs, core, core analyses, brine chemistry, and stratigraphy using the Java Profile app. A cross-section app allows for the display of log data for up to four wells at a time. 6. Integrated interpretations from the project’s interactive web-based mapping system to gain insights to aid in assessing the efficacy of geologic CO 2 storage in Kansas and insights toward understanding recent seismicity to aid in evaluating induced vs. naturally occurring earthquakes. 7. Developed a digital type-log system, including web-based software to modify and refine stratigraphic nomenclature to provide stakeholders a common means for communication about the subsurface. 8. Contracted use of a nuclear magnetic resonance (NMR) log and ran it slowly to capture response and characterize larger pores common for carbonate reservoirs. Used NMR to extend core analyses to apply permeability, relative permeability to CO 2 , and capillary pressure to the major rock types, each uniquely expressed as a reservoir quality index (RQI), present in the Mississippian and Arbuckle rocks. 9. Characterized and evaluated the possible role of microbes in dense brines. Used microbes to compliment H/O stable isotopes to fingerprint brine systems. Used perforation/swabbing to obtain samples from multiple hydrostratigraphic units and confirmed equivalent results using less expensive drill stem tests (DST). 10. Used an integrated approach from whole core, logs, tests, and seismic to verify and quantify properties of vuggy, brecciated, and fractured carbonate intervals. 11. Used complex geocellular static and dynamic models to evaluate regional storage capacity using large parallel processing. 12. Carbonates are complex reservoirs and CO 2 -EOR needs to move to the next generation to increase effectiveness of CO 2 and efficiency and safety of the injection.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2017

Characterization of the Triassic Newark Basin of New York and New Jersey for geologic storage of carbon dioxide

The Newark Basin is a Triassic-aged rift basin underlying densely populated, industrialized sections of New York, New Jersey and Pennsylvania. The Basin is an elongate half-graben encompassing an area of more than 7,510 square-kilometers (2,900 square-miles), and could represent a key storage component for commercial scale management of carbon dioxide emissions via geologic sequestration. The project team first acquired published reports, surface and subsurface maps, and seismic data, which formed the basis for a three-dimensional model framework for the northern end of the Basin incorporating stratigraphic, hydrologic, and water quality data. Field investigations included drilling, coring, and logging of two stratigraphic test borings in Clarkstown, NY (Exit 14 Tandem Lot Well No. 1), drilled to a depth of 2,099 meters (6,885 feet); and Palisades, NY (Lamont-Doherty Earth Observatory Test Well No. 4) drilled to a depth of 549 meters (1,802 feet). Two two-dimensional seismic reflection data lines arrayed perpendicularly were acquired by Schlumberger/WesternGeco to help characterize the structure and stratigraphy and as part of pre-drilling field screening activities for the deep stratigraphic borehole. A total of 47 meters (155 feet) of continuous whole core was recovered from the Tandem Lot boring from depths of 1,393 meters (4,570 feet) to 1,486 meters (4,877 feet). Twenty-five horizontal rotary cores were collected in mudstones and sandstones in the surface casing hole and fifty-two cores were taken in various lithologies in the deep borehole. Rotary core plugs were analyzed by Weatherford Laboratories for routine and advanced testing. Rotary core plug trim end thin sections were evaluated by the New York State Museum for mineralogical analysis and porosity estimation. Using core samples, Lawrence Berkley National Laboratory designed and completed laboratory experiments and numerical modeling analyses to characterize the dissolution and reaction of carbon dioxide with formation brine and minerals, and resulting effects on injection rate, pressure, effective storage volume, and carbon dioxide migration within a prospective sandstone reservoir. $$Three potential porous and permeable sandstone units were identified in the Passaic Formation at the New York State Thruway Exit 14 location. Potential Flow Unit 1, at a depth of 643 meters (2,110 feet) to 751 meters (2,465 feet); Potential Flow Unit 2 at a depth of 853 meters (2,798 feet) to 1,000 meters (3,280 feet); and Potential Flow Unit 3, at a depth of 1,114 meters (3,655 feet) to 1,294 meters (4,250 feet). Reactive transport simulations of interactions between carbon dioxide, brine and formation minerals were carried out to evaluate changes in formation water chemistry, mineral precipitation and dissolution reactions, and any potential resulting effects on formation permeability. The experimental and modeling analyses suggest that mineral precipitation and dissolution reactions (within the target formation) are not expected to lead to significant changes to the underground hydrologic system over time frames (~30 years) typically relevant for carbon dioxide injection operations. Key findings of this basin characterization study include an estimate of carbon dioxide storage capacity in the Newark Basin. Assuming an average porosity of twelve percent and an aquifer volume of 6.1E+12 meters3, calculated ranges of likely storage capacity range from 1.9 – 20.2 gigatonnes under high temperature (low carbon dioxide density) conditions; and 2.9 – 30.2 gigatonnes under low temperature (low carbon dioxide density) conditions. Intra-basin faulting, geometry of the Palisades Sill, and the presence of altered meta-sediments above and below the Sill, increase potential compartmentalization within the basin. A structural/stratigraphic trap type may occur where porous/permeable sediments are cross-cut by the Palisades Sill. Potential injection intervals are present within the Stockton Formation of the Newark Basin. Additional porous/permeable intervals may be present within sandstones of the Passaic Formation, increasing projected storage capacity. Deeper wedges of strata are likely present in the deeper portions of the basin in southern New York and into northern New Jersey. Abundant mudstones are present in the Passaic, Lockatong, and Stockton Formations. These intervals have the requisite petrophysical properties to form effective primary and secondary containment intervals to industrial-scale sequestration of carbon dioxide in the Newark Basin. Hydro-thermally altered meta-sediments in the region immediately surrounding the top and base of the Palisades Sill is devoid of porosity/permeability and forms an additional effective lateral/vertical sealing cap rock.

58 GEOSCIENCES↗

DOE OSTI2022

Predicting Rare Earth Element Potential in Produced and Geothermal Waters of the United States via Emergent Self-Organizing Maps

This work applies emergent self-organizing map (ESOM) techniques, a form of machine learning, in the multidimensional interpretation and prediction of rare earth element (REE) abundance in produced and geothermal waters in the United States. Visualization of the variables in the ESOM trained using the input data shows that each REE, with the exception of Eu, follows the same distribution patterns and that no single parameter appears to control their distribution. Cross-validation, using a random subsample of the starting data and only using major ions, shows that predictions are generally accurate to within an order of magnitude. Using the same approach, an abridged version of the U.S. Geological Survey Produced Waters Database, Version 2.3 (which includes both data from produced and geothermal waters) was mapped to the ESOM and predicted values were generated for samples that contained enough variables to be effectively mapped. Results show that in general, produced and geothermal waters are predicted to be enriched in REEs by an order of magnitude or more relative to seawater, with maximum predicted enrichments in excess of 1000-fold. Cartographic mapping of the resulting predictions indicates that maximum REE concentrations exceed values in seawater across the majority of geologic basins investigated and that REEs are typically spatially co-associated. The factors causing this co-association were not determined from ESOM analysis, but based on the information currently available, REE content in produced and geothermal waters is not directly controlled by lithology, reservoir temperature, or salinity.

Engle, Mark A. (ORCID:0000000152587374)↗

DOE OSTI2021

Femtosecond Laser Desorption Postionization MS vs ToF-SIMS Imaging for Uncovering Biomarkers Buried in Geological Samples

The study of lipid molecular fossils by traditional biomarker analysis requires bulk sample crushing, followed by solvent extraction, and then the analysis of the extract by gas chromatography-mass spectrometry (GC–MS). This traditional analysis mixes all organic compounds in the sample regardless of their origins, with a loss of information on the spatial distribution of organic molecules within the sample. These shortcomings can be overcome using the chemical mapping of intact samples. Spectroscopic techniques such as UV fluorescence or Raman spectroscopy, laser ablation inductively coupled plasma mass spectrometry, and time-of-flight secondary ion mass spectrometry (ToF-SIMS) are among those elemental and molecular mapping techniques. This study employed femtosecond (fs) laser ablation combined with single-photon ionization, a method called fs-laser desorption postionization mass spectrometry (fs-LDPI-MS). In this work, a pulsed ~75 fs, 800 nm laser was used to ablate the geological sample, which was then photoionized after a few microseconds by a pulsed 7.9 eV vacuum ultraviolet laser. An organic carbon-rich geological sample was used for this study to map hydrocarbon biomarkers in sediments that were previously studied by GC–MS. The petrography of this sample was examined by optical and fluorescence microscopy. It is demonstrated here that fs-LDPI-MS combined with petrography for multimodal imaging can expose buried compounds within the sample via in situ layer removal. When used in conjunction with traditional organic geochemical analysis, this method has the potential to determine the spatial distribution of organic biomarkers in geological material. Finally, fs-LDPI-MS imaging data are compared with ToF-SIMS imaging that is commonly used for such studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Hydrologic Impacts of a Strike-Slip Fault Zone: Insights from Joint 3D Body-Wave Tomography of Rock Valley

The Rock Valley fault zone (RVFZ), an intraplate strike-slip fault zone in the southern Nevada National Security Site (NNSS), hosted a series of very shallow (<3 km) earthquakes in 1993. The RVFZ may also have hydrological significance within the NNSS, potentially playing a role in regional groundwater flow, but there is a lack of local hydrological data. In the Spring of 2021, we collected active-source accelerated weight drop seismic data over part of the RVFZ to better characterize the shallow subsurface. Here, we manually picked ~17,000 P-wave travel times and over 14,000 S-wave travel times, which were inverted for P-wave velocity (V P ), S-wave velocity (V S ), and V P /V S ratio in a 3D joint tomographic inversion scheme. Seismic velocities are imaged as deep as ~700 m in areas and generally align with geologic and structural expectations. V P and V S are relatively reduced near mapped and inferred faults, with the most prominent lower V P and V S zone around the densest collection of faults. We image V P /V S ratios ranging from ~1.5 to ~2.4, the extremes of which occur at a depth of ~100 m and are juxtaposed across a fault. One possible interpretation of the imaged seismic velocities is enhanced fault damage near the densest collection of faults with relatively higher porosity and/or crack density at ~100 m depth, with patches of semiperched groundwater present in the sedimentary rock in higher V P /V S areas and drier rock in lower V P /V S areas. A relatively higher V P /V S area beneath the densest faults persists at depth, which suggests percolation of groundwater via the fault damage zone to the regionally connected lower carbonate aquifer. Potentially, the presence and movement of groundwater may have played a role in the 1993 earthquake aftershocks.

58 GEOSCIENCES↗

DOE OSTI2015

Discovering Geothermal Systems in the Great Basin Region: An Integrated Geologic, Geochemical, and Geophysical Approach for Establishing Geothermal Play Fairways

The Great Basin currently hosts ~24 geothermal power plants with more than 600 MW of capacity, but estimates suggest that the region is capable of producing much greater amounts of geothermal energy. Recent studies indicate that a combination of geological, geochemical, and geophysical features characterize the most robust geothermal systems. Considering that most of the geothermal resources in this region are blind, it is imperative that the favorable characteristics for geothermal activity be synthesized and techniques perfected for the discovery of new viable systems. We therefore utilized a multi-disciplinary approach and synthesized 9 geologic, geochemical, and geophysical parameters to produce a new detailed geothermal play fairway map of a large transect across the Great Basin region (96,000 km2) from west-central to eastern Nevada, with the primary objective of facilitating discovery of commercial-grade, blind geothermal fields (i.e. systems with no surface hot springs or fumaroles). Our analysis included: 1) structural settings (i.e., patterns of faulting), 2) age of recent faulting, 3) slip rates on recent faults, 4) regional-scale strain rates, 5) the tendency of faults to slip or dilate based on their orientation in the regional stress field, 6) earthquake density, 7) gravity data, 8) temperature at 3 km depth, and 9) geochemistry from springs and wells. These parameters were grouped into subsets to delineate rankings for local permeability, regional permeability, and heat, which collectively defined the geothermal play fairways (i.e., most likely locations for significant geothermal fluid flow). A major challenge was developing appropriate weightings of individual data types to best predict permeability and overall geothermal potential. Rigorous statistical methods, utilizing 34 benchmarks of known relatively high-temperature (>130°C) geothermal systems in the region, were employed to determine the hierarchal weights of each parameter. The new geothermal potential maps will help to reduce the risks in prospecting for new geothermal systems throughout the region.

15 GEOTHERMAL ENERGY↗

DOE OSTI2020

Regional Assessment for the CO 2 Storage Potential in the Northern Niagaran Pinnacle Reef Trend

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. The report describes the regional analysis of the Northern Niagaran Pinnacle Reef Trend (NNPRT) in Northern Michigan as part of the MRCSP programs Regional Assessment goals. The goal of this task was to perform an initial assessment of the geologic storage capacity and injectivity of the Northern Fairway of Michigan’s Niagaran Reef Trend. This was achieved through three main objectives: 1) development of a reef atlas, 2) understanding the regional trends in geology, and 3) estimating CO 2 and CO 2 -EOR resources. Data collected across multiple subtasks were integrated into a mappable database that could be accessed from many of the leading software. The reef atlas was used to map the reefs and associated characteristics across the entire NNPRT. Several methodologies were explored to estimate the CO 2 and CO 2 -EOR resources across the trend which included three scenarios: 1) storage only using fluid substitution and volumetric estimations, 2) CO 2 -EOR which applied proximity analysis concepts and measured performance metrics from the Core Energy reefs to predict performance at all oil reefs, and 3) enhanced storage scenario which combined CO 2 -EOR with fluid substitution to represent maximized storage after completion of CO 2 -EOR.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2015

Estimating and mapping ecological processes influencing microbial community assembly

Ecological community assembly is governed by a combination of (i) selection resulting from among-taxa differences in performance; (ii) dispersal resulting from organismal movement; and (iii) ecological drift resulting from stochastic changes in population sizes. The relative importance and nature of these processes can vary across environments. Selection can be homogeneous or variable, and while dispersal is a rate, we conceptualize extreme dispersal rates as two categories; dispersal limitation results from limited exchange of organisms among communities, and homogenizing dispersal results from high levels of organism exchange. To estimate the influence and spatial variation of each process we extend a recently developed statistical framework, use a simulation model to evaluate the accuracy of the extended framework, and use the framework to examine subsurface microbial communities over two geologic formations. For each subsurface community we estimate the degree to which it is influenced by homogeneous selection, variable selection, dispersal limitation, and homogenizing dispersal. Our analyses revealed that the relative influences of these ecological processes vary substantially across communities even within a geologic formation. We further identify environmental and spatial features associated with each ecological process, which allowed mapping of spatial variation in ecological-process-influences. The resulting maps provide a new lens through which ecological systems can be understood; in the subsurface system investigated here they revealed that the influence of variable selection was associated with the rate at which redox conditions change with subsurface depth.

54 ENVIRONMENTAL SCIENCES↗