CN104113084B - 风电-制氢并网发电***的控制方法 - Google Patents

Abstract

本发明是一种风电-制氢并网发电***的控制方法，包括建立基于双馈感应发电机的风电机组模型、建立电解槽模型和对并网发电***进行控制等内容，能够充分反应对风电-制氢并网发电***建模与控制的有效性；提高了电网对新能源的接纳能力，适应性强，具有较高的实际应用价值。

Description

风电-制氢并网发电***的控制方法

技术领域

本发明是一种风电-制氢并网发电***的控制方法，应用于风电-制氢并网发电***建模仿真、并网控制策略研究、及其与电网相互作用机理研究与应用。

背景技术

由于风电功率输出具有随机性、间歇性和波动性，对电网安全稳定运行造成不利的影响。为有效解决这一问题，提高风电输出功率品质，许多科研工作者积极研究和探索风电-制氢并网发电***的理论意义和工程应用价值。通过采用科学创新的技术手段，实现了风力发电***与电解槽协调优化控制策略达到输出功率平稳，进而解决了风力发电***并网运行所产生的功率波动等技术难题，为风电大规模开发利用提供了有力的理论支撑和技术指导。风电-制氢并网发电***对提高了风能利用率，推进我国风能、氢资源的大规模开发利用，实现能源可持续发展有着非常重要的实际意义。

发明内容

本发明的目的是，提供一种电网对新能源吸纳能力强、适应性好、具有较高实用价值的风电-制氢并网发电***的控制方法。

本发明的目的是由以下技术方案来实现的：一种风电-制氢并网发电***的控制方法，其特征是，它包括以下步骤：

1)建立基于双馈感应发电机的风电机组模型

风电机组空气动力学数学模型为：

$P_M={\mathrm{\ρ}}_{\mathrm{air}}{C}_p(\mathrm{\λ},\mathrm{\β})\mathrm{\π}{R}^2{V}_w^3/2---\left(1\right)$

其中：P_M为风轮机捕获的风能转化成风电机组的机械功率，ρ_air为风电接入点的空气密度，C_p为叶片的风能转换效率系数，是风轮机叶尖速比和桨距控制角的函数，λ为风轮机的叶尖速比，β为风电机组的桨距控制角，π为圆周率，一般取3.1415926，R为风轮机叶轮半径，V_w为风电接入点的时时风速；

风轮机与发电机的两质块轴系数学模型方程为：

