[发明专利]一种倾斜层状粘弹性介质中波场正演模拟方法

说明书

技术领域

本发明属于地震勘探技术领域，涉及一种模拟方法，尤其是一种倾斜层状粘弹性介质中波场正演模拟方法，其是一种针对含倾斜界面的层状粘弹性介质中地震波场高效精确的正演模拟技术。

背景技术

波场模拟在油气勘探中占有重要地位，尽管已存在多种类型的波场模拟方法，然而，针对要解决的问题的要求，发展波场模拟理论及方法一直是很活跃的研究领域。另外，在油气勘探等领域，为了解释地震图，或进行数据反演，常常需要仅模拟指定类型的波，譬如仅模拟一次反射波，或仅模拟转换波等。常用的波场模拟方法有多种类型，如基于射线理论的方法，数值解方法，广义反射率法以及其推广而得到的反/透射矩阵法等。这些方法各有其优缺点。

射线方法能适用于较复杂的介质模型，可根据事先给定的射线码模拟感兴趣的类型的波，且计算速度快；不足之处是，它只是波动方程的高频近似，用于模拟复杂储层的地震响应具有局限性。采用数值法求解波动方程，是复杂介质中波场模拟的有效途径。常用的数值方法有：有限差分法，有限元法，伪谱法和谱元法等，这类方法可模拟地震波全波场，精度高。但计算量巨大，不能像射线类方法，有选择性的仅模拟我们感性趣的类型的波。

层状介质中波传播理论及波场模拟技术一直受到人们的重视。Aki和Müller等人发展的反射率法，能快速精确地模拟水平层状粘弹性介质的全波场或指定类型的波。Berkhout提出了波场外推法可用于带倾斜界面的弹性介质波场模拟。Chen等人发展的全局广义反/透射矩阵法和边界元法，能模拟带不规则界面的层状弹性介质中全波场。但这些方法不适用于含倾斜界面的层状粘弹性介质。

发明内容

针对上述问题，本发明的目的在于针对含倾斜波阻抗界面的层状粘弹性介质，提出一种倾斜层状粘弹性介质中波场正演模拟方法，该方法可以选择性模拟感兴趣的各种类型的波，如一次反射P波、S波和转换波，且模拟的波场非常精确。

本发明的目的是通过以下技术方案来解决的：

该种倾斜层状粘弹性介质中波场正演模拟方法，包括以下步骤：

1）将点源激励的球面波分解为平面波；

2）将平面波在倾斜层状粘弹性介质中传播，遇到界面时进行如下处理：以该界面为轴，界面与全局坐标系Z轴的交点为原点建立局部坐标系，垂直于界面向下，将全局坐标系中的入射波转换到局部坐标系中；在局部坐标系下求得反射波和透射波；再将局部坐标系下的反射波和透射波转换到全局坐标系下，得到全局坐标系下的反射系数和透射系数，以及反射波和透射波；在每个界面上根据需要选取波的类型，P波入射时能够选取的波包括反射PP波、反射PS波、透射PP波、透射PS波；S波入射时能够选取的波包括反射SP波、反射SS波、透射SP波、透射SS波；从而模拟所选取的波类型，得到模拟的平面波；

3）由步骤2）模拟得到的平面波根据下式合成点源激发的波场：