山西阳高～河北容城剖面大地电磁资料的初步解释──阻抗张量分解技术及其应用

介绍了阻抗张量分解技术及其在阳高～容城剖面的大地电磁资料解释中的应用。在某些测点局部电流畸变的影响较严重，用张量分解技术对资料分析后表明，沿剖面的区域异常总体上具二维性质，分解技术处理后的视电阻率和相位资料与地质构造有较好的对应性。事实表明，阻抗张量分解技术是对大地电磁资料进行分析解释的有力工具     

表面局部三维大地电磁曲线畸变校正——MT畸变校正阻抗张量分解方法

针对三维／二维地质结构开发了大地电磁阻抗张量分解技术,利用该技术仅需对各个独立测点的大地电磁资料进行分析即可分离和消除由三维局部异常体引起的其附近测点的ＭＴ响应函数的形态扭曲畸变。     

大地电磁阻抗张量正则分解的地质含义分析

由于常规的张量旋转分析方法的局限性,在处理三维阻抗问题时不能得到令人满意的结果.最新发展的三维阻抗张量的正则分解方法由于不需约束条件,可以自然的从三维阻抗张量中提取完整的参数体系,因而具有很好的应用前景.本文对正则分解的方法从理论上进行了分析和证明,对原有参数计算公式做适当的修改,将正则分解与坐标旋转结合起来,分析二者之间的联系与区别,并探讨参数的地质含义.     

一种简单的大地电磁阻抗张量畸变分解方法

提出了一种简单的大地电磁阻抗张量畸变分解方法，引入3d/2d模型，使得全畸变的大地电磁阻抗张量被生意人经，通过对比Swift旋转和Bahr分解后的资料，用联系二的畸变矩阵表示浅层三维电性异常体对区域二维构造的局部畸变影响，通过研究畸变张量矩阵，并对其进行分解，可将局部三维电性不均匀体产生的畸变与区域构造分离开来，从而了解区域构造的维数特征和局部畸变类型及其对区域构造的影响，最后，应用3个测点的实际资料，对这种简单的畸变分解方法进行了检验，结果表明，该方法简单且应用方便，解释效果较好。     

大地电磁阻抗张量正则分解的地质含义分析

由于常规的张量旋转分析方法的局限性，在处理三维阻抗问题时不能得到令人满意的结果.最新发展的三维阻抗张量的正则分解方法由于不需约束条件，可以自然的从三维阻抗张量中提取完整的参数体系，因而具有很好的应用前景.本文对正则分解的方法从理论上进行了分析和证明，对原有参数计算公式做适当的修改，将正则分解与坐标旋转结合起来，分析二者之间的联系与区别，并探讨参数的地质含义.     

相位张量分析约束下的大地电磁测深阻抗张量分解方法研究

由于浅部、局部不均匀三维异常体的影响,大地电磁测深(Magnetotelluric,MT)观测数据往往需要进行畸变校正.本文针对相位张量不受局部电场畸变影响的特点,研究相位张量分析约束下的大地电磁测深(MT)阻抗张量GB(Groom-Bailey)分解方法.分析由相位张量得到的相关参数,结合初步电性结构维性判断及二维或近似二维电性结构条件下电性主轴方向等信息,设计较合理的GB阻抗张量分解初始模型,提高GB算法的稳定性及可靠性.同时,在传统GB分解求解模型的基础上,将各向异性参数引入目标函数,并采用共轭梯度法进行模型的最优化求解.理论数据计算结果表明,MT阻抗张量GB分解方法的计算结果,依赖于初始模型中电性主轴角、阻抗张量值等参数的选择,相位张量分析约束下的GB分解方法能有效改善算法的稳定性和计算结果的可靠性.此外,研究表明,在合理的初始模型选择下,较之传统的GB分解方法,尽管各向异性参数仍旧无法求出,但带各向异性参数的GB分解方法在计算效率和结果可靠性方面均具有一定的优势.选择实际大地电磁测深点数据进行方法验证,取得了理想的电性分析结果.     

大地电磁GB张量分解法的改进

在对Groom-Baily的大地电磁(MT)张量分解法(简称GB分解法)修改之后,采用对单个频点进行张量分解的方法,并将参加分解的阻抗张量元素改用视电阻率和相位,使得分解方程数目及求解未知数数目减少,大大改善了方程的病态性能并减少计算时间,从而使张量分解更少地受分解条件的限制而稳定地进行.在此基础上,编写了单频点张量分解程序,并应用于实际资料.     

山西阳高──河北容城剖面大地电磁资料的二维反演解释

介绍了一种快速二维大地电磁自动反演方法及其在阳高－容城剖面大地电磁资料解释中的应用．二维反演结果显示，研究区浅层电性结构与地表地质有较好的对应性；在中、下地壳范围内发育着断断续续的低阻带，并沿剖面大致分为三段；上地幔低阻带由北西向南东方向逐渐变浅，上地幔低阻带深度陡变带与巨型重力异常高梯度带对应．与一此较大型的断裂带对应的区段在电性上有较清楚的显示．     

山西阳高──河北容城剖面大地电磁资料的二维反演解释

介绍了一种快速二维大地电磁自动反演方法及其在阳高－容城剖面大地电磁资料解释中的应用．二维反演结果显示,研究区浅层电性结构与地表地质有较好的对应性;在中、下地壳范围内发育着断断续续的低阻带,并沿剖面大致分为三段;上地幔低阻带由北西向南东方向逐渐变浅,上地幔低阻带深度陡变带与巨型重力异常高梯度带对应．与一此较大型的断裂带对应的区段在电性上有较清楚的显示．     

大地电磁(MT)阻抗张量的正则分解及其初步应用

鉴于大地电磁阻抗张量传统分析方法的局限性,提出了普遍适用的正则分解方法。在用正则分解方法分析时,主坐标系提供的８个标量参数,描述了地球系统传输特性和场源极化状态及其联系,因此是物理意义明确、易于解释的参数体系。讨论了正则分解方法在各种构造条件下的特点和场源的极化问题,并给出了简单的应用实例     

基于遗传算法的大地电磁阻抗张量分解方法研究

本文针对Groom-Bailey分解法在消除地表电性不均匀体对MT数据的畸变效应时存在的问题引入遗传算法,对基于传统线性最优化方法的GB分解算法进行改进,提出基于遗传算法的大地电磁阻抗张量分解方法.通过对理论合成数据以及三维/二维模型正演数据的分解试验,并且对青藏高原北缘阿尔金断裂地区实际MT数据进行分解处理,证明了基于遗传算法的GB分解能够更加有效地校正三维/二维情况下近地表三维电性不均匀体所造成的畸变影响.最后,在已有算法基础上,研究了基于遗传算法的多频率GB分解算法和MT数据静校正算法,并通过实际MT数据的处理证明了这些算法的有效性.     

大地电磁阻抗张量不变量及其Mohr圆分析

论述了Ｍｏｈｒ圆及其在大地电磁阻抗张量分析中的应用。并指出，大地电磁阻抗张量不变量由于提供信息的多样性和稳定性，在应用中具有独特的优点，而Ｍｏｈｒ圆又将各种大地电磁阻抗张量不变量有机地联系在一起，使其成为分析研究大地电磁阻抗张量性质的重要图示工具，能直观、清晰、简便地给出构造维数、二维偏离度、各向异性等重要信息，可用来研究测区电性结构沿纵向和横向的变化特征     

川西—藏东大地电磁资料的阻抗张量畸变分解

根据Swift和Bahr的分解结果,得到准三维/二维模型的区域构造和各个畸变因子, 应用于川西—藏东MT剖面资料的分析,结果表明：该地区的构造特点为浅部比较简单,深部相对复杂;以雅安为界,东部构造比西部构造简单,而且比较稳定;构造主轴方位基本在近南北向,其偏差很少超过30°;分解参数对几条大的断裂有较好的显示;测线1横穿3个主要构造单元,并被几条大断裂切割,构造最复杂;测线2接近3D/2D模型;测线3具有较强的静位移畸变.     

大地电磁测深数据处理——4元素和6元素阻抗张量元素的比较研究(英文)

阻抗张量元素的计算是在大地电磁测深数据处理的重要一步。按照常规,阻抗张量被定义为以Zxx,Zxy,Zyx,和Zyy为元素的2×2矩阵。在本次研究中,6个元素的阻抗张量的计算使用了一个含有Zxx,Zxy,Zyx,Zyy,Zxz和Zyz分量的2x3矩阵。对上述两类阻抗张量元素的属性进行了分析。利用由印度古吉拉特邦卡奇沉积盆地采集的5个分量大地电磁数据测试了文中的方法。从视电阻率和相位的计算中我们是观察到在大部分的频带范围内4个元素阻抗和6个元素阻抗Zxy和Zyx两类元素区别不大。然而,较长周期时间的数据,如超过100秒,观察到视电阻率的增加和相位的减少。我们还注意到,倾子幅度在大部分时间几乎是零,但较长周期(超过100秒),逐渐呈增加的趋势。卡奇沉积盆地的地电断面表明在较长的周期内浅层近水平层和深层异常高电导性的不均质层都可能是引起大的Hz分量的原因。这表明,磁场垂直分量Hz对在大的2D/3D结构区域内的电场参数估计发挥的重要作用。     

区域三维大地电磁阻抗张量畸变效应的Mohr圆分析

