蔡柏松

中國(guó)石化天然氣川氣東送管道分公司

大型天然氣管網(wǎng)有如下特點(diǎn)：①規(guī)模龐大。該系統(tǒng)包含的子系統(tǒng)﹑設(shè)備﹑進(jìn)氣點(diǎn)多，占有的空間大（距離長(zhǎng)﹑區(qū)域?qū)挬p地形起伏劇烈），投資巨大，運(yùn)行壽命期長(zhǎng)。若運(yùn)行管理正確，經(jīng)濟(jì)效益顯著；若運(yùn)行管理決策失誤，造成的經(jīng)濟(jì)損失難于彌補(bǔ)。②結(jié)構(gòu)復(fù)雜。該系統(tǒng)中各子系統(tǒng)﹑設(shè)備﹑節(jié)點(diǎn)之間的相互關(guān)系復(fù)雜，是多個(gè)相互制約和相互影響子系統(tǒng)有機(jī)結(jié)合的統(tǒng)一水動(dòng)力場(chǎng)。要使系統(tǒng)效益全面發(fā)揮，必須對(duì)該系統(tǒng)運(yùn)行管理進(jìn)行整體考慮。③功能綜合。該系統(tǒng)的目標(biāo)多樣（技術(shù)﹑經(jīng)濟(jì)﹑安全等等），因而，該系統(tǒng)的功能也是多方面的（經(jīng)營(yíng)管理﹑質(zhì)量控制﹑社會(huì)效益等）。④信息量大。該系統(tǒng)是多變量﹑多輸入﹑多輸出﹑多目標(biāo)﹑多參數(shù)﹑多干擾的系統(tǒng)，是必須借助計(jì)算機(jī)技術(shù)才能駕馭的系統(tǒng)，其最本質(zhì)的東西是工程海量信息的處理和工程活動(dòng)的輔助決策［1］。特別是在調(diào)峰調(diào)度管理中，如何對(duì)大型天然氣管網(wǎng)多個(gè)調(diào)峰方案進(jìn)行快速﹑準(zhǔn)確分析﹑評(píng)價(jià)以優(yōu)化該系統(tǒng)的運(yùn)行，提高經(jīng)濟(jì)效益，是天然氣管網(wǎng)調(diào)度人員關(guān)心的重大技術(shù)問題［1－2］。因此，筆者在大系統(tǒng)控制論［3］和多目標(biāo)決策思想［4］的啟發(fā)下，綜合應(yīng)用天然氣管道輸送理論［5］、綜合評(píng)價(jià)理論［6－7］和決策支持系統(tǒng)（Decision Support System，簡(jiǎn)稱DSS）理論［8－9］，開發(fā)研究《天然氣管網(wǎng)調(diào)峰方案評(píng)價(jià)系統(tǒng)》（Gas Pipeline Network Peak Shaving Project Evaluation System，簡(jiǎn)稱GPPES），對(duì)于更新工程管理觀念、優(yōu)化調(diào)峰調(diào)度運(yùn)行方案、節(jié)省經(jīng)營(yíng)管理費(fèi)用、縮短調(diào)度決策周期、提高天然氣管網(wǎng)運(yùn)行質(zhì)量具有重要的意義。

GPPES應(yīng)實(shí)現(xiàn)以下功能：在天然氣管網(wǎng)管轄區(qū)內(nèi)地形起伏、有若干進(jìn)氣點(diǎn)及已知計(jì)劃進(jìn)氣量的情況下，動(dòng)態(tài)掌握管網(wǎng)供氣狀況和各門站用氣狀況、若干進(jìn)氣點(diǎn)的最高（最低）壓力、各管段的儲(chǔ)氣能力、各壓氣站機(jī)組組合方式、壓氣站的進(jìn)、出站壓力、功率消耗［9］等；技術(shù)經(jīng)濟(jì)指標(biāo)包括動(dòng)力費(fèi)、管輸費(fèi)用等。當(dāng)具有多個(gè)預(yù)選方案的工藝參數(shù)和經(jīng)濟(jì)指標(biāo)后，用系統(tǒng)工程的思想方法整體考慮對(duì)各預(yù)選調(diào)峰方案進(jìn)行優(yōu)選排序。

預(yù)選調(diào)峰方案的生成必須根據(jù)天然氣管道輸送理論中氣體運(yùn)動(dòng)方程、連續(xù)性方程、能量方程和氣體狀態(tài)方程，采用隱式差分法或顯式差分法將上述偏微分方程離散為初等代數(shù)方程組，再根據(jù)管道具體結(jié)構(gòu)及運(yùn)行條件，引入相應(yīng)的邊界條件和初始條件進(jìn)行聯(lián)解才能得到。由于大型天然氣管網(wǎng)結(jié)構(gòu)的復(fù)雜性和計(jì)算方法的通用性，20世紀(jì)70年代開始，國(guó)外相關(guān)軟件公司、高等院校、科研院所相繼開發(fā)了天然氣管道仿真軟件，經(jīng)過近40年的應(yīng)用和改進(jìn)，目前比較成熟和可靠的商業(yè)化仿真軟件是 SPS （Stoner Pipeline Simulator，斯通納管道模擬器），它是目前國(guó)際上較先進(jìn)的瞬態(tài)流體仿真應(yīng)用程序，可以模擬任何現(xiàn)有的或規(guī)劃設(shè)計(jì)中的管道，對(duì)正?；蚍钦l件下各種不同控制策略的結(jié)果做出預(yù)測(cè)，可以模擬不同對(duì)象如用戶的用氣周期性變化、氣源斷流、壓縮機(jī)故障、管線出現(xiàn)泄漏以及管網(wǎng)系統(tǒng)的最長(zhǎng)工作時(shí)間等。SPS已被眾多國(guó)際知名工程公司所采用，用于管道設(shè)計(jì)及分析，一些公司已將SPS視為管道分析的標(biāo)準(zhǔn)。我國(guó)大型油氣管道管理和設(shè)計(jì)企業(yè)已用昂貴外匯引進(jìn)了該軟件，因此，如何發(fā)揮該軟件的使用效率，也是油氣輸送企業(yè)需要深入研究的技術(shù)問題之一［10－11］。筆者研究的GPPES，其預(yù)選調(diào)峰方案集的生成由SPS仿真軟件實(shí)現(xiàn)。以下分別介紹GPPES的功能設(shè)計(jì)和結(jié)構(gòu)設(shè)計(jì)。

1 GPPES的設(shè)計(jì)

1.1 功能設(shè)計(jì)

大型天然氣管網(wǎng)調(diào)峰方案評(píng)價(jià)系統(tǒng)應(yīng)具有如下功能（圖1）。

圖1 大型天然氣管網(wǎng)調(diào)峰方案評(píng)價(jià)系統(tǒng)的功能框圖

1）數(shù)據(jù)管理。對(duì)天然氣物性、用氣不均勻系數(shù)、沿線地溫傳熱系數(shù)等批量原始數(shù)據(jù)進(jìn)行追加、輸入、編輯、修改、計(jì)算、傳遞、打印及預(yù)處理時(shí)，具有數(shù)據(jù)管理功能。

2）SPS模擬管網(wǎng)工藝參數(shù)的調(diào)入。其中包括壓力、流量、溫度和功率的模擬（當(dāng)有壓縮機(jī)站、地下儲(chǔ)氣庫(kù)運(yùn)行時(shí)）。

3）儲(chǔ)氣調(diào)峰能力工藝計(jì)算。該模塊的功能是計(jì)算SPS模擬過程中各管段的儲(chǔ)氣能力，以及各干支線末段的儲(chǔ)氣能力，并將計(jì)算結(jié)果傳遞到“儲(chǔ)氣調(diào)峰能力評(píng)價(jià)方案”模塊中。

4）儲(chǔ)氣調(diào)峰能力評(píng)價(jià)。該模塊的功能是對(duì)系統(tǒng)建立的儲(chǔ)氣調(diào)峰評(píng)價(jià)模型分別用4種方法進(jìn)行評(píng)價(jià)。儲(chǔ)氣調(diào)峰能力的評(píng)價(jià)指標(biāo)是：管道儲(chǔ)氣量、儲(chǔ)氣最高壓力、儲(chǔ)氣最低壓力、增壓消耗功率、供需失衡戶數(shù)、實(shí)際供氣量、管輸費(fèi)用。

5）工況圖像顯示。為使調(diào)峰工藝參數(shù)和壓氣站工作特性的變化趨勢(shì)更加直觀，可將相關(guān)參數(shù)傳到繪圖系統(tǒng)中去，實(shí)現(xiàn)數(shù)據(jù)的可視化。

6）中間與評(píng)價(jià)成果輸出。各階段的數(shù)據(jù)處理結(jié)果都有相應(yīng)的報(bào)表輸出，報(bào)表輸出形式設(shè)計(jì)為.txt數(shù)據(jù)文檔，由用戶將數(shù)據(jù)文檔轉(zhuǎn)換為Word數(shù)據(jù)文檔，或轉(zhuǎn)換為Excel數(shù)據(jù)文檔，并連接到自己的生產(chǎn)報(bào)告文檔中去，成為生產(chǎn)報(bào)告的組成部分。

7）評(píng)價(jià)系統(tǒng)幫助瀏覽。為初學(xué)者提供操作指南，內(nèi)容包括：系統(tǒng)安裝、系統(tǒng)啟動(dòng)、運(yùn)行步驟、出錯(cuò)處理的詳細(xì)說明。

1.2 結(jié)構(gòu)設(shè)計(jì)

為實(shí)現(xiàn)上述功能設(shè)計(jì)提出的要求、決策人員操作的簡(jiǎn)捷性及儲(chǔ)氣調(diào)峰能力評(píng)價(jià)的自動(dòng)化和半自動(dòng)化，本系統(tǒng)結(jié)構(gòu)按DSS的理論與方法進(jìn)行設(shè)計(jì)（圖2）。

圖2 GPPES的結(jié)構(gòu)設(shè)計(jì)圖

GPPES的結(jié)構(gòu)由4大部分組成：評(píng)價(jià)數(shù)據(jù)庫(kù)及其管理系統(tǒng)、評(píng)價(jià)方法庫(kù)及其管理系統(tǒng)、評(píng)價(jià)模型庫(kù)及其管理系統(tǒng)、人機(jī)界面系統(tǒng)，如圖2虛框所示。

1）數(shù)據(jù)庫(kù)及其管理系統(tǒng)。GPPES中的數(shù)據(jù)庫(kù)結(jié)構(gòu)與一般的數(shù)據(jù)庫(kù)結(jié)構(gòu)不同，它所存儲(chǔ)的數(shù)據(jù)是調(diào)峰方案所需的數(shù)據(jù)。調(diào)峰方案評(píng)價(jià)的各模塊應(yīng)能方便有效地調(diào)用各種數(shù)據(jù)。處理后的數(shù)據(jù)又能方便地送入相應(yīng)的數(shù)據(jù)庫(kù)中存儲(chǔ)起來。本系統(tǒng)的數(shù)據(jù)庫(kù)按邏輯結(jié)構(gòu)分為8大類：天然氣物性數(shù)據(jù)庫(kù)、線路縱斷面數(shù)據(jù)庫(kù)、管網(wǎng)結(jié)構(gòu)數(shù)據(jù)庫(kù)、壓氣機(jī)組數(shù)據(jù)庫(kù)、技術(shù)經(jīng)濟(jì)指標(biāo)數(shù)據(jù)庫(kù)、儲(chǔ)氣調(diào)峰能力數(shù)據(jù)庫(kù)、儲(chǔ)氣調(diào)峰方案集評(píng)價(jià)數(shù)據(jù)庫(kù)。數(shù)據(jù)庫(kù)管理系統(tǒng)負(fù)責(zé)管理和維護(hù)GPPES所需要的各類數(shù)據(jù)，實(shí)現(xiàn)與方法庫(kù)、模型庫(kù)及人機(jī)界面的連接，使GPPES結(jié)構(gòu)中的各個(gè)模塊能方便有效地調(diào)用各種數(shù)據(jù)，完成對(duì)各類數(shù)據(jù)的處理和分析，達(dá)到有效決策的目的。

2）模型庫(kù)及其管理系統(tǒng)。模型庫(kù)是用來存儲(chǔ)輔助決策所需要的各類模型。本模型庫(kù)中有儲(chǔ)氣調(diào)峰能力評(píng)價(jià)模型、輸氣站影響力評(píng)價(jià)模型、管網(wǎng)動(dòng)態(tài)儲(chǔ)氣模型、輸氣站儲(chǔ)氣調(diào)峰模型［12］、曲線擬合模型［13］、表格模型［9］等。模型在 GPPES中采用2種存放形式：①把1個(gè)模型作為1個(gè)子程序，這樣模型庫(kù)實(shí)際上就成為程序庫(kù)；②有些模型也可以看作一組用數(shù)據(jù)集表示的關(guān)系，如壓氣機(jī)的特性參數(shù)［14］。這樣就可以作為數(shù)據(jù)存放，從而使得模型庫(kù)和數(shù)據(jù)庫(kù)能用統(tǒng)一的方法管理，便于模型的修改和更新；模型庫(kù)管理系統(tǒng)的功能，就是實(shí)現(xiàn)對(duì)各個(gè)模型的維護(hù)、修改和連接，以及實(shí)現(xiàn)模型庫(kù)與方法庫(kù)、數(shù)據(jù)庫(kù)及人機(jī)界面的連接。

3）方法庫(kù)及其管理系統(tǒng)。評(píng)價(jià)方法庫(kù)有多種方法；存儲(chǔ)有工藝計(jì)算方法的單元模塊和組合模塊，包括：天然氣物性參數(shù)和狀態(tài)參數(shù)、壓氣機(jī)特性參數(shù)、輸氣管熱力、水力等計(jì)算方法；方法庫(kù)管理系統(tǒng)的功能，就是實(shí)現(xiàn)對(duì)各方案評(píng)價(jià)、工藝計(jì)算方案的組合和連接，以及實(shí)現(xiàn)方法庫(kù)與數(shù)據(jù)庫(kù)、模型庫(kù)及人機(jī)界面的連接。

4）人機(jī)界面。人機(jī)界面是用于實(shí)現(xiàn)決策人員和GPPES連接的程序接口，其主要功能是：①檢查和接收決策人員輸入的數(shù)據(jù)和操作步驟，并輸出系統(tǒng)運(yùn)行后的結(jié)果，具體體現(xiàn)在數(shù)據(jù)輸入與輸出、儲(chǔ)氣調(diào)峰工藝計(jì)算、方案評(píng)價(jià)與優(yōu)選等菜單中；②協(xié)調(diào)和管理數(shù)據(jù)庫(kù)、模型庫(kù)、方法庫(kù)3大管理系統(tǒng)之間的連接和數(shù)據(jù)傳遞，保證各子系統(tǒng)之間友好協(xié)調(diào)地工作；③為決策人員提供良好而使用方便的接口功能，如操作簡(jiǎn)捷的按鈕、菜單式屏幕提示、出錯(cuò)處理和錯(cuò)誤類型提示，使GPPES成為無師自通的高級(jí)“傻瓜”軟件。

1.3 GPPES的規(guī)模

本系統(tǒng)用 MATLAB、VFP、SPS、Office聯(lián)合編程實(shí)現(xiàn)上述功能，系統(tǒng)由675個(gè)文件組成。系統(tǒng)全部運(yùn)行一遍可獲天然氣管網(wǎng)儲(chǔ)氣調(diào)峰能力方案評(píng)價(jià)成果20類報(bào)表。

2 工程實(shí)例

正在運(yùn)行的大型天然氣管網(wǎng)，主干網(wǎng)絡(luò)全長(zhǎng)2 182 km，管網(wǎng)經(jīng)過地區(qū)地形起伏較大；管網(wǎng)系統(tǒng)有25座輸氣站，3座壓氣站（A、B、C站），2座地下儲(chǔ)氣庫(kù)（D、E庫(kù)）。試?yán)肎PPES評(píng)價(jià)輸氣量為91×108m3／a的儲(chǔ)氣調(diào)峰方案。

某天然氣管網(wǎng)向輸氣站計(jì)劃供氣情況見表1。部分壓氣站投入運(yùn)行，地下儲(chǔ)氣庫(kù)不投入運(yùn)行；壓氣站出站最高壓力為9.88MPa；管網(wǎng)最低壓力點(diǎn)壓力為4.00MPa。

按照方案評(píng)價(jià)步驟，將原始數(shù)據(jù)準(zhǔn)備就緒后，送入SPS進(jìn)行仿真，仿真結(jié)果輸出后，生成儲(chǔ)氣調(diào)峰方案，再重新調(diào)整管網(wǎng)工藝參數(shù)，更新SPS軟件的數(shù)據(jù)接口文件，再次啟動(dòng)SPS，獲得另一組儲(chǔ)氣調(diào)峰方案。本工程實(shí)例中反復(fù)調(diào)整管網(wǎng)工藝參數(shù)6次，獲得6個(gè)儲(chǔ)氣調(diào)峰方案，結(jié)果見表2（篇幅所限，評(píng)價(jià)過程和各方案恕不一一列出），經(jīng)綜合評(píng)價(jià)后獲得各種評(píng)價(jià)方法的排序結(jié)果（表3）。

表1 某天然氣管網(wǎng)向各輸氣站計(jì)劃供氣情況一覽表

表2 某天然氣管網(wǎng)調(diào)峰評(píng)價(jià)預(yù)選方案一覽表

表3 某天然氣管網(wǎng)調(diào)峰方案各種方法綜合評(píng)價(jià)結(jié)果表

本系統(tǒng)中選取管網(wǎng)儲(chǔ)氣量、儲(chǔ)氣最高、最低壓力、增壓消耗功率、管輸動(dòng)力費(fèi)用、供需失衡戶數(shù)、實(shí)際供氣量、管輸收入作為儲(chǔ)氣調(diào)峰方案的評(píng)價(jià)指標(biāo)。其中，管網(wǎng)儲(chǔ)氣量、儲(chǔ)氣最高壓力、實(shí)際供氣量、管輸收入是高優(yōu)指標(biāo)；而儲(chǔ)氣最低壓力、增壓消耗功率、管輸費(fèi)用、供需失衡戶數(shù)是低優(yōu)指標(biāo)。

評(píng)價(jià)系統(tǒng)中分別采用秩和比法［15］（方法1）、累積值法［4］（方法2）、優(yōu)異度法［4］（方法3）對(duì)各儲(chǔ)氣調(diào)峰方案進(jìn)行評(píng)價(jià)排序。由于不同的評(píng)價(jià)方法獲得的評(píng)價(jià)結(jié)果有所差異，為了解決這個(gè)問題，再用最大兼容度法［16］（方法4）對(duì)前面3種方法進(jìn)行綜合評(píng)價(jià)。因此，方法4的排序結(jié)果即是綜合評(píng)價(jià)的結(jié)果。

經(jīng)評(píng)價(jià)排序表明，原第1方案為最佳儲(chǔ)氣調(diào)峰方案，管網(wǎng)儲(chǔ)氣量為1 362.476×104m3／d，儲(chǔ)氣最高壓力為9.834MPa，最低壓力為4.566MPa，增壓消耗功率為23.375MW／h、管輸動(dòng)力費(fèi)用為0.005 10元／（km3·km）、供需失衡戶數(shù)為14戶、實(shí)際供氣量為92.418×108m3／a、管輸收入為50.351億元／a。增壓所消耗的功率是啟動(dòng)的A、B號(hào)壓氣站，管網(wǎng)儲(chǔ)氣量按B號(hào)壓氣站后的管網(wǎng)長(zhǎng)度計(jì)算。

3 結(jié)論

1）全面、系統(tǒng)地提出了大型天然氣管網(wǎng)儲(chǔ)氣調(diào)峰能力評(píng)價(jià)的模式和方法。

2）通過綜合評(píng)價(jià)，實(shí)現(xiàn)了用簡(jiǎn)捷的評(píng)價(jià)模型獲得大型復(fù)雜天然氣管網(wǎng)的最佳調(diào)峰調(diào)度方案，突破了用經(jīng)典優(yōu)化方法建立、求解該管網(wǎng)優(yōu)化調(diào)度數(shù)學(xué)模型難度很大的瓶頸問題。

3）該系統(tǒng)有效地改善了輸氣調(diào)度管理決策的環(huán)境，提高了調(diào)峰調(diào)度管理決策的及時(shí)性和科學(xué)性。

4）方案生成與方案選擇自動(dòng)化。每運(yùn)行一次，自動(dòng)生成一種評(píng)價(jià)方案；對(duì)多次運(yùn)行的各方案進(jìn)行總體評(píng)價(jià)排序，并能挑選出既滿足工藝要求、滿足用戶用氣需求，又經(jīng)濟(jì)合理的調(diào)峰調(diào)度方案。

5）分別采用4種不同的評(píng)價(jià)方法實(shí)現(xiàn)了儲(chǔ)氣調(diào)峰方案的組合評(píng)價(jià)。

6）離散數(shù)據(jù)的自動(dòng)建模。參加方案評(píng)價(jià)的數(shù)據(jù)均能自動(dòng)建立方案評(píng)價(jià)的數(shù)學(xué)模型，供方案評(píng)價(jià)使用。

7）對(duì)更新工程技術(shù)人員在天然氣管網(wǎng)調(diào)峰調(diào)度管理中的觀念和模式，對(duì)優(yōu)化調(diào)峰調(diào)度運(yùn)行方案、節(jié)省經(jīng)營(yíng)管理費(fèi)，縮短調(diào)度決策周期，提高天然氣管網(wǎng)運(yùn)行質(zhì)量起到了促進(jìn)作用。

8）SPS仿真軟件利用率大大提高，在儲(chǔ)氣調(diào)峰方案生成中發(fā)揮了重要作用。

GPPES的開發(fā)是一項(xiàng)復(fù)雜的系統(tǒng)工程，有關(guān)工作剛剛起步，缺乏成熟經(jīng)驗(yàn)，本系統(tǒng)項(xiàng)目組將在實(shí)踐應(yīng)用的過程中對(duì)系統(tǒng)進(jìn)行的進(jìn)一步檢驗(yàn)、改進(jìn)和完善。

［1］姚安林，趙忠剛，張錦偉.油氣管道風(fēng)險(xiǎn)評(píng)估質(zhì)量評(píng)價(jià)技術(shù)［J］.天然氣工業(yè)，2013，33（12）：111－116.YAO Anlin，ZHAO Zhonggang，ZHANG Jinwei.Quality evaluation for oil and gas pipeline risk assessment［J］.Natural Gas Industry，2013，33（12）：111－116.

［2］吳志平，蔣宏業(yè)，李又綠，等.油氣管道完整性管理效能評(píng)價(jià)技術(shù)研究［J］.天然氣工業(yè)，2013，33（12）：131－137.Wu Zhiping，Jiang Hongye，Li Youlü，et al.Efficiency appraisal of oil ＆gas pipeline integrity management［J］.Natural Gas Industry，2013，33（12）：131－137.

［3］涂序彥，王樅，郭燕慧.大系統(tǒng)控制論［M］.北京：國(guó)防工業(yè)出版社，1994.TU Xuyan，WANG Cong，GUO Yanhui.Large systems cybernetics［M］.Beijing：National Defense University Press，1994.

［4］劉寧，王登瀛.實(shí)用多目標(biāo)分析與優(yōu)選［M］.武漢：武漢出版社，1998.LIU Ning，WANG Dengying.Practical multi－objective analysis and optimization［M］.Wuhan：Wuhan Publishing House，1998.

［5］LURIE M V.Modeling of oil product and gas pipeline transportation［M］.SINAISKI E G，Trans.Weinheim：WILEY－VCH Verlag GmbH ＆ Co.KGaA，2008.

［6］王新華，李堂軍，丁黎黎.復(fù)雜大系統(tǒng)評(píng)價(jià)理論與技術(shù)［M］.濟(jì)南：山東大學(xué)出版社，2010.WANG Xinhua，LI Tangjun，DING Lili.Evaluation theory and technology of large complex systems［M］.Jinan：Shandong University Press，2010.

［7］遲國(guó)泰，王衛(wèi).基于科學(xué)發(fā)展的綜合評(píng)價(jià)理論、方法與應(yīng)用［M］.北京：科學(xué)出版社，2009.CHI Guotai，WANG Wei.Comprehensive evaluation theory，method and application based on scientific development［M］.Beijing：Science Press，2009.

［8］高洪深.決策支持系統(tǒng)（DSS）理論與方法［M］.4版.北京：清華大學(xué)出版社，2009.GAO Hongshen.Decision support system（DSS）：Theory and method［M］.4thed.Beijing：Tsinghua University Press，2009.

［9］胡于進(jìn)，凌玲.決策支持系統(tǒng)的開發(fā)與應(yīng)用［M］.北京：機(jī)械工業(yè)出版社，2006.HU Yujin，LING Ling.Development and application of decision support systems［M］.Beijing：China Machine Press，2006.

［10］艾慕陽(yáng).大型油氣管網(wǎng)系統(tǒng)可靠性若干問題探討［J］.油氣儲(chǔ)運(yùn)，2013，32（12）：1265－1270.AI Muyang.Discussion on issues regarding the reliability of large－scale oil and gas pipeline network systems［J］.Oil ＆Gas Storage and Transportation，2013，32（12）：1265－1270.

［11］尤秋菊，樊建春，朱偉，等.天然氣管網(wǎng)系統(tǒng)風(fēng)險(xiǎn)評(píng)估［J］.油氣儲(chǔ)運(yùn)，2013，32（8）：834－839.YOU Qiuju，F(xiàn)AN Jianchun，ZHU Wei，et al.Risk assessment of natural gas pipeline networks［J］.Oil ＆ Gas Storage and Transportation，2013，32（8）：834－839.

［12］嚴(yán)銘卿，廉樂明.天然氣輸配工程［M］.北京：中國(guó)建筑工業(yè)出版社，2005.YAN Mingqing，LIAN Leming.Natural gas transmission and distribution engineering［M］.Beijing：China Architecture ＆Building Press，2005.

［13］宋葉志，賈東永.MATLAB數(shù)值分析與應(yīng)用［M］.北京：機(jī)械工業(yè)出版社，2009.SONG Yezhi，JIA Dongyong.MATLAB numerical analysis and application［M］.Beijing：China Machine Press，2009.

［14］ABBASPOUR M，KRISHNASWAMI P，CHAPMAN K S.Transient optimization in natural gas compressor stations for linepack operation［J］.Journal of Energy Resources Technology，2007，129（4）：314－324.

［15］陳正偉.綜合評(píng)價(jià)技術(shù)及應(yīng)用［M］.成都：西南財(cái)經(jīng)大學(xué)出版社，2013.CHEN Zhengwei.Technology and applications of comprehensive evaluation［M］.Chengdu：Press of Southwest Finance ＆ Economics University，2013.

［16］秦壽康.綜合評(píng)價(jià)原理與應(yīng)用［M］.北京：電子工業(yè)出版社，2003.QIN Shoukang.Principles and applications of comprehensive evaluation［M］.Beijing：Publishing House of Electronics Industry，2003.