宋 剛,張 彤
(浙江省電力設計院,浙江 杭州 310012)
溪洛渡左岸—浙江金華±800 kV特高壓直流輸電線路工程是國內(nèi)又一條高等級、大容量、遠距離特高壓直流線路,路徑總長度為1670.8 km,途經(jīng)四川、貴州、湖南、江西、浙江5省,沿線經(jīng)過地形中,山地占45.43%、高山占19.78%,地勢起伏大,海拔高,全線中、重冰區(qū)長度占到70%以上,屬于覆冰極為嚴重的輸電線路,重冰區(qū)設計是該工程設計的重點問題之一。
在經(jīng)歷了2008年輸電線路冰災后,設計人員對事故倒塔段進行計算分析,發(fā)現(xiàn)倒塔主要發(fā)生在大檔距、大高差、覆冰不均勻的情況下,覆冰造成桿塔縱向張力差大于設計耐受值[1-7]。因此桿塔設計時覆冰不平衡張力取值關系線路能否安全穩(wěn)定運行,及能否經(jīng)受住再一次冰災考驗。目前重冰區(qū)設計主要參考DL/T5440—2009《重覆冰架空輸電線路設計技術規(guī)程》(以下簡稱規(guī)程),該規(guī)程中對不同冰厚下的不平衡張力取值做了規(guī)定。為驗證該規(guī)定的不平衡張力取值能否承受百年一遇的冰災,本文以該工程浙江段實際耐張段和假設斷面為例,計算嚴重情況下的導、地線覆冰不平衡張力,探討直線塔的覆冰不平衡張力取值,給出嚴重情況下的特高壓直流線路重冰區(qū)直線塔不平衡張力建議值。
規(guī)程規(guī)定產(chǎn)出不平衡張力的不均勻冰荷載情況按未斷線、溫度-5℃、有不均勻冰、風速10 m/s的工況計算[8]。不平衡張力覆冰率計算條件和不平衡張力取值(以不平衡張力占最大使用張力的百分比表示)分別見表1和表2。
上述縱向不平衡張力的取值是參考等連續(xù)檔、連續(xù)等高差上山模型計算結果和實際觀測、運行經(jīng)驗得到的。文獻[9]提出了“等連續(xù)檔,僅中間一基桿塔前后有高差,其余塔前后檔無高差”的計算模型,文獻[10]利用該模型計算20~50 mm冰區(qū)不平衡張力,得出結論為:30 mm及以下重覆冰區(qū),不平衡張力百分數(shù)按照現(xiàn)行重覆冰區(qū)規(guī)程規(guī)定取值;40、50 mm重覆冰區(qū),應將不平衡張力百分數(shù)分別提高至35%和41%。
表1 直線塔不平衡張力覆冰率計算條件Tab.1 Icing rate for suspension tower unbalanced tension calculation
表2 直線塔不平衡張力取值Tab.2 Settings of suspension tower unbalanced tension
不平衡張力與檔距、高差、覆冰厚度、不均勻冰情況等眾多因素相關[10-13],以往理論計算只分析了幾種假設斷面,沒有對不均勻檔距、復雜高差、驗算覆冰等情況進行綜合分析,得出的計算值適用于平地、丘陵地區(qū)。對于山地復雜地形,特別是大高差、大檔距地形,未有文獻進行過實際分析。本文將實際耐張段和在嚴重斷面、不同覆冰厚度情況下的不平衡張力計算結果與已有文獻和規(guī)程進行比較,得出適用于本工程的實際不平衡張力系數(shù),指導特高壓直流重冰區(qū)線路設計。同時針對不平衡張力經(jīng)典算法中存在的收斂速度慢、收斂于非正確解等問題提出解決辦法。
本文采用經(jīng)典算法進行不平衡張力計算[14-17]。該算法認為出現(xiàn)最大不平衡張力的情況一般發(fā)生在檔距、高差相差懸殊,線路一側覆冰而另一側少冰或無冰的交界桿塔處。
假設一個耐張段內(nèi)有n個連續(xù)檔、n-1基直線塔,架線氣溫為tm,導、地線初伸長尚未放出時架線應力為σm,各直線桿塔上懸垂絕緣子串均處于中垂位置。當出現(xiàn)不均勻冰時,各檔應力不一,懸垂串發(fā)生偏移,檔距發(fā)生變化,如圖1所示。
圖1 連續(xù)檔不平衡張力示意圖Fig.1 Schematic diagram of unbalanced tension of continuous spans
第i檔檔距增量Δli與檔內(nèi)應力σi間的關系為:
其中,li、βi分別為耐張段內(nèi)懸垂串處于中垂位置時第 i檔的檔距(m)和高差角(°);α、E 分別為導線的溫度線膨脹系數(shù)(℃-1)和彈性系數(shù)(N/mm2);tm、σm、Δte、γm分別為導線架線時氣溫(℃)、相應氣溫下耐張段內(nèi)的架線水平應力(N/mm2)、架線時考慮初伸長降低的等效溫度(取正值,單位為℃)、架線時導線和地線的自重比載(N/(mm2·m));t、σi、γi、Δli分別為計算不平衡張力時的氣溫(℃)、第i檔的水平應力(N/mm2)、比載(N/(mm2·m))和檔距增量(縮短時為負值,單位為m)。
由于待求情況下各檔水平應力不同,導致相鄰檔間懸垂串出現(xiàn)不平衡張力差,它使懸垂串出現(xiàn)偏斜,如圖2所示。
圖2 懸垂絕緣子串受力偏斜圖Fig.2 Pendantinsulator tilted by unbalanced tension
第i基直線塔上懸垂串末端導線懸掛點的順線路水平偏移δi與兩側導線張力差的關系為:
通過式(2)、(3)解出 σi+1的顯函數(shù)為:
其中,σi為第 i檔的水平應力(N/mm2);δi為第 i基直線塔懸垂串導線懸掛點順線路水平偏距(m),偏向大編號為正,反之為負;λi、Gi分別為第 i基直線塔的懸垂串串長(m)及荷載(N);Wi為第 i基直線塔導線垂直荷載(N/mm2);hi、hi+1分別為懸垂串處于中垂位置時,第i基直線塔對第i-1基直線塔和第i+1基直線塔對第i基直線塔上導線懸掛點間的高差(m),較大編號的桿塔比較小編號的桿塔懸掛點高時為正值,反之為負值;βi為懸垂串處于中垂位置時,第i檔導線懸掛點間的高差角(°);li為懸垂串處于中垂位置時,第 i檔的檔距(m);A 為導線截面積(mm2)。
對于整個耐張段內(nèi),各檔檔距增量之和應為0,即第n基耐張塔上導線懸掛點的偏距應為0,即:
由式(1)可以列出n個方程,由式(4)可以列出n-1個方程,由式(5)可以列出1個方程,共2n個方程,有Δli、σi共2n個未知數(shù),方程是可以求解的。較為方便的是采用迭代求解法,通過迭代直至獲得最終解,迭代步驟如下:
a.令i=1,δ0=0;
b.將σi代入式(1)得 Δli;
c.令δi=δi-1+Δli;
d.將σi、δi代入式(2)得 δi+1;
e.將δi+1代入式(4)得到σi+1;
f.將σi+1代入式(1)得 Δli+1;
g.δi+1=δi+Δli+1,i=i+1;
h.如果i=n,則轉步驟i,否則轉步驟d;
i.如果δi=0,結束迭代;否則令 σ1=σ1+Δσ,轉步驟a。
經(jīng)典算法存在的主要問題是初始值σ1、迭代步長Δσ、收斂精度的設置和收斂結果判斷。初始值設置過小,需要迭代較多次數(shù),收斂時間長,且可能導致第一個收斂結果為非正確解,此時不平衡張力百分比將超過100%,這是由于經(jīng)典算法可能存在雙解現(xiàn)象;初始值設置過大,則無法得到解。迭代步長需和迭代精度配合,較大的迭代步長可能無法達到所需的迭代精度,轉而跳過計算解。
針對以上問題,在編程過程中先取較小初值σ0、較大的迭代步長(如 Δσ=100 N/mm2)和較大的判斷精度(如δn=0.1 m),初次收斂后取一稍小于上次計算收斂值σ1的值(如0.99σ1),減小迭代步長(如Δσ=10 N/mm2)和判斷精度(如δn=0.01m),如此再迭代一次即滿足工程設計精度要求。針對第一次收斂結果為非正確解的情況,判斷不平衡張力百分比是否大于100%,如果大于則將初值適當增大(如增加為σ1+1),繼續(xù)進行迭代計算直至得到最終解。通過該改進措施能保證一次得到正確解,且無需人工干預,減少收斂時間,提高收斂效率。
溪洛渡左岸—浙江金華±800 kV特高壓直流輸電線路(簡稱溪浙線)20 mm冰區(qū)的導線采用6分裂 JL/G2A-900/75鋼芯鋁絞線,地線采用JLB20A-150鋁包鋼絞線,具體參數(shù)見表3。
表3 導線和地線參數(shù)Tab.3 Parameters of conductors and ground wires
表4列出了某一實際耐張段的前側塔、后側塔的檔距及高差(后側塔相對前側塔,后同)情況,表中編號4145和4149的塔為耐張塔,其余為直線塔。
表4 實際耐張段地形情況Tab.4 Terrain of an actual strained segment
根據(jù)第2節(jié)不平衡張力計算方法,利用MATLAB編制程序,得到覆冰20 mm情況下的不平衡張力計算結果如表5、6所示。
由表5可見,4148號直線塔導線和地線的最大不平衡張力百分比分別達到29.96%、51.45%,超過設計規(guī)程規(guī)定的20 mm冰區(qū)直線塔導線不平衡張力百分比取值(25%)和地線不平衡張力百分比取值(46%)。
如果按照25 mm覆冰驗算(直線塔一側25 mm覆冰,另一側根據(jù)不均勻覆冰率取值),導線最大不平衡張力百分比將達到40.81%,地線最大不平衡張力百分比達到71.88%。
表5 直線塔20 mm覆冰不平衡張力計算結果Tab.5 Calculated unbalanced tension of suspension tower covered with 20 mm ice
對同一耐張段,假設設計覆冰分別為30 mm和40 mm,導線采用 6 分裂 JLHA1/G1A-800/55 鋼芯鋁合金絞線,地線采用JLB20A-240鋁包鋼絞線(導地線參數(shù)見表3),計算導、地線最大不平衡張力如表6所示。
由表6可見,根據(jù)實際地形計算的直線塔導線不平衡張力百分比較規(guī)范取值高5%~16%,地線高5%~9%。如果再考慮線路附近可能存在水汽增加覆冰的影響,應該在設計規(guī)范基礎上更大幅度提高不平衡張力系數(shù),以提高線路抗冰能力。
以上實例是選取了溪浙線特高壓直流工程浙江境內(nèi)一環(huán)境較惡劣的耐張段進行計算。另外,文獻[11]中已經(jīng)說明,不平衡張力隨著檔數(shù)增加而增加,文獻[10,12-13]論證了不平衡張力隨著檔距、高差增大而增大。特高壓直流線路一般都達到1500 km以上,經(jīng)過地形復雜,不排除在重覆冰地區(qū)存在一些連續(xù)檔數(shù)量多、大檔距、大高差同時存在、靠近水汽豐富地區(qū)等更嚴峻情況出現(xiàn),不平衡張力值無法完全參考規(guī)程取值。
另外通過實例計算,可以看出兩側承受大檔距的“一肩挑”式直線塔受到的不平衡張力最大。本文計算了一種假設的五檔斷面500-500-500-X-Y在各覆冰條件下的不平衡張力情況,其中X和Y可取700~1100 m之間的值,X和Y的高差h可取100 m、150 m和200 m,其余各檔無高差。計算結果如表7所示。
從表7可以看出,重冰區(qū)直線塔兩側存在大檔距、大高差時,不平衡張力容易超過規(guī)程值,地線不平衡張力超過規(guī)程值更明顯。近幾年覆冰事故中,大部分為地線支架倒塌或者地線支架先倒導致整塔折斷,除了地線覆冰較導線厚外,地線懸垂串較短導致不平衡張力也是引發(fā)事故的重要原因。另外可以看出覆冰越厚,越容易出現(xiàn)直線塔不平衡張力超規(guī)程值的情況。
表7 重冰區(qū)直線塔最大不平衡張力百分數(shù)Tab.7 Maximum unbalanced tension percentage of suspension tower located in area with heavy ice accretion
因此規(guī)程取值適用于大部分耐張段,具有一定的經(jīng)濟合理性。當直線塔存在大檔距(大于900 m)、大高差(超過100 m)配合或者靠近水汽豐富地區(qū)容易局部結冰時,建議按表8對不平衡張力取值,其中各冰區(qū)導線、地線選型參考溪浙線重冰區(qū)導線選型結論。
表8 重冰區(qū)直線塔最大不平衡張力百分數(shù)建議值Tab.8 Suggestive value of maximum unbalanced tension of suspension tower located at heavy ice accretion area
按表8取值基本能涵蓋目前重冰區(qū)大檔距、大高差情況下的直線塔不平衡張力范圍。
a.重冰區(qū)規(guī)程中對直線塔不平衡張力取值適用于一般地形,對存在大檔距、大高差情況下規(guī)程取值偏小。特高壓直流線路重冰區(qū)需根據(jù)實際地形、覆冰情況校核計算。
b.特高壓直流輸電線路重要性高,同時投資也較大,特別是鐵塔投資占40%左右。以往桿塔計算結果表明,不平衡張力主要控制桿塔塔身側面斜材和塔身橫隔面的部分桿件,對鐵塔重量影響較大,建議設計1個或2個大塔型適當提高不平衡張力系數(shù),滿足工程中特殊地形直線塔需求,而無需提高所有塔型的不平衡張力系數(shù),以節(jié)約工程投資。
c.針對重覆冰地區(qū)可能出現(xiàn)的大檔距(大于900 m)、大高差(超過100 m)配合或者靠近水汽豐富地區(qū)容易局部結冰的情況,建議直線塔導線不平衡張力取值如下:20 mm重冰區(qū)取30%,30 mm冰區(qū)取35%,40 mm冰區(qū)取50%;地線不平衡張力建議取值:20 mm重冰區(qū)取50%,30 mm冰區(qū)取55%,40 mm冰區(qū)取65%。
d.根據(jù)實際情況,計算驗算覆冰工況下的直線塔不平衡張力。
e.特高壓直流線路重覆冰地區(qū)應盡量避讓大檔距、大高差、水汽豐富地區(qū)。