冷庫氨泄漏擴(kuò)散規(guī)律探析

甄，王磊，郭靖，王長江

（哈爾濱商業(yè)大學(xué)能源與建筑工程學(xué)院，哈爾濱150028）

冷庫氨泄漏擴(kuò)散規(guī)律探析

（哈爾濱商業(yè)大學(xué)能源與建筑工程學(xué)院，哈爾濱150028）

本文介紹了冷庫氨泄漏的基本模式，以及氨氣泄漏與氨液泄漏的流量計(jì)算方法。以高貯器為例，利用Fluent分別模擬氨氣與氨液泄漏后濃度場分布，并分析氨擴(kuò)散規(guī)律，為氨泄漏事故預(yù)防及應(yīng)急處置提供參考。

冷庫；氨氣；氨液；泄漏；Fluent

氨作為一種天然制冷劑具有環(huán)境友好、價格低廉、制取方便的特點(diǎn)，占據(jù)著我國80%以上的冷庫制冷劑使用份額。同時氨也具有毒性、爆炸性等危險特性，一座萬噸級的氨冷庫其氨的沖注量可達(dá)20t以上，從而使冷庫成為較大的危險源。

據(jù)統(tǒng)計(jì)，目前我國庫齡在30年以上的冷庫占全國冷庫總量的1/2，這些高齡冷庫由于年久失修，設(shè)備老化從而造成泄漏事故頻發(fā)，近些年的幾起氨泄漏事故更是造成大量人員傷亡，使人們談氨色變。表1列舉出僅2015上半年我國發(fā)生的氨泄漏安全事故。

表1 2015年部分氨泄漏事故統(tǒng)計(jì)表Tab.1 part of ammonia leakage accidents statistics in 2015

其實(shí)事故的罪魁禍?zhǔn)撞⒉皇前北旧恚驹蛟谟谌藛T操作不規(guī)范，設(shè)備疏于維護(hù)，以及泄漏時應(yīng)急處置不當(dāng)。以氨作為制冷劑的冷庫，其制冷工藝流程非常復(fù)雜，整個系統(tǒng)中高低壓設(shè)備縱橫交錯，閥門、接頭數(shù)不勝數(shù)，氨在系統(tǒng)中的狀態(tài)更是變幻多端，任何一個部位都可能成為氨泄漏事故的源頭。因此，不同狀態(tài)的氨在不同部位泄漏時具有不同的泄漏規(guī)律，對不同泄漏規(guī)律的研究對于事故的預(yù)防及應(yīng)急處置具有重要意義。

1 氨的泄漏模式

圖1 常見泄漏情況Fig.1 Common ammonia leakage

氨泄漏大致可分為大面積泄漏、小孔泄漏及管道泄漏。大面積泄漏是指在短時間內(nèi)有大量的氨泄漏出來，貯罐的超壓爆炸就屬于大面積泄漏，這種情況在冷庫事故中并不常見。小孔泄漏是指氨以非常緩慢的速度持續(xù)地從小孔泄漏到外界。管道泄漏是指管道截面斷裂或泄漏孔徑與管道直徑相當(dāng)?shù)臅r候發(fā)生的泄漏。圖1顯示了常見的泄漏情況。

1.1 氨氣泄漏流量計(jì)算

1.1.1 氨氣經(jīng)小孔泄漏

氣體含有的能量從小孔泄漏或擴(kuò)散出去時，在壓力的作用下轉(zhuǎn)化為動能。氣體經(jīng)孔流出的過程中，其密度、壓力和溫度都發(fā)生著變化。若泄漏為自由擴(kuò)散，絕大部分壓力能轉(zhuǎn)化為動能，可假設(shè)為等熵過程。此時僅需獲取孔洞直徑參數(shù)。

圖2 氨氣自由擴(kuò)散泄漏Fig.2 Free diffusion leakage of ammonia gas

質(zhì)量流率的表達(dá)式：

式（1）描述了在等熵膨脹中任意點(diǎn)處的質(zhì)量流率。

表2 孔口形狀與泄漏系數(shù)和當(dāng)量直徑的關(guān)系Tab.2 The relationship between the shape of orifices,leakage factor and equivalent diameter

對于許多安全性研究，都需要計(jì)算通過小孔泄漏的最大流量，引起最大流速的壓力比為：

塞壓Pchoked為流量最大時下游最大壓力。當(dāng)下游壓力小于Pchoked時：（1）對于大多數(shù)情況，在小孔處氣體流速是聲速；（2）降低下游壓力，亦不能進(jìn)一步增加其流速及質(zhì)量流量，他們獨(dú)立于下游環(huán)境。通常稱這種流動為塞流、聲速流或臨界流，如3圖所示。

圖3 氨氣通過小孔的塞流Fig.3 Plug flow of ammonia gas pass through the small hole

確定最大流量為：

式中，γ是熱融比；M是泄漏氣體的摩爾質(zhì)量；T0是泄漏源的溫度（K）；Rg是理想氣體常數(shù)。

在雷諾數(shù)大于30 000的鋒利的小孔情況下，流出系數(shù)C0取0.61。然而，對于塞流和C0不確定的情況，推薦使用1.0。

1.1.2 氨氣經(jīng)管道泄漏

真實(shí)氣體流動介于絕熱和等溫之間，根據(jù)Crane實(shí)驗(yàn)，可按絕熱法建立流動模型。

有氣體流動的絕熱管道如圖4所示。對這一情況，出口處流速低于聲速，沿管道的壓力梯度驅(qū)動流動。當(dāng)氣體流經(jīng)管道時，因壓力下降而膨脹，導(dǎo)致氣體速度及動能增加，溫度升高。然而，氣體與管壁之間的摩擦力會使氣體溫度降低。氣體最終溫度的升降，取決于動能和摩擦力做功的大小。

圖4 氨氣通過管道的絕熱非臨界流動Fig.4 Adiabatic noncritical flow of ammonia gas pass through the pipe

式（4）即為氨氣經(jīng)管道泄漏的質(zhì)量流量表達(dá)式。

1.2 氨液泄漏流量計(jì)算

存在于制冷系統(tǒng)中的液氨往往高于其飽和蒸氣壓，若儲罐、管道或其他設(shè)備出現(xiàn)孔洞，部分液體會閃蒸為蒸氣，有時會發(fā)生爆炸。

閃蒸發(fā)生的速度非常快，過程可設(shè)為絕熱。過熱液體中的額外能量使液體氣化，并使溫度降到新的沸點(diǎn)。如果m是初始液體的質(zhì)量，cp是液體的比熱容，T0是降壓后液體的沸點(diǎn)，則包含在過熱液體中額外的能量為：

該能量使液體蒸發(fā)。若△Hv是液體的蒸發(fā)熱，則蒸發(fā)的液體質(zhì)量mv為：

液體蒸發(fā)比例是：

式（7）是基于液體的物理特性在T0到Tb的溫度范圍內(nèi)不變得到的。更一般的表達(dá)形式如下所述。

溫度T的變化導(dǎo)致液體質(zhì)量m的變化為：

在初始溫度T0（液體質(zhì)量為m）與最終沸點(diǎn)溫度Tb（液體質(zhì)量為m-mv）區(qū)間內(nèi)，對式（8）進(jìn)行積分，得到：

因?yàn)閮上嗔鞯拇嬖?，通過孔洞和管道泄漏出的閃蒸液體需要被考慮進(jìn)來。如果泄漏的沿程很短，則存在不平衡條件，液體在孔洞外閃蒸，流出方程應(yīng)選用不可壓流體。

若泄漏的沿程大于10cm，即可達(dá)到平衡閃蒸條件，且流動為塞流。假設(shè)塞壓與閃蒸液體的飽和蒸氣壓相等，質(zhì)量流率計(jì)算式為：

式中，A是釋放面積（m2）；C0是流出系數(shù)；ρf是液體密度（kg/m3）；p是儲罐內(nèi)壓力（Pa）；ρsat是閃蒸液體處于周圍溫度情況下的飽和蒸氣壓（Pa）。

2 基于Fluent的模擬舉例

Fluent采用計(jì)算流體力學(xué)方法模擬擴(kuò)散湍流的流動過程。通過建立質(zhì)量、動量、能量及組分等基本守恒方程，并設(shè)置適當(dāng)初始條件和邊界條件，運(yùn)用數(shù)值計(jì)算理論和方法，實(shí)現(xiàn)各種場的分布和預(yù)測，借以描述擴(kuò)散過程。其描述湍流動能的運(yùn)輸和湍流黏性系數(shù)的k-ε模型被廣泛應(yīng)用。

本文以某氨冷庫高貯器發(fā)生管道泄漏為例，模擬氨氣與氨液發(fā)生泄漏時的不同擴(kuò)散規(guī)律。

2.1 建立模型及模擬條件設(shè)置

假設(shè)液氨存儲在溫度25℃、壓力2.16Mpa、容積為

2.43 m3（液氨占80%）的高貯器中，發(fā)生泄漏部位是與高貯器相連接的管道，管道直徑100mm。由于資源有限，選用二維網(wǎng)格進(jìn)行模擬，泄漏空間為長度10m高5m制冷機(jī)房，機(jī)房一側(cè)有高為2m的門作為模擬出口，進(jìn)口位于模擬空間底邊的中央，泄漏方向垂直于底邊，泄漏過程中進(jìn)口上游條件不發(fā)生變化。

設(shè)定邊界參數(shù)時，選取湍流強(qiáng)度和水力半徑項(xiàng)作為TurbulenceSpecification Method的方法，可避免用二維網(wǎng)格模擬出來的結(jié)果非狹縫泄漏。采用標(biāo)準(zhǔn)κ-ε模型，開啟組分輸運(yùn)模型和氣化冷凝模型。

2.2 模擬結(jié)果及分析

根據(jù)上述所建模型和條件設(shè)置，分別模擬氨氣和氨液通過管道泄漏，模擬結(jié)果如下所示：

圖5 氨氣泄漏1s時濃度分布Fig.5 Concentration distribution at the 1s of ammonia gas leakage

圖6 氨氣泄漏5s時濃度分布Fig.6 Concentration distribution at the 5s of ammonia leakage

圖7 氨液泄漏1s時濃度分布Fig.7 Concentration distribution at the 1s of ammonia liquor leakage

圖8 氨液泄漏5s時濃度分布Fig.8 Concentration distribution at the 5s of ammonia liquor leakage

從模擬的濃度分布圖中可以看出，對于高貯器來說，無論是氨氣泄漏還是氨液泄漏，氨的擴(kuò)散速度都十分迅速，發(fā)生泄漏5s后整個模擬機(jī)房空間已全部被氨氣占據(jù)，可見氨泄漏事故的可怕程度，若泄漏時機(jī)房有作業(yè)人員留給其逃生時間極其短暫。但從泄漏開始到充滿整個空間，氨擴(kuò)散的方式并不同。氨氣泄漏時，氨氣濃度逐漸從低處，遠(yuǎn)離出口方向，向高處接近出口方向增加。而氨液泄漏則不然，氨氣濃度從高處，接近出口方向，向低處遠(yuǎn)離出口方向增加。這是由于泄漏發(fā)生在相對受限的空間內(nèi)且氨氣的密度小于空氣，泄漏發(fā)生時，受限空間內(nèi)的總壓升高，氣體會向出口流動，并在空間內(nèi)形成渦旋，圖5可以清晰地看出氣體的流動狀態(tài)，泄漏的氨氣會隨著渦旋擴(kuò)展開來。而氨液泄漏時，高速的氨液噴射到頂部并開始迅速氣化，急速向出口下方方向流動擴(kuò)散，從圖7可以看出剛泄漏出的氨液還沒有完全氣化，這使得泄漏處反而低于頂部。

3 結(jié)語

俗話說，知己知彼，百戰(zhàn)不殆。目前，人們對冷庫不同設(shè)備發(fā)生氨泄漏擴(kuò)散的規(guī)律仍知之甚少，本文僅以高貯器發(fā)生氨氣及氨液泄漏為例初步探索了氨泄漏規(guī)律，至于其他設(shè)備及不同泄漏條件仍有待進(jìn)一步研究。只要我們能夠知“氨”善用，一定能夠防范事故的發(fā)生，若事故真的發(fā)生，也一定能夠正確、及時地做出處置措施。

［1］楊一凡.氨制冷技術(shù)的應(yīng)用現(xiàn)狀及發(fā)展趨勢［J］.制冷學(xué)報(bào)，2007，（04）：12-19.

［2］DanietA.Crowl等.化工過程安全理論及應(yīng)用［M］.蔣軍成，等，譯.北京：化學(xué)工業(yè)出版社，2006.

［3］張偉.氨制冷系統(tǒng)壓力管道檢驗(yàn)及其安全評價分析［J］.石化技術(shù)，2015，（04）：44.

［4］王治華.受限空間內(nèi)氣體擴(kuò)散的數(shù)值模擬及分析［D］.大連：大連理工大學(xué)，2009.

［5］余運(yùn)波.沸騰液體膨脹蒸汽爆炸模型修正探討［J］.湖南安全與防災(zāi)，2011，（03）：46-47.

Study on the Diffusion Law of Ammonia Leakage in Cold Storage

ZHENBing,WANGLei,GUOJing,WANGChang-jiang
(School ofEnergyand Architectural Engineering,Harbin UniversityofCommerce,Harbin,150028,China)

In this paper,the basic mode of ammonia leakage in cold storage and the method for calculating the flowrate of ammonia leakage and ammonia liquid leakage was demonstrated.Taking the high pressure reservoir as an example,the concentration field distribution of ammonia vapor and ammonia liquid leakage was simulated by Fluent,and the ammonia diffusion law was analyzed.This work provides reference for the prevention and emergencydisposal ofammonia leakage accident。

Cold storage；Ammonia vapor；Ammonia liquid；Leakage；Fluent

x928.5

A

1674-8646（2015）10-0004-04

2015-08-19

黑龍江省教育廳科研項(xiàng)目（12531163）；哈爾濱商業(yè)大學(xué)博士科研啟動金（12DL003）

甄（1973-），男，山東莒南人，副教授，博士，在站博士后，從事制冷系統(tǒng)優(yōu)化設(shè)計(jì)、食品冷凍冷藏研究，e-mail:zb730812@sina.com。