4. STATCOM数学模型建立

STATCOM 的控制策略、参数整定和稳定性分析,都建立在它的数学模型之上。如果只把 STATCOM 当成一个“能发无功的盒子”,就无法解释为什么 dq 轴电流可以分别控制有功和无功、为什么直流电压会随有功电流变化、为什么锁相环动态会影响输出。本节从三相 abc 时域模型出发,逐步推导到 dq 旋转坐标系模型,并给出有功无功解耦表示和直流电压动态方程,形成后续控制器设计的统一基础。

一、参照约定与电路结构

先明确符号约定,全文及图示保持一致。设电网三相电压为 vsa、vsb、vsc,STATCOM 交流侧三相输出电压为 vca、vcb、vcc,并网电流为 ia、ib、ic,方向取从电网流向 STATCOM 为正。连接电抗为 L,等效电阻为 R,直流侧电容为 C,直流电压为 vdc。三相电压采用对称正弦、相序 a→b→c 为正序,dq 变换采用幅值不变派克变换,d 轴与电网电压空间矢量对齐,q 轴超前 d 轴 90°。

二、abc 坐标系下的时域模型

对每一相列写基尔霍夫电压方程,可得三相时域模型:

L·di_a/dt = v_sa − v_ca − R·i_a
L·di_b/dt = v_sb − v_cb − R·i_b
L·di_c/dt = v_sc − v_cc − R·i_c

三相电流之和满足 ia+ib+ic=0(三线制无中线),因此只有两个电流状态独立。直流侧由三相桥臂开关函数 sa、sb、sc(取值 ±1/2 或 0/1,取决于建模约定)与交流电流耦合,电容电流为三相桥臂电流之和:

C·dv_dc/dt = s_a·i_a + s_b·i_b + s_c·i_c − i_loss

其中 iloss 表示开关损耗、导通损耗等等效直流侧损耗电流。abc 模型物理直观,但三相变量都是时变量,且相互耦合,直接用于控制器设计非常不便,因此需要变换到旋转坐标系。

三、dq 旋转坐标系变换

采用幅值不变派克变换,将 abc 三相量投影到以角频率 ω 旋转的 dq 坐标系。设 θ=ωt 为电网电压相位,由锁相环提供,则变换矩阵为:

[x_d]  2 ⎡ cosθ        cos(θ−120°)  cos(θ+120°) ⎤ [x_a]
[x_q] = ─ ⎢ −sinθ      −sin(θ−120°) −sin(θ+120°)⎥ [x_b]
[x_0]  3 ⎣  1/2         1/2          1/2         ⎦ [x_c]

三线制系统中零序分量 x0 为零,可忽略。将变换作用于 abc 时域方程,并利用 dθ/dt=ω,得到 dq 坐标系下的 STATCOM 交流侧模型:

L·di_d/dt = v_sd − v_cd − R·i_d + ωL·i_q
L·di_q/dt = v_sq − v_cq − R·i_q − ωL·i_d

这里出现了两个关键耦合项:+ωL·iq 出现在 d 轴方程,−ωL·id 出现在 q 轴方程。它们表示旋转坐标系下 d、q 轴之间的交叉耦合,是后续前馈解耦控制要抵消的对象。若 d 轴与电网电压矢量对齐,则 vsd=Vm、vsq=0,模型进一步简化。

