成慶林 王 雪 孟 嵐 孫 巍 劉鶴皋

1. 東北石油大學(xué)提高油氣采收率教育部重點(diǎn)實(shí)驗(yàn)室, 黑龍江 大慶 163000;

2. 大慶油田工程有限責(zé)任公司, 黑龍江 大慶 163000

0 前言

中國(guó)屬于能源消耗大國(guó),一次能源占有率低,油氣消費(fèi)對(duì)外依存高,現(xiàn)有的能源結(jié)構(gòu)對(duì)生態(tài)環(huán)境造成過(guò)大壓力,迫切需要能源轉(zhuǎn)型及能耗優(yōu)化[1]。未來(lái)30年天然氣和非化石能源的地位將會(huì)顯著提升,中國(guó)能源消費(fèi)結(jié)構(gòu)將發(fā)生質(zhì)的變化[2]。由于天然氣資源充足而且潔凈有效,是可獲取、可接受、可消費(fèi)的優(yōu)質(zhì)能源,所以促進(jìn)天然氣發(fā)展、進(jìn)行天然氣改革是促進(jìn)中國(guó)燃料生產(chǎn)與消費(fèi)革命的重大措施[3-4]。

天然氣在凈化儲(chǔ)存和運(yùn)輸過(guò)程中,對(duì)天然氣含水量的要求比較嚴(yán)格。如果天然氣中有水分存在,在輸送過(guò)程中,一定程度上會(huì)降低天然氣的熱值和管道輸送能力,還會(huì)縮短管道使用年限,這與生產(chǎn)經(jīng)濟(jì)與安全息息相關(guān),故天然氣脫水工藝是凈化處理過(guò)程中的重要環(huán)節(jié)。三甘醇(TEG)脫水工藝是有效控制管道天然氣水露點(diǎn)的常用脫水方法,Mohamadbeigy K H[5]研究了甘醇循環(huán)量、汽提氣量與吸收塔塔板數(shù)對(duì)TEG脫水工藝的影響;李天斌[6]使用Aspen HYSYS軟件對(duì)TEG脫水及再生系統(tǒng)進(jìn)行模擬計(jì)算,優(yōu)化工藝參數(shù),確定各變量之間的相互關(guān)系,從而找到最佳操作點(diǎn);蔣洪等人[7]采用二次回歸正交實(shí)驗(yàn)設(shè)計(jì)能耗優(yōu)化方法,有效優(yōu)化了三甘醇脫水裝置的能耗;周軍等人[8]將BP神經(jīng)網(wǎng)絡(luò)與GA遺傳算法結(jié)合,得出脫水工藝過(guò)程優(yōu)化操作參數(shù),單位能耗下降18.13%。但在以往的優(yōu)化研究中，側(cè)重各操作參數(shù)的敏感度分析以影響三甘醇脫水效果,聚焦于脫水裝置的能耗優(yōu)化而忽略了整個(gè)系統(tǒng)的碳排放水平。

本文利用Aspen HYSYS軟件,根據(jù)西北某氣田天然氣TEG脫水工藝流程建立數(shù)學(xué)模型，根據(jù)實(shí)際生產(chǎn)情況,選取影響生產(chǎn)效率的實(shí)際參數(shù),對(duì)重沸器負(fù)荷進(jìn)行敏感度分析;引入碳排放指標(biāo),建立碳排放優(yōu)化模型;采用Aspen HYSYS軟件優(yōu)化器進(jìn)行優(yōu)化,求解出滿足處理工藝要求的最優(yōu)結(jié)果,為油氣田生產(chǎn)TEG脫水處理工藝的繼續(xù)優(yōu)化提供參考。

1 TEG脫水系統(tǒng)模型建立

1.1 TEG脫水系統(tǒng)工藝流程簡(jiǎn)介

TEG脫水系統(tǒng)主要包括天然氣常壓脫水和低壓高溫TEG富溶液(濃縮)再生兩部分。在三甘醇脫水裝置中,三甘醇溶液被用作此設(shè)計(jì)工藝流程脫水劑,作用是去除凈化天然氣中大部分的水分。經(jīng)三甘醇吸收器(工作壓力下水露點(diǎn)<-5 ℃)脫水的干燥凈化天然氣作為商業(yè)天然氣出口[9]。TEG脫水系統(tǒng)工藝流程見(jiàn)圖1。

圖1 TEG脫水系統(tǒng)工藝流程圖Fig.1 Process flow diagram of TEG dehydration system

