我做電力系統(tǒng)仿真這些年一個(gè)特別深的感受是微電網(wǎng)模型好搭但搭一個(gè)能說(shuō)服人的微電網(wǎng)模型很難。很多人一開(kāi)始圖省事用手畫(huà)三五節(jié)點(diǎn)的輻射狀拓?fù)淠靡粋€(gè)等效電源當(dāng)大電網(wǎng)跑通是跑通了可一到論文評(píng)審、項(xiàng)目驗(yàn)收別人問(wèn)一句你這個(gè)網(wǎng)絡(luò)結(jié)構(gòu)能不能反映真實(shí)系統(tǒng)的運(yùn)行特性立刻就卡住了。這次要聊的復(fù)合微電網(wǎng)模型基礎(chǔ)骨架選的是IEEE 14節(jié)點(diǎn)標(biāo)準(zhǔn)模型在Simulink里集成了柴油發(fā)電機(jī)、光伏陣列、電池儲(chǔ)能系統(tǒng)、電弧爐非線性負(fù)載。它最大的價(jià)值不在于單個(gè)模塊有多復(fù)雜而在于把標(biāo)準(zhǔn)輸電網(wǎng)骨架和分布式電源非線性負(fù)載組合到了一個(gè)仿真環(huán)境里既保留了標(biāo)準(zhǔn)算例的通用性和說(shuō)服力又覆蓋了微電網(wǎng)的核心特征。這篇文章我會(huì)從選型邏輯、各單元建模細(xì)節(jié)、電弧爐的坑、系統(tǒng)集成與工況設(shè)計(jì)、調(diào)試排錯(cuò)幾個(gè)維度展開(kāi)全程基于Simulink/Simscape Electrical環(huán)境適合正在做微電網(wǎng)仿真、分布式發(fā)電研究或者準(zhǔn)備把微電網(wǎng)當(dāng)作畢設(shè)/課題方向的讀者參考。1. 為什么微電網(wǎng)仿真要套IEEE 14節(jié)點(diǎn)從能跑到可信的選型邏輯1.1 標(biāo)準(zhǔn)節(jié)點(diǎn)模型在微電網(wǎng)仿真里的真正價(jià)值很多初學(xué)者會(huì)問(wèn)一個(gè)問(wèn)題微電網(wǎng)本身規(guī)模不大為什么要用IEEE 14節(jié)點(diǎn)這種輸電網(wǎng)級(jí)別的標(biāo)準(zhǔn)模型直接用一個(gè)小型輻射狀網(wǎng)絡(luò)不是更簡(jiǎn)單嗎這個(gè)問(wèn)題的答案得從**仿真結(jié)果的可信度**說(shuō)起。IEEE 14節(jié)點(diǎn)是IEEE發(fā)布的公開(kāi)標(biāo)準(zhǔn)算例拓?fù)浣Y(jié)構(gòu)、支路阻抗、變壓器變比、負(fù)荷分布全部有據(jù)可查全球研究者都在用同一個(gè)基準(zhǔn)做對(duì)比。你在這個(gè)骨架上做的控制策略、接入方案、非線性負(fù)載影響分析別人可以復(fù)現(xiàn)、可以對(duì)比、可以引用。這是自定義輻射狀網(wǎng)絡(luò)做不到的——你的拓?fù)浜蛥?shù)完全是自己定的別人無(wú)法橫向驗(yàn)證結(jié)論的說(shuō)服力天然矮一截。另外從技術(shù)角度講IEEE 14節(jié)點(diǎn)帶有環(huán)網(wǎng)結(jié)構(gòu)、多電壓等級(jí)和并聯(lián)支路比普通輻射狀微電網(wǎng)復(fù)雜得多。環(huán)網(wǎng)意味著功率流動(dòng)方向不是單一的電壓支撐和潮流分布需要真正算出來(lái)而不是看出來(lái)。這對(duì)微電網(wǎng)的控制策略設(shè)計(jì)是一個(gè)更嚴(yán)格也更真實(shí)的考驗(yàn)。用通俗的話說(shuō)你在簡(jiǎn)單網(wǎng)絡(luò)上驗(yàn)證的策略到了環(huán)網(wǎng)上不一定成立反過(guò)來(lái)能在14節(jié)點(diǎn)上成立的策略放到簡(jiǎn)單網(wǎng)絡(luò)通常沒(méi)問(wèn)題。1.2 IEEE 14節(jié)點(diǎn)關(guān)鍵參數(shù)與微電網(wǎng)改造思路IEEE 14節(jié)點(diǎn)標(biāo)準(zhǔn)模型的基本參數(shù)如下表所示這套數(shù)據(jù)在Matlab/Simulink的powerlib、Matpower等工具里都有現(xiàn)成算例可以直接導(dǎo)入?yún)⒖肌?shù)項(xiàng)數(shù)值說(shuō)明母線數(shù)量14包含1條平衡母線、多條PQ母線支路數(shù)量20含線路和變壓器支路同步發(fā)電機(jī)5臺(tái)接于母線1、2、3、6、8基準(zhǔn)功率100 MVA標(biāo)幺值系統(tǒng)統(tǒng)一基準(zhǔn)基準(zhǔn)電壓等級(jí)138kV / 13.8kV等原始輸電網(wǎng)參數(shù)無(wú)功補(bǔ)償并聯(lián)電容器組母線3、6、8配置直接用原始參數(shù)搭微電網(wǎng)有個(gè)問(wèn)題——輸電網(wǎng)的電壓等級(jí)和微電網(wǎng)不匹配。所以實(shí)際建模時(shí)要做兩點(diǎn)改造電壓等級(jí)縮放把微電網(wǎng)的基準(zhǔn)電壓降到中壓配電網(wǎng)級(jí)別常見(jiàn)做法是設(shè)為10kV或20kV與微電網(wǎng)工程實(shí)際吻合。注意改基準(zhǔn)電壓后支路阻抗和負(fù)荷功率要同步做標(biāo)幺值換算否則潮流結(jié)果完全失真。發(fā)電機(jī)節(jié)點(diǎn)改造原模型里母線1是平衡節(jié)點(diǎn)由同步發(fā)電機(jī)承擔(dān)系統(tǒng)功率平衡改造時(shí)把部分同步發(fā)電機(jī)替換為柴油發(fā)電機(jī)組的同步電機(jī)模型把光伏和儲(chǔ)能接到指定母線上。平衡節(jié)點(diǎn)的支撐角色可以由柴油發(fā)電機(jī)和電池儲(chǔ)能聯(lián)合承擔(dān)或者根據(jù)運(yùn)行模式切換。我的經(jīng)驗(yàn)是改造時(shí)先保留原模型的網(wǎng)絡(luò)骨架和負(fù)荷分布只替換發(fā)電源類型跑通之后再逐步調(diào)整控制策略。一上來(lái)就大改拓?fù)涑隽藛?wèn)題你都分不清是網(wǎng)絡(luò)參數(shù)錯(cuò)了還是控制不對(duì)。2. 四大核心單元在Simulink里的落地方式我把模型拆開(kāi)看微電網(wǎng)核心就四個(gè)部分柴油發(fā)電機(jī)、光伏、電池儲(chǔ)能、電弧爐。前三者是發(fā)電源/儲(chǔ)能電弧爐是非線性負(fù)載。它們的建模難度和側(cè)重點(diǎn)完全不同下面逐個(gè)過(guò)一遍。2.1 柴油發(fā)電機(jī)組骨架里的主力電源柴油發(fā)電機(jī)組由柴油機(jī)、同步發(fā)電機(jī)、調(diào)速器和勵(lì)磁系統(tǒng)四部分組成。在Simulink里模塊選型我建議這樣柴油機(jī)模型Simscape Electrical中沒(méi)有直接的柴油機(jī)模塊常規(guī)做法是用傳輸延遲Transport Delay加一階慣性環(huán)節(jié)模擬柴油機(jī)的油門響應(yīng)特性。點(diǎn)火延遲通常設(shè)為20~50ms對(duì)應(yīng)600~1500rpm下的幾個(gè)工作循環(huán)慣性時(shí)間常數(shù)按照機(jī)組轉(zhuǎn)動(dòng)慣量折算。同步發(fā)電機(jī)用標(biāo)準(zhǔn)的同步電機(jī)模型按機(jī)組容量設(shè)置定子電阻、電樞電感、勵(lì)磁參數(shù)。14節(jié)點(diǎn)模型各發(fā)電機(jī)容量不同母線1上的機(jī)組最大微網(wǎng)改造時(shí)柴油發(fā)電機(jī)組容量通常選擇2~5MVA。調(diào)速器Governor實(shí)現(xiàn)有功-頻率下垂控制即當(dāng)系統(tǒng)頻率下降時(shí)調(diào)速器增加油門開(kāi)度提高有功輸出。下垂系數(shù)一般設(shè)在4%~5%也就是頻率偏差1%對(duì)應(yīng)有功變化20%~25%。勵(lì)磁系統(tǒng)AVR維持機(jī)端電壓使用PI控制器目標(biāo)是無(wú)功與電壓的協(xié)調(diào)穩(wěn)定。有一個(gè)常見(jiàn)誤區(qū)要提醒很多人覺(jué)得柴油機(jī)就是一個(gè)恒功率源直接用可控電壓源代替。這在純穩(wěn)態(tài)分析里勉強(qiáng)能用但一旦涉及孤島運(yùn)行、負(fù)荷突變你就能看到頻率跌落、有功分配不均這些真實(shí)物理現(xiàn)象。柴油機(jī)的慣性響應(yīng)和調(diào)速器動(dòng)態(tài)是整個(gè)微電網(wǎng)頻率支撐的核心省不得。2.2 光伏模型MPPT和逆變器控制各管一攤光伏部分在Simulink里分三級(jí)PV陣列我用的是Simscape Electrical里的PV Array模塊可以直接設(shè)置串并聯(lián)數(shù)目。仿真環(huán)境里的光照強(qiáng)度W/m2和溫度是可變輸入方便做光照突變工況測(cè)試。在選擇PV面板參數(shù)時(shí)用廠商數(shù)據(jù)比如SunPower SPR-315E可以省去反復(fù)調(diào)參數(shù)的時(shí)間。Boost變換器與MPPT光伏輸出電壓一般遠(yuǎn)低于直流母線電壓需要Boost升壓。MPPT算法我用的擾動(dòng)觀察法PO步長(zhǎng)按仿真速度取舍——步長(zhǎng)太大會(huì)在最大功率點(diǎn)附近持續(xù)振蕩太小則追蹤速度慢。實(shí)現(xiàn)方式是通過(guò)S函數(shù)或MATLAB Function模塊把算法寫(xiě)進(jìn)去。逆變器與并網(wǎng)控制逆變器接受直流母線電能通過(guò)PWM調(diào)制輸出三相交流。并網(wǎng)控制采用PQ控制有功/無(wú)功解耦控制是目前微電網(wǎng)中最常規(guī)的并網(wǎng)策略——光伏按最大功率輸出有功無(wú)功按需求設(shè)定不參與頻率和電壓支撐。這里背后邏輯是光伏出力天然波動(dòng)讓光伏承擔(dān)調(diào)頻調(diào)壓任務(wù)反而會(huì)降低系統(tǒng)可靠性所以并網(wǎng)模式下優(yōu)先讓光伏最大化利用太陽(yáng)能把支撐任務(wù)交給柴油機(jī)和儲(chǔ)能。這也是微電網(wǎng)控制領(lǐng)域一個(gè)基本共識(shí)。2.3 電池儲(chǔ)能系統(tǒng)模式切換的關(guān)鍵角色電池儲(chǔ)能在微電網(wǎng)里的角色取決于系統(tǒng)處在并網(wǎng)模式還是離網(wǎng)模式并網(wǎng)模式儲(chǔ)能采用PQ控制此時(shí)它相當(dāng)于一個(gè)可調(diào)度的有功/無(wú)功源可以平滑光伏波動(dòng)、參與削峰填谷。離網(wǎng)模式儲(chǔ)能切換為V/f控制恒壓恒頻承擔(dān)起建立系統(tǒng)電壓和頻率基準(zhǔn)的任務(wù)。這個(gè)瞬間的控制模式切換是整個(gè)仿真里最容易出問(wèn)題的地方之一我后面會(huì)專門講調(diào)試經(jīng)驗(yàn)。電池本體我用的Simscape Electrical里的鋰電池模塊設(shè)置容量、SOC初值、內(nèi)阻等參數(shù)。雙向DC-DC變換器用來(lái)控制充放電電池側(cè)電流環(huán)加功率外環(huán)形成雙閉環(huán)。儲(chǔ)能容量要根據(jù)微電網(wǎng)負(fù)荷水平和光伏滲透率來(lái)確定不能隨便拍腦袋。一個(gè)粗略的估算邏輯是儲(chǔ)能額定功率應(yīng)覆蓋光伏最大波動(dòng)量的1.2~1.5倍能量容量則應(yīng)支撐系統(tǒng)在柴油機(jī)故障后能維持至少15~30分鐘的關(guān)鍵負(fù)荷。2.4 電弧爐非線性負(fù)載必須單獨(dú)說(shuō)的存在電弧爐是整個(gè)模型里最容易翻車的部分原因是它非線性太強(qiáng)。電弧爐的電阻不是固定的而是隨電弧電流和電弧長(zhǎng)度劇烈變化導(dǎo)致它從電網(wǎng)吸收的電流波形嚴(yán)重畸變產(chǎn)生大量諧波和電壓閃變。Simulink里搭建電弧爐主要有三種方式我在下表做了對(duì)比建模方式實(shí)現(xiàn)難度諧波還原度建議場(chǎng)景時(shí)變電阻模型低一般快速演示、整體系統(tǒng)級(jí)仿真電弧黑盒模型Cassie/Mayr中較高需要分析電弧動(dòng)態(tài)特性時(shí)電弧爐精煉爐組合模型高高完整工況模擬、電能質(zhì)量深入研究我在這個(gè)模型里采用的是帶諧波畸變的時(shí)變非線性電阻方案。它比用純電阻模擬要真實(shí)得多電弧電壓-電流特性存在明顯的非線性區(qū)每半個(gè)周波內(nèi)電弧都要經(jīng)歷引燃→穩(wěn)定燃燒→熄滅的過(guò)程。隨著燃弧過(guò)程的變化電阻值按特定規(guī)律周期性變化正是這種周期性變化導(dǎo)致了電壓跌落、諧波畸變等問(wèn)題。3. 電弧爐非線性負(fù)載這個(gè)模型最值得深挖的坑3.1 電弧爐到底在仿真里動(dòng)了什么電弧爐的非線性通俗講就是它不是一個(gè)線性電阻而是一個(gè)隨著電流時(shí)刻變化的可變電阻。電弧的物理特性決定了它的伏安特性曲線是負(fù)阻特性——電弧電流增大時(shí)電弧電阻反而下降。這個(gè)負(fù)阻特性讓系統(tǒng)的電流波形產(chǎn)生嚴(yán)重畸變不再是我們熟悉的按正弦規(guī)律變化的理想波形。這種畸變?cè)陬l域上表現(xiàn)為大量的諧波分量主要是2~7次諧波這里特指以工頻50Hz為基波的各次諧波。諧波電流注入系統(tǒng)后會(huì)在系統(tǒng)阻抗上產(chǎn)生諧波電壓導(dǎo)致母線電壓發(fā)生畸變電壓波形變得不平滑。如果你在Simulink里用FFT工具分析電弧爐接入母線的電壓電流會(huì)很直觀地看到50Hz基波旁邊多出一系列諧波尖峰。在實(shí)際的微電網(wǎng)運(yùn)行中這種非線性負(fù)載帶來(lái)的主要危害是電能質(zhì)量惡化諧波導(dǎo)致電壓/電流波形畸變率THD上升影響敏感負(fù)荷正常工作。設(shè)備損耗增加諧波電流使變壓器、電纜等設(shè)備額外發(fā)熱縮短絕緣壽命?？刂聘蓴_諧波可能干擾微電網(wǎng)中逆變器的鎖相環(huán)PLL和控制器導(dǎo)致控制精度下降甚至誤動(dòng)作。保護(hù)誤動(dòng)諧波可能引發(fā)過(guò)流保護(hù)、接地保護(hù)等誤動(dòng)作影響系統(tǒng)可靠性。3.2 建模方式選型我為什么推薦模塊微分方程的組合方案搭建電弧爐模型我推薦的思路是模塊組合微分方程的方法用Simscape Electrical里的基本元件搭一個(gè)電弧等效電路電弧的等效電阻和電抗用微分方程描述其動(dòng)態(tài)特性。具體實(shí)現(xiàn)上將電弧的動(dòng)態(tài)方程寫(xiě)入MATLAB Function模塊輸出一個(gè)實(shí)時(shí)變化的電阻值然后通過(guò)可控電阻模塊接入電路。這樣電弧電阻不是固定值而是隨電流/電壓狀態(tài)實(shí)時(shí)變化實(shí)現(xiàn)強(qiáng)非線性。微分方程采用簡(jiǎn)化的電弧模型核心關(guān)系是電弧電阻的變化率與當(dāng)前電阻、電弧功率損耗和電弧時(shí)間常數(shù)相關(guān)。電弧時(shí)間常數(shù)這個(gè)參數(shù)很關(guān)鍵它決定了電弧對(duì)電流變化的響應(yīng)速度通常取值在微秒到毫秒量級(jí)。時(shí)間常數(shù)越小電弧的動(dòng)態(tài)響應(yīng)越快諧波含量越豐富但仿真步長(zhǎng)要求也越高。我在模型里取的是0.5ms左右既能還原諧波特性的主要特征又不會(huì)讓仿真速度慢到無(wú)法接受。3.3 接入點(diǎn)和容量選擇非線性負(fù)載對(duì)系統(tǒng)影響的放大器電弧爐接在哪個(gè)母線、容量多大效果天差地別。原因在于非線性負(fù)載產(chǎn)生的諧波電流會(huì)沿著阻抗最小的路徑流動(dòng)而不同節(jié)點(diǎn)的系統(tǒng)阻抗不同。我的處理邏輯是這樣的母線選擇電弧爐作為工業(yè)負(fù)載通常接在負(fù)荷較集中的母線同時(shí)考慮該母線的短路容量。短路容量越大系統(tǒng)承受非線性負(fù)載的能力越強(qiáng)電壓畸變?cè)叫　Ｋ匀绻阆肟吹矫黠@的諧波影響和電壓波動(dòng)就應(yīng)該選短路容量相對(duì)較小的節(jié)點(diǎn)接入。容量設(shè)置電弧爐容量與系統(tǒng)基準(zhǔn)容量100MVA之間的比例要選好。我試過(guò)從3MVA到20MVA多組參數(shù)3MVA時(shí)畸變不明顯15MVA以上就會(huì)有顯著的電壓波形畸變。理想的仿真觀察區(qū)間是8~12MVA既能看出非線性影響又不會(huì)讓系統(tǒng)徹底崩潰。另外電弧爐接入的升壓變壓器要考慮漏抗參數(shù)。變壓器漏抗對(duì)電弧爐回路有兩個(gè)作用一是限制短路電流二是與電弧非線性特性相互作用形成電壓畸變。漏抗太小短路電流過(guò)大數(shù)值仿真容易出現(xiàn)震蕩漏抗太大諧波電流被過(guò)度抑制又看不到非線性效果。我一般取8%~12%的短路阻抗。3.4 電弧爐帶來(lái)的連鎖反應(yīng)對(duì)整個(gè)微電網(wǎng)的擾動(dòng)測(cè)試電弧爐在Simulink里不只影響它自己所在的母線。諧波通過(guò)線路傳播會(huì)影響到整個(gè)微電網(wǎng)的電壓質(zhì)量。我實(shí)測(cè)下來(lái)主要有三個(gè)連鎖反應(yīng)PCC公共連接點(diǎn)電壓畸變率升高即使電弧爐不接在PCC諧波也會(huì)通過(guò)線路傳導(dǎo)到PCC導(dǎo)致并網(wǎng)點(diǎn)的電能質(zhì)量指標(biāo)惡化。這在做并網(wǎng)符合性評(píng)估時(shí)是必須關(guān)注的指標(biāo)。逆變器交互影響光伏和儲(chǔ)能的逆變器采用PLL鎖相諧波電壓會(huì)干擾PLL的鎖定精度導(dǎo)致逆變器輸出電流附加波形畸變。這種諧波-控制-諧波的閉環(huán)交互是電能質(zhì)量惡化的重要推手。功率振蕩電弧爐的快速波動(dòng)會(huì)引起系統(tǒng)有功和無(wú)功的同步波動(dòng)在弱網(wǎng)條件下可能誘發(fā)功率振蕩。觀察柴油發(fā)電機(jī)的出力曲線你會(huì)看到頻率波動(dòng)明顯加大。做微電網(wǎng)仿真的人常常只關(guān)注源-荷平衡和控制策略但引入電弧爐這類非線性負(fù)載后電能質(zhì)量變成了不可忽視的約束條件。這也是為什么這個(gè)模型值得認(rèn)真對(duì)待。4. 系統(tǒng)集成、控制協(xié)調(diào)與運(yùn)行工況設(shè)計(jì)4.1 整體架構(gòu)四個(gè)單元怎么接成一個(gè)微電網(wǎng)在Simulink里的整體架構(gòu)我有兩條經(jīng)驗(yàn)可以分享一是電源接入方式。柴油發(fā)電機(jī)作為同步電源直接并接在交流母線上光伏和儲(chǔ)能通過(guò)電力電子變換器Boost逆變器/Buck-Boost逆變器并接指定母線上電弧爐通過(guò)電弧爐變壓器接到特定負(fù)荷母線。并網(wǎng)開(kāi)關(guān)將微電網(wǎng)與大電網(wǎng)連接。二是控制架構(gòu)。微電網(wǎng)的經(jīng)典控制架構(gòu)分為三層底層是設(shè)備級(jí)控制本地控制器中間層是微電網(wǎng)中央控制MGCC頂層是能量管理系統(tǒng)EMS。在仿真模型里底層控制是核心——柴油機(jī)的調(diào)速器和勵(lì)磁、光伏的MPPT和PQ控制、儲(chǔ)能的V/f或PQ控制都屬于底層。MGCC和EMS在模型中簡(jiǎn)化為運(yùn)行模式切換和功率指令分配。架構(gòu)上最容易犯的錯(cuò)誤是電源類型與母線性質(zhì)不匹配。比如把逆變器型電源直接接到?jīng)]有同步電源支撐的孤立母線上電壓和頻率沒(méi)有參考仿真必然發(fā)散。必須保證任意一個(gè)孤島子網(wǎng)內(nèi)至少有一個(gè)能提供電壓和頻率支撐的電源柴油發(fā)電機(jī)或儲(chǔ)能V/f控制。4.2 三套運(yùn)行工況仿真模型的價(jià)值要看它能測(cè)什么一個(gè)只有穩(wěn)態(tài)運(yùn)行曲線的微電網(wǎng)模型說(shuō)服力是不夠的。我建議至少設(shè)計(jì)三套運(yùn)行工況來(lái)驗(yàn)證模型工況一并網(wǎng)模式下光伏波動(dòng)光照突變?cè)O(shè)定光照強(qiáng)度在t2s時(shí)從1000W/m2突降到400W/m2觀察光伏出力驟降后柴油發(fā)電機(jī)和儲(chǔ)能如何響應(yīng)。預(yù)期結(jié)果是柴油機(jī)功率增加、儲(chǔ)能放電系統(tǒng)母線電壓維持在允許范圍內(nèi)。這個(gè)工況驗(yàn)證的是功率平衡調(diào)節(jié)能力。工況二孤島切換瞬間微電網(wǎng)從大電網(wǎng)解列t3s時(shí)斷開(kāi)并網(wǎng)開(kāi)關(guān)微電網(wǎng)轉(zhuǎn)為孤島運(yùn)行。在這個(gè)瞬間系統(tǒng)失去大電網(wǎng)支撐儲(chǔ)能必須從PQ控制切換到V/f控制。如果切換邏輯有問(wèn)題系統(tǒng)頻率和電壓會(huì)出現(xiàn)嚴(yán)重的暫態(tài)振蕩甚至崩潰。這個(gè)工況驗(yàn)證的是控制策略魯棒性和模式切換平滑性。工況三電弧爐非線性負(fù)載投切t4s時(shí)投入電弧爐負(fù)載觀察非線性負(fù)載對(duì)母線電壓、系統(tǒng)頻率、柴油發(fā)電機(jī)出力的影響。投入前和投入后分別用Powergui的FFT工具分析電壓波形和THD含量量化非線性負(fù)載對(duì)系統(tǒng)電能質(zhì)量的惡化程度。這個(gè)工況驗(yàn)證的是非線性源荷交互特性。4.3 仿真參數(shù)配置步長(zhǎng)和算法選不好模型跑不動(dòng)也跑不準(zhǔn)Simulink的仿真參數(shù)設(shè)置是很多人后期卡殼的重災(zāi)區(qū)。我在這套模型里的配置如下求解器ode23t變步長(zhǎng)梯形法適合電力電子和剛性系統(tǒng)混合的模型。ode45在接入電力電子變換器后常常效率很低。最大步長(zhǎng)設(shè)為1e-4秒0.1ms考慮到逆變器PWM開(kāi)關(guān)頻率通常在10kHz量級(jí)最大步長(zhǎng)必須能分辨PWM的開(kāi)關(guān)動(dòng)作。如果只關(guān)心機(jī)電暫態(tài)和能量管理可以適當(dāng)放到1e-3秒但涉及電弧爐時(shí)建議收窄。仿真時(shí)長(zhǎng)根據(jù)工況需求設(shè)置穩(wěn)態(tài)分析1~5秒即可暫態(tài)分析和孤島切換需要10秒以上才能看到完整動(dòng)態(tài)過(guò)程。Powergui設(shè)置選擇Continuous模式如果你用Phasor模式雖然仿真速度快但無(wú)法看到諧波波形電弧爐的非線性特性就白搭了。這里要特別說(shuō)明步長(zhǎng)和精度是矛盾的步長(zhǎng)越小精度越高但時(shí)間成本成倍增長(zhǎng)。我的做法是分階段調(diào)試時(shí)用大步長(zhǎng)快速驗(yàn)證邏輯最后跑正式工況再切小步長(zhǎng)。5. 調(diào)試階段的高頻報(bào)錯(cuò)與實(shí)測(cè)優(yōu)化經(jīng)驗(yàn)5.1 代數(shù)環(huán)問(wèn)題電弧爐模型最常見(jiàn)的翻車點(diǎn)我第一次把電弧爐的動(dòng)態(tài)模型接入電氣系統(tǒng)時(shí)仿真直接報(bào)Algebraic Loop或者干脆卡死不動(dòng)。原因很典型電弧電阻的計(jì)算依賴于當(dāng)前電流而電流的計(jì)算又依賴于電弧電阻二者形成代數(shù)閉環(huán)。解決這個(gè)問(wèn)題我用過(guò)的有效方案有三種優(yōu)先推薦第一種在反饋回路中插入Memory模塊或Unit Delay模塊用上一步的電弧電阻值計(jì)算當(dāng)前電流打破代數(shù)環(huán)。代價(jià)是電弧動(dòng)態(tài)引入了一步延遲但步長(zhǎng)足夠小時(shí)影響微弱。用Simscape物理建模方式替代信號(hào)流方式Simscape的物理連接天然支持微分方程描述不容易出現(xiàn)代數(shù)環(huán)。改用離散求解器離散求解器天然切斷代數(shù)環(huán)但對(duì)非線性系統(tǒng)的精度有損失。5.2 初始化失敗與變壓器勵(lì)磁涌流另一個(gè)高頻報(bào)錯(cuò)是初始化階段就發(fā)散或者仿真剛開(kāi)始的0.1秒內(nèi)電流尖峰異常大。這個(gè)現(xiàn)象十有八九是變壓器空載合閘勵(lì)磁涌流。原因是變壓器的鐵芯磁通在初始時(shí)刻不能突變但模型初始化時(shí)沒(méi)有設(shè)置正確的初始磁通導(dǎo)致合閘瞬間出現(xiàn)很大的勵(lì)磁電流。處理方法在變壓器模塊參數(shù)中啟用初始磁通設(shè)置讓變壓器在仿真起點(diǎn)就處于穩(wěn)態(tài)磁通狀態(tài)?；蛘咴谀Ｐ椭麟娐返臄嗦菲鏖]合時(shí)間上做文章起步階段先讓部分負(fù)荷退出等系統(tǒng)穩(wěn)定后再逐步投入。5.3 步長(zhǎng)與PWM開(kāi)關(guān)頻率的配合問(wèn)題電力電子變換器的PWM開(kāi)關(guān)頻率高通常5~20kHz如果步長(zhǎng)選太大PWM波形就擬合不出來(lái)仿真出的波形全是毛刺。如果步長(zhǎng)選太小仿真速度又太慢跑完一個(gè)10秒的工況要等半天。我的做法是把模型拆成機(jī)電暫態(tài)和電磁暫態(tài)兩個(gè)視角分別對(duì)待。做整體控制策略研究時(shí)用開(kāi)關(guān)頻率平均模型Average Model替代詳細(xì)的PWM模型速度可以提高一個(gè)數(shù)量級(jí)。做諧波分析時(shí)再用詳細(xì)PWM模型。Simulink里逆變器的平均值模型就是一個(gè)受控電壓源用調(diào)制信號(hào)直接控制輸出電壓省去開(kāi)關(guān)過(guò)程。5.4 功率不平衡導(dǎo)致的系統(tǒng)發(fā)散排查有一次我調(diào)參數(shù)時(shí)發(fā)現(xiàn)仿真總是跑著跑著頻率越偏越大最后發(fā)散。排查到最后發(fā)現(xiàn)是光伏和儲(chǔ)能的容量設(shè)置超出了系統(tǒng)負(fù)荷太多導(dǎo)致柴油機(jī)無(wú)法在調(diào)速器調(diào)節(jié)范圍內(nèi)抵消多余功率。功率平衡是微電網(wǎng)仿真的底線約束。仿真前先做靜態(tài)潮流校核柴油機(jī)出力光伏最大出力儲(chǔ)能最大放電功率要大于系統(tǒng)最大負(fù)荷同時(shí)各母線電壓在標(biāo)幺值0.95~1.05之間。這一點(diǎn)在改造IEEE 14節(jié)點(diǎn)模型時(shí)尤其要注意——原模型是輸電網(wǎng)算例負(fù)荷分布和微電網(wǎng)場(chǎng)景不完全匹配改造后功率平衡需要重新校核。5.5 調(diào)試順序的建議根據(jù)我用這套模型的經(jīng)驗(yàn)調(diào)參順序往往決定了調(diào)試效率。如果你一上來(lái)就全系統(tǒng)投入任何一個(gè)小錯(cuò)誤都會(huì)被放大成發(fā)散很難定位。我的調(diào)試順序是先單獨(dú)跑通柴油發(fā)電機(jī)帶基礎(chǔ)負(fù)荷驗(yàn)證網(wǎng)絡(luò)骨架和發(fā)電機(jī)控制正常。再加入光伏觀察并網(wǎng)控制和功率平衡。再加入儲(chǔ)能驗(yàn)證模式切換邏輯。最后加入電弧爐調(diào)整非線性負(fù)載參數(shù)。每加一級(jí)都用FFT工具和Scope觀察關(guān)鍵母線電壓電流波形確認(rèn)沒(méi)有異常再進(jìn)入下一級(jí)。這套方法看著笨實(shí)際是最省時(shí)間的。5.6 常見(jiàn)報(bào)錯(cuò)與排查對(duì)照表最后整理一份我在調(diào)試過(guò)程中遇到的高頻報(bào)錯(cuò)表方便大家直接對(duì)照排查報(bào)錯(cuò)/問(wèn)題現(xiàn)象核心原因排查方向Algebraic Loop報(bào)錯(cuò)電弧爐/逆變器控制存在代數(shù)閉環(huán)插入Memory/Unit Delay改用Simscape物理建模初始0.01s發(fā)散變壓器初始磁通未設(shè)置設(shè)置變壓器初始磁通斷路器延遲合閘波形毛刺嚴(yán)重最大步長(zhǎng)過(guò)大無(wú)法分辨PWM開(kāi)關(guān)步長(zhǎng)調(diào)到1e-4秒以下或使用平均值模型頻率一直下跌系統(tǒng)總有功不足校核負(fù)載功率與電源功率平衡頻率振蕩發(fā)散調(diào)速器/下垂參數(shù)過(guò)大降低下垂系數(shù)檢查PI參數(shù)穩(wěn)定性電壓波形嚴(yán)重畸變電弧爐容量過(guò)大或接入點(diǎn)短路容量過(guò)小調(diào)整電弧爐容量更換接入母線孤島切換失敗儲(chǔ)能V/f切換邏輯時(shí)序錯(cuò)誤先切換控制模式再斷開(kāi)并網(wǎng)開(kāi)關(guān)最后再說(shuō)一點(diǎn)我個(gè)人的實(shí)操體會(huì)這套復(fù)合微電網(wǎng)模型調(diào)試的重點(diǎn)其實(shí)不在電氣參數(shù)本身而在于控制策略與仿真環(huán)境的配合。很多人花大量時(shí)間調(diào)變壓器漏抗、調(diào)線路阻抗效果卻不明顯真正影響仿真成敗的往往是步長(zhǎng)設(shè)置、代數(shù)環(huán)處理、初始化狀態(tài)這些仿真工程層面的細(xì)節(jié)。另外仿真模型的真實(shí)感很大程度上取決于你對(duì)物理過(guò)程的敬畏——比如電弧爐的非線性特性不要為了省事就用線性電阻替代那樣出來(lái)的結(jié)果不僅平滑得不可信還會(huì)讓你完全錯(cuò)過(guò)電能質(zhì)量分析這個(gè)微電網(wǎng)研究的重要方向。