大地电磁资料精细处理和二维反演解释技术研究一——阻抗张量分解与构造维性分析(1)

合集下载

大地电磁阻抗计算方法进展

大地电磁阻抗计算方法进展

大地电磁阻抗计算方法进展马海荣;王世彪;郭荣文;柳建新【摘要】大地电磁法主要依据视电阻率及相位曲线的形态特征来分析大地电性分布,因此阻抗张量的计算十分重要.在实际探测中,噪声的存在使得计算得到的阻抗张量往往存在一定的偏差,干扰了数据解释的准确性.本文概述了大地电磁阻抗张量计算方法的发展,首先介绍了现在常用的几种大地电磁阻抗计算方法及其不足之处,然后介绍了如何利用两种Jackknife方法来判断阻抗估计的稳定性.【期刊名称】《工程地球物理学报》【年(卷),期】2016(013)004【总页数】6页(P423-428)【关键词】大地电磁法;阻抗张量;Jackknife【作者】马海荣;王世彪;郭荣文;柳建新【作者单位】中南大学地球科学与信息物理学院,湖南长沙410083;中南大学有色资源与地质灾害探查湖南省重点实验室,湖南长沙410083;广东省有色地质环境中心,广东广州510080;中南大学地球科学与信息物理学院,湖南长沙410083;中南大学有色资源与地质灾害探查湖南省重点实验室,湖南长沙410083;中南大学地球科学与信息物理学院,湖南长沙410083;中南大学有色资源与地质灾害探查湖南省重点实验室,湖南长沙410083【正文语种】中文【中图分类】P631.3大地电磁测深法(Magnetotelluric, MT)诞生于20世纪50年代,随着技术及仪器的发展,已成为油气田、煤炭普查、岩石圈结构探测、地下水与地热资源寻找、工程勘探等领域的一种重要手段,是一种应用较广的地球物理方法。

在实际探测中,大地电磁往往受到较强的电磁噪声干扰,且现在的勘探环境相对复杂,勘探难度相对加大,如何获得质量较好的原始数据及对原始数据进行有效的处理,对后续反演等工作有着积极的作用。

在工程物探中,大地电磁阻抗计算采用传统的最小二乘法,但是因为大地电磁测深应用的是天然场源,易受到外界噪声的干扰,为了提高信噪比,获得无偏的阻抗估计,提出了两种改进方法:一是Robust稳健阻抗估计法[1-3],它的核心是减少较大异常值对整体数据的影响;二是远参考法[4],它的核心是设置参考点,消除输入场(磁场)的噪声带来的偏移。

大地电磁阻抗张量旋转方法和曲线圆滑方法的比较

大地电磁阻抗张量旋转方法和曲线圆滑方法的比较

大地电磁阻抗张量旋转方法和曲线圆滑方法的比较赵维俊;孙中任【摘要】随着三维深部地质填图和含页岩气盆地潜力的评价项目展开,大地电磁方法的重要性越来越突显.大地电磁资料处理和解释中,至今还存在几个挑战性难题,例如阻抗张量旋转、TE和TM模式识别、静态矫正、TE和TM模式选择二维反演问题.从WinGlink软件大地电磁模块出发,使用5种阻抗旋转方法和3种数据圆滑方法处理在内蒙古扎鲁特旗采集的大地电磁数据.以数据质量不同的3个测点的处理结果为依据,对各种旋转方法和圆滑方法进行详细评价,数据极化图显示地下结构模型是复杂的三维结构.对这些方法的比较研究,对三维结构下大地电磁数据处理,特别是对使用WinGlink软件处理大地电磁数据具有重要现实意义.【期刊名称】《物探与化探》【年(卷),期】2013(037)006【总页数】8页(P1125-1132)【关键词】大地电磁;阻抗张量旋转;曲线圆滑;WinGlink软件【作者】赵维俊;孙中任【作者单位】沈阳地质矿产研究所,辽宁沈阳 110034;沈阳地质矿产研究所,辽宁沈阳 110034【正文语种】中文【中图分类】P631大地电磁测深[1-2]是一种被动源电磁测深方法。

它利用在地表测量随着时间变化的电场和磁场来推断地下电阻率分布的电磁方法。

大地电磁(MT)方法在石油天然气勘探[3]、矿产勘查[4]和地热勘查[5]中得到了广泛应用。

近年来,中国地质调查局开展了三维地质填图和深部地质调查试点工作,MT这种球物理工作的重要作用越来越突出,而随着对能源需求的日益增加,非常规能源页岩气的勘探与开发已经提上日程,MT在页岩气勘探中也必将起到重要作用。

但是,在MT 资料处理与反演中,还存在一些难题没有很好地解决,例如大地电磁数据二维反演的数据旋转方向选择[6]、反演数据极化模式选择[7-8]等。

在中国地质调查局试点项目“松辽盆地外围深部地质调查”支持下,2012年在内蒙古扎鲁特地区采集了5条剖面的大地电磁数据。

大地电磁资料精细处理和二维反演解释技术研究一——阻抗张量分解与构造维性分析(1)

大地电磁资料精细处理和二维反演解释技术研究一——阻抗张量分解与构造维性分析(1)
表明观测阻抗和区域阻抗之间存在特定的关系这种关系不受局部畸变的影响而且不需要区域结构维性的假设故可用来消除局部畸变为叙述方便我们将式称为大地电磁阻抗张量的共轭阻抗变换按照上述共轭阻抗变换规则对于区域结构是一维模型有犣1d犣1d而对于区域结构是二维模型假设观测坐标系为构造主轴方向阻抗用表示转换阻抗用犣tp表示则有可以发现转换阻抗只包含相位的信息亦即对于一维二维模型而言经过共轭阻抗变换后我们可以提取完整的阻抗张量元素的相位信息而其幅度信息则全部遗失在这一点上与bibby等提出的相位张量技术一样此外对于一维二维模型而言当观测坐标系与构造主轴一致时转换阻抗具有反对角元素为利用这个特征我们可以采用与经典swift方法完全相同的方法非常方便地求得区域最佳主轴方位畸变因子区域阻抗等重要参数完成整个张量分解任务22大地电磁阻抗张量分解的共轭变换阻抗法以上研究已经为我们提出了一种消除局部畸变的大地电磁阻抗张量分解新方法只要把实际观测的阻抗张量按照式实行共轭阻抗变换在新的转换阻抗的基础上按照与经典的swift完全相似的办法求得原阻抗张量的一系列重要参数消除电场局部畸变的影响我们把这种新的张量分解方法称为共轭变换阻抗法221确定最佳主轴方位对于转换阻抗按照阻抗旋转公式也可获得平面上任意坐标系的转换阻抗现借鉴swift旋转准则使得旋转后目标函数min则可得到最佳主轴方位的求解公式为10可得到个解按照与swift最佳主轴完全相似的方式通过比较目标函数值的大小来确定正确的结果222确定区域阻抗相位求出区域结构主轴方向后将实测的转换阻抗min最终求得区域相位的表达联立起来共有16种组合区域阻抗相位为其中使目标函数223确定畸变因子和区域阻抗张量前人的研究已经表明畸变矩阵可以分解为可确定部分和不可确定部分18由于在主轴观测坐标系下受静位移影响的二维阻抗张量与区域阻抗张量具有形式上的等效性故畸变矩阵中所包含的静位移因子无法通过张量分解来进行参数分离导致的对角元素不能严格确定21从这个意义上看静位移效应是局部畸变效应的一部分但这部分效应采用阻抗张量分解技术无法解决由于静位移因子无法确定阻抗张量分解技术的应用得到了极大的限制在实际资料处理中如果要采用张量分解得到的区域阻抗张量进行反演计算还必须极其仔细地对每个频点的视电阻率曲线进行静位移校正在静位移因子不能确定的前提下为了求得畸变矩阵和区域阻抗张量需要引入一定的约束条件我们假定如果不存在局部三维小异常

基于二维小波分解的MT二维反演方法

基于二维小波分解的MT二维反演方法

基于二维小波分解的MT二维反演方法刘苗;吴健生;于鹏;张罗磊【摘要】对于大地电磁反演,期望得到快速稳定的反演结果并且能够反映清晰的、未发生畸变的地质界面.为了解决这一问题,提出了基于二维小波分解的MT二维反演技术,对小波分解得到的背景场采用光滑反演,细节部分采用聚焦反演.通过模型试验,证明了该方法的有效性,同时对某区的实际资料进行了处理,表明了该方法的适用性和效果.【期刊名称】《物探与化探》【年(卷),期】2009(033)003【总页数】5页(P337-341)【关键词】大地电磁测深;二维小波分解;分频反演【作者】刘苗;吴健生;于鹏;张罗磊【作者单位】同济大学,海洋地质教育部重点实验室,上海,200092;中国石化集团,上海海洋石油局,上海,200092;同济大学,海洋地质教育部重点实验室,上海,200092;同济大学,海洋地质教育部重点实验室,上海,200092;同济大学,海洋地质教育部重点实验室,上海,200092【正文语种】中文【中图分类】P6312D大地电磁反演方法有RRI(rapid relaxation inversion,快速松弛反演),OCCCAM(极小可能构造意义下符合数据的模型反演), NLCG(nonlinearconjugate gradients,非线性共轭梯度法),REBOCG(rduced basis occam's inversion,简化的OCCAM反演)和SBI(sharp boundary inversion,尖边界反演)等[1-3]。

这些反演方法各有侧重点,如OCCAM反演偏重于最小结构的光滑模型反演,RRI反演偏重于反演速度,SBI反演偏重于反演具有确切地质意义的目标体及其界面。

另外还有相关的非线性反演方法,如蒙特卡洛法、人工神经元反演法、模拟退火反演(SA)和基于小波分析的多尺度反演法等,并且开始用于实际MT资料的处理解释[4-6]。

以往的多尺度反演是把目标函数分解成不同尺度的分量,根据不同尺度上目标函数的特征逐步缩小全局最小点。

基于大地电磁三维响应的二维反演解释研究

基于大地电磁三维响应的二维反演解释研究

基于大地电磁三维响应的二维反演解释研究作者简介:周汝峰,男,1992年102017年6获得工学硕士学位。

摘要实际的地电构造是三维的,故应在大地电磁中采取三维的反演方法来得到地下的电性结构,但现在大地电磁三维的反演因为存在一些问题而仍未得到广泛使用。

故研究大地电磁三维响应的二维反演解释是较有意义的。

本文首先介绍分析大地电磁三维正演基本的一些理论和公式、三维响应特征,然后简要介绍了二维反演的一些基本理论和公式并对比介绍了不同二维反演方法的特点及其在不同极化模式下的特征;其后对大地电磁的阻抗张量进行了理解和学习,重点学习了相位张量分解的特点和用法;然后建立了一个三维空间中存在不同电阻率块体的典型三维模型,以分析其在不同极化模式下的三维响应特征,深刻研究其特征差异,即不同极化模式对电性分界面分辨的灵敏度,并探究了其原因;建立了另外一个空间中存在一个低阻块体的三维模型,以分析其在不同极化模式下的二维反演结果,探究了极化模式的选择对二维反演的影响;最后综合利用上述的研究成果对两个较复杂的模型进行了试算,采取阻抗相位张量分解,进行构造维性分析、寻求最佳旋转角,并将旋转后的数据作二维反演,探求了一种大地电磁资料处理解释的合理流程。

发现对于同一异常体,不同极化模式的正演响应是不一样的,其不同极化模式的反演结果也会出现一些差异,有时差异较大,也就是说不同的极化模式对电性分界面的分辨能力是不一样的;在野外工作,布设测线时,应当尽量垂直于构造的走向,并且采取合适的极化模式来进行反演。

并且对于三维响应尤其是实测资料很容易出现受浅表低阻体等影响而引起畸变的情况,或者遇到采样方向不与构造走向垂直的情况,假设不对观测资料进行任何的校正,而直接做反演解释,那么所得到的结果将会出现偏离实际构造的情况;另外,假如采用的校正方法不恰当的话,也会使得后续的反演结果得到错误的结论;故应首先选取合适的地球物理模型、张量分解方式,进行阻抗张量分解、构造维性分析、求取最佳旋转角,并将数据旋转到最佳角度,再进行合适极化模式的二维反演,最终将可能消除冗余构造、达到合理的反演效果。

大地电磁三维模型的一维二维反演近似问题研究

大地电磁三维模型的一维二维反演近似问题研究

文章编号:1672—7940(2006)01—0009—07大地电磁三维模型的一维二维反演近似问题研究陈小斌1,赵国泽1,马 霄2(1.中国地震局地质研究所,北京100029;2.大港油田集团石油工程监理中心,天津300280)作者简介:陈小斌(1972—),男,中国地震局地质研究所副研究员,主要从事电磁测深正反演方法及地球动力学数值模拟方面的研究。

E -mail :cxb @pku 1edu 1cn摘 要:设计了一个简单的三维模型———均匀半空间的三维低阻异常体模型。

采用51×51×38的网格计算得到了331333Hz 至010071Hz 之间共34个频率的三维正演响应。

选取其中5条具有代表性的测线,比较了不同测线上自适应正则化(ARIA )一维反演结果和非线性共轭梯度法(NL CG )二维反演结果相互之间及其与原三维理论模型之间的异同,研究大地电磁三维模型的一维、二维反演的近似情况。

发现:对于三维模型响应,一维、二维反演都能得到反映模型真实信息的结果,但二维模型对三维模型的近似程度比一维模型要高得多;在二维模型反演中,采用TM 模式更易得到真实的模型信息。

研究中,采用蒲兴华等(1994)开发的M T 三维正演程序进行正演计算,采用“大地电磁数据管理器(M TDA TAMN G )”软件进行数据处理和一维、二维反演。

关键词:大地电磁;三维;一维;二维;反演中图分类号:P631124文献标识码:A收稿日期:2005—12—18RESEARCH ON INVERSION OF MT3D MODE L APPR OXIMATE LY B Y 1D ,2D INVERSION METH ODC H EN Xiao 2bin 1,ZHAO Guo 2ze 1,MA Xiao 2(11I nstitute of Geolog y ,China Seismological B ureau ,B ei j ing 100029,China;21Dagang Pet roleum S u pervision Cent re ,Tianj ing 300280,China )Abstract :In t his paper ,a simple 3d model ,low -resistivity abnormity 3d model in halfspace ,is designed 1Responses of 34frequencies from 331333Hz to 010071Hz are obtained in3d modeling wit h 515138grids 15rep resentatio nal lines are adopted to st udy approxima 2ting 3d model wit h 1d ,2d inversion result s 1Result s of 1d inversion ,2d inversionfor each line are comp uted ,and are compared wit h each ot her and wit h t he original model 1St udies show t hat bot h 1d and 2d inversion result s contain t rue 3d model information ,but t he 2d inversion is far more app roached to 3d model t han 1d inversio n 1Furt hermore ,when app roaching 3d models wit h 2d inversion ,it is t he TM polarization modedata rat her t han t he TE polarization mode data t hat sho uld be adopted 1The 3d modelingprogram used in t his paper was developed by Dr 1PU Xinghua 1The new visual integratedsystem ,M TDA TAMN G ,is used in data processing and inversion 1The 1d inversion met h 2第3卷第1期2006年2月 工程地球物理学报CHIN ESE J OU RNAL OF EN GIN EERIN G GEOP H YSICS Vol 13,No 11Feb 1,2006od is adaptive regularized inversion algorit hm(A RIA),and2d inversion met hod is nonlin2 ear conjugate gradient s algorit hm(NL C G)1K ey w ords:M T;3d;1d;2d;inversion1 引 言当前,大地电磁测深的一维和二维反演技术已相对成熟,并成为实际资料处理和解释中的主要手段。

基于多测点多频点的大地电磁一维区域阻抗恢复

基于多测点多频点的大地电磁一维区域阻抗恢复

基于多测点多频点的大地电磁一维区域阻抗恢复王琪琪;马涛;王红兵;谭鑫;高鑫;汪强强【期刊名称】《工程地球物理学报》【年(卷),期】2022(19)4【摘要】在实际勘探中,浅部三维电性不均匀体的存在会导致大地电磁响应发生畸变,使阻抗张量偏离真实区域阻抗,影响后续反演和解释工作。

本文主要利用1D(One-dimensional)区域构造中静态因子的统计学性质对多测点多频点畸变校正方法进行了改进。

首先通过粒子群全局优化算法得到GB(Groom-Bailey,GB)分解参数中确定参数的全局最优解空间,将其作为高斯牛顿法的初值来进行局部优化,然后通过ssq(the Sum of the Squared Elements,SSQ)旋转不变量的统计特性得到不确定参数,最后利用正演模型验证了该算法的有效性。

结果表明:该方法结合了粒子群算法和高斯牛顿法,避免结果陷入局部最优,通过正演模型发现不确定参数随测点位置呈“正态分布”,在异常体中心位置取最大,并且不确定参数可以指示静态效应的影响频段。

同时静态校正后视电阻率和相位曲线与未受畸变影响的数据具有较高的重合度,说明该方法很好地消除了1D区域构造中多个浅部异常体的影响,恢复了原始区域阻抗信息。

【总页数】8页(P566-573)【作者】王琪琪;马涛;王红兵;谭鑫;高鑫;汪强强【作者单位】中铁第一勘察设计院集团有限公司陕西铁道工程勘察有限公司【正文语种】中文【中图分类】P631.3【相关文献】1.大地电磁资料精细处理和二维反演解释技术研究(四)——阻抗张量分解的多测点-多频点统计成像分析2.区域三维大地电磁阻抗张量畸变效应的Mohr圆分析3.基于K中心点聚类分析的大地电磁阻抗识别4.大地电磁测深首枝频点统计求平均法在连片静位移区的校正实验研究5.对大地电磁测深阻抗张量的解释:区域感应及局部大地电流畸变因版权原因,仅展示原文概要,查看原文内容请购买。

大地电磁资料精细处理和二维反演解释技术研究(二)——反演数据极化模式选择

大地电磁资料精细处理和二维反演解释技术研究(二)——反演数据极化模式选择

Absr c Ai n tt hoc fd t o a ia i n mo e i g t t Iu i D n e so ta t mi g a he c ieo a a p l rz to d n ma ne o e I rc 2 i v r i n,s v r l e e a 2 a D D nd 3 mod l r e i ne o s u y t if r n ebe we n 3 r s o e a D e po s ,t e sa e d sg d t t d he d fe e c t e D e p ns nd 2 r s n e he
用T E模 式 或利 用 T +TM 模式 联 合 反 演 都更 合 理 , 演 结 果 中 的虚 假 结 构 明显 减 少 ; 于 T 模 式 , 位 受 三维 E 反 对 M 相
畸变 影 响 较小 , 电阻 率较 大 , 以二 维 反演 中可 适 当加 大 相 位 的 权 ; 于 实测 数 据 的 二 维 反 演 , 优 先 考 虑 采 用 视 所 对 应 TM 模式 数 据 进行 二 维 反 演 , 次是 TM+T 其 E模 式 , 般 不 要单 独采 用 T 一 E模 式 . 关键 词 大地 电磁 , 化模 式 , 极 三维 模 拟 , 维 反 演 二
W hi h da a p l r z to o e s o l e u e n 2 i e s o c t o a i a i n m d h u d b s d i D nv r i n
CAIJ n Ta u — o ~.CHEN XioBi a— n
( 1 : 7 3 2 1 DOI 1 . 9 9 jis . 0 15 3 . 0 0 1 . 1 1 ) 2 0  ̄ 7 4, :0 3 6 /.sn 0 0 — 7 3 2 1 . 1 0 8 C i a T,Ch n X . Reie e h iu sfrd t r c sig a d t — i n in liv rin i g eo elrc I: hc aa J e B f d tc nq e o aa p o esn n wodme so a n eso n ma n t tl i n u 1 W ih d t

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

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

一种简单的大地电磁阻抗张量畸变分解方法王立凤;晋光文;孙洁;白登海【期刊名称】《地震工程学报》【年(卷),期】2001(023)002【摘要】提出了一种简单的大地电磁阻抗张量畸变分解方法.引入3d/2d模型,使得全畸变的大地电磁阻抗张量被简化.通过对比Swift旋转和Bahr分解后的资料,用联系二者的畸变矩阵表示浅层三维电性异常体对区域二维构造的局部畸变影响.通过研究畸变张量矩阵,并对其进行分解,可将局部三维电性不均匀体产生的畸变与区域构造分离开来,从而了解区域构造的维数特征和局部畸变类型及其对区域构造的影响.最后,应用3个测点的实际资料,对这种简单的畸变分解方法进行了检验.结果表明,该方法简单且应用方便,解释效果较好.%The necessity of establishing a 3d/2d model is discussed. Owing to introducing the model, the full electric and magnetic field galvanic distortion tensor is simplified. The effects of a surface 3-d conductivity anomaly on a regional 2-d structure is derived by comparing data through Swift rotation with data through Bahr decomposition, then is expressed by a distortion matrix. The effect of a local anomaly can be separated from a regional structure by researching and decomposing the distortion tensor matrix. So it is available to find out the regional structure dimensionality, classification of the local structure and the effect of it on the regional structure. Eventually, the simple distortion decomposition method is tested through three sites of practical data.【总页数】9页(P172-180)【作者】王立凤;晋光文;孙洁;白登海【作者单位】中国地震局地质研究所,;中国地震局地质研究所,;中国地震局地质研究所,;中国地震局地质研究所,【正文语种】中文【中图分类】P631.3+25【相关文献】1.大地电磁阻抗张量格林—贝利分解方法及其应用 [J], 李伟;欧阳明明;周成涛2.基于遗传算法的大地电磁阻抗张量分解方法研究 [J], 尹曜田;魏文博;叶高峰;金胜;董浩3.表面局部三维大地电磁曲线变校正:MT畸变校正阻抗张量分解方法 [J], 王书明4.基于改进粒子群算法的大地电磁阻抗张量分解方法 [J], 陈先洁;王绪本;李德伟;谢卓良;乃国茹5.相位张量分析约束下的大地电磁测深阻抗张量分解方法研究 [J], 谢成良;魏文博;金胜;叶高峰;景建恩;张乐天;董浩;张帆;王辉;姚硕因版权原因,仅展示原文概要,查看原文内容请购买。

二维与三维大地电磁反演的敏感度研究

二维与三维大地电磁反演的敏感度研究

二维与三维大地电磁反演的敏感度研究孟绿汀;张慧茜;黄清华【期刊名称】《地球物理学报》【年(卷),期】2022(65)8【摘要】大地电磁法可应用于岩石圈结构探测、金属固体矿产以及地热资源勘察等方面,目前三维大地电磁反演技术日趋成熟,但二维反演依然应用广泛.受台网布设及复杂地质构造影响,大地电磁数据往往结合二维与三维反演技术来减轻计算负担和改善反演效果,而评估二维与三维反演模型及其差异对准确获取地下电性结构信息至关重要.本文提出了一种基于敏感度矩阵的电性异常体边界划定以及反演模型可靠性评估方法,该方法通过求解大地电磁正演的伴随问题来构建敏感度矩阵,无需设计额外模型即可分析反演结果的可靠性.此外,本文还利用合成数据讨论了大地电磁二维与三维反演的有效探测深度.当地下结构三维性不强时,二维反演可以有效恢复深部低阻异常体的电阻率值,为深部地质解释提供定量约束.但当电性结构三维性较强时,三维反演重建的边界与真实位置更加吻合,而二维反演会产生虚假异常,不利于后续地质解释.【总页数】10页(P3078-3087)【作者】孟绿汀;张慧茜;黄清华【作者单位】北京大学地球与空间科学学院地球物理系;河北红山巨厚沉积与地震灾害国家野外科学观测站(北京大学)【正文语种】中文【中图分类】P319【相关文献】1.大地电磁三维模型的一维二维反演近似问题研究2.大地电磁资料精细处理和二维反演解释技术研究(六)——交错模型的大地电磁二维反演3.大地电磁资料精细处理和二维反演解释技术研究(三)——构建二维反演初始模型的印模法4.利用大地电磁三维反演方法获得二维剖面附近三维电阻率结构的可行性5.大地电磁三维反演在土屋铜矿床探测中的应用研究--三维反演和解释因版权原因,仅展示原文概要,查看原文内容请购买。

MT阻抗张量在二维与三维介质中的分析

MT阻抗张量在二维与三维介质中的分析

MT阻抗张量在二维与三维介质中的分析胡玮哲;熊高君【摘要】大地电磁测深法主要研究岩石圈的电性结构,以及对高阻覆盖层下的低阻矿物进行探测.大地电磁测深法利用天然交变的平面电磁波(10-4Hz~104Hz)为场源,依据在地表观测到的包含地下地质体信息的电场值与磁场值计算出阻抗张量,并由阻抗张量计算出的视电阻率以及相位曲线特征来分析大地电性结构特征.通过对阻抗张量的旋转分析、SVD 分解、通用分析以及阻抗张量扩展,来获得阻抗张量元素的地质信息,并提出了解决三维介质阻抗张量的计算方法.%Magnetotelluric sounding method has high resolving power in studying the electrical structure of lithosphere and low resistivity minerals under high resistivity cover.The use of natural alternating electromagnetic wave(10-4~104Hz)field source,based on electric field contains the underground geological information on the surface resistivity observation value calcu-lated by impedance tensor and magnetic field values,and the apparent resistivity and phase curves are calculated to analyze the electrical structure of the earth.The present research and analysis of the impedance tensor is of great significance.This paper through the rotation analysis,SVD decomposition,general geological information analysis and extension of impedance tensor magnetotelluric impedance tensor element,calculation method is proposed to solve the three-dimensional dielectric impedance tensor.【期刊名称】《物探化探计算技术》【年(卷),期】2018(040)002【总页数】7页(P205-211)【关键词】大地电磁;阻抗张量;电阻率【作者】胡玮哲;熊高君【作者单位】成都理工大学地球物理学院成都 610059;成都理工大学地球物理学院成都 610059【正文语种】中文【中图分类】P631.20 引言大地电磁测深法(Magnetotelluric ,MT)是20世纪50年代初,由前苏联学者A.T.Tikhonov[1]与法国地球物理学家L.cagnird[2]分别提出,以天然交变电磁场为场源来获取有关地下地质体的电性结构特征的信息,属于频率域测深的一种方法。

二维大地电磁资料的偏移

二维大地电磁资料的偏移

二维大地电磁资料的偏移
魏胜;王家映
【期刊名称】《地球物理学报》
【年(卷),期】1993(000)002
【摘要】本文采用反射映像(即U/D成像)原理对大地电磁资料进行成像,将地面观测的波场分解成上行波和下行波,并分别向下延拓,利用上、下行波的时间一致性原理确定地下电性界面的位置.对TE和TM模式的MT响应分别处理可以得到两个深度剖面,两剖面的一致性便为实际反射界面的最佳估计.理论和实际资料的计算表明,MT偏移技术是比较有效的,它可获得地下界面直观的映像.和常规MT反演方法相比较,该方法具有算法简单等优点,并且还能获得表示地下地质构造真实映像的MT深度剖面.
【总页数】1页(P256)
【作者】魏胜;王家映
【作者单位】不详;不详
【正文语种】中文
【中图分类】P318.6
【相关文献】
1.大地电磁资料三维傅里叶变换偏移方法及其应用 [J], 宋维琪
2.大地电磁资料偏移成像 [J], 宋维琪;邹文勇
3.大地电磁资料精细处理和二维反演解释技术研究(六)——交错模型的大地电磁二
维反演 [J], 郭春玲;陈小斌
4.大地电磁资料精细处理和二维反演解释技术研究(三)——构建二维反演初始模型的印模法 [J], 叶涛;陈小斌;严良俊
5.对二维偏移资料进行图偏移 [J], Whit.,DN;谭凡
因版权原因,仅展示原文概要,查看原文内容请购买。

电性可控源音频大地电磁法二维反演方法技术研究

电性可控源音频大地电磁法二维反演方法技术研究

收稿 日期 :O O 0 — 1 2 l 一 5 2 作者简介 : 汪英宏( 94 , , 15 一)男 安徽舒城人 , 工程师 , 现主要从事地球物理勘探工作 。
第 2 卷 第 3期 0
汗英宏, 电性可控源音频大地电磁法二维反演方法技术研究 等:
以及 地 磁 资 料 的分 析 ,网 格 密度 较 小 的反 演结 果 图 之 间的高阻异常 以及 1 1 号点 至大号 点处的高阻异 5
中图分类号 :P 3 _ 6 l2
文献标志码 :A
电阻率 P 由 电场 振 幅E 、磁场 振 幅 H 计 算得 之 。 x y
0 引 言
于 具 有 发射 功 率 大 、信 噪 比高 、快速 高 效 等 优点 ,
卡 尼 亚 电阻 率 P 只在 远 区 是 正确 的 ,在 过 渡 区和 近区 p 公 式无 效 ,所计 算 的视 电阻率 不能真 实反 映 图 1( b)是 C AMT 一 条 孔 旁 测 深 卡 尼 亚 电阻 率 S 法
地 电断 面 的反 演 方 法 技术 的研 究 尚不够 。本 文 用美 它 使 连续 的求 解 区域 离 散 化 ,即按 一定 的规 则将 计 国Z NG 公 司S 2 维 反演 软 件对 G .2 型仪 算 区域 剖 分 为 一些 节 点 处相 互 连 接 的 网格 单元 ,因 O E CS DZ. DP3 “
电性 可 控源 音频 大 地 电磁法 ( 简称 C AMT)由 地 下 电阻 率 变化 ,在 近 区 电阻 率 曲线 呈4 o 直线 。 S 5
近几年来在深部找矿工作 中扮演着重要角色 。它通 常 以 剖 面性 工 作 方式 为 主 ,因此所 获 得 原 始 资料 是 二 维形 式 的卡 尼 亚 电 阻率 断 面 。这种 断 面 虽 然反 应 了地 中介 质 电性 结 构分 布 特 征 ,但 含 有 地 表 不均 匀 地 质 体 的影 响—— 静 态 效 应 以及 地形 起 伏 的影 响 ;

用大地电磁法研究构造走向及维性特征

用大地电磁法研究构造走向及维性特征

用大地电磁法研究构造走向及维性特征
杨长福;林长佑
【期刊名称】《地震工程学报》
【年(卷),期】2002(024)001
【摘要】介绍了大地电磁GB张量分解法及其对它的改进法,可确定出更可靠更真实的区域构造走向.将分解结果结合传统的座标旋转法所确定的视电阻率、相位、走向、偏离度等响应函数及维权参数进行分析,可得到更详细的电性结构维性质信息.对兰州地区的实测资料研究表明,区域电性结构主体呈2-D结构,走向方向大致为南北或东西向.
【总页数】7页(P65-71)
【作者】杨长福;林长佑
【作者单位】南京大学地球科学系,江苏,南京,210093;中国地震局兰州地震研究所,甘肃,兰州,730000
【正文语种】中文
【中图分类】P319
【相关文献】
1.天然源大地电磁法在构造勘察中的应用研究 [J], 祁晓雨;许广春;李志华
2.有限差分法二维地电构造模型大地电磁法正演模拟 [J], 张弘宇; 温建亮; 张富明
3.有限差分法二维地电构造模型大地电磁法正演模拟 [J], 张弘宇; 温建亮; 张富明
4.基于非结构网格的三维大地电磁法有限内存拟牛顿反演研究(英文) [J], 曹晓
月;殷长春;张博;黄鑫;刘云鹤;蔡晶
5.基于大地电磁法的深部地质构造勘探研究 [J], 潘双彦;郭锋;甘宇升;张波
因版权原因,仅展示原文概要,查看原文内容请购买。

用大地电磁法研究构造走向及维性特征

用大地电磁法研究构造走向及维性特征
解释, 乃是 当前 最 流 行 而 又 实 际 的研 究 途 径 ] 基 于 目前 的 这 种 现 状 , 量 地 确 定 电 性 构 造 的 . 定 维 参 数 和 区域 构 造 走 向也 显得 越 来 越重 要 . 如 人 们 必 须 在 事 先 确 定 在 误 差 范 围 内 是 否 表 现 例 强 烈 的 3D 性 ,一 一 3 D性 对 MT 响应 为 lD或 2 D假 设 的 有 效 性 的 影 响 程 度 . 果 3D性 很 强 的 — 一 如 一
西





第 2 4卷
地 表 示 区 域 3 D程 度 . 它 和 传 统 方 法 确 定 的 偏 离 度 进 行 比 较 , 以 帮 助 我 们 判 断 局 部 和 区 域 一 将 可
构 造 的维 数 . 在 模 型 较 为 简 单 的 情 况 下 , 所 确 定 的 走 向 和 维 参 数 很 不 稳 定 . 近 P ae 但 它 最 rtr和
资 料 用 2D 反 演 方 法 进 行 解 释 , 必 造 成 很 大 误 差 . 存 在 非 均 匀 性 的 情 况 下 , 传 统 的 座 标 一 势 在 用
旋 转 法 求 出 的 区 域 性 构 造 走 向 , 于表 面 非 均 匀 性 的 强 烈 影 响 , 时 是 凌 乱 或 不 正 确 的. 外 , 由 有 另 由 于 响 应 函 数 随 座 标 旋 转 角 ( 向角 ) 走 而变 化 , 此走 向角 的偏 差 , 必 影 响 到 用 于 反 演 的 响应 因 势 函数值 , 而也必将 影响反演 结果. 从
维普资讯
第 2 4卷 第 1期 20 0 2年 3月
西



大地电磁法数据模式识别与阻抗旋转在铁路地质勘察中应用

大地电磁法数据模式识别与阻抗旋转在铁路地质勘察中应用

大地电磁法数据模式识别与阻抗旋转在铁路地质勘察中应用马涛【摘要】在音频大地电磁法(AMT)的数据处理中,经常会遇到诸如模式识别、阻抗旋转等问题.重点讨论大地电磁数据处理过程中的模式识别和阻抗旋转并结合案例说明处理效果.【期刊名称】《铁道勘察》【年(卷),期】2017(043)005【总页数】4页(P90-93)【关键词】大地电磁;模式识别;阻抗旋转【作者】马涛【作者单位】中铁第一勘察设计院集团有限公司,陕西西安 710043【正文语种】中文【中图分类】P631在大地电磁法的处理解释中,采用不同的解释手段所得到的处理结果往往有很大差异,产生问题的原因在于:测区构造异常复杂,很难找到与其构造相适应的资料处理、分析与解释方法。

目前,音频大地电磁数据的后续处理步骤主要为曲线编辑、静态分析,很少做大地电磁曲线模式识别工作,大地电磁阻抗旋转处理的应用更是少之又少。

为使大地电磁解释结果更为可靠、合理,将重点讨论大地电磁数据处理过程中的模式识别和阻抗旋转,并结合案例说明处理效果。

大地电磁测深中表面阻抗Z定义为地表电场和磁场水平分量的比值Z=E/H。

实际的地质模型比较复杂,一维介质模型仅具有理论研究价值,如地质研究中经常用到的背斜或向斜构造,其电阻率分布沿走向和倾向两个方向变化。

此时,需要采用二维模型来研究。

当介质电性为多方向变化时,需要采用三维模型。

阻抗张量是地质模型的响应函数。

对于一般的三维不均匀大地体,阻抗是空间坐标的函数,阻抗张量是普通二维复数矩阵,完整的描述应当是含有四个元素的张量,每个元素与场的正交分量有关在地质模型为二维情况下,阻抗张量Z2D为其中式中,α是地质构造走向顺时针旋转到观测电场EY方向的角度。

2.1 大地电磁极化模式的定义大地电磁极化模式是根据二维模型条件下的地质体走向进行定义的[4-8]。

在工程勘察中,将构造倾向定义为X方向,将构造走向定义为Y方向,这时,Rxy、Pxy 代表TM模式,Ryx、Pyx代表TE模式(如图1)。

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

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

大地电磁GB张量分解法的改进
杨长福;林长佑;徐世浙
【期刊名称】《地球物理学报》
【年(卷),期】2002(045)0z1
【摘要】在对Groom-Baily的大地电磁(MT)张量分解法(简称GB分解法)修改之后,采用对单个频点进行张量分解的方法,并将参加分解的阻抗张量元素改用视电阻率和相位,使得分解方程数目及求解未知数数目减少,大大改善了方程的病态性能并减少计算时间,从而使张量分解更少地受分解条件的限制而稳定地进行.在此基础上,编写了单频点张量分解程序,并应用于实际资料.
【总页数】9页(P356-364)
【作者】杨长福;林长佑;徐世浙
【作者单位】南京大学地球科学系,南京,210093;中国地震局兰州地震研究所,兰州,730000;浙江大学地球科学系,杭州,310027
【正文语种】中文
【中图分类】P631
【相关文献】
1.大地电磁相位张量分解特性 [J], 陈江平;周亚东
2.大地电磁资料精细处理和二维反演解释技术研究(四)——阻抗张量分解的多测点-多频点统计成像分析 [J], 陈小斌;蔡军涛;王立凤;叶涛
3.大地电磁GB张量分解法的改进 [J], 杨长福;林长佑;等
4.关于大地电磁阻抗张量分解技术之Bahr分解稳定性与实用性的简要分析 [J], 张

5.基于改进粒子群算法的大地电磁阻抗张量分解方法 [J], 陈先洁;王绪本;李德伟;谢卓良;乃国茹
因版权原因,仅展示原文概要,查看原文内容请购买。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

犃 犫 狊 狋 狉 犪 犮 狋 h es e c i a lr e l a t i o n s h i e t w e e nt h eo b s e r v e da n dr e i o n a l i m e d a n c et e n s o r sc a nb e T p pb g p d e r i v e df r o mam a t h e m a t i c a lt r a n s f o r m a t i o nc a l l e dc o n u a t ei m e d a n c et r a n s f o r m a t i o no ft h e j g p m a n e t o t e l l u r i c i m e d a n c et e n s o r .T h i sr e l a t i o n s h i si n d e e n d e n to fg a l v a n i cd i s t o r t i o n.N o g p pi p a s s u m t i o na b o u t t h ed i m e n s i o n a l i t f t h eu n d e r l i n e i o n a l c o n d u c t i v i t t r u c t u r e i sr e u i r e d p yo y gr g ys q a n dt h i sr e r e s e n t a t i o ni sa l i c a b l et o3 Dr e i o n a lc o n d u c t i v i t t r u c t u r e .T h r o u ht h e p p p g ys g ,f t r a n s f o r m a t i o no ft h eo b s e r v e di m e d a n c et e n s o r i r s t l h er e i o n a ls t r i k ea n l ec a nb e p yt g g

通讯作者 陈小斌 , 男, 博士 , 副研究员 , 主要从事电磁测深正反演及其应用 、 地球动力学模拟等方面的研究 . 1 9 7 2 年生 , : c x b@p k u. e d u . c n E m a i l
0期 1
蔡军涛等 : 大地电磁资料精细处理和二维反演解释技术研究 ( 一) — — — 阻抗张量分解与构造维性分析
: 犚 犲 犳 犻 狀 犲 犱 狋 犲 犮 犺 狀 犻 狌 犲 狊 犳 狅 狉犱 犪 狋 犪狆 狉 狅 犮 犲 狊 狊 犻 狀 狀 犱 狋 狑 狅 犱 犻 犿 犲 狀 狊 犻 狅 狀 犪 犾 犻 狀 狏 犲 狉 狊 犻 狅 狀 犻 狀犿 犪 狀 犲 狋 狅 狋 犲 犾 犾 狌 狉 犻 犮 Ⅰ 狇 犵犪 犵 犜 犲 狀 狊 狅 狉犱 犲 犮 狅 犿 狅 狊 犻 狋 犻 狅 狀犪 狀 犱犱 犻 犿 犲 狀 狊 犻 狅 狀 犪 犾 犻 狋 狀 犪 犾 狊 犻 狊 狆 狔犪 狔
大地电磁资料精细处理和二维反演解释技术研究 ( 一)
— — — 阻抗张量分解与构造维性分析
2 蔡军涛1, ,陈小斌1 , 赵国泽1 北京 1 1 中国地震局地质研究所 , 0 0 0 2 9 南京 2 1 0 0 1 4 2 江苏省地震局 ,
摘 要 本文对大地电磁观测阻抗实施一种数学变换 — — —“ 共轭阻抗变换 ” , 发现 转 换 后 的 观 测 阻 抗 和 区 域 阻 抗 之 这种特定关系不受电场局部畸变的影响 , 而且同样不需要关 于 地 下 区 域 结 构 维 性 的 假 设 , 在区 间存在特定的关系 , 域结构是三维的情况下也是成立的 . 对转换后的观测阻抗采用已有的 S 然 w i f t旋转方法即可求得区域主轴 方 位 角 , 后采用最优化方法求取区域阻抗相位 、 振幅以及畸变 因 子 . 同时根据转换后的观测阻抗重新定义了不受电场局部 畸变影响的构造维性参数 . 采用合成理论数据验证 了 新 算 法 的 正 确 性 , 和S 、 、 相 位 张 量、 w i f t B a h r G B、 WA L方法进 并将新方法应用于实测资料的解释 , 发展 了 一 套 消 除 局 部 畸 变 和 进 行 构 造 维 性 分 析 的 MT 精 细 资 行了对比分析 , 料处理技术 . 进一步的工作是将新算法推广到多点 多 频 , 并 采 用 最 优 化 技 术 分 解, 发 展 一 种 稳 定 性 好、 多测点多频 点的阻抗张量分解技术 . 关键词 大地电磁 , 局部畸变 , 共轭阻抗变换 , 张量分解 , 构造维性分析 : / D O I 1 0 . 3 9 6 9 . i s s n . 0 0 0 1 5 7 3 3 . 2 0 1 0 . 1 0 . 0 2 5 j 中图分类号 P 6 3 1 收稿日期 2 , 0 1 0 0 8 2 4 2 0 1 0 0 9 3 0收修定稿
第5 3卷 第1 0期 2 0 1 0年1 0月
地 球 物 理 学 报
CH I N E S E J OUR NA L O F G E O P HY S I C S
V o l . 5 3,N o . 1 0 , O c t . 2 0 1 0
蔡军涛 ,陈小斌 , 赵国泽 . 大地电磁资料精细处理和二维反演解释技术研究 ( 一) — — — 阻抗张量分解与构造维性 分 析 . 地球物理 学报 , ( ) : : / 2 0 1 0, 5 3 1 0 2 5 1 6~2 5 2 6, D O I 1 0. 3 9 6 9 . i s s n . 0 0 0 1 5 7 3 3. 2 0 1 0. 1 0. 0 2 5 j C a i JT, C h e nX B, Z h a oG Z.R e f i n e dt e c h n i u e sf o rd a t ap r o c e s s i n n dt w o d i m e n s i o n a l i n v e r s i o ni nm a n e t o t e l l u r i c Ⅰ: q ga g ) , : : T e n s o rd e c o m o s i t i o na n dd i m e n s i o n a l i t n a l s i s . 犆 犺 犻 狀 犲 狊 犲犑. 犌 犲 狅 犺 狊 .( i nC h i n e s e 2 0 1 0, 5 3( 1 0) 2 5 1 6~2 5 2 6, D O I 1 0. p ya y 狆 狔 / 3 9 6 9 . i s s n . 0 0 0 1 5 7 3 3. 2 0 1 0. 1 0. 0 2 5 j
基金项目 国家自然科学基金项 目 ( ; 国家重点基础研究发展计划项目“ 项目( ; 国家自然科学基金重点项目 4 0 6 7 4 0 3 6) 9 7 3” 2 0 0 4 C B 4 1 8 4 0 2) ( ) 资助 . 4 0 5 3 4 0 2 3 作者简介 蔡军涛 , 男, 博士 , 助理研究员 , 主要研究方向为电磁方法与应用 . : 1 9 7 7 年生 , E m a i l c u n t a o o h u . c o m @s j
2 5 1 7
,t ,a d e t e r m i n e du s i n o n v e n t i o n a lS w i f tm e t h o d .S e c o n d l h ep h a s e s m l i t u d eo fr e i o n a l gc y p g i m e d a n c et e n s o ra n d d i s t o r t i o n p a r a m e t e r sa r eo b t a i n e d b t i m i z a t i o n m e t h o d .N e w p y o p d i m e n s i o n a l i t a r a m e t e r so ft h eg e o e l e c t r i cs t r u c t u r ew h i c ha r ef r e eo fg a l v a n i cd i s t o r t i o na r e yp d e f i n e db a s e do n t h e t r a n s f o r m a t i o no f t h eo b s e r v e d i m e d a n c e t e n s o r s . S n t h e t i cd a t a a r eu s e d t o p y ,B ,G a h r B,p h a s e c h e c kt h ev a l i d i t fn e w m e t h o d sa n df o rc o m a r a t i v ea n a l s i so fS w i f t yo p y ,WA , t e n s o r L, a n dn e wa r o a c h. F i n a l l t h en e wm e t h o d s a r e a l i e d t o t h e f i e l dd a t a a n da s e t p p y p p o f r e f i n e dp r o c e s s i n e c h n i u e sf o rm a n e t o t e l l u r i cd a t a i n c l u d i n e m o v i n o c a ld i s t o r t i o na n d gt q g gr gl d i m e n s i o n a l i t n a l s i sa r ed e v e l o e d .T h ef u t u r ew o r ki st oe x t e n dt h i sn e ws t r a t e o ya y p g yt ,m m u l t i s i t e u l t i f r e u e n c e c o m o s i t i o nb s i n t i m i z a t i o nm e t h o da n dd e v e l o o r e s t a b l e q yd p yu go p pam ,m u l t i f r e u e n c e n s o rd e c o m o s i t i o nf o rm a n e t o t e l l u r i cd a t ap r o c e s s i n . m u l t i s i t e q yt p g g ,L ,C 犓 犲 狑 狅 狉 犱 狊 M a n e t o t e l l u r i c o c a ld i s t o r t i o n o n u a t ei m e d a n c et r a n s f o r m a t i o n, T e n s o r g j g p 狔 , d e c o m o s i t i o n D i m e n s i o n a l i t n a l s i s p ya y 轴方向 , 根据其构造 维 性 参 数 也 不 能 获 得 正 确 的 区
相关文档
最新文档