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

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

0 引言

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

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

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

1 LNG卸船工藝過程

1.1 卸船工藝

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

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

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

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

1.2 操作變量

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

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

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

1.3 操作費(fèi)用分析

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

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

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

在LNG接收站卸船工藝流程中，壓縮機(jī)的作用是將產(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，此時(shí)由于泄壓過程中保冷循環(huán)LNG流量的減少，卸船總管的溫度必然升高，同時(shí)伴隨著各個(gè)通過閥的節(jié)流現(xiàn)象，管道中會產(chǎn)生大量的BOG。因此，壓縮機(jī)在整個(gè)卸船工藝流程持續(xù)工作，其功耗與BOG產(chǎn)生量有關(guān)。壓縮機(jī)功率用下式計(jì)算[14]。

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

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

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

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

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

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

2.1 運(yùn)動方程

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

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

2.2 傳熱方程

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

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

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

2.3 狀態(tài)方程

Stryjek和Vera[20]于1986年對PR方程進(jìn)行了修正，提出了新溫度函數(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表示每個(gè)純組分的可調(diào)整參數(shù)；ω表示偏心因子。

2.4 動態(tài)模擬

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

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

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

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

目標(biāo)函數(shù)：

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

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

邊界條件：

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

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

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

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

基于江蘇如東LNG接收站實(shí)際運(yùn)行參數(shù)進(jìn)行實(shí)例分析。計(jì)算實(shí)例中天然氣組分的物性參數(shù)如表1所示。接收站現(xiàn)有卸船工藝的運(yù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)有卸船工藝進(jìn)行動態(tài)模擬，得到接收站現(xiàn)有卸船工藝的年總功耗為2 010 MW·h。

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

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

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

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

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

應(yīng)用卸船工藝優(yōu)化運(yùn)行模型對接收站卸船工藝進(jìn)行優(yōu)化，同時(shí)進(jìn)行非穩(wěn)態(tài)傳熱過程分析，得到接收站卸船期間最小年功耗及對應(yīng)運(yùn)行工況。圖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，此時(shí)接收站年卸船總功耗為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)傳熱過程進(jìn)行了動態(tài)模擬，建立了卸船工藝優(yōu)化運(yùn)行模型。與傳統(tǒng)運(yùn)行模式相比，卸船工藝優(yōu)化運(yùn)行模型的提出使得LNG接收站在卸船工藝操作期間能夠通過對操作變量的調(diào)整，實(shí)現(xiàn)最小化接收站操作費(fèi)用的目的。

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