成品油管道泄漏计算新模型研究.pdfVIP

维普资讯 第 10卷 第 4期 重庆科技学院学报 自然科学版 2008年 8月 成品油管道泄漏计算新模型研究 刘广文 张刘樯z陈次昌 张 涛。董耀文。 1．中国石油管道科技研究中心，廊坊065000；2．西南石油大学，成都610500 摘 要 ：以油滴为研究对象，利用经典力学的基本理论对油滴 的运动过程和受力情况进行分析，建立 了计算成 品 油管道小孔泄漏喷射高度的数学模型，利用该模型可直接根据油品喷射高度和角度反算泄漏孔径。实例计算表 明模型准确性较高。 关键词：成品油管道；泄漏孔径；喷射高度；反算 中图分类号 ：TE832 文献标识码 ：A 文章编号 ：1673-1980 2008 04-0039-03 成品油管道 的输送压力通常较高，这是 由于管 道的腐蚀穿孔或第三方破坏 如打孔盗油等 ，会在 管道上形成泄漏小孔 ，从而造成油品的喷射泄漏 。泄 漏油品在地面和大气中扩散 ，如果遇到火源则会引 起火灾或爆炸 ，危害巨大。为了准确评估油品发生泄 漏后的扩散规律 ，从而对泄漏危害后果进行评价 ，首 先需要对油品的泄漏速率进行计算 。而考虑到安全 因素 ，在油品‘发生喷射泄漏时不可能近距离测量泄 图 1 油 品管化运动过程示意图 漏孔径，只能利用油柱的喷射高度和角度来反算其 孔径 。李大全利用紊动射流原理对成 品油管道的喷 作出如下假设翻： 射泄漏进行研究，提出了喷射高度的计算公式n】。但 1 油舌 由一群大小各异的油滴组成 ； 2 油滴大小在泄漏孔 13处已经形成 ，油滴在运 由于其中的修正系数难 以确定，因此其适用性有限。 本文所建立的模型方便于实际计算。 动过程中不发生破碎 ，并保持球形 ； 3 油滴之间不发生相互作用 。 1 模型的建立 1 ． 2 油滴受力分析 1．1 油滴运动的基本过程及基本假设 油滴在大气中运动 ，受到空气阻力、重力和空气 油品从泄漏孔 口喷出到形成各种直径的油滴落 浮力的共同作用 ．如图2所示。 到地上是一个极其复杂的过程。油品从孔 口高速喷 出，形成油舌 。部分油 品形成油滴落在孔 口近处，未 受影响的油品在初速度作用下 向前运动 。受气流反 作用而逐层剥落，当油滴表面张力不足 以抵抗外力 时，最初生成的较大直径的油滴会再次分裂，形成更 小的油滴翻。油品简化运动过程为油品从A点喷出。 达到最高点 后作平抛运动 ，最后落到 C点，在地 面形成 圆形油池 。油品不断从 C点注入油池 ，使得 油池半径不断扩大 ，其过程如图 1所示。 要精确分析上述变化过程相当困难 ．因此需要 图 2 油滴受力示意 图 收稿 日期 ：2008—03—21 作者简介 ：刘广文 1963一 ，男，河北廊坊人 ，高级工程师。 · 39 · 维普资讯 刘广文，张刘樯，陈次昌，张涛，董耀文：成品油管道泄漏计算新模型研究 1 空气阻力 l 0 根据流体力学理论 ，当运动物体 的雷诺数 e j 很大时．阻力与速度的2次方成正 比，惯性作用远大 于粘性作用 ；只有 当雷诺数 e很小时，粘性效应才 空气浮力在 、Z轴上的投影分别为 ： 会大于惯性效应 ，阻力与速度 的 1次方成正 比。一 l 0 般情况下 ．阻力与速度 的 1～2次方成正 比。 喷射油舌的雷诺数都很大 ，处于阻力平方 区，因 将式 6 ～ 8 代入式 5 并化简得到X和 Y方 此空气阻力可用牛顿定律计算： 1 ’ 向的加速度 ： FK 1CDpSV 1 式 中： 为空气 阻力 ，N；C。为惯性 阻力系数 ，取 一 净 0．44；p为空气密度 ，kg／m ；V为 t时刻油滴运动速 _-_3CDP‘a 度，m]s；S为油滴投影面积，m2；S 1 ，其 中do为 厶 1．4 喷射高度的计算 油滴直径 ，由于没有关于油滴直径的计算公式 ，此处 由于式 9 是二阶微分方程 ，很难求得解析解 ， 近似采用水滴直径的计算公式，由式 1 计算3[1。 ng2o699 —0203 R’ 只能求其数值解。将时间离散，利用差分格式得到： 0．0699282D H 2 、-一 式 中：D为泄漏孔 口孔径 ，mm；日为油品压力 ，kPa。 Vx ￡ Vxo-At22ax t 2 重力 G 1 Vz ￡ Vzo-At∑ ￡ G 1T g 3 式中： 为油滴直径，m；p。为油滴密度 ，kg／m ；g为重 f ：、／ f +2 f 1o 力加速度 ，9．8m／s。 O, arctan- 3 空气浮力 1 1 1 g 4 H, Ho+At∑Vz ￡ 式中，P为空气密度 ，k m。 Vz t 10 1．3 油滴运动方程 式中：Vx t 为 t时刻油滴竖直方 向分速度 ，rrds；V