原料天然氣首先通過(guò)分離裝置進(jìn)入TEG吸收塔底部,天然氣中的水分在TEG吸收塔頂部通過(guò)與TEG溶液逆流接觸后被脫除,干天然氣途經(jīng)TEG吸收塔的頂部到達(dá)分離器,經(jīng)分離裝置分離的產(chǎn)品氣被輸送出去。富TEG溶液從吸收塔出來(lái)后與再生塔頂部之間進(jìn)行熱量交換,然后進(jìn)入TEG富液換熱器,再進(jìn)入閃蒸罐,溶解在閃蒸罐中的碳?xì)浠衔飼?huì)很快出現(xiàn),交換之后依次進(jìn)入TEG富液換熱器和閃蒸罐。TEG脫水處理工藝再生過(guò)程中,產(chǎn)生的廢氣被廢氣焚燒爐焚燒,再生的TEG貧液從再生系統(tǒng)出來(lái)后依次進(jìn)入重沸器和緩沖罐。進(jìn)行熱交換后,通過(guò)將貧TEG溶液送入TEG循環(huán)泵,再送入換熱器,最后送至TEG吸收塔頂部,實(shí)現(xiàn)TEG的循環(huán)利用。

1.2 設(shè)計(jì)物料參數(shù)及組分

西北某氣田的天然氣水露點(diǎn)要求<-5 ℃,濕天然氣的實(shí)際參數(shù)及產(chǎn)品氣的設(shè)計(jì)物料參數(shù)見(jiàn)表1,物料組分見(jiàn)表2。

表2 設(shè)計(jì)物料摩爾組分表Tab.2 Design material molar components

1.3 Aspen HYSYS模擬軟件建立

Aspen HYSYS模擬軟件是一種常見(jiàn)的大型油氣加工模擬軟件,包含動(dòng)態(tài)和穩(wěn)態(tài)兩個(gè)部分。Aspen HYSYS模擬軟件可以將穩(wěn)態(tài)系統(tǒng)、控制理論、動(dòng)態(tài)化工和熱力學(xué)模型、動(dòng)態(tài)數(shù)據(jù)處理等有機(jī)地結(jié)合起來(lái),通過(guò)求解巨型常微分方程組進(jìn)行動(dòng)態(tài)仿真模擬[10-12]。本文選用Aspen HYSYS模擬軟件進(jìn)行工藝流程建模、參數(shù)提取、參數(shù)分析和優(yōu)化過(guò)程。

Aspen HYSYS模擬軟件中,Peng-Robinson狀態(tài)方程(即P-R狀態(tài)方程)具有很多優(yōu)點(diǎn),比如精確度較高、適用范圍廣等,尤其是在TEG脫水工藝流程的氣液平衡模型中表現(xiàn)出良好的收斂性和較高的精度[13]。為使脫水工藝流程快速收斂,同時(shí)減少誤差,本文選P-R狀態(tài)方程[14-15]，該軟件初步建立的TEG脫水工藝流程模擬見(jiàn)圖2。

圖2 TEG脫水工藝流程模擬圖Fig.2 Simulation diagram of TEG dehydration process flow

2 TEG系統(tǒng)優(yōu)化方法建立

碳排放的衡量離不開(kāi)能耗,能耗是導(dǎo)致碳排放的必要因素之一。根據(jù)文獻(xiàn)調(diào)研和實(shí)際生產(chǎn)情況,當(dāng)富甘醇進(jìn)入再生塔時(shí),入塔時(shí)的溫度對(duì)產(chǎn)品氣的質(zhì)量有相對(duì)較小的影響。在工藝流程中,汽提氣的用量也會(huì)影響其熱值和汽提氣加熱器的負(fù)荷。在TEG脫水處理工藝流程中,如果入塔時(shí)貧甘醇流量、重沸器溫度及汽提氣的用量被嚴(yán)格控制住,可以減少蒸汽的消耗量,從而大大減少TEG脫水裝置的能量損耗[16-17],導(dǎo)致碳排放量增加。

2.1 因素敏感度分析

通過(guò)敏感性分析,可以確定項(xiàng)目的敏感性因素,并確定這些因素在評(píng)估指標(biāo)中改變時(shí)的影響程度[18]。利用HYSYS軟件分別對(duì)脫水系統(tǒng)內(nèi)TEG循環(huán)量、汽提氣用量、重沸器溫度對(duì)重沸器負(fù)荷的影響進(jìn)行敏感度分析,根據(jù)重沸器負(fù)荷變化得出關(guān)鍵因素。重沸器負(fù)荷變化見(jiàn)圖3。

圖3 重沸器負(fù)荷變化圖Fig.3 Reboiler load variation

