CN110649577B - 一种基于纯电流特征的直流输电线路保护方法和*** - Google Patents

一种基于纯电流特征的直流输电线路保护方法和*** Download PDF

Info

Publication number
CN110649577B
CN110649577B CN201910599291.3A CN201910599291A CN110649577B CN 110649577 B CN110649577 B CN 110649577B CN 201910599291 A CN201910599291 A CN 201910599291A CN 110649577 B CN110649577 B CN 110649577B
Authority
CN
China
Prior art keywords
current
transmission line
direct current
current transmission
criterion
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910599291.3A
Other languages
English (en)
Other versions
CN110649577A (zh
Inventor
周泽昕
柳焕章
吕鹏飞
曹虹
王兴国
阮思烨
杜丁香
李勇
杨国生
金明亮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI filed Critical State Grid Corp of China SGCC
Priority to CN201910599291.3A priority Critical patent/CN110649577B/zh
Publication of CN110649577A publication Critical patent/CN110649577A/zh
Priority to PCT/CN2020/096763 priority patent/WO2021000738A1/zh
Priority to US17/596,826 priority patent/US11973340B2/en
Application granted granted Critical
Publication of CN110649577B publication Critical patent/CN110649577B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02H—EMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H7/00—Emergency protective circuit arrangements specially adapted for specific types of electric machines or apparatus or for sectionalised protection of cable or line systems, and effecting automatic switching in the event of an undesired change from normal working conditions
    • H02H7/26—Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured
    • H02H7/268—Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured for dc systems
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
    • G01R19/165—Indicating that current or voltage is either above or below a predetermined value or within or outside a predetermined range of values
    • G01R19/16533—Indicating that current or voltage is either above or below a predetermined value or within or outside a predetermined range of values characterised by the application
    • G01R19/16538—Indicating that current or voltage is either above or below a predetermined value or within or outside a predetermined range of values characterised by the application in AC or DC supplies
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02H—EMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H1/00—Details of emergency protective circuit arrangements
    • H02H1/0007—Details of emergency protective circuit arrangements concerning the detecting means
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02H—EMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H1/00—Details of emergency protective circuit arrangements
    • H02H1/0092—Details of emergency protective circuit arrangements concerning the data processing means, e.g. expert systems, neural networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02H—EMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H3/00—Emergency protective circuit arrangements for automatic disconnection directly responsive to an undesired change from normal electric working condition with or without subsequent reconnection ; integrated protection
    • H02H3/26—Emergency protective circuit arrangements for automatic disconnection directly responsive to an undesired change from normal electric working condition with or without subsequent reconnection ; integrated protection responsive to difference between voltages or between currents; responsive to phase angle between voltages or between currents
    • H02H3/28—Emergency protective circuit arrangements for automatic disconnection directly responsive to an undesired change from normal electric working condition with or without subsequent reconnection ; integrated protection responsive to difference between voltages or between currents; responsive to phase angle between voltages or between currents involving comparison of the voltage or current values at two spaced portions of a single system, e.g. at opposite ends of one line, at input and output of apparatus
    • H02H3/286—Emergency protective circuit arrangements for automatic disconnection directly responsive to an undesired change from normal electric working condition with or without subsequent reconnection ; integrated protection responsive to difference between voltages or between currents; responsive to phase angle between voltages or between currents involving comparison of the voltage or current values at two spaced portions of a single system, e.g. at opposite ends of one line, at input and output of apparatus involving comparison of similar homopolar quantities

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Emergency Protection Circuit Devices (AREA)

Abstract

本发明提供一种基于纯电流特征的直流输电线路保护方法和***。所述方法和***通过采集直流输电线路电流确定线路的电流差分值和电流突变量，并根据所述线路电流差分值和电流突变量确定预先设置的保护判据是否成立，来确定直流保护动作是否出口。本发明所述的基于纯电流特征的直流输电线路保护方法和***能够进一步加快金属性故障保护动作速度，大幅提升高阻故障保护动作速度，雷击或异常大数时保护可靠不动作，单极接地故障工况下非故障极线路保护可靠不误动，全面提升直流输电线路保护性能。

Description

一种基于纯电流特征的直流输电线路保护方法和***
技术领域
本发明涉及电力***继电保护领域,并且更具体地，涉及一种基于纯电流特征的直流输电线路保护方法和***。
背景技术
在与传统交流输电***相比，高压直流输电***具有输送容量大、传输距离远、损耗低等优点，在远距离输电、大区域的电网互联、地下电缆输电等方面得到了广泛应用。高压直流输电线路肩负能源产地和负荷中心时间段电能传输的重任，其输电距离远，运行条件恶劣，故障率相对直流***其他部分更高，约占直流***故障的50％。因此，高性能的高压直流输电线路保护对提高整个电网的安全稳定性具有重要意义。
目前，直流输电线路保护主要包括行波保护、微分欠压保护、纵联差动保护等。行波保护、欠压微分保护作为直流输电线路主保护，可以快速响应直流线路故障(3～5ms)，但耐受过渡电阻能力差，动作性能易受到雷击、异常大数、噪声干扰的影响，受极间互感影响非故障极线路保护易误动。纵联差动保护作为直流输电线路后备保护，具备良好的反映高阻接地故障的能力，但在外部短路时，长距离输电线路的分布电容电流会产生较大差流，为防止保护误动，其动作延时设置较长(百毫秒甚至秒级)。
发明内容
为了解决现在技术中行波保护、欠压微分保护耐受过渡电阻能力差，动作性能易受到雷击、异常大数、噪声干扰的影响，而纵联差动保护动作延时较长的技术问题，本发明提供一种基于纯电流特征的直流输电线路保护方法，所述方法包括：
采集第一直流输电线路M侧在k0时刻的线路电流iM(k0)和k0-ts时刻的线路电流iM(k0-ts)，并根据所述线路电流iM(k0)和iM(k0-ts)计算k0时刻第一直流输电线路M侧电流的差分值diM(k0)，并确定所述电流的差分值diM(k0)是否满足保护启动判据，其中，所述第一直流输电线路为输电线路中的任意一条，所述第一直流输电线路M侧是所述第一直流输电线路整流侧和逆变侧中的一侧，ts为采样时间间隔；
当k0时刻第一直流输电线路M侧电流的差分值diM(k0)满足保护启动判据时，启动所述M侧的直流保护,并采集第一直流输电线路M侧在j时刻的线路电流iM(j)和j-ts时刻的线路电流iM(j-ts)，并根据所述线路电流iM(j)和iM(j-ts)计算j时刻第一直流输电线路M侧电流的差分值diM(j)，采集第一直流输电线路M侧在k时刻的线路电流iM(k)和k-ts时刻的线路电流iM(k-ts)，并根据所述线路电流iM(k)和iM(k-ts)计算k时刻第一直流输电线路M侧电流的差分值diM(k)，采集第一直流输电线路N侧在j-Ttran时刻的线路电流iN(j-Ttran)和j-ts-Ttran时刻的线路电流iN(j-ts-Ttran)，并根据所述线路电流iN(j-Ttran)和iN(j-ts-Ttran)计算j-Ttran时刻第一直流输电线路N侧电流的差分值diN(j-Ttran)，采集第二直流输电线路M侧在j时刻的线路电流iM'(j)和j-ts时刻的线路电流iM'(j-ts)，并根据所述线路电流iM'(j)和iM'(j-ts)计算j时刻第二直流输电线路M侧电流的差分值diM'(j)，以及根据所述j时刻第一直流输电线路M侧电流的差分值diM(j)计算k时刻第一直流输电线路M侧电流突变量△iM(k)，根据所述j-Ttran时刻第一直流输电线路N侧电流的差分值diN(j-Ttran)计算k-Ttran时刻第一直流输电线路N侧电流突变量△iN(k-Ttran)，其中，将所述启动保护的时刻记为t0，所述第二直流输电线路为除第一直流输电线路外的其他直流输电线路，第二直流输电线路M侧与第一直流输电线路M侧位于同一侧，所述第一直流输电线路N侧为所述直流输电线路中与M侧相对的另一侧，Ttran为直流线路传输通道延时，t0≤j≤k，t0≤k≤t；
根据第一直流输电线路M侧电流的差分值diM(k)和N侧电流突变量△iN(k-Ttran)计算第一直流输电线路差分积聚动作量iΣ(t)，根据第一直流输电线路M侧电流的差分值diM(k)计算第一直流输电线路差分积聚门槛值isetz，并确定所述第一直流输电线路差分积聚动作量iΣ(t)和差分积聚门槛值isetz是否满足预先设置的对端助增差分电流积聚量判据，其中，
t0≤k≤t；
根据第一直流输电线路M侧电流突变量△iM(k)和N侧电流突变量
△iN(k-Ttran)计算第一直流输电线路方向动作量iΣ△(t)，并确定所述第一直流输电线路方向动作量iΣ△(t)是否满足预先设置的对端助增突变量电流的积聚量高定值方向判据、对端助增突变量电流的积聚量低定值方向判据，其中，t0≤k≤t；
在保护启动后的tlimit时间内，确定所述j时刻第一直流输电线路M侧电流的差分值diM(j)和第二直流输电线路M侧电流的差分值diM'(j)是否满足预先设置的第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据，在保护启动t0～tf时间内，确定所述j时刻第一直流输电线路M侧电流的差分值diM(j)是否满足预先设置的第一极比幅式突变量电流选极判据，在保护启动tf时间后，确定第一直流输电线路方向动作量iΣ△(t)是否满足比率制动电流选极判据，其中，当在t0～tlimit时间段持续同时不满足第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据、在t0～tf中某一个时刻满足第一极比幅式突变量电流选极判据或者tf时间后任一时刻满足比率制动电流选极判据中的任意一个成立都表示阶段式电流选极判据成立，所述第一极是第一直流输电线路M侧，所述第二极是第二直流输电线路M侧，t0≤j≤k，t0≤k≤t，tf>tlimit；
根据k时刻第一直流输电线路M侧电流的差分值diM(k)和第一直流输电线路N侧电流突变量△iN(k-Ttran)计算k时刻第一直流输电线路M侧的差分量i△(k)，基于所述差分量i△(k)确定k时刻第一直流输电线路M侧的置反差分量i△'(k)，并确定所述k时刻第一直流输电线路M侧的置反差分量i△'(k)是否满足预先设置的利用差分积聚量的防大数保护判据；
根据第一直流输电线路M侧电流突变量△iM(k)和N侧电流突变量△iN(k-Ttran)计算k时刻第一直流输电线路防大数突变量iz(k)，并确定k时刻第一直流输电线路防大数突变量iz(k)是否满足利用电流突变量的防大数保护判据；
当所述t时刻第一直流输电线路的对端助增差分电流积聚量判据、对端助增突变量电流的积聚量高定值方向判据和利用差分积聚量的防大数保护判据同时成立，在t0～t时间段满足利用电流突变量的防大数保护判据的采样点数大于预先设置的点数阈值，且满足利用电流突变量的防大数保护判据的采样点数大于不满足利用电流突变量的防大数保护判据的采样点数时，或者所述t时刻第一直流输电线路对端助增差分电流积聚量判据、对端助增突变量电流的积聚量低定值方向判据、阶段式电流选极判据和利用差分积聚量的防大数保护判据同时成立，在t0～t时间段满足利用电流突变量的防大数保护判据的采样点数大于预先设置的点数阈值，且满足利用电流突变量的防大数保护判据的采样点数大于不满足利用电流突变量的防大数保护判据的采样点数时，第一直流输电线路M侧保护动作出口，其中，t0≤k≤t。
进一步地，所述根据k0时刻第一直流输电线路M侧的线路电流iM(k0)和iM(k0-ts)计算k0时刻第一直流输电线路M侧电流的差分值diM(k0)，当所述电流的差分值diM(k0)满足保护启动判据时，启动所述M侧的直流保护包括：
根据k0时刻第一直流输电线路M侧的线路电流iM(k0)和iM(k0-ts)计算k0时刻第一直流输电线路M侧电流的差分值diM(k0)，其计算公式为：
diM(k0)＝iM(k0)-iM(k0-ts)
式中，ts为采样时间间隔；
确定所述电流的差分值diM(k0)是否满足保护启动判据，其中，所述保护启动判据的计算公式为：
|diM(k0)|>iset0
式中，iset0为启动门槛，按照本线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，计算|diM(k0)|的值，利用所述|diM(k0)|的值再除以一个大于1的系数，将其结果作为启动门槛iset0。
进一步地，根据所述线路电流iM(j)和iM(j-ts)计算j时刻第一直流输电线路M侧电流的差分值diM(j)，其计算公式为：
diM(j)＝iM(j)-iM(j-ts)
根据所述线路电流iM(k)和iM(k-ts)计算k时刻第一直流输电线路M侧电流的差分值diM(k)，其计算公式为：
diM(k)＝iM(k)-iM(k-ts)
根据所述线路电流iN(j-Ttran)和iN(j-ts-Ttran)计算j-Ttran时刻第一直流输电线路N侧电流的差分值diN(j-Ttran)，其计算公式为：
diN(j-Ttran)＝iN(j-Ttran)-iN(j-ts-Ttran)
根据所述线路电流iM'(j)和iM'(j-ts)计算j时刻第二直流输电线路M侧电流的差分值diM'(j)，其计算公式为：
diM′(j)＝iM′(j)-iM′(j-ts)
根据所述j时刻第一直流输电线路M侧电流的差分值diM(j)计算k时刻第一直流输电线路M侧电流突变量△iM(k)，其计算公式为：
Figure GDA0002273247850000051
根据所述j-Ttran时刻第一直流输电线路N侧电流的差分值
diN(j-Ttran)计算k-Ttran时刻第一直流输电线路N侧电流突变量△iN(k-Ttran)，其计算公式为：
Figure GDA0002273247850000052
式中，t0≤j≤k，t0≤k≤t，ts为采样时间间隔，Ttran为直流线路传输通道延时。
进一步地，所述根据k时刻第一直流输电线路M侧电流的差分值diM(k)和N侧电流突变量△iN(k-Ttran)计算第一直流输电线路差分积聚动作量iΣ(t)，根据k时刻第一直流输电线路M侧电流的差分值diM(k)计算第一直流输电线路差分积聚门槛值isetz，并确定所述t时刻第一直流输电线路差分积聚动作量iΣ(t)和差分积聚门槛值isetz是否满足预先设置的对端助增差分电流积聚量判据包括：
根据k时刻第一直流输电线路M侧电流的差分值diM(k)和N侧电流突变量△iN(k-Ttran)计算t时刻第一直流输电线路差分积聚动作量iΣ(t)，其计算公式为：
Figure GDA0002273247850000061
式中，n为t0～t时间段的采样点数，|diM(k)|为k时刻第一直流输电线路M侧电流的差分值diM(k)的绝对值；
根据第一直流输电线路M侧电流的差分值diM(k)计算第一直流输电线路差分积聚门槛值isetz，其计算公式为：
isetz＝kk(iset1+γ)
Figure GDA0002273247850000062
式中，kk为可靠系数，iset1为固定门槛，确定为区外金属性故障时的差分积聚动作量在故障后T0时间内的最大值，γ为负斜率电流浮动门槛，
Figure GDA0002273247850000063
为负斜率电流，T为浮动门槛计算窗长，ρ为浮动门槛系数，ρmax为比例常数，ρmax>1，T'为定值计算窗长，T'>T，λ为裕度系数，λ>1；
确定所述第一直流输电线路差分积聚动作量iΣ(t)和差分积聚门槛值isetz是否满足预先设置的对端助增差分电流积聚量判据，其中，所述对端助增差分电流积聚量判据的计算公式为：
iΣ(t)>isetz。
进一步地，所述根据第一直流输电线路M侧电流突变量△iM(k)和N侧电流突变量△iN(k-Ttran)计算第一直流输电线路方向动作量iΣ△(t)，并确定所述第一直流输电线路方向动作量iΣ△(t)是否满足预先设置的对端助增突变量电流的积聚量高定值方向判据、对端助增突变量电流的积聚量低定值方向判据包括：
根据第一直流输电线路M侧电流突变量△iM(k)和N侧电流突变量△iN(k-Ttran)计算第一直流输电线路方向动作量iΣ△(t)，其计算公式为：
Figure GDA0002273247850000071
式中，n为t0～t时间段的采样点数；
确定所述第一直流输电线路方向动作量iΣ△(t)是否满足预先设置的对端助增突变量电流的积聚量高定值方向判据、对端助增突变量电流的积聚量低定值方向判据，其中，所述对端助增突变量电流的积聚量高定值方向判据的计算公式为：
iΣ△(t)>iset2H
所述对端助增突变量电流的积聚量低定值方向判据的计算公式为：
iΣ△(t)>iset2L
式中，iset2H为对端助增突变量电流的积聚量方向判别高定值，iset2L为方向判别低定值，iset2H按照第二直流输电线路末端金属性故障整定，即在第二直流输电线路末端生成一个金属性故障，并在所述故障时计算第一直流输电线路方向动作量iΣ△(t)的值，并将所述iΣ△(t)的值乘以一个大于1的常数的结果作为iset2H，iset2L按照第一直流输电线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，并在所述故障时计算第一直流输电线路方向动作量iΣ△(t)的值，并将所述iΣ△(t)的值除以一个大于1的常数的结果作为iset2L。
进一步地，所述在保护启动后的tlimit时间内，确定所述j时刻第一直流输电线路M侧电流的差分值diM(j)和第二直流输电线路M侧电流的差分值diM'(j)是否满足预先设置的第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据，在保护启动t0～tf时间内，确定所述j时刻第一直流输电线路M侧电流的差分值diM(j)是否满足预先设置的第一极比幅式突变量电流选极判据，在保护启动tf时间后，确定第一直流输电线路方向动作量iΣ△(t)是否满足比率制动电流选极判据包括：
在保护启动后的tlimit时间内，确定所述j时刻第一直流输电线路M侧电流的差分值diM(j)和第二直流输电线路M侧电流的差分值diM'(j)是否满足预先设置的第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据，其中，所述第一极比幅式突变量电流选极判据的计算公式为：
Figure GDA0002273247850000081
所述第二极比幅式突变量电流选极判据的计算公式为：
Figure GDA0002273247850000082
当在t0～tf时间段内任一时刻第一极比幅式突变量电流选极判据成立，则阶段式电流选极判据成立，当在t0～tlimit时间段内第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据持续不成立时，在tlimit时刻及tlimit～tf时间段内判定发生双极故障，所述阶段式电流选极判据成立；
在保护启动tf时间后，确定第一直流输电线路方向动作量iΣ△(t)是否满足比率制动电流选极判据，其中，所述比率制动电流选极判据的计算公式为：
iΣ△(t)>αiD△(t)
Figure GDA0002273247850000083
式中，σ为选极系数，σ>1，α为比率制动系数，α>1，αiD△(t)为方向制动量，保护启动tf时间后，以比率制动电流选极判据的结果作为阶段式电流选极判据的结果，保护启动tf时间后任意一个时刻比率制动电流选极判据成立则阶段式电流选极判据成立。
进一步地，所述根据k时刻第一直流输电线路M侧电流的差分值diM(k)和第一直流输电线路N侧电流突变量△iN(k-Ttran)计算k时刻第一直流输电线路M侧的差分量i△(k)，基于所述差分量i△(k)确定k时刻第一直流输电线路M侧的置反差分量i△'(k)，并确定所述k时刻第一直流输电线路M侧的置反差分量i△'(k)是否满足预先设置的利用差分积聚量的防大数保护判据包括：
根据k时刻第一直流输电线路M侧电流的差分值diM(k)和第一直流输电线路N侧电流突变量△iN(k-Ttran)计算k时刻第一直流输电线路M侧的差分量i△(k)，其计算公式为：
i△(k)＝|diM(k)|+△iN(k-Ttran)
基于所述差分量i△(k)确定k时刻第一直流输电线路M侧的置反差分量i△'(k)，其计算公式为：
i'△(k)＝f(i△(k))
式中，f为置反处理方法，是将t0～t时间段绝对值最大的i△(k)强制乘以-ν作为k时刻的置反差分量，而将其他时刻的i△(k)作为对应时刻的置反差分量，ν>1；
确定所述k时刻第一直流输电线路M侧的置反差分量i△'(k)是否满足预先设置的利用差分积聚量的防大数保护判据，其中，所述利用差分积聚量的防大数保护判据的计算公式为：
if(t)>iset3
Figure GDA0002273247850000091
式中，if(t)为防大数差分动作量，iset3为防大数差分定值，按照第一直流输电线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，并在所述故障时计算第一直流输电线路防大数差分动作量if(t)的值，并将所述if(t)的值除以一个大于1的常数的结果作为iset3。
进一步地，所述根据第一直流输电线路M侧电流突变量△iM(k)和N侧电流突变量△iN(k-Ttran)计算k时刻第一直流输电线路防大数突变量iz(k)，并确定k时刻第一直流输电线路防大数突变量iz(k)是否满足利用电流突变量的防大数保护判据包括：
根据第一直流输电线路M侧电流突变量△iM(k)和N侧电流突变量
△iN(k-Ttran)计算k时刻第一直流输电线路防大数突变量iz(k)，其计算公式为：
iZ(k)＝△iM(k)+△iN(k-Ttran)
确定k时刻第一直流输电线路防大数突变量iz(k)是否满足利用电流突变量的防大数保护判据，其中，所述利用电流突变量的防大数保护判据的计算公式为：
|iZ(k)|>iset4
式中，iset4为防大数突变量定值，按照第一直流输电线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，并在所述故障时计算第一直流输电线路防大数突变量iz(k)的值，并将所述iz(k)的值除以一个大于1的常数的结果作为iset4。
根据本发明的另一方面，本发明提供一种基于纯电流特征的直流输电线路保护***，所述***包括：
初始设置单元，其用于设置第一直流输电线路和第二直流输电线路，以及第一直流输电线路的M侧和N侧，以及和第二直流输电线路的M侧，设置启动保护判据、对端助增差分电流积聚量判据、对端助增突变量电流的积聚量高定值方向判据、对端助增突变量电流的积聚量低定值方向判据、第一极比幅式突变量电流选极判据、第二极比幅式突变量电流选极判据、比率制动电流选极判据、利用差分积聚量的防大数保护判据和利用电流突变量的防大数保护判据的计算公式，以及为上述判据中具备定值的参数和采样点数阈值赋值，其中，所述第一直流输电线路为输电线路中的任意一条，第二直流输电线路为除第一直流输电线路外的其他直流输电线路，第一直流输电线路M侧是所述第一直流输电线路整流侧和逆变侧中的一侧，第二直流输电线路M侧与第一直流输电线路M侧位于同一侧，所述第一直流输电线路N侧为所述直流输电线路中与M侧相对的另一侧；
数据采集单元，其用于实时采集第一直流输电线路M侧和N侧的线路电流，以及第二直流输电线路M侧的线路电流；
数据处理单元，其用于根据第一直流输电线路M侧的线路电流计算第一直流输电线路M侧电流的差分值，根据第一直流输电线路N侧的线路电流计算第一直流输电线路N侧电流的差分值，根据第二直流输电线路M侧的线路电流计算第二直流输电线路M侧电流的差分值，根据第一直流输电线路M侧电流的差分值计算第一直流输电线路M侧的电流突变量，根据第一直流输电线路N侧电流的差分值计算第一直流输电线路N侧的电流突变量，并根据所述计算确定的结果判定初始设置单元中的判据是否成立；
保护启动单元，其用于在保护启动判据成立时，启动第一直流输电线路M侧的直流保护,并将所述启动保护的时刻记为t0；
保护出口单元，其用于在启动第一直流输电线路M侧的直流保护后，当所述t时刻第一直流输电线路的对端助增差分电流积聚量判据、对端助增突变量电流的积聚量高定值方向判据和利用差分积聚量的防大数保护判据同时成立，在t0～t时间段满足利用电流突变量的防大数保护判据的采样点数大于预先设置的点数阈值，且满足利用电流突变量的防大数保护判据的采样点数大于不满足利用电流突变量的防大数保护判据的采样点数时，或者所述t时刻第一直流输电线路对端助增差分电流积聚量判据、对端助增突变量电流的积聚量低定值方向判据、阶段式电流选极判据和利用差分积聚量的防大数保护判据同时成立，在t0～t时间段满足利用电流突变量的防大数保护判据的采样点数大于预先设置的点数阈值，且满足利用电流突变量的防大数保护判据的采样点数大于不满足利用电流突变量的防大数保护判据的采样点数时，第一直流输电线路M侧保护动作出口，其中，当在t0～tlimit时间段内持续同时不满足第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据、在t0～tf时间段中任一个时刻满足第一极比幅式突变量电流选极判据或者在tf时间以后任一时刻满足比率制动电流选极判据中的任意一个成立都表示阶段式电流选极判据成立，所述第一极是第一直流输电线路M侧，所述第二极是第二直流输电线路M侧，所述t为第一直流输电线路M侧的直流保护启动后，数据采集单元采样结束的时刻，tlimit和tf为保护启动后的两个时刻，t0≤j≤k，t0≤k≤t，tf>tlimit。
进一步地，所述数据处理单元根据第一直流输电线路M侧的线路电流计算第一直流输电线路M侧电流的差分值，根据第一直流输电线路N侧的线路电流计算第一直流输电线路N侧电流的差分值，根据第二直流输电线路M侧的线路电流计算第二直流输电线路M侧电流的差分值，根据第一直流输电线路M侧电流的差分值计算第一直流输电线路M侧的电流突变量，根据第一直流输电线路N侧电流的差分值计算第一直流输电线路N侧的电流突变量，其计算公式分别为：
diM(j)＝iM(j)-iM(j-ts)
式中，diM(j)为j时刻第一直流输电线路M侧电流的差分值，iM(j)为第一直流输电线路M侧在j时刻的线路电流，iM(j-ts)为j-ts时刻的线路电流，ts为采样时间间隔；
diN(j-Ttran)＝iN(j-Ttran)-iN(j-ts-Ttran)
式中，diN(j-Ttran)为j-Ttran时刻第一直流输电线路N侧电流的差分值，
iN(j-Ttran)为第一直流输电线路N侧在j-Ttran时刻的线路电流，iN(j-ts-Ttran)为j-ts-Ttran时刻的线路电流，ts为采样时间间隔，Ttran为直流线路传输通道延时，t0≤j≤t；
diM′(j)＝iM′(j)-iM′(j-ts)
式中，diM'(j)为j时刻第二直流输电线路M侧电流的差分值，iM'(j)为第二直流输电线路M侧在j时刻的线路电流，iM'(j-ts)为j-ts时刻的线路电流，ts为采样时间间隔；
Figure GDA0002273247850000131
式中，△iM(k)为k时刻第一直流输电线路M侧电流突变量，diM(j)为j时刻第一直流输电线路M侧电流的差分值，△iN(k-Ttran)为k-Ttran时刻第一直流输电线路N侧电流突变量△iN(k-Ttran)，diN(j-Ttran)为j-Ttran时刻第一直流输电线路N侧电流的差分值，t0≤j≤k，t0≤k≤t，t0为保护启动时刻，ts为采样时间间隔，Ttran为直流线路传输通道延时。
进一步地，所述初始设置单元设置的启动保护判据、对端助增差分电流积聚量判据、对端助增突变量电流的积聚量高定值方向判据、对端助增突变量电流的积聚量低定值方向判据、第一极比幅式突变量电流选极判据、第二极比幅式突变量电流选极判据、比率制动电流选极判据、利用差分积聚量的防大数保护判据和利用电流突变量的防大数保护判据的计算公式包括：
第一直流输电线路保护启动判据的计算公式为：
|diM(k0)|>iset0
式中，|diM(k0)|为k0j时刻第一直流输电线路M侧电流的差分值diM(k0)的绝对值，iset0为启动门槛，按照本线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，计算|diM(k0)|的值，利用所述|diM(k0)|的值再除以一个大于1的系数，将其结果作为启动门槛iset0；
第一直流输电线路对端助增差分电流积聚量判据的计算公式为：
iΣ(t)>isetz
Figure GDA0002273247850000132
isetz＝kk(iset1+γ)
Figure GDA0002273247850000141
式中，iΣ(t)为t时刻第一直流输电线路差分积聚动作量，isetz为差分积聚门槛值，n为t0～t时间段的采样点数，|diM(k)|为k时刻第一直流输电线路M侧电流的差分值diM(k)的绝对值，kk为可靠系数，iset1为固定门槛，确定为区外金属性故障时的差分积聚动作量在故障后T0时间内的最大值，γ为负斜率电流浮动门槛，
Figure GDA0002273247850000143
为负斜率电流，T为浮动门槛计算窗长，ρ为浮动门槛系数，ρmax为比例常数，ρmax>1，T'为定值计算窗长，T'>T，λ为裕度系数，λ>1；
第一直流输电线路对端助增突变量电流的积聚量高定值方向判据的计算公式为：
iΣ△(t)>iset2H
Figure GDA0002273247850000142
式中，iΣ△(t)为t时刻第一直流输电线路方向动作量，iset2H为对端助增突变量电流的积聚量方向判别高定值，iset2H按照第二直流输电线路末端金属性故障整定，即在第二直流输电线路末端生成一个金属性故障，并在所述故障时计算第一直流输电线路方向动作量iΣ△(t)的值，并将所述iΣ△(t)的值乘以一个大于1的常数的结果作为iset2H，△iM(k)为k时刻第一直流输电线路M侧电流突变量，△iN(k-Ttran)为k时刻第一直流输电线路N侧电流突变量，n为t0～t时间段的采样点数；
第一直流输电线路对端助增突变量电流的积聚量低定值方向判据的计算公式为：
iΣ△(t)>iset2L
Figure GDA0002273247850000151
式中，iΣ△(t)为t时刻第一直流输电线路方向动作量，iset2L为方向判别低定值，iset2L按照第一直流输电线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，并在所述故障时计算第一直流输电线路方向动作量iΣ△(t)的值，并将所述iΣ△(t)的值除以一个大于1的常数的结果作为iset2L，△iM(k)为k时刻第一直流输电线路M侧电流突变量，△iN(k-Ttran)为k时刻第一直流输电线路N侧电流突变量，n为t0～t时间段的采样点数，iset2H>iset2L；
第一极比幅式突变量电流选极判据的计算公式为：
Figure GDA0002273247850000152
第二极比幅式突变量电流选极判据的计算公式为：
Figure GDA0002273247850000153
第一直流输电线路比率制动电流选极判据的计算公式为：
iΣ△(t)>αiD△(t)
Figure GDA0002273247850000154
式中，diM(j)为j时刻第一直流输电线路M侧电流的差分值，diM'(j)为j时刻第二直流输电线路M侧电流的差分值，σ为选极系数，σ>1，iΣ(t)为t时刻第一直流输电线路差分积聚动作量，α为比率制动系数，α>1，
αiD△(t)为方向制动量，t0≤j≤k，t0≤k≤t，ts为采样时间间隔，Ttran为直流线路传输通道延时，当在t0～tf时间段内任一时刻第一极比幅式突变量电流选极判据成立，则阶段式电流选极判据成立，当在t0～tlimit时间段内第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据持续不成立时，在tlimit时刻及tlimit～tf时间段内判定发生双极故障，所述阶段式电流选极判据成立，保护启动tf时间后任意一个时刻，以比率制动电流选极判据的结果作为阶段式电流选极判据的结果，比率制动电流选极判据成立则阶段式电流选极判据成立；
利用差分积聚量的防大数保护判据的计算公式为：
if(t)>iset3
Figure GDA0002273247850000161
i'△(k)＝f(i△(k))
i△(k)＝|diM(k)|+△iN(k-Ttran)
式中，if(t)为防大数差分动作量，iset3为防大数差分定值，按照第一直流输电线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，并在所述故障时计算第一直流输电线路防大数差分动作量if(t)的值，并将所述if(t)的值除以一个大于1的常数的结果作为iset3，i△'(k)为k时刻第一直流输电线路M侧的置反差分量，f为置反处理方法，是将t0～t时间段绝对值最大的i△(k)强制乘以-ν作为k时刻的置反差分量，而将其他时刻的i△(k)作为对应时刻的置反差分量，ν>1，i△(k)为k时刻第一直流输电线路M侧的差分量，diM(k)为k时刻第一直流输电线路M侧电流的差分值，△iN(k-Ttran)为k-Ttran时刻第一直流输电线路N侧电流突变量，t0≤k≤t，ts为采样时间间隔，Ttran为直流线路传输通道延时；
第一直流输电线路利用电流突变量的防大数保护判据的计算公式为：
|iZ(k)|>iset4
iZ(k)＝△iM(k)+△iN(k-Ttran)
式中，iz(k)为k时刻第一直流输电线路防大数突变量，iset4为防大数突变量定值，按照第一直流输电线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，并在所述故障时计算第一直流输电线路防大数突变量|iz(k)|的值，并将所述|iz(k)|的值除以一个大于1的常数的结果作为iset4，△iM(k)为k-Ttran时刻第一直流输电线路M侧电流突变量，△iN(k-Ttran)为第一直流输电线路时刻N侧电流突变量。
本发明技术方案提供的基于纯电流特征的直流输电线路保护方法和***通过采集直流输电线路电流确定线路的电流差分值和电流突变量，并根据所述线路电流差分值和电流突变量确定预先设置的保护判据是否成立，来确定直流保护动作是否出口。本发明所述的基于纯电流特征的直流输电线路保护方法和***能够进一步加快金属性故障保护动作速度，大幅提升高阻故障保护动作速度，雷击或异常大数时保护可靠不动作，单极接地故障工况下非故障极线路保护可靠不误动，全面提升直流输电线路保护性能。
附图说明
通过参考下面的附图，可以更为完整地理解本发明的示例性实施方式：
图1为根据本发明优选实施方式的基于纯电流特征的直流输电线路保护的方法的流程图；
图2为根据本发明优选实施方式的直流输电线路的示意图；
图3为根据本发明优选实施方式的第一直流输电线路的对端助增时序示意图；
图4为根据本发明优选实施方式的第一直流输电线路的阶段式电流选极判据时序示意图；
图5为根据本发明优选实施方式的基于纯电流特征的直流输电线路保护逻辑示意图；
图6(a)为根据本发明优选实施方式的第一直流输电线路末端发生高阻故障时差分积聚动作量和差分积聚门槛值随时间变化的波形示意图；
图6(b)为根据本发明优选实施方式的第一直流输电线路末端发生高阻故障时方向动作量和方向判别低定值随时间变化的波形示意图；
图6(c)为根据本发明优选实施方式的第一直流输电线路末端发生高阻故障时第一极比幅式突变量电流幅值和第二极比幅式突变量电流幅值随时间变化的波形示意图；
图6(d)为根据本发明优选实施方式的第一直流输电线路末端发生高阻故障时防大数差分动作量和防大数差分定值随时间变化的波形示意图；
图6(e)为根据本发明优选实施方式的第一直流输电线路末端发生高阻故障时防大数突变量绝对值和防大数突变量定值随时间变化的波形示意图；
图6(f)为根据本发明优选实施方式的第一直流输电线路末端发生高阻故障时满足利用电流突变量的防大数保护判据的点数和不满足利用电流突变量的防大数保护判据的点数随时间变化的波形示意图；
图7(a)为根据本发明优选实施方式的第一直流输电线路末端发生双极金属性故障时差分积聚动作量和差分积聚门槛值随时间变化的波形示意图；
图7b)为根据本发明优选实施方式的第一直流输电线路末端发生双极金属性故障时方向动作量和方向判别高定值随时间变化的波形示意图；
图8(a)为根据本发明优选实施方式的第二直流输电线路末端发生金属性故障时方向动作量和方向判别高定值随时间变化的波形示意图；
图8(b)为根据本发明优选实施方式的第二直流输电线路末端发生金属性故障时方向动作量和方向判别低定值随时间变化的波形示意图；
图8(c)为根据本发明优选实施方式的第二直流输电线路末端发生金属性故障时第一极比幅式突变量电流幅值和第二极比幅式突变量电流幅值随时间变化的波形示意图；
图8(d)为根据本发明优选实施方式的第二直流输电线路末端发生故障时方向动作量和方向制动量随时间变化的波形示意图；
图9为根据本发明优选实施方式的基于纯电流特征的直流输电线路保护***的结构示意图。
具体实施方式
现在参考附图介绍本发明的示例性实施方式，然而，本发明可以用许多不同的形式来实施，并且不局限于此处描述的实施例，提供这些实施例是为了详尽地且完全地公开本发明，并且向所属技术领域的技术人员充分传达本发明的范围。对于表示在附图中的示例性实施方式中的术语并不是对本发明的限定。在附图中，相同的单元/元件使用相同的附图标记。
除非另有说明，此处使用的术语(包括科技术语)对所属技术领域的技术人员具有通常的理解含义。另外，可以理解的是，以通常使用的词典限定的术语，应当被理解为与其相关领域的语境具有一致的含义，而不应该被理解为理想化的或过于正式的意义。
图1为根据本发明优选实施方式的基于纯电流特征的直流输电线路保护的方法的流程图。如图1所示，本优选实施方式所述的基于纯电流特征的直流输电线路保护的方法100从步骤101开始。
在步骤101，采集第一直流输电线路M侧在k0时刻的线路电流iM(k0)和k0-ts时刻的线路电流iM(k0-ts)，并根据所述线路电流iM(k0)和iM(k0-ts)计算k0时刻第一直流输电线路M侧电流的差分值diM(k0)，并确定所述电流的差分值diM(k0)是否满足保护启动判据，其中，所述第一直流输电线路为输电线路中的任意一条，所述第一直流输电线路M侧是所述第一直流输电线路整流侧和逆变侧中的一侧，ts为采样时间间隔。
图2为根据本发明优选实施方式的直流输电线路的示意图。如图2所示，在本优选实施方式中，所述第一直流输电线路为极Ⅰ线路，所述第一直流输电线路M侧是极Ⅰ线路的整流侧，所述第一直流输电线路N侧是极Ⅰ线路的逆变侧，iM为直流线路M侧电流值，iN为直流线路N侧电流值，电流正方向为极母线指向线路。
优选地，所述根据k0时刻第一直流输电线路M侧的线路电流iM(k0)和iM(k0-ts)计算k0时刻第一直流输电线路M侧电流的差分值diM(k0)，当所述电流的差分值diM(k0)满足保护启动判据时，启动所述M侧的直流保护包括：
根据k0时刻第一直流输电线路M侧的线路电流iM(k0)和iM(k0-ts)计算k0时刻第一直流输电线路M侧电流的差分值diM(k0)，其计算公式为：
diM(k0)＝iM(k0)-iM(k0-ts)
式中，ts为采样时间间隔；
确定所述电流的差分值diM(k0)是否满足保护启动判据，其中，所述保护启动判据的计算公式为：
|diM(k0)|>iset0
式中，iset0为启动门槛，按照本线路末端高阻故障有灵敏度整定，例如±500kV***按照线路末端800欧姆过渡电阻有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，计算|diM(k0)|的值，利用所述|diM(k0)|的值再除以一个大于1的系数，将其结果作为启动门槛iset0。
在步骤102，当k0时刻第一直流输电线路M侧电流的差分值diM(k0)满足保护启动判据时，启动所述M侧的直流保护,并采集第一直流输电线路M侧在j时刻的线路电流iM(j)和j-ts时刻的线路电流iM(j-ts)，并根据所述线路电流iM(j)和iM(j-ts)计算j时刻第一直流输电线路M侧电流的差分值diM(j)，采集第一直流输电线路M侧在k时刻的线路电流iM(k)和k-ts时刻的线路电流iM(k-ts)，并根据所述线路电流iM(k)和iM(k-ts)计算k时刻第一直流输电线路M侧电流的差分值diM(k)，采集第一直流输电线路N侧在j-Ttran时刻的线路电流iN(j-Ttran)和j-ts-Ttran时刻的线路电流iN(j-ts-Ttran)，并根据所述线路电流iN(j-Ttran)和iN(j-ts-Ttran)计算j-Ttran时刻第一直流输电线路N侧电流的差分值diN(j-Ttran)，采集第二直流输电线路M侧在j时刻的线路电流iM'(j)和j-ts时刻的线路电流iM'(j-ts)，并根据所述线路电流iM'(j)和iM'(j-ts)计算j时刻第二直流输电线路M侧电流的差分值diM'(j)，以及根据所述j时刻第一直流输电线路M侧电流的差分值diM(j)计算k时刻第一直流输电线路M侧电流突变量△iM(k)，根据所述j-Ttran时刻第一直流输电线路N侧电流的差分值diN(j-Ttran)计算k-Ttran时刻第一直流输电线路N侧电流突变量△iN(k-Ttran)，其中，将所述启动保护的时刻记为t0，所述第二直流输电线路为除第一直流输电线路外的其他直流输电线路，第二直流输电线路M侧与第一直流输电线路M侧位于同一侧，所述第一直流输电线路N侧为所述直流输电线路中与M侧相对的另一侧，Ttran为直流线路传输通道延时，t0≤j≤k，t0≤k≤t。
优选地，根据所述线路电流iM(j)和iM(j-ts)计算j时刻第一直流输电线路M侧电流的差分值diM(j)，其计算公式为：
diM(j)＝iM(j)-iM(j-ts)
根据所述线路电流iM(k)和iM(k-ts)计算k时刻第一直流输电线路M侧电流的差分值diM(k)，其计算公式为：
diM(k)＝iM(k)-iM(k-ts)
根据所述线路电流iN(j-Ttran)和iN(j-ts-Ttran)计算j-Ttran时刻第一直流输电线路N侧电流的差分值diN(j-Ttran)，其计算公式为：
diN(j-Ttran)＝iN(j-Ttran)-iN(j-ts-Ttran)
根据所述线路电流iM'(j)和iM'(j-ts)计算j时刻第二直流输电线路M侧电流的差分值diM'(j)，其计算公式为：
diM′(j)＝iM′(j)-iM′(j-ts)
根据所述j时刻第一直流输电线路M侧电流的差分值diM(j)计算k时刻第一直流输电线路M侧电流突变量△iM(k)，其计算公式为：
Figure GDA0002273247850000211
根据所述j-Ttran时刻第一直流输电线路N侧电流的差分值diN(j-Ttran)计算k-Ttran时刻第一直流输电线路N侧电流突变量△iN(k-Ttran)，其计算公式为：
Figure GDA0002273247850000212
式中，t0≤j≤k，t0≤k≤t，ts为采样时间间隔，Ttran为直流线路传输通道延时。
在步骤103，根据第一直流输电线路M侧电流的差分值diM(k)和N侧电流突变量△iN(k-Ttran)计算第一直流输电线路差分积聚动作量iΣ(t)，根据第一直流输电线路M侧电流的差分值diM(k)计算第一直流输电线路差分积聚门槛值isetz，并确定所述第一直流输电线路差分积聚动作量iΣ(t)和差分积聚门槛值isetz是否满足预先设置的对端助增差分电流积聚量判据，其中，t0≤k≤t。
优选地，所述根据k时刻第一直流输电线路M侧电流的差分值diM(k)和N侧电流突变量△iN(k-Ttran)计算第一直流输电线路差分积聚动作量iΣ(t)，根据k时刻第一直流输电线路M侧电流的差分值diM(k)计算第一直流输电线路差分积聚门槛值isetz，并确定所述t时刻第一直流输电线路差分积聚动作量iΣ(t)和差分积聚门槛值isetz是否满足预先设置的对端助增差分电流积聚量判据包括：
根据k时刻第一直流输电线路M侧电流的差分值diM(k)和N侧电流突变量△iN(k-Ttran)计算t时刻第一直流输电线路差分积聚动作量iΣ(t)，其计算公式为：
Figure GDA0002273247850000221
式中，n为t0～t时间段的采样点数，|diM(k)|为k时刻第一直流输电线路M侧电流的差分值diM(k)的绝对值；
根据第一直流输电线路M侧电流的差分值diM(k)计算第一直流输电线路差分积聚门槛值isetz，其计算公式为：
isetz＝kk(iset1+γ)
Figure GDA0002273247850000222
Figure GDA0002273247850000231
式中，kk为可靠系数，iset1为固定门槛，确定为区外金属性故障时的差分积聚动作量在故障后T0时间内的最大值，建议T0<50ms，γ为负斜率电流浮动门槛，
Figure GDA0002273247850000232
为负斜率电流，T为浮动门槛计算窗长，建议T<5ms，ρ为浮动门槛系数，ρmax为比例常数，ρmax>1，T'为定值计算窗长，T'>T，λ为裕度系数，λ>1；
确定所述第一直流输电线路差分积聚动作量iΣ(t)和差分积聚门槛值isetz是否满足预先设置的对端助增差分电流积聚量判据，其中，所述对端助增差分电流积聚量判据的计算公式为：
iΣ(t)>isetz
图3为根据本发明优选实施方式的第一直流输电线路的对端助增时序示意图。如图3所示，第一直流输电线路M侧直流保护启动时刻为t0，由于第一直流输电线路N侧电流突变量是采用的k-Ttran时刻的量，对于M侧判据，保护启动至Ttran时间段内，N侧电流突变量△iN(k-Ttran)为0，判据仅采用了M侧单端的电流量，称为单端量过程；在保护启动后Ttran时刻之后，N侧突变量传递过来，当区内故障时，N侧电流突变量△iN(k-Ttran)大于0，叠加了电流突变量△iN(k-Ttran)之后，动作量更大，判据更容易满足。当末端区外发生故障时，N侧电流突变量△iN(k-Ttran)小于0，叠加了△iN(k-Ttran)之后，动作量更小，判据更不容易满足，由此可见，对于第一直流输电线路M侧直流保护而言，N侧突变量传递过来后对判据的整体性能起到提升的作用，因此保护启动Ttran时刻之后称为对端助增过程。
在步骤104，根据第一直流输电线路M侧电流突变量△iM(k)和N侧电流突变量△iN(k-Ttran)计算第一直流输电线路方向动作量iΣ△(t)，并确定所述第一直流输电线路方向动作量iΣ△(t)是否满足预先设置的对端助增突变量电流的积聚量高定值方向判据、对端助增突变量电流的积聚量低定值方向判据。
优选地，所述根据第一直流输电线路M侧电流突变量△iM(k)和N侧电流突变量△iN(k-Ttran)计算第一直流输电线路方向动作量iΣ△(t)，并确定所述第一直流输电线路方向动作量iΣ△(t)是否满足预先设置的对端助增突变量电流的积聚量高定值方向判据、对端助增突变量电流的积聚量低定值方向判据包括：
根据第一直流输电线路M侧电流突变量△iM(k)和N侧电流突变量△iN(k-Ttran)计算第一直流输电线路方向动作量iΣ△(t)，其计算公式为：
Figure GDA0002273247850000241
式中，n为t0～t时间段的采样点数；
确定所述第一直流输电线路方向动作量iΣ△(t)是否满足预先设置的对端助增突变量电流的积聚量高定值方向判据、对端助增突变量电流的积聚量低定值方向判据，其中，所述对端助增突变量电流的积聚量高定值方向判据的计算公式为：
iΣ△(t)>iset2H
所述对端助增突变量电流的积聚量低定值方向判据的计算公式为：
iΣ△(t)>iset2L
式中，iset2H为方向判别高定值，iset2L为方向判别低定值，iset2H按照第二直流输电线路末端金属性故障整定，即在第二直流输电线路末端生成一个金属性故障，并在所述故障时计算第一直流输电线路方向动作量iΣ△(t)的值，并将所述iΣ△(t)的值乘以一个大于1的常数的结果作为iset2H，iset2L按照第一直流输电线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，并在所述故障时计算第一直流输电线路方向动作量iΣ△(t)的值，并将所述iΣ△(t)的值除以一个大于1的常数的结果作为iset2L。
在步骤105，在保护启动后的tlimit时间内，确定所述j时刻第一直流输电线路M侧电流的差分值diM(j)和第二直流输电线路M侧电流的差分值diM'(j)是否满足预先设置的第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据，在保护启动t0～tf时间内，确定所述j时刻第一直流输电线路M侧电流的差分值diM(j)是否满足预先设置的第一极比幅式突变量电流选极判据，在保护启动tf时间后，确定第一直流输电线路方向动作量iΣ△(t)是否满足比率制动电流选极判据，其中，当在t0～tlimit时间段持续同时不满足第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据、在t0～tf中某一个时刻满足第一极比幅式突变量电流选极判据或者tf时间后某一时刻满足比率制动电流选极判据中的任意一个成立都表示阶段式电流选极判据成立，所述第一极是第一直流输电线路M侧，所述第二极是第二直流输电线路M侧，t0≤j≤k，t0≤k≤t，tf>tlimit。
优选地，所述在保护启动后的tlimit时间内，确定所述j时刻第一直流输电线路M侧电流的差分值diM(j)和第二直流输电线路M侧电流的差分值diM'(j)是否满足预先设置的第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据，在保护启动t0～tf时间内，确定所述j时刻第一直流输电线路M侧电流的差分值diM(j)是否满足预先设置的第一极比幅式突变量电流选极判据，在保护启动tf时间后，确定第一直流输电线路方向动作量iΣ△(t)是否满足比率制动电流选极判据包括：
在保护启动后的tlimit时间内，确定所述j时刻第一直流输电线路M侧电流的差分值diM(j)和第二直流输电线路M侧电流的差分值diM'(j)是否满足预先设置的第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据，其中，所述第一极比幅式突变量电流选极判据的计算公式为：
Figure GDA0002273247850000251
所述第二极比幅式突变量电流选极判据的计算公式为：
Figure GDA0002273247850000252
当在t0～tf时间段内任一时刻第一极比幅式突变量电流选极判据成立，则阶段式电流选极判据成立，当在t0～tlimit时间段内第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据持续不成立时，在tlimit时刻及tlimit～tf时间段内判定发生双极故障，所述阶段式电流选极判据成立；
在保护启动tf时间后，确定第一直流输电线路方向动作量iΣ△(t)是否满足比率制动电流选极判据，其中，所述比率制动电流选极判据的计算公式为：
iΣ△(t)>αiD△(t)
Figure GDA0002273247850000261
式中，σ为选极系数，σ>1，α为比率制动系数，α>1，αiD△(t)为方向制动量，保护启动tf时间后，以比率制动电流选极判据的结果作为阶段式电流选极判据的结果，比率制动电流选极判据成立则阶段式电流选极判据成立。图4为根据本发明优选实施方式的第一直流输电线路的阶段式电流选极判据时序示意图。如图4所示，第一直流输电线路M侧保护启动时刻为t0，当在t0～tf时间段内任一时刻第一极比幅式突变量电流选极判据成立，则阶段式电流选极判据成立，当在t0～tlimit时间段内第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据持续不成立时，在tlimit时刻及tlimit～tf时间段内判定发生双极故障，所述阶段式电流选极判据成立，保护启动tf时间后，以比率制动电流选极判据的结果作为阶段式电流选极选极判据的结果，比率制动电流选极判据成立则阶段式电流选极选极判据成立。
在步骤106，根据k时刻第一直流输电线路M侧电流的差分值diM(k)和第一直流输电线路N侧电流突变量△iN(k-Ttran)计算k时刻第一直流输电线路M侧的差分量i△(k)，基于所述差分量i△(k)确定k时刻第一直流输电线路M侧的置反差分量i△'(k)，并确定所述k时刻第一直流输电线路M侧的置反差分量i△'(k)是否满足预先设置的利用差分积聚量的防大数保护判据。
优选地，所述根据k时刻第一直流输电线路M侧电流的差分值diM(k)和第一直流输电线路N侧电流突变量△iN(k-Ttran)计算k时刻第一直流输电线路M侧的差分量i△(k)，基于所述差分量i△(k)确定k时刻第一直流输电线路M侧的置反差分量i△'(k)，并确定所述k时刻第一直流输电线路M侧的置反差分量i△'(k)是否满足预先设置的利用差分积聚量的防大数保护判据包括：
根据k时刻第一直流输电线路M侧电流的差分值diM(k)和第一直流输电线路N侧电流突变量△iN(k-Ttran)计算k时刻第一直流输电线路M侧的差分量i△(k)，其计算公式为：
i△(k)＝|diM(k)|+△iN(k-Ttran)
基于所述差分量i△(k)确定k时刻第一直流输电线路M侧的置反差分量i△'(k)，其计算公式为：
i'△(k)＝f(i△(k))
式中，f为置反处理方法，是将t0～t时间段绝对值最大的i△(k)强制乘以-ν作为k时刻的置反差分量，而将其他时刻的i△(k)作为对应时刻的置反差分量，ν>1；
确定所述k时刻第一直流输电线路M侧的置反差分量i△'(k)是否满足预先设置的利用差分积聚量的防大数保护判据，其中，所述利用差分积聚量的防大数保护判据的计算公式为：
if(t)>iset3
Figure GDA0002273247850000271
式中，if(t)为防大数差分动作量，iset3为防大数差分定值，按照第一直流输电线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，并在所述故障时计算第一直流输电线路防大数差分动作量if(t)的值，并将所述if(t)的值除以一个大于1的常数的结果作为iset3。
在步骤107，根据第一直流输电线路M侧电流突变量△iM(k)和N侧电流突变量△iN(k-Ttran)计算k时刻第一直流输电线路防大数突变量iz(k)，并确定k时刻第一直流输电线路防大数突变量iz(k)是否满足利用电流突变量的防大数保护判据。
优选地，所述根据第一直流输电线路M侧电流突变量△iM(k)和N侧电流突变量△iN(k-Ttran)计算k时刻第一直流输电线路防大数突变量iz(k)，并确定k时刻第一直流输电线路防大数突变量iz(k)是否满足利用电流突变量的防大数保护判据包括：
根据第一直流输电线路M侧电流突变量△iM(k)和N侧电流突变量
△iN(k-Ttran)计算k时刻第一直流输电线路防大数突变量iz(k)，其计算公式为：
iZ(k)＝△iM(k)+△iN(k-Ttran)
确定k时刻第一直流输电线路防大数突变量iz(k)是否满足利用电流突变量的防大数保护判据，其中，所述利用电流突变量的防大数保护判据的计算公式为：
|iZ(k)|>iset4
式中，iset4为防大数突变量定值，按照第一直流输电线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，并在所述故障时计算第一直流输电线路防大数突变量|iz(k)|的值，并将所述|iz(k)|的值除以一个大于1的常数的结果作为iset4。
在步骤108，当所述t时刻第一直流输电线路的对端助增差分电流积聚量判据、对端助增突变量电流的积聚量高定值方向判据和利用差分积聚量的防大数保护判据同时成立，在t0～t时间段满足利用电流突变量的防大数保护判据的采样点数大于预先设置的点数阈值，且满足利用电流突变量的防大数保护判据的采样点数大于不满足利用电流突变量的防大数保护判据的采样点数时，或者所述t时刻第一直流输电线路对端助增差分电流积聚量判据、对端助增突变量电流的积聚量低定值方向判据、阶段式电流选极判据和利用差分积聚量的防大数保护判据同时成立，在t0～t时间段满足利用电流突变量的防大数保护判据的采样点数大于预先设置的点数阈值，且满足利用电流突变量的防大数保护判据的采样点数大于不满足利用电流突变量的防大数保护判据的采样点数时，第一直流输电线路M侧保护动作出口。
图5为根据本发明优选实施方式的基于纯电流特征的直流输电线路保护逻辑示意图。如图5所示，本优选实施方式所述的基于纯电流特征的直流输电线路保护逻辑为对端助增突变量电流的积聚量高定值方向判据与对端助增差分电流积聚量判据经“与”门逻辑，对端助增突变量电流的积聚量低定值方向判据、对端助增差分电流积聚量保护判据与阶段式电流选极方法经“与”门逻辑，两个“与”门逻辑结果再经“或”门逻辑，利用差分积聚量的防大数保护判据与利用电流突变量的防大数保护判据经“与”门逻辑，“或”门逻辑与“与”门逻辑再经“与”门逻辑出口。
在本优选实施方式中，当如图2所示的第一直流输电线路末端F3发生800Ω高阻故障时：
图6(a)为根据本发明优选实施方式的第一直流输电线路末端发生高阻故障时差分积聚动作量和差分积聚门槛值随时间变化的波形示意图。如图6(a)所示，虚线为差分积聚动作量随时间变化波形，实线为差分积聚门槛值随时间变化波形。由波形可知，在第一直流输电线路M侧保护启动后0ms，差分积聚动作量大于差分积聚门槛值，对端助增差分电流积聚量判据满足。
图6(b)为根据本发明优选实施方式的第一直流输电线路末端发生高阻故障时方向动作量和方向判别低定值随时间变化的波形示意图。如图6(b)所示，虚线为方向动作量随时间变化波形，实线为方向判别低定值随时间变化波形，由于其是根据直流输电线路末端发生高阻故障时进行灵敏度整定得到的值，故为常数，因此随时间变化波形是一条横线。由波形可知，在第一直流输电线路M侧保护启动后3ms，方向动作量大于方向判别低定值，对端助增突变量电流的积聚量低定值方向判据满足。
图6(c)为根据本发明优选实施方式的第一直流输电线路末端发生高阻故障时第一极比幅式突变量电流幅值和第二极比幅式突变量电流幅值随时间变化的波形示意图。如图6(c)所示，虚线为第一极比幅式突变量电流幅值随时间变化波形，实线为第二极比幅式突变量电流幅值随时间变化波形。由波形可知，在第一直流输电线路M侧保护启动后1.6ms，第一极比幅式突变量电流幅值大于第二极比幅式突变量电流幅值，第一极比幅式突变量电流选极判据满足，故阶段式电流选极判据满足。
图6(d)为根据本发明优选实施方式的第一直流输电线路末端发生高阻故障时防大数差分动作量和防大数差分定值随时间变化的波形示意图。如图6(d)所示，虚线为防大数差分动作量随时间变化波形，实线为防大数差分定值随时间变化波形，由于其是根据直流输电线路末端发生高阻故障时进行灵敏度整定得到的值，故为常数，因此随时间变化波形是一条横线。由波形可知，在第一直流输电线路M侧保护启动后0.2ms，防大数差分动作量大于防大数差分定值，利用差分积聚量的防大数保护判据满足。
图6(e)为根据本发明优选实施方式的第一直流输电线路末端发生高阻故障时防大数突变量绝对值和防大数突变量定值随时间变化的波形示意图。如图6(e)所示，虚线为防大数突变量绝对值随时间变化波形，实线为防大数突变量定值随时间变化波形，由于其是根据直流输电线路末端发生高阻故障时进行灵敏度整定得到的值，故为常数，因此随时间变化波形是一条横线。由波形可知，在第一直流输电线路M侧保护启动后0ms，防大数突变量绝对值大于防大数突变量定值，利用电流突变量的防大数保护判据满足。
图6(f)为根据本发明优选实施方式的第一直流输电线路末端发生高阻故障时满足利用电流突变量的防大数保护判据的点数和不满足利用电流突变量的防大数保护判据的点数随时间变化的波形示意图。如图6(f)所示，实线为满足利用电流突变量的防大数保护判据的点数随时间变化波形，虚线为不满足利用电流突变量的防大数保护判据的点数随时间变化波形，点划线为常数m，m为预先设置的利用电流突变量的防大数保护判据的点数阈值，其随时间变化的波形是一条横线。由波形可知，在第一直流输电线路M侧保护启动后0.7ms，满足利用电流突变量的防大数保护判据的点数大于常数m且大于不满足利用电流突变量的防大数保护判据的点数，利用电流突变量的防大数保护判据满足。
根据图5所示的基于纯电流特征的直流输电线路保护逻辑可知，对端助增突变量电流的积聚量低定值方向判据、对端助增差分电流积聚量保护判据与阶段式电流选极方法经“与”门逻辑，利用差分积聚量的防大数保护判据与利用电流突变量的防大数保护判据经“与”门逻辑后，两组判据再经与门后，第一直流输电线路的保护出口，因此，在本实施例中，在第一直流输电线路M侧保护启动3ms后，上述所有判据均满足，第一直流输电线路M例的直流保护出口。
在本优选实施方式中，当如图2所示的第一直流输电线路末端F3发生双极金属性故障时：
图7(a)为根据本发明优选实施方式的第一直流输电线路末端发生双极故障时差分积聚动作量和差分积聚门槛值随时间变化的波形示意图。如图7(a)所示，虚线为差分积聚动作量随时间变化波形，实线为差分积聚门槛值随时间变化波形。由波形可知，在第一直流输电线路M侧保护启动后0ms，差分积聚动作量大于差分积聚门槛值，对端助增差分电流积聚量判据满足。
图7(b)为根据本发明优选实施方式的第一直流输电线路末端发生双极故障时方向动作量和方向判别高定值随时间变化的波形示意图。如图7(b)所示，虚线为方向动作量随时间变化波形，实线为方向判别高定值随时间变化波形，由于其是根据第二直流输电线路末端发生金属性故障时进行整定得到的值，故为常数，因此随时间变化波形是一条横线。由波形可知，在第一直流输电线路M侧保护启动后0ms，方向动作量大于方向判别高定值，方向判别对端助增突变量电流的积聚量高定值判据满足。
由于在本优选实施方式中，在启动后0.2ms，防大数差分动作量大于防大数差分定值，利用差分积聚量的防大数保护判据满足，在启动后0.7ms，当满足利用电流突变量的防大数保护判据的点数大于常数m且大于不满足利用电流突变量的防大数保护判据的点数时，，根据图5所示的基于纯电流特征的直流输电线路保护逻辑可知，对端助增突变量电流的积聚量高定值方向判据与对端助增差分电流积聚量判据经“与”门逻辑，利用差分积聚量的防大数保护判据与利用电流突变量的防大数保护判据经“与”门逻辑，两组判据再经“与”门后，第一直流输电线路的保护出口，因此，在本实施例中，在第一直流输电线路M侧保护启动0.7ms后，上述所有判据均满足，第一直流输电线路M例的直流保护出口。
在本优选实施方式中，当如图2所示的第二直流输电线路末端发生金属性故障时：
第一直流输电线路M侧保护启动后0ms，差分积聚动作量大于差分积聚门槛值，对端助增差分电流积聚量判据满足。
图8(a)为根据本发明优选实施方式的第二直流输电线路末端发生金属性故障时方向动作量和方向判别高定值随时间变化的波形示意图。如图8(a)所示，虚线为方向动作量随时间变化波形，实线为方向判别高定值随时间变化波形，其为常数，因此随时间变化波形为一条横线。由波形可知，在第一直流输电线路M侧保护启动后500ms，方向动作量始终小于方向判别高定值，方向判别对端助增突变量电流的积聚量高定值方向判据不满足。
图8(b)为根据本发明优选实施方式的第二直流输电线路末端发生金属性故障时方向动作量和方向判别低定值随时间变化的波形示意图。如图8(b)所示，虚线为方向动作量随时间变化波形，实线为方向判别低定值随时间变化波形，其为常数，因此随时间变化波形为一条横线。由波形可知，在第一直流输电线路M侧保护启动后0.1ms，方向动作量大于方向判别低定值，对端助增突变量电流的积聚量低定值方向判据满足。
图8(c)为根据本发明优选实施方式的第二直流输电线路末端发生金属性故障时第一极比幅式突变量电流幅值和第二极比幅式突变量电流幅值随时间变化的波形示意图。如图8(c)所示，虚线为第一极比幅式突变量电流幅值随时间变化波形，实线为第二极比幅式突变量电流幅值随时间变化波形。由波形可知，在第一直流输电线路M侧保护启动后tf时间内，第一极比幅式突变量电流幅值小于第二极比幅式突变量电流幅值，第一极比幅式突变量电流选极判据不满足，第二极比幅式突变量电流选极判据满足，故阶段式电流选极判据不满足。
图8(d)为根据本发明优选实施方式的第二直流输电线路末端发生故障时方向动作量和方向制动量随时间变化的波形示意图。如图8(d)所示，虚线为方向动作量随时间变化波形，实线为方向制动量随时间变化波形。由波形可知，在启动tf时间后，方向动作量小于方向制动量，比率制动电流选极判据不满足，故阶段式电流选极判据不满足。
根据图5所示的基于纯电流特征的直流输电线路保护逻辑可知，在对端助增差分电流积聚量判据成立，而对端助增突变量电流的积聚量高定值方向判据不成立时，两者的“与”门逻辑不成立，或者在对端助增差分电流积聚量判据和对端助增突变量电流的积聚量低定值方向判据的“与”门逻辑成立，但阶段式电流选极判据不成立时，即使利用差分积累量的防大数保护判据和利用电流突变量的防大数保护判据成立，所述逻辑示意图中的最后一个“与”门逻辑也不成立，故第一直流输电线路M侧的保护启动后，其直流保护也不出口。
除本优选实施方式中的第一直流线路发生高阻故障、双极故障以及第二直流线路末端金属性故障，所述基于纯电流特征的直流输电线路保护方法还可适用于直流输电线路单极金属性故障、直流线路末端区外故障以及直流线路发生雷击或采样出现异常大数时，直流输电线路保护可靠出口或者不出口。
对于直流输电***单极运行时，将对端电流的电流值设为0，仍可按照上述基于纯电流特征的直流保护判据进行直流保护动作出口。
图9为根据本发明优选实施方式的基于纯电流特征的直流输电线路保护***的结构示意图。如图9所示，本优选实施方式所述的基于纯电流特征的直流输电线路保护***900包括：
初始设置单元901，其用于设置第一直流输电线路和第二直流输电线路，以及第一直流输电线路的M侧和N侧，以及和第二直流输电线路的M侧，设置启动保护判据、对端助增差分电流积聚量判据、对端助增突变量电流的积聚量高定值方向判据、对端助增突变量电流的积聚量低定值方向判据、第一极比幅式突变量电流选极判据、第二极比幅式突变量电流选极判据、比率制动电流选极判据、利用差分积聚量的防大数保护判据和利用电流突变量的防大数保护判据的计算公式，以及为上述判据中具备定值的参数和采样点数阈值赋值，其中，所述第一直流输电线路为输电线路中的任意一条，第二直流输电线路为除第一直流输电线路外的其他直流输电线路，第一直流输电线路M侧是所述第一直流输电线路整流侧和逆变侧中的一侧，第二直流输电线路M侧与第一直流输电线路M侧位于同一侧，所述第一直流输电线路N侧为所述直流输电线路中与M侧相对的另一侧；
数据采集单元902，其用于实时采集第一直流输电线路M侧和N侧的线路电流，以及第二直流输电线路M侧的线路电流；
数据处理单元903，其用于根据第一直流输电线路M侧的线路电流计算第一直流输电线路M侧电流的差分值，根据第一直流输电线路N侧的线路电流计算第一直流输电线路N侧电流的差分值，根据第二直流输电线路M侧的线路电流计算第二直流输电线路M侧电流的差分值，根据第一直流输电线路M侧电流的差分值计算第一直流输电线路M侧的电流突变量，根据第一直流输电线路N侧电流的差分值计算第一直流输电线路N侧的电流突变量，并根据所述计算确定的结果判定初始设置单元中的判据是否成立；
保护启动单元904，其用于在保护启动判据成立时，启动第一直流输电线路M侧的直流保护,并将所述启动保护的时刻记为t0；
保护出口单元905，其用于在启动第一直流输电线路M侧的直流保护后，当所述t时刻第一直流输电线路的对端助增差分电流积聚量判据、对端助增突变量电流的积聚量高定值方向判据和利用差分积聚量的防大数保护判据同时成立，在t0～t时间段满足利用电流突变量的防大数保护判据的采样点数大于预先设置的点数阈值，且满足利用电流突变量的防大数保护判据的采样点数大于不满足利用电流突变量的防大数保护判据的采样点数时，或者所述t时刻第一直流输电线路对端助增差分电流积聚量判据、对端助增突变量电流的积聚量低定值方向判据、阶段式电流选极判据和利用差分积聚量的防大数保护判据同时成立，在t0～t时间段满足利用电流突变量的防大数保护判据的采样点数大于预先设置的点数阈值，且满足利用电流突变量的防大数保护判据的采样点数大于不满足利用电流突变量的防大数保护判据的采样点数时，第一直流输电线路M侧保护动作出口，其中，当在t0～tlimit时间段内持续同时不满足第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据、在t0～tf时间段中任一个时刻满足第一极比幅式突变量电流选极判据或者在tf时间后满足比率制动电流选极判据中的任意一个成立都表示阶段式电流选极判据成立，所述第一极是第一直流输电线路M侧，所述第二极是第二直流输电线路M侧，所述t为第一直流输电线路M侧的直流保护启动后，数据采集单元采样结束的时刻，tlimit和tf为保护启动后的两个时刻，tf>tlimit。
优选地，所述数据处理单元903根据第一直流输电线路M侧的线路电流计算第一直流输电线路M侧电流的差分值，根据第一直流输电线路N侧的线路电流计算第一直流输电线路N侧电流的差分值，根据第二直流输电线路M侧的线路电流计算第二直流输电线路M侧电流的差分值，根据第一直流输电线路M侧电流的差分值计算第一直流输电线路M侧的电流突变量，根据第一直流输电线路N侧电流的差分值计算第一直流输电线路N侧的电流突变量，其计算公式分别为：
diM(j)＝iM(j)-iM(j-ts)
式中，diM(j)为j时刻第一直流输电线路M侧电流的差分值，iM(j)为第一直流输电线路M侧在j时刻的线路电流，iM(j-ts)为j-ts时刻的线路电流，ts为采样时间间隔；
diN(j-Ttran)＝iN(j-Ttran)-iN(j-ts-Ttran)
式中，diN(j-Ttran)为j-Ttran时刻第一直流输电线路N侧电流的差分值，
iN(j-Ttran)为第一直流输电线路N侧在j-Ttran时刻的线路电流，iN(j-ts-Ttran)为j-ts-Ttran时刻的线路电流，ts为采样时间间隔，Ttran为直流线路传输通道延时，t0≤j≤t；
diM′(j)＝iM′(j)-iM′(j-ts)
式中，diM'(j)为j时刻第二直流输电线路M侧电流的差分值，iM'(j)为第二直流输电线路M侧在j时刻的线路电流，iM'(j-ts)为j-ts时刻的线路电流，ts为采样时间间隔；
Figure GDA0002273247850000361
Figure GDA0002273247850000362
式中，△iM(k)为k时刻第一直流输电线路M侧电流突变量，diM(j)为j时刻第一直流输电线路M侧电流的差分值，△iN(k-Ttran)为k-Ttran时刻第一直流输电线路N侧电流突变量△iN(k-Ttran)，diN(j-Ttran)为j-Ttran时刻第一直流输电线路N侧电流的差分值，t0≤j≤k，t0≤k≤t，t0为保护启动时刻，ts为采样时间间隔，Ttran为直流线路传输通道延时。
优选地，所述初始设置单元901设置的启动保护判据、对端助增差分电流积聚量判据、对端助增突变量电流的积聚量高定值方向判据、对端助增突变量电流的积聚量低定值方向判据、第一极比幅式突变量电流选极判据、第二极比幅式突变量电流选极判据、比率制动电流选极判据、利用差分积聚量的防大数保护判据和利用电流突变量的防大数保护判据的计算公式包括：
第一直流输电线路保护启动判据的计算公式为：
|diM(k0)|>iset0
式中，|diM(k0)|为k0时刻第一直流输电线路M侧电流的差分值diM(k0)的绝对值，iset0为启动门槛，按照本线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，并在所述故障时计算第一直流输电线路M侧电流的差分值diM(k0)的绝对值|diM(k0)|，并将所述|diM(k0)|的值除以一个大于1的常数的结果作为iset4；
第一直流输电线路对端助增差分电流积聚量判据的计算公式为：
iΣ(t)>isetz
Figure GDA0002273247850000363
isetz＝kk(iset1+γ)
Figure GDA0002273247850000371
式中，iΣ(t)为t时刻第一直流输电线路差分积聚动作量，isetz为差分积聚门槛值，n为t0～t时间段的采样点数，|diM(k)|为k时刻第一直流输电线路M侧电流的差分值diM(k)的绝对值，kk为可靠系数，iset1为固定门槛，确定为区外金属性故障时的差分积聚动作量在故障后T0时间内的最大值，γ为负斜率电流浮动门槛，
Figure GDA0002273247850000373
为负斜率电流，T为浮动门槛计算窗长，ρ为浮动门槛系数，ρmax为比例常数，ρmax>1，T'为定值计算窗长，T'>T，λ为裕度系数，λ>1；
第一直流输电线路对端助增突变量电流的积聚量高定值方向判据的计算公式为：
iΣ△(t)>iset2H
Figure GDA0002273247850000372
式中，iΣ△(t)为t时刻第一直流输电线路方向动作量，iset2H为方向判别对端助增突变量电流的积聚量高定值，iset2H按照第二直流输电线路末端金属性故障整定，即在第二直流输电线路末端生成一个金属性故障，并在所述故障时计算第一直流输电线路方向动作量iΣ△(t)的值，并将所述iΣ△(t)的值乘以一个大于1的常数的结果作为iset2H，△iM(k)为k时刻第一直流输电线路M侧电流突变量，△iN(k-Ttran)为k时刻第一直流输电线路N侧电流突变量，n为t0～t时间段的采样点数；
第一直流输电线路对端助增突变量电流的积聚量低定值方向判据的计算公式为：
iΣ△(t)>iset2L
Figure GDA0002273247850000381
式中，iΣ△(t)为t时刻第一直流输电线路方向动作量，iset2L为方向判别对端助增突变量电流的积聚量低定值，iset2L按照第一直流输电线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，并在所述故障时计算第一直流输电线路方向动作量iΣ△(t)，并将所述iΣ△(t)的值除以一个大于1的常数的结果作为iset2L，△iM(k)为k时刻第一直流输电线路M侧电流突变量，△iN(k-Ttran)为k时刻第一直流输电线路N侧电流突变量，n为t0～t时间段的采样点数，iset2H>iset2L；
第一极比幅式突变量电流选极判据的计算公式为：
Figure GDA0002273247850000382
第二极比幅式突变量电流选极判据的计算公式为：
Figure GDA0002273247850000383
第一直流输电线路比率制动电流选极判据的计算公式为：
iΣ△(t)>αiD△(t)
Figure GDA0002273247850000384
式中，diM(j)为j时刻第一直流输电线路M侧电流的差分值，diM'(j)为j时刻第二直流输电线路M侧电流的差分值，σ为选极系数，σ>1，iΣ(t)为t时刻第一直流输电线路差分积聚动作量，α为比率制动系数，α>1，αiD△(t)为方向制动量，t0≤j≤k，t0≤k≤t，ts为采样时间间隔，Ttran为直流线路传输通道延时，当在t0～tf时间段内任一时刻第一极比幅式突变量电流选极判据成立，则阶段式电流选极判据成立，当在t0～tlimit时间段内第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据持续不成立时，在tlimit时刻及tlimit～tf时间段内判定发生双极故障，所述阶段式电流选极判据成立，保护启动tf时间后，以比率制动电流选极判据的结果作为阶段式电流选极判据的结果，比率制动电流选极判据成立则阶段式电流选极判据成立；
利用差分积聚量的防大数保护判据的计算公式为：
if(t)>iset3
Figure GDA0002273247850000391
i'△(k)＝f(i△(k))
i△(k)＝|diM(k)|+△iN(k-Ttran)
式中，if(t)为防大数差分动作量，iset3为防大数差分定值，按照第一直流输电线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，并在所述故障时计算第一直流输电线路防大数差分动作量if(t)，并将所述if(t)的值除以一个大于1的常数的结果作为iset3，i△'(k)为k时刻第一直流输电线路M侧的置反差分积聚量，f为置反处理方法，是将t0～t时间段绝对值最大的i△(k)强制乘以-ν作为k时刻的置反差分量，而将其他时刻的i△(k)作为对应时刻的置反差分量，ν>1，i△(k)为k时刻第一直流输电线路M侧的差分量，diM(k)为k时刻第一直流输电线路M侧电流的差分值，△iN(k-Ttran)为k-Ttran时刻第一直流输电线路N侧电流突变量，t0≤k≤t，ts为采样时间间隔，Ttran为直流线路传输通道延时；
第一直流输电线路利用电流突变量的防大数保护判据的计算公式为：
|iZ(k)|>iset4
iZ(k)＝△iM(k)+△iN(k-Ttran)
式中，iz(k)为k时刻第一直流输电线路防大数突变量，iset4为防大数突变量定值，按照第一直流输电线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，并在所述故障时计算第一直流输电线路防大数突变量|iz(k)|，并将所述|iz(k)|的值除以一个大于1的常数的结果作为iset4，△iM(k)为k-Ttran时刻第一直流输电线路M侧电流突变量，△iN(k-Ttran)为第一直流输电线路时刻N侧电流突变量。
本发明所述基于纯电流特征的直流输电线路保护***对直流输电线路进行保护的方法与本发明所述直流输电线路保护判据的步骤相同，并且达到的技术效果也相同，此处不再赘述。
已经通过参考少量实施方式描述了本发明。然而，本领域技术人员所公知的，正如附带的专利权利要求所限定的，除了本发明以上公开的其他的实施例等同地落在本发明的范围内。
通常地，在权利要求中使用的所有术语都根据他们在技术领域的通常含义被解释，除非在其中被另外明确地定义。所有的参考“一个/所述/该[装置、组件等]”都被开放地解释为所述装置、组件等中的至少一个实例，除非另外明确地说明。这里公开的任何方法的步骤都没必要以公开的准确的顺序运行，除非明确地说明。

Claims (11)

1.一种基于纯电流特征的直流输电线路保护方法，其特征在于，所述方法包括：
采集第一直流输电线路M侧在k0时刻的线路电流iM(k0)和k0-ts时刻的线路电流iM(k0-ts)，并根据所述线路电流iM(k0)和iM(k0-ts)计算k0时刻第一直流输电线路M侧电流的差分值diM(k0)，并确定所述电流的差分值diM(k0)是否满足保护启动判据，其中，所述第一直流输电线路为输电线路中的任意一条，所述第一直流输电线路M侧是所述第一直流输电线路整流侧和逆变侧中的一侧，ts为采样时间间隔；
当k0时刻第一直流输电线路M侧电流的差分值diM(k0)满足保护启动判据时，启动所述M侧的直流保护,并采集第一直流输电线路M侧在j时刻的线路电流iM(j)和j-ts时刻的线路电流iM(j-ts)，并根据所述线路电流iM(j)和iM(j-ts)计算j时刻第一直流输电线路M侧电流的差分值diM(j)，采集第一直流输电线路M侧在k时刻的线路电流iM(k)和k-ts时刻的线路电流iM(k-ts)，并根据所述线路电流iM(k)和iM(k-ts)计算k时刻第一直流输电线路M侧电流的差分值diM(k)，采集第一直流输电线路N侧在j-Ttran时刻的线路电流iN(j-Ttran)和j-ts-Ttran时刻的线路电流iN(j-ts-Ttran)，并根据所述线路电流iN(j-Ttran)和iN(j-ts-Ttran)计算j-Ttran时刻第一直流输电线路N侧电流的差分值diN(j-Ttran)，采集第二直流输电线路M侧在j时刻的线路电流iM'(j)和j-ts时刻的线路电流iM'(j-ts)，并根据所述线路电流iM'(j)和iM'(j-ts)计算j时刻第二直流输电线路M侧电流的差分值diM'(j)，以及根据所述j时刻第一直流输电线路M侧电流的差分值diM(j)计算k时刻第一直流输电线路M侧电流突变量△iM(k)，根据所述j-Ttran时刻第一直流输电线路N侧电流的差分值diN(j-Ttran)计算k-Ttran时刻第一直流输电线路N侧电流突变量△iN(k-Ttran)，其中，将所述启动保护的时刻记为t0，所述第二直流输电线路为除第一直流输电线路外的其他直流输电线路，第二直流输电线路M侧与第一直流输电线路M侧位于同一侧，所述第一直流输电线路N侧为所述直流输电线路中与M侧相对的另一侧，Ttran为直流线路传输通道延时，t0≤j≤k，t0≤k≤t，t为启动保护后的当前采样时刻；
根据第一直流输电线路M侧电流的差分值diM(k)和N侧电流突变量△iN(k-Ttran)计算第一直流输电线路差分积聚动作量iΣ(t)，根据第一直流输电线路M侧电流的差分值diM(k)计算第一直流输电线路差分积聚门槛值isetz，并确定所述第一直流输电线路差分积聚动作量iΣ(t)和差分积聚门槛值isetz是否满足预先设置的对端助增差分电流积聚量判据；
根据第一直流输电线路M侧电流突变量△iM(k)和N侧电流突变量△iN(k-Ttran)计算第一直流输电线路方向动作量iΣ△(t)，并确定所述第一直流输电线路方向动作量iΣ△(t)是否满足预先设置的对端助增突变量电流的积聚量高定值方向判据、对端助增突变量电流的积聚量低定值方向判据；
在保护启动后的tlimit时间内，确定所述j时刻第一直流输电线路M侧电流的差分值diM(j)和第二直流输电线路M侧电流的差分值diM'(j)是否满足预先设置的第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据，在保护启动t0～tf时间内，确定所述j时刻第一直流输电线路M侧电流的差分值diM(j)是否满足预先设置的第一极比幅式突变量电流选极判据，在保护启动tf时间后，确定第一直流输电线路方向动作量iΣ△(t)是否满足比率制动电流选极判据，其中，当在t0～tlimit时间段持续同时不满足第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据、在t0～tf中某一个时刻满足第一极比幅式突变量电流选极判据或者tf时间后某一时刻满足比率制动电流选极判据中的任意一个成立都表示阶段式电流选极判据成立，所述第一极是第一直流输电线路M侧，所述第二极是第二直流输电线路M侧，tf＞tlimit；
根据k时刻第一直流输电线路M侧电流的差分值diM(k)和第一直流输电线路N侧电流突变量△iN(k-Ttran)计算k时刻第一直流输电线路M侧的差分量i△(k)，基于所述差分量i△(k)确定k时刻第一直流输电线路M侧的置反差分量i△'(k)，并确定所述k时刻第一直流输电线路M侧的置反差分量i△'(k)是否满足预先设置的利用差分积聚量的防大数保护判据；
根据第一直流输电线路M侧电流突变量△iM(k)和N侧电流突变量△iN(k-Ttran)计算k时刻第一直流输电线路防大数突变量iz(k)，并确定k时刻第一直流输电线路防大数突变量iz(k)是否满足利用电流突变量的防大数保护判据；
当所述t时刻第一直流输电线路的对端助增差分电流积聚量判据、对端助增突变量电流的积聚量高定值方向判据和利用差分积聚量的防大数保护判据同时成立，在t0～t时间段满足利用电流突变量的防大数保护判据的采样点数大于预先设置的点数阈值，且满足利用电流突变量的防大数保护判据的采样点数大于不满足利用电流突变量的防大数保护判据的采样点数时，或者所述t时刻第一直流输电线路对端助增差分电流积聚量判据、对端助增突变量电流的积聚量低定值方向判据、阶段式电流选极判据和利用差分积聚量的防大数保护判据同时成立，在t0～t时间段满足利用电流突变量的防大数保护判据的采样点数大于预先设置的点数阈值，且满足利用电流突变量的防大数保护判据的采样点数大于不满足利用电流突变量的防大数保护判据的采样点数时，第一直流输电线路M侧保护动作出口。
2.根据权利要求1所述的方法，其特征在于，所述根据k0时刻第一直流输电线路M侧的线路电流iM(k0)和iM(k0-ts)计算k0时刻第一直流输电线路M侧电流的差分值diM(k0)，当所述电流的差分值diM(k0)满足保护启动判据时，启动所述M侧的直流保护包括：
根据k0时刻第一直流输电线路M侧的线路电流iM(k0)和iM(k0-ts)计算k0时刻第一直流输电线路M侧电流的差分值diM(k0)，其计算公式为：
diM(k0)＝iM(k0)-iM(k0-ts)
式中，ts为采样时间间隔；
确定所述电流的差分值diM(k0)是否满足保护启动判据，其中，所述保护启动判据的计算公式为：
|diM(k0)|＞iset0
式中，iset0为启动门槛，按照本线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，计算|diM(k0)|的值，利用所述|diM(k0)|的值再除以一个大于1的系数，将其结果作为启动门槛iset0。
3.根据权利要求1所述的方法，其特征在于，根据所述线路电流iM(j)和iM(j-ts)计算j时刻第一直流输电线路M侧电流的差分值diM(j)，其计算公式为：
diM(j)＝iM(j)-iM(j-ts)
根据所述线路电流iM(k)和iM(k-ts)计算k时刻第一直流输电线路M侧电流的差分值diM(k)，其计算公式为：
diM(k)＝iM(k)-iM(k-ts)
根据所述线路电流iN(j-Ttran)和iN(j-ts-Ttran)计算j-Ttran时刻第一直流输电线路N侧电流的差分值diN(j-Ttran)，其计算公式为：
diN(j-Ttran)＝iN(j-Ttran)-iN(j-ts-Ttran)
根据所述线路电流iM'(j)和iM'(j-ts)计算j时刻第二直流输电线路M侧电流的差分值diM'(j)，其计算公式为：
diM′(j)＝iM′(j)-iM′(j-ts)
根据所述j时刻第一直流输电线路M侧电流的差分值diM(j)计算k时刻第一直流输电线路M侧电流突变量△iM(k)，其计算公式为：
Figure FDA0003209259660000041
根据所述j-Ttran时刻第一直流输电线路N侧电流的差分值diN(j-Ttran)计算k-Ttran时刻第一直流输电线路N侧电流突变量△iN(k-Ttran)，其计算公式为：
Figure FDA0003209259660000051
式中，t0≤j≤k，t0≤k≤t，ts为采样时间间隔，Ttran为直流线路传输通道延时。
4.根据权利要求3所述的方法，其特征在于，所述根据k时刻第一直流输电线路M侧电流的差分值diM(k)和N侧电流突变量△iN(k-Ttran)计算第一直流输电线路差分积聚动作量iΣ(t)，根据k时刻第一直流输电线路M侧电流的差分值diM(k)计算第一直流输电线路差分积聚门槛值isetz，并确定所述t时刻第一直流输电线路差分积聚动作量iΣ(t)和差分积聚门槛值isetz是否满足预先设置的对端助增差分电流积聚量判据包括：
根据k时刻第一直流输电线路M侧电流的差分值diM(k)和N侧电流突变量△iN(k-Ttran)计算t时刻第一直流输电线路差分积聚动作量iΣ(t)，其计算公式为：
Figure FDA0003209259660000052
式中，n为t0～t时间段的采样点数，t0≤k≤t，|diM(k)|为k时刻第一直流输电线路M侧电流的差分值diM(k)的绝对值；
根据第一直流输电线路M侧电流的差分值diM(k)计算第一直流输电线路差分积聚门槛值isetz，其计算公式为：
isetz＝kk(iset1+γ)
Figure FDA0003209259660000053
Figure FDA0003209259660000054
Figure FDA0003209259660000061
Figure FDA0003209259660000062
式中，kk为可靠系数，iset1为固定门槛，确定为区外金属性故障时的差分积聚动作量在故障后T0时间内的最大值，γ为负斜率电流浮动门槛，
Figure FDA0003209259660000063
为负斜率电流，T为浮动门槛计算窗长，ρ为浮动门槛系数，ρmax为比例常数，ρmax＞1，T'为定值计算窗长，T'＞T，λ为裕度系数，λ＞1；
确定所述第一直流输电线路差分积聚动作量iΣ(t)和差分积聚门槛值isetz是否满足预先设置的对端助增差分电流积聚量判据，其中，所述对端助增差分电流积聚量判据的计算公式为：
iΣ(t)＞isetz。
5.根据权利要求3所述的方法，其特征在于，所述根据第一直流输电线路M侧电流突变量△iM(k)和N侧电流突变量△iN(k-Ttran)计算第一直流输电线路方向动作量iΣ△(t)，并确定所述第一直流输电线路方向动作量iΣ△(t)是否满足预先设置的对端助增突变量电流的积聚量高定值方向判据、对端助增突变量电流的积聚量低定值方向判据包括：
根据第一直流输电线路M侧电流突变量△iM(k)和N侧电流突变量△iN(k-Ttran)计算第一直流输电线路方向动作量iΣ△(t)，其计算公式为：
Figure FDA0003209259660000064
式中，n为t0～t时间段的采样点数，t0≤k≤t；
确定所述第一直流输电线路方向动作量iΣ△(t)是否满足预先设置的对端助增突变量电流的积聚量高定值方向判据、对端助增突变量电流的积聚量低定值方向判据，其中，所述对端助增突变量电流的积聚量高定值方向判据的计算公式为：
iΣ△(t)＞iset2H
所述对端助增突变量电流的积聚量低定值方向判据的计算公式为：
iΣ△(t)＞iset2L
式中，iset2H为方向判别高定值，iset2L为方向判别低定值，iset2H按照第二直流输电线路末端金属性故障整定，即在第二直流输电线路末端生成一个金属性故障，并在所述故障时计算第一直流输电线路方向动作量iΣ△(t)的值，并将所述iΣ△(t)的值乘以一个大于1的常数的结果作为iset2H，iset2L按照第一直流输电线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，并在所述故障时计算第一直流输电线路方向动作量iΣ△(t)的值，并将所述iΣ△(t)的值除以一个大于1的常数的结果作为iset2L。
6.根据权利要求5所述的方法，其特征在于，所述在保护启动后的tlimit时间内，确定所述j时刻第一直流输电线路M侧电流的差分值diM(j)和第二直流输电线路M侧电流的差分值diM'(j)是否满足预先设置的第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据，在保护启动t0～tf时间内，确定所述j时刻第一直流输电线路M侧电流的差分值diM(j)是否满足预先设置的第一极比幅式突变量电流选极判据；在保护启动tf时间后，确定第一直流输电线路方向动作量iΣ△(t)是否满足比率制动电流选极判据包括：
在保护启动后的tlimit时间内，确定所述j时刻第一直流输电线路M侧电流的差分值diM(j)和第二直流输电线路M侧电流的差分值diM'(j)是否满足预先设置的第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据，其中，所述第一极比幅式突变量电流选极判据的计算公式为：
Figure FDA0003209259660000071
所述第二极比幅式突变量电流选极判据的计算公式为：
Figure FDA0003209259660000081
当在t0～tf时间段内任一时刻第一极比幅式突变量电流选极判据成立，则阶段式电流选极判据成立，当在t0～tlimit时间段内第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据持续不成立时，在tlimit时刻及tlimit～tf时间段内判定发生双极故障，所述阶段式电流选极判据成立；
在保护启动tf时间后，确定第一直流输电线路方向动作量iΣ△(t)是否满足比率制动电流选极判据，其中，所述比率制动电流选极判据的计算公式为：
iΣ△(t)＞αiD△(t)
Figure FDA0003209259660000082
式中，σ为选极系数，σ＞1，α为比率制动系数，α＞1，αiD△(t)为方向制动量，保护启动tf时间后，以比率制动电流选极判据的结果作为阶段式电流选极判据的结果，保护启动tf时间后任意一个时刻比率制动电流选极判据成立则阶段式电流选极判据成立。
7.根据权利要求3所述的方法，其特征在于，所述根据k时刻第一直流输电线路M侧电流的差分值diM(k)和第一直流输电线路N侧电流突变量△iN(k-Ttran)计算k时刻第一直流输电线路M侧的差分量i△(k)，基于所述差分量i△(k)确定k时刻第一直流输电线路M侧的置反差分量i△'(k)，并确定所述k时刻第一直流输电线路M侧的置反差分量i△'(k)是否满足预先设置的利用差分积聚量的防大数保护判据包括：
根据k时刻第一直流输电线路M侧电流的差分值diM(k)和第一直流输电线路N侧电流突变量△iN(k-Ttran)计算k时刻第一直流输电线路M侧的差分量i△(k)，其计算公式为：
i△(k)＝|diM(k)|+△iN(k-Ttran)
基于所述差分量i△(k)确定k时刻第一直流输电线路M侧的置反差分量i△'(k)，其计算公式为：
i'△(k)＝f(i△(k))
式中，f为置反处理方法，是将t0～t时间段绝对值最大的i△(k)强制乘以-ν作为k时刻的置反差分量，而将其他时刻的i△(k)作为对应时刻的置反差分量，ν＞1；
确定所述k时刻第一直流输电线路M侧的置反差分量i△'(k)是否满足预先设置的利用差分积聚量的防大数保护判据，其中，所述利用差分积聚量的防大数保护判据的计算公式为：
if(t)＞iset3
Figure FDA0003209259660000091
式中，if(t)为防大数差分动作量，iset3为防大数差分定值，按照第一直流输电线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，并在所述故障时计算第一直流输电线路防大数差分动作量if(t)的值，并将所述if(t)的值除以一个大于1的常数的结果作为iset3，t0≤k≤t。
8.根据权利要求3所述的方法，其特征在于，所述根据第一直流输电线路M侧电流突变量△iM(k)和N侧电流突变量△iN(k-Ttran)计算k时刻第一直流输电线路防大数突变量iz(k)，并确定k时刻第一直流输电线路防大数突变量iz(k)是否满足利用电流突变量的防大数保护判据包括：
根据第一直流输电线路M侧电流突变量△iM(k)和N侧电流突变量△iN(k-Ttran)计算k时刻第一直流输电线路防大数突变量iz(k)，其计算公式为：
iZ(k)＝△iM(k)+△iN(k-Ttran)
确定k时刻第一直流输电线路防大数突变量iz(k)是否满足利用电流突变量的防大数保护判据，其中，所述利用电流突变量的防大数保护判据的计算公式为：
|iZ(k)|＞iset4
式中，iset4为防大数突变量定值，按照第一直流输电线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，并在所述故障时计算第一直流输电线路防大数突变量|iz(k)|的值，并将所述|iz(k)|的值除以一个大于1的常数的结果作为iset4，t0≤k≤t。
9.一种基于纯电流特征的直流输电线路保护***，其特征在于，所述***包括：
初始设置单元，其用于确定第一直流输电线路和第二直流输电线路，以及第一直流输电线路的M侧和N侧，以及第二直流输电线路的M侧，设置启动保护判据、对端助增差分电流积聚量判据、对端助增突变量电流的积聚量高定值方向判据、对端助增突变量电流的积聚量低定值方向判据、第一极比幅式突变量电流选极判据、第二极比幅式突变量电流选极判据、比率制动电流选极判据、利用差分积聚量的防大数保护判据和利用电流突变量的防大数保护判据的计算公式，以及为上述判据中具备定值的参数和采样点数阈值赋值，其中，所述第一直流输电线路为输电线路中的任意一条，第二直流输电线路为除第一直流输电线路外的其他直流输电线路，第一直流输电线路M侧是所述第一直流输电线路整流侧和逆变侧中的一侧，第二直流输电线路M侧与第一直流输电线路M侧位于同一侧，所述第一直流输电线路N侧为所述直流输电线路中与M侧相对的另一侧；
数据采集单元，其用于实时采集第一直流输电线路M侧和N侧的线路电流，以及第二直流输电线路M侧的线路电流；
数据处理单元，其用于根据第一直流输电线路M侧的线路电流计算第一直流输电线路M侧电流的差分值，根据第一直流输电线路N侧的线路电流计算第一直流输电线路N侧电流的差分值，根据第二直流输电线路M侧的线路电流计算第二直流输电线路M侧电流的差分值，根据第一直流输电线路M侧电流的差分值计算第一直流输电线路M侧的电流突变量，根据第一直流输电线路N侧电流的差分值计算第一直流输电线路N侧的电流突变量，并根据所述计算确定的结果判定初始设置单元中的判据是否成立；
保护启动单元，其用于在保护启动判据成立时，启动第一直流输电线路M侧的直流保护,并将所述启动保护的时刻记为t0；
保护出口单元，其用于在启动第一直流输电线路M侧的直流保护后，当t时刻第一直流输电线路的对端助增差分电流积聚量判据、对端助增突变量电流的积聚量高定值方向判据和利用差分积聚量的防大数保护判据同时成立，在t0～t时间段满足利用电流突变量的防大数保护判据的采样点数大于预先设置的点数阈值，且满足利用电流突变量的防大数保护判据的采样点数大于不满足利用电流突变量的防大数保护判据的采样点数时，或者所述t时刻第一直流输电线路对端助增差分电流积聚量判据、对端助增突变量电流的积聚量低定值方向判据、阶段式电流选极判据和利用差分积聚量的防大数保护判据同时成立，在t0～t时间段满足利用电流突变量的防大数保护判据的采样点数大于预先设置的点数阈值，且满足利用电流突变量的防大数保护判据的采样点数大于不满足利用电流突变量的防大数保护判据的采样点数时，第一直流输电线路M侧保护动作出口，其中，当在t0～tlimit时间段内持续同时不满足第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据、在t0～tf时间段中任一个时刻满足第一极比幅式突变量电流选极判据或者tf时间后任一个时刻满足比率制动电流选极判据中的任意一个成立都表示阶段式电流选极判据成立，所述第一极是第一直流输电线路M侧，所述第二极是第二直流输电线路M侧，所述t为第一直流输电线路M侧的直流保护启动后，数据采集单元采样结束的时刻，tlimit和tf为保护启动后的两个时刻，tf＞tlimit。
10.根据权利要求9所述的***，其特征在于，所述数据处理单元根据第一直流输电线路M侧的线路电流计算第一直流输电线路M侧电流的差分值，根据第一直流输电线路N侧的线路电流计算第一直流输电线路N侧电流的差分值，根据第二直流输电线路M侧的线路电流计算第二直流输电线路M侧电流的差分值，根据第一直流输电线路M侧电流的差分值计算第一直流输电线路M侧的电流突变量，根据第一直流输电线路N侧电流的差分值计算第一直流输电线路N侧的电流突变量，其计算公式分别为：
diM(j)＝iM(j)-iM(j-ts)
式中，diM(j)为j时刻第一直流输电线路M侧电流的差分值，iM(j)为第一直流输电线路M侧在j时刻的线路电流，iM(j-ts)为j-ts时刻的线路电流，ts为采样时间间隔；
diN(j-Ttran)＝iN(j-Ttran)-iN(j-ts-Ttran)
式中，diN(j-Ttran)为j-Ttran时刻第一直流输电线路N侧电流的差分值，iN(j-Ttran)为第一直流输电线路N侧在j-Ttran时刻的线路电流，iN(j-ts-Ttran)为j-ts-Ttran时刻的线路电流，ts为采样时间间隔，Ttran为直流线路传输通道延时，t0≤j≤t；
diM′(j)＝iM′(j)-iM′(j-ts)
式中，diM'(j)为j时刻第二直流输电线路M侧电流的差分值，iM'(j)为第二直流输电线路M侧在j时刻的线路电流，iM'(j-ts)为j-ts时刻的线路电流，ts为采样时间间隔；
Figure FDA0003209259660000121
Figure FDA0003209259660000122
式中，△iM(k)为k时刻第一直流输电线路M侧电流突变量，diM(j)为j时刻第一直流输电线路M侧电流的差分值，△iN(k-Ttran)为k时刻第一直流输电线路N侧电流突变量△iN(k-Ttran)，diN(j-Ttran)为j时刻第一直流输电线路N侧电流的差分值，t0≤j≤k，t0≤k≤t，t0为保护启动时刻，ts为采样时间间隔，Ttran为直流线路传输通道延时。
11.根据权利要求9所述的***，其特征在于，所述初始设置单元设置的启动保护判据、对端助增差分电流积聚量判据、对端助增突变量电流的积聚量高定值方向判据、对端助增突变量电流的积聚量低定值方向判据、第一极比幅式突变量电流选极判据、第二极比幅式突变量电流选极判据、比率制动电流选极判据、利用差分积聚量的防大数保护判据和利用电流突变量的防大数保护判据的计算公式包括：
第一直流输电线路保护启动判据的计算公式为：
|diM(k0)|＞iset0
式中，|diM(k0)|为k0时刻第一直流输电线路M侧电流的差分值diM(k0)的绝对值，iset0为启动门槛，按照本线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，计算|diM(k0)|的值，利用所述|diM(k0)|的值再除以一个大于1的系数，将其结果作为启动门槛iset0；
第一直流输电线路对端助增差分电流积聚量判据的计算公式为：
iΣ(t)＞isetz
Figure FDA0003209259660000131
isetz＝kk(iset1+γ)
Figure FDA0003209259660000132
Figure FDA0003209259660000133
Figure FDA0003209259660000134
Figure FDA0003209259660000135
式中，iΣ(t)为t时刻第一直流输电线路差分积聚动作量，isetz为差分积聚门槛值，n为t0～t时间段的采样点数，|diM(k)|为k时刻第一直流输电线路M侧电流的差分值diM(k)的绝对值，kk为可靠系数，iset1为固定门槛，确定为区外金属性故障时的差分积聚动作量在故障后T0时间内的最大值，γ为负斜率电流浮动门槛，
Figure FDA0003209259660000141
为负斜率电流，T为浮动门槛计算窗长，ρ为浮动门槛系数，ρmax为比例常数，ρmax＞1，T'为定值计算窗长，T'＞T，λ为裕度系数，λ＞1；
第一直流输电线路对端助增突变量电流的积聚量高定值方向判据的计算公式为：
iΣ△(t)＞iset2H
Figure FDA0003209259660000142
式中，iΣ△(t)为t时刻第一直流输电线路方向动作量，iset2H为方向判别高定值，iset2H按照第二直流输电线路末端金属性故障整定，即在第一直流输电线路末端生成一个金属性故障，并在所述故障时计算第一直流输电线路方向动作量iΣ△(t)的值，并将所述iΣ△(t)的值乘以一个大于1的常数的结果作为iset2H，△iM(k)为k时刻第一直流输电线路M侧电流突变量，△iN(k-Ttran)为k-Ttran时刻第一直流输电线路N侧电流突变量，n为t0～t时间段的采样点数；
第一直流输电线路对端助增突变量电流的积聚量低定值方向判据的计算公式为：
iΣ△(t)＞iset2L
Figure FDA0003209259660000143
式中，iΣ△(t)为t时刻第一直流输电线路方向动作量，iset2L为方向判别低定值，iset2L按照第一直流输电线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，并在所述故障时计算第一直流输电线路方向动作量iΣ△(t)的值，并将所述iΣ△(t)的值除以一个大于1的常数的结果作为iset2L，△iM(k)为k时刻第一直流输电线路M侧电流突变量，△iN(k-Ttran)为k时刻第一直流输电线路N侧电流突变量，n为t0～t时间段的采样点数，iset2H＞iset2L；
第一极比幅式突变量电流选极判据的计算公式为：
Figure FDA0003209259660000151
第二极比幅式突变量电流选极判据的计算公式为：
Figure FDA0003209259660000152
第一直流输电线路比率制动电流选极判据的计算公式为：
iΣ△(t)＞αiD△(t)
Figure FDA0003209259660000153
式中，diM(j)为j时刻第一直流输电线路M侧电流的差分值，diM'(j)为j时刻第二直流输电线路M侧电流的差分值，σ为选极系数，σ＞1，iΣ(t)为t时刻第一直流输电线路差分积聚动作量，α为比率制动系数，α＞1，αiD△(t)为方向制动量，t0≤j≤k，t0≤k≤t，ts为采样时间间隔，Ttran为直流线路传输通道延时，当在t0～tf时间段内任一时刻第一极比幅式突变量电流选极判据成立，则阶段式电流选极判据成立，当在t0～tlimit时间段内第一极比幅式突变量电流选极判据和第二极比幅式突变量电流选极判据持续不成立时，在tlimit时刻及tlimit～tf时间段内判定发生双极故障，所述阶段式电流选极判据成立，保护启动tf时间后，以比率制动电流选极判据的结果作为阶段式电流选极判据的结果，比率制动电流选极判据成立则阶段式电流选极判据成立；
利用差分积聚量的防大数保护判据的计算公式为：
if(t)＞iset3
Figure FDA0003209259660000154
i'△(k)＝f(i△(k))
i△(k)＝|diM(k)|+△iN(k-Ttran)
式中，if(t)为防大数差分动作量，iset3为防大数差分定值，按照第一直流输电线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，并在所述故障时计算第一直流输电线路防大数差分动作量if(t)的值，并将所述if(t)的值除以一个大于1的常数的结果作为iset3，i△'(k)为k时刻第一直流输电线路M侧的置反差分量，f为置反处理方法，是将t0～t时间段绝对值最大的i△(k)强制乘以-ν作为k时刻的置反差分量，而将其他时刻的i△(k)作为对应时刻的置反差分量，ν＞1，i△(k)为k时刻第一直流输电线路M侧的差分量，diM(k)为k时刻第一直流输电线路M侧电流的差分值，△iN(k-Ttran)为k时刻第一直流输电线路N侧电流突变量，t0≤k≤t，ts为采样时间间隔，Ttran为直流线路传输通道延时；
第一直流输电线路利用电流突变量的防大数保护判据的计算公式为：
|iZ(k)|＞iset4
iZ(k)＝△iM(k)+△iN(k-Ttran)
式中，iz(k)为k时刻第一直流输电线路防大数突变量，iset4为防大数突变量定值，按照第一直流输电线路末端高阻故障有灵敏度整定，即在第一直流输电线路末端生成一个高阻故障，并在所述故障时计算第一直流输电线路防大数突变量iz(k)的值，并将所述iz(k)的值除以一个大于1的常数的结果作为iset4，△iM(k)为k时刻第一直流输电线路M侧电流突变量，△iN(k-Ttran)为第一直流输电线路k-Ttran时刻N侧电流突变量，t0≤k≤t。
CN201910599291.3A 2019-07-04 2019-07-04 一种基于纯电流特征的直流输电线路保护方法和*** Active CN110649577B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201910599291.3A CN110649577B (zh) 2019-07-04 2019-07-04 一种基于纯电流特征的直流输电线路保护方法和***
PCT/CN2020/096763 WO2021000738A1 (zh) 2019-07-04 2020-06-18 基于纯电流特征的直流输电线路保护方法和***
US17/596,826 US11973340B2 (en) 2019-07-04 2020-06-18 Direct-current transmission line protection method and system based on pure current characteristics

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910599291.3A CN110649577B (zh) 2019-07-04 2019-07-04 一种基于纯电流特征的直流输电线路保护方法和***

