风电‑氢储能与煤化工多能耦合系统气体分配控制方法与流程

本发明涉及一种风电-氢储能与煤化工多能耦合系统的气体分配控制方法。

背景技术：

在当今中国国民日益增强的健康观念和现实环境“霾影重重”的严峻背景下，国家发改委和能源局在“十三五”规划中将煤炭的清洁高效利用与储能技术、可再生能源研究等一同定为能源发展八大重点工程，以此致力于实现“北煤南运、西煤东调”的总体战略性布局。借此，提出利用氢清洁储能，重组并发挥我国新疆等北部地区兼备风、煤资源的优势，消除风电大规模并网消纳瓶颈，提升煤基化工产业的环境友好度，在切实实现煤炭优化目标的同时，利用煤、电辅助氢储能系统打破传统氢气存储与输送限局。在此过程中，如何依据风、氢、煤各子系统实际运行工况和需求，对多能耦合系统内部能量进行协调控制，是发挥各子系统优势，切实提升大规模风电并网消纳能力，优化储能工艺，改善煤化工工艺清洁度的重要环节。

现有的国内外研究主要致力于利用氢储能技术平滑电源功率波动，实现电力削峰填谷，提高系统备用容量。或避开风电上网的瓶颈问题，将煤化工产业作为非并网风电的一条消纳路径，间接提高风电利用率，减少弃风。在风、氢、煤集成系统控制研究方面，文献[袁铁江,段青熙,秦艳辉,等.风电–氢储能与煤化工多能耦合系统能量广域协调控制架构[j].高电压技术,2016,42(9):2748-2755.]基于分层递阶控制原理，对风电-氢储能与煤化工多能耦合系统广域协调控制方案进行了探讨。专利2015102890549提出了一种分散式接入电场的风-氢储能耦合系统控制方法，基于风电场实际出力与本地负荷差值，协同上级电网调度计划，以氢储能等效soc状态为判断依据，给出了可弃风情况下，风-氢储能耦合系统的风电功率分配方法。专利2015105419016整合了氢、氧储气罐状态给出氢储能系统等效soc状态，并以此为前提，详细讨论了风电-氢储能与煤化工多能耦合系统内电能需求情况，并以提高分散式风电消纳能力为目标，煤化工系统不间断稳定运行为约束，对系统内部电能分配进行控制研究。鉴于风电–氢储能与煤化工多能耦合系统内部存在化学能链路，因此，为保证多能耦合系统内部各子系统协调运行，还需基于上述电能分配控制方案，结合其内部电能与化学能的转换关系，对多能耦合系统内气体分配进行控制研究。

技术实现要素：

为克服现有技术的缺点，完善上述多能耦合系统控制方法，本发明提出一种风电-氢储能与煤化工多能耦合系统的气体分配控制方法。

本发明基于由风力发电系统、氢储能系统、煤化工系统和电网构成的风电-氢储能与煤化工多能耦合系统，所述的氢储能系统由制氢、储氢的等效充电环节和氢氧燃料电池等效放电环节两部分组成。本发明基于该多能耦合系统的实际运行情况，以气体产、耗过程动态平衡为原则，以灵活调整优先级顺序为主导思想，以风电消纳最大化为目标，以煤化工系统额定不间断运行、氢储能系统储能状态处于正常区间为约束，协调控制多能耦合系统内部气体的流向和流量，实现多能耦合系统内部气体合理分配。

本发明的具体步骤为：

1、以制得满足煤化工系统以额定功率运行时所需氢气的功率pmhg和煤化工系统以额定功率运行时消耗氢气的流率nmhg为基准，依据氢储能系统储能状态esocs划分氢储能系统控制区间，制定各区间内控制的优先级顺序，并据此给出风电出力用于电解水的功率公式

2、基于构建的的气体分配基本控制公式，分析不同的优先级情况下，多能耦合系统内部可能出现的气体流动情况；

3、在步骤1和步骤2的基础上，通过研究氢氧燃料电池等效充放电控制情况，构建完整的多能耦合系统气体分配控制公式；

4、考虑多能耦合系统可能出现的极限情况，确定多能耦合系统控制时间间隔取值。

以上步骤具体如下：

1、所述步骤1中，

氢储能系统储能状态esocs＝(esoch×vhcap+esoco×vocap)/(vhcap+vocap)，

其中：esoch为氢储能系统中储氢罐的状态：

esoch＝phre/phcap

式中：phre为储氢罐剩余压强，phcap为储氢罐额定压强。

esoco为氢储能系统中储氧罐的状态：

esoco＝pore/pocap

式中：pore为储氧罐剩余压强，pocap为储氧罐额定压强，vhcap为储氢罐容积，vocap为储氧罐容积。

得到氢储能系统控制区间划分及各区间内控制的优先级顺序如表1所示：

表1氢储能系统控制区间划分及各区间内控制依据

表中，esocs_max为氢储能系统储能状态的上限；esocs_min为氢储能系统储能状态的下限；pjh为集群功率调度计划；pwind为与集群功率调度计划pjh时间尺度匹配的风电场实际出力；pload为与集群功率调度计划pjh时间尺度匹配的本地负荷；i为区间定义标号，取值为1、2、3；为氢储能系统储能状态esocs处于区间i时风电出力满足对应情况下并网的必须需求后剩余的功率。

则风电出力用于电解水的功率为：

其中，

式中：pmhg为以制得满足煤化工系统以额定功率运行时所需氢气的功率，为氢储能系统储能状态esocs处于区间i时，风电出力满足对应情况下并网的必须需求后剩余的功率。

2、所述步骤2中，由气体产、耗过程动态平衡可知，在任意时间尺度下，多能耦合系统运行过程中电解水产生气体的物质的量和消耗储气罐内原有气体的物质的量之和等于煤化工系统和氢氧燃料电池运行所需气体的物质的量之和，据此得到气体分配基本控制公式为：

式中：为δt时间内电解水产生氢气的物质的量；为δt时间内电解水产生氧气的物质的量；nδh为δt时间内应消耗储气罐内原有氢气的物质的量，nδh大于0表示储存的氢气减少、nδh小于0表示储存的氢气增加；nδo为δt时间内应消耗储气罐内原有氧气的物质的量，nδo大于0表示储存的氧气减少、nδo小于0表示储存的氧气增加；nmhg_h为δt时间内煤化工系统消耗氢气的物质的量；nmhg_o为δt时间内煤化工系统消耗氧气的物质的量；nrh为δt时间内需通入氢氧燃料电池的氢气的物质的量；nro为δt时间内需通入氢氧燃料电池的氧气的物质的量；δt为系统控制时间间隔。

其中，δt时间内电解水产生氢气的物质的量氧气的物质的量为：

式中：为电解水制氢气摩尔流率，δt为系统控制时间间隔，ηe为电解池电解效率，为风电出力用于电解水的功率，ue为电解槽电压，f为法拉第常数。

进而，视煤化工系统进气流速、流量为定值，且整个过程为理想状态，则有δt时间内煤化工系统消耗氢气的物质的量nmhg_h、氧气的物质的量nmhg_o为：

式中：nmhg煤化工系统以额定功率运行时消耗氢气的流率，mw为风电场容量，单位：mw，vm为气体摩尔体积22.4l/mol。

据此，电解水制氢气量、氧气量和煤化工系统消耗氢气量、氧气量均能可直接求得，但因储气罐内原有气体的消耗情况和氢氧燃料电池放电的用气情况依据实际运行工况的变化而变化，所以，为实现多能耦合系统气体分配控制意愿，应重点分析氢储能系统储能状态esocs与氢氧燃料电池运行状态之间的关系，即为对氢储能系统等效充放电情况分析如下：

不考虑氢氧燃料电池出力的上下阈值，依次就3种不同优先级情况下氢储能系统运行状态进行分析，以下分析以氢气路径为例，氧气路径同理可得：

(1)esocs>esocs_max，时，氢氧燃料电池启动：

①当时，有氢氧燃料电池放电，补足本地负荷pload与集群功率调度计划pjh之和，即氢氧燃料电池放电功率进而可知nrh，则

②当时，有氢氧燃料电池放电，调节esocs降至esocs_max，即nδh和nmhg_h已知，则并可推导出pfc。

(2)esocs_min≤esocs≤esocs_max，时：

为保证进气量为0，出气量最大的恶劣情况下，在控制时间内，esocs不低于下限，则设氢氧燃料电池启停临界值esocsa为：

式中：nrh_max为氢氧燃料电池以最大功率输出时，通入氢气的流率，r为理想气体常数，th为储氢罐内温度，to为储氧罐内温度，γh-o为通入氢氧燃料电池入口的氢气和氧气摩尔流率比，phcap为储氢罐额定压强，pocap为储氧罐额定压强。

①当时，氢氧燃料电池放电保证本地负荷需求，进而可得nrh，则

②当或时，氢氧燃料电池停止，pfc＝0，nrh＝0，则

(3)esocs＜esocs_min，时，氢氧燃料电池停止，pfc＝0，nrh＝0，则

上述过程中：esocs为氢储能系统储能状态；esocs_max为氢储能系统储能状态的上限；esocs_min为氢储能系统储能状态的下限；vhcap为储氢罐容积；vocap为储氧罐容积；pjh为集群功率调度计划；pwind为与集群功率调度计划pjh时间尺度匹配的风电场实际出力；pload为与集群功率调度计划pjh时间尺度匹配的本地负荷；pmhg为以制得满足煤化工系统以额定功率运行时所需氢气的功率；δt为控制时间间隔；i为区间定义标号，取值为1、2、3；为氢储能系统储能状态esocs处于区间i时，风电出力满足对应情况下并网的必须需求后剩余的功率；为风电出力用于电解水的功率；pfc为氢氧燃料电池放电功率；为δt时间内电解水产生氢气的物质的量；为δt时间内电解水产生氧气的物质的量；nδh为δt时间内应消耗储气罐内原有氢气的物质的量，nδh大于0表示储存的氢气减少、nδh小于0表示储存的氢气增加；nδo为δt时间内应消耗储气罐内原有氧气的物质的量，nδo大于0表示储存的氧气减少、nδo小于0表示储存的氧气增加；nmhg_h为δt时间内煤化工系统消耗氢气的物质的量；nmhg_o为δt时间内煤化工系统消耗氧气的物质的量；nrh为δt时间内需通入氢氧燃料电池的氢气的物质的量；nro为δt时间内需通入氢氧燃料电池的氧气的物质的量。

由于氢氧燃料电池作为一种能量转换装置，通过电化学反应将化学能直接转换为电能，并不能存储能量，这使其仅具备启动和停止两种模式，分别记为flag1＝1和flag1＝0。继而，通过对氢氧燃料电池充放电状态的建模，构建完整气体分配控制模型。

3、所述步骤3依据步骤1和步骤2的综合分析，得到三种氢氧燃料电池等效充放电控制情况如下：

(1)当时，消耗储氢罐和储氧罐中气体调节esocs，等效为降低esocs值。此时flag1＝1。逻辑顺序为：已知δt时间内应消耗储气罐内原有氢气的物质的量nδh、氧气的物质的量nδo，可推导出δt时间内需通入氢氧燃料电池的氢气的物质的量nrh、氧气的物质的量nro，进而得到氢氧燃料电池放电功率pfc。

此时，氢储能系统储能状态esocs超出上限，这表明氢储能系统中储氢罐的状态超出上限，或氢储能系统中储氧罐的状态超出上限，或两者均超出上限。该时刻表征氢储能系统中储氢罐和储氧罐状态超出上限的部分的物质的量公式为：

式中：phcap为储氢罐额定压强，pocap为储氧罐额定压强，vhcap为储氢罐容积，vocap为储氧罐容积，esoch为氢储能系统中储氢罐的状态，esoco为氢储能系统中储氧罐的状态，esoch_max为氢储能系统中储氢罐状态的上限，esoco_max氢储能系统中储氧罐状态的上限，r为理想气体常数，th为储氢罐内温度，to为储氧罐内温度。

则此时氢储能系统储能环节可能出现的情况及计算该时刻应消耗气罐中氢气、氧气的物质的量需使用的参数如表2所示：

表2氢储能系统超上限情况及对应情况选取的计算参数

表中：n'δh为氢储能系统中储氢罐状态超出上限的部分的物质的量，n'δo为氢储能系统中储氧罐状态超出上限的部分的物质的量，nδh为δt时间内应消耗储气罐内原有氢气的物质的量，nδo为δt时间内应消耗储气罐内原有氧气的物质的量。

因此，当分别选取储氢超上限部分和储氧超上限部分作为对象进行计算时，对应δt时间内需通入氢氧燃料电池的氢气的物质的量nrh、氧气的物质的量nro分别为：

a、选取储氢超上限部分进行计算时：

式中：为δt时间内电解水产生氢气的物质的量；为δt时间内电解水产生氧气的物质的量；nδh为δt时间内应消耗储气罐内原有氢气的物质的量，nδh大于0表示储存的氢气减少，nδh小于0表示储存的氢气增加；nδo为δt时间内应消耗储气罐内原有氧气的物质的量，nδo大于0表示储存的氧气减少、nδo小于0表示储存的氧气增加；nmhg_h为δt时间内煤化工系统消耗氢气的物质的量；nmhg_o为δt时间内煤化工系统消耗氧气的物质的量；γh-o为通入氢氧燃料电池入口的氢气和氧气摩尔流率比。

b、选取储氧超上限部分进行计算时：

由于本发明所述的氢氧燃料电池由多个单元氢氧燃料电池组成，则根据物料守恒，得到单元氢氧燃料电池内氢气分压pfc_h、氧气分压pfc_o，单位atm，分别为：

式中：n为氢氧燃料电池单元个数，va为单元氢氧燃料电池阳极体积，vca为单元氢氧燃料电池阴极体积，t为氢氧燃料电池的工作温度，ifc为氢氧燃料电池负载电流，为氢气阀摩尔常数，为氧气阀摩尔常数。

进而，由氢氧燃料电池模型可得氢氧燃料电池放电功率pfc为：

式中：ηf为氢氧燃料电池实际发电效率，ufc为氢氧燃料电池输出电压，pfc_min为氢氧燃料电池出力范围的上限，pfc_max为氢氧燃料电池出力范围的下限。

其中，氢氧燃料电池输出电压ufc：

ufc＝n·ucell＝n·(enernst-uohmic-uc)

式中：ucell为单节氢氧燃料电池输出电压，enernst为热力学电动势，uohmic为欧姆极化过电压，uc为等效电压。

热力学电动势：

enernst＝1.229-8.5×10^-4(t-298.15)+4.3085×10^-5t(lnpfc_h+0.5lnpfc_o)

欧姆极化过电压：

式中：zm为等效膜阻抗，zc为阻碍质子过膜阻抗，l为质子交换膜厚度，a为质子交换膜的活化面积，为质子交换膜含水量。

电池等效电压：

式中：b为方程系数，j为实际电流密度，jmax为能达到的最大电流密度，c为电池等效电容，ε1、ε2、ε3、ε4为氢氧燃料电池系统经验参数。

(2)当或时，氢储能系统储能状态esocs允许，但风电并网功率不能满足预计需求，氢氧燃料电池放电辅助风电并网，flag1＝1。逻辑顺序为：已知pfc，推导出nrh、nro，进而得到nδh、nδo。此时，氢氧燃料电池放电功率pfc为：

进而，由氢氧燃料电池放电功率pfc＝ηf(ufc·ifc)可反推出δt时间内需通入氢氧燃料电池的氢气的物质的量nrh和氧气的物质的量nro。

(3)当或时，氢氧燃料电池停止，flag1＝0，氢氧燃料电池放电功率pfc＝0，nrh＝0，

至此，多能耦合系统气体分配控制公式构建完毕。进而，得到储氢罐剩余压强，即下一时刻控制开始时储氢罐实际压强phre，单位，pa：

phre＝ph_0+δph_1+δph_2-δph_3-δph_4

其中，

式中：ph_0为储氢罐初始压强，δph_1为电解水制氢气增压，δph_2为回收氢氧燃料电池未完全反应的氢氧气致储气罐增压，δph_3为用于煤化工反应的氢气输出减压，δph_4为通入氢氧燃料电池的氢气输出减压，k'为气体回收率。

则储氧罐剩余压强pore，单位，pa：

式中：po_0为储氧罐初始压强。

4、所述步骤4中，为使气体分配控制策略合理、有效，保证氢储能系统储能状态在控制时间间隔δt内不会出现骤降或激增，所述的骤降为δt内氢储能系统储能状态esocs从超出氢储能系统储能状态的上限esocs_max，变为低于氢储能系统储能状态的下限esocs_min，所述的激增为δt内氢储能系统储能状态esocs从低于氢储能系统储能状态的下限esocs_min，变为超出氢储能系统储能状态的上限esocs_max，若在这两种氢储能系统运行的极端情况下，控制时间间隔δt仍能满足系统运行需求，确保气体控制策略正常实现控制效果，则能确保多能耦合系统在运行过程中，不会出现因控制时间不合理造成的系统运行失稳情况。因此，基于以上两种情况，确定控制时间间隔δt的取值。

a、以储气罐为研究对象，当以最大流率充气，以最小流率放气时，符合的情况。此时，esocs＜esocs_min，风电出力除了取4倍的制得满足煤化工系统以额定功率运行时所需氢气的功率pmhg值用于电解水之外，全部并网。为保证在δta时间内氢储能系统中储氢气罐的状态esoch和氧罐的状态esoco均不超出上限，δta应取δtah、δtao中的最小值：

式中：δta为当以最大流率给储气罐充气，以最小流率给储气罐放气时，应保证的多能耦合系统气体分配控制时间间隔；δtah为当以最大流率给储氢罐充气，以最小流率给储氢罐放气时，应保证的控制时间间隔；δtao为当以最大流率给储氧罐充气，以最小流率给储氧罐放气时，应保证的控制时间间隔；esoch_max为氢储能系统中储氢罐状态的上限；esoch_min为氢储能系统中储氢罐状态的下限；esoco_max为氢储能系统中储氧罐状态的上限；esoco_min为氢储能系统中储氧罐状态的下限；phcap为储氢罐额定压强；pocap为储氧罐额定压强；vhcap为储氢罐容积；vocap为储氧罐容积；nmhg为煤化工系统以额定功率运行时消耗氢气的流率；r为理想气体常数；th为储氢罐内温度；to为储氧罐内温度。

b、以储气罐为研究对象，当以最小流率充气，以最大流率放气时，符合步骤3中的情况(2)。鉴于氢氧燃料电池启停临界值esocsa是依据δt推出的，因此在此仅需讨论情况，此时esocs>esocs_max，为保证在δtb时间内氢储能系统中储氢罐的状态esoch和储氧罐的状态esoco均不低于其下限，δtb应取δtbh、δtbo中最小值：

式中：δtb为当以最小流率给储气罐充气，以最大流率给储气罐放气时，应保证的多能耦合系统气体分配控制时间间隔；δtah为当以最小流率给储氢罐充气，以最大流率给储氢罐放气时，应保证的控制时间间隔；δtao为当以最小流率给储氧罐充气，以最大流率给储氧罐放气时，应保证的控制时间间隔；nrh_max为氢氧燃料电池以最大功率输出时，通入氢气的流率。

则控制时间间隔δt取值应满足：

δt＜min(δta，δtb)。

附图说明

图1风电-氢储能与煤化工多能耦合系统气体分配控制方法流程图；

图2氢储能系统运行极端情况示意图。

具体实施方式

下面结合附图及具体实施方式对本发明做进一步说明。

本发明风电-氢储能与煤化工多能耦合系统气体分配控制方法流程图如图1所示，图1基于具体步骤1、步骤2和步骤3，反映了多能耦合系统气体分配控制方法。为控制多能耦合系统内部电能与非电能量依需流动、转换，协调各子系统运行状态，确保在煤化工系统额定不间断运行及氢储能系统储能状态处于正常区间的前提下，实现风电消纳最大化的目标，本发明依据氢储能系统储能状态，以灵活调整优先级顺序为主导思想，梳理集群功率调度计划pjh、与集群功率调度计划pjh时间尺度匹配的风电场实际出力pwind、与集群功率调度计划pjh时间尺度匹配的本地负荷pload三者之间的关系，基于风电出力用于电解水的功率公式，以气体产、耗过程动态平衡的原则，给出风电-氢储能与煤化工多能耦合系统气体分配控制方法。

