共查询到16条相似文献，搜索用时 171 毫秒
1.
有限差分法和有限单元法在大地电磁场数值模拟中已经得到了广泛的应用,但其数值结果的精度在很大程度上依赖于网格的离散程度.当模拟起伏地形、弯曲界面等复杂地电模型大地电磁场响应时,常常需要花费大量的时间以便得到较合理的离散网格.无网格局部Petrov-Galerkin法(MLPG)不同于有限差分法和有限元法,其形函数和权函数脱离了网格的束缚.本文详细推导了二维大地电磁场边值问题的弱式形式,并将其离散为局部积分域内的表达形式.通过模拟二维海洋地电模型大地电磁场响应,并与结构网格有限元结果进行对比,验证了本文算法和程序的正确性及精度.设计了一个含有弯曲界面的二维地电模型,讨论了不同离散网格对MLPG无网格法模拟结果的影响,并与结构有限元法结果进行了比较,结果表明MLPG无网格法模拟结果受离散网格影响较小.最后利用MLPG无网格法计算了两个海洋起伏地形模型的大地电磁响应,讨论了海底起伏地形对大地电磁响应的影响. 相似文献
2.
本文实现了2.5维电导率正交各向异性海洋可控源电磁等参有限元数值模拟.利用傅里叶变换导出了电导率正交各向异性2.5维海洋可控源电磁法波数域电磁场耦合方程,采用伽里金加权余量法推导了相应的有限元方程;采用任意四边形单元对研究区域进行剖分,在单元中进行双二次插值,将有限元方程化为线性代数方程组;最后,求解线性方程组并进行反傅里叶变换获得空间域电磁场值.这个方法可以模拟海底起伏地形条件下地下任意形状电导率正交各向异性的复杂模型.与一维模型的数值模拟结果对比表明,电磁场数值解与解析解吻合.二维模型的计算结果与二维自适应非结构有限元模拟结果也吻合.水平海底二维地电模型考察了不同各向异性系数对海洋可控源电磁响应的影响特征.海底起伏地形地电模型的数值结果表明,电导率各向异性对海洋可控源电磁响应影响明显,有可能淹没海底地形和高阻油气藏引起的异常. 相似文献
3.
为了计算带任意地形的各向异性介质中二维大地电磁响应,本文在非结构化网格的基础上,采用有限体积法,开发了二维大地电磁各向异性正演模拟的新算法.首先,从Maxwell方程出发,推导二维各向异性介质中大地电磁场的边值问题;然后,采用三角网格自动生成技术对求解区域进行非结构化网格剖分,进而构建节点中心控制体积单元,利用有限体积方法,得到求解边值问题的大型稀疏线性方程组;最后,利用Pardiso精确地计算了大地电磁响应值.三个各向异性模型的计算结果表明,本文开发的有限体积算法,不仅能够高精度求解带任意地形的大地电磁电导率各向异性问题,而且对于同一模型,该方法的计算消耗和精度都与有限单元法相当.因此,有限体积法是处理电磁法各向异性问题的一种有效方法. 相似文献
4.
目前,对于可控源电磁法各向异性介质2.5维问题,主要采用一次场、二次场分离的方法消除场源奇异性并降低截断边界对计算区域的影响.该方法数值计算精度高,但是很难适用于复杂地形条件下的数值模拟.针对复杂地形问题,基于总场的有限元方法表现出一定的优越性,然而,这种方法存在场源奇异性问题和截断边界问题.本文采用基于总场计算的方法对带地形的可控源电磁法2.5维各向异性介质进行模拟研究,推导了考虑电导率和介电常数各向异性的2.5维控制方程;引入网格加密-收缩算法降低场源奇异性的影响范围,提升数值计算效率;引入行波分解边界条件降低截断边界的影响;提出任意采样反傅里叶变换方法,快速、高精度地计算出空间域电磁场分量.理论模型数值算例中:首先,验证了本文算法的有效性;其次,对任意各向异性倾角产生的可控源电磁响应规律进行研究;最后,采用山丘模型对各向异性介质电磁场的响应规律进行了模拟和分析. 相似文献
5.
带地形的大地电磁二维正演数值模拟多数基于电性各向同性理论,由于地球内部电性各向异性现象的普遍存在,基于电性各向异性理论研究地形起伏情况下大地电磁二维正演数值模拟就显得非常迫切.本文首先由麦克斯韦方程出发,引入张量电导率,求得一组关于平行走向的电场分量Ex和磁场分量Hx的二阶偏微分方程,使用有限差分法求解出Ex和Hx的近似解,并以此求得其他场分量;其次,引入地形因素,改变变量在网格节点中的排列方式,选择交错排列方式从而给有限差分系数矩阵的最大带宽分配合理的存储空间;最后,使用Weaver的方法解决TM模式下,在地-空分界面垂直于构造走向的一些区域存在不同电导率的问题.通过对带地形的二维电性各向异性结构做正演模拟,研究地形因素对大地电磁响应的影响;以电性各向异性理论为基础,将地形因素引入对实测大地电磁资料的处理中,通过做二维正演拟合和未引入地形因素的结果做对比,说明电性各向异性现象的普遍存在,认识地形因素对观测大地电磁场的影响,为今后分析解释实测大地电磁资料包含地形因素和电性各向异性情况提供理论基础和技术指导. 相似文献
6.
在前人工作的基础上,本文推导了电导率任意各向异性介质的海洋可控源电磁三维谱元法正演方程.采用一次场/二次场分离算法结合混合阶矢量基函数,可以有效避免源点的奇异性的影响,从而提高数值解的精度.采用任意六面体单元离散研究区域,有利于模拟复杂地形和地电结构.利用不完全LU分解的Induced Dimension Reduction(IDR(s))迭代算法求解线性方程组,有效地提高了求解的效率.设计典型的地电模型进行正演计算,并将计算结果与有限元解进行对比,对比结果表明本文提出的基于混合阶矢量基函数的海洋可控源电磁三维谱元数值模拟算法是正确的、有效的.本文算法具有良好的通用性,可推广用于电导率呈任意各向异性的陆地电磁、井中电磁等数值模拟研究. 相似文献
7.
8.
传统三维大地电磁各向异性模拟均是基于规则六面体网格,计算精度有限且较难拟合复杂地质条件.本文采用面向目标自适应非结构矢量有限元法,对三维大地电磁各向异性介质进行模拟.首先从电场双旋度方程出发,利用伽辽金方法建立变分方程;然后利用电流密度连续性条件构建适合大地电磁各向异性问题的加权后验误差估计方法,实现面向目标的网格自适应正演;最后通过典型算例分析各向异性对网格自适应和大地电磁响应的影响特征以及各向异性的识别方法.本文算法能够高精度地拟合起伏地表和任意各向异性介质,适用于分析复杂地电条件大地电磁响应特征,为提高大地电磁资料解释水平提供了理论基础. 相似文献
9.
为提高大地电磁正演计算速度,开展了基于多重网格有限元法的大地电磁二维正演模拟计算研究.将稳定双共轭梯度算法作为多重网格法的细网格松弛迭代算法,插值算子采用完全加权算子,限制算子设计基于网格单元面积率,使多重网格法更适于求解大型复系数方程组.二维均匀半空间模型、低阻体模型和高阻体模型的大地电磁正演模拟结果表明:当计算量较小时(网格剖分数量少),多重网格法在计算效率方面并未有优势,网格剖分数量较大时,多重网格有限元算法在收敛速度方面的优势明显,多重网格有限元法的大地电磁正演精度优于一般数值算法.这为三维多重网格有限元的大地电磁正演研究奠定了基础. 相似文献
10.
为了易于模拟野外复杂地形和地下任意形状地电体模型,将有限元单元网格设计为三角单元;并考虑到野外实际勘探中,地球介质的电性参数均是连续变化的情况,单元内的场值和电性参数被设计为双线性变化;推导出二维起伏地形条件下大地电磁法有限元数值模拟算法;根据单元节点主场值和线性插值形函数间的关系,计算出单元节点的辅助场值;在二维起伏地形情况下,定义TE、TM模式视电阻率和阻抗相位.4个模型的计算的结果与解析法的均方根误差小于1%,地形模拟与前人的计算结果相符,模拟倾斜界面异常体,能有效的反映出其异常形态. 相似文献
11.
经过半个多世纪的发展,国内外利用大地电磁法研究地球内部电性结构取得了令人瞩目的成就,这些研究成果多数是基于电性各向同性理论.然而地球内部普遍存在电性各向异性现象,地壳和上地幔中存在的电性各向异性是地电模型、地下结构和构造模型间一个重要的联系因素.本文首先由麦克斯韦方程出发,引入张量电导率,根据二维电性各向异性结构的特点,得到一组关于平行走向的电场分量Ex和磁场分量Hx的偏微分方程.使用有限差分法求解偏微分方程,求出Ex和Hx的近似解,并以此求得其它场分量;随后,通过对普通及特殊的二维电性各向异性结构做正演模拟,研究其对观测大地电磁场的影响,从而认识在普通及某种特定地质条件下的电磁传播特性,为其后对大地电磁实测资料的处理解释奠定理论基础;最后,以本文的研究成果为基础,将电性各向异性理论引入对实测大地电磁资料的处理解释中,通过对新疆某地的大地电磁资料做二维正演拟合解释,说明了电性各向异性现象的普遍存在,也验证了理论的正确性及算法的实用性,为今后分析解释大地电磁资料中的电性各向异性现象提供理论依据和技术指导,并开拓了对大地电磁实测资料处理的思路和方法. 相似文献
12.
本文以基于改进BISQ模型的二维双相各向同性介质一阶速度-应力方程为基础,推导出了曲线坐标系下对应的方程,然后采用低频散、低耗散的同位网格MacCormack有限差分法来离散方程,并采用紧致的单边MacCormack差分格式结合牵引力镜像法来施加自由地表边界条件,实现了地震波场数值模拟.曲线网格有限差分法采用贴体网格来描述自由表面,地表的网格线紧贴地形,避免了台阶近似造成的数值散射.数值模拟结果表明,在双相介质起伏自由地表和分界面处,各类波型复杂的反射透射规律可以清晰展现,曲线网格有限差分法可以精确地解决地震波在含起伏地表的双相各向同性介质中的传播问题. 相似文献
13.
本文研究并建立了一种模拟各向异性地层中三维电性异常体电磁响应的积分方程算法.首先讨论了并矢Green函数及其相关积分的计算,将水平层状各向异性地层中的电场并矢Green函数分解成含有奇异项的直达波与非奇异的来自各个层界面的反射和透射波两个部分,再应用等效体积单元和表面积分技术对积分方程的奇异核进行离散化处理以便提高离散方程的精度.然后为了节省计算机内存以及计算时间,引入基于Krylov子空间的迭代算法求解积分方程的离散化矩阵方程.最后通过与现有文献中的结果作对比从而检验了所述算法的有效性,并结合具体算例考察分析了地层的各向异性对三维电性异常体电磁响应的影响特征和规律. 相似文献
14.
A magnetotelluric finite-element modelling algorithm is developed, which is capable of handling three-dimensional conductive and magnetic anisotropic anomalies. Different from earlier three-dimensional magnetotelluric anisotropic modelling methods, the algorithm we presented has taken the magnetic anisotropy into consideration. The variational equations are produced by the Galerkin method and the governing equations are solved using a hexahedral vector edge finite-element method. The accuracy of this algorithm is firstly validated by comparing its solutions with the results of finite-difference method for a three-dimensional conductive arbitrary anisotropic model, and then validated by comparing with analytical solutions for a one-dimensional magnetic model. The responses of four kinds of models under different conditions are studied, and some conclusions are obtained. It shows that for materials with a high magnetic permeability, its influence on magnetotelluric responses cannot be ignored in some circumstances. Especially, if the magnetic susceptibility is exceptionally high, it may really distort the apparent resistivities of lower resistive anomalies. These conclusions are also beneficial for magnetotelluric survey. 相似文献
15.
In this paper, we deduced the corresponding first-order velocity–stress equation for curvilinear coordinates from the first-order velocity–stress equation based on the modified Biot/squirt model for a two-dimensional two-phase medium. The equations are then numerically solved by an optimized high-order non-staggered finite difference scheme, that is, the dispersion relation preserving/optimization MacCormack scheme. To implement undulating free-surface topography, we derive an analytical relationship between the derivatives of the particle velocity components and use the compact finite-difference scheme plus a traction-image method. In the undulating free surface and the undulating subsurface interface of two-phase medium, the complex reflected wave and transmitted wave can be clearly recognized in the numerical simulation results. The simulation results show that the curvilinear-grid finite-difference method, which uses a body-conforming grid to describe the undulating surface, can accurately reduce the numerical scattering effect of seismic wave propagation caused by the use of ladder-shaped grid to fit the surfaces when undulating topography is present in a two-phase isotropic medium. 相似文献
16.
瞬变电磁法应用广泛,三维数值模拟是研究复杂地质模型异常响应规律的重要技术手段之一,也是反演的基础.目前瞬变电磁数值模拟的不足主要有两个方面:第一,场源是在地表水平、浅层介质均匀的条件下计算的,限制了应用范围;第二,地下边界采用Dirichlet边界条件,导致计算空间很大,耗时较长.针对上述问题,在三维正演时,场源采用有限长细导线模型,在Maxwell有源差分方程中直接加入电流密度进行计算.在地表面加入空气层,避免了复杂的向上延拓计算,也可以对地形影响下的响应规律进行分析.在空气边界和地下边界均采用CPML吸收边界条件,并改进了CPML的参数分布,能够吸收空气介质和大地介质中的低频电磁波而反射误差极小,在满足计算精度的条件下可以有效减小节点数量.对循环迭代方法进行优化,将计算域、CPML区域和场源的空间循环统一转化为矩阵方式,加快了计算速度,但是空间消耗增大了约4~5倍.采用三维有限差分正演算法对均匀半空间模型、层状模型和地形模型进行了计算,并与解析解进行了对比验证. 相似文献