CN103969552B - 一种分布式发电系统的谐波源定位分析方法 - Google Patents

一种分布式发电系统的谐波源定位分析方法 Download PDF

Info

Publication number
CN103969552B
CN103969552B CN201410220507.8A CN201410220507A CN103969552B CN 103969552 B CN103969552 B CN 103969552B CN 201410220507 A CN201410220507 A CN 201410220507A CN 103969552 B CN103969552 B CN 103969552B
Authority
CN
China
Prior art keywords
harmonic
node
generation system
current
power generation
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
Application number
CN201410220507.8A
Other languages
English (en)
Other versions
CN103969552A (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.)
Hohai University HHU
Original Assignee
Hohai University HHU
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 Hohai University HHU filed Critical Hohai University HHU
Priority to CN201410220507.8A priority Critical patent/CN103969552B/zh
Publication of CN103969552A publication Critical patent/CN103969552A/zh
Application granted granted Critical
Publication of CN103969552B publication Critical patent/CN103969552B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种分布式发电系统的谐波源定位分析方法,包含以下内容,首先,根据分布式发电系统特点,获取系统参数以及各种基波和谐波量测数据;然后,通过基波和谐波状态估计方法计算出分布式发电系统所有节点的基波和各次谐波电压;接着,根据各支路两端节点的基波和各次谐波电压及其支路参数,分别计算出系统所有支路的基波和各次谐波电流;最后,通过支路电流与节点注入电流的关系公式计算出各节点的节点注入基波电流和各次谐波电流,并进一步计算出系统所有节点的电流总谐波畸变率分布情况,便可准确定位各谐波源分布的位置和分析其特性。本发明具有对分布式发电系统的多谐波源定位准确、计算过程清晰、易于编程等特点。

Description

一种分布式发电系统的谐波源定位分析方法
技术领域
本发明属于电力系统分析与计算领域,具体涉及一种分布式发电系统的谐波源定位分析方法。
背景技术
谐波源定位技术是谐波研究领域的重要分支,是了解电网谐波源分布特性和进行谐波治理的关键技术,可计算出谐波源的准确位置及其关键的谐波指标,进而提出抑制谐波的有效措施。随着分布式发电技术和主动配电网(ADN)技术迅速发展,越来越多的分布式电源(DG)接入配电网,形成了一个个分布式发电系统。而这些带有大量电力电子器件的分布式电源接入配电网,必将给配电网的电能质量带来新的问题和挑战。在分布式发电系统快速发展的背景下,谐波源定位技术是对分布式发电系统进行谐波治理的一项关键技术,是分析系统谐波源分布和特性的有效手段,有利于快速有效的实现谐波治理及改善电能质量。因此,研究分布式发电系统的多谐波源定位分析方法对分布式发电系统的进一步发展和应用具有重要意义。本发明从谐波源定位分析方法着手,提出一种分布式发电系统的谐波源定位分析方法,给分布式发电系统中的多谐波源定位和治理提供一种有效的技术方法。
发明内容
发明目的:针对上述现有技术存在的问题和不足,本发明的目的是提供一种分布式发电系统的谐波源定位分析方法,该方法能够准确定位分布式发电系统中的多谐波源位置,对研究分布式发电系统中各谐波源的特性和分布及其有效治理具有重要作用,从根本上解决了分布式发电系统中多谐波源准确定位的难题,极大提高了多谐波源定位的准确性和对其进行有效分析和治理的效率,具有对分布式发电系统的多谐波源定位准确、计算过程清晰、易于编程等特点。
技术方案:为实现上述发明目的,本发明采用的技术方案为一种分布式发电系统的谐波源定位分析方法,该方法包括以下步骤:
步骤A、获取分布式发电系统的网络参数以及各种基波和谐波量测数据,包括系统总节点数,节点序号,支路序号,独立节点数,支路数,参考节点,平衡节点,支路阻抗,网络拓扑结构,基波量测数据,以及各次谐波量测数据,如针对具有N个节点和含M个分布式电源的辐射型或树状分布式发电系统,假设系统的首节点作为参考节点和平衡节点,则系统独立节点数为n=N-1,支路数为b=n;
步骤B、通过常规基波和谐波状态估计方法计算出分布式发电系统各节点的基波电压和各次谐波电压其中,为节点i的基波电压,为节点i的第h次谐波电压,h为谐波次数,i为节点序号,i=1,2,…,n,n为系统独立节点数,上标“(1)”表示基波分量,上标“(h)”表示第h次谐波分量;
步骤C、根据公式计算分布式发电系统各支路的基波电流和各次谐波电流其中,分别为分布式发电系统支路r的基波电流和第h次谐波电流,分别为分布式发电系统支路r的基波阻抗和第h次谐波阻抗,分别为分布式发电系统支路r的首端节点j和末端节点i的基波电压,分别为分布式发电系统支路r的首端节点j和末端节点i的第h次谐波电压,j和i分别为支路r的首端节点和末端节点序号,r为支路序号,j、i=1,2,…,n,r=1,2,…,n;
步骤D、根据公式计算分布式发电系统各节点的节点注入基波电流和各次谐波电流其中,分别为分布式发电系统独立节点i的节点注入基波电流和第h次谐波电流,为由各独立节点i注入基波电流组成的n×1阶相量矩阵,即 为由各支路r基波电流组成的n×1阶相量矩阵,即 为由各独立节点i注入第h次谐波电流组成的n×1阶相量矩阵,即 为由各支路r第h次谐波电流组成的n×1阶相量矩阵,即A为分布式发电系统的n×n阶节-支关联矩阵,i=1,2,…,n,r=1,2,…,n;
步骤E、根据不同谐波次数h,重复步骤B、步骤C和步骤D即可求出分布式发电系统所有节点的节点注入全部谐波电流,并计算出分布式发电系统各谐波源向系统输入的各次谐波电流为其中,为各谐波源向分布式发电系统输入的第h次谐波电流,i为谐波源节点序号,i=1,2,…,n;
步骤F、根据电流总谐波畸变率THD的计算公式和电流单个谐波畸变率IHD的计算公式分别计算出分布式发电系统中各谐波源的电流总谐波畸变率和单个谐波畸变率的分布情况,从而准确定位分布式发电系统中各谐波源的节点序号和位置,并根据谐波源的各次谐波含量以及各谐波源的总谐波畸变率分布水平采取相应的治理措施,最终实现分布式发电系统多谐波源准确定位及有效治理,其中,THDi(%)表示谐波源的电流总谐波畸变率,表示谐波源的电流单个谐波畸变率,h为谐波次数,H为谐波源的最高谐波次数,为各谐波源输出的基波电流,且有i为谐波源节点序号,i=1,2,…,n。
有益效果:本发明针对现有谐波源定位技术存在的不足,提出了一种分布式发电系统的谐波源定位分析方法,该方法根据分布式发电系统特点,获取系统参数以及各种基波和谐波量测数据;然后,通过基波和谐波状态估计方法计算出分布式发电系统所有节点的基波和各次谐波电压;接着,根据各支路两端节点的基波和各次谐波电压及其支路参数,分别计算出系统所有支路的基波和各次谐波电流;最后,通过支路电流与节点注入电流的关系公式计算出各节点的节点注入基波电流和各次谐波电流,并进一步计算出系统所有节点的电流总谐波畸变率分布情况,便可准确定位各谐波源分布的位置和分析其特性。该方法能够准确定位分布式发电系统中的多谐波源位置,对研究分布式发电系统中各谐波源的特性和分布及其有效治理具有重要作用,从根本上解决了分布式发电系统中多谐波源准确定位的难题,极大提高了多谐波源定位的准确性和对其进行有效分析和治理的效率,具有对分布式发电系统的多谐波源定位准确、计算过程清晰、易于编程等特点。由此可见,本发明具有很好的工程应用价值和借鉴意义。
附图说明
图1为本发明的总体流程图。
具体实施方式
下面结合附图和具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
图1为本发明的总体流程图,具体方法如下:
针对具有N个节点和含M个分布式电源(DG)的辐射型或树状分布式发电系统,假设系统的首节点作为参考节点和平衡节点,则系统独立节点数为n=N-1,支路数(也即树支支路数)为b=n。
辐射型分布式发电系统的有向连通图有N个节点,b条支路,独立节点数为n=N-1,支路数等于独立节点数,即b=n,在不考虑参考节点的情况下,则其节点-支路(下面简称“节-支”)关联矩阵A描述独立节点和支路之间的联系。其中节-支关联矩阵A是一个n×n阶矩阵,规定每条支路的正方向都是从小号节点指向大号节点,如果节点i是支路r的首端节点,则A(i,r)=1,如果节点i是支路r的末端节点,则A(i,r)=-1,如果节点i不是支路r的端节点,则A(i,r)=0。
1)获取分布式发电系统的网络参数以及各种基波和谐波量测数据,包括系统总节点数,节点序号,支路序号,独立节点数,支路数,参考节点,平衡节点,支路阻抗,网络拓扑结构,基波量测数据,以及各次谐波量测数据。
2)通过常规基波和谐波状态估计方法计算出分布式发电系统各节点的基波电压和各次谐波电压其中,为节点i的基波电压,为节点i的第h次谐波电压,h为谐波次数,i为节点序号,i=1,2,…,n,n为系统独立节点数,上标“(1)”表示基波分量,上标“(h)”表示第h次谐波分量。
3)根据公式计算分布式发电系统各支路的基波电流和各次谐波电流其中,分别为分布式发电系统支路r的基波电流和第h次谐波电流,分别为分布式发电系统支路r的基波阻抗和第h次谐波阻抗,分别为分布式发电系统支路r的首端节点j和末端节点i的基波电压,分别为分布式发电系统支路r的首端节点j和末端节点i的第h次谐波电压,j和i分别为支路r的首端节点和末端节点序号,r为支路序号,j、i=1,2,…,n,r=1,2,…,n。
4)根据公式计算分布式发电系统各节点的节点注入基波电流和各次谐波电流其中,分别为分布式发电系统独立节点i的节点注入基波电流和第h次谐波电流,为由各独立节点i注入基波电流组成的n×1阶相量矩阵,即 为由各支路r基波电流组成的n×1阶相量矩阵,即 为由各独立节点i注入第h次谐波电流组成的n×1阶相量矩阵,即 为由各支路r第h次谐波电流组成的n×1阶相量矩阵,即A为分布式发电系统的n×n阶节-支关联矩阵,i=1,2,…,n,r=1,2,…,n。
5)根据不同谐波次数h,重复步骤2)、步骤3)和步骤4)即可求出分布式发电系统所有节点的节点注入全部谐波电流,并计算出分布式发电系统各谐波源向系统输入的各次谐波电流为其中,为各谐波源向分布式发电系统输入的第h次谐波电流,i为谐波源节点序号,i=1,2,…,n。
6)根据电流总谐波畸变率THD的计算公式和电流单个谐波畸变率IHD的计算公式分别计算出分布式发电系统中各谐波源的电流总谐波畸变率和单个谐波畸变率的分布情况,从而准确定位分布式发电系统中各谐波源的节点序号和位置,并根据谐波源的各次谐波含量以及各谐波源的总谐波畸变率分布水平采取相应的治理措施,最终实现分布式发电系统多谐波源准确定位及有效治理。
步骤6)中电流总谐波畸变率THD的计算公式为:
THD i ( % ) = Σ h = 2 H | I · i ( h ) | 2 | I · i ( 1 ) | × 100 % - - - ( 1 )
电流单个谐波畸变率IHD的计算公式为:
IHD i ( h ) ( % ) = | I · i ( h ) | | I · i ( 1 ) | × 100 % - - - ( 2 )
其中,THDi(%)表示谐波源的电流总谐波畸变率,表示谐波源的电流单个谐波畸变率,h为谐波次数,H为谐波源的最高谐波次数,为各谐波源输出的基波电流,且有i为谐波源节点序号,i=1,2,…,n。

