基于故障樹-蒙特卡洛的輸電可靠性設(shè)計

劉 洋，王 佳，江 泳，高 源

（1.北京智芯半導(dǎo)體科技有限公司，北京 102202；2.中核控制系統(tǒng)工程有限公司，北京 102401）

隨著電網(wǎng)規(guī)模的不斷擴(kuò)大，架空輸電線路也隨之增多，然而，由于雨雪、結(jié)冰、雷擊等造成架空輸電線路的故障也在不斷考驗著電網(wǎng)的輸電能力[1-3]。作為電網(wǎng)輸電系統(tǒng)的命脈，為避免對電力系統(tǒng)和居民生活帶來重大損失，分析并提高線路的可靠性對電網(wǎng)的穩(wěn)定運行意義重大[4-5]。

架空輸電線路的可靠性是指線路在規(guī)定的時間和條件下輸送電能的能力，目前研究架空輸電線路的可靠性多采用定性分析或?qū)＜掖蚍址绞絒6-8]。該文依據(jù)故障樹分析方法建立架空輸電線路故障樹模型，通過故障樹-蒙特卡洛方法對系統(tǒng)執(zhí)行定量可靠性分析，明確系統(tǒng)各部件的重要度，從而確定系統(tǒng)薄弱環(huán)節(jié)[9]。改進(jìn)或加強這些環(huán)節(jié)的可靠性，對架空輸電線路的穩(wěn)定運行具有現(xiàn)實指導(dǎo)意義[10]。

1 架空輸電線路可靠性設(shè)計

根據(jù)架空輸電線路的結(jié)構(gòu)特點、運行數(shù)據(jù)和檢修記錄，采用故障樹分析方法建立架空輸電線路的故障樹模型。將架空輸電線路故障作為故障樹頂事件，研究各部件故障的概率、相互關(guān)系以及導(dǎo)致的結(jié)果。

架空輸電線路主要由圖1所示的8個部分組成[11]。

圖1 架空輸電線路組成

1.1 故障樹分析

故障樹本質(zhì)上是一種邏輯因果關(guān)系圖，它將系統(tǒng)的故障狀態(tài)設(shè)為頂層事件，然后找出導(dǎo)致這種故障可能的直接原因（中間事件），如此層層追溯和跟蹤，直至找出導(dǎo)致系統(tǒng)故障的所有基本原因（底事件）。利用邏輯門符號將圖形化模型轉(zhuǎn)化為數(shù)學(xué)模型，方便后續(xù)計算失效概率和系統(tǒng)中各部件的重要度[12]。

故障樹頂事件命名為架空輸電線路故障S，故障樹中間事件命名為Mj(j=1,2,…,14)，故障樹底事件命名為Xi(i=1,2,…,37)。將底事件進(jìn)行細(xì)分并根據(jù)其邏輯關(guān)系建立的故障樹如圖2 所示。

圖2 架空輸電線路故障樹

故障樹分析法是根據(jù)基本部件的故障概率來計算頂事件的故障率，同時根據(jù)系統(tǒng)結(jié)構(gòu)函數(shù)計算各基本部件的重要度，從而對系統(tǒng)的安全性及可靠性做出定量評價。然而，大型復(fù)雜系統(tǒng)搭建的故障樹規(guī)模及其龐大，理論計算費時費力，開展較為困難，急需結(jié)合其他方法來提高計算速度。

1.2 蒙特卡洛方法

蒙特卡洛方法又被稱為隨機抽樣統(tǒng)計試驗方法或概率模擬方法。對架空輸電線路而言，其真實的故障樹結(jié)構(gòu)通常較為復(fù)雜，往往會出現(xiàn)多級中間事件，僅僅采用數(shù)學(xué)公式計算難以得到頂事件的故障率及各基本部件的重要度。而蒙特卡洛法可有效解決此類問題，利用計算機來快速準(zhǔn)確計算故障樹頂事件的故障率，獲得各基本部件的重要度排序，確定系統(tǒng)的薄弱環(huán)節(jié)和頂事件的故障情況，為定量分析故障樹提供有力支持。

