1.本发明涉及多能流系统调度技术领域，尤其涉及一种考虑柔性氢需求的电氢能源系统调度方法。
背景技术：2.电氢能源系统是一种基于电转气技术的新型能源系统，通过电转气技术将多余的电能转化成氢能并进行储存，解决了电能不能长期、大规模存储的痛点。电氢能源系统的这些特点为季节性显著、间歇性突出、不确定性强的可再生能源大规模接入电力系统提供了技术支撑，因此在工程应用上具有光明的前景，正得到越来越多的关注。
3.电氢能源系统的经济运行取决于电、氢能量的协同优化。一方面，电
‑
氢耦合影响使系统运行问题变得十分复杂。另一方面，系统中电、氢负荷具有一定的柔性，即可在一定时段内对负荷进行转移，为系统运行带来了额外灵活性的同时，也使得电氢能源系统的调度决策成为一个极具挑战的难题。如果不对电、氢能流进行有效的协同调度，一方面可能导致调度计划在物理上是不可行的，给系统的安全稳定运行造成威胁，另一方面将影响系统运行的经济性。
4.鉴于上述问题，本设计人基于从事此类产品工程应用多年丰富的实务经验及专业知识，并配合学理的运用，以期设计考虑柔性氢需求的电氢能源系统调度方法，保证系统安全、高效、经济的运行。
技术实现要素：5.本发明的目的是：提供一种考虑柔性氢需求的电氢能源系统调度方法，能够有效解决电、氢能流协同问题，同时考虑了电、氢负荷的柔性，进一步为系统运行提供额外的灵活性，通过所提出的调度方法获得经济可行的电氢综合能源系统调度计划，保证系统安全、高效、经济的运行。
6.为了达到上述目的，本发明提供一种考虑柔性氢需求的电氢能源系统调度方法，包括：
7.根据电力系统的运行参数，建立电负荷柔性方程、购电售电约束方程和可再生能源出力约束方程；
8.根据电制氢系统的运行参数，建立氢负荷柔性方程和电制氢安全运行约束方程；
9.利用电力系统和电制氢系统的电功率平衡关系，建立电功率平衡约束方程；
10.以调度周期内电氢能源系统的运行成本最低为目标函数，以所述电负荷柔性方程、购电售电约束方程、可再生能源出力约束方程、氢负荷柔性方程、电制氢安全运行约束方程和电功率平衡约束方程为约束条件，建立电氢能源系统调度模型；
11.求解所述电氢能源系统调度模型，得到优化调度结果。
12.优选地，所述电力系统的运行参数包括电力系统的网络参数、拓扑、电负荷数据和可再生能源出力数据，所述电制氢系统的运行参数包括电制氢系统的网络参数、拓扑、氢负
荷数据、设备类型和设备容量。
13.优选地，所述目标函数的表达式为：
[0014][0015]
式中，t代表调度时段，t代表调度周期，代表电力系统在t时段向电网购电的单位成本，代表电力系统在t时段向电网购电的功率，代表电力系统在t时段向电网售电的单位成本，p
tgrid
‑
代表电力系统在t时段向电网售电的功率，c
shift
代表电力系统电负荷转移需要向用户支付的单位惩罚成本，p
te,in
代表电力系统在t时段从其它时段转入的电负荷。
[0016]
优选地，所述电负荷柔性方程的表达式为：
[0017][0018]
式中，p
te
代表电力系统在t时段的总电负荷，p
te,fix
代表电力系统在t时段的固定电负荷，p
te,out
代表电力系统在t时段转出到其他时段的电负荷，p
te,in,max
代表电力系统在t时段从其它时段可转入的电负荷的最大值，p
te,out,max
代表电力系统在t时段可转出到其他时段的电负荷的最大值。
[0019]
优选地，所述购电售电约束方程的表达式为：
[0020][0021]
式中，代表电力系统在t时段向电网购电时的变量，代表电力系统在t时段向电网售电时的变量，代表电力系统与电网所允许交换电功率的最大值。
[0022]
优选地，所述可再生能源出力约束方程的表达式为：
[0023][0024]
式中，代表电力系统在t时段的光伏出力调度值，代表电力系统在t时段的光伏出力预测值，σ
t,fore
代表电力系统在t时段的光伏出力预测值的标准差，代表标准正态分布n(0,1)的逆累计分布函数，η代表光伏预测值的置信水平。
[0025]
优选地，所述氢负荷柔性方程的表达式为：
[0026][0027]
式中，t
h
代表电制氢系统的氢负荷所需的供需平衡周期长度，k代表电制氢系统的氢负荷的供需平衡周期序列，p
th
代表电制氢系统在t时段用户的氢负荷；代表电制氢系统在第k个供需平衡周期内用户所需的氢需求。
[0028]
优选地，所述电制氢安全运行约束方程的表达式为：
[0029][0030]
式中，η
p2h
代表电制氢设备的转化效率，p
tp2h
代表在t时段电制氢设备消耗的电功率，λ代表电制氢设备的最小负载水平，代表在t时段电制氢设备的运行状态，c
p2h
代表电制氢设备的容量。
[0031]
优选地，所述电功率平衡约束方程的表达式为：
[0032]
p
tgrid+
‑
p
tgrid
‑
+p
tpv
‑
p
tp2h
＝p
te
。
[0033]
优选地，所述求解所述电氢能源系统调度模型，得到优化调度结果，包括：
[0034]
利用混合整数线性规划法求解所述电氢能源系统调度模型，得到优化调度结果。
[0035]
本发明至少具有以下有益效果：
[0036]
本发明通过建立电负荷柔性方程、购电售电约束方程、可再生能源出力约束方程、氢负荷柔性方程、电制氢安全运行约束方程和电功率平衡约束方程，并以调度周期内电氢能源系统的运行成本最低为目标函数，构建电氢能源系统调度模型，最后求解该模型可以得到电氢能源系统的优化调度结果，由于充分考虑了电负荷和氢负荷的柔性，可以为电氢综合能源系统运行提供额外的灵活性，以提高系统运行经济性，与现有技术相比，能够有效解决电、氢能流协同问题，同时考虑了电、氢负荷的柔性，进一步为系统运行提供额外的灵活性，通过所提出的调度方法获得经济可行的电氢综合能源系统调度计划，保证系统安全、高效、经济的运行。
附图说明
[0037]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0038]
图1为本发明实施例中考虑柔性氢需求的电氢能源系统调度方法的流程示意图；
[0039]
图2为本发明实施例中典型电氢能源系统中的电负荷和单位光伏出力的示意图；
[0040]
图3为本发明实施例中采用考虑柔性氢需求的电氢能源系统调度方法在不同氢负荷柔性下系统经济性的对比图；
[0041]
图4为本发明实施例中采用考虑柔性氢需求的电氢能源系统调度方法在不同氢负
荷柔性下电制氢设备输入电功率的示意图；
[0042]
图5为本发明实施例中采用考虑柔性氢需求的电氢能源系统调度方法在不同氢负荷柔性下系统电负荷转移的示意图。
具体实施方式
[0043]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。
[0044]
需要说明的是，当元件被称为“固定于”另一个元件，它可以直接在另一个元件上或者也可以存在居中的元件。当一个元件被认为是“连接”另一个元件，它可以是直接连接到另一个元件或者可能同时存在居中元件。本文所使用的术语“垂直的”、“水平的”、“左”、“右”以及类似的表述只是为了说明的目的，并不表示是唯一的实施方式。
[0045]
除非另有定义，本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。在本发明的说明书中所使用的术语只是为了描述具体的实施例的目的，不是旨在于限制本发明。本文所使用的术语“及/或”包括一个或多个相关的所列项目的任意的和所有的组合。
[0046]
请参阅图1，本发明实施例提供一种考虑柔性氢需求的电氢能源系统调度方法，包括：
[0047]
s110、根据电力系统的运行参数，建立电负荷柔性方程、购电售电约束方程和可再生能源出力约束方程。
[0048]
本发明实施例中，电氢能源系统中的电负荷具有一定的柔性，可以在一定时段内对电负荷进行转移，将当前时段的电负荷转出到其他时段，或者从其他时段转入电负荷到当前时段，根据电力系统的运行参数可以建立电负荷柔性方程。同时，电氢能源系统可以与上级电网进行电功率交换，向上级电网购电或售电，根据电力系统的运行参数还可以建立购电售电约束方程。另外，电氢能源系统可以通过风电和光伏等可再生能源进行发电，根据电力系统的运行参数还可以建立可再生能源出力约束方程。
[0049]
具体实施时，电力系统的运行参数包括电力系统的网络参数、拓扑、电负荷数据和可再生能源出力数据。
[0050]
s120、根据电制氢系统的运行参数，建立氢负荷柔性方程和电制氢安全运行约束方程。
[0051]
本发明实施例中，电氢能源系统中的氢负荷同样具有一定的柔性，在供需平衡周期内用户在不同时段所需的氢需求有所不同，根据电制氢系统的运行参数可以建立氢负荷柔性方程。同时，电制氢设备在运行时受到电制氢系统的安全约束，根据电制氢系统的运行参数可以建立电制氢安全运行约束方程。
[0052]
具体实施时，电制氢系统的运行参数包括电制氢系统的网络参数、拓扑、氢负荷数据、设备类型和设备容量。
[0053]
s130、利用电力系统和电制氢系统的电功率平衡关系，建立电功率平衡约束方程。
[0054]
本发明实施例中，电氢能源系统中的电力系统和电制氢系统的电功率具有一定的平衡关系，可以建立相应的电功率平衡约束方程。
[0055]
s140、以调度周期内电氢能源系统的运行成本最低为目标函数，以电负荷柔性方
程、购电售电约束方程、可再生能源出力约束方程、氢负荷柔性方程、电制氢安全运行约束方程和电功率平衡约束方程为约束条件，建立电氢能源系统调度模型。
[0056]
本发明实施例中，在调度周期内，以电氢能源系统的运行成本最低为优化目标，并将上述步骤中建立各个约束方程作为约束条件，建立电氢能源系统调度模型。具体实施时，所建立的电氢能源系统调度模型可以是混合整数线性规划模型。
[0057]
s150、求解电氢能源系统调度模型，得到优化调度结果。
[0058]
本发明实施例中，在满足上述步骤中建立的约束条件前提下，求解该电氢能源系统调度模型的目标函数，最终可以得到电氢能源系统的优化调度结果。
[0059]
以上可知，本发明实施例提供的考虑柔性氢需求的电氢能源系统调度方法，通过建立电负荷柔性方程、购电售电约束方程、可再生能源出力约束方程、氢负荷柔性方程、电制氢安全运行约束方程和电功率平衡约束方程，并以调度周期内电氢能源系统的运行成本最低为目标函数，构建电氢能源系统调度模型，最后求解该模型可以得到电氢能源系统的优化调度结果，由于充分考虑了电负荷和氢负荷的柔性，可以为电氢综合能源系统运行提供额外的灵活性，以提高系统运行经济性，与现有技术相比，能够有效解决电、氢能流协同问题，同时考虑了电、氢负荷的柔性，进一步为系统运行提供额外的灵活性，通过所提出的调度方法获得经济可行的电氢综合能源系统调度计划，保证系统安全、高效、经济的运行。
[0060]
具体的，上述实施例中，目标函数的表达式为：
[0061][0062]
式中，t代表调度时段，t代表调度周期，代表电力系统在t时段向电网购电的单位成本，p
tgrid+
代表电力系统在t时段向电网购电的功率，代表电力系统在t时段向电网售电的单位成本，p
tgrid
‑
代表电力系统在t时段向电网售电的功率，c
shift
代表电力系统电负荷转移需要向用户支付的单位惩罚成本，p
te,in
代表电力系统在t时段从其它时段转入的电负荷。
[0063]
本发明实施例中，电氢能源系统运行调度的优化目标为调度周期内运行成本最低，优化变量为系统在每个时段向电网购电的功率p
tgrid+
、系统在每个时段向电网售电的功率p
tgrid
‑
和系统在每个时段从其它时段转入的电负荷p
te,in
。
[0064]
进一步地，上述实施例中，电负荷柔性方程的表达式为：
[0065][0066]
式中，p
te
代表电力系统在t时段的总电负荷，p
te,fix
代表电力系统在t时段的固定电负荷，p
te,out
代表电力系统在t时段转出到其他时段的电负荷，p
te,in,max
代表电力系统在t时段从其它时段可转入的电负荷的最大值，p
te,out,max
代表电力系统在t时段可转出到其他时
段的电负荷的最大值。
[0067]
本发明实施例中，电氢能源系统在一定时段内对电负荷进行转移时，从其它时段转入的电负荷p
te,in
和转出到其他时段的电负荷p
te,out
均受到一定的限制，即电氢能源系统中电负荷具有的柔性受到上述电负荷柔性方程的约束。
[0068]
更进一步地，上述实施例中，购电售电约束方程的表达式为：
[0069][0070]
式中，代表电力系统在t时段向电网购电时的变量，代表电力系统在t时段向电网售电时的变量，代表电力系统与电网所允许交换电功率的最大值。
[0071]
本发明实施例中，电氢能源系统在与上级电网进行电功率交换时，在一定时段内向上级电网购电的功率p
tgrid+
和售电的功率p
tgrid
‑
受到系统与电网所允许交换电功率的最大值的限制，即电氢能源系统中向上级电网的购电售电受到上述购电售电约束方程的约束。
[0072]
更进一步地，上述实施例中，可再生能源出力约束方程的表达式为：
[0073][0074]
式中，代表电力系统在t时段的光伏出力调度值，代表电力系统在t时段的光伏出力预测值，σ
t,fore
代表电力系统在t时段的光伏出力预测值的标准差，代表标准正态分布n(0,1)的逆累计分布函数，η代表光伏预测值的置信水平。
[0075]
本发明实施例中，电氢能源系统中的可再生能源主要为光伏，且光伏在一定时段内出力调度值受到光伏发电本身的限制，即电氢能源系统中可再生能源出力受到上述可再生能源出力约束方程的约束。
[0076]
更进一步地，上述实施例中，氢负荷柔性方程的表达式为：
[0077][0078]
式中，t
h
代表电制氢系统的氢负荷所需的供需平衡周期长度，k代表电制氢系统的氢负荷的供需平衡周期序列，p
th
代表电制氢系统在t时段用户的氢负荷；代表电制氢系统在第k个供需平衡周期内用户所需的氢需求。
[0079]
本发明实施例中，电氢能源系统在一定时段内用户的氢负荷p
th
在供需平衡周期内的波动受到一定的限制，即电氢能源系统中氢负荷具有的柔性受到上述氢负荷柔性方程的约束。
[0080]
更进一步地，上述实施例中，电制氢安全运行约束方程的表达式为：
[0081][0082]
式中，η
p2h
代表电制氢设备的转化效率，p
tp2h
代表在t时段电制氢设备消耗的电功率，λ代表电制氢设备的最小负载水平，代表在t时段电制氢设备的运行状态，c
p2h
代表电制氢设备的容量。
[0083]
本发明实施例中，电氢能源系统在一定时段内电制氢设备消耗的电功率p
tp2h
受到电制氢系统的安全约束，即电氢能源系统安全运行受到上述电制氢安全运行约束方程的约束。
[0084]
更进一步地，上述实施例中，电功率平衡约束方程的表达式为：
[0085]
p
tgrid+
‑
p
tgrid
‑
+p
tpv
‑
p
tp2h
＝p
te
。
[0086]
本发明实施例中，电氢能源系统在一定时段内向上级电网购电的功率p
tgrid+
、向上级电网售电的功率p
tgrid
‑
、光伏出力调度值和电制氢设备消耗的电功率p
tp2h
与电氢能源系统的总电负荷p
te
之间具有平衡关系，受到上述电功率平衡约束方程的约束。
[0087]
作为本发明一种优选的实施例，步骤s150包括：
[0088]
利用混合整数线性规划法求解电氢能源系统调度模型，得到优化调度结果。
[0089]
需要说明的是，优化问题的求解方法主要分为两大类：即启发式方法和数学优化方法。本发明实施例中，所建立的电氢能源系统调度模型为混合整数线性规划模型，属于混合整数线性规划问题，可以采用cplex、gurobi等求解器进行求解，得到电氢能源系统的优化调度结果。
[0090]
下面，结合具体的案例对本发明实施例提供的考虑柔性氢需求的电氢能源系统调度方法在经济性方面的优点进行验证。
[0091]
设置调度周期为1年，时间间隔为1h。为使仿真结果更具实用性，通过采集我国东部沿海某地区的小时级年电力负荷曲线作为电负荷数据，利用美国国家可再生能源实验室开发的pvwatts calculator计算年小时级单位光伏发电数据。电负荷和单位光伏出力如图2所示。
[0092]
此外，设电制氢设备的电氢转换效率为0.613(即每生产1kg氢需消耗54.3kmh电)，设备容量为2mw，最小负载水平为0.1。
[0093]
将柔性氢负荷的供氢时间尺度分别设置为1h、6h、12h、24h和168h，对应的氢需求分别为12.5kg/h、75.0kg/6h、150.0kg/12h、300.0kg/24h和2100.0kg/168h，研究从小时级到周级尺度下氢需求对电氢能源系统的影响。为了验证本发明所提出的方法在经济性方面的优势，光伏配置容量为20mw下不同柔性氢需求的系统经济性对比情况如图3所示。
[0094]
从图3中可以看出，不同光伏配置容量下，系统年运行成本均随着供氢时间尺度的增加而降低。具体地，随着时间尺度从1h到6h，系统的购电成本在不同光伏配置容量下均大幅降低，进而使系统总费用降低。随着时间尺度从6h到168h，系统的成本降低逐渐放缓，后期变化很小。以上结果说明，随着柔性氢需求供氢时间尺度的增加，系统购电费用、售电费用和补贴费用均逐渐降低，这有利于缓解对上级电网和需求侧响应的依赖作用，同时验证
了本发明实施例所提出的调度方法在运行经济性方面的优势。
[0095]
请参阅图4，图4展示了不同氢负荷柔性下电制氢设备的输入电功率结果。可以看出，在1h氢需求约束中，电制氢设备需保持恒功率运行来满足实时的用氢需求。当氢负荷平衡周期增加到6h时，电制氢设备运行灵活性迅速提升，在一天的调度周期内主要工作于3：00
‑
4：00,11：00
‑
12：00，14：00
‑
15：00和23：00
‑
24:00电价低谷和光伏出力高峰时段，因此系统经济性有所提升，但电制氢设备启停过于频繁(一天内启停三次)。而随着氢负荷平衡周期的进一步增加，电制氢设备进一步集中于在中午时段以额定功率进行工作，且在一个典型日内只启停一次。
[0096]
请参阅图5，图5展示了不同氢负荷柔性下的电负荷转移结果。可以看出，不同氢负荷平衡周期下系统主要将其他时段的电负荷需求转移到8：00
‑
17:00时段光伏出力高峰时段，且通过需求侧响应曲线可以发现，氢负荷平衡周期为1h时波动最大，随后随着时间尺度增加波动越来越小。以上结果可以发现，随着氢负荷平衡周期逐渐增大(终端存储能力逐渐增强)，系统在设备运行、与上级电网交易和需求侧响应等方面效果显著。
[0097]
本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。