Y6 ~ m) G' z b.按经验取TTC的一个固定百分比(比如4%)。这种方法简单方便,较为常用。 6 w8 L3 O6 e; X7 f, ]5 Q1 V2 o5 y* P
c.利用统计学或概率方法求解[6]。这类方法的计算过程较为繁琐,也不太成熟。 / B" o1 i e; W" t7 a* s$ `/ [! l9 K# R1 @( s& f. B- w8 K
d.降低额定值法:系统中不确定因素对所有电力设备的影响是相对一致的,可以通过降低设备的计算额定值来计及这些不确定因素。典型的额定值降低为2-5%,当然降低的幅值会随着时间跨度的增大而增大[4]。 2 N9 F" m/ s# K 5 F! w, }' s0 s" h, B 在实际应用中,并不需要实时计算TRM值,欧共体内部电力市场只要求每半年刷新一次TRM值。计算方法更多是取TTC的一个固定百分比,对于电气连接紧密,无功支持充裕的电网,TRM所占TTC的百分比可以取的小一些。6 w+ ]6 B r+ h, w8 A0 Y$ W
, v$ |3 r' J& F5 t
2.2 容量效益裕度CBM& t( k, x% J+ R, B: P
6 _" b9 D( _+ Z, m
CBM指的是:为了能够从其它互联系统中获得电力来满足发电可靠性需求而预留的输电容量裕度。 + q# Y, U+ j5 \$ ?$ q( e " V1 k b, t; N7 \; t 备用发电容量对保证供电可靠性至关重要,在发电机停运或因其它设备故障而失去部分电源时,需要启动备用电源以保证向负荷供电。对于互联的电力系统,当发生故障引起电源短缺时可以从其它子系统得到紧急援助,因而每个系统内部都可以减少发电备用容量,以提高经济性。为了保证系统在任何时刻都有能力将系统外的电能输送到负荷中心,需要预留部分输电容量,这个容量就是容量效益裕度CBM,这里的效益指电网互联的效益。 ( H( \# [3 C. b+ L7 o7 O: P/ y: F' G; x% u8 e2 ~3 n
负荷服务公司LSE(Load Serving Entity)通过共享互联系统中其它区域的备用容量而受益,最终应该反馈给用户。在确定CBM时,需要考虑所有相关的用户,并且保证所有的LSE都能享有这一裕度带来的效益。保留的CBM减少了发电机备用容量,只有在紧急发电不足时才可以被使用。7 l$ p+ Y4 w& g7 V- l0 _. g
) q, I& Z* {$ Y S+ v+ a8 d 目前计算CBM没有很好的方法,一般取为系统内最大发电单元出力的一个倍数,或者取为TTC的一个固定的百分比,这两种方法可统称为确定性方法。文献[7]引用发电可靠性指数LOLE(Loss ofLoad Expected)来计算CBM的值。 : E7 _8 K% V3 ?8 m4 W$ I! X7 n + k3 }# _3 Y9 ?$ r( T! J 一般假设发电可靠性指数为LOLE<1天/10年,或LOLE<2.4小时/年,也就是说在10年中,负荷超过最大发电容量的时间累计最多不超过1天。如果某区域自身不能满足该标准,就需要预留输电裕度以便从其它区域得到功率。这些需要的功率在其它区域的分配是按照其它各区域的LOLE指数相反的比例进行的,即LOLE指数的区域将提供更多的功率。显然当区域内部发电备用充裕时所需CBM可能为零。 / E9 S1 x& g1 {2 r $ |! l- I4 ~; e$ _' P7 l) O 无论用哪一种方法,输电提供商对所有的用户服务公司的要求准则必须同等。而计算方法的选取往往取决于时间跨度的大小。举例来说,研究很近时间断面的情况,发电机强迫停运或保养停运的不确定度很低,确定性的方法就比较适用。而在研究长时间跨度的情况时,由于不确定度的原因,随机方法更适用。 : s- w" y. R: K6 y: E 1 N, ~! S" u! J! R# @ 通常来说,CBM不需要实时计算,一个比较长的时段里的CBM可以取其中的最大值。) B8 M/ C0 D8 ]" k$ R. B
2 Z9 _- V0 Q6 W' K$ O% b
3 ATC的计算方法- R' `3 G; k6 p# J
: Q. a9 [+ H* v ATC计算方法大体上可以分为确定性的求解方法和基于概率的求解方法两大类。确定性的方法主要包括线性分布因子法(LDF:Linear Distribution Fac-tor)、重复潮流法(RPF:Repeated Power Flow)、连续潮流法(CPF:Continuation Power Flow)、最优潮流法(OPF:Optimal Power Flow)、灵敏度分析法等。另外不少文献采用概率性的方法计算ATC时,主要基于很多电力系统参数所具有的各种随机特性,因而期望在概率框架下能得到与实际情况更为接近的解。事实上,一方面这些随机不确定因素在TRM中已被计及,另一方面在求解ATC时基于概率的方法与确定性方法的基本原则是一致的,因此这里只评述确定性的求解方法。 8 o' L1 M, k( t% ]9 v0 E; v ; ?6 e, X5 L7 L, \* O3.1 线性分布因子法) {9 a E4 ?4 y0 {! x2 r7 v
' ^& o3 ?, K0 S% _ L- p/ W
线性分布因子法也叫直流灵敏度系数法,是基于直流潮流分析实际网络响应系数的方法,一般用到多种线性分布因子[8][9]。( b% @, j+ c/ G0 M1 K2 ~) M
( ?3 ?2 k. A8 |! v
直流潮流假设节点电压幅值为常数,计及支路电抗而忽略支路电阻,因而不存在线损。直流潮流模型是线性的,不需要迭代,因而计算速度快,目前在电力系统各个领域中得到广泛的应用。& f1 ]1 w: G; O$ B
5 a- A. _5 W1 h3 ^ 文献[10]介绍了支路停运分布因子(LODF)、功率传输分布因子(PTDF)和发电机停运分布因子(GODF),给出了在基准状态、支路停运和发电机停运情况下的ATC计算模型。其中,LODF描述了当电网中发生单条支路停运时其它支路上有功潮流的变化;PTDF描述了在指定的送受端间多传输单位有功功率时各支路潮流的变化;GODF描述了某一发电机停运后,系统各支路有功潮流的变化。这样在基准状态、单支路停运和单发电机停运下各条支路的有功潮流的增量都与假想的功率传输的增量成线性关系。在已知各条支路过负荷极限的情况下,可以方便地计算最大的输送功率增量。这个数量实际上是TRM、CBM与ATC之和。. x0 d: a8 G4 ], x
5 k U; Q; m' O
线性分布因子法能很方便地考虑“N-1”静态安全约束和支路过负荷约束。在计算过程中不需迭代,求解速度快,可以满足在线应用。但无法计及电压约束和其它稳定约束,并且由于忽略电压和无功因素,在电网结构不紧密、无功支持不充足的系统中将存在较大误差。文献[2]中提及的网络响应法(Network Response Method)和额定系统路径法(RatedSystemPath Method)计算ATC也可以看成是属于线性分布因子的计算方法。 8 w* }- m4 c: h 9 I) ~$ \: X5 N( c f! s+ n3.2 重复潮流计算法 3 D5 t+ D) x+ I h ; G( f. A2 P* f8 p) A J 重复潮流法又叫常规潮流法[7]。这种方法基于常规交流潮流,计算中考虑节点电压限值约束,支路过负荷约束以及其它可能的稳定约束。其要点是逐渐增加负荷侧的负荷,同时相应增加发电侧的出力,直到某一约束生效为止,此时通过所研究断面的有功潮流之和即为最大输电能力TTC。 2 L/ K( M* j0 O" Y8 `3 P% }4 t: F9 P# L- ]/ R
典型的电力系统交流潮流方程可以表示为:4 B0 w: p" ~" r$ N
" G u! j% i- k0 v. V% j. A
9 D. }6 J. P4 |. D" w, ~2 k 3 b) N1 C+ L+ ~% ?7 o 式中x=[V,θ]。取步长为ΔP(ΔQ),逐步增加P(Q)的值,重复求解状态变量x,检查是否越界,并计算出各条支路的潮流,检验是否满足过负荷要求。若满足,继续增加步长,直到出现某一越限时,将步长减半,重复上述过程,直至最后步长满足误差要求为止。& Z$ ]) u, D2 p
. K R* Y! T- p$ i5 U 为了加快计算速度,可以对重复潮流计算做一定的改进[11]。首先采用线性分布因子法,逐步增加受电侧的负荷和发电侧的出力,直至有支路过负荷为止。然后采用交流潮流计算各节点电压,检查是否有电压越限的现象发生,如果有,则逐步减少受电侧的负荷和发电侧的出力,直到满足所有节点电压限值约束为止。- W, F# x8 R- b2 t$ O0 g6 [
$ ?$ A; K9 z' b' ?& z/ N9 X7 ] 重复潮流计算ATC的方法原理简单,可以计及系统的电压和无功对ATC的影响,计算结果能较好的反映实际运行状况,但需要重复计算交流潮流,计算时间长,不适合在线应用。 - K! w( Q" \ s! b: A " o+ S- E$ N+ ~7 p! y& r3.3 连续潮流法 $ y/ M* g# \( s+ k2 u ) G0 r; Y5 L0 B( ^- S% z! o 作为一种求解非线性代数方程的数值方法,连续潮流法在上个世纪90年代引起人们关注,是因为该方法在电压稳定性研究方面有其独特的优越性。CPF广泛应用于计算静态电压稳定的P-V曲线中的极限功率点(NP:Nose-Point)[12]。由于牛顿法在电压稳定极限点附近因雅可比矩阵奇异,引起潮流不收敛。而CPF法可从当前潮流解出发,逐步增加指定送端母线功率,通过迭代求解,沿P-V曲线准确得到NP点相应的发电功率,因而可以被方便地用来计算静态电压稳定约束下的ATC。( Z4 O% D4 Y2 R& i! ^6 x# i
8 x, o* O1 A, `. }+ [4 q
当系统中发电机功率或负荷发生缓慢变化,如果用P0和Q0表示对应于系统当前状态下的节点有功和无功向量,则可将系统方程参数化为, 0 Y) N" P( o, @5 B ( e) V4 w, U4 G4 z* o
, c. \7 i, }9 {* Y+ _ 其中,u为控制变量,包括发电机有功输出,机端电压,变压器变比等;x为状态变量,包括节点电压幅值和相角;等式约束g(x,u)=0为潮流方程;h(x,u)≤0表示所有不等式约束。 6 T! G: t: p/ p 9 g3 q/ T5 N& _" f8 y! u 优化潮流方法的关键瓶颈在于计算速度。在实际使用中,往往只能作为离线计算工具。当对计算速度要求较高时,可以采用简化的直流潮流方程作为等式约束。但若要考虑“N-1”静态安全约束,将使问题的求解规模大大增加。文献[16]采用的Benders分解法有效地解决了这一问题,Benders分解法将问题分为主从两层,主问题处理基态潮流及相应约束,而每一个预想事故则形成一个子问题,每个子问题可单独求解,其起作用的约束以Benders割集的形式返回主问题,主子问题反复迭代直至全部约束满足,求得最优解。$ `7 o F/ e2 f. A% @1 q7 c
/ D' Y; j& q0 {2 O: m6 ~" X: M; Y 基于最优潮流的ATC计算方法同连续潮流法相比,对约束条件有更强的处理能力,理论上可以处理各种约束,还可以进行有功和无功优化,计算结果更准确。但这类方法有下述问题:很难考虑系统稳定性这样的动态约束条件;所获得的最优运行点是一个理想的结果,实际上难以达到;需要的计算时间长,难以满足在线要求;目前尚难应用于超大规模电力系统。1 C1 `/ E8 S" O6 y; f. L) U
* Z$ Q1 u0 e- d1 q& O/ K6 t! _ 美国电力科学研究院EPRI开发的可用来计算ATC的商用软件TRACE(Transfer Capability Evalu-ation)就是基于优化潮流方法的,目前已经得到广泛的应用。但计算规模有限,无法计及整个互联大区的相互影响,只能用于离线的分析。因此NERC建议各大区电网采用北美整个电网的结构,使用线性分布因子法,使得计算的ATC能有效的跟踪电网和市场的变化。' G7 E0 p. a) E
* o/ `9 G1 o( Z3 }3.5 灵敏度分析法6 |; B0 V1 p, m) C/ H
+ t+ P: y8 ]6 T" M1 R/ y) o 由以上分析可知,线性分布因子法计算速度快,但精度不高;重复潮流法、连续潮流法和最优潮流法计算时间长,难以在线应用,而且重复潮流法和连续潮流法计算结果往往偏保守。灵敏度分析法就是在这样的背景下被提出来的,其计算速度快,满足在线计算要求,同时能保证一定的精度。但是需要指出,它不是一个独立的ATC计算方法,它需要以某种ATC计算结果为基础[11]。8 X0 {" e" V' p& `) s, Y( d
( i* S* \. H) Q6 H& U1 G P& [
基于交流潮流的灵敏度分析法的最大优点是当电力系统中某些运行参数发生变化后,它可以快速计算出其对ATC的影响,而不需要重新进行潮流计算,从而在系统运行状态改变后,获得特定断面的ATC。该方法从某一运行点下的已知的ATC值出发,只分析当系统参数在此基础上发生微小变化时对ATC值的影响。ATC对某参数变化的灵敏度可以是一阶的,也可以是二阶或高阶的,实际使用中,一阶模型已经足够。显然该方法的快速性是以牺牲一定的准确度为代价的,因为一阶灵敏度系数只能反映各变化因素与ATC之间的线性关系,不能计及它们之间的非线性关系。 3 i! l0 L* H+ l3 y$ t$ p0 _1 L, H9 a( Y x3 B: B
对于一组给定的系统参数,可以利用重复潮流法、连续潮流法或最优潮流法计算某一断面的ATC值T0,并计算ATC值对各种可能变化参数P的一阶灵敏度系数Tp。当P发生变化时ΔP,新的ATC值可以由下式直接得到[19]: 9 R- q9 |6 y- t, V : D3 V, H5 F3 ? k3 I( T; j+ K: i
6 b% N }1 F; B" U* t7 Q2 U0 A9 g+ p 基于一阶灵敏度分析法计算ATC速度快,满足在线应用的要求。不过,对于给定的系统运行状态和参数,需要预先计算较多的灵敏度系数TP,且当运行点发生较大变化后,这些灵敏度系数需要重新计算。在系统运行情况变化不大时,求得的ATC还是比较准确的,但当系统运行情况发生较大变化时,如支路或发电机停运,可能会存在较大误差。但总的说来,灵敏度分析法与前述某种方法,如OPF法,结合在一起使用是一种非常实用的综合在线计算方案。 & O2 {; d! q. B2 t4 p# v2 {& R0 o3 |- R; c: C8 W9 H* ], S
4 结论 $ z3 J% {4 z4 o; s/ C * t& r, b8 \/ [* N7 f0 p0 p4 d 本文介绍了电网可用输电能力的计算框架,解 释并给出了TRM和CBM这两种裕度的内容和计算方法。分析了计算ATC的各种算法,分析比较其优缺点,得出如下结论: _4 _3 C: P8 l- m2 V7 u * o( O2 }* V/ E. ^7 {% W6 X a.电网可用传输能力是反映系统安全性能的一个重要指标,可以为电力市场参与者提供电网的使用情况,是电力市场正常运作的基础,因此要求ATC的计算结果准确、快速、全面。 + S% U8 |$ D) }9 D; Y ) S! e9 \( U- [ b.基于线性分布因子法的ATC计算速度快,可有效考虑输电设备过负荷约束和“N-1”静态安全约束。 9 K$ j' m( |; Q9 l+ `/ L, r, t/ c0 X! |7 F* v
c.基于重复潮流和连续潮流的ATC计算可以方便考虑系统设备过负荷约束,能计及系统的电压和无功对ATC的影响。但计算结果偏保守,计算时间较长,考虑“N-1”静态安全约束时无法满足实际在线计算要求,可用于离线的ATC计算。( u/ m/ V5 p, C0 J: c! ^3 \
+ f; V5 m* A4 c% P( s! y% d
d.基于最优潮流的ATC计算方法对约束条件有更强的处理能力,能够计及暂态稳定和动态稳定约束,可以实现更大范围内的发电和负荷分布优化,计算精度更高;但计算耗时更长,无法满足实际在线计算要求,是一种很好的ATC离线计算工具。 M! S% g z u
- V! Q; L* {8 H; ^ e.基于灵敏度分析的ATC计算速度快,计算中无需任何迭代。与其它精度较高的ATC计算方法配合使用,能在系统状态发生变化时快速得到新的ATC,特别适合应用在系统运行状态变化不大时的ATC计算,但当系统参数变化较大时精度稍差。1 z7 ^2 u1 Z, g
! u# n$ {8 I: p( q3 [0 U5 p! U6 e
我国电网已形成东北、华北、华中、西北、华东和南方六大区域电网。计算各大区域之间的ATC有很大的实际意义。在我国电力市场运行初期,大区间联网网架比较薄弱,建议先使用线性分布因子法实时计算ATC,并可采用重复潮流法或连续潮流法来离线效验。随着市场的成熟及基于最优潮流的ATC算法的发展和完善,再逐步过渡为用最优潮流加上灵敏度分析方法,将会提供实时准确的ATC。2 [9 A6 v- E, M+ z
4 m/ u9 t* u. T# N& X2 d7 e8 R作者: 6 {8 r \) _1 Z3 u" [3 ?5 G刘皓明1,倪以信2,吴军基1,邹云1 % ~5 n. B( ^; W+ x% ^# G5 V
(1.南京理工大学动力学院,江苏南京210094; " M: S. ~* ^9 h* F9 ~* E" s2.香港大学电机电子工程系香港); {' h+ S& N, s
+ @9 {; K/ d6 {6 R6 `; \" ^参考文献; \7 k P( B% T
% _& y; f9 Y; Y
[1]李国庆,王成山,余贻鑫.大型互联电力系统区域间功率交换能力研究综述[J].中国电机工程学报,2001,21(4):2025. ! V C$ n- V/ J% I0 {9 M$ K9 ~[2]Transmission Transfer Capability Task Force.AvailableTransmission Capability Definitions and Determination.NorthAmerican Electric Reliability Council,Princeton,NewJer-sey,June 1996.; v3 w. Y. ^& |. B& P1 G5 i
[3]Federal Energy Regulatory Commission.Open Access Same-Time Information Systemand Standards of Conduct.Docket No.RM95-9-000,Order 889,April 1996.; K7 s: I: H3 ]& d! ]9 ~1 L$ @7 `$ P
[4]Transmission Capability Marginsand Their Usein ATCDeter-mination-White Paper.Prepared by the North AmericanElectric Reliability Council Available Transfer Capability Working Group.June,1999. ; E$ p; Q4 P j& W7 S, L, ][5]P.W.Sauer.Alternativesfor Calculating Transmission Reli-ability Margin(TRM)in Available Transfer Capability(ATC).Proceedingsofthe Thirty-First HawaiiInternational Conference on SystemSciences,1998(3).( C$ t! M1 e% y, _: r9 v
[6]I.Dobson,etc.Electric powertransfer capability:concepts,applications,sensitivity,uncertainty.Power Systems Engi-neering Research Center publication 00-34,November 2001.availablefromhttp://www.pserc.wisc.edu/. 4 _- C3 u& d4 B/ a[7]Ou.Yan,C.Singh.Assessmentof Available Transfer Capa-bility and Margins.IEEE Transactions on Power Systems,2002,17(2):463~468.4 K4 J- ]& N) G9 X9 G4 r2 \2 ^
[8]M.K.Enns,J.J.Quada,B.Sackett.Fast Linear Contin-gency Analysis.IEEETransactions on Power Apparatus and System,1982,PAS-101(4):783~791. & J1 a' S2 T' N% u4 v! @3 h7 y[9]Peter W.Sauer,Karl E.Reinhard,Thomas J.Overbye.Extended Factors for Linear Contingency Analysis.Proceed-ings ofthe 34th Hawaiiinternational conference on systemsci-ences-2001:697~703.. {. e( X4 z) W; l! `! s- _4 H. g) g
[10]Ejebe G.C.,WaightJ.G.,Sanots-Nieto M.,Tinney W.F.Fast Calculation of Linear Available Transfer Capability. B* D7 @7 t1 d2 z0 P3 U$ C0 G/ q
IEEETransactions on Power Systems,2000,15(3):1112~1116. ) }3 z) C+ _# C) d2 @( t i X# C[11]M.H.Gravener,C.Nwankpa.Availabletransfer capabilityand first order sensitivity.IEEE Transactions on Power Sys-tems,1999,14(2):512~518. 8 b3 ?# }5 r" G) N: {[12]H.D.Chiang,A.Flueck,K.S.Shah,N.Balu.CP-FLOW:A practical tool for tracing power system steady-state stationary behavior due to load and generation varia-tions.IEEETransactions on Power Systems.1995,10(2):623~634. 1 e( Z. a, |) @7 z' e$ O ~[13]A.J.Flueck,H.D.Chiang,K.S.Shah.Investigatingtheinstalled real power transfer capability of a large scale powersystemunder a proposed multiareainterchange schedule usingCPFLOW.IEEE Transactions on Power Systems,1996,11(2):883~889. 4 b$ Q, p- L, z0 R[14]丁晓莺,王锡凡.最优潮流在电力市场环境下的最新发展[J].电力系统自动化,2002,26(13):1~7.6 _9 U4 ^5 [# {! e* J% m
[15]刁勤华,Mohamed Shaaban,倪以信.运用连续二次规划法计算区域间极限传输容量[J].电力系统自动化,2000,24(24):5~8.) W* [* ^) e9 D6 j9 O
[16]Youjie Dai,J.D McCalley,V.Vittal.Simplification,Ex-pansion and Enhancement of Direct Interior Point Algorithm for Power System Maximum Loadability.IEEE Transactions on Power Systems,2000,15(3):1014~1021. * z+ P, f* T+ f! S: Z[17]X.Luo,A.D.Patton,C.Singh.Real powertransfer ca-pability calculations using multi-layer feed-forward neuralnetworks.IEEE Transactions on Power Systems,2000,15(2):903~908.& A9 D W8 l# b' e# F1 n
[18]焦连伟,Mohamed Shaaban,倪以信.可用传输容量计算的分层算法研究[J].中国电力,2002,35(1):33~36.6 B: g. W. c2 g0 G; A
[19]S.Greene,I.Dobson,F.L.Alvarado.Sensitivity ofTransfer Capability Margins with a Fast Formula.IEEETransactions on Power Systems,2002,17(1):34~40.-%-9142-%--%-10637-%--%-10622-%-