高 龍，王 迪，車(chē)德勇，李少華，白 章

(1.東北電力大學(xué)能源與動(dòng)力工程學(xué)院，吉林吉林132012;2.中國(guó)大唐集團(tuán)科學(xué)技術(shù)研究院，北京100033;3.中國(guó)科學(xué)院工程熱物理研究所，北京100190)

火電廠是復(fù)雜能量系統(tǒng)，運(yùn)用傳統(tǒng)的能量計(jì)算及分析方法［1－2］，不但計(jì)算量大而且不能準(zhǔn)確分析系統(tǒng)變量之間的關(guān)系。近年來(lái)，以定量和定性相結(jié)合的方法對(duì)大系統(tǒng)作全面的模擬、分析取得了較好的效果［3－5］。其中，應(yīng)用過(guò)程系統(tǒng)工程的模擬方法［6］，確定系統(tǒng)各部位的信息并推定系統(tǒng)特性，在復(fù)雜的化工系統(tǒng)中應(yīng)用廣泛。本文以某200 MW熱電聯(lián)產(chǎn)系統(tǒng)為研究對(duì)象，應(yīng)用通用型模擬軟件Aspen Plus對(duì)發(fā)電系統(tǒng)中多循環(huán)嵌套的復(fù)雜網(wǎng)絡(luò)進(jìn)行建模研究，并對(duì)系統(tǒng)性能進(jìn)行計(jì)算分析。

1 系統(tǒng)模擬原理及方法

1.1 系統(tǒng)結(jié)構(gòu)分析

設(shè)備數(shù)量巨大、多循環(huán)結(jié)構(gòu)復(fù)雜的網(wǎng)絡(luò)系統(tǒng)，需要數(shù)以萬(wàn)計(jì)的各類(lèi)偏微分、線性及其它形式的方程組進(jìn)行數(shù)學(xué)描述，方程組共同求解對(duì)于常規(guī)方法難以實(shí)現(xiàn)，并且存在計(jì)算效率及精度等問(wèn)題。因此，通過(guò)系統(tǒng)結(jié)構(gòu)分析將復(fù)雜網(wǎng)絡(luò)系統(tǒng)分解成多個(gè)相對(duì)獨(dú)立的子系統(tǒng)，以降低系統(tǒng)計(jì)算維度及收斂迭代次數(shù)，通常包括系統(tǒng)結(jié)構(gòu)的數(shù)學(xué)描述、系統(tǒng)分隔、不可再分塊的切斷以及計(jì)算次序的確定四個(gè)主要步驟［7］，如圖1所示。

圖1 系統(tǒng)結(jié)構(gòu)分析過(guò)程

(1)系統(tǒng)結(jié)構(gòu)描述

實(shí)際的大型復(fù)雜系統(tǒng)都是由眾多設(shè)備和管路組成的精密系統(tǒng)，但其中許多輔助設(shè)備及管路對(duì)系統(tǒng)整體影響不大，系統(tǒng)越復(fù)雜越不利于計(jì)算機(jī)對(duì)數(shù)據(jù)的讀取及整體的運(yùn)算，因此對(duì)復(fù)雜多循環(huán)網(wǎng)絡(luò)系統(tǒng)進(jìn)行歸納、簡(jiǎn)化和數(shù)學(xué)描述是必要的。鑒于此，運(yùn)用離散數(shù)學(xué)中的圖論方法，將復(fù)雜的網(wǎng)絡(luò)結(jié)構(gòu)、拓?fù)潢P(guān)系抽象成圖，并作為一定邏輯關(guān)系的特定表達(dá)形式［8］。以單元設(shè)備或虛擬單元為結(jié)點(diǎn)，將對(duì)系統(tǒng)影響不大的設(shè)備簡(jiǎn)化，對(duì)代表工質(zhì)流動(dòng)和能量交換的過(guò)程以符號(hào)標(biāo)示，利用標(biāo)有符號(hào)點(diǎn)和線的對(duì)系統(tǒng)設(shè)備、工質(zhì)及能量關(guān)系進(jìn)行數(shù)學(xué)描述。

(2)系統(tǒng)分隔

系統(tǒng)分割過(guò)程中相對(duì)獨(dú)立的分離復(fù)雜的循環(huán)結(jié)構(gòu)，并將其分解至不可繼續(xù)拆分的程度，該過(guò)程不僅劃分了系統(tǒng)計(jì)算的層次性，也是數(shù)學(xué)計(jì)算中的一個(gè)降階過(guò)程。升冪法、鄰接矩陣通路搜索等方法是目前系統(tǒng)分隔過(guò)程常用的方法，其中鄰接矩陣通路搜索法在通用型模擬系統(tǒng)中應(yīng)用廣泛，通過(guò)相對(duì)復(fù)雜的算法編程實(shí)現(xiàn)復(fù)雜系統(tǒng)的自動(dòng)分隔計(jì)算，可在計(jì)算過(guò)程中節(jié)省大量的存儲(chǔ)空間，詳細(xì)計(jì)算框圖如圖2所示。

(3)不可再分環(huán)路的切斷與排序

圖2 鄰接矩陣通路搜索法邏輯框圖

若不可再分的子系統(tǒng)是由若干單元和環(huán)路組成，需斷開(kāi)部分外部變量，以方便整體的線性化，此時(shí)迭代變量維數(shù)大大降低，在得到切斷變量的計(jì)算之后，應(yīng)用迭代算法進(jìn)行迭代計(jì)算，從數(shù)學(xué)角度出發(fā)，該過(guò)程也是對(duì)系統(tǒng)整體計(jì)算的進(jìn)一步降階和簡(jiǎn)化。在不可再分環(huán)路的計(jì)算順序上，被切斷流股視為已知，輸入部分全部已知的單元及其輸出部分為已知，并且優(yōu)先計(jì)算該單元，然后使用上述原則繼續(xù)尋找下一個(gè)可計(jì)算單元，如此重復(fù)至到回到切斷位置，不可在分塊M1及計(jì)算次序如圖3－圖4所示。

圖3 不可再分塊M1

圖4 M1的計(jì)算次序

1.2 單元過(guò)程模型化

系統(tǒng)模擬的數(shù)學(xué)模型由兩部分組成，即系統(tǒng)所包含的單元數(shù)學(xué)模型和描述系統(tǒng)中各單元之間拓?fù)潢P(guān)系的結(jié)構(gòu)數(shù)學(xué)模型。系統(tǒng)流程的單元數(shù)學(xué)模擬屬于混合型模型，在建模過(guò)程中通過(guò)做出適當(dāng)?shù)募俣ê秃?jiǎn)化，不僅易于數(shù)學(xué)描述、求解方便有效，又能反映單元過(guò)程的基本規(guī)律。在對(duì)客觀過(guò)程充分理解的基礎(chǔ)上，對(duì)單元過(guò)程建立物理模型，然后通過(guò)推導(dǎo)加以數(shù)學(xué)描述，其基本步驟如表1所示。

表1 單元過(guò)程數(shù)學(xué)建模方法

1.3 系統(tǒng)模擬求解方法

系統(tǒng)模擬的實(shí)質(zhì)是求解系統(tǒng)中全部的非線性方程組，若采用一般性算法對(duì)該組方程進(jìn)行求解會(huì)造成代數(shù)和微分模型混合、高度非線性等復(fù)雜情況，不利于高效準(zhǔn)確的獲得系統(tǒng)參數(shù)，而序貫?zāi)K法的應(yīng)用能夠有效的解決以上問(wèn)題，通過(guò)對(duì)單元過(guò)程數(shù)學(xué)模型的整合，再將各單元過(guò)程以一定的次序聯(lián)接，序貫計(jì)算系統(tǒng)數(shù)學(xué)模型。系統(tǒng)模擬的描述性方程組包括:單元模型方程、流股連接方程、物性估算方程和設(shè)計(jì)規(guī)定方程，序貫?zāi)K法求解系統(tǒng)模擬問(wèn)題和設(shè)計(jì)型問(wèn)題時(shí)的迭代收斂過(guò)程包括:基本參數(shù)下的工質(zhì)物性估算迭代收斂過(guò)程，單元模塊變化過(guò)程計(jì)算的迭代收斂過(guò)程，子系統(tǒng)循環(huán)計(jì)算過(guò)程的迭代收斂過(guò)程和計(jì)算符合設(shè)計(jì)規(guī)定值的迭代收斂過(guò)程［9］。其中，子系統(tǒng)循環(huán)迭代收斂過(guò)程為系統(tǒng)模擬中的重要環(huán)節(jié)，對(duì)迭代變量進(jìn)行修正的基本方法包括:直接迭代法、Wegstein法、Newton法和綜合算法等。

2 火電機(jī)組熱力系統(tǒng)建模與仿真

火電廠的汽輪機(jī)熱力系統(tǒng)是具有代表的復(fù)雜能量轉(zhuǎn)換系統(tǒng)［10］。本文以某熱電聯(lián)產(chǎn)機(jī)組的汽輪機(jī)熱力系統(tǒng)為研究對(duì)象，該系統(tǒng)包括一次中間再熱的汽輪機(jī)本體、凝汽器、低壓加熱器、除氧器和高壓加熱器等主要設(shè)備。

2.1 系統(tǒng)建模

(1)系統(tǒng)結(jié)構(gòu)分析

在該系統(tǒng)中將鍋爐系統(tǒng)以?xún)膳_(tái)換熱器的形式簡(jiǎn)化，系統(tǒng)模擬完成后通過(guò)熱流與鍋爐系統(tǒng)耦合;該機(jī)組汽輪機(jī)為三缸結(jié)構(gòu)共32級(jí)、8段抽汽，考慮到缸體內(nèi)抽汽級(jí)處蒸汽的狀態(tài)變化，在模擬過(guò)程中將實(shí)際缸體按抽汽級(jí)數(shù)繼續(xù)按模塊劃分;將多臺(tái)凝結(jié)水泵或給水泵分別以一個(gè)模塊進(jìn)行模擬。經(jīng)過(guò)結(jié)構(gòu)分析及簡(jiǎn)化后的熱力系統(tǒng)如圖5所示。

圖5 熱力系統(tǒng)基本結(jié)構(gòu)信息圖

(2)物性方法及計(jì)算模型選擇

由于汽輪機(jī)熱力系統(tǒng)中工質(zhì)只有水和水蒸汽，因此在蒸汽循環(huán)的物性方法中Aspen Plus提供了兩種方法，分別是STEAMNBS和STEAM－TA，本文經(jīng)過(guò)對(duì)比選擇了STEAMNBS方法，該方法采用了1984年美國(guó)國(guó)家標(biāo)準(zhǔn)局(NBS)和科學(xué)研究委員會(huì)(NRC)制定的水蒸汽相關(guān)熱力學(xué)性質(zhì)，水蒸汽性能?chē)?guó)際學(xué)會(huì)(IAPS)制定的水蒸汽相關(guān)傳遞特性，適應(yīng)的溫度范圍約是0℃至1 700℃，最大壓力為1 000 MPa。

(3)單元模塊選擇

根據(jù)系統(tǒng)工程理論的模塊化建模方法，可以將該系統(tǒng)中的設(shè)備分為三大類(lèi):壓力變送設(shè)備、換熱設(shè)備及相關(guān)輔助設(shè)備，具體模塊及功能如表2所示。

表2 汽輪機(jī)熱力系統(tǒng)單元模塊分類(lèi)及應(yīng)用

(4)多循環(huán)系統(tǒng)建模方法

