Simulated CO2 storage efficiency in sandstone and carbonate reservoirs: CO2-SCREEN tool upgrade
For presentation at the 2023 FECM/NETL Carbon Management Research Project Review Meeting, Pittsburgh, PA, August 28-September 1, 2023.
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.
For presentation at the 2023 FECM/NETL Carbon Management Research Project Review Meeting, Pittsburgh, PA, August 28-September 1, 2023.
Carbonated brine storage is a low-risk geologic CO2 sequestration strategy compared to traditional supercritical CO2 injection and storage. Our study provides an evaluation methodology for CBI that can be used for other potential pilot injections in suitable fields and opens the opportunity to work with industry for CBI using existing Class II wells to gain CO2 tax credits.
During CO 2 injection into geological formations, petrophysical and geomechanical properties of host formations can be altered due to mineral dissolution and precipitation. Field and laboratory results have shown that sandstone and siltstone can be altered by CO 2 -water mixtures, but few quantitative studies have been performed to fully investigate underlying mechanisms. Based on the hypothesis that CO 2 -water mixtures alter the integrity of rock structure by attacking cements rather than grains, we attempt to explain the degradation of cementation due to long-term contact with CO 2 and water and mechanisms for changes in rock mechanical properties. Many sandstones, including calcite-cemented quartzitic sandstone, chlorite-cemented quartzitic sandstone, and hematite-cemented quartzitic sandstone, contain interparticle cements that are more readily affected by CO 2 -water mixtures than grains. A model that couples the discrete element method and the bonded-particle model is used to perform simulations of indentation tests on synthetic rocks with crystal and random packings. The model is verified against the analytical cavity expansion model and validated against laboratory indentation tests on Entrada sandstone with and without CO 2 alteration. Sensitivity analysis is performed for cementation microscopic parameters including stiffness, size, axial, and shear strength. The simulation results indicate that the CO 2 -related degradation of mechanical properties in bleached Entrada sandstone can be attributed to the reduction of cement size rather than cement strength. Finally, our study indicates that it is possible to describe the CO 2 -related rock alteration through particle-scale mechanisms.
The reactivity of sandstones was studied under conditions relevant to the injection of supercritical carbon dioxide in the context of carbon geosequestration. The emphasis of the study was on the reactivity of iron-bearing minerals when exposed to supercritical CO 2 (scCO 2 ) and scCO 2 with commingled aqueous solutions containing H 2 S and/or SO 2 . Flow through and batch experiments were conducted. Results indicate that sandstones, irrespective of their mineralogy, are not reactive when exposed to pure scCO2 or scCO 2 with commingled aqueous solutions containing H 2 S and/or SO 2 under conditions simulating the environment near the injection point (flow through experiments). However, sandstones are reactive under conditions simulating the edge of the injected CO 2 plume or ahead of the plume (batch experiments). Sandstones containing hematite (red sandstone) are particularly reactive. The composition of the reaction products is strongly dependent on the composition of the aqueous phase. The presence of dissolved sulfide leads to the conversion of hematite into pyrite and siderite. The relative amount of the pyrite and siderite is influenced by the ionic strength of the solution. Little reactivity is observed when sulfite is present in the aqueous phase. Sandstones without hematite (grey sandstones) show little reactivity regardless of the solution composition.
Utah has a wealth of potential geological reservoirs for carbon dioxide storage (CS) and a long history of geologic research resulting in an abundance of available subsurface data to evaluate CS potential. Reservoirs may include sandstone, carbonate, and basalt; these rock types are plentiful in Utah’s subsurface and the complex Phanerozoic history throughout the state requires evaluating each geologic region individually for promising reservoir-seal pairs for CO2 storage. Classifying Utah by geologic provinces (or “geo-regions”) allows for customized thinking about suitable CS reservoir and seal distribution, CO2 point sources, land use, and existing infrastructure. Preliminary results from this study highlight the geologic CS potential across 15 geo-regions. Four regions stand out as having high CS potential: the Uinta Basin, San Rafael Swell, Paradox Basin, and the southern Basin and Range Province. The Uinta Basin and San Rafael Swell geo-regions are well suited for CS and have several projects ongoing to evaluate Cretaceous Frontier and Naturita Formations, Jurassic Navajo Sandstone and Entrada Sandstone, and Permian Weber Sandstone reservoir units that lie beneath robust sealing units like the ~5000-ft-thick Mancos Shale and Carmel Formation. Reservoirs such as the Navajo and Weber Sandstones have been demonstrated to be suitable reservoirs through a long history of oil and gas exploration in Utah. New areas of interest include the southern Basin and Range in southwest Utah, where the Jurassic Navajo Sandstone is overlain by the sealing Carmel Formation at suitable depths (>3000 ft), and has good porosities based on outcrop analogue data. Just to the north (e.g., central Basin and Range), legacy wells and 2D seismic data show possible salt and subsurface basalt flows that may provide additional possible CS reservoirs and seals. This geo-region also has the advantage of being coupled with geothermal energy resources that may be used to power burgeoning direct air capture technologies. In the Paradox Basin of southeastern Utah, the Leadville Limestone is a potential storage reservoir beneath the thick (4000–5000 ft), salt-bearing Pennsylvanian Paradox Formation. Although the northern and western parts of Utah offer CS potential, these areas typically contain less infrastructure and subsurface penetrations, creating geologic uncertainty associated with subsurface seals and reservoirs due to a lack of data. Overall, this statewide assessment and ranking is the first step to aid in evaluating CS potential across Utah and provides a foundation for future research in the most favorable locations.
A microindentation scanner is constructed that measures the spatial fluctuation in the elastic properties of natural rocks. This novel instrument performs automated indentation tests on the surface of a rock slab and outputs 2-D maps of the indentation modulus at submillimeter resolution. Maps obtained for clean, well-consolidated, sandstone are presented and demonstrate the capabilities of the instrument. We postulate that the elastic structure of sandstones correlates well with their visual appearance. Moreover, we observe that the probability distribution of the indentation modulus fluctuations across the slab surfaces can be modeled using a lognormal probability density function. To illustrate possible use of the data obtained with the microindentation scanner, we use roughly 10 cm × 10 cm scans with millimeter resolution over four sandstone planar slabs to numerically compute the overall drained elastic moduli for each sandstone sample. We demonstrate that such numerically computed moduli are well modeled using the multicomponent form of the Hashin-Shtrikman lower bound that employs the observed lognormal probability distribution for the mesoscopic-scale moduli (the geometric mean works almost the same). We also compute the seismic attenuation versus frequency associated with wave-induced fluid flow between the heterogeneities in the scanned sandstones and report relatively small values for the inverse quality factor (Q -1 <10 -2 ) in the seismic frequency band 10 2 Hz 4 Hz. The numerically computed frequency dependence in the attenuation varies from one type of sandstone to another, and we determine vital anisotropy in the attenuation associated with waves propagating in different directions.
Well logs and 3D seismic reflection data are integrated to image Lower Miocene depocenters of the upper Texas and westernmost Louisiana coastal and offshore areas. Although previously interpreted at a large scale (3rd order cycles), the detailed stratigraphy and depositional history of the early Miocene succession has not been fully developed. The upper lower Miocene interval from Robulus L (MFS10) to Amphistegina B (MFS9) has been further subdivided into five 4th order cycles (each 200 to 300 m thick) to provide finer scale stratigraphic detail and interpretation of depositional environments. Results of the finer scale correlations are presented in a series of sandstone percent maps. Overall, the maps display a strike-elongate (sub-parallel to the present day coastline) orientation for the sandstone bodies which thin to the southeast (basinward). Net sandstone thickness is relatively low in older cycles and increases in younger cycles suggesting an overall progradation of the deltaic system into the mapped area through time. Continued marine regression and southward movement of the shoreline are also confirmed by seismic amplitude 2 maps. After a major transgression associated with Robulus L deltaic progradation occurred under rising sea level conditions. Deltaic sediments prograded southward and sandstone brought to the upper Texas coast was distributed laterally by longshore currents to form strandplain coast lines. Maximum regression occurred during the interval MFS9_2 to MFS9_3 when a deltaic depocenter formed offshore Texas in east High Island area. However, sandstone thickness progressively decreased laterally to the east (offshore Louisiana) where more shale was present and marine processes reworked deltaic derived sandstones into shore parallel bars. There was a clear retreat of the shoreline during the youngest cycle marking the beginning of transgression associated with Amphistegina B . Detailed stratigraphic interpretation (at a 4th order scale) shows more variability in the dominant shoreline processes during early Miocene than previously thought. Recognition of sedimentary systems variability at a 4th order scale is critical for improved hydrocarbon exploration and in understanding potential future CO 2 storage in the area.
Oxygen isotope analyses of diagenetic cements can provide detailed evidence of sedimentary burial processes and conditions, as the δ 18 O values of precipitating minerals reflect contemporaneous local δ 18 O water and temperature conditions. Uncertainties in the timing and rates of pore water δ 18 O evolution in sedimentary basins can complicate interpretation of these records. Fracture-bridging (0.5–1 mm) quartz cements observed in sandstones of the Cretaceous Travis Peak Formation in the East Texas basin show clear growth-zoning by cathodoluminescence and contain detailed fluid inclusion records of temperature that make them excellent candidates for interrogating prolonged histories of basin temperature and the evolution of δ 18 O in basin pore water. New secondary ion mass spectrometer (SIMS) δ 18 O quartz isotopic data from fluid inclusion-rich quartz bridges in Travis Peak sandstones record a steady increase of pore water δ 18 O values from ~5 to 7‰ (VSMOW; Vienna Standard Mean Ocean Water) as the sandstone warms from ~130 to 150 °C. To help evaluate whether this trend could be generated solely from local water-rock interactions in response to burial compaction, a one-dimensional closed system isotopic burial model was created to simulate how δ 18 O water values change in a quartz-dominated sandstone during diagenesis. Using both directly measured and inferred rates of Travis Peak compaction, the magnitude of change in δ 18 O water that we calculate from quartz bridge geochemistry cannot be reasonably modeled solely by local quartz mechanical compaction, pressure solution, and cementation processes, necessitating significant fluxes of silica and high-δ 18 O water from outside of the sandstones prior to maximum burial. In conclusion, this indicates that even units which appear surrounded by significant barriers to fluid flow (i.e., mudrock-bounded channel sandstones) may have been infiltrated and diagenetically modified by large fluxes of fluid on geologic time scales.
A Portsmouth On-Site Waste Disposal Facility (OSWDF) All-Pathways analysis has been conducted that considers the radiological impacts to a resident farmer. It is assumed that the resident farmer utilizes a farm pond contaminated by the OSWDF to irrigate a garden and pasture and water livestock from which food for the resident farmer is obtained, and that the farmer utilizes groundwater from the Berea sandstone aquifer for domestic purposes (i.e. drinking water and showering). As described by FBP 2014b the Hydrologic Evaluation of Landfill Performance (HELP) model (Schroeder et al. 1994) and the Surface Transport Over Multiple Phases (STOMP) model (White and Oostrom 2000, 2006) were used to model the flow and transport from the OSWDF to the Points of Assessment (POAs) associated with the 680-ft elevation sandstone layer (680 SSL) and the Berea sandstone aquifer. From this modeling the activity concentrations radionuclides were projected over time at the POAs. The activity concentrations were utilized as input to a GoldSimTM (GTG 2010) dose model, described herein, in order to project the dose to a resident farmer over time. A base case and five sensitivity cases were analyzed. The sensitivity cases included an evaluation of the impacts of using a conservative inventory, an uncased well to the Berea sandstone aquifer, a low waste zone uranium distribution coefficient (Kd), different transfer factors, and reference person exposure parameters (i.e. at 95 percentile). The maximum base case dose within the 1,000 year assessment period was projected to be 1.5E-14 mrem/yr, and the maximum base case dose at any time less than 10,000 years was projected to be 0.002 mrem/yr. The maximum projected dose of any sensitivity case was approximately 2.6 mrem/yr associated with the use of an uncased well to the Berea sandstone aquifer. This sensitivity case is considered very unlikely because it assumes leakage from the location of greatest concentration in the 680 SSL in to the Berea sandstone aquiver over time and does not conform to standard private water well construction practices. The bottom-line is that all predicted doses from the base case and five sensitivity cases fall well below the DOE all-pathways 25 mrem/yr Performance Objective.
To establish the potential of safe, long-term commercial-scale carbon capture & storage (CCS) in southwestern Wyoming, the University of Wyoming formed a CCS coordination team (CCT) to complete a Phase I CarbonSAFE pre-feasibility study adjacent to the region’s largest coal-fired power plant, the Jim Bridger Plant (JBP). The study’s primary objective was to evaluate the potential of storing 50+ million metric tons of anthropogenic carbon dioxide (CO 2 ) in stacked reservoirs over a period of 25 years. In addition, this integrated study evaluated related aspects of carbon capture utilization & storage (CCUS) at the study site, including: (1) the CO 2 character of post-combustion flue gas from JBP and the technical capacity for retrofitting proven commercial-scale CO 2 capture technology relative to CO 2 source character and technological constraints; (2) utilization of the existing CO 2 pipeline network in the immediate vicinity of JBP for enhanced oil recovery (CO 2 -EOR) and utilizing that network as part of a regional hub for CO 2 storage from other anthropogenic sources; (3) assessment of the challenges and benefits of meeting all of CarbonSAFE programmatic Phase I goals relative to Wyoming’s carbon management regulatory framework; and (4) assessment of the community and environmental factors that may impact CCS at JBP. With respect to storage, we analyzed the potential of stacked Mesozoic reservoirs at the study site, focusing on the Entrada and Nugget sandstones. The reservoir response to injection simulations varied greatly between formations, though they both were capable of CO 2 storage. The Entrada Sandstone is relative thin (~40’) and heterogeneous; storage using one injection and production well was limited to 1.2 Mt over 25 years at the study site. The Nugget Sandstone is much thicker (> 400’) and generally homogeneous; storage using one injection and production well approached 15.0 Mt over 25 years; this scenario suggests that the storage potential in the Nugget Sandstone is close to the estimated P90 calculation of 9.6 Mt/mi 2 as models indicate a storage potential of 8.3 Mt/mi 2 at saturation. Both reservoirs’ storage capacity benefited greatly from coupling an injection well with an offset production well; total storage capacity increased and reservoir pressure stayed below critical thresholds. This indicates that implementing a pressure management strategy at the study site would optimize storage and decrease risk. Seals associated with the targeted reservoirs were shown to be capable of retaining the potential storage volumes, and the site benefits from multiple seal redundancies (i.e. > 7,000’ of seal formations). Furthermore, assessment of different reservoir variables (e.g., fluid composition, pressure, etc.) indicates all reservoirs at the study site are isolated. Storage assessments from Mesozoic reservoirs were coupled with simulations from the site’s Paleozoic reservoirs (the Weber Sandstone and Madison Limestone), which were the focus of previous study. The fully integrated stacked reservoir pre-feasibility studies suggest that the Nugget and Madison formations alone could safely store 50 million metric tons of CO 2 within the study site. Furthermore, both formations are true saline reservoirs near JBP, making them ideal target reservoirs for CO 2 sequestration. Analysis of Wyoming’s laws and regulations confirm that the State has a favorable policy environment for CCS/CCUS. The project’s economic models estimated that the following revenues collectively are sufficient to finance the project over its lifespan if CO 2 utilization is coupled with CCS via the following strategies: (1) sales of CO 2 for EOR (approximately $69 million/year); (2) use of CO 2 tax credits such as amended §45Q and §48A (approximately $484 million total); and (3) sales of low-carbon electricity and marketable carbon offset/credits into carbon-constrained markets (approximately $11-$17 million/year) to West Coast consumers. These revenue estimates are broadly consistent with other studies that concluded arbitrage in JBP retail electricity sales between California and Wyoming could help to support the cost of deploying CCS at JBP. The models also indicate that maintaining a successful business case for long-term CCS -- wherein associated capture, transportation and operational costs are fully addressed at JBP -- is reliant on numerous external factors that complicate predictions. This study suggests that this site meets the CarbonSAFE program’s requirements of being able to feasibly store 50+ million metric tons of CO 2 over 25 years within the site’s stacked reservoirs, especially if coupled with pressure management strategies. This study suggests that the Nugget Sandstone and the Madison Limestone are two of the best reservoirs in the State with respect to overall storage capacity. In addition, the site’s proximity to existing CO 2 transportation networks, CO 2 -EOR opportunities, ability to sell to markets that value low-carbon electricity, Wyoming’s existing carbon regulatory framework, and a public that is well-educated with respect to energy markets and issues collectively increase the potential for implementing commercial-scale CCS adjacent to the JBP. With respect to the available technical and non-technical site data, and the reservoir quality at the study site, we suggest that this site is capable of moving directly into Phase III of the CarbonSAFE program.
For much of the Mesozoic an extensive desert environment extended across the Western Interior of the USA. These desert systems deposited vast dunefields, including one of the largest ergs preserved in the rock record. The Glen Canyon Group in Utah is a Triassic-Jurassic aged aeolian succession, consisting of thick, laterally extensive aeolian dune deposits, interdune lakes and oases, and dryland fluvial systems. These aeolian deposits have been identified as potential CO2 sequestration reservoirs. The Jurassic-aged Navajo Sandstone has been the subject of research into its potential as a CO2 reservoir, as well as naturally occurring CO2 seeps in the Green River area. It has excellent reservoir properties, consisting of thick sandstones with high porosity and permeability, and has industry data including well logs and core. It occurs in both outcrop and subcrop, allowing for comparison and sense checking of interpretations across small and large scales. The Triassic Wingate Sandstone has received much less attention, however it too is composed of thick aeolian dune deposits. Like the Navajo Sandstone, the Wingate Sandstone has industry data, and occurs in both outcrop and subcrop. Using extensive industry data which exists across Utah is an effective way to pivot from a hydrocarbon focus to CO2 injection, and so contribute to green energy and carbon neutral emission goals. Here we use historic industry data in conjunction with field work, to study the Navajo and Wingate Sandstones and their potential as CO2 reservoirs, in addition to increasing our understanding of lithologic and stratigraphic complexity in one of the most significant aeolian systems in the world.
The Beaverhead impact structure, in southwestern Montana, was identified recently by the presence of shatter cones and impactites in outcrops of Proterozoic sandstones of the Belt Supergroup. The cones occur over an area greater than 100 sq km. Because the geologic and tectonic history of this region is long and complex, the outline of the original impact crater is no longer identifiable. The extent of the area over which shatter cones occur suggests, however, that the feature may have been at least 60 km in diameter. The absence of shatter cones in younger sedimentary units suggests that the impact event occurred in late Precambrian or early Paleozoic time. We have collected samples of shocked sandstone from the so-called 'Main Site' of dark-matrix breccias, and of impact breccias and melts from the south end of Island Butte. The melts, occurring often as veins through brecciated sandstone, exhibit a distinctive fluidal texture, a greenish color, and a cryptocrystalline matrix, with small inclusions of deformed sandstone. Samples of the same type, along with country rock, were analyzed previously for major- and trace-element abundances. It was found that, although the major-element composition as relatively uniform, trace-element composition showed variations between the melt material and the adjacent sandstone. These variations were attributed to extensive weathering and hydrothermal alteration. In a more specific search for a possible meteoritic signature in the breccia and the melt material we have conducted a new series of trace-element analyses on powders of our own samples by thermal neutron activation analysis. Our results indicate that Ir abundances in the breccia, the melts, and the adjacent sandstone clasts are no greater than about 0.1 ppb, suggesting no Ir enrichment of the breccia or the melts relative to the country rock. However, both the breccia and the melt material exhibit notable enrichments in Cr (8- and 10-fold), in U (9- and 5-fold), and in the heavy REE's (1.5- and 3-fold), respectively.
The Knox Project objective is to evaluate the potential of formations within the Cambrian-Ordovician strata above the Mt. Simon Sandstone (St. Peter Sandstone and Potosi Dolomite) as potential targets for carbon dioxide (CO2) sequestration in the Illinois and Michigan Basins. The suitability of the St. Peter Sandstone and Potosi Dolomite to serve as reservoirs for CO2 sequestration is discussed in separate reports. In this report the data gathered from the Knox project, the Illinois Basin – Decatur Project (IBDP) and Illinois Industrial Carbon Capture and Sequestration project (IL-ICCS) are used to make some conclusions about the suitability of the Maquoketa shale as a confining layer for CO2 sequestration. These conclusions are then upscaled to basin-wide inferences based on regional knowledge. Data and interpretations (stratigraphic, petrophysical, fractures, geochemical, risk, seismic) applicable to the Maquoketa Shale from the above mentioned projects was inventoried and summarized. Based on the analysis of these data and interpretations, the Maquoketa Shale is considered to be an effective caprock for a CO2 injection project in either the Potosi Dolomite or St. Peter Sandstone because it has a suitable thickness (~200ft. ~61m), advantageous petrophysical properties (low effective porosity and low permeability), favorable geomechanical properties, an absence of observable fractures and is regionally extensive. Because it is unlikely that CO2 would migrate upward through the Maquoketa Shale, CO2, impact to above lying fresh water aquifers is unlikely. Furthermore, the observations indicate that CO2 injected into the St. Peter Sandstone or Potosi Dolomite may never even migrate up into the Maquoketa Shale at a high enough concentrations or pressure to threaten the integrity of the caprock. Site specific conclusions were reached by unifying the data and conclusions from the IBDP, ICCS and the Knox projects. In the Illinois Basin, as one looks further away from these sites, the formation characteristics are expected to vary. The degree of how well this data can be extrapolated throughout the Basins (regionalized) is difficult to quantify because of the limited amount of data collected on the Maquoketa Shale away from IBDP, IL-ICCS and the Knox projects. Data gathered from the IBDP/IL-ICCS/Knox projects were used to make conclusions about the suitability of the Maquoketa shale as a confining layer for CO2 sequestration. This study indicates that the Maquoketa Shale would be a suitable caprock for a CO2 injection program in either the Potosi Dolomite or St. Peter Sandstone.
Fracture growth processes under subsurface conditions are known to be affected by fluid chemical environment. Chemical effects on fracture growth are likely to be relevant during host rock and organic matter diagenesis and catagenesis, hydraulic fracturing, in systems undergoing enhanced oil recovery, in geothermal systems, and in CO 2 sequestration reservoirs. To address the chemical effects on opening-mode fracture growth under a wide range in fluid chemical conditions and for different rock types we conducted double-torsion fracture mechanics tests to determine mode-I fracture toughness (KIC), subcritical index (SCI), and the stress-intensity factor vs fracture velocity (K-V) behavior. We tested sandstones, shales and, altered crystalline rocks. Samples were tested in aqueous conditions with variable salinity, solution pH, and at temperatures of up to 70 °C. Fracture growth in shales is highly sensitive to the environmental conditions, with a reduction of KIC of up to 47% from ambient air to water-saturated conditions. The primary interactions between shales and reactive solutions are clay-water and carbonate dissolution. For clay-rich Woodford and Mancos shales, solutions with higher salinity reduce subcritical fracture growth, while acidic conditions favor fracture growth in carbonate-rich Marcellus shale. Increasing the solution temperature promotes fracture growth. In sandstone, chemical effects on fracture parameters are most pronounced for clay- and carbonate-rich samples. Both KIC and SCI are lower in water than in air for nearly all sandstone lithologies tested. SCI is further reduced in chlorite-cemented Tuscaloosa Sandstone in acidic solutions. Mount Simon Sandstone is also sensitive to pH, most significantly in diagenetically bleached samples. Compositionally mature Aztec Sandstone is affected by water saturation, but shows little sensitivity to chemical environment. Silicified gabbro shows a slight reduction in KIC but a >60% reduction in SCI in aqueous environments compared to ambient conditions. The reduction in SCI is most pronounced in alkaline and saline solutions, suggesting silica dissolution and cation exchange at the fracture tip facilitate subcritical fracture growth. Limited reduction in KIC is attributed to low permeability. Overall, our results demonstrate that fracture properties depend strongly on water saturation, on fluid composition and pH, and on temperature, with an increase in temperature lowering fracture resistance.
Geopressured geothermal reservoirs are characterized by high temperatures and high pressures with correspondingly large quantities of dissolved methane. Two methods utilizing reservoir simulation techniques have been used to estimate recoverability factors of geothermal brine and methane based on well log data from a specific reservoir in Texas. The first assumes a simplified reservoir that hasthree layers: upper shale, sandstone, and lower shale with the sandstone layer thickness equal to the net sandstone in the reservoir interval. The second method uses a detailed reservoir model that accounts for multiple sandstone and shale layering. This method includes 12 layers of sandstone or shale with the layer depths determined from a well log. These two methods are used to answer thequestion on the sensitivity of the results to the level of detail that is included in the reservoir model. It was found that incorporating multiple thin layers of lower permeability sandstone can noticeably impact the results of the reservoir simulation. The heterogeneous model resulted in greater flow rates of both geothermal brine and total methane. Both models demonstrate that thegeopressured geothermal reservoir is capable of producing hot geothermal fluid at flow rates over a long duration that are sufficient for electricity production from binary power plants. The results indicate that simplified models of geopressured geothermal reservoirs that approximate actual reservoir details can be applied to give a reasonable, albeit conservative, estimate of the recoverableresource over broad areas using generalized data sets.
The Wabash #1 well, drilled for the Wabash CarbonSAFE Project and located in Vigo County, Indiana, USA, was drilled in early 2020 as a stratigraphic test well to characterize and evaluate the basal Cambrian Mt. Simon Sandstone for carbon dioxide storage (TD=8750 ft; 2667 m). The Wabash #1 well is located along the eastern flank of a newly interpreted Cambrian aulacogen that occurs in western Indiana and eastern Illinois. Here we present 938 new detrital zircon U-Pb ages (LA-ICPMS) from early Cambrian sandstones sampled near the base of the well. A basalt lava flow was penetrated at ~8530 ft (2600 m) and has an 40 Ar/ 39 Ar age of 525.03 +/- 1.10 Ma, which represents the first known Cambrian crust in the Illinois Basin. The two sandstone samples from beneath the basalt are dominated by zircons derived from the Midcontinent Granite-Rhyolite terrane. The sandstone samples from above the basalt reflect a mixture of these locally derived Mazatzal and Granite-Rhyolite terrane zircons, but also distal Archean, Grenville, and Yavapai zircons. Each sample has small numbers of Cambrian zircons, which is consistent with those in basal Cambrian sandstones in other deep wells to the west. These early Cambrian detrital zircons and early Cambrian age of the basalt, combined with sediment thickness patterns permit the interpretation of the Illinois aulacogen, which formed during the final stage of Rodinian rifting.
We focus on potential reservoirs in the Logan Canyon Formation. Based on lithostratigraphy, Libby-French originally divided the Logan Canyon sandstones into lower and upper Logan Canyon sandstone units separated by the Sable Shale equivalent. Miller et al. reinterpreted the well-log correlations with biostratigraphy using a sequence stratigraphic approach and found three distinct Logan Canyon highstand sandstones associated with three Aptian to lower Cenomanian stratigraphic sequences, the LC3, LC2, and LC1, respectively. In addition, Miller et al. found that the Sable Shale Formation identified by Libby-French was placed above and below the LC1 basal sequence boundary in some wells, highlighting issues with correlating units based only on lithology. The Upper Cretaceous strata in the BCT are termed the Dawson Canyon Equivalent Formation, correlate to the DCx composite sequence, and cap the Logan Canyon sequences below made up of a very thick, calcareous Cretaceous sequence stratigraphy mudstone, divided by a single, thick, laterally continuous Coniacian to Campanian sandstone called the Middle Sandstone.
The presence of compliant pores in rocks is important for understanding the stress–strain behaviors under different stress conditions. This paper describes findings on the effect of compliant pores on the mechanical behavior of a reservoir sandstone under hydrostatic and triaxial compression. Laboratory experiments were conducted at reservoir temperature on Weber Sandstone samples from the Rock Springs Uplift, Wyoming. Each experiment was conducted at three sequential stages: (stage 1) increase in the confining pressure while maintaining the pore pressure, (stage 2) increase in the pore pressure while maintaining the confining pressure, and (stage 3) application of the deviatoric load to failure. The nonlinear pore pressure – volumetric strain relationship governed by compliant pores under low confining pressure changes to a linear behavior governed by stiff pores under higher confining pressure. The estimated compressibilities of the matrix material in sandstone samples are close to the typical compressibility of quartz. Because of the change in pore structures during stage 1 and stage 2 loadings, the estimated bulk compressibilities of the sandstone sample under the lowest confining pressure decrease with increasing differential pressure. The increase in crack initiation stress is limited with increasing differential pressure because of similar total crack length governed by initial compliant porosity in sandstone samples.