燃气管道泄漏流量的计算.pdf

书书书 燃气管道泄漏流量的计算 黄小美 彭世尼 徐海东 杨茂华 重庆大学 城市建设与环境工程学院 重庆 潍坊港华燃气有限公司 山东 潍坊 重庆燃气设计研究院有限责任公司 重庆 摘 要 根据一元气体流动基本方程式 推导了孔口泄漏在绝热过程下的小孔模型和适合管 道完全断裂的多变过程泄漏流量计算的管道模型 联合两种模型计算任何泄漏孔口直径下的泄漏 流量 讨论了燃气最大泄漏流量的限制 进行了实例计算并对比了不同模型的计算结果 关键词 泄漏流量计算 管道模型 小孔模型 管道小孔综合模型 流量限制 中图分类号 文献标识码 文章编号 概述 在燃气管道事故定量风险评价 事故抢险预案 制定和漏气损失评估时 首先要计算泄漏流量 燃 气管道在事故破损时 燃气可通过两种途径进入到 大气中 一种是燃气直接泄漏到大气环境中 另一种 是泄漏到土壤中 通过土壤渗透进入大气环境 前 者可以通过理论推导得出泄漏流量的计算公式 后 者理论计算比较复杂且不确定性很大 本文主要分 析和讨论前一种情况下的泄漏流量计算 第三方破 坏是城市燃气管道泄漏的主要原因之一 其主要表 现是挖掘机器 钻孔机器破坏管道 在这种情况下 燃气通常直接泄漏到大气中 此外 架空管道泄漏 也是直接泄漏到大气中 小孔模型的推导 管道泄漏示意图见图 小孔模型是将泄漏孔 口当作孔径很小的小孔 从而建立泄漏流量计算的 模型 图 管道泄漏示意图 图中 点 管道起点 点 泄漏口入口点 点 泄漏口出口截面上的点 点 点 上游附近的某点 泄漏点至管道起点的距离 泄漏点上游管道体积流量 泄漏体积流量 图 中 点 通常为该管道上游的调压器出口 第 卷 第 期 年 月 煤 气 与 热 力 其压力通常保持不变 假设点 的断面流量及其平 均流速方向不受泄漏影响 而点 下游至泄漏口处 的任何点管道断面平均流速由于受到泄漏影响而不 再沿管道轴线方向 点 至点 的距离非常小 可以 忽略不计 因而点 的压力近似等于点 的压力 小孔模型假设管内燃气全部从该小孔泄漏 即 管道上游无支管或支管燃气流量为 这样假设是 为了保证从小孔泄漏的燃气流量是最大值 由于泄 漏小孔孔径较小 泄漏流量有限 因而忽略管道沿程 阻力 认为泄漏处的管内压力等于管道起点压力 即 式中 图 中点 的绝对压力 图 中点 的绝对压力 在泄漏孔处 燃气流速一般较快 燃气没有足够 的时间与环境进行热量交换 因此燃气泄漏过程 即 从点 到点 的燃气流动过程可被视为可压缩气体 绝热流动过程 可见泄漏孔口与喷嘴相似 孔口泄 漏瞬间的流动可以看作是一维流动 气体的一元流 动欧拉运动微分方程为 式中 燃气绝对压力 燃气密度 燃气断面平均流速 因为泄漏过程为绝热过程 所以有 式中 燃气的等熵指数 常量 是温度的函数 在常温下理想气体的 可近 似当作定值 对于天然气等由多原子分子组成的 气体 取 因为城市燃气压力大多在 以下 燃气 温度为常温 所以燃气可以看作理想气体 因此 有 式中 燃气的气体常数 燃气的温度 将式 代入式 并在小孔入口点 和小孔 出口点 进行积分 然后将式 代入整理可得 式中 点 的燃气断面平均流速 即为燃气 的泄漏出口流速 点 沿管道泄漏口轴线方向的流速 点 处的燃气温度 点 的绝对压力 可按式 计算 式中 环境压力 临界压力 临界压力比 在小孔模型下 一般很小或近似为 因而 远远小于 可以忽略不计 因此式 可简化 为 槡 泄漏孔出口压力与入口压力的比值等于临界压 力比时的泄漏出口燃气流速称为临界流速 将式 代入式 有 槡 式中 泄漏出口燃气临界流速 当 达到临界压力 后 点 的流速 就等于 临界流速 根据孔口出流的质量流量公式 可得出泄漏质 量流量为 式中 泄漏质量流量 流量系数 可取 泄漏孔口当量直径 点 处的燃气密度 不规则孔口当量直径按下式计算 槡 第 卷 第 期 煤 气 与 热 力 式中 泄漏孔口面积 将式 代入式 并整理可 得 当 时 槡 槡 当 时 槡 槡 在小孔模型下 认为式 成立 且由于管内流 速小 管内流动可视为等温过程 因而对于任一管 道 只要知道管道起点压力 管内燃气温度和泄漏孔 口当量直径 便可按照式 和 计算泄漏质量 流量 管道模型的推导 管道模型适合于管道完全断裂的情形 即泄 漏当量直径等于管道内径 点 即为管道末端 点 和点 重合 管道泄漏流量就等于管输流量 此时可 按管道水力计算公式来计算管道流量 燃气泄漏过程中 管道上游阀门关闭之前 管内 燃气流动可以视为稳定流动 在燃气管道完全断裂 泄漏情况下 燃气流速较大 管内燃气可能没有足够 的时间和周边环境进行充分的热交换 管内流动也 不能看作等温过程 而只能看作多变过程 因而对于 管内燃气流动过程有 式中 多变指数 常量 仍然假设从管道起点至泄漏点之间的管段上没 有支管 或所有支管的流量都为 同时考虑到管内 燃气流速较大时 不能忽略对流项 将文献 提 供的管内流动基本微分方程组联合式 可以推 导得到式 式中 泄漏点上游管道内燃气的质量流 量 管道内径 摩擦阻力系数 可按文献 或 提 供的公式计算 式 中 时 表示管内流动为等温过程 此时认为管内燃气与周边环境有充分的热交换 管 内温度等于环境温度且保持不变 当 时 表示 管内流动为绝热过程 此时认为管内流速太快或者 管道太短 管内燃气完全没有和环境进行热交换 实 际上这两种理想状况都不存在 的值在 和 之 间 为了简化计算 通常在管内燃气流速较小 管道 较长时 取 而在管道流速很大或者管道很短时 取 在管道模型下 点 为管道末端且暴露于大气 环境中 因而泄漏口处的压力等于大气环境压力 即 此时只要知道管径 泄漏处至管道起点的距离 管道起点的压力 选取适当的值 便可按照式 计算管道的质量流量 管道 小孔综合模型 当燃气管道完全断裂时 按管道模型计算泄漏 流量比较准确 当燃气管道只有孔径很小的破损孔 时 按小孔模型计算泄漏流量比较准确 但在实际 情况下 特别是由于施工开挖导致的断裂 泄漏口既 不是小孔 也不是完全断裂 因此用这两种模型都不 准确 管道 小孔综合模型就是将管道模型和小孔 模型结合起来进行泄漏流量计算 事实上 在小孔模型下 只要有燃气泄漏 管道 内燃气就会流动 由于摩擦的存在 管道内必然会有 沿程阻力 因此式 并不成立 实际的 还需要根 据式 计算 前面已经假设管输燃气全部从小 孔泄漏 这种假设对于风险评价和事故应急抢险来 说是保守的 因而有 式 中 可根据式 和 计算 式 和 中的 可由式 和 联立求得 可以根据式 计算 因此联立式 便可计算出泄漏质量流量 泄漏处的管内压力 及管内温度 该模型对 于任何泄漏口尺寸的稳定泄漏过程都比较准确 计 黄小美 等 燃气管道泄漏流量的计算 第 卷 第 期 算时需要利用计算机语言编程求解 为了符合工程习惯 通常需要将泄漏质量流量 转换成泄漏体积流量 式中 标准状态下燃气密度 最大泄漏流量限制的讨论 燃气在管道内的最大流量有限制 并不是按照 管道模型或者管道 小孔模型计算得到的结果就一 定是实际流量 而必须对计算结果进行两个方面的 校验 一是调压器最大通过能力的限制 二是管内临 界流量的限制 调压器最大通过能力的限制 城市燃气管道都与上游调压器相连 管道内的 流量和泄漏流量都不能大于上游调压器的最大通过 能力 调压器说明书通常都给出了该调压器在 各进口压力下的 值 因此无论用哪种模型计 算 泄漏流量的计算结果最大只能取 管内燃气临界流量的限制 与泄漏孔处的绝热流动相似 燃气管内流动也 存在临界流量问题 由于沿程摩擦阻力对燃气运动 参数的影响 实际气体一元运动微分方程为 将式 代入式 参照文献 推导可得 式中 马赫数 马赫数按式 计算 槡 式中 当地声速 管道起始端流速一般较小 马赫数远小于 因 而有 由式 和 可知 随着流动距离的增长 压力不断减小 速度不断增加 但压力不可能无限降 低 流速不可能无限增大 即 不可能为 因此式 恒成立 槡 又根据管内稳定流动质量流量方程有 将式 代入式 整理可得 槡 在等内径稳定管道流动中 式 中参数除断 面平均流速 外均为常数 因而马赫数随 增大而 增大 图 中点 处马赫数最大 有 槡 式中 点 处的马赫数 将式 代入式 整理可得 式中 点 处的燃气绝对压力 点 处的燃气温度 根据图 中的假设有 将式 代入式 并转化成体积流量 则 有 槡 槡 令管内临界体积流量 为 槡 槡 式中 管内燃气临界体积流量 无论按照何种模型计算所得的体积流量必须校 验是否满足式 若计算得到的泄漏体积流量大 于 则只能取 实例计算 某管道受施工破坏而连续泄漏 泄漏口近似为 圆形 该管道上游调压器出口绝对压力为 调压器最大通过能力为 泄漏点 距离调压器 管道内径为 已知天然 气密度为 管道周边环境温度为 分析在不同泄漏当量直径下的泄漏流量 最大流量限制 在管道完全断裂时泄漏流量最大 此时 由于燃气流速较大 假设管内流动为绝 第 卷 第 期 煤 气 与 热 力 热过程 取 将 代入式 计算得到管道 内燃气临界体积流量为 而该管道上游调压器最大通过能力为 因此最大泄漏体积流量限制为 管道模型计算 管道完全断裂情况下 根据式 编程计算 可得燃气泄漏的体积流量为 小于最大泄漏体积流量限制值 因此该值为 管道在此位置的最大可能泄漏体积流量 管道 小孔综合模型与小孔模型计算比较 采用管道 小孔综合模型时 可联立式 编程计算 其中 按 式 确定 式 假设多变指数与泄漏孔口面积成线性 关系 当管道完全断裂时 管内流动为绝热过程 而 当泄漏孔口非常小的时候 为等温过程 经计算 泄漏孔口当量直径从 增加到 管道 小孔综合模型泄漏体积流量 泄漏处 的管内压力 泄漏处的管内温度 随泄漏孔口当 量直径的变化规律分别见图 中的曲线 小 孔模型计算结果则分别见图 中的曲线 图 泄漏体积流量与泄漏孔口当量直径的关系 不同泄漏位置最大泄漏体积流量计算 本文推导的管道模型与文献 或 中提供 的管道水力计算公式相比 考虑了管内流动的温度 变化和对流项的影响 为了比较这两个公式计算结 果的差别 在本例中假设泄漏当量直径等于管道内 径 分别用两个公式计算管道不同泄漏位置的泄漏 体积流量 结果对比见图 当泄漏处至管道起点 较远时 两者计算结果基本一致 前者计算结果比 后者略高 而当泄漏处距管道起点较近时 前者计 图 泄漏处管内压力与泄漏孔口当量直径的关系 图 泄漏处管内温度与泄漏孔口当量直径的关系 算结果比后者低 曲线 前段为直线是由于考虑了 最大流量的限制而形成的 图 泄漏点位置与该点最大泄漏体积流量的关系 结论及建议 燃气管道泄漏流量计算是燃气管道风险评价的 重要环节 也是事故应急抢险的重要依据 利用合理 假设推导的管道 小孔综合模型计算泄漏流量 同 时考虑泄漏流量的限制条件 能够得到比较精确的 结果 能够满足风险评价和应急抢险的需要 本文的模型没有考虑非稳定过程 实际工程中 在管道开始泄漏瞬间和当泄漏管道的上游阀门被关 闭后 燃气流动是非稳定流动 但在燃气持续泄漏期 间 燃气流动可看作稳定流动 在风险评价和事故应 急抢险过程中 往往关心的是这种持续泄漏的影响 黄小美 等 燃气管道泄漏流量的计算 第 卷 第 期 因此可不必进行复杂的非稳定流动计算 若需要更 为精确的计算结果 则需要建立管内燃气与周边环 境之间的传热关系 建立非稳定流动和传热方程组 确定边界条件 采用计算流体力学方法计算 参考文献 周谟仁 流体力学 泵与风机 第 版 北京 中国建筑工业出版社 邱信立 廉乐明 李力能 工程热力学 第 版 北京 中国建筑工业出版社 段常贵 燃气输配 第 版 北京 中国建筑工 业出版社 城镇燃气设计规范 作者简介 黄