电工技术学报

参数	数值	参数	数值
直流参考电压U_dc^ref/V	800	滤波电感1/mH	0.5
无功电流参考值i_c_,i_,q,ref/A	0	滤波电感2/mH	0.3
开关频率/kHz	10	滤波电容/μF	10
交流电网电压u_g/V	415	采样时间/μs	1
比例系数k_p,_i_,dc	0.5	积分系数k_in,_i_,dc	25
比例系数k_p,_i_,c	3	积分系数k_in,_i_,c	40
比例系数k_p,_i_,pll	3	积分系数k_in,_i_,pll	24
电网电阻/μΩ	1	系统容量/(kV·A)	1 000
电网电感/mH	1.5	光照强度/(W/m²)	1 000

参数

数值

参数

数值

直流参考电压U_dc^ref/V

800

滤波电感1/mH

0.5

无功电流参考值i_c_,i_,q,ref/A

滤波电感2/mH

0.3

开关频率/kHz

滤波电容/μF

交流电网电压u_g/V

415

采样时间/μs

比例系数k_p,_i_,dc

0.5

积分系数k_in,_i_,dc

比例系数k_p,_i_,c

积分系数k_in,_i_,c

比例系数k_p,_i_,pll

积分系数k_in,_i_,pll

电网电阻/μΩ

系统容量/(kV·A)

1 000

电网电感/mH

1.5

光照强度/(W/m²)

1 000

参数	数值	参数	数值
直流参考电压U_dc^ref/V	800	滤波电感1/mH	0.5
无功电流参考值i_c_,i_,q,ref/A	0	滤波电感2/mH	0.3
开关频率/kHz	10	滤波电容/μF	10
交流电网电压u_g/V	415	采样时间/μs	1
比例系数k_p,_i_,dc	0.5	积分系数k_in,_i_,dc	25
比例系数k_p,_i_,c	3	积分系数k_in,_i_,c	40
比例系数k_p,_i_,pll	3	积分系数k_in,_i_,pll	24
电网电阻/μΩ	1	系统容量/(kV·A)	1 000
电网电感/mH	1.5	光照强度/(W/m²)	1 000

参数

数值

参数

数值

直流参考电压U_dc^ref/V

800

滤波电感1/mH

0.5

无功电流参考值i_c_,i_,q,ref/A

滤波电感2/mH

0.3

开关频率/kHz

滤波电容/μF

交流电网电压u_g/V

415

采样时间/μs

比例系数k_p,_i_,dc

0.5

积分系数k_in,_i_,dc

比例系数k_p,_i_,c

积分系数k_in,_i_,c

比例系数k_p,_i_,pll

积分系数k_in,_i_,pll

电网电阻/μΩ

系统容量/(kV·A)

1 000

电网电感/mH

1.5

光照强度/(W/m²)

1 000

参数	数值	参数	数值
直流参考电压 $U dc ref /V$	800	滤波电感1/mH	0.5
无功参考电流i_c_,s_,q,ref/A	50	滤波电感2/mH	0.3
开关频率/kHz	10	滤波电容/μF	10
交流电网线电压u_g/V	415	采样时间/μs	1
比例系数k_p,_s_,dc	5	积分系数k_in,_s_,dc	20
比例系数k_p,_s_,c	0.3	积分系数k_in,_s_,c	40
比例系数k_p,_s_,pll	3	积分系数k_in,_s_,pll	40
电网电阻/μΩ	1	系统容量/(kV·A)	1 000
电网电感/mH	1.5	光照强度/(W/m²)	1 000

参数

数值

参数

数值

直流参考电压

U dc ref /V

800

滤波电感1/mH

0.5

无功参考电流i_c_,s_,q,ref/A

滤波电感2/mH

0.3

开关频率/kHz

滤波电容/μF

交流电网线电压u_g/V

415

采样时间/μs

比例系数k_p,_s_,dc

积分系数k_in,_s_,dc

比例系数k_p,_s_,c

0.3

积分系数k_in,_s_,c

比例系数k_p,_s_,pll

积分系数k_in,_s_,pll

电网电阻/μΩ

系统容量/(kV·A)

1 000

电网电感/mH

1.5

光照强度/(W/m²)

1 000

参数	数值	参数	数值
直流参考电压 $U dc ref /V$	800	滤波电感1/mH	0.5
无功参考电流i_c_,s_,q,ref/A	50	滤波电感2/mH	0.3
开关频率/kHz	10	滤波电容/μF	10
交流电网线电压u_g/V	415	采样时间/μs	1
比例系数k_p,_s_,dc	5	积分系数k_in,_s_,dc	20
比例系数k_p,_s_,c	0.3	积分系数k_in,_s_,c	40
比例系数k_p,_s_,pll	3	积分系数k_in,_s_,pll	40
电网电阻/μΩ	1	系统容量/(kV·A)	1 000
电网电感/mH	1.5	光照强度/(W/m²)	1 000

参数

数值

参数

数值

直流参考电压

U dc ref /V

800

滤波电感1/mH

0.5

无功参考电流i_c_,s_,q,ref/A

滤波电感2/mH

0.3

开关频率/kHz

滤波电容/μF

交流电网线电压u_g/V

415

采样时间/μs

比例系数k_p,_s_,dc

积分系数k_in,_s_,dc

比例系数k_p,_s_,c

0.3

积分系数k_in,_s_,c

比例系数k_p,_s_,pll

积分系数k_in,_s_,pll

电网电阻/μΩ

系统容量/(kV·A)

1 000

电网电感/mH

1.5

光照强度/(W/m²)

1 000

d通道	q通道
p₁	-8.032+0.003i	p₂₂	-8.032+0.006i
p₂₃	-8.032-0.006i
p₂	-8.032-0.003i	p₅₀	-8.038
p₅₁	-8.025
p₂₁	-3457.421	p₅₂	-4.812
p₅₃	-4.811

d通道

q通道

p₁

-8.032+0.003i

p₂₂

-8.032+0.006i

p₂₃

-8.032-0.006i

p₂

-8.032-0.003i

p₅₀

-8.038

p₅₁

-8.025

p₂₁

-3457.421

p₅₂

-4.812

p₅₃

-4.811

d通道	q通道
p₁	-8.032+0.003i	p₂₂	-8.032+0.006i
p₂₃	-8.032-0.006i
p₂	-8.032-0.003i	p₅₀	-8.038
p₅₁	-8.025
p₂₁	-3457.421	p₅₂	-4.812
p₅₃	-4.811

d通道

q通道

p₁

-8.032+0.003i

p₂₂

-8.032+0.006i

p₂₃

-8.032-0.006i

p₂

-8.032-0.003i

p₅₀

-8.038

p₅₁

-8.025

p₂₁

-3457.421

p₅₂

-4.812

p₅₃

-4.811

通道	参与因子
d	p_k1	1.28×10^-12
p_k2	1.28×10^-12
p_k17	3.32×10^-13
q	p_k22	6.25×10^-12
p_k23	6.25×10^-12
p_k29	9.05×10^-12
p_k30	9.05×10^-12
p_k31	3.32×10^-13
p_k32	3.35×10^-13

通道

参与因子

p_k1

1.28×10^-12

p_k2

1.28×10^-12

p_k17

3.32×10^-13

p_k22

6.25×10^-12

p_k23

6.25×10^-12

p_k29

9.05×10^-12

p_k30

9.05×10^-12

p_k31

3.32×10^-13

p_k32

3.35×10^-13

通道	参与因子
d	p_k1	1.28×10^-12
p_k2	1.28×10^-12
p_k17	3.32×10^-13
q	p_k22	6.25×10^-12
p_k23	6.25×10^-12
p_k29	9.05×10^-12
p_k30	9.05×10^-12
p_k31	3.32×10^-13
p_k32	3.35×10^-13

通道

参与因子

p_k1

1.28×10^-12

p_k2

1.28×10^-12

p_k17

3.32×10^-13

p_k22

6.25×10^-12

p_k23

6.25×10^-12

p_k29

9.05×10^-12

p_k30

9.05×10^-12

p_k31

3.32×10^-13

p_k32

3.35×10^-13

含静止同步补偿器的光伏并网系统静态同步稳定分析

PDF下载

马睿聪 ¹ , 曹永吉 ² , 张恒旭 ¹ , 李常刚 ¹

电工技术学报 | 电力系统与综合能源 2025,40(11): 3410-3426

收起

电工技术学报 | 电力系统与综合能源 2025, 40(11): 3410-3426

含静止同步补偿器的光伏并网系统静态同步稳定分析

全屏

马睿聪¹, 曹永吉², 张恒旭¹, 李常刚¹

作者信息

¹ 山东大学电气工程学院济南 250061

² 山东大学智能创新研究院济南 250101

