CN108802491B - 基于谐波源独立性的谐波发射水平估计求解新方法 - Google Patents

基于谐波源独立性的谐波发射水平估计求解新方法 Download PDF

Info

Publication number
CN108802491B
CN108802491B CN201810735975.7A CN201810735975A CN108802491B CN 108802491 B CN108802491 B CN 108802491B CN 201810735975 A CN201810735975 A CN 201810735975A CN 108802491 B CN108802491 B CN 108802491B
Authority
CN
China
Prior art keywords
pcc
harmonic
covariance matrix
matrix
equation
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.)
Expired - Fee Related
Application number
CN201810735975.7A
Other languages
English (en)
Other versions
CN108802491A (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.)
Sichuan University
Original Assignee
Sichuan University
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 Sichuan University filed Critical Sichuan University
Priority to CN201810735975.7A priority Critical patent/CN108802491B/zh
Publication of CN108802491A publication Critical patent/CN108802491A/zh
Application granted granted Critical
Publication of CN108802491B publication Critical patent/CN108802491B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16Spectrum analysis; Fourier analysis

Landscapes

  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明公开了基于谐波源独立性的谐波发射水平估计求解新方法,所述方法包括以下步骤:步骤(A):对一定时间段内的公共耦合点电压、电流值进行采样,使用采用的量测数据建立量测数据矩阵;步骤(B):对步骤(A)中的量测数据矩阵,计算其协方差矩阵;同时根据已知的电路方程,建立其协方差矩阵的另一个表达式,根据两个协方差矩阵表达式的对应关系,建立阻抗计算方程;步骤(C):基于PCC点两侧谐波源独立的特性,对步骤(B)中的协方差矩阵的另一个表达式进行简化,由两协方差矩阵表达式对应相等建立方程组求解谐波阻抗。

Description