对应步骤2中，气体分配基本控制公式由任意时间尺度下，氢储能系统电解水产生气体的物质的量和消耗储气罐内原有气体的物质的量之和等于煤化工系统和氢氧燃料电池运行所需气体的物质的量之和，可得：

其中：δt时间内电解水产生氢气的物质的量氧气的物质的量为：

式中：为电解水制氢气摩尔流率，δt为系统控制时间间隔，ηe为电解池电解效率，为风电出力用于电解水的功率，ue为电解槽电压，f为法拉第常数

对应步骤1中，风电出力用于电解水的功率为：

其中，

视煤化工系统进气流速、流量为定值，且整个过程为理想状态，则有δt时间内煤化工系统消耗氢气的物质的量nmhg_h和消耗氧气的物质的量nmhg_o为：

式中：nmhg为煤化工系统以额定功率运行时消耗氢气的流率，mw为风电场容量，单位：mw，vm为气体摩尔体积22.4l/mol。

且所述风电-氢储能与煤化工多能耦合系统气体分配情况如下：

不考虑氢氧燃料电池出力上下阈值，依次分析3种不同优先级情况下氢储能系统运行状态：

(1)esocs>esocs_max，时，氢氧燃料电池启动：

①当时，有氢氧燃料电池放电补足本地负荷与集群功率调度计划之和，即氢氧燃料电池放电功率进而可知nrh、nro，则

②当时，有氢氧燃料电池放电调节esocs降至esocs_max，此时nδh、nδo已知，则nrh＝nδh+nh2o_h-nmhg_h、可推出pfc。

(2)esocs_min≤esocs≤esocs_max，时：

为保证进气量为0，出气量最大的恶劣情况下，在控制时间内，esocs不低于下限，则设氢氧燃料电池启停临界值esocsa为：

①当时，氢氧燃料电池放电保证本地负荷需求，进而可得nrh、nro，则

②当或时，氢氧燃料电池停止，pfc＝0，nrh＝0、nro＝0，则

(3)esocs＜esocs_min，时，氢氧燃料电池停止，pfc＝0，nrh＝0、nro＝0，则

上述过程中，氢储能系统储能状态esocs为：

其中：phre为储氢罐剩余压强，phcap为储氢罐额定压强；esoch为氢储能系统中储氢罐的状态，esoco为氢储能系统中储氧罐的状态；pore为储氧罐剩余压强，pocap为储氧罐额定压强；vhcap为储氢罐容积，vocap为储氧罐容积；esocs_max为氢储能系统储能状态的上限、esocs_min为氢储能系统储能状态的下限；pmhg为以制得满足煤化工系统以额定功率运行时所需氢气的功率；δt为控制时间间隔；i为区间定义标号，取值为1、2、3；为氢储能系统储能状态esocs处于区间i时风电出力满足对应情况下并网的必须需求后剩余的功率；为风电出力用于电解水的功率；pfc为氢氧燃料电池放电功率；esocsa为氢氧燃料电池启停临界值，γh-o为通入氢氧燃料电池入口的氢气和氧气摩尔流率比，nrh_max为氢氧燃料电池以最大功率输出时通入氢气的流率，r为理想气体常数，th为储氢罐内温度，to为储氧罐内温度；为δt时间内电解水产生氢气的物质的量，为δt时间内电解水产生氧气的物质的量，nδh为δt间内应消耗储气罐内原有氢气的物质的量，nδh大于0表示储存的氢气减少、nδh小于0表示储存的氢气增加；nδo为δt时间内应消耗储气罐内原有氧气的物质的量，nδo大于0表示储存的氧气减少，nδo小于0表示储存的氧气增加；nmhg_h为δt时间内煤化工系统消耗氢气的物质的量，nmhg_o为δt时间内煤化工系统消耗氧气的物质的量；nrh为δt时间内需通入氢氧燃料电池的氢气的物质的量，nro为δt时间内需通入氢氧燃料电池的氧气的物质的量。

则由气体分配基本控制公式可知，在和nmhg_h、nmhg_o均能计算确定的情况下，只要知道nδh、nδo和nrh、nro任意一方，便能得到另一方的值，进而可以求得一个周期结束时氢储能系统的状态。据此，对应步骤3，将氢氧燃料电池充放电归纳为以下三种情况，依次构建对应情况下氢氧燃料电池充放电模型，给出nδh、nδo和nrh、nro求解方法：

