含气输油管道的不稳定流动.docxVIP

含气输油管道的不稳定流动 油挤出道中含有一定的轻组分。如果管道输送压力降低到一定值，轻组分将消失。逸出的气体在管壁上或在油品中的悬浮颗粒周围形成气泡。含气输油管道不稳定流动时,一方面可能会随着空泡增大,使管道压瘪;另一方面可能会随着压力升高,空泡溃灭,使管道破坏。含有气泡的输油管道不稳定流动与不含气泡的输油管道的不稳定流动是有区别的,特别是当气泡集聚在一起产生液柱分离时,如果随着空泡进一步增大,就会使管道压瘪;或者随着压力再度升高时,空泡发生溃灭,使管内压力发生跃升。因此,研究含气管道不稳定流动具有重要意义。 一、 气体密度l 若液体中含有少量的气泡,在进行波速计算时,可以忽略管壁弹性的影响,此时波速公式可表示为: a=√Κ/ρ(1)a=K/ρ????√(1) 式中a——流体波速,m/s; ρ——流体密度,kg/m3; K——流体体积弹性模量,Pa。 式(1)中应该考虑气泡对液体密度ρ和弹性模量K的影响。 ρ=ρg?g?+ρl?l?(2)ρ=ρg?g?+ρl?l?(2) Κ=-ΔΡΔ?/?=-ΔΡΔ?l+Δ?g?(3)K=?ΔPΔ?/?=?ΔPΔ?l+Δ?g?(3) 设: Κl=-ΔΡΔ?l/?l?Κg=-ΔΡΔ?g/?(4)Κ=11Κl?l?+1Κg?g?=Κl1+?g?(ΚlΚg-1)(5)a=√Κl/ρ1+ΚlDEδ+mRΤΡ(ΚlΚg-1)(6)Kl=?ΔPΔ?l/?l?Kg=?ΔPΔ?g/?(4)K=11Kl?l?+1Kg?g?=Kl1+?g?(KlKg?1)(5)a=Kl/ρ1+KlDEδ+mRTP(KlKg?1)???????????????√(6) 式中ρg——气体密度,kg/m3; ρl——液体密度,kg/m3; ?g——单位管长内气相体积,m3; ?l——单位管长内液相体积,m3; ? ——单位管长内流体体积,m3; Δ? ——单位管长内流体体积变化率,m3; Δ?g——单位管长内气相体积变化率,m3; Δ?l——单位管长内液相体积变化率,m3; Kg——气体体积弹性模量,Pa; Kl——液体体积弹性模量,Pa; E——管材体积弹性模量,Pa; D——管子内径,m; δ——管子壁厚,m; P——流体压力,Pa;T——流体温度,K; m——单位体积内气体摩尔数,kmol/m3; R——气体常数,取8 314.3 J/(kmol·K)。 若认为气体压缩过程是等温过程,则有Kg=P,同时,由于ΚlΚg?1,式(6)可简化为: a=√Κl/ρ1+ΚlDEδ+mRΤΚlΡ2(7) 或 a=a0√1+C1m/Ρ2(8)a0=√Κl/ρl1+ΚlDEδ(9)C1=RΤΚl1+ΚlD/Eδ(10) 式中a0——没有自由气体时的液体波速,m/s。 若以绝对压头H代替压力P,即P=ρgH,则式(8)可以表示为: 二、 气泡色散配置模型 1、 等温规律 气泡在液体中分布是一个复杂问题,为了计算方便,气泡离散布置模型将气泡集中在各计算截面上,且其压力和体积的变化遵循等温规律。而计算截面之间的液体则假定是纯液体,不含任何自由气体。文献在介绍计算方法时,没有考虑波速的变化,为此本研究在这一方面作了改进。 2、 数学模型及其处理 (1) 计算变量及方程 一维不稳定流动的基本方程可描述为: ?U?t+B?U?x=C(12) 式(12)中,U为自变量列向量,U=[H,ν]T,系数矩阵B及列向量C表达式分别为: B=(υa2ggυ)?C=(υsinθ-fυ|υ|2D)(13) 式中υ——流体速度,m/s; H——流体压头,m; g——重力加速度,m/s2; θ——管道与水平面之间的夹角,rad; f——管道摩阻系数。 (2) 节点处气泡体积随加权因子的变化 对于方程,令|B-λE| =0,其中E为单位矩阵。解得方程组的特征根为: {λ1=υ+aλ2=υ-a(14) 因此,可得到方程的两条特征线和相容方程为: 将波速方程式(11)代入特征方程式(15)和相容方程式(16)得到: 按照模型的假设,气泡是集中在节点上,但气泡大小会随时发生变化的,因此描述其变化规律的方程应该作为必须满足的条件予以考虑,即应将每一节点作为内边界来处理。 按图1,沿着C+及C-特征线积分相容性方程,应分别有: g(√C2+Η2Ρ+√C2+Η2R)a0(ΗΡ+ΗR)ΗΡ-ΗRΔt+QΡu-QRAΔt+f|QΡu+QR|(QΡu+QR)8Ad-g(√C2+Η2Ρ+√C2+Η2R)a0(ΗΡ+ΗR)QΡu+QRAsinθ=0(19)g(√C2+Η2Ρ+√C2+Η2S)a0(ΗΡ+ΗS)ΗΡ-ΗSΔt-QΡd-QSAΔt-f|QΡu+QS|(QΡu+QS)8Ad-g(√C2+Η2Ρ+√C2+Η2S)a0(ΗΡ+ΗS)QΡu+QSAsinθ=0(20) 节点处气泡体积的变化应满