基于FDTD混合算法的雷擊桿塔的浪涌分析

谷建宏+廖成+葉志紅+楊丹

摘 要： 由于傳統的電磁暫態分析程序ATP?EMTP在對雷電直擊桿塔的浪涌進行估計時不夠準確，故將傳統FDTD方法與細導線FDTD方法結合起來，對雷擊桿塔的過程進行建模，模擬和分析雷擊桿塔引起的浪涌信號在塔上三相電線上耦合產生的過電流響應，同時分析了絕緣子的擊穿電壓和擊穿時間。研究結果表明，基于FDTD方法雷擊桿塔的線上電流比ATP?EMTP要小，可以有效地估計輸電線路和絕緣子的耐雷水平。

關鍵詞： 細導線FDTD； 桿塔結構； 過電流響應； 絕緣子； 雷擊桿塔； 浪涌分析

中圖分類號： TN99?34； TM863 文獻標識碼： A 文章編號： 1004?373X（2018）05?0116?04

Abstract： Since the traditional electromagnetic transient analysis program ATP?EMTP can′t accurately estimate the surges of the lightning striking tower， the traditional FDTD method and thin?wire FDTD method are combined to model the process of the lightning tower， simulate and analyze the overcurrent response generated by three?phase electric wires coupling of surge signal induced by the lightning tower， and analyze the breakdown voltage and breakdown time of the insulator. The research results show that the FDTD method has smaller line current of lightning tower than ATP?EMTP， and can effectively estimate the lightning tolerance level of the transmission line and insulator.

Keywords： thin?wire FDTD； tower structure； overcurrent response； insulator； lightning tower； surge analysis

0 引 言

電力設施中產生的浪涌信號估計是在電力系統的雷電保護設計中的關鍵步驟。在電力設備的雷擊浪涌分析中，使用電磁暫態分析軟件已經大大提高了準確性。但是電磁暫態分析軟件仿真是基于一個粗略的簡易模型，雷擊的空間傳播過程分析仍有待研究，比如說，在地面上浪涌沿著導體不均勻傳播，垂直于或傾向于導體傳播，仍亟待解決。這就衍生了沿著傳輸桿塔、沿著輸電線路、沿著地線傳播的浪涌分析。

在國內的雷電保護設計中的雷擊分析大部分是通過電磁暫態分析軟件分析桿塔受到雷擊并且產生回擊的情況[1?2]。由于雷擊桿塔時傳播的電磁波為非橫電磁波[3]，對尺寸較大、結構相對復雜的輸電線路桿塔來說，使用EMTP軟件仿真結果會有一定誤差。采用FDTD方法從電磁場的角度對雷擊桿塔進行仿真，可以比較精確地模擬出桿塔內部的瞬態過程，可以更加有效地估計輸電線路和桿塔絕緣子的耐雷水平。

本文對500 kV同塔雙回線路桿塔的簡化模型在標準波形的雷擊條件下的雷電浪涌現象使用傳統FDTD與細導線FDTD結合的方法進行數值分析。分析了回閃時刻絕緣子的擊穿電壓值以及傳輸線上的電流響應，同時，與在相同條件下電路分析的結果進行對比，在分析結果的基礎上從雷電耐受電壓的角度進行評估。

1 細導線FDTD基本理論

FDTD方法是一種在空間和時間的推演中進行差分迭代來求解Maxwell方程的一種方法。考慮到桿塔結構和輸電線路均可近似為細導線結構，本文結合準靜態場推導的細導線處理方法對雷擊桿塔輸電線上過電流的問題進行建模和仿真。由于雷電流源的頻率相對較低[4]，為滿足FDTD算法的數值穩定性和色散性，可以將仿真區域的網格設置得比較大。桿塔的金屬結構相對網格的大小滿足細導線近似的使用條件，同時，輸電線路和避雷線的粗細相對網格大小同樣可以近似為細導線結構進行處理。

由于桿塔結構、輸電線和避雷線都為細導線結構，在雷擊桿塔的仿真過程中，最重要的就是細導線等效模型的正確性。

