Skip to main content

Introduction

A preliminary assessment of the geologic CO2 sequestration potential was made for an area surrounding the LG&E-KU Cane Run Generation Station in Jefferson County, Kentucky. A circular area with a 20-km radius around the plant was defined as the primary focus of the evaluation, but data from beyond 20 km were also used because of limited data from the primary area (Figure 1).

The following data were compiled for the evaluation:

  1. Locations of all petroleum-exploration and waste-disposal wells penetrating the Cambrian-Ordovician Knox Group or deeper formations (Kentucky and Indiana Geological Surveys)
  2. Formation tops for geologic units from the top of the Ordovician to the Precambrian
  3. Available digital geophysical logs for Knox and deeper wells (Kentucky and Indiana Geological Surveys)
  4. Core analyses (porosity and permeability) for Mount Simon Sandstone and Eau Claire Formation

The data extracted or interpreted from these sources were input into National Energy Technology Laboratory’s CO2-SCREEN modeling software package for subsurface volumetric calculations (Haeri et al., 2024). See the Executive Summary for detailed modeling results.

Geologic Setting Surface Technology

The Cane Run Generation Station lies on the west flank of the Cincinnati Arch, a broad anticline that separates the deep Illinois Basin in western Kentucky from the Appalachian Basin in eastern Kentucky. The arch developed in Middle Ordovician time, and rock units deposited prior to this time have been tilted to the west toward the Illinois Basin. In general, subsurface data dips westward in the study area on the arch, deepening to the west and shallowing to the east.

Figure 1
Figure 1 - Project location map. The bold gray line indicates the Indiana (to NW) and Kentucky (to SE) state border, counties are outlined in green, mapped surface faults are in red, and the blue circle represents a 12.4-mi (20-km) radius around the Cane Run Generation Station. Locations for oil and gas wells that were used for data analysis are marked by their completion code symbols. Location of shotpoints along proprietary seismic line XM-7 displayed in purple.

The Cane Run facility is located on unconsolidated sediments deposited along the Ohio River in the West Louisville 7.5’ Quadrangle (Kepferle, 1974). These sediments are Quaternary (Pleistocene) in age and are interpreted to be glacial outwash deposits (Wisconsin age). The Quaternary outwash was deposited unconformably onto the Mississippian bedrock (siltstones and shales of the Borden Formation), which is exposed in the hills and bluffs to the SSE of the station (Figure 2).

Figure 2
Figure 2 - Surface geologic map of area surrounding Cane Run Station (red spot). Yellow areas are Quaternary alluvium and loess, green and orange areas are Mississippian Borden Fm., and brown areas are artificial fill.

Surface geology does not have a direct impact on carbon sequestration potential, since CO2 injection will occur at much deeper depths. However, the abundance of low-permeability shales in the near-surface Upper Devonian and Lower Mississippian rocks could serve as secondary confining layers (where unfractured) in the unlikely event CO2 was to migrate through the deeper primary seals.

Only one pre-existing well penetrates the reservoir interval within a 12.4-mi (20-km) radius of the Cane Run facility, the nearby (~ 3.5 mi NE) DuPont #1WAD DuPont Fee well (Figure 1). Fortunately, extensive data exists for this 1971 waste-disposal well, including high-quality geophysical logs and lab analyses from core samples extracted during drilling. Although the mapped surfaces in this report were 8 produced using all the available wells in the region (~ 80-mi radius around Cane Run), the reservoir properties utilized in modeling were all based upon the DuPont #1WAD well.

Stratigraphy and Structure

In this part of the central US, efficient geologic storage of CO2 is generally confined to depths greater than 2,600 ft below the surface, where the temperatures and pressures are such that CO2 should be in the supercritical, or dense, phase (Frailey et al., 2005; Wickstrom et al., 2003; Parris et al., 2010). Supercritical CO2 has properties of both a liquid and a gas but much higher density than gaseous CO2, which dramatically reduces the pore volume needed to contain the injected gas (i.e., more CO2 fits into a smaller volume). In western Jefferson County, this 2,600-ft depth would likely occur within the Cambrian Eminence Dolomite of the lower Knox Supergroup. Geologic formations in this area below 2,600 ft include the basal part of the Knox, the Upper/Middle Cambrian Eau Claire Formation, the Middle Cambrian Mount Simon Sandstone, and Precambrian Eastern Granite-Rhyolite Province basement (Figure 2). These formations are briefly described below, from oldest to youngest.

Precambrian Basement (Low-Permeability Lower Sealing Unit)

The Precambrian basement in the study area consists of mostly igneous rocks (felsic plutonic rocks with scattered mafic intrusions) assigned to the Eastern Granite-Rhyolite Province (1.42–1.50 Ga) (Bickford et al., 1986; Van Schmus et al., 1996). However, some wells in the region have penetrated clastic deposits of the 1.0-Ga-old Middle Run sandstones (Drahovzal et al., 1992). The Middle Run was deposited in nonmarine fluvial environments and consists of fine-grained, red lithic sandstones and minor siltstone and shale (Drahovzal et al., 1992). However, due to about 1.0 billion years of compaction and diagenesis, the sandstone (where present) is well cemented and lacks any measurable porosity or permeability and therefore is unlikely to have any potential for carbon storage in the study area. In practical terms, regardless of whether the Precambrian lithology is a lithic arenite or igneous crystalline rock, it forms the lower confining layer for potential injection into the overlying Mount Simon Sandstone.

