CN103967482A - 一种多相流体的饱和度解释方法 - Google Patents

Abstract

本发明公开了一种多相流体的饱和度解释方法，按流体的流动特征对流体的体积系数进行加权，将多相流体的等效粘度与多相流体的等效流动体积系数以及等效累积体积系数进行有效分离，不仅更符合流动本质，而且提高了计算精度，同时在确定井底计算压力的过程中引入多相流体的等效累积体积系数，使等效粘度的变化与绝对渗透率的变化特征异化，从而可同时解释出绝对渗透率，扩大数据解释的范围。

Description

一种多相流体的饱和度解释方法

技术领域

本发明属于油藏开发技术领域，尤其涉及一种多相流体的饱和度解释方法。

背景技术

在石油工业领域，流体的饱和度是指储层岩石孔隙中某种流体所占的体积百分数，它表示了孔隙空间为某种流体所占据的程度，岩石中由多相流体充满其孔隙，这几相流体的饱和度之和为1。

申请号为201110209177.9的发明专利申请公开了《一种多相流体的饱和度解释方法》，无需借助除压力检测装置之外的其他专用检测设备，也无需使用放射性物质，可以对油气井更大范围内多相流体的饱和度进行解释。

在石油工业中，多相流体的饱和度是储量计算最重要的参数，若地层中各相流体的饱和度发生误差时，会错失油层。因此，如何提高多相流体饱和度的解释精度，是本领域技术人员亟待解决的问题。

发明内容

有鉴于此，本发明的目的在于提供一种多相流体的饱和度解释方法，相对现有技术可以提高多相流体饱和度的解释精度。

为实现上述目的，本发明提供如下技术方案：

一种多相流体的饱和度解释方法，包括：

确定地层中一个区域内多相流体的等效粘度，包括：设定油藏参数，所述油藏参数包括油气井的地质模型、井筒类型、油藏边界、地层中多个区域内多相流体的预设等效粘度、原始地层压力、绝对渗透率、所述油气井的预设井储常数和预设井筒表皮；确定所述预设等效粘度为等效粘度因数，利用所述等效粘度因数解释所述多相流体中各相流体的饱和度；确定所述各相流体的体积系数和粘度；利用所述各相流体的饱和度、体积系数、粘度以及与各相流体的饱和度对应的相对渗透率，计算多相流体的等效流动体积系数和等效累积体积系数；利用所述油藏参数、等效流动体积系数和等效累积体积系数进行渗流方程求解，确定井底计算压力；获取井底实测压力，对所述井底计算压力和井底实测压力进行拟合；判断所述拟合的精度是否满足预设要求，若是，确定当前多个区域内多相流体的预设等效粘度为多相流体的等效粘度，否则，重新设定油藏参数；

确定所述等效粘度为等效粘度因数，利用所述等效粘度因数解释所述多相流体中各相流体的饱和度。

优选的，

在上述方法中，所述利用所述等效粘度因数解释所述多相流体中各相流体的饱和度，包括：

设定所述各相流体的预设饱和度，所述各相流体的预设饱和度之和为1；

分别依据所述各相流体的相对渗透率曲线，确定与所述各相流体的预设饱和度对应的相对渗透率；

计算所述各相流体的相对渗透率与粘度比值之和，并进行记录；

按照预设策略调整所述各相流体的预设饱和度；

当各相流体调整后的预设饱和度位于预设范围时，返回执行确定与所述各相流体的预设饱和度对应的相对渗透率的步骤；

当各相流体调整后的预设饱和度超出预设范围时，分别比较记录的各比值之和与所述等效粘度因数的倒数是否相同，当仅有一个比值之和与所述等效粘度因数的倒数相同时，确定与该比值之和所对应的预设饱和度为地层中该区域内各相流体的饱和度，当有两个比值之和与所述等效粘度因数的倒数相同时，分别确定与所述两个比值之和所对应的预设饱和度为各相流体的备用饱和度，并按照预设规则确定所述备用饱和度中的一组为各相流体的饱和度。

优选的，

在上述方法中，所述利用所述等效粘度因数解释所述多相流体中各相流体的饱和度，包括：

设定所述各相流体的预设饱和度，所述各相流体的预设饱和度之和为1；

分别依据所述各相流体的相对渗透率曲线，确定与所述各相流体的预设饱和度对应的相对渗透率；

计算所述各相流体的相对渗透率与粘度比值之和；

判断所述比值之和与所述等效粘度因数的倒数是否相同；

当所述比值之和与所述等效粘度因数的倒数相同时，确定与该比值之和对应的预设饱和度为各相流体的备用饱和度并记录，若当前记录有两组备用饱和度，则按照预设规则确定所述备用饱和度中的一组为各相流体的饱和度，若当前记录有一组备用饱和度，则转向执行按照预设策略调整所述各相流体的预设饱和度的步骤；

当所述比值之和与所述等效粘度因数的倒数不同时，执行按照预设策略调整所述各相流体的预设饱和度的步骤；

判断调整后的各相流体的预设饱和度是否处于预设范围内，若是，则返回执行确定与所述各相流体的预设饱和度对应的相对渗透率的步骤，若否，则确定当前记录的备用饱和度为各相流体的饱和度。

优选的，

在上述方法中，所述利用所述各相流体的饱和度、体积系数、粘度以及与各相流体的饱和度对应的相对渗透率，计算多相流体的等效流动体积系数和等效累积体积系数，包括：

分别计算各相流体的流动系数，每相流体的流动系数为该相流体的相对渗透率和粘度的比值；

确定各相流体的流动系数之和作为总的流动系数；

分别计算各相流体的分流系数，每相流体的分流系数为该相流体的流动系数和所述总的流动系数的比值；

计算所述多相流体的等效流动体积系数：当所述多相流体为油水两相流体或气水两相流体时，所述多相流体的等效流动体积系数为各相流体的分流系数与体积系数之比的和值的倒数；当所述多相流体为油气两相流体或油气水三相流体时，所述多相流体的等效流动体积系数为各相流体的分流系数与体积系数之比、以及油相流体的分流系数和溶解气油比的乘积与油相流体的体积系数之比的和值的倒数；；

计算所述多相流体的等效累积体积系数：当所述多相流体为油水两相流体或气水两相流体时，所述多相流体的等效累积体积系数为各相流体的饱和度与体积系数之比的和值的倒数；当所述多相流体为油气两相流体或油气水三相流体时，所述多相流体的等效累积体积系数为各相流体的饱和度与体积系数之比、以及油相流体的饱和度和溶解气油比的乘积与油相流体的体积系数之比的和值的倒数。

优选的，

在上述方法中，按照预设规则确定所述备用饱和度中的一组为各相流体的饱和度，包括：

分别确定两个所述备用饱和度对应的理论含水率；

获取实测含水率；

比对两个所述理论含水率与所述实测含水率，确定两个理论含水率中与所述实测含水率之间差值较小的理论含水率；

确定与具有较小差值的理论含水率对应的备用饱和度为各相流体的饱和度。

优选的，

在上述方法中，按照预设规则确定所述备用饱和度中的一组为各相流体的饱和度，包括：

获取所述多相流体中任意一相流体的最大饱和度和最小饱和度；

将两个所述备用饱和度中该相流体的饱和度分别与所述最大饱和度、最小饱和度进行比较；

将位于所述最大饱和度和最小饱和度之间的该相流体饱和度所在的备用饱和度作为所述各相流体的饱和度。

优选的，

在上述方法中，所述按照预设策略调整所述各相流体的预设饱和度，包括：

当所述多相流体为两相流体时，以第一预设间隔增大所述两相流体中一相流体的预设饱和度，相应的以所述第一预设间隔减小另一相流体的预设饱和度；

当所述多相流体为三相流体时，以第二预设间隔增大所述三相流体中一相流体的预设饱和度，以第三预设间隔减小另一相流体的预设饱和度。

优选的，

在上述方法中，所述拟合包括压力史拟合、压力降落或恢复拟合、以及压力导数拟合。

由此可见，本发明公开的多相流体的饱和度解释方法，按流体的流动特征对流体的体积系数进行加权，将多相流体的等效粘度与多相流体的等效流动体积系数以及等效累积体积系数进行有效分离，不仅更符合流动本质，而且提高了计算精度。在确定井底计算压力的过程中引入多相流体的等效累积体积系数，使等效粘度的变化与绝对渗透率的变化特征异化，从而可同时解释出绝对渗透率，扩大数据解释的范围。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明公开的一种多相流体的饱和度解释方法的流程图；

图2为本发明公开的一种利用等效粘度因数解释各相流体的饱和度的方法的流程图；

图3为本发明公开的另一种利用等效粘度因数解释各相流体的饱和度的方法的流程图；

图4为本发明公开的一种计算多相流体的等效流动体积系数和等效累积体积系数的方法的流程图；

图5为本发明公开的一种按照预设规则确定备用饱和度中的一组为各相流体的饱和度的方法的流程图；

图6为本发明公开的另一种按照预设规则确定备用饱和度中的一组为各相流体的饱和度的方法的流程图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明公开了一种多相流体的饱和度解释方法，具有较高的解释精度。

参见图1，图1为本发明公开的一种多相流体的饱和度解释方法的流程图。包括：

步骤S1：设定油藏参数。

该油藏参数包括油气井的地质模型、井筒类型、油藏边界、地层中多个区域内多相流体的预设等效粘度、原始地层压力、绝对渗透率、油气井的预设井储常数和预设井筒表皮。

步骤S2：确定该预设等效粘度为等效粘度因数，利用等效粘度因数解释多相流体中各相流体的饱和度。

本发明公开了多种利用等效粘度因数解释各相流体的饱和度的方法。参见图2，图2示出了一种利用等效粘度因数解释各相流体的饱和度的方法，包括：

步骤S21：设定各相流体的预设饱和度。

油气井中的多相流体可能为油气两相流体、油水两相流体、气水两相流体或者油气水三相流体，多相流体中各相流体的饱和度之和为1，因此设定的各相流体的预设饱和度之和为1。

步骤S22：分别依据各相流体的相对渗透率曲线，确定与各相流体的预设饱和度对应的相对渗透率。

在试井过程中，可以确定各相流体的相对渗透率曲线，相对渗透率曲线反映了流体的饱和度和渗透率之间的关系。分别在各相流体的预设饱和度在该相流体的相对渗透率曲线上进行插值，就可以确定与各相流体的预设饱和度对应的相对渗透率。

步骤S23：计算各相流体的相对渗透率与流体粘度比值之和，并记录。

步骤S24：按照预设策略调整各相流体的预设饱和度。

油气井中的多相流体可能为两相流体，也可能为三相流体。当油气井中的多相流体为两相流体（即油气两相流体、油水两相流体、气水两相流体）时，以第一预设间隔增大两相流体中一相流体的预设饱和度，相应的以第一预设间隔减小另一相流体的预设饱和度；当油气井中的多相流体为三相流体（即油气水三相流体）时，以第二预设间隔增大三相流体中一相流体的预设饱和度，以第三预设间隔减小另一相流体的预设饱和度，则第三相流体的预设饱和度相应发生变化。

以油水两相流体为例，油相流体的初始饱和度为S_o1、残余饱和度为S_o2，相应的，水相流体的最大饱和度S_w1为1-S_o2，最小饱和度S_w2为1-S_o1，随着油气井的开发，油相流体的饱和度会逐渐减小，水相流体的饱和度会逐渐增加。如果在执行步骤S21时，设定油相流体的预设饱和度为S_o1，设定水相流体的预设饱和度为1-S_o1，那么在步骤S24中可以第一预设间隔增大水相流体的预设饱和度，以第一预设间隔减小油相流体的预设饱和度；如果在执行步骤S21时，设定油相流体的预设饱和度为S_o2，设定水相流体的预设饱和度为1-S_o2，那么在步骤S24中可以第一预设间隔增大油相流体的预设饱和度，以第一预设间隔减小水相流体的预设饱和度。

步骤S25：判断各相流体调整后的预设饱和度是否位于预设范围内，若是，则转向步骤S22，若否，则转向步骤S26。

地层的多相流体中各相流体的饱和度在油气井开发过程会持续发生规律性的变化，各相流体均存在一个最大饱和度和最小饱和度，地层中多相流体中各相流体的饱和度在整个开发过程中不会超出最大饱和度和最小饱和度。在步骤S25中，可以将各相流体饱和度的预设范围确定为由各相流体的最大饱和度和最小饱和度构成的区间。当按照预设策略调整后的各相流体的预设饱和度处于预设范围时，根据当前各相流体的预设饱和度确定的各比值之和可能与等效粘度因数的倒数相同，也就是说当前各相流体的预设饱和度可能是各相流体的饱和度；当按照预设策略调整后的各相流体的预设饱和度超出预设范围时，此时的预设饱和度肯定不是各相流体的饱和度，因此不需要对其进行后续操作。

步骤S26：分别比较记录的各比值之和是否与等效粘度因数的倒数相同，并确定与等效粘度因数的倒数相同的比值之和的个数。

步骤S27：判断与等效粘度因数的倒数相同的比值之和是否仅有一个，若是，转向步骤S28，若否，转向步骤S29。

步骤S28：确定与该比值之和对应的预设饱和度为各相流体的饱和度。

步骤S29：分别确定与两个比值之和对应的预设饱和度为各相流体的备用饱和度，按照预设规则确定备用饱和度中的一组为各相流体的饱和度。

由于地层中各相流体的相对渗透率曲线有一个交叉点，所以与等效粘度因数的倒数相同的各比值之和可能存在两组，相应的会得到两组饱和度值，其中一组饱和度值为假值，要进一步根据产量中的含水率、初始开发时的初始饱和度、油层的残余饱和度等信息舍弃备用饱和度中不合理的值。

在步骤S2中将多相流体的预设等效粘度作为等效粘度因数，利用等效粘度因数解释出的各相流体的饱和度是与该预设等效粘度对应的各相流体的饱和度。

参见图3，图3示出了另一种利用等效粘度因数解释各相流体的饱和度的方法，包括：

步骤S31：设定所述各相流体的预设饱和度。

需要说明的是，各相流体的预设饱和度之和为1。

步骤S32：分别依据所述各相流体的相对渗透率曲线，确定与所述各相流体的预设饱和度对应的相对渗透率。

步骤S33：计算各相流体的相对渗透率与粘度比值之和。

步骤S34：判断该比值之和与等效粘度因数的倒数是否相同，当该比值之和与等效粘度因数的倒数相同时，执行步骤S35，当该比值之和与等效粘度因数的倒数不同时，执行步骤S38。

步骤S35：确定与该比值之和对应的预设饱和度为各相流体的备用饱和度，并记录。

步骤S36：判断当前是否记录有两组备用饱和度，若是，则执行步骤S37，否则执行步骤S38。

步骤S37：按照预设规则确定所述备用饱和度中的一组为各相流体的饱和度。

步骤S38：按照预设规则调整各相流体的预设饱和度。

步骤S39：判断调整后的各相流体的预设饱和度是否处于预设范围内，若是，则返回执行步骤S32，若否，则执行步骤S310。

步骤S310：确定当前记录的备用饱和度为各相流体的饱和度。

与图2所示方法相比，图3所示的利用等效粘度因数解释各相流体的饱和度的方法中，当获取两组备用饱和度之后就不再对其他处于预设范围内的预设饱和度进行遍历，降低了运算量。

步骤S3：确定各相流体的体积系数和粘度。

各相流体的体积系数可以通过多种方式确定。例如：通过多次实验确定各相流体的体积系数列表，通过在体积系数列表中插值确定各相流体的体积系数。

或者，当多相流体为油水两相流体时，按照下述公式确定油相流体或水相流体的体积系数，

$B_l=\frac{B_l^{\mathrm{ref}}}{1+C_l(p-p^{\mathrm{ref}})}(1=i,o)$

在该公式中，B_l为油相流体或水相流体的体积系数，为油相流体或水相流体的参考体积系数，C_l为油相流体或水相流体的压缩系数，p为压力，p^ref为参考体积压力，p^ref、和C_l均为实测数据。

或者，当多相流体为油气水三相流体时，油相流体和气相流体的体积系数可以通过在体积系数列表中插值确定。并按照下述公式确定水相流体的体积系数，

$B_w=\frac{B_w^{\mathrm{ref}}}{1+C_w(p-p^{\mathrm{ref}})}$

在该公式中，B_w为水相流体的体积系数，为水相流体的参考体积系数，C_w水相流体的压缩系数，p为压力，p^ref为参考体积压力，p^ref、和C_w均为实测数据。

各相流体的的粘度可以通过多种方式确定。例如：通过多次实验确定各相流体的粘度列表，通过在粘度表中插值确定各相流体的粘度。或者，也可以利用Beggs-Robinson公式、Standing公式、GlasO公式、Khan公式或Ng-Egbogah公式计算各相流体的粘度。

步骤S4：计算多相流体的等效流动体积系数和等效累积体积系数。

利用各相流体的饱和度、体积系数、粘度以及与各相流体的饱和度对应的相对渗透率，计算多相流体的等效流动体积系数和等效累积体积系数。

步骤S5：利用油藏参数、等效流动体积系数和等效累积体积系数进行渗流方程求解，确定井底计算压力。

将渗流方程中流动项的体积系数用多相流体的等效流动体积系数替换，将渗流方程中累积项的体积系数用多相流体的等效累积体积系数替换，并利用油藏参数，即地质模型、井筒类型、油藏边界、原始地层压力、预设井储常数、预设井筒表皮和多个区域内多相流体的预设等效粘度进行渗流方程求解，确定井底计算压力。

步骤S6：获取井底实测压力，对井底计算压力和井底实测压力进行拟合。

拟合的过程包括压力史拟合、压力降落或恢复拟合、以及压力导数拟合。

步骤S7：判断拟合的精度是否满足预设要求，若满足，则转向步骤S8，若不满足，则转向步骤S9。

步骤S8：确定当前多个区域内多相流体的预设等效粘度为多相流体的等效粘度，执行步骤S10。

步骤S9：重新设定油藏参数，转向步骤S2。

根据油气井的地质模型、井筒类型、油藏边界、原始地层压力、预设井储常数、预设井筒表皮、地层中多个区域内多相流体的预设等效粘度、多相流体的等效流动体积系数和多相流体的等效累积体积系数计算井底压力，即井底计算压力。对各油藏参数进行调整，当根据各油藏参数确定的井底计算压力与实测压力的拟合精度满足预设要求时，则认为当前的各油藏参数为油气井的实测油藏参数，当前的多个区域内多相流体的预设等效粘度为多相流体的等效粘度。

步骤S10：确定该等效粘度为等效粘度因数，利用等效粘度因数解释多相流体中各相流体的饱和度。

在步骤S10中，将地层中一个区域内多相流体的等效粘度作为等效粘度因数，之后利用如图2或图3所示的方法解释多相流体中各相流体的饱和度，将解释出的各相流体的饱和度作为最终用于储量计算的最终饱和度。

以油水两相流体为例，对本发明公开的多相流体的饱和度解释方法进行说明。

确定地层中一个区域内油水两相流体的等效粘度μ_e，获取油相流体的相对渗透率曲线S_w～K_ro和粘度μ_o，获取水相流体的相对渗透率曲线S_w～k_rw和粘度μ_w；设定水相流体的预设饱和度为S_w，则油相流体的预设饱和度S_o为1-S_w；在油相流体的相对渗透率曲线上通过插值获得与S_o对应的相对渗透率K_ro，在水相流体的相对渗透率曲线上通过插值获得与S_w对应的相对渗透率K_rw；计算油相流体的相对渗透率与粘度的比值K_ro/μ_o，计算水相流体的相对渗透率与粘度的比值K_rw/μ_w，求取两比值之和并进行记录；按照预设策略调整水相流体和油相流体的饱和度，当调整后的水相流体和油相流体的饱和度仍处于预设范围时，重新执行上述确定相对渗透率的步骤以及后续步骤，当调整后的水相流体或者油相流体的饱和度超出预设范围时，获得多个两比值之和分别对多个两比值之和与等效粘度μ_e的倒数进行比较，并确定与等效粘度μ_e的倒数一致的两比值之和有几个，即判断满足这一等式的两比值之和有几个，若仅有一个，则与该两比值之和对应的预设饱和度就是地层中当前区域内油水两相流体的饱和度，若有两个，则确定与这两个比值之和对应的预设饱和度为油水两相流体的备用饱和度，进一步根据产量中的含水率、初始开发时的初始饱和度、油层的残余饱和度等信息舍弃备用饱和度中不合理的值。

在本发明上述公开的多相流体的饱和度解释方法中，按流体的流动特征对流体的体积系数进行加权，将多相流体的等效粘度与多相流体的等效流动体积系数以及等效累积体积系数进行有效分离，不仅更符合流动本质，而且提高了计算精度。在确定井底计算压力的过程中引入多相流体的等效累积体积系数，使等效粘度的变化与绝对渗透率的变化特征异化，从而可同时解释出绝对渗透率，扩大数据解释的范围。

实施中，计算多相流体的等效流动体积系数和等效累积体积系数的方法，如图4所示。包括：