Publications (2)

Publication Number Publication Date
CN110649577A CN110649577A (zh) 2020-01-03
CN110649577B true CN110649577B (zh) 2021-10-01

Family

ID=69009400

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910599291.3A Active CN110649577B (zh) 2019-07-04 2019-07-04 一种基于纯电流特征的直流输电线路保护方法和***

Country Status (3)

Country Link
US (1) US11973340B2 (zh)
CN (1) CN110649577B (zh)
WO (1) WO2021000738A1 (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110649577B (zh) * 2019-07-04 2021-10-01 中国电力科学研究院有限公司 一种基于纯电流特征的直流输电线路保护方法和***
CN112564066B (zh) * 2020-12-24 2023-08-22 中国电力科学研究院有限公司 一种直流输电线路孪生电流差分保护方法及***
CN112564067B (zh) * 2020-12-24 2023-08-18 中国电力科学研究院有限公司 一种基于暂态电流变化的直流输电线路保护方法及***
CN113113899B (zh) * 2021-03-12 2022-06-14 南方电网科学研究院有限责任公司 一种直流输电线路保护方法、装置及存储介质
CN113488972B (zh) * 2021-08-27 2023-10-13 中国电力科学研究院有限公司 基于电流时域累积量的高压直流输电线路保护方法和***
CN114142442B (zh) * 2021-11-05 2022-06-28 昆明理工大学 一种直流输电线路保护方法及***
CN117277234B (zh) * 2023-08-29 2024-06-14 天津大学 适用于并联型多端高压直流输电***的单端量保护方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103597687A (zh) * 2011-04-04 2014-02-19 Abb技术有限公司 用于hvdc断路器的快速断路器故障检测
WO2014121438A1 (en) * 2013-02-05 2014-08-14 Alstom Technology Ltd. Method and apparatus for current differential protection for uhvdc transmission line
CN104113045A (zh) * 2014-04-29 2014-10-22 南京南瑞继保电气有限公司 Lc并联回路失谐故障的继电保护方法和装置
CN106208128A (zh) * 2016-08-24 2016-12-07 许继集团有限公司 一种混合三端高压直流输电***的功率反转方法
CN107112742A (zh) * 2014-03-21 2017-08-29 Abb瑞士股份有限公司 用于dc配电***的保护的方法和装置
CN107681644A (zh) * 2017-10-31 2018-02-09 华北电力大学（保定） 一种直流母线故障的比率制动式电流差动保护方法
CN109061397A (zh) * 2018-10-11 2018-12-21 南方电网科学研究院有限责任公司 一种线路故障区域识别方法
CN109787196A (zh) * 2019-01-09 2019-05-21 国网江苏省电力有限公司 低压直流配网保护控制方法和***

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009076769A1 (en) 2007-12-17 2009-06-25 The University Of British Columbia A method and system for protecting an electrical power transmission network
CN101662143B (zh) 2009-09-25 2012-04-25 国电南京自动化股份有限公司 一种防止在电流消失时误动的差动保护方法
CN101694936B (zh) 2009-10-22 2011-12-14 国电南京自动化股份有限公司 基于变化量带量方向的母线保护方法
CN106253244B (zh) 2016-08-10 2018-08-07 三峡大学 一种基于电流自构参考量的电流方向纵联保护方法
CN106684842B (zh) 2017-02-15 2018-08-21 东南大学 双端直流配电网保护方法
US10923907B2 (en) * 2018-03-07 2021-02-16 Southern States, Llc Direction-to-fault and zone-based distance-to-fault electric power sectionalizer systems
US11474139B2 (en) * 2019-04-05 2022-10-18 Schweitzer Engineering Laboratories, Inc. Fault direction calculation during current transformer saturation
CN110649577B (zh) 2019-07-04 2021-10-01 中国电力科学研究院有限公司 一种基于纯电流特征的直流输电线路保护方法和***
EP4027078A4 (en) * 2019-09-30 2023-10-25 Daikin Industries, Ltd. SOLID STATE COOLING MODULE
EP4030573A1 (en) * 2021-01-19 2022-07-20 Hitachi Energy Switzerland AG Device, system, and method for performing an online-update of a two-port equivalent

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103597687A (zh) * 2011-04-04 2014-02-19 Abb技术有限公司 用于hvdc断路器的快速断路器故障检测
WO2014121438A1 (en) * 2013-02-05 2014-08-14 Alstom Technology Ltd. Method and apparatus for current differential protection for uhvdc transmission line
CN107112742A (zh) * 2014-03-21 2017-08-29 Abb瑞士股份有限公司 用于dc配电***的保护的方法和装置
CN104113045A (zh) * 2014-04-29 2014-10-22 南京南瑞继保电气有限公司 Lc并联回路失谐故障的继电保护方法和装置
CN106208128A (zh) * 2016-08-24 2016-12-07 许继集团有限公司 一种混合三端高压直流输电***的功率反转方法
CN107681644A (zh) * 2017-10-31 2018-02-09 华北电力大学（保定） 一种直流母线故障的比率制动式电流差动保护方法
CN109061397A (zh) * 2018-10-11 2018-12-21 南方电网科学研究院有限责任公司 一种线路故障区域识别方法
CN109787196A (zh) * 2019-01-09 2019-05-21 国网江苏省电力有限公司 低压直流配网保护控制方法和***

Also Published As

Publication number Publication date
US11973340B2 (en) 2024-04-30
WO2021000738A1 (zh) 2021-01-07
CN110649577A (zh) 2020-01-03
US20220285931A1 (en) 2022-09-08

Similar Documents

Publication Publication Date Title
CN110649577B (zh) 一种基于纯电流特征的直流输电线路保护方法和***
DE69414328T2 (de) Lichtbogenerfassung über die Stromänderung
RU2422964C1 (ru) Устройство токовой защиты электрических сетей от однофазных замыканий на землю (варианты)
CN103872667B (zh) 一种防合并单元异常大数的线路电流差动保护方法
CN110336255B (zh) 一种基于电流积聚量的直流输电线路保护方法和***
CN111224386B (zh) 纵联保护方法、***及终端设备
CN110364997B (zh) 一种基于单端电流差分量的直流输电线路后备保护方法及***
CN110445103B (zh) 一种基于单端差分电流积聚量的直流输电线路保护方法和***
CN109217270B (zh) 一种直流线路分段式定时限过流保护方法与装置
EP2466710A1 (de) Verfahren und Vorrichtung zum Detektieren eines Störlichtbogenereignisses in einem Wechselspannungsnetz
CN110336258A (zh) 一种对端电流突变量助增的直流输电线路双端保护方法和***
CN109088425B (zh) 一种基于交流母线电压扰动量的换相失败预防方法
CN110336257B (zh) 一种防雷击和异常大数的直流输电线路保护方法和***
CN112564067B (zh) 一种基于暂态电流变化的直流输电线路保护方法及***
CN1614435A (zh) 电力***线路故障方向检测保护方法
CN110112715B (zh) 基于线路两侧电流采样值离散Fréchet距离的纵联保护方法
CN113488972B (zh) 基于电流时域累积量的高压直流输电线路保护方法和***
CN110429571B (zh) 一种可选极的直流输电线路保护方法和***
CN109103846B (zh) 一种基于比相和比幅原理的抽能电抗器保护方法及***
CN109752630B (zh) 一种同塔双回直流线路故障选线方法
CN113917366B (zh) 用于三相充电桩墙盒的相间短路检测电路及检测方法
DE2903533A1 (de) Einrichtung fuer den fehlerrichtungs- vergleichsschutz
CN110336256B (zh) 一种直流输电线路比率制动选极方法和***
RU162402U1 (ru) Устройство адаптивной резервной защиты трансформаторов на ответвлениях воздушной линии
CN112564065A (zh) 一种直流输电线路伪同步差动保护方法及***

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant