# 2D Model Study of CO2 Plumes in Saline Reservoirs by Borehole Resistivity Tomography

Authors: Said A. al Hagrey
Journal: International Journal of Geophysics
 (2011)

---

## Abstract

The performance of electrical resistivity tomography (ERT) in boreholes is studied numerically regarding changes induced by CO2 sequestration in deep saline reservoirs. The new optimization approach is applied to generate an optimized data set of only 4% of the comprehensive set but of almost similar best possible resolution. Diverse electrode configurations (mainly tripotential α and β) are investigated with current flows and potential measurements in different directions. An extensive 2.5D modeling (>100,000 models) is conducted systematically as a function of multiparameters related to hydrogeology, CO2 plume, data acquisition and methodology. ERT techniques generally are capable to resolve storage targets (CO2 plume, saline host reservoir, and impermeable cap rock), however with the common smearing effects and artefacts. Reconstructed tomograms show that the optimized and multiply oriented configurations have a better-spatial resolution than the lateral arrays with splitting of potential and current electrode pairs between boreholes. The later arrays are also more susceptible to telluric noise but have a lower level of measurement errors. The resolution advance of optimized and multiply oriented configurations is confirmed by lower values for ROI (region of index) and residual (relative model difference). The technique acceptably resolves targets with an aspect ratio down to 0.5.


---

## Body

## 1. Introduction

The need to manage the global CO2 emissions for mitigating the greenhouse effect has led to a world wide research to reduce atmospheric CO2. Techniques of carbon capture and storage (CCS) must (1) be effective and cost-competitive, (2) provide stable, long-term storage, and (3) be environmentally benign. Potential terrestrial media for CO2 storage include depleted oil and gas reservoirs, unmineable coal seams, and deep saline water reservoirs capped by impermeable rock to prevent upward leakage.CO2 exists in the gas phase at standard atmospheric temperature and pressure. Above the dynamic critical point (>31.1°C, >7.38 MP, density >0.469 g/cm3), CO2 changes to a supercritical fluid phase; it diffuses through solids like a gas and dissolves material like a liquid. CO2 has long been injected in the subsurface to enhance oil, gas, and coal-bed methane recovery and storage. This injection has been mainly monitored using seismic time-lapse imaging (e.g., Sleipner oil field in North Sea, e.g., [1]). However, investigations on brine-saturated sandstones showed that the electrical resistivity (ρ) is more sensitive to CO2 saturation than is seismic velocity (Figure 1). This may justify application of electrical resistivity tomography (ERT), particularly in boreholes, for monitoring resistive supercritical CO2 plumes in a deep saline reservoir (e.g., [2]). This reservoir formation normally consists of a highly resistive matrix (e.g., sandstone and limestone) and a conductive pore brine. Here, CO2 saturation can be predicted using the law of Archie [3]:(1)ρ=aρwΦ-mSw-n′,(2)ρCO2=aρwΦ-m(1-SCO2)-n′,
where ρ, ρw, ρCO2 = bulk, fluid, and CO2 resistivity, respectively, Φ = porosity, Sw, SCO2 = water, and CO2 saturation, a, m, n′ = constants.Figure 1
Experimental P-wave velocity (vp) and electrical resistivity (ρ) of sandstone reservoir versus brine saturation (S) showing ρ far more sensitive to S than vp [4].Recent developments have enabled installing deep boreholes with coated (insulating) casing and fixed electrode arrays for ERT monitoring (e.g., [5]). Forward modeling and inversion algorithms have also been developed for better monitoring CO2 plume scenarios in deep saline reservoirs and coal seams (e.g., [6, 7]). In 2008 we started the interdisciplinary project “CO2 MoPa” (modeling and parameterization of CO2 storage in deep saline formations for dimensions and risk analysis). It aims at studying long-term CO2 attenuation and migration in deep and shallow layers (including saline and fresh water aquifers), along with assessing storage capacity and analyzing risk. Various synthetic, almost realistic, storage scenarios are simulated for formations of the North German Basin that seem suitable for CO2 storage. Our main task is to develop optimized, constrained monitoring strategy techniques for CCS using a combined seismic and ERT approach. This approach focuses on (1) developing a constrained electrical resistivity-modeling strategy based on a priori subsurface knowledge from seismic time-lapse imaging (some years) and logging data, and (2) reliably inverting continuous ERT time-lapse imaging to yield spatiotemporal developments in the intrinsic physico-chemical properties of CO2 reservoirs and cap rocks with time.


### 1.1. Optimized and Reliable Borehole ERT

