专利名称：一种精确测量气液两相混合流体中气相流量和液相流量的方法
技术领域：
本发明涉及流体流量测量领域，更具体地，涉及一种测量气液两相混合流体中的气相流量和液相流量的方法。
背景技术：
采油工业中，经常从油井中采出包含液相和气相的气液混合流体，业内常称之为“湿气”。其中所述气相包括例如空气、油田气或任何在常温下不凝的气体，其中该油田气一般为较轻的烷烃甲烷、乙烷、丙烷、丁烷等。所述液相可包括原油本身以及采用过程中注入油井中的水以及其它液体添加剂。实时准确地测量从油井中采出的气液混合流体中气体的体积流量和液体的体积流量，是生产管理和生产优化所必需的基础数据。当前，已经有一些方法来实现这种气液两相混合流体中气体体积流量和液体体积流量的在线测量。传统方法是通过分离器将气液两相混合流体分离成气相和水相，然后分别计量它们的体积流量。但由于分离器以及相关附属设施重达数十吨，占地上百平方米，且控制环节多，使得分离器的维护和管理十分复杂，不利于生产过程管理的自动化，尤其不利于在沙漠和海上油田中使用。另一类方法是不对气液两相进行分离，通过测量气液两相混合流体的总体积流量Q和气相分率GVF来测量气相体积流量和液相体积流量。其中气相分率GVF是指气相流量占气液两相混合流体总流量的百分比，亦可表示为某一横截面处，气相所占据的面积与整个截面面积之比。理想情况下的处理方法是气相体积流Qg = QXGVF,液相体积流量Q1 =QX (I-GVF)0但这种测量方法要基于一个前提气液两相在该截面上具有相同的速度。但是，实际上，由于气液之间在密度、粘度等性质方面差别巨大，在实际管路中，上述前提并不成立。在实际管路中，由于气体密度较小，气相速度常常大于液相速度，这使得上述理想情况下测量的气相体积流量和液相体积流量出现了一定误差，故需要一种改进的方法来对上述理想情况进行校正。本领域中将气液两相混合流体在同一管路中的气相速度与液相速度之差称为滑差，而将气液两相混合流体在同一管路中的气相速度与液相速度之比称为气液滑差因子，下文中将对气液滑差以及气液滑差因子进行详细介绍。目前气液两相混合流体的计量技术中对于由于气液滑差引起的计量误差的校正主要有两种方法—种方法是基于气液两相的Lockhart-Martinelli参数进行迭代计算,主要用于极高含气的湿气计量中的气量虚高修正，主要代表有ISO湿气模型，例如参见出版物ISO/TRl 1583:2012,其英文题目是 “Measurement of wet gas flow by means of pressuredifferential devices inserted incircular cross-section conduits，，。这种方法都是使用单相仪表将湿气当做单相进行计量，不会使用多相仪表，所以液量需要通过别的途径预先知道，或知道个概值，所以修正都是针对气量的修正。这种方法有三个主要缺点一是没有明确的动力学机制；二是只修正了气量，没有对液量计量进行修正；第三，该方法的适用范围仅限于极高含气的很窄含气率量程。另一种方法基于经验模型，比如实验数据拟合模型、或者信号处理模型的两相流流量气液滑差修正。这种方法有两个缺点经验模型对实验数据和测量条件的依赖性很强，无法做到通用性和精度兼顾。另外因为多相流的流态复杂性，使得经验模型通过对建模参量选择和数学关系作出更合理的假定非常困难。信号处理完全把多相流当作黑箱处理，计量精度和范围都有很大的人为特征和不确定性。因此，本领域需要一种能够简单地且精确地测量气液两相混合流体中气相体积流量和液相体积流量的方法，且还要能够对气液滑差引起的测量误差进行校正。发明概述本发明的第一方面涉及一种精确测量气液两相混合流体中气相体积流量和液相体积流量的方法，包括以下步骤·
a.使气液两相混合流体流过竖直设置的圆管，以形成环状气液两相流；所述环状气液两相流是指这样的一种流型气体以圆柱状围绕圆管轴线，而液体以环状分布在气体与管壁之间；b.测量该气液两相混合流体的总体积流量Q、温度T、压力P和气相分率GVF ；c.基于以下公式计算气液滑差因子S和气体体积流量Qg和液体体积流量Q1 S = 2-GVF/(u kX (GVF— I))Qg = QXGVFXS/(1-GVF+S*GVF)Q1 = QX (I-GVF)/(1_GVF+S*GVF)其中，Q为总体积流量，单位为m3/s ；S为以m/s计的气相速度与以m/s计的液相速度之比，无量纲；GVF为气相分率，无量纲；u k为在气液两相混合流的测量温度T和测量压力P下以任何粘度单位例如Pa-S计的气体粘度U g与以任何粘度单位例如Pa s计的液体粘度U !之比，无量纲。本发明的第二方面涉及一种精确测量气液两相混合流体中气相体积流量和液相体积流量的方法，包括以下步骤a.使气液两相混合流体流过竖直设置的圆管，以形成环状气液两相流；所述环状气液两相流是指这样的一种流型气体以圆柱状围绕圆管轴线，而液体以环状分布在气体与管壁之间；b.测量该气液两相混合流体中的液相平均速度VI、温度T、压力P和气相分率GVF ；c.基于以下公式计算气液滑差因子S和气体体积流量Qg和液体体积流量Q1 S = 2-GVF/(u kX (GVF-1))Qg = VlXAXGVFXSQ1 = VlXAX (I-GVF)其中，Vl为液相平均速度，单位为m/s ；
S为以m/s计的气相速度与以m/s计的液相速度之比，无量纲；GVF为气相分率，无量纲；U k为在气液两相混合流的测量温度T和测量压力P下以任何粘度单位例如Pa-S计的气体粘度U g与以任何粘度单位例如Pa s计的液体粘度U I之比，无量纲；A为圆管的横截面积，单位为m2。本发明的第三方面涉及一种精确测量气液两相混合流体中气相体积流量和液相体积流量的方法，包括以下步骤a.使气液两相混合流体流过竖直设置的圆管，以形成环状气液两相流；所述环状气液两相流是指这样的一种流型气体以圆柱状围绕圆管轴线，而液体以环状分布在气体与管壁之间；b.测量该气液两相混合流体中的气相平均速度Vg、温度T、压力P和气相分率·GVF ；c.基于以下公式计算气液滑差因子S和气体体积流量Qg和液体体积流量Q1 S = 2-GVF/(u kX (GVF— I))Qg = VgXAXGVFQ1 = VgXAX (I-GVF)/S其中，Vg为气相平均速度，单位为m/sS为以m/s计的气相速度与以m/s计的液相速度之比，无量纲；GVF为气相分率，无量纲；u k为在气液两相混合流的测量温度T和测量压力P下以任何粘度单位例如Pa-S计的气体粘度U g与以任何粘度单位例如Pa s计的液体粘度U I之比，无量纲；A为圆管的横截面积，单位为m2。附图
简述图I是气液两相混合流体在竖直管路内形成环状流后的截面示意图。其中，中心白色部分为气相，环状阴影部分为液相。图2是气液两相混合流体在竖直管路内形成环状流后，在某一截面处，气体速度和液体速度沿径向分布的示意图。图3是气液两相混合流体在竖直管路内形成环雾状流后的截面示意图。其中，中心白色部分为雾，由气相和散布其中的液滴组成，环状阴影部分为液相。图4为本发明的实施例中使用的测量方法的流程示意图。图5是采用本发明环状流模型修正后的气相流量与真值相比的相对误差以及未进行修正的气相流量与真值相比的相对误差图。图6是采用本发明环状流模型修正后的液相流量与真值相比的绝对误差以及未进行修正的液相流量与真值相比的绝对误差图。以上附图仅仅是为了举例说明而提供，不打算以任何意图来限制本发明。本发明的保护范围仅由权利要求确定。发明详述本发明中使用的术语的定义如下
“气液两相混合流体”是指运动中的气液混合物，其中气相与液相可以以任何方式混合。例如，气液两相可以具有连续的边界，形成所谓的“分层流”;或者液相悬浮于气体中，形成所谓的“液滴流”;或者气相悬浮于液相中，形成所谓的“气泡流”;或者管路内某一段是气体，另某一段为液体，形成所谓的“栓塞流”，等等。“流量”是指单位时间内流过某截面的流体的体积。流量有多种计算方式，按时间分，有瞬间流量和积分流量之分，瞬间流量是指在某一时刻即时间点时的流量，积分流量是指在某一时间段内的流量，可以认为是瞬间流量在该时间段内的积分值，故此得名。生产上关心的是积分流量。故本发明中，除非另有专门说明，所提及的“流量”均是指积分流量。“速度”，也称作“流速”，是指流体单位时间内流过的路程。“速度”这一术语在物理学和流体力学中本来是矢量，但在本发明中，将其作为标量来处理，认为管路中所有的流体的速度都是平行于管路轴线方向并指向流体前进方向的，故仅考虑其大小即可，而忽略其方向。这也是符合采油工业的实际需要的，因为人们通常只关心速度的大小，不关心其方向如何。“速度”又可细分为在截面上某一点处的速度和在整个截面上的平均速度，前者称为点速度，后者一般简称为截面平均速度。从概念上讲，截面平均速度是点速度在截面上的面积分值。本发明在使用速度的概念时，除非专门说明是点速度，否则都是指截面平均速度。·在多相流中，某一相的速度，通常应该是指该相在其自己截面上的平均速度，本专利中的速度比也是各相的“平均速度”之比。“气液滑差”是指气液两相混合流体在同一管路中的气相速度与液相速度之差。“气液滑差因子”是指气液两相混合流体在同一管路中的气相速度与液相速度之比，无量纲。“相分率”，是指在多相流的某一流通横截面(在管道中流动时，即为管道横截面)上，某一相所占据的面积与该流通横截面面积的比例。根据所指称的相的不同，该术语又可分为“气相分率”(简称为GVF，又称截面含气率)和“液相分率”(又称截面含液率)，分别指气相和液相在某一流通横截面上的所占据的面积与该流通横截面面积的比例。相分率可以用相分率计来测量。“稳态流”是指流体的流型在宏观上不随时间变化，即达到所谓的“稳态”。在本发明的第一方面的步骤a中，使气液两相混合流体流过竖直设置的圆管，以形成环状气液两相流；所述环状气液两相流是指这样的一种流型气体以圆柱状围绕圆管轴线，而液体以环状分布在气体与管壁之间。一般认为，当气液两相混合流体在水平或倾斜管道中流动时，由于液相的密度要大于气相的密度，故液相会更多地分布在管道的中下部位置，这给气液两相各自的体积流量的计量以及校正气液滑差代来的误差造成了困难。为避免这种情况，本发明使使气液两相混合流体流过竖直设置的圆管，以形成气体在管中心流动而且液体在管壁和气体之间流动的流型，本发明中将其称作环状气液两相流。这样的环状气液两相流在测量上的优势是在稳态下气液速度场是围绕圆形管道的轴线对称分布的，这给气液滑差的校正带来了便利。本发明人通过大量的实验发现，通过控制以下条件表观气速> 15m/s，可以在竖直圆管内稳定且可重现地实现这种环状气液两相流，其中表观气速的测量可通过其它传统的测量方式加以粗略地测定，例如通过将气相和液相进行分离后测量气相的体积流量，再除以管路横截面积得到表观气速。在本发明的第一方面的步骤b中，测量该气液两相混合流体的总体积流量Q、温度T、压力P和气相分率GVF。其中，通过使用相分率计来在线测量所述气相分率GVF，所述相分率计选自单能伽马射线相分率计或双能伽马射线相分率计。这些相分率计都是可商购得到的成熟产品，本领域技术人员知晓去何处购买以及如何使用，故在此不再赘述。通常认为，所测得的气相分率GVF就是某一截面处气体所占据的面积与截面总面积之比。在本发明的环状流的情况下，气相分率GVF = r(l2/R2，其中r°是如图I所示的圆柱状气体的半径，R是指圆形管的内径。温度T可以通过任何常规测温手段获得，例如通过热电偶或温度计来测量。压力P可以使用任何常规测压设备例如压力传感器或压力表来测得。总体积流量Q可以通过使用任何常规体积流量测量方法来测量，例如使用文丘里流量计、孔板流量计、转子流量计来测量，这些设备和方法都是已知的，且所使用的设备均可商购得到。具体来说，在本方面的一种实施方法中，使用文丘里流量计来计量总体积流量。该文丘里流量计可安装在管路任何位置，但优选安装在该竖直圆管中的某一段上，最优选安装在该竖直圆管中部，并距离该竖直圆管的入口和出口都有一定距离，例如大于该圆管直径的6倍以上的距离，以使得测得的流量比较准确。一般地，文丘里流量计都配有差压变送器，用于测量该流量计入口处与其喉部的压差，该压差可用于计算总体积流量Q，具体计算公式可从任何一本流体力学教科书或从流量计的使用说明书中查找。使用直接测量体积流量的流量计也是可以的。·
在本发明的第一方面的步骤c中，基于以下公式计算气液滑差S和气体体积流量Qg和液体体积流量Q1 S = 2-GVF/(u kX (GVF - I))Qg = QXGVFX S/(I _GVF+S*GVF)Q1 = QX (I-GVF)/(I_GVF+S*GVF)其中，Q为总体积流量，单位为m3/s ；已经由步骤b测得。S为以m/s计的气相速度与以m/s计的液相速度之比，无量纲；GVF为气相分率，无量纲；已经由步骤b测得。u k为在气液两相混合流的温度T下气体粘度U g与液体粘度U I之比，无量纲。认为气体粘度Ug与液体粘度U1在本发明中是容易获得的物理量，且可以以任何粘度单位例如Pa *s计来表示。例如，测定温度T后，通过常规的气液组成通过查表或数值计算获得气体粘度y g和液体粘度y:，或者，通过任何常规方法将液体和气体分开后，分别用相应的仪器例如粘度计来进行测定。值得一提的是，本发明中对气液滑差因子S的计算采用了现有技术中从未采用的方法。其具体计算过程如下本发明充分利用了步骤a中在圆形管道中得到的气液两相环状流的对称性和稳态性，并结合牛顿流体等特征(采油工业中所遇到的流体均为牛顿流体，关于牛顿流体的概念，可参见任何一本流体力学教科书)，简化描气液述两相流体的普适公式Navier-Stokes动量方程，求出管道截面中的气液两相的速度分布的解析解，然后对气液各自的截面域积分后求出气液两相的流量，再利用气液两相体积流量和气液两相面积计算出气液两相各自在某一截面上的真实平均速度，二者的比值即为气液滑差因子的精确解析解表达式S。具体推导过程如下I.求解环状流中气液两相在柱坐标系下轴向速度沿径向分布的速度流场解析解形式，如图I所示首先，可以列出沿直的圆形管流动的环状流的Navier-Stokes方程在柱坐标系下的表示形式如下(这里只列出动量守恒方程，对所有流体状态及物性量通过下标g或I以区分气相或液相，其余各符号的含义见说明书末页的符号含义说明)
「 n (dur dur u0 dur dti 0 V .p —^ + ^ + ―—^ +~^ (I)
v d ' or rcz r J
dp |~1 5 f dur) I d2u d2u u 2 diia·—--￡_+"--尸-L. -J---L-f--1---1-----f 02
Br [r5rV dr J r2 d02 d=2 r2 r2 o0 J rpi^ + ^ +++
、dt dr r 50 dz r J
dp Idf durA) I d2n0 d2u0 u& 2 dur
=—i + # 一7T r+~~^0~ + ~^2---~ ++ PS-0
a0 r dr Sk dr J r 50' dz r r 80
「 I foil du_ H0 du7 dtt'),、p + ++(3)
、ot or r o0 dz J
dp Ildf 一) I S2U7 d2tL
—---h U--T-— H—r--— H--+ 02.
dz _ r dr v Sr J r' 30' d: _根据气液两相的动量方程、对称性、牛顿流体假设、稳态、绝热、相间界面的连续性、管壁处的无滑移，并且忽略重力，对N-S方程进行简化，上述的方程退化为如下的常微分方程组，边界值问题也转化成了初值问题，初始条件、求解域由嵌套圆环退化为线性I维区间套，如图I所示
d2u, I (du,、I dp/ +—=--—f A \
dr r 乂 dr j JU1 ct^ 4 jr G [r0, R]
Ci2Ugz ^ ifI dp
dr2 r Ik dr J jug dz( 5 )r G
ulz (r) =0 I r = E(6)
^- = ^ = 01 (7)
dr dr
"“厂)=〃、“))=0“( 8)对上述方程组进行积分求解，可得到环状流中气液两相在柱坐标系下轴向速度沿径向分布的速度流场解析解形式如下
权利要求
1.一种精确测量气液两相混合流体中气相体积流量和液相体积流量的方法,包括以下步骤 a.使气液两相混合流体流过竖直设置的圆管，以形成环状气液两相流；所述环状气液两相流是指这样的一种流型气体以圆柱状围绕圆管轴线，而液体以环状分布在气体与管壁之间； b.测量该气液两相混合流体的总体积流量Q、温度T、压力P和气相分率GVF； c.基于以下公式计算气液滑差因子S和气体体积流量Qg和液体体积流量Q1 S = 2-GVF/(u kX (GVF — I))Qg = Q X GVF X S/(1-GVF+S*GVF)Q1 = QX (I-GVF)/(1-GVF+S*GVF) 其中， Q为总体积流量，单位为m3/s ； S为以m/s计的气相速度与以m/s计的液相速度之比，无量纲； GVF为气相分率，无量纲； U k为在气液两相混合流的测量温度T和测量压力P下以任何粘度单位例如Pa s计的气体粘度U g与以任何粘度单位例如Pa s计的液体粘度U I之比，无量纲。
2.—种精确测量气液两相混合流体中气相体积流量和液相体积流量的方法,包括以下步骤 a.使气液两相混合流体流过竖直设置的圆管，以形成环状气液两相流；所述环状气液两相流是指这样的一种流型气体以圆柱状围绕圆管轴线，而液体以环状分布在气体与管壁之间； b.测量该气液两相混合流体中的液相平均速度V1、温度T、压力P和气相分率GVF； c.基于以下公式计算气液滑差因子S和气体体积流量Qg和液体体积流量Q1 S = 2-GVF/(u kX (GVF — I)) Qg = V1XaXGVFXS Q1 = V1XAX (I-GVF) 其中， V1为液相平均速度，单位为m/s ； S为以m/s计的气相速度与以m/s计的液相速度之比，无量纲； GVF为气相分率，无量纲； U k为在气液两相混合流的测量温度T和测量压力P下以任何粘度单位例如Pa s计的气体粘度U g与以任何粘度单位例如Pa s计的液体粘度U I之比，无量纲； A为圆管的横截面积，单位为m2。
3.—种精确测量气液两相混合流体中气相体积流量和液相体积流量的方法,包括以下步骤 a.使气液两相混合流体流过竖直设置的圆管，以形成环状气液两相流；所述环状气液两相流是指这样的一种流型气体以圆柱状围绕圆管轴线，而液体以环状分布在气体与管壁之间； b.测量该气液两相混合流体中的气相平均速度Vg、温度T、压力P和气相分率GVF；C.基于以下公式计算气液滑差因子S和气体体积流量Qg和液体体积流量Q1 S = 2-GVF/(u kX (GVF — I)) Qg = VgXAXGVFQ1 = VgXAX (I-GVF)/S 其中， Vg为气相平均速度，单位为m/s S为以m/s计的气相速度与以m/s计的液相速度之比，无量纲； GVF为气相分率，无量纲； U k为在气液两相混合流的测量温度T和测量压力P下以任何粘度单位例如Pa s计的气体粘度U g与以任何粘度单位例如Pa s计的液体粘度U I之比，无量纲； A为圆管的横截面积，单位为m2。
4.权利要求I一 3中任一项的方法，其中使用相分率计来在线测量所述相分率GVF，所述相分率计选自单能伽马射线相分率计或双能伽马射线相分率计。
5.权利要求I的方法，其中使用文丘里流量计、孔板流量计、转子流量计来测量所述总体积流量Q。
6.权利要求2的方法，其中使用互相关法来测量液相平均速度
7.权利要求3的方法，其中使用互相关法来测量气相平均速度Vg。
8.权利要求I一 3中任一项的方法，其中所述气液两相混合流体以大于15m/s的表观气速流过所述竖直设置的圆管。
9.一种精确测量气液两相混合流体中气相体积流量和液相体积流量的方法,包括以下步骤 a.使气液两相混合流体流过竖直设置的圆管，以形成环雾状气液两相流；所述环雾状气液两相流是指这样的一种流型其中悬浮有液滴的气体以圆柱状围绕圆管轴线，而液体以环状分布在气体与管壁之间； b.测量该气液两相混合流体的总体积流量Q、温度T、压力P和气相分率GVF； c.基于以下公式计算气液滑差因子S和气体体积流量Qg和液体体积流量Q全文摘要
本发明涉及一种用于环状气液两相流的流量计量方法。此方法首先使气液两相混合流体流过竖直设置的圆管，以形成环状气液两相流。然后测量该气液两相混合流体的总流量或液相平均速度或气相平均速度，并测量温度、压力、气相分率等物理量，最后基于气液滑差因子的解析解来精确计算出液体积流量和气体积流量。
文档编号G01F1/00GK102749104SQ20121025786
公开日2012年10月24日 申请日期2012年7月24日 优先权日2012年7月24日
发明者陈继革 申请人:兰州海默科技股份有限公司