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

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

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

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

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

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

     

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

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

     

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

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

     

基于区域分解法的海洋可控源电磁三维正演模拟研究

为准确刻画复杂海底地质结构中电磁场变化特征,有效处理大规模海洋可控源电磁三维正演模拟问题,本文将基于矢量有限元的 FETI-DP(Dual-Primal Finite-Element Tearing and Interconnecting)区域分解法引入到海洋电磁三维正演模拟中 该方法首先利用非结构四面体网格对海洋地电...     

基于Lanczos-模型降阶算法的三维频率域可控源电磁快速正演

快速的电磁法正演模拟是反演的核心要素之一.本文针对陆地可控源电磁(CSEM)问题, 实现了三维多频的快速正演.首先基于电场的微分控制方程, 利用非结构化网格和伽辽金方法建立关于电场矢量的有限元线性方程组.为快速求解大规模复系数有限元线性方程组, 本文引入了基于有理Krylov子空间的单极点模型降阶算法.该方法只需求解实系数矩阵方程, 然后使用Lanczos算法快速构建有理Krylov正交基子空间, 达到系数矩阵维数降阶的目的, 从而加快正演求解速度.层状模型和低阻块体模型的计算模拟表明: 在精度满足条件的情况下, 同常规有限元方法相比, 加速比达到10~20多倍.对于百万未知数级、几十至上百个频率的三维可控源电磁正演问题, 在个人PC机上只需几分钟便可快速求解.

     

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

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

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

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

基于谱元法的频率域三维海洋可控源电磁正演模拟

高精度、快速有效的正演模拟算法是三维电磁正反演的前提.为了提高海洋电磁三维数值模拟的精度和效率,本文提出利用基于Gauss-Lobatto-Chebyshev（GLC）基函数的谱元法进行海洋可控源三维电磁正演模拟.谱元法结合有限元法和谱方法的优点.我们通过应用伽辽金加权残差法离散二次电场矢量亥姆赫兹方程,在单元内选择混合阶GLC多项式的张量积作为高阶矢量插值基函数,在求解大型稀疏线性方程组时利用直接求解器进行快速求解,从而实现了三维海洋可控源电磁快速高精度正演模拟.一维和三维模型正演结果验证了本文算法的有效性和准确性.典型模型的数值结果表明谱元法是一种有效的三维海洋可控源电磁正演数值方法,能在稀疏网格剖分情况下获得精确的海洋电磁正演模拟响应.     

三维CSAMT法非结构化网格有限元数值模拟

考虑到可控源音频大地电磁法（CSAMT）电偶极发射源与地下介质的三维结构特点,本文采用非结构化网格剖分技术,开展了三维CSAMT方法有限元数值模拟研究,将三维电磁场的背景场和异常场分别求解,避免了电偶极发射源的奇异性问题,并减小了计算区域.推导了三维异常电场遵循的有限元方程,加入散度条件进行约束以消除电场伪解;对非结构化网格单元采用高斯加权平均算法,得到了精度较高的异常磁场.针对层状介质模型,与积分方程法对比,验证了有限元算法的正确性;计算分析了典型三维地质模型的电磁响应,异常体反映明显.结果表明本文算法正确、可靠,适用于三维地质模型的CSAMT方法正反演研究.

     

基于三维模拟的海洋CSEM资料处理

海洋可控源电磁法已经成为海洋油气勘探一个重要工具,但是其资料处理和解释还处于定性和一维模拟阶段.在积分方程三维模拟的基础上对Troll油田实测数据进行了处理,采用人机交互三维模拟寻找背景模型和异常体初始模型,最后对异常体电阻率采用准线性近似快速反演,取得了定量的结果.同时,说明对于二维测线和二维模型依然可以用三维来模拟,其结果优于二维反演.在电子计算机技术快速发展的今天,可以预计三维反演将成为资料处理解释的主流.     

基于微增模型的海洋可控源电磁法三维非结构化矢量有限元数值模拟

为了分析海洋中含金属矿的储油气砂层的成分、结构对于海洋可控源电磁场的影响特征,引入多相微增等效介质模型,结合非结构化网格实现了海洋可控源电磁三维矢量有限元正演.首先对多相微增模型进行了介绍,并分析不同参数影响下的等效电导率变化特征.而后,详细推导了频率域可控源电磁法非结构化矢量有限元方程.接着,利用预处理的IDR（s）迭代算法求解线性方程组.最后设计典型的海洋地电模型,通过正演计算,验证了本文算法的正确性和有效性,同时研究了海洋储油砂层中的金属矿的含量、孔隙度、含水饱和度等参数对于可控源电磁场响应特征的影响.

     

起伏海底地形时间域海洋电磁三维自适应正演模拟

本文基于自适应非结构有限元算法实现海洋电磁起伏海底地形三维正演模拟.通过采用隐式后推欧拉时间离散技术,保证在较大的时间步长条件下获得正确结果.为获得多时间道海洋电磁正演模拟的有效网格,我们采用基于法向电流连续的后验误差估计的自适应方法和网格融合技术;同时为了控制网格数量和保证正演模拟稳定性,我们还在网格融合过程中应用了随机网格挑选技术.对于方程组求解我们使用MUMPUS直接求解器.当时间步长不变时,只需对系数矩阵进行一次分解,大大提高计算速度.将本文计算结果与半空间模型解析解进行对比,验证了本文算法精度.针对海洋电磁半拖曳式和双船拖曳式工作方式,我们通过典型模型的模拟计算,研究海底地形影响及海底高阻层识别特征.

     

电导率各向异性的海洋电磁三维有限单元法正演

本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性, 文中选用二次场算法, 将背景电阻率设置为水平层状且各向异性, 场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组, 通过Galerkin加权余量法形成大型稀疏有限元方程, 采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位; 进而, 利用滑动平均方法得到二次矢量位和标量位在空间的导数, 由此得到二次电磁场; 通过一维模型对算法的可靠性进行验证, 与此同时, 针对实际复杂海洋电磁模型, 比较有限元模拟结果与积分方程模拟结果, 进一步验证算法精度.若干计算结果均表明, 文中算法具有良好的通用性, 适用于井中电磁、航空电磁, 环境地球物理等非均匀且各向异性介质中的电磁感应基础研究.     

三维三分量CSAMT法有限元正演模拟研究初探

首先从麦克斯韦方程出发,用伽里金方法推导了三维三分量CSAMT法的有限元方程.在研究过程中,认识到加入散度条件的必要性,在公式中强加了散度条件,提高了解的完备性.其次将成功应用于二维线源频率域电磁法有限元模拟中的两种技术推广到三维中,一是边界条件统一采用一阶吸收边界,使线源产生的电磁波在边界上按波的传播规律被吸收,以降低平面波假设造成的影响;二是总体系数矩阵的存储,用两个二维数组分别记录总体系数矩阵的非零元素及其在总体结点编号中所处的位置,使总体系数矩阵的存储量达到最小的同时,物理意义明确,迭代求解时迅速简便.最后用均匀半空间模型进行了验证.     

基于Coulomb规范势的电导率呈任意各向异性海洋可控源电磁三维非结构化有限元数值模拟

本文以非结构化网格剖分为基础,推导了电导率呈任意各向异性条件下的海洋可控源电磁法二次场磁矢量位、标量位所满足的有限元方程.将不完全LU分解预条件因子（ILU）与Induced dimension reduction（IDR（s））迭代算法相结合对有限元方程进行求解.之后,采用加权移动最小二乘法对二次场矢量位、标量位进行求导得到相应的电磁场各分量.最后,开展了两个地电模型的电磁场计算与分析.结果表明：本文所提出的电导率呈任意各向异性条件下的海洋可控源电磁法三维正演算法正确;不完全LU分解预处理的IDR（s）算法计算效率明显优于常用的ILU-QMR、ILU-BICGSTAB算法;算法具有良好的通用性,可用于陆地电磁、航空电磁、井中电磁等地球物理勘探方法在电导率任意各向异性方面的研究.

     

基于多尺度有限元法的频率域航空电磁三维正演研究

正演模拟是电磁数据反演的基础, 其计算速度与精度一直是制约电磁反演的两大核心问题.在三维电磁正反演中, 传统方法通过加密网格或增加插值基函数阶数提高计算精度, 但由此也降低了计算效率, 制约了三维电磁反演的实用化.因此, 如何实现大尺度模型高精度快速正演是目前电磁三维正反演中亟需解决的问题.本文将多尺度有限元法应用到麦克斯韦方程求解中.我们首先在粗网格尺度上构建满足局部特性微分算子的多尺度基函数, 进而在粗网格尺度上对原问题进行求解, 通过建立粗细两套网格间场的映射关系, 在未知数较少的粗网格上实现电磁问题求解之后, 利用粗细两套网格间场的映射关系获取细网格上电磁场响应, 由此可以在保证计算精度前提下快速获取不同尺度电磁场正演响应, 计算速度得到很大提高.此外, 本文还基于八叉树思想进行网格优化, 进一步改善三维正演效率.我们通过对典型地电结构进行多尺度有限元正演模拟并与传统有限元结果对比验证算法的有效性.最后, 我们通过模拟加拿大Voisey's Bay卵形体镍铜硫化矿区航空电磁响应以检验本文算法模拟地下复杂异常体的能力.