簡介本資源是一套面向電力系統(tǒng)仿真初學(xué)者與新能源方向研究生的雙饋風(fēng)力發(fā)電機DFIG與混合儲能系統(tǒng)聯(lián)合建模實戰(zhàn)資料聚焦風(fēng)電并網(wǎng)穩(wěn)定性提升這一工程痛點提供從原理理解到Simulink實操的完整閉環(huán)。壓縮包共77個文件含3個核心Simulink模型.mdl、9個MATLAB數(shù)據(jù)文件.mat、13個HTML格式建模說明文檔、9個C語言代碼文件用于RTW代碼生成及1個PDF版《雙饋風(fēng)機仿真詳細說明研二》全面覆蓋電機建模、變流器控制、多源儲能協(xié)調(diào)策略與實時仿真部署等關(guān)鍵環(huán)節(jié)4MB體積精煉實用無冗余素材。已有743人學(xué)習(xí)下載配套文檔逐模塊解析系統(tǒng)架構(gòu)、參數(shù)設(shè)定邏輯與仿真結(jié)果判據(jù)并附帶bat腳本與r2016b專屬編譯配置顯著降低MATLAB 2016b環(huán)境下模型編譯與運行門檻。 雙饋風(fēng)機加混合儲能這套仿真我前前后后折騰了大半個月從搭第一版單風(fēng)機模型到最后并網(wǎng)功率波動被削平到讓人滿意的程度中間踩的坑比想象中多得多。這個項目本身并不算罕見但真正能讓人照著搭起來、理解透的完整案例其實不多尤其是“雙饋風(fēng)機混合儲能”這個組合大多數(shù)人要么只做過風(fēng)機側(cè)要么只做過儲能側(cè)合在一起的時候又是另一回事。這篇文章就圍繞這個仿真項目展開把我搭建過程中的方案選擇、模型細節(jié)、參數(shù)計算、踩坑記錄和排查思路全部拆開來講適合正在做風(fēng)電并網(wǎng)仿真、儲能平抑波動研究或者準(zhǔn)備復(fù)現(xiàn)這個方向的畢業(yè)設(shè)計/課題的同學(xué)參考。1. 項目整體設(shè)計與方案選型思路1.1 為什么選雙饋風(fēng)機而不是直驅(qū)風(fēng)機風(fēng)力發(fā)電仿真里風(fēng)機模型主要就兩種路線雙饋異步風(fēng)機DFIG和永磁直驅(qū)風(fēng)機PMSG。在咱們這個項目場景下選雙饋風(fēng)機不是隨手拍腦袋而是有考量的。雙饋風(fēng)機的核心特點是轉(zhuǎn)子側(cè)通過背靠背變流器接入電網(wǎng)定子側(cè)直接并網(wǎng)。變流器容量只需要承擔(dān)轉(zhuǎn)差功率通常是風(fēng)機額定功率的30%左右這意味著變流器的成本、損耗都比直驅(qū)全功率變流器低很多。在仿真里這個特性也帶來一個挺有意思的建模差異——你不需要像直驅(qū)那樣把整個定子側(cè)功率都經(jīng)過變流器控制策略的重心落在轉(zhuǎn)子勵磁的調(diào)節(jié)上。從學(xué)術(shù)和工程價值角度說雙饋風(fēng)機到現(xiàn)在仍是陸上風(fēng)電的主流裝機機型之一研究它的并網(wǎng)特性、故障穿越和功率控制比只做直驅(qū)更有現(xiàn)實意義。在Simulink里搭DFIG模型工具箱里有現(xiàn)成的“Wind Farm (DFIG)”示例但是直接用現(xiàn)成模型和自己搭一套的區(qū)別在于——現(xiàn)成模型是黑盒參數(shù)改起來麻煩尤其你要在這個基礎(chǔ)上加儲能協(xié)調(diào)控制很多內(nèi)部信號根本引不出來。所以我最后選的是用Simscape Electrical Specialized Power Systems里的雙饋異步電機本體模型加上外部風(fēng)機氣動模型和轉(zhuǎn)子/網(wǎng)側(cè)變流器控制這樣內(nèi)部dq軸電流、磁鏈、轉(zhuǎn)矩信號全部可見可改加混合儲能的時候才靈活。1.2 混合儲能方案電池超級電容的互補邏輯儲能這塊很多人的第一反應(yīng)是用電池一組打天下但真做過風(fēng)電波動平抑仿真的話你會發(fā)現(xiàn)純電池方案挺難受的。電池的能量密度高能扛住持續(xù)性的功率缺額但它響應(yīng)速度慢循環(huán)壽命受頻繁充放電影響很大而風(fēng)電功率波動恰恰是既有持續(xù)性的低頻分量又有陣風(fēng)引起的瞬時高頻分量。超級電容的特性正好補上電池的短板功率密度大、響應(yīng)速度毫秒級、循環(huán)壽命幾十萬次但能量密度小扛不住長時間持續(xù)出力。所以在混合儲能系統(tǒng)里主流思路就是讓兩者打配合低頻波動分量給電池高頻瞬時分量給超級電容。這個邏輯用一句大白話說就是“大馬拉長活快馬接急活”。具體到Simulink實現(xiàn)我采用的是功率分頻協(xié)調(diào)控制策略。核心思路是先把并網(wǎng)功率參考值P_ref通過一階低通濾波器得出一條平滑的功率基準(zhǔn)線然后把風(fēng)機原始輸出功率與基準(zhǔn)線的差值作為儲能系統(tǒng)需要補償?shù)目偣β手噶钤儆靡粋€高通/低通分解結(jié)構(gòu)把總指令拆成電池和超級電容各自的出力指令。這個結(jié)構(gòu)不復(fù)雜但效果很直觀后面第3章會給出具體參數(shù)和實現(xiàn)步驟。1.3 系統(tǒng)拓?fù)渑c控制框架總覽整個仿真系統(tǒng)的拓?fù)浣Y(jié)構(gòu)是這樣的風(fēng)輪機通過傳動軸驅(qū)動DFIG轉(zhuǎn)子DFIG定子直接接入690V交流母線轉(zhuǎn)子通過背靠背變流器RSC轉(zhuǎn)子側(cè)變流器 GSC網(wǎng)側(cè)變流器接入同一個交流母線混合儲能系統(tǒng)通過雙向DC/DC變換器掛在直流母線上再通過一個三相逆變器并聯(lián)到交流母線。交流母線再通過升壓變壓器和輸電線路接入無窮大電網(wǎng)?？刂茖用娣譃槿龑语L(fēng)機層風(fēng)機氣動模型Cp-λ曲線槳距角控制RSC/GSC矢量控制實現(xiàn)最大功率跟蹤和直流母線電壓穩(wěn)定儲能層電池和超級電容各自的充放電控制雙向DC/DC電流內(nèi)環(huán)功率外環(huán)SOC管理與限幅保護協(xié)調(diào)層功率分頻分配算法負(fù)責(zé)計算儲能總出力指令并分解到兩種儲能單元。這個框架的好處是模塊之間耦合度低我調(diào)試風(fēng)機的時候可以把儲能模塊整個斷掉調(diào)試儲能的時候也可以把風(fēng)機輸出固定成一組實測風(fēng)速數(shù)據(jù)方便對照。仿真項目做到后面你會發(fā)現(xiàn)模塊化設(shè)計的價值在排障時體現(xiàn)得淋漓盡致。2. 核心模塊建模細節(jié)剖析2.1 雙饋風(fēng)機氣動與機械傳動建模雙饋風(fēng)機的氣動模型是整個仿真的源頭風(fēng)速進來之后風(fēng)輪捕獲的機械功率是P_m 0.5 · ρ · π · R2 · C_p(λ, β) · v3公式本身不難難在C_p函數(shù)的擬合。C_p是葉尖速比λ和槳距角β的函數(shù)很多論文里直接給一個經(jīng)驗公式C_p 0.22 · (116/λ? ? 0.4·β ? 5) · exp(?12.5/λ?)其中 1/λ? 1/(λ 0.08β) ? 0.035/(β3 1)在Simulink里實現(xiàn)這個函數(shù)我建議用Fcn模塊或者MATLAB Function模塊寫直接按公式搭就行。做的時候要注意一點λ的計算要用風(fēng)速和風(fēng)輪轉(zhuǎn)速的實時值風(fēng)速是輸入風(fēng)輪轉(zhuǎn)速來自發(fā)電機反饋兩者之間在啟動階段容易算出很大的λ值需要加一個飽和限幅不然C_p曲線數(shù)值會發(fā)瘋仿真直接發(fā)散。傳動系統(tǒng)我用的是兩質(zhì)量塊模型高速軸和低速軸分開建模中間用阻尼和剛度系數(shù)連接。這樣做的好處是能觀察軸系扭振現(xiàn)象尤其是在風(fēng)速突變時兩質(zhì)量塊模型比把傳動系統(tǒng)簡化為剛性單質(zhì)量塊更真實對后續(xù)儲能平抑功率波動時的動態(tài)響應(yīng)評估也更可靠。2.2 DFIG電機本體與轉(zhuǎn)子側(cè)控制電機本體我直接用的Simscape Electrical里的Asynchronous Machine Doubly-Fed模型參數(shù)就按1.5MW典型機組設(shè)置額定功率1.5MW、定子電壓690V、頻率50Hz、極對數(shù)2轉(zhuǎn)子側(cè)參數(shù)按折算后的標(biāo)幺值給定。轉(zhuǎn)子側(cè)變流器RSC的控制目標(biāo)是實現(xiàn)最大功率跟蹤和機組無功調(diào)節(jié)?？刂平Y(jié)構(gòu)是典型的雙閉環(huán)外環(huán)功率環(huán)有功功率指令來自MPPT查表無功功率指令按需求給定內(nèi)環(huán)轉(zhuǎn)子dq軸電流環(huán)。轉(zhuǎn)子d軸電流控制有功q軸電流控制無功這個分配關(guān)系跟傳統(tǒng)同步電機正好相反剛上手的人很容易搞混我也在這上面栽過跟頭。關(guān)鍵點在定子磁鏈定向。轉(zhuǎn)子側(cè)控制需要把轉(zhuǎn)子電流變換到以定子磁鏈為定向的dq坐標(biāo)系里所以必須用PLL鎖相環(huán)實時檢測定子電壓相位。我之前試過直接用工頻50Hz的固定相位代替鎖相電網(wǎng)電壓理想的時候問題不大但只要電壓有點跌落或者相位抖動控制就崩了。后來老老實實把PLL加上穩(wěn)定多了。鎖相環(huán)參數(shù)我用的是MATLAB自帶的PLL模塊默認(rèn)值微調(diào)比例系數(shù)0.2積分系數(shù)8跟蹤500Hz以內(nèi)的相位波動都沒問題。2.3 儲能單元建模電池與超級電容的差異化處理電池模型在Simulink里有好幾種選擇我推薦用Simscape Battery工具箱里的等效電路模型帶SOC、溫度特性、開路電壓?SOC曲線這三個維度。如果沒用這個工具箱用SimPowerSystems老版里的Battery模型也行把額定容量、初始SOC、內(nèi)阻等參數(shù)設(shè)好就能跑。比較多人踩坑的地方是電池容量設(shè)置。按照平抑功率波動的需求電池需要承擔(dān)的補償能量大約是分鐘級別的假設(shè)平均補償功率為100kW、持續(xù)平抑5分鐘那電池可用容量至少需要E P·t 100kW × 300s ≈ 8.33kWh考慮SOC不能從0干到100%實際運行一般在20%~90%區(qū)間有效容量系數(shù)取0.7所以電池額定容量至少配置到E_rated 8.33 / 0.7 ≈ 11.9kWh這個數(shù)量級已經(jīng)超過1C充放電倍率約束了所以還要反推電池功率等級。如果SOC初始設(shè)定在50%允許最大放電深度70%額定容量12kWh的電池最多能放出8.4kWh平均放電功率能達到100kW的話就需要電池能承受約12C的短時倍率這顯然不現(xiàn)實。最后我把電池容量提升到50kWh這樣短時功率倍率在2C左右靠譜多了。仿真里不用真算壽命但參數(shù)別離譜否則審稿人一眼就能看出來。超級電容相對簡單我用的是串聯(lián)RC等效電路模型電容容值C15F等效串聯(lián)內(nèi)阻R0.01Ω初始電壓設(shè)為250V。超級電容的端電壓在充放電過程中是線性變化的這在變換器控制里要注意它的電壓環(huán)沒法像電池那樣恒定所以超級電容側(cè)DC/DC一般采用電流控制外加電壓上限保護避免過充。2.4 雙向DC/DC變換器與協(xié)調(diào)控制實現(xiàn)儲能單元要掛到直流母線上必須有雙向DC/DC變換器做電壓變換和功率方向切換。這里用的是雙向Buck-Boost拓?fù)潆姵貍?cè)低壓、母線側(cè)高壓升壓模式放電、降壓模式充電絕緣柵雙極晶體管用PWM驅(qū)動頻率我設(shè)在5kHz可以兼顧精度和速度?？刂撇呗允请娏鲀?nèi)環(huán)功率外環(huán)。外環(huán)接收協(xié)調(diào)層給出的有功指令P_bat和P_sc除以各自端電壓得到電流參考值I_ref經(jīng)過PI調(diào)節(jié)器輸出占空比再進PWM發(fā)生器。這里的核心技巧是電流環(huán)的PI參數(shù)整定。我用工程整定法先根據(jù)電感L和開關(guān)周期估算穿越頻率再調(diào)Kp和Ki。電流環(huán)Kp初始值取0.5Ki取20實際調(diào)試時逐步增大到振蕩臨界值的0.6倍作為最終參數(shù)效果很好。協(xié)調(diào)控制層的功率分頻邏輯用兩個一階低通濾波器實現(xiàn)。設(shè)總儲能出力指令為P_store那么P_bat P_store / (T_bat·s 1)P_sc P_store ? P_bat濾波器時間常數(shù)T_bat的選擇影響分頻點頻率。我用的T_bat30s對應(yīng)分頻頻率大約是0.005Hz也就是說周期大于200秒的低頻分量歸電池管高頻快速波動歸超級電容這樣一個組合下來并網(wǎng)功率曲線既平滑又不會讓電池頻繁動作。3. 實操過程與關(guān)鍵步驟實現(xiàn)3.1 從零搭建Simulink仿真模型的完整流程下面這一步一步的操作流程是我實際搭建時反復(fù)調(diào)整后固化下來的你照著走基本能跑通。第一步新建模型并配置仿真環(huán)境。Simulink模型里先設(shè)置求解器為變步長、ode23tb算法最大步長0.001s相對誤差1e-4。這是電力電子仿真比較穩(wěn)的組合既能捕捉PWM開關(guān)細節(jié)又不至于太慢。電源頻率50Hz、采樣時間基值用離散化后的1e-5s做PWM波形生成。第二步搭風(fēng)機氣動模型。從Simulink庫拖入Constant風(fēng)速輸入、FcnCp函數(shù)、Product、Gain等模塊先算出機械轉(zhuǎn)矩輸出到DFIG電機的輸入端。建議把Cp-λ曲線先畫出來驗算一下最大Cp應(yīng)該能到0.438左右達不到說明公式或參數(shù)有問題。第三步搭DFIG電機與背靠背變流器主回路。電機模型放到SimPowerSystems域里轉(zhuǎn)子側(cè)和網(wǎng)側(cè)各配一個兩電平VSC變流器直流母線電容22000μF母線電壓預(yù)設(shè)1150V690V交流對應(yīng)直流母線電壓的典型值。變流器要用Universal Bridge模塊搭并配好PWM脈沖輸入口。第四步搭RSC和GSC控制器?？刂破髟谄胀⊿imulink域里搭用電壓/電流測量模塊把信號引入經(jīng)過坐標(biāo)變換、PI調(diào)節(jié)、PWM發(fā)生再輸出到變流器開關(guān)管。這一步是最費時間的建議先從空載調(diào)通變流器再接電機。第五步搭混合儲能系統(tǒng)。電池和超級電容各接一個雙向Buck-Boost變換器變換器高壓側(cè)掛在直流母線上低壓側(cè)接儲能單元?？刂菩盘杹碜詤f(xié)調(diào)層P_bat和P_sc指令。第六步搭協(xié)調(diào)控制層。采集風(fēng)機輸出功率P_wind通過低通濾波器得到P_smooth用P_storeP_smooth?P_wind得到儲能總指令再經(jīng)分頻分解成P_bat和P_sc送出去。第七步加電網(wǎng)與測量模塊。交流母線經(jīng)0.69/30kV升壓變壓器、30km輸電線路接無窮大電源并網(wǎng)點加Three-Phase V-I Measurement模塊測功率、電壓、電流數(shù)據(jù)用To Workspace導(dǎo)出方便后續(xù)畫圖。3.2 關(guān)鍵參數(shù)一覽表親測可用配置參數(shù)名數(shù)值說明風(fēng)機額定功率1.5 MWDFIG容量定子額定電壓690 V線電壓有效值直流母線電壓1150 V背靠背變流器公共母線直流母線電容22000 μF母線穩(wěn)壓電容電流環(huán)Kp/Ki0.5 / 20RSC和GSC電流環(huán)初值功率環(huán)Kp/Ki2 / 20外環(huán)功率調(diào)節(jié)低通濾波時間常數(shù)T_smooth10 s并網(wǎng)功率平滑基準(zhǔn)分頻時間常數(shù)T_bat30 s電池低頻通道開關(guān)頻率5 kHz全部變流器PWM電池額定容量50 kWh磷酸鐵鋰等效模型超級電容容值15 FRC簡化模型這套參數(shù)我不能保證你直接抄作業(yè)就完美因為控制器參數(shù)跟模型版本、求解器設(shè)置都有關(guān)系但作為起點非常夠用大部分情況微調(diào)Kp就能收斂。3.3 如何記錄和分析仿真結(jié)果仿真數(shù)據(jù)的記錄方法直接影響你分析問題的效率。用Scope看波形現(xiàn)場調(diào)試還行但要做后處理和生成圖片就麻煩了我建議所有關(guān)鍵信號都并行接一份To Workspace變量名用結(jié)構(gòu)體帶時間的形式保存。風(fēng)速序列、風(fēng)機輸出功率、儲能總指令、電池功率、超級電容功率、并網(wǎng)點功率這六個信號必存有條件還可以把SOC、直流母線電壓、轉(zhuǎn)子dq電流都存下來。在MATLAB腳本里把Simulink模型跑完再用plot命令把這些變量畫在一張圖里對比比在Scope里截圖直觀太多。比如你要驗證混合儲能的平抑效果就畫三根曲線原始風(fēng)機功率、平滑參考值、實際并網(wǎng)點功率一眼就能看出高頻毛刺是不是被超級電容吃掉了持續(xù)偏差是不是被電池補上了。數(shù)據(jù)分析時還有一個指標(biāo)可以算并網(wǎng)功率波動率。定義是最大功率?最小功率/額定功率用平滑前后的兩組數(shù)據(jù)分別算就能量化儲能介入的效果。4. 常見問題與排查技巧實錄4.1 仿真發(fā)散和超調(diào)問題我碰到最多的bug就是仿真跑著跑著直接發(fā)散或者某些變量出現(xiàn)NaN。排查順序建議是這樣的第一步檢查初始化。DFIG電機模型的初始狀態(tài)必須與工況匹配如果風(fēng)速突然從5m/s跳到12m/s而風(fēng)機轉(zhuǎn)速還在低速區(qū)軸系扭矩會瞬間爆表。解決方法是給風(fēng)速加斜坡變化或者把風(fēng)機初始轉(zhuǎn)速設(shè)置到匹配風(fēng)速的穩(wěn)態(tài)值。第二步檢查PI參數(shù)。電流環(huán)Kp過大會導(dǎo)致控制量震蕩發(fā)散表現(xiàn)為PWM占空比在上下限之間來回跳。我遇到這種情況時先把Kp降一半看看波形是否平復(fù)再按0.5~0.8倍臨界值回去這個方法屢試不爽。第三步檢查代數(shù)環(huán)。Simulink里如果在控制回路中出現(xiàn)直通路徑會報代數(shù)環(huán)錯誤表現(xiàn)為仿真卡死或結(jié)果異常。解決方法是在測量信號后面加一個很小的延遲比如1e-6s的Memory模塊把代數(shù)環(huán)打破。這個技巧在雙閉環(huán)控制里特別常用。4.2 超級電容SOC和電壓越限超級電容容值小特性類似一個大電容最容易出現(xiàn)的問題就是端電壓越限。當(dāng)風(fēng)電功率劇烈波動時高頻補償功率指令可能超出超級電容的承載能力電壓沖到上限或跌到下限。我的處理策略是在協(xié)調(diào)控制層加限制環(huán)節(jié)超級電容的功率指令通過一個飽和模塊限幅限幅值根據(jù)當(dāng)前端電壓實時計算電壓越接近上限允許的充電功率越小越接近下限允許的放電功率越小。這個邏輯用MATLAB Function寫也就十幾行代碼但加了之后整個系統(tǒng)魯棒性提升非常明顯。4.3 電池SOC管理策略電池SOC的問題主要是仿真跑久了SOC漂移出合理區(qū)間然后整個儲能系統(tǒng)輸出能力失效?；A(chǔ)的解決手段是SOC滯環(huán)控制SOC低于20%時限制放電功率高于90%時限制充電功率中間區(qū)間正常工作。還有一個容易被忽略的點是初始SOC的設(shè)置。如果你的仿真場景是中午風(fēng)速持續(xù)偏高、風(fēng)機出力大于平滑參考值那儲能大部分時間在充電初始SOC就不要設(shè)到80%否則跑一會兒就頂?shù)缴舷蘖?。反過來風(fēng)速持續(xù)偏低時初始SOC不要低于40%。這個要靠對場景的理解來預(yù)判沒有萬能初值。4.4 仿真速度優(yōu)化變步長求解器加PWM開關(guān)切換仿真速度經(jīng)常慢得讓人抓狂。我的經(jīng)驗是以下三個優(yōu)化手段疊加使用速度能提升3~5倍第一把PWM開關(guān)頻率從5kHz降到2kHz動態(tài)特性差別不大但仿真步數(shù)顯著減少第二把測量環(huán)節(jié)的采樣時間從連續(xù)改成離散設(shè)置合適的采樣步長比如20μs避免連續(xù)狀態(tài)過多第三在仿真參數(shù)里把“Signal logging”只保留關(guān)鍵信號不要勾選所有輸出信號不然內(nèi)存和計算都吃緊。如果你后面要做多工況掃參這套優(yōu)化是必須的。5. 建模說明文檔的整理思路這個項目標(biāo)題里專門提到“附詳細建模說明文檔”我就多聊幾句文檔怎么寫才有價值。很多人寫建模說明就是貼一堆模型截圖加參數(shù)表這其實沒什么用。有參考價值的建模說明文檔應(yīng)該包含四塊內(nèi)容第一塊是系統(tǒng)拓?fù)鋱D和控制框架圖讓人不看模型就知道你整體結(jié)構(gòu)是啥樣的第二塊是每個模塊的建模方法選擇和理由比如為什么電池用等效電路模型而不用電化學(xué)模型這體現(xiàn)建模者的思考深度第三塊是參數(shù)設(shè)置表和計算過程尤其那些需要通過公式算出來的參數(shù)必須給出從物理量到標(biāo)幺值或者實際值的推導(dǎo)步驟第四塊是調(diào)試過程中遇到的問題和解決辦法這部分是最珍貴的。在寫這個項目文檔時我把自己調(diào)試過程中的關(guān)鍵截圖、每次參數(shù)改動后的波形對比、踩坑記錄都整理進去了花了一整天時間專門整理文檔后來回頭看非常值。仿真模型的代碼和模型文件本身會過時但建模思路和排查經(jīng)驗不會過期這才是文檔的核心價值。最后再分享一個我自己的體會這個項目調(diào)通后后續(xù)擴展方向非常多。把風(fēng)機模型換成永磁直驅(qū)就是另一個課題把混合儲能的控制策略改成模糊控制或者模型預(yù)測控制可以寫論文把并網(wǎng)模型從無窮大電源改成實際電網(wǎng)阻抗模型可以做弱電網(wǎng)適應(yīng)性分析。也就是說這套仿真框架搭好了就是一個可以反復(fù)利用的研究平臺遠比單獨跑通一個模型有價值。本文還有配套的精品資源點擊獲取