微電網(wǎng)穩(wěn)定性這個領(lǐng)域我入行前五年其實一直處在“會用工具但講不清原理”的狀態(tài)。仿真軟件里點幾下能算出特征值、能畫出功角曲線但別人一問“這個特征值實部為什么是負的”“等面積法則到底在算哪塊面積”我就答不上來了。直到后來被一個做微電網(wǎng)示范工程的調(diào)試現(xiàn)場逼著啃了小半年的數(shù)學(xué)基礎(chǔ)才算是把穩(wěn)定性分析這根鏈條打通。這篇文章就把我梳理出來的小信號分析與暫態(tài)穩(wěn)定兩條主線做個系統(tǒng)整理適合剛接觸微電網(wǎng)控制的研究生、工程師以及想從“看仿真結(jié)果”進階到“理解仿真過程”的從業(yè)者。1. 為什么微電網(wǎng)穩(wěn)定性不能只會“看結(jié)論”1.1 微電網(wǎng)穩(wěn)定性到底在討論什么微電網(wǎng)和我們熟悉的傳統(tǒng)電力系統(tǒng)有一個核心區(qū)別傳統(tǒng)電網(wǎng)依賴大型同步發(fā)電機的轉(zhuǎn)動慣量來維持頻率和電壓而微電網(wǎng)里絕大多數(shù)電源是逆變器接口的分布式電源比如光伏、儲能、風電。逆變器沒有物理上的旋轉(zhuǎn)部件它輸出的電壓和頻率完全由控制算法決定。也就是說微電網(wǎng)的穩(wěn)定性和控制器參數(shù)、控制結(jié)構(gòu)、通信延遲甚至濾波器的阻抗特性都直接掛鉤。這就帶來了一個很實際的問題你用仿真軟件搭一個微電網(wǎng)模型改一個下垂系數(shù)可能系統(tǒng)就振蕩了再改一個低通濾波器時間常數(shù)振蕩又消失了。平臺上看這個加一個低通濾波器看起來只是一個小改動背后實際影響的是整個閉環(huán)系統(tǒng)的動態(tài)特性。如果不懂小信號分析的基本原理就只能靠試錯來調(diào)參數(shù)試驗效率很低而且一旦系統(tǒng)規(guī)模變大、電源數(shù)量變多試錯的空間幾乎不存在。我常跟身邊的人說微電網(wǎng)穩(wěn)定性研究核心要回答三個問題第一給定一個工作點系統(tǒng)在小擾動下能不能回到原來的平衡狀態(tài)第二如果發(fā)生大擾動比如負荷突變、逆變器切機系統(tǒng)能不能過渡到新的穩(wěn)定狀態(tài)第三從數(shù)學(xué)上能不能定量給出穩(wěn)定的邊界以及哪些參數(shù)在推動邊界移動。前兩個問題分別對應(yīng)小信號穩(wěn)定和暫態(tài)穩(wěn)定第三個問題才是“數(shù)學(xué)基礎(chǔ)”四個字的真正含義。1.2 兩條主線小尺度與大尺度穩(wěn)定性分析按照擾動規(guī)模和時間尺度分成兩條主線這個劃分不是學(xué)術(shù)上拿著玩兒的而是直接決定你用哪套數(shù)學(xué)工具。小信號分析研究的是系統(tǒng)在某個穩(wěn)態(tài)工作點附近受到微小擾動后的動態(tài)行為。所謂“微小”指的是擾動幅度小到可以讓非線性模型在平衡點附近做線性化處理。這種分析的時間尺度通常是毫秒到秒級別關(guān)心的是振蕩模態(tài)的阻尼、頻率以及哪些狀態(tài)變量參與了這個振蕩。小信號分析的數(shù)學(xué)核心是狀態(tài)空間模型、雅可比矩陣和特征值分析本質(zhì)上使用線性代數(shù)和常微分方程穩(wěn)定性理論也就是李雅普諾夫第一法的線性化形式。暫態(tài)穩(wěn)定研究的是系統(tǒng)經(jīng)歷大擾動后的行為例如三相短路、大負荷投切這類讓系統(tǒng)結(jié)構(gòu)或運行點發(fā)生劇烈改變的事件。此時線性化假設(shè)不再成立因為系統(tǒng)的狀態(tài)可能被推到了很遠的位置非線性項主導(dǎo)了動態(tài)過程。暫態(tài)穩(wěn)定分析的時間尺度一般是秒級甚至幾十秒級關(guān)心的問題是功角或者逆變器輸出相位是否還能保持同步電壓和頻率是否還能恢復(fù)到可接受的范圍。數(shù)學(xué)核心是非線性微分方程、能量函數(shù)和數(shù)值積分方法。這兩條線不是互相獨立的。實際工程中我遇到過不少情況小信號分析顯示系統(tǒng)有很好的穩(wěn)定裕度但一做大擾動仿真系統(tǒng)卻在膝跳式地振蕩幾下后失穩(wěn)了。原因在于大擾動激發(fā)了非線性區(qū)域里的動力學(xué)行為小信號分析根本覆蓋不到那個區(qū)域。所以成熟的分析流程一定是小信號分析和暫態(tài)仿真配合使用前者幫你理解模態(tài)、調(diào)參數(shù)后者幫你驗證大擾動下的生存能力。2. 小信號分析的數(shù)學(xué)骨架從狀態(tài)空間到特征值2.1 建狀態(tài)空間模型的常規(guī)流程小信號分析的第一步是建立微電網(wǎng)的狀態(tài)空間模型。說完這句話很多新手會問微電網(wǎng)模型里有逆變器、濾波器、線路、負荷還有控制器狀態(tài)變量怎么選選狀態(tài)變量的原則是凡是含記憶元件或積分環(huán)節(jié)的地方都要有狀態(tài)變量。對逆變器來說LC濾波器中的電感電流和電容電壓是自然的狀態(tài)變量控制環(huán)節(jié)里的積分器輸出、低通濾波器的輸出、鎖相環(huán)內(nèi)部的積分狀態(tài)也都要納入狀態(tài)變量。舉個具體例子一臺電壓控制型逆變器比如采用下垂控制的逆變器加上LC濾波器加上本地負荷狀態(tài)變量至少包括兩個電感電流、兩個電容電壓、下垂控制器中兩個功率濾波器的輸出、電壓控制環(huán)中兩個積分器的輸出、電流控制環(huán)中兩個積分器的輸出。這還沒算鎖相環(huán)的狀態(tài)。一臺逆變器輕松到十來個狀態(tài)變量一個含十臺逆變器的微電網(wǎng)就有上百個狀態(tài)變量所以微電網(wǎng)小信號建模是典型的“多狀態(tài)變量”問題。建模型的常用方法有兩大類一種是對每臺逆變器單獨建立dq坐標系下的狀態(tài)空間方程然后通過坐標變換統(tǒng)一到公共坐標系一種是直接建立整個系統(tǒng)的微分代數(shù)方程組然后在平衡點線性化。實際中大部分文獻和工具采用前一種因為分模塊建模思路清晰便于檢查模型錯誤。狀態(tài)空間方程的標準形式一般寫為dx/dt A x B u y C x D u其中x是狀態(tài)向量u是輸入向量一般是參考值或擾動輸入A矩陣就是系統(tǒng)動態(tài)特性的核心載體。小信號穩(wěn)定性只看A矩陣的特征值B、C、D主要用于研究輸入對狀態(tài)的影響和輸出觀測不直接影響穩(wěn)定性結(jié)論。2.2 平衡點與小擾動線性化的關(guān)鍵細節(jié)很多人把“在平衡點線性化”這句話輕輕帶過了但這其實是整個小信號分析最容易被工程化帶偏的一步。所謂平衡點就是讓所有狀態(tài)變量導(dǎo)數(shù)為零的工作點。找平衡點看起來只需要解一組代數(shù)方程實際上麻煩在于微電網(wǎng)的潮流方程和控制器參考值之間互相耦合。一個典型的場景你用下垂控制的逆變器給一個恒功率負荷供電。下垂控制的目標是讓輸出頻率和電壓隨有功無功變化所以平衡點的電壓和頻率不是預(yù)先給定的而是整個系統(tǒng)的穩(wěn)態(tài)條件決定的。你需要先做潮流計算得到每個節(jié)點的電壓幅值和相角再反推逆變器輸出功率進而確定控制器穩(wěn)態(tài)輸出值最后才能把全部狀態(tài)變量在平衡點的數(shù)值算出來。線性化的過程則是把一個形如dx/dt f(x, u)的非線性系統(tǒng)在平衡點附近展開。設(shè)平衡點為x0、u0令x x0 Δx、u u0 Δu于是f(x, u) ≈ f(x0, u0) ?f/?x |(x0,u0) · Δx ?f/?u |(x0,u0) · Δu因為平衡點上有f(x0, u0) 0所以得到dx/dt A Δx B ΔuA矩陣就是雅可比矩陣。這個展開丟掉的是Δx的二階及以上項所以適用前提是Δx足夠小。這里有一個工程上非常容易踩的坑dq坐標系下坐標變換本身是角度相關(guān)的非線性環(huán)節(jié)。如果你把坐標變換方程也做了線性化那么A矩陣里必須包含相角的偏差量。很多初學(xué)模型的人會忽略相角和鎖相環(huán)動態(tài)在擾動后的作用導(dǎo)致特征值分析結(jié)果偏樂觀。我的經(jīng)驗是鎖相環(huán)動態(tài)和線路電阻在微電網(wǎng)小信號模型中的影響往往被低估它們恰恰是阻尼不足振蕩的常見來源。2.3 特征值、阻尼與參與因子怎么讀拿到A矩陣之后核心工作就是對A做特征值分解。A的特征值λ是方程det(A - λI) 0的根。如果狀態(tài)變量是實數(shù)特征值要么是實數(shù)要么是以共軛對形式出現(xiàn)的復(fù)數(shù)。對于連續(xù)的線性時不變系統(tǒng)穩(wěn)定性判據(jù)很簡單全部特征值實部都小于零系統(tǒng)漸進穩(wěn)定只要存在一個實部大于零的特征值系統(tǒng)不穩(wěn)定實部為零則處于臨界穩(wěn)定工程上視為不安全。為什么實部正負有這種決定性的意義可以從解的形態(tài)來理解。對線性系統(tǒng)來說每個特征值對應(yīng)的運動模式隨時間演化是e^(λt)。如果λ σ jω那么e^(λt) e^(σt) · e^(jωt)前一個是幅度增減因子后一個是振蕩因子。σ小于零時這個模式隨時間衰減最終消失系統(tǒng)回到平衡點σ大于零時這個模式不斷增長振鈴越來越大一直到非線性約束生效系統(tǒng)失穩(wěn)。σ為零時振蕩不衰減系統(tǒng)會持續(xù)擺動這種臨界狀態(tài)在工程上也是不可接受的。特征值是共軛復(fù)數(shù)時振蕩頻率就是ω/(2π)阻尼比ζ定義為ζ -σ / sqrt(σ2 ω2)阻尼比是很有工程意義的量。微電網(wǎng)中小信號振蕩的阻尼比如果小于0.05你會看到功率波形衰減得很慢要振蕩好幾個周期才能穩(wěn)定下來。如果阻尼比為負那就是增幅振蕩對應(yīng)特征值實部為正。你直覺上覺得逆變器的響應(yīng)好像很“拖泥帶水”其實大概率就是某個振蕩模態(tài)阻尼不足。特征值本身只能告訴你系統(tǒng)有哪些振蕩模式還不能讓你知道哪個模式是哪臺逆變器惹的禍。這時要用參與因子矩陣。參與因子是一個把特征值和狀態(tài)變量關(guān)聯(lián)起來的指標數(shù)值上是左特征向量和右特征向量對應(yīng)元素的乘積歸一化結(jié)果。實際閱讀方法是對于某個關(guān)心的振蕩模態(tài)某對共軛特征值看參與因子表里哪些狀態(tài)變量數(shù)值最大這些狀態(tài)變量就是該模態(tài)的主要參與者。我調(diào)試一個含四臺逆變器的微電網(wǎng)仿真時發(fā)現(xiàn)有一個頻率大約11 Hz的振蕩模態(tài)阻尼比只有0.02。參與因子分析結(jié)果顯示這個模態(tài)主要參與的是兩臺并聯(lián)儲能逆變器的有功下垂濾波環(huán)節(jié)和LC濾波器電流。最后定位到問題是兩臺逆變器的下垂系數(shù)差異過大導(dǎo)致有功分配不均加上功率濾波時間常數(shù)設(shè)置太小有功環(huán)路由產(chǎn)生了負阻尼。花了一個上午調(diào)整下垂系數(shù)和濾波時間常數(shù)把阻尼比拉到了0.15。沒有參與因子分析我只能瞎猜各臺機器的貢獻。2.4 為什么根軌跡和參數(shù)靈敏度很重要特征值分析給的是“當前一組參數(shù)下系統(tǒng)穩(wěn)不穩(wěn)定”的快照。工程上更關(guān)心的往往是“參數(shù)在什么范圍內(nèi)變化時系統(tǒng)還能穩(wěn)定”這就是根軌跡分析的意義。根軌跡的做法是選定一個關(guān)心的參數(shù)比如下垂系數(shù)m_p在合理范圍內(nèi)連續(xù)變化畫出系統(tǒng)特征值在復(fù)平面上的移動軌跡。如果某條軌跡隨著參數(shù)增加從左半平面穿過虛軸進入右半平面那這個參數(shù)的臨界值就是系統(tǒng)穩(wěn)定邊界。靈敏度分析則是定量看雅可比矩陣元素變化對特征值的影響。數(shù)學(xué)上特征值對某個參數(shù)的靈敏度可以通過左右特征向量計算?λ_i / ?k (w_i^T · ?A/?k · v_i) / (w_i^T · v_i)其中w_i和v_i分別是對應(yīng)特征值的左、右特征向量。這個公式不復(fù)雜但意義很大它直接告訴你最有控制效益的調(diào)節(jié)方向。我在分析一個微電網(wǎng)方案時發(fā)現(xiàn)特征值靈敏度最高的是逆變器輸出濾波器電容值而不是控制器增益。這意味著如果設(shè)計階段能優(yōu)化電容選型遠比費勁調(diào)PI參數(shù)效率高。這種判斷沒有靈敏度分析做支撐就只能是玄學(xué)。根軌跡和靈敏度分析本質(zhì)上都是在回答“穩(wěn)定裕度”的問題系統(tǒng)離失穩(wěn)有多遠。行業(yè)里常用“阻尼比不小于0.1”作為微電網(wǎng)設(shè)計的一個經(jīng)驗準則而完整的小信號分析流程能幫你量化這個裕度到底夠不夠而不是等真機跑起來振蕩了才后知后覺。3. 暫態(tài)穩(wěn)定的數(shù)學(xué)底色非線性、能量與臨界切除3.1 暫態(tài)過程為什么要“換一套”數(shù)學(xué)工具大擾動發(fā)生后微電網(wǎng)各個換流器的輸出相位和功率會出現(xiàn)大幅擺動。此時系統(tǒng)的數(shù)學(xué)模型依然是那組微分方程但你不能在某個平衡點附近做線性化因為狀態(tài)變量的偏離幅度太大了。假設(shè)偏離量小這件事一旦不成立線性化模型給出的運動軌跡就完全不是系統(tǒng)真實的運動軌跡。數(shù)學(xué)上暫態(tài)穩(wěn)定分析面對的是一個非線性自治或非自治系統(tǒng)形如dx/dt f(x)現(xiàn)在首要任務(wù)是判斷系統(tǒng)在大的初始擾動下是否收斂到某個穩(wěn)定平衡點以及會經(jīng)過怎樣的大幅擺動過程。這里有一個和線性系統(tǒng)完全不同的關(guān)鍵現(xiàn)實非線性系統(tǒng)可能存在多個平衡點有些穩(wěn)定、有些不穩(wěn)定初始條件不同收斂到的平衡點也可能不同。這意味著對大擾動穩(wěn)定性來說初始狀態(tài)在哪里、擾動要去哪這兩個信息必須同時掌握。傳統(tǒng)電力系統(tǒng)的暫態(tài)穩(wěn)定分析一直用“單機無窮大母線”模型來起步。對于微電網(wǎng)來說這種簡化更復(fù)雜一些因為微電網(wǎng)里沒有明確的無窮大母線多臺逆變器之間相互影響。但分析思想是一樣的通過等面積法則EEAC方法、或者從能量函數(shù)角度看暫態(tài)過程中的加速面積和減速面積來評估穩(wěn)定性。3.2 等面積法則的推導(dǎo)與使用邊界等面積法則是理解暫態(tài)穩(wěn)定的第一道門檻。它從搖擺方程出發(fā)本質(zhì)上是在說一個同步發(fā)電機對于微電網(wǎng)我們可以將一臺逆變器或聚合的分布式電源看作一個等效的擺動體的轉(zhuǎn)子運動方程近似為M · d2δ/dt2 P_m - P_e其中M是慣性時間常數(shù)δ是功角P_m是機械對于逆變器來說是輸入?yún)⒖脊β使β蔖_e是電磁輸出功率。這個方程描述的是輸入功率和輸出功率不平衡時功角隨時間的變化。暫態(tài)過程可以分解為幾個階段。以最經(jīng)典的輸電線路三相短路為例短路前系統(tǒng)運行在某個功角δ0輸出功率為P_e0短路期間線路無法傳輸正常功率P_e大幅下降輸入功率不變所以出現(xiàn)了加速功率P_m - P_e。功角開始增大。如果繼電保護在某個時間切除故障系統(tǒng)拓撲改變輸出功率特性曲線也變了此時如果功角已經(jīng)增大太多減速面積小于加速面積系統(tǒng)就會失去同步。等面積法則就是在說暫態(tài)過程中功角從初始值增大到切除時刻值這一段系統(tǒng)吸收了多少加速能量對應(yīng)加速面積切除后系統(tǒng)能吸收多少減速能量對應(yīng)減速面積。如果減速面積大于或等于加速面積那么功角在達到最大值后開始回落系統(tǒng)能恢復(fù)穩(wěn)定如果加速面積更大功角將持續(xù)增大甚至越過不穩(wěn)定平衡點系統(tǒng)失穩(wěn)。數(shù)學(xué)上加速面積和減速面積分別是對應(yīng)時間段內(nèi)不平衡功率對功角的積分。雖然這里討論起來看著像古典電機學(xué)內(nèi)容但它在微電網(wǎng)暫態(tài)分析中依然有直接對應(yīng)的應(yīng)用比如逆變器在低電壓穿越期間的功率傳輸特性變化、故障期間有功輸出的跌落、故障切除后恢復(fù)控制參數(shù)的動作都可以用類似“加速面積-減速面積”的框架來判斷臨界切除時間的量級。當然微電網(wǎng)的暫態(tài)分析比單機無窮大復(fù)雜得多但核心思想“擾動后的能量積累與耗散能力對比”是非常通用的。值得強調(diào)的是等面積法則的前提是系統(tǒng)可以化簡為一個等效搖擺方程。微電網(wǎng)多機系統(tǒng)的暫態(tài)穩(wěn)定問題需要擴展的等面積法則EEAC也就是把多機系統(tǒng)在某個時段內(nèi)聚合成一個單機等效系統(tǒng)來做近似判斷。EEAC在某些場景下很準在某些場景下又會偏樂觀或偏保守。工程上我會把它當快速篩選工具用最終結(jié)論還是靠全系統(tǒng)的時域仿真來驗證。等面積法則的價值不是精確判穩(wěn)而是給你一個直觀的物理圖像和臨界切除時間的初始估計。3.3 李雅普諾夫直接法判斷穩(wěn)定但不給軌跡除了數(shù)值積分和等面積法則暫態(tài)穩(wěn)定還有一條理論路線李雅普諾夫直接法。這個方法的思想是不直接求解微分方程而是構(gòu)造一個能量函數(shù)V(x)然后看這個函數(shù)沿系統(tǒng)軌跡的導(dǎo)數(shù)變化。對于微電網(wǎng)模型如果存在一個正定函數(shù)V(x)沿著系統(tǒng)軌跡滿足dV/dt ≤ 0那么系統(tǒng)在平衡點附近是穩(wěn)定的。如果進一步d2V/dt2也滿足一些條件還能給出漸近穩(wěn)定結(jié)論。這個思路在物理上很像“水往低處流”找到一個能量度量系統(tǒng)狀態(tài)只會沿著能量減小的方向運動最后停在能量最低點。李雅普諾夫直接法在微電網(wǎng)暫態(tài)穩(wěn)定分析中的應(yīng)用最有代表性的方向是基于結(jié)構(gòu)保持模型的能量函數(shù)。這類方法需要構(gòu)造包含線路儲能、電容器儲能和控制器狀態(tài)的李氏函數(shù)。實際應(yīng)用中最大的問題在于微電網(wǎng)含有大量逆變器控制環(huán)這些控制環(huán)不是簡單的保守系統(tǒng)純能量函數(shù)經(jīng)常無法滿足嚴格的微分約束。所以工程上直接法主要用于給出保守的穩(wěn)定域估計或者配合時域仿真判斷穩(wěn)定裕度的相對大小。我在實際項目里的做法是先用李雅普諾夫能量函數(shù)估算一個偏保守的臨界切除時間下限再用時域仿真掃出準確的臨界切除時間。兩者對比之后如果仿真值比理論值大很多說明系統(tǒng)還有優(yōu)化空間如果兩者接近說明系統(tǒng)已經(jīng)工作在比較極限的狀態(tài)這時候設(shè)計保護定值要留足裕度。4. 實操記錄一個微電網(wǎng)穩(wěn)定性分析項目的完整流程4.1 案例背景與系統(tǒng)參數(shù)這一節(jié)我拿一個自己做過的小型微電網(wǎng)仿真項目來串一遍完整流程。硬件配置是一個380 V交流母線上的微電網(wǎng)包含兩臺儲能逆變器、一臺光伏逆變器、一組恒阻抗負荷和一組恒功率負荷。儲能逆變器采用電壓源型的下垂控制光伏逆變器采用電流源型的MPPT控制。分析方法上小信號部分用狀態(tài)空間建模加特征值分析暫態(tài)部分用電磁暫態(tài)仿真。系統(tǒng)的基本參數(shù)如下額定容量每臺逆變器30 kWLC濾波器電感0.5 mH、電容50 μF開關(guān)頻率10 kHz下垂控制中有功下垂系數(shù)2%即滿載時頻率偏差1 Hz無功下垂系數(shù)5%即滿載時電壓偏差5%功率濾波時間常數(shù)0.1 s。線路阻抗從幾十毫歐到幾百毫歐不等。選這個規(guī)模是有講究的。臺數(shù)不多你可以手工檢查每一步模型但狀態(tài)變量已經(jīng)超過60個也足夠展示特征值分析、參與因子和根軌跡的實際操作流程。如果你一上來就做十臺機的模型光排查模型錯誤就夠你熬夜的。4.2 小信號分析的實施步驟與判讀具體實施分為以下幾步第一步是建立每臺逆變器的dq坐標系下狀態(tài)空間模型。這一步的工作量最大因為逆變器的五個部分都有各自的狀態(tài)方程。以儲能逆變器為例下垂控制部分有兩個狀態(tài)有功和無功濾波值電壓環(huán)有兩個積分狀態(tài)電流環(huán)有兩個積分狀態(tài)LC濾波器有四個狀態(tài)兩個電感電流、兩個電容電壓。把這些方程按“控制環(huán)給指令、主電路執(zhí)行指令”的邏輯串起來。第二步是把各臺逆變器、線路和負荷方程組裝成一個整體系統(tǒng)。這個組裝過程最容易出錯的就是坐標變換。每臺逆變器的dq坐標系是相對于自己鎖相得到的角度定義的必須把所有狀態(tài)同步到公共參考坐標系。我習(xí)慣先把整個系統(tǒng)的方程寫成矩陣形式再用符號計算工具做集總求雅可比人工推導(dǎo)容易錯且很難檢查。第三步是計算平衡點。首先用潮流計算求節(jié)點電壓幅值和相角然后根據(jù)功率平衡條件求每臺逆變器的輸出功率再反推下垂控制器的穩(wěn)態(tài)輸出頻率和電壓幅值進而求所有控制環(huán)的穩(wěn)態(tài)狀態(tài)。這一串計算下來你會得到一組完整的平衡點狀態(tài)向量。第四步才輪到雅可比矩陣和特征值計算。得到特征值后先分群看哪些是實特征值對應(yīng)非振蕩衰減模式哪些是共軛特征值對應(yīng)振蕩模式。然后對振蕩模態(tài)逐個算阻尼比重點關(guān)注阻尼比低于0.05的模態(tài)。最后用參與因子分析定位每個模態(tài)的物理來源。我在這個案例里發(fā)現(xiàn)一個有意思的現(xiàn)象兩臺儲能逆變器之間有一個振蕩模態(tài)頻率大概8 Hz阻尼比只有0.03。參與因子顯示主導(dǎo)這個模態(tài)的是一號儲能逆變器的有功下垂濾波狀態(tài)和二號儲能逆變器的電容電壓狀態(tài)兩個設(shè)備的LC濾波器耦合形成了功率振蕩。繼續(xù)看靈敏度分析發(fā)現(xiàn)改變兩臺逆變器之間的線路阻抗對阻尼比影響最大而改變下垂系數(shù)的影響次之。后來我把線路感抗從0.3 mH增加到0.5 mH阻尼比從0.03升到了0.09。這個結(jié)果初看起來反直覺因為線路電抗增加通常會讓無功分配能力變差但在這個場景下它提高了兩臺機器的動態(tài)解耦程度反而改善了阻尼。這就是靈敏度分析比拍腦袋好用的地方。4.3 暫態(tài)仿真與臨界切除時間求解小信號分析做完后我習(xí)慣緊接著做一輪大擾動測試檢驗小信號結(jié)論的魯棒性。做法是在電磁暫態(tài)仿真平臺上搭建同樣的系統(tǒng)設(shè)置三相短路故障發(fā)生在公共耦合點故障持續(xù)時間從50 ms開始逐步增加每次增加20 ms觀察系統(tǒng)是否還能恢復(fù)同步。第一步是觀察擾動后的有功功率、頻率和電壓波形故障切除后功率波形如果出現(xiàn)振蕩收斂系統(tǒng)是穩(wěn)定的如果振蕩不收斂而是幅度越來越大就說明到了穩(wěn)定邊界。用二分法或者逐步掃描法逼近臨界切除時間。在我的案例中臨界切除時間大約在160 ms附近。第二步是回到小信號結(jié)果去解釋大擾動行為的合理性。這個案例中故障切除后頻率能夠恢復(fù)但出現(xiàn)了一段持續(xù)約1秒的功率振蕩。這個振蕩頻率和小信號分析里8 Hz模態(tài)對應(yīng)不上說明不是線性模態(tài)主導(dǎo)而是大擾動把系統(tǒng)推入了非線性區(qū)域。這時再用等面積法則估算臨界切除時間作為交叉驗證。等面積法則給出的估計值是120 ms左右比仿真值保守一些處于安全側(cè)。通過這個案例你應(yīng)該能看到小信號分析和暫態(tài)穩(wěn)定并不是完全獨立的兩套分析它們互相印證、互相補充。小信號分析沒有直接給出“耐大擾動能力”的結(jié)論但提供了模態(tài)信息和穩(wěn)定裕度參考暫態(tài)分析雖然直接回答“能不能活下來”但定位問題根源仍要靠小信號工具。4.4 分析工具與腳本組織經(jīng)驗我在微電網(wǎng)穩(wěn)定性分析中用到的工具整體圍繞“模型清晰、結(jié)果可追溯”來組織。建模過程我用符號計算確認方程系統(tǒng)數(shù)值計算采用腳本編程最終仿真驗證用電磁暫態(tài)工具。實際工作中更重要的是腳本的組織方式。我的習(xí)慣是一個案例一個目錄目錄下分成模型定義腳本每臺逆變器一個模塊、系統(tǒng)組裝腳本線路和負荷連接關(guān)系、平衡點計算腳本、特征值分析和參與因子輸出腳本、參數(shù)掃描腳本根軌跡和靈敏度、結(jié)果輸出目錄特征值表、參與因子表、波形數(shù)據(jù)。全部基于文本的腳本格式組織好處是版本管理方便改了參數(shù)可以立刻看到哪一步結(jié)果變了也方便給同事解釋問題出在哪。還有一個細節(jié)求解特征值之前一定要先做歸一化檢查。狀態(tài)變量里有的量綱是安培、有的是伏特、有的是標幺值功率如果直接對原始單位的A矩陣求特征值數(shù)值尺度差異可能引發(fā)病態(tài)問題。一般建議把所有狀態(tài)變量做標幺化處理再做雅可比計算。這個習(xí)慣能減少很多數(shù)值噪聲導(dǎo)致的假特征值。5. 最容易踩坑的幾個地方5.1 小信號穩(wěn)定不等于暫態(tài)穩(wěn)定這是我在實際項目中見過最多的認知誤區(qū)。小信號穩(wěn)定只保證系統(tǒng)在小擾動下能回到平衡點但對大擾動沒有任何承諾。一個在特征值分析中阻尼比很高的小信號穩(wěn)定系統(tǒng)在大擾動下完全可能因超出控制環(huán)的線性工作范圍而失穩(wěn)。控制環(huán)的限幅、PWM的飽和、鎖相環(huán)在大頻率偏移下的非線性行為都會在小信號模型中被忽略或簡化而在大擾動中恰恰是這些約束能決定生死。所以正確的工程流程一定是“小信號分析暫態(tài)仿真”雙軌并行。小信號分析負責參數(shù)設(shè)計和模態(tài)分析暫態(tài)仿真負責安全驗證。任何拿小信號穩(wěn)定裕度直接推斷“系統(tǒng)大擾動也穩(wěn)定”的做法都值得警惕。5.2 內(nèi)部節(jié)點消除的誤差為了降低系統(tǒng)階數(shù)很多人會在建模過程中把一些代數(shù)變量消掉比如消除純代數(shù)節(jié)點電壓或者消去非重點關(guān)注的線路內(nèi)部節(jié)點。這種降階往往省事但不一定安全。以輻射狀微電網(wǎng)為例饋線末端的節(jié)點電壓受負荷波動影響很大。如果在建模時把這些節(jié)點當作代數(shù)變量直接消去相當于默認這些節(jié)點電壓在動態(tài)過程中滿足瞬時代數(shù)關(guān)系這對準穩(wěn)態(tài)分析問題不大但對暫態(tài)過程來說可能引入較大誤差。尤其是帶有動態(tài)負荷或恒功率負荷的節(jié)點電壓動態(tài)和負荷動態(tài)是耦合的消去節(jié)點相當于丟掉了耦合信息。我的建議是分兩步走第一階段可以用降階模型做快速的參數(shù)掃描和根軌跡分析等鎖定了候選參數(shù)組后切換到全階模型做復(fù)核。兩條腿走路既保證效率又避免降階誤差影響最終結(jié)論。5.3 判斷矩陣正定的三件事小信號穩(wěn)定分析中使用李雅普諾夫函數(shù)來判穩(wěn)可以構(gòu)造某些二次型能量函數(shù)并要求對應(yīng)的矩陣正定。這里有幾個工程中常被忽略的前提第一矩陣正定的結(jié)論依賴于你選的坐標基底。線性變換不改變特征值結(jié)論但會改變矩陣正定性判別的表達形式所以任何判穩(wěn)結(jié)論都要注明坐標系。第二系統(tǒng)矩陣A的特征值判據(jù)是充分必要條件而李雅普諾夫方程P的求解結(jié)果則依賴你選的Q矩陣。實際工程里用李雅普諾夫方程時結(jié)果只能說明你找到了某個能量函數(shù)是下降的不代表所有可能的能量函數(shù)都是下降的。也就是說有時候給了“充分但非必要”的結(jié)論。第三求解李雅普諾夫方程前要檢查A是否滿足可控性和可觀測性條件。如果系統(tǒng)存在不可觀的模態(tài)該模態(tài)的特征值不影響外部輸出你會誤以為它對系統(tǒng)穩(wěn)定性無關(guān)緊要實際上它仍然參與內(nèi)部動態(tài)。這三件事看著是理論潔癖但在工程項目評審或者調(diào)試排故時能幫你避免很多無謂的爭論。5.4 模型維度與計算成本的取舍微電網(wǎng)狀態(tài)空間模型的維度膨脹速度遠超多數(shù)人的直覺。一臺含完整控制器模型的逆變器動輒十幾個狀態(tài)變量一個含二十臺逆變器加線路動態(tài)的系統(tǒng)就可能逼近上千階。高維度帶來的問題不僅僅是計算慢更是數(shù)值算法穩(wěn)定性。對上千階的矩陣求特征值時特征值算法的數(shù)值誤差可能產(chǎn)生虛假的弱阻尼模態(tài)。我覺得有三條實用經(jīng)驗可以分享第一在建模階段就對“哪些狀態(tài)必須保留”做取舍。對暫態(tài)穩(wěn)定分析鎖相環(huán)動態(tài)必須保留對小信號分析功率濾波狀態(tài)和LC濾波器動態(tài)必須保留但如果只是為了觀察機電振蕩對本微電網(wǎng)來說即逆變器輸出功率的低頻擺動逆變器的電流內(nèi)環(huán)因為帶寬高可以適當簡化。第二預(yù)先做量綱歸一化處理避免數(shù)值病態(tài)。第三特征值計算完成后用逆冪迭代法驗證阻尼比較小的模態(tài)確認真實性。6. 最后說點實操上的體會做微電網(wǎng)穩(wěn)定性分析這幾年我最大的體會是越是覺得“仿真軟件已經(jīng)幫我算好了”的地方越值得回去翻數(shù)學(xué)推導(dǎo)。特征值是誰的特征值、線性化時丟掉了什么、暫態(tài)過程的能量從哪來又去哪了這些問題的答案全在數(shù)學(xué)里但表達出來的卻是一個個真實的物理過程。如果非要給一條最快的上手路徑我的建議是挑一個三到五臺逆變器的小系統(tǒng)手工推導(dǎo)完整的dq模型寫腳本算出特征值然后用電磁暫態(tài)仿真去驗證每個關(guān)鍵模態(tài)在大擾動下的表現(xiàn)。這個過程看起來費時間但能一次性把模型、算法、物理意義全部串起來之后再看任何微電網(wǎng)穩(wěn)定性分析報告你都能夠快速抓住哪些是湊數(shù)的結(jié)論、哪些才是真正決定系統(tǒng)安全的判斷。最后再分享一個小技巧在你的分析腳本里把每一次計算的特征值表、參與因子表和暫態(tài)響應(yīng)的截圖統(tǒng)一歸檔成帶備注的文本記錄。穩(wěn)定性的判斷不是只看一次計算結(jié)果而是看參數(shù)變化時系統(tǒng)行為的趨勢。手邊如果沒有歷史對比數(shù)據(jù)你很難說清楚這個穩(wěn)定裕度到底是提高了還是下降了。系統(tǒng)的記錄習(xí)慣在這個領(lǐng)域的重要性不亞于任何數(shù)學(xué)推導(dǎo)。