變電站接地網參數的計算與分析

秦大海 戴玉松 徐 聞 黃滿意 姜聿涵

（1.西華大學電氣信息學院，成都 610039；2.國網成都供電公司，成都 610041）

隨著電力輸送向著長距離、高電壓、大容量的發展，簡單的接地網已經不能滿足電力系統安全性和可靠性的要求，因而發變電站的地網規模不斷擴大，結構越來越復雜。故障或雷擊電流在地網中的泄流也越來越復雜，過電流引起接地網參數的變化也越來越復雜。與此同時，隨著變電站自動化水平的不斷提高，越來越多的微電子裝置應用于變電站的測量、監控、保護等中，而這些微電子裝置的耐壓水平很低、抗干擾的能力較弱[1-9]。綜上所述，接地網的好壞不僅僅危及工作人員的安全，也與設備的正常運行息息相關，因此對接地網的相關參數進行評估具有現實的意義。

接地電阻一直是接地網評估的一個重要指標，因為接地電阻與地面的電位分布息息相關。但是即使接地電阻滿足規程要求，也有可能導致與地網相連的設備受到反擊的影響，同時在地面上產生危險的跨步電壓和接觸電壓。因此，對于接地網的評估，應該綜合考慮接地電阻、地面電位分布、電位梯度分布、接觸電壓、跨步電壓[10-13]。

本文采用基于恒定電流場理論的數值計算方法，首先根據地網結構參數得到接地網的數據文件，然后將數據文件用編寫的程序計算，得到接地網的電阻、電位文件，最后使用數據處理軟件繪制接地網電位分布和梯度分布的三維圖、色階圖和灰度圖。根據繪制的圖形對地網的安全性進行評估，最后提出相應的防護措施。

1 接地網參數計算的原理和方法

1）恒定電流場理論

設一恒定電流I由接地電極注入地網，由恒定電流場理論，使用格林函數的原理，設無窮遠處為零電位點，則電極泄漏電流在任意一點P的電位為

式中，J(Q) 是接地體電極表面S上的Q點處的泄漏電流密度，G(P,Q)是相應電極幾何形狀的格林函數，此處代表單位電流經過電極表面Q點在P電產生的電位值。

經過接地網所有導體在土壤中的泄漏電流等于注入接地電極的電流I：

設電極表面等電位，則邊界條件為

2）表面電荷法

取無窮遠處為零電位點，注入接地極的電流為恒定值I，均勻土壤電阻率為ρ。設電極總長為L，通過電極總的泄流電流為I，將L分成n個微段，第j段長為Lj，其中心為Oj，第j段的泄流電流為Ij，則有

G(P,Oj)是以Oj為中心的微段向地中注入單位 電流時在P點產生的電位。若P點在電極微段i上，每個微段通過泄流均要在i點產生電位，可用互電 阻Rij表示，當i=j時，定義為i段的自電阻Rii，則式（5）可以表示為

為簡化計算，我們忽略金屬導體上的電壓降，即認為電極表面等電位，設電極電壓為VG，則

由式（4）和式（5），得

由式（5）和式（6）組成n+ 1個線性方程組，寫成矩陣形式為

式中，E為單位陣，VG是列矢量，即

則接地體接地電阻為

可見，只要求出了電阻系數，然后根據式（5）就可以求得任一點P的電位，最后由式（11）就可以求出接地電阻。

3）Rij和G(P,Oj)的計算

由前面分析可知，格林函數代表一個微段單元注入土壤一單位電流時，在P點產生的電位。

當P點距離電極微段單元很遠時，即微段電極的尺寸與距離相比可以忽略不計，則電極可以看成一個點模型，則有

圖1中r和r'分別是點P到微段中心點和到微段鏡像的中心點的距離。

圖1 鏡像法計算示意圖

當P點與電極微段很近時，特別是P點在微段表面時，采用點電源公式會產生很大的誤差。這時可以把電極微段進行二次分段處理，即把電極微段按長度等分為N段。設dkj為微分段Kj 流入地中的電流，由j段流入地中電流為單位電流時，則有

當單位電流由微分段Kj注入土壤時，在P點產生的電位為φkj(P)，則有

如圖2，a為導體半徑，b為微分段長度；r和z分別為P點的橫坐標和縱坐標，k0和k1分別是零階與一階的第二類修正的貝塞爾函數，0+表示積分下限不能為0。

圖2 微分段坐標及幾何尺寸示意圖

式（14）是三維空間均勻介子情況的計算，對于半無限均勻介質，計算應該由鏡像法求解。先求出φkj(r',z')，則P點的電位為

當P點與微段很接近時，只要求得各就可得到P點與該微段較近時的格林函數值G(P,Oj)。

設有單位電流由第j微段注入地中，這時第j微段的導體電位升高在數值上等于Rij。把P點分別固定于各微段上，考慮到電極表面等電位，得到如下N+ 1個方程組：

微段i和微段j的互電阻，即在j段上注入單位電流而在i段上產生的平均電位，有

式中，Li為第i微段導體總長度；Pdis為第i段上微分單元ds中心點；G(Pd is,Oj)為單位電流注入j段在Pdis點產生的電位，有Rij=Rji。

為使式（15）計算簡化，當i、j段較遠時，可以將Li、Lj微段作為點電源處理，這時Rij就可用式（9）計算。當i、j段相距很近時，對分布在i段上的三個點P1i、P2i、P3i，用式（14）分別計算G(P1i,Oj)、G(P2i,Oj)和G(P3i,Oj)，取它們的平均值作為Rij的近似值：

根據地面上任一P點的電位Vp，可算出P點的電位梯度值 ?VP，以及P點的接觸電壓UTP，有

式中，Vg為接地設備離地面1.8m 處的對地電位升。

設人腳如果接觸地表點P和距離點P為0.8m的點M兩點，則跨步電壓US為

式中，VM為M點的電位值。

由上述可知，給定了一個接地網的原始數據，就可以計算出地網接地電阻；確定接地電流之后，就可以計算地面上任意一點電位；最后計算出地面電位梯度、跨步電壓和接觸電壓。

2 算例計算與分析

220kV 保國變電站的接地網結構，如圖3所示。圖中黑色小圓點表示垂直接地體的布置位置。該接地網的相關參數如下：接地網的長為174m，寬為101.3m，土壤電阻率為200Ω·m，接地網的埋深為0.8m。該接地網一共設置了三處帽檐式均壓帶，內外弧的埋深分別為：1m、1.5m。接地網使用的材料為圓鋼，垂直接地體的半徑為0.01m，長度為2.5m；水平導體的半徑為0.0061804m。規則導體的數目為101 根，不規則導體的數目為826 根。注入地網的電流為10kA 的恒定電流。為了使計算結果更加精確，本文選取計算步長為 0.25m，計算的點數取 768×660。

圖3 變電站接地網結構圖

根據程序計算的結果，將電位分布數據矩陣應用數據處理軟件（Origin8.0）繪制得到電位分布色階圖，如圖4所示。在該圖中，橫縱坐標代表不同電位點，不同的顏色代表不同的電位值，電位的單位為伏特。由圖可知，電位分布基本按網格狀分布，并且在不同的電位之間通過黑線來隔開。電位分布在地網所在區域的地面上電位較高，離開地網后電位迅速下降，因此在地網邊緣處可能會出現較高的電位差，危及工作人員的人身安全和設備的正常運行。在地網中心部分有些密集的黑色重疊區，說明此處有多個電位等級，而在地網邊緣出現黑色重疊區是因為靠近地網邊緣電位迅速下降的緣故，在這些黑色重疊區電位梯度可能較大，若最大電壓差超出允許值，應采取適當的均壓措施。從該圖還可以看出在變電站左下角設置有均壓帶的地方電位變化比較緩慢，電位梯度較小。

圖4 地網地面電位分布色階圖

圖5 地網地面電位分布灰度圖

圖5所示為通過Origin 畫出的電位分布灰度圖，該圖不同亮度代表電位的高低，亮度越高的地方代表電位越高。由該圖可以看出，接地網所在區域的地面上的亮度除少數地方外亮度相當，即整個地網的電位均壓效果較好。同時，在地網的邊緣處亮度變化很大，對應的電位變化較快。通過觀察各點的亮度，可以得出電位在地面的分布情況，在靠近中心的網格節點處亮度最高，說明此處是整個地網地面電位最高的點，對照圖形的橫縱坐標找到大致位置，具體數值可以通過地網地面電位分布數值矩陣得出。

