[发明专利]一种高精度地震波层速度反演方法

权利要求书

1.一种高精度地震波层速度反演方法，特点是采用以下步骤：

1)地面震源激发，井下检波器接收得到VSP数据；

2)对地震记录进行傅里叶变换得到振幅谱，将高频部分的振幅谱清零，然后反傅里叶变换，得到只保留低频信号的地震记录，然后拾取每道地震记录的初至；

3)对步骤2)中的地震记录进行傅里叶变换得到振幅谱，然后在深度方向再进行一次傅里叶变换得到F-K谱，并将上行波对应的F-K谱清零，然后分别在波数方向和频率方向进行反傅里叶变换得到下行波；

4)从每道的下行波中沿初至往后开一时窗，对该时窗内的信号进行傅里叶变换，得到每个频率的相位谱；

5)重复步骤4)，得到所有道下行波每个频率的相位谱；

6)求取相邻道中相同频率的相位谱之差；

7)对步骤6)中的相位谱之差进行最小二乘拟合，得到相邻道初至差；

8)相邻道的深度差除以相邻道初至差得到高精度的层速度。

2.根据权利要求1所述的方法，特点是步骤6)所述的相邻道相位谱之差计算如下式：

Si+1(ω)=∫-∞+∞si(t-Δti+1)e-j2πftdt=∫-∞+∞si(t)e-j2πf(t-Δti+1)dt]]>

则有：θ_i+1(ω)＝2πfΔt_i+1

式中，s_i(t-Δt_i+1)为深度为h_i的检波器得到的下行波为s_i(t)，深度为h_i+1的检波器得到的下行波；

S_i+1(ω)为s_i+1(t)的傅里叶变换；

θ_i+1(ω)为s_i+1(t)的傅里叶变换的相位谱

Δt_i+1为第i+1道和第i道初至差；

j为复数单位。