分裂步傅里叶法地震波场数值模拟.pdfVIP

2008年第2期 内蒙古石油化工 l25 分裂步傅里叶法地震波场数值模拟 熊 忠 (成都理工大学信息工程学院，四川成都 610059) 摘 要 ：通过用分裂步傅 里叶方法在频率波数域作波场延拓 ，实现 了二维复杂地质模型的地震波场 数值模拟 ，得到 了合理的正演地震记录和偏移图像。结果表明该方法计算效率高，精度高，稳定性好 ，不 受倾 角限制 ，能适应速度纵横 向任意变化 ，是一种比较好的地震波场模拟算法。 关键词：分裂步傅里叶法；波场延拓；数值模拟 1 引言 △0P+∞ s5P一一∞ (2s0△s+△s0)P (3) 地震波场数值模拟无论是在天然地震监测还是 略去扰动的二阶项△s ，得： 在人工地震勘探工程中都十分重要，一直以来科学 △ P+∞ s5P一一2∞。s0△sP (4) 工作者们对其研究热情始终不减，发展出了很多方 对非齐次方程(4)，令为对的傅氏变换，使用参 法。波场模拟包括正演和偏移两个方面，正演是假设 考慢度作相移延拓 地下地质情况已知，如地质体的空间分布、介质的速 P1(k ，k ，Zi+△z，∞)==P(k ，k ，zi，∞)eikzo (5) 度、密度等，应用地震波运动学和动力学的基本原 其中：K o—sqrt(∞ s5一k：一k；) (6) 理，计算出所给地质模型的地震响应。而在速度模型 在频率一空间域中，对慢度的扰动△s(x，y，z)， 已知的情况下，偏移可视为地震波在地下传播过程 使用第二次相移延拓得到： 的计算机逆模拟。采用波动方程理论模拟地震波的 P(x，y，zi+△z，∞)一e’ i P1(x，y，zI+△z，∞) 传播，不仅能保持地震波的运动学特征，而且还能保 (7) 持地震波的动力学特征。求解波动方程的方法有以 其中P1(x，y，z +△z，∞)是 由P；(k ，ky，zi+△z，∞) 下三种：有限差分法、Kirchhoff积分法和频率一波 对k ，k 作反傅氏变换得到的。 数域法。用于地质模型正演与偏移的波场数值模拟 对有限频带积分有： 方法需满足保振幅、高精度、高稳定性、能适应大倾 1 f 一 角和纵横向速度剧烈变化等特殊要求，而频率波数 P(x，y，zi+△z)一 J P(x，y，zi+△z，∞)d∞ (8) 域方法因其快速简洁、不受倾角限制、精度高、稳定 式(5)和式(7)即分裂步傅里叶法波场延拓的基 性好而特别受到人们的青睐，因此本文采用 P．L． 本公式，其中式(5)在频率一波数域实现，式(7)在频 2 Stoffa提出的分裂步傅里叶法频率波数域波场延拓 ～ 算子对地震波场进行延拓，并进一步利用爆炸反射 率一空间域实现。在计算中要注意，(k：+k；) 的绝 原理对其进行正演和偏移。 对值要小于等于1，若不成立就令k 一0。正演和偏移 2 分裂步傅里叶法波场模拟的原理 相互之间是共轭的，计算时方向和符号是