1.计及需求侧响应和共享储能的多微电网综合能源系统优化调度方法，其特征在于包括以下步骤：步骤一：将共享储能系统加入到含电能交互的多微电网综合能源系统中，构成含共享储能和电力交互的多微电网综合能源系统模型；

步骤二：在含共享储能和电力交互的多微电网综合能源系统模型上加入微电网内用户的自主响应行为，考虑用户侧自主响应行为对含共享储能和电力交互的多微电网综合能源系统模型带来的影响；

所述步骤二中，用户的自主响应行为考虑可移动的电负荷、可移动的气负荷、灵活的热负荷和灵活的冷负荷四种负荷类型：可移动的电负荷是指时间可进行平移的电负荷，可移动电负荷包括可移动可中断电负荷和可移动不可中断电负荷两类；可移动可中断电负荷指的是时间可平移且使用中途可中断电负荷；可移动不可中断电负荷指的是时间可平移而使用中途不可中断电负荷；

可移动可中断电负荷包括洗衣机、洗碗机电气设备；t时刻可移动不可中断电负荷的表达形式为式(5)所示：式(5)中，LEDR，s(t)表示t时刻可移动不可中断电负荷量；n表示微电网内用户户数；np表示微电网内参与可移动不可中断电负荷的响应的用户数； 表示np户用户在t时刻使用第x种设备产生的负荷；X表示可用设备的集合； 可表示为式(6)所示：式(6)中，np表示微电网内参与可移动不可中断电负荷的响应的用户数；Px,e表示设备x在调度周期内所使用的功率；

可移动可中断电负荷考虑微电网中的电动汽车充电负荷，设一个微电网中有电动汽车的用户数为nc，其充电时长为Tsp,将其等效为用户数为 充电时长为1小时的电动汽车负荷；t时刻可转移可中断负荷模型可表示为式(7)所示：式(7)中LEDR，d(t)表示t时刻可移动可中断电负荷量；Px,e表示充电汽车在一个调度周期内的充电功率；nc’(t)表示t时刻参与充电的电动汽车户数；

可转移不可中断气负荷考虑燃气壁挂炉、燃气热水器设备；t时刻可移动不可中断气负荷的表达形式为式(8)所示：式(8)中，Lg，mob(t)表示t时刻可移动不可中断气负荷量；n表示微电网内用户户数；np1表示微电网内参与可移动不可中断气负荷的响应的用户数； 表示np1户用户在t时刻使用第g种设备产生的负荷；G表示可用设备的集合； 可表示为式(9)所示：式(9)中，np1表示微电网内参与可移动不可中断电负荷的响应的用户数；Pg,e表示设备g在调度周期内所使用的功率；

灵活的热负荷考虑热水负荷，设微电网内用户对热水水温的接受范围为[Th,min,Th,max]，用Hh,min和Hh,max分别表示最小热水负荷功率和最大热水负荷功率，其表达式如式(10)和式(11)表示：式(10)‑(11)中，Th,in表示t时刻加入水的温度；Cw和ρW分别表示水的比热容和水的密度；Vc(t)表示在t时刻加入冷水的体积；Δt表示时间步长；Hh(t)表示t时刻的热负荷；灵活的热负荷满足约束条件式(12)：Hh,min(t)≤Hh(t)≤Hh,max(t)                             (12)式(12)中，Hh,min和Hh,max分别表示最小热水负荷功率和最大热水负荷功率；Hh(t)表示t时刻灵活的热负荷值；

灵活的冷负荷考虑用户对室内供冷温度的可接受范围，设用户对室内温度的可接受范围表示为[Tc,min,Tc,max]，t时刻最小制冷负荷和最大制冷负荷分别表示为式(13)和式(14)：式(13)‑(14)中，Cc,min(t)和Cc,max(t)分别表示t时刻最小制冷负荷和最大制冷负荷；

Tc,min和Tc,max分别为用户对室内温度的可接受范围上下限；Ris表示用户房屋热阻；To(t)表示t时刻用户房屋室外温度；灵活的冷负荷表示为：Cc,min(t)≤Cc(t)≤Cc,max(t)  (15)

式(15)中，Cc,min(t)和Cc,max(t)分别表示t时刻最小制冷负荷和最大制冷负荷；Cc(t)表示t时刻灵活的冷负荷；

步骤三：提出基于主从‑合作博弈的两阶段优化模型，对整个优化过程及共享储能装置最佳容量进行求解，对联盟微电网贡献度的利润进行分配；

所述步骤三中，利用主从‑合作博弈的两阶段优化模型对整个优化过程及共享储能装置最佳容量进行求解，第一阶段优化模型可以引入主从博弈模型对其进行求解，第一阶段优化模型求出各个微电网在与各微电网用户达成主从博弈均衡的情况下的能源转换情况及用户负荷使用情况，并将各个微电网的电能最佳产销情况传输给第二阶段优化模型；