马睿聪男，1999年生，博士研究生，研究方向为电力系统同步稳定分析与紧急控制。E-mail：202114573@mail.sdu.edu.cn

通讯作者:

曹永吉男，1992年生，副研究员，硕士生导师，研究方向为电力系统稳定分析与控制、可再生能源并网及储能技术应用。E-mail：yongji@sdu.edu.cn

Static Synchronization Stability Analysis of PV Grid-Connected System with Static Synchronous Compensator

Ruicong Ma¹, Yongji Cao², Hengxu Zhang¹, Changgang Li¹

Affiliations

¹ School of Electrical Engineering Shandong University Jinan 250061 China

² Academy of Intelligent Innovation Shandong University Jinan 250101 China

出版时间: 2025-06-10 doi: 10.19595/j.cnki.1000-6753.tces.240754

文章导航

摘要

收起

光伏电站并网时需要具备一定的无功支撑能力，通常需要配备无功补偿装置，二者之间的相互作用增加了静态同步稳定问题的分析难度。该文首先提出一种含静止同步补偿器（STATCOM）的光伏并网系统静态同步稳定分析方法；然后考虑不同并网变换器dq坐标系间的变换，构建多并网变换器并联系统的等效小信号模型，在阻抗分析方法基础上改进稳定判据，并利用衰减系数指标定量刻画系统的静态同步稳定性能，构建静态同步稳定量化评估指标体系；最后结合改进的判据和量化指标，详细地分析了STATCOM对光伏并网系统静态同步稳定性能的影响，总结了各控制参数对于系统运行零极点分布的影响规律。

关键词

光伏系统 / 静止同步补偿器（STATCOM） / 静态同步稳定 / 阻抗分析方法 / 衰减系数

Abstract

收起

PV systems are typically equipped with reactive power compensation devices when connected to the grid, and static synchronous compensator (STATCOM) devices are widely employed due to their flexible control capabilities. The increasing utilization of power electronic devices in the power grid has resulted in a shift from physical synchronization to control synchronization as the dominant mode of system operation. Analyzing static synchronization stability problem is more challenging for these systems compared to conventional power systems, as converter output characteristics are influenced by control strategies. Therefore, it is imperative to urgently address the problem of static synchronization stabilization under control strategy dominance.

First, this paper establishes the static synchronous stability analysis model of the grid-connected converter based on the control loop and circuit structure of each converter in a parallel system under the respective dq reference frame. The dq reference frame of the converter is determined by the phase information provided by the control loop in a multiple converter parallel system, thus enabling a unified coordinate system for static synchronization stability analysis. Subsequently, an equivalent small signal model is developed for analyzing multiple grid-connected converters. In comparison with existing coordinate conversion methods, Kirchhoff's current law is incorporated to enhance accuracy and reduce errors.

Then, the stability criterion for impedance analysis is enhance, and the static grid-synchronization performance indices are created. The small perturbation oscillation characteristics are measured using overshooting and regulation time, while the participation factor is employed to analyze the impact of each pole of the system. Finally, the attenuation coefficient is utilized to assess static synchronization stability performance.

The model is developed in Matlab/Simulink for simulation verification. Subsequently, an analysis is conducted on the impact of parameters such as the control loop parameters and STATCOM capacity on static synchronization stability. The main conclusions are summarized as follows:

(1) Optimizing the reactive power output of the grid-connected converter based on known control parameters can significantly enhance static synchronization stability performance, with the dominant influence of small perturbations after oscillation mode being attributed to poles generated by phase-locked loop control.

(2) The attenuation coefficient of the system exhibits a rapid increase in proximity to the critical stability region. Hence, it is imperative for the system to possess a certain margin of attenuation coefficient during operation. Based on the simulation analysis results presented in this study, static synchronous instability phenomena occur when the attenuation coefficient of the PV system exceeds 300. Conversely, when the attenuation coefficient falls below 250, the system remains in a state of static synchronous stability. These findings establish a criterion for analyzing and assessing static synchronization stability within such systems.

(3) The addition of STATCOM to the PV system primarily impacts the conductance matrix transfer function of the q-coupled channel. The phase-locked-loop coupling oscillations between the grid-connected converters do not affect the dd channel. Within the stable operating region, an increase in bandwidth for the DC voltage control loop, active current control loop, and reactive current control loop results in an amplification of both attenuation coefficient and system oscillation amplitude.

(4) When the phase-locked loop parameters of the STATCOM are the same as the PV system, the stability performance of the system is mainly affected by the grid impedance and the equivalent conductance transfer function of each grid-connected converter. And each grid-connected converter can independently connect to the grid and achieve stable operation to ensure that the system achieves static synchronous stability in this case.

(5) In cases where the active output of the PV system is low, STATCOM typically adjusts its capacitive or inductive reactive power provision to improve static synchronization stability performance. Conversely, when compensating for capacitive reactive power, utilizing STATCOM may yield superior results compared to using the PV grid-connected converter alone. Hence, allocating an optimal capacity for STATCOM can significantly enhance static synchronization stability performance.

Key words

PV system / static synchronous compensator (STATCOM) / static synchronization stability / impedance analysis method / attenuation coefficient

引用本文

马睿聪, 曹永吉, 张恒旭, 李常刚. 含静止同步补偿器的光伏并网系统静态同步稳定分析. 电工技术学报, 2025 , 40 (11) : 3410 -3426 . DOI: 10.19595/j.cnki.1000-6753.tces.240754

Ruicong Ma, Yongji Cao, Hengxu Zhang, Changgang Li. Static Synchronization Stability Analysis of PV Grid-Connected System with Static Synchronous Compensator[J]. Transactions of China Electrotechnical Society, 2025 , 40 (11) : 3410 -3426 . DOI: 10.19595/j.cnki.1000-6753.tces.240754

正文

收起

0 引言

收起

电力系统内各动态元件处于相对静止的安全稳定运行状态，被称为电力系统处于同步状态^[1]。风电和光伏等新能源采用电力电子装置作为并网接口，具有非线性、异构性和低惯量等明显特征，导致系统同步方式从以同步发电机为主导的物理同步转变为以并网变换器为主导的控制同步，引起了新型电力系统中的同步机理转变^[2-3]。随着新能源装机规模的增加，电网传统物理同步方式的主导特性持续减弱，并网变换器多控制回路间的强耦合特性对电网影响日益显著，新型电力系统静态同步稳定问题机理更加复杂^[4-5]。相较于传统电力系统，新型电力系统静态同步稳定问题分析的难点在于并网变换器的应用使得输出特性与控制策略密切相关，亟须对控制策略主导下的新能源静态同步稳定问题开展分析^[6-7]。

现有研究大多采用阻抗分析方法评估新能源发电系统的静态同步稳定特性^[8]。文献[9]针对多并网变换器系统，给出详细的阻抗分析建模过程，分析系统的稳定性，采用仿真验证了所提方法的准确性。文献[10]针对线性变参数（Linear Parameter-Varying, LPV）控制的功率同步跟网型变流器（Power-Synchronized Grid Following Inverter, PSGFLI），构建输出的阻抗矩阵，分析控制参数运行对系统稳定性能的影响。文献[11]构建含储能多馈入系统的交流电网动态阻抗模型，基于多变量反馈理论实现系统的模态解耦。文献[12]构建模块化多电平变换器（Modular Multilevel Converter, MMC）多输入多输出模型，采用块对角占优理论实现模型降维。阻抗分析方法的难点在于稳定判据的选取，要进一步拓展阻抗分析稳定判据^[13-16]。

影响系统静态稳定的因素主要有以下几点：

1）系统的静态同步稳定性受电压影响，在电压较低时，并网变换器输出的电流值上限较低，限制并网变换器输出功率及静态储备系数，不利于系统静态稳定；在发生电压跌落时，并网变换器输出功率急剧下降，系统易同步失稳^[17-18]。

2）系统的静态同步稳定性受电网阻抗的影响，在并网变换器接入系统后，较高的电网阻抗值会导致系统整体等效阻尼较小，甚至呈现负阻尼，引发低频振荡^[19]。利用电网短路比（Short Circuit Ratio, SCR）表征并网变换器接入电网的相对强度，当构网型（Grid Forming, GFM）装置接入点SCR在1附近时，系统不一定存在稳定工作点，因此短路比较小的系统更容易发生静态同步稳定失稳^[20-21]。

