可控源音频大地电磁三维共轭梯度反演研究

可控源音频大地电磁法在资源勘探等领域中发挥着重要的作用.我们把有限差分数值模拟方法用于可控源音频大地电磁三维正演,结合正则化反演方案和共轭梯度反演的思路,将反演中的雅可比矩阵计算问题转为求解两次"拟正演"问题,得到模型参数的更新步长,形成反演迭代,实现了可控源音频大地电磁三维共轭梯度反演算法.该反演算法可用于对有限长度电偶源激发下采集到的可控源音频大地电磁全区(近区、过渡区和远区)视电阻率和相位资料进行三维反演定量解释,获得地下三维模型的电阻率结构.理论模型合成数据的反演算例验证了所实现的可控源音频大地电磁三维共轭梯度反演算法的有效性和稳定性.     

CSAMT单分量数据解释方法

可控源音频大地电磁法(CSAMT)一直沿用大地电磁法(MT)的办法,通过计算电场分量与磁场分量的比值,求取卡尼亚视电阻率.而CSAMT场源已知,电场分量和磁场分量都与地下电阻率存在一定的关系,可以单独采用CSAMT电场分量或者磁场分量提取地下介质的视电阻率.本文通过分析电场分量与磁场分量的数据特性,提出利用CSAMT电场单分量数据进行视电阻率的计算,用改进的广义逆矩阵反演方法,使初始模型中的地电层数等于频道个数,克服了以往反演计算中层数较少的问题;实现全场区电场分量视电阻率曲线的拟合反演.同时对单分量视相位计算方法进行分析,结合山西大同地区积水采空区探测及数据解释结果,论证本文提出的单分量解释方法的有效性.     

利用CSAMT电磁场响应提取激电效应的方法初探

在传统的CSAMT电磁场一雏正演中加入柯尔-柯尔模型,模拟可极化地电断面的CSAMT电磁场响应.通过对均匀半空间、H型地电断面的理论计算,分析了柯尔-柯尔模型参数对CSAMT的Ex、Hy影响特征.基于电场相对磁场对激电效应的反映更加灵敏,提出分别利用Hy与Ex分量反演电阻率、厚度以及激电参数的方法.反演计算表明,该方法克服了电阻率、厚度、激电参数一同反演时受初始模型影响太大,容易陷入局部极小值的问题,可以在一定程度上提取出合适的激电参数.     

三维电阻率测深数据Zohdy近似反演方法

Zohdy方法近似反演三维电阻率测深数据。正演计算采用有限单元法。反演初始模型由测量视电阻率数据给出。通过比较实测视电阻率值和预测模型计算的视电阻率值对数差来修改模型网格电阻率．为了解决任意电极距测深数据的反演,采用大、小双网格剖分。大网格反映地下电性分布情况。小网格用于实际有限元正演计算．在电阻率调整公式中加入一个迭代系数,能够加快收敛速度．并对加5％随机噪声的模型理论视电阻率测深断面数据进行反演,得到的电阻率分布与模型电阻率基本一致．     

尖锐边界反演在临江市六道沟可控源音频大地电磁测深数据反演中的应用

在实际大地电磁测深数据反演过程中,为了得到更加清晰的电性单元分界面,传统的Occam平滑反演方法遇到了挑战,本文讨论了sharp boundary inversion(SBI)反演的理论及其应用.通过模型试算,对比Occam和SBI的反演效果,证明了SBI反演划分层状介质的可靠性.本文对临江市六道沟地区实测可控源音频大地电磁测深(CSAMT)数据进行了全区视电阻率转换,然后对全区视电阻率进行Occam反演和SBI反演.最后对比已有地质资料,表明SBI反演能够用于CSAMT数据的处理中.     

高密度电阻率法逐级反演与实践——以湖南省常德市鼎城区某地质剖面为例

为提高高密度电阻率法的探测精度和效果,在常德市鼎城开发区开展了基于比照印证已有地质钻孔数据进行约束模型的一维到二维的逐级反演试验.文中首先阐述了一维层状模型的正反演解析解的计算方法,并完成了高密度电阻率法二维模型的正演响应公式和基于高斯牛顿法的反演理论推导;其次设计层状模型利用有限单元法进行了正演模拟,从其响应数据中抽取了单点测深数据进行层状反演,并利用一维反演结果结合已有地层电阻率信息构建了二维初始模型进行反演,实现了高密度电阻率法逐级反演的理论研究;最后,选取常德市鼎城区跨越常德—益阳区域构造的高密度电法剖面进行试验,并根据地质情况和钻孔信息对野外采集的连续剖面的每个测点的数据进行了层状划分,构建二维模型对其进行反演,完成了高密度电阻率法逐级反演的应用研究.结果表明:在城市开发区利用逐级反演的思想对高密度电阻率法探测数据进行处理,有助于提高反演精度和地质解释的准确度.     

