[发明专利]一种基于等级相关分析的多分量转换波裂缝预测方法

摘要

本发明提供了一种基于等级相关分析的多分量转换波裂缝预测方法，所述预测方法包括以下步骤：选取计算时窗，将待处理采样点的径向分量和横向分量按照计算时窗进行截取；对待处理采样点进行等级相关分析，得到最大等级相关系数；判断待处理采样点处是否具有裂缝发育，确定裂缝发育方位；对判定有裂缝发育的待处理采样点径向分量和横向分量进行分离，得到快横波和慢横波；根据快横波和慢横波时差计算各向异性系数；重复所述选取计算时窗步骤至所述计算各向异性系数步骤，循环至下一个待处理采样点直到地震数据计算完毕。根据本发明的方法采用了判别方法对计算结果进行判别，剔除裂缝明显不发育区域，在对裂缝发育进行统计时更为准确。

主权项

1.一种基于等级相关分析的多分量转换波裂缝预测方法，其特征在于，所述预测方法包括以下步骤：/n选取一个包含目标深度范围的计算时窗，将待处理采样点的径向分量地震数据和横向分量地震数据按照计算时窗进行截取；/n对待处理采样点进行等级相关分析，得到等级相关系数最大值；/n判断待处理采样点处是否具有裂缝发育，若有裂缝发育，确定裂缝发育方位；/n对判定有裂缝发育的待处理采样点进行快横波和慢横波分离；/n根据快横波和慢横波时差计算各向异性系数；/n重复所述选取计算时窗的步骤至所述计算各向异性系数的步骤，循环至下一个待处理采样点直到地震数据计算完毕，其中，/n所述对待处理采样点进行等级相关分析的步骤包括：/n将截取的待处理采样点的径向分量地震数据和横向分量地震数据分别记录到数组R和数组T中；/n对数组R和数组T进行旋转角度θ_j和时移τ_j的变换，得到变换后的分量LR(θ_j,t)和LT(θ_j,t+τ_j)，其中，旋转角度变换式为：/n /nLT(θ_j,t+τ_j)＝LT(θ_j,t)*τ_j，其中，θ_j＝θ_k+Δθ，τ_j＝τ_k+Δτ，θ_k为起始旋转角度，τ_k为起始时移，Δθ为旋转角度循环步长，Δτ为时移循环步长，t为时间；/n分别求取LR(θ_j,t)分量和LT(θ_j,t+τ_j)分量的秩，记为ZR_i和ZT_i；/n利用式1计算等级相关系数，得到等级相关系数最大值r_max＝max{r(θ_j,τ_j)}，其中，式1为/n 其中，D_i＝ZR_i-ZT_i，n为分量中数据个数；/n所述判断待处理采样点处是否具有裂缝发育的步骤包括：/n将待处理采样点的最大值r_max与判别值T₁进行比较，若r_max＞T₁，则计算所述对待处理采样点进行等级相关分析步骤中的待处理采样点所有r(θ_j,τ_j)值的标准差D_r；若D_r＞T₂，则判定该待处理采样点处有裂缝发育，并记录最大值r_max下对应的旋转角度和时移，其中，判别值T₁，T₂为经验值或给定值。/n