光伏-储能并网系统仿真(MATLAB/SIMULINK)-part1

光伏-储能并网系统仿真(MATLAB/SIMULINK)-part1

如果不想看原理，可以直接点击文件下载连接

本文介绍如何利用MATLAB/SIMULINK完成一个光伏-储能并网系统的搭建。
一篇文章就能读懂光伏并网系统的控制，以及PI控制器的设计方法。
本文将会介绍四个部分：

1. 整体系统控制逻辑
2. 光伏并网逆变器的控制理论（dq解耦控制）
3. 光伏MPPT控制以及恒功率控制
4. 电池的功率控制

每个部分都会包括：

1. 控制器原理
2. PI控制器参数理论计算

1. 整体系统控制逻辑

1.1 系统要求

图1 系统整体结构 图1展示了系统的基本结构包括各个模块的功率和连接关系。

1. 储能DC/DC变换器与光伏DC/DC变换器共直流母线（母线电压750V）运行，母线电压纹波需要小于5%。
2. 储能电池系统双向功率控制，达到平滑系统功率输出的目的

1.2 系统控制方法

该系统需要考虑了两种场景分别是：

1. 场景一：电池SOC较为健康（40%-90%），此时可以帮助平抑光伏的功率波动
2. 场景二：电池SOC过低或者过高，此时电池需要恒功率的充电或者放电，以保护电池

针对这两种场景，系统需要采取不同的控制方式。
对于场景一，按照表1所示的控制方式可以实现系统有功输出和无功输出都保持恒定，从而可以满足系统的要求。
对于场景二，按照表1所示的控制方式，进行光伏弃光（即恒功率控制），同时对电池进行恒功率充电或者放电以保持系统总功率输出不变。

场景	光伏模块	电池	逆变器
场景一（电池SOC健康）	MPPT	输出系统有功输出和光伏功率输出的差值	控制母线电压+控制无功输出
场景二	恒功率控制	SOC过低（高）：恒功率充（放）电	控制母线电压+控制无功输出

表1 两种模式下的控制方法 注：MPPT=最大功率点跟踪，具体原理会在后续阐释

1.3 系统控制原理解释

控制的基本原理是：维持母线电容电压的稳定等价于光伏输出功率加上电池输出功率（充电记为负功率）与逆变器输出功率平衡。

如果记光伏出力为$P_1$，电池出力为$P_2$，逆变器出力为$P_3$。母线电压稳定就意味着电容的能量保持稳定，因此有如下关系：
$$P_1+P_2=P_3$$
为了更清晰地解释，我们假设系统运行在场景一下。此时，光伏采用MPPT控制，因此光照和温度的变化都会引起光伏出力（$P_1$）的变化。假设系统要求出力为$P^{\ast }$，那么应该控制电池输出功率为$P_2$（$P_2=P^{\ast }-P_1$）。又因为逆变器控制了母线电压稳定，因此整个系统的输出就是$P_3$（$P_3=P_1+P_2=P^{\ast }$）。因此也就满足了系统处理要求。

2. 光伏并网逆变器的控制理论（dq解耦控制）

2.1 逆变器有功无功控制

