：從暫態(tài)原理到Simulink仿真建模)
簡介基于MATLAB的三相短路分析資料包適合電力系統(tǒng)專業(yè)學生與工程技術(shù)人員用于掌握故障建模、仿真與結(jié)果分析的方法。資源從三相短路基本概念出發(fā)結(jié)合Simulink平臺完整演示了電力系統(tǒng)模型搭建、故障類型與位置設置、仿真運行及數(shù)據(jù)解讀的全過程并給出可復用的test.m腳本和詳細課設報告。壓縮包共12個文件大小僅1.09MB包含m腳本、PDF報告、4個txt數(shù)據(jù)文件與6張PNG結(jié)果圖txt文件保存了電路原始數(shù)據(jù)、節(jié)點阻抗/導納矩陣及輸出結(jié)果PNG圖清晰呈現(xiàn)發(fā)電機等效前后的矩陣變化便于對照報告逐項理解。目前已有2576人學習尤其適合正在完成電力系統(tǒng)課設或MATLAB仿真作業(yè)的讀者。通過研讀代碼與報告可以快速掌握用Simulink分析三相短路電流、驗證保護設備動作特性并能直接修改參數(shù)擴展至其他故障類型是一份實用的電力仿真學習資料。 做短路計算校驗的時候三相短路幾乎是每個電力系統(tǒng)從業(yè)者都躲不開的基本功。我最早接觸三相短路分析是在學暫態(tài)過程的時候當時全靠手算公式、查沖擊系數(shù)曲線算一條輻射線路的短路電流要折騰大半天后來真正參與設備選型和保護整定發(fā)現(xiàn)短路電流校核幾乎貫穿所有項目才開始系統(tǒng)性用MATLAB來搭仿真模型。MATLAB里的Simscape Electrical模塊老版本叫SimPowerSystems支持直接拖拽電源、變壓器、線路和故障模塊設置一個故障時間就能看到三相短路后的電壓電流波形比手算直觀太多。這篇把整套流程和踩過的坑整理出來給準備用MATLAB做三相短路分析的同學和工程師一個可以照做的參考。1. 三相短路分析的本質(zhì)算的是暫態(tài)電流不是穩(wěn)態(tài)電流1.1 為什么會短路短路電流為什么可怕三相短路是指供電系統(tǒng)中三相之間發(fā)生金屬性或非金屬性的直接連接相當于負載被短路掉回路阻抗驟減系統(tǒng)電壓被瞬間拉低短路電流可以高達額定電流的數(shù)倍甚至數(shù)十倍。為什么都拿三相短路作為最嚴重的故障工況來校核因為對稱短路時正序阻抗最小短路電流通常最大設備動穩(wěn)定、熱穩(wěn)定、斷路器遮斷容量、保護定值整定幾乎都要圍繞這個最大工況來做。真正讓三相短路分析變得棘手的是暫態(tài)過程。短路瞬間電流不能突變?yōu)榱司S持磁鏈守恒回路里會感應出一個非周期分量電流這個直流分量和周期分量疊加導致短路電流波形嚴重偏置。用公式來表達故障相的短路電流大約為i(t) sqrt(2) I * [sin(ωt α - φ) - sin(α - φ) * e^(-t/Ta)]其中I是次暫態(tài)電流有效值α是故障瞬間的電壓相位角φ是短路回路阻抗角Ta L/R為直流分量衰減時間常數(shù)。初始的沖擊電流峰值在半個周波附近出現(xiàn)50Hz系統(tǒng)約0.01s其幅值可以近似為 Kimp * sqrt(2) * IKimp就是沖擊系數(shù)取值通常在1.8左右。如果不理解這個暫態(tài)過程仿真跑出來的波形就只是“看到了一個波形”沒辦法從里面提取設備校核需要的關(guān)鍵參數(shù)。1.2 MATLAB做短路分析的兩條路線用MATLAB做三相短路分析我試過兩條路線各有適用場景。第一條是純腳本編程把系統(tǒng)等效成R-L阻抗用解析公式直接算短路電流周期分量、非周期分量和沖擊電流。好處是快速、透明適合簡單輻射網(wǎng)和批量參數(shù)掃描比如我要比較不同線路長度對短路電流的影響寫一個循環(huán)就出來了。缺點是處理復雜拓撲和多節(jié)點網(wǎng)絡時手寫節(jié)點阻抗矩陣非常費勁很容易出錯。第二條是用Simulink里的Simscape Electrical模塊搭建數(shù)學模型這也是我后期項目里用得最多的一條路。它把電源、變壓器、線路、負載、故障模塊做成可視化元件搭完模型就能看到各點電壓電流波形還能方便地切換故障類型、故障時刻和故障位置。缺點是模型自由度多求解器、步長、阻抗參數(shù)設置不當結(jié)果可能歪到離譜。我個人的建議是做教學演示或驗證理論公式優(yōu)先用腳本做工程系統(tǒng)分析、保護整定或需要給同事看波形時用Simulink建模。兩者結(jié)合是最好的先用腳本算一個基準值再拿Simulink去驗證和擴展。2. 建模前的參數(shù)整理從單線圖到Simulink元件參數(shù)2.1 系統(tǒng)阻抗與標幺值換算不管用什么工具三相短路分析的輸入基礎(chǔ)都是系統(tǒng)阻抗。模型里的每個元件都要有明確的R和X這些參數(shù)可以從設備銘牌、線路單位阻抗表和系統(tǒng)短路容量中獲取。Simulink里可以直接填歐姆Ω和亨H但我建議在做理論校核時先換算成標幺值這樣對比手算結(jié)果更方便。基準值的選取遵循 S_B 100 MVA或10 MVA基準電壓取電網(wǎng)額定電壓例如10kV系統(tǒng)選U_B 10.5kV基準阻抗為 Z_B U_B^2 / S_B。舉個例子我常用的一個驗證模型是10kV輻射網(wǎng)系統(tǒng)電源側(cè)短路容量為200MVA折算到10.5kV基準下的電源等效阻抗約 Z_source U_B^2 / S_k 10.5^2 / 200 0.551Ω。變壓器容量2MVA短路電壓百分比 U_k% 6%折算到10.5kV側(cè)變壓器阻抗約為 Z_T U_k% * U_B^2 / (100 * S_T) 6 * 10.5^2 / (100 * 2) 3.307Ω。線路用10km的10kV架空線每公里電抗0.35Ω電阻0.17Ω線路阻抗 Z_L 1.7 j3.5Ω。這些值最后填到Simulink元件里就是有名值。2.2 在Simscape Electrical庫中選型與搭建新版MATLAB的Simulink庫瀏覽器里路徑一般是Simscape Electrical Specialized Power Systems Fundamental Blocks。老版本直接叫SimPowerSystems位置類似。需要用到的基本元件有Three-Phase Source三相電壓源設置相電壓/線電壓有效值、頻率、內(nèi)阻和電感。Three-Phase Transformer (Two Windings)雙繞組三相變壓器設置額定容量、頻率、一次二次額定電壓、短路阻抗百分比。Three-Phase Series RLC Branch線路阻抗可以直接填R、L。Three-Phase Fault三相故障模塊可選故障類型和切換時間。Three-Phase V-I Measurement測量三相電壓電流輸出到Scope或Workspace。Powergui電力系統(tǒng)圖形用戶界面模塊選擇連續(xù)模式或相量模式。我一般選Continuous才能看到暫態(tài)波形。搭建順序不復雜三相電壓源接到變壓器一次側(cè)變壓器二次側(cè)經(jīng)過線路阻抗到故障模塊再接到負載。故障模塊放在短路點測量模塊放在電源側(cè)和故障點兩側(cè)各一個方便對比電壓跌落和電流變化。別忘了在模型中放置Powergui模塊否則整個模型會報錯。把上節(jié)計算出的參數(shù)填入元件對話框Three-Phase Source的Phase-to-phase RMS voltage填10.5kV如果電源側(cè)視為無窮大系統(tǒng)內(nèi)阻為0但為了仿真接近實際建議填0.5Ω和0.01H變壓器一次10.5kV、二次0.4kV也可以不降壓用于等值短路計算線路RLC填1.7Ω和0.35H對應10km注意電抗換算到電感需要除以2πf。2.3 故障模塊參數(shù)設置細節(jié)Three-Phase Fault模塊是所有仿真里的核心開關(guān)。雙擊模塊Fault type選Three-phase faultABC這時三相同時短路。Fault resistance填0.01Ω不要填0數(shù)值上可以更穩(wěn)定如果做金屬性短路填一個很小的值即可。切換方式我習慣用外部控制External control用Step模塊從0跳到1觸發(fā)故障模塊導通也可以用內(nèi)部定時直接設置Transition time為0.1s故障初始狀態(tài)選Open經(jīng)過0.1s后切換到Closed。這里有個容易被忽略的點Transition time是故障開始時間不是故障持續(xù)時間。如果希望故障持續(xù)0.2s后恢復需要設置兩個Transition time比如0.1和0.3s狀態(tài)依次為Open→Closed→Open。只填一個時間故障會一直保持到仿真結(jié)束。對于短路分析我們通常只關(guān)心故障后前幾個周波所以故障持續(xù)0.2s足夠。我還習慣在故障模塊的External Fault Breaker控制端口接一個Signal Generator或者Step這樣可以精確控制故障觸發(fā)相位。想研究短路初始角的影響就通過Step的延遲時間調(diào)整短路發(fā)生的合閘角。3. 仿真運行與波形解讀短路的真正信息藏在暫態(tài)過程里3.1 從波形中讀出周期分量和非周期分量模型搭好設置仿真時長0.4s最大步長0.1ms運行后雙擊Scope看到故障點電流波形。穩(wěn)定情況下看到故障前是規(guī)則正弦故障發(fā)生后電流跳變前半段波形上下嚴重不對稱包絡線一高一低隨后逐漸恢復到對稱正弦。這個不對稱就是非周期分量在衰減而對稱部分就是短路電流的周期分量。我在實際讀波形時會把數(shù)據(jù)導到MATLAB工作區(qū)再畫圖方便標注峰值。比如用To Workspace模塊輸出三相電流然后在命令窗口提取數(shù)據(jù)t out.tout; iA out.iabc(:,1); plot(t, iA); xline(0.1, --); grid on;如果仿真的故障開始時間是0.1s波形圖里0.1s之前的正弦就是負荷電流或空載電流0.1s之后才是短路電流。用數(shù)據(jù)游標讀出故障后第一個周期波峰大約15kA再往后幾個周波峰值回落到約8kA前者包含了直流分量后者才接近周期分量幅值。3.2 沖擊電流與衰減時間常數(shù)的實測方法沖擊電流是短路后第一個峰值是設備動穩(wěn)定校驗的核心參數(shù)。我在模型中讀取A相電流最大值為15kA左右與理論估算的沖擊電流對照。理論值怎么算先由短路回路總阻抗算出次暫態(tài)周期分量有效值I然后乘以sqrt(2)再乘以沖擊系數(shù)Kimp。Kimp可以用1.8的典型值也可以按Kimp 1 e^(-0.01/Ta)計算。T_a L/R (總電抗/2πf)/總電阻。如果L/R算下來是0.03sKimp 1 e^(-0.333) ≈ 1.717理論峰值約為 sqrt(2) * I * 1.717。衰減時間常數(shù)還能直接從波形里測出來把故障后的電流包絡線提取出來直流分量衰減到初始值的36.8%所需的時間就是Ta。由于非周期分量在示波器里表現(xiàn)為包絡線的不對稱差可以取上包絡線減下包絡線再除以2得到直流衰減曲線用對數(shù)坐標擬合得出時間常數(shù)。這個方法適合驗證實際模型中的回路時間常數(shù)是否與手算一致。3.3 不同故障時刻對波形的影響短路初始角對沖擊電流影響很大。這是因為公式里sin(α-φ)決定了非周期分量的初始幅值。如果短路瞬間電壓過零α≈0回路阻抗角φ約80度sin(α-φ)很大直流分量幾乎以最大初始值出現(xiàn)沖擊電流最大如果短路瞬間電壓剛好在峰值附近非周期分量初始值很小波形基本對稱沖擊系數(shù)接近1。用MATLAB模擬這個現(xiàn)象很簡單在Three-Phase Source模塊里設置Phase angle of phase A或者用外部Step控制故障觸發(fā)時間把故障時間從0.1s改成0.105s等觀察第一個峰值的不同。我實測發(fā)現(xiàn)同一個模型只改故障時刻沖擊電流能從14kA變到16kA差別超過了10%。所以在做保護整定時不能光看一組仿真結(jié)果要找最嚴重的合閘角。4. 理論手算與仿真結(jié)果交叉驗證差距從哪來4.1 手算短路電流的簡化模型對于輻射網(wǎng)手算三相短路電流的簡化步驟是忽略負荷把電源等效為次暫態(tài)電勢E短路點電壓在故障前按額定電壓考慮短路電流有效值I_sc c U_N / (sqrt(3) Z_total)。電壓系數(shù)c在高壓系統(tǒng)通常取1.1低壓系統(tǒng)取1.05。Z_total是短路點到電源的全部等值阻抗之和包括電源內(nèi)阻抗、變壓器阻抗和線路阻抗。我用前面10kV模型來算Z_source 0.551ΩZ_T 3.307ΩZ_L 1.7 j3.5Ω。要注意電抗和電阻要分別加再求模。總電阻R 0.5 0 1.7 2.2Ω電源電阻很小忽略或按前面填的0.5Ω總電抗X 0.551 3.307 3.5 7.358Ω假設電源內(nèi)阻電抗0.551。Z_total sqrt(2.2^2 7.358^2) ≈ 7.68Ω。取c 1.1短路電流有效值 I_sc 1.1 * 10.5k / (sqrt(3) * 7.68) ≈ 0.868kA 868A。這里因為模型用了10kV/0.4kV變壓器和0.4kV側(cè)故障阻抗折算會不同手算時應注意基準側(cè)選擇為了簡單這里是在10kV側(cè)等值短路點。把仿真結(jié)果的短路電流周期分量有效值從波形里提煉出來兩者通常能對上偏差在5%以內(nèi)。如果偏差超過10%不要急著說仿真不對先去檢查參數(shù)折算和基準值。4.2 仿真和手算不一致的典型原因我在幾次對比中總結(jié)出幾個最常見的不一致來源變壓器模型是否包含勵磁支路和飽和非線性。手算基本忽略勵磁電流Simulink的變壓器模型如果選擇飽和飽和特性短路電流會略有偏差但對短路電流影響不大如果變壓器模型沒有選擇模擬飽和則和手算接近。故障電阻。手算假設金屬性短路R_fault0仿真為了數(shù)值收斂填了0.01Ω多了一個微小的電阻測出來的電流會偏小一點可以忽略。負荷的影響。手算短路電流通常忽略負荷但Simulink模型里如果接了額定負載故障前的電流不是0故障瞬間的初始相位條件會變進而影響非周期分量但周期分量有效值基本不受負荷影響。線路的分布電容。手算用集中阻抗Simulink的分布式參數(shù)線路模型Distributed Parameters Line在短路瞬間會有暫態(tài)行波看起來波形高頻振蕩這和簡化手算差別很大。對于短線路建議用Series RLC Branch或PI型線路模型波形更干凈結(jié)果也更接近手算。4.3 用腳本批量對比參數(shù)我做參數(shù)校核時會用腳本循環(huán)修改線路長度把仿真后的周期分量有效值和手算值放在一起畫曲線這樣能直觀看出趨勢和偏差。做法是把模型參數(shù)用變量L_km替代在仿真腳本里用set_param改參數(shù)并運行sim再把結(jié)果存到數(shù)組。這里提供一段簡化腳本思路L_vec 2:2:20; I_sim zeros(size(L_vec)); for i 1:length(L_vec) set_param(short_circuit/Line RLC, R, num2str(0.17*L_vec(i))); set_param(short_circuit/Line RLC, L, num2str(0.35*L_vec(i)/(2*pi*50))); sim(short_circuit); % 提取故障后電流穩(wěn)態(tài)峰值 I_sim(i) max(abs(out.iabc(300:end,1))) / sqrt(2); end這里set_param只改阻抗模型會自動按新參數(shù)重新仿真非常方便掃參。注意L的單位是亨利需要把線路電抗除以314。5. 實操中踩過的坑與調(diào)試經(jīng)驗5.1 仿真步長與求解器設置三相短路仿真最怕步長太大把沖擊電流峰值給抹平了。默認的變步長求解器在快速暫態(tài)過程中會自動減步長但也可能不夠密。我通常在模型配置參數(shù)里把最大步長固定為1e-4s0.1ms求解器選ode23tb或ode15s這兩種對電力電子和電磁暫態(tài)問題比較友好。如果只關(guān)心波形宏觀趨勢可以放開到0.5ms但讀沖擊電流時必須用0.1ms左右否則讀出的峰值可能偏小5%以上。另一個坑是仿真總時長。故障時間設0.1s仿真總時長至少要0.2s才能看到故障后的暫態(tài)衰減過程我一般設0.4s如果還想看故障清除后的恢復過程就設0.6s并在故障模塊里配置第二個Transition time。5.2 測量環(huán)節(jié)的讀數(shù)陷阱用Three-Phase V-I Measurement輸出到工作區(qū)后數(shù)據(jù)是結(jié)構(gòu)體不是數(shù)組。我一開始直接plot出錯了后來用out.iabc.signals.values提取才行。如果測量模塊里選擇的是Use a label那么會被命名成變量需要在模型里用Signal Label把標簽對應好。Scope里看到的電流是瞬時值最大值是峰值不是有效值。很多新手直接把Scope里的峰值當作短路電流有效值去套公式結(jié)果對不上。要取有效值應該把故障后達到穩(wěn)態(tài)的波形做RMS計算或者用峰值除以sqrt(2)。還要注意區(qū)分是故障相電流還是線電流三相故障時三相電流等大但如果是單相故障讀故障相和非故障相意義完全不同。另外故障模塊如果設置了外部控制Step模塊從0跳到1之后要確保Step的Final value是1若Final value是0故障模塊可能一直在導通或斷開狀態(tài)檢查仿真波形有無異常通常一眼能看出來。5.3 模型擴展從三相短路到其他故障類型三相短路分析是基礎(chǔ)但實際系統(tǒng)中單相接地、兩相短路也很常見。Simulink里切換故障類型很快Three-Phase Fault模塊選擇Phase A fault、Phase B fault或AB fault即可。如果是接地故障記得把Fault resistance和接地方案配好要觀察零序電流和三相不對稱波形。故障位置也能隨便挪我把故障模塊從線路末端移到母線側(cè)、變壓器出口、發(fā)電機出口測不同位置的短路電流能直觀理解短路點越靠近電源阻抗越小短路電流越大。這種擴展對一個校核項目來說很有價值比如斷路器能不能裝在當前母線位置就看該點的短路電流是否在遮斷容量以內(nèi)。5.4 關(guān)于模型收斂的最后提醒有時候模型運行報錯顯示“Solver encountered an error”十有八九是初始條件或參數(shù)突變造成的。常見解決辦法把故障電阻從0.01Ω提高到0.1Ω看能否運行把Powergui的Simulation type切到PhasorPhasor模式不看暫態(tài)波形但能快速算穩(wěn)態(tài)短路電流把負荷從模型中斷開用純阻抗網(wǎng)絡測試。我在調(diào)試時習慣把故障模塊的Transition time設置成0.1s但先不接入線路直接測空載模型是否穩(wěn)定再加故障一層層排查能省很多時間。總的來說基于MATLAB做三相短路分析本質(zhì)上不是學會拖模塊而是能把系統(tǒng)阻抗、暫態(tài)過程、沖擊電流這些物理概念用數(shù)字曲線表達出來并反復驗證。這個流程跑通之后后面不管換什么系統(tǒng)接線我都是先用腳本估一個數(shù)再搭Simulink模型精細看波形最后回到手算公式解釋偏差。整個過程里最有價值的一點就是不要盲信仿真波形任何一組數(shù)據(jù)都要能回到物理原理上講得通才算真可靠。本文還有配套的精品資源點擊獲取 MATLAB三相短路分析實戰(zhàn)：從暫態(tài)原理到Simulink仿真建模
簡介基于MATLAB的三相短路分析資料包適合電力系統(tǒng)專業(yè)學生與工程技術(shù)人員用于掌握故障建模、仿真與結(jié)果分析的方法。資源從三相短路基本概念出發(fā)結(jié)合Simulink平臺完整演示了電力系統(tǒng)模型搭建、故障類型與位置設置、仿真運行及數(shù)據(jù)解讀的全過程并給出可復用的test.m腳本和詳細課設報告。壓縮包共12個文件大小僅1.09MB包含m腳本、PDF報告、4個txt數(shù)據(jù)文件與6張PNG結(jié)果圖txt文件保存了電路原始數(shù)據(jù)、節(jié)點阻抗/導納矩陣及輸出結(jié)果PNG圖清晰呈現(xiàn)發(fā)電機等效前后的矩陣變化便于對照報告逐項理解。目前已有2576人學習尤其適合正在完成電力系統(tǒng)課設或MATLAB仿真作業(yè)的讀者。通過研讀代碼與報告可以快速掌握用Simulink分析三相短路電流、驗證保護設備動作特性并能直接修改參數(shù)擴展至其他故障類型是一份實用的電力仿真學習資料。 做短路計算校驗的時候三相短路幾乎是每個電力系統(tǒng)從業(yè)者都躲不開的基本功。我最早接觸三相短路分析是在學暫態(tài)過程的時候當時全靠手算公式、查沖擊系數(shù)曲線算一條輻射線路的短路電流要折騰大半天后來真正參與設備選型和保護整定發(fā)現(xiàn)短路電流校核幾乎貫穿所有項目才開始系統(tǒng)性用MATLAB來搭仿真模型。MATLAB里的Simscape Electrical模塊老版本叫SimPowerSystems支持直接拖拽電源、變壓器、線路和故障模塊設置一個故障時間就能看到三相短路后的電壓電流波形比手算直觀太多。這篇把整套流程和踩過的坑整理出來給準備用MATLAB做三相短路分析的同學和工程師一個可以照做的參考。1. 三相短路分析的本質(zhì)算的是暫態(tài)電流不是穩(wěn)態(tài)電流1.1 為什么會短路短路電流為什么可怕三相短路是指供電系統(tǒng)中三相之間發(fā)生金屬性或非金屬性的直接連接相當于負載被短路掉回路阻抗驟減系統(tǒng)電壓被瞬間拉低短路電流可以高達額定電流的數(shù)倍甚至數(shù)十倍。為什么都拿三相短路作為最嚴重的故障工況來校核因為對稱短路時正序阻抗最小短路電流通常最大設備動穩(wěn)定、熱穩(wěn)定、斷路器遮斷容量、保護定值整定幾乎都要圍繞這個最大工況來做。真正讓三相短路分析變得棘手的是暫態(tài)過程。短路瞬間電流不能突變?yōu)榱司S持磁鏈守恒回路里會感應出一個非周期分量電流這個直流分量和周期分量疊加導致短路電流波形嚴重偏置。用公式來表達故障相的短路電流大約為i(t) sqrt(2) I * [sin(ωt α - φ) - sin(α - φ) * e^(-t/Ta)]其中I是次暫態(tài)電流有效值α是故障瞬間的電壓相位角φ是短路回路阻抗角Ta L/R為直流分量衰減時間常數(shù)。初始的沖擊電流峰值在半個周波附近出現(xiàn)50Hz系統(tǒng)約0.01s其幅值可以近似為 Kimp * sqrt(2) * IKimp就是沖擊系數(shù)取值通常在1.8左右。如果不理解這個暫態(tài)過程仿真跑出來的波形就只是“看到了一個波形”沒辦法從里面提取設備校核需要的關(guān)鍵參數(shù)。1.2 MATLAB做短路分析的兩條路線用MATLAB做三相短路分析我試過兩條路線各有適用場景。第一條是純腳本編程把系統(tǒng)等效成R-L阻抗用解析公式直接算短路電流周期分量、非周期分量和沖擊電流。好處是快速、透明適合簡單輻射網(wǎng)和批量參數(shù)掃描比如我要比較不同線路長度對短路電流的影響寫一個循環(huán)就出來了。缺點是處理復雜拓撲和多節(jié)點網(wǎng)絡時手寫節(jié)點阻抗矩陣非常費勁很容易出錯。第二條是用Simulink里的Simscape Electrical模塊搭建數(shù)學模型這也是我后期項目里用得最多的一條路。它把電源、變壓器、線路、負載、故障模塊做成可視化元件搭完模型就能看到各點電壓電流波形還能方便地切換故障類型、故障時刻和故障位置。缺點是模型自由度多求解器、步長、阻抗參數(shù)設置不當結(jié)果可能歪到離譜。我個人的建議是做教學演示或驗證理論公式優(yōu)先用腳本做工程系統(tǒng)分析、保護整定或需要給同事看波形時用Simulink建模。兩者結(jié)合是最好的先用腳本算一個基準值再拿Simulink去驗證和擴展。2. 建模前的參數(shù)整理從單線圖到Simulink元件參數(shù)2.1 系統(tǒng)阻抗與標幺值換算不管用什么工具三相短路分析的輸入基礎(chǔ)都是系統(tǒng)阻抗。模型里的每個元件都要有明確的R和X這些參數(shù)可以從設備銘牌、線路單位阻抗表和系統(tǒng)短路容量中獲取。Simulink里可以直接填歐姆Ω和亨H但我建議在做理論校核時先換算成標幺值這樣對比手算結(jié)果更方便。基準值的選取遵循 S_B 100 MVA或10 MVA基準電壓取電網(wǎng)額定電壓例如10kV系統(tǒng)選U_B 10.5kV基準阻抗為 Z_B U_B^2 / S_B。舉個例子我常用的一個驗證模型是10kV輻射網(wǎng)系統(tǒng)電源側(cè)短路容量為200MVA折算到10.5kV基準下的電源等效阻抗約 Z_source U_B^2 / S_k 10.5^2 / 200 0.551Ω。變壓器容量2MVA短路電壓百分比 U_k% 6%折算到10.5kV側(cè)變壓器阻抗約為 Z_T U_k% * U_B^2 / (100 * S_T) 6 * 10.5^2 / (100 * 2) 3.307Ω。線路用10km的10kV架空線每公里電抗0.35Ω電阻0.17Ω線路阻抗 Z_L 1.7 j3.5Ω。這些值最后填到Simulink元件里就是有名值。2.2 在Simscape Electrical庫中選型與搭建新版MATLAB的Simulink庫瀏覽器里路徑一般是Simscape Electrical Specialized Power Systems Fundamental Blocks。老版本直接叫SimPowerSystems位置類似。需要用到的基本元件有Three-Phase Source三相電壓源設置相電壓/線電壓有效值、頻率、內(nèi)阻和電感。Three-Phase Transformer (Two Windings)雙繞組三相變壓器設置額定容量、頻率、一次二次額定電壓、短路阻抗百分比。Three-Phase Series RLC Branch線路阻抗可以直接填R、L。Three-Phase Fault三相故障模塊可選故障類型和切換時間。Three-Phase V-I Measurement測量三相電壓電流輸出到Scope或Workspace。Powergui電力系統(tǒng)圖形用戶界面模塊選擇連續(xù)模式或相量模式。我一般選Continuous才能看到暫態(tài)波形。搭建順序不復雜三相電壓源接到變壓器一次側(cè)變壓器二次側(cè)經(jīng)過線路阻抗到故障模塊再接到負載。故障模塊放在短路點測量模塊放在電源側(cè)和故障點兩側(cè)各一個方便對比電壓跌落和電流變化。別忘了在模型中放置Powergui模塊否則整個模型會報錯。把上節(jié)計算出的參數(shù)填入元件對話框Three-Phase Source的Phase-to-phase RMS voltage填10.5kV如果電源側(cè)視為無窮大系統(tǒng)內(nèi)阻為0但為了仿真接近實際建議填0.5Ω和0.01H變壓器一次10.5kV、二次0.4kV也可以不降壓用于等值短路計算線路RLC填1.7Ω和0.35H對應10km注意電抗換算到電感需要除以2πf。2.3 故障模塊參數(shù)設置細節(jié)Three-Phase Fault模塊是所有仿真里的核心開關(guān)。雙擊模塊Fault type選Three-phase faultABC這時三相同時短路。Fault resistance填0.01Ω不要填0數(shù)值上可以更穩(wěn)定如果做金屬性短路填一個很小的值即可。切換方式我習慣用外部控制External control用Step模塊從0跳到1觸發(fā)故障模塊導通也可以用內(nèi)部定時直接設置Transition time為0.1s故障初始狀態(tài)選Open經(jīng)過0.1s后切換到Closed。這里有個容易被忽略的點Transition time是故障開始時間不是故障持續(xù)時間。如果希望故障持續(xù)0.2s后恢復需要設置兩個Transition time比如0.1和0.3s狀態(tài)依次為Open→Closed→Open。只填一個時間故障會一直保持到仿真結(jié)束。對于短路分析我們通常只關(guān)心故障后前幾個周波所以故障持續(xù)0.2s足夠。我還習慣在故障模塊的External Fault Breaker控制端口接一個Signal Generator或者Step這樣可以精確控制故障觸發(fā)相位。想研究短路初始角的影響就通過Step的延遲時間調(diào)整短路發(fā)生的合閘角。3. 仿真運行與波形解讀短路的真正信息藏在暫態(tài)過程里3.1 從波形中讀出周期分量和非周期分量模型搭好設置仿真時長0.4s最大步長0.1ms運行后雙擊Scope看到故障點電流波形。穩(wěn)定情況下看到故障前是規(guī)則正弦故障發(fā)生后電流跳變前半段波形上下嚴重不對稱包絡線一高一低隨后逐漸恢復到對稱正弦。這個不對稱就是非周期分量在衰減而對稱部分就是短路電流的周期分量。我在實際讀波形時會把數(shù)據(jù)導到MATLAB工作區(qū)再畫圖方便標注峰值。比如用To Workspace模塊輸出三相電流然后在命令窗口提取數(shù)據(jù)t out.tout; iA out.iabc(:,1); plot(t, iA); xline(0.1, --); grid on;如果仿真的故障開始時間是0.1s波形圖里0.1s之前的正弦就是負荷電流或空載電流0.1s之后才是短路電流。用數(shù)據(jù)游標讀出故障后第一個周期波峰大約15kA再往后幾個周波峰值回落到約8kA前者包含了直流分量后者才接近周期分量幅值。3.2 沖擊電流與衰減時間常數(shù)的實測方法沖擊電流是短路后第一個峰值是設備動穩(wěn)定校驗的核心參數(shù)。我在模型中讀取A相電流最大值為15kA左右與理論估算的沖擊電流對照。理論值怎么算先由短路回路總阻抗算出次暫態(tài)周期分量有效值I然后乘以sqrt(2)再乘以沖擊系數(shù)Kimp。Kimp可以用1.8的典型值也可以按Kimp 1 e^(-0.01/Ta)計算。T_a L/R (總電抗/2πf)/總電阻。如果L/R算下來是0.03sKimp 1 e^(-0.333) ≈ 1.717理論峰值約為 sqrt(2) * I * 1.717。衰減時間常數(shù)還能直接從波形里測出來把故障后的電流包絡線提取出來直流分量衰減到初始值的36.8%所需的時間就是Ta。由于非周期分量在示波器里表現(xiàn)為包絡線的不對稱差可以取上包絡線減下包絡線再除以2得到直流衰減曲線用對數(shù)坐標擬合得出時間常數(shù)。這個方法適合驗證實際模型中的回路時間常數(shù)是否與手算一致。3.3 不同故障時刻對波形的影響短路初始角對沖擊電流影響很大。這是因為公式里sin(α-φ)決定了非周期分量的初始幅值。如果短路瞬間電壓過零α≈0回路阻抗角φ約80度sin(α-φ)很大直流分量幾乎以最大初始值出現(xiàn)沖擊電流最大如果短路瞬間電壓剛好在峰值附近非周期分量初始值很小波形基本對稱沖擊系數(shù)接近1。用MATLAB模擬這個現(xiàn)象很簡單在Three-Phase Source模塊里設置Phase angle of phase A或者用外部Step控制故障觸發(fā)時間把故障時間從0.1s改成0.105s等觀察第一個峰值的不同。我實測發(fā)現(xiàn)同一個模型只改故障時刻沖擊電流能從14kA變到16kA差別超過了10%。所以在做保護整定時不能光看一組仿真結(jié)果要找最嚴重的合閘角。4. 理論手算與仿真結(jié)果交叉驗證差距從哪來4.1 手算短路電流的簡化模型對于輻射網(wǎng)手算三相短路電流的簡化步驟是忽略負荷把電源等效為次暫態(tài)電勢E短路點電壓在故障前按額定電壓考慮短路電流有效值I_sc c U_N / (sqrt(3) Z_total)。電壓系數(shù)c在高壓系統(tǒng)通常取1.1低壓系統(tǒng)取1.05。Z_total是短路點到電源的全部等值阻抗之和包括電源內(nèi)阻抗、變壓器阻抗和線路阻抗。我用前面10kV模型來算Z_source 0.551ΩZ_T 3.307ΩZ_L 1.7 j3.5Ω。要注意電抗和電阻要分別加再求模。總電阻R 0.5 0 1.7 2.2Ω電源電阻很小忽略或按前面填的0.5Ω總電抗X 0.551 3.307 3.5 7.358Ω假設電源內(nèi)阻電抗0.551。Z_total sqrt(2.2^2 7.358^2) ≈ 7.68Ω。取c 1.1短路電流有效值 I_sc 1.1 * 10.5k / (sqrt(3) * 7.68) ≈ 0.868kA 868A。這里因為模型用了10kV/0.4kV變壓器和0.4kV側(cè)故障阻抗折算會不同手算時應注意基準側(cè)選擇為了簡單這里是在10kV側(cè)等值短路點。把仿真結(jié)果的短路電流周期分量有效值從波形里提煉出來兩者通常能對上偏差在5%以內(nèi)。如果偏差超過10%不要急著說仿真不對先去檢查參數(shù)折算和基準值。4.2 仿真和手算不一致的典型原因我在幾次對比中總結(jié)出幾個最常見的不一致來源變壓器模型是否包含勵磁支路和飽和非線性。手算基本忽略勵磁電流Simulink的變壓器模型如果選擇飽和飽和特性短路電流會略有偏差但對短路電流影響不大如果變壓器模型沒有選擇模擬飽和則和手算接近。故障電阻。手算假設金屬性短路R_fault0仿真為了數(shù)值收斂填了0.01Ω多了一個微小的電阻測出來的電流會偏小一點可以忽略。負荷的影響。手算短路電流通常忽略負荷但Simulink模型里如果接了額定負載故障前的電流不是0故障瞬間的初始相位條件會變進而影響非周期分量但周期分量有效值基本不受負荷影響。線路的分布電容。手算用集中阻抗Simulink的分布式參數(shù)線路模型Distributed Parameters Line在短路瞬間會有暫態(tài)行波看起來波形高頻振蕩這和簡化手算差別很大。對于短線路建議用Series RLC Branch或PI型線路模型波形更干凈結(jié)果也更接近手算。4.3 用腳本批量對比參數(shù)我做參數(shù)校核時會用腳本循環(huán)修改線路長度把仿真后的周期分量有效值和手算值放在一起畫曲線這樣能直觀看出趨勢和偏差。做法是把模型參數(shù)用變量L_km替代在仿真腳本里用set_param改參數(shù)并運行sim再把結(jié)果存到數(shù)組。這里提供一段簡化腳本思路L_vec 2:2:20; I_sim zeros(size(L_vec)); for i 1:length(L_vec) set_param(short_circuit/Line RLC, R, num2str(0.17*L_vec(i))); set_param(short_circuit/Line RLC, L, num2str(0.35*L_vec(i)/(2*pi*50))); sim(short_circuit); % 提取故障后電流穩(wěn)態(tài)峰值 I_sim(i) max(abs(out.iabc(300:end,1))) / sqrt(2); end這里set_param只改阻抗模型會自動按新參數(shù)重新仿真非常方便掃參。注意L的單位是亨利需要把線路電抗除以314。5. 實操中踩過的坑與調(diào)試經(jīng)驗5.1 仿真步長與求解器設置三相短路仿真最怕步長太大把沖擊電流峰值給抹平了。默認的變步長求解器在快速暫態(tài)過程中會自動減步長但也可能不夠密。我通常在模型配置參數(shù)里把最大步長固定為1e-4s0.1ms求解器選ode23tb或ode15s這兩種對電力電子和電磁暫態(tài)問題比較友好。如果只關(guān)心波形宏觀趨勢可以放開到0.5ms但讀沖擊電流時必須用0.1ms左右否則讀出的峰值可能偏小5%以上。另一個坑是仿真總時長。故障時間設0.1s仿真總時長至少要0.2s才能看到故障后的暫態(tài)衰減過程我一般設0.4s如果還想看故障清除后的恢復過程就設0.6s并在故障模塊里配置第二個Transition time。5.2 測量環(huán)節(jié)的讀數(shù)陷阱用Three-Phase V-I Measurement輸出到工作區(qū)后數(shù)據(jù)是結(jié)構(gòu)體不是數(shù)組。我一開始直接plot出錯了后來用out.iabc.signals.values提取才行。如果測量模塊里選擇的是Use a label那么會被命名成變量需要在模型里用Signal Label把標簽對應好。Scope里看到的電流是瞬時值最大值是峰值不是有效值。很多新手直接把Scope里的峰值當作短路電流有效值去套公式結(jié)果對不上。要取有效值應該把故障后達到穩(wěn)態(tài)的波形做RMS計算或者用峰值除以sqrt(2)。還要注意區(qū)分是故障相電流還是線電流三相故障時三相電流等大但如果是單相故障讀故障相和非故障相意義完全不同。另外故障模塊如果設置了外部控制Step模塊從0跳到1之后要確保Step的Final value是1若Final value是0故障模塊可能一直在導通或斷開狀態(tài)檢查仿真波形有無異常通常一眼能看出來。5.3 模型擴展從三相短路到其他故障類型三相短路分析是基礎(chǔ)但實際系統(tǒng)中單相接地、兩相短路也很常見。Simulink里切換故障類型很快Three-Phase Fault模塊選擇Phase A fault、Phase B fault或AB fault即可。如果是接地故障記得把Fault resistance和接地方案配好要觀察零序電流和三相不對稱波形。故障位置也能隨便挪我把故障模塊從線路末端移到母線側(cè)、變壓器出口、發(fā)電機出口測不同位置的短路電流能直觀理解短路點越靠近電源阻抗越小短路電流越大。這種擴展對一個校核項目來說很有價值比如斷路器能不能裝在當前母線位置就看該點的短路電流是否在遮斷容量以內(nèi)。5.4 關(guān)于模型收斂的最后提醒有時候模型運行報錯顯示“Solver encountered an error”十有八九是初始條件或參數(shù)突變造成的。常見解決辦法把故障電阻從0.01Ω提高到0.1Ω看能否運行把Powergui的Simulation type切到PhasorPhasor模式不看暫態(tài)波形但能快速算穩(wěn)態(tài)短路電流把負荷從模型中斷開用純阻抗網(wǎng)絡測試。我在調(diào)試時習慣把故障模塊的Transition time設置成0.1s但先不接入線路直接測空載模型是否穩(wěn)定再加故障一層層排查能省很多時間。總的來說基于MATLAB做三相短路分析本質(zhì)上不是學會拖模塊而是能把系統(tǒng)阻抗、暫態(tài)過程、沖擊電流這些物理概念用數(shù)字曲線表達出來并反復驗證。這個流程跑通之后后面不管換什么系統(tǒng)接線我都是先用腳本估一個數(shù)再搭Simulink模型精細看波形最后回到手算公式解釋偏差。整個過程里最有價值的一點就是不要盲信仿真波形任何一組數(shù)據(jù)都要能回到物理原理上講得通才算真可靠。本文還有配套的精品資源點擊獲取