基于多馈入交直流系统关键直流电压稳定性的判别方法
【技术领域】
[0001] 本发明属于电力系统的稳定计算和分析领域，具体涉及一种基于多馈入交直流系 统关键直流电压稳定性的判别方法。
【背景技术】
[0002] 随着直流输电技术的发展和特高压电网的建设，我国已经出现并将继续出现多个 多馈入交直流系统。直流输电系统在为受端交流系统提供有功功率的同时，逆变站在换流 过程中需要吸收大量的无功功率，正常运行时，直流系统所需无功可由补偿装置完全提供。 但当交流系统发生故障时，特别是逆变侧交流系统发生故障及故障恢复期间，直流输电系 统需要从交流系统吸收大量无功功率，这可能成为交流系统电压失稳或崩溃的诱因。
[0003] 与纯交流系统和单馈入交直流相比，多馈入交直流系统有更多的动态元件与非线 性环节，引起系统电压稳定问题的因素更多也更复杂，如多回直流系统间的相互影响、直流 控制方式与控制参数的变化、无功补偿方式、发电机及励磁系统、负荷模型和有载调压变压 器(on-load tap changer，0LTC)等等。这些动态元件和非线性环节对电压稳定都起着重要 的作用。
[0004] 小干扰稳定分析方法是电力系统稳定性分析的一般性方法，同样适用于电压稳定 分析。小干扰电压稳定实际上是一种李雅普诺夫意义下的渐近稳定，它可以计及与电压稳 定问题有关的各元件的动态，其实质在于将所考虑的动态元件的微分方程在运行点处线性 化，通过分析状态方程特征矩阵的特征根来判断系统的稳定性和各元件的作用。许多文献 在小干扰电压稳定研究中考虑了发电机及励磁系统、0LTC、无功补偿设备及负荷的动态特 性，但考虑直流输电系统(HVDC)的动态研究非常少。并且，在基于小干扰电压稳定分析方法 判断多馈入交直流系统的关键直流方面，目前国内外尚未见相关的报道。
【发明内容】
[0005] 针对现有技术的不足，本发明基于小干扰电压稳定分析方法，提出基于多馈入交 直流系统关键直流电压稳定性的判别方法。
[0006] 为了实现上述发明目的，本发明采取如下技术方案：
[0007] 基于多馈入交直流系统关键直流电压稳定性的判别方法，所述方法包括如下步 骤：
[0008] 步骤1:采集多馈入交直流系统网络参数、动态元件模型及模型参数，构建多馈入 交直流系统的微分-代数方程组DAE;在系统稳态点(X〇，Y Q)处执行线性化处理，获取线性动 态系统模型；
[0009] 步骤2:基于小干扰稳定性分析方法求解系统特征值及相应的左、右特征向量，通 过熄弧角相关比Ργ?获取多馈入交直流系统中与稳定性相关的特征值λυ?( i = 1，2，…，m)，并 确定最小特征值λυ min;
[0010] 步骤3:利用所述最小特征值λυ_判断多馈入交直流系统的小干扰电压稳定性，根 据λυ min对应的参与因子Pj min确定影响多馈入交直流系统电压稳定性的关键直流。
[0011]优选的，所述步骤1中，构建多馈入交直流系统的微分-代数方程组DAE的表达式 为：
[0013 ]式中，X e Rn为系统状态变量;γ e Rm为代数变量;μ e R为变化参量;f和g分别为系统 的微分方程组和代数方程组;其中，微分方程表示元件、负荷节点和直流系统的动态特性； 代数方程表示由网络决定的潮流约束；
[0014]在系统平衡点(Χο，Υο)处将式（1)线性化，并消去变化参量μ，获取线性动态系统模 型：
[0017]如果D非奇异，消去式中的代数变量增量ΔΥ，获得标准状态矩阵方程如下描述：
[0019] 式中，X = A_BD-iG为ηΧη维系数矩阵，即该线性动态系统模型的状态矩阵。
[0020] 优选的，所述步骤2中，通过分析状态矩阵1的特征值，判断线性动态系统模型的 小扰动稳定性；
[0021 ]如果多馈入交直流系统规模较小，即状态变量低于1000，则采用QR算法计算状态 矩阵Α的全部特征值及相应的左、右特征向量;如果多馈入交直流系统规模较大，即状态变 量高于1000,则采用隐式重启动Arnoldi算法计算自定义区域内的特征值及其相应的左、右 特征向量。
[0022]优选的，所述步骤2的熄弧角相关比PYi表示状态矩阵A的第i个特征值h与直流熄 弧角γ的相关程度，将其定义为：
[0024]式中，XY j (j = 1，2，…，η)表示与逆变器熄弧角γ相关的状态变量;η为直流回数； Pki为量度第k个状态量Xk与第i个特征值λ?相关物理量；
[0025]将pklS义为参与因子，其表达式为：
[0027] 式中，Vi和Ui分别为与λ?对应的左、右特征向量;vki和Uki分别为左、右特征向量的 第k个元素;p kl的模| pkl |的大小反映 Xk和\的相关性大小；| Pkl |越大，则相关性越强。
[0028] 优选的，所述步骤2中，通过熄弧角相关比ργι获取多馈入交直流系统中与电压稳 定性相关的特征值的方法如下：
[0029]若根据步骤2求得的某个特征值Μ满足熄弧角相关比ργι>1，且虚部为〇,则认为该 特征值与多馈入交直流系统的小干扰电压稳定性相关;令λυ^ληλυ,表示多馈入交直流 系统的小干扰电压稳定程度。
[0030] 优选的，所述步骤2中，确定最小特征值λυ min包括，选取Aui(i = l，2，···，m)中距离 虚轴最近的特征值，即模I Aui |最小的特征值λυ min。
[0031] 优选的，所述步骤3中，利用最小特征值λυ min判断多馈入交直流系统的小干扰电压 稳定性的方法包括：
[0032] 若最小特征值λυ min>〇,则该多馈入交直流系统小干扰电压不稳定;若λυ min<〇， 且距离虚轴较远，则该多馈入交直流系统小干扰电压稳定;若AUmin<0,则该系统处于小干 扰电压稳定运行临界点。
[0033] 优选的，所述步骤3中，根据最小特征值λυ min对应的参与因子Pj min确定影响多馈 入交直流系统电压稳定性关键直流的方法包括：
[0034] 若系统中第j回直流逆变器熄弧角相关的状态变量XYj的最小特征值AUmin对应的 参与因子最大，则第j回直流为所述多馈入交直流系统的关键直流。
[0035] 与最接近的现有技术比，本发明的优异效果为：
[0036] 本发明可运用于电力系统的理论和计算分析，基于小干扰电压稳定分析方法，提 出基于多馈入交直流系统关键直流电压稳定性的判别方法，该方法能有效地找出影响多馈 入交直流系统电压稳定性的关键直流，从而判别系统中最易导致电压失稳的区域;对预防 大电网发生电压崩溃事故提供了新的技术支撑。
【附图说明】
[0037] 图1为本发明实施例中两馈入交直流算例等值系统示意图；
[0038] 图2为本发明实施例中等值系统稳态运行时的系统潮流图；
[0039] 图3为本发明提供的基于多馈入交直流系统关键直流电压稳定性的判别方法流程 图。
【具体实施方式】
[0040] 下面结合附图对本发明做进一步详细说明。
[0041] 基于多馈入交直流系统关键直流电压稳定性的判别方法，如图3所示，所述方法包 括如下步骤：
[0042] 实施例：以一个含两回直流输电线路的算例系统为例，研究判断多馈入交直流系 统中影响电压稳定性的关键直流的方法。两馈入交直流算例系统的等值电路如图1所示，两 回直流的逆变侧连接于相邻的同侧，逆变侧换流母线通过一回交流线路连接，直流系统整 流侧相互独立。直流系统1的主控制模式为恒功率控制，直流系统2的主控制模式为恒电流 控制。送端系统为无穷大机组，受端交流系统模型用系统等值电势源串联等值阻抗来表示。 稳态运行时的系统潮流图如图2所示。
[0043 ]步骤1:采集多馈入交直流系统网络参数、动态元件模型及模型参数，构建多馈入 交直流系统的微分-代数方程组DAE;在系统稳态点(Xo，Yo)处执行线性化处理，获取线性动 态系统模型；
[0044]步骤1中，构建多馈入交直流系统的微分-代数方程组DAE的表达式为：
[0046] 式中，X e Rn为系统状态变量;Y e Rm为代数变量;μ e R为变化参量;f和g分别为系统 的微分方程组和代数方程组;其中，微分方程表示元件、负荷节点和直流系统的动态特性； 代数方程表示由网络决定的潮流约束；
[0047] 在系统平衡点(X〇，Y〇)处将式（1)线性化，并消去变化参量μ，获取线性动态系统模 型：
[0049] 其中，系统状态变量增量△ X总共有20个，如下所示：
[0050] Δ X - [ Δ 0 rl，Α ，Λ CO r2，$2，Δ CO r3，$3，Δ Xa-reel，Α Χγ-invl，A Vmesl，
[0051 ] Λ Vdi-mesl，A Idr-mesl，A Idi-mesl，A 7" i-mesl, Δ Xa-rec2 ,
[0052] Δ Χγ-inv2，Δ Vmes2，Δ Vdi-mes2，A Idr-mes2，A Idi-mes2，A y i-mes2]
[0054]如果D非奇异，消去式中的代数变量增量ΔΥ，获得标准状态矩阵方程如下描述：
[0056]式中，A = A - BD C为η Xη维系数矩阵，即该线性动态系统模型的状态矩阵。 [