各微电网综合能源系统的第一阶段目标函数表示为各微电网利润最大化函数，如式(23)所示：式(23)中，Pr表示微电网日经营利润；Lj,e(t)表示微电网j在t时刻能源e的实际负荷量；γe(t)表示能源e在t时刻的价格；CM表示微电网运行经营成本，可由式(24)表示：CM＝Cepe+Ceql                    (24)

式(24)中，Cepe表示微电网j全天从配电网和天然气网购买电能和天然气所花成本之和，其表达式为式(25)：式(25)中，E’＝{ele,gas}； 表示t时刻微电网从外界购买能源e的功率值；ζe(t)表示t时刻外界能源e的价格，由配电网和天然气网制定；

式(24)中，Ceqi表示微电网内所有设备的日维护成本，其表达式如式(26)所示：式(24)中， 表示t时刻设备b的输出功率；vb表示设备b的损耗系数；B表示所有设备的集合；

由上述可知，由各微电网运营商组成联盟微电网并先向联盟微电网内用户公布能源价格，随后用户根据能源价格进行需求响应；在这个过程中，联盟微电网运营商为先决策方，微电网用户为后决策方，因此将联盟微电网运营商称为领导者，微电网用户称为跟随者；整个第一阶段优化模型可用两个步骤进行表示：步骤1：跟随者根据领导者公布的能源价格进行需求响应，目标函数为式(22)；将跟随者负荷使用情况反馈给领导者；

步骤2：领导者根据跟随者所反馈的负荷使用情况，目标函数为式(23)；对微电网内设备的出力情况进行调整优化；

重复上述步骤，直到联盟微电网得到最优日利润，则认为达到博弈均衡；

通过第一阶段的求解，可得到各个微电网与微电网用户达到主从博弈均衡情况下的日电能产销情况；由于第二阶段中需要对共享储能装置容量进行求解，因此需对各微电网全年电能产销情况进行分析；此外，第二阶段优化模型中需考虑各微电网之间的电能交互以及各微电网对共享储能装置的出资建设情况，由于各个微电网属于不同的运营商所有，各微电网之间既有合作的关系也有竞争的关系，同时包含着复杂的利益交互关系，因此引入合作博弈对第二阶段优化模型进行求解；第二阶段目标函数可表示为式(27):式(27)中，Ccost表示联盟微电网全年最低购电成本与共享储能年建设成本之和；M表示典型日个数；W表示典型日天数；J表示联盟微电网内微电网个数；Pj,w,m(t)表示t时刻联盟微电网向配电网购买的电功率；ζele(t)表示t时刻配电网电价；Cinv,y表示共享储能建设成本。

2.根据权利要求1所述计及需求侧响应和共享储能的多微电网综合能源系统优化调度方法，其特征在于：所述步骤一中，多个微电网使用同一个储能电站对电能进行存储或释放，这个储能电站就称为共享储能系统，储能电站的年平均投资成本和维护成本表示为式(1)：式(1)中，Cinv,y表示储能电站年平均投资费用；μs表示储能电站的容量成本，单位为：元/(kWh)；μp表示储能电站的功率成本，单位为：元/kW； 和 分别表示储能电站的最大容量和最大充放电功率；Ys表示储能电站的预期使用年数；Mess为储能电站的年维护成本。

3.根据权利要求2所述计及需求侧响应和共享储能的多微电网综合能源系统优化调度方法，其特征在于：所述步骤一中，共享储能系统根据各微电网的电能剩余情况，有充电和放电两种工作状态：若t时刻微电网j多余功率Pj,s(t)>0，则说明此时微电网j有剩余电能，向共享储能系统发出储能信号；若Pj,s(t)<0，则说明此时微电网j缺少电能，向共享储能系统发出放能信号；

共享储能系统通过收集各微电网的供需信号进行供需匹配；总需求供给表达式和共享储能的总充放电需求量表达式分别如式(2)和式(3)所示：PD(t)＝|Dsum(t)|‑|Ssum(t)|                           (3)式(2)和式(3)中，Dsum为所有微电网t时刻向共享储能系统上报的总能量需求量；Ssum为所有微电网t时刻向共享储能系统上报的总能量供给量；Pj,s(t)表示t时刻微电网j多余功率；PD(t)表示共享储能的总充放电需求量，若PD(t)>0,说明在t时刻接入共享储能电站的所有微电网总放电需求为放电，共享储能系统调控中心使用储能电站放电来满足微电网群的需求；若PD(t)＜0,说明在t时刻接入共享储能电站的所有微电网总放电需求为充电，共享储能系统调控中心使用储能电站充电来满足微电网群的需求。

