[0094] 为使本发明的动态评估方法易于明白,以下结合附图和具体实例,对本发明的实施方式进行阐述。下面以S市核心城区2016-6-5日7:55~8:00的管网DMA区为例,从中选取17个节点在不同指标下进行动态综合评估。
[0095] 1.初级动态评估及权重确定
[0096] 运行状态初级评估需要根据SCADA历史数据和行业标准确定相应的评价标准,部分评价指标需要通过惩罚曲线进行转换,以此求出水力性能和水质性能中评估指标的评价向量。
[0097] (1)评价指标的选取
[0098] 在供水管网运行的实际情况中,可以选择一个甚至多个指标进行动态评估。本发明采用因子分析法对当前时刻前5个五分钟中影响供水管网动态评估的节点压力、节点压力波动、管道流量、余氯、浊度进行分析,确定下一个五分钟对管网动态评估具有代表性的指标。如表3所示。
[0099] 表3在7:30-7:55时间断内各指标的值
[0100] 节点压力x1 节点压力波动x2 管道流量x3 余氯x4 浊度x57:30-7:35 28.3949 0.0316 1331.05 0.486 0.087
7:35-7:40 28.4265 0.0873 1347.14 0.460 0.076
7:40-7:45 28.5108 0.0843 1366.95 0.403 0.069
7:45-7:50 28.6427 0.1319 1391.96 0.464 0.084
7:50-7:55 28.8320 0.1893 1419.76 0.459 0.076
[0101] 分析步骤如下:
[0102] 1)选取前一时间段内的待选择的管网评价指标原始数据进行标准化处理,以消除量纲的影响。
[0103]
[0104] 其中 n为样本数,代表不同时刻或者某一时间段内的时刻数,j=1,2,…,p,p为指标数量。
[0105] 计算得到的标准化矩阵X为:
[0106]
[0107] 2)计算指标变量的相关系数矩阵A,求相关系数矩阵的前P个特征值,如表4所示。
[0108]
[0109] 表4对应成份的特征值和累积贡献率
[0110]成份 特征值(λ) 累积贡献率z/%
1 3.092 61.842
2 1.755 35.107
3 0.130 2.604
4 0.022 0.446
5 2.5e-16 5.028e-15
[0111] λ1=3.092,λ2=1.755,λ3=0.130,λ4=0.022,λ5=2.5e-16[0112] 3)确定公共因子的个数,设为m个,一般按累积贡献率不小于85%来确定。
[0113] 成份1和2累积贡献率大于85%。则选择成份1和2作为公共因子。m=2。
[0114] 4)估计初始因子载荷矩阵A=(aij)p×m,如表5所示。
[0115] 表5初始因子载荷矩阵
[0116]
[0117] 5)对初始因子载荷矩阵进行方差最大正交旋转变换,得到比较理想的新的因子载荷矩阵(如表6)。并对各公共因子进行合理命名和解释。
[0118] 表6变换后因子载荷矩阵
[0119]
[0120] 由因子载荷矩阵看出成份1在节点压力、节点压力波动、管道流量上的载荷值比较大,所以命名成份1为水力因子;命名成份2在余氯和浊度上的载荷因子比较大,所以成份2为水质因子。
[0121] 6)确定不同的公共因子指标系数t以及主成分的线性组合。
[0122] Fm=t1mx1+t2mx2+…+tpmxp
[0123] 其中, Fpm为因子载荷矩阵数据(如表6),x1,x2,…,xp为参与的指标,t1,t2,…,tp为参与指标的权重系数,p为指标的数量,m为公共因子数。
[0124] 计算得到:
[0125] t11=0.564,t21=0.561,t31=0.558,t41=-0.015,t51=-0.082[0126] t12=-0.021,t22=-0.098,t32=-0.075,t42=-0.735,t52=0.726[0127] F1=0.564x1+0.561x2+0.558x3-0.015x4-0.082x5
[0128] F2=-0.021x1-0.098x2-0.075x3+0.735x4+0.726x5
[0129] 7)以各因子的方差贡献率作为权重,对因子的得分加权平均得到综合评价函数以及指标权重系数排名。通过各指标权重系数的排名,确定选择前三个指标进行下一时刻的管网的综合评价。
[0130] 综合得分F=w1x1+w2x2+…+wpxp,
[0131] 其中, z为各主成份的方差,p为参与分析的指标数,m为主成份数。
[0132] 计算得到:F=0.352x1+0.322x2+0.329x3+0.257x4+0.211x5
[0133] 表7指标权重系数及排名
[0134] 指标系数 排名
w1 0.352 1
w2 0.322 3
w3 0.329 2
w4 0.257 4
w5 0.211 5
[0135] 综上所述,根据前五个五分钟的分析,指标权重系数排名前三的为节点压力、管道流量、节点压力波动,即对供水管网的影响相对较大的三个指标,所以在7:55-8:00时选择节点压力、节点压力波动、管道流量进行评价。
[0136] (2)构建供水管网评价指标模糊子集:
[0137] H={h1,h2,…,hn};
[0138] DH={dh1,dh2,…,dhn};
[0139] F={f1,f2,…,fm};
[0140] 其中,H为节点压力的模糊子集,DH代表节点压力波动子集,F代表管道流量子集,n为参与水力评估的节点数,m代表参与水力评估的管道数。
[0141] 当n=m=17时,得到:
[0142] H={h1,h2,…,h17};
[0143] DH={dh1,dh2,…,dh17};
[0144] F={f1,f2,…,f17};
[0145] (3)确定使用的评语集为C={优秀,良好,合格,不合格}。如表8所示。
[0146] 根据惩罚曲线(如图1-3)将节点压力值、节点压力波动值、管道流量值分别转换到0~4之间的数值,转换完成后根据统一的评估标准对节点压力、节点压力波动、管道流量进行模糊评估。
[0147] 表8节点压力评估标准
[0148] 优秀 良好 合格 不合格
压力转换值 3.5-4 2-3.5 1-2 0-1
[0149] 同时,令转换值为x,并且求出转换值在评估标准中四种分级隶属度[0150]
[0151] 可以得到转换值(如表9)。
[0152] 表9各指标的转换值
[0153]
[0154]
[0155] (4)求解不同因素子集的评判矩阵
[0156]
[0157] (5)求不同选取指标的权重系数
[0158] 对不同指标的相关性分析,相关性越高代表两个相关的密切程度,说明与之相关指标因素对此指标更具有代表性。
[0159] 1)节点压力相关权重系数
[0160] 根据节点之间的平均相关性来确定。将每个压力节点与其他压力节点分别进行相关性分析并求解平均相关性。平均相关性越高的压力节点,表明该节点更为准确地反映整个管网情况。
[0161]
[0162] 上式中,ω1i表示节点i的权重, 表示节点i的平均相关性,n为参与评估的节点数。
[0163] 当n=17可得到:
[0164] ωH=(0.058 0.062 0.064 0.063 0.060 0.057 0.057 0.049 0.060 0.059 0.063 0.057 0.055 0.056 0.054 0.063 0.062)
[0165] ωDH=(0.058 0.062 0.064 0.063 0.060 0.058 0.057 0.049 0.060 0.059 0.063 0.057 0.055 0.056 0.054 0.063 0.062)
[0166] 2)管道流量相关权重系数根据管道直径确定。
[0167]
[0168] 上式中,ω2i表示管道j的权重,Dj表示管道j的管径,m为参与评估的管道数。
[0169] m=17时,计算得到:
[0170] ωF=(0.087 0.087 0.070 0.087 0.087 0.034 0.026 0.053 0.053 0.053 0.070 0.061 0.061 0.018 0.054 0.027 0.07)(6)求出评价向量
[0171] 根据不同因素评价矩阵和权重子集,采用 算子,计算评价向量:
[0172]
[0173] 其中,B为评价向量,ω为指标权重子集,R为指标评价矩阵,y为某个需要评价的指标。
[0174] 计算可以得到:
[0175] 节点压力评价向量BH为:
[0176] 节点压力波动评价向量BDH为:
[0177] 管道流量波动评价向量BF为:
[0178] 根据最大隶属度原则,从单一指标进行评估的结果分别是“良好”、“及格”、“良好”。根据不同节点和管道指标相关性的不同,可以选取多个或单个指标进行综合的动态评估。
[0179] 2.二级动态评估及权重确定
[0180] (1)将一级评估得到的节点压力、节点压力波动、管道流量评价向量合并成综合评价矩阵B1为:
[0181]
[0182] (2)使用层次分析法进行分析及权重确定,步骤如下:
[0183] 1)建立层次结构模型
[0184] 根据评价目标及问题建立目标层,准则层,指标层的三层结构模型[0185] 2)建立判断矩阵
[0186] 将一级评估中的指标根据层次分析法,确定这些指标的权重。指标之间的重要性通过“1~9”标度法得到评估指标间重要性表10。
[0187] 表10评估指标间的重要性表
[0188]
[0189] 根据上述表11得到判断矩阵A,
[0190]
[0191] 3)一致性检验及计算权重向量、权重系数
[0192] 定义评估因素V={v1,v2,…,vn}中各个目标的权重为w1,w2,…,wn则首先将判断矩阵按行相加。
[0193]
[0194] 计算得到:M1=6, M3=6
[0195] 经过规范化后得到权重向量wi:
[0196]
[0197] w1=0.38,w2=0.26,w3=0.38
[0198] 设λmax是判断矩阵A的最大特征值,计算得到λmax=3,将最大特征值代入下面公式中
[0199]
[0200] 其中,n为判断矩阵的阶数。
[0201] 得到C.I.=0,因为判断矩阵阶数大于2阶时,难以构造满足一致性的矩阵,因此使用平均随机一致性指标R.I.修正C.I.,即公式
[0202]
[0203] 当C.R<0.1,满足判断矩阵的一致性检验结果。所以将C.I.=0代入上述公式中,得到C.R=0<0.1由此得到,评价因素的权重向量W:
[0204] W=(0.38 0.26 0.38)
[0205] 4)根据 算子得到二级评估向量B’:
[0206]
[0207] (3)确定评估结果
[0208] 根据最大隶属原则,确定2016-6-5日7:55~8:00最终评估结果为“良好”。