一种直流输电系统的离散特征值分析方法及装置与流程

1.本发明涉及电力系统分析技术领域，尤其涉及一种直流输电系统的离散特征值分析方法及装置。


背景技术：

2.高压直流输电因其线路造价较低、线路有功损耗较小等优势，现下被广泛应用于远距离、大容量电流输送中。但随着电力系统“双高”特性日益凸显，直流输电系统宽频带振荡的风险加剧，会对电力运输造成影响，因此有必要对直流输电系统进行稳定性分析，评估系统振荡风险。
3.特征值法因其严谨的理论基础、较高的准确性、信息的丰富性，被广泛应用于直流输电系统稳定性分析中。但是，建立直流输电系统连续状态空间模型时没有分元件模块化建模，而是将全系统视为一个整体考虑，使得系统各组成部分之间缺乏独立性；此外，为消除连续状态空间模型的中间变量，需要针对每一个中间变量补充列写相关代数方程，该过程亦较为繁琐，且消元过程为纯代数运算，没有实际物理意义。


技术实现要素：

4.本发明目的在于，提供一种直流输电系统的离散特征值分析方法及装置，使得离散特征值分析在直流输电系统中也能得到运用。
5.为实现上述目的，本发明提供一种直流输电系统的离散特征值分析方法，所述直流输电系统，包括直流输电元件、换流变和交流系统等效模型，所述直流输电元件包括换流器、直流支撑电容、直流电缆、换流器出口等效电感和电阻，所述离散特征值分析方法包括：
6.以整流侧公共连接点处的电压和逆变侧公共连接点处的电压为输入变量，以整流侧公共连接点处的电流和逆变侧公共连接点处的电流为输出变量，得到非线性动态方程；
7.对所述非线性动态方程做全微分线性化处理，得到稳态运行点处的小信号状态空间模型；其中，所述小信号状态空间模型包括直流输电元件的小信号连续状态空间方程、直流输电元件整流侧和逆变侧端口输出电流表达式；
8.将所述小信号状态空间模型基于自身dq坐标系的端口量转换到系统基准两相同步旋转xy坐标系下，得到旋转坐标状态空间模型；
9.采用梯形积分规则对所述旋转坐标状态空间模型进行离散化处理，得到直流元件的离散状态空间方程和离散等效电路模型。
10.优选地，所述非线性动态方程的表达式为：
[0011][0012]
其中，x
dc
为连续状态变量构成的矩阵，u
dcrdq
为整流侧公共连接点处的电压，u
dcidq
为逆变侧公共连接点处的电压，i
dcrdq
整流侧公共连接点处的电流，i
dcidq
为逆变侧公共连接点处的电流，p为微分算子，t为方程所描述的时刻，f(x)表示px
dc
与x
dc
、u
dcrdq
、u
dcidq
之间的函
数关系，g(x)表示i
dcrdq
、i
dcidq
与x
dc
、u
dcrdq
、u
dcidq
之间的函数关系。
[0013]
优选地，所述小信号状态空间模型的表达式为：
[0014][0015]
其中，δ表示微增量，a
dc1
、b
dcr1
、b
dci1
、c
dc1
、d
dcr1
、d
dci1
、c
dc2
、d
dcr2
、d
dci2
为系数矩阵，pδx
dc
表示状态变量微增量的微分，δx
dc
为状态变量微增量，δu
dcrdq
为整流侧公共连接点处的电压微增量，δu
dcidq
为逆变侧公共连接点处的电压微增量，δi
dcrdq
为整流侧公共连接点处的电流微增量，δi
dcidq
为逆变侧公共连接点处的电流微增量。
[0016]
优选地，所述采用梯形积分规则对所述旋转坐标状态空间模型进行离散化处理，得到直流元件的离散状态空间方程和离散等效电路模型，包括：
[0017]
将t时刻直流元件离散状态变量用前一时刻连续状态变量、整流侧端口输入电压和逆变侧端口输入电压线性表示，迭代得到t时刻端口电流、离散状态变量和端口电压的线性关系；
[0018]
根据所述线性关系，得到直流元件的等效电导和两个电流源构成的等效电路，根据直流元件的等效电导和两个电流源构成的等效电路得到离散等效电路模型；
[0019]
根据所述线性关系，得离散状态变量的迭代表达式，根据离散状态变量的迭代表达式得到离散状态空间方程。
[0020]
本发明还提供一种直流输电系统的离散特征值分析装置，所述直流输电系统，包括直流输电元件、换流变和交流系统等效模型，所述直流输电元件包括换流器、直流支撑电容、直流电缆、换流器出口等效电感和电阻，其特征在于，所述离散特征值分析装置包括：
[0021]
电压分析模块，用于以整流侧公共连接点处的电压和逆变侧公共连接点处的电压为输入变量，以整流侧公共连接点处的电流和逆变侧公共连接点处的电流为输出变量，得到非线性动态方程；
[0022]
线性化处理模块，用于对所述非线性动态方程做全微分线性化处理，得到稳态运行点处的小信号状态空间模型；其中，所述小信号状态空间模型包括直流输电元件的小信号连续状态空间方程、直流输电元件整流侧和逆变侧端口输出电流表达式；
[0023]
坐标系旋转模块，用于将所述小信号状态空间模型基于自身dq坐标系的端口量转换到系统基准两相同步旋转xy坐标系下，得到旋转坐标状态空间模型；
[0024]
离散化处理模块，用于采用梯形积分规则对所述旋转坐标状态空间模型进行离散化处理，得到直流元件的离散状态空间方程和离散等效电路模型。
[0025]
优选地，所述非线性动态方程的表达式为：
[0026][0027]
其中，x
dc
为连续状态变量构成的矩阵，u
dcrdq
为整流侧公共连接点处的电压，u
dcidq
为逆变侧公共连接点处的电压，i
dcrdq
整流侧公共连接点处的电流，i
dcidq
为逆变侧公共连接点处的电流，p为微分算子，t为方程所描述的时刻，f(x)表示px
dc
与x
dc
、u
dcrdq
、u
dcidq
之间的函数关系，g(x)表示i
dcrdq
、i
dcidq
与x
dc
、u
dcrdq
、u
dcidq
之间的函数关系。
[0028]
优选地，所述小信号状态空间模型的表达式为：
[0029][0030]
其中，δ表示微增量，a
dc1
、b
dcr1
、b
dci1
、c
dc1
、d
dcr1
、d
dci1
、c
dc2
、d
dcr2
、d
dci2
为系数矩阵，pδx
dc
表示状态变量微增量的微分，δx
dc
为状态变量微增量，δu
dcrdq
为整流侧公共连接点处的电压微增量，δu
dcidq
为逆变侧公共连接点处的电压微增量，δi
dcrdq
为整流侧公共连接点处的电流微增量，δi
dcidq
为逆变侧公共连接点处的电流微增量。
[0031]
优选地，所述离散化处理模块还包括：
[0032]
线性关系分析单元，用于将t时刻直流元件离散状态变量用前一时刻连续状态变量、整流侧端口输入电压和逆变侧端口输入电压线性表示，迭代得到t时刻端口电流、离散状态变量和端口电压的线性关系；
[0033]
等效电路模型获取单元，用于根据所述线性关系，得到直流元件的等效电导和两个电流源构成的等效电路，根据直流元件的等效电导和两个电流源构成的等效电路得到离散等效电路模型；
[0034]
状态空间方程获取单元，用于根据所述线性关系，得离散状态变量的迭代表达式，根据离散状态变量的迭代表达式得到离散状态空间方程。
[0035]
本发明实施例还提供一种计算机终端设备，包括一个或多个处理器和存储器。存储器与所述处理器耦接，用于存储一个或多个程序；当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如上述任一实施例所述的直流输电系统的离散特征值分析方法。
[0036]
本发明实施例还提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现如上述任一实施例所述的直流输电系统的离散特征值分析方法。
[0037]
与现有技术相比，本发明具有以下有益效果：
[0038]
本发明提供的直流输电系统的离散特征值分析方法及装置，适用于系统离散特征值分析的直流输电系统的离散状态空间模型和离散等效电路模型的建模方法，能够使已有的离散特征值分析方法的适用范围从纯交流系统延伸到了交直流混联系统中；
[0039]
本发明提供的直流输电系统的离散特征值分析方法及装置，可以嵌入到原有离散特征值分析程序中，采用与交流离散等效电路中相同的节点法消去中间变量，使得该过程具有明确的物理意义，仅需列写节点电压方程和注入电流方程即可，简单易行。
附图说明
[0040]
为了更清楚地说明本发明的技术方案，下面将对实施方式中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0041]
图1是本发明某一实施例提供的直流输电系统的结构示意图；
[0042]
图2是本发明某一实施例提供的直流输电系统的整流侧控制框图；
[0043]
图3是本发明某一实施例提供的直流输电系统的逆变侧控制框图；
[0044]
图4是本发明某一实施例提供的直流输电系统的离散特征值分析方法的流程示意图；
[0045]
图5是本发明另一实施例提供的直流输电系统的离散特征值分析方法的流程示意图；
[0046]
图6是本发明某一实施例提供的直流输电系统的直流输电元件电路图简化示意图；
[0047]
图7是本发明另一实施例提供的直流输电系统的直流输电元件离散等效电路示意图；
[0048]
图8是本发明某一实施例提供的直流输电系统的离散特征值分析装置的结构示意图；
[0049]
图9是本发明另一实施例提供的直流输电系统的离散特征值分析装置的结构示意图；
[0050]
图10是本发明某一实施例提供的计算机终端设备的结构示意图；
[0051]
图11是本发明某一实施例提供的直流输电系统的离散等效电路示意图。
具体实施方式
[0052]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0053]
应当理解，文中所使用的步骤编号仅是为了方便描述，不作为对步骤执行先后顺序的限定。
[0054]
应当理解，在本发明说明书中所使用的术语仅仅是出于描述特定实施例的目的而并不意在限制本发明。如在本发明说明书和所附权利要求书中所使用的那样，除非上下文清楚地指明其它情况，否则单数形式的“一”、“一个”及“该”意在包括复数形式。
[0055]
术语“包括”和“包含”指示所描述特征、整体、步骤、操作、元素和/或组件的存在，但并不排除一个或多个其它特征、整体、步骤、操作、元素、组件和/或其集合的存在或添加。
[0056]
术语“和/或”是指相关联列出的项中的一个或多个的任何组合以及所有可能组合，并且包括这些组合。
[0057]
为建立用于离散特征值分析的直流输电系统的离散状态空间模型及其离散等效电路模型，首先，需要对系统所含元件按类别建立连续状态空间模型，再做离散化处理，得到各类元件的离散状态空间模型和离散等效电路；在此基础上进行融合建模，构建全系统的离散状态空间模型和离散等效电路，基于后者消去状态空间方程的中间变量，得到离散状态变量的迭代表达式；随后便可提取离散状态矩阵，进行离散特征值分析，根据离散特征根位置判断系统稳定性。
[0058]
请参阅图1，图1是本发明某一实施例提供的直流输电系统的结构示意图。图中下标中dc表示直流输电元件，dq为整流侧或逆变侧锁相环对应的两相同步旋转坐标系，所标注的电气量为该坐标系下的量值，r和i分别代表整流侧和逆变侧，v代表阀侧，s代表系统
侧，t代表换流变。直流输电系统由直流输电元件、换流变及交流系统等效模型构成，其中直流输电元件包括整流侧和逆变侧的换流器及相关控制、直流支撑电容、直流电缆、换流器出口等效电感和电阻。整流侧和逆变侧均采用双闭环控制，其中整流侧外环采用定有功功率控制和定无功控制，逆变侧外环采用定直流电压控制和定无功功率控制，内环均为电流跟踪控制，整流侧控制框图如图2所示，逆变侧控制框图如图3所示。
[0059]
请参阅图4，图4是本发明某一实施例提供的直流输电系统的离散特征值分析方法的流程示意图。本实施例提供的直流输电系统的离散特征值分析方法包括以下步骤：
[0060]
s110，以整流侧公共连接点处的电压和逆变侧公共连接点处的电压为输入变量，以整流侧公共连接点处的电流和逆变侧公共连接点处的电流为输出变量，得到非线性动态方程；
[0061]
作为示例，结合图1进行说明，以整流侧公共连接点pcc1和逆变侧公共连接点pcc2处电压u
dcrdq
和u
dcidq
为直流输电元件的输入变量、电流i
dcrdq
和i
dcidq
作为输出变量，则可用一组非线性动态方程描述直流输电元件特性：
[0062][0063]
其中，x
dc
为连续状态变量构成的矩阵，u
dcrdq
为整流侧公共连接点处的电压，u
dcidq
为逆变侧公共连接点处的电压，i
dcrdq
整流侧公共连接点处的电流，i
dcidq
为逆变侧公共连接点处的电流，p为微分算子，t为方程所描述的时刻，f(x)表示px
dc
与x
dc
、u
dcrdq
、u
dcidq
之间的函数关系，g(x)表示i
dcrdq
、i
dcidq
与x
dc
、u
dcrdq
、u
dcidq
之间的函数关系。
[0064]
s120，对所述非线性动态方程做全微分线性化处理，得到稳态运行点处的小信号状态空间模型；其中，所述小信号状态空间模型包括直流输电元件的小信号连续状态空间方程、直流输电元件整流侧和逆变侧端口输出电流表达式；
[0065]
在一具体实施例中，对方程组(1)做全微分线性化处理，得到稳态运行点处的小信号状态空间模型：
[0066][0067]
其中，δ表示微增量，a
dc1
、b
dcr1
、b
dci1
、c
dc1
、d
dcr1
、d
dci1
、c
dc2
、d
dcr2
、d
dci2
为系数矩阵，pδx
dc
表示状态变量微增量的微分，δx
dc
为状态变量微增量，δu
dcrdq
为整流侧公共连接点处的电压微增量，δu
dcidq
为逆变侧公共连接点处的电压微增量，δi
dcrdq
为整流侧公共连接点处的电流微增量，δi
dcidq
为逆变侧公共连接点处的电流微增量。式(2)中第一式为直流输电元件的小信号连续状态空间方程，即状态变量微增量的微分pδx
dc
可以用状态变量微增量δx
dc
与端口输入电压微增量δu
dcrdq
、δu
dcidq
的线性组合表示；第二和第三式分别为直流输电元件整流侧和逆变侧端口输出电流表达式，即输出电流微增量δi
dcrdq
、δi
dcidq
可以用状态变量微增量δx
dc
和端口输入电压微增量δu
dcrdq
、δu
dcidq
的线性组合表示。
[0068]
s130，将所述小信号状态空间模型基于自身dq坐标系的端口量转换到系统基准两相同步旋转xy坐标系下，得到旋转坐标状态空间模型；
[0069]
在一具体实施例中，将方程(2)基于自身dq坐标系的端口量转换到系统基准两相
同步旋转xy坐标系下，需要说明的是：由于直流解耦，整流侧和逆变侧的两相同步旋转基准坐标系相互独立，建模者可自行规定xy坐标位置。采用q轴超前于d轴90
°
、y轴超前于x轴90
°
的惯例，记d轴领先x轴的电角度为γ，则有dq
‑
xy坐标变换矩阵：
[0070][0071]
采用该变换关系对端口量进行处理，将所得端口量转换关系式代入式(2)即可得到基于全系统基准两相同步旋转xy坐标系的小信号连续状态空间模型，即旋转坐标状态空间模型：
[0072][0073]
其中，a
dc2
、b
dcr2
、b
dci2
、d
dc3
、d
dcr3
、d
dci3
、c
dc4
、d
dcr4
、d
dci4
均为系数矩阵，δu
dcrdq
为整流侧公共连接点处的电压微增量，δu
dcidq
为逆变侧公共连接点处的电压微增量，δi
dcrdq
为整流侧公共连接点处的电流微增量，δi
dcidq
为逆变侧公共连接点处的电流微增量。
[0074]
s140，采用梯形积分规则对所述旋转坐标状态空间模型进行离散化处理，得到直流元件的离散状态空间方程和离散等效电路模型。
[0075]
在一具体实施例中，如图5所示，步骤s140，采用梯形积分规则对所述旋转坐标状态空间模型进行离散化处理，得到直流元件的离散状态空间方程和离散等效电路模型，包括：
[0076]
s141，将t时刻直流元件离散状态变量用前一时刻连续状态变量、整流侧端口输入电压和逆变侧端口输入电压线性表示，迭代得到t时刻端口电流、离散状态变量和端口电压的线性关系；
[0077]
由信号与系统相关知识可知，线性时不变(lti)系统可以用离散状态空间模型表示。对目标系统选定一组离散状态变量，则状态变量在离散时间点t处的值可以以状态变量和输入变量在(t
‑
δt)时刻值的线性组合表示，输出变量在t时刻的值可以用状态变量和输入变量在t时刻的值线性组合表示：
[0078][0079]
其中，h为离散状态变量，y
in
为输入变量，y
out
为输出变量。
[0080]
以下将给出直流输电元件离散状态空间模型的推导过程，基于离散状态空间模型可生成直流输电元件的离散等效电路。
[0081]
采用步长δt的梯形积分规则对小信号连续状态空间模型(式(4))做离散化处理，则对于(4)的第一式有：
[0082][0083]
定义离散状态变量h
dc
，将t时刻直流元件离散状态变量h
dc
用(t
‑
δt)时刻连续状态变量x
dc
、整流和逆变侧端口输入电压u
rxy
、u
ixy
线性表示：
[0084]
h
dc
(t)＝c
dc5
x
dc
(t
‑
δt) d
dcr5
u
dcrxy
(t
‑
δt) d
dci5
u
dcixy
(t
‑
δt)
ꢀꢀ
(7)
[0085]
其中，
[0086][0087]
再把式(7)代入式(6)，整理可得：
[0088]
x
dc
(t)＝a
dc3
h
dc
(t) b
dcr3
u
dcrxy
(t) b
dci3
u
dcixy
(t)
ꢀꢀ
(8)
[0089]
其中，a
dc3
＝i
n
(n为直流输电元件阶数,下同)，
[0090][0091]
将式(8)代入(2)的后两式(即直流输电元件端口输出电流表达式)，可得t时刻端口电流i
rxy
、i
ixy
和离散状态变量h
dc
、端口电压u
rxy
、u
ixy
的线性关系：
[0092][0093]
其中，c
dcdr
＝c
dc3
a
dc3
、d
dcdrr
＝c
dc3
b
dcr3
d
dcr3
、d
dcdri
＝c
dc3
b
dci3
d
dci3
、c
dcdi
＝
‑
c
dc4
a
dc3
、d
dcdir
＝c
dc4
b
dcr3
d
dcr4
、d
dcdii
＝c
dc4
b
dci3
d
dci4
，u
dcrdq
为整流侧公共连接点处的电压，u
dcidq
为逆变侧公共连接点处的电压，i
dcrdq
整流侧公共连接点处的电流，i
dcidq
为逆变侧公共连接点处的电流。
[0094]
s142，根据线性关系，得到直流元件的等效电导和两个电流源构成的等效电路，根据直流元件的等效电导和两个电流源构成的等效电路得到离散等效电路模型；
[0095]
根据线性关系，可将直流元件表示为由等效电导(包含了关联节点处的自电导和节点间的互电导)和两个电流源构成的等效电路如图6
‑
7所示，该电路为后续构建的全系统离散等效电路的重要组成部分。
[0096]
s143，根据线性关系，得离散状态变量的迭代表达式，根据离散状态变量的迭代表达式得到离散状态空间方程。
[0097]
整流侧节点自电导g
dcrrd
＝d
dcrrd
，逆变侧自电导g
dciid
＝
‑
d
dciid
，整流侧
‑
逆变侧互电导g
dcrid
＝
‑
d
dcrid
，逆变侧
‑
整流侧互电导g
dcird
＝d
dcird
，电流源值为c
dcdr
h
dc
和c
dcdi
h
dc
电流(当等效电路模型中所有量均为同一时刻值时，为方便表示，省略了时间点的标注，下同)。
[0098]
此外，联立(7)(9)两式可得离散状态变量h
dc
的迭代表达式，即直流元件的离散状态空间方程：
[0099][0100]
其中a
dc
＝c
dc5
a
dc3
，b
dcr
＝c
dc5
b
dcr3
d
dcr5
，b
dci
＝c
dc5
b
dci3
d
dci5
。
[0101]
式(10)、(9)共同构成了直流输电元件的离散状态空间模型：
[0102][0103]
请参阅图8，图8是本发明某一实施例提供的直流输电系统的离散特征值分析装置的结构示意图。本实施例提供的直流输电系统的离散特征值分析装置，所述直流输电系统，
包括直流输电元件、换流变和交流系统等效模型，所述直流输电元件包括换流器、直流支撑电容、直流电缆、换流器出口等效电感和电阻，其特征在于，所述离散特征值分析装置包括：
[0104]
电压分析模块210，用于以整流侧公共连接点处的电压和逆变侧公共连接点处的电压为输入变量，以整流侧公共连接点处的电流和逆变侧公共连接点处的电流为输出变量，得到非线性动态方程；
[0105]
线性化处理模块220，用于对所述非线性动态方程做全微分线性化处理，得到稳态运行点处的小信号状态空间模型；其中，所述小信号状态空间模型包括直流输电元件的小信号连续状态空间方程、直流输电元件整流侧和逆变侧端口输出电流表达式；
[0106]
坐标系旋转模块230，用于将所述小信号状态空间模型基于自身dq坐标系的端口量转换到系统基准两相同步旋转xy坐标系下，得到旋转坐标状态空间模型；
[0107]
离散化处理模块240，用于采用梯形积分规则对所述旋转坐标状态空间模型进行离散化处理，得到直流元件的离散状态空间方程和离散等效电路模型。
[0108]
在某一实施例中，如图9所示，离散化处理模块240还包括：
[0109]
线性关系分析单元241，用于将t时刻直流元件离散状态变量用前一时刻连续状态变量、整流侧端口输入电压和逆变侧端口输入电压线性表示，迭代得到t时刻端口电流、离散状态变量和端口电压的线性关系；
[0110]
等效电路模型获取单元242，用于根据所述线性关系，得到直流元件的等效电导和两个电流源构成的等效电路，根据直流元件的等效电导和两个电流源构成的等效电路得到离散等效电路模型；
[0111]
状态空间方程获取单元243，用于根据所述线性关系，得离散状态变量的迭代表达式，根据离散状态变量的迭代表达式得到离散状态空间方程。
[0112]
关于直流输电系统的离散特征值分析装置的具体限定可以参见上文中对于直流输电系统的离散特征值分析方法的限定，在此不再赘述。上述直流输电系统的离散特征值分析装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
[0113]
请参阅图10，本发明实施例提供一种计算机终端设备，包括一个或多个处理器和存储器。存储器与所述处理器耦接，用于存储一个或多个程序，当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如上述任意一个实施例中的直流输电系统的离散特征值分析方法。
[0114]
处理器用于控制该计算机终端设备的整体操作，以完成上述的直流输电系统的离散特征值分析方法的全部或部分步骤。存储器用于存储各种类型的数据以支持在该计算机终端设备的操作，这些数据例如可以包括用于在该计算机终端设备上操作的任何应用程序或方法的指令，以及应用程序相关的数据。该存储器可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，例如静态随机存取存储器(static random access memory，简称sram)，电可擦除可编程只读存储器(electrically erasable programmable read
‑
only memory，简称eeprom)，可擦除可编程只读存储器(erasable programmable read
‑
only memory，简称eprom)，可编程只读存储器(programmable read
‑
only memory，简称prom)，只读存储器(read
‑
only memory，简称rom)，磁存储器，快闪存储器，磁盘或光盘。
[0115]
在一示例性实施例中，计算机终端设备可以被一个或多个应用专用集成电路(application specific 1ntegrated circuit，简称as1c)、数字信号处理器(digital signal processor，简称dsp)、数字信号处理设备(digital signal processing device,简称dspd)、可编程逻辑器件(programmable logic device，简称pld)、现场可编程门阵列(field programmable gate array，简称fpga)、控制器、微控制器、微处理器或其他电子元件实现，用于执行上述的直流输电系统的离散特征值分析方法，并达到如上述方法一致的技术效果。
[0116]
在另一示例性实施例中，还提供了一种包括程序指令的计算机可读存储介质，该程序指令被处理器执行时实现上述任意一个实施例中的直流输电系统的离散特征值分析方法的步骤。例如，该计算机可读存储介质可以为上述包括程序指令的存储器，上述程序指令可由计算机终端设备的处理器执行以完成上述的直流输电系统的离散特征值分析方法，并达到如上述方法一致的技术效果。
[0117]
在某一具体实施例中，将各元件的离散等效电路按照系统拓扑连接得到直流输电系统的离散等效电路，如图11所示；另一方面联立各元件的离散状态空间模型构建全系统的离散状态空间模型，可以发现系统离散状态空间模型中的状态空间方程含有中间变量(元件输入电压)待消除。
[0118]
为消除中间变量，可采用与已有的交流系统的离散特征值分析类似的思路，即基于含直流输电的目标系统的离散等效电路，一方面列写其节点电压方程，得到节点注入电流i
xy
和元件输入电压u
xy
的关系式；另一方面根据系拓扑结构获得节点注入电流i与各元件离散状态变量h的关系式。然后将它们一起代入到系统的离散状态空间方程中，便可消去方程中的中间电压变量，得到系统全部离散状态变量h的差分迭代表达式，即系统不含中间变量的离散状态空间方程。
[0119]
在某一具体实施例中，直流输电系统的直流输电系统电气参数如表1所示，直流输电系统控制参数如表2所示。
[0120]
表1直流输电系统电气参数
[0121][0122]
表2直流输电系统控制参数
[0123][0124][0125]
提取离散状态空间方程中的离散状态矩阵后即可进行特征值求解，若求得的离散特征根全部位于单位圆内，则系统稳定；若有一个或一个以上离散特征根位于单位圆外，则系统不稳定，存在振荡风险。另外，根据梯形积分规则下的离散特征根与连续特征根的转换
关系，可以将离散状态空间特征值变换回连续状态空间特征值，从而获取模态的振荡频率等信息。
[0126]
解得图1系统的部分模态和对应离散特征值如表1第一、二列所示，可见模态λ
11,12
的离散特征根于单位圆外，由离散状态空间稳定性判据知该模态不稳定，将导致系统振荡发散。将离散特征值变换为连续特征值，如表1第三列所示，可见模态λ
11,12
特征值虚部对应振荡频率为11.2981hz。
[0127]
为验证上述结果，搭建该算例系统的pscad/emtdc电磁暂态仿真模型进行时域仿真验证，设定第5s后的系统参数与特征值分析的参数一致(表1～2)。时域仿真后观察逆变侧公共连接点pcc2有功波形(如图8所示)得到：系统在表1～2的参数下逐渐振荡发散，其振频率为11.11hz，与特征值求解结果11.2981hz一致，从而验证了系统离散状态空间模型的正确性，同时佐证了将直流输电系统离散等效电路用于离散特征值分析的可行性。
[0128]
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。