(1)当时，消耗储氢罐和储氧罐内气体调节esocs，等效为降低esocs值。此时flag1＝1。逻辑顺序为：已知nδh、nδo，推导出nrh、nro，进而得到pfc。

此时，氢储能系统储能状态esocs超出上限，这表明氢储能系统中储氢罐的状态esoch超出上限或储氧罐的状态esoco超出上限或两者均超出上限。该时刻表征氢储能系统中氢和储氧罐状态超出上限的部分的物质的量公式为：

则计算该时刻应消耗气罐内原有氢气、氧气的物质的量nδh、nδo的参考值选取方法存在以下两种情况：

a、当(n'δh>0)&(n'δo≤0)或(n'δh≥2.3n'δo＞0)时，以n'δh为计算参考值，有nδh＝n'δh，则nro＝nrh/γh-o；

b、当(n'δo>0)&(n'δh≤0)或(2.3n'δo＞n'δh＞0)时，以n'δo为计算参考值，有nδo＝n'δo，则nrh＝γh-onro。

上述情况中：n'δh为氢储能系统中储氢罐状态超出上限的部分的物质的量，n'δo为氢储能系统中储氧罐状态超出上限的部分的物质的量，为δt时间内电解水产生氢气的物质的量，为δt时间内电解水产生氧气的物质的量，nδh为δt时间内应消耗储气罐内原有氢气的物质的量，nδo为δt时间内应消耗储气罐内原有氧气的物质的量，nmhg_h为δt时间内煤化工系统消耗氢气的物质的量，nmhg_o为δt时间内煤化工系统消耗氧气的物质的量；nrh为δt时间内需通入氢氧燃料电池的氢气的物质的量，nro为δt时间内需通入氢氧燃料电池的氧气的物质的量，γh-o为通入氢氧燃料电池入口的氢气和氧气摩尔流率比。

由于本发明所述的氢氧燃料电池由多个单元氢氧燃料电池组成，根据物料守恒，得到单元氢氧燃料电池内氢气分压pfc_h、氧气分压pfc_o，单位atm，分别为：

式中：n为氢氧燃料电池单元个数，va为单元氢氧燃料电池阳极体积，vca为单元氢氧燃料电池阴极体积，t为氢氧燃料电池的工作温度，ifc为氢氧燃料电池负载电流，为氢气阀摩尔常数，为氧气阀摩尔常数。得到该情况下氢氧燃料电池放电功率pfc为：

式中：ηf为氢氧燃料电池实际发电效率，ufc为氢氧燃料电池输出电压，pfc_min为氢氧燃料电池出力范围的上限，pfc_max为氢氧燃料电池出力范围的下限。

其中，氢氧燃料电池输出电压ufc：

ufc＝n·ucell＝n·(enernst-uohmic-uc)

式中：ucell为单节氢氧燃料电池输出电压，enernst为热力学电动势，uohmic为欧姆极化过电压，uc为等效电压。

热力学电动势：

enernst＝1.229-8.5×10^-4(t-298.15)+4.3085×10^-5t(lnpfc_h+0.5lnpfc_o)

欧姆极化过电压：

式中：zm为等效膜阻抗，zc为阻碍质子过膜阻抗，l为质子交换膜厚度，a为质子交换膜的活化面积，为质子交换膜含水量。

电池等效电压：

式中：b为方程系数，j为实际电流密度，jmax为能达到的最大电流密度，c为电池等效电容，ε1、ε2、ε3、ε4为氢氧燃料电池系统经验参数。

(2)当或时，氢储能系统储能状态esocs允许，但风电并网功率不能满足预计需求，氢氧燃料电池放电辅助风电并网，flag1＝1。逻辑顺序为：已知pfc，推导出nrh、nro，进而得到nδh、nδo。

进而，由前述氢氧燃料电池放电功率pfc＝ηf(ufc·ifc)可反推出δt时间内需通入氢氧燃料电池的氢气的物质的量nrh、氧气的物质的量nro，则此时

(3)当或时，氢氧燃料电池停止flag1＝0，pfc＝0，nrh＝0、nro＝0，此时

