跳到主要内容

跟网型储能变流器-平均模型-v1

案例介绍

跟网型储能变流器-平均模型由平均化电气拓扑、锁相环、功率外环、电流内环、调制环节等模块构成,还附加了标幺化、倍乘等值、量测与输出、多短路比测试实现,直流侧为储能电池,后续可拓展斩波控制、穿越控制、脱网控制等模块,模型整体架构如下图所示。

跟网/构网型变流器并网系统整体架构示意图
跟网/构网型变流器并网系统整体架构示意图

模型已封装为元件,具备单元测试回路,可单独使用,也可在其他算例中调用。

元件图形
元件图形

使用方法说明

适用场景

跟网型储能变流器-平均模型支持单机或接入大规模电力系统算例的仿真测试,适用于以下分析场景:

  • 小信号稳定性分析
  • 不同短路比下的变流器并网适应性测试
  • 考虑系统功角、频率、电压的电网稳定性分析
  • 高比例新能源接入下的系统支撑能力评估
  • 跟网型与构网型控制策略的对比研究

适用范围

  • 建议步长范围:1-50μs
  • 当前参数设置下建议短路比范围:>3

功能概述

  • 变流器并网系统的发电单元除储能外,还可以接入风机、光伏或SVG,以满足不同的研究需求
  • 跟网型变流器控制可以附加频率-有功和电压-无功的支撑控制,提升系统的稳定性和可靠性

算例介绍

跟网型储能变流器-平均模型由电气主拓扑、标幺制、倍乘等值、变流器控制、量测与输出等五个模块组成。

电气主拓扑由蓄电池、直流电容、受控电流源组成的直流侧和受控电压源、交流滤波器、电压源组成的交流侧,以及功率设定和单元测试构成,实现平均化变流器等效建模。

  • 功率设定输出由用户设定的输出有功、无功功率参考值,并包含输出小信号扰动功率的功能

  • 单元测试中设置经线路阻抗与理想电压源相连,连接阻抗的大小由用户设置,系统的短路比在量测与输出模块中计算并输出。

    电气主拓扑
    电气主拓扑
    功率设定与单元测试
    功率设定与单元测试

标幺制处理使仿真中除时间外的电气量均为标幺值,便于参数的通用化设置。为使得所构建的变流器模型适用于各种电压等级及功率送出场景,需要构造一个物理概念清晰、使用方便的变流器模型标幺制系统,将控制环节中的量测值、参考值、元件参数、控制参数等均转换为标幺值,经过调试和优化可以得到一套通用的控制参数,实现变流器模型标幺化。

标幺制系统
标幺制系统

倍乘等值的基本原理为:按设备额定电压、额定容量、并联数量变化倍数线性缩放电阻、电感、电容等参数,并保持电压、电流标幺值不变(有名值需要折算),可在简化模型规模的同时,使等效模型在并网点呈现与真实多机系统相似的功率输出和动态响应,从而满足系统级稳定性分析的需求。

倍乘等值
倍乘等值

变流器控制由锁相环、Park变换、功率外环控制、电流内环控制、调制模块、平均化控制等部分组成,实现输出功率的控制,此外还可附加支撑控制,提升系统的稳定性和可靠性。平均化变流器的控制信号通过功率平衡法计算得到。

变流器控制
变流器控制

量测与输出模块量测控制环节所需电气量,计算并输出用户关注的电气量。

量测与输出  

算例仿真测试

针对跟网型储能变流器-平均模型分别进行了小信号状态空间模型验证和并网适应性测试。

小信号状态空间模型验证

依据跟网型储能变流器-平均模型构建了全阶小信号状态空间模型,如下式所示。

Δx˙=AΔx+BΔu\Delta \dot x = A\Delta x + B\Delta u

其中:

x=[ζ,  δ,  γd,  γq,  ξd,  ξq,  id,  iq,  iLd,  iLq,  vd,  vq]Tx = {\left[ {\zeta ,\;\delta ,\;{\gamma _d},\;{\gamma _q},\;{\xi _d},\;{\xi _q},\;{i_d},\;{i_q},\;{i_{Ld}},\;{i_{Lq}},\;{v_d},\;{v_q}} \right]^{\rm{T}}} u=[Pref,Qref]Tu = {\left[ {{P_{{\rm{ref}}}},{Q_{{\rm{ref}}}}} \right]^{\rm{T}}} A=[0000000000011Tp0000000000Kp000000vd0vq000id0iq0000000vq0vd000iq0id0001To1000Ko1vd0Ko1vq010Ko1id0Ko1iq00001To100Ko1vq0Ko1vd001Ko1iq0Ko1id0iq0Tpvg0sin(δ0)Lg0000RgLgω0001LgKpiq0id0Tpvg0cos(δ0)Lg0000ω0RgLg000KpLgid01Lg00Ki1LfTo101LfTi10Ki1Ko1vd0LfKi1Ko1vq0LfKi1+RfLf0Ki1Ko1id0LfKi1Ko1iq0Lf000Ki1LfTo101LfTi1Ki1Ko1vq0LfKi1Ko1vd0Lf0Ki1+RfLfKi1Ko1iq0LfKi1Ko1id0Lfvq0Tp000001Cf01Cf00Cfω0+CfKpvq0Cfvd0Tp0000001Cf01Cfω0Kpvd0]\tiny{A} = \begin{bmatrix} 0&0&0&0&0&0&0&0&0&0&0&1\\[4mm] {\dfrac{1}{{{T_p}}}}&0&0&0&0&0&0&0&0&0&0&{{K_p}}\\[4mm] 0&0&0&0&0&0&{ - {v_{d0}}}&{ - {v_{q0}}}&0&0&{ - {i_{d0}}}&{ - {i_{q0}}}\\[4mm] 0&0&0&0&0&0&{{v_{q0}}}&{ - {v_{d0}}}&0&0&{ - {i_{q0}}}&{{i_{d0}}}\\[4mm] 0&0&{\dfrac{1}{{{T_{o1}}}}}&0&0&0&{ - {K_{o1}}{v_{d0}}}&{ - {K_{o1}}{v_{q0}}}&{ - 1}&0&{ - {K_{o1}}{i_{d0}}}&{ - {K_{o1}}{i_{q0}}}\\[4mm] 0&0&0&{\dfrac{1}{{{T_{o1}}}}}&0&0&{{K_{o1}}{v_{q0}}}&{ - {K_{o1}}{v_{d0}}}&0&{ - 1}&{ - {K_{o1}}{i_{q0}}}&{{K_{o1}}{i_{d0}}}\\[4mm] {\dfrac{{{i_{q0}}}}{{{T_p}}}}&{\dfrac{{{v_{g0}}\sin ({\delta _0})}}{{{L_g}}}}&0&0&0&0&{ - \dfrac{{{R_g}}}{{{L_g}}}}&{{\omega _0}}&0&0&{\dfrac{1}{{{L_g}}}}&{{K_p}{i_{q0}}}\\[4mm] { - \dfrac{{{i_{d0}}}}{{{T_p}}}}&{\dfrac{{{v_{g0}}\cos ({\delta _0})}}{{{L_g}}}}&0&0&0&0&{ - {\omega _0}}&{ - \dfrac{{{R_g}}}{{{L_g}}}}&0&0&0&{ - \dfrac{{{K_p}{L_g}{i_{d0}} - 1}}{{{L_g}}}}\\[4mm] 0&0&{\dfrac{{{K_{i1}}}}{{{L_f}{T_{o1}}}}}&0&{\dfrac{1}{{{L_f}{T_{i1}}}}}&0&{ - \dfrac{{{K_{i1}}{K_{o1}}{v_{d0}}}}{{{L_f}}}}&{ - \dfrac{{{K_{i1}}{K_{o1}}{v_{q0}}}}{{{L_f}}}}&{ - \dfrac{{{K_{i1}} + {R_f}}}{{{L_f}}}}&0&{ - \dfrac{{{K_{i1}}{K_{o1}}{i_{d0}}}}{{{L_f}}}}&{ - \dfrac{{{K_{i1}}{K_{o1}}{i_{q0}}}}{{{L_f}}}}\\[4mm] 0&0&0&{\dfrac{{{K_{i1}}}}{{{L_f}{T_{o1}}}}}&0&{\dfrac{1}{{{L_f}{T_{i1}}}}}&{\dfrac{{{K_{i1}}{K_{o1}}{v_{q0}}}}{{{L_f}}}}&{ - \dfrac{{{K_{i1}}{K_{o1}}{v_{d0}}}}{{{L_f}}}}&0&{ - \dfrac{{{K_{i1}} + {R_f}}}{{{L_f}}}}&{ - \dfrac{{{K_{i1}}{K_{o1}}{i_{q0}}}}{{{L_f}}}}&{\dfrac{{{K_{i1}}{K_{o1}}{i_{d0}}}}{{{L_f}}}}\\[4mm] {\dfrac{{{v_{q0}}}}{{{T_p}}}}&0&0&0&0&0&{ - \dfrac{1}{{{C_f}}}}&0&{\dfrac{1}{{{C_f}}}}&0&0&{\dfrac{{{C_f}{\omega _0} + {C_f}{K_p}{v_{q0}}}}{{{C_f}}}}\\[4mm] { - \dfrac{{{v_{d0}}}}{{{T_p}}}}&0&0&0&0&0&0&{ - \dfrac{1}{{{C_f}}}}&0&{\dfrac{1}{{{C_f}}}}&{ - {\omega _0}}&{ - {K_p}{v_{d0}}} \end{bmatrix} B=[0010Ko1000Ki1Ko1Lf00000010Ko1000Ki1Ko1Lf00]TB = \begin{bmatrix} 0&0&{1{\kern 1pt} }&{0{\kern 1pt} }&{{K_{o1}}{\kern 1pt} }&{0{\kern 1pt} }&{0{\kern 1pt} }&{0{\kern 1pt} }&{\dfrac{{{K_{i1}}{K_{o1}}}}{{{L_f}}}{\kern 1pt} }&{0{\kern 1pt} }&{0{\kern 1pt} }&0\\ 0&0&{0{\kern 1pt} }&{ - 1{\kern 1pt} }&{0{\kern 1pt} }&{ - {K_{o1}}{\kern 1pt} }&{0{\kern 1pt} }&{0{\kern 1pt} }&{0{\kern 1pt} }&{ - \dfrac{{{K_{i1}}{K_{o1}}}}{{{L_f}}}{\kern 1pt} }&{0{\kern 1pt} }&0 \end{bmatrix}^{\rm{T}}

