CN106383295B - 一种供用电系统谐波源识别方法及模型 - Google Patents
一种供用电系统谐波源识别方法及模型 Download PDFInfo
- Publication number
- CN106383295B CN106383295B CN201610946158.7A CN201610946158A CN106383295B CN 106383295 B CN106383295 B CN 106383295B CN 201610946158 A CN201610946158 A CN 201610946158A CN 106383295 B CN106383295 B CN 106383295B
- Authority
- CN
- China
- Prior art keywords
- harmonic
- power
- load
- branch
- source
- 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
- 230000005611 electricity Effects 0.000 title claims abstract description 83
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000012544 monitoring process Methods 0.000 claims abstract description 7
- 238000005070 sampling Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 6
- 240000002853 Nelumbo nucifera Species 0.000 claims description 5
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 5
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 5
- DMBHHRLKUKUOEG-UHFFFAOYSA-N diphenylamine Chemical compound C=1C=CC=CC=1NC1=CC=CC=C1 DMBHHRLKUKUOEG-UHFFFAOYSA-N 0.000 claims description 3
- 230000007935 neutral effect Effects 0.000 claims description 3
- 238000001179 sorption measurement Methods 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 3
- 230000001360 synchronised effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 7
- 238000010521 absorption reaction Methods 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000000047 product Substances 0.000 description 2
- MYUPPPVXZSAJQS-VQTJNVASSA-N (4aR,9bS)-2,6-diacetyl-4a-ethoxy-1,7,9-trihydroxy-8,9b-dimethyl-4H-dibenzofuran-3-one Chemical compound CCO[C@@]12CC(=O)C(C(C)=O)=C(O)[C@]1(C)c1c(O2)c(C(C)=O)c(O)c(C)c1O MYUPPPVXZSAJQS-VQTJNVASSA-N 0.000 description 1
- 208000032368 Device malfunction Diseases 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/088—Aspects of digital computing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/086—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种供用电系统谐波源识别方法及模型,包括:同步实时采集供用电系统电源和所有负荷的电压和电流全信息,并标注时标存储;根据所采集的带时标的信息,计算该时刻的基波和各次谐波的电压、电流和阻抗的幅值和相位;计算和确定供用电系统相应时刻的电源和各用户的基波有功功率及谐波有功功率的大小和方向;根据供用电系统相应时刻的电源和各负荷的谐波功率方向识别谐波源;根据功率平衡原则,验证谐波源识别的正确性。本发明用于供用电系统谐波源识别、计算和分析谐波源对供用电系统产生的附加功率损耗,为供用电系统的谐波源监测、谐波功率计量和谐波污染责任的划分提供理论依据和有效技术手段,实时、客观、简单、实用和易推广。
Description
技术领域
本发明涉及供用电系统电能质量技术领域,特别涉及一种供用电系统谐波源识别方法及模型。
背景技术
随着电网规模的不断扩大,越来越多的新型用电设备和高精度用电设备接入到电力系统,向电网注入大量的谐波功率,使得输、供和用电产生额外的电能损耗,影响电气设备以及邻近的通讯设备的正常工作,同时也使得电力系统的继电保护装置误动作,可靠性降低。因此,供用电系统中进行谐波源识别,能够提高电网运行的可靠性、灵活性、安全性和经济性。
电网采用的电能计量仪表大多数是基于基波进行设计和计量,在计量时并没有考虑到电网存在谐波的情况。但由于谐波源的存在,将会直接影响电网和用电户,出现一些电能质量问题。本发明根据在供用电系统中的谐波源特点,提出了供用电系统谐波源识别方法及模型,为解决供用电系统中准确识别谐波源问题提供理论依据和有效技术手段。
发明内容
本发明提供一种供用电系统谐波源识别方法及模型,旨在解决供用电系统中,由于谐波源的存在,其产生的谐波电流会反馈到电网和电力用电户,直接导致电网和电力用电户的能量质量降低以及电能计量不公平的问题。
本发明的目的是这样实现的:
供用电系统谐波源识别方法,其特征在于,包括下列步骤:
步骤1、同步实时采集供用电系统电压源和所有负荷的电压和电流全信息,并标注时标后存储,具体是:分别通过电压互感器、电流互感器、高速高精度模拟量采集器和全球卫星定位系统同步实时采集供用电系统电压源和所有负荷的电压和电流全信息;当同步实时采集供用电系统电压源和所有负荷的电压和电流全信息时,其最低采样速率和采样精度满足谐波测量的国家标准和行业标准;同时,设置北斗定位授时系统,保证规定时间和精度内各采样点采样的同步性和实时性;
步骤2、根据步骤1所采集的带时标的电压和电流全信息,计算相应时刻的基波和各次谐波的电压、电流和阻抗的幅值和相位,基于以下公式:基波阻抗如式(1)所示;
谐波阻抗如式(2)所示;
Zk为谐波阻抗,Rk为谐波电阻,Xk为谐波电抗,k为谐波次数,nmax为最大谐波次数;
步骤3、计算和确定供用电系统相应时刻的电源和各用户的基波有功功率及谐波有功功率的大小和方向,其中,
供用电系统相应时刻的电源和各用户的基波和谐波的有功功率、无功功率和功率因数的大小计算如式(3)所示;
供用电系统相应时刻含各次谐波的总电压、电流、有功功率、无功功率和功率因数的大小计算如式(4)所示;
步骤4、根据步骤3得到的供用电系统相应时刻的电源和各负荷的谐波功率方向识别谐波源,具体是:
当电源的谐波功率为负,该电源产生谐波,即:被识别为谐波源电源或非线性电源;
当负荷的谐波功率为负时,该负荷产生谐波,即:被识别为谐波源负荷或非线性负荷;
步骤5、根据功率平衡原则,验证谐波源识别的正确性,即:
条件一:所有负荷和电网电源的各次对应谐波功率之和为0;
条件二:所有负荷和电网电源的基波功率之和也为0;
同时满足条件一和条件二的所有负功率,均为谐波源,即:谐波源的识别正确,否则,谐波源的识别有误。
在上述的供用电系统谐波源识别方法,所述步骤1中,三相供用电系统的电压互感器和电流互感器采用Y-N型接线。三个电压互感器的一次部分的同名端分别接A、B和C相电压端,异名端共同接入中性点N端;二次部分的各同名端分别定义为A、B和C相电压的正极性、各异名端分别定义为A、B和C相电压的负极性。三个电流互感器的一次部分分别串联接入A、B和C三相线路,各同名端分别接入A、B和C相电源侧、各异名端分别接入A、B和C相负荷侧;二次部分的各同名端分别定义为A、B和C相电流的正极性、各异名端分别定义为A、B和C相电流的负极性。
在上述的供用电系统谐波源识别方法,所述步骤2中,以时标为检索信息,调用的相应的电压和电流采样值,采样快速傅里叶算法,计算该时刻电压和电流的基波和各次谐波的电压和电流的幅值和相位,定义被采样的模拟信号是一个周期性时间函数x(t),除基波外还含有不衰减的直流分量和各次谐波,表示为:
式中:
式中n为自然数,n=0,1,2,……;an和bn则分别为各次谐波的正弦项和余弦项的振幅,a0为直流分量的值;
x(t)中的各次谐波分量xn(t)=ancos nω1t+bnsin nω1t的实部、虚部幅值an,bn分别为:
各次谐波分量xn(t)的有效值Xn和相角θn计算,各次谐波分量xn(t)的有效值Xn和相角θn为:
各次谐波分量xn(t)的有效值Xn和相角θn与电压和电流的对应关系为:
Xn→Vn,In;
(功率因数角或阻抗角)
然后根据基波和各次谐波的电压和电流的幅值和相位,计算基波和各次谐波阻抗的幅值和相位。
在上述的供用电系统谐波源识别方法,所述步骤3中,根据负荷功率定义各节点的电压和电流参考方向,即电压和电流同方向;电源基波功率因数角绝对值小于90°时,电源基波有功功率为负;当电源谐波功率因数角绝对值小于90°时,电源的谐波功率为负;当电源谐波功率因数角绝对值大于90°时,该电源的谐波功率为正,该电源吸收谐波;当负荷谐波功率因数角绝对值小于90°时,该负荷的谐波功率为正,该负荷吸收谐波;当负荷谐波功率因数角绝对值大于90°时,负荷的谐波功率为负,该负荷产生谐波。
在上述的供用电系统谐波源识别方法,所述步骤5具体方法包括:在基波电源作用下,各基波有功功率之和计算如公式(5)所示;
l为供用电系统的支路数,I1i为第i条支路的基波电流有效值,基波功率为负值时,取“-”;基波功率为正值时,取“+”;
在谐波源作用下,各次谐波有功功率之和计算如公式(6)所示;
k为谐波次数,l为供用电系统的支路数,Iki为第i条支路第k次谐波电流的有效值,R1i为第i条支路的基波电阻,nmax为最大谐波次数;当谐波功率为负值时,取“-”;当谐波功率为正值时,取“+”;
非线性负荷谐波源所发出的能量来自于系统电压源所发出能量的吸收转化;非线性负荷会消耗系统电压源的基波功率,反馈给系统和线性负荷谐波功率,系统和线性负荷被迫吸收谐波功率;系统和线性负荷吸收的谐波功率来自于系统的各个谐波源;所以,根据功率平衡的原则,所有负荷和电网电源消耗的各次对应谐波功率之和为0,即:公式(6)等于零;同时,消耗的基波功率之和也为0,即:公式(5)也等于零;满足这些条件的所有负功率,均为谐波源,即:谐波源的识别正确;若公式(6)不等于零,则则返回至步骤1重新计算。
一种供用电系统谐波源识别模型,包括系统电压源us和系统内阻抗Zs串联支路、线性负荷阻抗支路和非线性负荷谐波电流源与非线性负荷等效阻抗并联支路构成的按系统支路、线性负荷支路和非线性支路顺序并联的供用电系统谐波源识别模型。
线性负荷用支路阻抗表示,并且线性支路为第1条支路到第i条支路;将非线性负荷用谐波电流源和支路阻抗的并联来表示,并且,非线性负荷支路为第i+1条支路到第m条支路。因此,线性负荷支路和非线性负荷支路共同组成m条支路,系统电压源和系统阻抗支路向线性负荷支路和非线性负荷支路供电。非线性负荷支路的谐波电流源可以通过电流检测和傅里叶分解获得,它对于各次谐波电流之和;线性和非线性负荷支路的支路阻抗可以通过检测的负荷电压和负荷电流计算获得,负荷支路的基波支路阻抗对于基波电压与基波电流的比值;负荷支路的某次谐波支路阻抗对于该次谐波电压与该次谐波电流的比值。
非线性负荷将谐波电流源和支路阻抗并联支路可以变换为谐波电压源和支路阻抗串联支路来表示,谐波电压源等于谐波电流源与其支路阻抗的乘积,原支路阻抗值不变。
按系统电压源和各次谐波电压源分别单独作用,计算各支路的基波和各次谐波的电流、电压、有功功率、无功功率和功率因数。
负荷支路为第1条支路到第m条支路,其中:线性负荷i条支路,非线性负荷m-i条支路。先计算系统电压源单独作用在各负荷支路产生基波电流Ii和基波阻抗(电阻Ri和电抗Xi),根据计算公式计算出各负荷支路的基波功率的大小;利用负荷基波电流与监测点电压的相位关系,得到负荷基波功率的方向。再计算各次谐波源单独作用在各负荷支路产生谐波电流Ii(k)和谐波阻抗(电阻Ri(k)和电抗Xi(k)),根据计算公式计算出各支路负荷的k次谐波功率的大小;利用负荷k次谐波电流与监测点电压的相位关系,得到负荷k次谐波功率的方向。当支路谐波功率小于零时,该支路为谐波源支路。
因此,本发明的特点是:1、可以采用电能计量系统识别供用电系统谐波源,无需另建数据采集系统;2、计算和分析谐波源对供用电系统产生的附加功率损耗;3、谐波源识别模型简单、实用和易推广。
附图说明
图1是本发明供用电系统谐波源识别方法实施例流程示意图。
图2是本发明供用电系统谐波电流源识别模型示意图。
图3是本发明供用电系统谐波电压源识别模型示意图。
图4是本发明谐波电压源单独作用电路示意图。
图5是本发明系统电压源单独作用电路示意图。
图6是本发明线性负荷和非线性负荷网络模型简化示意图。
具体实施方式
下面结合附图和实施例详细说明:
本发明提供供用电系统谐波源识别方法,流程图如图1所示,具体包括:
步骤1、同步实时采集供用电系统电压源和所有负荷的电压和电流全信息,并标注时标后存储,具体是:分别通过电压互感器、电流互感器、高速高精度模拟量采集器和全球卫星定位系统同步实时采集供用电系统电压源和所有负荷的电压和电流全信息;当同步实时采集供用电系统电压源和所有负荷的电压和电流全信息时,其最低采样速率和采样精度满足谐波测量的国家标准和行业标准;同时,设置北斗定位授时系统,保证规定时间和精度内各采样点采样的同步性和实时性;
步骤2、根据步骤1所采集的带时标的电压和电流全信息,计算相应时刻的基波和各次谐波的电压、电流和阻抗的幅值和相位,基于以下公式:基波阻抗如式(1)所示;
谐波阻抗如式(2)所示;
Zk为谐波阻抗,Rk为谐波电阻,Xk为谐波电抗,k为谐波次数,nmax为最大谐波次数;
步骤3、计算和确定供用电系统相应时刻的电源和各用户的基波有功功率及谐波有功功率的大小和方向,其中,
供用电系统相应时刻的电源和各用户的基波有功功率及谐波有功功率、功率因数的大小计算如式(3)所示;
供用电系统相应时刻含各次谐波的总电压、电流、有功功率、无功功率和功率因数的大小计算如式(4)所示;
步骤4、根据步骤3得到的供用电系统相应时刻的电源和各负荷的谐波功率方向识别谐波源,具体是:
当电源的谐波功率为负,该电源产生谐波,即:被识别为谐波源电源或非线性电源;
当负荷的谐波功率为负时,该负荷产生谐波,即:被识别为谐波源负荷或非线性负荷;
步骤5、根据功率平衡原则,验证谐波源识别的正确性,即:
条件一:所有负荷和电网电源的各次对应谐波功率之和为0;
条件二:所有负荷和电网电源的基波功率之和也为0;
同时满足条件一和条件二的所有负功率,均为谐波源,即:谐波源的识别正确,否则,谐波源的识别错误。
其中,
所述步骤1中,三相供用电系统的电压互感器和电流互感器采用Y-N型接线。三个电压互感器的一次部分的同名端分别接A、B和C相电压端,异名端共同接入中性点N端;二次部分的各同名端分别定义为A、B和C相电压的正极性、各异名端分别定义为A、B和C相电压的负极性。三个电流互感器的一次部分分别串联接入A、B和C三相线路,各同名端分别接入A、B和C相电源侧、各异名端分别接入A、B和C相负荷侧;二次部分的各同名端分别定义为A、B和C相电流的正极性、各异名端分别定义为A、B和C相电流的负极性。
步骤2中,以时标为检索信息,调用的相应的电压和电流采样值,采样快速傅里叶算法,计算该时刻电压和电流的基波和各次谐波的电压和电流的幅值和相位,定义被采样的模拟信号是一个周期性时间函数x(t),除基波外还含有不衰减的直流分量和各次谐波,表示为:
式中:
式中n为自然数,n=0,1,2,……;an和bn则分别为各次谐波的正弦项和余弦项的振幅,a0为直流分量的值;
x(t)中的各次谐波分量xn(t)=ancos nω1t+bnsin nω1t的实部、虚部幅值an,bn分别为:
各次谐波分量xn(t)的有效值Xn和相角θn计算,各次谐波分量xn(t)的有效值Xn和相角θn为:
各次谐波分量xn(t)的有效值Xn和相角θn与电压和电流的对应关系为:
Xn→Vn,In;
(功率因数角或阻抗角)
然后根据基波和各次谐波的电压和电流的幅值和相位,计算基波和各次谐波阻抗的幅值和相位。
其中,步骤3中,根据负荷功率定义各节点的电压和电流参考方向,即电压和电流同方向;电源基波功率因数角绝对值小于90°时,电源基波有功功率为负;当电源谐波功率因数角绝对值小于90°时,电源的谐波功率为负;当电源谐波功率因数角绝对值大于90°时,该电源的谐波功率为正,该电源吸收谐波;当负荷谐波功率因数角绝对值小于90°时,该负荷的谐波功率为正,该负荷吸收谐波;当负荷谐波功率因数角绝对值大于90°时,负荷的谐波功率为负,该负荷产生谐波。
其中,步骤5具体包括:
在基波电源作用下,各基波有功功率之和计算如公式(5)所示;
l为供用电系统的支路数,I1i为第i条支路的基波电流有效值,基波功率为负值时,取“-”;基波功率为正值时,取“+”;
在谐波源作用下,各次谐波有功功率之和计算如公式(6)所示;
k为谐波次数,l为供用电系统的支路数,Iki为第i条支路第k次谐波电流的有效值,R1i为第i条支路的基波电阻,nmax为最大谐波次数;当谐波功率为负值时,取“-”;当谐波功率为正值时,取“+”;
非线性负荷谐波源所发出的能量来自于系统电压源所发出能量的吸收转化;非线性负荷会消耗系统电压源的基波功率,反馈给系统和线性负荷谐波功率,系统和线性负荷被迫吸收谐波功率;系统和线性负荷吸收的谐波功率来自于系统的各个谐波源;所以,根据功率平衡的原则,所有负荷和电网电源消耗的各次对应谐波功率之和为0,即:公式(6)等于零;同时,消耗的基波功率之和也为0,即:公式(5)也等于零;满足这些条件的所有负功率,均为谐波源,即:谐波源的识别正确;若公式(6)不等于零,则需要重新计算。
具体实施案例:
如图2所示,本发明供用电系统谐波电流源识别模型示意图,由系统电压源为网络模型提供基波电能。其中,将线性负荷用支路阻抗表示,并且线性支路为第1条支路到第i条支路;将非线性负荷用谐波电流源和支路阻抗的并联来表示,并且,非线性负荷支路为第i+1条支路到第m条支路。因此,线性负荷支路和非线性负荷支路共同组成m条支路,系统电压源和系统阻抗支路向线性负荷支路和非线性负荷支路供电。非线性负荷支路的谐波电流源可以通过电流检测和傅里叶分解获得,它对于各次谐波电流之和;线性和非线性负荷支路的支路阻抗可以通过检测的负荷电压和负荷电流计算获得,负荷支路的基波支路阻抗对于基波电压与基波电流的比值;负荷支路的某次谐波支路阻抗对于该次谐波电压与该次谐波电流的比值。
图3所示,本发明供用电系统谐波电压源识别模型示意图,图2中非线性负荷将谐波电流源和支路阻抗并联支路变换为谐波电压源和支路阻抗串联支路来表示,谐波电压源等于谐波电流源与其支路阻抗的乘积。
图4所示,将图3电路中谐波源单独作用时提取出来,非线性负荷支路为第i+1条支路到第m条支路,共m-i条非线性支路上的k次谐波源共同作用,在电网中产生k次谐波电流,再根据计算公式即可计算出各支路负荷的k次谐波功率的大小;利用负荷k次谐波电流与监测点电压的相位关系,得到负荷k次谐波功率的方向。
图5所示,将图3电路中系统电压源单独作用时提取出来,负荷支路为第1条支路到第m条支路,共m条支路。系统电压源在电网中产生基波电流,再根据计算公式即可计算出各负荷支路的基波功率的大小;利用负荷基波电流与监测点电压的相位关系,得到负荷基波功率的方向。
图6所示,将线性负荷和非线性负荷的电路网络模型进行简化,其中us=100cos(2πft-90°),f=50Hz,ω0=2πf,Z1=0.5Ω,Z2=10Ω,Z3=15Ω,R4=1Ω,L=0.01H。将Us的一个周期内进行12倍基波频率的采样,则最高谐波次数为5次。通过快速傅里叶变换,参考电压uab为:
uab=7.53+80.44cos(ω0t-89.4°)+5.32cos(2ω0t-15.5°)+0.62cos(3ω0t-176°)+1.13cos(4ω0t-14.3°)+0.61cos(5ω0t-164°);
流经Z1支路的电流为:
iZ1=15.06+39.2cos(ω0t-89.4°)+10.63cos(2ω0t-15.5°)+1.24cos(3ω0t-176°)+2.25cos(4ω0t-14.3°)+1.23cos(5ω0t-164°);
流经Z2支路的电流为:
iZ2=0.75+8.04cos(ω0t-89.4°)+0.53cos(2ω0t-15.5°)+0.06cos(3ω0t-176°)+0.11cos(4ω0t-14.3°)+0.06cos(5ω0t-164°);
流经Z3支路的电流为:
iZ3=0.5+5.36cos(ω0t-89.4°)+0.35cos(2ω0t-15.5°)+0.04cos(3ω0t-176°)+0.08cos(4ω0t-14.3°)+0.04cos(5ω0t-164°);流经Z4支路的电流为:
iZ4=16.31+25.85cos(ω0t-95.7°)+11.52cos(2ω0t+164.5°)+1.34cos(3ω0t+4°)+2.44cos(4ω0t+165.7°)+1.33cos(5ω0t+16°)。对比,流经Z1,Z2,Z3,Z4所在支路电流与参考电压uab,可以看出,在基波下,uab相位与流经Z1的电流相位之差则为电阻负荷,高次谐波下,则相位同相,说明Z1为线性电阻负荷;在基波和高次谐波下,uab相位与流经Z2和Z3的电流相位之差为零,则相位同相,说明Z2和Z3均为线性电阻负荷;在基波下,uab相位与流经Z4的电流相位之差为6.3°∈(-90°,-90°),则为感性负荷,在高次谐波下,uab相位与流经Z4的电流相位之差为-180°∈[-180°,-90°]∪[90°,180°],相位均反相,说明Z4为非线性负荷,即谐波源。
根据功率平衡的原则,来验证谐波源识别的正确性,根据功率计算式:其中,u0,i0为谐波源导致产生的直流分量,Uk,Ik为k次谐波有效值。则
PZ1=113.36+384.21+28.25+0.39+1.27+0.38=527.85(W);PZ2=5.67+323.56+1.41+0.02+0.06+0.02=330.74(W);PZ3=3.78+215.7+0.94+0.01+0.04+0.01=220.5(W);PZ4=-122.81+1033.02-30.6-0.42-1.37-0.41=877.41(W)。所有负荷和电网电源消耗的各次对应谐波功率之和为0,同时,消耗的基波功率之和也为0,并且,可以看出Z1,Z2,Z3所消耗的有功功率均为正值,所以Z1,Z2,Z3均为线性负荷;而Z4除基波下所消耗的有功功率为正,在高次谐波下,吸收的有功功率均为负。所以说明利用功率平衡的验证是正确的。本文中所描述的具体实施例仅仅是对本发明精神作举例说明。
本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (5)
1.一种供用电系统谐波源识别方法,其特征在于,基于一种供用电系统谐波源识别模型,该模型包括系统电压源us和系统内阻抗Zs串联支路、线性负荷阻抗支路和非线性负荷谐波电流源与非线性负荷等效阻抗并联支路构成的按系统支路、线性负荷支路和非线性支路顺序并联的供用电系统谐波源识别模型;
线性负荷用支路阻抗表示,并且线性支路为第1条支路到第i条支路;将非线性负荷用谐波电流源和支路阻抗的并联来表示,并且,非线性负荷支路为第i+1条支路到第m条支路;因此,线性负荷支路和非线性负荷支路共同组成m条支路,系统电压源和系统阻抗支路向线性负荷支路和非线性负荷支路供电;非线性负荷支路的谐波电流源可以通过电流检测和傅里叶分解获得,它对于各次谐波电流之和;线性和非线性负荷支路的支路阻抗可以通过检测的负荷电压和负荷电流计算获得,负荷支路的基波支路阻抗对于基波电压与基波电流的比值;负荷支路的某次谐波支路阻抗对于该次谐波电压与该次谐波电流的比值;
非线性负荷将谐波电流源和支路阻抗并联支路可以变换为谐波电压源和支路阻抗串联支路来表示,谐波电压源等于谐波电流源与其支路阻抗的乘积,原支路阻抗值不变;
按系统电压源和各次谐波电压源分别单独作用,计算各支路的基波和各次谐波的电流、电压、有功功率、无功功率和功率因数;
负荷支路为第1条支路到第m条支路,其中:线性负荷i条支路,非线性负荷m-i条支路;先计算系统电压源单独作用在各负荷支路产生基波电流Ii和基波阻抗,根据计算公式计算出各负荷支路的基波功率的大小;利用负荷基波电流与监测点电压的相位关系,得到负荷基波功率的方向;再计算各次谐波源单独作用在各负荷支路产生谐波电流Ii(k)和谐波阻抗,根据计算公式计算出各支路负荷的k次谐波功率的大小;利用负荷k次谐波电流与监测点电压的相位关系,得到负荷k次谐波功率的方向;当支路谐波功率小于零时,该支路为谐波源支路;
识别方法包括下列步骤:
步骤1、同步实时采集供用电系统电源和所有负荷的电压和电流全信息,并标注时标后存储,具体是:分别通过电压互感器、电流互感器、高速高精度模拟量采集器和全球卫星定位系统同步实时采集供用电系统电源和所有负荷的电压和电流全信息;当同步实时采集供用电系统电源和所有负荷的电压和电流全信息时,其最低采样速率和采样精度满足谐波测量的国家标准和行业标准;同时,设置北斗定位授时系统,保证规定时间和精度内各采样点采样的同步性和实时性;
步骤2、根据步骤1所采集的带时标的电压和电流全信息,计算相应时刻的基波和各次谐波的电压、电流和阻抗的幅值和相位,基于以下公式:基波阻抗如式(1)所示;
谐波阻抗如式(2)所示;
k=2,3,4,...,nmax
Zk为谐波阻抗,Rk为谐波电阻,Xk为谐波电抗,k为谐波次数,nmax为最大谐波次数;
步骤3、计算和确定供用电系统相应时刻的电源和各用户的基波有功功率及谐波有功功率的大小和方向,其中,
供用电系统相应时刻的电源和各用户的基波和谐波的有功功率、无功功率和功率因数的大小计算如式(3)所示;
供用电系统相应时刻含各次谐波的总电压、电流、有功功率、无功功率和功率因数的大小计算如式(4)所示;
步骤4、根据步骤3得到的供用电系统相应时刻的电源和各负荷的谐波功率方向识别谐波源,具体是:
当电源的谐波功率为负,该电源产生谐波,即:被识别为谐波源电源或非线性电源;
当负荷的谐波功率为负时,该负荷产生谐波,即:被识别为谐波源负荷或非线性负荷;
步骤5、根据功率平衡原则,验证谐波源识别的正确性,即:
条件一:所有负荷和电网电源的各次对应谐波功率之和为0;
条件二:所有负荷和电网电源的基波功率之和也为0;
同时满足条件一和条件二的所有负功率,均为谐波源,即:谐波源的识别正确,否则,谐波源的识别有误。
2.根据权利要求1所述的一种供用电系统谐波源识别方法,其特征在于,所述步骤1中,三相供用电系统的电压互感器和电流互感器采用Y-N型接线;三个电压互感器的一次部分的同名端分别接A、B和C相电压端,异名端共同接入中性点N端;二次部分的各同名端分别定义为A、B和C相电压的正极性、各异名端分别定义为A、B和C相电压的负极性;三个电流互感器的一次部分分别串联接入A、B和C三相线路,各同名端分别接入A、B和C相电源侧、各异名端分别接入A、B和C相负荷侧;二次部分的各同名端分别定义为A、B和C相电流的正极性、各异名端分别定义为A、B和C相电流的负极性。
3.根据权利要求1所述的一种供用电系统谐波源识别方法,其特征在于,所述步骤2中,以时标为检索信息,调用的相应的电压和电流采样值,采样快速傅里叶算法,计算该时刻电压和电流的基波和各次谐波的电压和电流的幅值和相位,定义被采样的模拟信号是一个周期性时间函数x(t),除基波外还含有不衰减的直流分量和各次谐波,表示为:
式中:
式中n为自然数,n=1,2,……;an和bn则分别为各次谐波的正弦项和余弦项的振幅,a0为直流分量的值;
x(t)中的各次谐波分量xn(t)=ancos nω1t+bnsin nω1t的实部、虚部幅值an,bn分别为:
各次谐波分量xn(t)的有效值Xn和相角θn计算,各次谐波分量xn(t)的有效值Xn和相角θn为:
各次谐波分量xn(t)的有效值Xn和相角θn与电压和电流的对应关系为:
Xn→Vn,In;
θn→θVn,θIn
φn=θVn-θIn
然后根据基波和各次谐波的电压和电流的幅值和相位,计算基波和各次谐波阻抗的幅值和相位。
4.根据权利要求1所述的一种供用电系统谐波源识别方法,其特征在于,所述步骤3中,根据负荷功率定义各节点的电压和电流参考方向,即电压和电流同方向;电源基波功率因数角绝对值小于90°时,电源基波有功功率为负;当电源谐波功率因数角绝对值小于90°时,电源的谐波功率为负;当电源谐波功率因数角绝对值大于90°时,该电源的谐波功率为正,该电源吸收谐波;当负荷谐波功率因数角绝对值小于90°时,该负荷的谐波功率为正,该负荷吸收谐波;当负荷谐波功率因数角绝对值大于90°时,负荷的谐波功率为负,该负荷产生谐波。
5.根据权利要求1所述的一种供用电系统谐波源识别方法,其特征在于,所述步骤5具体方法包括:在基波电源作用下,各基波有功功率之和计算如公式(5)所示;
l为供用电系统的支路数,I1i为第i条支路的基波电流有效值,基波功率为负值时,取“-”;基波功率为正值时,取“+”;
在谐波源作用下,各次谐波有功功率之和计算如公式(6)所示;
k为谐波次数,l为供用电系统的支路数,Iki为第i条支路第k次谐波电流的有效值,R1i为第i条支路的基波电阻,nmax为最大谐波次数;当谐波功率为负值时,取“-”;当谐波功率为正值时,取“+”;
非线性负荷谐波源所发出的能量来自于系统电源所发出能量的吸收转化;非线性负荷会消耗系统电源的基波功率,反馈给系统和线性负荷谐波功率,系统和线性负荷被迫吸收谐波功率;系统和线性负荷吸收的谐波功率来自于系统的各个谐波源;所以,根据功率平衡的原则,所有负荷和电网电源消耗的各次对应谐波功率之和为0,即:公式(6)等于零;同时,消耗的基波功率之和也为0,即:公式(5)也等于零;满足这些条件的所有负功率,均为谐波源,即:谐波源的识别正确;若公式(6)不等于零,则返回至步骤1重新计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610946158.7A CN106383295B (zh) | 2016-10-26 | 2016-10-26 | 一种供用电系统谐波源识别方法及模型 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610946158.7A CN106383295B (zh) | 2016-10-26 | 2016-10-26 | 一种供用电系统谐波源识别方法及模型 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106383295A CN106383295A (zh) | 2017-02-08 |
CN106383295B true CN106383295B (zh) | 2018-06-19 |
Family
ID=57957100
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610946158.7A Active CN106383295B (zh) | 2016-10-26 | 2016-10-26 | 一种供用电系统谐波源识别方法及模型 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106383295B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107238756B (zh) * | 2017-06-14 | 2019-10-25 | 扬州万泰电子科技有限公司 | 一种智能配变终端冲击性负荷谐波源识别方法 |
CN107144732B (zh) * | 2017-06-29 | 2019-08-30 | 西安科技大学 | 基于集总功率的用户侧谐波源定位方法 |
CN109100568A (zh) * | 2018-06-25 | 2018-12-28 | 南方电网科学研究院有限责任公司 | 电网谐波普测分析方法、装置、设备及介质 |
CN109061300B (zh) * | 2018-09-19 | 2020-11-10 | 江苏中凌高科技股份有限公司 | 一种pcc点的特征谐波源定位方法 |
CN110927426B (zh) * | 2018-09-20 | 2021-09-07 | 中国石油化工股份有限公司 | 一种配电网的谐波源区段定位方法 |
CN110927425B (zh) * | 2018-09-20 | 2021-09-07 | 中国石油化工股份有限公司 | 一种基于谐波电流监测的配电网的谐波源区段定位方法 |
CN109301827B (zh) * | 2018-10-10 | 2022-05-31 | 国电南瑞科技股份有限公司 | 基于谐波源自动辨识及分级治理的谐波控制方法和系统 |
CN109342815B (zh) * | 2018-10-29 | 2022-06-21 | 中国电力科学研究院有限公司 | 一种变电站谐波源定位方法与系统 |
CN110007147A (zh) * | 2018-11-29 | 2019-07-12 | 广西电网有限责任公司电力科学研究院 | 一种单母线多用户的谐波责任计算方法 |
CN109870615B (zh) * | 2019-03-12 | 2020-05-26 | 西安科技大学 | 适用于工矿企业的谐波责任评定方法 |
CN112904083A (zh) * | 2021-01-29 | 2021-06-04 | 许继集团有限公司 | 一种基于宽频测量的功率平衡校验方法及装置 |
CN117216721B (zh) * | 2023-11-09 | 2024-02-02 | 安徽安大清能电气科技有限公司 | 一种基于多判据融合的谐振状态及谐振路径识别方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7746084B2 (en) * | 2006-10-05 | 2010-06-29 | Newport News Shipbuilding And Dry Dock Company | Impedance measurement using line-to-line current injection |
CN101227092B (zh) * | 2007-12-28 | 2010-10-06 | 华北电力大学 | 一种配电网络中的谐波源识别方法 |
CN102323494B (zh) * | 2011-05-24 | 2013-05-15 | 山东电力研究院 | 一种多谐波源辨识方法 |
CN104280612B (zh) * | 2014-10-28 | 2017-05-24 | 清华大学 | 一种基于单频电流传输特性的分布式谐波源辨识方法 |
CN104375026B (zh) * | 2014-11-04 | 2017-12-01 | 国家电网公司 | 一种基于特征量分析辨识谐波源类型的方法 |
-
2016
- 2016-10-26 CN CN201610946158.7A patent/CN106383295B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106383295A (zh) | 2017-02-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106383295B (zh) | 一种供用电系统谐波源识别方法及模型 | |
CN104793082A (zh) | 基于谐波相关分析的电力系统谐波源识别装置 | |
WO2012067747A2 (en) | Data alignment in large scale electrical system applications | |
CN106501612B (zh) | 基于非均匀传输线的双回部分同塔线路零序参数测量方法 | |
CN103901273B (zh) | 电网谐波检测方法和电网谐波检测装置 | |
CN102096019A (zh) | 小电流接地系统单相接地故障测距方法和装置 | |
CN106405470A (zh) | 电子式电流互感器谐波准确度整体检测方法和检测装置 | |
CN108918973A (zh) | 基于upqc的高精度阻抗测量装置及其切换控制方法 | |
CN106771901A (zh) | 基于定频投切电阻原理的直流系统绝缘监测方法 | |
CN107144732B (zh) | 基于集总功率的用户侧谐波源定位方法 | |
CN106646136B (zh) | 带并联阻尼调匝式消弧线圈单相接地故障测距方法及系统 | |
Herrera et al. | Harmonic disturbance identification in electrical systems with capacitor banks | |
Han et al. | A new assessment for the total harmonic contributions at the point of common coupling | |
CN109557494A (zh) | 一种基于电容式电压互感器的谐波测量装置和系统 | |
CN110133563B (zh) | 一种电流传感器抗直流性能检测方法及装置 | |
CN205594091U (zh) | 一种电流互感器误差测量装置 | |
CN110146780A (zh) | 中性点不接地柔性配电网系统铁磁谐振判别方法 | |
Dong et al. | Research and analysis of signal injection method for measuring capacitive current in distribution network | |
CN110516692A (zh) | 一种基于超平面聚类的pmu数据检测方法 | |
Rai et al. | Algorithms for Synchrophasor-Enabled Digital Relay in Differential Protection Schemes | |
Zhang et al. | Small-current grounding fault location method based on transient main resonance frequency analysis | |
Ren et al. | A novel self-powered smart current sensor for power equipment | |
CN107490721A (zh) | 一种基于NB‑IoT的智能电表电能计量装置 | |
Filianin | Identification of harmonic distortion sources in distribution systems using the discrete fourier transform on periods | |
Blazic et al. | A modified harmonic current vector method for harmonic contribution determination |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20180601 Address after: 430072 Wuhan University, Luojia mountain, Wuchang District, Wuhan, Hubei Applicant after: Wuhan University Address before: 430072 Wuhan University, Luojia mountain, Wuchang District, Wuhan, Hubei Applicant before: Wuhan University Applicant before: XIANGYANG POWER SUPPLY COMPANY OF STATE GRID HUBEI ELECTRONIC POWER COMPANY Applicant before: State Grid Corporation of China |
|
TA01 | Transfer of patent application right |