[发明专利]基于剪切破裂共轭准则的叠后裂缝储层表征方法

权利要求书

1.基于剪切破裂共轭准则的叠后裂缝储层表征方法，其特征在于，包括以下步骤：

步骤一、通过识别目标数据体的断层目标特征，分析沿地震轴变化趋势方向/倾向的相似特征属性，获得成果断层保护数据体；

步骤二、针对步骤一的成果断层保护数据体，进行脊线断裂增强处理，通过寻找沿地震轴倾向方向上的产状突变特征，获得断裂增强成果数据体；

步骤三、在步骤二的断裂增强成果数据体之上，按照蚁群种子点追踪的方法，识别微小断裂；

步骤四、根据微小断裂微观“共轭”的线性特征，结合其成带分布统计规律，对步骤三成果进行更小的裂缝特征的进一步强化，根据各功能需求进行差异特征的提取/共性特征的增强，依据剪切破裂共轭准则进行X型裂缝储层预测；

所述步骤四中根据各功能需求进行差异特征的提取/共性特征的增强，依据剪切破裂共轭准则进行X型裂缝储层预测过程具体包括以下子步骤：

S401，以目标模拟点张量场为基础，计算邻域内各方向的权值核函数K，使用最优化聚类方法加权拟合模拟点主要方向的裂缝属性值，如下式所示：

式中p_z表示数据窗口内第z点位置；p₀表示窗口中心位置；h_p表示核函数K_p的方差；表示方向模型m的z点与中心点的距离；h_m表示方向模型m的方差；K_m(M_z,M₀)表示方向模型m的z点位置处的核函数值；

S402，建立张量梯度方向的和，进行种子点路径拟合，如下式所示：

式中STGO_(i,j)表示窗口w′内(i，j)点的张量梯度方向的和；G_w′(m,n)表示窗口w′内(m，n)点处的张量梯度方向；G_w′(i,j)表示使用高斯核函数滤波后的梯度；表示窗口w′内(i，j)点的梯度，该梯度使用高斯核函数平滑后计算得到，以降低噪音影响；T表示梯度的判断阈值，绝对值大于这个值作为有效方向，小于此值，方向不加入统计；

S403，合成正交张量梯度方向和，该过程如下式所示：

ISTGO_h(i,j)＝STGO_ru(i,j)+STGO_rd(i,j)

ISTGO_v(i,j)＝STGO_ld(i,j)+STGO_rd(i,j)

S_I＝max{|ISTGO₁|...|ISTGO_i″|...|ISTGO₈|}

式中ru表示右上方向，rd表示右下方向，右上和右下可以合成水平方向的梯度标识ISTGO_h；ld表示左下方向，左下和右下可以合成垂直方向的梯度标识ISTGO_v；S_I表示主方向标识；ISTGO_i″表示第i″个方向上的合成正交张量梯度，取梯度最大的方向作为主方向；

S404，进行张量梯度方向的判定，基于裂缝模拟结果，统计单位半径R内的相对裂缝密度与裂缝发育主方位。

2.基于权利要求1所述的剪切破裂共轭准则的叠后裂缝储层表征方法，其特征在于，所述步骤一中获得成果断层保护数据体过程具体包括；通过分析目标位置线、道测线位置的不同地震同向轴的产状和趋势，保留现有的构造特征及趋势，避免裂缝特征被滤波所模糊化，从而得到一个去除噪声背景并且保留原始地震产状及构造特征的成果断层保护数据体。

3.基于权利要求1所述的剪切破裂共轭准则的叠后裂缝储层表征方法，其特征在于，所述步骤二中针对步骤一的成果断层保护数据体，进行脊线断裂增强处理过程包括：

S201，通过基于种子区域增长的脊线方向估计，该过程如下式所示：

式中RB表示待计算脊线方向的区域，RB_adj表示该区域的四/八联通区域；θ表示搜索区域内用正弦波建模拟合的方向角度；f表示搜索区域内正弦波建模的频率；A表示搜索区域内所有样点的振幅值；w表示该区域内的高斯拟合函数，用于加权平滑区域内的振幅以降低区域内噪音水平；x和y表示区域RB的中心位置坐标；当两个区域连通 时，O_f值为1，不连通时值为0；

S202，通过非线性脊线加强方法，进行脊线的相干增强处理，该过程如下式所示：

式中ρ是相干增强因子，取值范围在0～1之间；AA'和BB'分别表示(i′，j′)位置处的断层主方向及其法向；g(n_AA'-k)表示AA'方向上，以第k个CDP点为中心的高斯函数；表示第k个CDP点与(i′，j′)位置点相干值的梯度；ΔC_BB'表示BB'方向上的相干梯度；C_i′j′是(i′，j′)的CDP位置处的相干属性值。