CN106837271A - 一种确定泡沫油高压物性参数的计算方法 - Google Patents

Abstract

一种确定泡沫油高压物性参数的计算方法，用于表征泡沫油油藏降压生产过程中稠油捕获气体的能力，反映稠油产生泡沫油能力的大小。本发明形成了一系列确定平衡状态下泡沫油溶解气油比、压缩系数、体积系数、密度和粘度的计算公式，并通过建立的平衡状态下高压物性参数、平衡时间与非平衡状态下高压物性参数、不同静止时间的函数关系，确定泡点压力下非平衡泡沫油高压物性参数。最后，确定不同静止时间和压力下泡沫油溶解气油比、压缩系数、体积系数、密度和粘度。从而实现了快速、准确确定平衡和非平衡状态下泡沫油高压物性参数的目的。本发明对于该类稠油油藏的储量计算、开发方式选择、油气田地面工程集输设计、开发方案制定等具有重要的意义。

Description

一种确定泡沫油高压物性参数的计算方法

技术领域

本发明涉及一种确定泡沫油高压物性参数的计算方法，属于稠油油藏开发的技术领域。

背景技术

世界范围内稠油资源十分丰富，占油气资源总量的三分之一以上。我国探明和控制储量19亿吨以上，主要分布于辽河、新疆、胜利等十几个油田。此外，随着我国油气勘探开发的国际化，大量国外稠油资源亟待开发。因此，在常规油气勘探开发难度越来越大，能源矛盾日益突出的当下，大力开发稠油资源符合国家经济发展战略。

在加拿大、委内瑞拉以及中国部分稠油油藏溶解气驱过程中，出现了异于常规稠油油藏的生产特征：①油藏产出油呈现连续的泡沫状态，原油中含有大量稳定气泡。②油藏生产气油比上升速度缓慢。③油藏采收率与采油速度较高(较常规稠油油藏溶解气驱采收率高出5-25％，采油速度高出10-30倍，有的甚至高达100倍)。研究表明，泡沫油的产生是上述现象的主要原因。泡沫油的产生是由于地层压力降至泡点压力后，溶解气从稠油中析出，由于特殊稠油较高的粘滞力和压力梯度，较低的扩散速度，溶解气以小气泡的形式分散在油相中，形成泡沫油。

原油高压物性主要用于表征不同油藏压力和温度下原油体积及其相态特性的变化规律。原油高压物性研究结果广泛应用于油田开发的储量计算、开发方式选择、油气田地面工程集输设计、开发方案制定(油藏数值模拟和动态分析)等多个方面，是油田开发人员了解原油性质，优化油田生产过程的重要依据。由于泡沫油是一种具有分散气泡的非平衡体系，分散在原油中的气泡随着时间的延长和压力的降低逐渐脱离原油，这使得泡沫油压缩系数、粘度、密度等高压物性参数不仅与地层压力、温度和原油组分有关，还与时间有关。因此，除了通过常规高压物性实验方法测量平衡状态下的泡沫油高压物性参数外，还应确定不同时间非平衡状态下的泡沫油高压物性参数。目前，国内外专家学者采用非常规实验方法研究非平衡泡沫油高压物性特征。该实验过程与常规高压物性实验相似，不同之处在于压力低于泡点压力时，每级降压后不对原油体系进行搅拌，避免在人为作用下气泡聚并形成气相，从而使得溶解气分散在原油中保持泡沫油状态，并通过改变每一级降压的时间(静止时间)，通过多组重复实验，研究不同时间下(静止时间不同导致油相中溶气量不同)对泡沫油高压物性特征的影响。但是，通过上述实验过程确定泡沫油高压物性特征存在以点不足：

①由于实验样品为稠油，气体难以溶解，配制实验用地层样品时需要较长的搅拌时间。此外，由于需要改变每一级降压的时间进行多组重复实验，同时确定平衡和非平衡状态下泡沫油高压物性参数，也大大增加了实验时间；②原油高压物性实验装置包括PVT仪、配样仪、粘度计和密度计等高精度、耐高温、高压设备。因此，开展室内实验研究价格昂贵，操作复杂；③具有泡沫油特征的特殊稠油样品主要来自加拿大、委内瑞拉以及中国部分稠油油藏，十分稀少，难以大量获得。

