为了提高广域电磁法三维正演精度和效率, 本文提出采用基于Gauss-Lobatto-Legendre(GLL)多项式的谱元法进行广域电磁法三维正演模拟.首先从麦克斯韦方程组出发, 推导了二次场满足的电场双旋度控制方程, 利用伽辽金加权残差法将微分形式的边值问题转换为积分弱形式, 再通过单元剖分和高阶正交基函数插值对全域问题进行离散, 最后通过求解大型线性方程组得到全局的数值解.层状介质模型验证了算法的正确性和精度, 三维地电模型分析了算法稳定性以及广域电磁法响应特征.研究表明谱元法是进行广域电磁法三维正演模拟的有效方法, 具有高精度、低网格依赖性等优势.
相似文献大地电磁正演常采用电场双旋度方程,由于工作频率低、空气电导率为零,电场双旋度方程含有空解空间,导致离散化的线性方程组严重病态、具有高的条件数.为降低大地电磁场双旋度正演方程的病态性,本文提出了一种基于空气与地球分离策略的大地电磁混合正演公式.首先,在空气层中将双旋度算子转换为Laplace算子,在地下利用电导率大于零的特征保留双旋度算子,建立空气与地下分离的混合公式系统,以便实现混合公式系统对应的离散化线性方程组条件数的大幅度降低.然后,在空气与地下分界面上,施加电磁场的切向与法向连续性条件实现空气与地下区域的耦合与信息交换,以保证电磁场分布的唯一性.再则,利用非结构化四面体网格离散技术来高精度模拟任意起伏地形和复杂地下三维地电结构,并结合节点有限元与棱边有限元法,高精度求解该混合公式系统.最后,利用国际标准大地电磁测试模型来验证该空气与地球分离算法的正确性,并展示其有效降低原有电场双旋度方程的病态性的能力.本文提出的空气与地球分离思想为寻求大地电磁最优化控制方程开启了新的求解思路,预计会产生更多、性能更佳的大地电磁混合公式,以便服务于大地电磁正反演技术的快速发展.
相似文献正演模拟是电磁数据反演的基础, 其计算速度与精度一直是制约电磁反演的两大核心问题.在三维电磁正反演中, 传统方法通过加密网格或增加插值基函数阶数提高计算精度, 但由此也降低了计算效率, 制约了三维电磁反演的实用化.因此, 如何实现大尺度模型高精度快速正演是目前电磁三维正反演中亟需解决的问题.本文将多尺度有限元法应用到麦克斯韦方程求解中.我们首先在粗网格尺度上构建满足局部特性微分算子的多尺度基函数, 进而在粗网格尺度上对原问题进行求解, 通过建立粗细两套网格间场的映射关系, 在未知数较少的粗网格上实现电磁问题求解之后, 利用粗细两套网格间场的映射关系获取细网格上电磁场响应, 由此可以在保证计算精度前提下快速获取不同尺度电磁场正演响应, 计算速度得到很大提高.此外, 本文还基于八叉树思想进行网格优化, 进一步改善三维正演效率.我们通过对典型地电结构进行多尺度有限元正演模拟并与传统有限元结果对比验证算法的有效性.最后, 我们通过模拟加拿大Voisey's Bay卵形体镍铜硫化矿区航空电磁响应以检验本文算法模拟地下复杂异常体的能力.
相似文献在经典位场理论中,许多简单形体位场异常难以通过积分得到全空间的解析式.圆柱体是一类很重要的理论模型体,常用于模拟圆柱状地质体或非地质体(如管线),但目前还不能用解析公式正演有限长圆柱体在三维空间里的磁异常,而多是采用近似简化为有限长磁偶极子或线模型代替.对于有限长圆柱体,特别是半径相对于上顶埋深较大时,这种近似的误差不可忽略.本文利用共轭复数变量替换法,推导出有限长圆柱体在全空间的引力位一阶、二阶导数,利用Poisson关系得到磁异常正演公式,进而利用有限长圆柱体磁异常正演公式求解管状体的磁异常,得到不同磁化方向、不同大小的管线产生的磁场的特征,并将其推广到截面为椭圆的情况.最后通过模拟计算定量给出了将圆柱体近似为线模型的条件.
相似文献三维磁异常快速反演是磁法勘探进行地质解释的一项重要内容.本文在方向tilt梯度基础上进行导数处理,并进行合理改进,得到了6个改进的方向tilt梯度导数,丰富了不同方向上的磁异常信息.对常规欧拉反褶积方程进行三个方向求导,所构成的方程组与6个改进的方向tilt梯度导数进行结合,从而构建出了两个基于方向tilt梯度导数的快速反演方程组(方向tilt-Euler法),用于反演场源的位置参数.另外,将三个方向导数下的欧拉反褶积进行平方和及均方根处理,得到了场源构造指数的求解公式.单一模型试验表明,6个改进的方向tilt梯度导数异常基本不受磁化方向影响,可从不同方向上反映磁异常信息,方向tilt-Euler能够准确地反演出场源的位置及构造指数.叠加模型试验证实了方向tilt-Euler法相对于tilt-Euler法,具有反演解聚集度更强、连续性更好、反演准确度更高、受噪声影响更小及无虚假反演解等优点.将本文方法应用于内蒙古塔木素航磁资料处理中,获得了丰富的地下磁源分布信息,为了解该地区隐伏断裂构造及岩体分布等提供了有效依据.
相似文献本文提出一种空间波数混合域磁异常场三维数值模拟方法.该方法利用磁位三维空间域积分为卷积的特点,沿水平方向进行二维傅里叶变换,把空间域磁位满足的三维积分问题转化为不同波数之间相互独立的垂向一维积分问题.保留垂向为空间域,优势之一在于便于浅层单元剖分可适当加密,随着深度增加,单元剖分适当稀疏,可以准确模拟任意复杂地形和磁性体的磁异常,兼顾了计算精度与计算效率;优势之二在于一维积分垂向可离散为多个单元积分之和,每个单元采用二次形函数表征磁化强度,可得出单元积分的解析表达式,计算精度高、效率高.该方法充分利用一维形函数积分的高效和高精度、快速傅里叶变换的高效性及算法高度并行性,实现了磁异常场高效、高精度的数值模拟.设计棱柱体模型,将模型解析解与空间波数混合域法的数值解对比,结果表明该方法计算精度高、效率高.设计了组合棱柱体复杂模型,对比分析了标准FFT扩边法与Gauss-FFT法的计算精度与计算效率,总结了标准FFT的扩边系数选取策略.针对任意复杂地形条件下的磁异常模拟问题,本文提出一种适用于起伏地形条件下的磁异常场快速计算方法,并对其有效性进行了验证.
相似文献磁异常反演是获取地下场源磁化率分布的重要手段之一,在地球勘探中扮演着重要角色.在磁异常反演中,对比光滑反演,稀疏反演的结果具有边界分明,物性参数分布集中的特点,更符合实际情况.针对稀疏反演,本文首先构建了具有代表性的基于L1范数目标函数,利用交替方向乘子算法可分离凸函数的特点,将极小化L1范数的优化问题分解为一系列的子问题,通过对子问题求解获得原问题的解;为了增强交替方向乘子算法的适应性,本文结合广义软阈值函数将交替方向乘子法推广于Lp(0 < p < 1)范数的反演中.为了验证本文提出的算法的有效性,采用了三种常规模型进行模拟实验.与基于L2范数的反演算法进行实验对比,结果表明,本文算法得到了边界清晰,磁化率分布更集中的反演结果.最后,将基于交替方向乘子算法的L1和Lp(0 < p < 1)范数的反演应用到青海省尕林格铁矿保护区获得的实际磁异常数据中,获得了较为符合实际地质情况的稀疏反演结果.
相似文献