直流电阻率三维正演的代数多重网格方法

多重网格方法在求解由偏微分方程的边值问题离散所得线性系统时,具有非常高的计算效率.但常用的几何多重网格法在处理带跃变系数的偏微分方程时存在一定缺陷,限制了其应用.本文应用代数多重网格（AMG）方法求解三维直流电阻率法正演模拟形成的有限差分线性方程组,通过求解二次场的方法消除了总场中由点电源导致的奇异性,从而获得快速、精确的三维电阻率数值模拟.对两个存在大的电性差异的模型进行了模拟计算,以验证代数多重网格法的收敛效率.计算结果表明,与不完全Cholesky共轭梯度（ICCG）方法相比,代数多重网格方法具有更高的计算效率及稳定性.而且,随着三维网格节点数的增加,代数多重网格方法计算的高效性更加明显.     

频率域电磁法三维有限元正演线性方程组迭代算法

在三维频率域电磁法的正演模拟方法中,有限元方法具有计算精度高、适应性强的优点,近年来来得到了越来越多的关注.在正演过程中,主要的计算量集中在求解由偏微分方程组离散得到的线性方程组上,因此求解线性方程组关系着正演计算速度以及模拟精度.由于由有限元方法离散得到的复系数线性方程组条件数非常大,使用常规的迭代法和预条件很难收敛.目前大多数的研究工作采用直接解法,需要大量的计算机内存,限制了可求解问题的规模.本文研究了线性方程组的迭代解法,通过将复系数线性方程组转化为其实对称形式,构造分块对角预条件.在应用预条件的过程中,需要求解两个较小的实数方程,通过辅助空间解法求解.本文的算法适用于可控源电磁法和大地电磁法,对一系列的数值算例的模拟结果证明了迭代算法的效率,结果表明迭代算法可以在小于20次迭代内收敛,同时迭代次数与模型电阻率、问题规模和频率无关.     

基于几何多重网格算法的海洋可控源电磁法三轴各向异性介质的三维正演研究

海洋可控源电磁法(MCSEM)三维正反演理论现如今已经成为地球物理学研究的热点和难点之一,准确、高效、稳定的正演计算是实现快速反演计算的基础.三维正演数值模拟技术的发展已相对成熟,一些学者已将研究如何提高正演计算效率的目光转移到研究如何提高线性方场组的计算速度.为了提高MCSEM的三维正演问题的计算效率,本文首先从频域三维海洋电磁控制方程出发,然后利用Yee氏交错网格有限体积法在三维空间离散方程组,并施以第一类Dirichlet边界条件获得大型稀疏复系数线性方程组,最后引入3种不同几何多重网格迭代算法求解该线性方程组.为了检验GMG算法的正确性,通过建立一维层状油气模型,将3种GMG算法计算结果与Kerry Key等开发的二维开源程序MARE2DEM计算结果进行对比,两种程序求解电场分布的曲线能够很好的吻合,表明GMG算法能正确求解海洋电磁正演问题,且两种程序求解的相对误差数量级在1以下,表明GMG算法具有较高的求解精度.为分析GMG算法的计算效率,我们首先想要模拟出一个更加真实的海洋地下环境,将沉积(背景)层电阻率设计为三轴各向异性,然后在此环境中建立三维海洋油气油气模型,实现MCSEM三维正演计算.通过改变网格数,实施3种GMG迭代算法与GCROT迭代算法求解,结果表明:GMG算法求解三维海洋可控源电磁正演问题算法稳定,计算效率高.GMG算法作为Krylov子空间迭代算法的预条件器求解三维海洋可控源电磁正演问题,不仅能加快求解速度,而且能提高算法的稳定性.     

复杂地形三维直流电阻率有限元数值模拟

系统地论述了用有限单元法研究复杂地形条件下三维直流电阻率的正演计算技术.首先给出了三维构造中点源电场的边值问题以及相应的变分问题;然后利用有限单元法求解变分问题,采用四面体单元对研究区域进行剖分,在单元中进行三线性函数插值,将变分方程化为线性代数方程组;最后,考虑到节约计算时间,利用对称超松弛顸条件共轭梯度迭代算法求解大型线性方程组,得到了各节点的电位值,进而计算出地表的视电阻率.通过理论模型的计算检验了算法的可行性之后,给出了几种常见纯地形异常的数值模拟结果和一个组合模型的计算结果,其研究工作为研究三维直流电阻率反演奠定了基础.     

直流电阻率法2.5维正演的外推瀑布式多重网格法

引入外推瀑布式多重网格法(EXCMG)求解2.5维直流电阻率有限元计算形成的大型稀疏线性方程组,结合基于地址矩阵的压缩存贮方式以及最优化离散波数,使得2.5维电阻率正演程序的计算速度大大提高而内存需求大大减小.研究结果表明:EXCMG法的收敛速度与网格尺寸无关,计算速度明显优于不完全Cholesky共轭梯度(ICCG)...     

基于多重网格有限元法的大地电磁二维正演计算研究

为提高大地电磁正演计算速度,开展了基于多重网格有限元法的大地电磁二维正演模拟计算研究.将稳定双共轭梯度算法作为多重网格法的细网格松弛迭代算法,插值算子采用完全加权算子,限制算子设计基于网格单元面积率,使多重网格法更适于求解大型复系数方程组.二维均匀半空间模型、低阻体模型和高阻体模型的大地电磁正演模拟结果表明:当计算量较小时(网格剖分数量少),多重网格法在计算效率方面并未有优势,网格剖分数量较大时,多重网格有限元算法在收敛速度方面的优势明显,多重网格有限元法的大地电磁正演精度优于一般数值算法.这为三维多重网格有限元的大地电磁正演研究奠定了基础.     

基于多重网格的Helmholtz方程解算及简单地质模型的大地电磁正演模拟

为了提高大地电磁正演模拟的计算效率,开展了基于多重网格算法的Helmholtz方程解算研究;对比Gauss-Seidel、Richardson迭代算法的收敛性,在达到同样的误差限值下,多重网格算法仅需十余次套迭代即可满足误差精度要求,而GaussSeidel和Richardson算法分别需要近500次和5000次迭代;多重网格算法的收敛速度极快,多重网格算法的收敛速度比一般数值算法快近3个数量级;由于其是一种基于套迭代技术,该方法完成一次循环的迭代效率明显比一般数值算法低;多重网格算法极快的收敛速度为将其用于提高大地电磁正演模拟效率提供了可能.     

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

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

复杂二维/三维大地电磁的有限单元法正演模拟策略

复杂二维和三维大地电磁模型的正演数值模拟具有一定的挑战性。对于复杂的二维和三维大地电磁正演问题,我们采用有限单元法进行求解。有限单元法最后形成一个线性方程组,系数矩阵是大型稀疏的带状对称复系数矩阵,并且其条件数远大于1,为严重病态矩阵,求解其对应方程组会遇到很多困难。不完全LU分解处理的Bi-CGSTAB迭代方法可用于该线性方程组的求解,并且具有速度快、精度高和稳定性好等优点;为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分;在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间。通过对二维和三维模型电磁响应的计算,验证了算法的正确性。     

弹性波场数值模拟的隐式差分多重网格算法

为了得到稳定的弹 性波数值模拟,而不得不选择隐式差分方程;为了提高解的精度,又不得不增加节点数目, 但同时也降低了隐式迭代求解的收敛速度. 为此,本文使用隐式差分的多重网格算法进行弹 性波数值模拟,多重网格算法通过粗网格收敛较快的迭代过程求出近似解,以近似解为初值 使用细网格进行精确的迭代求解,从而加速了隐式迭代求解的过程,能够以较高计算速度、 精度、稳定性完成弹性波传播过程的数值模拟.     

基于改进Krylov子空间算法的井中激电反演

井中激电是二次找矿重要的地球物理勘探手段,快速而稳定的正反演算法有助于方法的推广和应用.本文在正演模拟中,给出了考虑井眼影响下的网格剖分方式;用右端项校正技术减小边界效应和源点奇异性引起的模拟误差;并采用循环Krylov子空间算法提高多线性方程组的求解效率.反演用Gauss-Newton法结合Jacobian-free Krylov迭代求解技术,给出了Jacobian矩阵向量积的简化计算方法;用不精确预处理共轭梯度法对模型修正量方程近似求解以减少计算量;采用不同于正演的反演网格剖分降低不适定性.数值算例验证了相关算法的有效性和可靠性.     

利用共轭梯度算法的电阻率三维有限元正演

引入对称超松弛预条件共轭梯度(SSOR PCG)迭代算法求解电阻率三维有限元计算形成的大型线性方程组，并有机结合系数矩阵的稀疏存储模式，使得三维有限元正演计算的速度大大提高而内存需求则大大减少. 该算法可方便地应用于求三维异常电位，在保持快速计算的基础上，正演计算的精度显著提高.     

基于混合阶矢量基函数的海洋可控源电磁三维谱元法数值模拟

在前人工作的基础上,本文推导了电导率任意各向异性介质的海洋可控源电磁三维谱元法正演方程.采用一次场/二次场分离算法结合混合阶矢量基函数,可以有效避免源点的奇异性的影响,从而提高数值解的精度.采用任意六面体单元离散研究区域,有利于模拟复杂地形和地电结构.利用不完全LU分解的Induced Dimension Reduction(IDR(s))迭代算法求解线性方程组,有效地提高了求解的效率.设计典型的地电模型进行正演计算,并将计算结果与有限元解进行对比,对比结果表明本文提出的基于混合阶矢量基函数的海洋可控源电磁三维谱元数值模拟算法是正确的、有效的.本文算法具有良好的通用性,可推广用于电导率呈任意各向异性的陆地电磁、井中电磁等数值模拟研究.     

