CN109446643B - 一种基于实测数据建立家用电器负荷谐波模型的方法 - Google Patents
一种基于实测数据建立家用电器负荷谐波模型的方法 Download PDFInfo
- Publication number
- CN109446643B CN109446643B CN201811258954.7A CN201811258954A CN109446643B CN 109446643 B CN109446643 B CN 109446643B CN 201811258954 A CN201811258954 A CN 201811258954A CN 109446643 B CN109446643 B CN 109446643B
- Authority
- CN
- China
- Prior art keywords
- matrix
- household appliance
- load
- harmonic
- voltage
- 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 57
- 239000011159 matrix material Substances 0.000 claims abstract description 70
- 239000013598 vector Substances 0.000 claims description 21
- 230000001419 dependent effect Effects 0.000 claims description 16
- 238000005259 measurement Methods 0.000 claims description 10
- 230000003190 augmentative effect Effects 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 241000282465 Canis Species 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000011156 evaluation Methods 0.000 abstract description 2
- 238000013488 ordinary least square regression Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 230000005611 electricity Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000003416 augmentation Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于实测数据建立家用电器负荷谐波模型的方法,包括:获取家用电器负荷的实测电压和电流数据,对测量得到的数据进行傅里叶分析得到电压和电流的各次谐波的相量值;建立家用电器负荷的交叉频率导纳矩阵谐波模型;将得到的电压和电流的各次谐波的相量值作为输入,基于广义线性复偏最小二乘法计算交叉频率导纳矩阵模型参数,得到家用电器交叉频率导纳谐波模型。本发明有益效果:该方法可在不需知道家用负荷内部电路结构与参数的基础上,具有准确、快速且方法简单实用,对于保证准确分析与评估配电网谐波产生情况有重大意义。
Description
技术领域
本发明属于负荷谐波建模领域,尤其涉及一种基于实测数据建立家用电器负荷谐波模型的方法。
背景技术
当今电力***中,居民用电占电能总消耗的30%以上,而且随着智能配电网的发展,居民用电比例还在持续增长。然而,由于家用负荷中存在着大量的电力电子装置,使得大多数家用电器的谐波电流畸变非常严重,有些电流谐波总畸变率甚至超过了100%,实测数据表明,某些居民配电馈线上的谐波电流总需求畸变率可达12%,使得居民负荷成为非常重要的一类谐波源。由此而产生的谐波会导致居民配电网中谐波损耗增加,变压器和中性线过载,电子设备出现故障的概率增加等一系列问题。因此为了准确的评估和分析家用负荷的大量接入对民用配电网的影响,并相应地提出适当的治理措施来提高电能质量,非常有必要对家用负荷建立准确的谐波模型。
目前关于家用电器负荷的谐波建模方法可以分为基于等值电路分析的方法和基于实测数据的建模方法。基于等值电路分析的方法需要已知家用电器负荷的等值电路以及电路参数的具体数值。然而实际上,很多家用电器负荷的供电电路结构复杂,难以获得其简化等值电路,而且电路的参数不恒定或难以估算。基于实测数据的建模方法仅需要对负荷进行实际测量获取其电压和电流数据即可。以往分析中常采用恒流源模型计算家用电器负荷的谐波电流,但恒流源模型简化了供电电压对负荷谐波电流的影响,使得所计算的谐波电流不跟随供电电压变化,因此对供电电压敏感的谐波负荷会产生较大误差。为此提出利用交叉频率导纳矩阵模型对家用电器负荷进行谐波建模,然而有利用普通最小二乘法求解交叉频率导纳矩阵模型参数过程中,在计算过程中会出现矩阵病态的问题,导致求解的模型参数不准确,从而影响建模精度。
发明内容
为了解决上述问题,本发明提出了一种基于实测数据建立家用电器负荷谐波模型的方法,能够准确的评估家用负荷大量接入对民用配电网所产生的影响,有效地提升配电网的电能质量。
为了实现上述目的,本发明采用如下技术方案:
在一种或多种实施方式中公开的一种基于实测数据建立家用电器负荷谐波模型的方法,包括:
(1)获取家用电器负荷的实测电压和电流数据,对测量得到的数据进行傅里叶分析得到电压和电流的各次谐波的相量值;
(2)建立家用电器负荷的交叉频率导纳矩阵谐波模型;
(3)将步骤(1)中得到的电压和电流的各次谐波的相量值作为输入,基于广义线性复偏最小二乘法计算交叉频率导纳矩阵模型参数,得到家用电器交叉频率导纳谐波模型。
进一步地,所述获取家用电器负荷的实测电压和电流数据的方法,具体为:
建立实验测量平台,家用负荷的供电电压为实际电网电压,通过电压钳和电流钳测量其供电电压uM和实际电流iM,通过数据采集装置获取实际测量的供电电压和实际电流数据。
进一步地,建立家用电器负荷的交叉频率导纳矩阵谐波模型,具体为:
进一步地,将家用电器负荷的交叉频率导纳矩阵谐波模型进行扩展;
将扩展后的家用电器负荷的交叉频率导纳矩阵谐波模型变换为紧凑形式,所述紧凑形式为因变量矩阵I等于自变量矩阵U与交叉频率导纳矩阵Y的乘积;
对自变量矩阵U和因变量矩阵I分别进行增广,通过迭代方法提取最大相关的变量,基于广义线性复偏最小二乘法提取主成分,得到自变量矩阵U和因变量矩阵I的估测值;
计算得到交叉频率导纳矩阵。
进一步地,将家用电器负荷的交叉频率导纳矩阵谐波模型进行扩展,具体为:
进一步地,对自变量矩阵U和因变量矩阵I分别进行增广,分别为U=[U,U*],I=[I,I*],其中上标(·)*为取共轭运算;
通过迭代方法提取最大相关的变量,具体为:
其中,t m是U和I的第m个主成分相量,p m和q m分别是U和I的第m个主成分所对应的载荷向量,T是对应的主成分矩阵,P和Q是对应的载荷矩阵。U M和I M是U和I的残差矩阵,M表示最大主成分数,通过交叉检验法确定。
进一步地,基于广义线性复偏最小二乘法提取主成分,具体过程如下:
对自变量矩阵U和因变量矩阵I分别进行增广,可以得到初始输入矩阵分别为U 1=[U,U*],I 1=[I,I*],U1=U,I1=I;其中上标(·)*为取共轭运算;
假设i=1,…,M;
重复上述步骤,提取M个主成分向量,其中M的大小通过交叉检验法确定。
进一步地,计算得到交叉频率导纳矩阵,具体为:
Y=(M+N×e-j2θ)。
与现有技术相比,本发明的有益效果是:
针对家用电器负荷,所提出的利用广义线性复偏最小二乘法计算其交叉频率导纳矩阵模型参数,该方法可在不需知道家用负荷内部电路结构与参数的基础上,具有准确、快速且方法简单实用,对于保证准确分析与评估配电网谐波产生情况有重大意义。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。
图1为一种基于实测数据建立家用电器负荷谐波模型的方法流程图;
图2为用于实验测量的实验平台电路图;
图3(a)为紧凑型荧光灯基于本发明方法、基于普通最小二乘法计算的和实际电流波形的电流波形图以及频谱对比图;
图3(b)为电动自行车充电器基于本发明方法、基于普通最小二乘法计算的和实际电流波形的电流波形图以及频谱对比图;
图3(c)为台式机电脑基于本发明方法、基于普通最小二乘法计算的和实际电流波形的电流波形图以及频谱对比图。
具体实施方式
应该指出,以下详细说明都是例示性的,旨在对本申请提供进一步的说明。除非另有指明,本发明使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
实施例一
在一个或多个实施方式中公开的一种基于实测数据建立家用电器负荷谐波模型的方法,具体流程如下:
1、建立实验测量平台,获取家用电器负荷的实测电压和电流数据。
如图2所示,建立实验测量平台,为了考虑实际电网中的电压对家用电器负荷的影响,家用负荷的供电电压为实际电网电压,通过电压钳和电流钳测量其供电电压uM和实际电流iM,结合数据采集板和数据采集软件LABVIEW获取其实际电压和电流数据,之后对所采集的数据进行傅里叶分析,获取其电压相量值和电流的相量值。
2、建立家用电器负荷的交叉频率导纳矩阵谐波模型。
针对家用电器负荷建立其交叉频率导纳矩阵谐波模型,具体表达式为:
为了利用测量所得到的电压和电流相量值计算交叉频率导纳矩阵谐波模型的参数,需要将式(1)扩展为:
式(2)的紧凑形式为:
I=YU(3)
3、基于广义线性复偏最小二乘法计算其交叉频率导纳矩阵谐波模型参数。
对步骤(3)中的式(3)所得到的自变量矩阵U和因变量矩阵I,分别进行增广,可以得到其增广矩阵,分别为U=[U,U*],I=[I,I*],其中上标(·)*为取共轭运算。通过迭代方法提取最大相关的变量,如式(4)所示:
其中t m是U和I的第m个主成分相量,p m和q m分别是U和I的第m个主成分所对应的载荷向量,T是对应的主成分矩阵,P和Q是对应的载荷矩阵。U M和I M是U和I的残差矩阵,其中M表示最大主成分数,通过交叉检验法确定。
基于广义线性复偏最小二乘法提取主成分,首先对自变量矩阵U和因变量矩阵I分别进行增广,可以得到初始输入矩阵分别为U 1=[U,U*],I 1=[I,I*],U1=U,I1=I,其中上标(·)*为取共轭运算;当i=1时,求取矩阵最大特征值所对应的特征向量w1,其中Eigmax{·}表示取最大特征值的特征向量运算;然后再求第一个主成分向量和其增广向量之后再求取第一个主成分向量对应的载荷向量最后通过主成分向量和载荷向量可以得到i=2时的输入矩阵为再重复之前的循环步骤,可以提取M个主成分向量,其中M的大小可以通过交叉检验法确定。
基于广义线性复偏最小二乘法提取主成分步骤如下:
1)初始化:U 1=[U,U*],I 1=[I,I*],U1=U,I1=I
2)循环:假设i=1,…,M
7)存储wi,ti,p i和q i
8)结束循环
矩阵M和N可以计算:
因此利用广义线性复偏最小二乘法计算出的交叉频率导纳矩阵为:
Y=(M+N×e-j2θ) (7)
为了验证本发明所提计算方法,采用了普通最小二乘法作为对比的计算方法,各家用负荷的供电电压均是实际电网电压,对如下表所示的家用负荷进行实验验证:
下表所示的为以上三种家用负荷的模型重构电流和实际电流的各次谐波电流幅值和相角之间的误差:
以上三种家用负荷通过本发明所提的方法计算、实测结果与普通最小二乘法比较的结果如图3(a)-(c)所示,由图中结果可以知道基于本发明所提出的方法用于家用负荷谐波建模与实际结果的电流波形图基本吻合,且误差基本保持在10%以内,而基于普通最小二乘法进行谐波建模所得到的电流波形与实际电流波形偏差较大,误差基本在30%以上。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (5)
1.一种基于实测数据建立家用电器负荷谐波模型的方法,其特征在于,包括:
(1)获取家用电器负荷的实测电压和电流数据,对测量得到的数据进行傅里叶分析得到电压和电流的各次谐波的相量值;
(2)建立家用电器负荷的交叉频率导纳矩阵谐波模型,
(3)将步骤(1)中得到的电压和电流的各次谐波的相量值作为输入,基于广义线性复偏最小二乘法计算交叉频率导纳矩阵模型参数,得到家用电器交叉频率导纳谐波模型;
将家用电器负荷的交叉频率导纳矩阵谐波模型进行扩展;
将扩展后的家用电器负荷的交叉频率导纳矩阵谐波模型变换为紧凑形式,所述紧凑形式为因变量矩阵I等于自变量矩阵U与交叉频率导纳矩阵Y的乘积;
对自变量矩阵U和因变量矩阵I分别进行增广,通过迭代方法提取最大相关的变量,基于广义线性复偏最小二乘法提取主成分,得到自变量矩阵U和因变量矩阵I的估测值;
计算得到交叉频率导纳矩阵;
计算得到交叉频率导纳矩阵具体为:
Y=(M+N×e-j2θ)。
2.如权利要求1所述的一种基于实测数据建立家用电器负荷谐波模型的方法,其特征在于,所述获取家用电器负荷的实测电压和电流数据的方法,具体为:
建立实验测量平台,家用负荷的供电电压为实际电网电压,通过电压钳和电流钳测量其供电电压uM和实际电流iM,通过数据采集装置获取实际测量的供电电压和实际电流数据。
5.如权利要求4所述的一种基于实测数据建立家用电器负荷谐波模型的方法,其特征在于,基于广义线性复偏最小二乘法提取主成分,具体过程如下:
对自变量矩阵U和因变量矩阵I分别进行增广,可以得到初始输入矩阵分别为U 1=[U,U*],I 1=[I,I*],U1=U,I1=I;其中上标(·)*为取共轭运算;
假设i=1,…,M;
重复上述步骤,提取M个主成分向量,其中M的大小通过交叉检验法确定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811258954.7A CN109446643B (zh) | 2018-10-26 | 2018-10-26 | 一种基于实测数据建立家用电器负荷谐波模型的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811258954.7A CN109446643B (zh) | 2018-10-26 | 2018-10-26 | 一种基于实测数据建立家用电器负荷谐波模型的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109446643A CN109446643A (zh) | 2019-03-08 |
CN109446643B true CN109446643B (zh) | 2020-10-30 |
Family
ID=65547563
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811258954.7A Active CN109446643B (zh) | 2018-10-26 | 2018-10-26 | 一种基于实测数据建立家用电器负荷谐波模型的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109446643B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110502804B (zh) * | 2019-07-29 | 2020-04-21 | 山东大学 | 一种电动汽车充电过程时变谐波电流评估方法及*** |
CN110896218B (zh) * | 2019-11-08 | 2021-09-03 | 山东大学 | 一种用于建立集合性居民负荷的谐波建模方法及*** |
CN114935688B (zh) * | 2022-07-25 | 2022-10-14 | 山东大学 | 基于功率分段的电弧炉供电***谐波评估方法与*** |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103258103A (zh) * | 2013-05-29 | 2013-08-21 | 国家电网公司 | 基于偏最小二乘回归的戴维南等值参数辨识方法 |
CN103532205A (zh) * | 2013-10-31 | 2014-01-22 | 重庆大学 | 一种电动汽车三相充电机谐波模型的建模方法 |
CN103544378A (zh) * | 2013-09-28 | 2014-01-29 | 南方电网科学研究院有限责任公司 | 一种直流输电用交流***谐波阻抗计算方法 |
CN104123451A (zh) * | 2014-07-16 | 2014-10-29 | 河海大学常州校区 | 基于偏最小二乘回归的疏浚作业产量预测模型建立方法 |
CN104198575A (zh) * | 2014-09-23 | 2014-12-10 | 华中科技大学 | 具有溯源功能的电感耦合等离子质谱信号分析方法 |
CN104569528A (zh) * | 2014-12-25 | 2015-04-29 | 国网重庆市电力公司永川供电分公司 | 一种计及负荷的谐波电能损耗的电能计费方法 |
CN105911334A (zh) * | 2016-04-19 | 2016-08-31 | 国网上海市电力公司 | 一种谐波阻抗未知条件下的谐波电流估计方法 |
CN109116183A (zh) * | 2018-06-28 | 2019-01-01 | 全球能源互联网研究院有限公司 | 谐波模型参数辨识方法、装置、存储介质及电子设备 |
-
2018
- 2018-10-26 CN CN201811258954.7A patent/CN109446643B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103258103A (zh) * | 2013-05-29 | 2013-08-21 | 国家电网公司 | 基于偏最小二乘回归的戴维南等值参数辨识方法 |
CN103544378A (zh) * | 2013-09-28 | 2014-01-29 | 南方电网科学研究院有限责任公司 | 一种直流输电用交流***谐波阻抗计算方法 |
CN103532205A (zh) * | 2013-10-31 | 2014-01-22 | 重庆大学 | 一种电动汽车三相充电机谐波模型的建模方法 |
CN104123451A (zh) * | 2014-07-16 | 2014-10-29 | 河海大学常州校区 | 基于偏最小二乘回归的疏浚作业产量预测模型建立方法 |
CN104198575A (zh) * | 2014-09-23 | 2014-12-10 | 华中科技大学 | 具有溯源功能的电感耦合等离子质谱信号分析方法 |
CN104569528A (zh) * | 2014-12-25 | 2015-04-29 | 国网重庆市电力公司永川供电分公司 | 一种计及负荷的谐波电能损耗的电能计费方法 |
CN105911334A (zh) * | 2016-04-19 | 2016-08-31 | 国网上海市电力公司 | 一种谐波阻抗未知条件下的谐波电流估计方法 |
CN109116183A (zh) * | 2018-06-28 | 2019-01-01 | 全球能源互联网研究院有限公司 | 谐波模型参数辨识方法、装置、存储介质及电子设备 |
Non-Patent Citations (2)
Title |
---|
主电路串联电抗器的无源滤波器研究;王华军;《中国优秀硕士学位论文全文数据库.信息科技辑》;20100515;第2010年卷(第5期);摘要、正文第5章 * |
基于复数域偏最小二乘法与等值法的多谐波源责任划分;王辉 等;《电力***自动化》;20170225;第41卷(第4期);摘要、正文第1-4节 * |
Also Published As
Publication number | Publication date |
---|---|
CN109446643A (zh) | 2019-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Roinila et al. | Online grid-impedance measurement using ternary-sequence injection | |
CN109446643B (zh) | 一种基于实测数据建立家用电器负荷谐波模型的方法 | |
Jo et al. | On improving distortion power quality index in distributed power grids | |
CN107121609A (zh) | 一种基于prbs扰动注入的电网阻抗在线辨识方法及装置 | |
CN109030941A (zh) | Hanning自乘卷积窗FFT三谱线插值谐波分析方法 | |
Nascimento et al. | Harmonic distortion monitoring for nonlinear loads using neural-network-method | |
CN101807795A (zh) | 一种电能计量仿真***构成方法及其装置 | |
Xie et al. | Data-driven dynamic harmonic model for modern household appliances | |
CN108318852B (zh) | 一种智能电能表方波影响试验方法 | |
Caicedo et al. | Frequency domain modeling of nonlinear loads, considering harmonic interaction | |
Wang et al. | Data-driven probabilistic small signal stability analysis for grid-connected PV systems | |
Chen et al. | Robust data-driven linearization for distribution three-phase power flow | |
CN114301055B (zh) | 基于宽频量测的电力***间谐波潮流获取方法与*** | |
Zolfaghari et al. | Evaluation of windowed ESPRIT virtual instrument for estimating Power Quality Indices | |
Gerber et al. | Parameter estimation of a ferro-resonance damping circuit using pseudo-random impulse sequence perturbations | |
CN106356867B (zh) | 一种节点无功电压灵敏度的确定方法 | |
CN110221168B (zh) | 一种主导谐波源定位和谐波污染传播路径追踪方法 | |
CN106953331A (zh) | 一种考虑不确定性和三相不平衡的电力***谐波潮流算法 | |
CN109840717A (zh) | 一种含分布式电源的配电网线性潮流计算方法 | |
CN105429143B (zh) | 一种针对电力***光伏专供线的谐波质量监测点选址方法 | |
Pigazo et al. | Accurate and computationally efficient implementation of the IEEE 1459-2000 standard in three-phase three-wire power systems | |
Uhl et al. | Non-linear behavioral X-Parameters model of single-phase rectifier in the frequency domain | |
CN106503290B (zh) | 一种基于pls-svm的通用单相谐波负荷建模方法 | |
CN110896218B (zh) | 一种用于建立集合性居民负荷的谐波建模方法及*** | |
CN107561380A (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 |