CN113820566B - 一种煤矿电网选漏方法 - Google Patents

一种煤矿电网选漏方法 Download PDF

Info

Publication number
CN113820566B
CN113820566B CN202111178390.8A CN202111178390A CN113820566B CN 113820566 B CN113820566 B CN 113820566B CN 202111178390 A CN202111178390 A CN 202111178390A CN 113820566 B CN113820566 B CN 113820566B
Authority
CN
China
Prior art keywords
line
transient zero
poles
iteration
coal mine
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
CN202111178390.8A
Other languages
English (en)
Other versions
CN113820566A (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.)
Xian University of Science and Technology
Original Assignee
Xian University of Science and Technology
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 Xian University of Science and Technology filed Critical Xian University of Science and Technology
Priority to CN202111178390.8A priority Critical patent/CN113820566B/zh
Publication of CN113820566A publication Critical patent/CN113820566A/zh
Application granted granted Critical
Publication of CN113820566B publication Critical patent/CN113820566B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08—Locating faults in cables, transmission lines, or networks
    • G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
    • G01R31/086—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08—Locating faults in cables, transmission lines, or networks
    • G01R31/088—Aspects of digital computing
    • 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/32—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 corresponding points in different conductors of a single system, e.g. of currents in go and return conductors

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Testing Of Short-Circuits, Discontinuities, Leakage, Or Incorrect Line Connections (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开了一种煤矿电网选漏方法，包括步骤：一、获取煤矿电网中各条线路的暂态零序电流；二、采用牛顿插值法求取各条线路暂态零序电流的包络线；三、采用可变遗忘因子递推最小二乘法进行各条线路暂态零序电流的包络线辨识，计算出煤矿电网中各条线路的拟合斜率；四、采用健全线路与故障线路包络线拟合斜率值正负相反的特点选择出故障线路。本发明能够在各种单相接地故障条件下准确选出故障线路，有利于提高漏电保护的灵敏性、可靠性和快速性，便于应用于实际漏电保护中。

Description

一种煤矿电网选漏方法
技术领域
本发明属于煤矿电网漏电保护技术领域，具体涉及一种煤矿电网选漏方法。
背景技术
我国是煤炭生产大国，煤炭生产为我国经济建设做出巨大贡献，因此需要保证煤煤矿下的安全。煤煤矿下空气潮湿，湿度普遍在95％左右，散热条件差，线路会受到井下矿石的挤压，使得对电气设备与电缆参数的要求更加苛刻，存在易爆的瓦斯煤尘且工作人员比较集中。在此恶劣环境下，电网的漏电有可能导致人身伤害、瓦斯煤尘爆炸和电雷管的先期爆发。因此，一旦漏电发生，保护装置必须迅速判断出漏电线路，动作于跳闸。这一要求比中压电网单相接地保护的要求严格，一般煤矿电网的漏电保护必须满足人身触电安全电流30mA·s的速动要求。
目前，应用于煤煤矿下的漏电保护主要有附加直流源式和零序功率方向保护。附加直流源式保护无选择性，它能快速地检测出漏电故障，但漏电时跳闸造成的停电范围大，从而造成更大的损失，且动作时间较长(一般大于50ms)。零序功率方向保护具有选择性，是利用基频零序电流电压的稳态特征保护法，稳态漏电信号一般比较微弱，又由于受到漏电过渡电阻和线路参数的影响，这些原本微弱的稳态信号幅值和相位的测量就更加困难，直接影响选漏的灵敏性和可靠性。现有技术还无法很好的保障煤矿电网供电的安全性、可靠性和稳定性。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种煤矿电网选漏方法，其能够在各种单相接地故障条件下准确选出故障线路，有利于提高漏电保护的灵敏性、可靠性和快速性，便于应用于实际漏电保护中。
为解决上述技术问题，本发明采用的技术方案是：一种煤矿电网选漏方法，该方法包括以下步骤：
步骤一、获取煤矿电网中各条线路的暂态零序电流；
步骤二、采用牛顿插值法求取各条线路暂态零序电流的包络线；
步骤三、采用可变遗忘因子递推最小二乘法进行各条线路暂态零序电流的包络线辨识，计算出煤矿电网中各条线路的拟合斜率；
步骤四、采用健全线路与故障线路包络线拟合斜率值正负相反的特点选择出故障线路。
上述的一种煤矿电网选漏方法，步骤一中所述获取煤矿电网中各条线路的暂态零序电流，是获取煤矿电网中各条线路
Figure BDA0003296318330000023
周波的暂态零序电流。
上述的一种煤矿电网选漏方法，步骤二中所述采用牛顿插值法求取各条线路暂态零序电流的包络线时，其中采用牛顿插值法求取任意一条线路暂态零序电流的包络线的具体过程为：
步骤201、对该条线路的暂态零序电流进行挑选，挑选出S个暂态零序电流的极值点，并将第s个暂态零序电流的极值点表示为i(ts)，s的取值为1～S的自然数；
步骤202、对S个暂态零序电流的极值点进行筛选，求取每个极值点两侧相邻两个极值的均值，并对该均值设置裕度d，按照中间极值点在该均值裕度d的范围内则保留的原则进行筛选,得到P个插值极点，并将第p个插值极点表示为i(tp)，p的取值为1～P的自然数；
步骤203、根据牛顿插值法得到插值极点关于采样时刻的各阶均差为：
Figure BDA0003296318330000021
Figure BDA0003296318330000022
Figure BDA0003296318330000031
其中，i[t1,t2]为插值极点关于采样时刻t1、t2的一阶均差，i(t1)为采样时刻t1的插值极点，i(t2)为采样时刻t2的插值极点；i[t1,t2,t3]为插值极点关于采样时刻t1、t2、t3的二阶均差，i[t2,t3]为插值极点关于采样时刻t2、t3的一阶均差，且
Figure BDA0003296318330000032
i(t3)为采样时刻t3的插值极点；i[t1,…,tP]为插值极点关于采样时刻t1、…、tP的P-1阶均差，i[t2,…,tP]为插值极点关于采样时刻t2、…、tP的P-2阶均差，且
Figure BDA0003296318330000033
i[t1,…,tP-1]为插值极点关于采样时刻t1、t2、…、tP-1的P-2阶均差且
Figure BDA0003296318330000034
i[tx,tx+1,…,tP]为插值极点关于采样时刻tx、tx+1、…、tP的P-x阶均差，且
Figure BDA0003296318330000035
i[tx,tx+1,…,tP-1]为插值极点关于采样时刻tx、tx+1、…、tP-1的P-1-x阶均差，且
Figure BDA0003296318330000036
i[t1,…,tP-y]为插值极点关于采样时刻t1、t2、…、tP-y的P-y-1阶均差且
Figure BDA0003296318330000037
x和y的取值为1～P-1的自然数，P的取值为大于等于2的自然数；
步骤204、根据插值极点关于采样时刻的各阶均差得到暂态零序电流包络线的表达式为：
i'(t)＝i(t1)+i[t1,t2](t-t1)+i[t1,t2,t3](t-t1)(t-t2)+…+i[t1,…,tP](t-t1)…(t-tP-1)(t-tP) (F4)。
上述的一种煤矿电网选漏方法，步骤三中所述采用可变遗忘因子递推最小二乘法进行各条线路暂态零序电流的包络线辨识，计算出煤矿电网中各条线路的拟合斜率的具体过程为：
步骤301、将可变遗忘因子的递推最小二乘法目标函数W表示为：
Figure BDA0003296318330000038
其中，L为总迭代次数，r为迭代次数且r的取值为1～L的自然数；tr′为第r次迭代时数据窗内的时间样本向量且t′r＝[t′r,1,t′r,2,…t′r,J]，t′r,j为第r次迭代时数据窗内的第j个时间样本，j的取值为1～J的自然数，J为第r次迭代时数据窗内的时间样本总数量，
Figure BDA0003296318330000041
为t'r的转置；i'(t′r)为第r次迭代时数据窗内的包络线曲线函数值向量；k′r为第r次迭代时M条线路被辨识的斜率参数向量且k′r＝[k′r,1,k′r,2,…,k′r,M]，k′r,m为第r次迭代时第m条线路被辨识的斜率参数，m的取值为1～M的自然数，M为线路总条数；λ为遗忘因子；每迭代一次数据窗向后推移一次；
具体实施时，J与M的取值相等；
步骤302、带可变遗忘因子的递推最小二乘法参数估计迭代公式为公式(F6)～(F9)：
Figure BDA0003296318330000042
Figure BDA0003296318330000043
Figure BDA0003296318330000044
Figure BDA0003296318330000045
其中，Qr为第r次迭代时的卡尔曼滤波增益矩阵，Pr为迭代至第r次时的协方差矩阵，Pr-1为迭代至第r-1次时的协方差矩阵，I为单位矩阵；k′r-1为第r-1次迭代时M条线路被辨识的斜率参数向量且k′r-1＝[k′r-1,1,k′r-1,2,…,k′r-1,M]，k′r-1,m为第r-1次迭代时第m条线路被辨识的斜率参数；er为先验误差，
Figure BDA0003296318330000046
为k′r-1的转置；
步骤303、在0<λ≤1的取值范围内，将λ的取值公式表示为：
Figure BDA0003296318330000047
Figure BDA0003296318330000048
Figure BDA0003296318330000049
Figure BDA00032963183300000410
Figure BDA00032963183300000411
Figure BDA00032963183300000412
其中，λmax为λ的最大取值，λmax＝1；ε(r)为后验误差，E{ε2(r)}为ε(r)的能量，
Figure BDA00032963183300000413
为系统噪声的功率；E{e2(r)}为er的能量，
Figure BDA00032963183300000414
为误差信号的功率；E{q2(r)}为q(r)的能量，q(r)为中间变量，
Figure BDA0003296318330000051
为系统功率；ξ为常数；
步骤304、给带可变遗忘因子的递推最小二乘法参数估计迭代公式赋初值，
Figure BDA0003296318330000052
e1＝0；根据初值和λ的取值公式对公式(F6)～(F9)进行迭代，得到煤矿电网中各条线路的拟合斜率。
上述的一种煤矿电网选漏方法，步骤四中所述采用健全线路与故障线路包络线拟合斜率值正负相反的特点选择出故障线路的具体过程为：
步骤401、将各条线路的拟合斜率构造出故障选线向量表示为：
K＝[k1,k2,…,kM] (F16)
其中，km为第m条线路的斜率参数，m的取值为1～M的自然数，M为线路总条数；
步骤402、当故障选线向量K中只有第f条线路对应的暂态零序电流包络线拟合斜率kf与其他线路的暂态零序电流包络线拟合斜率正负相反时，判断为第f条线路为故障线路；当所有线路的暂态零序电流包络线拟合斜率正负都相同时，判断为母线发生了单相接地故障。
本发明与现有技术相比具有以下优点：
1、本发明提出的煤矿电网选漏方法，对故障零序电流利用牛顿插值得到包络线，采用可变遗忘因子的递推最小二乘法计算包络线的拟合斜率作为特征值，再根据各条线路拟合斜率的正负判断故障线路，能够在各种单相接地故障条件下准确选出故障线路。
2、本发明提出的煤矿电网选漏方法，利用暂态过程特征选漏；煤矿电网发生漏电时，煤矿电网要经暂态过渡过程进入故障稳态过程，漏电暂态过程的信号由于含有比稳态分量大十几倍到几十倍的暂态分量，所以暂态过程的信号瞬时值很大；本发明利用暂态过程特征选漏，有利于提高漏电保护的灵敏性和可靠性，而且，利用暂态过程特征选漏，可能在暂态过程尚未结束之前就切除漏电故障，提高漏电保护的快速性。
3、本发明提出的煤矿电网选漏方法，能够解决已有方法在高阻及小故障合闸角条件下选线困难的问题，且方法实施简单，便于应用于实际漏电保护中。
综上所述，本发明能够在各种单相接地故障条件下准确选出故障线路，有利于提高漏电保护的灵敏性、可靠性和快速性，便于应用于实际漏电保护中。
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明的方法流程框图；
图2为本发明进行实验验证时故障选线实验系统的结构示意图。
具体实施方式
如图1所示，本发明的煤矿电网选漏方法，该方法包括以下步骤：
步骤一、获取煤矿电网中各条线路的暂态零序电流；
步骤二、采用牛顿插值法求取各条线路暂态零序电流的包络线；
步骤三、采用可变遗忘因子递推最小二乘法进行各条线路暂态零序电流的包络线辨识，计算出煤矿电网中各条线路的拟合斜率；
步骤四、采用健全线路与故障线路包络线拟合斜率值正负相反的特点选择出故障线路。
本实施例中，步骤一中所述获取煤矿电网中各条线路的暂态零序电流，是获取煤矿电网中各条线路
Figure BDA0003296318330000061
周波的暂态零序电流。
本实施例中，步骤二中所述采用牛顿插值法求取各条线路暂态零序电流的包络线时，其中采用牛顿插值法求取任意一条线路暂态零序电流的包络线的具体过程为：
步骤201、对该条线路的暂态零序电流进行挑选，挑选出S个暂态零序电流的极值点，并将第s个暂态零序电流的极值点表示为i(ts)，s的取值为1～S的自然数；
步骤202、对S个暂态零序电流的极值点进行筛选，求取每个极值点两侧相邻两个极值的均值，并对该均值设置裕度d，按照中间极值点在该均值裕度d的范围内则保留的原则进行筛选,得到P个插值极点，并将第p个插值极点表示为i(tp)，p的取值为1～P的自然数；
具体实施时，每条线路插值极点的个数P不一定相等。
步骤203、根据牛顿插值法得到插值极点关于采样时刻的各阶均差为：
Figure BDA0003296318330000071
Figure BDA0003296318330000072
Figure BDA0003296318330000073
其中，i[t1,t2]为插值极点关于采样时刻t1、t2的一阶均差，i(t1)为采样时刻t1的插值极点，i(t2)为采样时刻t2的插值极点；i[t1,t2,t3]为插值极点关于采样时刻t1、t2、t3的二阶均差，i[t2,t3]为插值极点关于采样时刻t2、t3的一阶均差，且
Figure BDA0003296318330000074
i(t3)为采样时刻t3的插值极点；i[t1,…,tP]为插值极点关于采样时刻t1、…、tP的P-1阶均差，i[t2,…,tP]为插值极点关于采样时刻t2、…、tP的P-2阶均差，且
Figure BDA0003296318330000075
i[t1,…,tP-1]为插值极点关于采样时刻t1、t2、…、tP-1的P-2阶均差且
Figure BDA0003296318330000076
i[tx,tx+1,…,tP]为插值极点关于采样时刻tx、tx+1、…、tP的P-x阶均差，且
Figure BDA0003296318330000077
i[tx,tx+1,…,tP-1]为插值极点关于采样时刻tx、tx+1、…、tP-1的P-1-x阶均差，且
Figure BDA0003296318330000078
i[t1,…,tP-y]为插值极点关于采样时刻t1、t2、…、tP-y的P-y-1阶均差且
Figure BDA0003296318330000079
x和y的取值为1～P-1的自然数，P的取值为大于等于2的自然数；
步骤204、根据插值极点关于采样时刻的各阶均差得到暂态零序电流包络线的表达式为：
i'(t)＝i(t1)+i[t1,t2](t-t1)+i[t1,t2,t3](t-t1)(t-t2)+…+i[t1,…,tP](t-t1)…(t-tP-1)(t-tP) (F4)。
本实施例中，步骤三中所述采用可变遗忘因子递推最小二乘法进行各条线路暂态零序电流的包络线辨识，计算出煤矿电网中各条线路的拟合斜率的具体过程为：
步骤301、将可变遗忘因子的递推最小二乘法目标函数W表示为：
Figure BDA0003296318330000081
其中，L为总迭代次数，r为迭代次数且r的取值为1～L的自然数；t′r为第r次迭代时数据窗内的时间样本向量且t′r＝[t′r,1,t′r,2,…t′r,J，t′r,j为第r次迭代时数据窗内的第j个时间样本，j的取值为1～J的自然数，J为第r次迭代时数据窗内的时间样本总数量，
Figure BDA0003296318330000082
为t'r的转置；i'(t′r)为第r次迭代时数据窗内的包络线曲线函数值向量；k′r为第r次迭代时M条线路被辨识的斜率参数向量且k′r＝[k′r,1,k′r,2,…,k′r,M]，k′r,m为第r次迭代时第m条线路被辨识的斜率参数，m的取值为1～M的自然数，M为线路总条数；λ为遗忘因子；每迭代一次数据窗向后推移一次；
具体实施时，J与M的取值相等；
步骤302、带可变遗忘因子的递推最小二乘法参数估计迭代公式为公式(F6)～(F9)：
Figure BDA0003296318330000083
Figure BDA0003296318330000084
Figure BDA0003296318330000085
Figure BDA0003296318330000086
其中，Qr为第r次迭代时的卡尔曼滤波增益矩阵，Pr为迭代至第r次时的协方差矩阵，Pr-1为迭代至第r-1次时的协方差矩阵，I为单位矩阵；k′r-1为第r-1次迭代时M条线路被辨识的斜率参数向量且k′r-1＝[k′r-1,1,k′r-1,2,…,kr′-1,M]，k′r-1,m为第r-1次迭代时第m条线路被辨识的斜率参数；er为先验误差，
Figure BDA0003296318330000087
为k′r-1的转置；
步骤303、在0<λ≤1的取值范围内，将λ的取值公式表示为：
Figure BDA0003296318330000091
Figure BDA0003296318330000092
Figure BDA0003296318330000093
Figure BDA0003296318330000094
Figure BDA0003296318330000095
Figure BDA0003296318330000096
其中，λmax为λ的最大取值，λmax＝1；ε(r)为后验误差，E{ε2(r)}为ε(r)的能量，
Figure BDA0003296318330000097
为系统噪声的功率；E{e2(r)}为er的能量，
Figure BDA0003296318330000098
为误差信号的功率；E{q2(r)}为q(r)的能量，q(r)为中间变量，
Figure BDA0003296318330000099
为系统功率；ξ为常数；
具体实施时，ξ为一个取值很小的常数，防止分母为零；
步骤304、给带可变遗忘因子的递推最小二乘法参数估计迭代公式赋初值，
Figure BDA00032963183300000910
e1＝0；根据初值和λ的取值公式对公式(F6)～(F9)进行迭代，得到煤矿电网中各条线路的拟合斜率。
具体实施时，λ的初值可以取为0.98；
本实施例中，步骤四中所述采用健全线路与故障线路包络线拟合斜率值正负相反的特点选择出故障线路的具体过程为：
步骤401、将各条线路的拟合斜率构造出故障选线向量表示为：
K＝[k1,k2,…,kM] (F16)
其中，km为第m条线路的斜率参数，m的取值为1～M的自然数，M为线路总条数；
步骤402、当故障选线向量K中只有第f条线路对应的暂态零序电流包络线拟合斜率kf与其他线路的暂态零序电流包络线拟合斜率正负相反时，判断为第f条线路为故障线路；当所有线路的暂态零序电流包络线拟合斜率正负都相同时，判断为母线发生了单相接地故障。
为了验证本发明的方法应用于煤矿电网漏电保护中的技术效果，进行了如下实验：
利用380V煤矿低压配电系统漏模拟平台，模拟一具有四条馈线路路的配电系统并设置单相接地故障，搭建的故障选线实验系统如图2所示；图2中，ES表示电源，Lk表示消弧线圈，TA1～TA4表示4个零序电流采样传感器，线路l1-l4对地电容分别为C1＝0.44μF、C2＝0.2μF、C3＝0.32μF、C4＝0.1μF；线路l1-l4的等效电感分别为L1＝6mH、L2＝6mH、L3＝6mH、L4＝6mH；线路l1-l4的等效电阻分别为R1＝2Ω、R2＝2Ω、R3＝2Ω、R4＝2Ω；得到的不同过渡电阻故障选漏结果如表1所示：
表1不同过渡电阻故障选漏结果(P＝10％,l3)
Figure BDA0003296318330000101
得到的不同运行方式故障选漏结果如表2所示：
表2不同运行方式故障选漏结果(Rg＝1000Ω,l3)
Figure BDA0003296318330000102
得到的不同故障位置故障选漏结果如表3所示：
表3不同故障位置故障选漏结果(Rg＝1000Ω,P＝10％)
Figure BDA0003296318330000103
表1和表3中，P为补偿度，Rg为故障过渡电阻；
从表1、表2和表3可以看出，本发明的煤矿电网选漏方法，能够在不同过渡电阻、不同运行方式、不同故障位置情况下，准确地进行选漏。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
前述对本发明的具体示例性实施方案的描述是为了说明和例证的目的。这些描述并非想将本发明限定为所公开的精确形式，并且很显然，根据上述教导，可以进行很多改变和变化。对示例性实施例进行选择和描述的目的在于解释本发明的特定原理及其实际应用，从而使得本领域的技术人员能够实现并利用本发明的各种不同的示例性实施方案以及各种不同的选择和改变。本发明的范围意在由权利要求书及其等同形式所限定。

Claims (3)

1.一种煤矿电网选漏方法，其特征在于，该方法包括以下步骤：
步骤一、获取煤矿电网中各条线路的暂态零序电流；
步骤二、采用牛顿插值法求取各条线路暂态零序电流的包络线；
步骤二中所述采用牛顿插值法求取各条线路暂态零序电流的包络线时，其中采用牛顿插值法求取任意一条线路暂态零序电流的包络线的具体过程为：
步骤201、对该条线路的暂态零序电流进行挑选，挑选出S个暂态零序电流的极值点，并将第s个暂态零序电流的极值点表示为i(ts)，s的取值为1～S的自然数；
步骤202、对S个暂态零序电流的极值点进行筛选，求取每个极值点两侧相邻两个极值的均值，并对该均值设置裕度d，按照中间极值点在该均值裕度d的范围内则保留的原则进行筛选,得到P个插值极点，并将第p个插值极点表示为i(tp)，p的取值为1～P的自然数；
步骤203、根据牛顿插值法得到插值极点关于采样时刻的各阶均差为：
Figure FDA0004137970900000011
Figure FDA0004137970900000012
Figure FDA0004137970900000013
其中，i[t1,t2]为插值极点关于采样时刻t1、t2的一阶均差，i(t1)为采样时刻t1的插值极点，i(t2)为采样时刻t2的插值极点；i[t1,t2,t3]为插值极点关于采样时刻t1、t2、t3的二阶均差，i[t2,t3]为插值极点关于采样时刻t2、t3的一阶均差，且
Figure FDA0004137970900000014
i(t3)为采样时刻t3的插值极点；i[t1,…,tP]为插值极点关于采样时刻t1、…、tP的P-1阶均差，i[t2,…,tP]为插值极点关于采样时刻t2、…、tP的P-2阶均差，且
Figure FDA0004137970900000015
i[t1,…,tP-1]为插值极点关于采样时刻t1、t2、…、tP-1的P-2阶均差且
Figure FDA0004137970900000016
i[tx,tx+1,…,tP]为插值极点关于采样时刻tx、tx+1、…、tP的P-x阶均差，且
Figure FDA0004137970900000021
i[tx,tx+1,…,tP-1]为插值极点关于采样时刻tx、tx+1、…、tP-1的P-1-x阶均差，且
Figure FDA0004137970900000022
i[t1,…,tP-y]为插值极点关于采样时刻t1、t2、…、tP-y的P-y-1阶均差且
Figure FDA0004137970900000023
x和y的取值为1～P-1的自然数，P的取值为大于等于2的自然数；
步骤204、根据插值极点关于采样时刻的各阶均差得到暂态零序电流包络线的表达式为：
i'(t)＝i(t1)+i[t1,t2](t-t1)+i[t1,t2,t3](t-t1)(t-t2)+…+i[t1,…,tP](t-t1)…(t-tP-1)(t-tP) (F4)；
步骤三、采用可变遗忘因子递推最小二乘法进行各条线路暂态零序电流的包络线辨识，计算出煤矿电网中各条线路的拟合斜率；
步骤三中所述采用可变遗忘因子递推最小二乘法进行各条线路暂态零序电流的包络线辨识，计算出煤矿电网中各条线路的拟合斜率的具体过程为：
步骤301、将可变遗忘因子的递推最小二乘法目标函数W表示为：
Figure FDA0004137970900000024
其中，L为总迭代次数，r为迭代次数且r的取值为1～L的自然数；t′r为第r次迭代时数据窗内的时间样本向量且t′r＝[t′r,1,t′r,2,…t′r,J]，t′r,j为第r次迭代时数据窗内的第j个时间样本，j的取值为1～J的自然数，J为第r次迭代时数据窗内的时间样本总数量，t'r T为t'r的转置；i'(t′r)为第r次迭代时数据窗内的包络线曲线函数值向量；k′r为第r次迭代时M条线路被辨识的斜率参数向量且k′r＝[k′r,1,k′r,2,…,kr,M]，k′r,m为第r次迭代时第m条线路被辨识的斜率参数，m的取值为1～M的自然数，M为线路总条数；λ为遗忘因子；每迭代一次数据窗向后推移一次；
步骤302、带可变遗忘因子的递推最小二乘法参数估计迭代公式为公式(F6)～(F9)：
Figure FDA0004137970900000031
Figure FDA0004137970900000032
k'r＝k'r-1+Qr[i′(t′r)-t'r Tk'r-1] (F8)
er＝i′(t′r)-k'r-1 Tt'r (F9)
其中，Qr为第r次迭代时的卡尔曼滤波增益矩阵，Pr为迭代至第r次时的协方差矩阵，Pr-1为迭代至第r-1次时的协方差矩阵，I为单位矩阵；k′r-1为第r-1次迭代时M条线路被辨识的斜率参数向量且k′r-1＝[k′r-1,1,k′r-1,2,…,k′r-1,M]，k′r-1,m为第r-1次迭代时第m条线路被辨识的斜率参数；er为先验误差，k'r-1 T为k′r-1的转置；
步骤303、在0<λ≤1的取值范围内，将λ的取值公式表示为：
Figure FDA0004137970900000033
ε(r)＝i′(tr)-k'r Tt'r (F11)
Figure FDA0004137970900000034
Figure FDA0004137970900000035
q(r)＝t'r TPrt'r (F14)
Figure FDA0004137970900000036
其中，λmax为λ的最大取值，λmax＝1；ε(r)为后验误差，E{ε2(r)}为ε(r)的能量，
Figure FDA0004137970900000037
为系统噪声的功率；E{e2(r)}为er的能量，
Figure FDA0004137970900000038
为误差信号的功率；E{q2(r)}为q(r)的能量，q(r)为中间变量，
Figure FDA0004137970900000039
为系统功率；ξ为常数；
步骤304、给带可变遗忘因子的递推最小二乘法参数估计迭代公式赋初值，P1＝(t'1 Tt'1)-1，e1＝0；根据初值和λ的取值公式对公式(F6)～(F9)进行迭代，得到煤矿电网中各条线路的拟合斜率；
步骤四、采用健全线路与故障线路包络线拟合斜率值正负相反的特点选择出故障线路。
2.按照权利要求1所述的一种煤矿电网选漏方法，其特征在于：步骤一中所述获取煤矿电网中各条线路的暂态零序电流，是获取煤矿电网中各条线路
Figure FDA0004137970900000041
周波的暂态零序电流。
3.按照权利要求1所述的一种煤矿电网选漏方法，其特征在于：步骤四中所述采用健全线路与故障线路包络线拟合斜率值正负相反的特点选择出故障线路的具体过程为：
步骤401、将各条线路的拟合斜率构造出故障选线向量表示为：
K＝[k1,k2,…,kM] (F16)
其中，km为第m条线路的斜率参数，m的取值为1～M的自然数，M为线路总条数；
步骤402、当故障选线向量K中只有第f条线路对应的暂态零序电流包络线拟合斜率kf与其他线路的暂态零序电流包络线拟合斜率正负相反时，判断为第f条线路为故障线路；当所有线路的暂态零序电流包络线拟合斜率正负都相同时，判断为母线发生了单相接地故障。
CN202111178390.8A 2021-10-10 2021-10-10 一种煤矿电网选漏方法 Active CN113820566B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111178390.8A CN113820566B (zh) 2021-10-10 2021-10-10 一种煤矿电网选漏方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111178390.8A CN113820566B (zh) 2021-10-10 2021-10-10 一种煤矿电网选漏方法

Publications (2)

Publication Number Publication Date
CN113820566A CN113820566A (zh) 2021-12-21
CN113820566B true CN113820566B (zh) 2023-05-09

Family

ID=78920007

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111178390.8A Active CN113820566B (zh) 2021-10-10 2021-10-10 一种煤矿电网选漏方法

Country Status (1)

Country Link
CN (1) CN113820566B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116466189A (zh) * 2023-05-09 2023-07-21 国网江苏省电力有限公司宿迁供电分公司 基于粒子群算法优化支持向量机的配电网故障选线方法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2764071A1 (fr) * 1997-06-03 1998-12-04 Electricite De France Procede d'identification de defauts monophases et dispositif pour la mise en oeuvre d'un tel procede
CN1604418A (zh) * 2004-11-04 2005-04-06 河北工业大学 漏电自动跟踪保护器及其运行方法
DE102009020908A1 (de) * 2008-10-17 2010-04-22 Edc Gmbh Erweitertes Verfahren zur Erkennung der Richtung eines Erdschlusses
CN102129012A (zh) * 2011-01-21 2011-07-20 昆明理工大学 一种利用形态峰谷检测的配电网故障选线方法
CN102832601A (zh) * 2012-08-29 2012-12-19 义马煤业集团股份有限公司 基于相电流分时比值的煤矿高压漏电保护方法
CN105606955A (zh) * 2015-12-04 2016-05-25 云南电网公司楚雄供电局 一种基于数值微分与经验模态分解的故障线路判别方法
CN106199342A (zh) * 2016-09-20 2016-12-07 西安科技大学 一种配电网单相接地故障选线方法
CN106370960A (zh) * 2016-08-15 2017-02-01 西安科技大学 基于故障瞬时正序电流实部分量的矿井电网漏电辨识方法
CN106597223A (zh) * 2016-12-19 2017-04-26 山东科技大学 一种基于改进hht的单相接地故障选线方法
RU2637378C1 (ru) * 2016-06-28 2017-12-04 Галина Андреевна Филатова Способ дистанционного определения места однофазного замыкания на землю
CN110879330A (zh) * 2019-12-02 2020-03-13 昆明理工大学 一种基于零序伏安曲线面积的配电网单相接地故障发展态势判别方法
CN110940894A (zh) * 2019-12-16 2020-03-31 国网北京市电力公司 一种基于泄漏电流瞬时幅相积和特征的触电故障检测方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11680977B2 (en) * 2020-03-18 2023-06-20 Mitsubishi Electric Research Laboratories, Inc. Transient based fault location method for ungrounded power distribution systems

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2764071A1 (fr) * 1997-06-03 1998-12-04 Electricite De France Procede d'identification de defauts monophases et dispositif pour la mise en oeuvre d'un tel procede
CN1604418A (zh) * 2004-11-04 2005-04-06 河北工业大学 漏电自动跟踪保护器及其运行方法
DE102009020908A1 (de) * 2008-10-17 2010-04-22 Edc Gmbh Erweitertes Verfahren zur Erkennung der Richtung eines Erdschlusses
CN102129012A (zh) * 2011-01-21 2011-07-20 昆明理工大学 一种利用形态峰谷检测的配电网故障选线方法
CN102832601A (zh) * 2012-08-29 2012-12-19 义马煤业集团股份有限公司 基于相电流分时比值的煤矿高压漏电保护方法
CN105606955A (zh) * 2015-12-04 2016-05-25 云南电网公司楚雄供电局 一种基于数值微分与经验模态分解的故障线路判别方法
RU2637378C1 (ru) * 2016-06-28 2017-12-04 Галина Андреевна Филатова Способ дистанционного определения места однофазного замыкания на землю
CN106370960A (zh) * 2016-08-15 2017-02-01 西安科技大学 基于故障瞬时正序电流实部分量的矿井电网漏电辨识方法
CN106199342A (zh) * 2016-09-20 2016-12-07 西安科技大学 一种配电网单相接地故障选线方法
CN106597223A (zh) * 2016-12-19 2017-04-26 山东科技大学 一种基于改进hht的单相接地故障选线方法
CN110879330A (zh) * 2019-12-02 2020-03-13 昆明理工大学 一种基于零序伏安曲线面积的配电网单相接地故障发展态势判别方法
CN110940894A (zh) * 2019-12-16 2020-03-31 国网北京市电力公司 一种基于泄漏电流瞬时幅相积和特征的触电故障检测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Shangbin Sun等.Leakage protection scheme of coal mine low voltage distribution network based on parameter identification. 《2011 International Conference on Electrical and Control Engineering》.2011,全文. *

Also Published As

Publication number Publication date
CN113820566A (zh) 2021-12-21

Similar Documents

Publication Publication Date Title
Kasztenny et al. Practical setting considerations for protective relays that use incremental quantities and traveling waves
CN111141995B (zh) 一种基于比幅原理的线路双端稳态量测距方法和系统
CN113820566B (zh) 一种煤矿电网选漏方法
Datsios et al. Implementation of leader development models in ATP-EMTP using a type-94 circuit component
CN105785229A (zh) 中性点不接地系统的故障识别法
Idris et al. Teaching distance relay using Matlab/Simulink graphical user interface
CN106487293A (zh) 一种非对称故障下双馈风机短路电流计算方法及装置
CN103543376A (zh) 用于小电流接地系统故障选线的径向基神经网络方法
Li et al. Protection challenges in future converter dominated power systems: Demonstration through simulation and hardware tests
Ji et al. Improved fault location through analysis of system parameters during autoreclose operations on transmission lines
Khaleghi et al. Transient fault area location and fault classification for distribution systems based on wavelet transform and adaptive neuro-fuzzy inference system (ANFIS)
Chaffey et al. Requirements for functional testing of HVDC protection IEDs
CN104135038A (zh) 一种交直流混联系统不对称故障分析方法
Almalki et al. Capacitor bank switching transient analysis using frequency dependent network equivalents
CN108539716B (zh) 一种基于暂态量的保护启动方法和装置
CN104852363B (zh) 一种基于故障录波数据的暂态串补电容电压计算方法
CN113437732A (zh) 一种光伏发电并网联络线纵联保护方法及系统
Hilbrich et al. Development of a hybrid platform for automated type and online application testing of protection & control schemes
Dahane et al. A novel algorithm for differential protection of untransposed transmission line using synchronized measurements
Lyamets et al. Description of 6–35 kV lines in recognition of earth fault in a distribution network
Garg et al. Fault classification, location in a series compensated power transmission network using Online Sequential Extreme Learning Machine
CN105447216A (zh) 一种提高电路仿真精度的方法及装置
CN104316791A (zh) 一种输电线路高阻接地故障原因分析查找方法
Zalitis et al. Experimental testing of distance protection performance in transient fault path resistance environment
Heo et al. Realization of distance relay algorithm using EMTP MODELS

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