三维地形大地电磁数值模拟方法
摘要
大地电磁场地形影响数值模拟国内外已有不少研究.如ChouteauandBouchard(1988)和Wannamaker,StodtandRiio(1986)用有限元法模拟了二维地形对大地电磁测深信号的影响;XuS.Z.andZhouH.(1977)用边界元法模拟了二维地形影响;1997年,徐世浙、阮百尧和周辉等人又用边界元法模拟了大地电磁测深的三维地形影响.本文在徐世浙、阮百尧和周辉等人研究基础上,对原来需要在地上空气介质和地下导电介质2个区域的边界面积分通过迭加方法变为只对地形界面的面积分,对原来只考虑磁场方程变为同时考虑磁场方程和电场方程,对原来三角单元电磁场使用零次插值改为仅地形影响项采用零次插值.由于这些改进,不仅使方法的计算精度得到了提高,而且使最后形成的线性方程组具有明显的对角占优特征,可用SSOR方法求解,因此与原方法相比大大提高了计算速度。