高壓直流輸電系統(tǒng)接地極周圍地表電位分析

李 強，葉紅楓，山江川，劉 念，周運鴻

（1.中核核電運行管理有限公司，浙江 海鹽 314300；2.四川大學(xué)電氣信息學(xué)院，四川 成都 610065）

高壓直流輸電系統(tǒng)接地極周圍地表電位分析

李 強1，葉紅楓1，山江川1，劉 念2，周運鴻2

（1.中核核電運行管理有限公司，浙江 海鹽 314300；2.四川大學(xué)電氣信息學(xué)院，四川 成都 610065）

高壓直流輸電系統(tǒng)大地回路運行時地表電位分布的研究，是評估地下金屬管道腐蝕、變壓器直流偏磁、通信干擾等不良效應(yīng)受地電流影響的重要依據(jù)。利用鏡像法和電磁波的折射、反射規(guī)律，推導(dǎo)出多層垂直和水平復(fù)合分層的地表電位的計算公式，并通過對比不同模型的計算結(jié)果，建立華東地區(qū)的3層水平帶海洋大地模型。最后，通過不同模型仿真計算華東地區(qū)高壓直流輸電系統(tǒng)大地回路單極運行時接地極地表電位的分布情況，得到華東地區(qū)理想的仿真計算模型。

高壓直流輸電系統(tǒng)；地表電位模型；接地極模型；地表電位分布

0 引 言

直流輸電系統(tǒng)接地極周圍地表電位不僅受到高壓直流輸電系統(tǒng)的影響（包括系統(tǒng)的運行方式、輸送功率、傳輸電壓等），而且受接地極的影響（包括接地極的尺寸、形狀、埋設(shè)方式等），其地表電位高低還取決于大地的參數(shù)，包括電阻率的空間分布以及能夠引起空間電阻率改變的濕度、溫度、含鹽比例等[1-2]。特別是由復(fù)雜的地形地貌帶來多變的大地參數(shù)，是地表電位研究的難點。

單極大地回路方式、功率或電壓不相等的雙極不對稱或不平衡方式和同極性方式都有電流通過接地極流向大地，也就是高壓直流輸電系統(tǒng)大地回路運行方式[3]。強大的電流使大地不再為理論上的零電勢，地表出現(xiàn)電位分布。地表不為零的電勢將造成很多負面影響，如金屬管道腐蝕、變壓器直流偏磁、無線電干擾等。地表電位的量化對評估以上負面影響的嚴重程度非常重要，受限于人力、物力，實地測量的方法不能實施，所以模型仿真計算的方法受到廣泛的重視。文獻[4-6]采用單一水平分層、2層水平分層或是采用2層垂直分層的大地模型，具有算法簡單、計算快等優(yōu)點，但沒有對這3種模型的精度加以比較。

本文采用鏡像法和電磁波的折射反射規(guī)律，推導(dǎo)出多層垂直和水平復(fù)合分層的地表電位的計算公式,并通過對比不同模型的計算結(jié)果，分析建立華東地區(qū)的3層水平帶海洋大地模型。最后，通過不同模型仿真計算華東地區(qū)高壓直流輸電系統(tǒng)大地回路單極運行時地表電位的分布情況。

1 與接地極相關(guān)聯(lián)的大地結(jié)構(gòu)

直流接地極的作用：1）作為工作接地，長時間為系統(tǒng)輸送電力，提高系統(tǒng)運行的可靠性；2）作為保護接地，鉗制換流站中性點電位[7]。直流接地極的形狀、尺寸、埋設(shè)的方式等因素同樣對地表電位的計算產(chǎn)生影響。文獻[8-9]給出了接地極電抗和接地極上下表面電流密度不同不影響電位計算的證明，將接地極用線電源模擬。

由于大地的電阻率受許多因素影響，如地理條件和氣候環(huán)境等，所以量化電阻率時要考慮許多變量。大地是高壓直流輸電系統(tǒng)中重要的組成部分，因為系統(tǒng)中的單極電流或者同極的不平衡和不對稱電流都是通過接地極流入大地。地震波對大地進行了粗略的水平劃分，如圖1所示，它必須與具體研究的地區(qū)相結(jié)合，例如華東地區(qū)必須考慮海洋的垂直分層地況。

圖1 大地的結(jié)構(gòu)

地表電位的計算精度隨著劃分空間電阻率的精細程度而改變，除了上述水平劃分和垂直劃分，常見的還有斜面和柱面劃分。當然，用面對大地分層適用于計算遠距離電位。對近距離而言，電流使接地極附近的土壤發(fā)生物理或化學(xué)變化，從而使電阻率發(fā)生變化，變化規(guī)律必須通過函數(shù)關(guān)系式或其他方法來表征，電阻率不能看作是常數(shù)[10]。

不論采用遠距離還是近距離大地模型，分析大地的電場，常用的方法還是求解拉普拉斯方程，輔以狄里赫利邊值條件、紐曼邊值條件或是混合邊值條件。恒定電場中，由電流連續(xù)性定理得到拉普拉斯方程：

式中：γ——電導(dǎo)率標量，S/m；

φ——電位標量，V。

求解拉普拉斯方程可以用解析方法中的鏡像法，也可以用數(shù)值方法中的有限元法或邊界元法等[11-13]。河流、山川對地表電位有影響，但影響不大，但海洋對地表電位的影響較大，解決華東地區(qū)地表電位的分布，就得考慮海洋對地表電位的影響[14]。

2 求解不同分層土壤模型

2.1 無限大均勻大地模型

根據(jù)電流連續(xù)性定理，得到無限大均勻大地模型中點電源周圍的電位分布，電位φ的計算式：

式中：x0、y0、z0——點電源的坐標；

x、y、z——任意點的坐標；

ρ——介質(zhì)的電阻率，Ω·cm；

i——點電源的電流，A。

下面將采用鏡像法和電磁波的折射與反射規(guī)律，推導(dǎo)多層土壤電位分布計算的解析公式，其中ρm（m=0，1，2，3，4，5，6）是相應(yīng)土壤的電阻率，Γm（m=0，1，2，3，4，5）和 Tm（m=1，2）分別是分界面上的反射率和折射率，Qm（m=1，2）和 Pm（m=1，2，3）分別是垂直和水平分界面。

2.2 垂直分層大地模型

圖2 3層垂直大地模型

圖3 2層水平大地模型

如圖2所示，這是一個用2個垂直平面將大地分成3層的大地模型。電流波i在傳播的過程中遇到電阻率不同的土壤介質(zhì)分界面會發(fā)生反射，如果認為在反射波的延長線上出現(xiàn)了一個鏡像電流源，這個鏡像電流源的大小等于反射電流波的大小，那么這個鏡像電源就等效為分界面對電流波的影響，3層垂直不均勻大地模型轉(zhuǎn)化到了單層均勻大地模型中求解。簡化后的單層均勻模型又可以做一次關(guān)于地表的鏡像，最后只需在無窮大均勻大地模型中求解。當然，所求電位的點必須限定在Q2和Q3之間的開區(qū)間，如果要求解左右兩側(cè)土壤的電位分布，需要重新建立鏡像電流源系統(tǒng)。

圖4 3層水平大地模型

只發(fā)生一次反射得到的電位計算公式為

根據(jù)式（3）推出無限次反射后的電位計算解析式為

式（3）和式（4）是針對垂直分界面的反射情況，如果還要考慮水平地面的反射情況，只需要將計算出來的φV替換掉等式（5）中的φ即可。

2.3 2層水平大地模型

如圖3所示，與垂直分層大地模型類似，得到P1和P2之間的開區(qū)間的點電位計算公式（5）。

2層水平分層與3層垂直分層的區(qū)別在于電流波在水平分層模型中傳播時遇到地表會發(fā)生全反射，反射率Γ3=1，也就是沒有折射波。2層水平分層大地模型與3層垂直大地模型一樣沒有折射波。

2.4 3層水平大地模型

大地被平面分為多層水平或垂直層，折射波對電位的影響必須考慮。這里以3層水平分層大地模型轉(zhuǎn)化到無限大均勻大地模型為例，解釋多層模型的求解最終也是轉(zhuǎn)化到求解無限大均勻大地模型。

