含地形起伏的MT模型非规则四边形FEM正演与反演

从大地电磁(MT)边值问题满足的变分方程出发,采用非规则四边形网格、双线性插值有限单元法(FEM)开展复杂起伏地形MT模型正演,探讨二维Jacobian变换行列式的计算方法,推导任意非规则四边形单元的插值方式及单元系数矩阵表达式,实现起伏地形MT模型的高精度正演。然后,介绍光滑约束的Tikhonov正则化反演算法,针对反演中正则化参数选取困难的问题,将广泛应用的L-curve法引入反演的正则化参数选取中。研究结果表明:L-curve法的曲线中曲率最大的拐点准确地指示了最优正则化参数;L-curve法选取的最优正则化参数对应的反演结果与原模型所示结果吻合度最高,反演效果最好。