范開明 王成江 李 光 李紅艷 李如鋒

（1.三峽大學電氣與新能源學院，湖北宜昌 443002；2.臨沂市供電公司，山東臨沂 276000）

變電站的接地網(wǎng)是確保電力系統(tǒng)穩(wěn)定，設備和運行人員安全的一個重要環(huán)節(jié).在我國，普遍使用鍍鋅的碳鋼作為地網(wǎng)的敷設材料.由這種材料連接而成的接地網(wǎng)經(jīng)常會因為腐蝕而降低其電氣性能，成為設備和人身安全事故的隱患，因此接地網(wǎng)故障診斷具有重要意義［1］.

在接地網(wǎng)故障診斷領域已經(jīng)取得了很多成就.文獻［2］從地網(wǎng)拓撲簡約的角度探討接地網(wǎng)故障診斷的可測性問題.文獻［3］探討了減小測量誤差對診斷結(jié)果影響的措施以及怎樣通過優(yōu)化測試方案來提高診斷結(jié)果的精確度.文獻［4］研究了接地網(wǎng)腐蝕診斷初步診斷時的選線原則，提出了大分塊結(jié)合不動點測量的診斷測量方案.文獻［5－7］將非線性規(guī)劃理論應用到接地網(wǎng)故障診斷中來，用迭代的方法求得滿足工程要求的最優(yōu)解.

本文在前人研究的基礎上，運用電網(wǎng)絡理論，在接地網(wǎng)拓撲分層簡約理論的基礎上，建立基于本征接地網(wǎng)層上的故障診斷方程，再結(jié)合物理學中的能量最低原理建立接地網(wǎng)故障診斷數(shù)學模型，成功解決了欠定方程組的求解問題.最后運用最優(yōu)化技術迭代求解，得到了能反應接地網(wǎng)各段導體腐蝕狀況的最優(yōu)解.

1 接地網(wǎng)的拓撲簡約［2］

對于一個具有n個節(jié)點和b條支路的網(wǎng)格形接地網(wǎng)，由于在測量時采用恒定直流電流源進行激勵，故可以忽略其支路的電感和對地電容，將接地網(wǎng)看成是一個無源純電阻網(wǎng)絡.若n個節(jié)點中有m（m＜n）個節(jié)點是連接有接地引上線露出地表可接觸到的，稱其為可及節(jié)點，其余埋于地下的節(jié)點則稱為不可及節(jié)點.其電路模型如圖1所示，其中標記有黑色實心圓點的節(jié)點為可及節(jié)點，其余的為不可及節(jié)點.

圖1 接地網(wǎng)電路模型圖

文獻［2］將一個接地網(wǎng)在拓撲簡約的過程中分別命名為接地網(wǎng)、準元板塊，元板塊、元網(wǎng)絡、可及接地網(wǎng)和本征接地網(wǎng).詳細的拓撲簡約方法和步驟在文獻［2］中有詳細的論述，本文不再贅述.依據(jù)該文獻的結(jié)論可知，本征接地網(wǎng)中的節(jié)點是原始地網(wǎng)中的全部可及節(jié)點，并且任意兩個節(jié)點之間最多只有一條支路，稱其為本征支路.顯然這些支路都可以通過一定的激勵和測量來唯一確定其阻值，然后通過逐層上溯的方法就可以求得接地網(wǎng)層上的支路電阻值或者把故障支路限定在很小的范圍之內(nèi).

2 故障診斷方程的建立

由上節(jié)拓撲簡約分層理論可知，對于具有n個節(jié)點b條支路的本征接地網(wǎng)，因為其中不存在不可及節(jié)點，故各可及節(jié)點相對于公共節(jié)點之間的電壓都是可以實地測量出來的，只要測試儀器的精度滿足要求，再利用一定的措施就可以得到誤差滿足工程要求的可及節(jié)點電壓值，這為下文故障診斷方程的建立提供了可行性保證.

由電網(wǎng)絡理論可知，純電阻網(wǎng)絡的節(jié)點電壓方程的矩陣形式為［8］

其中，A為節(jié)點關聯(lián)矩陣，In為節(jié)點的電流源激勵列向量（非激勵節(jié)點的電流值為零），Yb為支路導納矩陣（為對角方陣，主對角線上的元素是各支路的電阻Rj（j＝1，2，3，…，b），非對角元素為零），Vn為節(jié)點電壓列向量，Yn為節(jié)點導納矩陣.由式（1）可得：

上式左右兩邊對Rj求偏導得：

又由式（2）可知：

將式（3）、（5）代入式（4）可得

