裂缝介质中瑞雷面波传播的渐变非均匀交错网格数值模拟

提出了渐变非均匀交错网格高阶差分方法,模拟裂缝介质中瑞雷面波的传播过程,并引入具有二阶时间精度和六阶空间精度的附加差分公式解决非均匀交错网格的不对称问题。这种差分方法可描述为：在裂缝处采用细网格,在其他地方采用粗网格,细网格到粗网格之间在空间上采用渐变步长,达到既细致地描述了局部细微结构、又保持了计算量适当的目的。在变步长的选取中提出了采用汉宁窗函数,这样既可以较快地从粗网格向细网格过渡,又避免了粗细网格过渡时发生频散,算例表明该方法完全适用于模拟裂缝介质中瑞雷面波的传播,计算精度较高,且能节约时间和内存,数值模拟结果非常直观地揭示了瑞雷面波在裂缝介质中的波场传播特征。     

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

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