3）系统的静态同步稳定性受并网变换器控制方式及参数的影响，利用复转矩系数法分析可得到阻尼传递函数D(s)，增加新的控制环节能够增加系统的静态同步稳定裕度^[22-23]。文献[24]揭示了各控制参数与特征根间的关系，改变系统的控制参数能够提高扰动振荡过程的衰减速度。文献[25]针对锁相环（Phase Locked Loop, PLL）控制引入线性反馈补偿器，使得跟网型（Grid Following, GFL）装置的静态稳定能力不依赖初始的PLL参数设计，通过改变线性反馈补偿器，实现GFL装置在弱电网下的静态运行能力的提升。文献[26]提出一种基于相位补偿器的改进PLL控制策略，有效地改善了GFM装置在弱电网下的静态稳定性能。

虽然业界已对并网变换器与电网交互系统的静态同步稳定问题开展了系统研究，但对并网变换器间交互机理的研究还不够深入，未区分系统控制环节产生的主导特征值。在实际工程中，新能源机组并网运行时会增设无功补偿装置，以确保并网变换器输出电压稳定^[27-28]。但是现有研究尚未分析静止同步补偿器（Static Synchronous Compensator, STATCOM）对光伏发电系统静态稳定机理的影响，缺乏含静止同步补偿器STATCOM的光伏并网系统静态同步稳定判据。相较于现有的研究，本文考虑STATCOM锁相环及电流控制环节对光伏系统静态同步稳定的影响，提出了适用于多并联变换器系统的静态同步稳定分析方法，构建了静态同步稳定判别。相较于现有研究中线性化的参考坐标系确定方法，本文结合基尔霍夫电流定律，建立了电压参考坐标系转换传递函数，减少了分析时的旋转误差。

1 光伏并网系统和STATCOM模型

收起

1.1 光伏并网系统模型

在变换器同步稳定分析方面，现有研究一般假设并网输出电压等于参考信号，主要考虑控制环节的影响，忽略变换器本身特性对系统同步稳定的影响^[8,29]。因此，级联H桥电路、三电平和两电平等不同类型的并网变换器的分析方法和结果是相同的。本文以两电平并网逆变器为例进行分析，光伏发电系统模型如图1所示。其中，光伏由侧Boost两级式升压电路组成，光伏功率输出模式为最大功率跟踪（Maximum Power Point Tracking, MPPT）控制，并网逆变器采用电压电流内外双环控制。

新能源机组并网控制信号输出主要由电流环和锁相环调节，其余环节通过影响电流环参考值来改变并网变换器的控制信号。多并网变换器系统如图2所示，假设系统中光伏并网变换器数量为M，STATCOM变换器数量为N，总变换器数量为K。

在忽略端口电阻的情况下，并网变换器i端口处输出电压u_t_,i,_abc和电流i_t_,i,_abc间的关系为

（1）

u t, i, abc = L t, i, 1 d i c, i, abc d t + L t, i, 1 C t, i, 1 d 2 u c, i, abc d t + u c, i, abc

（2）

u t, s, abc = L t, s, 1 d i c, s, abc d t + L t, s, 1 C t, s, 1 d 2 u c, s, abc d t + u c, s, abc

（3）

u c, s, abc = u c, i, abc

（4）

i g, abc = i c, i, abc + i c, s, abc

式中，下标i为光伏变换器的序号，i=1,2,…,N；下标s为STATCOM装置变换器序号，s=1,2, …,M；L_t_,i,₁和L_t_,s,₁分别为并网变换器i和s网侧滤波电感；C_t,_i_,1和C_t,_s_,1分别为光伏并网变换器i和s滤波电容；u_c_,i,_abc为并网变换器i网侧滤波电容端电压；u_c_,s,_abc为STATCOM并网变换器s滤波电容端电压；i_g_,_abc为电网输入电流；i_c_,i,_abc为光伏并网变换器i输出电流；i_f_,i,_abc为光伏并网变换器i电容电流；i_c_,s,_abc为STATCOM并网变换器s输出电流。

光伏变换器i采用双闭环电压电流矢量控制，实现输出功率的解耦。其中，d轴电流控制采用定直流电压控制，q轴电流控制采用定无功电流控制，并网变换器控制策略和锁相环控制环节分别如图3和图4所示。图中，u_t,_i_,d,ref、u_t,_i_,q,ref分别为光伏系统并网变换器i输出d轴、q轴电压控制信号；i_i_,d,ref、i_i_,q,ref分别为变换器i输入d轴、q轴电流控制信号；ω₀为电网稳定运行状态下角频率；k_p,_i_,c为并网变换器i电流环比例系数；k_in,_i_,c为并网变换器i电流环积分系数；u_i_,dc,ref为并网变换器i直流侧电压参考值；u_i_,dc为并网变换器i直流侧电压值；k_p,_i_,dc为并网变换器i电压环比例系数；k_in,_i_,dc为并网变换器i电压环积分系数；θ_pll,_i为并网变换器i输出相位；ω_pll,_i为并网变换器i锁相环输出频率。k_p,_i_,pll为并网变换器i锁相环比例系数；k_in,_i_,pll为并网变换器i锁相环积分系数。

1.2 STATCOM模型

STATCOM具备快速调节性能，其拓扑结构如图5所示。STATCOM并网侧的结构与新能源机组相同，其输出电压u_t_,s,_abc和电流i_t_,s,_abc间的关系同样与新能源机组相同。

STATCOM锁相环控制方程与光伏并网变换器锁相环控制分析相同，假设在运行过程中无功输出不超过容量限制，其控制策略如图6所示。图中u_t,_s_,d,ref、u_t,_s_,q,ref分别为STATCOM的变换器s输出d轴、q轴电压控制信号；i_s_,d,ref、i_s_,q,ref分别为变换器s输入d轴、q轴电流控制信号；k_p,_s_,c、k_in,_s_,c、u_s_,dc,ref、u_s_,dc、k_p,_s_,dc和k_in,_s_,dc分别为STATCOM的电流环比例系数、电流环积分系数、直流电压参考值、直流电压值、电压环比例系数和电压环积分系数。

2 静态同步稳定建模及分析

收起

2.1 静态同步稳定分析模型

2.1.1 光伏系统信号建模

当锁相环输出功角发生变化时，Park变换的线性化公式可表示为

（5）

T park = T abc / dq + Δ θ pll, i ω pll, i T ˙ abc / dq

（6）

T abc / dq = 23 cos θ cos θ − 2 π 3 cos θ + 2 π 3 − sin θ − sin θ − 2 π 3 − sin θ + 2 π 3 121212

将式（1）和式（2）转换至dq坐标系下，线性化可得

（7）

Δ u t, i, d = Δ u c, i, d + s L t, i, 1 (Δ i c, i, d + Δ i f, i, d) − ω pll, i, 0 L t, i, 1 (Δ i c, i, q + Δ i f, i, q) + ω pll, i, 0 L t, i, 1 (i c, i, d, 0 + i f, i, d, 0) Δ θ pll, i − (u t, i, q, 0 − u c, i, q, 0 ︸ 近 似 为 0) Δ θ pll, i

（8）

Δ u t, i, q = Δ u c, i, q + s L t, i, 1 (Δ i c, i, q + Δ i f, i, q) + ω pll, i, 0 L t, i, 1 (Δ i c, i, d + Δ i f, i, d) + ω pll, i, 0 L t, i, 1 (i c, i, q, 0 + i f, i, q, 0) Δ θ pll, i +(u t, i, d, 0 − u c, i, d, 0 ︸ 近 似 为 0) Δ θ pll, i

式中，ω_pll,_i_,0为并网变换器i锁相环稳态输出频率；u_c,_i_,d,0和u_c,_i_,q,0分别为并网点处d轴和q轴稳态电压；u_t,_i_,d,0和u_t,_i_,q,0分别为变换器i端口处d轴和q轴稳态电压；i_c,_i_,d,0和i_c,_i_,q,0分别为并网点处d轴和q轴稳态电流；i_f,_i_,d,0和i_f,_i_,q,0分别为滤波电容的d轴和q轴稳态电流；Δu_c,_i_,d和Δu_c,_i_,q分别为并网点处d轴和q轴线性化电压；Δi_c,_i_,d和Δi_c,_i_,q分别为并网点处d轴和q轴线性化电流；Δu_t,_i_,d和Δu_t,_i_,q分别为变换器i端口处d轴和q轴线性化电压；Δi_f,_i_,d和Δi_f,_i_,q分别为滤波电容的d轴和q轴线性化电流。

对滤波电容电流进行分析，可得

（9）

$Δ i f, i, d = C t, i, 1 s Δ u c, i, d − ω pll, i, 0 C t, i, 1 Δ u c, i, q + (ω pll, i, 0 C t, i, 1 u c, i, d, 0 − i i, q, 0) Δ θ pll, i$

（10）

$Δ i f, i, q = C t, i, 1 s Δ u c, i, q + ω pll, i, 0 C t, i, 1 Δ u c, i, d + (ω pll, i, 0 C t, i, 1 u c, i, q, 0 − i i, d, 0) Δ θ pll, i$

