馮華然，庾力維，何敬開，梁錫亮，蔡昭群

(1.廣東電網(wǎng)有限責(zé)任公司東莞供電局，廣東 東莞 523000；2.廣州依萊科電力科技有限公司，廣東 廣州 510630)

隨著電網(wǎng)發(fā)展，配電網(wǎng)規(guī)劃、供電部門日常負(fù)荷管理等工作日益精細(xì)化?？臻g負(fù)荷預(yù)測(cè)[1-2]作為配電網(wǎng)規(guī)劃和負(fù)荷管理的基礎(chǔ)，其預(yù)測(cè)的精細(xì)度、準(zhǔn)確性要求越來越高。目前不少地區(qū)已逐漸普及以網(wǎng)格[3-4]、地塊[5-6]為單位進(jìn)行空間負(fù)荷預(yù)測(cè)。而現(xiàn)狀負(fù)荷統(tǒng)計(jì)是負(fù)荷預(yù)測(cè)的基礎(chǔ)，無論是短期負(fù)荷預(yù)測(cè)、中期負(fù)荷預(yù)測(cè)、遠(yuǎn)期飽和負(fù)荷預(yù)測(cè)均需對(duì)比參照現(xiàn)狀負(fù)荷數(shù)據(jù)。因此，準(zhǔn)確統(tǒng)計(jì)地塊現(xiàn)狀負(fù)荷，對(duì)配電網(wǎng)規(guī)劃和負(fù)荷管理意義重大。

目前供電部門現(xiàn)有各數(shù)據(jù)系統(tǒng)可方便地查詢、統(tǒng)計(jì)省、市、鎮(zhèn)等各級(jí)行政單位負(fù)荷、各變電站負(fù)荷、線路負(fù)荷、變壓器負(fù)荷、電表負(fù)荷，但無記錄、查詢地塊負(fù)荷。目前規(guī)劃管理人員統(tǒng)計(jì)地塊現(xiàn)狀負(fù)荷往往依賴于人工，常用的主要有以下2種估算方法：

a)配用變壓器(以下簡(jiǎn)稱“配變”)負(fù)載率法。通過統(tǒng)計(jì)現(xiàn)狀地塊內(nèi)各行業(yè)配變?nèi)萘?，乘各行業(yè)配變實(shí)際平均負(fù)載率[7]，結(jié)果求和再乘同時(shí)率得出地塊現(xiàn)狀負(fù)荷。該方法需人工逐個(gè)地塊統(tǒng)計(jì)配變?nèi)萘?，?dāng)?shù)貕K數(shù)量多時(shí)，該工作量巨大。該方法默認(rèn)同一行業(yè)的配變采用相同的負(fù)載率，同時(shí)率取值受主觀影響較大。

b)線路負(fù)荷分配法。通過統(tǒng)計(jì)現(xiàn)狀地塊內(nèi)各配變?nèi)萘考捌渌鶎?0 kV線路，按各回路10 kV線路位于該地塊內(nèi)的配變?nèi)萘空荚摼€路總配變?nèi)萘康谋壤?，?0 kV線路負(fù)荷分配到各地塊，從而得出各地塊現(xiàn)狀負(fù)荷。該方法需人工逐個(gè)地塊統(tǒng)計(jì)配變?nèi)萘亢退鶎倬€路，當(dāng)?shù)貕K數(shù)量多時(shí)，該工作量巨大。該方法默認(rèn)同一線路的配變采用相同的負(fù)載率，同時(shí)率取值受主觀影響較大。

傳統(tǒng)統(tǒng)計(jì)方法較多依賴于人工，不適合大規(guī)模統(tǒng)計(jì)。為此，本文提出一種基于地理信息系統(tǒng)(geographic information system，GIS)[8-9]技術(shù)的統(tǒng)計(jì)方法。該方法通過GIS平臺(tái)確定配變與地塊的耦合關(guān)系，無需人工逐個(gè)地塊統(tǒng)計(jì)，適用于大規(guī)模地塊現(xiàn)狀負(fù)荷統(tǒng)計(jì)，可有效提高統(tǒng)計(jì)效率。該方法顆粒度精細(xì)至配變，各配變最大負(fù)荷均由計(jì)量自動(dòng)化系統(tǒng)[10-11]導(dǎo)出，通過負(fù)荷曲線疊加計(jì)算地塊負(fù)荷，擺脫了傳統(tǒng)統(tǒng)計(jì)方法中配變負(fù)載率一刀切和同時(shí)率主觀取值影響，可有效提高統(tǒng)計(jì)結(jié)果可靠性。

