CN104320157A - 一种电力线双向工频通信上行信号检测方法 - Google Patents

一种电力线双向工频通信上行信号检测方法 Download PDF

Info

Publication number
CN104320157A
CN104320157A CN201410597083.7A CN201410597083A CN104320157A CN 104320157 A CN104320157 A CN 104320157A CN 201410597083 A CN201410597083 A CN 201410597083A CN 104320157 A CN104320157 A CN 104320157A
Authority
CN
China
Prior art keywords
signal
imf
frequency
power
noise
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
Application number
CN201410597083.7A
Other languages
English (en)
Other versions
CN104320157B (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.)
State Grid Corp of China SGCC
Northeast Electric Power University
Information and Telecommunication Branch of State Grid Jilin Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Northeast Dianli University
Information and Telecommunication Branch of State Grid Jilin Electric Power 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 State Grid Corp of China SGCC, Northeast Dianli University, Information and Telecommunication Branch of State Grid Jilin Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201410597083.7A priority Critical patent/CN104320157B/zh
Publication of CN104320157A publication Critical patent/CN104320157A/zh
Application granted granted Critical
Publication of CN104320157B publication Critical patent/CN104320157B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of Resistance Or Impedance (AREA)
  • Monitoring And Testing Of Transmission In General (AREA)

Abstract

本发明是一种电力线工频通信上行信号检测方法,其特点是:包括在双向工频通信上行信号检测时,首先通过加权求和法使调制信号的背景工频信号相互抵消并且增强调制信号,再对增强后的调制信号进行求导,利用集合经验模分解对调制信号进行经验模分解;通过相关运算找到n阶本征模函数中150HZ~500HZ所在的阶数,对所在阶数的本征模函数找到对应阶数的阈值并进行阈值去噪,对去噪处理后的信号利用希尔伯特变换计算信号的希尔伯特边际谱,从而检测调制信号所在时域;能有效的去除大部分噪声及谐波干扰,准确检测调制信号所在时域。具有科学合理,容易实现,检测精度高等优点。

Description

