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

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

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

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

Occam反演在西准噶尔包古图地区某铜矿AMT资料处理中的应用

Occam反演方法经常应用在音频大地电磁测深数据中.为了了解Occam反演方法在金属矿的使用效果,本文设计了层状和脉状异常体的理论模型,利用音频大地电磁有限单元法做正演响应计算,并对加高斯噪声后的理论数据使用Occam反演进行反演.实验表明,Occam反演对音频大地电磁法的理论数据反演具有稳定收敛性,反演结果和理论模型能够很好地吻合.并以实例验证了Occam反演的有效性,为以后AMT方法寻找金属矿产提供了理论依据和实际经验.     

基于散射理论与动态时间规整的时移地震波动方程差异反演

时移地震勘探作为储层动态监测的有效手段,对获取地下介质的弹性或物性参数变化起到了至关重要的作用.由于多期地震勘探之间介质的传播速度存在差异,同时环境条件、采集方式及处理手段等方面的不同,都会使地震记录产生走时差异.而这种走时差异对时移地震反演存在着严重的影响,尤其是基于波动方程的反演方法.为了消除时移差异,我们引入动态时间规整(DTW)方法实现角道集优化匹配,以获取地震记录的有效变化信息.此外,本文提出根据散射理论,通过弹性波动方程进行时移地震差异反演,得到地层的差异弹性参数,实现储层差异定量表征.该方法采用WKBJ近似模拟非均匀背景,具有更高的反演精度,尤其是在复杂地质构造条件下优势更加明显.本文对所提出的方法进行了数值模拟实验与实际数据应用,证明该方法得到的反演结果优于传统方法,可以实现对地层时移差异的准确刻画.     

多种频率域电磁法在冷水坑矿集区的应用效果对比

江西省冷水坑矿集区具有浅部斑岩型矿体与深部层控叠生型矿体的复杂组合.利用现代地球物理方法对矿集区地下深部结构进行了探测,对圈定地下容矿、示矿地质信息和认识地下深部成矿过程具有重要意义.本文以位于江西省贵溪市境内的冷水坑银铅锌矿田为试验研究区,开展了混场源音频大地电磁测深法(EH4)、音频大地电磁测深法(AMT)以及可控源音频大地电磁测深法(CSAMT)三种不同频率域的探测试验,研究和总结了频率域方法以及不同仪器电磁法数据采集特点和应用效果,为频率域电磁法勘探提供借鉴和参考.通过对冷水坑矿集区典型勘探线上三种电磁方法的数据处理,绘制了视电阻率和相位单点测深曲线、原始拟断面及二维反演断面图.发现单点曲线形态基本一样.由于AMT在低频段受矿区工业干扰,视电阻率曲线呈45°上升;视电阻率、相位拟断面对地下电性分层结构基本一样,总体可分为三层,且分层结构较清晰;二维反演断面能够较好反映地下高低阻异常响应,但响应效果有所差异且CSAMT法效果较好.因此在干扰严重地区,建议选用人工源发射的可控源音频大地电磁测深法(CSAMT).     

可控源音频大地电磁测深法勘查深度研究

可控源音频大地电磁测深法的勘查深度问题在理论研究和实际工作中都得到了广泛的关注,本文分析了水平层状大地电磁正反演方法的基本原理,说明了可控源音频大地电磁测深法在深部勘查方面具有可靠性.通过建立两个三层地电模型,并对吉林某地区的CSAMT实测数据进行了带源反演成像,与已知钻井资料对比分析表明,可控源音频大地电磁测深法对地下1000~2000 m深度的反演较为准确,能够反映地下电性的垂向变化.     

海岸效应对近海地区大地电磁测深数据畸变作用研究

在近海地区采集的大地电磁测深数据通常受到海岸效应的影响,使得大地电磁测深数据发生畸变,因而很难利用大地电磁测深资料较为可靠地获得地下深部的电性结构.本文通过正演模拟方法,分析和总结海水深度变化和海底地形变化对近海地区大地电磁测深数据的畸变影响.当测区与海岸线的距离小于目标频率的大地电磁场趋肤深度时,高导海洋的存在会严重影响测区内电磁场的分布.由于海岸效应的影响,大地电磁测深视电阻率曲线和相位曲线均会发生不同程度的畸变,在低频部分,这种畸变作用尤为明显.大地电磁测深一维Occam反演方法和二维非线性共轭梯度反演方法,对近海地区浅部地层具有较好的反演效果.随着海水深度的增加和海底地形的复杂变化,两种反演方法均会出现不同程度的假异常,为地质解释工作造成了影响.近渤海地区的实测大地电磁测深数据在低频部分可能受到海岸效应的影响而导致视电阻率曲线的严重畸变.     