通過(guò)敏感性分析,可以找出其關(guān)鍵敏感因素。由圖3可知:當(dāng)改變汽提氣用量時(shí),重沸器負(fù)荷變化幅度是最小的;當(dāng)改變重沸器溫度時(shí),重沸器負(fù)荷變化幅度較為明顯;當(dāng)改變TEG循環(huán)量時(shí),重沸器負(fù)荷變化幅度最為明顯。

2.2 TEG系統(tǒng)碳排放源識(shí)別

“雙碳”目標(biāo)的提出,將引領(lǐng)石油化工行業(yè)進(jìn)入前所未有的時(shí)代,全球石油化工行業(yè)競(jìng)爭(zhēng)格局也將發(fā)生重大轉(zhuǎn)變[19]?！吨袊?guó)石油天然氣生產(chǎn)企業(yè)溫室氣體排放核算方法與報(bào)告指南(試行)—2015》主要包含燃料燃燒、逸散排放、工藝放空的直接排放源,以及用電消耗的間接碳排放源。此次設(shè)計(jì)對(duì)TEG脫水系統(tǒng)進(jìn)行碳排放源識(shí)別分析。在此次設(shè)計(jì)工藝流程中,TEG循環(huán)泵設(shè)備為主要能耗設(shè)備,屬于間接碳排放源,而重沸器設(shè)備為耗熱設(shè)備,再生塔設(shè)備也為耗熱設(shè)備,重沸器設(shè)備與再生塔設(shè)備都屬于燃燒源。此外,由于工藝流程中的緩沖罐和閥門都會(huì)產(chǎn)生少量的逸散碳排放量,因此廢氣焚燒爐中存在大量放空污染源。但由于選擇優(yōu)化的參數(shù)并不涉及緩沖罐、閥門和廢氣焚燒爐,因此TEG系統(tǒng)的碳排放量?jī)?yōu)化邊界并不涉及逸散排放源和工藝放空源。燃料燃燒源碳排放核算見(jiàn)式(1)。

(1)

間接排放源碳排放核算見(jiàn)式(2)。

E=AD2×EF

(2)

西北某氣田TEG脫水系統(tǒng)的熱能由天然氣燃燒供給,天然氣低位熱值為38 931 kJ/m3,含碳量為0.000 596 tCO2/m3,碳氧化率為0.99,西北區(qū)域電網(wǎng)平均CO2排放因子為0.667 1 tCO2/(MW·h)。

2.3 TEG系統(tǒng)碳排放量?jī)?yōu)化模型建立

TEG脫水工藝流程中,TEG循環(huán)泵、重沸器和再生塔為主要耗能裝置,碳排放節(jié)點(diǎn)為TEG循環(huán)泵的電耗，以及重沸器和汽提氣的熱耗。根據(jù)關(guān)鍵參數(shù)的控制選取,建立該系統(tǒng)碳排放量?jī)?yōu)化模型見(jiàn)式(3)～(5)。

(3)

s.t.l≤g(x1,x2,x3)≤u

(4)

z≤-5

(5)

2.4 Aspen HYSYS優(yōu)化方法選用

TEG循環(huán)量、濕天然氣進(jìn)塔時(shí)溫度、汽提氣的體積流量和重沸器溫度等參數(shù)對(duì)天然氣脫水效果和重沸器碳排放量影響相對(duì)較大[3,20]。由于濕天然氣的溫度主要由上游流程決定,該系統(tǒng)無(wú)法進(jìn)行調(diào)節(jié);此外,Optimizer優(yōu)化器不能將再生塔內(nèi)部的約束條件設(shè)置為自變量,所以針對(duì)脫水流程可優(yōu)化的變量為TEG循環(huán)量和汽提氣的體積流量。令TEG循環(huán)量x1,汽提氣用量x2為自變量,以單位用能E為目標(biāo)函數(shù),建立優(yōu)化模型[21],其約束方程設(shè)定見(jiàn)表3。

表3 Aspen HYSYS優(yōu)化約束條件表Tab.3 Aspen HYSYS optimization constraints

Aspen HYSYS軟件自帶優(yōu)化方法包括BOX、SQP、MIX、Fletcher Reeves、Quasi-Newton,見(jiàn)表4。

表4 Aspen HYSYS軟件自帶優(yōu)化方法表Tab.4 Aspen HYSYS built-in optimizer optional optimization method

根據(jù)TEG系統(tǒng)建立的碳排放量?jī)?yōu)化模型,結(jié)合Aspen HYSYS軟件自帶優(yōu)化方法的特點(diǎn),選擇BOX、SQP、MIX方法進(jìn)行優(yōu)化求解。優(yōu)化流程見(jiàn)圖4。

