Influence of noise estimation on Electrical Resistivity Tomography Data Inversion

Document Type : Research Article

Authors

1 Department of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: yosraazadi@ut.ac.ir

2 Corresponding Author, Department of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: rghanati@ut.ac.ir

Abstract

Electrical resistivity tomography is a simple, cost-effective, and highly practical method for surveying near-surface properties. Today, this method is widely used in the discovery and exploitation of water resources, archeology, and environmental and hydro-geophysical studies (such as estimating the hydrogeological parameters of the aquifer). In electrical resistivity imaging, according to the purpose and location of data collection, the electrodes are placed in specific arrays, and data collection is performed. The collected data (potential distribution or apparent resistivity) is then transformed into a distribution of actual electrical resistivity values using inverse modeling methods. Imaging requires defining and solving a nonlinear inverse problem. In this strategy, we optimize the objective function, which consists of fitting field and theoretical data. First, the physics of the problem (forward model) is presented by solving Poisson's equation with the finite difference numerical solution method. An accurate and efficient forward calculation is the basis of most processes of the inversion. Calculation of resistivity forward responses is carried out using simulation of the current flow into the earth’s surface through solving Poisson’s equation. In this contribution, a finite-difference algorithm is applied to discretize the simulated models, restricted by a mixed boundary condition. One of the merits of the finite-difference method over the other methods is its well-known ability to quickly approximate the solutions for any arbitrary and complex structure models. The finite-difference method is relatively fast compared with the finite-element method. However, to include a general topography, the finite-element method becomes a better selection despite being computationally expensive. The partial differential equations governing the resistivity problem are obtained by using the principle of conservation of charge and the continuity equation.
The inverse problem is then solved by linearizing the problem in different iterations. A significant part of this research is how to perform inverse modeling of electrical resistivity data. The formulation and solution of the forward and inverse problem in this dissertation have been programmed in MATLAB and part of the program has been written in the C language to increase the computing speed. The field data is noisy due to the non-ideal measuring instruments, improperly filed conditions, operator errors, and geological conditions. Noise values can play a pivotal role in the inversion of electrical resistivity due to the special properties of the inverse problem. A proper estimation of field measurements noise level prevents over- or under-fitting of the calculated data and field data during inversion. Improper fitting (i.e., fitting where the value of the parameter  is much more or less than one) leads to creating an artifact or loss of important details in the final inverted model. In this paper, to deal with the effect of noise level on the ERT inversion results, two methods of reciprocity error method and stacking error method have been used. The results of numerical modeling show that the appropriate estimation of the noise level leads to the estimation of subsurface resistivity models close to the ground reality. We also provide a comparison between the inversion results obtained with the presence of noise level and those derived without including the weighting matrix into the objective function.

Keywords

Main Subjects

References