根据二维大地电磁阻抗张量分解原理，推导出剪切、扭曲和各向异性因素影响区域性三维响应的关系式。这既适用于任意复杂的三维构造，也包含了二维和一维情况。并用Ｍｏｈｒ圆分析了受畸变因素影响时的三维响应及其在微弱畸变与强烈畸变影响下不同的性质和特征     

大地电磁资料精细处理和二维反演解释技术研究（四）——阻抗张量分解的多测点-多频点统计成像分析

基于大地电磁阻抗张量分解技术,本文提出了两种电性主轴方位的统计描述图像：随频率变化的统计分布成像（频率分布云图）和随测点序列变化的统计分布成像（测点分布云图）.这两种图像与传统的统计玫瑰图一起,较全面地描述了最佳主轴的分布特征.在进行构造维性分析过程中,通过定义二维有效因子e2d,来压制一维结构和三维结构、突出纯二维结构的影响.e2d被用于电性主轴的统计加权,有效地起到了滤波的作用;同时,统计成像中还考虑了数据质量的影响.为了得到稳定、高质量的区域阻抗张量数据,提出并实现了共主轴的多测点-多频点阻抗张量分解新算法.最终,完成了以上各项处理手段的可视化实现.本文通过两个理论模型和一个实测算例,以共轭阻抗法（CCZ法）为基础,展示了这一新技术的有效性.     

大地电磁资料精细处理和二维反演解释技术研究(五)--利用阻抗张量成像识别大地线性构造

大地线性构造往往与块体或物性的边界、断裂的几何延伸相关联,是地学研究的重要目标构造.在大地电磁法中,传统的方法是通过反演得到的电性结构的差异来解译线性构造,是一种间接分析方法.由于线性构造发育处走向和倾向明显,二维性较其他地方一般会更强,故而在阻抗张量成像分析中,相应的二维有效因子就会比较大,最佳主轴统计分布会更加明显.据此,本文利用大地电磁阻抗张量成像（包括最佳主轴统计成像、构造维性成像）来探测识别大地线性构造.本文首先通过三维理论模型的研究结果来验证这一技术在线性结构探测中的有效性,然后将其应用于郯庐断裂带南段实测剖面的具体分析中,展示了这一技术在线性构造探测中的实用性及其数据处理与分析过程.本文研究表明,大地电磁阻抗张量成像技术可独立于反演为地质解释尤其是线性构造的解释提供诸多参考信息.

     

长白山天池火山区大地电磁测深初步观测

在长白山火山区沿着两条剖面进行了１５个测点的大地电磁测深观测，对地磁感应矢量和地下电性分布维数特征进行了分析，并使用阻抗张量分解和二维自动反演等新技术对资料进行了解释。结果表明，火山区及其周围地区的地下电性结构具有二维特征，但在火山口及其附近地下存在局部三维异常体，异常体埋深约１２ｋｍ，初步推测可能是火山岩浆囊     

基于混合优化算法的MT阻抗张量畸变分解方法

在前人研究基础上,对Groom-Bailey(GB)张量分解畸变因子和区域阻抗的求解方法进行了改进.首先,通过Swift旋转与GB分解的扭变和剪切矩阵的求逆变换,利用变换后区域阻抗主对角元素为0的条件获得关于扭变因子和剪切因子的超定方程组,采用模拟退火全局优化算法进行求解.其次,由得到的扭变因子和剪切因子,结合Swift旋转确定的走向角和区域阻抗元素的估计,作为非线性最小二乘局部优化算法的初始值,对GB分解定义式的超定方程组进行求解,得到各畸变参数和区域阻抗的解.通过模型试验验证了方法的正确,对方法的稳定性进行了比较与评价,并通过与已有结果的对比和实际资料的应用,表明了方法实际应用的效果.     

网式大地电磁阻抗张量及其影响因素分析

网式大地电磁（Network-MT,N-MT）法采用长数公里至数十公里的电话线为电极线测量电场,很难形成两条笔直且相互垂直的电极线,因此阻抗张量的计算不如大地电磁法中直接.本文依据阻抗张量的旋转规则提出了一种计算N-MT阻抗张量的简便算法.依据该算法计算了中国东北地区5个N-MT测站的阻抗张量,获得了基于阻抗张量的视电阻率、相位曲线和最佳主轴方位角分布图像,为我国东北N-MT资料的进一步处理和解释提供了基础性数据.此外,本文对比分析了朝阳测站中6条N-MT观测电极线上阻抗张量旋转值与观测值之间的差异,重点讨论了产生这种系统性偏差的各种因素,提出其主要因素可能来自“电场等效各向异性”效应,即测站附近的地壳内部存在与观测电极线尺度相比拟的横向非均匀构造,而测站各电极线沿不同方向跨越不均匀构造,此时各电极线上的电场分量不遵循同一电场矢量的分解准则,导致地表观测三角形内阻抗分量不满足统一的阻抗张量旋转规则.