摘要：對于LNG空溫式氣化器的單根翅片管進(jìn)行模擬換熱計算，計算時將翅片管分為單項液、兩項以及單項氣3區(qū)，得出了管外壁溫度、管內(nèi)流體溫度以及管內(nèi)流體表面?zhèn)鳠嵯禂?shù)沿管長的分布。

關(guān)鍵詞：液化天然氣；空溫式氣化器；換熱計算；三區(qū)模型

Abstract：The heat transfer calculation of a single finned-tube in ambient air Vaporinzer is simulated．The sing1e finned-tube is divided into 3 areas including single-phase fluid,two phases and single-phase gas. The distributions of external wall temperature，intemal fluid temperature and the surface heat transfer coefficient of the internal fluid along the length are obtained.

Key words：LNG；ambient air vaporizer；heat transfer calculation；three—area model

1概述

隨著經(jīng)濟(jì)的持續(xù)穩(wěn)步增長，我國的能源需求量不斷增大。據(jù)BP世界能源統(tǒng)計結(jié)果，2000-2008年我國各類一次能源的年均增長率中，天然氣的增幅達(dá)到了l5％，明顯高于煤炭和石油^[1]。常壓下將天然氣冷卻至-162℃時，氣態(tài)天然氣轉(zhuǎn)變?yōu)橐夯烊粴?span lang="EN-US">(LNG)，其體積約為同質(zhì)量氣態(tài)天然氣的1／600。天然氣以LNG的形式儲存和運輸，具有成本低、使用方便、安全可靠等諸多優(yōu)點。近幾年來，

我國陸續(xù)規(guī)劃和建造了多座大型LNG接收終端、LNG生產(chǎn)廠以及LNG衛(wèi)星站，我國LNG產(chǎn)業(yè)進(jìn)入快速發(fā)展階段，LNG產(chǎn)業(yè)鏈日趨成熟。

對各種LNG氣化器的經(jīng)濟(jì)以及環(huán)境對比顯示，空溫式氣化器(ambient air vaporizers，簡稱AAVs)以其節(jié)能環(huán)保的優(yōu)勢得到了大力提倡^[2]，在LNG氣化

站中AAVs已經(jīng)得到了廣泛的應(yīng)用，近年來大型的LNG接收終端也開始采用AAVs來實現(xiàn)LNG的氣化。

目前國內(nèi)外學(xué)者對AAVs的研究已經(jīng)取得一定的成果。Lee^[3]與Kong^[4]分別通過實驗與數(shù)方法研究了翅片管形式、翅片角、空氣側(cè)參數(shù)等對低溫領(lǐng)域AAVs換熱性能的影響，Hyo-Min Jeong^[5]等用數(shù)值模擬方法研究了LNG空溫式氣化器外表面霜層厚度與翅片厚度以及翅片角之間的關(guān)系，Hyo-Min Jeong和Han-Shik Chung^[6]對目前應(yīng)用較廣泛的8fin501e型與4fin751e型翅片管在結(jié)霜工況下?lián)Q熱特性進(jìn)行了實驗比較，國內(nèi)的劉小川^[7]利用FLUENT軟件對換熱器結(jié)霜工況下空氣側(cè)的傳熱傳質(zhì)狀況進(jìn)行了模擬，陳瑞球等^[8]通過實驗分析了風(fēng)速、肋片間距以及空氣濕度等因素對結(jié)霜工況翅片管空氣側(cè)換熱性能的影響。另外，高華偉^[9]等針對LNG空溫式氣化器的換熱過程，將管內(nèi)的氣化分為單相液、單相氣與兩相區(qū)3段，并建立了單相段天然氣的換熱模型，進(jìn)行了數(shù)值模擬。

由上述可以看出，以往對AAVs的研究主要集中在空氣側(cè)特性的研究，而對LNG這一特定介質(zhì)的管內(nèi)氣化過程的換熱研究較少。本文對LNG空溫式氣化器單根翅片管的換熱進(jìn)行計算，并分析管內(nèi)的相變吸熱特點。

2 LNG空溫式氣化器結(jié)構(gòu)

AAVs是由翅片管按照一定的間距并聯(lián)而成的，一般是單程式，實際運行中的LNG空溫式氣化器見圖1^[10]。換熱管一般較長，為了增大空氣側(cè)的換熱面積，在換熱管的外側(cè)加裝翅片。8翅片星型翅片管見圖2^[10]，目前最常用的是8翅片結(jié)構(gòu)，另外還有l2翅片與4翅片結(jié)構(gòu)。

本文以實際運行中的某空溫式氣化器為例，進(jìn)行單根翅片管的換熱計算，計算時將翅片管分為單相液、兩相以及單相氣3區(qū)。氣化器的工作壓力為1．6 MPa，每根翅片管均為l2翅片結(jié)構(gòu)，單根翅片管長度為13．316 m，內(nèi)徑d₁為10 mm，外徑d₂為l4mm，翅片高度h為80 mm，翅片厚度δ為2 mm，翅片管的材質(zhì)為鋁合金。LNG中甲烷、乙烷與丙烷的摩爾分?jǐn)?shù)分別為90．07％、7．36％與2．57％。LNG進(jìn)入氣化器的溫度為111 K。

3 LNG管內(nèi)氣化過程分析

根據(jù)文獻(xiàn)[11]，LNG管內(nèi)氣化機(jī)理為：當(dāng)混合物進(jìn)入氣化器后，沸點低、蒸氣壓高的組分先氣化，但由于液相斷斷續(xù)續(xù)和傳熱面接觸以及液體被氣泡破裂時分裂成的氣沫所氣化，導(dǎo)致各種組分都得到了氣化。

LNG是由不同烴類組成的非共沸混合物，開始?xì)饣臏囟?span lang="EN-US">(即泡點)與完全氣化的溫度(即露點)不同，氣化過程中溫度不斷發(fā)生變化。LNG管內(nèi)氣化過程可分為過冷區(qū)、兩相區(qū)和過熱區(qū)。當(dāng)溫度低于泡點時，處于單相液體換熱區(qū)；高于泡點而低于露點時，處于兩相區(qū)；高于露點時，處于單相氣體換熱區(qū)^[12]。

4翅片管換熱計算

4．1  翅片管換熱計算簡化假設(shè)

①剛剛進(jìn)入翅片管的低溫LNG與氣化器中原有LNG的混合是在一瞬間完成的，即翅片管中LNG的溫度與各組分的組成在液體內(nèi)部處處均勻。

②  由于氣化器工作時間比較長，故采用穩(wěn)態(tài)分區(qū)計算方法。

③認(rèn)為氣化器運行中，壓力是恒定的^[8]。

④管流采用一維近似。

⑤空氣按干空氣處理。

⑥管外空氣側(cè)按大空間自然對流換熱處理。

4．2泡點與露點計算

要將翅片管進(jìn)行分區(qū)，首先需要計算LNG在工作壓力下的泡點與露點。在LNG的吸熱過程中，開始產(chǎn)生第一個氣泡時的溫度稱為泡點^[13]，定義式為：

最后一個液滴消失時的溫度稱為露點^[13]，定義式為：

利用式(1)～(3)編制Matlab程序，計算得到該算例下LNG的泡點為l63.7 K，露點為211.2 K。

4．3單相液(氣)區(qū)換熱計算

按從管內(nèi)到管外的順序，翅片管的熱量傳遞依次為管內(nèi)流體與管壁的強(qiáng)制對流換熱、通過管壁的導(dǎo)熱、管外空氣與翅片和管壁的對流換熱。

①管內(nèi)流體的對流換熱

管內(nèi)LNG或氣態(tài)天然氣的對流換熱可以按照內(nèi)部強(qiáng)制對流換熱的實驗準(zhǔn)則關(guān)聯(lián)式進(jìn)行計算。

由管內(nèi)對流換熱的準(zhǔn)則關(guān)系式(8)、(9)確定出管內(nèi)流體側(cè)的表面?zhèn)鳠嵯禂?shù)α₁：

由此可以確定出管內(nèi)單位長度的對流換熱量q_l：

②通過管壁的導(dǎo)熱

利用通過圓筒壁的導(dǎo)熱公式，計算沒有翅片部分管外壁的溫度T_w，此溫度即為翅片的根部溫度。

③管外空氣側(cè)的對流換熱

空氣側(cè)換熱可看作大空間自然對流換熱，單個翅片的換熱可以看作空氣與等截面直肋的換熱。

首先由空氣的定性溫度T_m確定其物性參數(shù)，T_m計算公式為：

由于翅片管為立式安裝，故此處自然對流換熱計算選用豎壁準(zhǔn)則關(guān)聯(lián)式：

利用式(4)～(19)編制Matlab程序進(jìn)行換熱計算，計算思路為：首先假設(shè)管內(nèi)壁溫度T_n，利用式(4)～(10)計算出微元段管內(nèi)單位長度的對流換熱量q₁，然后利用式(11)計算出管外壁溫度T_w，再利用式(12)～(19)計算出空氣側(cè)單位長度的自然對流換熱量q₂。若，說明所假設(shè)的管內(nèi)壁溫度T_n不正確，需重新假設(shè)計算；若，說明假設(shè)的管內(nèi)壁溫度T_n正確，即可以進(jìn)行下一步長的計算。

取空氣溫度為l0℃，步長(進(jìn)行一次計算所取的管道長度)L=0.1 m，同時根據(jù)LNG達(dá)到泡點(163.7 K)時代表單相液區(qū)結(jié)束、天然氣溫度處于露點(211.2 K)與263 K之間時為單相氣區(qū)，也可以計算出單相液區(qū)和單相氣區(qū)的長度。

4．4兩相區(qū)換熱計算

由于LNG是多元組分混合物，而混合物的沸騰換熱很復(fù)雜，因而是較難處理的，目前對二元混合物有一些沸騰傳熱系數(shù)的計算公式，但其中均含有與特定二元物質(zhì)組合有關(guān)的實驗系數(shù)，其流動換熱關(guān)系式不能直接應(yīng)用于本課題的LNG介質(zhì)。因此，本文在進(jìn)行兩相區(qū)的換熱計算時，采用簡化假設(shè)。由于LNG中甲烷的摩爾分?jǐn)?shù)高達(dá)90％以上，故將此區(qū)的LNG看作甲烷。將兩相區(qū)分為流動沸騰區(qū)和缺液區(qū)，在流動沸騰區(qū)溫度保持泡點不變，在缺液區(qū)溫度由泡點升高至露點。

①流動沸騰區(qū)換熱計算

常壓下沸點小于120 K的流體稱為低溫流體，其熱物性相似性判據(jù)J≥3.5(J=100 p_r，p_r指對比溫度為0.625時的對比壓力)^[14]。根據(jù)此定義，甲烷屬于低溫流體。Klimenko的方法是目前計算低溫流體流動沸騰換熱最精確的關(guān)系式^[15]，該方法的具體表述如下：

②缺液區(qū)換熱計算

由于在兩相區(qū)熱流密度變化不大，所以在缺液區(qū)取與流動沸騰區(qū)相同的熱流密度值和管內(nèi)流體表面?zhèn)鳠嵯禂?shù)值。缺液區(qū)管內(nèi)流體的溫度計算公式為：

流動沸騰區(qū)和缺液區(qū)的管外壁溫度計算與單相液(氣)區(qū)計算相似。

5計算結(jié)果分析

根據(jù)以上3個區(qū)的計算結(jié)果，可得到以下結(jié)論：

①經(jīng)計算得翅片管的總長度為13.45m，與翅片管的實際長度l3.316m接近。其中單相液區(qū)最短，只占總長度的14％，兩相區(qū)占45％，單相氣區(qū)占41％。

②管內(nèi)流體表面?zhèn)鳠嵯禂?shù)沿翅片管的分布見圖3。3區(qū)中，兩相區(qū)的管內(nèi)流體表面?zhèn)鳠嵯禂?shù)最大，接近2 200 W／(m²·K)，單相液區(qū)居中，單相氣區(qū)最小。由此可見，相變的發(fā)生使得換熱強(qiáng)化，在較小的傳熱溫差下即可產(chǎn)生較大的傳熱熱流。并且由兩相區(qū)的計算得到，在本文的熱流強(qiáng)度、工作壓力及流量下，N_CB始終小于l.2×10⁴，說明甲烷始終處于核態(tài)沸騰換熱區(qū)，沒有出現(xiàn)液膜強(qiáng)制對流。

③管外壁溫度與流體溫度沿翅片管的分布見圖4。流體溫度在單相液區(qū)急劇升高；在兩相區(qū)的核態(tài)沸騰區(qū)保持不變，到達(dá)缺液區(qū)后不斷升高；在單相氣區(qū)，相對于單相液區(qū)緩慢升高。

管外壁溫度在單相液區(qū)急劇升高，到達(dá)兩相區(qū)后有一個突降，而后保持不變，到達(dá)缺液區(qū)后逐漸升高，在單相氣區(qū)也是緩慢升高。由管外壁溫度沿翅片管的分布可以看出，在單相液區(qū)和兩相區(qū)的核態(tài)沸騰區(qū)，即從人口至管長大約6m處，管外壁溫度較低，容易發(fā)生結(jié)露或結(jié)霜現(xiàn)象。