3 優(yōu)化結(jié)果分析

模擬結(jié)果見(jiàn)表5。對(duì)上述影響碳排放量的三個(gè)重要參數(shù)TEG循環(huán)量、汽提氣用量和重沸器溫度進(jìn)行分析后認(rèn)為:當(dāng)重沸器溫度降低時(shí),貧甘醇溶液的含量也相應(yīng)降低,同時(shí)產(chǎn)品氣的質(zhì)量也隨著貧甘醇溶液含量的降低而相應(yīng)下降,同時(shí)重沸器的加熱負(fù)荷和重沸器消耗的蒸汽量也隨之降低。與貧甘醇流量和重沸器溫度這兩個(gè)參數(shù)相比較,汽提氣用量的變化對(duì)重沸器蒸汽消耗量的影響相對(duì)較小。

圖4 HYSYS優(yōu)化流程圖Fig.4 Aspen HYSYS optimization flow chart

表5 Aspen HYSYS優(yōu)化前后結(jié)果對(duì)比表Tab.5 Comparison of results before and after Aspen HYSYS optimization

模擬結(jié)果表明,由于經(jīng)過(guò)優(yōu)化,TEG再生重沸器的加熱負(fù)荷明顯降低,從而導(dǎo)致了整體的單位綜合能耗顯著下降,節(jié)電效益突出,在實(shí)際工業(yè)生產(chǎn)過(guò)程中值得推廣應(yīng)用。通過(guò)對(duì)TEG脫水工藝參數(shù)和工藝過(guò)程優(yōu)化的研究,確定了TEG脫水工藝參數(shù),由于Aspen HYSYS優(yōu)化模型為理想狀態(tài)下的最優(yōu)結(jié)果,故在實(shí)際現(xiàn)場(chǎng)中選取SQP模型的模擬優(yōu)化結(jié)果作為操作參考。當(dāng)汽提氣用量、入塔貧甘醇流量這兩個(gè)主要參數(shù)的數(shù)值均滿足設(shè)計(jì)要求且重沸器的溫度設(shè)定為188 ℃時(shí),利用Aspen HYSYS軟件模擬三甘醇脫水工藝流程,通過(guò)優(yōu)化前后結(jié)果對(duì)比可知,優(yōu)化后能耗降低了約6.7%。

通過(guò)TEG系統(tǒng)優(yōu)化前后碳排放核算對(duì)比,以SQP模型下的結(jié)果為參考,該系統(tǒng)內(nèi)燃料燃燒源碳排放可達(dá)到0.050 43 tCO2/d,間接排放源排放可達(dá)到0.005 066 tCO2/d,總碳排放量可達(dá)到0.055 496 tCO2/d,優(yōu)化前的系統(tǒng)原始碳排放量為0.057 567 565 tCO2/d,優(yōu)化后的脫水系統(tǒng)可減少碳排放量0.002 07 tCO2/d,每年可減少0.756 0 tCO2的排放量。該優(yōu)化研究對(duì)“雙碳”的實(shí)施有著一定的經(jīng)濟(jì)參考價(jià)值。

4 建議

1)在TEG脫水工藝模擬過(guò)程中,將TEG溶液的循環(huán)量提高,其它參數(shù)保持恒定不變時(shí),露點(diǎn)溫降會(huì)變大;但是當(dāng)TEG溶液循環(huán)量提高到一定程度時(shí),重沸器設(shè)備的負(fù)荷將會(huì)增加,同時(shí)也會(huì)增加能耗。在實(shí)際生產(chǎn)過(guò)程中,通過(guò)合理的工藝流程優(yōu)化,選擇合適的工藝參數(shù),不僅可以提高生產(chǎn)效率,還可以達(dá)到節(jié)能環(huán)保減排的效果。

2)在使用Aspen HYSYS數(shù)值模擬軟件進(jìn)行工藝參數(shù)敏感性分析時(shí),結(jié)果顯示TEG循環(huán)量對(duì)重沸器熱負(fù)荷敏感度影響最大,其次是重沸器溫度,影響較小的是汽提氣用量。

3)基于碳核算的TEG系統(tǒng)優(yōu)化研究,建立碳排放模型,找出最佳工藝參數(shù),對(duì)企業(yè)碳排放核算清單的建立具有一定的參考價(jià)值。對(duì)于油氣生產(chǎn)系統(tǒng)的優(yōu)化,碳排放是未來(lái)參考的重要指標(biāo)之一。