# Interpretation¶

## Apparent Resistivities¶

Plotting apparent resistivities, as pseudosections or in plan-view maps, is informative. The images are useful for recognizing data outliers, bad electrodes, and validating normalizations that might have been applied to the data. In addition, they sometimes provide valuable information about earth structure.

As an illustration, Fig. 180 (a) shows an earth model consisting of a single prism (\(\rho=10\; \Omega \cdot m\)) buried in a uniform halfspace (\(\rho= 100\; \Omega \cdot m\)). A dipole-dipole survey is carried out along a line that passes directly above the conductive prism. The resulting pseudosection is shown in Fig. 180 (b). The prism is manifested as a region of lower apparent resistivity in the center of the image and there are wings extending outwards and downward. The apex of the image can be used to estimate the horizontal location of the prism but the depth to the body is less evident since the vertical scale of the pseudosection is in n-values and not in meters.

Despite the above success, the situation worsens if the earth is more complex. This is illustrated in Fig. 180 (c) where some near-surface inhomogeneities are added. The same dipole-dipole survey is carried out and resultant pseudosection is shown in Fig. 180 (d). The response of the prism is masked and attempting to infer existence and location of the prism is extremely challenging.

This example can be downloaded here.

Gradient array surveys are often used in reconnaissance modes and it is insightful to repeat the above analysis with a representative example. A plan view of the resistivity model and electrode geometry is shown in Fig. 181 (a). The survey consists of a grid of 13 x 13 receivers located between a 450 meter dipole current source. Each receiver is a 20 meter dipole. The corresponding apparent resistivity map is shown in Fig. 181 (b). An estimate of the horizontal location of the center of the prism can be obtained, but again there is no quantitative information about the depth.

Contaminating the model by adding some conductive and resistive features (Fig. 181 (c)) leads to an apparent resistivity map that is very complicated and in which the signal of the prism is masked (Fig. 181 (d)). To obtain information about the electrical conductivity we need to invert the data.

## Inversion of DCR data¶

The data are the measured voltages and the goal is to find an electrical conductivity that satisfactorily reproduces these data and agrees with a priori geologic structure and petrophysical constraints.

To illustrate the importance of inverting the data, we take a two spheres setup, one resistive, one conductive. Although the geology is 3D, we first invert the data using a 2D inversion algorithm. Parameters used for the inversion of the dipole-dipole data (Fig. 182 (b)) are provided in Table 13.

Number of sources | 43 |
---|---|

Number of data | 540 |

Data uncertainties | \(2\%|d| (percentage) + 2 \times 10^{-5} V\) (floor) |

Mesh Size | \(10 \times 10 \times 10\) meters |

Reference conductivity | \(0.01\) S/m |

Regularization Scales ( \(\alpha_s, \alpha_x,\alpha_y,\alpha_z\) ) | \(0.01, 1, 1, 1\) |

Fig. 182 (c) presents the recovered 2D conductivity model after convergence of the algorithm.

**Important comments:**

- Even though there are no contaminating near-surface blocks the pseudosection does not clearly indicate two bodies. This is in contrast to Fig. 180 (a) where a single prism was clearly identified in the pseudosection.
- The two spheres are recovered but they have lower conductivity contrasts with respect to the halfspace than do the true spheres. This occurs for three reasons: (i) the inversion generates smooth models and this extends structures and reduces amplitudes; (ii) the true spheres extend into regions where there is limited depth of investigation; and (iii) the 2D inversion assumes that the structures are cylindrical.

Similar to the prism model example (Fig. 180), we repeat the experiment with the same survey setup but use a more complicated resistivity model that has near-surface inhomogeneities (Fig. 182 (d)). The resulting pseudosection (Fig. 182 (e)) is challenging to interpret visually. The 2D resistivity model recovered from the inversion ( Fig. 182 (f)) unravels the data complexity.

**Important comments:**

- The pseudosection of data is complicated and dominated by the near-surface conductors.
- The inversion recovers the contaminating surface conductors. It also recovers the two spheres with about the same fidelity as in the simple case.

This example can be downloaded here.

## Depth of Investigation¶

The inverse problem is nonunique and the DCR data are sensitive to conductivity only in a region in the vicinity of the electrode array. Conductivity structures that exist outside this region are unreliable and likely artifacts of the inversion. There are several methods proposed in the literature to quantify this limits of this region for a specific DCR survey. The following example uses the Depth of Investigation (DOI) method proposed by LiDWO1999.

In its simplest form, the DOI analysis requires the data to be inverted twice with all parameters the same except for the background conductivity. For the two-sphere example shown in Fig. 182 (c), the synthetic data is inverted a second time with a reference halfspace conductivity of \(10\; \Omega \cdot m\). Fig. 183 (a) shows the recovered 2D resistivity model. Note that the region away from the electrode locations returns to a value close to the reference model.

We now have a discretized volume of the Earth and two conductivity models that
can equally reproduced the observed data. Let \(\sigma_1, \sigma_2\) be
the conductivity values recovered at some location (*x,z*), a DOI index is calculated as:

where the DOI index will approach 1 for similar model values obtained with both inversions regardless of the chosen reference models \(\sigma_1^{ref}, \sigma_2^{ref}\). Conversely, the DOI will approach 0 where the recovered models return to their respective reference conductivities. Fig. 183 (b) presents the calculated DOI index for the two- sphere problem, showing a lower confidence over the bottom half of the domain. The DOI mask is applied to our preferred 2D model, presented in Fig. 183 (c), with transparency applied proportionally to the DOI index.

### Downloads¶

Data, model and inversion files used in this page can be downloaded below:

Utilities: UBC-DC2D data viewer and model viewer

### References¶

Oldenburg, Douglas W., and Yaoguo Li. “Estimating depth of investigation in DC resistivity and IP surveys.” Geophysics 64.2 (1999): 403-416.