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