复杂介质地震波场模拟的单程波方法.pdfVIP

第 48卷第 5期 石 油 物 探 Vo1．48，No．5 2009年 9月 GEoPHYSICAL PRoSPECTING F0R PETROLEUM Sep．，2009 文章编号：1000—1441(2009)05 0479—04 复杂介质地震波场模拟的单程波方法 首 皓 ，孙鲁平 ，刘红伟 ，赵晓龙。 (1．中国石油勘探开发研究院，北京 100083；2．中国科学院地质与地球物理研究所，北京 100029； 3．中国石油天然气集团公司川I庆钻探工程有限公司工程技术研究院，陕西西安 710018) 摘要：单程波方法是常用的地震正演模拟方法之一。与双程波方法相比，单程波方法的计算量成倍减少，且能保 持较高的精度。但由于单程波方法不能获得准确的入射角信息 ，因而多采用垂直入射的反射系数，使得模拟波 场的振幅特征不是实际地层的真实反映。为此，提出了结合频率一波数域反射系数近似公式和相移加插值的改 进单程波方法。分别采用双程波虚谱法、常规单程波方法和改进单程波方法，基于Marmousi模型进行了数值模 拟，对比模拟结果发现，改进的单程波方法可以准确地模拟振幅随入射角的变化，且具备与双程波方法相同的 精度 。 关键词：地震波正演模拟；波动方程；单程波 ；相移加插值；反射系数随入射角变化 DOI：10．3969／j．issn．1000—1441．2009．05．009 中图分类号：P631．4 文献标识码 ：A 地震波正演模拟不仅是研究地震波传播规律 值模拟 中。 的重要手段，也是地震资料偏移成像的基础。地震 利用单程波方法进行波场延拓时无法获得波 波正演模拟通常通过求解波动方程来获得走时和 场的入射角，因此不能利用随入射角变化的反射系 振幅信息，但在横向非均匀介质情况下难 以获得波 数近似公式获得准确的振幅特征 。我们将 Gazdag 场准确的解析解 ，由此产生了高频近似射线法、双 等_l](1984)提 出 的相移 加插值 法和 deBruin 程波方程解法和单程波方程解法等，以使其适应非 等E](1990)提出的频率一波数域反射系数近似公 均匀介质。对 于高频近似解法，Cerveny等口] 式与常规单程波正演方法相结合，相移加插值法可 (1977)提出了射线渐进近似解法；Chanpman等 ¨2] 以对复杂介质中的地震波进行准确模拟，频率一波 (1982)提出了Maslov方法；CervenyE。](1983)提出 数域反射系数包含了角度信息，因此能够获得准确 了高斯射线束方法。基于射线理论的方法一般仅 的振幅特征 。基于 Marmousi模型和双程波方法 适用于弱非均匀介质条件 。在双程波方程解法中， 模拟的地震记 录，对本文方法 的正确性进行 了 BooreE4](1972)和 KellyEI_(1976)：~]1用有 限差分法 、 验证。 SmithE](1974)利用有限元法、Wong等l7](1975) 利用边界积分法、Orzag[8](1980)利用虚谱法直接 1 方法原理 求解双程波方程 。双程波方法适用于横 向非均匀 介质，但这类方法计算时间长 ，在实际中难 以得到 1．1 单程波方程 广泛应用 。 只考虑二维纵波情况，地震波的传播可用以下 单程波方法是将双程波方程分解为上行波波 声波方程描述 ： 广 1 ] 场和下行波波场 ，通过反射系数来实现上行波波场 01(，) ‘1． P(z，z，J+ 和下行波波场的转换 。单程波算法可以正确模拟