电导率任意各向异性海洋可控源电磁三维矢量有限元数值模拟

现有海洋可控源电磁三维数值模拟方法大多基于电导率各向同性介质理论,不能模拟海底地层电导率各向异性的实际情况.本文给出了电导率各向异性三维介质中电性源海洋可控源电磁二次电场的边值问题以及相应的变分问题,采用长方体单元对研究区域剖分,将场分量定义在剖分单元的边上,利用矢量有限单元法求解变分问题,实现了电导率任意各向异性海洋可控源电磁三维矢量有限元数值模拟.这个新的正演方法可以计算电导率任意各向异性三维地电模型的海洋可控源电磁响应,基于二次场矢量有限元法直接求解电磁场,避免了传统有限元方法可能遇到的伪解问题和难于处理电场法向分量不连续的问题,提高了数值模拟计算精度.一维电导率各向异性模型电磁场数值解与解析解吻合得相当好,无论在源附近还是远离源处相对误差均不超过1%.电导率各向异性二维模型的计算结果与已有文献采用的非结构有限元模拟结果十分吻合.三维地电模型数值模拟结果显示,电导率各向异性张量电导率主轴分量和欧拉角对不同装置海洋可控源电磁响应均有着明显的影响.     

海洋可控源电磁三维非结构矢量有限元数值模拟

本文实现了海洋可控源电磁三维矢量有限元数值模拟.由于采用非结构四面体单元进行三维网格剖分,该方法可以模拟复杂电性异常体和海底地形.一维模型的数值模拟结果表明,电场实、虚部均与解析解吻合得相当好,计算误差基本小于1%.二维模型的计算结果与已有的二维自适应非结构有限元模拟结果吻合很好.带地形模型的数值模拟结果显示,海底地形对电场影响很大,有可能掩盖海底油气藏产生的异常.     

面向目标自适应海洋可控源电磁三维矢量有限元正演

本文实现了一种面向目标自适应海洋可控源电磁三维矢量有限元方法.为满足三维复杂电性结构模拟的需求,网格剖分采用非结构化六面体.在组装刚度矩阵之后,形成的大型复数线性方程组分解为等价的实数形式,利用带预条件的广义最小残差法进行求解.在获得微分方程的解之后,为提高解的准确性,通过面向目标的自适应误差估计来指示网格细化,重点加密能使观测点数值模拟精度提高的网格.对于大规模三维数据,为了使模型空间的并行计算达到均衡负载的效果,我们使用METIS函数库来进行网格计算任务量的划分.最后,通过对比一维解析解与三维自适应矢量有限元计算结果,验证了程序的正确性;通过自适应过程中误差指示子的分布,验证了面向目标自适应的有效性;通过对三维复杂模型进行均衡负载下的并行计算,测试了程序的可扩展性.     

2.5维电导率正交各向异性海洋可控源电磁等参有限元数值模拟

本文实现了2.5维电导率正交各向异性海洋可控源电磁等参有限元数值模拟.利用傅里叶变换导出了电导率正交各向异性2.5维海洋可控源电磁法波数域电磁场耦合方程,采用伽里金加权余量法推导了相应的有限元方程;采用任意四边形单元对研究区域进行剖分,在单元中进行双二次插值,将有限元方程化为线性代数方程组;最后,求解线性方程组并进行反傅里叶变换获得空间域电磁场值.这个方法可以模拟海底起伏地形条件下地下任意形状电导率正交各向异性的复杂模型.与一维模型的数值模拟结果对比表明,电磁场数值解与解析解吻合.二维模型的计算结果与二维自适应非结构有限元模拟结果也吻合.水平海底二维地电模型考察了不同各向异性系数对海洋可控源电磁响应的影响特征.海底起伏地形地电模型的数值结果表明,电导率各向异性对海洋可控源电磁响应影响明显,有可能淹没海底地形和高阻油气藏引起的异常.     

基于二次场的海洋可控源电磁法三维有限元正演

根据库伦规范势的定义,推导出关于磁矢量势和电标量势的偏微分方程,为了克服由电流源引起的奇异性和数值模拟计算困难,将电磁总场分解为一次场和二次场,一次场由基于Schelkunoff势函数的一维正演算法得到,二次场由有限元法计算得到,实现了海洋可控源电磁法三维有限元正演算法。通过一维数值模拟实例,验证该算法的计算精度。然后,利用该算法对带海底地形的三层储层模型进行正演,分析了海底地形对海洋控源电磁场各分量产生的影响。      

