共查询到18条相似文献，搜索用时 156 毫秒
1.
交错网格有限差分算法以其高效、精确、实用等优点在地震波数值模拟中得到广泛应用。目前交错网格有限差分的精度已达到时间4阶、空间2N阶;然而在求空间三次导数时,差分格式实际上并未达到所谓的2N阶精度,而是采用了低阶的差分格式及差分系数,这样有利于提高大尺度空间正演时的计算效率;但从计算精度的角度考虑,有必要推导出准确的满足2N阶精度的交错网格有限差分格式及差分系数,以得到更高精度的正演结果。笔者利用Taylor公式展开首次推导出了可导函数任意次导数的任意偶数阶精度的差分近似式及相应的差分系数,从而完善了常规高精度交错网格有限差分算法。采用新推导的交错网格有限差分格式得到的正演波形与解析解进行了对比,证明了新推导的差分格式的正确性,并与常规差分格式的正演波形进行了比较,结果显示,新推导出的交错网格有限差分格式模拟结果稳定性好,精度更高。 相似文献
2.
基于一阶速度-应力波动方程,采用高阶交错网格有限差分数值模拟方法,对弹性及黏弹性TTI介质进行正演数值模拟。模拟时采用完全匹配层吸收边界条件(PML)消除边界反射。同时设计了层状介质模型、断层模型,通过模型的正演计算,得到了不同时刻的地震波波场快照及合成地震记录,分析其波场运动学及动力学特征。模拟结果表明,交错网格有限差分法可以很好地完成对复杂介质的波场模拟,具有较高的精度和可靠性。 相似文献
3.
殷文 《吉林大学学报(地球科学版)》2008,38(1):144-151
在弹性波频率空间域有限差分数值模拟方面,差分网格及边界条件是影响弹性波模拟成功与否的关键,为了压制数值模拟中的网格频散,采用25点有限差分算子,建立了有限差分矩阵方程,且借鉴匹配层衰减边界条件思想,设计了弹性波频率空间域有限差分数值模拟算法。由于采用高阶有限差分法来提高差分格式的精度,将会导致计算量显著增加,为此,对频率空间域有限差分弹性波数值模拟方法,采用流水线技术与分治策略进行了并行算法研究,提高了计算效率,使得在合理的计算时间内更精确地模拟弹性波在弹性介质中的传播过程。 相似文献
4.
波动方程数值模拟是研究地震波传播机理的重要工具,有限差分求解波动方程是当前地震波数值模拟的主要方法之一。当地下介质中的地震波速度较低或地震波高频成分丰富时,常规有限差分技术常常产生严重的数值频散误差,这种误差会降低数值模拟的精度,影响对地震波传播机理的分析。为压制地震波数值模拟时产生的数值频散误差,提高波场模拟精度,提出了基于NAD算子的时间四阶精度波动方程差分格式。根据对应的差分格式,分析了该差分格式的数值频散关系。与常规四阶精度差分算法的频散曲线相比,基于NAD时间四阶精度差分方法不但能够实现时间频散的有效压制,同时其基于更多网格点的位移分量和位移梯度分量空间微分求解方法还能够实现空间频散的有效压制。另外在相同模型条件下,基于NAD算法的声波方程时间四阶差分解法可采用大网格对模拟空间进行差分离散,减少网格数,提高计算效率。 相似文献
5.
井间地震正演模拟技术是研究地震波在井间传播规律的重要手段之一,可以帮助认识井间地震的复杂波场.从二维井间地震波传播波动方程出发,结合初始、边界条件,推导出了交错网格任意偶阶精度差分格式,阐述了非均匀差分网格的实现方法.在此基础上,编制了跨孔波场交错网格有限差分正演程序,并应用该程序对半空间粘弹性跨孔模型进行了正演计算,得到了震源在不同位置时的波场快照及地震记录图.这对于跨孔模型波传播规律的研究,具有重要意义. 相似文献
6.
在井间地震有限差分数值模拟中,用离散化的高阶差分方程近似连续导数的波动方程时,不可避免地会产生数值频散,而数值频散程度则直接影响到地震波数值模拟精度,因此为了得到清晰准确的地震波场记录,必须尽可能地压制数值频散。这里在一阶速度应力弹性波方程的基础上,利用两个约束条件构造拉格朗日函数获取优化差分系数,与泰勒展开差分系数下的交错网格高阶差分模拟结果比较,发现改进的优化交错网格差分算子的高阶差分数值模拟能更有效地压制数值频散,进一步提高交错网格高阶差分数值模拟的精度,为高精度井间地震数据的波场成像、纵横波联合解释等提供可靠依据。 相似文献
7.
非均匀介质中交错网格高阶有限差分数值模拟 总被引：5，自引：0，他引：5
地震波场的数值模拟一直是地球物理学的一个重要的研究领域,而在数值正演模拟方法的研究中,计算精度和计算效率是评价该方法有效性及优越性的二个关键问题。这里从一阶速度—应力弹性波动方程出发,着重介绍如何构造离散化模型的网格,如何求解空间导数,如何选取边界条件等内容,从而更有效地提高数值计算的精度与计算效率。文中构造了不同类型的介质模型,并在交错网格中,利用高阶有限差分模拟非均匀介质的波场传播。模拟结果表明,该方法实现简单,具有很好地稳定性和较高的精度,能够直观、高效地反映出介质中波场的传播规律。 相似文献
8.
弹性波模拟或逆时偏移时,对空间偏导数采用高阶差分格式可提高计算精度,但这种算法的稳定性条件过于严格,要求差分离散的时间步长必须足够小以确保算法稳定。在常规空间高阶差分格式的基础上,将速度(应力)对时间的高阶导数转化为不同精度的应力(速度)对空间的差分,得到了一种新的基于交错网格的时间高阶、空间高阶差分格式。通过对交错网格时间高阶差分格式稳定性的分析,认为该算法的稳定性条件较常规算法宽松,在弹性波场的求解过程中可以采用更大的时间步长。 相似文献
9.
10.
从一阶速度—应力弹性波方程出发,基于旋转交错网格,推导了时间二阶精度空间2M阶精度的有限差分离散格式。阐述了递归卷积复频移完全匹配层(CPML)边界条件的原理,建立了一阶速度—应力弹性波高阶差分CPML边界条件的递推公式。开展了CPML边界中关键参数m、κ和α的选取实验,通过分析反射误差分布图,选取了CPML边界条件中最优参数。全局反射误差与波场快照都说明,CPML较PML对隐失波具有更优的吸收性能。基于Matlab平台,编写了基于CPML边界的旋转交错网格弹性波正演模拟程序,应用该程序对各向异性介质及随机介质进行了模拟,得到了弹性波正演剖面记录及波场快照,通过对正演剖面记录及波场快照的分析,可以更清楚地了解弹性波在各向异性介质及随机介质的传播特性,指导非均匀介质中地震勘探资料解释。 相似文献
11.
研究井间地震波场的形成过程以及波场的传播机理、规律,对于指导实际井间地震勘探有着重要的意义.基于具有垂直对称轴的横向各向同性(VTI)介质中的一阶准P波方程,应用具有无条件稳定性质的紧致交错网格隐式差分格式求解该方程.重点研究了紧致交错网格求解该方程的完全匹配层(PML)吸收边界条件,在此基础上实现了VTI介质中一阶准P波方程的井间地震波场的正演模拟.数值算例表明:紧致交错网格能精准模拟VTI介质中准P波的传播过程,得到高精度的正演结果.一阶准P波方程能以足够的精度描述VTI介质中准P波特征.完全匹配层吸收边界能有效地解决人工边界问题,是一种高效的边界吸收算法. 相似文献
12.
地震波传播的褶积微分算子法数值模拟 总被引：2，自引：0，他引：2
为了解决地震波场数值模拟中的速度与精度的匹配及局部信息与全局信息的耦合等问题, 该文借助新推出的基于Forsyte广义正交多项式的褶积微分算子, 将计算数学中的Forsyte多项式褶积微分算子应用到地震波传播的数值模拟中.复杂非均匀介质模型数值模拟结果说明了该方法的可行性和优越性.该方法同时具有广义正交多项式褶积微分算子的高精度和有限差分短算子算法的高速度.通过对算子长度的调节及算子系数的优化, 可同时兼顾波场解的全局信息与局部信息. 相似文献
13.
采用交错网格技术将速度-应力方程中的速度对时间的导数转化为应力对空间的导数, 将弹性波动方程表示为与二阶双曲方程等价的一阶应力-速度公式, 以实现三维三分量地震波场模拟。对获得的波场快照、VSP记录图像和地面记录图像进行波场分析, 发现在波场图像中存在明显的拟P波、拟快横波和慢横波, 还出现了横波分裂、横波分裂盲区及波面三分叉等地震波在方位各向异性介质传播时产生的特殊现象。研究结果表明, 用交错网格三维模拟方法研究方位各向异性介质对方位各向异性介质中地震波传播反演能起到很强的辅助作用。 相似文献
14.
地震波混合阶褶积算法模拟 总被引：6，自引：0，他引：6
基于正反傅立叶变换,提出了地震波模拟的混合阶褶积算法。该方法原理简单、易于实现,结合了有限差分与伪谱法的优点,具有较高的精度和计算效率,适用于地震波场的正演计算,同时还给出了二维地震波场的理论计算实例。计算结果表明,此算法模拟结果正确、精度高、速度快、能适应较为复杂地质模型,并且易于推广到各向异性介质中去。 相似文献
15.
地震勘探广泛应用于油气、煤田勘探。地震波场数值模拟是整个地震勘探数据处理技术的基石。将径向基函数（RBF）引入地震声波波场数值模拟中,在空间上用径向基函数无网格法来构造二阶导数,而在时间上采用简单的二阶差分公式,并重点讨论了形状参数c对该方法精度的影响,总结c经验取值范围为2~4倍平均数据点间距。设计不同模型,利用径向基函数无网格法进行声波波场模拟,并与空间四阶时间二阶的有限差分计算结果进行对比,结果表明:同样精度下,径向基函数每个波长所取的数据点数远小于空间四阶矩形网格有限差分每个波长所取的网格点数,即径向基函数的空间采样率更低,这表明径向基函数具有更小的数值频散。 相似文献
16.
Waveform modelling is essential for seismic imaging and inversion. Because including more physical characteristics can potentially yield more accurate Earth models, we analyse strategies for elastic seismic wave propagation modelling including topography. We focus on using finite differences on modified staggered grids. Computational grids can be curved to fit the topography using distribution functions. With the chain rule, the elasto-dynamic formulation is adapted to be solved directly on curved staggered grids. The chain-rule approach is computationally less expensive than the tensorial approach for finite differences below the 6th order, but more expensive than the classical approach for flat topography (i.e. rectangular staggered grids). Free-surface conditions are evaluated and implemented according to the stress image method. Non-reflective boundary conditions are simulated via a Convolutional Perfect Matching Layer. This implementation does not generate spurious diffractions when the free-surface topography is not horizontal, as long as the topography is smoothly curved. Optimal results are obtained when the angle between grid lines at the free surface is orthogonal. The chain-rule implementation shows high accuracy when compared to the analytical solution in the case of the Lamb’s problem, Garvin’s problem and elastic interface. 相似文献
17.
Zengli Du Jianjun Liu Jinping Liu Feng Xu Yongzhang Li 《Arabian Journal of Geosciences》2017,10(21):473
Since the precision of spatial derivative of wavefield calculated by staggered difference is superior to that by central difference, the formula of second-order staggered differential coefficient with 2N-order precision is constructed via mathematical deduction on combination of first-order differential coefficient on each grid point, and its validity is verified by numerical simulation in constant velocity media. With the staggered mesh, the maximum peak frequency of Ricker wavelet which does not generate numerical dispersion may be increased by 10–15 Hz, thus the vertical resolution of simulated seismograms may be enhanced. In order to eliminate the boundary reflection which is characterized as high frequency, the absorbing boundary condition of perfectly matched layer with staggered mesh is utilized in numerical simulation of scalar-wave equation, and the suppression effect is superior to that of the central difference. The precision of simulated records may be improved when staggered mesh is adopted in calculation of the second-order spatial derivative, and the ringing phenomena appeared before the first break may be avoided. The computational cost of second-order staggered difference increased by only 25% on 3.60 GHz Intel Core i7–4790 processor when the precision of simulated records is same as that simulated by second-order central difference. 相似文献
18.
针对起伏地表条件下的地震波数值模拟问题,提出将起伏的地表映射到一个规则的长方形网格坐标系中,并在此基础上推导出变换域中的波动方程。根据导出来的波动方程应用有限余弦变换有限差分方法进行地震波场的数值模拟,以解决起伏地表条件下的数值模拟。 相似文献