[发明专利]一种基于倾角成像增强的火山岩裂缝预测方法

权利要求书

1.一种基于倾角成像增强的火山岩裂缝预测方法，其特征在于，包括步骤：

S1倾角成像增强：对原始地震资料进行多尺度多方位倾角成像增强，利用倾角及方位角的变化计算相邻道的相似性，以提高地震资料横向信噪比；

S2断层加强处理：对原始地震资料进行断层加强处理，以提高地震资料垂向分辨率，提高断层可识别度；

S3基于张量方向场约束的裂缝预测：在倾角成像增强和断层加强处理后的地震资料基础上，基于张量场方向约束加权拟合进行火山岩裂缝预测，预测出裂缝发育强度，结合区域地质分析最终确定火山岩储层分布区。

2.如权利要求1所述的一种基于倾角成像增强的火山岩裂缝预测方法，其特征在于，所述步骤S1包括：

S11：根据地震三维数据体，定义地震数据的分方位分尺度特征为ψ_u,v(x)，

式(1)中，i为复数算子，σ为滤波器带宽，k_u,v＝k_v(cosθ,sinθ)^T，k_v＝2^(2-v)/2π，θ＝u·π/K，v对应于滤波器的尺度，u对应滤波器的方向，||*||表示模；

对于给定的三维地震数据I(x),对于给定的任意一个空间点x＝(in,cl,t)，定义不同方位不同尺度的地震数据为D_u,v(x)：

D_u,v(x)＝ψ_u,v(x)*I(x) (2),

式(2)中，D_u,v(x)为根据不同的(u,v)决定的滤波器得到的不同方位不同尺度的地震数据；

采用3个尺度和8个方向对地震数据进行联合分析，得到D_u,v(x)的24中表达形式：

式(3)中，变量m和n为下标索引变量；

S12：对D_u,v(x)进行编码，

式(4)、(5)中，和为D_u,v(x)在(u,v)参数下多尺度多方位滤波数据的二进制编码，并计算合成数据和

式(6)、(7)中，和为通过不同尺度和不同方位的滤波处理；

S13：根据分解数据的相关性，定义相关矩阵为：

将用特征值和特征向量表示为：

式(10)中，v是尺度下标，为的三个非负特征值，且满足分别为对应的特征向量，构成一个局部的正交坐标系，为局部领域内对比度最大的方向，成一个垂直于的局部平面，通过求取倾角属性：

将用特征值和特征向量表示为：

式(12)中,v是尺度下标，为的三个非负特征值，且满足分别为对应的特征向量，构成一个局部的正交坐标系，为局部领域内对比度最大的方向，成一个垂直于的局部平面，通过求取倾角属性：

最佳倾角dip(x)为：