投递文章
投稿指南赵新华1，陈春芳1,郑毅2
(1．天津大学 环境工程系，天津 300072；2.天津市自来水集团公司，天津 300040）
|
摘 要：介绍了给水管网可靠度计算的两种方法，重点介绍了最小割集法。首次提出了给水管 网计算中管件可用度和平均可用度的概念，并建立了相应的数学模型。通过管网故障状态下 水力模拟仿真模块，计算节点和系统的最小割集，利用可用度计算模块，计算出组件的可用 度和不可用度，进而求得节点和系统的可靠度。 关键词： 给水管网； 最小割集； 可靠度； 可用度 中图分类号： TU991.33 文献标识码： C 文章编号： 1000-4602(2000)01-0057-04 给水系统的可靠性是指在一定时间内，在一定的运行条件下给水系统完成预定功能的性质。给水系统的功能是保证用水对象获得所需的水量、水压和水质。如果用水过程中发生恶化 , 使这三个指标中任意一个下降到允许的限度以下，就说系统发生了故障(本文研究的是水 压 和水量恶化引起故障的情况)。按照系统发生故障的原因，可以将可靠性分成两大类：机械 可靠性（Mechanic Reliability）和水力可靠性（Hydraulic Reliability）。前者主要是 由于管件发生损坏或被隔断引起的供水能力的下降，通常通过维修或更新来处理；而后者主 要是由于长期使用，管件内部摩阻增大，造成水头损失增大而引起, 通常可采取水泵调节或 刮管涂衬来处理。本文主要研究基于机械可靠性的给水管网可靠度的计算。 1可用度模型及可靠度算法 对所研究的对象有四点假设： ① 部件和系统都只有正常和失效两种状态； 设系统S由n个部件组成，用二值变量Xi来表示第i个部件的状态
系统的状态可以用下述函数ø(X)=ø(x1,x2,…,xn)来 表示，其中X是n维向量øX=(x1,x2,…,xn)，ø(X)是n维向量的二值函数， 即n元二值函数，有
把这种n维二值变量的二值函数称为n维结构函数。 计算给水管网可靠度可以采用以下两种方法。 ① 不交最小路算法 a. 对于简单小型管网，可以采用邻接矩阵法求出所有最小路。 b. 对于大型复杂管网，节点数很大，邻接矩阵往往是稀疏矩阵，因此用邻接矩阵法求 最小路要大容量去存贮，可以采用计算机算法求最小路。 若已知全部最小路集，由初等概率论知
得可靠度：
② 对于多水源管网，应用枚举法求系统及节点的最小割集较为方便。 本文采用最小割集法（Minimum Cut Set）求节点和系统的可靠度。 按图论的理论，所谓割集，就是系统部件的状态变量集合（x1,x2,…,xn）中满足下列条件的子集：设该子集为{xj1,…,xja},j=1,…,m，{ xj1,…,xja}∈{ x1,x2,…,xn}， 当xj1=…=xja=0时 ,Ø(x)=0，即该子集所对应的全体部件故障时，系统s必定故障。 其中m为割集数，为结构函数。那些去掉一个部件即不构成割集的割集称为最小割集 。 城市供水管网实际运作是个动态过程，用管网图形的拓扑结构生成最小割集进而计算节点和系统的可靠度，并不能很好地反映供水管网运行状态的实时变化。因此，我们可以采用管 网故障状态水力模拟的方法确定实际状态下的最小割集。 在实际管网中，当管段发生故障，可以采用关闸方式关闭有关管段进行维修，这时认为管段 被从管网中除去, 就可以用水力计算方法进行管段（管段组合）故障状态下的状态模拟。如 果节点水压高程(水量)值小于所要求的节点水压高程(水量)的最低值, 则认为节点发生故障 ，该管段（管段组合）就是这个节点的一个最小割集，同时也是系统的一个最小割集。对管 网中每根管段及管段的组合进行上述模拟，就可以求出节点和系统的全部最小割集。对于一般城市管网，枚举各种管段组合的数量是相当大的，又因为两个以上管段同时发生 故障的概率很小，所以可以首先只对单个管段进行故障状态下的模拟，然后对优化结果进行单个管段和两管段同时发生故障的可靠度的计算。如果结果不甚令人满意，可重新考虑对单 个管段和两个管段同时故障时的模拟。 在以往的文献中，一般都是以可靠度这个指标来表示给水管网的可靠性。它没有反映出管段 、节点和系统的各自特点以及它们之间的相互关系。 根据定义，部件寿命分布函数为 一般R(t)称为该部件的可靠度函数或可靠度。 由定义可知，可靠度表示了部件寿命这个随机变量，对于初始状态部件为完好的情形，F(t) 是部 件首次出现故障前的时间分布，R(t)为从某个初始状态出发，在时刻t仍处于正常状态的概率，这里t是指寿命时间。而给水管网属于可修复系统，对于可修复系统，系统可由正常和故障 两个状态不断交替的过程来描述。一般来说， 给水系统组件的故障和修复率都服从指数分布。由指数分布的无记忆性可知，初始时刻t以后系统发展的概率规律完全由时刻t是正常还 是故障所决定，而与该系统在时刻t已经工作了多长时间或已修理了多长时间无关。因此，组件在任意时刻正常工作的概率应该用可用度A(t)来表示。可用度的定义为：组件在规定的 条件下，在任意时刻上正常工作的概率。它的数学式表示为：
式中 λ-- 故障强度，也叫失效率，表示单位时间单位长度管线发生故障的次数 ，a-1·km-1 μ-- 修复强度，也叫维修率，表示单位时间管段维修的次数，h-1 这里研究初始时刻部件处于完好状态的情形，同理也可以求得该时刻部件处于故障状态的情 形。考虑［0,T］时间段内组件的平均可用度，即为求A(t)的数学期望E［A(t)］。由函数的数学 期望公式得： 相对于可用度A(t)是不可用度Q(t)，定义为：组件在规定的工作条件下使用时，在任意时刻 上处于故障状态的概率。 Q(t) =1-A(t) (12) 研究节点和系统的可靠性，引入可靠度这一概念。节点可靠度的定义为： 在某一时 间段内，节点的水压高程大于节点所要求的最小水压高程的概率。系统可靠度的定义为： 在某一时间段内，所有节点都可靠的概率。 设某节点j或系统s的一个最小割集MCi中包含n条管段, 则该最小割集的故障概率为：
Rj=1-Pj (17) 本文提出的算法可以用流程图（图1）来表示。
① 搜集原始数据，获得管网结构图，平差基本数据，组件故障率、修复率资料及节点和系统可靠度指标值。 2 算例 本文采用［1］的管网（图2）对以上提出的算法进行了验算，并编制了相应的程序 。
对于一些原始数据，如管径、管长、节点流量、节点初始水压、节点压力水头的约束均 采用 了该篇文章的数据，对于另外一些参数，如故障率、修复率则是通过查阅大量文献资料并结 合我国国情现状确定的（见表1）。建议在实际应用中，这些参数从历史资料的统计分析中 获得。算例中研究的时间区间为2年，由于该管网规模较小，枚举的运算量不是很大，所以 讨论了一条或两条管段故障的情况，在586PⅡ300机上占CPU时间约3min。 计算结果见表1。 从表1可以看出：水源点附近的点比距离水源点较远的点的可靠性高；当降低故障 率或提高修复率时，节点或系统的可靠度都会提高。这一结果对运行的管理和控制有一定的 实际指导意义。 3 结论 综上所述，讨论给水管网中节点及系统可靠度，首先应在掌握历史资料的基础上 ，通过统计回归的方法得到各管段的故障及维修参数λ和μ，求出各管段的可用度A(t)和不可用度 Q(t) ，然后调用水力模拟仿真模块计算出各个节点在故障状态下的水压Hj，并与允许的水压最 低值Hjmin进行比较，如果Hj＜Hjmin，则模拟的故障管段即为最小 割集 。在得到各节点及系统的最小割集后，运用可靠度计算模块计算节点及系统的可靠度。
参考文献： ［1］ 徐祖信.水分配系统以可靠性为基础的线性优化模式［J］.同济大学学报, 1 996，5. ［2］布拉莫夫.给水系统可靠性［M］.北京：中国建筑工业出版社，1990 |
(16)
