4. 并网逆变器小信号阻抗建模

小信号阻抗建模的核心目的,是把并网逆变器在某一稳态工作点附近“线性化”为一个等效的诺顿/戴维南电路,从而用阻抗比判据分析弱电网下的稳定性。与直接列写状态空间方程相比,阻抗模型更贴近工程测量:它把逆变器看作一个受控电流源并联输出导纳,把电网看作电压源串联线路阻抗,二者在公共耦合点(PCC)交互。本节围绕“坐标系选择—开环导纳—闭环导纳—关键传递函数”这条主线展开,并给出可复现的推导与示例参数。

一、坐标系约定与建模前提

建模前必须先固定参照系,否则正负号与交叉耦合项会前后矛盾。本节统一采用以下约定:

  • 三相电压、电流采用幅值不变(amplitude-invariant)的 Clarke/Park 变换,即 v_d = (2/3)(v_a cosθ + v_b cos(θ-120°) + v_c cos(θ+120°)),αβ 与 dq 之间为旋转矩阵关系。
  • dq 坐标系由 PCC 电压矢量定向,稳态时 V_q = 0,V_d = V_m(相电压幅值)。
  • 小信号量用前缀 Δ 表示,稳态量用大写字母表示,如 I_d、D_d。
  • 电流正方向定义为逆变器流向电网。

选择 dq 还是 αβ,取决于控制结构与关注频段。dq 坐标系下,基波变为直流量,PI 调节器可对直流量实现零稳态误差,因此适合分析低频段(数 Hz 至数百 Hz)的稳定性;但 dq 阻抗是 2×2 矩阵,存在 dq 交叉耦合。αβ 坐标系下,基波仍是 50 Hz 交流量,PR 谐振调节器可直接对交流量增益,阻抗矩阵在静止坐标系下为 2×2 但通常近似解耦,适合分析中高频段(数百 Hz 至数 kHz)及谐波谐振问题。两者可通过旋转矩阵相互转换,工程上常按“低频用 dq、高频用 αβ”分工。

并网逆变器小信号阻抗建模:从开环导纳到闭环输出导纳 物理对象 直流源 Vdc 三相桥 + LC/LCL 滤波 PCC 电压 Vd, Vq 输出电流 Id, Iq 电网阻抗 Zg 稳态工作点 (Vd, Id, Dd) 控制与线性化 电流环 PI/PR + 前馈解耦 PLL 小信号模型 Gpll(s) 调制与采样延迟 e^(-1.5Ts·s) 稳态点线性化 → 小信号 Δd = f(Δx, Δv) 坐标变换矩阵 T(θ) 输出导纳模型 开环导纳 Yo(s) 闭环导纳 Yc(s) = Yo/(1+T) dq 域 2×2 矩阵 αβ 域近似对角化 Zinv = Yc⁻¹ 与 Zg 比较判稳定 闭环输出导纳推导链(dq 域) ① 状态方程 L di/dt = v - Ri - e ② 稳态点 Vd, Id, Dd 恒定 ③ 小信号扰动 Δd, Δv, Δi ④ 消去 Δd Δd = Gc·Δi + Gff·Δv ⑤ 得 Yo(s) Δi = Yo·Δv 闭环导纳:Yc(s) = Yo(s) / [1 + T(s)],其中 T(s) = Gc(s)·Gd(s)·Yo(s) 为回路增益 稳定性判据:Zg(s)/Zinv(s) 满足 Nyquist 判据,弱电网下 Zg 增大易引发谐振

二、开环输出导纳 Yo(s) 的建立

以 L 型滤波并网逆变器为例,忽略寄生电阻,在 dq 坐标系下三相电感方程可写为:

L·d(Id)/dt = Vd - Ed + ωL·Iq
L·d(Iq)/dt = Vq - Eq - ωL·Id

其中 Ed、Eq 为桥臂输出电压的 dq 分量,ω 为基波角频率。加入小信号扰动并消去稳态量,得到线性化方程:

L·s·ΔId = ΔVd - ΔEd + ωL·ΔIq
L·s·ΔIq = ΔVq - ΔEq - ωL·ΔId

桥臂电压由调制波决定:Ed = (Vdc/2)·Dd,Eq = (Vdc/2)·Dq。将 ΔEd、ΔEq 代入并整理,可得到开环导纳矩阵:

Yo(s) = [ sL,  -ωL ;  ωL,  sL ]⁻¹
      = 1/(s²L² + ω²L²) · [ sL,  ωL ; -ωL, sL ]

