CN102420426B - 基于注入功率灵敏度的电网局部降损方法 - Google Patents

基于注入功率灵敏度的电网局部降损方法 Download PDF

Info

Publication number
CN102420426B
CN102420426B CN 201110393178 CN201110393178A CN102420426B CN 102420426 B CN102420426 B CN 102420426B CN 201110393178 CN201110393178 CN 201110393178 CN 201110393178 A CN201110393178 A CN 201110393178A CN 102420426 B CN102420426 B CN 102420426B
Authority
CN
China
Prior art keywords
line
power
node
partiald
output
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
CN 201110393178
Other languages
English (en)
Other versions
CN102420426A (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.)
State Grid Corp of China SGCC
North China Electric Power University
Electric Power Research Institute of State Grid Gansu Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
North China Electric Power University
Electric Power Research Institute of State Grid Gansu Electric Power Co Ltd
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 State Grid Corp of China SGCC, North China Electric Power University, Electric Power Research Institute of State Grid Gansu Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN 201110393178 priority Critical patent/CN102420426B/zh
Publication of CN102420426A publication Critical patent/CN102420426A/zh
Application granted granted Critical
Publication of CN102420426B publication Critical patent/CN102420426B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/30Reactive power compensation

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了电网管理控制技术领域中的一种基于注入功率灵敏度的局部降损方法。该方法包括:计算不同电压等级下各个线路的线损及线损率,并按照线损率排序;确定待降损的线路,设定降损后期望得到的线损率;计算待降损的线路的输送功率的变化量;计算待降损的线路的输送功率对各个第一类节点有功出力的灵敏度和各个第二类节点注入无功功率的灵敏度;其中,第一类节点为连接发电机组的节点,第二类节点为连接无功补偿设备的节点;利用反向等量配对原则调整发电机有功出力和节点的无功补偿量。本发明解决了现有的电力系统降损方法存在的问题。

Description

基于注入功率灵敏度的电网局部降损方法
技术领域
本发明属于电网管理控制技术领域,尤其涉及一种基于注入功率灵敏度的局部降损方法。
背景技术
在电力系统中,输电、变电、配电和营销各环节中都会产生电能损失。如何对线损进行管理,是供电企业成本管理的一项重要内容。作为衡量电网运行经济性的一项重要技术指标,线损率的降低对于提高电能传输效率、贯彻国家节能减排政策具有重要意义。
常用电力系统中传统的降损方法主要包括两类:一是根据理论计算的数学模型,通过改变线路、变压器的相关运行参数来实现降低损耗的目的,但是该方法主要应用于电网的设计领域;二是采用最优潮流法,以全网线损最小为目标函数,对电网的潮流分布进行优化。但此方法具有受初值影响大、收敛速度慢且不易得到最优解的缺陷。显然,这两种方法都不能满足电网实时经济调度中快速、有效降低线损率的要求,故不适宜在线制定降损策略。
发明内容
本发明的目的在于,提出一种基于注入功率灵敏度的电网局部降损方法,用以解决现有的电力系统降损方法存在的问题。
为实现上述目的,本发明提供的技术方案是,一种基于注入功率灵敏度的电网局部降损方法,其特征是所述方法包括:
步骤1:计算不同电压等级下各个线路的线损及线损率,并按照线损率排序;
步骤2:确定待降损的线路,设定降损后期望得到的线损率;
步骤3:计算待降损的线路的输送功率的变化量;
步骤4:计算待降损的线路的输送功率对各个第一类节点有功出力的灵敏度和各个第二类节点注入无功功率的灵敏度;
所述第一类节点为连接发电机组的节点;
所述第二类节点为连接无功补偿设备的节点;
步骤5:利用反向等量配对原则调整发电机有功出力和节点的无功补偿量。
所述计算不同电压等级下各个线路的线损具体利用公式
Figure BDA0000115122280000021
其中,PLOSS为线路(i,j)的线损,i和j分别为线路(i,j)两端的节点,Pi为线路(i,j)所输送的有功功率,Qi为线路(i,j)所输送的无功功率,Ui为线路(i,j)的节点i端的电压,R为线路(i,j)的电阻。
所述计算不同电压等级下线路的线损率具体利用公式
Figure BDA0000115122280000022
其中,r为线路(i,j)的线损率,Pi为线路(i,j)所输送的有功功率,PLOSS是线路(i,j)的线损,i和j分别为线路(i,j)两端的节点。
所述计算待降损的线路的输送功率的变化量具体是:
当待降损的线路满足
Figure BDA0000115122280000023
时,则仅计及待降损的线路的无功变化,待降损的线路的输送功率的变化量为 ΔQ i = Q i - Q i ′ = ( 1 - b a ) Q i ;
当待降损的线路满足
Figure BDA0000115122280000025
时,则仅计及待降损的线路的有功变化,待降损的线路的输送功率的变化量为 ΔP i = P i - P i ′ = ( 1 - b a ) P i ;
当待降损的线路满足
Figure BDA0000115122280000032
时,待降损的线路的输送功率的变化量为 ΔP i = P i - P i ′ = ( 1 - b a ) P i , ΔQ i = Q i - Q i ′ = ( 1 - b a ) Q i ;
其中,Pi为待降损的线路(i,j)所输送的有功功率,Qi为待降损的线路(i,j)所输送的无功功率,P′i为降损后线路(i,j)所输送的有功功率,Q′i为降损后线路(i,j)所输送的无功功率,a为线损率,b为降损后期望得到的线损率,D1为第一设定阈值,D2为第二设定阈值,i和j分别为线路(i,j)两端的节点。
所述第一设定阈值D1=100。
所述第二设定阈值D2=0.01。
所述计算待降损的线路的输送功率对各个第一类节点有功出力的灵敏度采用公式
S PG i = ∂ P i ∂ P G = ∂ P i ∂ δ i ∂ δ i ∂ P G + ∂ P i ∂ δ j ∂ δ j ∂ P G + ∂ P i ∂ U i U i · ∂ U i ∂ P G 1 U i + ∂ P i ∂ U i U j · ∂ U j ∂ P G 1 U j
其中,
Figure BDA0000115122280000036
为待降损的线路(i,j)的输送功率对各个第一类节点有功出力的灵敏度,Pi为待降损的线路(i,j)的输送的有功功率,PG为发电机有功出力,Ui为节点i的电压幅值,δi为节点i的电压相角,Uj为节点j的电压幅值,δj为节点j的电压相角。
所述计算待降损的线路的输送功率对各个第二类节点注入无功功率的灵敏度采用公式
S QN i = ∂ Q i ∂ Q N = ∂ Q i ∂ δ i ∂ δ i ∂ Q N + ∂ Q i ∂ δ j ∂ δ j ∂ Q N + ∂ Q i ∂ U i U i · ∂ U i ∂ Q N 1 U i + ∂ Q i ∂ U j U j · ∂ U j ∂ Q N 1 U j
其中,
Figure BDA0000115122280000038
为待降损的线路(i,j)的输送功率对各个第二类节点注入无功功率的灵敏度,Qi为待降损的线路(i,j)的输送的无功功率,QN为第二类节点的注入无功功率,Ui为节点i的电压幅值,δi为节点i的电压相角,Uj为节点j的电压幅值,δj为节点j的电压相角。
所述利用反向等量配对原则调整发电机有功出力具体包括下列步骤:
步骤101:当待降损的线路的输送功率对第一类节点有功出力的灵敏度大于零,即
Figure BDA0000115122280000041
时,将与该第一类节点对应的发电机组等效为减出力发电机组,并将所述减出力发电机组加入减出力发电机组集合GPLUS
当待降损的线路的输送功率对第一类节点有功出力的灵敏度小于零,即
Figure BDA0000115122280000042
时,将与该第一类节点对应的发电机组等效为加出力发电机组,并将所述加出力发电机组加入加出力发电机组集合GMINUS
当待降损的线路的输送功率对发电机节点有功出力的灵敏度等于零,即
Figure BDA0000115122280000043
时,将与该第一类节点对应的发电机组等效为平衡功率发电机组,并将所述平衡功率发电机组加入平衡功率发电机组集合GZERO
步骤102:按照有功出力的灵敏度的绝对值的大小,分别对减出力发电机组集合GPLUS中的各个减出力发电机组和加出力发电机组集合GMINUS中的各个加出力发电机组进行排序;排序后减出力发电机组集合中的元素记为
Figure BDA0000115122280000044
Figure BDA0000115122280000045
为排序后第k个减出力发电机组有功出力的灵敏度,且有
Figure BDA0000115122280000046
k=1,2,...,n-1,n为减出力发电机组的个数;排序后加出力发电机组集合记为
Figure BDA0000115122280000047
为排序后第l个加出力发电机组有功出力的灵敏度,且有
Figure BDA0000115122280000049
l=1,2,...,m-1,m为加出力发电机组的个数;
步骤103:令k=1且l=1;
步骤104:判断k>n是否成立,如果成立,则从平衡功率发电机组集合GZERO中取出一个未使用的元素,将其作为第k个减出力发电机组有功出力的灵敏度
Figure BDA0000115122280000051
并将该取出的未使用的元素标识为已使用元素;
步骤105:判断l>m是否成立,如果成立,则从平衡功率发电机组集合GZERO中取出一个未使用的元素,将其作为第l个加出力发电机组有功出力的灵敏度
Figure BDA0000115122280000052
并将该取出的未使用的元素标识为已使用元素;
步骤106:减少第k个减出力发电机组的有功出力,减少量为
Figure BDA0000115122280000053
增加第l个加出力发电机组的有功出力,增加量为
Figure BDA0000115122280000054
其中,
ΔP PLUS k = ΔP MINUS l = min { ΔP i S PG , PLUS k - S PG , MINUS l , ΔP PLUS , k max , ΔP MINUS , l max }
Figure BDA0000115122280000056
为第k个减出力发电机组有功出力与该发电机组有功出力的下限的差值;
Figure BDA0000115122280000057
为第l个加出力发电机组有功出力的上限与该发电机组有功出力的差值;
步骤107:判断待降损的线路是否达到降损后期望得到的线损率,或者是否达到减出力发电机组集合的调节下限,或者是否达到加出力发电机组集合的调节上限,如果达到,则执行步骤108;否则,令k=k+1且l=l+1,返回步骤104;
所述减出力发电机组集合的调节下限为
Figure BDA0000115122280000058
其中,
Figure BDA0000115122280000059
为减出力发电机组集合的调节下限,n为减出力发电机组的个数;
所述加出力发电机组集合的调节上限为
Figure BDA00001151222800000510
其中,为加出力发电机组集合的调节上限,m为加出力发电机组的个数;
步骤108:结束。
所述利用反向等量配对原则调整节点的无功补偿量具体包括下列步骤:
步骤201:当待降损的线路的输送功率对第二类节点注入无功功率的灵敏度大于零,即时,将其等效为减出力节点,并将其加入减出力节点集合QPLUS
当待降损的线路的输送功率对第二类节点注入无功功率的灵敏度小于零,即时,将其等效为加出力节点,并将其加入加出力节点集合QMINUS
当待降损的线路的输送功率对第二类节点注入无功功率的灵敏度等于零,即
Figure BDA0000115122280000063
时,将其等效为平衡功率节点,并将其加入平衡功率节点集合QZERO
步骤202:按照节点注入无功功率的灵敏度的绝对值的大小,分别对减出力节点集合GPLUS中的各个减出力节点和加出力节点集合GMINUS中的各个加出力节点进行排序;排序后减出力节点集合中的元素记为
Figure BDA0000115122280000064
为排序后第k个减出力节点注入无功功率的灵敏度,且有
Figure BDA0000115122280000065
k=1,2,...,n-1,n为减出力节点的个数;排序后加出力节点集合记为
Figure BDA0000115122280000066
Figure BDA0000115122280000067
为排序后第l个加出力节点注入无功功率的灵敏度,且有
Figure BDA0000115122280000068
l=1,2,...,m-1,m为加出力节点的个数;
步骤203:令k=1且l=1;
步骤204:判断k>n是否成立,如果成立,则从平衡功率节点集合QZERO中取出一个未使用的元素,将其作为第k个减出力节点注入无功功率的灵敏度
Figure BDA0000115122280000069
并将该取出的未使用的元素标识为已使用元素;
步骤205:判断l>m是否成立,如果成立,则从平衡功率节点集合QZERO中取出一个未使用的元素,将其作为第l个加出力节点注入无功功率的灵敏度
Figure BDA0000115122280000071
并将该取出的未使用的元素标识为已使用元素;
步骤206:减少第k个减出力节点的无功注入,减少量为
Figure BDA0000115122280000072
增加第l个加出力节点的无功注入,增加量为
Figure BDA0000115122280000073
其中,
ΔQ PLUS k = ΔQ MINUS l = min { ΔQ i S QN , PLUS k - S QN , MINUS l , ΔQ PLUS , k max , ΔQ MINUS , l max }
Figure BDA0000115122280000075
为第k个减出力节点注入无功功率与该节点无功补偿量的下限的差值;
Figure BDA0000115122280000076
为第l个加出力节点无功补偿量的上限与该节点注入无功功率的差值;
步骤207:判断待降损的线路是否达到降损后期望得到的线损率,或者是否达到减出力节点集合的调节下限,或者是否达到加出力节点集合的调节上限,如果达到,则执行步骤208;否则,令k=k+1且l=l+1,返回步骤204;
所述减出力节点集合的调节下限为
Figure BDA0000115122280000077
其中,
Figure BDA0000115122280000078
为减出力节点集合的调节下限,n为减出力节点的个数;
所述加出力节点集合的调节上限为
Figure BDA0000115122280000079
其中,
Figure BDA00001151222800000710
为加出力节点集合的调节上限,m为加出力节点的个数;
步骤208:结束。
本发明提供的方法,在保证全网线损不提高的前提下,基于电网实时数据,以线损率最小为目标函数,根据交流灵敏度法和反向等量配对原则调整发电机出力或节点无功补偿量,实现了对某一线损过高支路的降损,有效解决了现有的电力系统降损方法存在的问题。
附图说明
图1是基于注入功率灵敏度的电网局部降损方法流程图;
图2是IEEE 118节点系统接线图。
具体实施方式
下面结合附图,对优选实施例作详细说明。应该强调的是,下述说明仅仅是示例性的,而不是为了限制本发明的范围及其应用。
实施例1
图1是基于注入功率灵敏度的电网局部降损方法流程图,图1中,本发明提供的基于注入功率灵敏度的电网局部降损方法包括:
步骤1:计算不同电压等级下各个线路的线损及线损率,并按照线损率排序。
计算不同电压等级下各个线路的线损具体利用公式
Figure BDA0000115122280000081
其中,PLOSS为线路(i,j)的线损,i和j分别为线路(i,j)两端的节点,Pi为线路(i,j)所输送的有功功率,Qi为线路(i,j)所输送的无功功率,Ui为线路(i,j)的节点i端的电压,R为线路(i,j)的电阻。
计算不同电压等级下线路的线损率具体利用公式
Figure BDA0000115122280000082
其中,r为线路(i,j)的线损率,Pi为线路(i,j)所输送的有功功率,PLOSS是线路(i,j)的线损,i和j分别为线路(i,j)两端的节点。
步骤2:确定待降损的线路,设定降损后期望得到的线损率。
步骤3:计算待降损的线路的输送功率的变化量。具体是:
当待降损的线路满足
Figure BDA0000115122280000083
时,则仅计及待降损的线路的无功变化,待降损的线路的输送功率的变化量为
Figure BDA0000115122280000091
当待降损的线路满足
Figure BDA0000115122280000092
时,则仅计及待降损的线路的有功变化,待降损的线路的输送功率的变化量为
Figure BDA0000115122280000093
当待降损的线路满足
Figure BDA0000115122280000094
时,待降损的线路的输送功率的变化量为 ΔP i = P i - P i ′ = ( 1 - b a ) P i , ΔQ i = Q i - Q i ′ = ( 1 - b a ) Q i . 其中,Pi为待降损的线路(i,j)所输送的有功功率,Qi为待降损的线路(i,j)所输送的无功功率,P′i为降损后线路(i,j)所输送的有功功率,Q′i为降损后线路(i,j)所输送的无功功率,a为线损率,b为降损后期望得到的线损率。
步骤4:计算待降损的线路的输送功率对各个第一类节点有功出力的灵敏度和各个第二类节点注入无功功率的灵敏度。其中,第一类节点为连接发电机组的节点,第二类节点为连接无功补偿设备的节点。
计算待降损的线路的输送功率对各个第一类节点有功出力的灵敏度采用公式
S PG i = ∂ P i ∂ P G = ∂ P i ∂ δ i ∂ δ i ∂ P G + ∂ P i ∂ δ j ∂ δ j ∂ P G + ∂ P i ∂ U i U i · ∂ U i ∂ P G 1 U i + ∂ P i ∂ U i U j · ∂ U j ∂ P G 1 U j
其中,
Figure BDA0000115122280000098
为待降损的线路(i,j)的输送功率对各个第一类节点有功出力的灵敏度,Pi为待降损的线路(i,j)输送的有功功率,PG为发电机有功出力,Ui为节点i的电压幅值,δi为节点i的电压相角,Uj为节点j的电压幅值,δj为节点j的电压相角。
计算待降损的线路的输送功率对各个第二类节点注入无功功率的灵敏度采用公式
S QN i = ∂ Q i ∂ Q N = ∂ Q i ∂ δ i ∂ δ i ∂ Q N + ∂ Q i ∂ δ j ∂ δ j ∂ Q N + ∂ Q i ∂ U i U i · ∂ U i ∂ Q N 1 U i + ∂ Q i ∂ U j U j · ∂ U j ∂ Q N 1 U j
其中,
Figure BDA00001151222800000910
为待降损的线路(i,j)的输送功率对各个第二类节点注入无功功率的灵敏度,Qi为待降损的线路(i,j)输送的无功功率,QN为第二类节点的注入无功功率,Ui为节点i的电压幅值,δi为节点i的电压相角,Uj为节点j的电压幅值,δj为节点j的电压相角。
待降损的线路(i,j)输送的有功功率Pi对待降损的线路(i,j)两端电压相角与幅值的偏导数为:
∂ P i ∂ δ i = U i U j ( G ij sin δ ij - B ij cos δ ij ) ∂ P i ∂ δ j = - U i U j ( G ij sin δ ij - B ij cos δ ij ) ∂ P i ∂ U i U i = 2 U i 2 G ij - U i U j ( G ij cos δ ij + B ij sin δ ij ) ∂ P i ∂ U j U j = - U i U j ( G ij cos δ ij + B ij sin δ ij )
待降损的线路(i,j)输送的无功功率Qi对待降损的线路(i,j)两端电压相角与幅值的偏导数为:
∂ Q i ∂ δ i = - U i U j ( B ij sin δ ij + G ij cos δ ij ) ∂ Q i ∂ δ j = U i U j ( B ij sin δ ij + G ij cos δ ij ) ∂ Q i ∂ U i U i = - 2 U i 2 B ij + U i U j ( B ij cos δ ij - G ij sin δ ij ) ∂ Q i ∂ U j U j = U i U j ( B ij cos δ ij - G ij sin δ ij )
这样,由上式可从潮流计算中间结果雅克比矩阵J中获得
Figure BDA0000115122280000103
Figure BDA0000115122280000104
右端的四个偏导数的值。又由于雅可比矩阵逆阵J-1中的元素就是各节点电压的相角与幅值对各节点的P,Q导数,可以用J-1来计算
Figure BDA0000115122280000105
Figure BDA0000115122280000106
中的其他4个偏导数。
其中,G=(G1,G2,…Gm)为第一类节点,N=(N1,N2,…,Nn)为第二类节点,Pi、Qi为待降损的线路(i,j)输送的有功功率和无功功率,PG、QN分别为第一类节点有功出力和第二类节点注入无功功率,Ui、Uj、δi、δj分别为节点i和j电压的幅值与相角,Bij为待降损的线路(i,j)的电纳,Gij为待降损线路的电导、δij为节点i和j电压的相角差。
步骤5:利用反向等量配对原则调整发电机有功出力和节点的无功补偿量。
利用反向等量配对原则调整发电机有功出力具体包括下列步骤:
步骤101:当待降损的线路的输送功率对第一类节点有功出力的灵敏度大于零,即
Figure BDA0000115122280000111
时,将与该第一类节点对应的发电机组等效为减出力发电机组,并将所述减出力发电机组加入减出力发电机组集合GPLUS
当待降损的线路的输送功率对第一类节点有功出力的灵敏度小于零,即
Figure BDA0000115122280000112
时,将与该第一类节点对应的发电机组等效为加出力发电机组,并将所述加出力发电机组加入加出力发电机组集合GMINUS
当待降损的线路的输送功率对发电机节点有功出力的灵敏度等于零,即
Figure BDA0000115122280000113
时,将与该第一类节点对应的发电机组等效为平衡功率发电机组,并将所述平衡功率发电机组加入平衡功率发电机组集合GZERO
步骤102:按照有功出力的灵敏度的绝对值的大小,分别对减出力发电机组集合GPLUS中的各个减出力发电机组和加出力发电机组集合GMINUS中的各个加出力发电机组进行排序;排序后减出力发电机组集合中的元素记为
Figure BDA0000115122280000114
Figure BDA0000115122280000115
为排序后第k个减出力发电机组有功出力的灵敏度,且有
Figure BDA0000115122280000116
k=1,2,...,n-1,n为减出力发电机组的个数;排序后加出力发电机组集合记为
Figure BDA0000115122280000121
为排序后第l个加出力发电机组有功出力的灵敏度,且有
Figure BDA0000115122280000122
l=1,2,...,m-1,m为加出力发电机组的个数。
步骤103:令k=1且l=1。
步骤104:判断k>n是否成立,如果成立,则从平衡功率发电机组集合GZERO中取出一个未使用的元素,将其作为第k个减出力发电机组有功出力的灵敏度
Figure BDA0000115122280000123
并将该取出的未使用的元素标识为已使用元素。
步骤105:判断l>m是否成立,如果成立,则从平衡功率发电机组集合GZERO中取出一个未使用的元素,将其作为第l个加出力发电机组有功出力的灵敏度并将该取出的未使用的元素标识为已使用元素;。
步骤106:减少第k个减出力发电机组的有功出力,减少量为
Figure BDA0000115122280000125
增加第l个加出力发电机组的有功出力,增加量为
Figure BDA0000115122280000126
其中,
ΔP PLUS k = ΔP MINUS l = min { ΔP i S PG , PLUS k - S PG , MINUS l , ΔP PLUS , k max , ΔP MINUS , l max }
为第k个减出力发电机组有功出力与该发电机组有功出力的下限的差值;为第l个加出力发电机组有功出力的上限与该发电机组有功出力的差值。
步骤107:判断待降损的线路是否达到降损后期望得到的线损率,或者是否达到减出力发电机组集合的调节下限,或者是否达到加出力发电机组集合的调节上限,如果达到,则执行步骤108;否则,令k=k+1且l=l+1,返回步骤104。减出力发电机组集合的调节下限为
Figure BDA00001151222800001210
其中,
Figure BDA00001151222800001211
为减出力发电机组集合的调节下限,n为减出力发电机组的个数。加出力发电机组集合的调节上限为
Figure BDA0000115122280000131
其中,
Figure BDA0000115122280000132
为加出力发电机组集合的调节上限,m为加出力发电机组的个数。
步骤108:结束。
利用反向等量配对原则调整节点的无功补偿量具体包括下列步骤:
步骤201:当待降损的线路的输送功率对第二类节点注入无功功率的灵敏度大于零,即
Figure BDA0000115122280000133
时,将其等效为减出力节点,并将其加入减出力节点集合QPLUS
当待降损的线路的输送功率对第二类节点注入无功功率的灵敏度小于零,即
Figure BDA0000115122280000134
时,将其等效为加出力节点,并将其加入加出力节点集合QMINUS
当待降损的线路的输送功率对第二类节点注入无功功率的灵敏度等于零,即
Figure BDA0000115122280000135
时,将其等效为平衡功率节点,并将其加入平衡功率节点集合QZERO
步骤202:按照节点注入无功功率的灵敏度的绝对值的大小,分别对减出力节点集合GPLUS中的各个减出力节点和加出力节点集合GMINUS中的各个加出力节点进行排序;排序后减出力节点集合中的元素记为
Figure BDA0000115122280000136
为排序后第k个减出力节点注入无功功率的灵敏度,且有
Figure BDA0000115122280000137
k=1,2,...,n-1,n为减出力节点的个数;排序后加出力节点集合记为
Figure BDA0000115122280000139
为排序后第l个加出力节点注入无功功率的灵敏度,且有l=1,2,...,m-1,m为加出力节点的个数。
步骤203:令k=1且l=1。
步骤204:判断k>n是否成立,如果成立,则从平衡功率节点集合QZERO中取出一个未使用的元素,将其作为第k个减出力节点注入无功功率的灵敏度
Figure BDA0000115122280000141
并将该取出的未使用的元素标识为已使用元素。
步骤205:判断l>m是否成立,如果成立,则从平衡功率节点集合QZERO中取出一个未使用的元素,将其作为第l个加出力节点注入无功功率的灵敏度
Figure BDA0000115122280000142
并将该取出的未使用的元素标识为已使用元素。
步骤206:减少第k个减出力节点的无功注入,减少量为
Figure BDA0000115122280000143
增加第l个加出力节点的无功注入,增加量为其中,
ΔQ PLUS k = ΔQ MINUS l = min { ΔQ i S QN , PLUS k - S QN , MINUS l , ΔQ PLUS , k max , ΔQ MINUS , l max }
为第k个减出力节点注入无功功率与该节点无功补偿量的下限的差值;
Figure BDA0000115122280000147
为第l个加出力节点无功补偿量的上限与该节点注入无功功率的差值。
步骤207:判断待降损的线路是否达到降损后期望得到的线损率,或者是否达到减出力节点集合的调节下限,或者是否达到加出力节点集合的调节上限,如果达到,则执行步骤208;否则,令k=k+1且l=l+1,返回步骤204。减出力节点集合的调节下限为
Figure BDA0000115122280000148
其中,为减出力节点集合的调节下限,n为减出力节点的个数。加出力节点集合的调节上限为
Figure BDA00001151222800001410
其中,为加出力节点集合的调节上限,m为加出力节点的个数。
步骤208:结束。
实施例2
按照上述方法,实例2采用IEEE 118节点系统(图2为该系统接线图)作为校验模型,分析如下:
降损前,在330kV的线路中,线路38的线损率最高,为1.7%。因此,对该线路进行降损。降损前线路的潮流计算结果如表1所示。
Figure BDA0000115122280000151
表1降损前线路的潮流计算结果
设降损后线路的线损率降低为1%,则线路功率的变化量为ΔPi=97.9026MW,ΔQi=-5.0153MW。
通过灵敏度计算,节点26是最适合减出力的发电机节点,其灵敏度系数为-0.556494,节点8是最适合加出力的发电机节点,其灵敏度系数为0.123906。节点27是最适合减少无功补偿容量的节点,其灵敏度系数为0.119631,节点17是最适合增加无功补偿容量的节点,其灵敏度系数为0.115739。因此,发电机节点和无功补偿节点所需的调整量为:
ΔP G = ΔP i S G B - S G A = 97.9026 0.123906 - ( - 0.556494 ) = 143.89 MW
ΔQ G = ΔQ i S Q B - S Q A = - 5.0153 0.119634 - ( - 0.115739 ) = - 21.31 MW
改变发电机节点的注入功率和无功补偿节点的补偿容量后,进行潮流计算,结果如表2所示。
Figure BDA0000115122280000154
表2降损后线路的潮流计算结果
实例分析表明:本方法克服了传统方法中需要改造电网设备来实现降损的缺陷,仅利用电网中发电机的可调出力和节点的无功补偿裕量就实现了对线损过高支路的降损。另外,由于该方法采用了稀疏矢量法等手段,大大提高了计算速度,因此更适合于在线应用,进而辅助调度人员快速制定降损策略。
本发明所提供的方法是为在线降损提供辅助策略的一种有效手段。该方法是在保证全网线损不升高的前提下,以线损率最小为目标函数,通过电网运行方式的调整,实现对某一线损过高的支路进行降损。与传统的通过电网改造来实现降损的方法相比,该方法基于电网实时数据,根据交流灵敏度法和反向等量配对原则调整发电机出力或节点无功补偿量,改变电网潮流分布,从而有效地降低设备线损。基于注入功率灵敏度的电网局部降损方法具有以下特点:
(1)局部性,针对电网中某一线损过高的支路进行局部降损,而非对全网的优化;
(2)实时性,通过改变潮流分布,即改变发电机出力和节点无功补偿量来降低线损,而非对设备进行改造;
(3)快速性,采用了稀疏矢量法及多种提高计算速度的手段,能够快速计算分析并给出最终的调整方案;
(4)准确性,采用了交流灵敏度法及反向等量配对原则,能够反映系统的实际情况,所得结果准确可信。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。

Claims (9)

1.一种基于注入功率灵敏度的电网局部降损方法,其特征是所述方法包括:
步骤1:计算不同电压等级下各个线路的线损及线损率,并按照线损率排序;
步骤2:确定待降损的线路,设定降损后期望得到的线损率;
步骤3:计算待降损的线路的输送功率的变化量;
所述计算待降损的线路的输送功率的变化量具体是:
当待降损的线路满足
Figure FDA0000364106450000011
时,则仅计及待降损的线路的无功变化,待降损的线路的输送功率的变化量为
Figure FDA0000364106450000012
当待降损的线路满足
Figure FDA0000364106450000013
时,则仅计及待降损的线路的有功变化,待降损的线路的输送功率的变化量为
Figure FDA0000364106450000014
当待降损的线路满足
Figure FDA0000364106450000015
时,待降损的线路的输送功率的变化量为 Δ P i = P i - P i ′ = ( 1 - b a ) P i , Δ Q i = Q i - Q i ′ = ( 1 - b a ) Q i ;
其中,Pi为待降损的线路(i,j)所输送的有功功率,Qi为待降损的线路(i,j)所输送的无功功率,Pi′为降损后线路(i,j)所输送的有功功率,Qi′为降损后线路(i,j)所输送的无功功率,a为线损率,b为降损后期望得到的线损率,D1为第一设定阈值,D2为第二设定阈值,i和j分别为线路(i,j)两端的节点;
步骤4:计算待降损的线路的输送功率对各个第一类节点有功出力的灵敏度和各个第二类节点注入无功功率的灵敏度;
所述第一类节点为连接发电机组的节点;
所述第二类节点为连接无功补偿设备的节点;
步骤5:利用反向等量配对原则调整发电机有功出力和节点的无功补偿量。
2.根据权利要求1所述的基于注入功率灵敏度的电网局部降损方法,其特征是所述计算不同电压等级下各个线路的线损具体利用公式
Figure FDA0000364106450000021
其中,PLOSS为线路(i,j)的线损,i和j分别为线路(i,j)两端的节点,Pi为线路(i,j)所输送的有功功率,Qi为线路(i,j)所输送的无功功率,Ui为线路(i,j)的节点i端的电压,R为线路(i,j)的电阻。
3.根据权利要求1所述的基于注入功率灵敏度的电网局部降损方法,其特征是所述计算不同电压等级下线路的线损率具体利用公式
Figure FDA0000364106450000022
其中,r为线路(i,j)的线损率,Pi为线路(i,j)所输送的有功功率,PLOSS是线路(i,j)的线损,i和j分别为线路(i,j)两端的节点。
4.根据权利要求1所述的基于注入功率灵敏度的电网局部降损方法,其特征是所述第一设定阈值D1=100。
5.根据权利要求1所述的基于注入功率灵敏度的电网局部降损方法,其特征是所述第二设定阈值D2=0.01。
6.根据权利要求1所述的基于注入功率灵敏度的电网局部降损方法,其特征是所述计算待降损的线路的输送功率对各个第一类节点有功出力的灵敏度采用公式
S PG i = ∂ P i ∂ P G = ∂ P i ∂ δ i ∂ δ i ∂ P G + ∂ P i ∂ δ j ∂ δ j ∂ P G + ∂ P i ∂ U i U i · ∂ U i ∂ P G 1 U i + ∂ P i ∂ U j U j · ∂ U j ∂ P G 1 U j
其中,
Figure FDA0000364106450000024
为待降损的线路(i,j)的输送功率对各个第一类节点有功出力的灵敏度,Pi为待降损的线路(i,j)的输送的有功功率,PG为发电机有功出力,Ui为节点i的电压幅值,δi为节点i的电压相角,Uj为节点j的电压幅值,δj为节点j的电压相角。
7.根据权利要求1所述的基于注入功率灵敏度的电网局部降损方法,其特征是所述计算待降损的线路的输送功率对各个第二类节点注入无功功率的灵敏度采用公式
S QN i = ∂ Q i ∂ Q N = ∂ Q i ∂ δ i ∂ δ i ∂ Q N + ∂ Q i ∂ δ j ∂ δ j ∂ Q N + ∂ Q i ∂ U i U i · ∂ U i ∂ Q N 1 U i + ∂ Q i ∂ U j U j · ∂ U j ∂ Q N 1 U j
其中,为待降损的线路(i,j)的输送功率对各个第二类节点注入无功功率的灵敏度,Qi为待降损的线路(i,j)的输送的无功功率,QN为第二类节点的注入无功功率,Ui为节点i的电压幅值,δi为节点i的电压相角,Uj为节点j的电压幅值,δj为节点j的电压相角。
8.根据权利要求1所述的基于注入功率灵敏度的电网局部降损方法,其特征是所述利用反向等量配对原则调整发电机有功出力具体包括下列步骤:
步骤101:当待降损的线路的输送功率对第一类节点有功出力的灵敏度大于零,即
Figure FDA0000364106450000033
时,将与该第一类节点对应的发电机组等效为减出力发电机组,并将所述减出力发电机组加入减出力发电机组集合GPLUS
当待降损的线路的输送功率对第一类节点有功出力的灵敏度小于零,即
Figure FDA0000364106450000034
时,将与该第一类节点对应的发电机组等效为加出力发电机组,并将所述加出力发电机组加入加出力发电机组集合GMINUS
当待降损的线路的输送功率对发电机节点有功出力的灵敏度等于零,即
Figure FDA0000364106450000035
时,将与该第一类节点对应的发电机组等效为平衡功率发电机组,并将所述平衡功率发电机组加入平衡功率发电机组集合GZERO
步骤102:按照有功出力的灵敏度的绝对值的大小,分别对减出力发电机组集合GPLUS中的各个减出力发电机组和加出力发电机组集合GMINUS中的各个加出力发电机组进行排序;排序后减出力发电机组集合中的元素记为
Figure FDA0000364106450000041
为排序后第k个减出力发电机组有功出力的灵敏度,且有
Figure FDA0000364106450000042
k=1,2,...,n-1,n为减出力发电机组的个数;排序后加出力发电机组集合记为
Figure FDA0000364106450000043
为排序后第l个加出力发电机组有功出力的灵敏度,且有l=1,2,...,m-1,m为加出力发电机组的个数;
步骤103:令k=1且l=1;
步骤104:判断k>n是否成立,如果成立,则从平衡功率发电机组集合GZERO中取出一个未使用的元素,将其作为第k个减出力发电机组有功出力的灵敏度
Figure FDA0000364106450000045
并将该取出的未使用的元素标识为已使用元素;
步骤105:判断l>m是否成立,如果成立,则从平衡功率发电机组集合GZERO中取出一个未使用的元素,将其作为第l个加出力发电机组有功出力的灵敏度
Figure FDA0000364106450000046
并将该取出的未使用的元素标识为已使用元素;
步骤106:减少第k个减出力发电机组的有功出力,减少量为
Figure FDA0000364106450000047
增加第l个加出力发电机组的有功出力,增加量为其中,
Δ P PLUS k = Δ P MINUS l = min { Δ P i S PG , PLUS k - S PG , MINUS l , Δ P PLUS , k max , Δ P MINUS , l max }
Figure FDA00003641064500000410
为第k个减出力发电机组有功出力与该发电机组有功出力的下限的差值;为第l个加出力发电机组有功出力的上限与该发电机组有功出力的差值;
步骤107:判断待降损的线路是否达到降损后期望得到的线损率,或者是否达到减出力发电机组集合的调节下限,或者是否达到加出力发电机组集合的调节上限,如果达到,则执行步骤108;否则,令k=k+1且l=l+1,返回步骤104;
所述减出力发电机组集合的调节下限为
Figure FDA0000364106450000051
其中,
Figure FDA0000364106450000052
为减出力发电机组集合的调节下限,n为减出力发电机组的个数;
所述加出力发电机组集合的调节上限为
Figure FDA0000364106450000053
其中,
Figure FDA0000364106450000054
为加出力发电机组集合的调节上限,m为加出力发电机组的个数;
步骤108:结束。
9.根据权利要求1所述的基于注入功率灵敏度的电网局部降损方法,其特征是所述利用反向等量配对原则调整节点的无功补偿量具体包括下列步骤:
步骤201:当待降损的线路的输送功率对第二类节点注入无功功率的灵敏度大于零,即
Figure FDA0000364106450000055
时,将其等效为减出力节点,并将其加入减出力节点集合QPLUS
当待降损的线路的输送功率对第二类节点注入无功功率的灵敏度小于零,即
Figure FDA0000364106450000056
时,将其等效为加出力节点,并将其加入加出力节点集合QMINUS
当待降损的线路的输送功率对第二类节点注入无功功率的灵敏度等于零,即
Figure FDA0000364106450000057
时,将其等效为平衡功率节点,并将其加入平衡功率节点集合QZERO
步骤202:按照节点注入无功功率的灵敏度的绝对值的大小,分别对减出力节点集合GPLUS中的各个减出力节点和加出力节点集合GMINUS中的各个加出力节点进行排序;排序后减出力节点集合中的元素记为
Figure FDA0000364106450000058
为排序后第k个减出力节点注入无功功率的灵敏度,且有
Figure FDA0000364106450000059
k=1,2,...,n-1,n为减出力节点的个数;排序后加出力节点集合记为
Figure FDA00003641064500000510
为排序后第l个加出力节点注入无功功率的灵敏度,且有
Figure FDA00003641064500000511
l=1,2,...,m-1,m为加出力节点的个数;
步骤203:令k=1且l=1;
步骤204:判断k>n是否成立,如果成立,则从平衡功率节点集合QZERO中取出一个未使用的元素,将其作为第k个减出力节点注入无功功率的灵敏度
Figure FDA00003641064500000611
并将该取出的未使用的元素标识为已使用元素;
步骤205:判断l>m是否成立,如果成立,则从平衡功率节点集合QZERO中取出一个未使用的元素,将其作为第l个加出力节点注入无功功率的灵敏度
Figure FDA0000364106450000061
并将该取出的未使用的元素标识为已使用元素;
步骤206:减少第k个减出力节点的无功注入,减少量为
Figure FDA0000364106450000062
增加第l个加出力节点的无功注入,增加量为
Figure FDA0000364106450000063
其中,
Δ Q PLUS k = Δ Q MINUS l = min { Δ Q i S QG , PLUS k - S QG , MINUS l , Δ Q PLUS , k max , Δ Q MINUS , l max }
Figure FDA0000364106450000065
为第k个减出力节点注入无功功率与该节点无功补偿量的下限的差值;
Figure FDA0000364106450000066
为第l个加出力节点无功补偿量的上限与该节点注入无功功率的差值;
步骤207:判断待降损的线路是否达到降损后期望得到的线损率,或者是否达到减出力节点集合的调节下限,或者是否达到加出力节点集合的调节上限,如果达到,则执行步骤208;否则,令k=k+1且l=l+1,返回步骤204;
所述减出力节点集合的调节下限为
Figure FDA0000364106450000067
其中,
Figure FDA0000364106450000068
为减出力节点集合的调节下限,n为减出力节点的个数;
所述加出力节点集合的调节上限为
Figure FDA0000364106450000069
其中,
Figure FDA00003641064500000610
为加出力节点集合的调节上限,m为加出力节点的个数;
步骤208:结束。
CN 201110393178 2011-12-01 2011-12-01 基于注入功率灵敏度的电网局部降损方法 Active CN102420426B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110393178 CN102420426B (zh) 2011-12-01 2011-12-01 基于注入功率灵敏度的电网局部降损方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110393178 CN102420426B (zh) 2011-12-01 2011-12-01 基于注入功率灵敏度的电网局部降损方法

Publications (2)

Publication Number Publication Date
CN102420426A CN102420426A (zh) 2012-04-18
CN102420426B true CN102420426B (zh) 2013-10-30

Family

ID=45944709

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110393178 Active CN102420426B (zh) 2011-12-01 2011-12-01 基于注入功率灵敏度的电网局部降损方法

Country Status (1)

Country Link
CN (1) CN102420426B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106410784B (zh) * 2016-06-02 2020-03-17 国网江西省电力公司赣东北供电分公司 变电站有功负荷对地区电网有功网损的灵敏度计算方法
CN112260399B (zh) * 2020-09-18 2024-01-26 国网河北省电力有限公司行唐县供电分公司 一种利用配电监测系统的中压线损管控方法
CN113162041A (zh) * 2021-04-29 2021-07-23 南方电网电力科技股份有限公司 基于无功相关性优化线路损耗的监测方法、装置及设备

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101039031A (zh) * 2007-04-17 2007-09-19 清华大学 一种大区电网与省级电网的协调电压控制方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101039031A (zh) * 2007-04-17 2007-09-19 清华大学 一种大区电网与省级电网的协调电压控制方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
《基于灵敏度的电力系统无功优化算法研究》;葛莲等;《变频器世界》;20110731;第89页,第109-111页 *
《电力系统有功安全校正策略的反向等量配对调整法》;邓佑满;《电力系统自动化》;19990930;第23卷(第18期);5-8 *
葛莲等.《基于灵敏度的电力系统无功优化算法研究》.《变频器世界》.2011,第89页,第109-111页.
虞忠年.线路功率损耗计算.《电力网电能损耗》.中国电力出版社,2000,第5页,第36-37页. *
邓佑满.《电力系统有功安全校正策略的反向等量配对调整法》.《电力系统自动化》.1999,第23卷(第18期),5-8.

Also Published As

Publication number Publication date
CN102420426A (zh) 2012-04-18

Similar Documents

Publication Publication Date Title
CN106655227B (zh) 一种基于智能软开关的有源配电网馈线负载平衡方法
CN107316125A (zh) 一种基于经济运行域的主动配电网经济运行评价方法
Sreejith et al. Analysis of FACTS devices on security constrained unit commitment problem
Skarvelis-Kazakos et al. Implementing agent-based emissions trading for controlling Virtual Power Plant emissions
CN109523303A (zh) 一种基于节点边际电价的低压有源配电网阻塞管理方法
CN101609989B (zh) 一种城市电网供电充裕度计算系统
CN102522756B (zh) 一种避免电压越限风险的电网感性无功补偿方法
CN103050985B (zh) 一种风储系统广域优化配置的方法
CN106549402B (zh) Tcsc抑制交直流并联运行系统次同步振荡的方法
CN109066694A (zh) 含线间潮流控制器的电力系统的多目标潮流优化方法
CN105046584A (zh) 一种基于k-means算法的理想线损率的计算方法
CN104794531A (zh) 基于网损等值负荷模型的改进直流最优潮流方法
CN104701858A (zh) 计及分区动态无功储备的无功电压控制方法
CN103824124B (zh) 一种针对电网公司的节能潜力评估方法
CN102420426B (zh) 基于注入功率灵敏度的电网局部降损方法
CN106786550A (zh) 一种微电网成本优化的分布式控制方法及装置
CN103887823A (zh) 基于模糊层次分析的微电网并网位置选取方法
CN105279297B (zh) 一种基于工程应用的pss/e电力系统等值方法
CN105896547A (zh) 一种风电接入下的大电网分级电压控制方法
CN103178551B (zh) 海上油田群电网的功率优化控制方法
CN112600244A (zh) 一种基于神经网络的光伏电站电压控制方法
Zhang et al. Day-ahead stochastic optimal dispatch of LCC-HVDC interconnected power system considering flexibility improvement measures of sending system
CN109802448A (zh) 一种可再生能源最大消纳容量分析计算方法
CN115528670A (zh) 基于多主体博弈的分布式电源及配套电网投资决策建模方法
Esmaeili et al. A new multiobjective optimal allocation of multitype FACTS devices for total transfer capability enhancement and improving line congestion using the harmony search algorithm

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: NORTH CHINA ELECTRIC POWER UNIVERSITY STATE GRID C

Free format text: FORMER OWNER: NORTH CHINA ELECTRIC POWER UNIVERSITY

Effective date: 20121126

C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20121126

Address after: Qilihe District of Gansu city in Lanzhou Province, 730050 West East Road No. 648

Applicant after: Gansu Electric Power Research Institute

Applicant after: North China Electric Power University

Applicant after: State Grid Corporation of China

Address before: Qilihe District of Gansu city in Lanzhou Province, 730050 West East Road No. 648

Applicant before: Gansu Electric Power Research Institute

Applicant before: North China Electric Power University

C14 Grant of patent or utility model
GR01 Patent grant