简介:瑞利波勘探主要是建立在弹性介质分层半空间模型基础上的。当实际地层中包含孔隙介质层时,需要将孔隙介质简化为弹性介质来进行分析和处理。这种简化处理究竟对瑞利波勘探造成怎样的影响是本文所研究的问题。本文基于弹性介质与孔隙介质共同构成的分层半空间模型,首先推导了分层半空间中两种介质处于不同相对位置时的瑞利波频散方程,解决了不同阶数矩阵之间的变量传递问题然后,针对传递矩阵法在求解频散函数时可能出现的溢出问题,给出了一种可以有效提高计算范围的解决方案;同时,提出了一套新的数值算法用于复频散方程的快速求解。数值计算结果表明:当孔隙介质位于半空间表面时对低频瑞利波频散特性的影响最为显著,而在其它情况下对瑞利波频散特性的影响相对较小。
简介:频率空间域地震波数值模拟具有独特的优势:可以同时模拟多源的波传播、每个频率之间独立并行地计算、计算频带选择灵活、不存在累计误差、容易模拟粘弹性介质中地震波传播。但是该方法的最大瓶颈是对于计算机内存的巨大需求。我们使用压缩存储系数矩阵的方法,极大地减少了计算机内存的需求量。同时为了减少短差分算子的数值频散,引用了频率空间域25点弹性波波动方程的差分格式,并使用了最小二乘意义下求出的优化差分系数。为了克服边界反射,采用了最佳匹配层吸收边界条件。数值模拟试验证明:用压缩存储系数矩阵及优化差分系数的频率空间域25点差分格式进行弹性波正演模拟,可以减少数值频散,提高计算精度。使用较大的网格间距,降低计算机内存需求,并保持较高的计算效率。该正演方法为后续弹性波偏移和弹性参数反演提供较好的基础。