基于有限元法的同走廊交直流線路混合電場的計算研究

王建林，靳 斌，李博文，李 興，熊 明

(西華大學 電氣與電子信息學院，四川 成都 610039)

基于有限元法的同走廊交直流線路混合電場的計算研究

王建林，靳 斌，李博文，李 興，熊 明

(西華大學 電氣與電子信息學院，四川 成都 610039)

特高壓直流線路與超高壓交流線路同走廊架設將會有效地減少架設輸電線路的土地成本，但卻帶來了新的電場計算問題.為了研究二者的混合電場，首先使用有限元分析軟件Comsol Multiphysics計算了±800kV特高壓直流線路下方的合成電場及500kV超高壓交流輸電線路下方的工頻電場,進而計算了混合電場，最后重點探討了導線布局對混合電場的影響.結果表明，基于Comsol Multiphysics的混合電場計算是合理可靠的；合理的接近距離和導線最小對地高度能分別有效地減少混合區域場強和混合電場最大值.

特高壓直流線路；Comsol；同走廊；混合電場；有限元

近年來，為滿足持續增長的用電負荷需求，需要建設更多的輸電線路.與此同時，原有的土地資源變得越來越緊缺，尤其是在經濟發達地區，新輸電走廊的獲取變得越來越困難，輸電工程中勢必會出現交直流線路同走廊或同塔架設的情況,雖然這樣可以節約土地資源并提高單位走廊的輸送容量，但也帶來了新的電磁環境問題[1].因此對其混合電場的計算就變得很有必要了.

近年來，中外學者基于ANSYS、Ansoft等有限元分析軟件對輸電線路的電場分布作了不少的研究.但使用Comsol Multi- physics來研究輸電線路電場分布的還不多.文獻[2]基于Comsol研究了特高壓直流輸電線路的電場分布，對使用Comsol求解輸電線路的電場分布有一定的借鑒意義.本文中，筆者首先基于Comsol分別計算了特高壓直流輸電線路和超高壓交流輸電線路的電場，并在此基礎上計算了二者的混合電場，最后分析了不同的交直流線路接近距離和導線最小對地高度對混合電場分布的影響規律，以期為交直流線路同走廊建設和運行提供參考.

1 直流線路電場的計算

1.1 幾何模型

以同塔雙回同相序交流輸電線及一條雙極直流輸電線路同走廊并行為例(圖1)，直流導線型號為6×LGJ-720/50，分裂間距為0.45m，子導線半徑為1.68cm，極間距22m，高度為21m.交流導線型號為4×LGJ-400/35，分裂間距為0.45m，弧垂最低點對地高度為15m.交、直流輸電線路接近距離D為70m，為簡化計算不考慮避雷線的影響.

圖1 同走廊交、直流輸電線路的幾何位置示意圖

1.2 直流線路電場的數學控制方程及邊界條件

高壓直流輸電線路正常運行時，導線表面電場強度高于電暈起始電場強度會產生電暈[3],由電暈產生的帶電離子所形成空間電場將會與輸電線路上的電荷產生的標稱電場疊加并形成合成電場.目前，計算分析直流輸電線路的合成場強時主要使用的方法有：半經驗公式法，基于Deustch假設的解析法和二維數值計算方法[4].本文采取在Comsol上建立直流輸電線路模型[5-6]，并結合經驗公式的方法來計算直流輸電線路的合成場強.其控制方程為

(1)

(2)

(3)

(4)

(5)

(6)

式中:φ為標量電位；ρ+、ρ-分別為正、負空間電荷密度；E為合成電場的電場強度；J+和J-分別為正、負離子流密度；K+、K-分別為正、負離子遷移率；R為離子復合系數;ε0為空氣介電常數；e為電子電量.

邊界條件：

導線表面電壓為其運行電壓φ±=±U;

地面上的電位φ0=0.

1.3 網格剖分

在分析分裂導線的離子流場時通常將將分裂導線用等效半徑導線代替，等效公式為[7]

(7)

式中：req為等效半徑；Rf為分裂間距；n為分裂數；r為子導線半徑.

以本文的直流輸電線為例，采用公式(7)計算可知req為0.38m.

考慮到精度要求和盡量減少系統運算量，選擇2D場進行處理，使用靜電穩態求解器對模型進行求解，所設置的求解域為半徑為50m的半圓形域(圖2).

