實現(xiàn)操作費用最小化的LNG接收站卸船工藝優(yōu)化

吳 明 朱祚良 孫東旭 何俊男 唐 凱 胡本源 田士章

1.遼寧石油化工大學(xué)石油天然氣工程學(xué)院 2.中國石油大學(xué)（華東）儲運與建筑工程學(xué)院 3.中石油江蘇液化天然氣有限公司

0 引言

LNG接收站最首要的過程是接收由LNG船運送至接收站的LNG[1]，稱為卸船工藝流程。卸船過程中的蒸發(fā)氣（Boil Off Gas，以下簡稱BOG）產(chǎn)生量在接收站運行過程中占有較大比重。因此，系統(tǒng)地對卸船工藝流程進行優(yōu)化對減少接收站運行成本具有現(xiàn)實意義。何淼[2]通過項目實際討論了LNG接收站卸船過程中影響B(tài)OG產(chǎn)生量的各種因素。賈士棟等[3]通過建立一維流動傳熱模型對LNG接收站卸料管線的預(yù)冷過程進行了模擬，并分析了不同因素對預(yù)冷過程的影響。Park等[4]考慮了BOG產(chǎn)生量對直接外輸和再冷凝兩種BOG處理工藝之間切換的影響，研究了使操作費用最小化的保冷循環(huán)LNG流量。Lee等[5]針對地上儲罐和地下儲罐混合卸船情況，提出了優(yōu)化的卸船流程。Srikanth等[6]開發(fā)了卸載操作的動態(tài)模型，實現(xiàn)了預(yù)冷和卸載操作策略之間的平衡。然而，針對操作變量對卸船工藝功耗的影響及其相關(guān)優(yōu)化的研究卻少有報道。

筆者基于卸船工藝操作費用的構(gòu)成，建立了卸船工藝流程動態(tài)仿真模型，進而構(gòu)建了以最小年總功耗為目標函數(shù)的卸船工藝優(yōu)化運行模型，并應(yīng)用優(yōu)化運行模型對計算實例進行經(jīng)濟效益對比分析。

圖1 LNG接收站卸船工藝流程圖

1 LNG卸船工藝過程

1.1 卸船工藝

LNG接收站卸船工藝流程如圖1所示。卸船工藝流程可分為3個階段：保冷循環(huán)階段、泄壓階段和卸船階段。以下分述之。

1）保冷循環(huán)階段。由于LNG本身低溫、易揮發(fā)的特性，一旦有外界的熱量流入[7]，LNG將不可避免地蒸發(fā)產(chǎn)生BOG[8]。因此當LNG接收站內(nèi)無船只進行卸船時，接收站需要持續(xù)對卸船總管進行保冷循環(huán)，以減少卸船階段BOG的產(chǎn)生[9]。具體操作為：LNG儲罐內(nèi)低壓泵出口分出一部分LNG通過保冷循環(huán)管線流入卸船總管進行保冷循環(huán)，循環(huán)后的LNG一小部分返回罐內(nèi)，大部分輸送至低壓外輸總管，以避免罐內(nèi)因流入過量高溫LNG而產(chǎn)生大量BOG[10]。

2）泄壓階段。由于LNG船內(nèi)壓力遠小于保冷循環(huán)時卸船總管內(nèi)壓力，因此在船舶進港卸船前需要降低卸船總管的壓力。在此過程中，關(guān)閉保冷循環(huán)進、出口閥，使得管道內(nèi)的LNG通過旁通閥流入儲罐，直至卸船總管中的壓力降至略小于碼頭卸料臂的壓力，以便于卸船階段的順利進行。

3）卸船階段。在對LNG進行全流量卸載前，需將卸料臂冷卻至正常操作溫度[11]。隨后打開卸料閥和卸料臂閥，將船內(nèi)LNG卸載至儲罐內(nèi)。在LNG船完成卸船操作后，保冷循環(huán)系統(tǒng)將重新啟動，直到下一艘LNG船到達接收站[12]。

1.2 操作變量

有兩個主要操作變量可用于調(diào)控LNG卸船工藝流程，分別是保冷循環(huán)LNG體積流量（Qrec）和保冷循環(huán)期間支路管線LNG流量（qbr）。

Qrec的作用是維持管線溫度介于－160～－158℃。由于循環(huán)冷卻，循環(huán)出口溫度介于－157～－155 ℃，比入口溫度（－160 ℃）稍高一些。根據(jù)接收站安全運行規(guī)定，循環(huán)進出口溫度差推薦值為3～5 ℃[13]。通過改變循環(huán)流量，可以調(diào)節(jié)溫度差值。當減少循環(huán)流量時，由于從環(huán)境傳熱，循環(huán)進出口溫度差增加，反之亦然。

qbr的作用是保冷支路管線。支路管線流量一般比總循環(huán)流量小得多。在冷卻支路管線后，LNG流入儲罐。可以通過調(diào)節(jié)支路管線流量控制冷卻工況。

1.3 操作費用分析

操作費用是接收站功耗的體現(xiàn)，卸船工藝流程的功耗分為泵功耗和壓縮機功耗兩部分。

在保冷循環(huán)階段，泵的作用是通過保冷循環(huán)管線將LNG輸送至卸船總管，從而起到保冷卸船管線的作用；而在泄壓階段和卸船階段，由于保冷循環(huán)系統(tǒng)關(guān)閉，泵停止向保冷循環(huán)管線輸送LNG。因此，泵僅在保冷循環(huán)階段工作，功耗與參與保冷循環(huán)的LNG流量有關(guān)。泵功率用下式計算[14]。

式中Lpump（Qrec）表示受Qrec影響的泵功率，W；Δp（Qrec）表示受Qrec影響的保冷管線進出口壓力損失，Pa；S表示安全系數(shù)；ηpump表示泵效率；ηmotor表示電動機效率。

在LNG接收站卸船工藝流程中，壓縮機的作用是將產(chǎn)生的BOG輸送至下游再冷凝器。在保冷循環(huán)階段，由于LNG經(jīng)保冷循環(huán)后溫度升高，流入儲罐的部分LNG必然會產(chǎn)生BOG；且由于閥門的節(jié)流現(xiàn)象[15]，LNG通過旁通閥后壓力下降，將產(chǎn)生額外的BOG。在泄壓階段，由于管道內(nèi)壓力逐漸降低和閥門的節(jié)流現(xiàn)象，管道內(nèi)會產(chǎn)生一定BOG，但產(chǎn)量較低。在卸船階段，LNG船卸載了大量的低溫LNG，此時由于泄壓過程中保冷循環(huán)LNG流量的減少，卸船總管的溫度必然升高，同時伴隨著各個通過閥的節(jié)流現(xiàn)象，管道中會產(chǎn)生大量的BOG。因此，壓縮機在整個卸船工藝流程持續(xù)工作，其功耗與BOG產(chǎn)生量有關(guān)。壓縮機功率用下式計算[14]。

式中Lcomp（Qrec,qbr）表示受Qrec和qbr影響的壓縮機功率，W；k表示比熱容比；QI（Qrec,qbr）表示受Qrec和qbr影響的BOG體積流量，m3/h；pI表示壓縮機入口壓力，Pa；pO表示壓縮機出口壓力，Pa；ηB表示機械效率。

從工藝功耗的角度出發(fā)，引入卸船周期的概念。一個卸船周期可分為保冷循環(huán)和卸船兩個階段。由于泄壓階段停止了LNG保冷循環(huán)，泵不產(chǎn)生功耗，同時因泄壓時間較短，壓縮機的功耗可忽略，所以卸船周期未包含泄壓階段。

因此，一個卸船周期內(nèi)的功耗可分為3部分：用于輸送保冷循環(huán)LNG的低壓泵功耗、保冷循環(huán)階段，用于將BOG排出儲罐的壓縮機功耗和卸船階段，用于將BOG排出儲罐的壓縮機功耗?？梢耘袛嗟蛪罕霉呐c卸船階段壓縮機的功耗是相互對立的：當保冷循環(huán)LNG流量增大時，低壓泵功耗增加，管道保冷效果改善，管壁溫度降低，卸船時產(chǎn)生的BOG氣體量隨之減少，BOG壓縮機功耗減小；反之，當保冷循環(huán)LNG流量減小時，雖然泵提供的動力功耗減小，但由于管道溫度較高，使卸船時產(chǎn)生的BOG量增大，即增加了壓縮機功耗。因此可以說泵貢獻了卸船工藝的動力功耗，而壓縮機貢獻了卸船工藝因熱力問題而產(chǎn)生的功耗，兩者聯(lián)系緊密。泵與壓縮機的功耗之和即為卸船工藝總功耗，一個卸船周期內(nèi)總功耗計算式為：

式中WT（Qrec,qbr）表示受Qrec和qbr影響的卸船周期總功耗，kW·h；Lreccomp（Qrec,qbr）表示受Qrec和qbr影響的保冷循環(huán)階段壓縮機功率，W；Lunlcomp（Qrec,qbr）表示受Qrec和qbr影響的卸船階段壓縮機功率，W；T表示卸船周期，h；Tunl表示卸船時間，h。

2 卸船工藝流程控制方程及動態(tài)模擬

為了準確模擬實際的卸船過程，需要解決3方面問題：LNG卸船過程中的流動、傳熱與相變規(guī)律。運動方程描述了卸船系統(tǒng)流量與進出口壓差、管道阻力之間的關(guān)系；傳熱方程描述了管道內(nèi)LNG溫度場分布；狀態(tài)方程解答了流體相態(tài)與溫度、壓力和密度之間的變化規(guī)律。

2.1 運動方程

LNG在管道內(nèi)流動時由于溫度升高和壓力降低，可能蒸發(fā)為BOG氣體，即存在多相流動問題?；诰嗔鲃幽Ｐ偷贸龅腂eggs-Brill方程適用于水平管、垂直管和傾斜管的多相流問題求解[16]，壓降梯度方程如式（4）所示。

式（4）中等號右側(cè)分子第一項與第二項分別表示消耗于位差的壓力損失和摩擦阻力引起的壓力損失，分母表示動能變化引起的壓力損失[17]。

2.2 傳熱方程

導(dǎo)熱微分方程描述了管道溫度場的分布規(guī)律，對于管道溫度場計算適合使用柱坐標系形式的導(dǎo)熱微分方程。當忽略管道周向溫度變化、無內(nèi)熱源且導(dǎo)熱系數(shù)為常數(shù)時的導(dǎo)熱微分方程如式（5）所示[18-19]。

式中ρ表示導(dǎo)熱介質(zhì)密度，kg/m3；c表示導(dǎo)熱介質(zhì)比熱容，J/（ kg·K）； t表示溫度，℃；τ表示時間，s；K表示介質(zhì)的傳熱系數(shù)；r表示距離管道中心軸線的徑向距離，m；z表示管道軸向距離，m。

圖2 LNG接收站卸船工藝流程動態(tài)模型圖

2.3 狀態(tài)方程

Stryjek和Vera[20]于1986年對PR方程進行了修正，提出了新溫度函數(shù)形式的PR方程，稱為PRSV方程[21]。此方程可用于極性和非極性物質(zhì)，采用普通的混合規(guī)則就可以精確地預(yù)測氣液平衡數(shù)據(jù)[22]。

PRSV方程的形式為[23-24]：

式中p表示壓力，MPa；R表示氣體常數(shù)，8.314 J/（mol·K）；t表示溫度，K；V表示摩爾體積，m3/mol；a、b表示狀態(tài)方程參數(shù)。

式中tc表示臨界溫度，K；pc表示臨界壓力，Pa；α（t）表示溫度函數(shù)；k、k0表示狀態(tài)參數(shù)；tr表示對比溫度；k1表示每個純組分的可調(diào)整參數(shù)；ω表示偏心因子。

2.4 動態(tài)模擬

LNG卸船過程每一階段都有各自的邊界條件，下一階段的初始條件正是由上一階段建立起的溫度場、壓力場和物性條件構(gòu)成的。研究采用HYSYS V8.4為工具，使用Beggs-Brill方程完成多相流問題求解，選用PRSV方程作為物性計算方法，建立LNG卸船工藝流程動態(tài)模型并求解。模型如圖2所示?；诮邮照镜膶嶋H運行工況輸入相關(guān)的環(huán)境變量參數(shù)值以及LNG物性參數(shù)，并通過設(shè)置模型中閥門的開關(guān)狀態(tài)依次模擬保冷循環(huán)、泄壓和卸船3個過程。在模擬卸船過程時，由于瞬間產(chǎn)生大量的BOG與LNG混合，在管道內(nèi)形成兩相流動，可能會形成不穩(wěn)定的模擬結(jié)果。若模擬結(jié)果不穩(wěn)定，應(yīng)在調(diào)整保冷循環(huán)流量和支路管線流量后重新模擬。在每一階段都得到穩(wěn)定的模擬結(jié)果后，即可計算泵、壓縮機和卸船周期總功耗。動態(tài)模擬流程如圖3所示。

3 卸船工藝優(yōu)化運行模型及應(yīng)用

3.1 卸船工藝優(yōu)化運行模型

結(jié)合式（1）～（6），卸船工藝優(yōu)化運行模型由目標函數(shù)、控制方程、邊界條件和初始條件構(gòu)成。

目標函數(shù)：

式中WA（Qrec,qbr）表示受Qrec和qbr影響的年總功耗，kW·h；A表示年卸船周期數(shù)。

控制方程：選用式（4）～（6）為優(yōu)化模型的控制方程。

邊界條件：

初始條件：保冷循環(huán)達到穩(wěn)態(tài)時各變量值。

模型中有兩個變量：Qrec和qbr。優(yōu)化模型的目標函數(shù)、控制方程是非線性的，研究用如下方法求解優(yōu)化模型：將保冷循環(huán)LNG流量作為外層循環(huán)，將支路管線流量作為內(nèi)層循環(huán)。根據(jù)目標函數(shù)、控制方程及約束條件計算得到可行運行方案的集合，集合中最小值即為最優(yōu)解，計算流程如圖4所示。

圖4 求最優(yōu)解流程圖

3.2 實例計算與經(jīng)濟效益分析

基于江蘇如東LNG接收站實際運行參數(shù)進行實例分析。計算實例中天然氣組分的物性參數(shù)如表1所示。接收站現(xiàn)有卸船工藝的運行參數(shù)不隨環(huán)境參數(shù)變化，保冷循環(huán)流量（Qrec）和支路管線流量（qbr）分別為262 m3/h和12 m3/h。根據(jù)接收站下游用氣量波動與環(huán)境參數(shù)的變化，可將卸船工況分為11月—次年3月的用氣高峰期與4—10月的用氣低谷期，涉及到的環(huán)境變量參數(shù)如表2所示。使用如圖2所示動態(tài)仿真模型對現(xiàn)有卸船工藝進行動態(tài)模擬，得到接收站現(xiàn)有卸船工藝的年總功耗為2 010 MW·h。

表1 天然氣組分的物性參數(shù)表

表2 環(huán)境變量及相關(guān)工藝參數(shù)表

圖5為接收站用氣高峰期與低谷期泄壓階段與部分卸船階段的LNG與BOG流量曲線，其中0～10 min為泄壓階段，無LNG進罐且無BOG產(chǎn)生。卸船操作在第10 min開始，在卸船開始的瞬間，BOG大量產(chǎn)生并流入儲罐，隨后BOG流量逐漸下降，直至趨于平穩(wěn)。這是因為在卸船開始的瞬間，由于管道與所卸載LNG的溫差較大，同時伴隨著閥門的節(jié)流效應(yīng)，BOG流量會瞬間達到峰值；隨著閥門完全開啟且由于卸載低溫LNG導(dǎo)致管道溫度下降，BOG流量會逐漸減小直至趨于平穩(wěn)。LNG流量在開始卸船后較短時間內(nèi)增至最大卸船流量，并持續(xù)至卸船結(jié)束。

圖5 用氣高峰期與低谷期接收站泄壓階段和部分卸船階段的LNG與BOG流量曲線圖

圖6 用氣高峰期接收站年卸船總功耗隨保冷循環(huán)LNG流量和支路管線流量變化趨勢圖

應(yīng)用卸船工藝優(yōu)化運行模型對接收站卸船工藝進行優(yōu)化，同時進行非穩(wěn)態(tài)傳熱過程分析，得到接收站卸船期間最小年功耗及對應(yīng)運行工況。圖6、7展示了用氣高峰期與低谷期接收站卸船總功耗隨操作變量的變化情況。從圖6、7可知，隨著Qrec的增加，用氣高峰期卸船總功耗隨之先減小后增加，用氣低谷期卸船總功耗逐漸增加；在這兩種工況下，qbr對總功耗的影響相對較小。經(jīng)優(yōu)化后的用氣高峰期操作變量為Qrec=178.5 m3/h，qbr=9 m3/h，用氣低谷期為Qrec=156.7 m3/h，qbr=9 m3/h，此時接收站年卸船總功耗為1 749 MW·h，相對于優(yōu)化前的2 010 MW·h節(jié)能13.0%。

圖7 用氣低谷期接收站年卸船總功耗隨保冷循環(huán)LNG流量和支路管線流量變化趨勢圖

4 結(jié)論

1）通過對LNG接收站卸船工藝非穩(wěn)態(tài)傳熱過程的分析，建立了動態(tài)仿真模型，對卸船期間非穩(wěn)態(tài)傳熱過程進行了動態(tài)模擬，建立了卸船工藝優(yōu)化運行模型。與傳統(tǒng)運行模式相比，卸船工藝優(yōu)化運行模型的提出使得LNG接收站在卸船工藝操作期間能夠通過對操作變量的調(diào)整，實現(xiàn)最小化接收站操作費用的目的。

2）通過將研究成果應(yīng)用于江蘇如東LNG接收站，在使用卸船工藝優(yōu)化運行模型進行優(yōu)化后，通過調(diào)節(jié)保冷循環(huán)LNG流量和支路管線流量可使年卸船總功耗降低13.0%。