专利名称：一种获得碳酸盐岩储层储集空间体积的方法
技术领域：
本发明涉及一种获得碳酸盐岩储层储集空间体积的方法，尤其是涉及获得碳酸盐岩储层的油气储量体积的方法。
背景技术：
近年来，碳酸盐岩逐渐成为油气勘探的重要新领域，国内碳酸盐岩油气田通常埋藏深，储层成因复杂，包括同生期和准同生期岩溶作用、表生岩溶作用和埋藏岩溶作用。不同地质背景下和不同成因形成的储层类型多种多样，主要包括洞穴型、孔洞型、裂缝型和裂缝-孔洞型这四类。各类储层的生产措施和产量差别较大，目前经济效益最佳的储层类型是洞穴型储层，所谓的洞穴型碳酸盐岩储层并不是无充填的洞穴，而是地下深部岩层经过挤压垮塌以后被各种杂乱堆积的角砾碎屑和流体充填的洞穴状岩体。如图I所示，在致密 灰岩I中，在箭头2所示位置发生塌陷和断裂，挤压垮塌的角砾碎屑和流体堆积形成岩溶洞穴3。勘探初期，应用常规地震属性预测技术确定储层的大致发育位置即可满足生产需求。但碳酸盐岩储层非均质性强，勘探实践证明，同一区块的相邻不同油井的储层类型和流体类型都可能有巨大的差异，实际开发工作更是提出了 “一井一藏”的高要求，因此需要针对每口井进行储集空间计算，以计算单井和区块的储量，满足地质评估需要。计算油气藏的储量的方法通常分为静态法和动态法两大类，对于均值性较强的砂岩油气藏，孔隙度等参数容易确定，这两类方法均可获得较准确的储量。但碳酸盐储层类型和物性参数都较难确定，传统的基于均值介质模型的静态储量计算方法难以获得准确的储量。动态法计算碳酸盐岩油藏储量，无论是油藏数值模拟法还是经验统计法，均需要油藏稳产相当一段时期后，才可根据生产时期取得的动态资料进行历史匹配和预测，这种时间滞后性不能满足勘探早期单井的储量评估和生产措施决策的工作需求。

发明内容
针对碳酸盐储层的非均值性，测井资料能够准确地刻画出井眼附近储层物性的纵向变化。研究表明对于单井小范围的储层而言，各种弹性参数测井(声波和密度等)之间具备确定的关系。而应用这种定量关系，就可以将地震资料和测井资料建立联系，从而可以有效地发挥地震资料横向预测范围大的优势，求得地质体在三维空间的物性参数分布情况。基于上述地震参数和测井参数存在的定量关系，本发明提出了一种获得碳酸盐岩储层储集空间的方法，以解决现有技术中存在的缺陷，从而获得更准确的单井和区块储量，满足地质评估需要。本发明提供的方法步骤包括I)在待测工区布设信号激发装置和接收装置，用所述信号激发装置产生测量用地震纵波，用接收装置采集地震纵波传输地震信号。2)对所述地震信号进行校正和保幅处理，并进行叠前时间偏移。
3)进行叠前纵波速度的反演，获得三维地震速度体，在所述速度体中根据速度值进行颜色填充的方法显示速度剖面，结合所述工区前期钻井资料进行目的层位标定和解释。 4)对所述速度剖面进行逐条观察和检视，寻找目的层段地震纵波速度异常区，确定储层发育的大致位置。其中所述的地震纵波速度异常区优选包括纵波速度轴叠置成串珠状、纵波速度轴排列成丘状或速度值突减等。5)应用测井曲线中的密度值DEN计算出孔隙度值P0R。优选通过下列公式进行计算，公式为POR = (2. 71-DEN)/1. 71其中2. 71是灰岩骨架密度值，I. 71是岩石骨架值与底层水密度值的差值。6)根据步骤5)算出的孔隙度值拟合出其与纵波速度Vp之间的定量关系；根据该定量关系将速度体计算转换为孔隙度体。所述拟合优选为投点法拟合，获得孔隙度值POR和纵波速度Vp之间的定量关系满足公式POR = a*Vp+b其中a、b为拟合直线的对应系数。7)在叠前反演的速度剖面上解释出储层的顶底界面作为计算边界。优选为人工解释确定。8)根据上述计算所得的孔隙度和对应边界计算得到待测工区储层的体积V。优选采用积分方法进行计算，优选积分公式为
NV = ^ajHiPi
i=l其中N为储层顶底内包含的体素个数，Bi是每个体素的面积，hi是每个体素的高度，Pi是每个体素的孔隙度。本发明提供的方法综合地震资料和测井资料确定碳酸盐储层的物性参数，将地球物理反演应用到油藏储量计算中，进而确立碳酸盐岩储层物性的空间分布情况，应用体素积分计算总体积的方法则克服了储层非均质性和不规则边界的影响。