基于谐波源独立性的谐波发射水平估计求解新方法
技术领域
本发明涉及谐波源定位中的谐波发射水平研究领域,具体地,涉及一种基于谐波源独立性的谐波发射水平估计求解新方法。
背景技术
随着电网的不断发展,越来越多的电力电子元件接入了电网中,从而引发的谐波问题越来越引起重视,如何定量划分接入系统中各谐波源的谐波责任是实现电能质量智能管理、控制的前提。
通常采用等效电路的方法。即选定一个公共连接点(PCC),将PCC两端的系统侧电路和用户侧电路等效成谐波电压源和谐波电流源,并定义系统侧和用户侧谐波参考阻抗,将谐波责任表示成PCC处量测电压、电流和参考阻抗,并划分两侧的谐波责任。
目前进行谐波责任划分的方法主要有波动量法、线性回归方法、独立分量法等。
综上所述,本申请发明人在实现本申请发明技术方案的过程中,发现上述技术至少存在如下技术问题:
在现有技术中,现有的波动量方法、线性回归方法存在估计误差随着两侧谐波源发射水平逐渐接近,误差增长速度较快的技术问题。波动量法根据波动量比值特征符号进行估计,在某些场景下估计精度不能保证;线性回归方法适用于系统侧和用户侧平稳的条件,在波动量较大时拟合效果不好。
发明内容
本发明的目的在于提供一种根据谐波源独立性计算谐波阻抗的计算方法,依据谐波源之间的独立性,通过PCC处测量得到的电压、电流波动量,建立量测数据矩阵并计算其协方差矩阵;同时通过等效电路建立PCC处测量电压、电流与谐波源电流波动量、参考阻抗的方程;由系统侧谐波源电流波动量远小于用户侧谐波源电流波动量的特性,可对方程进行简化,并求解谐波阻抗和谐波发射水平。
为实现上述发明目的,本申请提供了基于谐波源独立性的谐波发射水平估计求解新方法,包括以下步骤:
步骤(A):对一定时间段内的公共耦合点电压、电流值进行采样,使用采用的量测数据建立量测数据矩阵;
步骤(B):对步骤(A)中的量测数据矩阵,计算其协方差矩阵;同时根据已知的电路方程,建立其协方差矩阵的另一个表达式,根据两个协方差矩阵表达式的对应关系,建立阻抗计算方程;
步骤(C):基于PCC点两侧谐波源独立的特性,对步骤(B)中的协方差矩阵的另一个表达式进行简化,由两协方差矩阵表达式对应相等建立方程组求解谐波阻抗。
进一步的,所述步骤(A)中的量测数据矩阵建立按照以下步骤进行:
步骤(A1):假设采样电压、电流为Vpcc、Ipcc,采样点数为N;对采样值去均值后,建立量测数据矩阵A;
上式中ΔVpcc为PCC点电压波动量,ΔIpcc为PCC点电流波动量;
步骤(A2):上式中的ΔVpcc、ΔIpcc可由下式计算得到:
ΔVpcc=Vpcc(n+1)-Vpcc(n)n=1,2,...,N-1;
ΔIpcc=Ipcc(n+1)-Ipcc(n)n=1,2,...,N-1。
进一步的,所述步骤(B)中的协方差矩阵的两个表达式的建立按照以下步骤进行:
步骤(B1):量测数据矩阵A的协方差矩阵(1)建立如下:
上式中上注符H表示共轭转置;上注符*表示共轭;
步骤(B2):根据PCC诺顿等效电路可建立方程:
A=ZI
上式中,Zu、Zc分别为PCC诺顿等效电路系统侧和用户侧的等效阻抗;ΔIu和ΔIc分别为系统侧和用户侧电流波动量,Z为等效阻抗矩阵,I为PCC两侧电流波动量矩阵;
步骤(B3):由电路方程改写量测数据矩阵A的协方差矩阵的另一个表达式(2)如下:
AAH=(ZI)(ZI)H=Z(IIH)ZH (2)。
进一步的,所述步骤(C)中的协方差矩阵简化和谐波阻抗计算按照以下步骤进行:
步骤(C1):IIH可简化如下:
其中,λ1、λ2为PCC系统侧和用户侧的电流波动量幅值平方和,可由下式表示:
步骤(C2):由(B2)计算协方差矩阵AAH如下:
步骤(C3):联立(B1)和(C2)中协方差矩阵AAH计算方程可建立方程组:
步骤(C4):(C3)中方程可简化如下,求解简化后方程即可求解谐波阻抗并计算谐波发射水平:
本发明与现有技术相比,具有如下的优点和有益效果:
本发明基于谐波源独立性的谐波发射水平估计求解新方法,在建立量测数据协方差矩阵时可忽略谐波源电流波动量之间交叉相乘项,减少了PCC两侧谐波源的相互影响对计算结果的影响;2本发明基于谐波源独立性的谐波发射水平估计求解新方法,可减少谐波责任计算的计算量,加快计算速度;3本发明基于谐波源独立性的谐波发射水平估计求解新方法,算法计算结果较为准确快速,误差较小,具有良好的实用性。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解,构成本申请的一部分,并不构成对本发明实施例的限定;
图1为本发明流程示意图;
图2为实施例一的PCC点等效电路模型图。
具体实施方式
本发明提供了一种基于双目立体视觉的接触网几何参数检测测量系统,解决了现有的基于单目视觉的检测方法存在先天的误差的技术问题,实现了测量系统设计合理,几何参数检测测量结果准确的技术效果。
为了能够更清楚地理解本发明的上述目的、特征和优点,下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是,在相互不冲突的情况下,本申请的实施例及实施例中的特征可以相互组合。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是,本发明还可以采用其他不同于在此描述范围内的其他方式来实施,因此,本发明的保护范围并不受下面公开的具体实施例的限制。
请参考图1,本申请提供了基于谐波源独立性的谐波发射水平估计求解新方法,所述方法包括以下步骤:
步骤(A):对一定时间段内的公共耦合点电压、电流值进行采样,使用采用的量测数据建立量测数据矩阵;
步骤(B):对步骤(A)中的量测数据矩阵,计算其协方差矩阵;同时根据已知的电路方程,建立其协方差矩阵的另一个表达式,根据两个协方差矩阵表达式的对应关系,建立阻抗计算方程;
步骤(C):基于PCC点两侧谐波源独立的特性,对步骤(B)中的协方差矩阵的另一个表达式进行简化,由两协方差矩阵表达式对应相等建立方程组求解谐波阻抗;
实施例一:在Matlab中搭建如图2所示的PCC点等效电路,各元件参数表1所示,使用本发明对谐波阻抗进行求解。并将计算结果、计算误差与其他方法进行比较如表2-表3所示。
表1
表2
表3
其中,在本申请实施例中,所述步骤(A)中的量测数据矩阵建立按照以下步骤进行:
步骤(A1):假设采样电压、电流为Vpcc、Ipcc,采样点数为N;对采样值去均值后,建立量测数据矩阵A;
式中ΔVpcc为PCC点电压波动量,ΔIpcc为PCC点电流波动量。
步骤(A2):式(1)中的ΔVpcc、ΔIpcc可由下式计算得到:
ΔVpcc=Vpcc(n+1)-Vpcc(n)n=1,2,...,N-1;
ΔIpcc=Ipcc(n+1)-Ipcc(n)n=1,2,...,N-1;
其中,在本申请实施例中,所述步骤(B)中的协方差矩阵的两个表达式的建立按照以下步骤进行:
步骤(B1):量测数据矩阵A的协方差矩阵(1)建立如下:
式中上注符H表示共轭转置;上注符*表示共轭。
步骤(B2):根据PCC诺顿等效电路可建立方程:
A=ZI
式中,Zu、Zc分别为PCC诺顿等效电路系统侧和用户侧的等效阻抗;ΔIu和ΔIc分别为系统侧和用户侧电流波动量,Z为等效阻抗矩阵,I为PCC两侧电流波动量矩阵。
步骤(B3):由电路方程改写量测数据矩阵A的协方差矩阵的另一个表达式(2)如下:
AAH=(ZI)(ZI)H=Z(IIH)ZH (2)
其中,在本申请实施例中,所述步骤(C)中的协方差矩阵简化和谐波阻抗计算按照以下步骤进行:
步骤(C1):IIH可简化如下:
其中,λ1、λ2为PCC系统侧和用户侧的电流波动量幅值平方和,可由下式表示:
步骤(C2):由(B2)计算协方差矩阵AAH如下:
步骤(C3):联立(B1)和(C2)中协方差矩阵AAH计算方程可建立方程组:
步骤(C4):(C3)中方程可简化如下,求解简化后方程即可求解谐波阻抗并计算谐波发射水平:
尽管已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。

