光儲(chǔ)互補(bǔ)微電網(wǎng)Simulink建模與仿真實(shí)戰(zhàn)詳解)
1. 風(fēng)光儲(chǔ)互補(bǔ)微電網(wǎng)建模的全局思路1.1 為什么需要把光伏、風(fēng)電、儲(chǔ)能放在同一個(gè)系統(tǒng)里建模仿真先說一個(gè)我自己的判斷很多人一上來(lái)就奔著“怎么用Simulink搭模塊”去但真正把風(fēng)光儲(chǔ)互補(bǔ)微電網(wǎng)跑起來(lái)、跑穩(wěn)、跑出能寫進(jìn)論文里的曲線容易卡住的地方恰恰不是某個(gè)模塊怎么點(diǎn)而是對(duì)整個(gè)系統(tǒng)的理解。這個(gè)項(xiàng)目標(biāo)題里的關(guān)鍵詞是“風(fēng)光儲(chǔ)互補(bǔ)微電網(wǎng)”拆開看就是三個(gè)電源和一個(gè)電網(wǎng)架構(gòu)的問題光伏發(fā)電、風(fēng)力發(fā)電、儲(chǔ)能單元再加上把它們連起來(lái)的母線、負(fù)荷和控制系統(tǒng)。所謂“互補(bǔ)”本質(zhì)上是利用光伏和風(fēng)電在時(shí)間尺度上的天然錯(cuò)峰特性配合儲(chǔ)能來(lái)平抑功率波動(dòng)讓微電網(wǎng)在并網(wǎng)和孤島兩種模式下都能穩(wěn)定運(yùn)行。那為什么要建模說白了你不能為了驗(yàn)證一個(gè)控制策略真的去搭一套光伏板、一架風(fēng)機(jī)和一個(gè)儲(chǔ)能電站——成本太高而且很多邊界條件和故障工況在實(shí)物上是很難安全再現(xiàn)的。Simulink建模的價(jià)值就在這里用數(shù)學(xué)模型去等效物理過程把輻照度變化、風(fēng)速波動(dòng)、負(fù)荷投切這些“外部擾動(dòng)”注入系統(tǒng)看控制策略能不能扛得住扛到什么程度。我在做這類項(xiàng)目時(shí)習(xí)慣先畫一張系統(tǒng)級(jí)的功率流程圖。光伏和風(fēng)電是源儲(chǔ)能是“緩沖池”負(fù)荷是消費(fèi)端而控制策略負(fù)責(zé)告訴整個(gè)系統(tǒng)“什么時(shí)候該充電、什么時(shí)候該放電、什么時(shí)候切并網(wǎng)模式、什么時(shí)候轉(zhuǎn)孤島”。這張圖畫清楚之后Simulink里怎么搭模塊就只是時(shí)間問題了。1.2 系統(tǒng)架構(gòu)和母線電壓等級(jí)怎么選接下來(lái)一個(gè)重要決策是母線電壓等級(jí)。風(fēng)光儲(chǔ)微電網(wǎng)里常見的母線電壓有三種48V低壓直流、400V直流、或者690V交流。我實(shí)測(cè)的體驗(yàn)是如果是課程設(shè)計(jì)或畢業(yè)設(shè)計(jì)性質(zhì)的建模建議用400V直流母線原因有兩個(gè)。第一400V直流母線可以直接對(duì)應(yīng)市面上成熟的光伏逆變器、儲(chǔ)能變流器PCS的直流側(cè)電壓范圍后續(xù)如果做硬件在環(huán)測(cè)試或者真機(jī)實(shí)驗(yàn)參數(shù)可以平順平移過去。第二Boost變換器在400V母線下的占空比不至于太極端調(diào)制起來(lái)比較平滑。48V母線雖然更安全、更適合小功率實(shí)驗(yàn)室但Boost占空比太高甚至接近極限很多控制規(guī)律的效果會(huì)被壓縮。交流母線的情況要復(fù)雜一些。如果微電網(wǎng)里有交流負(fù)載或者需要并網(wǎng)你就必須引入三相逆變器或者三相整流器然后面對(duì)交流側(cè)電壓幅值、頻率、相位的同步問題。Simulink里這個(gè)同步是靠PLL鎖相環(huán)實(shí)現(xiàn)的初學(xué)者在這里栽跟頭的概率非常高——我后面會(huì)在排查章節(jié)里專門展開。1.3 建模仿真的工具鏈選型工具鏈上我用的是Matlab/Simulink本體加Simscape Electrical Toolbox以前叫SimPowerSystems這是目前2025年前后做微電網(wǎng)仿真比較主流的選擇。Simulink負(fù)責(zé)控制策略的搭建Simscape Electrical負(fù)責(zé)主電路光伏陣列、電池、變換器、逆變器、電網(wǎng)的物理建模。有兩條路線可以選擇一是全部用Simscape自帶的庫(kù)——比如光伏用PV Array模塊、電池用Battery模塊這種方案建模速度快模塊內(nèi)部已經(jīng)包含了工程上常見的等效電路但缺點(diǎn)是有些內(nèi)部參數(shù)被封裝起來(lái)想精確控制某個(gè)細(xì)節(jié)時(shí)不太順手。二是用基本的電阻、電感、電容、開關(guān)管自己搭主電路再用方程寫控制邏輯這種方案自定義空間大、模型透明但工程量大且容易出數(shù)值問題。我個(gè)人的建議是如果目標(biāo)是快速出一個(gè)“風(fēng)光儲(chǔ)互補(bǔ)微電網(wǎng)”的運(yùn)行仿真并重點(diǎn)展現(xiàn)控制策略和模式切換邏輯就選第一條路線如果目標(biāo)是做學(xué)術(shù)論文、想深度分析某個(gè)特定的電力電子拓?fù)浠騍VPWM調(diào)制細(xì)節(jié)就選第二條路線。這篇文章后面以第一條路線為主線展開。2. 三大核心電源模塊的Simulink實(shí)現(xiàn)細(xì)節(jié)2.1 光伏組件建模從等效電路到PV Array參數(shù)換算光伏建模這塊很多人直接默認(rèn)用Simscape里的PV Array模塊這沒問題關(guān)鍵是模塊參數(shù)怎么填得合理。PV Array模塊內(nèi)部實(shí)際上是一個(gè)基于單二極管五參數(shù)模型的等效電路——光生電流源、一個(gè)二極管、一個(gè)串聯(lián)電阻和一個(gè)并聯(lián)電阻。其I-V特性表達(dá)式是I Iph - I0 × [exp((V I × Rs) / (a × Vt)) - 1] - (V I × Rs) / Rsh其中Iph是光生電流I0是飽和電流Rs和Rsh是串聯(lián)和并聯(lián)等效電阻a是二極管理想因子Vt是熱電壓。這個(gè)方程的五個(gè)參數(shù)在模塊里通常不會(huì)直接暴露而是讓你填“廠家數(shù)據(jù)表里的標(biāo)稱值”比如開路電壓Voc、短路電流Isc、最大功率點(diǎn)電壓Vmp、最大功率點(diǎn)電流Imp。填參數(shù)有一個(gè)最重要的換算點(diǎn)數(shù)據(jù)表給的標(biāo)稱值是在標(biāo)準(zhǔn)測(cè)試條件STC輻照度1000W/m2溫度25°C下的而你的場(chǎng)景可能是400W/m2或者35°CSimulink靠?jī)?nèi)部的溫度和輻照度修正公式自動(dòng)算出來(lái)。不過你需要在PNel最大輸出功率這一欄確認(rèn)一下填的是峰值功率還是模塊串并聯(lián)后的陣列總功率——我見過不少人在這一步算錯(cuò)導(dǎo)致后續(xù)MPPT曲線完全不對(duì)。我一般先設(shè)置單個(gè)組件Module的參數(shù)再通過串聯(lián)數(shù)Ns和并聯(lián)數(shù)Np組成陣列。比如選一塊峰值功率約280W的組件Voc38V、Isc9.2A、Vmp31V、Imp8.8A然后用4塊串聯(lián)組建一串48V左右電壓、8串并聯(lián)組建一個(gè)陣列總功率約9kW。這個(gè)9kW的光伏陣列配一個(gè)10kW級(jí)的Boost變換器基本符合“互補(bǔ)微電網(wǎng)里光伏作為主力電源之一”的設(shè)定。MPPT最大功率點(diǎn)跟蹤算法我優(yōu)先推薦擾動(dòng)觀察法Simulink里實(shí)現(xiàn)起來(lái)就是“采樣光伏輸出電壓和電流 → 算當(dāng)前功率 → 和上一拍功率比 → 決定下一拍的參考電壓是加還是減”。我給一個(gè)在實(shí)際項(xiàng)目里驗(yàn)證過的擾動(dòng)觀察法參考代碼風(fēng)格以函數(shù)方式封裝用持久變量或Unit Delay實(shí)現(xiàn)上一拍記憶function Vref mppt_pvo(P_sim, P_prev, V_sim, V_prev, dV) % P_sim 當(dāng)前功率P_prev 上一拍功率 % V_sim 當(dāng)前電壓V_prev 上一拍電壓 if P_sim P_prev if V_sim V_prev Vref V_sim dV; % 功率上升且電壓上升 → 繼續(xù)加電壓 else Vref V_sim - dV; % 功率上升但電壓下降 → 繼續(xù)減電壓 end else if V_sim V_prev Vref V_sim - dV; % 功率下降且電壓上升 → 反向減電壓 else Vref V_sim dV; % 功率下降且電壓下降 → 反向加電壓 end end end步長(zhǎng)dV的選擇直接影響MPPT精度和動(dòng)態(tài)速度。我實(shí)測(cè)下來(lái)對(duì)于400V母線、10kW級(jí)光伏陣列dV取0.5V左右比較合適太小了在輻照突變時(shí)半天追不上最大功率點(diǎn)太大了又在穩(wěn)態(tài)時(shí)來(lái)回震蕩造成輸出功率紋波。2.2 風(fēng)電模型的簡(jiǎn)化與等效風(fēng)速輸入風(fēng)電建模比光伏麻煩因?yàn)楣夥恍枰拜椪斩葴囟取眱蓚€(gè)外部變量而風(fēng)機(jī)涉及空氣動(dòng)力學(xué)、傳動(dòng)鏈、發(fā)電機(jī)、變流器好幾個(gè)環(huán)節(jié)。但在Simulink里做微電網(wǎng)級(jí)仿真時(shí)我們關(guān)心的不是風(fēng)機(jī)內(nèi)部某一片槳葉的受力而是它反映到電網(wǎng)側(cè)的功率輸出特性和波動(dòng)特性。所以完全可以做合理簡(jiǎn)化。我推薦用“風(fēng)速-功率曲線”模型。核心公式是空氣動(dòng)力學(xué)捕獲的風(fēng)功率P 0.5 × ρ × π × R2 × v3 × Cp(λ, β)ρ是空氣密度R是風(fēng)輪半徑v是風(fēng)速Cp是風(fēng)能利用系數(shù)——它是葉尖速比λ和槳距角β的函數(shù)。在Simulink里實(shí)現(xiàn)時(shí)Cp可以用查表或解析公式來(lái)算如果不做槳距角控制β固定為0°Cp函數(shù)可以直接用一個(gè)多項(xiàng)式或者三角函數(shù)擬合。一個(gè)常用近似公式是Cp(λ, β) 0.22 × (116/λi - 0.4β - 5) × exp(-12.5/λi)其中λi 1 / (1/(λ 0.08β) - 0.035/(β3 1))。這個(gè)公式是學(xué)術(shù)界流行的一個(gè)近似模型實(shí)測(cè)下來(lái)在額定風(fēng)速附近的擬合效果還不錯(cuò)。然后加上傳動(dòng)鏈和發(fā)電機(jī)模型。如果你不想把模型搞得太重傳動(dòng)鏈可以用單質(zhì)量塊模型代替兩質(zhì)量塊模型——兩質(zhì)量塊的扭振模態(tài)在微網(wǎng)運(yùn)行分析里影響不大但會(huì)顯著拉低仿真速度。發(fā)電機(jī)部分我習(xí)慣直接用一個(gè)“三相電壓源功率控制”的組合代替完整的異步發(fā)電機(jī)/永磁同步發(fā)電機(jī)FOC矢量控制這樣系統(tǒng)主電路不會(huì)過于復(fù)雜風(fēng)電側(cè)就以一個(gè)“可控電流源/可控功率源”的形式體現(xiàn)。風(fēng)速輸入也別圖省事用恒定值——那樣子根本看不到“風(fēng)光互補(bǔ)”的優(yōu)勢(shì)。你需要一個(gè)隨時(shí)間變化的風(fēng)速序列。可以是實(shí)測(cè)數(shù)據(jù)導(dǎo)入From Workspace或From File也可以用幾個(gè)正弦波疊加自噪聲模塊來(lái)合成模擬風(fēng)速。我常用的是組合信號(hào)一個(gè)平均風(fēng)速8m/s、幅值2m/s、周期約50s的慢變動(dòng)分量再疊加一個(gè)幅值較小的隨機(jī)噪聲分量。這樣風(fēng)機(jī)的輸出功率會(huì)呈現(xiàn)明顯的波動(dòng)特性儲(chǔ)能的平抑功能就“有事可做”了。2.3 儲(chǔ)能蓄電池建模與SOC估算儲(chǔ)能是微電網(wǎng)里真正意義上的“壓艙石”。在孤島模式下儲(chǔ)能要維持母線電壓和頻率相當(dāng)于一個(gè)可充電的電壓源在并網(wǎng)模式下儲(chǔ)能要做削峰填谷和功率平抑。Simulink里的Battery模塊已經(jīng)封裝了多種電池類型的等效電路我常用的設(shè)置是“Lithium-Ion Battery”額定電壓按母線要求設(shè)計(jì)。電池模塊內(nèi)部本質(zhì)上是“受控電壓源內(nèi)阻”結(jié)構(gòu)開路電壓是SOC的函數(shù)而內(nèi)阻在充放電和老化過程中會(huì)有變化。參數(shù)設(shè)置上重點(diǎn)關(guān)注額定容量Ah、初始SOC百分比、響應(yīng)時(shí)間常數(shù)。響應(yīng)時(shí)間常數(shù)這個(gè)參數(shù)挺關(guān)鍵的——它決定了電池在一個(gè)功率階躍之后多久能達(dá)到新的穩(wěn)態(tài)電壓如果設(shè)得太短比如0.001s模型會(huì)表現(xiàn)得像一個(gè)剛性電壓源在動(dòng)態(tài)仿真里容易引起電流突變?nèi)绻O(shè)得太長(zhǎng)又會(huì)顯得儲(chǔ)能反應(yīng)遲鈍母線電壓在負(fù)載突變時(shí)可能壓不住。我建議在1s到10s之間選取具體按容量和場(chǎng)景微調(diào)。SOC荷電狀態(tài)的估算和控制邏輯深度綁定。Simulink的Battery模塊里會(huì)實(shí)時(shí)給出SOC信號(hào)這個(gè)信號(hào)是基于安時(shí)積分法算出來(lái)的SOC(t) SOC? - (1/Q)∫i(t)dtQ是電池總?cè)萘?。如果你需要看SOC的變化曲線直接把這個(gè)引出即可如果你要自己設(shè)計(jì)SOC的滯回區(qū)間比如低于20%強(qiáng)制切負(fù)載、高于95%限制充電功率需要寫一點(diǎn)邏輯判斷這個(gè)時(shí)候可以用Relay模塊或者Stateflow Chart來(lái)做。這里有一個(gè)我在實(shí)驗(yàn)里踩過的坑Battery模塊在SOC接近100%或接近0%時(shí)等效電路里的開路電壓和非線性內(nèi)阻會(huì)發(fā)生劇烈變化很容易導(dǎo)致仿真步長(zhǎng)突然縮得非常小整段仿真像蝸牛一樣慢。我在遇到這種情況時(shí)會(huì)把SOC的工作區(qū)間約束在20%~95%同時(shí)設(shè)定仿真終止條件“SOC超出該區(qū)間即觸發(fā)控制策略切換”既保護(hù)模型數(shù)值穩(wěn)定性又貼近工程實(shí)際。3. 控制策略設(shè)計(jì)與系統(tǒng)聯(lián)調(diào)3.1 主從控制與下垂控制微電網(wǎng)的核心控制框架在風(fēng)光儲(chǔ)互補(bǔ)微電網(wǎng)中控制策略通常不是單一模塊搞定而是分層分布式的。我實(shí)際項(xiàng)目中用得最順的是主從控制Master-Slave架構(gòu)儲(chǔ)能作為主電源Master負(fù)責(zé)維持母線電壓和頻率的穩(wěn)定光伏和風(fēng)電作為從電源Slave只負(fù)責(zé)按最大功率或限功率輸出不參與頻率調(diào)節(jié)。孤島模式下儲(chǔ)能系統(tǒng)采用V/f控制。所謂V/f控制說白了就是“我不管功率來(lái)自哪里都是通過調(diào)節(jié)變換器的開關(guān)來(lái)讓輸出電壓幅值和頻率保持恒定”。具體的Simulink實(shí)現(xiàn)方式是把母線電壓的d軸參考值設(shè)為母線電壓峰值的√2倍或按你的坐標(biāo)定義q軸參考值設(shè)為0頻率設(shè)為50Hz然后用PI控制器分別調(diào)節(jié)d軸和q軸的電流再經(jīng)PWM/PFM生成開關(guān)信號(hào)。PI參數(shù)我一般先按系統(tǒng)時(shí)間常數(shù)估算再通過試湊法微調(diào)——比例項(xiàng)先把穩(wěn)態(tài)偏差砍掉積分項(xiàng)再慢慢消掉殘差。并網(wǎng)模式下儲(chǔ)能就不再是嚴(yán)格的V/f主源了而是轉(zhuǎn)為PQ控制恒定有功/無(wú)功控制跟電網(wǎng)電壓走。這個(gè)過程里的“模式切換”是整個(gè)建模里最容易仿真出錯(cuò)的地方因?yàn)閂/f控制切換到PQ控制時(shí)控制器的初始狀態(tài)和積分器“記憶”是連續(xù)性的問題——如果切換瞬間控制器的積分初始值沒有對(duì)齊就會(huì)出現(xiàn)電流突變甚至振蕩發(fā)散。解決辦法是在切換邏輯里加一個(gè)狀態(tài)初始化比如把PQ控制器的積分項(xiàng)初值設(shè)為切換前一刻的V/f控制輸出或者干脆用平滑切換的滯回/斜坡過渡。下垂控制Droop Control是另一種我經(jīng)常用的方案特別適合多個(gè)儲(chǔ)能單元并聯(lián)或者微電網(wǎng)之間有組網(wǎng)需求的場(chǎng)景相當(dāng)于從邏輯上“模擬同步發(fā)電機(jī)的自調(diào)節(jié)特性”。頻率下垂是f f? - m×P電壓下垂是V V? - n×Q其中m是頻率有功下垂系數(shù)n是電壓無(wú)功下垂系數(shù)。下垂系數(shù)的設(shè)計(jì)原則是輸出功率達(dá)到滿額時(shí)頻率偏移不超過±0.5Hz比如額定功率100kW、額定頻率50Hz那么m (50-49.5)/100k 5×10?? Hz/W。這類參數(shù)在仿真里很好驗(yàn)證突加負(fù)載后頻率降多少理論計(jì)算和仿真曲線一對(duì)比就知道垂感系數(shù)設(shè)合理沒有。3.2 Boost變換器與PWM驅(qū)動(dòng)參數(shù)的計(jì)算過程儲(chǔ)能變流器或者光伏Boost變換器的電感電容參數(shù)直接決定電流紋波和電壓紋波是否合格。這里我以光伏Boost為例給出一個(gè)完整的計(jì)算流程。假設(shè)光伏陣列MPPT電壓近似為Boost輸入電壓Vin為120V母線電壓Vout是400V開關(guān)頻率f_sw是10kHz額定輸出功率P是10kW。Boost變換器的占空比D 1 - Vin/Vout 1 - 120/400 0.7。電感量按下式計(jì)算L Vin×D / (2×f_sw×ΔIL)。ΔIL是允許的電流紋波通常取輸入電流的20%~30%。輸入電流平均值I P/Vin 10000/120 ≈ 83.3A取紋波系數(shù)25%那么ΔIL 0.25×83.3 ≈ 20.8A。代入公式L 120×0.7 / (2×10000×20.8) ≈ 202μH。實(shí)際我會(huì)稍微取大一點(diǎn)到220μH留一點(diǎn)余量。輸出電容按紋波電壓要求算ΔVout取母線電壓的1%也就是4V那么C P×D / (f_sw×ΔVout×Vout) 10000×0.7 / (10000×4×400) ≈ 437μF工程上取470μF即可。這兩個(gè)參數(shù)在Simulink里直接輸入到Buck/Boost變換器模塊或者你自建的開關(guān)管模型里。PWM生成用Simulink自帶的PWM Generator——注意載波頻率要跟你算的f_sw一致仿真步長(zhǎng)至少設(shè)成載波周期的1/20以上否則PWM波形會(huì)很粗糙。具體來(lái)說10kHz的PWM意味著周期是100μs仿真步長(zhǎng)最好在5μs以下最好用變步長(zhǎng)求解器把最大步長(zhǎng)限制在10μs以內(nèi)不然“開關(guān)事件”在步長(zhǎng)太大時(shí)會(huì)被跳過波形失真。3.3 從空白模型到完整系統(tǒng)的搭建流程我平時(shí)從零開始搭這樣一個(gè)系統(tǒng)大致分四步走。第一步是搭主電路骨架放直流母線、光伏Boost、風(fēng)機(jī)等效功率源、儲(chǔ)能電池Buck/Boost雙向變換器、兩個(gè)可控負(fù)載一個(gè)固定阻性負(fù)載、一個(gè)動(dòng)態(tài)負(fù)載用于測(cè)試擾動(dòng)。這一步只看波形是否正常先不接任何閉環(huán)控制光伏給固定占空比儲(chǔ)能直接并聯(lián)在母線上做“硬充電”目的是確認(rèn)電氣連接沒有錯(cuò)誤母線電壓能撐起來(lái)。第二步是逐個(gè)模塊接入控制先給光伏接MPPT觀察母線電壓和光伏輸出曲線是否隨輻照度變化合理追蹤再給儲(chǔ)能接V/f控制或下垂控制看看負(fù)載突變時(shí)電壓能不能穩(wěn)住最后給風(fēng)機(jī)接功率控制同時(shí)把風(fēng)速擾動(dòng)加進(jìn)去觀察系統(tǒng)功率波動(dòng)情況。第三步是聯(lián)調(diào)互補(bǔ)邏輯把風(fēng)光儲(chǔ)三個(gè)源都接上固定輻照度、固定風(fēng)速、固定負(fù)載先跑一段穩(wěn)態(tài)確認(rèn)“三源一荷”平衡。然后突然把輻照度從1000W/m2降到300W/m2看看儲(chǔ)能是補(bǔ)功率缺口還是母線電壓瞬間跌落觀察動(dòng)態(tài)響應(yīng)。這是整個(gè)項(xiàng)目最核心的一個(gè)驗(yàn)證場(chǎng)景。第四步是模式切換測(cè)試設(shè)置一個(gè)并網(wǎng)變流器令微電網(wǎng)在并網(wǎng)和孤島兩種模式之間切換。模式切換的動(dòng)作可以用手動(dòng)斷路器Breaker模塊模擬也可以設(shè)定自動(dòng)判據(jù)比如檢測(cè)到大電網(wǎng)異常失壓時(shí)自動(dòng)切換。我建議先在穩(wěn)態(tài)下手動(dòng)切換出現(xiàn)哨叫再逐步加入自動(dòng)判據(jù)這樣排查問題會(huì)輕松很多。4. 仿真調(diào)試中的常見問題與實(shí)戰(zhàn)排查4.1 仿真速度慢從求解器到離散化的一整套優(yōu)化這是所有人在做微電網(wǎng)仿真時(shí)都會(huì)碰到的問題。我剛開始做10kW風(fēng)光儲(chǔ)微電網(wǎng)時(shí)用一個(gè)普通ODE45變步長(zhǎng)求解器設(shè)置仿真時(shí)長(zhǎng)60秒結(jié)果跑了十幾分鐘還沒跑完。問題出在好幾個(gè)點(diǎn)上。第一PWM開關(guān)頻率附近的詳細(xì)電力電子模型會(huì)產(chǎn)生非常多的事件點(diǎn)每次開關(guān)動(dòng)作必須把步長(zhǎng)縮小到能精確捕捉邊沿的量級(jí)結(jié)果就是大量微小步長(zhǎng)累積。第二Simscape Electrical的物理模型求解是隱式迭代每一次都涉及矩陣求逆模型節(jié)點(diǎn)一多更慢。第三MPPT和V/f控制里如果用了過快的采樣時(shí)間比如1μs但意義不大純屬浪費(fèi)計(jì)算。我的優(yōu)化方案是三管齊下。首先是換求解器把ODED45換成ode23tb或ode15s。這類剛性求解器對(duì)電路模型更友好。然后是最大步長(zhǎng)控制把Max Step Size設(shè)為50μs或100μs保證足夠捕捉10kHz~20kHz的開關(guān)頻率又不至于沿事件無(wú)限細(xì)分。最后是離散化改造如果控制目標(biāo)是驗(yàn)證系統(tǒng)級(jí)功率平衡可以把Boost、雙向DC/DC這些平均化——用受控電壓源/電流源加傳遞函數(shù)來(lái)等效拓?fù)鋭?dòng)態(tài)保留電壓電流包絡(luò)而不是鋸齒波開關(guān)紋波。平均化以后仿真速度能提升1~2個(gè)數(shù)量級(jí)代價(jià)是需要你自己寫一點(diǎn)平均模型等效公式。4.2 代數(shù)環(huán)導(dǎo)致的仿真發(fā)散與處理代數(shù)環(huán)Algebraic Loop是Simulink仿真里最折磨人的問題之一在微電網(wǎng)模型里尤其常見。它本質(zhì)上是一個(gè)信號(hào)回路在沒有狀態(tài)量緩沖的情況下直接形成了瞬時(shí)依賴閉環(huán)比如“電壓 → 控制器 → 占空比 → 平均模型 → 電壓”這個(gè)回路里如果控制器只是一個(gè)比例增益而沒有狀態(tài)沒有積分或一階慣性Simulink就必須在每一拍迭代求解這個(gè)代數(shù)環(huán)一旦迭代不收斂就直接報(bào)錯(cuò)或者發(fā)散。處理代數(shù)環(huán)的思路很直接在環(huán)上插入一個(gè)離散的狀態(tài)保持環(huán)節(jié)比如Unit Delay單位延遲、Memory模塊或一個(gè)時(shí)間常數(shù)極小的低通濾波器。我實(shí)踐里最常用的是在“母線電壓采樣 → 控制器”之間串一個(gè)Transfer Function模塊形式是1/(Ts1)T取1e-4s左右。這相當(dāng)于給反饋信號(hào)加了非常輕微的低通濾波物理上不影響控制效果但能有效打破代數(shù)環(huán)讓Simulink順利推進(jìn)。需要留一張“排查問題速查表”可以大幅壓縮調(diào)參時(shí)間。把報(bào)錯(cuò)類型、可能原因、定位方法、解決方案四欄列清楚每次遇到問題直接查表比翻文檔高效太多。4.3 各階段典型異?,F(xiàn)象速查表異常現(xiàn)象可能原因定位思路解決方案仿真報(bào)錯(cuò)Derivative of state not finite電感/電容初始值不合理或控制輸出出現(xiàn)NaN斷開控制環(huán)先跑開環(huán)看主電路波形修正初始條件給控制器輸出加限幅Saturation母線電壓起不來(lái)幾乎為0Boost占空比給定不正確或光伏組件參數(shù)矛盾檢查PV模塊的Vmp是否大于Boost輸出電壓的1/1-D校正占空比范圍檢查Voc與Vmp一致性MPPT后功率一直不增加擾動(dòng)步長(zhǎng)過大導(dǎo)致來(lái)回震蕩噪聲觀察P-V曲線是否反復(fù)跳變減小dV或改用增量電導(dǎo)法孤島切換瞬間電壓暴跌或飛升模式切換時(shí)控制器積分初值未對(duì)齊加日志記錄切換前后參考值差異切換時(shí)給PI積分器初始化設(shè)置SOC先不變后突然跳變安時(shí)積分初始化不對(duì)或容量單位填錯(cuò)檢查電池容量單位Ah是否匹配電流單位A重設(shè)容量參數(shù)SOC范圍約束20%~95%風(fēng)速不變但風(fēng)機(jī)輸出功率振蕩Cp函數(shù)查表沒有覆蓋當(dāng)前葉尖速比范圍畫Cp-λ曲線看工作點(diǎn)是否落在有效區(qū)擴(kuò)展查表范圍或改用分段函數(shù)插值仿真中途變慢到走不動(dòng)開關(guān)事件過密或進(jìn)入高剛性數(shù)值區(qū)查看步長(zhǎng)監(jiān)視器是否頻繁觸發(fā)小步長(zhǎng)換ode23tb采用平均化模型限Max Step Size4.4 幾個(gè)來(lái)自實(shí)測(cè)現(xiàn)場(chǎng)的調(diào)參經(jīng)驗(yàn)最后分享幾個(gè)只靠文檔絕對(duì)學(xué)不到的細(xì)節(jié)。第一風(fēng)機(jī)模型的風(fēng)速擾動(dòng)信號(hào)不要直接和直流母線功率模型相乘——因?yàn)轱L(fēng)速含有高頻成分功率目標(biāo)是三次方關(guān)系直接乘會(huì)讓功率信號(hào)劇烈跳變。我會(huì)讓風(fēng)速信號(hào)先經(jīng)過一個(gè)一階低通濾波時(shí)間常數(shù)2s到5s模擬風(fēng)輪和傳動(dòng)系統(tǒng)的慣性輸出功率曲線會(huì)平滑很多也更真實(shí)。第二光伏MPPT在輻照度突變時(shí)有一個(gè)經(jīng)典問題PO算法在輻照度快速變化時(shí)可能“誤判方向”因?yàn)镻O假設(shè)功率變化主要由電壓調(diào)整引起但輻照突變本身也會(huì)改變功率算法可能把兩個(gè)原因混在一起。如果輻照突變很頻繁比如你要模擬多云天建議改用增量電導(dǎo)法它對(duì)輻照變化更魯棒。Simulink里增量電導(dǎo)的實(shí)現(xiàn)其實(shí)就是算dP/dV并判斷電流變化方向代碼量比PO多幾行但值得的。第三電池模型的“Intial SOC”別隨便設(shè)成100%。很多人在建模時(shí)圖省事直接默認(rèn)滿電結(jié)果一跑仿真就進(jìn)入電池過壓保護(hù)控制表現(xiàn)成“充不進(jìn)電、母線卻缺電”的怪現(xiàn)象。我給儲(chǔ)能一個(gè)初始SOC60%留足上下調(diào)節(jié)空間微電網(wǎng)的動(dòng)態(tài)過程才能真正體現(xiàn)“互補(bǔ)”兩個(gè)字的價(jià)值。結(jié)尾我在實(shí)際做這個(gè)風(fēng)光儲(chǔ)互補(bǔ)微電網(wǎng)Simulink仿真的過程中最深的一點(diǎn)體會(huì)是模型復(fù)雜不等于模型正確很多看起來(lái)“高大上”的多模態(tài)控制跑出來(lái)的曲線不一定比一個(gè)結(jié)構(gòu)清晰的簡(jiǎn)化模型更有說服力。做微電網(wǎng)仿真最好的路徑是先想透系統(tǒng)功率流的方向再把每個(gè)模塊的物理邊界條件算清楚最后才去Simulink里搭建和調(diào)參。遇到仿真發(fā)散先別急著加阻尼回頭看看是不是初始條件、電容電感參數(shù)或者模式切換邏輯本身出了問題。這套流程走完你的Simulink模型不僅是一個(gè)能出圖的仿真工程更是一套可以用來(lái)指導(dǎo)實(shí)際系統(tǒng)設(shè)計(jì)的可靠工具。