一种基于随机矩阵的柔性直流电网线路纵联保护方法与流程

1.本发明属于直流电网继电保护技术领域，尤其涉及一种基于随机矩阵的柔性直流电网线路纵联保护方法。


背景技术：

2.基于模块化多电平换流器的多端直流电网具有谐波含量低、没有换相失败问题和能够独立控制有功功率和无功功率等优点，被认为是实现未来能源输送的理想方式。多端直流电网的直流线路故障将产生巨大的故障电流，使得多端直流电网的换流器在极短时间内闭锁、停运，影响整个电网的供电可靠性。快速准确识别出故障直流线路并切除是保证多端直流电网安全稳定运行的基础。已有的多端直流电网保护方法大多利用电压变化率或数字信号处理方法进行设计，当发生高阻故障时，可能不能正确识别出故障线路，并且当采样数据受到噪声和异常数据干扰时，可能会使方法输出错误的故障识别结果。


技术实现要素：

3.针对现有技术的不足，为解决现有多端直流电网保护方法不能准确识别高阻故障，且识别结果易受噪声和异常数据影响的问题。本发明提供一种基于随机矩阵的柔性直流电网线路纵联保护方法。
4.本发明的一种基于随机矩阵的柔性直流电网线路纵联保护方法，包括以下步骤：
5.步骤1：获取直流线路、直流母线和相邻直流线路电压的n个采样值，根据采样值计算1模电压故障分量和0模电压反行波，将直流线路和直流母线的1模电压故障分量、直流母线和相邻直流线路1模电压故障分量组成2个2
×
n的矩阵，对两个矩阵分别进行复制、平移和镜像处理，扩展成为两个原始矩阵。
6.步骤2：对直流线路的0模电压反行波积分，根据积分值大小判定故障类型，再根据故障类型和系统拓扑参数求取噪声矩阵，并分别与原始矩阵叠加形成2个随机矩阵，根据随机矩阵理论，将2个随机矩阵转换为2个标准矩阵积。
7.步骤3：构造柔性直流电网线路的故障方向判据。
8.确定阈值，计算2个标准矩阵积的平均谱半径，如果2个平均谱半径同时小于给定阈值，为正向故障，输出本地检测信号tpm＝1，并向线路对侧保护装置发送。
9.步骤4：构造柔性直流电网的故障线路判据。
10.等待对侧检测信号tpn，如果两端检测信号tpm和tpn都为1，则该线路为故障线路。
11.步骤1具体为：
12.和分别为直流线路、直流母线和相邻直流直线路1模电压故障分量。
13.将组成2个2
×
tn的矩阵和
14.[0015][0016]
其中，k表示当前采样时刻，和分别代表直流线路、直流母线和相邻直流线路当前时刻的1模电压故障分量采样值，tn表示时间窗宽度，即为步骤1所述的n，k-tn 1，k-tn 2也分别为采样时刻。
[0017]
柔性直流线路要求其保护方法的识别速度快，其保护算法的时间窗非常短，但随机矩阵理论要求原始数据矩阵具有一定的维数。因此将和各行复制w-1次并向下平移后，再镜像一次，形成2w
×
2tn的原始矩阵和
[0018][0019][0020]
步骤2具体为：
[0021]
将直流线路的0模电压反行波积分，根据积分值大小判定故障类型，判据如下：
[0022][0023]
其中，为直流线路的0模电压行波，k-i为采样时刻，fn为采样频率，s
set
为判定阈值。
[0024]
根据故障类型选择相应的1模电压初始行波幅值柔性直流电网在对称双极接线方式下，不同故障类型的如下：
[0025][0026]
其中，u
dc
为柔性直流电网直流侧的额定电压，和分别为直流线路的1模和0模波阻抗，rf为过渡电阻。
[0027]
根据电路理论，和与在复频域有如下幅值关系：
[0028][0029]
其中，s为复频率，和为复频域下直流线路、母线和相邻线路1模电压幅值，为复频域下的1模电压初始反行波幅值，l
dc
为直流线路两侧的限流电抗器电抗值，z
mmc1
为换流器1在复频域下等效阻抗，由下式求得：
[0030][0031]
其中，z
mmc
为换流器复频域下的等效阻抗，l
arm
和c
sub
为换流器桥臂电抗值和子模块电容值，h为桥臂上子模块的个数。
[0032]
定义复频域下的系数：
[0033][0033][0033][0034]
其中，γ
α
(s)、γ
β
(s)和γ
γ
(s)分别为复频域下直流线路、母线和相邻线路1模电压故障分量和1模电压初始反行波的幅值比；对γ
α
(s)、γ
β
(s)和γ
γ
(s)进行逆拉普拉斯变换，得到γ
α
(t)、γ
β
(t)和γ
γ
(t)，定义0时刻γ
β
(0)和γ
γ
(0)的差值的0.5倍为最大差异系数；则噪声矩阵的幅值am为：
[0035][0036]
其中，kf为故障系数；建议kf∈(0.6,0.8)。
[0037]
构造与原始矩阵维数相同的噪声矩阵，并与原始矩阵叠加，形成每个矩阵元素均为随机数的随机矩阵x
loc
和x
adj
：
[0038][0039][0040]
其中，e
gas
为2w
×
2tn的高斯白噪声矩阵，其每个元素为服从标准正态分布的随机数。
[0041]
将随机矩阵x
loc
和x
adj
转换为标准矩阵积：
[0042]
首先对x
loc
和x
adj
中的每个元素按下式进行处理得到归一化矩阵和
[0043][0044]
其中，x
i,j
和分别为随机矩阵和归一化矩阵的每个元素，下标i和j分别代表矩阵的行和列，μ(xj)和σ(xj)分别为随机矩阵第j列元素的期望和方差，和为归一化
矩阵第j列元素的期望和方差；归一化矩阵的元素满足期望和i＝1,
···
2w；j＝1,
···
2tn。
[0045]
然后求取归一化矩阵和的奇异值等价矩阵和
[0046][0047][0048]
其中，u为2w
×
2w的哈尔酉矩阵，上标h代表矩阵的转置；按下式直接对奇异值等价矩阵进行归一化处理得到标准矩阵积和
[0049][0050]
其中，为标准矩阵积的第i行数据，为奇异值等价矩阵的第i行数据，为奇异值等价矩阵第i行数据的方差。
[0051]
步骤3具体为：
[0052]
根据非渐进随机矩阵理论中的单环定理，当标准矩阵积的元素为期望μ＝0、方差σ2＝1的独立同分布高斯随机变量时，如果的行数m、列数n足够大时，且保持行列比不变时，的特征值的经验谱分布满足概率密度函数f(λ)：
[0053][0054]
其中，π为圆周率，d为原始随机矩阵的个数，c为标准矩阵积的行列比，即c＝m/n,c∈(0,1]，m和n分别的行数和列数，λ为标准矩阵积的特征值，|λ|为特征值的模值。
[0055]
故障直流线路的1模电压故障分量远大于直流母线和相邻线路的1模电压故障分量，由故障直流线路1模电压故障分量构造的两个标准矩阵积不符合单环定理中各元素满足独立同分布的高斯随机变量这一条件，的特征值的经验谱分布将不满足概率密度函数f(λ)。本发明的故障方向判据阈值r
set
由概率密度函数f(λ)的下界(1-c)
0.5d
确定：
[0056]rset
＝k
rel
(1-c)
0.5d
[0057]
其中，k
rel
为一个小于1的可靠系数，本发明的k
rel
＝0.95。d的取值为1。
[0058]
求取标准矩阵积和的所有特征值，由下式计算平均谱半径r
loc
和r
adj
：
[0059][0060]
其中，|λi|表示第i个特征值的模值，c为所有特征值的个数。
[0061]
柔性直流电网线路的故障方向判据：如果平均谱半径r
loc
和r
adj
同时小于r
set
，则为正向故障，本地检测信号tpm＝1，并向线路对侧保护发送：
[0062]
[0063]
步骤4具体为：
[0064]
柔性直流电网的故障线路判据：线路对侧保护采用同样的方式计算了检测信号tpn，如果tpm和tpn都为1，判定直流线路为故障线路：
[0065][0066]
根据步骤2得到故障类型，对相对应故障极线路发送跳闸命令。
[0067]
本发明的有益技术效果为：
[0068]
1)本发明利用直流线路和直流母线、直流母线和相邻直流线路的1模电压故障分量组成原始矩阵，叠加由故障类型和系统拓扑参数求取的噪声矩阵后，通过矩阵变换得到了两个标准矩阵积，将两个标准矩阵积的平均谱半径与阈值比较来判断故障方向，并向线路对侧发送检测信号。本发明的故障方向信息为布尔量，对数据的同步性要求不高，能在发生高阻故障时准确判断故障方向。
[0069]
2)本发明的纵联保护方法不受故障位置、系统拓扑结构变化等影响，当过渡电阻为600欧姆时也能准确识别故障。
[0070]
3)本发明的纵联保护方法，利用数据窗内采样数据的统计信息进行故障识别，具有良好的抗异常数据能力，当采样数据受到行业标准规定的最大噪声或占20％数据窗宽度的异常数据干扰下，本发明所提的纵联保护方法仍能准确识别故障线路。
[0071]
4)当系统参数出现10％的误差，本发明仍能准确识别故障线路，本发明还可应用于不同接线方式下的多端直流电网。
附图说明
[0072]
图1为本发明基于随机矩阵的柔性直流电网线路纵联保护方法的流程图。
[0073]
图2为四端直流电网测试系统。
[0074]
图3为线路l1的1模故障附加网络。
[0075]
图4为γ
α
(t)、γ
β
(t)和γ
γ
(t)的变化图。
[0076]
图5为直流线路杆塔图。
[0077]
图6为f1发生经300ω正极接地故障时的平均谱半径。
[0078]
图7为发生区外故障(f4处双极故障)时的平均谱半径。
[0079]
图8为发生区外故障(f6处双极故障)时的平均谱半径。
[0080]
图9为收到异常数据干扰时的电压波形图和平均谱半径。
[0081]
图10为30db干扰下的电压波形图。
[0082]
图11为30db干扰下的平均谱半径。
[0083]
图12为系统参数存在误差时的平均谱半径。
具体实施方式
[0084]
下面结合附图和具体实施例对本发明做进一步详细说明。
[0085]
本发明所提一种基于随机矩阵的柔性直流电网线路纵联保护方法如图1所示，首先获取直流线路、直流母线和相邻直流线路电压的n个采样值，根据采样值计算其1模电压
故障分量和0模电压反行波，将直流线路和直流母线的1模电压故障分量、直流母线和相邻直流线路1模电压故障分量组成2个2
×
n的矩阵，对两个矩阵分别进行复制、平移和镜像处理，扩展成为两个原始矩阵；对直流线路的0模电压反行波积分，根据积分值大小判定故障类型，然后根据故障类型和系统拓扑参数计算噪声矩阵，并分别与原始矩阵叠加形成2个随机矩阵，根据随机矩阵理论，将2个随机矩阵转换为2个标准矩阵积；确定阈值，计算2个标准矩阵积的平均谱半径，如果2个平均谱半径同时小于给定阈值，为正向故障，输出本地检测信号tpm＝1，并向线路对侧保护装置发送；同时等待对侧检测信号tpn，如果两端检测信号tpm和tpn都为1，则该线路为故障线路。具体包括如下步骤：
[0086]
1、对于如图2所示的四端直流电网，其中，mmc1～mmc4代表四个换流站；l1，l2，l3和l4代表四条直流线路，l
13
代表没投入的备用直流线路；b1，b2，b3和b4分别代表4个直流母线；l
dc
为限流电抗器，在实际直流输电工程中通常设置为100-200mh，本发明算例中为150mh；p
ij
为保护装置与直流断路器，下标i,j为线路连接换流站的编号；f1～f8为本发明算例中设置的故障点。
[0087]
以线路l1的p
12
为例阐述2个原始矩阵的形成方法，如图3所示为线路l1故障时的1模故障附加网络。其中m，b1和l分别代表直流线路，直流母线和相邻直流线路的测量点，f为故障点，为1模附加电压分量，rf为过渡电阻；为故障点的1模电压，z
mmc1
和z
mmc2
为换流站1和2的等效阻抗。
[0088]
选择常用的直流电网故障启动判据——电压梯度算法作为本发明的启动判据，电压梯度算法定义如下：
[0089][0090]
其中，为直流线路的1模电压故障分量，k-j为采样值时刻，为第k个采样值的电压梯度，为对电压梯度取模值；为启动阈值，本发明当检测到电网发生故障时，将直流线路和直流母线的1模电压故障分量和直流母线和相邻直流线路1模电压故障分量和组成2个2
×
tn的矩阵：
[0091][0092][0093]
其中，k表示当前采样时刻，和分别代表直流线路、直流母线和相邻直流线路当前时刻的1模电压故障分量采样值，tn表示时间窗宽度，即为本发明权利要求书中的n，k-tn 1,k-tn 2也分别为采样时刻。
[0094]
柔性直流线路要求其保护方法的识别速度快，其保护算法的时间窗非常短，但随机矩阵理论要求原始数据矩阵具有一定的维数。因此将和各行复制w-1次并向下平
移后，再镜像一次，形成2w
×
2tn的原始矩阵和
[0095][0096][0097]
2、将直流线路的0模电压反行波积分，根据积分值大小判定故障类型，判据如下：
[0098][0099]
其中，为直流线路的0模电压行波，k-i为采样时刻，fn为采样频率，s
set
为判定阈值，本发明s
set
＝5；根据故障类型选择相应的1模电压初始行波幅值柔性直流电网在对称双极接线方式下，不同故障类型的如下：
[0100][0101]
其中，u
dc
为柔性直流电网直流侧的额定电压，和分别为直流线路的1模和0模波阻抗，rf为过渡电阻。根据电路理论，和与在复频域有如下幅值关系：
[0102][0103]
其中，s为复频率，和为复频域下直流线路、母线和相邻线路1模电压幅值，为复频域下的1模电压初始反行波幅值，l
dc
为直流线路两侧的限流电抗器
电抗值，z
mmc1
为换流器1在复频域下的等效阻抗，可由下式求得：
[0104][0105]
其中，z
mmc
为换流器复频域下的等效阻抗，l
arm
和c
sub
为换流器桥臂电抗值和子模块电容值，h为桥臂上子模块的个数。定义复频域下的系数：
[0106][0106][0106][0107]
其中，γ
α
(s)、γ
β
(s)和γ
γ
(s)分别表示复频域下直流线路、母线和相邻线路1模电压和1模电压初始反行波的幅值比。对γ
α
(s)、γ
β
(s)和γ
γ
(s)进行逆拉普拉斯变换，可以得到γ
α
(t)、γ
β
(t)和γ
γ
(t)，其变化趋势图如附图4所示。定义0时刻γ
β
(0)和γ
γ
(0)的差值的0.5倍为最大差异系数。则噪声矩阵的幅值am为：
[0108][0109]
其中，kf为故障系数，本发明建议kf∈(0.6,0.8)。构造与原始数据矩阵维数相同的噪声矩阵，并与原始矩阵叠加，形成每个矩阵元素均为随机数的随机矩阵x
loc
和x
adj
[0110][0111][0112]
其中，e
gas
为2w
×
2tn的高斯白噪声矩阵，其每个元素为服从标准正态分布的随机数。
[0113]
由随机矩阵理论可知，可以将随机矩阵x
loc
和x
adj
按如下步骤转换为标准矩阵积。首先对x
loc
和x
adj
中的每个元素按下式进行处理得到归一化矩阵和
[0114][0115]
其中，x
i,j
和分别为随机矩阵和归一化矩阵的每个元素，下标i和j分别代表矩阵的行和列，μ(xj)和σ(xj)分别为随机矩阵第j列元素的期望和方差，和为归一化矩阵第j列元素的期望和方差。归一化矩阵的元素满足期望和i＝1,
···
2w；j＝1,
···
2tn。然后求取归一化矩阵和的奇异值等价矩阵和
[0116][0117][0118]
其中，u为2w
×
2w的哈尔酉矩阵，上标h代表矩阵的转置。本发明求得的奇异值等价矩阵和中只包含1个随机矩阵，所以按下式直接对奇异值等价矩阵进行归一化处理得到标准矩阵积和
[0119][0120]
其中，为标准矩阵积的第i行数据，为奇异值等价矩阵的第i行数据，为奇异值等价矩阵第i行数据的方差。
[0121]
3、根据非渐进随机矩阵理论中的单环定理，当标准矩阵积的元素为期望μ＝0、方差σ2＝1的独立同分布高斯随机变量，如果的行数m、列数n足够大时，且保持行列比不变时，的特征值的经验谱分布满足概率密度函数f(λ)：
[0122][0123]
其中，π为圆周率，d为原始随机矩阵的个数，c为标准矩阵积的行列比，即c＝m/n,c∈(0,1]，m和n分别的行数和列数，λ为标准矩阵积的特征值，|λ|为特征值的模值。
[0124]
故障直流线路的1模电压故障分量远大于直流母线和相邻线路的1模电压故障分量，由故障直流线路1模电压故障分量构造的两个标准矩阵积不符合单环定理中各元素为独立同分布的高斯随机变量这一条件，的特征值的经验谱分布将不满足概率密度函数f(λ)。故本发明故障方向判据阈值r
set
由概率密度函数f(λ)的下界(1-c)
0.5d
确定：
[0125]rset
＝k
rel
(1-c)
0.5d
[0126]
其中，k
rel
为可一个小于1的可靠系数，本发明的k
rel
＝0.95。本发明求得的奇异值等价矩阵和中只包含1个随机矩阵，故d＝1。
[0127]
求取标准矩阵积和的所有特征值λ，按下式计算平均谱半径(msr)r
loc
和r
adj
：
[0128][0129]
其中，r为平均谱半径，|λi|表示第i个特征值的模值，c为所有特征值的个数。
[0130]
柔性直流电网线路的故障方向判据：如果平均谱半径r
loc
和r
adj
同时小于r
set
，则为正向故障，本地检测信号tpm＝1，并向线路对侧保护发送：
[0131][0132]
4、柔性直流电网的故障线路判据：线路对侧采用同样方式计算了检测信号tpn，如果tpm和tpn都为1，判定直流线路为故障线路。
[0133][0134]
根据步骤2判断的故障类型，向对应故障极线路发送跳闸命令。
[0135]
实施例
[0136]
以下是利用电磁暂态软件pscad/emtdc搭建四端直流电网采集实验数据进行的实
验算例，电网拓扑结构图如图2所示。
[0137]
其中，换流站子模块为半桥型，换流站mmc1采用定直流电压控制，换流站mmc2～mmc4采用定有功功率控制，测试系统的主要参数如表1所示；直流线路采用依频模型，其参数和杆塔结构如表2和附图5所示；系统采样频率fn＝50khz，时间窗为1ms，结合采样率和时间窗可得tn＝50；复制次数w＝25，故原始矩阵和的维数为50
×
100。将行列比50/100带入阈值计算公式，可得阈值r
set
＝k
rel
(1-50/100)
0.5
＝0.672。本实施例中所有算例以线路l1两侧的保护装置p
12
和p
21
为例，所有故障时间均设置在2.3s。其中，r
loc_12
和r
adj_12
、r
loc_21
和r
adj_21
分别代表保护装置p
12
、p
21
计算的两个平均谱半径；s
12
和s
21
代表p
12
和p
21
计算的电压反行波积分值，由式计算；tp
12
和tp
21
为两侧保护算法得到的检测信号。
[0138]
表1 四端直流电网参数
[0139][0140][0141]
表2 直流线路参数
[0142][0143]
算例1：
[0144]
为验证本发明对区内各种类型故障的识别效果，在附图2中l1的5％处(f1)，50％处(f2)和95％处(f3)分别设置正极金属性接地故障(pgf)和双极故障(ppf)，p
12
和p
21
的识别结果如表3所示。
[0145]
表3 区内故障识别结果
[0146][0147]
由表3中数据可知，两侧保护装置的两个平均谱半径r
loc_12
和r
adj_12
、r
loc_21
和r
adj_21
均同时小于阈值r
set
＝0.672，本发明可以准确识别出区内不同故障位置，不同故障类型的故障。
[0148]
算例2：
[0149]
为验证本发明识别高阻故障的能力，当f2发生经300ω正极接地故障时，保护装置p
12
和p
21
计算的平均谱半径r
loc_12
和r
adj_12
、r
loc_21
和r
adj_21
的变化如附图6所示。由附图6可知，本发明能在故障发生后短时间内快速输出检测信号。在附图2中l1的5％处(f1)，50％处(f2)和95％处(f3)分别设置经300ω和600ω接地的正极故障，识别结果表4所示。
[0150]
表4 经不同过渡电阻正极故障
[0151][0152][0153]
由表4中数据可知，相比于正极金属性接地故障，经过渡电阻接地时，两侧保护装置的两个平均谱半径略有上升，但r
loc_12
和r
adj_12
、r
loc_21
和r
adj_21
均同时小于阈值r
set
＝0.672，本发明所提方法可以识别过渡电阻高达600ω的区内故障。
[0154]
算例3：
[0155]
纵联保护方法要求方法不受各种类型区外故障的影响，为验证本发明在区外故障发生时能可靠不误动，在f4、f5和f6处设置双极故障。同时为了验证本发明在换流器交流侧故障时也能可靠不误动，在f7处设置三相接地故障。仿真结果如表5所示。
[0156]
表5 区外故障识别结果
[0157][0158]
由表中数据可知，两侧保护装置计算的两个平均谱半径(msr)不会同时小于给定阈值，在区外故障时可靠不动作。f4和f6故障时，平均谱半径的变化附图7和附图8所示。由附
图7、附图8可知，当发生区外故障时，本发明所提出的故障特征量r
loc_12
、r
adj_12
、r
loc_21
和r
adj_21
不会同时低于给定阈值，可以保证所提方法在区外故障时可靠不误动。
[0159]
算例4：
[0160]
为验证本发明所提方法在受到异常数据和噪声干扰时的性能，以保护装置p
12
为例，对1模电压故障分量测量值在2.29s和2.3006s添加异常数据，异常数据幅值为直流电压额定值的1倍标幺值。添加异常数据个数为1和10，分别占整个数据窗宽度的2％和20％。当f2处发生正极接地时且直流线路1模电压故障分量存在20％数据窗宽度的噪声干扰时，的波形和p
12
计算的两个平均谱半径变化如附图9所示，可知，当受到20％数据窗宽度的异常数据干扰时，本发明所提方法也不会对非故障线路输出错误识别结果，并且仍能正确识别故障线路。
[0161]
当受到不同数据窗宽度的异常数据干扰时，不同故障类型不同故障位置的识别结果如表6所示，由表6中数据可知，在异常数据干扰下，本发明能可靠正确识别出故障线路。
[0162]
表6 异常数据干扰时的识别结果
[0163][0164]
根据行业标准，采样装置的最低信噪比(snr)应大于30db。为验证本发明所提方法在受到噪声干扰时的性能，当f1点发生金属性接地故障，对1模电压故障分量叠加30db的白噪声。叠加噪声后的直流线路，直流母线和相邻直流线路的1模电压故障分量如附图10所示，两侧保护装置计算的平均谱半径(msr)如附图11所示。从附图11可得，相比于未叠加30db噪声的附图6，当叠加了30db白噪声后，正常状态的msr较未受30db噪声干扰时略有下降，从0.84变为了0.77，但仍大于所给阈值，即当受到噪声干扰时，不会向正常直流线路发出错误检测信号。由附图11可知，在数据受到30db噪声干扰时，本发明所提方法仍能可靠识别出故障线路。
[0165]
算例5：
[0166]
为验证当系统参数存在误差时，本发明所提方法的性能，当f1和f2点发生金属性接地故障时，对系统参数叠加了
±
10％的标准正态分布随机误差，进行了1000次仿真实验。实验结果如附图12所示。图中的每一个点代表p
12
和p
21
计算的一个平均谱半径值(msr)，图中右边柱状图代表了每一点落在图中区域的概率。可以看出，当f1和f2点发生金属性接地故障，p
12
和p
21
计算的每一个平均谱半径都小于给定阈值，说明即便系统参数存在误差时，本发明所提方法仍能正确识别故障线路。
[0167]
算例6：
[0168]
为了验证当系统拓扑结构发生变化时本发明所提方法的性能，设置了如下实验。
[0169]
工况1：当一条直流线路因为检修而退出运行时。具体为通过直流断路器使l4线路退出运行，并在此时设置f1和f2发生金属性接地故障。
[0170]
工况2：当一换流站因为检修而退出运行时。具体为当换流站mmc4因检修而退出运行，通过备用直流线路l
13
连接换流站mmc1和mmc3，并在此时设置f1和f2发生金属性接地故障。工况1和2的实验结果如表7所示。
[0171]
表7 不同工况下的实验结果
[0172][0173]
由表7中数据可知，当四端直流电网在不同工况运行下，p
12
和p
21
的平均谱半径值(msr)r
loc_12
和r
adj_12
、r
loc_21
和r
adj_21
均小于给定阈值，本发明所提方法能在不同工况下正确识别故障线路。
[0174]
算例7
[0175]
为验证本发明所提方法在对称单极系统下的性能，将附图2所示的对称双极四端电网改造为对称单极的四端电网，具体为：除直流电压额定值和功率额定值外，系统和线路参数不变。功率额定值变为1500mva，直流电压额定值变为500kv。换流器交流侧变压器采用yn/d接法，且二次侧经高阻接地。
[0176]
当对称单极系统发生直流线路单极接地故障时，故障电流回路与对称双极不同，主要区别在于换流器等效阻抗不同，此时换流器等效阻抗变为当本发明所提方法应用于对称单极系统单极故障时，需要将等效阻抗修改后带入，并重新计算最大差异系数。当对称单极系统发生直流线路双极故障时，故障电流回路与对称双极系统的单极故障相同，换流器等效阻抗不变。需要说明的是，在对称单极系统中，故障极识别判据和启动判据不变。此时在f1和f2点分别设置金属性接地故障和双极故障，仿真结果表8所示。
[0177]
表8 对称单极系统仿真结果
[0178][0179]
由表8中数据可知，当本发明所提方法应用于对称单极系统时，仍可利用0模电压反行波判定故障极，再根据故障极和最大差异系数求取相应的噪声矩阵，进行标准矩阵积的构造，最后仍利用两个标准矩阵积的平均谱半径识别故障方向，并根据双端的检测信号识别故障线路。拓扑结构变化时，本发明只需稍加改变参数数值便可应用于不同接线方式
的直流电网。