砂砾岩体无效储层扣除方法与流程

本发明涉及油田开发技术领域，特别是涉及到一种砂砾岩体无效储层扣除方法。

背景技术：

砂砾岩体储层广泛发育，是目前及今后一段时间内勘探与增储上产的重要领域，在我国第二大油田——胜利油田中，砂砾岩体油气资源量占有相当大的比例。因此，确定砂砾岩有效储层分布范围及规模，不仅是此类领域勘探目标优选的依据，也是准确计算油气储量的基础。受沉积分异、成岩作用等的影响，同一砂砾岩体的不同部位储集物性存在很大差异，局部物性低于成藏下限，即为无效储层，以往的方法主要就是经验法、测试法、含油产状法，这些方法主观性强，缺乏明确、可靠、科学的操作步骤。为此我们发明了一种新的砂砾岩体无效储层扣除方法，解决了以上技术问题。

技术实现要素：

本发明的目的是提供一种砂砾岩体无效储层扣除方法，解决砂砾岩体油气勘探过程中，确定砂砾岩有效储层分布范围的问题。

本发明的目的可通过如下技术措施来实现：砂砾岩体无效储层扣除方法，该砂砾岩体无效储层扣除方法包括：(1)获取钻井的常规测井资料、录井岩性资料、试油资料；(2)计算钻井钻遇的砂砾岩体的孔隙度、渗透率；(3)将步骤(2)得到的孔隙度、渗透率的计算结果绘制在半对数坐标系内；(4)将砂砾岩体中的非有效储层和有效储层分界处对应的孔隙度和渗透率值读出即为有效储层物性下限值；(5)统计钻井钻遇的砂砾岩体厚度结合地震相特征绘制砂砾岩体平面分布图；(6)将钻井的孔隙度、渗透率计算数据标注在步骤(5)得到砂砾岩体平面分布图上并绘制砂砾岩体物性等值线图；(7)以步骤(4)确定的有效储层物性下限值为界限，在砂砾岩体物性等值线图上扣除物性下限值以下的砂砾岩体无效储层确定有效储层平面展布。

本发明的目的还可通过如下技术措施来实现：

在步骤2中，根据泥岩含量SH，通过测井资料中的密度测井DEN、声波测井AC、补偿中子测井CNL三种测井资料中的任意一种计算砂砾岩体的孔隙 度POR、渗透率PERM。

在步骤2中，通过密度测井DEN计算砂砾岩体的孔隙度POR、渗透率PERM，计算公式为：

$\mathrm{POR}=\frac{\mathrm{DEN}-\mathrm{DG}}{\mathrm{DF}-\mathrm{DG}}-\mathrm{SF}\×\frac{\mathrm{DSH}-\mathrm{DG}}{\mathrm{DF}-\mathrm{DG}}$

式中DEN为密度测井曲线值；DG为岩石骨架密度，取值2.65；DF为流体密度，取值1；DSH为泥岩密度，取值2.65；SH为泥岩含量；

$\mathrm{PERM}=0.136\×\frac{{(100\×\mathrm{POR})}^{4.4}}{{\mathrm{SIRR}}^2}$

式中SIRR为束缚水饱和度，取值40。

在步骤2中，通过声波测井AC计算砂砾岩体的孔隙度POR、渗透率PERM，计算公式为：

$\mathrm{POR}=\frac{\mathrm{AC}-\mathrm{TM}}{\mathrm{TF}-\mathrm{TM}}\×\frac{1}{\mathrm{CP}1}-\mathrm{SH}\×\frac{\mathrm{TSH}-\mathrm{TM}}{\mathrm{TF}-\mathrm{TM}}$

式中AC为声波测井曲线值；TM为岩石骨架声波时差，取值180；TF为流体声波时差，取值620；CP1为压实系数，取值1；TSH为泥岩声波时差，取值180；SH为泥岩含量；

