CN114236235B - 一种采用相似特征匹配的局部不可观系统谐波源定位方法 - Google Patents
一种采用相似特征匹配的局部不可观系统谐波源定位方法 Download PDFInfo
- Publication number
- CN114236235B CN114236235B CN202111554443.1A CN202111554443A CN114236235B CN 114236235 B CN114236235 B CN 114236235B CN 202111554443 A CN202111554443 A CN 202111554443A CN 114236235 B CN114236235 B CN 114236235B
- Authority
- CN
- China
- Prior art keywords
- harmonic
- equivalent
- brh
- node
- harmonic current
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000004364 calculation method Methods 0.000 claims abstract description 11
- 238000001914 filtration Methods 0.000 claims abstract description 8
- 230000009471 action Effects 0.000 claims abstract description 7
- 238000012417 linear regression Methods 0.000 claims abstract description 7
- 238000011549 displacement method Methods 0.000 claims abstract description 4
- 238000012216 screening Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 36
- 238000005259 measurement Methods 0.000 claims description 22
- 230000008569 process Effects 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
- G01R23/163—Spectrum analysis; Fourier analysis adapted for measuring in circuits having distributed constants
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/40—Arrangements for reducing harmonics
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种采用相似特征匹配的局部不可观系统谐波源定位方法,包括以下步骤:基于多端口网络等效原理,采用谐波电流移置法,将局部不可观系统转化为等效可观系统;建立谐波源定位的等值模型,通过卡尔曼滤波算法求解谐波电流动态信息,以此确定谐波源疑似存在区域;通过相关性计算,筛选不可观区域内主谐波源作用下的谐波电流状态量,并采用线性回归辨识边界节点状态量之间的回归系数,将其与区域内节点对应的等效参数进行匹配,将匹配度最高的节点视为主谐波源所在位置。该方法有利于对局部不可观系统中的谐波源进行准确定位。
Description
技术领域
本发明属于谐波状态估计技术领域,具体涉及一种采用相似特征匹配的局部不可观系统谐波源定位方法。
背景技术
近些年,随着大量分布式电源接入电网中,谐波污染问题日益严重。谐波源定位是划分谐波源责任和采用经济手段惩治谐波问题的基础,具有重要的研究意义。然而,分布式电源具有强随机性和波动性,仅通过单一时刻断面下的量测量难以充分反映出谐波特性;此外,量测装置因成本限制,无法做到对电网信息的全面监测,系统的不完全可观性进一步加剧了谐波源定位的难度。因此,有必要研究在有限量测装置数量的前提下,如何有效利用谐波电流的动态信息,实现谐波源的准确定位。
对局部不可观系统中的谐波源进行定位,其模型通常为欠定方程,难以通过常规方法求解。目前,大多谐波源定位方法主要围绕谐波源的稀疏特性,采用压缩感知等算法求解欠定模型,重构出某一时间断面下的谐波源信号,但该类方法容易忽略谐波的波动特征,且易受测点位置等因素的影响,难以保证重构结果的准确度,容易对谐波源位置造成误判。另外,目前基于动态谐波状态估计的谐波源定位方法,通常需要保证系统完全可观的前提下才能够顺利完成对谐波源定位的目标,因此无法直接适用于局部不可观系统。
发明内容
本发明的目的在于提供一种采用相似特征匹配的局部不可观系统谐波源定位方法,该方法有利于对局部不可观系统中的谐波源进行准确定位。
为实现上述目的,本发明采用的技术方案是:一种采用相似特征匹配的局部不可观系统谐波源定位方法,包括以下步骤:
基于多端口网络等效原理,采用谐波电流移置法,将局部不可观系统转化为等效可观系统;
建立谐波源定位的等值模型,通过卡尔曼滤波算法求解谐波电流动态信息,以此确定谐波源疑似存在区域;
通过相关性计算,筛选不可观区域内主谐波源作用下的谐波电流状态量,并采用线性回归辨识边界节点状态量之间的回归系数,将其与区域内节点对应的等效参数进行匹配,将匹配度最高的节点视为主谐波源所在位置。
进一步地,建立谐波源定位的等值模型的具体方法如下:
化简电力系统内的不可观节点,保留相关可观节点,将局部不可观系统转化为等效可观系统;
由此得到所述局部不可观系统化简后的表达式,记为:
其中,为等效后的节点导纳矩阵,/>为保留节点的电压,/>为原注入电流在边界节点上产生的等效电流;
进一步得到谐波源定位的等值模型如下:
式中,h表示谐波频次,为谐波量测量,含节点谐波电压与支路谐波电流;为系统等效后的谐波量测矩阵;/>为谐波电流;/>为谐波量测的噪声;/>表示复数域,上标m表示量测量数,上标n表示状态量数;
求解上述模型,当不可观区域内存在谐波源时,将在边界节点产生等效谐波电流。
进一步地,求解谐波电流动态信息的具体方法如下:
根据谐波源定位等值模型中连续多个时刻的量测量,采用卡尔曼滤波算法对谐波的状态量进行迭代更新计算;
在预测阶段,谐波电流的预测值和先验误差通过下式获取:
其中,上标-表示先验估计,上标+表示后验估计,下标k表示k时刻,下标k+1表示k+1时刻,即下一时刻,上标^表示估计值;A为状态转移矩阵;为k+1时刻的先验误差矩阵;为k时刻的后验误差矩阵;Q为系统噪声;/>为k时刻的后验状态估计值;/>表示k+1时刻的先验状态估计值,指谐波电流的预测值;
在校正阶段,谐波电流的最优估计值和后验误差通过下式获取:
其中,上标-1为矩阵求逆;R为测量噪声;H为量测矩阵;Kk+1为k+1时刻的增益矩阵;zk+1为k+1时刻的观测值;G为单位矩阵;为k+1时刻的后验误差矩阵;/>为k+1时刻的后验状态估计值,指校正后的谐波电流最优估计值;
谐波电流的动态信息围绕式(3)至(7)进行迭代更新;式(3)预测k+1时刻状态量;式(4)则表示真实值与预测状态量之间的误差;式(5)通过计算增益矩阵来校正更新式(6)中k+1时刻最优的状态量;式(7)更新最优状态量与真实值之间的误差,以执行下一次迭代。
进一步地,基于相似特征匹配定位谐波源,其具体方法如下:
定义不可观区域内存在主谐波源,即其余谐波源保持相对稳定,此时边界节点的等效谐波电流受主谐波源影响,边界节点谐波电流的波动特征同时取决于该主谐波源的波动特征,与该主谐波源的波动特征具有相似性,进而可知边界节点谐波电流之间存在一定的联系,因此进一步描述边界节点之间的等效谐波电流一般线性关系如下:
式中,为边界节点l的等效谐波电流;/>为边界节点r的等效谐波电流变量;为回归系数,对应主谐波源节点对应的等效参数;/>为常量;
对边界节点的谐波电流进行窗口划分,设定该划分的窗宽为n,窗口数为m,对每个窗口进行皮尔逊相关性计算,并选择相关性大于阈值α的谐波电流窗口,此时为主谐波源作用下的窗口;
对窗口内谐波电流进行线性回归,并辨识其回归系数,回归过程如下所示:
Ibrh,l,real=kbrh,r,real,1Ibrh,r,real+kbrh,r,imag,1Ibrh,r,imag+b1 (9)
Ibrh,l,imag=kbrh,r,imag,2Ibrh,r,real+kbrh,r,real,2Ibrh,r,imag+b2 (10)
式中,Ibrh,r,real为边界节点r的等效谐波电流实部,Ibrh,r,imag为边界节点r的等效谐波电流虚部,Ibrh,l,real为边界节点l的等效谐波电流实部,Ibrh,l,imag为边界节点l的等效谐波电流虚部,b为常量;
从上式总共得到四个变量系数,即kbrh,r,real,1、kbrh,r,imag,1、kbrh,r,real,2、kbrh,r,imag,2;将上述四个变量系数与区域内节点对应的等效参数进行匹配对比,定义匹配度评估指标gmatch如下所示:
式中,ktrue,real、ktrue,imag即式(8)节点对应的等效参数实部与虚部;
gmatch越大,表明回归系数越接近节点对应的等效参数,主谐波源处于该节点的可能性越大;选取匹配度最高的节点,该节点即为主谐波源所处位置;
当相关性大于阈值α的谐波电流窗口不存在时,判定该局部不可观区域不存在主谐波源,只进行谐波源区域定位。
与现有技术相比,本发明具有以下有益效果:提供了一种采用相似特征匹配的局部不可观系统谐波源定位方法,该方法利用网络化简手段,将局部不可观系统转化为等效可观系统,从而使欠定模型转化为正定或超定模型,然后采用卡尔曼滤波算法获取谐波动态信息,通过相关性计算筛选出高度线性相关的状态量,并根据边界节点之间的回归系数确定不可观区域内部主谐波源位置。该方法实现了在局部不可观系统中利用谐波动态特性来对谐波源进行准确定位。
附图说明
图1是本发明实施例的方法实现流程图。
图2是本发明实施例中局部不可观系统转化为等效可观系统的示意图。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
应该指出,以下详细说明都是示例性的,旨在对本申请提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
如图1所示,本实施例提供了一种采用相似特征匹配的局部不可观系统谐波源定位方法,包括以下步骤:
1)基于多端口网络等效原理,采用谐波电流移置法,将局部不可观系统转化为等效可观系统。
2)建立谐波源定位的等值模型,通过卡尔曼滤波算法求解谐波电流动态信息,以此确定谐波源疑似存在区域。
3)通过皮尔逊相关性计算,筛选不可观区域内主谐波源作用下的谐波电流状态量,并采用线性回归辨识边界节点状态量之间的回归系数,将其与区域内节点对应的等效参数进行匹配,将匹配度最高的节点视为主谐波源所在位置。
边界节点谐波电流的波动特征与不可观区域内部主谐波源的波动特征存在相似特征,基于该相似特征的匹配,从边界节点的谐波电流获取回归系数与区域内部节点对应的等效参数,反推不可观区域内部谐波源的位置。
1、建立谐波源定位的等值模型
化简电力系统内的不可观节点,保留相关可观节点,将局部不可观系统转化为等效可观系统;转化过程如图2所示。
由此得到所述局部不可观系统化简后的表达式,记为:
其中,为等效后的节点导纳矩阵,/>为保留节点的电压,/>为原注入电流在边界节点上产生的等效电流。
进一步得到谐波源定位的等值模型如下:
式中,h表示谐波频次,为谐波量测量,含节点谐波电压与支路谐波电流;为系统等效后的谐波量测矩阵;/>为谐波电流;/>为谐波量测的噪声;/>表示复数域,上标m表示量测数,上标n表示状态量数。
求解上述模型,当不可观区域内存在谐波源时,将在边界节点产生等效谐波电流。
2、求解谐波电流动态信息
根据谐波源定位等值模型中连续多个时刻的量测量,采用卡尔曼滤波算法对谐波的状态量进行迭代更新计算。
在预测阶段,谐波电流的预测值和先验误差通过下式获取:
其中,上标-表示先验估计,上标+表示后验估计,下标k表示k时刻,下标k+1表示k+1时刻,即下一时刻,上标^表示估计值;A为状态转移矩阵;为k+1时刻的先验误差矩阵;为k时刻的后验误差矩阵;Q为系统噪声;/>为k时刻的后验状态估计值;/>表示k+1时刻的先验状态估计值,指谐波电流的预测值。
在校正阶段,谐波电流的最优估计值和后验误差通过下式获取:
其中,上标-1为矩阵求逆;R为测量噪声;H为量测矩阵;Kk+1为k+1时刻的增益矩阵;zk+1为k+1时刻的观测值;G为单位矩阵;为k+1时刻的后验误差矩阵;/>为k+1时刻的后验状态估计值,指校正后的谐波电流最优估计值。
谐波电流的动态信息围绕式(3)至(7)进行迭代更新;式(3)预测k+1时刻状态量;式(4)则表示真实值与预测状态量之间的误差;式(5)通过计算增益矩阵来校正更新式(6)中k+1时刻最优的状态量;式(7)更新最优状态量与真实值之间的误差,以执行下一次迭代。
3、基于相似特征匹配定位谐波源
定义不可观区域内存在主谐波源,即其余谐波源基本维持不变,此时边界节点的等效谐波电流受主谐波源影响,边界节点谐波电流的波动特征同时取决于该主谐波源的波动特征,与该主谐波源的波动特征具有一定的相似性,进而可知边界节点谐波电流之间存在一定的联系,因此可进一步描述边界节点之间的等效谐波电流一般线性关系如下:
式中,为边界节点l的等效谐波电流;/>为边界节点r的等效谐波电流变量;为回归系数,对应主谐波源节点对应的等效参数;/>为常量。
对边界节点的谐波电流进行窗口划分,设定该划分的窗宽为n,窗口数为m,对每个窗口进行皮尔逊相关性计算,并选择相关性大于阈值α的谐波电流窗口,此时为主谐波源作用下的窗口。
对窗口内谐波电流进行线性回归,并辨识其回归系数,回归过程如下所示:
Ibrh,l,real=kbrh,r,real,1Ibrh,r,real+kbrh,r,imag,1Ibrh,r,imag+b1 (9)
Ibrh,l,imag=kbrh,r,imag,2Ibrh,r,real+kbrh,r,real,2Ibrh,r,imag+b2 (10)
式中,Ibrh,r,real为边界节点r的等效谐波电流实部,Ibrh,r,imag为边界节点r的等效谐波电流虚部,Ibrh,l,real为边界节点l的等效谐波电流实部,Ibrh,l,imag为边界节点l的等效谐波电流虚部,b为常量。
从上式总共得到四个变量系数,即kbrh,r,real,1、kbrh,r,imag,1、kbrh,r,real,2、kbrh,r,imag,2;将上述四个变量系数与区域内节点对应的等效参数进行匹配对比,定义匹配度评估指标gmatch如下所示:
式中,ktrue,real、ktrue,imag即式(8)节点对应的等效参数实部与虚部。
gmatch越大,表明回归系数越接近节点对应的等效参数,主谐波源处于该节点的可能性越大;选取匹配度最高的节点,该节点即为主谐波源所处位置。
当相关性大于阈值α的谐波电流窗口不存在时,判定该局部不可观区域不存在主谐波源,只进行谐波源区域定位。
以上所述,仅是本发明的较佳实施例而已,并非是对本发明作其它形式的限制,任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例。但是凡是未脱离本发明技术方案内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型,仍属于本发明技术方案的保护范围。
Claims (3)
1.一种采用相似特征匹配的局部不可观系统谐波源定位方法,其特征在于,包括以下步骤:
基于多端口网络等效原理,采用谐波电流移置法,将局部不可观系统转化为等效可观系统;
建立谐波源定位的等值模型,通过卡尔曼滤波算法求解谐波电流动态信息,以此确定谐波源疑似存在区域;
通过相关性计算,筛选不可观区域内主谐波源作用下的谐波电流状态量,并采用线性回归辨识边界节点状态量之间的回归系数,将其与区域内节点对应的等效参数进行匹配,将匹配度最高的节点视为主谐波源所在位置;
基于相似特征匹配定位谐波源,其具体方法如下:
定义不可观区域内存在主谐波源,即其余谐波源保持相对稳定,此时边界节点的等效谐波电流受主谐波源影响,边界节点谐波电流的波动特征同时取决于该主谐波源的波动特征,与该主谐波源的波动特征具有相似性,进而可知边界节点谐波电流之间存在一定的联系,因此进一步描述边界节点之间的等效谐波电流一般线性关系如下:
式中,为边界节点l的等效谐波电流;/>为边界节点r的等效谐波电流变量;/>为回归系数,对应主谐波源节点对应的等效参数;/>为常量;
对边界节点的谐波电流进行窗口划分,设定该划分的窗宽为n,窗口数为m,对每个窗口进行皮尔逊相关性计算,并选择相关性大于阈值α的谐波电流窗口,此时为主谐波源作用下的窗口;
对窗口内谐波电流进行线性回归,并辨识其回归系数,回归过程如下所示:
Ibrh,l,real=kbrh,r,real,1Ibrh,r,real+kbrh,r,imag,1Ibrh,r,imag+b1 (9)
Ibrh,l,imag=kbrh,r,imag,2Ibrh,r,real+kbrh,r,real,2Ibrh,r,imag+b2 (10)
式中,Ibrh,r,real为边界节点r的等效谐波电流实部,Ibrh,r,imag为边界节点r的等效谐波电流虚部,Ibrh,l,real为边界节点l的等效谐波电流实部,Ibrh,l,imag为边界节点l的等效谐波电流虚部,b为常量;
从上式总共得到四个变量系数,即kbrh,r,real,1、kbrh,r,imag,1、kbrh,r,real,2、kbrh,r,imag,2;将上述四个变量系数与区域内节点对应的等效参数进行匹配对比,定义匹配度评估指标gmatch如下所示:
式中,ktrue,real、ktrue,imag即式(8)节点对应的等效参数实部与虚部;
gmatch越大,表明回归系数越接近节点对应的等效参数,主谐波源处于该节点的可能性越大;选取匹配度最高的节点,该节点即为主谐波源所处位置;
当相关性大于阈值α的谐波电流窗口不存在时,判定该局部不可观区域不存在主谐波源,只进行谐波源区域定位。
2.根据权利要求1所述的一种采用相似特征匹配的局部不可观系统谐波源定位方法,其特征在于,建立谐波源定位的等值模型的具体方法如下:
化简电力系统内的不可观节点,保留相关可观节点,将局部不可观系统转化为等效可观系统;
由此得到所述局部不可观系统化简后的表达式,记为:
其中,为等效后的节点导纳矩阵,/>为保留节点的电压,/>为原注入电流在边界节点上产生的等效电流;
进一步得到谐波源定位的等值模型如下:
式中,h表示谐波频次,为谐波量测量,含节点谐波电压与支路谐波电流;为系统等效后的谐波量测矩阵;/>为谐波电流;/>为谐波量测的噪声;/>表示复数域,上标m表示量测量数,上标n表示状态量数;
求解上述模型,当不可观区域内存在谐波源时,将在边界节点产生等效谐波电流。
3.根据权利要求1所述的一种采用相似特征匹配的局部不可观系统谐波源定位方法,其特征在于,求解谐波电流动态信息的具体方法如下:
根据谐波源定位等值模型中连续多个时刻的量测量,采用卡尔曼滤波算法对谐波的状态量进行迭代更新计算;
在预测阶段,谐波电流的预测值和先验误差通过下式获取:
其中,上标-表示先验估计,上标+表示后验估计,下标k表示k时刻,下标k+1表示k+1时刻,即下一时刻,上标^表示估计值;A为状态转移矩阵;为k+1时刻的先验误差矩阵;/>为k时刻的后验误差矩阵;Q为系统噪声;/>为k时刻的后验状态估计值;/>表示k+1时刻的先验状态估计值,指谐波电流的预测值;
在校正阶段,谐波电流的最优估计值和后验误差通过下式获取:
其中,上标-1为矩阵求逆;R为测量噪声;H为量测矩阵;Kk+1为k+1时刻的增益矩阵;zk+1为k+1时刻的观测值;G为单位矩阵;为k+1时刻的后验误差矩阵;/>为k+1时刻的后验状态估计值,指校正后的谐波电流最优估计值;
谐波电流的动态信息围绕式(3)至(7)进行迭代更新;式(3)预测k+1时刻状态量;式(4)则表示真实值与预测状态量之间的误差;式(5)通过计算增益矩阵来校正更新式(6)中k+1时刻最优的状态量;式(7)更新最优状态量与真实值之间的误差,以执行下一次迭代。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111554443.1A CN114236235B (zh) | 2021-12-17 | 2021-12-17 | 一种采用相似特征匹配的局部不可观系统谐波源定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111554443.1A CN114236235B (zh) | 2021-12-17 | 2021-12-17 | 一种采用相似特征匹配的局部不可观系统谐波源定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114236235A CN114236235A (zh) | 2022-03-25 |
CN114236235B true CN114236235B (zh) | 2024-02-06 |
Family
ID=80758460
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111554443.1A Active CN114236235B (zh) | 2021-12-17 | 2021-12-17 | 一种采用相似特征匹配的局部不可观系统谐波源定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114236235B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0826972A2 (en) * | 1996-08-28 | 1998-03-04 | Matsushita Electric Industrial Co., Ltd. | Harmonic measurement system |
CN103675446A (zh) * | 2013-12-23 | 2014-03-26 | 重庆大学 | 一种实时精确的谐波源定位系统和方法 |
CN103969552A (zh) * | 2014-05-22 | 2014-08-06 | 河海大学 | 一种分布式发电系统的谐波源定位分析方法 |
KR20150028450A (ko) * | 2013-09-06 | 2015-03-16 | 한국전력공사 | 고조파 발생원 추정 장치 및 방법 |
CN104793082A (zh) * | 2015-04-23 | 2015-07-22 | 江苏中凌高科技股份有限公司 | 基于谐波相关分析的电力系统谐波源识别装置 |
CN105911334A (zh) * | 2016-04-19 | 2016-08-31 | 国网上海市电力公司 | 一种谐波阻抗未知条件下的谐波电流估计方法 |
CN109581103A (zh) * | 2018-11-21 | 2019-04-05 | 上海交通大学 | 基于广域监测的电网谐波源定位方法 |
CN109839538A (zh) * | 2019-03-29 | 2019-06-04 | 云南电网有限责任公司电力科学研究院 | 一种基于相关性分析的谐波源辨识方法和系统 |
CN110308366A (zh) * | 2019-06-26 | 2019-10-08 | 国网山西省电力公司电力科学研究院 | 一种基于正交匹配追踪算法的谐波源定位方法 |
CN113051826A (zh) * | 2021-03-29 | 2021-06-29 | 福州大学 | 基于高斯过程回归的谐波源通用不确定性建模方法 |
-
2021
- 2021-12-17 CN CN202111554443.1A patent/CN114236235B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0826972A2 (en) * | 1996-08-28 | 1998-03-04 | Matsushita Electric Industrial Co., Ltd. | Harmonic measurement system |
KR20150028450A (ko) * | 2013-09-06 | 2015-03-16 | 한국전력공사 | 고조파 발생원 추정 장치 및 방법 |
CN103675446A (zh) * | 2013-12-23 | 2014-03-26 | 重庆大学 | 一种实时精确的谐波源定位系统和方法 |
CN103969552A (zh) * | 2014-05-22 | 2014-08-06 | 河海大学 | 一种分布式发电系统的谐波源定位分析方法 |
CN104793082A (zh) * | 2015-04-23 | 2015-07-22 | 江苏中凌高科技股份有限公司 | 基于谐波相关分析的电力系统谐波源识别装置 |
CN105911334A (zh) * | 2016-04-19 | 2016-08-31 | 国网上海市电力公司 | 一种谐波阻抗未知条件下的谐波电流估计方法 |
CN109581103A (zh) * | 2018-11-21 | 2019-04-05 | 上海交通大学 | 基于广域监测的电网谐波源定位方法 |
CN109839538A (zh) * | 2019-03-29 | 2019-06-04 | 云南电网有限责任公司电力科学研究院 | 一种基于相关性分析的谐波源辨识方法和系统 |
CN110308366A (zh) * | 2019-06-26 | 2019-10-08 | 国网山西省电力公司电力科学研究院 | 一种基于正交匹配追踪算法的谐波源定位方法 |
CN113051826A (zh) * | 2021-03-29 | 2021-06-29 | 福州大学 | 基于高斯过程回归的谐波源通用不确定性建模方法 |
Non-Patent Citations (1)
Title |
---|
电能质量谐波监测与评估综述;严静等;《电气技术》 * |
Also Published As
Publication number | Publication date |
---|---|
CN114236235A (zh) | 2022-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jayawardena et al. | Noise reduction and prediction of hydrometeorological time series: dynamical systems approach vs. stochastic approach | |
CN107092744B (zh) | 基于emd-svr的地表沉降量预测方法 | |
CN107658881A (zh) | 基于戴维南等值方法的电压稳定临界点判断方法 | |
CN111046327B (zh) | 适用于低频振荡与次同步振荡辨识的Prony分析方法 | |
CN109946518B (zh) | 基于贝叶斯方法的电力谐波信号分析方法与分析设备 | |
CN108647807B (zh) | 河流流量的预测方法 | |
CN102288843A (zh) | 一种电能质量扰动信号检测方法 | |
CN110119588A (zh) | 基于扩展卡尔曼滤波状态估计值的在线优化设计方法 | |
CN109754013A (zh) | 一种基于无迹卡尔曼滤波的电力系统混合量测融合方法 | |
CN111541255B (zh) | 基于动力学系统的低频振荡模态识别方法及系统 | |
Song et al. | Iterative interpolation for parameter estimation of LFM signal based on fractional Fourier transform | |
CN115201563A (zh) | 一种基于联合熵的多谐波源定位方法及系统 | |
CN114236235B (zh) | 一种采用相似特征匹配的局部不可观系统谐波源定位方法 | |
CN112697215B (zh) | 一种用于超声波水表数据滤波的卡尔曼滤波参数调试方法 | |
CN107204616B (zh) | 基于自适应稀疏伪谱法的电力系统随机状态估计方法 | |
CN113514743A (zh) | 一种基于多维特征的gis局部放电模式识别系统构建方法 | |
CN116644583A (zh) | 一种基于vmd和弗雷歇距离相结合的油气管道泄漏信号预处理方法及系统 | |
CN116565856A (zh) | 一种考虑未知拓扑变化的配电网状态估计方法 | |
CN115293090A (zh) | 基于重构数据处理的多谐波源责任量化方法 | |
CN114236236B (zh) | 一种基于区间动态状态估计的谐波源定位方法 | |
CN113705335B (zh) | 一种低频稀疏线谱信号时频超分辨率提取方法 | |
CN107576866A (zh) | 一种基于光滑逼近稀疏重构与内点法的多谐波源辨识方法 | |
CN104881077A (zh) | 一种光伏系统中最大功率点跟踪控制方法 | |
CN117630487B (zh) | 基于lof筛选-启发式分割算法的系统谐波阻抗估计方法 | |
CN113204831B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |