排序方式： 共有50条查询结果，搜索用时 15 毫秒
1.
数值求解波动方程是大尺度正演波场模拟、基于波动方程的地震偏移和反演成像的关键.本文针对求解二维声波方程的Runge-Kutta 间断有限元（RKDG）方法的数值频散问题,从理论推导和数值分析的角度进行了深入研究,并将其与近似解析离散化方法（Optimal Nearly Analytic Discrete Method,简称ONAD 方法）、Lax-Wendroff 修正方法、交错网格（Staggered-Grid,简称SG）方法的数值频散进行了比较研究.结果表明：RKDG方法以及近似解析离散化方法在压制数值频散方面要好于上述其他方法,特别是空间精度为3阶的RKDG方法,即使当空间步长取波长的一半,即一个波长内取2个网格点时,最大的频散误差也不超过1.67%.同时,我们也通过波场模拟对比研究了不同数值方法的数值频散问题,进一步直观地验证了数值频散的理论分析结果. 相似文献
2.
本文定义了各向异性黏弹性参数修正因子,并将其引入到黏弹性模型中以体现泥质含量对黏弹性机制的影响,同时将波传播过程中孔隙介质骨架黏弹性力学机制与两种孔隙流体流动力学机制（Biot流动和喷射流动机制）有机地统一起来处理,从而给出了描述含泥质低孔渗孔隙各向异性介质中波传播规律的黏弹性Biot/squirt (BISQ)模型.数值计算结果表明,入射波的方位角、各向异性渗透率以及泥质含量等对含流体复杂孔隙介质中波频散和衰减的影响具有显著的方位各向异性特征,在低频范围内（地震波勘探频率）黏弹性力学机制对波传播能量的衰减起主导作用. 相似文献
3.
To simultaneously take into account the Biot-flow mechanism, the squirt-flow mechanism, and the frame-viscoelasticity mechanism, a generalized viscoelastic BISQ （Biot/squirt） model is developed for wave propagation in clay-bearing sandstones based on Dvorkin＇s elastic BISQ model. The present model is extended to a wide range of permeabilities （k 〉 0.05 mD） by introducing a dimensionless correction factor for viscoelastic parameters, defined as a function of the permeability and the clay content. We describe the frame＇s stress-strain relationship of the clay-bearing sandstones by the differential constitutive equations of generalized viscoelasticity and then derive the viscoelastic-wave dynamic equations. With the assumption of a plane-wave solution, we finally yield the phase velocities and the attenuation coefficients by solving the dynamic wave equations in the frequency and wave number domain. The comparison of numerical results and experimental data shows that the generalized viscoelastic BISQ model is applicable for modeling the wave propagation in most of the sandstones mainly bearing kaolinite clay. 相似文献
4.
地震层析成像是研究地球内部结构的主要方法.文章利用首都圈区域内3231个地震挑选出来的29839条直达P波走时和29972条直达S波走时数据,首次使用拟波形层析成像方法对首都圈地区地壳的P波和S波速度,以及泊松比进行反演成像,获得了该区域三维高分辨率VP、VS以及泊松比结构.成像结果表明,首都圈地壳结构非均匀性强烈,地表地质构造特征与浅层的速度异常具有明显的一致性.中地壳康拉德面深度附近广泛存在的相对低速异常层可能是多期构造运动叠加的结果.该区域地震多发生于高低速及泊松比异常交界区域.唐山地震和历史上的三河-平谷地震震源下方都存在明显的相对低速异常,可能与流体的存在有关.深部速度结构异常的一致性可能暗示唐山地震和三河-平谷地震的深部发震构造具有一定的相似性.高分辨率精细速度结构的获得,为深入认识首都圈地区地壳结构的精细非均匀性和动力学机制提供了重要的科学依据. 相似文献
5.
ＡｎａｌｔｅｒｎａｔｅａｌｇｏｒｉｔｈｍｏｆｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｉｎａｎｉｓｏｔｒｏｐｉｃｍｅｄｉａＤＩＮＧ－ＨＵＩＹＡＮＧ（杨顶辉），ＪＩ－ＷＥＮＴＥＮＧ（滕吉文）ａｎｄＺＨＯＮＧ－ＪＩＥＺＨＡＮＧ（张中杰）Ｉｎｓｔｉｔｕｔｅｏ... 相似文献
6.
模拟地震波传播的三维逐元并行谱元法 总被引：2，自引：0，他引：2
高效地震波场正演模拟对于复杂模型中地震波传播与成像研究至关重要.本文在谱元法原理框架内,对已有逐元谱元法改进,提出一种新的逐元并行谱元法求解三维地震波运动方程,并得到地震波场.逐元并行谱元法的核心思想在于在单元上进行质量矩阵与解向量的乘积运算,并将此运算平均分配至每一个CPU计算核心,此处理有利提升谱元法的并行计算效率.同时,根据Gauss-Lobatto-Legendre(GLL)数值积分点与插值点重合的特点,将稠密单元刚度矩阵的存储转化成单元雅克比矩阵行列式的值及其逆的存储,大幅减少谱元法计算内存开销.此外,在模型边界上利用逐元并行谱元法求解二阶位移形式完美匹配层(PML)吸收边界条件,消除边界截断而引入的虚假反射.通过逐元并行谱元法得到的数值解与解析解对比,以及实际地震波场模拟,数值结果证实了逐元并行谱元法用于地震波场模拟的高效性. 相似文献
7.
间断有限元（Discontinuous Galerkin：DG）方法具有低数值频散、网格剖分灵活、能模拟地震波在复杂介质中传播等优点.因此，本文将一种新的DG方法推广到双相和黏弹性等复杂介质的地震波场模拟，发展了求解Biot弹性波方程和D'Alembert介质波动方程的DG方法.首先通过引入辅助变量将Biot双相介质弹性波方程和D'Alembert介质波动方程转化为关于时间-空间的一阶偏微分方程组，然后对该方程组进行DG空间离散，得到半离散化的常微分方程组.最后，对此常微分方程组，应用加权的Runge-Kutta格式进行时间推进计算.数值结果表明，DG方法可以有效地求解Biot双相介质弹性波方程和D'Alembert介质波动方程，并能很好地压制因离散求解波动方程而产生的数值频散，获得清晰的各种地震波震相. 相似文献
8.
9.
谱元法已成为区域性乃至大陆性尺度地震波场模拟的重要工具.对于区域或大陆尺度层析成像而言,地球曲率不可忽略,此时模拟地震波传播采用球坐标系更为合适.本文从球坐标系下弹性波动方程弱形式出发,基于球坐标系变分原理给出了球坐标系下求解三维地震波方程的谱元法.另一方面,计算Fréchet敏感核是进行全波形反演的关键,本文借助伴随原理,推导了全波走时层析成像三维Fréchet敏感核表达式.为了验证球坐标系下谱元法的精度,我们将数值模拟结果与normal mode方法得到的解析解在1-D PREM模型下进行了对比.同时,我们将此方法应用到华北克拉通区域,以期获得地球内部结构精确成像.基于3-D全球径向各向异性地幔模型S362ANI和3-D地壳模型Crust1.0,我们建立了华北克拉通初始3-D背景模型,并将数值模拟结果与实际观测台站记录波形资料进行对比分析,利用互相关方法提取走时残差,最后给出了Fréchet敏感核在3-D空间中的分布,这些工作为下一步开展球坐标系下三维大尺度全波形反演奠定了基础. 相似文献
10.
本文基于声波方程的哈密尔顿系统,构造了一种新的保辛数值格式,简称NSPRK方法.该方法在时间上采用二阶辛可分Runge-Kutta方法,空间上采用近似解析离散算子进行离散逼近.针对本文发展的新方法,我们给出了NSPRK方法在一维和二维情况下的稳定性条件、一维数值频散关系以及二维数值误差,并在计算效率方面与传统辛格式和四阶LWC方法进行了比较.最后,我们将本文方法应用于声波在三层各向同性介质和异常体模型中的波传播数值模拟.数值结果表明,本文发展的NSPRK方法能有效压制粗网格或具有强间断情况下数值方法所存在的数值频散,从而极大地提高了计算效率,节省了计算机内存. 相似文献