殷长春,,
任秀艳,
刘云鹤,
惠哲剑,
谷宇
吉林大学地球探测科学与技术学院, 长春 130026
基金项目: 北京市科技计划"地球深部探测技术攻关"专项（Z181100005718001），国家自然科学基金重点项目（42030806），国家自然科学基金面上项目（41774125，42174080）及广西自然科学基金创新研究团队项目（GXNSFGA380004）联合资助
详细信息
作者简介: 朱姣, 女, 1992年生, 博士, 主要从事电磁勘探正反演方法技术研究.E-mail: zhujiaojlu@163.com
通讯作者: 殷长春, 男, 1965年生, 教授, 国家特聘专家, 主要从事电磁勘探理论, 特别是航空和海洋电磁方面的研究.E-mail: yinchangchun@jlu.edu.cn
中图分类号: P631收稿日期:2021-03-29
修回日期:2021-06-08
上线日期:2021-12-10
DC resistivity forward modelling for arbitrarily anisotropic media using the unstructured spectral element method
ZHU Jiao,YIN ChangChun,,
REN XiuYan,
LIU YunHe,
HUI ZheJian,
GU Yu
College of Geoexploration Science and Technology, Jilin University, Changchun 130026, China
More Information
Corresponding author: YIN ChangChun,E-mail:yinchangchun@jlu.edu.cn
MSC: P631--> Received Date: 29 March 2021
Revised Date: 08 June 2021
Available Online: 10 December 2021
摘要
摘要:各向异性是地电异常解释中不可忽视的因素，广泛存在于裂隙或层理发育的地质环境中.本文针对任意各向异性条件下直流电阻率法三维正演问题进行研究，结合非结构谱元法建立模拟算法，充分利用谱方法的指数收敛性以及非结构有限元对地形和复杂异常体刻画能力，提高计算精度和效率.通过灵活的四面体网格剖分和高阶谱插值，实现了复杂介质任意各向异性模型电阻率响应的高精度数值模拟.我们首先通过层状各向异性模型验证本文非结构谱元法的计算精度，进而我们以半空间中立方体模型为例分析各向异性对电阻率响应的影响特征，并通过计算针对不同各向异性参数的视电阻率极性图，探究地下介质各向异性特征识别方法.最后，我们针对典型的山脊模型计算和分析存在地形效应条件下各向异性电流场分布及视电阻率特征.模型计算结果表明基于四面体网格的谱元法模拟带有复杂地形和异常体的任意各向异性模型具有很高的计算精度.本文的研究成果将在推进电阻率方法用于解决裂隙及层理等环境和工程地质问题中发挥积极作用.
关键词: 谱元法/
非结构网格/
任意各向异性/
直流电阻率法/
正演模拟
Abstract:Anisotropy widely exists in the geological environment with well-developed fractures or bedding, which is a key factor that cannot be ignored when interpreting DC resistivity data. In this paper, we study three-dimensional forward modelling of the DC resistivity method for arbitrarily anisotropic media, and establish a simulation algorithm based on the unstructured spectral element method. This method makes full use of exponential convergence of the spectral method and the ability of the unstructured finite-element to delineate terrain and complex underground structures, so that the modelling accuracy and efficiency can be improved simultaneously. Via flexible tetrahedral meshing and high-order spectral interpolation, a high-precision numerical simulation of resistivity response for arbitrarily anisotropic media can be achieved. To do so, we first verify the accuracy of the unstructured spectral element method on a layered anisotropic model, and then take the model of a cube buried in a half space as an example to analyze the influence of anisotropy on resistivity response. We compare the polar plots of apparent resistivity for different anisotropy parameters to explore the means of identifying the anisotropy of underground media. Finally, we calculate and analyze the distribution of the current and apparent resistivity affected by the topography in a typical anisotropic ridge model. Calculation results show that the spectral element method based on tetrahedral grids has high calculation accuracy for simulating arbitrary anisotropic models with complex terrain and abnormal bodies. Research results of this paper would contribute to the advancement of the DC resistivity method applied to the solution of the environmental and engineering problems.
Key words:Spectral element method/
Unstructured grid/
Arbitrarily anisotropic medium/
DC resistivity method/
Forward modelling
PDF全文下载地址:
http://www.geophy.cn/data/article/export-pdf?id=16d0e40b-3543-4975-90f7-6dbc754c7219
