基于一致性算法的“源-網-荷-儲”協同優化方法

張鑫，鄧莉榮，李敬光，王康，陳威洪，黎嘉樂

(1．廣東電網有限責任公司東莞供電局，廣東省東莞市 523000；2．清華-伯克利深圳學院，廣東省深圳市 518057)

0 引 言

隨著新一代智能電網技術的發展，可再生能源和儲能裝置大規模并網，需求響應技術不斷深化，電力系統由電源單邊傳輸轉向“源荷”雙邊互動，進一步轉變為“源-網-荷-儲”多參與主體的復雜系統，給現有的調度模式帶來了巨大的挑戰[1]。為了應對可再生能源的高波動性和間歇性，負荷側的不確定性[2]，基于“源-網-荷-儲”協調互動的調度方案，成為近年來學者們關注的熱點[1,3-4]。

通過充分利用電源、電網、負荷和儲能之間的靈活和互補特性，“源-網-荷-儲”協調互動模式可以提高系統應對風險的能力，增強動態功率平衡的控制能力[5]。盡管目前源-源[6]、源-網[7]、網-荷-儲[8]等方面已經開展了較為成熟的研究，但是全面考慮“源-網-荷-儲”4個要素的有功調度策略還沒有得到充分研究。

目前，“源-網-荷-儲”協同優化大多采用集中式優化方法[9]，如拉格朗日乘子法[10]，或者智能算法，如粒子群算法[1,11]。

隨著分布式算法的迅速發展，“源-網-荷-儲”協同優化調度問題可以通過使用分布式算法，減小模型的復雜度，提高求解效率。分布式算法具有保護用戶隱私，減少與中央控制器的通信量，避免單點故障造成系統崩潰等優點。常見的分布式算法有拉格朗日松弛法、交替乘子迭代法、輔助問題原理法、最優性條件分解、一致性算法等[12]。其中，一致性算法具有較高的靈活性，可以應用于任何級別的分區，例如，一個單獨的節點或者一個大的區域都可以成為一個分區[12]。因此，這一特性應用于“源-網-荷-儲”這種多元結構具有天然的優勢。系統中的每一個主體：電源、負荷、儲能等，都可以看作是一個代理，代理之間相互協作，可以大大提升系統總體運行效率。但目前，應用于“源-網-荷-儲”協同系統的一致性全分布式算法還不多見。

針對以上空白，本文考慮電源、負荷、儲能的響應特性以及風電的波動性，建立“源-網-荷-儲”最優資源分配模型，然后提出基于一致性算法的協同優化算法，最后通過算例驗證協同的必要性：促進可再生能源的消納吸收，引導彈性負荷參與需求響應。

列寧關于批評要有內容的思想是列寧主義的重要組成部分，是在俄國革命實踐中產生的黨的建設思想，有著自身鮮明的特點。它在俄共（布）黨內曾經產生積極影響，有著重要歷史意義。

1 一致性算法

一致性算法是指每個代理各自收集本地信息，并且和相鄰代理之間交互邊界信息，將問題分布式求解，同時保證與集中式求解算法解集的一致性[13-14]。根據一致性算法，代理i的信息更新過程表示如下：

(1)

式中：ξi(t+1)表示代理i在第t+1次迭代時的局部信息；ξj(t)表示代理j在第t次迭代時的局部信息；wij為代理i和j之間的通信系數；Ω為參與信息傳遞過程的總代理數。

我們采用平均metropolis算法確定系數wij：

彈性負荷參與協同，即考慮負荷具有響應特性，可隨著系統邊際成本的變化調整自身用電情況，從而增大風電消納，提高社會總體效益。

式中：ni和nj分別為連接到i和j的代理數；ε設為一個較小的數值；Ωi為與代理i相連通的代理索引。

2 “源-網-荷-儲”最優資源分配模型

電網中有功功率平衡約束為

施術準備階段幫助患者清潔皮膚進行消毒處理，可問詢患者疼痛耐受度，大部分患者無需麻醉，其他怕疼的患者進行心理疏導后放松心態接受診療，操作中將火針的針尖部位置于酒精燈之上逐漸燒紅，穩定其熱量快速的置于扁平疣部位，針柄用布包裹，以不導熱為宜，而后刺入患者的病患區域，待針稍有冷卻后繼續燒紅反復操作，一次將患者的扁平疣清除干凈，注意針要燒透，這樣滅菌效果好痛感小，垂直要直刺正中一下。刺入程度不宜太深，入皮最多1～2 mm。皮損較大的的情況下，利用三頭火針直刺，小部分的皮損可利用粗火針斜刺[2] 。

(3)

式中：PG,k，PW,r，PB,m和PD,l分別為傳統發電機k，可再生能源r，儲能m，負荷l的有功功率，假設向電網注入功率為正，向電網吸收功率為負；NG，NW，NB，ND分別為傳統機組、可再生能源、儲能、負荷的索引。

實際系統中，電源包括傳統燃煤發電機，風能、太陽能等可再生能源。可再生能源的運行成本可以忽略不計。由于其零邊際成本導致負荷優先由可再生能源供給，為了保證傳統發電機和儲能的可調度性，本文假設可再生能源出力不能完全供給負荷。傳統發電機的成本函數為

(4)

式中ak，bk，ck分別為機組k的燃料費用系數。

一般認為，用戶消耗的電量越大，完成的任務越多，效用越高[15]。因此，假設負荷的負利潤函數為

歧義詞通常可分為交集型、組合型和真歧義三類，不同的消除方法對歧義詞的類型有所偏重。消除歧義詞主要方法有兩種：一是基于規則的歧義消解算法，它主要采用語義、語法、詞性等規則對歧義字段進行消除；一是基于統計的歧義消解算法，它是統計相鄰字同時出現的頻次，確定組成詞的可能性，具體有：單純以詞頻為依據的詞頻法；多個單字之間相連趨勢為依據的t-測試法；兩字之間結合緊密程度為依據的互信息法。

(6)

彈性負荷具有一定的調節范圍：

失眠是指不能正常獲得睡眠的一種病癥,目前約有10%～20%的人群受失眠所困擾[1]。長期失眠不僅會妨礙人們的身體健康、降低人們的生活質量,同時還常伴發注意力不集中、精神障礙、記憶力減退等多種精神和軀體疾病[2-3]。目前臨床上用于治療失眠癥的藥物主要包括苯二氮卓類和巴比妥類,但這些藥物往往有不同程度的不良反應,且容易形成藥物依賴[4]。已有研究結果表明一些芳香植物精油具有改善睡眠的功效,且不良反應少,如薰衣草、纈草[5-6]。使用具有鎮靜催眠效果的芳香植物揮發性化合物來改善睡眠,是上述亞健康狀態的重要輔助療法,因此芳香植物精油的鎮靜催眠功效具有很大的研發價值。

Ploss,l=slPD,l

(8)

通過首次制作招平斷裂北段跨4礦區長達5000m的地質剖面圖，系統研究招平斷裂北段深部構造行跡，招平斷裂北段玲瓏金礦田內帚狀構造群和招平斷裂北段主干九曲蔣家208斷裂以及招平斷裂北段分支破頭青斷裂之間的空間關系得以清晰呈現(圖6)。

相應地，式(3)需要考慮損耗項：

(10)

根據文獻[16，19]，若σ足夠小，由式(16)—(19)構成的系統會逐步收斂到穩態。因此，λi(t)收斂到最優解λ*，失配量的局部估計值Pmis,i(t)會消失，即功率平衡約束滿足。

呂凌子的意思是，既然索賠沒了指望，就是借錢也得把這個事情給結了。肇事者是她的丈夫，丈夫雖然不是故意所為，但畢竟對一個無辜女人造成了難以啟齒的傷害，這種傷害也許是無辜女人一輩子的心理陰影。錢不是萬能的，然而除了錢，呂凌子實在想不出比這更好的解決辦法。

線路傳輸損耗通常占總負荷的3%～7%。為了描述損耗對總成本的影響，本文認為損耗可近似為負荷的線性函數[16]，用傳輸損耗因子sl表示：

雨心坐在病床邊，把手伸過去，捏握著蔣浩德干枯的手，看見他晚境如此凄涼，不禁感嘆：“人老了，確實要個伴兒，你出院后，搬到我家去住吧！”

民辦高校學生普遍存在學習基礎薄弱問題，并且上課注意力集中度較差，課堂時間利用率不高，時間管理能力有待加強。同時，整理和匯總理論知識，構建知識體系能力有待提高。另外，傳統的教育方式會統一每個班級學生的教學時間，教學場地仍然以教室、實訓室等實體場所為主，限制了學生自由學習的時間和地點。

圖1顯示了各參與主體：電源、負荷、儲能的成本函數或者負利潤函數。x表示輸出功率，y表示成本。作為生產者，電源的成本函數滿足產出與邊際成本成正比；對于負荷(輸出功率設為負)，從數值而言，負荷越大，負利潤越大，逐漸趨于飽和；對于儲能而言，儲能可以充電，可以放電，因此既是源，也是負荷。

圖1 系統各決策主體的成本函數/負利潤函數示意圖Fig.1 Cost function or minus utility function of each participant in the system

對于源側而言，效益函數為總售電收益與總成本之差；對于負荷側而言，效益函數為總利潤與總購電費用之差。從全局來看，“源-網-荷-儲”的優化調度模型需要協調各參與主體使得社會總效益最大，購售電費用在統一市場出清價格機制和電網功率平衡約束下相互抵消，則優化問題可表述為最大化源側負成本函數和負荷利潤函數：

假設各決策主體：發電機k、可再生能源r、負荷l、儲能m都有一個本地代理。為方便起見，代理i管理的功率統一寫成Pi，其邊際成本λi為

λi=aiPi+bi

(14)

若式(13)求解的最優邊際成本為λ*，根據等耗量微增量準則[18]，各代理可以通過式(15)安排出力計劃：

每一個代理只使用本地的信息以及和相連的領域代理交換信息。各代理的優化步驟如下：

3 基于一致性算法的協同優化策略

模型(13)可以通過集中式算法求解。集中式算法需要一個中央控制器；中央控制器和各代理之間進行信息交互[16]。中央控制器首先向各代理征集信息，包括電廠、負荷、儲能的運行情況和出力界限，成本/利潤函數等。然后控制器求解最優邊際成本，并下達給各代理。高滲透率的可再生能源給系統運行控制增加了難度：控制器需要更加頻繁地動作以響應動態變化的源端供給。這種運行工況的快速變化和不可預期性，集中式算法也許并不能及時響應。而全分布式算法以其靈活性，可靠性和可擴展性，不需要與中央控制器通信，可能成為“源-網-荷-儲”協同調度的有效解決方案。

1.2 試驗方法 試驗地點設置在長興仙山湖庫區，根據庫區及入庫河流水系特點，選取水深較淺、透明度較高的入庫河流近岸區域設置樣方，分別在庫區的北湖設置9個樣方N1～N9，平均水深1.5 m，透明度常年可見底；南湖設置3個樣方S1～S3，水平水深0.8 m，透明度常年可見底。樣方設置為5 m×2 m，四周立桿標記邊界。

式中：λi(t+1)表示第(t+1)次迭代時代理i的邊際成本；σ表示可調節收斂速度的步長；Pmis,i(t)為全局供給與需求失配量的局部估計。

對于儲能設備，成本函數和功率約束可以寫成如下形式[17]：

本地代理的算法框圖如圖2所示。代理首先按照式(17)和(18)更新本地信息，然后與鄰居節點交換局部功率的估計值；本地代理整合本地和鄰居信息中的邊際成本和局部功率估計值，根據式(16)和(19)求出更新后的邊際成本和失配量的局部估計。

圖2 代理的算法框圖Fig.2 Algorithm diagram of an agent

4 算例與分析

仿真系統(調整的IEEE-14節點系統)如圖3所示，其中包含4個傳統發電機，1個儲能裝置，1個風力發電機和9個負荷。風力發電機的可用功率曲線見圖4。各參與主體的詳細參數見表1。

圖3 4機-14節點的“源-網-荷-儲”協同系統及其通信拓撲Fig.3 4-generator 14-node “generation-grid-load-storage” system

圖4 風電功率曲線Fig.4 Wind power generation

系統的通信拓撲通常與當地的通信設施建設相關，可以與實際物理系統拓撲不同。本文假設通信節點i(即代理i)與節點i-1，i-2，i+1，i+2互聯(若節點i=1，則與系統中編號最后的2個節點通信，依此類推)。

4.1 彈性負荷參與協同

圖5給出了各機組的出力變化。節點邊際成本如圖6所示。可以看出，風電出力與系統邊際成本呈負相關。由于沒有阻塞，機組出力搭界等情況，各節點的邊際成本相同。傳統發電機組為獲得更多的收益，在高邊際成本時段會增加相應的電出力，因此傳統機組出力曲線與邊際成本變動曲線趨勢相同。

由式(8)和式(9)可知，一個量子位數量為m的量子個體可以表示2m個狀態，因此小種群的量子個體就可表示傳統方法下的大數量個體。量子門的操作使得量子進化能夠擁有很強的全局搜索能力，隨著概率幅的收斂，搜索的結果也會自動變為局部搜索。

表1算例系統的參數設置
Table1Parametersofthetestsystem

注：G—發電機；B—儲能裝置；L—負荷。

而從負荷側考慮(見圖7，負荷按照數值畫出，省略負號)，系統邊際成本直接影響用戶用電成本，在相對低邊際成本時段，可調負荷會相對增加以補償在高邊際成本時段減少的負荷。從數值上分析，負荷在邊際成本減少時將增大消耗量，反之，在邊際成本增加時，選擇減少使用量。

圖8顯示，基于源網荷協同調度模型，最終實現協同系統的功率平衡，即源側總功率等于負荷側總功率。

圖5 負荷參與協同時各決策主體的出力變化Fig.5 Each decision marker output update with demand response

圖6 負荷參與協同時的節點邊際成本Fig.6 Nodal incremental cost with demand response

圖7 負荷的功率變化Fig.7 Local demand update

圖8 負荷參與協同時的系統供求不平衡量估計Fig.8 System mismatch estimation with demand response

4.2 負荷不參與協同

負荷不參與協同時，可看成固定負荷。此時，系統吸納風電的能力減弱，系統總體靈活性降低。

經過大量的實驗表明，在螺栓連接較為緊密時，兩傳感器所采樣的振動信號基本以振動激勵的振動函數為主，此時兩信號除振動幅度以外，基本不存在其他差異，也反映出兩被聯件連接相當緊密受激勵部件完全跟隨主振器件發生振動。然而，當螺栓連接不緊密，出現松動時，其振動波形包含高次諧波，而高次諧波在本系統中不作為判定條件，屬于沒有意義的引入雜波，需對其進行濾除。本文選用FPGA進行數字濾波，其主要特點是對不同的應用場合可通過改變FPGA內濾波器的方式來適應不同的應用場合，例如高頻振動環境下的濾波并不適用于低頻振動環境下的濾波情況，而不需要更換硬件設施，只需將FPGA的代碼做更改即可適應不同工作環境[11]。

圖9給出了負荷不參與協同時各機組和儲能的出力變化。在60～300 min時，各機組以及儲能的出力達到下限，此時系統只能采取棄風限電措施，見圖10。同理，在第1 320～1 440 min時，系統也發生了棄風現象。說明負荷不參與需求響應時，系統吸納風電能力減弱。

圖9 負荷不參與協同時各機組和儲能的出力變化Fig.9 Local generator and storage output update without demand response

圖10 負荷不參與協同時的風電功率輸出Fig. 10 Wind power generation without demand response

棄風現象發生時，各節點的邊際成本不再完全一致(見圖11)。因為此時各機組出力達到下限，邊際成本按照各機組最小功率求取。圖12顯示系統供求平衡，算法收斂。

圖11 負荷不參與協同時的節點邊際成本Fig. 11 Nodal incremental cost without demand response

為了評估本文所提方法的有效性，表2和圖13給出了本文算法和粒子群算法的結果對比。粒子群算法作為一種智能算法，在凸問題中可以較好地搜索出全局最優解。從表2可以看出，本文所提的分布式算法與粒子群算法所得結果基本一致。從圖13中可以看出，本文所提方法在迭代92次后收斂，粒子群算法在迭代112次后收斂，因此在收斂精度相同的情況下，本文所提方法收斂速度更快。

圖12 負荷不參與協同時的系統供求不平衡量估計Fig. 12 System mismatch estimation without demand response

圖13 算法收斂次數對比Fig. 13 Methods comparison in convergence iterations

5 結 論

本文計及了可再生能源出力的波動性，彈性負荷的響應特性，考慮了網絡中的線損，建立了“源-網-荷-儲”的最優資源分配模型，并基于一致性算法對資源分配模型進行了求解。算例表明：考慮“源-網-荷-儲”協同的模型可以提高可再生能源的消納率，通過價格信號促進負荷進行需求響應；所采用的一致性算法計算速度快，收斂性好，適用于對實時性要求較高的動態電力系統控制，對多主體參與的實時市場運行也有一定的幫助。

2004年拍《歷史的天空》，殷桃初生牛犢，第一天第一場戲便是跟李雪健、張豐毅搭。這得多重要。殷桃就不會演了，很簡約的一個打招呼，NG了20多遍。