各向异性介质对大地电磁观测数据的影响往往不可忽略,因此需要提高大地电磁各向异性三维反演的可靠性和有效性.为了满足大地电磁各向异性三维反演的需求,本文研究了一种基于交叉梯度结构约束的大地电磁主轴各向异性并行三维反演算法.根据大地电磁平面波理论假设,正演方程采用背景场与二次场分离的计算方式,二次场利用交错网格有限差分法求解.由于各向异性反演的多解性,本文将各向异性介质简化为主轴各向异性,并在此基础上进一步采用有限内存拟牛顿LBFGS法实现三维各向异性反演.为了提高各向异性反演的分辨率,反演目标函数中引入交叉梯度项,利用先验的结构信息,对三个方向的电阻率参数进行结构约束,最终的反演进一步利用MPI(Message Passing Interface,消息传递接口)技术实现分频并行计算,测试结果显示并行接近线性加速比.
相似文献大地电磁正演常采用电场双旋度方程,由于工作频率低、空气电导率为零,电场双旋度方程含有空解空间,导致离散化的线性方程组严重病态、具有高的条件数.为降低大地电磁场双旋度正演方程的病态性,本文提出了一种基于空气与地球分离策略的大地电磁混合正演公式.首先,在空气层中将双旋度算子转换为Laplace算子,在地下利用电导率大于零的特征保留双旋度算子,建立空气与地下分离的混合公式系统,以便实现混合公式系统对应的离散化线性方程组条件数的大幅度降低.然后,在空气与地下分界面上,施加电磁场的切向与法向连续性条件实现空气与地下区域的耦合与信息交换,以保证电磁场分布的唯一性.再则,利用非结构化四面体网格离散技术来高精度模拟任意起伏地形和复杂地下三维地电结构,并结合节点有限元与棱边有限元法,高精度求解该混合公式系统.最后,利用国际标准大地电磁测试模型来验证该空气与地球分离算法的正确性,并展示其有效降低原有电场双旋度方程的病态性的能力.本文提出的空气与地球分离思想为寻求大地电磁最优化控制方程开启了新的求解思路,预计会产生更多、性能更佳的大地电磁混合公式,以便服务于大地电磁正反演技术的快速发展.
相似文献本文将大地电磁场分解为一次场和二次场,应用交错网格有限差分法模拟计算大地电磁二次场,并引入各向异性最佳匹配层(PML)吸收边界条件作为二次场边界条件,实现了耦合PML吸收边界条件的三维大地电磁二次场有限差分正演模拟.为了确保正演的稳定性和效率,QMR求解器和磁感应矢量散度校正技术被用于PML吸收边界条件下系数矩阵的快速求解.三维模型正演响应表明,基于二次场的三维大地电磁有限差分算法具有较高的计算精度和可靠性.通过计算分析不同PML吸收因子条件的大地电磁正演结果,显示在适当的吸收因子下,PML吸收边界条件可较大幅度的减小外边界距离,从而有效的压缩模型求解空间,最终提高三维大地电磁正演模拟的效率.
相似文献几何多重网格法(GMG)将细网格上的大型稀疏矩阵的求解转化为较粗网格上的更容易求解的问题, 从而快速求解大型稀疏方程组.但是由于大地电磁法(MT)正演模拟中涉及双旋度算子, 传统GMG无法有效平滑高频误差导致其收敛慢甚至发散.为此, 我们引入了四色分块高斯赛德尔法(GS)作为平滑算法, 该算法局部满足电流散度为零的条件, 无需额外的散度校正且具有高度并行性, 可以显著提高GMG的收敛效率.但是随着系数矩阵各向异性(比如电导率的强烈变化等)增加, GMG收敛速度会变慢.Krylov子空间求解法如稳定双共轭梯度法(BiCGstab)可以改善这种收敛变慢的问题.因此, 在本文中针对交错网格有限差分法(FDM)提出了一种结合四色分块GS平滑算法GMG和BiCGstab的MT高效正演模拟方法.在该方法中, 将四色分块GS平滑算法GMG作为BiCGstab求解器的预条件技术, 从而显著提高正演效率.我们设计了一个层状电阻率模型, 通过与其解析解对比验证本文所提算法的正确性.然后设计了一个双异常体电阻率模型和一个Dublin模型1(DTM1), 基于BiCGstab, 对比了GMG预条件技术与其他传统预条件技术的数值表现, 如超松弛预条件技术(SSOR)、分块不完全LU分解预条件技术(block ILU)和高斯赛德尔预条件技术(GS).结果显示本文提出的算法在迭代次数, 计算时间和稳定性方面都远远优于传统预条件技术.对于所有例子, GMG预条件技术均能在10次以内达到收敛, 计算时间比传统预条件技术减少70%以上, 显示了本方法的稳定性和高效性.
相似文献大地电磁正演线性方程组求解的主要加速手段有并行技术和多重网格技术.传统的几何多重网格(GMG)方法依赖于均匀嵌套的正交网格,处理带跃变系数的问题时存在一定缺陷,且对各向异性问题需采用半粗化、线/面磨光等策略特殊处理,限制了GMG方法的应用.本文提出一种基于非均匀网格的外推瀑布式多重网格方法(EXCMG),快速求解三维大地电磁各向异性有限元正演形成的大型复线性方程组.首先,对库仑规范下电磁耦合势满足的向量Helmholtz方程,由Galerkin加权余量法推导有限元离散系统,得到大规模、稀疏、复线性方程组.然后,从最密的非均匀正交网格出发,逐层粗化得到一系列嵌套网格;借助Richardson外推及Lagrange二次插值技术,设计非均匀正交网格下全新的多网格延拓算子;利用前两层网格上的数值解,构造下层密网上有限元解的高精度逼近,作为多网格磨光算子——复数域不完全LU分解预处理稳定双共轭梯度算法(BiCGStab)的迭代初值,加速收敛.最后,通过典型地电模型验证算法的精度和效率.数值实验表明,本文提出的算法适用于较宽频带,加速效果明显,相较于预条件BiCGStab方法,求解效率提升数十倍.本文算法能够求解上亿自由度的超大规模、强各向异性问题,且求解问题的规模越大,算法效率优势越明显.EXCMG算法有望在其他地球物理领域得到应用.
相似文献大地电磁法可应用于岩石圈结构探测、金属固体矿产以及地热资源勘察等方面,目前三维大地电磁反演技术日趋成熟,但二维反演依然应用广泛.受台网布设及复杂地质构造影响,大地电磁数据往往结合二维与三维反演技术来减轻计算负担和改善反演效果,而评估二维与三维反演模型及其差异对准确获取地下电性结构信息至关重要.本文提出了一种基于敏感度矩阵的电性异常体边界划定以及反演模型可靠性评估方法,该方法通过求解大地电磁正演的伴随问题来构建敏感度矩阵,无需设计额外模型即可分析反演结果的可靠性.此外,本文还利用合成数据讨论了大地电磁二维与三维反演的有效探测深度.当地下结构三维性不强时,二维反演可以有效恢复深部低阻异常体的电阻率值,为深部地质解释提供定量约束.但当电性结构三维性较强时,三维反演重建的边界与真实位置更加吻合,而二维反演会产生虚假异常,不利于后续地质解释.
相似文献有别于传统基于梯度信息的反演方法在正则化约束中用总梯度逼近海塞逆矩阵的技术,本文将正则化约束问题的数据拟合项和模型光滑项分开考虑,只利用数据拟合函数的梯度信息对数据拟合项的海塞矩阵进行逼近,通过求解类高斯牛顿下降方向方程得到不依赖前几次迭代正则化因子的更精确下降方向,在求解当前迭代下降方向的过程中,通过保证右端项中两个向量的二范数在同一数量级的原则,实现了正则化因子的自动更新.对理论模型的试算表明这种自适应正则化反演方案可以在拟牛顿反演框架下基本达到OCCAM的算法稳定性,反演结果对初始模型依赖性较小,同时又无需在一次迭代中多次搜索最佳正则化因子.本文还基于此算法讨论了大地电磁各参数对于反演结果的影响,由于本文的反演结果能得到充分的正则化约束,因而在此框架下讨论阻抗和倾子在反演中的作用相对更为客观.
相似文献本文将面向目标的自适应算法应用于三维大地电磁数值模拟.使用基于非结构网格的矢量有限单元法对起伏地表大地电磁正演模拟问题进行求解.使用利用垂向电流密度在物性界面上的连续性对后验误差进行估算的算法指导网格优化.由于全局自适应算法针对观测点优化网格的能力较差,本文通过求解正演问题的对偶问题计算后验误差的加权系数,并对相关加权系数进行改进,从而实现了面向目标的自适应算法.与传统基于结构化网格的电磁正演算法相比,采用非结构网格能够更好地拟合起伏地表和地下不规则异常体.由于使用了面向目标的自适应算法,本文能够使用更少的网格达到较高的计算精度.通过对比本文模拟结果与半空间响应和全局自适应算法计算结果,并通过对比使用改进前和改进后加权系数得到的网格剖分结果验证了本文算法的有效性.
相似文献大地电磁勘探方法由于其成本低、施工简单、探测深度广等优点,广泛应用于矿产资源普查、油气勘探和深部构造研究等领域.如何提高大规模三维大地电磁数值模拟的精度和效率一直是研究热点.本文基于空间-波数域方法,实现了基于Lorenz规范的空间-波数域三维大地电磁数值模拟.基于二次场计算原理,引入Lorenz规范,将Maxwell方程组转化为关于二次场矢量位的亥姆霍兹方程;利用水平方向二维傅里叶变换,将空间域三维偏微分方程转换为多个波数下相互独立的常微分方程,方程采用二次插值有限单元法计算,得到定带宽线性方程组,方程计算量小、并行性好,采用追赶法求解,提高了算法效率;引入压缩算子,用迭代法逐次逼近真实解.充分利用了空间-波数域方法数值精度高、内存需求少、效率高的特点.设计棱柱体模型验证了算法的正确性、分析了算法的收敛性,说明算法对不同频率、不同电导率对比度模型均具有很好的适应性.利用Dublin(DTM1)模型进行三维大地电磁数值模拟,结果表明:在满足精度要求的前提下,空间-波数域算法比空间域算法占用内存少、耗时短;相比基于Coulomb规范的空间-波数域算法,基于Lorenz规范的空间-波数域方法耗时更短、占用内存更少,效率提高至少3倍以上,体现了新方法的优势.
相似文献