《二维各向同性弹性介质瑞雷波波场数值模拟》.pptVIP

中国地质大学(武汉) 地球物理与空间信息学院 交流提纲 一、引 言 二、弹性介质瑞雷波波场数值模拟 弹性波波动方程的建立 数值模拟各模块说明 数值模拟算法实现 三、波场数值模拟中的思考 最小偏移距的优化问题 地质模型泊松比分析 Rayleigh Wave Method具有快速、无损、衰减小、抗干扰能力强、浅层分辨率高，且不受各地层速度关系的影响等优点，目前已经成为浅地表地球物理工程的一个研究热点，并且在环境与工程地球物理、能源勘探地球物理等领域都有了广泛的应用。通过波场数值模拟的方法来研究瑞雷波的传播特性和频散特性，将大大提高瑞雷波勘探推断地下复杂介质结构的能力，并能进一步拓宽瑞雷波勘探应用领域。 本文的主要内容：通过交错网格高阶有限差分的数值模拟方法，实现了对二维各向同性弹性介质中的瑞雷波的波场数值模拟，并将正演的地震记录转化为频散能量图进行分析。 一、引 言 二、弹性介质瑞雷波波场数值模拟 弹性波波动方程是研究地震波传播的出发点，根据地震学理论，弹性波波动方程可由柯西方程、本构方程和运动平衡微分方程来确定。本文只考虑二维P-Sv波情况，即不考虑y方向。 弹性波波动方程的建立 柯西方程(Cauchy Equation)（描述应变和位移之间的关系） 其中 和 为正应变分量， 为剪应变分量， 和 为质点在 ， 方向的位移分量。 本构方程(广义虎克定律)（描述应力和应变之间的关系） 其中， 为弹性系数张量。 在均匀各向同性介质中，弹性系数可用拉梅常数 和 表示。其矩阵表示形式为： 弹性波波动方程的建立(续) 其中 和 为质点在 和 方向的正应力， 为质点在 平面内的剪切应力。 运动平衡微分方程(Navier Equation)（描述位移和应力的关系） 二、弹性介质瑞雷波波场数值模拟 弹性波波动方程的建立(续) 上式是不考虑外力作用影响的方程。其中， 为时间， 为各向同性介质的密度。 联立柯西方程、本构方程和运动平衡微分方程，可以得到弹性各向同性介质中的一阶速度-应力波动方程： 其中， 和 为拉梅常数，则拉梅常数可以表示为： 式中， 和 分别为纵波和横波速度， 为各向同性介质的密度。 二、弹性介质瑞雷波波场数值模拟 数值模拟各模块说明 数值模拟方法：采用差分近似为时间上2阶、空间上12阶(2,12)的交错网格有限差分； 自由边界条件：采用Xu等(2007)提出的声学边界替代法(AEA)作为自由边界条件； 人工边界条件：采用完全匹配层(PML)吸收边界条件； 震源条件：采用雷克子波作为震源函数，采用垂直集中力震源作为加载方式； 稳定性条件：采用董良国等(2000)在一阶弹性波方程交错网格高阶有限差分解法的稳定性的讨论中提出的2N阶差分稳定性条件。 二、弹性介质瑞雷波波场数值模拟 数值模拟算法实现 参数说明：网格大小为 ，空间步长为1m，时间间隔为0.1ms；震源采用雷克子波，主频为20Hz，时间延迟为30ms。PML吸收层网格数10，吸收层厚度10m，吸收系数0.01；道集数48，最小偏移距10m，道间距1m，排列长度47m。采用差分近似为时间2阶、空间12阶的交错网格有限差分。 层序号 横波速度 (m/s) 纵波速度 (m/s) 密度 (g/cm3) 层厚度 (m) 1 200 346 1.9 Infinite 表1 均匀半空间地质模型物性参数表 均匀半空间模型 二、弹性介质瑞雷波波场数值模拟 均匀半空间模型(续) 图1 300ms时垂直分量波场快照 瑞雷波传播速度分析：根据图1的波场快照，可以看出纵波速度最快，横波速度其次，瑞雷波速度最慢，且瑞雷波速度跟横波速度很接近。从图2中可以计算出瑞雷波速度约为183m/s，而模型给出的横波速度为200m/s，由此可知，瑞雷波传播速度略小于横波速度。 图2 均匀半空间模型地震记录 二、弹性介质瑞雷波波场数值模拟 均匀半空间模型(续) 图3 均匀半空间模型频散能量图 瑞雷波频散分析：图3是由高分辨率线性拉东变换得到的频散能量图，图中低频段出现喇叭口形状，这是因为在进行高分辨率线性拉东变换过程中，进行速度扫描时由于低频段数据对应的速度范围较大造成的，其他频段频散曲线基本成一条直线，不随频率的改变而改变，这表明，在均匀半空间中瑞雷波不发生频散。 瑞雷波传播深度分析：瑞雷波只沿着地表传播，是一种面波，由图1可看出，瑞雷波在模型空间中的传播深度约为10m。在本次数值模拟中，所用震源子波主频为 ，