共查询到16条相似文献，搜索用时 202 毫秒
1.
基于精细积分法的三维弹性波数值模拟(英文) 总被引：1，自引：0，他引：1
波动方程有限差分法是地震数值模拟中的一种重要的方法,对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想,产生了不稳定性。精细积分法在有限差分法的基础上,在时间域采用解析解的表达形式,在空间域保留任意差分格式,发展成为半解析的数值方法。本文结合并发展了以往学者的成果,推导了任意精细积分法的三维弹性波正演模拟计算公式,并对其稳定性进行了数值分析。在计算实例中,实现了精细积分法二维和三维弹性波模型的地震正演模拟,对计算结果的分析表明,精细积分法反射信号走时准确,稳定性好,弹性波场相较于声波波场,弹性波波场成分更为丰富,包含了更多波型成分(PP-和PS-反射波、透射波和绕射波),这对实际地震资料的解释和储层分析有重要的意义。实践证明,该方法可直接应用到弹性波的地质模型的数值模拟中。 相似文献
2.
《地球物理学进展》2016,(2)
地震波场数值模拟不仅是研究复杂地区地震资料采集、处理和解释的有效辅助手段,而且是研究地球深部精细构造和地球深部探测的有效工具.射线理论和波动方程理论是地震波场数值模拟的理论基础.射线理论主要刻画地震波在介质中传播的走时场、地震射线等运动学属性;波动方程理论通过求解波动方程来描述地震波在介质中传播的弹性动力学响应(能量衰减、相位特征、偏振属性、以及全波形等).基于波动方程理论的波场数值模拟由于能够引入丰富的波场信息,使得人们对不同介质中地震波的传播过程有了较全面的了解.本文以二维层状均匀介质模型为例,通过射线追踪法和交错网格有限差分法模拟得到的波场快照图、单炮地震记录剖面、合成理论地震图的分析比较,不但对地震波在各向同性层状均匀介质中的传播规律和特点有了深刻的理解和认识,同时又可以相互验证两种不同方法的正确性和有效性. 相似文献
3.
4.
求解波动方程的任意差分精细积分法 总被引：5，自引：0，他引：5
在地震偏移成像技术中，常常要反向外推波场，因此涉及到波动方程的数值求解问题。本文提出了一种求解偏微分方程的新的半解析方法——任意差分精细积分(ADPI)法。其大体思路是：空间域上作坐标离散，但不采用古典差分法的等分离散方式，而是一种相对自由的、任意的离散法；时间域上则采用子域精细积分的方法，既保留了精细积分法的高精度，又克服了工作量大、占用存储大等缺点。该方法具有精度高、带宽小、稳定性好等多项突出优点，并且可以灵括处理各类边界条件。本文从简单一维、较普遍一维、以及二维3种情形讨论该方法对波动方程的具体应用，结合实际的算例，详细分析了各种算法的可行性和精度特点。 相似文献
5.
6.
单程波近似实际上是一种多次前向散射和单次后向散射近似.利用单程波近似来描述波传播可以极大地节省地震数值模拟的计算时间和内存,实现地震波长距离传播模拟和三维地震模拟快速计算.本文基于单程波近似和波动积分方程的分离变量逼近,从广义Lippmann-Schwinger波动积分方程推导出耦合反射/透射系数的单程波传播算子.该算子由两部分构成：分离变量Fourier单程波传播算子和薄板间的反射/透射系数表达.前者将常规的Fourier分裂步单程波传播算子(SSF)推广适应横向强速度变化介质和大角度传播波场.后者是利用垂直波数来表示反射/透射系数,自然耦合到波场传播的计算过程中,其为地质界面倾角的隐式表达,精确描述振幅随入射角的变化,能适应任意复杂的模型.通过两个数值算例和一个实际地质模型的计算,本文将该方法和边界元法进行了比较,结果表明：在算例给出的介质横向速度变化情况下,本文提出的方法在相位和振幅方面与全波数值方法基本吻合. 相似文献
7.
本文针对声-弹耦合介质,为尽可能的减少频率域正演模拟的计算内存,提高计算效率,在一阶非均质位移-应力波动方程的基础上,借助等效交错网格思想并充分考虑密度参数空间变化对地震波传播的影响,推导了声-弹耦合地震波波动方程.在流相介质和固相介质中分别采用非均质情况频率域二阶声压标量波、二阶纯位移控制方程,为保证流、固相介质间地震波能量的稳定传输和有效交换,提出了声-弹耦合界面转换过渡层方法,并详细阐述了过渡层与上下介质空间差分具体耦合方法.在与非均质纯位移波动方程正演结果对比分析的基础上,首先采用各向同性单层流相介质模型进行正演模拟验证了声-弹耦合方程数值模拟中过渡层策略的有效性和准确性,随后又数值模拟了地震波在声-弹耦合介质简单模型和复杂Marmousi2模型中的传播,验证了本文方法稳定性和准确性,同时该方法可以简单的推广到三维情况. 相似文献
8.
9.
10.
吸收衰减是地震波在实际地球介质中传播的固有特征.在实际应用中,通常假设表征吸收衰减特征的品质因子Q在地震频带范围内不随频率变化.高阶广义流变模型能够在时间域内精确的表征品质因子Q不随频率变化的特征,为黏弹性介质波动方程精细模拟奠定了基础.基于广义标准线性体模型理论,采用最小二乘拟合方法对Q值不随频率变化特征进行拟合,分析了不同阶次广义标准线性体模型对黏弹性介质Q值特征的拟合程度,在权衡计算精度和三维计算量的基础上,确定了五阶广义标准线性体模型并建立了相应的三维黏弹性波的速度-应力方程,结合CFS-PML边界条件开展了高精度三维黏弹性波正演模拟.通过均匀介质正演模拟,验证算法的正确性,明确了地震波的传播时的吸收衰减特征,对三维盐丘模型进行数值模拟,表明了五阶广义标准线性体可以精确的模拟黏弹性介质地震波吸收衰减特征. 相似文献
11.
Jia Xiaofeng 《中国地震研究》2004,18(1):36-41
Wave equation migration is often applied to solve seismic imaging problems. Usually, the finite difference method is used to obtain the numerical solution of the wave equation. In this paper, the arbitrary difference precise integration (ADPI) method is discussed and applied in seismic migration. The ADPI method has its own distinctive idea. When dispersing coordinates in the space domain, it employs a relatively unrestrained form instead of the one used by the conventional finite difference method. Moreover, in the time domain it adopts the sub-domain precise integration method. As a result, it not only takes the merits of high precision and narrow bandwidth, but also can process various boundary conditions and describe the feature of an inhomogeneous medium better. Numerical results show the benefit of the presented algorithm using the ADPI method. 相似文献
12.
13.
14.
波动方程数值模拟方法是研究地震波场传播的一种重要手段,本文采用交错网格高阶有限差分方法分别对双程声波方程和双程弹性波方程进行了波场数值模拟,并且根据定位原理采用傅立叶有限差分算子进行了单程波方程数值模拟,在分析定位原理的基础上,对其计算过程稍作修改,将延拓到地面的波场直接由每个检波点接收,无需横向叠加过程,得到了单程声波方程共炮记录.基于不同波动方程的数值模拟结果表明,双程波方程结果包含直达波、多次波等干扰波,信噪比低;单程波数值模拟结果只包含了介质分界面的一次反射波,信噪比高,但对于大角度入射波误差较大,并且对于同一个地质模型而言,双程弹性波方程计算速度最慢,双程声波方程次之,单程声波方程计算速度最快.因此对于复杂地质模型,三种模拟方法可以取长补短,综合应用. 相似文献
15.
正演模拟技术在地震采集设计中的应用 总被引：1，自引：0，他引：1
随着地震勘探开发的不断深入和发展,地震勘探的主战场逐渐向复杂地区转移.复杂地区既指地表条件复杂的地区,也指地下地质构造和地层条件复杂的地区,这些都对地震勘探提出了新的挑战和更高的要求.地震正演模拟正是开展此类问题研究的一个重要手段和方法.目前市场上具有正演模拟功能的软件大多是根据射线理论采用射线追踪的方法来完成正演模拟的,这种方法不能很好地反映地震波的动力学特征,特别是在复杂地区难以得到正确的结果.本文利用高阶有限差分有效克服了常规有限差分算法求解波动方程的频散问题,并以高效的OpenMP并行计算模式进行了并行优化,较大程度上提高了正演计算的速度和精度;同时实现了二阶Higdon边界条件,改善了边界吸收效果;也在一定程度上提高了计算速度.塔中地区主要目的层埋藏深,逆断层发育,地震反射特征复杂,增加了地震勘探的难度.本文依据该地区的地质模型,利用波动方程正演技术论证了该地区的地震采集观测系统,为该区地震采集观测系统的设计提供了科学依据. 相似文献
16.
模拟地震波传播的三维逐元并行谱元法 总被引：2，自引：0，他引：2
高效地震波场正演模拟对于复杂模型中地震波传播与成像研究至关重要.本文在谱元法原理框架内,对已有逐元谱元法改进,提出一种新的逐元并行谱元法求解三维地震波运动方程,并得到地震波场.逐元并行谱元法的核心思想在于在单元上进行质量矩阵与解向量的乘积运算,并将此运算平均分配至每一个CPU计算核心,此处理有利提升谱元法的并行计算效率.同时,根据Gauss-Lobatto-Legendre(GLL)数值积分点与插值点重合的特点,将稠密单元刚度矩阵的存储转化成单元雅克比矩阵行列式的值及其逆的存储,大幅减少谱元法计算内存开销.此外,在模型边界上利用逐元并行谱元法求解二阶位移形式完美匹配层(PML)吸收边界条件,消除边界截断而引入的虚假反射.通过逐元并行谱元法得到的数值解与解析解对比,以及实际地震波场模拟,数值结果证实了逐元并行谱元法用于地震波场模拟的高效性. 相似文献