由于該系統(tǒng)多循環(huán)嵌套的復(fù)雜特性，模擬過(guò)程不宜采用整體建模計(jì)算的方法。通過(guò)結(jié)構(gòu)分析將汽輪機(jī)高中低壓缸按抽汽位置被劃分為9個(gè)模塊，汽輪機(jī)本體模塊為純序貫結(jié)構(gòu)，不存在循環(huán)及收斂過(guò)程，因此首先應(yīng)用Comper、Valve、FSplit等模塊對(duì)汽輪機(jī)本體進(jìn)行建模，其中Comper模塊提供了包括標(biāo)準(zhǔn)等熵模型(Isentropic)、ASME模型、GPSA模型等8種計(jì)算模型，排氣壓力(Discharge Pressure)、壓力比率(Pressure Ratio)等模型參數(shù)。本文選擇了標(biāo)準(zhǔn)等熵模型，通過(guò)設(shè)置排氣壓力、等熵效率和機(jī)械效率對(duì)排氣溫度、輸出功率等參數(shù)進(jìn)行計(jì)算，在汽輪機(jī)本體模型構(gòu)建完成后計(jì)算發(fā)現(xiàn)在功率和抽汽參數(shù)方面略有偏差，原因是尚未進(jìn)行回?zé)峒訜嵯到y(tǒng)的耦合，耦合后通過(guò)各層循環(huán)的收斂計(jì)算會(huì)對(duì)汽輪機(jī)本體參數(shù)進(jìn)行修正;另一方面，由于汽輪機(jī)調(diào)節(jié)級(jí)與高壓缸內(nèi)其它級(jí)組在等熵?zé)峤怠⑺俣缺?、反?dòng)度等方面都存在很大差異，因此將調(diào)節(jié)級(jí)單獨(dú)劃分模塊進(jìn)行模擬，避免影響高壓缸整體精度;最后，在低壓缸排汽位置的模擬過(guò)程中，需要在收斂計(jì)算的有效相態(tài)(Valid Phases)選項(xiàng)中選擇汽液共存(Vapor－Liquid)的狀態(tài)。

通過(guò)鄰接矩陣通路搜索法確定系統(tǒng)分隔流股為各回?zé)峒訜崞鞯氖杷鞴?，因此在?yīng)用MHeatX模塊對(duì)回?zé)峒訜崞鬟M(jìn)行模擬時(shí)，采用逆向單循環(huán)的方法，如圖6所示，即H1+JG2+B1→H2+JG1+JG2+B1+H1……以此類(lèi)推。在應(yīng)用MHeatX對(duì)加熱器加熱過(guò)程進(jìn)行單元模擬時(shí)，需要應(yīng)用設(shè)計(jì)規(guī)定(Design Spec)方法，在保證所求流股在允許偏差內(nèi)的情況下對(duì)該級(jí)抽汽量進(jìn)行修正計(jì)算，計(jì)算下一循環(huán)時(shí)先將上一循環(huán)的設(shè)計(jì)規(guī)定設(shè)置為隱藏(Hide)狀態(tài)，本循環(huán)收斂后再將上一循環(huán)的設(shè)計(jì)規(guī)定恢復(fù)進(jìn)行二次核算，以保證計(jì)算精度和計(jì)算效率。以汽輪機(jī)本體與JG2耦合為例，在模擬JG2的過(guò)程中將疏水S1與下一循環(huán)斷開(kāi)，給水GS1作為所求的目標(biāo)流股，通過(guò)對(duì)GS2的詳細(xì)參數(shù)和S1目標(biāo)參數(shù)的設(shè)定，所求得的GS1溫度與設(shè)計(jì)溫度相差0.1℃。建立設(shè)計(jì)規(guī)定DS－1，設(shè)定溫度允許偏差為0.05℃，調(diào)整變量為C1流量，通過(guò)系統(tǒng)迭代計(jì)算將C1流量修正0.002 t/h。

圖6 逆向單循環(huán)法示意圖

2.2 模擬誤差分析

高壓缸中調(diào)節(jié)級(jí)內(nèi)效率與2至12級(jí)變化幅度很大，因此單獨(dú)模擬。模擬結(jié)果如圖7所示，內(nèi)效率趨勢(shì)與設(shè)計(jì)數(shù)據(jù)吻合，各模塊平均相對(duì)誤差為0.004 87，其中，由于27級(jí)為汽輪機(jī)排氣，壓力由第26級(jí)的0.0185 MPa降至背壓狀態(tài)的0.004 9 MPa，并且工質(zhì)處于汽液兩相狀態(tài)，造成模擬誤差增大到0.041，但對(duì)整體趨勢(shì)及功率影響不大;在汽輪機(jī)功率的模擬結(jié)果中，各模塊功率平均相對(duì)誤差為0.006 52，如圖8所示，調(diào)節(jié)級(jí)結(jié)構(gòu)功能特殊，工質(zhì)流量、壓力、溫度等在調(diào)節(jié)級(jí)變化都很大，設(shè)計(jì)等熵?zé)峤禐?9.28 kJ/kg，其余級(jí)均為36 kJ/kg左右，導(dǎo)致相對(duì)誤差達(dá)到0.032，但對(duì)整體趨勢(shì)影響不大。

圖7 額定工況汽輪機(jī)內(nèi)效率模擬誤差分析

圖8 額定工況汽輪機(jī)各模塊功率模擬誤差分析

本系統(tǒng)為定壓模擬，涉及到循環(huán)收斂計(jì)算的抽汽系統(tǒng)與加熱器疏水系統(tǒng)溫度與流量會(huì)存在誤差，設(shè)定汽輪機(jī)模塊的排汽壓力，通過(guò)調(diào)整內(nèi)效率調(diào)整各模塊的抽汽溫度，各模塊內(nèi)效率發(fā)生改變的同時(shí)引起汽輪機(jī)功率的變化，對(duì)各工況進(jìn)行模擬后，汽輪機(jī)功率的相對(duì)誤差最大值為0.019 86，平均相對(duì)誤差為0.007 43，如圖9所示。

圖9 不同工況汽輪機(jī)整體功率模擬誤差分析

在回?zé)峒訜嵯到y(tǒng)中，為了調(diào)整回?zé)峒訜崞鞯恼羝肟跍囟?、疏水溫度、上下端差等參?shù)更接近于設(shè)計(jì)值，采用設(shè)計(jì)規(guī)定的方法依次調(diào)節(jié)各子系統(tǒng)的循環(huán)流量，通過(guò)整體調(diào)整后抽汽溫度平均相對(duì)誤差0.000 38，抽汽流量平均相對(duì)誤差0.013 11，疏水溫度平均相對(duì)誤差0.005 34，疏水流量平均相對(duì)誤差0.010 07，不同工況模擬抽汽溫度及流量如表3所示。

表3 不同工況抽汽溫度及流量模擬參數(shù)

3 模擬結(jié)果分析

(1)純凝工況機(jī)組經(jīng)濟(jì)性評(píng)價(jià)

通過(guò)系統(tǒng)流程模擬、耦合及收斂計(jì)算，得出不同工況條件下系統(tǒng)各點(diǎn)的模擬參數(shù)。以熱力學(xué)第一定律為基礎(chǔ)，在Aspen Plus中運(yùn)用Fortran語(yǔ)言編程，用熱量法定量計(jì)算該熱電廠在純凝工況下熱經(jīng)濟(jì)性指標(biāo)以評(píng)價(jià)其熱經(jīng)濟(jì)性。如圖10所示，選取6組純凝工況，模擬計(jì)算其汽輪機(jī)絕對(duì)內(nèi)效率ηi、鍋爐效率ηb、總效率ηcp等熱經(jīng)濟(jì)指標(biāo)，其中選取部分管道計(jì)算管道效率，平均管道效率均在0.98－0.99之間;汽輪機(jī)絕對(duì)內(nèi)效率直接影響總效率的變化情況，并隨負(fù)荷變化較大，負(fù)荷位于75%以上時(shí)變化趨于平穩(wěn);鍋爐效率在純凝工況隨負(fù)荷降低而略微減小，最低效率不低于0.9。

圖10 機(jī)組主要效率隨負(fù)荷變化規(guī)律

汽輪發(fā)電機(jī)組的汽耗率d和熱耗率q為評(píng)價(jià)機(jī)組運(yùn)行經(jīng)濟(jì)性的重要指標(biāo)，其中q的大小主要取決于汽輪機(jī)絕對(duì)內(nèi)效率ηi，這兩項(xiàng)指標(biāo)均是發(fā)電廠重要熱經(jīng)濟(jì)性指標(biāo)，從圖11中可以看出該機(jī)組在純凝工況下隨著機(jī)組負(fù)荷的增加汽耗率先降低后升高，最高汽耗率在30%負(fù)荷時(shí)，汽輪發(fā)電機(jī)組汽耗率為3.1375 kg/(Kw·h)，最低汽耗率在大約75% 負(fù)荷2.8 726 kg/(Kw·h);汽輪發(fā)電機(jī)組的熱耗率q隨負(fù)荷的增加逐漸下降，最高熱耗率在30%負(fù)荷時(shí)，汽輪發(fā)電機(jī)組熱耗率為9442.51 kJ/(Kw·h)，最低熱耗率為額定狀態(tài)下的8 221.76 kJ/(Kw·h)。

圖11 純凝工況汽耗率和熱耗率隨負(fù)荷變化規(guī)律

圖12 發(fā)電廠熱耗量與熱耗率隨負(fù)荷變化規(guī)律

標(biāo)準(zhǔn)煤的低位發(fā)熱量為29 270 KJ/kg，該機(jī)組設(shè)計(jì)煤種低位發(fā)熱量17 997 KJ/kg，發(fā)電廠純凝工況熱耗量及熱耗率隨負(fù)荷變化趨勢(shì)如圖12所示，隨著機(jī)組負(fù)荷的增加全廠熱耗量逐漸增大，熱耗率逐漸降低;在煤耗量與煤耗率的計(jì)算中，隨著機(jī)組負(fù)荷的增加煤耗量逐漸增加、煤耗率逐漸下降，全廠供電標(biāo)準(zhǔn)煤耗率為0.3 kg/(Kw·h)左右，該機(jī)組在該煤種下的煤耗率與供電標(biāo)準(zhǔn)煤耗率相差約0.13 kg/(Kw·h)，如圖13所示。

(2)采暖工況機(jī)組經(jīng)濟(jì)性分析

圖13 發(fā)電廠煤耗量與煤耗率隨負(fù)荷變化規(guī)律

在供熱狀態(tài)下，本文模擬6個(gè)采暖抽汽工況包括:額定工業(yè)、采暖抽汽工況，最大工業(yè)、采暖抽汽工況，最大采暖抽汽工況，最大工業(yè)抽汽工況，額定工業(yè)、采暖抽汽(0.23 MPa)工況，額定工業(yè)、采暖抽汽(0.32 MPa)工況。在各工況下汽輪機(jī)一部分蒸汽作為工業(yè)抽汽和采暖抽汽，新汽量增加到640 t/h至670 t/h，汽耗率均高于純凝發(fā)電狀態(tài)，在額定工業(yè)、采暖抽汽工況下汽耗率為4.351 8 kg/(Kw·h);由于熱電聯(lián)產(chǎn)蒸汽流股對(duì)熱能的充分利用，在熱耗率方面均低于純凝工況，在額定工業(yè)、采暖抽汽工況下熱耗率為5 117.09 KJ/(Kw·h)，其它采暖抽汽工況主要參數(shù)如表4所示，其中除4號(hào)(最大工業(yè)抽汽工況)熱耗率與汽耗率與其它工況差別較大，其余工況的熱耗率及汽耗率均在5 000 KJ/(Kw·h)和4.2 Kg/(Kw·h)左右，如圖14所示;各采暖工況供給熱用戶的熱耗約為750 GJ/h，燃料利用系數(shù)為0.75，其中不同工況的對(duì)比情況如圖15所示。

圖14 不同采暖抽汽工況熱耗率及汽耗率情況

圖15 不同采暖工況供熱熱量及燃料利用系數(shù)

表4 采暖抽汽工況機(jī)組主要參數(shù)對(duì)比

熱電聯(lián)產(chǎn)系統(tǒng)總熱耗的傳統(tǒng)分配方法主要包括熱量法、實(shí)際焓降法和火用方法三類(lèi)方法。其中熱量法以只考慮能量的數(shù)量，不考慮能量的質(zhì)量差別的特點(diǎn)，使發(fā)電部分的節(jié)能效果更為明顯;實(shí)際焓降法是按聯(lián)產(chǎn)汽流在汽輪機(jī)中做功的比例來(lái)分配，由發(fā)電方面承擔(dān)系統(tǒng)的損失，供熱方面的節(jié)能效果更為明顯;火用方法，同時(shí)考慮熱能的質(zhì)量和數(shù)量上的差別，將熱耗量按最大做功能力比例分配，具有較為完善的熱力學(xué)理論基礎(chǔ)，以三種分配方式的計(jì)算方法計(jì)算最大采暖抽汽工況如表5所示，其中從三種方法對(duì)比，火用方法折中并更為接近實(shí)際情況。

表5 熱電聯(lián)產(chǎn)系統(tǒng)總熱耗的分配方法對(duì)比

4 結(jié) 論

本文針對(duì)復(fù)雜能量系統(tǒng)的全局計(jì)算和模擬問(wèn)題進(jìn)行了研究，以某200 MW熱電聯(lián)產(chǎn)機(jī)組熱力系統(tǒng)為研究對(duì)象，通過(guò)應(yīng)用多循環(huán)系統(tǒng)的建模方法和整體求解方法，實(shí)現(xiàn)了復(fù)雜系統(tǒng)的整體定量計(jì)算，與傳統(tǒng)方法相比提高了計(jì)算效率和計(jì)算精度。通過(guò)對(duì)研究對(duì)象的誤差分析和模擬計(jì)算分析，實(shí)現(xiàn)了系統(tǒng)局部關(guān)鍵參數(shù)的判斷與預(yù)測(cè)以及系統(tǒng)整體綜合狀態(tài)評(píng)價(jià)，該模擬方法在火電機(jī)組和其它復(fù)雜能量轉(zhuǎn)換系統(tǒng)中的應(yīng)用，可在節(jié)能改造、經(jīng)濟(jì)性分析、故障診斷等方面起到重要作用。

［1］李勇，李仁杰，曹麗華，等.凝汽式汽輪機(jī)低壓缸相對(duì)內(nèi)效率在線檢測(cè)方法［J］.化工機(jī)械，2014，41(3):316－322.

［2］呂國(guó)強(qiáng)，王華，馬文會(huì)，等.小龍?zhí)峨姀S300MW機(jī)組熱力系統(tǒng)火用分析［J］.動(dòng)力工程學(xué)報(bào)，2011，31(2):85－89.

［3］韓中合，王營(yíng)營(yíng)，王繼選，等.碳補(bǔ)集系統(tǒng)與燃煤機(jī)組熱力系統(tǒng)耦合的熱經(jīng)濟(jì)性分析［J］.化工進(jìn)展，2014，33(6):1616－1623.

［4］陳巒，黃琦，陳池.基于目標(biāo)參數(shù)和儲(chǔ)能容量約束的風(fēng)電功率平滑策略［J］.中國(guó)電力，2013，46(6):63－69.

［5］胡劍輝，林汝謀.總能系統(tǒng)模塊化建模及新算法的探討［J］.燃?xì)廨啓C(jī)技術(shù)，1996，9(4):14－20.

［6］趙淵，周念成，謝開(kāi)貴，等.大電力系統(tǒng)可靠性評(píng)估的系統(tǒng)狀態(tài)抽取方法研究［J］.中國(guó)電力，2006，39(6):5－9.

［7］陸恩錫，張慧娟.化工過(guò)程模擬—原理與應(yīng)用［M］.北京:化學(xué)工業(yè)出版社，2011(3):5－6.

［8］冉鵬，李庚生，張樹(shù)芳，等.基于圖論的火電機(jī)組熱經(jīng)濟(jì)性定量分析方法［J］.中國(guó)電機(jī)工程學(xué)報(bào)，2012(14):20－26.

［9］陸恩錫，張慧娟.化工過(guò)程模擬—原理與應(yīng)用［M］.北京:化學(xué)工業(yè)出版社，2011(3):5－6.

［10］金建國(guó)，商建波，劉歡，曹麗華.汽輪機(jī)相對(duì)內(nèi)效率測(cè)量及計(jì)算方法研究［J］.東北電力大學(xué)學(xué)報(bào)，2012，32(6):72－76.