线源时间域电测深全期全区视电阻率求解

由于线源一维的正演国内外研究尚少,文中主要利用水平电场分量求取全区视电阻率。通过计算发现水平电场与电阻率之间是一一对应的关系,这使得定义全区视电阻率具有惟一性和可靠性。文中以水平电场定义视电阻率,然后利用二分搜索算法进行全区视电阻率求解。模型计算表明,该方法定义的全区视电阻率能较准确地反映地下电阻率信息。进一步研究发现:对时间域线源而言,在近区算不好晚期的数据,因为电磁波的传播、衰减很快,同样在极远区算不好早期的数据,而且探测深度是随收发距的增加而增大。最后,加入5%的随机噪声后与理论结果进行对比,得到了收发距与最大探测深度的关系。     

海洋可控源电磁勘探中MT噪声降噪效果评估方法研究

信号去噪是海洋可控源电磁勘探(Marine Controlled-Source Electromagnetic,MCSEM)数据处理过程中的重要环节.在中、远收发距处,大地天然电磁信号(Magneto Telluric,MT)是影响MCSEM反演质量的主要噪声之一.尽管已有学者提出一些MT噪声降噪方法,但由于缺乏降噪效果的量化评价,导致这些方法的有效性亟待验证.针对上述问题,本文将信号和噪声一维正演、噪声降噪处理、电阻率一维反演相结合,提出一种MT噪声降噪效果评估方法.该方法利用电阻率反演结果与电阻率输入值的均方相对误差作为降噪效果的量化评价指标.为提高灵活性,既可以利用正演的手段生成MT噪声,也可以直接输入MCSEM发射机不开机时间段的实测数据作为MT噪声.利用该评估方法,验证了短时窗结合预白化滤波降噪方法的有效性.在所开展的实验中,实测数据的均方相对误差为231%.相比于不进行降噪处理的直接反演结果,均方相对误差减小了126%.本文所研究的评估方法可应用于MT噪声降噪方法有效性的量化评估,也可推广至其他类型噪声降噪效果的量化评估.     

地震资料约束下大地电磁资料反演

应用地震资料约束下大地电磁资料反演古潜山或基底内部结构的新方法. 首先利用地震及井资料,反演得到古潜山或基底之上各层的厚度及电阻率,建立了古潜山或基底以上的地电模型;并利用大地电磁一维反演资料,建立了古潜山或基底及其以下的地电模型. 在初始模型建立的基础上,应用高效快速的二维大地电磁正演算法——模式匹配方法,进行正演计算. 利用共扼梯度迭代方法,进行正反演迭代计算. 通过对胜利油田南北618剖面反演, 基底以下在地震资料没有显示的构造信息,在反演结果中得到了较好地反映.      

复杂介质有限元法2.5维可控源音频大地电磁法数值模拟

利用有限元2.5维可控源音频大地电磁法（简称CSAMT）数值模拟方法,对100Ωm均匀半空间介质中有限长度的电偶极源产生的电场、磁场及视电阻率、相位特征进行了数值模拟,研究了场的空间变化规律. 在一个象限中,场的特征存在双叶现象,当收发距大于4个趋肤深度时,电阻率较接近介质真实的电阻率,这些结果为观测系统和收发距的选择提供了依据. 波数域场的特征表明,低波数对源的贡献占较大的比例. 有限元法2.5维CSAMT数值模拟的优势在于能较准确地获得复杂介质结构的波场特征. 本文结合直立异常体、倾斜异常体及断陷模型对CSAMT电阻率、相位剖面特征及频率曲线特征的可靠性进行了研究. 数值模拟结果直观地给出了异常体的剖面异常形态. 通过对比研究异常体的剖面异常形态和半空间场的特征进一步说明本文方法和软件在模拟复杂介质结构场特征时是可靠的. 这为认识观测数据,指导反演解释提供了较好的依据.     

基于主轴各向异性介质的三维张量CSAMT正演模拟

在地质情况复杂化以及地质勘探精细化的要求下,有必要开展张量CSAMT正反演研究.考虑到各向异性对正反演结果影响较大,本文开展基于主轴各向异性介质的张量CSAMT正演研究.根据含源的电场控制方程,引入主轴各向异性介质的张量电阻率,将电场总场分解为一次场和二次场,一次场由Key的算法得到,二次场利用有限差分法计算,即实现了主轴各向异性介质的张量CSAMT三维正演.通过与二维有限元算法结果比较来验证算法的正确性;通过设计几组主轴各向异性三维模型,分析卡尼亚电阻率和相位的响应特征.结果表明卡尼亚电阻率ρxy和ρyx分别对异常体X方向和Y方向的电阻率变化敏感,说明张量CSAMT方法对识别各向异性具有优势.本研究为开展任意各向异性介质的张量CSAMT正反演研究提供基础.     

具有干预机制的CSAMT数据一维最优化反演

在CSAMT测深中,测深曲线往往受到电偶极源的影响而出现畸变,这给CSAMT数据的反演解释带来很大的困惑。文中依据CSAMT法在波区与MT理论相似的原理,提出用Bostick半定量结果作为CSAMT一维反演的初始模型,采用传统最小二乘算法,结合多种常识性物性变化特点,优化层参数并控制反演过程的迭代方向,使纯粹的数值迭代转化为按地球物理特征自动迭代,其结果可更好地满足物探解释需要。算例和实测数据反演表明,优化后的CSAMT一维反演方法,反演精度较高、计算速度较快,其反演结果的正演响应曲线和实测曲线拟合得很好。     

基于牛顿迭代法和遗传算法的CSAMT近场校正

由于可控源音频大地电磁法（CSAMT）采用人工场源在大大增加信号强度的同时也带来了在近区产生非平面波效应的问题,其表现特征为在双对数坐标中卡尼亚电阻率曲线呈45°上升,即所谓的近场效应.本文首先提出了基于牛顿迭代法求解视电阻率方程的近场效应校正方法,通过对二层、三层理论模型试算验证了该方法的有效性;为进一步增强对噪声的压制能力,本文进一步引入了遗传反演的思想,将全频段误差最小作为目标函数,采用自适应正则化方法引入模型参数梯度最小作为稳定因子,由遗传算法求解得到校正后视电阻率,通过对理论模拟数据加随机噪声的校正结果,表明该方法在实现近场校正的同时能剔除噪声干扰,得到的结果与理论模拟曲线吻合度较高.通过对内蒙古曹四夭钼矿区实测CSAMT数据处理结果表明,牛顿法和遗传算法均能明显校正CSAMT因近场效应引起的假高阻异常,能更好地反映地电结构特征.     

可控源音频大地电磁数据的反演方法

从反演方程、构造目标函数和求解三方面对用于可控源音频大地电磁法(CSAMT)的实用反演方法中的四种进行了描述．水平层状地层CSAMT法资料的直接反演法首次尝试了一维空间的全资料CSAMT反演，效果较好，但该方法尚难应用于2D、3D复杂介质中；奥克姆反演方法既考虑了横向的光滑函数，又考虑了纵向的光滑函数，得到比较光滑的横向、纵向变化的背景电性结果，但有可能把一些小构造光滑掉．快速松驰反演算法和共轭梯度算法由于计算速度快，占内存少而被用于三维反演中，二者相比，快速松驰算法在求解雅可比矩阵时只做一次正演计算，在更新模型时解小型方程组，所以在速度上更胜一筹．在后三种算法中，由于复杂电性结构无解析解，正演计算都采用数值计算．数值计算的可靠性、速度影响着反演算法的有效性，这方面的研究也将是2D、3D复杂电性结构反演的研究方向之一.     

基于矢量有限元的带地形大地电磁三维反演研究

研究了基于矢量有限元方法的大地电磁带地形三维反演算法并开发了三维反演计算程序代码.在大地电磁场正演数值模拟方面,采用并行直接稀疏求解器PARDISO且无需进行散度校正的快速正演方案,对典型地形模型,在中等规模计算条件下,与双共轭梯度法(BICG)计算结果比较,发现PARDISO比BICG快10倍以上;通过理论模型试算,并与前人的有限元法计算结果对比,验证了带地形三维正演计算程序的正确性.在反演方面,本研究基于共轭梯度方法编写了大地电磁带地形三维反演代码,为了避免直接求取雅可比矩阵,将反演中的雅可比矩阵计算问题转为求解两次“拟正演”问题,进而将PARDISO的快速正演方案应用于“拟正演”问题的求解,以提高反演计算效率.利用开发的反演算法对多个带地形地电模型的合成数据进行了三维反演,反演结果能很好地重现理论模型的电性结构,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了某矿区大地电磁实测数据,反演得到的三维电性结构清晰地反映了研究区的地电特征,将反演结果与该区已有地质资料结合进行解释,应用效果明显,进一步验证了本文算法的有效性.     

