利用一个六层模型(图7-1)分析品质因子QP和QS对瑞雷波衰减系数的贡献。令VS取VP的25%～50%的值,当VS/VP增加时,QP对瑞雷波衰减系数的贡献增加,QS对瑞雷波衰减系数的贡献减小(图7-2)。当VS/VP接近0.45时,QP的贡献对所有频率来说都变得显著...[继续阅读]

    1.示例模型示例模型的合成炮集记录见图4-20。模型参数见表4-1。图4-20　示例模型的合成炮集记录:60道,最小偏移距1m,道间距1m;采样率1000Hz,记录时间1.024s表4-1　二层理论模型参数2.程序部分代码%%%%%%相关参数%%%%%%%%%%seis　　　%地震数据...[继续阅读]

    传统单轴PML技术对于高泊松比(>0.38)浅地表地质模型,是不稳定的。泊松比越高,传统PML算法发散越快。自由表面的存在造成了数值模拟时的不稳定。自由表面复杂的波现象,是影响PML内数值误差快速累积的重要因素。对泊松比在0.1...[继续阅读]

    因为层数的多变性,多层模型情况有些复杂。采用较少层数的模型来代替多层模型来简化问题,表5-3就用来验证该想法的正确性。用三层模型来代替多层模型,将多层模型的频散曲线作为实测数据,并用Xia等(1999)讨论的反演方法将其反演...[继续阅读]

    瑞雷波具有频散的特性,且在同一频率有多个不同的相速度。在某一频率最小传播速度被称为基阶模式(一阶模式)(图6-27),相邻的高于基阶模式的相速度被称为一阶高阶模式(二阶模式),以此类推。所有高于基阶模式传播速度的模式均称...[继续阅读]

    与瑞雷波频散方程推导方式类似,我们可以通过连续两层上界面之间的关系,应力位移连续边界条件、自由边界条件及无穷远辐射条件得到勒夫波的频散函数:其中,当1≤m≤n-2时,式中:Qm＝krβmhm,rβm＝-i。当m=n-1时且半空间的上界面为...[继续阅读]

    当表层的横波速度高于某些下伏地层时,由目前算法得到的瑞雷波相速度在其高频部分会趋近除半空间以外的最低横波层速度,而不是与表层横波速度相关的某个值。基于波动方程的数值模拟结果表明,瑞雷波的相速度在其波长远小于...[继续阅读]

    1.反演方程的建立在瑞雷波的反演算法中,首先定义目标函数和权重矩阵,利用L-M方法和奇异值分解方法使得目标函数最小化。假设S波速度是一个n维长度的向量x,x=[VS1,VS2,VS3,…VSn]T,同样地,瑞雷波的相速度也可以用一个m维长度的向量...[继续阅读]

    前一节的实例告诉我们,通过反演基阶模式瑞雷波相速度通常能得出相对误差小于15%的S波速度。然而,当基阶波信号被体波或高阶模式波“污染”后,反演过程将会变得不稳定,求得的浅地表横波速度的准确性将降低。众所周知,地球物...[继续阅读]

    王周等(2012)总结求解弹性波有限差分法中自由边界常用的处理方法:直接法、应力镜像法、改进应力镜像法、横向各向同性介质替换法和声学边界替换法等。如图3-3,自由边界处垂直于边界的正应力为零,这是空气与地下介质的分界面...[继续阅读]