瑞利波泄漏模式的模拟研究

瑞利导波相速度频散曲线有时会出现截止或截断现象,此时必须考虑泄漏模式波的存在,如何快速求取泄漏模式波的相速度大小是研究的核心问题．采用添加附加层的方法,使瑞利波实频散方程求根的范围包含所有的层速度区域,从而可用Abo—zena,Menke和Bixing Zhang等学者研究过的传递矩阵理论,利用二分法求解相速度大小,即用导波来模拟计算泄漏模式波相速度大小．从模拟计算结果来看,附加层法不仅能模拟出原模型的实数导波,而且还能计算出泄漏模式波的变化．由此可见,用附加层法来计算瑞利波的泄漏模式是可行的;同时,若将该方法应用于瑞利波勘探的反演,它不会增加未知数的个数。     

计算瑞利波频散曲线的快速矢量传递算法

基于轴对称柱面瑞利面波，得到了一种计算层状介质中瑞利面波频散曲线的快速矢量传递算法，频散方程类似Menke方法的F矢量上传形式，上传矩阵F为3个五阶矩阵的乘积形式，提高了计算速度，而且各矩阵的元素均为无量纲量且为实数值，避免了以往方法中同一矩阵中各元素的数量级相差较大、出现复数运算的缺陷，提高了计算的精度及稳定性。此算法方法也避免了高频数值精度丢失问题和高频数值溢出问题。     

基于起伏地形的瑞利面波有限差分数值模拟

利用瑞利面波的频散特性及其不同频率瑞利面波能量随探测深度衰减变化规律,能够用来进行浅地表的地质勘探和地球内部构造的研究。利用有限差分数值求解的算法,模拟起伏条件下的地下介质中的瑞利面波;建立速度-应力弹性波动方程,利用交错网格的差分方式,引入完全匹配层(Perfectly matched layer,简称PML)的人工边界条件,对几种简单的起伏地表模型进行正演模拟。结果表明交错网格有限差分方法有比较好的精度,PML边界条件能够对来自边界的虚反射有较好的吸收处理。     

正弦型波纹板导波传播的理论模型与数值计算

为了探究正弦型波纹板中超声导波的传播特性,基于经典的自由平板中导波频散方程的推导方法,引入波纹板的结构参数,对波纹板中SH波和Lamb波的频散方程进行了推导.采用二分法数值求解频散方程,得到了波纹板相速度和群速度频散曲线.将波纹板和平板的频散曲线进行了比较,研究了导波不同模态的群速度随结构参数的变化规律.同时对波纹板中超声导波传播特性进行了有限元仿真,分析了正弦型波纹板中导波的实际传播距离随结构参数的变化规律.     

弹性波传播的低数值频散波场模拟方法

压制数值频散以提高计算精度是检验地震波数值模拟方法的一个重要标准。基于弹性波传播方程,建立了低数值频散波场模拟的八阶FNRK方法。该方法以Runge-Kutta方法对时间导数进行三阶离散,以近似解析离散算子替代差分算子对空间偏导数进行八阶离散,结合通量校正传输技术消除离散后的数值频散。弹性波场模拟结果表明,与高阶有限差分方法相比,该方法能在压制数值频散方面具有明显的优势,计算精度提高,且适应于地震波在大规模复杂介质中传播的波场模拟。     

求解声波方程的高精度NAD方法及其波场模拟

通过引入高阶精度的近似解析离散算子,给出了一种求解声波方程的八阶NAD方法.数值误差分析和计算效率结果显示,与四阶LWC方法和八阶LWC方法相比,八阶NAD方法具有高数值精度、高计算效率和低存储量.应用NAD方法模拟地震波在复杂非均匀3层介质和Marmousi模型中的传播,数值结果表明该方法能有效压制数值频散,具有较强的地震波模拟适用能力.     

改进BISQ模型的双相介质地震波场数值模拟及频散校正

从改进BlSQ模型的双相介质所对应的速度-应力运动方程出发,构建2×2N阶交错网格有限差分模拟算法:同时,为压制模拟过程中的数值频散现象,采用通量校正传输(FCT)技术获得带FCT修正的交错网格有限差分模拟算法,对各向同性双相介质进行地震波场数值模拟.研究结果表明:(1)波场存在快纵波、慢纵波和横波等波场特征,并与模型的理论响应相符,说明交错网格模拟算法具有正确性和可行性;(2)采用FCT修正的交错网模拟算法能够有效压制数值频散,并保留真实的波场特征.     

面波有限频率法3-D灵敏度核函数的计算

高分辨率面波成像应克服射线理论的高频近似假设,考虑有限频率效应.该文给出了一种计算面没有限频率3-D灵敏度核函数的方法.该方法可直接计算球壳层状地球模型中面波振型特征值和特征函数,而不再需要采用球化平变换来近似计算球壳层状地球模型中面波频散参数.作者计算了相位延迟、振幅扰动和群到时延迟的3-D灵敏度核函数,并对面波有限频率效应进行了讨论.分析结果显示:面波相位延迟的3-D灵敏度核函数同2-D体波走时的核函数形状类似,说明面波具有二维传播特性.对相速度数据而言,Fresnel带方法适用于低频面波的反演,射线方法适用于高频面波的反演,而对于中等频率的面波反演应采用有限频率方法.频率的商低是相对于所要取得的分辨率而言的.在反演地壳上地幔细结构方面,群频散数据更具有优越性;然而,群到时延迟核函数的极强的旁瓣可能使得反演变得不稳定.     

岩石边坡失稳的瑞利面波频散响应特征

根据边坡失稳的地球物理模型,系统研究了岩石边坡失稳的瑞利面波基阶模的相速度、波数响应特征,总结了边坡从稳定→失稳演化过程的瑞利面波基阶模频散特征变化规律,从中提取边坡失稳的瑞利面波基阶模相速度、波数、频散曲线等特征标志,使瑞利面波法不仅能用于探测边坡失稳的滑动面或结构面,而且还能评价边坡的稳定性,对边坡失稳进行综合预测预报.     

非饱和多孔介质层中Love波频散特性

为建立更准确的地层反演模型,基于非饱和多孔介质本构关系,用连通率处理上面两层间的边界条件,理论推导出三层非饱和多孔介质中Love波的频散方程,计算分析多阶模态Love波频散曲线。结果表明,Love波第1阶模态没有截止频率,对2阶以上模态,阶数越高截止频率越大。各阶模态的Love波最大速度趋近半空间层的横波速度,而最小速度趋近第1层横波速度。Love波速曲线以第2层横波速度为界分为上下两簇,各模态均在第2层横波波速处衔接。模态不同,衔接处曲线的光滑程度不同,模态越高,光滑程度越好。对各参数的敏感性分析表明,孔隙度、饱和度、连通率、黏滞系数和渗透率的变化对Love波相速度影响敏感性依次减小。