Precambrian strata are overlain by Cambrian strata above an erosional surface called the Precambrian unconformity. The Precambrian unconformity surface dips to the southwest in the study area (~ 0.7° at an azimuth of 250°), consistent with the trend of the Cincinnati Arch (Figure 3). This structure map is based on the few wells that penetrate the Precambrian surface in the area. As such, it should be considered a general representation of the structure of the area. Drilling would be needed to confirm the actual depth to Precambrian basement. The estimated depth to basement at Cane Run is approximately 5,950 ft.

Figure 3
Figure 3 - Geophysical well logs of the DuPont #1WAD DuPont well in Jefferson County, Kentucky in the left column, the green curve represents passive gamma-ray response (both increase to the right), and the orange curve is the caliper reading. In the right column, the red curve is bulk density log (increases to the right), and the blue curve is the neutron-porosity log (decreases to the right).

 

Figure 4
Figure 4 - Structure map of the top of the Precambrian basement. Elevations are in feet above/below sea level

Cambrian Mount Simon Sandstone (Primary Injection Zone)

The Cambrian Mount Simon Sandstone unconformably overlies the Precambrian basement in most of the study area. The quartz-rich Mount Simon Sandstone is the primary CO2 injection zone in the study because of its depth and porosity. The Mount Simon is the primary reservoir target in much of the Illinois Basin to the west (Frailey et al., 2005; Leetaru and McBride, 2009; Gray, 2015; Greenberg, 2017) and was the target of a carbon storage demonstration well to the east at Duke Energy’s East Bend Generating Station (Battelle, 2011; Sminchak, 2012). The Mount Simon has been encountered in five wells in the study area. Cores from the Mount Simon Sandstone are available from two of these wells: the Battelle Duke Energy well and the DuPont #1WAD waste-injection well in Louisville. Porosity and permeability data measured in these cores are described further in the Reservoir Quality section.

Using well data from the area, structure and thickness maps for the Mount Simon were constructed. Other studies have used data from seismic lines outside this study area to map the extent of the Mount Simon Sandstone across Kentucky. The broader regional data show the Mount Simon thickens to the north and northwest and pinches out toward the southeast (Figure 4). The Mount Simon is known to be absent in several wells in central Kentucky, but the exact southern limit of the pinchout is uncertain, and the limit shown in the map should be considered a preliminary boundary that may be revised with more data in the future.

Figure 5
Figure 5 - Thickness map of the Mount Simon Sandstone in Kentucky. Bold lines are interpreted basement faults, and thin lines represent county boundaries. Modified from Greb and Solis, 2009. Blue circle indicates the 20-km radius around the Cane Run Facility, defining the project study area.

The top of the Mount Simon is at 5,103-ft depth in the DuPont well in Louisville (Figure 5), and ranges in thickness from ~ 400 to 750 ft across the radius of the study area (Figure 4). Approximately 78 mi to the northeast in Boone County, Kentucky, Duke Energy was able to successfully inject 1,000 tons of CO2 into the Mount Simon Sandstone at their East Bend facility in 2009 (Battelle, 2011; Dillon et al., 2020). At the Cane Run Facility, the top of the Mount Simon is estimated to be at a subsea elevation of 12 approximately −4,770 ft with an estimated thickness of approximately 820 ft. Drilling would be needed to confirm thicknesses and depths.

Figure 6
Figure 6 - Structure on the top of the Cambrian Mount Simon Sandstone. Elevations are in ft above/below sea level.

The Eau Claire Formation directly overlies the Mount Simon Sandstone. The Eau Claire Formation was cored in the DuPont #1WAD well, at depths ranging from 4,409 to 4,459 ft and 4,842 to 4,871 ft. In core, the Eau Claire is composed of green and gray marine shale, thin siltstones, and interbedded dolomite. Based on regional data, the Eau Claire should be 600 to 700 ft thick (Figure 6) in the study area, and the top of the unit should be from approximately 3,700 to 4,700 ft below sea level (Figure 7). The Eau Claire deepens to the southwest into the southern Illinois Basin and the Rough Creek Graben (Hickman, 2011). The top of Eau Claire Formation at Cane Run is at a subsea elevation of approximately -4,160 ft at Cane Run.

Thick shales in the Eau Claire Formation are generally considered the primary confining layer (seal) for gas and liquid injection into underlying Mount Simon Sandstone reservoirs in the region (Frailey et al., 2005; Wickstrom et al., 2005; Leetaru and McBride, 2009; Battelle, 2011). Where the formation contains thick shales, shaly facies commonly exhibit permeabilities of 10-5 mD (Sminchak, 2012), which is adequate for the unit to be considered a primary confining layer (seal). The top of this confining layer is 2,000 ft deeper than the minimum depth required for stable supercritical CO2 conditions in the study area.

