首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
研究了各向异性介质中弹性波动方程交错网格高阶差分法求解的数值相速度色散和数值相速度误差。从而表明了交错网格差分法一般情况下要求每个最短波长离散网格点数约为6。同时表明了高阶差分(4,6)阶的数值各向异性要比低阶的(2,2)阶小得多。  相似文献   

2.
对气煤煤样孔隙率测试的基础上,讨论了各向异性系数与孔隙率的关系,表明各向异性的大小与孔隙率之间有一定的关系。应用高阶交错网格差分算法,对一阶速度-应力弹性波方程进行模拟,分析了其稳定性和边界条件。通过对Thomsen各向异性参数的模拟表明,主要控制的是P波的各向异性强弱,δ主要控制水平垂直正交方向上的各向异性强弱。对煤田各向异性研究具有重要意义。  相似文献   

3.
王宁  吕希华  闫磊  苗长盛 《世界地质》2010,29(1):130-137
采用交错网格技术将速度-应力方程中的速度对时间的导数转化为应力对空间的导数, 将弹性波动方程表示为与二阶双曲方程等价的一阶应力-速度公式, 以实现三维三分量地震波场模拟。对获得的波场快照、VSP记录图像和地面记录图像进行波场分析, 发现在波场图像中存在明显的拟P波、拟快横波和慢横波, 还出现了横波分裂、横波分裂盲区及波面三分叉等地震波在方位各向异性介质传播时产生的特殊现象。研究结果表明, 用交错网格三维模拟方法研究方位各向异性介质对方位各向异性介质中地震波传播反演能起到很强的辅助作用。  相似文献   

4.
在地震勘探中,粘弹性各向异性介质相比弹性介质更能实际表征地下介质的性质,这里主要针对Kelvin模型从本构方程、几何方程和运动微分方程,推导其一阶应力速度方程,运用交错网格有限差分法对其进行数值模拟,研究该介质中地震波的波场特征与传播规律。同时分析了品质因子Q对地震波在振幅以及频率上的衰减、吸收作用的影响。在数值模拟中,使用交错网格算法必然会造成数值频散或假频现象,从而干扰数值模拟结果的正确性。为了提高数值模拟效果,对常规FCT(通量校正传输法)方法进行优化来压制频散,模拟实例证明,优化后的FCT方法要比常规FCT方法在消除频散方面更加有效,能有效地改善数值模拟精度。  相似文献   

5.
相比于VTI介质q P波数值模拟方法,考虑倾角因素的TTI介质q P波数值模拟方法能够更加准确地描述各向异性介质中波场的传播特征。常规的声学近似方法往往会造成TTI介质中倾角急剧变化的区域出现数值不稳定,笔者首先引入一个各向异性控制参数,推导了稳定形式的TTI介质二阶q P波方程,保证q P波场的稳定传播;其次,通过引入波场的伪速度分量,推导了其等价的一阶应力—速度形式及相应的PML边界条件,并应用优化的旋转交错网格有限差分(RSGFD)方法实现了精确的数值模拟。数值结果表明:TTI介质一阶q P波方程能够稳定、有效地模拟q P波的运动学特征,利用优化的RSGFD方法可以得到精确的合成地震记录,同时可以相对地提高计算效率。  相似文献   

6.
井间地震方法观测到的波场信息丰富且复杂,其难点是如何识别和分离这些波场信息。这里从纵波波场为无旋场,横波波场为无散场出发,利用一阶速度~应力方程的交错网格高阶有限差分法,对井间弹性波纵、横波场分离问题进行了数值模拟。经数值模拟试验表明,同二阶的位移方程相比,使用等价的波场分离交错网格的一阶速度应力方程,无论是模拟精度还是模拟的稳定性方面都有了较大的提高,同时该横波资料的获取方法为井间地震资料的解释应用,提供了一种新的技术思路。  相似文献   

7.
应用高阶交错网格有限差分算法,并加入吸收边界条件和衰减带,对弹性波方程进行模拟,分析了其稳定性和收敛性。通过对各向同性和各向异性介质模型的模拟表明,高阶差分波动方程模拟的网格频散较小、精度较高、效果较好,可为地震勘探及其资料解译提供技术手段。  相似文献   