CSAMT三维正演数值模拟研究进展

可控源大地电磁(CSAMT)资料的三维正、反演问题已成为国际地球内部电磁感应领域研究的前沿课题.本文介绍了可控源音频大地电磁法中三种主要的数值模拟正演方法,即有限元法、有限差分法,积分方程法,比较了这三种方法的优缺点.积分方程法只需对异常区进行剖分,仅需计算小体积异常区的场,计算速度快,但只适合模拟简单模型,精确度也较低;有限元法与有限差分法虽然精确度较高,但要求对全部区进行离散化,占用的计算机容量较大,计算时间长,因此基于并行算法的三维电磁场正演研究可能是可控源音频大地电磁法未来的发展趋势.     

浅部频率域电磁勘探方法综述

适用于近地表(2000m以内)勘探的频率域电磁法主要有音频大地电磁法(audio-frequency magnetotellurics,AMT),无线电大地电磁法(radio-magnetotellurics,RMT),可控源音频大地电磁法(controlled source audiofrequency magnetotellurics,CSAMT),广域电磁法(Wide Field Electromagnetic Method,WFEM).本文拟从最新的数据采集技术、数据处理技术、正反演算法、实例等四个方面,论述适用于浅部勘探的AMT,RMT,CSAMT和WFEM方法的国内外最新进展,总结目前AMT,RMT,CSAMT和WFEM方法遇到的困难,对潜在的发展方向提出建议.综述表明:(1)张量测量、多站阵列、多站叠加可提高AMT、RMT和CSAMT数据的质量.利用近区数据WFEM法可获得良好的效果.国产与国外仪器在质量方面的差距正在逐步缩小.(2)数学形态滤波技术、Hilbert-Huang变换等可有效分离出有用的数据,局部畸变仍然是亟待解决的难题,需要更为深入的研究.(3)矢量有限元与非结构网格的出现大幅度提高了有限元处理复杂电磁问题模拟的精度与应用范围,成为目前电磁正演的首选工具.完全非线性反演算法仍然局限于1D、2D问题,共轭梯度法和高斯牛顿算法等为解决3D问题的发展趋势.地质约束的引入和多数据联合反演可以减小反演的非唯一性.各向异性的反演为目前反演研究的热点之一.(4)野外数据解释的正确性严重依赖于对地下结构先期的维性判别,在2D特性不明显、3D特性明显时,需要采用3D进行反演解释.     

2.5维可控源音频大地电磁法Occam反演理论及应用

Occam反演是一种正则化的光滑模型反演方法,它在寻找最小拟合差的同时追求最光滑模型.因地下介质的电性通常是连续变化的,Occam反演考虑了地下介质横向和纵向的光滑情况,应用Occam反演避免了模型电性参数被随意间断.通过建立层状和存在两个异常体的模型,利用2.5维可控源音频大地电磁有限元法作正演响应计算,对响应数据基...     

起伏地形下CSAMT二维正反演研究与应用

CSAMT在山区金属矿勘查中,采用各种滤波和相位积分之类的处理方法,校正因地形起伏和局部电性不均匀引起的静态效应,往往难保奏效,开发消除静态效应的新方法是提高CSAMT资料处理与解释水平和方法应用效果的重要研究课题. 本文以如何消除地形影响为重点,对起伏地形下CSAMT二维大地三维源地电模型,采用加权余弦数值积分法,进行波数域电磁场二维有限单元法正演. 为模拟复杂地形地电模型,选取交叉对称网格三角形剖分法,实现了在国内常用赤道电偶极装置的CSAMT二维正演计算;在二维正演的基础上,开发了基于奥克姆反演法的CSAMT二维反演技术,研制出一套起伏地形下CSAMT二维正反演处理与解释方法技术系统. 通过理论模型试算和实测数据处理证实,本系统能有效地削减起伏地形影响. 在找矿应用中,该系统反演的电阻率断面,极大地消除了起伏地形影响和静态效应,突显出清晰的控矿构造和矿体的异常,取得了重要成效.     

频率域双极源全区视电阻率的计算及分析

频率域可控源电磁法由于其发射源的非偶极性,收发距受到发射场源功率的限制,在近区和过渡区,场源不能看做偶极子源,而是双极源.对于双极源数据在近区和过渡区的校正,若采用偶极予源校正,存在很大的误差,收发距一定时,偶极距越大,误差越大.双极源全区视电阻率利用半空间模型视电阻率等效地下复杂模型视电阻率计算全区视电阻率,对比偶极...