張健,劉剛,廖勝明
中南大學(xué) 能源科學(xué)與工程學(xué)院,湖南 長(zhǎng)沙 410083
地?zé)崮苡捎趦?chǔ)量大、分布廣等優(yōu)勢(shì),在眾多清潔可持續(xù)能源中備受關(guān)注[1?2]。增強(qiáng)型地?zé)嵯到y(tǒng)(enhanced geothermal system,EGS)是指用人工工程形成的裂隙,從低滲透性的高溫?zé)釒r中,經(jīng)濟(jì)地采出熱能進(jìn)行發(fā)電[3]。EGS產(chǎn)熱溫度通??蛇_(dá)150~350 ℃,用于發(fā)電具有高能量轉(zhuǎn)換效率的特點(diǎn)[4?5],開(kāi)發(fā)潛力巨大[6?8]。
目前地?zé)岚l(fā)電技術(shù)主要包括干蒸汽發(fā)電、閃蒸發(fā)電、雙工質(zhì)循環(huán)發(fā)電和卡琳娜循環(huán)發(fā)電等,其中發(fā)電效率較高的Kalina循環(huán),受到國(guó)內(nèi)外學(xué)者的廣泛關(guān)注[9?11],很多學(xué)者嘗試研究其改進(jìn)方向和衍生系統(tǒng),以進(jìn)一步提升發(fā)電效率和系統(tǒng)性能。牛成珂等[12]將吸收式熱泵耦合到Kalina系統(tǒng)中,提出吸收升溫Kalina循環(huán)系統(tǒng),并搭建了吸收升溫試驗(yàn)臺(tái),使系統(tǒng)凈發(fā)電量提升了6.8%;郭占偉[13]提出了一種可實(shí)現(xiàn)熱電聯(lián)供的氨水卡琳娜?朗肯循環(huán)組合系統(tǒng)和余熱下的雙壓力蒸發(fā)卡琳娜循環(huán)系統(tǒng),并分別對(duì)系統(tǒng)進(jìn)行了研究分析;婁聚偉等[14]建立了帶一個(gè)回?zé)崞鞯腒alina循環(huán)系統(tǒng)穩(wěn)態(tài)數(shù)學(xué)模型,通過(guò)模擬計(jì)算,分析了熱源溫度等參數(shù)對(duì)系統(tǒng)循環(huán)凈功和熱效率的影響。
在EGS高熱源溫度的條件下,直接采用Kalina循環(huán)發(fā)電,地?zé)崴毓鄿囟雀?,地?zé)崮芾貌怀浞侄鴰?lái)能源的浪費(fèi)。本文提出了兩級(jí)閃蒸?Kalina 循環(huán)聯(lián)合發(fā)電系統(tǒng)(double flash?Kalina cycle system,DFKCS),將兩級(jí)閃蒸系統(tǒng)和Kalina循環(huán)系統(tǒng)有機(jī)耦合,有效發(fā)揮閃蒸系統(tǒng)適于高溫?zé)嵩春蚄alina循環(huán)適于中低溫?zé)嵩吹膬?yōu)勢(shì),實(shí)現(xiàn)地?zé)崮艿亩嗉?jí)利用。本文將會(huì)對(duì)聯(lián)合發(fā)電系統(tǒng)建模并計(jì)算,并研究不同參數(shù)對(duì)系統(tǒng)熱力性能的影響。
DFKCS系統(tǒng)的兩級(jí)閃蒸部分包括依次連接的汽水分離器、二級(jí)閃蒸器、混壓式汽輪機(jī)和第一冷凝系統(tǒng)。Kalina循環(huán)部分包括連為一體的發(fā)生器、分離器、高溫回?zé)崞?、?jié)流閥、背壓式透平、低溫回?zé)崞骱偷诙淠到y(tǒng),Kalina循環(huán)的發(fā)生器串接在二級(jí)閃蒸器之后,再利用后的地?zé)崤潘c兩級(jí)閃蒸部分經(jīng)冷凝后的乏汽溶液混合后回灌到注入井中,如圖1所示。
圖1 兩級(jí)閃蒸?Kalina 循環(huán)聯(lián)合系統(tǒng)原理
聯(lián)合系統(tǒng)熱力循環(huán)過(guò)程如圖2所示,從生產(chǎn)井抽出的地?zé)崴?經(jīng)過(guò)2次閃蒸分離過(guò)程,汽水分離器產(chǎn)生的地?zé)嵴羝?推動(dòng)高壓汽輪機(jī)做功,二級(jí)閃蒸器產(chǎn)生的地?zé)嵴羝?與高壓汽輪機(jī)乏汽5混合形成混合蒸汽9,用于推動(dòng)低壓汽輪機(jī)做功發(fā)電,剩下的飽和地?zé)崴?進(jìn)入到發(fā)生器,用來(lái)加熱氨水溶液工質(zhì),使之蒸發(fā)成蒸汽14推動(dòng)透平做功發(fā)電,從而實(shí)現(xiàn)能量的多級(jí)利用。
圖2 兩級(jí)閃蒸?Kalina 循環(huán)聯(lián)合系統(tǒng) T?S 圖
為了分析DFKCS系統(tǒng)的影響參數(shù)和熱力性能,本文將循環(huán)中的部分參數(shù)設(shè)定為常數(shù),需要做優(yōu)化研究的參數(shù)在合理范圍內(nèi)取值,初始參數(shù)如表1所示。
本文以地?zé)崃黧w和氨水工質(zhì)為研究對(duì)象,為了簡(jiǎn)化計(jì)算,模型建立過(guò)程中作出如下假設(shè):1)系統(tǒng)在穩(wěn)定條件下運(yùn)行,所有部件均處于熱力學(xué)穩(wěn)態(tài)平衡狀態(tài);2)忽略循環(huán)中管道以及兩系統(tǒng)連接部分的熱損和壓損;3)忽略地?zé)崴畯牡叵碌降厣线^(guò)程中的焓值變化;4)忽略地?zé)崴诜蛛x和冷凝過(guò)程中的壓力和溫度損失;5)混壓汽輪機(jī)內(nèi)蒸汽的混合為等壓過(guò)程;6)發(fā)生器與冷凝器的出口假設(shè)為飽和氨水溶液,不存在過(guò)冷度;7)工質(zhì)泵前后的溶液溫度相等,節(jié)流閥前后的溶液焓值相等;8)忽略地?zé)崴┧眉袄鋮s塔水泵的耗功。
針對(duì)兩級(jí)閃蒸?Kalina循環(huán)聯(lián)合發(fā)電系統(tǒng)的循環(huán)特點(diǎn),在建模過(guò)程中,主要采用地?zé)崴|(zhì)量守恒、能量守恒以及氨水質(zhì)量守恒、能量守恒、氨質(zhì)量分?jǐn)?shù)守恒等控制方程,來(lái)描述循環(huán)狀態(tài)點(diǎn)和各部件的運(yùn)行情況。
穩(wěn)態(tài)運(yùn)行時(shí),系統(tǒng)主要的控制方程如下:質(zhì)量守恒方程:
根據(jù)以上守恒定律和對(duì)循環(huán)系統(tǒng)的理論分析,列出系統(tǒng)主要換熱部件的控制方程如表2、3所示。在系統(tǒng)模擬中,閃蒸器中的閃蒸和節(jié)流閥中的節(jié)流都是等焓過(guò)程,分離器中的分離是等壓過(guò)程,混壓式汽輪機(jī)和透平做功的理想狀態(tài)均為等熵過(guò)程。
表2 兩級(jí)閃蒸部分控制方程
表3 Kalina 循環(huán)部分控制方程
公式(1)~(19)中,m為質(zhì)量流量,kg/s;h為比焓,kJ/kg;x為氨水質(zhì)量濃度;WHPT、WLPT和 WKT分別為兩級(jí)閃蒸系統(tǒng)高壓汽輪機(jī)、低壓汽輪機(jī)和Kalina系統(tǒng)透平做功,kW;QDF,con和 QK,con分別為兩級(jí)閃蒸系統(tǒng)和Kalina系統(tǒng)的冷凝熱量,kJ;Qgeo為地?zé)崴峁┙oKalina系統(tǒng)的熱量,kW;Wpump為工質(zhì)泵消耗的功率,kW;數(shù)字下標(biāo)與圖1所示狀態(tài)點(diǎn)相對(duì)應(yīng)。
DFKCS系統(tǒng)在穩(wěn)定運(yùn)行條件下,根據(jù)系統(tǒng)各部件的控制方程,以及氨水混合物的物性數(shù)據(jù)庫(kù),組成非線性方程組,通過(guò)迭代計(jì)算得出模擬計(jì)算結(jié)果。
由于兩級(jí)閃蒸系統(tǒng)和Kalina系統(tǒng)的運(yùn)行相對(duì)獨(dú)立,因此在模型驗(yàn)證中,可以分別驗(yàn)證2個(gè)原有系統(tǒng)建模的正確性。本文通過(guò)采用文獻(xiàn)[15]和[16]的研究成果,設(shè)置相同的環(huán)境條件和初始參數(shù),將系統(tǒng)模擬結(jié)果分別與相應(yīng)數(shù)據(jù)進(jìn)行對(duì)比分析,以驗(yàn)證DFKCS系統(tǒng)模型的正確性。具體數(shù)據(jù)列于表4、5中。
表4 兩級(jí)閃蒸系統(tǒng)模型計(jì)算結(jié)果的驗(yàn)證
表5 Kalina 循環(huán)系統(tǒng)模型計(jì)算結(jié)果的驗(yàn)證
表4中狀態(tài)點(diǎn)1的溫度值與文獻(xiàn)結(jié)果有較大的誤差,主要是因?yàn)楸疚暮雎粤说責(zé)崴畯牡叵碌降厣线^(guò)程中的溫度損失,誤差由假設(shè)條件產(chǎn)生,因此并不影響對(duì)模型正確性的驗(yàn)證。通過(guò)表4和5可以看出,2個(gè)系統(tǒng)的計(jì)算結(jié)果與文獻(xiàn)值基本吻合,從而較好地驗(yàn)證了模型的正確性。
兩級(jí)閃蒸?Kalina循環(huán)聯(lián)合發(fā)電系統(tǒng)的熱力性能受不同因素的影響,包括地?zé)崴肟跍囟?、分離溫度、閃蒸溫度、冷凝溫度、氨質(zhì)量濃度、蒸發(fā)壓力等參數(shù)。其中,閃蒸溫度作為連接兩級(jí)閃蒸系統(tǒng)和Kalina循環(huán)系統(tǒng)的重要參數(shù),對(duì)聯(lián)合發(fā)電系統(tǒng)有較為明顯的影響。本文選用系統(tǒng)凈功率、循環(huán)效率、?效率等指標(biāo),分析不同情況下閃蒸溫度對(duì)系統(tǒng)熱力性能的影響,在地?zé)崴肟跍囟?50~300 ℃條件下,合理設(shè)定閃蒸溫度范圍為90~150 ℃,選用評(píng)價(jià)指標(biāo)列于表6中。
表6 系統(tǒng)所選用的評(píng)價(jià)指標(biāo)
取地?zé)崴肟跍囟葹?00 ℃,一級(jí)閃蒸后分離溫度取160 ℃,冷卻水溫度為5 ℃,透平入口壓力為 3 000 kPa,氨質(zhì)量濃度為 0.8。熱力性能指標(biāo)隨閃蒸溫度變化如圖3所示。圖3(a)為系統(tǒng)凈功率隨閃蒸溫度的變化規(guī)律曲線。從圖中可以看出,系統(tǒng)總凈功率和兩級(jí)閃蒸部分凈功率隨閃蒸溫度的提高先增大后減小,呈拋物線趨勢(shì),系統(tǒng)總凈功率最大值對(duì)應(yīng)的閃蒸溫度要高于兩級(jí)閃蒸系統(tǒng),而Kalina部分凈功率呈線性趨勢(shì)增大。原因是在此研究范圍內(nèi),閃蒸溫度的提高會(huì)增加氨蒸氣質(zhì)量流量,對(duì)Kalina循環(huán)是有利的。系統(tǒng)最大凈功率對(duì)應(yīng)閃蒸溫度為120 ℃,當(dāng)閃蒸溫度低于120 ℃時(shí),Kalina循環(huán)對(duì)系統(tǒng)影響更大;當(dāng)閃蒸溫度高于120 ℃時(shí),兩級(jí)閃蒸循環(huán)對(duì)系統(tǒng)影響更大。圖3(b)給出了系統(tǒng)循環(huán)效率、動(dòng)力回收效率和?效率隨閃蒸溫度的變化規(guī)律曲線。從圖中可以看出,系統(tǒng)循環(huán)效率隨閃蒸溫度的提高而增大,當(dāng)閃蒸溫度大于130 ℃時(shí),循環(huán)效率的增長(zhǎng)趨勢(shì)較緩慢,這是因?yàn)榘l(fā)生器換熱量和排放溫度在不斷提高;動(dòng)力回收效率和?效率隨閃蒸溫度的提高先增大后減小,呈拋物線趨勢(shì),最大值對(duì)應(yīng)的閃蒸溫度均為120 ℃,與凈功率變化規(guī)律一致。因此在設(shè)定條件下,認(rèn)為存在最佳閃蒸溫度為 120 ℃。
圖3 熱力性能指標(biāo)隨閃蒸溫度的變化曲線
圖4為系統(tǒng)余熱回收率和地?zé)崴欧艤囟入S閃蒸溫度的變化規(guī)律曲線。隨著閃蒸溫度的提高,地?zé)崴欧艤囟纫苍诓粩嗌撸到y(tǒng)余熱回收率在不斷下降。這是因?yàn)橛酂峄厥章逝c熱源排放溫度有關(guān);同時(shí),閃蒸溫度的提高導(dǎo)致了二級(jí)閃蒸階段單位質(zhì)量吸熱量的減小,使得系統(tǒng)的總吸熱量降低。
圖4 余熱回收率和回灌溫度隨閃蒸溫度的變化曲線
當(dāng)閃蒸溫度為120 ℃時(shí),取4種不同地?zé)崴肟跍囟?,?jì)算得到相應(yīng)的最佳分離溫度如表7所示。將4組參數(shù)分別做模擬計(jì)算,分析在不同地?zé)釡囟认?,閃蒸溫度對(duì)DFKCS系統(tǒng)熱力性能的影響。
表 7 不同熱源溫度對(duì)應(yīng)的最佳分離溫度 ℃
圖5為熱力性能指標(biāo)隨熱源溫度與閃蒸溫度的變化曲線。
圖5 熱力性能指標(biāo)隨熱源溫度與閃蒸溫度的變化曲線
圖5(a)表示的是不同熱源溫度下系統(tǒng)凈功率隨閃蒸溫度的變化規(guī)律。不同熱源溫度下,凈功率隨閃蒸溫度的變化規(guī)律基本一致,均為先增大后減小的拋物線趨勢(shì),但系統(tǒng)凈功率達(dá)到最大值時(shí)所對(duì)應(yīng)的閃蒸溫度,呈線性增大趨勢(shì),熱源溫度每提升50 ℃,對(duì)應(yīng)最佳閃蒸溫度提高10 ℃。在同一閃蒸溫度下,隨著熱源溫度的提高,系統(tǒng)凈功率增幅明顯提升。因此,從凈功率指標(biāo)來(lái)看,提升系統(tǒng)的熱源溫度,對(duì)系統(tǒng)發(fā)電是有利的。
不同熱源溫度下系統(tǒng)循環(huán)效率隨閃蒸溫度的變化規(guī)律如圖5(b)所示。可以看出,不同熱源溫度下,系統(tǒng)循環(huán)效率隨閃蒸溫度的提升而增大,均在接近臨界溫度時(shí)達(dá)到最大值,同時(shí),熱源溫度越高,系統(tǒng)循環(huán)效率越大。圖5(c)為不同熱源度下系統(tǒng)?效率隨閃蒸溫度的變化規(guī)律曲線,熱源溫度提高,?效率隨閃蒸溫度的變化趨勢(shì)越平緩,當(dāng)熱源溫度分別為150 ℃和300 ℃時(shí),隨著閃蒸溫度的提高,?效率增幅分別為15.66%和11.75%。
冷卻水溫度會(huì)隨著季節(jié)的變化而變化,對(duì)DFKCS系統(tǒng)工質(zhì)冷凝效果的影響非常顯著。圖6給出了不同冷卻水溫度下系統(tǒng)凈功率和循環(huán)效率的變化規(guī)律曲線??梢钥闯觯W蒸溫度提高,系統(tǒng)凈功率先增大后減小,循環(huán)效率則以不斷變緩的趨勢(shì)逐漸增大。隨著冷卻水溫度的提高,系統(tǒng)凈功率和循環(huán)效率都隨之減小,主要原因在于氨工質(zhì)冷凝溫度提高,造成透平的出口壓力增大,系統(tǒng)的凈輸出功減小,這說(shuō)明冷卻水溫度越高,對(duì)系統(tǒng)越不利。
圖6 熱力性能指標(biāo)隨冷卻水溫度與閃蒸溫度的變化曲線
不同冷卻水溫度下,系統(tǒng)凈功率對(duì)應(yīng)的最佳閃蒸溫度基本相同,在閃蒸溫度較低范圍內(nèi),系統(tǒng)凈功率和循環(huán)效率的增幅隨著冷卻水溫度的降低而明顯提升,可以看出Kalina循環(huán)在中低溫?zé)嵩礂l件下有較好的熱力性能。
DFKCS系統(tǒng)中Kalina循環(huán)部分的循環(huán)工質(zhì)為氨水溶液,工質(zhì)中氨的質(zhì)量濃度會(huì)直接影響到發(fā)生器蒸發(fā)過(guò)程氣液兩相的比例,從而影響氨蒸氣質(zhì)量流量和透平做功。圖7表示的是不同氨質(zhì)量濃度下系統(tǒng)凈功率和循環(huán)效率的變化規(guī)律曲線。同樣,當(dāng)閃蒸溫度提高時(shí),系統(tǒng)凈功率先增大后減小,而循環(huán)效率則以不斷變緩的趨勢(shì)逐漸增大。
圖7 熱力性能指標(biāo)隨氨質(zhì)量濃度與閃蒸溫度的變化曲線
圖7(a)中,氨質(zhì)量濃度越高,凈功率就不斷提高,同時(shí)隨閃蒸溫度的變化越平緩。當(dāng)氨質(zhì)量濃度為0.5、0.6、0.7和0.8時(shí),凈功率最大值時(shí)的閃蒸溫度分別為 129、126、123、120 ℃,呈線性減小趨勢(shì)。圖7(b)中,當(dāng)閃蒸溫度低于120 ℃時(shí),隨著氨質(zhì)量濃度提升,系統(tǒng)的循環(huán)效率不斷提高。當(dāng)閃蒸溫度高于120 ℃時(shí),系統(tǒng)循環(huán)效率受氨質(zhì)量濃度的影響較小,這說(shuō)明系統(tǒng)受多種因素的綜合影響,氨質(zhì)量濃度在一定條件下,并不是越高越好。
針對(duì)增強(qiáng)型地?zé)嵯到y(tǒng)產(chǎn)熱溫度高的特點(diǎn),本文提出了兩級(jí)閃蒸?Kalina循環(huán)聯(lián)合發(fā)電系統(tǒng),并通過(guò)模擬計(jì)算,有效驗(yàn)證了系統(tǒng)模型的正確性。得到結(jié)論如下:
1)分析了系統(tǒng)運(yùn)行中的閃蒸溫度對(duì)DFKCS系統(tǒng)的熱力性能影響,在合理設(shè)定其他參數(shù)的情況下,隨著閃蒸溫度的提高,系統(tǒng)凈功率、動(dòng)力回收效率和?效率呈先增大后減小的趨勢(shì),而系統(tǒng)循環(huán)效率呈緩慢上升趨勢(shì),余熱回收率降低,說(shuō)明系統(tǒng)存在最佳閃蒸溫度,本文中最佳閃蒸溫度為 120 ℃。
2)地?zé)崴肟跍囟让刻嵘?0 ℃,系統(tǒng)凈功率達(dá)到最大值時(shí)所對(duì)應(yīng)的閃蒸溫度提高10 ℃;同一閃蒸溫度下,地?zé)崴肟跍囟仍礁?,DFKCS系統(tǒng)發(fā)電功率越高,熱力性能越好。
3)冷卻水溫度提高時(shí),系統(tǒng)凈功率和循環(huán)效率隨之減小,說(shuō)明冷卻水溫度越高,對(duì)DFKCS系統(tǒng)越不利;氨質(zhì)量濃度提高時(shí),系統(tǒng)凈功率緩慢提高,凈功率最大值對(duì)應(yīng)的閃蒸溫度呈線性減小趨勢(shì),結(jié)合循環(huán)效率的變化,發(fā)現(xiàn)氨質(zhì)量濃度卻并非越高越好。因此在實(shí)際應(yīng)用中,要考慮冷卻水溫度、氨質(zhì)量濃度等參數(shù)的合理選擇,保證系統(tǒng)有更好地?zé)崃π阅堋?/p>