CN115792507A - 一种基于短时窗斜率单调性的多重雷击判别方法及*** - Google Patents
一种基于短时窗斜率单调性的多重雷击判别方法及*** Download PDFInfo
- Publication number
- CN115792507A CN115792507A CN202310086723.7A CN202310086723A CN115792507A CN 115792507 A CN115792507 A CN 115792507A CN 202310086723 A CN202310086723 A CN 202310086723A CN 115792507 A CN115792507 A CN 115792507A
- Authority
- CN
- China
- Prior art keywords
- time window
- slope
- lightning stroke
- lightning
- judging
- 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.)
- Granted
Links
- 208000025274 Lightning injury Diseases 0.000 title claims abstract description 125
- 238000012850 discrimination method Methods 0.000 title description 5
- 238000000034 method Methods 0.000 claims abstract description 51
- 230000000630 rising effect Effects 0.000 claims abstract description 36
- 230000005540 biological transmission Effects 0.000 claims abstract description 29
- 230000003121 nonmonotonic effect Effects 0.000 claims abstract description 25
- 238000004364 calculation method Methods 0.000 claims abstract description 20
- 230000035772 mutation Effects 0.000 claims abstract description 16
- 238000002347 injection Methods 0.000 claims abstract description 8
- 239000007924 injection Substances 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims description 19
- 238000004458 analytical method Methods 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 6
- 238000012216 screening Methods 0.000 claims description 5
- 238000012544 monitoring process Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000000007 visual effect Effects 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 abstract description 2
- BULVZWIRKLYCBC-UHFFFAOYSA-N phorate Chemical compound CCOP(=S)(OCC)SCSCC BULVZWIRKLYCBC-UHFFFAOYSA-N 0.000 description 11
- 238000002474 experimental method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 4
- 238000009434 installation Methods 0.000 description 3
- 238000012806 monitoring device Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage or fault management, e.g. fault detection or location
Landscapes
- Locating Faults (AREA)
Abstract
本发明涉及一种基于短时窗斜率单调性的多重雷击判别方法及***。本发明包括步骤：当被保护输电线路遭受雷击产生故障时，选取雷击故障相电流波形的突变量；对所选取的突变量的上升沿处进行微分计算，求取相电流的短时窗斜率；判断各个时窗斜率的单调性；若短时窗斜率值满足单调变化，则判断为单次雷击；若短时窗斜率值满足存在非单调点，则判断该点为后续雷击注入时刻；统计录波时窗内所有非单调点个数记为n，得出在该段录波时窗内发生n+1重雷击。本发明大量不同幅值的雷电流与不同雷击时刻的仿真结果表明此方法有效且可靠，其物理概念直观清晰，实现方法简单，可广泛运用于输电线路的防雷保护。具有判断准确率高、速度快、适应性强的优点。
Description
技术领域
本发明涉及一种基于短时窗斜率单调性的多重雷击判别方法及***，属于新型电力***智能制造与继电保护技术领域。
背景技术
110kV架空输电线路作为国内电力***中的重要部分，确保电能稳定安全的输送具有重大意义。雷电流是一种超强的电磁脉冲电流，雷云放电时会在输电线路上形成暂态过电压，引起雷击跳闸事故。国内外的诸多学者主要的研究方向是基于首次雷击而引起的故障分量，对其进行计算分析，却忽略了后续雷击的影响。据统计数据分析，输电线路雷击跳闸等故障数占总故障数的60%以上，而根据雷电监测***的数据显示，70%以上的地面落雷均为多重雷击，但是由于缺少对应的监测装置，目前对后续雷电流的识别与分析的难度较大，实际操作相对困难，难以确定架空输电线路是由单重雷击还是多重雷击造成的。因此亟需一种简单可靠的方法，能够在雷击故障后精确地判别出是单重雷击还是多重雷击，为输电线路的雷击防护针对性分析起到可靠的基础性作用。
发明内容
本发明提供了一种基于短时窗斜率单调性的多重雷击判别方法及***，以用于从新型电力***智能制造与继电保护的角度出发，解决现有缺少对应的监测装置，目前对后续雷电流的识别与分析的难度较大，实际操作相对困难，难以确定架空输电线路是由单重雷击还是多重雷击造成的等不足之处。本发明基于保护装置对雷击故障识别的基础之上，提出通过在故障相电流波形的突变量处选取短时窗，计算短时窗斜率是否存在非单调点的方法，识别单重雷击和多重雷击，其原理清晰直观、所用电气量简单、信号处理过程简单、判别结果可靠，可为输电线路防雷研究提供简单易行的手段，为输电线路新型继电保护提供可靠的雷击辨识方法。
本发明的技术方案是：第一方面，本发明提供一种基于短时窗斜率单调性的多重雷击判别方法，本发明的识别方法总流程如图10所示，采集故障后输电线路三相电流信号，所述判别方法的具体步骤如下：
Step1、当被保护输电线路遭受雷击产生故障时，选取雷击故障相电流波形的突变量；
Step2、对上述所选取的突变量的上升沿处进行微分计算，求取相电流的短时窗斜率；
Step3、判断各个时窗斜率的单调性；
Step4、若短时窗斜率值满足单调变化，则判断为单次雷击；若短时窗斜率值满足存在非单调点，则判断该点为后续雷击注入时刻；
Step5、统计录波时窗内所有非单调点个数记为n，得出在该段录波时窗内发生n+1重雷击。
作为本发明的进一步方案，所述Step1中，选取雷击故障相电流波形的突变量，具体包括：
Step1.1、记时间t i 所对应的数据为data[t i ]，t i-1为t i 的上一个时刻，data[t i-1]为data[t i ]的上一个数据；
若满足：
则初步判断：故障相电流波形在时间t i 处发生突变；
Step1.2、记上述初步选取的突变量时刻为t i ，选取t i 之后第一个极值点；
Step1.3、记极值点所对应时刻为t im ；
所选取的突变量时刻t i 与极值点所对应时刻t im 之间应满足：
|t i －t im |＜t set
则得到突变量时刻t i ，t set 为整定值。
作为本发明的进一步方案，所述Step1.2中，选取t i 之后第一个极值点，具体包括：
Step1.2.1、记时间t m 所对应的数据为data[t m ]
若满足：
data[t m ]＞data[t m-1]
data[t m ]＞data[t m-2]
data[t m ]＞data[t m+1]
data[t m ]＞data[t m+2]
则初步判断：t m 为极值点所对应的时刻。
Step1.2.2、需要对初步所选极值进行优化处理，即并不是对所有的极值进行选取，需要进行优化筛选，需要剔除非保护线路的母线反射、折射波所带来的影响。
作为本发明的进一步方案，所述Step2中，对上述所选取的突变量的上升沿处进行微分计算，求取相电流的短时窗斜率，具体包括：
Step2.1、通过采用上述步骤得到优化筛选后的极值点，确定对应的波头；
Step2.2、选取由筛选过后的极值点所对应的突变量的波头，以在波前取达到峰值的10%的时刻t f 与在波尾取达到峰值的90%所对应的时刻t b ，将t f ～t b 取作短时窗；
Step2.3、记采样率为f s ，波头上升沿处的数据为data[t j ]，data[t j ]均为波头上升沿处的值；
通过计算短时窗内突变量的波头上升沿处的斜率：
选取其最大值，记为该短时窗的斜率值。
作为本发明的进一步方案，所述判断各个时窗斜率的单调性，具体包括：
记短时窗的斜率为k i (i=1,2,···) ；
若k i －k i+1＞0对所有i值均成立，则判断短时窗斜率单调变换，进而判断为单次雷击；
若k i －k i+1＞0存在i值不成立，则判断短时窗斜率非单调变换，进而判断为存在后续雷击注入时刻。
作为本发明的进一步方案，统计录波时窗内所有非单调点个数记为n，得出在该段录波时窗内发生n+1重雷击，具体包括：
记短时窗的斜率为k i (i=1,2,···)，在上述步骤中，若k i －k i+1＜0，则将第i+1个时窗所对应的值记作非单调点；
统计所有时窗内k i －k i+1＜0发生的次数，记为n，则在该段录波时窗内发生n+1重雷击。
第二方面，本发明提供一种基于短时窗斜率单调性的多重雷击判别***，包括：
信号采集与处理模块：用于采集线路的电流信号并进行处理与存储；具体包括：
电流变送单元、A/D转换单元、数据存储单元。
故障启动模块：用于比较电流变化率d i /d t 与整定值的大小，控制高频装置的启动；具体包括：
电流变化率监测单元、高频装置启动控制单元；
故障判别模块：用于判断故障是否为雷击故障；具体包括：
雷击故障判别单元：用于判断线路故障是否为雷击故障而非短路或断线等其他故障。
数据计算模块：对雷击故障相电流波形的突变量上升沿处进行微分计算，求取相电流的短时窗斜率；具体包括：
数据计算单元：对雷击故障相电流波形的突变量上升沿处进行微分计算，求取相电流的短时窗斜率；
图像处理分析与雷击识别模块：判断各时窗斜率的单调性，识别单次雷击与多重雷击。具体包括：
图像处理分析单元：判断各时窗斜率的单调性；
雷击类型识别单元：判断短时窗斜率值是否存在非单调点，识别单次雷击与多重雷击。
判别结果输出单元：用于输出雷击类型判别结果。
本发明的原理是，雷击会对输电线路产生较大的感应过电压和过电流，会引起线路跳闸等故障，根据雷击瞬时内在输电线路上雷击的次数分为单重雷击和多重雷击。单重雷击和多重雷击对输电线路的影响有显著区别，需要对其进行正确的识别。由于二者对输电线路雷击次数的不同，电流录波的传输形式也不同，两种情况下在站端采集的电流录波数据存在明显的区别；
由于能量守恒定律，输电线路在遭受单重雷击后，能量是逐渐消耗；而输电线路在遭受多重雷击后，能量在逐渐消耗时，后续的雷击注入到输电线路的能量会打破原先能量逐渐消耗的趋势，即会在后续的雷击注入时刻使得能量有一个显著的上升趋势。利用线路上安装的监测装置所采集的故障电流波形，将上述波形的突变量上升沿处微分后求取相电流的短时窗斜率，判断各时窗斜率是否存在非单调点，来判别出是单重雷击还是多重雷击。本发明可集成于保护***中，大量不同幅值的雷电流与不同雷击时刻的仿真结果表明此方法有效且可靠，其物理概念直观清晰，实现方法简单，可广泛运用于输电线路的防雷保护。具有判断准确率高、速度快、适应性强的优点。
本发明的有益效果是：
1、原理直观简单，物理概念清晰，易于实现，适合于实际工程应用；
2、仅处理和比较雷击后线路的故障相电流信息，无需考虑电流的其余信息以及其余电气量的信息，判据简单所以结果可靠；
3、信号处理简单，无需使用复杂算法，只需提取电流极值所在即可，信号处理过程简单因此出错概率极低。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对其所需要使用的附图作简单地介绍，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，可以根据这些附图获得所需的信息。
图1为交流***模型；
图2为杆塔模型；
图3为距首端20km处单次雷击干扰首端处测得电流原始波形情况；
图4为距首端20km处单次雷击干扰短时窗斜率值变化情况；
图5为距首端20km处多重雷击干扰短时窗斜率值变化情况；
图6为距首端20km处单次雷击故障短时窗斜率值变化情况；
图7为距首端20km处多重雷击故障短时窗斜率值变化情况；
图8为距首端100km处单次雷击干扰短时窗斜率值变化情况；
图9为距首端100km处多重雷击干扰短时窗斜率值变化情况；
图10为输电线路多重雷识别流程图。
具体实施方式
实施例1：如图1-图10所示，当输电线路发生雷击故障时，利用本发明所述方法和***可以准确识别是单重雷击与多重雷击，输电线路单重雷击与多重雷击识别流程如图10所示。第一方面，本发明提供一种基于短时窗斜率单调性的多重雷击判别方法，具体步骤如下：
Step1、当被保护输电线路遭受雷击产生故障时，选取雷击故障相电流波形的突变量；
Step2、对上述所选取的突变量的上升沿处进行微分计算，求取相电流的短时窗斜率；
Step3、判断各个时窗斜率的单调性；
Step4、若短时窗斜率值满足单调变化，则判断为单次雷击；若短时窗斜率值满足存在非单调点，则判断该点为后续雷击注入时刻；
Step5、统计录波时窗内所有非单调点个数记为n，得出在该段录波时窗内发生n+1重雷击。
作为本发明的进一步方案，所述Step1中，选取雷击故障相电流波形的突变量，具体包括：
Step1.1、记时间t i 所对应的数据为data[t i ]，t i-1为t i 的上一个时刻，data[t i-1]为data[t i ]的上一个数据；
若满足：
则初步判断：故障相电流波形在时间t i 处发生突变；
Step1.2、记上述初步选取的突变量时刻为t i ，选取t i 之后第一个极值点；
Step1.3、记极值点所对应时刻为t im ；
所选取的突变量时刻t i 与极值点所对应时刻t im 之间应满足
|t i －t im |＜t set
则得到突变量时刻t i ，t set 为整定值。
作为本发明的进一步方案，所述Step1.2中，选取t i 之后第一个极值点，具体包括：
Step1.2.1、记时间t m 所对应的数据为data[t m ]
若满足：
data[t m ]＞data[t m-1]
data[t m ]＞data[t m-2]
data[t m ]＞data[t m+1]
data[t m ]＞data[t m+2]
则初步判断：t m 为极值点所对应的时刻。
Step1.2.2、需要对初步所选极值进行优化处理，即并不是对所有的极值进行选取，需要进行优化筛选，需要剔除非保护线路的母线反射、折射波所带来的影响。
作为本发明的进一步方案，所述Step2中，对上述所选取的突变量的上升沿处进行微分计算，求取相电流的短时窗斜率，具体包括：
Step2.1、通过采用上述步骤得到优化筛选后的极值点，确定对应的波头；
Step2.2、选取由筛选过后的极值点所对应的突变量的波头，以在波前取达到峰值的10%的时刻t f 与在波尾取达到峰值的90%所对应的时刻t b ，将t f ～t b 取作短时窗；
Step2.3、记采样率为f s ，波头上升沿处的数据为data[t j ]，data[t j ]均为波头上升沿处的值；
通过计算短时窗内突变量的波头上升沿处的斜率k：
选取其最大值，记为该短时窗的斜率值。
作为本发明的进一步方案，所述判断各个时窗斜率的单调性，具体包括：
记短时窗的斜率为k i (i=1,2,···)；
若k i －k i+1＞0对所有i值均成立，则判断短时窗斜率单调变换，进而判断为单次雷击；
若k i －k i+1＞0存在i值不成立，则判断短时窗斜率非单调变换，进而判断为存在后续雷击注入时刻。
作为本发明的进一步方案，统计录波时窗内所有非单调点个数记为n，得出在该段录波时窗内发生n+1重雷击，具体包括：
记短时窗的斜率为k i (i=1,2,···)，在上述步骤中，若k i －k i+1＜0，则将第i+1个时窗所对应的值记作非单调点；
统计所有时窗内k i －k i+1＜0发生的次数，记为n，则在该段录波时窗内发生n+1重雷击。
第二方面，本发明提供一种基于短时窗斜率单调性的多重雷击判别***，包括：
信号采集与处理模块：用于采集线路的电流信号并进行处理与存储；具体包括：
电流变送单元、A/D转换单元、数据存储单元。
故障启动模块：用于比较电流变化率d i /d t 与整定值的大小，控制高频装置的启动；具体包括：
电流变化率监测单元、高频装置启动控制单元；
故障判别模块：用于判断故障是否为雷击故障；具体包括：
雷击故障判别单元：用于判断线路故障是否为雷击故障而非短路或断线等其他故障。
数据计算模块：对雷击故障相电流波形的突变量上升沿处进行微分计算，求取相电流的短时窗斜率；具体包括：
数据计算单元：对雷击故障相电流波形的突变量上升沿处进行微分计算，求取相电流的短时窗斜率；
图像处理分析与雷击识别模块：判断各时窗斜率的单调性，识别单次雷击与多重雷击。具体包括：
图像处理分析单元：判断各时窗斜率的单调性；
雷击类型识别单元：判断短时窗斜率值是否存在非单调点，识别单次雷击与多重雷击。
判别结果输出单元：用于输出雷击类型判别结果。
上述为本发明具体的识别步骤，实施例说明如下：本发明具体的识别步骤如上所述，按照图1所示交流***模型，通过仿真软件，以两组实验验证本发明的可靠性：
实验一：110kV线路，300km，在距首端20km处分别发生单重雷击幅值为30kA和多重雷击，其雷击幅值分别为30kA、8kA、5kA的雷击干扰，于站端的保护安装处测量所需电流行波并加以数据处理分析，判断雷击重数；
为保证实验的可靠性，这里直搭建一条母线，确保不存在非保护线路的母线反射、折射波所带来的影响，采样率为100kHz，t set =100μs。
由图3可知，在时间t 1.00036 时刻满足：
同理可得，在时间t 1.00236 ，时间t 1.00305 亦满足。
初步判断在时间t 1.00036 ，时间t 1.00236 ，时间t 1.00305 为突变点。
由图3可知，时间t 1.00039 所对应的数据data[t 1.00039 ]满足
data[t 1.00039 ]＞data[t 1.00038 ]
data[t 1.00039 ]＞data[t 1.00037 ]
data[t 1.00039 ]＞data[t 1.00040 ]
data[t 1.00039 ]＞data[t 1.00041 ]
可以判断：在时间t 1.00039 为时间t 1.00036 之后第一个极值点所对应时刻。
由图3可知，时间t 1.00036 之后第一个极值点所对应时刻为t 1.00039 ，满足：
|t 1.00039— t 1.00036 |＜100μs
可以判断时间t 1.00036 为突变时刻，时间t 1.00039 为到达极值点时刻；
同理可得，时间t 1.00236 为突变时刻，时间t 1.00241 为到达极值点时刻；
时间t 1.00305 为突变时刻，时间t 1.00309 为到达极值点时刻；
记第一个短时窗取作t 1.00035 ～t 1.00040 ，计算短时窗内突变量的波头上升沿处的斜率，最大值为600732.1；
记第二个短时窗取作t 1.00236 ～t 1.00244 ，计算短时窗内突变量的波头上升沿处的斜率，最大值为268187；
记第三个短时窗取作t 1.00305 ～t 1.00309 ，计算短时窗内突变量的波头上升沿处的斜率，最大值为96564；
由图4可知，单次雷击时，短时窗斜率值单调变化，判断为单重雷击，判断准确；
同理由图5可知，多重雷击时，短时窗斜率值变化存在非单调点，判断为多重雷击，且非单调点有两个，判断为三重雷击，判断准确；
实验二：110kV线路，300km，在距首端20km处分别发生单重雷击幅值为30kA和多重雷击，其雷击幅值分别为30kA、8kA的雷击故障，于站端的保护安装处测量所需电流行波并加以数据处理分析，判断雷击重数。
为保证实验的可靠性，这里直搭建一条母线，确保不存在非保护线路的母线反射、折射波所带来的影响，采样率为100kHz, t set =100μs。
通过上述步骤计算可得：
记第一个短时窗取作t 1.00006 ～t 1.00008 ，计算短时窗内突变量的波头上升沿处的斜率，最大值为738485；
记第二个短时窗取作t 1.00020 ～t 1.00024 ，计算短时窗内突变量的波头上升沿处的斜率，最大值为203872；
记第三个短时窗取作t 1.00209 ～t 1.00211 ，计算短时窗内突变量的波头上升沿处的斜率，最大值为110140；
记第四个短时窗取作t 1.00402 ～t 1.00405 ，计算短时窗内突变量的波头上升沿处的斜率，最大值为26609.4；
由图6可知，单次雷击时，短时窗斜率值单调变化，判断为单重雷击，判断准确；
同理，由图7可知，多重雷击时，短时窗斜率值变化存在非单调点，判断为多重雷击，且非单调点有一个，判断为二重雷击，判断准确；
实验三：110kV线路，300km，在距首端100km处分别发生单重雷击幅值为30kA和多重雷击，其雷击幅值分别为30kA、8kA、5kA的雷击干扰，于站端的保护安装处测量所需电流行波并加以数据处理分析，判断雷击重数；
为保证实验的可靠性，这里直搭建一条母线，确保不存在非保护线路的母线反射、折射波所带来的影响，采样率为100kHz, t set =100μs。
通过上述步骤计算可得：
记第一个短时窗取作t 1.00033 ～t 1.00035 ，计算短时窗内突变量的波头上升沿处的斜率，最大值为950570；
记第二个短时窗取作t 1.00169 ～t 1.00172 ，计算短时窗内突变量的波头上升沿处的斜率，最大值为258453；
记第三个短时窗取作t 1.00234 ～t 1.00238 ，计算短时窗内突变量的波头上升沿处的斜率，最大值为170914；
记第四个短时窗取作t 1.00369 ～t 1.00374 ，计算短时窗内突变量的波头上升沿处的斜率，最大值为107817；
记第五个短时窗取作t 1.00435 ～t 1.00440 ，计算短时窗内突变量的波头上升沿处的斜率，最大值为94684.9；
由图8可知，单次雷击时，短时窗斜率值单调变化，判断为单重雷击，判断准确；
同理，由图9可知，多重雷击时，短时窗斜率值变化存在非单调点，判断为多重雷击，且非单调点有两个，判断为三重雷击，判断准确。
上面结合附图对本发明的具体实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化。
Claims (10)
1.一种基于短时窗斜率单调性的多重雷击判别方法，其特征在于：采集故障后输电线路三相电流信号，所述判别方法的具体步骤如下：
Step1、当被保护输电线路遭受雷击产生故障时，选取雷击故障相电流波形的突变量；
Step2、对上述所选取的突变量的上升沿处进行微分计算，求取相电流的短时窗斜率；
Step3、判断各个时窗斜率的单调性；
Step4、若短时窗斜率值满足单调变化，则判断为单次雷击；若短时窗斜率值满足存在非单调点，则判断该点为后续雷击注入时刻；
Step5、统计录波时窗内所有非单调点个数记为n，得出在该段录波时窗内发生n+1重雷击。
2.根据权利要求1所述的基于短时窗斜率单调性的多重雷击判别方法，其特征在于：所述Step1中，选取雷击故障相电流波形的突变量，具体包括：
Step1.1、记时间
所对应的数据为
，
为
的上一个时刻，
为
的上一个数据；
若满足：
则初步判断：故障相电流波形在时间
处发生突变；
Step1.2、记上述初步选取的突变量时刻为
，选取
之后第一个极值点；
Step1.3、记极值点所对应时刻为
；
所选取的突变量时刻
与极值点所对应时刻
之间应满足
3.根据权利要求2所述的基于短时窗斜率单调性的多重雷击判别方法，其特征在于：所述Step1.2中，选取
之后第一个极值点，具体包括：
Step1.2.1、记时间
所对应的数据为
，若满足：
则初步判断：
为极值点所对应的时刻；
Step1.2.2、需要对初步所选极值进行优化处理，即并不是对所有的极值进行选取，需要进行优化筛选，需要剔除非保护线路的母线反射、折射波所带来的影响。
4.根据权利要求2所述的基于短时窗斜率单调性的多重雷击判别方法，其特征在于：所述Step2中，对上述所选取的突变量的上升沿处进行微分计算，求取相电流的短时窗斜率，具体包括：
Step2.1、通过采用上述步骤得到优化筛选后的极值点，确定对应的波头；
Step2.2、选取由筛选过后的极值点所对应的突变量的波头，以在波前取达到峰值的10%的时刻
与在波尾取达到峰值的90%所对应的时刻
，将
取作短时窗；
Step2.3、记采样率为
，波头上升沿处的数据为
，
均为波头上升沿处的值；
通过计算短时窗内突变量的波头上升沿处的斜率
：
选取其最大值，记为该短时窗的斜率值。
5.根据权利要求1所述的基于短时窗斜率单调性的多重雷击判别方法，其特征在于：所述判断各个时窗斜率的单调性，具体包括：
记短时窗的斜率为
；
若
对所有
值均成立，则判断短时窗斜率单调变换，进而判断为单次雷击；
若
存在
值不成立，则判断短时窗斜率非单调变换，进而判断为存在后续雷击注入时刻。
6.根据权利要求1所述的基于短时窗斜率单调性的多重雷击判别方法，其特征在于：统计录波时窗内所有非单调点个数记为n，得出在该段录波时窗内发生n+1重雷击，具体包括：
记短时窗的斜率为
，在上述步骤中，若
，则将第
个时窗所对应的值记作非单调点；
统计所有时窗内
发生的次数，记为n，则在该段录波时窗内发生n+1重雷击。
7.一种基于短时窗斜率单调性的多重雷击判别***，其特征在于，包括：
信号采集与处理模块：用于采集线路的电流信号并进行处理与存储；
故障启动模块：用于比较电流变化率
与整定值的大小，控制高频装置的启动；
故障判别模块：用于判断故障是否为雷击故障；
数据计算模块：对雷击故障相电流波形的突变量上升沿处进行微分计算，求取相电流的短时窗斜率；
图像处理分析与雷击识别模块：判断各时窗斜率的单调性，识别单次雷击与多重雷击。
8.根据权利要求7中所述的基于短时窗斜率单调性的多重雷击判别***，其特征在于：所述信号采集与处理模块具体包括：
电流变送单元、A/D转换单元、数据存储单元。
9.根据权利要求7中所述的基于短时窗斜率单调性的多重雷击判别***，其特征在于：所述故障启动模块具体包括：
电流变化率监测单元、高频装置启动控制单元；
所述故障判别模块具体包括：
雷击故障判别单元：用于判断线路故障是否为雷击故障而非短路或断线等其他故障。
10.根据权利要求7中所述的基于短时窗斜率单调性的多重雷击判别***，其特征在于：所述数值计算模块具体包括：
数据计算单元：对雷击故障相电流波形的突变量上升沿处进行微分计算，求取相电流的短时窗斜率；
所述图像处理分析与雷击识别模块具体包括：
图像处理分析单元：判断各时窗斜率的单调性；
雷击类型识别单元：判断短时窗斜率值是否存在非单调点，识别单次雷击与多重雷击；判别结果输出单元：用于输出雷击类型判别结果。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310086723.7A CN115792507B (zh) | 2023-02-09 | 2023-02-09 | 一种基于短时窗斜率单调性的多重雷击判别方法及*** |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310086723.7A CN115792507B (zh) | 2023-02-09 | 2023-02-09 | 一种基于短时窗斜率单调性的多重雷击判别方法及*** |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115792507A true CN115792507A (zh) | 2023-03-14 |
| CN115792507B CN115792507B (zh) | 2023-05-23 |
Family
ID=85430664
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310086723.7A Active CN115792507B (zh) | 2023-02-09 | 2023-02-09 | 一种基于短时窗斜率单调性的多重雷击判别方法及*** |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115792507B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117970037A (zh) * | 2024-03-29 | 2024-05-03 | 昆明理工大学 | 一种基于sod变化的输电线路多重雷识别方法 |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU95107817A (ru) * | 1995-05-15 | 1997-02-10 | Серпуховское высшее военное командно-инженерное училище ракетных войск им.Ленинского комсомола | Способ определения местоположения источника импульсного сигнала |
| CN103675602A (zh) * | 2013-12-03 | 2014-03-26 | 昆明理工大学 | 一种输电线路雷击闪络与未闪络的判别方法 |
| CN105318940A (zh) * | 2015-10-08 | 2016-02-10 | 华南理工大学 | 一种多路流量计检定装置脉冲计数信号重构方法 |
| CN107064677A (zh) * | 2017-04-05 | 2017-08-18 | 昆明理工大学 | 一种判别输电线路雷击闪络与未闪络的方法 |
| CN109471003A (zh) * | 2018-09-30 | 2019-03-15 | 昆明理工大学 | 一种基于110kV同塔双回输电线路雷击闪落与未闪落辨识方法 |
| CN110297146A (zh) * | 2019-07-30 | 2019-10-01 | 华北电力大学 | 基于暂态波形特征的输电线路雷击干扰与故障识别方法 |
| CN112415334A (zh) * | 2020-11-16 | 2021-02-26 | 国网上海市电力公司 | 一种多合一站接地***的数字孪生***及建模方法 |
| CN114520500A (zh) * | 2022-02-15 | 2022-05-20 | 山东大学 | 柔性直流电网输电线路保护及雷击识别方法及*** |
| CN114527332A (zh) * | 2022-01-21 | 2022-05-24 | 中国测试技术研究院电子研究所 | 雷击浪涌脉冲参数的自动测量方法、装置、介质 |
| WO2022236064A2 (en) * | 2021-05-06 | 2022-11-10 | Strong Force Iot Portfolio 2016, Llc | Quantum, biological, computer vision, and neural network systems for industrial internet of things |
| CN115616350A (zh) * | 2022-12-19 | 2023-01-17 | 昆明理工大学 | 一种基于曲率非单调性的雷击故障识别方法及*** |
-
2023
- 2023-02-09 CN CN202310086723.7A patent/CN115792507B/zh active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU95107817A (ru) * | 1995-05-15 | 1997-02-10 | Серпуховское высшее военное командно-инженерное училище ракетных войск им.Ленинского комсомола | Способ определения местоположения источника импульсного сигнала |
| CN103675602A (zh) * | 2013-12-03 | 2014-03-26 | 昆明理工大学 | 一种输电线路雷击闪络与未闪络的判别方法 |
| CN105318940A (zh) * | 2015-10-08 | 2016-02-10 | 华南理工大学 | 一种多路流量计检定装置脉冲计数信号重构方法 |
| CN107064677A (zh) * | 2017-04-05 | 2017-08-18 | 昆明理工大学 | 一种判别输电线路雷击闪络与未闪络的方法 |
| CN109471003A (zh) * | 2018-09-30 | 2019-03-15 | 昆明理工大学 | 一种基于110kV同塔双回输电线路雷击闪落与未闪落辨识方法 |
| CN110297146A (zh) * | 2019-07-30 | 2019-10-01 | 华北电力大学 | 基于暂态波形特征的输电线路雷击干扰与故障识别方法 |
| CN112415334A (zh) * | 2020-11-16 | 2021-02-26 | 国网上海市电力公司 | 一种多合一站接地***的数字孪生***及建模方法 |
| WO2022236064A2 (en) * | 2021-05-06 | 2022-11-10 | Strong Force Iot Portfolio 2016, Llc | Quantum, biological, computer vision, and neural network systems for industrial internet of things |
| CN114527332A (zh) * | 2022-01-21 | 2022-05-24 | 中国测试技术研究院电子研究所 | 雷击浪涌脉冲参数的自动测量方法、装置、介质 |
| CN114520500A (zh) * | 2022-02-15 | 2022-05-20 | 山东大学 | 柔性直流电网输电线路保护及雷击识别方法及*** |
| CN115616350A (zh) * | 2022-12-19 | 2023-01-17 | 昆明理工大学 | 一种基于曲率非单调性的雷击故障识别方法及*** |
Non-Patent Citations (3)
| Title |
|---|
| TAO LU 等: "Current distribution within a new meteorological gradient observation tower under direct lightning strike in Shenzhen and spatial electrical field around the tower" * |
| 向强铭: "基于奇异值分解理论的输电线路雷击故障测距的研究" * |
| 张广斌 等: "基于Hough变换直线检测的行波波头标定" * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117970037A (zh) * | 2024-03-29 | 2024-05-03 | 昆明理工大学 | 一种基于sod变化的输电线路多重雷识别方法 |
| CN117970037B (zh) * | 2024-03-29 | 2024-06-07 | 昆明理工大学 | 一种基于sod变化的输电线路多重雷识别方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115792507B (zh) | 2023-05-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104242267B (zh) | 一种风力发电送出输电线路距离保护方法 | |
| CN101345415B (zh) | 直流输电线路雷电绕击与反击分辨的行波分析识别方法 | |
| CN102590708B (zh) | 一种输电线路雷击故障与非雷击故障的辨别方法 | |
| CN103728538A (zh) | 一种小电流接地***的接地故障选线方法 | |
| CN103323747A (zh) | 一种小电流接地***的单相接地故障选线方法 | |
| CN103675602A (zh) | 一种输电线路雷击闪络与未闪络的判别方法 | |
| CN109375051A (zh) | 基于频谱密度衰减的雷电暂态信号识别方法及*** | |
| CN105138843A (zh) | 电力***采样飞点检测及其修复方法 | |
| CN115792507A (zh) | 一种基于短时窗斜率单调性的多重雷击判别方法及*** | |
| CN112083276A (zh) | 一种弧光接地故障的识别方法及*** | |
| CN105092997A (zh) | 一种特高压输电线路雷电绕击与反击的识别方法 | |
| CN100395932C (zh) | 一种暂态量保护的雷击故障识别方法 | |
| CN115616330B (zh) | 一种基于波形相似度的输电线路多重雷识别方法及*** | |
| CN108594156B (zh) | 一种改进的电流互感器饱和特性识别方法 | |
| CN103558484B (zh) | 一种变电站电缆受雷击干扰检测装置及方法 | |
| CN106771928A (zh) | 一种局部放电脉冲初至时刻在线拾取方法 | |
| CN115494350B (zh) | 一种交流输电线路雷击故障识别方法及*** | |
| CN115542087B (zh) | 一种适用于行波录波一体化设备的多重雷击识别方法 | |
| CN115616349B (zh) | 一种基于多元数据融合的多重雷击识别方法及*** | |
| CN117970037B (zh) | 一种基于sod变化的输电线路多重雷识别方法 | |
| CN114062880A (zh) | 基于自适应信号处理及cart树集成学习的直流故障电弧检测方法 | |
| CN114142442A (zh) | 一种直流输电线路保护方法及*** | |
| CN115575772B (zh) | 一种输电线路雷电绕击与反击的暂态零模量识别方法 | |
| CN118244041B (zh) | 一种高压直流输电线路雷电绕击与反击的识别方法与*** | |
| CN110943447A (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 |