CN108152679A - 一种基于hht与闪变功率的闪变源定位方法 - Google Patents
一种基于hht与闪变功率的闪变源定位方法 Download PDFInfo
- Publication number
- CN108152679A CN108152679A CN201711447962.1A CN201711447962A CN108152679A CN 108152679 A CN108152679 A CN 108152679A CN 201711447962 A CN201711447962 A CN 201711447962A CN 108152679 A CN108152679 A CN 108152679A
- Authority
- CN
- China
- Prior art keywords
- flicker
- signal
- voltage
- power
- flickering
- 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
Links
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/088—Aspects of digital computing
-
- 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/086—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于HHT与闪变功率的闪变源定位方法,包括以下步骤:1)对从PCC点处检测到的含有闪变的电压、电流信号,使用Hilbert变换将含有闪变的电压、电流信号的包络信号提取出来;2)使用EMD分解方法对电压、电流信号的包络信号进行处理,得到一组仅含单一分量的闪变信号;3)将分解出的相同频率分量的电压、电流闪变信号依据闪变功率的数学公式计算出该频率下的闪变功率值;4)按照闪变功率法的闪变源定位标准对各个频率的闪变源进行定位识别。本发明有效克服了闪变功率法仅能对PCC点处的主导闪变源进行定位识别的问题,可以对不同频率的闪变源同时开展定位识别工作。
Description
技术领域
本发明涉及一种基于HHT与闪变功率的闪变源定位方法,属于信号处理技术领域。
背景技术
随着电网中新能源的接入以及大量冲击性负荷的并网应用,致使电压波动现象频繁发生,由此引起的电压闪变问题对电网中的电能质量造成了较大影响。闪变源定位识别的研究,既有助于发现电压闪变产生的根本原因,又可以有效推动电能质量问题的责任认定工作。
闪变源定位识别研究,目前主要采用的方法为间谐波功率流向法,电压电流比率法,闪变功率流向法。间谐波功率流向法通过间谐波功率流向对闪变源进行判别,但是谐波源有时会向系统中注入谐波功率,有时也会从系统中吸收谐波功率,因此,依据间谐波功率流向法对闪变源进行定位识别存在准确性不高的问题。电压电流比率法通过电压电流的比值正负情况对闪变源进行定位识别,由于改方法采用的电压和电流值是有效值,因此对系统中出现的随机干扰具有一定的抑制作用,但由于电力系统中的电压与电流存在一定的暂态性,因此无法确保可靠性。闪变功率流向法类似间谐波功率流向法,但是该方法只能对公共连接点处的主导闪变源进行定位识别,无法对所有频率的闪变源进行定位识别。
发明内容
本发明的主要目的在于,克服现有技术中的不足,提供一种基于HHT与闪变功率的闪变源定位方法,可以对公共连接点处多种闪变频率的闪变源进行定位识别。
本发明所采用的技术方案是:
一种基于HHT(Hilbert-Huang Transform,希尔伯特-黄变换)与闪变功率的闪变源定位方法,包括以下步骤:
步骤1),在公共连接点(Point of Common Coupling,PCC)处截取信号序列,从截取的信号序列中检测出含有闪变的电压信号和电流信号,使用Hilbert变换将含有闪变的电压信号和电流信号的包络信号提取出来;
步骤2),运用经验模态分解(Empirical Mode Decomposition,EMD)方法对电压信号和电流信号的包络信号进行处理,分解出一组仅含单一频率分量的电压闪变信号和电流闪变信号;
步骤3),将分解出的相同频率分量的电压闪变信号和电流闪变信号依据闪变功率的数学公式计算出所述频率下的闪变功率值;
步骤4),按照闪变功率法的闪变源定位标准对各个频率的闪变源进行定位识别。
较优地,步骤1)Hilbert变换具体为,
式中:X(t)为从PCC点处检测的电压信号或者电流信号(X(t)可以为电压信号也可以为电流信号),Hilbert变换能够将含有闪变的电压信号和电流信号的包络信号提取出来,当X(t)输入的为电压信号时,Hilbert变换将含有闪变的电压信号的包络信号提取出来,当X(t)输入的为电流信号时,Hilbert变换将含有闪变的电流信号的包络信号提取出来;Y(t)为X(t)的Hilbert变换,*表示卷积;t表示时间;τ表示积分分量;a(t)为信号的幅值;hp(t)为系统的脉冲响应。
较优地,步骤2)中具体包括以下步骤:
令h0(t)=X(t),在截取的信号序列X(t)中找出所有的极值点,采用三次样条插值算法进行极值点的拟合,得到极值点对应的信号包络线,信号包络线包括上包络线和下包络线,将上包络线和下包络线进行相加后除以2,得到信号包络线的均值mk(t),将原信号hk(t)减去均值mk(t),得到hk+1(t),依据停止条件对hk+1(t)进行判断,如果不满足停止条件,则令k=k+1,k=0,1,2…,循环从信号hk+1(t)中找出所有的极值点,重新进行信号包络线提取;如果满足停止条件则分解结束,分解出了一组仅含单一频率分量的电压闪变信号或者电流闪变信号。
较优地,闪变功率值的计算公式为:
式中,Π(t)为闪变功率值,t表示时间,T为测量闪变的电压幅值与电流幅值的积分周期,mu(τ)、mi(τ)分别表示某种频率下闪变电压、闪变电流的实时波动分量。
较优地,步骤4)具体包括以下步骤,
当PCC点处的闪变功率值Π(t)为正时,闪变源处于以PCC点为基准点的系统侧,当PCC点处的闪变功率值Π(t)为负时,闪变源处于以PCC点为基准点的负荷侧。
较优地,停止条件具体为:令
式中,T为信号时间长度;当SD小于设定阈值时,则满足停止条件。
较优地,设定阈值为0.2~0.3。
与现有技术相比,本发明具有的有益效果是:
本文提出了一种基于HHT与闪变功率的闪变源定位方法,在闪变功率法能够对PCC点处主导闪变源的进行定位的前提下,利用HHT分解对闪变功率计算方法进行改进,通过将检测到的节点闪变电压和支路闪变电流进行分解,得到含单一频率的闪变电压信号和闪变电流信号,计算每个频率下对应的闪变功率值,使用闪变功率的闪变源定位规则对各个频率的闪变源进行定位识别。
本发明解决了闪变功率法只能对主导闪变源进行定位的缺点,使现在可以对所有频率的闪变源进行定位识别了。
附图说明
图1为本发明一种基于HHT与闪变功率的闪变源定位方法流程图;
图2为本发明的一种含有单闪变源的配电网系统图;
图3为本发明经验模态分解流程图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
如图1所示,本实施例公开一种基于HHT与闪变功率的闪变源定位方法,包括以下步骤:
步骤1)、公共连接点(Point of Common Coupling,PCC)处截取信号序列,从截取的信号序列中检测出含有闪变的电压信号和电流信号,使用Hilbert变换将含有闪变的电压信号和电流信号的包络信号提取出来;
Hilbert变换具体为,
式中:X(t)为从PCC点处检测的电压信号或者电流信号,Y(t)为X(t)的Hilbert变换,*表示卷积;t表示时间;τ表示积分分量;a(t)为信号的幅值;hp(t)为系统的脉冲响应。Hilbert变换能够将含有闪变的电压信号和电流信号的包络信号提取出来,当X(t)输入的为电压信号时,Hilbert变换将含有闪变的电压信号的包络信号提取出来,当X(t)输入的为电流信号时,Hilbert变换将含有闪变的电流信号的包络信号提取出来。
步骤2)、运用经验模态分解(Empirical Mode Decomposition,EMD)方法对电压信号和电流信号的包络信号进行处理,分解出一组仅含单一频率分量的电压闪变信号和电流闪变信号;
步骤2)中具体包括以下步骤:
令h0(t)=X(t)(初始化),在截取的信号序列X(t)中找出所有的极值点,采用三次样条插值算法进行极值点的拟合,得到极值点对应的信号包络线,信号包络线包括上包络线和下包络线,将上包络线和下包络线进行相加后除以2,得到信号包络线的均值mk(t),将原信号hk(t)减去均值mk(t),得到hk+1(t),依据停止条件对hk+1(t)进行判断,如果不满足停止条件,则令k=k+1,k=0,1,2…,循环从信号hk+1(t)中找出所有的极值点,重新进行信号包络线提取;如果满足停止条件则分解结束,分解出了一组仅含单一频率分量的电压闪变信号和电流闪变信号。
停止条件具体为:令
式中,T为信号时间长度;当SD小于设定阈值时,则满足停止条件。
设定阈值为0.2~0.3,本实施例,设定阈值为0.25。
步骤3)、将分解出的相同频率分量的电压闪变信号和电流闪变信号依据闪变功率的数学公式计算出所述频率下的闪变功率值;
闪变功率值的计算公式为:
式中,Π(t)为闪变功率值,t表示时间,T为测量闪变的电压幅值与电流幅值的积分周期,mu(τ)、mi(τ)分别表示某种频率下闪变电压、闪变电流的实时波动分量。
步骤4)、按照闪变功率法的闪变源定位标准对各个频率的闪变源进行定位识别。
当PCC(Point of Common Coupling,PCC)点处的闪变功率值Π(t)为正时,闪变源处于以PCC点为基准点的系统侧,当PCC点处的闪变功率值Π(t)为负时,闪变源处于以PCC点为基准点的负荷侧。
本实施例如图2为一种含有单闪变源的配电网系统图,表1为图2中含多闪变源的配网系统参数,表2为本实施例计算出含有多闪变源的配网各节点功率。
表1为本实施例含多闪变源的配网系统参数
表2为本实施例计算出含有多闪变源的配网各节点功率
表2为节点1,节点2,节点3处计算得到的闪变功率值,从中可以看出,节点1,节点3处闪变频率为10hz与15hz的闪变功率为正值,而节点2处闪变频率为10hz与15hz的闪变功率为负值,从而可以判断节点2处有闪变频率为10hz与15hz的闪变源。
本发明的创新点在于,针对应用闪变功率流向法仅能对公共连接点处的主导闪变源进行定位识别,从而采用HHT算法对从公共连接点处检测到的电压与分支处检测到的电流进行分解处理,从而得到各个频率的电压闪变和电流闪变信号,随后通过闪变功率数学公式计算各个频率下的闪变功率值。最后使用闪变功率的闪变源定位规则对各个频率的闪变源进行定位识别。
以上仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (7)
1.一种基于HHT与闪变功率的闪变源定位方法,其特征在于,包括以下步骤:
步骤1)、公共连接点处截取信号序列,从截取的信号序列中检测出含有闪变的电压信号和电流信号,使用Hilbert变换将含有闪变的电压信号和电流信号的包络信号提取出来;
步骤2)、运用经验模态分解方法对电压信号和电流信号的包络信号进行处理,分解出一组仅含单一频率分量的电压闪变信号和电流闪变信号;
步骤3)、将分解出的相同频率分量的电压闪变信号和电流闪变信号依据闪变功率的数学公式计算出所述频率下的闪变功率值;
步骤4)、按照闪变功率法的闪变源定位标准对各个频率的闪变源进行定位识别。
2.根据权利要求1所述的一种基于HHT与闪变功率的闪变源定位方法,其特征在于:
步骤1)Hilbert变换具体为,
式中:X(t)为从PCC点处检测的电压信号或者电流信号,Y(t)为X(t)的Hilbert变换,*表示卷积,t表示时间,τ表示积分分量,a(t)为信号的幅值,hp(t)为系统的脉冲响应。
3.根据权利要求1所述的一种基于HHT与闪变功率的闪变源定位方法,其特征在于:
步骤2)具体包括以下步骤:
令hk(t)=X(t),在截取的信号序列X(t)中找出所有的极值点,采用三次样条插值算法进行极值点的拟合,得到极值点对应的信号包络线,信号包络线包括上包络线和下包络线,将上包络线和下包络线进行相加后除以2,得到信号包络线的均值mk(t),将原信号hk(t)减去均值mk(t),得到hk+1(t),依据停止条件对hk+1(t)进行判断,如果不满足停止条件,则令k=k+1,k=0,1,2…,循环从信号hk+1(t)中找出所有的极值点,重新进行信号包络线提取;如果满足停止条件则分解结束,分解出了一组仅含单一频率分量的电压闪变信号或者电流闪变信号。
4.根据权利要求1所述的一种基于HHT与闪变功率的闪变源定位方法,其特征在于:
闪变功率值的计算公式为:
式中,Π(t)为闪变功率值,t表示时间,T为测量闪变的电压幅值与电流幅值的积分周期,mu(τ)、mi(τ)分别表示闪变电压、闪变电流的实时波动分量。
5.根据权利要求1所述的一种基于HHT与闪变功率的闪变源定位方法,其特征在于:
所述步骤4)具体包括以下步骤,
当PCC点处的闪变功率值Π(t)为正时,闪变源处于以PCC点为基准点的系统侧,当PCC点处的闪变功率值Π(t)为负时,闪变源处于以PCC点为基准点的负荷侧。
6.根据权利要求3所述的一种基于HHT与闪变功率的闪变源定位方法,其特征在于:
所述停止条件具体为:令
式中,T为信号时间长度;当SD小于设定阈值时,则满足停止条件。
7.根据权利要求6所述的一种基于HHT与闪变功率的闪变源定位方法,其特征在于:
设定阈值为0.2~0.3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711447962.1A CN108152679A (zh) | 2017-12-27 | 2017-12-27 | 一种基于hht与闪变功率的闪变源定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711447962.1A CN108152679A (zh) | 2017-12-27 | 2017-12-27 | 一种基于hht与闪变功率的闪变源定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108152679A true CN108152679A (zh) | 2018-06-12 |
Family
ID=62462228
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711447962.1A Pending CN108152679A (zh) | 2017-12-27 | 2017-12-27 | 一种基于hht与闪变功率的闪变源定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108152679A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI668454B (zh) * | 2018-08-06 | 2019-08-11 | 國立中正大學 | 閃爍分析裝置及方法 |
WO2020087556A1 (zh) * | 2018-10-30 | 2020-05-07 | 重庆先进光电显示技术研究院 | 一种快速寻找显示面板的公共电压的方法及测试机台 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1742208A (zh) * | 2002-12-23 | 2006-03-01 | 联合动力股份公司 | 用于确定闪变源方向的测量方法 |
CN101055299A (zh) * | 2007-01-31 | 2007-10-17 | 湖南大学 | 一种冲击性负荷电能计量方法 |
CN101598760A (zh) * | 2009-07-01 | 2009-12-09 | 华北电力大学(保定) | 一种闪变源定向的分形分析方法 |
JP4531636B2 (ja) * | 2005-06-13 | 2010-08-25 | 東芝三菱電機産業システム株式会社 | フリッカ抑制効果シミュレーション装置 |
CN102621454A (zh) * | 2012-04-06 | 2012-08-01 | 河海大学 | 谐波源闪变源实时定位仪及确定谐波源闪变源方位的方法 |
CN106779470A (zh) * | 2017-01-04 | 2017-05-31 | 南京工程学院 | 一种基于改进hht的电压闪变检测算法 |
-
2017
- 2017-12-27 CN CN201711447962.1A patent/CN108152679A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1742208A (zh) * | 2002-12-23 | 2006-03-01 | 联合动力股份公司 | 用于确定闪变源方向的测量方法 |
JP4531636B2 (ja) * | 2005-06-13 | 2010-08-25 | 東芝三菱電機産業システム株式会社 | フリッカ抑制効果シミュレーション装置 |
CN101055299A (zh) * | 2007-01-31 | 2007-10-17 | 湖南大学 | 一种冲击性负荷电能计量方法 |
CN101598760A (zh) * | 2009-07-01 | 2009-12-09 | 华北电力大学(保定) | 一种闪变源定向的分形分析方法 |
CN102621454A (zh) * | 2012-04-06 | 2012-08-01 | 河海大学 | 谐波源闪变源实时定位仪及确定谐波源闪变源方位的方法 |
CN106779470A (zh) * | 2017-01-04 | 2017-05-31 | 南京工程学院 | 一种基于改进hht的电压闪变检测算法 |
Non-Patent Citations (5)
Title |
---|
HADI MOGHADAM BANAYEM等: "Flicker Source Tracing byWavelet Transform", 《ELECTRIC POWER COMPONENTS AND SYSTEMS》 * |
刘稳坚 等: "经验模态分解及Hilbert谱分析在电压闪变检测中的应用", 《计算机测量与控制》 * |
李天云 等: "基于HHT的电能质量检测新方法", 《中国电机工程学报》 * |
王志群 等: "Hilbert变换求取电压闪变有关参数", 《电力系统自动化》 * |
王畅: "基于HHT的微网电能质量检测与分析", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI668454B (zh) * | 2018-08-06 | 2019-08-11 | 國立中正大學 | 閃爍分析裝置及方法 |
WO2020087556A1 (zh) * | 2018-10-30 | 2020-05-07 | 重庆先进光电显示技术研究院 | 一种快速寻找显示面板的公共电压的方法及测试机台 |
US10891885B2 (en) | 2018-10-30 | 2021-01-12 | Chongqing Advance Display Technology Research | Method and test machine platform for quickly searching for common voltage of display panel |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Multiscale adaptive fault diagnosis based on signal symmetry reconstitution preprocessing for microgrid inverter under changing load condition | |
CN104502700A (zh) | 电力信号的正弦参数测量方法和系统 | |
TW201720002A (zh) | 紀錄供應網路中電壓之方法及裝置 | |
CN104362628B (zh) | 一种避免大型风电场谐振的控制方法及其装置 | |
CN104319794B (zh) | 一种提高系统稳定性的交直流协调控制方法 | |
CN105811408B (zh) | 面向mmc换流站的谐波能量传递函数建模方法 | |
CN108152679A (zh) | 一种基于hht与闪变功率的闪变源定位方法 | |
Liu et al. | A new classification method for transient power quality combining spectral kurtosis with neural network | |
CN104360235A (zh) | 含dg复杂配电网的电压暂降源定位方法 | |
CN108418238A (zh) | 基于谐波综合换相系数评估换相失败发生的方法 | |
CN102998527A (zh) | 通频型基波、谐波、直流分量检测方法 | |
CN108008187A (zh) | 基于变分模态分解的电网谐波检测方法 | |
CN107045081B (zh) | 一种保护装置中实时频率的测量方法 | |
CN110749769A (zh) | 一种三相电压跌落的快速检测方法 | |
Meyur et al. | A wavelet-adaptive network based fuzzy inference system for location of faults in parallel transmission lines | |
Dubey et al. | A novel out-of-step detection algorithm using angle of power signal | |
CN104868499A (zh) | 三相不平衡情况下的基于自适应陷波器的三相pq控制方法及系统 | |
CN110098620A (zh) | 一种用于优化换流站电压的控制方法、装置及系统 | |
CN114825415A (zh) | 一种柔性直流高频谐振的自适应抑制方法及系统 | |
Otchere et al. | A 2 nd order LPF wavelet based control scheme for shunt active power filter | |
CN107359839A (zh) | 一种静止变频器启动同步电机过程中机组磁链计算方法 | |
Satyanarayana et al. | Identification of Faults in HVDC System using Wavelet Analysis | |
Wu et al. | A novel active islanding fault detection based on even harmonics injection and set-membership filtering | |
CN111707856A (zh) | 一种利用波形叠加原理设计的电力系统故障电流测量方法 | |
CN203396841U (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180612 |