劉思賢

（國網(wǎng)山東省電力公司濟寧供電公司，山東 濟寧 272100）

0 引言

近年來，竊電現(xiàn)象高發(fā)，一直困擾著供電企業(yè)，嚴(yán)重影響著經(jīng)營效益。一方面竊電現(xiàn)象擾亂了正常的供用電秩序，另一方面竊電也給電網(wǎng)安全帶來嚴(yán)重挑戰(zhàn)，尤其是高科技、隱蔽化、大額化的竊電手段給國家造成了嚴(yán)重的經(jīng)濟損失。從竊電數(shù)據(jù)上看，高壓專用變壓器用戶竊電造成的損失占總竊電量的48.3%［1］。因此，研究高壓專用變壓器用戶竊電識別和損失評估具有一定的經(jīng)濟價值和社會價值。

竊電的手段多種多樣，但無外乎改變電能計量的三要素：電壓、電流及兩者間的相位關(guān)系，改變?nèi)刂械娜我庖粋€，均可使電能表反轉(zhuǎn)、停轉(zhuǎn)或者慢轉(zhuǎn)。竊電的基本方式有欠壓、失流、移相、擴大誤差以及無表、動表等高科技竊電方法［2］。針對用電側(cè)數(shù)據(jù)異常，有學(xué)者依托用電信息采集系統(tǒng)、一體化電量和線損管理系統(tǒng)開展了竊電監(jiān)測和分析工作［3-5］，對欠壓法和失流法竊電取得了較好的識別效果。但單純針對用電數(shù)據(jù)開展竊電監(jiān)測工作存在局限性： 一是無法對全失壓法竊電與無表法竊電進行監(jiān)測；二是對擴大誤差法等竊電方式下的損失電量追補計算缺乏依據(jù)，單純的更正系數(shù)法因電能表殘壓、三相不對稱負(fù)載造成較大誤差。

文獻(xiàn)［6-8］推導(dǎo)出竊電發(fā)生后，配電網(wǎng)絡(luò)電壓、電流、支路功率和竊電負(fù)荷之間的數(shù)據(jù)關(guān)系，提出了基于支路電流、支路電阻以及支路功率損耗作為判據(jù)，從配電網(wǎng)絡(luò)分析的角度實現(xiàn)竊電的識別，給出了損失電量的評估方法。針對低壓配電網(wǎng)絡(luò)，利用智能電表采集電壓數(shù)據(jù)和低壓拓?fù)浣Y(jié)構(gòu)，逐級迭代，實現(xiàn)中間節(jié)點電壓計算值比較分析，若小于0.3 V，則對應(yīng)下游節(jié)點存在竊電嫌疑［9］，該算法在損失電量的計算上，仍然依賴于臺區(qū)線路損耗數(shù)據(jù)。Berrisford 基于智能電能表電流和電壓數(shù)據(jù)，引入低壓配電網(wǎng)絡(luò)狀態(tài)估計的方法實現(xiàn)了竊電識別［10］。文獻(xiàn)［11］建立了支路電阻隨溫度變化的仿真模型，實現(xiàn)了非技術(shù)損耗的計算和估計，提高了竊電偵測成功率。文獻(xiàn)［12-16］基于用電大數(shù)據(jù)分析的方法，實現(xiàn)了竊電的監(jiān)測。

現(xiàn)有竊電監(jiān)測方法均沒有突破營配業(yè)務(wù)邊界，方法缺乏普遍適用性，不能對竊電產(chǎn)生的損失做出有效評估。因此，提出基于營配融合數(shù)據(jù)分析的高壓專用變壓器用戶竊電識別判據(jù)，從配電網(wǎng)絡(luò)的角度建立實際應(yīng)用中更為常見的支路竊電分析模型。在實現(xiàn)竊電識別的基礎(chǔ)上，基于二分迭代法，實現(xiàn)了損失負(fù)荷的精準(zhǔn)評估；針對量測數(shù)據(jù)不充分的配電網(wǎng)絡(luò)，結(jié)合了智能電能表量測數(shù)據(jù)，實現(xiàn)了損失電量的估算。

1 竊電判據(jù)

竊電會導(dǎo)致網(wǎng)絡(luò)運行方式變化，反映在局部區(qū)域或者局部元件上為自導(dǎo)納、互導(dǎo)納的變化。改變一條支路的狀態(tài)或者參數(shù)僅僅影響該支路兩端節(jié)點的自導(dǎo)納和不同節(jié)點之間互導(dǎo)納，因此，發(fā)生一條支路竊電后不必重新形成導(dǎo)納矩陣，只需要對原有導(dǎo)納矩陣相應(yīng)修改。

在實際網(wǎng)絡(luò)中，難以單獨增加一個節(jié)點實施竊電，通常增加一條支路和一個對地節(jié)點實施竊電，所增加支路和對地節(jié)點在網(wǎng)絡(luò)量測系統(tǒng)中是未知的。如網(wǎng)絡(luò)中一條農(nóng)用排灌分支線路實施竊電，量測系統(tǒng)遠(yuǎn)程采集表碼無法正常獲取，造成長期未進行抄表算費，這種情況下，不僅會影響10 kV 同期線損的計算，長期不抄表還會給供電企業(yè)帶來經(jīng)濟損失。

1.1 支路竊電時的節(jié)點電壓和支路損耗網(wǎng)絡(luò)方程

在系統(tǒng)斷面潮流已知的前提下，可將配電網(wǎng)絡(luò)等效為電路網(wǎng)絡(luò)［7］。假設(shè)一個配電網(wǎng)絡(luò)具有n 個節(jié)點，q 個電源，電源為已知量，以電流源描述，記為竊電發(fā)生后，電力系統(tǒng)發(fā)電機出力變化小，可近似為出力不變，即發(fā)電機的等效電流源不變。網(wǎng)絡(luò)如圖1 所示。

圖1 網(wǎng)絡(luò)中增加一條非法用電支路

而圖1 中，從節(jié)點k 引出一條支路，增加不可測量節(jié)點k′，變化后的節(jié)點導(dǎo)納矩陣增加一階。新增加支路與原網(wǎng)絡(luò)節(jié)點k 相連接。原節(jié)點導(dǎo)納矩陣的元素Ykk改變，新增節(jié)點k′只通過支路Lkk′與原網(wǎng)絡(luò)中節(jié)點k 相連，而與其他節(jié)點不直接相連。因此，新節(jié)點導(dǎo)納矩陣中第n+1 列與第n+1 行非對角元素中除Yk（n+1）和Y（n+1）k外，其余均為0。變化后的節(jié)點電壓方程如式（1）所示。

式（1）可寫為

式中：Δe為竊電發(fā)生前網(wǎng)絡(luò)節(jié)點導(dǎo)納矩陣對應(yīng)行列式的值；]。由式（4）—式（6）可知，在網(wǎng)絡(luò)中其他參數(shù)不變的前提下，支路竊電發(fā)生后，節(jié)點k和k′的電壓與竊電支路導(dǎo)納和竊電負(fù)荷有關(guān)。

1）節(jié)點電壓變化。

由式（4）可知，對于網(wǎng)絡(luò)中任一節(jié)點，由于竊電支路和竊電節(jié)點的存在，導(dǎo)致了節(jié)點電壓實際量測值小于竊電發(fā)生前量測值，且大于同一時刻由下游節(jié)點計算電壓值。竊電發(fā)生后節(jié)點i 電壓為

2）支路電流變化。

3）支路功率變化。

4）支路功率損耗變化。

1.2 支路竊電時中間節(jié)點電壓的變化

當(dāng)網(wǎng)絡(luò)中無竊電發(fā)生時，網(wǎng)絡(luò)中各節(jié)點電壓就保持一種穩(wěn)定狀態(tài)，由于電壓數(shù)據(jù)可通過式（11）快速計算，各節(jié)點至末端計量點線路阻抗可通過節(jié)點電壓方程進行求解。通過監(jiān)測節(jié)點電壓變化，實現(xiàn)竊電點識別。

式中：ΔU 為電壓降落的縱分量；δU 為電壓降落橫分量。

取網(wǎng)絡(luò)中任一支路i—j 為例，說明中間節(jié)點電壓迭代求解過程。

1）假定0 號節(jié)點的初始電壓為額定電壓Ubus。

2）由末端向首端推算各個支路損耗以及功率分布。已知節(jié)點j 注入功率為Ploadj+jQloadj，初始電壓為Ubus_j，支路i—j 阻抗為Ri+jXi。

3）支路損耗為：

考慮節(jié)點i 的負(fù)荷PLi+jQLi，注入節(jié)點的功率為：

由于末端注入功率已知，可利用上述方法逐級計算首端節(jié)點的注入功率。

4）假定Ubus_i=Ubus_j，則支路電壓降落分量為：

由支路首端電壓求得支路末端電壓為

在配電網(wǎng)絡(luò)中，中間節(jié)點多為饋線終端（Feeder Terminal Unit，F(xiàn)TU）、站所終端（Distribution Terminal Unit，DTU）設(shè)備，可實現(xiàn)節(jié)點電壓監(jiān)測。由末端計量點電壓、電流、有功功率與無功功率可求得末端節(jié)點到中間節(jié)點的支路阻抗，根據(jù)支路阻抗與末端計量點電壓求出中間節(jié)點的電壓。將中間節(jié)點電壓計算值與FTU 設(shè)備監(jiān)測值比較，若電壓計算值小于監(jiān)測值，則末端計量點或其附近節(jié)點存在異常，從而導(dǎo)致末端節(jié)點電壓降低。圖2 給出了由中間節(jié)點變化實現(xiàn)的竊電識別流程。

圖2 由中間節(jié)點電壓變化實現(xiàn)竊電識別流程

2 IEEE 33 節(jié)點模型的竊電支路識別

在仿真中，網(wǎng)絡(luò)開環(huán)運行，首端基準(zhǔn)電壓12.66 kV。模擬支路竊電如圖3 所示，圖3 中紅色節(jié)點和支路為竊電支路。設(shè)定竊電負(fù)荷為60 kW，阻抗為0.341 0+j0.536 2 Ω。

圖3 改進的IEEE 33 節(jié)點配電網(wǎng)絡(luò)示意

假設(shè)配電網(wǎng)絡(luò)中所有量測裝置可以實時采集、準(zhǔn)確計量，每條支路首端和末端均裝有量測裝置，所有支路的統(tǒng)計損耗可準(zhǔn)確計算，支路統(tǒng)計損耗為輸入電量減去輸出電量。若待計算支路和其他線路有聯(lián)絡(luò)開關(guān)聯(lián)絡(luò)，須對輸入電量進行修正，此時聯(lián)絡(luò)線路的輸入電量為總輸入電量減去聯(lián)絡(luò)線路輸入電量，而聯(lián)絡(luò)線路輸入電量為聯(lián)絡(luò)線路所有計量點電量之和。仿真結(jié)果如表1 所示。

表1 中支路損耗的增加率是計算的每條支路有功功率增加率，更好地模擬了實際情況。統(tǒng)計損耗增加率最大的支路附近存在竊電的可能性最大。并且，距離該支路越遠(yuǎn)，有功功率損耗的增加越不明顯。支路31—32 統(tǒng)計損耗的增加率為454 848.48%，更直觀地反映了網(wǎng)絡(luò)中的異常情況。

將節(jié)點31 作為中間節(jié)點，給出正常網(wǎng)絡(luò)中間節(jié)點在MATLAB 仿真環(huán)境下的計算結(jié)果，如表2 所示。

仿真計算的結(jié)果表明，竊電的發(fā)生會使節(jié)點32電壓降低，與上文理論推導(dǎo)結(jié)果一致。通過節(jié)點32反推計算節(jié)點31 的電壓值和量測值相比有所減小，與上文分析結(jié)果一致。

表1 支路發(fā)生竊電前后支路損耗對比

表2 支路竊電前后中間節(jié)點電壓變化

在實際校驗中，通過同期線損監(jiān)測數(shù)據(jù)分析，出現(xiàn)線損異常之后，逐一計算中間節(jié)點電壓量值。通過計算，發(fā)現(xiàn)一些農(nóng)業(yè)排灌計量點運行狀態(tài)與系統(tǒng)中不一致，實際是運行狀態(tài)，而系統(tǒng)中對應(yīng)停運狀態(tài)，導(dǎo)致了配電設(shè)備監(jiān)測到的電壓比計算出的節(jié)點電壓值偏大，完成運行狀態(tài)整改之后，配電線路同期日線損降至正常范圍。

3 支路竊電損失電量追補計算

3.1 二分迭代法計算損失電量

找到竊電支路后，采用以下方法估算竊電負(fù)荷的大小。

1）根據(jù)配電網(wǎng)絡(luò)原始參數(shù)，應(yīng)用潮流計算程序進行計算，計算出電力網(wǎng)絡(luò)各線路流經(jīng)的復(fù)功率值和理論線損；

2）利用配電自動化設(shè)備得到電網(wǎng)線路功率測量值，計算各個支路統(tǒng)計線損值；

3）比較這些支路統(tǒng)計線損值，如果比理論值明顯偏大，判斷該分支線路存在竊電或者其他異常行為；

4）假設(shè)一點k′為計算分析得到的竊電懷疑點，給定一個竊電負(fù)荷初始值，這個負(fù)荷初始值往往由統(tǒng)計線損的增加量估算得到。由該值歸算到上一個分支線路，計算該分支線路的功率及功率損耗；

