中国科学技术大学学报 ›› 2020, Vol. 50 ›› Issue (5): 705-714.DOI: 10.3969/j.issn.0253-2778.2020.05.019

• 论著 • 上一篇    

起伏地形下时间域激发极化法三维非结构有限元数值模拟

黄渡,吴小平   

  1. 1.中国科学技术大学地球和空间科学学院地震与地球内部物理实验室,安徽合肥 230026; 2.蒙城地球物理国家野外科学观测研究站,安徽蒙城 233500
  • 收稿日期:2019-08-26 修回日期:2020-05-08 接受日期:2020-05-08 出版日期:2020-05-31 发布日期:2020-05-08
  • 通讯作者: 吴小平
  • 作者简介:黄渡,男,1994年生,硕士.研究方向:直流电法.E-mail:huangdu@mail.ustc.edu.cn
  • 基金资助:
    国家自然科学基金(41874084,U2039206)和国家重点研发计划项目(2018YFE0208300,2018YFC0603500)资助.

3D time-domain IP forward modeling using unstructured finite element method under uneven terrain

HUANG Du, WU Xiaoping   

  1. 1. Laboratory of Seismology and Physics of Earth’s Interior, School of Earth and Space Sciences, University of Science and Technology of China, Hefei 230026, China; 2. National Geophysical Observatory at Mengcheng, Mengcheng 233500, China
  • Received:2019-08-26 Revised:2020-05-08 Accepted:2020-05-08 Online:2020-05-31 Published:2020-05-08

摘要: 一般认为,起伏地形对激发极化法勘探影响不大,这是因为在地下介质极化率均匀的理想情况下起伏地形没有影响.当存在矿体的复杂介质模型时,起伏地形对激发极化法预期会存在一定影响.因为起伏地形的模拟,特别是对于三维模型具有较大的挑战性,所以目前还缺乏地形影响的定量化研究.本文实现了时间域激发极化三维非结构有限元数值模拟,由于采用了非结构网格剖分,特别适合任意起伏地形三维地电模型的数值模拟及其对激发极化影响的研究.在此基础上本研究对起伏地形下的板状、球状异常体进行了三维激发极化数值模拟,探讨了起伏地形对激发极化响应的影响.研究结果表明:非对称地形会使激发极化法的激电异常中心位置发生偏离,斜坡坡度越大异常的中心位置偏离越大,而且异常幅值也减小;地形影响的幅度与极化异常体非对称性存在关联性.

关键词: 激发极化法, 起伏地形, 等效电阻率, 有限元, 非结构, 三维正演

Abstract: It is generally believed that the undulating terrain has little influence on the induced polarization (IP) exploration, because it has no influence under the relatively ideal condition of uniform polarization of the underground medium. However, for the complex medium model with ore body under the condition of undulating terrain, the terrain will certainly have an impact on the induced polarization method. Three have been very few studies on the quantitative characterization of terrain effect on the induced polarization method,mainly because the undulating terrain is not easy to simulate, especially for the three-dimensional model. In this paper, the three-dimensional unstructured finite element numerical simulation of induced polarization is developed. The unstructured mesh is particularly suitable for the numerical simulation of three-dimensional geoelectric model under arbitrary undulating terrain and the study of its influence on induced polarization. With the newly developed method, the three-dimensional induced polarization numerical simulation for plate-shaped and spherical anomaly bodies in rolling terrain is carried out, and the influence of rolling terrain on induced polarization response is discussed. The results show that the asymmetric topography will make the center position of IP anomaly deviate. The larger the slope gradient is, the larger the deviation will be, and the IP amplitude becomes smaller with increasing slope gradient. The asymmetry of polarization anomaly also plays an important role in the topography impact for its IP response.

Key words: induced polarization method, uneven terrain, equivalent resistivity, finite element method, unstructured, 3D forward modeling

中图分类号: