簡介電力市場節(jié)點邊際電價出清優(yōu)化方案圍繞史新紅《機組運行約束對機組節(jié)點邊際電價的影響分析》進行復現(xiàn)以單時段出清模型為核心在忽略爬坡約束的前提下完整考慮機組運行約束適合電力市場、電力經(jīng)濟相關專業(yè)學生及研究者用于學習節(jié)點邊際電價LMP與系統(tǒng)邊際電價的形成機理。采用KKT條件將對偶問題轉(zhuǎn)化為矩陣形式求解拉格朗日乘子影子價格程序結構清晰、注釋完整便于對照報告理解影子價格在出清優(yōu)化中的經(jīng)濟含義。壓縮包共5個文件包含3個Matlab程序m、1篇參考文獻caj和1份完整報告doc總大小僅383KB輕量易下載。目前已有4668人學習資源內(nèi)附出清過程詳細報告并針對運行問題提供答疑支持可幫助讀者快速理解邊際電價出清的規(guī)范做法并復現(xiàn)算例。1. 打開這個zip之前先弄懂節(jié)點邊際電價為什么是一堆影子價格我做電力市場出清仿真這幾年最常被問的一句話是節(jié)點邊際電價LMP到底是怎么算出來的很多人第一反應是“用報價加權平均”或者“按負荷比例分攤成本”這兩個直覺在單節(jié)點系統(tǒng)里能湊合一旦把電網(wǎng)拓撲和線路阻塞放進來誤差能到幾十美元。真實情況是LMP不是市場規(guī)則里拍腦袋定的而是安全約束經(jīng)濟調(diào)度SCED這個線性規(guī)劃問題里功率平衡約束和線路阻塞約束的拉格朗日乘子組合出來的結果。換句話說價格不是算出來的是從優(yōu)化問題的對偶問題里“問”出來的。這個名為“基于yalmipcplex的電力市場-節(jié)點邊際電價出清優(yōu)化”的zip本質(zhì)上是把這條鏈路完整打通用YALMIP在MATLAB里把出清模型寫成約束和目標交給CPLEX求解再從求解結果里抽出每個節(jié)點的邊際電價。它適合三類人正在做電力市場仿真的研究生、售電或發(fā)電側(cè)做策略的算法工程師、想從通用優(yōu)化轉(zhuǎn)行電力調(diào)度的開發(fā)者。讀完你不僅能跑通一個三節(jié)點系統(tǒng)的出清還能自己拆LMP、判斷結果對不對、避開一堆讓人撓頭的坑。2. SCED與LMP直流潮流下的安全約束經(jīng)濟調(diào)度怎么定價2.1 沒有線路阻塞時電網(wǎng)其實是一個“大池子”先看最簡單的情況假設整個電網(wǎng)是一個點所有發(fā)電機和負荷接在同一個母線上不考慮線路容量。那么出清問題就是在滿足總負荷的前提下讓發(fā)電總成本最小。這是一個再簡單不過的線性規(guī)劃目標是最小化所有機組的發(fā)電成本之和約束只有一個——總出力等于總負荷。這個問題的對偶乘子拉格朗日乘子就是統(tǒng)一出清價所有節(jié)點完全一樣。為什么因為在這個假設下不管電從哪個機組來、送到哪個節(jié)點去成本只跟機組自身的報價有關跟“位置”無關。這就是電力市場里常說的“全網(wǎng)一口價”。但現(xiàn)實電網(wǎng)不是這樣。輸電線路有容量上限電能必須遵循基爾霍夫定律流動。于是同樣1MW負荷加在送端節(jié)點和加在受端節(jié)點引起的發(fā)電成本增量完全不同。這個“不同”就是節(jié)點邊際電價存在的原因。我一般會給新人打一個比方LMP是什么是“在這個節(jié)點多接1MW負荷整個系統(tǒng)最優(yōu)發(fā)電成本會增加多少錢”。這個增量包含了機組邊際報價也包含了線路阻塞帶來的機會成本。2.2 線路越過容量價格就開始“撕裂”用一個最小例子說明。兩個節(jié)點一條線路連著。送端節(jié)點1有一臺便宜機組報價10美元/MWh受端節(jié)點2有一臺貴機組報價20美元/MWh負荷全在節(jié)點2。線路容量很大時便宜機組直接送電全網(wǎng)價格10美元?，F(xiàn)在把線路容量調(diào)小比如只能送60MW而負荷是100MW。便宜機組最多往節(jié)點2送60MW剩下的40MW必須讓節(jié)點2的貴機組發(fā)。這時候全網(wǎng)不再是統(tǒng)一價。節(jié)點1的LMP是10美元——因為在節(jié)點1多1MW負荷本地便宜機組直接多發(fā)1MW就行。節(jié)點2的LMP變成20美元——因為在節(jié)點2多1MW負荷線路已經(jīng)滿了只能讓貴機組加出力。兩個節(jié)點的價差正好等于阻塞線路的“稀缺租金”。這個租金不是任何人主觀定的是線路潮流約束的對偶乘子。安全約束經(jīng)濟調(diào)度SCED的核心就是把這些約束全部放進優(yōu)化問題里讓價格自己“長出來”。2.3 為什么偏偏是YALMIP CPLEX理論上你自己用MATLAB寫一個單純形法或者調(diào)linprog也能解這個LP。但出清模型有個特殊要求你得知道每個約束對應的影子價格而且要把影子價格準確地映射到“功率平衡乘子”“線路阻塞乘子”這些業(yè)務概念上。自己寫求解器變量和約束多了以后對偶信息的跟蹤極其容易出錯。YALMIP的價值在于建模層和求解層分離。你只需要聲明變量是sdpvar然后像寫數(shù)學公式一樣寫約束YALMIP負責把模型轉(zhuǎn)換成標準LP/QP形式傳給底層求解器。CPLEX則在求解層提供了成熟的單純形和對偶單純形算法數(shù)值穩(wěn)定大規(guī)模稀疏問題也能扛。更關鍵的是YALMIP的dual函數(shù)能直接取回約束對應的對偶乘子這比從CPLEX原生接口里翻對偶解要友好得多。我見過不少直接從CPLEX Java/C接口取乘子的工程代碼里全是矩陣索引和符號轉(zhuǎn)換維護成本極高。另外YALMIP的好處是求解器可替換。同一個人樣模型把sdpsettings里的solver改成其他商業(yè)求解器或開源求解器也能跑只是性能和對偶信息獲取的方便程度不同。我一般在教學和驗證階段固定用CPLEX因為它的LP對偶解默認返回完整且預求解不會輕易“吃掉”你關注的約束乘子。3. YALMIP寫三節(jié)點出清建模、求解和第一個LMP3.1 直流潮流的近似邏輯與適用邊界市場出清仿真里最常用的潮流模型是直流潮流DCOPF不是交流潮流。交流潮流要同時算有功、無功、電壓幅值和相角問題是非凸的求解慢、可能不收斂而且LMP的定義在交流模型下要復雜得多——無功價格、電壓約束的影子價格都會摻進來。直流潮流做的事情很粗暴忽略無功、忽略網(wǎng)損、假設電壓幅值都是1只保留有功和相角的關系。這個關系寫出來很簡單線路有功潮流P_ij約等于兩端相角差除以電抗也就是P_ij(θ_i-θ_j)/x_ij。x是線路電抗單位用標幺值。這個線性關系讓整個出清問題變成一個純線性規(guī)劃凸、可收斂、對偶乘子有清晰的價格含義。代價是精度——直流潮流不適用于電壓問題、無功優(yōu)化場景但在節(jié)點邊際電價仿真、市場力分析、阻塞管理這些場景里它是絕對的主流選擇。你去看各地的現(xiàn)貨市場出清系統(tǒng)核心引擎基本都是某種形式的DCOPF或加入了網(wǎng)損近似的擴展版。需要注意一個邊界x_ij不能為0。電抗為零意味線路導納無窮大B矩陣奇異潮流方程無解。我見過有人在測試數(shù)據(jù)里為了省事把所有線路電抗寫成0結果YALMIP一跑就報“NaN”或者不可行最后查了半天發(fā)現(xiàn)是電抗單位寫錯了。3.2 三節(jié)點系統(tǒng)的參數(shù)與完整求解代碼為了把原理落到代碼上我構造一個最簡單但能體現(xiàn)阻塞定價的三節(jié)點鏈式系統(tǒng)。節(jié)點1是送端節(jié)點2是中轉(zhuǎn)節(jié)點3是負荷中心。節(jié)點1有一臺便宜機組節(jié)點2有一臺貴機組負荷100MW全部在節(jié)點3。線路1-2容量設為60MW線路2-3容量設為150MW。這樣一定會出現(xiàn)阻塞——便宜機組的電最多送60MW到節(jié)點2節(jié)點3的缺口只能讓節(jié)點2的貴機組頂上。系統(tǒng)參數(shù)如下對象參數(shù)數(shù)值線路1-2電抗(p.u.) / 容量(MW)0.1 / 60線路2-3電抗(p.u.) / 容量(MW)0.1 / 150機組G1所在節(jié)點 / 報價($/MWh) / 出力上下限(MW)1 / 10 / 0~100機組G2所在節(jié)點 / 報價($/MWh) / 出力上下限(MW)2 / 20 / 0~100負荷D3所在節(jié)點 / 大小(MW)3 / 100下面是完整的YALMIP CPLEX求解代碼。這段代碼我建議你原樣跑通一次再逐步改參數(shù)觀察價格變化。%% 三節(jié)點直流潮流出清YALMIP CPLEX % 節(jié)點1送端2中轉(zhuǎn)3受端(負荷) % 線路L12連接1-2L23連接2-3 clear; clc; %% 系統(tǒng)數(shù)據(jù) mpc.nb 3; % 節(jié)點數(shù) mpc.nl 2; % 線路數(shù) % 線路: [起點, 終點, 電抗(p.u.), 容量(MW)] mpc.line [1 2 0.1 60; 2 3 0.1 150]; % 機組: [節(jié)點, 報價($/MWh), 出力下限(MW), 出力上限(MW)] mpc.gen [1 10 0 100; 2 20 0 100]; mpc.load [0 0 100]; % 節(jié)點3負荷100MW %% 構建直流潮流B矩陣和線路-節(jié)點關聯(lián)矩陣Bf B zeros(mpc.nb, mpc.nb); for k 1:mpc.nl f mpc.line(k,1); t mpc.line(k,2); x mpc.line(k,3); B(f,f) B(f,f) 1/x; B(f,t) B(f,t) - 1/x; B(t,f) B(t,f) - 1/x; B(t,t) B(t,t) 1/x; end Bf zeros(mpc.nl, mpc.nb); for k 1:mpc.nl f mpc.line(k,1); t mpc.line(k,2); x mpc.line(k,3); Bf(k,f) 1/x; Bf(k,t) -1/x; end ref 3; % 參考節(jié)點選3相角置0 %% 定義優(yōu)化變量 P sdpvar(2, 1); % 兩臺機組出力 theta sdpvar(3, 1); % 節(jié)點相角 %% 約束 cons_plim (mpc.gen(:,3) P mpc.gen(:,4)); % 機組出力上下限 cons_ref (theta(ref) 0); % 參考節(jié)點相角為0 F Bf * theta; % 線路有功潮流表達式 cons_fline (F mpc.line(:,4)); % 線路潮流上限 cons_fdown (-mpc.line(:,4) F); % 線路潮流下限 cons_balance (sum(P) sum(mpc.load)); % 功率平衡 % 合并所有約束傳給求解器 all_cons [cons_plim, cons_ref, cons_fline, cons_fdown, cons_balance]; %% 目標函數(shù)與求解 obj mpc.gen(:,2) * P; % 總發(fā)電成本線性成本 options sdpsettings(solver, cplex, verbose, 1); sol optimize(all_cons, obj, options); %% 讀取結果 P_opt value(P); F_opt value(F); theta_opt value(theta); fprintf(G1出力: %.2f MW, G2出力: %.2f MW\n, P_opt(1), P_opt(2)); fprintf(線路L12潮流: %.2f MW, 線路L23潮流: %.2f MW\n, F_opt(1), F_opt(2)); fprintf(總發(fā)電成本: %.2f 美元\n, value(obj));邏輯上說得很直白先構建B矩陣和Bf矩陣B矩陣用來算節(jié)點相角和功率注入的關系Bf矩陣用來從相角直接得到每條線路的潮流。然后聲明兩臺機組的出力P和三個節(jié)點的相角theta為優(yōu)化變量。約束分四組機組出力上下限、參考節(jié)點相角、線路潮流上下限、功率平衡。目標函數(shù)是線性成本之和。參數(shù)說明要強調(diào)幾個地方。第一參考節(jié)點必須固定否則所有相角同時加一個常數(shù)潮流不變但矩陣求解不唯一YALMIP不會報錯但結果看著像隨機數(shù)。第二F Bf * theta不是約束是表達式它把相角和潮流聯(lián)系起來后續(xù)的線路容量約束都作用在這個表達式的值上。第三sdpsettings里solver指定用cplexverbose1讓求解器打印迭代日志如果CPLEX沒有正確安裝或授權YALMIP會報錯這時可以把solver改成其他支持的求解器先驗證模型邏輯。跑完這段代碼你會得到G1出力60MW、G2出力40MW、總成本1400美元。這個結果本身就是阻塞定價的起點——便宜機組只發(fā)到線路容量允許的上限剩下的全靠貴機組。3.3 第一次看結果哪些約束在“卡脖子”求解完成后不要只看出力數(shù)字要看哪條約束是積極約束active constraint也就是dual值不等于0的約束。積極約束就是“卡脖子”的位置線路L12的潮流上限約束和功率平衡約束都會有非零對偶乘子。機組G1的出力上限60MW不會觸發(fā)它的上限是100MW觸發(fā)的是線路容量60MW。我習慣把結果分成兩種場景做對比。場景A線路L12容量60MW價格分層場景B把線路L12容量改成150MW也就是不阻塞。兩張場景的LMP對比如下節(jié)點場景AL12容量60MW場景BL12容量150MW節(jié)點110 美元/MWh10 美元/MWh節(jié)點220 美元/MWh10 美元/MWh節(jié)點320 美元/MWh10 美元/MWh機組出力G160, G240G1100, G20場景B里全網(wǎng)一口價10美元跟2.1節(jié)的理論預測完全一致。場景A里節(jié)點1和節(jié)點2/3差了10美元差的這10美元就是線路L12阻塞的影子價格。這個“一堵就分層”的現(xiàn)象是所有節(jié)點邊際電價仿真的核心。后面第四章我要講的LMP分解就是把這10美元從數(shù)學上拆出來。4. LMP的提取與分解從對偶乘子到節(jié)點電價4.1 取對偶乘子YALMIP的dual()與符號約定YALMIP里取對偶乘子非常簡單直接調(diào)用dual(約束句柄)。但要命的是符號約定對于形如“表達式 0”的約束YALMIP返回的乘子是非負的對于等式約束乘子可以是任意實數(shù)符號由等式方向決定。我初學時在這里栽過跟頭——把功率平衡約束的dual取出來直接當成LMP結果算出來的價格全反了。標準做法是把功率平衡等式約束寫成sum(P) sum(mpc.load)然后dual(cons_balance)就是能量價格基礎分量λ。線路潮流約束F Fmax返回的是非負乘子μ_up-Fmax F返回的也是非負乘子μ_dn。不需要手動取負號只需要在匯總公式里按KKT條件的方向組合。4.2 拆LMP能量分量、阻塞分量與參考節(jié)點有了對偶乘子還需要一個矩陣節(jié)點轉(zhuǎn)移因子PTDF。PTDF的含義是某個節(jié)點注入1MW功率參考節(jié)點吸收每條線路上會增加多少潮流。直流潮流下PTDF矩陣可以由Bf和B矩陣直接算出來。下面這段代碼接著上一節(jié)的系統(tǒng)數(shù)據(jù)計算PTDF并把三個節(jié)點的LMP拆成能量分量和阻塞分量。%% 計算PTDF矩陣 nonref [1 2]; % 非參考節(jié)點參考節(jié)點3的列保持0 Bred B(nonref, nonref); % 去掉參考節(jié)點行列的B矩陣 PTDF zeros(mpc.nl, mpc.nb); PTDF(:, nonref) Bf(:, nonref) / Bred; % 線路注入靈敏度參考節(jié)點列為0 %% 取對偶乘子 lambda dual(cons_balance); % 功率平衡乘子能量價格基礎 mu_up dual(cons_fline); % 線路FFmax的乘子非負 mu_dn dual(cons_fdown); % 線路-FmaxF的乘子非負 %% 組合LMP LMP zeros(mpc.nb, 1); for i 1:mpc.nb LMP(i) lambda - sum((mu_up - mu_dn) .* PTDF(:, i)); end %% 打印結果 fprintf(能量分量(全網(wǎng)統(tǒng)一): %.2f $/MWh\n, lambda); for i 1:mpc.nb blocking lambda - LMP(i); % 阻塞分量含符號 fprintf(節(jié)點%d LMP%.2f $/MWh, 阻塞分量%.2f $/MWh\n, ... i, LMP(i), blocking); end這里PTDF(:,i)是一個列向量表示節(jié)點i的注入對各條線路潮流的靈敏度。用mu_up減去mu_dn再和PTDF逐元素相乘后求和得到該節(jié)點阻塞分量的修正值。公式為什么是“減去”而不是“加上”回到第二章的雙節(jié)點例子送端節(jié)點1的LMP10全網(wǎng)能量分量為20阻塞修正值是-10方向必須能讓送端更便宜、受端更貴。如果你拿到的數(shù)據(jù)算出來方向相反優(yōu)先檢查PTDF的參考節(jié)點是否選對了以及線路約束寫成“F Fmax”時YALMIP返回的dual是否真的非負。在三節(jié)點例子里這個代碼會輸出lambda20美元節(jié)點1的LMP10美元阻塞分量-10節(jié)點2和節(jié)點3的LMP20美元阻塞分量0。這跟3.3節(jié)表格完全一致。阻塞分量只作用于阻塞線路的“送端”這是直流潮流的典型特征。注意阻塞盈余的概念也由此而來節(jié)點3的用戶按20美元買電節(jié)點1的機組按10美元賣電中間每度電的10美元差額就是阻塞盈余在真實市場里會通過阻塞費用結算機制分配。關于參考節(jié)點再補充一句LMP的數(shù)值總量不依賴參考節(jié)點的選擇三個節(jié)點算出來都一致。但分解出來的能量分量和阻塞分量會隨參考節(jié)點變化——參考節(jié)點的阻塞分量恒為0能量分量會被“吸收”一部分。所以做分解時不要在報告里只說“能量價是多少”一定要注明參考節(jié)點是哪個。4.3 擾動校核用電價的定義反推你的結果這是我最推薦的驗證習慣。LMP的嚴格定義是“某節(jié)點增加單位負荷引起的最優(yōu)成本增量”。與其信公式不如直接讓計算機“擾動”一下。做法是把求解過程封裝成一個函數(shù)然后對每個節(jié)點做一次1MW負荷擾動重新求解比較成本差。function [P, F, obj_val] solve_sced(mpc) % 把第三章的建模和求解包成函數(shù)返回出力、潮流和目標值 % mpc結構里有l(wèi)ine/gen/load/nb/nl等字段 % 具體建模代碼同第三章略 end %% 擾動校核對每個節(jié)點加1MW負荷求成本增量 c0 0; % 先求基礎場景的總成本 [~, ~, c0] solve_sced(mpc); for i 1:mpc.nb mpc2 mpc; mpc2.load(i) mpc2.load(i) 1; % 只在節(jié)點i加1MW [~, ~, c1] solve_sced(mpc2); fprintf(節(jié)點%d擾動法LMP%.2f $/MWh\n, i, c1 - c0); end在三節(jié)點例子里這個擾動校核的結果會和4.2節(jié)的LMP完全吻合節(jié)點1加1MW負荷成本增加10美元節(jié)點2和節(jié)點3加1MW負荷成本增加20美元。邏輯也很好驗證——節(jié)點1的負荷增加不會加重線路L12的負擔本地負荷吸收本地出力便宜機組直接加發(fā)節(jié)點2和節(jié)點3的負荷增加則必須靠貴機組。擾動法看著笨但它是唯一不考慮對偶理論也能解釋得通的價格驗證手段。我每次寫新模型都會保留這個校核腳本作為LMP結果的“金標準”。5. 出清仿真的五類翻車現(xiàn)場與排查路徑5.1 INFEASIBLE電抗為零、參考缺失與負荷不平衡現(xiàn)象optimize返回求解狀態(tài)為不可行或者YALMIP直接報NaN。原因有三大類按我實際遇到頻次排序第一線路電抗x寫成0B矩陣奇異第二參考節(jié)點相角約束沒有加模型自由度過大導致數(shù)值病態(tài)第三總發(fā)電出力上限小于總負荷數(shù)學上無解。還有一類隱蔽的負荷數(shù)據(jù)寫在某個節(jié)點但該節(jié)點沒有發(fā)電機也沒有外來線路形成孤島。解決先查B矩陣的行列式行列式接近0立刻定位電抗異常檢查所有節(jié)點是否至少有發(fā)電機或線路連通到參考節(jié)點?？焖俣ㄎ徊豢尚屑s束的辦法是給功率平衡加一個人工松弛變量比如把平衡約束改成sum(P)slacksum(load)讓slack無界求解后看slack落在哪個節(jié)點就知道哪里功率缺了。5.2 dual()全為零對偶信息被預求解消掉了現(xiàn)象模型求解成功出力、潮流都對但dual(cons_balance)返回0。LMP算出來全是0怎么看都不對。原因CPLEX預求解階段可能判定某條約束冗余并直接刪掉對偶乘子隨之丟失。尤其當功率平衡約束的表達式寫法不夠“標準”時容易觸發(fā)。另外YALMIP里多個約束合并成復合約束后dual返回的向量順序和約束順序如果不匹配也容易讓人誤讀成零。解決把需要取乘子的約束單獨定義成獨立的句柄不要和其他約束揉在一起。求解前在sdpsettings里關閉部分預求解options.cplex.preprocessing.reduce 0。如果還是取不到用擾動法4.3節(jié)交叉驗證——不需要信dual也能得到價格。這是最保險的兜底手段。5.3 全網(wǎng)一口價是不是哪里沒設容量現(xiàn)象明明加了線路容量約束三節(jié)點系統(tǒng)還是全網(wǎng)統(tǒng)一價跟第二章“大池子”一樣。原因線路容量設得太大了。容量上限在200MW線路實際潮流才100MW約束雖然存在但從未被激活乘子當然為0價格不會分層。很多人第一次跑通代碼后發(fā)現(xiàn)價格都一樣以為自己優(yōu)化錯了其實是數(shù)據(jù)沒卡到阻塞點。解決逐步縮小關鍵線路容量觀察LMP分叉的臨界點。這個臨界值其實有物理含義線路容量緊到一定程度市場上的便宜機組送不過去受端價格開始爬升。我在3.3節(jié)里把L12容量從150降到60就是故意制造這個分叉。調(diào)參時每次都對比“潮流值vs容量上限”兩個數(shù)字相等的那條線就是當前最擁堵的線。5.4 LMP是負的或高于機組報價稀缺定價不是bug現(xiàn)象某節(jié)點LMP出現(xiàn)負值或者高于全網(wǎng)最高機組報價20美元達到35美元。原因這是線性規(guī)劃里約束影子價格的正常表現(xiàn)不是求解器抽風。負價格通常出現(xiàn)在機組有最小技術出力比如最小出力50MW但負荷只有30MW或者必須開機維持電壓支撐的場景價格高于最高報價則說明某條受約束的線路把稀缺性放大受端愿意付出的邊際代價遠超本地最貴機組的報價。解決不要“修”結果先理解是哪條約束在定價。把該節(jié)點的擾動法LMP手動算一遍你會發(fā)現(xiàn)它和影子價格完全一致。如果業(yè)務上不能接受負價格那是市場規(guī)則需要加價格下限不是優(yōu)化模型的問題。仿真時保留這個現(xiàn)象即可。5.5 多解與退化同一份數(shù)據(jù)不同的價格現(xiàn)象代碼沒變、數(shù)據(jù)沒變不同機器上跑出來G1/G2出力的分配有差異甚至LMP在某些節(jié)點上有微小波動。原因LP存在多個最優(yōu)解尤其是機組報價都是線性沒有二次項時成本曲線是一條直線多個出力組合可以達到相同的總成本。單純形法會沿著頂點走起點和預求解差異會讓它落在不同的最優(yōu)頂點上。價格本身不一定變但個別節(jié)點的阻塞分量可能因為參考約束的取舍而不同。解決給目標函數(shù)加一個極小的二次項比如0.0001×P^2讓問題變成凸二次規(guī)劃最優(yōu)解唯一。這在實際現(xiàn)貨市場出清里也是常用手段——報價段內(nèi)部按線性處理但出清程序里會加微小的平滑項避免解漂移。代價是價格含義從純線性變成帶微弱二次修正但量級完全可以忽略。6. 從單時段到多時段SCED擴展成時序價格的進階技巧6.1 給模型加時間軸爬坡約束與24小時價格單時段模型解決好“某一時刻的電價”但現(xiàn)貨市場價格是時序的調(diào)度員關心的是機組能不能在15分鐘或1小時內(nèi)從50MW爬升到80MW。給模型加時間軸的套路是固定的所有變量加一個時間維度P變成2×T的矩陣theta變成3×T的矩陣功率平衡、線路容量約束對每個時段單獨寫目標函數(shù)變成所有時段的成本之和。額外要加的約束是爬坡約束T 24; % 24個時段 P sdpvar(2, T); % 每臺機組每個時段的出力 R_up 30; % 每時段最大爬坡上限MW R_dn 30; % 每時段最大爬坡下限MW cons_ramp []; for t 2:T cons_ramp [cons_ramp, P(:,t) - P(:,t-1) R_up]; cons_ramp [cons_ramp, P(:,t-1) - P(:,t) R_dn]; end這里P(:,t)-P(:,t-1)用矩陣形式一次性約束兩臺機組比寫兩個標量約束干凈。爬坡約束會對價格產(chǎn)生實質(zhì)影響低谷時段便宜機組被強制壓低出力高峰時段又爬不上去受端節(jié)點在爬坡受限的時段會看到價格尖峰。每個時段的功率平衡約束dual取出來就是該時段的能量價格序列把它們連成曲線就是典型的24小時出清價。注意多時段模型里LMP的取值要區(qū)分“該時段的lambda”和“跨時段爬坡約束的影子價格”后者會影響多個時段的價格形成分析時不要混在一起。6.2 MIP與LP的兩階段出清機組啟停后怎么取LMP這是最容易被忽視的進階問題。一旦引入機組啟停變量變成0/1整數(shù)模型從LP變成MIP。從數(shù)學上嚴格說MIP的對偶乘子不再有經(jīng)典LP影子價格的解釋直接從混合整數(shù)解里取dual得到的“價格”是失真且沒有經(jīng)濟含義的。行業(yè)里標準的做法是兩階段出清第一階段用MIP求解機組啟停和出力計劃目標函數(shù)包含啟停成本第二階段把第一步解出來的0/1狀態(tài)變量固定成常數(shù)重新求解一次純LP從第二次LP里取對偶乘子作為LMP。YALMIP里做這個操作也很方便——binvar聲明啟停變量第一階段optimize之后把啟停變量的值用value提取出來替換成約束值再二次optimize。第二階段模型不包含整數(shù)變量CPLEX自動走LP路徑dual()就能正常取回乘子。我去年在某模擬項目X里第一次做帶啟停的出清時直接對MIP結果取dual算出來的價格低得離譜后來查文獻才發(fā)現(xiàn)這個坑。從此固定了“先MIP定狀態(tài)、再LP定價”的習慣。6.3 校驗腳本與工程習慣最后說一個長期受益的習慣每次建新模型第一件事不是跑結果而是寫一個自動校核腳本包含三類檢查。一是潮流檢查把解出來的P和theta代回Bf*theta確認潮流方程成立。二是平衡檢查每時段的發(fā)電機總出力減去負荷誤差不超過1e-6。三是價格檢查用4.3節(jié)的擾動法隨機抽三個節(jié)點驗證LMP。這三關過了再談報告和曲線。這套習慣幫我擋掉了無數(shù)個“看起來很美但實際算錯”的版本。說到底出清優(yōu)化最玄學的地方不在YALMIP語法也不在CPLEX配置而在“你有沒有真正理解什么是影子價格”。我踩過的坑從符號搞反到預求解吞乘子沒有一個能靠調(diào)參解決全是回到定義本身才找到答案。希望這篇筆記能幫你少走這些彎路把時間花在分析和業(yè)務上。本文還有配套的精品資源點擊獲取