$\mathrm{PERM}=0.136\×\frac{{(100\×\mathrm{POR})}^{4.4}}{{\mathrm{SIRR}}^2}$

式中SIRR为束缚水饱和度，取值40。

在步骤2中，通过补偿中子测井CNL计算砂砾岩体的孔隙度POR、渗透率PERM，计算公式为：

POR＝ＣＮＬ－ＳＨ×ＮＳＨ

式中CNL为补偿中子测井曲线值；NSH为泥岩中子，取值28；SH为泥岩含量；

$\mathrm{PERM}=0.136\×\frac{{(100\×\mathrm{POR})}^{4.4}}{{\mathrm{SIRR}}^2}$

式中SIRR为束缚水饱和度，取值40。

在步骤2中，泥岩含量SH通过测井资料中自然伽马测井GR、补偿中子测井CNL、自然电位测井SP、中子寿命测井NLL和深探测电阻率测井RT五种测井资 料中的任意一种计算，由于可以使用5种测井资料中任意一种进行计算，因此在公式中所选择计算的测井曲线的最大值用GMXi代表；所选择计算的测井曲线的最小值用GMNi代表，公式为：

$\mathrm{SH}=\frac{{2.0}^{(\mathrm{GCUR}\×\mathrm{VSH})}-1.0}{{2.0}^{\mathrm{GCUR}}-1.0}$

$\mathrm{VSH}=\frac{\mathrm{SHLG}-\mathrm{GMNi}}{\mathrm{GMXi}-\mathrm{GMNi}}$

式中GCUR为计算泥质含量经验系数，对于新生界取值1，对于中生界、古生界和太古界取值2；SHLG为所选择计算的测井曲线值。

在步骤4中，非有效储层和有效储层判断方法是根据试油结果，干层为非有效储储层，油层、低产油层、油水同层、含油水层、水层为有效储层。

在步骤5中，砂砾岩体地震相特征是，在地震剖面上砂砾岩体呈现同向轴不连续的杂乱反射或空白反射特征，与一般正常地层同向轴连续成层性好的反射特征有明显区别。

该砂砾岩体无效储层扣除方法涉及油气勘探领域钻探目标选取及储量计算的基础工作——有效储层分布范围的圈定，主要针对盆地陡坡带近岸水下扇砂砾岩储层，提供了一种扣除无效储层、圈定有效储层分布范围的方法。本发明中的砂砾岩体无效储层扣除方法，通过对(近岸水下扇)砂砾岩含油气性、物性、沉积展布的研究，用试油法确定有效储层界限明确有效储层的物性下限，该方法比测试法简单、方便，同时又比含油产状法准确可靠。最后结合钻井与地震相特征圈定无效储层的分布范围。该方法结合有效储层物性界线和扇体展布分析结果，可以确定有效储层平面展布，扣除无效储层，从而可以更好地指导油气勘探和储量申报。

附图说明

图1为本发明的砂砾岩体无效储层扣除方法的一具体实施例的流程图；

图2为本发明的一具体实施例的济阳坳陷渤南洼陷北部砂砾岩体孔隙度和渗透率交汇图(半对数坐标)；

图3为本发明的一具体实施例的济阳坳陷渤南洼陷北部砂砾岩体地震剖面图；

图4为本发明的一具体实施例中渤南洼陷北部砂砾岩体有效储层分布 图。

具体实施方式

为使本发明的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施例，并配合所附图式，作详细说明如下。

如图1所示，图1为本发明的砂砾岩体无效储层扣除方法的流程图。

在步骤101，获取钻井的常规测井资料、录井岩性资料、试油资料。流程进入到步骤102。

在步骤102，计算钻井钻遇的砂砾岩体的孔隙度、渗透率；

计算砂砾岩体的孔隙度POR、渗透率PERM，是通过测井资料中的密度测井DEN、声波测井AC、补偿中子测井CNL三种测井资料中的任意一种计算出来的，公式分别为：

