大地电磁相位张量和倾子的三维主轴各向异性反演

地壳和上地幔普遍存在的电性各向异性是大地电磁测深方法在实践应用中的关键和难点问题,而由近地表局部地电异常引起的电流型畸变可进一步加剧各向异性反演的多解性和不确定性.为压制电流型畸变影响,本文实现了一种利用大地电磁相位张量和倾子数据获取地下介质主轴电性各向异性结构的三维反演算法.该算法正演问题的求解采用交错采样有限差分法,反演问题的求解采用非线性共轭梯度法.理论模型测试表明,相位张量和倾子的三维各向异性反演不仅恢复了真实的各向异性异常分布特征,且压制了为拟合电流型畸变出现的假异常.由于相位张量和倾子数据不含幅值信息,先验模型选取对反演模型电阻率值的恢复效果影响较大,但初始模型的改变基本不影响反演结果.为进一步验证该算法的有效性,将其应用于USArray项目在美国西北部采集的109个长周期大地电磁实测数据.反演结果表明,相对于阻抗数据,相位张量和倾子数据的各向异性反演结果对浅层结构的刻画更为准确,与地表地质特征更为吻合.

     

基于有限内存拟牛顿法的电阻率法三维主轴各向异性反演研究

基于各向同性假设的电阻率法三维正反演方法成熟并得到广泛应用.由于地下介质普遍存在电阻率各向异性,为提高反演精度,各向异性介质的三维正反演问题有待深入研究.本文将三维正演中总场分为一次场和二次场计算,二次场使用非结构网格有限单元法求解,基于有限内存拟牛顿法,系统研究了电阻率法三维主轴各向异性反演.通过设计低阻异常体与双异常体模型,利用地表、地井、井地与井间多种观测方式的合成数据进行了反演试算.反演结果表明：三维主轴各向异性反演能够可靠恢复异常体的位置与形状.与仅在地表采集数据的反演结果对比,地井、井地与井间多种观测数据的加入,反演结果能更好地恢复异常体的各向异性特征.证明了本文基于有限内存拟牛顿法电阻率三维反演方法的可行性与有效性.     

基于交叉梯度约束的重力、磁法和大地电磁三维联合反演

为了降低单一地球物理方法反演的多解性及受噪声的影响程度,本文围绕重力、磁法和大地电磁法开展了三维联合反演的研究.重、磁采用基于对数障碍法的正则化反演算法,大地电磁使用limited-memory Broyden-Fletcher-Goldfarb-Shanno（L-BFGS）反演算法,引入交叉梯度函数实现了三种物性结构的相互耦合,最终开发出一套重磁电三维联合反演算法,并实现MPI并行加速计算.通过理论模型算例验证了算法的准确性,结果表明：不论是单棱柱体模型还是组合棱柱体模型,联合反演结果相较单独反演对于异常体的空间形态刻画以及物性数值恢复具有较好的提升;单棱柱体模型算例使得异常体的物性参数（密度、磁化率和电阻率）更加接近于真实的物性参数;组合棱柱体模型的联合反演结果不仅仅消除了围岩物性参数的假异常,而且还增强了异常体边界结构的恢复程度.