1 概述

LNG儲罐的保冷性能直接關系到儲罐的安全運營，靜態(tài)蒸發(fā)率是衡量儲罐保冷性能的重要指標

。文獻[2]建立了預測LNG儲罐靜態(tài)蒸發(fā)率和罐內壓力變化的數(shù)學模型，分析了靜態(tài)蒸發(fā)率的影響因素。文獻[3]給出了多種LNG儲罐靜態(tài)蒸發(fā)率的測量和計算方法，并且詳細論述了各種測量方法的適用情況。以往文獻對大型LNG儲罐保冷性能的研究主要集中在熱流量和靜態(tài)蒸發(fā)率計算方法的分析上，對儲罐熱流量影響因素的研究較少。本文以某3×10

m

雙金屬LNG儲罐為例，給出了詳細的熱流量和靜態(tài)蒸發(fā)率計算方法，并計算分析了環(huán)境溫度、環(huán)境風速、充滿率對儲罐熱流量的影響。

2 雙金屬LNG儲罐基本結構

雙金屬LNG儲罐的主體結構通常由不銹鋼外罐、不銹鋼內罐、熱角保護層和保冷層等組成。本文以濟南某雙金屬LNG儲罐為例進行計算分析，圖1和圖2為該LNG儲罐的保冷結構，內外罐體材質均為S30408不銹鋼，內外罐之間填充保冷材料。本文約定，外罐罐壁板、內罐罐壁板以及它們之間的夾層空間為罐壁，罐壁的保冷主要依靠膨脹珍珠巖和彈性玻璃纖維氈的雙層組合保冷材料；在外罐罐壁板內側底部設有一定厚度的泡沫玻璃磚作為熱角保護層；內罐罐底板以下的儲罐部分為罐底，罐底的保冷材料主要為泡沫玻璃磚；穹頂、吊頂板以及它們之間的玻璃棉和BOG儲存空間為罐頂，吊頂板和玻璃棉合稱為吊頂，罐頂?shù)谋＠渲饕揽康蹴數(shù)牟Ａ薇＠鋵印?/p>

3 儲罐熱流量和靜態(tài)蒸發(fā)率的計算

根據(jù)LNG儲罐的結構特點，分析建立儲罐熱流量和靜態(tài)蒸發(fā)率計算的數(shù)學模型。為了簡化計算，對模型做出如下5個假設：儲罐各層材料之間不存在空隙，忽略鄰層之間的接觸熱阻。罐內處于氣液平衡狀態(tài)，罐內溫度分布均勻且相等。將儲罐與外界環(huán)境的傳熱過程簡化為一維穩(wěn)態(tài)傳熱過程。 LNG的組成為純甲烷。大氣壓按101.325 kPa考慮。

3.1 無太陽輻射的罐壁熱流量計算

① 儲罐罐壁熱流量的計算

在無太陽輻射時，罐壁的熱流量主要是空氣傳遞給儲罐的熱流量。由于LNG與BOG對內罐罐壁板內表面的對流傳熱強度不同，需要分別考慮BOG和LNG對應罐壁的熱流量。吊頂以上罐壁內外溫差很小，忽略熱流量。根據(jù)圓筒壁的熱流量計算公式

可推導得到無太陽輻射時儲罐罐壁的熱流量計算公式，見式(1)、(2)。

使用直接請求策略，主要有兩種狀態(tài)，其一：直接表達會話核心，只用中心行為語。其二：使用稱呼語引起注意，然后表達核心內容，即起始行為語+中心行為語。如（1）“幫我拿那件衣服看看?！薄拔乙@個。”（2）“爸爸，給我打點錢?！薄盃敔?，您能給我講講你們時代的事嗎？”“xxx，清明一起去旅游吧？”采訪發(fā)現(xiàn)，留學生直接請求策略使用頻率相對不高，表意較委婉。

(1)

(2)

式中

——無太陽輻射時內罐LNG對應壁面的熱流量，W

桔皮書是實施專利鏈接制度的基礎，只有列入桔皮書中的專利才能適用這一制度。列入桔皮書的專利情況能夠反映原研藥企業(yè)對這一制度的利用情況，桔皮書中登記的專利的量和類型變化也反映了藥品專利的發(fā)展情況和趨勢。

