摘要:在电力市场环境下,电网可用输电能力(ATC)是反映输电系统可用于传输电力的剩余容量的重要指标。文中根据国际公认的NERC关于ATC的定义,讨论了TRM和CBM两种裕度的含义及几种计算方法;详述了现有的多种ATC的计算方法,包括线性分布因子法、重复潮流法、连续潮流法、最优潮流法和灵敏度分析法,并分析比较了各种方法的优缺点。最后对我国电力市场环境下大型互联电网之间的ATC计算提出了一些建议。 ) \0 n* ^$ ?3 ]0 ^6 X2 p0 y 4 ]0 u" @7 }3 }+ L' R) @1 引言6 e7 }7 o; L0 J' P) M
7 z4 z% d& I+ R0 H) `
电力系统区域间功率输送能力计算的研究始于20世纪70年代,至今已有30年的历史。在传统的电力工业运行模式下,它主要用于评估系统互联强度及比较不同输电系统结构的优劣[1]。而在电力市场环境下,系统中往往存在着大量的频繁变化的电力交易,这使得输电系统负荷增加、环流增大、容量裕度降低,进而稳定裕度减少,电力系统的安全与稳定问题更加突出。为保证电力系统的安全运行,需要实时评估电力系统运行的安全水平。在这样的背景下,可用输电能力(ATC:Available Transfer Capabili-ty)的计算就显得十分重要,它不仅可以显示电网运行的安全与稳定裕度,从而减少阻塞发生的概率,也可以为市场参与者提供电网使用状况的详细信息,以指导他们参与市场的行为。( E0 C4 \. E" d( |) F9 w9 V N
' i: t1 S& z% }: N4 O# R/ b/ b% N1 d
北美电力可靠性委员会NERC(North AmericanElectric Reliability Council)在上个世纪九十年代统一了有关输电极限的概念,提出了可用输电能力ATC的详细定义与计算框架[2],这一工作在国际上已得到了更为广泛的认可。ATC是在已有的协议基础上,在实际输电系统中可以用于进一步商业活动的富余输电能力。数学上,ATC可以表示为线路最大输电能力TTC(Total Transfer Capability)减去输电可靠性裕度TRM(Transmission Reliability Margin),减去容量效益裕度CBM(Capacity Benefit Margin),再减去现存输送协议ETC(Existing Transmission Commit-ments)。ATC的计算值应实时发布在电力市场公共信息系统网页上,以进一步促进输电网络的开放,进而促进电力市场的交易[3]。7 I$ P5 N: d7 Q; b" c7 k
! M# m/ m: ]8 S5 r" ~) k
按照NERC建议的计算框架,ATC的计算需要考虑三种物理约束:电压约束;设备过负荷约束;包括静态、动态和暂态等在内的各种稳定极限约束[2]。这里静态安全约束为“N-1”准则。0 I; l9 e0 M2 A$ T1 H, ^
5 X( m) H% F |. m {7 Z
计算ATC必须要计算TRM和CBM这两种输电容量裕度。本文将先分析TRM和CBM两种裕度的含义,介绍几种计算两种裕度的方法;然后详述了现有的多种ATC的计算方法,包括线性分布因子法、重复潮流法、连续潮流法、最优潮流法和灵敏度分析法等,分析比较了这些方法的优缺点。在文章最后对我国电力市场环境下大型互联电网的ATC的计算提出了一些的建议。" X+ k; B& h! I0 L9 c
7 l3 ^9 e) k5 w4 K5 ~ ~9 \
2 两种容量裕度1 W6 s8 e3 `, q
/ y7 o, ]* j$ @( E5 V! n% t
计算ATC时,一般以一个给定的基本运行状态为基础,计算可进一步用以传输电能的最大裕度,在裕度中减去TRM和CBM,即得ATC值。目前大多数文献在计算ATC时都忽略这两个裕度是不正确的。本节将详细的解释它们的含义,并介绍几种可用的计算方法[4]。 1 b* l/ h. O2 Q6 K: a ( Q# ?& @3 U3 D9 V# @2.1 输电可靠性裕度TRM[5]0 X3 v1 Z; i" s5 _4 a- _' j# W
6 P0 m8 ^- ]3 P9 e
TRM指的是预留的必要的电网输电能力,以确保当系统运行参数在合理范围内发生变化时,整个系统能够安全稳定的运行。" C2 ~( ^) f" I
! h( _3 `' g0 C: Q
这些不确定的运行参数可能包括支路停运、发电调度、负荷预测、并行潮流等。一般来说,时间跨度越大,它们的不确定性越大。因此,TRM与考虑的时间断面有关,时间跨度越大,往往需要预留的TRM越多。5 W5 o1 ~) h+ V' {# ~, H- M
! V$ A& f. h$ ^( ?
计算TRM时需要全面考虑到各种不确定因素,并将它们合理地组合到一起,而不是仅仅代数叠加。TRM的计算方法主要包括:% E# b& K( X: |0 H% _$ c* t
- d1 k) M) V/ p6 x/ p a.在基准情况下不断改变假设和预测的参量,重复计算TTC,所得到的TTC中最大值与最小值之差即为TRM。理论上需要取遍所有可能的参数变化组合。7 X0 {& {' x, Q; N
3 C6 N6 w1 N3 G8 T9 \/ P* f% D
b.按经验取TTC的一个固定百分比(比如4%)。这种方法简单方便,较为常用。 ( @* Y' l4 j( g0 r& `; q' k# {' t2 _+ A
c.利用统计学或概率方法求解[6]。这类方法的计算过程较为繁琐,也不太成熟。$ v& C6 l+ Q/ d3 C H- g2 u0 |/ W/ n
; ^& Y5 c( Z X) e$ ^& U( e: S
d.降低额定值法:系统中不确定因素对所有电力设备的影响是相对一致的,可以通过降低设备的计算额定值来计及这些不确定因素。典型的额定值降低为2-5%,当然降低的幅值会随着时间跨度的增大而增大[4]。, e) G9 X: o& N3 l
5 V( x' T4 q* d" [. s. j/ M9 a 在实际应用中,并不需要实时计算TRM值,欧共体内部电力市场只要求每半年刷新一次TRM值。计算方法更多是取TTC的一个固定百分比,对于电气连接紧密,无功支持充裕的电网,TRM所占TTC的百分比可以取的小一些。3 D* V3 q. @% E- F
) P1 b0 h( I3 j- q6 }# R$ i$ B2.2 容量效益裕度CBM + [1 G7 A5 c m" B o- X4 L4 T' Q0 b7 f/ p' ^% l
CBM指的是:为了能够从其它互联系统中获得电力来满足发电可靠性需求而预留的输电容量裕度。' y* `5 K5 U+ J9 c
0 s/ P" _9 N. |. e1 b; D4 l! b
备用发电容量对保证供电可靠性至关重要,在发电机停运或因其它设备故障而失去部分电源时,需要启动备用电源以保证向负荷供电。对于互联的电力系统,当发生故障引起电源短缺时可以从其它子系统得到紧急援助,因而每个系统内部都可以减少发电备用容量,以提高经济性。为了保证系统在任何时刻都有能力将系统外的电能输送到负荷中心,需要预留部分输电容量,这个容量就是容量效益裕度CBM,这里的效益指电网互联的效益。 ( }, w, z8 ~+ G: Y9 y- s + l. J' Q3 @7 D2 S 负荷服务公司LSE(Load Serving Entity)通过共享互联系统中其它区域的备用容量而受益,最终应该反馈给用户。在确定CBM时,需要考虑所有相关的用户,并且保证所有的LSE都能享有这一裕度带来的效益。保留的CBM减少了发电机备用容量,只有在紧急发电不足时才可以被使用。# c2 N/ |. E) G
! f1 H* f0 b, `7 ~, `
目前计算CBM没有很好的方法,一般取为系统内最大发电单元出力的一个倍数,或者取为TTC的一个固定的百分比,这两种方法可统称为确定性方法。文献[7]引用发电可靠性指数LOLE(Loss ofLoad Expected)来计算CBM的值。 / N! T; r$ d y7 s' w+ ~/ F* f! G 9 O3 O$ d/ \! ]5 p# ~7 ^ 一般假设发电可靠性指数为LOLE<1天/10年,或LOLE<2.4小时/年,也就是说在10年中,负荷超过最大发电容量的时间累计最多不超过1天。如果某区域自身不能满足该标准,就需要预留输电裕度以便从其它区域得到功率。这些需要的功率在其它区域的分配是按照其它各区域的LOLE指数相反的比例进行的,即LOLE指数的区域将提供更多的功率。显然当区域内部发电备用充裕时所需CBM可能为零。0 @6 z) L' O- j( Z6 q& b
0 F% Q' B- S" [- O- Q
无论用哪一种方法,输电提供商对所有的用户服务公司的要求准则必须同等。而计算方法的选取往往取决于时间跨度的大小。举例来说,研究很近时间断面的情况,发电机强迫停运或保养停运的不确定度很低,确定性的方法就比较适用。而在研究长时间跨度的情况时,由于不确定度的原因,随机方法更适用。2 Y! o# n* J/ M# h
5 q2 i" A% A0 h% A% h. n! `
通常来说,CBM不需要实时计算,一个比较长的时段里的CBM可以取其中的最大值。 , r8 w+ L2 W. g8 e. P U - C }) z* l4 \6 P* h3 ATC的计算方法 $ m, I) p* N0 r( f+ ?2 N" S; r' Z. s: I* N" L! }5 D' K+ _
ATC计算方法大体上可以分为确定性的求解方法和基于概率的求解方法两大类。确定性的方法主要包括线性分布因子法(LDF:Linear Distribution Fac-tor)、重复潮流法(RPF:Repeated Power Flow)、连续潮流法(CPF:Continuation Power Flow)、最优潮流法(OPF:Optimal Power Flow)、灵敏度分析法等。另外不少文献采用概率性的方法计算ATC时,主要基于很多电力系统参数所具有的各种随机特性,因而期望在概率框架下能得到与实际情况更为接近的解。事实上,一方面这些随机不确定因素在TRM中已被计及,另一方面在求解ATC时基于概率的方法与确定性方法的基本原则是一致的,因此这里只评述确定性的求解方法。 V# i$ P) |2 x. |- W2 { % }# H3 S4 @; B* J& K' J3.1 线性分布因子法; [, C. [ {, I8 U! Z. v
6 ^* i! E; c ^/ T {0 b" U 线性分布因子法也叫直流灵敏度系数法,是基于直流潮流分析实际网络响应系数的方法,一般用到多种线性分布因子[8][9]。 5 Y2 I. u/ e& c) g $ Q. u0 } o2 \0 x 直流潮流假设节点电压幅值为常数,计及支路电抗而忽略支路电阻,因而不存在线损。直流潮流模型是线性的,不需要迭代,因而计算速度快,目前在电力系统各个领域中得到广泛的应用。+ l2 G' q; ~1 W% c7 P. s0 h
% c5 N8 m/ I- i- \4 [ p+ z3 p* m% u 文献[10]介绍了支路停运分布因子(LODF)、功率传输分布因子(PTDF)和发电机停运分布因子(GODF),给出了在基准状态、支路停运和发电机停运情况下的ATC计算模型。其中,LODF描述了当电网中发生单条支路停运时其它支路上有功潮流的变化;PTDF描述了在指定的送受端间多传输单位有功功率时各支路潮流的变化;GODF描述了某一发电机停运后,系统各支路有功潮流的变化。这样在基准状态、单支路停运和单发电机停运下各条支路的有功潮流的增量都与假想的功率传输的增量成线性关系。在已知各条支路过负荷极限的情况下,可以方便地计算最大的输送功率增量。这个数量实际上是TRM、CBM与ATC之和。 , Z3 _* C' L! N & g: |! X4 S, N6 q- c' K! e 线性分布因子法能很方便地考虑“N-1”静态安全约束和支路过负荷约束。在计算过程中不需迭代,求解速度快,可以满足在线应用。但无法计及电压约束和其它稳定约束,并且由于忽略电压和无功因素,在电网结构不紧密、无功支持不充足的系统中将存在较大误差。文献[2]中提及的网络响应法(Network Response Method)和额定系统路径法(RatedSystemPath Method)计算ATC也可以看成是属于线性分布因子的计算方法。9 f T/ j' n4 M+ X. _ j3 T, H p
% [5 j/ ?$ l; X) f, N$ u3.2 重复潮流计算法 , W+ G' n: l; X6 G: m " @4 a2 A$ y. b% Y! M7 j) f" f6 Z 重复潮流法又叫常规潮流法[7]。这种方法基于常规交流潮流,计算中考虑节点电压限值约束,支路过负荷约束以及其它可能的稳定约束。其要点是逐渐增加负荷侧的负荷,同时相应增加发电侧的出力,直到某一约束生效为止,此时通过所研究断面的有功潮流之和即为最大输电能力TTC。 $ F2 w# {+ a! m. ~/ R: C# l' k i$ a; r5 R3 f8 Y" c% u
典型的电力系统交流潮流方程可以表示为: 4 [0 n8 p2 x$ ? * e6 {7 D% T# j% D8 M: R3 m6 V! C+ ^
4 C' s/ x. e2 h" ` n) c
. H2 z, g4 G; t( i' P( `9 u
其中,u为控制变量,包括发电机有功输出,机端电压,变压器变比等;x为状态变量,包括节点电压幅值和相角;等式约束g(x,u)=0为潮流方程;h(x,u)≤0表示所有不等式约束。+ X& g% |* V7 v
`: D; W6 L1 x Z# U 优化潮流方法的关键瓶颈在于计算速度。在实际使用中,往往只能作为离线计算工具。当对计算速度要求较高时,可以采用简化的直流潮流方程作为等式约束。但若要考虑“N-1”静态安全约束,将使问题的求解规模大大增加。文献[16]采用的Benders分解法有效地解决了这一问题,Benders分解法将问题分为主从两层,主问题处理基态潮流及相应约束,而每一个预想事故则形成一个子问题,每个子问题可单独求解,其起作用的约束以Benders割集的形式返回主问题,主子问题反复迭代直至全部约束满足,求得最优解。! ~. k5 E3 S. C; ]4 }+ s8 a
# k9 n( l$ x0 \1 V
基于最优潮流的ATC计算方法同连续潮流法相比,对约束条件有更强的处理能力,理论上可以处理各种约束,还可以进行有功和无功优化,计算结果更准确。但这类方法有下述问题:很难考虑系统稳定性这样的动态约束条件;所获得的最优运行点是一个理想的结果,实际上难以达到;需要的计算时间长,难以满足在线要求;目前尚难应用于超大规模电力系统。4 g" Y& d" I3 e% U- a) d
7 X' A' `8 n# C; e, W- v( |; Y 基于一阶灵敏度分析法计算ATC速度快,满足在线应用的要求。不过,对于给定的系统运行状态和参数,需要预先计算较多的灵敏度系数TP,且当运行点发生较大变化后,这些灵敏度系数需要重新计算。在系统运行情况变化不大时,求得的ATC还是比较准确的,但当系统运行情况发生较大变化时,如支路或发电机停运,可能会存在较大误差。但总的说来,灵敏度分析法与前述某种方法,如OPF法,结合在一起使用是一种非常实用的综合在线计算方案。 / P* y' m, C# e) A |2 c0 k! r$ l, y9 z, T9 L4 Y% e* R
4 结论0 M& l/ X/ p* O5 W
! K' ]# v( ?# W& ~( |/ z
本文介绍了电网可用输电能力的计算框架,解 释并给出了TRM和CBM这两种裕度的内容和计算方法。分析了计算ATC的各种算法,分析比较其优缺点,得出如下结论: / t! v# Q9 W: e; x1 z3 O: |$ L * C8 `/ N W) Y+ a |* V a.电网可用传输能力是反映系统安全性能的一个重要指标,可以为电力市场参与者提供电网的使用情况,是电力市场正常运作的基础,因此要求ATC的计算结果准确、快速、全面。 + ?( n0 C+ s" l4 m3 o: T ; D0 `" {$ I8 s# K% w" ~$ ^( A b.基于线性分布因子法的ATC计算速度快,可有效考虑输电设备过负荷约束和“N-1”静态安全约束。4 E% ^' n" U6 O7 C3 L$ q0 B5 h
; |$ b8 o$ X/ `5 u$ _8 S
c.基于重复潮流和连续潮流的ATC计算可以方便考虑系统设备过负荷约束,能计及系统的电压和无功对ATC的影响。但计算结果偏保守,计算时间较长,考虑“N-1”静态安全约束时无法满足实际在线计算要求,可用于离线的ATC计算。 c: Y7 `2 D# c& ^- T1 s$ ?5 o5 q$ z
: V# Y9 e7 S5 R' ~ d.基于最优潮流的ATC计算方法对约束条件有更强的处理能力,能够计及暂态稳定和动态稳定约束,可以实现更大范围内的发电和负荷分布优化,计算精度更高;但计算耗时更长,无法满足实际在线计算要求,是一种很好的ATC离线计算工具。 1 I- J; e/ |' Y4 o, U 5 f$ q, c" b0 O$ ? e.基于灵敏度分析的ATC计算速度快,计算中无需任何迭代。与其它精度较高的ATC计算方法配合使用,能在系统状态发生变化时快速得到新的ATC,特别适合应用在系统运行状态变化不大时的ATC计算,但当系统参数变化较大时精度稍差。3 o5 D0 X: O+ k% A
0 R+ {7 r8 G5 K( J' r3 n
我国电网已形成东北、华北、华中、西北、华东和南方六大区域电网。计算各大区域之间的ATC有很大的实际意义。在我国电力市场运行初期,大区间联网网架比较薄弱,建议先使用线性分布因子法实时计算ATC,并可采用重复潮流法或连续潮流法来离线效验。随着市场的成熟及基于最优潮流的ATC算法的发展和完善,再逐步过渡为用最优潮流加上灵敏度分析方法,将会提供实时准确的ATC。! Y! b- D- {5 j5 d: y
! p9 D2 y4 x$ k" P$ h7 _2 g# ]0 a# E
作者: - R1 w* C; n' O1 m" u9 w7 l' G刘皓明1,倪以信2,吴军基1,邹云1 1 d" K2 G) E* ?; m
(1.南京理工大学动力学院,江苏南京210094; 9 }+ ]4 t' e+ \7 \! v4 v2.香港大学电机电子工程系香港) 7 b# F! W- B9 j$ K/ R* N5 l; o) A1 L1 r5 X; \9 y
参考文献5 H% H0 r6 ?7 U6 [! ~& `9 J
3 K8 c: ~' \# O
[1]李国庆,王成山,余贻鑫.大型互联电力系统区域间功率交换能力研究综述[J].中国电机工程学报,2001,21(4):2025.* r% c) T x3 X5 r6 |7 d5 O. \
[2]Transmission Transfer Capability Task Force.AvailableTransmission Capability Definitions and Determination.NorthAmerican Electric Reliability Council,Princeton,NewJer-sey,June 1996. , `- t3 l! J* K2 D$ V+ B* v[3]Federal Energy Regulatory Commission.Open Access Same-Time Information Systemand Standards of Conduct.Docket No.RM95-9-000,Order 889,April 1996. 8 Y4 j5 C7 s* I# U1 x[4]Transmission Capability Marginsand Their Usein ATCDeter-mination-White Paper.Prepared by the North AmericanElectric Reliability Council Available Transfer Capability Working Group.June,1999. # D8 d5 F% r! F0 a' J[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). , e# I" ^! ?6 M# y8 y2 i[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/.2 o" B$ U" _' n9 y& l8 f5 F& N
[7]Ou.Yan,C.Singh.Assessmentof Available Transfer Capa-bility and Margins.IEEE Transactions on Power Systems,2002,17(2):463~468. " @3 Z7 b+ X ]$ R0 i ^6 e% T$ ?[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. : g& {+ I1 ~- i+ _: f[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. 7 [1 ^# |2 @6 T. R! A[10]Ejebe G.C.,WaightJ.G.,Sanots-Nieto M.,Tinney W.F.Fast Calculation of Linear Available Transfer Capability.# [9 t$ ^+ x$ a! o* k7 }
IEEETransactions on Power Systems,2000,15(3):1112~1116.$ c6 i) P) ?: P2 X. d0 o
[11]M.H.Gravener,C.Nwankpa.Availabletransfer capabilityand first order sensitivity.IEEE Transactions on Power Sys-tems,1999,14(2):512~518., H, x0 z, L9 H" j: y. j$ t* W
[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.: H t+ G: P, j; ~
[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. ' E L; R9 n D9 z7 y8 H1 \[14]丁晓莺,王锡凡.最优潮流在电力市场环境下的最新发展[J].电力系统自动化,2002,26(13):1~7.) |# f1 ~/ d1 Y( P0 u. y5 Z8 h
[15]刁勤华,Mohamed Shaaban,倪以信.运用连续二次规划法计算区域间极限传输容量[J].电力系统自动化,2000,24(24):5~8. P+ X! N8 E+ I! @
[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. 4 t7 f! r* L( A9 F% ^[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. * z% |4 R" R, f( B* P[18]焦连伟,Mohamed Shaaban,倪以信.可用传输容量计算的分层算法研究[J].中国电力,2002,35(1):33~36.; B8 W# ?2 I& L2 b8 {
[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-%-