不依赖线路参数及边界元件的柔性直流线路方向纵联保护方法

1.本发明涉及电力系统继电保护领域，特别是涉及一种适用于柔性直流线路的新型方向纵联保护方法。


背景技术：

2.直流线路快速保护是柔性直流电网安全可靠运行的核心关键技术。目前，柔直线路保护原理一般以单端量保护为主保护，以纵联保护作为后备。其中，直流单端量保护利用线路两端安装的限流电抗器对故障行波的阻滞作用，可靠区分区内外故障，在几个毫秒之内即可超快速动作。纵联保护依赖双端信息判别区内外故障，由于通信的引入，动作时间在十几毫秒至几十毫秒等级，因此一般考虑作为直流线路后备保护。
3.根据工作原理的不同，直流线路纵联保护一般可以分为两种类型：一种类型基于线路模型(依赖线路参数)构建纵联判据，另一种类型则利用线路边界导致的边界特性构建纵联判据。常规直流输电线路后备保护时域电流差动保护属于典型的基于线路模型的纵联保护方法。为了消除分布电容电流影响，提高保护动作速度，文献《ultra-fast current differential protection with high-sensitivity for hvdc transmission lines》提出的基于贝瑞隆模型的快速电流差动保护利用精确线路模型计算差电流，能够自动免疫线路分布电容电流影响。此外，行波方向纵联保护利用正、反向行波幅值比判断故障方向，同样是利用了输电线路分布参数模型的行波传播特征。文献《apilotprotection scheme ofdc lines for multi-terminal hvdc grid》等提出的方向纵联保护方法则是利用线路边界元件(滤波器、电抗器)产生的边界特性判断故障方向，在此基础上构造方向纵联保护识别故障区段。
4.总结而言，现有的直流线路方向纵联保护主要存在两方面问题：1)基于线路模型的纵联保护方法需要提前获取精确的线路参数，否则将难以保证动作灵敏性和可靠性。由于不同工程的容量、电压等级不同，线路参数也会有所不同。因此，该类方法在不同工程实现时需要分别获取对应线路参数，对保护装置参数进行调整。显然，在工程实际中会大幅增加工程实施难度。2)基于线路边界元件的保护方法高度依赖线路两端装设的限流电抗器等边界元件。但是，随着直流电网拓扑结构的不断复杂化、规模的不断扩大，线路两端均存在边界元件的条件将无法保障。该类方法将难以再保证可靠动作。因此，有必要研究既不依赖线路参数、又不依赖边界元件的直流线路保护新原理，以降低保护的工程实现难度、提升工程普适性、应对线路边界缺失的场景。


技术实现要素：

