CN1678815A - 用于解释在油气生产井中测得数据的方法 - Google Patents

Abstract

一种计算在井中流动的多相流出物的至少一种相的相对体积流率的方法包括获取在特定深度处所述相横跨井某个截面的局部体积分数和/或速度的第一个步骤。该方法还包括：对所述局部体积分数和/或速度测量值进行修正，以便使其相互一致和/或与流出物流动状况一致；选择一个合适的流动模型从数学上代表流出物流动；对所述局部体积分数的测量值和所述局部体积速度测量值通过所选择的流动模型进行内插，以便获得井内在所述深度横跨所述截面的流出物的至少一种相的体积分数分布线和速度分布线。通过对井内在所述深度横跨所述截面的所述的至少一种相的体积分数分布线和速度分布线的积分计算所述的至少一种相的相对体积流率。

Description

用于解释在油气生产井中测得数据的方法

技术领域

本发明涉及一种用于解释从处于生产中的油气井(hydrocarbon well)中获得的数据的方法。

更确切地说，设计本发明的方法是为了确保在生产过程中在井中收集的数据可以得到正确的解释，以便精确地了解井内部的流动状况，也就是相对体积流率(relative volumetric flow-rate)。

背景技术

为了保证处于生产中的油气井的监测和诊断函数，必须获得一定数量的数据，主要是物理数据。这些数据实质上与在井中流动的多相流体有关(流量、不同相的各自比例、温度、压力等)。它们也能够涉及井眼本身的某些特性(成椭圆形、倾斜)。例如，这些数据可以定量地确定所有深度的流体进入区和流体退出区，以便能够在需要的时候对相应的区域重新堵塞(re-plug)或重新射孔(re-perforate)。这种监测因此可以把不想要的流体比如水的地面产量降到最低，把有价值的流体比如油的产量提到最高。

为了获得所述的数据，尤其如文献FR-A-2732068所示，传统的方法包括下面的步骤：首先借助于井轴上的纺锤形流量计(spinner)对井中的流体流动速度进行全面的测量，第二步，进行测量(也可能是局部的)使得能够确定井中特定区域的流体的各种相的比例。速度测量和持有率测量是在不同深度进行的。持有率测量是通过不同的传感器进行的，这些传感器可能是电阻率传感器、光学传感器等。

为了确定在井中流动的不同相的流体的流率(flow rate)，横跨井中某个截面的流体的流率是从由所述的现有仪器通过将总的速度测量值(包括在井中心测量的速度乘以一个通常在0.6至1之间的系数)乘以在进行测量的位置处的井截面的面积而得到的测量结果计算出来的。而后将由传感器确定的与所考虑的相有关的比例应用到所述的流率中去。

但是，众所周知，在油井中流动的流体的不同相的分布是随着井是竖直、倾斜或水平的不同而变化。由于流体不同相的密度的差异，所述的相随着井的倾斜程度的增加而变得更为循序渐进地分层。因此，在包括水、油和气的三相流体的情况下，当井是水平的或非常倾斜时，这三种相趋向于一种相在另一种相的上部流动。因此，井中相的分布(相的持有率)和每种相的速度在所述井中的一个截面上是不一致的：需要对井剖面的这些函数进行更为明确的描述以计算每种相的流率。

借助于如文献WO 01/11190所述的设置在井截面的不同的已知位置的传感器(例如局部探针和微型纺锤形流量计)将可以获得这些描述。文献WO 01/11190是基于如下认识的：流体的任何一种相的流率不等于流体的总的(平均)速度与井的截面积和与在该流动流体中所述相的体积分数的乘积，而是所考虑相的平均速度与截面积和与所述相的体积分数的乘积。

因此，为计算所有深度的相对体积流率在每一个局部纺锤形流量计和探针上收集的数据的解释成为一种非常重要的程序，以便准确估计组成井流出物的每种流体的动态。这种解释实际上需要以应用于流出物流动(emuent flow)的流体模型为基础的值的内插(interpolation of values)，所述的模型根据当前流体状况而改变。

发明内容

因此，本发明的一个目的是提供一种用于计算在井中流动的多相流出物的至少一种相的相对体积流率的方法，所述方法包括首先获取横跨处在特定深度处的井的某个截面的所述相的局部的体积分数和速度的第一步骤，该方法还包括：

对所述的局部体积分数和/或速度的测量值进行修正，以便使其互相一致和/或与流出物流动状况一致；

选择一个从数学上代表流出物流动的合适的流动模型；

通过所选择的流动模型对所述的局部体积含量测量结果和所述的局部体积速度测量结果进行内插，以便获得在所述深度通过井的所述截面的流出物的至少一种相的体积分数分布线(profile)和速度分布线。

通过在所述深度在井的所述截面上对所述的至少一种相的体积分数分布线和速度分布线的积分计算该所述的至少一种相的相对体积流率。

这种计算流出物的相的体积流率的方法允许考虑能够改进或改变局部测量结果的任何与测量装置或流动状况有关的因素。

附图说明

在参考了详细描述和提供的图后，本发明的另外的目的和优势将会对本领域中的普通技术人员变得明显。

图1A和图1B代表了可以对速度测量采取的第一校正步骤的示意性例子；

图2A和图2B代表了可以对该组局部速度和持有率测量值采取的第二可能的校正步骤的示意例子；

图3A和图3B是第三可能的校正步骤的示意例子；

图4代表了一个根据选定的流体模型在局部持有率测量结果上实现的内插法的例子；

图5是一个采用本发明的方法的一个例子的流程示意图；

图6A、6B、6C给出了根据本发明测量和确定的持有率值的曲线图；

图7A、7B、7C给出了根据本发明测量和确定的速度值的曲线图；

图8给出了根据泊肃叶模型(Poiseuille model)的管道中流速的表达；和

图9给出了根据EP 0519809中的模型的管道中流速的表达。

具体实施方式

正如在WO 01/11190中申请人已经解释的那样，为了确定在井的给定深度处的多相流出物的流体的体积分数和所述流体的相对体积流率，我们需要确定每种流体(通常为油、水和气体)的局部体积分数(通常叫持有率(hold-up))和这些流出物中每种流体的局部速度。正如这篇文献里所描述的那样，给定流体(i)(在这里i可以是w，即水，o，即油，g，即气体)在某个时间(t)的体积流率Q_i(z，t)是通过对当前深度(z)的一个局部流体持有率(h_i)与井截面(s)处的局部流体速度(v_i)的乘积的简单积分计算的，计算过程遵循下述原则：

v_i和h_i是三维位置和时间(x，y，z，t)的两个函数。

令z是与当前深度的井截面垂直的方向，x是与垂直平面和管道截面相交的方向。那么，在时刻t，而且i＝w，o，g，局部体积流率为：

$Q_i(z,t)=\underset{s}{\∫}{v}_{i.}{h}_{i.}\mathrm{ds}$

如果井是垂直的，体系的动力学性质使得速度和持有率仍然是流体的函数和深度(z)的函数，但是在井的横截面S上是一致的(处于第一次近似)。那么上述关系式简化成了：

               Q_i(z，t)＝Vi(z).Hi(z).S(z)

其中，V_i和H_i在横截面S的范围内分别是v_i和h_i的恒定值。

如果井不是垂直的，由于在该横截面的范围内v_i和h_i不再一致，就不能进行前面的简化。因此，为了计算每个深度处的Q_i，在该井横截面的范围内就需要对速度和持有率分布线进行更为精确的描述。在根据本发明方法中的第一步骤中，通过设置在该横截面的不同已知位置上的局部传感器采集每种相的局部持有率和局部速度的测量值。

因此，在多相流体在斜井中流动的情况下，Q_i的计算就成为一个更为复杂的过程，这是因为它涉及基于流动模型的测量值的内插法，而流动模型可以随着当前的流动状况而改变的：在这种类型的内插法中，井斜(deviation)是一个主要的影响参数，而且，在那些井中能遇到的流动模式的范围是非常宽的。因此，考虑到任何能够改变流动模式的因素而且仔细检查收集的任何局部测量值是非常重要的，这是因为很小的缺陷或错误在确定流动模型并且因此确定局部体积流率时也能够导致大错误。

