Efficient 2D inversion of long ERT sections

被引:9
|
作者
Tsourlos, Panagiotis [1 ]
Papadopoulos, Nikos [2 ]
Papazachos, Costas [1 ]
Yi, Myeong-Jong [3 ]
Kim, Jung-Ho [3 ]
机构
[1] Aristotle Univ Thessaloniki, Dept Geophys, Sch Geol, Thessaloniki 54124, Greece
[2] Fdn Res & Technol Hellas FORTH, Inst Mediterranean Studies, Lab Geophys Satellite Remote Sensing & Archaeoenv, Rethimnon 74100, Crete, Greece
[3] Korean Inst Geosci & Mineral Resources KIGAM, Mineral Resources Res Div, Dept Explorat Geophys & Min Engn, Taejon, South Korea
关键词
Fast 2D inversion; Jacobian matrix; Least-squares; LSQR; RESISTIVITY DATA; DC RESISTIVITY; LEAST-SQUARES; ALGORITHM;
D O I
10.1016/j.jappgeo.2014.03.022
中图分类号
P [天文学、地球科学];
学科分类号
07 ;
摘要
In this work a new algorithm for the efficient and fast two dimensional (2D) inversion of long electrical resistivity tomography (ERT) sections is introduced. The algorithm is based on 2.5D finite element method (FEM) scheme to solve Poisson's equation that describes the current flow into the earth's subsurface. The adjoint equation technique was incorporated into the FEM framework to estimate the sensitivity values. The reconstructed 20 resistivity models are recovered through an iterative, non-linear smoothness constrained least-squares approach. The algorithm further incorporates an experimental procedure to avoid the calculation and storage of the entire Jacobian matrix. The basic concept of this new algorithm relies on the fact that for every measurement there is a number of model parameters which are located in parts of the 2-D model at distant locations from potential and current electrodes. The corresponding absolute Jacobian matrix values in such cases are very small (almost zero) and can be omitted by the Jacobian calculation. Around every measurement a fixed rectangular threshold region is defined a-priori based on geometrical criteria. The algorithm calculates only Jacobian matrix values for the model parameters that are included in this threshold area omitting the calculation of the Jacobian entries related to model parameters outside this region. This approach speeds up the Jacobian matrix calculations while the efficient storage of the sparse Jacobian and Smoothness matrices and the inversion using an iterative routine like LSQR method increase significantly the inversion speed and reduce the memory requirements. The new algorithm is almost more than one order of magnitude (similar to 30 times) faster and consumes one order of magnitude (similar to 90%) less storage memory than the original one based on full Jacobian calculations for typical applications. The application of the new algorithm to synthetic and real data sets shows that the reconstructed models exhibit comparable accuracy to the standard inversion approach. (C) 2014 Elsevier B.V. All rights reserved.
引用
收藏
页码:213 / 224
页数:12
相关论文
共 50 条
  • [41] Estimating free parameters in 2D inversion: example of gravity inversion in a rifted basement
    Silva, Joao B. C.
    Santos, Darciclea F.
    Monteiro, Daniele P.
    GEOPHYSICAL PROSPECTING, 2020, 68 (03) : 1057 - 1071
  • [42] 2D or not 2D? Adaptive 3D Convolution Selection for Efficient Video Recognition
    Li, Hengduo
    Wu, Zuxuan
    Shrivastava, Abhinav
    Davis, Larry S.
    2021 IEEE/CVF CONFERENCE ON COMPUTER VISION AND PATTERN RECOGNITION, CVPR 2021, 2021, : 6151 - 6160
  • [43] 2D resistivity inversion of perturbation in the layered-earth
    Yan, YL
    Ma, XB
    Di, QY
    Kong, XR
    Gao, LB
    Chen, G
    CHINESE JOURNAL OF GEOPHYSICS-CHINESE EDITION, 2004, 47 (06): : 1139 - 1144
  • [44] 2D inversion processing of geoelectric measurements with archaeological aim
    E. Turai
    L. Hursán
    Acta Geodaetica et Geophysica Hungarica, 2012, 47 : 245 - 255
  • [45] Real time 2D inversion of induction logging data
    Tabarovsky, LA
    Rabinovich, MB
    JOURNAL OF APPLIED GEOPHYSICS, 1998, 38 (04) : 251 - 275
  • [46] LANDSLIDE STUDY WITH 2D ELECTRICAL RESISTIVITY TOMOGRAPHY (ERT): A CASE STUDY FROM TURKEY
    Alpaslan, Nuray
    Bayram, Mehmet
    CARPATHIAN JOURNAL OF EARTH AND ENVIRONMENTAL SCIENCES, 2020, 15 (02): : 391 - 403
  • [47] 2D gravity inversion technique in the study of Cheshire basin
    1600, E-Journal of Geotechnical Engineering (19 S):
  • [48] 2D APPROACH TO QUANTITATION OF INVERSION-RECOVERY DATA
    VANORMONDT, D
    DEBEER, R
    MARIEN, AJH
    DENHOLLANDER, JA
    LUYTEN, PR
    VERMEULEN, JWAH
    JOURNAL OF MAGNETIC RESONANCE, 1990, 88 (03): : 652 - 659
  • [49] Fault detection by 2D resistivity inversion on a topographic area
    Yang, CH
    Chen, PH
    Chen, JJ
    TERRESTRIAL ATMOSPHERIC AND OCEANIC SCIENCES, 2000, 11 (04): : 909 - 928
  • [50] Layered and laterally constrained 2D inversion of resistivity data
    Auken, Esben
    Christiansen, Anders Vest
    Geophysics, 3 (752-761):