共查询到18条相似文献，搜索用时 914 毫秒
1.
TTI介质qP波数值模拟方法因为考虑了倾角因素,可以比VTI介质qP波数值模拟方法更加准确地描述各向异性介质中地震波场的传播规律。文中用拟声波方程对TTI介质中的地震波场进行了高阶有限差分数值模拟,在改进衰减函数分布方式后,通过坐标变换,利用改进的完全匹配层(perfectly matched layer,PML)边界控制方程对波场边界进行吸收处理,取得了良好的效果;然后分析了拟声波方程数值模拟中的稳定性问题,并对波场中的伪横波进行压制。通过对不同模型的数值模拟,验证了文中使用的TTI介质拟声波波动方程的稳定性以及所采用的PML边界控制方程的可靠性和适用性。 相似文献
2.
常规变加密网格有限差分波动方程数值模拟方法采用水平分层加密网格,该网格剖分策略在适应地形起伏和近地表速度结构变化特征方面效果较差。针对该问题,提出一种起伏多重变加密网格有限差分波动方程数值模拟方法。该方法根据地形起伏和近地表低速层到高速层的速度分布特征进行网格加密;采用不同网格中的变系数差分格式离散声波方程,在保证波场模拟精度的同时兼顾计算效率;同时,为了进一步保障地表附近波场模拟的精度,地表附近最细网格中的差分格式不做降阶处理,针对位于地表以上的虚像点的波场值,提出一种融入自由地表边界条件的法向虚像外推法。算例分析验证该算法对速度模型不同区域进行的网格多次加密显著提高了计算效率,以黄土塬实际模型为例,耗时为常规1 m×1 m网格耗时的43.3%,并可达到和细网格基本一致的模拟精度,模拟误差控制在10-12范围内,同时表现出很好的近地表散射压制和边界吸收效果,且算法能稳定地适应实际复杂地表介质。 相似文献
3.
4.
通过对声波方程采用二阶有限差分与四阶有限差分近似,进行数值模拟及波场分析,边界条件利用吸收边界和增加衰减带混合使用的方法。计算结果表明,四阶有限差分的精度明显高于二阶有限差分,并且边界吸收效果好。对于单炮记录难于模拟复杂模型的情况下,实现了单次叠加剖面和爆炸反射面的零炮检距剖面,并进行了偏移处理。地震剖面更加清晰地反映了复杂构造的信息,更好地实现了对复杂地质模型的正演模拟。 相似文献
5.
常规地震勘探理论无法完全解释复杂多变的地质特征,较难达到高精度地震资料处理的目的,因此必须发展适应复杂非均匀介质的地震波场传播理论,以更好地刻画真实的复杂构造、岩性、裂隙油气藏的特点。根据实际火山岩样的随机特征,构建了相应的多尺度组合型随机孔洞介质模型,并进行了高精度地震波正演数值模拟和逆时成像研究,制作了在不同频率条件下的自激自收剖面和逆时偏移剖面,计算结果表明,构建的多尺度随机孔洞介质模型可以有效模拟真实的地质变化特征,同时指出准确区分由地层随机非均匀性引起的散射波场响应和干扰信号是实现高分辨率地震勘探的前提。 相似文献
6.
地震波斜入射情形下水平成层半空间自由场的时域计算 总被引：8，自引：1，他引：8
在应用有限元模拟地震工程中的波动散射问题时，为了计算边界节点的运动和实现波动能量的输入，必须首先确定边界节点的内行波场即背景水平成层半空间的自由场。利用显式有限差分法的内节点位移计算公式和基于水平成层介质波动传播特点建立的相邻节点间自由场运动的关系式，给出了计算水平成层半空间自由场的时城方程(以入射侧边界节点为例)。本方法是显式有限元意义上的数值精确解，可应用于近场波动数值模拟中内行波场与散射波场的时城同步模拟。 相似文献
7.
深水区地震波传播特性研究 总被引：1，自引：0，他引：1
根据层速度与Q值的经验公式和对射线参数进行泰勒展开的方法，以及二维有限差分声波波动方程波场数值模拟技术，为海上长排列大容量震源地震资料采集及高分辨率地震勘探野外施工的观测系统设计及数据处理方案等提供参考意见。 相似文献
8.
起伏地表条件下的声波散射数值模拟的积分方程法 总被引：2，自引：1，他引：2
从散射理论的角度来看,起伏地表可以看作是一种特殊的扰动介质,因此应用散射积分方程求解起伏地表条件下的散射场在理论上是可行的。从三维频率域声波方程出发,由格林函数定理,得出起伏地表条件下的散射积分方程。散射积分方程为关于起伏地表的面积分和与速度扰动体有关的体积分之和,同时给出了格林函数在奇异点的积分方法。由于数值离散求解积分方程存在着计算时间太长和存储内存不足的问题,采用电磁散射积分方程的拟解析近似的方法。在假设反射函数为缓变函数的基础上,最终得到其近似表达式,因此散射场的数值求解不必再借助于代数方程组,只要进行数值积分即可。这种方法避开了传统数值计算方法存在的问题,为地震散射波场快速正演模拟打下了基础。理论分析表明,这种方法适用于小扰动的问题。当扰动较大时,拟解析近似会产生较大的误差。 相似文献
9.
基于波动方程的地震波正演模拟能够准确揭示波在介质中传播的振幅、频率与相位变化,真实地反映出波的动力学特征。有限差分法因其精度高、计算效率高而成为波场数值模拟的一个重要方法。根据地震波的近似方程———声波方程,采用二阶有限差分网格数值算法,利用透明边界条件,设计典型地质模型,进行地面波场的地震记录模拟。通过稳定性条件分析、震源子波优选以及频散抑制,成功合成了倾斜地层模型与正断层模型的地面共炮点波场记录,直观地显示出直达波、反射波、绕射波等丰富的波场信息,获取各模型的波场快照。结果表明,有限差分法结合透明边界条件对于获取介质丰富的波场信息具有很大的优势。 相似文献
10.
波动方程数值模拟是研究地震波传播机理的重要工具,有限差分求解波动方程是当前地震波数值模拟的主要方法之一。当地下介质中的地震波速度较低或地震波高频成分丰富时,常规有限差分技术常常产生严重的数值频散误差,这种误差会降低数值模拟的精度,影响对地震波传播机理的分析。为压制地震波数值模拟时产生的数值频散误差,提高波场模拟精度,提出了基于NAD算子的时间四阶精度波动方程差分格式。根据对应的差分格式,分析了该差分格式的数值频散关系。与常规四阶精度差分算法的频散曲线相比,基于NAD时间四阶精度差分方法不但能够实现时间频散的有效压制,同时其基于更多网格点的位移分量和位移梯度分量空间微分求解方法还能够实现空间频散的有效压制。另外在相同模型条件下,基于NAD算法的声波方程时间四阶差分解法可采用大网格对模拟空间进行差分离散,减少网格数,提高计算效率。 相似文献
11.
等时叠加波动方程叠前正演 总被引：6，自引：3，他引：6
非零炮检距地震正演模拟在复杂区叠前地震资料处理和解释中有着重要意义。通常使用射线追踪法和全程波动方程法来实现这一模拟。通过提出等时叠加原理和数学检波器表达式，形成了非零炮检距地震正演的新一类方法——单程波动方程法。该方法与射线追踪法相比能保持波的动力学特征，正演结果与实际地震记录相似性好。并且在计算过程中，人工干预的程度低，操作简便。与全程波动方程法相比，该方法计算效率高，且不产生直达波、多次波等干扰，信噪比高，与目的层相对应的一次波突出，波型识别难度不大。因此，单程波动方程法特别适合于复杂油气藏反射地震资料的正、反演研究。 相似文献
12.
利用定位原理进行三维叠前正演,仅靠单程声波方程的延拓即可形成复杂地质条件下的正演记录,运算效率高无多次反射出现。但速度模型中对于过大或过小的速度常使波场在延拓时有频率和波数成分的损失,这种损失引起波场畸变并严重影响正演记录的质量。本文采用波场替换技术,在每次延拓时用零相位子波形成一个同时保留了原波场的动力学特征与空间位置信息的新波场,然后用新波场取代原畸变波场并参与正演记录合成。在频率域单程方程三维叠前正演中特别适合做波场替换,因为震源脉冲与接收响应(均可视为点脉冲)的延拓波场在形态上相似,而且在延拓波场中每道记录之振幅包络只有一个极值点。理论试算结果表明,采用波场替换技术后最终正演记录质量比替换前有了很大改善。 相似文献
13.
利用Liner的2.5-D时间域声波方程和高阶有限差分近似,研究了2.5维声波方程在时间-空间域中的衰减特性。简单地质模型(水平二层介质)的数值计算结果证实了下列事实:在Liner的声波方程中的最后一项,即含有1/t2的一项对整个计算结果没有影响,因此可以忽略。这一结论与文献中已经发表的数值结果相吻合。同样的结论对于由Williamson和Pratt提出的2.5-D声波方程也适用。 相似文献
14.
15.
南极和北极地区在全球变化中的作用研究 总被引：15，自引：0，他引：15
多次波问题是海洋地震勘探中最突出的问题之一 ,如何有效地压制多次波是解决问题的关键。多次波数据处理的压制方法和技术分为运动学的滤波方法和动力学的波动方程方法两大类。波动方程方法由于几乎不需要先验信息和具有较好的压制效果 ,已经成为多次波压制方法和技术发展的主要趋势。文章首先对波动方程方法进行了较全面的综述 ;然后 ,对波动方程方法中的反馈迭代反演方法做了深入的讨论 ;最后 ,通过理论模型的正反演计算 ,验证了运用波动方程反馈迭代反演法不需要知道震源函数、地下岩层结构和岩性等先验信息。此方法适用于压制多种类型的多次波 ,具有优异的振幅保真性等独特的优点。 相似文献
16.
二维声波方程稳健迭代速度反演 总被引：1，自引：0，他引：1
从二维声波方程初、边值问题出发,通过引进一个广义幂指数误差分布函数,以及对波动方程的Ｌｉｐｐｍａｎ－Ｓｃｈｗｉｎｇｅｒ方程的解进行Ｂｏｎ近似,导出一个残差加权迭代最小二乘（稳健迭代）算法,进而实现对地下介质剖面的速度结构以演计算。模拟计算表明,本方法具有较高的精度和较强的抗干扰能力,是求解声波方程反演问题的一种有效方法,同时也为地下介质的速度分布成像提供一种新的技术。 相似文献
17.
零偏移距地震记录是叠后偏移的基础,常用射线追踪法来制作,但当介质复杂时,射线会出现交叉现象,对三维情况计算量更大。从二维标量波动方程出发,用伪谱法研究了零偏移距地震记录的合成,给出了不同情况下的数值计算公式,并对边界吸收问题作了有效处理。数值计算结果表明了方法的正确性和有效性。该方法可用于偏移中正演模拟的快速计算,且易推广到三维情形。 相似文献
18.
分离P-P波和P-SV波最常用的是τ-p变换法, 但是在τ-p域中, P-P波和P-SV波常有重叠部分, 不能很好分离。地震勘探中, 当炮检距小于界面深度时, P-P波和P-SV波的时距曲线方程都可近似为抛物线方程。针对这种情况, 提出了一种基于抛物线方程的坐标拉伸τ-p变换的P-P、P-SV 波分离方法, 该方法首先对地震记录沿空间方向上进行坐标拉伸, 再进行τ-p正变换, 在τ-p域中分离P-P波和P-SV波, 其次对P-P波和P-SV波分别进行τ-p反变换, 实现P-P波和P-SV波的分离。与常规τ-p变换法对比分析表明, 利用该方法能很好地将合成地震记录中的P-P波和P-SV波分离。 相似文献