直驱风机接入弱电网宽频带振荡机理与抑制方法(一)：阻抗建模

原创已于 2025-06-20 07:12:02 修改 · 848 阅读

18 ·

CC 4.0 BY-SA版权

文章标签：

#matlab #阻抗建模 #永磁直驱风机 #谐波线性化

于 2025-06-19 16:31:17 首次发布

永磁直驱风机宽频振荡建模和分析专栏收录该内容

1 篇文章

订阅专栏

一、研究背景

二新能源机组的频率耦合阻抗/导纳建模方法

三直驱风机网侧变换器序阻抗建模

3.1 永磁直驱风力发电装备的拓扑结构与控制方法

3.2 永磁直驱风力发电装备的序阻抗建模

四、阻抗特性曲线绘制

五、参考文献

一、研究背景

受风、光资源的制约，新能源发电基地大多分布在草原、山区、荒漠等偏远地区，长距离的输电线路使得这些地区的电网较为薄弱，新能源发电并网变换器往往运行在弱交流电网的条件下。在并网变换器发展的初期，由于并网变换器的容量小，数量少。从变换器看，电网就是一个无穷大电源，传统并网变换器的控制与参数设计方法是将电网阻抗看作为零来进行设计。因此，当新能源并网变换器渗透率提高，而交流电网的特性变弱时，电网已不能再看作是理想电网，新能源并网变换器接入弱电网容易引发宽频带振荡问题。2009年美国德州发生了第一起公开报道的新能源并网次同步振荡事故；2012 年以来，我国张北地区也发生了多起双馈风电机组经串补并网次同步振荡事故；2015 年 7 月，我国新疆哈密地区发生了永磁直驱机组并入弱交流电网次同步振荡事故，造成了哈密地区 3 台火电机组轴系扭振切机；多处新能源发电场站通过长距离输电线路接入电网引发了高频振荡问题。这些宽频带振荡问题已经成为影响我国电力系统安全稳定运行的重要威胁之一，研究其产生机理及新型抑制技术，以解决我国新能源电力系统中现实存在的宽频带振荡新问题，对新能源的消纳和电力系统的安全稳定运行具有重要的理论价值和工程意义。

二新能源机组的频率耦合阻抗/导纳建模方法

这部分可以参考我在csdn上面的另一篇文章：

VSG宽频振荡建模与机理分析(一)：考虑频率耦合效应和电压电流双环的VSG阻抗建模_频率耦合阻抗模型-CSDN博客文章浏览阅读1.2k次，点赞29次，收藏28次。摘要：本文研究了虚拟同步发电机(VSG)的阻抗建模方法，重点分析考虑频率耦合效应的序阻抗建模。针对新能源发电装备中逆变器的非线性特性，比较了dq阻抗建模和序阻抗建模两种方法，指出后者具有物理意义清晰、阻抗形式简单等优势。通过建立VSG的拓扑结构和控制模型，推导考虑频率耦合效应的频域导纳矩阵表达式，并绘制阻抗特性曲线。研究表明，VSG系统中的P/Q控制和dq控制会产生显著的镜频耦合效应，需采用2×2导纳矩阵准确描述其阻抗特性。该建模方法为分析VSG并网系统的稳定性问题提供了理论基础。_频率耦合阻抗模型https://blog.csdn.net/zqjx0113/article/details/148192106?spm=1011.2415.3001.10575&sharefrom=mp_manage_link

三直驱风机网侧变换器序阻抗建模

3.1 永磁直驱风力发电装备的拓扑结构与控制方法

永磁直驱风力发电装备的拓扑结构图如图 1 所示，通过并网逆变器将直流电转换为交流电并入电网。图 1 中：ea、eb 和 ec 为并网逆变器桥臂中点电压；va、vb 和 vc 为 PCC 点电压；ia、ib 和 ic 为并网逆变器电感电流；Lf、Cf和 Rf分别为并网逆变器的滤波电感、滤波电容和阻尼电阻；Lg 和 Rg 为电网的等效线路电感和电阻；vga、vgb 与 vgc 和 iga、igb 与 igc 分别为电网电压和入网电流。

永磁直驱风力发电装备并网逆变器的控制方法如图 2 所示，采用电压、电流双闭环控制和同步旋转坐标系下锁相环（phase locked loop, PLL）控制。图 2中：vdcr 为并网逆变器直流侧电压指令；vdc 为并网逆变器直流侧电压；Idr、Iqr 分别为同步旋转坐标系下 d 轴、q 轴输出电流指令；Kd 为前馈解耦增益系数；Kf为电压前馈系数；cd 和 cq 为同步旋转坐标系下的调制波；ca、cb 和 cc 为静止坐标系下的调制波；kp_I 和 ki_I 分别是电流 PI 控制的比例与积分系数；kp_PLL 和 ki_PLL 分
别是锁相环 PI 控制的比例与积分系数；Q1~Q6 为并网逆变器中各个 IGBT 的驱动信号。

图1 永磁直驱风力发电装备拓扑结构

图2 永磁直驱并网逆变器控制框图

3.2 永磁直驱风力发电装备的序阻抗建模

对永磁直驱风力/光伏发电装备建模时，由于风机 PWM 整流器或者光伏板与并网逆变器之间经过容值较大的直流电容相连接，在稳定运行时，直流侧的电压波动很小，可将直流侧看作是稳定值，故并网逆变器的电压控制外环在阻抗建模时可以忽略不计，可降低建模难度。因此，并网逆变器的电流控制指令 Idr、Iqr可看作为给定常量。

假设在并网逆变器与电网连接的 PCC 点处注入正、负序扰动电压，电压幅值的大小不超过工频电压的 10%，即小于系统工频电压一个数量级，因此不会对系统造成太大的波动而改变系统的稳定工作点。此时，A 相的电压和电流为：

$v_a(t) = V_1 \cos(2\pi f_1 t) + V_p \cos(2\pi f_p t + \varphi_{vp}) + V_n \cos(2\pi f_n t + \varphi_{vn}) (1)$

$i_a(t) = I_1 \cos(2\pi f_1 t + \varphi_{i1}) + I_p \cos(2\pi f_p t + \varphi_{ip}) + I_n \cos(2\pi f_n t + \varphi_{in}) (2)$

式中，V1、Vp 和 Vn 分别为基波电压、正序扰动电压和负序扰动电压的幅值；I1、Ip 和 In 分别为基波电流、正序扰动电流响应和负序扰动电流响应的幅值；f1、fp和 fn 分别为基波频率、正序扰动频率和负序扰动频率；φvp 和 φvn 分别为正序扰动电压和负序扰动电压的初相角；φi1、φip 和 φin 分别为基波电流、正序扰动电流响应和负序扰动电流响应的初相角。在频域中，电压 va 和电流 ia 可写成：

$\mathbf{V}_a[f] = \begin{cases} \mathbf{V}_1, & f = \pm f_1 \\ \mathbf{V}_p, & f = \pm f_p \\ \mathbf{V}_n, & f = \pm f_n \end{cases} \quad \mathbf{I}_a[f] = \begin{cases} \mathbf{I}_1, & f = \pm f_1 \\ \mathbf{I}_p, & f = \pm f_p \\ \mathbf{I}_n, & f = \pm f_n \end{cases}(3)$

式中， ${\text{V}_1}=V_1/2$ ； ${\text{V}_{\text{p}}}=(V_{\text{p}}/2)e^{\pm{j\varphi_{\text{vp}} }}$ ； ${\text{V}_{\text{n}}}=(V_{\text{n}}/2)e^{\pm{j\varphi_{\text{vn}} }}$ ； ${\text{I}_{\text{1}}}=(I_{\text{1}}/2)e^{\pm{j\varphi_{\text{i1}} }}$ ； ${\text{I}_{\text{p}}}=(I_{\text{p}}/2)e^{\pm{j\varphi_{\text{ip}} }}$ ； ${\text{I}_{\text{n}}}=(I_{\text{n}}/2)e^{\pm{j\varphi_{\text{in}} }}$ 。

由图1可得到并网逆变器的电感电流、桥臂中点电压和 PCC 点电压的关系式为：

$sL_f \begin{bmatrix} i_a \\ i_b \\ i_c \end{bmatrix} = \begin{bmatrix} m_a \\ m_b \\ m_c \end{bmatrix} \frac{V_{dc}}{2} - \begin{bmatrix} v_a \\ v_b \\ v_c \end{bmatrix} (4)$

随后，分别对锁相环、电流内环等控制电路进行小信号建模，推导出ma在频域中关于ia和va的关系式（平衡工况以单相为例即可），即可推导出相应的正负序阻抗，详细推导过程见如下论文，不再详述：

新能源发电接入弱电网的宽频带振荡机理及抑制方法研究 - 中国知网 (cnki.net)https://kns.cnki.net/kcms2/article/abstract?v=i7m15r_oBqqWLoQOSWqq5cJeS8zjpvaQIf0OoY43MuQZ_iYE4iz4S0HY21_RnRTRc9rRy6FBIpdh9Kx8Yz2_OK8_lLK6nJ93frgN39AvTxe52T4I0GDi-pQOgvlR5mMgBpmZ2wkBZglyE0ZJPPdsQY3rfuApFQ4sdyfO_uU1IJnIkkzzJLMIEht4S36hjg-0&uniplatform=NZKPT&language=CHS永磁直驱/光伏并网逆变器的正、负序阻抗表达式分别为：

$Z_p(s) = -\frac{V_p(s)}{I_p(s)} = \frac{ sL_f + \frac{V_{dc}}{2} \left( H_i(s - j2\pi f_1) - jK_d \right) G_i(s) }{ 1 - \frac{V_{dc}}{2} K_f G_v(s) - \frac{V_{dc}}{4} \frac{T_{PLL}(s - j2\pi f_1)}{V_1} G_v(s) H_i(s - j2\pi f_1)(I_{dr} + jI_{qr}) } (5)$

$Z_n(s) = -\frac{V_n(s)}{I_n(s)} = \frac{ sL_f + \frac{V_{dc}}{2} \left( H_i(s + j2\pi f_1) + jK_d \right) G_i(s) }{ 1 - \frac{V_{dc}}{2} K_f G_v(s) - \frac{V_{dc}}{4} \frac{T_{PLL}(s + j2\pi f_1)}{V_1} G_v(s) H_i(s + j2\pi f_1)(I_{dr} - jI_{qr}) }(6)$

这里Zp(s)=-Vp(s)/Ip(s)是因为直驱风机并网电流是从逆变器流向电网的，而我们要计算逆变器的阻抗，电压除以电流，电流是反方向，所以加个负号。

四、阻抗特性曲线绘制

在MATLAB中编程，绘制690V永磁直驱风机的理论阻抗bode图，表1为永磁直驱风机的参数

表1 永磁直驱风机的参数

参数	数值
kp_PLL	0.2659
ki_PLL	10.9988
kp	0.01715
ki	22.85715
Lg	0.1mH
Cf	20e-6 F
Rf	1.5Ω
Vg	690V
f1	50Hz

在MATLAB中，通过扫频法绘制690V永磁直驱风机的理论阻抗bode图，即s=1i*2*pi*fp，fp为频率，通过给fp赋值，一般范围为1Hz~10000Hz，如下面代码所示：

f_begin=1;
f_end=10000;
Point=100000;
fsample= (logspace(log10(f_begin),log10(f_end),Point));

h=1;
for n=1:1:Point
    fp=fsample(n);
    w=2*pi*fp;
    s=w*1i;


    ……


end

690V永磁直驱风机理论正负序阻抗伯德图如下：