共查询到18条相似文献，搜索用时 140 毫秒
1.
电磁法的三维数值模拟是一个对数值算法和计算机硬件要求都非常高的问题.对常用的微分类方法如有限单元法和有限差分法而言,求解最后所得的大型线性方程组是至关重要的一步,直接影响到正演算法的实用性.如何高效、稳定且准确地解线性方程长期以来一直是被探讨的问题.本文实现了基于线性系统直接求解技术的频率域可控源电磁(CSEM)三维正演.使用交错网格有限体积法(FV)来离散化关于二次电场的Helmholtz方程;使用直接解法取代传统的迭代解法来求解离散线性系统,即对系统矩阵进行完全LU分解,具体通过调用大规模并行矩阵直接求解器(MUMPS)来实现.基于理论模型做了一系列数值实验,首先证明了直接解法的高精度和稳定性,并考察了其内存需求、计算时间和并行可伸缩性等主要计算性能,最后检验了所开发的算法快速模拟多场源CSEM问题的能力以及对常规海洋和陆地CSEM模拟的有效性. 相似文献
2.
本文采用有理函数Krylov子空间模型降阶算法实现了同时求解多频可控源电磁法三维正演响应的快速计算.首先采用基于Yee氏交错网格的拟态有限体积法实现控制方程的空间离散，将任意频率的电场响应表示为关于频率参数的传递函数.采用有理函数Krylov子空间算法求解该传递函数.针对构建m维有理函数Krylov子空间需要求解m次（几十到上百）关于有理函数极点和离散控制方程系数矩阵的线性方程组的问题，本文提出采用单个重复极点的有理函数Krylov子空间模型降阶算法，结合直接法求解器PARDISO，采用Gram-Schmidt方法，只需要1次系数矩阵分解和m次矩阵回代即可实现有理函数Krylov子空间的构建，极大地减少了计算量.针对最优化有理函数极点选取问题，本文根据传递函数的有理函数Krylov子空间投影算法的误差分析理论，引入关于单个重复极点的收敛率函数，通过求解有理函数的最大收敛率直接给出最优化的单个重复极点公式.最终实现了不同发射频率的可控源电磁法三维正演响应的快速计算.分别计算了典型层状模型多发射频率的CSAMT和海洋CSEM的正演响应，通过与解析解的对比验证了本文算法在多发射频率正演的计算精度和计算效率；并通过一个三维海洋CSEM勘探设计最优化发射频率和接收区域选取的例子进一步说明本文算法的优点. 相似文献
3.
本文实现了可控源电磁(CSEM)场三维自适应矢量有限元正演算法,该算法采用非结构四面体单元进行三维网格剖分,能够真实模拟地形起伏和复杂电性异常体.采用一次场和二次场分离的方式计算电磁场响应,能够有效解决有限元模拟中的源点奇异性,提高场源附近电磁场数值精度,其中一次场利用CSEM一维正演算法解析求得,二次场采用矢量有限元方法求得.并利用基于后验误差估计的自适应网格细化算法指导网格细化,以减少人为设计网格导致的误差.通过一维和三维模型的数值模拟,验证了本文算法的有效性:一维模型有限元解与解析解吻合得很好,电磁场振幅相对误差在1%左右,相位差整体小于1°;三维模型有限元解与有限体积解吻合得也很好.模拟了一个含三维倾斜板状异常体的可控源电磁响应,表明了本文算法模拟复杂地电结构电磁场的能力和有效性. 相似文献
4.
在前人工作的基础上,本文推导了电导率任意各向异性介质的海洋可控源电磁三维谱元法正演方程.采用一次场/二次场分离算法结合混合阶矢量基函数,可以有效避免源点的奇异性的影响,从而提高数值解的精度.采用任意六面体单元离散研究区域,有利于模拟复杂地形和地电结构.利用不完全LU分解的Induced Dimension Reduction(IDR(s))迭代算法求解线性方程组,有效地提高了求解的效率.设计典型的地电模型进行正演计算,并将计算结果与有限元解进行对比,对比结果表明本文提出的基于混合阶矢量基函数的海洋可控源电磁三维谱元数值模拟算法是正确的、有效的.本文算法具有良好的通用性,可推广用于电导率呈任意各向异性的陆地电磁、井中电磁等数值模拟研究. 相似文献
5.
大规模地球物理电磁数据的定量解释需要发展高效、稳定的三维正反演算法.本文通过求解离散化的三维电场矢量Helmholtz方程,实现了基于有限体积法的频率域可控源电磁(CSEM)三维正演算法.为模拟具有强电性差异的三维电性介质,该算法采用拟态有限体积法(MFV)对Maxwell方程组进行离散化;另外,为获得稳定、高精度的正演数值结果,采用直接矩阵分解技术来求解离散所得到的大型稀疏线性方程组.对于具有多个发射源的CSEM测量来说,一次矩阵分解结果能够用于同频率下所有场源的正演计算.为降低场源奇异性及边界条件对数值精度的影响,采用虚拟场源校正技术,避免了散射场公式中在构建场源项时所需的大量时间.对于具有多个频率的CSEM的模拟计算,采用分频并行策略来加快三维正演计算.最后,通过与一维层状模型及三维模型的数值结果的对比验证了本文所开发的正演算法对频率域CSEM模拟计算的准确性及有效性,表明该正演算法能够有效应用于三维介质的数值计算.另外,对于多频率CSEM的并行测试结果表明基于分频并行策略的并行计算能够显著地降低正演计算时间. 相似文献
6.
三维反演解释是电磁法勘探发展的重要趋势,而如何提高三维反演的可靠性、稳定性和计算效率是算法开发者们目前的研究重点.本文实现了一种频率域可控源电磁(CSEM)三维反演算法.其中正演基于拟态有限体积法离散化,利用直接矩阵分解技术来求解大型线性系统方程,不仅准确、稳定,而且特别有利于含有大量发射场源位置的CSEM勘探情况;对目标函数的最优化采用高斯牛顿法(GN),具有近似二次的收敛性;使用预条件共轭梯度法(PCG)求解每次GN迭代所得到的法方程,避免了显式求解和存储灵敏度矩阵,减小了计算量.以上这些方法的结合应用,使得本文的三维反演算法准确、稳定且高效.通过陆地和海洋CSEM勘探场景中的典型理论模型的反演测试,验证了本文算法的有效性. 相似文献
7.
在使用电偶极发射源的可控源电磁法(CSEM)勘探中,发射源的方位、长度、形状等对观测数据有重要的影响,然而现有的大部分三维数值模拟方法没有全面地将这些因素考虑进来,很多都只能应对非常简单的场源形态,例如单一方位的点电偶极子,这有可能显著降低模拟结果的准确性.本文实现了基于交错网格有限体积(FV)离散的海洋CSEM三维正演算法,能够模拟形态相对复杂的场源,包括任意方位的有限长直导线和弯曲导线发射源.该算法使用一次场/二次场方法,只需对二次场使用FV法求解,避免了场源的奇异性问题;一次场的计算为一维正演问题,使用准解析法求解,并且只要在计算一次场时考虑复杂的场源形态便可以实现同样场源的三维正演.通过与一维理论模型的解析解对比验证了三维程序的准确性,并针对三维理论模型进行了一系列正演测试,初步考察了场源形态对三维正演结果的影响. 相似文献
8.
直流电阻率数值模拟对于实测电阻率数据反演、解释地下真实电阻率结构具有重要意义.为解决复杂地电模型剖分以及处理截断边界时面临的困难，本文提出一种新的三维自然单元-无限元耦合法.该算法在自然单元区域灵活布设自然节点以填充复杂模型，并引入无限单元替代截断边界条件以促使刚度矩阵与电源位置无关.首先推导了稳定电流场基本方程；然后介绍了二、三维自然单元法及其形函数的构建过程；再引入三维多向映射无限单元并实现自然单元-无限元的有效耦合；最后通过4个数值模型，结合传统有限单元法、有限元-无限元耦合法以及COMSOL Multiphysics软件的计算结果，验证了算法的正确性及有效性.数值结果表明，本文所提算法满足精度要求，能有效减小计算区域，无需二次更新刚度矩阵，对复杂模型适应性强. 相似文献
9.
高精度、快速有效的正演模拟算法是三维电磁正反演的前提.为了提高海洋电磁三维数值模拟的精度和效率，本文提出利用基于Gauss-Lobatto-Chebyshev（GLC）基函数的谱元法进行海洋可控源三维电磁正演模拟.谱元法结合有限元法和谱方法的优点.我们通过应用伽辽金加权残差法离散二次电场矢量亥姆赫兹方程，在单元内选择混合阶GLC多项式的张量积作为高阶矢量插值基函数，在求解大型稀疏线性方程组时利用直接求解器进行快速求解，从而实现了三维海洋可控源电磁快速高精度正演模拟.一维和三维模型正演结果验证了本文算法的有效性和准确性.典型模型的数值结果表明谱元法是一种有效的三维海洋可控源电磁正演数值方法，能在稀疏网格剖分情况下获得精确的海洋电磁正演模拟响应. 相似文献
10.
三维陆地可控源电磁法有限元快速正演的主要瓶颈在于多频率大型稀疏方程组求解问题.本文引入一种基于模型降阶的Krylov子空间投影算法,推导了有限元刚度矩阵的模型降阶形式,构建了频率域传递函数;采用标准正交向量序列,构建一个远远小于有限元刚度矩阵维度的矩阵,该矩阵与频率无关,通过一次模型降阶即可实现多频点有限元方程快速求解.采用基于电场的变分方程,加入散度校正条件,以消除伪解;引入伪δ函数,消除了源点的奇异性,可适用于复杂背景模型三维有限元数值模拟,并为多源的求解奠定了基础;以层状介质模型解析解为标准,通过和基于Pardiso直接求解器的有限元算法（3DFEM）进行比较,模型降阶法计算时间小于前者的1/10,平均相对误差在1.72%,在满足精度要求下,实现了高效率三维有限元数值求解;分别设计了横向高低阻模型和纵向高低阻模型,分析了从近区到远区电场和卡尼亚视电阻率的变化规律,假极值的表现特征,阴影效应的影响等,从而也验证了该算法的正确性.最后,建立了一个地层陷落柱模型,通过模型降阶有限元正演模拟,发现视电阻率断面图在陷落柱上方出现"凹陷",与模型设计吻合,表明该算法对复杂地层模拟具有同样的适用性. 相似文献
11.
本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性,文中选用二次场算法,将背景电阻率设置为水平层状且各向异性,场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组,通过Galerkin加权余量法形成大型稀疏有限元方程,采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位;进而,利用滑动平均方法得到二次矢量位和标量位在空间的导数,由此得到二次电磁场;通过一维模型对算法的可靠性进行验证,与此同时,针对实际复杂海洋电磁模型,比较有限元模拟结果与积分方程模拟结果,进一步验证算法精度.若干计算结果均表明,文中算法具有良好的通用性,适用于井中电磁、航空电磁,环境地球物理等非均匀且各向异性介质中的电磁感应基础研究. 相似文献
12.
从电偶源三维地电断面可控源电磁法的二次电场边值问题及其变分问题出发,采用任意六面体单元对研究区域进行剖分,并且在单元分析中同时对电导率及二次电场进行三线性插值,实现电导率分块连续变化情况下,基于二次场的可控源电磁三维有限元数值模拟.这个新的可控源电磁三维正演方法可以模拟实际勘探中地下任意形状及电性参数连续变化的复杂模型.理论模型的计算结果表明,均匀大地计算的视电阻率误差和相位误差分别为0.002%和0.0005°.分层连续变化模型的有限元计算结果表明,其与对应的分层均匀模型解析结果有明显差异.三维异常体组合模型以及倾斜异常体等复杂模型的有限元计算结果也有效地反映了异常形态. 相似文献
13.
14.
可控源电磁法具有分辨率高及抗干扰能力强等特点，是一种重要的地电磁勘探方法.目前，可控源电磁法的高精度正演计算一直是其核心研究问题之一.传统积分方程法一般采用近似积分公式、简单矩形网格和近似的奇异性体积分计算技术，制约了体积分方程法处理复杂地下异常体的能力，降低了计算精度.针对上述问题，本文基于完全积分公式、四面体非结构化网格和奇异体积分的精确解析解来高精度求解复杂可控源电磁模型的正演响应.首先，从电场积分公式出发，推导了可控源电磁问题满足的积分方程；其次，借助于非结构化四面体网格离散技术，实现了地下复杂异常体的有效模拟.最后，利用散度定理把强奇异值体积分转换为一系列弱奇异性的面积分公式，并通过推导获得了这些弱奇异性的面积分公式的解析解，从而最终实现三维可控源电磁问题的高精度积分求解.以块状低阻体地电模型为测试模型，采用本文提出的积分方程方法获得的数值解与其他公开数值算法解进行对比分析，其对比结果具有高度的吻合性，验证了算法的正确性；同时，设计了球状及复杂地电模型进行算法收敛性测试，进一步验证算法的正确性以及能够处理地下复杂模型的能力. 相似文献
15.
频率域地空电磁探测方法是指在地面布设人工场源,在空中测量电磁场的一种高效的地球物理勘探技术.该方法具有大范围、高通过性、快速测量的优势,尤其适合崎岖山地、沙漠、沼泽、海陆交互带等复杂地貌区域的资源勘查.但是这些地区的地形起伏通常较大,因此分析地形对地空电磁响应的影响具有重要意义.本文利用有限元法对频率域地空电磁响应进行了正演计算,分析了起伏地表条件下的频率域地空电磁响应特征.首先利用傅里叶变换将2.5维问题转化成二维问题,利用伽辽金加权余量法推导了相应的离散有限元方程组.采用任意四边形单元对区域进行不均匀网格剖分,源和异常体附近网格加密处理,保证计算精度,远离目标区域网格逐渐稀疏,模拟无穷远边界,降低对计算资源的要求.在单元内进行插值,将有限元方程组变换为线性方程组,采用总场算法,利用具有一定面积的伪δ函数表达源电流分布,源项近似为分布在以电偶极源为中心的25个节点上.通过求解线性方程组得到波数域电磁响应,再对波数域电磁场响应进行反傅里叶变换从而获得空间域2.5维频率域电磁场值.通过对比2.5维正演结果与均匀半空间解析解,验证了本文算法的精度,同时本文还对地空电磁场与地面电磁场的响应特性进行了对比. 相似文献
16.
根据库伦规范势的定义，推导出关于磁矢量势和电标量势的偏微分方程，为了克服由电流源引起的奇异性和数值模拟计算困难，将电磁总场分解为一次场和二次场，一次场由基于Schelkunoff势函数的一维正演算法得到，二次场由有限元法计算得到，实现了海洋可控源电磁法三维有限元正演算法。通过一维数值模拟实例，验证该算法的计算精度。然后，利用该算法对带海底地形的三层储层模型进行正演，分析了海底地形对海洋控源电磁场各分量产生的影响。 相似文献
17.
为了分析海洋中含金属矿的储油气砂层的成分、结构对于海洋可控源电磁场的影响特征，引入多相微增等效介质模型，结合非结构化网格实现了海洋可控源电磁三维矢量有限元正演.首先对多相微增模型进行了介绍，并分析不同参数影响下的等效电导率变化特征.而后，详细推导了频率域可控源电磁法非结构化矢量有限元方程.接着，利用预处理的IDR （s）迭代算法求解线性方程组.最后设计典型的海洋地电模型，通过正演计算，验证了本文算法的正确性和有效性，同时研究了海洋储油砂层中的金属矿的含量、孔隙度、含水饱和度等参数对于可控源电磁场响应特征的影响. 相似文献
18.
地层宏观电性各向异性会对可控源电磁响应产生重要影响.由于海底地层电性结构常表现为电导率各向异性，若仅对海洋可控源电磁（MCSEM）数据进行常规各向同性反演，有可能无法获得准确的反演解释结果，从而削弱MCSEM技术的可靠性.本文实现了电导率垂直各向异性（VTI）条件下频率域海洋可控源电磁数据三维反演算法.其中，三维正演采用基于二次场控制方程的交错网格有限体积法，并利用直接矩阵分解技术来求解离散所得的大型线性方程组，有利于快速计算多场源的响应.反演采用具有近似二次收敛性的高斯牛顿算法对目标函数进行最优化.最后，对具有VTI电性各向异性特征的盐丘构造模型的MCSEM合成数据分别进行了电导率各向同性和垂直各向异性三维反演，结果表明：各向同性三维反演算法无法对受VTI介质影响的MCSEM数据进行正确的反演解释，而垂直各向异性三维反演能够获得更为可靠的地下电阻率结构和异常体分布，展现出对海底电性各向异性结构更为优良的反演解释能力. 相似文献