紀馨，陳雙雙，宋陽，林文勝

（上海交通大學制冷與低溫工程研究所，上海200240）

引 言

中間流體氣化器 （intermediate fluid vaporizer，IFV）是 FSRU［1］（floating storage re－gasification unit）再氣化模塊中經常使用的一種氣化器。現(xiàn)有的IFV［2］通常由3個單元組成：蒸發(fā)器、冷凝器和調溫器，一般采用丙烷為中間流體。

丙烷的相變換熱流程有效利用了丙烷的相變潛熱，減少了丙烷的儲存量和循環(huán)量。但由于氣態(tài)丙烷的比容遠大于液態(tài)，因此，雖然這種流程的丙烷循環(huán)量較小，但換熱器的整體體積仍然較大。根據FSRU特點，其系泊在海上，船體空間有限，為了探討減小換熱器體積的可能性，使系統(tǒng)盡可能安全，提出了過冷丙烷換熱流程［3］，其換熱流程如圖1所示。

圖1 過冷中間流體換熱流程Fig.1 Heat transfer process of super cooled liquid

在E－100換熱器中，丙烷提供足夠的熱量，LNG被液態(tài)丙烷加熱至常溫后輸出。丙烷與LNG換熱后溫度降低，隨后經過泵P－101輸送至E－101換熱器，被海水加熱至常溫后，回到E－100，完成丙烷的循環(huán)。丙烷在整個循環(huán)中均在一定壓力下工作，從而維持在過冷的液體狀態(tài)。

以過冷中間流體換熱器的理論計算為基礎，對該換熱過程進行實驗研究，研究換熱流程中流體的換熱特性，及結合實際FSRU運行背景討論對該換熱過程的影響因素?？紤]到LNG－IFV在海上工作，海面晃動會造成設備傾斜，本文進行了靜態(tài)條件下傾斜角度對換熱過程影響的實驗研究；研究了海水質量流量對換熱過程的影響；考慮到海水溫度隨地理位置和季節(jié)變化的不同，研究了海水入口溫度對換熱過程的影響。

1 實驗裝置

1.1 實驗系統(tǒng)

對丙烷作為制冷劑的換熱特性已經開展了較多研究［4－10］，雖然這些研究基本上都是考察蒸發(fā)或冷凝的相變傳熱過程的，很少涉及到無相變傳熱特性。參考現(xiàn)有文獻，設計了本文過冷中間流體換熱過程實驗裝置。

從安全角度出發(fā)，本實驗系統(tǒng)中用液氮代替LNG作為冷源，實驗系統(tǒng)如圖2所示，丙烷由丙烷泵8從儲罐9中抽出后流經流量計5后進入換熱段6，吸收水的熱量，由于其壓力較高，丙烷吸收熱量后不會發(fā)生蒸發(fā)相變過程，其仍然為液態(tài)，換熱后溫度較高的丙烷流入換熱段4與液氮進行換熱，丙烷熱量被液氮吸收，其溫度降低，經流量調節(jié)閥后流入丙烷儲罐9中完成其循環(huán)回路。

水經過冷水機組15進行溫度調控，經過流量計10后流入換熱段6與丙烷換熱，將其熱量傳遞給丙烷后流入冷水機組，完成其循環(huán)回路。

液氮由儲罐1經過閥2調節(jié)流量后進入換熱段4，吸收丙烷的熱量，降低丙烷溫度，換熱后直接排空。

結合FSRU海上運行背景，必須考慮到海平面的晃動對換熱的影響，故研究在靜態(tài)條件下，不同傾斜角度0°、5°、10°下的換熱特性。

圖2 實驗裝置Fig.2 Diagram of experimental set up

1.2 實驗原理

過冷丙烷換熱流程實驗中，將重點研究丙烷／水換熱段換熱特性，丙烷／水換熱段的示意圖如圖3所示。

圖3 換熱段示意圖Fig.3 Diagram of heat transfer process

丙烷與水為逆流換熱，前后端分別留有430 mm作為充分混合流動段，換熱測試段每段為228 mm，則每一段的換熱量由式 （1）和式 （2）求解

第i段換熱段的水與內壁面的局部傳熱系數的求解按照式 （3）～式 （5）求解

整個換熱段內水與內壁面的平均傳熱系數hw，avg由式 （6）求解所得

對于丙烷側傳熱系數的求解原理與求解水側傳熱系數一致，按照式 （7）～式 （9）求解

對于整個換熱過程中的總傳熱系數Kavg由式（10）求解

2 實驗結果及分析

2.1 流體換熱特性

圓管內水的換熱關聯(lián)式的研究較為成熟，實驗中將對殼程丙烷側換熱特性進行研究，通過整理實驗數據得到丙烷側傳熱系數，將該實驗值與Dittus－Boelter、Sieder－Tate、Miheev公式［11］和文獻 ［12］推薦的適用于管殼式換熱器的關聯(lián)式進行對比。

2.1.1 水平換熱段的換熱特性 結合實際FSRU冬天的運行背景，海水的水溫較低，故實驗中入口水溫設置為7℃，在不同的丙烷流速下進行實驗，其Nusselt數與Reynolds數的關系曲線如圖4所示，并將實驗值與各個經驗換熱關聯(lián)式進行對比。

從圖中可見，D－B、S－T、Miheev、適用于管殼式換熱器的關聯(lián)式計算值曲線隨著流速的增加均線性增長，實驗數值曲線也幾乎隨著流速的增加而線性增加，這是由于在整個換熱過程中，為了維持丙烷一直處于過冷液體，則丙烷的壓力較高，根據其物性，在不同的溫度下，其熱力參數幾乎不變化。

圖4 丙烷側換熱實驗值與經驗公式值對比Fig.4 Comparison between experimental data and empirical formula calculations

實驗值與D－B、S－T、Miheev換熱關聯(lián)式之間的誤差較大，在高流速時相對誤差更大。原因為丙烷側的換熱為殼程的換熱，其特征長度不是管內徑，而D－B、S－T、Miheev方程均適用于管內流體的換熱，其特征長度即為管內徑，故前文中的換熱器理論計算不能使用這3種經驗關聯(lián)式。

實驗值與管殼式換熱器計算值之間的誤差較小，而且隨著流速的增大反而減小，在低流速時相對誤差也在容許的誤差范圍之內，故前文中的換熱器理論計算應采用該關聯(lián)式進行計算。

2.1.2 傾斜角度下?lián)Q熱特性 在傾斜角度θ＝10°時，進行丙烷流動速度變化實驗，丙烷側Nusselt數與Reynolds數的關系曲線如圖5所示，并將實驗值與各個經驗換熱關聯(lián)式進行對比。

從圖中可見，D－B、S－T、Miheev方程計算值隨著流速的增加而線性增長。實驗值曲線幾乎平行于適用于管殼式換熱器的換熱關聯(lián)式曲線，與其相對誤差最小，而在θ＝0°實驗工況下，實驗值與管殼式關聯(lián)式計算值曲線相互交叉幾乎重合。

在θ＝10°時，實驗值與各個關聯(lián)式之間的相對誤差均有所增大，而且隨著流速的增大而增大，在小流動速度下與Miheev方程之間的相對誤差最小，而在θ＝0°實驗工況下，實驗值與管殼式關聯(lián)式之間的誤差最小，與各個關聯(lián)式計算值之間的相對誤差較小。

圖5 丙烷側換熱實驗值與經驗公式值對比Fig.5 Comparison between experimental data and empirical formula calculations

2.2 過冷中間流體換熱過程的影響因素

2.2.1 傾斜角度對換熱的影響 在Tw＝7℃、uw＝2.8m·s－1、up＝2.2m·s－1時，不同傾斜角度0°、5°、10°下，水側Nusselt數沿著管長方向的變化規(guī)律如圖6所示。

圖6 不同傾斜角度下水側Nusselt數沿著管長的變化Fig.6 Water’s local Nuunder differentθ

從圖中可以看出，當傾斜角度從0°增加到5°時，傳熱系數略微增大，但增大幅度不大，當傾斜角度從5°增加到10°時，傳熱系數幾乎沒有發(fā)生改變，兩條曲線幾乎重合。

相同實驗工況下，圖7為丙烷側Nusselt數沿著管長方向的變化規(guī)律。

圖中可以看出，隨著傾斜角度的增大，丙烷側的傳熱系數有所增加，從0°增加到5°時，傳熱系數的增幅較大，在5°到10°時增加的幅度減小。相對于水側換熱情況，丙烷側傳熱系數隨著傾斜角度的變大其增加幅度較大，說明傾斜對丙烷側的換熱影響較大。