1 總體方法流程

本文所提出的基于GIS技術(shù)的大規(guī)模地塊現(xiàn)狀負(fù)荷統(tǒng)計(jì)方法流程如圖1所示。

圖1 總體方法流程Fig.1 Overall method flow

基于GIS技術(shù)的大規(guī)模地塊現(xiàn)狀負(fù)荷統(tǒng)計(jì)方法流程為：

a)根據(jù)政府規(guī)劃部門提供的總規(guī)和控規(guī)資料，利用MapInfo軟件，建立地塊GIS信息數(shù)據(jù)庫。

b)將配變坐標(biāo)轉(zhuǎn)換[12-13]為GIS信息數(shù)據(jù)庫所用坐標(biāo)，利用MapBasic編程對(duì)MapInfo軟件進(jìn)行二次開發(fā)[14-15]，增加配變與地塊耦合功能，實(shí)現(xiàn)配變與地塊的耦合。

c)確定地區(qū)和配變的行業(yè)，提取主要行業(yè)日負(fù)荷曲線[16-17]和配變最大負(fù)荷，進(jìn)而模擬出各配變?nèi)肇?fù)荷曲線。

d)將各地塊內(nèi)配變?nèi)肇?fù)荷曲線進(jìn)行疊加得出各地塊日負(fù)荷曲線，從中提取地塊最大負(fù)荷。

2 建立地塊GIS信息數(shù)據(jù)庫

2.1 MapInfo平臺(tái)

建立地塊GIS信息數(shù)據(jù)首先需要GIS平臺(tái)，如本文采用的MapInfo軟件。MapInfo可以將地塊面積、用地屬性、容積率、地塊配變?nèi)萘俊⒇?fù)荷水平等數(shù)據(jù)在GIS平臺(tái)上記錄、匯總,形成各種地理分布專題圖，方便規(guī)劃人員查看、統(tǒng)計(jì)任意范圍內(nèi)的數(shù)據(jù)，專題圖如圖2所示。

圖2 MapInfo地塊信息統(tǒng)計(jì)Fig.2 Information statistics of land parcel in MapInfo

2.2 建立地塊GIS信息數(shù)據(jù)庫

政府規(guī)劃部門提供的總規(guī)和控規(guī)地塊信息一般為CAD文件或mdb數(shù)據(jù)庫文件，需要先整理原始資料，提取有用數(shù)據(jù)，再轉(zhuǎn)換為MapInfo TAB文件，最后處理TAB文件字段、圖層，導(dǎo)入地塊信息形成地塊GIS數(shù)據(jù)庫，如圖3所示。

3 配變與地塊耦合

3.1 配變坐標(biāo)轉(zhuǎn)換和導(dǎo)入GIS數(shù)據(jù)庫

首先從供電部門營(yíng)銷系統(tǒng)[18]提取各配變坐標(biāo)數(shù)據(jù)，該系統(tǒng)坐標(biāo)數(shù)據(jù)通常為WGS84大地坐標(biāo)，需統(tǒng)一轉(zhuǎn)換為GIS數(shù)據(jù)庫所用坐標(biāo)(一般為平面坐標(biāo))。

圖3 地塊GIS信息數(shù)據(jù)庫Fig.3 GIS information database of land parcel

對(duì)于范圍較廣的地區(qū)坐標(biāo)轉(zhuǎn)換，如某市配變坐標(biāo)轉(zhuǎn)換，可采用七參數(shù)法轉(zhuǎn)換，以降低誤差?？蓮膰?guó)土部門直接獲取轉(zhuǎn)換參數(shù)，或者獲取若干控制點(diǎn)坐標(biāo)自行計(jì)算七參數(shù)。

利用MapInfo的創(chuàng)建點(diǎn)功能，將配變導(dǎo)入MapInfo形成配變點(diǎn)，1個(gè)配變1個(gè)點(diǎn)，如圖4所示。

圖4 配變導(dǎo)入GIS數(shù)據(jù)庫Fig.4 Import of distribution transformer into GIS database

3.2 配變與地塊耦合

GIS數(shù)據(jù)庫中可清楚顯示配變位置、所在地塊，但要查看所在地塊編號(hào)還需要單獨(dú)點(diǎn)開查看地塊信息，這樣顯然效率太低。

通過MapBasic編程對(duì)MapInfo軟件進(jìn)行二次開發(fā)，增加配變與地塊耦合功能，可批量實(shí)現(xiàn)配變與地塊的耦合，精準(zhǔn)大批量確定配變所在地塊。

3.2.1 MapBasic與二次開發(fā)

MapBasic是MapInfo軟件的二次開發(fā)工具，運(yùn)行于MapInfo環(huán)境下的編程語言，語法與Basic相似，可用于擴(kuò)展MapInfo軟件功能，實(shí)現(xiàn)程序的自動(dòng)重復(fù)操作。MapBasic具有功能很強(qiáng)的內(nèi)含地理特性的語句和函數(shù)，利用它可以完成很多有關(guān)圖元地理特性的復(fù)雜運(yùn)算和操作，并且支持多種技術(shù)標(biāo)準(zhǔn)，易與其他應(yīng)用軟件互通、集成，便于集成二次開發(fā)。

目前GIS二次開發(fā)主要有獨(dú)立開發(fā)、單純二次開發(fā)和集成二次開發(fā)3種方式。其中單純二次開發(fā)是指利用GIS平臺(tái)提供的二次開發(fā)工具，直接在該GIS平臺(tái)上進(jìn)行開發(fā)，開發(fā)出的功能、插件直接運(yùn)行于該GIS平臺(tái)。此方法開發(fā)周期短，省時(shí)省力，本文采用單純二次開發(fā)方式實(shí)現(xiàn)配變與地塊的耦合。

3.2.2 配變與地塊耦合功能開發(fā)

地理運(yùn)算符[19]能夠根據(jù)某些對(duì)象與其他對(duì)象的空間關(guān)系來選中。MapInfo有1個(gè)能用于地理運(yùn)算符的特殊列名：“obj”或“object”。此列名引用與表相連接的圖形對(duì)象，地理運(yùn)算符作用于所指定的對(duì)象。可據(jù)此判斷配變所屬地塊，具體實(shí)現(xiàn)步驟如下：

a)判斷打開的圖層是否為2個(gè)(地塊圖層、配變圖層)，如果不是則提示用戶打開這2個(gè)圖層，程序停止往下執(zhí)行。

b)檢測(cè)地塊圖層是否含有線、弧等非區(qū)域圖元，如果是，則提示用戶，程序停止往下執(zhí)行。

c)設(shè)置內(nèi)部坐標(biāo)系，以便程序能夠在MapInfo上準(zhǔn)確運(yùn)行。

d)將配變圖層的“所屬地塊”字段清空。

e)逐個(gè)計(jì)算配變所屬地塊。利用MapBasic的地理運(yùn)算符intersects計(jì)算配變所屬地塊。

f)若配變不在任何地塊中，則建立緩沖區(qū)放寬條件再查找。以配變所在位置為中心建立緩沖區(qū)，緩沖區(qū)半徑從2 m開始，搜索該緩沖區(qū)是否與某地塊相交，如果搜索不到任何地塊則緩沖區(qū)半徑增加2 m，直到半徑為200 m為止。

g)如果配變或其緩沖區(qū)與某一地塊相交，則將該地塊的編碼寫入該配變的“所屬地塊”字段。

4 模擬配變?nèi)肇?fù)荷曲線并統(tǒng)計(jì)地塊負(fù)荷

理論上要統(tǒng)計(jì)地塊現(xiàn)狀最大負(fù)荷，只需將地塊內(nèi)全部配變負(fù)荷曲線疊加取最大值即可,但配變負(fù)荷數(shù)據(jù)來源于計(jì)量自動(dòng)化系統(tǒng)，其采樣周期為15 min，每臺(tái)配變每年采樣52.56萬次。以某市為例，全市共約8萬臺(tái)配變，則合共采樣約420億次，且該市共規(guī)劃5萬多個(gè)地塊。對(duì)于規(guī)劃人員，要從計(jì)量自動(dòng)化系統(tǒng)提取如此巨量的數(shù)據(jù)，并按地塊進(jìn)行疊加分析顯然是不現(xiàn)實(shí)的。

本文提出利用各配變最大負(fù)荷和各行業(yè)日負(fù)荷曲線，模擬各配變?nèi)肇?fù)荷曲線進(jìn)行疊加統(tǒng)計(jì)地塊負(fù)荷。

根據(jù)計(jì)量自動(dòng)化系統(tǒng)的電量和負(fù)荷統(tǒng)計(jì)數(shù)據(jù)，確定地區(qū)主要行業(yè)分類。系統(tǒng)現(xiàn)有行業(yè)分類[20]200多個(gè)，理論上行業(yè)分類越精細(xì)，則統(tǒng)計(jì)結(jié)果越準(zhǔn)確；但考慮到實(shí)際可操作性問題，可以將性質(zhì)相近的行業(yè)合并，最終合并出少量主要的行業(yè)分類。

