共查询到20条相似文献,搜索用时 194 毫秒
1.
2.
基于一阶速度-应力波动方程,采用高阶交错网格有限差分数值模拟方法,对弹性及黏弹性TTI介质进行正演数值模拟。模拟时采用完全匹配层吸收边界条件(PML)消除边界反射。同时设计了层状介质模型、断层模型,通过模型的正演计算,得到了不同时刻的地震波波场快照及合成地震记录,分析其波场运动学及动力学特征。模拟结果表明,交错网格有限差分法可以很好地完成对复杂介质的波场模拟,具有较高的精度和可靠性。 相似文献
3.
在井间地震有限差分数值模拟中,用离散化的高阶差分方程近似连续导数的波动方程时,不可避免地会产生数值频散,而数值频散程度则直接影响到地震波数值模拟精度,因此为了得到清晰准确的地震波场记录,必须尽可能地压制数值频散。这里在一阶速度应力弹性波方程的基础上,利用两个约束条件构造拉格朗日函数获取优化差分系数,与泰勒展开差分系数下的交错网格高阶差分模拟结果比较,发现改进的优化交错网格差分算子的高阶差分数值模拟能更有效地压制数值频散,进一步提高交错网格高阶差分数值模拟的精度,为高精度井间地震数据的波场成像、纵横波联合解释等提供可靠依据。 相似文献
4.
5.
《物探化探计算技术》2015,(6)
利用Taylor公式展开导出了交错网格中可导函数任意次导数的任意偶数阶精度的差分近似式及相应的差分系数,从而改进了常规高精度交错网格有限差分算法中对各个空间导数采用不一致精度的问题。采用推导的交错网格高阶差分算法对一阶弹性波动方程进行数值模拟,得到了精确的数值模拟结果,证明了推导的交错网格高阶差分算法的正确性。 相似文献
6.
殷文 《吉林大学学报(地球科学版)》2008,38(1):144-151
在弹性波频率空间域有限差分数值模拟方面,差分网格及边界条件是影响弹性波模拟成功与否的关键,为了压制数值模拟中的网格频散,采用25点有限差分算子,建立了有限差分矩阵方程,且借鉴匹配层衰减边界条件思想,设计了弹性波频率空间域有限差分数值模拟算法。由于采用高阶有限差分法来提高差分格式的精度,将会导致计算量显著增加,为此,对频率空间域有限差分弹性波数值模拟方法,采用流水线技术与分治策略进行了并行算法研究,提高了计算效率,使得在合理的计算时间内更精确地模拟弹性波在弹性介质中的传播过程。 相似文献
7.
笔者给出了一种能够模拟弹性波在任意各向异性介质中传播的二维三分量高阶有限差分算法.相对于常规交错网格有限差分方法,旋转交错网格有限差分方法在介质具有强差异性时能更精确地模拟地震波的传播,避免常规交错网格中因对弹性系数进行插值而带来的误差.采用高阶旋转交错网格有限差分方法模拟并分析了零偏移距横波分裂现象随裂缝介质方位角和倾角变化的响应特征.结果表明:结合完全匹配层(PML)吸收边界条件的高阶旋转交错网格有限差分方法能获得高精度的地震波场模拟数据,并且在边界具有良好的吸收效果;横波分裂现象主要受裂缝走向与波的极化方向之间的夹角影响,受裂缝倾角影响较小,且快慢横波的能量也跟裂缝走向与波极化方向间的夹角有关.具有倾斜对称轴的横向各向同性(TTI)介质倾角的变化可能会导致记录中波到达时的变化,影响快慢横波的时差.利用横波分裂的能量分布和方位各向异性特征,可以帮助检测裂缝的方位角和倾角.横波在多层TTI介质中传播时会发生多次分裂的现象. 相似文献
8.
波动方程数值模拟是研究地震波传播机理的重要工具,有限差分求解波动方程是当前地震波数值模拟的主要方法之一。当地下介质中的地震波速度较低或地震波高频成分丰富时,常规有限差分技术常常产生严重的数值频散误差,这种误差会降低数值模拟的精度,影响对地震波传播机理的分析。为压制地震波数值模拟时产生的数值频散误差,提高波场模拟精度,提出了基于NAD算子的时间四阶精度波动方程差分格式。根据对应的差分格式,分析了该差分格式的数值频散关系。与常规四阶精度差分算法的频散曲线相比,基于NAD时间四阶精度差分方法不但能够实现时间频散的有效压制,同时其基于更多网格点的位移分量和位移梯度分量空间微分求解方法还能够实现空间频散的有效压制。另外在相同模型条件下,基于NAD算法的声波方程时间四阶差分解法可采用大网格对模拟空间进行差分离散,减少网格数,提高计算效率。 相似文献
9.
煤巷地震超前探测合成记录的数值模拟方法研究 总被引:2,自引:0,他引:2
地震数值模拟技术是地震勘探方法研究的重要手段之一。应用地震波反射法探测技术,对巷道前方地质异常进行预报预测,可以更加有效地保证煤矿井下生产的安全。这里对弹性波动方程应用交错网格高阶有限差分技术,对巷道中的地震勘探技术进行数值模拟。采用的交错网格高阶差分技术网格频散小,PML使边界条件的处理简单而高效。通过不同的模型,模拟和研究井下地震波的传播情况,并对转换波的情况进行了对比。结果表明,该方法具有实现简单,精度高,稳定性好的特点。 相似文献
10.
《物探化探计算技术》2015,(2)
将真空法和旋转交错网格相结合,实现了地下空洞高阶旋转交错网格波场数值模拟,克服了旋转交错网格在模拟自由边界条件时只在空间二阶差分时稳定的不足。通过和已有普通交错网格下改进的真空法对比,验证了该方法的正确性和稳定性。为了减弱模型人工边界产生的反射,采用了吸收效果较好的完全匹配层技术。模拟了不同模型参数组合情况下的地震波传播以及地震波与地下空洞的相互作用,通过对比不同参数组合的地震记录和波场快照,以及包含空洞和不含空洞地震记录和波场快照,分析了地震波在地下存在空洞时的传播特征以及震源子波频率、地下空洞大小和埋深对地表地震记录的影响。 相似文献
11.
12.
从具有水平对称轴的横向各向同性(HTI)介质中的弹性波动方程出发,在交错网格空间中采用高阶差分算子对弹性波动方程进行差分离散,得到了HTI介质中地震波正演的高阶有限差分格式,研究并实现了PML吸收边界条件。在此基础上实现了HTI介质中弹性波方程的多波正演。数值算例表明,该方法能够精确模拟弹性波在复杂各向异性介质中的传播过程,得到高精度的正演记录。 相似文献
13.
为了研究巷道内地震波的波场特征,识别出波场中的各种波,利用交错网格高阶有限差分正演方法模拟了煤矿巷道内点震源激发及平面波震源激发产生的弹性波场。通过波场分析得出煤矿巷道内的共炮点道集中除了有常规地面地震中的各种波外,还有在巷道内上下界面的多次反射波和多次折射转换波,是巷道内地震勘探的主要干扰波。利用平面波震源在巷道内激发产生的波场中,巷道下方的空洞在平面波震源激发的地震记录上,其多次反射绕射波互相叠合在一起,表现为“糖葫芦”状的明显异常特征。这些都为巷道内地震勘探资料处理和解释提供了重要依据。 相似文献
14.
利用高阶交错网格有限差分模拟Kelvin-Voigt黏弹性介质中传播的地震波,同时将完全匹配层吸收边界条件引入到其边界处理中。数值模拟结果表明,完全匹配层吸收边界效果好,高阶有限差分能模拟得到的黏弹性介质波场精度较高。对模拟的黏弹性波场进行分析,表明介质的粘滞性使地震反射波的能量变弱,高频衰减明显,并比低频衰减得快,主频向低频方向移动,有效频带变窄,即降低了地震波的分辨率;并且反射转换波比反射纵波要衰减得快;而且还随着传播距离的增加,其峰值频率也逐渐降低。通过数值模拟分析具有不同的粘滞系数介质对地震波的吸收和衰减,结果表明随着粘滞系数的增大,地下介质对地震波的吸收衰减更明显。 相似文献
15.
目前,场地影响评价方法可以分为理论法和经验法。地脉动水平分量与竖向分量谱比法属经验法。这种方法将地脉动水平分量与竖向分量谱比的最大值对应的频率和幅值分别看作是观测场地的卓越频率和放大因子。虽然该法在国内外很多工程实例中得到了理想的结果。但由于这种方法建立过程中应用了许多假设,所以很多地震工程学家对这种方法持怀疑态度。本文讨论了基岩地脉动水平分量与竖向分量谱比为1这一假定的合理性。并采用泊松弹性半空间模型,应用概率分析方法考虑地脉动的随机性,推导了地下体波斜入射弹性半空间自由表面形成地脉动的水平分量与竖向分量谱比表达式。通过计算论证体波斜入射泊松弹性半空间形成的地脉动水平分量与竖向分量谱比特性。结果表明:当地脉动主要成分为体波时,基岩表面地脉动水平分量与竖向分量谱比接近于1;当地脉动主要成分为面波时,基岩表面地脉动水平分量与竖向分量谱比介于0.54~0.79之间。对泊松基岩,该比值是0.68。 相似文献
16.
从双程声波方程出发,在交错网格空间中推导了地震波逆时延拓的高阶有限差分算子,依据最佳匹配层(PML)的方程分裂思路,得到了一阶声波方程的PML边界条件及其高阶差分格式,采用零时间成像条件和上行、下行波场互相关成像条件,实现了声波方程的叠后与叠前逆时深度偏移。逆时偏移对sigsbee_2b模型理论数据的偏移成像得到了满意效果。 相似文献
17.
18.
研究井间地震波场的形成过程以及波场的传播机理、规律,对于指导实际井间地震勘探有着重要的意义.基于具有垂直对称轴的横向各向同性(VTI)介质中的一阶准P波方程,应用具有无条件稳定性质的紧致交错网格隐式差分格式求解该方程.重点研究了紧致交错网格求解该方程的完全匹配层(PML)吸收边界条件,在此基础上实现了VTI介质中一阶准P波方程的井间地震波场的正演模拟.数值算例表明:紧致交错网格能精准模拟VTI介质中准P波的传播过程,得到高精度的正演结果.一阶准P波方程能以足够的精度描述VTI介质中准P波特征.完全匹配层吸收边界能有效地解决人工边界问题,是一种高效的边界吸收算法. 相似文献
19.
For the case of seismic waves in a more complex architecture porous medium wave field occurs slightly changed with changes in the structure, several parameters characteristics of the structure of the media which influence on various types of seismic wave propagation were studied. Firstly, the article establishes the isotropic elastic porous medium model, derive the corresponding elastic wave equation, and uses high-order staggered-grid finite difference method for forward modeling, the article also analyze the pore structure parameters such as porosity, viscosity and penetration influence on the wave field characteristics were also analyzed. At the same time, the study analyze the influence of porosity, permeability and viscosity on phase velocity and attenuation coefficient was analyzed in the stady. The results showed that the influence on attenuation coefficient was more sensitive than that on phase velocity. This study helps to deepen the understanding of seismic wave propagation in the practical system of complex medium. 相似文献
20.
基于BISQ模型双相各向同性介质中地震波数值模拟 总被引:9,自引:2,他引:9
基于BISQ机制,推导了双相各向同性介质中弹性波数值模拟方程,并借助交错网格有限差分方法进行了波场模拟。证实了三种波(快纵波、慢纵波、横波SV)的存在和波场特征;在合成地震记录中,观测列三种波的反射,由于各种波相互之间的转换,转换波在地震记录中显示也很清楚,这使得波场变得更为复杂化。 相似文献