5.针对目前直流线路纵联保护高度依赖线路参数或线路边界元件的问题，本发明提出一种不依赖线路参数及边界元件的柔性直流线路方向纵联保护方法，利用高低频测量阻抗幅值比判断故障方向，通过方向纵联实现区内外故障判别。
6.本发明提出了一种不依赖线路参数及边界元件的柔性直流线路方向纵联保护方法，利用高低频测量阻抗幅值比判断故障方向，实现区内故障及区外故障的判别；该方法具体步骤如下：
7.步骤1、设定保护的启动时刻为t
start
，确定采样时窗为(t
start-t1,t
start
+t2)；其中，t1为启动时刻以前的采样时长，由经验值给出；t2为启动时刻以后的采样时长，由故障电磁波在被保护线路上进行一次往返传播所需的时间与故障电磁波分别在保护背侧相邻线路1～线路n上进行一次往返传播所需的时间中的最小值来确定；
8.步骤2、获取保护安装处直流线路正极电压在步骤1确定的采样时窗内的采样数据获取保护安装处直流线路负极电压在步骤1确定的采样时窗内的采样数据在步骤1确定的采样时窗内的采样数据获取保护安装处直流线路正极电流在步骤1确定的采样时窗内的采样数据获取保护安装处直流线路负极电流在步骤1确定的采样时窗内的采样数据k-＝t1/ts，k
+
＝t2/ts，ts为采样周期；
9.步骤3、计算电压、电流的线模分量步骤3、计算电压、电流的线模分量计算公式如下所示：
[0010][0011]
其中，k-＝t1/ts，k
+
＝t2/ts，ts为采样周期，g为采样的数据点，k-为启动时刻以前的采样数据点，k
+
为启动时刻以后的采样数据点；
[0012]
步骤4、对线模电压、线模电流分别进行小波变换，获取第1层小波变换细节系数u
d_1
(k)和i
d_1
(k)，获取第j层小波变换细节系数u
d_j
(k)和i
d_j
(k)，其中k＝1,2,3,
…
,k-+k
+
+1；
[0013]
步骤5、计算高、低频线模测量阻抗幅值|z
m_high_f
|、|z
m_low_f
|，如下式所示：
[0014][0015]
步骤6、计算高、低频测量阻抗幅值比，如下式所示：
[0016][0017]
步骤7、判断故障方向：当kz》k
set
，判断故障为正方向故障，方向信号r
local
设为1；当kz≤k
set
，判断故障为反方向故障，r
local
设为0；
[0018]
步骤8、向对端发送方向信号r
local
，接收对端保护的方向判断信号r
opposite
；
[0019]
步骤9、判断故障区段：当r
local
&r
opposite
＝1，判断故障为区内故障；当r
local
&r
opposite
＝0，判断故障为区外故障。
[0020]
与现有技术相比，本发明能够达成以下技术效果：
[0021]
1)保护方法不依赖线路参数，在不同工程中进行应用时，无需根据线路参数的变化进行具体调整，因此能够大幅降低工程实施难度、提升工程普适性；
[0022]
2)保护方法不依赖线路边界元件，因此能够很好地适用于线路两端不存在边界元件的场景，很好地弥补了现有直流线路保护方法高度依赖线路边界元件的缺陷。
附图说明
[0023]
图1为柔性直流电网的典型系统结构图；
[0024]
图2为本发明中利用的柔直电网正、反向故障时测量阻抗幅值随频率变化的曲线图；
[0025]
图3为本发明提出的不依赖线路参数及边界元件的柔性直流线路方向纵联保护方法流程图。
具体实施方式
[0026]
下面结合附图及实施例对本发明作进一步详细说明。该实施例根据本发明技术方案进行实施，给出了详细具体的实施方式和操作过程，但本发明的应用范围不限于下述实施例。
[0027]
如图1所示，为柔性直流电网的典型系统结构图。图中以该四端柔性直流电网为例，line1～line4为直流输电线路，bus1～bus4为直流汇流母线，s1～s4为直流换流站。m、n分别代表线路line1两侧的保护。保护m处测得的正负极直流电压、电流分别为u
m_dcp
、u
m_dcn
、i
m_dcp
、i
m_dcn
，保护n处测得的正负极直流电压、电流分别为u
n_dcp
、u
n_dcn
、i
n_dcp
、i
n_dcn
。
[0028]
以保护m为例：
[0029]
当正方向发生故障时，线模电压、电流的测量阻抗表示为式(1)：
[0030][0031]
当反方向发生故障时，线模电压、电流的测量阻抗表示为式(2)：
[0032][0033]
其中，u
m_lm
(ω)代表保护m处线模(line mode)电压(由正负极电压计算得到)在某一频率ω(频段)下的分量、i
m_lm
(ω)代表保护m处线模电流(由正负极电流计算得到)在某一频率ω(频段)下的分量，z
c1
(ω)、z
c4
(ω)分别代表线路line1和line4的线模波阻抗。
[0034]
如图2所示，为正、反向故障测量阻抗幅值随频率变化的曲线图。根据上述公式(1)、(2)推导可得，当反方向发生故障以后，测量阻抗幅值在低频段(0～10hz左右)快速下降，随着频率的上升，测量阻抗快速稳定到一个固定值附近(非常轻微的下降)。整体而言，测量阻抗单调下降，最后近似稳定到一个固定值。因此，在反方向发生故障时，低频测量阻抗幅值大于高频测量阻抗幅值，即|z
m_high_f
|/|z
m_low_f
|≤1。
[0035]
当正方向发生故障以后，测量阻抗幅值在临界频率f
critical
之前逐渐下降，将临界频率f
critical
表示为式(3)：
[0036][0037]
而在临界频率之后，测量阻抗幅值则单调上升。因此在临界频率之后，|z
m_high_f
|/|z
m_low_f
|》1。
[0038]
根据上述特性差异，即可可靠判断故障方向。
[0039]
如图3所示，为本发明的一种不依赖线路参数及边界元件的柔性直流线路方向纵联保护方法整体流程图。该方法具体实施步骤如下：
[0040]
步骤1、设保护的启动时刻为t
start
，利用故障电磁波在被保护线路的传输时间、故障电磁波在保护背侧相邻线路1～线路n的传输时间中的最小值来确定采样时窗为(t
start-t1,t
start
+t2)。其中，t1为经验值，建议选取为1ms；t2＝min(2l
protected
/v、2l
backwar1
/v、2l
backwar2
/v、
…
2l
backwarn
/v)/k，k为可靠系数，一般取略大于1的数(例如1.2～1.3)，l
protected
为被保护线路长度，l
backward1
～l
backwardn
为保护背侧相邻线路1～线路n的线路长度，v为故障电磁波的传播速度。取故障电磁波在被保护线路上进行一次往返传播所需的时间2l
protected
/v、故障电磁波在分别在保护背侧线路1～n上进行一次往返传播所需的时间2l
backwar1
/v、2l
backwar2
/v、
…
2l
backwarn
/v的最小值，除以可靠系数，即为t2。以图2中保护m为例时，被保护线路为line1，背侧相邻线路为line4(该实施例中只有1条背侧相邻线路)；
[0041]
步骤2、获取保护安装处直流线路正极电压在步骤1确定的采样时窗内的采样数据获取保护安装处直流线路负极电压在步骤1确定的采样时窗内的采样数据在步骤1确定的采样时窗内的采样数据获取保护安装处直流线路正极电流在步骤1确定的采样时窗内的采样数据获取保护安装处直流线路负极电流在步骤1确定的采样时窗内的采样数据k-＝t1/ts，k
+
＝t2/ts，ts为采样周期；
[0042]
步骤3、计算电压、电流的线模分量步骤3、计算电压、电流的线模分量步骤3、计算电压、电流的线模分量计算公式如下所示：
[0043][0044]
其中，k-＝t1/ts，k
+
＝t2/ts，ts为采样周期，g为采样的数据点，k-为启动时刻以前的采样数据点，k
+
为启动时刻以后的采样数据点；
[0045]
步骤4、对线模电压、线模电流分别进行小波变换，获取第1层小波变换细节系数u
d_1
(k)和i
d_1
(k)；获取第j层小波变换细节系数u
d_j
(k)和i
d_j
(k)，其中k＝1,2,3,
…
,k-+k
+
+1，变量j的取值原则为确保式(5)成立的最大值：
[0046][0047]
其中，fs代表保护装置采样率，l
arm
代表保护所在换流站mmc桥臂电抗电感值，c
sm
代表子模块电容值，nv代表mmc每个桥臂的子模块数量；
[0048]
步骤5、计算高、低频线模测量阻抗幅值|z
m_high_f
|、|z
m_low_f
|，如式(6)所示：
[0049][0050]
步骤6、计算高、低频测量阻抗幅值比，如式(7)所示：
[0051][0052]
步骤7、判断故障方向：当kz》k
set
，判断故障为正方向故障，方向信号r
local
设为1；当kz≤k
set
，判断故障为反方向故障，r
local
设为0；k
set
为大于1的整定值，建议取1.5～2；
[0053]
步骤8、向对端发送方向信号r
local
，接收对端保护的方向判断信号r
opposite
；
[0054]
步骤9、判断故障区段：当r
local
&r
opposite
＝1，判断故障为区内故障；当r
local
&r
opposite
＝0，判断故障为区外故障。