参数列表如下:

符号物理意义取值
Kp锁相环比例系数1.00p.u.
Tp锁相环积分时间常数0.01s
Ko功率外环比例系数0.01p.u.
To功率外环积分时间常数0.01s
Ki电流内环比例系数0.10p.u.
Ti电流内环积分时间常数0.10s
Kω调频系数0p.u.
Ku调压系数0p.u.
Rf0原始滤波电阻2.00mΩ
Lf0原始滤波电感5.00mH
Cf0原始滤波电容73.10μF
V0原始额定电压0.381kV
SB0原始单机额定容量0.05MVA
N并联数量10
VB额定电压0.40kV
SB1单机额定容量2.00MVA
Pref有功功率参考值0.90p.u.
Qref无功功率参考值0.10p.u.
SCR系统短路比5.0
X/R系统电抗电阻比5.0

由小信号状态空间模型计算得到并网系统特征根分布如下图所示,验证了系统的稳定性。

特征根分布
特征根分布

小扰动下跟网型储能变流器-平均模型状态空间方程理论计算结果与电磁暂态仿真结果基本相同,各电气量均方根误差均小于5×1045×10^{-4},验证了模型的准确性。

小扰动验证
小扰动验证

并网适应性测试

跟网型储能变流器-平均模型进行了多短路比下的并网适应性测试。由下图可知,随着系统短路比增大,系统稳定性单调增强,而系统稳定性随电抗电阻比变化非单调,同一短路比下电抗电阻比过小或过大系统稳定性均减弱,不同短路比下的最佳电抗电阻比不同,例如该图中系统最稳定工况为SCR=10,X/R=9,而当短路比小于等于3时系统均失稳。

