趙思越，克紅娟，林文勝*，徐世龍

（1-上海交通大學(xué)制冷與低溫工程研究所，上海 200240；2-中國(guó)石化銷售有限公司河南分公司，河南鄭州 450016）

基于L-CNG加氣站冷能利用的蓄冰槽性能模擬與分析

趙思越1，克紅娟2，林文勝*1，徐世龍1

（1-上海交通大學(xué)制冷與低溫工程研究所，上海 200240；2-中國(guó)石化銷售有限公司河南分公司，河南鄭州 450016）

本文根據(jù)L-CNG加氣站的特點(diǎn)提出了基于L-CNG加氣站冷能利用的蓄冷系統(tǒng)，并根據(jù)蓄冰槽管外結(jié)冰的特點(diǎn)提出了蓄冰槽管外結(jié)冰模型并對(duì)模型進(jìn)行簡(jiǎn)化。本文結(jié)合影響蓄冰槽管外結(jié)冰的因素對(duì)管外結(jié)冰特性進(jìn)行了分析計(jì)算，并研究了不同的載冷劑進(jìn)口溫度、流速、濃度以及蓄冰槽內(nèi)水的初溫等因素對(duì)蓄冰槽蓄冷特性的影響。根據(jù)計(jì)算分析結(jié)果，本文提出了適合 L-CNG冷能利用系統(tǒng)的蓄冰槽最佳運(yùn)行工況，為L(zhǎng)-CNG冷能利用提供理論基礎(chǔ)。

L-CNG；冷能利用；蓄冰槽；結(jié)冰；載冷劑

0 引言

L-CNG加氣站結(jié)合了液化天然氣成本低、燃燒熱值高、環(huán)保安全、運(yùn)輸方便、靈活建站[1]和壓縮天然氣可車載氣瓶裝載的優(yōu)點(diǎn)[2]，并具有加壓效率高、占地面積小、運(yùn)行費(fèi)用低等特點(diǎn)[3]。LNG在LCNG加氣站氣化器氣化放出830kJ/kg的冷量[4]，在常規(guī)加氣站中這部分冷量通常被舍棄，造成了很大的能源浪費(fèi)。隨著 L-CNG加氣站數(shù)量逐漸的增加，以及人們對(duì)于提高能源利用效率的追求，L-CNG加氣站冷能利用問(wèn)題逐漸得到了人們的重視。與傳統(tǒng)的LNG冷能利用不同，L-CNG的冷能利用具有冷量較小且冷能不連續(xù)的特點(diǎn)，可以通過(guò)蓄冷裝置將 L-CNG氣化放出的冷能收集起來(lái)，供給 L-CNG加氣站周邊辦公、住宅等使用，這樣可以提高 L-CNG加氣站的能源利用效率，也可緩解夏季高峰期用電壓力，起到節(jié)約能源的作用[5]。

蓄冷系統(tǒng)可以分為水蓄冷系統(tǒng)、冰蓄冷系統(tǒng)和相變蓄冷3種[6]，其中冰蓄冷系統(tǒng)由于其蓄冷密度大、設(shè)備體積小以及成本較低的特點(diǎn)使其得到最為廣泛的應(yīng)用。在冰蓄冷系統(tǒng)中目前使用較多的形式為盤(pán)管式蓄冰槽[7]。相對(duì)于內(nèi)融冰蓄冰槽，外融冰式蓄冰槽可以提高蓄冷的經(jīng)濟(jì)性[8]，因此本文對(duì)采用載冷劑間接蓄冰的盤(pán)管式外融冰蓄冰槽的蓄冰過(guò)程進(jìn)行蓄冰的數(shù)值模擬和分析，了解蓄冰槽傳熱的動(dòng)態(tài)過(guò)程。

1 L-CNG冷能利用系統(tǒng)

將 L-CNG冷能利用應(yīng)用于冰蓄冷空調(diào)系統(tǒng)主要流程主要分為3個(gè)部分：

1）將儲(chǔ)罐的 LNG經(jīng)低溫高壓泵加壓達(dá)到25 MPa，使其滿足CNG對(duì)于壓力的要求，然后LNG經(jīng)過(guò)氣化器氣化為 CNG同時(shí)與載冷劑換熱，將來(lái)自LNG的冷量傳遞給氣化器另一端的載冷劑；

2）選用適當(dāng)?shù)妮d冷劑，通過(guò)載冷劑循環(huán)將從氣化器得到的冷量傳遞給蓄冰槽內(nèi)的水；

3）蓄冰槽內(nèi)水接收來(lái)自載冷劑的冷量后結(jié)冰，完成蓄冷過(guò)程的同時(shí)也要保證隨時(shí)將冷量供給下游空調(diào)系統(tǒng)。

基于 L-CNG冷能利用的蓄冷系統(tǒng)在冷源的提供上與傳統(tǒng)的冰蓄冷空調(diào)系統(tǒng)有很大的不同，前者依靠液化天然氣氣化提供冷量，后者利用制冷機(jī)組在用電低谷時(shí)期蓄冷，起到“移峰填谷”的作用[9]，因此兩者最佳結(jié)冰工況存在較大差異。

2 蓄冷模型的建立及求解

2.1 蓄冷模型建立

蓄冰槽管外結(jié)冰的過(guò)程是一個(gè)三維非穩(wěn)態(tài)的過(guò)程，結(jié)冰過(guò)程中冰層的厚度不斷發(fā)生變化，蓄冰槽內(nèi)水的溫度不斷變化，是一個(gè)非常復(fù)雜的過(guò)程。載冷劑在管內(nèi)流動(dòng)，與管外的水或者冰換熱主要經(jīng)歷以下3個(gè)過(guò)程：

1）載冷劑在管內(nèi)流動(dòng)與盤(pán)管內(nèi)壁的強(qiáng)制對(duì)流換熱；

2）通過(guò)盤(pán)管壁的導(dǎo)熱；

3）盤(pán)管外壁與冰的導(dǎo)熱以及水的自然對(duì)流換熱。

為了簡(jiǎn)化模型，對(duì)結(jié)冰過(guò)程做出如下假設(shè)[10]：

1）盤(pán)管的材料均一，各處的導(dǎo)熱系數(shù)相等；

2）冰的物理性質(zhì)視為常數(shù)，即冰的密度、導(dǎo)熱系數(shù)、比熱容、結(jié)冰潛熱視為常數(shù)，不隨溫度變化；

3）盤(pán)管在管外結(jié)冰可視為與盤(pán)管同心的圓環(huán)，可以按照純導(dǎo)熱計(jì)算其熱阻；

4）忽略軸向?qū)?，將其化為徑向的一維流動(dòng)問(wèn)題；

5）認(rèn)為固液界面的溫度保持不變，為0 ℃；

6）忽略結(jié)冰過(guò)程中水的顯熱，認(rèn)為開(kāi)始結(jié)冰后水的溫度保持0 ℃。

通過(guò)簡(jiǎn)化，可以將復(fù)雜三維結(jié)冰過(guò)程模型簡(jiǎn)化為一維傳熱模型。

如圖1所示，將蓄冰槽盤(pán)管沿載冷劑流動(dòng)方向分為 N段等長(zhǎng)的微元段，對(duì)每一小段進(jìn)行能量分析。當(dāng)分段足夠小時(shí)，可以認(rèn)為該段載冷劑進(jìn)口溫度Tin,i與該微元段載冷劑出口溫度Tout,i相差很小，此段載冷劑溫度Ts,i≈Tin,i。對(duì)于第i微元段來(lái)說(shuō)，該段的入口溫度等于上一微元段的出口溫度，即

蓄冷過(guò)程以是否結(jié)冰為界限，分為水的顯熱蓄冷和冰的潛熱蓄冷兩部分。

圖1 蓄冰槽分段示意圖

2.2 水的顯熱蓄冷

顯熱蓄冷階段載冷劑和水的換熱量等于載冷劑的吸熱量也等于水的放熱量，即：

Tin,i、Tout,i——第i段內(nèi)乙二醇溶液的入口溫度和出口溫度，℃；

qm——乙二醇溶液的質(zhì)量流量 kg/s；

Cp,s——乙二醇溶液的定壓熱容J/(kg·℃)；

kAi——傳熱系數(shù)與傳熱面積的乘積，有：

式中：

Rs,i——乙二醇對(duì)流換熱熱阻，℃/W，其中乙二醇管內(nèi)強(qiáng)制對(duì)流換熱系數(shù)由Dittus-Boelter公式[11]求得：

Rtube,i——管壁熱阻，℃/W；

Rwater,i——水對(duì)流換熱熱阻，℃/W，其中水自然對(duì)流換熱系數(shù)由Churchill and Chu[12]經(jīng)驗(yàn)式求得：

同樣，在水側(cè)有：

式中：

mwater.i——第i段水的質(zhì)量，kg；

Cp,water——水的比熱容，J/(kg·℃)；——上一時(shí)刻第i段水槽內(nèi)水的溫度和本時(shí)刻內(nèi)水槽的溫度，K。

2.3 冰的潛熱蓄冷

在潛熱蓄冷階段水和載冷劑的換熱量等于載冷劑的吸熱量也等于結(jié)冰的放熱量，即：

式中：

hice——冰的熔解熱，取3.35×105J/kg；

ρice——冰的密度，取900 kg/m3。

式中，Rs,i，Rtube,i，Rwater,i與顯熱蓄冷中相同，Rice,i是冰層的導(dǎo)熱熱阻。

由于冰層厚度 dice未知，可以通過(guò)迭代的方法求得冰層厚度，再由冰層厚度求得載冷劑出口溫度。由于采用的是外融冰式蓄冰槽，在其蓄冰率（IPF）達(dá)到40%時(shí)停止蓄冰。

3 蓄冷模擬結(jié)果及分析

3.1 載冷劑的選擇

在常規(guī)的蓄冷裝置中，載冷劑一般選擇質(zhì)量分?jǐn)?shù)為 25%的乙二醇溶液，其凝固點(diǎn)為-10.7 ℃。在基于L-CNG冷能利用的蓄冰槽設(shè)計(jì)中，由于LNG的氣化溫度低，因此為了防止載冷劑被凍結(jié)，應(yīng)該選擇凝固點(diǎn)更低的載冷劑，同時(shí)還應(yīng)該滿足載冷劑無(wú)毒、不可燃、性質(zhì)穩(wěn)定、易于獲得、價(jià)格低廉以及對(duì)環(huán)境不會(huì)產(chǎn)生污染等特點(diǎn)[13]。

由圖2可以看出，在一定濃度范圍內(nèi)，乙二醇水溶液凝固點(diǎn)隨著溶液濃度增加而降低，因此可以通過(guò)增加乙二醇溶液濃度的方法降低其凝固點(diǎn)，使載冷劑滿足要求。

圖2 乙二醇溶液濃度-冰點(diǎn)曲線圖

圖3給出了不同質(zhì)量濃度的乙二醇溶液粘度的變化[14]。乙二醇溶液的濃度越高，在同一溫度下粘度越高，越不利于輸送，而且濃度越高會(huì)對(duì)管道產(chǎn)生一定腐蝕。雖然乙二醇溶液的凝固點(diǎn)可以達(dá)到-50 ℃，但由于低溫下其粘度上升迅速，其工業(yè)應(yīng)用價(jià)值在-20 ℃以上，因此在選擇乙二醇溶液濃度是應(yīng)該綜合其凝固點(diǎn)、粘度和腐蝕性進(jìn)行考慮[15]。在下面章節(jié)的分析中，除了第3.6節(jié)的分析計(jì)算外，其余小節(jié)的分析計(jì)算中選擇的乙二醇溶液質(zhì)量分?jǐn)?shù)均為 60%，其凝固點(diǎn)為-48.3 ℃；選擇此濃度是因?yàn)檩d冷劑的最低溫度為-25 ℃，60%濃度乙二醇溶液在此溫度下不會(huì)結(jié)冰，且粘度較低，滿足對(duì)載冷劑的要求。

3.2 蓄冷結(jié)果分析

將8 m長(zhǎng)的蓄冰管段等分為80份，對(duì)每一微元段進(jìn)行能量分析。取壁面溫度為 4 ℃，時(shí)間間隔為30 s，在這一時(shí)間間隔里，認(rèn)為其蓄冰槽內(nèi)每一微元段的溫度保持恒定且每一微元段結(jié)冰厚度相同。

利用Matlab編程進(jìn)行計(jì)算，在計(jì)算中由于載冷劑的溫度和水的溫度不斷變化，其物性也不斷變化，由表查得不同濃度不同溫度下乙二醇物性后進(jìn)行擬合并載入分析，可以得到比較精確的解。

1）載冷劑出口溫度的變化

