共查询到20条相似文献，搜索用时 24 毫秒
1.
波动方程数值模拟是研究地震波传播机理的重要工具,有限差分求解波动方程是当前地震波数值模拟的主要方法之一。当地下介质中的地震波速度较低或地震波高频成分丰富时,常规有限差分技术常常产生严重的数值频散误差,这种误差会降低数值模拟的精度,影响对地震波传播机理的分析。为压制地震波数值模拟时产生的数值频散误差,提高波场模拟精度,提出了基于NAD算子的时间四阶精度波动方程差分格式。根据对应的差分格式,分析了该差分格式的数值频散关系。与常规四阶精度差分算法的频散曲线相比,基于NAD时间四阶精度差分方法不但能够实现时间频散的有效压制,同时其基于更多网格点的位移分量和位移梯度分量空间微分求解方法还能够实现空间频散的有效压制。另外在相同模型条件下,基于NAD算法的声波方程时间四阶差分解法可采用大网格对模拟空间进行差分离散,减少网格数,提高计算效率。 相似文献
2.
3.
煤巷地震超前探测合成记录的数值模拟方法研究 总被引：2，自引：0，他引：2
地震数值模拟技术是地震勘探方法研究的重要手段之一。应用地震波反射法探测技术,对巷道前方地质异常进行预报预测,可以更加有效地保证煤矿井下生产的安全。这里对弹性波动方程应用交错网格高阶有限差分技术,对巷道中的地震勘探技术进行数值模拟。采用的交错网格高阶差分技术网格频散小,PML使边界条件的处理简单而高效。通过不同的模型,模拟和研究井下地震波的传播情况,并对转换波的情况进行了对比。结果表明,该方法具有实现简单,精度高,稳定性好的特点。 相似文献
4.
《物探化探计算技术》2016,(1)
在地震勘探中,粘弹性各向异性介质相比弹性介质更能实际表征地下介质的性质,这里主要针对Kelvin模型从本构方程、几何方程和运动微分方程,推导其一阶应力速度方程,运用交错网格有限差分法对其进行数值模拟,研究该介质中地震波的波场特征与传播规律。同时分析了品质因子Q对地震波在振幅以及频率上的衰减、吸收作用的影响。在数值模拟中,使用交错网格算法必然会造成数值频散或假频现象,从而干扰数值模拟结果的正确性。为了提高数值模拟效果,对常规FCT(通量校正传输法)方法进行优化来压制频散,模拟实例证明,优化后的FCT方法要比常规FCT方法在消除频散方面更加有效,能有效地改善数值模拟精度。 相似文献
5.
6.
《物探化探计算技术》2015,(6)
利用Taylor公式展开导出了交错网格中可导函数任意次导数的任意偶数阶精度的差分近似式及相应的差分系数,从而改进了常规高精度交错网格有限差分算法中对各个空间导数采用不一致精度的问题。采用推导的交错网格高阶差分算法对一阶弹性波动方程进行数值模拟,得到了精确的数值模拟结果,证明了推导的交错网格高阶差分算法的正确性。 相似文献
7.
Matlab环境下瑞利波有限差分正演与曲线绘制 总被引：1，自引：0，他引：1
以瑞利波频散方程为出发点的Abo-Zena传递矩阵等方法,只能研究层状介质中瑞利波的传播特性,对于非层状介质,传递矩阵无能为力。因此,不得不考虑有限元、有限差分等方法。根据弹性动力学方程,采用交错网格有限差分方法对均匀弹性半空间介质进行全波场模拟,并在Matlab环境下实现编程计算,再现了瑞利波在近地表的传播状态,实现了地震剖面曲线的正振幅充填问题,从而获得了类似野外地震勘探的剖面记录。模拟结果表明,采用一阶差分格式所得到的地震记录存在较严重的数值频散,在采用有限差分法研究瑞利波"之"字形成因时,须尽量压制数值频散现象,才能获得更接近实际情况的频散曲线。 相似文献
8.
殷文 《吉林大学学报(地球科学版)》2008,38(1):144-151
在弹性波频率空间域有限差分数值模拟方面,差分网格及边界条件是影响弹性波模拟成功与否的关键,为了压制数值模拟中的网格频散,采用25点有限差分算子,建立了有限差分矩阵方程,且借鉴匹配层衰减边界条件思想,设计了弹性波频率空间域有限差分数值模拟算法。由于采用高阶有限差分法来提高差分格式的精度,将会导致计算量显著增加,为此,对频率空间域有限差分弹性波数值模拟方法,采用流水线技术与分治策略进行了并行算法研究,提高了计算效率,使得在合理的计算时间内更精确地模拟弹性波在弹性介质中的传播过程。 相似文献
9.
从具有水平对称轴的横向各向同性（HTI）介质中的弹性波动方程出发,在交错网格空间中采用高阶差分算子对弹性波动方程进行差分离散,得到了HTI介质中地震波正演的高阶有限差分格式,研究并实现了PML吸收边界条件。在此基础上实现了HTI介质中弹性波方程的多波正演。数值算例表明,该方法能够精确模拟弹性波在复杂各向异性介质中的传播过程,得到高精度的正演记录。 相似文献
10.
基于波动方程的地震波正演模拟能够准确揭示波在介质中传播的振幅、频率与相位变化,真实地反映出波的动力学特征。有限差分法因其精度高、计算效率高而成为波场数值模拟的一个重要方法。根据地震波的近似方程———声波方程,采用二阶有限差分网格数值算法,利用透明边界条件,设计典型地质模型,进行地面波场的地震记录模拟。通过稳定性条件分析、震源子波优选以及频散抑制,成功合成了倾斜地层模型与正断层模型的地面共炮点波场记录,直观地显示出直达波、反射波、绕射波等丰富的波场信息,获取各模型的波场快照。结果表明,有限差分法结合透明边界条件对于获取介质丰富的波场信息具有很大的优势。 相似文献
11.
井间地震正演模拟技术是研究地震波在井间传播规律的重要手段之一,可以帮助认识井间地震的复杂波场.从二维井间地震波传播波动方程出发,结合初始、边界条件,推导出了交错网格任意偶阶精度差分格式,阐述了非均匀差分网格的实现方法.在此基础上,编制了跨孔波场交错网格有限差分正演程序,并应用该程序对半空间粘弹性跨孔模型进行了正演计算,得到了震源在不同位置时的波场快照及地震记录图.这对于跨孔模型波传播规律的研究,具有重要意义. 相似文献
12.
潘海滨 《物探化探计算技术》2009,31(4):369-373
交错网格高阶差分解法是地震波场模拟的一种有效方法。较大的时间延拓步长可以提高效率,但可能引起算法不稳定,使模拟无法进行。因此,对于空间网格大小和时间步长的选择要适中。虽然震源子波的主频对模拟精度的影响不大,但适当提高震源子波的主频,可以提高地震波对于薄层的分辨率。FCT方法能够有效地压制数值频散,将FCT方法与交错网格有限差分法相结合,可以提高波场模拟的精度和运算效率。 相似文献
13.
交错网格有限差分算法以其高效、精确、实用等优点在地震波数值模拟中得到广泛应用。目前交错网格有限差分的精度已达到时间4阶、空间2N阶;然而在求空间三次导数时,差分格式实际上并未达到所谓的2N阶精度,而是采用了低阶的差分格式及差分系数,这样有利于提高大尺度空间正演时的计算效率;但从计算精度的角度考虑,有必要推导出准确的满足2N阶精度的交错网格有限差分格式及差分系数,以得到更高精度的正演结果。笔者利用Taylor公式展开首次推导出了可导函数任意次导数的任意偶数阶精度的差分近似式及相应的差分系数,从而完善了常规高精度交错网格有限差分算法。采用新推导的交错网格有限差分格式得到的正演波形与解析解进行了对比,证明了新推导的差分格式的正确性,并与常规差分格式的正演波形进行了比较,结果显示,新推导出的交错网格有限差分格式模拟结果稳定性好,精度更高。 相似文献
14.
利用高阶交错网格有限差分模拟Kelvin-Voigt黏弹性介质中传播的地震波,同时将完全匹配层吸收边界条件引入到其边界处理中。数值模拟结果表明,完全匹配层吸收边界效果好,高阶有限差分能模拟得到的黏弹性介质波场精度较高。对模拟的黏弹性波场进行分析,表明介质的粘滞性使地震反射波的能量变弱,高频衰减明显,并比低频衰减得快,主频向低频方向移动,有效频带变窄,即降低了地震波的分辨率;并且反射转换波比反射纵波要衰减得快;而且还随着传播距离的增加,其峰值频率也逐渐降低。通过数值模拟分析具有不同的粘滞系数介质对地震波的吸收和衰减,结果表明随着粘滞系数的增大,地下介质对地震波的吸收衰减更明显。 相似文献
15.
在基于波动方程的有限差分数值模拟中,会不可避免地出现数值频散（也称网格发散）问题。数值频散问题通常会给数值模拟的结果造成严重影响,因此在数值模拟中,应尽量设法消除这种现象。这里在前人的基础上,在基于各向同性介质的弹性波方程的数值模拟过程中,通过引入通量校正方法来解决数值频散问题。由数值模拟结果表明,该方法可以有效地消除数值频散现象,大大改善数值模拟的结果。 相似文献
16.
基于一阶速度-应力波动方程,采用高阶交错网格有限差分数值模拟方法,对弹性及黏弹性TTI介质进行正演数值模拟。模拟时采用完全匹配层吸收边界条件(PML)消除边界反射。同时设计了层状介质模型、断层模型,通过模型的正演计算,得到了不同时刻的地震波波场快照及合成地震记录,分析其波场运动学及动力学特征。模拟结果表明,交错网格有限差分法可以很好地完成对复杂介质的波场模拟,具有较高的精度和可靠性。 相似文献
17.
18.
TTI介质弹性波频率-空间域有限差分数值模拟 总被引：3，自引：0，他引：3
由周期性薄互层引起的VTI介质是研究比较广泛的一类各向异性介质。当VTI介质对称轴偏离垂向,本构坐标系与观测坐标系不重合时,会形成观测坐标系下的TTI介质。引入25点优化差分算子,推导出二维TTI介质频率域弹性波动方程;为压制边界反射,采用完全匹配层法吸收边界条件,并计算出优化差分系数;最后采用集中力源,模拟了弹性波在TTI介质中的传播过程。从波场快照和地面共炮记录可以看出,笔者采用的数值模拟算法能有效压制数值频散。TTI介质中的波场传播比较复杂,纵波传播相对稳定,横波波前的三分叉现象比较明显,并存在振幅奇异性。当VTI介质的对称轴偏转后,还会增加地面地震记录的复杂性。 相似文献
19.
研究井间地震波场的形成过程以及波场的传播机理、规律,对于指导实际井间地震勘探有着重要的意义.基于具有垂直对称轴的横向各向同性(VTI)介质中的一阶准P波方程,应用具有无条件稳定性质的紧致交错网格隐式差分格式求解该方程.重点研究了紧致交错网格求解该方程的完全匹配层(PML)吸收边界条件,在此基础上实现了VTI介质中一阶准P波方程的井间地震波场的正演模拟.数值算例表明:紧致交错网格能精准模拟VTI介质中准P波的传播过程,得到高精度的正演结果.一阶准P波方程能以足够的精度描述VTI介质中准P波特征.完全匹配层吸收边界能有效地解决人工边界问题,是一种高效的边界吸收算法. 相似文献
20.
非均匀介质中交错网格高阶有限差分数值模拟 总被引：5，自引：0，他引：5
地震波场的数值模拟一直是地球物理学的一个重要的研究领域,而在数值正演模拟方法的研究中,计算精度和计算效率是评价该方法有效性及优越性的二个关键问题。这里从一阶速度—应力弹性波动方程出发,着重介绍如何构造离散化模型的网格,如何求解空间导数,如何选取边界条件等内容,从而更有效地提高数值计算的精度与计算效率。文中构造了不同类型的介质模型,并在交错网格中,利用高阶有限差分模拟非均匀介质的波场传播。模拟结果表明,该方法实现简单,具有很好地稳定性和较高的精度,能够直观、高效地反映出介质中波场的传播规律。 相似文献