胡 云， 師 濤， 王文林， 王 斌， 覃日升

(1. 云南電網(wǎng)有限責(zé)任公司 玉溪供電局 玉溪 653100；2. 云南電力試驗(yàn)研究院(集團(tuán))有限公司 電力研究院 昆明 650000)

0 引言

架空線路是電力系統(tǒng)中重要的能源輸送通道，其所處自然環(huán)境復(fù)雜多樣，因而極易遭受雷害威脅。運(yùn)行經(jīng)驗(yàn)表明，線路因雷擊引起的線路跳閘事故占總跳閘事故的40%～70%[1]。雷擊引起的沖擊電暈會造成過電壓波形的衰減和畸變，從宏觀角度看，沖擊電暈增加導(dǎo)線周圍的空間電荷，進(jìn)而造成導(dǎo)線對地電容的增加。Wagner和Lloyd[2]通過實(shí)驗(yàn)對導(dǎo)線的庫伏(q-v)曲線進(jìn)行了實(shí)測，通過q-v曲線的斜率提取等值電容，在仿真中將傳輸線分成幾段π型等值電路，并在每個π型等值電路的連接點(diǎn)添加該等效電容來分析沖擊電暈對雷電沖擊下暫態(tài)過程的影響。

國內(nèi)外學(xué)者采用了不同方法研究了雷擊導(dǎo)線情況下沖擊電暈對輸電線路過電壓水平的影響，一致認(rèn)為沖擊電暈效應(yīng)會造成電壓波形的衰減和畸變，且沖擊電暈效應(yīng)引起的過電壓波前衰減和畸變遠(yuǎn)大于因?qū)Ь€自身電阻引起的相應(yīng)損耗，在仿真計算中可以忽略導(dǎo)線的電阻[3～5]。對于直擊雷而言，判斷電暈的起暈條件是以導(dǎo)線電壓u超過電暈臨界起暈電壓u0，且導(dǎo)線電壓未達(dá)到峰值。

相對于高壓輸電線路，配電線路的高度較低，絕緣等級不高，大多數(shù)的雷電事故是因雷擊在線路附近大地或建筑物等產(chǎn)生的感應(yīng)雷過電壓造成的[6～8]。之前人們在計算感應(yīng)雷過電壓時很少考慮電暈的影響，大多僅考慮了有損大地的影響。后來一些學(xué)者發(fā)現(xiàn)沖擊電暈會導(dǎo)致過電壓波形的畸變，它的影響與有損大地同樣重要。2000年和2010年，Nucci等人[9,10]基于Agrawal場線耦合模型[11]對考慮沖擊電暈的單導(dǎo)體架空線路過電壓進(jìn)行了計算，模型中將電暈等效成一個隨導(dǎo)線電壓變化的動態(tài)電容，分析了沖擊電暈對感應(yīng)雷過電壓的影響。在仿真中，Nucci等人沿用了直擊雷情況下對電暈的處理方法，以導(dǎo)線上的全電壓作為起暈依據(jù)，實(shí)現(xiàn)對雷電感應(yīng)過電壓的計算與分析。而由電暈機(jī)理可知，電暈產(chǎn)生的根本原因是電場，2014年，Thang[12]從電磁場的角度基于3維時域有限差分(FDTD)法開展了考慮沖擊電暈的架空線感應(yīng)雷過電壓計算方法的研究，但此種方法計算量龐大，對計算機(jī)的內(nèi)存和計算能力要求較高。2017年，王吉文[13]借助Agrawal場線耦合模型，采用FDTD方法分析了雷擊沖擊電暈對架空線路不同位置處感應(yīng)雷過電壓的影響。

為了綜合分析沖擊電暈對架空線路直擊雷和感應(yīng)雷過電壓的影響，本文基于沖擊電暈的等效電容模型，結(jié)合多導(dǎo)體傳輸線電報方程，建立了考慮沖擊電暈的雷電過電壓計算模型，并基于時域有限元法(TDFEM)對模型進(jìn)行了求解。此外，重點(diǎn)分析了沖擊電暈對感應(yīng)雷過電壓的影響機(jī)理，結(jié)果表明，對于感應(yīng)雷，導(dǎo)線周圍場強(qiáng)主要受散射電場決定，而不應(yīng)采用與直擊雷過電壓相同的判據(jù)，需采用散射電場作為判據(jù)并以此確定電暈套半徑。因此，沖擊電暈對感應(yīng)雷過電壓的影響與直擊雷過電壓并不相同，在電壓達(dá)到峰值后，線路仍將產(chǎn)生沖擊電暈，進(jìn)而造成感應(yīng)雷過電壓波形的波尾畸變。本文方法所得結(jié)果與電磁場方法的分析結(jié)果一致，同時由于采用電路模型，可有效降低仿真復(fù)雜度，提高計算效率。

1 建模方法

1.1 沖擊電暈的伏庫特性

導(dǎo)線上發(fā)生暈后的電氣參數(shù)的變化可以通過q-v曲線特性來表征[14]。q-v曲線特性反映了導(dǎo)線起暈后，導(dǎo)線上的電壓瞬時值u與導(dǎo)線上及其周圍電暈套內(nèi)的總電荷q的關(guān)系，用函數(shù)關(guān)系表示為

q=f(u)

(1)

圖1為導(dǎo)線發(fā)生沖擊電暈時的q-v特性曲線，曲線的斜率近似表征導(dǎo)線的電容，即C=dq/du。

圖1 導(dǎo)線的沖擊電暈q-v曲線

如圖所示，q-v曲線大體可以分為三部分，即：OA段：波頭部分(uC0；BC段：u

從能量守恒角度分析，沖擊電暈伴隨著能量消耗，過電壓波形會發(fā)生明顯的衰減。由于導(dǎo)線等效電容增大，導(dǎo)線波阻抗及波的傳播速度都將減小，進(jìn)而造成波形的畸變。電壓等級越高的線路，沖擊電暈對沿線傳播的雷電波的影響也越大。

1.2 考慮沖擊電暈輸電線路模型

考慮大地?fù)p耗的N+1線多導(dǎo)體傳輸線時域電報方程為

(2)

式中：u和i分別表示傳輸線上沿線各點(diǎn)的電壓和電流列向量；R0，G0，L0，C0分別為傳輸線的單位長電阻、電感、電導(dǎo)、電容矩陣；Zf(t)表示大地?fù)p耗和導(dǎo)體集膚效應(yīng)對應(yīng)的總的阻抗矩陣[15]，其維數(shù)為N×N。

考慮到?jīng)_擊電暈相當(dāng)于擴(kuò)大了導(dǎo)線的等效半徑，對電感參數(shù)的影響很小，僅是增大了導(dǎo)線的電容和電導(dǎo)，同時由于單位長度電納參數(shù)對計算結(jié)果影響甚微，所以可以用附加電容等效沖擊電暈的影響。

(1)考慮沖擊電暈的直擊雷過電壓計算模型

結(jié)合電暈等效電路模型，可得到考慮沖擊電暈多導(dǎo)體傳輸線電報方程

(3)

式中：電容矩陣C(u)為描述電暈效應(yīng)的單位長度等效電容，該電容可通過計算電暈等效半徑結(jié)合同軸導(dǎo)體電容計算公式進(jìn)行計算。關(guān)于電暈等效半徑rc的計算，可依據(jù)如下公式迭代計算

(4)

式中：U為導(dǎo)體電位；r0為導(dǎo)體半徑；h為導(dǎo)線高度；Ec為臨界起暈場強(qiáng)；對于正極性電暈取Ec+=0.4～0.5 MV/m，負(fù)極性電暈取Ec-=1.1～1.8 MV/m。可以看出，由于沖擊電暈的區(qū)域受線路周圍電場影響，因此線路電容參數(shù)為電壓的非線性函數(shù)，因此所建模型為非線性模型。

近幾年，TDFEM[16]廣泛應(yīng)用于多導(dǎo)體傳輸線的求解，相比于時域有限差分法(FDTD)[17]，TDFEM具有邊界條件方程易于處理和無條件穩(wěn)定的優(yōu)點(diǎn)，為此，本文選取TDFEM對所建立的模型進(jìn)行求解。

求解過程中對非線性電容的可通過如下公式進(jìn)行計算

(5)

(6)

(7)

式中：Cl-corona為電暈區(qū)域的等效電容；Ccorona-g為電暈套與大地之間的電容；C(u)為導(dǎo)線與大地間的等效非線性電容。

(2)考慮沖擊電暈感應(yīng)雷過電壓計算模型

Agrawal場線耦合模型被廣泛應(yīng)用于求解外場激勵下的多導(dǎo)體傳輸線暫態(tài)響應(yīng)，其與傳統(tǒng)電報方程之間的區(qū)別在于需在電壓方程右側(cè)附加入射水平電場，具體形式為一個非齊次偏微分方程，即

(8)

同樣地，采用TDFEM計算考慮沖擊電暈的感應(yīng)雷過電壓，電暈?zāi)Ｐ瓦x取與直擊雷情況的相同，但區(qū)別在于通過導(dǎo)線上的散射電壓而非全電壓求取沖擊電暈等效電容，其原因?qū)⒃谙挛慕o出。

2 算例驗(yàn)證

2.1 直擊雷過電壓算例分析

本部分選取文獻(xiàn)[18]中的實(shí)驗(yàn)為例，以驗(yàn)證本文所提模型及求解方法的準(zhǔn)確性。該實(shí)驗(yàn)由兩組導(dǎo)線構(gòu)成，一組為單導(dǎo)體，另一組為4分裂導(dǎo)線。其中單導(dǎo)體導(dǎo)線長度為1 410 m，離地高度為22.2 m，4分裂導(dǎo)線距離地面14 m，與單導(dǎo)體導(dǎo)線距離2 m，分裂導(dǎo)線間距0.4 m。在線路首端施加負(fù)極性沖擊電壓，波形如圖3所示，分別記錄不同位置處的電壓波形。利用本文方法，計算距離線路首端352.5 m、705 m和1057.5 m處的電壓波形，如圖2所示。通過對比仿真結(jié)果和實(shí)驗(yàn)結(jié)果，可以看出兩者吻合良好，驗(yàn)證了本文所提方法的準(zhǔn)確性和有效性。

圖2 仿真結(jié)果與實(shí)測結(jié)果對比圖

2.2 感應(yīng)雷過電壓算例分析

選取某10 kV配網(wǎng)線路為例，對考慮沖擊電暈的感應(yīng)雷過電壓進(jìn)行分析計算。架空線路長度為1 000 m，離地高度為10 m，三相導(dǎo)線水平排列，導(dǎo)線間距為3.7 m，導(dǎo)線半徑為9.14 mm。雷擊點(diǎn)在線路中點(diǎn)處，距離A相導(dǎo)線距離為R，如圖3所示。算例中選取的雷電流和雷電通道模型分別為Heilder雷電流模型和MTLE回?fù)裟Ｐ?，具體選取參數(shù)與文獻(xiàn)[19]一致。

圖3 水平排列導(dǎo)線感應(yīng)雷計算示意圖

選取大地電導(dǎo)率為σg=0.01 S/m，雷擊點(diǎn)距離A相導(dǎo)線為R=50 m，激勵選取幅值60 kA的負(fù)極性雷，計算仿真結(jié)果如圖4和圖5所示。

圖4 A相導(dǎo)線中點(diǎn)①處的電壓波形

圖5 A相導(dǎo)線端點(diǎn)②處的電壓波形

可以看出，感應(yīng)雷過電壓的數(shù)值比不考慮電暈時的數(shù)值要高4%～7%，這與直擊雷恰好相反。因此，對于感應(yīng)雷同樣需要考慮沖擊電暈的影響。

眾所周知，電暈放電存在能量消耗，會降低過電壓峰值，而沖擊電暈造成了感應(yīng)雷過電壓峰值的升高，這是因?yàn)闆_擊電暈會造成導(dǎo)線上散射電壓的峰值下降，從圖6中可以看到散射電壓的值為負(fù)值，根據(jù)全電壓的計算公式u(x)=us(x)+uinc(x)，全電壓的值就會因電暈而增大。

圖6 全電壓和散射電壓波形

此外，從本文的仿真結(jié)果可以發(fā)現(xiàn)，當(dāng)電壓過波峰之后考慮電暈與不考慮電暈的電壓波形曲線一直呈現(xiàn)分離的狀態(tài)，這與學(xué)者Nucci仿真結(jié)果存在很大的差異，這是由于學(xué)者Nucci選用全電壓u作為起暈判據(jù)，而本文選用導(dǎo)線周圍電場強(qiáng)度E作為起暈判據(jù)。

圖7 t=5 μs時刻A相導(dǎo)線觀測點(diǎn)①正下方電場分布圖

在實(shí)際計算中應(yīng)以散射電壓為依據(jù)確定電暈等效電容，這恰恰是造成沖擊電暈對感應(yīng)雷過電壓影響規(guī)律與直擊雷相反的原因。本文方法彌補(bǔ)了Nucci等人在考慮沖擊電暈感應(yīng)雷計算模型的不足，同時也避免了采用3維全場方法復(fù)雜度大的缺點(diǎn)。因此本文方法與現(xiàn)有方法相比，具有準(zhǔn)確高效的優(yōu)點(diǎn)。

3 結(jié)論

本文在電報方程的基礎(chǔ)上，結(jié)合沖擊電暈的等效模型，建立了考慮沖擊電暈的雷電過電壓計算模型?；谒⒌哪Ｐ?，對線路直擊雷過電壓和感應(yīng)雷過電壓進(jìn)行了計算分析。仿真結(jié)果表明：

(1)對于直擊雷過電壓，沖擊電暈會造成過電壓波形的衰減和畸變，在雷電過電壓計算中應(yīng)充分考慮其影響。

(2)對于感應(yīng)雷過電壓，沖擊電暈會造成過電壓波形的畸變。但其影響規(guī)律與直擊雷情況存在區(qū)別，考慮沖擊電暈感應(yīng)雷過電壓的數(shù)值會比不考慮電暈時的數(shù)值要高4%～7%。

(3)與直擊雷不同，在考慮沖擊電暈的感應(yīng)雷過電壓計算時，起暈的判別條件應(yīng)以導(dǎo)線周圍電場強(qiáng)度，而不是以全電壓作為起暈依據(jù)。

參考文獻(xiàn)：

[1] VITA V, EKONOMOU L, CHRISTODOULOU C A. The impact of distributed generation to the lightning protection of modern distribution lines[J]. Energy Systems, 2016, 7(2):357-364.

[2] WAGNER C F, LLOYD B L. Effects of corona on traveling waves[J]. IEEE Transactions on Power Apparatus & Systems, 1955, 74(3):858-872.

[3] 黃克捷. 沖擊電暈特性及其對線路雷電瞬態(tài)計算的影響[D]. 北京：北京交通大學(xué), 2017.

[4] 王榮珠. 雷擊高壓輸電線路沖擊電暈特性及其影響研究[D]. 南京：南京信息工程大學(xué), 2015.

[5] 陳則煌, 王榮珠, 周中山,等. 沖擊電暈對輸電線路中雷電波傳播的影響研究[J]. 電瓷避雷器, 2015(2):94-98.

[6] YU Z, ZHU T, WANG Z, et al. Calculation and experiment of induced lightning overvoltage on power distribution line[J]. Electric Power Systems Research, 2016, 139:52-59.

[7] RIZK M E M, MAHMOOD F, LEHTONEN M, et al. Computation of peak lightning-induced voltages due to the typical first and subsequent strokes considering high ground resistivity[J]. IEEE Transactions on Power Delivery, 2017, 32(4):1861-1871.

[8] PIANTINI A. Lightning-induced overvoltages on overhead medium-voltage lines[C]// IEEE International Conference on High Voltage Engineering and Application. IEEE, 2016:1-7.

[9] NUCCI C A, GUERRIERI S, BARROS M T C D, et al. Influence of corona on the voltages induced by nearby lightning on overhead distribution lines[J]. IEEE Transactions on Power Delivery, 2000, 15(4):1265-1273.

[10] DRAGAN G, FLOREA G, NUCCI C A, et al. On the influence of corona on lightning-induced overvoltages[J]. The Progressive Fish-Culturist, 2010, 33(4):187-192.

[11] AGRAWAL A K, PRICE H J, GURBAXANI S H. Transient response of multiconductor transmission lines excited by a nonuniform electromagnetic field[J]. IEEE Transactions on Electromagnetic Compatibility, 2007, 22(2):119-129.

[12] THANG T H, BABA Y, NAGAOKA N, et al. FDTD simulations of corona effect on lightning-induced voltages[J]. IEEE Transactions on Electromagnetic Compatibility, 2014, 56(1):168-176.

[13] 王吉文. 雷擊電暈對輸電線感應(yīng)電壓的影響[J]. 電瓷避雷器, 2017(3):133-137.

[14] 董博文. 考慮頻變參數(shù)的架空線分?jǐn)?shù)階建模及其雷電過電壓計算方法的研究[D]. 保定：華北電力大學(xué), 2016.

[15] LIU X, CUI X, QI L. Time-domain finite-element method for the transient response of multiconductor transmission lines excited by an electromagnetic field[J]. IEEE Transactions on Electromagnetic Compatibility, 2011, 53(2):462-474.

[16] LI J, HUANG Y. Time domain finite element methods for metamaterials[M]. Time-Domain Finite Element Methods for Maxwell's Equations in Metamaterials. Springer Berlin Heidelberg, 2013:53-125.

[17] BRIGNONE M, DELFINO F, PROCOPIO R, et al. Evaluation of power system lightning performance, part I: model and numerical solution using the PSCAD-EMTDC platform[J]. IEEE Transactions on Electromagnetic Compatibility, 2016, 59(1):137-145.

[18] INOUE, A. Study on propagation characteristics of high-voltage traveling waves with corona discharge[J]. Japanese CRIEPI, 1983, 114-123

[19] LIU X, YANG J, WANG L, et al. A new method for fast evaluation of the cooray-rubinstein formula in time domain[J]. IEEE Transactions on Electromagnetic Compatibility, 2017, 59(4):1188-1195.