$\left\{\begin{array}{c}2H_{T}d{\mathrm{\ω}}_T/\mathrm{dt}=T_M-K_S{\mathrm{\θ}}_S-D_T{\mathrm{\ω}}_T\\ 2H_{G}d{\mathrm{\ω}}_G/\mathrm{dt}=K_S{\mathrm{\θ}}_S-T_E-D_G{\mathrm{\ω}}_G\\ d{\mathrm{\θ}}_S/\mathrm{dt}={\mathrm{\ω}}_0({\mathrm{\ω}}_T-{\mathrm{\ω}}_G)\end{array}\right.---\left(2\right)$

其中：H_T为风轮机的惯性常数，ω_T为风轮机的电角速度，T_M为风轮机的机械转矩，K_S为风轮机和发电机轴的刚度系数，θ_S为两质块之间相对角位移，D_T为风轮机转子阻尼系数，H_G为发电机的惯性常数，ω_G为感应发电机的电角速度，T_E为发电机的电磁转矩，D_G为发电机转子阻尼系数，ω₀为电网的同步角速度，dω_T/dt为风轮机角速度对时间的导数，dω_G/dt为感应发电机电角速度对时间的导数，dθ_S/dt为两质块之间相对角位移对时间的导数；

同步旋转坐标系下双馈感应发电机的电压方程为：

$\left\{\begin{array}{c}{u}_{\mathrm{sd}}=d{\mathrm{\ψ}}_{\mathrm{sd}}/\mathrm{dt}-{\mathrm{\ω}}_s{\mathrm{\ψ}}_{\mathrm{sq}}+R_s{i}_{\mathrm{sd}}\\ u_{\mathrm{sq}}=d{\mathrm{\ψ}}_{\mathrm{sq}}/\mathrm{dt}-{\mathrm{\ω}}_s{\mathrm{\ψ}}_{\mathrm{sd}}+R_s{i}_{\mathrm{sq}}\\ u_{\mathrm{rd}}=d{\mathrm{\ψ}}_{\mathrm{rd}}/\mathrm{dt}-s{\mathrm{\ω}}_s{\mathrm{\ψ}}_{\mathrm{rq}}+R_r{i}_{\mathrm{rd}}\\ u_{\mathrm{rq}}=d{\mathrm{\ψ}}_{\mathrm{rq}}/\mathrm{dt}-s{\mathrm{\ω}}_s{\mathrm{\ψ}}_{\mathrm{rd}}+R_r{i}_{\mathrm{rq}}\end{array}\right.---\left(3\right)$

磁链方程为

$\left\{\begin{array}{c}{\mathrm{\ψ}}_{\mathrm{sd}}=L_s{i}_{\mathrm{sd}}+L_m{i}_{\mathrm{rd}}\\ {\mathrm{\ψ}}_{\mathrm{sq}}=L_s{i}_{\mathrm{sq}}+L_m{i}_{\mathrm{rq}}\\ {\mathrm{\ψ}}_{\mathrm{rd}}=L_r{i}_{\mathrm{rd}}+L_m{i}_{\mathrm{sd}}\\ {\mathrm{\ψ}}_{\mathrm{rq}}=L_r{i}_{\mathrm{rq}}+L_m{i}_{\mathrm{sq}}\end{array}\right.---\left(4\right)$

其中：u_sd与u_sq分别为发电机定子d轴和q轴绕组电压，u_rd与u_rq分别为发电机转子d轴和q轴绕组电压，ψ_sd与ψ_sq分别为发电机定子d轴和q轴绕组磁链，ψ_rd与ψ_rq分别为发电机转子d轴和q轴绕组磁链，i_sd与i_sq分别为发电机定子d轴和q轴绕组电流，i_rd与i_rq分别为发电机转子d轴和q轴绕组电流，ω_s为坐标系旋转角速度，R_s为发电机定子绕组的电阻，R_r为发电机转子绕组的电阻，s为发电机的滑差率，L_s为发电机定子绕组的自感抗，L_r为发电机转子绕组的自感抗，L_m为发电机转子与定子绕组之间的互感抗，dψ_sd/dt与dψ_sq/dt分别为发电机定子d轴和q轴绕组磁链对时间的导数，dψ_rd/dt与dψ_rq/dt分别为发电机转子d轴和q轴绕组磁链对时间的导数；

发电机定子电压矢量方向设定为d轴，故发电机定子d轴绕组电压等于发电机定子电压矢量，发电机定子q轴绕组电压等于0，因此，双馈感应发电机输出的有功功率与无功功率为：

$\left\{\begin{array}{c}{P}_s=-3U_s{L}_m{i}_{\mathrm{rd}}/\left(2L_s\right)\\ Q_s=3(U_s^2/{\mathrm{\ω}}_s{L}_s+U_s{L}_m{i}_{\mathrm{rq}}/L_s)/2\end{array}\right.---\left(5\right)$

其中：P_s与Q_s分别为双馈感应发电机输出的有功功率与无功功率，U_s为发电机定子电压矢量，L_m为发电机转子与定子绕组之间的互感抗，i_rd与i_rq分别为发电机转子d轴和q轴绕组电流，L_s为发电机定子绕组的自感抗，ω_s为坐标系旋转角速度；

2)建立电解槽模型

电解槽阳极和阴极电极反应为：

其中：H₂O为水，O₂为氧气，H₂为氢气，H⁺为氢离子，e为电子；

阳极平衡方程式为

$\left\{\begin{array}{c}\frac{{\mathrm{dN}}_{O_2}}{\mathrm{dt}}=N_{O_2}^{\mathrm{ina}}-N_{O_2}^{\mathrm{outa}}+N_{O_2}^{\mathrm{gen}}\\ \frac{d{N}_{H_2\mathrm{Oa}}}{\mathrm{dt}}=N_{H_{2}O}^{\mathrm{ina}}-N_{H_{2}O}^{\mathrm{outa}}-N_{H_{2}O}^{\mathrm{mem}}\end{array}\right.---\left(7\right)$

其中：与分别为阳极氧气和水的摩尔量，与分别为阳极流入和流出氧气的摩尔流速，与分别为阳极流入和流出水的摩尔流速，为阳极产生的氧气的流速，为电迁移与扩散流速，与分别为阳极氧气和水的摩尔流速；

阴极平衡方程式为：

$\left\{\begin{array}{c}\frac{{\mathrm{dN}}_{H_2}}{\mathrm{dt}}=N_{H_2}^{\mathrm{inc}}-N_{H_2}^{\mathrm{outc}}+N_{H_2}^{\mathrm{gen}}\\ \frac{d{N}_{H_2\mathrm{Oc}}}{\mathrm{dt}}=N_{H_{2}O}^{\mathrm{inc}}-N_{H_{2}O}^{\mathrm{outc}}-N_{H_{2}O}^{\mathrm{mem}}\end{array}\right.---\left(8\right)$

其中：与分别为阴极氢气和水的摩尔量，与分别为阴极流入和流出氧气的摩尔流速，与分别为阴极流入和流出水的摩尔流速，为阴极产生的氢气的流速，为电迁移与扩散流速，与分别为阴极氢气和水的摩尔流速；

电解槽整体电压为：

$V_{\mathrm{ele}}=E_{\mathrm{ele}}+V_{\mathrm{ele}}^{\mathrm{act}}+V_{\mathrm{ele}}^{\mathrm{ohm}}---\left(9\right)$

式中

$\left\{\begin{array}{c}{E}_{\mathrm{ele}}=\frac{1}{2F}(\mathrm{\Δ}{G}_{\mathrm{ele}}+{\mathrm{RT}}_{\mathrm{ele}}[\ln \left(\frac{{\mathrm{\ρ}}_{H2}^{\mathrm{ele}}\sqrt{{\mathrm{\ρ}}_{O2}^{\mathrm{ele}}}}{{\mathrm{\α}}_{H2O}^{\mathrm{ele}}}\right))\\ V_{\mathrm{ele}}^{\mathrm{act}}=\frac{{\mathrm{RT}}_{\mathrm{ele}}}{2\mathrm{\βF}}\ln \left(\frac{{\stackrel{\‾}{I}}_{\mathrm{ele}}}{{\stackrel{\‾}{I}}_{\mathrm{ele}0}}\right)\\ V_{\mathrm{ele}}^{\mathrm{ohm}}={\stackrel{\‾}{I}}_{\mathrm{ele}}{R}_{\mathrm{ele}}^{\mathrm{ohm}}\end{array}\right.---\left(10\right)$

其中：V_ele为电解槽整体电压，E_ele为开路电压，为活化极化电压，为欧姆极化电压，F为法拉第常数，ΔG_ele为电化学反应过程的Gibbs自由能变，R为气体常数，T_ele为电解槽温度，为阴极氢气分压，为阳极氧气分压，为阳极和电解质之间的水活度，β为传递系数，为电流密度，为交换电流密度，为膜电阻；

3)对并网发电***的控制

风电场标准差功率为：

$P_{W-\mathrm{\δ}}=\sqrt{\frac{\underset{t-T}{\overset{t}{\∫}}{(P_W-P_{W-\mathrm{ESM}})}^2\mathrm{dt}}{T}}---\left(11\right)$

其中：P_W-δ为风电场标准差功率，P_W为风电场输出功率，P_W-ESM为经指数平滑法预测出的风电场输出较平稳的功率，T为时间间隔；

并网发电***参考功率为：

P_REF＝P_W-ESM-P_W-δ(12)

其中：P_REF为并网发电***上网参考功率，P_W-ESM为经指数平滑法预测出的风电场输出较平稳的功率，P_W-δ为风电场标准差功率；

电解槽消耗功率为：

P_C＝P_W-P_REF(13)

其中：P_C为电解槽消耗的功率，P_W为风电场输出功率，P_REF为并网发电***上网参考功率；

应用“先进先出”算法，进行电解槽优化控制，具体流程为：

A初始化电解槽开关序号、AELN和LAEL，

B若TELN大于AELN进入C，否则进入I，

C电解槽开关打开和电解槽打开次数小于电解槽持续打开最大次数，进入D，否则进入E，

D电解槽打开次数等于上时刻电解槽打开的次数加1，进入E，

E该电解槽序号达到最大，进入F，否则判断下一个电解槽返回C，

F若LAEL的下一个电解槽开关是关闭的，进入G，否则进入H，

G将电解槽开关置于开状态将电解槽打开的次数置本时刻AELN，本时刻AELN等于上时刻AELN加本时刻LAEL，并且等于当前电解槽的序号，进入H，

H若AELN达到TELN，结束，否则判断下一个电解槽返回F，

I寻找打开次数最大的电解槽，进入J，

J该电解槽打开的次数达到电解槽持续打开的最大次数，进入K，否则进入M，

K将电解槽开关置于关状态，将电解槽的打开次数归零AELN等于AELN减1，进入L

M该电解槽序号达到最大，进入N，否则判断下一个电解槽返回J，

N电解槽已打开最大次数等于电解槽已打开最大次数减1，进入O，

O若AELN达到TELN，结束，否则判断下一个电解槽返回J，

其中：AELN为某一时刻已经激活的电解槽个数，LAEL为上一次打开的电解槽的序号，TELN为制氢装置中需要安装的电解槽个数。

本发明的风电-制氢并网发电***的控制方法，它包括建立基于双馈感应发电机的风电机组模型、建立电解槽模型和对并网发电***进行控制等内容，能够充分反应对风电-制氢并网发电***建模与控制的有效性；提高了电网对新能源的接纳能力，适应性强，具有较高的实际应用价值。

附图说明

图1是风电-制氢并网发电***整体结构示意图；

图2是风电-制氢并网发电各子***之间功率关系示意图；

图3是制氢***消耗的功率与对应电解槽数；

图4是电解槽开关状态示意图；

图5是电解槽连续打开次数示意图；

图6是风电-制氢并网发电***上网功率示意图。

图7是电解槽优化控制流程框图。

图中：1风电场，2风电场汇流母线，3外部***同步发电机，4外部***高压母线，5外部电网，6风电-制氢并网发电***高压母线，7风电-制氢并网发电***中压母线，8制氢子***高压交流母线，9制氢子***低压交流母线，10制氢装置交流/直流变换器，11制氢装置直流母线，12氢电解堆模块，21风电场输出功率，22经指数平滑法预测出的风电场输出较平稳的功率，23电解槽消耗的功率，24风电场标准差功率，25并网发电***上网参考功率。

具体实施方式

本发明的风电-制氢并网发电***的控制方法，包括以下步骤：

1)建立基于双馈感应发电机的风电机组模型

风电机组空气动力学数学模型为：

$P_M={\mathrm{\ρ}}_{\mathrm{air}}{C}_p(\mathrm{\λ},\mathrm{\β})\mathrm{\π}{R}^2{V}_w^3/2---\left(1\right)$

其中：P_M为风轮机捕获的风能转化成风电机组的机械功率，ρ_air为风电接入点的空气密度，C_p为叶片的风能转换效率系数，是风轮机叶尖速比和桨距控制角的函数，λ为风轮机的叶尖速比，β为风电机组的桨距控制角，π为圆周率，一般取3.1415926，R为风轮机叶轮半径，V_w为风电接入点的时时风速；

风轮机与发电机的两质块轴系数学模型方程为：

$\left\{\begin{array}{c}2H_{T}d{\mathrm{\ω}}_T/\mathrm{dt}=T_M-K_S{\mathrm{\θ}}_S-D_T{\mathrm{\ω}}_T\\ 2H_{G}d{\mathrm{\ω}}_G/\mathrm{dt}=K_S{\mathrm{\θ}}_S-T_E-D_G{\mathrm{\ω}}_G\\ d{\mathrm{\θ}}_S/\mathrm{dt}={\mathrm{\ω}}_0({\mathrm{\ω}}_T-{\mathrm{\ω}}_G)\end{array}\right.---\left(2\right)$

其中：H_T为风轮机的惯性常数，ω_T为风轮机的电角速度，T_M为风轮机的机械转矩，K_S为风轮机和发电机轴的刚度系数，θ_S为两质块之间相对角位移，D_T为风轮机转子阻尼系数，H_G为发电机的惯性常数，ω_G为感应发电机的电角速度，T_E为发电机的电磁转矩，D_G为发电机转子阻尼系数，ω₀为电网的同步角速度，dω_T/dt为风轮机角速度对时间的导数，dω_G/dt为感应发电机电角速度对时间的导数，dθ_S/dt为两质块之间相对角位移对时间的导数；

同步旋转坐标系下双馈感应发电机的电压方程为：

$\left\{\begin{array}{c}{u}_{\mathrm{sd}}=d{\mathrm{\ψ}}_{\mathrm{sd}}/\mathrm{dt}-{\mathrm{\ω}}_s{\mathrm{\ψ}}_{\mathrm{sq}}+R_s{i}_{\mathrm{sd}}\\ u_{\mathrm{sq}}=d{\mathrm{\ψ}}_{\mathrm{sq}}/\mathrm{dt}-{\mathrm{\ω}}_s{\mathrm{\ψ}}_{\mathrm{sd}}+R_s{i}_{\mathrm{sq}}\\ u_{\mathrm{rd}}=d{\mathrm{\ψ}}_{\mathrm{rd}}/\mathrm{dt}-s{\mathrm{\ω}}_s{\mathrm{\ψ}}_{\mathrm{rq}}+R_r{i}_{\mathrm{rd}}\\ u_{\mathrm{rq}}=d{\mathrm{\ψ}}_{\mathrm{rq}}/\mathrm{dt}-s{\mathrm{\ω}}_s{\mathrm{\ψ}}_{\mathrm{rd}}+R_r{i}_{\mathrm{rq}}\end{array}\right.---\left(3\right)$

磁链方程为

$\left\{\begin{array}{c}{\mathrm{\ψ}}_{\mathrm{sd}}=L_s{i}_{\mathrm{sd}}+L_m{i}_{\mathrm{rd}}\\ {\mathrm{\ψ}}_{\mathrm{sq}}=L_s{i}_{\mathrm{sq}}+L_m{i}_{\mathrm{rq}}\\ {\mathrm{\ψ}}_{\mathrm{rd}}=L_r{i}_{\mathrm{rd}}+L_m{i}_{\mathrm{sd}}\\ {\mathrm{\ψ}}_{\mathrm{rq}}=L_r{i}_{\mathrm{rq}}+L_m{i}_{\mathrm{sq}}\end{array}\right.---\left(4\right)$

其中：u_sd与u_sq分别为发电机定子d轴和q轴绕组电压，u_rd与u_rq分别为发电机转子d轴和q轴绕组电压，ψ_sd与ψ_sq分别为发电机定子d轴和q轴绕组磁链，ψ_rd与ψ_rq分别为发电机转子d轴和q轴绕组磁链，i_sd与i_sq分别为发电机定子d轴和q轴绕组电流，i_rd与i_rq分别为发电机转子d轴和_q轴绕组电流，ω_s为坐标系旋转角速度，R_s为发电机定子绕组的电阻，R_r为发电机转子绕组的电阻，s为发电机的滑差率，L_s为发电机定子绕组的自感抗，L_r为发电机转子绕组的自感抗，L_m为发电机转子与定子绕组之间的互感抗，dψ_sd/dt与dψ_sq/dt分别为发电机定子d轴和q轴绕组磁链对时间的导数，dψ_rd/dt与dψ_rq/dt分别为发电机转子d轴和q轴绕组磁链对时间的导数；

发电机定子电压矢量方向设定为d轴，故发电机定子d轴绕组电压等于发电机定子电压矢量，发电机定子q轴绕组电压等于0，因此，双馈感应发电机输出的有功功率与无功功率为：

$\left\{\begin{array}{c}{P}_s=-3U_s{L}_m{i}_{\mathrm{rd}}/\left(2L_s\right)\\ Q_s=3(U_s^2/{\mathrm{\ω}}_s{L}_s+U_s{L}_m{i}_{\mathrm{rq}}/L_s)/2\end{array}\right.---\left(5\right)$

其中：P_s与Q_s分别为双馈感应发电机输出的有功功率与无功功率，U_s为发电机定子电压矢量，L_m为发电机转子与定子绕组之间的互感抗，i_rd与i_rq分别为发电机转子d轴和q轴绕组电流，L_s为发电机定子绕组的自感抗，ω_s为坐标系旋转角速度；

2)建立电解槽模型

电解槽阳极和阴极电极反应为：

其中：H₂O为水，O₂为氧气，H₂为氢气，H⁺为氢离子，e为电子；

阳极平衡方程式为

$\left\{\begin{array}{c}\frac{{\mathrm{dN}}_{O_2}}{\mathrm{dt}}=N_{O_2}^{\mathrm{ina}}-N_{O_2}^{\mathrm{outa}}+N_{O_2}^{\mathrm{gen}}\\ \frac{d{N}_{H_2\mathrm{Oa}}}{\mathrm{dt}}=N_{H_{2}O}^{\mathrm{ina}}-N_{H_{2}O}^{\mathrm{outa}}-N_{H_{2}O}^{\mathrm{mem}}\end{array}\right.---\left(7\right)$

其中：与分别为阳极氧气和水的摩尔量，与分别为阳极流入和流出氧气的摩尔流速，与分别为阳极流入和流出水的摩尔流速，为阳极产生的氧气的流速，为电迁移与扩散流速，与分别为阳极氧气和水的摩尔流速；

阴极平衡方程式为：

$\left\{\begin{array}{c}\frac{{\mathrm{dN}}_{H_2}}{\mathrm{dt}}=N_{H_2}^{\mathrm{inc}}-N_{H_2}^{\mathrm{outc}}+N_{H_2}^{\mathrm{gen}}\\ \frac{d{N}_{H_2\mathrm{Oc}}}{\mathrm{dt}}=N_{H_{2}O}^{\mathrm{inc}}-N_{H_{2}O}^{\mathrm{outc}}-N_{H_{2}O}^{\mathrm{mem}}\end{array}\right.---\left(8\right)$

其中：与分别为阴极氢气和水的摩尔量，与分别为阴极流入和流出氧气的摩尔流速，与分别为阴极流入和流出水的摩尔流速，为阴极产生的氢气的流速，为电迁移与扩散流速，与分别为阴极氢气和水的摩尔流速；

电解槽整体电压为：

$V_{\mathrm{ele}}=E_{\mathrm{ele}}+V_{\mathrm{ele}}^{\mathrm{act}}+V_{\mathrm{ele}}^{\mathrm{ohm}}---\left(9\right)$

式中

$\left\{\begin{array}{c}{E}_{\mathrm{ele}}=\frac{1}{2F}(\mathrm{\Δ}{G}_{\mathrm{ele}}+{\mathrm{RT}}_{\mathrm{ele}}[\ln \left(\frac{{\mathrm{\ρ}}_{H2}^{\mathrm{ele}}\sqrt{{\mathrm{\ρ}}_{O2}^{\mathrm{ele}}}}{{\mathrm{\α}}_{H2O}^{\mathrm{ele}}}\right))\\ V_{\mathrm{ele}}^{\mathrm{act}}=\frac{{\mathrm{RT}}_{\mathrm{ele}}}{2\mathrm{\βF}}\ln \left(\frac{{\stackrel{\‾}{I}}_{\mathrm{ele}}}{{\stackrel{\‾}{I}}_{\mathrm{ele}0}}\right)\\ V_{\mathrm{ele}}^{\mathrm{ohm}}={\stackrel{\‾}{I}}_{\mathrm{ele}}{R}_{\mathrm{ele}}^{\mathrm{ohm}}\end{array}\right.---\left(10\right)$

其中：V_ele为电解槽整体电压，E_ele为开路电压，为活化极化电压，为欧姆极化电压，F为法拉第常数，ΔG_ele为电化学反应过程的Gibbs自由能变，R为气体常数，T_ele为电解槽温度，为阴极氢气分压，为阳极氧气分压，为阳极和电解质之间的水活度，β为传递系数，为电流密度，为交换电流密度，为膜电阻；

3)对并网发电***的控制

风电场标准差功率为：

$P_{W-\mathrm{\δ}}=\sqrt{\frac{\underset{t-T}{\overset{t}{\∫}}{(P_W-P_{W-\mathrm{ESM}})}^2\mathrm{dt}}{T}}---\left(11\right)$

其中：P_W-δ为风电场标准差功率，P_W为风电场输出功率，P_W-ESM为经指数平滑法预测出的风电场输出较平稳的功率，T为时间间隔；

并网发电***参考功率为：

P_REF＝P_W-ESM-P_W-δ(12)

其中：P_REF为并网发电***上网参考功率，P_W-ESM为经指数平滑法预测出的风电场输出较平稳的功率，P_W-δ为风电场标准差功率；

电解槽消耗功率为：

P_C＝P_W-P_REF(13)

其中：P_C为电解槽消耗的功率，P_W为风电场输出功率，P_REF为并网发电***上网参考功率；

应用“先进先出”算法，进行电解槽优化控制，具体流程见图7所示：

A初始化电解槽开关序号、AELN和LAEL，

B若TELN大于AELN进入C，否则进入I，

C电解槽开关打开和电解槽打开次数小于电解槽持续打开最大次数，进入D，否则进入E，

D电解槽打开次数等于上时刻电解槽打开的次数加1，进入E，

E该电解槽序号达到最大，进入F，否则判断下一个电解槽返回C，

F若LAEL的下一个电解槽开关是关闭的，进入G，否则进入H，

G将电解槽开关置于开状态将电解槽打开的次数置本时刻AELN，本时刻AELN等于上时刻AELN加本时刻LAEL，并且等于当前电解槽的序号，进入H，

H若AELN达到TELN，结束，否则判断下一个电解槽返回F，

I寻找打开次数最大的电解槽，进入J，

J该电解槽打开的次数达到电解槽持续打开的最大次数，进入K，否则进入M，

K将电解槽开关置于关状态，将电解槽的打开次数归零AELN等于AELN减1，进入L

M该电解槽序号达到最大，进入N，否则判断下一个电解槽返回J，

N电解槽已打开最大次数等于电解槽已打开最大次数减1，进入O，

O若AELN达到TELN，结束，否则判断下一个电解槽返回J，

其中：AELN为某一时刻已经激活的电解槽个数，LAEL为上一次打开的电解槽的序号，TELN为制氢装置中需要安装的电解槽个数。

下面利用附图和实施例对本发明作进一步说明。

具体实例：参照图1，本发明的风电-制氢并网发电***的风电场1功率输出具有随机性、间歇性和无规律性等特点，风电场汇流母线2将风电场中的各台风机功率汇集到一起统一外送，外部***同步发电机3通过升压变压器接入外部***高压母线4并入外部电网5，制氢***通过风电-制氢并网发电***高压母线6同外部电网连接，经过一次降压变压器与风电-制氢并网发电***中压母线7，再由二次降压变压器连接制氢子***高压交流母线8和制氢子***低压交流母线9，最后通过制氢装置AC/DC变换器10和制氢装置DC母线11与氢电解堆模块12相连，实现风电-制氢***与并网发电。

参照图2，为风电-制氢并网发电各子***之间功率关系示意图，风电场输出功率21经过指数平滑计算得到经指数平滑法预测出的风电场输出较平稳的功率22，通过标准差计算获得风电场标准差功率24，进而得到并网发电***上网参考功率25，由风电场输出功率与并网发电***上网参考功率差值计算，得电解槽消耗的功率23。电解槽是制氢装置。

采用图1所示***作为算例基础，风机额定功率2MW，额定电压0.69kV，切入风速6.2m/s，切出风速21.5m/s，额定风速11m/s。电解槽额定容量0.6MW，交流/直流变换器交流额定电压0.97kV，电解槽安装数目为4个。制氢***消耗的功率与对应电解槽数如图3所示，当风电出力接近2MW时，4个电解槽全开来吸收电网无法消纳的风电弃风电量，分析可知，电解槽全开时段很少，大部分时段是在1到2个之间。电解槽开关状态与电解槽连续打开次数分别如图4和图5所示，图4与图5验证了多个电解槽的控制策略，由图4可知，4个电解槽的运行时间是平均的，当有电解槽需要打开时，总是打开已经关闭时间最长的电解槽，即最早关闭的电解槽；当有电解槽需要关闭时，总是关闭已经运行时间最长的电解槽，即最先打开的电解槽。图5同样说明这个道理，图中的纵坐标表示电解槽持续打开的次数，电解槽持续打开的次数越大，表示其连续运行的时间最长，所以当有电解槽需要被关闭时，总是关闭纵坐标最大的那个电解槽。这样平衡了每个电解槽的运行时间。风电-制氢并网发电***上网功率如图6所示，图6中将大于参考功率之上的部分风电输出功率输送到制氢装置***中，将其用来作为电解水制氢的能源。因此，经过仿真验证表明，这种风电-制氢并网发电系可有效地降低风电场输出有功功率的波动，同时获得大量清洁的氢能源，风电-制氢并网发电***的控制方法是高效且实用的。

Claims (1)

1.一种风电-制氢并网发电***的控制方法，其特征是，它包括以下步骤：

1)建立基于双馈感应发电机的风电机组模型

风电机组空气动力学数学模型为：

$P_M={\mathrm{\ρ}}_{air}{C}_p(\mathrm{\λ},\mathrm{\β}){\mathrm{\πR}}^2{V}_w^3/2---\left(1\right)$

其中：P_M为风轮机捕获的风能转化成风电机组的机械功率，ρ_air为风电接入点的空气密度，C_p为叶片的风能转换效率系数，是风轮机叶尖速比和桨距控制角的函数，λ为风轮机的叶尖速比，β为风电机组的桨距控制角，π为圆周率，一般取3.1415926，R为风轮机叶轮半径，V_w为风电接入点的实时风速；

风轮机与发电机的两质块轴系数学模型方程为：

$\left\{\begin{array}{c}2H_T{\mathrm{d\ω}}_T/dt=T_M-K_S{\mathrm{\θ}}_S-D_T{\mathrm{\ω}}_T\\ 2H_G{\mathrm{d\ω}}_G/dt=K_S{\mathrm{\θ}}_S-T_E-D_G{\mathrm{\ω}}_G\\ {\mathrm{d\θ}}_S/dt={\mathrm{\ω}}_0({\mathrm{\ω}}_T-{\mathrm{\ω}}_G)\end{array}\right.---\left(2\right)$

其中：H_T为风轮机的惯性常数，ω_T为风轮机的电角速度，T_M为风轮机的机械转矩，K_S为风轮机和发电机轴的刚度系数，θ_S为两质块之间相对角位移，D_T为风轮机转子阻尼系数，H_G为发电机的惯性常数，ω_G为感应发电机的电角速度，T_E为发电机的电磁转矩，D_G为发电机转子阻尼系数，ω₀为电网的同步角速度，dω_T/dt为风轮机角速度对时间的导数，dω_G/dt为感应发电机电角速度对时间的导数，dθ_S/dt为两质块之间相对角位移对时间的导数；

同步旋转坐标系下双馈感应发电机的电压方程为：

$\left\{\begin{array}{c}{u}_{sd}={\mathrm{d\ψ}}_{sd}/dt-{\mathrm{\ω}}_s{\mathrm{\ψ}}_{sq}+R_s{i}_{sd}\\ u_{sq}={\mathrm{d\ψ}}_{sq}/dt-{\mathrm{\ω}}_s{\mathrm{\ψ}}_{sd}+R_s{i}_{sq}\\ u_{rd}={\mathrm{d\ψ}}_{rd}/dt-{\mathrm{s\ω}}_s{\mathrm{\ψ}}_{rq}+R_r{i}_{rd}\\ u_{rq}={\mathrm{d\ψ}}_{rq}/dt-{\mathrm{s\ω}}_s{\mathrm{\ψ}}_{rd}+R_r{i}_{rq}\end{array}\right.---\left(3\right)$

磁链方程为

$\left\{\begin{array}{c}{\mathrm{\ψ}}_{sd}=L_s{i}_{sd}+L_m{i}_{rd}\\ {\mathrm{\ψ}}_{sq}=L_s{i}_{sq}+L_m{i}_{rq}\\ {\mathrm{\ψ}}_{rd}=L_r{i}_{rd}+L_m{i}_{sd}\\ {\mathrm{\ψ}}_{rq}=L_r{i}_{rq}+L_m{i}_{sq}\end{array}\right.---\left(4\right)$

其中：u_sd与u_sq分别为发电机定子d轴和q轴绕组电压，u_rd与u_rq分别为发电机转子d轴和q轴绕组电压，ψ_sd与ψ_sq分别为发电机定子d轴和q轴绕组磁链，ψ_rd与ψ_rq分别为发电机转子d轴和q轴绕组磁链，i_sd与i_sq分别为发电机定子d轴和q轴绕组电流，i_rd与i_rq分别为发电机转子d轴和q轴绕组电流，ω_s为坐标系旋转角速度，R_s为发电机定子绕组的电阻，R_r为发电机转子绕组的电阻，s为发电机的滑差率，L_s为发电机定子绕组的自感抗，L_r为发电机转子绕组的自感抗，L_m为发电机转子与定子绕组之间的互感抗，dψ_sd/dt与dψ_sq/dt分别为发电机定子d轴和q轴绕组磁链对时间的导数，dψ_rd/dt与dψ_rq/dt分别为发电机转子d轴和q轴绕组磁链对时间的导数；

发电机定子电压矢量方向设定为d轴，故发电机定子d轴绕组电压等于发电机定子电压矢量，发电机定子q轴绕组电压等于0，因此，双馈感应发电机输出的有功功率与无功功率为：

$\left\{\begin{array}{c}{P}_s=-3U_s{L}_m{i}_{rd}/\left(2L_s\right)\\ Q_s=3(U_s^2/{\mathrm{\ω}}_s{L}_s+U_s{L}_m{i}_{rq}/L_s)/2\end{array}\right.---\left(5\right)$

其中：P_s与Q_s分别为双馈感应发电机输出的有功功率与无功功率，U_s为发电机定子电压矢量，L_m为发电机转子与定子绕组之间的互感抗，i_rd与i_rq分别为发电机转子d轴和q轴绕组电流，L_s为发电机定子绕组的自感抗，ω_s为坐标系旋转角速度；

2)建立电解槽模型

电解槽阳极和阴极电极反应为：

其中：H₂O为水，O₂为氧气，H₂为氢气，H⁺为氢离子，e为电子；

阳极平衡方程式为

$\left\{\begin{array}{c}\frac{{\mathrm{dN}}_{O_2}}{dt}=N_{O_2}^{ina}-N_{O_2}^{outa}+N_{O_2}^{gen}\\ \frac{{\mathrm{dN}}_{H_{2}Oa}}{dt}=N_{H_{2}O}^{ina}-N_{H_{2}O}^{outa}-N_{H_{2}O}^{mem}\end{array}\right.---\left(7\right)$

其中：与分别为阳极氧气和水的摩尔量，与分别为阳极流入和流出氧气的摩尔流速，与分别为阳极流入和流出水的摩尔流速，为阳极产生的氧气的流速，为电迁移与扩散流速，与分别为阳极氧气和水的摩尔流速；

阴极平衡方程式为：

$\left\{\begin{array}{c}\frac{{\mathrm{dN}}_{H_2}}{dt}=N_{H_2}^{inc}-N_{H_2}^{outc}+N_{H_2}^{gen}\\ \frac{{\mathrm{dN}}_{H_{2}Oc}}{dt}=N_{H_{2}O}^{inc}-N_{H_{2}O}^{outc}+N_{H_{2}O}^{mem}\end{array}\right.---\left(8\right)$

其中：与分别为阴极氢气和水的摩尔量，与分别为阴极流入和流出氧气的摩尔流速，与分别为阴极流入和流出水的摩尔流速，为阴极产生的氢气的流速，为电迁移与扩散流速，与分别为阴极氢气和水的摩尔流速；

电解槽整体电压为：

$V_{ele}=E_{ele}+V_{ele}^{act}+V_{ele}^{ohm}---\left(9\right)$

式中

$\left\{\begin{array}{c}{E}_{ele}=\frac{1}{2F}\left({\mathrm{\ΔG}}_{ele}+{\mathrm{RT}}_{ele}\[\ln \left(\frac{{\mathrm{\ρ}}_{H2}^{ele}\sqrt{{\mathrm{\ρ}}_{O2}^{ele}}}{{\mathrm{\α}}_{H2O}^{ele}}\right)\]\right)\\ V_{ele}^{act}=\frac{{\mathrm{RT}}_{ele}}{2\mathrm{\β}F}\ln \left(\frac{{\stackrel{\‾}{I}}_{ele}}{{\stackrel{\‾}{I}}_{ele0}}\right)\\ V_{ele}^{ohm}={\stackrel{\‾}{I}}_{ele}{R}_{ele}^{ohm}\end{array}\right.---\left(10\right)$

其中：V_ele为电解槽整体电压，E_ele为开路电压，为活化极化电压，为欧姆极化电压，F为法拉第常数，ΔG_ele为电化学反应过程的Gibbs自由能变，R为气体常数，T_ele为电解槽温度，为阴极氢气分压，为阳极氧气分压，为阳极和电解质之间的水活度，β为传递系数，为电流密度，为交换电流密度，为膜电阻；

3)对并网发电***的控制

风电场标准差功率为：

$P_{W-\mathrm{\δ}}=\sqrt{\frac{\underset{t-T}{\overset{t}{\∫}}{(P_W-P_{W-ESM})}^{2}dt}{T}}---\left(11\right)$

其中：P_W-δ为风电场标准差功率，P_W为风电场输出功率，P_W-ESM为经指数平滑法预测出的风电场输出较平稳的功率，T为时间间隔；

并网发电***参考功率为：

P_REF＝P_W-ESM-P_W-δ(12)

其中：P_REF为并网发电***上网参考功率，P_W-ESM为经指数平滑法预测出的风电场输出较平稳的功率，P_W-δ为风电场标准差功率；

电解槽消耗功率为：

P_C＝P_W-P_REF(13)

其中：P_C为电解槽消耗的功率，P_W为风电场输出功率，P_REF为并网发电***上网参考功率；

应用“先进先出”算法，进行电解槽优化控制，具体流程为：

(A)初始化电解槽开关序号、AELN和LAEL，

(B)当TELN大于AELN时，进入(C)，否则进入(I)，

(C)当满足电解槽打开，并且电解槽打开次数小于电解槽持续打开最大次数时，进入(D)，否则进入(E)，

(D)本时刻电解槽打开的次数，为上一时刻电解槽打开的次数加1，进入(E)，

(E)该电解槽序号达到最大，进入(F)，否则判断下一个电解槽返回(C)，

(F)LAEL的下一个电解槽开关是关闭的，进入(G)，否则进入(H)，

(G)将电解槽开关打开，同时将电解槽打开的次数存入AELN，本时刻AELN等于上时刻AELN加本时刻LAEL，并且等于当前电解槽的序号，进入(H)，

(H)AELN达到TELN，结束，否则判断下一个电解槽返回(F)，

(I)寻找打开次数最大的电解槽，进入(J)，

(J)当该电解槽打开的次数，达到电解槽持续打开的最大次数时，进入(K)，否则进入(M)，

(K)将电解槽开关置于关状态，将电解槽的打开次数归零AELN等于AELN减1，进入(L)，

(M)该电解槽序号达到最大，进入(N)，否则判断下一个电解槽返回(J)，

(N)当本时刻电解槽已打开的最大次数，为上一时刻电解槽已打开的最大次数减1时，进入(O)，

(O)AELN达到TELN，结束，否则判断下一个电解槽返回(J)，

其中：AELN为某一时刻已经激活的电解槽个数，LAEL为上一次打开的电解槽的序号，TELN为制氢装置中需要安装的电解槽个数。