逆变器与电网的连接可以用如下公式近似表示：
$$U_{abc}=i_{abc}R+L\frac{d{i}_{abc}}{dt}+e_{abc}（1）$$
其中，下表abc表示abc三相所对应量所构成的列向量；$U_{abc}$表示逆变器的电压；$i_{abc}$表示三相电流；$e_{abc}$表示网侧的三相电压。
对公式(1)进行park变换，可以得到如下的公式：
$$\{\begin{array}{c}{v}_d=R{i}_d+L\frac{d{i}_d}{dt}-\omega L{i}_q+e_d\\ v_q=R{i}_q+L\frac{d{i}_q}{dt}+\omega L{i}_d+e_q\end{array}(2)$$
其中， $\omega$为park变化中dq坐标系的旋转角速度（一般选取为电网的频率）。（park变换的相关知识见文章park变换是什么？）
经过变换后，输出给电网的有功功率和无功功率可以表示为如下公式：
$$\{\begin{array}{c}P=\frac{3}{2}(e_d{i}_d+e_q{i}_q)\\ Q=\frac{3}{2}(e_q{i}_d-e_d{i}_q)\end{array}(3)$$
其中，P，Q分别代表输出给电网的有功和无功功率。
更进一步，在进行park变换时可以定义电网侧的电压矢量与d轴平行。（所谓的电压矢量与d轴平行其实就是在进行park变换时，通过调整变换矩阵的相位使得$e_q=0$）。因此，我们可以看到当$e_q=0$时，上面的功率的公式变成了更加简单的形式：
$$\{\begin{array}{c}P=\frac{3}{2}(e_d{i}_d)\\ Q=\frac{3}{2}(-e_d{i}_q)\end{array}(4)$$
上面的公式也就是常说的q轴电流控制无功，d轴电流控制有功。（但其实如果令电压矢量与q轴平行，不就是q轴电流控制有功了？ 😃 ）

那么该如何控制dq轴电流呢？下面介绍电流环的设计思路和流程。

假设id和iq的参考值（$i_d^{\ast },i_q^{\ast }$）已知，那么我们可以通过控制逆变器的电压输出使得$i_d=i_d^{\ast },i_q=i_q^{\ast }$。
首先，对公式（2）进行laplace变换：
$$\{\begin{array}{c}{i}_d=\frac{V_d+\omega L{i}_q-e_d}{R+Ls}\\ i_q=\frac{V_q-\omega L{i}_d-e_q}{R+Ls}\end{array}(5)$$其中，s为Laplace算子。
$V_d$和$V_q$是逆变器的dq电压，因此这两个量是我们可以控制的。如果我们采取的是PI控制器，那么我们可以令这两个控制量为如下形式：
$$\{\begin{array}{c}{V}_d=e_d-\omega L{i}_q+(K_p+K_i/s)(i_d^{\ast }-i_d)\\ V_q=e_q+\omega L{i}_d+(K_p+K_i/s)(i_q^{\ast }-i_q)\end{array}(6)$$
将（6）带入（5），我们就可以整理得到如下的等式：
$$\{\begin{array}{c}{i}_d=\frac{\frac{K_p}{L}s+\frac{K_i}{L}}{s^2+(\frac{R+K_p}{L})s+\frac{K_i}{L}}{i}_d^{\ast }\\ i_q=\frac{\frac{K_p}{L}s+\frac{K_i}{L}}{s^2+(\frac{R+K_p}{L})s+\frac{K_i}{L}}{i}_q^{\ast }\end{array}(7)$$
式（7）表明如果我们确实可以让控制变量满足式（6）所示的关系，那么整个系统就可以表示为式（7）所示的一个二阶系统。并且我们可以看到如果$i_d^{\ast }$是一个阶跃信号，那么利用Laplace的终值定理可以得到如下关系：
$$\underset{t\to \infty }{\lim }{i}_d(t)=\underset{s\to 0}{\lim }\frac{\frac{K_p}{L}s+\frac{K_i}{L}}{s^2+(\frac{R+K_p}{L})s+\frac{K_i}{L}}\frac{i_d^{\ast }}{s}\ast s=i_d^{\ast }(8)$$
这也就意味着随着时间的推移，$i_d$就会趋向于我们设定的$i_d^{\ast }$。
可是，我们往往不希望需要等到海枯石烂才能趋向于我们设定的值，因此我们还需要对PI控制器参数进行合理的配置，才能让式（7）所示的二阶系统具有足够好的动态要求。（合理的超调量和满足要求的收敛时间）
配置的方式有很多，这里介绍一种。
我们可以写出二阶系统所对应的阻尼系数和自由振荡频率，并且我们可以规定我们所期望的系统收敛时间$t_r$：
$$\{\begin{array}{c}\xi =\frac{R+K_p}{2\sqrt{K_{i}L}}=0.8\\ {\omega }_n=\sqrt{\frac{K_i}{L}}=\frac{3}{\xi t_r}\\ t_r=0.01\end{array}(9)$$上面我们设定了系统的收敛时间为$t_r=0.01s$，阻尼系数为0.8。两个方程以及两个未知数（$K_p,K_i$），我们利用式（9）不难求得合适的PI参数值。

图2 电流环控制框图 图2展示了上面的理论推导所对应的控制框图。

以上内容解决了电流环的设计问题，那么剩下的问题便是：我们如何给出合适的dq电流参考值（$i_d^{\ast },i_q^{\ast }$）从而控制有功功率和无功功率呢？
下面介绍母线电压控制和无功控制环的思路

根据式（4），我们知道，如果$i_d$影响着逆变器的有功输出大小，因此如果母线电压较低时，我们就需要减小$i_d$，这样逆变器的功率输出（$P_3$）就会减小，导致光伏和电池输出的能量之和大于逆变器的功率输出（$P_3$），进而可以对电容进行充电也就提高了电容的电压。相反，母线电压低于参考值时，我们就应该提高$i_d$。
母线电压$U$的变化可以用如下公式近似表示：
$$\frac{d}{dt}(\frac{1}{2}C{U}^2)=CU\frac{dU}{dt}=\frac{3}{2}{e}_d{i}_d-P^{\ast }(10)$$其中，$P^{\ast }$表示系统所要求的有功处理。显然，这里做了一个假设：光伏的有功输出和电池的有功输出之和已经接近系统所要求的有功出力要求。我们进一步假设$U$已经接近设定的母线电压值$U_{bus}$。式（10）可以近似线性化为：
$$C{U}_{bus}\frac{dU}{dt}=\frac{3}{2}{e}_d{i}_d-P^{\ast }(11)$$其中，$C$为母线电容容值。
对线性化的式（11）进行Laplace变换，我们可以得到：
$$sU=\frac{3}{2CU}{e}_d(i_d-\frac{2}{3e_d}{P}_{ref})(12)$$
此时，$i_d$是我们的控制变量（请注意，借助电流环，我们可以控制$i_d$为我们的设定值，因此可以将它看作控制变量，也就是电流环是内环）。令$i_d$为如下形式：
$$i_d=(K_p+K_i/s)(U-U_{ref})+\frac{2}{3e_d}{P}_{ref}(13)$$
将（13）代入（12），我们可以得到如下系统：
$$U=\frac{C_0{K}_{p}s+C_0{K}_i}{s^2+C_0{K}_{p}s+C_0{K}_i}{U}_{ref}(14)$$其中，$C_0=\frac{3e_d}{2U_{bus}C}=\frac{3\times 220\sqrt{2}}{2\times 750\times C}$
类似于电流环的控制分析，我们可以知道式（14）代表，母线电压最终会趋向于$U_{ref}$。PI参数的配置方法和电流环的配置是类似的：
$$\{\begin{array}{c}\xi =\frac{1}{2}\sqrt{\frac{C_0}{K_i}}{K}_p=0.8\\ {\omega }_n=\sqrt{C_0{K}_i}=\frac{3}{\xi t_r}\\ t_r=0.05\end{array}(15)$$
这里我们希望外环的收敛速度更慢。这是因为我们假设了对于外环来说，内环的控制速度是最够快的。
无功控制环，的配置就不再赘述了。这里只给出其控制框图：

图3 无功控制环

逆变器控制的最后内容！——锁相环
显然，我们上面的分析都用到了dq值。那么我们该怎么得到dq值呢？现实中，我们只能测量得到三相的电流或者电压，因此我们需要一个器件得到dq值。这个器件就是锁相环。
我们希望park变换后，电网侧的电压的q轴电压为零，因此，我们输入的就是电网侧的三相电压。锁相环的具体原理见下图。

图4 锁相环原理 显然，原理图实现了将网侧电压的q轴控制为零。MATLAB中有现成的锁相环模块，因此可以直接使用。

剩下的部分见part2。