王振宇，蘇詠梅，易善明

（1.鄭州職業(yè)技術(shù)學(xué)院，鄭州450121；2.國網(wǎng)河南省電力公司，鄭州450018）

0 引言

對雷電回擊電磁場以及其對周圍建筑、架空線路感應(yīng)的過電壓研究已經(jīng)成為當(dāng)今雷電方向的重要研究課題。在電網(wǎng)方面，造成輸電線路故障的主要原因之一是由于雷擊感應(yīng)電壓引起的[1-2]。而雷擊輸電線的感應(yīng)過電壓主要為入射電壓和散射電壓構(gòu)成，且在架空線不同位置處的感應(yīng)電壓值不相同，因此，研究架空線不同位置處入射電壓和散射電壓變化特征，對線路防雷具有重要意義[3]。

國內(nèi)外一些學(xué)者對閃電引起的地磁場傳播機制以及耦合過電壓進行了研究。Bruce等[4]首次提出雷電回擊模型這一概念。Norton[5]在早期提出了“斜波法”，來研究閃電通道周圍電磁場變化特征。Cooray等[6]在研究地面土壤為有限電導(dǎo)率情況下，對閃電回擊電磁場的影響時，提出了一種近似解析的算法，該方法能夠計算出距離通道幾百米范圍內(nèi)的電磁場變化情況。在后期Cooray[7]通過根據(jù)大地的阻抗，研究了回擊電磁場變化特征。Rubinstein[8]在Cooray研究基礎(chǔ)上，提出了一種適用于更遠范圍的回擊電磁場計算方法。李福壽[9]對配電線路的雷擊過電壓進行了分析。文習(xí)山等[10]利用時域有限差分方法，對雷擊架空線路耦合過電壓峰值變化趨勢進行了分析。還有其他學(xué)者采用不同方法，對雷擊過電壓討論，分析[11-13]。

雷擊輸電線的感應(yīng)過電壓主要為入射電壓和散射電壓構(gòu)成，上述學(xué)者，在對架空線路雷擊過電壓研究工作中，均沒有研究入射電壓以及散射電壓變化情況。因此，筆者通過對感應(yīng)過電壓進行拆分分析，分析計算了架空輸電線中點和端點處感應(yīng)入射、散射電壓變化特征。

1 雷電電磁場數(shù)值計算方法

1.1 地閃回擊模型

閃電通道底部基電流主要有兩種：雙指數(shù)電流模型和Heidler模型。Heidler模型[14-15]是Heidler于1985年提出的一種較新的雷電電流的數(shù)學(xué)模型：

式中：τ1和τ2分別為波頭和波尾的時間常數(shù)；n為電流波形陡度因子，通常取n=2；η是峰值電流修正因子，其表達式為

通道底部電流采用的是n階（n=2）雙Heidler函數(shù)，其表達式為

式中：I01和I02為回擊電流的幅值；τ11和τ21為電流波形的上升沿；τ12和τ22為電流波形的下降沿。

1.2 算法介紹

筆者采用二維FDTD算法，計算通道周圍電磁場值，見圖1。假設(shè)閃電通道垂直于地面，由于閃電通道為軸對稱類型，因此采用閃電通道模型建立在二維柱坐標系下，此時的激勵源為包括H?、Er和Ez的TM波。

圖1 FDTD算法計算雷電通道周圍電磁場模型Fig.1 FDTD algorithm calculates the electromagnetic field model around the lightning channel

根據(jù)中心差分方法和電磁場各分量節(jié)點位置，TM波的Maxwell方程離散后整理為

1.3 耦合模型

在對雷擊線路，耦合過電壓的計算模型研究工作中，常常采用Agrawal耦合計算模型，該耦合模型能夠?qū)崿F(xiàn)將電場轉(zhuǎn)化成雷擊電壓值，該模型在考慮傳輸線和雷電通道周圍電磁場耦合時，只考慮了空間電場對耦合過電壓的影響，圖2為輸電線路Agrawal耦合計算模型。

Agrawal模型的方程如下：

式中：L′和C′為傳輸線單位長度電感和單位長度電容；I(x)為傳輸線電流；Vs(x)為散射電壓。散流電壓Vs(x)與總電壓V(x)的關(guān)系是：

式中，Ve(x)是激勵電壓。Ve(x)可表示為

求解上述模型方程的端部條件為

圖2 Agrawal模型等效電場激勵傳輸線Fig.2 Equivalent electric field excited transmission line of Agrawal model

2 模擬結(jié)果分析

為了得到雷電電磁脈沖在架空輸電線上耦合的感應(yīng)過電壓，必須要知道雷電回擊通道與架空輸電線之間的相對位置關(guān)系以及架空線長度、架高等幾何尺寸和雷電回擊放電參數(shù)，建立場-線一體化模型。見圖3。筆者計算架空輸電線的感應(yīng)過電壓時選取了一組參數(shù)對一體化模型中的所有參量進行賦值，即設(shè)架空輸電線長度L為1 200 m，架高h為10 m，輸電線半徑a為0.005 m，兩終端電阻R1=R2=498 Ω，閃電通道與線路中心位置處的垂直距離d為100 m，土壤電導(dǎo)率σ為0.01 S/m，介質(zhì)介電常數(shù)ε為10 F/m，雷電回擊通道高度H為7 500 m，回擊電流速度v為1.5×108m/s，閃電回擊模型選取的是MTLL模型，通道底部雷電采用雙Heidler基電流函數(shù)。

圖3 場-線一體化模型Fig.3 Integrated model of field-line

2.1 雷電電磁場的計算

因為我們計算感應(yīng)電壓采用的是Agrawal耦合模型，所以計算時只考慮回擊通道周圍電場，包括水平電場和垂直電場。計算結(jié)果如圖4所示。從圖4中可以看出，雷電回擊電流在架空線中點處產(chǎn)生的垂直電場最大。架空線中點的水平電場垂直于架空線，而水平電場的垂直分量對感應(yīng)電壓的計算沒有作用，我們計算感應(yīng)電壓時所用的是水平電場的沿架空輸電線的切向分量，程序中我們采用了線性插值的方法算出了水平電場的切向分量，所以在輸出結(jié)果中架空輸電線中點的水平電場為零。在線路端點位置處雷擊感應(yīng)水平電場主要呈雙極性；水平電場主要用來計算散射電壓，垂直電場主要用來計算入射電壓。

圖4 雷電回擊電流在架空線上產(chǎn)生的電場Fig.4 Electric field produced by lightning counterattack current on overhead lines

2.2 感應(yīng)過電壓結(jié)果分析

計算出水平電場和垂直電場之后，將其作為激勵源帶入Agrawal耦合模型中計算感應(yīng)過電壓。雷電在架空線上耦合的感應(yīng)過電壓是由散射電壓和入射電壓兩部分組成的，其中散射電壓主要由水平電場激勵產(chǎn)生，根據(jù)經(jīng)過傅里葉逆變換得到時域范圍Agrawal耦合模型得出的。

式中，Eer(x,h,t)是沿架空輸電線方向的水平電場值，也就是Agrawal耦合模型中的激勵源項，在程序中是根據(jù)線性插值的方法得到的。而入射電壓由垂直電場激勵產(chǎn)生，架空線上的任一點的垂直電場分量也可以根據(jù)線性插值的方法得到。圖5給出了架空輸電線中點和端點的感應(yīng)總電壓及其分量的波形圖。

圖5 架空線雷電感應(yīng)總電壓和其兩個分量波形圖Fig.5 The total overhead line lightning induced voltage and its two component waveforms

從圖5中可以看出，架空輸電線中點的感應(yīng)過電壓峰值最大，端點處的感應(yīng)過電壓峰值最小。從程序運行的其他計算結(jié)果來看，隨著觀測點從中點向兩端移動，感應(yīng)過電壓的峰值不斷的下降。此外，輸電線路中點處的電壓一直處于正極性，而線路上其他各點均出現(xiàn)雙極性，即在開始階段，電壓先為負值，隨后，電壓極性反轉(zhuǎn)為正值；并且隨著觀測點向兩端移動，雙極性現(xiàn)象更加明顯，即感應(yīng)過電壓初始負峰值越來越大。

從圖5中可以看出，在架空輸電線中點處，入射電壓和散射電壓的幅值相近，但是極性相反，對感應(yīng)過電壓的作用相近；而隨著觀測點向輸電線端點處移動，由垂直電場激勵產(chǎn)生的入射電壓對感應(yīng)過電壓的貢獻越來越微弱，在線路端點處，感應(yīng)過電壓主要是由散射電壓控制的。這是因為從傳輸線的中點到兩端點，入射電壓的衰減很大，到端點時，入射電壓的作用很小，感應(yīng)過電壓也就主要由散射電壓控制了。從圖中我們還可以發(fā)現(xiàn)，入射電壓始終是正極性的，而散射電壓可能是正極性的也可能是負極性的或者雙極性的。這是由于散射電壓由水平電場激勵產(chǎn)生，在有限土壤電導(dǎo)率情況下，離地面h高度處的水平電場表達為其中H∞是無限地壤電導(dǎo)率情況下的磁場，w是衰減因子，Z是阻抗，是無限土壤電導(dǎo)率情況下h高度處的水平電場，當(dāng)h=0時，項為0。將水平電場進行拆分，為正值，為負值，初始時刻，影響較大，后期則由控制，這樣，水平電場可能出現(xiàn)雙極性，散射電壓也就出現(xiàn)雙極性現(xiàn)象。

3 結(jié)語

筆者立足于雷電電磁脈沖對架空輸電線路的影響，選取雙Heidler函數(shù)作為雷電通道基電流，MTLL模型作為回擊模型，利用FDTD算法和Agraw?al耦合模型，分析計算了架空輸電線中點和端點處感應(yīng)入射、散射電壓變化特征。具體模擬結(jié)果表明：1）在線路中心位置處，模擬出的雷擊過電壓波形呈正極性變化趨勢，且中心位置處的過電壓峰值要大于線路端點位置峰值，在端點處雷擊感應(yīng)過電壓波形呈現(xiàn)雙極性，從中點位置處向斷點處移動，雙極性變化特征越顯著；2）在架空線中點位置處，散射、入射電壓峰值相接近，但波形極性相反。在線路端點位置處的耦合過電壓主要為散射電壓。這主要是由于從架空線中點到端點，入射電壓的衰減很大。

[1]歐陽雙.復(fù)雜地表下雷電電磁環(huán)境及架空輸電線感應(yīng)過電壓模擬研究[D].南京：南京信息工程大學(xué)，2012.OUYANG shuang.Simulation of lightning electromagnetic environment and induced overvoltages in overhead trans?mission lines under complex surface[D].Nanjing：Nanjing University of Information Science&Technology，2012.

[2]譚湘海.輸電線路的防雷設(shè)計[D].長沙：湖南大學(xué)，2004.TAN Xiangmei.Lightning protection design of transmis?sion line[D].Changsha：Hunan University，2004.

[3]方瑜.配電網(wǎng)過電壓[M].北京：水利電力出版社，1994.FANG Yu.Distribution network overvoltage，[M].Beijing：Hydraulic Power Press，1994.

[4] BRUCE C E R，GOLDE R H.The lightning discharge[J].Philosophical Magazine，1975，88（192）:487-505.

[5] NORTON K A.The Propagation of Radio Waves over the Surface of the Earth and in the Upper Atmosphere[J].1936，25（9）:1203-1236.

[6] COORAY V.Effects of propagation on the return stroke ra?diation fields[J].Radio Science，1987，22（5）:757-768.

[7] COORAY V.Horizontal fields generated by return strokes[J].Radio Science，2016，27（4）:529-537.

[8] RUBINSTEIN M.An approximate formula for the calcula?tion of the horizontal electric field from lightning at close，intermediate，and long range[J].IEEE Transactions on Electromagnetic Compatibility，1996，38（3）:531-535.

[9]李福壽.電力系統(tǒng)過電壓計算[M].北京：水利電力出版社，1986.LI Fushou.Calculation of overvoltage in power system[M].Beijing：Hydraulic Power Press，1986.

[10]文習(xí)山，彭向陽，解廣潤.架空配電線路感應(yīng)過電壓的數(shù)值計算[J].中國電機工程學(xué)報，1998，18（4）：299-301.WEN Xishan，PENG Xiangyang，JIE Guangrun.Numeri?cal calculation of induced overvoltages in overhead distri?bution lines[J].Proceedings of The Chinese Society for Electrical Engineering，1998，18（4）：299-301.

[11]莫付江，陳允平，阮江軍.架空輸電線感應(yīng)過電壓耦合機理及計算方法分析[J].電網(wǎng)技術(shù)，2005，29（6）:72-77.MO Fujiang，CHEN Yunping，RUAN Jiangjun.Analysis of coupling mechanism and calculation method of induced overvoltage in overhead transmission line[J].Power System Technology，2005.29（6）:72-77.

[12]任合明，周璧華，余同彬，等.雷電電磁脈沖對架空電力線的耦合效應(yīng)[J].強激光與粒子束，2005，17（10）:1539-1543.REN Heming，ZHOU Bihua，YU Tongbin，et al.Cou?pling effects of lightning electromagnetic pulse on over?head power lines[J].High Power Laser and Particle Beams，2005，17（10）:1539-1543.

[13]蘭海燕.傳輸線耦合雷電電磁波過程的分析方法[J].電瓷避雷器，2017（1）:7-13.LAN Haiyan.Analysis method of transmission line cou?pling lightning electromagnetic wave[J].Insulators and Surge Arresters，2017（1）:7-13.

[14]HEIDLER F.Analytische blitzstromfunktion zur LEMPBerechnung[J].paper presented at the 18th international Conference on lightning protection，1985.

[15]唐力，李海明，高晉文，等.雷電流建模方法及其對架空線路耐雷性能分析影響的研究[J].電瓷避雷器，2015（2）:37-43.TANG Li，LI Haiming，GAO Jinming，et al.Research on modeling method of lightning current and its influence on the lightning performance of overhead line[J].Insulators and Surge Arresters，2015（2）:37-43.