Backus, G., & Gilbert, F. (1968). The resolving power of gross earth data. Geophysical Journal International, 16(2), 169-205.
Claerbout, J.F., & Muir, F. (1973). Robust modeling with erratic data. Geophysics, 38(5), 826-844.
DeGroot-Hedlin, C., & Constable, S.C. (1990). Occam’s inversion to generate smooth, two-dimensional models from magneto-telluric data. Geophysics, 55, 1613-1624.
Dahlin, T. (1996). 2D resistivity surveying for environmental and engineering applications. First break, 14(7), 275-283.
Dey, A., & Morrison, H.F. (1979a). Resistivity modeling for arbitrarily shaped two‐dimensional structures. Geophysical Prospecting, 27(1), 106-136.
Dey, A., & Morrison, H.F. (1979b). Resistivity modeling for arbitrarily shaped three-dimensional structures. Geophysics, 44(4), 753-780.
Edwards L.S. (1977). A modified pseudosection for resistivity and induced-polarization. Geophysics, 42, 1020-1036.
Fallah Safari, M., & Ghanati, R. (2022). DC Electrical Resistance Tomography Inversion, Journal of the Earth and Space Physics, 47(4), 87-98.
Ghanati, R., & Fallahsafari, M. (2022). Fréchet Derivatives calculation for electrical resistivity imaging using forward matrix method, Iranian Journal of Geophysics, 15(4), 153-163.
Griffiths D.H., & and Barker R.D. (1993). Two-dimensional resistivity imaging and modelling in areas of complex geology. Journal of Applied Geophysics, 29, 211-226.
Habberjam, G.M. (1967). Short note: On the application of the reciprocity theorem in resistivity prospecting. Geophysics, 32, 918.
Jackson, D.D. (1972). Interpretation of inaccurate, insufficient and inconsistent data. Geophysical Journal International, 28(2), 97-109.
Kemna, A., Binley, A., Cassiani, G., Niederleithinger, E., Revil, A., Slater, L., Williams, K.H., Orozco, A .F., Haegel, F.H., Hoerdt, A., & Kruschwitz, S. (2012). An overview of the spectral induced polarization method for near-surface applications. Near Surface Geophysics, 10(6), 453-468.
LaBrecque, D.J., Mletto, M., Daily, W., Ramirez, A.L., & Owen, E. (1996). The effects of noise on Occam’s inversion of resistivity tomography data. Geophysics, 61, 538.
Loke, M.H. (1994). The inversion of two-dimensional resistivity data. Ph.D. thesis, University of Birmingham.
Loke, M.H., & Barker, R.D. (1994). Rapid least-squares inversion of apparent resistivity pseudo-sections. 54th EAEG Meeting, Vienna, Austria.
Loke, M.H., & Barker, R.D. (1995). Least-squares deconvolution of apparent resistivity pseudo-sections. Geophysics, 60, 1682-1690.
Loke, M.H., Acworth, I., & Dahlin, T. (2003). A comparison of smooth and blocky inversion methods in 2D electrical imaging surveys. Exploration Geophysics, 34(3), 182-187.
Loke, M.H., Chambers, J.E., Rucker, D.F., Kuras, O., & Wilkinson, P.B. (2013). Recent developments in the direct-current geo-electrical imaging method. Journal of applied geophysics, 95, 135-156.
McGillivray, P.R. (1992). Forward modeling and inversion of DC resistivity and MMR data. Ph.D. thesis, University of British Columbia.
Oldenburg, D. W., & Li, Y. (1994). Inversion of induced polarization data. Geophysics, 59(9), 1327-1341.
Oldenburg, D.W., & Li, Y. (1999). Estimating depth of investigation in dc resistivity and IP surveys. Geophysics, 64(2), 403-416.
Park, S.K., & Van, G.P. (1991). Inversion of pole-pole data for 3-D resistivity structure beneath arrays of electrodes. Geophysics, 56, 951-960.
Sasaki, Y. (1992). Resolution of resistivity tomography inferred from the numerical simulation. Geophysical Prospecting, 40, 453-464.
Schlumberger, C., & Schlumberger, M. (1929). Electrical Logs and correlations in Drill Holes. Mining Metallurgy, 10, 515-518.
Seigel, H., Nabighian, M., Parasnis, D.S., & Vozoff, K. (2007). The early history of the induced polarization method. The Leading Edge, 26(3), 312-321.
Smith, N.C., & Vozoff, K. (1984). Two-dimensional DC resistivity inversion for dipole-dipole data: IEEE Trans. Geosci. Remote Sensing, 22, 21-28.
Tripp, A.C., Hohmann, G.W., & Swift Jr., C.M. (1984). Two-dimensional resistivity inversion. Geophysics, 49, 1708-1717.
Tso, C.-H.M., Kuras, O., Wilkinson, P.B., Uhlemann, S., Chambers, J.E., Meldrum, P.I., Graham, J., Sherlock, E.F., & Binley, A. (2017). Improved characterisation andmodelling of measurement errors in electrical resistivity tomography (ERT) surveys. Journal of Applied Geophysics, 146, 103– 119.
Wilkinson, P., Chambers, J., Uhlemann, S., Meldrum, P., Smith, A., Dixon, N., & Loke, M.H. (2016). Reconstruction of landslide movements by inversion of 4-D electrical resistivity tomography monitoring data. Geophys. Res. Lett. 43, 1166–1174
Zhdanov, M., Endo, M., Cox, L., & Sunwall, D. (2018). Effective-medium inversion of induced polarization data for mineral exploration and mineral discrimination: Case study for the copper deposit in Mongolia. Minerals, 8(2), 68.
Zhou, B., & Dahlin, T. (2003). Properties and effects of measurement errors on. Near Surf. Geophysics, 1(3), 105-117.
Journal of the Earth and Space Physics
Volume 49, Issue 1
online publication date: 14 June 2023
June 2023
Pages 75-95

Files

Share

How to cite

Statistics