CN109829652B - 一种长时间尺度动态谐波责任划分方法 - Google Patents
一种长时间尺度动态谐波责任划分方法 Download PDFInfo
- Publication number
- CN109829652B CN109829652B CN201910110713.6A CN201910110713A CN109829652B CN 109829652 B CN109829652 B CN 109829652B CN 201910110713 A CN201910110713 A CN 201910110713A CN 109829652 B CN109829652 B CN 109829652B
- Authority
- CN
- China
- Prior art keywords
- responsibility
- harmonic
- dynamic
- index
- long
- 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 32
- 238000012544 monitoring process Methods 0.000 claims abstract description 21
- 238000010606 normalization Methods 0.000 claims abstract description 7
- 238000012545 processing Methods 0.000 claims abstract description 4
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000005611 electricity Effects 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 230000007774 longterm Effects 0.000 abstract description 10
- 238000011156 evaluation Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000013475 authorization Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- 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
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种长时间尺度动态谐波责任划分方法,包括以下步骤:步骤S1:采集需划分责任时段内谐波电压监测数据序列和各个用户在用电信息采集系统的平均有功功率计量数据;步骤S2:根据得到的谐波电压监测数据序列和平均有功功率计量数据,计算所有用户的平均有功功率与谐波电压数据的动态相关系数;步骤S3:根据得到的动态相关系数和评价有功功率值,构造动态谐波责任指标,并进行归一化处理;步骤S4:根据归一化后的动态谐波责任指标,构造长期谐波责任指标;步骤S5:根据得到的长期谐波责任指标,进行需划分责任时段内的谐波责任划分。本发明可用于解决用户难以接受短时谐波责任划分结果的问题。
Description
技术领域
本发明涉及一种长时间尺度动态谐波责任划分方法。
背景技术
电力电子负荷越来越多地接入电网,使得电网中的谐波情况更加复杂,影响负荷正常工作,严重时还会导致电力设备损坏或引发停电事故。为了有效的实施“奖惩性方案”,需要制定合理的谐波责任指标,以定量衡量每个用户发射谐波对电网的影响。实际电力系统中一条母线往往连接多个发出谐波的终端用户,这些用户间的谐波相互影响,如何根据已有数据制定合理的谐波责任指标,是有效治理谐波,并在未来执行“优质优价”方案,保证电网和用户安全高效运行的前提;同时,谐波是一种长期的稳态电能质量扰动,因此,在实际工程中应考虑长时间尺度的责任划分,利用现有监测系统数据,得出用户更易认可的结果。目前现有的谐波责任划分指标尚无法实现长时间尺度的动态责任划分。
目前多谐波源谐波责任指标可分为两种:谐波电压责任指标、谐波电流责任指标。谐波电压责任指标与谐波电流责任指标的基本原理相同,都是计算用户等效谐波源单独作用时,作用于公共连接点的电压或电流向量,在公共连接点电压或电流向量方向的投影,占公共连接点电压或电流向量的百分比作为指标。
目前此问题存在以下技术缺陷:
1、针对短时段的谐波责任划分难以被用户所接受。现有指标大多基于局部线性化的原理求解,这种方法要求分析时段较短。但是用户发射谐波是一个长期动态行为,用短时分析确定用户的谐波责任,难以被用户所接受。
2、针对馈线的谐波责任划分难以具体到每个用户。目前电能质量监测系统监测点有限,无法做到对每个专变用户均进行监测,仅靠安装在变电站母线上的电能质量监测终端提供的监测数据无法对每个用户的谐波责任指标进行准确的计算。
3、针对有效值的谐波责任划分指标难以运用于实际。从现有谐波责任指标使用的数据来看,需要PCC谐波电压与各条馈线谐波电流有效值数据以及相位数据;然而目前国网公司部署的电能质量监测系统中只能提供长时统计值,相位数据的统计值没有物理意义,且电压与电流统计值也并不对应,这导致传统方法在现有监测系统中均无法适用。
因此,实际工程上,迫切需要一种利用实际电能质量监测系统统计数据,融合其它电力系统数据(如用电信息采集系统的用户用电数据),制定一个新的长期谐波责任指标。
发明内容
有鉴于此,本发明的目的在于提供一种长时间尺度动态谐波责任划分方法,可用于解决用户难以接受短时谐波责任划分结果的问题。
为实现上述目的,本发明采用如下技术方案:
一种长时间尺度动态谐波责任划分方法,包括以下步骤:
步骤S1:采集需划分责任时段内谐波电压监测数据序列和各个用户在用电信息采集系统的平均有功功率计量数据;
步骤S2:根据得到的谐波电压监测数据序列和平均有功功率计量数据,计算所有用户的平均有功功率与谐波电压数据的动态相关系数,分别为D1=[d11,d12,…,d1n]、D2=[d21,d22,…,d2n]、……、Dm=[dm1,dm2,…,dmn];
步骤S3:根据得到的动态相关系数和评价有功功率值,构造动态谐波责任指标,并进行归一化处理;
步骤S4:根据归一化后的动态谐波责任指标,构造长期谐波责任指标;
步骤S5:根据得到的长期谐波责任指标,进行需划分责任时段内的谐波责任划分。
进一步的,所述步骤S1具体为:采集需划分责任时段内谐波电压监测数据序列U=[u1,u2,…,un],此时段内连接于该母线的各个用户在用电信息采集系统的平均有功功率计量数据,设共有m个用户,用户的平均有功功率序列为P1=[p11,p12,…,p1n]、P2=[p21,p22,…,p2n]、……、Pm=[pm1,pm2,…,pmn]。
进一步的,所述步骤S3具体为:
步骤S31:构造动态谐波责任指标Ak
步骤S32:对Ak使用归一化方法,使每个用户每一时刻谐波责任在区间[0,1]之内,且同一时刻所有用户谐波责任总和为100%,得到归一化后的动态谐波责任指标Fk。
进一步的,所述责任拟合函数具体为:
g(x)=a+bx+cx2+dx3+ex4 (2)
其中a、b、c、d、e为待定系数。
进一步的,采用网格搜索法确定责任拟合函数中的待定系数,具体为:
1)计算每个用户每个时间节点的动态谐波责任指标Ak,设第k个用户i时刻的动态谐波责任指标为Ak[1,i],可表示为
4)采用网格搜索法求取系数,确定搜索步长,搜索网格中每一组a、b、c、d、e式(5)的取值,当式(5)达到最大值时的a、b、c、d、e即为所求系数。
进一步的,所述采用网格搜索法求取系数具体为:
(1)对a、b、c、d、e、f赋初值,并且确定网格搜索下限为-100、上限为+100、步长为1;
(2)从网格下限开始,计算每个步长a、b、c、d、e的值,并代入式(5)计算f,比较本步f与前一步f的大小,若本步f大于前一步f,则记录本步a、b、c、d、e的值;若本步f小于前一步f,则记录前一步a、b、c、d、e的值,直到网格上限;
(3)返回a、b、c、d、e的值。
进一步的,所述步骤S32具体为:
步骤S321:将谐波责任指标Ak在不改变原本指标的相对大小关系的前提下去除负值,得到Ek[i,1]
式中,Amax[i,1]、Amin[i,1]分别表示所有用户在第i个分析时段的谐波责任指标Ak的最大值、最小值;
步骤S322:根据Ek[i,1],得到归一化后的谐波责任指标:
进一步的,所述步骤S5具体为:
设用户k长时段总谐波责任指标Hk;长时间谐波责任指标是对动态责任指标随时间变化曲线进行积分,再计算各积分结果占所有用户积分总和的百分比,计算方法如式(8)所示:
本发明与现有技术相比具有以下有益效果:
本发明可用于解决用户难以接受短时谐波责任划分结果的问题,不需要专门进行测试或增加安装专用设备,原理简单、使用方便,利用经国家授权检定的电表数据,其结果客观可信,更具工程实用性、具有较强的推广应用价值。
附图说明
图1是本发明方法流程图;
图2是本发明实施例中谐波责任工程问题示意图;
图3是本发明实施例中责任拟合函数示意图。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
请参照图1和图2,本发明提供一种长时间尺度动态谐波责任划分方法,包括以下步骤:
步骤S1:采集需划分责任时段内谐波电压监测数据序列和各个用户在用电信息采集系统的平均有功功率计量数据;
步骤S2:根据得到的谐波电压监测数据序列和平均有功功率计量数据,计算所有用户的平均有功功率与谐波电压数据的动态相关系数,分别为D1=[d11,d12,…,d1n]、D2=[d21,d22,…,d2n]、……、Dm=[dm1,dm2,…,dmn];
步骤S3:根据得到的动态相关系数和评价有功功率值,构造动态谐波责任指标,并进行归一化处理;
步骤S4:根据归一化后的动态谐波责任指标,构造长期谐波责任指标;
步骤S5:根据得到的长期谐波责任指标,进行需划分责任时段内的谐波责任划分。
在本实施例中,所述步骤S1具体为:采集需划分责任时段内谐波电压监测数据序列U=[u1,u2,…,un],此时段内连接于该母线的各个用户在用电信息采集系统的平均有功功率计量数据,设共有m个用户,用户的平均有功功率序列为P1=[p11,p12,…,p1n]、P2=[p21,p22,…,p2n]、……、Pm=[pm1,pm2,…,pmn]。
在本实施例中,所述步骤S3具体为:
步骤S31:构造动态谐波责任指标Ak
步骤S32:对Ak使用归一化方法,使每个用户每一时刻谐波责任在区间[0,1]之内,且同一时刻所有用户谐波责任总和为100%,得到归一化后的动态谐波责任指标Fk。
在本实施例中,所述责任拟合函数具体为:
g(x)=a+bx+cx2+dx3+ex4 (2)
其中a、b、c、d、e为待定系数。
进一步的,采用网格搜索法确定责任拟合函数中的待定系数,具体为:
1)计算每个用户每个时间节点的动态谐波责任指标Ak,设第k个用户i时刻的动态谐波责任指标为Ak[1,i],可表示为
4)采用网格搜索法求取系数,确定搜索步长,搜索网格中每一组a、b、c、d、e式(5)的取值,当式(5)达到最大值时的a、b、c、d、e即为所求系数。
在本实施例中,所述采用网格搜索法求取系数具体为:
(1)对a、b、c、d、e、f赋初值,并且确定网格搜索下限为-100、上限为+100、步长为1;
(2)从网格下限开始,计算每个步长a、b、c、d、e的值,并代入式(5)计算f,比较本步f与前一步f的大小,若本步f大于前一步f,则记录本步a、b、c、d、e的值;若本步f小于前一步f,则记录前一步a、b、c、d、e的值,直到网格上限;
(3)返回a、b、c、d、e的值。
在本实施例中,所述步骤S32具体为:
步骤S321:将谐波责任指标Ak在不改变原本指标的相对大小关系的前提下去除负值,得到Ek[i,1]
式中,Amax[i,1]、Amin[i,1]分别表示所有用户在第i个分析时段的谐波责任指标Ak的最大值、最小值;
步骤S322:根据Ek[i,1],得到归一化后的谐波责任指标:
在本实施例中,所述步骤S5具体为:
设用户k长时段总谐波责任指标Hk;长时间谐波责任指标是对动态责任指标随时间变化曲线进行积分,再计算各积分结果占所有用户积分总和的百分比,计算方法如式(8)所示:
以上所述仅为本发明的较佳实施例,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本发明的涵盖范围。
Claims (5)
1.一种长时间尺度动态谐波责任划分方法,其特征在于,包括以下步骤:
步骤S1:采集需划分责任时段内谐波电压监测数据序列和各个用户在用电信息采集系统的平均有功功率计量数据;
步骤S2:根据谐波电压监测数据序列和平均有功功率计量数据,计算所有用户的平均有功功率与谐波电压监测数据的动态相关系数,分别为D1=[d11,d12,…,d1n]、D2=[d21,d22,…,d2n]、……、Dm=[dm1,dm2,…,dmn];
步骤S3:根据动态相关系数和平均有功功率,构造动态谐波责任指标,并进行归一化处理;
步骤S4:根据归一化后的动态谐波责任指标,构造长时间谐波责任指标;
步骤S5:根据长时间谐波责任指标,进行需划分责任时段内的谐波责任划分;
所述步骤S3具体为:
步骤S31:构造动态谐波责任指标Ak
步骤S32:对Ak使用归一化方法,使每个用户每一时刻谐波责任在区间[0,1]之内,且同一时刻所有用户谐波责任总和为100%,得到归一化后的动态谐波责任指标Fk;
所述责任拟合函数具体为:
g(x)=a+bx+cx2+dx3+ex4 (2)
其中a、b、c、d、e为待定系数;
采用网格搜索法确定责任拟合函数中的待定系数,具体为:
1)计算每个用户每个时刻的动态谐波责任指标Ak,设第k个用户i时刻的动态谐波责任指标为Ak[1,i],表示为:
式中,x[i,i+T]表示从序列x的i时刻的数据到i+T时刻的数据,y[i,i+T]表示从序列y的i时刻的数据到i+T时刻的数据,则与动态相关系数D是一个只与责任拟合函数的系数取值有关的函数,用f表示,如式(5)所示:
4)采用网格搜索法求取系数,确定搜索步长,搜索网格中每一组a、b、c、d、e的取值,当式(5)达到最大值时的a、b、c、d、e即为所求系数。
2.根据权利要求1所述的一种长时间尺度动态谐波责任划分方法,其特征在于,
所述步骤S1具体为:采集需划分责任时段内谐波电压监测数据序列U=[u1,u2,…,un],根据此时段内连接于母线的各个用户在用电信息采集系统的平均有功功率计量数据,设共有m个用户,用户的平均有功功率序列为P1=[p11,p12,…,p1n]、P2=[p21,p22,…,p2n]、……、Pm=[pm1,pm2,…,pmn]。
3.根据权利要求1所述的一种长时间尺度动态谐波责任划分方法,其特征在于,所述采用网格搜索法求取系数具体为:
(1)对a、b、c、d、e、f赋初值,并且确定网格搜索下限为-100、网格搜索上限为+100、步长为1;
(2)从网格搜索下限开始,计算每个步长a、b、c、d、e的值,并代入式(5)计算f,比较此步f与前一步f的大小,若此步f大于前一步f,则记录此步a、b、c、d、e的值;若此步f小于前一步f,则记录前一步a、b、c、d、e的值,直到网格搜索上限;
(3)返回a、b、c、d、e的值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910110713.6A CN109829652B (zh) | 2019-02-12 | 2019-02-12 | 一种长时间尺度动态谐波责任划分方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910110713.6A CN109829652B (zh) | 2019-02-12 | 2019-02-12 | 一种长时间尺度动态谐波责任划分方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109829652A CN109829652A (zh) | 2019-05-31 |
CN109829652B true CN109829652B (zh) | 2023-04-07 |
Family
ID=66863472
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910110713.6A Active CN109829652B (zh) | 2019-02-12 | 2019-02-12 | 一种长时间尺度动态谐波责任划分方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109829652B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110618316B (zh) * | 2019-10-12 | 2021-09-03 | 南方电网科学研究院有限责任公司 | 一种谐波源辨识方法、装置、终端设备和存储介质 |
CN110927491B (zh) * | 2019-12-03 | 2021-07-20 | 南方电网科学研究院有限责任公司 | 基于无相位数据的多谐波源责任划分方法及装置 |
CN112230056B (zh) * | 2020-09-07 | 2022-04-26 | 国网河南省电力公司电力科学研究院 | 一种基于OFMMK-Means聚类和复合分位数回归多谐波源贡献计算方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130173189A1 (en) * | 2011-12-29 | 2013-07-04 | University of Connecticut, a public institution of higher education | Robust high resolution spectrum estimation method for accurate phasor, harmonic and interharmonic measurement in power systems |
CN105759117B (zh) * | 2015-10-15 | 2018-08-24 | 华北电力大学 | 一种基于谐波分析综合等值电路的谐波责任量化方法 |
CN108845194B (zh) * | 2018-03-30 | 2020-09-11 | 东南大学 | 基于幅值波动量峰度检测原理的谐波责任量化方法 |
CN109239456B (zh) * | 2018-08-03 | 2020-12-25 | 福州大学 | 一种基于动态规划时间序列相似性算法的谐波溯源方法 |
-
2019
- 2019-02-12 CN CN201910110713.6A patent/CN109829652B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109829652A (zh) | 2019-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111398885B (zh) | 一种结合线损分析的智能电表运行误差监控方法 | |
US10580096B2 (en) | Method for determining integrated network loss rate in UHV AC cross-regional electricity trading | |
CN109829652B (zh) | 一种长时间尺度动态谐波责任划分方法 | |
CN103278685B (zh) | 一种基于统计数据相关性分析的谐波责任定量划分方法 | |
CN106990324B (zh) | 一种配电网接地故障检测定位方法 | |
CN109239456B (zh) | 一种基于动态规划时间序列相似性算法的谐波溯源方法 | |
CN105098763A (zh) | 一种基于wams和ems的暂态电压稳定在线评估方法 | |
CN110579734B (zh) | 一种单相电能表电能质量分析的方法及单相电能表 | |
CN104156889A (zh) | 一种基于wams数据的风电场性能评估系统及其评估方法 | |
CN106443307A (zh) | 变电设备绝缘在线监测系统 | |
CN103926449A (zh) | 电力电缆接地电流的自适应监测方法 | |
CN108614192B (zh) | 一种配电网雷电监测装置布点决策方法 | |
CN117748507B (zh) | 基于高斯回归模型的配网谐波接入不确定性评估方法 | |
CN110047013B (zh) | 反专变用户间断型窃电方法 | |
Milanović et al. | Compound index for power quality evaluation and benchmarking | |
CN118091526B (zh) | 一种电流互感器测量误差在线监测方法及装置 | |
CN109146265A (zh) | 基于动态规划时间序列相似性的三相不平衡源溯源方法 | |
CN109613372B (zh) | 一种基于多元电网数据库的电网故障诊断方法 | |
Qiu et al. | Comprehensive evaluation of power quality based on improved TOPSIS-RSR method | |
Zhang et al. | A waveform-similarity-based protection scheme for the VSC-HVDC transmission lines | |
CN117856245A (zh) | 基于网格化管理的台区线损分析方法、装置及存储介质 | |
Li et al. | Estimation method for voltage sag occurrence frequency based on calculation of protection operating time characteristics | |
CN117013593A (zh) | 一种基于数字孪生的分布式电源评估方法及系统 | |
CN107064633B (zh) | 城市轨道交通负荷谐波电流迭加系数确定方法 | |
Gu et al. | Power quality early warning based on anomaly detection |
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 |