各向异性介质地震波场数值模拟及特征分析-ct理论与应用研究.pdfVIP

第26 卷 第6 期 CT 理论与应用研究 Vol.26, No.6 2017 年12 月 （661-668） CT Theory and Applications Dec., 2017 龚明平, 张军华, 刘文革, 等. 各向异性介质地震波场数值模拟及特征分析[J]. CT理论与应用研究, 2017, 26 (6): 661-668. doi:10.15953/j.1004-4140.2017.26.06.01. Gong MP, Zhang JH, Liu WG, et al. The Seismic wave numerical simulation and characteristics analysis in anisotropic medium [J]. CT Theory and Applications, 2017, 26(6): 661-668. (in Chinese).doi:10.15953/j.1004-4140.2017.26.06.01. 各向异性介质地震波场数值模拟及特征分析 1 1 2 1 1 龚明平 ，张军华 ，刘文革 ，刘杨 ，吴系源 1.中国石油大学（华东）地球科学与技术学院，山东 青岛266580 2.西南石油大学地球科学与技术学院，成都610500 摘要：由于各向异性广泛存在于地下岩石中，随着勘探精度的不断提高，对地下介质的各向同 性假设越来越不能够满足于现状，因此对各向异性介质的数值模拟显得更为重要。本文推导了 各向异性介质的弹性波动方程，总结了震源类型，通过 PML 方法处理了人工边界问题，通过快 照分析验证了数值频散、稳定性条件。研究结果表明：①PML 完全匹配层，可较好地解决人工 边界问题；②减小空间采样间隔压制数值频散比减小时间采样间隔效果要好得多，盲目减小时 间采样间隔会大大降低数值模拟的运算效率；③各向异性介质中弹性波场中除含有准纵波外， 还含有速度较慢的准横波；④准纵波波前能量要比由各向异性引起的准横波能量强，准纵波和 准横波的波前随着各向异性介质参数的变化而变化。 关键词：各向异性介质；数值模拟；有限差分；数值频散 doi:10.15953/j.1004-4140.2017.26.06.01 中图分类号：O 242；P 315 文献标志码：A 勘探地震学常将单元体最大尺度小于地震频率范围内的最小波长，且表现出弹性性质 随方向而改变的介质称为各向异性介质。通常各向异性介质中纵波的偏振（或极化）方向 并不与传播方向平行，横波的偏振方向也不与传播方向垂直，因此通常把横向各向同性介 [1] 质(TI)中的纵波称为qP 波，把两种横波分别称为qSV 波和SH 波 。 对各向异性的研究可以追溯到Kelvin[2] （1856）发表的、具有现代形式的关于弹性各向异 性的第一篇论文。该文最早提出了弹性张量的无坐标表示方法，第一次提出了各向异性介质的 表示方法。Helbig[3] （1863）研究了各向异性介质中格林函数的运用。Backus[4] （1965）提出 [5] 了弱各向异性介质中速度与岩性的关系式。宋海滨等 （1997）进行了层状横向各向异性介质 [6] 的反问题探究。董国良等 （1999）对横向各向同性介质弹性波的物理模型进行了模拟。程玖 [7] [8] 兵等 （2001）研究了波动方程共炮检距道集叠前深度偏移。裴正林等 （2006）模拟了交错 [9] 网格高阶有限差分的三维各向异性介质中的弹性波。殷学鑫等 （201