专利名称：利用无序度地震属性预测碳酸盐岩油气藏优质储层的方法
技术领域：
本发明涉及地球物理勘探技术领域，特别是涉及一种利用无序度地震属性预测碳酸盐岩油气藏优质储层的方法。
背景技术：
地震勘探是在地表激发人工震源，由震源所引起的震动以地震波的形式向地下传播，并在一定的条件下向上反射传回地表，然后由地表的仪器(检波器)记录反射回来的地震波，从而得到地震记录(也叫地震资料)。地震记录是经由地下介质而得到的，必然受到地下介质物理性质(如岩性、孔隙度、密度以及流体性质)的影响，这种影响继而被作为地下介质的综合信息保存在地震记录中，有针对性的对地震资料进行处理和解释便可以得到有用的地下信息，进而指导地质目标勘探。 地震属性是地震资料处理与解释的一个重要分支，是地震资料及其衍生资料的几何学、运动学、动力学及统计学特征的一种量度，通常是由叠前或者叠后地震数据经过各种数学变换得到的。地震属性最早出现于20世纪60年代末，以当时直接用于烃类检测的亮点技术为标志，随后，涌现出大量的属性定量提取方法，例如，70年代末到80年代末的复地震道分析、自相关分析、频谱分析等，90年代以后，又出现了以相干体、倾角、方位角等为代表的一批属性，它们在油气勘探与开发中发挥着越来越重要的作用。碳酸盐油气藏通常具有控制因素复杂、规律性差、非均质性强等特点，其储层预测工作一直是地球物理勘探的一个难点。目前，尽管地震属性种类多样，但由于它们的物理意义和适用条件不明确，储层预测结果的可靠性难以评价，很少能够直接用于指导碳酸盐岩油气藏优质储层预测。

发明内容
碳酸盐岩储层质量通常与其孔洞缝发育程度和空间结构密切相关。当地震波传播到碳酸盐岩储层发育区时，由于该区孔缝洞结构相对发育，导致地震波反射、绕射作用复杂多样，并最终表现为地震数据体中的某种无序性。因此，定量描述地震响应的无序程度是进行碳酸盐岩储层识别和物性预测一条新途径。本发明的目的是根据碳酸盐岩储层的上述特点，提出了一种快速、便捷的利用无序度地震属性预测碳酸盐岩油气藏优质储层的方法。本发明实施例提供了一种利用无序度地震属性预测碳酸盐岩油气藏优质储层的方法，所述方法包括I)设计地震记录采集方式(也称定义观测系统)，在地表激发人工震源，利用地表的仪器(检波器)获取地震记录(也称地震资料);对地震资料进行处理，得到三维叠后地震数据；利用测井资料进行井震标定，确定目的层段；在目的层段进行地震数据层位解释。步骤I :所述的三维叠后地震数据可以表示为{D(i, X，n) I imin〈i〈imax, xmin〈x〈xmax, nmin〈n〈nmax}，其中，i 代表主测线方向的坐标，x代表联络测线方向的坐标，n代表时间轴方向的坐标。步骤2 :对于三维叠后地震数据D(i，x，n)的每一道地震数据求取复地震道，得到三维复地震道数据体Z (i，x，n)。其中所述的每一道地震数据是指在D(i，x, n)中取i和X为定值Itl和Xtl时，由D(i0, x0, n)所表示的一维数据。所述的求取复地震道的方法是对于一个一维数据x(n)，它的复地震道为c (n) =x (n) +j y (n)，其中 j 为虚数单位，y (n) =IDFT [-j (Z (k) -X (k))],
权利要求
1.一种利用无序度地震属性预测碳酸盐岩油气藏优质储层的方法，其特征在于，所述方法包括 步骤I:在地表激发人工震源，利用地表的仪器(检波器)获取地震记录，对地震资料进行处理，得到三维叠后地震数据，利用测井资料进行井震标定，确定目的层段，在目的层段进行地震数据层位解释； 其中所述的三维叠后地震数据表示为{D(i, X，n) I imin〈i〈imax, xmin〈x〈xmax, nmin〈n〈nmax},其中，i 代表主测线方向的坐标，x 代表联络测线方向的坐标，n代表时间轴方向的坐标； 步骤2 :对于三维叠后地震数据D (i, X, n)的每一道地震数据求取复地震道，得到三维复地震道数据体Z(i，x，n)； 其中所述的每一道地震数据是指在D(i，x，n)中取i和X为定值Itl和Xtl时，由D(i0, x0, n)所表示的一维数据； 所述的求取复地震道的方法包括对于一个一维数据x(n)，它的复地震道为c (n) =x (n) +j y (n)，其中 j 为虚数单位，y (n) =IDFT [-j (Z (k) -X (k))], X{k)k =O Z(k) = <2X(k)灸=1，2，….警-I N 、 2 X (k) =DFT [x (n) ]，k=0, I,…，N-I DFT[]代表正离散傅里叶变换算子，IDFT[]代表反离散傅里叶变换算子。
2.如权利要求I所述的利用无序度地震属性预测碳酸盐岩油气藏优质储层的方法，其特征在于，所述方法还包括 步骤3:将三维复地震道数据体Z(i，x，n)的每个数据点作为Z(l，按迭代公式zn+i=zn_f (zn)/f' (Zn)，n=0, 1，2…迭代 30 次,将最终迭代结果记为 Ziter (i，x，n)； 其中，f(zn)表达式为f (Zn)=Zm-Pm,其中m为整数，一般取m=12、13、…、18，P =max (D (i, x, n)) /100。
3.如权利要求I所述的利用无序度地震属性预测碳酸盐岩油气藏优质储层的方法，其特征在于，所述方法还包括 步骤4 :对ZitOT(i, X，n)的每个数据点求取幅角，记为A (i, X，n)； 其中，求取幅角的方法包括对于任意复数c=x+i y，其幅角定义为a=arctan(y/x)。
4.如权利要求I所述的利用无序度地震属性预测碳酸盐岩油气藏优质储层的方法，其特征在于，所述方法还包括 步骤5 :选取计算无序度滑动窗体，窗体大小为NiXNxXNn,将该窗体在三维数据体A(i,x, n)上沿各个坐标轴方向逐点滑动，截取小数据体 fNN NN NN I Aub (h X, n) = \A{i,x,}i)+ .A- — ~~^-<n<n + ^- \，计算小数据体的无序度，并将计算结果作为小数据体中心点的无序度，当滑动窗体遍历整个三维数据体A(i，x，n)时，最终得到无序度地震属性体I (i，x，n)； 其中，滑动窗体大小通常取为Ni=Il, Nx=Il, Nn=Il ； 其中，小数据体Asub(i，x，n)的无序度的计算方法为
5.如权利要求I所述的利用无序度地震属性预测碳酸盐岩油气藏优质储层的方法，其特征在于，所述方法还包括 步骤6 :在目的层段，利用解释好的地震层位开时窗，统计时窗内无序度地震属性体I(i,x, n)沿时间轴方向的累计求和值，得到平面储层预测图。
全文摘要
本发明实施例提供了一种利用无序度地震属性预测碳酸盐岩油气藏优质储层的方法。所述方法包括在地表激发人工震源，利用地表的仪器(检波器)获取地震记录，对地震资料进行处理，得到三维叠后地震数据，利用测井资料进行井震标定，确定目的层段，在目的层段进行地震数据层位解释，提取无序度地震属性，利用无序度地震属性指导储层预测。与常规地震属性储层预测方法相比，具有物理意义明确，预测结果可靠等特点。
文档编号G01V1/28GK102707315SQ20121018337
公开日2012年10月3日 申请日期2012年6月5日 优先权日2012年6月5日
发明者张研, 李劲松, 杨昊, 郑晓东 申请人:中国石油天然气股份有限公司