做電力系統(tǒng)分析的同行應(yīng)該都有同感課本里把暫態(tài)穩(wěn)定講得明明白白等面積法則一畫(huà)、臨界切除時(shí)間一算好像就完事了?？烧嬉肧imulink把“單機(jī)無(wú)窮大系統(tǒng)暫態(tài)穩(wěn)定性仿真”完整跑起來(lái)盯著兩相接地短路故障時(shí)發(fā)電機(jī)轉(zhuǎn)速變化的波形一條條解讀往往會(huì)在模塊選型、參數(shù)配置、故障設(shè)置和波形判讀上反復(fù)折騰。這篇內(nèi)容就是我從零搭建這個(gè)仿真項(xiàng)目的全過(guò)程復(fù)盤(pán)把模型背后的物理依據(jù)、Simulink里每個(gè)關(guān)鍵模塊的配置方式、故障時(shí)序的設(shè)計(jì)邏輯、以及我實(shí)際調(diào)試中踩過(guò)的坑都掰開(kāi)講清楚。適合正在做電力系統(tǒng)暫態(tài)穩(wěn)定課程設(shè)計(jì)、畢業(yè)設(shè)計(jì)或者想用Simulink復(fù)現(xiàn)典型短路故障仿真的同學(xué)參考。坦白說(shuō)這類仿真在文獻(xiàn)里看起來(lái)就是一張系統(tǒng)圖加幾張波形圖但自己動(dòng)手建模時(shí)會(huì)發(fā)現(xiàn)任何一個(gè)參數(shù)不合理轉(zhuǎn)速波形都可能出現(xiàn)完全反物理的走勢(shì)。下面按我的實(shí)際操作順序來(lái)寫(xiě)。1. 為什么單機(jī)無(wú)窮大系統(tǒng)是研究暫態(tài)穩(wěn)定的最佳切入點(diǎn)1.1 暫態(tài)穩(wěn)定研究的物理本質(zhì)機(jī)械功率與電磁功率的失衡暫態(tài)穩(wěn)定研究的是系統(tǒng)在遭受大擾動(dòng)之后同步發(fā)電機(jī)能否繼續(xù)保持同步運(yùn)行的能力。這里的“大擾動(dòng)”最常見(jiàn)的就是短路故障。要理解整個(gè)仿真必須先抓住一個(gè)核心矛盾發(fā)電機(jī)轉(zhuǎn)子上同時(shí)存在兩個(gè)功率一個(gè)是原動(dòng)機(jī)輸入的機(jī)械功率Pm一個(gè)是發(fā)電機(jī)向系統(tǒng)輸出的電磁功率Pe。正常穩(wěn)態(tài)運(yùn)行時(shí)Pm和Pe嚴(yán)格相等轉(zhuǎn)子以同步轉(zhuǎn)速旋轉(zhuǎn)功角δ保持恒定。一旦線路上發(fā)生兩相接地短路網(wǎng)絡(luò)結(jié)構(gòu)發(fā)生變化發(fā)電機(jī)出口到無(wú)窮大母線的功率傳輸能力被嚴(yán)重削弱Pe會(huì)瞬間跌落。但原動(dòng)機(jī)側(cè)的汽門(mén)開(kāi)度或者水門(mén)開(kāi)度還來(lái)不及調(diào)整Pm基本維持原值。這就會(huì)出現(xiàn)“輸入多、輸出少”的局面多出來(lái)的功率全部轉(zhuǎn)化為轉(zhuǎn)子動(dòng)能發(fā)電機(jī)轉(zhuǎn)速開(kāi)始上升。轉(zhuǎn)速一旦偏離同步速功角就開(kāi)始擺動(dòng)系統(tǒng)能不能回到新的穩(wěn)定運(yùn)行點(diǎn)就是暫態(tài)穩(wěn)定要回答的問(wèn)題。這個(gè)邏輯是所有后續(xù)分析的根基。模型可以建得很復(fù)雜但你要盯住的永遠(yuǎn)是Pm和Pe之間的差值以及這個(gè)差值如何在時(shí)間軸上推動(dòng)轉(zhuǎn)速和功角變化。1.2 單機(jī)無(wú)窮大SMIB模型的簡(jiǎn)化邏輯與適用邊界單機(jī)無(wú)窮大系統(tǒng)由發(fā)電機(jī)、升壓變壓器、輸電線路通常是雙回線和無(wú)窮大母線構(gòu)成。無(wú)窮大母線指的是電壓幅值和頻率都恒定、可以吸收或發(fā)出任意功率的理想節(jié)點(diǎn)用它來(lái)等效外部大電網(wǎng)。引入SMIB模型最大的價(jià)值是把一個(gè)復(fù)雜的多節(jié)點(diǎn)電力網(wǎng)絡(luò)退化成一條單功率傳輸通路。發(fā)電機(jī)只有一個(gè)外部電網(wǎng)用一個(gè)恒壓源代表于是系統(tǒng)動(dòng)態(tài)的核心變量只剩下功角δ也就是發(fā)電機(jī)內(nèi)電勢(shì)相對(duì)無(wú)窮大母線電壓的相位差。教材里推導(dǎo)的功率特性方程 Pe(E·V/XΣ)·sinδ以及等面積法則都嚴(yán)格依賴這個(gè)模型。多機(jī)系統(tǒng)當(dāng)然更接近實(shí)際電網(wǎng)但多機(jī)系統(tǒng)的暫態(tài)穩(wěn)定本質(zhì)上是機(jī)群之間的相對(duì)功角搖擺問(wèn)題復(fù)雜度成倍上升。從學(xué)習(xí)角度不把單機(jī)模型的穩(wěn)態(tài)運(yùn)行點(diǎn)、故障期間功率特性變化、故障切除后的減速過(guò)程吃透直接上多機(jī)仿真只會(huì)一頭霧水。所以先做SMIB不是簡(jiǎn)化過(guò)度而是最為科學(xué)的遞進(jìn)路徑。1.3 為什么故障場(chǎng)景選兩相接地短路很多人做仿真時(shí)習(xí)慣一上來(lái)就搞三相短路因?yàn)槿喽搪纷顕?yán)重、最直觀。但工程實(shí)際中短路故障按發(fā)生概率從高到低依次是單相接地、兩相接地、兩相相間、三相短路。兩相接地屬于典型的不對(duì)稱短路在架空線路、電纜線路中都有較高的發(fā)生概率而且它比單相接地更嚴(yán)重又不像三相短路那樣把三相對(duì)稱性完全破壞分析邏輯剛好處于“需要?jiǎng)佑眯蚓W(wǎng)絡(luò)但又不會(huì)過(guò)于復(fù)雜”的中間位置。更關(guān)鍵的是兩相接地短路時(shí)故障點(diǎn)存在負(fù)序和零序分量。負(fù)序分量會(huì)在發(fā)電機(jī)轉(zhuǎn)子表面感應(yīng)出兩倍頻的負(fù)序磁場(chǎng)引起轉(zhuǎn)子附加損耗和發(fā)熱這在對(duì)稱的三相短路里看不到。所以用兩相接地做暫態(tài)穩(wěn)定分析不僅能看轉(zhuǎn)速搖擺行為還能順帶研究不對(duì)稱故障對(duì)發(fā)電機(jī)的額外影響信息量比三相短路大得多。2. Simulink模型搭建元件選型與參數(shù)配置全記錄2.1 用現(xiàn)有電力系統(tǒng)元件庫(kù)還是自己搭微分方程動(dòng)手之前先想清楚技術(shù)路線。Simulink里建電力系統(tǒng)暫態(tài)模型兩條路一是用Simscape Electrical庫(kù)里的Specialized Power Systems元件庫(kù)里面有現(xiàn)成的同步電機(jī)、變壓器、輸電線路、故障模塊二是完全用基礎(chǔ)Simulink模塊把轉(zhuǎn)子運(yùn)動(dòng)方程、派克方程、網(wǎng)絡(luò)方程全部手工搭建。我選擇的是元件庫(kù)路線。原因很直接這個(gè)項(xiàng)目的目標(biāo)是研究暫態(tài)穩(wěn)定性、分析轉(zhuǎn)速變化規(guī)律而不是重新發(fā)明電力系統(tǒng)模型。同步電機(jī)模塊內(nèi)部已經(jīng)封裝了定子繞組電磁暫態(tài)方程、轉(zhuǎn)子運(yùn)動(dòng)方程和磁鏈方程參數(shù)對(duì)話框里填上標(biāo)幺值就能工作省去大量推導(dǎo)和調(diào)試時(shí)間。元件庫(kù)的劣勢(shì)在于參數(shù)項(xiàng)比較多物理單位、標(biāo)幺制混在一起容易看花眼但只要理解了每個(gè)參數(shù)的含義反而比自己搭方程更可控。2.2 模型整體結(jié)構(gòu)從發(fā)電機(jī)到無(wú)窮大母線模型拓?fù)渚褪堑湫蚐MIB結(jié)構(gòu)按順序串聯(lián)以下模塊。模塊庫(kù)位置主要參數(shù)說(shuō)明同步發(fā)電機(jī)Specialized Power Systems / Machines200 MVA13.8 kV50 HzH4 s標(biāo)幺值凸極機(jī)模型升壓變壓器Specialized Power Systems / Elements200 MVA13.8/230 kV漏抗10%Yd11接線雙回輸電線路Specialized Power Systems / Elements100 km正序參數(shù)按典型值每回線各用一組PI Section無(wú)窮大母線Specialized Power Systems / Sources230 kV內(nèi)阻抗極小用Three-Phase Source等效故障模塊Specialized Power Systems / Elements兩相接地0.1 s投入0.2 s切除和線路并聯(lián)接入測(cè)量模塊Specialized Power Systems / Measurements電壓電流測(cè)量送入Scope觀察Powergui電力系統(tǒng)專用模塊Continuous仿真必須存在用于潮流初始化同步發(fā)電機(jī)選標(biāo)幺值模型輸出轉(zhuǎn)速直接是pu值1.0對(duì)應(yīng)同步轉(zhuǎn)速便于后續(xù)計(jì)算偏差。變壓器高壓側(cè)中性點(diǎn)接地方式要注意選擇接地是為了給零序電流提供回路這對(duì)兩相接地故障的仿真結(jié)果影響很大。輸電線路用PI Section模型足夠分布參數(shù)模型在100 km級(jí)別差異不大PI模型在暫態(tài)研究中是主流做法。無(wú)窮大母線直接用Three-Phase Source等效。內(nèi)部電感、電阻要設(shè)置得足夠小這樣外部系統(tǒng)才能表現(xiàn)出現(xiàn)電壓支撐能力。如果內(nèi)阻抗給大了母線電壓在故障期間會(huì)明顯跌落那就不是真正的無(wú)窮大母線了波形會(huì)失真。2.3 同步發(fā)電機(jī)參數(shù)設(shè)置的講究同步電機(jī)模塊參數(shù)對(duì)話框里的每一項(xiàng)都對(duì)應(yīng)真實(shí)物理量。項(xiàng)目里最需要反復(fù)調(diào)整的有三組慣性常數(shù)H、阻尼系數(shù)D、以及d/q軸電抗。慣性常數(shù)H直接決定轉(zhuǎn)子“輕重”標(biāo)準(zhǔn)定義是轉(zhuǎn)子在額定轉(zhuǎn)速下儲(chǔ)存的動(dòng)能與發(fā)電機(jī)額定容量之比。H越小轉(zhuǎn)子越輕同樣的功率差額引起的轉(zhuǎn)速變化率就越大暫態(tài)穩(wěn)定性越差。我在模型里取H4 s這是200 MVA級(jí)汽輪發(fā)電機(jī)的典型值域如果做靈敏度分析可以把它從2 s掃到9 s觀察轉(zhuǎn)速波形的上升斜率差異。阻尼系數(shù)D反映了轉(zhuǎn)子上的機(jī)械阻尼和電氣阻尼的綜合效果。有些模型允許直接填D有些則是通過(guò)阻尼繞組參數(shù)隱式體現(xiàn)。如果D填0故障切除后的轉(zhuǎn)速振蕩衰減會(huì)非常緩慢實(shí)際機(jī)組由于阻尼繞組、原動(dòng)機(jī)阻尼的存在不會(huì)這樣。建議首次仿真保留輕微阻尼等波形穩(wěn)定后再逐步調(diào)小觀察臨界狀態(tài)的失穩(wěn)模式。d軸、q軸電抗參數(shù)直接影響功率特性曲線的幅值。Xd1.5 pu左右、Xd0.35 pu左右、Xd0.25 pu左右是比較典型的凸極機(jī)參數(shù)。要注意的是暫態(tài)穩(wěn)定分析通常關(guān)心的是q軸暫態(tài)電抗Xd和d軸暫態(tài)時(shí)間常數(shù)Td0因?yàn)楣收习l(fā)生后電磁功率的恢復(fù)過(guò)程與它們強(qiáng)相關(guān)。3. 兩相接地短路故障的設(shè)置邏輯與序網(wǎng)分析3.1 Three-Phase Fault模塊里的兩相接地配置Specialized Power Systems里的Three-Phase Fault模塊是一個(gè)可投切的三相斷路器加接地開(kāi)關(guān)組合。要模擬兩相接地短路需要在Fault Type欄勾選Fault A和Fault B同時(shí)勾選Ground fault這樣A、B兩相在故障時(shí)刻同時(shí)接到地C相保持正常構(gòu)成典型的A-B兩相接地故障。時(shí)序參數(shù)Transition Times填[0.1 0.2]含義是0.1秒時(shí)故障投入0.2秒時(shí)故障切除。這里0.1秒之前的時(shí)段是系統(tǒng)空載/滿載穩(wěn)定運(yùn)行時(shí)間用來(lái)讓潮流計(jì)算和機(jī)電暫態(tài)進(jìn)入穩(wěn)態(tài)0.1到0.2秒之間模擬的是短路持續(xù)期0.2秒后的波形則是故障切除后的恢復(fù)和搖擺過(guò)程。0.1秒的故障持續(xù)時(shí)間對(duì)應(yīng)的是保護(hù)裝置的動(dòng)作時(shí)間現(xiàn)實(shí)中高壓線路主保護(hù)動(dòng)作時(shí)間在幾十毫秒到一百毫秒左右取0.1秒很典型。還有一個(gè)容易忽略的細(xì)節(jié)Fault Resistance不能填0。這個(gè)在第五章詳細(xì)說(shuō)這里只強(qiáng)調(diào)我用的是0.001 Ω既保證物理上接近金屬性短路又不會(huì)讓網(wǎng)絡(luò)方程出現(xiàn)數(shù)值奇異。3.2 序網(wǎng)絡(luò)視角兩相接地故障為什么會(huì)讓功率特性“塌”下去要理解兩相接地短路對(duì)發(fā)電機(jī)轉(zhuǎn)速的影響離不開(kāi)對(duì)稱分量法。A、B兩相接地后故障點(diǎn)在正序、負(fù)序、零序三個(gè)序網(wǎng)絡(luò)里同時(shí)存在邊界條件。對(duì)兩相接地短路正序、負(fù)序、零序網(wǎng)絡(luò)在故障點(diǎn)處是并聯(lián)關(guān)系忽略過(guò)渡電阻時(shí)先并聯(lián)再和故障阻抗串聯(lián)這決定了故障點(diǎn)的等效阻抗。從發(fā)電機(jī)端往系統(tǒng)看故障發(fā)生前的等效電抗是變壓器漏抗加上線路電抗故障發(fā)生后由于負(fù)序網(wǎng)絡(luò)和零序網(wǎng)絡(luò)的引入正序網(wǎng)絡(luò)中透入了一個(gè)附加阻抗使得從發(fā)電機(jī)看過(guò)去的系統(tǒng)總電抗顯著增大。功率特性方程Pe(E·V/XΣ)·sinδ中XΣ變大意味著同樣的功角δ下能傳輸?shù)碾姶殴β蚀蠓陆倒β侍匦郧€整體被壓低。這就是故障期間Pe驟降的根本原因。理解這一步就不需要死記結(jié)論只要看一下序網(wǎng)的并聯(lián)結(jié)構(gòu)就知道為什么兩相接地比三相短路的功率特性降得少一些但比正常運(yùn)行低得多。故障切除后負(fù)序零序通路消失XΣ恢復(fù)到故障前水平功率特性重新抬升減速面積開(kāi)始起作用。3.3 故障切除方式與保護(hù)時(shí)序的對(duì)應(yīng)關(guān)系故障切除在模型里有兩種做法結(jié)果差異很大。做法一是故障模塊在0.2秒斷開(kāi)但輸電線路仍然完整系統(tǒng)回到故障前的拓?fù)渥龇ǘ怯脭嗦菲髂K把故障線路整體跳開(kāi)相當(dāng)于永久性故障后的線路切除系統(tǒng)網(wǎng)絡(luò)結(jié)構(gòu)永久改變。我這次做得是第一種也就是等面積法則教材里的經(jīng)典設(shè)定故障點(diǎn)在0.2秒消失所有線路保持不變。這種情況下故障切除后的功率特性曲線完全恢復(fù)到故障前水平減速面積最大。如果要模擬“線路跳閘后單回線運(yùn)行”則需要用Breaker把其中一回線斷開(kāi)并考慮跳閘后XΣ增大、功率特性低于故障前暫態(tài)穩(wěn)定性會(huì)明顯變差。兩種設(shè)定沒(méi)有對(duì)錯(cuò)之分對(duì)應(yīng)的是不同的工程場(chǎng)景但解讀波形時(shí)一定要清楚自己用的是哪種。如果想把故障持續(xù)時(shí)間變成變量反復(fù)掃可以把Transition Times里的0.2改成工作區(qū)變量比如tf然后在腳本里用for循環(huán)調(diào)用sim命令批量仿真。我習(xí)慣的做法是先讓0.1秒作為固定故障投入時(shí)間然后把切除時(shí)間從0.11秒逐步增加到0.18秒觀察轉(zhuǎn)速波形的失穩(wěn)邊界這個(gè)思路第4章會(huì)展開(kāi)。4. 轉(zhuǎn)速波形的全過(guò)程解讀與穩(wěn)定性判據(jù)4.1 無(wú)故障基準(zhǔn)仿真先證明模型本身正確做故障仿真之前一定要先跑一遍無(wú)故障基準(zhǔn)工況。具體做法是把故障模塊的Transition Times改成[10 10]也就是仿真時(shí)間內(nèi)不投入或者干脆把故障模塊從通路中斷開(kāi)。此時(shí)觀察發(fā)電機(jī)轉(zhuǎn)速應(yīng)該是恒定的1.0 pu輸出功率和功角穩(wěn)定在由潮流計(jì)算決定的初始值上。這一步不是走過(guò)場(chǎng)而是排查模型問(wèn)題的最高效手段。我在這個(gè)環(huán)節(jié)遇到過(guò)功角緩慢漂移的情況后來(lái)發(fā)現(xiàn)是潮流初始化時(shí)發(fā)電機(jī)的機(jī)械功率初值和電磁功率沒(méi)對(duì)齊導(dǎo)致仿真一開(kāi)始轉(zhuǎn)子就有微小的不平衡轉(zhuǎn)矩?；鶞?zhǔn)工況能穩(wěn)定運(yùn)行后面故障仿真的所有異常才值得分析和解讀。4.2 故障發(fā)生后轉(zhuǎn)速變化的三個(gè)階段接上故障模塊后轉(zhuǎn)速波形可以清晰地分成三個(gè)階段來(lái)讀。第一階段是0到0.1秒的穩(wěn)態(tài)段轉(zhuǎn)速保持在1.0 pu附近波形是一條平滑直線。第二階段是0.1到0.2秒的故障段。故障瞬間電磁功率Pe突然下降機(jī)械功率Pm來(lái)不及變化凈加速轉(zhuǎn)矩出現(xiàn)轉(zhuǎn)速開(kāi)始從1.0 pu往上爬。爬升的斜率由兩個(gè)因素決定一是Pe跌落的深度也就是功率特性被壓低的程度二是慣性常數(shù)H的大小。H越小、加速功率越大波形就越陡。這個(gè)階段轉(zhuǎn)速不會(huì)立刻大幅跳變因?yàn)檗D(zhuǎn)子有慣性但0.1秒的故障時(shí)間里轉(zhuǎn)速已經(jīng)積累了明顯偏差。第三階段是0.2秒后的恢復(fù)段。故障切除后Pe恢復(fù)此時(shí)功角δ已經(jīng)比故障前增大了不少。如果此時(shí)電磁功率超過(guò)機(jī)械功率轉(zhuǎn)子開(kāi)始減速轉(zhuǎn)速先是一段上升慣性帶來(lái)的沖高然后回落越過(guò)同步速產(chǎn)生類似過(guò)阻尼/欠阻尼的往復(fù)振蕩最終在阻尼作用下穩(wěn)定在1.0 pu附近。如果故障時(shí)間過(guò)長(zhǎng)轉(zhuǎn)速則表現(xiàn)為單調(diào)上升不再回頭。讀這一段波形時(shí)要注意一個(gè)細(xì)節(jié)轉(zhuǎn)速曲線在故障投入和切除瞬間可能會(huì)出現(xiàn)非常窄的高頻毛刺這是電磁暫態(tài)和開(kāi)關(guān)切換引起的不是模型出錯(cuò)第5章會(huì)詳細(xì)說(shuō)怎么處理。4.3 轉(zhuǎn)速和功角兩種判斷穩(wěn)定性的角度轉(zhuǎn)速波形只能提供一個(gè)方向性判斷最終回到1.0 pu基本就是穩(wěn)定的持續(xù)偏離就是失穩(wěn)。但更精細(xì)的穩(wěn)定性邊界要靠功角δ的搖擺曲線來(lái)判斷。獲取功角的方式同步電機(jī)模塊輸出中有一個(gè)轉(zhuǎn)子角度信號(hào)另一個(gè)是無(wú)窮大母線電壓的相角兩者做差得到功角。有些版本的電機(jī)模塊直接輸出Machine Power Angle功率角可以直接引出。如果只能用電壓測(cè)量模塊獲取母線電壓相角需要注意相位是以哪個(gè)參考點(diǎn)為基準(zhǔn)的建議用PLL鎖相環(huán)鎖定母線電壓相位信號(hào)和機(jī)械轉(zhuǎn)子角信號(hào)一起送入示波器。功角判據(jù)的含義是SMIB系統(tǒng)在給定機(jī)械功率Pm下功率特性曲線與Pm水平線有兩個(gè)交點(diǎn)左側(cè)的是穩(wěn)定平衡點(diǎn)δs右側(cè)的是不穩(wěn)定平衡點(diǎn)δu。故障期間轉(zhuǎn)子加速功角向右擺動(dòng)如果功角越過(guò)了δu即使故障切除后電磁功率大于機(jī)械功率累計(jì)的減速面積也已經(jīng)不足以抵消加速面積系統(tǒng)必然失步。所以看功角曲線上是否逼近甚至穿越δu比單純看轉(zhuǎn)速更接近本質(zhì)。4.4 故障持續(xù)時(shí)間掃描估算臨界切除時(shí)間為了把這個(gè)仿真從定性變成定量可以做一組簡(jiǎn)單的參數(shù)掃描。固定故障投入時(shí)間0.1秒不變把切除時(shí)間從0.11秒逐步增加到0.18秒每次記錄轉(zhuǎn)速曲線的形態(tài)和功角最大值。故障切除時(shí)間轉(zhuǎn)速特征穩(wěn)定判斷0.11 s轉(zhuǎn)速小幅波動(dòng)后快速回到1.0 pu穩(wěn)定0.13 s振蕩明顯約2個(gè)周期后收斂穩(wěn)定0.15 s第一擺功角逼近不穩(wěn)定平衡點(diǎn)回落緩慢臨界0.17 s轉(zhuǎn)速持續(xù)上升不回到同步速失穩(wěn)在項(xiàng)目這套參數(shù)下臨界切除時(shí)間大致在0.14到0.16秒之間。這個(gè)數(shù)值和用等面積法手算的結(jié)果可以做對(duì)照。手動(dòng)算加速面積、減速面積應(yīng)該能落在同一個(gè)區(qū)間二者互相驗(yàn)證。這個(gè)過(guò)程做完對(duì)“臨界切除時(shí)間”這個(gè)概念就不再是背定義而是真正有了直觀認(rèn)識(shí)。5. 反復(fù)踩坑后的實(shí)用調(diào)試經(jīng)驗(yàn)5.1 初始化失敗多半在潮流計(jì)算這一步Simulink里電力系統(tǒng)模型運(yùn)行前Powergui會(huì)先做潮流初始化這一步經(jīng)常報(bào)錯(cuò)。最典型的情況是Infinity bus或Synchronous Machine的初值不給合理范圍潮流不收斂。我之前遇到過(guò)一次把發(fā)電機(jī)功率設(shè)為滿載、但無(wú)窮大母線電壓初值相位和變壓器聯(lián)結(jié)組號(hào)不匹配一運(yùn)行就初始化失敗。解決套路是先把發(fā)電機(jī)出力調(diào)小、把變壓器分接頭設(shè)為默認(rèn)讓系統(tǒng)在一個(gè)非?！拜p載”的狀態(tài)下完成初始化再逐步增加出力觀察穩(wěn)態(tài)波形。還有一個(gè)經(jīng)驗(yàn)線路長(zhǎng)度如果從100公里改成200公里末端電壓會(huì)明顯降低導(dǎo)致潮流初始化困難先短后長(zhǎng)的做法可以快速定位是線路還是電源的問(wèn)題。5.2 求解器與最大步長(zhǎng)不連續(xù)點(diǎn)怎么控制故障投入和切除本質(zhì)上是網(wǎng)絡(luò)拓?fù)涞耐蛔僑imulink里的變步長(zhǎng)求解器在開(kāi)關(guān)切換時(shí)刻容易出現(xiàn)步長(zhǎng)收縮和數(shù)值振蕩。我對(duì)比過(guò)幾種求解器配置最終穩(wěn)定使用的是ode23tb最大步長(zhǎng)限制在1e-4秒量級(jí)。如果用ode45這種顯式非剛性算法故障瞬間容易出現(xiàn)波形尖峰、甚至積分失敗。還有一點(diǎn)容易忽略Powergui的仿真類型選擇Continuous時(shí)不要同時(shí)在模型里混入離散模塊這會(huì)導(dǎo)致?tīng)顟B(tài)更新不同步。如果非用離散模式跑不可步長(zhǎng)選50微秒以內(nèi)比較穩(wěn)妥。轉(zhuǎn)速信號(hào)這種低速機(jī)電量用連續(xù)仿真是最省心的。5.3 故障電阻為什么不能填0不少第一次做故障仿真的同學(xué)會(huì)在Fault Resistance里填0結(jié)果一運(yùn)行就報(bào)矩陣奇異或者波形上出現(xiàn)超高幅值尖峰。原因很簡(jiǎn)單理想短路點(diǎn)阻抗為零故障支路電流趨向無(wú)窮大數(shù)值求解會(huì)碰到病態(tài)矩陣。工程上實(shí)際電弧總有一定電阻取0.001到0.01歐姆既不影響短路特性的物理規(guī)律又能保證數(shù)值穩(wěn)定。不要為了追求所謂“理想故障”而死磕0值這既不真實(shí)也不必要。5.4 轉(zhuǎn)速波形上的“毛刺”和數(shù)值振蕩的辨別和同步電機(jī)轉(zhuǎn)速信號(hào)打交道時(shí)會(huì)發(fā)現(xiàn)波形在故障投切瞬間出現(xiàn)一個(gè)極窄的高頻脈沖很多初學(xué)者把這個(gè)當(dāng)成了轉(zhuǎn)速真實(shí)突變甚至據(jù)此得出“發(fā)電機(jī)瞬間加速到失步”的錯(cuò)誤結(jié)論。這個(gè)脈沖主要來(lái)自兩部分一是開(kāi)關(guān)切換瞬間的電磁暫態(tài)二是數(shù)值積分在間斷點(diǎn)附近的暫態(tài)誤差。辨別方法很簡(jiǎn)單把最大步長(zhǎng)從1e-4改成1e-5重新仿真同一工況如果脈沖幅度變化不大說(shuō)明是真實(shí)的電磁暫態(tài)高頻分量如果幅度明顯變小那就是數(shù)值振蕩。機(jī)電暫態(tài)分析的頻帶在0.1赫茲到幾赫茲之間這些高頻分量本來(lái)就不在該研究的頻段內(nèi)可以在示波器前串聯(lián)一個(gè)一階慣性環(huán)節(jié)做低通濾波比如1/(0.01s1)就能看到干凈的轉(zhuǎn)速搖擺曲線。這套仿真完整走下來(lái)之后我最大的感受是課本上那些等面積法則、臨界切除時(shí)間第一次以“會(huì)變化的曲線”形式活生生出現(xiàn)在眼前。故障期間的加速、切除后的減速往復(fù)、越過(guò)臨界點(diǎn)后的失步過(guò)程全都在轉(zhuǎn)速波形上留下了清晰的痕跡。后來(lái)我又順手把故障類型改成三相短路、把雙回線改成單回線對(duì)比同一故障時(shí)長(zhǎng)下的轉(zhuǎn)速上升斜率對(duì)暫態(tài)穩(wěn)定的直觀理解比翻十遍書(shū)都管用。如果你也準(zhǔn)備做類似仿真別急著追求復(fù)雜的多機(jī)模型先把單機(jī)無(wú)窮大系統(tǒng)這一條鏈路跑透把波形里的每一個(gè)拐彎都解釋清楚再去擴(kuò)展你的系統(tǒng)規(guī)模。