李蒙贊

(國網(wǎng)山西省電力公司 電力科學(xué)研究院，山西 太原 030001)

工業(yè)園區(qū)屬于工業(yè)企業(yè)相對集中的城市功能區(qū)域，區(qū)域內(nèi)大負荷用電企業(yè)較為密集，部分企業(yè)采用多電源回路上的專用接入點接入，或在2條以上多電源回路上部署多個專用接入點接入；部分企業(yè)采用多路專線供電接入，或?qū)＞€供電配合多電源回路接入。此處多電源回路包括拉手回路、環(huán)形回路等，拉手回路一般擁有至少2個電源接入點，環(huán)形回路一般擁有至少3個電源接入點，即多電源回路上的電源性質(zhì)本身即為2級負荷電源[1]。

綜合能源信息管理管理系統(tǒng)是區(qū)域用電負荷調(diào)度的重要管理信息系統(tǒng)，通過將特定區(qū)域內(nèi)的用電負荷需求統(tǒng)一管理，對不同電壓級別，不同供電回路的電源狀態(tài)、衍生電源狀態(tài)進行統(tǒng)一管理，以確保重點區(qū)域的供電可靠性[2]。

研究采用平衡調(diào)度理念，將某個案工業(yè)園區(qū)的綜合能源管理信息系統(tǒng)的調(diào)度專家系統(tǒng)算法進行充分優(yōu)化，使其策略編制功能更貼近實際需求。

1 個案基本情況

該個案工業(yè)園區(qū)面積7.4平方km，為省級工業(yè)開發(fā)區(qū)，園區(qū)內(nèi)擁有用電企業(yè)12個，用電總?cè)萘繛?8.1 MW，其中35 kV專線用電容量12.4 MW，共2個用戶；10 kV專線用電容量6.5 MW，共3個用戶；10 kV多電源回路雙節(jié)點接入用戶用電容量6.0 MW，共3個用戶；10 kV多電源回路單節(jié)點接入用戶用電容量3.2 MW，共4個用戶。具體線路分布情況如表1所示。

表1 園區(qū)內(nèi)供電線路分布示意圖Tab.1 Distribution diagram of power supply lines in the park

由表1可知，園區(qū)內(nèi)供電線路包括針對5個用戶的10條供電專線，其中35 kV專線4條，10 kV專線6條，另布置3條10 kV拉手線路[3]。因為工業(yè)園區(qū)拉手線路中每個用電節(jié)點的用電負荷需求均較高，所以其拉手回路內(nèi)的接入點密度遠小于居民區(qū)接入點密度，3條拉手線路電源分別來自35 kV變電所的105站、108站、116站、109站共4個變電站。10 kV專線電源分別來自105站、108站、116站共3個變電站，專線電源站包含在拉手線路變電站內(nèi)。35 kV專線電源分別來自35 kV的303站、307站、313站共3個變電站。

涉及到3個35 kV電源站和4個10 kV電源站，該7個電源站均分布在工業(yè)園內(nèi)或工業(yè)園周邊。所有用戶分為4個級別，分別為：35 kV雙專線供電、10 kV雙專線供電、拉手線路雙節(jié)點供電、拉手下單節(jié)點供電。該4個級別均可實現(xiàn)2級負荷直接供電[4]。

2 平衡調(diào)度算法設(shè)計

2.1 平衡調(diào)度算法的數(shù)學(xué)模型

CLASS分析法是將用戶負荷和電源站負荷賦予數(shù)學(xué)屬性，將客觀實體轉(zhuǎn)化為數(shù)學(xué)實體的研究方法。在對該個案的研究中，分析12個用電用戶的最大用電容量，假定其最大用電容量為PU，則其可能產(chǎn)生負荷在Rand(0,PU)區(qū)間內(nèi)，即采用隨機數(shù)模擬法可以產(chǎn)生12個用電用戶分別用電的負荷量[5]。而電源站的負荷量為其向用戶供應(yīng)電能的回路通斷情況，即對于10 kV電源站，如式(1)：

PΣ=∑Rand(0,PUi)

(1)

式中：PUi為第i個接通負荷的最大負荷容量；Rand(*)為隨機數(shù)取值過程，自變量為隨機數(shù)取值空間。

對于35 kV電源站，如式(2)：

(2)

假定所有7個電源站的負荷量均為PΣi，其中i=7，要實現(xiàn)該7個電源站的負荷均衡，則需要對其標準偏差率進行計算，如式(3)：

(3)

2.2 平衡調(diào)度算法下的開關(guān)量提取

