摘 要: | 
结构化六面体网格对地形起伏的模拟效果不好, 并且会对网格外延区域进行不必要的加密.有限元法在对区域剖分的网格选择上比较灵活, 可以采用非结构化网格对研究区域进行划分, 更适合对起伏地形和复杂的地质构造进行模拟.因此, 本文基于矢量有限元算法, 采用非结构化变形六面体网格来实现大地电磁三维正演, 既能较好的模拟地形起伏, 也能对局部重点区域进行加密, 减少外扩区域的网格划分数量, 减轻计算量. 通过对理论模型的正演试算, 验证了程序的正确性.对于反演, 采用高斯-牛顿算法, 分别对平地、双峰和双谷地形的理论模型进行了反演试算, 对比了带地形反演和不带地形反演的结果差异, 最后对日本Kayabe地区的阻抗张量数据进行了反演.理论模型和实测数据的反演结果表明本文算法能够用于三维大地电磁反演.

|