本发明涉及砂砾岩油藏地质勘探领域中的测井评价技术领域，特别是涉及到一种砂砾岩有效储层常规测井定量识别方法。
背景技术：
砂砾岩油藏岩性复杂多变、储层非均质强，有效储层的判识制约了砂砾岩油藏的高效勘探开发。成像、核磁测井资料可直接对有效储层进行识别，但在实际的勘探实践中，这类资料较少，无法推广应用。因此，需要研究建立针对复杂砂砾岩油藏中有效储层的常规测井定量识别方法。
目前，已有国内外学者对砂砾岩有效储层的测井识别方法开展了较多的研究和论述。其中，大多数文献采用的仍是传统的砂砾岩测井解释评价方法，如《东营凹陷盐222块沙四段上亚段有效储层识别》(张清，2008)，《苏里格地区有效储层测井识别方法研究》(张海涛，2010)。该方法仅将储集层作为研究对象，通过建立“四性”(岩性、物性、电性、含油气性)关系的测井解释模型，采用交会图分析确定各项参数下限的半定量识别方法识别砂砾岩有效储层。但这种方法受人为因素影响和不同地区沉积条件变化的限制，对有效储层的判识准确率不高，往往造成误判，甚至会遗漏产油层段，延误勘探进程。也有学者注意到成像、核磁测井资料的优势，结合多种信息，建立了更为准确的有效储层识别模式，如《考虑岩石结构的砂砾岩有效储层测井判识方法》(闫建平等，2011)。但该方法仅限于具有成像、核磁测井资料的单井有效储层识别，无法推广至全区。为此我们发明了一种新的砂砾岩有效储层常规测井定量识别方法，解决了以上技术问题。
技术实现要素：
本发明的目的是提供一种适用于成像、核磁测井资料缺乏的情况下，利用多种常规测井资料建立计算模型，实现定量识别砂砾岩有效储层的方法。
本发明的目的可通过如下技术措施来实现：砂砾岩有效储层常规测井定量识别方法，该砂砾岩有效储层常规测井定量识别方法包括：步骤1，校正测井资料环境影响；步骤2，明确有效储层测井响应特征，优选敏感测井曲线；步骤3，计算储层识别参数r；步骤4，计算粒度指示参数g；步骤5，计算有效空间参数p；步骤6，构建储层有效性指数的计算模型e；步骤7，结合试油情况，建立有效储层定量划分标准。
本发明的目的还可通过如下技术措施来实现：
在步骤1中，以砂砾岩的实测密度、声波时差值为反演的约束条件，针对不同系列的测井仪器，优选相应的环境影响校正图版，将图版读值拟合成校正公式的基础上进行自动可视化校正处理，并考虑到图版法的局限性，采用非线性数学方法来迭代反演砂砾岩的真实密度、声波时差、自然伽马这些测井信息。
在步骤2中，明确了砂砾岩储层在常规测井曲线上具有两高两低一变化的特征，即中高电阻、高密度、低声波、低中子，自然伽马随砂砾岩母岩成分的不同而变化；优选自然伽马曲线识别储层发育程度，优选声波时差和电阻率曲线识别岩石粒度大小，选择声波时差和统计的孔隙度下限值判断有效储集空间发育程度。
在步骤3中，储层识别参数r反映了砂砾岩储层的发育程度，r值越大，砂砾岩储层越发育；反之，砂砾岩储层越不发育；采用自然伽马计算不同砾石成分的储层识别参数r。
在步骤3中，对于灰岩砾石，具有低自然伽马特征，采用公式(1)计算储层识别参数r：
式中，r为储层识别参数，无量纲；gr为实测的自然伽马值，api；grshale为纯泥岩的自然伽马值；grgravel-h为纯灰岩砾石的自然伽马值，api。
在步骤3中，对于片麻岩砾石，具有高自然伽马特征，采用公式(2)计算储层识别参数r：
式中，r为储层识别参数，无量纲；gr为实测的自然伽马值，api；grgravel-p为纯片麻岩砾石的自然伽马值，api。
在步骤4中，粒度指示参数g与岩石的粒度大小成正比，g值越大，砂砾岩的岩石粒度越大；受到岩石骨架声波时差和电阻率的影响，与岩石骨架声波时差成反比，与电阻率成正比。
在步骤4中，粒度指示参数g的计算公式(3)为：
其中：
式中，g为粒度指示参数，无量纲；rt为实测电阻率值，ω·m；δtma为岩石骨架声波时差，us/ft；△tf为流体声波时差，取180us/ft；△tac为实测声波时差，us/ft；φcnl为实测补偿中子孔隙度值，％。
在步骤5中，有效空间参数p定义为储层孔隙度与孔隙度下限的差值，p值越大，砂砾岩有效储集空间越发育；先采用试油法、束缚水饱和度法、最小有效孔喉半径法、分布函数曲线法，求取不同埋深下有效储层的孔隙度下限值，拟合有效储层孔隙度下限与埋深的关系。
在步骤5中，储层孔隙度下限值的计算公式(4)为：
φl＝-7.99×in(depth)+69.965(4)
式中，φl为储层孔隙度下限值；depth为埋深，m；
有效空间参数p定义为储层孔隙度与孔隙度下限的差值，即公式(5)：
p＝φ-φl(5)
其中，
式中，p为有效空间参数，无量纲；φ为储层孔隙度；δtma为岩石骨架声波时差，us/ft；△tf为流体声波时差，取180us/ft；△t为实测声波时差，us/ft；cp为压实校正系数，取值1；φl为储层孔隙度下限值。
在步骤6中，构建的储层有效性指数的计算模型，即公式(6)为：
e＝r×g×p(6)
储层有效性指数e越大，表明砂砾岩储层的各类沉积和储集条件有利，成藏及产能情况较好。
在步骤7中，将储层有效性指数计算模型应用于不同典型井中，得到典型井的储层有效性指数值，计算特定埋深的单井极限初产油流，确定优质储层日产油量，结合实际井段的试油情况，建立有效储层的定量划分标准。
在步骤7中，建立的有效储层的定量划分标准为：e值大于0.25时，成藏条件最好，不采取措施或压裂后稳产，为i类储层；e值介于0.25至0.1之间时，成藏条件较好，经酸化压裂后增产明显，为ii类储层；e值介于0.1至0.05之间时，成藏条件较差，经酸化压裂后效果较差，为iii类储层；e值介于0至0.05之间时，成藏条件最差，酸化压裂无效，为iv类储层；依据此标准实现对砂砾岩有效储层的定量识别。
本发明中的砂砾岩有效储层常规测井定量识别方法，在明确砂砾岩有效储层常规测井响应特征基础上，优选测井曲线分别计算了储层识别参数、粒度指示参数和有效空间参数，并构建了储层有效性指数的计算模型；最终结合试油情况，建立有效储层划分识别标准，达到了定量识别砂砾岩有效储层的目的。这种砂砾岩有效储层常规测井定量识别方法以典型井试油井段的测井响应特征为依据，融合多种常规测井曲线，建立了有效储层的定量公式和划分识别标准，提高了砂砾岩有效储层的定量识别精度，可在成像、核磁测井资料缺乏的情况下，为复杂砂砾岩油藏中有效储层的划分识别和试油层段优选提供重要参考依据，具有良好的应用效果和推广前景。
附图说明
图1为本发明的砂砾岩有效储层常规测井定量识别方法的一具体实施例的流程图；
图2为本发明的一具体实施例中储层识别参数r在d7井3880～4060m的识别效果图；
图3为本发明的一具体实施例中不同岩石粒度的电阻率与岩石骨架的交会分析图；
图4为本发明的一具体实施例中粒度指示参数g在c6井4340～4380m的识别效果与成像测井对比图；
图5为本发明的一具体实施例中实施例地区中有效储层孔隙度下限与深度关系图；
图6为本发明的一具体实施例中有效空间参数p在d7井3900～4070m的识别效果图；
图7为本发明的一具体实施例中ⅰ类储层的储层有效性指数e的储层定量识别成果图；
图8为本发明的一具体实施例中ⅱ类储层的储层有效性指数e的储层定量识别成果图；
图9为本发明的一具体实施例中ⅲ类储层的储层有效性指数e的储层定量识别成果图；
图10为本发明的一具体实施例中ⅳ类储层的储层有效性指数e的储层定量识别成果图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施例，并配合附图所示，作详细说明如下。
如图1所示，图1为本发明的砂砾岩有效储层常规测井定量识别方法的流程图。
在步骤101，对测井资料环境影响进行校正，保证用于有效储层定量识别的测井资料真实可靠。以砂砾岩的实测密度、声波时差值为反演的约束条件，针对不同系列的测井仪器，优选相应的环境影响校正图版，将图版读值拟合成校正公式的基础上进行计算机自动可视化校正处理，并考虑到图版法的局限性，采用非线性数学方法来迭代反演砂砾岩的真实密度、声波时差、自然伽马等测井信息，以保障用于有效储层定量识别参数计算的测井资料真实可靠。流程进入到步骤102。
在步骤102，明确砂砾岩有效储层的常规测井响应特征，优选敏感测井曲线。在此实施例中，砂砾岩储层在常规测井曲线上具有“两高两低一变化”的特征，即中高电阻、高密度、低声波、低中子，自然伽马随砂砾岩母岩成分的不同而变化。最终，优选自然伽马曲线识别储层发育程度，优选声波时差和电阻率曲线识别岩石粒度大小，选择声波时差和统计的孔隙度下限值判断有效储集空间发育程度。流程进入到步骤103。
在步骤103，计算储层识别参数r，该参数主要反映了砂砾岩储层的发育程度，r值越大，砂砾岩储层越发育；反之，砂砾岩储层越不发育。主要采用自然伽马计算不同砾石成分的储层识别参数r。对于灰岩砾石，具有低自然伽马特征，采用公式(1)计算储层识别参数r：
对于片麻岩砾石，具有高自然伽马特征，采用公式(2)计算储层识别参数r：
式中，r为储层识别参数，无量纲；gr为实测的自然伽马值，api；grshale为纯泥岩的自然伽马值；grgravel-h为纯灰岩砾石的自然伽马值，api；
grgravel-p为纯片麻岩砾石的自然伽马值，api。流程进入到步骤104。
在步骤104，计算粒度指示参数g，该参数与岩石的粒度大小成正比，g值越大，砂砾岩的岩石粒度越大。主要受到岩石骨架声波时差和电阻率的影响，与岩石骨架声波时差成反比，与电阻率成正比，即公式(3)：
其中：
式中，g为粒度指示参数，无量纲；rt为实测电阻率值，ω·m；δtma为岩石骨架声波时差，us/ft；△tf为流体声波时差，取180us/ft；△tac为实测声波时差，us/ft；φcnl为实测补偿中子孔隙度值，％。流程进入到步骤105。
在步骤105，计算有效空间参数p，该参数定义为储层孔隙度与孔隙度下限的差值，p值越大，砂砾岩有效储集空间越发育。需要先采用试油法、束缚水饱和度法、最小有效孔喉半径法、分布函数曲线法，求取不同埋深下有效储层的孔隙度下限值，拟合有效储层孔隙度下限与埋深的关系，即公式(4)：
φl＝-7.99×in(depth)+69.965(4)
式中，φl为储层孔隙度下限值；depth为埋深，m。
有效空间参数p定义为储层孔隙度与孔隙度下限的差值，即公式(5)：
p＝φ-φl(5)
其中，
式中，p为有效空间参数，无量纲；φ为储层孔隙度；δtma为岩石骨架声波时差，us/ft；△tf为流体声波时差，取180us/ft；△t为实测声波时差，us/ft；cp为压实校正系数，取值1；φl为储层孔隙度下限值。流程进入到步骤106。
在步骤106，构建储层有效性指数的计算模型，即公式(6)：
e＝r×g×p(6)
储层有效性指数e越大，表明砂砾岩储层的各类沉积和储集条件有利，成藏及产能情况较好。流程进入到步骤107。
在步骤107，结合试油情况，建立有效储层定量划分标准。将储层有效性指数计算模型应用于不同典型井中，得到典型井的储层有效性指数值。依据中华人民共和国地质矿产行业标准(dz/t0217-2005)石油天然气储量计算规范，依据实施例地区所用成本油价、钻井投资、地面投资等计算实施例地区埋深3000～4500m的单井极限初产油为2.5～3t/d，将初产油大于10t/d定为优质储层。据此，结合实际井段的试油情况，可建立有效储层的定量划分标准：
e值大于0.25时，成藏条件最好，试油日产大于10吨，往往不采取措施或压裂后稳产，为i类储层；e值介于0.25至0.1之间时，成藏条件较好，试油日产大于10吨，经酸化压裂后增产明显，为ii类储层；e值介于0.1至0.05之间时，成藏条件较差，试油日产在0.5～10吨，经酸化压裂后效果较差，为iii类储层；e值介于0至0.05之间时，成藏条件最差，试油日产小于0.5吨，酸化压裂无效，为iv类储层。依据此标准可实现对砂砾岩有效储层的定量识别。
图2是储层识别参数r在d7井3880～4060m的识别效果图，该参数可有效区分片麻岩砾石(硅质含量较多)和灰岩砾石(灰质含量较多)。图3是不同岩石粒度的电阻率与岩石骨架的交会分析图，该图表明岩石粒度的大小主要与岩石骨架声波时差成反比，与电阻率成正比。图4是粒度指示参数g在c6井4340～4380m的识别效果与成像测井对比图，可以看出岩石粒度定量识别结果与实测的成像测井具有较好的对应关系。图5是实施例地区中有效储层孔隙度下限与深度关系图，图6是有效空间参数p在d7井3900～4070m的识别效果图。图5和图6表明通过计算孔隙度和统计的孔隙度下限两者相结合，可以有效识别储层中的有效储集空间。图7-10是储层有效性指数e的四类储层定量识别成果图，利用建立的分类识别标准推广至全区中判别精度达到80％以上，实际应用效果较好，实现了砂砾岩有效储层的常规测井定量识别。
本说明书中采用的表示方法，是本领域技术人员习惯用法，本领域技术人员熟知，不做更详细解释。
如上所述，对本发明的实施例进行了详细地说明，并非用于限定本发明的保护范围。只要实质上没有脱离本发明的发明点及效果可以有很多的变形，这对本领域的技术人员来说是显而易见的。因此，这样的变形例也全部包含在本发明的保护范围之内。