沒有選沈從文的作品，主要原因恐怕是魯迅不認同沈從文的作品。但是，兩個人的矛盾影響了選本編選的可能是存在的。1935年11月，沈從文寫了《讀〈中國新文學大系〉》的書評，沒有談及自己的作品是否應該入選，但認為魯迅的選本“有抑彼揚此處”，“取舍之間不盡合理”，特別點明未選入王統(tǒng)照等幾人的作品，對沉鐘社、莽原社的評價過高，“皆與印行這套書籍的本意稍稍不合”。

——LNG液位高度，m

——環(huán)境溫度，K

——LNG儲罐運行溫度，K

——空氣對儲罐外表面的表面?zhèn)鳠嵯禂?shù)，W/(m

·K)

——外罐外直徑，m

——罐壁保冷材料的總層數(shù)

——罐壁保冷材料層數(shù)的編號

——罐壁第

層保冷材料的熱導率，W/(m·K)

——罐壁第

層保冷材料的外直徑，m

——罐壁第

層保冷材料的內直徑，m

——LNG對罐壁內表面的表面?zhèn)鳠嵯禂?shù)，W/(m

·K)

——內罐內直徑，m

二是成立國務院直屬的跨部門、跨行政區(qū)的管理協(xié)調機構，推進一體化管理。長江大保護要一體化推進已經(jīng)形成共識，建議成立國務院直屬的跨部門、跨行政區(qū)的管理協(xié)調機構。早在20世紀，部分發(fā)達國家開始建立強有力的流域管理機構，它們以流域為對象，建立流域管理機構，并通過法律賦予有效的行政管理權。如美國田納西流域管理局、法國塞納河流域管理局、保護萊茵河國際委員會等。此外，日本“琵琶湖綜合發(fā)展工程”也非常成功，均可借鑒。

——無太陽輻射時內罐BOG對應壁面的熱流量，W

——內罐BOG所占空間高度，m

——BOG對罐壁內表面的表面?zhèn)鳠嵯禂?shù)，W/(m

·K)

空氣對儲罐外表面的表面?zhèn)鳠嵯禂?shù)隨環(huán)境風速的增大而增大，可按式(3)

計算。

=6.02+3.46

（3）估計該家庭使用節(jié)水龍頭后，一年能節(jié)省多少水？（一年按365天計算，同一組中的數(shù)據(jù)以這組數(shù)據(jù)所在區(qū)間中點的值作代表．）

(3)

式中

——環(huán)境風速，m/s

② BOG對罐壁內表面的表面?zhèn)鳠嵯禂?shù)計算

內罐BOG與罐壁內表面的對流傳熱可以近似看作豎直熱面的大空間自然對流傳熱，設BOG對應罐壁內表面的溫度為

，由于

和

溫差較小，計算時定性溫度取

，特征長度取

，其實驗關聯(lián)式見式(4)

。

隨著社會的不斷發(fā)展，以及對不動產(chǎn)測繪要求的提升，給不動產(chǎn)測繪技術也提出了新的要求。GPS技術作為新興技術測繪手段，在不動產(chǎn)測繪工作中應用十分廣泛，它可以實現(xiàn)對目標區(qū)域的精準定位，一方面提高了不動產(chǎn)測繪結果的準確性，另一方面也有效提高了測繪的整體效率，縮短了測繪工作的周期。同時，GPS技術在進行定位測量中，還可以根據(jù)目標區(qū)域的實際情況開展測量分析，確保測得的數(shù)據(jù)結果具有更高的可靠性。另外，不同于其他測繪技術，GPS技術在實際的應用過程中不會受到氣候等外部因素的影響，可以實現(xiàn)全天候的測量，給不動產(chǎn)測繪工作的開展提供了巨大的便利。

(4)

式中

——內罐BOG與罐壁內表面對流傳熱過程的努塞爾數(shù)

——內罐BOG與罐壁內表面對流傳熱過程的格拉曉夫數(shù)

——內罐BOG的普朗特數(shù)

、

、

根據(jù)式(5)、(6)、(7)

進行計算。

2.3 4組患者治療前后才藤氏吞咽障礙評定結果比較 治療2周、4周后，觀察3組后一時間節(jié)點治愈率較前一時間節(jié)點均有顯著增加(均P<0.05)，同時，觀察3組治愈率較對照組、觀察1組、觀察2組相同時間節(jié)點有顯著增加(均P<0.05)；觀察1組與對照組比較，治愈率差異無統(tǒng)計學意義；觀察2組治療1周、治療2周后與對照組比較，治愈率差異無統(tǒng)計學意義，治療4周后治愈率較對照組明顯提高(P<0.05)。見表4。

(5)

(6)

(7)

(8)

(9)

式中

——內罐BOG的熱導率，W/(m·K)

在初中信息技術教學中，教師要做的不僅是教給學生辦公軟件的應用、網(wǎng)頁制造等方面的計算機知識，還要將更多的時間與精力放在學生網(wǎng)絡學習能力的培養(yǎng)上來。我們當前正處在一個信息爆炸的時代，而對于中學生來說，網(wǎng)絡無疑是他們學習知識、認識世界最有效、最便利的途徑。因此，教師要在日常教學之余教給學生一些正確的快速搜索信息、處理信息的方法，使學生能夠利用網(wǎng)絡這一途徑解決問題。這樣不僅能夠逐步改善信息技術教育的現(xiàn)狀，使其朝著正確的方向發(fā)展，還能逐步增強學生的網(wǎng)絡學習能力，促進學生綜合素質的提高。

——重力加速度，m/s

，取9.8 m/s

——內罐BOG的體脹系數(shù)，K

——BOG對應罐壁內表面的溫度，K

式中

——內罐LNG對罐底上表面對流傳熱過程的努塞爾數(shù)

——內罐BOG的密度，kg/m

——內罐BOG的黏度，Pa·s

g

——內罐BOG的比定壓熱容，J/(kg·K)

——BOG在溫度為273.15 K和壓力為0時的密度，kg/m

5.提升績效人員素質，推動績效審計不斷發(fā)展。油田企業(yè)內部審計部門應結合企業(yè)特點，一方面，積極吸納各專業(yè)技術人才加入審計隊伍，并注重其財經(jīng)審計知識、績效審計方法和技術的培訓；另一方面，盡快建立一套切實可行的培訓制度，進一步督促審計人員對新業(yè)務的學習，快速掌握績效審計的方法和技術，著力培養(yǎng)復合型人才，優(yōu)化審計隊伍的知識結構，提高審計人員素質，加快建設適應績效審計要求的審計隊伍，使之適應績效審計發(fā)展的需要，組建一支具有多學科、復合型的審計隊伍。

——與空氣接觸的承臺下表面面積，m

使用云天化鎂立硼復合肥+四全水溶復合肥的示范田與對照田比較，前期長勢明顯，中期病蟲害少、落果少，后期肥勁長、不脫肥、著色好。云天化示范田與其它品牌對照田相比，同等的氣候條件下，示范田櫻桃畝增產(chǎn)28.5kg，增產(chǎn)率達8%。這表明，使用云天化系列化肥可以顯著提高該地區(qū)櫻桃的產(chǎn)量。

外部環(huán)境通過BOG對應罐壁傳入儲罐內部的熱流量等于BOG從罐壁接收的熱流量，即：

=π

(

-

)

(10)

聯(lián)立式(4)和(10)，可求得

，然后代入式(2)可得到

。

③ LNG對罐壁內表面的表面?zhèn)鳠嵯禂?shù)計算

LNG與罐壁內表面的對流傳熱也可以近似看作豎直熱面的大空間自然對流傳熱，計算時定性溫度取

，特征長度取

，其實驗關聯(lián)式見式(11)

。

(11)

式中

——內罐LNG與罐壁內表面對流傳熱過程的努塞爾數(shù)

——內罐LNG與罐壁內表面對流傳熱過程的格拉曉夫數(shù)

聯(lián)立式(16)和(17)，可求得

，然后代入式(14)和(15)可得到

和

。

——LNG在溫度為

時的普朗特數(shù)

——LNG對應罐壁內表面的溫度，K

由于

與

差別較小，在計算中校正因子(

/

)

近似取1；

根據(jù)LNG熱導率計算的經(jīng)驗公式

進行計算；

按式(12)計算

，在

計算中LNG近似按純甲烷處理，甲烷在不同溫度下的密度可通過甲烷狀態(tài)圖

查得；LNG在溫度為

時的黏度

和比定壓熱容

L

根據(jù)文獻[12]中LNG黏度和比定壓熱容的計算方法進行計算。

(12)

式中

——LNG在

～

溫度范圍內的平均體脹系數(shù)，K

——LNG在溫度為

時的密度，kg/m

——LNG在溫度為

時的密度，kg/m

——地面溫度，K

——LNG的溫度范圍下限，K

外部環(huán)境通過LNG對應罐壁傳入儲罐內部的熱流量等于LNG從罐壁接收的熱流量，即：

=π

(

-

)

(13)

聯(lián)立式(11)和(13)，可求得

，然后代入式(1)可得到

。

3.2 罐底熱流量計算

儲罐罐底不受太陽輻射的影響。罐底的熱流量主要包括2部分：空氣以對流傳熱的形式傳遞給儲罐的熱流量，土壤以導熱的形式傳遞給儲罐的熱流量。由于承重柱與周圍環(huán)境溫差較小，忽略其對流換熱的熱流量。地面處的承重柱溫度可按地面溫度進行取值，在已知地面溫度的情況下可只對地面以上罐底部分進行分析。由平壁的熱流量計算公式

可推導得到儲罐罐底熱流量的計算公式。

(14)

(15)

式中

——與空氣接觸的承臺下表面對應罐底的熱流量，W

——LNG對罐底上表面的表面?zhèn)鳠嵯禂?shù)，W/(m

·K)

LNG儲罐內壓力通常處于微正壓狀態(tài)，罐內LNG和BOG與壓力有關的參數(shù)均按壓力為0進行取值。

根據(jù)文獻[6]

中天然氣的熱導率計算方法進行計算；

按式(8)計算

；BOG在273.15 K和壓力為0時的密度

在文獻[7]中查得；

根據(jù)文獻[8]中天然氣黏度的計算方法計算；BOG在溫度為

時的比定壓熱容

g

在文獻[9]中查得。

——罐底保冷材料的總層數(shù)

作為前明尼蘇達州參議員(1964－1976)、卡特時期副總統(tǒng)(1977－1981)，蒙代爾僅贏得明尼蘇達州和哥倫比亞特區(qū)的13張選舉人票，創(chuàng)下選舉人票最大差距紀錄——13：525。敗選后的蒙代爾離開政壇，重操舊業(yè)，加入明尼蘇達州一家律所，做起了律師的老本行?？肆诸D政府時期重返政壇，擔任美國駐日本大使(1993-1996)。之后從日本回到明尼蘇達老家，重新在律所做事，當普通公民。

——罐底保冷材料層數(shù)的編號

——罐底第

層保冷材料的厚度，m

——罐底第

層保冷材料的熱導率，W/(m·K)

10月10日，五師八十九團“棉思眠”農(nóng)業(yè)合作社來了一位成都老板，他與“棉思眠”農(nóng)業(yè)合作社社長陳建輝簽訂了2000套千層無網(wǎng)被合同，這是今年新棉下來后合作社簽訂的第一份訂單合同。

——承重柱對應罐底的熱流量，W

——LNG的溫度范圍上限，K

——承重柱總截面積，m

LNG與罐底上表面的對流傳熱可以看作水平熱面向上的大空間自然對流傳熱，計算時定性溫度取

，特征長度取0.25

，其實驗關聯(lián)式見式(16)

。

(16)

上莊鄉(xiāng)正在向全省乃至省外大力推廣本鄉(xiāng)的舍飼養(yǎng)羊技術，希望各地都能通過該模式解決羊規(guī)?；B(yǎng)殖帶來的生態(tài)平衡問題，同時加速羊只的生長速度，提高養(yǎng)殖戶的整體經(jīng)濟效益，為新野縣乃至河南省的畜牧產(chǎn)業(yè)長期可持續(xù)穩(wěn)定發(fā)展創(chuàng)造有利條件。

——內罐LNG對罐底上表面對流傳熱過程的格拉曉夫數(shù)

外部環(huán)境通過罐底傳入儲罐內部的熱流量等于LNG通過對流傳熱從罐底接收的熱流量，即：

+

=(

+

)

(

-

)

(17)

式中

——罐底上表面的平均溫度，K

——LNG在溫度為

時的普朗特數(shù)

3.3 無太陽輻射的罐頂熱流量計算

對于雙金屬LNG儲罐，穹頂不具備保冷性能，在計算分析中可以近似認為穹頂內外表面溫度相同且內外表面積相同。

① 傳熱熱流量的關系

根據(jù)對流傳熱熱流量基本計算公式

可得無太陽輻射時空氣向穹頂傳遞的熱流量為：

=

(

-

)

(18)

式中

——無太陽輻射時空氣向穹頂傳遞的熱流量，W

——穹頂外表面面積，m

——無太陽輻射時穹頂溫度，K

同理可得穹頂向罐頂BOG傳遞的熱流量為：

=

(

-

)

(19)

式中

——無太陽輻射時穹頂向罐頂BOG傳遞的熱流量，W

——罐頂BOG對穹頂?shù)谋砻鎮(zhèn)鳠嵯禂?shù)，W/(m

·K)

——無太陽輻射時罐頂BOG溫度，K

同理可得罐頂BOG向吊頂傳遞的熱流量為：

=

(

-

)

(20)

式中

——無太陽輻射時罐頂BOG向吊頂傳遞的熱流量，W

——吊頂上表面面積，m

——罐頂BOG對吊頂上表面的表面?zhèn)鳠嵯禂?shù)，W/(m

·K)

——無太陽輻射時吊頂上表面溫度，K

根據(jù)兩表面封閉系統(tǒng)的輻射傳熱熱流量計算公式

，可得無太陽輻射時穹頂向吊頂輻射傳熱的熱流量為：

(21)

式中

——無太陽輻射時穹頂向吊頂輻射傳熱的熱流量，W

——斯忒藩-玻爾茲曼常量，W/(m

·K

)，取5.67×10

W/(m

·K

)

——穹頂內表面的發(fā)射率

——穹頂內表面對吊頂上表面的角系數(shù)

——吊頂上表面的發(fā)射率

根據(jù)平壁的熱流量計算公式得到無太陽輻射時吊頂上表面向儲罐內部傳遞的熱流量為：

(22)

式中

——無太陽輻射時吊頂上表面向儲罐內部傳遞的熱流量，W

——吊頂保冷材料的總層數(shù)

——吊頂保冷材料層數(shù)的編號

——吊頂?shù)?/p>

層保冷材料的厚度，m

——吊頂?shù)?/p>

層保冷材料的熱導率，W/(m·K)

——內罐BOG對吊頂下表面的表面?zhèn)鳠嵯禂?shù)，W/(m

·K)

根據(jù)對流傳熱熱流量基本計算公式，可得無太陽輻射時吊頂下表面向內罐BOG傳遞的熱流量為：

=

(

-

)

(23)

式中

——無太陽輻射時吊頂下表面向內罐BOG傳遞的熱流量，W

——吊頂下表面面積，m

——吊頂下表面溫度，K

② 對流傳熱系數(shù)的計算

穹頂近似按水平面處理，罐頂BOG與穹頂?shù)膶α鱾鳠峥梢钥醋魉綗崦嫦蛳碌淖匀粚α鱾鳠帷Ｓ捎隈讽敳痪哂斜＠湫阅埽?/p>

和

均接近

，計算時定性溫度取

，特征長度取

/(π

)，其實驗關聯(lián)式見式(24)

。

(24)

式中

——罐頂BOG對穹頂內表面對流傳熱過程的努塞爾數(shù)

——罐頂BOG對穹頂內表面對流傳熱過程的格拉曉夫數(shù)

——罐頂BOG在溫度為

時的普朗特數(shù)

罐頂BOG與吊頂上表面的對流傳熱可以看作為水平冷面向上的自然對流傳熱，吊頂上表面和罐頂BOG溫度雖然與環(huán)境溫度有一定溫差，但是溫差較小，計算時定性溫度近似取

，特征長度取

/4，其實驗關聯(lián)式見式(25)

。

(25)

式中

——罐頂BOG對吊頂上表面對流傳熱過程的努賽爾數(shù)

——罐頂BOG對吊頂上表面對流傳熱過程的格拉曉夫數(shù)

內罐BOG與吊頂下表面的對流傳熱可以看作為水平熱面向下的自然對流傳熱，計算時定性溫度近似取

，特征長度取

/4，其實驗關聯(lián)式見式(26)

。

(26)

式中

——罐頂BOG對吊頂上表面對流傳熱過程的努塞爾數(shù)

——罐頂BOG對吊頂上表面對流傳熱過程的格拉曉夫數(shù)

③ 吊頂向儲罐內部傳遞的熱流量求解

傳入與流出穹頂?shù)臒崃髁肯嗟?，即?/p>

=

+

(27)

穹頂向罐頂BOG傳遞的熱流量等于罐頂BOG向吊頂傳遞的熱流量，即：

=

(28)

從穹頂和罐頂BOG流入到吊頂上表面的熱流量等于從吊頂上表面經(jīng)吊頂各層流出到內罐BOG的熱流量，見式(29)。

+

=

(29)

從吊頂上表面經(jīng)吊頂各層流入到吊頂下表面的熱流量等于從吊頂下表面流出到內罐BOG的熱流量，見式(30)。

=

(30)

聯(lián)立式(24)～(30)，可求得

和

，然后代入式(22)可求得

。

3.4 有太陽輻射的罐壁和罐頂熱流量計算

在有太陽輻射時，儲罐罐壁外表面和穹頂(穹頂內、外表面溫度相同)在接收太陽輻射熱后溫度升高，增加儲罐罐壁和罐頂?shù)臒崃髁?。罐壁和穹頂外表面附近空氣溫度相比無太陽輻射時的溫升見公式(31)

。以

表示不同朝向的儲罐外表面，將儲罐罐壁外表面按8個朝向各占罐壁外表面積1/8進行處理。以東朝向罐壁面為例，正東左右各偏22.5°(共45°)的罐壁圓柱面為東朝向罐壁面(

=1)，其他方向以此類推，西、南、北、東南、東北、西南、西北朝向罐壁面分別為

=2～8，穹頂外表面近似按水平表面(

=9)進行處理。儲罐各朝向的日平均太陽輻射照度可在文獻[13]

中查得。罐外表面的太陽輻射熱吸收系數(shù)

參照GB 50176—2016《民用建筑熱工設計規(guī)范》附錄第B.5條選取。

(31)

式中 Δ

——太陽輻射導致的

朝向儲罐外表面附近空氣溫度相比無太陽輻射時的溫升，K

——儲罐不同朝向外表面的編號

——

朝向日平均太陽輻射照度，W/m

——儲罐外表面的太陽輻射熱吸收系數(shù)

儲罐罐壁外表面附近空氣的平均溫升為：

(32)

式中 Δ

——太陽輻射導致的罐壁外表面附近空氣溫度相比無太陽輻射時的平均溫升，K

有太陽輻射的罐壁熱流量計算方法與無輻射時的相同，只是計算過程中用到的環(huán)境溫度不同。太陽輻射下罐壁熱流量計算用到的環(huán)境溫度

等于

加上溫升Δ

。

同理，太陽輻射下罐頂熱流量計算用到的環(huán)境溫度

等于

加上溫升Δ

。

3.5 LNG蒸發(fā)量的計算

本文定義

為儲罐單日接收太陽輻射的時間占單日總時間(24 h)的比例，對儲罐有太陽輻射和無太陽輻射時的熱流量進行加權平均可得儲罐總熱流量的計算公式，見式(33)。

=

+

+(1-

)(

+

+

)+

(

+

+

)

(33)

式中

——儲罐的總熱流量，W

——儲罐單日接收太陽輻射的時間占單日總時間(24 h)的比例

——有太陽輻射時LNG對應罐壁的熱流量，W

——有太陽輻射時BOG對應罐壁的熱流量，W

——有太陽輻射時吊頂向儲罐內部傳遞的熱流量，W

儲罐24 h內LNG蒸發(fā)的質量可通過式(34)、(35)計算，其中LNG的氣化潛熱

可通過文獻[6]

查得。

=

(34)

式中

——24 h內環(huán)境流入儲罐的總熱量，J

——時間，s，取86 400 s

儲罐24 h內LNG蒸發(fā)的質量為：

(35)

式中

——24 h內LNG蒸發(fā)的質量，kg

——LNG的氣化潛熱，J/kg

3.6 靜態(tài)蒸發(fā)率的計算

LNG儲罐的靜態(tài)蒸發(fā)率是指LNG儲罐在額定充滿率(額定充液高度時儲存LNG的體積與LNG儲罐內罐總容積之比)下，靜置達到熱平衡后，折算成在儲罐外部環(huán)境為標準狀態(tài)(1.013 25×10

Pa，293.15 K)下儲罐24 h內自然蒸發(fā)產(chǎn)生的LNG蒸發(fā)氣(BOG)的質量與儲罐有效容積(額定充滿率時儲存LNG的體積)下LNG總質量的比例，其計算公式見式(36)。

(36)

式中

——LNG儲罐的靜態(tài)蒸發(fā)率

——

取293.15 K時儲罐24 h內LNG蒸發(fā)的質量，kg

——LNG的密度，kg/m

——LNG儲罐的有效容積，m

4 實例計算

4.1 熱流量和靜態(tài)蒸發(fā)率計算

應用上述計算方法對濟南市某3×10

m

雙金屬LNG儲罐的熱流量和靜態(tài)蒸發(fā)率進行計算分析。案例儲罐的基本參數(shù)和保冷材料的物性參數(shù)見表1、2，表2中的參數(shù)參照GB 50264—2013《工業(yè)設備及管道絕熱工程設計規(guī)范》附錄A和GB/T 13350—2017《絕熱用玻璃棉及其制品》第5.6.3條。

儲罐熱流量和靜態(tài)蒸發(fā)率計算用到的其他參數(shù)如下：

=293.15 K；

=111.15 K；

=23.6 m，

=1.9 m；濟南市全年平均風速為2.6 m/s，取

=2.6 m/s；

=0.717 4 kg/m

；

=3.24×10

Pa·s；

g

=2.202 kJ/(kg·K)；

=0.01 W/(m·K)；

=140 K時

=400 kg/m

，

=100 K時

=460 kg/m

；

=116.43×10

Pa·s；

L

=3.22 kJ/(kg·K)；

=0.48 W/(m·K)；

=293.15 K；罐底承重柱數(shù)量為144根，直徑1 m，高出地面1.5 m，

=113.09 m

；

=0.8，

=0.7；

=0.03 W/(m·K)；

=10.3×10

Pa·s；

vg

=2.202 kJ/(kg·K)；

=0.48；濟南不同朝向表面接收的日平均太陽輻射照度見表3

；假設儲罐單日接收太陽輻射的時間為12 h，即

=0.5；在溫度為111.15 K和壓力為0時，

=509.74 kJ/kg，

=450.74 kg/m

。

基于以上已知參數(shù)，計算可得儲罐各部分熱流量，見表4。利用式(33)～(36)，計算可得儲罐的靜態(tài)日蒸發(fā)率

=0.074 8%?？梢?，儲罐的靜態(tài)蒸發(fā)率較小，符合國際上對大型LNG儲罐的靜態(tài)蒸發(fā)率控制要求：罐容為(1～5)×10

m

的大型LNG儲罐的靜態(tài)蒸發(fā)率≤0.08%

。

4.2 儲罐熱流量影響因素分析

由于儲罐熱流量在有無太陽輻射工況下的差別較小，并且受外部環(huán)境條件和運行工況影響的變化趨勢相同，下面以無太陽輻射工況為例，對儲罐各部分熱流量的影響因素進行計算分析。

.

.

環(huán)境溫度的影響

儲罐在額定充滿率和環(huán)境風速為2.6 m/s情況下各部分熱流量隨環(huán)境溫度的變化見圖3。罐壁、罐頂和罐底的熱流量隨環(huán)境溫度的變化基本呈線性變化規(guī)律，其變化率分別為182.5、94.3、51.8 W/℃。儲罐各部分熱流量受環(huán)境溫度影響較大。

.

.

環(huán)境風速的影響

LNG儲罐外部的環(huán)境風速主要是通過改變儲罐外表面的表面?zhèn)鳠嵯禂?shù)來影響儲罐各部分的熱流量，儲罐在不同環(huán)境風速下的外表面表面?zhèn)鳠嵯禂?shù)可按式(3)計算。儲罐在額定充滿率和環(huán)境溫度取濟南全年平均氣溫15.15 ℃情況下各部分的熱流量隨環(huán)境風速的變化見圖4?？梢?，隨著環(huán)境風速的增大，儲罐的各部分熱流量都會增加，但是增加的幅度較小，環(huán)境風速對儲罐各部分熱流量的影響較小。

.

.

儲罐充滿率的影響

LNG儲罐充滿率主要影響罐壁的熱流量，圖5為LNG儲罐在環(huán)境風速為2.6 m/s和環(huán)境溫度為15.15 ℃時罐壁熱流量隨儲罐充滿率的變化曲線?？梢?，隨著儲罐充滿率的增加，儲罐罐壁熱流量呈線性變化規(guī)律，平均每增加10%的充滿率，BOG對應罐壁的熱流量減少3 278.7 W，LNG對應罐壁的熱流量增加3 325.9 W，儲罐罐壁總熱流量增加47.2 W，儲罐充滿率對儲罐罐壁熱流量的影響較小。

5 結論

根據(jù)LNG儲罐的結構，給出了儲罐熱流量和靜態(tài)蒸發(fā)率的計算方法，利用該方法對濟南某3×10

m

雙金屬LNG儲罐熱流量和靜態(tài)蒸發(fā)率進行計算，分析儲罐熱流量的影響因素，結果表明：

① 案例儲罐的靜態(tài)蒸發(fā)率為0.074 8%，靜態(tài)蒸發(fā)率較小，儲罐的保冷性能良好。

② 環(huán)境溫度對儲罐熱流量的影響較大，環(huán)境風速和充滿率對儲罐熱流量的影響較小。

[1] 彭明,丁乙. 全容式LNG儲罐絕熱性能及保冷系統(tǒng)研究[J]. 天然氣工業(yè),2012 (3):94-97.

[2] 王武昌,李玉星,孫法峰,等. 大型LNG儲罐內壓力及蒸發(fā)率的影響因素分析[J]. 天然氣工業(yè),2010 (7):87-92.

[3] 李曉明,許燕,姚淑婷,等. LNG儲罐日蒸發(fā)率測量方法及計算[J]. 低溫與超導,2013 (9):29-32.

[4] 楊世銘,陶文銓. 傳熱學[M]. 4版. 北京:高等教育出版社,2015.

[5] 李兆慈,郭保玲,嚴俊偉. LNG儲罐溫度場計算及影響因素分析[J]. 油氣儲運,2015 (3):244-247.

[6] 顧安忠，魯雪生，金國強，等. 液化天然氣技術手冊[M]. 北京:機械工業(yè)出版社，2010.

[7] 崔永章，史永征，陳彬劍. 燃氣氣源[M]. 北京: 機械工業(yè)出版社，2013: 285-286.

[8] 高華偉，張樹龍，李秀娟，等. 天然氣黏度和熱導率的計算[J]. 煤氣與熱力，2008 (5):B10-B13.

[ 9] 黃永華，陳國邦. 低溫流體熱物理性質[M]. 2版. 北京: 國防工業(yè)出版社，2014：428-432.

[10] 詹淑慧. 燃氣供應[M]. 2版. 北京: 中國建筑工業(yè)出版社，2011:31-32.

[11] 段常貴. 燃氣輸配[M]. 5版. 北京: 中國建筑工業(yè)出版社，2015:26-27.

[12] 劉璐，劉寶玉，李少華，等. LNG熱物性計算[J]. 當代化工，2010(6): 696-698.

[13] 電子工業(yè)部第十設計研究院. 空氣調節(jié)設計手冊[M]. 2版. 北京:中國建筑工業(yè)出版社,1995.

[14] 金明皇,許克軍,程松民,等. 大型LNG儲罐靜態(tài)日蒸發(fā)率的計算方法[J]. 油氣儲運,2016 (4): 386-390.