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

[0001]
本发明涉及输电线故障测距技术领域，尤其涉及一种输电线故障测距基波分量提取中量测误差抑制方法及系统。


背景技术：

[0002]
本部分的陈述仅仅是提供了与本发明相关的背景技术信息，不必然构成在先技术。
[0003]
现有基于量测阻抗的故障测距算法都需要利用故障后电压和电流的基波分量进行计算。基波分量的准确性会极大的影响故障测距精度。而量测误差的存在，会使得用于测距计算的电压电流基波分量存在误差，影响测距结果的精度。
[0004]
为克服量测误差对测距精度的影响，现有方法一般是通过滑动采样数据窗提取多组基波分量构建冗余测距方程，利用最小二乘法提升测距精度。虽然该方法取得了一定效果，但对于配电网中较为常见的t接线或多分支接线而言，求解的测距方程比较复杂，利用冗余数据增加方程组数较多，不利于非线性测距方程的数值求解。


技术实现要素：

[0005]
为了解决上述问题，本发明提出了一种输电线故障测距基波分量提取中量测误差抑制方法及系统，可以直接对采样数据进行滤波处理，得到更加精准的基波分量计算结果，从而提高故障测距结果的精度。
[0006]
在一些实施方式中，采用如下技术方案：
[0007]
一种输电线故障测距基波分量提取中量测误差抑制方法，包括：
[0008]
通过滑动采样数据窗口的起始时刻形成多组采样数据；
[0009]
使用离散傅里叶变换从每组采样数据中分别提取基波分量；
[0010]
对基波分量的幅值和初相分别进行线性拟合，得到准确的幅值和相位。
[0011]
在另一些实施方式中，采用如下技术方案：
[0012]
一种输电线故障测距基波分量提取中量测误差抑制系统，包括：
[0013]
数据采样模块，用于通过滑动采样数据窗口的起始时刻形成多组采样数据；
[0014]
基波分量提取模块，用于使用离散傅里叶变换从每组采样数据中分别提取基波分量；
[0015]
线性拟合模块，用于对基波分量的幅值和初相分别进行线性拟合，得到准确的幅值和相位。
[0016]
在另一些实施方式中，采用如下技术方案：
[0017]
一种终端设备，其包括处理器和计算机可读存储介质，处理器用于实现各指令；计算机可读存储介质用于存储多条指令，所述指令适于由处理器加载并执行上述的输电线故障测距基波分量提取中量测误差抑制方法。
[0018]
在另一些实施方式中，采用如下技术方案：
[0019]
一种计算机可读存储介质，其中存储有多条指令，所述指令适于由终端设备的处理器加载并执行上述的输电线故障测距基波分量提取中量测误差抑制方法。
[0020]
与现有技术相比，本发明的有益效果是：
[0021]
本发明根据故障后电压和电流基波分量幅值和初相位恒定不变特征，提出一种新的量测误差抑制方法，可以直接对采样数据进行滤波处理，得到更加精准的基波分量计算结果，从而提高故障测距结果的精度。
[0022]
本发明的其他特征和附加方面的优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本方面的实践了解到。
附图说明
[0023]
图1为本发明实施例中输电线故障测距基波分量提取中量测误差抑制方法流程图。
具体实施方式
[0024]
应该指出，以下详细说明都是例示性的，旨在对本申请提供进一步的说明。除非另有指明，本发明使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
[0025]
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
[0026]
实施例一
[0027]
在一个或多个实施方式中，公开了一种输电线故障测距基波分量提取中量测误差抑制方法，参照图1，具体包括如下步骤：
[0028]
(1)通过滑动采样数据窗口的起始时刻形成多组采样数据；
[0029]
(2)使用离散傅里叶变换从每组采样数据中分别提取基波分量；
[0030]
(3)对基波分量的幅值和初相分别进行线性拟合，得到准确的幅值和相位。
[0031]
具体地，假设故障后基波分量表达式为：
[0032]
u1(t)＝u1cos(w1t+θ0)
ꢀꢀꢀ
(1)
[0033]
式中，w1＝2πf1为基波角频率，f1为基波分量的频率，θ0为基波分量的初相位，简称初相，u1为基波分量的幅值。根据基波分量的表达式可知其幅值和初相不随故障后时间推移变化，是常量。相位则是随时间变化的线性函数。
[0034]
因此，利用故障后工频一周波采样点进行基波分量提取，如果结果准确，理论上其幅值不随采样点数据窗起点时刻变化而变化。
[0035]
假设故障起始时刻为0时刻，采样周期为δt，如按δt步长依次平移提取基波分量的一周波采样点数据窗，对于起始时刻为kδt的数据窗，则相位是步长δt的线性函数起始时刻的相角
[0036]
根据上述基波分量的幅值和初相恒定不变的特征规律，本实施例选用了线性拟合
滤波，通过滑动数据窗口起始时刻形成多组采样数据，使用离散傅里叶变换(dft)从每个数据窗口分别提取基波分量，并求出幅值和初相，这些基波分量的幅值和初相围绕对应的准确值上下波动，进而可以借助线性拟合获得准确的幅值和相位。
[0037]
下面以故障后量测点m处所得ab相间线电压基波分量的幅值u
mab
和初相θ0为例进行具体阐述。
[0038]
在故障后选取32个连续的采样数据窗口，每个数据窗的长度是基波分量的一个周期，使用离散傅里叶变换(dft)从每个采样数据窗口内提取基波分量并将幅值记为u
mab
(i)(i＝0,1,2,
……
,31)，待求的线性拟合结果为u
mab
，基于最小二乘法的线性拟合是指拟合一个近似直线y＝u
mab
,使得对于所有数据来说总的误差和err1最小，即求：
[0039][0040]
为使误差和最小，对误差函数求导，并且令导函数等于0：
[0041][0042]
进而有：
[0043][0044][0045]
再对初相进行拟合，由于在数据窗内相角随时间一次线性变化，所以相角满足关系式基波分量频率f1＝50hz，w1＝2πf1＝100πrad/s,则由于初相θ0与时间无关，故是个恒定不变的常数。
[0046]
将32个数据窗中提取的基波分量的相角记为为α
k
(k＝0,1,2
……
,31)，则计算出的初相θ
k
＝α
k-w1kδt(k＝0,1,2
……
,31),待求的线性拟合结果为θ0，基于最小二乘法的线性拟合是指拟合一个近似直线y＝θ0，使得对于所有θ
k
，总的误差和err2最小，即求：
[0047][0048]
为使误差和最小，误差函数对θ0求导，并且令导函数等于0：
[0049][0050]
进而解出：
[0051][0052]
为了验证所提出的量测误差抑制方法的自洽性和有效性，在matlab中对进行仿真验证，首先构造信号
[0053][0054]
式(9)中基波分量的幅值为100，初相位为再给信号添加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％。数据表明本发明所提基于线性拟合的量测噪声抑制方法能够很好地抑制随机噪音的干扰，基波分量幅值和初相计算结果相比直接提取结果更接近基波分量的准确值。
[0055]
表1 20组含随机噪声的线性拟合结果和未经线性拟合基本分量提取结果比较
[0056][0057]
[0058]
实施例三
[0059]
在一个或多个实施方式中，公开了一种主端可穿戴遥操作控制装置的遥操作控制器，包括服务器，所述服务器包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现实施例一中的方法。为了简洁，在此不再赘述。
[0060]
应理解，本实施例中，处理器可以是中央处理单元cpu，处理器还可以是其他通用处理器、数字信号处理器dsp、专用集成电路asic，现成可编程门阵列fpga或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
[0061]
存储器可以包括只读存储器和随机存取存储器，并向处理器提供指令和数据、存储器的一部分还可以包括非易失性随机存储器。例如，存储器还可以存储设备类型的信息。
[0062]
在实现过程中，上述方法的各步骤可以通过处理器中的硬件的集成逻辑电路或者软件形式的指令完成。
[0063]
实施例一中的方法可以直接体现为硬件处理器执行完成，或者用处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器、闪存、只读存储器、可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器，处理器读取存储器中的信息，结合其硬件完成上述方法的步骤。为避免重复，这里不再详细描述。
[0064]
实施例四
[0065]
在一个或多个实施方式中，公开了一种计算机可读存储介质，其中存储有多条指令，所述指令适于由终端设备的处理器加载并实施例一中所述的一种输电线故障测距基波分量提取中量测误差抑制方法。
[0066]
上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。