6結(jié)論

本文對LNG空溫式氣化器單根翅片管的換熱進(jìn)行了模擬計算，計算時沿管長方向?qū)⒊崞芊譃閱蜗嘁骸上嗯c單相氣3區(qū)，得出了管外壁溫度、LNG(或氣態(tài)天然氣)溫度以及管內(nèi)流體表面?zhèn)鳠嵯禂?shù)沿管長的分布。翅片管的換熱計算是進(jìn)行氣化器設(shè)計的關(guān)鍵，通過進(jìn)行換熱過程的模擬計算，可以確定設(shè)備是否能達(dá)到預(yù)定的氣化能力，并且可以分析LNG的組成變化或環(huán)境條件改變對換熱器氣化效果的影響。本文的計算結(jié)果可以為LNG空溫式氣化器的設(shè)計與優(yōu)化提供依據(jù)。

參考文獻(xiàn)：

[1]  付國忠，陳超．我國天然氣供需現(xiàn)狀及煤制天然氣工藝技術(shù)和經(jīng)濟(jì)性分析[J]．中外能源，2010，15(6)：28-34．

[2] 陳雪．LNG接收終端工藝對比及選擇[J]．石油規(guī)劃設(shè)計，2008，19(2)：44-47．

[3] LEE K S，KIM W S．The effects of design and operating factors on the frost growth and thermal performance of a flat plate fin-tube heat exchanger under the frosting condition『J]．KSME International Journal，1999，13(12)：973-981．

[4] KONG T W．LEE S C，LEE Y H，et al．A study on the air vaporizer for liquefied natural gas with super low temperature[C]／／Proceedings of the 3rd Asian Conference on Refrigeration and Air-conditioning. Gyeongju (South Korea)：[s.n.]，2006．

[5] JEONG H．CHUNG H．Optimum design of vaporizer fin with liquefied natural gas by numerical analysis[J].KSME International Journal，2006，20(4)：545-553．

[6] JEONG H，CHUNG H．Experimental study on the characteristics of longitudinal fin air-heating vaporizers in different seasons[J]．Journal of Mechanical Science and Technology，2008，22(5)：981-990．

[7]  劉小川．結(jié)霜工況下翅片管換熱器傳熱傳質(zhì)的數(shù)值模擬(碩士學(xué)位論文)[D]．上海：上海交通大學(xué)，2007：9-68．

[8]  陳瑞球，周廣．霜形成對翅片管換熱器空氣側(cè)表面?zhèn)鳠嵯禂?shù)影響的實驗研究[J]．制冷技術(shù)，2010，38(4)：54-56．

[9]  高華偉，段常貴，解東來，等．LNG空溫式氣化器氣化過程的數(shù)值分析[J]．煤氣與熱力，2008，28(2)：Bl9-B22

[10] GAVELLI F．Computational fluid dynamics simulation of fog clouds due to ambient air vaporizers[J]．Journal of Loss Prevention in the Process Industries，2010，23(6)：773-，780．

[11]  姜正侯．燃?xì)夤こ碳夹g(shù)手冊[M]．上海：同濟(jì)大學(xué)出版社，1993：230-254．

[12]  田貫三，朱斌．LPG氣化壓力與溫度用CSD方程的分析[J]．煤氣與熱力，2000，20(6)：448-450.

[13]  寇虎，嚴(yán)銘卿，廉樂明，等．液化石油氣泡點的直接計算及其應(yīng)用[J]．煤氣與熱力，2001，21(5)：443-445．

[14] KLIMNEKO V V．Heat transfer intensity at forced flow boiling of cryogenic liquids[J]．Cryogenic，1982，22 (6)：569-576．

[15]  KIIMNEKO V V．Generalized correlation for two-phase forced flow heat transfer-second assessment[J]．  International Journal of Heat and Mass Transfer，1990，33 (10)：2073-2088．

本文作者：楊聰聰  焦文玲  楊 光  趙元偉  馬 建

作者單位：大慶高新區(qū)規(guī)劃建筑設(shè)計院  哈爾濱工業(yè)大學(xué)市政環(huán)境工程學(xué)院  深圳市燃?xì)饧瘓F(tuán)股份有限公司  大慶油田有限責(zé)任公司