激发极化法是矿产勘查、水文地质和环境监测中应用广泛的技术之一，尤其是在金属矿产勘查中不可或缺。现有工作表明，各向异性介质中的电磁场与各向同性介质场有较大不同，这导致在使用现有各向同性介质理论时，数据解释出现较大偏差甚至严重错误。迄今为止，以往文献中尚未发表各向异性极化介质的三维激电正演模拟。本文使用有限元法实现了具有任意各向异性的三维激电正演模拟。通过对一个经典的两层各向异性模型的数值解和解析解的比较，证明了该算法的有效性。通过非结构化网格在平坦地形条件下实现了各向异性介质的三维激电正演模拟，并研究了地形和各向异性对激电数据解释的影响，这是三维激电法领域的新发展，具有重要的理论意义。
The induced polarization method is one of the most widely used technologies in mineral exploration, hydrogeology, and environmental monitoring, especially in metal mineral exploration. The existing study shows that the electromagnetic field of anisotropic mediums is very different from that of isotropic mediums, which leads to large deviations or even serious errors in data interpretation whenusing the existing isotropic theory. Up to now, three-dimensional induced polarization forward modeling of anisotropic mediums has not been published in previous literature. In this paper, a three-dimensional induced polarization forward modeling with arbitrary anisotropy is established by using the finite element method. The effectiveness of the algorithm is proved by comparing the numerical and analytical solutions of a classical two-layer anisotropic model. Three-dimensional induced polarization forward modeling for arbitrary anisotropic mediums is achieved through unstructured grids under flat terrain conditions, and the influence of terrain and anisotropy on induced polarization data interpretation is studied. This is a new development in the field of the three-dimensional induced polarization method, which has important theoretical significance.
2023,44（3）: 40-48 收稿日期：2023-01-30
DOI：10.3969/j.issn.1003-3246.2023.03.006
基金项目：国家重点研发计划（项目编号：2018YFE0208300）
作者简介：胡代明(1992-),男,博士,工程师,主要从事电磁法勘探的相关研究。E-mail:hdm0206@mail.ustc.edu.cn
*通讯作者：刘昊,男,高级工程师,主要从事水利水电、地质工程管理工作。E-mail:330633026@qq.com
参考文献：
傅良魁. 电法勘探教程[M]. 北京:地质出版社,1983:50-60.
景强. 高密度电阻率法和激发极化法在青山水厂找水中的应用[J]. 煤矿现代化,2018,(1):89-91.
李帝铨,王光杰,底青云,等. 大功率激发极化法在额尔古纳成矿带中段找矿中的应用[J]. 地球物理学进展,2007,22(5):1 621-1 626.
李金铭. 地电场与电法勘探[M]. 北京:地质出版社,2005:25-35.
孙仁斌,楚丽霞,王宁,等. 时间域激电测深二维反演在内蒙古兴和盆地贫水区找水勘查中若干案例研究[J]. 地球物理学进展,2017,32(1):387-394.
向葵,胡文宝,严良俊,等. 页岩气储层特征及地球物理预测技术[J]. 特种油气藏,2016,23(2):5-8.
赵后越. 大功率激电在张家口地区石墨矿勘查中的应用[J]. 工程地球物理学报,2017,14(5):546-551.
赵荣春,吕玉增,凌嘉宣,等. 激电中梯和对称四极测深在广西某铅锌矿区的应用[J]. 工程地球物理学报,2016,13(3):271-276.
郑冰. 频谱激电法在某铅锌银矿的应用[J]. 工程地球物理学报,2015,12(6):750-754.
Boaga J. The use of FDEM in hydrogeophysics:A review[J]. Journal of Applied Geophysics, 2017, 139:36-46.
Ciarlet P G. The finite element method for elliptic problems[M]. Philadelphia:SIAM, 2002:301-311.
Dias C A. A non-grounded method for measuring electrical induced polarization and conductivity[M]. Berkeley:University of California, 1968:17-30.
Herwanger J V, Pain C C, Binley A, et al. Anisotropic resistivity tomography[J]. Geophys J Int, 2004, 158(2):409-425.
Li Y G, Luo M, Pei J X. Adaptive finite element modeling of marine controlled-source electromagnetic fields in two-dimensional general anisotropic media[J]. J Ocean Univ China, 2013, 12(1):1-5.
Linde N, Pedersen L B. Evidence of electrical anisotropy in limestone formations using the rmt technique[J]. Geophysics, 2004, 69(4):909-916.
Liu W Q, Lin P R, Lü Q T, et al. Time domain and frequency domain induced polarization modeling for three-dimensional anisotropic medium[J]. Journal of Environmental and Engineering Geophysics, 2017, 22(4):435-439.
Okay G, Cosenza P, Ghorbani A, et al. Localization and characterization of cracks in clay-rocks using frequency and time-domain induced polarization[J]. Geophysical Prospecting, 2013, 61(1):134-152.
Pekşen E, Yas T, Kıyak A. 1-D DC resistivity modeling and interpretation in anisotropic media using particle swarm optimization[J]. Pure and Applied Geophysics, 2014, 171(9):2 371-2 389.
Pelton W H, Ward S H, Hallof P G, et al. Mineral discrimination and removal of inductive coupling with multifrequency IP[J]. Geophysics, 1978, 43(3):588-609.
Revil A, Florsch N. Determination of permeability from spectral induced polarization in granular media[J]. Geophysical Journal International, 2010, 181(3):1 480-1 498.
Schmutz M, Albouy Y, Guérin R, et al. Joint electrical and time domain electromagnetism (TDEM) data inversion applied to the Super Sauze earthflow (France)[J]. Surveys in Geophysics, 2000, 21(4):371-390.
Schwartz N, Shalem T, Furman A. The effect of organic acid on the spectral-induced polarization response of soil[J]. Geophysical Journal International, 2014, 197(1):269-276.
Seigel H O. Mathematical formulation and type curves for induced polarization[J]. Geophysics, 1959, 24(3):547-565.
Soueid Ahmed A, Revil A. 3-D time-domain induced polarization tomography:a new approach based on a source current density formulation[J]. Geophysical Journal International, 2018, 213(1):244-260.
Yang X H, He J S, Tong X Z. Numerical simulation of frequency-domain IP with FEM[J]. Progress in Geophysics, 2008, 23(4):1 186-1 189.