跟网型变流器并网系统主导特征根实部热力图
跟网型变流器并网系统主导特征根实部热力图
机端功率
机端功率

模型地址

点击打开模型地址:跟网型储能变流器-平均模型

附录

参数

基本设置

基本设置

参数名键名类型 [单位]描述
并联数量num实数并联储能设备数量
并网点电压Vpcc实数 [kV]交流端口线电压有效值
电压初始相位Init_Phase实数 [Deg]电压初始相位
单机额定容量Sbase实数 [MVA]单机额定容量(1~3MVA)
有功功率参考值Ppu_ref_PCC实数 [p.u.]有功功率参考值
无功功率参考值Qpu_ref_PCC实数 [p.u.]无功功率参考值

跟网控制参数设置

跟网控制参数设置

参数名键名类型 [单位]描述
频率-有功支撑系数Kf_VSC实数 [p.u.]频率-有功支撑系数(取零则为无支撑)
电压-无功支撑系数Ku_VSC实数 [p.u.]电压-无功支撑系数(取零则为无支撑)
功率外环比例系数Ko_VSC实数VSC 控制中,功率外环 PI 环节的比例系数
功率外环积分时间常数To_VSC实数 [s]VSC 控制中,功率外环 PI 环节的积分时间常数
功率外环限幅Olim_VSC实数 [p.u.]VSC 控制中,功率外环 PI 环节的限幅
电流内环比例系数Ki_VSC实数VSC 控制中,电流内环 PI 环节的比例系数
电流内环积分时间常数Ti_VSC实数 [s]VSC 控制中,电流内环 PI 环节的积分时间常数
电流内环限幅Ilim_VSC实数 [p.u.]VSC 控制中,电流内环 PI 环节的限幅

潮流设置

潮流设置

参数名键名类型 [单位]描述
Bus TypeBusType选择节点类型
Injected Active Powerpf_P实数 [MW]节点注入有功功率
Injected Reactive Powerpf_Q实数 [MVar]节点注入无功功率
Bus Voltage Magnitudepf_V实数 [p.u.]母线电压幅值
Bus Voltage Anglepf_Theta实数 [Deg]母线电压相位
Lower Voltage Limitpf_Vmin实数 [p.u.]母线电压下限
Upper Voltage Limitpf_Vmax实数 [p.u.]母线电压上限
Lower Reactive Power Limitpf_Qmin实数 [MVar]无功功率下限
Upper Reactive Power Limitpf_Qmax实数 [MVar]无功功率上限

事件设置

事件设置

参数名键名类型 [单位]描述
并网时间OnTime实数 [s]储能并网运行时刻
启动时间Startup_Time实数 [s]启动时间
功率设定改变时间ChargeSwitch实数 [s]小扰动时间

量测设置

量测设置

参数名键名类型 [单位]描述
有功量测(MW)P_PCC虚拟引脚(输出)有功量测(MW)
无功量测(MVar)Q_PCC虚拟引脚(输出)无功量测(MVar)

单元测试

单元测试

参数名键名类型 [单位]描述
UnitTestUnitTest布尔单元测试开关,打开时接入外部测试电路。
线路电阻R实数 [Ω]线路等效电阻
线路电感L实数 [H]线路等效电感

引脚

引脚名键名类型维度描述
AC0电气3 × 1AC