從圖 4中可以看出，當(dāng)盤(pán)管進(jìn)口為-20 ℃的60%乙二醇溶液時(shí)，出口溫度的最高值發(fā)生在初始時(shí)刻，因?yàn)榇藭r(shí)載冷劑與水的換熱溫差最大，且此時(shí)沒(méi)有結(jié)冰換熱熱阻最小，換熱量最大，此后隨著水和載冷劑的換熱溫差的減小，載冷劑進(jìn)出口溫度逐漸降低。在11 min時(shí)，出口溫度達(dá)到最低值，這是因?yàn)殡S著管外水溫度的下降，在11 min時(shí)，水的溫度為4 ℃左右，水與壁面溫度的溫差最小，水的自然換熱非常弱，因此導(dǎo)致熱阻增加。之后水的溫度逐漸降低，自然對(duì)流換熱溫差又有所增大，因此載冷劑出口溫度略有回升。在18 min左右，盤(pán)管外開(kāi)始結(jié)冰；結(jié)冰初期，由于盤(pán)管為光管，換熱熱阻較小，導(dǎo)致載冷劑出口溫度會(huì)有短暫的升高；隨著冰層厚度的逐漸增加，冰層的導(dǎo)熱熱阻成為主要熱阻，導(dǎo)致載冷劑進(jìn)出口溫差逐漸減小，但變化逐漸趨向平緩。

圖3 不同濃度乙二醇溶液粘度

圖4 載冷劑出口溫度隨時(shí)間變化

2）蓄冰槽內(nèi)水的溫度的變化

蓄冰槽內(nèi)水的初溫為 15 ℃時(shí)，且進(jìn)口載冷劑為-20 ℃的60%乙二醇溶液時(shí)蓄冰槽內(nèi)水的平均溫度的變化見(jiàn)圖5。從圖中可以看到在18 min時(shí)蓄冰槽溫度達(dá)到0 ℃；隨著時(shí)間的增加，蓄冰槽內(nèi)水溫逐漸降低，但開(kāi)始階段降低比較快，主要是因?yàn)轱@熱蓄冷的初始階段水和載冷劑的溫差較大，換熱比較強(qiáng)；11 min左右時(shí)水溫降低的速度有所減慢，隨后降溫速度又有所提升。主要是因?yàn)?1 min時(shí)槽內(nèi)水的溫度為4 ℃左右，水和壁面的換熱溫差較小導(dǎo)致自然對(duì)流的換熱量減小，水的溫度下降比較慢。

圖5 蓄冰槽內(nèi)平均水溫隨時(shí)間變化

3）冰層厚度隨時(shí)間的變化

當(dāng)載冷劑進(jìn)口溫度為-20 ℃，蓄冰槽內(nèi)水的溫度為15 ℃時(shí)，在18 min時(shí)開(kāi)始結(jié)冰，由圖6可以看到，結(jié)冰速度先逐漸增加后趨于平緩，這是由于剛開(kāi)始結(jié)冰時(shí)，盤(pán)管為光管，熱阻較小，載冷劑與水換熱強(qiáng)烈，冰層增長(zhǎng)率較快，結(jié)冰過(guò)程中冰層厚度逐漸增加，冰層的導(dǎo)熱熱阻逐漸增加，并且成為主要的換熱熱阻，換熱逐漸減小，冰層厚度的增加逐漸減緩。當(dāng)結(jié)冰厚度為 20 mm（即 IPF≈40%左右時(shí)）結(jié)冰停止。

3.3 載冷劑進(jìn)口溫度的影響

改變載冷劑進(jìn)口溫度，取-25 ℃、-20 ℃、-15 ℃、-10 ℃和-5 ℃，觀察載冷劑進(jìn)口溫度對(duì)于結(jié)冰的影響，如圖7～圖9所示。

由圖7可以看出，載冷劑溫度變化時(shí)，蓄冰槽內(nèi)水溫變化趨勢(shì)基本相同。載冷劑進(jìn)口溫度越低，水溫下降得越快，水溫降到0 ℃所需的時(shí)間越短，載冷劑進(jìn)口溫度為-25 ℃、-20 ℃、-15 ℃、-10 ℃和-5 ℃時(shí)蓄冰槽內(nèi)水下降到 0 ℃需要的時(shí)間分別為為16 min、19 min、24 min、31 min和48 min。

由圖8可以看出，載冷劑進(jìn)口溫度對(duì)于結(jié)冰的影響。載冷劑進(jìn)口溫度越低，結(jié)冰厚度達(dá)到20 mm所需要的時(shí)間越短，并且隨著載冷劑溫度的降低，結(jié)冰初期冰層厚度的增長(zhǎng)率越快。

圖6 結(jié)冰厚度隨時(shí)間的變化

圖7 載冷劑溫度對(duì)蓄冰槽內(nèi)水的平均溫度的影響

圖8 載冷劑溫度對(duì)于結(jié)冰的影響

由圖 9可以看到，載冷劑溫度由-5℃降低到-25 ℃時(shí)，結(jié)冰時(shí)間變?yōu)樵瓉?lái)的20%，特別是載冷劑溫度低于-10 ℃時(shí)對(duì)于結(jié)冰時(shí)間有很大的減小。傳統(tǒng)的冰蓄冷空調(diào)需要制冷主機(jī)提供冷源，載冷劑進(jìn)口溫度的降低可以減少結(jié)冰時(shí)間，但是過(guò)低的載冷劑溫度會(huì)降低制冷主機(jī)的效率，一般制冷主機(jī)的蒸發(fā)溫度每下降 1 ℃，主機(jī)的效率下降 3%[16]。因此設(shè)計(jì)時(shí)應(yīng)該綜合考慮，不宜選取過(guò)低的載冷劑溫度，一般溫度選擇在-5 ℃左右，此時(shí)制冷機(jī)組COP下降比較少且結(jié)冰量比較多[17]?；?L-CNG冷能利用的蓄冷系統(tǒng)的冷源與傳統(tǒng)冰蓄冷空調(diào)有較大的不同。對(duì)于基于 L-CNG加氣站冷能利用的蓄冰槽的設(shè)計(jì)，LNG的最低溫度為-162 ℃，選取合適的載冷劑，合理設(shè)計(jì) LNG氣化器，完全可以為蓄冰槽提供-20 ℃甚至更低溫度的冷量。但是如果蓄冷溫度過(guò)低，LNG冷能利用率將下降，對(duì)于氣化器盤(pán)管材料也會(huì)提出更高的要求，會(huì)導(dǎo)致制作成本的增加，因此綜合考慮，-15 ℃～-20 ℃是較為合適的載冷劑溫度區(qū)間。

圖9 不同載冷劑進(jìn)口溫度蓄冷時(shí)間比較

3.4 蓄冰槽內(nèi)水的初溫的影響

由圖10和圖11可以看到，不同的蓄冰槽水的初溫下，蓄冰槽內(nèi)載冷劑的出口溫度以及水的溫度變化趨勢(shì)基本是相同的。由圖 10可以看到，在同樣的載冷劑進(jìn)口溫度-20 ℃下，蓄冰槽內(nèi)水的初溫越高，載冷劑與水的溫差越大，換熱量越大，載冷劑初始時(shí)刻出口溫度越高。

由圖11可以看到，蓄冰槽內(nèi)水的溫度越高，蓄冰槽內(nèi)水降到 0 ℃所需要的時(shí)間越長(zhǎng)，而且蓄冰槽內(nèi)水的初溫越高，水的溫度變化率越高，這也是因?yàn)槌鯗卦礁咚洼d冷劑的溫差越大，單位時(shí)間內(nèi)水溫下降的越多。但是蓄冰槽的初始水溫的高低對(duì)于結(jié)冰并沒(méi)有影響，因?yàn)榻Y(jié)冰時(shí)水已經(jīng)全部降低到0 ℃，已經(jīng)和水的初溫沒(méi)有關(guān)系。因此，在不同的蓄冰槽初始溫度下，蓄冰槽內(nèi)水的結(jié)冰時(shí)間和結(jié)冰厚度以及結(jié)冰過(guò)程中載冷劑出口溫度變化是相同的。

由圖10和圖11還可以看到，當(dāng)蓄冰槽內(nèi)水的初溫低于4 ℃時(shí)，載冷劑出口溫度和水槽平均溫度與其他溫度有所不同，載冷劑出口溫度不會(huì)出現(xiàn)進(jìn)出口溫差極小值，這是因?yàn)榈陀? ℃時(shí)不會(huì)出現(xiàn)水和壁面自然對(duì)流換熱的極小值。

圖10 蓄冰槽內(nèi)水的初溫對(duì)載冷劑出口溫度影響

圖11 蓄冰槽內(nèi)水的初溫對(duì)槽內(nèi)水平均溫度的影響

3.5 載冷劑流速的影響

改變載冷劑的流速，取 0.3 m/s～2 m/s，圖 12顯示了流速對(duì)于結(jié)冰時(shí)間的影響；流速在1 m/s～2 m/s時(shí)，結(jié)冰時(shí)間幾乎沒(méi)有太大的變化，而當(dāng)流速小于1m/s時(shí)結(jié)冰時(shí)間明顯增加。

圖 13為載冷劑流速對(duì)于載冷劑溫度的影響?？梢钥吹剑淖冚d冷劑的流速對(duì)于初始的載冷劑出口溫度的影響較大，主要是因?yàn)樵趽Q熱初期換熱不均勻，載冷劑流速對(duì)于換熱的影響較大；因此，流速降低時(shí)，載冷劑和水有充分的時(shí)間換熱，載冷劑出口溫差提高。隨著換熱的進(jìn)行，特別是隨著結(jié)冰的進(jìn)行，影響換熱的主要因素為冰的導(dǎo)熱熱阻；在結(jié)冰結(jié)束時(shí)載冷劑的流速對(duì)于載冷劑出口溫度影響不大。

圖12 不同載冷劑流速對(duì)結(jié)冰時(shí)間的影響

圖13 不同載冷劑流速對(duì)載冷劑出口溫度的影響

3.6 載冷劑濃度的影響

選取質(zhì)量分?jǐn)?shù)為30%、40%和60%的乙二醇溶液，分析在載冷劑進(jìn)口溫度為-10 ℃時(shí)的結(jié)冰情況。圖14和圖15為不同乙二醇濃度載冷劑出口溫度和結(jié)冰厚度隨時(shí)間的變化圖。由圖中可以看到，對(duì)于不同濃度相同溫度的乙二醇溶液，濃度對(duì)于結(jié)冰的影響很小，因此只需要考慮濃度滿足其凝固溫度的影響即可，乙二醇濃度過(guò)高反而會(huì)腐蝕管道，對(duì)傳熱造成影響。

圖14 不同濃度乙二醇溶液下載冷劑出口溫度變化

圖15 不同濃度乙二醇溶液對(duì)結(jié)冰厚度的影響

4 總結(jié)與展望

本文對(duì)基于 L-CNG冷能利用的蓄冰槽管外結(jié)冰的過(guò)程進(jìn)行了分析，通過(guò)數(shù)值分析可以觀察到載冷劑初溫、載冷劑流速以及槽內(nèi)水的初溫以及載冷劑的種類對(duì)于結(jié)冰的影響，得出如下結(jié)論：

1）降低載冷劑的溫度可以提高蓄冷的效率，減少結(jié)冰時(shí)間，但載冷劑溫度不宜過(guò)低，因?yàn)闇囟冗^(guò)低后對(duì)于結(jié)冰效率的提升并不明顯，反而會(huì)增加氣化器的材料成本，降低了其經(jīng)濟(jì)性；

2）蓄冰槽內(nèi)水的初溫對(duì)于顯熱蓄冷階段的影響較大，對(duì)于潛熱蓄冷階段沒(méi)有很大的影響；

3）載冷劑的流速小于0.5 m/s時(shí)，結(jié)冰時(shí)間明顯增加且載冷劑進(jìn)出口溫差也有明顯升高，而載冷劑流速在1 m/s～2 m/s時(shí)對(duì)于結(jié)冰的影響不大；可以在此流速區(qū)間內(nèi)選擇載冷劑的流速；

4）乙二醇溶液可以作為 L-CNG冷能利用系統(tǒng)的中間流體。不同濃度的乙二醇對(duì)結(jié)冰的影響不大，只需要根據(jù)其凝固點(diǎn)選擇合適濃度乙二醇即可。

[1]楊文剛, 陳杰. 無(wú)動(dòng)力 LNG汽車加注裝置[J]. 制冷技術(shù), 2013, 33(3): 27-29.

[2]錢文斌, 章陽(yáng), 李丹華, 等. L-CNG與 CNG加氣站技術(shù)經(jīng)濟(jì)比較[J]. 煤氣與熱力, 2012, 32(4): B14-B16.

[3]祝家新, 林文勝. 天然氣汽車加氣站發(fā)展趨勢(shì)及L-CNG加氣站技術(shù)特點(diǎn)[J]. 能源技術(shù), 2007, 28(1):32-35.

[4]顧安忠, 魯雪生, 汪順榮, 等. 液化天然氣技術(shù)[M]. 2版. 北京: 機(jī)械工業(yè)出版社, 2015.

[5] 王玉偉, 后磊. L-CNG 加氣站及其冷能回收用于冰蓄冷系統(tǒng)的設(shè)計(jì)[J]. 石油庫(kù)與加油站, 2015, 24(12):14-17.

[6]李浙, 黃銳遠(yuǎn). 冰蓄冷空調(diào)系統(tǒng)設(shè)計(jì)方法[J]. 制冷技術(shù), 2001, 21(2): 11-15.

[7]方貴銀. 蓄能空調(diào)技術(shù)[M]. 北京: 機(jī)械工業(yè)出版社,2006: 80-120.

[8]朱學(xué)錦. 某金融中心冰蓄冷系統(tǒng)的應(yīng)用與經(jīng)濟(jì)性分析[J]. 制冷技術(shù), 2016, 36(4): 72-78.

[9]盧建波. 論冰蓄冷空調(diào)系統(tǒng)的優(yōu)劣[J]. 建筑安全,2013(5): 39-41.

[10]鄒同華, 申江, 成天及. 水平管外結(jié)冰特性的理論與實(shí)踐研究[J]. 制冷技術(shù), 2001, 21(3): 32-35.

[11]楊世銘, 陶文銓. 傳熱學(xué)[M]. 北京: 高等教育出版社,2007: 243-255.

[12]呂曉艷. 管殼式外融冰槽熱工性能及應(yīng)用研究[D]. 北京: 清華大學(xué), 2004.

[13]于航. 空調(diào)蓄冷技術(shù)與設(shè)計(jì)[M]. 北京: 機(jī)械工業(yè)出版社, 2007: 50-123.

[14]ASHRAE. ASHRAE handbook 2009 Fundamentals[M].Atlanta: American of Society of Heating, Refrigerating and Air-conditioning Engineers, Inc., 2009.

[15]李曼, 張哲, 田津津. LNG冷能用于冷藏車的系統(tǒng)設(shè)計(jì)及分析[J]. 制冷技術(shù), 2014, 34(4): 38-52.

[16]陳林. 不完全凍結(jié)方式的導(dǎo)熱塑料蓄冰盤(pán)管傳熱性能研究[D]. 杭州: 浙江大學(xué), 2002.

[17]張奕, 張小松. 冷媒進(jìn)出口溫度對(duì)管外結(jié)冰過(guò)程的影響[J]. 暖通空調(diào), 2004, 34(10): 110-113.

Simulation and Analysis of Performance of Ice Storage Tank Based on Cold Energy Utilization in L-CNG Fueling Station

ZHAO Siyue1, KE Hongjuan2, LIN Wensheng*1, XU Shilong1
(1- Institute of Refrigeration and Cryogenics, Shanghai Jiao Tong University, Shanghai 200240, China;2-Henan Petroleum Branch, SinoPec Sales Co. Ltd., Zhengzhou, Henan 450016, China)

In this paper, the ice storage system based on the cold energy utilization in L-CNG fueling station has been proposed. On the basis of the characteristics of the ice-forming on the outside of the coils, the model is proposed and simplified, and the analysis and calculation to the icing process under the influence of different factors are performed. Meanwhile the influence factors, such as the temperature, the flow rate and the concentration of the secondary refrigerant and so on, are also studied. The optimal operational states of the ice storage tank are suggested, which can provide theoretical basis for the study of the cold energy utilization in L-CNG fueling station.

L-CNG; Cold energy utilization; Ice storage tank; Ice-forming; Secondary refrigerant

10.3969/j.issn.2095-4468.2017.04.200

*林文勝（1967-），男，博士，副教授。研究方向：液化天然氣。聯(lián)系地址：上海市閔行區(qū)東川路800號(hào)上海交通大學(xué)，郵編：200240。聯(lián)系電話：021-34206533。E-mail：linwsh@sjtu.edu.cn。

本論文選自2016年第九屆全國(guó)制冷空調(diào)新技術(shù)研討會(huì)。