综上，储氢罐剩余压强，即下一时刻控制开始时储氢罐实际压强phre，单位，pa：

phre＝ph_0+δph_1+δph_2-δph_3-δph_4(12)

其中，

式中：ph_0为储氢罐初始压强，δph_1为电解水制氢气增压，δph_2为回收氢氧燃料电池未完全反应的氢氧气致储气罐增压，δph_3为用于煤化工反应的氢气输出减压，δph_4为通入氢氧燃料电池的氢气输出减压，k'为气体回收率。

则储氧罐剩余压强pore，单位，pa：

式中：po_0为储氧罐初始压强。

图2所示为氢储能系统运行极端情况，是对具体步骤4的反应。为确保氢储能系统储能状态在控制时间间隔δt内不会出现骤降或激增，所述的骤降为δt内氢储能系统储能状态esocs从超出氢储能系统储能状态的上限esocs_max，变为低于氢储能系统储能状态的下限esocs_min，所述的激增为δt内氢储能系统储能状态esocs从低于氢储能系统储能状态的下限esocs_min，变为超出氢储能系统储能状态的上限esocs_max，保证控制策略的合理性，δt取值应综合考虑如图2所示的当储气罐以最大流率充气、最小流率放气和最小流率充气、最大流率放气两种极端运行情况，当这两种情况下控制时间间隔δt的取值仍能保证气体控制策略正常发挥控制效果，则能确保在整个控制时间内系统不会出现因控制策略制定不合理造成系统运行失稳的情况。

因此，基于以上两种情况，讨论控制时间间隔δt的取值。

a、如图2中(a)、(c)两种情况，以储气罐为研究对象，当以最大流率充气，以最小流率放气时，符合的情况。此时，氢储能系统中储氢罐的状态esoch或储氧罐的状态esoco低于下限，或两者均低于下限，风电出力除了取4pmhg值用于电解水之外，全部并网。为保证在δta时间内氢储能系统中储氢罐的状态esoch和氧罐的状态esoco均不超出其上限，δta应取δtah、δtao中最小值：

式中：δta为当以最大流率给储气罐充气，以最小流率给储气罐放气时，应保证的多能耦合系统气体分配控制时间间隔；δtah为当以最大流率给储氢罐充气，以最小流率给储氢罐放气时，应保证的控制时间间隔；δtao为当以最大流率给储氧罐充气，以最小流率给储氧罐放气时，应保证的控制时间间隔；esoch_max为氢储能系统中储氢罐状态的上限；esoch_min为氢储能系统中储氢罐状态的下限；esoco_max为氢储能系统中储氧罐状态的上限；esoco_min为氢储能系统中储氧罐状态的下限；phcap为储氢罐额定压强；pocap为储氧罐额定压强；vhcap为储氢罐容积；vocap为储氧罐容积；nmhg为煤化工系统以额定功率运行时消耗氢气的流率；r为理想气体常数；th为储氢罐内温度；to为储氧罐内温度。

b、如图2中(b)、(d)两种情况所示，以储气罐为研究对象，当以最小流率充气，以最大流率放气时，符合氢储能系统储能状态esocs允许，但风电并网功率不能满足预计需求，氢氧燃料电池放电辅助风电并网的情况，此时或鉴于氢氧燃料电池启停临界值esocsa是基于δt推出得到，因此在此仅需讨论情况。此时，氢储能系统中储氢罐的状态esoch或储氧罐的状态esoco超出上限，或两者均超出上限，为保证在δtb时间内氢储能系统中储氢罐的状态esoch和储氧罐的状态esoco均不低于其下限，则δtb应取δtbh、δtbo中最小值：

式中：δtb为当以最小流率给储气罐充气，以最大流率给储气罐放气时，应保证的多能耦合系统气体分配控制时间间隔；δtah为当以最小流率给储氢罐充气，以最大流率给储氢罐放气时，应保证的控制时间间隔；δtao为当以最小流率给储氧罐充气，以最大流率给储氧罐放气时，应保证的控制时间间隔；nrh_max为氢氧燃料电池以最大功率输出时，通入氢气的流率。

综上，控制时间间隔δt取值应满足：

δt＜min(δta，δtb)。