一种用于简化柔性交直流输电系统中RGA计算的方法

摘要

本发明公开了一种用于简化柔性交直流输电系统中RGA计算的方法，通过先对柔性交直流输电系统模型中的微分代数方程组进行线性化，再通过简明的矩阵运算得到传递函数，再在此基础上进行RGA计算。虽然形成的线性化系数矩阵元素较多，但避免了多次矩阵求逆导致的矩阵元素表达式过于复杂的问题，易于利用程序实现，提高了运算速度，降低了技术成本和对运行环境的要求，适用于高维数柔性交直流输电系统。

权利要求

1.一种用于简化柔性交直流输电系统RGA计算的方法，其包括以下步骤：步骤A：构建柔性交直流输电系统数学模型：

g(x，y，u)＝0

其中式 为描述系统各元件动态的微分方程，包括发电机及其励磁系统和FACTS装置的动态；式g(x，y，u)＝0为表示网络各个节点电压电流关系；

步骤B：将描述发电机和励磁系统动态的微分方程写成直角坐标形式并进行线性化，得电力系统的线性化模型：步骤C：对步骤B中的电力系统的线性化模型再次进行线性化可得：其中，Yphzt为2n×4m阶矩阵，Yphds为2n×2n阶矩阵，Yphkz为2n×2阶矩阵；

步骤D：对用于RGA计算的传递函数进行线性化，可得：式中，Yscsc为2×2阶矩阵，Ysczt为2×4m阶矩阵，Yscds为2×2n阶矩阵，Ysckz为2×2阶矩阵。

步骤E：联合步骤B、步骤C和步骤D中的三个线性化模型得出传递函数：-1 -1 -1 -1

其中，A＝Ywfzt-Ywfds*(Yphds Yphzt)；B＝Ywfkz-Ywfds*(Yphds Yphkz)C＝(Yscsc Yscds)*(Yphds Yp-1 -1 -1hzt)-Yscsc Ysczt；D＝(Yscsc Yscds)*(Yphds Yphkz)。

步骤F：将步骤E得到的传递函数转化为频域传递函数：

2.根据权利要求1所述的用于简化柔性交直流输电系统RGA计算的方法，其特征在于，发电机系统采用的是三阶实用模型：其中，δ为发电机功角；w为角速度；Pm为原动机输入的机械功率；Eq为发电机的暂态电势； 为发电机的暂态电抗；Xd为发电机的d轴同步电抗；Xq为发电机的q轴同步电抗；

Tj为发电机的惯性时间常数； 为发电机的直轴暂态开路时间常数；Efd为励磁电压；D为阻尼系数。

3.根据权利要求1所述的用于简化柔性交直流输电系统RGA计算的方法，其特征在于，励磁系统的数学模型为：

说明书

一种用于简化柔性交直流输电系统中RGA计算的方法

技术领域

[0001] 本发明涉及一种柔性交直流输电系统中RGA计算的方法，尤其涉及的是一种用于简化柔性交直流输电系统中RGA计算的方法。

背景技术

[0002] 柔性交流输电系统即将现代电力电子技术和现代自动控制技术引入到交直流输电系统中，以提高原有交直流输电系统的可控性、灵活性、运行的稳定性和经济性。该现代电力电子技术和现代自动控制技术即为柔性交流输电技术，即在传统的输电系统中安装FACTS(柔性交流输电)装置。然而，研究表明FACTS装置的不同控制回路之间存在着交互影响，这种交互影响对其本身的控制效果和系统的稳定性都有着较大的负面作用。因此，需要用到相对增益矩阵(RGA)方法对FACTS多个控制回路之间的交互影响进行量化分析，从而采取相应的技术手段以减少其负面作用。现有技术通常是对该柔性交直流输电系统模型的平衡方程进行节点消去从而得到传递函数，再在此基础上对柔性交直流输电系统模型进行RGA计算。该方法容易多次矩阵求逆导致的矩阵元素表达式过于复杂的问题，用程序实现比较复杂，计算成本高，尤其不适用于高维数电网模型的RGA计算。

[0003] 因此，现有技术还有待于改进和发展。

发明内容

[0004] 本发明的目的在于提供一种用于简化柔性交直流输电系统中RGA计算的方法，以解决现有技术中对该柔性交直流输电系统模型的平衡方程进行 节点消去从而得到传递函数，再在此基础上对柔性交直流输电系统模型进行RGA计算，从而导致的计算成本高，实现程序复杂，不适于高维数电网模型RGA计算的问题。

