VSG多机并联系统小信号建模

参考文献：曾德银等：虚拟同步发电机多机并联系统的频率小信号稳定性分析研究

多个VSG并联系统拓扑

看单条线路上的情况，从左到右：

VSG端口电压：edqi（i=1、2、3）是第i台VSG输出的端口电压

滤波器参数：Rfi是滤波器等效内阻，Lfi是滤波器电感，Cfi是滤波器电容

考虑漏感：Rei是VSG到PCC处线路电阻，Lei是漏感等效电阻

电网端：Rg是电网等效内阻，Lg是电网等效电感

uodqi和iodqi是第i太VSG输出电压和电流

ubdq是PCC点电压

在仿真和建模的时候，忽略了和频率响应不在一个时间尺度的一些变化和因素，就是说这部分发生的太快或太慢，和频率响应影响因素不一致，比如开关延时、电压电流内环控制和功率计算的低通滤波。

我们主要研究在机电时间尺度下通过分析系统外特性来研究VSG多机并联的频率稳定性。

参考坐标系和坐标变换

先把所有VSG变换到同一个坐标系下：这里选的是第一台VSG的dq坐标系（不懂dq坐标系的看电力系统分析，同步发电机那节里面讲的有，派克变换）作为参考，那么第i台和第一台的角度差为

坐标系中表现为：

 

解释：攻角上面打点就是表示它的微分，这个其实就是上面那个式子两边取微分得到的，表现为变化量，就是它的小信号建模，这个在模电那本书里面有，可以去看看，另外，这是一种常见的小信号表示，上面打点表示微分。

在建模的时候，要把第i个VSG控制的逆变器输出电流iodqi转变到DQ坐标系下，方法如下：

 线性化以后的小信号模型如下：

 解释：这里是什么意思呢？为什么dq坐标系下的i还要变换到另一个DQ坐标系呢？其实DQ是第一台VSG的旋转坐标系，而我们前面分析的和这个iodqi是第i台VSG的旋转坐标系，这两个坐标系的差值是 δ i，也就是我们把第一台VSG角度作为0，参考角度，那么第i台就是 δ i。这里这个公式是将原本在第i台VSG的dq坐标系下的输出电流变到第一台DQ坐标系下面去了，而Tsi和Tci具体，这是一个复杂一点的矩阵。

同理：

把位于第i台VSG处的uob变换到第一台VSG坐标系下的小信号建模是这样的。

功率控制环小信号建模

单台VSG功率外环控制的数学表达式是这样的：

 对其线性化，得到小信号模型:

然后通过输出电压和输出电流求瞬时功率（dq轴下功率计算）：

对其进行小信号线性化得到功率反馈小信号模型：

 结合上面的式子，把瞬时功率带入功率外环小信号模型里面去，得到VSG功率外环小信号模型标准形式：

 ABCD矩阵参数：

滤波器和连接线路小信号建模

如果假定我们功率外环输出的参考电压就是逆变器的输出端电压（指令等于实际值），计算在LC滤波阶段和线路上的数学模型:

 对其线性化，得到小信号模型：

 AB参数矩阵:

VSG单机和多机并联系统小信号建模

单机

以第1台VSG作为公共参考坐标系，结合第i台VSG和第一台VSG在坐标系上的换算方法，单台VSG功率外环计算方法，单台VSG输出滤波电路小信号建模，可以得到多台并联的情况下，第i台VSG在公共DQ坐标系下的小信号模型：

 单台VSG系统状态变量如下：

 参数矩阵：

多机

上面的是单台的情况，对于多台，可以直接写为下面的形式：

 其中参数矩阵：

输电线路和无穷大电网的小信号模型

既然是视为无穷大电网，容量是无穷大的，那么电压幅值和频率都是不变的：设为Eg和ωg，其与第一台参考VSG的相位差δg表示为：

 结合拓扑，可以得到PCC点和无穷大电网的电路方程如下：其中的igDQ表示公共DQ坐标系下流经输电线路的总电流（总汇入电网电流）：

 结合前面的式子，可以得到无穷大电网的小信号建模：

VSG多机并联并网系统小信号建模

相当于是一个汇总，把之前的VSG并网（包括滤波电路）、输电线路、无穷大电网的小信号建模结合一起，得到VSG的并联并网系统小信号建模，为了消除并网点的节点电压δubDQ（就是不想那这个表达式带进去算，为什么呢？可能是为了消除电网电压处的攻角计算吧），引入了对第电阻RVN，如果RVN足够大，那么可以忽略Rvn对系统的影响，同时可以将上面几个部分模型统一到一起来，引入虚拟对地电阻后PCC点处的节点电压方程如下：

 做线性化，得到小信号模型：

 上面的所有小信号模型最后的结果，得到总的小信号模型，所谓的VSG多机并联并网系统的通用小信号模型：

 其中的[Δxsys]为

一些参数：