Inverse ERT algorithms, however, tend to smear resistivity values from any given voxel to adjacent voxels (e.g., [8]). Based on the sensitivity functions, ERT in boreholes performs near the electrodes (boreholes) better than in the interwell region. Electrical sensitivity is used to select an array for a certain target but one that does not necessarily has the best possible resolution. Thus, a new approach of array optimization was recently developed to search for electrode configurations that maximize survey resolution (e.g., [9]). The optimization algorithms take into account the trade-off between the spatial and temporal (measurement time) resolution. They select measurements based on the contribution to the cumulative sensitivity of the array (e.g., [10]) or the model resolution matrix, R (e.g., [11]). R depends on sensitivities of all configurations plus regularization types used in the inversion [12]. For an arbitrary electrode array, the algorithms generate optimized (opt) data sets that have far less size than the comprehensive one and almost the same resolution of the targets. This comprehensive set includes all possible viable electrode configurations conducted within this array and possesses the maximum possible resolution, see next sections. The application of this 2D array optimization was recently extended into borehole-borehole and surface-borehole surveys [13–15]. The last algorithm is applied here. This algorithm strongly improves the ERT resolution in the interwell region (commonly low) to approach the resolution of the highly sensitive region close the boreholes.The sensitivity accounts only for data sampling and model heterogeneities. As an alternative, the region of investigation index (ROI) is used to assess the whole 2D inversion procedure such as the data sampling and noise, model discretization and regularization, and nonlinearity [16, 17]. Thus, the reliability of ERT 2D tomograms will be evaluated here by ROI in addition to the relative model difference (residual) between each input and inverted output model.For ERT in borehole surveys, the topographical aspect ratio (AR) is defined by the vertical length of the electrode array divided by the horizontal crosshole offset. Thus, resolution is enhanced by increasing the density (i.e., number and thus costs) of expensive monitoring wells in the area of plume migration. Newmark et al. [18] studied AR values of 2, 1.5, and 1 and found that they (in this order) show the best, intermediate, and worst resolution, respectively. In this study the lower boundary of AR is extended down to 0.25 which leads to a further decrease in the number (and thus costs) of monitoring wells. A broad AR range (0.25–2) is tested to determine its optimum value between the highest and lowest resolution (of AR = 2 and 0.25, resp.) that corresponds to the highest and lowest number of monitoring wells (i.e., costs), respectively.


### 1.2. Problem and Objectives

Until now ERT is rarely applied for the CCS problematic in deep saline reservoirs. Only few recent studies partly treated this problem, for example, feasibility studies by Christensen et al. [19] and sensitivity investigations for some specified CO2 plume forms using point and long (metal-cased) borehole electrodes by Ramirez et al. [20]. Also there is a deficit of field sites for CO2 sequestration that are equipped by adequate acquisition infrastructures for ERT surveys. All these lead to a strong demand for systematic ERT modeling investigations. In this study, extensive, systematic numerical ERT 2.5D modeling is carried out, and the results are analyzed for different virtual scenarios of injected wedge-like CO2 plumes (dimensions, SCO2 or ρ) as a function of electrode configuration, burial depth, AR, data noise, and setup parameters of modeling constraints (mainly regularization parameters, see next sections). Moreover, ROI analyses and residuals are applied to evaluate the resolving capability of various electrode configurations and inversion procedures. The technique’s robustness in the field is tested by adding three different random errors to data sets. These studies aim to test the capability of (non-)standard and optimized ERT techniques (partly developed here) to resolve the subsurface CO2 storage targets as a function of diverse parameters related to hydro-/geologic and geochemical subsurface properties (mainly of saline reservoir and cap rock), CO2 plume, survey design, data acquisition, and modeling techniques.In the next sections, I describe the applied borehole electrode configurations, the experiment setup of the subsurface model scenarios, and modeling varieties. I then discuss, summarize, and conclude the results of the different numerical simulations. Optimization algorithms, noises, depth effect, and modeling constraints are not shown here. They are contained in Hagrey [2, 14, 15] and Hagrey and Petersen [21].


## 1.1. Optimized and Reliable Borehole ERT

Inverse ERT algorithms, however, tend to smear resistivity values from any given voxel to adjacent voxels (e.g., [8]). Based on the sensitivity functions, ERT in boreholes performs near the electrodes (boreholes) better than in the interwell region. Electrical sensitivity is used to select an array for a certain target but one that does not necessarily has the best possible resolution. Thus, a new approach of array optimization was recently developed to search for electrode configurations that maximize survey resolution (e.g., [9]). The optimization algorithms take into account the trade-off between the spatial and temporal (measurement time) resolution. They select measurements based on the contribution to the cumulative sensitivity of the array (e.g., [10]) or the model resolution matrix, R (e.g., [11]). R depends on sensitivities of all configurations plus regularization types used in the inversion [12]. For an arbitrary electrode array, the algorithms generate optimized (opt) data sets that have far less size than the comprehensive one and almost the same resolution of the targets. This comprehensive set includes all possible viable electrode configurations conducted within this array and possesses the maximum possible resolution, see next sections. The application of this 2D array optimization was recently extended into borehole-borehole and surface-borehole surveys [13–15]. The last algorithm is applied here. This algorithm strongly improves the ERT resolution in the interwell region (commonly low) to approach the resolution of the highly sensitive region close the boreholes.The sensitivity accounts only for data sampling and model heterogeneities. As an alternative, the region of investigation index (ROI) is used to assess the whole 2D inversion procedure such as the data sampling and noise, model discretization and regularization, and nonlinearity [16, 17]. Thus, the reliability of ERT 2D tomograms will be evaluated here by ROI in addition to the relative model difference (residual) between each input and inverted output model.For ERT in borehole surveys, the topographical aspect ratio (AR) is defined by the vertical length of the electrode array divided by the horizontal crosshole offset. Thus, resolution is enhanced by increasing the density (i.e., number and thus costs) of expensive monitoring wells in the area of plume migration. Newmark et al. [18] studied AR values of 2, 1.5, and 1 and found that they (in this order) show the best, intermediate, and worst resolution, respectively. In this study the lower boundary of AR is extended down to 0.25 which leads to a further decrease in the number (and thus costs) of monitoring wells. A broad AR range (0.25–2) is tested to determine its optimum value between the highest and lowest resolution (of AR = 2 and 0.25, resp.) that corresponds to the highest and lowest number of monitoring wells (i.e., costs), respectively.


## 1.2. Problem and Objectives

Until now ERT is rarely applied for the CCS problematic in deep saline reservoirs. Only few recent studies partly treated this problem, for example, feasibility studies by Christensen et al. [19] and sensitivity investigations for some specified CO2 plume forms using point and long (metal-cased) borehole electrodes by Ramirez et al. [20]. Also there is a deficit of field sites for CO2 sequestration that are equipped by adequate acquisition infrastructures for ERT surveys. All these lead to a strong demand for systematic ERT modeling investigations. In this study, extensive, systematic numerical ERT 2.5D modeling is carried out, and the results are analyzed for different virtual scenarios of injected wedge-like CO2 plumes (dimensions, SCO2 or ρ) as a function of electrode configuration, burial depth, AR, data noise, and setup parameters of modeling constraints (mainly regularization parameters, see next sections). Moreover, ROI analyses and residuals are applied to evaluate the resolving capability of various electrode configurations and inversion procedures. The technique’s robustness in the field is tested by adding three different random errors to data sets. These studies aim to test the capability of (non-)standard and optimized ERT techniques (partly developed here) to resolve the subsurface CO2 storage targets as a function of diverse parameters related to hydro-/geologic and geochemical subsurface properties (mainly of saline reservoir and cap rock), CO2 plume, survey design, data acquisition, and modeling techniques.In the next sections, I describe the applied borehole electrode configurations, the experiment setup of the subsurface model scenarios, and modeling varieties. I then discuss, summarize, and conclude the results of the different numerical simulations. Optimization algorithms, noises, depth effect, and modeling constraints are not shown here. They are contained in Hagrey [2, 14, 15] and Hagrey and Petersen [21].


## 2. Borehole Electrode Configurations

Similar to surface surveys, ERT data acquisition between two borehole electrode arrays can be conducted in the tripotential 4-pole configurationsα (CPPC, C = current electrode, P = potential electrode), β (CCPP) and γ (CPCP), and their reciprocals. The γ configurations can be derived from α and β measurements, that is, they are not independent and are usually excluded from the data [22]. The rest of α and β measurements are accomplished in vertical (v, at 90°), horizontal (h, 0°), and lateral (l, >0–<90°) modes (Figure 2). These 4-pole modes are carried out within the same borehole (inhole) or distributed between the two boreholes (crosshole). Table 1 shows all possible configurations α and β for the inhole and crosshole modes [23]. Each configuration consists of two inhole and three crosshole arrangements.Table 1
All possible 4-pole tripotential configurationsα and β (nonequivalent, nonreciprocal) of a survey between two borehole electrode arrays.
GroupConfigureInholeCrosshole4-00-43-11-32-2αCPPC-00-CPPCCPP-CC-PPCCP-PCβCCPP-00-CCPPCCP-PC-CPPCC-PPFor example, group 4-0 denotes for number of electrodes (current, C/potential, P) in the first (4) and second (0) boreholes, respectively.Tripotential electrode configurationsα and β applied for electric resistivity tomography in boreholes (inhole and crosshole) using 4-pole of current (C) and potential (P) electrode pairs. The survey can be conducted in circulating (c) vertical, v (a), lateral, l and horizontal, h (b) modes and tripole-pole, tp-p (c) which is a special type of circulating vertical configuration with fixed C electrodes [24].
(a)(b)(c)For an array ofN electrodes, the whole comprehensive data set contains [N(N-1)(N-2)(N-3)/8] independent nonreciprocal 4-pole configurations (Table 2, [25]). Excluding the less stable inversion configurations from this whole set results in a more effective data set, simply called here comprehensive data set. These redundant configurations include γ configurations and those with geometric factors larger than that produced by the dipole-dipole array where the maximum dipole separation is 6𝔞 (𝔞 = electrode spacing, Figure 2). The comprehensive data set with all viable configurations should provide the best possible resolution [26]. It contains all subsurface information that can be gathered by an N-electrode array. Table 2 displays the size of the whole and effective comprehensive data sets in comparison with the standard ones for a collinear 32 electrode array and circulating arrays between two boreholes, each of 16 electrodes. The effective comprehensive set in boreholes contains more than 80,000 data points and is even about 0.65 of its whole set but more than 450 times each of the standard dipole-dipole and Wenner sets. This justifies applying the new approach of electrode optimization to generate optimized data sets of far lower size and almost the same resolution as the comprehensive ones, that is, of highly spatiotemporal resolution (for more details, see [14, 15]). Briefly this opt array improves the ERT resolution particularly in the interwell region which is commonly low compared to that of the region close the boreholes.Table 2
Data set size of 4-pole configurations for a collinearN (32) electrode array and circulating crosshole array of 16 electrodes in each borehole.
Data setRelationCollinear  (N=32)Circulating crosshole array (N = 16↔16)Comprehensive, whole (independent, nonreciprocals)N(N-1)(N-2)(N-3)/8  (α + β + γ) data107,880122,760Comprehensive, effective  (whole compreh.—redundant)(α + β) data  (no γ + no noisy)68,62180,738Dipole-dipoleβ (𝔞=1, n = 1–6)((N-3)+(N-8))×3159165Wenner,α (CPPC)(N-1)(N-2)/6155165Redundant data include configurations ofγ and noisy data of very low voltage.Noisy data are those with geometric factor larger than that of dipole-dipole configuration of maximum dipole factor of 6.Ten different borehole surveys of standard, nonstandard, and optimized configurations are investigated here (Table3, Figure 2). The vertical (v) circulating (c) configurations αvc and βvc represent the comprehensive data sets with all possible α and β electrode combinations, respectively. Their corresponding subconfigurations αvcs (s = symmetrical around the midpoint) and βvcs include only the conventional symmetrical arrangements of Wenner and Schlumberger, and dipole-dipole, respectively. The configurations βl and βh represent lateral (l) and horizontal (h) bipole-bipole configurations, where their C- and P-pairs are split between the two boreholes (CP-CP). The configuration βl represents the comprehensive data set acquired in all possible electrode combinations and orientations of this bipole-bipole (CP-CP) arrangement. Its subset βh is conducted with horizontal current electrodes (flows) only.Table 3
Applied configurations for two vertical arrays, each of 16 electrodes (see Figure2). The corresponding electrodes in each borehole are set at equal depths with unit interval.
ConfigurationData no.Characteristics(1) tp-p813Vertical current flow, horizontal resolution(2)βh1,240Horizontal/lateral current flow, better data quality, better vertical (less lateral) resolution(3)βl14,440(4)βhtp2,053Sum ofβh and tp-p ((1) + (2))(5)αvcs1,796Sum of Wenner and Schlumberger(6)βvcs1,177Dipole-dipole(7)αβvcs2,973Sum ofαvcs and βvcs ((5) + (6))(8)αvc34,788Comprehensive including nonstandard and standard (αvcs and βvcs in αvc and βvc, resp.) configurations(9)βvc31,510(10) opt3,000Optimized data set with less than 4% size but 97% resolution relative to that of comprehensive setc: circulating, v: vertical,α, β: tripotential configuration, l: lateral, s: symmetrical, tp-p: tripole-pole, h: horizontal, opt: optimized.The tripole-pole (tp-p) is a special type of vertical circulating crosshole configurations [24]. It has fixed C electrodes (the respective upper- and lowermost electrode of the first and second hole, and vice versa; Figure 2) and moving P electrodes between any other possible combinations of electrode pairs in each borehole separately (i.e., CPP-C, C-PPC, PPC-C, and C-CPP of Table 1). It shows near-vertical current flows that are strongly influenced by horizontal layers in the interwell region. As opposed to tp-p, bipole-bipole configurations (βl and βh) have mostly lateral current flows, that is, a low resolution for horizontal structures, but detect better vertical structures. The new complex data sets βhtp (sum of βh and tp-p) and αβvcs (sum of αvcs and βvcs) are introduced for the first time in this study. Each should reflect the advantages of its constituting arrays. Their combined current flows and potential measurements in various possible combinations and directions should be able to resolve targets of varying orientations. The applied opt data set (3000 data points) is less than 4% of the comprehensive set but has an average relative resolution of 0.96% (Table 3).


## 3. Subsurface Model Scenarios

To keep the virtual CO2 sequestration modeling more realistic, the formation parameters of the starting subsurface scenarios used to generate the synthetic data have been taken from published data, for example, CO2 SINK test site of Ketzin, near Berlin (e.g., [5, 27, 28]). The applied single subsurface models consist of the electrically almost insulating, supercritical CO2 plume (ρCO2≈∞) sequestrated at the top of a conductive saline sandstone reservoir (ρreservoir = 3 Ωm, ρbrine = 0.20 Ωm, salinity 35–55 g/L, Φ = 20–25%), (Figure 3). This reservoir is capped by an impermeable siltstone (ρ = 8 Ωm with varying thicknesses in the range 2𝔞−9.5𝔞 (𝔞 = electrode spacing)) to prevent upward CO2 leakages. Model dimensions in this study are given in unit of 𝔞 which is often assumed to be 1 m. The CO2 plume is simulated by the common wedge shape with bulk resistivities, ρplume of 100, 30, 15, and 10 Ωm (corresponding to saturations, SCO2 of 80, 60, 40, and 30%, resp., as calculated from (2), and varying thicknesses (0.5𝔞–13𝔞) and widths (0.5𝔞–13𝔞)) [29].Figure 3
Subsurface model scenarios of saline sandstone reservoir (resistivity,ρ = 3 Ωm) with CO2 wedge-like plume of varying width (d), thickness (h), depth and saturation (80, 60, 40, and 30%, i.e., ρ = 100, 30, 15, and 10 Ωm, resp.) capped by siltstone (8 Ωm) together with borehole electrode arrays (●) used in the survey. d1, d2, h1, h2 are given in units of electrode spacing 𝔞.


## 4. Applied Procedures

