CN112464421B - 基于联合信息熵的供水管网漏损识别传感器优化布置方法 - Google Patents

基于联合信息熵的供水管网漏损识别传感器优化布置方法 Download PDF

Info

Publication number
CN112464421B
CN112464421B CN202011318087.9A CN202011318087A CN112464421B CN 112464421 B CN112464421 B CN 112464421B CN 202011318087 A CN202011318087 A CN 202011318087A CN 112464421 B CN112464421 B CN 112464421B
Authority
CN
China
Prior art keywords
sensor
leakage
node
entropy
candidate
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011318087.9A
Other languages
English (en)
Other versions
CN112464421A (zh
Inventor
胡祖康
陈文龙
沈定涛
苏正猛
陈蓓青
谭德宝
顾世祥
张煜
叶松
夏煜
陈喆
王珺珂
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changjiang River Scientific Research Institute Changjiang Water Resources Commission
Original Assignee
Changjiang River Scientific Research Institute Changjiang Water Resources Commission
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Changjiang River Scientific Research Institute Changjiang Water Resources Commission filed Critical Changjiang River Scientific Research Institute Changjiang Water Resources Commission
Priority to CN202011318087.9A priority Critical patent/CN112464421B/zh
Publication of CN112464421A publication Critical patent/CN112464421A/zh
Application granted granted Critical
Publication of CN112464421B publication Critical patent/CN112464421B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/18Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Geometry (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Computing Systems (AREA)
  • Pipeline Systems (AREA)
  • Examining Or Testing Airtightness (AREA)

Abstract

本发明提供一种基于联合信息熵的供水管网漏损识别传感器优化布置方法,其基于联合信息熵对供水管网传感器进行优化布置,得到确定性条件下的传感器最优布局;考虑到传感器发生突然中断的情况,对确定性条件下的传感器布局方案进行改进,得到具有鲁棒性的传感器优化布局方案。该方法不仅能够确保所得到布局方案在确定性条件下具有良好的漏损识别率,同时,在传感器发生故障时所获得的传感器布局方案仍然具有良好的鲁棒性,对提高供水管网漏损识别尤其是在发生故障时的鲁棒性具有重要意义。

Description

基于联合信息熵的供水管网漏损识别传感器优化布置方法
技术领域
本发明属于城市水务信息化领域,具体涉及一种基于联合信息熵的供水管网漏损识别传感器优化布置方法。
背景技术
供水管网是现代城市必需的基础设施,对保障人民生活安定和城市建设发展起着至关重要的作用,被形象地称作城市的“水动脉”。然而,随着城市供水管网整体规模的不断扩大以及在役时间的逐年增加,对供水管网进行监测成为人们关注的焦点。
“漏损”是供水管网维护中常见的问题之一。管网漏损是由设备老化、腐蚀等原因引起的,会造成巨大的经济、社会和环境成本。管网一方面会造成大量的水流失,造成资源浪费;另一方面也会使污染物进入供水管网,影响供水安全。对漏损的快速检测与识别已成为当前城市管理中一个关键问题。
对管网的运行状况进行实时监测是管网漏损检测与定位的重要基础。由于传感器的安装与维护成本,不可能在管网每个节点布设传感器,需要对传感器进行优化布置。传感器的优化布置涉及在传感器数量有限的情况下实现传感器的最优布置,使传感器网络的信息增益最大。
传感器优化布置对漏损检测和识别至关重要,基于信息论的传感器布置方法得到了广泛应用。信息熵度量允许对每个配置中涉及不同数量的传感器进行比较,信息熵能够对传感器的数量和传感器获取的信息做出权衡,从而获得经济效益高的传感器优化布置方案。基于信息熵的方法在传感器优化布置中取得了良好效果,但是之前的方法没有考虑系统不确定性和传感器之间互信息的影响。单独计算每个位置最大的熵而没有考虑各个位置之间的相似性,传感器优化布置应可能减少传感器之间的信息冗余。一些方法利用条件熵考虑传感器中信息冗余的影响,或者对模型参数后验分布的信息熵进行修正以及减少候选模型的数量,但是这些方法计算效率低下,计算成本高。
此外,以往的方法大多将传感器优化布置看作是确定性问题,传感器完全可靠。实际上,由于各种因素的影响传感器可能会发生突然中断,由于传感器数量有限,单个传感器发生故障时可能会造成严重后果。因此,传感器优化布置时需要考虑鲁棒性,确保在传感器发生故障时仍然具有较好地漏损识别性能。
鉴于此,本发明针对管网漏损识别,提出了一种基于联合信息熵的传感器优化布置方法,并考虑了传感器突然发生中断的情况,确保在传感器发生故障时整个监控系统仍然具有较高的漏损识别率。
发明内容
本发明旨在提出一种基于联合信息熵的供水管网传感器优化布置方法,对管网漏损进行识别,并考虑传感器突然中断的情况,得到具有鲁棒性的传感器优化布置方案。
为实现以上目的,本发明提供一种基于联合信息熵的供水管网传感器优化布置方法,包括以下步骤:
步骤(1):利用管网水力模型对各种漏损情形进行模拟,得到各种漏损情形下各候选节点压力模拟测量值;
步骤(2):利用管网水力模型对正常工况进行模拟,得到各候选节点的压力预测值,并与压力模拟测量值进行比较,得到各候选节点各漏损情形下的压力残差矩阵;
步骤(3):基于压力残差矩阵对各种漏损情形进行区间划分,使其分布在不同的区间内,并计算各个候选节点的信息熵,选择信息熵最大的候选节点添加到传感器列表中,在该候选节点划分区间基础上,继续对区间进行细分,并计算各候选节点的信息熵,选择熵最大的候选节点,重复对区间进行细分和信息熵计算,依次对候选节点进行选择,直到添加新的候选节点后信息熵不再继续增大,得到确定性条件下的传感器最优布置方案,所述传感器最优布置方案包括传感器数量N、位置和次序;
步骤(4):使用穷举法对N个传感器依次进行“冷藏”,并重新进行区间划分和信息熵计算,依次对传感器进行选择,得到新的传感器列表,并将“冷藏”的候选节点“解冻”后将其添加到新的传感器列表中,从而得到N种改进方案,并对N种改进方案进行比较,选择传感器发生故障时漏损识别率最高的方案作为具有鲁棒性的传感器最优布置方案。
进一步的,所述步骤(1)具体包括:
步骤(1.1)、利用EPANET建立供水管网水力模型,并根据实测的节点压力和管道流量对水力模型的初始参数进行校正,使管网水力模型模拟供水系统的实际运行状况;
步骤(1.2)、定义各种漏损情形Si
Si=[Θi1,…,Θik];
Θik表示节点i处发生漏损,漏损强度为k;
步骤(1.3)、对各种漏损情形进行模拟,得到各种漏损情形下各候选节点压力模拟测量值g;
Figure GDA0003644220790000031
gij表示漏损情形为Θjk时候选节点i处的压力模拟测量值,即假设节点j发生漏损,漏损强度为k。
进一步的,所述步骤(2)具体包括:
步骤(2.1)、利用校核好的管网水力模型对管网正常运行工况进行模拟,得到各候选节点压力模拟预测值y;
y=[y1,…yi,…yn]
yi表示在正常工况下节点i处的压力模拟预测值;
步骤(2.2)、将正常工况下的压力模拟预测值与各种漏损情形下的压力模拟测量值进行比较,得到压力残差R=y-g;
Figure GDA0003644220790000032
其中rij为节点j处发生漏损,节点i处的水力模型预测值与测量值之差,R中的每一行对应于一个候选节点在各种漏损情形下的压力残差,R中的每一列对应于一种漏损情形在各候选节点的压力残差。
进一步的,所述步骤(3)具体如下:
(3.1)对于各候选节点i,根据区间界限(Iw,i)n将rij划分到不同的区间内,满足
Figure GDA0003644220790000033
Figure GDA0003644220790000034
其中NI是第i个候选节点划分得到的最大区间数;
(3.2)计算各个候选节点位置的熵H(y)i
Figure GDA0003644220790000041
其中,H(R)i是压力残差R在候选节点i处的熵,P(Rn)i是变量分布在第i个传感器第n个区间内的概率,NI为第i个传感器总的区间数;
(3.3)选择熵最大的节点放入传感器节点列表中,并保存当前的区间划分结果;
(3.4)在当前区间划分结果基础上,利用各候选节点对区间进一步划分,并根据步骤(3.2)中的方法计算各候选节点处的熵;
(3.5)对各候选节点的熵进行比较,选择熵最大的候选节点,并比较该候选节点的熵与传感器列表中最后一个节点的熵的大小,如果该候选节点的熵大于传感器列表最后一个节点的熵,则进入步骤(3.3),否则,进入步骤(3.6);
(3.6)将当前的传感器列表作为传感器优化布置方案。
进一步的,所述步骤(4)具体包括:
(4.1)根据所述步骤(3)得到确定性条件下传感器优化布局方案;
(4.2)得到N种具有鲁棒性的传感器布局方案,对于第i种方案,在分层算法的第i步将信息熵最高的传感器“冷藏”,选择次优的候选节点添加到传感器列表中;
(4.3)继续利用分层算法对依次对后续传感器进行选择,得到传感器列表;
(4.4)将第i步“冷藏”的候选节点“解冻”并将其添加到传感器列表中,并将其作为具有鲁棒性的传感器优化布局方案;
(4.5)对各种方案在传感器发生故障时的漏损识别率进行比较,得到在传感器发生故障时漏损识别率最高的方案。
进一步的,所述步骤(4.5)具体如下:
(4.5.1)对于第i种传感器优化布置方案,传感器数量为m;
(4.5.2)对于第i种方案,重复进行m次漏损识别实验,在第j次将第j个传感器“冷藏”,利用剩余的m-1个传感器对各种漏损情形进行识别,得到各种漏损情形下的识别率,并计算m次实验的平均漏损识别率;
(4.5.3)比较N种方案在发生故障时的漏损识别率,选择漏损识别率最高的方案作为具有鲁棒性的传感器优化布置方案。
本发明基于联合信息熵对确定性条件下供水管网漏损识别传感器进行优化布置,在传感器数量有限的情况下,能够得到具有良好信息增益的传感器优化布置方案;同时,基于分层策略的传感器依次选择方法考虑了传感器之间互信息,避免了信息冗余;此外,考虑到实际中传感器可能发生故障的情况,对确定性条件下的传感器优化布置方案做出改进,获得了具有鲁棒性的传感器优化布局方案,确保在传感器发生故障时仍然能够保持良好的漏损检测性能,对提高供水管网漏损识别尤其是在发生故障时的鲁棒性具有重要意义。
附图说明
图1为本发明基于联合信息熵的供水管网传感器优化布置方法的原理图;图1主要分为4个部分:(a)利用管网水力模型对各种漏损情形进行模拟,得到各种漏损情形下各候选节点压力模拟测量值;(b)利用管网水力模型对正常工况进行模拟,得到各候选节点的压力模拟预测值,并与压力模拟测量值进行比较,得到各候选节点各漏损情形下的压力残差矩阵;(c)采用联合信息熵对传感器进行优化布置,获得确定性条件下的传感器最优布置;(d)采用联合信息熵对传感器进行优化布置,获得具有鲁棒性的传感器布置方式;
图2为本发明实施例EPANET示例管网Net3示意图;
图3为本发明实施例供水管网漏损识别流程图;
图4为本发明实施例确定性条件下传感器优化布置示意图;
图5为本发明实施例基于鲁棒性的传感器优化布置示意图;
图6为比较各方案在传感器发生故障时漏损识别性能流程图;
图7为本发明实施例各传感器依次添加传感器网络信息熵以及漏损识别率变化曲线图;
图8为本发明实施例确定性条件下各传感器发生故障时信息熵及漏损识别率变化曲线图;
图9为本发明实施例具有鲁棒性传感器布局各方案在传感器发生故障时漏损识别性能的比较。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例提出一种基于联合信息熵的供水管网漏损识别的传感器优化布置方法,用于供水管网漏损识别的传感器优化布置,在传感器数量有限的情况下,采用联合信息熵选取各候选节点对传感器进行布置,获得确定性条件下的传感器最优布置,并考虑传感器发生故障的情况,获得具有鲁棒性的传感器配置。所述方法具体包括如下步骤:
步骤1、利用管网水力模型对各种漏损情形进行模拟,得到各种漏损情形下各候选节点的压力模拟测量值。所述步骤1具体包括:
步骤1.1、利用EPANET建立供水管网水力模型,并根据实测的节点压力和管道流量对水力模型的初始参数(节点需求和阻力系数)进行校正,使管网水力模型模拟供水系统的实际运行状况;
步骤1.2、定义各种漏损情形Si
Si=[Θi1,…,Θik];
Θik表示节点i处发生漏损,漏损强度为k;
步骤1.3、对各种漏损情形进行模拟,得到各种漏损情形下各候选节点压力模拟测量值g;
Figure GDA0003644220790000061
gij表示漏损情形为Θjk时候选节点i处的压力模拟测量值,即假设节点j发生漏损,漏损强度为k。
如图2所示,Net3示例管网共有92个节点,假设每个节点发生漏损,共有92种漏损情形;
对92种漏损情形进行模拟,得到各候选节点的压力模拟测量值g,g为92×92矩阵;
Figure GDA0003644220790000071
gij为节点j发生漏损时候选节点i处的压力模拟测量值;
步骤2、利用水力模型对正常工况进行模拟,得到正常工况下各候选节点的压力模拟预测值,并与各漏损情形下的各候选节点的压力模拟测量值进行比较,得到各种漏损情形下的压力残差矩阵R。所述步骤2具体包括:
步骤2.1、利用校核好的管网水力模型对管网正常运行工况进行模拟,得到各候选节点压力模拟预测值y;
y=[y1,…yi,…yn]
yi表示在正常工况下节点i处的压力模拟预测值;
步骤2.2、如图1(b)所示,将正常工况下的压力模拟预测值与各种漏损情形下的压力模拟测量值进行比较,得到压力残差R=y-g;
Figure GDA0003644220790000072
其中rij为节点j处发生漏损,节点i处的水力模型预测值与测量值之差,R中的每一行对应于一个候选节点在各种漏损情形下的压力残差,R中的每一列对应于一种漏损情形在各候选节点的压力残差。
步骤3、基于压力残差矩阵对各种漏损情形进行区间划分,使其分布在不同的区间内,并计算各个候选节点的信息熵,选择信息熵最大的候选节点添加到传感器列表中,在该候选节点划分区间基础上,继续对区间进行细分,并计算各候选节点的信息熵,选择熵最大的候选节点,重复对区间进行细分和信息熵计算,依次对候选节点进行选择,直到添加新的候选节点后信息熵不再继续增大,得到确定性条件下的传感器最优布置方案,所述传感器最优布置方案包括传感器数量N、位置和次序。所述步骤3具体包括:
步骤3.1、对于各候选节点i,根据区间界限(Iw,i)n将rij划分到不同的区间内,满足
Figure GDA0003644220790000074
Figure GDA0003644220790000075
其中NI是第i个候选节点划分得到的最大区间数;
步骤3.2、计算各个候选节点位置的熵H(y)i
Figure GDA0003644220790000073
其中,H(R)i是压力残差R在候选节点i处的熵,P(Rn)i是变量分布在第i个传感器第n个区间内的概率,NI为第i个传感器总的区间数;
为了计算熵,计算每个区间mi内的子集数,并计算区间的概率(mi/N);
步骤3.3、选择熵最大的节点放入传感器节点列表中,并保存当前的区间划分结果;
步骤3.4、在当前区间划分结果基础上,利用各候选节点对区间进一步划分,并根据步骤3.2中的方法计算各候选节点处的熵;
步骤3.5、对各候选节点的熵进行比较,选择熵最大的候选节点,并比较该候选节点的熵与传感器列表中最后一个节点的熵的大小,如果该候选节点的熵大于传感器列表最后一个节点的熵,则进入步骤(3.3),否则,进入步骤(3.6);
步骤3.6、将当前的传感器列表作为传感器优化布置方案。
如图4所示,采用分层策略对传感器依次进行选择,一开始所有漏损情形在一个区间内,随着添加传感器,各种漏损情形被划分到不同的区间内,各区间内的漏损情形数不断减少;
计算各个候选节点的熵,选择熵最大的节点添加到传感器列表中;
如图6所示,在第1次区间划分后,节点'113'处的熵最大,因此将节点'113'添加到传感器列表中;
添加节点'113'后,区间划分结果为:[['10','15','20','35','40','50','60','601','61','101'],['103','105','107'],['109','111'],['113','115','117','119','120','121','123','125'],['127','129','131','139','141','143','145','147','149','151','153','157','159','161','163'],['164','166','167','169','171'],['173','177','179','181','183'],['184','185','187'],['189','191','193'],['195','197','199'],['201','203','204'],['205','206','207','208','209'],['211','213','215','217','219'],['225','229','231','237','239','241','243','247','249'],['251','253','255','257','259'],['261','263','265'],['267','269','271'],['273','275']],92种漏损情形被划分为18个区间,传感器网络的熵为2.72,漏损识别率为0%;
添加第2个传感器,所选择节点为'145',区间进一步划分为:[['10'],['15','20','35','40','50','60','601','61'],['101'],['103','105','107'],['109','111'],['113','115','117','119','120'],['121','123','125'],['127','129','131'],['139'],['141'],['143'],['145'],['147'],['149'],['151','153','157','159','161','163'],['164','166'],['167','169','171'],['173','177','179','181','183'],['184','185','187'],['189','191','193'],['195','197'],['199'],['201','203','204'],['205','206','207','208','209'],['211','213','215','217','219'],['225','229'],['231','237','239','241','243','247','249'],['251','253','255','257','259'],['261','263','265'],['267','269'],['271'],['273','275']],92种漏损情形被划分为32个区间,传感器网络的熵为3.26,漏损识别率为10.87%;
添加第3个传感器,所选择节点为'275',区间进一步划分为:[['10'],['15','20','35','40','50','60','601','61'],['101'],['103','105','107'],['109'],['111'],['113','115','117','119','120'],['121','123','125'],['127','129','131'],['139'],['141'],['143'],['145'],['147'],['149'],['151','153','157','159'],['161','163'],['164','166'],['167'],['169','171'],['173','177','179'],['181','183'],['184'],['185','187'],['189','191'],['193'],['195','197'],['199'],['201','203'],['204'],['205'],['206','207'],['208','209'],['211'],['213','215'],['217','219'],['225'],['229'],['231','237','239'],['241','243','247','249'],['251','253'],['255','257','259'],['261','263'],['265'],['267','269'],['271'],['273'],['275']],92种漏损情形被划分为48个区间,传感器网络的熵为3.69,漏损识别率为29.35%;
添加第4个传感器,所选择节点为'107',区间划分为:[['10'],['15','20','35','40','50','60','601','61'],['101'],['103'],['105'],['107'],['109'],['111'],['113','115','117'],['119','120'],['121','123','125'],['127','129','131'],['139'],['141'],['143'],['145'],['147'],['149'],['151','153','157'],['159'],['161','163'],['164','166'],['167'],['169'],['171'],['173','177','179'],['181','183'],['184'],['185','187'],['189'],['191'],['193'],['195','197'],['199'],['201','203'],['204'],['205'],['206','207'],['208','209'],['211'],['213'],['215'],['217','219'],['225'],['229'],['231','237','239'],['241','243'],['247','249'],['251','253'],['255','257'],['259'],['261','263'],['265'],['267','269'],['271'],['273'],['275']],92种漏损情形被划分为57个区间,传感器网络的熵为3.88,漏损识别率为36.96%;
添加第5个传感器,所选择节点为'206',区间划分为:[['10'],['15','20','35','40','50','60','601','61'],['101'],['103'],['105'],['107'],['109'],['111'],['113'],['115','117'],['119','120'],['121','123','125'],['127','129','131'],['139'],['141'],['143'],['145'],['147'],['149'],['151','153','157'],['159'],['161','163'],['164','166'],['167'],['169'],['171'],['173'],['177','179'],['181','183'],['184'],['185','187'],['189'],['191'],['193'],['195','197'],['199'],['201'],['203'],['204'],['205'],['206'],['207'],['208'],['209'],['211'],['213'],['215'],['217'],['219'],['225'],['229'],['231'],['237','239'],['241','243'],['247'],['249'],['251','253'],['255'],['257'],['259'],['261','263'],['265'],['267','269'],['271'],['273'],['275']],92种漏损情形被划分为66个区间,传感器网络的熵为4.04,漏损识别率为53.26%;
添加第6个传感器,所选择节点为'204',区间划分为:[['10'],['15','20','35','40','50','60','601','61'],['101'],['103'],['105'],['107'],['109'],['111'],['113'],['115','117'],['119','120'],['121','123'],['125'],['127','129','131'],['139'],['141'],['143'],['145'],['147'],['149'],['151','153','157'],['159'],['161','163'],['164','166'],['167'],['169'],['171'],['173'],['177','179'],['181'],['183'],['184'],['185','187'],['189'],['191'],['193'],['195'],['197'],['199'],['201'],['203'],['204'],['205'],['206'],['207'],['208'],['209'],['211'],['213'],['215'],['217'],['219'],['225'],['229'],['231'],['237','239'],['241','243'],['247'],['249'],['251'],['253'],['255'],['257'],['259'],['261'],['263'],['265'],['267'],['269'],['271'],['273'],['275']],92种漏损情形被划分为72个区间,传感器网络的熵为4.13,漏损识别率为65.22%;
添加第7个传感器,所选择节点为'601',区间划分为:[['10'],['15','20','35','40','50','60'],['601'],['61'],['101'],['103'],['105'],['107'],['109'],['111'],['113'],['115','117'],['119','120'],['121','123'],['125'],['127','129','131'],['139'],['141'],['143'],['145'],['147'],['149'],['151','153'],['157'],['159'],['161','163'],['164','166'],['167'],['169'],['171'],['173'],['177','179'],['181'],['183'],['184'],['185','187'],['189'],['191'],['193'],['195'],['197'],['199'],['201'],['203'],['204'],['205'],['206'],['207'],['208'],['209'],['211'],['213'],['215'],['217'],['219'],['225'],['229'],['231'],['237','239'],['241','243'],['247'],['249'],['251'],['253'],['255'],['257'],['259'],['261'],['263'],['265'],['267'],['269'],['271'],['273'],['275']],92种漏损情形被划分为75个区间,传感器网络的熵为4.22,漏损识别率为68.48%;
添加第8个传感器,所选择节点为'117',区间划分为:[['10'],['15','20','35','40','50','60'],['601'],['61'],['101'],['103'],['105'],['107'],['109'],['111'],['113'],['115'],['117'],['119','120'],['121'],['123'],['125'],['127','129','131'],['139'],['141'],['143'],['145'],['147'],['149'],['151','153'],['157'],['159'],['161','163'],['164','166'],['167'],['169'],['171'],['173'],['177'],['179'],['181'],['183'],['184'],['185'],['187'],['189'],['191'],['193'],['195'],['197'],['199'],['201'],['203'],['204'],['205'],['206'],['207'],['208'],['209'],['211'],['213'],['215'],['217'],['219'],['225'],['229'],['231'],['237','239'],['241','243'],['247'],['249'],['251'],['253'],['255'],['257'],['259'],['261'],['263'],['265'],['267'],['269'],['271'],['273'],['275']],92种漏损情形被划分为79个区间,传感器网络的熵为4.28,漏损识别率为77.17%;
添加第9个传感器,所选择节点为'151',区间划分为:[['10'],['15','20','35','40','50','60'],['601'],['61'],['101'],['103'],['105'],['107'],['109'],['111'],['113'],['115'],['117'],['119'],['120'],['121'],['123'],['125'],['127','129','131'],['139'],['141'],['143'],['145'],['147'],['149'],['151'],['153'],['157'],['159'],['161'],['163'],['164','166'],['167'],['169'],['171'],['173'],['177'],['179'],['181'],['183'],['184'],['185'],['187'],['189'],['191'],['193'],['195'],['197'],['199'],['201'],['203'],['204'],['205'],['206'],['207'],['208'],['209'],['211'],['213'],['215'],['217'],['219'],['225'],['229'],['231'],['237','239'],['241','243'],['247'],['249'],['251'],['253'],['255'],['257'],['259'],['261'],['263'],['265'],['267'],['269'],['271'],['273'],['275']],92种漏损情形被划分为82个区间,传感器网络的熵为4.32,漏损识别率为83.70%;
添加第10个传感器,所选择节点为'171',区间划分为:[['10'],['15','20','35','40','50','60'],['601'],['61'],['101'],['103'],['105'],['107'],['109'],['111'],['113'],['115'],['117'],['119'],['120'],['121'],['123'],['125'],['127','129','131'],['139'],['141'],['143'],['145'],['147'],['149'],['151'],['153'],['157'],['159'],['161'],['163'],['164'],['166'],['167'],['169'],['171'],['173'],['177'],['179'],['181'],['183'],['184'],['185'],['187'],['189'],['191'],['193'],['195'],['197'],['199'],['201'],['203'],['204'],['205'],['206'],['207'],['208'],['209'],['211'],['213'],['215'],['217'],['219'],['225'],['229'],['231'],['237','239'],['241'],['243'],['247'],['249'],['251'],['253'],['255'],['257'],['259'],['261'],['263'],['265'],['267'],['269'],['271'],['273'],['275']],92种漏损情形被划分为84个区间,传感器网络的熵为4.35,漏损识别率为88.04%;
添加第11个传感器,所选择节点为'205',区间划分为:[['10'],['15','20','35','40','50','60'],['601'],['61'],['101'],['103'],['105'],['107'],['109'],['111'],['113'],['115'],['117'],['119'],['120'],['121'],['123'],['125'],['127','129','131'],['139'],['141'],['143'],['145'],['147'],['149'],['151'],['153'],['157'],['159'],['161'],['163'],['164'],['166'],['167'],['169'],['171'],['173'],['177'],['179'],['181'],['183'],['184'],['185'],['187'],['189'],['191'],['193'],['195'],['197'],['199'],['201'],['203'],['204'],['205'],['206'],['207'],['208'],['209'],['211'],['213'],['215'],['217'],['219'],['225'],['229'],['231'],['237'],['239'],['241'],['243'],['247'],['249'],['251'],['253'],['255'],['257'],['259'],['261'],['263'],['265'],['267'],['269'],['271'],['273'],['275']],92种漏损情形被划分为85个区间,传感器网络的熵为4.37,漏损识别率为90.22%;
在添加第11个传感器后,继续添加新的传感器,区间不能进一步被细分,漏损识别率也不再提高,因此,得到确定性条件下的传感器优化布局方案['113','145','275','107','206','204','601','117','151','171','205'];
如图7所示,为传感器网络的信息熵和漏损识别率随着添加的传感器数量的变化曲线,在添加第1个传感器后漏损识别率为0%,在添加第11个传感器后,漏损识别率升高至90.22%;
考虑到传感器发生故障的情况,利用图6中所示的流程,考虑各个传感器发生故障的情况,得到了在一次有1个传感器发生故障时的漏损识别率;
如图8所示,在各个传感器发生故障时,传感器网络的信息熵以及漏损识别率均显著下降,其中节点'145'处的传感器发生故障时信息熵以及漏损识别率最低。
步骤4、如图1(d)所示,考虑到传感器发生故障的情况,对步骤3中提出的方法做出改进,根据图5中所示的方法得到具有鲁棒性的传感器优化布置方案,具体的,使用穷举法对N个传感器依次进行“冷藏”,并重新进行区间划分和信息熵计算,依次对传感器进行选择,得到新的传感器列表,并将“冷藏”的候选节点“解冻”后将其添加到新的传感器列表中,从而得到N种改进方案,并对N种改进方案进行比较,选择传感器发生故障时漏损识别率最高的方案作为具有鲁棒性的传感器最优布置方案。
所述步骤4具体包括:
步骤4.1、根据步骤3得到确定性条件下传感器最优布局方案,得到传感器的数量11,以及相应的传感器节点位置['113','145','275','107','206','204','601','117','151','171','205'];
步骤4.2、由于确定性条件下传感器的数量为11,因此,基于改进的分层策略得到11种传感器优化布局方案,并对各种方案在传感器发生故障时的漏损识别率以及熵值进行比较,选择最优的方案作为具有鲁棒性的传感器优化布局方案。
(1)对于改进方案1,在分层策略的第1步,将最优候选节点'113'“冷藏”,而选择次优的候选节点'109'作为第1次选择的传感器节点;
继续利用步骤3中的分层策略选择候选的传感器,直到传感器网络的熵不再继续增加,得到传感器列表['109','209','143','187','199','105','237','191','15','601','120','111'];
将节点'113'“解冻”并将其添加到传感器列表中,将整个传感器列表['113','109','209','143','187','199','105','237','191','15','601','120','111']作为传感器优化布局方案1;
(2)对于改进方案2,首先利用分层策略选择第1个节点['113'],并在第2步将最优候选节点'145'“冷藏”,而选择次优的候选节点'208'作为传感器2的布设位置;
在添加第2个传感器后,继续利用分层策略依次对传感器进行选择,直到传感器网络的熵不再继续增加,得到传感器列表['113','208','143','107','207','105','111','189','15','601','10','101'];
将节点'145'“解冻”并将其添加到传感器列表中,将整个传感器列表['113','145','208','143','107','207','105','111','189','15','601','10','101']作为传感器优化布局方案2;
(3)对于改进方案3,首先利用分层策略选择前3个传感器布设位置['113','145'],在第3次选择时将最优候选节点'275'“冷藏”,选择次优的候选节点'243';
在添加第3个传感器后,继续利用分层策略依次对传感器进行选择,直到传感器网络的熵不再增加,得到传感器列表['113','145','243','267','191','151','601','101','111','15','120','149'];
将节点'275'“解冻”并将其添加到传感器列表中,将整个传感器列表['113','145','275','243','267','191','151','601','101','111','15','120','149']作为传感器优化布局方案3;
(4)对于改进方案4,首先利用分层策略选择前3个传感器布设位置['113','145','275'],在第4次选择时将最优候选节点'107'“冷藏”,选择次优的候选节点'243';
在添加第4个传感器后,继续利用分层策略依次对传感器进行选择,直到传感器网络的熵不再增加,得到传感器列表['113','145','275','243','117','171','601','105','204','15','111']
将节点'107'“解冻”并将其添加到传感器列表中,将整个传感器列表['113','145','275','107','243','117','171','601','105','204','15','111']作为传感器优化布局方案4;
(5)对于改进方案5,首先利用分层策略选择前4个传感器布设位置['113','145','275','107'],在第5次选择时将最优候选节点'206'“冷藏”,选择次优的候选节点'211';
在添加第5个传感器后,继续利用分层策略依次对传感器进行选择,直到传感器网络的熵不再增加,得到传感器列表['113','145','275','107','211','117','601','171','105','115','120'];
将节点'206'“解冻”并将其添加到传感器列表中,将整个传感器列表['113','145','275','107','206','211','117','601','171','105','115','120']作为传感器优化布局方案5;
(6)对于改进方案6,首先利用分层策略选择前5个传感器布设位置['113','145','275','107','206'],在第6次选择时将最优候选节点'204'“冷藏”,选择次优的候选节点'117';
在添加第6个传感器后,继续利用分层策略依次对传感器进行选择,直到传感器网络的熵不再增加,得到传感器列表['113','145','275','107','206','117','601','171','211','103','120'];
将节点'204'“解冻”并将其添加到传感器列表中,将整个传感器列表['113','145','275','107','206','204','117','601','171','211','103','120']作为传感器优化布局方案6;
(7)对于改进方案7,首先利用分层策略选择前6个传感器布设位置['113','145','275','107','206','204'],在第7次选择时将最优候选节点'601'“冷藏”,选择次优的候选节点'61';
在添加第7个传感器后,继续利用分层策略依次对传感器进行选择,直到传感器网络的熵不再继续增加,得到传感器列表['113','145','275','107','206','204','61','117','151','171','205'];
将节点'601'“解冻”并将其添加到传感器列表中,将整个传感器列表['113','145','275','107','206','204','601','61','117','151','171','205']作为传感器优化布置方案7;
(8)对于改进方案8,首先利用分层策略选择前7个传感器布设位置['113','145','275','107','206','204','601'],并在第8次选择时将最优候选节点'117'“冷藏”,选择次优的候选节点'103';
在添加第8个传感器后,继续利用分层策略依次对传感器进行选择,直到传感器网络的熵不再继续增加,得到传感器列表['113','145','275','107','206','204','601','103','171','157','111','205'];
将节点'117'“解冻”并将其添加到传感器列表中,将整个传感器列表['113','145','275','107','206','204','601','117','103','171','157','111','205']作为传感器优化布置方案8;
(9)对于改进方案9,首先利用分层策略选择前8个传感器布设位置['113','145','275','107','206','204','601','117'],并在第9次选择时将最优候选节点'151'“冷藏”,选择次优的候选节点'171';
在添加第9个传感器后,继续利用分层策略依次对传感器进行选择,直到传感器网络的熵不再继续增加,得到传感器列表['113','145','275','107','206','204','601','117','171','105','205'];
将节点'151'“解冻”并将其添加到传感器列表中,将整个传感器列表['113','145','275','107','206','204','601','117','151','171','105','205']作为传感器优化布置方案9;
(10)对于改进方案10,首先利用分层策略选择前9个传感器布设位置['113','145','275','107','206','204','601','117','151'],并在第10次选择时将最优候选节点'171'“冷藏”,选择次优的候选节点'173';
在添加第10个传感器后,继续利用分层策略依次对传感器进行选择,直到传感器网络的熵不再继续增加,得到传感器列表['113','145','275','107','206','204','601','117','151','173','205'];
将节点'171'“解冻”并将其添加到传感器列表中,将整个传感器列表['113','145','275','107','206','204','601','117','151','171','173','205']作为传感器优化布置方案10;
(11)对于改进方案11,首先利用分层策略选择前10个传感器布设位置['113','145','275','107','206','204','601','117','151','171'],并在第11次选择时将最优候选节点'205'“冷藏”,选择次优候选节点'208';
在添加第11个传感器后,继续利用分层策略依次对传感器进行选择,直到传感器网络的熵不再继续增加,得到传感器列表['113','145','275','107','206','204','601','117','151','171','208'];
将节点'205'“解冻”并将其添加到传感器列表中,将整个传感器列表['113','145','275','107','206','204','601','117','151','171','205','208']作为传感器优化布置方案11;
在确定性条件下所获得的每个候选节点“冷藏”1次之后,得到了11种改进方案以及确定性条件下的传感器优化布局方案0,如表1所示;
表1,各种传感器优化布局方案汇总;
Figure GDA0003644220790000161
步骤4.3、在得到各种传感器优化布局方案后,利用各方案节点的压力模拟测量值对92种漏损情形进行识别,主要考虑在传感器发生故障时(1次有一个传感器发生故障)各方案的漏损识别率;
(1)对11种改进方案以及确定性条件下传感器优化布局方案(方案12)在传感器发生故障时的漏损识别性能进行比较;
(2)根据图6所示的比较各方案在传感器发生故障时漏损识别性能流程图,对12方案在传感器发生故障时的漏损识别性能进行比较(漏损识别方法如图3所示),得到各方案各个传感器发生故障时的平均漏损识别率;
表2为各方案各个传感器发生故障时的漏损识别率以及平均漏损识别率;
Figure GDA0003644220790000171
(3)选择在传感器发生故障时漏损识别率最高的方案作为具有鲁棒性的传感器优化布局方案;
如图9所示,为各方案在传感器发生故障时的漏损识别率和信息熵,与方案12(确定性条件下传感器优化布局方案),各种改进方案在传感器发生故障时稳健性更好,具有更高的信息熵和漏损识别率,其中方案4(“冷藏”节点'107')时的信息熵和漏损识别率最高,因此,选择方案4作为具有鲁棒性的传感器优化布局方案。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何属于本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。

Claims (5)

1.一种基于联合信息熵的供水管网漏损识别的传感器优化布置方法,其特征在于包括如下步骤:
步骤(1):利用管网水力模型对各种漏损情形进行模拟,得到各种漏损情形下各候选节点压力模拟测量值;
步骤(2):利用管网水力模型对正常工况进行模拟,得到各候选节点的压力预测值,并与压力模拟测量值进行比较,得到各候选节点各漏损情形下的压力残差矩阵;
步骤(3):基于压力残差矩阵对各种漏损情形进行区间划分,使其分布在不同的区间内,并计算各个候选节点的信息熵,选择信息熵最大的候选节点添加到传感器列表中,在该候选节点划分区间基础上,继续对区间进行细分,并计算各候选节点的信息熵,选择熵最大的候选节点,重复对区间进行细分和信息熵计算,依次对候选节点进行选择,直到添加新的候选节点后信息熵不再继续增大,得到确定性条件下的传感器最优布置方案,所述传感器最优布置方案包括传感器数量N、位置和次序;
步骤(4):使用穷举法对N个传感器依次进行“冷藏”,并重新进行区间划分和信息熵计算,依次对传感器进行选择,得到新的传感器列表,并将“冷藏”的候选节点“解冻”后将其添加到新的传感器列表中,从而得到N种改进方案,并对N种改进方案进行比较,选择传感器发生故障时漏损识别率最高的方案作为具有鲁棒性的传感器最优布置方案;
所述步骤(4)具体包括:
(4.1)根据所述步骤(3)得到确定性条件下传感器优化布局方案;
(4.2)得到N种具有鲁棒性的传感器布局方案,对于第i种方案,在分层算法的第i步将信息熵最高的传感器“冷藏”,选择次优的候选节点添加到传感器列表中;
(4.3)继续利用分层算法依次对后续传感器进行选择,得到传感器列表;
(4.4)将第i步“冷藏”的候选节点“解冻”并将其添加到传感器列表中,并将其作为具有鲁棒性的传感器优化布局方案;
(4.5)对各种方案在传感器发生故障时的漏损识别率进行比较,得到在传感器发生故障时漏损识别率最高的方案。
2.如权利要求1所述的基于联合信息熵的供水管网漏损识别的传感器优化布置方法,其特征在于:所述步骤(1)具体包括:
步骤(1.1)、利用EPANET建立供水管网水力模型,并根据实测的节点压力和管道流量对水力模型的初始参数进行校正,使管网水力模型模拟供水系统的实际运行状况;
步骤(1.2)、定义各种漏损情形Si
Si=[Θi1,…,Θik];
Θik表示节点i处发生漏损,漏损强度为k;
步骤(1.3)、对各种漏损情形进行模拟,得到各种漏损情形下各候选节点压力模拟测量值g;
Figure FDA0003653619470000021
gij表示漏损情形为Θjk时候选节点i处的压力模拟测量值,即假设节点j发生漏损,漏损强度为k。
3.如权利要求2所述的基于联合信息熵的供水管网漏损识别的传感器优化布置方法,其特征在于:所述步骤(2)具体包括:
步骤(2.1)、利用校核好的管网水力模型对管网正常运行工况进行模拟,得到各候选节点压力模拟预测值y;
y=[y1,…,yi,…yn]
yi表示在正常工况下节点i处的压力模拟预测值;
步骤(2.2)、将正常工况下的压力模拟预测值与各种漏损情形下的压力模拟测量值进行比较,得到压力残差R=y-g;
Figure FDA0003653619470000022
其中rij为节点j处发生漏损,节点i处的水力模型预测值与测量值之差,R中的每一行对应于一个候选节点在各种漏损情形下的压力残差,R中的每一列对应于一种漏损情形在各候选节点的压力残差。
4.如权利要求3所述的基于联合信息熵的供水管网漏损识别的传感器优化布置方法,其特征在于所述步骤(3)具体如下:
(3.1)对于各候选节点i,根据区间界限(Iw,i)n将rij划分到不同的区间内,满足
Figure FDA0003653619470000031
Figure FDA0003653619470000032
其中NI是第i个候选节点划分得到的最大区间数;
(3.2)计算各个候选节点位置的熵H(R)i
Figure FDA0003653619470000033
其中,H(R)i是压力残差R在候选节点i处的熵,P(Rn)i是变量分布在第i个传感器第n个区间内的概率,NI为第i个传感器总的区间数;
(3.3)选择熵最大的节点放入传感器节点列表中,并保存当前的区间划分结果;
(3.4)在当前区间划分结果基础上,利用各候选节点对区间进一步划分,并根据步骤(3.2)中的方法计算各候选节点处的熵;
(3.5)对各候选节点的熵进行比较,选择熵最大的候选节点,并比较该候选节点的熵与传感器列表中最后一个节点的熵的大小,如果该候选节点的熵大于传感器列表最后一个节点的熵,则进入步骤(3.3),否则,进入步骤(3.6);
(3.6)将当前的传感器列表作为传感器优化布置方案。
5.如权利要求1所述的基于联合信息熵的供水管网漏损识别的传感器优化布置方法,其特征在于所述步骤(4.5)具体如下:
(4.5.1)对于第i种传感器优化布置方案,传感器数量为m;
(4.5.2)对于第i种方案,重复进行m次漏损识别实验,在第j次将第j个传感器“冷藏”,利用剩余的m-1个传感器对各种漏损情形进行识别,得到各种漏损情形下的识别率,并计算m次实验的平均漏损识别率;
(4.5.3)比较N种方案在发生故障时的漏损识别率,选择漏损识别率最高的方案作为具有鲁棒性的传感器优化布置方案。
CN202011318087.9A 2020-11-23 2020-11-23 基于联合信息熵的供水管网漏损识别传感器优化布置方法 Active CN112464421B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011318087.9A CN112464421B (zh) 2020-11-23 2020-11-23 基于联合信息熵的供水管网漏损识别传感器优化布置方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011318087.9A CN112464421B (zh) 2020-11-23 2020-11-23 基于联合信息熵的供水管网漏损识别传感器优化布置方法

Publications (2)

Publication Number Publication Date
CN112464421A CN112464421A (zh) 2021-03-09
CN112464421B true CN112464421B (zh) 2022-07-05

Family

ID=74799214

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011318087.9A Active CN112464421B (zh) 2020-11-23 2020-11-23 基于联合信息熵的供水管网漏损识别传感器优化布置方法

Country Status (1)

Country Link
CN (1) CN112464421B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113111513B (zh) * 2021-04-13 2024-04-12 上海商汤临港智能科技有限公司 传感器配置方案确定方法、装置、计算机设备及存储介质
CN113624915A (zh) * 2021-07-01 2021-11-09 同济大学 一种室内突发空气污染预警传感器布置优化方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111119282A (zh) * 2019-11-26 2020-05-08 中国地质大学(武汉) 一种针对供水管网的压力监测点优化布置方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE202008017243U1 (de) * 2008-05-02 2009-09-17 Rbs Wave Gmbh System zur Überwachung und/oder Ortung von Wasserverlusten in einem Wasserrohrleitungsnetz
CN106897679A (zh) * 2017-02-13 2017-06-27 长江水利委员会长江科学院 一种基于改进模糊c均值聚类的语义变化检测方法及系统
CN109783979B (zh) * 2019-02-15 2021-03-26 吉林大学 一种城市供水管网半监督条件下泄漏监测传感器布局优化方法
CN110043808B (zh) * 2019-05-29 2020-05-19 浙江大学 基于时间序列分析的供水管网漏损监测预警方法
CN110851774B (zh) * 2019-10-30 2024-05-14 杭州电子科技大学 一种城市供水管网的漏损估计方法
CN111853555B (zh) * 2020-07-07 2022-03-25 杭州电子科技大学 一种基于动态过程的供水管网暗漏识别方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111119282A (zh) * 2019-11-26 2020-05-08 中国地质大学(武汉) 一种针对供水管网的压力监测点优化布置方法

Also Published As

Publication number Publication date
CN112464421A (zh) 2021-03-09

Similar Documents

Publication Publication Date Title
CN112464421B (zh) 基于联合信息熵的供水管网漏损识别传感器优化布置方法
CN108984873B (zh) 供水管网实时漏损检测方法、设备、系统及存储介质
Bao et al. An analysis of availability for series Markov repairable system with neglected or delayed failures
KR101283828B1 (ko) 상수관망 진단 시스템
CN105912844A (zh) 一种天然气管网供气可靠性的评价方法及装置
CN111784116A (zh) 一种多能耦合下配电网弹性能力评估方法及系统
Liu et al. Lifecycle operational reliability assessment of water distribution networks based on the probability density evolution method
CN108615098A (zh) 基于贝叶斯生存分析的供水管网管道爆管风险预测方法
CN114462688A (zh) 一种基于lstm模型和动态阈值确定算法的爆管检测方法
Jun et al. Identifying pipes and valves of high importance for efficient operation and maintenance of water distribution systems
CN110848578A (zh) 一种基于pdd模型的城市供水管网既有漏损定位方法
CN110543108B (zh) 基于机理模型预测的供热系统测量数据的校正方法及系统
CN116307260B (zh) 一种面向缺陷路段扰动的城市路网韧性优化方法及系统
CN108536980B (zh) 一种考虑可靠性因素的气体探测器离散选址优化方法
CN116562012A (zh) 一种基于Kriging模型的动态可靠性分析方法
CN114169590B (zh) 一种基于多情景划分的水库入库径流预报校正方法和系统
CN113139584B (zh) 用于供水管网污染物入侵点识别的传感器优化布置方法
CN114154305A (zh) 一种用于供水管网漏损检测与识别的传感器优化布置多准则决策分析方法
CN114519259A (zh) 配电信息物理系统可靠性评估方法
Sun et al. Reliability prediction of long-lived linear assets with incomplete failure data
Ottenburger et al. MCDA-based genetic algorithms for developing disaster resilient designs of critical supply networks
CN110991839B (zh) 一种电-气综合能源系统整体脆弱性的评估方法
Cubillo et al. Water distribution system risk assessment method
CN114117790A (zh) 一种复杂混联系统的可靠性仿真评估方法和系统
CN114065937A (zh) 一种设备故障诊断方法、装置及终端设备

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant