1.砂砾岩体无效储层扣除方法,其特征在于,该砂砾岩体无效储层扣除方法包括:
(1)获取钻井的常规测井资料、录井岩性资料、试油资料;
(2)计算钻井钻遇的砂砾岩体的孔隙度、渗透率;
(3)将步骤(2)得到的孔隙度、渗透率的计算结果绘制在半对数坐标系内;
(4)将砂砾岩体中的非有效储层和有效储层分界处对应的孔隙度和渗透率值读出即为有效储层物性下限值;
(5)统计钻井钻遇的砂砾岩体厚度结合地震相特征绘制砂砾岩体平面分布图;
(6)将钻井的孔隙度、渗透率计算数据标注在步骤(5)得到砂砾岩体平面分布图上并绘制砂砾岩体物性等值线图;
(7)以步骤(4)确定的有效储层物性下限值为界限,在砂砾岩体物性等值线图上扣除物性下限值以下的砂砾岩体无效储层确定有效储层平面展布。
2.根据权利要求1所述的砂砾岩体无效储层扣除方法,其特征在于,在步骤2中,根据泥岩含量SH,通过测井资料中的密度测井DEN、声波测井AC、补偿中子测井CNL三种测井资料中的任意一种计算砂砾岩体的孔隙度POR、渗透率PERM。
3.根据权利要求2所述的砂砾岩体无效储层扣除方法,其特征在于,在步骤2中,通过密度测井DEN计算砂砾岩体的孔隙度POR、渗透率PERM,计算公式为:
式中DEN为密度测井曲线值;DG为岩石骨架密度,取值2.65;DF为流体密度,取值1;DSH为泥岩密度,取值2.65;SH为泥岩含量;
式中SIRR为束缚水饱和度,取值40。
4.根据权利要求2所述的砂砾岩体无效储层扣除方法,其特征在于,在步骤2中,通过声波测井AC计算砂砾岩体的孔隙度POR、渗透率PERM,计算公 式为:
式中AC为声波测井曲线值;TM为岩石骨架声波时差,取值180;TF为流体声波时差,取值620;CP1为压实系数,取值1;TSH为泥岩声波时差,取值180;SH为泥岩含量;
式中SIRR为束缚水饱和度,取值40。
5.根据权利要求2所述的砂砾岩体无效储层扣除方法,其特征在于,在步骤2中,通过补偿中子测井CNL计算砂砾岩体的孔隙度POR、渗透率PERM,计算公式为:
POR=CNL-SH×NSH
式中CNL为补偿中子测井曲线值;NSH为泥岩中子,取值28;SH为泥岩含量;
式中SIRR为束缚水饱和度,取值40。
6.根据权利要求2所述的砂砾岩体无效储层扣除方法,其特征在于,在步骤2中,泥岩含量SH通过测井资料中自然伽马测井GR、补偿中子测井CNL、自然电位测井SP、中子寿命测井NLL和深探测电阻率测井RT五种测井资料中的任意一种计算,由于可以使用5种测井资料中任意一种进行计算,因此在公式中所选择计算的测井曲线的最大值用GMXi代表;所选择计算的测井曲线的最小值用GMNi代表,公式为:
式中GCUR为计算泥质含量经验系数,对于新生界取值1,对于中生界、古生界和太古界取值2;SHLG为所选择计算的测井曲线值。
7.根据权利要求1所述的砂砾岩体无效储层扣除方法,其特征在于,在步骤4中,非有效储层和有效储层判断方法是根据试油结果,干层为非有效储储层,油层、低产油层、油水同层、含油水层、水层为有效储层。
8.根据权利要求1所述的砂砾岩体无效储层扣除方法,其特征在于,在步骤5中,砂砾岩体地震相特征是,在地震剖面上砂砾岩体呈现同向轴不连续的杂乱反射或空白反射特征,与一般正常地层同向轴连续成层性好的反射特征有明显区别。