8.
在井间地震有限差分数值模拟中,用离散化的高阶差分方程近似连续导数的波动方程时,不可避免地会产生数值频散,而数值频散程度则直接影响到地震波数值模拟精度,因此为了得到清晰准确的地震波场记录,必须尽可能地压制数值频散。这里在一阶速度应力弹性波方程的基础上,利用两个约束条件构造拉格朗日函数获取优化差分系数,与泰勒展开差分系数下的交错网格高阶差分模拟结果比较,发现改进的优化交错网格差分算子的高阶差分数值模拟能更有效地压制数值频散,进一步提高交错网格高阶差分数值模拟的精度,为高精度井间地震数据的波场成像、纵横波联合解释等提供可靠依据。  相似文献   

9.
利用Taylor公式展开导出了交错网格中可导函数任意次导数的任意偶数阶精度的差分近似式及相应的差分系数,从而改进了常规高精度交错网格有限差分算法中对各个空间导数采用不一致精度的问题。采用推导的交错网格高阶差分算法对一阶弹性波动方程进行数值模拟,得到了精确的数值模拟结果,证明了推导的交错网格高阶差分算法的正确性。  相似文献   

10.
从具有水平对称轴的横向各向同性(HTI)介质中的弹性波动方程出发,在交错网格空间中采用高阶差分算子对弹性波动方程进行差分离散,得到了HTI介质中地震波正演的高阶有限差分格式,研究并实现了PML吸收边界条件。在此基础上实现了HTI介质中弹性波方程的多波正演。数值算例表明,该方法能够精确模拟弹性波在复杂各向异性介质中的传播过程,得到高精度的正演记录。  相似文献   

11.
这里实现了在2.5维各向异性介质中地震波传播的数值模拟。首先从2.5维一阶速度-应力弹性各向异性波动方程出发,得出了该方程的拟谱法数值解法;然后通过对Marmousi速度模型进行数值模拟,表明了拟谱法的有效性;最后计算并分析了均匀和混合各向异性介质中波的传播快照和理论模拟的地震记录,进一步认识了波在各向异性介质中的传播规律。  相似文献   

12.
隋竞函  刘财 《世界地质》2018,37(4):1239-1249
基于一阶速度-应力波动方程,采用高阶交错网格有限差分数值模拟方法,对弹性及黏弹性TTI介质进行正演数值模拟。模拟时采用完全匹配层吸收边界条件(PML)消除边界反射。同时设计了层状介质模型、断层模型,通过模型的正演计算,得到了不同时刻的地震波波场快照及合成地震记录,分析其波场运动学及动力学特征。模拟结果表明,交错网格有限差分法可以很好地完成对复杂介质的波场模拟,具有较高的精度和可靠性。  相似文献   

13.
数值模拟是研究地震波传播机理的重要工具,也是地震资料采集、处理、和反演领域的重要内容。基于常规矩形网格或交错网格的差分算法存在模拟精度和计算量之间的矛盾,变网格技术可以缓解这一矛盾,但对于TTI介质,常规变网格、变系数差分算法无法实现各分量的交错对应,旋转网格与旋转交错网格技术存在同样的问题。针对上述问题,提出了一种适用于TTI介质正演的新方法:即菱形交错网格剖分法,并给出了基于菱形交错网格的变系数、变网格方法,实现了TTI介质中弹性波场的高精度模拟。  相似文献   

14.
根据各向异性介质的扰动理论来计算走时。首先简述各向异性介质的扰动理论,并推导非退化和退化情况的相速度扰动公式,而走时扰动可表示为沿未扰动介质中的射线路径的积分;当未扰动介质选为均匀各向性介质时,计算就较为简化。所给的数值例了证明了相速度扰动公式和走时计算公式的正确性。  相似文献   

15.
贴体网格有限差分正演模拟算法不仅能够精确模拟任意起伏地形下的波场特征,且计算效率较高,是一种很有应用前景的处理西部复杂地表问题的方法;然而,目前求解波动方程时常用的同位网格和标准交错网格,在处理贴体网格起伏地表正演模拟时存在诸多问题。为此,将全交错网格引入到曲线坐标系下,避免了标准交错网格的插值误差和同位网格中奇偶失联引起的高频振荡现象,提高了模拟精度,减小了算法实现的复杂度。在自由边界条件实施时,采用牵引力镜像法计算速度分量,速度自由边界条件配合紧致交错差分格式更新应力分量,得到了较好的效果。随后,重点研究了贴体全交错网格正演模拟算法的影响因素,考虑了网格正交性、网格间距和网格拼接等的影响,并取得了如下认识:算法对网格的正交性没有过分要求;网格间距的突变会引起虚假反射的产生;不同类型的网格拼接对模拟结果不会造成明显的影响。  相似文献   