STATCOM 数学模型:abc 时域 → dq 旋转 → 解耦与直流动态 abc 坐标系时域模型 L·di_a/dt = v_sa − v_ca − R·i_a L·di_b/dt = v_sb − v_cb − R·i_b L·di_c/dt = v_sc − v_cc − R·i_c C·dv_dc/dt = Σ s_k·i_k − i_loss 三相时变、相互耦合,不便于控制设计 Park 变换 θ = ωt (PLL) dq 旋转坐标系模型 L·di_d/dt = v_sd − v_cd − R·i_d + ωL·i_q L·di_q/dt = v_sq − v_cq − R·i_q − ωL·i_d d 轴对齐电网电压:v_sd = V_m, v_sq = 0 耦合项 ±ωL·i 需前馈抵消 稳态时 dq 量为直流量,便于 PI 控制 解耦 有功/无功解耦 P = 1.5·v_sd·i_d Q = −1.5·v_sd·i_q i_d 控有功/直流电压 i_q 控无功/交流电压 两通道近似独立 直流电压动态方程与功率平衡 C·dv_dc/dt = (P_ac − P_loss) / v_dc P_ac = 1.5·(v_cd·i_d + v_cq·i_q) 稳态:P_ac = P_loss,v_dc 恒定;i_d 偏差 → v_dc 变化 Δi_d ΔP Δv_dc → 外环 PI 调节 i_d* 注:图中参数为教学示例符号,实际取值由装置容量、电抗与开关频率决定

四、有功无功解耦表示

在 d 轴与电网电压对齐的条件下,STATCOM 与电网交换的瞬时功率可写成:

P = 1.5·(v_sd·i_d + v_sq·i_q) = 1.5·V_m·i_d
Q = 1.5·(v_sq·i_d − v_sd·i_q) = −1.5·V_m·i_q

这说明:d 轴电流主要对应有功功率,q 轴电流主要对应无功功率。当 iq 为负时,Q 为正,STATCOM 向电网注入感性无功(发出无功);当 iq 为正时,Q 为负,STATCOM 吸收无功。由于 Vm 近似恒定,只要分别控制 id 和 iq,就能独立调节有功和无功。这正是 dq 解耦控制的物理依据。

需要强调,解耦是“模型层面”的解耦,实际 d、q 轴方程中仍存在 ωL 交叉耦合项。工程上通常采用前馈解耦:在 d 轴电流控制器输出上叠加 −ωL·iq,在 q 轴控制器输出上叠加 +ωL·id,再叠加电网电压前馈 vsd、vsq,使被控对象近似为两个独立的一阶惯性环节。

五、直流电压动态方程

直流侧电容电压由交流侧有功功率与损耗功率之差决定。忽略开关谐波,用平均模型表示:

C·v_dc·dv_dc/dt = P_ac − P_loss
P_ac = 1.5·(v_cd·i_d + v_cq·i_q)

稳态时 Pac=Ploss,vdc 保持不变。若 id 出现偏差,交流侧有功功率变化,电容充放电,vdc 随之变化。因此直流电压外环通常以 vdc 为反馈,输出 id* 作为 d 轴电流内环指令;而 iq* 由无功或交流电压外环给出。这样,整个 STATCOM 就形成“直流电压外环 + 无功外环 + dq 电流内环”的双环控制结构。

教学示例:设 Vm=311 V(相电压峰值),L=2 mH,R=0.1 Ω,C=5000 μF,ω=314 rad/s。若要求发出 1 Mvar 感性无功,则 iq≈−1×106/(1.5×311)≈−2143 A(峰值),对应耦合项 ωL·iq≈314×0.002×(−2143)≈−1346 V,必须在 d 轴控制器中前馈补偿,否则 d 轴电流会被显著扰动。以上数值为教学示例,实际工程需按装置容量和电抗重新核算。

六、模型使用要点

  • 变换一致性:幅值不变与功率不变变换的系数不同,推导 P、Q 公式时必须与所用变换匹配,否则会出现 2/3 或 3/2 系数错误。
  • 锁相环影响:θ 由 PLL 提供,若 PLL 动态较慢或存在相位误差,d、q 轴将出现附加耦合,解耦效果变差。
  • 损耗等效:iloss 或 Ploss 可用等效电阻或恒流源近似,用于解释直流电压稳态偏差。
  • 离散化实现:数字控制器中 dq 变换、PI 调节和 PWM 更新均按采样周期 Ts 离散执行,模型推导时需注意采样延迟对电流环带宽的限制。