The 2.5D forward and inverse ERT modeling of deep CO2 plume was carried out using new codes based on algorithms for shallow surveys (e.g., [30]). All codes use the half-space solution with a fine mesh grid to accurately model the whole region. The Neumann and mixed boundary conditions are used for the top surface and the side/bottom boundaries, respectively. The program optimizes automatically a sufficient number of additional mesh grids far from the electrodes (5𝔞–10𝔞) at these side/bottom boundaries.An extensive numerical investigation was started by generating synthetic data sets of apparent resistivity (ρa) from 2.5D forward simulation as a function of (1) ten wedge-like models (Figure 2), (2) four SCO2 values, (3) ten electrode configurations (Table 3, Figure 1), (4) seven ARs, (5) two burial depths (1𝔞 and 101𝔞) and (6) three noise levels (1%, 3%, and 5%). This forward modeling has resulted in more than 8000 synthetic data sets for these diverse model scenarios with and without the wedge-like CO2 plume. Each synthetic data set was filtered to remove any outliers and any ρa values of high geometric factor which may cause a potential leakage. The filter ensures also the absence of any equivalent, reciprocal, or γ configurations from the data set. Each filtered data set has been inverted independently eight times using different setup constraint parameters. This results in more than 100,000 tomograms for all applied synthetic data sets. The diverse setup constraints applied in the inversion mainly include regularizations with the minimization methods of least squares (L2) or robust blocky normalization (L1), and initial models of a constant homogeneous resistivity or an approximate inverse model.The reliability of the forward modeling using finite-element method (with the more accurate trapezoidal elements instead of triangle ones) was confirmed before starting inversions using a homogeneous medium with a constant resistivity. Compared with this constant resistivity, the resulting deviations of the single apparent resistivities (ρa) for each applied synthetic data set generally do not exceed 3%. These deviations are similar to the normal error level in the real data and are considered here as a noise in the synthetic data sets.Among the 10 examined wedge-like scenarios, this study focuses mainly on models W1, W42, and W8 with the largest (13𝔞), intermediate (5.5𝔞), and least (0.5𝔞) thicknesses for the CO2 plume (Figure 3). The electrode coverage for the cap rock is least for W1 (only 1𝔞), intermediate for W42 (5𝔞), and best for W8 (11𝔞). These CO2 storage scenarios represent examples with optimum target dimensions (W42) and problems of thin layers (plume in W8 and reservoir below the plume in W1) and thin widths (lowermost triangle apex in W1). Among the eight independent inversions for each data set, only the one with the best-fitting model (least root mean square error, rms-error) is considered for further studies.In the following, any inverted tomogram with its single target anomalies (CO2 plume, host reservoir rock, and cap rock) is evaluated relative to the starting input model according to the following criteria:(1)
the reconstructed geometry, shape, and position,(2)
the recovered resistivity magnitude,(3)
the sharpness of the boundaries.


## 5. Results and Discussions

In the following, the reconstructed ERT tomograms for W1, W42, and W8 scenarios will be described and discussed only as a function of the applied CO2 plume scenarios (dimensions, SCO2 or ρ), electrode configurations, ARs and noises.The overall average rms-error for all inverted data sets approaches 0.9% with a nearly similar distribution between tomograms inverted by the robust L1 norm and their corresponding models of the L2 norm. For all studied data sets, every best-fitting tomogram (least rms-error almost of <0.5% and average iteration number of 5), among its independent eight inversions, was optimized with the L1 norm. This low rms-error value is explained by the good convergence of the synthetic data sets toward the final solution. The average rms-error values are least for the lateral/horizontal configurations (βl/βh), intermediate for the opt and vertical circulating symmetrical configurations (αvcs, βvcs, αβvcs), and highest for the comprehensive (αvc, βvc), tp-p and βhtp. One may note that the misfit distribution in these arrays is directly proportional to the data number with singularity problems and/or high geometric factors (i.e., of voltage leakages).In all studied cases, the L1 tomograms always show sharper boundaries and better magnitude recovery for the single targets than the L2 tomograms, that is, L1 inversion fits better for resolving sequestration targets with sharp boundaries. These results are in accordance with the evidence that the L1 norm tends to produce models that are piecewise constant, whereas the L2 norm tends to smear out the sharp boundaries [19].Briefly every best resolved output tomogram (among its eight independent inversions) resulted from the set of inversion constraint parameters which incorporates the use of (1) accurate calculation options at the expense of computation time (e.g., the standard Gauss-Newton algorithm to (re-)calculate the Jacobian matrix and optimization, 4 nodes per𝔞 and mesh grid to reduce singularity errors, the model refinement of half width and crosshole model of half size, etc.), (2) a robust blocky L1 norm for sharp boundaries instead of smooth L2 norm for gradual changes, (3) sufficient boundary meshes (>5𝔞), and (4) low damping parameters for these synthetic data (almost noise free). Based on these discussions, all (best-fitting) tomograms considered throughout the next sections were inverted using the robust L1 norm.


## 6. Effect of Configurations and Model Scenarios

Figure4 shows examples of the reconstructed 2D tomograms (with the best-fitting least rms-error) for the subsurface scenarios W1, W42, and W8 of different CO2 plume dimensions and only the two extreme SCO2 of 80 and 30 g/L (corresponding ρ of 100 and 10 Ωm, resp., Figure 3). These scenarios are situated at a depth of 101𝔞 (to the uppermost electrode). They are reconstructed for the applied 10 electrode configurations (Table 3, Figure 2) using the aspect ratio (AR) of 1.Figure 4
Numerical resistivity tomograms (3rd–12th column) inverted for different wedge-like CO2 plumes (W1, W42, W8, 1st-2nd column, see Figure 3) of varying saturations—resistivity of 100 (a) and 10 Ωm (b)—and 4-pole electrode configurations (1st row, see Figure 2 and Table 3 for symbols) between two borehole electrode arrays (●). All models show root mean square errors of less than 1%. Solid lines refer to target boundaries.In most applied cases, the resistive wedge-like CO2 plume together with the conductive saline reservoir of hosting sandstone and the impermeable siltstone cap rock are generally mapped directly by the inverted absolute ρ tomograms (no model differencing between tomograms with and without plume). Most configurations generally reconstruct these three targets with varying degrees of smearing and artifacts. This smearing is reflected on the single tomograms by target anomalies of lower magnitude, larger volume, and blurred boundary relative to that of their input models. This smearing is evaluated quantitatively by the ROI and residual analyses, see next section. The mapping capability generally decreases with decreasing dimensions (thickness, width) and resistivity of the targets. With the exception of the lateral/horizontal (βl/βh), all other configurations can reconstruct even the worst CO2 plume scenario (of the least thickness of W8 and resistivity of 10 Ωm). Based on the previously mentioned criteria for evaluating the inverted tomogram anomalies (geometry, position, ρ amplitude, and boundary sharpness), the mapping capability generally is better for configurations with multiply orientated current flows and potential measurements (αvc and βvc, their subsets αvcs, βvcs, and αβvcs, as well as opt and βhtp) than those of only lateral (βl), horizontal (βh), and vertical (tp-p) current injection. Tomographic anomalies of the later configurations generally show lower ρ magnitude than that of the former configurations, that is, configurations βl, βh and tp-p underestimate the ρ magnitude compared with the input models, see next section.Among the applied scenarios, the mapping capability for the sequestration targets (cap rock, plume, and reservoir) is best for W42, intermediate for W1, and least for W8. Tomogram W1 shows strong smearing effects for the lowermost apex of plume triangle and the cap rock with their surroundings. Both targets have the thin layer problem with electrode coverage of ≤1 which causes this smearing governed by the equivalence principle (e.g., [15]). For a thin resistive CO2 plume target, the inversion code is not able to resolve thickness and resistivity parameters of certain resistive model features individually but only their product (resistance). The horizontal boundaries of the upper reservoir/plume and lower plume in W42 and partly W1 are generally resolved slightly better than the inclined boundaries of the plume.Regarding the applied configurations, the vertical circulatingαvc and αvcs generally resolve inclined boundaries better than their corresponding βvc and βvcs. This may explain the occasional moderate resolution of αβvcs tomograms (containing βvcs data) for the lower triangle apex of the plume in W1. It is clear that the individual complex configurations βhtp (sum of βh and tp-p) and αβvcs (sum of αvcs and βvcs) combine the features (mostly of better resolution) of their corresponding single constituents. Each complex configuration contains more data measured in more orientations and thus carries more information than its constituent data sets (cf. [31]). On the other hand, models of the vertical circulating β (βvc and βvcs) tomograms generally reflect higher (better) resistivity magnitudes of the targets than their corresponding α (αvc and αvcs) models. The resistivity magnitude approaches its real value in the target centre and deviates increasingly from this with distance toward the contact with the next target. As opposed to real data, the inversion models do not show the usual poor resolution with strong artifacts around the boreholes. These are due to the heterogeneities caused by boring and electrode installation and/or the distorting effects of the conductive borehole fluid relative to the resistive host rock [32].Regarding the data set size, the resolution of opt tomograms (nearly 3,000 data points) is almost similar to the best possible resolution of each of the comprehensiveαvc and βvc models (each of >30,000 points) and far better than that of βl tomograms (14,440). Obviously, the resolution for the noncomprehensive data sets is best for opt, above moderate for the complex αβvcs and βhtp, moderate for αvcs and βvcs, and least for βh and tp-p.Briefly, the lateral/horizontal configurations (βl and βh) are more robust against measurement errors due to current and voltage leakages. In these configurations, splitting the current and potential electrode pair between boreholes (i.e., CP-CP) maximizes the measured voltage (e.g., [23]). Large borehole offset (i.e., P1-P2 distance), however, may include the telluric noise in the data (e.g., [33]). This noise can be minimized by periodically reversing the current flow in the current electrodes. These arrays result in poorly resolved anomalies, especially for lateral boundaries, often with underestimated ρ magnitudes. This is due to the predominance of current flows and potential measurements in the lateral directions with poor data coverage in the vertical one. Configurations having either a C or P electrode pair or both in the same borehole (inhole) such as most of the vertical circulating α and β (αvc, βvc, αvcs, and βvcs, Table 1) configurations have the singularity problem of very low voltages. However, their filtered data sets after removing this noise can resolve well the subsurface targets of CO2 plume, reservoir, and cap rock assuming enough coverage of more than one electrode spacing. These conclusions are in accordance with that of, for example, Oldenburg and Li [16] and Oldenborger et al. [34] and opposed to that of Zhou and Greenhalgh [23]. In solute transport experiments, the former authors found that the vertical βvcs tomograms show more reliable subsurface structures than the horizontal β tomograms. Bing and Greenhalgh stated that the acquisition data for vertical inhole configurations are easily obscured by background noise and yield images inferior to those from lateral/horizontal configurations βl and βh.


## 7. Evaluation of Reconstructed Tomograms

This evaluation is carried out by the ROI analysis and the model difference (residual) relative to the input model. The ROI analysis is started by conducting two independent inversions of the same dataset using different resistivity values (qA and qb) for homogeneous reference models. Values of qA and qb are calculated from the average logarithmic ρa using two different multiplication factors. The inversion results from these two starting models are used to calculate the ROI value for each pixel, defined as [35]:
(3)ROI(x,z)=[qA(x,z)-qB(x,z)](qA-qB)-1.The ROI approaches zero when the model is well constrained by the data (the two inversions reproduce very similarρ values) and one when the model has a very poor data coverage. The two inversions of each ROI analysis were conducted using the same set of inversion setup parameters (resulted in the best fitting tomograms, see before) and only three iterations. Higher iteration numbers were tested and resulted in artifacts as the algorithm tries to reduce the data misfit by modeling the noise as well. Different pairs of multiplication factors (0.1 and 10, 0.1 and 10, and 0.2 and 5) were applied here and yield similar results (see also [34, 36]).On the other hand, the relative resistivity difference or residual (Δρ) between each corresponding pixel of the input (ρinput) and output (ρoutput) 2D models is calculated by
(4)Δρ=(ρoutput-ρinput)(ρoutput+ρinput)-1.Moreover, a random noise at 1, 2, and 5% levels was added to the synthetic data sets, in addition to their forward modeling errors (up to 3%), to evaluate noisy field effects on the inversion. Adding noise in this order generally increases the rms-error values by a factor of 2 to 9 but slightly decreases the mapping capability, particularly for W1 and W42. Ramirez et al. [20] obtained similar results; the effect of the random error is insignificant for anomalies of a large size and magnitude.Figure5 shows examples of the resulting ROI and Δρ tomograms as a function of the applied electrode configurations carried out for the subsurface scenarios W1, W42, and W8 at 101𝔞 depth only. Constraining of inverted tomograms by data coverage is best (lowest values for ROI and partly Δρ) for the comprehensive αvc and βvc, and opt arrays, intermediate for the complex αβvcs and βhtp, and poor (highest ROI values) for the others. This confirms the effectiveness of our optimization approach applied here to generate a practical opt dataset of high resolution.Figure 5
Evaluation of the reconstructed electrical resistivity tomograms (depth of 101𝔞, 𝔞 = electrode spacing) using methods of region of index, ROI (top, (3)) and relative model difference, Δρ (bottom, (4)) for survey between two borehole electrode arrays (●) as a function of subsurface scenario (W1, W24, and W8, 1st column), and electrode configuration (see Table 3 for symbols). Solid lines refer to target boundaries.Unlike real field data, this synthetic modeling analysis does not show the usual disadvantageous low resolution of high values of ROI andΔρ (of bad data and thus resolution) around the boreholes due to the heterogeneities resulting from boring and electrode installation. The highest values of ROI and Δρ with the least data coverage and resolution are concentrated in the central interwell region which is a common disadvantage for ERT results. Also ROI and Δρ analyses reflect similar results between the corresponding single tomograms of shallow (1𝔞 depth, not shown here) and deep scenarios (101𝔞 depth). This similarity reflects the reliability of the applied techniques. A comparison with published results (e.g., [34]) shows that the vertical configurations reflect better resolution (of lower ROI values) than the lateral ones which confirms the results obtained here.Briefly, tomograms ofΔρ and ROI generally show similar results regarding the mapping capability of the single applied arrays. They reflect well the common smearing effects of varying degrees. These effects lead to overpredicted volumes, underpredicted magnitudes, and blur boundaries of the target anomalies. This study shows clearly that the complex (βhtp and αβvcs) arrays with multiply oriented measurements in addition to opt arrays with practical data sizes are recommended for highly spatiotemporal resolution and will be considered further in this study.


## 8. Effect of Tomographic Aspect Ratio (AR)

Extensive tests for seven different AR steps within 0.25–2.0 range were performed as a function of the applied subsurface scenarios (Figure2), SCO2, burial depths, electrode configurations (Figure 1, Table 3), noises, and modeling setup parameters. At the AR values (0.5, 0.7, 1.0, 1.3 and 2), Figure 6 shows some best-fitting (least rms-error) tomograms resulting from unconstrained inversions only for the highly resolving arrays (βhtp, αβvcs and opt) and the subsurface scenario W42 with SCO2 of 60% (ρCO2=30 Ωm).Figure 6
Effect of aspect ratios, AR (first row) on ERT tomograms, resulting from wedge-like CO2 plume scenario W42 (Figure 3, second row) with 30 Ωm plume resistivity for 4-pole configurations βhtp, αβvcs and opt (see Figure 2 and Table 3 for symbols) between two borehole electrode arrays (●). All models show root mean square errors of less than 1%. Solid lines refer to target boundaries.Figure6 shows that the obtained mapping capability (including resolution) generally increases with increasing AR from 0.5 to 2. This occurs although the lateral boundaries of the CO2 wedge plume become more vertical with increasing AR. This mapping improvement with increasing AR is also associated with a slight increase of the recovered ρ magnitude of the CO2 plume relative to that of the starting input model. Most other results of the previous sections are manifested here such as the poor mapping resolution of thin layers (of W1 and W8 scenarios, not shown here). At AR = 1, the lowermost triangle apex of W1 is predominated either from the reservoir anomaly or from a wide smeared zone of the plume. At AR>1, the apex resolution increases with increasing AR and approaches best results at AR of 2. Based on the applied electrode configurations, the resolution for tomograms of the vertical circulating and opt arrays is better than for those of lateral/horizontal βl/βh and partly tp-p and βhtp configurations (partly not shown here). Compared with W1 and W8, W42 models are generally better resolved due to the better coverage of their target. The ability to detect and often map the three sequestration targets (CO2 plume, reservoir, and cap rock) by unconstrained inversions is still possible with AR values down to 0.5 for the most studied scenarios (even those with the worst scenario of least thickness and ρ). This result is superior to that of published studies (e.g., [18]). These authors applied ERT techniques for site characterization and process monitoring and determined a minimum AR with acceptable output resolution of 1. In comparison, the minimum AR value (0.5) currently obtained in this study can lead to a decrease of the number of the expensive monitoring wells and the costs by a factor of 4 (i.e., n(1/AR), n: well number for AR = 1). The reconstructed output tomograms for lower AR values (<0.5) achieve a satisfactory resolution only for constrained inversions with an a priori fixing of boundaries and/or resistivities of the targets. The resolution increases with increasing the number of constraints.


## 9. Summary and Conclusions

Electrical resistivity tomography (ERT) techniques in boreholes are powerful in monitoring intrinsic property changes for storing the resistive (supercritical) CO2 in conductive saline reservoirs. In this study, the mapping capability of various ERT techniques is studied for diverse wedge-like CO2 plumes in a deep saline aquifer capped by an impermeable rock. Extensive, systematic 2.5D modeling studies (>100,000 models) were performed to test the ERT sensitivity for a multitude of parameters related to the subsurface setting (hydrogeology and geochemistry of reservoir and cap rock), CO2 plume reservoir, survey design, data acquisition, and modeling techniques. The new array optimization approach is applied to generate optimized data sets (opt) of only 4% of the comprehensive set but of almost similar resolution. Forward simulations were carried out to generate diverse synthetic data sets (>8000) as a function of plume scenarios (different dimensions and CO2 saturations SCO2 or resistivity, ρ), burial depths, electrode configurations, random noise, and aspect ratios (AR). The data quality (<3% noises) is confirmed by results of tests on a homogeneous model with constant ρ. This numerical study principally reveals the capability of ERT techniques to resolve the various deep subsurface scenarios with the CO2 sequestration targets (plume, host reservoir, and cap rock). Most important results may be summarized as follows:(1)
Most applied ERT configurations can generally map the sequestration targets of sharp boundaries directly by the absoluteρ tomograms (no model differencing) using L1 robust inversions. All models, however, reflect smeared anomalies of lower magnitude, larger area, and blurred boundary.(2)
Superior to published studies, the detection of CO2 targets is possible even for the worst scenario of 0.5𝔞 thicknesses (𝔞 = unit electrode length), 30%  SCO2, and 0.5 AR. At lower AR values (<0.5), a satisfactory resolution can result only from constrained inversions with an a priori fixing of boundaries and/or resistivities of targets.(3)
The developed opt and complex (αβvcs, βhtp) arrays (nearly 3000 data points) are recommended for surveys of highly spatiotemporal resolution. Their resolution is the second best after the comprehensive vertical arrays (αvc and βvc, each of >30,000 data) and far better than the comprehensive lateral array (βl, 14,400 data). These arrays improve the common low resolution of the ERT technique in the interwell region by combining configurations with current flows and potential measurements in all possible orientations.(4)
Lateral arrays (βl and βh) are more robust against measurement noise; their synthetic data sets result in tomograms with the least rms-error misfit, but their real field data may include natural telluric noise at large borehole offset.(5)
Analyses of the region of index, ROI, and residual model difference relative to the input model show that inverted tomograms of comprehensive (αvc and βvc), opt, and complex (αβvcs) configurations are better constrained by the data than those of the other applied ones. This tomogram reliability generally increases by increasing the data size.(6)
Sometimes the configuration tripole-pole (tp-p) is able to detect horizontal structures due to its near-vertical current flows. Contrary to this, lateral arrays (βl and βh) have lateral current flows and thus better resolution for vertical structures.(7)
The vertical circulatingα and β data sets are collected partly with either or both of the current and potential electrode pair inhole and thus may show singularity problems. Filtering these data sets results in tomograms of better resolution than those of lateral arrays.(8)
Adding random noise (1%, 3%, and 5%) to the synthetic data (in addition to its forward modeling errors, up to 3%) increases the rms-error values (according to the error continuation law) by a factor of 2–9 but slightly decreases the mapping capability of the techniques, particularly for large targets.


## 10. Outlook and Recommendation

The current results give answers to some studied problems of 2D ERT mapping for CO2 sequestration in deep saline aquifers. Many other problems related to the 2D/3D mapping and 4D monitoring are currently studied within the research activities of our MoPa project. These tasks include the following:(i)
Continuing the systematic 2.5D modeling studies for other CO2 plume scenarios and extending these to 3D investigations with varying parameters related to modeling, data acquisition and methodology, geological setting and plume reservoir developments based on petrophysical approaches.(ii)
Applying approaches of time-lapse imaging for 4D monitoring of any CO2 migration either laterally and downwards within the saline reservoir or upwards through probable postinjection fracturing of the cap rock caused by the high injection pressure. This includes monitoring of pre- and postinjection scenarios, developments of the new CO2 reservoir, and any change in the porosity/permeability with time.(iii)
Risk analyses for any possible postinjection fracturing in the cap rock and upward CO2 seepages into the near surface zone and especially freshwater aquifers.(iv)
Evaluating the spatiotemporal resolution of output tomograms using techniques of residual, 2D ROI, and 3D VOI. Quantitative analyses ofρ magnitude and spatial extension in the single pixels/voxels of the recovered tomograms with respect to the pre-injection model.(v)
Applying developed modeling techniques on more realistic subsurface scenarios by using available (meta)data of the North German Basin andρa (hard data) inversion constrained by seismic and log information (soft data) to reduce the problems of nonuniqueness. Constrained inversion schemes should be able to fix resistivity regions and include boundaries in the output models.(vi)
Using petrophysical approaches to quantify inverted resistivity models in saturation, and the quality of the pore filler (brine and CO2).(vii)
Refining and verifying the use of optimized configurations for enhancing the spatial resolution without a significant decline of temporal resolution.(viii)
Testing the developed techniques on real field data, for example, the sequestration site of Ketzin in collaboration with partners of the project CO2 sink.



---

*Source: 805059-2011-12-15.xml*