多相管流理论与案例概述.docx

PAGE   气液多相管流压力计算程序编制 要 求： 1.按专业做与本专业相关的多相管流压力计算程序； 2.写出计算原理及计算方法，手算一段压力； 3.写出程序编制语言、运行环境、程序功能及程序框图； 4.程序编写及结果输出。 1.压降计算模型（Beggs-Brill方法） 1.1Beggs-Brill方法原理 Beggs-Brill方法是可用于水平、垂直和任意倾斜气液两相管流动计算的方法。从能量守恒方程出发，推导了考虑起伏影响的两相管路压降梯度计算式。 在稳定流动下，建立气液混合物沿管路流动时的能量守恒方程： （1-1） 式中各项分别代表单位质量气液混合物的压降梯度、单位管长动能、位能变化和摩阻损失之和。 由于： （1-2） 由公式（1-1）（1-2）得出 （1-3） 式中：-管段内流动介质的平均压力（绝对）， -管段长度， -截面含液率； -两相混输水力摩阻系数； -液相密度，； -气相密度，； -气液混合物的质量流量，； -气液混合物速度，； -气相折算速度，； -管道内径，； -管段倾角，度或弧度 -重力加速度，。 1.2截面含液率的计算 引入修正系数，表示倾斜管截面含液率与水平管截面含液率之比值，即 (1-4) 式中 -倾角为时的截面含液率 -水平管的截面含液率 1.2.1的计算 ， （1-5） 其中： -体积含液率 -液相速度准数 -富劳德准数 -与流型有关的系数，其值见表1-1. 对于过渡流型，则先分别用分离流和间歇流计算出之后采取内插法确定其持液率。 ， 1.2.2 水平管截面含液率的计算 ， （1-6） 其中-数值取决于流型其值见表1-2 表1-2 流型  分离流0.98 0.48460.0868间歇流 0.845 0.53510.0173 分散流1.0650.59290.0609 , ，式中下表分别表示过渡流、分离流和间歇流。 1.3水力摩阻系数计算 （1-7） 其中 1.4贝格斯-布里尔法判断流型见表1-3 表1-3 流型判别准则分离流过渡流间歇图分散流 1.5计算步骤 2计算实例 一上倾的油气混输管路，已知倾角，管内径，长，在管路平均压力和温度下，原油流量，密度，粘度，天然气流量，密度，粘度。原油的表面张力，管路平均压力。计算倾斜管压降梯度。 MATLAB运行结果： Beggs-Brill法计算上倾管道压降间歇流 dfdl =1.517185410566795e+002 用MATLAB软件编的程序如下 %Beggs-Brill法计算上倾管道压降 clc clear all format long fprintf(Beggs-Brill法计算上倾管道压降); %已知参数 qj=pi/90;%管道倾角，rad d=0.2;%管道半径，m L=500;%管长，m Ql=1.2*10^(-2);%液体流量，m3/s pl=810;%液体密度，kg/m3 nl=5*10^(-3);%液体粘度Pa.s Qg=3.4*10^(-2);%气体流量 pg=6.2;%气体密度 ng=1.2*10^(-5);%气体粘度 bmzl=25*10^(-3);%表面张力N/m p=0.6*10^6;%平均压力MPa g=9.8; %计算相关参数 A=pi*d^2/4;%面积 u=(Ql+Qg)/A;%混合物流速，m/s usg=Qg/A;%气相折算速度，m/s usl=Ql/A;%液相折算速度，m/s RL=Ql/(Ql+Qg);%体积含液率 nd=ng*(1-RL)+nl*RL;%两相混合物粘度 pf=RL*pl+(1-RL)*pg;%两相混合物流动密度 Fr=u^2/g/d;%富劳德准数 Nlw=usl*(pl/g/bmzl)^(1/4); L1=316*RL^0.302; L2=9.252*10^(-4)*RL^(-2.4684); L3=0.10*RL^(-1.4516); L4=0.5*RL^(-6.738); %流型判断，计算截面含液率 if RL0.01FrL1 fprintf(分离流) a=0.98;b=0.4846;c=0.0868; d1=0.011; e=-3.768