基于电场Helmholtz方程的回线源瞬变电磁法三维正演.pdfVIP

第56卷第12期 地 球 物 理 学 报 Voi．56，No．12 2013年12月 CHINESEJOURNALOFGEOPHYSICS Dec．，2013 I doi：10．6038／cj920131228． Li X S J H，Hu H，eta1．Three—dimensionalforwardcalculationfor sourcetransient method Y，Zeng loop electromagnetic based onelectricfieldHelmholtz equation．ChineseJ．Geophys．(in cj920131228． 基于电场Helmholtz方程的回线源 瞬变电磁法三维正演 李建慧1，胡祥云h，曾思红1，路金阁1，霍光谱2，韩波1，彭荣华1 1中国地质大学地球物理与空间信息学院，武汉430074 2河南省地质调查院，郑卅I 450001 摘要正演是电磁法勘探野外工作参数选取、室内资料处理与解释的基础，精确、稳定、高效的三维正演算法尤 为重要．本文采取先求解拉普拉斯域电场、再由Gaver—Stehfest算法获得时间域磁场的思路，基于电场异常场 Helmhohz方程实现了交错网格有限差分法和有限体积法对回线源瞬变电磁法的三维正演．通过对比低阻块状体 的积分方程法、时域有限差分法、矢量有限单元法和SLDM法的数值解，验证了交错网格有限差分法和有限体积法 的正确性．由于交错网格有限差分法、有限体积法和基于矩形块单元的矢量有限单元法将待求电场均定义在矩形 块单元棱边上，因此三种数值算法可采用相同方法进行电场待求量编码、计算背景场和后处理．然而，与矢量有限 单元法相比，交错网格有限差分法和有限体积法的系数矩阵更加稀疏，求解效率更高．通过对水平低阻板状体三维 模型的数值模拟，我们发现本研究中交错网格有限差分法比有限体积法精度更高；再利用一维解析法求解相应三 层层状地电模型的感应电动势，我们还发现两种数值算法和一维解析法计算的感应电动势等值线形状吻合程度 高，只是数值范围略有差异． 关键词Helmholtz方程，回线源，瞬变电磁法，有限差分法，有限体积法 doi：10．6038／cj920131228中图分类号P631 收稿13期2013-05—20，2013—08—27收修定稿 Three—dimensionalforwardcalculationfor sourcetransient loop methodbasedonelectricfieldHelmholtz equation LIJian—Huil，HU Jin—Gel， Xiang—Yunh，ZENGSi—Hon91，LU HUO B01，PENG Guang—Pu2，HANRong—Hual 1 Institute andGeometics，ChinaGeoscience，Wuhan430074，China ofGeophysics Universityof 2 HenanInstitute 450001，China