4.根据权利要求3所述计及需求侧响应和共享储能的多微电网综合能源系统优化调度方法，其特征在于：所述步骤一中，储能系统的荷电状态以及相应的功率应满足一定的约束条件，如式(4)所示：abs relea

式(4)中，Soc(t)代表共享储能t时刻存储的电能；η 和η 分别表示充电效率和放电效率； 和 分别表示共享储能t时刻的充电功率和放电功率；Socmax和Socmin分别表示储能系统荷电状态的上下限； 和 表示储能电站的充电状态变量和放电状态变量； 表示共享储能装置的最大充放电功率。

5.根据权利要求1所述计及需求侧响应和共享储能的多微电网综合能源系统优化调度方法，其特征在于：所述步骤二中，微电网j内用户t时刻的实际负荷量可表示为式(16)‑(19)：

Lj,e(t)＝Lj,BE(t)+Lj,e,mob(t)   (16)Lj,g(t)＝Lj,GE(t)+Lj,g,mob(t)                       (17)Lj,h(t)＝Lj,HE(t)+Hj,h(t)                         (18)Lj,c(t)＝Lj,CE(t)+Cj,c(t)                          (19)式(16)‑(19)中，Lj,e(t)、Lj,g(t)、Lj,h(t)、Lj,c(t)分别表示t时刻微电网j内用户的实际电负荷、实际气负荷、实际热负荷和实际冷负荷；Lj,BE(t)、Lj,GE(t)、Lj,HE(t)、Lj,CE(t)分别表示微电网j内t时刻的基本电负荷、基本气负荷、基本热负荷和基本冷负荷；Lj,e,mob(t)、Lj,g,mob(t)、Hj,h(t)、Cj,c(t)分别表示微电网j内的可移动电负荷、可移动不可中断气负荷、灵活的热负荷和灵活的冷负荷。

6.根据权利要求5所述计及需求侧响应和共享储能的多微电网综合能源系统优化调度方法，其特征在于：所述步骤二中，将效用函数也加入到用户目标函数中，效用函数采用二次函数来表示，用户在每个时间段内有着最佳用能量，当用户用能偏离最佳用能量时会产生满意度损失；

微电网j内用户的效用函数及满意度损失分别如式(20)、式(21)表示：

式(20)‑(21)中，Uj,ut(t)表示t时刻用户的效用函数；Uj,loss(t)表示t时刻用户满意度损失；αj,e和βj,e均表示微电网j内用户的用能偏好系数；λj,e和θj,e表示微电网j内能源e的满意度损失参数；Lj,e(t)和Lj,e，B(t)表示微电网j在t时刻能源e的实际负荷量与基线负荷；其中E＝{ele,gas,cold,heat}，表示用户购能种类；

微电网j内用户全年整体利益诉求可表示为式(22)形式：

式(22)中，Uj,eu(t)表示微电网j内用户的全年综合效益值；Uj,ut(t)表示t时刻用户的效用函数；Uj,loss(t)表示t时刻用户满意度损失；γe(t)表示能源e在t时刻的价格，由微电网运营商制定公布；d表示天数；t表示时刻数。

7.根据权利要求1所述计及需求侧响应和共享储能的多微电网综合能源系统优化调度方法，其特征在于：所述步骤三中，联盟微电网最终年利润函数如式(28)所示：式(28)中，Pr’表示年利润函数；M表示典型日个数；W表示典型日天数；E表示能源类型集合；Lj,e,w,m(t)表示微电网j在第m种典型日第w天t时刻能源e的实际负荷量；γe(t)表示能源e在t时刻的价格； 表示微电网j在第m种典型日第w天t时刻从外界购气量；ζgas(t)表示外界天然气价格； 表示微电网j内t时刻设备b的输出功率；vb表示设备b的损耗系数；B表示所有设备的集合；Ccost表示联盟微电网全年购电成本与共享储能年建设成本之和；

联盟总额外利润为各微电网联盟后的利润减去各微电网未联盟前的利润之和，如式(29)所示：式(29)中，Uext为联盟所得额外利润；Pr’为各微电网联盟后的利润， 为各微电网未联盟前的利润之和；利润分配模型如式(30)‑(32)所示：Aj＝Cvj·Uext   (30)

式(30)‑(32)中，Aj表示微电网j所分得的额外利润；Cj表示微电网j的总贡献量；Esj(t)表示微电网j在t时刻传输给其他微电网的电能总和；Pj,s(t)表示t时刻微电网j从共享储能中获取的功率；Cvj表示微电网j在联盟中的贡献度。