CN102104252A - 一种适用于电磁暂态分析的电力系统动态等值方法 - Google Patents

一种适用于电磁暂态分析的电力系统动态等值方法 Download PDF

Info

Publication number
CN102104252A
CN102104252A CN2011100469576A CN201110046957A CN102104252A CN 102104252 A CN102104252 A CN 102104252A CN 2011100469576 A CN2011100469576 A CN 2011100469576A CN 201110046957 A CN201110046957 A CN 201110046957A CN 102104252 A CN102104252 A CN 102104252A
Authority
CN
China
Prior art keywords
bus
equivalent
class
generator
equivalence
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.)
Granted
Application number
CN2011100469576A
Other languages
English (en)
Other versions
CN102104252B (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.)
Zhejiang University ZJU
Zhejiang Electric Power Test and Research Insititute
Original Assignee
Zhejiang University ZJU
Zhejiang Electric Power Test and Research Insititute
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 Zhejiang University ZJU, Zhejiang Electric Power Test and Research Insititute filed Critical Zhejiang University ZJU
Priority to CN2011100469576A priority Critical patent/CN102104252B/zh
Publication of CN102104252A publication Critical patent/CN102104252A/zh
Application granted granted Critical
Publication of CN102104252B publication Critical patent/CN102104252B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种适用于电磁暂态分析的电力系统动态等值方法,它保留电力系统中500kV及以上电压等级的主干网架,并将500kV母线分为两类,分别对第一类500kV母线和第二类500kV母线分解进行等值,且两者采用不同的等值模型。动态仿真分析结果表明利用本发明等值方法进行等值电力系统等值前后的动态特性基本一致,说明了该等值方法的可行性和有效性。

Description

一种适用于电磁暂态分析的电力系统动态等值方法
技术领域
本发明涉及电力系统动态等值方法,尤其涉及一种适用于电磁暂态分析的电力系统动态等值方法。
背景技术
随着我国交直流大电网的发展,使得交直流输电系统越来越庞大。现有的交直流电力系统机电暂态仿真程序只能模拟基波对称情况下的系统,对偏离理想假设条件的运行状态则不能给出精确和严格的分析;交直流混合系统不对称故障情况下,再按照机电暂态仿真程序的正序分量计算会得出与实际情况不符合的结果,所以目前应用于电力系统中的机电暂态仿真程序无法满足系统不对称故障期间的仿真要求。为了准确模拟含有多回大容量直流输电线路的交直流大电网在不对称故障期间的暂态特性,必须采用计算量非常大的电磁暂态仿真分析方法。对于大规模系统,直接用电磁暂态仿真计算,计算量太大,一般计算前都要对被仿真网络进行简化等值。
电力系统动态等值方法主要有同调等值法、模态等值法、辨识等值法三类,但这三类等值方法主要适合于机电暂态分析,在等值过程中会在连接内外部系统节点间引入等值的阻抗支路,有些支路会出现负电阻。但在电磁暂态仿真中无法处理负电阻支路,因此需要开发一套行之有效的适用于电磁暂态仿真的动态等值方法。
发明内容
本发明提供了一种适用于电磁暂态分析的电力系统动态等值方法,避免传统等值方法在等值过程中出现负阻抗支路的问题。
一种适用于电磁暂态分析的电力系统动态等值方法,保留电力系统中500kV及以上电压等级的主干网架,将500kV母线分为两类,其中第一类:500kV母线之间只通过500kV或以上电压等级的主干线路相连;第二类:500kV母线之间除了通过500kV或以上电压等级的主干线路相连之外,还通过低于500kV电压等级的线路相连;
所述的等值方法包括对第一类500kV母线进行等值和对第二类500kV母线进行等值,其中对第一类500kV母线进行等值为对与每条第一类500kV母线连接的低于500kV电压等级的网络逐一进行等值,包括构建等值模型和确定模型参数:
构建等值模型:将与每条第一类500kV母线连接的低于500kV电压等级的网络等值为一个等值发电机和一个等值负荷,其中等值发电机与该条第一类500kV母线通过一个等值升压变压器相连,等值负荷与该条第一类500kV母线通过一个等值降压变压器相连;
确定模型参数:确定等值发电机的容量、该条第一类500kV母线电压、从低于500kV电压等级的网络注入到该条第一类500kV母线的短路电流以及等值发电机的次暂态电抗,根据等值前后注入第一类500kV母线短路电流保持不变的原则,确定等值发电机的参数;根据等值前后电力系统潮流不变的原则,确定与该条第一类500kV母线相连的等值负荷的参数;
所述的对第二类500kV母线进行等值为将所有第二类500kV母线根据是否通过同一低于500kV电压等级的网络相连划分为不同集合,将与同一低于500kV电压等级的网络相连的第二类500kV母线划分在同一集合内,对与每个集合相关的低于500kV电压等级的网络逐一进行等值,包括构建等值模型和确定模型参数:
构建等值模型:将与每个集合内所有第二类500kV母线连接的低于500kV电压等级的网络等值为一个等值发电机和数量与该集合内第二类500kV母线相同的等值负荷,其中等值发电机与该集合内每条第二类500kV母线之间分别通过各自的等值升压变压器相连,每个等值负荷与和其对应的第二类500kV母线通过各自的等值降压变压器相连;
确定模型参数:确定等值发电机的容量、有功功率、无功功率;该集合内所有第二类500kV母线的母线电压;从低于500kV电压等级的网络注入到所有该集合内每条第二类500kV母线的短路电流以及等值发电机的次暂态电抗,根据注入所有第二类500kV母线的短路电流保持不变的原则,确定的等值发电机的参数;根据等值前后电力系统潮流不变的原则,确定每个等值负荷的参数。
所述的对第一类500kV母线进行等值中确定模型参数的方法为:
(1)获取与单条第一类500kV母线通过低于500kV电压等级的网络相连的所有发电机的容量和惯性时间常数,计算得到相应的等值发电机的容量和惯性时间常数;
(2)获取该条第一类500kV母线的电压幅值和相角,得到该条第一类500kV母线电压;
(3)根据等值前电力系统的短路电流计算得到从低于500kV电压等级的网络注入该条第一类500kV母线的短路电流;
(4)以等值发电机的容量、等值发电机次暂态电抗、该条第一类500kV母线电压、从低于500kV电压等级的网络注入到该条第一类500kV母线的短路电流为已知量,在等值发电机功率因数一定的条件下,通过电力系统正常时的功率平衡方程以及电力系统短路情况下的基尔霍夫定律方程,求解等值升压变压器的阻抗和等值发电机的机端电压相量;通过欧姆定律方程和功率平衡方程,求解等值负荷的值和等值降压变压器的漏抗。
所述的对第二类500kV母线进行等值中确定模型参数的方法为:
(1)获取与同一集合内所有第二类500kV母线通过低于500kV电压等级的网络相连的发电机的容量、惯性时间常数、有功功率及无功功率,计算得到等值发电机的容量、惯性时间常数、有功功率及无功功率;
(2)获取该集合内第二类500kV母线的电压幅值和相角,得到该集合内所有第二类500kV母线的母线电压;
(3)根据等值前电力系统的短路电流计算分别得到从低于500kV电压等级的网络注入该集合内每条第二类500kV母线的短路电流;
(4)以等值发电机的容量、有功功率、无功功率、次暂态电抗;所有第二类500kV母线的母线电压;从低于500kV电压等级的网络注入到所有第二类500kV母线的短路电流;在等值发电机发出有功功率条件下,通过电力系统正常时的功率平衡方程以及电力系统短路情况下的基尔霍夫定律方程,求解等值升压变压器的阻抗和等值发电机的机端电压相量;通过欧姆定律方程和功率平衡方程,求解等值负荷的值和等值降压变压器的漏抗。
使用本发明等值方法对南方电网2013年规划数据中的全部500kV母线进行等值。等值前整个电网中含有2981条母线、555台发电机、7条直流线路,等值后只有1017条母线、240台发电机、7条直流线路,等值前后500kV主干网架上的潮流完全一致。对原型系统和等值系统进行500kV母线三相短路故障分析,图8-11为等值前后云广直流逆变侧附近主干网架上三相短路0.1秒后切除故障各个地区机组的功角曲线,动态仿真分析结果表明等值前后系统的动态特性基本一致,从而验证了该等值方法的可行性和有效性。
附图说明
图1为第一类500kV母线等值后的结构;
图2为正常潮流情况下第一类500kV母线等值后的结构的简化等效电路图;
图3为三相短路情况下第一类500kV母线等值后结构的简化等效电路图;
图4为第二类500kV母线等值前的结构;
图5为第二类500kV母线等值后的结构;
图6为第二类500kV母线等值后短路电流计算时的等效电路图;
图7为第二类500kV母线等值后忽略
Figure BDA0000048144850000041
情况下的等效电路图;
图8为云广直流逆变侧附近主干网架三相短路0.1秒后切除故障广东某机组等值前后的功角曲线;
图9为云广直流逆变侧附近主干网架三相短路0.1秒后切除故障广西某机组等值前后的功角曲线;
图10为云广直流逆变侧附近主干网架三相短路0.1秒后切除故障贵州某机组等值前后的功角曲线;
图11为云广直流逆变侧附近主干网架三相短路0.1秒后切除故障云南某机组等值前后的功角曲线。
具体实施方式
一种用于电磁暂态分析的电力系统动态等值方法,该方法将系统按照电压等级进行分层,认为交直流互联大系统的动态特性主要取决大容量直流线路接入层的系统强度,即所谓的发电机对主干网架的支撑强度。
等值前后需要满足以下几个条件:(1)等值前后系统的发电机总容量保持不变;(2)等值前后主干网架上各条母线的短路电流保持不变;(3)等值前后系统的潮流保持一致;(4)等值前后系统的动态特性保持一致,从而达到动态等值的目的。
进行动态等值时应该对以下系统进行保留:(1)系统500kV及以上电压等级的主干网架;(2)各回直流输电线路的详细模型;(3)与直流并联运行的交流通道;(4)与500kV主干网架直接相连的发电容量较大的发电机组。
根据上述等值原理,该等值方法如下:对原型系统中500kV母线进行分类。系统主干网架上500kV母线可以按要求分为两类。第一类:500kV母线之间只通过500kV或以上电压等级的主干线路相连;第二类:500kV母线之间除了通过500kV或以上电压等级的主干线路相连之外,还通过低于500kV电压等级的线路相连。
对于第一类500kV母线,等值后系统主干网架上500kV母线标准结构如图1所示。它将与每条第一类500kV母线连接的低于500kV电压等级的网络等值为一个等值发电机和一个等值负荷,其中等值发电机与该条第一类500kV母线通过一个等值升压变压器相连,等值负荷与该条第一类500kV母线通过一个等值降压变压器相连;
以下重点研究500kV母线到等值机之间升压变压器以及等值机的参数计算和调整问题,原来直接连接在500kV母线上的发电机组给予保留。为了方便建立计算公式,现在设发电机母线电压Ut∠θt是未知量,500kV母线电压Us∠θs是已知量可由正常潮流计算得到,升压变压器阻抗设为Zt=jXt(忽略电阻Rt)是未知量,等值发电机次暂态电抗设为
Figure BDA0000048144850000051
,等值发电机的容量SG可以根据潮流数据计算得出为已知量,从低于500kV电压等级的网络注入到500kV母线的短路电流I″可以根据短路电流计算得出为已知量(I″是指原型系统中该500kV母线下属的所有低于500kV的线路注入该母线的短路电流,与该500kV母线的总短路电流是不一样的)。
结合工程实际,把等值机的次暂态电抗
Figure BDA0000048144850000052
用典型值,
Figure BDA0000048144850000053
(以等值机容量为基准的标幺值),现在已知等值发电机的容量SG,等值发电机次暂态电抗设为
Figure BDA0000048144850000061
,500kV母线电压为Us∠θs,从低于500kV电压等级的网络注入到500kV母线的短路电流
Figure BDA0000048144850000062
在等值发电机功率因数一定的条件下,可以根据功率平衡方程以及短路情况下的基尔霍夫定律方程联立求解升压变压器阻抗和等值发电机的机端电压相量,图2和图3给出了计算时的等效电路图。以下推导过程中的各物理量的参考方向参见图1~图3,并且各物理量均归算至标么值。
三相短路时:
E → ′ ′ = j ( X d ′ ′ + X t ) I → ′ ′ - - - ( 1 )
上式中,
Figure BDA0000048144850000064
是等值发电机的次暂态电势,是等值发电机次暂态电抗,Xt是升压变压器的短路电抗,是从低于500kV电压等级的网络注入到被等值500kV母线的短路电流。
稳态情况下,根据图2:
E → ′ ′ = U → s + j ( X d ′ ′ + X t ) I → t - - - ( 2 )
U → t = U → s + j X t I → t - - - ( 3 )
上面两式中,
Figure BDA0000048144850000069
是被等值的500kV母线的电压相量,
Figure BDA00000481448500000610
是流过升压变压器的稳态电流,
Figure BDA00000481448500000611
是等值发电机的机端电压相量。
忽略变压器对地的等值导纳:
P + jQ = U → t I ^ t - - - ( 4 )
上式中,P,Q分别是等值发电机发出的有功和无功.
将式(1)代入式(2),可得:
I → t = I → ′ ′ - U → s j ( X d ′ ′ + X t ) - - - ( 5 )
由式(4)知:
U → t = P + jQ I ^ t - - - ( 6 )
由式(3)、(6)消去
Figure BDA00000481448500000615
可得:
P + jQ = U → s I ^ t + j | I → t | 2 Xt - - - ( 7 )
将式(5)取共轭:
I ^ t = I ^ ′ ′ - U ^ s - j ( X d ′ ′ + X t ) - - - ( 8 )
将式(8)代入式(7)得:
P + jQ = U → s ( I ^ ′ ′ - U ^ s - j ( X d ′ ′ + X t ) ) + j | I → t | 2 Xt - - - ( 9 )
取式(9)的实部:
P = real ( U → s I ^ ′ ′ ) - - - ( 10 )
取式(9)的虚部并结合式(5):
Q = imag ( U → s I ^ ′ ′ ) - | U ^ s | 2 ( X d ′ ′ + X t ) ) + | I → ′ ′ - U → s j ( X d ′ ′ + X t ) | 2 Xt - - - ( 11 )
等值后发电机出力P可以通过
Figure BDA0000048144850000074
Figure BDA0000048144850000075
根据式(10)求解得到。在等值发电机功率因数一定的条件下,可以求解得到Q,则式(11)是一个只含有一个未知数Xt的一元三次方程,可以求解得到发电机升压变压器电抗参数Xt。根据式(6)和式(8)可以求解得到机端电压相量
Figure BDA0000048144850000076
求解出发电机升压变压器电抗参数Xt和机端电压相量
Figure BDA0000048144850000077
之后,只需调整负荷大小就可以保证等值后主干网架上线路与原型系统主干网架上线路潮流一致。调整负荷时需要从潮流数据中获取主干网架上流出等值节点的潮流,记为Fout,计算过程用到的变量及参考方向参见图1。
Figure BDA0000048144850000078
求出以后,可以求得稳态时等值机注入500kV母线的电流:
I → t = ( U → t - U → s ) Z t - - - ( 12 )
稳态时等值机注入500kV母线的功率:
P 2 = jQ 2 = U → s I ^ t - - - ( 13 )
等值负荷流出500kV母线的功率:
P1+jQ1=P2+jQ2-Fout    (14)
在与等值负荷直接相连的降压变压器的容量就等于流过变压器的视在功率的情况下,并设此变压器的Uk%为某一定值,则可以求得变压器的漏抗XL,则流过变压器的稳态电流和负荷节点的电压:
I → L = ( P 1 - jQ 1 ) U ^ s - - - ( 15 )
U → L = U → s - j I → l X L - - - ( 16 )
等值负荷大小为:
P 3 + j Q 3 = U → L I ^ L - - - ( 17 )
到此,第一类500kV母线的等值过程所需要的所有参数求解完毕。
对于上述第二类500kV母线的等值,不妨设某个第二类500kV母线集合S内有m条第二类500kV母线通过同一低于500kV电压等级的网络相连,结构图如图4所示,这种结构采用如图5的模型来进行等值。它将与每个集合内所有第二类500kV母线连接的低于500kV电压等级的网络等值为一个等值发电机和数量与该集合内第二类500kV母线相同的等值负荷,其中等值发电机与该集合内每条第二类500kV母线之间分别通过各自的等值升压变压器相连,每个等值负荷与和其对应的第二类500kV母线通过各自的等值降压变压器相连;
与第一类500kV母线等值方法类似,重点研究的是500kV母线到等值机之间升压变压器以及等值机的参数计算和调整,原来直接连接在500kV母线的发电机组给予保留。现在设发电机机端母线电压Ute∠θte是未知量,可以从潮流数据中获取集合S内所有第二类500kV母线的电压幅值和相角,得到集合S内500kV各母线电压{U1∠θ1,U2∠θ2,...},各升压变压器阻抗设为{Z1=R1+jX1,Z2=R2+jX2,...}是未知量,结合工程实际,等值发电机次暂态电抗设为
Figure BDA0000048144850000081
(以自身容量为基准的标幺值),等值发电机的容量SG、有功出力P、发出无功功率Q可以根据潮流数据计算得出作为已知量,集合S内500kV各母线短路时从低于500kV电压等级的网络注入到集合内每条第二类500kV母线的短路电流
Figure BDA0000048144850000082
可以根据短路电流计算得出。求解思路是根据功率平衡方程及基尔霍夫定律方程,求解得到各升压变压器阻抗和等值发电机的机端电压相量。以下给出具体求解过程:
不失一般性,设集合S内共有m条第二类500kV母线,并且近似地认为等值机的次暂态电势
Figure BDA0000048144850000083
且假定升压变压器阻抗Zi=Ri+jXi中Xi是Ri的10倍。因为等值发电机相对于基准值容量很大,把
Figure BDA0000048144850000084
(以自身容量为基准的标幺值)转化成系统容量时等值机的次暂态电抗很小,相对于升压变压器阻抗要小得多。于是图6中当第i个母线短路时,
Figure BDA0000048144850000085
要比
Figure BDA0000048144850000086
小得多,等值机经过第i个升压变压器注入到短路点的短路电流大小主要取决于第i个升压变压器的阻抗,所以在近似计算第i个升压变压器的阻抗时,可以忽略
Figure BDA0000048144850000087
的影响,见图7。等值发电机经过升压变压器连接到被等值500kV母线,可以作以下近似:
X di ≈ | E → ′ ′ | | I → i ′ ′ | - - - ( 18 )
R di = X di 10 - - - ( 19 )
上面两式中,
Figure BDA0000048144850000093
是原系统中从低于500kV电压等级的网络注入到第i个第二类500kV母线的短路电流,Rdi,Xdi分别是第i个升压变压器的短路电阻和短路电抗,并且Zi=Ri+jXi
另一方面,等值机的次暂态电抗相对于升压变压器阻抗要小得多,所以第i个升压变压器的阻抗参数可近似地看成
Z i = R i + jX i ≈ Z di - j X d ′ ′ - - - ( 20 )
设正常运行时等值机机端输出的电流为
Figure BDA0000048144850000096
则:
I → te = Σ i = 1 m ( U → te - U → i ) Z i - - - ( 21 )
等值机输出的复数功率:
P + jQ = U → te I ^ te - - - ( 22 )
Figure BDA0000048144850000099
将(21)式代入(22)式,并分开实部和虚部可得:
P=a(U2 tex+U2 tey)+bUtex+cUtey    (23)
Q=d(U2 tex+U2 tey)-cUtex+bUtey    (24)
其中:
a = Σ i = 1 m R i R i 2 + X i 2 - - - ( 25 )
b = - Σ i = 1 m U ix R i + U iy X i R i 2 + X i 2 - - - ( 26 )
c = Σ i = 1 m U ix X i - U iy R i R i 2 + X i 2 - - - ( 27 )
d = Σ i = 1 m X i R i 2 + X i 2 - - - ( 28 )
式(23)、(24)是一个含有两个未知数Utex和Utey的二元二次方程组,由于两个方程相似,可以消去二次项,再代入到式(15),得到一个一元二次方程,进而求解出Utex和Utey
求解出发电机升压变压器电抗参数Xi和机端电压相量之后,只需调整负荷大小就可以保证等值后主干网架上线路与原型系统主干网架上线路潮流一致。求解负荷大小的方法与第一类500kV母线等值时求解负荷的方法一样。
所有参数求解完之后,保留原型系统中的主干网架,删除被等值掉的网架,并在此基础上搭建好等值后的网架,添加等值机的典型参数,便完成了整个等值过程。
为了更好地说明等值过程,以下举两个例子详细叙述各参数的计算、选择过程,以便公众更好地掌握本发明的具体实施手段。
实施例1
选取南方电网2013年规划数据中广东东莞500kV母线的等值作为实例。
(1)确定母线类型
搜索是否有等于或者高于500kV电压等级的母线通过低于500kV电压等级的线路与东莞500kV母线相连,发现没有这样的母线,所以东莞500kV母线属于第一类母线。
(2)确定等值发电机的容量和惯性时间常数
遍历与东莞500kV母线通过低于500kV电压等级的线路相连的发电机,获取这些发电机的容量Si和惯性时间常数Hi(i=1,2,...N,N为发电机的台数)。等值的发电机的容量按照容量相加的原则,惯性时间常数HG的计算思路是将各个发电机本身的Hi值换算成系统基准容量SB下的值后再相加,最后再转换成以等值机容量SG为基准的值,系统参考容量基准值取100MVA,等值机惯性时间常数和容量计算如下:
S G = Σ i = 1 N S i - - - ( 29 )
H G = ( H 1 * S 1 S B + H 2 * S 2 S B + . . . + H N * S N S B ) * S B S G (30)
= H 1 * S 1 S G + H 2 * S 2 S G + . . . + H N * S N S G
求解得到SG=22.35(pu),HG=2.8814。
(3)获取东莞500kV母线的电压幅值和相角
从潮流数据中获取东莞500kV母线的电压幅值和相角得到电压相量Us∠θs=1.056609∠5.445°(pu)。
(4)短路电流计算
通过短路电流计算得到东莞500kV母线从低于500kV电压等级的网络注入的短路电流。东莞500kV母线下共有4条低于500kV的线路与之相连,4条线路注入到东莞500kV母线的短路电流都等于1.76-10.51j(pu),相加得到
Figure BDA0000048144850000111
(5)发电机升压变压器电抗参数Xt和机端电压相量的计算
根据式(10)计算得到
Figure BDA0000048144850000113
如果假设等值机按照功率因素0.9发出有功,则Q=1.549(pu),结合式(11),得出Xt=0.01667(pu)。再根据式(6)和式(8)可以得到
(6)确定等值负荷及负荷所连变压器参数
从潮流数据中获取主干网架上流出等值节点的潮流,主干网架上与东莞500kV母线相连的线路有两条双回线。将4条支路上流出东莞500kV母线的潮流相加,得到Fout=18.474-2.126j(pu)。根据式(13)和式(14)求得P1+jQ1=21.672-0.757j(pu)。若变压器短路电压百分数Uk%=12,变压器容量按
Figure BDA0000048144850000115
计算,则可得到变压器漏抗XL=0.005533(pu),再根据式(15)、(16)和(17)解得P3+jQ3=21.672-3.088j(pu)。
(7)用求解出来的参数搭建出等值系统,添加等值机的典型参数。
实施例2
选取南方电网2013年规划数据中云南厂口500kV母线的等值作为实例。
(1)确定母线类型
搜索是否有等于或者高于500kV电压等级的母线通过低于500kV电压等级的线路与厂口500kV母线相连,发现只有曲靖500kV母线与厂口500kV母线通过低于500kV电压等级的线路间接相连,所以厂口和曲靖500kV母线属于第二类母线。
(2)确定等值发电机的容量、惯性时间常数、有功出力和发出的无功功率
遍历与厂口及曲靖500kV母线通过低于500kV电压等级的线路相连的发电机,获取这些发电机的容量Si和惯性时间常数Hi(i=1,2,...N,N为发电机的台数)、有功出力和发出的无功功率。根据式(29)和(30)求解得到SG=20.108(pu),HG=3.2379。等值机的有功出力和发出无功功率根据累加的原则,求得P=15.04(pu),Q=16.84(pu)。
(3)获取厂口和曲靖500kV母线的电压幅值和相角
从潮流数据中获取厂口和曲靖500kV母线的电压幅值和相角得到{U1∠θ1=1.0068∠47.1°(pu),U2∠θ2=1.0156∠44.9°(pu)}。
(4)短路电流计算
通过短路电流计算得到厂口和曲靖500kV母线从低于500kV电压等级的网络注入的短路电流。厂口500kV母线下共有2条低于500kV的线路与之相连,2条线路注入到厂口500kV母线的短路电流都等于9.2481-8.7161j(pu),相加得到
Figure BDA0000048144850000121
曲靖500kV母线下共有1条低于500kV的线路与之相连,故
Figure BDA0000048144850000122
(5)发电机升压变压器电抗参数Xt和机端电压相量
Figure BDA0000048144850000123
的计算
根据式(18)计算得到
Figure BDA0000048144850000124
根据式(19)和式(20)可以得到升压变压器阻抗Z1=0.00469+0.03694j,Z2=0.00583+0.04835j。联立二元一次方程组式(23)和式(24),可以解得等值机机端电压相量
Figure BDA0000048144850000126
(6)确定等值负荷及负荷所连变压器参数
从潮流数据中获取主干网架上流出等值节点厂口500kV母线的潮流,主干网架上与厂口500kV母线相连的线路有四条双回线。将8条支路上流出厂口500kV母线的潮流相加,得到Fout1=7.956-2.558j(pu)。根据式(13)和式(14)求得P1+jQ1=15.706+2.965j。若变压器短路电压百分数Uk%=12,变压器容量按
Figure BDA0000048144850000127
计算,则可得到变压器漏抗XL=0.007508(pu),再根据式(15)、(16)和(17)解得P3+jQ3=15.706+1.073j(pu)。
同理从潮流数据中获取主干网架上流出等值节点曲靖500kV母线的潮流Fout2=6.507-0.768j(pu)。与上面类似地可以求解得到P′1+jQ′1=13.261+0.600j(pu),X′L=0.009040(pu),P′3+jQ′3=13.261-0.944j(pu)。
(7)用求解出来的参数搭建出等值系统,添加等值机的典型参数。
依照上述对方法整个电力系统进行等值,图8~10为云广直流逆变侧附近主干网架三相短路0.1秒后切除故障某机组等值前后的功角曲线,表明系统的动态特性基本一致,从而验证了该等值方法的可行性和有效性。

Claims (3)

1.一种适用于电磁暂态分析的电力系统动态等值方法,保留电力系统中500kV及以上电压等级的主干网架,将500kV母线分为两类,其中第一类:500kV母线之间只通过500kV或以上电压等级的主干线路相连;第二类:500kV母线之间除了通过500kV或以上电压等级的主干线路相连之外,还通过低于500kV电压等级的线路相连;
所述的等值方法包括对第一类500kV母线进行等值和对第二类500kV母线进行等值,其中对第一类500kV母线进行等值为对与每条第一类500kV母线连接的低于500kV电压等级的网络逐一进行等值,包括构建等值模型和确定模型参数:
构建等值模型:将与每条第一类500kV母线连接的低于500kV电压等级的网络等值为一个等值发电机和一个等值负荷,其中等值发电机与该条第一类500kV母线通过一个等值升压变压器相连,等值负荷与该条第一类500kV母线通过一个等值降压变压器相连;
确定模型参数:确定等值发电机的容量、该条第一类500kV母线电压、从低于500kV电压等级的网络注入到该条第一类500kV母线的短路电流以及等值发电机的次暂态电抗,根据等值前后注入第一类500kV母线短路电流保持不变的原则,确定等值发电机的参数;根据等值前后电力系统潮流不变的原则,确定与该条第一类500kV母线相连的等值负荷的参数;
所述的对第二类500kV母线进行等值为将所有第二类500kV母线根据是否通过同一低于500kV电压等级的网络相连划分为不同集合,将与同一低于500kV电压等级的网络相连的第二类500kV母线划分在同一集合内,对与每个集合相关的低于500kV电压等级的网络逐一进行等值,包括构建等值模型和确定模型参数:
构建等值模型:将与每个集合内所有第二类500kV母线连接的低于500kV电压等级的网络等值为一个等值发电机和数量与该集合内第二类500kV母线相同的等值负荷,其中等值发电机与该集合内每条第二类500kV母线之间分别通过各自的等值升压变压器相连,每个等值负荷与和其对应的第二类500kV母线通过各自的等值降压变压器相连;
确定模型参数:确定等值发电机的容量、有功功率、无功功率;该集合内所有第二类500kV母线的母线电压;从低于500kV电压等级的网络注入到所有该集合内每条第二类500kV母线的短路电流以及等值发电机的次暂态电抗,根据注入所有第二类500kV母线的短路电流保持不变的原则,确定的等值发电机的参数;根据等值前后电力系统潮流不变的原则,确定每个等值负荷的参数。
2.根据权利要求1所述的电力系统动态等值方法,其特征在于,所述的对第一类500kV母线进行等值中确定模型参数的方法为:
(1)获取与单条第一类500kV母线通过低于500kV电压等级的网络相连的所有发电机的容量和惯性时间常数,计算得到相应的等值发电机的容量和惯性时间常数;
(2)获取该条第一类500kV母线的电压幅值和相角,得到该条第一类500kV母线电压;
(3)根据等值前电力系统的短路电流计算得到从低于500kV电压等级的网络注入该条第一类500kV母线的短路电流;
(4)以等值发电机的容量、等值发电机次暂态电抗、该条第一类500kV母线电压、从低于500kV电压等级的网络注入到该条第一类500kV母线的短路电流为已知量,在等值发电机功率因数一定的条件下,通过电力系统正常时的功率平衡方程以及电力系统短路情况下的基尔霍夫定律方程,求解等值升压变压器的阻抗和等值发电机的机端电压相量;通过欧姆定律方程和功率平衡方程,求解等值负荷的值和等值降压变压器的漏抗。
3.根据权利要求1所述的电力系统动态等值方法,其特征在于,所述的对第二类500kV母线进行等值中确定模型参数的方法为:
(1)获取与同一集合内所有第二类500kV母线通过低于500kV电压等级的网络相连的发电机的容量、惯性时间常数、有功功率及无功功率,计算得到等值发电机的容量、惯性时间常数、有功功率及无功功率;
(2)获取该集合内第二类500kV母线的电压幅值和相角,得到该集合内所有第二类500kV母线的母线电压;
(3)根据等值前电力系统的短路电流计算分别得到从低于500kV电压等级的网络注入该集合内每条第二类500kV母线的短路电流;
(4)以等值发电机的容量、有功功率、无功功率、次暂态电抗;所有第二类500kV母线的母线电压;从低于500kV电压等级的网络注入到所有第二类500kV母线的短路电流;在等值发电机发出有功功率条件下,通过电力系统正常时的功率平衡方程以及电力系统短路情况下的基尔霍夫定律方程,求解等值升压变压器的阻抗和等值发电机的机端电压相量;通过欧姆定律方程和功率平衡方程,求解等值负荷的值和等值降压变压器的漏抗。
CN2011100469576A 2011-02-28 2011-02-28 一种适用于电磁暂态分析的电力系统动态等值方法 Expired - Fee Related CN102104252B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011100469576A CN102104252B (zh) 2011-02-28 2011-02-28 一种适用于电磁暂态分析的电力系统动态等值方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011100469576A CN102104252B (zh) 2011-02-28 2011-02-28 一种适用于电磁暂态分析的电力系统动态等值方法

