摘 要：超級(jí)開架式氣化器(SuperORV)是一種以海水為熱源的新型氣化器，它采用雙層結(jié)構(gòu)的傳熱管，可有效改善傳統(tǒng)開架式氣化器管束外結(jié)冰的狀況并提高換熱效率，而目前國(guó)內(nèi)對(duì)該裝置傳熱性能的研究還較少。為此，對(duì)SuperORV關(guān)鍵傳熱單元——傳熱管的換熱過程進(jìn)行了模擬研究，建立了SuperORV傳熱管整體換熱過程的傳熱計(jì)算模型。該模型利用兩組離散化方程組分別描述了SuperORV傳熱管氣化段和加熱段的傳熱過程，在給定的尺寸和邊界條件下對(duì)傳熱管的整體換熱性能進(jìn)行了數(shù)值模擬，得到了傳熱管各個(gè)局部的表面換熱系數(shù)和溫度分布曲線，進(jìn)而推導(dǎo)出了傳熱管總換熱系數(shù)和熱流密度的分布曲線。海水和外翅片管上的溫度分布曲線可用于預(yù)測(cè)傳熱管外表面易結(jié)冰的位置，傳熱管總換熱系數(shù)和熱流密度的分布曲線則可為傳熱管的整體換熱性能描述提供幫助。該模型及相關(guān)模擬分析可望為該類氣化器的設(shè)計(jì)、選型和運(yùn)行管理提供參考。
關(guān)鍵詞：超級(jí)開架式氣化器 LNG接收終端 氣化 低溫傳熱 沸騰傳熱 數(shù)值模擬 表面換熱系數(shù) 溫度分布曲線
Numerical simulation analysis of a heat transfer tube in SuperORV
Abstract：Super Open-Rack LNG Vaporizer(SuperORV)，a new type of vaporizer using sea water as a heat source，has a doubledtube structure which can so markedly weaken icing conditions outside the tube as to improve the overall heat transfer efficiency．However，such related research in China has not yet been found on this subject．Therefore，in this paper,a simulation studv was conducted to investigate the LNG(Liquefied Natural Gas)evaporating process in a SuperORV heat transfer tube．At first，calculation models were established of the overall heat transfer process．Two groups of discrete equations were then utillzed to describe the heat transfer process in vaporizing and heating sections of this equipment．Under the given size and boundary conditions，a numerical study was carried out of the performance of the heat transfer tube and the heat exchange coefficient and temt)crature distribution curves were obtained for the surface of each zone，based on which the distribution curves of heat flux density and cumulative heat exchange coefficient were thus derived．These curves would be helpful for the description of overall heat exchanging performance of the tube；on the other hand，the temperature distribution curve of sea water vs．the exterior fin tube could be used to predict the freeze position on the external surface of the heat transfer tube．This study provides reference for the design，type selection and operation management of LNG SuperORV facilities．
Key words：SuperORV，LNG receiving terminal，LNG evaporation，cryogenic heat transfer，boiling heat transfer，simulation
開架式氣化器(ORV)是一種以海水為熱源的氣化器，是液化天然氣(LNG)接收終端中用于基本負(fù)荷的大型氣化裝置，但由于海水容易在該氣化裝置底部結(jié)冰，使得裝置盼傳熱性能下降，導(dǎo)致其氣化性能無(wú)法得到進(jìn)一步的提升。超級(jí)開架式氣化器(SuperORV)是在ORV基礎(chǔ)上二的改進(jìn)版本，其傳熱管可分為氣化段和加熱段兩部分(圖1)。加熱段中采用雙層結(jié)構(gòu)的傳熱管，使得管外結(jié)冰的狀況得到了進(jìn)有效的改善。據(jù)有關(guān)技術(shù)數(shù)據(jù)[1-2]，SuperORV可以使氣化裝置的氣化能力提高3倍，海水流量和安裝空間可分別減小15％和40％。
我國(guó)先前建成的粵東LNG接收終端以及其他一些在建的LNG接收終端項(xiàng)目多采用的是傳統(tǒng)的ORV[3-6]，而日本和歐洲的LNG接收終端早在2005年以前都已大量采用SuperORV作為基本負(fù)荷型氣化設(shè)備(表1) [7]。未來(lái)中國(guó)將會(huì)大量建設(shè)LNG接收終端[8-9]，但國(guó)內(nèi)對(duì)該裝置傳熱性能的研究目前還不是很多。為此，筆者擬對(duì)SuperORV關(guān)鍵傳熱單元傳熱管的換熱過程進(jìn)行數(shù)值模擬，并對(duì)模擬結(jié)果反映的氣化器性能進(jìn)行分析和討論，期望能為該類型氣化器的設(shè)計(jì)、選型和運(yùn)行管理提供參考。
2 SuperORV傳熱管傳熱計(jì)算模型
2．1 傳熱管的傳熱過程分析
SuperORV上部裝有海水噴淋裝置，將海水自管束板外自上而下噴淋，海水經(jīng)分布器分配后，形成薄膜均勻沿管束下降，使管內(nèi)的LNG受熱氣化。在傳熱管的加熱段，LNG從底部的分配器先后進(jìn)入內(nèi)套管和外管之間的環(huán)狀間隙豎直向上流動(dòng)。間隙中的LNG流量較小，在流動(dòng)過程中被翅片管外的海水加熱后立即氣化。內(nèi)套管中流動(dòng)的LNG被間隙里已經(jīng)氣化的天然氣(NG)加熱，氣化逐漸進(jìn)行。隨著含氣率的增加，內(nèi)套管中的LNG流體先后經(jīng)歷不同的流態(tài)，可被劃分為不同的換熱區(qū)間：?jiǎn)蜗嘁后w對(duì)流換熱區(qū)、欠熱沸騰換熱區(qū)、飽和沸騰換熱區(qū)和缺液區(qū)[10-11]。待內(nèi)套管中的氣化過程完畢后，內(nèi)外管的NG混合后進(jìn)入傳熱管的加熱段；在加熱段，翅片管內(nèi)的低溫NG單相氣體通過管壁，被海水液膜加熱到設(shè)計(jì)要求的出口溫度，最后離開傳熱管。
為了方便傳熱區(qū)間的劃分，根據(jù)LNG流體的焓值變化，沿?fù)Q熱管將管段劃分成若干個(gè)換熱單元，每個(gè)小單元里的流體進(jìn)行流型和狀態(tài)的判斷，作為選擇換熱經(jīng)驗(yàn)公式的依據(jù)。對(duì)小單元進(jìn)行能量守衡的恒算，根據(jù)已知工況，計(jì)算出單元出口(單元人口或出口的叫法取作與管內(nèi)LNG流動(dòng)方向相同，與海水流動(dòng)方向相反，下同)的溫度、焓值等數(shù)據(jù)，每個(gè)單元出口的數(shù)據(jù)作為下個(gè)單元入口的已知量，海水和LNG進(jìn)出口焓值作為邊界條件，依次計(jì)算，直到天然氣出口滿足出口溫度要求，從而得到每個(gè)單元的各種參數(shù)以及其他相應(yīng)參數(shù)沿傳熱管的分布情況。
由于流體的氣化過程比較復(fù)雜，為了更好地分析傳熱管的傳熱性能，筆者僅考慮翅片管壁不結(jié)冰狀況下的傳熱過程，并對(duì)傳熱過程做如下簡(jiǎn)化：
1)傳熱管是光管，不考慮強(qiáng)化傳熱措施和污垢熱阻的影響。
2)管道內(nèi)壓力恒定，LNG在恒壓條件下沸騰。
3)LNG沿管道長(zhǎng)度方向上的飽和溫度不變。
4)不考慮流體和管道厚度，即忽略它們?cè)趶较蛏系臏囟忍荻葘?duì)傳熱的影響。
5)進(jìn)人環(huán)狀間隙的LNG在進(jìn)入瞬間即開始?xì)饣瑑?nèi)套管中的LNG離開氣化段時(shí)剛好被全部氣化。
6)LNG氣化時(shí)，不考慮管壁的過熱度，即假設(shè)壁面過熱度為0。
7)環(huán)狀間隙和內(nèi)套管中的兩股流體混合后的焓值與混合前兩股流體的焓值之和相等。
2．2 傳熱計(jì)算模型
由于傳熱管氣化段和加熱段結(jié)構(gòu)是不同的，對(duì)這兩部分的傳熱過程分別建立傳熱模型(圖2、3)，其中表示與相鄰壁面的換熱量。
2．2．1氣化段
對(duì)管外海水側(cè)能量恒算得到離散化方程：
式中m表示質(zhì)量流量，kg／s；h表示進(jìn)出口流體的比焓值，J／kg；D代表傳熱管直徑，m；Dz表示單元長(zhǎng)度，m；a為對(duì)流傳熱表面換熱系數(shù)，W／(m2·K)；T為溫度，K；上標(biāo)1、2分別代表單元進(jìn)出口截面；下標(biāo)中的w代表海水，I表示翅片管壁，o表示其外側(cè)。
對(duì)環(huán)狀間隙里NG的能量恒算得到離散化方程為：
式中下標(biāo)g代表NG，Ⅱ表示翅片內(nèi)側(cè)和套管管壁，i表示管內(nèi)側(cè)，其余同上。
對(duì)內(nèi)套管LNG能量恒算得到離散化方程為：
式中下標(biāo)m代表LNG兩相流體，其余同上。
上述3個(gè)能量守恒方程與通過2個(gè)管壁的能量分別相等列出的方程聯(lián)立，可得到下面的方程組：
以上方程組中，由于工藝參數(shù)mw、mg、mm、DⅠ,o、DⅠ,i和DⅡ,i都系已知，由邊界條件和上一單元的計(jì)算結(jié)果，也可知hw、hg和hm，計(jì)算壓力下對(duì)應(yīng)的Tw、Tg和Tm也已知，根據(jù)單元?jiǎng)澐值囊罁?jù)，出口的LNG比焓值hm也為已知量，而各個(gè)表面對(duì)流換熱系數(shù)依賴于流體的物性參數(shù)，而物性參數(shù)的取值取決于此處的溫度和壓力，由于入口溫度已知，將其取作單元入口溫度對(duì)應(yīng)的物性。所以，方程組中獨(dú)立的未知量?jī)H有Dz、TⅠ、TⅡ、hw和hg這5個(gè)未知量，與方程數(shù)目相等，方程組封閉，可解。
2．2．2加熱段
仿照上述方法，得到加熱段各單元的傳熱模型為：
式中下標(biāo)mg代表該段的NG氣體，wall代表管壁，o代表管外側(cè)，i表示管內(nèi)側(cè)，其余同上。
方程中的未知量為單元長(zhǎng)度Dz、壁溫Twall和單元出口的海水比焓值hw，方程同樣封閉，可解。
2．3 傳熱經(jīng)驗(yàn)公式的選擇
2．3．1海水側(cè)液膜的換熱系數(shù)
水膜下降過程中的換熱系數(shù)比較高，但目前對(duì)于豎直降膜的傳熱經(jīng)驗(yàn)公式多集中在對(duì)液膜冷卻換熱設(shè)備的研究，很少有適用于液膜加熱低溫流體工況的公式。考慮到液膜換熱機(jī)理復(fù)雜，尤其當(dāng)液膜的流道表面比較復(fù)雜的情況(如外翅片管)下，液膜傳熱關(guān)聯(lián)式的適用條件比較苛刻，模擬選用的換熱系數(shù)為常數(shù)5800W／(m2·K)，該數(shù)據(jù)源于跟該類氣化器有關(guān)的本文參考文獻(xiàn)[12]。
2．3．2單相流體的對(duì)流換熱
單相流體的對(duì)流換熱適用于在氣化段內(nèi)套管中過冷LNG液體、環(huán)狀間隙以及加熱段中的NG氣體。單相強(qiáng)迫對(duì)流傳熱系數(shù)用aspl表示。初步計(jì)算顯示，SuperORV中的單相流動(dòng)主要處于湍流狀態(tài)，考慮采用經(jīng)典的Dittus-Boelter關(guān)系式[13]：
Nu=0.023Re0.8Pr0.4 (6)
2．3．3欠熱沸騰換熱區(qū)的對(duì)流換熱系數(shù)
欠熱沸騰換熱區(qū)的傳熱公式采用Rohesnow方法[10-11]。該換熱區(qū)的傳熱量被看作是單相對(duì)流貢獻(xiàn)部分qspl和壁面過熱的泡核沸騰傳熱貢獻(xiàn)部分qsub之和。其中，qspl=aspl(Tw-Tf)，可根據(jù)對(duì)流傳熱Dittus-Boelter關(guān)系式計(jì)算aspl，qsub則可使用Rohesnow的池內(nèi)沸騰換熱公式獲得。具體公式和參數(shù)的選擇見本文參考文獻(xiàn)[10-11]。
2．3．4飽和沸騰換熱區(qū)的對(duì)流換熱系數(shù)
選用Chen關(guān)系式來(lái)描述飽和沸騰換熱區(qū)的傳熱特征[10-11]。Chen關(guān)系式最適合描述非金屬純流體的，飽和沸騰，其對(duì)有機(jī)流體的飽和沸騰描述適用性也很強(qiáng)。該方法認(rèn)為，飽和泡核沸騰區(qū)內(nèi)存在2種基本傳熱模式：泡核沸騰傳熱和強(qiáng)制對(duì)流傳熱。在飽和沸騰下，兩相流的傳熱系數(shù)(aTP)可由下式得到：
aTP＝amac +amic (7)
式中amac描述對(duì)流傳熱，稱宏觀傳熱分量；amic描述泡核沸騰傳熱，稱微觀傳熱分量。Chen公式中參數(shù)的選擇和計(jì)算方法見本文參考文獻(xiàn)[10-11]。
2．3．5缺液區(qū)的對(duì)流換熱系數(shù)
缺液區(qū)以液滴在氣體的彌狀流動(dòng)為主。關(guān)于彌狀流的傳熱公式選擇下式[11]：
式中的選擇見本文參考文獻(xiàn)[11]。
在加熱段，模型根據(jù)條件判斷出不同的傳熱區(qū)間，從而選擇特定傳熱區(qū)間的公式進(jìn)行計(jì)算。
3 計(jì)算結(jié)果分析
3．1主要技術(shù)參數(shù)及邊界條件
翅片管外徑為40mm，內(nèi)徑為20mm；內(nèi)套管外徑為l8mm，內(nèi)徑為l4mm。根據(jù)本文參考文獻(xiàn)[1]中整臺(tái)氣化器的蒸發(fā)能力，換算得到海水噴淋的流量為2.5kg／s；LNG總流量為0.05kg／s(內(nèi)套管內(nèi)LNG流量為0.036kg／s，夾層里的氣態(tài)NG流量為0.014kg／s)。根據(jù)氣化工藝要求，LNG入口溫度為170K，出口溫度為275K，管內(nèi)操作壓力為4MPa；海水不結(jié)冰，其主流出口溫度為280K。
海水和LNG都是混合物，為方便起見，在模擬中直接采用它們主要成分(水和甲烷)的物性參與計(jì)算。在軟件Matlab中直接引用美國(guó)國(guó)家標(biāo)準(zhǔn)研究院(NIST)物性軟件REFPROP中的數(shù)據(jù)，進(jìn)而計(jì)算得到相應(yīng)的結(jié)果。
3．2 結(jié)果及分析
計(jì)算得到傳熱管氣化段為2.80m，加熱段為6.56m，傳熱管總長(zhǎng)為9.36m，該長(zhǎng)度比實(shí)際管長(zhǎng)偏大17％，這主要是由于在計(jì)算中并未考慮強(qiáng)化傳熱措施的影響。計(jì)算得到海水入口處的溫度為282.5K，工程案例中人口海水的溫度為283K[7]。
圖4給出了氣化段內(nèi)套管中LNG表面對(duì)流換熱系數(shù)的分布曲線，圖5給出了內(nèi)套管壁及其兩側(cè)流體的溫度分布曲線。從圖4可以看出，飽和沸騰區(qū)的傳熱系數(shù)有一個(gè)很大的提升，正是南于這個(gè)原因，圖5中內(nèi)套管的壁溫會(huì)有所下降，而偏向低溫流體。在缺液區(qū)，圖4曲線上出現(xiàn)了較大的下降臺(tái)階，這是由擘面被蒸干所導(dǎo)致的傳熱性能下降而引起的，可以看到，與其對(duì)應(yīng)的圖5中的內(nèi)套管壁溫也有一個(gè)明顯的上升。環(huán)狀間隙里的NG流體的溫度在0.5m之后保持在了一個(gè)較穩(wěn)定的溫度位上，這也是SuperORV利用環(huán)隙中的NG進(jìn)行保溫，緩解了結(jié)冰狀況的依據(jù)。
圖6給出了整個(gè)換熱管上翅片管壁和海水的溫度分布曲線。海水流量較之LNG流量大得多，所以其溫度變化不大。在1m以下的位置，管壁的溫度急劇降低，甚至低于冰點(diǎn)，很容易導(dǎo)致靠近壁面的海水結(jié)冰。該區(qū)域確實(shí)也是結(jié)冰情況較嚴(yán)重的區(qū)域。
圖7、8分別給出了反映SuperORV整體換熱性能的總換熱系數(shù)和熱流密度的分布曲線?？梢钥吹剑跉饣魏图訜岫芜^渡區(qū)，2條曲線都出現(xiàn)了跳躍。這是由于SuperORV結(jié)構(gòu)上的變化引起的。圖7中，氣化段的總換熱系數(shù)總體優(yōu)于加熱段的總換熱系數(shù)。圖8中，加熱段開頭的熱流密度先高于氣化段末尾的熱流密度而后逐漸下降，這是由NG流體混合前后的溫度變化引起的，從氣化段出來(lái)的2個(gè)流道里的流體混合后導(dǎo)致溫度有大的下降，而后隨著加熱溫度又逐漸上升。
4 結(jié)束語(yǔ)
建立了SuperORV傳熱管的整體換熱過程傳熱計(jì)算模型，該模型利用兩組離散化方程組分別描述了SuperORV傳熱管氣化段和加熱段的傳熱過程，并在給定的尺寸和邊界條件下對(duì)傳熱管的整體換熱性能進(jìn)行了數(shù)值模擬，得到了傳熱管各個(gè)局部的表面換熱系數(shù)和溫度分布曲線，并利用它們推導(dǎo)出了傳熱管總換熱系數(shù)和熱流密度的分布曲線。根據(jù)海水和外翅片管上的溫度分布曲線可預(yù)測(cè)傳熱管外表面易結(jié)冰的位置(如本例中的外翅片管1m以下的位置)。傳熱管總換熱系數(shù)和熱流密度的分布曲線則可為傳熱管的整體換熱性能描述提供幫助。該模型及相關(guān)模擬分析可望為該類氣化器的設(shè)計(jì)、選型和運(yùn)行管理提供參考。
參考文獻(xiàn)
[1]MORIMOTO N，YAMAMOTO S，YAMASAKI Y，et al．Development and practical application of a high performance open-rack LNG vaporizer(SuperORV)[R]．Tokyo：International Gas Union，2003．
[2]陳永東，陳學(xué)東．LNG成套裝置換熱器關(guān)鍵技術(shù)分析[J]．天然氣工業(yè)，2010，30(1)：96-100．
CHEN Yongdong，CHEN Xuedong．A technical analysis of heat exchangers in LNG plants and terminals[J]．Natural Gas Industry，2010，30(1)：96-100．
[3]張韶．粵東LNG汽化器ORV的方案選擇及技術(shù)要素的分析[J]．中國(guó)新技術(shù)新產(chǎn)品，2010，24(1)：4-5．
ZHANG Shao．A technical analysis of heal exchangers in LNG plants and terminals[J]．China New Technologies and Products，2010，24(1)：4-5．
[4]杜光能．LNG終端接收站工藝及設(shè)備[J]．天然氣工業(yè)，1999，19(5)：82-86．
DU Guangneng．Technology and equipment for LNG receiving terminal station[J]．Natural Gas Industrv，l999，19(5)：82-86．
[5]呂俊，王蕾．浙江LNG接收站項(xiàng)目氣化器選型及系統(tǒng)優(yōu)化[J]．天然氣工業(yè)，2008，28(2)：132-135．
LU Jun，WANG Lei．Selection of vaporizer types and optimization of vaporizer system in LNG receiving terminal project of Zhejiang Province[J]．Natural Gas Industry，2008，28(2)：132-135．
[6]張勇，常宏崗，艾志久，等．冷板式換熱反應(yīng)器結(jié)構(gòu)設(shè)計(jì)[J]．石油與天然氣化工，2011，41(2)：ll6-118．
ZHANG Yong，CHANG Honggang，AI Zhijiu，et al．Cold Plate Heat Transfer Reactor Structure Design[J]．Chemical Engineering of Oil&Gas，20ll，41(2)：ll6-ll8．
[7]陳永東．大型LNG汽化器的選材和結(jié)構(gòu)研究[J]．壓力容器，2007，24(11)：40-47．
CHEN Yongdong．Material choice and structure research of large LNG vaporizer[J]．Pressure Vessel Technology．2007，24(11)：40-47．
[8]邢云，劉淼兒．中國(guó)液化天然氣產(chǎn)業(yè)現(xiàn)狀及前景分析[J]．天然氣工業(yè)，2009，29(1)：l20-123．
XING Yun，LIU Miao'er．Status quo and prospect analysis on LNG industry in China[J]．Natural Gas Industry．2009．29(1)：120-123．
[9]LIN Wensheng，ZHANG Na，GU Anzhong．LNG(1iquefied natural gas)：A necessary part in China's ruturc energy infrastructure[J]．Energy，2010，35(11)：4383-4391．
[10]徐濟(jì)望．沸騰傳熱和氣液兩相流[M]．北京：原子能出版社，2001：273-319．
XU Jijun．Boiling heat transfer and gels liquid two phase flow[M]．Beijing：Atomic Energy Press，2001：273-319．
[11]CAREY V P．Liquid vapor phase change phenomena：An introduction to the thermophysics of vaporization and condensation processes in heat transfer equipmenl[M]．Washington：Taylor＆，Francis，2007．
[12]顧安忠，魯學(xué)生，汪榮順，等．液化滅然氣技術(shù)[M]．北京：機(jī)械工業(yè)出版社，2004．
GU Anzhong，LU Xuesheng，WANG Rongshun,et al．Liquefied natural gas(LNG)technology[M]．Beijing：China Machine Press，2004．
[13]CAO E．Heat transfer in process engineering[M]．Chicago：McGraw Hill Companies，Incorporated，2009．
本文作者：王萌 金滔 湯珂 陳國(guó)邦
作者單位：浙江大學(xué)制冷與低溫研究所
您可以選擇一種方式贊助本站
支付寶轉(zhuǎn)賬贊助
微信轉(zhuǎn)賬贊助
