Well Intergity
Well Intergity
Well Intergity
[email protected] 614-424-7392
Recipient: Battelle
505 King Avenue
Columbus, OH 43201-2696
DUNS Number: 00 790 1598
Project Grant Period: October 1, 2015, and ending September 30, 2018
U.S. Department of Energy Disclaimer
Battelle Disclaimer
Support for this project was provided by the United States Department of
Energy National Energy Technology Laboratory (NETL) under award
#DE-FE0026585. The project team very much thanks NETL for supporting
this effort, which is an important issue for CO2 storage projects. Project
guidance was provided by Bill O’Dowd (NETL). The project team thanks
Core Energy and the Petroleum Technology Research Council for
providing access to CO2 wells for testing and other support information.
Support for the Appalachian Basin site characterization was provided by
the West Virginia Geologic and Economic Survey, with special thanks to
Phil Dinterman, Gary Daft, Jessica Moore, and Ron McDowell.
Project results reflect contributions from many people on the project team.
The project lead was Battelle, and J.R. Sminchak was the project manager.
Mark Moody and Dr. Neeraj Gupta were technical advisors. Task leads
were Dr. Andrew Duguid, Matt Place, and Jared Hawkins. Other
contributions were made by Autumn Haagsma, Desiree Padgett, Kelly
Quigley, Leo Chiques, Jen O’Brian, Jacob Markiewicz, and others. Carol
Brantley, Colin Kegel, Theresa Schaff, and Jackie Gerst assisted with
project administration and management.
Abstract/Plain Language Summary
This Final Technical Report presents the findings for the project Integrated Wellbore Integrity Analysis
Program for CO2 Storage Applications (FE0026585). The project is part of a U.S. Department of Energy
(DOE) National Energy Technology Laboratory (NETL) program to develop and advance technologies
that will significantly improve the effectiveness and reduce the cost of implementing carbon storage. The
project is designed to assess well integrity for wells exposed to carbon dioxide (CO2) environments in the
subsurface through a combination of field testing and record analysis.
1.2 Objectives
The project goal was to develop an integrated program to identify, survey, measure, and analyze CO2
migration in wellbores. After a well has been constructed and/or plugged, the only indication of migration
through the outer well materials may be pressure buildup on the well, referred to as sustained casing
pressure (SCP). The impact of CO2 on wellbore integrity was determined by integrating field casing
pressure test results with analysis of cement sealing potential, well construction details, well logs, cement
bond logs (CBLs), and well history. In addition, the types of well defects (micro-annulus, cracks, porous
cement, and incomplete cement coverage) were explored by analyzing casing pressure buildup curves
measured in the field on CO2 wells. Meta-modeling methods were used on CO2 storage test fields to
investigate the impact of pressure, gas saturation, and chemistry on well integrity.
The project was organized into eight main tasks with more detailed subtasks. Table 1-1 lists the task
objectives, technical approach, and corresponding deliverables. A key objective of the project was to
complete field measurements of casing pressure at several locations with existing boreholes that have
wells exposed to CO2 in the Michigan Basin, Appalachian Basin, and Williston Basin. The casing
pressure testing results were integrated with analysis of cement sealing to predict well integrity problems
in wells exposed to CO2 in terms of leak location, nature, and severity. Based on SCP test results at the
field sites, Tasks 6 and 7 were modified to evaluate cement sealing conditions.
1.4 Summary of Appalachian Basin, Michigan Basin, and Williston Basin Field
Sites
A key part of the project was testing and analysis at three field test sites. The three field test sites
identified for wellhead SCP testing were located in the Appalachian, Michigan, and Williston
sedimentary basins. These sites have wells exposed to CO2 at depths of 1,000 to 7,000 ft and 5 to 50+
years of age (Table 1-3). Therefore, they provided an excellent opportunity to examine CO2 storage
wellbore integrity at field sites. The wells were surveyed for indications of SCP, and a subsample of wells
were tested for wellhead casing pressure buildup. The field test sites were characterized for geologic
setting, field history, well construction specifications, and hydrologic conditions. Field site characteristics
are summarized in Table 1-3 and are discussed in the following subsections.
The Appalachian Basin Indian Creek site is a natural CO2 and methane field located in Kanawha County,
West Virginia (Figure 1-2). The field contains approximately 58 wells at total depths between 6,200 ft
and 6,700 ft. The Indian Creek field produces in the Tuscarora sandstone, where the percent of CO2 in
some wells ranged from 44% to 83% and nitrogen ranged from 13.9% to 35% (Avary, 1996). A
completion report for the discovery well, API number 4703901684, listed CO2 at 65%. Two scout cards,
for wells 4703902718 and 4703902719, both listed CO2 at 60%. Hamak & Sigler (1991) and Hamak &
Gage (1992) reported produced gas with an average CO2 at 65.8%, 305 British thermal units per cubic
foot, and 1.214 grams per centimeter average gas gravity. One gas sample in the field taken by Jenden et
al. (1993) reported a CO2 content of over 61%. At depth, pressure conditions in the Indian Creek field
were likely >2,800 pounds per square inch (psi) and temperatures were >110 °F. Thus, the CO2 would be
in supercritical state.
Figure 1-2. Location of the Indian Creek field in Kanawha County, West Virginia.
The Michigan Basin study site is located in the northern portion of the Niagaran reef trend in Otsego
County, Michigan (Figure 1-3). The fields have been developed since the 1960s in the region, and
selected reefs have been subject to CO2 enhanced oil recovery (EOR) since the 1990s. There are several
hundred Antrim gas wells at depths between 1,000 ft and 1,500 ft, and there are dozens of Niagaran reef
wells at depths ranging from 5,000 ft to 7,000 ft in Otsego County. The Antrim shale wells produce
methane and CO2. The CO2 volume in the produced gas is between 5% and 30%. Some of the Niagaran
reefs in the area have been subject to CO2 EOR over the past 10 to 20 years. Overall, the Antrim shale
wells and Niagaran reef CO2 EOR wells are attractive candidates for well integrity analysis.
Williston Basin Site
The Williston Basin is a large sedimentary basin which spans the southern portion of Saskatchewan,
Canada, and the north-central United States. The Weyburn field is located on the northwestern edge of the
Williston Basin geologic feature and was a major oil play which is now used for CO2 EOR (Figure 1-4).
The main reservoir for the Weyburn oil field consists of the Marly Midale and Vuggy Midale beds, which
are part of the Mississippian-aged Madison group (Wilson & Monea, 2004). Approximately 3,000 wells
are located in the Williston Basin testing site. The testing site is a mature oil field that began production in
1954. Prior to the start of CO2 EOR, the Williston Basin testing site produced roughly 340 million barrels
of oil, or an estimated 25% of the field’s total reserve. As such, conventional methods for oil production
are no longer viable for profit. To keep the oil field active, CO2 EOR was started at the Weyburn field in
2000 using CO2 transported via pipeline from a gasification synfuels plant in North Dakota.
Well integrity issues may arise from the materials, construction, operations, and subsurface conditions
specific to a well. The objective of this task was to develop a wellbore integrity registry that describes
actual or potential well integrity issues that may affect carbon storage projects. The registry was based on
existing research and on experience related to well construction methods, well casing integrity issues,
well cement issues, and geologic processes for CO2 environments.
The carbonate species from Equation 1 interact with the cement, C-S-H, and CH to create CaCO3
(Equations 2 through 5):
Ca(OH)2(s) + 2H+ + CO32- → CaCO3(s) + 2H2O [Equation 2]
Ca3Si2O7H•4H2O(s) + 2H+ + CO32- → CaCO3(s) + SiOxOHx(s) [Equation 3]
Ca(OH)2(s) + H+ + HCO3- → CaCO3(s) + 2H2O [Equation 4]
Ca3Si2O7H•4H2O(s) + H + HCO3 → CaCO3(s) + SiOxOHx(s)
+ -
[Equation 5]
For cement exposed to wet supercritical CO2 and carbonated brine that is not refreshed often (diffusion
conditions with a small volume of brine as compared to the volume of cement), the reaction with CO2
stops with CaCO3 being created within the cement pores. However, if cement is exposed to flowing
carbonated brine or the carbonated brine is refreshed often, the carbonate that formed in the pores can be
reacted away (Equations 6 and 7), leaving a soft silicate-hydrate gel (SiOxOHx in Equations 3 and 5).
CO2 + H2O + CaCO3(s) ↔ Ca2+ + 2HCO3- [Equation 6]
2H+ + CaCO3(s) ↔ CO2 + Ca2+ + H2O [Equation 7]
Figure 2-2. Photos showing the results of well cement exposed to flowing carbonated brine
at 50°C and pH 3. Left photo shows the outer orange and brown silicate-hydrate zone
outside of a CH-depleted cement zone, with an unreacted cement in the center.
Right photo shows the outer orange and brown silicate-hydrate zones, the
white CaCO3-rich zone, and a CH-depleted cement zone in the center.
Multiple field investigations have been conducted on wells exposed to CO2 or generally on well
construction as it relates to carbon capture and storage (CCS) leakage risk. Carey et al. (2007), Crow et al.
(2010), and Duguid et al. (2014) have all studied wells to identify well integrity defects and identify
changes brought on by exposure to CO2. Carey et al., Crow et al., and Duguid et al. all found evidence of
carbonation in cements in and above the CO2 reservoir.
• Carey et al. (2007) looked at a 50-year-old well that was exposed to CO2 for 30 years in the
SACROC oil field in Scurry County, Texas. The cement was collected at the surface during a
sidetracking operation. The cement in the well was a neat Portland Type 1. The authors saw
carbonation and discoloration of the samples and measured a permeability of an air-dried sample
of 0.1 millidarcy.
Task 3 was divided into three subtasks to compile available information on existing wells at the field test
sites in Otsego County, Michigan (Michigan Basin site), Kanawha County, West Virginia (Appalachian
Basin site), and Saskatchewan, Canada (Williston Basin site). The well data collection task was
completed in June 2017. Information was compiled and evaluated for well cementing/drilling, operation,
and well workover records. This dataset included over 1,000 items related to wellbore construction in the
three study sites. The information from the well collection task will be used to evaluate the overall
condition of boreholes in the study areas and as input for future tasks (Figure 3-1).
Oil and gas records were acquired from publicly available resources such as the Michigan Department of
Environmental Quality (DEQ), the West Virginia Geologic Survey (WVGES), and the Government of
Saskatchewan. The primary focus was collecting data on active and producing wells. Operational and
workover records were also acquired from the same sources. Additional, proprietary information was
shared by select operators.
Figure 3-2. Typical Antrim (left) and Niagaran (right) well construction.
CBLs were available for 87 wells in the Michigan Basin field site area. The types of cement were not
typically entered in well records. The surface casing in Antrim wells was mostly driven (95%) with only
5% having recorded using cement. The amount of cement used ranged from 350 to 460 sacks. The cement
used in the intermediate casing ranged from 160 to 700 sacks of cement, with an average of 360 sacks.
Some wells did not have cement data available but were recorded to have had cement circulate to the
surface. The deep casing had recorded cement ranging from 55 to 660 sacks with an average of 243 sacks.
Occasionally, wells were recorded to have had cement circulate to the surface.
The conductor casing in Niagaran wells was mostly driven (70%) with 18% having recorded using
cement. The amount of cement ranged from 100 to 200 sacks. The remaining 12% of wells did not have
recorded cement. The cement used in the second surface casing string ranged from 200 to 1,130 sacks,
with an average of around 600 sacks. A subset of wells (9%) did not have a second casing string. The
Figure 3-3. Monthly and cumulative oil production from Niagaran wells showing
greatest monthly production in the early 1980s and a steady monthly decline through 2016.
In addition to construction data, information on perforations, treatments, and treatment volumes were
recorded. These data provide information on well stimulations and additional stresses applied to the
wellbores. Only 22% of the wells were recorded to have perforations, while many wells were completed
open hole. Nearly 60% of the wells were treated with acid and fractured. The volumes and types of
materials used for treatments varied greatly by well.
CBLs were available for 15 wells in the Indian Creek field. The CBL analyses were conducted, and are
summarized, under Task 4. Results showed moderate thickness of good quality cement with frequent
intervals of poor cement (Figure 3-5). Additionally, caliper logs were collected to be used in Task 4 with
CBLs and cement sacks to determine the volume and thickness of a cement column. The number of sacks
of cement was recorded for each casing string; cement type was not readily available. Eight wells did not
have cement data available. The conductor casing was set with 10 to 826 sacks of cement with an average
of 183 sacks of cement. Wells were frequently recorded to have had cement circulate to the surface. The
surface casing string used 150 to 1,121 sacks of cement with an average of 548 sacks. It was also
commonly recorded to have been cemented to the surface. The intermediate casing used 148 to 1,226
sacks of cement with an average of 544 sacks. The fourth casing and deep casing strings used similar
amounts of cement, ranging from 75 to 280 sacks.
Production data were available for 52 wells in the Indian Creek field on a monthly basis since 1985. Over
40 billion standard cubic feet (Bcf) of gas were produced since 1985 between the 52 production wells.
The cumulative production by well ranged from 0 MCF to over 4.3 million MCF. On average, a well had
a cumulative gas production of 780,000 MCF. The Indian Creek field experienced two major periods
without production: from January 1987 to October 1987 and again from January 2012 to October 2013.
Production ended in January 2015, and the field has recently been sold to new ownership.
Figure 3-6 shows the annual and cumulative production of gas for each well in the Indian Creek field; in
the figure, the periods of no production are visible. The production of gas was not consistent across the
field. The greatest amount of gas was produced in the southernmost wells, while the least amount was
produced in the northern wells.
No well workover or leakages were recorded for wells in the Indian Creek field. A search on the West
Virginia DEP oil and gas database showed no violations listed for Indian Creek wells. Detailed operator
records were unavailable for the field. Informal discussions with the field technician suggest that the wells
had few problems with no significant cause for frequent repairs or other corrective actions.
Casing details were available for 1,328 wells at the Williston Basin site. The wellbores were consistently
constructed using two casing strings: a surface casing and a production casing. The diameter of the
surface casing ranged from 7.0 to 13.375 inches. The most common surface casing diameter was 8.625 or
9.625 inches, with 96% of the wells drilled no deeper than 1,000 feet. The production casing ranged in
diameter from 2.875 to 10.75 inches. Approximately 50% of the wells had 7-inch-diameter production
casing, with 30% of wells having production casing of 5.5 inches in diameter and 18% with 4.5-inch
casing. The production casing was set at depths up to 5,600 ft, with most wells set between 4,000 and
5,000 ft deep to target the Midale.
Cementing data were recorded as cement volumes rather than the amount of materials used. Cement
volumes were only available for 135 wells. The volumes ranged from 170 cubic feet (ft3) to 1,560 ft3,
with an average of 545 ft3. Most wells had less than 400 ft3 of cement. Figure 3-8 shows a histogram of
cement volumes for available wells at the Williston Basin site.
Full operational data and workover histories were not publicly available for Williston Basin wells.
However, numerous wells (80%) were tested for SCP as part of a monitoring program. The level of SCP
was categorized by leakage severity. Most of the wells tested did not show any SCP or leakage through
the wellbore; a small percentage had minor SCP associated with non-serious leakage. Some wells (34%)
showed significant SCP and were categorized as having serious leakage. Figure 3-9 illustrates the number
of wells that fell into each leakage category.
Figure 3-9. Casing pressure category for wells at the Williston Basin site.
The objective of Task 4 was to analyze well integrity for the three field sites based on logs and well
records. Task 4 was divided into three subtasks: Log Analysis, Well Record Analysis, and Well Integrity
Evaluation. The well risk analysis was focused on field test sites in Otsego County, Michigan (the
Michigan Basin site), Kanawha County, West Virginia (the Appalachian Basin site), and Saskatchewan,
Canada (the Williston Basin site). Task 4 was completed in October 2017. Well log data for the Michigan
Basin and Appalachian Basin field sites were collected in Task 3. Log data for Weyburn Field were
collected as part of the SaskCO2USER project (Duguid et al., 2011) and were provided to Battelle already
interpreted for risk assessment. The available log data are summarized in Table 4-1.
Table 4-1. Number of logs used for Otsego, Michigan, and Indian Creek, West Virginia, fields.
Log Type Otsego, MI Indian Creek, WV
Cement Bond Log (CBL) 1 12
Gamma Ray (GR) 54 11
Formation Bond (FBL) -- 7
Neutron Log (NEU) 43 1
Density (DEN) 38 3
Acoustic (XMAC) 26 --
Resistivity (RT) 49 --
Pulsed Neutron Capture (PNC) 3 --
Caliper Log (CAL) 53 36
Bit Size Log (BIT) 6 --
Table 4-2. Number of wells in each field with known data for
each category ranked in the well integrity evaluation.
Risk Michigan Appalachian Williston
Field
Type Basin Basin Basin
Total number of wells (records collected) 54 47 1,391
Cemented through caprock 51 41 176
Likelihood of
leakage
Figure 4-1. Histograms of Total Likelihood and Total Severity Scores (top) and
Total Risk Scores (bottom), Northern reef trend, Otsego County, Michigan.
Total Risk Scores for the Appalachian Basin wells ranged from 55 to 165 (the minimum possible Total
Risk Score was 15, while the maximum possible Total Risk Score was 375). The average Total Risk
Score was 98.2 with a standard deviation of 45.3. The histogram for risk scores was a bell-shaped curve
centered between the bins 76-100 and 101-125 that skewed toward lower values (Figure 4-3). More than
half of the wells had a Total Risk Score of 100 or less (25 locations). The wells with lower Total Risk
Scores are located in the eastern and central portions of the field, away from population centers
(Figure 4-4), suggesting that distance to population was a determining factor for total risk.
Histograms for the Total Likelihood Scores, Total Severity Scores, and Total Risk Scores for Williston
Basin wells included in this study are presented in Figure 4-5. Maps are not presented for the Weyburn
Field at the request of the operator that provided the data. Total Risk Scores for the Williston Basin wells
included in this study range from 27 to 300 (the minimum possible Risk Score is 15, while the maximum
possible Risk Score is 375). The average Total Risk Score was 101.7 with a standard deviation of 25.6.
The histogram of Total Risk Scores for the wells included in this study is skewed toward lower values.
The largest number of wells has a Total Risk Score between 101 and 125. More than half of the wells
included in this study (756 locations) had Total Risk Scores of 100 or less.
Figure 4-4. Total Risk Scores, Indian Creek Field, Kanawha County, West Virginia.
The three field test sites were assessed using a Features, Events, and Processes (FEP)-based risk ranking
methodology. The fields varied in size, geography, and geology. The assessment showed that Total
Severity was probably more important in determining Total Risk. Correlation analyses of the Total
Likelihood data showed that cement through caprock, well age, and well status had the largest effect on
the Total Severity rank, with cement through caprock being the most important for each field. The CBL
analysis conducted for the Indian Creek and Otsego wells shows that some wells may have primary
cement that may pose a risk with values of less than 80% bond index.
The effect of unknown data was also most pronounced in the Alberta dataset, with only 176 of 1,391
wells with data for the cement through caprock category. However, for each of the other fields, the most
missing data were in the cement through caprock category as well. The strong correlation may be due to
giving the worst-case score to categories with missing data, but it also highlights the need to collect
additional data to ensure that these categories can be fairly and fully evaluated. Correlation analyses of
the data show that Total Risk follows more closely the Severity categories than the Likelihood categories.
This finding has implications for risk mitigation. Mitigation of severity is most likely to come from early
warnings from monitoring tools, implying that the receptors of concern (the local population,
groundwater, and surface water and environmentally sensitive areas) may need to be monitored.
4.4 Key Findings of Log- and Testing-Based Well Integrity Assessment
The three field test sites in Otsego, Michigan, Indian Creek, West Virginia, and Alberta, Saskatchewan,
were assessed using a risk ranking methodology. The fields varied in size, geography, and geology. The
assessment showed that Total Severity was more important in determining the Total Risk.
Correlation analyses of the Total Likelihood data showed that cement through caprock, well age, and well
status had the largest effect on the Total Severity rank, with cement through caprock being the most
important for each field. This finding may indicate that well isolation for each CCS project may need to
be better understood in order to properly assess leakage risk. The CBL analysis conducted for the Indian
A key part of the project was testing SCP in CO2 wells at the field test sites. The tests provide a means to
measure the nature and severity of well defects in CO2 wells. The procedures, results, and analysis of SCP
field testing are provided in Sections 5.1 (Michigan Basin) and 5.2 (Williston Basin).
5.1 Michigan Basin SCP Field Testing
Niagaran reefs at the Michigan Basin field test site have been subject to CO2 EOR since 1996, with
operations expanding to a total of 10 reefs. The Niagaran fields were developed since the 1960s in the
region. The production history of the tested reefs is outlined in Figure 5-1. Overall, the Antrim
shale/Niagaran reef is an attractive area for examining wells exposed to different types of CO2
environments in the subsurface. The site also has other hydrocarbon wells and injection wells for
comparison testing (Figure 5-2).
Multiple wells in the Michigan Basin were measured for casing pressure. If well conditions indicated that
pressure might continue in the well, the well was selected to be part of enhanced testing. Six selected
wells were tested using a pressure bleed-down/buildup test, and gas samples were collected from the
wells. In addition, the gas chamber volume was measured in select wells. Table 5-1 summarizes the well
construction specifications for the 23 wells surveyed for indications of SCP. The majority of the tested
wells were drilled as primary production wells in the 1970s; however, a few of the wells were drilled as
part of recent CO2 EOR operations.
Monitoring was performed by recording the temperature and pressure to have a record of pre-test
conditions. The process was continued through depleting the pressure in the annulus to atmospheric
conditions. In cases where gas chamber volume was needed, a flow meter was attached to the correct
valve, and the volume of flow out of the annulus space was recorded. Later, these data were converted
from atmospheric conditions (pressure and volume coming out of the wellhead into the atmosphere) using
the annulus pressure and a form of the Ideal Gas Law formula to calculate the gas volume released.
After opening the valve, venting the gas, and properly closing the valve, a memory gauge was attached.
This was done to monitor the pressure buildup in the annulus. The gauge was programmed to record a
pressure/temperature measurement at a preselected interval of time, such as a minute-by-minute
measurement. The well was then left alone for a minimum of 2.5 weeks to see if the pressure would build
up in the annulus. After this time period had fully passed, the gauge was removed and the pressure
buildup data were collected. The data revealed how long it took for the pressure to return to previous
levels, as well as the influence of thermal changes on the gas in the annulus.
If there was SCP in the annulus of the well, a gas sample was taken. This sample would have been
procured after the initial drawdown of the well to ensure that the sample was representative of the gas that
was making its way into the annulus. It was important to ensure that the sampling cylinder was visually
A pressure bleed-down/build-up test was performed on each well involved in this study to evaluate the
pressure-response curve related to the SCP. The casing valve associated with the annulus being tested was
then opened to allow the pressure to bleed down to near-atmospheric conditions. Next, a data-recording
pressure/temperature gauge was connected to the annular space to monitor the pressure recovery curve.
Figure 5-4 shows the comparison of pressure and temperature versus time for the Michigan Basin wells.
For MB-1, MB-4, and MB-6, it is apparent that temperature is driving the changes in pressure. MB-3
seems to be building pressure over time. MB-2 and MB-5 are inconclusive; these tests, for example, may
have been affected by the formation of ice, hydrate, or another substance. The beginning and end points
show where the gauge is reading standard conditions before and after being connected to the well.
In two of the wells, the annular space was bled down through a flow meter to determine the volume of the
gas released; a valve at the outlet of the wellhead maintained a constant flow rate. Using the time required
to blow the well down to near-atmospheric levels and the initial pressure, the gas volume released was
calculated. The total volume and the initial pressure were used with the Ideal Gas Law to determine the
gas volume released. The volume of gas released from the annulus of well MB-1 was calculated to be
about 18 ft3, which was calculated using a well pressure of 17 psia, and bleeding down to a pressure of
15 psia. The total open annulus volume was calculated to be about 420 ft3, using standard casing and hole
sizes. The volume of gas released from the annulus of well MB-3 was calculated to be about 66 ft3, which
was calculated using a well pressure of 32 psia, and bleeding down to a pressure of 18 psia. The total
open annulus volume was calculated to be about 865 ft3, using standard casing sizes. The open annulus
Figure 5-5. Map of studied wells in the Williston Basin field test site.
Table 5-4 summarizes the well construction specifications for the 17 Williston Basin wells tested for
SCP. Figure 5-6 shows the wellbore diagrams for WB-1, WB-4, WB-5, WB-10, WB-13, and WB-15
(other well diagrams are provided in the full field testing report). Testing showed that these six wells had
casing pressure buildup. The tested wells were between 1,381 and 1,582 meters deep (measured depth
[MD]). The tested wells in the field are being used for production, WAG, or water injection. Eight of the
wells were drilled in the 1950s, with the rest being drilled in the 1990s or 2000s. CO2 EOR operations
began in 2000, and the wells have been in CO2 environments in the subsurface since the early 2000s,
some of which have likely been exposed to CO2 since being drilled. Detailed operational data, well
histories, and certain details pertaining to well construction, such as cementing records, were not available
from the operator.
After this process, the gauge was removed, and the pressure buildup data were collected. The data
revealed how long it took for the pressure to return to previous levels, as well as how much influence
thermal changes throughout the course of the day had upon the gas in the annulus. Ideally, if there was
SCP in the annulus of the well, then a gas sample would have been collected after the initial drawdown of
the well to ensure that the sample was representative of the gas making its way into the annulus. Each
sample would have been tracked with a tag that included the well name, API, date/time, and all other
pertinent identifiers. These samples would then be sent to a laboratory for composition testing. However,
for this site, gas samples could not be collected due to low pressure and flow at the wellhead.
Consequently, historical operator records were examined to determine the gas chemistry and the source of
the gas. These operator records were compared to available simple gas samples taken by the operator,
which were used for a rough indication of where the gas was coming from. Low amounts of CO2 indicate
that CO2 is not originating from the injection zone.
Seventeen (17) of the 60 wells considered were tested with the apparatus containing the pressure and
temperature gauge. The 17 wells were ultimately chosen by the operator because of the operator’s
extensive knowledge of the field and well characteristics. Battelle was informed by the operator that a gas
A pressure bleed-down/buildup test was performed on each of the 17 wells involved in this study to
evaluate the pressure-response curve related to the SCP. The casing valve associated with the annulus
being tested was then opened to allow the pressure to bleed down to near-atmospheric conditions. Next, a
data-recording pressure/temperature gauge was connected to the annular space to monitor the pressure
recovery curve.
Figures 5-8 and 5-9 show the comparison of pressure and temperature versus time. For WB-1, WB-4,
WB-5, WB-10, WB-13, and WB-15, buildup can be seen. WB-1 shows a moderate dip in the pressure
during the middle of the test period, which may have been a device error. For WB-1, the test
administrator had noted that the testing apparatus ‘froze up’ during the test. WB-2 and WB-3 show no
significant pressure increase, with the temperature in WB-3 having a slight effect on the pressure. WB-5
shows a possible mechanical defect, indicated by the sharp logarithmic increase in pressure. WB-10
shows an unexplained dip in the pressure buildup, perhaps due to a small release of pressure, or possibly
another factor. The other wells without buildup have varying amounts of pressure dependence on
temperature, with a handful of the wells experiencing a negligible decline in pressure.
Eleven (11) of the 17 wells tested at the Williston Basin site showed no significant pressure rebound
pattern. However, six of the wells did exhibit a minor pressure buildup pattern that might be useful for the
operator to evaluate potential well defects present in shallow zones. Well construction information,
subsurface conditions, and pressure buildup data were input into the MATLAB script designed to
estimate well defect factors based on the methodology by Moody and Dotson (2015). This methodology
examines the gas influx rate change over time to determine the nature of a well defect based on a defect
model curve, similar to pressure transient well testing methods that examine pressure derivative curves.
There was some uncertainty on several input parameters related to height of liquid volume in the annulus,
gas volume due to low flow conditions, and gas properties. It was assumed that the calculations used
500 ft true vertical depth (TVD) of liquid length and a gas chamber length of 30 ft.
The SCP analysis for the six wells suggests mainly a pressure-limited orifice type defect factor. None of
the wells correlated to a permeable cement flow model, which would suggest cement degradation and gas
migration from the reservoir zone. Figure 5-11 shows the results for the wellbore integrity analysis. For
each well, the analysis provides a curve of pressure buildup over time and gas influx over time. The gas
influx curve shows a wellhead model curve of the type of cement defect as detailed by Moody and Dotson
(2015). For the six wells analyzed, the most suitable gas influx curve was determined to be a pressure-
limited orifice (O) or non-pressure-limited vein (V) model. The wells did not show a strong match to the
gas influx rate type curve for orifice type defect, so the understanding of defect type is uncertain.
Table 5-6. SCP analysis results for Williston Basin field site.
Well Asymptotic Defect Instantaneous Sustained Flow Factor
Gas
Well Depth Pressure Response Leakage Metric Leakage Metric (square
Source
(ft) (psi) Type (MSCFD) (MSCFD) microns)
Unclear
WB-1 4,777 Shallow 50 0.027 0.0076 564
(Orifice)
Unclear
WB-4 4,600 Shallow 28 0.026 0.0077 627
(Orifice)
Unclear
WB-5 4,626 Shallow 65 0.033 0.1500 >1,000
(Orifice)
Unclear
WB-10 4,826 Shallow 75 0.036 0.0082 555
(Orifice)
Unclear
WB-13 4,784 Shallow 45 0.022 0.0050 379
(Orifice)
Unclear
WB-15 4,784 Shallow 27 0.220 0.0200 >1,000
(Orifice)
The objective of Task 6 was to analyze field data on mineralogy, fluids, cement, hydrologic conditions,
and CO2 exposure for the three field test sites to determine cement sealing and well integrity relationship.
Data from previous tasks were analyzed for indicators of wellbore integrity factors that may predict
cement sealing conditions in relation to CO2 leakage potential in legacy oil and gas wells. These results
provide a better understanding of wellbore integrity effects for CO2 storage applications.
Appalachian Basin. Mineralogy of reservoir and caprocks for the Appalachian Basin Indian Creek
field was characterized by the WVGES based on analysis of available core samples from the Tuscarora
sandstone. Cores were inspected, and 45 thin sections were prepared for core samples from Clay
County well 4701500513 and Kanawha County well 4703903914. Rock core from the well was in
poor condition, and exact sample depths were not available. The Tuscarora sandstone samples were
>90% quartz with a small fraction of accessory minerals. Rose Hill shale caprock samples were not
readily available, so mineralogy was based on research by Jin et al. (2010). The Rose Hill is mostly
Table 6-2. Summary of brine geochemistry for reservoir samples from the field test sites.
Appalachian Basina Michigan Basin Siteb Williston Basinc
Constituent
(mg/L) Baseline Post-CO2 Injection
Low High Low High
Low High Low High
Ca 10,000 40,000 67,500 110,000 84,900 99,400 1,100 2,000
Mg 800 8,000 7,985 12,100 8,060 11,200 320 490
Na 40,000 110,000 15,300 22,500 14,400 21,300 12,000 26,000
K 400 1,700 11,000 16,600 16,200 18,400 120 600
SO4 300 1,300 81 97.1 44 150 3,180 3,620
Cl 73,000 183,000 251,000 274,000 188,000 270,000 18,000 43,000
HCO3- 10 160 361 468 NDd 956 NA NA
Br 879 4,650 2,280 3,030 NDe 3,250 NA NA
Al 0.54 3.51 ND 1.0 NDf 1.0 NA NA
Fe 186 --- 10 129 52.4 654 NA NA
SiO2 0.005 0.006 6.90x10-4 1.46x10-2 2.04x10-3 2.35x10-2
pH 5.1 6.4 4.83 5.88 4.10 4.87 6.3 --
Alkalinity NA NA
7 104 296 384 NDg 785
(as HCO3-)
TOC NA NA NDf 79 27 343 NA NA
DOC NA NA NDf 66 16 295 NA NA
Salinity 120,000 300,000 348,000 450,000 380,000 450,000 35,000 110,000
Water Sat. 40% 50% 30% 40%
a: Based on produced water for Medina-Tuscarora sandstones in WV & southeast OH from Battelle (2015) and Breen et al. (1985).
b: Baseline (pre-injection) data obtained from samples from five wells in fields (or lobes) that have not seen CO2.
c: Based on data from Cantucci et al. (2009); Mills et al. (2011); and Hassani et al. (2014).
d: Minimum detected concentration: 316 mg/L.
e: Minimum detected concentration: 2,720 mg/L.
f: One detected sample.
g: Minimum detected concentration: 259 mg/L.
Appalachian Basin. Well completion records were reviewed for 53 wells in the Appalachian Basin
field site to determine cement composition. The majority of records did not include details on cement
composition because the well completion forms did not require the information. Records with cement
information show generally a neat Portland Class A cement for the production casing. Records indicate
a 50/50 Pozmix with calcium chloride (CaCl2) and latex additive in the shallow and intermediate
casing strings to approximately 2,300 ft and 5,000 ft. Since the same operator installed most of the
The geochemical model is an equilibrium speciation model. Kinetic reactions for cement are not well-
defined, so the model cannot demonstrate the timing of reactions. The LLNL geochemical model is
calculated based on extended Debye-Hückel, which is not typically used for waters with high ion activity
(e.g., brine). Pitzer calculations are used for slightly higher activity brines; however, the activities of the
brines of interest in this study are beyond the conditions that dictate when Pitzer calculations are usually
used. The specific model setup is based on brine chemistry across each field test site combined with solid
mineral complexes from different locations. Site-specific cement mineral phases were not available, so a
cement value from another site was used.
The midpoint brine sample was used in the model to determine the SIs, mineral precipitation, and
matrix/brine/CO2 interactions (Tables 6-6 and 6-7 for the Appalachian and Michigan Basins,
respectively). Batch reaction modeling was done using a range of parameters guided by the minimum and
maximum values for each parameter for the Michigan Basin and Appalachian Basin brines. Furthermore,
each parameter was extended beyond the maximum value to determine the range of possibilities for other
potential brines. The batch reaction parameters are defined in Table 6-8.
Equilibrium Model
The results of the Appalachian Basin equilibrium model are shown in Figure 6-1. The stability of most
cement minerals, as measured by log SIs, increases from the minimum brine value to the midpoint brine
value. The increase in the log SI value from the midpoint to the maximum value is not as large, even
though the change in parameters was the same between the minimum and midpoint values as the midpoint
and maximum values. This indicates that the change begins to level off as the parameters increase in
value. The most stable cement mineral phases are calcite, aragonite, and vaterite, which are at or near
saturation in the midpoint and maximum values and have log SI values less than -3 in the minimum
sample. The cement mineral SIs calculated using the summary data for the Michigan Basin are presented
in Figure 6-2. Calcite, aragonite, and vaterite are all saturated in the minimum, midpoint, and maximum
values. This is likely due to the buffering capacity provided by the defined alkalinity.
10
0
Log Saturation Indices
-10
-20
-30
-40
-50
-60
10
Log Saturation Indices
-10
-20
-30
-40
Table 6-11. Amount of each mineral phase that precipitates from Appalachian Basin summary
solutions, with and without 1.1 mol of CO2. Precipitates are reported in mol/kg water. Columns
labeled “Δ with CO2” show the change in precipitation when the solution is reacted with CO2.
Dolomite-
Solution Anhydrite Diaspore Hematite
Anhydrite Diaspore ordered Dol-ord Hematite
Minimum --- Δ with 2.00E-5 Δ with --- Δ with 3.34E-8 Δ with
Midpoint --- CO2 7.83E-5 CO2 --- CO2 8.51E-4 CO2
Maximum 3.64E-3 1.30E-4 3.91E-4 8.15E-5
Min. w. CO2 --- --- --- -2.00E-5 --- --- 1.19E-8 -2.15E-8
Mid. w. CO2 --- --- --- -7.83E-5 --- --- 8.51E-4 -0-
Max. w. CO2 3.57E-3 -6.5E-5 --- -1.30E-4 --- -3.91E-4 8.14E-5 -1.00E-7
The mineral phases that precipitated in the Michigan Basin summary solutions are listed in Table 6-12.
Three mineral phases precipitate from all three summary solutions: diaspore (α-AlO(OH)), dolomite
(CaMg(CO3)2), and hematite (Fe2O3). When 1.1 mol of CO2 is added to the solutions, the amount of each
of diaspore and dolomite that precipitates from each solution decreases. The amount of hematite that
precipitates did not change when CO2 was added. The amount of diaspore that precipitated after the
addition of CO2 decreased by about a third in the minimum solution and by less than 3% in the midpoint
and maximum samples. Dolomite did not precipitate in the minimum solution after the addition of CO2
and decreased by around an order of magnitude in the midpoint and maximum samples.
Table 6-12. Amount of each mineral phase that precipitates from Michigan Basin summary
solutions, with and without 1.1 mol of CO2. Precipitates are reported in mol/kg water. Columns
labeled “Δ with CO2” show the change in precipitation when the solution is reacted with CO2.
Dolomite-
Solution Diaspore Hematite
Diaspore ordered Dol-ord Hematite
Minimum 2.92E-6 Δ with 1.22E-3 Δ with 6.79E-5 Δ with
Midpoint 1.80E-5 CO2 2.24E-3 CO2 2.33E-3 CO2
Maximum 3.31E-5 3.79E-3 3.90E-3
Min. w. CO2 --- -2.92E-6 --- -1.22E-3 6.77E-5 -2.00E-7
Mid. w. CO2 --- -1.80E-5 --- -2.24E-3 2.33E-3 -0-
Max. w. CO2 --- -3.31E-5 --- -3.79E-3 3.90E-3 -0-
6.3 Review of Well Defects and Pore Network Dimensions for Field Sites
Field data from the Michigan Basin site were evaluated to determine the range of well defects and pore
network dimensions in relation to potential for CO2 migration. Several researchers have analyzed the
relationship of fracture aperture diameter and cement sealing potential (Carroll et al., 2016). Pore throat
radius of rock core samples, fracture width of reservoir/caprock, and cement features were evaluated to
determine potential for gas migration.
Image Log Fracture Characterization. Image logs were evaluated to estimate the dimensions of induced
and natural fractures around the borehole at the Michigan Basin site. These fractures may be pathways for
CO2 migration and carbonation. Electrical and acoustic image logs are geophysical wireline tools for
describing features around the borehole. Fracture width was measured by hand from a physical copy of
the resistivity formation image log (Baker Hughes, Star tool) from the Lawnichak 9-33 well of the Dover
33 reef field. When fractures with well-defined dimensions were identified, the apparent width of the
fractures was measured by hand in centimeters and converted into degrees to calculate the scaled apparent
aperture. With track width spanning the full 360 degrees of the borehole, features may be measured by
hand or electronically and converted into degrees. Degrees may then be converted into apparent aperture
using the following equation:
𝐴⁰(2𝜋𝑟)
𝐴′ = 360
Where A' is the apparent fracture aperture, A⁰ is the measured feature width in degrees, and r is the
borehole radius (obtained from caliper).
CBL Features. Defects in cement were also evaluated as pathways for CO2 migration. In 2009, a CBL
mapping tool (the Schlumberger Isolation Scanner) was run through the intermediate casing section of the
State Charlton 4-30 well of the Charlton 30/31 fields in Otsego County, Michigan. The isolation scanner
utilizes ultrasonic imaging technologies in combination with flexural wave imaging to yield a full
azimuthal profile of the borehole. This profile can be used to identify cement bond integrity, possible
channels or fractures within cement, casing integrity, and degree of zonal isolation between the formation
and a wellbore. This study utilized the Isolation Scanner to measure and calculate the apparent aperture of
several cement channels observed within the State Charlton 4-30 well.
Cement features like channels or fractures were identified within the Charlton 4-30 well by noting sharp
changes in flexural attenuation where a low value of attenuation appeared. The presence of channels or
fractures was verified using the Solid Liquid Gas Map, where interpreted cement features would be
saturated with either gas or water. The apparent width of the features was measured by hand using a ruler
and a physical copy of the log. With track width spanning the full 360 degrees of the borehole, features
Where A' is the apparent width, A⁰ is the measured feature width in degrees, and r is the production
casing radius.
Four channels/fractures were interpreted to be present within the Charlton 4-30 well, located at 2,552 ft
MD, 2,802 ft MD, 3,014 ft MD, and 3,052 ft MD. A cement channel was identified by its low attenuation
and borehole-parallel shape from 2,552 ft to 2,562 ft. This feature has a measured maximum width of
61 degrees and a calculated maximum width of 106 mm. A large irregular feature interpreted to be a
cement channel was identified from 2,802 ft to 2,828 ft based on its irregular shape, low attenuation, and
gas saturation. The feature’s maximum width was measured as 247 degrees, and its maximum width was
calculated to be 425 mm. A small cement channel/feature was identified from 3,014 ft to 3,017 ft based
on its low attenuation, partial gas-water saturation, and distinct ellipsoid shape elongated perpendicular to
the borehole. The feature’s measured maximum width was 154 degrees, and its calculated maximum
aperture was 270 mm. The fourth and final cement fracture/channel was identified from 3,052 ft to
3,080 ft; it was identified by its low flexural attenuation, distinct shape, and gas saturation (Figure 6-3).
The feature’s measured maximum width was 123 degrees, and its calculated maximum aperture was
210 mm. Over the entire intermediate casing interval, the mean aperture size of cement features was
calculated to be 252 mm, with a standard deviation of 6.65 mm and 95% confidence upper and lower
bounds of 464 mm and 41 mm, respectively (Figure 6-4).
Figure 6-3. Cement channel/fracture identified within the Charlton 4-30 well from
3,052 to 3,080 ft, highlighted by the yellow circle. Note the low attenuation (blue color)
and its gas saturation (red) illustrated on the Solid Liquid Gas Map (right-most track).
6.4 Key Findings of Field Analysis of CO2 Cement Sealing and Well Integrity
Appalachian Basin, Michigan Basin, and Williston Basin field test sites were evaluated for cement sealing
conditions based on mineralogy, brine geochemistry, cement chemistry, hydrologic conditions, and CO2
exposure. The analysis was completed because SCP testing and well history review on CO2 wells
suggested that well defects were not a common problem at these field sites. Therefore, the sites were
investigated to determine if cement carbonation would potentially seal well defects, limiting gas
migration and casing pressure buildup in the CO2 wells. The results provide a better understanding of
wellbore integrity effects for CO2 storage applications as listed below:
Analysis of Subsurface Setting for Cement Sealing. Thin sections, well materials, cement properties,
and hydrologic, temperature, pressure, brine chemistry, and saturation conditions were summarized for
the pertinent CO2 reservoir zones and caprocks to provide a basis for further analysis of CO2 sealing
conditions in the Appalachian Basin, Michigan Basin, and Williston Basin field sites. Key findings of the
subsurface parameters are as follows:
• The field sites have wells exposed to CO2 at depths of 1,000 to 7,000 ft, pressures of 2,000 to
3,000 psi, temperatures of 105 °F to 145 °F, and CO2 exposure durations of 5 to 50 years.
• Reservoirs were sandstone and carbonates with shale and evaporite caprocks. The reservoir zones
contain a fairly large portion of quartz in the sandstone and dolomite in the carbonate reservoirs,
with a small amount of trace minerals.
• Brine geochemistry includes high salinity at the Appalachian Basin and Michigan Basin sites.
• Standard Portland Class A cement and carbon steel well casing were used for well materials at all
three field sites.
Together, these findings help define the subsurface system for the field sites in terms of interactions
between cement, brine, minerals, and CO2 in the subsurface.
Geochemical Analysis to Predict Cement Sealing Conditions. Information compiled on the field sites
was utilized to analyze geochemical interactions in relation to cement sealing conditions and well defects
for the three field sites with geochemical model PHREEQ-C. Key conclusions of the geochemical
analysis are as follows:
Ranges of formation geothermal gradients, water chemistry, subsurface pressure, and mineralogy were
evaluated with geochemical models to determine if there are suitable indicators for cement sealing
conditions in the subsurface. Potential indicators were analyzed using a CO2 batch reaction model and the
solid equilibrium phase model. Geochemical sealing conditions in the subsurface were also evaluated for
four test study areas with meta-modeling techniques.
Figure 7-1. pH and pe of the Appalachian Basin summary solutions, CO2 batch reactions.
Reference pH and pe are simulated for CO2 added to pure water for comparison.
Cement Phases
Figure 7-3. CO2 batch delta molar mass (moles of mineral phases that precipitated or dissolved
from previous CO2 batch step) and moles of mineral phases with midpoint summary solution
for Tuscarora sandstone, Rose Hill formation, and cement phases.
Michigan Basin
The change in precipitation of mineral phases with the addition of CO2 to the Michigan Basin midpoint
brine summary for the State Chester limestone, Dover 33 dolomite, Dover 33 A2 Carbonate, and cement
mineral solid equilibrium phases is shown in Figure 7-4. K-feldspar dissolves and quartz precipitates as
CO2 is added. The remaining minerals are at or near equilibrium, although, at a lower rate, dolomite
dissolves and calcite precipitates. In the Dover 33 dolomite, the dolomite precipitates and calcite and
halite dissolve with the addition of CO2. The minerals are then at or near equilibrium as CO2 nears
saturation. The remaining calcite then begins to dissolve (around 0.2 mol calcite dissolve in the final step
of the batch reaction). The Michigan Basin midpoint brine summary and Dover 33 A2 Carbonate solid
equilibrium phases shows that illite and quartz precipitate and K-feldspar dissolves. All minerals then
reach a state of equilibrium when 6.0 moles of CO2 added. The changes in the Michigan Basin cement
with the addition of CO2 are similar to the cement solid equilibrium cement phases analysis for the
Appalachian Basin, indicating that the initial brine solution is not an important factor when comparing the
midpoint values of the Michigan and Appalachian Basins.
Dover 33 Limestone
Dover 33 A2 Carbonate
Cement Phases
Figure 7-4. CO2 batch delta molar mass (moles of mineral phases that precipitated or dissolved
from previous CO2 batch step) (top) and moles of mineral phases (bottom), State Chester
limestone, Dover 33 limestone, Dover 33 A2 carbonate, and cement mineral phases
with midpoint summary solution. The reference pH and pe (Ref. pH and Ref. pe)
were found by reacting the CO2 batch with pure water alone.
Figure 7-5. Log SI values for cement minerals batch reaction for
pH, temperature, and pressure, Appalachian Basin.
Test Study Area 1. Test study area 1, located in the vicinity of Calhoun County, Michigan, is 7 x 7
km based on the parameters listed in Table 7-1, which were selected using wireline data and literature.
The study area encompasses 22 oil and gas wells that primarily target the Trenton-Black River
formations in the Albion-Scipio play. The potential storage zone in this study area is the Mt. Simon
sandstone, with the Eau Claire shale as the primary confining layer.
Test Study Area 2. Test study area 2, located in the vicinity of Saint Clair County, Michigan, was
estimated as 6 x 6 km based on the parameters listed in Table 7-1, which were selected using wireline
data and literature. The test study area encompasses 155 oil and gas wells that primarily target the
Niagara reef system. These Niagaran Reefs are the primary target for CO2, and this site would be a
CO2 EOR scenario where the depleted reefs would be filled with CO2. The reefs are fairly shallow at
approximately 975 m depth.
Test Study Area 3. Test study area 3, located in the vicinity of Muskingum County, Ohio, was
estimated as 15 x 15 km based on the parameters listed in Table 7-1, which were selected using
wireline data and literature. The test study area has 1,221 oil and gas wells that primarily target the
Clinton-Cataract group. There are 12 large CO2-emitting facilities nearby. The targeted storage
formations are the Copper Ridge dolomite down to the basal sandstone, with the Queenston shale to
the Black River group as the confining layers.
Test Study Area 4. Test study area four, located in the vicinity of Trumbull County, Ohio, was
estimated as 15 x 15 km based on the parameters listed in Table 7-1, which were selected using
wireline data and literature. The study area encompasses 357 oil and gas wells that primarily target the
Cataract-“Clinton” sandstone group. The potential storage zone in this study area is the Copper Ridge
dolomite down to the basal sandstone, with the Queenston shale to Black River group as the confining
layers.
Well construction specifications were also tabulated for wells at the four test study areas. Well locations,
completion date, total depth, casing schedule, cementing information, plugging details, and well status
details were compiled. The thickness of the cement in the production casing was also mapped out for test
study areas, because this directly relates to potential CO2 migration from the storage zone along the
boring. In addition, the thickness of the plugs was mapped out for plugged and abandoned wells.
Figure 7-7 shows maps of cement thickness outside of the production casing string. As shown, all sites
generally have 100 to 200 ft production casing cement, and many wells were cemented more than 500 ft
into overlying casing string. Production zone plugs are less thick and have a greater variation.
Subsurface Conditions. Subsurface conditions for the test study areas were determined for
caprock/reservoir mineralogy, pressure, temperature, salinity, and pH. This provides a range of conditions
to be examined with the uncertainty analysis and meta-modeling. Mineralogy was based on regional
trends in oil and gas reservoirs (Roan et al. 1996), oil and gas drilling records, and research on CO2
storage. Table 7-3 summarizes the general mineralogy for the caprock and reservoir zones for each study
area. Overall, the mineralogy would be expected to be fairly consistent across these local test study areas.
Brine geochemistry was based on regional studies on brine samples from the reservoir zones for each site
(Table 7-4). As shown, the sites have high salinity and moderate pH, which is typical for Paleozoic-age
rocks in the Midwest United States.
To help visualize the variations in subsurface conditions for wells at the test study areas, maps were
prepared for temperature and pressure. Temperature and pressure were based on well depths in each study
area. Subsurface temperatures were estimated with mean ambient temperature and 1° F per 100-ft
gradient. Pressures were based on a 0.444 psi/ft pressure gradient. Figure 7-8 shows the estimated
subsurface pressure distribution in wells. Since the pressure and temperatures were based on well depth,
the maps mainly reflect depths of the wells.
Meta-Model Development
To evaluate the cement sealing potential in wells at the four test study areas, meta-models were applied to
the sites based on subsurface pressure, temperature, and pH conditions present at the sites. The meta-
models portray how cement sealing conditions in the subsurface may vary across CO2 storage areas with
legacy wells. A meta-model is a model of a model, or a systematic method to portray a problem based on
metadata or input parameters. In this case, the metadata are the pressure, temperature, pH, fluid
chemistry, mineralogy, and cement makeup for the legacy wells. The meta-model approach assumes that
CO2 would be injected, contact the legacy wells near the well total depth, and introduce changes in
pressure, temperature, and pH (Figure 7-9). The meta-modeling predictions were constrained to
parameters in the space of interest. For example, the pressure ranges were constrained to initial pressure
and maximum injection pressure anticipated in the subsurface for each test study area. In this case, the
PHREEQ-C model output was represented by a proxy-model estimating CaCO3 saturation indices (SIs) at
various pressure, depth, and pH conditions. These scenarios represent conditions likely to be present in
the subsurface due to the injection of CO2.
The PHREEQC geochemical simulations results were processed to develop a proxy model to estimate
geochemical changes as a function of the depth, temperature, pressure, and introduced pH. The results
focused on the calcite SIs as an indicator of cement sealing potential. The proxy model was used to
develop a database of CaCO3 SIs for the wells at the four test study areas for different combinations of
pressure, temperature, and pH. For example, Figure 7-10 shows the meta-model output for the 357 wells
at test study area 4 at temperature +10 °C, pressure + 0%, and pH -2 from baseline conditions. The user
may adjust the ‘sliders’ for pressure, temperature, and pH to see how they affect the SI for the wells.
Thus, the meta-model is a way to visualize how differences in subsurface conditions may affect cement
sealing potential. In general, the meta-models showed only minor amount of variation in SI across the test
study areas. For example, the scenario shown in Figure 7-10 ranged from -4.9 to -5.2, a fairly small range.
Meta-models for the other test study areas show small changes for a combination of metadata. This
suggests that CO2-related cement sealing processes would not be sensitive to subsurface conditions.
The objective of the reporting and technology transfer task was to document project results and provide
project data to other CO2 research projects. Reporting and technology transfer activities included
preparation of task reports, technical meetings, project review meetings, and synergistic activities with
other carbon storage research projects.
8.1 Reporting
Summary reports were prepared for the major technical tasks and field work. Table 8-1 summarizes
reports generated and submitted to DOE-NETL. In addition, routine quarterly research performance
progress reports and financial reports were submitted to document technical and financial progress.
Table 8-1. Summary of deliverables.
Planned Submission
Task Milestone Description Deliverable
Due Date Date
1 Update Project Management Project Management Plan 30 days
October 6,
Plan after initial
2015
award
2 Complete Wellbore Integrity Well Integrity Registry Summary
June 2016 June 30, 2016
Registry Report
3 Collect Well Record Data Well Record Data Summary Report
June 2017 June 29, 2017
4 Complete Log & Testing Log & Testing Based Well Integrity
September November 1,
Based Well Integrity Assessment Summary Report
2017 2017
Assessment
5 Collect All SCP Analysis Compiled database of SCP
Data Analysis Data (uploaded to EDX
September 2018)
December
Appalachian Basin Field Testing March 2018
2017
Summary Report (March 2017)
Williston Basin Field Testing
Summary Report (February 2018)
6-7 Wellbore Integrity Sealing Wellbore Integrity Sealing Factor
June 2018 June 29, 2018
Factor Analysis Summary Report
8 Final Technical Report Final report with all project September September 27,
methods, results, and conclusions 2018 2018
Subsurface exposure to CO2 is a concern for wellbore integrity at CO2 storage sites, because CO2 can
corrode well materials and migrate along defects around the borehole. These processes may affect new
wells and legacy oil and gas wells. Consequently, this project completed a program to evaluate well
integrity in CO2 wells with a combination of direct field testing and analysis of well records. Project
results were used with geochemical analysis to identify trends that lead to better understanding and
prediction of well integrity issues for CO2 storage applications.
In this project, approximately 1,500 wells at three sites were reviewed in terms of well construction,
history of exposure to CO2, geochemistry, mineralogy, and well materials. The field sites included the
following locations:
• Appalachian Basin. The Appalachian Basin Indian Creek site is a natural CO2 and methane field
located in Kanawha County, West Virginia. The field contains approximately 58 wells at total
depths between 6,200 ft and 6,700 ft. The Indian Creek field produces in the Tuscarora sandstone.
The field has CO2 levels up to 60% in some areas.
• Michigan Basin. The Michigan Basin study site is located in the northern portion of the Niagaran
reef trend in Otsego County, Michigan. The fields were developed since the 1960s in the region,
and selected reefs have been subject to CO2 EOR since the 1990s. CO2 is also present in the
Antrim gas wells between 5% and 30% at depths between 1,000 ft and 1,500 ft. The Niagaran
reef EOR wells are completed at depths ranging from 5,000 ft to 7,000 ft.
• Williston Basin. The Williston Basin field is located on the northwestern edge of the Williston
Basin geologic feature and was a major oil play which is used for CO2 EOR. Approximately
3,000 wells are located in the Williston Basin testing site, completed at depth of 6,000 to 7,000 ft.
The Williston Basin testing site is a mature oil field that began production in 1954. CO2 EOR was
started in 2000 at the site, expanding to additional areas over time.
A total of 83 CO2 wells were surveyed at the Michigan Basin (23 wells) and Williston Basin (60 wells)
sites for wellhead casing pressures that may indicate well defects. The Appalachian Basin site was not
available for field testing, because the asset was sold to a new operator. Detailed SCP testing was
completed on 23 wells that had some indication of significant SCP.
The test results did not show significant well defects, with casing pressures less than 1 MPa and minor
pressure buildup patterns. The wells demonstrated zonal isolation, with no CO2 gas migrating to the
wellhead b-annulus. Analysis of wells with minor pressure buildup was inconclusive in terms of nature
and severity of the defects. These results were surprising, given that 15-20% of typical oil and gas wells
develop casing pressure. It was expected that more defects would be present in CO2 wells, because of
carbonic acid evolution and corrosive conditions in the subsurface.
Additional geochemical modeling and meta-modeling for the three field sites and four test study areas
were completed to determine if subsurface conditions at the field sites were suitable for cement sealing of
gas migration pathways via CaCO3 precipitation. Well construction and/or geochemical conditions for
cement carbonation appear to have contributed to well integrity. Results support management of CO2
storage applications in areas with many legacy oil and gas wells.
Some key conclusions of the integrated wellbore integrity analysis program include the following:
• The wellbore integrity registry developed in this project provides a catalog of the well
component, integrity issues, causes, timing, and leakage pathways that may occur in wells. Most
wellbore integrity problems are in the casing, cement, or interface between the two components
or arise due to geological processes such as formation lithology and geomechanical stresses.
Aadnoy, B. S. (1990). In-Situ Stress Directions from Borehole Facture Traces. Journal of Petroleum
Science and Engineering, 4: 143-153
Avary, K. L. (1996). The Lower Silurian Tuscarora sandstone fractured anticlinal play. The atlas of major
Appalachian gas plays, Roen, J. B., & Walker, B. J. (Eds), West Virginia Geological and
Economic Survey Publication, 25: 151-155.
Azaroual, M., Durst, P., Czernichowski-Lauriol, I., Olsen, D., Stentoft, N., Springer, N., Rochelle, C. A.,
Pearce, J., Bateman, K., & Birchall, D. (2004). The geochemical reactions resulting from CO2
injection into the midale formation, Weyburn oilfield; A laboratory experimental and modelling
study. GHGT7: 7th International Conference on Greenhouse Gas Control Technologies,
Vancouver, Canada, 5-9 September 2004, 4 p.
Battelle. (2015). Development of Subsurface Brine Disposal Framework in the Northern Appalachian
Basin. Report prepared for Research Partnership to Secure Energy for America under contract
#11123-73, 412 p.
Battelle. (2016). Final Technical Report: Systematic Assessment of Wellbore Integrity for Geologic
Carbon Storage Projects Using Regulatory and Industry Information. Prepared for U.S.DOE-
NETL under contract DE-FE009367.
Bell, J. S., & Gough, D. I. (1979). Northeast-southwest compressive stress in Alberta: Evidence from oil
wells, Earth Planet, Sci. Lett., 45: 475-482.
Braunberger, J., Bremer, J., Liu, G., Gorecki, C., Peck, W., Steadman, E., & Harju, J. (2012).
Characterization, Petrography, and Static Modeling of an Unconventional Carbonate Reservoir:
Intervals of the Midale and Rival “Nessons” Beds in the Mississippian Madison Group, Burke
County, North Dakota. 012 AAPG Annual Convention and Exhibition, 22-25 April 2012, Long
Beach, California.
Breen, K., Angelo, C., Masters, R., & Sedam, A. (1985). Chemical and isotopic characteristics of brines
from three oil- and gas-producing sandstones in eastern Ohio, with applications to the
geochemical tracing of brine sources. Water-Resources Investigations Report 84-4314. Prepared
in cooperation with the Ohio Department of Natural Resources, Division of Oil and Gas.
Cantucci, B., Montegrossi, G., Vaselli, O., Tassi, F., Quattrocchi, F., & Perkins, E. H. (2009).
Geochemical modeling of CO2 storage in deep reservoirs: The Weyburn Project (Canada) case
study. Chemical Geology, 265: 181-197.
Carey, J. W. (2013). Geochemistry of Wellbore Integrity in CO2 Sequestration: Portland Cement-Steel-
Brine-CO2 Interactions. Reviews in Mineralogy and Geochemistry, 77(1): 505-539.
Carey, J. W., Wigand, M., Chipera, S., WoldeGabriel, G., Pawar, R., Lichtner, P., Wehner, S.,
Raines, M., & Guthrie Jr., G. D. (2007). Analysis and performance of oil well cement with
30 years of CO2 exposure from the SACROC Unit, West Texas, USA. International Journal of
Greenhouse Gas, 1(1): 75-85.
Carroll S., Carey. J. W., Dzombak, D., Huerta, N., Li, L., Richard, T., Um, W., Walsh, S., & Zhang, L.
(2016). Review: Role of Chemistry, Mechanics, and Transport on Well Integrity in CO2 Storage
Environments. International Journal of Greenhouse Gas Control, 49: 149-160.
Chamberlain, B., Chiesl, N., Day, J., Dowd, L., Overocker, B., Pape, D., Petrus, M., Swanson, M., &
Toles, J. (1995). Concrete - A Material for the New Stone Age. University of Illinois, Materials
Science and Technology, A MAST Module. Retrieved from <illinois.edu>.