电气传动2021年第51卷第2期
Stability Analysis and Parameter Designing for Virtual Synchronous Generator
CHENG Shiyao ，WANG Fang ，LI Xuecheng ，SHI Pengbo ，CHI Yuan
（Electric Power Research Institute ，State Grid Beijing Electric Power Company ，Beijing 100075，China ）Abstract:The virtual synchronous generator （VSG ）control scheme formed by the cascaded outer power control loop in voltage and current double closed loop core can introduce inertia into the inverter ，so as to enhance the stability of the system after grid connection.Aiming at the parameter design and stability of this VSG ，a VSG parameter tuning method guided by the eigenvalue sensitivity matrix of the linearized system model was proposed.Due to the interactions between the cascaded control loops and the complex dependence of the system dynamics on the controller parameters ，the traditional parameter tuning method is not effective at low switching frequencies.The new solution was implemented in an iteratively optimized form to ensure system stability and remove system eigenvalues from critical locations.Finally ，a time-domain comparison simulation was carried out to verify that the parameters designed with the new scheme significantly and improve the dynamic performance of the VSG compared to the traditional parameter design scheme.
Key words:inverter ；virtual synchronous generator （VSG ）；parameters tuning ；stability
基金项目：国网山东省电力公司科技资助项目（52061416000v ）
作者简介：程诗尧（1990—），男，硕士，高级工程师，Email ：*****************
摘要：在电压电流双闭环内核外级联功率控制外环构成的虚拟同步发电机（VSG ）控制方案可引入惯性到逆变器中，从而增强并网后对系统稳定性的支持。针对此种VSG 的参数设计和稳定性问题，提出了一种由线性化系统模型特征值灵敏度矩阵引导的VSG 参数整定方法。由于级联控制回路之间的相互作用以及系统动态对控制器参数的复杂依赖性，使传统的参数整定方法在低开关频率下效果欠佳。而新方案以迭代优化的形式实现，可确保系统稳定性，并使系统特性值从关键位置处移开。最后，开展了时域对比仿真，验证了相对于传统参数设计方案，采用新方案设计的参数明显改善了VSG 的动态性能。
关键词：逆变器；虚拟同步发电机；参数整定；稳定性中图分类号：TM46
文献标识码：A
DOI ：10.19457/j.1001-2095.dqcd20394
虚拟同步发电机稳定性分析与参数设计
程诗尧，王芳，李雪城，史鹏博，迟源
（国网北京市电力公司电力科学研究院，北京100075）
虚拟同步发电机（virtual synchronous genera⁃tor ，VSG ）是一种模拟同步发电机特性的控制策略[1]。VSG 可工作在并网或孤岛模式，并且能够像不间断电源一样利用下垂特性均分负载[2]。
VSG 的惯性对电力系统稳定性尤为重要[3]。
VSG 控制器的实现目前已有几种方案，VSG
首次由Beck 和Hesse 在文献[4]中提出，其实现方
案中将传统发电机的全阶模型作为控制系统的一部分。控制中电机模型输入为所测量的并网点电压，从而计算出电机电流，然后将这些电流作为电流滞环控制器的参考值来实现发电机特性模拟。文献[5]和文献[6]中设计了全阶模型的VSG 方案，
且改为电压源型设计，更准确地模拟了发电机特性。后续为了降低系统复杂度，将降阶模型代替全阶模型，不再考虑对同步发电机的电磁暂态模拟，而只是引入了惯性和阻尼[7]。文献[8]和文献[9]中的VSG 最为简单，直接采用转子方程来实现惯性，即由转子方程和无功功率控制器结合一起为PWM 生
成模块提供电压参考。该设计确保了控制器结构简单，惯性模拟和无功功率控制确定的电压参考值可以输入到电压电流双闭环级联控制器的内部。该控制方案可方便地实现电流和电压限制，以保护逆变器免受过流影响并避免意外的过调制。尽管使用电压电流双闭环级联控制器可
ELECTRIC DRIVE 2021Vol.51No.2
51
程诗尧，等：虚拟同步发电机稳定性分析与参数设计
电气传动2021年第51卷第2期方便地实现逆变器控制，但相对较多的控制器参数和各个环路之间的相互作用使参数整定变难，同时还需详细分析系统稳定性[10-11]。
对此，先介绍传统的控制器参数整定。然后，导出VSG 线性化状态空间模型，利用特征值评估小信号稳定性。分析结果表明，如果逆变器以低开关频率（如2～3kHz ）运行，则内环电流控制器的带宽受到限制，从而导致系统失稳。因此，提出了一种由线性化模型的特征值的灵敏度矩阵引导的VSG 参数调整方法。利用特征值的灵敏度可识别对临界模态影响最大的参数，同时利用迭代算法将系统特征值从临界位置移开，从而增强系统稳定性。最后，通过对比仿真验证了所提出参数整定方法的有效性。
1VSG 控制策略概览
本文的研究对象为并入电网或者带负载的
三相逆变器，逆变器输出端包含有LC 滤波器。图1为所分析的逆变器配置示意图。
本文设计的系统主要参数如下：额定容量S b =
1MV·A ，额定输出线电压有效值U S_LL_RMS =690V ，LC 滤波器电阻R 1=0.003（标幺值），LC 滤波器电容C 1=0.074（标幺值），LC 滤波器电感L 1=0.08（标幺值），VSG 虚拟惯性时间常数T a =2s ，VSG 阻尼
系数k d =10000，额定电网频率f n =50Hz ，电源电阻
R s =0.003（标幺值），电源电感L s =0.1（标幺值）。以下分析将以这些参数为基础。为了简化分析，假设直流链路上的储能单元所蕴含的能量等效于VSG 模拟的电机旋转惯性。1.1
同步参考坐标系中的电气模型
逆变器在忽略开关暂态后可由连续时间平均模型描述：ìí
îï
ï
ïïL 1ωb d i L d
d t =u L d -R 1⋅i L d +ωpu ⋅L 1⋅i L q -u O d L 1ωb d i L q
d t =u L q -R 1⋅i L q -ωpu ⋅L 1⋅i L d -u O q （1）ìíîïïïïC 1
ωb du O d
d t
=i L d +ωpu ⋅C 1⋅u O q -i O d
C 1ωb
d u O q
d t
=i L q -ωpu ⋅C 1⋅u O d -i O q
（2）
式中：L 1，C 1，R 1分别为滤波电感、滤波电容和寄生电阻；i L d ，i L q 分别为电感电流d ，q 轴分量；u L d ，u L q 分别为电感前端电压d ，q 轴分量；i O d ，i O q 分别为逆变
器输出电流的d ，q 轴分量；u O d ，u O q 分别为逆变器输出电压的d ，q 轴分量；ωb ，ωpu 分别为频率基值、
标幺值。
逆变器侧滤波电感方程如式（1）所示，滤波电容方程由式（2）给出。网侧或负载侧可由RL 电路描述，相应模型形式和式（1）一样。1.2
基于VSG 的控制系统
图2为基于VSG 的逆变器控制器结构。
从图2中可看出，控制系统含有五部分，其中控制器末端为PWM 调制模块，然后依次是基于同步参考坐标系的电流、电压控制环，无功功率下垂控制模块和VSG 惯性和阻尼模块。q 轴电压参考值等于0，同步参考坐标系定向到滤波电感电压矢量上，d 轴电压参考则是由传统的无功功率下垂控制器提供的。同步参考坐标系的相位角是根据下式给出的转子方程对惯性和阻尼模拟后输出的频率进行积分得到的。
T v ⋅s ⋅ωpu ≈P *-P -k d (ωpu -ωgpu )
（3）
式中：T v 为虚拟转矩；ωpu ，ωgpu 分别为角频率标幺
值、电网角频率标幺值；P *，P 分别为有功功率设定值和实际测量值；k d 为下垂系数；s 为算子。
如果逆变器在并网模式下运行，则VSG 通过
转子方程的功率平衡与电网电压相角同步。因此，VSG 的正常运行时不依赖于任何电网同步算
图1
所分析的逆变器配置示意图
Fig.1
Schematic of the analyzed inverter configuration
图2基于VSG 的逆变器控制器结构Fig.2
Controller structure of the VSG -based
inverter
52
程诗尧，等：虚拟同步发电机稳定性分析与参数设计电气传动2021年第51卷第2期
法，如锁相环等。与传统的发电机一样，VSG 本身也具有孤岛运行的能力。
2传统的VSG 控制器参数整定
2.1
电流闭环参数设计
为了对图2中所示电流闭环参数进行设计，
需考虑逆变器PWM 工作模式的延迟效应。采用一阶传递函数近似PWM 的影响，并假设测量电压的前馈和d ，q 轴电压理想解耦，则电流控制器的开环传递函数h cc ，dq 可定义如下：
h cc,dq ≈(k pc +k ic s )⋅11+T d ⋅s ⋅1
r (1+T 1⋅s )（4）
其中
T 1=
L 1R 1=l 1r 1⋅ωb T d ≈12⋅f sw
式中：k pc ，k ic 分别为电流闭环PI 调节器的比例、积
分系数；L 1，R 1分别为滤波电感、滤波电阻标幺值；l 1，r 1分别为滤波电感、滤波电阻实际值；T 1为滤波电感时间常数；T d 为延迟效应近似时间常数；f sw
为开关频率。
基于式（4），按照极点抵消原则来选择增益，从而产生临界阻尼闭环传递函数，对应k pc 和k ic 的计算式如下：
ì
í
î
ïï
ï
ïk pc =
l 12⋅ωb ⋅T d
k ic ≈r 12⋅T
d （5）由式（5）可知，PI 参数和带宽直接受逆变器
开关频率限制。以开环传递函数极点抵消原则设置参数得到的系统相角裕度约为65°，可进一步降低相角裕度以获取更多的控制器增益和带宽。2.2
电压闭环参数设计
从图2中的级联控制器结构可看出，电压闭环参数受到电流控制带宽的限制。为了简化，将电流控制器用一阶传递函数近似。假设d ，q 轴理想化去耦，则电压开环传递函数为下式：
h vc,dq ≈(k pu +k iu s )⋅11+T eqcc ⋅s ⋅
1
T C ⋅s （6）
其中
T eqcc =2⋅T d
T C =
c 1
ωb
式中：k pu ，k iu 分别为电压闭环PI 调节器的比例、积
分系数；C 1，c 1分别为滤波电容标幺值、实际值；T eqcc 为电流闭环时间常数；T C 为滤波电容积分时间常数。
电压闭环参数是根据确保电压开环传递函数的穿越频率处最大相角裕度的原则设计的。增加设计参数a 将电压闭环PI 调节器的参数整定和闭环传递函数阻尼因子ζ关联，如下所示：
ì
í
î
ïïï
ïk pu
≈T C a ⋅T eqcc
k iu ≈T C a 3⋅T 2
eqcc （7）
其中
a ≈2ζ+1
传统参数整定要求电压环和电流环间带宽保持较大差异。这是通过增大a 的值来保证的，直到电压闭环带宽比电流控制带宽低至少十倍。
3稳定性和参数敏感度分析
对于图2中所示的VSG 控制方案，必须显著
降低电压控制闭环的带宽，否则无法确保无功功率控制和虚拟惯性控制的预期性能以及整个控制系统的稳定性。而如果电流控制闭环带宽较低，则可能发生不同控制回路之间的相互作用。故下面对控制
系统稳定性和动态性能开展分析。3.1
状态空间系统模型
由图1和图2，可导出VSG 小信号数学模型如下式所示：
x =[Δu O d
Δu O q Δξd Δξq
Δi d
Δi q Δσd Δσq
Δi O d
Δi O q Δq m Δωpu
Δγ]T
x =éëêêêêêêêêêùû
úúúúúúúúúA 0
B 0C D 0E 0000F 0G H I J K L 0M
N E 0O P 0Q 000R S T U 000V W 00
00
00
ωb
0x （8）
其中
A =é
ëêùûú0ωb ωpu,0-ωb ωpu,00B =éëêùûúωb /c 1
00s parameter
ωb /c 1C =é
ëêùûú-ωb /c 100-ωb /c 1D =éëêùûú000-ωb u O d ,0E =éëêùûú-100-1F =é
ëêùûú0-m q 0
0G =éëêêêê
êùû
úúúúú-1-k pc k pu +k FFu l 1/ωb -k pc c 1ωpu,0l 1/ωb k pc c 1ωpu,0l 1/ωb -1-k pc k pu +k FFu l 1/ωb 53
程诗尧，等：虚拟同步发电机稳定性分析与参数设计
电气传动2021年第51卷第2期H =éëêêêêêùûúúúúúk pc k iu l 1/ωb
00k pc k iu l 1/ωb I =éëêêê
êêù
û
úúúúú
-k pc +r 1
l 1/ωb 0
0-k pc +r 1l 1/ωb J =éëêêêêêùûúúúúúk ic l 1/ωb 00k ic l 1/ωb K =éë
êêê
ê
êù
û
ú
úúúúk FF i k pc
l 1/ω
b
00k FFi k pc l 1/ωb L =éëêêêê
êù
ûú
úúúú-k pc k pu m q l 1/ωb 00k pc c 1u O d ,0l 1/ωb M =éëê
ùûú-k pu -c 1ωpu,0c 1ωpu,0-k pu N =é
ëêùûúk iu 00k iu O =éëêùûúk FF i 00k FF i P =éëêùûú-k pu m q 00c 1ωO d ,0Q =éë
êêêê
êêù
ûúúúúúú1
l 1/ωb
001l 1/ωb R =éëêêê
êêùû
úúúúú-r s l 1/ω
b ωb ωpu,0-ωb ωpu,0-r s l 1/ωb S =éëêùûúωb i O q ,0
0-ωb i O d ,00T =éëêêêêêùû
úúú
úúu S,0sin δ0l 1/ωb ωb u S sin δ0l 1/ωb U =éëêêêêùûúúúúωb ωf,pu i O q ,0-ωb ωf,pu i O d ,0-i O d ,0T a ωpu,0-i O q ,0T a ωpu,0V =éëêêêêùû
ú
ú
úú0ωb ωf,pu u O d ,0-u O d ,0T a ωpu,00W =éë
êêêêù
ûúúúú-ωb ωf,pu 00-k d T a ωg,pu 式中：“Δ”代表状态变量在工作点附近的小信号扰动；k FF u 为电流闭环中电压反馈系数；k FF i 为电压闭环中电流反馈系数；ξ，σ分别为电流和电压
控制器中积分器的状态辅助变量。下标“0”为工作点的变量值；γ为逆变器控制系统和电源相关的两个同步坐标系之间的相位差，有γ0=θVSG ，0
−ωg t 0。3.2
特征值分析和系统稳定性
基于式（8）中状态矩阵的特征值可分析VSG 稳定性。
图3为调节a 从0.5变化至50时，系统特征值在复平面中的轨迹。图3为前述传统参数整定的3种情况，但具有不同的开关频率。当系统不稳
定时，特征值为黑，图3中结果表明，当开关频率为2kHz 时，系统稳定工作范围非常窄，当开关频率增加到5kHz 时，系统稳定区间更宽，但仍对控制闭环的参数整定非常敏感。进一步将开关频率提高到10kHz ，系统具有相当宽的稳定工作范围，这意味着传统参数整定方法适用于高开关频率的低功率逆变器设计，如kW 级小规模微电网应用VSG 技术时。
图3所示的特征值轨迹分析结果表明，系统稳定性受电流闭环带宽的强烈影响，即受到逆变器开关频率影响。此外，当VSG 方案中采用相对较低的开关频率时，传统的参数整定方法可导致较差的动态响应性和低稳定裕度，因为各个级联控制闭环不能充分解耦。实际上，尽管各个控制闭环设计均为稳定的，
但它们之间的相互作用会对系统整体稳定性产生负面影响，从而导致稳定工作范围变窄。因此，有必要探索新的参数整定方法。3.3
参数灵敏度矩阵
定义系统极点的参数灵敏度为系统特征值对于系统参数的导数，从而可通过参数灵敏度辨识对系统临界特征值影响最大的参数，并可以揭示修改这些参数影响特征值位置的规律。
假设n 阶系统具有k 个可调参数，灵敏度矩阵为n ×k 阶。参数ρk 相对于特征值λn 的灵敏度αn ，k 可由下式表示：
αn ,k =
∂λn
∂ρk
=ΦT n
∂Γ
∂ρk
Ψn ΦT n Ψn
（9）
式中：Ψn ，Φn T 分别为与特征值λn 相关联的左、右特征向量；Г为n 阶系统的状态矩阵。
图3对应参数a 变化的系统特性值轨迹Fig.3
The trajectories of the system eigenvalues with variation of parameter
a
54
程诗尧，等：虚拟同步发电机稳定性分析与参数设计电气传动2021年第51卷第2期
灵敏度的实部和虚部分别与沿着实轴、虚轴的极点位置的导数相关联。由于极点位置的实部表征了系统的稳定性，故下面仅需研究灵敏度矩阵的实部。
参数灵敏度的两个例子如图4所示，系统配置见第1节中系统主要参数设置，开关频率为2kHz，设定参数a=4。
图4a为参数整定导致系统失稳的情况下不稳定特征值的分布，其中深点代表稳定特征根，浅点代表不稳定特征根。图4b为不稳定的特征值对13个参数的灵敏度，该不稳定的特征值位置主要取决于控制器参数k pc和k pu，但LC滤波器参数和电网电感参数也有一定影响，如k pc或k pu 增大，则该特征值将向左移动，从而稳定了系统。类似地，图4c、图4d绘出了接近于虚轴的另一个特征值的分析结果，图4c中浅点代表接近于虚轴的特征值，反之则为深点。在这种情况下，电网电感参数有显著影响，但特征值也受k i u和k d 的影响。
4VSG控制器参数整定算法
利用3.3节的灵敏度矩阵可整定系统参数。下面将引入VSG控制器的新型参数整定算法。参数整定基于一个迭代过程，其中每一步都将修改控制器参数将系统极点配置到一个更合适的位置。
图5为参数整定算法流程图。
由图5可知，VSG控制器的新型参数整定算法主要步骤如下：首先，将系统参数和一组初始控制器参数配置作为迭代起点，初始参数通常是采用传统整定方法得到的。然后，计算出稳态工作点。值得注意
的是，计算得到的工作点可能不稳定，但是不会影响到后续迭代过程和最终结果。接着，计算围绕该工作点的小扰动线性化系统的状态矩阵Г及其特征值。同时，辨识出具有最大实部的极点并将其选择为当前优化步骤的目标。继续计算灵敏度矩阵以得到该关键极点相对于每个可调参数的灵敏度，并到最高灵敏度对应的参数，在每次迭代中，略微改变该参数。灵敏度的符号确定了参数是增加还是减小，以便迫使目标极点向左移动。
通常，可以根据预设的终止条件终止迭代过程。为简单起见，在该示例中，迭代过程采用固定数量的迭代。示例中整定后的参数将确保最关键的极点移动到具有较小实部的位置。尽管可以考虑更高级的改变系统参数的标准，但该算法采用简单条件通常是有益的，因为造成不稳定的极点将逐渐移动到稳定区，而极点特性将通过迫使特征值向左移动而得到改善。然而，值得注意的是，某些参数的修改可能导致其他极点右移，导致算法可能表现出在每个迭代周期的不同临界极点之间不停交替振荡。还要注意的是，可调参数不仅仅限于控制器增益，还包括物理参数
图4参数敏感度分析Fig.4Parameter sensitivity analysis
图5参数整定算法流程图Fig.5Flow chart of parameter tuning
algorithm
55

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。