[发明专利]一种台盆区走滑断裂空腔储层刻画方法

权利要求书

1.一种台盆区走滑断裂空腔储层刻画方法，其中，所述台盆区走滑断裂空腔储层刻画方法包括以下步骤：

S100、获取台盆区走滑断裂空腔储层的钻井资料、测井资料和地震资料；

S200、针对断裂空腔储层内部的孔洞型储层进行波阻抗反演，采用基于子波分解的降频速度反演方法，刻画串珠空间分布；

S300、计算三维地震数据体的梯度构形属性，以识别缝洞异常体；

S400、使用已钻井标定构形属性门槛值范围，以刻画断裂空腔储层的断裂带边界；

S500、计算三维地震的不连续属性，以刻画断裂空腔储层内部的裂缝结构；

S600、在孔洞储层的波阻抗反演基础上进行地质学统计趋势反演，以提高对小缝洞储层的识别和刻画能力，结合S500，完成对断裂空腔储层内部的刻画；

S300中所述计算梯度构形属性的过程包括：

S301、建立梯度构形属性：

计算三维地震数据体每一点梯度矢量：

式中x、y、z分别为三维地震体素的线、道位置及双程旅行时间方向；g为梯度矢量；g₁、g₂、g₃分别为沿x、y、z方向的方向导数；u为在对应方向的导数矩阵；

利用三维高斯窗对三维数据体梯度进行平滑：

式中G为高斯核函数；g为数据体计算出的振幅梯度场；g₁、g₂、g₃分别为沿x、y、z方向的方向导数；1、2、3分别指代x、y、z方向；σ_g为对应i方向的长度；

S302、计算梯度构形属性T：

S303、利用高斯窗对梯度构形属性的每一个成分进行平滑：

S304、计算梯度构形属性的特征值和特征向量：

其中λ和v分别是特征值和特征向量；λ₁≥λ₂≥λ₃＞0，三个特征向量v₁、v₂、v₃两两正交，1、2、3分别指代x、y、z方向；包括以下几种情况：

1)λ₁＝λ₂＝λ₃＞0，这种情况下是各向同性的；

2)λ₁＞λ₂＝λ₃＝0，这时的λ₁对应的矢量方向接近Z方向；

3)λ₁＞λ₂＞λ₃＝0，这种情况下根据特征值和特征向量准确计算滤波向量，同时对有效的断裂进行保存；

S500中计算地震不连续属性的步骤包括：

设时间窗口内由N个样点元素组成的一个分析道与以道为单位的窗口内的J道的数据排列在一起可以形成数据矩阵：

其中，u_ij表示第j道上第i个地震波的振幅值；矩阵U的n行是分析体内每个地震道上的第n个振幅值，则第n个采样的方差矩阵为：

从而相干分析时窗内道中各元素协方差构成的协方差矩阵为：

设协方差矩阵M有J个从大到小排序的特征值λ₁，λ₂，...，λ_J，对应的特征向量为v₁，v₂，...，v_J，则基于特征结构的相干体计算公式为：

式中，λ_j为计算时窗内地震数据协方差矩阵的特征值；

S600中所述地质学统计趋势反演的过程包括：

S601、建立趋势模型；

设某空间地理要素的实际观测数据为：Z_i(x_i|y_i)(i＝1，2，...，n)，趋势拟合值为则有式中：ε_i即为剩余值，即残差值；

从实际观测值出发推算趋势面或体，采用回归分析方法，使得残差平方和趋于最小，即：

根据观测值x_i，y_i，z_i(i＝1，2，...，n)确定多项式的系数a₀，a₁，...，a_p，使残差平方和最小；

S602、趋势面模型的适度检验，包括趋势面与实际面的拟合度系数R2检验和F检验；

所述地质学统计趋势反演中的趋势体生成方法包括：

对于绝对阻抗属性体进行大网格滤波，形成一个渐变的趋势阻抗趋势体，再将原绝对阻抗属性体与此渐变趋势体进行相减，即得到了一个剩余阻抗体将剩余组抗体以常规地震属性的形式表现；

研究区根据趋势面方法原理，用绝对阻抗与趋势体相减得到一个相对低的地质学统计趋势波阻抗反演体。