大地电磁有限差分数值解对比

摘要：

基于网格精度、形成系统方程的方式、边界条件以及预条件线性算子等,文中对大地电磁(MT)有限差分数值解作了对比.对不同网格剖分方式下的三个均匀半空间模型的一维MT响应对比显示,在降低首层厚度的同时保持层间厚度变化在合理范围可以同时提高主场和辅助场的精度.在利用正常中心网格法(主场和辅助场都定义在单元顶面的中心)计算二维(2-D)TM模式响应时,应该从Maxwell一次差分方程开始组建二次差分方程,这样可以更充分考虑模型电阻率的变化.在对边界值如何影响数值解的测试表明,仅仅提高一维(1-D)边界值的精度对提高2-D MT有限差分数值解的精度是有限的.线性算子对提高MT解的效率十分重要,简单的对比进一步表明合适的预条件再配合好的线性算子(如文中求解2-D MT时所采用的DILU-BICGSTAB方法)不仅可以加速收敛,而且可以降低迭代次数.

关键词：大地电磁测深、有限差分、数值解、对比

所属期刊栏目：53

分类号：P631

资助基金：This study has been supported by the Director Foundation of Institute of Geology,China Earth-quake AdministrationGrant DF-IGCEA-0608-2-16;the National Natural Science Foundation of ChinaGrant 40534023 and Grant 40974042

在线出版日期：2010-05-19（万方平台首次上网日期，不代表论文的发表时间）

页数：共9页

页码：622-630