电法勘探正演数值模拟的若干结果

本文前一部分叙述了点源二维电阻率法数值模拟正演计算的有限单元法,文中采用混合边界条件、用LL~T分解解线性方程组等优化措施,使二维有限单元法电阻率法正演计算的速度和精度都比目前国外流行的L.Rijo方法和程序有了进一步的提高。文中介绍了有关的方法并引入了若干计算结果。本文后一部分叙述了三维电阻率法数值模拟的积分方程近似解法,由于对K.Dieter等人提出的方法作了一些近似处理,并用迭代法求解积分方程表示式,实算结果说明在保证计算精度的条件下,提高了计算速度,从而提高了三维电场正演计算的有效性和实用性。所提出的方法容易推广到激发极化法的正演计算中。     

基于并行化直接解法的频率域可控源电磁三维正演

电磁法的三维数值模拟是一个对数值算法和计算机硬件要求都非常高的问题.对常用的微分类方法如有限单元法和有限差分法而言,求解最后所得的大型线性方程组是至关重要的一步,直接影响到正演算法的实用性.如何高效、稳定且准确地解线性方程长期以来一直是被探讨的问题.本文实现了基于线性系统直接求解技术的频率域可控源电磁(CSEM)三维正演.使用交错网格有限体积法(FV)来离散化关于二次电场的Helmholtz方程;使用直接解法取代传统的迭代解法来求解离散线性系统,即对系统矩阵进行完全LU分解,具体通过调用大规模并行矩阵直接求解器(MUMPS)来实现.基于理论模型做了一系列数值实验,首先证明了直接解法的高精度和稳定性,并考察了其内存需求、计算时间和并行可伸缩性等主要计算性能,最后检验了所开发的算法快速模拟多场源CSEM问题的能力以及对常规海洋和陆地CSEM模拟的有效性.     

三维电阻率测深数据Zohdy近似反演方法

Zohdy方法近似反演三维电阻率测深数据。正演计算采用有限单元法。反演初始模型由测量视电阻率数据给出。通过比较实测视电阻率值和预测模型计算的视电阻率值对数差来修改模型网格电阻率．为了解决任意电极距测深数据的反演,采用大、小双网格剖分。大网格反映地下电性分布情况。小网格用于实际有限元正演计算．在电阻率调整公式中加入一个迭代系数,能够加快收敛速度．并对加5％随机噪声的模型理论视电阻率测深断面数据进行反演,得到的电阻率分布与模型电阻率基本一致．     

电阻率任意各向异性三维有限元快速正演

电阻率各向异性由于其可能导致观测数据解释的很大偏差,受到越来越多的关注.由于其复杂性,三维各向异性研究国内尚未开展.本文先推导了各向异性介质中直流电法满足的控制方程并建立其变分问题,然后用有限元方法将问题转化为求解二次场电位的大型线性方程组,结合系数矩阵一维压缩存储技术和解大型线性方程组的对称超松弛预条件共轭梯度(SSOR-PCG)迭代算法,大大提高了计算速度,减少了内存需求.本文提出的算法可以任意调整介质的电阻率取向,从而实现了直流电阻率三维任意各向异性的有限元快速、精确正演计算.模型计算结果表明,各向异性介质响应相对于各向同性介质有明显差异,其介质取向对于视电阻率观测有明显影响.     

基于几何多重网格预条件技术的三维大地电磁高效正演模拟

几何多重网格法(GMG)将细网格上的大型稀疏矩阵的求解转化为较粗网格上的更容易求解的问题,从而快速求解大型稀疏方程组.但是由于大地电磁法(MT)正演模拟中涉及双旋度算子,传统GMG无法有效平滑高频误差导致其收敛慢甚至发散.为此,我们引入了四色分块高斯赛德尔法(GS)作为平滑算法,该算法局部满足电流散度为零的条件,无需额...     

基于Lorenz规范条件下磁矢势和标势耦合方程的频率域电磁法三维正演

为了克服空气层和地表耦合以及避免一次场计算,开发适合不同类型场源、不同应用范围的频率域三维正演模拟统一平台,本文从麦克斯韦基本方程出发,推导基于Lorenz规范条件的磁矢势和标势耦合方程;通过将不同类型场源分解成一系列短导线(电性)源组合,采用交错网格采样和有限体积技术对方程进行离散得到对称大型稀疏线性方程组,并采用Jacobi迭代预处理QMR(Quasi-Minimum-Residual,拟最小残差)算法进行求解,我们成功实现不同类型场源、不同应用范围的频率域电磁法三维正演模拟.通过层状模型下大地电磁法以及有限长接地导线和大回线磁性源激发下的电磁场响应模拟,并与一维解析解对比验证算法的有效性.进而,我们利用该算法平台的模拟结果对典型地电模型在不同场源激发下频率域电磁法响应特征进行对比分析.本文算法研究及实现为建立频率域电磁法三维正反演统一框架打下基础.     

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

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