本文中對細導線的等效處理參考Taku Noda等人在2002年提出的基于準靜態場推導的細導線處理方法[3，5]。Taku Noda修正了細導線橫截面平面方向4個元胞中的介電常數，乘上一個常數將周圍4個元胞的磁導率除了一個常數同時將細導線徑向位置上的電場強制賦值為0，其中為細導線半徑，為空間步長。基于準靜態場推導的細導線等效示意圖見圖1。

2 建模過程

2.1 雷電通道建模

本文用上文所述的細導線結構來模擬雷電通道[6]，雷電通道建模側視圖見圖2。在通道底部添加一個集總電流源，在雷電通道的周圍一個網格設置相對介電常數為9的電介質來模擬雷電流的回擊速度。

在雷電通道的底部添加集總電流源，需要對電流源位置周圍4個方向的磁場進行強制賦值來滿足：

式中：為空間步長，單位為m；單位為A。

在仿真過程中，通常簡單的將4個磁場分量設置為絕對值相等：

2.2 桿塔建模

在FDTD仿真中，最重要的就是對桿塔進行建模。輸電桿塔結構由很多細導線結構組成，這些細導線結構中除了小部分垂直和平行于地面的，大部分并非水平豎直，而是傾斜的。而FDTD無法對傾斜的細導線進行直接處理，所以首先要對桿塔結構進行階梯近似處理再在FDTD中進行建模。endprint

首先將端點位置不在YEE元胞角點上的直線經過近似處理使其兩端點位于元胞角點處。然后對傾斜直線進行階梯近似處理，在三維空間條件下選擇階梯近似路徑需要遵循以下兩個原則：路徑必需沿從起點開始到終點的方向行進；所選取的路徑必須要最靠近原本的直線段。根據這兩個條件在Matlab軟件中進行編程，提取CAD圖中桿塔的所有線段始終點坐標進行處理，最后得到所有階梯近似后的線段始終點的坐標，如圖3所示。

2.3 三相線及避雷線

在圖3中的四個橫擔上的線路從上到下依次為避雷線、上相線、中相線、下相線。避雷線直接與桿塔相連接，而各相線則通過絕緣子與桿塔相接。

在高壓輸電線路中，通常使用分裂導線形式的輸電線。由于FDTD分析雷擊桿塔仿真的空間步長通常比較大，而分裂導線中各個導線的距離比較小，通常對分裂導線進行等效處理為半徑為的單根導線[3]。計算公式為：

式中：為分裂導線的根數；為實際導線半徑；為分裂導線間距。

3 仿真條件設置及結果

在FDTD仿真過程中，設置的網格大小為1 m×1 m×1 m，整個空間大小為100 m×100 m×500 m，整個空間中有100×100×500個元胞，邊界條件設置為PML邊界條件，時間步長為1 ns。在空間的下方為厚度50 m的土壤層，其相對介電常數為10，相對磁導率為1。桿塔模型上方設置為長369 m，半徑為0.1 m的雷電通道，在通道周圍的1 m空間內設置相對介電常數為9的無損介質，以此來模擬雷電流回擊的速度。桿塔上設置橫向長為100 m的輸電線路，根據分裂導線參數的等效原則，可等效為單根半徑為0.21 m的導線。在桿塔與雷電通道之間加入電流源，波形為雙指函數，幅值為120 000 A，波頭時間為2.6 μs，波尾時間為50 μs。

FDTD仿真中絕緣子兩端電壓的變化如圖4所示，其中上相線和中相線的絕緣子分別在322 ns和375 ns發生擊穿，擊穿電壓均約為1 400 kV，而下相線的絕緣子則在824 ns才發生擊穿，擊穿電壓約為1 100 kV。圖5為FDTD仿真中避雷線及各相線上的電流變化曲線，可以看出各相線上的電流在擊穿后顯著地增加并在短時間內達到峰值。

4 ATP仿真結果對比分析

本節通過在EMTP軟件[7]中建立基本的桿塔模型和線路在相同的條件下進行仿真，將得到的結果與FDTD仿真結果進行對比分析，如圖6，圖7所示。

