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

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

直流电阻率法与大地电磁法的二维联合反演

本文对电性联合反演进行了深入研究,以减少地球物理反演的多解性.将直流电阻率(DC)与大地电磁(MT)数据加入到同一反演数据集中.引入Tikhonov正则化思想建立反演目标函数,使反演过程更加高效稳定.在解决正则化反演问题过程中,分别采用了二阶最大平滑稳定因子和改进的L-curve法,提高了反演结果的稳定性和正则化因子的求取精度;最后运用非线性共轭梯度法(NLCG)对反演目标函数实现最优化求解.经研究表明:联合反演方法与单一反演方法相比,能够更加有效的约束反演模型范围;反演算法快速稳定,提高了反演精度,减少了对地下地质结构认识的模糊性.