海洋可控源电磁法数值滤波解及算例

讨论了海洋可控源电磁法水平电偶极子源发射时电磁场的正演计算问题,采用Guptasarma线性数值滤波公式求解。对比验证了两组Guptasarma滤波系数的计算精度,分析了电场表达式核函数的变化规律。模型算例表明,高阻薄层的存在能使电场响应增强,若要提高海洋可控源电磁法对浅部高阻薄层的分辨能力,应增大发射电流的频率并采用短偏移距测量。振幅随偏移距变化曲线同样表明该方法对高阻水平薄层有较高的分辨率,适用于高阻油气藏的识别,这些都说明开展一维正演问题的研究,可以了解地质体的电磁响应特征及变化规律,为野外资料的反演解释奠定基础。     

关于频率电磁测深几个问题的探讨(五)——频率电磁测深中电磁场波型及意义

引出水平分层大地面上水平电偶极源及垂直磁偶极源的电磁场分量表达式,对式中核函数进行了变换,分析所含的波型及其与场区的关系,指出了TE型场、TM型场和TEM型场穿透地层的能力。根据波型与地层作用特点,优化电磁场分量,以提高频率电磁测深的探测效果。      

用混合算法计算电偶源瞬变电磁响应

从水平层状介质表面上垂直磁偶源瞬变电磁响应的计算公式出发,根据互换原理,得出目前常用的以电偶极源发射、在中垂线两侧张角为θ范围内用线圈测量感应电动势装置的正演计算公式,用M. M. Goldman提出的混合算法计算水平层状介质电偶极源的瞬变电磁响应,计算结果证明了混合算法的优越性。     

各向异性介质中的浅海海洋可控源电磁响应特征

由于受到空气波的影响,浅海海洋可控源电磁数据对海底储油层的反映较弱,如何对浅海数据进行处理和解释一直是海洋电磁理论研究的热点。随着近年来海洋电磁理论的不断完善,浅海数据已经可以被较好地处理与反演,但是其解释水平仍然受基本理论研究不足的制约。针对这一现状,本文开展了海底电各向异性对浅海数据影响规律的研究。具体方法为:利用欧拉旋转建立不同的海底电性各向异性模型,然后采用交错网格有限差分法计算浅海海洋可控源电磁响应,最后通过分析同线情况下电场E_x分量的振幅和相位曲线特征以及海底电场及电流密度分布规律,分析各向异性对浅海海洋可控源电磁响应影响的物理机制,并讨论浅海各向异性情况下海洋电磁对高阻储油层的识别能力。得出的结论为各向异性介质中的浅海海洋电磁响应特征与深海有较大区别,在进行数据的处理、反演和解释时应区别于深海情况。     

海洋可控源电磁法多参数正演响应特征分析

鉴于介质的色散现象,充分考虑介质电磁参数(ε、μ、σ),获得各层介质中的电磁场和势函数。为了考察岩、矿石电磁参数、海水深度和海底地层厚度等变化时海洋电磁响应特征规律,给出典型地电模型和观测系统参数,借助高精度快速汉克尔滤波系数,采用30点高斯勒让德数值求积方法,计算水平电性源频率域可控源电磁法(CSEM)在海底各观测点电场Ez和磁场Bx的正演响应。结果表明,随海水深度变浅,电场Ez和磁场Bx振幅曲线变化幅度增大,空气波逐渐占据主导地位,当水深超过3km时,观测区内可以忽略空气波的影响;电场Ez分量对海底高阻层引起的异常大,而中低阻层引起的异常小,磁场Bx分量幅值对目标层电导率变化与Ez相同,因此,水平电性源频率域CSEM法不适宜探测海底低阻目标。并且对厚层目标层的探测能力高于薄层。针对岩层中磁导率变化的情况,海底沉积物电性均匀或者呈层状分布时,电场Ez和磁场Bx幅值均受到影响。此外,电场Ez和磁场Bx的幅值基本不受介电常数变化而影响。海洋可控源的探测效果与偏移距的选择有密切关系,有利偏移距范围为3000m^12000m。     

基于场延拓的海洋可控源电磁正演模拟及各向异性特征识别

