城市给水管网可靠度的计算

《城市给水管网可靠度的计算》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、城市给水管网可靠度的计算赵新华1，刘英梅1，乔宇2(1.天津大学环境科学与工程学院，天津300072；2.天津大学机械学院，天津300072)摘要：研究了城市给水管网可靠度的计算，首先从可靠度与图论原理出发，阐释了算法设计的依据，然后从管网结构的角度，以最小路原理为依据，应用最小路不交化方法设计了给水管网可靠度的算法并进行了实例验证。关键词：城市给水管网；可靠度；最小路；不交化中图分类号：TU991.33文献标识码：C文章编号：1000-4602(2002)04-0053-03以最小路不交化算法为核心，从管网的网络系统结构入手，分析管段的关联关系对管网可靠度的影响，采用遍历计算的方式，

2、从单个管段的可靠度计算出管网中水源点与任意非水源点的可靠度，取最小值作为管网系统可靠性的评价。1计算依据如果给定一个网络，将其整体作为一个系统、每条弧作为一个组件。系统可靠性理论的中心问题是要确定系统可靠度与其组件可靠度之间的关系。1.1管网组件可用度A给水管网在实际运行中经常出现管段故障的情况，经维修后仍投入正常运行，那么称这样的管段为可维修组件，其可靠度用组件可用度A表示。组件在规定的条件下，在任意时刻上正常工作的概率称为组件瞬时可用度A(t)；与瞬时可用度相对应的是瞬时不可用度，表示组件在规定的条件下使用时任意时刻故障的概率用Q(t)表示，显然Q(t)=1-A(t)。采用组件瞬时

3、可用度公式：式中λ——管网的故障率μ——管网的修复率将组件可用度A定义为管段从开始到使用寿命T这一时间段内的平均可用度，即求A(t)的数学期望E［A(t)］。由函数的数学期望公式得：又因为f(t)在［0，T］时间段内服从均匀分布，f(t)=1/T、t∈［0，T］，故组件可用度的数学式表示为：相应的组件不可用度Q=1-A。1.2管网系统可靠度在研究系统可靠性时，首先对所研究的对象作4点假设：①组件和系统都只有正常和失效两种状态；②系统的状态完全由系统的逻辑关系和组件的状态决定；③组件的状态转移率即故障率λ和修复率μ均为常数；④组件或系统的故障和修复都相互独立。采用最小路不交化算法求解。设

4、G是一个给定的网络，V1、V2是指定的两个节点，则基本问题是求时刻T由输入节点V1可以到达输出节点V2的概率，亦即求时刻T系统正常的概率R。对网络的路及最小路的定义如下：定义1对于给定网络G，从指定的节点V1经过一串弧序列(或其中的一部分弧)可以到达节点V2，则称这个弧序列为从V1到V2的一条路。定义2对于给定网络G，从节点V1到V2的弧序列称作一条最小路，若满足①它是一条路，②最小性，即从这个弧序列中除去任意一条弧后它即不是从V1到V2的路。由于系统正常这一事件可表示为S={V1可以到达V2}={由V1至少有一条最小路通到V2}，这里最小路是指组成这条最小路的弧都正常这一事件。因此，

5、若记A1，A2，…，Am为网络的所有最小路，则，，即将求可靠度问题转化为求随机事件和的问题。由概率论可知：称为容斥定理。若能把系统正常这一事件表示成不交(即相斥)事件之和称为相斥求和公式。当最小路数目不多时，用直观方法可以得到全部最小路；若系统结构复杂，产生的最小路数目较多时，应采用图论中的邻接矩阵法［1、2］求全部最小路。用最小路不交化方法将相依事件转化为相斥事件，然后利用相斥求和公式求解系统的可靠度，从而简化了计算过程。定理1设A1，KAm为网络G的所有最小路，A1，L，Ar的长度都小于n-1，Ar+1，L，Am的长度为n-1，且r＜m，则：其中右端、Br+1、L、Bm之间不交，而

6、Br+1、L、Bm为Ar+1、L、Am中分别添上在其中不出现的(l-n+1)条弧的逆而成之事件。2算法设计对于已知结构的管网系统，每个水源点与每个非水源点之间的系统可靠度都可以采用最小路不交化算法求得，即：将水源点处理为输入节点，非水源点处理为输出节点，采用邻接矩阵法求得网络的所有最小路，并对所求得的所有最小路进行不交化处理，就可以结合组件的可用度进行此时的系统可靠度计算。假设管网的管段可用度Ax与管段不可用度Qx均已求出，且通过邻接矩阵分析已知A1，KAm为管网G中水源点i与非水源点j之间的所有最小路，A1，L，Ar的长度都小于n-1，Ar+1、L、Am的长度为n-1，则其系统可靠度

7、R的计算过程如下：①挑出长度为n-1的最小路全体Ar+1、KAm，利用公式将求集合并的问题简化为求代数和的形式，Bi是指在Ai中添上在其中不出现的(l-n+1)条弧的逆而成之事件，Bi之间不交，其概率之和作为系统可靠度的一部分，即最小路长度为n-1子集的可靠度R1。循环计算，其中每一步使Aj与A1、L、Aj-1、不交，j=2,……，r(r为集合元素数)。d.F中还有项(即rk+1≥1)吗?若有，则K←K+1，转b步骤)，否则结果R2就作为系统可