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

详细信息   

• 责任者：肖骑彬^c中国地震局地质研究所 ; 赵国泽
• 刊名：地球物理学报
• 出版年：2010
• 卷期：53^b3
• 关键词：大地电磁测深 ; 有限差分法20^a大地电磁法 ; 电法勘探 ; 一维模型 ; 电阻率 ; 数值模拟 ; 边界元法
• 页码：622～630
• 索书号：P260.66/141
• 图表：^c参30插图
• 分类号：2100
• 文摘期号：2010/09

摘要

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