一种基于密闭取心的原始油藏含水饱和度计算方法与流程

本发明属于石油勘探中的复杂油气藏储层流体性质测井评价技术领域，涉及一种基于密闭取心的原始油藏含水饱和度计算方法。

背景技术：

毛管压力曲线是评价储层原始含油饱和度的主要方法之一，通过对目的层段选取有限的岩心开展压汞实验，采用j函数法或者基于流动单元分析法将毛管压力资料进行分类得到相同类型的平均毛管压力曲线，利用平均毛管压力曲线求取油藏原始含油饱和度；或者建立不同进汞压力下进汞饱和度与砂岩渗透率之间函数模型,对整个储层毛管压力曲线逐点构造拟合，利用构造的毛管压力曲线逐点评价储层原始含油饱和度。另外，通过岩石物理实验获得相渗数据，采用实验数据与生产数据相互结合的方法，在油水渗流理论的基础上计算油藏原始含水饱和度。其它获得储层原始含水饱和度的方法是通过井壁取心或密闭取心，利用实验测试工艺求取单个样品的饱和度，不能连续计算储层原始油藏含水饱和度。同时，通过压汞实验确定油藏原始含水饱和度是一种间接的方法，计算的精度相对较低。

本发明根据油藏的成藏特性，在构建原始含水饱和度计算模型时引入反映成藏过程中油藏驱替压力的油柱高度；其次引入能反映储层排驱压力大小的孔隙结构特征参数最后引入反映储层储集空间大小的孔隙度。原始含水饱和度计算模型直接利用密闭取心实验测量的含水饱和度直接进行刻度，使模型的计算精度较高。

技术实现要素：

为解决现有技术中存在的上述缺陷，本发明的目的在于提供一种基于密闭取心的原始油藏含水饱和度计算方法，在密闭取心岩石物理实验分析的基础上，通过测井资料连续求取储层原始油藏含水饱和度，能够指导测井解释，识别储层流体性质。

本发明是通过下述技术方案来实现的。

一种基于密闭取心的原始油藏含水饱和度计算方法，包括以下步骤：

步骤101，首先利用区块已有的密闭取心井，由实验分析得到岩心的孔隙度、渗透率和饱和度参数；

步骤102，结合区块的测井与试油资料确定研究区油藏原始油水界面；

步骤103，构建油藏原始含水饱和度计算模型，结合岩石物理实验分析数据，通过最小二乘法确定模型中的系数；

步骤104，利用岩心实验分析的孔隙度和渗透率，结合测井曲线，通过最小二乘法建立孔隙度和渗透率计算模型；

步骤105，利用孔隙度、渗透率模型、原始含水饱和度计算模型及已确定的油藏原始油水界面，通过区块的测井资料计算油藏原始含水饱和度。

进一步地，所述步骤101，利用区块已有的密闭取心井，由岩心实验分析得到岩心的孔隙度、渗透率和饱和度参数。

进一步地，所述的岩心实验分析按照《岩心分析方法sy/t5336-2006》标准规定的流程进行。

进一步地，所述步骤103构建油藏原始含水饱和度计算模型如下：

式中：swy为储层原始油藏含水饱和度；dep为储层当前测井深度；owl为油藏原始油水界面；k为储层渗透率；por为储层孔隙度；a、b、c、d为模型公式中的系数。

进一步地，模型公式(1)中的系数a、b、c、d利用密闭取心岩心分析的孔隙度、渗透率和饱和度，通过最小二乘法拟合获得。

进一步地，所述步骤104，建立的孔隙度、渗透率计算模型如下：

por＝a1×den+b1(2)

式中：k为储层渗透率；por为储层孔隙度；den为密度测井值；δgr为自然伽马相对值；gr为自然伽马测井值；grmin为纯砂岩自然伽马测井值；grmax为纯泥岩段自然伽马测井值；a1、b1、a2、b2、c2为模型公式中的系数，通过密闭取心岩心实验分析数据与测井数据，采用最小二乘法拟合获得。

本发明的有益效果在于，实现了在密闭取心岩心实验分析饱和度的基础上，建立原始油藏含水饱和度计算模型，通过测井资料求取储层原始含水饱和度，提供了一种不通过储层电阻率确定饱和度的方法，能够指导测井解释，进行储层流体性质识别。

附图说明

图1为本发明提供的一种基于密闭取心原始油藏含水饱和度计算方法流程图；

图2为本发明实施例中提供的xx11井油藏原始油水界面确定图；

图3为本发明实施例中提供的原始含水饱和度计算结果验证图；

图4为本发明实施例中提供的水淹层原始含水饱和度计算成果图。

具体实施方式

下面通过对本实施例的具体实施情况做进一步详细说明，以支持本发明所要解决的技术问题，但并不作为对发明做任何限制的依据。

如图1所示，本发明实施例提供的一种基于密闭取心的原始油藏含水饱和度计算方法，包括如下步骤：

步骤101：首先利用区块已有的密闭取心井，由实验分析得到岩心的孔隙度、渗透率和饱和度参数；

步骤102：结合区块的测井与试油资料确定研究区油藏原始油水界面；

步骤103：构建油藏原始含水饱和度计算模型，结合岩石物理实验分析数据，通过最小二乘法确定模型中的系数；

步骤104：利用岩心实验分析的孔隙度和渗透率，结合测井曲线，通过最小二乘法建立孔隙度和渗透率计算模型；

步骤105：利用孔隙度、渗透率模型、原始含水饱和度计算模型及已确定的油藏原始油水界面，通过区块的测井资料计算油藏原始含水饱和度。

下面，通过对本实施例的具体实施情况做进一步详细说明，以支持本发明所要解决的技术问题，按照以下步骤进行操作：

步骤一，选取区域内进行过从上到下连续密闭取心的井xx1-6-9井，选取不同取心深度、不同粒度下的岩心按照《岩心分析方法sy/t5345-2007》标准规定，开展岩心孔隙度、渗透率、含水饱和度测量，获取不同深度变化下岩心的孔隙度、渗透率、含水饱和度，为后续的建模准备数据。

步骤二，结合区块的xx11测井与试油资料综合确定研究区油藏原始油水界面为5827m，如图2所示。

步骤三，利用步骤一中密闭取心岩心分析的孔隙度、渗透率、饱和度，通过最小二乘法拟合得到模型公式(1)中的系数a、b、c、d，确定油藏原始含水饱和度计算模型如下：

式中：swy为储层原始油藏含水饱和度；dep为储层当前测井深度；owl为油藏原始油水界面；k为储层渗透率；por为储层孔隙度；a、b、c、d为模型公式中的系数，系数a、b、c、d利用密闭取心岩心分析的孔隙度、渗透率和饱和度，通过最小二乘法拟合获得。

步骤四，利用步骤一中岩心实验分析的孔隙度、渗透率，结合测井曲线，通过最小二乘法拟合确定孔隙度、渗透率计算模型如下：

por＝71.43×den+147.86(2)

k＝10^-2.42×por^2.34/δgr^1.12(3)

式中：k为储层渗透率；por为储层孔隙度；den为密度测井值；δgr为自然伽马相对值；gr为自然伽马测井值；grmin为纯砂岩自然伽马测井值；grmax为纯泥岩段自然伽马测井值。

步骤五，通过密度测井与伽马测井，利用步骤四中式(2)、(3)、(4)确定的孔隙度、渗透率模型计算储层的孔隙度与渗透率；在储层孔隙度、渗透率计算的基础上，结合步骤二中确定原始油水界面，通过步骤三中式(1)确定的原始含水饱和度计算模型，逐点求取随深度变化的储层原始含水饱和度，如图3、图4所示。

图3为本发明提供的一口井原始油藏含水饱和度计算结果，求取的原始油藏含水饱和度与密闭取心分析结果基本一致，验证了本方法的可靠性；图4为在新井中的应用，通过原始油藏含水饱和度与目前计算含水饱和度进行对比，可以有效识别水淹层，指导测井资料解释评价。

最后所应说明的是，以上具体实施方式仅用以说明本发明的技术方案而非限制，尽管参照实例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。