Claims (1)

1.基于谐波源独立性的谐波发射水平估计求解新方法,其特征在于,所述方法包括以下步骤:
步骤(A):对一定时间段内的公共耦合点电压、电流值进行采样,使用采用的量测数据建立量测数据矩阵;
步骤(B):对步骤(A)中的量测数据矩阵,计算其协方差矩阵;同时根据已知的电路方程,建立其协方差矩阵的另一个表达式,根据两个协方差矩阵表达式的对应关系,建立阻抗计算方程;
步骤(C):基于PCC点两侧谐波源独立的特性,对步骤(B)中的协方差矩阵的另一个表达式进行简化,由两协方差矩阵表达式对应相等建立方程组求解谐波阻抗;
所述步骤(A)中的量测数据矩阵建立按照以下步骤进行:
步骤(A1):假设采样电压、电流为Vpcc、Ipcc,采样点数为N;对采样值去均值后,建立量测数据矩阵A;
上式中ΔVpcc为PCC点电压波动量,ΔIpcc为PCC点电流波动量;
步骤(A2):上式中的ΔVpcc、ΔIpcc可由下式计算得到:
ΔVpcc=Vpcc(n+1)-Vpcc(n)n=1,2,…,N-1;
ΔIpcc=Ipcc(n+1)-Ipcc(n)n=1,2,…,N-1;
所述步骤(B)中的协方差矩阵的两个表达式的建立按照以下步骤进行:
步骤(B1):量测数据矩阵A的协方差矩阵(1)建立如下:
上式中上注符H表示共轭转置;上注符*表示共轭;
步骤(B2):根据PCC诺顿等效电路可建立方程:
A=ZI
上式中,Zu、Zc分别为PCC诺顿等效电路系统侧和用户侧的等效阻抗;ΔIu和ΔIc分别为系统侧和用户侧电流波动量,Z为等效阻抗矩阵,I为PCC两侧电流波动量矩阵;
步骤(B3):由电路方程改写量测数据矩阵A的协方差矩阵的另一个表达式(2)如下:
AAH=(ZI)(ZI)H=Z(IIH)ZH (2);
所述步骤(C)中的协方差矩阵简化和谐波阻抗计算按照以下步骤进行:
步骤(C1):IIH可简化如下:
其中,λ1、λ2为PCC系统侧和用户侧的电流波动量幅值平方和,可由下式表示:
步骤(C2):由(B2)计算协方差矩阵AAH如下:
步骤(C3):联立(B1)和(C2)中协方差矩阵AAH计算方程可建立方程组:
步骤(C4):(C3)中方程可简化如下,求解简化后方程即可求解谐波阻抗并计算谐波发射水平:
CN201810735975.7A 2018-07-06 2018-07-06 基于谐波源独立性的谐波发射水平估计求解新方法 Expired - Fee Related CN108802491B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810735975.7A CN108802491B (zh) 2018-07-06 2018-07-06 基于谐波源独立性的谐波发射水平估计求解新方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810735975.7A CN108802491B (zh) 2018-07-06 2018-07-06 基于谐波源独立性的谐波发射水平估计求解新方法