圖7 不同傾斜角度下丙烷側Nusselt數沿著管長的變化Fig.7 Propane’s local Nuunder differentθ

2.2.2 水流動速度對換熱的影響 整個過冷中間流體換熱流程中，海水作為熱源，提供中間流體熱量，中間流體再將這部分熱量傳遞給LNG，而每天處理的LNG氣化量不一樣，則換熱流程中水的流量也不一樣，則研究水的流動速度對換熱過程的影響。

在不同傾斜角度θ＝0°、10°下，Nusselt數隨著流動速度的變化規(guī)律如圖8所示。

圖8 流動速度對換熱的影響Fig.8 Effect of velocity of water on heat transfer

從圖8中的兩條曲線可以看出，不管在水平還是傾斜換熱段中，水側Nusselt數均隨著水流動速度的增加而增加。流動速度的增加導致黏性底層的厚度變小，從而使得換熱強度變大，傳熱系數增加。

2.2.3 水入口溫度對換熱的影響 結合實際FSRU應用背景，LNG的氣化過程中所采用的熱源為海水，海水的溫度隨著季節(jié)的變化而變化，考慮到實際情況下，海水的溫度變化在7～22℃下進行變化，實驗中研究水溫對換熱段的影響。

不同傾斜角度θ＝0°、10°下，不同入口水溫時水的進出口溫度差如圖9所示。

圖9 進口水溫對水溫差的影響Fig.9 Effect of water’s temperature on water

從圖中可以看出，水的進出口溫度差隨著進口水溫的變化幾乎不發(fā)生改變，當進口水溫從7～22℃變化時，在傾斜角度下的水溫變化幅度差略微比水平角度下水溫差變化幅度大一些，但是增幅非常非常小，幾乎可以忽略不計。在同一個水溫下，發(fā)現(xiàn)傾斜角度10°下的水溫差明顯高于水平換熱段，增加幅度小于1℃。根據前文中傾斜角度對換熱的影響可知，隨著傾斜角度的增加，水側的換熱強度增強，傳熱系數增大，則相應的進出口的水溫差也比水平換熱段大。

丙烷作為中間流體，吸收海水的熱量后將這部分熱量傳遞給LNG使其氣化，丙烷／水換熱段中丙烷進出口溫差代表了丙烷所能夠提供LNG氣化所需的熱量，圖10為不同傾斜角度下，不同進口水溫時丙烷進出口溫差。

圖10 進口水溫對丙烷溫差的影響Fig.10 Effect of water’s temperature on propane

從圖中可以看出，當水溫度在7～22℃變化時，丙烷進出口溫差隨著進口水溫的升高而變大。結合水進出口溫差變化曲線，丙烷進出口溫差變化趨勢和水的進出口溫差變化趨勢相符合，因為丙烷／水換熱段中兩種流體的換熱量應該相等。同一進口水溫下，傾斜角度下丙烷進出口溫差比水平換熱段中的大，根據傾斜角度對換熱的影響來看，隨著傾斜角度的增大，換熱增強，則丙烷進出口溫度差也相應變大。

根據丙烷／水換熱段中換熱量相等

利用水側的進出口溫度差，求解得到丙烷進出口理論溫度差，圖11和圖12分別為θ＝0°和θ＝10°下，丙烷進出口理論計算值與實驗值的對比。

圖11 θ＝0°實驗值與理論值的對比Fig.11 Comparison of experimental data and calculated data underθ＝0°

圖12 θ＝10°實驗值與理論值的對比Fig.12 Comparison of experimental data and calculated data underθ＝10°

從圖11和圖12中可以看出，在水平換熱段中，根據換熱公式計算出的丙烷進出口溫差和實驗值相差不大，兩條曲線幾乎重合，說明實驗數據測量較為準確，而在傾斜角度θ＝10°下所計算的理論值與實驗值相差較為明顯，兩條曲線之間的差距在1.3℃左右，可能是由于換熱段傾斜，從而造成丙烷的入口溫度偏低，丙烷進出口溫度差變大，也可能是由于換熱段的傾斜，使得測量丙烷溫度的鉑電阻誤差變大，也可能是由于傾斜角度下，流體的混合更為充分，對溫度有所影響。

3 結 論

（1）進行過冷中間流體換熱實驗，利用實驗數據驗證換熱器理論計算時殼程流體換熱特性，實驗結果表明，適用于管殼式換熱器設計的換熱關聯(lián)式能夠較好地描述其換熱特性。

（2）結合FSRU實際應用背景，充分考慮海平面的晃動，從而進行靜態(tài)條件下，傾斜角度對換熱過程的影響實驗，實驗結果表明，傾斜角度的增加有利于換熱。

（3）結合FSRU實際應用背景，研究海水的質量流量對換熱過程的影響，實驗結果表明隨著海水流速的增大，換熱增強，海水的進出口溫度差減小。

（4）結合FSRU實際應用背景，一年四季海水溫度不同，研究海水入口溫度對換熱過程的影響，實驗結果表明，隨著海水入口溫度的升高，丙烷進出口溫度差變大，在傾斜角度下更為明顯，則說明提供給LNG的熱量更多。

符 號 說 明

cp——比定壓熱容，kJ·kg－1·K－1

G——質量流量，kg·s－1

h——傳熱系數，W·m－2·K－1

Q——換熱段的換熱量，kW

q——熱通量，W·m－2

T——溫度，K

Δz——換熱段長度，m

λ——熱導率，W·m－1·K－1

下角標

avg——平均量

i——第i換熱段

in——入口處

loc——局部換熱段

out——出口處

p——丙烷

w——水

［1］ Zellouf Y，Portannier B.First step in optimizing LNG storages for offshore terminals ［J］.NaturalGasScienceand Engineering，2011，3 （5）：582－590.

［2］ Park C，Lee C J，Lim Y.Optimization of recirculation operating in liquefied natural gas receiving terminal ［J］.TaiwanInstituteofChemicalEngineers，2010，41 （4）：482－491.

［3］ Ji Xin （紀 馨 ），Chen Shuangshuang （陳 雙 雙 ），Lin Wensheng（林文勝）.Research on single－phase heat transfer in intermediate fluid vaporizer／／Cryogenic Engineering Conference ［C］.Zunyi，2013.

［4］ Zou X，Gong M Q，Chen G F，Sun Z H，Zhang Y，Wu J F.Experimental study on saturated flow boiling heat transfer of R170／R290mixtures in a horizontal tube ［J］.InternationalJournalofRefrigeration，2010，33：371－380.

［5］ Lee H S，Yoon J I，Kim J D，Bansal P K.Characteristics of condensing and evaporating heat transfer using hydrocarbon refrigerants ［J］.AppliedThermalEngineering，2006，26：1054－1062.

［6］ Gu H F，Li H Z，Wang H J.Experimental investigation on convective heat transfer from a horizontal miniature tube to methane at supercritical pressures ［J］.AppliedThermal Engineering，2013，58 （1）：490－498.

［7］ Kim S Y，Kang B H，Hyun J M.Heat transfer in the thermally developing region of a pulsating channel flow ［J］.InternationalJournalofHeatTransferandMassTransfer，1997，7 （10）：2461－2466.

［8］ Han Xiaoxia （韓 曉 霞），Nan Xiaohong （南 曉 紅），Liu Xianding （劉 咸 定 ），F(xiàn)erreira C A I.In－tube evaporation pressure－drop characteristics of R290and R404A ［J］.JournalofRefrigeration（制 冷 學 報），2006，27 （2）：17－20.

［9］ Ning Jinghong （寧靜紅），Liu Shengchun （劉圣春），Ye Qingyin （葉慶銀）.Analysis on heat transfer characteristic of dropwise condensation of natural refrigerant R290and R22［J］.CIESCJournal（化 工 學 報 ），2013，64 （8）：2827－2832.

［10］ Ning Jinghong （寧靜紅），Liu Shengchun （劉圣春），Li Huiyu （李 慧 宇）.Simulation on condensation heat transfer characteristic of natural refrigerant R290inside tube ［J］.CIESCJournal（化工學報），2012，63 （4）：1038－1043.

［11］ Yang Shiming （楊 世 銘），Tao Wenquan （陶 文 銓）.Heat Transfer（傳熱學） ［M］.Beijing：Higher Education Press，2006：246－322.

［12］ Qian Songwen （錢頌文）.Handbook of Heat Transfer Design（換熱器設計手冊）［M］.Beijing：Chemical Industry Press，2002：59－67.