摘要

为了解决传统有限元截断边界所引起的问题,提出一种新的2.5D直流电阻率有限元-无限元耦合数值模拟方法。首先推导无限元2.5D单元映射函数,然后提出一种无限元形函数,将其与传统有限元相结合,取代传统的混合边界,使得电位在无限域内连续并在无限远处衰减为0 V,最终形成的刚度矩阵稀疏对称并与场源位置无关。研究结果表明:在相同的网格剖分下,有限元无限元耦合方法比传统有限元法能够在边界测点处得到更高的计算精度,能够在较小的计算范围内得到更优的计算结果,从而有利于减少节点数,提高计算速度;由于系数矩阵不随场源位置改变,有利于加速反演计算。