多微電網(wǎng)的拓撲設(shè)計說到底是決定微電網(wǎng)之間可交換功率的連接結(jié)構(gòu)。用約束差分進化算法做大規(guī)模矩陣優(yōu)化是目前在Matlab里比較實用的一套拓撲設(shè)計技術(shù)路線。這篇文章把建模、編碼、約束處理和代碼實現(xiàn)完整梳理一遍也記錄一些實戰(zhàn)中才會暴露的坑。這套方案適合正在做微電網(wǎng)群規(guī)劃、園區(qū)多微電網(wǎng)互聯(lián)設(shè)計的工程師或研究生。如果你已經(jīng)跑過一些啟發(fā)式優(yōu)化算法但被高維矩陣編碼、約束不可行解、收斂振蕩這些問題卡住那這篇內(nèi)容應該能直接派上用場。我會把矩陣變量怎么構(gòu)造、約束怎么處理、參數(shù)怎么調(diào)、代碼怎么組織都展開講。1. 多微電網(wǎng)拓撲優(yōu)化到底在解決什么問題1.1 從孤島微電網(wǎng)到微電網(wǎng)群拓撲設(shè)計的價值一個微電網(wǎng)自給自足不難難的是在極端天氣、機組檢修或負荷突增時仍然可靠。多個微電網(wǎng)通過聯(lián)絡線互聯(lián)之后就可以在空間上共享備用容量某個區(qū)域光伏大發(fā)、負荷偏低時可以把富余功率送到隔壁隔壁負荷高峰、本地機組不夠時又能從其他區(qū)域接收支援。這種互聯(lián)結(jié)構(gòu)就是多微電網(wǎng)的拓撲。拓撲不是越密越好。每增加一條聯(lián)絡線意味著一次投資、一份運維成本和一部分網(wǎng)損。但線路太少又會導致功率互濟能力不足一個節(jié)點故障可能讓整片區(qū)域失去支援。所以拓撲設(shè)計本質(zhì)上是在經(jīng)濟性、可靠性和運行靈活性之間做最優(yōu)折中。工程上這類問題經(jīng)常被簡化成“手動方案比選”靠經(jīng)驗列幾個備選方案算潮流、算可靠性然后挑一個。當微電網(wǎng)數(shù)量達到幾十個、候選聯(lián)絡線數(shù)量達到數(shù)百條時手算完全不現(xiàn)實。必須把拓撲設(shè)計轉(zhuǎn)成一個數(shù)學優(yōu)化問題交給算法去搜。1.2 為什么優(yōu)化問題會變成矩陣優(yōu)化假設(shè)系統(tǒng)里有N個微電網(wǎng)節(jié)點節(jié)點之間的連接關(guān)系天然可以用一個N×N的矩陣來表示。比如矩陣第i行第j列的元素x_ij代表節(jié)點i和節(jié)點j之間聯(lián)絡線的狀態(tài)或功率。對角線元素表示本地信息非對角線元素表示互聯(lián)關(guān)系。直接把拓撲作為優(yōu)化變量時決策矩陣就是鄰接矩陣或功率傳輸矩陣。這個矩陣一般是對稱的因為微電網(wǎng)之間的聯(lián)絡線是無向的x_ij和x_ji描述的是同一條線路。實際操作時只需要優(yōu)化上三角部分的N(N-1)/2個元素下三角通過對稱性補全。這就引出了“大規(guī)模矩陣優(yōu)化”的問題。當N20時決策維度是190維N30時維度變成435維N50時直接沖到1225維。如果再考慮多個典型運行場景、多種線路容量選項決策變量和約束數(shù)量還會進一步膨脹。普通的窮舉法、枚舉法到這里基本都是死路。1.3 這類問題難在什么地方第一個難點是非凸性。拓撲變量里包含離散的“建/不建”決策功率分配又是連續(xù)變量二者耦合在一起。目標函數(shù)里還有建設(shè)成本這種帶階躍性質(zhì)的項整個優(yōu)化曲面既不平滑也不凸。第二個難點是強約束。每個節(jié)點都要滿足功率平衡每條聯(lián)絡線的傳輸功率不能超過容量上限整個系統(tǒng)還要滿足一定的可靠性要求。這些約束不是軟性的違反任何一條方案在工程上就是廢的。第三個難點是維度高。前面算過決策變量隨節(jié)點數(shù)平方增長。差分進化這類群體算法在高維空間里很容易失去多樣性搜索效率斷崖式下降。如果編碼方式再不合理算法甚至會花大量時間在“全不可行解”的區(qū)域內(nèi)打轉(zhuǎn)。2. 約束差分進化算法為什么選它核心機制是什么2.1 差分進化算法的基本流程差分進化算法Differential Evolution, DE是經(jīng)典進化算法里比較皮實的一種。它的核心思想其實很樸素讓種群個體之間互相“借鑒”差異信息生成新個體然后跟舊個體競爭上崗。標準流程分四步。第一步是初始化隨機生成NP個個體每個個體代表一個候選解。第二步是變異從種群中隨機挑三個不同個體用其中兩個的差乘以縮放因子F加到第三個個體上得到變異個體。第三步是交叉變異個體與原個體按位混合生成試驗個體混合程度由交叉概率CR控制。第四步是選擇比較試驗個體和原個體留下更好的一個。整個過程不需要梯度信息對目標函數(shù)是否可導完全不敏感。這一特性對多微電網(wǎng)拓撲優(yōu)化這種非凸、離散連續(xù)混合的問題非常合適。只要能把方案編碼成向量DE就能搜。用生活化的方式理解每個個體就是一張“連線方案圖”。變異操作相當于參考幾個已有方案拍腦袋組合出一個新方案交叉操作相當于把兩個方案拼拼湊湊你取一段我取一段選擇操作就是把好方案留下壞方案淘汰。反復迭代下去種群整體越來越接近最優(yōu)連線方案。2.2 約束處理分步篩選比一票否決更靠譜多微電網(wǎng)拓撲優(yōu)化最大的敵人不是目標函數(shù)復雜而是約束太多導致可行解稀少。如果約束處理方式不對算法一開始生成的大多數(shù)個體都是不可行的進化過程很快會退化。常見的做法有三種罰函數(shù)法、可行性規(guī)則、ε約束法。罰函數(shù)法最直觀把約束違反量乘一個懲罰系數(shù)加到目標函數(shù)上但懲罰系數(shù)非常難調(diào)。調(diào)小了不可行解輕松勝出調(diào)大了搜索又被過早牽進可行域邊緣找不到真正好的方案?？尚行砸?guī)則是更穩(wěn)妥的選擇。它按三條準則比較兩個個體可行解永遠優(yōu)于不可行解同為可行解目標函數(shù)值小的勝同為不可行解約束違反總量小的勝。這套準則不需要任何懲罰系數(shù)簡單有效是約束差分進化算法里最常用的處理框架。ε約束法可以看成可行性規(guī)則的升級版。它的想法是在前期允許個體有一定程度的約束違反讓種群在早期能保持足夠的探索活力而不是一上來就被嚴格約束鎖死。隨著迭代代數(shù)的增加ε值逐漸減小約束一步步收緊最后收斂到完全可行區(qū)域。2.3 面向矩陣編碼的算法調(diào)整直接把標準DE套在矩陣問題上很容易翻車。原因在于DE的變異算子是用實數(shù)向量的線性加減來生成新個體如果決策向量里混著0和1的拓撲決策差分變異產(chǎn)生的中間值會大量落在(0,1)之間導致拓撲處于“半連接不連接”的模糊狀態(tài)。一個實用的調(diào)整是決策變量先用連續(xù)值表示聯(lián)絡線的功率容量或傳輸功率水平用一個閾值τ來判斷線路是否實際建設(shè)。比如x_ij小于τ就視為不連接大于等于τ就視為連接。這樣矩陣既保留了連續(xù)變量的優(yōu)化特性又能落地成工程可用的離散拓撲。另一個調(diào)整是對稱性修復。矩陣優(yōu)化中每次調(diào)用適應度函數(shù)之前都要把決策向量還原成對稱矩陣強制對角線為零。不然變異和交叉會把原本對稱的矩陣破壞掉x_ij和x_ji出現(xiàn)不一致目標函數(shù)計算直接出錯。3. Matlab代碼實現(xiàn)從建模到收斂的完整拆解3.1 仿真場景與模型參數(shù)我先用一個20節(jié)點多微電網(wǎng)算例說明整個實現(xiàn)流程。每個節(jié)點內(nèi)部包含一臺燃氣機組、一個光伏電站和一片負荷區(qū)域。節(jié)點之間通過候選聯(lián)絡線連接每條聯(lián)絡線有傳輸容量上限。優(yōu)化目標是在滿足全天不同時段功率平衡的前提下最小化年綜合成本。綜合成本包括三部分聯(lián)絡線建設(shè)成本的等年值、燃氣機組的發(fā)電成本、線路傳輸過程中的能量損耗成本。發(fā)電成本用二次函數(shù)近似網(wǎng)損成本跟傳輸功率的平方掛鉤建設(shè)成本用Sigmoid函數(shù)平滑逼近避免階躍項給搜索過程引入不連續(xù)點。參考參數(shù)設(shè)置如下表所示。參數(shù)取值說明節(jié)點數(shù)N20微電網(wǎng)節(jié)點數(shù)量單節(jié)點負荷峰值300~800 kW按區(qū)域差異隨機生成光伏裝機200~600 kW午后出力明顯聯(lián)絡線候選容量200~500 kW不同位置上限不同線路等年建設(shè)成本8~15 萬元/km含折損和維護燃氣發(fā)電成本系數(shù)0.35~0.55 元/kWh二次函數(shù)線性項3.2 種群初始化與矩陣解碼差分進化的種群個體在算法內(nèi)部是一維向量但工程問題是一個對稱矩陣。所以第一件事就是把向量和矩陣之間的轉(zhuǎn)換寫好。function X vec2sym(v) D length(v); n round((1 sqrt(1 8 * D)) / 2); X zeros(n, n); k 1; for i 1:n-1 for j i1:n X(i,j) v(k); X(j,i) v(k); k k 1; end end end初始化時要注意如果所有元素都在[0, P_max]之間隨機均勻分布那初始矩陣會非常稠密幾乎每個節(jié)點都互相“弱連接”拓撲結(jié)構(gòu)完全沒有可區(qū)分度。我建議初始值大部分落在[0, τ]區(qū)間內(nèi)即默認不連接只讓少部分元素進入建設(shè)區(qū)間。function [X0] init_population(NP, D, Pmax_upper) X0 zeros(D, NP); for i 1:NP % 稀疏初始化絕大多數(shù)值在[0, 0.02*Pmax_upper] v rand(D, 1) * 0.02 * Pmax_upper; % 隨機選 2~5 個位置給出真實線路候選值 idx randperm(D, randi([2, 5])); v(idx) Pmax_upper * (0.3 0.7 * rand(length(idx), 1)); X0(:, i) v; end end這樣做的好處是初始種群里的拓撲結(jié)構(gòu)天然稀疏不會一開始就把搜索空間浪費在“全連接”這種工程無效方向上。3.3 目標函數(shù)與約束計算目標函數(shù)負責把決策向量轉(zhuǎn)化成成本數(shù)值。實現(xiàn)時先把向量恢復成矩陣再根據(jù)矩陣中非零線路計算建設(shè)成本、運行成本和網(wǎng)損。function [cost] objective_func(v, model) X vec2sym(v); n size(X, 1); % 建設(shè)成本項Sigmoid平滑逼近階躍 inv_cost 0; for i 1:n-1 for j i1:n if X(i,j) model.tau inv_cost inv_cost model.inv_unit(i,j) / (1 exp(-10*(X(i,j)/model.Pmax(i,j) - 0.5))); end end end % 運行成本項 op_cost sum(model.alpha .* model.Pg.^2 model.beta .* model.Pg model.gamma); % 網(wǎng)損成本項 loss_cost model.loss_coef * sum(sum(X.^2)); cost inv_cost op_cost loss_cost; end約束計算需要返回每個個體的約束違反總量。這里要特別小心每條約束不是簡單判斷“是否違反”而是要量化“違反了多少”否則可行性規(guī)則沒法比較兩個不可行解的優(yōu)劣。function [viol] constraint_violation(v, model) X vec2sym(v); n size(X, 1); viol 0; % 線路容量約束 for i 1:n-1 for j i1:n if X(i,j) model.Pmax(i,j) viol viol (X(i,j) - model.Pmax(i,j)) / model.Pmax(i,j); end end end % 節(jié)點功率平衡約束簡化版本用功率不平衡量衡量 net model.Pg model.Ppv - model.Pload; for i 1:n balance net(i); for j 1:n if i ~ j balance balance model.incidence(i,j) * X(i,j); end end if abs(balance) model.balance_tol viol viol abs(balance) / model.balance_tol; end end end每個約束的偏差量做歸一化后累加這樣不同量綱的約束才能放到同一個標量里比較。否則線路容量的千瓦級偏差會直接把功率平衡的微小偏差淹沒。3.4 主循環(huán)與約束支配選擇主循環(huán)按“變異-交叉-選擇”的順序迭代。變異采用經(jīng)典的DE/rand/1策略交叉用二項式交叉。選擇階段采用可行性規(guī)則。for gen 1:MaxGen for i 1:NP % 隨機選取三個互不相同的個體 r randperm(NP, 4); r(r i) []; r1 r(1); r2 r(2); r3 r(3); % 變異 V X(:, r1) F * (X(:, r2) - X(:, r3)); % 邊界處理 V min(max(V, LB), UB); % 二項式交叉 jrand randi(D); U X(:, i); mask rand(D, 1) CR; mask(jrand) true; U(mask) V(mask); % 可行性規(guī)則選擇 f_new objective_func(U, model); v_new constraint_violation(U, model); f_old objective_func(X(:, i), model); v_old constraint_violation(X(:, i), model); if (v_new model.eps v_old model.eps f_new f_old) || ... (v_new model.eps v_old model.eps) || ... (v_new model.eps v_old model.eps v_new v_old) X(:, i) U; f_pop(i) f_new; v_pop(i) v_new; end end % 記錄最優(yōu)解 [best_f, best_idx] min(f_pop 1e6 * (v_pop model.eps)); best_v X(:, best_idx); curve(gen) best_f; % ε逐漸收緊 model.eps model.eps0 * (1 - gen / MaxGen) ^ 2; end選擇邏輯說明兩條都可行時取成本低的新個體可行而舊個體不可行新個體獲勝兩條都不可行時取約束違反量小的。這個邏輯正是Deb可行性規(guī)則的標準實現(xiàn)不需要任何懲罰系數(shù)。3.5 參數(shù)表與推薦取值差分進化雖然皮實但參數(shù)設(shè)置不當還是會出問題。下表是我在多次測試中比較穩(wěn)定的取值范圍按20~50節(jié)點規(guī)模調(diào)整。參數(shù)推薦范圍備注種群規(guī)模NP80~150維度高時可適當增大縮放因子F0.4~0.7過大易振蕩過小易早熟交叉概率CR0.7~0.9影響局部探索強度最大迭代代數(shù)800~1500看收斂曲線趨勢拓撲閾值τ0.05~0.1相對線路容量ε初始值0.1~0.3前期允許一定違反我的經(jīng)驗是NP不要低于維度D的1/5否則種群多樣性不夠。F取0.5左右比較穩(wěn)CR取0.8能讓新個體保留更多變異信息。ε的衰減曲線用二次下降比線性下降效果更好前中期探索空間更充分。4. 典型算例結(jié)果與關(guān)鍵觀察4.1 小規(guī)模算例驗證正確性先做5節(jié)點小算例這個規(guī)模下可以把所有可能的拓撲組合枚舉出來用窮舉結(jié)果校驗差分進化算法是否找得到最優(yōu)解。我把候選線路限制在10條以內(nèi)這樣2^101024種組合枚舉很快。對比結(jié)果顯示約束差分進化算法在多數(shù)隨機種子下都能找到理論最優(yōu)解少數(shù)幾次會卡在次優(yōu)解。這說明算法實現(xiàn)本身的正確性沒有問題。這一步一定不要跳。直接上大算例一旦結(jié)果不合理根本說不清是模型錯還是算法錯。先用小規(guī)模把代碼邏輯驗證干凈再放大節(jié)點規(guī)模能省去大量排查時間。4.2 大規(guī)模算例性能表現(xiàn)20節(jié)點算例運行800代收斂曲線呈現(xiàn)出明顯的前期陡降、后期趨平趨勢。前200代成本下降最快400代之后基本進入平臺期。最終優(yōu)化拓撲相比等容量隨機布線的方案綜合成本下降約15%到20%。30節(jié)點算例的決策維度是435維種群規(guī)模取120運行1000代在普通配置的臺式機上耗時約40分鐘。從最終拓撲看算法傾向于保留一組“主干聯(lián)絡線”形成兩到三個環(huán)狀支撐結(jié)構(gòu)而不是把所有節(jié)點全部互聯(lián)。這種結(jié)果符合工程直覺少量關(guān)鍵線路承擔主要互濟功能冗余線路的投資收益比不高。從約束滿足情況看最終解的功率平衡約束和線路容量約束全部滿足說明可行性規(guī)則配合ε衰減機制成功把種群引向可行區(qū)域。4.3 參數(shù)敏感性關(guān)鍵結(jié)論我用控制變量法重點測了F、CR和閾值τ三個參數(shù)?？s放因子F從0.3調(diào)到0.9的過程中0.5附近收斂速度和最終精度最均衡。F0.3時早期收斂快但150代之后容易陷入局部最優(yōu)F0.8以上振蕩明顯收斂曲線經(jīng)常出現(xiàn)回跳。交叉率CR在0.8附近表現(xiàn)最好。CR太低試驗個體跟原個體太像缺乏探索性CR太高矩陣結(jié)構(gòu)的連續(xù)性被破壞拓撲經(jīng)?！爸щx破碎”。后者在矩陣編碼問題上尤其明顯。閾值τ的影響容易被低估。τ太小很多微弱的“偽連接”會被當成有效線路最終方案拓撲復雜、建設(shè)成本虛高τ太大真實有用的線路被砍掉可靠性下降。建議τ取線路容量上限的5%左右并在優(yōu)化完成后做一次連通性校驗。5. 實戰(zhàn)中踩過的坑與排查方法5.1 收斂曲線不下降或者振蕩劇烈這個問題的常見原因是種群中大部分個體不可行算法只能在不可行域里互相比較約束違反度。表現(xiàn)就是cost曲線長時間橫盤或者上下亂跳。排查時先看每個個體進入選擇階段前的約束違反總量分布。如果90%以上的個體viol都大于0說明約束太嚴或初始化太激進。這時調(diào)大ε初始值或者改進初始化策略讓初始種群里多生成一些結(jié)構(gòu)簡單、容易滿足約束的方案。5.2 優(yōu)化結(jié)果矩陣幾乎全零矩陣全零意味著算法認為“一條線路都不建最劃算”。這通常是建設(shè)成本系數(shù)設(shè)置過高或者目標函數(shù)中缺少可靠性維度的約束。多微電網(wǎng)拓撲優(yōu)化的意義本來就建立在“必須保證一定互聯(lián)水平”的工程前提上。解決辦法有兩個方向一是給系統(tǒng)增加可靠性約束比如要求每個節(jié)點至少有一條外聯(lián)路徑或孤島負荷覆蓋率不低于某個閾值二是在目標函數(shù)中加入線路數(shù)量的懲罰系數(shù)讓算法在“省投資”和“保互濟”之間自然權(quán)衡。5.3 矩陣幾乎全連接拓撲失去稀疏性這個坑跟5.2相反出現(xiàn)在省略建設(shè)成本或建設(shè)成本系數(shù)過低時。另一個誘因是閾值τ設(shè)得太小導致大量本該歸零的低功率元素被判定為有效線路。我的處理方式是在目標函數(shù)里加一個“線路稀疏獎勵項”即統(tǒng)計有效線路數(shù)量再乘一個較小的單位成本系數(shù)。同時在最終方案輸出前對所有線路做一次容量利用率篩選如果某條線路在多個典型場景下的平均傳輸功率都低于其容量的10%就把它剔除重新驗算。5.4 運行耗時過長大規(guī)模矩陣優(yōu)化的計算量主要集中在目標函數(shù)和約束函數(shù)上。如果每次計算都做雙層for循環(huán)掃描矩陣20節(jié)點可能還能忍50節(jié)點就非常吃力。改進方法有三個。第一核心矩陣操作全部向量化用矩陣運算替代逐元素for循環(huán)。第二用parfor對種群個體的適應度進行并行計算多核CPU直接跑滿。第三設(shè)置最大有效線路數(shù)量解碼時只統(tǒng)計前若干個最大元素其余直接視作零減少不必要的計算。5.5 結(jié)果復現(xiàn)性不好多次運行差異大差分進化算法本身帶有隨機性直接對比多輪結(jié)果時最優(yōu)值和目標函數(shù)都會出現(xiàn)波動。這屬于算法固有特性但工程上需要穩(wěn)定可復現(xiàn)的結(jié)果。建議在每次運行前固定隨機種子并記錄隨機數(shù)生成器的狀態(tài)。正式匯報時保留三個關(guān)鍵指標最優(yōu)解、多次運行的中位數(shù)解、最優(yōu)解對應的完整矩陣。同時采用“多次運行取最優(yōu)”的策略比如連續(xù)運行5次每次800代把最優(yōu)結(jié)果作為最終方案。我個人的體會是多微電網(wǎng)拓撲設(shè)計模型建模比算法調(diào)參重要得多。把工程需求翻譯成矩陣變量和約束條件這一步做踏實了剩下交給差分進化算法去搜基本能在合理時間內(nèi)給出工程上可用的方案。如果你想在這套代碼上繼續(xù)擴展可以考慮加入分布式協(xié)同進化把大矩陣按區(qū)域拆分后并行優(yōu)化再通過邊界協(xié)調(diào)收斂到全局一致解這會進一步提升大規(guī)模問題的求解能力。