$\mathrm{POR}=\frac{\mathrm{DEN}-\mathrm{DG}}{\mathrm{DF}-\mathrm{DG}}-\mathrm{SF}\×\frac{\mathrm{DSH}-\mathrm{DG}}{\mathrm{DF}-\mathrm{DG}}$

式中DEN为密度测井曲线值；DG为岩石骨架密度，取值2.65；DF为流体密度，取值1；DSH为泥岩密度，取值2.65；SH为泥岩含量。

$\mathrm{POR}=\frac{\mathrm{AC}-\mathrm{TM}}{\mathrm{TF}-\mathrm{TM}}\×\frac{1}{\mathrm{CP}1}-\mathrm{SH}\×\frac{\mathrm{TSH}-\mathrm{TM}}{\mathrm{TF}-\mathrm{TM}}$

式中AC为声波测井曲线值；TM为岩石骨架声波时差，取值180；TF为流体声波时差，取值620；CP1为压实系数，取值1；TSH为泥岩声波时差，取值180；SH为泥岩含量。

POR＝GNL-SH×NSH

式中CNL为补偿中子测井曲线值；NSH为泥岩中子，取值28；SH为泥岩含量。

$\mathrm{PERM}=0.136\×\frac{{(100\×\mathrm{POR})}^{4.4}}{{\mathrm{SIRR}}^2}$

式中SIRR为束缚水饱和度，取值40。

三种计算公式中都有泥岩含量SH需要计算，由于各种测井资料可能不全，因此为了保证本专利的可行性，泥岩含量SH可以通过测井资料中自然伽马测井GR、补偿中子测井CNL、自然电位测井SP、中子寿命测井NLL和深探测电阻率测井RT五种测井资料中的任意一种计算，由于可以使用5种测井资料中任 意一种进行计算，因此在公式中所选择计算的测井曲线的最大值用GMXi代表；所选择计算的测井曲线的最小值用GMNi代表。公式为：

$\mathrm{SH}=\frac{{2.0}^{(\mathrm{GCUR}\×\mathrm{VSH})}-1.0}{{2.0}^{\mathrm{GCUR}}-1.0}$

$\mathrm{VSH}=\frac{\mathrm{SHLG}-\mathrm{GMNi}}{\mathrm{GMXi}-\mathrm{GMNi}}$

式中GCUR为计算泥质含量经验系数，对于新生界取值1，对于中生界、古生界和太古界取值2；SHLG为所选择计算的测井曲线值；GMXi为所选择计算的测井曲线的最大值；GMNi为所选择计算的测井曲线的最小值。流程进入到步骤103。

在步骤103，将步骤102得到的孔隙度、渗透率的计算结果绘制在半对数坐标系内，如图2所示。流程进入到步骤104。

在步骤104，将砂砾岩体中的非有效储层和有效储层分界处对应的孔隙度和渗透率值读出即为有效储层物性下限值，如图2所示；

非有效储层和有效储层判断方法是根据试油结果，干层为非有效储储层，油层、低产油层、油水同层、含油水层、水层为有效储层。流程进入到步骤105。

在步骤105，统计钻井钻遇的砂砾岩体厚度结合地震相特征绘制砂砾岩体平面分布图；

砂砾岩体地震相特征是，在地震剖面上砂砾岩体呈现同向轴不连续的杂乱反射或空白反射特征。与一般正常地层同向轴连续成层性好的反射特征有明显区别，如图3所示，图3是说明砂砾岩体的地震相特征，用以明确砂砾岩体的展布范围，并在平面成图。流程进入到步骤106。

在步骤106，将钻井的孔隙度、渗透率计算数据标注在步骤105得到砂砾岩体平面分布图上并绘制砂砾岩体物性等值线图，如图4所示。流程进入到步骤107。

在步骤107，以步骤104确定的有效储层物性下限值为界限，在砂砾岩体物性等值线图上扣除物性下限值以下的砂砾岩体无效储层确定有效储层平面展布，如图4所示。