共 12 条
三维地形直流电阻率有限元法模拟
被引：52
作者：
强建科
[1
]
罗延钟
[2
]
机构：
[1] 中南大学信息物理工程学院
[2] 中国地质大学地球物理与空间信息学院
来源：
关键词：
电阻率法;
有限单元法;
三维数值模拟;
起伏地形;
D O I：
暂无
中图分类号：
P631.3 [电法勘探];
学科分类号：
摘要：
基于稳定电流场的基本方程、三维区域满足的边值问题以及相应的变分问题,研究了三维起伏地形条件下电阻率的有限单元数值模拟算法.离散积分区域时,以三棱柱为最小研究单元,推导了含有地形特征信息的三线性插值型函数以及单元刚度矩阵.采用变带宽、一维数组方式只存储稀疏刚度矩阵中非零元素,能够节约内存.利用Cholesky分解法只分解一次大型稀疏矩阵,通过回代可以求出方程组的全部解,当求解有多个供电点的测深问题时可以缩短计算时间.模型计算表明,在水平层状介质模型上,三维计算结果与解析解或二维数值解十分吻合,计算精度满足误差要求.在二维山脊上的二极剖面或三维山谷上的中间梯度剖面上,其三维计算结果与相应模型的土槽实验结果或边界元法计算结果也非常接近.
引用
收藏
页码：1606 / 1613
页数：8
相关论文