該模型下涉及到不同電壓等級的多種負荷，其電氣連接方式包括各種變壓器、互感器等線圈繞組設(shè)備，隔離開關(guān)、接地開關(guān)、斷路器等開關(guān)設(shè)備，且母線結(jié)構(gòu)較為復(fù)雜，如果直接研究上述系統(tǒng)，勢必造成研究對象的復(fù)雜化導(dǎo)致無法提取對應(yīng)的數(shù)學(xué)模型[6]。而該研究的本質(zhì)只考慮到上級電源站對下級電源站的供電回路和對各負荷的供電回路，對任何供電回路來說，其回路兩側(cè)斷路器以及邊界回路的對應(yīng)斷路器，可以看做是其唯一功能節(jié)點。所以，為了研究平衡調(diào)度算法的相關(guān)屬性，將上述系統(tǒng)充分簡化，即形成一套最簡開關(guān)量構(gòu)成拓撲[7]。對網(wǎng)絡(luò)拓撲進行拓撲學(xué)模型提取，具體如圖1所示。

圖1 園區(qū)平衡調(diào)度所用開關(guān)拓撲關(guān)系圖Fig.1 Topological diagram of switches used for balanced dispatching in the park

從圖1可以看出，共涉及到33個開關(guān)，其中互斥開關(guān)共12對，即自由控制開關(guān)共21個。根據(jù)組合定律，上述開關(guān)量共可能形成開關(guān)組合模式為221=2.097×106種倒閘狀態(tài)，而其中經(jīng)常使用的倒閘狀態(tài)遠小于該值，多數(shù)倒閘狀態(tài)會造成該網(wǎng)絡(luò)內(nèi)大面積失電，屬于故障倒閘狀態(tài)[8]。所以，通過機器學(xué)習(xí)算法，在2.097×106種倒閘狀態(tài)中選擇出最適合當(dāng)前平衡調(diào)度需求的倒閘狀態(tài)，避免可能發(fā)生的故障倒閘狀態(tài)，是本文研究的重點[9]。

2.3 神經(jīng)網(wǎng)絡(luò)的設(shè)計

將上述21個開關(guān)節(jié)點設(shè)定為1個一維矩陣，如式(4)：

S={Si=(1,0),i∈1,2,…,n|n=21}

(4)

對上述7個母線負荷設(shè)定為1個一維矩陣，如式(5)：

P={Pi=Rand(0,PUi),i=1,2,…,n|n=7}

(5)

通過向神經(jīng)網(wǎng)絡(luò)中輸入S與P矩陣，得到最終的S矩陣判斷結(jié)果，將S矩陣判斷結(jié)果輸入到后置仿真模塊中，判斷其產(chǎn)生的P矩陣結(jié)果，比較P矩陣結(jié)果的標準偏差率，當(dāng)P矩陣標準偏差率小于閾值時，則輸出該結(jié)果到倒閘專家系統(tǒng)；具體結(jié)果如圖3所示。

圖2 神經(jīng)網(wǎng)絡(luò)及外圍模塊數(shù)據(jù)流圖Fig.2 Data flow diagram of neural network and peripheral modules

該神經(jīng)網(wǎng)絡(luò)中，輸入模塊共21+7=28個，為21個二值化變量(Logical格式)和7個雙精度浮點變量(Double格式)，輸出模塊共21個，為21個二值化變量(Logical格式)。21個輸出模塊分別直接從所有28個輸入模塊讀取信息，經(jīng)過數(shù)據(jù)處理后輸出二值化分析結(jié)果，即該神經(jīng)網(wǎng)絡(luò)是一族21列分列單模塊神經(jīng)網(wǎng)絡(luò)架構(gòu)[10]。

分析21個神經(jīng)網(wǎng)絡(luò)模塊屬性，在統(tǒng)計學(xué)上需要對28個輸入數(shù)據(jù)進行充分降維融合，7個雙精度浮點變量即P數(shù)列因子，輸入數(shù)據(jù)之間需要將數(shù)據(jù)投影到[0,1]空間，以避免其量綱差異和投影區(qū)間差異對[0,1]區(qū)間上的二值化數(shù)據(jù)帶來影響，其重投影函數(shù)如式(6)：

(6)

式中：Pi為數(shù)列中第i個母線的負荷值；maxP為7個母線負荷值中的最大值；Yi為重投影的輸出值；當(dāng)Pi=0時，Yi=0。

該神經(jīng)網(wǎng)絡(luò)模塊前3層隱藏層的作用為數(shù)據(jù)降維，且在降維時將損失信息量保存到數(shù)據(jù)節(jié)點中，所以數(shù)據(jù)節(jié)點的待回歸系數(shù)應(yīng)足夠豐富且其曲線拐點豐度應(yīng)滿足復(fù)雜曲線的擬合要求，故應(yīng)采用多項式函數(shù)進行節(jié)點設(shè)計，其基函數(shù)如式(7)：

(7)

式中：Xi為第i個輸入數(shù)據(jù)；Y為輸出數(shù)據(jù)；j為多項式階數(shù)；Aj為第j階多項式的待回歸系數(shù)。

為實現(xiàn)充分降維，前3層隱藏層的節(jié)點數(shù)量應(yīng)分別布置為29、11和5個，該45個節(jié)點均采用式(5)的多項式回歸函數(shù)進行節(jié)點設(shè)計，因為每個節(jié)點包含6個待回歸系數(shù)，則該45個節(jié)點共包含270個待回歸系數(shù)[11]。

第4、5層隱藏層的統(tǒng)計學(xué)意義是將上述數(shù)據(jù)充分二值化，以實現(xiàn)最后形成1個充分二值化的輸出值；故第4層布置3個二值化節(jié)點，第5層布置1個二值化節(jié)點，二值化回歸函數(shù)的基函數(shù)可寫為式(8)：

(8)

式中：Xi為第i個輸入數(shù)據(jù)；Y為輸出數(shù)據(jù)；e為自然常數(shù)，此處取近似值e=2.718 281 828；A為待回歸系數(shù)[12]。

經(jīng)過2層二值化輸出的1個二值化結(jié)果，因為存在多個二值化值的疊加，所以其二值化效果仍不可靠，所以在數(shù)據(jù)訓(xùn)練過程中，確保連續(xù)50次未出現(xiàn)所有數(shù)據(jù)二值化結(jié)果未落點在[0.05,0.95]后，方認為該神經(jīng)網(wǎng)絡(luò)已經(jīng)充分收斂。

3 仿真測試結(jié)果

數(shù)據(jù)訓(xùn)練過程中，選擇該工業(yè)園區(qū)自2018年1月1日至2019年12月31日實際運行數(shù)據(jù)作為訓(xùn)練數(shù)據(jù)，在該時間戳范圍內(nèi)隨機抽取真實運行數(shù)據(jù)對神經(jīng)網(wǎng)絡(luò)進行訓(xùn)練，使神經(jīng)網(wǎng)絡(luò)得到充分收斂[13]。

仿真測試選用電力CAE平臺下的SimuWorks仿真控件，采用隨機負荷法對上述模型進行仿真分析。參照組選用該工業(yè)園之前使用的早期平衡調(diào)度平臺算法，采用同一數(shù)據(jù)對早期平衡調(diào)度平臺算法和革新平衡調(diào)度算法進行比較分析[14]。

2種算法平衡負荷后，實際負荷變化率曲線的標準偏差率σ如式(3)和最大偏差比如式(9)：

(9)

式中：maxP、minP分別為7個母線負荷值的最大值和最小值。

該σ值和ΔP值仿真比較結(jié)果如表2所示。

表2 系統(tǒng)控制精度仿真比較結(jié)果表Tab.2 Comparison results of system control accuracy simulation

由表2可知，革新算法平衡后的7個個案變電站負荷均得到了更為有效的平衡，革新算法較早期算法的標準偏差率提升5.1倍，最大偏差比提升3.8倍。標準偏差率和最大偏差比取值減小，標志著7個個案變電站的電源母線負荷的均一性得到更大程度的統(tǒng)一，電網(wǎng)負荷的平衡性更佳[15-16]。標準偏差率和最大偏差比關(guān)系到電網(wǎng)節(jié)點的調(diào)度控制精度，早期全人工編制倒閘策略，人工執(zhí)行倒閘的過程中，因為倒閘頻率較低，倒閘過程的復(fù)雜性和現(xiàn)場管理風(fēng)險導(dǎo)致現(xiàn)場無法對變電站節(jié)點的負荷比給出精密控制。

為考察該目標的實際控制效果，使用電力CAE平臺下的SimuWorks仿真控件產(chǎn)生12個用電用戶負荷的連續(xù)時域變化量，比較電網(wǎng)內(nèi)的電能質(zhì)量，結(jié)果如表3所示。

表3 電能質(zhì)量控制效果比較表Tab.3 Comparison of power quality control effects

由表3可知，無功增壓比標志著工業(yè)園配電網(wǎng)絡(luò)中無功功率對電能質(zhì)量的影響；電壓峰值波動比標志著工業(yè)園內(nèi)的沖擊電壓對電能質(zhì)量的影響，該2項指標的評價結(jié)果中，革新算法的表現(xiàn)顯著優(yōu)于早期算法。而在革新算法的驅(qū)動下，開關(guān)倒閘頻率較早期算法略有增加，無功補償器投切頻率基本保持不變[17]。

4 結(jié)語

革新算法將倒閘頻率提升43.9%的前提下，實現(xiàn)無功增壓比和電壓峰值波動比分別下降71.2%和50.6%，同時實現(xiàn)7個電源站負荷情況最大偏差比提升3.8倍的控制效果。經(jīng)過該技術(shù)革新，變電站的負荷管理控制過程更加精細化，對用電側(cè)提供的電能波形更加接近于理論波形，工業(yè)園企業(yè)獲得的電能質(zhì)量服務(wù)更加負荷智能電網(wǎng)高質(zhì)量純凈電能供應(yīng)的技術(shù)要求，即，該算法對實現(xiàn)工業(yè)園區(qū)用電環(huán)境的平衡調(diào)度有積極意義。