该矩阵的非对角项 ±ωL 正是 dq 交叉耦合的来源。若在 αβ 坐标系下重写,电感方程变为 L·s·Δiα = Δvα - Δeα、L·s·Δiβ = Δvβ - Δeβ,导纳矩阵退化为对角形式 Yo_αβ = diag(1/(sL), 1/(sL)),这也是 αβ 域分析更直观的原因。

三、闭环输出导纳 Yc(s) 的推导

闭环导纳的关键在于把控制环“折合”进导纳。以 dq 域电流环为例,控制律为:

Dd = Gc(s)·(Id* - Id) + Vd/(Vdc/2) - ωL·Iq/(Vdc/2)
Dq = Gc(s)·(Iq* - Iq) + Vq/(Vdc/2) + ωL·Id/(Vdc/2)

其中 Gc(s) 为 PI 调节器,Gc(s) = Kp + Ki/s;后两项为前馈解耦项。对控制律取小信号,并考虑采样与调制延迟 Gd(s) = e^(-1.5Ts·s),可得:

ΔDd = -Gc(s)·Gd(s)·ΔId + ΔVd/(Vdc/2) - ωL·ΔIq/(Vdc/2)
ΔDq = -Gc(s)·Gd(s)·ΔIq + ΔVq/(Vdc/2) + ωL·ΔId/(Vdc/2)

将 ΔDd、ΔDq 代入桥臂电压小信号表达式,再代回电感线性化方程,整理后得到闭环导纳:

Yc(s) = Yo(s) / [1 + T(s)]
T(s) = Gc(s)·Gd(s)·(Vdc/2)·Yo(s)

这里 T(s) 是回路增益,Vdc/2 来自调制增益。若忽略交叉耦合,Yc(s) 可近似为标量形式 Yc(s) ≈ 1/(sL + Gc(s)Gd(s)Vdc/2),物理意义清晰:闭环后逆变器输出阻抗在低频段被控制环“压低”,呈现近似电流源特性;在高频段控制环增益衰减,输出阻抗回升至 sL,呈现电感特性。

四、PLL 与坐标变换的附加影响

上述推导假设 dq 坐标系与 PCC 电压严格同步。实际中 PLL 会引入额外的小信号耦合:当 PCC 电压发生扰动时,PLL 输出角度 Δθ 会随之波动,进而通过旋转矩阵影响 dq 电流的测量值。其小信号关系可写为:

Δθ = Gpll(s)·ΔVq
ΔId_meas = ΔId + Iq·Δθ
ΔIq_meas = ΔIq - Id·Δθ

把 Δθ 代入控制律后,闭环导纳矩阵会新增与 Iq·Gpll(s)、Id·Gpll(s) 相关的项。在单位功率因数运行时 Iq≈0,PLL 影响主要体现在 q 轴;但在弱电网下 PCC 电压畸变严重,PLL 带宽若过高,会显著改变中低频段导纳相位,甚至引入负阻尼。因此工程上常将 PLL 带宽限制在 20~50 Hz,并在阻抗模型中显式保留 Gpll(s)。

五、示例参数与结果判读

取一组教学示例参数(非实测):L = 2 mH,Vdc = 700 V,Vd = 311 V,Id = 20 A,Iq = 0 A,开关频率 fs = 10 kHz,采样周期 Ts = 100 μs,电流环 Kp = 8,Ki = 400,PLL 带宽 30 Hz。按上述公式在 MATLAB/PLECS 中扫频,可得到以下典型特征:

频段Yc 幅值趋势相位趋势物理含义
10 Hz~100 Hz较高,近似 1/(Kp·Vdc/2)接近 0°电流源特性,控制环主导
100 Hz~1 kHz随频率下降逐渐滞后PI 增益衰减,过渡区
1 kHz~5 kHz近似 1/(sL)接近 -90°电感特性,滤波电感主导
>5 kHz受延迟影响出现相位跌落可低于 -180°采样与调制延迟主导,易谐振

判读时重点关注两个量:一是 Zg/Zinv 的幅值交点频率处相位裕度是否大于 30°;二是 Yc 的相位是否在弱电网阻抗峰值附近出现负阻尼。若交点处相位裕度不足,可通过降低电流环带宽、增大 PLL 阻尼比或加入有源阻尼项来修正。需要强调的是,上述数值仅为教学示例,实际系统必须结合具体滤波器参数、死区与采样延迟重新扫频验证。

建模时务必保持“稳态量—小信号量—坐标系”三者一致。若在 dq 域推导却用 αβ 域的导纳矩阵代入,交叉耦合项会丢失,导致弱电网谐振频率预测偏差可达数十赫兹。