Publications (2)

Publication Number Publication Date
CN108802491A CN108802491A (zh) 2018-11-13
CN108802491B true CN108802491B (zh) 2019-07-19

Family

ID=64075323

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810735975.7A Expired - Fee Related CN108802491B (zh) 2018-07-06 2018-07-06 基于谐波源独立性的谐波发射水平估计求解新方法

Country Status (1)

Country Link
CN (1) CN108802491B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110045186B (zh) * 2019-04-25 2020-06-23 四川大学 电力系统公共连接点处的系统谐波阻抗估计方法和系统
CN110261682A (zh) * 2019-08-14 2019-09-20 四川大学 实测数据相位缺失下的系统侧谐波阻抗获取方法和系统
CN110441601B (zh) * 2019-10-10 2020-01-07 四川大学 基于独立分量分析的多谐波源责任量化方法
CN115032458B (zh) * 2022-06-02 2023-08-11 北京妙微科技有限公司 一种基于电网系统的谐波阻抗估计方法及计算机存储介质

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102135569A (zh) * 2011-01-21 2011-07-27 清华大学 基于波动量法的用户侧谐波发射水平实用化的估计方法
CN102998535A (zh) * 2012-12-10 2013-03-27 华北电力大学(保定) 一种基于极大似然估计理论的系统谐波阻抗计算方法
CN103344828A (zh) * 2013-07-12 2013-10-09 重庆大学 一种新的谐波发射水平评估方法
CN104407214A (zh) * 2014-11-27 2015-03-11 云南电网公司电力科学研究院 一种谐波源辨识方法
CN104698273A (zh) * 2015-03-19 2015-06-10 广州供电局有限公司 谐波责任划分方法和系统
CN104898019A (zh) * 2015-05-11 2015-09-09 国家电网公司 一种应用于主动配电网的谐波源定位和谐波责任划分方法
CN105548809A (zh) * 2015-12-18 2016-05-04 南瑞(武汉)电气设备与工程能效测评中心 一种电力系统谐波责任的确定方法
CN107576853A (zh) * 2017-07-12 2018-01-12 国网上海市电力公司 基于典型相关性分析的配电网谐波阻抗计算方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102135569A (zh) * 2011-01-21 2011-07-27 清华大学 基于波动量法的用户侧谐波发射水平实用化的估计方法
CN102998535A (zh) * 2012-12-10 2013-03-27 华北电力大学(保定) 一种基于极大似然估计理论的系统谐波阻抗计算方法
CN103344828A (zh) * 2013-07-12 2013-10-09 重庆大学 一种新的谐波发射水平评估方法
CN104407214A (zh) * 2014-11-27 2015-03-11 云南电网公司电力科学研究院 一种谐波源辨识方法
CN104698273A (zh) * 2015-03-19 2015-06-10 广州供电局有限公司 谐波责任划分方法和系统
CN104898019A (zh) * 2015-05-11 2015-09-09 国家电网公司 一种应用于主动配电网的谐波源定位和谐波责任划分方法
CN105548809A (zh) * 2015-12-18 2016-05-04 南瑞(武汉)电气设备与工程能效测评中心 一种电力系统谐波责任的确定方法
CN107576853A (zh) * 2017-07-12 2018-01-12 国网上海市电力公司 基于典型相关性分析的配电网谐波阻抗计算方法

Also Published As

Publication number Publication date
CN108802491A (zh) 2018-11-13

Similar Documents

Publication Publication Date Title
CN108802491B (zh) 基于谐波源独立性的谐波发射水平估计求解新方法
Chan et al. The available capacity computation model based on artificial neural network for lead–acid batteries in electric vehicles
CN105912799B (zh) 一种液态或半液态金属电池的建模方法
CN103197251B (zh) 一种动力锂电池二阶rc等效模型的辨识方法
CN104698273B (zh) 谐波责任划分方法和系统
CN107658881A (zh) 基于戴维南等值方法的电压稳定临界点判断方法
CN106483462A (zh) 一种电池充电电量的测量方法和装置
CN109633479A (zh) 基于嵌入式容积卡尔曼滤波的锂电池soc在线估算方法
CN102798760A (zh) 一种光伏阵列对地绝缘阻抗检测方法及电路
CN106093571B (zh) 基于阻抗约束的判断用户侧为主谐波源的谐波溯源方法
CN110441601A (zh) 基于独立分量分析的多谐波源责任量化方法
CN106300331B (zh) 一种精确计算配电网支路阻抗的方法
CN102095934A (zh) 一种交流电信号的相位差测量方法
CN102545177A (zh) 一种基于贝杰龙模型的交流输电线路故障选相的测后模拟方法
CN110350517A (zh) 一种基于运行风险的电动汽车并网配电网络重构方法
CN109888773A (zh) 一种电力系统多区域分布式状态评估方法
CN106443496A (zh) 一种带改进型噪声估计器的电池荷电状态估计方法
CN111141995A (zh) 一种基于比幅原理的线路双端稳态量测距方法和系统
CN106990275A (zh) 对地绝缘电阻电压采样方法及装置
CN104979829A (zh) 一种基于改进的FastICA算法的谐波电流估计方法
CN115000947A (zh) 基于智能电表量测的配电网拓扑结构与线路参数辨识方法
CN104362638A (zh) 基于相量量测单元量测的电网聚合的关键节点电压调控方法
CN108318819A (zh) 一种估算电池荷电状态的方法
CN107064808B (zh) 电池的soc估计方法及装置
CN111123027A (zh) 基于wams系统的多端传输线路故障测距方法

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190719

Termination date: 20210706