共查询到19条相似文献，搜索用时 156 毫秒
1.
2.
笔者给出了一种能够模拟弹性波在任意各向异性介质中传播的二维三分量高阶有限差分算法.相对于常规交错网格有限差分方法,旋转交错网格有限差分方法在介质具有强差异性时能更精确地模拟地震波的传播,避免常规交错网格中因对弹性系数进行插值而带来的误差.采用高阶旋转交错网格有限差分方法模拟并分析了零偏移距横波分裂现象随裂缝介质方位角和倾角变化的响应特征.结果表明:结合完全匹配层(PML)吸收边界条件的高阶旋转交错网格有限差分方法能获得高精度的地震波场模拟数据,并且在边界具有良好的吸收效果;横波分裂现象主要受裂缝走向与波的极化方向之间的夹角影响,受裂缝倾角影响较小,且快慢横波的能量也跟裂缝走向与波极化方向间的夹角有关.具有倾斜对称轴的横向各向同性(TTI)介质倾角的变化可能会导致记录中波到达时的变化,影响快慢横波的时差.利用横波分裂的能量分布和方位各向异性特征,可以帮助检测裂缝的方位角和倾角.横波在多层TTI介质中传播时会发生多次分裂的现象. 相似文献
3.
这里实现了在2.5维各向异性介质中地震波传播的数值模拟。首先从2.5维一阶速度-应力弹性各向异性波动方程出发,得出了该方程的拟谱法数值解法;然后通过对Marmousi速度模型进行数值模拟,表明了拟谱法的有效性;最后计算并分析了均匀和混合各向异性介质中波的传播快照和理论模拟的地震记录,进一步认识了波在各向异性介质中的传播规律。 相似文献
4.
5.
《物探化探计算技术》2015,(2)
通过旋转交错网格高阶有限差分法对一阶速度应力弹性波动方程进行了波场模拟,将耦合的纵、横波场分别投影到纵波、横波的极化向量上纵、横波的分离,然后将分离之后的纵波、横波场按照矢量极化向量的比例再分配,达到了将纵、横波场分解,分解之后矢量波场的和与分离分解之前的矢量波场的和相等。由于分离分解是在旋转交错网格下进行的,因此不需要对波场进行近似或插值处理,提高了波场分解的精度。通过对均匀介质模型和复杂层状介质模型分解前、后波场快照与地震记录的对比,验证了该方法的正确性。波场分离分解技术能够加强对地震波传播规律的认识,具有重要的理论意义。 相似文献
6.
TTI介质qP波数值模拟方法因为考虑了倾角因素,可以比VTI介质qP波数值模拟方法更加准确地描述各向异性介质中地震波场的传播规律。文中用拟声波方程对TTI介质中的地震波场进行了高阶有限差分数值模拟,在改进衰减函数分布方式后,通过坐标变换,利用改进的完全匹配层(perfectly matched layer,PML)边界控制方程对波场边界进行吸收处理,取得了良好的效果;然后分析了拟声波方程数值模拟中的稳定性问题,并对波场中的伪横波进行压制。通过对不同模型的数值模拟,验证了文中使用的TTI介质拟声波波动方程的稳定性以及所采用的PML边界控制方程的可靠性和适用性。 相似文献
7.
井间地震方法观测到的波场信息丰富且复杂,其难点是如何识别和分离这些波场信息。这里从纵波波场为无旋场,横波波场为无散场出发,利用一阶速度～应力方程的交错网格高阶有限差分法,对井间弹性波纵、横波场分离问题进行了数值模拟。经数值模拟试验表明,同二阶的位移方程相比,使用等价的波场分离交错网格的一阶速度应力方程,无论是模拟精度还是模拟的稳定性方面都有了较大的提高,同时该横波资料的获取方法为井间地震资料的解释应用,提供了一种新的技术思路。 相似文献
8.
9.
完全3维弹性波数值模拟计算时间长,并且占用庞大的计算资源,这不利于在计算机配置不高的情况下进行科学研究,而二维弹性波数值模拟又达不到三维模拟的精度;同时,当模型、波场空间分布比较复杂时,传统的3维波动方程拟谱法模拟结果比较差.因此,在较高数值精度的一阶应力-速度弹性波动方程的基础上,采用傅氏变换仅计算y方向的偏导数,利用有限差分方法计算x、z方向和时间的偏导数,即利用2.5维数值模拟方法,实现在二维介质中计算三维弹性波场.最后通过数值模拟实现了在各向异性介质中多波波场的数值模拟,验证了2.5维方法是一种高精度、高效率、且能适应复杂模型的正演模拟方法,通过波场分析进一步认识了波在各向异性介质中的传播规律. 相似文献
10.
《物探与化探》2016,(1)
相比于VTI介质q P波数值模拟方法,考虑倾角因素的TTI介质q P波数值模拟方法能够更加准确地描述各向异性介质中波场的传播特征。常规的声学近似方法往往会造成TTI介质中倾角急剧变化的区域出现数值不稳定,笔者首先引入一个各向异性控制参数,推导了稳定形式的TTI介质二阶q P波方程,保证q P波场的稳定传播;其次,通过引入波场的伪速度分量,推导了其等价的一阶应力—速度形式及相应的PML边界条件,并应用优化的旋转交错网格有限差分(RSGFD)方法实现了精确的数值模拟。数值结果表明:TTI介质一阶q P波方程能够稳定、有效地模拟q P波的运动学特征,利用优化的RSGFD方法可以得到精确的合成地震记录,同时可以相对地提高计算效率。 相似文献
11.
基于一阶速度-应力波动方程,采用高阶交错网格有限差分数值模拟方法,对弹性及黏弹性TTI介质进行正演数值模拟。模拟时采用完全匹配层吸收边界条件(PML)消除边界反射。同时设计了层状介质模型、断层模型,通过模型的正演计算,得到了不同时刻的地震波波场快照及合成地震记录,分析其波场运动学及动力学特征。模拟结果表明,交错网格有限差分法可以很好地完成对复杂介质的波场模拟,具有较高的精度和可靠性。 相似文献
12.
13.
14.
研究井间地震波场的形成过程以及波场的传播机理、规律,对于指导实际井间地震勘探有着重要的意义.基于具有垂直对称轴的横向各向同性(VTI)介质中的一阶准P波方程,应用具有无条件稳定性质的紧致交错网格隐式差分格式求解该方程.重点研究了紧致交错网格求解该方程的完全匹配层(PML)吸收边界条件,在此基础上实现了VTI介质中一阶准P波方程的井间地震波场的正演模拟.数值算例表明:紧致交错网格能精准模拟VTI介质中准P波的传播过程,得到高精度的正演结果.一阶准P波方程能以足够的精度描述VTI介质中准P波特征.完全匹配层吸收边界能有效地解决人工边界问题,是一种高效的边界吸收算法. 相似文献
15.
非均匀介质中交错网格高阶有限差分数值模拟 总被引：5，自引：0，他引：5
地震波场的数值模拟一直是地球物理学的一个重要的研究领域,而在数值正演模拟方法的研究中,计算精度和计算效率是评价该方法有效性及优越性的二个关键问题。这里从一阶速度—应力弹性波动方程出发,着重介绍如何构造离散化模型的网格,如何求解空间导数,如何选取边界条件等内容,从而更有效地提高数值计算的精度与计算效率。文中构造了不同类型的介质模型,并在交错网格中,利用高阶有限差分模拟非均匀介质的波场传播。模拟结果表明,该方法实现简单,具有很好地稳定性和较高的精度,能够直观、高效地反映出介质中波场的传播规律。 相似文献
16.
17.
岩石介质具有黏滞性,黏弹性介质模型相对于弹性介质模型更为接近岩石介质的真实情况。根据Kelvin-Voigt黏弹性单元体模型本构关系,推导了柱坐标系下各向同性黏弹性介质时间域交错网格有限差分方程,进行了黏弹性介质井孔声场的数值模拟。模拟结果表明：品质因子和声源中心频率对黏弹性介质中的井内外声场均有影响：井外声场和井内声场的衰减均随品质因子的增大而减小,井外声波和井内反射波振幅随之增大;井外声场的衰减随声源中心频率的增大而减小,井外声波振幅随之增大,井内声场的衰减随声源中心频率的增大而增大,井内反射波振幅随之减小。 相似文献
18.
贴体网格有限差分正演模拟算法不仅能够精确模拟任意起伏地形下的波场特征,且计算效率较高,是一种很有应用前景的处理西部复杂地表问题的方法;然而,目前求解波动方程时常用的同位网格和标准交错网格,在处理贴体网格起伏地表正演模拟时存在诸多问题。为此,将全交错网格引入到曲线坐标系下,避免了标准交错网格的插值误差和同位网格中奇偶失联引起的高频振荡现象,提高了模拟精度,减小了算法实现的复杂度。在自由边界条件实施时,采用牵引力镜像法计算速度分量,速度自由边界条件配合紧致交错差分格式更新应力分量,得到了较好的效果。随后,重点研究了贴体全交错网格正演模拟算法的影响因素,考虑了网格正交性、网格间距和网格拼接等的影响,并取得了如下认识：算法对网格的正交性没有过分要求;网格间距的突变会引起虚假反射的产生;不同类型的网格拼接对模拟结果不会造成明显的影响。 相似文献
19.
利用时间分裂的错格伪谱法模拟地震波在基于改进BISQ模型的双相介质中的传播 总被引：1，自引：0，他引：1
改进的BISQ(Biot-Squirt)模型中各参数具有明确的物理意义和可实现性,在不引入特征喷流长度的情况下可将Biot流动和喷射流动两种力学机制有机地结合起来;而高精度的地震波场数值模拟技术是研究双相介质地震波传播规律的重要手段。本文从本构方程、动力学方程和动力学达西定律出发,推导了基于改进BISQ模型的双相各向同性介质的一阶速度--应力方程组;采用时间分裂错格伪谱法求该方程组的数值解,模拟半空间及层状双相介质中的地震波场。数值模拟结果表明:①与传统方法相比,时间分裂错格伪谱法波场数值模拟的精度更高,压制网格频散效果更好;②在非黏滞相界情况下,慢纵波呈传播性,而在黏滞相界情况下,慢纵波呈扩散性,以静态模式出现在震源位置;③双相介质分界面处,各类波型复杂的反射透射规律可由数值模拟结果清晰展现。 相似文献