一种电力线双向工频通信上行信号检测方法
技术领域
本发明涉及一种电力线载波通信的信号检测方法,更具体地说,是一种电力线双向工频通信上行信号检测方法。
背景技术
双向工频通信系统是近年来新兴的一种基于电力配电网络的通信系统,该系统以电力线为数据传输载体,利用过零调制原理实现信号的调制,工频通信上行信号是通过晶闸管的开关控制在基波电压信号过零点附近处产生的畸变电流脉冲信号,是一种频谱随时间变化的时变信号。双向工频通信上行信号检测是在大背景下检测小信号,传统的检测方法是时域检测法,这种方法会损失大量有用的信息,对电网有一定要求,具有一定的局限性,现有的时域检测方法使用最普遍的小波变换信号检测算法,此方法对调制信号进行时频分析时会受到小波基选取的影响,导致无法进行精确的时频分析。
发明内容
本发明的目的是提供一种科学合理,简便实用,检测准确的电力线双向工频通信上行信号检测方法,其特征在于,它包括如下步骤:
1)对信号加权求和消除背景信号:
工频通信上行信号的检测是在大背景电流中提取微弱电流,所以利用加权求和抵消工频信号,,利用其中IA,IB,IC是配电网中的三相电流,能消除绝大部分的基波电流影响,同时使信号增强1倍;
2)对加权求和后的通信信号进行求导:
双向工频通信上行调制信号是经晶闸管的开关在工频电压过零点处产生的电流脉冲,每个脉冲含有一个极大值或一个极小值以及两个过零点,对其求导后,每个脉冲含有一个极大值,一个极小值和一个过零点,满足EEMD的分解条件,所述EEMD为集合经验模分解的英文缩写,本领域通常采用其英文缩写EEMD;
3)对处理后的信号进行EEMD:
在原始信号中加入均匀分布的高斯白噪声用于填补原始调制畸变电流中部分尺度缺失,使待测的调制电流在不同程度上连续;选定进行EMD的分解次数和白噪声的幅值系数,分解次数取100,幅值系数取0.2,所述EMD为经验模分解的英文缩写,本领域通常采用其英文缩写EMD;
4)测量各本征模函数频率:
通过傅里叶变换测量各阶IMF的频率,找出150HZ~500HZ频率所在的阶数进行保留,所述IMF为本征模函数的英文缩写,本领域通常采用其英文缩写IMF;
5)对保留的IMF进行阈值去噪:
由于电网中的噪声都是随机的,其能量大部分是均匀分布,有用信号的能量相对集中,根据有用信号和噪声信号能量的分布不同选择合适的阈值对保留的IMF进行降噪处理;
6)计算各阶IMF希尔伯特边际谱:
获得IMF分量后,对于其中一个IMF分量记为imf(t),经过希尔伯变换后的H[imf(t)]与原分量imf(t)组成解析信号,记为其中 瞬时角频率为H(ω,t)是振幅I(t)以灰度形式在频率-时间平面上显示的希尔伯特谱,对H(ω,t)在时间上积分,进而求得希尔伯特边际谱 h ( ω ) = ∫ 0 T H ( ω , t ) dt .
所述步骤3)对处理后的信号进行EEMD步骤如下:
(1)在原调制畸变电流信号中加入幅值系数为0.2的高斯白噪声,对所得信号进行EMD,得到n个IMF记为imfj
(2)对原调制畸变电流信号添加的下一个白噪声序列重复上述步骤,获得各自的IMF;
(3)计算剩余信号Ri,1=xi(t)-imfi,1(t),xi(t)为原始信号,imfi,1(t)为混入第i个白噪声序列分解得到的第一阶IMF,把剩余信号作为待处理信号进行EMD分解,依次获得imfi,1,imfi,2,....imfi,j
(4)对所得结果求平均其中j=1,2,…n,n为分解的阶数,去除添加的白噪声,得到待测信号。
本发明与现有技术相比,本发明的有益效果是:
在双向工频通信上行信号检测时,首先通过加权求和法使调制信号的背景工频信号相互抵消并且增强调制信号,然后对增强后的调制信号进行求导,求导目的就是使原始信号满足EEMD的条件,利用EEMD对调制信号进行EMD,大大减小了EEMD的计算量,再通过相关运算找到n阶IMF中150HZ~500HZ的阶数,对此类IMF进行阈值降噪处理,对降噪后的信号进行利用希尔伯特变换计算希尔伯特边际谱,从而检测调制信号所在时域;本发明不存在小波基选取的问题,能有效的去除大部分噪声及谐波干扰,准确检测调制信号所在时域。其科学合理,容易实现,检测精确度高。
附图说明:
图1为本发明的一种电力线双向工频通信上行信号检测方法的流程框图;
图2为双向工频通信上行信号调制等效电路原理图;
图3为双向工频通信上行信号调制信号调制波形图。
具体实施方式
下面结合附图对本发明的一种电力线双向工频通信上行信号检测方法作进一步具体说明。
本发明的一种电力线双向工频通信上行信号检测方法,包括以下步骤:
1.对信号加权求和消除背景信号
工频通信上行信号检测是在大背景中提取微弱电流,所以利用加权求和抵消工频信号,利用其中IA,IB,IC是配电网中的三相电流,能消除绝大部分的基波电流影响,同时使信号增强1倍;
2.对加权求和后的通信信号进行求导
如图2双向工频通信上行调制信号调制等效电路原理图所示,上行信号是通过对电流基波波形的调制来传送信息,信号由用户终端发出,携带相关的信息传回子站,其调制信号的频率在150Hz~500Hz范围内变化;
如图3上行调制信号的调制波形图所示,经晶闸管的开关在工频电压过零点处产生的电流脉冲,每个脉冲含有一个极大值或一个极小值以及两个过零点,对这样的信号进行EEMD分解计算量会非常大,对其求导后,每个脉冲含有一个极大值,一个极小值和一个过零点,满足EEMD的分解条件,大大减小了计算量;对其求导的计算过程为本领域技术人员所熟悉的简单数学运算;
3.对处理后的信号进行EEMD分解
在原信号中加入均匀分布的高斯白噪声,填补原始调制畸变电流中部分尺度缺失,使待测的调制电流在不同程度上连续;,减小模态混叠的程度;选定进行EMD的总次数n与白噪声的幅值系数k,这里n取100,k取0.2,分解步骤如下:
1)在原始信号中加入幅值系数为0.2的高斯白噪声,对所得待分解信号进行EMD,得到n个本征模函数记为imfj
2)对原始信号加入的下一个混入白噪声序列重复上述步骤,获得各自的IMF;
3)计算剩余信号Ri,1=xi(t)-imfi,1(t),xi(t)为原始信号,imfi,1(t)为混入第i个白噪声序列分解得到的第一阶IMF,把剩余信号作为待处理信号进行EMD分解,依次获得imfi,1,imfi,2,....imfi,j
4)对所得结果求平均其中j=1,2,...n,n为分解的阶数,去除添加的白噪声,得到待测信号。
4.测量各IMF的频率:
测量各阶IMF的频率,上行信号是通过对电流基波波形的调制来传送信息,信号由用户终端发出,携带相关的信息传回子站,其调制信号的频率在150Hz~600Hz范围内变化,找出150HZ~500HZ频率所在的阶数进行保留,保留的150HZ~500HZ所在的范围为调制信号所在阶数;
5.对保留的IMF进行阈值去噪:
由于电网中的噪声都是随机的,其能量大部分是均匀分布,有用信号的能量相对集中,根据有用信号和噪声信号能量的分布不同选择合适的阈值对保留的IMF进行降噪处理;
6.计算各阶IMF希尔伯特边际谱:
获得IMF分量后,对于其中一个IMF分量记为imf(t),经过希尔伯变换后的H[imf(t)]与原分量imf(t)组成解析信号,记为其中 瞬时角频率为以灰度形式在频率-时间平面上显示的希尔伯特谱,对H(ω,t)在时间上积分,进而求得希尔伯特边际谱
根据得到的希尔伯特边际谱可以清晰的判断出调制信号所在的时域。
本发明采用的EEMD是基于任何信号均由不同时间尺度振荡构成,所以能够对非线性、非平稳信号自适应的分解成为一系列线性、平稳的IMF信号,在对信号进行时频分析时与经常使用的小波变换相比避免了受小波基选择的影响,不仅能消除大部分噪声干扰而且能很好的抑制谐波干扰,在检测波动信号时具有较好的动态特性,由于双向工频通信上行调制信号本身特性不适合利用EEMD,双向工频通信上行调制信号是一种经晶闸管的开关在工频电压过零点处产生的电流脉冲,每个脉冲含有一个极值点和两个过零点,对这样的信号进行EEMD计算量会非常大,对其求导后满足IMF的两个条件大大减小了计算量;对其进行求导使之更好的适应EEMD,这种检测算法更加精确。
本发明所描述的具体方式仅为一个实例,并非穷举,本领域技术人员可以对所描述的具体实例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离发明的精神或者超越权利要求所记载的保护范围。

Claims (2)

1.一种电力线双向工频通信上行信号检测方法,其特征在于,它包括如下步骤:
1)对信号加权求和消除背景信号:
工频通信上行信号检测是在大背景中提取微弱电流,所以利用加权求和抵消工频信号,利用其中IA,IB,IC是配电网中的三相电流,能消除绝大部分的基波电流影响,同时使信号增强1倍;
2)对加权求和后的通信信号求导:
双向工频通信上行调制信号是经晶闸管的开关在工频电压过零点处产生的电流脉冲,每个脉冲含有一个极大值或一个极小值以及两个过零点,对其求导后,每个脉冲含有一个极大值,一个极小值和一个过零点,满足EEMD的分解条件,所述EEMD为集合经验模分解;
3)对处理后的信号进行EEMD:
在原始信号中加入均匀分布的高斯白噪声用于填补原始调制畸变电流中部分尺度缺失,使待测的调制电流在不同程度上连续;选定进行EMD的分解次数和白噪声的幅值系数,分解次数取100,幅值系数取0.2,所述EMD为经验模分解;
4)测量各本征模函数频率:
通过傅里叶变换测量各阶IMF的频率,找出150HZ~500HZ频率所在的阶数进行保留,所述IMF为本征模函数;
5)对保留的IMF进行阈值去噪:
由于电网中的噪声都是随机的,其能量大部分是均匀分布,有用信号的能量相对集中,根据有用信号和噪声信号能量的分布不同选择合适的阈值对保留的IMF进行降噪处理;
6)计算各阶IMF希尔伯特边际谱:
对得到的IMF分量其中一个本征模函数分量记为imf(t),经过希尔伯变换后的H[imf(t)]与原imf(t)分量组成解析信号,记为其中 瞬时角频率为H(ω,t)是振幅I(t)以灰度形式在频率-时间平面上显示的希尔伯特谱,对H(ω,t)在时间上积分,进而求得希尔伯特边际谱 h ( ω ) = ∫ 0 T H ( ω , t ) dt .
2.按照权利要求1所述的一种电力线工频通信上行信号检测方法,其特征在于,所述步骤3)对处理后的信号进行EEMD步骤如下:
(1)在原调制畸变电流信号中加入幅值系数为0.2的高斯白噪声,对所得信号进行EMD,得到n个IMF记为imfj
(2)对原调制畸变电流信号添加的下一个白噪声序列重复上述步骤,获得各自的IMF;
(3)计算剩余信号Ri,1=xi(t)-imfi,1(t),xi(t)为原始信号,imfi,1(t)为混入第i个白噪声序列分解得到的第一阶IMF,把剩余信号作为待处理信号进行EMD分解,依次获得imfi,1,imfi,2,....imfi,j
(4)对所得结果求平均其中j=1,2,…n,n为分解的阶数,去除添加的白噪声,得到待测信号。
CN201410597083.7A 2014-10-30 2014-10-30 一种电力线双向工频通信上行信号检测方法 Expired - Fee Related CN104320157B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410597083.7A CN104320157B (zh) 2014-10-30 2014-10-30 一种电力线双向工频通信上行信号检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410597083.7A CN104320157B (zh) 2014-10-30 2014-10-30 一种电力线双向工频通信上行信号检测方法

Publications (2)

Publication Number Publication Date
CN104320157A true CN104320157A (zh) 2015-01-28
CN104320157B CN104320157B (zh) 2016-08-24

Family

ID=52375351

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410597083.7A Expired - Fee Related CN104320157B (zh) 2014-10-30 2014-10-30 一种电力线双向工频通信上行信号检测方法

Country Status (1)

Country Link
CN (1) CN104320157B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108267638A (zh) * 2016-12-30 2018-07-10 钜泉光电科技(上海)股份有限公司 基于工频畸变的工频同步通信信号检测装置、系统和方法
CN108768452A (zh) * 2018-05-30 2018-11-06 国网山东省电力公司经济技术研究院 一种电力线双向工频通信的检测方法
CN109257156A (zh) * 2018-09-04 2019-01-22 国网山东省电力公司经济技术研究院 基于时频分析的提取twacs通信信号特征的方法
CN110061762A (zh) * 2019-03-27 2019-07-26 国网山东省电力公司青岛供电公司 基于决策融合的电力线工频通讯信号检测算法
CN111103489A (zh) * 2019-12-02 2020-05-05 国网吉林省电力有限公司吉林供电公司 一种基于工频通信技术的配电变压器运行工况在线监测系统
CN114019211A (zh) * 2021-10-14 2022-02-08 国网福建省电力有限公司福州供电公司 一种畸变电流采样方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080055067A1 (en) * 2006-09-01 2008-03-06 Walter Curt Method and apparatus for a power line communications test system
CN101980451A (zh) * 2010-10-09 2011-02-23 华北电力大学 基于时段自适应的电力线工频通信系统及方法
CN203071952U (zh) * 2012-12-07 2013-07-17 四川电力科学研究院 一种用电信息采集终端的上下行通信测试装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080055067A1 (en) * 2006-09-01 2008-03-06 Walter Curt Method and apparatus for a power line communications test system
CN101980451A (zh) * 2010-10-09 2011-02-23 华北电力大学 基于时段自适应的电力线工频通信系统及方法
CN203071952U (zh) * 2012-12-07 2013-07-17 四川电力科学研究院 一种用电信息采集终端的上下行通信测试装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈晓娟等: "一种实用的工频通信上行信号检测方法", 《电测与仪表》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108267638A (zh) * 2016-12-30 2018-07-10 钜泉光电科技(上海)股份有限公司 基于工频畸变的工频同步通信信号检测装置、系统和方法
CN108267638B (zh) * 2016-12-30 2021-08-10 钜泉光电科技(上海)股份有限公司 基于工频畸变的工频同步通信信号检测装置、系统和方法
CN108768452A (zh) * 2018-05-30 2018-11-06 国网山东省电力公司经济技术研究院 一种电力线双向工频通信的检测方法
CN108768452B (zh) * 2018-05-30 2021-08-24 国网山东省电力公司经济技术研究院 一种电力线双向工频通信的检测方法
CN109257156A (zh) * 2018-09-04 2019-01-22 国网山东省电力公司经济技术研究院 基于时频分析的提取twacs通信信号特征的方法
CN110061762A (zh) * 2019-03-27 2019-07-26 国网山东省电力公司青岛供电公司 基于决策融合的电力线工频通讯信号检测算法
CN110061762B (zh) * 2019-03-27 2021-12-03 国网山东省电力公司青岛供电公司 基于决策融合的电力线工频通讯信号检测算法
CN111103489A (zh) * 2019-12-02 2020-05-05 国网吉林省电力有限公司吉林供电公司 一种基于工频通信技术的配电变压器运行工况在线监测系统
CN114019211A (zh) * 2021-10-14 2022-02-08 国网福建省电力有限公司福州供电公司 一种畸变电流采样方法及系统

Also Published As

Publication number Publication date
CN104320157B (zh) 2016-08-24

Similar Documents

Publication Publication Date Title
CN104320157A (zh) 一种电力线双向工频通信上行信号检测方法
Khamis et al. Islanding detection in a distributed generation integrated power system using phase space technique and probabilistic neural network
Achlerkar et al. Variational mode decomposition and decision tree based detection and classification of power quality disturbances in grid-connected distributed generation system
Ye et al. An improved fault-location method for distribution system using wavelets and support vector regression
Zhang et al. EMTR‐based fault location for DC line in VSC‐MTDC system using high‐frequency currents
CN103344875A (zh) 谐振接地系统单相接地故障分类选线方法
Soto et al. Electromagnetic field due to lightning striking on top of a cone-shaped mountain using the FDTD
CN106533591B (zh) 一种小型基站电磁干扰识别方法
Huai et al. Single-ended line fault location method for multi-terminal HVDC system based on optimized variational mode decomposition
Guo et al. Hilbert–Huang transform‐based transient busbar protection algorithm
CN102095929B (zh) 一种快速测量交流电信号频率的方法
CN108333564A (zh) 用于谐波雷达频谱感知与频率选择的方法
He et al. The improved fault location method based on natural frequency in MMC-HVDC grid by combining FFT and MUSIC algorithms
CN105137181A (zh) 基于Nuttall-Kaiser组合窗双谱线插值的谐波分析方法
CN102901855B (zh) 一种特高压直流电晕电流信号去噪方法
CN107390022A (zh) 基于离散频谱校正的电能计量方法
Guo et al. Calculation of lightning induced voltages on overhead lines using an analytical fitting representation of electric fields
CN108920760A (zh) 基于多项式混沌展开的天线高空电磁脉冲耦合响应性能统计方法
CN106771586A (zh) 一种直流控制保护板卡的回路信号分析方法及装置
Da Silva et al. A hybrid fault locator for three-terminal lines based on wavelet transforms
CN102087311B (zh) 一种提高电力互感器测量精度的方法
CN108008187A (zh) 基于变分模态分解的电网谐波检测方法
Xie et al. Altitude correction of radio interference of HVDC transmission lines part I: converting method of measured data
Li et al. Research on lightning electromagnetic fields associated with first and subsequent return strokes based on Laplace wavelet
CN103618287A (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
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160824

Termination date: 20171030

CF01 Termination of patent right due to non-payment of annual fee