基于局部加密非结构网格的海洋可控源电磁法三维有限元正演

本文基于非结构网格实现了海洋可控源电磁法三维有限元正演模拟.该算法采用完全非结构网格剖分,可以模拟任意起伏地形和复杂地电模型.为了避免场源的奇异性,采用一次场/二次场分解算法,一次场由基于Schelkunoff势函数的一维解析公式得到.为了提高算法的精度和效率,采用对测点附近单元和异常体区域进行体积约束加密的方法,实现了非结构网格的局部加密.一、二维模型计算和分析表明,本文采用的局部加密方法能够明显地改善算法的精度,最大相对误差基本在1%以内.对三维模型计算及对比分析,说明了该算法对三维可控源电磁正演的实用性.复杂海底地形模型的正演模拟表明,海底地形对电磁场的影响很大,在进行海洋可控源电磁资料解释时,地形的影响有必要考虑在内.     

基于改进的接收点插值算法的频率域海洋可控源电磁法2.5维正演

在海洋可控源电磁法勘探中,接收站常置于海底.在进行海洋电磁场模拟时,由于海水和海底介质存在显著电性差异,这给海底接收点处场值的求取带来困难.本文提出一种新的接收点插值算法,该算法考虑到海底电场法向分量不连续性问题,用法向电流分量进行插值以准确求取海底任意接收点处电磁场值.本文利用交错网格有限差分法实现了二维介质中频率域海洋可控源法(CSEM)正演.对构造走向做傅里叶变换,将三维电磁模拟问题转换为波数域2.5维问题,即三维场源激励下针对二维地电模型的电磁模拟问题.使用交错网格有限差分法,基于一次场/二次场分离方法导出波数域二次电场离散形式,并进一步求得波数域电磁场.采用本文提出的改进的插值算法可求得海底任意接收点处波数域电磁场,采用傅里叶逆变换对波数域电磁场进行积分可得到接收点处空间域电磁场.模型算例表明,与常规的线性插值和严格插值算法相比,本文提出的改进的插值算法具有更高的精度.     

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

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

时间域航空电磁各向异性大地三维自适应有限元正演研究

各向异性介质模型电性结构复杂,如何进行合理的网格剖分成为获得高精度正演结果的关键,为此本文开展时间域航空电磁各向异性大地三维自适应有限元正演算法研究.通过结合非结构时间域有限元算法和自适应网格优化技术,实现各向异性介质条件下三维时间域航空电磁自适应正演.考虑到时间域航空电磁响应随时间的衰减特性,为了综合评价不同时刻的后验误差,本文将时间作为加权因子,调整各个时刻后验误差的相对权重,进而实现对浅部和深部网格的同步优化.通过与一维解析结果进行对比验证了本文算法的可靠性.数值实验结果显示电导率各向异性对自适应网格影响严重,其最大主轴电导率的数值及其分布特征直接决定了网格加密效果.此外,各向异性对时间域航空电磁三分量响应的分布形态和异常幅值也会产生严重影响,利用全域视电阻率极性图,可以很好地识别各向异性主轴方向.     

基于拟态有限体积法的频率域可控源三维正演计算

大规模地球物理电磁数据的定量解释需要发展高效、稳定的三维正反演算法.本文通过求解离散化的三维电场矢量Helmholtz方程,实现了基于有限体积法的频率域可控源电磁(CSEM)三维正演算法.为模拟具有强电性差异的三维电性介质,该算法采用拟态有限体积法(MFV)对Maxwell方程组进行离散化;另外,为获得稳定、高精度的正演数值结果,采用直接矩阵分解技术来求解离散所得到的大型稀疏线性方程组.对于具有多个发射源的CSEM测量来说,一次矩阵分解结果能够用于同频率下所有场源的正演计算.为降低场源奇异性及边界条件对数值精度的影响,采用虚拟场源校正技术,避免了散射场公式中在构建场源项时所需的大量时间.对于具有多个频率的CSEM的模拟计算,采用分频并行策略来加快三维正演计算.最后,通过与一维层状模型及三维模型的数值结果的对比验证了本文所开发的正演算法对频率域CSEM模拟计算的准确性及有效性,表明该正演算法能够有效应用于三维介质的数值计算.另外,对于多频率CSEM的并行测试结果表明基于分频并行策略的并行计算能够显著地降低正演计算时间.     

