[发明专利]一种基于加权匹配追踪的三维VSP地震资料内插方法

说明书

本发明公开了一种基于加权匹配追踪的三维VSP地震资料内插方法，包括以下步骤：S1、对野外采集的原始地震数据分别进行静校正、反褶积、VSP动校正处理，消除不同偏移距的道间时差，得到内插前的三维VSP地震数据；S2、根据内插前的三维VSP地震数据的炮点坐标，得到三维VSP地震数据的空间采样位置；S3、用空白数据对内插前的三维VSP地震数据在空间方向上进行加密，得到空间采样规则的数据体；S4、将步骤S3得到的空间采样规则的数据体进行F‑K变换，得到F‑K谱；S5、将F‑K谱进行线性拉冬变换，并根据拉冬变换结果求取加权算子；S6、用步骤S5得到的加权算子处理F‑K谱得到新的F‑K谱；S7、对新的F‑K谱做F‑K反变换得到内插后的三维VSP地震数据。

技术领域

本发明涉及地质勘探领域，尤其涉及一种基于加权匹配追踪的三维VSP地震资料内插方法。

背景技术

目前，三维VSP技术是一项高精度地震物理技术，可以获得高精度的三维成像资料，以其具有更直接的油藏描述(吸收衰减、速度、各向异性)能力、更好的时深储层标定能力和高于地面地震的分辨能力，已广泛应用于精细刻划井周复杂的构造、储层和流体的空间展布，被誉为地面、井间地震、测井和地质相结合的“桥梁”。

由于VSP数据采集在空间上通常是不规则的，影响对后续的资料的处理和解释。因此，在偏移前进行内插处理，将其加密成空间采样规则的数据是必要的。常规的VSP数据内插方法容易产生空间假频，反而不利于后续的资料处理和解释；用于常规地震资料的内插方法，由于观测系统的差异，不能直接对VSP数据进行内插处理。

发明内容

本发明目的是针对现有技术的不足，提供一种基于加权匹配追踪的三维VSP地震资料内插方法。

为了实现上述目的，本发明的技术方案是：

用三维VSP地震资料在空间方向的稀疏性作为约束条件，在压缩感知的理论框架下将内插处理转换为去噪问题；然后，分析了变换域结果中噪音和有效信号的分布特征，并根据该分布特征构造加权算子；最后，利用加权算子设计加权匹配追踪算法求解该压缩感知问题，得到内插后的三维VSP地震数据。

具体包括如下步骤：

S1、对野外采集的原始地震数据分别进行静校正、反褶积、VSP动校正处理，消除不同偏移距的道间时差，得到内插前的三维VSP地震数据；

S2、根据内插前的三维VSP地震数据的炮点坐标，得到三维VSP地震数据的空间采样位置；

S3、用空白数据对内插前的三维VSP地震数据在空间方向上进行加密，得到空间采样规则的数据体；

S4、将步骤S3得到的空间采样规则的数据体进行F-K变换，得到F-K谱；

S5、将F-K谱进行线性拉冬变换，并根据拉冬变换结果求取加权算子；

S6、用步骤S5得到的加权算子处理F-K谱得到新的F-K谱；

S7、对新的F-K谱做F-K反变换得到内插后的三维VSP地震数据。

进一步的，所述步骤S1中野外采集的原始地震数据是地面放炮井中接收或井中放炮地面接收得到的地震数据。

进一步的，所述步骤S4中F-K变换公式为：

其中，t为时间，x是空间，d(t，x)是原始地震数据，f为频率，k为波数，是转换到频率空间域得到F-K谱。

进一步的，所述步骤S5中的线性拉冬变换公式为：

其中，为斜率，R为拉冬变换结果。

进一步的，所述步骤S5中加权算子的计算公式为：