CN110889594A - 一种城市管廊综合能源系统的脆弱性评估方法 - Google Patents
一种城市管廊综合能源系统的脆弱性评估方法 Download PDFInfo
- Publication number
- CN110889594A CN110889594A CN201911067108.1A CN201911067108A CN110889594A CN 110889594 A CN110889594 A CN 110889594A CN 201911067108 A CN201911067108 A CN 201911067108A CN 110889594 A CN110889594 A CN 110889594A
- Authority
- CN
- China
- Prior art keywords
- energy system
- comprehensive
- node
- demand
- comprehensive energy
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 96
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 71
- 239000003345 natural gas Substances 0.000 claims description 48
- 238000003860 storage Methods 0.000 claims description 41
- 230000008878 coupling Effects 0.000 claims description 22
- 238000010168 coupling process Methods 0.000 claims description 22
- 238000005859 coupling reaction Methods 0.000 claims description 22
- 239000007789 gas Substances 0.000 claims description 17
- 239000000126 substance Substances 0.000 claims description 14
- 238000004891 communication Methods 0.000 claims description 12
- 238000010248 power generation Methods 0.000 claims description 12
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000003287 optical effect Effects 0.000 claims 1
- 230000007123 defense Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 8
- 238000007726 management method Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 230000006378 damage Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000013486 operation strategy Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000006722 reduction reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/067—Enterprise or organisation modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/80—Management or planning
- Y02P90/82—Energy audits or management systems therefor
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种城市管廊综合能源系统的脆弱性评估方法,包括:S1、结合综合管廊和综合能源系统以及二者的关联性,建立城市管廊综合能源系统的脆弱性评估模型;S2、实时采集综合能源系统的中各能源需求节点的实际能源需求量、攻击成本和保护成本,输入到城市管廊综合能源系统的脆弱性评估模型中,得到城市管廊综合能源系统的脆弱性评估结果。采用通过结合综合管廊和综合能源系统,并考虑综合管廊和综合能源系统间的地理关联所建立的城市管廊综合能源系统的脆弱性评估模型,准确度较高。另外,本发明所提供的方法能够准确求解出城市管廊综合能源系统在蓄意攻击下最优保护策略,帮助城市管理者在有限的成本下做出最优的保护决策,减少系统的损失。
Description
技术领域
本发明属于关键基础设施韧性领域,更具体地,涉及一种城市管廊综合能源系统的脆弱性评估方法。
背景技术
综合管廊是城市管道的地下综合走廊,即在城市地下建造一个隧道空间,将各种工程管线集中于该空间内敷设,实施统一规划、统一设计、统一建设和统一管理,是保障城市运行的重要基础设施和生命线。如今,综合管廊在许多现代城市广泛建造。然而,修建综合管廊带来了新的威胁,城市综合管廊网络是城市综合能源系统的载体,由于综合管廊增加了城市基础设施系统之间的地理相互依赖性,当城市管廊遭受到破坏时,城市综合能源系统便会遭到破坏,能源是为人类一切生产活动提供能量的物质资源,因此综合管廊的破坏可能导致所有携带系统的破坏。
现有的城市管廊综合能源系统的脆弱性评估方法主要关注综合管廊的修建为城市带来的便捷,如方便集中管理和维修、减少线路管道的裸露、使得城市更加美观等。然而,现有的方法都忽略了城市综合管廊作为城市综合能源系统的载体,与城市基础设施系统之间存在地理相互依赖关系,当城市管廊的被破坏时,综合能源系统也会被极大破坏,脆弱性评估结果的准确率较低。
综上所述,提供一种准确率较高的城市管廊综合能源系统的脆弱性评估方法是亟待解决的问题。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种城市管廊综合能源系统的脆弱性评估方法,用于解决现有技术由于忽略城市综合管廊与城市基础设施系统之间存在地理相互依赖关系而导致的脆弱性评估准确率较低的问题。
为实现上述目的,本发明提供了一种城市管廊综合能源系统的脆弱性评估方法,包括以下步骤:
S1、结合综合管廊和综合能源系统以及二者的关联性,建立城市管廊综合能源系统的脆弱性评估模型;
S2、实时采集综合能源系统的中各能源需求节点的实际能源需求量、攻击成本和保护成本,输入到城市管廊综合能源系统的脆弱性评估模型中,得到城市管廊综合能源系统的脆弱性评估结果。
进一步优选地,步骤S1中建立城市管廊综合能源系统的脆弱性评估模型的方法,包括以下步骤:
S11、实时采集综合管廊的结构信息数据、综合能源系统的结构数据和供给需求数据以及综合管廊和能源系统的耦合信息;
S12、基于所采集的综合管廊的结构信息数据、综合能源系统的结构数据和供给需求数据,构建综合管廊和综合能源系统的网络模型;
S13、根据综合管廊和综合能源系统的网络模型和耦合信息、实时攻击成本和保护成本,构建攻击约束和保护约束;
S14、根据综合管廊和综合能源系统的网络模型和耦合信息,确定综合管廊和综合能源系统的地理关联约束;
S15、根据所得综合能源系统的网络结构,确定综合能源系统中各能源系统的运作模型,构建各能源系统的运作约束;
S16、根据所得综合能源系统的网络结构,以及综合能源系统中各能源系统的运作模型,构建综合能源系统中各能源系统间的关联关系约束;
S17、根据系统管理者的目标需求,并结合所得综合管廊和综合能源系统的地理关联约束、攻击约束、保护约束、关联关系约束以及各能源系统的运作约束,构建城市管廊综合能源系统的脆弱性评估模型。
进一步优选地,综合管廊的结构信息数据包括综合管廊中每个节点的位置、综合管廊中每一管廊段的起始节点和终止节点位置;综合能源系统的结构信息数据包括各能源系统中每个节点的位置以及每一条边的起始节点和终止节点位置。
进一步优选地,综合能源系统包括电力系统、天然气系统、水利系统、通讯系统。
进一步优选地,当综合能源系统中存在电力系统时,该能源系统的节点表示发电节点和电力需求节点,能源系统中的边表示每一段电线;此时,综合能源系统的结构信息数据包括各发电节点和电力需求节点的位置,每一段电线的起始节点和终止节点位置,每一段电线的电容、电阻、电抗和电感;综合能源系统的供给需求数据包括各发电节点的发电能力和成本,各需求节点的实际电力需求量;
当综合能源系统中存在天然气系统时,该能源系统的节点表示气井、天然气需求节点和天然气存储节点,能源系统中的边表示每一段天然气管道;此时,综合能源系统的结构数据包括各气井、天然气需求节点和天然气存储节点的位置,各加压站和阀门的位置,每一段天然气管道的起始节点和终止节点位置,各气井的出气能力和成本,每一段天然气管道的容量,各加压站的加压能力,各天然气存储节点的存储能力;综合能源系统的供给需求数据包括各气井、天然气需求节点和天然气存储节点的气压,各天然气需求节点的实际天然气需求量,各天然气存储节点的实际天然气存储量;
当综合能源系统中存在水利系统时,该能源系统的节点表示各水泵、水利需求节点和水利存储节点,能源系统中的边表示每一段水管;此时,综合能源系统的结构数据包括各水泵、水利需求节点和水利存储节点的位置,各加压站和阀门的位置,每一段水管的容量及其起始节点和终止节点位置,各水泵的出水能力和成本,各水利存储节点的存储能力,每一段水管的容量,各加压站的加压能力;综合能源系统的供给需求数据包括各水泵、水利需求节点和水利存储节点的水压、各水利需求节点的实际水利需求量,各水利存储节点的实际水利存储量;
当综合能源系统中存在通讯系统时,能源系统的节点表示信息站和通讯需求点,能源系统中的边表示每一段信道;此时,综合能源系统的结构数据包括各信息站和通讯需求点的位置,各信道的起始节点和终止节点位置。
进一步优选地,综合管廊和综合能源系统的耦合信息包括综合管廊和能源系统之间的位置耦合关系和综合能源系统中各能源系统之间的能源转换关系。
进一步优选地,综合管廊和综合能源系统的地理关联约束为:
其中,为能源系统k中的边e的故障状态,为与能源系统k中的边e存在耦合关系的综合管廊节点n的被保护状态,为与能源系统k中的边e存在耦合关系的综合管廊节点n的被攻击状态,为与能源系统k中的边e完全耦合的综合管廊的边e'的被保护状态,为与能源系统k中的边e完全耦合的综合管廊的边e'的被攻击状态,为与综合管廊的边e'部分耦合的能源系统k中的边e的被保护状态,为与综合管廊的边e'部分耦合的能源系统k中的边e的被攻击状态,其中,取值为1表示被保护,取值为0表示未被保护,取值为0表示被攻击,取值为1表示未被攻击。
进一步优选地,上述城市管廊综合能源系统的脆弱性评估模型为:
其中,vulnerability为城市管廊综合能源系统的脆弱性评估结果,为能源系统k中的节点n的被保护状态,为与综合管廊的边e'部分耦合的能源系统k中的边e的被保护状态,为综合管廊的中的节点n'的被保护状态,为综合管廊的中的节点e'的被保护状态,为能源系统k中节点n实际能源供给量,为能源系统k中第n个节点的需求,为能源系统k中边e上边的流量,为能源系统k中的节点i和能源系统l中的节点j之间的关联关系,ηk为能源供应系统k中需求未满足的单位惩罚成本,为对能源系统k中第n个节点的最大需求量,上述参数满足综合管廊和综合能源系统的地理关联约束、攻击约束、保护约束、关联关系约束以及各能源系统的运作约束。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有益效果:
1、本发明提出了一种城市管廊综合能源系统的脆弱性评估方法,由于城市综合管廊作为城市综合能源系统的载体,与城市基础设施系统之间存在地理相互依赖关系,当城市管廊的被破坏时,综合能源系统也会被极大破坏,通过结合综合管廊和综合能源系统,并考虑综合管廊和综合能源系统间的地理关联,建立城市管廊综合能源系统的脆弱性评估模型,相比现有的管廊综合能源系统的脆弱性评估方法,脆弱性评估结果准确率较高。
2、本发明所提出的一种城市管廊综合能源系统的脆弱性评估方法,将博弈策略应用到城市管廊网络综合能源系统,能够准确求解出城市管廊网络综合能源系统在蓄意攻击下最优保护策略,能够帮助城市管理者在有限的成本下做出最优的保护决策,减少系统的损失。
3、本发明所提出的一种城市管廊综合能源系统的脆弱性评估方法,可以在不同的攻击成本、保护成本和不同目标需求下制定出不同的攻击策略和保护策略,是一种在不同成本预算和目标需求下都普适的模型。
4、本发明所提出的一种城市管廊综合能源系统的脆弱性评估方法,考虑了不同能源系统间的物理关联,对综合能源系统进行整体优化,与不考虑能源系统间的物理关联的模型相比,这种方式可实现综合能源系统整体损失的最小化。
附图说明
图1是本发明所提供的一种城市管廊综合能源系统的脆弱性评估方法流程图;
图2是本发明实施例所提供的天津市综合管廊和能源系统结构图;
图3是本发明实施例所提供的天津市综合能源系统结构图;其中,图(a)为电力系统结构图,图(b)为水利系统结构图;
图4是分别在不考虑综合管廊和能源系统的地理关联关系和考虑综合管廊和能源系统的地理关联关系的情况下,所得脆弱性评估结果随攻击成本和防御成本的变化曲线;其中,图(a)为不考虑综合管廊和能源系统的地理关联关系的情况下,所得脆弱性评估结果随攻击成本和防御成本的变化曲线,图(b)为考虑综合管廊和能源系统的地理关联关系的情况下,所得脆弱性评估结果随攻击成本和防御成本的变化曲线。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
为实现上述目的,本发明提供了一种城市管廊综合能源系统的脆弱性评估方法,如图1所示,包括以下步骤:
S1、结合综合管廊和综合能源系统以及二者的关联性,建立城市管廊综合能源系统的脆弱性评估模型;
具体的,包括以下步骤:
S11、实时采集综合管廊的结构信息数据、综合能源系统的结构数据和供给需求数据以及综合管廊和能源系统的耦合信息;
具体的,综合管廊的结构信息数据包括综合管廊中每个节点的位置、综合管廊中每一管廊段的起始节点和终止节点位置。综合能源系统的结构信息数据包括各能源系统中每个节点的位置以及每一条边的起始节点和终止节点位置;综合管廊和综合能源系统的耦合信息包括综合管廊和能源系统之间的位置耦合关系和综合能源系统中各能源系统之间的能源转换关系。
具体的,综合能源系统包括电力系统、天然气系统、水利系统、通讯系统。进一步地,当综合能源系统中存在电力系统时,该能源系统的节点表示发电节点和电力需求节点,能源系统中的边表示每一段电线;此时,综合能源系统的结构信息数据包括各发电节点和电力需求节点的位置,每一段电线的起始节点和终止节点位置,每一段电线的电容、电阻、电抗和电感;综合能源系统的供给需求数据包括各发电节点的发电能力和成本,各需求节点的实际电力需求量;
当综合能源系统中存在天然气系统时,该能源系统的节点表示气井、天然气需求节点和天然气存储节点,能源系统中的边表示每一段天然气管道;此时,综合能源系统的结构数据包括各气井、天然气需求节点和天然气存储节点的位置,各加压站和阀门的位置,每一段天然气管道的起始节点和终止节点位置,各气井的出气能力和成本,每一段天然气管道的容量,各加压站的加压能力,各天然气存储节点的存储能力;综合能源系统的供给需求数据包括各气井、天然气需求节点和天然气存储节点的气压,各天然气需求节点的实际天然气需求量,各天然气存储节点的实际天然气存储量;
当综合能源系统中存在水利系统时,该能源系统的节点表示各水泵、水利需求节点和水利存储节点,能源系统中的边表示每一段水管;此时,综合能源系统的结构数据包括各水泵、水利需求节点和水利存储节点的位置,各加压站和阀门的位置,每一段水管的容量及其起始节点和终止节点位置,各水泵的出水能力和成本,各水利存储节点的存储能力,每一段水管的容量,各加压站的加压能力;综合能源系统的供给需求数据包括各水泵、水利需求节点和水利存储节点的水压、各水利需求节点的实际水利需求量,各水利存储节点的实际水利存储量。
当综合能源系统中存在通讯系统时,能源系统的节点表示信息站和通讯需求点,能源系统中的边表示每一段信道;此时,综合能源系统的结构数据包括各信息站和通讯需求点的位置,各信道的起始节点和终止节点位置。
S12、基于所采集的综合管廊的结构信息数据、综合能源系统的结构数据和供给需求数据,构建综合管廊和综合能源系统的网络模型;
具体的,构建综合管廊网络为GT(NT,ET),其中NT为管廊节点点集,ET为管廊段边集。
具体的,综合能源系统κ由一组能源系统构成,对每个能源系统k∈κ,构建综合能源系统网络为Gk(Nk,Ek),其中,Nk为能源系统k的节点点集,Ek为能源系统k的边集。对于每个能源系统k∈κ,存在一组源节点和一组需求节点
S13、根据综合管廊和综合能源系统的网络模型和耦合信息、实时攻击成本和保护成本,构建攻击约束和保护约束;
具体的,根据综合管廊和综合能源系统的网络模型和耦合信息确定得到攻击元件和保护元件,均包括四种类型的组件,分别对应为每个能源系统中的可破坏节点未完全敷设在综合管廊中的可破坏边(其中为每个能源系统中的可破坏的边)、可破坏的综合管廊节点以及可破坏的管廊边对于以上四种组件,分别定义二进制变量来表示每种类型的组件是否被保护,若被保护,则取值为1,否则取值为0。分别定义二进制变量来表示每种类型的组件是否被攻击,若被攻击,则取值为0,否则取值为1。
具体的,为了实现攻击成本和保护成本之间的博弈,所得保护约束为:
其中,为保护能源系统k中的节点n的成本,为能源系统k中的可破坏节点的被保护状态,为保护能源系统k中的未完全敷设于管廊中的边e成本,为能源系统k中未完全敷设在综合管廊中的可破坏边e的被保护状态,为保护可破坏的综合管廊节点n'的成本,为可破坏的综合管廊节点n'的被保护状态,为保护可破坏的管廊边e'的成本,为可破坏的管廊边e'的被保护状态,BD为总保护成本,其中,被保护状态w∈{0,1},当w取值为1时,表示被保护状态,当w取值为0时,表示未被保护状态。
所得攻击约束为:
其中,为攻击能源系统k中的节点n的成本,为能源系统k中的可破坏节点的被攻击状态,为攻击能源系统k中的未完全敷设于管廊中的边e成本,为能源系统k中未完全敷设在综合管廊中的可破坏边e的被攻击状态,为攻击可破坏的综合管廊节点n'的成本,为可破坏的综合管廊节点n'的被攻击状态,为攻击可破坏的管廊边e'的成本,为可破坏的管廊边e'的被攻击状态,BA为总攻击成本,其中,被攻击状态y∈{0,1},当y取值为0时,表示被攻击状态,当y取值为1时,表示未被攻击状态。另外
S14、根据综合管廊和综合能源系统的网络模型和耦合信息,确定综合管廊和综合能源系统的地理关联约束;
具体的,综合管廊和综合能源系统的地理关联表示为对于每个能源系统k∈κ,若能源系统k中的边e经过了综合管廊的节点n',则的值为1,否则就为0;若能源系统k中的边e有一部分敷设在综合管廊的边e'中,则的值为1,否则就为0。由于在实际案例中,存在能源系统中的某些边的一部分铺设于综合管廊中,一部分位于综合管廊之外的情况,如果敷设在综合管廊内部的这部分边得到了保护,而管廊外部的部分边未被保护,则这个保护是无用的,这条边仍然可能被损坏且和未保护的情况具有相同的风险等级。为了解决这个问题,需要为能源系统k定义一个向量若能源系统k中的边e没有完全敷设在综合管廊的边e'内,则的值为1,否则为0。
具体的,对于能源系统的边,其故障状态取决于该边通过的所有综合管廊节点或边的故障状态以及未完全敷设在综合管廊中的能源系统的边的故障状态。如果此管廊节点或边或未完全敷设在综合管廊中的能源系统的边受到攻击且未被保护,则所有与其相关的能源系统的边都将被破坏。例如,对于能源系统k中的边e,当与其关联的任一管廊节点被攻击且未被保护时,能源系统k中的边e处于被攻击状态,即综上,能源系统的边的故障状态由该边经过的管廊节点、与该边完全耦合的管廊边和未完全敷设在综合管廊中的能源系统的边的故障状态共同决定。故综合管廊和综合能源系统的地理关联约束为:
其中,为能源系统k中的边e的故障状态,为与能源系统k中的边e存在耦合关系的综合管廊节点n的被保护状态,为与能源系统k中的边e存在耦合关系的综合管廊节点n的被攻击状态,为与能源系统k中的边e完全耦合的综合管廊的边e'的被保护状态,为与能源系统k中的边e完全耦合的综合管廊的边e'的被攻击状态,为与综合管廊的边e'部分耦合的能源系统k中的边e的被保护状态,为与综合管廊的边e'部分耦合的能源系统k中的边e的被攻击状态,其中,取值为1表示被保护,取值为0表示未被保护,取值为0表示被攻击,取值为1表示未被攻击。
S15、根据所得综合能源系统的网络结构,确定综合能源系统中各能源系统的运作模型,构建各能源系统的运作约束;
具体的,综合能源系统中各能源系统的运作模型均选择网络流模型,该模型的核心思想是,能源系统中每个节点都要满足流量守恒。具体的,综合能源系统中各能源系统的运作需满足以下约束:
(1)每个节点的流量守恒,表示如下:
(2)每个边的流量需在该边的最大流量范围内,表示如下:
其中,为能源系统k中边e上边的流量,表示能源系统k中边e上的容量,为能源系统k中边e的故障状态,为能源系统k中边e的起点节点的故障状态,为能源系统k中边e的终点节点的故障状态,ok(e)为能源系统k中边e的起点节点,dk(e)为能源系统k中边e的终点节点。
(3)每个节点的能量供应量需在该节点的最大能源供应量范围内,表示为:
(4)用户的需求满足量不能超过每个负载节点的所需需求量,表示为:
S16、根据所得综合能源系统的网络结构,以及综合能源系统中各能源系统的运作模型,构建综合能源系统中各能源系统间的关联关系约束;
其中,为能源系统k中第i个节点的需求,为能源系统k中的节点i和能源系统l中的节点j之间的关联关系,若取值为1,则表示从能源系统k中的节点i到能源系统l中的节点j的物理关联可以正常工作,若取值为0,则表示从能源系统k中的节点i到能源系统l中的节点j的物理关联不能正常工作,为对节点i的能源最大需求量,表示能源系统k中的节点i和能源系统l中的节点j形成的物理关联节点对的集合。
此外,对于节点对中的每个节点j,根据其类型,应有不同的约束。
如果节点j是源节点,则其输出能量受零或最大输出量的限制,需满足如下条件:
如果节点j是需求节点,其需求水平由零或实际需求限制,需满足如下条件:
如果节点j损坏,其所有连接的边都将损坏,并且这些边上的流量也应该为零,如下所示:
S17、根据系统管理者的目标需求,并结合所得综合管廊和综合能源系统的地理关联约束、攻击约束、保护约束、关联关系约束以及各能源系统的运作约束,构建城市管廊综合能源系统的脆弱性评估模型。
具体的,城市管廊综合能源系统的脆弱性评估模型为:
其中,vulnerability为城市管廊综合能源系统的脆弱性评估结果,为能源系统k中的节点n的被保护状态,为与综合管廊的边e'部分耦合的能源系统k中的边e的被保护状态,为综合管廊的中的节点n'的被保护状态,为综合管廊的中的节点e'的被保护状态,为能源系统k中节点n实际能源供给量,为能源系统k中第n个节点的需求,为能源系统k中边e上边的流量,为能源系统k中的节点i和能源系统l中的节点j之间的关联关系,ηk为能源供应系统k中需求未满足的单位惩罚成本,为对能源系统k中第n个节点的最大需求量,上述参数满足所得综合管廊和综合能源系统的地理关联约束、攻击约束、保护约束、关联关系约束以及各能源系统的运作约束。
具体的,此模型是一个三阶防御者-攻击者-防御者模型,引入虚拟的攻击者和防御者的概念。其中,攻击者试图在一定攻击成本内采取最有效的攻击策略来破坏系统,而防御者则通过采取事前防御策略(保护关键元件)或在事后通过寻求最优的系统运作策略来最小化攻击者的破坏所带来的弹性损失。在模型的第一阶段,防御者选择防御策略,而后第二阶段中,攻击者在了解了第一阶段防御者的防御策略后选择自己的攻击策略,最后,防御者在清楚了前期的防御策略和攻击策略后,调整综合能源供应系统的运作模式,来使得系统损失V最小化。另外,当系统管理者的目标需求不同,能源供应系统k中需求未满足的单位惩罚成本ηk不同,
S2、实时采集综合能源系统的中各能源需求节点的实际能源需求量、攻击成本和保护成本,输入到城市管廊综合能源系统的脆弱性评估模型中,得到城市管廊综合能源系统的脆弱性评估结果。
为了更清楚的说明本发明所提出的城市管廊综合能源系统的脆弱性评估方法,下面结合实施例进行详述:
以天津市生态城内部管廊综合能源系统为例进行说明,如图2所示为天津市综合管廊和综合能源系统结构图,其中综合能源系统包含两种能源系统,电力系统和水利系统,其系统结构图分别如图3中的图(a)和图(b)所示。
构建城市管廊综合能源系统的脆弱性评估模型如下:
保护约束:
攻击约束:
地理约束:
运作约束:
关联关系约束:
上述中各参数的物理含义同上所述,不做赘述。
分别在不考虑综合管廊和能源系统的地理关联关系和考虑综合管廊和能源系统的地理关联关系的情况下对城市管廊综合能源系统进行脆弱性评估,所得脆弱性评估结果随攻击成本和防御成本的变化曲线如图4所示,其中,图(a)为不考虑综合管廊和能源系统的地理关联关系的情况下,所得脆弱性评估结果随攻击成本和防御成本的变化曲线,图(b)为考虑综合管廊和能源系统的地理关联关系的情况下所得脆弱性评估结果随攻击成本和防御成本的变化曲线。其中,横坐标代表攻击成本(分别表示为1,2,3,4,5,6),不同曲线代表不同的防御成本(分别表示为0,1,2,5,9),纵坐标代表在相应的攻击成本和防御成本下求解模型所得脆弱性的值。对比两图可以看出,当考虑综合管廊结构带来的地理关联时,综合能源系统的脆弱性要明显大于不考虑此地理关联的情况,这说明了综合管廊的构建增加了城市基础设施系统之间的地理相互依赖性,使得城市综合能源系统在攻击下更加损失更大,这也说明了此模型在求解管廊综合能源系统脆弱性的准确性要比以往的模型提高了很多。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种城市管廊综合能源系统的脆弱性评估方法,其特征在于,包括以下步骤:
S1、结合综合管廊和综合能源系统以及二者的关联性,建立城市管廊综合能源系统的脆弱性评估模型;
S2、实时采集综合能源系统的中各能源需求节点的实际能源需求量、攻击成本和保护成本,输入到城市管廊综合能源系统的脆弱性评估模型中,得到城市管廊综合能源系统的脆弱性评估结果。
2.根据权利要求1所述的城市管廊综合能源系统的脆弱性评估方法,其特征在于,步骤S1所述的方法包括以下步骤:
S11、实时采集综合管廊的结构信息数据、综合能源系统的结构数据和供给需求数据以及综合管廊和能源系统的耦合信息;
S12、基于所采集的综合管廊的结构信息数据、综合能源系统的结构数据和供给需求数据,构建综合管廊和综合能源系统的网络模型;
S13、根据综合管廊和综合能源系统的网络模型和耦合信息、实时攻击成本和保护成本,构建攻击约束和保护约束;
S14、根据综合管廊和综合能源系统的网络模型和耦合信息,确定综合管廊和综合能源系统的地理关联约束;
S15、根据所得综合能源系统的网络结构,确定综合能源系统中各能源系统的运作模型,构建各能源系统的运作约束;
S16、根据所得综合能源系统的网络结构,以及综合能源系统中各能源系统的运作模型,构建综合能源系统中各能源系统间的关联关系约束;
S17、根据系统管理者的目标需求,并结合所得综合管廊和综合能源系统的地理关联约束、攻击约束、保护约束、关联关系约束以及各能源系统的运作约束,构建城市管廊综合能源系统的脆弱性评估模型。
3.根据权利要求2所述的城市管廊综合能源系统的脆弱性评估方法,其特征在于,所述综合管廊的结构信息数据包括综合管廊中每个节点的位置、综合管廊中每一管廊段的起始节点和终止节点位置;所述综合能源系统的结构信息数据包括各能源系统中每个节点的位置以及每一条边的起始节点和终止节点位置。
4.根据权利要求3所述的城市管廊综合能源系统的脆弱性评估方法,其特征在于,所述综合能源系统包括电力系统、天然气系统、水利系统、通讯系统。
5.根据权利要求4所述的城市管廊综合能源系统的脆弱性评估方法,其特征在于,当综合能源系统中存在电力系统时,该能源系统的节点表示发电节点和电力需求节点,能源系统中的边表示每一段电线;此时,综合能源系统的结构信息数据包括各发电节点和电力需求节点的位置,每一段电线的起始节点和终止节点位置,每一段电线的电容、电阻、电抗和电感;综合能源系统的供给需求数据包括各发电节点的发电能力和成本,各需求节点的实际电力需求量;
当综合能源系统中存在天然气系统时,该能源系统的节点表示气井、天然气需求节点和天然气存储节点,能源系统中的边表示每一段天然气管道;此时,综合能源系统的结构数据包括各气井、天然气需求节点和天然气存储节点的位置,各加压站和阀门的位置,每一段天然气管道的起始节点和终止节点位置,各气井的出气能力和成本,每一段天然气管道的容量,各加压站的加压能力,各天然气存储节点的存储能力;综合能源系统的供给需求数据包括各气井、天然气需求节点和天然气存储节点的气压,各天然气需求节点的实际天然气需求量,各天然气存储节点的实际天然气存储量;
当综合能源系统中存在水利系统时,该能源系统的节点表示各水泵、水利需求节点和水利存储节点,能源系统中的边表示每一段水管;此时,综合能源系统的结构数据包括各水泵、水利需求节点和水利存储节点的位置,各加压站和阀门的位置,每一段水管的容量及其起始节点和终止节点位置,各水泵的出水能力和成本,各水利存储节点的存储能力,每一段水管的容量,各加压站的加压能力;综合能源系统的供给需求数据包括各水泵、水利需求节点和水利存储节点的水压、各水利需求节点的实际水利需求量,各水利存储节点的实际水利存储量;
当综合能源系统中存在通讯系统时,能源系统的节点表示信息站和通讯需求点,能源系统中的边表示每一段信道;此时,综合能源系统的结构数据包括各信息站和通讯需求点的位置,各信道的起始节点和终止节点位置。
6.根据权利要求2所述的城市管廊综合能源系统的脆弱性评估方法,其特征在于,所述综合管廊和综合能源系统的耦合信息包括综合管廊和能源系统之间的位置耦合关系和综合能源系统中各能源系统之间的能源转换关系。
7.根据权利要求2所述的城市管廊综合能源系统的脆弱性评估方法,其特征在于,所述综合管廊和综合能源系统的地理关联约束为:
8.根据权利要求2所述的城市管廊综合能源系统的脆弱性评估方法,其特征在于,所述城市管廊综合能源系统的脆弱性评估模型为:
其中,vulnerability为城市管廊综合能源系统的脆弱性评估结果,为能源系统k中的节点n的被保护状态,为与综合管廊的边e'部分耦合的能源系统k中的边e的被保护状态,为综合管廊的中的节点n'的被保护状态,为综合管廊的中的节点e'的被保护状态,为能源系统k中节点n实际能源供给量,为能源系统k中第n个节点的需求,为能源系统k中边e上边的流量,为能源系统k中的节点i和能源系统l中的节点j之间的关联关系,ηk为能源供应系统k中需求未满足的单位惩罚成本,为对能源系统k中第n个节点的最大需求量,上述参数满足综合管廊和综合能源系统的地理关联约束、攻击约束、保护约束、关联关系约束以及各能源系统的运作约束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911067108.1A CN110889594B (zh) | 2019-11-04 | 2019-11-04 | 一种城市管廊综合能源系统的脆弱性评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911067108.1A CN110889594B (zh) | 2019-11-04 | 2019-11-04 | 一种城市管廊综合能源系统的脆弱性评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110889594A true CN110889594A (zh) | 2020-03-17 |
CN110889594B CN110889594B (zh) | 2022-04-22 |
Family
ID=69746858
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911067108.1A Active CN110889594B (zh) | 2019-11-04 | 2019-11-04 | 一种城市管廊综合能源系统的脆弱性评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110889594B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111738457A (zh) * | 2020-05-12 | 2020-10-02 | 中国兵器科学研究院 | 确定对城市供水管网进行维护的方法和设备 |
CN111832907A (zh) * | 2020-06-22 | 2020-10-27 | 华中科技大学 | 关联电力–天然气系统在不同故障下的脆弱性评估方法 |
CN113420937A (zh) * | 2021-07-05 | 2021-09-21 | 山西大学 | 一种考虑输送路径脆弱性的综合能源切换方法及系统 |
CN114338426A (zh) * | 2021-12-31 | 2022-04-12 | 电子科技大学 | 一种多状态关联基础设施网络弹性优化方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103413015A (zh) * | 2013-04-24 | 2013-11-27 | 重庆科技学院 | 一种城市燃气管网弱性评估模型建立的方法 |
US20170185698A1 (en) * | 2015-12-23 | 2017-06-29 | Gas Technology Institute | Utility situational awareness system |
CN107909276A (zh) * | 2017-11-20 | 2018-04-13 | 广东电网有限责任公司电力调度控制中心 | 一种电力信息物理融合系统的脆弱性评估方法 |
EP3429049A1 (en) * | 2017-07-11 | 2019-01-16 | Siemens Aktiengesellschaft | A method and apparatus for controlling of utilities in a multimodal utility network |
CN109409729A (zh) * | 2018-10-19 | 2019-03-01 | 中国环境科学研究院 | 油气管道周边城市脆弱性评价方法 |
-
2019
- 2019-11-04 CN CN201911067108.1A patent/CN110889594B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103413015A (zh) * | 2013-04-24 | 2013-11-27 | 重庆科技学院 | 一种城市燃气管网弱性评估模型建立的方法 |
US20170185698A1 (en) * | 2015-12-23 | 2017-06-29 | Gas Technology Institute | Utility situational awareness system |
EP3429049A1 (en) * | 2017-07-11 | 2019-01-16 | Siemens Aktiengesellschaft | A method and apparatus for controlling of utilities in a multimodal utility network |
CN107909276A (zh) * | 2017-11-20 | 2018-04-13 | 广东电网有限责任公司电力调度控制中心 | 一种电力信息物理融合系统的脆弱性评估方法 |
CN109409729A (zh) * | 2018-10-19 | 2019-03-01 | 中国环境科学研究院 | 油气管道周边城市脆弱性评价方法 |
Non-Patent Citations (3)
Title |
---|
HUI ZHANG: ""Simplified operation models of integrated power and gas systems for vulnerability analysis"", 《PHYSICA A: STATISTICAL MECHANICS AND ITS APPLICATIONS》 * |
刘菁等: ""城市供水管网脆弱性研究进展"", 《灾害学》 * |
张志霞等: ""基于级联失效的燃气管网节点重要度评估"", 《消防科学与技术》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111738457A (zh) * | 2020-05-12 | 2020-10-02 | 中国兵器科学研究院 | 确定对城市供水管网进行维护的方法和设备 |
CN111738457B (zh) * | 2020-05-12 | 2021-03-23 | 中国兵器科学研究院 | 确定对城市供水管网进行维护的方法和设备 |
CN111832907A (zh) * | 2020-06-22 | 2020-10-27 | 华中科技大学 | 关联电力–天然气系统在不同故障下的脆弱性评估方法 |
CN113420937A (zh) * | 2021-07-05 | 2021-09-21 | 山西大学 | 一种考虑输送路径脆弱性的综合能源切换方法及系统 |
CN113420937B (zh) * | 2021-07-05 | 2022-12-30 | 山西大学 | 一种考虑输送路径脆弱性的综合能源切换方法及系统 |
CN114338426A (zh) * | 2021-12-31 | 2022-04-12 | 电子科技大学 | 一种多状态关联基础设施网络弹性优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110889594B (zh) | 2022-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110889594B (zh) | 一种城市管廊综合能源系统的脆弱性评估方法 | |
CN106959608B (zh) | 一种基于聚类粒子群算法的供水管网渗漏优化控制方法 | |
Rowell et al. | Obtaining layout of water distribution systems | |
US9762473B2 (en) | Network optimization | |
CN111027730B (zh) | 一种基于阀门操作和在线用水计量的供水管网漏失高效定位方法 | |
Gavanelli et al. | An ASP approach for the valves positioning optimization in a water distribution system | |
CN113449993B (zh) | 一种城市水源供水调度方法 | |
CN113803646A (zh) | 适用于缺陷管道排查的排水管网监测布点方法 | |
CN108708428A (zh) | 一种带有监测系统的水管路径规划方法 | |
Suzuki et al. | Water supply planning and operation in the Metropolitan region of Saõ Paulo: worst drought in history, conflicts, response, and resilience | |
CN115795997A (zh) | 一种市政消火栓取水组合优化方法 | |
CN103020739B (zh) | 一种基于网络算法的报警预案处理方法及系统 | |
CN106058907A (zh) | 一种多端柔性直流输电系统中线路投入直流网络的方法 | |
CN113361054B (zh) | 电力信息物理系统路由优化方法及系统 | |
CN215053230U (zh) | 综合管廊 | |
Bobylev | Sustainability and vulnerability analysis of critical underground infrastructure | |
CN109918795A (zh) | 一种兼顾物理性能和可靠性的能量枢纽标准矩阵稳态建模方法 | |
Belle et al. | A mathematical framework for the optimal coupling of interdependent critical infrastructures | |
CN104348668A (zh) | 一种基于快速近非迭代的双树形可靠pon规划方法 | |
CN110984302B (zh) | 一种输水管网传感器部署定位系统及定位方法 | |
Giustolisi | Development of a modularity index for reliability assessment of isolation valve systems | |
CN114120630B (zh) | 一种考虑道路交通动态特征的应急电源配置方法 | |
Shihu et al. | Decision support system of water distribution network expansion | |
CN216286198U (zh) | 一种智慧水务管网地理信息管理平台系统 | |
Grebenyuk et al. | Search and Selection of Blocking Cross-sections in the Analysis of Vulnerability and Efficiency of Engineering Networks |
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 |