DOE OSTI2014
As a part of a larger project co-funded by the United States Department of Energy (US DOE) to evaluate the potential of formations within the Cambro-Ordovician strata above the Mt. Simon as potential targets for carbon sequestration in the Illinois and Michigan Basins, the Illinois Clean Coal Institute (ICCI) requested Schlumberger to evaluate the potential injectivity and carbon dioxide (CO2) plume size of the Cambrian Potosi Formation. The evaluation of this formation was accomplished using wireline data, core data, pressure data, and seismic data from the US DOE-funded Illinois Basin–Decatur Project (IBDP) being conducted by the Midwest Geological Sequestration Consortium in Macon County, Illinois. In 2010, technical performance evaluations on the Cambrian Potosi Formation were performed through reservoir modeling. The data included formation tops from mud logs, well logs from the VW1 and the CCS1 wells, structural and stratigraphic formation from three dimensional (3D) seismic data, and field data from several waste water injection wells for Potosi Formation. Intention was for two million tons per annum (MTPA) of CO2 to be injected for 20 years. In the preceding, the 2010 Potosi heterogeneous model (referred to as the "Potosi Dynamic Model 2010" in this topical report) was re-run using a new injection scenario; 3.2 MTPA for 30 years. The extent of the Potosi Dynamic Model 2010, however, appeared too small for the new injection target. It was not sufficiently large enough to accommodate the evolution of the plume. The new model, Potosi Dynamic Model 2013a, was built by extending the Potosi Dynamic Model 2010 grid to 30 miles x 30 miles (48.3km x48.3km), while preserving all property modeling workflows and layering. This model was retained as the base case of Potosi Dynamic Model 2013a. The Potosi reservoir model was updated to take into account the new data from the verification well VW2 which was drilled in 2012. The new porosity and permeability modeling was performed to take into account the log data from the new well. Revisions of the 2010 modeling assumptions were also done on relative permeability, capillary pressures, formation water salinity, and the maximum allowable well bottomhole pressure. Dynamic simulations were run using the injection target of 3.2 MTPA for 30 years. This new dynamic model was named Potosi Dynamic Model 2013b. Due to the major uncertainties on the vugs permeability, two models were built; the Pessimistic and Optimistic Cases. The Optimistic Case assumes vugs permeability of 9,000 mD, which is analog to the vugs permeability identified in the pressure fall off test of a waste water injector in the Tuscola site, approx. 40 miles (64.4km) away from the IBDP area. The Pessimistic Case assumes that the vugs permeability is equal to the log data, which does not take into account the permeability from secondary porosity. The probability of such case is deemed low and could be treated as the worst case scenario, since the contribution of secondary porosity to the permeability is neglected and the loss circulation events might correspond to a much higher permeability. It is considered important, however, to identify the range of possible reservoir performance since there are no rigorous data available for the vugs permeability. The Optimistic Case gives an average CO2 injection rate of 0.8 MTPA and cumulative injection of 26 MT in 30 years, which corresponds to 27% of the injection target. The injection rate is approx. 3.2 MTPA in the first year as the well is injecting into the surrounding vugs, and declines rapidly to 0.8 MTPA in year 4 once the surrounding vugs are full and the CO2 start to reach the matrix. This implies that according to this preliminary model, a minimum of four (4) wells could be required to achieve the injection target. This result is lower than the injectivity estimated in the Potosi Dynamic Model 2013a (43 MT in 30 years), since the permeability model applied in the Potosi Dynamic Model 2013b is more conservative. This revision was deemed necessary to treat the uncertainty in a more appropriate manner. As the CO2 follows the paths where vugs interconnection exists, a reasonably large and irregular plume extent was created. For the Optimistic Case, the plume extends 17 miles (27.4km) in E-W and 14 miles (22.5km) in N-S directions after 30 years. After injection is completed, the plume continues to migrate laterally, mainly driven by the remaining pressure gradient. After 100 years post injection, the plume extends 20 miles (32.2km) in E-W and 15.5 miles (24.9km) in N-S directions. Should the targeted cumulative injection of 96 MT be achieved; a much larger plume extent could be expected. For the Optimistic Case, the increase of reservoir pressure at the end of injection is approximately 1200 psia (8,274 kPa) around the injector and gradually decreases away from the well. The reservoir pressure increase is less than 30 psia (206.8 kPa) beyond 14 miles (22.5km) away from injector. Should the targeted cumulative injection of 96 MT be achieved; a much larger areal pressure increase could be expected. The initial reservoir pressure is nearly restored after approximately 100 years post injection. The presence of matrix slows down the pressure dissipations. The Pessimistic Case gives an average CO2 injection rate of 0.2 MTPA and cumulative injection of 7 MT in 30 years, which corresponds to 7% of the injection target. This implies that in the worst case scenario, a minimum of sixteen (16) wells could be required to achieve the injection target. The present evaluation is mainly associated with uncertainty on the vugs permeability, distribution, and interconnectivity. The different results indicated by the Optimistic and Pessimistic Cases signify the importance of vugs permeability characterization. Therefore, injection test and pressure interference test among the wells could be considered to evaluate the local vugs permeability, extent, and interconnectivity. Porosity mapping derived from the seismic inversion could also be used in the succeeding task to characterize the lateral porosity distribution within the reservoir. With or without seismic inversion porosity mapping, it is worth exploring whether increased lateral heterogeneity plays a significant role in Potosi injectivity. Investigations on vugular, dolomitic outcrops suggest that there may be significantly greater lateral heterogeneity than what has been modeled here. Facies modeling within the Potosi has yet to be thoroughly addressed. The carbonates during the time of deposition are believed to be regionally extensive. However, it may be worth delineating the reservoir with other regional wells or modern day analogues to understand the extent of the Potosi. More specifically, the model could incorporate lateral changes or trends if deemed necessary to represent facies transition. Data acquisitions to characterize the fracture pressure gradient, the formation water properties, the relative permeability, and the capillary pressure could also be considered in order to allow a more rigorous evaluation of the Potosi storage performance. A simulation using several injectors could also be considered to determine the required number of wells to achieve the injection target while taking into account the pressure interference.