如圖4所示，電流波在到達P2平面時會發(fā)生折射，折射波在P3平面上發(fā)生反射后折射到第1層土壤中，折射進來的電流波在P1平面全反射后又可以開始第2個周期的傳播。這里忽略了電流波在第1層和第 2 層土壤中經(jīng)過 m 次（m=1，2，3，…）反射，再折射到第1層土壤中的情況。

圖4中注明了第1次和第2次折射波的等效鏡像電流源的大小和坐標。鏡像電流源的大小是通過反射和折射的次數(shù)決定的，如第1個鏡像電流源大小為T1Γ5T2i；坐標是通過幾何圖形推導(dǎo)的，如第1個鏡像電流源z軸的坐標如下：

式中，α1、α2是入射角和反射角，ε1、ε2是土壤的介電常數(shù)，不難發(fā)現(xiàn)鏡像的坐標只與土壤的性質(zhì)、接地極的埋設(shè)深度和大地每一層的深度相關(guān)，而與電流波的大小、入射角和折射角無關(guān)。其他鏡像電流源的大小和坐標可以按上述方法求得，這里就沒有在圖上標識了。

P3平面的等效鏡像電流源取代了折射波對第1層土壤中電位計算的影響，3層水平分層模型轉(zhuǎn)化到2層水平分層模型，最后轉(zhuǎn)化到無限大均勻大地模型中求解。在計算電位的過程中，3層水平分層只是比2層水平分層多了2n（n是折射的次數(shù)）個鏡像電流源，計算方法一樣。

表1 華東地區(qū)地質(zhì)結(jié)構(gòu)

表2 仿真模型參數(shù)

圖5 不同模型仿真計算接地極表面由位分布

3 華東地區(qū)地表電位的仿真計算

3.1 仿真模型參數(shù)

結(jié)合圖1和華東地區(qū)土壤的實際特性得到華東地區(qū)的地質(zhì)結(jié)構(gòu)，如表1所示。選用表2的4種仿真模型，采用Matlab對地表電位計算得到如圖5所示的地表電位曲線。

3.2 仿真分析

比較圖5中不同模型的接地極地表電位曲線，得到如下結(jié)論：

1）地表電位最大值出現(xiàn)在接地極的表面，3000A的入地電流，接地極表面電位能達到100～200 V，并且電位的大小與所選模型無明顯關(guān)系，說明接地極表面的電位不受大范圍土壤的影響。

2）所有模型計算的地表電位都是在0～20km這一段發(fā)生銳減，說明接地極的選址最好離城市或者是110kV以上的變壓器有20km的距離。

3）比較水平分層的2層、3層土壤模型發(fā)現(xiàn)3層土壤模型的地表電位一直比2層的大，而且在20km左右達到最大差值，說明第3層土壤的反射波對地表電位的影響不能忽視。

4）比較水平分層的3層、4層模型發(fā)現(xiàn)兩條電位曲線幾乎重合，說明第4層土壤的反射波對地表電位的影響可以忽略；并且海洋的存在將近海端的整個區(qū)域的地表電位都拉低了，而且在海洋部分的地表電位逐漸為0。

因此，華東地區(qū)的電位分布可采用3層水平帶海洋的復(fù)合土壤模型仿真計算接地極地表電位，這樣既保證了計算的準確性，又能節(jié)省計算機的內(nèi)存和減少耗時。

4 結(jié)束語

本文利用復(fù)合鏡像法和電磁波的反射、折射規(guī)律推導(dǎo)出多層大地模型中地表電位的計算方法。這種方法不僅適合于分析地表電位，對任何電磁波在不均勻介質(zhì)中傳播要求分析介質(zhì)中的某些電磁量都是實用的。結(jié)合華東地區(qū)3層水平帶海洋的復(fù)合土壤模型的仿真結(jié)果，最后提出兩點建議：

1）位于沿海地區(qū)中性點接地的變壓器或者是金屬管道要特別注意高壓直流輸電系統(tǒng)大地回路運行。由于大海的低阻抗的作用，使靠近海洋的城市出現(xiàn)了低電位。大電位差使得該地區(qū)的直流變大，造成更嚴重的金屬腐蝕或者是變壓器的直流偏磁等不良影響。

