第八章瞬变流1

输油管道设计与管理第八章输油管道瞬变流第一节管道的瞬变流动过程第二节管道瞬变流动过程的数学描述与特征线解法第三节输油管道中的水力瞬变流动第四节输油管道瞬变流动过程的控制稳定流动：压力管道中任一点的流速和压力仅与该点的位置有关，而与时间无关的流动称为稳定流动。不稳定流动瞬变流动（Transientflow）：流体从一种稳定流动状态过渡到另一种新的稳定流动状态的过程。水击（WaterHammer）：又称“水锤”。有压管道中因流速、压力急剧变化而引起压力波在水中沿管道传播的现象。第一节管道的瞬变流动过程水击现象（Water-hammerPhenomena）：在有压管道系统中，由于某一管路中的工作状态的突然改变，引起管内液体流速的急剧变化，同时液体压强大幅度波动，这种现象称为水击现象。正水击（PositiveWater-hammer）：当管道阀门迅速关闭时，管中流速迅速减小，压强显著增大，这种水击称为正水击。负水击（SuctionWater-hammer）：当管道阀门迅速开启时，管中流速迅速增大，压强显著减小，这种水击称为负水击。第一节管道的瞬变流动过程一、水击的研究历史1897年，茹可夫斯基（俄）发表了他的经典报告，全面而正确的阐述了水击的机理，给出了惯性水击压力和水击波传播速度的计算公式。1902年，阿列维（意）建立了不稳定流动的微分方程，奠定了水击分析的理论基础。20世纪60年代以后，借助于电子计算机的发展，数值法求解阿列维方程成为可能。斯特里特（美）1963年创建的特征线法，开创了水击分析的新局面。二、水击产生的原因—流量突变计划地调整输量或改变输送流程开关阀门：阀门最后关闭的0.1～0.3秒，水击最强。油品的注入与分输泵站节流、调速、启停泵改变输送压力、流量油品切换混油通过中间泵站二、水击产生的原因—流量突变事故引起的流量变化停电阀门失灵管道泄漏三、水击形成过程水击过程四个阶段：增压逆波减压顺波减压逆波增压顺波四、管道瞬变流的基本公式直接瞬变压力公式根据儒可夫斯基的水击理论，由流速瞬间变化直接产生的压力脉动值，可用下式计算如果阀门阶段性关闭，则四、管道瞬变流的基本公式压力波传播速度管内压力波的传播速度取决于液体的可压缩性和管子的弹性。含气液体中的波速计算（略）五、水击危害1、正水击—超压，管道破裂2、负水击—管道失稳、液柱分离、气体逸出3、振动六、水击分析的意义在优化设计中的应用在编制运行方案中的应用预测混油界面及检漏在操作培训中的应用七、水击的控制对于管道瞬变流动过程，控制的目的是：避免管道超压（包括超高压和超低压）；减轻管道运行参数的脉动，维持管道的平稳运行。用于控制管内瞬变流动过程的装置和措施很多，根据其作用原理可分为两类：改变流速变化过程（幅度和时间）使用各种压力保护措施

压力保护装置和措施调节阀压力保护泄压阀回流保护泵机组自动停运超前保护七、水击的控制第二节

管道瞬变流的数学描述与特征线法分析管道的瞬变流动过程，使用弹性理论比较符合客观实际。描述管道瞬变流动过程的弹性理论以液体可以压缩和管壁产生变形为基础，分析受压力波作用的微元体的受力与质量守恒，建立运动方程和连续方程的微分形式。将此两方程联立，对应相应的边界条件，是用数值方法即可分析管道的瞬变流动过程。一、管道瞬变流动过程的数学描述运动方程连续方程（一）运动方程运动方程是从处于瞬变流动的管内流体中选取隔离体，应用牛顿第二定律建立的。推导方程是有两点假设：（1）管内流体为均质一维流动，这就意味着管道横截面上的压力、流速和密度是均匀分布的；（2）管子和流体变形均在线性弹性变形范围内。

把压力P换算成液柱高度：由于：隔离体内视为均匀密度。由于隔离体足够小，（8-9）（二）连续方程连续方程描述了时刻t流入及流出隔离体的流体净质量等于隔离体内的质量随时间的变化率。对隔离体，考虑时间的质量守恒，应存在

展开，忽略高阶小量，整理后可得

由于（8-10）二、特征线解法分析液体管道不稳定流的两个基本方程一阶拟线性双曲型偏微分方程组由于二次摩擦项的存在，目前用分析方法还无法得到它们的解析解，只能用数值方法求解。求解液体管内不稳定流问题的最常用的一种数值方法——特征线解法（一）特征方程特征方程是用一对特征值，，将两个偏微分方程线性组合，即

（一）特征方程（二）有限差分方程特征方程虽然是常微分方程的形式，但由于摩阻项是非线性的，仍不能用微分方法得到解析解。下面介绍用有限差分方法求数值解。使用数值方法求解特征方程，必须把管道沿线长度方向离散成若干管段，把瞬变过程离散成若干个时间步长，在小范围内进行近似计算。沿管道（简单管道）长度等间距分成N段，每段长，称为空间步长，管道的结点数为N+1；压力波传播距离所需时间，称为时间步长。这样，把x-t平面画成网格。在一个网格范围内，沿特征线，在i-1和i之间，对方程分离变量后积分由于无法知道瞬变过程中V与x之间的函数关系，后面一项积分只能作数值计算。对于摩阻较小的管道，后面一项积分可采用一阶近似方法计算。差分方程为在进行管道瞬变流动计算时，往往用流量代替流速更方便。重新整理成下面的形式：（8-18）（8-19）带有p下标的表示当前时刻的计算值，不带p下标的表示已知参数。显然，是分别由i-1和i+1点已知的参数计算的常系数。联立求解（8-18）和（8-19）式可得对于i结点，时间终了时的参数和可由时间开始时i-1和i+1结点的参数计算求得。（8-20）（8-21）【例】站间管道如图所示，管材，长25.5km，稳态下流量Q=0.025m3/s，（V＝1.47m/s），m=0.125，f=14.13，a=1214m/s，稳态时末端进站压头为434米。末端阀门在0s时瞬间关闭，上游罐液位恒定。

0s1s2s3s4sQ=0.025502485468451434i=0.0141819202122（三）扰动的影响区沿管道不同的位置，受到扰动的影响的时间是不同的。发生扰动的位置不同，管道沿线受到扰动影响的过程也不同。以P点为起点，两特征线上方包络的范围内的状态受到影响。因此，该区间也被称为P点扰动的影响区。而两特征线下方区间管内没有受到P点扰动的影响，只取绝于初始状态，所以该区间也称为初始条件支配区。（四）