CN105763170B - 一种电力信号数字滤波方法 - Google Patents

一种电力信号数字滤波方法 Download PDF

Info

Publication number
CN105763170B
CN105763170B CN201610094330.0A CN201610094330A CN105763170B CN 105763170 B CN105763170 B CN 105763170B CN 201610094330 A CN201610094330 A CN 201610094330A CN 105763170 B CN105763170 B CN 105763170B
Authority
CN
China
Prior art keywords
value
sampling
filtering
sampled value
fluctuation
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
CN201610094330.0A
Other languages
English (en)
Other versions
CN105763170A (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.)
Zhejiang Harnpu Power Technology Co Ltd
Original Assignee
Zhejiang Harnpu Power Technology Co Ltd
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 Zhejiang Harnpu Power Technology Co Ltd filed Critical Zhejiang Harnpu Power Technology Co Ltd
Priority to CN201610094330.0A priority Critical patent/CN105763170B/zh
Publication of CN105763170A publication Critical patent/CN105763170A/zh
Application granted granted Critical
Publication of CN105763170B publication Critical patent/CN105763170B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00—Networks using digital techniques
    • H03H17/02—Frequency selective networks
    • H03H17/0201—Wave digital filters

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Physics (AREA)
  • Measurement Of Current Or Voltage (AREA)

Abstract

本发明公开一种电力信号数字滤波方法，该方法包括以下步骤：步骤1：获取保护级采样值WS1、采样波动限定值Wk、采样额定值W、采样波动连续超限次数上限J1、采样连续超额次数上限J2、和测量级采样值WS2；步骤2：对保护级采样值WS1进行分析：结合步骤1给定的参数，对检测到采样波动和采样超额进行计数来滤除偶然干扰；根据分析结果选择保护级采样值WS1或者测量级采样值WS2作为本次滤波采样值；步骤3：将步骤2获得的滤波采样值进行累加、求平均计算，得出累加平均值，作为测量输出值。

Description

一种电力信号数字滤波方法
技术领域
本发明涉及电力***中的电力信号的处理技术，具体涉及一种具备测量级准确度与保护级动态范围的电力信号数字滤波方法。
背景技术
随着现代化工业的不断发展，电力需求越来越大，电力输送的稳定性越来越重要，人们对电力信号监控的测量精度、响应时间、干扰滤波提出了更高要求。
本发明中提及的数字滤波主要针对电力***中信号采样值在输出前进行的处理。通常而言，在电力信号的变送输出中，***稳定、测量精度、响应时间和干扰滤波并重是目前解决的核心和难点。电力信号的一般数字滤波中，电信号被采样电路转化成采样值后，将连续N次转化所得采样值组成一个采样值队列，对采样值队列进行累加平均计算，用这个计算值相对上次输出值的波动被均分为平均了波动幅度，但对测量输出的动态响应和准确性并无作用，抗干扰能力也有待提高。
在分布式发电***中，有需要测量准确性的监控要求，有需要很好的动态响应保证实时性的场合，也有同时兼顾准确度和动态响应以适应高要求的环境，这样就需要一种即保证测量级准确度又保证保护级动态响应的测量装置。因此，本文提出了一种具备测量级准确度与保护级动态范围的电力信号数字滤波方法。
发明内容
有鉴于此，本发明的目的在于提供一种具备测量级准确度与保护级动态范围的电力信号数字滤波方法，其利用对检测到采样波动和采样超额进行计数来滤除偶然干扰。另外将相应的参数都作为可配置量，来适应不同的应用中对环境的适应性。
为达到以上目的，本发明是通过以下技术方案实现的：
一种电力信号数字滤波方法，该方法包括以下步骤：
步骤1：获取保护级采样值WS1、采样波动限定值Wk、采样额定值W、采样波动连续超限次数上限J1、采样连续超额次数上限J2、和测量级采样值WS2；
步骤2：对保护级采样值WS1进行分析，结合步骤1给定的参数，对检测到采样波动和采样超额进行计数来滤除偶然干扰；根据分析结果选择保护级采样值WS1或者测量级采样值WS2作为本次滤波采样值；也就是说，按分析结果选择保护级采样值WS1或者测量级采样值WS2作为有效的滤波采样值，并将有效的滤波采样值替换采样队列的最旧值；
步骤3：将步骤2获得的滤波采样值进行累加、求平均计算，得出累加平均值，作为测量输出值，将测量输出值记为WPn。然后重复步骤1-3，获得下一次的测量输出值。
与现有技术相比，本发明的有益效果在于，本发明通过各项可配置参数的分析，选择是测量级采样值参数或者保护级采样值参数进入采样队列，再进行算术平均滤波并输出。
进一步的，步骤1具体包括，通过保护级传感器和测量级传感器获取电力信号，并将该信号调制为一定范围的弱电信号，经过换算得到本次采样的保护级采样值WS1和测量级采样值WS2；同时设置好各项参考量：采样波动限定值Wk、采样额定值W、采样波动连续超限次数上限J1、采样连续超额次数上限J2，获取上一次测量输出值WPn。
进一步的，步骤2中，本步骤根据步骤1给定的参数分析，判断出我们需要的是保护级采样值还是测量级采样值作为本次滤波采样值，或者丢弃采样值。具体的，步骤2包括如下过程：
步骤(2.A)：如果|W1-WPn|≥Wk成立，则计算|W1-WPn|超出采样波动限定值Wk的次数M1；如果M1≥J1，即本次WS1波动大于采样波动限定值Wk次数足够多，排除偶然的脉冲干扰，由于一般保护级测量的动态范围更大，取WS1作为本次滤波采样值Wn+1；如果M1<J1，即本次WS1波动大于采样波动限定值Wk次数还少，可能是偶然的脉冲干扰，丢弃当前电力实时采样值WS1、WS2；
步骤(2.B)：如果|W1-WPn|≥Wk不成立，继而判断W1≥W是否成立，如果成立则计算WS1超出采样额定值W的次数M2；如果M2≥J2，本次WS1大于采样额定值W次数足够多，排除偶然的脉冲干扰，由于一般保护级测量的动态范围更大，取WS1为本次滤波采样值Wn+1；如果M2<J2，本次WS1大于采样额定值W次数很少，可能是偶然的脉冲干扰，由于测量级测量的准确度好，动态范围小，取WS2作为本次滤波采样值Wn+1；
步骤(2.C)：如果|W1-WPn|≥Wk不成立，继而当W1≥W也不成立时，说明此时测量采样值处于稳态，并在采样额定值W之内，此时取测量准确度高的WS2作为本次滤波采样值Wn+1。
进一步的，步骤3是进行滑动累加平均计算。其具体包括如下过程：将步骤2获得的本次滤波采样值记为Wn+1，将Wn+1加入长度为N的采样值信号队列[W1,W2,W3…Wn-1,Wn]，信号队列变成[Wn+1,W2,W3…Wn-1,Wn]，然后计算出累加平均值Wpn+1作为输出信号。N可配置，其大小决定滑动累加平均的速度。N＝1的时候，采样值直接作为测量输出值，响应实时快速。
本发明的有益效果在于：本发明所述的方法，利用保护级采样值WS1，采样波动限定值Wk，采样额定值W，采样波动连续超限次数上限J1，采样连续超额次数上限J2，测量级采样值WS2和采样值信号队列长度N，调整测量输出信号的平滑度、响应速度和准确度，改进了普通的数字滤波方法单一的滤波方式，进一步提高了滤波的灵活性。这对于电力信号各种不同环境的应用具有重要意义。
附图说明
图1为本发明的处理流程示意图；
图2是本发明的信号流图。
具体实施方式
现结合附图和具体实施方式对本发明进一步说明。
作为一个具体的实施例，图1为本发明的处理流程示意图，图2为信号流图，首先获取保护级采样值WS1，采样波动限定值Wk，采样额定值W，采样波动连续超限次数上限J1，采样连续超额次数上限J2，测量级采样值WS2；然后对WS1参考各设置参数分析，选择合适的采样值(保护级采样值WS1或者测量级采样值WS2)作为本次滤波采样值；然后将获得的多次滤波采样值进行累加、求平均计算，得出累加平均值，作为测量输出值。具体的，参见图1和图2，该电力信号的数字滤波方法包括四个步骤：
步骤2、电力信号采样信号有效值前64次的数据存入长度为64的RAM中，此队列标记为[W1,W2…W63,W64]，其累加平均值是上一次的输出信号WPn，即 计算本次保护级采样值和上次测量,输出值的差值|W1-WPn|，用这个差值的绝对值|W1-WPn|代表本次采样的波动值，并对此波动值进行判断，是否超过采样波动限定值Wk。判断公式是同时，根据需要判断本次滤波采样值是否超过采样额定值W，即W1≥W。同时对于和W1≥W的判断成功计数M1和M2，对M1和M2各自判断波动超限次数上限M1≥J1和额定超限次数上限M2≥J2，根据这四个判断标准有以下的情形：
步骤3、上一步是丢弃采样值WS1和WS2时，将上一次的输出信号WPn作为本次的输出信号WPn+1，这样偶然的脉冲干扰被滤除，电力信号变送输出稳定，记WPn+1＝WPn。在采用WS1或WS2作为滤波采样值时，本次采样值为有效数据，将本次滤波采样值Wn+1替换RAM中长度为64的采样信号队列[W1,W2…W65,W64]中的最旧值WS1，第一次替换，即用W65替换W1，最旧值的位置随着采样次数K变化为变化；将这个采样值队[W2…W63,W64,W65]做累加平均计算，得出本次输出信号 这一步中，采样信号队列长度64，采样波动限定值Wk，采样额定值W，采样波动连续超限次数上限J1，采样连续超额次数上限J2，存在ROM中，一般都在E2PROM中。配置的参数大小可决定真实的数据经过多少判断再反应到测量输出值。这一步的求和、除法运算，耗费一定的数字芯片时钟周期，根据不同的求和次数，时钟周期的消耗时间不能确定。按10000个时钟周期计算，在100M的芯片速率下，仅有0.1ms，延时时间在采样时间以内，忽略不计。在实际的应用环境中，输出的情况具体为以下三种：
最终输出信号比平滑较稳定。
2)丢弃偶然脉冲即最终输出WP(n+1)＝WPn。
3)超额功率电力信号上升较大，即最终输出经过几次循环，WPn+1接近稳定，并达到的稳态条件。在以上过程中，有M1从0增长到J1的过程，有M2从0增长到J2的过程，有更新长度为N的采样值队列的过程，此过程丢弃部分采样值，不能及时响应有效的采样值，造成响应延时。因为响应实时性要求高时，首先将J1、J2、N设置的更小。
步骤4、本次的输出信号WPn+1，常见的一种是，在数字芯片中作为PWM波输出，再经过数模转换电路，变送输出为4-20mA的电路。当信号WPn+1为0，对应输出4mA电流；当信号WPn+1到一个规定的满度值，对应输出20mA电流。
本发明说明书中未做详细描述的内容属于本领域专业技术人员的公知知识。
尽管结合优选实施方案具体展示和介绍了本发明，但所属领域的技术人员应该明白，在不脱离所附权利要求书所限定的本发明的精神和范围内，在形式上和细节上可以对本发明做出各种变化，均为本发明的保护范围。

Claims (4)

1.一种电力信号数字滤波方法，该方法包括以下步骤：
步骤1：获取保护级采样值WS1、采样波动限定值Wk、采样额定值W、采样波动连续超限次数上限J1、采样连续超额次数上限J2、和测量级采样值WS2；其中，保护级采样值WS1定义如下：
其中，Pt为t时刻的实时功率，P1为t＝1时的实时功率，P257为t＝257时的实时功率；
步骤2：对保护级采样值WS1进行分析：结合步骤1给定的参数，对检测到采样波动和采样超额进行计数来滤除偶然干扰；根据分析结果选择保护级采样值WS1或者测量级采样值WS2作为本次滤波采样值；
步骤3：将步骤2获得的滤波采样值进行累加、求平均计算，得出累加平均值，作为测量输出值Wpn；
所述步骤2具体包括如下过程：
如果|W1-Wpn|≥Wk成立，则计算|W1-Wpn|超出采样波动限定值Wk的次数M1；如果M1≥J1，则取WS1作为本次滤波采样值Wn+1；如果M1＜J1，则丢弃当前电力实时采样值WS1、WS2；其中，W1是第一次采样的信号有效值；Wpn是上一次的滤波采样值的累加平均值；
如果|W1-Wpn|≥Wk不成立，同时W1≥W成立，则计算WS1超出采样额定值W的次数M2；如果M2≥J2，则取WS1为本次滤波采样值Wn+1；如果M2＜J2，则取WS2作为本次滤波采样值Wn+1；
如果|W1-Wpn|≥Wk不成立，同时当W1≥W也不成立时，则取测量准确度高的WS2作为本次滤波采样值Wn+1。
2.根据权利要求1所述的一种电力信号数字滤波方法，其特征在于，所述步骤1具体包括，通过保护级传感器和测量级传感器获取电力信号，并将该信号调制为一定范围的弱电信号，经过换算得到本次采样的保护级采样值WS1和测量级采样值WS2；获取上一次测量输出值Wpn；同时设置好各项参考量：采样波动限定值Wk、采样额定值W、采样波动连续超限次数上限J1和采样连续超额次数上限J2。
3.根据权利要求1或2所述的一种电力信号数字滤波方法，其特征在于：所述步骤3是进行滑动累加平均计算。
4.根据权利要求3所述的一种电力信号数字滤波方法，其特征在于：步骤3具体包括如下过程：将步骤2获得的本次滤波采样值记为Wn+1，将Wn+1加入长度为N的采样值信号队列[W1,W2,W3…Wn-1,Wn]，信号队列变成[Wn+1,W2,W3…Wn-1,Wn]，然后计算出累加平均值Wpn+1作为输出信号。
CN201610094330.0A 2016-03-28 2016-03-28 一种电力信号数字滤波方法 Active CN105763170B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610094330.0A CN105763170B (zh) 2016-03-28 2016-03-28 一种电力信号数字滤波方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610094330.0A CN105763170B (zh) 2016-03-28 2016-03-28 一种电力信号数字滤波方法

Publications (2)

Publication Number Publication Date
CN105763170A CN105763170A (zh) 2016-07-13
CN105763170B true CN105763170B (zh) 2018-09-18

Family

ID=56330986

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610094330.0A Active CN105763170B (zh) 2016-03-28 2016-03-28 一种电力信号数字滤波方法

Country Status (1)

Country Link
CN (1) CN105763170B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106918735B (zh) * 2017-03-28 2019-07-02 辽宁大学 一种电流测量数据处理方法
CN106992661A (zh) * 2017-04-17 2017-07-28 广西大学 用于pwm数字控制的过采样即时信号处理方法
CN110221544A (zh) * 2019-06-14 2019-09-10 四川烟叶复烤有限责任公司 一种烟叶复烤入口水分数据滤波方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101895274A (zh) * 2010-07-21 2010-11-24 珠海天威技术开发有限公司 数字滤波电路及其滤波方法、耗材芯片
CN102118145A (zh) * 2010-12-16 2011-07-06 中国兵器工业第二○六研究所 改进的限速滤波方法
CN102684465A (zh) * 2012-04-28 2012-09-19 北方工业大学 基于电流指令生成的电网不平衡时pwm整流器控制方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10106403B4 (de) * 2001-02-12 2007-01-18 Rohde & Schwarz Gmbh & Co. Kg Vorrichtung und Verfahren zur Abtastratenumsetzung
CN101351967B (zh) * 2005-12-28 2012-05-30 美国亚德诺半导体公司 针对数模转换器和低通滤波器组合连续时间级与开关电容级的架构

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101895274A (zh) * 2010-07-21 2010-11-24 珠海天威技术开发有限公司 数字滤波电路及其滤波方法、耗材芯片
CN102118145A (zh) * 2010-12-16 2011-07-06 中国兵器工业第二○六研究所 改进的限速滤波方法
CN102684465A (zh) * 2012-04-28 2012-09-19 北方工业大学 基于电流指令生成的电网不平衡时pwm整流器控制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"基于STM32与AD7606的高精度和快速响应数字多功能表的设计";徐国明等;《电测与仪表》;20150625;第52卷(第12期);102-106 *

Also Published As

Publication number Publication date
CN105763170A (zh) 2016-07-13

Similar Documents

Publication Publication Date Title
CN116026389A (zh) 一种基于数据分析的智能传感器运行检测***
CN105763170B (zh) 一种电力信号数字滤波方法
CN105743462A (zh) 测量响应时间可调的电力信号数字滤波方法
CN109284933B (zh) 一种基于数理统计的电子式互感器状态评估***及方法
CN115330000A (zh) 一种工业自动化控制仪表运行智能监测管理***
CN115687983A (zh) 一种桥梁健康状态监测方法、***及电子设备
CN111812703B (zh) 一种伽玛能谱背景计数方法
CN111998918A (zh) 一种误差校正方法、误差校正装置及流量传感***
CN101988957B (zh) 变电站合并单元电量采样准确性的动态检测方法及***
CN114608678A (zh) 一种基于脉冲法的水表校准方法及装置
CN117293758B (zh) 基于故障识别的数字化配电馈线监测终端自动化保护方法
CN110488188A (zh) 基于动态阈值的机组三维健康量化评价方法
CN111707867B (zh) 一种交流电流测量方法、装置、电子设备及存储介质
CN108801320B (zh) 一种天然气测量***的诊断方法
CN105865606A (zh) 一种汽车衡故障诊断装置及诊断方法
CN103605323B (zh) 化工生产的离散控制方法及装置
CN109193639B (zh) 一种电力***抗差估计方法
CN112198208A (zh) 一种甲醛检测方法及甲醛检测装置
CN115950495A (zh) 一种超声水表测量频率调节方法
CN105115998A (zh) 一种粉末衍射择优取向参数确定方法
CN109974873B (zh) 一种移动平均的温度监测方法
CN112147409A (zh) 低压配电线路在线阻抗检测方法及装置
CN108982908B (zh) 一种导线顺线向积分尺度测量方法和***
CN117895660B (zh) 一种电力终端能耗数据采集分析处理方法
CN110849796A (zh) 腐蚀监测***

Legal Events

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