基于复杂网络的长距离输水系统节点重要性评估

杨元媛，朱昊彧，杜云彬，李 港

（1.南京水利科学研究院水文水资源与水利工程科学国家重点实验室，南京210029；2.长江科学院流域水环境研究所，武汉430010；3.四川发展国润水务投资有限公司，成都610047）

建设长距离输水系统是解决区域水资源紧缺问题的重要工程措施［1，2］，而闸阀、泵站等控配水建筑物作为长距离输水系统的重要组成部分［3］，对其开展重要性评估，根据其重要程度进行相应的水力设计和防护资源分配，可保障长距离输水系统的安全可靠运行，有效发挥其输水效益。目前的相关研究主要是针对流经控配水建筑物的水流特性对其破坏程度的大小来确定其重要性，梁圣辰等［4］模拟了粤东三江引水工程泵站的水锤效应，提出了针对重要泵站的联合防护方案；张焱炜［5］以贵州某管道输水工程为背景，进行了输水系统中泵站水击防护方案的设计和优化。然而，由输水管线和控配水建筑物构成的长距离输水系统同时还具有复杂的网络拓扑结构［6］，目前从复杂网络的角度进行控配水建筑物重要性评价的研究还较为匮乏。因此，本文以输水管线为边，控配水建筑物为节点，构建长距离输水系统拓扑网络，从节点的局部属性、全局属性和位置属性三个方面筛选指标，进行节点重要性的综合评估，可避免复杂的水力计算，并利用节点移除后的输水网络综合效率验证评估结果，为长距离输水系统的重要节点识别和监测保护提供理论指导。

1 节点重要性评估方法

1.1 输水系统拓扑网络构建

长距离输水系统的输水线路可能横跨几个甚至十几个省市，输水线路长度可达数百公里，而输水系统往往受到功能需求以及环境条件的制约，需要设置数量众多的输水建筑物（渠道、管道等）和控配水建筑物（泵站、闸阀等），形成连通发达的网络系统［7］。因此，运用复杂网络理论，将长距离输水系统输水设施抽象为边，控配水设施抽象为节点，以水流流向为方向，构建长距离输水网络拓扑结构图，记为G：

1.2 节点重要度指标确定

由于输水节点所处的位置以及与其连接的输水设施数量不同，使得节点的重要程度存在一定的差异［8］，可从输水网络中节点与邻居节点的联系（局部属性）、与全网节点间的联系（全局属性）以及节点所处的位置（位置属性）3 个方面来进行节点重要度指标的选取。

1.2.1 局部属性

输水节点的邻居节点数目可表征该节点与邻居节点建立联系的能力，邻居节点越多，输水节点失效后对输水系统的影响越大，引入节点的度［9］来表示节点的邻居节点数目：

式中：Ki为节点的度，表征与该节点连接的上下级节点数量；为节点的入度；为节点的出度，节点的度越大，与之相连节点越多，节点失效对整个输水网络影响越大。

对于度相同的输水节点，邻居节点间的相互联系程度不同，节点失效后对整个输水网络影响也存在一定的差异，邻居节点联系越紧密，节点失效后水流传输的可替代路径越多，对输水系统的影响越低。因此，引入输水节点的聚类系数［10］来表示邻居节点间的联系紧密程度，表示为邻居节点间实际相连边数与可能相连边数的比值。

式中：ci为输水节点的聚类系数；ekj为节点vi与其邻居节点的实际相连边数；为节点vi邻居节点的可能相连边数，聚类系数取值范围介于0～1，当聚类系数接近1 时，表明节点的邻居节点之间联系越紧密，节点重要性越低。

节点的度考虑了节点与相邻节点的联系程度，而聚类系数考虑了节点的邻居节点间联系紧密程度，本研究综合考虑上述两方面对节点局部重要性的影响，定义节点的局部中心性为：

式中：LCi为输水节点的局部中心性；α为重要度影响分配系数，取值为0.5。

1.2.2 全局属性

输水节点与全网中除邻居节点外其他位置节点的平均最短距离越小，则表明该节点与网络中其他节点越接近，水流传输到其他节点越容易，引入节点接近中心性［11］来表征节点与全网中除邻居节点外其他输水节点的联系程度：

式中：CCi为输水节点的接近中心性；dij为节点vi与节点vj之间的最短路径，若节点vi与节点vj之间没有路径可达，则dij= ∞。CC(i)值越大，表明节点vi与网络中其他节点距离越近，联系越紧密，失效后对输水系统影响越大。

1.2.3 位置属性

在输水网络中，一对节点之间通常存在多条最短路径，经过某个节点的最短路径数量越多，则该节点越重要。引入节点的介数中心性［12］为网络中节点对之间经过某节点的最短路径条数与节点对所有最短路径条数之比：

式中：gst为节点vs到节点vt的最短路径数目；为节点vs到节点vt的gst条路径中经过节点vi的最短路径数目，BC越大，节点vi失效后对输水系统影响越大。

处于网络中不同层级的节点其重要程度存在一定的差异，位于网络边缘的节点重要性相较于处于网络中心的节点其重要程度相对较低，利用K-shell分解方法［13］递归的移除网络中所有度值小于或等于kmin的节点来逐步确定节点在网络中的层级及其重要性，考虑到节点度指标中出度入度差异，定义度指标影响分配系数λ，提出节点交叉度概念：

式中：ki为输水节点的交叉度；λ取值为0.5，同时考虑出度入度为0 时计算的交叉度与实际存在偏差情况，引入一个虚拟的全局节点［14］，使该节点与网络中全部节点双向连接，则网络中所有节点的出度和入度都增加1，最终得到改进后的节点交叉度：

