光伏逆变器boost电路阻抗无源化重塑研究
光伏逆变器boost电路阻抗无源化重塑研究
edit:LPEI
gitee仓库 [ matlab 仿真文件]
目录
1、前言:
当光伏逆变器运行在光伏阵列 PV 曲线左半区域时,光伏阵列呈现出高内阻特性,
易诱发光伏逆变器前级 boost 电路振荡。
2、建模:
光伏逆变器前级boost电路可以等效为如下电路图:

其平均状态空间方程为:
传递函数为:
当光伏阵列等效内阻 Zpv 非常大时,上述系统中的阻尼项几乎消失,系统变为容易失稳
的 LC 谐振系统,其波特图如下:

一般可通过状态变量反馈实现有源阻尼,等效为电容并联虚拟电阻,如下图所示。
然而,引用虚拟电阻后,系统仍可能不稳定,依旧需要通过环路分析系统的稳定性。
参考LCL型逆变器阻抗无源化重塑方法,可以引入boost电路的等效输出阻抗,
利用阻抗稳定性分析LC谐振。
3、稳定性分析:
将光伏电池的输出电流视作外部输入,控制框图如下图所示:

合并Upv反馈支路:
其中Zo为红框内的闭环系统,其传递函数为:
定义:
根据上述,boost 电路可等效为理想受控电流源 Gg(s)Uref(s) 和恒流源 is(s) ,
并联输出阻抗 Zo(s) ,结合图1可得到如图6所示的等效电路。
根据上图,可以推出光伏电池输出电压 Upv 的表达式:
由上式可知,光伏电池输出电压的稳定性由N1(s)和N2(s)决定。
其中N1(s)等价于Zo(s)为无穷大时,光伏电池输出电压稳定性问题。
对应下面电路图:
对应控制框图如下:
这是一个开环系统,只需满足控制器自身稳定,系统就稳定。因此,可以认为 N1(s) 是
稳定的,光伏电池输出电压 Up v稳定性问题等价于 N2(s )的稳定性问题,N2(s) 等效为前
向通道为 1 ,反馈通道为 Zpv(s)/Zo(s) 的负反馈控制系统的闭环传递函数。为了保证 N2
(s) 稳定,Zpv(s)/Zo(s) 需满足 Nyquist 稳定性判据。根据上述阐述,可得到如下基
于阻抗分析法光伏电池输出电压稳定条件:
- (1) 控制器自身稳定;
- (2) 阻抗比 Zpv(s)/Zo(s)满足 Nyquist 稳定性判据。
这里对比一下跟网型逆变器阻抗分析方法的区别:

注意:boost电路等效的输出阻抗 Zo(s) 是非最小相位系统,通过 Nyquist 判断稳定性需要先确定 Zo 的不稳定极点数。
4、临界稳定:
从上述分析可知,光伏电池输出电压稳定问题等价于 N2(s) 的稳定性问题,
N2(s) 的开环传递函数为 Zpv(s)/Zo(s),考虑 Zpv(s) 为电阻 Rpv 情况,
那么 Rpv 可视为开环传递函数 1/Zo(s) 的增益,利用根轨迹,
可以知道内阻 Rpv 对系统稳定性影响。忽略PWM发波延迟,
以理想传递函数分析。参数取:
% matlab 脚本
fs = 20000;
Ts = 1/fs;
gain = 10;
GpiP = 1e-5 * gain;
GpiI = 0.05 * gain;
L = 0.6e-3;
C =20e-6;
Udc =700;
% case1: stable
% E = 4000;
% Rpv = 150;
% case2: unstable
E = 4000;
Rpv = 155;
% LC 谐振频率
wr = sqrt(1/(L*C));
fr = wr/2/pi
下图是根轨迹图,当 Rpv>152 时,极点出现在右半平面,系统不稳定。
% 分析脚本
s = tf([1 0],1);
Gpi = GpiP + GpiI/s;
Zo = s*L/(1 + Gpi*Udc + L*C*s*s);
%% 方法1:根轨迹,Zpv 可以视为开环系统的增益,因此可以分析 Zpv 临界稳定值
rlocus(1/Zo)
%% 方法2:直接判断 Zpv = Rpv
Rpv = 160;
N2 = Rpv/Zo;
margin(N2)
%% 方法3:求出闭环系统,分析极点
Gcl= feedback(N2,1);
step(Gcl)
pole(Gcl) % 极点在右半平面不稳定
simulink仿真,Rpv=150 时渐进收敛,Rpv=155 振荡。
5、阻抗无源化重塑:
若输出电阻 Zo(s) 在全频段都满足 Real(Zo(jw))>=0 ,则称输出阻抗 Zo(s) 是无源
的,基于无源性的稳定条件是光伏电池输出电压稳定的充分条件,据此对输出阻抗
进行无源化重塑,可保证宽范围的等效内阻下系统的稳定性。
从上式可以看出,输出阻抗 Zo(s) 是电感、电容和附加阻抗 Zd(s) 的并联,定义:
则等效电路如图所示:

若控制器中不含有积分环节,仅包含比例环节,在忽略延迟作用下,Zd(s)则表现为纯电感。
理论上,可以在任意 Zpv 下保持稳定。
Zo(s)阻抗重塑可以从两个方面入手:
- 1、修改控制器结构直接改变 Zd(s);
- 2、在Zo(s)中并联或串联 Gx。
从图5可知,对 Zo(s) 串并联阻抗等效控制框图的变化如下所示。
Zo 串联 Gx 实现:
Zo 并联 1/Gx 实现:
由于逆变器中通常没有 ipv 电流采样传感器,所以只能实现实现 Zo 并联 1/Gx 方式
的阻抗重塑。考虑一般情况,控制器采用 PI 控制器,由于时并联阻抗重塑,
将 Zd(s) 写为导纳形式 Yd(s)=1/Zd(s),具体表达式为:
不考虑延迟,以及 PWM 发波近似为1,Yd(s) 中的负电阻仅和控制器积分项有关,
代入 s=jw 可知,积分项完全表现为负电阻,因此进行阻抗重塑需要抵消该部分。
并联阻抗 1/Gx(s) ,如下图所示。

从表达式可以看出,这样实现前馈等价于在控制器中直接删去积分环节。
本质上,Upv 电压前馈和控制器修改是等效的,因为两者的引入点都是 Upv电压,
在控制框图中可以合并到同一个支路上。因此,必须采样其他状态量进行前馈,
才可以区别于直接修改控制器进行阻抗重塑。
由于 Upv 前馈方案与修改控制器方案近似,下面就直接修改控制器进行阻抗重塑。
方案1:用惯性环节替代积分环节,在 PI 控制中并联负的惯性环节。
阻抗重塑后,高频部分 Zd(s) 相位在 90 度以内,为正电阻,右边是根轨迹图,
可以看出极点都在左半平面,不管Rpv如何变化,系统都稳定。
% matlab 阻抗重塑脚本
s = tf([1 0],1);
%% 原始阻抗
Gcc = GpiP + GpiI/s;
Zo = s*L/(1 + Gcc*Udc + L*C*s*s);
Zd = s*L/(Gcc*Udc);
subplot(1,2,1)
bode(Zd); % 正电阻,相角应在 -90-90度之间
subplot(1,2,2)
rlocus(1/Zo)
%% case1: 负惯性环节
Glow = 1.2e-4/(s/10000 + 1);
Gcc1 = GpiP + GpiI/s - Glow;
Zd1 = s*L/Udc/Gcc1;
Zo1 = s*L/(1 + Gcc1*Udc + L*C*s*s);
subplot(1,2,1)
bode(Zd1); % 正电阻,相角应在 -90-90度之间
subplot(1,2,2)
rlocus(1/Zo1)
参考:
《LCL 型并网逆变器谐振阻尼和阻抗无源化重塑研究》
后续:
点赞超1000,更新逆变器阻抗建模。阻抗建模,可以分析强网下LCL谐振和弱网振荡问题。
参考:《新能源发电接入弱电网的宽频带振荡机理及抑制方法研究_伍文华》
改进:使用控制框图推导,将原本的方程组转换到控制框图中分析。
步骤:
1、恒压源下逆变器dq坐标系建模
2、转换到abc坐标系下,只考虑a相的正序和负序控制框图
3、锁相环建模,合并到上述控制框图中
4、在上述控制框图中,对电压外环进行小信号建模,等效对电感和Udc修正,合并到上述控制框图中
5、直流侧光伏以boost建模,等效到直流母线上,对Udc表达式进行修正,最终合并到上述控制框图中
AtomGit 是由开放原子开源基金会联合 CSDN 等生态伙伴共同推出的新一代开源与人工智能协作平台。平台坚持“开放、中立、公益”的理念,把代码托管、模型共享、数据集托管、智能体开发体验和算力服务整合在一起,为开发者提供从开发、训练到部署的一站式体验。
更多推荐



所有评论(0)