由于上述不足的存在，只依靠室内实验无法满足现场技术人员快速、全面了解泡沫油高压物性特征，高效开发该类稠油油藏的需求。因此，通过数学计算方法，快速、准确、系统地预测平衡和非平衡状态下泡沫油高压物性参数，是解决室内实验研究不足的重要途径。目前，对于常规稠油平衡状态下高压物性参数存在多种计算公式。但是，上述公式对特殊稠油(能产生泡沫油现象)的适用性有待进一步评价，更缺少计算非平衡状态下泡沫油高压物性参数的系统计算方法。因此，如果能够形成一种确定平衡和非平衡状态下泡沫油高压物性参数的计算方法，快速、准确地确定泡沫油高压物性参数，将大大减少室内实验所产生的时间和人力成本，对于该类特殊稠油油藏的高效开发具有重要的意义。

发明内容

针对现有技术的不足，本发明提供一种确定泡沫油高压物性参数的计算方法。

本发明的具体技术方案如下：

一种确定泡沫油高压物性参数的计算方法，包括如下步骤：

步骤1.根据不同静止时间下泡沫油溶解气油比值，确定不同静止时间下的气体分散系数α；

步骤2.确定所研究特殊稠油的物性参数，包括：泡点压力下溶解气油比R_sb、油藏温度T、死油密度ρ_STO；所述特殊稠油是指能产生泡沫油现象的稠油；

步骤3.确定特殊稠油的泡点压力p_b和每一级降压后油、气完全分离，达到平衡的时间t；

步骤4.确定平衡状态下泡沫油溶解气油比R_s(p)、压缩系数C_o(p)、体积系数B_o(p)、粘度μ_o(p)和密度ρ_o(p)；

步骤5.从油田开发实际需要出发，确定一个或者多个静止时间t’_i，i＝1，2，3，…，n，用于确定该静止时间下的非平衡泡沫油高压物性参数；

步骤6.重复步骤4，确定平衡状态下泡点压力p_b下泡沫油溶解气油比Rs(p_b)、压缩系数Co(p_b)、体积系数Bo(p_b)、粘度μo(p_b)和密度ρo(p_b)；

步骤7.根据建立的平衡状态下泡沫油高压物性参数、平衡时间与非平衡状态下泡沫油高压物性参数、不同静止时间的函数关系(7)～(11)，计算不同静止时间非平衡状态下拟泡点压力p_sb下泡沫油溶解气油比Rs’(p_sb)、压缩系数Co’(p_sb)、体积系数Bo’(p_sb)、密度ρo’(p_sb)和粘度μo’(p_sb)；

步骤8.给出一系列需要求取的压力点p，根据建立的公式(12)～(16)，计算不同静止时间和压力下非平衡泡沫油溶解气油比Rs(p)、压缩系数Co(p)、体积系数Bo(p)、粘度μo(p)和密度ρo(p)：

R_s'(p)＝R_s(p)(p≥p_b,p＝p_sc) (12)

R_s'(p)＝R_s(p_b) (p_psb≤p＜p_b)

C_o'(p)＝C_o(p) (p≥p_b,p＝p_sc) (13)

B_o'(p)＝B_o(p) (p≥p_b,p＝p_sc) (14)

ρ_o'(p)＝ρ_o(p)(p≥p_b,p＝p_sc) (15)

μ_o'(p)＝μ_o(p) (p≥p_b,p＝p_sc) (16)

根据本发明优选的，所述不同静止时间下的气体分散系数α的计算方法见公式(1)：

α＝(x-y)/x (1)

x为某一压力下特殊稠油中溶解气总体积，m³；由于该体积为平衡状态下特殊稠油逸出溶解气的最大值，因此该参数可由常规差异分离实验所确定的溶解气油比确定。常规差异分离实验是用于研究油藏降压过程中平衡状态下原油的溶解气油比、密度和粘度等高压物性参数的一种常用实验方法；y为非平衡状态下泡沫油中溶解气体积，m³。由于该体积为非平衡状态下泡沫油中溶解的气量，因此该参数可由非常规差异分离实验所确定的非平衡溶解气油比确定。常规差异分离实验过程与常规差异分离实验相似，不同之处在于压力低于泡点压力时，每级降压后不对稠油进行搅拌，避免在人为作用下气泡聚并形成气相，导致泡沫油现象消失。

根据本发明优选的，在所述步骤4中，根据公式(2)～(6)确定平衡状态下泡沫油溶解气油比R_s(p)、压缩系数C_o(p)、体积系数B_o(p)、粘度μ_o(p)和密度ρ_o(p)；

式中，p为压力,MPa；p_b为泡点压力,MPa；R_sb为泡点压力下的溶解气油比,m³/m³；

式中，T为温度,℃；γ_gs为分离器压力为0.6895MPa下，分离气体的相对密度；γ_o死油相对密度；

B_o(p)＝B_ob[1-C_o(p-p_b)]p≥p_b (4)

式中，B_ob为泡点压力下特殊稠油体积系数，m³/m³；

ρ_o(p)＝ρ_STO-1.22048×10^-3R_sγ_gs/B_o p≤p_b (5)

式中，ρ_STO为死油密度，g/cm³；

式中，μ_oD为死油粘度，mPa·s；μ_ob为泡点压力下特殊稠油粘度，mPa·s。

本发明的优势在于：

1.本发明利用气体分散系数计算方法可以表征泡沫油油藏降压生产过程中稠油捕获气体的能力，反映稠油产生泡沫油能力的大小。

2.本发明可以准确地确定不同静止时间和压力下平衡和非平衡泡沫油溶解气油比、压缩系数、体积系数、密度和粘度，明确泡沫油高压物性特征，大大降低了室内实验所产生的人力和时间成本，解决了只能依靠室内实验确定泡沫油高压物性的问题，满足了现场技术人员全面了解泡沫油高压物性的需求。

3.本发明所述计算方法所需稠油基础参数较少，计算过程简单、易懂，可以满足现场技术人员在参数较少情况下快速计算的需求。

4.本发明所述计算方法是在中国、加拿大、委内瑞拉和阿尔巴尼亚等多个国家多个特殊稠油样品实验研究基础上形成的，因此，该计算方法具有较强的普遍适用性。

5.本发明所述计算方法得到的泡沫油高压物性参数可以用于特殊稠油油藏的储量计算、开发方式选择、油气田地面工程集输设计、开发方案制定(油藏数值模拟和动态分析)等多个方面，具有较广的使用范围。

附图说明

图1是本发明所述泡沫油高压物性参数计算流程图；

图2a是本发明实施例1中平衡状态下泡沫油溶解气油比计算与实验结果对比图；

图2b是本发明实施例1中平衡状态下泡沫油压缩系数计算与实验结果对比图；

图2c是本发明实施例1中平衡状态下泡沫油体积系数计算与实验结果对比图；

图2d是本发明实施例1中平衡状态下泡沫油密度计算与实验结果对比图；

图3a是本发明实施例1中非平衡状态下泡沫油溶解气油比计算与实验结果对比图；

图3b是本发明实施例1中非平衡状态下泡沫油体积系数计算与实验结果对比图；

图3c是本发明实施例1中非平衡状态下泡沫油密度计算与实验结果对比图；

图3d是本发明实施例1中非平衡状态下泡沫油粘度计算与实验结果对比图；

图4是本发明实施例2中不同静止时间下气体分散系数图。

具体实施方式

下面结合实施例和说明书附图对本发明做详细的说明，但不限于此。实施例中共选择了表1所示10个能够形成泡沫油现象的特殊稠油样品。通过对比本发明所述方法的计算结果与常规和非常规差异分离实验所得平衡和非平衡泡沫油高压物性实验结果，全面验证本发明所述方法计算结果的可靠性。此外，由于上述样品来自委内瑞拉、加拿大和阿尔巴尼亚的多个稠油油藏，因此可以验证本发明所述方法是否具有广泛的适用性。

表1计算方法验证用实验数据

实施例1、

一种确定泡沫油高压物性参数的计算方法，计算流程图见图1所示，计算步骤如下：

步骤1.根据不同静止时间下泡沫油溶解气油比值，确定不同静止时间下的气体分散系数α；

步骤2.确定所研究特殊稠油的物性参数，包括：泡点压力下溶解气油比R_sb、油藏温度T、死油密度ρ_STO；

本实施例1中，确定所研究的特殊稠油为稠油样品9，其泡点压力下溶解气油比R_sb为31m³/m³，油藏温度T为43℃，死油密度ρ_STO为0.998g/cm³。

步骤3.确定特殊稠油的泡点压力p_b和每一级降压后油、气完全分离，达到平衡的时间t；

本实施例中，特殊稠油样品9的泡点压力p_b为17MPa，平衡时间t为900h。

步骤4.根据公式(2)～(6)确定平衡状态下泡沫油溶解气油比R_s(p)、压缩系数C_o(p)、体积系数B_o(p)、粘度μ_o(p)和密度ρ_o(p)；

式中，p为压力,MPa；p_b为泡点压力,MPa；R_sb为泡点压力下的溶解气油比,m³/m³；

式中，T为温度,℃；γ_gs为分离器压力为0.6895MPa下，分离气体的相对密度；γ_o死油相对密度；

B_o(p)＝B_ob[1-C_o(p-p_b)]p≥p_b (4)

式中，B_ob为泡点压力下特殊稠油体积系数，m³/m³；

ρ_o(p)＝ρ_STO-1.22048×10^-3R_sγ_gs/B_o p≤p_b (5)

式中，ρ_STO为死油密度，g/cm³；

式中，μ_oD为死油粘度，mPa·s；μ_ob为泡点压力下特殊稠油粘度，mPa·s。

本实施例中，通过公式(2)～(5)确定的平衡状态下泡沫油溶解气油比、压缩系数、体积系数和密度如图2a-2d所示。计算结果与实验结果误差表如表2所示。由图2a-2d可知，本发明所述方法计算得到的平衡状态下泡沫油溶解气油比、压缩系数、体积系数和密度值比现有其他方法所得结果更接近实验值。此外，由表2可知，通过本发明所述方法(公式(2)～(5))计算所得平衡状态下泡沫油溶解气油比、压缩系数、体积系数和密度计算误差分别为4.8％，10.20％、0.47％和0.34％，均小于其他现有计算方法。综上可知，本发明所述方法计算结果接近真实实验值，计算误差较小，可以满足现场工程应用的需要。

表2稠油样品9平衡状态下高压物性参数计算结果与实验结果误差表

步骤5.从油田开发实际需要出发，确定一个或者多个静止时间t’_i，i＝1，2，3，…，n，用于确定该静止时间下的非平衡泡沫油高压物性参数。

本实施例中，确定3个静止时间，分别为2小时、24小时和168小时。

步骤6.重复步骤4，确定平衡状态下泡点压力p_b下泡沫油溶解气油比Rs(p_b)、压缩系数Co(p_b)、体积系数Bo(p_b)、粘度μo(p_b)和密度ρo(p_b)。

本实施例中，通过公式(2)、(4)～(6)计算平衡状态泡点压力p_b下泡沫油溶解气油比Rs(p_b)、体积系数Bo(p_b)、粘度μo(p_b)和密度ρo(p_b)分别为30.0679m³/m³、1.0705m³/m³、1150mPa·s和0.9746g/cm³。

步骤7.根据建立的平衡状态下泡沫油高压物性参数、平衡时间与非平衡状态下泡沫油高压物性参数、不同静止时间的函数关系(7)～(11)，计算不同静止时间非平衡状态下拟泡点压力p_sb下泡沫油溶解气油比Rs’(p_sb)、压缩系数Co’(p_sb)、体积系数Bo’(p_sb)、密度ρo’(p_sb)和粘度μo’(p_sb)；。

本实施例中，通过公式(7)、(9)～(11)，确定不同静止时间非平衡状态下拟泡点压力p_sb下泡沫油溶解气油比Rs’(p_sb)、体积系数Bo’(p_sb)、密度ρo’(p_sb)和粘度μo’(p_sb)如表3所示：

表3稠油样品9拟泡点压力下泡沫油高压物性参数表

步骤8.给出一系列需要求取的压力点p，根据以下公式(12)～(16)，计算不同静止时间下非平衡泡沫油溶解气油比Rs(p)、压缩系数Co(p)、体积系数Bo(p)、粘度μo(p)和密度ρo(p)。

R_s'(p)＝R_s(p) (p≥p_b,p＝p_sc) (12)

R_s'(p)＝R_s(p_b) (p_psb≤p＜p_b)

C_o'(p)＝C_o(p) (p≥p_b,p＝p_sc) (13)

B_o'(p)＝B_o(p) (p≥p_b,p＝p_sc) (14)

ρ_o'(p)＝ρ_o(p) (p≥p_b,p＝p_sc) (15)

μ_o'(p)＝μ_o(p) (p≥p_b,p＝p_sc) (16)

本实施例中，根据以下公式(12)、(14)～(16)，计算不同静止时间下非平衡泡沫油溶解气油比Rs(p)、体积系数Bo(p)、粘度μo(p)和密度ρo(p)如图3a-3d所示。稠油样品9非平衡状态下高压物性参数计算与实验结果误差见表4。

由图3a-3d可知，本发明所述方法计算得到的非平衡状态下泡沫油溶解气油比、体积系数、密度和粘度值接近实验值。此外，由表4可知，通过本发明所述方法(公式(12)、(14)～(16))计算所得非平衡状态下泡沫油溶解气油比、体积系数、密度和粘度平均误差分别为7.45％，0.65％，0.58％和8.75％，该计算方法与实验计算误差小于10％。综上可知，本发明所述方法计算结果接近真实实验值，计算误差较小，从而可以满足快速、全面的得到不同静止时间和压力下非平衡状态下泡沫油高压物性参数的现场工程应用需要。

表4稠油样品9非平衡状态下高压物性参数计算与实验结果误差表

通过上述步骤，计算表1所示其他9个能够形成泡沫油现象的特殊稠油样品平衡和非平衡状态下的高压物性参数。计算结果如表5和表6所示。

表5平衡状态下泡沫油高压物性参数计算结果与实验结果误差表

表6非平衡状态下泡沫油高压物性参数计算结果与实验结果误差表

由表5和表6可知，通过本发明所述计算方法得到的10个特殊稠油样品平衡状态下溶解气油比、压缩系数、体积系数、密度和粘度平均误差分别为6.55％，9.05％，0.84％、0.46％和4.5％。非平衡状态下溶解气油比、压缩系数、体积系数、密度和粘度平均误差分别为7.07％，16.07％，0.73％、0.63％和10.68％。由此可知，该方法对于平衡和非平衡状态下泡沫油高压物性计算误差较小，对不同地区的多个样品具有广泛的适用性，可以满足现场工程应用的需要。

实施例2、

如实施例1所述的一种确定泡沫油高压物性参数的计算方法，其区别在于，通过该方法计算稠油样品1气体分散系数α。根据室内实验得到的不同静止时间下泡沫油溶解气油比值，确定不同静止时间下的气体分散系数α，计算方法见公式(1)：

α＝(x-y)/x (1)

x为某一压力下特殊稠油中溶解气总体积，m³；y为非平衡状态下泡沫油中溶解气体积，m³。

本实施例中，根据公式(1)计算所得不同静止时间下气体分散系数α见图4。由图4可知，不论何种静止时间下，随着压力的降低自由气分散系数先升高后降低。泡点压力之上，分散系数为0，即不存在分散气体(泡沫油)，当压力在泡点压力与拟泡点压力之间时，分散系数逐渐增加，拟泡点压力时达到最大值。例如，静止时间为2h时，气体分散系数最大值为0.72，表明该油藏原油具有较强的泡沫油特性，而对于常规原油或者轻质原油，气体分散系数为0或者低值。当压力在拟泡点压力与大气压力之间时，气体分散系数逐渐降低，表明分散气体逐渐脱离原油，原油的泡沫油现象逐渐消失。此外，静止时间对气体分散系数影响较大，当压力为2MPa时，静止时间为1h的气体分散系数为0.64，而静止时间为12h的气体分散系数为0.35，即随着静止时间的减小，气体分散系数逐渐增加。

Claims (3)

1.一种确定泡沫油高压物性参数的计算方法，其特征在于，所述计算方法包括如下步骤：

步骤1.根据不同静止时间下泡沫油溶解气油比值，确定不同静止时间下的气体分散系数α；

步骤2.确定所研究特殊稠油的物性参数，包括：泡点压力下溶解气油比R_sb、油藏温度T、死油密度ρ_STO；