Figure 7
Figure 7 - Thickness of the Cambrian Eau Claire Formation, in ft.
Figure 8
Figure 8 - Structure of the top of the Cambrian Eau Claire Formation. Elevations are in ft above/below sea level.

Cambrian-Ordovician Knox Supergroup

The Knox Supergroup overlies the Eau Claire Formation. In this part of the Illinois Basin, the Knox is divided into a series of dolomite units, from top to bottom: the Shakopee Dolomite, Oneota Dolomite, and Potosi Dolomite (Treworgy and Whitaker, 1990; Greb and Solis, 2009; Greb, 2017). A sandstone or quartzose dolomite unit equivalent to the Gunter Sandstone or an unnamed Shakopee sandstone occurs locally at or near the base of the Shakopee Dolomite but is generally thin. This sandstone was 16 encountered in the East Bend Duke Energy well (Battelle, 2011). The top of the Knox is a regional erosional unconformity that formed when the region was uplifted above sea level during the Early Ordovician. Hence, the upper part of the Shakopee Dolomite is progressively truncated eastward in the study area. The Knox is approximately 2,800 ft thick in the study area (Figure 10).

Knox dolomites are commonly tight, low-permeability units. However, they may also contain thin intervals with fracture and/or vuggy porosity which can act as local reservoirs (Greb and Solis, 2009; Lasemi and Askari, 2024; Leetaru et al., 2014). The Knox is used in other parts of the basin for brine disposal (Keller and Abdulkareem, 1980; Greb et al., 2009). Also, sandstones in the Knox can have reservoir potential. The Gunter or Shakopee sandstone near the base of the Shakopee Dolomite was successfully tested for CO2 injection at the KGS No. 1 Blan well in Hancock County, Kentucky, approximately 52 mi west of Cane Run (Bowersox et al., 2016).

Porous zones in the Knox have also been used for natural-gas storage by LG&E near the study area, in the Ballardsville and Eagle Creek storage fields in Grant and Oldham Counties, respectively (Greb and Solis, 2009). These storage fields are now abandoned but are too shallow for CO2 storage.

In the study area, the upper limit for CO2 in a supercritical phase (2,600-ft depth) lies within the Shakopee Dolomite of the Knox Supergroup. The unit would have to be cored and tested to determine if it would be a secondary confining interval above the Eau Claire or if it has local reservoirs. For this preliminary assessment, we did not consider the lower part of the Knox (below 2,600-ft depth) to be a viable injection target, because there is not a thick shale that would be considered a reliable sealing interval (containment zone) above the Eau Claire that is also below the 2,600-ft depth limit (although low-permeability Knox carbonates could offer secondary confining properties).

Because the elevations at the top of the Knox Supergroup are well above the minimum depth limit, it is the shallowest interval mapped in this evaluation. Many more wells have been drilled to the top of the Knox than to the deeper horizons, and thus more data are available for the Knox structure map (Figure 8). The Knox dips to the west, with the projected top of the Knox at about −1,350 ft subsea (1,770-ft depth) at the Cane Run facility.

The Knox isopach map (Figure 9) shows that the unit thins by more than 300 ft from southwest to northeast across the radius of the study-area. This thinning is primarily caused by erosional truncation at the top of the Knox section (sub-St. Peter Sandstone unconformity). This thinning is also illustrated in the regional cross section (Figure 11). The Knox Supergroup is interpreted to be around 2,820 ft thick at Cane Run.

Near-surface Formations

A series of Ordovician through Mississippian strata lie above the Knox Supergroup rocks at Cane Run. The Ordovician Black River Group (Kepferle, 1974) consists of limestones, minor dolomite, and interbedded shales that generally have low permeability (unless fractured) and could form a shallow secondary confining zone for CO2 injected into the deeper Mount Simon Sandstone (Harris and Hickman, 2013). However, all are above the 2,600-ft CO2 storage depth. In addition, above the Lexington Limestone lies ~ 250 ft of Upper Ordovician shale and 37 ft of Devonian shale (Figure 11). These clayrich shales have low porosity and permeability, creating shallow, tertiary seals for any potentially leaked CO2, although shallow units may fall within the zone of near-surface fracturing.

Deep Faults and Available Seismic Data

The only seismic data in the Kentucky Geological Survey inventory for the study area is one proprietary line running east-west in Floyd and Harrison Counties, Indiana, which was acquired by Vastar Energy in 1994 (Figure 1). This line shows flat-lying strata with no resolvable faults in the study area around Cane Run. Three faults have been mapped at the surface within the radius of the study area: one in Bullitt County, Kentucky, one in Floyd County, Indiana, and one in Harrison County, Indiana (Figure 10, Figure 11). All three of these unnamed faults are located near the edge of the study area and are all less than 4.6 mi in length. Because of these short lateral lengths, it is highly unlikely that their fault offsets or deformations penetrate either the seal (Eau Claire Formation shale) or the intended reservoir interval (Mount Simon Sandstone). Seismic data would likely have to be obtained at any future injection site to confirm that there are not unmapped faults in the study area.

