Fréchet Derivatives calculation for electrical resistivity imaging using forward matrix method

Document Type : Research Article

Authors

1 Assistant Professor, Institute of Geophysics, University of Tehran, Tehran, Iran

2 Ph.D Graduated, Institute of Geophysics, University of Tehran, Tehran, Iran

Abstract

Fréchet derivatives calculation or sensitivity matrix is an integral part of every non-linear inversion process. The sensitivity values indicate the variation of the forward response with respect to the variation of model parameters. Sensitivity patterns are also a criterion to assess the reliability of inverted models and to design optimum resistivity surveys. In this study, a numerical approach based on the forward matrix calculation in the framework of the 2.5D finite difference electrical resistivity forward modeling is presented. First, using the potential distribution in the Fourier space obtained from the forward calculation and the derivatives of the coupling coefficients with respect to the conductivity distribution, the sensitivity values in the wavenumber domain are computed. Then, these values are transformed into the space domain using an inverse Fourier technique. To verify and analyze the proposed numerical method, the sensitivity distributions assuming the homogeneous and inhomogeneous media for commonly used electrical resistivity tomography configurations (e.g. pole-pole, pole-dipole, dipole-dipole, and the Wenner arrays) are computed. The numerical experiments reveal that the sensitivity patterns vary spatially throughout the model depending not only on the resistivity distribution but also on the electrode configuration. It is also concluded that the sensitivity analysis can be used as a supplementary tool for any optimum electrical tomography survey design. 

Keywords

Main Subjects


Dey, A., and Morrison, H. F., 1979, Resistivity modelling for arbitrarily shaped two-dimensional structures: Geophysical Prospecting, 27(1), 106-136.
Edwards L. S., 1977, A modified pseudosection for resistivity and IP. Geophysics, 42, 1020–1036.
Geselowitz, D. B., 1971, An application of electrocardiographic lead theory to impedance plethysmography: IEEE Transactions on Biomedical Engineering, 18(1), 38–41.
Ghanati, R., Azadi, Y., and Fakhimi, R., 2020, RESIP2DMODE: A MATLAB-based 2D resistivity and induced polarization forward modeling software: Iranian Journal of Geophysics, 13(4), 60-78.
Greenhalgh, S. A., Zhou, B., Greenhalgh, M., Marescot, L., and Wiese, T., 2009, Explicit expressions for the Fréchet derivatives in 3D anisotropic resistivity inversion: Geophysics, 74, F31-F43.
Günther, T., Rücker, C., and Spitzer, K., 2006., Three-dimensional modelling and inversion of DC resistivity data incorporating topography – II, Inversion: Geophys. J. I., 166, 506–517.
Loke, M. H., and Barker, R. D., 1995, Least-squares deconvolution of apparent resistivity pseudosections: Geophysics, 60, 1682–1690.
McGillivray, P. R., and Oldenburg, D. W., 1990, Methods for calculating Fréchet derivatives and sensitivities for the non-linear inverse problem: a comparative study: Geophysical Prospecting, 38(5), 499–524.
Park, S. K., and Van, G. P., 1991, Inversion of pole-pole data for 3-d resistivity structure beneath arrays of electrodes: Geophysics, 56, 951–960.
Spitzer, K. 1998, The three-dimensional DC sensitivity for surface and subsurface sources: Geophys. J. I., 134(3), 736–746.
Szalai, S., and Szarka, L., 2008, Parameter sensitivity maps of surface geo-electric arrays I. Linear arrays: Acta Geodaetica et Geophysica, 43, 419–437.
Wiese, T., Greenhalgh, S., and Marescot, L., 2009, DC resistivity sensitivity patterns for tilted transversely isotropic media: Near Surface Geophysics, 7(2), 125–139.
Zhou, B., and Greenhalgh, S. A., 1999, Explicit expressions and numerical computation of the Fréchet and second derivatives in 2.5D Helmholtz equation inversion: Geophysical Prospecting, 47, 443–468.