CN105763170A - 一种电力信号数字滤波方法 - Google Patents
一种电力信号数字滤波方法 Download PDFInfo
- Publication number
- CN105763170A CN105763170A CN201610094330.0A CN201610094330A CN105763170A CN 105763170 A CN105763170 A CN 105763170A CN 201610094330 A CN201610094330 A CN 201610094330A CN 105763170 A CN105763170 A CN 105763170A
- Authority
- CN
- China
- Prior art keywords
- sampling
- value
- filtering
- fluctuation
- power signal
- 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
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,该电力信号的数字滤波方法包括四个步骤:
步骤1、电力信号中,最基本的就是电流电压功率信号。通过传感器获取实时电压电流信号Vt和It,模型Vt=A1*sin(ω1t+δ1),It=A2*sin(ω2t+δ2),将该信号在时间t经过采样电路的电压调制为一定范围的弱电信号Vt和It,此信号输入到芯片采样脚,经过计算得出的实时功率Pt为一个采样点。从而在数字芯片中利用一个周期所有采样点的均方根得到本次功率采样信号Wt。在此过程中,一个周波的采样分为256个点,256个点的均方根值作为本次采样的有效值。当Wst经过一个周期256个点的采集[P1,P2,P3…P255,P256],经过均方根计算得出被采集功率的有效值Wn。下一个采样点替换256点采样点中最旧值,新的采集点[P2,P3…P256,P257]采样信号有效值Wn+1。这一步延时约为采集时间为f在市电频率频率50hz时延时t为0.078125ms。这个采样值就作为保护级的采样值WS1, 测量级的采样值WS2一般有高精度的电能计量专用集成芯片获得。根据实际应用设置好的保护级采样值WS1,采样波动限定值Wk,采样额定值W,采样波动连续超限次数上限J1,采样连续超额次数上限J2,测量级采样值WS2,采样值队列长度N。在本次采样前已有上次输出测量值WPn。
步骤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,根据这四个判断标准有以下的情形:
判断公式成立,计算超出采样波动限定值Wk的次数M1,在M1≥J1的情况下,即本次WS1波动大于采样波动限定值Wk次数足够多,排除偶然的脉冲干扰,由于一般保护级测量的动态范围更大,采用测量采样值WS1作为滤波采样值,组合判断公式判断公式成立,计算WS1超出采样波动限定值Wk的次数M1,在M1<J1的情况下,即本次WS1波动大于采样波动限定值Wk次数少,可能是偶然的脉冲干扰,丢弃当前电力实时采样值WS1、WS2,直接按原测量值输出,组合判断公式为判断不成立,继而判断W1≥W成立,则计算WS1超出采样额定值W的次数M2,在M2≥J2的情况下,本次WS1大于采样额定值W次数足够多,排除偶然的脉冲干扰,取动态范围大的保护级WS1值,组合判断公式 判断不成立,继而判断W1≥W成立,则计算WS1超出采样额定值W的次数M2,在M2<J2的情况下,本次WS1大于采样额定值W次数很少,可能是偶然的脉冲干扰,由于测量级测量的准确度好,动态范围小,取WS2作为本次滤波采样值,组合判断公式为 判断不成立,继而判断W1≥W也不成立,这就是保护级采样值判断到波动较小,也没有超过采样额定值,则取测量级采样值WS2为滤波采样值,可保证滤波采样值取得的精确度和稳定度,组合判断公式为
步骤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,延时时间在采样时间以内,忽略不计。在实际的应用环境中,输出的情况具体为以下三种:
1)稳态即
最终输出信号比平滑较稳定。
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 (5)
1.一种电力信号数字滤波方法,该方法包括以下步骤:
步骤1:获取保护级采样值WS1、采样波动限定值Wk、采样额定值W、采样波动连续超限次数上限J1、采样连续超额次数上限J2、和测量级采样值WS2;
步骤2:对保护级采样值WS1进行分析:结合步骤1给定的参数,对检测到采样波动和采样超额进行计数来滤除偶然干扰;根据分析结果选择保护级采样值WS1或者测量级采样值WS2作为本次滤波采样值;
步骤3:将步骤2获得的滤波采样值进行累加、求平均计算,得出累加平均值,作为测量输出值WPn。
2.根据权利要求1所述的一种电力信号数字滤波方法,其特征在于,所述步骤2具体包括如下过程:
如果|W1-WPn|≥Wk成立,则计算|W1-WPn|超出采样波动限定值Wk的次数M1;
如果M1≥J1,则取WS1作为本次滤波采样值Wn+1;如果M1<J1,则丢弃当前电力实时采样值WS1、WS2;
如果|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。
3.根据权利要求1或2所述的一种电力信号数字滤波方法,其特征在于,所述步骤1具体包括,通过保护级传感器和测量级传感器获取电力信号,并将该信号调制为一定范围的弱电信号,经过换算得到本次采样的保护级采样值WS1和测量级采样值WS2;获取上一次测量输出值WPn;同时设置好各项参考量:采样波动限定值Wk、采样额定值W、采样波动连续超限次数上限J1和采样连续超额次数上限J2。
4.根据权利要求1或2所述的一种电力信号数字滤波方法,其特征在于:所述步骤3是进行滑动累加平均计算。
5.根据权利要求4所述的一种电力信号数字滤波方法,其特征在于:步骤3具体包括如下过程:将步骤2获得的本次滤波采样值记为Wn+1,将Wn+1加入长度为N的采样值信号队列[W1,W2,W3…Wn-1,Wn],信号队列变成[Wn+1,W2,W3…Wn-1,Wn],然后计算出累加平均值WPn+1作为输出信号。
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 true CN105763170A (zh) | 2016-07-13 |
CN105763170B 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) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106918735A (zh) * | 2017-03-28 | 2017-07-04 | 辽宁大学 | 一种电流测量数据处理方法 |
CN106992661A (zh) * | 2017-04-17 | 2017-07-28 | 广西大学 | 用于pwm数字控制的过采样即时信号处理方法 |
CN110221544A (zh) * | 2019-06-14 | 2019-09-10 | 四川烟叶复烤有限责任公司 | 一种烟叶复烤入口水分数据滤波方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020109617A1 (en) * | 2001-02-12 | 2002-08-15 | Markus Freidhof | Device and method for sampling rate conversion |
WO2007079097A2 (en) * | 2005-12-28 | 2007-07-12 | Analog Devices, Inc. | Architecture combining a continuous-time stage with a switched-capacitor stage for digital-to-analog converters and low-pass filters |
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整流器控制方法 |
-
2016
- 2016-03-28 CN CN201610094330.0A patent/CN105763170B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020109617A1 (en) * | 2001-02-12 | 2002-08-15 | Markus Freidhof | Device and method for sampling rate conversion |
WO2007079097A2 (en) * | 2005-12-28 | 2007-07-12 | Analog Devices, Inc. | Architecture combining a continuous-time stage with a switched-capacitor stage for digital-to-analog converters and low-pass filters |
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)
Title |
---|
徐国明等: ""基于STM32与AD7606的高精度和快速响应数字多功能表的设计"", 《电测与仪表》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106918735A (zh) * | 2017-03-28 | 2017-07-04 | 辽宁大学 | 一种电流测量数据处理方法 |
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 | 四川烟叶复烤有限责任公司 | 一种烟叶复烤入口水分数据滤波方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105763170B (zh) | 2018-09-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105743462A (zh) | 测量响应时间可调的电力信号数字滤波方法 | |
CN104646565B (zh) | 一种全自动冲压生产线实时监控系统 | |
CN107886171B (zh) | 一种基于pmu数据的断路器状态在线诊断方法及系统 | |
CN103034170B (zh) | 一种基于区间的数控机床性能预测方法 | |
CA2764082A1 (en) | Method and monitoring apparatus for automated surveillance of a wind turbine and a method for creating a linear model | |
CN105763170A (zh) | 一种电力信号数字滤波方法 | |
CN103336255B (zh) | 一种高精度故障指示器的校准方法 | |
EP1085534A3 (en) | Intelligent analysis system and method for fluid-filled electrical equipment | |
CN109284933B (zh) | 一种基于数理统计的电子式互感器状态评估系统及方法 | |
CN114065875A (zh) | 基于大数据的电力电网故障识别系统 | |
CN111382789A (zh) | 基于机器学习的电力负荷识别方法及系统 | |
EP1085635A3 (en) | Fluid-filled electrical equipment intelligent analysis system and method | |
CN203942803U (zh) | 用于检测谷物产量的装置 | |
CN115979350A (zh) | 一种海洋监测设备数据采集系统 | |
CN202900232U (zh) | 油水井工况数据远传计量管理系统 | |
CN207487708U (zh) | 一种磁性编码器校准装置 | |
CN112650068B (zh) | 一种基于智能电流传感器的家居故障检测方法及装置 | |
CN117957396A (zh) | 用于公共服务分配系统中的声学泄漏检测的系统和方法 | |
CN105526888A (zh) | 测厚方法及测厚系统 | |
CN104569675A (zh) | 一种电能计量芯片中的闪变检测电路及检测方法 | |
CN115455689A (zh) | 一种充电桩电能计量误差的高精度分析方法 | |
CN208270142U (zh) | 一种压力/差压变送器响应时间测试系统 | |
CN115219000A (zh) | 一种药箱水位自动标定方法及装置 | |
CN105022012A (zh) | 数字式电能质量监测终端的闪变精度检测方法及系统 | |
CN110729732B (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 |