步骤3.确定特殊稠油的泡点压力p_b和每一级降压后油、气完全分离，达到平衡的时间t；

步骤4.确定平衡状态下泡沫油溶解气油比R_s(p)、压缩系数C_o(p)、体积系数B_o(p)、粘度μ_o(p)和密度ρ_o(p)；

步骤5.从油田开发实际需要出发，确定一个或者多个静止时间t’_i，i＝1，2，3，…，n，用于确定该静止时间下的非平衡泡沫油高压物性参数；

步骤6.重复步骤4，确定平衡状态下泡点压力p_b下泡沫油溶解气油比Rs(p_b)、压缩系数Co(p_b)、体积系数Bo(p_b)、粘度μo(p_b)和密度ρo(p_b)；

步骤7.根据建立的平衡状态下泡沫油高压物性参数、平衡时间与非平衡状态下泡沫油高压物性参数、不同静止时间的函数关系(7)～(11)，计算不同静止时间非平衡状态下拟泡点压力p_sb下泡沫油溶解气油比Rs’(p_sb)、压缩系数Co’(p_sb)、体积系数Bo’(p_sb)、密度ρo’(p_sb)和粘度μo’(p_sb)；

$\frac{{R_s}^{\′}\left(p_{psb}\right)}{R_s\left(p_b\right)}=0.6901\frac{{t^{\′}}_i}{t}+0.7232---\left(7\right)$

$\frac{{C_o}^{\′}\left(p_{psb}\right)}{C_o\left(p_b\right)}=-11.909\frac{{t^{\′}}_i}{t}+5.1714---\left(8\right)$

$\frac{{B_o}^{\′}\left(p_{psb}\right)}{B_o\left(p_b\right)}=-0.2622\frac{{t^{\′}}_i}{t}+1.0613---\left(9\right)$

$\frac{{{\mathrm{\ρ}}_o}^{\′}\left(p_{psb}\right)}{{\mathrm{\ρ}}_o\left(p_b\right)}=0.2353\frac{{t^{\′}}_i}{t}+0.9451---\left(10\right)$

$\frac{{{\mathrm{\μ}}_o}^{\′}\left(p_{psb}\right)}{{\mathrm{\μ}}_o\left(p_b\right)}=-0.8945\frac{{t^{\′}}_i}{t}+1.4002---\left(11\right)$

步骤8.给出一系列需要求取的压力点p，根据建立的公式(12)～(16)，计算不同静止时间和压力下非平衡泡沫油溶解气油比Rs(p)、压缩系数Co(p)、体积系数Bo(p)、粘度μo(p)和密度ρo(p)：

R_s'(p)＝R_s(p)(p≥p_b,p＝p_sc) (12)

R_s'(p)＝R_s(p_b)(p_psb≤p＜p_b)

${R_s}^{\&prime;}\left(p\right)=R_s\left(p_b\right)-\frac{R_s\left(p_b\right)-R_s\left(p_{sc}\right)}{p_{psb}-p_{sc}}(p-p_{sc}),(p_{sc}<p<p_{psb})$

C_o'(p)＝C_o(p)(p≥p_b,p＝p_sc) (13)

${C_o}^{\&prime;}\left(p\right)=C_o\left(p\right)+\frac{\&lsqb;{C_o}^{\&prime;}\left(p_{psb}\right)-C_0\left(p_{psb}\right)\&rsqb;(p_b-p)}{p_b-p_{psb}},(p_{psb}\&le;p<p_b)$

${C_o}^{\&prime;}\left(p\right)=C_o\left(p\right)+\frac{{C_o}^{\&prime;}\left(p_{psb}\right)-C_o\left(p_{psb}\right)}{p_{psb}-p_{sc}}(p-p_{sc}),(p_{sc}<p<p_{psb})$

B_o'(p)＝B_o(p)(p≥p_b,p＝p_sc) (14)

ρ_o'(p)＝ρ_o(p)(p≥p_b,p＝p_sc) (15)

μ_o'(p)＝μ_o(p)(p≥p_b,p＝p_sc) (16)

${{\mathrm{\&mu;}}_o}^{\&prime;}\left(p\right)={\mathrm{\&mu;}}_o\left(p\right)+\frac{\&lsqb;{\mathrm{\&mu;}}_o\left(p_{psb}\right)-{\mathrm{\&mu;}}_0\left(p_{psb}\right)\&rsqb;(p_b-p)}{p_b-p_{psb}},(p_{psb}\&le;p<p_b)$

${{\mathrm{\&mu;}}_o}^{\&prime;}\left(p\right)={\mathrm{\&mu;}}_o\left(p\right)+\frac{{{\mathrm{\&mu;}}_o}^{\&prime;}\left(p_{psb}\right)-{\mathrm{\&mu;}}_o\left(p_{psb}\right)}{p_{psb}-p_{sc}}(p-p_{sc}),(p_{sc}<p<p_{psb}).$

2.根据权利要求1所述的一种确定泡沫油高压物性参数的计算方法，其特征在于，所述不同静止时间下的气体分散系数α的计算方法见公式(1)：

α＝(x-y)/x (1)

x为某一压力下特殊稠油中溶解气总体积，m³；y为非平衡状态下泡沫油中溶解气体积，m³。

3.根据权利要求1所述的一种确定泡沫油高压物性参数的计算方法，其特征在于，在所述步骤4中，根据公式(2)～(6)确定平衡状态下泡沫油溶解气油比R_s(p)、压缩系数C_o(p)、体积系数B_o(p)、粘度μ_o(p)和密度ρ_o(p)；

$R_s\left(p\right)=\frac{R_{sb}\&times;{(p/p_b)}^{0.83}}{1.031}---\left(2\right)$

式中，p为压力,MPa；p_b为泡点压力,MPa；R_sb为泡点压力下的溶解气油比,m³/m³；

$C_o\left(p\right)=\frac{-2000+4.968R_s+1.39\&times;{10}^{-7}T-429.56{\mathrm{\&gamma;}}_{gs}+\frac{2223.141}{{\mathrm{\&gamma;}}_o}}{10231.89p}---\left(3\right)$

式中，T为温度,℃；γ_gs为分离器压力为0.6895MPa下，分离气体的相对密度；γ_o死油相对密度；

B_o(p)＝B_ob[1-C_o(p-p_b)]p≥p_b (4)

$B_o\left(p\right)=1+2.6261\&times;{10}^{-3}{R}_s+\frac{1}{{\mathrm{\&gamma;}}_{gs}}(6.447\&times;{10}^{-2}-3.7441\&times;{10}^{-4}{R}_s)(6.4286\&times;{10}^{-2}T-1)\&times;(\frac{1.076}{{\mathrm{\&gamma;}}_o}-1),p<p_b$

式中，B_ob为泡点压力下特殊稠油体积系数，m³/m³；

ρ_o(p)＝ρ_STO-1.22048×10^-3R_sγ_gs/B_o p≤p_b (5)

${\mathrm{\&rho;}}_o\left(p\right)={\mathrm{\&rho;}}_{ob}\&times;e^{\&lsqb;C_o(p-p_b)\&rsqb;},p>p_b$

式中，ρ_STO为死油密度，g/cm³；

${\mathrm{\&mu;}}_o\left(p\right)={(5.615\&times;{10}^{-2}{R}_s+1)}^{-0.515}{{\mathrm{\&mu;}}_{oD}}^{{(3.7433\&times;{10}^{-2}{R}_s+1)}^{-0.338}},p<p_b---\left(6\right)$

${\mathrm{\&mu;}}_o\left(p\right)={\mathrm{\&mu;}}_{ob}{\left(\frac{p}{p_b}\right)}^{956.4295p^{1.187}\exp \&lsqb;-(11.513+1.3024\&times;{10}^{-2}p)\&rsqb;},p\&GreaterEqual;p_b$

${\mathrm{\&mu;}}_{ob}=\frac{0.09\sqrt{{\mathrm{\&gamma;}}_g}}{\sqrt[3]{\frac{R_s}{0.178}}{\left(\frac{1.8T+32+459.67}{459.67}\right)}^{4.5}{\left(1-{\mathrm{\&gamma;}}_o\right)}^3},p=p_b$

式中，μ_oD为死油粘度，mPa·s；μ_ob为泡点压力下特殊稠油粘度，mPa·s。