, 胡家兵1, 汪海蛟21. 华中科技大学 电气与电子工程学院, 强电磁工程与新技术国家重点实验室, 武汉 430074;
2. 中国电力科学研究院有限公司 新能源与储能运行控制国家重点实验室, 北京 100192
收稿日期：2020-11-30
基金项目：国家重点研发计划项目（2016YFB0900100）；国家电网科技项目（XTB17202000284）
作者简介：江克证(1995—), 男, 硕士研究生
通讯作者：朱建行, 博士后, E-mail：jh_zhu@hust.edu.cn
摘要：电力电子装备在电力系统中的广泛应用使得现代电力系统动态行为发生了显著变化，基于电网换相整流器的高压直流输电（LCC-HVDC）作为电力系统重要的组成部分，对电力系统的安全稳定运行具有重要影响。由于LCC-HVDC开关过程呈现断续时变特征，且开关频率与电磁尺度时间常数相当，因此断续时变开关过程的动态特性描述是研究LCC-HVDC电磁尺度动态稳定问题的关键。为此，该文提出了适用于计及开关过程的LCC-HVDC电磁尺度动态稳定分析的小信号模型。首先讨论了LCC-HVDC原始关系及其非线性断续周期时变特征，指出基本问题和挑战；然后基于线性周期时变理论建立了LCC-HVDC电磁尺度下的小信号模型。最后，通过与现有不考虑时变开关过程模型的对比分析，探究了开关过程对系统稳定性的影响，并通过时域仿真验证了分析的正确性。
关键词：基于电网换相整流器的高压直流输电(LCC-HVDC)开关过程周期时变谐波状态空间(HSS)电磁尺度
Small-signal modeling of LCC-HVDC systems with switching for electromagnetic timescale stability analyses of power systems
JIANG Kezheng1, ZHU Jianhang1
, HU Jiabing1, WANG Haijiao21. State Key Laboratory of Advanced Electromagnetic Engineering Technology, School of Electrical and Electronic Engineering, Huazhong University of Science and Technology, Wuhan 430074, China;
2. State Key Laboratory of Operation and Control of Renewable Energy & Storage Systems, China Electric Power Research Institute, Beijing 100192, China
Abstract: The wide application of power electronic equipment in power systems has caused significant changes in the dynamic behavior of modern power systems. As an important power system component, line commutated converter based high voltage direct current (LCC-HVDC) systems have an important impact on the safe, stable operation of power systems. Since the switching process in LCC-HVDC systems has intermittent time-varying characteristics and the switching frequency is equivalent to the electromagnetic timescale constant, these intermittent time-varying switching characteristics are the key to the stability of LCC-HVDC systems. This paper presents an electromagnetic timescale, small-signal model for LCC-HVDC systems that includes the switching process with a stability analysis. The model includes the nonlinear characteristics and the inherent periodic time-varying characteristics of LCC-HVDC systems for modeling the key challenges in the control system. Then, linear time-periodic theory is used to develop an electromagnetic timescale, small-signal model of the LCC-HVDC system. Finally, this model is compared with an existing model that does not consider the time-varying characteristics of the switching process to show the influence of the switching on the system stability with the results verified by time-domain simulations.
Key words: line commutated converter based high valtage direct current (LCC-HVDC)switching processperiodic time-varyingharmonic state space (HSS)electromagnetic timescale
随着电力电子技术在电力系统中的广泛应用，现代电力系统日益呈现出电力电子化的特征[1-2]。常规高压直流输电(LCC-HVDC)作为现代电力系统的重要组成部分，凭借高电压、大容量、远距离输送的经济性与灵活性高等优势被广泛应用到实际工程中[3-4]。对于中国的直流输电工程，南方电网已投运10回3 000万kW直流输电工程[5]，输送容量占远距离电能总容量的82%[6]，预计2030年将投运19回8 000万kW高压直流输电工程；而国家电网在“十四五”规划期间，华东、华中电网直流输电工程也将快速增加，其中华中电网中陕北至湖北、雅中至江西等特高压直流输电工程将于2021年相继投运[7]。LCC-HVDC的大量接入使得现代电力系统增加了大量响应时间常数不同的储能元件，进而导致现代电力系统的动态呈现多时间尺度的特征[8]，因此LCC-HVDC多时间尺度的动态分析是研究其对现代电力系统小扰动稳定影响的重要内容。
对于LCC-HVDC系统的建模与稳定性分析，文[9-11]基于交直流间的准稳态非线性关系，建立了LCC-HVDC的小信号模型，但忽略了开关动态和直流电流的动态，因此多用来分析LCC-HVDC对系统低频振荡的影响。文[12-13]利用基于Fourier级数展开和谐波平衡原理的动态相量方法，将原始的非线性周期时变系统转换到非线性时不变系统，从而在直流工作点上进行线性化处理，进而可以利用经典线性时不变理论进行小信号分析。当考虑高次分量时，由于控制器的非线性及开关函数的分段周期特性，动态相量建模的复杂性将急剧增加，因此该模型在实际应用时仅考虑基频和直流分量的动态，且多用来分析系统次同步振荡的相关问题。由于LCC-HVDC采用半控型开关器件，且开关频率与电磁尺度时间常数相当，因此LCC-HVDC开关过程的描述是分析该装备对电力系统电磁尺度动态稳定问题的关键，而现有文献并未对此深入分析。
根据LCC-HVDC的基本工作原理，半控型开关器件按照控制器输出触发角周期地开通、关断各开关管，且不同开关组合下的电压电流方程不同，因此该装备的原始方程呈现出断续周期时变的特征。另外，该装备原始断续方程的分段点受控制器输出影响，且控制器内部同时包含坐标变换等非线性环节，因此LCC-HVDC的原始方程又呈现出非线性特征。综上所述，对LCC-HVDC原始关系中非线性与周期时变性的理解、非线性关系的线性化和周期时变性的处理是该装备并网下系统电磁尺度小信号建模的基本挑战。因此本文在系统原始断续非线性时变方程的稳态周期轨迹上直接进行线性化处理，从而得到LCC-HVD C线性周期时变模型，进而利用文[14]所提的谐波状态空间方法，即利用Fourier级数展开和谐波平衡原理将线性周期时变模型变换为线性时不变模型，该方法得到的模型又称为谐波状态空间(HSS)模型。通过与现有仅考虑基频分量的动态相量模型的对比分析，探究了控制器参数对系统稳定性的影响，进而表征了LCC-HVDC开关过程对系统电磁尺度动态稳定问题的影响。
1 LCC-HVDC的基本结构和原始模型1.1 LCC-HVDC换流器的基本结构本文主要以LCC-HVDC整流站(见图 1)为例，主要包括换流器、控制器、直流线路，交流主回路和交流滤波器5部分，其中ucr和udcr分别为换流器交流母线和直流端口电压，ir和idcr分别为换流器交流和直流端口电流，Ldcr和Ldci为直流侧平波电抗器，Rdcr和Rdci为直流T型等效线路电阻，Cdc为T型等效线路电容，各部分详细拓扑结构及参数参考Benchmark标准模型[15]。图 2为控制器拓扑结构，控制方案采用定电流控制，锁相环控制为脉冲触发提供参考相位，与电流控制的输出共同产生触发序列[16-17]。
 |
| 图 1 LCC-HVDC拓扑结构图 |
| 图选项 |
 |
| 图 2 LCC-HVDC控制器结构图 |
| 图选项 |
1.2 LCC-HVDC原始数学关系1) 换流器的原始数学关系。
换流器是LCC-HVDC系统的重要组成部分，工程中多采用十二脉波换流器，输入是三相交流母线电压及直流电流，输出为三相交流电流与直流电压，因此基于开关函数的交直流关系可表示为
| $\left\{\begin{array}{l}u_{\mathrm{dcr}}=s_{\mathrm{rva}} u_{\mathrm{cra}}+s_{\mathrm{rvb}} u_{\mathrm{crb}}+s_{\mathrm{rvc}} u_{\mathrm{crc}}, \\i_{\mathrm{ra}}=s_{\mathrm{ria}} i_{\mathrm{dcr}}, \\i_{\mathrm{rb}}=s_{\mathrm{rib}} i_{\mathrm{dcr}}, \\i_{\mathrm{rc}}=s_{\mathrm{ric}} i_{\mathrm{der}} .\end{array}\right.$ |
由于开关器件的周期分段导通特性，使得用于描述开关器件导通状态的开关函数是分段周期函数，本文用0和1分别表示开关器件晶闸管的关断和导通状态，用0.5表示晶闸管的换相过程。值得注意的是，当晶闸管导通状态发生改变时，由于换相电感的存在，换相电流不能发生突变，故在换相过程中，换流器交直流电流具体表达式与交直流电压具体表达式不同。换相过程中交直流电流之间的关系可以通过斜率为±1/μ的一次函数进行线性近似处理[18]。以十二脉波换流器阀1电流开关函数为例进行说明，其具体表达式为
| ${s_{i1}} = \left\{ {\begin{array}{*{20}{l}}{(\theta - \alpha )/\mu ,}\\{\theta \in [\alpha + 2k{\rm{ \mathsf{ π} }},\alpha + \mu + 2k{\rm{ \mathsf{ π} }}];}\\{1,}\\{\theta \in [\alpha + \mu + 2k{\rm{ \mathsf{ π} }},\alpha + 2{\rm{ \mathsf{ π} }}/3 + 2k{\rm{ \mathsf{ π} }}];}\\{1 - (\theta - \alpha - n)/\mu ,}\\{\theta \in [\alpha + 2{\rm{ \mathsf{ π} }}/3 + 2k{\rm{ \mathsf{ π} }},\alpha + 2{\rm{ \mathsf{ π} }}/3 + 2k{\rm{ \mathsf{ π} }} + \mu ];}\\{0,}\\{\theta \in [\alpha + 2{\rm{ \mathsf{ π} }}/3 + 2k{\rm{ \mathsf{ π} }} + \mu ,\alpha + 2(k + 1){\rm{ \mathsf{ π} }}].}\end{array}} \right.$ |
图 3中，阀1电流开关函数对应的波形以2π为周期，且开关函数是关于触发角的函数。由于触发角是直流电流控制器的输出，因此电流开关函数本质上是关于直流电流的函数即si1(idcr)，从而导致开关过程的原始关系具有非线性特征。交直流电压表达式虽然不同于交直流电流表达式，但交直流电压关系同样体现出非线性特征，在此不再详细展开。
 |
| 图 3 电流开关函数波形示意图 |
| 图选项 |
2) 控制器的原始数学关系。
图 2中，控制器包括电流控制器和锁相环控制器2部分。电流控制器根据直流电流实际值与指令值之差进行PI调节，输出触发控制角，其输入输出关系为：
| $\begin{array}{*{20}{c}}{{i_{{\rm{drf}}}} = \frac{1}{{1 + {T_{\rm{r}}}s}}{i_{{\rm{dcr}}}},}\\{{\alpha _{{\rm{ord}}}} = {\rm{ \mathsf{ π} }} - \left( {{k_{{\rm{pcc}}}} + \frac{{{k_{{\rm{icc}}}}}}{s}} \right)\left( {{i_{{\rm{dref}}}} - {i_{{\rm{drr}}}}} \right).}\end{array}$ |
锁相环控制器通过实时监测并网点相位为控制和脉冲生成提供同步信号基准，其中锁相环输入和输出分别为端电压相位和锁相环输出相位，其原始关系为
| $\left\{\begin{array}{l}\boldsymbol{u}_{\alpha \beta}=\boldsymbol{T} \boldsymbol{u}_{\mathrm{cr}}, \\\boldsymbol{u}_{d q}=\boldsymbol{P}_{\theta} \boldsymbol{u}_{\alpha \beta}, \\\theta_{\mathrm{pll}}=\left(k_{\mathrm{ppll}} / s+k_{\mathrm{ipll}} / s^{2}\right) u_{\mathrm{q}} .\end{array}\right.$ |
当系统受到干扰时，电网电压相位会有一定的动态变化，锁相环的输出相位也会有一定的动态变化。由于锁相环输出相位为LCC-HVDC脉冲触发提供了参考相位，因此其动态变化将影响实际触发角。图 4给出了锁相环输出相位对实际触发角的影响，其中ucra′和ucrc′分别为ucra和ucrc受扰之后的量。
 |
| 图 4 锁相环对实际触发角的影响 |
| 图选项 |
由此可得LCC-HVDC实际触发角与锁相环输出相位及定电流控制输出触发角指令值之间的关系为
| $\alpha=\alpha_{\text {ord }}+\theta_{\text {pcc }}-\theta_{\text {pll }} \text { . }$ |
3) 直流线路的数学模型。
LCC-HVDC直流线路采用集中参数模型，其拓扑结构图如图 1所示，因此直流侧电压电流关系可以表示为：
| $\begin{array}{c}L_{\mathrm{dcr}} \frac{\mathrm{d} i_{\mathrm{dcr}}}{\mathrm{d} t}=u_{\mathrm{dcr}}-u_{\mathrm{dc}}-R_{\mathrm{dcr}} i_{\mathrm{dcr}}, \\C_{\mathrm{dc}} \frac{\mathrm{d} u_{\mathrm{dc}}}{\mathrm{d} t}=i_{\mathrm{dcr}}-i_{\mathrm{dci}}, \\L_{\mathrm{dci}} \frac{\mathrm{d} i_{\mathrm{dci}}}{\mathrm{d} t}=u_{\mathrm{dc}}-u_{\mathrm{dci}}-R_{\mathrm{dc}} i_{\mathrm{dci}}.\end{array}$ |
LCC-HVDC交流系统主要由交流侧主回路和交流滤波器组成，图 5以a相为例给出了整流站交流系统的单相结构图，分别选取电感电流及电容电压为状态变量，交流侧主回路与交流滤波器状态方程为：
| $\begin{array}{*{20}{c}}{\frac{{\rm{d}}}{{{\rm{d}}t}}{u_{{\rm{cra}}}} = {K_{11}}{u_{{\rm{cra}}}} + {K_{13}}{u_{{\rm{c}}3{\rm{ra}}}} + {K_{14}}{u_{{\rm{ctra}}}} + }\\{{K_{15}}{i_{2{\rm{ra}}}} + {K_{16}}{i_{3{\rm{ra}}}} + {K_{17}}{i_{1{\rm{ra}}}} + {K_{18}}{s_{{\rm{ria}}}}{i_{{\rm{dcr}}}},}\\{\frac{{\rm{d}}}{{{\rm{d}}t}}{u_{{\rm{c}}2{\rm{ra}}}} = {K_{25}}{i_{2{\rm{ra}}}},}\\{\frac{{\rm{d}}}{{{\rm{d}}t}}{u_{{\rm{c}}3{\rm{ra}}}} = {K_{31}}{u_{{\rm{cra}}}} + {K_{33}}{u_{{\rm{c}}3{\rm{ra}}}} + {K_{35}}{i_{2{\rm{ra}}}},}\\{\frac{{\rm{d}}}{{{\rm{d}}t}}{u_{{\rm{c4ra}}}} = {K_{41}}{u_{{\rm{cra}}}} + {K_{44}}{u_{{\rm{c4ra}}}},}\\{\frac{{\rm{d}}}{{{\rm{d}}t}}{i_{2{\rm{ra}}}} = {K_{51}}{u_{{\rm{cra}}}} + {K_{52}}{u_{{\rm{c}}2{\rm{ra}}}} + {K_{53}}{u_{{\rm{c}}3{\rm{ra}}}} + {K_{55}}{i_{2{\rm{ra}}}},}\\{\frac{{\rm{d}}}{{{\rm{d}}t}}{i_{3{\rm{ra}}}} = {K_{61}}{u_{{\rm{cra}}}} + {K_{64}}{u_{{\rm{c4ra}}}},}\\{\frac{{\rm{d}}}{{{\rm{d}}t}}{i_{1{\rm{ra}}}} = {K_{71}}{u_{{\rm{cra}}}} + {K_{77}}{i_{1{\rm{ra}}}} + {K_{78}}{u_{{\rm{sra}}}}.}\end{array}$ |
 |
| 图 5 整流站交流系统单相结构图 |
| 图选项 |
其中K11~K78均为常数(见表 1)。由于系数K18项引入了开关函数，使得LCC-HVDC原始关系状态方程中含有周期时变量。
表 1 系统常数
| 参数 | 数值 | 参数 | 数值 | |
| K11 | -(R22ra+R31ra)/(R22raR31raC1ra) | K44 | -1/ (R31raC4ra) | |
| K13 | 1/ (R22raC1ra) | K46 | 1/C4ra | |
| K14 | 1/ (R31raC1ra) | K51 | 1/L2ra | |
| K15 | -1/C1ra | K52 | -1/L2ra | |
| K16 | -1/C1ra | K53 | -1/L2ra | |
| K17 | 1/C1ra | K55 | -R21ra/L2ra | |
| K18 | -1/C1ra | K61 | 1/L3ra | |
| K25 | 1/C2ra | K64 | -1/L3ra | |
| K31 | 1/ (R22raC3ra) | K71 | -1/L1ra | |
| K33 | -1/ (R22raC3ra) | K77 | -R12ra/L1ra | |
| K35 | 1/C3ra | K78 | 1/L1ra | |
| K41 | 1/ (R31raC4ra) | K81 | -1/C1ra |
表选项
2 基于线性周期时变理论的小信号建模基于LCC-HVDC原始关系的认识可知该系统具有非线性周期时变特征，传统的将其变换到时不变系统下的建模思路将面临巨大的挑战，因此本文提出在系统原始方程的稳态周期轨迹上直接进行线性化处理，从而在更为一般化的线性体系下建立LCC-HVDC小信号模型即线性周期时变模型。现有针对分段线性周期时变模型的分析仍面临巨大的挑战，因此本文利用Fourier级数展开和谐波平衡原理，将LCC-HVDC线性周期时变(LTP)模型变换为无穷阶线性时不变模型即HSS模型，从而可以借助经典的线性时不变系统的稳定性分析方法对系统进行分析[19-21]，基本思路如图 6所示。
 |
| 图 6 线性周期时变理论分析思路 |
| 图选项 |
2.1 LCC-HVDC的HSS小信号模型基于上述线性周期时变理论的建模方法对LCC-HVDC原始数学模型进行处理，得到LCC-HVDC的HSS小信号模型。
换流器部分HSS小信号模型可以表示为
| $\left\{\begin{aligned}\Delta \boldsymbol{u}_{\mathrm{dc}}=& \boldsymbol{\varGamma}\left(u_{\mathrm{cra0}}\right) \Delta \boldsymbol{s}_{\mathrm{ra}}+\boldsymbol{\varGamma}\left(s_{\mathrm{rr0}}\right) \Delta \boldsymbol{u}_{\mathrm{cra}}+\boldsymbol{\varGamma}\left(u_{\mathrm{crb} 0}\right) \Delta \boldsymbol{s}_{\mathrm{rb}}+\\& \boldsymbol{\varGamma}\left(s_{\mathrm{rb0}}\right) \Delta \boldsymbol{u}_{\mathrm{crb}}+\boldsymbol{\varGamma}\left(u_{\mathrm{crc0}}\right) \Delta \boldsymbol{s}_{\mathrm{rc}}+\boldsymbol{\varGamma}\left(s_{\mathrm{rc} 0}\right) \Delta \boldsymbol{u}_{\mathrm{crc}}, \\\Delta \boldsymbol{i}_{\mathrm{ra}}=& \boldsymbol{\varGamma}\left(i_{\mathrm{dc0}}\right) \Delta \boldsymbol{s}_{\mathrm{ra}}+\boldsymbol{\varGamma}\left(s_{\mathrm{ra} 0}\right) \Delta \boldsymbol{i}_{\mathrm{dc}}, \\\Delta \boldsymbol{i}_{\mathrm{rb}}=& \boldsymbol{\varGamma}\left(i_{\mathrm{de0}}\right) \Delta \boldsymbol{s}_{\mathrm{rb}}+\boldsymbol{\varGamma}\left(s_{\mathrm{rb} 0}\right) \Delta \boldsymbol{i}_{\mathrm{dc}}, \\\Delta \boldsymbol{i}_{\mathrm{rc}}=& \boldsymbol{\varGamma}\left(i_{\mathrm{d} c 0}\right) \Delta \boldsymbol{s}_{\mathrm{rc}}+\boldsymbol{\varGamma}\left(s_{\mathrm{rc} 0}\right) \Delta \boldsymbol{i}_{\mathrm{dc}} .\end{aligned}\right.$ |
| $\begin{array}{l}\mathit{\boldsymbol{ \boldsymbol{\varGamma} }}\left( {{i_{{\rm{dco}}}}} \right) = \left[ {\begin{array}{*{20}{c}}{}& \vdots & \vdots & \vdots &{}\\ \cdots &{{i_{{\rm{dc}}0(0)}}}&{{i_{{\rm{dc}}0( - 1)}}}&{{i_{{\rm{dc}}0( - 2)}}}& \cdots \\ \cdots &{{i_{{\rm{dc}}0(1)}}}&{{i_{{\rm{dc}}0(0)}}}&{{i_{{\rm{dc}}0( - 1)}}}& \cdots \\ \cdots &{{i_{{\rm{dc0(2)}}}}}&{{i_{{\rm{dc}}0(1)}}}&{{i_{{\rm{dc}}0(0)}}}& \cdots \\{}& \vdots & \vdots & \vdots &{}\end{array}} \right],\\\begin{array}{*{20}{c}}{\Delta {\mathit{\boldsymbol{u}}_{{\rm{dc}}}} = }\\{{{\left[ {\Delta {u_{{\rm{dc}}( - k)}}, \cdots ,\Delta {u_{{\rm{dc}}( - 1)}},\Delta {u_{{\rm{dc}}(0)}},\Delta {u_{{\rm{dc}}(1)}}, \cdots ,\Delta {u_{{\rm{dc}}(k)}}} \right]}^{\rm{T}}}.}\end{array}\end{array}$ |
定直流电流控制器环节的HSS小信号模型为:
| $\begin{array}{c}(s+\boldsymbol{N}) \Delta \boldsymbol{i}_{\mathrm{drf}}=\left(\Delta \boldsymbol{i}_{\mathrm{dc}}-\Delta \boldsymbol{i}_{\mathrm{drf}}\right) / T_{\mathrm{r}}, \\(s+\boldsymbol{N}) \Delta \boldsymbol{x}_{\mathrm{CC}}=k_{\mathrm{icc}}\left(\Delta \boldsymbol{i}_{\mathrm{dref}}-\Delta \boldsymbol{i}_{\mathrm{drf}}\right) ,\\\Delta \boldsymbol{\alpha}_{\text {ord }}=-\left(\Delta \boldsymbol{x}_{\mathrm{CC}}+k_{\mathrm{pcc}}\left(\Delta \boldsymbol{i}_{\mathrm{dref}}-\Delta \boldsymbol{i}_{\mathrm{drf}}\right)\right) .\end{array}$ |
| $\boldsymbol{N}=\operatorname{diag}(-\mathrm{j} k \omega, \cdots,-\mathrm{j} \omega, 0, \mathrm{j} \omega, \cdots, \mathrm{j} k \omega).$ |
| $\begin{array}{c}\left\{\begin{aligned}\Delta \boldsymbol{u}_{\alpha} &=\boldsymbol{\Gamma}(T[1,1]) \Delta \boldsymbol{u}_{\mathrm{cra}}+\boldsymbol{\Gamma}(T[1,2]) \Delta \boldsymbol{u}_{\mathrm{crb}} \\&+\boldsymbol{\Gamma}(T[1,3]) \Delta \boldsymbol{u}_{\mathrm{crc}} ,\\\Delta \boldsymbol{u}_{\beta} &=\boldsymbol{\Gamma}(T[2,1]) \Delta \boldsymbol{u}_{\mathrm{cra}}+\boldsymbol{\Gamma}(T[2,2]) \Delta \boldsymbol{u}_{\mathrm{crb}} \\&+\boldsymbol{\Gamma}(T[2,3]) \Delta \boldsymbol{u}_{\mathrm{crc}}, \\\Delta \boldsymbol{u}_{\mathrm{d}} &=\boldsymbol{\Gamma}\left(P_{\theta}[1,1]\right) \Delta \boldsymbol{u}_{\alpha}+\boldsymbol{\Gamma}\left(P_{\theta}[1,2]\right) \Delta \boldsymbol{u}_{\beta} \\&+\boldsymbol{\Gamma}\left(u_{\mathrm{q} 0}\right) \Delta \boldsymbol{\theta}_{\mathrm{pll}}, \\\Delta \boldsymbol{u}_{\mathrm{q}} &=\boldsymbol{\Gamma}\left(P_{\theta}[2,1]\right) \Delta \boldsymbol{u}_{\alpha}+\boldsymbol{\Gamma}\left(P_{\theta}[2,2]\right) \Delta \boldsymbol{u}_{\beta} \\&-\boldsymbol{\Gamma}\left(u_{\mathrm{d} 0}\right) \Delta \boldsymbol{\theta}_{\mathrm{pll}} .\end{aligned}\right.\\\left\{\begin{array}{l}(s+\boldsymbol{N}) \Delta \boldsymbol{\theta}_{\mathrm{pll}}=\Delta \boldsymbol{\omega}, \\(s+\boldsymbol{N}) \Delta \boldsymbol{x}_{\mathrm{pll}}=k_{\mathrm{ipll}} \Delta \boldsymbol{u}_{\mathrm{q}} ,\\\Delta \boldsymbol{\omega}=\Delta \boldsymbol{x}_{\mathrm{pll}}+k_{\mathrm{ppll}} \Delta \boldsymbol{u}_{\mathrm{q}}.\end{array}\right.\end{array}$ |
| $\begin{array}{c}(s+\boldsymbol{N}) \Delta \boldsymbol{i}_{\mathrm{dcr}}=\frac{1}{L_{\mathrm{dcr}}} \Delta \boldsymbol{u}_{\mathrm{dcr}}-\frac{1}{L_{\mathrm{dcr}}} \Delta \boldsymbol{u}_{\mathrm{dc}}-\frac{R_{\mathrm{dcr}}}{L_{\mathrm{dcr}}} \Delta \boldsymbol{i}_{\mathrm{dcr}} ,\\(s+\boldsymbol{N}) \Delta \boldsymbol{u}_{\mathrm{dc}}=\frac{1}{C_{\mathrm{dc}}} \Delta \boldsymbol{i}_{\mathrm{dcr}}-\frac{1}{C_{\mathrm{dc}}} \Delta \boldsymbol{i}_{\mathrm{dci}} ,\\(s+\boldsymbol{N}) \Delta \boldsymbol{i}_{\mathrm{dci}}=\frac{1}{L_{\mathrm{dci}}} \Delta \boldsymbol{u}_{\mathrm{dc}}-\frac{1}{L_{\mathrm{dci}}} \Delta \boldsymbol{u}_{\mathrm{dci}}-\frac{R_{\mathrm{dci}}}{L_{\mathrm{dci}}} \Delta \boldsymbol{i}_{\mathrm{dci}}.\end{array}$ |
| $\begin{array}{*{20}{c}}{(s + \mathit{\boldsymbol{N}})\Delta {\mathit{\boldsymbol{u}}_{{\rm{cra}}}} = {K_{11}}\Delta {\mathit{\boldsymbol{u}}_{{\rm{cra}}}} + {K_{13}}\Delta {\mathit{\boldsymbol{u}}_{c3{\rm{ra}}}} + {K_{14}}\Delta {\mathit{\boldsymbol{u}}_{{\rm{c}}4{\rm{ra}}}} + }\\{{K_{15}}\Delta {\mathit{\boldsymbol{i}}_{2{\rm{ra}}}} + {K_{16}}\Delta {\mathit{\boldsymbol{i}}_{3{\rm{ra}}}} + {K_{17}}\Delta {\mathit{\boldsymbol{i}}_{1{\rm{ra}}}} + {K_{18}}\Delta {\mathit{\boldsymbol{i}}_{{\rm{ra}}}},}\\{(s + \mathit{\boldsymbol{N}}){\mathit{\boldsymbol{u}}_{{\rm{c}}2{\rm{ra}}}} = {K_{25}}\Delta {i_{2{\rm{ra}}}},}\\{(s + N){\mathit{\boldsymbol{u}}_{{\rm{c}}3{\rm{ra}}}} = {K_{31}}\Delta {\mathit{\boldsymbol{u}}_{{\rm{cra}}}} + {K_{33}}\Delta {\mathit{\boldsymbol{u}}_{{\rm{c}}3{\rm{ra}}}} + {K_{35}}\Delta {\mathit{\boldsymbol{i}}_{2{\rm{ra}}}},}\\{(s + \mathit{\boldsymbol{N}}){\mathit{\boldsymbol{u}}_{{\rm{c}}4{\rm{ra}}}} = {K_{41}}\Delta {\mathit{\boldsymbol{u}}_{{\rm{cra}}}} + {K_{44}}\Delta {\mathit{\boldsymbol{u}}_{{\rm{c}}4{\rm{ra}}}},}\\{(s + \mathit{\boldsymbol{N}}){\mathit{\boldsymbol{i}}_{2{\rm{ra}}}} = {K_{51}}\Delta {\mathit{\boldsymbol{u}}_{{\rm{cra}}}} + {K_{52}}\Delta {\mathit{\boldsymbol{u}}_{{\rm{c}}2{\rm{ra}}}} + }\\{{K_{53}}\Delta {\mathit{\boldsymbol{u}}_{{\rm{c}}3{\rm{ra}}}} + {K_{55}}\Delta {\mathit{\boldsymbol{i}}_{2{\rm{ra}}}},}\\{(s + \mathit{\boldsymbol{N}}){\mathit{\boldsymbol{i}}_{3{\rm{ra}}}} = {K_{61}}\Delta {\mathit{\boldsymbol{u}}_{{\rm{cra}}}} + {K_{64}}\Delta {\mathit{\boldsymbol{u}}_{{\rm{ctra}}}},}\\{(s + \mathit{\boldsymbol{N}}){\mathit{\boldsymbol{i}}_{1{\rm{ra}}}} = {K_{71}}\Delta {\mathit{\boldsymbol{u}}_{{\rm{cra}}}} + {K_{77}}\Delta {\mathit{\boldsymbol{i}}_{1{\rm{ra}}}} + {K_{78}}\Delta {\mathit{\boldsymbol{u}}_{{\rm{sra}}}}.}\end{array}$ |
表 2 LCC-HVDC系统参数
| 参数 | 数值 |
| 三相交流电压 | 345 kV |
| 交流主回路电感 | 0.02 H |
| 交流主回路电阻 | 0.056 Ω |
| 变压器变比 | 380/213 |
| 直流平波电抗器 | 0.596 8 H |
| 直流线路电阻 | 2.5 Ω |
| 直流线路电容 | 26 μF |
| 直流电流控制器kpcc | 1.098 9 |
| 直流电流控制器kicc | 91.575 |
| 锁相环控制器kppll | 10 |
| 锁相环控制器kipll | 50 |
表选项
与现有的HSS理论建模类似，LCC-HVDC的HSS小信号模型同样具有无穷阶的特性，不利于具体的分析和计算，因此首先需要进行截断处理，本文对HSS小信号模型进行13阶截断处理。由于LCC-HVDC系统中各变量在HSS模型中以各变量Fourier级数幅值的形式存在，因此各变量实际的时域波形可以通过Fourier级数进行还原：
| $\begin{array}{c}\Delta x(t)=\boldsymbol{F}(t) \Delta \boldsymbol{x}_{\mathrm{hss}}, \\\boldsymbol{F}(t)=\left[\mathrm{e}^{-\mathrm{j} k \omega t}, \cdots, \mathrm{e}^{-\mathrm{j} \omega t}, 0, \mathrm{e}^{\mathrm{j} \omega t}, \cdots, \mathrm{e}^{\mathrm{j} k \omega t}\right] .\end{array}$ |
 |
| 图 7 HSS小信号模型与详细非线性时域模型触发角时域波形 |
| 图选项 |
3 LCC-HVDC小信号稳定性分析由于在对HSS模型进行分析之前，需要进行高阶截断以保证模型的准确性，因此HSS模型的特征根具有特征根列的特征(见图 8)，其中特征根列的实部相同，虚部相差基频的整数倍。由于截断的影响，位于特征根列两端的特征根会部分偏移出特征根列，因此本文选择特征根列的中心特征根s0来表征系统特征。
 |
| 图 8 HSS模型系统特征值轨迹示意图 |
| 图选项 |
3.1 电流控制器带宽对系统稳定性的影响电流控制器是实现LCC-HVDC正常运行的主要控制器之一，其参数对换流器的调节性能具有决定性作用，是影响系统动态性能的关键因素之一。
在阻尼比保持不变的情况下，电流控制器带宽变化对HSS小信号模型稳定性的影响如图 9所示。
 |
| 图 9 LCC-HVDC系统HSS模型特征值轨迹图 |
| 图选项 |
随着带宽的增加，λ3、λ4和λ5、λ6单调向左半平面运动，表现为稳定性增强的趋势；λ7、λ8和λ9、λ10基本不受带宽变化的影响；λ1、λ2和λ11、λ12单调向右半平面移动，表现为稳定性逐步降低的趋势，但未到达复平面右半平面；λ13、λ14单调向右半平面移动并穿越虚轴到达右半平面，系统由稳定状态过渡到不稳定状态。为了探究开关过程在电流控制器带宽变化下对LCC-HVDC系统稳定性的影响，对比分析了电流控制器带宽同等条件变化下，基频动态相量模型的特征值变化，如图 10所示。随着带宽的增加，γ1、γ2和γ3、γ4基本不受带宽变化的影响，且系统在带宽变化过程中未出现不稳定状态。2种模型分析结果的不同，初步表明了电流控制器带宽变化下计及开关过程的LCC-HVDC建模对系统电磁尺度稳定性分析非常必要。
 |
| 图 10 LCC-HVDC系统基频动态相量模型特征值轨迹图 |
| 图选项 |
为了验证上述分析的正确性，在详细非线性时域模型中，保持阻尼比不变，当带宽由上述范围的最小带宽值切换到最大带宽值时，进行直流电流的时域仿真，结果如图 11所示。此时LCC-HVDC系统由稳定状态逐步发散，与上述HSS小信号模型分析结果相一致，而图 10基频动态相量模型的分析结果与图 11仿真结果不符，证明了考虑开关过程建模的必要性。
 |
| 图 11 LCC-HVDC详细非线性时域模型直流电流时域波形 |
| 图选项 |
3.2 锁相环控制器带宽对系统稳定性的影响锁相环为LCC-HVDC系统的导通触发提供参考相位，也是影响系统动态性能的关键因素之一。
在阻尼比保持不变的情况下，锁相环控制器带宽变化对LCC-HVDC系统稳定性的影响如图 12所示。随着带宽的增加，λ11、λ12和λ15、λ16单调向左半平面运动，表现为稳定性增强的趋势；λ5、λ6和λ7、λ8基本不受带宽变化的影响；λ1、λ2和λ3、λ4单调向右半平面移动，表现为稳定性逐步降低的趋势，但未到达复平面右半平面；λ9、λ10的运动轨迹先向右移动再向左移动，始终在左半平面；λ13、λ14单调向左半平面移动并由右半平面穿越虚轴到达左半平面，系统由不稳定状态过渡到稳定状态。为了探究开关过程在锁相环控制器带宽变化下对LCC-HVDC系统稳定性的影响，对比分析了锁相环控制器带宽同等条件变化下基频动态相量模型的特征值变化，结果如图 13所示。随着带宽的增加，γ1、γ2和γ3、γ4单调向左半平面运动，表现为稳定性增强的趋势；γ6、γ7和γ8基本不受带宽变化的影响；γ5单调向右半平面移动，表现为稳定性逐步降低的趋势，但未到达复平面右半平面；锁相环带宽变化过程中未出现不稳定状态。由此初步表明了锁相环带宽变化下计及开关过程的LCC-HVDC建模对系统电磁尺度稳定性分析具有重要影响。
 |
| 图 12 LCC-HVDC系统HSS模型特征值轨迹图 |
| 图选项 |
 |
| 图 13 LCC-HVDC系统基频动态相量模型特征值轨迹图 |
| 图选项 |
为了验证上述分析的正确性，在详细非线性时域模型中，保持阻尼比不变，当控制器带宽由上述范围的最大带宽值切换到最小带宽值时，进行直流电流的时域仿真，结果如图 14所示。此时LCC-HVDC系统由稳定状态逐步发散，与上述HSS小信号模型分析结果一致，而图 13基频动态相量模型的分析结果与图 14仿真结果不符，证明了计及开关过程的建模的必要性。
 |
| 图 14 LCC-HVDC详细非线性时域模型直流电流时域波形 |
| 图选项 |
4 结论本文讨论了LCC-HVDC输入、输出原始关系及其非线性断续周期时变特征，指出该装备电磁尺度的基本问题和挑战，并基于线性周期时变理论建立了计及LCC-HVDC开关过程的电磁尺度HSS小信号模型。将HSS小信号模型分析结果与传统基于基频动态相量的小信号模型分析结果进行了对比，探究了控制器参数对系统稳定性的影响，证实了LCC-HVDC开关过程的动态特性对电力系统电磁尺度稳定性具有重要影响；并进一步说明对于电力系统电磁尺度稳定性问题，传统LCC-HVDC装备模型将导致对系统稳定性的揭示发生本征错误；该结论与时域仿真结果吻合。
参考文献
| [1] | 周孝信, 鲁宗相, 刘应梅, 等. 中国未来电网的发展模式和关键技术[J]. 中国电机工程学报, 2014, 34(29): 4999-5008. ZHOU X X, LU Z X, LIU Y M, et al. Development models and key technologies of future grid in China[J]. Proceedings of the CSEE, 2014, 34(29): 4999-5008. (in Chinese) |
| [2] | 周孝信, 陈树勇, 鲁宗相. 电网和电网技术发展的回顾与展望——试论三代电网[J]. 中国电机工程学报, 2013, 36(22): 1-11. ZHOU X X, CHEN S Y, LU Z X. Review and prospect for power system development and related technologies: A concept of three-generation power systems[J]. Proceedings of the CSEE, 2013, 36(22): 1-11. (in Chinese) |
| [3] | 赵畹君. 高压直流输电工程技术[M]. 2版. 北京: 中国电力出版社, 2011. ZHAO W J. High voltage direct current transmission engineering technology[M]. 2nd ed. Beijing: China Electric Power Press, 2011. (in Chinese) |
| [4] | KIMBARK E W, Direct current transmission[M]. Direct current transmission[M]. New York: Wiley, 1971. |
| [5] | 陈国平, 李明节, 许涛, 等. 我国电网支撑可再生能源发展的实践与挑战[J]. 电网技术, 2017, 41(10): 3095-3103. CHEN G P, LI M J, XU T, et al. Practice and challenge of renewable energy development based on interconnected power grids[J]. Power System Technology, 2017, 41(10): 3095-3103. (in Chinese) |
| [6] | 国家电网公司电力调度控制中心. 2017年国家电网公司系统电网运行方式[R]. 北京: 国家电网公司, 2017. State Grid Corporation of China Power Dispatching Control Center. State grid corporation of China system grid operation mode in 2017[R]. Beijing: State Grid Corporation of China, 2017. (in Chinese) |
| [7] | 全球能源互联网发展合作组织. 中国"十四五"电力发展规划研究[R]. 北京: 全球能源互联网发展合作组, 2020. Global Energy Internet Development Cooperation Organization. Research on China's "14th Five-Year" electric power development plan[R]. Beijing: Global Energy Internet Development Cooperation Organization, 2020. (in Chinese) |
| [8] | 袁小明, 程时杰, 胡家兵. 电力电子化电力系统多尺度电压功角动态稳定问题[J]. 中国电机工程学报, 2016, 36(19): 5145-5154, 5395. YUAN X M, CHENG S J, HU J B. Multi-time scale voltage and power angle dynamics in power electronics dominated large power systems[J]. Proceedings of the CSEE, 2016, 36(19): 5145-5154, 5395. (in Chinese) |
| [9] | 郭春义, 宁琳如, 王虹富, 等. 基于开关函数的LCC-HVDC换流站动态模型及小干扰稳定性[J]. 电网技术, 2017, 41(12): 3862-3870. GUO C Y, NING L R, WANG H F, et al. Switching-function based dynamic model of LCC-HVDC station and small signal stability analysis[J]. Power System Technology, 2017, 41(12): 3862-3870. (in Chinese) |
| [10] | OSAUSKAS C M, WOOD A R. Small-signal dynamic modeling of HVDC systems[J]. IEEE Transactions on Power Delivery, 2003, 18(1): 220-225. DOI:10.1109/TPWRD.2002.803843 |
| [11] | KARAWITA C, ANNAKKAGE U D. Multi-infeed HVDC interaction studies using small-signal stability assessment[J]. IEEE Transactions on Power Delivery, 2009, 24(2): 910-918. DOI:10.1109/TPWRD.2008.2002986 |
| [12] | 戚庆茹, 焦连伟, 严正, 等. 高压直流输电动态相量建模与仿真[J]. 中国电机工程学报, 2003, 23(12): 31-35. QI Q R, JIAO L W, YAN Z, et al. Modeling and simulation of HVDC with dynamic phasors[J]. Proceedings of the CSEE, 2003, 23(12): 31-35. (in Chinese) |
| [13] | DARYABAK M, FILIZADEH S, JATSKEVICH J, et al. Modeling of LCC-HVDC systems using dynamic phasors[J]. IEEE Transactions on Power Delivery, 2014, 29(4): 1989-1998. DOI:10.1109/TPWRD.2014.2308431 |
| [14] | WERELEY N. Analysis and control of linear periodically time varying systems[D]. Boston, MA: Massachusetts Institute of Technology, 1991. |
| [15] | SZECHTMAN M. First benchmark model for HVDC control studies[J]. Electra, 1991, 135(4): 54-73. |
| [16] | LU J, YUAN X M, HU J B, et al. Motion equation modeling of LCC-HVDC stations for analyzing DC and AC network interactions[J]. IEEE Transactions on Power Delivery, 2020, 35(3): 1563-1574. DOI:10.1109/TPWRD.2019.2947572 |
| [17] | HU Y P, LU J, ZHOU H F. Commutation failure analysis considering direct current dynamics in LCC-HVDC systems[C]//12th IET International Conference on AC and DC Power Transmission (ACDC 2016). Beijing: IET, 2016: 1-5. |
| [18] | HU L, YACAMINI R. Harmonic transfer through converters and HVDC links[J]. IEEE Transactions on Power Electronics, 1992, 7(3): 514-525. DOI:10.1109/63.145139 |
| [19] | ZHU J H, HU J B, MA S C. Modeling and analysis of MMC in AC current control timescale considering PLL dynamics[C]//2019 IEEE Power & Energy Society General Meeting (PESGM). Atlanta, GA: IEEE, 2019: 1-5. |
| [20] | ZHU J H, MA S C, GONG X, et al. Mechanism underlying constant amplitude oscillation of grid-tied voltage source converters[C]//The 8th Renewable Power Generation Conference (RPG 2019). Shanghai: IET, 2019: 1-6. |
| [21] | ZHU J H, HU J B, WANG S C, et al. Small-signal modeling and analysis of MMC under unbalanced grid conditions based on linear time-periodic (LTP) method[J]. IEEE Transactions on Power Delivery, 2021, 36(1): 205-214. DOI:10.1109/TPWRD.2020.2976776 |