由于滤波器电容通常较小，L_t,_i_,1、C_t,_i_,1的值远小于其他参数。因此，将式（9）和式（10）分别代入式（7）和式（8），并转换至频域，整理可得

（11）

$Δ u t, i, d = (1 − ω pll, i, 0 2 L t, i, 1 C t, i, 1) Δ u c, i, d + L t, i, 1 s Δ i c, i, d − ω pll, i, 0 L t, i, 1 Δ i c, i, q − (2 ω pll, i, 0 L t, i, 1 C t, i, 1 s) Δ u c, i, q − (L t, i, 1 i c, i, q, 0 + L t, i, 1 C t, i, 1 u c, i, q, 0 s) Δ ω pll, i$

（12）

$Δ u t, i, q = (1 − ω pll, i, 0 2 L t, i, 1 C t, i, 1) Δ u c, i, q + L t, i, 1 s Δ i c, i, q + ω pll, i, 0 L t, i, 1 Δ i c, i, d + (2 ω pll, i, 0 L t, i, 1 C t, i, 1 s) Δ u c, i, d + (L t, i, 1 i c, i, d, 0 + L t, i, 1 C t, i, 1 u c, i, d, 0 s) Δ ω pll, i$

光伏并网变换器的功率关系可表示为

（13）

12 C i, dc U i, dc 2 = P out − P in

式中，C_i_,dc为稳压电容；P_out为新能源机组输出功率；P_in为输入电网功率。

在小扰动过程中，假设P_out不发生变化，可得

（14）

C i, dc U i, dc, 0 Δ U i, dc = 32 (u c, i, d, 0 Δ i c, i, d + i c, i, d, 0 Δ u c, i, d + u c, i, q, 0 ︸ 此 项 为 0 Δ i c, i, q + i c, i, q, 0 Δ u c, i, q)

对光伏并网系统的电流控制环节进行线性化，并转换至频域，整理可得

（15）

$Δ u t, i, d, ref = − Δ U i, dc G i, c (s) k p, i, dc + k in, i, dc s − Δ i c, i, d G i, c (s) − ω 0 L t, i, 1 Δ i c, i, q + Δ u c, i, d$

（16）

Δ u t, i, q, ref = − k p, i, c + k in, i, c s Δ i c, i, q + ω 0 L t, i, 1 Δ i c, i, d + Δ u c, i, q

式中，G_i_,c(s)和G_i_,dc(s)分别为电流控制环节和直流电压控制环节PI控制函数，表达式为

$G i, c (s) = k p, i, c + k in, i, c s$

$G i, dc (s) = k p, i, dc + k in, i, dc s$

对锁相环控制环节进行线性化，并转换至频域，整理可得

（17）

Δ θ pll, i = 1 s G i, pll (s) Δ u c, i, q

（18）

Δ ω pll, i = G i, pll (s) Δ u c, i, q

式中，G_i_,pll(s)为锁相环PI控制函数。

$G i, pll (s) = k p, i, pll + k in, i, pll s$

整理式（11）、式（12）、式（15）、式（16）和式（18）可得

（19）

Δ i c, i, d = G i, 1 (s) Δ u c, i, d + G i, 2 (s) Δ u c, i, q

（20）

Δ i c, i, q = − G i, 3 (s) Δ u c, i, d + G i, 4 (s) Δ u c, i, q

式中，G_i_,1(s)、G_i_,2(s)、G_i_,3(s)、G_i_,4(s)、g_i_,1(s)、g_i_,2(s)、g_i_,3(s)、g_i_,4(s)、g_i_,5(s)和g_i_,6(s)为中间变量，且

$G i, 1 (s) = ω pll, i, 0 2 L t, i, 1 C t, i, 1 − g i, 2 (s) L t, i, 1 s + g i, 1 (s)$

$G i, 2 (s) = g i, 3 (s) L t, i, 1 s + g i, 1 (s)$

$G i, 3 (s) = 2 ω pll, i, 0 L t, i, 1 C t, i, 1 s L t, i, 1 s + g i, 5 (s)$

$G i, 4 (s) = − g i, 6 (s) G pll, i (s) − ω pll, i, 0 2 L t, i, 1 C t, i, 1 L t, i, 1 s + g i, 5 (s)$

$g i, 1 (s) = G i, c (s) + 3 2 C i, dc U i, dc, 0 G i, dc (s) G i, c (s) u c, i, d, 0$

$g i, 2 (s) = 3 2 C i, dc U i, dc, 0 G i, dc (s) G i, c (s) i c, i, d, 0$

$g i, 3 (s) = g i, 4 (s) G i, pll (s) + 2 ω pll, i, 0 L t, i, 1 C t, i, 1 s − 3 2 C i, dc U i, dc, 0 G i, dc (s) G i, c (s) i c, i, q, 0$

$g i, 4 (s) = L t, i, 1 i c, i, q, 0 + L t, i, 1 C t, i, 1 u c, i, q, 0 s$

$g i, 5 (s) = G i, c (s)$

$g i, 6 (s) = L t, i, 1 i c, i, d, 0 + L t, i, 1 C t, i, 1 u c, i, d, 0 s$

2.1.2 STATCOM小信号建模

假设系统无功功率由各STATCOM平均补偿，且STATCOM在稳态时不提供有功功率，线性化整理可得

（21）

C s, dc U s, dc, 0 Δ U s, dc = 32 [u c, s, d, 0 Δ i c, s, d + i c, s, d, 0 ︸ 此 项 为 0 Δ u c, s, d + u c, s, q, 0 ︸ 此 项 为 0 Δ i c, s, q + i c, s, q, 0 Δ u c, s, q]

对STATCOM电流控制环节进行线性化，并转换至频域，整理可得

（22）

$Δ u t, s, d, ref = G s, c (s) G s, dc (s) [(− Δ U s, dc) − Δ i c, s, d] − ω 0 L t, s, 1 Δ i c, s, q + Δ u c, s, d$

（23）

Δ u t, s, q, ref = 1 N G s, c (s) ∑ Δ i c, i, q, s − Δ i s, q + ω 0 L t, s, 1 Δ i c, s, d + Δ u c, s, q

式中，Δi_c,_i_,q_,s为并网变换器i转换至并网变换器s的dq坐标系下的q轴并网点电流。

STATCOM的电路及锁相环控制方程的线性化方式与光伏并网系统相同，结合式（22）和式（23），整理可得d轴电流表达式为

（24）

Δ i c, s, d = G s, 1 (s) Δ u c, s, d + G s, 2 (s) u c, s, q

式中，G_s_,1(s)和G_s_,2(s)为中间变量，且

$G s, 1 (s) = ω pll, s, 0 2 L t, s, 1 C t, s, 1 − g s, 2 (s) L t, s, 1 s + g s, 1 (s)$

$G s, 2 (s) = g s, 3 (s) L t, s, 1 s + g s, 1 (s)$

在不同旋转坐标系下的系统电流可表示为

（25）

$Δ i c, i, d, r = Δ i c, i, d cos (θ pll, i, 0 − θ pll, r, 0) − Δ i c, i, q sin (θ pll, i, 0 − θ pll, r, 0) − i c, i, d, 0 sin (θ pll, r, 0 − θ pll, r, 0) (Δ θ pll, i − Δ θ pll, r) − i c, i, q, 0 cos (θ pll, i, 0 − θ pll, r, 0) (Δ θ pll, i − Δ θ pll, r)$

式中，Δi_c,_i_,d_,r为并网变换器i转换至并网变换器r的dq坐标系下的d轴并网点电流值。

在稳态运行状态下，各并网变换器锁相环输入参考电压相同，输出相角相同，因此有θ_i_,0=θ_r_,0，代入式（25）整理可得

（26）

Δ i c, i, d, r = Δ i c, i, d − i c, i, q, 0 s G i, pll (s) Δ u c, i, q + i c, i, q, 0 s G r, pll (s) Δ u c, r, q

同理可得

（27）

Δ i c, i, q, r = Δ i c, i, q + i c, i, d, 0 s G i, pll (s) Δ u c, i, q − i c, i, d, 0 s G r, pll (s) Δ u c, r, q

式中，Δi_c,_i_,q,_r为并网变换器i转换至并网变换器r的dq坐标系下的d轴并网点电流值。

因此，在STATCOM旋转坐标系下，整理可得

（28）

$Δ i c, s, q = − H s, 1 (s) Δ u c, i, d + H s, 2 (s) Δ u c, i, q + T 1 (s) ∑ i = 1 M Δ i c, i, q + ∑ i = 1 M i c, i, d, 0 s G i, pll (s) Δ u c, i, q$

式中，T₁(s)、H_s_,1(s)、H_s_,2(s)为中间变量，且

$T 1 (s) = 1 N G s, c (s) L t, s, 1 s + 1 N G s, c (s)$

$H s, 1 (s) = 2 ω pll, s, 0 C t, s, 1 L t, s, 1 s L t, s, 1 s + 1 N G s, c (s)$

$H s, 2 (s) = − 1 N G s, c (s) ∑ i = 1 M i c, i, d, 0 s G s, pll (s) L t, s, 1 s + 1 N G s, c (s) + ω pll . s .0 2 L t, s, 1 C t, s, 1 − g s, 7 (s) L t, s, 1 s + 1 N G s, c (s)$ $g s, 7 (s) = (L t, s, 1 i c, s, d, 0 + L t, s, 1 C t, s, 1 u c, s, d, 0 s) G s, pll (s)$

2.1.3 阻抗分析模型

基于基尔霍夫电流定律，对并网点处的电路方程分析可得电压旋转表达式为

（29）

Δ u c, i, q = N r (s) N i (s) Δ u c, r, q

（30）

Δ u c, i, d = Δ u c, r, d + M r (s) − M i (s) N r (s) N i (s) Δ u c, r, q

式中，Δu_c,_i_,d_,r为并网变换器i转换至并网变换器r在dq坐标系下的d轴并网点电压值；M_r(s)和N_r(s)为中间变量，且

$M r (s) = (s L g + R g) ∑ i = 1 K i c, i, q, 0 s G r, pll (s) + ω pll, r, 0 L g ∑ i = 1 K i c, i, d, 0 s G r, pll (s) − u g sin (δ 0, r) s + L g i g, q, 0, r G r, pll (s)$

$N r (s) = (s L g + R g) ∑ i = 1 K i c, i, d, 0 G r, pll (s) s + ∑ i = 1 K ω pll, r, 0 L g i c, i, q, 0 G r, pll (s) s − u g cos (δ 0, r) s − L g i g, d, 0, r G r, pll (s) − 1$

由式（28）~式（30）可得，STATCOM的q轴输出电流可表示为

（31）

Δ i s, q = − G s, 3 (s) Δ u c, s, d + G s, 4 (s) Δ u c, s, q

式中，G_s_,3(s)和G_s_,4(s)为中间变量，且

$G s, 3 (s) = H s, 1 (s) + T 1 (s) ∑ i = 1 M G i, 3 (s) G s, 4 (s) = H s, 2 (s) − T 1 (s) ∑ i = 1 M G i, 3 (s) M s (s) − M i (s) N s (s) N i (s) + T 1 (s) ∑ i = 1 M G i, 4 N s (s) N i (s) + ∑ i = 1 M i i, d, 0 s G i, pll (s) N s (s) N i (s)$

以并网变换器所在dq轴为参考系，可得并网点输出电流可表示为

（32）

Δ i c, d, 1 = W 11 (s) Δ u c, 1, d + W 12 (s) Δ u c, 1, q

（33）

Δ i c, q, 1 = W 21 (s) Δ u c, 1, d + W 22 (s) Δ u c, 1, q

式中，W₁₁(s)、W₁₂(s)、W₂₁(s)和W₂₂(s)分别为光伏并网系统的dd、dq、qd和qq通道阻抗模型，且

$W 11 (s) = ∑ i = 1 K G i, 1 (s) W 12 (s) = ∑ i = 1 K G i, 2 (s) N 1 (s) N i (s) + ∑ i = 1 K i i, q, 0 s G 1, pll (s) − ∑ i = 1 K i i, q, 0 N 1 (s) s N i (s) G i, pll (s) + ∑ i = 1 K G i, 1 (s) M 1 (s) − M i (s) N 1 (s) N i (s)$

$W 21 (s) = − ∑ i = 1 K G i, 3 (s) W 22 (s) = ∑ i = 1 K G i, 4 (s) N 1 (s) N i (s) − ∑ i = 1 K i i, d, 0 s G 1, pll (s) + ∑ i = 1 K i i, d, 0 N 1 (s) s N i (s) G i, pll (s) − ∑ i = 1 K G i, 3 (s) M 1 (s) − M i (s) N 1 (s) N i (s)$

系统的输出阻抗阵可表示为

（34）

Y inv = W 11 (s) W 12 (s) W 21 (s) W 22 (s)

同理可得，电网等效输出矩阵可表示为

（35）

Z g = − s L g + R g ω 0 L g − ω 0 L g − s L g + R g

式中，L_g为电网等效电感；R_g为电网等效阻抗；ω₀为并网点处电网的稳态频率。

2.2 静态同步稳定分析判据

GFL变流器并网系统的等效电路如图7所示。

由图7可得

（36）

I g (s) = I c (s) − U g (s) Z inv (s) 1 1 + Z g (s) / Z inv (s)

假设并网变换器和电网两个子系统都能够独立稳定运行，由式（36）可得，系统的稳定性取决于1+Z_g(s)/Z_inv(s)，该环节可等效为负反馈传递函数，其中，G(s)=1，H(s)=Z_g(s)/Z_inv(s)。

在dq坐标系下，GFL变流器并网系统还应满足

（37）

Z L (s) I c, dq (s) − U c, dq (s) Z inv (s) = U c, dq (s) − U g, dq (s)

因此，以等效电压源分析得到的稳定判据可表示为

（38）

U c, dq (s) = E − Z g (s) Y inv (s) − 1 ︸ K (s) U g, dq (s)

式中，E为单位矩阵。

由式（38）可得，系统满足静态同步稳定运行的条件为等效传递函数K(s)的特征根的实数部均为负数。

2.3 静态同步稳定量化指标

在对系统的静态同步稳定进行分析时，还需要分析系统在小扰动下的暂态性能，引入超调量σ和调节时间t_s两项暂态性能指标，可表示为

（39）

σ = c max − c (∞) c (∞) × 100 %

（40）

c (t) − c (∞) c (∞) ＜ 0.9 t ＞ t s

式中，c(t)为信号函数t时刻的取值；c_max为信号函数中的最大值；

c (∞)

为信号函数在时间无限大时的取值，即信号的稳态值，当信号不稳定时该值不存在。

为量化分析各极点对系统小扰动振荡过程的影响，引入参与因子指标

p k i

，可表示为

（41）

p k i = | Re s[K (s), p i] |

式中，Res[K(s), p_i]为在p_i处的留数；p_k_i为留数的模值，反映极点对系统振荡的影响程度。

此外，引入衰减系数β，以量化系统的小扰动稳定性能，可表示为

（42）

β = ∑ p k i e 0.01 p i

3 算例分析

收起

3.1 单光伏并网系统分析

在Matlab/Simulink中搭建模型进行仿真验证，其中，单光伏并网变换器系统参数见表1，电流基准值为100 A。

采用扫频法，通过在变换器端口串联电压源注入不同频率扰动并测量其响应，以获取系统的幅频特性与相频特性，与理论分析结果对比如图8所示。上述结果表明，扫频结果与理论分析模型基本一致，验证了所提模型的准确性。

在去掉传递函数的零极点对消后，系统的极点分布如图9所示。各通道的极点相同且实部均为负数，表明系统当前位于稳定域中。设置i_c,_i_,q,ref =0，且0.4 s时刻发生扰动故障，0.5 s后故障切除，并网变换器的输出电压与电流如图10所示。

由图10可得，系统在发生小扰动故障时能够保持稳定运行，与理论分析结果相符。设置i_c_,i_,q,ref的变化范围为[-2.5,2.5](pu)，以靠近虚轴近极点轨迹为例进行分析，临界极点变化轨迹如图11所示。随着并网变换器输出无功功率的增加，系统的p₂极点向右半平面偏移。这是因为并网变换器系统输出的无功功率逐渐增加，导致系统静态安全裕度降低，从而使得系统由稳定运行状态转变为不稳定运行状态。同时，计算系统参与因子与衰减系数，其变化轨迹分别如图12和图13所示。

由图12和图13可得，当系统处于安全域内时，锁相环控制产生的p₃极点占据主导地位，并且通过合理配置锁相环控制参数能够有效提升系统的稳定性能。随着输出无功功率的增加，极点

p 2

逐渐接近虚轴，同时

p k2

迅速上升，从而使系统的稳定性能主要受到有功功率控制环节的限制。通过对图中结果进行分析可以发现，在系统接近临界稳定运行边界时，衰减系数β会快速上升。因此，β可以作为衡量系统稳定性能的一个重要指标之一，并且应避免其过大。

图14为光伏并网变换器输出电流为0.6(pu)时小扰动波形。对比图10可得，当光伏并网逆变器向电网输送无功电流时，系统的超调量σ%从7.13%增加为12.76%，调节时间从0.017 s增加为0.038 s，静态同步稳定性能下降。上述分析结果与图10相同，光伏逆变器向电网输送无功电流不利于系统的静态同步稳定性能的提升。图15给出了系统小扰动不稳定时输出的电压电流波形，系统发生静态同步失稳现象，此时系统存在实部为正的极点p₂=6.27。仿真结果与理论分析相符，验证了本文提出的静态同步稳定分析方法的正确性。

3.2 含STATCOM的光伏并网系统分析

在3.1节光伏并网系统模型的基础上增加STATCOM，其参数见表2。采用扫频法分析模型的准确度，具体结果如图16所示，验证了所提模型的准确性。并网变换器等效导纳矩阵传递函数伯德图对比结果如图17所示。在增加STATCOM后，主要影响并网逆变器q耦合通道，原因在于跟网型并网变换器的同步环节使用锁相环控制策略，同步稳定性能受到并网点q轴电压影响。当系统发生振荡时，由于各并网变换器所属的dq参考坐标系不相同，导致其q轴电压振荡方式不同，在i_q控制通道上会额外产生耦合作用效果。

系统dd通道传递函数不受影响，因此在部分文献中，将dd通道的导纳值作为跟网型并网变换器同步稳定判据之一^[10]。系统极点分布如图18所示，部分关键极点见表3。

对比图9和图18可得，在配置STATCOM后，系统中在dd通道引入了p₁、p₂和p₁₇极点，在q耦合通道引入了p₂₂、p₂₃、p₅₀、p₅₁、p₅₂和p₅₃极点，其分布更加接近临界稳定运行边界。若按照传统的稳定分析标准，系统的稳定性能应降低。在发生小扰动故障后，并网点处电压和电流如图19所示。

对比图14和图19可得，系统的超调量σ%从12.76%减少为5.76%，调节时间从0.038 s减少为0.010 s，系统暂态稳定性能反而得到提升。因此，在并网变换器控制参数配置合理的情况下，配置一定容量的无功补偿装置，能够有效提升系统的静态同步稳定性能。仅以极点的实部幅值作为稳定性能评价会致使分析结果偏保守。部分关键极点的参与因子分布见表4。结果表明，STATCOM引入新极点的同时也引入了附近的零点，自身参与因子小且降低了原有临界极点的参与因子，提高了系统的静态稳定性能。

系统的衰减系数β从3.15×10^-4变为1.27×10^-5，降低了95.97%，其原因为配置STATCOM后改变了系统内零点和极点的分布位置，零点与极点分布更加相近，降低极点在振荡中的影响。图20为闭环传递函数K(s)的伯德图，配置STATCOM能够消除原并网变换器dq通道与qq通道中谐振频率点，从而降低系统在低频至中频段的幅值，并且与衰减系数分析结果相一致。

3.3 控制参数影响分析

本节主要分析STATCOM外环控制参数对于并网变换器静态同步稳定性的影响，并以光照强度为2 000、1 000、500 W/m²下的光伏典型运行工况为例进行分析。

3.3.1 比例控制系数k_p_,s,_dc影响

设置直流电压控制环节比例系数k_p_,s,_dc变化范围为[0,7]，其余参数不变。在上述三种运行工况下，当k_p_,s,_dc增大时，q通道极点受到影响较小，几乎不发生变化，原因在于STATCOM并网变换器与光伏并网变换器的锁相环参数配置相近，在旋转坐标系下耦合影响进一步削弱。d通道主导极点变化如图21所示，在三种光伏典型运行工况下，随着k_p_,s,_dc的增大，d通道主导极点向远离虚轴移动，光照强度较低运行工况需要设置较大的k_p_,s,_dc才能够保证系统的静态同步稳定性。

参与因子分析过程与上述分析过程相同，结果过于庞大，此处给出分析结论，结合图21极点分布有如下结论：在k_p_,s,_dc较小时，系统稳定性由极点p₃、p₄、p₁₇主导影响；当k_p,_s_,dc较大时，系统稳定性主要由中频段极点主导影响。

系统的衰减系数β如图22所示。在稳定运行区域，随着k_p,_s_,dc的增大，系统的衰减系数不断增大，系统静态同步稳定的暂态性能下降，不利于系统的稳定运行；光照强度较高的运行工况设定较小的k_p,_s_,dc便能够满足静态同步稳定性能。由上述分析可得知，在实际工程中，当控制参数不可调整时，可通过调整光伏机组接受的光照强度大小实现静态同步稳定性能的调整。

3.3.2 积分控制系数k_in_,s,_dc影响

设置直流电压控制环节积分系数k_in_,s,_dc变化范围为[0,30]，其余参数不变，可得系统的衰减系数β如图23所示。在光照强度较高的光伏运行工况下，系统的缩减系数β随着积分系数k_in_,s,_dc的增加而增加，较小的积分系数k_in_,s,_dc就能够实现较好的稳定性能。在其他光照强度下，系统的缩减系数β随着积分系数k_in_,s,_dc的增加先减小后增加，其变化轨迹上存在一个极值点。

结合图22直流电压控制环节比例系数k_p_,s,_dc衰减系数变化轨迹，可得：在稳定运行区域内，逐渐减小外环控制的带宽有助于提高系统的静态同步稳定性能；而当外环控制的带宽持续减小至阈值时，关键控制信号可能丢失，从而不利于系统的静态同步稳定性能提升。因此，在配置系统各控制环节的带宽值时需要合理考虑。

本文参数分析验证得到的主要结论如下：

1）在稳定运行区域内部，直流电压控制环节、有功电流控制环节及无功电流控制环节的带宽增大会导致系统的衰减系数增大，系统振荡幅值增大，衰减速度降低。

2）当STATCOM装置的锁相环参数与光伏逆变系统相同时，系统的稳定性能主要受到电网阻抗和各并网逆变器等效导纳传递函数影响，不存在d轴对q轴的耦合作用影响。此时只需满足各并网逆变器均可独立接入电网稳定运行即可保证系统静态同步稳定。

3）在实际工程中，当控制参数不可调整时，可通过调整光伏机组接受的光照强度大小实现静态同步稳定性能的调整。

4）根据本文仿真分析得到的结果，总结规律如下：当光伏系统的衰减系数大于300时，可判定系统出现静态同步失稳现象；当衰减系数指标小于250时，则表明系统处于静态同步稳定状态；当衰减系数指标介于二者间需进一步采用稳定判据进行判别。

3.4 STATCOM占比影响分析

由于STATCOM无功容量的输出是通过改变i_c_,s_,q,ref的数值实现的，所以将i_c_,s_,q,ref的变化范围设置为[-2.5 pu, 2.5 pu]，以分析不同STATCOM容量占比的影响，如图24所示。对比图13可得，在光照强度为1 000 W/m²的光伏典型运行工况下，由STATCOM向电网提供无功电流相较于光伏并网变换器提供无功电流的静态同步稳定性能更好。

由图24可得，当光照强度增加时，适当提高STATCOM容量以向光伏系统提供容性无功支撑，能够有效地提升系统的静态同步稳定性能，原因在于STATCOM向光伏系统提供容性无功支撑可以有效抑制光伏系统有功出力增加导致的电压下降，并为并网点电压稳定性提供保障。在光照强度较低的情况下，过高的无功补偿量会降低系统的静态同步稳定性能，反而不利于系统的稳定运行。因此，在不同的光伏运行工况下需要及时调整STATCOM的无功补偿量。

4 结论

收起

针对含有STATCOM的光伏并网系统，建立了静态同步稳定分析模型，改进了稳定判据和静态同步稳定评估指标，详细分析了增加STATCOM对光伏并网系统静态同步稳定性影响，主要结论如下：

1）在已知并网变换器控制参数的情况下，合理地优化并网变换器输出的无功功率可以有效提升其静态同步稳定性能；小扰动后振荡模式主要由锁相环控制产生的极点所主导。

2）在并网变换器运行工况接近临界稳定区域时，系统的衰减系数β会迅速增加。因此，在运行过程中，系统需要具备一定的衰减系数裕度。根据本文仿真分析结果表明，当光伏系统的衰减系数大于300时，可以判定系统出现静态同步失稳现象；而当衰减系数小于250时，则表示系统处于静态同步稳定状态。这为系统的静态同步稳定性分析提供了判别标准。

3）光伏并网系统在增加STATCOM后，主要影响了q耦合通道输出的导纳矩阵传递函数，而dd通道则不会受到并网变换器间锁相环耦合振荡的影响。在稳定运行区域内部，直流电压控制环节、有功电流控制环节以及无功电流控制环节的带宽增大，将导致系统衰减系数增加、系统振荡幅值增大，并且降低系统暂态性能中的衰减速度。

