[发明专利]地下裂缝预测方法

权利要求书

1.一种地下裂缝预测方法，该方法包括：

获得横波的径向分量R(t)和横向分量T(t)，其中t为获得地震数据的时间；

根据所述径向分量R(t)和横向分量T(t)求取快横波分量S₁(t)和慢横波分量S₂(t)之间的时差函数Δt(t)；

对所述快横波分量S₁(t)和慢横波分量S₂(t)之间的时差函数Δt(t)求导，得到该时差函数的导数函数dt(t)；以及

利用所述导数函数dt(t)对裂缝的发育强度进行预测。

2.根据权利要求1所述的方法，其中，在所述获得横波的径向分量R(t)和横向分量T(t)的步骤中，由地震勘探中的数字三分量检波器的两个水平分量接收到的地震数据得到所述横波的径向分量R(t)和横向分量T(t)。

3.根据权利要求1所述的方法，其中，在求取快横波分量S₁(t)和慢横波分量S₂(t)之间的时差函数Δt(t)的步骤中，通过以下步骤来求取所述时差函数Δt(t)：

a)根据下列方程计算所述快横波分量S₁(t)和慢横波分量S₂(t)之间的剩余误差E(θ，Δt，t)：

E(θ，Δt，t)＝(R(t)cosθ+T(t)sinθ)-(R(t-Δt)sinθ-T(t-Δt)cosθ)

-(N_R(t)cosθ+N_T(t)sinθ)-(N_R(t-Δt)sinθ-N_T(t-Δt)cosθ)

其中，θ为裂缝的方位角，Δt为所述快横波分量S₁(t)和所述慢横波分量S₂(t)之间的时差，N_R(t)和N_T(t)分别为沿径向和横向的噪声；

b)对于一时窗内的n个采样点，对所述快横波S₁(t)和慢横波S₂(t)之间的剩余误差E(θ，Δt，t)进行累加，得到如下的累加函数F(θ，Δt)，

F(θ,Δt)=(Σk=1nE(θ,Δt,tk)p)1p]]>

其中，t_k为第k个采样点的时间，p为累加指数，p＝2、3、4或5；

c)对方位角θ与时差Δt的值进行扫描获得的二维曲面，并得到使该二维曲面达到最小值点的极值方位角θ_min与极值时差Δt_min，该极值方位角θ_min与极值时差Δt_min作为所述时窗的时间中点处的裂缝方位角θ与时差Δt；

d)在时间轴上连续移动所述时窗，得到以时间t为变量的时差函数Δt(t)。

4.根据权利要求3所述的方法，其中，对时差函数求导的步骤还包括对得到的导数函数dt(t)进行下列处理中的至少一项处理：去除异常值处理、插值处理、以及平滑处理。

5.根据权利要求3所述的方法，其中所述采样点的个数n大于1，所述时窗的大小为采样点的个数n乘以采样间隔；

6.根据权利要求5所述的方法，其中所述采样间隔是相邻的两个采样点间的时间间隔，该时间间隔为0.1ms-20ms。