从CSAMT信号中提取IP信息探讨

CSAMT信号数据进行反演时,认为电阻率是一个与频率无关的实数,而实际上因为激电效应,地下可极化体的电阻率是一个与频率有关的复数.从CSAMT信号中提取激电信息可以加大激电法的使用范围,并可以提高CSAMT法勘探的精度.本文总结了国内外在这方面研究现状,介绍了从CSAMT法信号中提取激电信息的基本原理,并提出了未来值得深入研究的几个方面.     

三维大地电磁激电效应特征研究

在积分方程法的大地电磁三维正演模拟中引入Cole-Cole模型研究激电效应特征,对均匀半空间中存在极化体时的三维大地电磁测深响应进行了理论计算,分析了极化参数的影响规律.结果表明:激电效应使观测视电阻率值变低;地下高阻极化体比低阻极化体的激电效应强,特别是极化率大的高阻极化体对大地电磁响应影响大,使观测结果出现假异常.     

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

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

时间域航空电磁激发极化参数三维反演研究

时间域航空电磁中心回线(或重叠回线)装置晚期道数据受激电效应影响常出现符号反转现象.这类数据与多个激电参数相关,并且各参数之间灵敏度差异较大,导致反演存在严重的非唯一性.本文提出一种基于Pearson相关性约束和深度学习算法相结合的时间域航空电磁激发极化参数反演策略.该反演策略首先基于深度学习预测时间域航空电磁激电参数,进而给时间常数和频率相关系数一个较小的约束范围后再反演电阻率和极化率,由此大大减少反演的多解性.针对电阻率和极化率的反演,我们采用统计学中Pearson相关系数构建两种物性参数的相关性约束,进一步减少反演多解性.为验证反演策略的有效性,我们对双棱柱模型和拱形模型分别进行反演试算.理论测试结果表明,基于Pearson相关性约束的电阻率和极化率的反演结果比传统的高斯-牛顿反演结果更接近真实模型,而基于深度学习预测时间常数和频率相关系数后的电阻率和极化率反演结果与给定真实时间常数和频率相关系数后的反演结果效果相当.最后,我们对来自澳大利亚的带激电效应的航空电磁实测数据在考虑和不考虑激电效应条件下进行反演,结果表明考虑激电效应的反演无论数据拟合还是地电断面的连续性均得到明显改善...     

基于Cole-Cole模型的中间极化水平层大地电磁IP效应研究

利用Cole-Cole模型的复电阻率公式表征的介质激电效应,对一维层状介质O型、H型、A型、K型和Q型地电断面中间极化的大地电磁模型进行了正演计算.利用视电阻率比值和相位比值参数来表征IP效应,详细讨论了各地电类型的激发极化效应.对大地电磁场激发极化效应与极化体厚度和埋深的关系进行了讨论.     

基于自适应有限元的2.5D海洋可控源电磁场激电效应特征研究

油气藏自身会产生激电效应,明确海洋可控源电磁场激电效应特征非常重要。本文基于自适应有限元方法实现了带激电效应的2.5D海洋可控源频率域电磁场正演。采用非结构化三角网格单元剖分地电模型,可以模拟任意起伏地形和复杂地电构造。利用Cole-Cole模型表述油气藏的激电效应,实现了考虑激电效应的海洋电磁复电阻率正演算法。将本文的计算结果与已发表的一维模型结果对比,检验了本文算法的正确性。最后,设计了不同的极化模型,针对不同的场源方位角、激电参数和地形,计算了的电磁场响应结果,分析了激电效应对海洋可控源电磁场的影响规律,结果表明油气藏的激电效应和海底地形对海洋可控源电磁场有一定的影响。     

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

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

三维复电阻率模型电磁场正演模拟研究

复电阻率法是上世纪70年代发展起来的一种激电分支方法,自提出伊始就受到了国内外学者的广泛关注.很多学者开展了复电阻率电磁场正演计算,但缺少对正演结果的深入分析.本文采用体积分方程法,进行了复电阻率模型的正演模拟,通过对正演结果分析,得出三维复电阻率体的电磁响应是由累积电荷,电磁感应和激电效应共同产生的,并总结了三种效应的电磁响应规律,为准确分析异常提供了依据.     

基于油气藏两种极化模式的海洋可控源电磁法激电效应研究

海洋可控源电磁法(MCSEM)已日益成为海洋油气勘探的重要方法,但在传统MCSEM的数值模拟和反演解释中,通常考虑的是非极化海底介质模型,这与油气藏会产生激电效应的实际情况并不相符.在前人的研究基础上,本文总结了油气藏产生激电效应的两种模式:外生黄铁矿极化模式和油气藏自身极化模式.在典型的海洋4层地电模型基础上,本文设计了包含外生黄铁矿的海洋6层地电模型.基于这两种模型,通过推导MCSEM一维正演公式,引入Cole-Cole模型,本文分别计算了上述两种极化模式下MCSEM的激电响应,分析了各参数对两种模式下激电响应的影响特征,并对两种极化模式共同作用所产生的激电效应特征做了分析.结果表明,外生黄铁矿极化和油气层自身极化都会产生较明显的激电效应,而激电效应的强弱与观测频率、模型参数和激电参数有关.在相同的激电参数下,外生黄铁矿极化所产生激电效应的影响主要集中在小收发距处,而油气层自身极化所产生的激电效应则主要影响较大收发距处的场值.两种极化模式共同作用的效果相当于两种极化模式分别作用时的叠加.     

编码电磁测深

利用逆重复M序列伪随机信号良好的自相关特性,提出相关检测地电传输特性的编码电磁勘探法.分析了编码电磁测深原理及参数提取方法.同时记录发送电流信号和多收发距电磁场响应,对源信号和场信号进行相关运算,解卷积分离接收系统响应后可得到大地的频率特性或时间特性,以此实现地电断面的精细探测.在时间域,大地冲激响应和阶跃响应含有丰富的地电信息.通过冲激响应的峰值时间或阶跃响应的晚期渐近值可估计地电阻率分布.基于层状模型的大地冲激响应和阶跃响应正演计算结果表明,编码电磁测深法对大埋深薄层目标体有精细的分辨能力.可应用于能源、矿产资源、水资源、环境地质及工程地质勘察,有良好的应用前景.     

编码源地电阻率观测试验

地电阻率的高精度观测是实现地震预测预报的前提之一.面对日益严重的电磁干扰,提出了基于编码源循环互相关辨识技术的地电阻率观测方法,其实质是将待测地质体视为待辨识系统,利用编码源信号激励供电电极A和B产生的电流信号作为系统输入信号,测量电极M与N之间的电压信号作为系统响应输出,将输入和输出信号严格同步采样为时间序列,分别与参考信号进行循环互相关运算并转换至频率域计算获得待探测地电阻率谱(幅度和相位).由于系统环境的干扰和随机噪声与编码源信号不相关,通过循环互相关运算可以达到抑制环境随机噪声和干扰的目的.这种地电阻率观测体系在环境干扰较大的甘肃省兰州观象台和陇南汉王地震台站利用现有的观测场地和线路进行了观测试验,测量结果显示,数据的一致性好、均方差小,说明该方法在强干扰环境下具有较好的抗电磁干扰能力,观测频带较以往直流电法测量有较大的扩展.该方法为现有地面地电阻率台站持续发展提供了技术保障,可为地震预报与科学研究提供高质量的地电阻率观测数据.      

陆地可控源电磁法探测效果的频率响应

陆地可控源电磁法的观测资料可依据频段范围近似地划分为近区场、中间区场及远区场,但采用测量相互正交电、磁分量,并计算视电阻率的资料处理方式只适用于远区场数据.为更有效地利用陆地可控源电磁法不同区间场的观测资料,本文结合三维数值模拟技术并采用电场分量直接进行反演的策略,对不同区间电场的响应特征与探测效果进行了分析.数值模拟结果表明:近区电场的异常响应最明显,异常响应不随频率发生显著变化,但纵向分辨能力差;远区电场异常响应随频率发生显著变化,其探测深度取决于频率的高低;中间区场较为复杂,地表电场异常响应的等值线中心并不是位于异常体中心正上方,而是在沿场传播方向上向异常体与围岩的分界面处偏移,并且发现中间区场资料的加入会影响反演结果的准确性.综合合成数据和野外实测资料的反演结果,发现结合近区场和远区场资料而舍弃中间区场资料的反演效果更佳,这为陆地可控源电磁法资料的反演解释提供了一种有效途径.     

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

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

电性源瞬变电磁地空逆合成孔径成像

电性源地空瞬变电磁法具有工作效率高、勘探深度大、采集信号信噪比高、适用于地形地质条件复杂地区等优点.但是,到目前为止,由于尚未建立起该方法的解释系统,大大制约了该方法的发展.本文旨在建立起完整的地空电磁探测系统,丰富整个探测系统的理论.本文围绕地空瞬变电磁法全域视电阻率定义、瞬变电磁虚拟波场的克希霍夫偏移成像、逆合成孔径成像方法三个科学问题进行了系统研究.提出了用磁场强度定义全域视电阻率的迭代算法,理论模型试验结果表明计算出的视电阻率曲线首支趋于第一层电阻率,尾支趋于最后一层电阻率,实现了全空域、全时域视电阻率的计算;在先前研究的基础上,实现了适合电性源地空装置的瞬变电磁虚拟波场的克希霍夫偏移成像;采用相关迭加技术,实现孔径内多测点数据合成,将传统的单点处理方式发展成为逐点推移多次覆盖的逆合成孔径处理方法.层状模型试验表明:(1)全域视电阻率能够光滑、完整、渐变地反映出模型的电性信息变化;(2)当改变三层模型中间层电阻率时,全域视电阻率曲线随着参数的改变分异明显,对电性层的识别容易且直观;(3)由于在电阻率计算中同时考虑了接收机高度、偏移距、时间等各参数的影响,全域视电阻率可实现全空域、全时域的视电阻率计算.含水采空区的复杂模型算例表明:(1)根据不同测线的全域视电阻率结果可以看出,在靠近采空区的位置,全域视电阻率断面可以清晰地反映出采空区的空间位置,随着测线离采空区越来越远,采空区异常越来越弱直至消失;(2)波场变换和偏移成像的结果显示存在两个电性差异较大的界面,上界面指示地表,由于空气和大地之间的电性差异较大,故该界面波场信号反映强烈,遍布整个区域,下界面异常信号则主要集中在中部,向外逐渐减弱,指示采空区;(3)逆合成孔径成像结果表明地表界面在合成前后没有变化,而采空区异常合成后范围明显变小,且异常边界清晰,指示的采空区位置与模型吻合很好.本文借助于逆合成孔径雷达成像的基本思想,建立了一套电性源瞬变电磁地空逆合成孔径成像方法.基于反函数思想结合迭代算法提出的电性源地空瞬变电磁法的全域视电阻率定义方法,实现了全空域、全时域的视电阻率计算;借鉴瞬变电磁拟地震偏移成像算法,实现了瞬变场的三维成像;借鉴逆合成孔径雷达的思想,提出电性源瞬变电磁地空逆合成孔径算法,进一步提高了成像的分辨率.采空区模型算例表明相关叠加合成确实具有增强有用信号、提高信噪比、提高分辨率的诸多优点,证实了瞬变电磁地空逆合成孔径成像方法的有效性.     

多辐射场源地空瞬变电磁法多分量全域视电阻率定义

本文针对多辐射场源地空瞬变电磁法理论,建立了相应的多分量全域视电阻率定义方法.分析了利用磁场强度进行多辐射场源地空系统全域视电阻率定义的优点,针对磁场强度的各个分量,提出了各自的全域视电阻率算法,实现了多分量、全时域、全空域视电阻率计算,并分析了偏移距对全域视电阻率的影响.通过调整源的相对位置及电流方向等参数,多辐射源瞬变电磁地空系统不仅可以加强不同分量信号强度,削弱随机干扰,还可更好地分辨地下异常体的位置.通过对所设计模型的处理,证实了多辐射场源地空系统多分量全域视电阻率算法的有效性,也验证了多辐射场源地空系统的优势.     

多通道瞬变电磁法2D有限元模拟

多通道瞬变电磁法的反演工作是在大地脉冲响应的基础上进行的,因此本文用2D有限元正演来求取大地脉冲响应,实现多道瞬变电磁法的正演模拟研究.研究工作首先在有源大地电磁法2D有限元模拟的基础上进行,原有2D有限元法的频带相对于多通道瞬变电磁法的频带是窄频的,为此采用校正法将其扩展到所需的宽度,实现了频率域宽频带有源电磁勘探方法的正演模拟;然后通过频时变换变换到时间域,得到瞬变电磁法的阶跃响应;最后通过求取阶跃响应的时间导数,得到大地脉冲响应.针对多通道瞬变电磁方法的装置特点,用有一定埋深的顺层成矿模型分析了大地脉冲响应特征,发现对有一定埋深的顺层产出的矿脉模型来说,频率域电磁场结果以及时间域大地脉冲响应结果均需移除背景场才能突出异常体的存在,然而,在实际工作中,未知的背景场限制了其实用性.为此,本文利用多通道瞬变电磁方法数据量大的特点,提出用不同偏移距的等时曲线与等时断面来展示脉冲响应模拟结果,事实表明这两种展示方式可以清晰地分辨出矿体引起的异常,且能准确定位矿体的中心在地面的投影位置,从而说明多通道瞬变电磁方法相对于频率域有源电磁方法,对有一定埋深的顺层矿床有更高的分辨能力.     

海底油气藏及天然气水合物的时频电磁辨识

提出了伪随机多频海洋电磁法观测方案.采用伪随机多频信号作为激励场源,多偏移距同线偶极-偶极同时观测,相关辨识海底地电系统的频率特性和冲激响应,可以在时间域和频率域同时辨识海底高阻薄层.在时间域,瞬变冲激时刻可以直接指示海底地层电导率的变化;在频率域,利用多个频率的电场响应计算的频散率及其道闻变化量,相对相位道间变化量对高阻薄层有很好的反映.从而实现对海底油气及天然气水合物的多参数辨识.     

水平多层大地上垂直磁偶极频率测深的全波视电阻率

本文根据视电阻率定义的原则,以及用不同的场量定义的视电阻率效果不同这一事实,提出一种新的全波视电阻率定义.在全区同时用均匀大地上电磁场的三个分量来分区定义祝电阻率.在远区视电阻率由磁场的水平分量求出,在近区由磁场的垂直分量或其实分量定义,而在过渡区则由电场的水平分量确定.用这种方法定义的视电阻率为电磁响应的单值函数,它随频率变化的曲线显著改善,能直观地反映地层电阻率随深度的变化,数值比较接近地层的真电阻率值,假极值效应明显压低.在计算中用切比雪夫多项式分段拟合均匀大地电磁响应的反函数,并给出一套系数,由此算出的视电阻率误差小于1％.     

基于频谱分析识别电磁波信号干扰与地震异常

处理北京地区2009-2012年5个高采样率数字CNEM08-Ⅰ型电扰动台站数据,通过频谱分析,研究数字电磁波正常本底信号,地铁、地电供电、电磁感应等于扰信号,以及地震异常信号的频谱特征.结果表明：在时间域,数字电磁波的持续型干扰信号与震前电磁波信号表现相似,为连续高值脉冲信号;而在频域,数字电磁波干扰信号往往影响频带较窄,优势频率能量突出,优势频率频带窄到近乎以谱线的形式出现,谱线具有离散性、谐波性、衰减性特征;而震前信号具有连续的较宽频带,衰减较慢.