如果你正在研究光储并网系统,特别是关注如何提升系统在弱电网条件下的稳定性和电能质量,那么“谐波含量闭环调控”这个技术点很可能已经进入了你的视野。传统的光储虚拟同步发电机(VSG)控制策略主要关注基波功率的模拟与支撑,但在实际电网中,尤其是存在大量非线性负载和分布式电源的场合,谐波问题日益突出。谐波不仅会降低电能质量,还可能引发谐振,威胁系统安全。
一个常见的误区是,认为只要加入了LCL滤波器就能解决所有谐波问题。实际上,电网阻抗并非固定不变,特别是在弱电网场景下,电网阻抗的宽范围变化会显著影响LCL滤波器的谐振特性,传统的固定参数虚拟阻抗策略往往力不从心。这时,系统可能面临谐波放大甚至失稳的风险。
本文要深入解析的,正是一种能够应对这一挑战的自适应虚拟谐波阻抗策略。我们将通过Simulink仿真,完整复现一个包含该策略的光储VSG并网模型。这篇文章不会停留在理论公式的罗列,而是聚焦于解决一个核心工程问题:如何在Simulink中构建一个能够实时感知电网阻抗变化,并自动调整虚拟谐波阻抗,从而实现对并网电流谐波含量进行闭环、精准调控的系统?
读完本文,你将获得:
- 对自适应虚拟谐波阻抗策略原理的清晰理解,明白它为何比固定策略更有效。
- 一套可复用的Simulink仿真模型搭建指南,从主电路到控制回路,步步拆解。
- 关键的实现细节与参数整定方法,了解哪些地方容易出错,以及如何调整。
- 完整的仿真结果分析与验证,直观看到策略如何抑制谐波、提升稳定性。
我们将从问题本质出发,逐步深入到建模与仿真实现,最终让你不仅能“跑通”这个仿真,更能“吃透”其背后的控制思想。
1. 核心问题:为什么固定虚拟阻抗在弱电网下会“失灵”?
要理解自适应策略的必要性,首先要看清固定策略的局限性。光储VSG通过逆变器接入电网,中间通常会加入LCL滤波器来滤除开关频率附近的高频谐波。然而,LCL滤波器本身在特定频率下存在谐振峰。为了阻尼这个谐振,传统方法是在控制中引入“虚拟阻抗”,通常是一个虚拟电阻,来模拟实际电阻的阻尼效果,从而稳定系统。
这一切在“强电网”(电网阻抗接近于0)的假设下工作良好。因为电网阻抗很小,整个系统的谐振频率主要由LCL滤波器自身的参数决定,是固定不变的。因此,针对这个固定谐振点设计的固定虚拟阻抗,效果显著。
但现实是,随着分布式电源高比例接入,电网末端的“弱电网”特性越来越明显。所谓“弱电网”,其核心特征就是电网等效阻抗较大且可能随时间、运行方式变化。这时,系统的总阻抗(滤波器阻抗+电网阻抗)发生了改变,导致实际的谐振频率点发生了偏移。
如果仍然使用针对强电网设计的固定虚拟阻抗,就会出现两个严重问题:
- 失谐:虚拟阻抗的阻尼效果在新的谐振频率点上大大减弱,无法有效抑制谐振,可能导致并网电流振荡甚至发散。
- 谐波放大:对于电网背景谐波(如5次、7次等),固定的阻抗特性可能在某些频段形成非预期的阻抗匹配,反而放大这些谐波,恶化电能质量。
因此,问题的核心从“如何设计阻尼”转变为“如何让阻尼策略跟上电网阻抗的变化”。自适应虚拟谐波阻抗策略,正是为了解决这个动态匹配问题而生的。它的目标是:实时检测或估计系统谐波频段的阻抗特性,并动态调整虚拟阻抗的数值或频率特性,始终将系统的谐波阻抗维持在一个理想的范围内,从而实现对谐波含量的闭环调控。
2. 自适应虚拟谐波阻抗策略原理剖析
自适应策略的核心思想可以概括为“感知-决策-执行”闭环。下面我们拆解其技术原理。
2.1 系统结构与问题定义
一个典型的光储VSG并网系统简化控制框图如下(重点关注谐波环):
[光伏/储能] -> [DC/AC逆变器] -> [LCL滤波器] -> [电网] ^ | ^ | | (并网点电压、电流反馈) [VSG控制层] | | | [功率环(基波)] [谐波抑制环] | | +----------+VSG控制层通常包含功率(有功-频率、无功-电压)下垂控制、虚拟惯量等,这部分模拟同步发电机的基波外特性。而我们关注的是谐波抑制环,它通常并联在基波控制指令之后,共同构成逆变器的调制信号。
设并网点的电压为v_g,电流为i_g。在谐波频率下,我们希望逆变器输出的谐波电压v_h与谐波电流i_h满足一个理想的“谐波阻抗”关系:Z_desired(s) = v_h(s) / i_h(s)。这个Z_desired通常设计为一个纯电阻或带一定频率特性的阻抗,用于提供阻尼并塑造谐波频段的输出阻抗。
2.2 自适应律的核心:模型参考自适应控制(MRAC)思想
自适应策略广泛采用模型参考自适应控制(MRAC)或其变种。其基本思路是:
- 可调模型:将包含待调整虚拟阻抗
Z_v(s, θ)的系统视为可调模型。其中θ是可调参数(如虚拟电阻值)。 - 参考模型:定义一个理想的参考模型,其传递函数代表了期望的系统谐波阻抗
Z_desired(s)。 - 误差驱动:通过比较可调模型输出(实际谐波电流
i_h)与参考模型输出(期望谐波电流i_h_ref)之间的误差e。 - 参数调整:利用这个误差
e,通过一个设计好的自适应律(如梯度法、Lyapunov稳定性理论推导的律)来在线调整参数θ,使得误差e趋近于零。
通俗地说,就是系统不断问自己:“我现在的谐波阻抗表现和期望的‘模范生’(参考模型)差多少?”然后根据这个差距,自动修改自己的“虚拟阻抗”这个内部参数,努力向“模范生”看齐。
2.3 谐波提取与阻抗辨识
实现自适应的前提是能准确获取谐波分量。通常采用:
- 多谐振控制器:如针对特定次谐波(5th, 7th, 11th, 13th)的PR控制器或准PR控制器,在指定频率处提供高增益,从而分离出该次谐波指令或反馈。
- 带通滤波器组:设置中心频率为各次谐波的带通滤波器,提取相应频段的信号。
- 基于FFT/DFT的在线辨识:在控制周期内计算特定谐波分量的幅值与相位,但计算量较大。
对于电网阻抗的感知,有时不需要直接测量阻抗值,而是通过监测谐波环路的稳定性指标(如特定频段的相位裕度)或谐波电流的放大情况来间接驱动自适应过程。
3. Simulink仿真环境搭建与关键模块
现在,我们进入实战环节,在Simulink中搭建仿真模型。假设你使用的是MATLAB/Simulink R2021a或更新版本。
3.1 主电路与滤波器参数设计
首先建立主功率电路。
- 直流源:使用
DC Voltage Source模块模拟光伏阵列或储能电池的直流母线,电压设为800V。 - 三相逆变桥:使用
Universal Bridge模块,选择IGBT/Diodes,载波频率设置为10kHz。 - LCL滤波器:这是关键物理部件。参数设计需要权衡滤波效果、体积成本和系统稳定性。一组参考参数如下:
- 逆变器侧电感
L1: 1.2 mH - 网侧电感
L2: 0.3 mH - 滤波电容
Cf: 30 μF(通常串联阻尼电阻,如0.5 Ohm) - 计算谐振频率:
f_res = 1/(2π) * sqrt((L1+L2)/(L1*L2*Cf))。用上述参数计算约 1.8kHz,应远离基波(50Hz)和开关频率(10kHz),处于合理范围。
- 逆变器侧电感
- 电网模型:使用
Three-Phase Programmable Voltage Source模块模拟电网。可以设置基波电压(如311V峰值),并添加谐波电压源来模拟背景谐波,例如加入3%的5次谐波和2%的7次谐波。为了模拟弱电网,需要串联一个Three-Phase Series RLC Branch作为电网阻抗,设置其电阻Rg和电感Lg可调(例如Rg=0.1 Ohm,Lg=2 mH)。
3.2 基础VSG控制层实现
在控制部分,先搭建标准的VSG算法。
- 锁相环(PLL):使用
PLL (3ph)模块或SRF-PLL子系统,获取电网电压相位θ和频率ω_g。 - 功率计算:使用
abc to dq变换模块,将并网点三相电压v_g_abc和电流i_g_abc变换到同步旋转坐标系(dq轴)。有功功率P和无功功率Q通过P = 1.5*(v_d*i_d + v_q*i_q),Q = 1.5*(v_q*i_d - v_d*i_q)计算(通常v_q设为0)。 - VSG机械方程:模拟转子运动方程,生成虚拟频率
ω_vsg和相位θ_vsg。
其中% 伪代码体现在Simulink中用基本运算模块搭建 % J * d(ω_vsg)/dt = P_set - P - D*(ω_vsg - ω_g) // 有功-频率下垂 % d(θ_vsg)/dt = ω_vsg % K * (V_set - V_g) = Q_set - Q // 无功-电压下垂J为虚拟惯量,D为阻尼系数,K为电压调差系数。P_set,Q_set,V_set为设定值。 - 电压电流双闭环控制:
- 外环(电压环):以VSG输出的电压幅值
V_vsg和相位θ_vsg作为参考,与电容电压反馈比较,通过PI控制器生成电流参考。 - 内环(电流环):以电压环输出的电流为参考,与逆变器侧电感电流反馈比较,通过PI控制器生成调制电压。电流环带宽通常设计在1kHz左右以保证动态性能。
- 外环(电压环):以VSG输出的电压幅值
3.3 自适应虚拟谐波阻抗模块实现
这是本文的核心。我们以针对主要低次谐波(如5次、7次)的自适应电阻为例。
- 谐波提取:对并网电流
i_g_abc进行谐波分离。- 方法一:使用多个
Parallel RLC Branch作为带通滤波器,中心频率分别设为250Hz和350Hz。 - 方法二(更精确):在dq坐标系下,使用基于二阶广义积分器(SOGI)的谐波提取器,或直接使用
Discrete 3-Phase PLL并配置其谐波提取功能。 - 我们以方法一为例,在Simulink库中找到
Analog Filter Design或Discrete Filter模块,设计两个带通滤波器。
- 方法一:使用多个
- 构建自适应律:
- 定义可调参数:设虚拟谐波电阻为
R_vh(初始值可设为1 Ohm)。 - 定义参考模型:期望的谐波阻抗就是一个纯电阻
R_desired(例如2 Ohm)。这意味着希望谐波电压与谐波电流同相位,且比例为R_desired。 - 计算误差:对于提取出的5次谐波电流
i_h5,期望的谐波电压降为v_h_ref5 = R_desired * i_h5。实际中,我们通过测量并网点谐波电压v_h5和电流i_h5,可以计算当前等效阻抗R_actual = v_h5 / i_h5(需处理幅值)。误差可以定义为e = R_desired - R_actual或基于功率的误差。 - 实现参数更新:采用积分自适应律,
d(R_vh)/dt = γ * e * i_h5^2(或类似形式,γ为自适应增益)。在Simulink中用Integrator模块实现。
% 在Simulink中用基本模块搭建的自适应律示意 % 输入:i_h5 (5次谐波电流), v_h5 (5次谐波电压) % 输出:R_vh (自适应虚拟电阻) % 步骤: % 1. 计算当前阻抗幅值(简化):R_actual = abs(v_h5) / abs(i_h5) (需经过低通滤波平滑) % 2. 计算误差:e = R_desired - R_actual % 3. 积分更新:R_vh = (γ * e) 的积分 + R_vh_initial % 注意:为防止积分饱和,需设置上下限,如 [0.1, 10] Ohm。 - 定义可调参数:设虚拟谐波电阻为
- 注入虚拟阻抗:将自适应得到的
R_vh乘以提取的谐波电流i_h,得到一个谐波电压补偿量v_h_comp = R_vh * i_h。将这个补偿量叠加到电流环的参考电压上,或者作为额外的调制信号注入。这样,控制器就等效在谐波频段增加了一个值为R_vh的虚拟电阻。
4. 完整仿真模型集成与参数设置
将上述所有模块集成到一个完整的Simulink模型中。
- 子系统封装:将VSG控制、谐波提取、自适应律分别封装成子系统,使模型结构清晰。
- 关键参数设置:
- 仿真器:选择
ode23tb或ode15s等适用于电力电子仿真的变步长求解器,相对容差设为1e-4。 - 采样时间:控制部分使用固定的离散采样时间,例如
1e-5秒(100kHz),与实际的数字控制器一致。 - 自适应增益 γ:这是一个关键且敏感的参数。太小则自适应速度慢,太大则可能引起振荡。通常从很小值(如0.01)开始尝试,观察
R_vh的收敛过程。 - 期望电阻
R_desired:根据系统稳定性和谐波衰减要求选择,通常在1-5 Ohm之间。 - 电网阻抗
Lg:设置一个变化场景。例如,在0.4秒时将电网电感从1 mH突变到4 mH,以模拟电网强度变弱。
- 仿真器:选择
5. 仿真运行与结果分析
运行仿真,重点观察以下信号,并与采用固定虚拟阻抗的对照组进行对比。
场景一:电网阻抗阶跃变化
- 观察
R_vh的变化曲线:自适应电阻应能跟随电网阻抗的变化而自动调整,最终稳定在新的平衡点附近。 - 观察并网电流总谐波畸变率(THD):使用
Powergui模块的FFT Analysis工具。在电网阻抗突变后,自适应策略应能比固定策略更快地将THD恢复到较低水平(如<5%)。 - 观察特定次谐波(如5次)电流幅值:自适应策略应能更好地抑制该次谐波。
- 观察
场景二:背景谐波电压扰动
- 在电网电压源中设置持续的5次、7次谐波。
- 对比分析:固定阻抗下,谐波电流可能被放大;而自适应策略下,由于虚拟电阻调整到合适值,谐波电流应得到有效抑制。
关键波形对比(建议用Scope分组显示)
- Group 1: 自适应参数
R_vh(自适应虚拟电阻)R_actual(估算的实际谐波阻抗)e(误差)
- Group 2: 电能质量
i_g_a(A相并网电流)v_g_a(A相电网电压)THD(实时或周期计算的THD)
- Group 3: 谐波分量
i_h5(5次谐波电流)i_h7(7次谐波电流)
- Group 1: 自适应参数
预期效果:在自适应策略下,当电网阻抗增大导致系统谐波阻抗特性恶化时,R_vh会自动增大,提供更强的阻尼,从而稳定谐波环,抑制谐振峰,降低谐波电流。你会看到并网电流波形更正弦,THD值更低且更稳定。
6. 常见问题与排查思路
在复现过程中,你可能会遇到以下问题:
| 问题现象 | 可能原因 | 排查方式 | 解决方案 |
|---|---|---|---|
| 仿真发散或不收敛 | 1. 自适应增益γ过大。2. 电流环或电压环PI参数不合理。 3. 求解器设置不当。 | 1. 检查R_vh是否剧烈振荡。2. 先断开自适应环,只运行基础VSG,确保稳定。 3. 尝试更小的仿真步长或更严格的求解器。 | 1. 大幅减小γ,如设为0.001或更小。2. 重新整定双闭环PI参数,确保基波系统稳定。 3. 使用 ode23tb求解器,减小最大步长。 |
自适应效果不明显,R_vh几乎不变 | 1. 自适应增益γ过小。2. 谐波提取模块失效, i_h信号接近为零。3. 误差计算逻辑有误。 | 1. 检查i_h信号是否有足够的幅值。2. 检查带通滤波器的中心频率和带宽设置是否正确。 3. 检查 R_actual的计算是否合理(需取幅值、滤波)。 | 1. 逐步增大γ。2. 校准谐波提取模块,注入已知谐波测试。 3. 简化误差,尝试直接用 e = R_desired * i_h - v_h。 |
| 并网电流THD反而变差 | 1. 虚拟阻抗注入点错误,引起相位问题。 2. 自适应律引入正反馈。 | 1. 检查v_h_comp的注入是负反馈(与谐波电流方向相反以抑制它)。2. 在频率域分析系统环路,检查谐波频段相位。 | 1. 确保补偿电压为-R_vh * i_h(负号)。2. 在自适应更新律中引入泄漏项(leakage),防止积分漂移。 |
| 特定次谐波放大 | 1. 虚拟阻抗在该次谐波频率处呈现容性或感性,与系统阻抗发生谐振。 2. 滤波器参数不匹配。 | 1. 分析系统在谐波频率处的开环伯德图。 2. 检查LCL滤波器谐振频率是否偏移。 | 1. 考虑采用带频率特性的虚拟阻抗(如R-L或R-C串联),而不仅是纯电阻。 2. 重新核算LCL参数,确保谐振点远离主要谐波频率。 |
7. 工程实践建议与扩展方向
成功复现仿真只是第一步,要将此策略应用于实际工程或深入研究,还需考虑以下方面:
- 离散化实现:Simulink中的连续模型需转化为离散算法才能在DSP/FPGA中运行。注意对自适应律中的积分器进行离散化(如采用前向欧拉法),并考虑计算延时。
- 稳定性与鲁棒性分析:自适应系统可能存在非线性和时变特性,需采用李雅普诺夫稳定性理论或超稳定性理论严格证明所设计自适应律的稳定性,并分析其对参数变化、测量噪声的鲁棒性。
- 多谐波协同治理:本文示例针对单一谐波(如5次)。实际中需同时治理多个谐波。可以扩展为多个并行的自适应环,每个环针对一个主要谐波频段,但要注意各环之间的耦合影响。
- 与有源滤波功能的结合:光储VSG不仅可以抑制自身产生的谐波,还可以通过控制补偿电网中的背景谐波。自适应虚拟阻抗策略可以扩展为谐波电流补偿策略,使VSG具备有源滤波功能。
- 参数自整定:自适应增益
γ和期望阻抗R_desired本身也可以根据系统运行状态进行优化,实现二级自适应。
通过本文的详解与Simulink仿真实践,你应当已经掌握了自适应虚拟谐波阻抗策略从原理到实现的关键路径。这项技术的价值在于,它使光储并网系统从“被动耐受”电网环境,转向“主动适应”电网变化,是提升高比例新能源接入背景下系统韧性的重要技术手段。建议你在理解本模型的基础上,尝试修改电网阻抗变化规律、增加更多次谐波、或尝试不同的自适应律,以深化理解并探索其性能边界。模型中的具体参数需要根据你的实际系统规格进行调整,但核心的构建思路和控制逻辑是相通的。