CN116128110A - 一种电热综合能源系统的多阶段韧性提升方法及终端 - Google Patents

一种电热综合能源系统的多阶段韧性提升方法及终端 Download PDF

Info

Publication number
CN116128110A
CN116128110A CN202211643557.8A CN202211643557A CN116128110A CN 116128110 A CN116128110 A CN 116128110A CN 202211643557 A CN202211643557 A CN 202211643557A CN 116128110 A CN116128110 A CN 116128110A
Authority
CN
China
Prior art keywords
fault
stage
electric heating
toughness
constraint
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.)
Pending
Application number
CN202211643557.8A
Other languages
English (en)
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.)
Taiyuan University of Technology
Original Assignee
Taiyuan University of Technology
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 Taiyuan University of Technology filed Critical Taiyuan University of Technology
Priority to CN202211643557.8A priority Critical patent/CN116128110A/zh
Publication of CN116128110A publication Critical patent/CN116128110A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • 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
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Economics (AREA)
  • Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • General Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Health & Medical Sciences (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Quality & Reliability (AREA)
  • Development Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Water Supply & Treatment (AREA)
  • Game Theory and Decision Science (AREA)
  • Primary Health Care (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种电热综合能源系统的多阶段韧性提升方法及终端,以失负荷量最小为目标建立目标函数,结合电热综合能源系统在多阶段的网络拓扑约束条件、配电系统运行模型的约束条件,确定区域供热系统运行模型的约束条件,求解电热综合能源系统的多阶段韧性提升方案,实现了配电系统与区域供热系统的多阶段协同韧性提升。相较于现有技术中孤立的考虑配电网重构,本发明中的协同重构不仅通过在热源之间重新分配热负荷,提升区域供热系统韧性,还通过优化调整区域供热系统供热结构以匹配配电系统网络拓扑变化,有效控制系统间的故障蔓延,提升园区级电热综合能源系统韧性,有利于减少系统整体失负荷量,保证系统安全稳定运行。

Description

一种电热综合能源系统的多阶段韧性提升方法及终端
技术领域
本发明涉及综合能源系统的运行控制技术领域,特别涉及一种电热综合能源系统的多阶段韧性提升方法及终端。
背景技术
飓风、暴风雪和冰灾等极端自然灾害的发生,会使得关键供能设施,如输电线、天然气和供热管道遭到大范围破坏,供能中断事故屡有发生。
极端自然灾害下,大面积供能中断的恶劣事故频繁发生,严重威胁综合能源系统安全可靠地运行。因此提升能源系统的韧性引起了研究者的广泛关注。在电力系统的韧性研究方面,韧性提升策略可按极端事件发生的时序划分为灾前防御策略、灾中响应策略和灾后恢复策略。在灾前,通过加固风险水平较高的元件和部署灵活性发电资源及检修资源等技术措施提升韧性;在灾中,通过对分布式电源、储能装置等灵活性资源进行调度,在确保非故障区域的正常运行,降低系统失负荷损失;在灾后,通过对检修、运行人员调度以及供电恢复等多重恢复任务的协同优化,实现故障元件的修复和系统韧性的提升。如何提升灾后系统韧性是当前重点关注的问题。
近年来,热电联产机组(CHP)、电锅炉(EB)等耦合设备的广泛应用促进了配电系统(PDS)和区域供热系统(DHS)的深度耦合。子系统间的复杂耦合特性已成为电-热综合能源系统韧性提升的关键影响因素,具体原因如下:1)PDS/DHS故障通过耦合元件传播到另一个系统,单一系统故障可能引发连锁故障;2)仅依靠单一系统故障隔离可能引发故障传播,例如,不合理的PDS开关操作可能会导致CHP机组热输出功率降低,造成不必要的热负荷损失;3)当PDS和DHS协同运行时,才能充分挖掘灵活性资源紧急控制的能力(如CHP机组快速调节能力)。因此,有必要进行电-热耦合系统协同故障恢复。
发明内容
本发明所要解决的技术问题是:提供一种电热综合能源系统的多阶段韧性提升方法及终端,能够提升电热综合能源系统的韧性,保证系统安全稳定地运行。
为了解决上述技术问题,本发明采用的技术方案为:
一种电热综合能源系统的多阶段韧性提升方法,包括步骤:
获取电热综合能源的失负荷量,以所述失负荷量最小为目标建立目标函数;
建立电热综合能源系统在退化阶段、隔离阶段以及恢复阶段的网络拓扑约束条件,确定配电系统运行模型的约束条件,确定区域供热系统运行模型的约束条件;
结合所述目标函数,所述电热综合能源系统、所述配电系统以及所述区域供热系统的约束条件,计算得到所述电热综合能源系统的多阶段韧性提升方案。
为了解决上述技术问题,本发明采用的另一种技术方案为:
一种电热综合能源系统的多阶段韧性提升终端,包括存储器、处理器以及存储在所述存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述的一种电热综合能源系统的多阶段韧性提升方法的各个步骤。
本发明的有益效果在于:以失负荷量最小为目标建立目标函数,结合电热综合能源系统在多阶段的网络拓扑约束条件、配电系统运行模型的约束条件,确定区域供热系统运行模型的约束条件,求解电热综合能源系统的多阶段韧性提升方案,实现了配电系统与区域供热系统的多阶段协同韧性提升。相较于现有技术中孤立的考虑配电网重构,本发明中的协同重构不仅通过在热源之间重新分配热负荷,提升区域供热系统韧性,还通过优化调整区域供热系统供热结构以匹配配电系统网络拓扑变化,有效控制系统间的故障蔓延,提升园区级电热综合能源系统韧性,有利于减少系统整体失负荷量,保证系统安全稳定运行。
附图说明
图1为本发明实施例的一种电热综合能源系统的多阶段韧性提升方法的流程图;
图2为本发明实施例的一种电热综合能源系统的多阶段韧性提升终端的示意图;
标号说明:
1、一种电热综合能源系统的多阶段韧性提升终端;2、存储器;3、处理器。
具体实施方式
为详细说明本发明的技术内容、所实现目的及效果,以下结合实施方式并配合附图予以说明。
请参照图1,本发明实施例提供了一种电热综合能源系统的多阶段韧性提升方法,包括步骤:
获取电热综合能源的失负荷量,以所述失负荷量最小为目标建立目标函数;
建立电热综合能源系统在退化阶段、隔离阶段以及恢复阶段的网络拓扑约束条件,确定配电系统运行模型的约束条件,确定区域供热系统运行模型的约束条件;
结合所述目标函数,所述电热综合能源系统、所述配电系统以及所述区域供热系统的约束条件,计算得到所述电热综合能源系统的多阶段韧性提升方案。
从上述描述可知,本发明的有益效果在于:以失负荷量最小为目标建立目标函数,结合电热综合能源系统在多阶段的网络拓扑约束条件、配电系统运行模型的约束条件,确定区域供热系统运行模型的约束条件,求解电热综合能源系统的多阶段韧性提升方案,实现了配电系统与区域供热系统的多阶段协同韧性提升。相较于现有技术中孤立的考虑配电网重构,本发明中的协同重构不仅通过在热源之间重新分配热负荷,提升区域供热系统韧性,还通过优化调整区域供热系统供热结构以匹配配电系统网络拓扑变化,有效控制系统间的故障蔓延,提升园区级电热综合能源系统韧性,有利于减少系统整体失负荷量,保证系统安全稳定运行。
进一步地,以所述失负荷量最小为目标建立目标函数包括:
建立电热综合能源的失负荷量最小的目标函数:
Figure BDA0004008750540000031
式中,
Figure BDA0004008750540000032
表示在故障阶段t电负荷的损失量,
Figure BDA0004008750540000033
表示在故障阶段t热负荷的损失量,aT表示电负荷权重,bT表示热负荷权重,Tt表示故障阶段t的开始时间,pc表示故障场景c出现的可能性,C表示故障场景集合。
由上述描述可知,建立电-热综合能源失负荷量最小的目标函数,便于后续电热综合能源系统的多阶段韧性提升方案的计算。
进一步地,以所述失负荷量最小为目标建立目标函数之后还包括:
使用第一韧性指标计算恢复阶段的总体失负荷率,使用第二韧性指标计算多阶段故障恢复过程的总体失负荷率;
所述第一韧性指标Rr,c为:
Figure BDA0004008750540000041
所述第二韧性指标Rc为:
Figure BDA0004008750540000042
由上述描述可知,通过计算韧性指标,能够验证协同重构对于系统韧性提升的有效性。
进一步地,所述建立电热综合能源系统在退化阶段的网络拓扑约束条件包括:
当闭合管道或者线路发生故障时,将闭合管道或者线路两端的节点划分为故障区;
将闭合管道或者线路两端的节点同时划分为故障区和非故障区;
当热电联产机组在配电系统被划分为故障区时,将区域供热系统中的热电联产机组划分为故障区,当热电联产机组在区域供热系统被划分为故障区时,将配电系统中的热电联产机组划分为故障区。
由上述描述可知,针对不同的工况进行故障区和非故障区的划分,能够得到退化阶段的网络拓扑约束条件,便于后续生成多阶段的韧性提升方案。
进一步地,所述建立电热综合能源系统在隔离阶段的网络拓扑约束条件包括:
当闭合管道或者线路发生故障时,将闭合管道或者线路断开;
当故障区的管道或者线路检修后,将开关和阀门用于故障隔离,非故障区的开关和阀门直接用于故障隔离;
将闭合管道或者线路两端的节点同时划分为故障区和非故障区;
当热电联产机组在配电系统被划分为故障区时,将区域供热系统中的热电联产机组划分为故障区,当热电联产机组在区域供热系统被划分为故障区时,将配电系统中的热电联产机组划分为故障区。
由上述描述可知,针对不同的工况进行故障隔离,能够得到隔离阶段的网络拓扑约束条件,便于后续生成多阶段的韧性提升方案。
进一步地,所述建立电热综合能源系统在恢复阶段的网络拓扑约束条件包括:
当故障区的管道或者线路检修后,将开关和阀门用于故障恢复,非故障区的开关和阀门直接用于故障恢复;
将配电系统和区域供热系统的网络拓扑设置为辐射状约束;
筛除隔离阶段划分的故障区和非故障区。
由上述描述可知,针对不同的工况进行故障恢复,且在隔离阶段划分的故障区-非故障区不允许再次重新相连,能够得到恢复阶段的网络拓扑约束条件,便于后续生成多阶段的韧性提升方案。
进一步地,所述确定配电系统运行模型的约束条件包括:
确定配电系统运行模型的节点功率平衡方程约束、支路容量约束、热电联产机组功率约束、分布式电源功率约束、节点失负荷约束以及节点电压约束。
进一步地,所述确定区域供热系统运行模型的约束条件包括:
确定区域供热系统中热电联产机组的运行方程约束、热电联产机组的热功率约束、热泵的热功率约束、热网管道热量损失方程约束、管道传输容量约束、节点热平衡约束以及节点失负荷约束。
由上述描述可知,通过上述约束,能够解决热网重构问题。
进一步地,结合所述目标函数,所述电热综合能源系统、所述配电系统以及所述区域供热系统的约束条件,计算得到所述电热综合能源系统的多阶段韧性提升方案包括:
利用内点法,在所述目标函数的基础上,根据所述电热综合能源系统、所述配电系统以及所述区域供热系统的约束条件,计算得到所述电热综合能源系统的多阶段韧性提升方案;
所述电热综合能源系统的多阶段韧性提升方案包括电热综合能源系统内配电系统、区域供热系统网络拓扑的输出功率。
由上述描述可知,利用内点法能够求解得到电-热综合能源系统协同韧性提升策略,实现配电系统与区域供热系统的多阶段协同韧性提升。
请参照图2,本发明另一实施例提供了一种电热综合能源系统的多阶段韧性提升终端,包括存储器、处理器以及存储在所述存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述的一种电热综合能源系统的多阶段韧性提升方法的各个步骤。
本发明上述的一种电热综合能源系统的多阶段韧性提升方法及终端,适用于电热综合能源系统韧性提升策略制定,适配于原本的电力系统与区域供热系统能量管理系统,有利于减少系统整体失负荷量,保证系统安全稳定运行,以下通过具体的实施方式进行说明:
实施例一
请参照图1,一种电热综合能源系统的多阶段韧性提升方法,包括步骤:
S1、获取电热综合能源的失负荷量,以所述失负荷量最小为目标建立目标函数。
S11、以所述失负荷量最小为目标建立目标函数包括:
建立电热综合能源的失负荷量最小的目标函数:
Figure BDA0004008750540000061
式中,
Figure BDA0004008750540000062
表示在故障阶段t电负荷的损失量,
Figure BDA0004008750540000063
表示在故障阶段t热负荷的损失量,aT表示电负荷权重,bT表示热负荷权重,Tt表示故障阶段t的开始时间,pc表示故障场景c出现的可能性,C表示故障场景集合。
S12、使用第一韧性指标计算恢复阶段的总体失负荷率,使用第二韧性指标计算多阶段故障恢复过程的总体失负荷率。
其中,所述第一韧性指标Rr,c为:
Figure BDA0004008750540000071
所述第二韧性指标Rc为:
Figure BDA0004008750540000072
式中,pL表示电负荷的损失量,hL表示热负荷的损失量。
通过上述韧性指标能够验证协同重构对于系统韧性提升的有效性。
S2、建立电热综合能源系统在退化阶段、隔离阶段以及恢复阶段的网络拓扑约束条件,确定配电系统运行模型的约束条件,确定区域供热系统运行模型的约束条件。
S21、建立电热综合能源系统在退化阶段的网络拓扑约束条件包括:
S211、当闭合管道或者线路发生故障时,将闭合管道或者线路两端的节点划分为故障区:
Figure BDA0004008750540000073
Figure BDA0004008750540000074
其中,kpipe表示配电系统线路集合,kline表示区域供热系统管道集合,C表示故障场景集合,t表示故障阶段,fij,c表示线路/管道(i,j)故障状态,μij,0表示线路/管道(i,j)初始连接状态,ni,c,t表示在故障阶段t配电/供热系统节点i故障状态,nj,c,t表表示在故障阶段t配电/供热系统节点j故障状态,i表示线路的起点,j表示线路的终点。
S212、将闭合管道或者线路两端的节点同时划分为故障区和非故障区:
Figure BDA0004008750540000075
Figure BDA0004008750540000076
S213、当热电联产机组在配电系统被划分为故障区时,将区域供热系统中的热电联产机组划分为故障区,当热电联产机组在区域供热系统被划分为故障区时,将配电系统中的热电联产机组划分为故障区:
Figure BDA0004008750540000077
式中,nm,c,t在故障阶段t配电网中CHP机组节点m故障状态,nn,c,t在故障阶段t热网中CHP机组节点n故障状态,
Figure BDA0004008750540000078
表示配电系统中热电联产机组CHP集合,
Figure BDA0004008750540000081
表示供热系统中热电联产机组CHP集合。
S22、建立电热综合能源系统在隔离阶段的网络拓扑约束条件包括:
S221、当闭合管道或者线路发生故障时,将闭合管道或者线路断开:
Figure BDA0004008750540000082
式中,sij,0表示线路/管道(i,j)开关配置状态,μij,c,t表示在故障阶段t线路/管道(i,j)连接状态。
S222、当故障区的管道或者线路检修后,将开关和阀门用于故障隔离,非故障区的开关和阀门直接用于故障隔离:
为了保证检修人员的安全,仅当故障区内管道/线路检修后,其配置开关/阀门可用于故障隔离,非故障区配置的开关/阀门可直接用于故障隔离:
Figure BDA0004008750540000083
Figure BDA0004008750540000084
S223、将闭合管道或者线路两端的节点同时划分为故障区和非故障区:
Figure BDA0004008750540000085
Figure BDA0004008750540000086
S224、当热电联产机组在配电系统被划分为故障区时,将区域供热系统中的热电联产机组划分为故障区,当热电联产机组在区域供热系统被划分为故障区时,将配电系统中的热电联产机组划分为故障区:
Figure BDA0004008750540000087
S23、建立电热综合能源系统在恢复阶段的网络拓扑约束条件包括:
S231、当故障区的管道或者线路检修后,将开关和阀门用于故障恢复,非故障区的开关和阀门直接用于故障恢复:
Figure BDA0004008750540000088
S232、将配电系统和区域供热系统的网络拓扑设置为辐射状约束:
Figure BDA0004008750540000089
Figure BDA00040087505400000810
式中,kbus表示配电系统节点集合,knd表示区域供热系统节点集合,ωij,c,t、ωji,c,t均表示辅助变量,表示在故障恢复阶段t线路(i,j)两端节点的父子关系,ωij,c,t=1表示节点i为节点j的父节点,ωsj,c,t表示在故障恢复阶段t线路(s,j)两端节点的父子关系,ωsj,c,t=1表示节点s为节点j的父节点,gj、pj分别表示节点j热站/变电站配置状态,γj,c,t、dj分别表示分布式电源的启用状态及配置状态,π(j)表示节点j为末端节点的线路首段节点集、δ(j)表示以节点j为首端节点的线路末端节点集。
S233、筛除隔离阶段划分的故障区和非故障区:
Figure BDA0004008750540000091
Figure BDA0004008750540000092
S24、确定配电系统运行模型的约束条件包括:
确定配电系统运行模型的节点功率平衡方程约束、支路容量约束、热电联产机组功率约束、分布式电源功率约束、节点失负荷约束以及节点电压约束。
具体的,节点功率平衡方程为:
Figure BDA0004008750540000093
Figure BDA0004008750540000094
Figure BDA0004008750540000095
式中,A表示电网的节点-支路矩阵,
Figure BDA0004008750540000096
分别表示修正后节点注入有功、无功功率,
Figure BDA0004008750540000097
分别表示修正后线路有功、无功潮流,
Figure BDA0004008750540000098
分别表示修正后CHP和DG输出有功、无功功率,
Figure BDA0004008750540000099
分别表示修正后节点负荷有功和无功功率需求及失负荷量,修正后的功率值为实际功率和节点电压的比值。
支路容量约束:断开/故障输电线路传输功率为零,闭合输电线路传输功率不允许超过其限值。
Figure BDA00040087505400000910
Figure BDA00040087505400000911
Figure BDA00040087505400000912
式中,AoB表示矩阵Hadamard乘积,即AoB=(aij*bij),F表示电网的送端节点-支路矩阵,
Figure BDA00040087505400000913
表示输电线路传输容量限值,
Figure BDA00040087505400000914
表示修正后节点电压值,为实际节点电压的导数,M表示足够大的常数。
CHP功率约束:故障区内CHP机组有功、无功功率为零,非故障区内CHP机组的有功、无功功率在设定的安全运行上、下限值之间:
Figure BDA0004008750540000101
Figure BDA0004008750540000102
Figure BDA0004008750540000103
式中,
Figure BDA0004008750540000104
表示CHP-电网节点矩阵,pCHP,
Figure BDA0004008750540000105
表示CHP有功功率上下限,q CHP,
Figure BDA0004008750540000106
表示CHP无功功率上下限,n表示热网节点所属区域,1表示该节点属于故障区,0表示该节点属于非故障区。
DG(分布式电源)功率约束:故障区内DG有功、无功功率为零,非故障区内DG的有功、无功功率在设定的安全运行上、下限值之间:
Figure BDA0004008750540000107
Figure BDA0004008750540000108
Figure BDA0004008750540000109
式中,
Figure BDA00040087505400001010
表示DG-电网节点矩阵,pDG,
Figure BDA00040087505400001011
表示DG有功功率上下限,q DG,
Figure BDA00040087505400001012
表示DG无功功率上下限。
节点失负荷约束:故障区内电负荷将全部失去,非故障区内节点失负荷量小于或等于节点负荷量:
Figure BDA00040087505400001013
Figure BDA00040087505400001014
Figure BDA00040087505400001015
Figure BDA00040087505400001016
节点电压约束:故障区内修正后节点电压幅值设定为1:
Figure BDA00040087505400001017
Figure BDA00040087505400001018
式中,
Figure BDA00040087505400001019
表示节点位于故障区时修正后电压幅值;
非故障区内节点电压幅值在设定的电力系统安全运行电压的上下限值
Figure BDA00040087505400001020
u之间:
Figure BDA0004008750540000111
Figure BDA0004008750540000112
Figure BDA0004008750540000113
式中,
Figure BDA0004008750540000114
表示节点位于非故障区时修正后电压幅值,
Figure BDA0004008750540000115
表示修正后节点电压幅值,
Figure BDA0004008750540000116
u表示节点允许电压上限和下限值。
闭合输电线路节点电压需满足支路潮流约束:
Figure BDA0004008750540000117
式中,Ae为表示电网节点-支路矩阵,r、x表示输电线路的电阻、电抗。
S25、确定配电系统运行模型的约束条件包括:
确定配电系统运行模型的节点功率平衡方程约束、支路容量约束、热电联产机组功率约束、分布式电源功率约束、节点失负荷约束以及节点电压约束。
具体的,设定区域供热系统安全运行的等式和不等式约束。由于热网重构所带来的水力工况变化,供热管网中的循环水流向为变量,因此无法确定混温节点,进而无法建立混温约束模型;由于考虑水力工况变化(流量可变),该模型包含大量双线性及指数项,非线性强;而热网重构本身会引入阀门开闭等0,1变量,因此该模型为混合整数强非线性(非凸)模型,这会给电-热综合能源系统多阶段协同故障恢复问题的求解和工程应用带来极大挑战。因此考虑能流(Energy Flow,EF)模型约束,包括:
电力系统与区域供热系统的耦合元件——区域供热系统中热电联产机组(CHP)的运行特性方程约束:
pCHP≤diag(χ 1χ 2,...χ n)hCHP,
pCHP≥diag(χ 1χ 2,...χ n)hCHP,
式中,pCHP表示CHP机组的有功功率,hCHP表示CHP机组的热功率,χ i表示第i台CHP机组热电比的倒数,CHP机组热电比可从CHP机组的出厂说明书中获取。
区域供热系统中CHP机组的热功率约束:
故障区内CHP机组热功率为零,非故障区内CHP机组的热功率在设定的安全运行上、下限值h CHP
Figure BDA0004008750540000118
之间:
Figure BDA0004008750540000121
式中,
Figure BDA0004008750540000122
表示CHP机组-热网节点关联矩阵,hCHP表示CHP机组的热功率安全运行的下限,
Figure BDA0004008750540000123
表示CHP机组的热功率安全运行的上限,n表示热网节点所属区域,1表示该节点属于故障区,0表示该节点属于非故障区。
区域供热系统中热泵(HB)热功率约束:
故障区内HB热功率为零,非故障区内HB的热功率在设定的安全运行上、下限值hHB
Figure BDA0004008750540000124
之间:
Figure BDA0004008750540000125
式中,
Figure BDA0004008750540000126
表示HB-热网节点关联矩阵,hHB表示HB热功率,hHB表示HB的热功率安全运行的下限,
Figure BDA0004008750540000127
表示HB的热功率安全运行的上限。
HB热功率和燃料消耗之间的关系:hHB=ηofHB
式中,fHB表示HB的燃料消耗量,η表示HB热功率与燃料消耗量的比值,η可从HB的出厂说明书中获取。
区域供热系统中热网管道热量损失方程约束:
Figure BDA0004008750540000128
式中,
Figure BDA0004008750540000129
表示热网管道始端循环水的可利用热功率,即供水管网中该管道始端工质流所含热功率与回水管网中对应工质流所含热功率之差,
Figure BDA00040087505400001210
表示管道始端及末端循环水的热功率损失。
管道传输容量约束:断开/故障管道传输热功率为零,闭合管道传输热功率不允许超过其限值:
Figure BDA00040087505400001211
式中,
Figure BDA00040087505400001212
表示管道始端及末端循环水的可利用热功率的上限,μij表示管道通断状态。
节点热平衡约束:
Figure BDA00040087505400001213
hHS=hCHP+hHB,
式中,Fh,Th表示区域供热系统节点-支路矩阵,
Figure BDA00040087505400001214
表示热站节点矩阵,hHS表示热站输出热功率,hL、hLoss分别表示节点热负荷需求量及失负荷量。
节点失负荷约束:故障区内电负荷将全部失去,非故障区内节点失负荷量小于或等于节点负荷量:
-M(1-n)≤hL-hLoss≤M(1-n),
0≤hLoss≤hL.
S3、结合所述目标函数,所述电热综合能源系统、所述配电系统以及所述区域供热系统的约束条件,计算得到所述电热综合能源系统的多阶段韧性提升方案。
具体的,利用内点法,在所述目标函数的基础上,根据所述电热综合能源系统、所述配电系统以及所述区域供热系统的约束条件,计算得到所述电热综合能源系统的多阶段韧性提升方案;
所述电热综合能源系统的多阶段韧性提升方案包括电热综合能源系统内配电系统、区域供热系统网络拓扑的输出功率以及变电站、DG、CHP机组、电锅炉输出功率。
因此,本实施例考虑电热系统的紧密耦合与相互影响,实现了配电系统与区域供热系统的多阶段协同韧性提升。相比于孤立的考虑配电网重构,协同重构不仅通过在热源之间重新分配热负荷,提升DHS韧性,还通过优化调整DHS供热结构以匹配PDS网络拓扑变化并有效系统间的故障蔓延,提升园区级电-热综合能源系统韧性。本实施例可以实际应用于电-热综合能源系统韧性提升策略制定,适配于原本的电力系统与区域供热系统能量管理系统,有利于减少系统整体失负荷量,保证系统安全稳定运行。
实施例二
请参照图2,一种电热综合能源系统的多阶段韧性提升终端1,包括存储器2、处理器3以及存储在所述存储器2上并可在处理器3上运行的计算机程序,所述处理器3执行所述计算机程序时实现实施例一的一种电热综合能源系统的多阶段韧性提升方法的各个步骤。
综上所述,本发明提供的一种电热综合能源系统的多阶段韧性提升方法及终端,以失负荷量最小为目标建立目标函数,结合电热综合能源系统在多阶段的网络拓扑约束条件、配电系统运行模型的约束条件,确定区域供热系统运行模型的约束条件,求解电热综合能源系统的多阶段韧性提升方案,实现了配电系统与区域供热系统的多阶段协同韧性提升。相较于现有技术中孤立的考虑配电网重构,本发明中的协同重构不仅通过在热源之间重新分配热负荷,提升区域供热系统韧性,还通过优化调整区域供热系统供热结构以匹配配电系统网络拓扑变化,有效控制系统间的故障蔓延,提升园区级电热综合能源系统韧性,有利于减少系统整体失负荷量,保证系统安全稳定运行。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等同变换,或直接或间接运用在相关的技术领域,均同理包括在本发明的专利保护范围内。

Claims (10)

1.一种电热综合能源系统的多阶段韧性提升方法,其特征在于,包括步骤:
获取电热综合能源的失负荷量,以所述失负荷量最小为目标建立目标函数;
建立电热综合能源系统在退化阶段、隔离阶段以及恢复阶段的网络拓扑约束条件,确定配电系统运行模型的约束条件,确定区域供热系统运行模型的约束条件;
结合所述目标函数,所述电热综合能源系统、所述配电系统以及所述区域供热系统的约束条件,计算得到所述电热综合能源系统的多阶段韧性提升方案。
2.根据权利要求1所述的一种电热综合能源系统的多阶段韧性提升方法,其特征在于,以所述失负荷量最小为目标建立目标函数包括:
建立电热综合能源的失负荷量最小的目标函数:
Figure FDA0004008750530000011
式中,
Figure FDA0004008750530000012
表示在故障阶段t电负荷的损失量,
Figure FDA0004008750530000013
表示在故障阶段t热负荷的损失量,aT表示电负荷权重,bT表示热负荷权重,Tt表示故障阶段t的开始时间,pc表示故障场景c出现的可能性,C表示故障场景集合。
3.根据权利要求2所述的一种电热综合能源系统的多阶段韧性提升方法,其特征在于,以所述失负荷量最小为目标建立目标函数之后还包括:
使用第一韧性指标计算恢复阶段的总体失负荷率,使用第二韧性指标计算多阶段故障恢复过程的总体失负荷率;
所述第一韧性指标Rrc为:
Figure FDA0004008750530000014
所述第二韧性指标Rc为:
Figure FDA0004008750530000015
4.根据权利要求1所述的一种电热综合能源系统的多阶段韧性提升方法,其特征在于,所述建立电热综合能源系统在退化阶段的网络拓扑约束条件包括:
当闭合管道或者线路发生故障时,将闭合管道或者线路两端的节点划分为故障区;
将闭合管道或者线路两端的节点同时划分为故障区和非故障区;
当热电联产机组在配电系统被划分为故障区时,将区域供热系统中的热电联产机组划分为故障区,当热电联产机组在区域供热系统被划分为故障区时,将配电系统中的热电联产机组划分为故障区。
5.根据权利要求1所述的一种电热综合能源系统的多阶段韧性提升方法,其特征在于,所述建立电热综合能源系统在隔离阶段的网络拓扑约束条件包括:
当闭合管道或者线路发生故障时,将闭合管道或者线路断开;
当故障区的管道或者线路检修后,将开关和阀门用于故障隔离,非故障区的开关和阀门直接用于故障隔离;
将闭合管道或者线路两端的节点同时划分为故障区和非故障区;
当热电联产机组在配电系统被划分为故障区时,将区域供热系统中的热电联产机组划分为故障区,当热电联产机组在区域供热系统被划分为故障区时,将配电系统中的热电联产机组划分为故障区。
6.根据权利要求5所述的一种电热综合能源系统的多阶段韧性提升方法,其特征在于,所述建立电热综合能源系统在恢复阶段的网络拓扑约束条件包括:
当故障区的管道或者线路检修后,将开关和阀门用于故障恢复,非故障区的开关和阀门直接用于故障恢复;
将配电系统和区域供热系统的网络拓扑设置为辐射状约束;
筛除隔离阶段划分的故障区和非故障区。
7.根据权利要求1所述的一种电热综合能源系统的多阶段韧性提升方法,其特征在于,所述确定配电系统运行模型的约束条件包括:
确定配电系统运行模型的节点功率平衡方程约束、支路容量约束、热电联产机组功率约束、分布式电源功率约束、节点失负荷约束以及节点电压约束。
8.根据权利要求1所述的一种电热综合能源系统的多阶段韧性提升方法,其特征在于,所述确定区域供热系统运行模型的约束条件包括:
确定区域供热系统中热电联产机组的运行方程约束、热电联产机组的热功率约束、热泵的热功率约束、热网管道热量损失方程约束、管道传输容量约束、节点热平衡约束以及节点失负荷约束。
9.根据权利要求1所述的一种电热综合能源系统的多阶段韧性提升方法,其特征在于,结合所述目标函数,所述电热综合能源系统、所述配电系统以及所述区域供热系统的约束条件,计算得到所述电热综合能源系统的多阶段韧性提升方案包括:
利用内点法,在所述目标函数的基础上,根据所述电热综合能源系统、所述配电系统以及所述区域供热系统的约束条件,计算得到所述电热综合能源系统的多阶段韧性提升方案;
所述电热综合能源系统的多阶段韧性提升方案包括电热综合能源系统内配电系统、区域供热系统网络拓扑的输出功率。
10.一种电热综合能源系统的多阶段韧性提升终端,包括存储器、处理器以及存储在所述存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现上述权利要求1至9任一项所述的一种电热综合能源系统的多阶段韧性提升方法的各个步骤。
CN202211643557.8A 2022-12-20 2022-12-20 一种电热综合能源系统的多阶段韧性提升方法及终端 Pending CN116128110A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211643557.8A CN116128110A (zh) 2022-12-20 2022-12-20 一种电热综合能源系统的多阶段韧性提升方法及终端

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211643557.8A CN116128110A (zh) 2022-12-20 2022-12-20 一种电热综合能源系统的多阶段韧性提升方法及终端

Publications (1)

Publication Number Publication Date
CN116128110A true CN116128110A (zh) 2023-05-16

Family

ID=86303685

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211643557.8A Pending CN116128110A (zh) 2022-12-20 2022-12-20 一种电热综合能源系统的多阶段韧性提升方法及终端

Country Status (1)

Country Link
CN (1) CN116128110A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116894342A (zh) * 2023-07-19 2023-10-17 天津大学 基于天然气网络管存的电-气综合能源系统韧性提升方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116894342A (zh) * 2023-07-19 2023-10-17 天津大学 基于天然气网络管存的电-气综合能源系统韧性提升方法
CN116894342B (zh) * 2023-07-19 2024-03-12 天津大学 基于天然气网络管存的电-气综合能源系统韧性提升方法

Similar Documents

Publication Publication Date Title
Romero Blackouts illuminate India's power problems
Kiprakis et al. Maximising energy capture from distributed generators in weak networks
CN105977934B (zh) 一种考虑负荷重要性的配电网智能软开关供电恢复方法
JP7296657B2 (ja) マルチエネルギー協同の配電網用アクティブ解列制御方法
Mohamed et al. Microgrid online management and balancing using multiobjective optimization
CN104934971A (zh) 基于潮流转移比的动态断面控制方法
CN112260261B (zh) 一种电-气综合能源系统脆弱性评估与改善方法
CN113659572B (zh) 考虑网络重构和需求响应的气电综合能源配网鲁棒优化方法
CN103199512B (zh) 配电网单相接地故障网络重构方法
WO2021203480A1 (zh) 基于可靠性约束的配电网自动化系统综合规划方法
CN112529277B (zh) 一种基于资源调配的热电耦合系统灾前预防方法
CN116128110A (zh) 一种电热综合能源系统的多阶段韧性提升方法及终端
CN103997035A (zh) 电力变压器及线路经济运行闭环控制系统及方法
CN113852084B (zh) 考虑智能软开关优化配置的配电网多阶段韧性提升方法
CN113991653A (zh) 一种考虑负荷重要程度的三阶段拓扑运行方式优化方法
CN116934021A (zh) 考虑n-1安全约束的电-热综合能源系统优化调度方法
CN104135028A (zh) 一种结合黑启动优化模型与fcb工况的机组布点方法
CN116169673A (zh) 配电网的弹性评估方法、系统、电子设备及可读存储介质
CN113486509B (zh) 一种综合能源系统多目标优化控制方法
CN110556824B (zh) 一种基于母线分裂的输电能力提升方法
CN103822400A (zh) 一种高设备利用率模块型冷热电三联供系统
CN117996739A (zh) 含智能软开关的配电系统与区域供热系统的负荷恢复方法
CN112528480B (zh) 一种提升负荷裕度的热电耦合系统灾前预防方法
CN117910769A (zh) 一种考虑电-热综合能源系统重构的韧性调度方法及系统
CN115396294B (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