耦合有限单元法扩边的直流电阻率勘探无单元Galerkin法正演

作者：麻昌英; 柳建新; 郭荣文; 孙娅; 崔益安; 刘嵘; 刘海飞

来源：地球物理学报, 2018, 61(06): 2578-2588.

• 无单元Galerkin-有限单元耦合法
• RPIM形函数
• 无单元Galerkin法
• 有限单元法
• 直流电阻率

本文分析了目前直流电阻率正演模拟中的无单元Galerkin法(EFGM)和有限单元法(FEM)的优缺点,针对采用第一类边界条件需要足够大的计算域时EFGM计算成本高的问题,在计算域外围区域采用FEM扩边,提出了直流电阻率的无单元Galerkin-有限单元耦合法(EFG-FE).采用具有Kronecker delta函数性质的径向基点插值法(RPIM)构造EFGM形函数,在外围区域将EFGM与FEM直接耦合,无需其他处理手段,消除了传统EFGM与FEM耦合中存在的界面耦合困难.EFG-FE将模型计算域分割为EFGM区域和FEM区域,模型核心区域采用EFGM计算,发挥EFGM灵活性、适应性强和高精度的优点,使得模型建立简单方便,对任意复杂地电模型适应性强,同时获得高精度模拟结果.在模型计算域外围采用快速扩展的FEM单元网格进行剖分,利用其数值稳定性和高效性,使用少量FEM节点和单元网格将计算域大范围扩大满足第一类边界条件,同时不大幅增加计算成本,进而提高计算效率.最后,通过不同正演方法的模型算例的模拟结果对比,验证了本文提出的EFG-FE有效可行,其模拟结果具有很高的模拟精度,且相比于采用第三类边界条件的EFGM提高了计算效率,具有更好的模拟性能.

• 单位
  中南大学; 有色资源与地质灾害探查湖南省重点实验室; 物理学院