复杂地下异常体的可控源电磁法积分方程正演

可控源电磁法具有分辨率高及抗干扰能力强等特点,是一种重要的地电磁勘探方法.目前,可控源电磁法的高精度正演计算一直是其核心研究问题之一.传统积分方程法一般采用近似积分公式、简单矩形网格和近似的奇异性体积分计算技术,制约了体积分方程法处理复杂地下异常体的能力,降低了计算精度.针对上述问题,本文基于完全积分公式、四面体非结构化网格和奇异体积分的精确解析解来高精度求解复杂可控源电磁模型的正演响应.首先,从电场积分公式出发,推导了可控源电磁问题满足的积分方程;其次,借助于非结构化四面体网格离散技术,实现了地下复杂异常体的有效模拟.最后,利用散度定理把强奇异值体积分转换为一系列弱奇异性的面积分公式,并通过推导获得了这些弱奇异性的面积分公式的解析解,从而最终实现三维可控源电磁问题的高精度积分求解.以块状低阻体地电模型为测试模型,采用本文提出的积分方程方法获得的数值解与其他公开数值算法解进行对比分析,其对比结果具有高度的吻合性,验证了算法的正确性;同时,设计了球状及复杂地电模型进行算法收敛性测试,进一步验证算法的正确性以及能够处理地下复杂模型的能力.

     

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

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

复杂场源形态的海洋可控源电磁三维正演

在使用电偶极发射源的可控源电磁法(CSEM)勘探中,发射源的方位、长度、形状等对观测数据有重要的影响,然而现有的大部分三维数值模拟方法没有全面地将这些因素考虑进来,很多都只能应对非常简单的场源形态,例如单一方位的点电偶极子,这有可能显著降低模拟结果的准确性.本文实现了基于交错网格有限体积(FV)离散的海洋CSEM三维正演算法,能够模拟形态相对复杂的场源,包括任意方位的有限长直导线和弯曲导线发射源.该算法使用一次场/二次场方法,只需对二次场使用FV法求解,避免了场源的奇异性问题;一次场的计算为一维正演问题,使用准解析法求解,并且只要在计算一次场时考虑复杂的场源形态便可以实现同样场源的三维正演.通过与一维理论模型的解析解对比验证了三维程序的准确性,并针对三维理论模型进行了一系列正演测试,初步考察了场源形态对三维正演结果的影响.     

三维任意各向异性介质中海洋可控源电磁法正演研究

由于海底介质受沉积环境的影响,层理发育呈现明显各向异性特征.对于海洋可控源电磁法各向异性的研究以往主要局限于一维和二维模型,为更深入了解复杂情况下海底各向异性对海洋可控源电磁响应的影响规律,本文开展三维任意各向异性介质中海洋可控源电磁法正演研究.采用交错网格有限差分技术,通过对任意各向异性介质电导率张量实行体积和空间电流密度平均,完成海洋可控源电磁二次散射电场的离散化,成功实现任意各向异性介质中海洋可控源电磁正演模拟.通过对几种典型各向异性电性模型条件下海洋电磁电场多分量响应及分布特征和各向同性情况的对比分析,总结电各向异性对海洋电磁响应的影响规律和识别方法.本文算法研究及算例可为海洋可控源电磁数据精细化处理解释提供技术支撑.     

利用积分方程法的大地电磁三维正演

利用积分方程法实现了均匀导电半空间三维大地电磁响应的数值模拟。求取张量格林函数积分时,采用二次剖分算法解决计算中奇异值问题,对于含有贝塞尔函数的积分项,利用结合连分式展开的高斯求积代替常规的快速汉克尔变换方法,确保了张量格林函数的正确计算并提高了计算精度。最后通过数值模拟结果的对比及模型试算验证了算法的正确性,所实现的三维大地电磁数值模拟算法为理论研究三维地电构造的大地电磁响应的分布规律提供了有效的工具,也为研究三维反演算法奠定了基础。     

任意各向异性介质中三维可控源音频大地电磁正演模拟

在一些地层层理发育的地区,地下介质存在显著的电各向异性,此时基于各向同性模型解释含各向异性效应的可控源音频大地电磁（CSAMT）测深观测数据会导致错误的结果.本文通过引入3×3的对称正定张量表征电导率各向异性,采用非结构四面体网格和矢量有限元方法离散电场满足的矢量Helmholtz方程,并将电磁场源等效为系列电偶极子,实现任意各向异性介质中CSAMT高效数值模拟.本文首先通过层状各向异性模型检验三维有限元算法的精度和有效性,进一步建立三维地电模型研究异常体各向异性和围岩各向异性对CSAMT响应的影响,最后使用视电阻率极性图来识别各向异性电导率主轴方向.数值模拟结果表明,各向异性电导率对CSAMT视电阻率幅值及分布规律都有很大影响,视电阻率极性图能够很好地识别各向异性主轴方向.

     

陆地可控源电磁法三维勘探的激电效应影响研究

可控源电磁法(CSEM)的正演模拟与反演解释中常忽略激电介质的影响.本文基于电磁三维正反演开源平台ModEM,分别开展了陆地CSEM三维正反演以及考虑激电介质的陆地CSEM三维正演研究,并与一维模型解析解对比验证了三维计算的准确性.数值模拟结果显示,模型含激电效应与不含激电效应的电场E_x频率响应差异明显,激电效应会导致电场E_x幅值减小、相位增加,且对E_x的影响程度与激电异常体的频谱特性及观测的收发距均相关.合成数据的陆地CSEM反演结果显示,数据中含激电效应会使得反演出的异常体在深度、规模和电阻率值上都与真实异常体存在较大偏差,忽略观测数据中的激电效应会明显影响观测的准确性.本文研究表明当陆地CSEM三维观测资料中包含激电效应时,有必要考虑其影响.     

带地形的可控源电磁法2.5维各向异性介质正演研究

目前,对于可控源电磁法各向异性介质2.5维问题,主要采用一次场、二次场分离的方法消除场源奇异性并降低截断边界对计算区域的影响.该方法数值计算精度高,但是很难适用于复杂地形条件下的数值模拟.针对复杂地形问题,基于总场的有限元方法表现出一定的优越性,然而,这种方法存在场源奇异性问题和截断边界问题.本文采用基于总场计算的方法对带地形的可控源电磁法2.5维各向异性介质进行模拟研究,推导了考虑电导率和介电常数各向异性的2.5维控制方程;引入网格加密-收缩算法降低场源奇异性的影响范围,提升数值计算效率;引入行波分解边界条件降低截断边界的影响;提出任意采样反傅里叶变换方法,快速、高精度地计算出空间域电磁场分量.理论模型数值算例中：首先,验证了本文算法的有效性;其次,对任意各向异性倾角产生的可控源电磁响应规律进行研究;最后,采用山丘模型对各向异性介质电磁场的响应规律进行了模拟和分析.

     

海洋天然气水合物可控源电磁法的响应特征:一维正演模拟

We discuss the feasibility of using controlled-source electromagnetic （CSEM） in the frequency domain for prospecting marine gas hydrates. Based on the Ocean Drilling Program （ODP） Leg 164 log data, we have established several 1-D resistivity models which have different gas hydrate concentrations. Meanwhile, we analyzed the electromagnetic response of marine gas hydrates in the frequency domain based on these models. We also studied the relationship between electrical field magnitude or phase and parameters such as receiver-transmitter distance and frequency. Our numerical modeling results provide us with a quantitative reference for exploration and resource evaluation of marine gas hydrates.     

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

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

基于边界积分方程的2维电阻率层析成像