海洋可控源电磁法近年来被广泛应用于海洋油气勘探之中,凭借对高阻含油储层的有效识别大大降低干井率,节约钻探成本。受沉积环境的影响,海底地层常常表现出宏观各向异性,而对海洋各向异性模拟算法进行研究可以更有效地进行海洋电磁数据解释。本文从频率域麦克斯韦方程出发,基于标量势分解理论,通过标量势函数与场的连续性条件,在波数域中将电磁场分别向海底深部和海水中延拓,并在海底耦合到发射源上完成电磁场的求解,进而利用汉克尔变换得到了空间域电磁场。通过对海底介质一维各向异性模型电磁响应特征的分析研究发现,覆盖层各向异性对电磁响应有明显的影响;然而,海底高阻层仅垂向电阻率和厚度对海洋电磁响应产生明显影响。通过极性图的方式可以对海底介质各向异性特征进行有效识别。     

时间域电磁场偏移成像方法研究

采用 Gaver-Stehfest( G-S)方法模拟水平电偶极源在不同介质模型中地表面的负阶跃波电磁场响应 ,然后引入逆时序偏移的观点 ,把这些电磁场响应作为实际观测值 ,利用 Stratton-Chu积分公式将其向下偏移 ,对偏移后的电磁场值进行归一化处理并对场值求偏导数 ,重建地下拟地震深度剖面图 ,实现时间域电磁场的偏移成像。选用了均匀介质、含良导和高阻薄层水平均匀层状介质模型以及高阻柱体、断层等模型分别进行计算 ,成像的结果表明能够得到较高质量的地下地电构造图 ,从重建的深度剖面图中可以识别出地下反射界面的位置 ,该图比时间剖面直观。     

电偶源频率电磁测深中的EX分量

在电偶源频率电磁测深的各场量中,电场ＥＸ分量具有分辨率同,观测信噪比高,不需作近场改正等优点。在有地表电性不均匀引起的静态偏移情况下,可通过转换相位消除静态偏移。因此在资料处理与解释中可不局限于比值电阻率的模式,使频率电磁测深具有更大的灵活性与适应性。     

海洋可控源电磁场视电阻率计算方法

这里提出了一种计算海洋可控源电磁场视电阻率的方法,分析了水平电场分量振幅与视电阻率之间的关系,并介绍了海洋可控源电磁场视电阻率解释方法的原理。利用1D模型和2D模型海洋可控源电磁场响应计算了视电阻率,验证了视电阻率解释方法的有效性,为海洋可控源电磁资料定性解释提供了一种新的方法。     

MVO曲线方法在海洋可控源电磁数值模拟中的应用研究

采用有限差分方法,对水平电偶极子源在海洋二维介质情况下电磁响应进行数值模拟。利用油气高阻异常体模型的数值模拟结果,研究MVO(Magnitude Versus Offset)曲线解释方法的有效性,进而分析高阻油气异常体的水平位置的变化,对电场和磁场幅值的MVO曲线产生的影响,得出了通过MVO曲线定性分析高阻油气异常体相对水平位置的方法。本研究可为海洋可控源电磁法数据的定性分析提供部份理论依据。     

关于频率电磁测深几个问题的探讨(六)——频率电磁测深的电磁场分布与观测参量

分析了水平电偶源与垂直磁偶源的电磁场在半空间大地中的分布特点,指出观测各场分量的合适地面部位。为了定性解释,一般将场分量计算为单分量视电阻率或比值视电阻率;做定量解释时,可直接采用场强进行反演,以减小换算误差。基于当前电偶源与磁偶源的发射功率和电磁场衰减规律,在实际应用中,水平电偶源可探测更大深度(<3 000 m),垂直磁偶源探测浅部(<500 m)。      

测量磁场水平分量H y的电性源广域电磁测深法

基于层状大地表面水平电偶极源电磁场公式系统,提出了比值法获取H y分量全区视电阻率计算公式。计算了层状地电模型多个收发距的电磁场各场分量E x、H y、H z,并用这三个场分量计算了全区视电阻率,以及E x/H y卡尼亚电阻率,并与模型的MT卡尼亚电阻率频谱曲线进行对比。发现H y全区视电阻率在赤道装置和轴向装置的广大区域对地电模型均有较好的响应,可进行单分量广域电磁测深。轴向装置H y全区视电阻率与赤道装置垂直磁场H z全区视电阻率相似,在低频段与层状模型大地电磁测深卡尼亚电阻率有相似的频率响应特征,在较小收发距条件下对大埋深基底就能有较好的响应。赤道装置和轴向装置测量的水平电场E x分量全区视电阻率在低频区均会进入“饱和”区。对于各电阻率定义方法,合适的收发距是较好地反映出地电特征深度变化的前提,多收发距测量有利于揭示深部电性特征。     

对瞬变电磁测深几个问题的思考(一)——瞬变电磁测深中偶极子源及其转换

偶极子源是电磁法研究的基本对象,是构成更复杂探测装置与资料解释的前提,也是电磁法勘探技术开发与创新的理论基础。从偶极子源的概念出发,导出了均匀大地表面上的主要偶极子源瞬变电磁场公式,并指明正垂直阶跃脉冲与负垂直阶跃脉冲激发的瞬变电磁场公式的异同,以加深认识。举例说明了相应于电性源与磁性源的场分量可根据电磁场的互易定理相互转换,便捷推导公式。重点阐述了由偶极子源的场强公式运用叠加原理可衍生更复杂的电磁场公式,诸如大回线源及其任意一点的瞬变电磁场公式,盘活了瞬变场的应用。      

TEM框内回线装置发射框边界影响及消除方法

探讨了大定源回线框内观测装置的视电阻率计算问题。中心回线晚期视电阻率计算由于受边界影响及不满足晚期条件,计算的视电阻率值有畸变,影响了煤田采空区及富水性的解释。利用电偶源时域垂直磁场求取全区、全期视电阻率,解决了发射线框的边界影响问题及早期道视电阻率计算问题。理论模型计算和实测资料的处理结果表明了该方法的有效性。      

瞬变测深视电阻率曲线的正演计算及其初步解释

本文以垂直磁偶极子发射,接地导线测量这样一种发一收装置为例,讨论了求解时间减电磁场问题的计算方法,具体思路是：用线性滤波法,先求解频率域问题,然后且Ｇａｖｅｒ－Ｓｔｅｈｆｅｓｔ逆拉氏变换法将频率测深正演问题转化为相对应瞬变测深（时间域）正演问题,得到了垂直磁极子场中,水平层状介质表面瞬变电磁场分量的近似表达式,由此计算了视电阻率曲线,并从理论上探讨了瞬变测深视电阻率曲线的垂向分辨能力,本文试图对中     

对瞬变电磁测深几个问题的思考(三)——瞬变电磁场关断效应及全期视电阻率的普适算法

首先给出瞬变电磁（TEM）场卷积计算式,并将关断电流变化率函数看成冲击函数,垂直阶跃脉冲激发的场分量看成输入信号,卷积结果就是任意形状关断电流的TEM响应。把电流激发脉冲看成是许多垂直阶跃函数的叠加,通过叠加这些阶跃函数的响应求得任意电流脉冲激发的瞬变电磁响应。当已知TEM场的实测值和关断电流形状,就可按编程的方法计算全期视电阻率。文中给出了2个斜坡关断电流激发的电偶极源TEM垂直磁场算例,一个为均匀半空间的响应及其全期视电阻率,另一个为三层K型地电断面的响应及其全期视电阻率。      

A circular electric dipole: a transmitter for TEM surveys

Much experience has been gained lately in the use of a radial current source, the so-called circular electric dipole (CED), as a transmitter in transient electromagnetic (TEM) surveys. CED is a source of alternating transverse magnetic (TM) polarized field, a surface analog of a vertical electrical line in VES. In the course of two recent decades, the method has been developed theoretically and tested through the field practice. The respective published results are expected to provide an idea of TEM soundings with the optional use of either TE or TM mode. In this paper we report some new theoretical aspects and share our field experience of surveys with an CED system.© 2014, V.S. Sobolev IGM, Siberian Branch of the RAS. Published by Elsevier B.V. All rights reserved.     

CSAMT一维全区反演

由汉克尔变换和连分式求解水平层状大地电偶源的水平电场和垂直磁场,然后由电场和磁场求得卡尼亚视电阻率,不进行近区和中间区校正,直接用阻尼最小法对近区和中间区发生畸变的卡尼亚视电阻率反演,得到地下电性参数。模型试算结果表明:虽然近区和中间区的视电阻率不能正确反映地下电性垂向变化,但由其反演仍能得到正确的结果。