三维主轴各向异性介质中张量CSAMT正反演研究（英文）

相较于常用的标量CSAMT而言,由于使用多发射源,张量可控源音频大地电磁(CSAMT)可以接收到更加丰富的电场与磁场信息。然而,目前大多数关于张量CSAMT的理论、数值模拟算法以及反演研究都是基于远区的测量和地下介质为电阻率各向同性的假设。本文采用三维交错网格有限差分数值模拟算法,研究电阻率主轴各向异性介质与各向同性介质的响应差异,并在此基础上采用有限内存拟牛顿法(LBFGS)实现了三维张量CSAMT主轴各向异性反演。为加速和稳定LBFGS反演算法,本文采用了模型参数转换策略并修改了LBFGS第一次迭代的初始步长与线搜索停止条件。反演试算结果表明,当地质体存在主轴各向异性特征时,三维主轴各向异性的反演能在一定程度上较好的反演出地下构造,而不考虑各向异性效应,直接使用各向同性的三维张量CSAMT反演程序,会使最终的反演结果出现严重的假异常。     

时间域航空电磁法激电效应对电磁扩散的影响

由于激发极化效应的影响,时间域航空电磁晚期道信号经常会出现变号现象.基于电阻率的传统反演方法无法对变号数据进行正确反演,因此通常在数据处理中予以剔除.为深入了解极化介质的电磁扩散特征,认识航空瞬变电磁负响应的产生机理,本文研究时间域航空电磁系统的电磁扩散特征.我们以均匀极化、非极化半空间及层状介质模型为例,通过直接积分的方法求解频率域电场响应,并由欧姆定律得到电流响应,再经过汉克尔变换得到时间域电流响应.通过研究电流随时间在地下极化介质中的传播特征研究电磁扩散过程;通过对比不同激电参数对电磁扩散的影响,研究极化介质中感应电流与极化电流的扩散规律,从而合理地解释极化介质中负响应的产生机理.基于本文研究和分析结果,可加深对时间域航空电磁法中激电效应的认识.     

基于瞬变电磁矩变换的快速三维反演方法

瞬变电磁法的严格三维反演计算复杂、占用资源多,在普通计算机上难以实现.本文引入瞬变电磁矩变换的概念,提出一种快速三维反演方法.该方法基于阻性限制(resistive limit)特性,建立包含异常体的三维大地的一阶矩响应正演算法,根据不同约束条件,选择优化的最速下降法实现瞬变电磁快速三维反演.文中通过含异常体的三维大地正演一阶矩与仿真数据一阶矩的对比,验证了快速三维正演算法的有效性,之后在不同约束条件下,利用优化的最速下降法实现了对含噪声的仿真瞬变电磁数据的快速三维反演.结果表明,该方法能够在普通计算机上短时间内较为准确地反演出地下异常体的体积和位置,在瞬变电磁数据的实时解释工作中具有良好的应用前景.     

沁水盆地南部煤层气藏三维音频大地电磁探测

音频大地电磁法以其快速、经济、地形适应性强等优势,在复杂环境下煤层气藏勘探方面具有巨大应用潜力.本文分析了沁水盆地南部煤层气成藏模式,设计了适合音频大地电磁法勘探尺度的三维薄低阻、高阻目标层模型.模型典型测线的二维反演以及三维阻抗张量反演结果表明:在浅表不均匀体存在时,二维TE、TM、TETM模式反演解释会受到三维畸变不同程度的影响,很难判断哪种模式解释最佳.相比而言,三维反演不仅可以提取浅表异常,而且能以更高的分辨率圈出深部薄低阻层分布.但是上述反演手段对薄高阻目标层提取效果不佳,只有对初始模型(甚至目标层)施加电性约束时,薄高阻层反演结果才会改善.综合考虑,音频大地电磁法对薄低阻目标层进行三维探测是最具可行性的.选取沁水盆地南部胡底区块作为研究区,通过音频大地电磁三维反演提取富水低阻目标层分布.推断煤层气成藏模式,识别煤储层并推测储层排采状态,与已有的地质、勘探和生产资料基本吻合.以上分析为音频大地电磁三维探测法在油气藏勘探与动态监测方面的应用提供了有力支撑.     

基于各向异性岩石物理的缝隙流体因子AVAZ反演

裂缝型储层表现出较强的各向异性特征.缝隙中充填不同流体时,裂缝储层的地震响应特征也不相同.本文从各向异性岩石物理模型出发,引入可有效识别缝隙流体的指示因子,并研究缝隙充填流体类型、饱和度以及缝隙纵横比与流体因子的相互关系,进而分析不同流体充填时介质的地震响应特征,并基于AVAZ反演方法估测缝隙流体指示因子.首先对缝隙流体因子的敏感性进行了分析,讨论当缝隙充填不同流体时,缝隙流体因子值的变化特征,同时研究了不同流体类型充填时裂缝储层反射系数随方位角和入射角的变化特征.某工区测井数据和复杂裂缝模型应用表明,基于各向异性岩石物理的缝隙流体因子AVAZ反演方法合理、可靠,且具有良好的抗噪性,即当对合成地震记录添加信噪比不小于1/2的随机噪声时,利用AVAZ反演方法估测所得流体因子值与真实值仍然吻合较好.     

不同参数条件下远参考大地电磁处理效果对比分析

研究了在不同噪声强度、不同参考道、噪声仅存在基站或参考站时,远参考处理对测深曲线的影响.通过对标准时间序列进行仿真模拟,发现基站含不同信噪比的非相关噪声,远参考大地电磁对测深曲线提升效果不一:高信噪比(SNR=10.58 d B)条件下,远参考大地电磁能够较完美的恢复测深曲线,低信噪比(SNR=0.2 d B)条件下,远参考大地电磁对测深曲线质量提升效果有限;相同噪声仅存在于基站或者参考站对资料处理效果有区别:噪声仅存在于参考站的处理效果比仅存在于基站的效果要好;参考站无噪声情况下,电道或者磁道作为参考道对处理结果影响不大;基站与参考站之间的磁场相关性、基站磁场噪声与参考站磁场的非相关性与远参考处理效果关系密切.通过对实际资料的计算也有上述结论.     

有效压力与饱和度变化的气藏多波时移地震AVO响应

天然气在开发过程中,储层有效压力和含气饱和度均会发生变化,研究有效压力和含气饱和度的变化对地震响应特征的影响,在基于时移地震的剩余气分布预测研究中具有重要意义。天然气和石油的声学性质有着明显的差异,油藏时移地震的研究成果不能直接应用于气藏,因此需要开展气藏的时移地震研究。利用Shapiro模型表征干岩石弹性模量随有效压力的变化,借助Batzle-Wang方程描述流体速度随压力的变化关系,联合Gassmann理论进行流体替代,表征饱和流体岩石速度随含气饱和度的变化,建立了饱和流体岩石速度随有效压力和饱和度变化的岩石物理模型。基于该模型,对不同含气饱和度和不同有效压力下的气藏储层模型进行了多波时移地震叠前振幅变化（AVO）模拟。结果表明多波时移地震AVO技术可以有效地区分有效压力变化和含气饱和度变化,为进一步开展气藏多波时移地震流体监测提供了理论参考依据。      

正则化等效层重力向下延拓方法

近年来,利用时移微重力技术进行储层开发监测受到国内外学者广泛关注.时移微重力观测数据存在信噪比低,信号弱的问题,难以实现储层内物质运移的定量解释.为压制数据噪声,增强有效弱信号,本文研究了利用Tikhonov正则化方法反演等效层（源）,并由等效源实现重力场向下延拓的方法;在此基础上,本文推导了波数域正则化等效源向下延拓算子.针对向下延拓场幅值衰减问题,提出了正则化等效源迭代补偿算法.通过模拟数据实验研究了不同深度正则化等效源滤波算子及向下延拓算子的波数响应;与波数域Tikhonov正则化向下延拓方法相比,正则化等效源向下延拓方法的延拓精度更高、更稳定.最后,将基于迭代补偿的正则化等效源向下延拓技术应用于实测时移微重力数据证实了该方法能够有效增强局部异常,实现时移微重力数据大深度稳定向下延拓.     