為了提高計算的精度，筆者截取求解域邊緣附近5m寬的半圓環狀帶作無限元處理，這樣可將被人工截斷的物理域作1 000倍放大，即實際物理寬度變為5×103m；有限元中，在進行網格剖分時，理論上網格數目越多計算精度越高，但是過密的網格剖分將會占用大量的計算內存并降低計算速度，甚至于導致內存溢出[8].綜合考慮，本文設置自由剖分三角形網格的尺寸大小為“細化”，最大單元尺寸6.7m，最小單元尺寸0.03m, 從較小的單元到較大單元的最大單元生長率為1.3.

圖2 網格剖分圖

1.4 計算結果分析

1.4.1 標稱場強計算

從圖3可以看出，兩極導線的標稱場強以線路中心對稱分布，最大值出現在輸電線路外側為24.2kV/m，最小值出現在兩極線路中心為0.圖3中，每極導線最中間的空白部分為輸電線路的內徑，其電場值為0，通過高斯通量定理進行計算也可得到相同結果.

圖3 等效直流輸電線路標稱電場分布

求解直流線路的電場時，一般將電場評價參考點設在地面上0m或者1.5m處，這兩處的電場值差別不大.考慮到混合電場的求解，在下文的合成場強計算中，筆者將距地1.5m處作為計算的參考點.

1.4.2 合成場強計算

對于使用半經驗法求解直流線路的合成場強，法國電力科學院和美國EPRI等科研機構都作過相應的研究并得到了不同形式的求解公式.其中，最具代表性的是美國EPRI提出的方法.基于大量模擬試驗，EPRI認為直流線路下的電場有兩種極限情況：一是沒有電暈時，僅有導線上電荷決定的靜電場(標稱電場)；二是飽和電暈時僅由空間電荷決定的電場.計算直流線路下的空間電場時，首先計算出上述兩種極限情況的電場分布，在此基礎上再計算未飽和電暈放電時的合成場強.并提出了計算合成電場的經驗公式[4,8-9]：

(8)

式中:E、E0、ED分別為實際線路合成電場、標稱電場、飽和電暈電場；Ke為常數，查相關曲線可知；gmax、g0分別為導線表面最大場強和起暈場強，gmax可通過仿真計算得出，g0查表可知.

圖4 直流輸電線路合成場強分布

圖4為距地1.5m處的標稱場強和合成場強的分布圖.從圖中可以看出，合成場強的最大值是標稱場強的2倍以上，其最大值約為27.2kV/m.可見，空間電荷的存在對于合成場強的增大有明顯的作用.

2 交流線路電場的計算

2.1 計算分析

目前計算交流線路工頻電場的方法主要有模擬電荷法和有限元法等.雖然交流輸電線路的工頻電場是交變電場，但由于工頻頻率很低，且本文研究范圍遠小于200m，在討論工頻電場時可忽略工頻磁場的影響而將其視作準靜電場處理[9-10]，這就為使用ComsolMultiphysics求解交流工頻電場提供了可行性.將輸電線路視為理想的三相對稱，根據圖1的幾何布局圖，在ComsolMultiphysics上建立同塔雙回同相序的交流輸電線路模型.與上文求解直流線路電場類似，對分裂導線進行等效，網格剖分選為“細化”，求解域為半徑50m的半圓形域，不同的是，選取的求解器為靜電瞬態求解器.在施加激勵源時，通過調用函數的方式實現，以回路一A相輸電線路為例：

(9)

相較于A相輸電線路，B相輸電線路和C相線路的激勵源幅值不變，僅在相位上依次滯后120°.

2.2 計算結果

圖5為輸電線路電場強度幅值的分布圖，從圖中可看出來，電場強度兩側呈對稱分布，最大值出現在導線外側為514kV/m，離導線越遠電場強度越小，當距地只有1.5m左右時，中心處的電場最大值大概為后文所提的5.7kV/m.圖6為回路一的A相輸電線路的電場分布，從圖中可看出，由于此時不存在其他導線的屏蔽作用，電場最大值為532kV/m，距離A相輸電線路越遠，電場強度下降幅度越大.

圖5 等效交流輸電線路電場強度幅值分布

圖6 回路一的A相輸電線周圍電場分布

圖7為對仿真結果處理后所作距地1.5m處的工頻電場分布圖.從圖中可以看出，隨著距線路中心距離的增加，電場強度先增加后減小，最大電場強度出現在線下距中心10m處左右，其最大值約為6.4kV/m，距離超過40m后，電場強度緩慢變化并接近0.

圖7 離地1.5m處電場強度

3 混合電場的計算

3.1 計算方法

文獻[11-12]指出，在求解混合電場時，由于三相交流輸電線路下方的空間各點電場是旋轉的矢量，且直流線路電場可以視作恒定值，計算時，選擇將交流線路工頻電場和直流線路合成場強直接疊加得到混合電場的方法在工程中是可行的.則交、直流線路同走廊架設時，對空間中任意一點的混合電場強度有：

Eh=Ed+Ea

(10)

式中:Eh為混合電場強度，Ed為直流線路電場，Ea為交流線路電場，其中Ed、Ea上文已經給出.

需要注意的是，求解混合電場時，直流線路電場強度取其絕對值.

3.2 計算結果

以接近距離D的中心作為作圖參考點得到圖8.從圖中可以看出，由于場的疊加作用，相較于直流線路單獨作用時，同走廊區域中靠近交流線路側的電場有明顯的增大，但是混合電場的最大值主要還是由直流線路決定的，最大值27.9kV/m出現在距走廊中心約20m處，且混合電場滿足我國電力行業規定的限值.

圖8 同走廊交、直流輸電線路混合電場分布圖

3.3 結果驗證

目前，對于交直流同走廊的實測數據較少，文獻[13]在試驗室架設兩條架空線來模擬輸電線路，給其中一條導線通以+50kV直流電壓，另一條則通以有效值為36kV的交流電壓，交直流導線距地高度h均為1m，交直流導線間距D為1.8m，直徑d均為3.08mm，將兩條導線的幾何中心作為參考點.從表1的計算結果可以看出，本文的混合電場計算值和試驗值基本一致.在交、直流線路共走廊的地方由于忽略了其相互影響，計算值比試驗值偏大，以距中心0.6m內的數據為例，誤差最大的計算值出現在0.3m處，其誤差為8.9%，屬于工程設計中允許的范圍，通過理論分析，筆者也得到了相似結論.可見本文基于有限元的交、直流輸電線路同走廊時混合電場的計算是合理可靠的，對工程設計和預測具有一定參考價值.

表1 混合電場計算值和試驗值比較

水平距離/m試驗值/kV·m-1本文方法/kV·m-1文獻[10]/kV·m-1-1．213．813．713．6-0．916．316．516．0-0．615．816．015．6-0．315．316．415．8016．818．217．20．318．019．618．50．619．120．018．30．919．519．819．11．216．516．515．9

4 導線布局對混合電場的影響

對于單獨架設的直流輸電線路或交流輸電線路，導線對地高度、導線布置方式，分裂子導線間距和半徑，大氣條件等都會影響其線路下方的電場強度.同走廊架設的交、直流線路輸電線路混合電場的影響因素很多.筆者采用上文所提方法重點研究了改變交、直流輸電線路的接近距離和導線對地高度對混合電場的影響規律.

4.1 不同接近距離對混合電場的影響

在上文的基礎上，改變交、直流輸電線的接近距離，分別計算當接近距離D為50m和90m時的混合電場.計算結果如圖9所示.從圖中可以看出，增加或是減少二者的接近距離，場強最大值基本沒有變化.增加兩者的接近距離，混合區域的場強有一定程度減少；減少兩者的接近距離，混合區域內的場強增加，但是其場強并未超出最大值.交、直流線路共用走廊時, 由于線路內側不必考慮居民區，此時考慮走廊寬度時的原則是：在該公用走廊內最大地面合成場強未超過迭加電場的最大允許值, 走廊邊緣處場強的允許值與交直流線路單獨存在時相同[14].因此在滿足絕緣要求的情況下，適當地減少接近距離，既滿足了工程設計的要求，又能減少走廊寬度，從而有效地提高土地的利用率.對于同走廊架設的特高壓直流線路和超高壓交流直流線路建議接近距離為50m～70m.

圖9 不同接近距離時的混合電場分布

4.2 不同對地高度對混合電場的影響

分別改變交流輸電線路和直流輸電線路的對地高度，分析其對混合電場的影響規律.接近距離D=70m,計算結果如圖10所示，從圖中可以看出，保持直流線路的最小對地高度h=21m不變，當交流線路的最小對地高度從15m增加到18m時，混合電場的最大值基本無變化，混合區域內的場強有小幅度的減少；保持交流線路的最小對地高度h=15m不變，當直流線路的高度從21m增加到24m時，混合區域內的場強基本無變化，混合電場的最大值從27.9kV/m變為25.5kV/m，減幅達8.6%.可見，適當地抬高直流線路線路的最小對地高度可以明顯的減少混合電場的最大值.考慮到輸電線路實際運行時的電場強度控制原則：交流工頻電場不超過7kV/m，直流電場不超過30kV/m.對于同走廊架壓直流線路和超高壓交流直流線路，建議交流導線最小對地高度大于12m，直流導線最小對地高度大于21m.

圖10 導線對地高度不同時的混合電場分布

5 結束語

本文基于ComsolMultiphysics分別計算了同走廊架設的直流線路和交流線路的電場分布，在此基礎上，使用疊加法計算了二者的混合電場.通過與相關試驗數據的對比，進一步說明了計算的準確性和可靠性；重點討論了交、直流線路的接近距離和導線最小對地高度對混合電場分布的影響.結果表明：合理的接近距離能限制混合區域場強，合理的導線最小對地高能減少混合電場最大值.

[1] 趙永生, 張文亮. 同走廊高壓交、直流輸電線路混合電場分析[J]. 高電壓技術, 2014,40(3): 923-929.

[2] 趙丹, 張彼德, 李明坤，等. 特高壓直流輸電線路合成電場計算方法及其影響因素分析[J].水電能源科學, 2015,33(5):182-185．

[3] 黎小林, 張波, 王琦,等. 特高壓直流線路的場分布影響因素分析[J]. 高壓電器, 2006,42(6): 407-409．

[4] 劉振亞. 特高壓直流輸電工程電磁環境[M]. 北京：中國電力出版社,2009.

[5]張玉寶.基于COMSOLMultiphysics的MEMS建模及應用[M].北京：冶金工業出版社,2007.

[6]ROGERW.Pryor.MultiphysicsModeling:UsingCOMSOL5andMATLAB(Engineering).ElectromagneticForcesonParallelCurrent-CarryingWires[M/CD].MercuryLearning&Information:Har/DVDedition,2015.

[7] 周愷．高壓直流輸電線路電磁環境的計算研究[D].武漢: 華中科技大學, 2007.

[8] 李海峰, 吳冀川, 劉建波，等. 有限元網格剖分與網格質量判定指標[J]. 中國機械工程, 2012, 23(3): 368-377.

[9] 國家環境保護總局.HJ/T24-1998 500kV超高壓送變電工程電磁輻射環境影響評價技術規范[S].

[10] 彭湃, 程漢湘, 陳杏燦,等. 基于有限元的超高壓輸電線路電場的數值分析與測量[J]. 電測與儀表, 2016, 53(3): 56-61.

[11] 吳桂芳, 袁春峰, 陸家榆,等.特高壓直流線路與交流線路同走廊時混合電磁環境的計算[J].電網技術,2010, 34(2): 14-19．

[12]楊勇, 陸家榆, 鞠勇. 交流線路與±800kV直流線路同走廊時的地面混合電場研究[J]. 電網技術, 2009, 33(15): 54-59．

[13] 魯成棟, 肖登明, 秦松林. 交直流并行輸電線路地面混合電場的計算及試驗研究[J]. 高壓電器, 2014, 50(1): 81-86，91.

[14] 陸國慶, 何宏明, 張軍. 交直流輸電線路相鄰架設或共用走廊的探討[J]. 高電壓技術, 1997, 23(4): 68-70.

(編輯：姚佳良)

Calculation of hybrid electric field of AC/DC transmission lines erected in a common corridor based on the finite element method

WANG Jian-lin, JIN Bin, LI Bo-wen, LI Xing, XIONG Ming

(School of Electrical and Electronic Information，Xihua University，Chengdu 610039，China)

The cost of land will be effectively reduced when UHVDC and EHVAC transmission lines are erected in a common corridor, but it brings a new calculation problem of the electric field. In order to study the hybrid electric field, Comsol Multiphysics, a kind of finite element analysis software was used to calculate total electric field below ±800kV UHVDC transmission lines and power frequency electric field below EHVAC transmission lines. Then the hybrid electric field was calculated,and the influence of lines layout on hybrid electric field is discussed lastly. The results show that the hybrid electric field calculation based on Comsol Multiphysics is reasonable and reliable.The reasonable approaching distance and minimum clearance to ground can effectively reduce field strength of mixed region and maximum value of hybrid electric field, respectively.

UHVDC transmission lines; Comsol;common corridor; hybrid electric field; finite element method

2016-07-28

王建林，男，1157304821@qq.com;通信作者：靳斌，男，jb123456@163.com

1672-6197(2017)03-0046-05

TM 711

A