Claims (2)

1.一种分布式发电系统的谐波源定位分析方法,其特征是该方法包括以下步骤:
步骤A、获取分布式发电系统的网络参数以及各种基波和谐波量测数据,针对具有N个节点和含M个分布式电源的辐射型或树状分布式发电系统,假设系统的首节点作为参考节点和平衡节点,则系统独立节点数为n=N-1,支路数为b=n;
步骤B、通过常规基波和谐波状态估计方法计算出分布式发电系统各节点的基波电压和各次谐波电压其中,为节点i的基波电压,为节点i的第h次谐波电压,h为谐波次数,i为节点序号,i=1,2,…,n,n为系统独立节点数,上标“(1)”表示基波分量,上标“(h)”表示第h次谐波分量;
步骤C、根据公式计算分布式发电系统各支路的基波电流和各次谐波电流其中,分别为分布式发电系统支路r的基波电流和第h次谐波电流,分别为分布式发电系统支路r的基波阻抗和第h次谐波阻抗,分别为分布式发电系统支路r的首端节点j和末端节点i的基波电压,分别为分布式发电系统支路r的首端节点j和末端节点i的第h次谐波电压,j和i分别为支路r的首端节点和末端节点序号,r为支路序号,j=1,2,…,n,i=1,2,…,n,r=1,2,…,n;
步骤D、根据公式计算分布式发电系统各节点的节点注入基波电流和各次谐波电流其中,分别为分布式发电系统独立节点i的节点注入基波电流和第h次谐波电流,为由各独立节点i注入基波电流组成的n×1阶相量矩阵,即 为由各支路r基波电流组成的n×1阶相量矩阵,即 为由各独立节点i注入第h次谐波电流组成的n×1阶相量矩阵,即 为由各支路r第h次谐波电流组成的n×1阶相量矩阵,即A为分布式发电系统的n×n阶节-支关联矩阵,i=1,2,…,n,r=1,2,…,n;
步骤E、根据不同谐波次数h,重复步骤B、步骤C和步骤D即可求出分布式发电系统所有的节点注入谐波电流,并计算出分布式发电系统各谐波源向系统输入的各次谐波电流为其中,为各谐波源向分布式发电系统输入的第h次谐波电流,i为谐波源节点序号,i=1,2,…,n;
步骤F、根据电流总谐波畸变率THD的计算公式和电流单个谐波畸变率IHD的计算公式分别计算出分布式发电系统中各谐波源的电流总谐波畸变率和单个谐波畸变率的分布情况,从而准确定位分布式发电系统中各谐波源的节点序号和位置,并根据谐波源的各次谐波含量以及各谐波源的总谐波畸变率分布水平采取相应的治理措施,最终实现分布式发电系统多谐波源准确定位及有效治理,其中,THDi(%)表示谐波源的电流总谐波畸变率,表示谐波源的电流单个谐波畸变率,h为谐波次数,H为谐波源的最高谐波次数,为各谐波源输出的基波电流,且有i为谐波源节点序号,i=1,2,…,n。
2.如权利要求1所述的分布式发电系统的谐波源定位分析方法,其特征在于,步骤A中所述分布式发电系统的网络参数包括系统总节点数、节点序号、支路序号、独立节点数、支路数、参考节点、平衡节点、支路阻抗和网络拓扑结构,谐波量测数据包括各次谐波量测数据。
CN201410220507.8A 2014-05-22 2014-05-22 一种分布式发电系统的谐波源定位分析方法 Active CN103969552B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410220507.8A CN103969552B (zh) 2014-05-22 2014-05-22 一种分布式发电系统的谐波源定位分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410220507.8A CN103969552B (zh) 2014-05-22 2014-05-22 一种分布式发电系统的谐波源定位分析方法

Publications (2)

Publication Number Publication Date
CN103969552A CN103969552A (zh) 2014-08-06
CN103969552B true CN103969552B (zh) 2016-08-24

Family

ID=51239293

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410220507.8A Active CN103969552B (zh) 2014-05-22 2014-05-22 一种分布式发电系统的谐波源定位分析方法

Country Status (1)

Country Link
CN (1) CN103969552B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104407214A (zh) * 2014-11-27 2015-03-11 云南电网公司电力科学研究院 一种谐波源辨识方法
CN104849619B (zh) * 2015-05-26 2018-01-19 中国电力科学研究院 一种动态谐波源定位方法
CN104897992A (zh) * 2015-06-08 2015-09-09 湖南城市学院 一种分布式电源扰动分类装置及其在线检测方法
CN105510719A (zh) * 2015-11-30 2016-04-20 云南电网有限责任公司玉溪供电局 三相电网谐波阻抗测量方法
CN105717360A (zh) * 2016-01-26 2016-06-29 云南电网有限责任公司电力科学研究院 一种低压单相谐波阻抗测量方法及装置
CN106597099B (zh) * 2016-12-15 2019-02-01 国网新疆电力公司电力科学研究院 一种基于参数估计的谐波分析系统
CN107968402A (zh) * 2017-12-11 2018-04-27 华北电力大学 一种多电流源型振荡源的定位方法
CN109541308B (zh) * 2019-01-25 2020-12-04 中国科学院电工研究所 基于分布式准同步电网测量系统的谐波分析方法
CN114236281A (zh) * 2021-12-14 2022-03-25 福州大学 工业用户谐波污染分布情况辨识方法
CN114236235B (zh) * 2021-12-17 2024-02-06 福州大学 一种采用相似特征匹配的局部不可观系统谐波源定位方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201926727U (zh) * 2011-01-12 2011-08-10 烟台伊科电气技术有限公司 电能质量监测仪
CN102279327A (zh) * 2011-05-06 2011-12-14 重庆大学 光伏并网发电在线监测与状态评估系统
CN102323494A (zh) * 2011-05-24 2012-01-18 山东电力研究院 一种多谐波源辨识方法
CN202230173U (zh) * 2011-10-14 2012-05-23 山东电力研究院 自适应多谐波源定位装置
CN102540009A (zh) * 2010-12-23 2012-07-04 上海蓝瑞电气有限公司 一种配电网故障定位系统
CN102621454A (zh) * 2012-04-06 2012-08-01 河海大学 谐波源闪变源实时定位仪及确定谐波源闪变源方位的方法
CN202512198U (zh) * 2012-04-06 2012-10-31 河海大学 谐波源闪变源实时定位仪
CN102866315A (zh) * 2012-09-17 2013-01-09 华南理工大学 含逆变型分布式电源配电网的对称故障分析方法
CN103424620A (zh) * 2013-08-08 2013-12-04 华南理工大学 一种基于时间序列相关度的10kV配网谐波源辨识方法
CN103675446A (zh) * 2013-12-23 2014-03-26 重庆大学 一种实时精确的谐波源定位系统和方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3347995B2 (ja) * 1997-12-11 2002-11-20 株式会社日立製作所 高調波計測装置
CN102570460A (zh) * 2010-12-14 2012-07-11 中国电力科学研究院 一种计算电气化铁路的谐波和负序电流的方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102540009A (zh) * 2010-12-23 2012-07-04 上海蓝瑞电气有限公司 一种配电网故障定位系统
CN201926727U (zh) * 2011-01-12 2011-08-10 烟台伊科电气技术有限公司 电能质量监测仪
CN102279327A (zh) * 2011-05-06 2011-12-14 重庆大学 光伏并网发电在线监测与状态评估系统
CN102323494A (zh) * 2011-05-24 2012-01-18 山东电力研究院 一种多谐波源辨识方法
CN202230173U (zh) * 2011-10-14 2012-05-23 山东电力研究院 自适应多谐波源定位装置
CN102621454A (zh) * 2012-04-06 2012-08-01 河海大学 谐波源闪变源实时定位仪及确定谐波源闪变源方位的方法
CN202512198U (zh) * 2012-04-06 2012-10-31 河海大学 谐波源闪变源实时定位仪
CN102866315A (zh) * 2012-09-17 2013-01-09 华南理工大学 含逆变型分布式电源配电网的对称故障分析方法
CN103424620A (zh) * 2013-08-08 2013-12-04 华南理工大学 一种基于时间序列相关度的10kV配网谐波源辨识方法
CN103675446A (zh) * 2013-12-23 2014-03-26 重庆大学 一种实时精确的谐波源定位系统和方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
电力系统谐波源定位方法评述;徐志向等;《电工电能新技术》;20060731;第25卷(第3期);第64-67页 *

