3D forward modeling and analysis of the loop-source transient electromagnetic method based on the finite-volume method for an arbitrarily anisotropic medium

被引:14
作者
Liu YaJun [1 ,2 ]
Hu XiangYun [1 ,3 ]
Peng RongHua [1 ]
Yogeshwar, Pritam [2 ]
机构
[1] China Univ Geosci, Inst Geophys & Geomat, Wuhan 430074, Hubei, Peoples R China
[2] Univ Cologne, Inst Geophys & Meteorol, D-50923 Cologne, Germany
[3] Hubei Subsurface Multiscale Imaging Key Lab, Wuhan 430074, Hubei, Peoples R China
来源
CHINESE JOURNAL OF GEOPHYSICS-CHINESE EDITION | 2019年 / 62卷 / 05期
关键词
Large loop-source; Transient electromagnetic method; Anisotropy; Finite-volume method; 3D forward modeling; GROUNDED WIRE SOURCE; FIELD; INVERSION; SOLVER;
D O I
10.6038/cjg2019M0532
中图分类号
P3 [地球物理学]; P59 [地球化学];
学科分类号
0708 ; 070902 ;
摘要
Electrical anisotropy of strata has been long recognized by field and laboratory observations. However, nearly all of interpretations of transient electromagnetic (TEM) data is based on the assumption of electric isotropy of media, which can cause misleading data interpretation in regions with strong electrical anisotropy. To clarify the influence of electrical anisotropy on loop-source TEM responses, we present a three dimensional (3D) robust finite-volume (FV) algorithm for simulating TEM responses in an arbitrarily anisotropic medium through solving Helmholtz equations in the time domain. The time domain Maxwell equations are discretized using the mimetic finite-volume method (MFV) on a conventional staggered grid in the space domain and discretized in the time domain using the backward Euler method. To reduce time steps required by computation, the modeling time is first divided into several intervals, each of which has a constant time step size. Then, the resulting system matrix in each interval is factored and solved by a direct solver, which allows the factored matrix to be used to calculate the TEM responses for subsequent time steps in the same time interval. The accuracy of our algorithm is validated against quasi-analytic solutions of a 1-D layered anisotropic model. Finally, by numerical experiments for 3D models with different types of electrical anisotropy, we analyze the influences of electrical anisotropy on TEM responses. The results demonstrate that TEM responses are mainly affected by the horizontal conductivity. The effect of dipping anisotropy on TEM responses is much greater than horizontal anisotropy. Besides, the horizontal principle axis direction of electrical anisotropy could be inferred from the TEM signal.
引用
收藏
页码:1954 / 1968
页数:15
相关论文
共 47 条
  • [1] A fully asynchronous multifrontal solver using distributed dynamic scheduling
    Amestoy, PR
    Duff, IS
    L'Excellent, JY
    Koster, J
    [J]. SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS, 2001, 23 (01) : 15 - 41
  • [2] Transient electromagnetic surveys for the measurement of near-surface electrical anisotropy
    Dennis, Z. R.
    Cull, J. P.
    [J]. JOURNAL OF APPLIED GEOPHYSICS, 2012, 76 : 64 - 73
  • [3] A parallel finite-element time-domain method for transient electromagnetic simulation
    Fu, Haohuan
    Wang, Yingqiao
    Um, Evan Schankee
    Fang, Jiarui
    Wei, Tengpeng
    Huang, Xiaomeng
    Yang, Guangwen
    [J]. GEOPHYSICS, 2015, 80 (04) : E213 - E224
  • [4] 3D forward modeling of magnetotelluric fields in general anisotropic media and its numerical implementation in Julia
    Han, Bo
    Li, Yuguo
    Li, Gang
    [J]. GEOPHYSICS, 2018, 83 (04) : F29 - F40
  • [5] The electromagnetic response in a layered vertical transverse isotropic medium: A new look at an old problem
    Hunziker, Jurg
    Thorbecke, Jan
    Slob, Evert
    [J]. GEOPHYSICS, 2015, 80 (01) : F1 - F18
  • [6] Mimetic discretizations for Maxwell's equations
    Hyman, JM
    Shashkov, M
    [J]. JOURNAL OF COMPUTATIONAL PHYSICS, 1999, 151 (02) : 881 - 909
  • [7] Key K, 2012, GEOPHYSICS, V77, pF21, DOI [10.1190/geo2011-0237.1, 10.1190/GEO2011-0237.1]
  • [8] Joint inversion of CMP gathers of multi-channel transient electromagnetic data
    Li, Hai
    Xue, Guo-Qiang
    Zhong, Hua-Sen
    Di, Qing-Yun
    [J]. CHINESE JOURNAL OF GEOPHYSICS-CHINESE EDITION, 2016, 59 (12): : 4439 - 4447
  • [9] 3D numerical simulation for the transient electromagnetic field excited by the central loop based on the vector finite-element method
    Li, J. H.
    Zhu, Z. Q.
    Liu, S. C.
    Zeng, S. H.
    [J]. JOURNAL OF GEOPHYSICS AND ENGINEERING, 2011, 8 (04) : 560 - 567
  • [10] A vector finite element solver of three-dimensional modelling for a long grounded wire source based on total electric field
    Li Jian-Hui
    Farquharson, Colin G.
    Hu Xiang-Yun
    Zeng Si-Hong
    [J]. CHINESE JOURNAL OF GEOPHYSICS-CHINESE EDITION, 2016, 59 (04): : 1521 - 1534