论文部分内容阅读
大地电磁测深法(Magnetotelluric Sounding,MT)是以地下介质的电性差异为基础,以天然电磁场为场源,研究地下电性结构的频率测深方法。本文主要研究考虑地球曲率的二维MT正反演算法,并应用到云南福贡-陆良大地电磁剖面的资料处理中。 二维介质情况下,长剖面在地球尺度下表现为一条弯曲的弧线,在进行大地电磁正反演的时候需要考虑地球曲率。本文首先推导了大地电磁传播的偏微分方程和二维变分方程,正演采用扇形剖分方法、双二次插值、高斯数值积分、预条件共轭梯度法解方程的有限单元法,通过和均匀半空间、H型地电模型的数值解对比验证程序的正确性,设计高、低阻异常体与平地面大地电磁方法进行对比,发现大地电磁对低阻体响应更好,且TE模式纵向分辨率高,TM模式横向分辨率高,考虑地球曲率的结果更合理。采用NLCG方法对高、低阻异常体进行反演,发现大地电磁反演对低阻异常体的反演效果更好,且异常体的横向分辨率高于纵向分辨率,纵向上异常体的顶界面分辨率高于底界面。剖面长度对小规模异常体的影响较小,对大规模异常体的影响较大,在长剖面情况下,异常体长度越长受地球曲率的影响就越大,其响应结果表现为与实际地面相符合的弧形状。 介绍了云南福贡-陆良实测大地电磁剖面的构造背景、研究进展及意义,从测线的布置、资料的采集和处理方法进行研究,重点研究了ARIMA模型在大地电磁缺失数据中的应用;分析Niblett-Bostick穿透深度可知研究区域TM模式深度大于TE模式,整体上平均深度80km;根据Swift分解计算电性主轴及二维偏离度,发现电性主轴方向与构造走向吻合,低频部分二维偏离度较高,三维性较强。反演得出福贡-陆良的初步电性结构,自西向东可划分以下电性结构类型:兰坪壳幔高-低-高层状电性结构、金沙江壳幔低阻带状电性结构、扬子陆块西南缘壳幔高阻块状电性结构、滇中壳幔高-低-中低阻层状电性结构和滇东壳幔高阻带状电性结构。金顶北衙矿集区下发现两个高阻体初步推断成因是超基性岩浆侵入到地壳浅部的残留体。大多数地震发生在电性结构分界或邻近低阻的高阻区域内,研究区域内可分为五个地震活跃的区域(金沙江断裂带地震区,宾川-程何断裂带地震区,姚安断裂带地震区,绿汁江断裂带地震区和小江断裂带地震区)。