各向异性介质空间—波数混合域直流电阻率法三维数值模拟研究 |
| |
引用本文: | 戴世坤,凌嘉宣,陈轻蕊,张莹,李昆.各向异性介质空间—波数混合域直流电阻率法三维数值模拟研究[J].地球物理学报,2022(7):2729-2740. |
| |
作者姓名: | 戴世坤 凌嘉宣 陈轻蕊 张莹 李昆 |
| |
作者单位: | 1. 中南大学有色金属成矿预测与地质环境监测教育部重点实验室;2. 中南大学有色资源与地质灾害探查湖南省重点实验室;3. 中南大学地球科学与信息物理学院;4. 西南石油大学地球科学与技术学院 |
| |
摘 要: | 将各向异性介质分成各向同性背景介质和各向异性异常介质,并提出一种空间-波数混合域方法实现了各向异性介质下直流电阻率法的三维数值模拟.不同于传统直流电阻率数值模拟方法,本文算法直接对空间域异常电位满足的偏微分方程沿水平方向进行二维傅里叶变换,使水平方向转换成波数域,保留垂直方向为空间域,可根据地下介质电流密度变化的快慢灵活剖分.这样可把空间域异常电位满足的三维偏微分方程转化成不同波数满足的一维常微分方程,把一个大规模三维数值模拟问题分解为多个一维数值模拟问题,利用一维有限单元法求解方程组,并通过采用压缩算子迭代计算,最终获得较为精确的数值解.与自适应有限单元法对比验证了本文算法的正确性;测试了算法的收敛性,结果表明在满足精度要求的情况下,算法的收敛性只与异常体和围岩之间的电导率差异相关,而与异常体大小和埋深无关;分析了算法计算效率,结果表明算法的计算效率与剖分网格节点成线性关系,算法可在微型计算机中较快计算出剖分节点总数超过千万的各向异性模型的结果;设计简单观测系统并验证其具有反映地下各向异性结构特性的能力;最后模拟异常体沿着不同方向旋转不同角度时的响应特征,对比分析可知异常体为各向异性...
|
关 键 词: | 空间波数混合域 二维正反傅里叶变换 有限单元法 压缩算子 |
|