CN112379216A - 输电线故障测距基波分量提取中量测误差抑制方法及系统 - Google Patents

输电线故障测距基波分量提取中量测误差抑制方法及系统 Download PDF

Info

Publication number
CN112379216A
CN112379216A CN202011225612.2A CN202011225612A CN112379216A CN 112379216 A CN112379216 A CN 112379216A CN 202011225612 A CN202011225612 A CN 202011225612A CN 112379216 A CN112379216 A CN 112379216A
Authority
CN
China
Prior art keywords
fundamental component
fundamental
data
amplitude
phase
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.)
Pending
Application number
CN202011225612.2A
Other languages
English (en)
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.)
Shandong University
Original Assignee
Shandong University
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 Shandong University filed Critical Shandong University
Priority to CN202011225612.2A priority Critical patent/CN112379216A/zh
Priority to PCT/CN2020/133638 priority patent/WO2022095185A1/zh
Publication of CN112379216A publication Critical patent/CN112379216A/zh
Pending legal-status Critical Current

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/085—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution lines, e.g. overhead
    • 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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Locating Faults (AREA)

Abstract

本发明公开了一种输电线故障测距基波分量提取中量测误差抑制方法及系统，包括：通过滑动采样数据窗口的起始时刻形成多组采样数据；使用离散傅里叶变换从每组采样数据中分别提取基波分量；对基波分量的幅值和初相分别进行线性拟合，得到准确的幅值和相位。本发明有益效果：可以直接对采样数据进行滤波处理，得到更加精准的基波分量计算结果，从而提高故障测距结果的精度。

Description

输电线故障测距基波分量提取中量测误差抑制方法及系统
技术领域
本发明涉及输电线故障测距技术领域，尤其涉及一种输电线故障测距基波分量提取中量测误差抑制方法及系统。
背景技术
本部分的陈述仅仅是提供了与本发明相关的背景技术信息，不必然构成在先技术。
现有基于量测阻抗的故障测距算法都需要利用故障后电压和电流的基波分量进行计算。基波分量的准确性会极大的影响故障测距精度。而量测误差的存在，会使得用于测距计算的电压电流基波分量存在误差，影响测距结果的精度。
为克服量测误差对测距精度的影响，现有方法一般是通过滑动采样数据窗提取多组基波分量构建冗余测距方程，利用最小二乘法提升测距精度。虽然该方法取得了一定效果，但对于配电网中较为常见的T接线或多分支接线而言，求解的测距方程比较复杂，利用冗余数据增加方程组数较多，不利于非线性测距方程的数值求解。
发明内容
为了解决上述问题，本发明提出了一种输电线故障测距基波分量提取中量测误差抑制方法及系统，可以直接对采样数据进行滤波处理，得到更加精准的基波分量计算结果，从而提高故障测距结果的精度。
在一些实施方式中，采用如下技术方案：
一种输电线故障测距基波分量提取中量测误差抑制方法，包括：
通过滑动采样数据窗口的起始时刻形成多组采样数据；
使用离散傅里叶变换从每组采样数据中分别提取基波分量；
对基波分量的幅值和初相分别进行线性拟合，得到准确的幅值和相位。
在另一些实施方式中，采用如下技术方案：
一种输电线故障测距基波分量提取中量测误差抑制系统，包括：
数据采样模块，用于通过滑动采样数据窗口的起始时刻形成多组采样数据；
基波分量提取模块，用于使用离散傅里叶变换从每组采样数据中分别提取基波分量；
线性拟合模块，用于对基波分量的幅值和初相分别进行线性拟合，得到准确的幅值和相位。
在另一些实施方式中，采用如下技术方案：
一种终端设备，其包括处理器和计算机可读存储介质，处理器用于实现各指令；计算机可读存储介质用于存储多条指令，所述指令适于由处理器加载并执行上述的输电线故障测距基波分量提取中量测误差抑制方法。
在另一些实施方式中，采用如下技术方案：
一种计算机可读存储介质，其中存储有多条指令，所述指令适于由终端设备的处理器加载并执行上述的输电线故障测距基波分量提取中量测误差抑制方法。
与现有技术相比，本发明的有益效果是：
本发明根据故障后电压和电流基波分量幅值和初相位恒定不变特征，提出一种新的量测误差抑制方法，可以直接对采样数据进行滤波处理，得到更加精准的基波分量计算结果，从而提高故障测距结果的精度。
本发明的其他特征和附加方面的优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本方面的实践了解到。
附图说明
图1为本发明实施例中输电线故障测距基波分量提取中量测误差抑制方法流程图。
具体实施方式
应该指出，以下详细说明都是例示性的，旨在对本申请提供进一步的说明。除非另有指明，本发明使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
实施例一
在一个或多个实施方式中，公开了一种输电线故障测距基波分量提取中量测误差抑制方法，参照图1，具体包括如下步骤：
(1)通过滑动采样数据窗口的起始时刻形成多组采样数据；
(2)使用离散傅里叶变换从每组采样数据中分别提取基波分量；
(3)对基波分量的幅值和初相分别进行线性拟合，得到准确的幅值和相位。
具体地，假设故障后基波分量表达式为：
u1(t)＝U1cos(w1t+θ0) (1)
式中，w1＝2πf1为基波角频率，f1为基波分量的频率，θ0为基波分量的初相位，简称初相，U1为基波分量的幅值。根据基波分量的表达式可知其幅值和初相不随故障后时间推移变化，是常量。相位
Figure BDA0002763561140000041
则是随时间变化的线性函数。
因此，利用故障后工频一周波采样点进行基波分量提取，如果结果准确，理论上其幅值不随采样点数据窗起点时刻变化而变化。
假设故障起始时刻为0时刻，采样周期为Δt，如按Δt步长依次平移提取基波分量的一周波采样点数据窗，对于起始时刻为kΔt的数据窗，则相位是步长Δt的线性函数
Figure BDA0002763561140000042
起始时刻的相角
Figure BDA0002763561140000043
根据上述基波分量的幅值和初相恒定不变的特征规律，本实施例选用了线性拟合滤波，通过滑动数据窗口起始时刻形成多组采样数据，使用离散傅里叶变换(DFT)从每个数据窗口分别提取基波分量，并求出幅值和初相，这些基波分量的幅值和初相围绕对应的准确值上下波动，进而可以借助线性拟合获得准确的幅值和相位。
下面以故障后量测点m处所得AB相间线电压基波分量的幅值Umab和初相θ0为例进行具体阐述。
在故障后选取32个连续的采样数据窗口，每个数据窗的长度是基波分量的一个周期，使用离散傅里叶变换(DFT)从每个采样数据窗口内提取基波分量并将幅值记为Umab(i)(i＝0,1,2,……,31)，待求的线性拟合结果为Umab，基于最小二乘法的线性拟合是指拟合一个近似直线y＝Umab,使得对于所有数据来说总的误差和err1最小，即求：
Figure BDA0002763561140000044
为使误差和最小，对误差函数求导，并且令导函数等于0：
Figure BDA0002763561140000051
进而有：
Figure BDA0002763561140000052
Figure BDA0002763561140000053
再对初相进行拟合，由于在数据窗内相角随时间一次线性变化，所以相角满足关系式
Figure BDA0002763561140000054
基波分量频率f1＝50Hz，w1＝2πf1＝100πrad/s,则
Figure BDA0002763561140000055
由于初相θ0与时间无关，故是个恒定不变的常数。
将32个数据窗中提取的基波分量的相角记为为αk(k＝0,1,2……,31)，则计算出的初相θk＝αk-w1kΔt(k＝0,1,2……,31),待求的线性拟合结果为θ0，基于最小二乘法的线性拟合是指拟合一个近似直线y＝θ0，使得对于所有θk，总的误差和err2最小，即求：
Figure BDA0002763561140000056
为使误差和最小，误差函数对θ0求导，并且令导函数等于0：
Figure BDA0002763561140000057
进而解出：
Figure BDA0002763561140000058
为了验证所提出的量测误差抑制方法的自洽性和有效性，在MATLAB中对进行仿真验证，首先构造信号
Figure BDA0002763561140000059
式(9)中基波分量的幅值为100，初相位为
Figure BDA00027635611400000510
再给信号添加20组10％的随机噪音，分别采用滑动数据窗提取基波分量并线性拟合与从t＝0初始时刻直接提取基波分量两种方法获得基波分量幅值和初相位，并将结果进行比较如表1所示，经过线性拟合提取的基波分量平均幅值为99.8412，相对误差为0.15％，最大相对误差为0.63％；未经过线性拟合直接提取的基波分量平均幅值为97.47382，相对误差为2.53％，最大相对误差为4.59％。经过线性拟合提取的基波分量平均相角为59.9536°，相对误差为0.08％，最大相对误差为1.06％；未经过线性拟合提取的基波分量平均相角为58.2766，相对误差为2.87％，最大相对误差为6.37％。数据表明本发明所提基于线性拟合的量测噪声抑制方法能够很好地抑制随机噪音的干扰，基波分量幅值和初相计算结果相比直接提取结果更接近基波分量的准确值。
表1 20组含随机噪声的线性拟合结果和未经线性拟合基本分量提取结果比较
Figure BDA0002763561140000061
Figure BDA0002763561140000071
实施例三
在一个或多个实施方式中，公开了一种主端可穿戴遥操作控制装置的遥操作控制器，包括服务器，所述服务器包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现实施例一中的方法。为了简洁，在此不再赘述。
应理解，本实施例中，处理器可以是中央处理单元CPU，处理器还可以是其他通用处理器、数字信号处理器DSP、专用集成电路ASIC，现成可编程门阵列FPGA或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
存储器可以包括只读存储器和随机存取存储器，并向处理器提供指令和数据、存储器的一部分还可以包括非易失性随机存储器。例如，存储器还可以存储设备类型的信息。
在实现过程中，上述方法的各步骤可以通过处理器中的硬件的集成逻辑电路或者软件形式的指令完成。
实施例一中的方法可以直接体现为硬件处理器执行完成，或者用处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器、闪存、只读存储器、可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器，处理器读取存储器中的信息，结合其硬件完成上述方法的步骤。为避免重复，这里不再详细描述。
实施例四
在一个或多个实施方式中，公开了一种计算机可读存储介质，其中存储有多条指令，所述指令适于由终端设备的处理器加载并实施例一中所述的一种输电线故障测距基波分量提取中量测误差抑制方法。
上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。

Claims (10)

1.一种输电线故障测距基波分量提取中量测误差抑制方法，其特征在于，包括：
通过滑动采样数据窗口的起始时刻形成多组采样数据；
使用离散傅里叶变换从每组采样数据中分别提取基波分量；
对基波分量的幅值和初相分别进行线性拟合，得到准确的幅值和相位。
2.如权利要求1所述的一种输电线故障测距基波分量提取中量测误差抑制方法，其特征在于，通过滑动采样数据窗口的起始时刻形成多组采样数据，具体包括：
假设故障起始时刻为0时刻，采样周期为Δt，如按Δt步长依次平移提取基波分量的一周波采样点数据窗，对于起始时刻为kΔt的数据窗，则相位是步长Δt的线性函数
Figure FDA0002763561130000011
起始时刻的相角
Figure FDA0002763561130000012
其中，w1为基波角频率，θ0为基波分量的初相位，简称初相；k为从故障初始时刻到当前数据窗起点所需平移步长的个数。
3.如权利要求1所述的一种输电线故障测距基波分量提取中量测误差抑制方法，其特征在于，所述基波分量的幅值和初相不随故障后时间推移而变化，为常量。
4.如权利要求1所述的一种输电线故障测距基波分量提取中量测误差抑制方法，其特征在于，使用离散傅里叶变换从每组采样数据中分别提取基波分量并将幅值记为Umab(i)，i＝0,1,2,……,N-1，基于最小二乘法拟合一个近似直线y＝Umab,使得对于所有数据来说总的误差和err1最小，从而求得线性拟合结果Umab。
5.如权利要求4所述的一种输电线故障测距基波分量提取中量测误差抑制方法，其特征在于，所述线性拟合结果Umab具体为：
Figure FDA0002763561130000021
其中，N为采样数据窗口的个数。
6.如权利要求1所述的一种输电线故障测距基波分量提取中量测误差抑制方法，其特征在于，使用离散傅里叶变换从每组采样数据中分别提取基波分量并将相角记为αk，k＝0,1,2……,N-1，分别计算出初相θk，基于最小二乘法拟合一个近似直线y＝θ0,使得对于所有数据来说总的误差和err2最小，从而求得线性拟合结果θ0。
7.如权利要求6所述的一种输电线故障测距基波分量提取中量测误差抑制方法，其特征在于，所述线性拟合结果θ0具体为：
Figure FDA0002763561130000022
其中，N为采样数据窗口的个数。
8.一种输电线故障测距基波分量提取中量测误差抑制系统，其特征在于，包括：
数据采样模块，用于通过滑动采样数据窗口的起始时刻形成多组采样数据；
基波分量提取模块，用于使用离散傅里叶变换从每组采样数据中分别提取基波分量；
线性拟合模块，用于对基波分量的幅值和初相分别进行线性拟合，得到准确的幅值和相位。
9.一种终端设备，其包括处理器和计算机可读存储介质，处理器用于实现各指令；计算机可读存储介质用于存储多条指令，其特征在于，所述指令适于由处理器加载并执行权利要求1-7任一项所述的输电线故障测距基波分量提取中量测误差抑制方法。
10.一种计算机可读存储介质，其中存储有多条指令，其特征在于，所述指令适于由终端设备的处理器加载并执行权利要求1-7任一项所述的输电线故障测距基波分量提取中量测误差抑制方法。
CN202011225612.2A 2020-11-05 2020-11-05 输电线故障测距基波分量提取中量测误差抑制方法及系统 Pending CN112379216A (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202011225612.2A CN112379216A (zh) 2020-11-05 2020-11-05 输电线故障测距基波分量提取中量测误差抑制方法及系统
PCT/CN2020/133638 WO2022095185A1 (zh) 2020-11-05 2020-12-03 输电线故障测距基波分量提取中量测误差抑制方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011225612.2A CN112379216A (zh) 2020-11-05 2020-11-05 输电线故障测距基波分量提取中量测误差抑制方法及系统

Publications (1)

Publication Number Publication Date
CN112379216A true CN112379216A (zh) 2021-02-19

Family

ID=74579377

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011225612.2A Pending CN112379216A (zh) 2020-11-05 2020-11-05 输电线故障测距基波分量提取中量测误差抑制方法及系统

Country Status (2)

Country Link
CN (1) CN112379216A (zh)
WO (1) WO2022095185A1 (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115277294B (zh) * 2022-07-27 2024-03-26 中车青岛四方车辆研究所有限公司 一种差分总线在线健康诊断方法及系统
CN117471366B (zh) * 2023-12-27 2024-04-09 国网福建省电力有限公司 一种基于最小二乘法的计量装置中性线接触不良研判方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3812433A1 (de) * 1988-04-14 1989-10-26 Siemens Ag Verfahren zur bestimmung des fehlerortes auf einer elektrischen leitung
CN101441233A (zh) * 2008-12-02 2009-05-27 湖南海兴电器有限责任公司 基于Kaiser窗双谱线插值FFT的基波与谐波检测方法
CN101587147A (zh) * 2009-06-25 2009-11-25 中国电力科学研究院 一种同步相量测量装置进行相量修正的方法
CN106501602A (zh) * 2016-09-28 2017-03-15 西南交通大学 一种基于滑窗频谱分离的基波参数测量方法
CN109490625A (zh) * 2018-11-14 2019-03-19 国网江苏省电力有限公司徐州供电分公司 一种基于滑动窗口和半定规划的谐波信号分析方法
CN109831184A (zh) * 2019-02-18 2019-05-31 山东大学 一种带通滤波方法、系统、存储介质及终端设备
CN110244144A (zh) * 2019-05-28 2019-09-17 许昌许继软件技术有限公司 一种避雷器状态监测方法及交流数据采集方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8680872B2 (en) * 2009-09-29 2014-03-25 International Business Machines Corporation Identification of false positives in high impedance fault detection
CN109256762B (zh) * 2017-07-13 2021-10-22 许继集团有限公司 配电网单相故障电流补偿控制方法、装置与系统
CN109085455A (zh) * 2017-12-26 2018-12-25 贵州电网有限责任公司 一种用于配电线路高阻接地故障的判定方法
CN109145476B (zh) * 2018-08-31 2023-05-26 重庆水利电力职业技术学院 用于电力系统信号处理的时域自适应分段复指数级数法
CN109802372B (zh) * 2019-03-05 2020-06-02 华北电力大学 基于暂态电流波形特征的快速差动保护方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3812433A1 (de) * 1988-04-14 1989-10-26 Siemens Ag Verfahren zur bestimmung des fehlerortes auf einer elektrischen leitung
CN101441233A (zh) * 2008-12-02 2009-05-27 湖南海兴电器有限责任公司 基于Kaiser窗双谱线插值FFT的基波与谐波检测方法
CN101587147A (zh) * 2009-06-25 2009-11-25 中国电力科学研究院 一种同步相量测量装置进行相量修正的方法
CN106501602A (zh) * 2016-09-28 2017-03-15 西南交通大学 一种基于滑窗频谱分离的基波参数测量方法
CN109490625A (zh) * 2018-11-14 2019-03-19 国网江苏省电力有限公司徐州供电分公司 一种基于滑动窗口和半定规划的谐波信号分析方法
CN109831184A (zh) * 2019-02-18 2019-05-31 山东大学 一种带通滤波方法、系统、存储介质及终端设备
CN110244144A (zh) * 2019-05-28 2019-09-17 许昌许继软件技术有限公司 一种避雷器状态监测方法及交流数据采集方法

Also Published As

Publication number Publication date
WO2022095185A1 (zh) 2022-05-12

Similar Documents

Publication Publication Date Title
CN105137185B (zh) 一种基于离散傅里叶变换的频域插值电力谐波分析方法
Evrenosoglu et al. Travelling wave based fault location for teed circuits
CN112379216A (zh) 输电线故障测距基波分量提取中量测误差抑制方法及系统
CN103576002B (zh) 一种容性绝缘设备介质损耗角的计算方法
Wen et al. FFT-based amplitude estimation of power distribution systems signal distorted by harmonics and noise
CN111222088B (zh) 一种改进的平顶自卷积窗加权电力谐波幅值估计方法
CN107247193B (zh) 天线近场测试方法及装置
CN110837003A (zh) 一种基于三角窗的双窗全相位dft同步相量测量方法及系统
CN105785124A (zh) 一种采用谱估计和互相关的电力系统谐波和间谐波测量方法
CN111984920A (zh) 次/超同步谐波参数识别方法、装置、设备和介质
CN109521273B (zh) 一种同步相量测量方法、系统及装置
CN109145476A (zh) 用于电力系统信号处理的时域自适应分段复指数级数法
CN115389877A (zh) 电缆绝缘故障的定位方法、装置、终端及存储介质
CN109270406B (zh) 一种基于双端不同步电气量计算过渡电阻的方法
CN109541304B (zh) 基于六项最小旁瓣窗插值的电网高次弱幅值谐波检测方法
CN103543331A (zh) 一种计算电信号谐波和间谐波的方法
CN112946515A (zh) 一种用于换流变网侧套管的在线监测方法和装置
CN112730982A (zh) 一种混合直流输电系统的谐波检测方法
CN113162002A (zh) 一种计及宽频测量环节的直流行波保护方法及系统
Parikh et al. Decision tree based fault classification scheme for protection of series compensated transmission lines
CN109342814B (zh) 用于确定电缆局放工频相位频率的方法、装置及存储介质
Tan et al. Harmonic analysis based on time domain mutual-multiplication window
CN114184838A (zh) 基于sn互卷积窗的电力系统谐波检测方法、系统及介质
Lessa et al. A travelling wave-based fault locator for radial distribution systems using decision trees to mitigate multiple estimations
CN113419133B (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20210219

RJ01 Rejection of invention patent application after publication