共查询到18条相似文献，搜索用时 93 毫秒
1.
这里实现了在2.5维各向异性介质中地震波传播的数值模拟。首先从2.5维一阶速度-应力弹性各向异性波动方程出发,得出了该方程的拟谱法数值解法;然后通过对Marmousi速度模型进行数值模拟,表明了拟谱法的有效性;最后计算并分析了均匀和混合各向异性介质中波的传播快照和理论模拟的地震记录,进一步认识了波在各向异性介质中的传播规律。 相似文献
2.
采用交错网格技术将速度-应力方程中的速度对时间的导数转化为应力对空间的导数, 将弹性波动方程表示为与二阶双曲方程等价的一阶应力-速度公式, 以实现三维三分量地震波场模拟。对获得的波场快照、VSP记录图像和地面记录图像进行波场分析, 发现在波场图像中存在明显的拟P波、拟快横波和慢横波, 还出现了横波分裂、横波分裂盲区及波面三分叉等地震波在方位各向异性介质传播时产生的特殊现象。研究结果表明, 用交错网格三维模拟方法研究方位各向异性介质对方位各向异性介质中地震波传播反演能起到很强的辅助作用。 相似文献
3.
非均匀介质中交错网格高阶有限差分数值模拟 总被引：5，自引：0，他引：5
地震波场的数值模拟一直是地球物理学的一个重要的研究领域,而在数值正演模拟方法的研究中,计算精度和计算效率是评价该方法有效性及优越性的二个关键问题。这里从一阶速度—应力弹性波动方程出发,着重介绍如何构造离散化模型的网格,如何求解空间导数,如何选取边界条件等内容,从而更有效地提高数值计算的精度与计算效率。文中构造了不同类型的介质模型,并在交错网格中,利用高阶有限差分模拟非均匀介质的波场传播。模拟结果表明,该方法实现简单,具有很好地稳定性和较高的精度,能够直观、高效地反映出介质中波场的传播规律。 相似文献
4.
《物探化探计算技术》2019,(6)
为了提高二维直流电阻率法的数值模拟精度,这里实现了2.5维各向异性介质点源直流电场异常场的数值模拟。首先推导出2.5维各向异性介质总场和异常场满足的变分问题;然后采用矩形内三角网格的方式进行单元剖分,三角单元内二次场值和异常电阻率均设计为线性变化。为了简化背景场的计算,将背景介质的电性主轴与观测坐标系的夹角设计为零,从而避免了直接计算任意各向异性介质电位的傅里叶变换步骤。同时对计算中点源波数的选择和适用情况进行了讨论。最后通过对简单层状模型的计算,与解析解进行对比,验证了本数值算法的正确性和可靠性。比较常规各向同性介质模拟方式,本数值模拟发现直流电阻率法对横向电阻率的变化更为明显。 相似文献
5.
6.
基于一阶速度-应力波动方程,采用高阶交错网格有限差分数值模拟方法,对弹性及黏弹性TTI介质进行正演数值模拟。模拟时采用完全匹配层吸收边界条件(PML)消除边界反射。同时设计了层状介质模型、断层模型,通过模型的正演计算,得到了不同时刻的地震波波场快照及合成地震记录,分析其波场运动学及动力学特征。模拟结果表明,交错网格有限差分法可以很好地完成对复杂介质的波场模拟,具有较高的精度和可靠性。 相似文献
7.
各向异性介质储层中波场特征的研究,有利于对多波资料的准确解释。在弹性波方程基础上,结合数值模拟技术与AVO技术,对多种介质模型和第二类含气砂岩储层模型进行数值模拟与波场特征描述。以表征介质各向异性程度的Thomsen参数为出发点,研究各参数对地震波场动力学和运动学特征的影响,并对第二类含气砂岩介质中模拟得到的正演记录,以及各向异性参数对波场的影响程度进行分析,进一步加深了对该类储层各向异性特征的认识,这为各向异性介质中的参数反演与地震资料解释提供了有力的参考依据。 相似文献
8.
在地震传播理论中,地震波场的正演数值模拟一直是研究的热点。而在正演数值模拟方法的研究中,计算精度和计算效率是评价此方法的有效性及优越性的二个关键问题。伪谱法在计算精度与计算效率方面优越性十分明显,常用的有Fourier变换法和Hartley变换法;虽然Hartley变换法在求取导数时较Fourier变换法复杂,但由于Fourier变换法的计算同时涉及复数的实部与虚部,在计算速度和占用内存方面不如Hartley变换法。这里详细地阐述了利用Hartley变换求解正交各向异性介质波动方程的数值方法,并模拟了多种正交各向异性介质模型,对比分析了地震波在各向同性介质与正交各向异性介质中的传播差异。结论表明,该方法能正确、高效、直观地反应地震波在正交各向异性介质中的传播规律。 相似文献
9.
二阶弹性波动方程高精度交错网格波场分离数值模拟 总被引：4，自引：0，他引：4
给出了一种等价的二阶弹性波动方程,以解决弹性波场中完全弹性波动方程不能完全分离耦合的纵、横波波场问题.应用高阶交错网格有限差分法求解该波动方程,并使用通量校正技术(FCT)进一步压制频散,采用均匀介质模型和层状介质模型进行波场分离数值试验,精确得到了混合波场、完全分离的纯纵波及纯横波波场.数值结果分析表明,本文方法在均匀介质情况下准确可靠,在分离后的纯纵、横波波场中可观察到较为丰富的能量转换信息,这对认识复杂弹性波的传播规律及弹性波理论具有重要意义. 相似文献
10.
煤储层各向异性地震波场模拟 总被引：2，自引：1，他引：1
对中国煤层的裂隙系统进行了系统分类;阐述了各类裂隙的主要成因。基于二维二分量各向异性弹性波方程交错网格高阶有限差分法数值模拟,对各向同性介质模型和高裂隙密度(含干裂隙和饱和水裂隙)煤层的二维层状模型进行了弹性波场数值模拟。结果表明:煤层含裂隙后与不含裂隙时相比,反射时差的变化已不足以识别,但由于波速的改变会使煤层顶底、界面的波阻抗发生变化,从而引起波的动力学特征的改变,AVO现象尤其明显。含气和含水裂隙煤层中在x分量和z分量表现为极性相反,而且含水裂隙煤层中的波速比含气裂隙煤层的波速要快。因此,应用各向异性理论研究煤层裂隙的特征是可能的。 相似文献
11.
基于散度和旋度的弹性波波场分离数值模拟方法 总被引：4，自引：1，他引：3
弹性波波场分离数值模拟方法是研究混合波场中纯纵波和纯横波波场传播规律的一种重要手段,通常采用算子分裂来实现,这不但增加了计
算量,而且还浪费了大量计算时间。因此提出了基于散度和旋度的弹性波波场分离数值模拟新方法,仅依赖于交错网格来计算出散度场（纯纵波）
和旋度场（纯横波）,从而提高计算效率。在人工截断边界处采用完全匹配层吸收边界条件（PML）。数值模拟结果表明,该方法能够从混合波场中
准确地分离出纯纵横波波场,且精度高,边界吸收效果好,计算速度快,计算过程稳定,可用于指导实际多波多分量地震资料处理,同时指出本文
算法对相互耦合的波型（如面波）无法实现准确分离。 相似文献
12.
应用双相介质波动方程,推导了双相横向各向同性介质(TI)中波动方程的有限差分格式,对双相TI介质中弹性波有限差分数值进行了模拟.结果表明,弹性波在双相TI介质中传播时,除了存在常规的快纵波(qP1)和横波以外,还存在慢纵波( qP2).并且慢纵波的速度明显小于快纵波,而且受耗散系数的影响衰减地很快,所以在实际中很难观测到慢纵波.快纵波在固相和流相中相位相同,而慢纵波在固相和流相中的相位相反.慢纵波在流相中振幅大,而在固相中的振幅较小. 相似文献
13.
地震波混合阶褶积算法模拟 总被引：6，自引：0，他引：6
基于正反傅立叶变换,提出了地震波模拟的混合阶褶积算法。该方法原理简单、易于实现,结合了有限差分与伪谱法的优点,具有较高的精度和计算效率,适用于地震波场的正演计算,同时还给出了二维地震波场的理论计算实例。计算结果表明,此算法模拟结果正确、精度高、速度快、能适应较为复杂地质模型,并且易于推广到各向异性介质中去。 相似文献
14.
TTI介质弹性波频率-空间域有限差分数值模拟 总被引：3，自引：0，他引：3
由周期性薄互层引起的VTI介质是研究比较广泛的一类各向异性介质。当VTI介质对称轴偏离垂向,本构坐标系与观测坐标系不重合时,会形成观测坐标系下的TTI介质。引入25点优化差分算子,推导出二维TTI介质频率域弹性波动方程;为压制边界反射,采用完全匹配层法吸收边界条件,并计算出优化差分系数;最后采用集中力源,模拟了弹性波在TTI介质中的传播过程。从波场快照和地面共炮记录可以看出,笔者采用的数值模拟算法能有效压制数值频散。TTI介质中的波场传播比较复杂,纵波传播相对稳定,横波波前的三分叉现象比较明显,并存在振幅奇异性。当VTI介质的对称轴偏转后,还会增加地面地震记录的复杂性。 相似文献
15.
榴辉岩的弹性波速评述 总被引：2，自引：0，他引：2
文中评述了榴辉岩的密度和高温高压下的纵波速度、速度各向异性、泊松比以及榴辉岩声软化现象的成因 ,着重介绍了榴辉岩的密度和波速对探讨岩石圈物质组成、莫霍界面性质、超高压岩石对实现壳幔物质交换的重要意义。榴辉岩的密度为 3 2～ 3 6 5g·cm-3 ,其中超高压榴辉岩具有更高的密度 (3 4～ 3 6 5g·cm-3 )。榴辉岩的纵波速度 (vp)在 1GPa时为 7 3～8 9km·s-1,其裂隙闭合压力可能高于 1GPa。榴辉岩的压力系数为 0 3～ 0 4km·s-1·GPa-1,温度系数为 - 3 4× 10 -4 km·s-1·℃ -1。它具有最弱的vP 各向异性 (<3% )。超高压榴辉岩的泊松比为 0 2 54～ 0 2 75。大别山榴辉岩的密度和波速研究表明 ,现今的大别山深部地壳可能依然存在榴辉岩 ,但数量应很少 ;大别山上地幔具有同超高压榴辉岩类似的弹性特征 ;拆沉作用是解释超高压榴辉岩折返机制的重要模式之一 ;榴辉岩的形成过程包含了壳幔物质循环作用 ,一部分榴辉岩已拆沉进入深部地幔 ,另一部分则快速折返至地壳内或通过其他构造作用进一步抬升、暴露地表。 相似文献
16.
扭转导波在锚固锚杆中传播的数值模拟 总被引：1，自引：0，他引：1
采用有限元数值模拟的方法研究了扭转导波在锚固锚杆中的传播性质。建立了自由锚杆和锚固锚杆的有限元模型,在锚杆顶端激发20～60 kHz扭转导波信号,计算得到导波在自由锚杆和锚杆锚固段的传播速度值与理论值吻合很好,证明了扭转导波数值模拟方法的有效性。数值模拟结果表明：随着激发波频率的增大,扭转导波在自由锚杆和锚杆锚固段中的衰减值均呈线性递增趋势,导波在锚杆锚固段上界面的反射回波逐渐减弱;扭转导波在锚杆锚固段的衰减值较大,无法在锚杆顶端采集到锚杆底端反射回波信号,所以扭转导波不适用于锚杆长度的检测;随着锚固介质弹性模量的增大,同一频率扭转导波在锚杆锚固段上界面的反射波逐渐增强。通过检测扭转导波在锚杆锚固段上界面的反射回波,可以确定锚固介质弹性模量的大小。 相似文献
17.
18.
弹性波CT技术在码头加固效果检测中的应用 总被引：3，自引：0，他引：3
青岛港是我国北方一个很重要的港口码头,但由于年久失修,码头的使用出现安全隐患。经多种修补方案比较,最终采用垂直注浆加固方案。为了检测码头高压注浆加固处理效果,采用弹性波CT方法并结合岩芯波速测试及强度试验对其进行检测。本文简要叙述了弹性波CT法的原理与测试方法,探讨了该法在灌浆效果检测中的应用。试验证实,弹性波CT技术在探测水下高压注浆效果检测方面有其独特的效果。 相似文献