设初始网络为G0，根据上式，计算网络中所有节点的kc值，定义kc值最小的节点Ks值为Ks1= 1，移除网络中kc值最小的节点（集）及与其相连的边得到子图G1，更新G1中节点的kin和kout，计算剩余节点kc值，kc值最小的节点Ks值为。重复上述步骤直至网络中所有节点都有其对应的Ks值。

1.3 节点重要度计算

为避免单一指标评价片面性，综合考虑四种指标，通过熵权法［15］计算各指标权重。熵权法根据各项指标提供的信息熵大小确定权重，可以避免确定权重时的人为主观影响。具体步骤如下：

式中：Rij为计算指标x′ij的比重；ej为第j项指标的熵值；gj为第j项指标的差异性系数，值越大表明相应指标在评价中的重要性越强；wj为各指标的权重。

指标权重越大，表示该指标对节点重要度评价时的作用越大，节点重要度表示为各节点重要度指标与其权重之积：

式中：I(i)为节点的综合重要度，I(i)越大，节点vi在输水网络中越重要，失效后对输水系统影响越大。

1.4 输水网络综合效率计算

长距离输水网络中节点失效后，输水网络的结构发生变化，使得水流在输水网络中的流动重新分布，本文引入复杂网络抗毁性测度指标网络综合效率［16］，用来表示节点失效后网络输水功能的破坏程度以及水流在输水网络中流通的难易程度［17］：

式中：E为网络综合效率值；n为初始网络输水节点数目。

当输水网络中有输水节点vk遭到破坏时，若该节点vk不在节点对vi和vj的路径上，水流从节点vi输送到vj的最短路径不发生改变；若节点对vi和vj之间存在多条输水路径，节点vk处在节点对vi和vj的非最短路径上，则节点vk失效时，水流仍可以通过节点vi和节点vj的最短路径进行输送，输水网络综合效率保持不变；当节点vk处于节点对vi和vj之间的路径中的最短路径时，节点vk一旦遭到破坏，节点vi和vj之间的最短路径必然会被其他边取代，使得水流从节点vi传输到节点vj的距离变长，进一步影响到节点vi和vj之间的输水效率，使得水网综合效率变低。

2 案例分析

2.1 工程概况

图1表示某长距离输水工程，该工程包括泵站、闸阀等控配水设施以及管道、渠道等输水设施，将控配水设施抽象为节点，输水设施抽象为边，得到输水网络拓扑结构图，共包括67 个输水节点和116条输水边，1、4、9、22为水源节点。

图1 输水系统网络拓扑结构Fig.1 The topology of water delivery system

2.2 输水节点重要度分析

根据式（3）～（9）计算得到4个节点重要度指标，由于不同的指标计算结果处于不同的量级，经归一化后各指标下节点的重要度排序如图2所示。从图2中可以看出LC和Ks的排序结果呈现为一定的阶梯型分布，同一阶梯上的节点重要性相同，使得部分节点的重要性无法有效区分；BC和CC的排序结果虽然较为平滑，但根据单一指标的计算结果进行节点的重要性评价存在较大的片面性，需要从不同的角度选取多个节点重要性指标进行综合评价。

图2 节点重要度指标评价结果Fig.2 Evaluation results of node importance indices

利用熵权法计算得到LC、BC、CC、Ks的权重分别为0.247、0.238、0.252、0.263，根据式（14）计算得到各节点重要度I(i)，结果如表1所示，重要度较高的节点所处位置如图3所示。结合表1 和图3 可以看出，重要度排序前10 的节点为39、27、43、32、44、45、30、10、19、25，主要位于输水网络的中部位置，在实际的长距离输水工程中应对这些重要输水节点进行重点监测和保护。

图3 重要节点分布图Fig.3 Distribution of important nodes

表1 节点重要度Tab.1 Node importance

2.3 输水网络综合效率分析

根据式（15）计算得到输水系统初始网络综合效率值为0.12，再依次按节点重要度顺序移除节点并计算每个节点移除后的输水网络综合效率值［图4（a）］，同时比较了重要度排序前10 的节点单独移除后输水网络综合效率值的降低程度［图4（b）］。移除23个输水节点后，输水网络综合效率值从0.12下降到0.02，34%的节点影响了输水网络83%的综合效率值，即少量重要节点失效可导致输水系统的迅速崩溃，同时，按节点重要度顺序单独移除节点时，输水网络综合效率值的降低程度总体上随着节点重要度的降低而降低，表明节点的重要度评价结果总体较为准确，能够有效识别长距离输水系统中的重要节点。

图4 依次移除节点时输水网络综合效率演变过程及单独移除重要节点时水综合效率降低程度Fig.4 The evolution process of water conveyance comprehensive efficiency when nodes are removed in sequence and degradation degree when important nodes are removed separately

3 结论

（1）将长距离输水系统控配水建筑物抽象为节点，基于复杂网络理论从局部属性、全局属性和位置属性评估了节点的重要性，提出的综合评估方法可有效识别关键节点，为其针对性监测及防护资源分配提供理论指导。

（2）长距离输水网络中少部分关键节点失效即可对长距离输水系统造成较大的破坏，本案例中34%的节点影响了输水网络83%的输水网络综合效率，同时，位于输水网络中心的节点重要性往往大于边缘节点，应对这些节点进行重点监测和防护。

本研究将输水管线抽象为复杂网络的边时未考虑到输水管线的性质差异（如管道流量大小、有压无压等）对节点的重要性影响，如何区分上述差异并将其引入到输水系统节点重要性评估中来是笔者的下一步研究方向。□