部分配變可能包含多種類型負(fù)荷，但這部分負(fù)荷沒有明確的數(shù)據(jù)來源，現(xiàn)有技術(shù)條件下無法大規(guī)模精確統(tǒng)計(jì)。按照運(yùn)行經(jīng)驗(yàn)，這部分負(fù)荷很小，對(duì)于大規(guī)模的統(tǒng)計(jì)可忽略不計(jì)。如果是局部區(qū)域少量地塊統(tǒng)計(jì)，可細(xì)化分析低壓負(fù)荷性質(zhì)，使統(tǒng)計(jì)結(jié)果更精確。

從計(jì)量自動(dòng)化系統(tǒng)可導(dǎo)出上文主要行業(yè)分類的日負(fù)荷曲線。該日負(fù)荷曲線既要能夠準(zhǔn)確反映該類負(fù)荷日負(fù)荷特性，又要盡量精簡(jiǎn)數(shù)據(jù)以提高后期數(shù)據(jù)處理效率，日負(fù)荷曲線的采樣周期需要在兩者之間取得平衡。負(fù)荷曲線縱軸統(tǒng)一以各點(diǎn)負(fù)荷相對(duì)于最大負(fù)荷的比值來表示。

從計(jì)量自動(dòng)化系統(tǒng)可提取各配變所屬行業(yè)，按上文分類原則歸類到各主要行業(yè)。從計(jì)量自動(dòng)化系統(tǒng)后臺(tái)可導(dǎo)出各配變最大負(fù)荷。各配變最大負(fù)荷乘所屬主要行業(yè)日負(fù)荷曲線對(duì)應(yīng)數(shù)據(jù)即可模擬得出各配變?nèi)肇?fù)荷曲線。

根據(jù)之前確定的配變與地塊的耦合關(guān)系，把同一地塊內(nèi)的配變?nèi)肇?fù)荷曲線疊加即可得出各地塊的日負(fù)荷曲線，取其中最大值即為地塊現(xiàn)狀最大負(fù)荷值。

5 實(shí)例分析

5.1 實(shí)例計(jì)算

以某市2018年配電網(wǎng)為例，采用本文所提出的方法統(tǒng)計(jì)全市各地塊負(fù)荷現(xiàn)狀。

首先通過MapInfo建立全市地塊GIS信息數(shù)據(jù)庫，把全部配變通過坐標(biāo)匹配導(dǎo)入MapInfo，然后使用二次開發(fā)的配變與地塊耦合功能，確定各配變所屬地塊，如圖5所示。

圖5 確定配變所屬地塊Fig.5 Determination of land parcel of distribution transformer

將該市各行業(yè)劃分為居民、商業(yè)、辦公、電氣機(jī)械及設(shè)備制造業(yè)、電子信息制造業(yè)、紡織服裝鞋帽制造業(yè)、食品飲料及煙草制造業(yè)、造紙及紙制品業(yè)、玩具及文體用品制造業(yè)、木材加工及制品和家具制品業(yè)、化工制品制造業(yè)、包裝印刷業(yè)、金屬制造業(yè)、非金屬制造業(yè)、其他制造業(yè)和其他共16個(gè)主要行業(yè)。將該市全部配變歸類到該16個(gè)主要行業(yè)中。

從計(jì)量自動(dòng)化系統(tǒng)導(dǎo)出該16個(gè)主要行業(yè)日負(fù)荷曲線，采樣周期為1 h，如圖6所示。

通過主要行業(yè)日負(fù)荷曲線和各配變最大負(fù)荷，模擬計(jì)算各配變?nèi)肇?fù)荷曲線，計(jì)算結(jié)果見表1。

利用配變與地塊的耦合關(guān)系，疊加各地塊內(nèi)的配變?nèi)肇?fù)荷曲線，得出各地塊的日負(fù)荷曲線，計(jì)算結(jié)果見表2。各地塊日負(fù)荷曲線的最大值即為該地塊的現(xiàn)狀最大負(fù)荷值。

利用該方法統(tǒng)計(jì)的各地塊現(xiàn)狀負(fù)荷，已經(jīng)結(jié)合GIS平臺(tái)，應(yīng)用于該市供電局規(guī)劃部門日常工作中，包括地塊和網(wǎng)格負(fù)荷管理、用戶報(bào)裝管理、10 kV線路及變電站供電范圍劃分、變電站布點(diǎn)規(guī)劃、網(wǎng)架規(guī)劃等。

5.2 效率分析

該市2018年共有10 kV配變79 918臺(tái)，從政府規(guī)劃部門獲取各地區(qū)總規(guī)和控規(guī)，合計(jì)地塊51 174個(gè)。約8萬臺(tái)配變分布在5萬多個(gè)地塊上，無論配變數(shù)量還是地塊數(shù)量均規(guī)模巨大，依靠傳統(tǒng)方法人工篩查各地塊所含配變?nèi)萘?，進(jìn)而統(tǒng)計(jì)地塊現(xiàn)狀負(fù)荷，顯然效率太低。對(duì)比傳統(tǒng)方法和本文方法統(tǒng)計(jì)該市全市地塊現(xiàn)狀負(fù)荷所需時(shí)間，結(jié)果見表3。其中配變負(fù)載率法所需工時(shí)根據(jù)以往經(jīng)驗(yàn)按每個(gè)地塊1 min估算，線路負(fù)荷分配法所需工時(shí)根據(jù)以往經(jīng)驗(yàn)按每個(gè)地塊1.2 min估算，本文方法所需工時(shí)為實(shí)際操作真實(shí)數(shù)據(jù)。

圖6 主要行業(yè)日負(fù)荷曲線Fig.6 Daily load curves of main industries

表1 配變?nèi)肇?fù)荷Tab.1 Daily load curves of distribution transformers

由表3對(duì)比結(jié)果可知：與傳統(tǒng)方法相比，本文基于GIS技術(shù)的方法對(duì)于該市地塊負(fù)荷統(tǒng)計(jì)可節(jié)約75%～79%工時(shí)，按每天8 h計(jì)算約節(jié)省工作日149～189 d，按3人項(xiàng)目小組計(jì)算則可將工期從67～80 d縮減為17 d，工作效率有明顯提升。

5.3 可靠性分析

以某供電區(qū)域?yàn)槔?，區(qū)域1以商業(yè)、辦公負(fù)荷為主，區(qū)域2以居民負(fù)荷為主，區(qū)域3以工業(yè)負(fù)荷為主，現(xiàn)狀由3座變電站共10回10 kV線路對(duì)該區(qū)域供電,各區(qū)域用地性質(zhì)如圖7所示。

將利用本文方法對(duì)各區(qū)域的計(jì)算結(jié)果導(dǎo)入MapInfo，統(tǒng)計(jì)如圖8所示。

比較本文方法與傳統(tǒng)方法的結(jié)果見表4，計(jì)算結(jié)果誤差在0.4 MW以內(nèi)，誤差極小。本文方法顆粒度精細(xì)至配變，各配變最大負(fù)荷均由計(jì)量自動(dòng)化系統(tǒng)導(dǎo)出，利用負(fù)荷曲線疊加計(jì)算地塊負(fù)荷，擺脫了傳統(tǒng)統(tǒng)計(jì)方法中配變負(fù)載率一刀切和同時(shí)率主觀取值影響，結(jié)果更加客觀，計(jì)算機(jī)技術(shù)代替大量人工統(tǒng)計(jì)可有效提高統(tǒng)計(jì)結(jié)果可靠性。

表2 地塊日負(fù)荷Tab.2 Daily load curves of land parcels

表3 統(tǒng)計(jì)地塊負(fù)荷效率分析Tab.3 Analysis of statistical land load efficiency

圖7 區(qū)域地塊用地屬性Fig.7 Land use properties of regional parcels

圖8 現(xiàn)狀負(fù)荷統(tǒng)計(jì)Fig.8 Current load statistics

表4 地塊負(fù)荷計(jì)算結(jié)果比較Tab.4 Comparisons of calculation results of land parcel load MW

6 結(jié)束語

本文提出一種基于GIS技術(shù)的大規(guī)模地塊現(xiàn)狀負(fù)荷統(tǒng)計(jì)方法，解決了傳統(tǒng)方法過度依賴人工統(tǒng)計(jì)和同時(shí)率主觀取值的缺點(diǎn)，彌補(bǔ)了現(xiàn)有供電部門各數(shù)據(jù)系統(tǒng)無統(tǒng)計(jì)地塊負(fù)荷的不足，有效提高大規(guī)模地塊現(xiàn)狀負(fù)荷統(tǒng)計(jì)效率和結(jié)果可靠性，為配電網(wǎng)精細(xì)化規(guī)劃和日常負(fù)荷精細(xì)化管理提供基礎(chǔ)，具有一定實(shí)用價(jià)值。