共查询到20条相似文献，搜索用时 171 毫秒
1.
现有海洋可控源电磁三维数值模拟方法大多基于电导率各向同性介质理论,不能模拟海底地层电导率各向异性的实际情况.本文给出了电导率各向异性三维介质中电性源海洋可控源电磁二次电场的边值问题以及相应的变分问题,采用长方体单元对研究区域剖分,将场分量定义在剖分单元的边上,利用矢量有限单元法求解变分问题,实现了电导率任意各向异性海洋可控源电磁三维矢量有限元数值模拟.这个新的正演方法可以计算电导率任意各向异性三维地电模型的海洋可控源电磁响应,基于二次场矢量有限元法直接求解电磁场,避免了传统有限元方法可能遇到的伪解问题和难于处理电场法向分量不连续的问题,提高了数值模拟计算精度.一维电导率各向异性模型电磁场数值解与解析解吻合得相当好,无论在源附近还是远离源处相对误差均不超过1%.电导率各向异性二维模型的计算结果与已有文献采用的非结构有限元模拟结果十分吻合.三维地电模型数值模拟结果显示,电导率各向异性张量电导率主轴分量和欧拉角对不同装置海洋可控源电磁响应均有着明显的影响. 相似文献
2.
本文实现了2.5维电导率正交各向异性海洋可控源电磁等参有限元数值模拟.利用傅里叶变换导出了电导率正交各向异性2.5维海洋可控源电磁法波数域电磁场耦合方程,采用伽里金加权余量法推导了相应的有限元方程;采用任意四边形单元对研究区域进行剖分,在单元中进行双二次插值,将有限元方程化为线性代数方程组;最后,求解线性方程组并进行反傅里叶变换获得空间域电磁场值.这个方法可以模拟海底起伏地形条件下地下任意形状电导率正交各向异性的复杂模型.与一维模型的数值模拟结果对比表明,电磁场数值解与解析解吻合.二维模型的计算结果与二维自适应非结构有限元模拟结果也吻合.水平海底二维地电模型考察了不同各向异性系数对海洋可控源电磁响应的影响特征.海底起伏地形地电模型的数值结果表明,电导率各向异性对海洋可控源电磁响应影响明显,有可能淹没海底地形和高阻油气藏引起的异常. 相似文献
3.
4.
三维层状场地的精确动力刚度矩阵及格林函数 总被引：3，自引：1，他引：3
本文对Wolf二维层状场地精确动力刚度矩阵进行推广,给出了三维层状场地的精确动力刚度矩阵。刚度矩阵具有对称的特点,且因刚度矩阵是精确的,计算结果不受土层单元厚度的影响,可以大大提高计算效率。文中对刚度矩阵进行了数值验证。利用三维层状场地动力刚度矩阵,计算分析了基岩上单一土层场地的动力响应。最后作为动力刚度矩阵的另一重要应用,给出了表面或埋置矩形均布荷载或集中荷载的动力格林函数计算方法。 相似文献
5.
直流电阻率数值模拟对于实测电阻率数据反演、解释地下真实电阻率结构具有重要意义.为解决复杂地电模型剖分以及处理截断边界时面临的困难,本文提出一种新的三维自然单元-无限元耦合法.该算法在自然单元区域灵活布设自然节点以填充复杂模型,并引入无限单元替代截断边界条件以促使刚度矩阵与电源位置无关.首先推导了稳定电流场基本方程;然后介绍了二、三维自然单元法及其形函数的构建过程;再引入三维多向映射无限单元并实现自然单元-无限元的有效耦合;最后通过4个数值模型,结合传统有限单元法、有限元-无限元耦合法以及COMSOL Multiphysics软件的计算结果,验证了算法的正确性及有效性.数值结果表明,本文所提算法满足精度要求,能有效减小计算区域,无需二次更新刚度矩阵,对复杂模型适应性强. 相似文献
6.
7.
系统地论述了用有限单元法研究复杂地形条件下三维直流电阻率的正演计算技术.首先给出了三维构造中点源电场的边值问题以及相应的变分问题;然后利用有限单元法求解变分问题,采用四面体单元对研究区域进行剖分,在单元中进行三线性函数插值,将变分方程化为线性代数方程组;最后,考虑到节约计算时间,利用对称超松弛顸条件共轭梯度迭代算法求解大型线性方程组,得到了各节点的电位值,进而计算出地表的视电阻率.通过理论模型的计算检验了算法的可行性之后,给出了几种常见纯地形异常的数值模拟结果和一个组合模型的计算结果,其研究工作为研究三维直流电阻率反演奠定了基础. 相似文献
8.
复杂二维/三维大地电磁的有限单元法正演模拟策略 总被引：1，自引：0，他引：1
复杂二维和三维大地电磁模型的正演数值模拟具有一定的挑战性。对于复杂的二维和三维大地电磁正演问题,我们采用有限单元法进行求解。有限单元法最后形成一个线性方程组,系数矩阵是大型稀疏的带状对称复系数矩阵,并且其条件数远大于1,为严重病态矩阵,求解其对应方程组会遇到很多困难。不完全LU分解处理的Bi-CGSTAB迭代方法可用于该线性方程组的求解,并且具有速度快、精度高和稳定性好等优点;为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分;在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间。通过对二维和三维模型电磁响应的计算,验证了算法的正确性。 相似文献
9.
本文提出一种空间波数混合域磁异常场三维数值模拟方法.该方法利用磁位三维空间域积分为卷积的特点,沿水平方向进行二维傅里叶变换,把空间域磁位满足的三维积分问题转化为不同波数之间相互独立的垂向一维积分问题.保留垂向为空间域,优势之一在于便于浅层单元剖分可适当加密,随着深度增加,单元剖分适当稀疏,可以准确模拟任意复杂地形和磁性体的磁异常,兼顾了计算精度与计算效率;优势之二在于一维积分垂向可离散为多个单元积分之和,每个单元采用二次形函数表征磁化强度,可得出单元积分的解析表达式,计算精度高、效率高.该方法充分利用一维形函数积分的高效和高精度、快速傅里叶变换的高效性及算法高度并行性,实现了磁异常场高效、高精度的数值模拟.设计棱柱体模型,将模型解析解与空间波数混合域法的数值解对比,结果表明该方法计算精度高、效率高.设计了组合棱柱体复杂模型,对比分析了标准FFT扩边法与Gauss-FFT法的计算精度与计算效率,总结了标准FFT的扩边系数选取策略.针对任意复杂地形条件下的磁异常模拟问题,本文提出一种适用于起伏地形条件下的磁异常场快速计算方法,并对其有效性进行了验证. 相似文献
10.
在多源直流电阻率法有限元三维数值模拟中,传统混合边界条件由于刚度矩阵与源点位置相关,无法实现线性方程组的快速回代求解,目前常用齐次边界条件或无限元边界进行替代,虽然实现了快速回代求解,但同时也降低了数值模拟的精度.为了实现快速回代求解,并确保数值模拟的计算精度,本文提出了一种近似边界条件方法.其核心思想是将与场源位置相关的边界系数矩阵从刚度矩阵中分离出来,使得分离后的刚度矩阵与场源位置无关.并将边界系数矩阵与边界处一次电场向量的乘积移到线性方程组右端源项中,当场源位置发生改变时,只需要重新计算右端源项就可以实现快速回代求解.理论模型数值计算表明,在水平地形条件下,本文边界条件数值精度优于混合边界条件;在起伏地形条件下,与齐次边界条件相比,本文边界条件数值结果与混合边界条件吻合度更高. 相似文献
11.
三维陆地可控源电磁法有限元快速正演的主要瓶颈在于多频率大型稀疏方程组求解问题.本文引入一种基于模型降阶的Krylov子空间投影算法,推导了有限元刚度矩阵的模型降阶形式,构建了频率域传递函数;采用标准正交向量序列,构建一个远远小于有限元刚度矩阵维度的矩阵,该矩阵与频率无关,通过一次模型降阶即可实现多频点有限元方程快速求解.采用基于电场的变分方程,加入散度校正条件,以消除伪解;引入伪δ函数,消除了源点的奇异性,可适用于复杂背景模型三维有限元数值模拟,并为多源的求解奠定了基础;以层状介质模型解析解为标准,通过和基于Pardiso直接求解器的有限元算法（3DFEM）进行比较,模型降阶法计算时间小于前者的1/10,平均相对误差在1.72%,在满足精度要求下,实现了高效率三维有限元数值求解;分别设计了横向高低阻模型和纵向高低阻模型,分析了从近区到远区电场和卡尼亚视电阻率的变化规律,假极值的表现特征,阴影效应的影响等,从而也验证了该算法的正确性.最后,建立了一个地层陷落柱模型,通过模型降阶有限元正演模拟,发现视电阻率断面图在陷落柱上方出现"凹陷",与模型设计吻合,表明该算法对复杂地层模拟具有同样的适用性. 相似文献
12.
13.
L.A. Tabarovsky M.M. Goldman M.B. Rabinovich K.-M. Strack 《Journal of Applied Geophysics》1996,35(4):261-284
Understanding, using, or eliminating three-dimensional (3-D) effects in electromagnetic methods of geophysics are critical requirements. Numerous achievements in 3-D modeling sometimes give the impression that they are widely available today in geophysical practice. This is not necessarily true. Existing 3-D modeling packages prove that we know how to perform 3-D modeling. However, the computer resources and costs involved make the practical application of 3-D EM modeling in geophysical applications very limited.A practical compromise, or even alternative, is represented by 2.5-D modeling characterized by the use of a 3-D source in a 2-D medium. This combination allows one to mathematically describe the related boundary value problem as a sequence of independent two-dimensional problems. The typical technique leading to such a split formulation is Fourier analysis. That is why the individual terms of a split solution are often referred to as harmonics.Although each independent problem is two-dimensional, the algorithmic implementation of finite differences or integral equations for the higher harmonics has some specific features not present in the classical 2-D cases. In this paper, a hybrid scheme consisting of a combination of the finite difference technique with the integral equation approach for transient fields is described. Evaluation of algorithm accuracy is presented and a transient logging technique application is considered. The algorithm is fast and easily implemented on a personal computer 相似文献
14.
15.
在前人工作的基础上,本文推导了电导率任意各向异性介质的海洋可控源电磁三维谱元法正演方程.采用一次场/二次场分离算法结合混合阶矢量基函数,可以有效避免源点的奇异性的影响,从而提高数值解的精度.采用任意六面体单元离散研究区域,有利于模拟复杂地形和地电结构.利用不完全LU分解的Induced Dimension Reduction(IDR(s))迭代算法求解线性方程组,有效地提高了求解的效率.设计典型的地电模型进行正演计算,并将计算结果与有限元解进行对比,对比结果表明本文提出的基于混合阶矢量基函数的海洋可控源电磁三维谱元数值模拟算法是正确的、有效的.本文算法具有良好的通用性,可推广用于电导率呈任意各向异性的陆地电磁、井中电磁等数值模拟研究. 相似文献
16.
高精度及高效频率域声波数值模拟的关键在于高效求解声波方程经离散化后得到的大型稀疏线性方程组.该方程组系数矩阵具有很强的稀疏性,非对称性和非正定性等特征,常用的迭代算法难以准确、高效地求解.为了改善数值模拟迭代算法的收敛性与稳定性,在算法基础上添加预条件算子是求解该类方程的常用方案.本文基于以上思路,引入正则化技术来构造合适的预条件算子,提出正则化预条件迭代算法,以加速求解方程组.通过包含有均匀介质和高非均匀度介质(Marmousi)模型的数值模拟实验结果表明:与单独使用迭代算法相比,本文提出的正则化预条件迭代算法在计算量方面仅多了一次矩阵-矢量相乘,内存消耗未增加;同时,基于该算法的数值模拟结果能够满足精度要求,较单独使用迭代法能够有效改善收敛性质,加快收敛速度;而且,在二维模型算例下,与LU分解算法相比,基于该算法的内存消耗大幅下降. 相似文献
17.
18.
Piecewise heterogeneous media that the earth presents are composed of large-scale boundary structures and small-scale volume heterogeneities. Wave propagation in such piecewise heterogeneous media can be accurately superposed through the generalized Lippmann–Schwinger integral equation (GLSIE). Two different Born series modeling schemes are formulated for the boundary–volume integral equation with 2-D antiplane motion (SH waves). Both schemes decompose the resulting boundary–volume integral equation matrix into two parts: the self-interaction operator handled with a fully implicit manner, and the extrapolation operator approximated by a Born series. The first scheme associates the self-interaction operator with each boundary itself and the volume itself, and interprets the extrapolation operator as the cross-interaction between each boundary and other boundaries/volume scatterers in a subregion. The second scheme relates the self-interaction operator to each boundary itself and its cross-interaction with the volume scatterers on both sides, and expresses the extrapolation operator as both the direct and indirect (through the volume scatterers) cross-interactions between different boundaries in a subregion. By eliminating the displacement field from the volume scatterers, the second scheme reduces the dimension of the resulting boundary-volume integral equation matrix, leading to a faster convergence than the first scheme. Both the numerical schemes are validated by dimensionless frequency responses to a heterogeneous alluvial valley with the velocity perturbed randomly in the range of ca 5–20 %. The schemes are applied to wave propagation simulation in a heterogeneous multilayered model by calculating synthetic seismograms. Numerical experiments, compared with the full-waveform numerical solution, indicate that the Born series modeling schemes significantly improve computational efficiency, especially for high frequencies. 相似文献
19.
本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性,文中选用二次场算法,将背景电阻率设置为水平层状且各向异性,场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组,通过Galerkin加权余量法形成大型稀疏有限元方程,采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位;进而,利用滑动平均方法得到二次矢量位和标量位在空间的导数,由此得到二次电磁场;通过一维模型对算法的可靠性进行验证,与此同时,针对实际复杂海洋电磁模型,比较有限元模拟结果与积分方程模拟结果,进一步验证算法精度.若干计算结果均表明,文中算法具有良好的通用性,适用于井中电磁、航空电磁,环境地球物理等非均匀且各向异性介质中的电磁感应基础研究. 相似文献
20.
有限差分法和有限单元法在大地电磁场数值模拟中已经得到了广泛的应用,但其数值结果的精度在很大程度上依赖于网格的离散程度.当模拟起伏地形、弯曲界面等复杂地电模型大地电磁场响应时,常常需要花费大量的时间以便得到较合理的离散网格.无网格局部Petrov-Galerkin法(MLPG)不同于有限差分法和有限元法,其形函数和权函数脱离了网格的束缚.本文详细推导了二维大地电磁场边值问题的弱式形式,并将其离散为局部积分域内的表达形式.通过模拟二维海洋地电模型大地电磁场响应,并与结构网格有限元结果进行对比,验证了本文算法和程序的正确性及精度.设计了一个含有弯曲界面的二维地电模型,讨论了不同离散网格对MLPG无网格法模拟结果的影响,并与结构有限元法结果进行了比较,结果表明MLPG无网格法模拟结果受离散网格影响较小.最后利用MLPG无网格法计算了两个海洋起伏地形模型的大地电磁响应,讨论了海底起伏地形对大地电磁响应的影响. 相似文献