16.
弱各向异性介质地震波传播特征分析   总被引:1,自引:1,他引:0  
随着油气田勘探精度要求越来越高,各向异性已成为不得不考虑的重要因素,因此对地震波在各向异性介质中传播特征的研究具有一定的价值。利用Thomsen提出的弱各向异性理论,通过理论模型数值计算,分析了Thomsen各向异性参数对弱各向异性VTI介质中三种体波的影响。分析结果表明:①参数ε、δ分别影响qP波大角度、中小角度传播的各向异性;②参数(ε-δ)影响qSV波中等角度传播时的各向异性;③参数γ影响qSH波大角度传播的各向异性,且群速度的各向异性要比相速度的各向异性更明显,并会出现复杂的交会区。通过这些研究分析,进一步认识了弱各向异性介质相速度及群速度随各向异性参数变化的规律。  相似文献   

17.
从一阶速度—应力弹性波方程出发,基于旋转交错网格,推导了时间二阶精度空间2M阶精度的有限差分离散格式。阐述了递归卷积复频移完全匹配层(CPML)边界条件的原理,建立了一阶速度—应力弹性波高阶差分CPML边界条件的递推公式。开展了CPML边界中关键参数m、κ和α的选取实验,通过分析反射误差分布图,选取了CPML边界条件中最优参数。全局反射误差与波场快照都说明,CPML较PML对隐失波具有更优的吸收性能。基于Matlab平台,编写了基于CPML边界的旋转交错网格弹性波正演模拟程序,应用该程序对各向异性介质及随机介质进行了模拟,得到了弹性波正演剖面记录及波场快照,通过对正演剖面记录及波场快照的分析,可以更清楚地了解弹性波在各向异性介质及随机介质的传播特性,指导非均匀介质中地震勘探资料解释。  相似文献   

18.
李敏  刘洋 《物探与化探》2012,36(6):934-940
笔者给出了一种能够模拟弹性波在任意各向异性介质中传播的二维三分量高阶有限差分算法.相对于常规交错网格有限差分方法,旋转交错网格有限差分方法在介质具有强差异性时能更精确地模拟地震波的传播,避免常规交错网格中因对弹性系数进行插值而带来的误差.采用高阶旋转交错网格有限差分方法模拟并分析了零偏移距横波分裂现象随裂缝介质方位角和倾角变化的响应特征.结果表明:结合完全匹配层(PML)吸收边界条件的高阶旋转交错网格有限差分方法能获得高精度的地震波场模拟数据,并且在边界具有良好的吸收效果;横波分裂现象主要受裂缝走向与波的极化方向之间的夹角影响,受裂缝倾角影响较小,且快慢横波的能量也跟裂缝走向与波极化方向间的夹角有关.具有倾斜对称轴的横向各向同性(TTI)介质倾角的变化可能会导致记录中波到达时的变化,影响快慢横波的时差.利用横波分裂的能量分布和方位各向异性特征,可以帮助检测裂缝的方位角和倾角.横波在多层TTI介质中传播时会发生多次分裂的现象.  相似文献   

19.
根据各向异性介质的扰动理论来计算走时。首先简述各向异性介质的扰动理论,并推导非退化和退化情况的相速度扰动公式,而走时扰动可表示为沿未扰动介质中射线路径的积分;当未扰动介质选为均匀各向同性介质时,计算就较为简化。所给的数值例子证明了相速度扰动公式和走时计算公式的正确性。  相似文献   

20.
以Biot双相介质模型为背景,笔者推导了双相各向同性介质二维三分量一阶速度——应力弹性波方程方程,建立了各向同性双相介质波动方程的二维三分量有限差分格式。分别采用传统交错网格有限差分技术和旋转交错网格有限差分技术对均匀和非均匀双相各向同性介质进行了波场模拟。结果表明,旋转交错网格有限差分技术能够有效模拟双相各向同性介质中弹性波的传播情况;通过传统和旋转交错网格有限差分技术的对比,说明了旋转交错网格有限差分算法的稳定性更强,避免了插值带来的误差,是一种有效的地震波场模拟方法。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号