此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中图I是本发明获得碳酸盐岩储层储集空间体积的方法一实施例的流程图；图2是本发明获得碳酸盐岩储层储集空间体积的方法一实施例的洞穴型碳酸盐岩储层模型示意图；图3是本发明获得碳酸盐岩储层储集空间体积的方法一实施例的叠前反演速度(Vp)剖面；图4是本发明获得碳酸盐岩储层储集空间体积的方法一实施例的测井孔隙度和测井声波时差拟合曲线。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对本发明做进一步详细说明。在此，本发明的示意性实施方式及其说明用于解释本发明，但并不作为对本发明的限定。
根据本发明提供的方法，以洞穴型碳酸盐岩储层为例进行说明。请参阅图1，图I是本发明获得碳酸盐岩储层储集空间体积的方法一实施例的流程图；步骤110 :在待测工区布设信号激发装置和接收装置，用所述信号激发装置产生测量用地震纵波，用接收装置采集地震纵波传输地震信号。步骤120 :对所述地震信号进行校正和保幅处理，并进行叠前时间偏移。步骤130 :进行叠前纵波速度的反演，获得三维地震速度体，在所述速度体中根据速度值进行颜色填充的方法显示速度剖面，结合所述工区前期钻井资料进行目的层位标定和解释。步骤140 :对所述速度剖面进行逐条观察和检视，寻找目的层段地震纵波速度异常区，确定储层发育的大致位置。其中所述的地震纵波速度异常区优选包括纵波速度轴叠置成串珠状、纵波速度轴排列成丘状或速度值突减等。步骤150 :应用测井曲线中的密度值DEN计算出孔隙度值P0R。优选通过下列公式进行计算，公式为POR = (2. 71-DEN)/1. 71其中2. 71是灰岩骨架密度值，I. 71是岩石骨架值与底层水密度值的差值。步骤160 :根据步骤150计算出的孔隙度值拟合出其与纵波速度Vp之间的定量关系；根据该定量关系将速度体计算转换为孔隙度体。所述拟合优选为投点法拟合，获得孔隙度值POR和纵波速度Vp之间的定量关系满足公式POR = a*Vp+b其中a、b为拟合直线的对应系数。步骤170 :在叠前反演的速度剖面上解释出储层的顶底界面作为计算边界。优选为人工解释确定。步骤180 :根据上述计算所得的孔隙度和对应边界计算得到待测工区储层的体积V。优选采用积分方法进行计算，优选积分公式为
权利要求
1.一种获得碳酸盐岩储层储集空间体积的方法，其特征在于，该方法包括 在待测储集工区布设信号激发装置和接收装置，用所述信号激发装置产生测量用地震纵波，用接收装置采集地震纵波传输地震信号； 对所述地震信号进行校正和保幅处理，并进行叠前时间偏移； 进行叠前纵波速度的反演，获得三维地震速度体，在所述速度体中根据速度值进行颜色填充的方法显示速度剖面，结合所述工区前期钻井资料进行目的层位标定和解释； 对所述速度剖面进行逐条观察和检视，寻找目的层段地震纵波速度异常区，确定储层发育的位置； 应用测井曲线中的密度值DEN计算出孔隙度值POR ； 根据计算出的孔隙度值拟合出其与纵波速度Vp之间的定量关系；根据该定量关系将速度体计算转换为孔隙度体； 在叠前反演的速度剖面上解释出储层的顶底界面作为计算边界； 根据上述计算所得的孔隙度和对应边界计算得到待测工区储层空间的体积。
2.根据权利要求I所述的方法，其特征在于所述的地震纵波速度异常区包括纵波速度轴叠置成串珠状、纵波速度轴排列成丘状或速度值突减。
3.根据权利要求I所述的方法，其特征在于所述密度值DEN和孔隙度值POR之间的关系通过公式POR = (2. 71-DEN)/1. 71 计算得到孔隙度值POR ;其中2. 71是灰岩骨架密度值，I. 71是岩石骨架值与底层水密度值的差值。
4.根据权利要求3所述的方法，其特征在于所述拟合为投点法拟合，获得孔隙度值POR和纵波速度Vp之间的定量关系满足公式POR = a*Vp+b 其中a、b为拟合直线的对应系数。
5.根据权利要求I所述的方法，其特征在于所述顶底界面是人工解释确定。
6.根据权利要求I所述的方法，其特征在于所述体积的计算方法为积分方法，通过公式
全文摘要
本发明提供了一种获得碳酸盐岩储层储集空间体积的方法，其具体步骤包括对待测工区的地震信号进行校正和保幅处理，进行叠前时间偏移，经过叠前纵波速度的反演，获得三维地震速度体，在所述速度体中根据速度值进行颜色填充的方法显示速度剖面，进行目的层位标定和解释；再寻找目的层段地震纵波速度异常区，确定储层发育的大致位置；根据计算出的孔隙度值拟合出其与纵波速度Vp之间的定量关系；根据该定量关系将速度体计算转换为孔隙度体；并结合解释出的储层顶底界面，计算所得的孔隙度和对应边界，最终计算得到待测工区储层空间的体积。从而可以在油田生产之前即时获得较准确的油气藏储量。
文档编号G01V1/28GK102636812SQ20121011502
公开日2012年8月15日 申请日期2012年4月18日 优先权日2012年4月18日
发明者周波, 唐俊伟, 李明, 李艳东, 陶小晚 申请人:中国石油天然气股份有限公司