提出了一种新的二维电阻率层析成像方案。与国内外已有成像方法相比,该方法无需计算雅可比矩阵,不必要作线性化近似,这不仅节省了计算时间,而且使成像迭代过程,稳定,可以对实际上广泛存在的大一对比度介质获得满意的成像效果,因而大大提高了电阻率层析成像技术的这用性。中给出的数值模拟结果显示了该方法的有效性。     

退化的Fourier偏移算子及其在复杂断块成像中的应用

波动方程宽角抛物逼近得到的通常是非常系数的单程波传播算子,其系数是速度横向变化的函数,因此需要利用有限差分(FD)进行数值实施. 通过对Lippmann-Schwinger单程波动积分方程的退化核逼近,本文研究了一类宽角退化算子的偏移成像. 这种退化偏移算子只用快速Fourier变换进行波场延拓,将常规的Fourier分裂步地震偏移方法(SSF)推广适应强速度横向变化介质和大角度传播波场. 退化的Fourier偏移算子通过在两个分裂步项之间作波数域线性插值来实现波场延拓,每延拓一层需要比常规的SSF地震偏移方法多一次快速Fourier变换(FFT). 通过SEG/EAGE盐丘模型和实际地震资料的应用表明,退化Fourier偏移算子能很好地对盐下的陡倾角断层和实际地震剖面上的复杂小断块和大断裂地质构造成像.     

虚拟波动域三维海洋可控源电磁场正演模拟

本文基于对应原理将似稳态条件下频率域电磁场扩散方程转换成虚拟波动域电磁场波动方程,采用高阶时域有限差分进行求解,引入复频移完全匹配层吸收边界条件,降低了内存需求,提高了计算效率,并在虚拟波动域用伪δ函数离散电偶极源,实现了虚拟波动域任意取向电偶极源三维海洋可控源电磁场高阶时域有限差分正演算法.通过与拟解析解和频率域三维可控源电磁场数值模拟结果的对比,验证了本文算法的正确性和高效性,且探讨了网格参数和边界条件对不同频率电磁场模拟结果的影响.

     

一维电阻率各向异性对海洋可控源电磁响应的影响研究

地下介质的电阻率常常表现为各向异性, 海底裂隙地层和层状沉积序列可能形成宏观电阻率各向异性.在解释海洋电磁资料时, 电阻率各向异性的影响不应该被忽略, 否则可能会得到错误的海底地电模型.作者编写了电阻率任意各向异性一维层状介质海洋可控源电磁场计算程序, 计算了电阻率各向异性层状模型的海洋可控源电磁响应, 讨论了覆盖层和高阻储层分别具有电阻率各向异性时的电磁场响应特征.     

海洋可控源电磁法发射源姿态影响研究

海洋电磁理论是建立在水平发射偶极的基础上.然而,由于海底洋流作用海洋可控源发射偶极可能产生水平摆动、倾斜和水平旋转,造成测量数据与理想状态的水平电偶极子的电磁响应有较大的偏差.为了研究发射源姿态变化对场的影响,本文提出利用欧拉旋转将发射偶极由源坐标系转换到地球坐标系,进而分析发射源姿态变化造成海洋电磁观测数据的误差分布特征.研究发现,对于发射源水平摆动和倾斜造成的海洋电磁观测数据的畸变,其分布规律取决于接收机所在区域.以发射频率1 Hz为例,对于由发射源水平摆动或倾斜导致的地球坐标系中的两个等效电偶极子分量产生的水平电场Ex,在旁线区域,二者的相位差约为360°,接收机接收到的电场信号约为两个偶极子场之和;而在同线区域,二者的相位差接近于180°,接收机接收到的信号为两者之差.在两个区域的交界处,发射源姿态对电磁信号影响最大.相比之下,发射源水平摆动只改变发射源的位置,不会改变源的性质,其影响主要发生在收发距较小的区域.本文的研究成果有助于海洋电磁系统设计时确定发射源姿态变化的参数指标,同时也为未来实际海洋电磁观测时分析和校正发射源姿态变化造成的观测数据误差提供理论基础.