Publications (2)

Publication Number Publication Date
CN102104252A true CN102104252A (zh) 2011-06-22
CN102104252B CN102104252B (zh) 2012-11-07

Family

ID=44156851

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011100469576A Expired - Fee Related CN102104252B (zh) 2011-02-28 2011-02-28 一种适用于电磁暂态分析的电力系统动态等值方法

Country Status (1)

Country Link
CN (1) CN102104252B (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102280884A (zh) * 2011-08-26 2011-12-14 华中科技大学 一种电网等值方法
CN102280883A (zh) * 2011-08-19 2011-12-14 东北电网有限公司 一种用于电力系统动态仿真验证的广域模式分析方法
CN102403716A (zh) * 2011-10-28 2012-04-04 上海电力学院 一种多馈入交直流电网动态等值的方法
CN102495262A (zh) * 2011-11-22 2012-06-13 上海电力学院 短路电流经典假设的改进计算方法
CN103051265A (zh) * 2012-12-26 2013-04-17 浙江大学 一种电力系统故障情况下的发电机强励控制方法
CN103336866A (zh) * 2013-06-27 2013-10-02 华南理工大学 一种电磁暂态仿真中含负电阻支路的处理方法
CN103559347A (zh) * 2013-10-31 2014-02-05 浙江大学 一种大规模交直流电力系统电磁暂态仿真模型的构建方法
CN105119273A (zh) * 2015-08-13 2015-12-02 国电南瑞科技股份有限公司 一种计及系统网架强度和发电机有效送出能力的电网等值简化方法
CN105160058A (zh) * 2015-07-10 2015-12-16 国网四川省电力公司电力科学研究院 一种用于电网电磁暂态计算的网络边界等值方法
CN108092266A (zh) * 2017-12-29 2018-05-29 国网浙江省电力有限公司电力科学研究院 一种考虑零序阻抗的适用于电磁暂态仿真的大电网等值方法
CN109698497A (zh) * 2017-10-23 2019-04-30 南方电网科学研究院有限责任公司 一种电力系统动态等值中母线负荷的处理方法和装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006065859A2 (en) * 2004-12-14 2006-06-22 Brian Costa Dynamic load control in a customer distribution network
CN101814740A (zh) * 2010-02-09 2010-08-25 上海交通大学 适用于次同步谐振的电力系统外网等值方法
CN101872975A (zh) * 2010-05-04 2010-10-27 国网电力科学研究院 电力系统暂态功角稳定在线分析的自适应动态等值方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006065859A2 (en) * 2004-12-14 2006-06-22 Brian Costa Dynamic load control in a customer distribution network
CN101814740A (zh) * 2010-02-09 2010-08-25 上海交通大学 适用于次同步谐振的电力系统外网等值方法
CN101872975A (zh) * 2010-05-04 2010-10-27 国网电力科学研究院 电力系统暂态功角稳定在线分析的自适应动态等值方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
徐政等: "一种适用于电磁暂态仿真的等值简约方法", 《南方电网技术》, vol. 1, no. 1, 31 October 2007 (2007-10-31), pages 37 - 40 *
杨靖萍等: "基于同调机群识别的动态等值方法的工程应用", 《电网技术》, vol. 29, no. 17, 30 September 2005 (2005-09-30), pages 68 - 71 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102280883B (zh) * 2011-08-19 2013-04-03 东北电网有限公司 一种用于电力系统动态仿真验证的广域模式分析方法
CN102280883A (zh) * 2011-08-19 2011-12-14 东北电网有限公司 一种用于电力系统动态仿真验证的广域模式分析方法
CN102280884A (zh) * 2011-08-26 2011-12-14 华中科技大学 一种电网等值方法
CN102403716A (zh) * 2011-10-28 2012-04-04 上海电力学院 一种多馈入交直流电网动态等值的方法
CN102495262B (zh) * 2011-11-22 2014-04-09 上海电力学院 短路电流经典假设的改进计算方法
CN102495262A (zh) * 2011-11-22 2012-06-13 上海电力学院 短路电流经典假设的改进计算方法
CN103051265A (zh) * 2012-12-26 2013-04-17 浙江大学 一种电力系统故障情况下的发电机强励控制方法
CN103336866A (zh) * 2013-06-27 2013-10-02 华南理工大学 一种电磁暂态仿真中含负电阻支路的处理方法
CN103336866B (zh) * 2013-06-27 2016-04-13 华南理工大学 一种电磁暂态仿真中含负电阻支路的处理方法
CN103559347A (zh) * 2013-10-31 2014-02-05 浙江大学 一种大规模交直流电力系统电磁暂态仿真模型的构建方法
CN103559347B (zh) * 2013-10-31 2016-04-06 浙江大学 一种大规模交直流电力系统电磁暂态仿真模型的构建方法
CN105160058A (zh) * 2015-07-10 2015-12-16 国网四川省电力公司电力科学研究院 一种用于电网电磁暂态计算的网络边界等值方法
CN105160058B (zh) * 2015-07-10 2018-12-25 国网四川省电力公司电力科学研究院 一种用于电网电磁暂态计算的网络边界等值方法
CN105119273A (zh) * 2015-08-13 2015-12-02 国电南瑞科技股份有限公司 一种计及系统网架强度和发电机有效送出能力的电网等值简化方法
CN109698497A (zh) * 2017-10-23 2019-04-30 南方电网科学研究院有限责任公司 一种电力系统动态等值中母线负荷的处理方法和装置
CN108092266A (zh) * 2017-12-29 2018-05-29 国网浙江省电力有限公司电力科学研究院 一种考虑零序阻抗的适用于电磁暂态仿真的大电网等值方法

Also Published As

Publication number Publication date
CN102104252B (zh) 2012-11-07

Similar Documents

Publication Publication Date Title
CN102104252B (zh) 一种适用于电磁暂态分析的电力系统动态等值方法
CN103559347B (zh) 一种大规模交直流电力系统电磁暂态仿真模型的构建方法
CN101505061B (zh) 一种基于时域仿真的戴维南等值参数跟踪的计算方法
CN102682358B (zh) 一种评估风电并网规模与电网网架适应性的规划仿真方法
CN103001214B (zh) 一种基于中性点偏移的配电网三相潮流计算方法
CN104362622B (zh) 一种带小水电机组地区电网的负荷模型和建模方法
CN103036230A (zh) 一种基于工程应用的交直流混联大电网的动态等值方法
CN101882167A (zh) 一种大规模风电集中接入电网的风电场等值建模方法
CN108092266A (zh) 一种考虑零序阻抗的适用于电磁暂态仿真的大电网等值方法
CN104716646B (zh) 一种基于注入电流的节点耦合度分析方法
CN103595040A (zh) 含直流分布式电源的配电网及其综合负荷特性的分析方法
CN101958543A (zh) 一种变电站三相无阻抗非线性多源状态估计方法
CN107121604A (zh) 一种不对称短路电流直流分量衰减时间常数获取方法
CN105184669A (zh) 一种基于节点集GN分裂算法的220kV城市环网分区方法
CN106295160A (zh) 交直流互联电网戴维南等值参数在线计算方法
CN106549402A (zh) Tcsc抑制交直流并联运行系统次同步振荡的方法
CN110333394A (zh) 一种低压配电网线路阻抗估计方法
CN106443276A (zh) 一种交流高压多回输电线路无线电干扰计算方法及系统
CN106356843B (zh) 一种反映电网电压特性的综合负荷建模方法
CN102709911B (zh) 一种混合仿真谐波特性的接口设计方法
CN105823929B (zh) 基于双端异步数据的长距离同塔四回输电线路零序参数精确测量方法
CN112966364B (zh) 一种用于特征值计算的光伏电站等值模型建模方法及装置
CN103560494B (zh) 配电网保护整定的短路电流获取方法
CN106972518A (zh) 小型地方电网与能源基地直流外送系统接入模式选择方法
CN103914738B (zh) 一种配电网评估与网架优化分析方法

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

Granted publication date: 20121107

Termination date: 20150228

EXPY Termination of patent right or utility model