Also Published As

Publication number Publication date
CN103969552A (zh) 2014-08-06

Similar Documents

Publication Publication Date Title
CN103969552B (zh) 一种分布式发电系统的谐波源定位分析方法
CN103956735B (zh) 一种分布式发电系统的谐波潮流分析方法
CN105656038B (zh) 一种计及电网接纳能力的概率潮流计算方法
US9627886B2 (en) State estimation for power system using hybrid measurements
CN104573330B (zh) 氮化镓高电子迁移率晶体管i‑v模型参数的提取方法
CN106980044B (zh) 一种适应风电接入的电力系统谐波电流估计方法
CN106908665B (zh) 一种单相交流系统dq分解阻抗测量方法
CN103983851B (zh) 主动配电网谐波源定位分析方法
CN106786608B (zh) 一种适用于分布式电源接入的不确定谐波潮流计算方法
WO2016061745A1 (zh) 交流电力网的最小相位线性有功潮流的获取方法
CN111339624B (zh) 基于psasp和emtp/atp短路电流直流分量计算方法
CN103605856B (zh) 基于分数阶线路模型的输电线路参数估计方法
CN110556854B (zh) 一种柔性直流输电系统并网点的静态谐波分析方法及系统
CN104167732B (zh) 一种基于联络线相角差的电网等值方法
Meerimatha et al. A novel distribution system power flow algorithm using forward backward matrix method
CN105277852A (zh) 配电网线路工况分类识别方法
CN106546824B (zh) 基于变系数回归的含光伏电网系统谐波阻抗估计方法
CN103972886B (zh) 一种主动配电网谐波潮流分析方法
CN104849551B (zh) 一种谐相角分析方法
CN103576120A (zh) 三次谐波分量准同步信息传输校验与自愈算法
CN104092213A (zh) 一种基于优化方法的不确定潮流支路功率分析方法
CN110896218B (zh) 一种用于建立集合性居民负荷的谐波建模方法及系统
CN104764930B (zh) 基于相关特性的谐波发射水平估计方法
CN103956774B (zh) 一种含弱环分布式发电系统的谐波潮流分析方法
Tan et al. Convergence and stability of the compensated split-step θ-method for stochastic differential equations with jumps

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant