考虑电流型畸变的大地电磁三维反演

大地电磁测深(MT)的观测数据易受到由近地表小尺度非均匀体或地形起伏引起的电流型畸变干扰,消除或压制这种干扰对获取可靠的深部电性结构至关重要.当区域结构为二维时,电流型畸变可采用张量分解等方法予以消除或压制.当区域结构为三维时,畸变问题更加复杂和严重,传统张量分解方法往往效果不佳或无效,严重地制约了MT三维反演技术的实用性.对此,本文提出一种考虑电流型畸变的MT三维反演算法,将完整的电流型畸变参数引入到目标函数,并采用非线性共轭梯度法与电阻率参数同时反演,从而达到压制畸变的目的.该算法有两个关键点:一是通过分析实测数据所遭受畸变的分布特征,在目标函数中对其进行有效约束;二是在迭代过程中,通过自适应地调整双正则化因子保障算法的稳定和效率.理论模型测试结果显示,常规三维反演算法不能合理解释数据中的畸变成分,而只能通过引入虚假异常体强制地拟合受畸变数据,从而造成电阻率模型严重失真.与之相比,本文算法能够在反演中自动求解各测点所受到的畸变,获得更接近真实的电阻率模型.

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

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

阻抗张量分解进行大地电磁静校正的研究

地表电性不均匀体对大地电磁测深的静态影响，会使大地电磁测深曲线产生十分严重的畸变，如何进行静校正以获得合理的区域构造响应，引起了国内外大地电磁专家的广泛关注。阻抗张量分解就是近几年发展起来的一种有效的方法。在全面论述阻抗张量分解进行大地电磁静校正原理的基础上，通过模型试算证明了该方法的有效性。