higher permeability (Fig. 2.4). The correlation of this high permeability feature with a northwest trending break in seismic topography suggested that the carbon dioxide injected into well KB‐501 might flow preferentially along a fault zone on the flank of the anticline defining the field (Vasco et al., 2008).
The range changes observed over injector KB‐502 displayed a distinct double‐lobed pattern seen in Figures 2.3 and 2.5 (Vasco et al., 2010). As noted by collaborators from Pinnacle Technologies, such patterns were often observed in tilt‐meter monitoring of hydrofracturing, suggesting the tensile opening of a steeply dipping planar feature (Davis, 1983). For a tensile feature, such as a fracture, the fluid pressure induces the fracture to open. Thus, the volume change is primarily due to displacement normal to the fracture plane, a change in fracture aperture (fracture width). On the basis of Pinnacle's suggestion, a vertical to subvertical fault/fracture zone model was constructed. The zone extended 80 m above and below the reservoir and approximately 6 km to the northwest and 6 km to the southeast of the injection point in the plane, trending on an azimuth of 135 degrees. In addition, volume change within the reservoir layer was also accounted for in the fashion described above (Vasco et al., 2010).
Figure 2.4 (a) Onset time of the most rapid change in reservoir volume. The curves indicate the flow paths from the well to points within the reservoir. (b) Logarithm of the permeability multiplier found by solving equation 2.20 along the trajectories.
Figure 2.5 (a) Detailed view of range change above well KB‐502 following 1,060 days of injection. The darker colors signify range decrease associated with uplift above the well. (b) Seismic horizon displaying pushdown, most likely due to velocity decreases associated with the migration of injected carbon dioxide.
A homogeneous elastic model of the overburden did not produce the correct range change estimates, and hence an inversion based upon a model with uniform properties led to erroneous depth estimates for a tensile source. Thus, elastic layering, derived from well logs, had to be included in the modeling. The total range change from 12 July 2003 through 19 March 2007 was used to infer the distribution of volume change within the reservoir and the distribution of aperture change over the fault/fracture zone. It was found that the range changes could be matched by a combination of reservoir volume change and cumulative tensile opening of a fault/fracture zone confined to lie at depths within 100 m of the reservoir (Vasco et al., 2010). The lateral extent of the fracture opening was much greater, extending over 3 km from the injection point. Aspects of the fault/fracture zone were subsequently supported in an analysis of three‐dimensional surface seismic data. In particular, Gibson‐Poole and Raikes (2010) noted that time shifts, thought to be due to the injection of carbon dioxide and the resulting seismic velocity change, followed a remarkably linear zone with near parallel boundaries located between the two lobes of range change (Fig. 2.5). The orientation of the seismic feature agreed rather well with the azimuth of 135 degrees required to fit the range‐change data. These conclusions were subsequently supported by the work of Zhang et al. (2015, 2016).
X‐Band InSAR and Multicomponent Displacement Data
The Envisat C‐band satellite was deorbited on October 2010 and was no longer available for monitoring the range changes over the injectors at In Salah. Fortunately, several other satellite systems were functional as replacements before this date. In particular, the COSMO‐SkyMed (CSK) X‐band satellite pair were launched and operational before the loss of the Envisat satellite. Similarly, the TerraSAR X‐band satellite was also put in orbit, with a repeat cycle of 11 days. Thus, three X‐band satellites were accessible for estimating range changes at the In Salah site due to the injection of carbon dioxide. The repeat time attainable for these three satellites was 8 to 11 days, much better than the best possible Envisat return time of 35 days. The down side of this transition to X‐band data was the increased cost of data acquisition because the satellites have to be tasked to gather data at a particular site. However, with the help of British Petroleum, we were able to acquire X‐band SAR scenes from both ascending and descending acquisition geometries from 2009 until 2012.
Due to the favorable geometry of the two‐satellite systems, it was possible to obtain both quasi‐east‐west and quasi‐vertical components for all time intervals during this period. An example is shown in Figure 2.6 for displacements above well KB‐501. Note that there is a hint of a double‐lobed pattern in the vertical component of displacement. As at KB‐502, this pattern suggests the possible opening of a subvertical tensile feature, such as a fault or damage zone. The existence of such a damage zone is supported by the narrow, linear, high‐permeability trend shown in Figure 2.4, found by the diffusive travel time inversion of Vasco et al. (2008). Furthermore, a recent inversion of a decomposition of C‐band data into quasi‐vertical and quasi–east‐west components found that is was not possible to fit both data sets simply with reservoir volume change. A contribution from a subvertical tensile feature provided an acceptable fit to the two components (Rucci et al., 2013). Based upon these results, a steeply dipping (86o) tensile damage‐zone model was adopted to explain the surface displacement observed in Figure 2.6. The subvertical damage zone, with an orientation of 135o, was divided into a 20 x15 grid of surface elements. Each grid element could undergo a distinct aperture change in order to fit the two components of data. The resulting model (Fig. 2.7) contains an elongated region of aperture change from the depth of injection to about 100 m above the well. Significant aperture change extends roughly 6 to 7 km to the northwest of the injection well and about 5 km to the southeast.
The aperture change estimates, shown in Figure 2.7, provide a reasonable fit to the InSAR observations, even reproducing the small secondary lobe at the eastern edge of the vertical displacement anomaly (Fig. 2.8).
2.4.2. Aquistore, Canada
The Aquistore CO2 Storage Project (Worth et al., 2014) began in 2012 in southeastern Saskatchewan, Canada (Fig. 2.9). The site serves as buffer storage for the world's first commercial, postcombustion, carbon capture plant that began operation in October 2014 at the Boundary Dam coal‐fired power station. Injection started in April 2015 with a maximum anticipated injection rate of up to 1,500 tonnes/day. The carbon dioxide is sequestered in the brine‐filled clastic strata of the Deadwood and Winnipeg formations at 3,150 – 3,350 m depth, the deepest sedimentary units in the Williston Basin (Norford et al., 1994). Similar deep saline formations are found throughout the world. These formations are well suited for the storage of carbon dioxide because they are porous and permeable, have