1. 分布式电源接入配电网的影响分析
1.1 配电网结构特点与运行挑战
现代配电网最初设计为单向供电网络,其拓扑结构、保护系统和电压调节机制都是基于传统集中式发电模式构建的。典型的10kV配电网采用辐射状结构,通过变电站变压器将110kV或35kV电压降至10kV,再经配电变压器降至380V/220V供给终端用户。这种结构下,潮流方向固定,短路电流水平相对稳定,保护整定值设置简单明确。
当分布式电源(如双馈感应发电机)接入配电网后,最直接的改变是潮流方向可能发生逆转。以光伏电站为例,在阳光充足时段,光伏系统的出力可能超过本地负荷需求,导致功率反送至上级电网。这种双向潮流的出现会带来三个主要问题:
-
电压调节困难:传统配电网依靠有载调压变压器和并联电容器组进行电压调节,这些设备基于单向潮流的假设设计。当分布式电源注入功率时,线路电压分布特性改变,可能导致某些节点电压越限。
-
保护系统误动作:过电流保护装置的动作值通常按最大短路电流整定。分布式电源的接入会改变短路电流的大小和方向,可能导致保护拒动或误动。特别是当主系统侧发生故障时,分布式电源的持续馈入会维持故障电流,使保护无法正确识别故障区域。
-
电能质量恶化:风力发电等间歇性电源的输出功率波动会引发电压闪变和频率偏差。双馈电机的电力电子接口还会产生谐波污染,影响敏感负荷的正常运行。
1.2 双馈感应发电机的工作原理
双馈感应发电机(DFIG)是目前风力发电的主流机型,其独特之处在于转子侧通过背靠背变流器与电网连接。这种结构使得发电机可以在同步转速±30%的范围内运行,大幅提高了风能捕获效率。
DFIG的定子直接并网,转子通过滑环连接变流器。变流器由网侧换流器(GSC)和转子侧换流器(RSC)组成,中间通过直流母线连接。RSC控制转子电流的频率和相位,实现最大功率点跟踪(MPPT);GSC维持直流母线电压稳定,并调节无功功率输出。
当风速变化导致机械转速偏离同步速时,转子电流频率自动调整,保持定子输出频率与电网同步。例如,当转速低于同步速时(亚同步运行),变流器向转子注入正序电流;当转速高于同步速时(超同步运行),则从转子吸收能量。这种灵活的能量双向流动特性,使DFIG特别适合需要宽转速范围的风力发电应用。
1.3 仿真建模的关键参数
建立准确的DFIG模型需要考虑以下核心参数:
-
电机参数:
- 定子电阻(Rs)和漏感(Lls)
- 转子电阻(Rr')和漏感(Llr')
- 互感(Lm)
- 极对数(p)
-
变流器参数:
- 直流母线电压(Vdc)
- 开关频率(fsw)
- 滤波器电感(Lf)和电容(Cf)
-
控制参数:
- 外环电压/功率控制器的PI参数
- 内环电流控制器的PI参数
- 锁相环(PLL)带宽
这些参数需要通过电机铭牌数据或实测值确定。例如,一台1.5MW的DFIG典型参数可能为:Rs=0.023pu, Rr'=0.016pu, Lm=4.5pu, Lls=Llr'=0.18pu。在参数不确定时,可参考同容量电机的典型值进行估算。
2. DFIG建模与仿真实现
2.1 MATLAB/Simulink建模框架
在Simulink中搭建DFIG模型通常包含以下几个子系统:
-
风力机模型:
matlab复制% 风能捕获计算公式 P_m = 0.5 * rho * A * Cp * v_wind^3; T_m = P_m / (w_m * pole_pairs);其中Cp为风能利用系数,通常用如下近似表达式:
matlab复制Cp = 0.22*(116/lambda - 0.4*beta - 5)*exp(-12.5/lambda) lambda = (blade_radius * w_m)/v_wind; -
轴系模型:
采用双质块模型表示风机与发电机之间的传动系统:matlab复制% 运动方程 H_t * dw_t/dt = T_m - K_s * (theta_t - theta_g) - D_s * (w_t - w_g) H_g * dw_g/dt = K_s * (theta_t - theta_g) + D_s * (w_t - w_g) - T_e -
DFIG电气模型:
在同步旋转dq坐标系下建立电压方程:matlab复制% 定子电压方程 Vds = Rs*Ids + d(psi_ds)/dt - omega_s*psi_qs; Vqs = Rs*Iqs + d(psi_qs)/dt + omega_s*psi_ds; % 转子电压方程 Vdr = Rr*Idr + d(psi_dr)/dt - (omega_s - omega_r)*psi_qr; Vqr = Rr*Iqr + d(psi_qr)/dt + (omega_s - omega_r)*psi_dr; -
变流器控制系统:
采用典型的双闭环控制结构:- 外环:直流电压控制(GSC)、有功无功控制(RSC)
- 内环:电流跟踪控制
2.2 关键控制策略实现
转子侧换流器(RSC)控制策略:
matlab复制% 有功功率控制环(通过转速调节)
omega_ref = interp1(power_table, speed_table, P_ref);
error_omega = omega_ref - omega_m;
Id_ref = Kp_omega * error_omega + Ki_omega * integral(error_omega);
% 无功功率控制环
error_Q = Q_ref - Q_meas;
Iq_ref = Kp_Q * error_Q + Ki_Q * integral(error_Q);
% 电流控制环
Vd_inner = Kp_i * (Id_ref - Id_meas) + Ki_i * integral(Id_ref - Id_meas) - omega_slip * Lsigma * Iq_meas;
Vq_inner = Kp_i * (Iq_ref - Iq_meas) + Ki_i * integral(Iq_ref - Iq_meas) + omega_slip * (Lsigma * Id_meas + Lm^2/Ls * Ims);
网侧换流器(GSC)控制策略:
matlab复制% 直流电压控制环
error_Vdc = Vdc_ref - Vdc_meas;
Id_ref = Kp_Vdc * error_Vdc + Ki_Vdc * integral(error_Vdc);
% 交流电压控制环
error_Vac = Vac_ref - Vac_meas;
Iq_ref = Kp_Vac * error_Vac + Ki_Vac * integral(error_Vac);
% 电流控制环
Vd_inner = Kp_i * (Id_ref - Id_meas) + Ki_i * integral(Id_ref - Id_meas) - omega_s * Lf * Iq_meas;
Vq_inner = Kp_i * (Iq_ref - Iq_meas) + Ki_i * integral(Iq_ref - Iq_meas) + omega_s * Lf * Id_meas;
2.3 仿真场景设置
典型的仿真案例应包括以下测试场景:
-
风速阶跃变化:
- 初始风速8m/s,5秒时跃升至12m/s
- 观察最大功率点跟踪效果和转速调节过程
-
电网电压跌落:
- 在3秒时施加30%的电压跌落,持续500ms
- 验证低电压穿越(LVRT)能力
-
无功功率调节:
- 设定Q_ref从0到0.4pu阶跃变化
- 检查动态响应时间和超调量
仿真时间步长选择至关重要。对于包含电力电子器件的系统,建议采用50μs或更小的步长以确保开关过程的准确模拟。对于仅关注机电动态的过程,可采用1ms左右的步长以提高仿真效率。
3. 仿真结果分析与问题排查
3.1 典型波形解读
-
正常并网运行波形:
- 定子电流THD应小于5%
- 直流母线电压波动范围±2%
- 功率因数控制在设定值±0.01范围内
-
故障穿越波形:
- 电压跌落期间,转子电流不应超过2pu
- 故障清除后,有功功率应在1秒内恢复90%以上
- 整个过程中直流电压不超过1.15倍额定值
-
功率调节响应:
- 有功功率阶跃响应的调节时间应小于0.5秒
- 无功功率阶跃响应的超调量应小于10%
3.2 常见问题与解决方案
-
仿真发散问题:
- 现象:仿真中途报错停止或波形异常发散
- 可能原因:
- 初始条件设置不当(如电机初始滑差过大)
- 控制参数过于激进导致振荡
- 代数环(Algebraic Loop)问题
- 解决方案:
- 检查所有积分器的初始状态
- 降低PI控制器的比例系数
- 在适当位置插入Unit Delay模块
-
电流波形畸变:
- 现象:定子电流含有明显谐波
- 可能原因:
- PWM开关频率设置过低(建议≥2kHz)
- 电流控制器带宽不足
- 电网阻抗设置不合理
- 解决方案:
- 提高变流器开关频率
- 重新设计电流控制器(带宽取1/10开关频率)
- 在电网侧添加适当阻抗
-
功率控制响应迟缓:
- 现象:功率指令跟踪速度慢
- 可能原因:
- 外环PI参数整定不当
- 内环电流限幅值设置过小
- 机械惯性时间常数过大
- 解决方案:
- 采用Ziegler-Nichols法重新整定PI参数
- 检查电流限幅值(通常设为1.2-1.5倍额定)
- 验证轴系参数准确性
3.3 模型验证技巧
-
稳态验证:
- 在额定工况下,检查以下关系是否成立:
matlab复制abs(P_mech - P_elec)/P_rated < 0.01 abs(Q_meas - Q_ref)/S_rated < 0.005
- 在额定工况下,检查以下关系是否成立:
-
动态验证:
- 对比阶跃响应与理论计算值:
matlab复制rise_time_theoretical = 2.2/(bandwidth * 2*pi); overshoot_theoretical = exp(-zeta*pi/sqrt(1-zeta^2))*100;
- 对比阶跃响应与理论计算值:
-
参数敏感性分析:
- 对关键参数(如Lm、Rr)进行±10%的扰动,观察输出变化:
matlab复制
sensitivity = (delta_output/output)/(delta_parameter/parameter);
- 对关键参数(如Lm、Rr)进行±10%的扰动,观察输出变化:
在实际项目中,我们通常会先用简化模型验证控制策略的正确性,再逐步增加模型复杂度。例如,可以先使用理想电源代替风力机模型,待电气部分调试完成后再接入完整的机械系统。这种分阶段验证方法能有效定位问题源头。