步骤S41：分别计算各相流体的流动系数。

每相流体的流动系数为该相流体的相对渗透率和粘度的比值。

步骤S42：确定各相流体的流动系数之和作为总的流动系数。

步骤S43：分别计算各相流体的分流系数。

每相流体的分流系数为该相流体的流动系数和总的流动系数的比值。

步骤S44：计算多相流体的等效流动体积系数。

当多相流体为油水两相流体或气水两相流体时，该多相流体的等效流动体积系数为各相流体的分流系数与体积系数之比的和值的倒数。

当多相流体为油气两相流体或油气水三相流体时，该多相流体的等效流动体积系数为各相流体的分流系数与体积系数之比、以及油相流体的分流系数和溶解气油比的乘积与油相流体的体积系数之比的和值的倒数。也就是，计算多相流体中各相流体的分流系数和体积系数的比值，计算油相流体的分流系数和溶解气油比的乘积与油相流体的体积系数的比值，将所有比值进行加法运算得到和值，该和值的倒数为等效流动体积系数。

步骤S45：计算多相流体的等效累积体积系数。

当多相流体为油水两相流体或气水两相流体时，该多相流体的等效累积体积系数为各相流体的饱和度与体积系数之比的和值的倒数。

当多相流体为油气两相流体或油气水三相流体时，该多相流体的等效累积体积系数为各相流体的饱和度与体积系数之比、以及油相流体的饱和度和溶解气油比的乘积与油相流体的体积系数之比的和值的倒数。也就是，计算多相流体中各相流体的饱和度和体积系数的比值，计算油相流体的饱和度和溶解气油比的乘积与油相流体的体积系数的比值，将所有比值进行加法运算得到和值，该和值的倒数即为等效累积体积系数。

需要说明的是，溶解气油比为试井过程中获得的实测数据。

下面以油水两相流体为例，对确定其等效流动体积系数和等效累积体积系数的过程进行说明。

根据预设等效粘度计算油相流体的饱和度S_o和水相流体的饱和度S_w＝1-S_o。计算当前压力下的油相流体的体积系数B_o和水相流体的体积系数B_w。计算当前压力下的油相流体的粘度μ_o和水相流体的粘度μ_w。利用水相流体的相对渗透率曲线，确定水相流体在当前饱和度下的相对渗透率K_rw，利用油相流体的相对渗透率曲线，确定油相流体在当前饱和度下的相对渗透率K_ro。分别计算水相流体的流动系数和油相流体的流动系数计算总的流动系数计算水相分流系数计算油相分流系数β＝1-α。利用公式计算油水两相流体的等效流动体积系数B_e，利用公式计算油水两相流体的等效累积体积系数B_a。

当然，图3仅示出了一种计算多相流体的等效流动体积系数和等效累积体积系数的方法，实施中还可以采用其他方式。

例如：当多相流体为两相流体（即油水两相流体、油气两相流体和气水两相流体）时，首先分别计算各相流体的流动系数，确定总的流动系数，之后计算多相流体中任一相流体的分流系数（为了便于描述，记为第一分流系数），则另一相流体的分流系数为1减去第一分流系数，之后计算多相流体的等效流动体积系数和等效累积体积系数。

例如：当多相流体为油气水三相流体时，首先分别计算各相流体的流动系数，确定总的流动系数，之后计算多相流体中任两相流体的分流系数（为了便于描述，分别记为第一分流系数和第二分流系数），则第三相流体的分流系数为1减去第一分流系数和第二分流系数，之后计算多相流体的等效流动体积系数和等效累积体积系数。

本发明不再对计算多相流体的等效流动体积系数和等效累积体积系数的过程一一进行列举。

本发明公开了一种按照预设规则确定备用饱和度中的一组为各相流体的饱和度的方法，如图5所示，包括：

步骤S51：分别确定与两个备用饱和度对应的理论含水率；

步骤S52：获取实测含水率；

步骤S53：比对两个理论含水率与实测含水率，确定与实测含水率之间差值最小的理论含水率；

步骤S54：确定与具有最小差值的理论含水率对应的备用饱和度为各相流体的饱和度。

图5所示的方法中，根据井产量公式计算与每组备用饱和度对应的理论含水率，确定各个理论含水率与实测含水率之间的差值，并将最接近实测含水率的理论含水率所对应的备用饱和度作为多相流体的饱和度。

本发明公开了另一种按照预设规则确定备用饱和度中的一组为各相流体的饱和度的方法，如图6所示，包括：

步骤S61：获取多相流体中任一相流体的最大饱和度和最小饱和度；

步骤S62：将两个备用饱和度中该相流体的饱和度分别与最大饱和度和最小饱和度进行比较；

步骤S63：将位于最大饱和度和最小饱和度之间的该相流体饱和度所在的备用饱和度作为所述各相流体的饱和度。

随着油气井的开发，地层中各相流体的饱和度会发生变化，通过判断备用饱和度中某一相流体的饱和度是否位于该相流体的最大饱和度和最小饱和度之间，就可以判断该备用饱和度是否合理。

下面以油水两相流体为例进行说明。

获取油水两相流体中油相流体的最大饱和度和最小饱和度，此处油相流体的最大饱和度即为初始饱和度S_o1、最小饱和度即为残余饱和度S_o2，在油气井的开发过程中，油相流体的饱和度会逐渐减小，即残余饱和度S_o2小于初始饱和度S_o1；判断各备用饱和度中的油相流体的饱和度S_o是否位于初始饱和度S_o1（最大饱和度）和残余饱和度S_o2（最小饱和度）构成的区间，若某一备用饱和度中的油相流体的饱和度位于该区间，则该备用饱和度即为地层中各相流体的饱和度，另一备用饱和度被舍弃。

或者，获取油水两相流体中水相流体的最大饱和度1-S_o2和最小饱和度1-S_o1，在油气井的开发过程中，水相流体的饱和度会逐渐增大；判断各备用饱和度中的水相流体的饱和度S_w是否位于最大饱和度和最小饱和度构成的区间，若某一备用饱和度中的水相流体的饱和度位于该区间，则该备用饱和度即为地层中各相流体的饱和度，另一备用饱和度被舍弃。

本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。

结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块，或者二者的结合来实施。软件模块可以置于随机存储器（RAM）、内存、只读存储器（ROM）、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。

对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (8)

1.一种多相流体的饱和度解释方法，其特征在于，包括：

确定地层中一个区域内多相流体的等效粘度，包括：设定油藏参数，所述油藏参数包括油气井的地质模型、井筒类型、油藏边界、地层中多个区域内多相流体的预设等效粘度、原始地层压力、绝对渗透率、所述油气井的预设井储常数和预设井筒表皮；确定所述预设等效粘度为等效粘度因数，利用所述等效粘度因数解释所述多相流体中各相流体的饱和度；确定所述各相流体的体积系数和粘度；利用所述各相流体的饱和度、体积系数、粘度以及与各相流体的饱和度对应的相对渗透率，计算多相流体的等效流动体积系数和等效累积体积系数；利用所述油藏参数、等效流动体积系数和等效累积体积系数进行渗流方程求解，确定井底计算压力；获取井底实测压力，对所述井底计算压力和井底实测压力进行拟合；判断所述拟合的精度是否满足预设要求，若是，确定当前多个区域内多相流体的预设等效粘度为多相流体的等效粘度，否则，重新设定油藏参数；

确定所述等效粘度为等效粘度因数，利用所述等效粘度因数解释所述多相流体中各相流体的饱和度。

2.根据权利要求1所述的方法，其特征在于，所述利用所述等效粘度因数解释所述多相流体中各相流体的饱和度，包括：

设定所述各相流体的预设饱和度，所述各相流体的预设饱和度之和为1；

分别依据所述各相流体的相对渗透率曲线，确定与所述各相流体的预设饱和度对应的相对渗透率；

计算所述各相流体的相对渗透率与粘度比值之和，并进行记录；

按照预设策略调整所述各相流体的预设饱和度；

当各相流体调整后的预设饱和度位于预设范围时，返回执行确定与所述各相流体的预设饱和度对应的相对渗透率的步骤；

当各相流体调整后的预设饱和度超出预设范围时，分别比较记录的各比值之和与所述等效粘度因数的倒数是否相同，当仅有一个比值之和与所述等效粘度因数的倒数相同时，确定与该比值之和所对应的预设饱和度为地层中该区域内各相流体的饱和度，当有两个比值之和与所述等效粘度因数的倒数相同时，分别确定与所述两个比值之和所对应的预设饱和度为各相流体的备用饱和度，并按照预设规则确定所述备用饱和度中的一组为各相流体的饱和度。

3.根据权利要求1所述的方法，其特征在于，所述利用所述等效粘度因数解释所述多相流体中各相流体的饱和度，包括：

设定所述各相流体的预设饱和度，所述各相流体的预设饱和度之和为1；

分别依据所述各相流体的相对渗透率曲线，确定与所述各相流体的预设饱和度对应的相对渗透率；

计算所述各相流体的相对渗透率与粘度比值之和；

判断所述比值之和与所述等效粘度因数的倒数是否相同；

当所述比值之和与所述等效粘度因数的倒数相同时，确定与该比值之和对应的预设饱和度为各相流体的备用饱和度并记录，若当前记录有两组备用饱和度，则按照预设规则确定所述备用饱和度中的一组为各相流体的饱和度，若当前记录有一组备用饱和度，则转向执行按照预设策略调整所述各相流体的预设饱和度的步骤；

当所述比值之和与所述等效粘度因数的倒数不同时，执行按照预设策略调整所述各相流体的预设饱和度的步骤；

判断调整后的各相流体的预设饱和度是否处于预设范围内，若是，则返回执行确定与所述各相流体的预设饱和度对应的相对渗透率的步骤，若否，则确定当前记录的备用饱和度为各相流体的饱和度。

4.根据权利要求1、2或3所述的方法，其特征在于，所述利用所述各相流体的饱和度、体积系数、粘度以及与各相流体的饱和度对应的相对渗透率，计算多相流体的等效流动体积系数和等效累积体积系数，包括：

分别计算各相流体的流动系数，每相流体的流动系数为该相流体的相对渗透率和粘度的比值；

确定各相流体的流动系数之和作为总的流动系数；

分别计算各相流体的分流系数，每相流体的分流系数为该相流体的流动系数和所述总的流动系数的比值；

计算所述多相流体的等效流动体积系数：当所述多相流体为油水两相流体或气水两相流体时，所述多相流体的等效流动体积系数为各相流体的分流系数与体积系数之比的和值的倒数；当所述多相流体为油气两相流体或油气水三相流体时，所述多相流体的等效流动体积系数为各相流体的分流系数与体积系数之比、以及油相流体的分流系数和溶解气油比的乘积与油相流体的体积系数之比的和值的倒数；

计算所述多相流体的等效累积体积系数：当所述多相流体为油水两相流体或气水两相流体时，所述多相流体的等效累积体积系数为各相流体的饱和度与体积系数之比的和值的倒数；当所述多相流体为油气两相流体或油气水三相流体时，所述多相流体的等效累积体积系数为各相流体的饱和度与体积系数之比、以及油相流体的饱和度和溶解气油比的乘积与油相流体的体积系数之比的和值的倒数。

5.根据权利要求2或3所述的饱和度解释方法，其特征在于，按照预设规则确定所述备用饱和度中的一组为各相流体的饱和度，包括：

分别确定两个所述备用饱和度对应的理论含水率；

获取实测含水率；

比对两个所述理论含水率与所述实测含水率，确定两个理论含水率中与所述实测含水率之间差值较小的理论含水率；

确定与具有较小差值的理论含水率对应的备用饱和度为各相流体的饱和度。

6.根据权利要求2或3所述的饱和度解释方法，其特征在于，按照预设规则确定所述备用饱和度中的一组为各相流体的饱和度，包括：

获取所述多相流体中任意一相流体的最大饱和度和最小饱和度；

将两个所述备用饱和度中该相流体的饱和度分别与所述最大饱和度、最小饱和度进行比较；

将位于所述最大饱和度和最小饱和度之间的该相流体饱和度所在的备用饱和度作为所述各相流体的饱和度。

7.根据权利要求2或3所述的饱和度解释方法，其特征在于，所述按照预设策略调整所述各相流体的预设饱和度，包括：

当所述多相流体为两相流体时，以第一预设间隔增大所述两相流体中一相流体的预设饱和度，相应的以所述第一预设间隔减小另一相流体的预设饱和度；

当所述多相流体为三相流体时，以第二预设间隔增大所述三相流体中一相流体的预设饱和度，以第三预设间隔减小另一相流体的预设饱和度。

8.根据权利要求1所述的方法，其特征在于，所述拟合包括压力史拟合、压力降落或恢复拟合、以及压力导数拟合。