簡(jiǎn)介這是一份面向電力系統(tǒng)與儲(chǔ)能領(lǐng)域研究者的 MATLAB 開(kāi)源代碼聚焦電價(jià)套利背景下儲(chǔ)能系統(tǒng)的最優(yōu)運(yùn)行優(yōu)化問(wèn)題。通過(guò)建立以經(jīng)濟(jì)效益最大化為目標(biāo)函數(shù)的優(yōu)化框架并考慮充放電約束、電價(jià)波動(dòng)等關(guān)鍵條件幫助用戶在分時(shí)電價(jià)機(jī)制下制定合理的儲(chǔ)能調(diào)度策略。代碼基于兩篇國(guó)外相關(guān)文獻(xiàn)的方法提供了三種獨(dú)立可運(yùn)行的求解路線簡(jiǎn)單直接的運(yùn)行策略搜索、能夠處理隨機(jī)性的蒙特卡洛模擬法以及基于 fmincon 函數(shù)的數(shù)學(xué)規(guī)劃優(yōu)化方法。三種方法相互對(duì)照便于理解不同思路在收斂性、求解速度與結(jié)果精度上的差異也方便使用者按需替換或擴(kuò)展。壓縮包僅約114KB共10個(gè)文件以 MATLAB 的 .m 腳本為主5個(gè)同時(shí)包含電價(jià)時(shí)間序列文本、spot_price_2013.mat 數(shù)據(jù)、README 說(shuō)明文檔、許可證等目錄精簡(jiǎn)、結(jié)構(gòu)清晰下載后即可按腳本逐步運(yùn)行。目前已有973人學(xué)習(xí)下載適合具有一定 MATLAB 基礎(chǔ)、關(guān)注儲(chǔ)能經(jīng)濟(jì)調(diào)度或電力市場(chǎng)套利的入門與進(jìn)階學(xué)習(xí)者用于算法復(fù)現(xiàn)、對(duì)比實(shí)驗(yàn)或課題參考。1. 考慮經(jīng)濟(jì)性的儲(chǔ)能運(yùn)行優(yōu)化這套開(kāi)源 MATLAB 代碼解決什么問(wèn)題把「考慮經(jīng)濟(jì)性的儲(chǔ)能運(yùn)行優(yōu)化」這套開(kāi)源 MATLAB 代碼跑通你會(huì)得到一個(gè)可以直接復(fù)用的調(diào)度決策模板。它的典型場(chǎng)景是某園區(qū)裝了 200 kWh 儲(chǔ)能柜分時(shí)電價(jià)峰谷價(jià)差 0.9 元/kWh看著每天倒騰一次就能賺一百多但真跑起來(lái)就發(fā)現(xiàn)沒(méi)那么簡(jiǎn)單——什么時(shí)候充、什么時(shí)候放、充放多快、要不要留容量給下一個(gè)峰段每個(gè)決定都直接影響最終收益而充放次數(shù)多了還會(huì)加速電池老化。這套代碼做的事就是把上述決策寫成一個(gè)線性規(guī)劃問(wèn)題用 MATLAB 內(nèi)置的 linprog 求解。輸入只要三樣24 小時(shí)分時(shí)電價(jià)、儲(chǔ)能電氣參數(shù)容量、功率、效率、SOC 上下限輸出是一天逐小時(shí)的充放電功率和 SOC 曲線目標(biāo)是在全部運(yùn)行約束下讓凈收益最大。它不需要額外安裝求解器有 MATLAB 基礎(chǔ)就能改能跑。適合三類人做儲(chǔ)能調(diào)度研究的學(xué)生拿它當(dāng)基準(zhǔn)結(jié)果做對(duì)比評(píng)估儲(chǔ)能項(xiàng)目收益的工程師用它算不同價(jià)差下的理論收益上限以及想學(xué)優(yōu)化建模、需要一份干凈模板的開(kāi)發(fā)者。下面從建模開(kāi)始一步步把來(lái)龍去脈講清楚。2. 儲(chǔ)能運(yùn)行優(yōu)化建模目標(biāo)函數(shù)、約束條件與求解器選型先講透優(yōu)化模型動(dòng)手寫代碼之前先花十分鐘把數(shù)學(xué)形式寫清楚這比調(diào)試省時(shí)間得多。儲(chǔ)能運(yùn)行優(yōu)化的標(biāo)準(zhǔn)形式是給定一個(gè)調(diào)度周期通常 24 小時(shí)、步長(zhǎng) 1 小時(shí)在滿足電池自身運(yùn)行約束的前提下最大化周期內(nèi)的凈收益。核心就三件事目標(biāo)函數(shù)里放什么、約束條件有哪些、用哪類求解器。2.1 目標(biāo)函數(shù)怎么定峰谷套利、老化成本與需量電費(fèi)的取舍最基礎(chǔ)的目標(biāo)函數(shù)是峰谷套利。每個(gè)時(shí)段的收益等于放電收入減去充電成本寫成max Σ ( p_t * q_dis,t - p_t * q_ch,t ) * dt其中 p_t 是該時(shí)段電價(jià)q_dis,t 是放電功率q_ch,t 是充電功率dt 是步長(zhǎng)。這個(gè)形式只反映了電費(fèi)差忽略了一個(gè)關(guān)鍵問(wèn)題電池每放一度電循環(huán)壽命都在縮短這部分損耗是真實(shí)成本。所以第二個(gè)層次是把老化成本折算進(jìn)目標(biāo)函數(shù)放電收益變成 (p_t - c_deg) * q_dis,t * dt其中 c_deg 是每放電 1 kWh 折算出的電池?fù)p耗費(fèi)用具體折算方法見(jiàn)第 4 章。第三個(gè)層次是需量電費(fèi)管理。國(guó)內(nèi)工商業(yè)電價(jià)常見(jiàn)兩部制除了電度電費(fèi)還有按最大需量計(jì)收的基本電費(fèi)。儲(chǔ)能削減峰值負(fù)荷能直接省下這部分錢目標(biāo)函數(shù)里就要加入「需量電費(fèi)單價(jià) × 削減的峰值功率」這一項(xiàng)。但注意這一項(xiàng)會(huì)讓模型從純線性變成帶 max 算子的分段線性問(wèn)題通常需要用 epigraph 變換或引入輔助變量才能保持線性。對(duì)絕大多數(shù)入門場(chǎng)景來(lái)說(shuō)先做前兩個(gè)層次就夠了套利加老化成本。第三個(gè)層次建議在基礎(chǔ)版本跑通后再擴(kuò)展否則一邊學(xué)模型一邊處理 max 算子的重構(gòu)很容易卡住。2.2 約束條件怎么寫SOC 遞推、功率邊界與充放電互斥約束條件里最重要的一條是 SOC 遞推方程它描述電池能量隨時(shí)間的變化SOC(t1) SOC(t) [ eta_ch * q_ch,t - q_dis,t / eta_dis ] * dt / E_ratedeta_ch 和 eta_dis 分別是充放電效率E_rated 是額定容量。這條約束把相鄰時(shí)段的 SOC 串起來(lái)是整個(gè)模型的骨架。除了遞推方程還有幾類邊界約束SOC 上下限比如 0.1 到 0.9防止過(guò)充過(guò)放、充放電功率上下限0 到 P_max、初始 SOC 固定、以及終端 SOC 約束——要求周期結(jié)束時(shí)回到初始值避免優(yōu)化器把初始存電當(dāng)免費(fèi)資源一次性放空。還有一個(gè)容易糾結(jié)的點(diǎn)充放電互斥約束 q_ch * q_dis 0。這是非線性約束但絕大多數(shù)情況下可以安全丟棄。原因是只要充電和放電發(fā)生在同一時(shí)段能量先充進(jìn)去再放出來(lái)要經(jīng)過(guò)一次充放效率損耗0.95 × 0.95 ≈ 0.9白白損失 10% 能量電價(jià)是正的這種自我抵消的行為永遠(yuǎn)不可能是最優(yōu)解。省掉這個(gè)約束問(wèn)題就從非線性規(guī)劃退化成純線性規(guī)劃求解速度和穩(wěn)定性都大幅提升。只有當(dāng)模型里加了最小充放電功率、啟停成本這類條款時(shí)才需要引入 0-1 二進(jìn)制變量重新變成 MILP。2.3 為什么用線性規(guī)劃linprog 與 YALMIP/啟發(fā)式算法的邊界儲(chǔ)能價(jià)格套利問(wèn)題本質(zhì)上是 LP目標(biāo)線性、約束線性、變量連續(xù)。線性規(guī)劃的優(yōu)點(diǎn)是能保證全局最優(yōu)、求解可復(fù)現(xiàn)、速度快不存在啟發(fā)式算法那種每次運(yùn)行結(jié)果不穩(wěn)定、需要反復(fù)調(diào)參數(shù)的玄學(xué)。常見(jiàn)做法是直接用 MATLAB 內(nèi)置 linprog或者用 YALMIP 搭模型再調(diào)用底層求解器。YALMIP 的好處是建模語(yǔ)法接近數(shù)學(xué)表達(dá)式后續(xù)要加二進(jìn)制變量、做 MILP 擴(kuò)展時(shí)不用重寫整套矩陣。三者的選擇邊界大致如下求解方案能處理的模型獲取難度適用場(chǎng)景MATLAB linprog連續(xù)變量 LP內(nèi)置 Optimization Toolbox純價(jià)格套利、無(wú)二進(jìn)制變量的場(chǎng)景YALMIP GLPK/SCSLP/MILP免費(fèi)求解器需下載并配置路徑學(xué)生做原型、驗(yàn)證算法YALMIP Gurobi/CPLEX大規(guī)模 LP/MILP商業(yè)授權(quán)或?qū)W術(shù)授權(quán)高維問(wèn)題、強(qiáng)耦合約束經(jīng)驗(yàn)是先別急著上 YALMIP 和商業(yè)求解器。用 linprog 把純 LP 版本跑通確認(rèn)模型邏輯沒(méi)問(wèn)題再按需升級(jí)。很多人一上來(lái)就用遺傳算法之類的啟發(fā)式方法反而把簡(jiǎn)單問(wèn)題復(fù)雜化——非線性效率曲線、滯回特性這類模型才值得動(dòng)用啟發(fā)式線性能解的問(wèn)題用線性方法拿全局最優(yōu)解才是正路。3. 用 MATLAB 跑通儲(chǔ)能經(jīng)濟(jì)性優(yōu)化linprog 建模代碼逐段拆解模型立住之后代碼就是把數(shù)學(xué)翻譯成矩陣。下面用一組常見(jiàn)參數(shù)做演示200 kWh / 50 kW 儲(chǔ)能系統(tǒng)充放電效率 0.95SOC 范圍 0.10.9三段式分時(shí)電價(jià)。代碼按數(shù)據(jù)準(zhǔn)備、模型構(gòu)建、結(jié)果處理三塊走。3.1 參數(shù)與分時(shí)電價(jià)數(shù)據(jù)準(zhǔn)備一份可直接改的輸入模板%% 儲(chǔ)能參數(shù)與分時(shí)電價(jià)定義 clc; clear; close all; E_rated 200; % 額定容量/kWh P_max 50; % 最大充放電功率/kW SOC_init 0.5; % 初始荷電狀態(tài) SOC_min 0.1; % SOC 下限 SOC_max 0.9; % SOC 上限 eta_ch 0.95; % 充電效率 eta_dis 0.95; % 放電效率 dt 1; % 時(shí)間步長(zhǎng)/h % 24小時(shí)分時(shí)電價(jià)單位元/kWh % 0-7點(diǎn)谷段0.38-17點(diǎn)平段0.818-21點(diǎn)峰段1.222-23點(diǎn)平段0.8 price [0.3*ones(1,8), 0.8*ones(1,10), 1.2*ones(1,4), 0.8*ones(1,2)]; N 24; % 優(yōu)化周期內(nèi)時(shí)段數(shù)這里幾個(gè)要點(diǎn)price 向量必須按時(shí)間順序排列且長(zhǎng)度等于 N效率值用的是小數(shù)而不是百分?jǐn)?shù)dt 的單位要和電價(jià)單位匹配電價(jià)是元/kWh、功率是 kW乘上 dt 小時(shí)得到元。如果要做多日連續(xù)仿真把 price 向量按天拼接即可但要注意 SOC 會(huì)在日與日之間傳遞不在本章范圍第 6 章會(huì)講滾動(dòng)處理。3.2 核心建模與求解代碼決策變量、約束矩陣與 linprog 調(diào)用%% 決策變量布局 % x [q_ch(1..N), q_dis(1..N), SOC(1..N1)] % SOC(1)是第1時(shí)段開(kāi)始時(shí)SOCSOC(N1)是周期結(jié)束時(shí)SOC n_var 3*N 1; % 目標(biāo)函數(shù)linprog 求最小值所以充電成本為正、放電收益取負(fù) f zeros(n_var, 1); f(1:N) price * dt; % 充電花費(fèi) f(N1:2*N) -price * dt; % 放電收入取負(fù)號(hào) %% 等式約束SOC 遞推 % SOC(t1) SOC(t) eta_ch*q_ch(t)*dt/E_rated - q_dis(t)*dt/(E_rated*eta_dis) Aeq zeros(N, n_var); beq zeros(N, 1); for t 1:N Aeq(t,t) eta_ch * dt / E_rated; % q_ch(t) 系數(shù) Aeq(t,Nt) -dt / (E_rated * eta_dis); % q_dis(t) 系數(shù) Aeq(t,2*Nt) -1; % -SOC(t) Aeq(t,2*Nt1) 1; % SOC(t1) end % 初始 SOC 固定為 SOC_init Aeq [Aeq; zeros(1, n_var)]; Aeq(end, 2*N1) 1; beq [beq; SOC_init]; % 終端 SOC 約束周期結(jié)束回到初始值 Aeq [Aeq; zeros(1, n_var)]; Aeq(end, 3*N1) 1; beq [beq; SOC_init]; %% 變量邊界 lb zeros(n_var, 1); ub inf(n_var, 1); lb(1:2*N) 0; ub(1:2*N) P_max; % 充放電功率邊界 lb(2*N1:end) SOC_min; ub(2*N1:end) SOC_max; % SOC 邊界 lb(2*N1) SOC_init; ub(2*N1) SOC_init; % 初始 SOC 雙邊界固定 %% 求解 options optimoptions(linprog,Algorithm,dual-simplex,Display,iter); [x_opt, fval, exitflag] linprog(f, [], [], Aeq, beq, lb, ub, options);逐段說(shuō)明一下。決策變量把充電功率、放電功率、SOC 排成一個(gè)長(zhǎng)向量長(zhǎng)度是 3N1。SOC 多一個(gè)點(diǎn)是用來(lái)表示周期結(jié)束狀態(tài)的這個(gè)布局決定了后面所有矩陣的列索引改的時(shí)候要非常小心。目標(biāo)函數(shù)里充電對(duì)應(yīng)正系數(shù)因?yàn)?linprog 是最小化放電對(duì)應(yīng)負(fù)系數(shù)最小化負(fù)的收益等價(jià)于最大化收益。SOC 變量的系數(shù)是 0它本身沒(méi)有直接成本。SOC 遞推約束用循環(huán)逐時(shí)段生成每個(gè)時(shí)段一行等式系數(shù)來(lái)自遞推公式的移項(xiàng)。這里最容易犯的錯(cuò)是 eta_dis 寫錯(cuò)位置——它是放在放電功率的分母上不是乘在分子上寫反了結(jié)果會(huì)差一個(gè)效率的平方。最后兩個(gè)等式分別固定初始 SOC 和終端 SOC。終端 SOC 約束是防止優(yōu)化器把初始電量當(dāng)免費(fèi)午餐的后悔藥如果去掉典型結(jié)果會(huì)是最后一個(gè)峰段把電池放到 SOC 下限第二天無(wú)電可用。提示初次調(diào)試時(shí)把Display設(shè)為iter能看到單純形迭代過(guò)程確認(rèn)模型穩(wěn)定后再改回final減少刷屏。3.3 結(jié)果提取與繪圖充放電曲線、SOC 曲線與收益統(tǒng)計(jì)%% 結(jié)果提取 q_ch x_opt(1:N); % 各時(shí)段充電功率/kW q_dis x_opt(N1:2*N); % 各時(shí)段放電功率/kW SOC x_opt(2*N1:end); % SOC序列長(zhǎng)度 N1 profit -fval; % 凈收益 -目標(biāo)函數(shù)值 E_ch sum(q_ch) * dt; % 總充電電量/kWh E_dis sum(q_dis) * dt; % 總放電電量/kWh fprintf(凈收益: %.2f 元\n, profit); fprintf(充電電量: %.2f kWh, 放電電量: %.2f kWh\n, E_ch, E_dis); %% 繪圖功率、SOC、電價(jià)三聯(lián)圖 t 1:N; figure(Color,w,Position,[100 100 900 650]); subplot(3,1,1); bar(t, q_ch, 0.6, FaceColor,[0.2 0.5 1]); hold on; bar(t, -q_dis, 0.6, FaceColor,[1 0.6 0.2]); legend(充電,放電,Location,best); ylabel(功率/kW); ylim([-1.2*P_max 1.2*P_max]); xlim([0.5 N0.5]); subplot(3,1,2); stairs([t, N1], SOC, LineWidth,1.5); ylabel(SOC); ylim([0 1]); xlim([0.5 N0.5]); grid on; subplot(3,1,3); stairs(t, price, k, LineWidth,1.5); ylabel(電價(jià)/(元/kWh)); xlabel(時(shí)間/h); xlim([0.5 N0.5]); grid on;結(jié)果提取的關(guān)鍵是理解變量索引x_opt 里前 N 個(gè)是充電功率接著 N 個(gè)是放電功率最后 N1 個(gè)是 SOC。profit 取 -fval因?yàn)?linprog 最小化的目標(biāo)函數(shù)值本身是負(fù)的凈收益。繪圖時(shí)把放電功率畫成負(fù)值這樣一眼能看出充放時(shí)段和電價(jià)的對(duì)應(yīng)關(guān)系。正常的優(yōu)化結(jié)果應(yīng)該是谷段充電功率為正、SOC 抬升峰段放電功率為負(fù)、SOC 回落平段基本不動(dòng)。如果看到谷段不充、峰段不放或者功率在某個(gè)平臺(tái)期反復(fù)跳動(dòng)就是第 5 章要講的典型問(wèn)題。4. 參數(shù)敏感性分析容量配比、老化成本與電價(jià)結(jié)構(gòu)怎么影響收益代碼跑通只是第一步真正決定這套優(yōu)化有沒(méi)有價(jià)值的是參數(shù)。同一個(gè)模型換一組容量配比或電價(jià)數(shù)據(jù)結(jié)論可能完全相反。這一章把三個(gè)最影響結(jié)果的參數(shù)講透。4.1 容量與功率配比C-rate決定策略可行域P_max 與 E_rated 的比值就是 C-rate。C-rate 0.25 意味著滿功率充放電需要 4 小時(shí)0.5C 需要 2 小時(shí)。這個(gè)指標(biāo)直接決定電池能不能在谷段窗口內(nèi)充滿。以本文參數(shù)為例谷段持續(xù) 8 小時(shí)50 kW 對(duì) 200 kWh 是 0.25C理論上 4 小時(shí)就能從空充到滿谷段綽綽有余。但如果谷段只有 2 小時(shí)0.25C 系統(tǒng)只能充 100 kWh浪費(fèi)一半容量這時(shí)把 P_max 提到 100 kW0.5C才能吃滿谷段電量。反過(guò)來(lái)P_max 也不是越大越好——變流器成本隨功率上升而且功率過(guò)大時(shí)充放電過(guò)程更激進(jìn)對(duì)電池壽命的影響在老化成本模型里會(huì)進(jìn)一步放大。C-rate滿充時(shí)間對(duì) 2 小時(shí)谷段的適配0.25C4 小時(shí)只能充一半容量收益受限0.5C2 小時(shí)剛好充滿1C1 小時(shí)能充滿但變流器成本高、壽命壓力大調(diào)參時(shí)建議固定 E_rated掃描 P_max 看凈收益曲線。收益隨 P_max 上升會(huì)先快速增長(zhǎng)、后趨于飽和飽和點(diǎn)就是經(jīng)濟(jì)上最優(yōu)的功率配比。曲線拐點(diǎn)附近再疊加變流器成本才是真正的工程最優(yōu)解。4.2 循環(huán)壽命與老化成本折算目標(biāo)函數(shù)里最容易被忽略的一項(xiàng)老化成本折算的常見(jiàn)做法是線性折舊把電池總購(gòu)置成本平攤到全生命周期放電量上。c_deg 電池單位容量成本 / (循環(huán)壽命 * 放電深度)假設(shè)電池 500 元/kWh循環(huán)壽命 4000 次、按 80% 放電深度計(jì)則 c_deg 500 / (4000 × 0.8) ≈ 0.156 元/kWh也就是每放 1 kWh 電要計(jì)提約 0.16 元損耗。不同循環(huán)壽命下的折算值如下循環(huán)壽命80% DoDc_deg元/kWh 放電2000 次0.314000 次0.166000 次0.10加入 c_deg 之后目標(biāo)函數(shù)里放電收益從 p_t 變成 p_t - c_deg。這個(gè)改動(dòng)會(huì)真實(shí)改變優(yōu)化結(jié)果某些價(jià)差小、原本勉強(qiáng)可套利的時(shí)段扣掉老化成本后變成虧損優(yōu)化器會(huì)主動(dòng)放棄在那段時(shí)間放電。這是模型從「看起來(lái)賺錢」走向「真實(shí)賺錢」的關(guān)鍵一步。注意線性折舊是個(gè)簡(jiǎn)化模型實(shí)際電池老化與放電深度、倍率、溫度呈非線性關(guān)系。做方案對(duì)比用線性折算足夠做精確壽命預(yù)測(cè)則需要換用雨流計(jì)數(shù)或電化學(xué)老化模型那個(gè)就要上 MILP 或啟發(fā)式算法了。4.3 電價(jià)結(jié)構(gòu)與盈虧平衡價(jià)差先算門檻再調(diào)參判斷一個(gè)場(chǎng)景值不值得做儲(chǔ)能套利先算盈虧平衡價(jià)差不用跑代碼p_peak_min p_valley / (eta_ch * eta_dis) c_deg以谷價(jià) 0.3 元/kWh、雙效率 0.95、c_deg 0.156 為例p_peak_min 0.3 / 0.9025 0.156 ≈ 0.49 元/kWh。也就是說(shuō)峰段電價(jià)超過(guò) 0.49 元時(shí)才有正收益超過(guò)越多賺得越多。如果峰價(jià)只有 0.5 元每度電毛利不到一分錢扣掉運(yùn)維成本基本白干如果把效率換成 0.9門檻直接升到 0.53。電價(jià)結(jié)構(gòu)里還有一個(gè)常被忽略的因素是谷段時(shí)長(zhǎng)。谷段只有 2 小時(shí)而 C-rate 是 0.25C 時(shí)電池充不滿實(shí)際套利電量只有理論的一半。所以評(píng)估一個(gè)項(xiàng)目時(shí)不能只看峰谷價(jià)差還要看「價(jià)差 × 可套利電量 × 年循環(huán)次數(shù)」三者乘積。這個(gè)視角能解釋為什么同一套代碼換到某些省份的電價(jià)數(shù)據(jù)后收益大幅縮水——不是代碼錯(cuò)了是電價(jià)結(jié)構(gòu)根本不匹配。5. 儲(chǔ)能運(yùn)行優(yōu)化避坑指南MATLAB 求解中的五個(gè)常見(jiàn)問(wèn)題排查理論和代碼都講完了下面是實(shí)際調(diào)試中踩過(guò)的坑按「現(xiàn)象 → 原因 → 解決」寫每條都是血淚經(jīng)驗(yàn)照著排查能省好幾個(gè)晚上。5.1 linprog 報(bào)錯(cuò)與求解器路徑問(wèn)題現(xiàn)象運(yùn)行到 linprog 那行報(bào)Undefined function or variable linprog或者optimoptions直接報(bào)錯(cuò)。原因MATLAB 沒(méi)裝 Optimization Toolbox。linprog 不是基礎(chǔ)函數(shù)它屬于工具箱optimoptions 是 R2013b 之后引入的寫法老版本也不認(rèn)。解決先在命令行執(zhí)行ver(optim)如果顯示未安裝需要安裝對(duì)應(yīng)工具箱。確認(rèn)裝了但還是報(bào)錯(cuò)檢查是不是用了老版本語(yǔ)法老版本用optimset而不是optimoptions。如果不想裝工具箱就改用 YALMIP 加免費(fèi)求解器 GLPK建模代碼要相應(yīng)改成 YALMIP 語(yǔ)法。5.2 充放電曲線震蕩退化最優(yōu)解的判定與處理現(xiàn)象結(jié)果里相鄰時(shí)段功率劇烈交替比如 19 點(diǎn)放 50 kW、20 點(diǎn)充 50 kW、21 點(diǎn)又放SOC 曲線像鋸齒。目標(biāo)函數(shù)值看著正常但曲線完全不像能實(shí)際執(zhí)行的策略。原因電價(jià)平臺(tái)期存在退化解。LP 保證目標(biāo)函數(shù)最優(yōu)但不保證解唯一當(dāng)多個(gè)方案取得相同收益時(shí)單純形法返回的頂點(diǎn)解可能毫無(wú)規(guī)律。常見(jiàn)于峰段與峰段之間、平段與谷段邊界處。解決最干凈的方案是在目標(biāo)函數(shù)里加切換懲罰項(xiàng)用輔助變量把相鄰時(shí)段功率差值的絕對(duì)值計(jì)入目標(biāo)懲罰系數(shù)取一個(gè)很小的值比如 1e-3 量級(jí)既不影響總收益又能壓住震蕩。工程上更快的做法是如果震蕩幅度不大直接對(duì)結(jié)果做后處理平滑再重新核算目標(biāo)函數(shù)值如果收益損失在 1% 以內(nèi)說(shuō)明震蕩只是退化解平滑后的策略可以放心用。記住先算收益差再?zèng)Q定要不要平滑別憑感覺(jué)改。5.3 終端 SOC 失控優(yōu)化結(jié)果把電池放空的后悔藥現(xiàn)象結(jié)果里最后一個(gè)峰段電池被放到 SOC 下限曲線看起來(lái)很「完美」地把每一度電都賺干凈了但第二天這個(gè)策略根本沒(méi)法繼續(xù)執(zhí)行——電池沒(méi)電了。原因沒(méi)加終端 SOC 約束。初始 SOC 是給定的、免費(fèi)的優(yōu)化器會(huì)把它當(dāng)成可消耗資源在最后一個(gè)高價(jià)時(shí)段全部變現(xiàn)。解決在等式約束里加上 SOC(N1) SOC_init。如果做的是滾動(dòng)優(yōu)化每一步的初始 SOC 來(lái)自實(shí)測(cè)不需要強(qiáng)綁終端值但至少加一個(gè) SOC(N1) ≥ SOC_min 的下限約束兜底。這個(gè)坑極其隱蔽因?yàn)閱稳战Y(jié)果看起來(lái)完全正常要連續(xù)跑多日才會(huì)暴露。5.4 最優(yōu)收益為負(fù)先查數(shù)據(jù)再查模型現(xiàn)象linprog 正常返回exitflag 是 1但凈收益算出來(lái)是負(fù)的。原因要么電價(jià)數(shù)據(jù)有問(wèn)題要么這個(gè)場(chǎng)景本身就不經(jīng)濟(jì)。常見(jiàn)數(shù)據(jù)錯(cuò)誤包括峰谷價(jià)差寫反、電價(jià)單位不一致比如誤用元/MWh、效率參數(shù)填成 0.95 但實(shí)際想表達(dá) 95% 卻寫成了 95。解決先用第 4.3 節(jié)的盈虧平衡公式手算一遍門檻價(jià)看峰段電價(jià)是否高于門檻。如果高于門檻但結(jié)果還是負(fù)逐個(gè)檢查效率、電價(jià)單位、老化成本取值。如果確實(shí)低于門檻那負(fù)收益就是正確結(jié)論——說(shuō)明這個(gè)場(chǎng)景不該上儲(chǔ)能套利模型給你省了一筆投資這本身就是價(jià)值。5.5 單位不一致導(dǎo)致的數(shù)值翻車現(xiàn)象收益大得離譜比如 200 kWh 的電池一天賺了幾萬(wàn)塊或者 SOC 曲線亂跳、約束怎么調(diào)都矛盾。原因單位混用了——容量用 MWh、功率用 kW、電價(jià)用元/MWh導(dǎo)致約束矩陣?yán)锔黜?xiàng)數(shù)量級(jí)差好幾個(gè)量級(jí)數(shù)值病態(tài)。解決統(tǒng)一到 kW / kWh / 元/kWh 三個(gè)基礎(chǔ)單位。改完數(shù)據(jù)后在模型構(gòu)建前加一句打印核對(duì)各數(shù)組的最大最小值尤其是 E_rated 和 P_max 的比例關(guān)系。單位錯(cuò)誤這種事排查半天往往最后發(fā)現(xiàn)就是一處分母寫錯(cuò)屬于最冤的翻車方式。6. 進(jìn)階從單日優(yōu)化到滾動(dòng)優(yōu)化與規(guī)則策略對(duì)比驗(yàn)證6.1 滾動(dòng)優(yōu)化只執(zhí)行第一步每小時(shí)重算單日優(yōu)化適合算理論收益上限但真拿去運(yùn)行會(huì)有兩個(gè)問(wèn)題一是電價(jià)預(yù)測(cè)不可能完全準(zhǔn)二是實(shí)際執(zhí)行的充放電量和計(jì)劃總會(huì)有偏差SOC 越偏越遠(yuǎn)。常見(jiàn)做法是把單日求解封裝成函數(shù)然后用滾動(dòng)優(yōu)化的框架循環(huán)調(diào)用每個(gè)小時(shí)基于當(dāng)前實(shí)測(cè) SOC 和未來(lái)一段時(shí)間的電價(jià)預(yù)測(cè)重新求解一次但只執(zhí)行第一步。%% 滾動(dòng)優(yōu)化主循環(huán) T_total 24 * 7; % 仿真一周 H 24; % 優(yōu)化窗口長(zhǎng)度 soc_now SOC_init; result_q_ch zeros(T_total, 1); result_q_dis zeros(T_total, 1); for k 1:T_total % 取未來(lái) H 小時(shí)電價(jià)此處用真實(shí)電價(jià)代替預(yù)測(cè) idx mod((k-1:kH-2), 24) 1; price_fc price(idx); % 調(diào)用封裝好的單次優(yōu)化函數(shù) [q_ch, q_dis, ~] solve_storage_one_day(price_fc, E_rated, P_max, ... soc_now, SOC_min, SOC_max, eta_ch, eta_dis, dt); % 只執(zhí)行第一步 result_q_ch(k) q_ch(1); result_q_dis(k) q_dis(1); % 按實(shí)際充放電量更新 SOC soc_now soc_now (eta_ch * q_ch(1) - q_dis(1) / eta_dis) * dt / E_rated; end這一步的意義在于讓優(yōu)化永遠(yuǎn)從真實(shí)狀態(tài)出發(fā)。單日優(yōu)化里終端 SOC 約束導(dǎo)致的可調(diào)度容量損失在滾動(dòng)框架下自然消失——每天每步都在重新初始化。窗口長(zhǎng)度 H 的選取有個(gè)折中太短看不到完整谷峰周期太長(zhǎng)對(duì)預(yù)測(cè)精度要求高且求解規(guī)模變大我一般取 24 小時(shí)。6.2 與規(guī)則策略對(duì)比用基準(zhǔn)線驗(yàn)證優(yōu)化增益滾動(dòng)優(yōu)化跑出來(lái)的結(jié)果好不好看需要跟一個(gè)規(guī)則基準(zhǔn)對(duì)比。最簡(jiǎn)單的規(guī)則策略是「谷段充滿、峰段放完」谷段按最大功率充電峰段按最大功率放電其他時(shí)段不動(dòng)。把同一組電價(jià)和參數(shù)喂給規(guī)則策略算出收益再和優(yōu)化結(jié)果對(duì)比。兩者的差值就是優(yōu)化帶來(lái)的增益。如果增益很小說(shuō)明電價(jià)結(jié)構(gòu)太簡(jiǎn)單、規(guī)則策略已經(jīng)接近最優(yōu)如果增益明顯說(shuō)明復(fù)雜電價(jià)結(jié)構(gòu)下優(yōu)化確實(shí)值得。這個(gè)對(duì)比在項(xiàng)目匯報(bào)里特別有用——評(píng)審不會(huì)只看絕對(duì)收益更關(guān)心你比拍腦袋策略多賺了多少。我最早跑這個(gè)題目時(shí)沒(méi)加終端 SOC 約束結(jié)果每日收益虛高跟規(guī)則策略一對(duì)比才發(fā)現(xiàn)自己把初始存電免費(fèi)放掉了。后來(lái)把約束補(bǔ)齊、單位統(tǒng)一收益數(shù)字才真正經(jīng)得起復(fù)核。這類代碼模型寫簡(jiǎn)單一點(diǎn)沒(méi)關(guān)系邊界條件和單位才是真正翻車的地方。希望幫到你。本文還有配套的精品資源點(diǎn)擊獲取