1. 引言
在全球气候变化和城市热岛效应的背景下,夏季极端高温天气越来越普遍。例如,2016年夏季(6~8月),重庆全市平均气温较常年同期偏高0.9℃,35℃以上的高温日数较常年同期偏多11.8 d,共出现6次区域高温过程,其中8月10~26日出现了连续17天的区域高温过程,综合强度达特重等级 [1] 。夏季极端高温天气下空调用电量大幅度增长,电网负荷屡创新高,据相关数据显示,8月25日重庆电网统调最高负荷创历史新高,达到1815万千瓦,同比增长37.1% [2] 。期间重庆某供电公司出现多条220 kV断面部分时段不满足N−1安全准则,7条110 kV线路重载或过载,36座110 kV变电站的负载率高于80%。连续极端高温天气下,空调电负荷大,而高温天气下线路载流能力下降,因此架空线路的运行安全问题备受电网调度人员关注。
目前广泛使用的输电线路静态热极限值,是按历史气象条件重现期内的最恶劣气象参数(例如40℃的气温,0.5 m/s的风速)求取的保守数值 [3] [4] [5] ,利用线路的静态热极限值无法提高夏季连续高温天气下架空线路的载流能力 [6] 。为尽可能提高架空线路的电能输送能力,应尽可能提高设备利用率,故可在不进行电网结构改造的情况下利用架空线路的动态热容量实现电能输送的效能最大化 [7] [8] [9] [10] 。目前已开展了大量关于架空线路动态增容的研究,文献 [11] 提出利用输电线路的动态载流量代替静态载流量,利用在线热定值测量系统实时动态地计算输电线路的输送容量:输入输电线路的实际环境参数,即环境温度以及风速风向等参数,利用输电线暂态热平衡方程计算输电线路的输送容量。结果表明,输电线路的动态热容量相比于静态热容量,可以大幅提高输电线路的输电能力。针对在线热定值测量系统设备成本太高,难以广泛应用的缺点,文献 [12] 提出利用参数估计的方法确定实时气象环境参数,动态计算线路的最大输送容量。文献 [13] 研究了不同气候因素对架空线路动态载流量的影响程度,仿真分析发现风速对架空线路的载流量影响程度最大,日照强度对架空线路的载流量影响最小,但其仅针对不同气候因素的影响进行了定性分析,并未进行定量分析。而文献 [14] 提出了一种架空线路能力在线评估的简化方法,通过线路温度的在线测量以及利用短时气候预报信息,采用辨识修正气象因子,实时获取线路的动态载流能力,在线预测架空线路的最大载流能力。
综上可知,基于实时气象环境参数的动态载流量,实现了架空线路的有效增容,有效地提升了线路的输送能力,但是并未在连续高温天气下实时兼顾线路的动态载流量和运行温度,缺少架空线路动态载流能力的评估与运行温度的安全告警。因此,本文利用实时气象环境参数,借助导线载流的热平衡方程,在线估计导线运行温度以及动态载流量,动态辨识架空线路的安全裕度,分析导线在供电断面N−1潮流转移后的温升响应时间。最后以重庆地区某供电断面为算例,对所提方法的有效性和准确性进行验证。评估结果可为连续高温天气下的电网运行调度提供指导。
2. 夏季连续高温天气及负荷特点
重庆地形以山地为主,地势由西向东逐渐升高,由南向北长江河谷倾斜,起伏较大,四周山地环绕,地面散热困难,夏季易形成高温区。其中沿江河谷地区是暑热天气的多发区,特别是东北部的巫山、巫溪一带,中部的涪陵、丰都附近地区,以綦江为中心的西部偏南的沿江地区,更是三个高温集中区 [15] [16] [17] 。重庆全年42%的高温都出现于8月,持续极端高温伴随伏旱,使得电力负荷激增,森林火险等级处于高位,给电网安全运行带来不利影响。2017年夏季(6~8月),重庆全市平均气温较常年同期偏高0.7℃,35℃以上的高温日数较常年同期偏多12.4 d,共出现5次区域高温过程,其中7月16日~8月7日的区域高温过程,综合强度达到特重等级。图1所示为重庆某气象站2017年7月24日到7月30日的环境温度数据,以及同时期某双回线路供电断面的负荷电流数据。
由图1可知,负荷电流大小与温度高低大致呈正相关的关系,负荷电流随着温度升高而增大。在连续高温天气下,空调负荷大幅度增长,部分线路将处于重载甚至过载的工作状态。在多回线路并联运行时,当出现断面N−1事故时,剩余健全线路的电流会超过额定最大载流量,可能不满足N−1安全准则。因此,作为电网调度以及变电运维部门,在迎峰度夏用电高峰期间必须采取相关措施进行线路的载流能力评估,提前掌握架空线路的运行状态,保证供电断面和线路的运行安全。
Figure 1. Ambient temperature and load current from July 24 to July 30, 2017
图1. 2017年7月24日~7月30日环境温度与负荷电流
3. 架空线路的安全裕度
3.1. 导线载流的热平衡方程
架空线路的运行温度取决于线路运行过程中的载流值及其周围的环境气象条件(环境温度、日照强度、风速风向),导线温升响应满足暂态热平衡 [18] :
(1)
式中,m是单位长度导线质量,kg/m;c是导线的等效比热容,J/kg∙℃;式中,Te是导线温度,℃;qs、qJ、qc、qr分别是单位长度导线的日照吸热功率、焦耳热功率、对流散热功率、辐射散热功率,W/m。
对上式中qs、qJ、qc和qr的计算公式,简要说明如下:
1) 日照吸热功率qs
每单位长度的导体获得的日照发热功率qs,正比于导线的吸热系数α,导体外径D0 (m),导线所处地区海拔高度上的太阳辐射功率密度Qs (W/m2)以及太阳光入射方向与导线走向之间夹角θs (rad)。可由下式计算得到:
(2)
2) 焦耳热功率qJ
焦耳热由流过导体的电流决定:
(3)
其中I为载流值,A;R (Tc)为导线温度为Tc时的单位长度交流电阻值,Ω/m,可通过线性插值获得。
3) 对流散热功率qc
对流散热功率可以表示为努塞尔数Nu的函数,如式(4)所示:
(4)
其中λf是空气导热系数,Tc为导线表面温度,℃;Ta为导线环境温度,℃。
对流散热是导体散热的重要方式,导线的对流散热分为两种类型:自然对流散热和强迫对流散热,分别对应于无风和有风的情况。强迫对流散热功率取决于风速及风向角,在风速较高的情况下,强迫对流散热占主导地位,可忽略自然对流散热;在低风速下,起主要作用的则是自然对流散热。
强迫对流散热功率可按式(5)计算:
(5)
其中Nuδw为风向角为δw时的努塞尔数。
自然对流散热功率可按式(6)计算:
(6)
式中,其中Nunat为风速为零时的努塞尔数。
在计算导线的对流散热时,取qcn和qcf的最大值作为其对流散热功率,即:
(7)
4) 辐射散热功率qr
架空导线的辐射散热主要与导线和周围环境的温度差有关,根据斯蒂芬–玻尔兹曼定律,导线辐射散热功率为:
(8)
式中,ε是导线材料的辐射系数,σB是斯蒂芬–玻尔兹曼常数,其值为5.6697 × 10−8。
当式(1)左端等于零时,则表示架空线路的载流值确定,且周围气象条件处于稳定状态,其发热量和散热量处于平衡状态。即
(9)
在极端气象环境参数条件下,架空线路的运行温度则达到允许的最高温度,按照中国国家标准,钢芯铝绞线的最高允许温度为70℃。
3.2. 动态安全裕度表征方法
线路实际运行中,线路电流以及周围的环境温度、日照强度以及风速风向是不断变化的,线路的热平衡状态将会受到破坏,但由于热惯性的存在,导体的温度变化总是滞后于线路电流的变化。因此,本文引入架空线路的实时载流和实时温度两个判断指标,共同表征架空线路的动态安全裕度 [19] 。
1) 线路载流安全裕度
线路的动态载流量是根据实时气象环境参数以及线路运行最高允许温度所确定的,相比于架空线路的静态载流量,可以极大程度地挖掘架空线路的输电潜力,实现线路的有效增容,解决夏季连续高温天气下架空线路载流能力不足的问题。为确保线路的运行安全,并感知架空线路的载流能力,引入载流安全裕度概念,具体如图2所示。
线路载流安全裕度指当前线路的实时载流量到线路动态载流量之间的变化区间,可用如下公式表示:
(10)
式中,Imargin表示线路的载流安全裕度,Idy表示动态载流能力,It表示实时电流。
2) 导线温度安全裕度
线路的载流能力受电流、运行温度等因素限制,虽然运行时导体的温度与电流之间存在密切关联,但电流在表征线路载流能力上并不等同于运行温度的限制,制约线路载流能力的本质应是运行温度,温度升高时,线长增加,弧垂变大,安全距离减小,可能不满足对地或其他跨越物的电气距离;当线路温度近一步上升时,导线会因强度降低而断线,仅以电流表征输电线路的动态安全裕度会带来保守或者冒
Figure 2. Diagram of current carrying safety margin
图2. 载流安全裕度示意图
进的结果。因此,在挖掘线路送电能力需求以及准确把握线路动态安全裕度的驱使下,还应充分考虑制约线路载流能力的物理本质,从运行温度的角度对线路的安全运行状态进行评估。故引入线路的温度安全裕度概念,即
(11)
式中,Tmargin表示线路温度安全裕度,Tt表示线路实时运行温度,Tmax表示线路运行最高允许温度(中国标准一般采用70℃)。
4. 导线安全裕度评估与告警
4.1. 导线载流裕度评估
在连续高温天气下,空调负荷大幅度增长,部分线路将处于重载甚至过载的工作状态,电网调度以及变电运维部门必须完成架空线路动态载流能力的评估,实时掌握架空线路的运行状态,在线辨识架空线路的动态载流量,分析供电断面的载流安全裕度,校核供电断面的N−1风险水平,进行超前预控,确保供电断面的安全稳定运行与电力平稳有序供应。因此,需要进行架空线路动态载流量的实时计算以及线路载流安全裕度的确定。
由于各个气象环境参数与导线动态载流量之间的非线性关系,对线路动态载流量的影响是不同的,故必须对各种气象影响因素均进行考虑。风速是对线路动态载流量影响最大的环境因素,线路动态载流量对风速变化的灵敏度最高,同时环境温度与日照强度对线路动态载流量也有不同程度的影响。架空线路周围的气象环境参数是实时变化的,因此线路动态载流量的估计需要利用负荷电流以及短期气象环境参数。线路动态载流能力评估的框架具体如下(图3)。
借助短时气象预报系统,等时间间隔采集架空线路周围的气象环境参数值以及负荷电流大小,利用式(1)的导体暂态热平衡方程实时在线计算导体的动态载流量,构成一个等时间间隔的载流量随机序列,以此确定架空线路的动态载流值,再依据式(10)确定架空线路的载流安全裕度,完成架空线路的载流能力评估。评估精度由时间间隔决定,时间隔越短,灵敏性越高;时间间隔越大,灵敏性越低。
4.2. 供电断面N−1事故时导线温升告警
供电断面N−1事故是电网运行方式校核必须考虑的情况。当两条或多条线路构成的供电断面中的一条线路因为故障或其他原因退出运行时,潮流会转移到其他健全线路上,相当于在原有电流的基础上叠加了一个阶跃电流。由于热惯性的存在,架空线路的温度变化滞后于电流的变化,因此在输电线路负荷电流越限一段时间之后,导线温度才会上升至最高允许运行温度,到达输电线路安全运行的边界。电网调度工作人员可利用这段时间对电网结构进行优化,对将要达到最高允许运行温度的架空线路进行负荷转供,确保供电断面的安全稳定运行和电力平稳有序供应。因导线的热时间常数为min级,故可以假设
Figure 3. Framework of dynamic current carrying calculation
图3. 线路动态载流能力评估框架
该段时间的气象环境参数以及负荷电流不改变,利用式(1)计算得到线路从当前运行温度Tc0,即电流越限时刻t1,上升至线路最大允许温度所需时间tcs:
(12)
对系数A、B、K详细说明如下:
(13)
(14)
(15)
其中:Tmax为钢芯铝绞线的最高允许运行温度,为70℃;β为电阻随温度的变化率。
除获取完成负荷转供的时间裕度以外,还应实时辨识线路的温度安全裕度,进行线路运行温度的安全预警。可以以越限时刻为初始时刻,等时间间隔采样架空线路的运行电流以及实时气象环境参数,利用式(10)计算得到线路从越限时刻的运行温度Tt到最大允许运行温度Tmax的安全裕度Tmargin。利用实时安全裕度Tmargin完成运行温度的安全预警:
1) 当Tmargin逐渐减小时,且小至5℃时,则发出告警信号,告知作业人员线路的安全运行裕度不足,需要采取相关措施保证配电网的稳定与安全;当Tmar逐渐增大时,且大于5℃时,则解除告警信号,告知作业人员输电线路的安全运行裕度充足。
2) 当Tmargin ≤ 0℃时,调度运行人员必须完成负荷切除或者负荷转供,确保整个配电网的安全稳定。
电网调度工作人员可借助上述架空线路的运行温度预警方法准确地感知架空线路的温度安全裕度,以及利用热惯性提供的时间裕度对将要达到安全运行边界的输电线路完成负荷转供或者负荷切除,实现连续高温下配电网的协调控制,确保整个配电网的安全性与稳定性。
5. 算例分析
以图4所示的电气接线图为例,对连续高温天气下架空线路的动态载流能力进行评估以及运行温度的安全预警。选取重庆地区某110 kV供电断面作为算例,该断面通过双回线路向本地负荷以及下一级线路供电,其中架空线路的型号是2 × LGJ-300/25,根据国家标准《GB/T1179-2017圆线同心绞架空导线》 [20] ,该线路的具体参数如表1所示。以图1中的2017年7月27日为例,其温度和风速数据如图5所示。
Table 1. Parameters of LGJ-300/25
表1. LGJ-300/25参数
Figure 4. Diagram of power flow section
图4. 供电断面示意图
Figure 5. Weather parameters on July 27, 2017
图5. 2017年7月27日气象参数
连续高温天气下,当系统通过双回线路向负荷供电时,由气象环境参数实时确定的线路动态载流量如图6所示,可知连续高温天气下的日负荷与温度的变化呈正相关,在下午13:00达到了峰值836 A,在架空线路采用静态热定值运行方式下,该时刻留有178 A的载流安全裕度,而在动态热定值运行方式下则留有451 A的载流安全裕度,提高了近153%。由实时气象环境参数确定的动态载流量,相比于极端气象环境参数确定的静态载流值,可充分地提高输电线路的送电能力,即使在气象环境最恶劣的情况下,也留有438 A的载流安全裕度,相比于静态热定值运行方式的210 A,也提高了108%。因此,动态载流量可有效地提高输电线路的载流安全裕度,深度挖掘输电线路的输电潜力,确保该供电断面在迎峰度夏期间电力的平稳供应。
当送出线L1退出运行时,供电断面仅通过送出线L2向本地负荷以及下一级线路供电,L2将会承受因潮流转移引起的接近于2倍的负荷电流,供电断面通过单回线路供电时的动态载流量如图6所示。可知当供电断面发生N−1潮流转移之后,8:00~9:00之间的某个时刻架空线路的载流值将会超过静态最大载流值,但是并未越过架空线路的动态载流量,仍留有一定的载流安全裕度,在动态热定值下仍可以安全可靠运行。但是12:00~18:00当供电断面发生N−1潮流转移之后,该供电断面不再满足N−1准则,电网调度人员需要在该时间段内采取负荷转移或者切负荷等措施,保证架空线路留有足够的安全裕度,保证供电断面的可靠运行。例如,在13:00时,电网调度人员应该切除或者转移至少832.59 A负荷电流,才可以保证架空线路运行在载流安全边界之内。
由于热惯性的存在,架空线路的温度变化滞后于电流的变化,而运行温度又是制约架空线路安全运行问的物理本质,故还应实时辨识架空线路的运行温度。架空线路的实时温度如图7所示。
由图7可知,当系统通过双回线路向负荷供电时,输电线路具有充足的温度安全裕度,不存在越限时刻。当系统仅通过单回线路向负荷供电时,针对0:00到11:00之间供电断面发生N−1潮流转移之后,单回线路在动态热定值运行方式下能确保电力的平稳供应;但部分时段供电断面发生N−1潮流转移后,线路不再满足N−1安全准则。例如,当线路L1在12:00发生跳闸时,L2将会承受因潮流转移引起的负荷电流1642.68 A,发生了电流越限,不再满足N−1准则,故需要确定当前时刻上升至线路最大允许温度时刻的时间长度tcs,为电网调度人员完成负荷转供或者负荷切除提供参考。
利用式(12)计算如下可知tcs = 4.52 min,采用IEEE提供的导线温度对阶跃电流的响应曲线对该时间裕度的准确性进行验证 [21] ,计算的温升响应曲线如图8所示。
Figure 6. Dynamic capacity of overload lines
图6. 架空线路的动态载流能力
Figure 7. Real-time temperature of the overhead line
图7. 架空线路的实时温度
Figure 8. Response curve of transmission line temperature to step current
图8. 导线温度对阶跃电流的响应曲线
由图8可知,利用IEEE提供的导线温度对阶跃电流的响应曲线计算所得的时间裕度为4.75 min。由于式(12)推导过程中假设λf与v为常数,但导线在实际运行中λf与v是随着导线温度提高而增大的,故式(12)计算的时间裕度相比于IEEE标准计算的时间裕度会偏小,但误差范围是可接受的,故本文提出的时间裕度简化估算方法是合理的。以12:00越限时刻为初始时刻,开始进行线路气象环境参数的等时间间隔采样,进而计算输电线路的温度安全裕度,实现输电线路运行温度的安全告警,确保配电网的稳定安全运行。
6. 结论
针对连续高温天气下高负荷带来的线路运行安全问题,本文提出以动态载流量和运行温度共同表征架空线路的安全裕度,结合运行气象环境在线辨识架空线路的实时运行温度和动态载流能力,当载流裕度不足或导线温升越限时,向电网调度人员发出安全告警。通过研究得出以下结论:
1) 在连续高温天气下高负荷时,部分线路将处于重载甚至过载的工作状态,在多回线路并联运行时,当出现断面N−1事故时,剩余健全线路的电流会超过额定最大载流量,不满足N−1安全要求,因此必须采取相关措施保证供电断面和线路的运行安全。
2) 以动态载流量和运行温度共同表征的线路动态安全裕度,相比于静态安全裕度,可实时掌握输电线路的运行状态,进行超前预控,实现临近安全边界线路的负荷转移或者负荷切除,确保配电网的安全稳定运行与电力平稳有序供应。
3) 本文所提方法以架空线路的安全运行为前提,充分地利用了线路动态载流能力,可为连续高温天气下的电网得运行调度提供指导。
附录
公式(12)推导:
(1)
其中,日照吸热功率、焦耳热功率、辐射散热功率分别如下:
(2)
(3)
(4)
强迫对流散热功率取决于风速及风向角,在风速较高的情况下,强迫对流散热占主导地位,可忽略自然对流散热,则
(5)
其中:
(6)
(7)
(8)
(9)
上式之中,V为垂直于导线的风速,单位为m/s。
将(2)、(3)、(4)、(5)代入(1)可得:
(10)
(10)式中,
,β为电阻随温度的变化率。令
,可知C为常数,且在动态变化中环境温度不发生改变。可得:
(11)
令
为Tc的函数。可知当温度变化时,K、Nuδ以及λf均随温度的变化而变化。如(12)、(13)所示:
(12)
(13)
可知λf与v的变化很小,故可以近似认为λf与v为常数,则努塞尔数也可近似Nuδ为常数。
令
为常数,
,进一步简化方程(11)可得:
(14)
解微分方程(10)的初始条件为越限时刻时t1的运行温度Tc0,当t = ∞时,根据静态导线载流平衡方程可以求得越限时刻的电流对应的导线稳态温度Tmax,并且按此稳态温度求取K、Nuδ以及λf的值进而得到B的值。导体的温度随时间变化轨迹的近似表达式如公式(15):
(15)
导线通过的电流变大时,导线温度是慢慢上升,并最终达到静态平衡。由公式(15)求解可得导线电流增加时,导线温度上升至允许温度Tmax所需的时间为:
(16)