考虑随机停电供能可靠性的综合能源微网优化控制方法与流程

1.本发明涉及综合能源微网技术领域，尤其涉及一种考虑随机停电供能可靠性的综合能源微网优化控制方法。


背景技术：

2.随着多能耦合、集成优化、分布式能源发电等技术的推广应用，综合能源微网(integrated energy micro grid,iemg)发展迅速。iemg可以通过能源生产、转换和储存单元的协调互补，促进能源的高效利用，比独立的能源供应系统具有更好的能源供应经济性和可靠性。
3.iemg的优化控制需立足于当地的能源基础设施、自然条件、可获得的能源资源等，通过合理的供能单元选择和设备配置，提高用户供能的经济性和可靠性。iemg通常以并网模式运行，从外部电网(external power grid，epg)购买电能，并将其转换为其他形式的能源，如热能和冷能，以满足各种类型的负荷。但是，由于故障或检修等原因，epg可能会发生停电，无法向iemg输送电力。例如：2019年8月9日英国大停电和2020年美国加利福尼亚州轮流停电事故，严重影响了iemg的能源供应可靠性，尤其是对重要负荷。由于iemg的多能耦合特性，epg停电不仅影响其供电环节，还影响供气、供热环节的运行。
4.当epg停电时，iemg可以切换到孤岛模式运行，利用本地的调控资源(例如：分布式发电、储能设备等)维持能源供应。例如，在2011年日本大地震中，仙台的iemg被迫进入孤岛模式，采用电池、燃气发动机等装置实现了能源供应的快速恢复
1.。但在目前的优化控制方法中，没用充分挖掘iemg多能互补、灵活调控等特性对于抵御或减轻epg停电影响的优势，削弱了iemg的供能可靠性。
5.目前，在iemg优化控制方面已有大量的研究。文献[2]考虑了设备的投资约束，在满足用户能源需求的同时，通过合理的容量配置降低了iemg的运行成本。为了应对负荷和可再生能源的不确定性，文献[3]采用了一种iemg鲁棒规划方法，可以减少不确定性带来的成本变化。文献[4]提出了一种融合综合需求响应的iemg协同优化配置方法，该方法考虑了负荷的不确定性，可以降低能源利用成本。为了减少污染物排放，文献[5]建立了考虑环境外部性的iemg容量规划和优化模型，以促进可再生能源的使用。文献[6]同时考虑了经济性和环保性，建立了包括容量配置和运行优化的两阶段优化控制模型，可以降低iemg的用能和环境成本。
[0006]
上述研究侧重于降低iemg的优化控制成本，减轻能源价格、负荷、可再生能源等不确定因素对iemg未来运行的影响。然而，大多数研究忽略了epg停电，在并网模式下得到优化控制方案，却很少考虑iemg的孤岛运行模式。由于有epg的支撑，iemg发电设备的配置容量较小，当epg出现故障时，可能会造成能源供应短缺，带来经济损失。此外，由于epg的故障具有很强的随机性，电网的检修也具有不确定性，epg停电是一个随机事件，具有相互关联的多个不确定参数，如：持续时长、开始时间、发生日期等。在iemg的优化控制中，需要刻画这些不确定参数并对epg随机停电进行建模。
[0007]
因此，急需一种考虑随机停电供能可靠性的综合能源微网优化控制方法。


技术实现要素：

[0008]
本发明提供了一种考虑随机停电供能可靠性的综合能源微网优化控制方法，本发明解决了上述研究忽略epg停电对iemg运行的影响，导致iemg控制方案不合理，使iemg出现能量供给不足，削弱iemg供能可靠性等问题，详见下文描述：
[0009]
一种考虑随机停电供能可靠性的综合能源微网优化控制方法，所述方法包括：
[0010]
建立epg随机停电导致的iemg孤岛模式的概率模型，采用场景树生成若干孤岛场景；
[0011]
建立iemg多目标优化控制模型，通过对iemg进行协同优化控制，得到最优的长期优化控制方案和两种运行模式下的短期运行方案，用于提升iemg并网模式的运行经济性及孤岛模式的供能可靠性。
[0012]
其中，所述建立epg随机停电导致的iemg孤岛模式的概率模型为：
[0013]
采用随机向量ξ表示孤岛模式，由孤岛类型i、持续时间t
dur
、发生日期d
occ
和开始时间t
start
组成；
[0014]
其中，随机向量ξ的概率密度函数的联合pdf，定义如下：
[0015][0016]
式中，f(ξ)表示ξ的pdf；表示i、t
dur
、d
occ
和t
start
的联合pdf；i、t
dur
、d
occ
和t
start
分别表示i、t
dur
、d
occ
和t
start
的取值；fi(i)表示i的pdf；表示t
dur
的条件pdf；表示d
occ
的条件pdf；表示t
start
的条件pdf。
[0017]
进一步地，所述采用场景树生成若干孤岛场景为：
[0018]
采用场景树生成孤岛场景，从场景树的l1层的根节点开始，根据3种孤岛模式，在l2层生成3个节点；对于l2的每个节点，根据孤岛类型确定后持续时间的条件概率，生成l3层的节点；确定孤岛类型和持续时间后，根据发生日期的条件概率，生成l4层的节点；同理，考虑开始时间的条件概率，生成l5层的叶节点；从根节点到叶节点的路径被定义为孤岛场景；
[0019]
最终，生成s
total
个孤岛场景，包括si个i型孤岛场景，s
ii
个ii型孤岛场景和s
iii
个iii型孤岛场景，场景s的发生概率为：
[0020][0021][0022]
式中，p(s)表示场景s的发生概率；λi表示i型孤岛的发生频率；和分别表示和的条件概率；tdi、doi和tsi表示持续时间集、发生日期集和开始时间集。
[0023]
其中，所述多目标优化控制模型的目标函数包括：
[0024]
第一目标函数表示为最小化iemg的优化控制总成本，包括投资费用、运行费用、维
护费用和碳税费用；
[0025]
第二目标函数表示为孤岛模式下iemg的最小经济损失期望，包括能量供给不足惩罚和由于设备出力过多产生的弃能惩罚；
[0026]
采用线性加权和法将上述两个目标函数转化为加权目标函数，通过调整权重系数，以协调经济性和可靠性的关系。
[0027]
进一步地，所述多目标优化控制模型的约束条件包括：
[0028]
优化控制模型中的约束条件包括设备配置容量约束、并网模式下的运行约束及孤岛模式下的“可靠运行策略”及运行约束；
[0029]
其中，并网模式下的运行约束包括：
[0030]
与epg交换电能约束、设备运行约束和功率平衡约束；
[0031]
其中，所述孤岛模式下的可靠运行策略如下：
[0032]
步骤1：初始化运行情况；
[0033]
步骤2：判断是否chp、pv或chp、gb和hp满足电或热负荷需求，若是，执行步骤2-1，若否，执行步骤2-2；
[0034]
步骤2-1：将多余电或热能存储在es或hs中，并判断es或hs是否达到能量储存上限，若是，产生弃电或热；若否，结束；
[0035]
步骤2-2：调用es或hs放电或热；
[0036]
步骤3：判断是否满足电或热负荷需求，若是，结束，若否，执行步骤4；
[0037]
步骤4：优先切除nel或nhl，再切除pel或phl；
[0038]
步骤5：结束。
[0039]
本发明提供的技术方案的有益效果是：
[0040]
1、本发明通过目标函数和约束条件设定，所提的多目标优化控制模型不仅优化了长时间尺度控制方案，还优化了不同模式下的短期运行策略，提高了并网模式下的运行经济性，并提高了孤岛模式的供能可靠性；
[0041]
2、本发明采用场景树刻画了iemg未来可能出现的孤岛场景，通过合理的优化控制进一步挖掘了多能互补的优势，以牺牲较低的经济性为代价，大幅降低了epg停电造成的经济损失期望，大大提高了iemg的可靠性；
[0042]
3、本发明相比于独立供能系统的控制方法，优化控制总成本和经济损失期望得到了大幅降低，提高了iemg的供能经济性和可靠性。
附图说明
[0043]
图1为一种考虑随机停电供能可靠性的综合能源微网优化控制方法示意图；
[0044]
图2为基于场景树的孤岛生成方法示意图；
[0045]
图3为iemg规划-运行多目标协同优化控制模型示意图；
[0046]
图4为iemg孤岛模式“可靠运行策略”流程图；
[0047]
图5为iemg结构图；
[0048]
图6为负荷及光照强度曲线示意图；
[0049]
图7为分时电价曲线示意图；
[0050]
图8为各算例费用对比图；
[0051]
图9为不同类型孤岛的ens曲线；
[0052]
图10为case i在ii型孤岛的运行情况曲线；
[0053]
图11为case ii在ii型孤岛的运行情况曲线；
[0054]
图12为case iii在ii型孤岛的运行情况曲线。
具体实施方式
[0055]
为使本发明的目的、技术方案和优点更加清楚，下面对本发明实施方式作进一步地详细描述。
[0056]
本发明实施例提供了一种考虑随机停电供能可靠性的综合能源微网优化控制方法，建立epg随机停电导致的iemg孤岛模式的概率模型，采用场景树生成若干孤岛场景，进而建立iemg多目标优化控制模型，得到系统整体的优化控制策略，用于提升iemg并网模式的运行经济性及孤岛模式的供能可靠性。
[0057]
首先，根据epg停电的原因和严重程度将孤岛模式分为3种。考虑孤岛类型、持续时间、发生日期和开始时间4个不确定参数及它们之间的相互限制及约束，建立孤岛模式的概率模型。在此基础上，将概率模型离散化，通过场景树生成孤岛场景。进而，建立iemg多目标优化控制模型，该模型有2个目标函数：第1个目标函数为最小化优化控制总成本，以优化规划方案及并网模式的运行方案；第2个目标函数是最小化孤岛模式下的经济损失期望，旨在提高能源供应的可靠性。优化控制模型中的约束条件包括设备配置容量约束、并网模式下的运行约束及孤岛模式下的“可靠运行策略”及运行约束。该优化控制模型融合了长期优化控制方案和两种运行模式下的短期运行方案，通过协同优化控制，得到最优的结果。
[0058]
如图1所示，本发明实施例提供的一种考虑随机停电供能可靠性的综合能源微网优化控制方法，包括如下步骤：
[0059]
1)对iemg孤岛模式建模，并采用场景树生成若干孤岛场景。
[0060]
根据epg停电的原因及严重程度，将iemg的孤岛模式分为3类：
[0061]
i型孤岛：i型孤岛是由epg易于修复的小规模故障引起，持续时间较短，一般为1～5小时，假设其发生的频率为λi(次/a)。i型孤岛的持续时间、发生日期和开始时间是不确定的。
[0062]
ii型孤岛：ii型孤岛是由epg非常规、规模较大的故障引起的，故障修复时间长，一般为十几个小时甚至几天。其发生频率较低，假设为λ
ii
(次/a)，其持续时间、发生日期和开始时间是不确定的。
[0063]
iii型孤岛：iii型孤岛是由epg预安排停电引起，其发生频率假设为λ
iii
(次/a)，持续时间、发生日期和开始时间是随机的，但会提前通知iemg。
[0064]
采用随机向量ξ表示孤岛模式，它由4个随机变量组成：孤岛类型i、持续时间t
dur
、发生日期d
occ
和开始时间t
start
。
[0065]
其中，ξ定义如下：
[0066]
ξ＝(i,t
dur
,d
occ
,t
start
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0067]
例如，ξ＝(i型,3小时,6月5日,6:00)表示i型孤岛发生在6月5日的6:00，持续3个小时。
[0068]
ξ的概率密度函数(probability density function,pdf)，即i、t
dur
、d
occ
和t
start
的
联合pdf，定义如下：
[0069][0070]
式中，f(ξ)表示ξ的pdf；表示i、t
dur
、d
occ
和t
start
的联合pdf；i、t
dur
、d
occ
和t
start
分别表示i、t
dur
、d
occ
和t
start
的取值；fi(i)表示i的pdf；表示t
dur
的条件pdf，已知i＝i；表示d
occ
的条件pdf，已知i＝i且t
dur
＝t
dur
；表示t
start
的条件pdf，已知i＝i,t
dur
＝t
dur
且d
occ
＝d
occ
。
[0071]
如式(2)所示，4个随机变量不是独立的，它们是关联的，条件pdf表明后一随机变量受前一随机变量的约束和影响。
[0072]
采用场景树生成孤岛场景，如图2所示，从场景树的l1层的根节点开始，根据3种孤岛模式，在l2层生成3个节点。对于l2的每个节点，根据孤岛类型确定后持续时间的条件概率，生成l3层的节点。确定孤岛类型和持续时间后，根据发生日期的条件概率，生成l4层的节点。同理，考虑开始时间的条件概率，生成l5层的叶节点。从根节点到叶节点的路径被定义为孤岛场景。
[0073]
最终，生成s
total
个孤岛场景，包括si个i型孤岛场景，s
ii
个ii型孤岛场景和s
iii
个iii型孤岛场景。以场景s为例，s∈{1,2,
…
,s
total
}，假设其对应孤岛类型为i，持续时间为发生日期为开始时间为场景s的发生概率为：
[0074][0075]
式中，p(s)表示场景s的发生概率；λi表示i型孤岛的发生频率；和分别表示和的条件概率；tdi、doi和tsi表示持续时间集、发生日期集和开始时间集，包括持续时间、发生日期和开始时间的可能值。事实上，式(3)是式(2)的离散形式。
[0076]
2)建立iemg多目标优化控制模型，通过对iemg进行协同优化控制，得到最优的长期优化控制方案和两种运行模式下的短期运行方案，如图3所示。
[0077]
(1)iemg多目标优化控制模型整体框架
[0078]
首先将设备参数、负荷、光照强度、能源价格、孤岛场景等基础数据输入到优化控制模型中。
[0079]
规划-运行协同优化控制模型有两个目标函数：目标函数1是最小化iemg的优化控制总成本，包括投资费用、运行费用、维护费用和碳税费用，其目的是优化规划方案及并网模式下的运行方案。目标函数2的目的是通过优化孤岛模式运行方案，确保iemg尽可能满足能源需求，为与目标函数1相匹配，目标函数2表示为孤岛模式下iemg的最小经济损失期望，包括能量供给不足(energy not supplied,ens)惩罚和由于设备出力过多产生的弃能(energy abandoned,ea)惩罚。采用线性加权和法将上述两个目标函数转化为加权目标函数。调整加权目标函数中的权重系数ω，以协调经济性和可靠性的关系。
[0080]
优化控制模型的变量和约束分为3个模块。在规划模块，变量用向量x表示，包括每
个设备的容量，它们应满足设配配置容量的上下限约束。x还为两个运行模块提供设备运行限制。在并网运行模块，变量用向量y表示，包括各设备的出力和从epg购买的电能，它们应满足设备运行约束和功率平衡约束。孤岛运行模块中的变量由向量u表示，包含每个设备的出力、ens和ea。为了区分输入的不同孤岛场景，场景s中的变量表示为向量us。为了减少ens和ea带来的经济损失期望，孤岛模式采用“可靠运行策略”，iemg应满足相应的设备运行约束、功率平衡约束和运行逻辑约束。
[0081]
最后，将变量传递到目标函数的计算中，目标函数1采用x和y计算，目标函数2采用各场景的us计算。求解规划-运行协同优化模型，输出优化结果，包括规划方案及并网模式和孤岛模式下的运行方案。
[0082]
该多目标优化控制模型是一个milp问题，用矩阵形式表示如下：
[0083][0084]
式中，f(x,y)表示优化控制总成本；e[g(u,ξ)]表示孤岛模式下的经济损失期望；e[]表示期望计算公式；ω是权重系数，ω∈(0,1]；c、d和q为目标函数的系数列向量；向量us表示孤岛模式在场景s中的运行变量；a、b、c、d、e和f为约束中变量对应的系数矩阵；a、b、m和n为约束的系数列向量；xu表示x的上限。
[0085]
在优化控制模型的约束条件中，式(a)表示规划中设备的配置容量限制；式(b)表示iemg并网运行时的不等式约束，包括设备出力的上下限；式(c)表示并网运行时的等式约束，包括各设备的能量转换约束和功率平衡设备；对应地，式(d)式(e)表示iemg在孤岛运行时的不等式约束和等式约束，也包括“可靠运行策略”需满足的运行逻辑约束。
[0086]
(2)iemg多目标优化控制模型目标函数
[0087]
优化控制模型中目标函数1是最小化iemg的优化控制总成本，将其表示为等年值：
[0088]
obj1:min c
total
＝c
inv
+c
ope
+c
mai
+c
co2
ꢀꢀꢀꢀꢀꢀ
(5)
[0089]
式中，c
total
是iemg的优化控制总成本；c
inv
是设备投资费用等年值；c
ope
是并网模式下的年运行费用；c
mai
是并网模式下的年维护费用；c
co2
是并网模式下的年碳税费用。
[0090]
假设在优化控制中有n种候选设备，c
inv
计算如下：
[0091][0092][0093]
式中，c
inv,m
表示设备m单位容量的投资费用；wm表示设备m的配置容量；nm表示设备m的寿命；r表示折现率；表示年值系数，它的作用是将设备m的投资费用换算成其寿命内
的等年值；m＝1,2
…
,n。
[0094]cope
包括并网模式下的购电费用和天然气消耗费用，计算公式如下：
[0095]cope
＝c
grid
p
gridt
+sum(c
gasggas
)
ꢀꢀꢀꢀꢀꢀꢀ
(8)
[0096]
式中，c
grid
＝[c
grid
(t)]1×
8760
表示电价向量；c
grid
(t)为t时刻的电价；p
grid
＝[p
grid
(t)]1×
8760
表示在并网模式下从epg购电功率向量，p
grid
(t)为并网模式下t时刻的购电功率；c
gas
表示天然气价格；g
gas
＝[g
gas
(t)]1×
8760
表示并网模式下天然气消耗功率向量，g
gas
(t)为并网模式下t时刻的天然气消耗功率；sum()表示求和函数。
[0097]cmai
的表达式为：
[0098]cmai
＝sum(c
mai
p)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(9)
[0099]
式中，c
mai
＝[c
mai,m
]1×m表示设备单位功率可变维护费用向量，c
mai,m
为设备m单位功率可变维护费用；p＝[pm(t)]m×
8760
表示在并网模式下的设备输出功率矩阵，pm(t)为在并网模式下t时刻设备m的输出功率。
[0100]cco2
包括从epg购电产生的碳税费用和消耗天然气产生的碳税费用：
[0101]cco2
＝c
co2
(sum(μ
grid
p
grid
)+sum(μ
gasggas
))
ꢀꢀꢀꢀꢀ
(10)
[0102]
式中，c
co2
是每kg co2排放的碳税费用；μ
grid
和μ
gas
分别表示从epg购电以及消耗天然气对应的co2排放系数。
[0103]
目标函数2是最小化孤岛模式下iemg的经济损失期望，尽可能减少ens和ea，计算公式如下：
[0104][0105]
式中，c
loss
表示孤岛模式下的预期经济期望；gs(us)表示场景s中的经济损失，计算公式如下：
[0106][0107][0108][0109][0110][0111]
式中，表示场景s的经济损失；c
ens
表示单位功率ens的惩罚费用向量；和分别表示pel、nel、phl和nhl的单位功率ens的惩罚费用；表示t时刻ens功率向量；和分别表示t时刻pel、nel、phl和nhl的ens功率；c
ea
表示单位功率ea的惩罚费用向量；和分别表示电和热的单位功率ea惩罚费用；表示t时刻ea功率向量；和分别表示t时刻电和热的ea功率。
[0112]
通过改写式(11)-(16)，目标函数2也可以表示成下列形式：
[0113]
obj2:min c
loss
＝c
ens
eens+c
ea
eea
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(17)
[0114]
eens＝[eens
e1
,eens
e2
,eens
h1
,eens
h2
]
t
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(18)
[0115]
eea＝[eeae,eeah]
t
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(19)
[0116][0117][0118]
式中，eens表示能量供给不足期望(expected energy not supplied,eens)向量，eens
e1
、eens
e2
、eens
h1
和eens
h2
分别表示pel、nel、phl和nhl的eens；eea表示弃能期望(expected energy abandoned,eea)向量，eeae和eeah分别表示电和热的eea。
[0119]
以上两个目标函数通过线性加权和法转化为如下加权目标函数：
[0120]
obj:minωc
total
+(1-ω)c
loss
ꢀꢀꢀꢀꢀꢀ
(22)
[0121]
式中，ω的值取决于决策者的偏好。
[0122]
(3)iemg多目标优化控制模型约束条件
[0123]
优化控制模型中并网模式下的运行约束具体为：
[0124]
在并网模式下，iemg可以与epg进行电能交换，由于epg的支持，iemg可以为所有负荷供电而不会出现能源短缺。iemg运行在并网模式应满足与epg交换电能约束、设备运行约束和功率平衡约束。
[0125]
a)与epg交互电能功率约束：
[0126][0127]
式中，表示epg可向iemg传输的最大功率。
[0128]
b)pv运行约束：
[0129][0130]
0≤p
pv,c
(t)≤p
pv
(t)(24b)
[0131]
式中，p
pv
(t)为t时刻pv的最大输出功率；η
pv
为pv发电转换效率系数；w
pv
表示pv的配置容量；hs(t)表示t时刻的光照辐射密度；h
stc
表示标准测试条件下的光照辐射密度；p
pv,c
(t)为t时刻pv的消纳功率。
[0132]
c)能量生产、转换设备约束(chp、hp和gb)：
[0133]
p
et,out
(t)＝p
et,in
(t)η
et
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(25a)
[0134]
0≤p
et,out
(t)≤w
et et∈et＝{chp,hp,gb}
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(25b)
[0135]
式中，p
et,out
(t)表示t时刻设备et的输出能量功率；p
et,in
(t)表示t时刻设备et消耗的输入能量功率；η
et
表示设备et的能量转换效率；w
et
表示设备et的配置容量；et表示能量生产、转换设备集合。
[0136]
d)储能设备约束(es和hs)：
[0137]
0≤pc(t)≤p
c,max
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(26a)
[0138]
0≤pd(t)≤p
d,max
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(26b)
[0139][0140]
soc
min
≤soc(t)≤soc
max
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(26d)
[0141]
pc(t)pd(t)＝0
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(26e)
[0142]
soc(0)＝soc(t)
ꢀꢀꢀꢀꢀꢀ
(26f)
[0143]
式中，pc(t)和pd(t)分别表示t时刻ess(包括es和hs)充能功率与放能功率；p
c,max
和p
d,max
为最大充能、放能功率；soc(t)为ess在t时刻的荷电状态；soc
max
、soc
min
为ess的soc上下限；ηc和ηd分别为充能和放能效率；w
ess
为ess的配置容量；soc(0)、soc(t)分别为ess调度期初和期末的soc；δt为调度周期，本发明实施例取1小时。
[0144]
e)电/热功率平衡约束：
[0145][0146][0147]
式中，p
chp
(t)表示t时刻chp的发电功率；和分别表示t时刻es的充放电功率；p
hp
(t)表示t时刻hp的耗电功率；p
l1
(t)和p
l2
(t)分别表示t时刻重要电负荷和一般电负荷功率；h
chp
(t)表示t时刻chp的供热功率；h
hp
(t)表示t时刻hp的供热功率；h
gb
(t)表示t时刻gb的供热功率；和分别表示hs在t时刻的充放热功率；h
l1
(t)和h
l2
(t)分别表示t时刻重要热负荷和一般热负荷功率。
[0148]
优化控制模型中孤岛模式下的“可靠运行策略”及运行约束具体为：
[0149]
在孤岛模式下，为降低ens，iemg采用“可靠运行策略”，如图4所示，依靠自身的能源供应能力，以尽可能满足电/热负荷需求为目标，通过各环节的协调互补，控制各设备的出力，具体步骤如下：
[0150]
步骤1：初始化运行情况；
[0151]
步骤2：判断是否chp、pv(或chp、gb和hp)满足电(或热)负荷需求，若是，执行步骤2-1，若否，执行步骤2-2；
[0152]
步骤2-1：将多余电(或热)能存储在es(或hs)中，并判断es(或hs)是否达到能量储存上限，若是，产生弃电(或热)；若否，结束；
[0153]
步骤2-2：调用es(或hs)放电(或热)；
[0154]
步骤3：判断是否满足电(或热)负荷需求，若是，结束，若否，执行步骤4；
[0155]
步骤4：优先切除nel(或nhl)，再切除pel(或phl)；
[0156]
步骤5：结束。
[0157]
以场景s为例，在的区间内，iemg处于孤岛模式，为实现上述运行策略，除需满足式(5-24)-(5-26)的设备运行约束外，供电环节还需满足的约束为：
[0158]
a)iemg与epg的交换功率为0：
[0159]
p
grid
(t)＝0
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(29)
[0160]
b)除式(5-26)外，es的运行约束：
[0161]
p
chp
(t)+p
pv
(t)-p
l1
(t)-p
l2
(t)≥-mzd,zd∈{0,1}
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(30a)
[0162]
p
chp
(t)+p
pv
(t)-p
l1
(t)-p
l2
(t)≤m(1-zd)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(30b)
[0163][0164][0165]
式中，m表示一个充分大的正数；zd为0-1变量，用于表示chp和pv是否不满足电负荷需求，即是否需要调用es放电。
[0166]
c)弃电约束：
[0167]
1-zd≥za,za∈{0,1}
ꢀꢀꢀꢀ
(31a)
[0168][0169][0170]
式中，za为0-1变量，用于表示是否需要产生弃电。
[0171]
d)切负荷优先级约束：
[0172]
zd≥zn,zn∈{0,1}
ꢀꢀꢀꢀꢀꢀꢀ
(32a)
[0173][0174][0175][0176]zn
≥z
p
,z
p
∈{0,1}
ꢀꢀꢀꢀ
(32e)
[0177][0178][0179][0180]
式中，zn和z
p
为0-1变量，分别表示是否需要nel和pel；式(32e)表示优先切除nel。
[0181]
e)功率平衡约束
[0182][0183]
供热环节的相应约束与供电环节的约束类似，此处不再赘述。
[0184]
考虑外部电网随机停电供能可靠性的iemg多目标优化控制方法通过场景树刻画了iemg未来可能出现的孤岛场景，并在优化控制模型中将控制总成本最低及epg随机停电导致的经济损失期望最低作为目标函数，设置并网和孤岛模式下的运行策略及约束，通过合理的优化控制进一步挖掘了iemg多能互补的优势，不仅优化了控制方案，还优化了不同模式下的运行策略，提高了并网模式下的运行经济性，并提高了孤岛模式的供能可靠性。较传统仅考虑经济性的单目标优化控制方法，该方法以牺牲较低的经济性为代价，大幅降低了epg停电造成的经济损失期望，大大提高了iemg的可靠性。
[0185]
下面给出实施例。
[0186]
(1)实施例及参数
[0187]
以中国北方的一个iemg为例，如图5所示。假设iemg的优化控制周期为15年，候选设备及其参数见表1和表2。考虑iemg孤岛模式下ens和ea造成的经济损失，假设为290元/kw、为70元/kw、为60元/kw、为20元/kw。另外，假设为为1.35元/
kw，则为1.2元/kw。为反映负荷和可再生能源的波动性和季节性差异，采用夏季、冬季和过渡季三个典型日和夏季和冬季两个极端典型日的数据来表示全年情况，一年中各典型日的天数为100、109、152、2和2。各典型日负荷曲线和光照强度曲线如图6所示。其中，夏季极端典型日用电负荷比正常高10％，而光伏出力仅为70％；冬季极端典型日热负荷比正常高10％。
[0188]
表1能量生产和转换设备参数
[0189][0190]
表2 ess参数
[0191][0192]
iemg采用分时电价，如图7所示。天然气采用固定气价，为2.71元/m3，天然气低热值为9.7kwh/m3，经过折算后的天然气价格为0.28元/kwh。从外部电网购电和消耗天然气对应的co2排放系数分别为0.137kg/kwh和0.198kg/kwh，单位kg的碳税费用为0.02元。折现率r取8％。
[0193]
epg停电数据和孤岛模式参数由当地供电公司提供，采用场景树生成孤岛场景，如表3所示。
[0194]
表3孤岛模式参数和场景生成
[0195][0196]
为对比本发明实施例的多目标优化控制方法相较于仅考虑iemg经济性的单目标优化控制方法的有效性，以及采用多能耦合供能方式相较于传统分供方式在提升系统运行经济性以及供能可靠性方面的有效性，设置3个算例如下：
[0197]
case i：采用多目标优化控制，将iemg优化控制总成本最低和孤岛模式经济损失期望最低作为目标函数，设权重系数ω为0.5；系统采用多能耦合供能方式，候选设备包括chp、pv、gb、hp、es和hs。
[0198]
case ii：采用单目标优化控制，仅将iemg优化控制总成本最低作为目标函数；系
统采用多能耦合供能方式。
[0199]
case iii：采用多目标优化控制，设权重系数ω为0.5；系统采用传统分供方式，即iemg并网模式下电能通过外部电网供应，热能通过gb供应；在iemg孤岛模式下采用备用的柴油发电机(diesel generator,dg)供电。
[0200]
(2)优化控制经济性及可靠性分析
[0201]
iemg优化控制总成本c
total
、孤岛模式的经济损失期望c
loss
及其构成如图8所示。case i和case ii相比，case ii的c
total
比case i高10600元，仅占0.17％。case ii的c
loss
大约是case i的2.6倍，比case i出33800元。虽然case i的c
total
略高，但c
loss
低很多。因此，case i采用的多目标优化控制方法考虑了iemg的孤岛模式，能够有效抵御epg停电的影响，增强iemg供能可靠性。
[0202]
对比case i和case iii，虽然case iii的投资成本远低于case i，但case iii并网模式运行时高度依赖epg，运营成本约为case i的2.5倍。case iii的总成本比case i增加约28,357万元，增加了31％。当epg出现故障时，case iii缺乏epg的支持，只能通过备用dg供电，造成能源供应严重不足，c
loss
达到40200元，约为case i的2倍。因此，采用多能耦合的iemg，不仅可以降低供能成本，还可以提高epg停电下iemg的供能可靠性。
[0203]
各算例的eens如表4所示，其显示了孤岛模式下iemg的供能可靠性。由于case i和case ii的多能互补效应以及case iii中dg的备用支撑，3个算例都具有一定的可靠性，都可以完全供应pel和phl。由于epg停电对供热影响较小，case i和case iii也可以完全供应nhl，case ii只有少量的nhl供能不足。对于nel，三种情况的eens
e2
大小不同，其中case i最低，case iii次之，case ii最高。这说明采用多目标优化控制可以进一步挖掘iemg多能源耦合及协调互补的潜力，降低eens
e2
，提高孤岛模式下的供能可靠性。
[0204]
表4eens结果
[0205][0206]
进一步，比较3个算例nel的ens，如图9(a)-(b)所示。图9(a)后半部分显示了iemg在冬季极端典型日17:00-19:00运行在i型孤岛的ens，在此期间电负荷较高，3个算例均出现了ens。case i和case ii只有18:00和19:00出现ens，且case i的ens小于case ii，说明epg停电对case i的影响较小。而case iii在整个时段都有ens，且ens值较高，受epg停电影响较大。
[0207]
图9(a)的前半部分显示了iemg运行在iii型孤岛的ens。由于ⅲ型孤岛发生在负荷较低的过渡季典型日，对iemg供电的影响较小。case i和case ii可以确保在孤岛模式下可靠供电，没有ens。case iii在11:00-14:00出现ens，但数值较低。
[0208]
图9(b)显示了iemg在整个夏季极端典型日运行在ii型孤岛时的ens，此时电负荷较高，但pv的输出较低。case ii和case iii出现了ens，case ii中ens出现的时间段为6:00-9:00和22:00-0:00，持续时间为6小时；而case iii中ens出现的时间相对集中，在9:00-23:00，持续时间为13小时。case ii和case iii的ens都很高，但case ii的总ens低于case iii。而此时对case i的影响较小，case i中没有ens。
[0209]
(3)孤岛模式运行情况对比
[0210]
本小结分析3个算例在孤岛模式下供能情况不同的原因，由于夏季极端典型日电负荷高而光伏出力低，如果发生持续一整天的ii型孤岛，将对iemg的电力供应构成巨大威胁，该场景可以很好地反映不同优化控制方案维持能源供应的能力。因此，以该孤岛场景为例，3个算例的设备出力、ens、ea、供电和供热情况如图10-图12所示。
[0211]
通过所提的多目标优化控制方法，多能互补在提高供能可靠性方面的优势得到了进一步挖掘。如图10(a)所示，case i中chp的配置容量充足，当iemg处于孤岛模式时，chp的发电量较大，以弥补pv的不足出力，并且通过es的充放电，可以完全满足iemg的用电需求，实现孤岛模式下的可靠供电。由于chp工作在“以电定热”模式下，为满足电负荷的供给，chp的负载率较高，产生多余的热能，由于夏季热负荷较低，iemg通过弃热来保证热能供需平衡，如图10(b)所示。
[0212]
case ii的chp容量低于case i，削弱了iemg的电力自给能力。如图11(a)所示，chp在1:00-6:00满载运行，但其输出功率仍低于电力需求，通过es放电以满足用电需求。在6:00-9:00，es达到了储量下限，无法继续放电，此时电能供应不足。随后，随着光照强度的增加，pv的输出功率逐渐增加，可以保证iemg充足的电能供应，并将多余的电能储存在es中。在17:00-22:00，随着pv出力的减少，iemg仍然可以依靠chp的出力和es的放电来保证充足的电能供应。但22:00后，随着es再次达到储量下限，chp的出力无法满足用电需求，导致ens，部分nel被切除。由于chp全天运行在高负载率状态，iemg采取弃热措施以保证热能供需平衡，如图11(b)所示。因此，虽然case ii具有一定的可靠性，但由于优化控制中仅考虑经济性，而忽略了epg停电的影响，所得方案不合理，iemg的多能互补作用受到了限制。
[0213]
case iii的能源供应模式较简单，如图12(a)所示。当epg停电时，系统利用备用的dg供电，可满足大部分用电需求。而当电负荷大于dg出力上限时，如9:00-20:00，出现了供电不足，必须削减高于dg出力的nel。对于供热环节，如图12(b)所示，根据供热需求调节gb的输出功率，可以保证供热环节实时供需平衡，不需要采取弃热措施。
[0214]
总之，孤岛模式下，case ii和case iii的设备出力和调控无法支撑iemg用电需求，导致电能供应不足。而case i在优化控制阶段考虑了未来可能发生的epg停电场景，通过优化各设备的容量，进一步挖掘了多能互补在可靠性提升方面的潜力，使iemg能够有效应对孤岛模式运行的挑战，确保能源的可靠供应。
[0215]
综上，通过实施例表明，本发明实施例的一种考虑随机停电供能可靠性的综合能源微网优化控制方，对epg停电及其随机性进行建模，并通过场景树生成多个孤岛场景。考虑优化控制经济性和运行可靠性要求，将iemg的优化控制总成本和epg停电造成的经济损失期望作为优化控制模型的目标函数，对iemg的规划决策及并网模式和孤岛模式的运行方案进行协同优化，所得结论如下：
[0216]
a)本发明通过目标函数和约束条件设定，不仅优化了控制方案，还优化了不同模式下的运行策略，提高了并网模式下的运行经济性，并提高了孤岛模式的供能可靠性。
[0217]
b)与只考虑经济性的iemg单目标优化控制方法不同，本发明通过场景树刻画了iemg未来可能出现的孤岛场景，通过合理的优化控制进一步挖掘了多能互补的优势。以牺牲较低的经济性为代价，大幅降低了epg停电造成的经济损失期望，改善了iemg的可靠性。
[0218]
c)本发明相比于独立供能系统的控制方法，优化控制总成本和经济损失期望得到
了大幅降低，提高了iemg的供能经济性和可靠性。
[0219]
参考文献：
[0220]
[1]marnay c,aki h,hirose k,et al.japan's pivot to resilience:how two microgrids fared after the2011 thquake[j].ieee power and energy magazine,2015,13(3):44-57.
[0221]
[2]wang y,li r,dong h,et al.capacity planning and optimization of business park-level integrated energy system based on investment constraints[j].energy,2019,189:116345.
[0222]
[3]yu n,kang j s,chang c c,et al.robust economic optimization and environmental policy analysis for microgrid planning:an application to taichung industrial park,taiwan[j].energy,2016,113:671-682.
[0223]
[4]liu w,huang y,li z,et al.optimal allocation for coupling device in an integrated energy system considering complex uncertainties of demand response[j].energy,2020,198:117279.
[0224]
[5]sun j,ruze n,zhang j,et al.capacity planning and optimization for integrated energy system in industrial park considering environmental externalities[j].renewable energy,2021,167:56-65.
[0225]
[6]wang y,wang y,huang y,et al.planning and operation method of the regional integrated energy system considering economy and environment[j].energy,2019,171:731-750.
[0226]
本发明实施例对各器件的型号除做特殊说明的以外，其他器件的型号不做限制，只要能完成上述功能的器件均可。
[0227]
本领域技术人员可以理解附图只是一个优选实施例的示意图，上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。
[0228]
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