5）采用二分法逐步計算，把竊電負(fù)荷加到統(tǒng)計線損中，進行迭代計算，直到修正后的統(tǒng)計線損和網(wǎng)絡(luò)初始線損統(tǒng)計值小于設(shè)定誤差上限，誤差一般取0.000 1。

圖3 中節(jié)點32 出現(xiàn)一條支路竊電，竊電負(fù)荷的大小為60 kW，竊電支路阻抗為0.341 0+j0.536 2 Ω。下面用二分法迭代計算損失負(fù)荷，結(jié)果如表3 所示。

迭代第19 次時，支路32—33 有功功率統(tǒng)計損耗值和實測值相同，竊電的有功功率為60.013 3 kW。迭代到第22 次，支路32—33 無功功率統(tǒng)計損耗值和實測值相同，竊電無功功率為40.020 7 kvar，與已有數(shù)據(jù)一致，驗證二分法進行竊電負(fù)荷估計的可行性。

3.2 實用化的竊電支路損耗計算

在實際的配電網(wǎng)絡(luò)中，支路首端開關(guān)設(shè)備常不具備功率測量精度條件。因此，提出一種營配數(shù)據(jù)融合下的損失負(fù)荷評估方法，可依靠網(wǎng)絡(luò)中的分界負(fù)荷開關(guān)記錄的線電壓和線電流值實現(xiàn)損失電量的估計。

以高壓計量三相三線計量用戶兩相電壓全失壓為例進行計算。計量點A、C 兩相線電流可由分界負(fù)荷開關(guān)電流互感器測得，每5 min 可獲得1 個數(shù)據(jù)。線電壓值僅能獲取A、B 兩相線電壓標(biāo)量值Uab或者C、B 兩相線電壓標(biāo)量值Ucb，用已獲取線電壓近似估計另一線電壓標(biāo)量值。選取可獲取A、C 兩相電流的用戶負(fù)荷開關(guān)，獲取線電壓Ucb共288 個采樣點數(shù)據(jù)。由于智能電表只可召測1 個月內(nèi)共96 個采樣點的電壓、電流和功率因數(shù)數(shù)據(jù)，因此，通過負(fù)荷開關(guān)獲取96 個采樣點數(shù)據(jù)。

式中：P 為電能表計量有功功率；Q 為電能表計量無功功率；Ia為A 相電流二次值；Ic為C 相電流二次值；φ 為功率因數(shù)角。

計算前，須考慮分界開關(guān)到計量柜電壓互感器處線路壓降。取電流曲線相似的工作日96 個采樣點的電壓降，近似替代待求工作日電壓降。式（19）中cos（30°+φ）、cos（30°-φ）估算方法和線電壓類似，查找和目標(biāo)日A、C 相電流曲線最為相似的工作日，取當(dāng)天的96 個采樣點cos（30°+φ）、cos（30°-φ）數(shù)據(jù)，計算有功電量和無功電量。估計結(jié)果如表4所示。

表4 評估結(jié)果和實際值對比

由表4 可知，有功電量估計值與實際值誤差為1.12%，和實測值接近，將損失電量代入到整條線路統(tǒng)計線損中進行校核，與正常線損值一致。說明本算法在計算負(fù)荷較為平穩(wěn)用戶的損失電量有效、可行。

4 結(jié)語

從配電網(wǎng)絡(luò)分析的角度，推導(dǎo)出線路統(tǒng)計損耗以及中間節(jié)點電壓變化2 個竊電判據(jù)，基于IEEE 33 節(jié)點開環(huán)配電網(wǎng)絡(luò)進行了仿真，仿真的結(jié)果表明兩種竊電判據(jù)均為有效。在完成竊電識別的基礎(chǔ)上，給出了竊電負(fù)荷二分迭代計算方法，在量測條件不充分條件下，給出了一種基于配電網(wǎng)絡(luò)量測數(shù)據(jù)的損失電量估算方法，實際應(yīng)用驗證了該方法的有效性。研究成果可為供電企業(yè)反竊電工作提供借鑒，具有推廣應(yīng)用價值。

由于實際的配電網(wǎng)絡(luò)量測條件甚至還滿足不了同期線損計算的要求，一些分支線路、農(nóng)村低壓配電網(wǎng)絡(luò)管理較為薄弱，還沒有實現(xiàn)監(jiān)測數(shù)據(jù)全覆蓋，今后需要在網(wǎng)絡(luò)分析的基礎(chǔ)上，開展不充分量測條件下的負(fù)荷變化較為頻繁用戶竊電識別和損失評估，滿足當(dāng)前反竊電管理的需要。