多相管流理论与案例研讨.docxVIP

气液多相管流压力计算程序编制要求：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法计算上倾管道压降clcclear allformat longfprintf(Beggs-Brill法计算上倾管道压降);%已知参数qj=pi/90;%管道倾角，radd=0.2;%管道半径，mL=500;%管长，mQl=1.2*10^(-2);%液体流量，m3/spl=810;%液体密度，kg/m3nl=5*10^(-3);%液体粘度Pa.sQg=3.4*10^(-2);%气体流量pg=6.2;%气体密度ng=1.2*10^(-5);%气体粘度bmzl=25*10^(-3);%表面张力N/mp=0.6*10^6;%平均压力MPag=9.8;%计算相关参数A=pi*d^2/4;%面积u=(Ql+Qg)/A;%混合物流速，m/susg=Qg/A;%气相折算速度，m/susl=Ql/A;%液相折算速度，m/sRL=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.01FrL1fprintf(分离流) a=0.98;b=0.4846;c=0.0868; d1=0.011; e=-3.768; f=3.539; g1=-1.614; HL0=a*RL^b/Fr^c;%水平管的截面含液率 C=(1-RL)*log(d1*RL^e*Nlw^f*Fr^g); X=1+C*(sin(1.8*qj)-(1/3)*(sin(1.8*qj))^3);%修正系数elseif RL=0.01FrL2fprintf(分离流) a=0.98;b=0.4846;c=0.0868; d=10.011; e=-3.768; f=3.539; g1=-1.614; HL0=a*RL^b/Fr^c;%水平管的截面含液率 C=(1-RL)*log(d1*RL^e*Nlw^f*Fr^g); X=1+C*(sin(1.8*qj)-(1/3)*(s