大地电磁相位张量分解特性.docVIP

精品论文 参考文献 大地电磁相位张量分解特性 陈江平CHEN Jiang-ping曰周亚东ZHOU Ya-dong（成都理工大学地球物理学院，成都610059）（College of Geophysics，Chengdu University of Technology，Chengdu 610059，China） 　　摘要院在大地电磁资料的分析中，其中相当大的一部分问题都是因为近地表电性结构水平非均匀性导致二维或三维的区域构造中水平均匀的正常电磁场发生畸变，从而对MT 阻抗张量产生影响[1]。故研究MT 阻抗张量畸变和分解方法对于反演解释的正确性有重要的意义。 　　Abstract: In magnetotelluric data analysis, a considerable part of the problems are due to the distortion of even level normalelectromagnetic field in the 2D or 3D regional structure led by the uneven level of near-surface electrical structure, thus influencing the MTimpedance tensor. Therefore, Researching the MT impedance tensor distortion and the decomposition methods is important for thecorrectness of the inverse interpretation.关键词院大地电磁；阻抗张量；局部畸变Key words: magnetotelluric；impedance tensor；local distortion中图分类号院P318 文献标识码院A 文章编号院1006-4311（2014）31-0314-02 　　0 引言Swift 分解是将大地电磁张量Z 在水平面内旋转，使之与区域构造主轴达到最佳匹配的方法。严格的说，Swift分解法适合用于标准的二维构造。对观测阻抗张量做该分解后，产生区域构造的两个主轴阻抗和构造的主轴方位角，即从观测阻抗张量四个复数元素（八个标量参数）进行旋转分析[1]。通常情况下，实际构造并不是标准的二维构造，故观测阻抗张量将含局部畸变的影响，分解结果并不是很精确。阻抗相位张量旋转可以消除观测阻抗张量电场局部畸变的影响，而且不需要关于区域结构维性的假设。 　　对转换后的观测阻抗采用已有的Swift 旋转优化方法即可求得区域主轴方位角[2]。 　　1 阻抗相位张量分解原理地表的测量坐标系xyz 中，非均匀介质的大地电磁阻抗张量为：Zxx，Zxy，Zyx，Zyy。 　　分离复杂的阻抗张量得到他们的实部（X） 和虚部（Y），我们可以写为：Z=X+Yi （1）因此，我们定义阻抗相位张量的关系：渍=X-1Y （2）大地电磁的观测阻抗可以表示为：Z= Zxx Zxy Zyx Zyy 蓘蓡（3）因此我们可以得到阻抗相位张量为：渍= 渍xx 渍xy 渍yx 渍yy 蓘蓡（4）只要把实际观测的阻抗张量（3）按照（2）式进行相位变换得到（4），在新的转换基础上，按照与经典的Swift 完全相似的办法，求得阻抗张量的一系列重要参数[3]。 　　2 确定主轴的最佳方位对于转换的相位阻抗张量，根据Swift 旋转公式可获得平面上任意坐标系的转换阻抗。现使得旋转后目标函数：渍xx 2+ 渍yy 2寅Max；或者渍xy 2+ 渍yx 2寅min；则可获得最佳主轴方位角为：兹=14arctan（渍xy+渍yx）*（渍xx-渍yy）+（渍xx-渍yy）*（渍xy+渍yx）（渍xx-渍yy）*（渍xx-渍yy）-（渍xy+渍yx）*（渍xy+渍yx）由上式可得到4 个解，按照与Swift 旋转相似的方式，通过比较目标函数值的大小来确定正确的结果。 　　3 理论模型分析为说明阻抗相位张量方法不受电场局部畸变的影响，我们使用典型的二维模型的阻抗数据进行检验。 　　模型参数：走向沿正北方向的二维区域结构模型，Zr=蓘 -190-8i 180+7i 蓡 （5）再把（5）中的二维模型加入电场畸变，畸变矩阵为：C=蓘 csoinsxx -csoisnyy 蓡 渍0x 渍0y 蓘 蓡=蓘 00..197845 -00.9.26569 蓡x=10毅，y=15毅，渍x=渍y=1；则Z=CZr=蓘 -41.89.2315+4-2.70.7722i8i 13.73.1723++16..281958ii 蓡 （6）我们把加入电场畸变后的二维模型进行Swift