首页 > 文献资料
-
复杂二维/三维大地电磁的有限单元法正演模拟策略
复杂二维和三维大地电磁模型的正演数值模拟具有一定的挑战性.对于复杂的二维和三维大地电磁正演问题,我们采用有限单元法进行求解.有限单元法后形成一个线性方程组,系数矩阵足大型稀疏的带状对称复系数矩阵,并且其条件数远大于1,为严重病态矩阵,求解其对应方程组会遇到很多困难.不完全LU分解处理的Bi-CGSTAB迭代方法可用于该线性方程组的求解,并且具有速度快、精度高和稳定性好等优点;为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分;在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间.通过对二维和三维模型电磁响应的讣算,验证了算法的正确性.
关键词： 大地电磁 有限单元法 正演模拟 Bi-CGSTAB算法 -
大地电磁局部畸变异常特征三维数值模拟研究
局部畸变问题曾经困扰大地电磁资料反演解释几十年,大地电磁三维数值模拟技术的发展为剖析局部畸变特点和得到可靠的反演成像结果提供了技术基础.本文采用三维数值模拟成像方法对典型三维局部畸变模型进行模拟分析.三维数值模拟结果显示:电场分量垂直电性分界面的极化模式视电阻率曲线(对应二维情况下TM模式)在穿越低阻异常体界面时,曲线会先上移后下移,而在穿越高阻异常体界面时,曲线会先下移后上移,这与电性分界面处积累面电荷产生的二次电场有关.三维模型中XY模式、YX模式视电阻率和相位在三维异常体附近的水平变化是呈现近似垂向对称的,该现象与电场垂直跨越电性界面时视电阻率的变化规律是吻合的,当测线分别沿X方向和Y方向展布时,三维情况下的XY和YX模式分别对应二维情况下的TM模式.低阻小异常体对区域构造响应的畸变影响比高阻小异常体要严重.低阻小异常体对二维区域响应的两种极化模式视电阻率和相位都有非常明显的畸变影响,相比较而言对TE模式的畸变要大于TM模式,因此我们在做二维反演解释时,可优先考虑拟合TM模式数据.位于小异常体中心上方测点的三维畸变响应虽然与对应真实二维区域响应的差异比较大,但可以等效于某种二维模型响应,这种由局部畸变造成的假二维响应在实际野外数据的解释中是需要注意的.