基于桿塔的有損多波阻抗模型進行建模的[8]EMTP軟件的雷擊桿塔瞬態響應的仿真結果在一定程度上比較接近FDTD的仿真結果。圖6為避雷線上雷擊過電流的變化對比圖，FDTD結果和ATP?EMTP程序的避雷線仿真結果比較接近，FDTD的峰值電流為12 933 A，ATP?EMTP的峰值電流為13 136 A。圖7為桿塔上相線的雷擊過電流變化對比圖，FDTD的峰值電流為8 362 A，ATP?EMTP的峰值電流為8 637 A，FDTD仿真結果比EMTP仿真結果低275 A，并且在波形之后的衰減中FDTD仿真結果始終要比ATP?EMTP仿真結果低900 A左右。產生這種現象的原因主要是流動在傳輸桿塔和架空線路的雷電流的電磁感應現象，ATP?EMTP軟件仿真并未考慮空間的電磁波傳播問題，所以部分能量逸散到空間中，造成FDTD分析結果偏小。

圖8為ATP?EMTP中的絕緣兩端的電壓變化圖，與FDTD仿真結果相比，ATP?EMTP中絕緣子的擊穿時間比FDTD仿真結果要早一些，FDTD仿真結果的絕緣子雷電耐受電壓更高一些，如表1所示。

5 結 語

本文將傳統FDTD方法與細導線FDTD方法結合起來，對雷擊桿塔的過程在三維空間中進行簡易的建模，仿真和分析了雷擊桿塔引起的浪涌信號在塔上三相電線上耦合產生的過電流響應，同時分析了絕緣子的擊穿電壓和擊穿時間。研究結果表明，通過FDTD分析得到桿塔上線上的電流波形，要比ATP?EMTP仿真得到的電流值要小，這歸結于流動在傳輸桿塔和架空線路的雷電流的電磁感應現象，500 kV同塔雙回路桿塔絕緣子在雷電回擊時有著很高的雷擊耐受能力。

參考文獻

[1] 張志勁.500 kV同桿雙回輸電線路耐雷性能研究[D].重慶：重慶大學，2002.

ZHANG Zhijin. Analysis on the lightning protection performance of 500 kV double?circuit transmission line [D]. Chong?qing： Chongqing University， 2002.

[2] 劉渝根，劉緯，陳先祿.500 kV變電站雷電侵入波研究[J].重慶大學學報（自然科學版），2000（3）：17?19.

LIU Yugen， LIU Wei， CHEN Xianlu. Research on lightning intruded wave for 500 kV substations [J]. Journal of Chongqing University （natural science edition）， 2000（3）： 17?19.

[3] NODA T. A tower model for lightning overvoltage studies based on the result of an FDTD simulation [J]. IEEE transactions on power and energy， 2007， 127： 379?388.

[4] 丁美新，李慧峰，朱子述，等.雷電流波形的數學模型及頻譜仿真[J].高電壓技術，2002（6）：8?10.endprint

DING Meixin， LI Huifeng， ZHU Zishu， et al. The mathematical model and frequency spectrum simulation of lightning current [J]. High voltage engineering， 2002（6）： 8?10.

[5] NODA Taku， TATEMATSU Akiyoshi， SHIGERU Yokoyama. Improvements of an FDTD?based surge simulation code and its application to the lightning overvoltage calculation of a transmission tower [J]. Electric power systems research， 2007， 77（11）： 1495?1500.

[6] 李韋霖，楊琳，李慧，等.雷電通道模型研究與應用.現代電子技術，2010，33（5）：174?178.

LI Weilin， YANG Lin， LI Hui， et al. Research and application of lightning channel model [J]. Modern electronics technique， 2010， 33（5）： 174?178.

[7] AMETANI A， KAWAMURA T. A method of a lightning surge analysis recommended in Japan using EMTP [J]. IEEE transactions on power delivery， 2005， 20（2）： 867?875.

[8] ISHII M， KAWAMURA T， KOUNO T， et al. Multistory transmission tower model for lightning surge analysis [J]. IEEE tran?sactions on power delivery， 1991， 6（3）： 1327?1335.endprint