利用无序度地震属性预测碳酸盐岩油气藏优质储层的方法

摘要

本发明实施例提供了一种利用无序度地震属性预测碳酸盐岩油气藏优质储层的方法。所述方法包括：在地表激发人工震源，利用地表的仪器(检波器)获取地震记录，对地震资料进行处理，得到三维叠后地震数据，利用测井资料进行井震标定，确定目的层段，在目的层段进行地震数据层位解释，提取无序度地震属性，利用无序度地震属性指导储层预测。与常规地震属性储层预测方法相比，具有物理意义明确，预测结果可靠等特点。

主权项

1.一种利用无序度地震属性预测碳酸盐岩油气藏优质储层的方法，其特征在于，所述方法包括：步骤1：在地表激发人工震源，利用检波器获取地震记录，对地震资料进行处理，得到三维叠后地震数据，利用测井资料进行井震标定，确定目的层段，在目的层段进行地震数据层位解释；其中：所述的三维叠后地震数据表示为：{D(i,x,n)|i_min<i<i_max,x_min<x<x_max,n_min<n<n_max}，其中，i代表主测线方向的坐标，x代表联络测线方向的坐标，n代表时间轴方向的坐标；步骤2：对于三维叠后地震数据D(i,x,n)的每一道地震数据求取复地震道，得到三维复地震道数据体Z(i,x,n)；其中：所述的每一道地震数据是指：在D(i,x,n)中取i和x为定值i₀和x₀时，由D(i₀,x₀,n)所表示的一维数据；所述的求取复地震道的方法包括：对于一个一维数据x(n)，它的复地震道为c(n)=x(n)+j·y(n)，其中j为虚数单位，y(n)=IDFT[-j·(Z(k)-X(k))]，$Z\left(k\right)=\left\{\begin{array}{cc}X\left(k\right)& k=0\\ 2X\left(k\right)& k=\mathrm{1,2},...,\frac{N}{2}\\ 0& k=\frac{N}{2},...,N-1\end{array}-1,\right.$X(k)=DFT[x(n)]，k=0,1,…,N-1DFT[]代表正离散傅里叶变换算子，IDFT[]代表反离散傅里叶变换算子；步骤3：将三维复地震道数据体Z(i,x,n)的每个数据点作为z₀，按迭代公式z_n+1=z_n-f(z_n)/f′(z_n),n=0,1,2…迭代30次，将最终迭代结果记为Z_iter(i,x,n)；其中，f(z_n)表达式为f(z_n)=z^m-ρ^m，其中m为整数，一般取m=12、13、…、18，ρ=max(D(i,x,n))/100；步骤4：对Z_iter(i,x,n)的每个数据点求取幅角，记为A(i,x,n)；其中，求取幅角的方法包括：对于任意复数c=x+j·y，其中j为虚数单位，其幅角定义为a=arctan(y/x)；步骤5：选取计算无序度滑动窗体，窗体大小为N_i×N_x×N_n，将该窗体在三维数据体A(i,x,n)上沿各个坐标轴方向逐点滑动，截取小数据体$A_{\mathrm{sub}}(i,x,n)=\left\{A(i,x,n)\right||i-\frac{N_i}{2}<i<i+\frac{N_i}{2},x-\frac{N_x}{2}<x<x+\frac{N_x}{2},n-\frac{N_n}{2}<n<n+\frac{N_n}{2}\},$计算小数据体的无序度，并将计算结果作为小数据体中心点的无序度，当滑动窗体遍历整个三维数据体A(i,x,n)时，最终得到无序度地震属性体I(i,x,n)；其中，滑动窗体大小通常取为N_i=11，N_x=11，N_n=11；其中，小数据体A_sub(i,x,n)的无序度的计算方法为：$I_{\alpha }(P_0,P_1,...,P_{m-1})=\left\{\begin{array}{cc}{(2^{1-\alpha }-1)}^{-1}[\underset{s=0}{\overset{m-1}{\Sigma }}{P_s}^{\alpha }-1]& ,\alpha \ne 1\\ -\underset{s=0}{\overset{m-1}{\Sigma }}{P}_s{\log }_2{P}_s& ,\alpha =1\end{array}\right.$其中K(l)为在A_sub(i,x,n)中数值等于l的数据点的个数，并通常取α的值为α=1或α=2；步骤6：在目的层段，利用解释好的地震层位开时窗，统计时窗内无序度地震属性体I(i,x,n)沿时间轴方向的累计求和值，得到平面储层预测图。