圖形中亮度差別較大的地方說明此處電壓差較大，根據此圖，再結合電位分布數值矩陣可以計算出地網中跨步電壓較大的位置以及大小。

圖6 地網電位分布三維圖

根據對程序計算結果的電位分布數據矩陣進行處理，得到如圖6所示的電位分布三維圖，其中x、y軸分別為不同電位點的橫坐標與縱坐標，z軸表示地面上任一點的電位值。由圖可知地面電位分布呈現以下特點：導體上的電位高于網格中心的電位，有垂直接地極的地方電位較高，由規則導體向不規則導體的過度區域電位變化很大，地網邊緣處電位落差較大。由圖可知電位的最大值為6250V 左右，平均值在4250 左右。

圖7 地面電位梯度三維圖

由圖7電位梯度三維圖，我們可以看到在地網的中間區域，地面電位梯度都比較小。而在靠近地網邊緣處，電位梯度變化迅速。尤其在X為700 步，Y為600 步左右，電位梯度達到90V/m 左右。因此，根據計算結果，我們應該加強地網邊緣的防護措施。同時，我們還可以清晰的看到，在地網的左下角變電站大門處，由于敷設了帽檐式均壓帶，該處的電位梯度只有25V/m。

圖8 地面電位梯度等高線圖

圖9 地面電位梯度灰度圖

圖8和圖9分別是電位梯度分布的等高線圖和灰度圖，由等高線圖可以清晰的看到相應點的電位梯度，電位梯度值大多分布在20V/m 左右。由灰度圖可以清晰的看到該地網的總體設計效果，該接地網設計效果總體較好，只需要加強地網邊緣的防護，以確保人身和設備的安全。

3 結論

1）根據數值計算的結果，該接地網的接地電阻為0.64683 Ω，基本滿足了地網安全性要求。

2）從電位分布三維圖看，接地網的電位分布在地網所在區域較均勻，只有極少地方有高的電位升，如果在地網中間區域接有對過電壓特敏感的設備，則應采取相應的降壓措施。同時，在地網邊緣處的電位落差較大，跨步電壓和接觸電壓相應較高，為了人身和設備的安全，可在地網邊緣處鋪一層厚度為3～10cm 的礫石或用瀝青混凝土加大地表電阻率就可以使跨步電壓和接觸電壓降低到安全范圍內。

3）從電位梯度圖可知，該地網總體的均壓效果很好，在地網所在區域沒有明顯的突刺，只在地網邊緣有較少的突刺。

4）針對圖中出現的黑色重疊區，我們可以改變導體的布置，使用帽檐式均壓環等措施來降低電位差。

[1] 李雷.接地參數計算與接地網頻率響應特性下的研究[D].長沙: 長沙理工大學碩士論文,2010.

[2] 李欣.接地網地面電位分布及梯度分布的計算研究[D].成都: 西華大學碩士論文,2010.

[3] 何金良,曾嶸.電力系統接地技術[M].科學出版社,2007(2).

[4] 張波,崔翔,趙志斌,等.大型變電站接地網的頻域分析方法[J].中國電機工程學報,2002,22(9): 59-63.

[5] 魯志偉,文習山,史艷鈴,等.大型變電站接地網工頻接地參數的數值計算[J].中國電機工程學報,2003,23(12)9: 89-92.

[6] 解廣潤.電力系統接地技術[M].北京: 科學出版社,2007.

[7] 李景祿,等.實用電力接地技術[M].中國電力出版社,2002(1).

[8] 竇鑫,何金良,曾嶸.變電所接地阻抗的感性分量分析[J].高壓電器,2004,40(4): 245-246.

[9] 魯志偉,常樹生,張久祿.大型變電站接地網的網內電位差[J].電力建設,2004(9): 39-40.

[10] LIN L.Influence of soil and conducter of ground grid on satey of the grounding system.IEEE 2000.

[11] American National Standards Institute.IEEE Guide for Measuring Earth Resistivity,Ground Impedance,and Earth Surface Potentials of a Ground System.1983.

[12] 錢金忠.220kV 接地網改造的技術問題研究[D].南京: 東南大學工程碩士學術論文,2006.

[13] 甘艷,阮江軍,陳允平.短路故障下接地網電位分布的討論[J].華北電力,2003(12): 7-10.