1. 风光混合储能系统概述
三机并联的风光混合储能系统是当前新能源并网领域的热门研究方向，它巧妙地将光伏发电、风力发电和储能装置结合在一起，形成一个稳定可靠的微电网系统。作为一名电力电子工程师，我在过去三年里搭建过不下二十种不同配置的风光储系统，今天要分享的这个三机并联方案，可以说是集实用性和教学性于一体的典型案例。
这个系统的核心价值在于解决了新能源发电的两个致命问题：间歇性和波动性。光伏发电完全依赖光照强度，一片云飘过功率就可能掉一半；风机输出功率更是随风速的三次方变化，说变就变。而超级电容和电池组成的混合储能系统，就像给整个系统装上了"减震器"，能够快速平抑功率波动。
2. 光伏阵列建模与MPPT实现
2.1 光伏电池数学模型解析
光伏电池的电气特性可以用单二极管模型精确描述，这个模型本质上是一个非线性电流源。我常用的Matlab实现代码如下：
matlab复制function I = PVcell(V, G, T)
q = 1.6e-19; % 电子电荷量
k = 1.38e-23; % 玻尔兹曼常数
Rs = 0.05; % 串联电阻(Ω)
Rsh = 100; % 并联电阻(Ω)
Iph = G*0.75; % 光生电流(A)，0.75是转换系数
Irs = 1e-6*((T/298)^3)*exp(1.12*q/(k*T)); % 反向饱和电流
Vt = k*T/q; % 热电压(V)
I = Iph - Irs*(exp((V+Rs*I)/Vt)-1) - (V + Rs*I)/Rsh;
end
几个关键参数需要特别注意：
- 温度T的影响：温度每升高1℃，开路电压下降约0.4%，这个特性在做高温环境仿真时不能忽略
- 光照强度G：直接决定光生电流Iph，但要注意实际应用中存在"效率衰减"现象
- 串联电阻Rs：这个参数虽小，但在大电流工况下会引发明显的功率损耗
实际工程经验：在搭建实物系统时，建议先用厂家提供的参数表校准模型，特别是Rs和Rsh这两个参数，不同厂家的电池片差异可能达到±15%
2.2 扰动观察法MPPT实现
最大功率点跟踪(MPPT)是光伏系统的核心技术，这里展示的扰动观察法(P&O)虽然简单，但在80%的实际应用中都是首选方案：
matlab复制duty_cycle = 0.5; % 初始占空比
step = 0.01; % 扰动步长
prev_power = 0; % 上一周期功率
for t = 1:100
V = get_voltage(); % 获取当前电压
I = get_current(); % 获取当前电流
current_power = V*I; % 计算瞬时功率
if current_power > prev_power
duty_cycle = duty_cycle + step; % 功率增加，保持方向
else
step = -step; % 功率减小，改变方向
duty_cycle = duty_cycle + step;
end
prev_power = current_power;
end
这个算法有几点实操建议：
- 步长选择：一般取最大功率点处dP/dV斜率的1/5到1/10，太大会振荡，太小响应慢
- 采样周期：建议为10-100ms，需要与DC-DC变换器的开关周期匹配
- 死区设置：当功率变化小于1%时，可以暂停扰动以避免不必要的振荡
3. 永磁同步风机控制策略
3.1 风机MPPT的特殊性
风机的MPPT原理与光伏类似，但有两个关键差异点：
- 控制对象是转速而非电压，需要先通过叶尖速比(TSR)换算
- 机械惯性大，响应速度比光伏系统慢10-100倍
典型的风机MPPT实现代码：
matlab复制function omega_opt = WindMPPT(P_mech)
persistent prev_P delta_omega;
if isempty(prev_P)
prev_P = 0;
delta_omega = 0.1; % 转速扰动步长(rad/s)
end
if P_mech > prev_P
omega_opt = omega_opt + delta_omega;
else
delta_omega = -delta_omega;
omega_opt = omega_opt + delta_omega;
end
prev_P = P_mech;
end
3.2 矢量控制实现要点
永磁同步电机的矢量控制需要特别注意：
- 锁相环(PLL)的带宽设置：一般为电网频率的1/10
- dq轴解耦：电流环带宽至少是速度环的5倍
- 过调制处理：当调制比超过1.15时，需要启用过调制算法
4. 混合储能系统设计
4.1 超级电容与电池的配合
混合储能系统的核心思想是：
- 超级电容应对高频、大功率的瞬时波动
- 电池处理低频、长时间的功率调节
这种组合的典型效率对比：
| 储能类型 | 充放电效率 | 响应时间 | 循环寿命 |
|---|---|---|---|
| 超级电容 | 95%-98% | <10ms | >50万次 |
| 锂电池 | 85%-93% | 50-200ms | 3000-5000次 |
4.2 状态机控制实现
matlab复制switch mode
case 'CHARGE'
if V_sc > 300 % 电压上限保护
mode = 'HOLD';
elseif I_grid > 1.2*I_nom % 过载情况
I_sc = (I_grid - I_nom)/0.95; % 考虑充电效率
end
case 'DISCHARGE'
if V_sc < 200 % 电压下限保护
mode = 'CHARGE';
elseif I_grid < 0.8*I_nom % 欠载情况
I_sc = (I_nom - I_grid)/0.9; % 考虑放电效率
end
case 'HOLD'
if abs(I_grid - I_nom) > 0.1*I_nom % 超出死区
mode = 'CHARGE';
end
end
5. 并网控制策略
5.1 PQ控制双环设计
matlab复制% 电压外环
V_ref = 311*sin(2*pi*50*t); % 50Hz电网参考电压
V_error = V_ref - V_actual;
I_ref = kp_v*V_error + ki_v*integral(V_error);
% 电流内环
I_error = I_ref - I_actual;
duty = kp_i*I_error + ki_i*integral(I_error);
参数整定经验：
- 电流环带宽至少是电压环的5倍
- 典型取值：kp_v=0.3-0.6，ki_v=5-10；kp_i=3-8，ki_i=50-100
- 遇到振荡时，先减小积分系数ki，再调整比例系数kp
5.2 三机并联的注意事项
- 环流抑制：各单元输出端需要加装均流电抗器
- 同步问题：建议采用主从控制结构，避免多机竞争
- 保护协调：过流保护阈值需要根据并联数量重新计算
6. 仿真与调试技巧
6.1 典型测试场景
- 光照阶跃变化：从1000W/m²突降到600W/m²
- 风速渐变：3m/s线性增加到9m/s
- 负载投切：50%负载突然增加到100%
6.2 波形分析要点
- 母线电压波动：正常应保持在±5%以内
- 频率偏差：并网模式下应小于0.2Hz
- 谐波含量：THD<3%为合格
调试心得：遇到不收敛的情况，先检查各个子系统的初始状态是否合理，特别是PLL的初始相位和DC bus的初始电压。有时候把仿真步长从auto改为固定步长(如50us)也能解决奇怪的问题。