一种激发极化法2.5维正演的自适应有限元方法（英文）

传统的基于结构化网格有限元法采用的单元比较规则如矩形等,且网格剖分和加密要靠手动实现,所以传统的基于结构化网格有限元法不能准确和灵活地模拟复杂介质。本文采用易于模拟复杂介质模型的非结构化三角形网格进行剖分,且利用对偶加权后验误差估计指导网格自动细化过程,然后在电位模拟的基础上计算雅可比偏导矩阵,并依据Seigel（1959）理论实现激发极化法2.5维自适应有限元正演模拟算法。通过对垂直接触面模型进行正演分析,接收点附近网格得到了明显加密,电位数值解平均相对误差收敛到0.4%,视极化率平均相对误差收敛到1.2%,表明经自适应网格细化后,该算法数值解最终能收敛到精确解附近。最后对两个较复杂模型进行了正演计算与分析,进一步验证了该算法的准确性和灵活性。     

基于局部加密非结构化网格的三维电阻率法有限元数值模拟

总结了目前常用的结构化网格及其局限性,分析了非结构化网格对复杂地质体边界的适应性和Delaunay三角化算法及其网格加密策略,提出了一种新的局部节点加密方法,实现了对复杂模型的完全非结构化四面体全自动剖分,给出了三维直流电阻率模拟中四面体网格的质量评价标准和最优指标.计算和分析表明,数值解在点电流源及附近的奇异区精度最低,网格加密策略可以有效地减少其影响,极大地提高数值解的精度.本文提出的局部节点加密策略计算量最小,对精度的改善也优于经典方法.在精度要求苛刻或模型十分复杂时,局部体积加密策略和二次单元是高精度模拟的可靠保证.     

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

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

面向目标自适应三维大地电磁正演模拟

本文将面向目标的自适应算法应用于三维大地电磁数值模拟.使用基于非结构网格的矢量有限单元法对起伏地表大地电磁正演模拟问题进行求解.使用利用垂向电流密度在物性界面上的连续性对后验误差进行估算的算法指导网格优化.由于全局自适应算法针对观测点优化网格的能力较差,本文通过求解正演问题的对偶问题计算后验误差的加权系数,并对相关加权系数进行改进,从而实现了面向目标的自适应算法.与传统基于结构化网格的电磁正演算法相比,采用非结构网格能够更好地拟合起伏地表和地下不规则异常体.由于使用了面向目标的自适应算法,本文能够使用更少的网格达到较高的计算精度.通过对比本文模拟结果与半空间响应和全局自适应算法计算结果,并通过对比使用改进前和改进后加权系数得到的网格剖分结果验证了本文算法的有效性.     

基于电流密度连续性条件的直流电阻率各向异性问题自适应有限元模拟

直流电阻率法被广泛应用在工程和环境及水文地球物理、野外采矿、地热探测等领域.地下岩石常具有层理面和裂缝等具有方向依赖性的结构,岩石电导率常常呈现各向异性特征,因此研究复杂直流电阻率各向异性问题的高精度正演算法具有迫切的理论和学术需求.本文利用面向目标的自适应有限元算法和非结构化网格相结合的方式,解决了带地形任意复杂直流电阻率各向异性问题的高精度正演这一难题.有别于前人的研究成果,本文提出了一种特别的二次虚拟场算法实现带源的任意起伏地形问题模拟;另外,本文第一次基于电流密度连续性条件构建适合直流电阻率各向异性问题的后验误差估计算法,有效地驱动面向目标有限元网格自适应加密过程.最后,通过三组电阻率各向异性模型验证本文提出算法的正确性和适应性,测试结果表明：对于任意复杂直流电阻率各向异性问题,本文提出的算法具有精度高、适应性强等特点;另外,我们还发现电流密度连续性条件可用于设计直流电阻率问题的有效后验误差估计算法.     

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

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