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

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

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

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

可控源电磁三维频率域有限元模拟

本文采用电磁场的磁矢量位和电标量势,将Maxwell方程组化为位势的类似于Helmholtz型方程,并引入罚项及稳定化方法克服了电磁三维有限元算法中的伪解及数值不稳定性;采用人工边界把计算区域局域化,将均匀半空间中水平电偶极子源产生的位势值作为人工边界上的第一类边界条件以表示源的作用,减少了实际的计算区域.理论模型和复杂模型的计算结果均表明,可控源电磁三维有限元数值模拟给出了稳定、可靠的电磁场分布.     

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

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

频率域海洋可控源电磁垂直各向异性三维反演

地层宏观电性各向异性会对可控源电磁响应产生重要影响.由于海底地层电性结构常表现为电导率各向异性,若仅对海洋可控源电磁(MCSEM)数据进行常规各向同性反演,有可能无法获得准确的反演解释结果,从而削弱MCSEM技术的可靠性.本文实现了电导率垂直各向异性(VTI)条件下频率域海洋可控源电磁数据三维反演算法.其中,三维正演采用基于二次场控制方程的交错网格有限体积法,并利用直接矩阵分解技术来求解离散所得的大型线性方程组,有利于快速计算多场源的响应.反演采用具有近似二次收敛性的高斯牛顿算法对目标函数进行最优化.最后,对具有VTI电性各向异性特征的盐丘构造模型的MCSEM合成数据分别进行了电导率各向同性和垂直各向异性三维反演,结果表明:各向同性三维反演算法无法对受VTI介质影响的MCSEM数据进行正确的反演解释,而垂直各向异性三维反演能够获得更为可靠的地下电阻率结构和异常体分布,展现出对海底电性各向异性结构更为优良的反演解释能力.     

三维频率域可控源电磁反演研究

对于三维可控源电磁,反演计算效率、张量测量、旁侧效应以及阴影效应是目前研究的热点.本文正演采用基于库仑规范条件的耦合势有限体积算法,反演采用有限内存BFGS(L-BFGS)算法.合成数据反演结果表明:(1)有限内存BFGS法比非线性共轭梯度法,在反演计算效率上具有一定的优势,更适合求解大规模三维可控源电磁反演问题.(2)张量可控源电磁法相对于标量可控源电磁法,前者在模型分辨率上优于后者.(3)在某个区域无法布置测网的情况下,我们可利用旁侧效应在异常体周围布置测网进行三维反演,从而获得真实异常体的信息.同时,为避免阴影效应,我们应在测网外增加可控源电磁控制点,使得三维反演的数据更加完备.     

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

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

基于二次耦合势的广域电磁法有限体积三维正演计算

与可控源音频大地电磁（CSAMT）相比,广域电磁法通过采用全区视电阻率定义,突破了卡尼亚视电阻率所需的远区条件限制,极大拓展了可控源电磁观测区域和探测深度.考虑到电偶源激发场的三维特征以及地下复杂三维结构,为提高广域电磁数据解释精度,本文实现了基于二次耦合势的广域电磁法三维正演计算.该算法利用Helmholtz定理将麦克斯韦方程转化为库伦规范下的磁矢势和电标势耦合方程,有效改善了离散所得大型线性方程组的谱性质,并通过强加散度条件来消除电场伪解的影响.此外,采用散射场方法,其中一次场使用准解析法求解,二次场使用有限体积法求解,克服了局部激发场源奇异性问题.通过与一维层状模型下电偶源产生的电磁场准解析解对比,验证了本文算法的正确性.在此基础上,利用本文的正演算法对比分析了广域电磁法与CSAMT对典型三维目标体的探测能力,结果表明在相同的观测条件下,广域电磁法能够更准确地反映地下目标体信息,拥有更优的分辨能力.     

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

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

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

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

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

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

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

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

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

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

基于非正交网格的带地形三维瞬变电磁场模拟

本文采用非正交网格时间域有限差分（FDTD）算法开展带地形瞬变电磁场的正演研究,并提出了一种采用非等权投影的改进非正交网格算法,提高了带地形瞬变电磁场三维正演算法的精度.对不同地形模型的数值模拟结果表明,地形的深度、宽度以及地形与源的相对位置等都可能对瞬变电磁场产生影响.     

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

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