特高壓輸電線路下建筑物鄰近區(qū)域電場的數(shù)值仿真

裴春明，方福歆，王 璇，孔海洋，錢 航，關煥梅

(1.中國電力科學研究院武漢分院，湖北 武漢 430074；2.武漢大學電氣工程學院，湖北 武漢 430072；3.武漢大學計算機學院，湖北 武漢 430072）

0 引言

隨著輸電線路電壓等級的不斷提高，輸電過程中的電磁問題愈發(fā)受到人們的關注[1-3]。近年來國內(nèi)外學者紛紛進行相關研究，運用實驗和數(shù)值仿真的方法，對超高壓及特高壓線路下方的工頻電場進行測量和估算[4-8]。就數(shù)值仿真而言，主流的分析方法包括模擬電荷法和有限元法。模擬電荷法多適用于線路附近無建筑物，空間介質(zhì)分布均勻的情形；線路附近存在建筑物時，采用模擬電荷法計算畸變電場將產(chǎn)生較大誤差，此時多使用有限元法[9]。而無論何種情形，研究關注的核心區(qū)域多集中于輸電線路弧垂最低點附近。文獻[10]同時分析了長直導線模型和考慮弧垂的線路模型，得到兩種模型下電場分布的異同，但沒有考慮輸電線路附近有建筑物時的情況；文獻[11]提出了改變超高壓輸電線路鄰近建筑物的棱邊形狀，繼而改變其曲率半徑的方法以減小畸變場強，卻沒有給出具體的建筑物改變形狀及改變程度。基于此，本文欲探討如下兩個問題：（1）在使用有限元法對特高壓輸電線路下建筑物鄰近區(qū)域的電場進行仿真計算時，若導線模型不采用實際的懸鏈線模型，而采用簡化后的直導線代替，是否會影響所關注區(qū)域電場強度的分布結果。（2）建筑物屋頂以及陽臺的棱角處若通過倒角處理進行鈍化，會否影響其周圍畸變電場的分布；改變倒角深度，會否改變畸變場強的大小。

1 有限元法簡述

有限元法(Finite Element Method，F(xiàn)EM)的核心思想是將一個復雜的連續(xù)場域，分割成眾多簡單的子場域，進而針對每一個子場域建立不同的方程，最后將它們累加組合得到原始場域的解。主要求解過程如下[12]：（1）列出與偏微分方程邊值問題等價的條件變分問題；（2）將區(qū)域作三角形或者四面體單元剖分，并在單元中構造出線性插值函數(shù)；（3）將能量泛函的極值問題轉化為能量函數(shù)的極值問題，建立線性代數(shù)方程組；（4）求解線性代數(shù)方程組。本文探討的問題為輸電線路的工頻電場，故而在求解結果過程中需用到泊松方程

而有限元法的插值基函數(shù)一般如下

利用變分法將(1)式離散化，可得到其離散矩陣方程組

將式(2)分別帶入其中可以得到：

將式(4)、(5)代入式(3)中可求得U，其中U是各點電位的列向量，進而可通過下式計算得到各點電場強度E=-?φ,其中φ為各點的電位。

2 模型的建立與求解

2.1 模型參數(shù)

文獻[13]對比了單回輸電線路四種排列方式下的場強，結果表明三角布置形式下的場強值最小。因此本文以正三角排列貓頭型塔為原型搭建模型，導線則選用8分裂LGJ-500/45。由于有限元計算軟件多以四面體結構進行三維模型的剖分，對曲面的處理需用到大量的剖分單元，這將導致運算時間大大增長。因此本文將分裂導線做等效處理，導線的等效半徑求解公式如下式所示[14]，R=(rNAN-1)1N，

式中：r為子導線半徑，N為子導線數(shù)，A為分裂導線半徑。

假設導線沿x軸排布，則其懸鏈線方程如下式所示，式中：σ0為導線最低點的水平應力；g為導線的比載；h為弧垂最低點距地面高度。

簡化的直導線模型則取懸鏈線弧垂的最低點作為其對地高度。從文獻[14]中可知，架空地線對地面上方工頻電場的影響大約為2%左右，故而在確保仿真結果精度的條件下，略去模型中的地線以縮短計算時間。當線路下方鄰近建筑物的材料為鋼筋混凝土時，關注區(qū)域的場強較之其它材料的建筑物最小[15]，因此本文中的模型均使用鋼混結構，相對介電常數(shù)為6.4。鋼筋框架部分均與大地等電位，大地設置為零電位。

圖1 考慮弧垂前后線路模型對比圖Fig.1 The contrast between the model consider sag and the model regardless of the sag

2.2 模型構建

本文使用商業(yè)有限元軟件Ansoft Maxwell進行仿真計算。基于上文設定的參數(shù)，搭建線路模型如圖2所示。直線模型共包含2 068 579個剖分單元，而考慮弧垂的實際導線共包含2 957 726個剖分單元，單元量增加約25%，且單元規(guī)則程度降低，導致運算速度大大下降。圖2所示為建筑物的具體模型，同時給出棱邊倒角之后的建筑物結構。

圖2 棱邊倒角前后建筑物模型對比圖Fig.2 The contrast between the model after chamfering and the model with no chamfer

待求解區(qū)域分別為一樓上方1.5 m處，二樓上方1.5 m處以及建筑物樓頂處。待求解路徑均在各待求解區(qū)域內(nèi)取得，以中相導線為原點，沿導線的法線方向延展100 m。

3 結果分析

3.1 考慮弧垂前后仿真結果對比和分析

分別計算考慮弧垂的線路模型和簡化的直導線模型，得到一樓上方1.5 m處的場強分布結果如圖3所示。從圖中可以看出，考慮弧垂的實際導線模型與簡化直導線模型相比，電場強度較小一些。兩模型下場強的差異主要體現(xiàn)在距線路中相導線20 m及50 m處，即建筑物兩側10 m處，場強值的最大誤差約為7.33%。而對于待求解區(qū)域的其它位置，特別是建筑物所在區(qū)域，兩模型的場強大小幾乎相同，場強的畸變程度也幾乎一致。簡化直導線場強偏大的原因在于，本文設置的直導線模型的對地高度為實際線路的弧垂最低點處，等價于降低了實際線路的對地高度，因此造成場強增大。

圖3 距地1.5m處場強分布對比圖Fig.3 Field intensity distribution on the height of 1.5 m above the ground

同樣可以得到二樓上方1.5 m處的場強分布結果如圖4所示。簡化直導線模型的場強值仍然整體偏大，與考慮弧垂的實際導線模型相比，最大誤差約為9.76%，整體趨勢同上文相似。值得注意的是，此求解區(qū)域覆蓋了二樓的陽臺部分。而對于陽臺處(距線路中相40 m處)的畸變場強，簡化直導線模型卻較之實際導線模型偏小。究其原因，主要是由于考慮弧垂的實際導線模型本身彎曲，從而造成空間電場分布不均勻，在曲率半徑較小的地方這種不均勻將會更加顯著，最終導致陽臺處場強的畸變較之簡化直導線模型更為嚴重。

圖4 二樓上方1.5 m處場強分布對比圖Fig.4 Field intensity distribution on the height of 1.5 m above the second floor

樓頂處場強分布結果（圖5）恰可以佐證以上的分析。由于建筑物樓頂處的場強值很大，兩模型的差異不如前述明顯，但仍可以看出簡化直導線模型的場強值依然整體偏大。與此同時，建筑物樓頂兩側棱邊處(距線路中相31 m及38 m處)的畸變場強對比表明，簡化直導線模型仍較之考慮弧垂的實際導線模型偏小。

圖5 樓頂處場強分布對比圖Fig.5 Field intensity distribution on the roof

3.2 棱邊倒角前后仿真結果對比和分析

對建筑物模型的樓頂處及陽臺處進行倒角處理，倒角深度分別為5 cm、10 cm、15 cm和20 cm。仿真結果如圖6和圖7所示。倒角深度分別為0 cm和20 cm時樓頂處的場強分布云圖如圖8所示。

從上述圖中可以看出，建筑物棱邊經(jīng)倒角處理后，場強的整體分布沒有發(fā)生變化，五種情況下的場強分布曲線近乎重疊在一起，區(qū)別僅在于各棱邊處的畸變場強發(fā)生了變化。樓頂及陽臺處的棱邊經(jīng)不同程度倒角處理后，最大畸變場強出現(xiàn)的位置均發(fā)生變化，且最大畸變場強的數(shù)值也各不相同。

圖6 陽臺處棱邊倒角前后場強分布圖Fig.6 Field intensity distribution of the balcony before and after the edge chamfer

圖7 樓頂處棱邊倒角前后場強分布圖Fig.7 Field intensity distribution of the roof before and after the edge chamfer

圖8 樓頂處棱邊倒角前后場強分布云圖Fig.8 Nephogram of Field intensity distribution of the roof before and after the edge chamfer

圖9指出了最大畸變場強值隨棱邊倒角深度不斷增加時的變化規(guī)律。可以看出，陽臺處的最大畸變場強值隨倒角深度的增加變化不大，僅倒角深度為20 cm時，下降幅度較大，同未倒角模型相比，下降約20%。而樓頂處最大畸變場強則出現(xiàn)較大波動，且與未倒角建筑物相比，倒角后的建筑物畸變場強反而均有所增加，最大增幅達20%。造成此種現(xiàn)象的原因可能為棱角經(jīng)倒角操作后雖已鈍化，但棱角數(shù)量增加，畸變電場的相互疊加將導致場強值增大。

圖9 最大畸變場強隨倒角深度的變化曲線Fig.9 The largest distortion field intensity curve along with the change of chamfering depth

4 結語

本文針對特高壓輸電線路下建筑物鄰近區(qū)域電場的數(shù)值仿真計算，探討了其中尚未透徹分析的兩個問題，得到如下結論：

(1)簡化直導線模型同考慮弧垂的實際線路模型相比，非畸變區(qū)域的場強值均偏大，樓頂棱邊及陽臺等畸變區(qū)域的場強值均偏小。在求解精度要求較高時，不應以直導線模型代替實際線路模型進行計算。而對于求解精度要求較低，著重關注場強分布的情況，因簡化直導線模型的計算較之實際線路模型可節(jié)省運算時間約60倍，此時可考慮使用簡化直導線模型進行計算，有效提升運算效率。

(2)倒角后的模型與未倒角模型相比，陽臺處的畸變場強有所降低，但仍沒有降至國標限值之內(nèi)；樓頂處的場強則均有所增加，畸變更加嚴重。因此通過倒角的方式鈍化建筑物棱邊不能有效改善建筑物特定區(qū)域的電場環(huán)境。此時應考慮采取架設屏蔽線等措施改善畸變電場。

[參考文獻]（References）

[1]鄔雄.特高壓輸電的電磁環(huán)境問題研究[J].電力系統(tǒng)通信，2006，27(6)：1-5.Wu Xiong.Research on electromagnetic environment of ultra-high voltage grid[J].Telecommunications for Electric Power System,2006，27(6)：1-5.

[2]邵方殷.1 000 kV特高壓輸電線路的電磁環(huán)境[J].電網(wǎng)技術，2007，31(22)：1-6.Shao Fangyin.Electromagneticenvironmentof 1 000 kV UHV transmission line[J].Power System Technology,2007，31(22)：1-6.

[3]黃道春，阮江軍，文武，等.特高壓交流輸電線路電磁環(huán)境研究[J].電網(wǎng)技術，2007，31(1)：6-11.Huang Daochun,Ruan Jiangjun,Wen Wu,et al.Study on electromagnetic environment of UHV AC transmission lines[J].Power System Technology,2007，31(1)：6-11.

[4]Krajewski W.BEM analysis of electric field excited by overhead HV lines erected in built-up areas[J].IEE Proc of Science，Measurement and Technology，1997,144∶81-86.

[5]黃道春，阮江軍，余世峰，等.特高壓緊湊型輸電線路工頻電場強度計算[J].高電壓技術，2006，32(7)：69-71.Huang Daochun,Ruan Jiangjun,Yu Shifeng，et al.Calculation of the power frequency electric field intensity of UHV compact transmission lines[J].High Voltage Engineering,2006，32(7)：69-71.

[6]許丹，馬星河，王晨輝，等.500 kV輸電線路電場強度測量與計算分析[J].高壓電器，2013，49(7)：25-28.Xu Dan,Ma Xinghe,Wang Chenhui,et al.Electric field measurement and computational analysis on 500 kV transmission lines[J].High Voltage Apparatus,2013，49(7)：25-28.

[7]彭迎，阮江軍.模擬電荷法計算特高壓架空線路3維工頻電場[J].高電壓技術，2006，32(12)：69-73.Peng Ying,Ruan Jiangjun.Calculation of three dimensionalharmonic electric field around ultra-high voltage overhead line based on the charge simulation method[J].High Voltage Engineering,2006，32(12)：69-73.

[8]黃峰，孫苓生.跨越建筑物的高壓輸電線三維電場仿真[J].華東電力，2007，35(12)∶1165-1167.Huang Feng,Sun Lingsheng.3-D simulation of electric fields for high voltage transmission lines over buildings[J].East China Electric Power,2007，35(12)∶1165-1167.

[9]趙志斌，董松昭，鄧新麗，等.特高壓同塔雙回交流輸電線路鄰近建筑物時的安全距離[J].高電壓技術，2013，39(12)∶2943-2950.Zhao Zhibin,Dong Songzhao,Deng Xinli,et al.Safety distance of UHV double-circuit AC transmission lines on one tower adjacent to buildings[J]. High Voltage Engineering, 2013，39(12)∶2943-2950.

[10]蘭生，張振興，原永濱.考慮弧垂的交流特高壓輸電線三維電磁場[J].電機與控制學報，2012，16(12)∶42-46.Lan Sheng,Zhang Zhenxing,Yuan Yongbin.The 3D electromagnetic field of AC EHV transmission lines with sag[J].Electric Machines and Control,2012，16(12)∶42-46.

[11]汪泉弟，羅揚，楊帆，等.超高壓輸電線下建筑物臨近區(qū)域三維電場仿真[J].高壓電器，2013，49(10)∶1-6.WangQuandi,LuoYang,YangFan,etal.Three-dimensionalsimulation ofelectric field around building under EHV transmission lines[J].High Voltage Apparatus,2013，49(10)∶1-6.

[12]趙博，張洪亮.Ansoft 12在工程電磁場中的應用[M].北京：中國水利水電出版社，2010∶53-55.Zhao Bo,Zhang Hongliang.Application of Ansoft 12 in engineering electromagnetic field[M].Beijing∶China Water Conservancy and Hydropower Press,2010∶53-55.

[13]程煒，劉黎剛，張艷芳，等.特高壓輸電線路工頻電場的數(shù)值仿真研究[J].高壓電器，2012，48(2)：1-6.Cheng Wei,Liu Ligang,Zhang Yanfang,et al.Numerical simulation study on power frequency electric field of UHV AC transmission lines[J].High Voltage Apparatus,2012，48(2)：1-6.

[14]曾慶禹.特高壓輸電線路地面最大工頻電場強度和導線最大弧垂特性[J].電網(wǎng)技術，2008，32(6)：1-7.Zeng Qingyu.Study on characteristics of maximum power frequency electric field gradient on ground and maximum sag of bundle conductor for UHV AC transmission line[J].Power System Technology,2008，32(6)：1-7.

[15]Yang Bin，Wang Shuhong，Wang Qiang，et al.Simulation and analysis forpowerfrequency electric field of building close to power transmission lines[C]//IEEE International Symposium on Electromagnetic Compatibility(EMC).New York∶IEEE,2014∶451-454.