Figure 9
Figure 9 - Structure on the top of the Cambrian-Ordovician Knox Supergroup. Elevations in ft above/below sea level.
Figure 10
Figure 10 - Thickness of the Cambrian-Ordovician Knox Supergroup, in ft.
Figure 11
Figure 11 - Locations of wells utilized in the cross sections in Figures 12 and 13. Red lines represent faults mapped at the surface, green lines are county boundaries, and the bold black line indicates the state border between Kentucky (SE) and Indiana (NW).
Figure 12
Figure 12 - South to north well-based diagram (lateral spacing not to scale). The Mount Simon deepens to the south. Cane Run Station is 3.5 mi south of the DuPont #1WAD well. Well logs include the gamma-ray in the left track, and density and neutron-porosity logs in the right track. Depth scales along edges are in ft above/below sea level. Bold red line represents the upper limit for stable CO2 in a supercritical phase.

Reservoir Quality and Injection Zone Thickness

In order to estimate an initial potential carbon sequestration capacity for the study area, average porosity and thickness of the storage zone need to be estimated. Reasonable estimates for porosity and net injection zone thickness were calculated from nearby wells. Data from the DuPont #1WAD well were especially helpful, since quality well logs and core data are available from this well, which was drilled just 3.5 mi northeast of the Cane Run facility.

The reservoir zone (Mount Simon Sandstone) below the Cane Run facility is estimated to be about 847 ft thick, laying at a depth of 5,088–5,935 ft below the surface (Figure 12, Figure 13). On the eastern edge of the Illinois Basin, along the Ohio River, the Mount Simon directly overlies the Precambrian Basement, (e.g., Greb and Solis, 2009). Although the entire Mount Simon 843-ft-thick interval appears to be mostly sandstone on available well logs, minor variabilities in porosity and permeability (related to grain size and facies changes) can exist at a local scale and are evident in the geophysical well logs acquired from the DuPont well. As in many sandstones, in addition to porosity variations due to local facies changes, porosity in the Mount Simon Sandstone also decreases with increasing burial depth. This is primarily because of cementation and compaction and is a result of increased temperature, pressure, and the amount of time the rocks have been buried. A substantial set of Mount Simon porosity and permeability data from across the Midwest was published by Medina et al. (2011), in which they reported a dramatic decrease in porosity at depths greater than 7,000 ft. Fortunately for the Cane Run location, the entire reservoir is at depths less than 6,000 ft, so these effects should be minimal. Site-specific data would be needed to confirm porosity and permeability for the reservoir.

This reservoir is capped by an estimated 589 ft of low permeability, clay-rich marine shale and siltstones of the Eau Claire Formation (Figure 12). Although only sampling 34 ft within the upper Eau Claire Formation interval, all 36 core plugs in this unit resulted in < 0.1 millidarcies of vertical permeability, indicating quality sealing properties (KGS Well Record #11169). The high percentage of clays in the Eau Claire Formation also make through-going fractures (breaching the seal) highly unlikely, so the ~ 70% seal-to-reservoir thickness ratio estimated to be present at the Cane Run location based on available data should be more than sufficient for an adequate seal of a potential Mount Simon Sandstone CO2 reservoir. Site-specific data would be needed to confirm the Eau Claire’s confining properties.

Figure 13
Figure 13 - South to north cross section from Figure 11, flattened on the top of the Eau Claire Formation (primary reservoir seal) to illustrate thickening of Mt. Simon to the north and west. The column below the blue star represents the depths and thicknesses at the location of the Cane Run facility. In the left tracks of the above logs are gamma ray responses (green) and caliper curves (orange). In the right tracks are bulk density (red), neutron porosity (blue), sonic velocity (orange), and photoelectric effect (grey).

Core measurements are the most accurate method of determining porosity and permeability for a subsurface rock unit, but geophysical porosity logs can also provide a continuous estimate of porosity for a unit. Then, permeability can be estimated from log porosity based on porosity vs. permeability relationships developed from known core data from the Mount Simon Sandstone. Core-derived porosity and permeability data for the Mount Simon are available from cores of the DuPont #1WAD well, but unfortunately, only 35 ft of the 843-ft-thick Mount Simon were cored by DuPont in 1971. Therefore, both well log and core porosity data were used to estimate porosity at Cane Run.

The log-based average porosity of the Mount Simon was calculated using both bulk-density logs and neutron-porosity logs. After correcting for a sandstone lithology with an elevated percentage of potassium feldspar silt content found in the Mount Simon (which reduces bulk-density measurements and elevates gamma-ray responses relative to a pure quartz sandstone), porosities from each method were calculated for every 0.5 ft of the Mount Simon interval using IHS Petra petrophysical software. This resulted in an average porosity across the 847-ft-thick interval of 7.4%, with the majority of the results within range of 3–12%. DuPont drilled one core in the Mount Simon interval (Core #15, 5,718–5,748-ft depth). There is some variation in the porosity measurements from the #1WAD core (< 3 to 9.4%), but they are within the same range and appear to agree with the calculated values from the log responses. Therefore, for the initial capacity calculations (Table 1) made herein, an average porosity of 7.4% with a standard deviation of 3.5% was used.

Table 1 - Porosity and permeability results from DuPont #1WAD well, Core #15. Lab data from Oilfield Research, Inc., April 10, 1971., Kentucky Geological Survey, Well Record # 11169
DuPont Sampe #Depth (ft)Permeability (md)Porosity
HorzVert
1385719.54.407.307.4
1395721.54.005.609.4
1405723.5<0.10<0.104.4
1415725.50.14<0.105.3
1425727.5<0.10<0.10<3.0
1435729.50.240.113.7
1445731.53.900.12<3.0
1455733.50.221.603.2
1465735.535.004.306.0
1475737.541.001.305.0
1485743.5<0.10<0.104.9
1495745.5<0.10<0.10<3.0
1505747.5<0.10<0.10<3.0

CO2 Capacity Calculations

Carbon dioxide capacity calculations in this study utilize the methods described by Goodman et al. (2016). Specifically, the software CO2-SCREEN (v.5.0) from the US Department of Energy’s National Energy Technology Laboratory was used.

The potential Mount Simon reservoir is not confined laterally within the study area and so is considered an “open system,” in that injection pressure can dissipate as injected liquid migrates away from the borehole. In a gross sense, the unit is fairly homogenous with no dramatic lithology changes or internal baffles, so it was treated as a single volume. No oil shows have been reported within the Mount Simon in the region, so the pore fluids within Mount Simon porosity will be modeled as saline brine. Therefore, in CO2-SCREEN, an “open system, saline 1-grid” model was run. Analyses of the Mount Simon from deeper in the Illinois Basin to the west (Freiburg et al., 2022) indicate that the Mount Simon deposition occurred in both shallow marine (beach) and fluvial (river) environments. Because of this, calculations within CO2-SCREEN were performed for both modeled situations. See Tables 1 and 2 for more information.

A center depth of 5,500 ft was used to calculate physical parameters for the Mount Simon. Using geothermal and pressure gradients appropriate for central Kentucky (Bowersox et al., 2016), a reservoir temperature of 118°F (48°C) and a reservoir pressure of 3,190 psi (22 MPa) were calculated. Due to variations in depth, lithology, etc., for the model inputs, we estimated standard deviations of 9°F (5°C) and 290 psi (2 MPa), respectively.

Volumetric calculations were performed using CO2-SCREEN for two different areal extents (local and regional scales) and five different injection durations (1, 5, 10, 20, and 30 years). For the regional scale evaluation, a study area with a 12.4-mi (20-km) radius centered on the Cane Run facility was used (see blue circle in Figure 1). Because of the larger study area (485 mi2 ), the anticipated uncertainty in the reservoir gross thickness and net-to-total area values are higher than for the smaller, local area assessment, but all other input values remained equal (Table 1). Note that the larger regional area includes parts of two Kentucky counties and three Indiana counties, possibly complicating pore-space ownership issues. This study did not evaluate any legal or financial aspects of utilizing this site for either a local or regional scale scenario.

The local area assessment estimated the capacity for carbon sequestration within the Cane Run Station’s property footprint. For its calculations, the CO2-SCREEN software assumes a simple, right cylinder reservoir shape. To approximate such a volume within the boundaries of the Cane Run property, a circular 26 area with a radius of ~ 2,297 ft (0.7 km) was defined to maximize the useable area without crossing surface-owner property, county, or state boundaries (Figure 14).

Figure 14
Figure 14 - Local area surrounding Cane Run. The yellow line represents the 2,300 ft (0.7-km) radius that was used to define the local study boundary. Map image from Google Earth.

The storage capacity equations used in CO2-SCREEN include efficiency factors, which reduce the CO2 storage capacity. This factor is applied because the available pore volume is never completely saturated with CO2 because of fluid characteristics and geologic variability within the reservoir. The software uses P10 and P90 values to define the statistical envelope of values for a Monte Carlo simulation. The P10 and P90 values are entered as a value between zero and one (where 0.01 represents 0%, and 0.99 represents a 100% value). Although there is limited data available at reservoir depths, geologically realistic values were interpreted for this interval and are displayed in Table 2. The application of efficiency factors 27 significantly reduces the calculated storage capacities but is necessary to determine reasonable volume estimates. The results for multiple scenarios are displayed in Figure 15, Figure 16, Table 3 and in the Executive Summary. Values range from only 670,000 tons of CO2 (P10 after 1 year of injection into shallow marine sandstones under the Cane Run facility’s 0.6-mi2 (1.54-km2) property footprint) to 4,616,000,000 tons (P90 after 30 years of injection into fluvial sandstones with 485 mi2 (1,256 km2). This range of nearly four orders of magnitude in predicted storage resource volumes is partly due to the lack of deep well data needed to accurately model the stratigraphic variability of the Mount Simon reservoir.

The potential for induced seismicity from large-scale injections must also be considered in estimations of carbon storage potential in Mount Simon sandstone reservoirs in the region. Induced seismicity is a concern in any carbon storage project. In recent Illinois Basin Mount Simon injections near Decatur, Illinois, the potential for induced seismicity was mitigated by not using the lower Mount Simon Sandstone that was in direct contact with underlying Precambrian rocks as a reservoir (Williams-Stroud et al., 2020). Injecting into only the upper part of the Mount Simon at the Cane Run study site would decrease the overall potential storage volume, but may reduce the risk of induced seismicity. Acquiring a 3D seismic survey of the reservoir at and around Cane Run to identify possible unknown and unmapped subsurface faults is recommended prior to industrial-scale injection activities to alleviate the concern for induced seismicity.

Figure 15
Figure 15 - Modeled cumulative injection volume using the "local" study area scenario.
Figure 16
Figure 16 - Modeled cumulative injection volume using the "regional" study area scenario.

Conclusions

Cane Run Station has potential for geologic storage of CO2 beneath the property. The Mount Simon Sandstone is likely the only formation beneath the surface in the study area that would have suitable porosity and permeability at the depths required for supercritical-phase sequestration. Excellent confinement for injected CO2 would likely be provided by the Eau Claire Formation, which is more than 500 ft thick and contains thick shale intervals in available downhole logs in the region. However, local geologic data control for the reservoir is limited, with only one well penetrating the reservoir interval within a 20-km radius (the DuPont #1WAD well). Because of the limited data, the estimates provided herein should be regarded as very preliminary. Site-specific downhole well data would be needed for more accurate assessments of the reservoirs and seals.

Volumetric statistics calculated by the CO2-SCREEN software for the two Cane Run Station study areas are listed in Table 2. Although the Mount Simon thins to the northeast towards Cincinnati (Figure 4), the unit is also closer to the surface at the top of the Cincinnati Arch. There is less overburden compressing 30 the reservoir at shallow depths, so average porosities will be higher, increasing injection potential. Because of this aspect, it may be possible to move upstream along the Ohio River towards Cincinnati and find alternate sites to drill injection wells down to the Mount Simon that would result in increased storage volumes from the increased porosity, despite the thinner interval of Mount Simon located in the region.

References

Battelle, 2011, CO2 injection test in the Cambrian-age Mt. Simon formation, Duke Energy East Bend Generating Station, Boone County, Kentucky: Columbus, OH, Battelle, final report to the U.S. Department of Energy, DOE-NETL cooperative agreement DE-FC26-05NT42589, 890 p.

Bickford, M.E., Van Schmus, W.R., and Zietz, I., 1986, Proterozoic history of the Midcontinent region of North America: Geology, v. 14, no. 6, p. 492–496.

Bowersox, J.R., Williams, D.A., and Harris, D.C., 2016, Phase 1 geologic evaluation of the Kentucky Geological Survey Marvin Blan No. 1 deep saline reservoir CO2 injection test well, Hancock County, Kentucky: Kentucky Geological Survey, ser. 12, Contract Report 63, 58 p.

Dillon, D., Chu, R., Lu, H., Freeman, B., Elliot, W., and McKaskle, R., 2020, Initial engineering design of a post-combustion CO2 capture (PCC) system for Duke Energy’s East Bend Station using membrane-based technology: EPRI, final report to the U.S. Department of Energy, DOE-EPRI31589, 274 p.

Drahovzal, J.A., Harris, D.C., Wickstrom, L.H., Walker, D., Baranoski, M.T., Keith, B.D., and Furer, L.C., 1992, The East Continent Rift Basin: A new discovery: Kentucky Geological Survey, ser. 11, Special Publication 18, 25 p.

Frailey, S.M., Leetaru, H.E., Finley, R.J., Gustison, S.R., Korose, C.P., Garner, D.A., Rupp, J., and Drahovzal, J., 2005, An assessment of geologic sequestration options in the Illinois Basin: University of Illinois, Illinois State Geological Survey, contract report to the U.S. Department of Energy, DE–FC26–03NT41994, 477 p.

Freiburg, J.T., Delpomdor, F., McBride, J.H., Leetaru, H.E., Malone, D.H., Grigsby, N., Khosravi, M., Frailey, S., Damico, J.R., Monson, C., and Askari, Z., 2022, The geology of the Mt. Simon sandstone storage complex at the Wabash #1 well, Vigo Co., Indiana: University of Illinois, 31 Urbana-Champaign, technical report to the U.S. Department of Energy, DOE-FE0031626‐9, 192 p.

Goodman, A., Sanguinito, S., and Levine, J.S., 2016, Prospective CO2 saline resource estimation methodology: Refinement of existing US-DOE-NETL methods based on data availability: International Journal of Greenhouse Gas Control, v. 54, p. 242–249.

Gray, K., 2015, Carbon storage atlas [5th ed.]: Peachtree Corners, GA, Southern States Energy Board, DOE-SSEB-42590-120, 114 p., doi:10.2172/1814017.

Greb, S.F., 2017, Kentucky stratigraphy with stage correlations: Kentucky Geological Survey, ser. 12, Map and Chart 210, 1 sheet. Greb, S., and Solis, M.P., 2009, Geologic carbon storage (sequestration) potential in Kentucky, chap. 4 of Parris, T.M., Greb, S.F., and Nuttall, B.C., Evaluation of geologic CO2 sequestration potential and CO2 enhanced oil recovery in Kentucky: Kentucky Geological Survey, ser. 12, Report of Investigations 21, p. 55–212.

Greb, S. F., D. C. Harris, M. P. Solis, W. H. Anderson, J. A. Drahovzal, B. C. Nuttall, R. A. Riley, W. Solano-Acosta, J. A. Rupp, and N. Gupta, 2009, Cambrian–Ordovician Knox carbonate section as integrated reservoirs and seals for carbon sequestration in the eastern mid-continent United States: in M. Grobe, J. C. Pashin, and R. L. Dodge, eds., Carbon dioxide sequestration in geological media—State of the science: AAPG Studies in Geology 59, p. 241 – 259.

Greenberg, S.E., Bauer, R., Will, R., Locke II, R., Carney, M., Leetaru, H., and Medler, J., 2017, Geologic carbon storage at a one million tonne demonstration project: Lessons learned from the Illinois Basin–Decatur Project: Energy Procedia, v. 114, p. 5529–5539. Haeri, F., Goodman, A., and Myshakin, E.M., 2024, CO2 storage prospective resource estimation excel analysis (CO2-SCREEN) user’s manual: Python_V5.0: Pittsburgh, PA, National Energy Technology Laboratory, DOE/NETL-2024/4801, 40 p., doi:10.2172/2331225.

Harris, D.C., and Hickman, J.B., 2013, Evaluation of geologic CO2 storage potential at LG&E and Kentucky Utilities power plant locations, central and western Kentucky: Kentucky Geological Survey, ser. 12, Contract Report 1, 106 p.

Hickman, J.B., 2011, Structural evolution of an intracratonic rift system; Mississippi Valley Graben, Rough Creek Graben, and Rome Trough of Kentucky, USA: University of Kentucky Doctoral 32 Dissertations, no.144, 186 p., 31 pl., https://uknowledge.uky.edu/gradschool_diss/144 [accessed 06/17/2024].

Keller, S. and Abdulkareem, T.F., 1980. Post-Knox Unconformity-Significance at Unionport Gas-Storage Project and Relationship to Petroleum Exploration in Indiana: Indiana Geological Survey.

Kentucky Geological Survey, Well Record #11169, https://kgs.uky.edu/OG_images/0/0/0/1/1/R00011169/R00011169.pdf [accessed 5/01/2024].

Kepferle, R.C., 1974, Geologic map of parts of the Louisville West and Lanesville Quadrangles, Jefferson County, Kentucky: U.S. Geological Survey Geologic Quadrangle Map 1202, scale 1:24,000, 1 sheet.

Lasemi, Y. and Askari, Z., 2024, Multiple Paleokarst Events in the Cambrian Potosi Dolomite, Illinois Basin: The Sedimentary Record, 22(1).

Leetaru, H.E., and McBride, J.H., 2009, Reservoir uncertainty, Precambrian topography, and carbon sequestration in the Mt. Simon Sandstone, Illinois Basin: Environmental Geosciences, v. 16, no. 4, p. 235–243. Leetaru, H.E., Smith, V., Will, R., Freiburg, J. and Brown, A.L., 2014, The application of an integrated earth model in reservoir management of a CO2 plume: Energy Procedia, 63, pp.2903-2910.

Medina, C.R., Rupp, J.A., and Barnes, D.A., 2011, Effects of reduction in porosity and permeability with depth on storage capacity and injectivity in deep saline aquifers: A case study from the Mount Simon Sandstone aquifer: International Journal of Greenhouse Gas Control, v. 5, p. 146–156.

Parris, T.M.; Greb, S.F.; and Nuttall, B.C., 2010, Evaluation of Geologic CO2 Sequestration Potential and CO2 Enhanced Oil Recovery in Kentucky: Kentucky Geological Survey, Report of Investigations 23. https://uknowledge.uky.edu/kgs_ri/23

Sminchak, J., 2012. Simulation framework for regional geologic CO2 storage along arches province of midwestern United States. Battelle Memorial Institute, Columbus, OH (United States).

Treworgy, J.D., and Whitaker, S.T., 1990, Stratigraphic and structural framework along the 3 o'clock cross section in the Illinois Basin: Wayne County, Illinois, to Switzerland County, Indiana: Illinois State Geological Survey, Open File Series 1990-03, 1 sheet.

Van Schmus, W.R., Bickford, M.E., and Turek, A., 1996, Proterozoic geology of the east-central midcontinent basement, in van der Pluijm, B.A., and Catacosinos, P.A., eds., Basement and basins of eastern North America: Boulder, CO, Geological Society of America, Special Paper 308, p. 7–23.

Wickstrom, L.H., McDonald, J., Riley, R.A., Carr, T.R., Nuttall, B. and Rupp, J.A., 2003, May. Online Tools to Evaluate Oil and Gas Fields for CO2 Sequestration. In AAPG Annual Meeting 2003.

Williams-Stroud, S., Bauer, R., Leetaru, H., Volker, O., Stanek, F., Greenberg, S., and Langet, N., 2020, Analysis of microseismicity and reactivated fault size to assess the potential for felt events by CO2 injection in the Illinois Basin: Bulletin of the Seismological Society of America, v. 110, p. 2188–2204.

Executive Summary - Geologic CO2 Sequestration

Power Plant: Cane Run

County: Jefferson, Kentucky

Geologic Basin: Cincinnati Arch

 

Data Quality:

Distance to nearest well control in reservoir: 3.5 mi

Wells to primary injection zone within 12.4-mi radius: 1

Distance to nearest core in injection zone: 3.5 mi

 

Reservoirs:

Primary injection zone: Cambrian Mount Simon Sandstone

Rock lithology: sandstone (quartz arenite)

Drilling depth at plant site: 5,088 ft

Trapping mechanism: regional dip (capillary and solution trapping)

Maximum reservoir pressure: 3,190 psi

Reservoir temperature: 48°C

Salinity of reservoir fluid: 200,000 ppm (estimated)

Reservoir thickness: 847 ft

Average porosity: 7.4 percent

Secondary injection zone: none at this site

 

Confinement and Integrity:

Primary confining zone: Cambrian Eau Claire Shale

Rock type: shale and dolomitic siltstone

Thickness of primary confining zone: 589 ft

Height above primary injection zone: 0 (overlies injection zone)

Penetrations of primary seal within 12.4-mi radius: 1

Secondary confining zone: Ordovician Black River Group

Rock lithology: micritic limestone

Thickness of secondary confining zone: 604 ft

Height above primary injection zone: 3,366 ft

Penetrations of secondary seal within 12.4-mi radius: 10

 

Number of faults cutting primary seal within 12.4-mi radius: 0

Distance to nearest mapped fault: 9.3 mi

 

Storage Capacity:

Calculated CO2 storage capacity, primary injection zone: Table 3

Data compiled and interpreted from well records maintained by the Kentucky Geological Survey.

 

Table 2 - Model input values used for study.
 Regional StudyLocal Study
Physical ParametersMeanStandard DeviationMeanStandard Deviation
Area (km2)1256.601.540
Gross Thickness (m)258.056258.07
Porosity (%)7.43.57.43.5
Pressure (MPa)22.0222.02
Temperature (oC)48.0548.05
     
Efficiency FactorsP10P90P10P90
Net-to-Total Area0.850.950.950.99
Net-to-Gross Thickness0.850.950.850.95
Effective-to-Total Porosity0.850.950.850.95

 

Table 3 - CO2-SCREEN modeling results using the CO2BRA (2022) efficiency Mmodel.
Lithology and Depositional EnvironmentStorage Efficiency (%)Injection DurationStorage Resource (Mt)
P10P50P90RadiusGridsP10P50P90
Sandstone, Fluvial7.5810.0012.711 year20 km1471.821017.942132.58
Sandstone, Fluvial9.5613.0716.545 years20 km1641.021308.092597.22
Sandstone, Fluvial10.3514.0317.5810 years20 km1724.611486.993010.04
Sandstone, Fluvial10.7614.7818.5520 years20 km1859.941836.433684.55
Sandstone, Fluvial10.3214.6018.6430 years20 km1982.872012.24070.83
Sandstone, Shallow Marine4.236.409.261 year0.7 km10.671.382.76
Sandstone, Shallow Marine6.128.2210.665 years0.7 km10.921.773.31
Sandstone, Shallow Marine6.649.3512.6510 years0.7 km11.021.983.86
Sandstone, Shallow Marine8.2911.4515.1220 years0.7 km11.262.454.64
Sandstone, Shallow Marine9.0812.7416.9030 years0.7 km11.412.725.14
Sandstone, Fluvial7.009.3211.891 year20 km1799.141578.163150.48
Sandstone, Fluvial8.9412.1215.425 years20 km11018.872040.774044.58
Sandstone, Fluvial9.5413.0316.4510 years20 km11082.842230.044389.95
Sandstone, Fluvial9.9313.6317.2520 years20 km11124.032304.834667.09
Sandstone, Fluvial9.6413.0717.4030 years20 km11133.192314.234615.67
Sandstone, Shallow Marine3.956.008.711 year0.7 km11.132.154.09
Sandstone, Shallow Marine5.697.6410.005 years0.7 km11.442.775.35
Sandstone, Shallow Marine6.168.7311.7810 years0.7 km11.552.975.54
Sandstone, Shallow Marine7.7110.6514.1220 years0.7 km11.613.115.86
Sandstone, Shallow Marine8.4211.7815.7130 years0.7 km11.593.075.91