光伏发电系统在电网故障时面临的核心挑战是电压骤降情况下的持续稳定运行。传统光伏逆变器在检测到电网电压异常时通常会采取保护性脱网措施,但这种做法在现代电力系统中已无法满足要求。以德国中压电网导则为例,明确规定当电网电压跌落到额定值的15%时,光伏电站仍需保持并网运行至少150毫秒。这种技术能力就是我们所说的低电压穿越(Low Voltage Ride Through, LVRT)功能。
两级式结构(DC-DC Boost+逆变器)因其灵活的电压适配能力和高效率特性,成为中功率光伏系统的首选拓扑。前级Boost电路负责将光伏阵列的直流电压提升到适合逆变并网的母线电压,后级逆变器则完成直流到交流的转换。这种结构在LVRT工况下展现出独特优势:Boost级可以通过调节占空比快速稳定直流母线电压,而逆变级则专注于网侧电流的精确控制。
完整的LVRT控制体系包含三个关键控制环:最外层的直流电压环(仅在正常模式下工作)、中间层的功率分配环(负责有功/无功电流参考值生成)以及最内层的电流跟踪环。当检测到电网电压跌落时,系统立即转入LVRT模式,此时电压环被禁用,功率分配环按照预设的LVRT曲线生成电流指令。
典型的LVRT曲线要求:
传统扰动观察法在LVRT期间会导致直流母线电压失控的根本原因在于其单边追求最大功率点的控制逻辑。我们提出的改进方案通过引入动态限幅机制来解决这个问题:
matlab复制function [Vref, Pmax] = AdvancedMPPT(Vpv, Ipv, Vdc, Vdc_max)
persistent Vprev Pprev dV direction;
% 初始化
if isempty(Vprev)
Vprev = Vpv;
Pprev = Vpv*Ipv;
dV = 0.5; % 初始扰动步长
direction = 1;
end
Pnow = Vpv*Ipv;
deltaP = Pnow - Pprev;
% LVRT模式判断
if Vdc > 0.95*Vdc_max
% 正常MPPT模式
if deltaP > 0
direction = sign(Vpv - Vprev);
else
direction = -sign(Vpv - Vprev);
end
Vref = Vpv + direction*dV;
else
% LVRT保护模式
Vref = min(Vpv, 0.8*Vdc_max);
dV = 0.1; % 减小扰动步长
end
% 更新历史值
Vprev = Vpv;
Pprev = Pnow;
Pmax = Pnow;
end
该算法具有两个显著特点:
二阶广义积分器(DSOGI)锁相环相比传统SRF-PLL在谐波抑制方面具有明显优势。我们采用的改进结构包含以下关键环节:
正交信号生成器:
matlab复制function [v_alpha, v_beta] = SOGI(v_in, w, k)
persistent integrator1 integrator2;
% 实现状态空间方程
dv_alpha = w*k*(v_in - v_alpha) - w*v_beta;
dv_beta = w*v_alpha;
v_alpha = integrator1.update(dv_alpha);
v_beta = integrator2.update(dv_beta);
end
序列分离网络:
matlab复制function [v_pos, v_neg] = SequenceSeparation(v_alpha, v_beta)
v_pos = 0.5*(v_alpha - 1j*v_beta);
v_neg = 0.5*(v_alpha + 1j*v_beta);
end
频率自适应机制:
matlab复制function w_est = FrequencyAdaptation(v_alpha, v_beta, w_nom, k_gamma)
error = v_alpha*cos(w_nom*t) - v_beta*sin(w_nom*t);
w_est = w_nom + k_gamma*error;
end
实测表明,这种结构在电压畸变率15%的情况下仍能保持相位误差小于1度,完全满足LVRT要求。
设计参数示例(2000W系统):
电感计算:
matlab复制function L = BoostInductorDesign(Vin_min, Vout, Pmax, fs, ripple_ratio)
D = 1 - Vin_min/Vout; % 最大占空比
Iavg = Pmax/Vin_min;
deltaI = ripple_ratio*Iavg;
L = Vin_min*D/(fs*deltaI);
end
计算得L=2.5mH(取标准值2.2mH)
输出电容选择:
matlab复制function C = BoostCapacitorDesign(Pmax, Vdc, hold_up_time, deltaV_max)
E = Pmax*hold_up_time;
C = 2*E/(Vdc^2 - (Vdc - deltaV_max)^2);
end
按保持时间10ms、电压跌落不超过5%计算得C=680μF
采用三阶Butterworth特性设计:
总电感量限制:
matlab复制function Ltotal = LCL_TotalInductance(Vgrid, Pn, fs, ripple_max)
Ipeak = sqrt(2)*Pn/Vgrid;
Ltotal = Vgrid/(sqrt(2)*Ipeak*fs*ripple_max);
end
计算得Ltotal=3.2mH
电容取值约束:
matlab复制function Cmax = LCL_CapacitanceLimit(Pn, Vgrid, fn, reactive_ratio)
Qmax = reactive_ratio*Pn;
Cmax = Qmax/(3*2*pi*fn*Vgrid^2);
end
按无功不超过2.5%计算得C<6.8μF
谐振频率验证:
matlab复制function fres = LCL_ResonantFrequency(L1, L2, C)
fres = 1/(2*pi)*sqrt((L1+L2)/(L1*L2*C));
end
应满足:10fn < fres < 0.5fs
电网电压前馈的核心是消除电网扰动对电流环的影响。考虑逆变器输出电压方程:
Vinv = Vgrid + (sL + R)*Igrid
前馈项应补偿Vgrid波动,因此前馈传递函数为:
Gff(s) = 1/(sL + R)
离散化实现(Tustin变换):
matlab复制function v_ff = CurrentFeedforward(v_grid, L, R, Ts)
persistent z1;
% 双线性变换离散化
a0 = 2*L + R*Ts;
b0 = Ts;
b1 = Ts;
v_ff = (b0*v_grid + b1*z1)/a0;
z1 = v_grid;
end
采用复系数PI调节器实现正负序统一控制:
matlab复制function [v_d, v_q] = ComplexPIController(i_d_ref, i_q_ref, i_d, i_q, Kp, Ki, Ts)
persistent err_d err_q int_d int_q;
% 误差计算
err_d = i_d_ref - i_d;
err_q = i_q_ref - i_q;
% 积分项
int_d = int_d + Ki*Ts*err_d;
int_q = int_q + Ki*Ts*err_q;
% 输出限幅
int_d = min(max(int_d, -1), 1);
int_q = min(max(int_q, -1), 1);
% 控制输出
v_d = Kp*err_d + int_d;
v_q = Kp*err_q + int_q;
end
参数整定经验:
标准电压跌落测试波形生成:
matlab复制function v_grid = GenerateLVRTWaveform(t, Vn, dip_ratio, dip_duration, fault_type)
% fault_type: 1-对称跌落, 2-不对称跌落
if t < 0.1 || t > 0.1 + dip_duration
v_grid = Vn*sin(2*pi*50*t);
else
if fault_type == 1
v_grid = Vn*dip_ratio*sin(2*pi*50*t);
else
% 典型两相跌落
v_grid = Vn*[sin(2*pi*50*t);
dip_ratio*sin(2*pi*50*t - 2*pi/3);
sin(2*pi*50*t + 2*pi/3)];
end
end
end
实测数据示例:
实际系统中存在的计算延迟(通常1-2个PWM周期)会显著影响高频性能。补偿方法:
预测控制:
matlab复制function i_pred = CurrentPredictor(i_meas, v_inv, v_grid, L, Ts, delay)
di = (v_inv - v_grid - R*i_meas)/L;
i_pred = i_meas + di*(delay*Ts);
end
状态观测器:
matlab复制function x_hat = CurrentObserver(u, y, A, B, C, L, Ts)
persistent x;
% 预测步
x = A*x + B*u;
% 校正步
x = x + L*(y - C*x);
x_hat = x;
end
多级保护协调策略:
第一级:软件保护(电流环限幅)
第二级:硬件保护(比较器)
第三级:熔断保护
保护复位策略: