汶川地震汶川地震2008年，中國汶川發生8.0級地震，管網破壞嚴重，直接經濟損失高達8451億元。2012年颶風桑迪在美國新澤西州大西洋城附近沿海登陸，造成18個州超過820萬用戶和商家停電，直接經濟損失約為740億美元。研究背景電/氣/熱綜合能源系統→多種能源綜合能源系統在實現多能互補互濟的同時，也可能造成事故在不同能源系統之間耦合傳播，擴大事故影響，甚至導致多種能源供應同時停止，嚴重威臺北市新回雪風雪風臺灣“815”大停電●未收波及臺灣“815”大停電風電廠韌性指標的計算包括兩部分參數：一、極端災害情況下出現某種元件故障情況的概率；二、極端災害導致的系統負荷削減的大小。Qit)Rys=E(S,)Qit)Rys=E(S,)·故障場景下源、儲可用度臺風場景1地震場景1冰災場景1臺風場景245一多類型極端災害場景概率建模一多類型極端災害場景概率建模故障場景下電氣負荷削減改進算法故障場景下電氣負荷削減改進算法基于復用影響增量庫的韌性評估方法基于復用影響增量庫的韌性評估方法行進方向，初始中心氣壓差和行進速度。臺風概率模型應用組合枚,其發生概率Pw行進方向θ服從雙正態分布6臺風通常難以破壞室內設施及地下元件，故本研究只考慮架空線路在臺風下的受損情況。臺風攻擊模型根據Batts風場理論求得直接參數，即各地實時風速vt),并基于微元法計算架空線臺風路徑00yaya)①v①海岸線臺風風場示意圖臺風基礎參數風眼實時坐標7各地實時風速v(t)7《中國地震動參數區劃圖》提出了三級潛在震源區劃分方案：在地震基礎參數為震級和震中點。地震概率模型首先枚舉震級檔，再枚舉該震級檔對應的潛在震中點，從而生成潛在地震概率集。MAMMAMucAMMucAM地震活動斷層構造源A構造源B背景源C背景源D構造源A構造源B背景源C背景源DM~MD■Mg~M~MD震級上限■Mac~MB■Mac~MB81.2典型極端災害建模——地震概率模型參考烈度橢圓衰減關系，地震攻擊模型將等烈度線視為參考烈度橢圓衰減關系，地震攻擊模型將等烈度線視為同心同向、層層嵌套的橢圓群，并由震級和震中點推算橢圓方向、半徑和中心，進而繪制地震烈度I的分布圖，計算架空線路、變壓器和燃氣管道基烈度6區外邊界44烈度6區外邊界44長軸烈度8區烈度7區烈度6區烈度8區烈度8區烈度I分布地震烈度分布示意圖地震烈度分布示意圖9二三2.1電氣綜合能源系統的解耦優化框架考慮燃氣電廠等電氣耦合裝置，構建基于解耦優化的電氣互聯系開始確定燃氣電廠有功出力Pa根據耗量方程計算燃氣電廠供氣負荷Gu氣負荷削減優化確定燃氣電廠的供氣負荷削減量G否修改燃氣電廠有功出力上限否修改燃氣電廠有功出力上限Pumu優化收斂u計算總負荷削減量結束針對氣網最優負荷求解難收斂，受初值影響大等問題，將簡化混合整數線性規劃模型為第一階段模型，第一階段實數優化解作為第二階段非線性連續模型的內點法初值，整數優化解則在第二階段模型中沿用，從而實現對簡化解的修正。第一階段混合整數線性模型優化結果實數型變量氣源出力節點氣壓壓縮比氣負荷削減量初值壓縮機狀態管道燃氣流向沿用第二階段非線性連續模型氣負荷削減優化非線性連續模型第二階段固定了燃氣流向和壓縮機狀態，并以第一階段優化解為目標函數氣負荷削減量最小min約束條件流量-氣壓Weymouth方程f2=K,|n2-π}|f2=K,(g,-g,)(π2-π2)氣壓縮機耗量方程第一階段(忽略)壓縮機約束兩算法的計算結果存在明顯偏差，說明氣網第二階段優化模型發兩算法的計算結果存在明顯偏差，說明氣網第二階段優化模型發所提兩階段模型在三階以內故障狀態下均實現收斂，這說明氣網優化第一階段能夠為第二階段模型提供一個足夠好的內點法初值。參數單階段兩階段Q?-數關系，可將復雜的多能流優化問題轉化為簡單的函數計算，減少故障狀態分析所需的優化計算次數。型轉化為方程組，并構建了數據驅動的相似性匹配技術與相似矩陣快速求逆最優解最優解基于狀態相似性的線性方程組求解可靠性評估狀態相似性特征集合基于數據驅動的相似性特征選取特征提取狀態相似性特征集合基于數據驅動的相似性特征選取特征提取特征1特征N特征1特征N一多類型極端災害場景概率建模一多類型極端災害場景概率建模故障場景下電氣負荷削減改進算法故障場景下電氣負荷削減改進算法基于復用影響增量庫的韌性評估方法基于復用影響增量庫的韌性評估方法3.1基于影響增量的韌性評估方法針對多類型設備故障導致的系統狀態組合爆炸問題，將針對多類型設備故障導致的系統狀態組合爆炸問題，將高階故障狀態的影響拆分折算至相應的低階狀態，實現了單次狀態分析覆蓋多個系統狀態，提升評估效率。設備故障逐個分析可能的故障狀態，效率低、精度差源荷波動單次分析覆蓋多個故障狀態，效率大大提升3.1基于影響增量的韌性評估方法利用影響增量△I,替換掉Is,是傳統狀態枚舉法的改進方法，通過將狀態可獲得較高精度，提升評估效率。傳統狀態枚舉法(SE)引入影響增量基于影響增量的狀態枚舉法(IISE)pAI+p?AI?+'P?P?△I12低階狀態高階狀態(保留)(忽略)不隨災害場景的變化并反復調取，實現計算速度的大幅提升。 樣本抽樣抽樣抽樣dKdKrsys,Krm,KRIIDDSDW1面向規劃的韌性指標基于影響增量公式推導全局靈敏度，能夠分析不同參數之間的相互作用對韌性指標的影響，準確反映元件強化策略的韌性提升效果，指導韌性的精準提升。靈敏度分析各元件失效概率模型結果模型結果模型結果模型結果模型結果將電氣互聯系統附在中國臺灣將電氣互聯系統附在中國臺灣臨海地震區上作為算例，并應用組合枚舉法生成了10000個臺風場景和13806個地震場景構成潛在災害場景集。以特定地震場景為例驗證所提韌性評估方法的可行性。場景震級屬于7~7.5震級檔，此災害場景同時影響電網和氣網的元件，有一定的代表性。構造源A構造源B烈度8區背景源C背景源D代表性100km烈度7區近海水域0烈度6區(a)電力子系統(b)天然氣子系統0000于SE法，而3階IISE法誤差僅為0.92%,且復用影響增量庫在確保精度的同時能夠顯著提高韌性評估速度。方法rsys(MW)誤差(%)計算時間(s)-46.785