文中搭建了架空輸電線路的故障樹模型，在結(jié)構(gòu)函數(shù)的基礎(chǔ)上，通過蒙特卡洛仿真方法獲得系統(tǒng)的平均故障率及各基本事件的重要度，確定故障處理的優(yōu)先次序。可靠性分析流程如圖3所示。

圖3 可靠性分析流程

該文結(jié)合故障樹和蒙特卡洛法，根據(jù)故障樹分析法獲得的結(jié)構(gòu)函數(shù)，搭建合理的概率模型，通過執(zhí)行數(shù)值化試驗獲得抽樣值，通過處理試驗數(shù)據(jù)，獲得問題解。程序結(jié)構(gòu)簡便靈活、求解容易、占用計算機內(nèi)存少，仿真流程如圖4 所示。

圖4 仿真流程

1）初始化系統(tǒng)

分析故障樹得出系統(tǒng)基本事件個數(shù)m、狀態(tài)量函數(shù)、系統(tǒng)結(jié)構(gòu)函數(shù)。

2）初始化仿真

采用事件調(diào)度法驅(qū)動仿真策略，確定仿真次數(shù)N、仿真時間Tmax，將Tmax等分為M個時間間隔。

3）仿真建模

根據(jù)系統(tǒng)邏輯搭建故障樹概率模型。

4）數(shù)據(jù)采集

利用MATLAB 軟件生成仿真所需的隨機數(shù)，通過反函數(shù)得出相應(yīng)的時間變量[13]：（第j次仿真時，基本事件xi的失效時間）。

第j次仿真時，基本事件的隨機數(shù)序列為Tj=(t1j,t2j,…,tij,…,tmj),j=1,…,N。

仿真N次后，得到N行m列的隨機數(shù)矩陣：

5）排序矩陣

對基本事件的隨機數(shù)序列按從小到大進(jìn)行排序，得出第j次仿真時，基本事件的時間排序序列：

其中，TTFij≤TTF(i+1)j,j=1,…,Ns。

仿真N次后，得到基本事件時間排序矩陣：

6）失效時間序列

假定各基本事件按排序依次失效，若在第j次仿真時（1≤j≤N），仿真時鐘t=TTFij與初始隨機數(shù)tij相同，則代表在第j次仿真時，第i個基本事件失效。計算此時系統(tǒng)狀態(tài)函數(shù)值，若系統(tǒng)故障則系統(tǒng)失效時間為tsj=TTFij，經(jīng)過N次仿真，得出失效時間序列：ts=(ts1,ts2,…,tsj,…,tsN)

7）失效次數(shù)

在M個時間間隔內(nèi)各放置計數(shù)器ΔNsr，統(tǒng)計上一步得出的系統(tǒng)失效時間落入各時間間隔的個數(shù)，得出系統(tǒng)失效次數(shù)序列：ΔNs=(ΔNs1,…,ΔNsr,…,ΔNsM)

8）可靠性指標(biāo)

在tr時刻，系統(tǒng)的不可靠度為，可靠度為。

1.3 重要度分析

重要度分析可用于對系統(tǒng)設(shè)計進(jìn)行科學(xué)評估，找出系統(tǒng)薄弱環(huán)節(jié)，為系統(tǒng)設(shè)計提供改進(jìn)建議。實踐經(jīng)驗證明，系統(tǒng)中各部件并非同樣重要的，例如串聯(lián)式結(jié)構(gòu)中各部件通常比并聯(lián)式結(jié)構(gòu)中的部件更重要。假如系統(tǒng)由若干部件串聯(lián)而成，一旦某部件出現(xiàn)故障，系統(tǒng)就無法正常工作[14]。

假如為了增強系統(tǒng)可靠性，將所有基本部件按冗余配置，那么隨著備件增多系統(tǒng)可靠性也在提高，系統(tǒng)建設(shè)費用也在提高。因此，引入重要度分析可以在一定的總費用下，研究如何配置各部件的備件使系統(tǒng)的可靠性達(dá)到最大[15]。

通常，將基本事件對頂事件發(fā)生的貢獻(xiàn)稱為重要度，在定量可靠性分析時，可按照分析對象和要求的不同來計算各基本事件的重要度。按基本事件對頂事件發(fā)生的影響程度大小進(jìn)行排序后，可以很容易找到系統(tǒng)的薄弱環(huán)節(jié)，對改進(jìn)設(shè)計、制定針對性維護(hù)策略、檢修保養(yǎng)起著非常重要的作用。

2 仿真過程設(shè)計

依據(jù)第一節(jié)研究，建立仿真流程如下：

1）依據(jù)已知的結(jié)構(gòu)函數(shù)，建立仿真抽樣模型；

2）規(guī)定仿真次數(shù)、系統(tǒng)最大運行時間和時間間隔等參數(shù)；

3）在每次仿真時，根據(jù)架空輸電線路的各基本部件的故障規(guī)律[16]，隨機抽樣產(chǎn)生對應(yīng)的失效時間，代入系統(tǒng)結(jié)構(gòu)函數(shù)判斷架空輸電線路頂事件是否失效，若沒失效，則繼續(xù)執(zhí)行仿真；若失效，則記錄下頂事件失效時的仿真時間、各基本部件的失效時間及部件的序號；

4）經(jīng)過大量重復(fù)抽樣，獲取并記錄步驟3）中的數(shù)值；

5）統(tǒng)計各時間區(qū)間內(nèi)的系統(tǒng)失效次數(shù)，描點繪制不可靠度和可靠度曲線，計算系統(tǒng)平均無故障工作時間；

6）計算引起頂事件失效的各基本部件失效次數(shù)及系統(tǒng)總失效次數(shù)，得出各基本部件的結(jié)構(gòu)重要度。

由大量工程實踐經(jīng)驗反饋可知，大型復(fù)雜系統(tǒng)的故障規(guī)律大多服從指數(shù)分布。因此，對基于故障樹-蒙特卡洛方法的架空輸電線路可靠性設(shè)計進(jìn)行仿真，流程如圖5 所示。

圖5 可靠性設(shè)計仿真流程

假定蒙特卡洛仿真中架空輸電線路總運行時間Tmax=15 000 h，將其等分為M個時間間隔，令M=15 000，則每個時間間隔即最小單位時間為1 h。仿真總次數(shù)定為Ns，通過多次試算，根據(jù)仿真結(jié)果是否趨于穩(wěn)定得出Ns=10 000，N為仿真次數(shù)的序號，故N=1,2,…,Ns。

利用MATLAB 軟件產(chǎn)生[0,1]區(qū)間的隨機數(shù)，根據(jù)反函數(shù)產(chǎn)生與之對應(yīng)的隨機變量，每次仿真運行中，按照自然排序從小到大的順序?qū)ij進(jìn)行排序，得到TTFij。

每次仿真，若基本部件發(fā)生故障，代入到結(jié)構(gòu)函數(shù)Φ(t)中，判斷架空輸電線路頂事件是否發(fā)生。若該基本部件的失效時間為TTFkj，此時頂事件發(fā)生，則系統(tǒng)失效時間tsj=TTFkj。判斷此時tsj屬于哪個區(qū)間，計數(shù)器加1，此次仿真結(jié)束。

重復(fù)上述操作Ns=10 000 次，得出架空輸電線路的不可靠度和可靠度曲線，如圖6 所示。

圖6 系統(tǒng)不可靠度和可靠度曲線

通過仿真，得出架空輸電線路平均無故障工作時間MTBF=1 354 h，平均故障率λ=7.384 5×10-4，各基本事件的結(jié)構(gòu)重要度如圖7 所示。

3 結(jié)論

分析架空輸電線路故障樹結(jié)構(gòu)對搭建電力傳輸設(shè)備狀態(tài)感知庫起著非常重要的作用，可以為后期設(shè)備健康狀態(tài)評價、檢修運維決策提供數(shù)據(jù)支撐。通過故障樹分析并借助蒙特卡洛仿真可以取代繁瑣復(fù)雜的公式計算方法，快速找到系統(tǒng)薄弱點，由圖7得知，架空輸電線路結(jié)構(gòu)重要度排序由大到小的前3個基本事件為X30、X31、X32，說明通道環(huán)境是架空輸電線路中最重要的事件，在架空輸電線路維護(hù)和檢修中應(yīng)重點設(shè)計，從而優(yōu)化架空輸電線路巡檢維護(hù)設(shè)計方案。

圖7 基本事件結(jié)構(gòu)重要度柱狀圖