技术特征：
1.一种考虑随机停电供能可靠性的综合能源微网优化控制方法，其特征在于，所述方法包括：建立epg随机停电导致的iemg孤岛模式的概率模型，采用场景树生成若干孤岛场景；建立iemg多目标优化控制模型，通过对iemg进行协同优化控制，得到最优的长期优化控制方案和两种运行模式下的短期运行方案，用于提升iemg并网模式的运行经济性及孤岛模式的供能可靠性。2.根据权利要求1所述的一种考虑随机停电供能可靠性的综合能源微网优化控制方法，其特征在于，所述建立epg随机停电导致的iemg孤岛模式的概率模型为：采用随机向量ξ表示孤岛模式，由孤岛类型i、持续时间t
dur
、发生日期d
occ
和开始时间t
start
组成；其中，随机向量ξ的概率密度函数的联合pdf，定义如下：式中，f(ξ)表示ξ的pdf；表示i、t
dur
、d
occ
和t
start
的联合pdf；i、t
dur
、d
occ
和t
start
分别表示i、t
dur
、d
occ
和t
start
的取值；f
i
(i)表示i的pdf；表示t
dur
的条件pdf；表示d
occ
的条件pdf；表示t
start
的条件pdf。3.根据权利要求1所述的一种考虑随机停电供能可靠性的综合能源微网优化控制方法，其特征在于，所述采用场景树生成若干孤岛场景为：采用场景树生成孤岛场景，从场景树的l1层的根节点开始，根据3种孤岛模式，在l2层生成3个节点；对于l2的每个节点，根据孤岛类型确定后持续时间的条件概率，生成l3层的节点；确定孤岛类型和持续时间后，根据发生日期的条件概率，生成l4层的节点；同理，考虑开始时间的条件概率，生成l5层的叶节点；从根节点到叶节点的路径被定义为孤岛场景；最终，生成s
total
个孤岛场景，包括s
i
个i型孤岛场景，s
ii
个ii型孤岛场景和s
iii
个iii型孤岛场景，场景s的发生概率为：孤岛场景，场景s的发生概率为：式中，p(s)表示场景s的发生概率；λ
i
表示i型孤岛的发生频率；和分别表示和的条件概率；td
i
、do
i
和ts
i
表示持续时间集、发生日期集和开始时间集。4.根据权利要求1所述的一种考虑随机停电供能可靠性的综合能源微网优化控制方法，其特征在于，所述多目标优化控制模型的目标函数包括：第一目标函数表示为最小化iemg的优化控制总成本，包括投资费用、运行费用、维护费用和碳税费用；第二目标函数表示为孤岛模式下iemg的最小经济损失期望，包括能量供给不足惩罚和由于设备出力过多产生的弃能惩罚；
采用线性加权和法将上述两个目标函数转化为加权目标函数，通过调整权重系数，以协调经济性和可靠性的关系。5.根据权利要求1所述的一种考虑随机停电供能可靠性的综合能源微网优化控制方法，其特征在于，所述多目标优化控制模型的约束条件包括：优化控制模型中的约束条件包括设备配置容量约束、并网模式下的运行约束及孤岛模式下的“可靠运行策略”及运行约束；其中，并网模式下的运行约束包括：与epg交换电能约束、设备运行约束和功率平衡约束；其中，所述孤岛模式下的可靠运行策略如下：步骤1：初始化运行情况；步骤2：判断是否chp、pv或chp、gb和hp满足电或热负荷需求，若是，执行步骤2-1，若否，执行步骤2-2；步骤2-1：将多余电或热能存储在es或hs中，并判断es或hs是否达到能量储存上限，若是，产生弃电或热；若否，结束；步骤2-2：调用es或hs放电或热；步骤3：判断是否满足电或热负荷需求，若是，结束，若否，执行步骤4；步骤4：优先切除nel或nhl，再切除pel或phl；步骤5：结束。

技术总结
本发明公开了一种考虑随机停电供能可靠性的综合能源微网优化控制方法，所述方法包括：建立EPG随机停电导致的IEMG孤岛模式的概率模型，采用场景树生成若干孤岛场景；建立IEMG多目标优化控制模型，通过对IEMG进行协同优化控制，得到最优的长期优化控制方案和两种运行模式下的短期运行方案，用于提升IEMG并网模式的运行经济性及孤岛模式的供能可靠性。本发明解决了目前研究忽略EPG停电对IEMG运行的影响，导致IEMG控制方案不合理，使IEMG出现能量供给不足，削弱IEMG供能可靠性等问题。削弱IEMG供能可靠性等问题。削弱IEMG供能可靠性等问题。


技术研发人员：穆云飞 曹严 贾宏杰 王慧媛 陈盛 王新迎 蒲天骄 武国良 祖光鑫
受保护的技术使用者：中国电力科学研究院有限公司 国网黑龙江省电力有限公司电力科学研究院
技术研发日：2022.02.08
技术公布日：2022/5/25