本发明涉及工业虚拟电厂调度，尤其涉及一种考虑制氢储能的工业虚拟电厂优化调度方法。
背景技术：
1、虚拟电厂（virtual power plant，vpp）作为新型电力系统的一种，通过前沿技术将分布式能源聚集在一起，实现协调优化，提高电力系统的灵活性且低成本拓展可调资源。
2、然而，虚拟电厂中工业负荷在生产过程中存在用电、用热间的耦合关系及热电联产机组具有“热电耦合”特性。尽管热电联产机组的强制输出功率保持在高水平，但由于热电耦合特性，它降低用户供热需求输电量的可调范围，降低机组的灵活性，影响虚拟电厂利用新能源的能力。
技术实现思路
1、针对现有技术的不足，本发明提供了一种考虑制氢储能的工业虚拟电厂优化调度方法，解决了目前虚拟电厂由于热电耦合特性所导致的可调范围小以及降低机组灵活性的技术问题。
2、为解决上述技术问题，本发明提供了如下技术方案：一种考虑制氢储能的工业虚拟电厂优化调度方法，该方法包括以下过程：
3、构建包含电解槽、氢燃料电池的氢储能系统模型；
4、基于氢储能系统模型并结合工业虚拟电厂的内部资源和历史数据，建立考虑电制氢工业虚拟电厂资源的约束条件；
5、在约束条件下构建工业虚拟电厂调度模型，并以年利润最大化为目标构建工业虚拟电厂调度模型的目标函数；
6、采用改进型ieee30节点测试系统作为测试场景，并将工业虚拟电厂调度模型带入改进型ieee30节点测试系统中用以模拟优化调度并生成结果。
7、进一步地，氢储能系统模型由两部分组成，分别为氢燃料电池模型和电解槽模型，氢燃料电池模型的表达式为：
8、
9、其中，表示第天的时刻hfc输出的能量；表示第天的时刻的输入hfc的氢能；表示hfc的转换效率；、分别为输入hrc氢能的上下限；、分别为hfc的爬坡上、下限；、分别为燃料电池放电功率的上下限；
10、电解槽模型的表达式为：
11、
12、其中，表示第天的时刻el输出的能量；表示第天的时刻输入el的电能；表示el的转换效率；、分别为输入el电能的上下限；、分别为el的爬坡上、下限；、分别为el充电功率的上下限。
13、进一步地，约束条件包括考虑工业负荷生产特性的需求响应约束、多能耦合特性约束、储能规划约束、热电平衡约束以及其余分布式资源约束，其中，
14、工业负荷生产特性的需求响应约束包括最大负荷削减量约束、最小负荷削减时间间隔约束、最大负荷爬坡速率约束和负荷需求响应耦合比例约束；
15、多能耦合特性约束包括多能负荷综合需求响应中负荷的用户满意度约束以及热电联产机组中的电出力约束、热出力约束、爬坡约束；
16、储能规划约束包括电储能、热储能、氢储能各自的建设容量约束以及建设投资成本约束两部分。
17、进一步地，所述最大负荷削减量约束、最小负荷削减时间间隔约束、最大负荷爬坡速率约束和负荷需求响应耦合比例约束分别为：
18、最大负荷削减量约束为：
19、
20、其中，表示电负荷的最大削减率；表示第天的时刻电负荷的需求响应削减量；
21、最小负荷削减时间间隔约束为：
22、
23、其中，表示电负荷到第天的时刻已经连续运行的时间；表示电负荷的最小负荷削减时间间隔；表示电负荷在第天的时刻的开关机状态；
24、最大负荷爬坡速率约束为：
25、
26、其中，、分别为电负荷的最大滑坡、削减爬坡速率；表示第天的时刻电负荷的需求响应削减量；
27、负荷需求响应耦合比例约束为：
28、
29、其中，表示上下游负荷间的耦合延时间隔，即上游负荷至下游负荷所需要的最小物料传递时间的总和；表示从上游负荷到下游负荷生产损耗系数；、分别表示上游负荷和下游负荷运行的0-1变量；、分别表示上游负荷和下游负荷的需求响应削减量。
30、进一步地，所述用户满意度约束、电出力约束、热出力约束、爬坡约束分别为：
31、用户满意度约束为：
32、
33、其中，表示第天的时刻电负荷削减量；表示电负荷再单位时间内最大负荷削减率；表示电负荷在相邻两个时间段内连续削减率之和的最大值；表示第天的时刻电负荷的需求响应削减量；
34、电出力约束为：
35、
36、其中，、分别表示第个机组的电输出功率的最小值和最大值；和分别表示机组在最小凝气运行和最小进气运行时电功率与热功率之间的电热系数；和分别表示第天的时刻第个机组的电、热输出功率；
37、热出力约束为：
38、
39、其中，、为第个机组的热输出功率的最小值和最大值；表示第天的时刻第个机组的热输出功率；
40、爬坡约束为：
41、
42、其中，、分别表示第个机组的滑坡和爬坡率；表示第天的时刻第个机组的电输出功率。
43、进一步地，所述电储能、热储能、氢储能的规划约束相同，则电储能的规划约束为：
44、
45、其中，为第个电储能的最小储存容量；表示第天的时刻第个电储能的实际储能；表示电储能的建设容量；为电储能的最大建设容量；表示电储能单位容量的建设成本；表示电储能的最大建设投资成本。
46、进一步地，以年利润最大化为目标构建工业虚拟电厂调度模型的目标函数，具体包括以下过程：
47、计算虚拟电厂中热电联产机组运行成本、常规机组运行成本、多能负荷综合需求响应成本以及电解槽运行成本；
48、基于热电联产机组运行成本、常规机组运行成本、多能负荷综合需求响应成本以及电解槽运行成本构建工业虚拟电厂调度模型的目标函数，即：
49、
50、其中，表示虚拟电厂的年利润；表示一年的天数；表示虚拟电厂向电力市场售电收益，表示第天的时刻的市场电价，表示第天的时刻虚拟电厂在电力市场中的交易电量，卖电为正，购电为负；表示虚拟电厂电负荷的售电收益，表示第天的时刻的电负荷销售价格，表示第天的时刻的电负荷量；表示虚拟电厂热负荷的售电收益，表示第天的时刻的热负荷销售价格，表示第天的时刻的热负荷量；表示虚拟电厂氢燃料电池的售电收益，表示第天的时刻的氢燃料电池发电的销售价格，表示第天的时刻氢燃料电池的放电功率；表示第天的热电联产机组运行成本；表示第天的常规机组运行成本；表示第天的多能负荷综合需求响应成本；表示第天的电解槽机组成本；表示电储能成本；表示热储能成本；表示氢储能成本。
51、进一步地，热电联产机组运行成本、常规机组运行成本、多能负荷综合需求响应成本以及电解槽运行成本的计算公式分别为：
52、热电联产机组运行成本为：
53、
54、其中，是热电联产机组的单元数目；、、、、、表示机组成本系数；、分别表示第天的时刻第个机组的电、热输出功率；
55、常规机组运行成本为：
56、
57、其中，是常规机组的单元数目；，，表示机组成本系数；表示第天的时刻第个机组输出功率；
58、多能负荷综合需求响应成本为：
59、
60、其中，、分别表示第天电、热负荷需求响应成本；
61、电解槽运行成本为：
62、
63、其中，为电解槽单位数目；表示电解槽成本价格；表示第天的时刻第个电解槽的充电功率。
64、借由上述技术方案，本发明提供了一种考虑制氢储能的工业虚拟电厂优化调度方法，至少具备以下有益效果：
65、本发明所提出的优化调度方法能够降低工业虚拟电厂多能耦合的程度，增强热电联产机组的调峰能力，还能解决风电机组的消纳问题；
66、本发明通过对原有的ieee30节点测试系统进行改进，以适应工业虚拟电厂的需求和加入新型的能源技术，确保在复杂的多能耦合场景下实现高效和准确的运行；
67、本发明通过在约束条件下并结合目标函数搭建考虑制氢储能的工业虚拟电厂优化调度模型，以解决目前虚拟电厂由于热电耦合特性所导致的可调范围小以及降低机组灵活性的问题。
1.一种考虑制氢储能的工业虚拟电厂优化调度方法，其特征在于，该方法包括以下过程：
2.根据权利要求1所述的工业虚拟电厂优化调度方法，其特征在于，所述氢储能系统模型由两部分组成，分别为氢燃料电池模型和电解槽模型，氢燃料电池模型的表达式为：
3.根据权利要求1所述的工业虚拟电厂优化调度方法，其特征在于，所述约束条件包括考虑工业负荷生产特性的需求响应约束、多能耦合特性约束、储能规划约束、热电平衡约束以及其余分布式资源约束，其中，
4.根据权利要求3所述的工业虚拟电厂优化调度方法，其特征在于，所述最大负荷削减量约束、最小负荷削减时间间隔约束、最大负荷爬坡速率约束和负荷需求响应耦合比例约束分别为：
5.根据权利要求3所述的工业虚拟电厂优化调度方法，其特征在于，所述用户满意度约束、电出力约束、热出力约束、爬坡约束分别为：
6.根据权利要求3所述的工业虚拟电厂优化调度方法，其特征在于，所述电储能、热储能、氢储能的规划约束相同，则电储能的规划约束为：
7.根据权利要求1所述的工业虚拟电厂优化调度方法，其特征在于，所述以年利润最大化为目标构建工业虚拟电厂调度模型的目标函数，具体包括以下过程：
8.根据权利要求7所述的工业虚拟电厂优化调度方法，其特征在于，在所述热电联产机组运行成本、常规机组运行成本、多能负荷综合需求响应成本以及电解槽运行成本的计算公式分别为：
9.一种计算机设备，其特征在于，包括处理器和存储器，所述存储器用于存储计算机程序，所述计算机程序被所述处理器执行时实现如权利要求1至8任一项所述的工业虚拟电厂优化调度方法。