因此，业已研究出本发明的方法，以克服影响从井中获得的局部体积分数和速度测量值的极大数量的错误。因此，基于非常接近真实流动状况的测量值，可以确定最适合的流动模型，该流动模型使得可以在井横截面的范围内获得可靠的相体积分数和速度分布线。为了确定相对体积流率，本发明的方法的主要步骤包括一个校正步骤、一个选择流动模型的步骤、一个内插步骤和一个计算步骤。然而，注意到这些步骤几乎可以按照几乎任何顺序实现是很重要的。例如，首先选择流动模型然后校正测量值是可能的。相同的步骤进行两次也是可能的。例如，选择一个流动模型、进行校正步骤，然后根据已经进行的校正，选择第二个流动模型，这个过程是可能的。对局部测量值进行部分校正、选择流动模型然后再对这些测量值进行校正也是可能的。最后，在这里详细介绍的校正在下面只有部分被实施，而且校正步骤可以包括变动数量的校正。

本发明的方法需要获得两个主要的局部测量值：局部持有率(不同相的体积分数)和局部速度，这正如前面已经解释的。然而，为了计算每种相的相对体积流率还需要其它数据：缆索速度(cable velocity)和井的几何形状。局部速度是通过向下放到井中的工具获得的(注意到那些数据或者在把工具向下放到井的过程或者在把工具撤回的过程中获得)，从所述的局部速度测量值中减掉缆索速度是很必要的。井的几何形状包括了解其偏离程度和直径(当被假定为圆形时)或者成椭圆形的程度。为了计算传感器的位置，了解这工具在井中的位置(例如相对方位(relative bearing)：工具居中间位置、非居中间位置或倾斜)也很重要。了解流体的流变性(粘度、密度)也将有助于优化数据的解释。

局部持有率的测量使得可以确定在井中特定深度出现的流体的局部体积分数。根据流出物的组成，测量方法可以或者允许三相的差异(油、水和气)或者两相的差异(例如用电探针测油和水)。例如，测量装置(measurement means)可以是测量水持有率(water hold-up)的一些局部电探针和/或测量气体持有率(gas hold-up)的一些局部光探针。如果流体是三相的，就需要把这两种探针组合起来使用：油持有率(oil hold-up)将很容易地从其它两种相的持有率推导出。这种为了推导出油持有率而把所有探针的测量组合起来的方法可以在本发明的方法中的任何一步进行。如果流出物是两相的，仅仅一种类型的探针是必需的，而且不需要任何组合。如果流出物是单相的，持有率问题是不涉及的。测量手段也可以包括一些将获得局部速度值的纺锤形流量计，参见WO 01/11190。

下面描述的是能够用于局部测量的不同的校正方法，这些校正方法中任何一种是不是包含在本发明方法的校正步骤中，依赖于测量装置和流动状况。

在第一个实施例中，本发明的校正步骤可以包括对已知的系统测量误差的校正。事实上，测量装置的精确度经常是受到系统误差的影响。通常，它们不能被考虑在传感器校准中，典型的原因是它们依赖于一些其它的测量参数，而这些参数只有在进行测量时才能得到。在本发明的一个优选实施例中，这个校正是校正步骤的第一步。例如，对系统误差的校正将涉及纺锤形流量计。纺锤形流量计的效率可能受这个传感器在工具上的位置的影响。根据这个位置，速度的测量能够受工具自身引起的流动扰动的影响。那么，将要进行的校正将是一个简单的校正模型中的固定系数。然而，它也可以是一个扩展的校正模型中的函数，例如是井的截面积、工具位置(相对方位)和实际测量速度的函数。其它的系统误差也能够在局部持有率传感器或者光学探针上产生。

在第二个实施例中，校正步骤也可以包括对相的体积分数和/或速度惯用的(acquired for)局部测量值的检查。实际上，业已通过测量装置获得的一组局部测量值(相体积分数或速度)可以显示出一些足以在进行任何其它处理之前就能被检测到的不相关性(incoherence)。因此，这个校正的目的是检查同样种类的所有测量数据之间是相关的(coherent)。如果确定了不相关性，将进行一次校正。图1a和图1b给出了对局部速度测量值进行这种校正的一个例子。图1a显示了在井截面上的位置0、1、2、3上垂直地重新分配速度传感器，在该井截面上流动是稳定的，图1b显示了对一组错误的测量装置的校正。

如果传感器1和传感器3测得一个正的(不是零)速度值，那么传感器2(设置在传感器1和传感器3之间)应该也测得一个正(不是零)的值。如果传感器2测得一个负的值或是零，那么它的测量就可以被认为是错误的。一个简单的校正是决定忽略掉这次测量而同时进行下面的解释。另一种解决办法是如图1b所示那样把传感器2的值设定成传感器1和传感器3的测量值的平均值。

现在假设图1a代表了一个水/油两相稳定流，4个电探针沿着井截面的垂直直径(the vertical diameter)设置，所述的探针从底部到顶部设置，以数字0到3表示。如果传感器0显示0％的水持有率而探针1或2显示严格地正的水持有率值，那么至少一个传感器的测量将不得不被认定是错误的，因为这组测量值不符合基本的物理原理。一个简单的校正例如就是或者忽略这个测量值或者给传感器0与传感器1的测量值相同的值。

校正步骤也可以包括检验数据间的相互依从性(mutual compliance)。一旦最简单的测量错误已经得到了修正，就可以进行一个更高级的检验。在比较了同种类型的测量的相关性后，就需要检验它们之间的相互关联性，当已经确定时，对它们进行纠正。典型的是，局部持有率测量值将与局部速度测量值进行比较。

图2a和图2b显示了这种校正的一个例子。在油/水两相层流动(假定是稳定的)的情况下，水的速度和油的速度大致相同。让我们假设持有率传感器和速度传感器是耦合的：数对由一个持有率传感器和一个速度传感器组成的测量元件共同沿着井截面的垂直直径设置，如图2a所示的，速度传感器以数字0到3表示，在流动部分的左手侧；水持有率传感器以数字4到7表示，设置在它的右手侧。注意在下面的部分《o》代表油，《w》代表水，《g》代表气体。

如果偏离度不是90°(水平井)，V_o和V_w是不同的。让我们以偏离度为89°的情况为例。那么根据已知的物理原理，V_o大于V_w。

如果速度测值是图2b右手侧的图所示的那样，那么速度测量值和在左手侧图上所显示的持有率测量值之间存在明显的不相关性。的确，持有率2表明油/水界面低于它当前的位置而速度传感器6趋向于表明这个界面高于所述的传感器。那么，根据对所测速度和所测持有率的相对信赖程度，就需要对两个数据中的一个(可能两个)进行修正。

也可能在校正步骤里加入旨在改进这组测量值的进一步的校正。这个校正是以对速度测量值的更进一步分析为基础的，以便把这个分析的结果通过一个将使这个持有率体系精确的虚拟点注入到持有率测量值中去。也可能运用相反的程序并且分析这组持有率测量值以便注入一个虚拟的速度测量值，这会使速度体系的分析结果更为准确。

一个这种改进校正的例子可以在图3a和图3b中看到。流动形式假定为油/水两相层流(稳定的)，并且水的速度和油的速度大致一致。在图3a中，持有率探针和速度传感器不是沿着油井截面的垂直直径共同设置的，而是交叉定位的，持有率传感器以数字0到3表示，在井的左手侧；速度传感器以数字4到7表示，在右手侧。如图3b所表示的，该组持有率测量是用来改善这组速度测量值的。持有率传感器2测得水体积分数为60％。这表明这个传感器是设置在油/水界面上。然而，速度传感器不能察觉到这个界面，因为它们或者设置在油里或者设置在水里。因此，从给出界面的位置的持有率测量值，可以把这个信息融合到速度测量值中。这可以通过增加恰好在界面上面和下面分别靠近上方的真实的测量点6和下方的真实测量点5的真实的速度点而实现。

如果选择的模型是在一个接近水平的井(例如90°的偏离度)中的液体/气体两相流动，持有率传感器探测不到任何气体也是可能的，因为气体是在井截面的最上部。如果最高的速度传感器在最高的持有率传感器上面，这个速度传感器可以探测到增加的速度。这个特性将表明气体的存在，而未经校正的持有率内插分布线将忽略它的存在。因此，能够把一个真实的气体持有率传感器固定在与最高的速度传感器同样的位置，它的测量值可以被设定为它探测到100％的气体(在90°的偏离度和当前流动速度的情况，出现层流是合理的)。

注意到校正步骤可以包括这些校正过程的一个、几个或全部是很重要的。此外，这些校正可以无关紧要地轮流地、按任何顺序完成，或者也可以把允许内插这些组测量值的流动模型的选择包括在它们的一个之间。

流动模型是一个相应于已知流动模式的预先确定的持有率分布线和速度分布线的组合。为了描述一个在井中可以碰到的任何类型的流动，创造出不同的流动模型。已知的流动模型例子如下：

图8所示的泊肃叶模型是代表在柱状垂直管道中层状单相流动的传统的模型。这个模型适合半径对速度测量值的一维两项式函数。不涉及持有率测量。然而，在很多涉及井流出物的情况下并不使用这个模型，这是因为速度通常对于维持层状的流动来说太高，在多数情况下流动是湍流状态。

另一个模型在申请人的专利EP0519809中被披露，它适合分离的油/水流动。在这个模型中，油和水的速度假设为相同的。这个模型是基于油的速度和水的速度是常数这样一个假设。它如图9所示。

V_o和V_w分别代表横跨柱状管道截面S的油的速度和水的速度。H_w和H_o分别代表在测量截面中水的体积分数和油的体积分数。S_w和S_o分别代表在截面中的被水覆盖的管道表面和在所述截面处被油覆盖的管道表面。

流动模型的选择可以由负责实施本发明方法的人进行，或者由一个自动决策链(automatic decision chain)实行。在“人工”决策的情况下，是由人来监控本发明的方法，他根据自己的关于流出物流动、井的偏离等知识来选择给定的流动模型。在流动模型的选择是自动的情况下，这个过程将检查几个预先设定的标准以便决定选择使用哪个流动模型。例如，它将检查该流动是否是油/水流动：如果电探针探测到油和水，而光学探针没有探测到气体，该流动是一个油/水两相流动。而后，这个过程将检查这个流动是不是分层：如果井偏离接近水平而测得速度均低于最大值(例如2m/s，在15.24厘米(6英寸)直径的管道中)，该流动可能是分层的。另一个使用层流模型的标准可以是井偏离接近水平而除一个探针以外的探针测得水的持有率为1或者0。在自动过程识别一个油/水层流的情况下，例如，它将使用专利EP0519809所披露的流动模型。

根据本发明的方法，定义流动模型的步骤可以在任何测量值校正之前进行，也可以在两个校正之间或者整个校正过程结束之后进行。但是，如果流动模型在任何校正之前已经选定或者在选定所述的流体模型之后已经对测量值进行其它校正，检查流动模型是不是仍然合适或许会有用。此外，如果这个模型改变了，检查所有已经进行的校正是不是仍然关联可能有用。另外，校正步骤和流动模型的选择是按照一种反复的方式实现的，直到体积分数测量值和速度测量值适合最后选择的流动模型。

在校正步骤和流动模型的选择实现后，将对局部速度和体积分数测量值进行内插处理。事实上，这些测量值没必要完美地符合所选的流动模型所需要的输入。因此，需要一个转换步骤把这些局部测量值转化成“可输入”的格式。这可以通过一个简单的单位转换实现，但是也可能是数据物理性质的改变。例如，在一个斜井中的液体/气体两相流动中，局部速度测量值代表了混合物的速度。因此，它既不是液体的速度，也不是气体的速度。内插步骤可以应用到局部混合速度，然后液体速度分布线和气体速度分布线就可以从局部混合速度中推导出来。但是，也可以在一侧对局部液体速度进行内插，而在另一侧对局部气体速度进行内插。在这种情况下，需要在利用液体和气体这两种测量值完成内插步骤之前推导出液体和气体局部速度。事实上，为了从测得的混合速度计算这些局部速度必须使用第二种滑动模型(slippage model)。任何适合这种类型流出物的已知模型都可以采用。那么，根据选择的流动模型，对速度测量值和持有率测量值都完成内插步骤。图4所示的例子表示了一个内插步骤的结果，这个内插步骤采用了一个基于最小平方法的流动模型，是在从一个直径为20.32厘米(8英寸)的管道获得的四个局部水持有率测量值完成的。

内插步骤的结果将给出一组分布线。对于一个三相流出物，这些分布线将是一个水体积分数(持有率)分布线、一个油体积分数分布线曲线、一个油速度分布线和一个气体速度分布线。

本发明的方法的最后一个步骤是计算相的相对体积流率，这通过在进行局部测量的那个深度处的所述井截面上对体积分数和/或速度分布线进行积分得到的。

这种内插法是根据下面的公式完成的：

$Q_i\left(z\right)=\underset{S}{\∫}{v}_i.h_i\mathrm{dS}$

其中z是沿着井轴的坐标(测量的深度)，S是当前截面的面积，v_i是相i(i＝w，o或g)的速度分布线，h_i是相i(i＝w，o或g)的体积分数分布线。因此，可以看到，在这个阶段为了计算在任何深度井截面的面积，与井的几何参数有关的信息很重要。这种类型的信息将由已知的井径仪(callipertool)给出。

如果两个h_i和v_i分布线都是仅仅一个变量x(也就是当流出物在一个倾斜的但是圆形的井中流动)的函数，那么计算式变成了：

$Q_i\left(z\right)={\∫}_{-D/2}^{D/2}{v}_i{.h}_i.y\left(x\right)\mathrm{dx}$

其中D是井的直径，v_i是相i(i＝w，o或g)的速度分布线，h_i是相i(i＝w，o或g)的持有率分布线，y(x)是与坐标轴x垂直的井在位置x处的宽度。

如果两个分布线h_i和v_I都是相同的，那么计算式变成了：

Q_i(z)＝V_i(z).H_i(z).S(z)

其中，V_i(z)是v_i在深度z处的平均值，H_i(z)是h_i在深度z处的平均值，S(z)是井在深度z处的截面的面积。

因此，本发明的方法可以根据从井中获得的局部体积分数和速度测量值来计算在给定深度下的流出物每种相的相对体积流率。在井的获得局部测量的任何深度处重新应用这种方法将给出一个流出物流动特性的意见。那么，将可能通过把重点放在在特定区域并降低一些其它区域的产量来管理所述的生产井。

图5是实施本发明的方法的一个例子的流程示意图。方框代表了本发明的方法的不同步骤。在本发明的方法中的这个例子中可以看出，系统误差的校正是在选择流动模型之前进行的。此外，如前所述，在选择流动模型后的校正可以导致所述模型的新选择。

现在，本发明的方法在包括几个数据的一组数据中进行了测试。在以下，一个三相流出物的流动状况是已知的，并且使得：管道内径是6英寸，管道的偏离度是88°，水体积流率是10m³/h，油体积流流率是10m³/h，气体体积流率是6.9m³/h。在进行任何测量前，流动状况是不变的并且是稳定的。然而，流动状态在这种状况下是不稳定的：长的气泡在部分分离的油/水流顶部上面流动的这种状况下。

测量装置有五个固定在已知位置的沿着管道的垂直直径分布开的速度传感器(微型纺锤形流量计)。按照同样的方式，六个光学探针和六个电探针分别测量通过管道的水持有率和气体持有率。

这个例子中使用的数据是用上部的传感器测得的数据，在接近1分钟的周期内取平均值，再加上一些其它参数(流体流变形、传感器位置......)。

流动模型的选择在本申请中是用户的一种选择。图6a到6c和7a到7c描述了得到相对体积流率的连续步骤。

图6a给出了通过六个电探针测得的持有率值◇和通过光学探针测得的持有率值△。图6b给出了通过本发明的方法获得的水持有率(曲线A)和气体持有率(曲线B)的内插分布线。为了使用本发明的方法而已经选择的流动模型是基于最小平方法。图6c另外示出通过简单地减去两种其它的分布线得到的油持有率分布线(曲线C)。

图7a描绘了从五个纺锤形流量计的旋转速度推导出的混合物的局部测量速度◇。图7b描绘了与气体速度分布线(profile)相应的曲线(curve)，而该气体速度分布线是通过使用经过内插的气体持有率分布线和利用根据本发明的方法的滑动模型得到的气体/液体滑动速度而推导出的。图7c另外分别描绘了油和水的速度分布线(曲线C和A)，这是利用本发明的方法中的液体速度分布线(通过油和水的持有率分布线计算得到)和一个滑动模型以及油水滑动速度确定的。该方法的最后一步是以如下分布线为基础给出的相的体积流率：Q_w＝9.64m³/h，Q_o＝9.92m³/h，Q_g＝6.52m³/h。这些值表明，根据一个三相流动测得的数据本发明的方法能够给出令人满意的结果。

因此，本发明的模型使得可以得出在井中流动的流出物的每种相的体积分数和速度分布线，所述的分布线最大可能地接近真实的流动状况：实际上，在保持所有的这些测量值相互关联的同时，大多数扭曲局部测量值的因素都被仔细地进行了确定和校正。

Claims (12)

1.一种计算在一井中流动的一多相流出物的至少一种相的相对体积流率的方法，所述方法包括获取所述相横跨处在一深度处的所述井的一截面的局部体积分数和/或速度的一第一步骤，其特征在于，所述方法还包括：

校正所述局部体积分数和/或速度测量值，以便使它们相互一致和/或与所述流出物的流动状况相一致；

选择以数学的方式代表所述流出物流动的一合适的流动模型；

通过所述选择的流动模型对所述局部体积流率的测量值和/或所述局部速度测量值进行内插，以便获得横跨处在所述深度处的所述井的截面的所述流出物的至少一种相的一体积分数分布线和/或一速度分布线；

通过在所述深度处的所述井的截面上对所述体积分数和/或速度分布线积分，计算所述至少一种相的所述相对体积流率。

2.根据权利要求1所述的方法，其中所述校正步骤包括对由测量装置引起的系统测量误差的校正。

3.根据权利要求1或2所述的方法，其中所述校正步骤包括检查所述局部体积分数测量值它们自身之间的相关性和/或局部速度测量值它们自身之间的相关性。

4.根据权利要求1到3之中任意一项所述的方法，其中所述校正步骤包括检查定所述局部体积分数测量值与所述局部速度测量值之间的相关性。

5.根据以上权利要求之中任意一项所述的方法，其中校正步骤包括根据所述局部速度测量值和所述局部体积分数测量值各自的分析结果，对所述局部体积分数测量值和所述局部速度测量值进行改进。

6.根据以上权利要求之中任意一项所述的方法，其中所述合适的流动模型的选择是通过分析成组的体积分数和/或速度测量值而自动进行的。

7.根据以上权利要求之中任意一项所述的方法，其中还包括下面的步骤：

测量把所述体积分数和速度测量装置向下放到所述井中所使用的缆索的速度；

确定所述井的几何特征。

8.根据以上权利要求之中任意一项所述的方法，其中还包括所述井的相对方位的确定。

9.根据以上权利要求之中任意一项所述的方法，其中所述流动模型的选择是在对所述局部体积分数和速度测量值进行了一些校正后被证实的。

10.根据以上权利要求之中任意一项所述的方法，其中按如下所述完成以下步骤：

校正由于测量装置引起的系统测量误差；

选择一合适的流动模型；

检查所述局部体积分数测量值它们自身之间的相关性和/或所述局部速度测量值它们自身之间的相关性；

检查所述局部体积分数测量值与所述局部速度测量值之间的相关性；

根据所述局部速度测量值和所述局部体积分数测量值各自的分析结果，对所述局部体积分数测量值和所述局部速度测量值进行改进；

重复前面的步骤，直到所述局部体积分数和/或速度测量值适合业已选择的所述流动模型；

通过所述选择的流动模型，对所述局部体积分数测量值和/或所述局部速度测量值结果进行内插，以获得横跨处在所述深度处的所述井的截面的所述流出物的至少一种相的体积分数分布线和/或速度分布线；

通过在所述深度处的所述井的截面上对所述体积分数和/或速度分布线积分，计算所述至少一种相的相对体积流率。

11.根据以上权利要求之中任意一项所述的方法，其中所述局部体积分数测量值是所述井的所述流出物中的水体积分数和/或气体体积分数，所述流出物中的油体积分数从前者推导得出。

12.根据以上权利要求之中任意一项所述的方法，其中所述局部速度测量值是通过一组纺锤形流量计获得的。

Images