將式（6）表示為增量的形式，可表示為

由上式求得的n－1維列向量表示接地網(wǎng)第j條支路的電阻值由于發(fā)生了腐蝕而在設計值的基礎上發(fā)生一個小的變化ΔRj后，其對每個可及節(jié)點電壓值的影響的大小，即ΔVn.如果再用Rj乘以該列向量的每一項，則可以得到一個表示支路電阻變化一倍對節(jié)點電壓影響大小的向量.將式（3）對所有支路電阻變量依次求偏導數(shù)后，將所得到的所有j（j＝1，2，3，…，b）個n－1維列向量合在一起構(gòu)成一個（n－1）×b維的矩陣，記為L.其中的元素Lij反映了在對接地網(wǎng)施加給定激勵的情況下，第j支路的電阻變化對第i個節(jié)點電壓的影響.假設在某種激勵下，接地網(wǎng)n個節(jié)點（除公共節(jié)點）的實際測量值為V′n，理論計算值為Vn，電壓增量ΔVn＝V′n－Vn，則

3 故障診斷數(shù)學模型的建立和求解

由前文拓撲簡約理論可以很容易看出，任何一個層次上的網(wǎng)絡的節(jié)點數(shù)n都少于其相應的支路數(shù)b，故所建立的故障診斷方程組（8）一定是一個欠定的方程組，其方程數(shù)少于未知數(shù)的個數(shù)，由線性方程組解的理論可知，該方程組在大多數(shù)情況下都有無限多個解，而在實際的故障診斷中要求解是唯一的，因此求解該方程必須結(jié)合一定的約束條件或者目標函數(shù)，從而從無限多個解中選出最能反映真實腐蝕情況的最優(yōu)解.

自然界中一個普遍的規(guī)律是“能量越低越穩(wěn)定”，一切自然變化進行的方向都是使能量降低，一個系統(tǒng)總是要調(diào)整自己，使系統(tǒng)的總能量達到最低，使自己處于穩(wěn)定的平衡狀態(tài)，此即能量最低原理（宏觀世界與微觀世界都適用）.因此任何網(wǎng)絡都應滿足這一原理，運用在地網(wǎng)中則是在故障電流通過地網(wǎng)向大地擴散的過程中，所有支路上消耗的能量總和最小.

地網(wǎng)腐蝕后，支路電流表達式為

其中，j＝1，2，…，b，Vj為腐蝕后第j條支路兩端的電壓差，Rj為腐蝕后第j條支路的實際阻值.將式（9）在Rj0點展開成泰勒級數(shù)并忽略非線性項可得：

將式（11）代入式（10）可得：

故障電流通過地網(wǎng)支路流向大地的過程中，各段導體吸收的總能量為

式（13）即為目標函數(shù)，再者，地網(wǎng)腐蝕后，各條支路的電阻都是朝著增大的方向變化，故在建立故障診斷數(shù)學模型的時候還可以加上電阻增量非負的不等式約束條件.至此，接地網(wǎng)故障診斷的數(shù)學模型可以描述如下：

目標函數(shù)： min Q（X）

求得一個盡量靠近真實解的迭代初值X0，這樣可以避免迭代過程收斂于局部最優(yōu)解，然后采用最優(yōu)化理論中的簡約梯度法［9］對其迭代求解，最終獲得一個能反應地網(wǎng)真實腐蝕情況的最優(yōu)解.

4 算例分析

4.1 單支路故障

對圖1所示的地網(wǎng)化簡后所對應的本征接地網(wǎng)進行仿真計算，該層網(wǎng)絡節(jié)點數(shù)n＝10，支路數(shù)b＝21，假設該地網(wǎng)的均壓導體是均勻分布的，各支路電阻的設計值都是0.01Ω，仿真時用0.1Ω的電阻來模擬地網(wǎng)支路的腐蝕.單支路故障設置位置如圖2所示.

圖2 單支路故障示意圖

測量時，選定節(jié)點4為公共節(jié)點，直流電流源可分別加在其余可及節(jié)點與公共節(jié)點之間，但是為了得到比較好的測試效果，本文分別選擇（14，4），（10，4），（8，4）為電流源激勵位置，每次激勵都測量除公共節(jié)點外的所有節(jié)點的電壓值，為了減小測量誤差對診斷結(jié)果的影響，每個節(jié)點電壓值取3次測量的平均值.這樣，一共可以獲得3組測量值，每一組數(shù)據(jù)都可以用來建立故障診斷方程，并求解出能表示故障位置的最優(yōu)解，將診斷結(jié)果中多次被診斷為故障的支路確定為真正的故障支路，否則為偽故障支路.診斷結(jié)果如圖3所示.

圖3 故障診斷結(jié)果

由圖3可以看出，支路（11－16）3次都被診斷為故障支路，而支路（11－15），（11－17）和（14－15）的電阻值都有所增大，但是相對故障支路的增大倍數(shù)來說還是比較小的，且不是每次都有很明顯的增大，所以根據(jù)前文所述的真?zhèn)喂收吓袆e準則，支路（11－16）可以認為是真實故障，其它的都是偽故障，診斷結(jié)果與之前設定的故障位置一致.

4.2 多支路故障

在現(xiàn)實中，變電站的接地網(wǎng)大多數(shù)情況下是有多處腐蝕故障的，因此對多支路故障的仿真診斷是很有意義的.仿真方法、步驟以及故障真?zhèn)蔚呐袆e方法與上節(jié)一樣，只是在設置故障支路位置和故障程度的時候，為了反映故障的多樣性，本文設置了支路（6－16）的故障為腐蝕（支路電阻值增大10倍，即0.1Ω），支路（11－16）的故障為斷裂（支路電阻實際值遠遠大于設計值，本文取增大100倍，即1Ω）.診斷結(jié)果如圖4所示.

圖4 故障診斷結(jié)果

由圖4可以看出，診斷結(jié)果可以定位出本征接地網(wǎng)中故障支路的位置，而且在一定程度上可以反映故障支路的腐蝕程度.這里需要說明的是，由于故障診斷方程是欠定方程，所求的解是滿足目標函數(shù)基礎上的最優(yōu)解，故圖中的縱坐標值并不代表支路電阻的實際增加倍數(shù)，但是能達到故障定位的目的.

在定位出了本征接地網(wǎng)層上的故障支路后，采用逆拓撲簡約理論，即逐層上溯的方法可以確定本征接地網(wǎng)以上各層的支路電阻值［10］.

5 結(jié) 語

現(xiàn)代變電站的接地網(wǎng)占地面積大，節(jié)點數(shù)目一般都在3位數(shù)以上，腐蝕情況很復雜，在初步診斷時先運用電路理論將地網(wǎng)進行拓撲簡約化簡，得到只含全部可及節(jié)點的本征接地網(wǎng)，建立本征接地網(wǎng)層上的故障診斷方程，可以有效地降低故障診斷方程的欠定程度，再結(jié)合一定的優(yōu)化測量措施，大大降低了診斷結(jié)果的誤差.在求解故障診斷方程時，引入物理學中的能量最小原理建立故障診斷數(shù)學模型，運用最優(yōu)化技術迭代求解，得到了能反映故障位置和故障程度的最優(yōu)解.通過對一個小型地網(wǎng)的模擬計算驗證了該方法的正確性和可行性.

［1］ 劉 建，王建新，王 森.一種改進的接地網(wǎng)故障診斷算法及測試方案評價［J］.中國電機工程學報，2005，25（3）：71－77.

［2］ 劉 建，王樹奇，李志忠，等.地網(wǎng)腐蝕故障診斷的可測性研究［J］.高電壓技術，2008，34（1）：64－69.

［3］ 劉 建，王建新，王 森.量誤差對接地網(wǎng)故障診斷影響的分析［J］.高電壓技術，2006，32（4）：95－110.

［4］ 劉渝根，王碩，田金虎，等.接地網(wǎng)腐蝕診斷優(yōu)化測量方法［J］.重慶大學學報，2008，31（11）：1303－1306.

［5］ 張曉玲，黃青陽.電力系統(tǒng)接地網(wǎng)故障診斷［J］.電力系統(tǒng)及其自動化學報，2002，14（1）：48－51.

［6］ 肖新華，劉 華，陳先祿.接地網(wǎng)腐蝕和斷點的診斷理論分析［J］.重慶大學學報，2001，24（3）：72－75.

［7］ 劉渝根，吳立香，王 碩.大中型接地網(wǎng)腐蝕優(yōu)化診斷實用化分析［J］.重慶大學學報，2008，31（4）：417－420.

［8］ 邱關源.電路［M］.北京：高等教育出版社，2000.

［9］ 韋鶴平.最優(yōu)化技術應用［M］.上海：同濟大學出版社，1987.

［10］劉 健，王樹奇，李志忠，等.基于網(wǎng)絡拓撲分層簡約的接地網(wǎng)腐蝕故障診斷［J］.中國電機工程學報，2008，28（16）：122－128.