微网多能源系统与电网系统的联合调度方法与流程

1.本发明属于电气自动化领域，具体涉及一种微网多能源系统与电网系统的联合调度方法。


背景技术：

2.随着经济技术的发展和人们生活水平的提高，电能已经成为了人们生产和生活中必不可少的二次能源，给人们的生产和生活带来了无尽的便利。因此，保障电力的稳定可靠供应，就成为了电力系统最重要的任务之一。
3.目前，人们对于电能的需求越来越大。能源短缺与环境问题也由此而生。在此背景下，通过将多种性质各异的能源有机结合在一起，对其生产和转换环节进行协调优化，从而实现能源梯级高效利用和可再生能源大幅消纳的多能源系统(integrated energy system，ies)应运而生。微网多能源系统中存在较多的灵活资源，例如：储能、可中断电负荷、热水负荷等等，微网多能源系统可以通过调节灵活资源的应用从而参与到电网系统的调节中来。然而，单个多能源系统的可调节容量有限，由多个多能源系统参与电网辅助服务，可以提高多能源系统群的可响应容量和调节效率。
4.近年来，研究人员针对微网多能源系统的优化运行做出了许多研究：一方面，综合需求侧响应被提出，并用于改变自身的能源转换方式来调整负荷需求，从而优化微网多能源系统的经济成本和电力负荷曲线的峰谷差；另一方面，非合作博弈论用于实现综合能源产能基地、系统管理商和用户等多个存在利益冲突关系的主体的多目标优化运行，形成多方利益主体的相互制约。由此可见，目前的研究多关注于微网多能源系统内部的协调优化运行以及有利益冲突关系的多主体优化运行，但是并不存在针对微网多能源系统参与电网辅助服务的联合优化运行的研究。


技术实现要素：

5.本发明的目的在于提供一种可靠性高、实用性好且科学合理的微网多能源系统与电网系统的联合调度方法。
6.本发明提供的这种微网多能源系统与电网系统的联合调度方法，包括如下步骤：
7.s1.获取电网系统的数据信息；
8.s2.获取各个微网多能源系统的调度数据和对应的约束条件；
9.s3.根据步骤s1和步骤s2获取的数据信息，以各个微网多能源系统的响应容量和各个微网多能源系统的运行成本为双目标，构建联合调度模型；
10.s4.在步骤s1和步骤s2获取的数据信息基础上，对步骤s3构建的联合调度模型进行求解，得到联合调度方案；
11.s5.按照步骤s4得到的联合调度方案，对若干个微网多能源系统与电网进行联合调度。
12.步骤s2所述的获取各个微网多能源系统的调度数据和对应的约束条件，包括如下
步骤：
13.a.采用如下算式计算第i个微网多能源系统的可转移电负荷：
[0014][0015][0016]
式中为第i个微网多能源系统的可转移负荷m在时段t的转移容量；为第i个微网多能源系统的负荷功率可转移阈值系数；为第i个微网多能源系统的可转移电负荷m在时段t参与可转移调节的状态变量，且表示第i个微网多能源系统的可转移负荷m在时段t处于可转移工作状态且参与柔性调节，否则可转移负荷m在时段t处于可转移工作状态且参与柔性调节，否则为第i个微网多能源系统的可转移负荷m在t时段的日前预测功率；为第i个微网多能源系统的可转移负荷m的可转移工作时段集合；为第i个微网多能源系统的可转移负荷m在时段t的移出功率；为第i个微网多能源系统的可转移负荷m在时段t的移出时长；为第i个微网多能源系统的可转移负荷m在时段t的移入功率；为第i个微网多能源系统的可转移负荷m在时段t的移入时长；t为时间段t的取值；为第i个微网多能源系统的可转移负荷m在调节周期内的最大转移次数；为第i个微网多能源系统在时段t的可转移电负荷；为第i个微网多能源系统在时段t的日前预测总功率；j1为时段可移入的负荷种类数；k1为时段可移出的负荷种类数；
[0017]
b.采用如下算式计算第i个微网多能源系统的可削减电负荷：
[0018]
[0019][0020]
式中为第i个微网多能源系统的可削减电负荷n在时段t的削减容量；为第i个微网多能源系统的可削减电负荷n在时段t参与可削减调节的状态变量，且表示第i个微网多能源系统的可削减电负荷n在时段t参与可削减调节，否则表示第i个微网多能源系统的可削减电负荷n在时段t参与可削减调节，否则为第i个微网多能源系统的可削减电负荷n的功率可削减阈值系数；为第i个微网多能源系统的可削减电负荷n在t时刻的原负荷功率；为第i个微网多能源系统的可削减电负荷n的可削减工作时段集合；为第i个微网多能源系统的可削减电负荷n在时段t的削减比例系数；为第i个微网多能源系统的可削减电负荷n允许的削减速率系数；为第i个微网多能源系统的可削减电负荷n的参与可削减调节的最短时长；为第i个微网多能源系统的可削减电负荷n的参与可削减调节的时长；为第i个微网多能源系统的可削减电负荷n的参与可削减调节的最长时长；为第i个微网多能源系统在t时刻削减之后的电负荷总功率；为第i个微网多能源系统在t时刻的可削减电负荷的原总功率；k2为可削减负荷的个数；
[0021]
c.采用如下算式计算第i个微网多能源系统的柔性热负荷：
[0022][0023]
式中h
hw,min
为第i个微网多能源系统在t时刻的最小热水负荷功率；c
w
为水的比热容；p
w
为水的密度；为第i个微网多能源系统在t时刻的新加入的冷水体积；t
h,min
为用户可接受的最小水温；t
ini
为初始水温；δt为时间步长；h
hw,max
为第i个微网多能源系统在t时刻的最大热水负荷功率；t
h,max
为用户可接受的最大水温；为第i个微网多能源系统在t时刻的热水负荷功率；
[0024]
d.采用如下算式计算第i个微网多能源系统的柔性冷负荷：
[0025]
[0026]
式中c
air,min
为最小制冷负荷；为第i个微网多能源系统的t时刻的室外温度；t
air,min
为第i个微网多能源系统的满足用户要求的最低温度；r为建筑热阻；c
air,max
为最大制冷负荷；t
air,max
为第i个微网多能源系统的满足用户要求的最高温度；为第i个微网多能源系统的t时刻的冷却功率；
[0027]
e.采用如下算式作为第i个微网多能源系统的功率平衡约束：
[0028][0029][0030][0031]
式中为第i个微网多能源系统在t时刻的不可变电负荷；q
i,t
为第i个微网多能源系统在t时刻的调度响应容量；为第i个微网多能源系统在t时刻的电网购电量；为第i个微网多能源系统在t时刻的电储能充能量；为第i个微网多能源系统的产能设备y在t时刻的耗电量；y为产能设备的总数量；为第i个微网多能源系统中产电设备l在t时刻的产电量；l为微网i中所有产电设备；为第i个微网多能源系统在t时刻的电储能放能量；为第i个微网多能源系统在t时刻的不可变热负荷；为第i个微网多能源系统在t时刻的热储能充能量；为第i个微网多能源系统在t时刻的产热量；h为微网i中所有产热设备；为第i个微网多能源系统在t时刻的热储能放能量；为第i个微网多能源系统在t时刻的不可变冷负荷；为第i个微网多能源系统在t时刻的冷储能充能量；为第i个微网多能源系统在t时刻的产冷量；c为微网i中所有产冷设备；为第i个微网多能源系统在t时刻的冷储能放能量；
[0032]
f.采用如下算式作为第i个微网多能源系统的储能设备约束：
[0033][0034]
式中e
i,x,t
为第i个微网多能源系统在t时刻的x储能装置的当前容量；β
x
为x储能装置自损耗系数；为第i个微网多能源系统在t时刻的x储能装置的充能功率；为第i个微网多能源系统的x储能装置的充能效率；为第i个微网多能源系统在t时刻的x储能装置的放能效率；为第i个微网多能源系统的x储能装置的放能效率；为第i个微网多能源系统在t时刻的x储能装置的充能状态；为第i个微网多能源系统的x储能装置的最大充能功率；为第i个微网多能源系统在t时刻的x储能装置的放能状态；为第i个微网多能源系统的x储能装置的最大放能功率；e
i,x,min
为第i个微网多能源系统的x储能装置的最小储能容量；e
i,x,max
为第i个微网多能源系统的x储能装置的最大储能容量；x为指示变量，取值为1、2或3，x＝1表示电储能装置，x＝2表示热储能装置，x＝3表示冷储能装置；
[0035]
g.采用如下算式作为第i个微网多能源系统的设备出力约束：
[0036][0037]
式中为第i个微网多能源系统的设备q在t时刻的输出功率；η
i,q
为第i个微网多能源系统的设备q的效用系数；为第i个微网多能源系统的设备q在t时刻的输入功率；为第i个微网多能源系统的设备q的额定功率。
[0038]
步骤s3所述的根据步骤s1和步骤s2获取的数据信息，以各个微网多能源系统的响应容量和各个微网多能源系统的运行成本为双目标，构建联合调度模型，具体包括如下步骤：
[0039]
a.采用如下算式作为各个微网多能源系统的响应容量最大化目标函数：
[0040][0041]
式中t
s
为微网多能源系统参与调度的时间长度；start为微网多能源系统参与调度的初始时间；m为参与调度的微网多能源系统的个数；q
i,t
为第i个微网多能源系统在t时刻的调度响应容量；
[0042]
b.采用如下算式作为各个微网多能源系统的运行成本最小化目标函数：
[0043][0044]
式中为t时刻的电网购电价格；为第i个微网多能源系统在t时刻的电网购电量；c
gas
为购气价格；为第i个微网多能源系统在t时刻的购气量；r
i
为第i个微网多能源系统参与调度所获得的收益。
[0045]
步骤s4所述的在步骤s1和步骤s2获取的数据信息基础上，对步骤s3构建的联合调度模型进行求解，得到联合调度方案，具体为采用混沌粒子群+gourbi算法进行求解。
[0046]
步骤s4所述的在步骤s1和步骤s2获取的数据信息基础上，对步骤s3构建的联合调度模型进行求解，得到联合调度方案，具体包括如下步骤：
[0047]
(1)初始化混沌粒子群算法的参数，生成各个微网多能源系统的响应容量；
[0048]
(2)计算各个微网多能源系统参与调度的收益分配；
[0049]
(3)根据步骤(1)生成的各个微网多能源系统的响应容量，优化各个微网多能源系统的运行调度，得到当前轮次的各个微网多能源系统的最小运行成本；
[0050]
(4)判断步骤(3)得到的当前轮次的各个微网多能源系统的最小运行成本是否优于之前轮次得到的各个微网多能源系统的最小运行成本：
[0051]
若当前轮次的各个微网多能源系统的最小运行成本不优于之前轮次得到的各个微网多能源系统的最小运行成本，则返回步骤(2)重新进行优化；
[0052]
若当前轮次的各个微网多能源系统的最小运行成本优于之前轮次得到的各个微网多能源系统的最小运行成本，则进行后续步骤；
[0053]
(5)更新各个微网多能源系统的运行成本；
[0054]
(6)判断是否达到设定的迭代次数上限：
[0055]
若达到设定的迭代次数上限，则输出最终的各个微网多能源系统的运行策略和各个微网多能源系统的响应容量，并作为联合调度方案；结束；
[0056]
若未达到设定的迭代次数上限，则返回步骤(2)再次进行优化。
[0057]
步骤s4所述的在步骤s1和步骤s2获取的数据信息基础上，对步骤s3构建的联合调度模型进行求解，得到联合调度方案，具体包括如下步骤：
[0058]
(1)初始化混沌粒子群算法的参数，生成各个微网多能源系统的响应容量；
[0059]
(2)计算各个微网多能源系统参与调度的收益分配；
[0060]
(3)根据步骤(1)生成的各个微网多能源系统的响应容量，优化各个微网多能源系统的运行调度，得到当前轮次的各个微网多能源系统的最小运行成本；
[0061]
(4)判断步骤(3)得到的当前轮次的各个微网多能源系统的最小运行成本是否优于之前轮次得到的各个微网多能源系统的最小运行成本：
[0062]
若当前轮次的各个微网多能源系统的最小运行成本不优于之前轮次得到的各个微网多能源系统的最小运行成本，则返回步骤(2)重新进行优化；
[0063]
若当前轮次的各个微网多能源系统的最小运行成本优于之前轮次得到的各个微网多能源系统的最小运行成本，则进行后续步骤；
[0064]
(5)更新各个微网多能源系统的运行成本；
[0065]
(6)判断是否达到设定的迭代次数上限：
[0066]
若达到设定的迭代次数上限，则输出最终的各个微网多能源系统的运行策略和各个微网多能源系统的响应容量，并作为联合调度方案；结束；
[0067]
若未达到设定的迭代次数上限，则返回步骤(2)再次进行优化；
[0068]
(7)将步骤(6)得到的优化的各个微网多能源系统的运行策略和各个微网多能源系统的响应容量，发送至电网系统和对应的各个微网多能源系统；
[0069]
(8)电网系统和各个微网多能源系统对接收到的各个微网多能源系统的运行策略和各个微网多能源系统的响应容量进行审核：
[0070]
若审核通过，则将审核通过的各个微网多能源系统的运行策略和各个微网多能源系统的响应容量，并作为联合调度方案。
[0071]
若审核不通过，则返回步骤(1)再次进行优化，直至审核通过，得到最终的联合调度方案。
[0072]
所述的计算各个微网多能源系统参与调度的收益分配，具体为采用shapley值法计算各个微网多能源系统参与调度的收益。
[0073]
所述的计算各个微网多能源系统参与调度的收益分配，具体包括如下步骤：
[0074]
采用如下算式计算第i个微网多能源系统参与调度的收益r
i
：
[0075][0076]
式中s
i
为包含微网多能源系统i的所有子集；|s|为子集s中的微网多能源系统的个数；m为微网多能源系统的数量；r
′
(s)为集合s产生的收益；r
′
(s/i)为微网集合s除去微网多能源系统i后产生的收益。
[0077]
本发明提供的这种微网多能源系统与电网系统的联合调度方法，通过对若干个微网多能源系统进行联合调度和综合运行，不仅实现了多个微网多能源系统参与电网系统的联合调度，而且提高了电网的运行可靠性，减少了微网多能源系统的运行成本，而且可靠性高、实用性好且科学合理。
附图说明
[0078]
图1为本发明的方法流程示意图。
具体实施方式
[0079]
如图1所示为本发明的方法流程示意图：本发明提供的这种微网多能源系统与电网系统的联合调度方法，包括如下步骤：
[0080]
s1.获取电网系统的数据信息；
[0081]
s2.获取各个微网多能源系统的调度数据和对应的约束条件；具体包括如下步骤：
[0082]
a.采用如下算式计算第i个微网多能源系统的可转移电负荷：
[0083][0084][0085]
式中为第i个微网多能源系统的可转移负荷m在时段t的转移容量；为第i个微网多能源系统的负荷功率可转移阈值系数，取值范围为0～100％；为第i个微网多能源系统的可转移电负荷m在时段t参与可转移调节的状态变量，且表示第i个微网多能源系统的可转移负荷m在时段t处于可转移工作状态且参与柔性调节，否则个微网多能源系统的可转移负荷m在时段t处于可转移工作状态且参与柔性调节，否则为第i个微网多能源系统的可转移负荷m在t时段的日前预测功率；为第i个微网多能源系统的可转移负荷m的可转移工作时段集合；为第i个微网多能源系统的可转移负荷m在时段t的移出功率；为第i个微网多能源系统的可转移负荷m在时段t的移出时长；为第i个微网多能源系统的可转移负荷m在时段t的移入功率；为第i个微网多能源系统的可转移负荷m在时段t的移入时长；t为时间段t的取值；为第i个微网多能源系统的可转移负荷m在调节周期内的最大转移次数；为第i个微网多能源系统在时段t的可转移电负荷；为第i个微网多能源系统在时段t的日前预测总功率；j1为时段可移入的负荷种类数；k1为时段可移出的负荷种类数；对于可转移负荷，移入
‑
移出的用电量应保持不变；
[0086]
b.采用如下算式计算第i个微网多能源系统的可削减电负荷：
[0087][0088]
[0089]
式中为第i个微网多能源系统的可削减电负荷n在时段t的削减容量；为第i个微网多能源系统的可削减电负荷n在时段t参与可削减调节的状态变量，且表示第i个微网多能源系统的可削减电负荷n在时段t参与可削减调节，否则表示第i个微网多能源系统的可削减电负荷n在时段t参与可削减调节，否则为第i个微网多能源系统的可削减电负荷n的功率可削减阈值系数，取值范围为0～100％，当时，负荷n等效于可中断负荷；为第i个微网多能源系统的可削减电负荷n在t时刻的原负荷功率；为第i个微网多能源系统的可削减电负荷n的可削减工作时段集合；为第i个微网多能源系统的可削减电负荷n在时段t的削减比例系数；为第i个微网多能源系统的可削减电负荷n允许的削减速率系数；为第i个微网多能源系统的可削减电负荷n的参与可削减调节的最短时长；为第i个微网多能源系统的可削减电负荷n的参与可削减调节的时长；为第i个微网多能源系统的可削减电负荷n的参与可削减调节的最长时长；为第i个微网多能源系统在t时刻削减之后的电负荷总功率；为第i个微网多能源系统在t时刻的可削减电负荷的原总功率；k2为可削减负荷的个数；
[0090]
c.采用如下算式计算第i个微网多能源系统的柔性热负荷：
[0091][0092]
式中h
hw,min
为第i个微网多能源系统在t时刻的最小热水负荷功率；c
w
为水的比热容；p
w
为水的密度；为第i个微网多能源系统在t时刻的新加入的冷水体积；t
h,min
为用户可接受的最小水温；t
ini
为初始水温；δt为时间步长；h
hw,max
为第i个微网多能源系统在t时刻的最大热水负荷功率；t
h,max
为用户可接受的最大水温；为第i个微网多能源系统在t时刻的热水负荷功率；
[0093]
d.采用如下算式计算第i个微网多能源系统的柔性冷负荷：
[0094][0095]
式中c
air,min
为最小制冷负荷；为第i个微网多能源系统的t时刻的室外温度；t
air,min
为第i个微网多能源系统的满足用户要求的最低温度；r为建筑热阻；c
air,max
为最大制
冷负荷；t
air,max
为第i个微网多能源系统的满足用户要求的最高温度；为第i个微网多能源系统的t时刻的冷却功率；
[0096]
e.采用如下算式作为第i个微网多能源系统的功率平衡约束：
[0097][0098][0099][0100]
式中为第i个微网多能源系统在t时刻的不可变电负荷；q
i,t
为第i个微网多能源系统在t时刻的调度响应容量；为第i个微网多能源系统在t时刻的电网购电量；为第i个微网多能源系统在t时刻的电储能充能量；为第i个微网多能源系统的产能设备y在t时刻的耗电量；y为产能设备的总数量；为第i个微网多能源系统中产电设备l在t时刻的产电量；l为微网i中所有产电设备；为第i个微网多能源系统在t时刻的电储能放能量；为第i个微网多能源系统在t时刻的不可变热负荷；为第i个微网多能源系统在t时刻的热储能充能量；为第i个微网多能源系统在t时刻的产热量；h为微网i中所有产热设备；为第i个微网多能源系统在t时刻的热储能放能量；为第i个微网多能源系统在t时刻的不可变冷负荷；为第i个微网多能源系统在t时刻的冷储能充能量；为第i个微网多能源系统在t时刻的产冷量；c为微网i中所有产冷设备；为第i个微网多能源系统在t时刻的冷储能放能量；
[0101]
f.采用如下算式作为第i个微网多能源系统的储能设备约束：
[0102][0103]
式中e
i,x,t
为第i个微网多能源系统在t时刻的x储能装置的当前容量；β
x
为x储能装
置自损耗系数；为第i个微网多能源系统在t时刻的x储能装置的充能功率；为第i个微网多能源系统的x储能装置的充能效率；为第i个微网多能源系统在t时刻的x储能装置的放能效率；为第i个微网多能源系统的x储能装置的放能效率；为第i个微网多能源系统在t时刻的x储能装置的充能状态；为第i个微网多能源系统的x储能装置的最大充能功率；为第i个微网多能源系统在t时刻的x储能装置的放能状态；为第i个微网多能源系统的x储能装置的最大放能功率；e
i,x,min
为第i个微网多能源系统的x储能装置的最小储能容量；e
i,x,max
为第i个微网多能源系统的x储能装置的最大储能容量；x为指示变量，取值为1、2或3，x＝1表示电储能装置，x＝2表示热储能装置，x＝3表示冷储能装置；
[0104]
g.采用如下算式作为第i个微网多能源系统的设备出力约束：
[0105][0106]
式中为第i个微网多能源系统的设备q在t时刻的输出功率；η
i,q
为第i个微网多能源系统的设备q的效用系数；对于离心机和热泵等机组，其效用系数均大于1；对于热水锅炉和燃气轮机等机组，其效用系数小于1；为第i个微网多能源系统的设备q在t时刻的输入功率；为第i个微网多能源系统的设备q的额定功率；
[0107]
s3.根据步骤s1和步骤s2获取的数据信息，以各个微网多能源系统的响应容量和各个微网多能源系统的运行成本为双目标，构建联合调度模型；具体包括如下步骤：
[0108]
a.采用如下算式作为各个微网多能源系统的响应容量最大化目标函数：
[0109][0110]
式中t
s
为微网多能源系统参与调度的时间长度；start为微网多能源系统参与调度的初始时间；m为参与调度的微网多能源系统的个数；q
i,t
为第i个微网多能源系统在t时刻的调度响应容量；
[0111]
b.采用如下算式作为各个微网多能源系统的运行成本最小化目标函数：
[0112][0113]
式中为t时刻的电网购电价格；为第i个微网多能源系统在t时刻的电网购电量；c
gas
为购气价格；为第i个微网多能源系统在t时刻的购气量；r
i
为第i个微网多能源系统参与调度所获得的收益；
[0114]
s4.在步骤s1和步骤s2获取的数据信息基础上，对步骤s3构建的联合调度模型进行求解，得到联合调度方案；具体为采用混沌粒子群+gourbi算法进行求解，得到联合调度方案；
[0115]
具体实施时，本步骤可以有两种实施方案：
[0116]
方案一：具体包括如下步骤：
[0117]
(1)初始化混沌粒子群算法的参数，生成各个微网多能源系统的响应容量；
[0118]
(2)计算各个微网多能源系统参与调度的收益分配；具体为采用shapley值法计算各个微网多能源系统参与调度的收益；具体包括如下步骤：
[0119]
采用如下算式计算第i个微网多能源系统参与调度的收益r
i
：
[0120][0121]
式中s
i
为包含微网多能源系统i的所有子集；|s|为子集s中的微网多能源系统的个数；m为微网多能源系统的数量；r
′
(s)为集合s产生的收益；r
′
(s/i)为微网集合s除去微网多能源系统i后产生的收益；
[0122]
(3)根据步骤(1)生成的各个微网多能源系统的响应容量，优化各个微网多能源系统的运行调度，得到当前轮次的各个微网多能源系统的最小运行成本；
[0123]
(4)判断步骤(3)得到的当前轮次的各个微网多能源系统的最小运行成本是否优于之前轮次得到的各个微网多能源系统的最小运行成本：
[0124]
若当前轮次的各个微网多能源系统的最小运行成本不优于之前轮次得到的各个微网多能源系统的最小运行成本，则返回步骤(2)重新进行优化；
[0125]
若当前轮次的各个微网多能源系统的最小运行成本优于之前轮次得到的各个微网多能源系统的最小运行成本，则进行后续步骤；
[0126]
(5)更新各个微网多能源系统的运行成本；
[0127]
(6)判断是否达到设定的迭代次数上限：
[0128]
若达到设定的迭代次数上限，则输出最终的各个微网多能源系统的运行策略和各个微网多能源系统的响应容量，并作为联合调度方案；结束；
[0129]
若未达到设定的迭代次数上限，则返回步骤(2)再次进行优化；
[0130]
方案二：具体包括如下步骤：
[0131]
(1)初始化混沌粒子群算法的参数，生成各个微网多能源系统的响应容量；
[0132]
(2)计算各个微网多能源系统参与调度的收益分配；具体为采用shapley值法计算各个微网多能源系统参与调度的收益；具体包括如下步骤：
[0133]
采用如下算式计算第i个微网多能源系统参与调度的收益r
i
：
[0134][0135]
式中s
i
为包含微网多能源系统i的所有子集；|s|为子集s中的微网多能源系统的个数；m为微网多能源系统的数量；r
′
(s)为集合s产生的收益；r
′
(s/i)为微网集合s除去微网多能源系统i后产生的收益；
[0136]
(3)根据步骤(1)生成的各个微网多能源系统的响应容量，优化各个微网多能源系统的运行调度，得到当前轮次的各个微网多能源系统的最小运行成本；
[0137]
(4)判断步骤(3)得到的当前轮次的各个微网多能源系统的最小运行成本是否优于之前轮次得到的各个微网多能源系统的最小运行成本：
[0138]
若当前轮次的各个微网多能源系统的最小运行成本不优于之前轮次得到的各个微网多能源系统的最小运行成本，则返回步骤(1)重新进行优化；
[0139]
若当前轮次的各个微网多能源系统的最小运行成本优于之前轮次得到的各个微网多能源系统的最小运行成本，则进行后续步骤；
[0140]
(5)更新各个微网多能源系统的运行成本；
[0141]
(6)判断是否达到设定的迭代次数上限：
[0142]
若达到设定的迭代次数上限，则输出优化的各个微网多能源系统的运行策略和各个微网多能源系统的响应容量；结束；
[0143]
若未达到设定的迭代次数上限，则返回步骤(2)再次进行优化；
[0144]
(7)将步骤(6)得到的优化的各个微网多能源系统的运行策略和各个微网多能源系统的响应容量，发送至电网系统和对应的各个微网多能源系统；
[0145]
(8)电网系统和各个微网多能源系统对接收到的各个微网多能源系统的运行策略和各个微网多能源系统的响应容量进行审核：
[0146]
若审核通过，则将审核通过的各个微网多能源系统的运行策略和各个微网多能源系统的响应容量，并作为联合调度方案；
[0147]
若审核不通过，则返回步骤(1)再次进行优化，直至审核通过，得到最终的联合调度方案；
[0148]
s5.按照步骤s4得到的联合调度方案，对若干个微网多能源系统与电网进行联合调度。