2）通過擬合函數(shù)可以找到任何兩個點的電位差，這個電位差可以判斷跨步電壓、接觸電壓是否在安全值得范圍內(nèi)，優(yōu)化高壓直流輸電系統(tǒng)接地極的選址。

[1]郝治國，余洋，張保會，等.高壓直流輸電單極大地方式運行時地表電位分布規(guī)律[J].電力自動化設(shè)備，2009，29（6）：10-14.

[2]曾蓮生.高壓直流輸電陸地接地極設(shè)計——關(guān)于地面電位和跨步電壓分布的計算[J].電力建設(shè)，1994，15（1）：24-27.

[3] 葉廷路.葛南直流輸電系統(tǒng)的運行與控制[J].電網(wǎng)技術(shù)，1994，18（1）：20-25.

[4]陸繼明，肖東，毛承雄，等.直流輸電接地極對地表電位分布的影響[J].高電壓技術(shù)，2006，32（9）：55-58.

[5]潘卓洪，張露，劉虎，等.多層水平土壤地表電位分布的仿真分析[J].高電壓技術(shù)，2012，38（1）：116-123.

[6] 劉曲，李立浧，鄭健超，等.復(fù)合土壤模型下HVDC系統(tǒng)單極大地運行時的電流分布[J].中國電機工程學(xué)報，2007，27（36）：8-13.

[7]王彪，王渝紅，丁理杰，等.高壓直流輸電接地電極及相關(guān)問題綜述[J].電力系統(tǒng)及其自動化學(xué)報，2012，24（1）：66-72.

[8]陳水明，施廣德，趙智天，等.直線形直流輸電接地極電流場分析[J].高電壓技術(shù)，1994，20（3）：8-13.

[9]何金良，曾嶸，高延慶，等.電力系統(tǒng)接地技術(shù)研究進展[J].電力建設(shè)，2004，25（6）：1-7.

[10]VILLAS J E T，PORTELA C M.Calculation of electric field and potential distribution into soil and air media for a ground electrode of a HVDC system[J].IEEE Transactions on Power Delivery，2003，18（3）：867-873.

[11]唐新靈，姜建國，喬樹通，等.基于有限元法的地表電位分析[J].現(xiàn)代電力，2013，30（3）：65-68.

[12]傅晨釗，蘇磊，趙丹丹，等.基于有限元分析的直流偏磁電流計算模型研究[J].華東電力，2012，40（10）：1739-1742.

[13]戴傳友，文習(xí)山，方瑜，等.垂直多層土壤接地電阻的計算[J].高電壓技術(shù)，1996，22（3）：47-50.

[14]劉曲，鄭健超，李立呈，等.直流輸電系統(tǒng)變壓器中性點電流分布的影響因素[J].高電壓技術(shù)，2008，34（4）：643-646.

（編輯：李妮）

Analysis on earth surface potential around the electrode of the high voltage DC power transmission system

LI Qiang1， YE Hongfeng1， SHAN Jiangchuan1， LIU Nian2， ZHOU Yunhong2
（1.CNNC Nuclear Operation Management Co.，Ltd.，Haiyan 314300，China；2.College of Electrial Engineering and Information，Sichuan University，Chengdu 610065，China）

The study on distribution of earth surface potential of the high voltage DC power transmission system when ground return is operating is an important evidence of evaluating the underground metal pipeline corrosion， transformer DC bias， communication interference and other adverse effects.The images method and the law of refraction of electromagnetic wave were used to deduce the computing method of earth surface potentials of horizontal and vertical planes.By comparing the calculation results of different models，the model of three layers of soil with a layer of ocean was presented for the East China area.Finally，the distribution of ground return of high voltage DC power transmission system in East China was calculated by different models，and the ideal simulation model for the East China area was obtained.

high voltage DC power transmission system； earth surface potential model； electrode model；distribution of earth surface potential

A

1674-5124（2017）06-0140-05

10.11857/j.issn.1674-5124.2017.06.029

2016-11-20；

2016-12-08

李 強（1983-），男，遼寧新民市人，工程師，主要從事電氣設(shè)備管理檢修。