摘要

从大地电磁(MT)二维正演所满足的偏微分方程边值问题出发,应用矩形网格剖分和单元内双二次插值推导有限单元法求解大地电磁TE与TM两种极化模式正问题详细算法。应用反演理论将病态问题求解的正则化方法应用到最小二乘优化方法中,获得最光滑约束最小二乘正则化反演目标函数,并利用Matlab编制了大地电磁二维正反演计算程序。应用该程序对高低阻地电模型和Sasaki模型开展了正反演计算,并绘制TE模式和TM模式、TE&TM联合反演模式的反演成果剖面图。将所得的反演剖面与初始模型对比可知,TE模式反演剖面纵向分辨率较高,TM模式横向分辨率较高,TE&TM联合反演优于单一极化模式的反演,并证明双二次插值有限元法MT正演及最小二乘正则化反演算法的有效性与可行性。

全文