4）当STATCOM的锁相环参数与光伏逆变系统相同时，系统的稳定性能主要受到电网阻抗和各并网逆变器等效导纳传递函数影响。在这种情况下，只需确保各并网逆变器能够独立接入电网并实现稳定运行，以保证系统达到静态同步稳定。

5）当光伏系统有功出力较小时，STATCOM倾向于提供较小的容性无功出力或者提供较大的感性无功出力，以提高系统的静态同步稳定性能。然而，在光伏系统有功出力较大时，采用STATCOM来提供容性无功补偿可以更优于光伏并网变换器自身所提供的无功补偿。因此，合理配置STATCOM容量，能够有效地增强系统的静态同步稳定性能。

本文主要针对含STATCOM的光伏并网系统静态同步稳定分析进行研究，下一步将研究更加精细化控制参数的设计方法。

基金

收起

国家自然科学基金智能电网联合基金集成项目(U22B6008)
山东省自然科学基金项目(ZR2021QE133)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

韩应生, 孙海顺, 秦世耀, 等. 电压源型双馈风电并网系统小扰动低频稳定性分析[J]. 电工技术学报, 2023, 38(5): 1312-1324, 1374.

Han

Yingsheng

, Sun

Haishun

, Qin

Shiyao

, et al. Low-frequency stability analysis of voltage-sourced doubly-fed wind power grid-connected system under small disturbance[J]. Transactions of China Electrotechnical Society, 2023, 38(5): 1312-1324, 1374.

[2]

汪春江, 孙建军, 宫金武, 等. 并网逆变器与电网阻抗交互失稳机理及阻尼策略[J]. 电工技术学报, 2020, 35(增刊2): 503-511.

Wang

Chunjiang

, Sun

Jianjun

, Gong

Jinwu

, et al. Mechanism and damping strategy of interactive instability between grid-connected inverter and grid impedance[J]. Transactions of China Electrotechnical Society, 2020, 35(S2): 503-511.

[3]

曹永吉, 吴秋伟, 张恒旭, 等. 考虑分时电价和最大暂态频率偏移的储能机会约束优化配置[J]. 电力系统自动化, 2023, 47(2): 61-68.

Cao

Yongji

, Wu

Qiuwei

, Zhang

Hengxu

, et al. Chance-constrained optimal configuration for energy storage system considering time-of-use price and maximum transient frequency deviation[J]. Automation of Electric Power Systems, 2023, 47(2): 61-68.

[4]

马睿聪, 邹世豪, 曹永吉, 等. 基于灵敏度分析的同步调相机及其励磁系统参数分类辨识方法[J]. 电力系统保护与控制, 2024, 52(4): 26-37.

Ruicong

, Zou

Shihao

, Cao

Yongji

, et al. Sensitivity analysis-based classification and identification method of a synchronous condenser and its excitation system parameter[J]. Power System Protection and Control, 2024, 52(4): 26-37.

[5]

罗魁, 郭剑波, 王伟胜, 等. 跟网型新能源附加频率控制对频率稳定及小扰动同步稳定影响分析综述[J]. 中国电机工程学报, 2023, 43(4): 1261-1280.

Luo

Kui

, Guo

Jianbo

, Wang

Weisheng

, et al. Review of impact of grid following variable renewable energy supplementary frequency control on frequency stability and small-disturbance synchronization stability[J]. Proceedings of the CSEE, 2023, 43(4): 1261-1280.

[6]

San

Guocheng

, Zhang

Wenlin

, Luo

Ran

, et al. Small-signal multi-frequency model for grid-connected inverter system with PWM effect[J]. CSEE Journal of Power and Energy Systems, 2020, 6(2): 307-317.

[7]

Jiale

, Qi

Yang

, Deng

Han

, et al. Evaluating small-signal synchronization stability of grid-forming converter: a geometrical approach[J]. IEEE Transactions on Industrial Electronics, 2022, 69(9): 9087-9098.

[8]

胡祺, 付立军, 马凡, 等. 弱电网下基于锁相控制并网变换器小扰动同步稳定分析[J]. 中国电机工程学报, 2021, 41(1): 98-108.

, Fu

Lijun

, Ma

Fan

, et al. Small signal synchronizing stability analysis of PLL-based VSC connected to weak AC grid[J]. Proceedings of the CSEE, 2021, 41(1): 98-108.

[9]

Wen

, Dong

Dong

, Boroyevich

, et al. Impedance-based analysis of grid-synchronization stability for three-phase paralleled converters[J]. IEEE Transactions on Power Electronics, 2016, 31(1): 26-38.

[10]

Mansour

M Z

, Mohammed

, Ravanji

M H

, et al. Output impedance derivation and small-signal stability analysis of a power-synchronized grid following inverter[J]. IEEE Transactions on Energy Conversion, 2022, 37(4): 2696-2707.

[11]

郝艺, 周瑀涵, 刘晨曦, 等. 含跟网型储能的新能源多馈入系统小扰动电压支撑强度分析[J]. 电工技术学报, 2024, 39(11): 3569-3580.

Hao

, Zhou

Yuhan

, Liu

Chenxi

, et al. Small-disturbance voltage support strength analysis for renewable multi-infeed system with grid-following energy storage[J]. Transactions of China Electro-technical Society, 2024, 39(11): 3569-3580.

[12]

宗皓翔, 张琛, 吕敬, 等. MMC多入多出阻抗及其在不对称小扰动稳定分析中的应用[J]. 中国电机工程学报, 2022, 42(15): 5649-5664.

Zong

Haoxiang

, Zhang

Chen

, Lü

Jing

, et al. MMC MIMO impedance and its application in the asymmetric small-signal stability analysis[J]. Proceedings of the CSEE, 2022, 42(15): 5649-5664.

[13]

艾力西尔·亚尔买买提, 辛焕海, 宫泽旭, 等. 锁相环主导下变流器振荡问题的阻抗分析方法:标称性与鲁棒性[J]. 中国电机工程学报, 2024, 44(15): 6035-6047.

Yaermaimaiti

Ailixier

, Xin

Huanhai

, Gong

Zexu

, et al. Impedance analysis method of oscillation issues in converter caused by phase-locked loop: nominal performance and robustness[J]. Proceedings of the CSEE, 2024, 44(15): 6035-6047.

[14]

Zhang

Chen

, Føyen

, Suul

J A

, et al. Modeling and analysis of SOGI-PLL/FLL-based synchronization units: stability impacts of different frequency-feedback paths[J]. IEEE Transactions on Energy Conversion, 2021, 36(3): 2047-2058.

[15]

宫泽旭, 艾力西尔·亚尔买买提, 辛焕海, 等. 新能源电力系统并网设备小扰动稳定分析(一): 机理模型与稳定判据适用性[J]. 中国电机工程学报, 2022, 42(12): 4405-4418.

Gong

Zexu

, Yaermaimaiti

Ailixier

, Xin

Huanhai

, et al. Small signal stability analysis of equipment in renewable energy power system (part Ⅰ): mechanism model and adaptation of stability criterion[J]. Proceedings of the CSEE, 2022, 42(12): 4405-4418.

[16]

宫泽旭, 艾力西尔·亚尔买买提, 辛焕海, 等. 新能源电力系统并网设备小扰动稳定分析(二): 导出机理与稳定性分类探讨[J]. 中国电机工程学报, 2022, 42(14): 5126-5139.

Gong

Zexu

, Yaermaimaiti

Ailixier

, Xin

Huanhai

, et al.. Small signal stability analysis of grid-connected equipment in power system (part Ⅱ): discussion on mechanism derivation and classification of stability[J]. Proceedings of the CSEE, 2022, 42(14): 5126-5139.

[17]

康勇, 林新春, 郑云, 等. 新能源并网变换器单机无穷大系统的静态稳定极限及静态稳定工作区[J]. 中国电机工程学报, 2020, 40(14): 4506-4515.

Kang

Yong

, Lin

Xinchun

, Zheng

Yun

, et al. The static stable-limit and static stable-working zone for single-machine infinite-bus system of renewable-energy grid-connected converter[J]. Proceedings of the CSEE, 2020, 40(14): 4506-4515.

[18]

曹永吉, 张恒旭, 施啸寒, 等. 规模化分布式能源参与大电网安全稳定控制的机制初探[J]. 电力系统自动化, 2021, 45(18): 1-8.

Cao

Yongji

, Zhang

Hengxu

, Shi

Xiaohan

, et al. Preliminary study on participation mechanism of large-scale distributed energy resource in security and stability control of large power grid[J]. Automation of Electric Power Systems, 2021, 45(18): 1-8.

[19]

王凯, 原熙博, 王晖, 等. 弱电网下逆变器并网系统失稳的两种诱因[J]. 电网技术, 2023, 47(3): 1230-1241.

Wang

Kai

, Yuan

Xibo

, Wang

Hui

, et al. Two causes of the inverter-based grid-connected system instability in weak grids[J]. Power System Technology, 2023, 47(3): 1230-1241.

[20]

詹长江, 吴恒, 王雄飞, 等. 构网型变流器稳定性研究综述[J]. 中国电机工程学报, 2023, 43(6): 2339-2358.

Zhan

Changjiang

, Wu

Heng

, Wang

Xiongfei

, et al. An overview of stability studies of grid-forming voltage source converters[J]. Proceedings of the CSEE, 2023, 43(6): 2339-2358.

[21]

王冠中, 董炜, 辛焕海, 等. 基于广义短路比的电力电子多馈入系统小干扰概率稳定评估[J]. 电力系统自动化, 2018, 42(18): 17-24.

Wang

Guanzhong

, Dong

Wei

, Xin

Huanhai

, et al. Probabilistic stability evaluation of small disturbances for power electronic multi-infeed system based on generalized short-circuit ratio[J]. Automation of Electric Power Systems, 2018, 42(18): 17-24.

[22]

Kai

, Gan

Chun

, Hu

Yihua

, et al. Input-output small-signal stability analysis of a PLL-free direct power-controlled partially power-decoupled more-electric shipboard propulsion system[J]. IEEE Transactions on Transportation Electrification, 2021, 7(3): 1672-1686.

[23]

Chen

Yahong

, Xu

Wei

, Liu

, et al. Small-signal system frequency stability analysis of the power grid integrated with type-II doubly-fed variable speed pumped storage[J]. IEEE Transactions on Energy Conversion, 2023, 38(1): 611-623.

[24]

曾德银, 姚骏, 张田, 等. 虚拟同步发电机多机并联系统的频率小信号稳定性分析研究[J]. 中国电机工程学报, 2020, 40(7): 2048-2061.

Zeng

Deyin

, Yao

Jun

, Zhang

Tian

, et al. Research on frequency small-signal stability analysis of multi-parallel virtual synchronous generator-based system[J]. Proceedings of the CSEE, 2020, 40(7): 2048-2061.

[25]

Mansour

M Z

, Ravanji

M H

, Karimi

, et al. Small-signal synchronization stability enhancement of grid-following inverters via a feedback linearization controller[J]. IEEE Transactions on Power Delivery, 2022, 37(5): 4335-4344.

[26]

李红, 梁军杨, 王振民, 等. 跟网型变换器的小扰动同步稳定机理分析与致稳控制[J]. 电工技术学报, 2024, 39(12): 3802-3815.

Hong

, Liang

Junyang

, Wang

Zhenmin

, et al. Small signal synchronization stability analysis and improved control strategy for grid following converter[J]. Transactions of China Electrotechnical Society, 2024, 39(12): 3802-3815.

[27]

Sajadi

, Zhao

Shuang

, Clark

, et al. Small-signal stability analysis of large-scale power systems in response to variability of offshore wind power plants[J]. IEEE Systems Journal, 2019, 13(3): 3070-3079.

[28]

马富艺龙, 辛焕海, 刘晨曦, 等. 新能源基地柔性直流送出系统小扰动电压支撑强度评估[J]. 电工技术学报, 2023, 38(21): 5758-5770, 5938.

Fuyilong

, Xin

Huanhai

, Liu

Chenxi

, et al. Small-disturbance system voltage support strength assessment method for renewables VSC-HVDC delivery system[J]. Transactions of China Electrotechnical Society, 2023, 38(21): 5758-5770, 5938.

[29]

张宇, 张琛, 蔡旭, 等. 并网变换器的暂态同步稳定性分析: 稳定域估计与镇定控制[J]. 中国电机工程学报, 2022, 42(21): 7871-7883.

Zhang

, Zhang

Chen

, Cai

, et al. Transient grid-synchronization stability analysis of grid-tied voltage source converters: stability region estimation and stabilization control[J]. Proceedings of the CSEE, 2022, 42(21): 7871-7883.

2025年第40卷第11期

PDF下载

348

142

引用本文

BibTeX

文章信息

doi: 10.19595/j.cnki.1000-6753.tces.240754

接收时间：2024-05-10
首发时间：2025-11-06
出版时间：2025-06-10

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2024-05-10
修回日期：2024-06-25

基金

国家自然科学基金智能电网联合基金集成项目(U22B6008)

山东省自然科学基金项目(ZR2021QE133)

作者信息

¹ 山东大学电气工程学院济南 250061

² 山东大学智能创新研究院济南 250101

通讯作者:

曹永吉男，1992年生，副研究员，硕士生导师，研究方向为电力系统稳定分析与控制、可再生能源并网及储能技术应用。E-mail：yongji@sdu.edu.cn

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/dgjsxb/CN/10.19595/j.cnki.1000-6753.tces.240754

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT

参数	数值	参数	数值
直流参考电压U_dc^ref/V	800	滤波电感1/mH	0.5
无功电流参考值i_c_,i_,q,ref/A	0	滤波电感2/mH	0.3
开关频率/kHz	10	滤波电容/μF	10
交流电网电压u_g/V	415	采样时间/μs	1
比例系数k_p,_i_,dc	0.5	积分系数k_in,_i_,dc	25
比例系数k_p,_i_,c	3	积分系数k_in,_i_,c	40
比例系数k_p,_i_,pll	3	积分系数k_in,_i_,pll	24
电网电阻/μΩ	1	系统容量/(kV·A)	1 000
电网电感/mH	1.5	光照强度/(W/m²)	1 000

参数

数值

参数

数值

直流参考电压U_dc^ref/V

800

滤波电感1/mH

0.5

无功电流参考值i_c_,i_,q,ref/A

滤波电感2/mH

0.3

开关频率/kHz

滤波电容/μF

交流电网电压u_g/V

415

采样时间/μs

比例系数k_p,_i_,dc

0.5

积分系数k_in,_i_,dc

比例系数k_p,_i_,c

积分系数k_in,_i_,c

比例系数k_p,_i_,pll

积分系数k_in,_i_,pll

电网电阻/μΩ

系统容量/(kV·A)

1 000

电网电感/mH

1.5

光照强度/(W/m²)

1 000

参数	数值	参数	数值
直流参考电压 $U dc ref /V$	800	滤波电感1/mH	0.5
无功参考电流i_c_,s_,q,ref/A	50	滤波电感2/mH	0.3
开关频率/kHz	10	滤波电容/μF	10
交流电网线电压u_g/V	415	采样时间/μs	1
比例系数k_p,_s_,dc	5	积分系数k_in,_s_,dc	20
比例系数k_p,_s_,c	0.3	积分系数k_in,_s_,c	40
比例系数k_p,_s_,pll	3	积分系数k_in,_s_,pll	40
电网电阻/μΩ	1	系统容量/(kV·A)	1 000
电网电感/mH	1.5	光照强度/(W/m²)	1 000

参数

数值

参数

数值

直流参考电压

U dc ref /V

800

滤波电感1/mH

0.5

无功参考电流i_c_,s_,q,ref/A

滤波电感2/mH

0.3

开关频率/kHz

滤波电容/μF

交流电网线电压u_g/V

415

采样时间/μs

比例系数k_p,_s_,dc

积分系数k_in,_s_,dc

比例系数k_p,_s_,c

0.3

积分系数k_in,_s_,c

比例系数k_p,_s_,pll

积分系数k_in,_s_,pll

电网电阻/μΩ

系统容量/(kV·A)

1 000

电网电感/mH

1.5

光照强度/(W/m²)

1 000

d通道	q通道
p₁	-8.032+0.003i	p₂₂	-8.032+0.006i
p₂₃	-8.032-0.006i
p₂	-8.032-0.003i	p₅₀	-8.038
p₅₁	-8.025
p₂₁	-3457.421	p₅₂	-4.812
p₅₃	-4.811

d通道

q通道

p₁

-8.032+0.003i

p₂₂

-8.032+0.006i

p₂₃

-8.032-0.006i

p₂

-8.032-0.003i

p₅₀

-8.038

p₅₁

-8.025

p₂₁

-3457.421

p₅₂

-4.812

p₅₃

-4.811

通道	参与因子
d	p_k1	1.28×10^-12
p_k2	1.28×10^-12
p_k17	3.32×10^-13
q	p_k22	6.25×10^-12
p_k23	6.25×10^-12
p_k29	9.05×10^-12
p_k30	9.05×10^-12
p_k31	3.32×10^-13
p_k32	3.35×10^-13

通道

参与因子

p_k1

1.28×10^-12

p_k2

1.28×10^-12

p_k17

3.32×10^-13

p_k22

6.25×10^-12

p_k23

6.25×10^-12

p_k29

9.05×10^-12

p_k30

9.05×10^-12

p_k31

3.32×10^-13

p_k32

3.35×10^-13