面向目标自适应有限元法的带地形三维大地电磁各向异性正演模拟

传统三维大地电磁各向异性模拟均是基于规则六面体网格,计算精度有限且较难拟合复杂地质条件.本文采用面向目标自适应非结构矢量有限元法,对三维大地电磁各向异性介质进行模拟.首先从电场双旋度方程出发,利用伽辽金方法建立变分方程;然后利用电流密度连续性条件构建适合大地电磁各向异性问题的加权后验误差估计方法,实现面向目标的网格自适应正演;最后通过典型算例分析各向异性对网格自适应和大地电磁响应的影响特征以及各向异性的识别方法.本文算法能够高精度地拟合起伏地表和任意各向异性介质,适用于分析复杂地电条件大地电磁响应特征,为提高大地电磁资料解释水平提供了理论基础.     

大地电磁资料处理的计算方法

在地球物理勘探资料处理中,压制噪声、提高信噪比是中心环节。在大地电磁估算张量阻抗时,噪声的大小决定了张量阻抗的精度.而张量阻抗精度又直接影响了反演质量(即推断地下电性结构的准确性)。噪声可以分为两种类型,一类是相关噪声,第二类是不相关噪声。所谓相关与不相关是指噪声与信号之间、噪声与噪声之间的相关与不     

大地电磁与地震联合反演研究现状与展望

由于大地电磁探测和地震勘探各自存在优势和局限,因此联合反演比单独一种地球物理资料反演更优越.本文剖析了联合反演的两类分类方法,指出大地电磁与地震方法具备联合反演的基础-理论上相似,物性上有联系,并重点分析国内外大地电磁与地震数据联合反演的研究现状与存在问题,总结层状均匀和非均匀两种介质结构的电震联合反演算法.对于如何耦合不同类型的物性参数这一联合反演的核心问题,指出目前广泛使用的岩石物理法和结构法两类方法的适用范围.最后提出大地电磁与地震联合反演的发展方向:层状均匀介质结构要优化线性迭代技术;非均匀的复杂地质结构,需加紧非线性反演理论向非线性反演发展,其中结构法反演前景广阔,要加强实际资料的运用;此外,电震联合反演也需与岩石物性等先验地质信息相结合.     

多通道瞬变电磁法激励源降噪性能对比研究

对比研究双极性方波与m序列两种激励源在多通道瞬变电磁法(Multi-Transient Electromagnetic,MTEM)中的降噪性能,将大地脉冲响应的提取等效为辨识滤波器输出,通过分析辨识滤波器的谱特性量化评估不同激励源的降噪性能;通过数值模拟讨论两种激励源信号在工频噪声与白噪声环境下的大地脉冲响应辨识精度,模拟结果表明使用m序列替代双极性方波作为激励源能够大幅提高辨识精度,且残留噪声集中在码元频率及其谐频处;设计野外对比试验,采集两种激励源在弱/强干扰下的响应信号进行大地脉冲响应辨识,试验结果表明弱干扰环境下两种激励源的辨识精度相近,强干扰环境下m序列的降噪性能高于双极性方波,且噪声残留特征与数值模拟吻合.     

AMT 、TEM 、VES地层响应特征模拟分析及其联合反演探讨

通过正演模拟方式分析音频大地电磁法(AMT)、中心回线瞬变电磁法(TEM)、直流电测深法(VES)对地层的响应特征,结果表明:TEM、AMT对地层的响应特征类似,即对低阻层反应灵敏,而对高阻层的分辩能力较差;VES对高阻层及低阻层的分辩能力相当,相对于TEM、AMT,VES对低阻层的分辩能力差,而对高阻层的分辨能力较高.基于AMT或TEM与VES对地层的响应具有互补的特点,探讨了对AMT或TEM与VES数据进行的联合反演.结果表明,AMT或TEM与VES进行联合反演,相比单一方法进行反演,结果能得到明显的改善.