气体非牛顿流体两相流动

1、气体/非牛顿流体两相流动与气/水两相流动阻力特性的比较劳力云吴应湘郑之初李东晖石在虹（中国科学院力学研究所应用流体力学研究室）提要本文基于常用的均相流动模型， 主要针对油/气/水混输管道中经常遇到的气体 /藉律液体水平管道两相流动，对其阻力特性进行了阐述，其中包括流动基本模型和方程、 常用的经验估算方法等方面，将它与气/水两相流动中对应的有关问题进行了比较。Summary In this paper, the pressure loss characteristics of the gas and power-law liquid two-phase flow thatwas often me

2、t within oil/gas/water concurrent transportation pipeline were presented. Based the homogeneous model-oneof the frequently adapted models, the issues of fundamental model and equations of the flow and the estimation methodsof pressure losses were included. These results were also compared with the c

3、orrespondents related to gas/watertwo-phase flows.关键词均相模型气液两相流非牛顿流体Key words homogeneous model, gas/liquid two-phase flow, non-Newtonian fluid1引言在海洋石油资源开发中，油气混输正在成为一种越来越重要的油气输运手段。 为了了解油气混输管道中的能耗，进行混输管线与增压系统的合理设计，在油气混输过程中对油、气、水等介质进行有效准确的计量， 以及对整个混输系统的运行状态进行监控等等，需要深入了解油气混输管道的流动阻力特性。另一方面，油气混输问题实际上属于气液两

4、相流动中的气体/非牛顿液体的流动问题。而目前在国内外，气液两相流研究主要集中于对气体/牛顿液体（尤其是气/水两相流）的流动特性上，对于气体/非牛顿液体流动问题的研究相对较少，尚无比较成熟的结果。因此，对气体/非牛顿液体的流动阻力特性进行研究，不但具有重要的工程价值，也具有较大的学术意义。实际测试结果表明，我国油田所生产的原油，在常温输送时，其流变特性表现为藉率流（du V E = K I体的特性一一流动剪应力与剪应变呈现为藉率关系（对于管流即I激）1 2。其主要的参数为稠度系数K和藉率指数n。本文中的非牛顿液体即指这一类流体。 下面根据两相 流的均相流动模型和基本方程，针对油/气/水混输管道中

5、经常遇到的气体 /藉律液体水平两相 流动的阻力特性进行了阐述，并讨论了对摩擦系数的经验估算方法，将它与气 /水两相流动 中对应的有关问题进行了比较。2均相流动模型在某些两相流态中，如果一相均匀地弥散在另一相中，两相间的动量传递和能量传递足够快，两相的局部平均速度和温度就基本相同。此时，可以应用均相模型对两相流的流动特性进行研究。均相模型的基本思路是用一等效可压缩流体代替两相流体。一般，在研究气泡流、雾状流等流态下的流动特性时，常采用均相流模型。均相模型的基本假设为：(1)( 1)气液两相具有相同的线速度，即UG=UL=UH(式中u表示介质线速度, 下标G, L和H分别表示气相、液相和均相，下同

6、) ；(2)(2)两相间处于热平衡状态。对于两相管流，设其轴线方向为z,根据连续性、动量守恒和能量守恒原理以及上述假 设，可以得到：dHUH-dz一duH-1 dA -H- :却UH0A dzdz(1)dp P2dvH-.0 G2 :gs iEdz Adz(2)dp r dF2dvH、.-=：-H-G - - ：Hg s i ndz dzdz(3)上列式中，P为介质密度，A为流道截面积，p为管道中的压力，Pr为流道周长，工。为管1_道壁面剪应力，G=PHUH为质量流速，P为介质比容，g为重力加速度，0为管道轴线的倾角，F为管道中摩擦力所引起的能量损失。比较(2)、(3)两式，易知由摩擦力引起的

7、管道压降、一c 2dpF= Pr4品HE J /% = De2dFdz(4)式中fTP两相摩擦系数，De为两相流动水力直径。dpF实际上，(2)式反映了这样一个事实， 管道的总压降由管道摩擦压降dzA0(10)一业=G2dvH一 虹=管道加速压降VdzJ dz和管道重力压降 【dzJPHgsin8三部分构成的。流体的流变特性不影响其中的加速压降和重力压降，但将通过影响两相摩擦系数而直接影响摩擦压降。3气/水两相流的摩擦系数计算均相模型在本质上将两相流动视为单相流动，其摩擦系数的确定方法与单相流动类似。在单相流动的摩擦阻力计算中，摩擦系数f与流体介质动力粘度 卜的关系比较简单，主要R GD决定于

8、流动的Reynolds数（式中D为单相流动的水力直径）。当Re低于2000时（层流），摩擦系数为f =16 / Re（5）当Re高于2000时（湍流），则应用Blasius公式：.0.25f = 0.079 Re对于气体/牛顿液体的气液两相流，常采用以下三种方法计算粘度3 :（1）取两相粘度为单相粘度。当含气率较低时，取两相粘度为液相粘度TP当含液率较低时，取两相粘度为气相粘度TP(8)1x1 - x- =-+-P-R口McAdams :TP G LCicchitti :*TP=*L(1-x)+ 叩Dukler:=人(1一)七-(2)低浓度悬浮小球粘度计算式。对于空隙率a低于0.05的气泡流，

9、取两相粘度为TP-L(3)(3)经验拟合关系式。常用下列几种上列式中，x为质量含气率（干度）,P为体积含气率。(11)(12)4气体/籍率流体流动的摩擦系数计算对于气体/藉率流体两相流，在进行摩擦系数的计算时，其基本思路类似于气/牛顿液体, 但是应根据液相的流变特性进行必要的修正。藉率液相的Renoalds数可以采用下式计算4：n 2 _p DuLReL=-,3n十1 V(13)其临界Renoalds数(Re)c也与藉率指数有关:610 (2n - 1)(3n - 2)( 4n - 3)(Re)c =-3-(3n - 1)(14)在气相层流、气相紊流和液相层流时，可以参照牛顿流体的Blasiu

10、s公式计算单相状态摩擦系数；对于液相紊流时，则采用修正的Blasius公式计算液相摩擦系数：(15)0.1330.217甘中 a =0.3164 n ,b =0.25n获得气、液两相的摩擦系数后，可以根据fTP =fL(1-x)+fGX计算出两相流摩擦系数fTP。而实际的均相摩擦系数还应进行修正。文献2根据三种管径的泡状流提出以下 的修正公式：扁=AfTP(16)参数A、B仅与藉率指数有关,J2.0255A =0.3804 n_0.8 4 4 8B =0.89 1 n需要指出，上述计算中，涉及介质的速度时应取介质的表观速度。5结论在气体/藉率液体两相流各种流态中，气泡流和雾状流流态表现出较好的介质均布性，如各参数沿流道的变化率不大，可以采用均相流模型。 对于流动阻力中的重力分量和加速分量，可以通过管道尺寸以及流动工况进行直接计算；对于流动摩擦阻力，需要根据实际数据fLbReL(17)(18)回归出经验关系进行估算。与气 /水两相流相比较，