共查询到20条相似文献，搜索用时 15 毫秒
1.
三角网格谱元法地震波场数值模拟 总被引：5，自引：0，他引：5
谱元法结合了有限元法的灵活性和谱方法的指数收敛性,高效且高精度,是近年来发展的一种重要的地震波场数值模拟方法.经典的谱元法采用四边形(六面体)网格,利用一维Gauss-Legendre-Lobatto(GLL)积分的张量积得到对角的质量矩阵,以大大提高计算效率,但是四边形(六面体)网格不能够灵活地刻画复杂的几何模型的弯曲界面.为此,在谱元法中引入三角形(四面体)网格到二维(三维)是十分必要的.不同于经典的谱元法,在非结构化网格中不能使用GLL积分的张量积,使得非结构化网格的谱元法的实现存在着诸多的困难.目前,比较流行的三角网格谱元法,通过使用KoornwinderDubiner(KD)正交多项式,并正交化这些KD多项式构建基函数,同时利用重合的插值节点和积分节点以获取对角的质量矩阵;它所使用的积分点为优化的点集——Fekete点,且这些积分点能与四边形网格完全耦合.相比于四边形,三角网格谱元法能显著提高复杂模型的描述能力,对起伏地表模型有很大优势.本文引入高效的最佳匹配层(PML)吸收边界条件,并通过数值试验将三角网格谱元法与经典的谱元法进行对比研究.相比于经典的谱元法,三角网格谱元法显著缺点为较低的计算精度.对于7阶谱元,为了能够精确地模拟面波,三角网格谱元法需要在每个最短的面波波长内至少有11个采样点,然而经典的谱元法仅需4个采样点,并且前者所需的内存量约为后者的5.5倍. 相似文献
2.
为了提高广域电磁法三维正演精度和效率,本文提出采用基于Gauss-Lobatto-Legendre(GLL)多项式的谱元法进行广域电磁法三维正演模拟.首先从麦克斯韦方程组出发,推导了二次场满足的电场双旋度控制方程,利用伽辽金加权残差法将微分形式的边值问题转换为积分弱形式,再通过单元剖分和高阶正交基函数插值对全域问题进行... 相似文献
3.
各向异性是地电异常解释中不可忽视的因素,广泛存在于裂隙或层理发育的地质环境中.本文针对任意各向异性条件下直流电阻率法三维正演问题进行研究,结合非结构谱元法建立模拟算法,充分利用谱方法的指数收敛性以及非结构有限元对地形和复杂异常体刻画能力,提高计算精度和效率.通过灵活的四面体网格剖分和高阶谱插值,实现了复杂介质任意各向异性模型电阻率响应的高精度数值模拟.我们首先通过层状各向异性模型验证本文非结构谱元法的计算精度,进而我们以半空间中立方体模型为例分析各向异性对电阻率响应的影响特征,并通过计算针对不同各向异性参数的视电阻率极性图,探究地下介质各向异性特征识别方法.最后,我们针对典型的山脊模型计算和分析存在地形效应条件下各向异性电流场分布及视电阻率特征.模型计算结果表明基于四面体网格的谱元法模拟带有复杂地形和异常体的任意各向异性模型具有很高的计算精度.本文的研究成果将在推进电阻率方法用于解决裂隙及层理等环境和工程地质问题中发挥积极作用. 相似文献
4.
网格剖分程度直接影响着地震波正演数值模拟的计算精度及其计算量.以均匀介质模型为例,分析不同网格大小对波场模拟精度和计算量的影响,得出精细化网格剖分是实现高精度地震波正演模拟的有效方法,然而其计算量较大.以均匀倾斜介质模型为例,探讨倾斜地层网格剖分问题,数值实例分析不同震源频率对不同网格剖分方案引起的波场传播精度的影响,... 相似文献
5.
介绍一种基于离散粒子理论地震波传播数值模拟的网格剖分计算方法.根据离散粒子理论,将研究区域划分为由一系列相互作用的粒子组成的正六边形网格,这些粒子在它们的接触点处发生相互作用,并用Hooke定律和Newton定律描述.为解决六边形网格带来的网格交错而难以计算以及波场输出问题,将横向网格进行加密,加密处赋予假想的粒子,输出波场时选取偶数行偶数列点或奇数行奇数列点的波场值.均匀介质和层状介质模型的数值模拟结果表明,该网格剖分计算方法能够将离散粒子理论用于模拟弹性波在非均匀各向同性介质中地震波的传播. 相似文献
6.
7.
波动方程的数值求解是地震波正反演的重要环节,而数值算法的计算精度直接关系到地震波的模拟结果和成像质量.当前,谱元法由于同时具备有限元法的网格灵活性与谱方法的高精度性已被成功应用于不同尺度模型中的地震波模拟.然而,常见的Legendre谱元法在求解地震波运动方程时采用Gauss-Lobatto-Legendre(GLL)数值积分计算质量矩阵所包含的积分项,由于GLL数值求积无法对积分项精确估计,从而造成谱元法精度损失.针对谱元法精度上的不足,本文提出一种优化算法用于提升其精度.首先构造关于GLL数值求积积分权与质量矩阵对角线元素精确值的最小二乘目标函数,然后利用共轭梯度法求解目标函数得到优化权系数,该权系数能减小质量矩阵的离散误差最终提高谱元法的计算精度.通过数值频散分析、数值算例证实了本文给出的优化算法用于提升谱元法数值模拟精度的可行性和有效性. 相似文献
8.
利用谱元法的规则六面体单元进行网格剖分时,界面起伏较大处会出现阶梯状网格而导致模拟时产生数值散射.为消除阶梯状网格对起伏界面地震动模拟计算的影响,本文基于TrueGrid软件编写了应用程序,提出了起伏界面处六面体单元网格剖分方式,通过该程序可快速建立起伏界面处均匀的六面体网格模型.本文采取了删除四纵列拐角单元、删除一纵... 相似文献
9.
在交错网格有限差分算法中,模型网格剖分原则与正演计算效率密切相关.当模型存在小型非均质体或者低降速层等情况,为保证精度,满足稳定性条件,需缩小网格步长,导致局部过采样,计算效率低下.为保证模拟精度的同时保持高计算效率,通常采用变化的空间网格与时间步长相结合的高阶有限差分模拟方法对波场进行模拟.然而,时空双变算法存在着交错网格固有缺点,在模拟非均匀性较强的复杂介质波场传播时,需对介质参数进行平均或内插.同时,该算法在空间与时间上的变网格实现极为复杂.为压制变网格引起的虚假反射,提升模拟精度和计算效率,本文在时空双变网格算法的基础上,采用旋转差分角度的方式,提出了旋转时空双变交错网格算法.该方法既保留了旋转网格和双变网格的优势,又简化了时空双变算法流程,更利于推广和应用. 相似文献
10.
《应用地球物理》2017,(3)
航空电磁法具有高采样率探测特征,传统数值模拟方法为保证航电响应计算精度,特别是针对复杂地质体(如电各向异性),需要对模拟域进行精细剖分,导致计算工作量庞大。因此本文采用谱元法对航空电磁各向异性响应进行三维正演模拟,该方法结合谱方法和有限元法双重优势,它与谱方法类似采样高阶基函数代替有限元中的线性插值,同时具有有限元拟合边界的灵活性。因此该方法能改善离散网格内部的数值模拟精度,减弱数值模拟结果对网格的依赖性,实现航空电磁各向异性响应的高精度计算。首先,本文将旋转张量的各向异性电导率引入Maxwell方程,并给出GLL谱元基函数给出电场表示形式,采用伽辽金加权余量法形成基于谱元法大型线性方程组获得电磁响应,以VMD发射源为例进行航空电磁响应模拟。其次,本文采用粗物理网格四阶谱元法对任意各向异性半空间的航空电磁响应进行计算,通过与一维半解析解对比验证了该算法的高精度性。再次,本文采用不同物理网格和谱元阶数讨论不同各向异性条件下三维异常体的航空电磁响应,得出了不同各向异性异常体的电磁响应的收敛条件。最后,本文首次基于谱元算法,采用Hr/Hz方式给出了三维各向异性异常体,各向异性围岩以及围岩和异常体同时存在各向异性时航空电磁各向异性的识别方式,这航空电磁各向异性的三维反演和各向异性地质区域测量具有指导意义。 相似文献
11.
《地球物理学进展》2015,(5)
基于时域有限差分(FDTD)的瞬变电磁三维正演已经基本实现,但是复杂三维体的建模还存在问题,从而影响瞬变电磁对于复杂模型响应规律的分析.本文针对复杂模型的几何建模及Yee元胞的网格剖分问题,基于GID软件平台进行了二次开发,编写了"FDTD问题类型",通过编写代码,获得了能够直接应用于FDTD三维正演的网格剖分节点数据.文中利用GID软件对实际采矿巷道、隧道内TBM机模型、以及矿体模型进行建模,并采用FDTD方法进行正演计算.计算结果表明,基于二次开发的GID的可视化建模能够满足目前复杂模型建模的需求,为瞬变电磁FDTD三维正演模拟提供有力的支持. 相似文献
12.
高精度、快速有效的正演模拟算法是三维电磁正反演的前提.为了提高海洋电磁三维数值模拟的精度和效率,本文提出利用基于Gauss-Lobatto-Chebyshev（GLC）基函数的谱元法进行海洋可控源三维电磁正演模拟.谱元法结合有限元法和谱方法的优点.我们通过应用伽辽金加权残差法离散二次电场矢量亥姆赫兹方程,在单元内选择混合阶GLC多项式的张量积作为高阶矢量插值基函数,在求解大型稀疏线性方程组时利用直接求解器进行快速求解,从而实现了三维海洋可控源电磁快速高精度正演模拟.一维和三维模型正演结果验证了本文算法的有效性和准确性.典型模型的数值结果表明谱元法是一种有效的三维海洋可控源电磁正演数值方法,能在稀疏网格剖分情况下获得精确的海洋电磁正演模拟响应. 相似文献
13.
《地球物理学报》2021,64(7)
区域和全球地震定位越来越多地基于更接近实际的横向非均匀速度模型.速度模型主要来自于地震体波层析成像结果,分辨率不是很高.这样,模型宜于以稀疏网格剖分以减少计算时间和计算机内存需求.当剖分的模型单元较大时,基于射线追踪技术计算复杂介质地震定位中震源轨迹的方法——选取震源轨迹所经过单元节点(位于单元中心,称为震源轨迹节点)为轨迹参考点,利用最小走时树射线追踪方法计算绝对残差场中连接轨迹参考点的射线路径作为震源轨迹——计算结果误差较大,难以满足精定位需要.针对该问题,本文对其进行了改进:不将震源轨迹节点作为轨迹参考点,而是基于插值技术计算每个轨迹节点其法线点对(即该节点与其周围残差正负极性不同的相邻节点组成的点对中梯度绝对值最大的那对)间残差为零的点作为震源轨迹参考点.算例表明:和原方法相比,改进方法计算的震源轨迹更为精细,计算精度提高数十(线性插值)至数百倍(非线性插值),而计算效率基本保持在同一数量级,使利用震源轨迹进行直观、快速和高精度的区域或全球地震事件定位成为可能;壳幔界面反射纵波(PmP)对震源的约束和直达纵波(Pg)相似;同一台站PmP-Pg波到时差约束的震源轨迹对震源深度有很好的约束. 相似文献
14.
弹性波逆时偏移不受倾角和偏移孔径的限制,能够实现任意复杂构造的高精度多波成像,是目前最精确的多分量资料偏移成像方法之一.逆时偏移算法的核心是波场延拓,传统波场延拓以水平基准面为边界条件,基于固定采样步长进行规则网格剖分,采用阶梯近似法处理起伏地表和复杂构造界面时会产生台阶散射,严重影响起伏地表复杂构造的成像精度.基于无网格节点模型,定量分析了弹性波模拟中径向基函数有限差分法的频散关系和稳定性条件.基于此,提出一种基于QR径向基函数的高精度有限差分方法,并提出一种优化的起伏地表自适应节点剖分方法,推导了精确的无网格自由边界条件和弹性波无网格混合吸收边界条件,形成了新的基于无网格的起伏地表弹性波数值模拟方法.此外,本文将此无网格径向基函数有限差分方法应用于精确的纵横波场矢量分解公式,实现了起伏地表弹性波逆时偏移成像.通过对高斯山丘模型,起伏凹陷模型和起伏地表Marmousi-2模型进行数值试算,验证了本文方法的有效性和可行性. 相似文献
15.
提出混合ADI-FDTD亚网格技术开展频散介质GPR正演,即在物性参数变化剧烈局部区域采用细网格剖分ADI-FDTD计算,其他的区域采用粗网格剖分常规FDTD计算,ADI-FDTD突破了CFL条件的限制,可选取与粗网格一致的大时间步长,有效地提高了计算效率.本文首先基于Debye方程,推导了粗网格FDTD及细网格ADI-FDTD频散介质差分格式,着重对粗细两种网格结合的场值交换方式进行了深入探讨,给出了该算法的计算流程.然后以一个薄层模型为例,分别应用粗网格、细网格、混合ADI-FDTD亚网格算法对该模型进行正演,计算资源的占用及模拟精度说明了混合ADI-FDTD亚网格算法的优势.最后,建立频散介质与非频散介质的组合模型,应用3种方法对该模型进行正演,对比3种方法优劣,分析雷达剖面中非频散介质及频散介质中波形特征,有效地指导雷达资料的精确解释. 相似文献
16.
17.
《地球物理学进展》2017,(3)
复杂地表条件下的有限差分地震波场的数值模拟,由于受到低速层和地表起伏的限制,模型速度分布范围变大,一般使用精细的差分网格来抑制频散,提高模拟分辨率,但精细网格会显著增加计算成本.为了能有效地解决这一问题,本文提出一种步长自适应有限差分波动方程数值模拟方法.(1)新方法根据模型中的介质速度分布,对不同的速度区域采用与该速度匹配的空间步长,实现对模型空间网格的步长自适应精细划分.对于速度分布范围大的复杂地表模型,新方法不仅能够极大地减少模型的网格节点数,同时又能提高波场的时间采样步长,减少时间采样数,提高计算效率.(2)推导了不同步长边界网格节点Laplace算子的二阶有限差分表达式,避免了在这些结点进行插值计算产生的假扰动和数值不稳定问题.(3)为了降低有限差分产生的数值频散,本文在常规的差分方程中增加了一频散校正项,能有效地衰减了高波数成分,抑制了数值频散.对复杂近地表的波场数值模拟结果表明,本文提出的步长自适应新方法能够有效减少网格节点数和时间采样数,极大地提高计算效率,计算量比常规粗网格增加一些,但效果能够达到了常规精细网格的模拟结果. 相似文献
18.
利用交错网格有限差分和旋转交错网格有限差分进行各向异性介质弹性波场数值模拟时, 质点振动速度分量与应力张量的网格节点定义方式均不相同, 从而对各向异性波场分离效果产生不同的影响. 针对这一问题, 本文以具有垂直对称轴的横向各向同性(VTI)介质的波场分离为例, 首先分析了两种网格的参数定义方式以及VTI介质波场的分离过程; 其次, 详细研究和分析了这两种网格的参数定义方式对各向异性介质波场分离的影响, 并依据波前面连续性以及波场分离效果等方面, 通过数值模拟实验对该影响进行分析验证. 结果表明, 旋转交错网格的参数定义方式更有利于进行各向异性介质波场数值模拟和波场分离. 相似文献
19.
《中国科学:地球科学》2021,(8)
基于L-S热弹性理论,采用旋转交错网格伪谱法,实现了均匀各向同性介质的一阶速度-应力-温度微分方程组的数值求解和波场模拟.其中,用时间分裂法解决方程组的刚性问题,用旋转交错伪谱算子计算空间一阶导数,用中心伪谱算子计算空间二阶导数;对于热导率变化比较大的双层介质模型和参考温度随深度按梯度分布的模型,用Crank-Nicolson显式方法取代旋转交错网格伪谱法进行计算;讨论了热耦合波场性质和传播规律,对比了常规伪谱法、交错网格伪谱法和旋转交错网格伪谱法热耦合波场模拟效果.数值模拟结果表明:基于L-S热弹性理论,用时间分裂法结合旋转交错网格伪谱法对均匀各向同性介质的波场计算,能够得到稳定的、高精度的模拟结果,但是在热导率变化剧烈的情况下不能用大时间步长进行求解,而且在参考温度不均匀分布的情况下算法不稳定.本文在网格剖分方式与数值算法的优化组合应用方面进行了探索,为将这些方法推广到孔隙热弹性、热黏弹性和各向异性研究奠定了基础. 相似文献
20.
地震波场数值模拟是理论地球物理学和勘探地球物理学的重要研究手段.在众多数值模拟方法中边界元法和有限差分法是两种典型的地震波传播模拟计算方法.边界元法是一种半解析-半数值的边界型方法,它显式地利用边界连续条件,沿着地层边界进行离散,具有降维、高精度和自动满足远场辐射条件的优点;有限差分法是一种典型的基于微分的区域型方法,它隐式地使用边界连续条件,以空间网格形式进行离散和数值逼近,具有高效、实用和容易数值实现的优点.本文以一个半圆形均匀Valley模型和两个非均匀断裂/断层模型为例,从计算精度、计算效率、频散特性以及适用性等方面对这两种方法进行了比较研究.数值计算结果表明:边界元法可以精确地几何描述有内部断点、断面的复杂构造,能够精确地模拟内部不规则界面之间波的反射/传播;有限差分法不能以足够的精度描述几何断点和内部不规则边界.边界元法在高频时计算量大于有限差分法,有限差分法则需要更小的网格间距以压制数值频散.因此,在处理内部非均质和高频计算时,有限差分法更有效;在处理内部不规则边界、断点、大尺度等问题时,边界元法比有限差分法更有优势. 相似文献