提升灵活性的电热综合能源系统两阶段随机优化调度方法与流程

1.本发明涉及电热综合能源系统优化运行调度方法技术领域，是一种提升灵活性的电热综合能源系统两阶段随机优化调度方法。


背景技术：

2.风电作为清洁和可再生能源对于缓解全球能源短缺和环境污染问题具有重要的作用，发展前景广阔。近年来我国风电产业发展迅速，装机容量位居世界第一，但弃风问题突出。在“三北”地区供热季，热电联产机组工作在“以热定电”运行模式，且具有较强的“热电耦合”特性，灵活调节能力受限，并且风电出力具有明显的随机性和波动性，使得电力系统面临严峻的调峰和风电消纳问题。传统的能源系统仅考虑电、热等独立单一的能源形式，割裂了不同能源系统的联系，难以发挥互补优势。因此，考虑风电的不确定性，构建电热综合能源系统协调优化调度模型对提升系统风电消纳能力具有重要意义。
3.一方面，当前关于电热综合能源系统协调优化的研究多侧重于储热、电锅炉、热泵等供给侧，或热网、需求侧的热动态特性，但计及风电不确定性，综合考虑供给侧、网侧和需求侧的协调优化的研究尚不多见。
4.另一方面，为将风电强随机性和波动性等特性融入调度决策中，现有研究多采用随机优化和鲁棒优化理论构建优化调度模型。两阶段随机优化在决策过程中充分考虑了风电不确定性，但鲜少应用于电热综合能源系统优化调度中，并且如何综合考虑热电联产机组的热-电解耦运行特性和提供备用的潜力也有待深入研究。


技术实现要素：

5.本发明提供了一种提升灵活性的电热综合能源系统两阶段随机优化调度方法，克服了上述现有技术之不足，其综合考虑了电力系统和供热系统的灵活性，利用热力系统的动态运行特性和用能灵活特性，实现热电解耦，在保证供热质量的前提下提升电力系统的运行灵活性，拓宽风电的消纳区间，实现大规模风电并网消纳。
6.本发明的技术方案是通过以下措施来实现的：一种提升灵活性的电热综合能源系统两阶段随机优化调度方法，包括以下步骤：
7.步骤1，电力系统和热力系统的运行参数和历史运行数据的收集梳理；
8.步骤2，构建计及储热能力的热网动态特性模型；
9.步骤3，生成风电动态随机场景表征风电出力的随机特性；
10.步骤4，构建提升灵活性的电热综合能源系统两阶段随机优化调度模型；
11.步骤5，求解所述电热综合能源系统两阶段随机优化调度模型，得到电热综合能源系统调度计划。
12.下面是对上述发明技术方案的进一步优化或/和改进：
13.上述步骤2中，供热系统(热网)由热电联产机组、储热装置等热源、供热管网、热交换站和热负荷等组成。热源产生的热量通过循环热水供应给用户，以维持室内温度，其中一
次管网和二次管网通过热交换站进行热量传递。一方面，热力系统传输较慢，具有较大的惯性，其延时效应可等效为系统储热能力，而供热区域具有一定的热惯性，也能提供一定的储热能力，两者结合可在不增加额外投资的情况下提升系统运行灵活性。由于二次管网相对较短、储能容量较小，因此针对一次管网进行建模分析。另一方面，在热电联产机组配置储热罐可在电负荷高峰期间储热，在电负荷低谷期间放热，同时建筑物对温度要求具有一定的弹性，具有一定的储热能力，综合作用下能够实现热电联产机组发电功率的灵活调节，拓宽风电的出力空间。热网采用广泛应用和研究的质调节运行方式，即保持热网中的热水质流速率恒定，仅通过调节热源处的供水温度实现供热功率的调整。
14.所述热网动态特性模型包括热源供热功率约束、热网动态运行约束和热负荷交换站运行约束；
15.步骤2.1，配置储热罐的热电联产机组作为热源共同向热网供热，热源供热功率约束(总供热功率约束)如式(1)所示，
[0016][0017]
式中，h
chp,h,t
、h
hs,s,t
分别为热电机组h及其配置的储热罐s在时段t的热功率；为热源与热网进行热量交换的热源换热站j在时段t的热水质流速率；和分别为热源交换站所在的管道节点n的供水管道入口和回水管道出口在时段t的热水温度；c
water
为水的比热容；i
chp
、i
hs
、t分别为热电联产机组、热源交换站、运行时段的集合；为与热源交换站相连的管道节点集合；
[0018]
步骤2.2，热网动态运行约束包括节点质流速率连续性约束、节点温度混合约束、管道供水及回水温度约束和管道温度动态及热损耗特性；
[0019]
步骤2.2.1，每个节点热水流入质流速率等于流出质流速率，供水和回水网络的节点质流速率连续约束如式(2)和式(3)所示，
[0020][0021][0022]
式中，分别为管道b在时段t的热水质流速率；分别为以管道节点n为起点和终点的管道集合；i
nd
为热网节点集合；
[0023]
步骤2.2.2，不同温度的热水从不同管道流经同一节点进行混合，混合后从该节点流入不同管道的热水温度保持相同，则节点温度混合约束(混合媒介的温度约束)如式(4)至式(7)所示，
[0024][0025][0026]
[0027][0028]
式中，分别为供水管道和回水管道在时段t的出口处的热水温度；分别为供水管道和回水管道在时段t的入口处的热水温度；分别为供水网络节点和回水网络节点的热水混合温度；
[0029]
步骤2.2.3，管道供水及回水温度约束如式(8)至式(11)所示，
[0030][0031][0032][0033][0034]
式中，为供水管道的热水温度的上、下限；为回水管道的热水温度的上、下限；s
pipe,s
、s
pipe,r
为供水管道集合、回水管道集合；
[0035]
步骤2.2.4，管道温度动态及热损耗特性为单个管道的热动态特性体现在管道供热媒介的温度传输时延和热损耗，其出口和入口热水温度的动态关系如(12)至式(14)所示，
[0036][0037][0038][0039]
式中，为管道周围的土壤温度；t
phr,b
为管道的热储备系数；δt为调度时间间隔；γb为管道的供热延迟时间，tb是处理后的供热延迟时间，等于γb按照δt离散化之后最近的整数，便于进行计算；ρ
water
为水的密度；lb、rb、μb、分别为管道长度、半径、热损失系数、热水的质流速率；
[0040]
步骤2.3，热负荷交换站从热网获得热量用于满足用户的热负荷需求，则热负荷交换站运行约束如式(15)所示，
[0041][0042]
式中，为热负荷交换站从热网获得的热功率；为热负荷交换站向用户提供的热功率；为热负荷交换站k在时段t的热水质流速率；i
hes
、分别为热负荷交换
站集合、与热负荷交换站相连的节点集合。
[0043]
上述步骤3中，为了准确刻画风电的随机性和波动性，采用计及随机变量相关性的动态场景法生成风电出力场景集，风电动态随机场景生成的步骤如下：
[0044]
步骤3.1，对风电功率预测值、实测值、预测误差组成的历史数据进行标幺化处理，将历史数据按照风电预测值划分为n个预测箱，每个预测箱含有n
l
个数据组，各预测箱误差样本从到大排序为，拟合得到预测误差e的累计经验分布函数，如式(16)所示，进而基于日前各时段的风电预测值确定各时段所在的预测箱的累积概率分布函数；
[0045][0046]
步骤3.2，生成ns个服从多元联合标准正态分布的随机向量z＝(z1,z2...z
t
)
′
，z～n(μ0,σ)，μ0为t维零向量，协方差矩阵σ如式(17)所示，
[0047][0048][0049]
式中，为变量zm和zn的协方差；t为预测时段数；ε为调节参数，用于控制不同时段随机变量间的相关性强度；
[0050]
步骤3.3，通过等概率逆变换，将ns个服从多元标准正态分布的随机向量转化为具有相关性的风电功率误差场景值e，进而将不同时段的误差场景值和预测值相加后并将其转换为有名值，得到ns个风电动态随机场景，
[0051][0052]
式中：φ为标准正态分布的累积概率分布函数。
[0053]
上述步骤4中，所述电热综合能源系统两阶段随机优化调度模型以第一阶段的日前调度成本f
da
和第二阶段的随机场景下的实时调整成本期望f
rt
之和f最小为目标函数，如式(20)所示；第一部分为日前阶段的调度成本f
da
，包括常规火电机组运行成本(开停机成本、发电成本、备用容量成本)，热电联产机组运行成本(发电成本、备用容量成本)，弃风惩罚成本，如式(21)所示；第二部分为实时阶段的调度成本f
rt
，包括常规火电机组备用调用成本期望值、热电联产机组备用调用成本期望值、非自愿切负荷和弃风成本期望值，如式(22)所示；
[0054]
minf＝min(f
da
+f
rt
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(20)
[0055][0056][0057]
式中，u
i,t
、y
i,t
、z
i,t
分别为常规火电机组i的运行状态变量、启动变量、停机变量；s
u,i
和s
d,i
分别为常规火电机组i单次启动成本和停机成本；a
1,i
、a
2,i
、a
3,i
为常规火电机组i的发电成本系数；p
g,i,t
为常规火电机组i的发电功率；和分别为常规火电机组i提供的向上备用和价格系数；和分别为常规火电机组i提供的向下备用和价格系数；p
chp,h,t
和h
chp,h,t
分别为热电联产机组h的发电功率和热功率；b
1,h
、b
2,h
、b
3,h
、b
4,h
、b
5,h
、b
6,h
为热电联产机组的运行成本系数；和分别为热电联产机组h提供的向上备用和价格系数；和分别为热电联产机组h提供的向下备用和价格系数；和p
w,t
为风电场w的预测功率和调度功率；c
wind
为弃风惩罚成本系数；和分别为常规火电机组i在场景s中调用的向上备用和调用能量价格系数；和分别为常规火电机组i在场景s中调用的向下备用和调用能量价格系数；和分别为热电联产机组h在场景s中调用的向上备用和调用能量价格系数；和分别为热电联产机组h在场景s中调用的向下备用和调用价格系数；和c
load
为场景s下的切负荷功率和切负荷惩罚成本系数；和p
w,t,s
为风电场w在场景s的实际功率和调度功率；n
chp
、ng、nw、n
load
、ns分别为热电联产机组、常规火电机组、风电场、用电负荷、随机场景的数量；ps为场景s的概率。
[0058]
上述电热综合能源系统两阶段随机优化调度模型第一阶段日前调度约束条件是指预测场景下的日前调度约束条件，包括耦合元件运行约束、电网运行约束和热网运行约束，
[0059]
(1)耦合元件运行约束即热电联产机组运行约束，热电联产机组运行约束包括电出力和热出力运行约束、热电联产机组爬坡约束以及热电机组调节灵活性备用约束；
[0060]
抽气式热电联产机组的运行域是一个多面体，其电出力和热出力运行约束如式(23)所示，
[0061][0062]
式中：分别为热电联产机组h纯凝工况下的最大、最小的发电功率；c
v,h
为热电联产机组h总进气量恒定时抽取单位供热量下发电功率的减小量；c
m,h
为热电联产机组h背压状态运行下的电功率和热功率的弹性系数；为热电联产机组最大热出力；c
k,h
为常数项；
[0063]
热电联产机组爬坡约束如式(24)所示，
[0064][0065]
式中，和分别为热电联产机组h发电功率的向上爬坡能力和向下爬坡能力；
[0066]
热电机组调节灵活性备用约束如式(25)至(26)所示，
[0067][0068][0069]
(2)电网运行约束包括电功率平衡约束、线路潮流约束和常规火电机组运行约束、常规火电机组爬坡约束、常规火电机组备用容量约束、常规火电机组最小运行时间约束、常规火电机组最小停机时间约束、常规火电机组开停机状态约束、风电场出力约束；
[0070]
1)电功率平衡约束如式(27)所示，
[0071][0072]
式中，p
d,t
为负荷d的用电功率；
[0073]
2)线路潮流约束如式(28)所示，
[0074][0075]
式中，λ
h,l
、λ
i,l
、λ
w,l
、λ
d,l
分别为热电联产机组h、常规火电机组i、风电场w、负荷d对网架输电线路l的潮流分布转移因子；f
l,max
为输电线路l的功率输送限值；
[0076]
3)常规火电机组运行约束如式(29)所示，
[0077]
p
g,i,min
≤p
g,i,t
≤p
g,i,max
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(29)
[0078]
式中，p
i,min
和p
i,max
分别为常规火电机组i的最小有功出力和最大有功出力；
[0079]
4)常规火电机组爬坡约束如式(30)所示，
[0080]
[0081]
式中，和分别为常规火电机组i的向上爬坡能力和向下爬坡能力；
[0082]
5)常规火电机组备用容量约束如式(31)所示，
[0083][0084]
6)常规火电机组最小运行时间约束如式(32)所示，
[0085][0086]
式中，t
u,i
和u
i,0
分别为常规机组i的最小启动运行时间和初始运行时间；
[0087]
7)常规火电机组最小停机时间约束如式(33)所示，
[0088][0089]
8)常规火电机组开停机状态约束如式(34)所示，
[0090]yi,t-z
i,t
＝u
i,t-u
i,t-1
，y
i,t
+z
i,t
≤1
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(34)
[0091]
9)风电场出力约束如式(35)所示，
[0092][0093]
(3)热网运行约束包括热源供热功率约束、热网动态运行约束、储热罐的热功率约束和建筑物热动态约束，热源供热功率约束如式(1)所示，热网动态运行约束如式(2)至(14)所示，
[0094]
在热电联产机组处配置储热罐提升其运行灵活性，储热罐的热功率约束如式(36)所示，
[0095]
[0096]
式中：和为分别为第s个储热罐在t时段的储热状态和放热状态；h
hs,s,t
、分别为第s个储热罐在t时段的热功率、储热功率和放热功率；分别为第s个储热罐在t时段的热功率、储热功率和放热功率；分别为第s个储热罐的最小储热功率、最大储热功率、最小放热功率、最大放热功率；e
hs,s,t
为第s个储热罐在t时段的储热量；e
hs,s,min
、e
hs,s,max
、e
hs,s,0
为第s个储热罐的最小储热量、最大储热量、初始储热量；
[0097]
考虑到建筑物对室内温度要求具有一定的弹性，通过优化调节可以提供一定的储热能力，建筑物热动态约束包括室内温度-建筑物的热量-调度时段耦合关系约束和室内温度范围约束；
[0098]
室内温度-建筑物的热量-调度时段耦合关系约束如式(37)所示，
[0099][0100]
式中，τ
od,k,t
、τ
id,k,t
分别为建筑物k的室外温度、室内温度；χ
bt,k
、t
bhr,k
分别为建筑物k的单位温差传热系数、热储备系数；
[0101]
室内温度范围约束如式(38)所示，为了保证热负荷用户的用热舒适度，建筑物的室内温度要保持在合理的范围内。
[0102][0103]
式中，分别不影响用户的用热舒适度下的最大室内温度、最小室内温度。
[0104]
上述电热综合能源系统两阶段随机优化调度模型第二阶段约束条件基于第一阶段确定的日前调度计划，第二阶段为应对实时风电随机场景的不确定性，对日前调度计划进行优化调整，最小化实时调整的期望成本，并且在日内实时调度阶段时，不再调整热网中热源的供热计划，仅对供电计划进行优化调整，优化调整调整过程中，场景s下的电功率平衡约束、常规火电机组爬坡约束、热电联产机组爬坡约束、风电功率、切负荷约束、备用容量约束、线路潮流约束如下：
[0105]
1)场景s下的电功率平衡约束如式(39)所示，
[0106][0107]
2)场景s下的常规火电机组爬坡约束如式(40)至(41)所示，
[0108][0109][0110]
3)场景s下的热电联产机组爬坡约束如式(42)至(43)所示，
[0111][0112]
[0113]
4)场景s下的风电功率、切负荷约束如式(44)所示，
[0114][0115]
5)备用容量约束如式(45)至(46)所示，
[0116][0117][0118]
线路潮流约束如式(47)所示，
[0119][0120]
上述步骤5中，基于所述电热综合能源系统两阶段随机优化调度模型，采用yalmip工具箱结合商业化优化软件cplex进行求解，得到日前调度计划中电-热综合能源系统中常规火电机组的运行状态u
i,t
、发电功率p
g,i,t
、提供的向上备用容量、向下备用容量热电联产机组的发电功率p
chp,h,t
、供热功率h
chp,h,t
、提供的向上备用容量、向下备用容量，储热罐的储热功率、放热功率，风电场的出力计划p w,t
，用户建筑物内的温度τ
id,k,t
，实现电-热综合能源系统的协调优化，提升系统运行的灵活性。
[0121]
本发明综合考虑了电力系统和供热系统的灵活性，利用热力系统的动态运行特性和用能灵活特性，实现热电解耦，在保证供热质量的前提下提升电力系统的运行灵活性，拓宽风电的消纳区间，实现大规模风电并网消纳。
[0122]
在电力系统方面，所构建的调度模型能够发挥常规火电机组的运行灵活性和提供备用的能力，热电联产机组通过与热力系统协调实现热电解耦，能够拓宽出力区间并发挥提供备用的能力。
[0123]
在热力系统方面，综合考虑储热罐运行约束、热网的动态运行特性、建筑物热动态特性，挖掘供给侧，热网侧、热需求侧的综合储热特性和灵活运行特性，实现供电和供热的有效解耦，进而提升电力系统运行灵活性。基于风电动态随机场景，通过构建两阶段随机优化模型，最优化预测场景运行成本和随机场景的期望调节成本，能够兼顾电热综合能源系统的运行灵活性和调节灵活性。
[0124]
本发明方法的优点在于考虑热力系统中储热罐运行模型、热网的动态运行特性、建筑物热动态特性，挖掘供给侧、热网侧、热需求侧的综合储热特性和灵活运行特性进而实现电-热供应解耦，同时考虑电力系统中常规火电机组和热电联产机组的运行灵活性和提供备用的能力，挖掘热、电系统的潜在灵活性，能够从整个综合能源系统范围内实现灵活性资源的协调优化，提升整体运行的灵活性，为调度运行人员提供更合理的供热和供电调度计划。
附图说明
[0125]
附图1为本发明的流程图。
具体实施方式
[0126]
本发明不受下述实施例的限制，可根据本发明的技术方案与实际情况来确定具体的实施方式。
[0127]
本发明中，如无特别说明，使用的工具、软件均为现有公知公用的工具、软件。
[0128]
下面结合实施例对本发明作进一步描述：
[0129]
实施例：如附图1所示，该提升灵活性的电热综合能源系统两阶段随机优化调度方法，包括以下步骤：
[0130]
步骤1，电力系统和热力系统的运行参数和历史运行数据的收集梳理；
[0131]
步骤2，构建计及储热能力的热网动态特性模型；
[0132]
步骤3，生成风电动态随机场景表征风电出力的随机特性；
[0133]
步骤4，构建提升灵活性的电热综合能源系统两阶段随机优化调度模型；
[0134]
步骤5，求解所述电热综合能源系统两阶段随机优化调度模型，得到电热综合能源系统调度计划。提升灵活性的电热综合能源系统两阶段随机优化调度模型中，第一阶段基于预测信息，满足电-热系统的运行约束，以系统运行成本最小为目标优化得到日前调度计划，第二阶段在第一阶段结果的基础上进行实时不确定场景下的再优化，即以风电随机场景下的备用调用成本、弃风成本和切负荷成本最小为目标进行调整优化，保证优化结果能够满足实时运行要求。
[0135]
步骤2中，所述热网动态特性模型包括热源供热功率约束、热网动态运行约束和热负荷交换站运行约束；
[0136]
步骤2.1，配置储热罐的热电联产机组作为热源共同向热网供热，热源供热功率约束(总供热功率约束)如式(1)所示，
[0137][0138]
式中，h
chp,h,t
、h
hs,s,t
分别为热电机组h及其配置的储热罐s在时段t的热功率；为热源与热网进行热量交换的热源换热站j在时段t的热水质流速率；和分别为热源交换站所在的管道节点n的供水管道入口和回水管道出口在时段t的热水温度；c
water
为水的比热容；i
chp
、i
hs
、t分别为热电联产机组、热源交换站、运行时段的集合；为与热源交换站相连的管道节点集合；
[0139]
步骤2.2，热网动态运行约束包括节点质流速率连续性约束、节点温度混合约束、管道供水及回水温度约束和管道温度动态及热损耗特性；
[0140]
步骤2.2.1，每个节点热水流入质流速率等于流出质流速率，供水和回水网络的节点质流速率连续约束如式(2)和式(3)所示，
[0141][0142][0143]
式中，分别为管道b在时段t的热水质流速率；分别为以管道节点n为起点和终点的管道集合；i
nd
为热网节点集合；
[0144]
步骤2.2.2，不同温度的热水从不同管道流经同一节点进行混合，混合后从该节点流入不同管道的热水温度保持相同，则节点温度混合约束(混合媒介的温度约束)如式(4)至式(7)所示，
[0145][0146][0147][0148][0149]
式中，分别为供水管道和回水管道在时段t的出口处的热水温度；分别为供水管道和回水管道在时段t的入口处的热水温度；分别为供水网络节点和回水网络节点的热水混合温度；
[0150]
步骤2.2.3，管道供水及回水温度约束如式(8)至式(11)所示，
[0151][0152][0153][0154][0155]
式中，为供水管道的热水温度的上、下限；为回水管道的热水温度的上、下限；s
pipe,s
、s
pipe,r
为供水管道集合、回水管道集合；
[0156]
步骤2.2.4，管道温度动态及热损耗特性为单个管道的热动态特性体现在管道供热媒介的温度传输时延和热损耗，其出口和入口热水温度的动态关系如(12)至式(14)所示，
[0157][0158][0159][0160]
式中，为管道周围的土壤温度；t
phr,b
为管道的热储备系数；δt为调度时间间隔；γb为管道的供热延迟时间，tb是处理后的供热延迟时间，等于γb按照δt离散化之后最近的整数，便于进行计算；ρ
water
为水的密度；lb、rb、μb、分别为管道长度、半径、热损失系数、热水的质流速率；
[0161]
步骤2.3，热负荷交换站从热网获得热量用于满足用户的热负荷需求，则热负荷交换站运行约束如式(15)所示，
[0162][0163]
式中，为热负荷交换站从热网获得的热功率；为热负荷交换站向用户提供的热功率；为热负荷交换站k在时段t的热水质流速率；i
hes
、分别为热负荷交换站集合、与热负荷交换站相连的节点集合。
[0164]
上述步骤3中，为了准确刻画风电的随机性和波动性，采用计及随机变量相关性的动态场景法生成风电出力场景集，风电动态随机场景生成的步骤如下：
[0165]
步骤3.1，对风电功率预测值、实测值、预测误差组成的历史数据进行标幺化处理，将历史数据按照风电预测值划分为n个预测箱，每个预测箱含有n
l
个数据组，各预测箱误差样本从到大排序为拟合得到预测误差e的累计经验分布函数，如式(16)所示，进而基于日前各时段的风电预测值确定各时段所在的预测箱的累积概率分布函数；
[0166][0167]
步骤3.2，生成ns个服从多元联合标准正态分布的随机向量z＝(z1,z2...z
t
)
′
，z～n(μ0,σ)，μ0为t维零向量，协方差矩阵σ如式(17)所示，
[0168][0169][0170]
式中，为变量zm和zn的协方差；t为预测时段数；ε为调节参数，用于控制不同时段随机变量间的相关性强度；
[0171]
步骤3.3，通过等概率逆变换，将ns个服从多元标准正态分布的随机向量转化为具有相关性的风电功率误差场景值e，进而将不同时段的误差场景值和预测值相加后并将其转换为有名值，得到ns个风电动态随机场景，
[0172][0173]
式中：φ为标准正态分布的累积概率分布函数。
[0174]
上述步骤4中，所述电热综合能源系统两阶段随机优化调度模型以第一阶段的日前调度成本f
da
和第二阶段的随机场景下的实时调整成本期望f
rt
之和f最小为目标函数，如式(20)所示；第一部分为日前阶段的调度成本f
da
，包括常规火电机组运行成本(开停机成本、发电成本、备用容量成本)，热电联产机组运行成本(发电成本、备用容量成本)，
弃风惩罚成本，如式(21)所示；第二部分为实时阶段的调度成本f
rt
，包括常规火电机组备用调用成本期望值、热电联产机组备用调用成本期望值、非自愿切负荷和弃风成本期望值，如式(22)所示；
[0175]
minf＝min(f
da
+f
rt
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(20)
[0176][0177][0178]
式中，u
i,t
、y
i,t
、z
i,t
分别为常规火电机组i的运行状态变量、启动变量、停机变量；s
u,i
和s
d,i
分别为常规火电机组i单次启动成本和停机成本；a
1,i
、a
2,i
、a
3,i
为常规火电机组i的发电成本系数；p
g,i,t
为常规火电机组i的发电功率；和分别为常规火电机组i提供的向上备用和价格系数；和分别为常规火电机组i提供的向下备用和价格系数；p
chp,h,t
和h
chp,h,t
分别为热电联产机组h的发电功率和热功率；b
1,h
、b
2,h
、b
3,h
、b
4,h
、b
5,h
、b
6,h
为热电联产机组的运行成本系数；和分别为热电联产机组h提供的向上备用和价格系数；和分别为热电联产机组h提供的向下备用和价格系数；和p
w,t
为风电场w的预测功率和调度功率；c
wind
为弃风惩罚成本系数；和分别为常规火电机组i在场景s中调用的向上备用和调用能量价格系数；和分别为常规火电机组i在场景s中调用的向下备用和调用能量价格系数；和分别为热电联产机组h在场景s中调用的向上备用和调用能量价格系数；和分别为热电联产机组h在场景s中调用的向下备用和调用价格系数；和为场景s下的切负荷功率和切负荷惩罚成本系数；和p
w,t,s
为风电场w在场景s的实际功率和调度功率；n
chp
、ng、nw、n
load
、ns分别为热电联产机组、常规火电机组、风电场、用电负荷、随机场景的数量；ps为场景s的概率；
[0179]
所述电热综合能源系统两阶段随机优化调度模型第一阶段日前调度约束条件是指预测场景下的日前调度约束条件，包括耦合元件运行约束、电网运行约束和热网运行约束，
[0180]
(1)耦合元件运行约束即热电联产机组运行约束，热电联产机组运行约束包括电出力和热出力运行约束、热电联产机组爬坡约束以及热电机组调节灵活性备用约束；
[0181]
抽气式热电联产机组的运行域是一个多面体，其电出力和热出力运行约束如式(23)所示，
[0182][0183]
式中：分别为热电联产机组h纯凝工况下的最大、最小的发电功率；c
v,h
为热电联产机组h总进气量恒定时抽取单位供热量下发电功率的减小量；c
m,h
为热电联产机组h背压状态运行下的电功率和热功率的弹性系数；为热电联产机组最大热出力；c
k,h
为常数项；
[0184]
热电联产机组爬坡约束如式(24)所示，
[0185][0186]
式中，和分别为热电联产机组h发电功率的向上爬坡能力和向下爬坡能力；
[0187]
热电机组调节灵活性备用约束如式(25)至(26)所示，
[0188][0189][0190]
(2)电网运行约束包括电功率平衡约束、线路潮流约束和常规火电机组运行约束、常规火电机组爬坡约束、常规火电机组备用容量约束、常规火电机组最小运行时间约束、常规火电机组最小停机时间约束、常规火电机组开停机状态约束、风电场出力约束；
[0191]
1)电功率平衡约束如式(27)所示，
[0192][0193]
式中，p
d,t
为负荷d的用电功率；
[0194]
2)线路潮流约束如式(28)所示，
[0195][0196]
式中，λ
h,l
、λ
i,l
、λ
w,l
、λ
d,l
分别为热电联产机组h、常规火电机组i、风电场w、负荷d对网架输电线路l的潮流分布转移因子；f
l,max
为输电线路l的功率输送限值；
[0197]
3)常规火电机组运行约束如式(29)所示，
[0198]
p
g,i,min
≤p
g,i,t
≤p
g,i,max
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(29)
[0199]
式中，p
i,min
和p
i,max
分别为常规火电机组i的最小有功出力和最大有功出力；
[0200]
4)常规火电机组爬坡约束如式(30)所示，
[0201][0202]
式中，和分别为常规火电机组i的向上爬坡能力和向下爬坡能力；
[0203]
5)常规火电机组备用容量约束如式(31)所示，
[0204][0205]
6)常规火电机组最小运行时间约束如式(32)所示，
[0206][0207]
式中，t
u,i
和u
i,0
分别为常规机组i的最小启动运行时间和初始运行时间；
[0208]
7)常规火电机组最小停机时间约束如式(33)所示，
[0209][0210]
8)常规火电机组开停机状态约束如式(34)所示，
[0211]yi,t-z
i,t
＝u
i,t-u
i,t-1
，y
i,t
+z
i,t
≤1
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(34)
[0212]
9)风电场出力约束如式(35)所示，
[0213][0214]
(3)热网运行约束包括热源供热功率约束、热网动态运行约束、储热罐的热功率约束和建筑物热动态约束，热源供热功率约束如式(1)所示，热网动态运行约束如式(2)至(14)所示，
[0215]
在热电联产机组处配置储热罐提升其运行灵活性，储热罐的热功率约束如式(36)所示，
[0216][0217]
式中：和为分别为第s个储热罐在t时段的储热状态和放热状态；h
hs,s,t
、分别为第s个储热罐在t时段的热功率、储热功率和放热功率；分别为第s个储热罐在t时段的热功率、储热功率和放热功率；分别为第s个储热罐的最小储热功率、最大储热功率、最小放热功率、最大放热功率；e
hs,s,t
为第s个储热罐在t时段的储热量；e
hs,s,min
、e
hs,s,max
、e
hs,s,0
为第s个储热罐的最小储热量、最大储热量、初始储热量；
[0218]
考虑到建筑物对室内温度要求具有一定的弹性，通过优化调节可以提供一定的储热能力，建筑物热动态约束包括室内温度-建筑物的热量-调度时段耦合关系约束和室内温度范围约束；
[0219]
室内温度-建筑物的热量-调度时段耦合关系约束如式(37)所示，
[0220][0221]
式中，τ
od,k,t
、τ
id,k,t
分别为建筑物k的室外温度、室内温度；χ
bt,k
、t
bhr,k
分别为建筑物k的单位温差传热系数、热储备系数；
[0222]
室内温度范围约束如式(38)所示，
[0223][0224]
式中，分别不影响用户的用热舒适度下的最大室内温度、最小室内温度。
[0225]
所述电热综合能源系统两阶段随机优化调度模型第二阶段约束条件基于第一阶段确定的日前调度计划，第二阶段为应对实时风电随机场景的不确定性，对日前调度计划进行优化调整，最小化实时调整的期望成本，并且在日内实时调度阶段时，不再调整热网中热源的供热计划，仅对供电计划进行优化调整，优化调整调整过程中，场景s下的电功率平衡约束、常规火电机组爬坡约束、热电联产机组爬坡约束、风电功率、切负荷约束、备用容量约束、线路潮流约束如下：
[0226]
1)场景s下的电功率平衡约束如式(39)所示，
[0227]
[0228]
2)场景s下的常规火电机组爬坡约束如式(40)至(41)所示，
[0229][0230][0231]
3)场景s下的热电联产机组爬坡约束如式(42)至(43)所示，
[0232][0233][0234]
4)场景s下的风电功率、切负荷约束如式(44)所示，
[0235][0236]
5)备用容量约束如式(45)至(46)所示，
[0237][0238][0239]
线路潮流约束如式(47)所示，
[0240][0241]
所述步骤5中，基于所述电热综合能源系统两阶段随机优化调度模型，采用yalmip工具箱结合商业化优化软件cplex进行求解，得到日前调度计划中电-热综合能源系统中常规火电机组的运行状态u
i,t
、发电功率p
g,i,t
、提供的向上备用容量、向下备用容量热电联产机组的发电功率p
chp,h,t
、供热功率h
chp,h,t
、提供的向上备用容量向下备用容量，储热罐的储热功率、放热功率，风电场的出力计划p
w,t
，用户建筑物内的温度τ
id,k,t
，实现电-热综合能源系统的协调优化，提升系统运行的灵活性。
[0242]
以上技术特征构成了本发明的实施例，其具有较强的适应性和实施效果，可根据实际需要增减非必要的技术特征，来满足不同情况的需求。

技术特征：
1.一种提升灵活性的电热综合能源系统两阶段随机优化调度方法，其特征在于包括以下步骤：步骤1，电力系统和热力系统的运行参数和历史运行数据的收集梳理；步骤2，构建计及储热能力的热网动态特性模型；步骤3，生成风电动态随机场景表征风电出力的随机特性；步骤4，构建提升灵活性的电热综合能源系统两阶段随机优化调度模型；步骤5，求解所述电热综合能源系统两阶段随机优化调度模型，得到电热综合能源系统调度计划。2.根据权利要求1所述的提升灵活性的电热综合能源系统两阶段随机优化调度方法，其特征在于步骤2中，所述热网动态特性模型包括热源供热功率约束、热网动态运行约束和热负荷交换站运行约束；步骤2.1，配置储热罐的热电联产机组作为热源共同向热网供热，热源供热功率约束如式(1)所示，式中，h
chp,h,t
、h
hs,s,t
分别为热电机组h及其配置的储热罐s在时段t的热功率；为热源与热网进行热量交换的热源换热站j在时段t的热水质流速率；和分别为热源交换站所在的管道节点n的供水管道入口和回水管道出口在时段t的热水温度；c
water
为水的比热容；i
chp
、i
hs
、t分别为热电联产机组、热源交换站、运行时段的集合；为与热源交换站相连的管道节点集合；步骤2.2，热网动态运行约束包括节点质流速率连续性约束、节点温度混合约束、管道供水及回水温度约束和管道温度动态及热损耗特性；步骤2.2.1，每个节点热水流入质流速率等于流出质流速率，供水和回水网络的节点质流速率连续约束如式(2)和式(3)所示，流速率连续约束如式(2)和式(3)所示，式中，分别为管道b在时段t的热水质流速率；分别为以管道节点n为起点和终点的管道集合；i
nd
为热网节点集合；步骤2.2.2，不同温度的热水从不同管道流经同一节点进行混合，混合后从该节点流入不同管道的热水温度保持相同，则节点温度混合约束如式(4)至式(7)所示，不同管道的热水温度保持相同，则节点温度混合约束如式(4)至式(7)所示，不同管道的热水温度保持相同，则节点温度混合约束如式(4)至式(7)所示，
式中，分别为供水管道和回水管道在时段t的出口处的热水温度；分别为供水管道和回水管道在时段t的入口处的热水温度；分别为供水网络节点和回水网络节点的热水混合温度；步骤2.2.3，管道供水及回水温度约束如式(8)至式(11)所示，步骤2.2.3，管道供水及回水温度约束如式(8)至式(11)所示，步骤2.2.3，管道供水及回水温度约束如式(8)至式(11)所示，步骤2.2.3，管道供水及回水温度约束如式(8)至式(11)所示，式中，为供水管道的热水温度的上、下限；为回水管道的热水温度的上、下限；s
pipe,s
、s
pipe,r
为供水管道集合、回水管道集合；步骤2.2.4，管道温度动态及热损耗特性为单个管道的热动态特性体现在管道供热媒介的温度传输时延和热损耗，其出口和入口热水温度的动态关系如(12)至式(14)所示，水温度的动态关系如(12)至式(14)所示，水温度的动态关系如(12)至式(14)所示，式中，为管道周围的土壤温度；t
phr,b
为管道的热储备系数；δt为调度时间间隔；γ
b
为管道的供热延迟时间，t
b
是处理后的供热延迟时间，等于γ
b
按照δt离散化之后最近的整数，便于进行计算；ρ
water
为水的密度；l
b
、r
b
、μ
b
、分别为管道长度、半径、热损失系数、热水的质流速率；步骤2.3，热负荷交换站从热网获得热量用于满足用户的热负荷需求，则热负荷交换站运行约束如式(15)所示，式中，为热负荷交换站从热网获得的热功率；为热负荷交换站向用户提供的热功率；为热负荷交换站k在时段t的热水质流速率；i
hes
、分别为热负荷交换站集合、与热负荷交换站相连的节点集合。3.根据权利要求1或2所述的提升灵活性的电热综合能源系统两阶段随机优化调度方
法，其特征在于步骤3中，风电动态随机场景生成的步骤如下：步骤3.1，对风电功率预测值、实测值、预测误差组成的历史数据进行标幺化处理，将历史数据按照风电预测值划分为n个预测箱，每个预测箱含有n
l
个数据组，各预测箱误差样本从到大排序为拟合得到预测误差e的累计经验分布函数如式(16)所示，进而基于日前各时段的风电预测值确定各时段所在的预测箱的累积概率分布函数；步骤3.2，生成n
s
个服从多元联合标准正态分布的随机向量z＝(z1,z2...z
t
)
′
，z～n(μ0,σ)，μ0为t维零向量，协方差矩阵σ如式(17)所示，为t维零向量，协方差矩阵σ如式(17)所示，式中，为变量z
m
和z
n
的协方差；t为预测时段数；ε为调节参数，用于控制不同时段随机变量间的相关性强度；步骤3.3，通过等概率逆变换，将n
s
个服从多元标准正态分布的随机向量转化为具有相关性的风电功率误差场景值e，进而将不同时段的误差场景值和预测值相加后并将其转换为有名值，得到n
s
个风电动态随机场景，式中：φ为标准正态分布的累积概率分布函数。4.根据权利要求3所述的提升灵活性的电热综合能源系统两阶段随机优化调度方法，其特征在于步骤4中，所述电热综合能源系统两阶段随机优化调度模型以第一阶段的日前调度成本f
da
和第二阶段的随机场景下的实时调整成本期望f
rt
之和f最小为目标函数，如式(20)所示；第一部分为日前阶段的调度成本f
da
，包括常规火电机组运行成本热电联产机组运行成本弃风惩罚成本如式(21)所示；第二部分为实时阶段的调度成本f
rt
，包括常规火电机组备用调用成本期望值热电联产机组备用调用成本期望值非自愿切负荷和弃风成本期望值如式(22)所示；min f＝min(f
da
+f
rt
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(20)
式中，u
i,t
、y
i,t
、z
i,t
分别为常规火电机组i的运行状态变量、启动变量、停机变量；s
u,i
和s
d,i
分别为常规火电机组i单次启动成本和停机成本；a
1,i
、a
2,i
、a
3,i
为常规火电机组i的发电成本系数；p
g,i,t
为常规火电机组i的发电功率；和分别为常规火电机组i提供的向上备用和价格系数；和分别为常规火电机组i提供的向下备用和价格系数；p
chp,h,t
和h
chp,h,t
分别为热电联产机组h的发电功率和热功率；b
1,h
、b
2,h
、b
3,h
、b
4,h
、b
5,h
、b
6,h
为热电联产机组的运行成本系数；和分别为热电联产机组h提供的向上备用和价格系数；和分别为热电联产机组h提供的向下备用和价格系数；和p
w,t
为风电场w的预测功率和调度功率；c
wind
为弃风惩罚成本系数；和分别为常规火电机组i在场景s中调用的向上备用和调用能量价格系数；和分别为常规火电机组i在场景s中调用的向下备用和调用能量价格系数；和分别为热电联产机组h在场景s中调用的向上备用和调用能量价格系数；和分别为热电联产机组h在场景s中调用的向下备用和调用价格系数；和c
load
为场景s下的切负荷功率和切负荷惩罚成本系数；和p
w,t,s
为风电场w在场景s的实际功率和调度功率；n
chp
、n
g
、n
w
、n
load
、n
s
分别为热电联产机组、常规火电机组、风电场、用电负荷、随机场景的数量；p
s
为场景s的概率。5.根据权利要求4所述的提升灵活性的电热综合能源系统两阶段随机优化调度方法，其特征在于电热综合能源系统两阶段随机优化调度模型第一阶段日前调度约束条件是指预测场景下的日前调度约束条件，包括耦合元件运行约束、电网运行约束和热网运行约束，(1)耦合元件运行约束即热电联产机组运行约束，热电联产机组运行约束包括电出力
和热出力运行约束、热电联产机组爬坡约束以及热电机组调节灵活性备用约束；抽气式热电联产机组的运行域是一个多面体，其电出力和热出力运行约束如式(23)所示，式中：分别为热电联产机组h纯凝工况下的最大、最小的发电功率；c
v,h
为热电联产机组h总进气量恒定时抽取单位供热量下发电功率的减小量；c
m,h
为热电联产机组h背压状态运行下的电功率和热功率的弹性系数；为热电联产机组最大热出力；c
k,h
为常数项；热电联产机组爬坡约束如式(24)所示，式中，和分别为热电联产机组h发电功率的向上爬坡能力和向下爬坡能力；热电机组调节灵活性备用约束如式(25)至(26)所示，热电机组调节灵活性备用约束如式(25)至(26)所示，(2)电网运行约束包括电功率平衡约束、线路潮流约束和常规火电机组运行约束、常规火电机组爬坡约束、常规火电机组备用容量约束、常规火电机组最小运行时间约束、常规火电机组最小停机时间约束、常规火电机组开停机状态约束、风电场出力约束；1)电功率平衡约束如式(27)所示，式中，p
d,t
为负荷d的用电功率；2)线路潮流约束如式(28)所示，式中，λ
h,l
、λ
i,l
、λ
w,l
、λ
d,l
分别为热电联产机组h、常规火电机组i、风电场w、负荷d对网架输电线路l的潮流分布转移因子；f
l,max
为输电线路l的功率输送限值；3)常规火电机组运行约束如式(29)所示，p
g,i,min
≤p
g,i,t
≤p
g,i,max
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(29)式中，p
i,min
和p
i,max
分别为常规火电机组i的最小有功出力和最大有功出力；4)常规火电机组爬坡约束如式(30)所示，
式中，和分别为常规火电机组i的向上爬坡能力和向下爬坡能力；5)常规火电机组备用容量约束如式(31)所示，6)常规火电机组最小运行时间约束如式(32)所示，式中，t
u,i
和u
i,0
分别为常规机组i的最小启动运行时间和初始运行时间；7)常规火电机组最小停机时间约束如式(33)所示，8)常规火电机组开停机状态约束如式(34)所示，y
i,t-z
i,t
＝u
i,t-u
i,t-1
，y
i,t
+z
i,t
≤1
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(34)9)风电场出力约束如式(35)所示，(3)热网运行约束包括热源供热功率约束、热网动态运行约束、储热罐的热功率约束和建筑物热动态约束，热源供热功率约束如式(1)所示，热网动态运行约束如式(2)至(14)所示，储热罐的热功率约束如式(36)所示，式中：和为分别为第s个储热罐在t时段的储热状态和放热状态；h
hs,s,t
、
分别为第s个储热罐在t时段的热功率、储热功率和放热功率；分别为第s个储热罐在t时段的热功率、储热功率和放热功率；分别为第s个储热罐的最小储热功率、最大储热功率、最小放热功率、最大放热功率；e
hs,s,t
为第s个储热罐在t时段的储热量；e
hs,s,min
、e
hs,s,max
、e
hs,s,0
为第s个储热罐的最小储热量、最大储热量、初始储热量；建筑物热动态约束包括室内温度-建筑物的热量-调度时段耦合关系约束和室内温度范围约束；室内温度-建筑物的热量-调度时段耦合关系约束如式(37)所示，式中，τ
od,k,t
、τ
id,k,t
分别为建筑物k的室外温度、室内温度；χ
bt,k
、t
bhr,k
分别为建筑物k的单位温差传热系数、热储备系数；室内温度范围约束如式(38)所示，式中，分别不影响用户的用热舒适度下的最大室内温度、最小室内温度。6.根据权利要求5所述的提升灵活性的电热综合能源系统两阶段随机优化调度方法，其特征在于所述电热综合能源系统两阶段随机优化调度模型第二阶段约束条件基于第一阶段确定的日前调度计划，第二阶段为应对实时风电随机场景的不确定性，对日前调度计划进行优化调整，最小化实时调整的期望成本，并且在日内实时调度阶段时，不再调整热网中热源的供热计划，仅对供电计划进行优化调整，优化调整调整过程中，场景s下的电功率平衡约束、常规火电机组爬坡约束、热电联产机组爬坡约束、风电功率、切负荷约束、备用容量约束、线路潮流约束如下：1)场景s下的电功率平衡约束如式(39)所示，2)场景s下的常规火电机组爬坡约束如式(40)至(41)所示，2)场景s下的常规火电机组爬坡约束如式(40)至(41)所示，3)场景s下的热电联产机组爬坡约束如式(42)至(43)所示，3)场景s下的热电联产机组爬坡约束如式(42)至(43)所示，4)场景s下的风电功率、切负荷约束如式(44)所示，
5)备用容量约束如式(45)至(46)所示，5)备用容量约束如式(45)至(46)所示，线路潮流约束如式(47)所示，7.根据权利要求5或6所述的提升灵活性的电热综合能源系统两阶段随机优化调度方法，其特征在于基于所述电热综合能源系统两阶段随机优化调度模型，得到日前调度计划中电-热综合能源系统中常规火电机组的运行状态u
i,t
、发电功率p
g,i,t
、提供的向上备用容量向下备用容量热电联产机组的发电功率p
chp,h,t
、供热功率h
chp,h,t
、提供的向上备用容量向下备用容量储热罐的储热功率放热功率风电场的出力计划p
w,t
，用户建筑物内的温度τ
id,k,t
。

技术总结
本发明涉及电热综合能源系统优化运行调度方法技术领域，是一种提升灵活性的电热综合能源系统两阶段随机优化调度方法，其包括构建计及储热能力的热网动态特性模型和提升灵活性的电热综合能源系统两阶段随机优化调度模型，求解所述模型，得到电热综合能源系统调度计划。本发明考虑热力系统中的综合储热特性和灵活运行特性，进而实现电-热供应解耦，同时考虑电力系统中常规火电机组和热电联产机组的运行灵活性和提供备用的能力，挖掘热、电系统的潜在灵活性，能够从整个综合能源系统范围内实现灵活性资源的协调优化，提升整体运行的灵活性，为调度运行人员提供更合理的供热和供电调度计划。调度计划。调度计划。


技术研发人员：唐君毅 刘震 董雪涛 秦艳辉 朱子民 李德存 南东亮 孙冰 段青熙 段玉 王小云 祁晓笑 张媛 马星 糟伟红 马健 杨琪
受保护的技术使用者：国家电网有限公司
技术研发日：2021.12.03
技术公布日：2022/3/8