油气储运工程外文翻译应用特征方法的天然气管道破裂问题的计算.doc

中文译文 应用特征方法的天然气管道破裂问题的计算 J.A. Olorunmaiye ，机械工程系，Ilorin大学，P.M.B. 1525, Ilorin(Nigeria)) N.E. Imide，机械工程系，Abubakar Tafawa Balewa大学，Bauchi (Nigeria) （完稿于1992年7月4日，校正后定稿于1992年12月2日） 摘 要 将在长距离、高压、伴随有突然断裂的天然气管道内的流动模拟为不稳定的一维等温流动。集合的双曲线偏微分方程已由数值特征方法所解答。应用具备二次插值的线性特性时，其数值计算的准确性是满足要求的。我们也发现，等温流动中特征曲率的作用并不像Flatt的绝热流的报道中所说的那样明显。为了评估天然气管道破裂的危险，有必要了解断点处气体的外泄速率。预测的断裂口处泄漏气体的质量流量要比应用绝热流理论预算的小18%，然而，这与早些时候利用等温流动理论的工作者得出的结果却很好的吻合，其理论计算方法是基于加权残差的。 说 明 天然气管道的瞬态压力波动通常是由阀门或者压缩机的突然开启与关闭引起的。管道系统内存在压力波的传输与反射，从而会因为压力过大导致管道破裂。天然气工业对长输管道在事故破裂后的气体不稳定流动是相当感兴趣的，因为管道内含有大量的可燃气体，它的破裂也就意味着相当大的危害。为计算气体扩散范围和空气—天然气混合气的地面浓度，有必要了解气体的泄漏速率。破裂后，由于在低压段的流动条件下，膨胀波进入管道后会造成流动逆转，低压段管道内的气体流动将比高压段更复杂。在图1中，假定了破裂相当于管道截面大小，并立即向左方延伸（其高压管段），阀门在断裂的瞬间完全关闭。 Fannelsp and Ryhming [1]研究了高压下发生破裂时天然气管道内的不稳定流动。他们将流动划分为三个时间状态。紧接着破裂发生的“早期”状态受波动进程影响，并且开口端的压力接近周围环境压力。“中期”状态伴随着一个内部压力峰值的出现，大约对应于流动逆转速度为0的位置。当压力峰顶到达低压尽端（截断处）时，“后期”状态便出现了。在这一时段内，压力从关闭端到开口端单调递减。 图 1 破裂后膨胀波在管道内传播的物理环境 在假定等温流动的条件下，他们用积分法分析了管道内的不稳定流动。并且忽略惯性条件得到了一组抛物线方程。不过，他们的方法无法解决“早期”状态的问题同时在内部流动问题上也毫无建树，因为压力和压力—流速积在x=0至x=lo 上的空间分布是被假定的，其目的是获得解决方法。 Lang and Fannelsp [2] 使用了一些加权残差族中近似的程序在前人的成果上有所提高。他们的流动没有以Fannelsp and Ryhming 的工作为先验，但是却成为了解决方法的一部分。他们能够获得更准确的结果，但是该方法在“早期”状态的预测中会产生因变量的震荡。 Ryhming [3] 利用了匹配渐进展开的方法来研究长输管道突然破裂后气体的不稳定等温流动。他证明了流动处于临界状态下，管壁的摩擦使得流速顺着管道在破裂处变的异常。他的复合第一阶解给出了不同时段的流速分布图。因为他的分析集中于预算流速分布，并未考虑压力或者密度的变化，所以此方法不能预算“早期”状态的外泄速率。 Flatt[4]研发了一种特征方法来估算长输管道破裂后天然气的不稳定绝热流动。因为使用了全套的双曲线偏微分阶级方程，他的方法能够处理三个时态的流动问题。 然而，像管道破裂产生的高压变化，流体与管道周围媒介之间的热传递是相当可观的，除非管道加有防冻保温层。事实上，流动过程会更接近等温而不是绝热状态，因为密度和临界流速的补偿误差发生在前者。目前的努力都是针对于开发一种数值计算方法来解决双曲线偏微分阶级方程以描述破裂发生时的不稳定等温流动。因此，这就要向Ryhming, Fannelep 和Lang提供了“早期”状态的分析中所缺少的数据表示祝贺。它对三个时态的流动分析也同样有利。 Flatt[4]并没有将绝热模型应用到预测直到“后期”状态时的流动，因为这需要很长的计算时间。由于等温流动方程（2个）要比绝热方程（3个）简单，这样用等温方程解决问题会花费较少的时间，从而使其在计算所有时态内的流动时更经济。 数学模型 对于等径管内的一维流动，其质量和动量守恒方程为： 式中： 将u2表示为u | u|以保证摩擦力方向总是与运动方向相反的。 理想气体等温流动的状态转移方程是： 联立式（1）和（2）,消掉密度ρ，带入式（4）中，则连续性方程和动量方程变为： 和 其中，等温条件下的声速表示为 尽管等温条件下的声速是恒定的，膨胀波仍然会分散传播，如图1。因为波头在流体粒子中的传播具有正速度，然而波尾以声速c的传播相对于流体粒子来说却是负速度。 以管道长度，等温声速，环境压力和温度为参考常量