Скачать книгу

current flow. EIT methods are highly sensitive to changes in electrical resistivity. In EIT methods, a known amount of electrical current is injected in the medium of the model and resulting electrical potential field is measured around the boundary points of the medium of the model. From these data, it is possible to perform electrical resistivity inversion of these potential field data to determine the electrical conductivity or resistivity inside the medium of the model that is being probed by the currents. In this way, the internal resistivity or electrical conductivity distribution is reconstructed using electrical measurement from the boundary of the medium of the model under investigation. Such electrical resistivity or conductivity distribution gives valuable information about the interior of the medium. Basically, all EIT methods deal with solving forward problem and inverse problem iteratively. EIT methods are used for various shallow depth archaeological prospecting, geothermal resource prospecting, geo-environmental monitoring, hydrogeological invitations, and geotechnical investigations etc. A good account of EIT methods and electrical resistivity imaging tools may be found in many geophysical literatures [13–20]. In brief, all these inverse problems are solved iteratively until a best possible solution is obtained.

      Hohmann and Raiche [21] have given a good description and classification of model properties in terms of relevance, non-relevance, or irrelevance. The importance and relevance of model properties can also be described theoretically using eigenvalues and its size. The “singular value decomposition” (SVD) of a matrix has been found useful in the least-squares electrical inversion methods. Tripp et al. [22] have given an excellent overview of this technique in electrical resistivity inversion. It is important to note that the SVD technique needs some special efforts that are computationally demanding and time consuming.

      Another computationally attractive technique is described in the electrical resistivity inversion, which uses a “damping factor”. Different criteria are described in the research literature for the selection of “damping factor”. Narayan et al. [9] have given a good account of the various electrical resistivity inversion schemes and methods. The main benefit to use a damping in inversion methods is that it produces a stable solution. Constable et al. [23] have introduced an Occam’s inversion, which is receiving attention in the electrical geophysics. Without knowing magnitudes of the eigenvalue, it is difficult to quantify an appropriate amount of damping. Oristaglio and Worthington [24] have described a method and criterion for selecting the amount of damping. It is important to note that the selection of the amount of damping is highly debatable. A different approach of selecting relatively high damping in the initial part and gradually low in the later part of the inversion is introduced by Eaton [25]. It may only be done based on a genuine experience.

      There are two important steps for the inversion of electrical data as shown below:

       A fast numerical modelling scheme, which computes electrical data theoretically for a given numerical/physical model (i.e., a fast “forward modelling scheme”), and

       An efficient numerical approach to compute first derivatives of the data (also called “coefficient matrix”).

      Unfortunately, the second requirement is not easily ready for 2-D or 3-D inverse geophysical and geological models. Furthermore, there is no claim for an optimum inverse technique or approach for the interpretation of 2-D and/or 3-D geophysical data. The purpose of this article is to introduce an effective solution for this problem and to propose multi-dimensional electrical resistivity inversion techniques.

      Based on “reciprocity theorem (1961)” and “perturbation analysis”, a new method for the multi-dimensional resistivity inversion is introduced herein. This approach for inverse formulation is entirely different from resistivity inversion methods published in the geophysical literature. The main benefits of the proposed method are listed here:

       The first derivatives of the data with respect to the model properties are computed efficiently.

       It provides a mathematical picture of sensitivity analysis. It is expressed as an amount of power that is lost in the anomalous zone.

      The governing equations (also known as “Poisson” equation) for potential field distribution due to a 3-D direct current (DC) point source Is(x, y, z) are given by [26]

      and

      (2.4)images

      where

      (2.6)images

      After expanding equation (2.5) completely and upon neglecting second-order terms here, it gives