[0005] 本发明的技术方案如下：

[0006] 一种用于简化柔性交直流输电系统RGA计算的方法，其包括以下步骤：

[0007] 步骤A：构建柔性交直流输电系统数学模型：

[0008]

[0009] g(x，y，u)＝0

[0010] 其中式 为描述系统各元件动态的微分方程，包括发电机及其励磁系统和FACTS装置的动态；式g(x，y，u)＝0为表示网络各个节点电压电流关系；本步骤主要目的是将柔性交直流输电系统用数学模型表示出来，包含各发电机及其励磁系统、FACTS装置的动态特性，各个节点的电压、电流之间的关系。

[0011] 步骤B：将描述发电机和励磁系统动态的微分方程写成直角坐标形式并进行线性化，得电力系统的线性化模型：

[0012]

[0013]

[0014] 本步骤主要解决了高维数微分方程的求解问题，将其变化为线性化模型，从而可以将该方法应用于含有柔性交直流输电装置的大电网的此类问题的求解。

[0015] 步骤C：对步骤B中的电力系统的线性化模型再次进行线性化可得：

[0016]

[0017] 其中，Yphzt为2n×4m阶矩阵，Yphds为2n×2n阶矩阵，Yphkz为2n×2阶矩阵； [0018] 本步骤主要是线性化模型做技术处理，提高求解效率。

[0019] 步骤D：对用于RGA计算的传递函数进行线性化，可得：

[0020]

[0021]

[0022] 式中，Yscsc为2×2阶矩阵，Ysczt为2×4m阶矩阵，Yscds为2×2n阶矩阵，Ysckz为2×2阶矩阵。

[0023] 本步骤主要是对用于RGA计算的传递函数进行线性化，解决传递函数高维数，难以求解的问题。

[0024] 步骤E：联合步骤B、步骤C和步骤D中的三个线性化模型得出传递函数：

[0025]

[0026] 其中，A＝Ywfzt-Ywfds*(Yphds-Yphzt)；B＝Ywfkz-Ywfds*(Yphds-Yphkz)；C＝(Yscsc-1Yscds)*(Y-1 -1 -1 -1phds Yphzt)-Yscsc Ysczt；D＝(Yscsc Yscds)*(Yphds Yphkz)。

[0027] 本步骤将上述三个线性化模型综合变换，得出传递函数。利用本传递函数即可求解并分析柔性交直流输电系统中电网各元件相互影响的问题。

[0028] 步骤F：将步骤E得到的传递函数转化为频域传递函数：

[0029]

[0030] 本步骤主要是将传递函数转化为频域传递函数，将其变换为RGA求解的形式，利用RGA算法求解。

[0031] 所述的用于简化柔性交直流输电系统RGA计算的方法，其中，发电机系统采用的是三阶实用模型：

[0032]

[0033]

[0034]

[0035] 其中，δ为发电机功角；w为角速度；Pm为原动机输入的机械功率； 为发电机的暂态电势； 为发电机的暂态电抗；Xd为发电机的d轴同步电抗；Xq为发电机的q轴同步电抗；Tj为发电机的惯性时间常数； 为发电机的直轴暂态开路时间常数；Efd为励磁电压；D为阻尼系数。

[0036] 所述的用于简化柔性交直流输电系统RGA计算的方法，其中，励磁系统的数学模型为：

[0037]

[0038] 本发明通过先对柔性交直流输电系统模型中的微分代数方程组进行线性化，再通过简明的矩阵运算得到传递函数，再在此基础上进行RGA计算。虽然形成的线性化系数矩阵元素较多，但避免了多次矩阵求逆导致的矩阵元素表达式过于复杂的问题，易于利用程序实现，提高了运算速度，降低了技术成本和对运行环境的要求，适用于高维数柔性交直流输电系统。

具体实施方式

[0039] 为使本发明的目的、技术方案及优点更加清楚、明确，以下举实施例对本发明进一步详细说明。

[0040] 本发明中方法的基本流程如下：

[0041] 步骤A：构建柔性交直流输电系统的数学模型：

[0042]

[0043] g(x，y，u)＝0…………(2)

[0044] 式(1)为描述系统各元件动态的微分方程，包括发电机及其励磁系统和FACTS装置的动态；式(2)为描述网络各个节点电压电流关系的代数方程。

[0045] 本步骤主要目的是将柔性交直流输电系统用数学模型表示出来，包含各发电机及其励磁系统、FACTS装置的动态特性，各个节点的电压、电流之间的关系。

[0046] 柔性交直流输电系统动态过程可用式(1)和(2)的模型描述，x为状态变量，即节点电压的幅值和相角；y为代数变量，一般是网络中节点电压向量；u为控制变量，包括发电机有功、无功输出功率、发电机机端电压和变压器变比等。

[0047] 本发明提供的发电机系统采用的是三阶实用模型，如式(3)～(5)所表示的： [0048]

[0049]

[0050]

[0051] 其中，δ为发电机功角，w为角速度，Pm为原动机输入的机械功率， 为发电机的暂态电势、 为发电机的暂态电抗、Xd为发电机的d轴同步电抗、Xq为发电机的q轴同步电抗、Tj为发电机的惯性时间常数、 为发电机的直轴暂态开路时间常数；Efd为励磁电压；D为阻尼系数。

[0052] 其中，id和iq可由下式求出：

[0053]

[0054] id、iq分别表示d轴、q轴的电流(A)，通过坐标变换将电力系统从abc三相坐标变换到dq0坐标，参数ua，ub，uc，ia，ib，ic，磁链a，磁链b，磁链c这些量都变换到dq0坐标中，如果有需要可以逆变换回来。从物理意义上讲，该变换就是将ia，ib，ic电流投影等效到d轴和q轴上，将定子上的电流都等效到直轴和交轴上去。对于稳态来说，等效之后，iq，id正好就是一个常数。

[0055] 构建励磁系统的数学模型，如式(7)：

[0056]

[0057] 步骤B：假设多级系统中共有n个节点，m台发电机。对式(1)中描述发电机和励磁系统动态的微分方程，即发电机系统的三阶实用模型、以及励磁系统的数学模型写成直角坐标形式并进行线性化，可得到：

[0058]

[0059]

[0060] 其 中，Δδ[Δδ1Δδ2…Δδm]T；Δω ＝ [Δω1Δω2…Δωm]T；ΔE ′ q ＝T T T[ΔE′q1ΔE′q2…ΔE′qm] ；ΔEfd＝[ΔEfd1ΔEfd2…ΔEfdm] ；Vi＝[Vix Viy] ；Ywfzt为4m×4m阶矩阵，Ywfds为4m×2n阶矩阵，Ywfkz为4m×2阶矩阵。

[0061] 将描述网络各个节点电压电流关系的代数方程式(2)，改写成为：

[0062] I-(Y0+Y′)V＝0 (9)

[0063] 其中，I＝[Ix1Iy1…IxiIyi…IxnIyn]T，为节点注入电流向量；V＝[Vx1Vy1…VxiVyi…VxnVyn]T，为节点电压向量；Y0为网络节点导纳矩阵，Y′＝diag(Y1′，…，Yi′，…，Yn′)，是由发电机、负荷和FACTS装置并入网络的导纳所构成的矩阵。

[0064] 若节点i为发电机的节点，则其导纳为：

[0065]

[0066]

[0067] 式(10)中，Gxi、Bxi、Byi、Gyi为发电机i及其凸极效应后并入网络的导纳，其中Gxi、Bxi、Byi、Gyi通过式(11)求得：

[0068]

[0069] 若节点i安装STATCOM，则该节点的注入电流为：

[0070] I′＝I+IST (12)

[0071] 其中IST可以表示为：

[0072]

[0073] 若节点i安装SVC，则该节点的自导纳为：

[0074]

[0075] 本步骤主要解决了高维数微分方程的求解问题，将其变化为线性化模型，从而可以将该方法应用于含有柔性交直流输电装置的大电网的此类问题的求解。

[0076] 步骤C：式(8)是含状态变量、代数变量和控制变量的代数方程组，对其线性化可得：

[0077]

[0078] 式中，Yphzt为2n×4m阶矩阵，Yphds为2n×2n阶矩阵，Yphkz为2n×2阶矩阵。 [0079] 由式(8)、(15)构成了含FACTS装置的多机电力系统线性化模型。用于RGA计算