砂岩储层的孔隙度的预测方法、装置及计算机存储介质制造方法及图纸

本说明书提供砂岩储层的孔隙度的预测方法、装置及计算机存储介质，包括：根据目标储层的原始沉积体的参数，计算目标储层的原始孔隙度；其中，所述原始沉积体用于表示在沉积末期，埋藏初期所形成的沉积物；根据目标储层的埋深，计算在指定埋深作用下目标储层的孔隙度；根据目标储层的构造应力，计算在指定构造应力作用下目标储层的孔隙度减孔量；根据实测砂岩储层的孔隙度减孔量计算埋深作用下的减孔系数和构造应力作用下的减孔系数；根据所述目标储层的原始孔隙度、埋深作用下目标储层的孔隙度、构造应力作用下目标储层的孔隙度减孔量和埋深作用下的减孔系数以及构造应力作用下的减孔系数预测所述目标储层的孔隙度。提高了孔隙度计算精度。

【技术实现步骤摘要】
砂岩储层的孔隙度的预测方法、装置及计算机存储介质
本说明书涉及石油天然气勘探领域，特别涉及砂岩储层的孔隙度的预测方法、装置及计算机存储介质。
技术介绍
孔隙度是指岩样中所有孔隙空间体积之和与该岩样体积的比值，储层的孔隙度越大，说明岩石中孔隙空间越大，孔隙不仅能够储存油气，而且可以允许油气在其中渗滤，孔隙度是储层表征的重要参数，储层孔隙度的大小直接代表了储层储集性能的好坏。致密砂岩储层是非常重要的油气储层，属于非常规油气储层，按照国内为对超深层储层的定义，超过3500-4000m的储层为深层储层。二十世纪七十年代，美国联邦能源管理委员会将储层基质渗透率小于0.1×10-3μm2的气藏定义为致密气藏。在我国，塔里木盆地白垩系巴什基奇克组砂岩储层埋深达6000-7000m，样品渗透率小于0.1×10-3μm2，属超深层致密砂岩储层。近年来我国非常规油气资源勘探开发取得了重大突破，新增探明储量逐年增加，占我国油气新增探明储量的60％以上。因此研究塔里木盆地白垩系巴什基奇克组超深层致密砂岩储层的储集能力，预测计算其孔隙度对指导塔里木盆地白垩系巴什基奇克组砂岩储层的油气资源勘探与开发具有十分重要的意义。在实现本说明书的过程中，专利技术人发现现有技术中存在如下问题。目前，国内外对储层孔隙度的计算预测方法大致可以分为深度指数法以及基于深度指数法的相应孔隙度预测方法；基于成岩物理实验的孔隙度预测方法；基于储层孔隙演化史孔隙度模拟方法；基于地震相的储层孔隙度预测方法以及基于局部的决策树算法的孔隙度预测方法。上述方法对储层孔隙度的预测大多基于储层成岩演化过程或是基于储层的地震属性，鉴于影响储层孔隙度发育因素的多样性，上述预测方法存在不足。
技术实现思路
本说明书实施方式提供砂岩储层的孔隙度的预测方法、装置及计算机存储介质。通过计算得到目标储层原始孔隙度、埋藏压实作用下目标储层的孔隙度、构造挤压作用下目标储层的孔隙度以及埋藏压实作用下的减孔系数和和构造挤压作用下的减孔系数，从而预测得到目标储层的孔隙度。本说明书实施方式提供一种砂岩储层的孔隙度的预测方法。包括，根据目标储层的原始沉积体的参数，计算目标储层的原始孔隙度；其中，所述原始沉积体用于表示在沉积末期，埋藏初期所形成的沉积物；所述原始孔隙度用于表示所述目标储层在沉积末期，埋藏初期时的孔隙度；根据目标储层的埋深，计算在指定埋深作用下目标储层的孔隙度；根据目标储层的构造应力，计算在指定构造应力作用下目标储层的孔隙度减孔量；根据实测砂岩储层的孔隙度减孔量计算埋深作用下的减孔系数和构造应力作用下的减孔系数；根据所述目标储层的原始孔隙度、埋深作用下目标储层的孔隙度、构造应力作用下目标储层的孔隙度减孔量和埋深作用下的减孔系数以及构造应力作用下的减孔系数预测所述目标储层的孔隙度。本说明书实施方式提供一种砂岩储层的孔隙度的预测装置。包括，原始孔隙度计算模块；根据目标储层的原始沉积体的参数，计算目标储层的原始孔隙度；其中，所述原始沉积体用于表示在沉积末期，埋藏初期所形成的沉积物；所述原始孔隙度用于表示所述目标储层在沉积末期，埋藏初期时的孔隙度；埋深作用下孔隙度计算模块；用于根据目标储层的埋深，计算在指定埋深作用下目标储层的孔隙度；构造应力作用下孔隙度减孔量计算模块；用于根据目标储层的构造应力，计算在指定构造应力作用下目标储层的孔隙度减孔量；减孔系数计算模块；根据实测砂岩储层的孔隙度减孔量计算埋深作用下的减孔系数和构造应力作用下的减孔系数；孔隙度计算模块；用于根据所述目标储层的原始孔隙度、埋深作用下目标储层的孔隙度、构造应力作用下目标储层的孔隙度减孔量和埋深作用下的减孔系数以及构造应力作用下的减孔系数预测所述目标储层的孔隙度。本说明书实施方式提供一种计算机存储介质。所述计算机存储介质存储有计算机程序指令，在所述计算机程序指令被执行时实现：根据目标储层的原始沉积体的参数，计算目标储层的原始孔隙度；其中，所述原始沉积体用于表示在沉积末期，埋藏初期所形成的沉积物；所述原始孔隙度用于表示所述目标储层在沉积末期，埋藏初期时的孔隙度；根据目标储层的埋深，计算在指定埋深作用下目标储层的孔隙度；根据目标储层的构造应力，计算在指定构造应力作用下目标储层的孔隙度减孔量；根据实测砂岩储层的孔隙度减孔量计算埋深作用下的减孔系数和构造应力作用下的减孔系数；根据所述目标储层的原始孔隙度、埋深作用下目标储层的孔隙度、构造应力作用下目标储层的孔隙度减孔量和埋深作用下的减孔系数以及构造应力作用下的减孔系数预测所述目标储层的孔隙度。本说明书的有益效果在于基于区域大量可靠的实验资料、砂岩动力成岩理论、数值分析和地质统计学方法，通过计算得到目标储层原始孔隙度、埋深作用下目标储层的孔隙度、构造应力作用下目标储层的孔隙度以及埋深作用下的减孔系数和构造应力作用下的减孔系数，综合计算目标储层的孔隙度。本说明书实施方式考虑了砂岩原始沉积体对储层原生孔隙及其演化的影响，还考虑了后期埋藏压实以及构造挤压作用对储层孔隙度的影响，并计算埋藏压实与构造挤压作用下的减孔系数，最终应用预测孔隙度，解决了现有技术中砂岩储层孔隙度预测过程中因为没有考虑到孔隙度发育因素多样性所带来的不足。附图说明为了更清楚地说明本说明书实施方式或现有技术中的技术方案，下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本说明书中记载的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1为本说明书实施方式中提供的砂岩储层的孔隙度的预测方法的流程图；图2为本说明书实施方式中提供的砂岩储层的孔隙度的预测装置的结构框图；图3为本说明书实施方式中提供的砂岩储层孔隙度与平均粒径拟合图；图4为本说明书实施方式中提供的砂岩储层孔隙度与泥质含量拟合图；图5为本说明书实施方式中提供的砂岩储层孔隙度与分选系数拟合图；图6为本说明书实施方式中提供的砂岩储层孔隙度与刚性颗粒含量拟合图；图7为本说明书实施方式中提供的实验条件下砂岩储层孔隙度随埋深变化图；图8为本说明书实施方式中提供的砂岩储层孔隙度随埋深变化图；图9为本说明书实施方式中提供的砂岩储层孔隙度随构造应力变化图；图10为本说明书实施方式中提供的砂岩储层的埋藏史线、孔隙史线、地温史线。具体实施方式为了使本
的人员更好地理解本说明书中的技术方案，下面将结合本说明书实施方式中的附图，对本说明书实施方式中的技术方案进行清楚、完整地描述，显然，所描述的实施方式仅仅是本说明书一部分实施方式，而不是全部的实施方式。基于本说明书中的实施方式，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施方式，都应当属于本说明书保护的范围。在一个场景实施方式中，已有库车坳陷秋探1井、东秋8井白垩系巴什基奇克组储集层孔隙度，通过使用本说明书提供的实施方式对库车坳陷秋探1井、东秋8井白垩系巴什基奇克组储集层的孔隙度进行计算，并与已有的储集层孔隙度进行比较，验证孔隙度计算精度。在本场景实施方式中，已有克拉2气田白垩系砂岩储层的地质资料，包括岩屑、岩心粒度分析结果，岩石矿物成分、埋深等资料。经过资料分析，该区域在埋深3600-4000米有残余原生孔隙型储层本文档来自技高网...

【技术保护点】
1.一种砂岩储层的孔隙度的预测方法，其特征在于，包括：根据目标储层的原始沉积体的参数，计算目标储层的原始孔隙度；其中，所述原始沉积体用于表示在沉积末期，埋藏初期所形成的沉积物；所述原始孔隙度用于表示所述目标储层在沉积末期，埋藏初期时的孔隙度；根据目标储层的埋深，计算在指定埋深作用下目标储层的孔隙度；根据目标储层的构造应力，计算在指定构造应力作用下目标储层的孔隙度减孔量；根据实测砂岩储层的孔隙度减孔量计算埋深作用下的减孔系数和构造应力作用下的减孔系数；根据所述目标储层的原始孔隙度、埋深作用下目标储层的孔隙度、构造应力作用下目标储层的孔隙度减孔量和埋深作用下的减孔系数以及构造应力作用下的减孔系数预测所述目标储层的孔隙度。

【技术特征摘要】
1.一种砂岩储层的孔隙度的预测方法，其特征在于，包括：根据目标储层的原始沉积体的参数，计算目标储层的原始孔隙度；其中，所述原始沉积体用于表示在沉积末期，埋藏初期所形成的沉积物；所述原始孔隙度用于表示所述目标储层在沉积末期，埋藏初期时的孔隙度；根据目标储层的埋深，计算在指定埋深作用下目标储层的孔隙度；根据目标储层的构造应力，计算在指定构造应力作用下目标储层的孔隙度减孔量；根据实测砂岩储层的孔隙度减孔量计算埋深作用下的减孔系数和构造应力作用下的减孔系数；根据所述目标储层的原始孔隙度、埋深作用下目标储层的孔隙度、构造应力作用下目标储层的孔隙度减孔量和埋深作用下的减孔系数以及构造应力作用下的减孔系数预测所述目标储层的孔隙度。2.根据权利要求1所述的方法，其特征在于，所述目标储层的原始沉积体的参数至少包括：所述原始沉积体的平均粒径、分选系数、泥质含量以及刚性颗粒含量数据。3.根据权利要求2所述的方法，其特征在于，计算所述原始孔隙度包括：根据实测砂岩储层原始沉积体的孔隙度、平均粒径、分选系数、泥质含量与刚性颗粒含量通过回归算法得到所述平均粒径、分选系数、泥质含量与刚性颗粒含量对应的回归系数；根据与所述目标储层关联的指定地层的平均粒径、分选系数、泥质含量与刚性颗粒含量以及所述回归系数计算所述原始孔隙度。4.根据权利要求3所述的方法，其特征在于，所述原始孔隙度根据公式：计算得出；其中，表示目标储层的原始孔隙度，So表示分选系数，D表示平均粒径、Vsh表示泥质含量、VQ表示刚性颗粒含量，a、b、c、d、e分别表示回归系数。5.根据权利要求1所述的方法，其特征在于，计算所述指定埋深作用下目标储层的孔隙度，包括：根据实测砂岩储层孔隙度和埋深数据，通过回归算法计算回归系数；根据所述目标储层的埋深和所述回归系数计算在指定埋深作用下目标储层的孔隙度。6.根据权利要求5所述的方法，其特征在于，所述方法还包括，根据不同的砂岩类型计算不同砂岩类型对应的回归系数。7.根据权利要求5所述方法，其特征在于，所述指定埋深作用下目标储层的孔隙度，根据公式：计算得出；其中f、g表示回归系数，H表示所述指定埋深。8.根据权利要求1所述的方法，其特征在于，所述计算在指定构造应力作用下目标储层的孔隙度减孔量，包括：根据实测砂岩储层孔隙度和构造应力数据，通过回归算法计算回归系数；根据所述目标储层所受的构造应力和所述回归系数计算在指定构造应力作用下目标储层的孔隙度。9.根据权利要求8所述的方法，其特征在于，所述方法还包括，根据不同砂岩类型计算不同砂岩类型对应的回归系数。10.根据权利要求8所述的方法，其特征在于，所述计算在指定构造应力作用下目标储层的孔隙度减孔量，根据公式计算得出；其中，h、i表示回归系数，T表示所述指定构造应力。11.根据权利要求1所述的方法，其特征在于，计算所述埋深作用下的减孔系数和构造挤压作用下的减孔系数，包括：根据实测砂岩储层的地史数据，生成所述实测砂岩储层的埋藏史线、孔隙史线、地温...

【专利技术属性】
技术研发人员：张荣虎，斯春松，曾庆鲁，刘春，王俊鹏，陈戈，张惠良，寿建峰，王波，赵继龙，王珂，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京,11