氨制冷機房氨氣泄漏擴散的數值模擬

董曉強 李 坤 劉玉婷 李俊明

（1清華大學熱能工程系 北京 100084；2國民核生化災害防護國家重點實驗室 北京 102205；3國內貿易工程設計研究院 北京 100069）

氨制冷機房氨氣泄漏擴散的數值模擬

董曉強1，2李 坤1，3劉玉婷1李俊明1

（1清華大學熱能工程系 北京 100084；2國民核生化災害防護國家重點實驗室 北京 102205；3國內貿易工程設計研究院 北京 100069）

本文選取氨制冷機房作為典型場景，應用計算流體力學方法研究氨氣泄漏擴散規律，采用Fluent軟件對無通風及上下通風情況下的制冷機房氨氣泄漏擴散進行了數值模擬。分析了氨氣泄漏后不同時刻制冷機房縱剖面和橫截面氨氣濃度分布，研究了不同泄漏量下氨氣濃度隨高度變化的特征，提出了氨氣報警器偏上安裝的建議。對比上進風下排風和下進風上排風兩種應急排風方案泄漏擴散的數值模擬結果，提出了應急排風進風口偏下設置的建議。

氨；泄漏；擴散；模擬；制冷機房

氨在世界上的使用已經有150多年歷史，在我國也是使用最早的一種制冷工質。氨具有制冷效率高、能效系數大、熱物性好、價格低廉、自然環保等優勢，作為大型制冷系統的制冷工質性價比非常高。因此，我國約80%的冷庫、冷凍加工項目均采用氨作為制冷工質［1－2］。為提高效率，目前還出現了氨與其他工質復合使用的趨勢［3］。在“臭氧層殺手”氟利昂逐步替代后，GWP和ODP均為零且天然存在的氨制冷劑重新受到各國的重視，又開始采用氨工質替代氟利昂。目前，在日本和美國，氨作為制冷工質的大型制冷系統占比大于80%。近年來，導致多起人員傷亡事故。因此在泄漏早期進行快速監測報警與妥善應急處置至關重要。深入研究氨泄漏傳質特性，科學確定氨氣報警器的安裝位置，實現氨泄漏后及時準確地監測氨氣濃度，合理制定應急排風及稀釋噴淋頭的布點方案，對于避免氨泄漏造成的生命和財產損失有重要意義。

歐美日等發達國家在泄漏事故研究方面一直位于世界前列，并取得了一系列的研究成果。由于國外氨泄漏造成亡人事故并不多見，所以研究氨泄漏事故的研究工作也較少。近年來，我國研究機構也相繼開展了危險化學品泄漏事故的研究。但2013年以前，國內室內氨泄漏引發的中毒死亡事故也相對較少，所以對氨泄漏擴散的研究相對較少，且主要集中于采用高斯模式研究氨氣泄漏擴散過程［4］。 王志鵬［5］對外場液氨儲罐區域的泄漏擴散進行了模擬，Ji Jie等［6］進行了外場液氨儲罐區域泄漏擴散的可視化模擬，Pang Bing等［7］通過數值模擬研究了液氨儲罐泄漏事故中不同泄漏孔大小、位置和泄漏速率的影響，熊立春等［8］對引入時間疊加的高斯模式在液氨泄漏擴散進行了模擬研究并指導人員疏散，杜娟麗等［9］研究了冷庫制冷管路的檢漏和定位模型，王如竹等［10－12］研究了氨吸收傳質特性和滲透率，陳小磚等［13－14］對垂直管、橫紋管氨水溶液降膜吸收特性進行了研究。上述研究很少涉及室內氨氣泄漏模擬，對氨氣報警器的安裝位置和應急排風進出風口位置選擇的研究更少。本文選取氨控制機房作為典型場景，采用計算流體力學方法研究典型的氨氣泄漏擴散規律，利用FLUENT軟件對室內氨氣在無通風情況進行了數值模擬，為確定最佳氨氣報警器安裝位置和應急排風進出風口位置提出了相應的建議。

1 模型與方法

1.1 模型基本假設

1）初始狀態時，假設儲罐內初始溫度與儲罐周邊環境溫度相同，泄漏發生后泄漏點初始溫度為液氨沸騰相變溫度240 K；2）氨氣泄漏速度不隨時間的改變而改變；3）泄漏過程中所涉及氣體均為理想氣體；4）泄漏口處液氨已轉變為氨氣，不涉及相變，在模擬區域擴散的也為氨氣。

1.2 控制方程

氨氣擴散控制方程包括連續性方程、能量守恒方程、動量守恒方程、組分運輸方程。

1）連續性方程

2）能量守恒方程

3）動量守恒方程

4）組分運輸方程

式中：Sm為質量源項，kg/（m3·s）；E為內能，J；k為導熱系數，W/（m·K）；hi為i組分的焓，J/kg；Ji為擴散通量，kg/（m2·s）；Φ為能量耗散項，J/（m3·s）；Sh為能量源項，J/（m3·s）；τ為切應力，N/m2；mi為質量分數；Si表示i組分的質量源項，kg/（m3·s）。

2 模型的構建

2.1 幾何模型

本文選擇北京某公司氨制冷機房作為模擬對象。兩個高壓儲罐單獨位于1個隔間；制冷機房中有2個中壓儲罐，4個低壓儲罐，9個壓縮機組；控制室位于制冷機房的門口左側單獨隔開，并通過一個小門與機房相通。整個區域尺寸為24 m×12 m×7.1 m，采用Solidworks建立三維幾何模型，如圖1所示。

圖1 氨制冷機房幾何模型Fig．1 Geometric model of ammonia refrigeration room

2.2 網格劃分與邊界條件

采用ICEM對模型進行網格劃分，泄漏孔所在位置設為INLET，壓縮機散熱面設為 HEAT，地面、屋頂、墻面、無泄漏孔的儲罐和無散熱面的壓縮機表面設為絕熱壁面，制冷機房中的空氣設為FLUID1，控制機房中的空氣設為FLUID2，網格最大單元為200，采用非結構化網格生成四面體網格，生成的網格如圖2所示。

邊界條件設置：1）熱源：壓縮機，壁面邊界類型，恒溫323 K；2）入口：儲罐泄漏孔，質量流量入口，0.1 kg/s或 2.4 kg，240 K；3）流體：房間內空氣；4）出口：房門，壓力出口；5）壁面：房間墻壁。

2.3 網格簡化

圖2 氨制冷機房網格劃分Fig．2 Mesh partition of ammonia refrigeration room

由于幾何模型復雜、網格數量巨大，難以進行擴散模擬，因此對上述幾何模型按“應急通風氨氣泄漏傳質模擬”及“制冷機房及控制室氨氣泄漏傳質模擬”兩種情況進行簡化處理。為了模擬應急通風對液氨泄漏擴散的影響，將系統簡化為一個泄漏孔、一個壓縮機和兩對進、出風口的的制冷機房，采用ICEM進行網格劃分，生成網格如圖3（a）所示。為了模擬制冷機房內液氨泄漏對控制室中操作人員的影響，刪除圖2中的高壓儲罐所在封閉區域，將制冷機房中的壓縮機和低壓儲罐均簡化為長方體，其中一個低壓儲罐的截面為泄漏面，所有壓縮機均有一個面為熱源面。為模擬控制室與制冷機房之間房門常處于關閉狀態，將房門四周的縫隙設置為一個環狀內部面，兩個空間中的流體可以自由進出，并設置一個上部窗戶為自然排風狀態，生成網格如圖3（b）所示。

圖3 氨氣泄漏傳質簡化網格Fig．3 Simplified mesh of ammonia leakage mass transfer

3 無通風情況下氨氣泄漏擴散模擬

按圖3所示的簡化網格，模擬了無通風情況下400 s內0.1 kg/s泄漏速度下的氨氣濃度分布。由于氨氣密度僅為空氣的0.6倍，而距地面1.5 m附近是人體受氨氣毒性傷害的視覺、呼吸等主要器官所處高度，故圖4主要給出了氨氣泄漏后機房內距地面高1.5 m處，氨氣濃度隨泄漏時間變化分布的數值模擬。由圖4（a）可以看出，泄漏3 s后人體呼吸區域內的氨氣濃度達到可察覺臭味的12 mg/kg；由圖4（b）可以看出，泄漏4 s后人眼高度區域內的氨氣濃度達到對眼睛造成輕度刺激的24 mg/kg；由圖4（c）可以看出，泄漏6 s后人員呼吸高度的氨氣濃度已經全部達到長期工作危害劑量60 mg/kg；由圖4（d）可以看出，泄漏9 s后人員工作區域內的氨氣濃度已經達到接觸后會造成嚴重的喉部、鼻孔和上呼吸道刺激的240 mg/kg；由圖4（e）可以看出，泄漏 11 s后人員呼吸高度的氨氣濃度達到如果接觸時間超過30 min會造成永久性傷害的420 mg/kg；由圖4（f）可以看出，泄漏16 s后人員呼吸高度的氨氣濃度達到造成嚴重咳嗽、影響中央神經系統、導致痙攣、可能致命的1 020 mg/kg；由圖4（g）可以看出，泄漏 35 s后人員呼吸高度的氨氣濃度達到會造成嚴重浮腫（體內液體異常累積，導致膨脹、窒息，喪失意識），立即致命的 3 000 mg/kg。

依據《氨氣檢測報警儀技術規范》［15］規定，氨氣檢測報警儀的響應時間為泵吸式不大于60 s，擴散式不大于90 s，通過上述泄漏過程的時間分析可知，氨氣泄漏發生后擴散極快。實際上大多數氨氣報警器平時為達到長期工作的可靠性大多數時間均處于擴散檢測模式，這種情況下報警限附近的報警響應時間非常慢。為確定氨氣報警器最佳安裝位置，本文研究了氨氣以0.1 kg/s速度泄漏90 s后1.5 m 高和離屋頂0.1 m（即 7.0 m 高）處的氨氣濃度分布，結果表明：屋頂（7.0 m）氨氣濃度遠高于人身高度（1.5 m）的氨氣濃度。由此可見，90 s響應時間過慢，應配備一定數量的傳感器外置非防爆型快速響應氨氣報警器。氨氣泄漏后最重要的就是事故現場隔離和周邊人員疏散，不僅要有最短的疏散路線，還要經常進行事故預案宣貫和演練，避免出現處置失誤而造成人身和財產損失。如果泄漏濃度較低，可以啟動應急排風。如果泄漏濃度達到420 mg/kg，會造成人身永久傷害，必須進行隔離疏散，應急處置只能采用自動噴淋等無人作業方式。當泄漏濃度低于420 mg/kg，且必須人員進入處置時，只有配帶防護設備專業人員才能進行處置作業，否則很容易造成人員傷亡。

圖4 氨氣泄漏后的1.5 m高處的濃度分布Fig．4 Ammonia concentration distribution of 1.5 m height after leakage

4 制冷機房及控制室氨氣泄漏傳質模擬

根據圖3所示的簡化網格，當泄漏量為0.1 kg/s，研究了復雜熱環境下氨氣泄漏時控制室內氨氣濃度分布，如圖5所示。

由圖5（a）～（h）可以看出：氨氣泄漏后氨制冷機房下部由于存在液氨罐和壓縮機等障礙物，容易出現氨氣濃度很低的擴散死角，因此制冷機房下部安裝氨氣報警器很容易造成漏報或遲報；由圖5（i）～（p）可以看出：氨氣泄漏后除泄漏點附近外，氨制冷機房頂部氨氣濃度高于房間下部濃度，控制室內濃度低于制冷機房內氨氣濃度，氨氣濃度隨高度的增加而不斷升高，關閉控制室的門窗能夠大幅降低氨氣向控制室內擴散。

圖5 氨氣泄漏后制冷機房的濃度分布Fig．5 Ammonia concentration distribution of refrigeration room after leakage

5 結論

通過對氨制冷機房氨氣泄漏擴散過程進行數值模擬，得出如下結論：

1）氨氣泄漏發生后擴散速度極快，以0.1 kg/s的小泄漏速率，11 s內氨氣濃度即可達到致使現場人員永久傷害的700 mg/kg，該時間遠低于國家標準規定的氨氣檢測報警儀的響應時間（泵吸式≤60 s以及擴散式≤90 s）。

2）氨氣泄漏后，氨制冷機房內房間頂部的氨氣濃度高于底部的氨氣濃度，氨制冷機房內下部容易出現氨氣濃度很低的擴散死角。

3）氨制冷機房與控制室之間門窗關閉，可延緩氨氣向控制室擴散，大幅降低控制室內氨氣濃度。

根據以上結論，提出如下建議：

1）氨氣泄漏發生后，最重要的就是現場隔離與人員疏散，制冷機房中要設置疏散通道，還應編制事故預案，重在專業培訓和處置演練，防止錯誤處置。

2）分析不同泄漏時間氨氣濃度的分布，在泄漏點位置不能事先確定的情況下，沿房頂均勻布點并兼顧可能泄漏點是最佳氨氣報警器安裝方案。為兼顧報警器長期工作的可靠性和火災爆炸事故后工作的可靠性，并快速實現氨氣泄漏報警，應同時使用傳感器外置的快速響應的非防爆型氨氣報警器（響應時間更快）和響應時間較慢的防爆型氨氣報警器（事故后更可靠）。

3）在應急排風裝置滿足換氣次數和排風量不小于國家相關規范要求值的前提下，采用下進風上排風的應急排風是一種優選方案，可以保證操作人員安全，且利于有毒氣體及時排出。

4）操作維護人員所在控制室必須與氨制冷機房和氨制冷系統完全隔離，制冷機房與控制室之間大門關閉，控制室中操作人員便有足夠疏散和應急處置準備時間。

［1］楊一凡．氨制冷技術的應用現狀及發展趨勢［J］．制冷學報，2007，28（4）：12?19．（YANG Yifan．Application and development of ammonia refrigeration technology［J］．Journal of Refrigeration，2007，28（4）：12?19．）

［2］張建一，徐穎．國內外大中型冷庫制冷劑的現狀和發展動向［J］．制冷學報，2009，30（4）：51?57．（ZHANG Jianyi，XU Ying．Development on refrigerants used in large and medium?size refrigerated warehouses［J］．Journal of Refrigeration，2009，30（4）：51?57．）

［3］卞宜峰，何國庚，蔡德華，等．吸收式制冷工質對的研究進展［J］．制冷學報，2015，36（6）：17?26．（BIAN Yifeng，HE Guogeng，CAI Dehua，et al．Research pro?gress of absorption refrigeration working pairs［J］．Journal of Refrigeration，2015，36（6）：17?26．）

［4］黃斌．高斯模式在液氨泄漏事故環境污染控制中的應用研究［D］．昆明：昆明理工大學，2007．（HUANG Bin．Research on controlling environmental pollution in liquid ammonia leakage accident by gaussian model［D］．Kun?ming： Kunming University of Science and Technology，2007．）

［5］王志鵬．氨氣泄漏數值模擬及應急響應措施研究［D］．北京：中國地質大學，2013．（WANG Zhipeng．Response measures and numerical simulation for ammonia leakage［D］．Beijing： China University of Geosciences，2013．）

［6］JI Jie，CHEN Xiaoxia，HAN Xuefeng．Rapid simulation and visualization analysis of liquid ammonia tank leakage risk［J］．Procedia Engineering，2014，84： 682?688．

［7］PANG Bing，LIANG Chen，MA Lan，et al．Numerical simulation of leakage and diffusion of liquid ammonia tank［C］//2014 Fourth International Conference on Instrumen?tation and Measurement，Computer，Communication and Control（IMCCC）．Harbin： IEEE，2014： 216?220．

［8］熊立春，陳建宏，石東平．引入時間疊加的高斯液氨泄漏擴散模擬及人員疏散研究［J］．中國安全生產科學技術，2015，11（11）：76?82．（XIONG Lichun，CHEN Jian?hong，SHI Dongping．Simulation on liquid ammonia leakage by Gaussian model introducing temporal superposition and personnel evacuation［J］．Journal of Safety Science and Technology，2015，11（11）：76?82．）

［9］杜娟麗，田紳，邵雙全，等．冷庫制冷劑管路檢漏與定位實時模型研究［J］．制冷學報，2015，36（5）：43?48．（DU Juanli，TIAN Shen，SHAO Shuangquan，et al．Real?time model for leak detection and positioning of liquid pipe?lines in the refrigeration system［J］．Journal of Refrigera?tion，2015，36（5）：43?48．）

［10］董曼，夏再忠，王如竹，等．大壓差下靜止氨水表面吸收特性研究［J］．制冷學報，2014，35（2）：1?5．（DONG Man，XIA Zaizhong，WANG Ruzhu，et al．Characteristic research of ammonia vapor absorption into the surface of stagnant ammonia?water under big pressure difference［J］．Journal of Refrigeration，2014，35（2）：1?5．）

［11］江龍，金哲權，王麗偉，等．不同吸附量條件下氯化鍶混合吸附劑—氨的導熱系數與滲透率研究［J］．制冷學報，2013，34（3）：35?39．（JIANG Long，JIN Zhequan，WANG Liwei，et al．Study on the thermal conductivity and permeability of composite SrCl2?NH3working pair under the condition of different values for adsorption quantity［J］．Journal of Refrigeration，2013，34（3）：35?39．）

［12］徐振中，夏再忠，林芃，等．垂直管氨水降膜吸收傳質研究［J］．制冷學報，2012，33（1）：1?5．（XU Zhenz?hong，XIA Zaizhong，LIN Peng，et al．Mass transfer study of ammonia?water falling?film absorption in vertical tube［J］．Journal of Refrigeration，2012，33（1）：1?5．）

［13］陳小磚，柳建華，張良，等．基于橫紋管氨水溶液垂直管外降膜吸收性能的變化規律［J］．制冷學報，2013，34（3）：50?55．（CHEN Xiaozhuan，LIU Jianhua，ZHANG Liang，et al．Performance variation of aqua ammonia falling film absorption outside vertical tube based on transversally grooved tube［J］．Journal of Refrigeration，2013，34（3）：50?55．）

［14］陳小磚，柳建華，張良，等．橫紋管氨水降膜吸收溫度場與濃度場計算與測量［J］．制冷學報，2013，34（6）：29?34．（CHEN Xiaozhuan，LIU Jianhua，ZHANG Liang，et al．Calculation and measurement of temperature field and mass fraction field about aqua ammonia falling film absorp?tion based on transversally grooved tube［J］．Journal of Re?frigeration，2013，34（6）：29?34．）

［15］氨氣檢測報警儀技術規范：AQ/T 3044—2013［S］．北京：中國標準出版社，2013．（Technical specifications for ammonia detection instruments： AQ/T 3044—2013［S］．Beijing： Standards Press of China，2013．）

Numerical Simulation of the Ammonia?leakage Diffusion in Ammonia Refrigeration Room

Dong Xiaoqiang1，2Li Kun1，3Liu Yuting1Li Junming1
（1．Department of Thermal Engineering，Tsinghua University，Beijing，100084，China；2．State Key Laboratory of National Nuclear and Biological Disaster Protection，Beijing，102205，China；3．Internal Trade Engineering Design and Research Institute，Beijing，100069，China）

In this study，an ammonia refrigeration room was selected as a typical scenario．Computational fluid dynamics（CFD） was applied to study the rules of the ammonia?leakage diffusion．The ammonia?leakage diffusion in a refrigeration room with up?down ventila?tion and without ventilation was numerically simulated using the Fluent CFD software package．The ammonia concentrations in the longi?tudinal profile and the cross?section of the refrigeration room were analyzed at different times after the ammonia leakage．The numerical simulation results show that the ammonia concentration increases with height in the refrigeration room，and a lower fresh?air inlet can quickly decrease the ammonia concentration．Based on the above results，it is proposed that an ammonia?gas alarm should be installed near the ceiling to detect the leakage，and an emergency exhaust inlet should be installed near the floor to decrease the ammonia con?centration．

ammonia；leakage；diffusion；simulation；refrigeration room

Li Junming，male，Ph．D．/professor，Department of Thermal En?gineering，Tsinghua University，＋86 10?62771001，E?mail： lijm＠ tsinghua．edu．cn．Research fields： micro flow and heat trans?fer，condensation heat transfer，energy conservation of air condi?tioning system，heat and mass transfer of hazardous chemicals leakage．

TB651；TB64；TP391．9

A

0253－4339（2017）06－0012－08

10.3969 /j.issn.0253 － 4339.2017.06.012

科技部中英政府間合作項目（2015DFG62060）和國家973項目（2011CB706900）資助。 （The project was supported by the Innovate UK?China Most Project（No．2015DFG62060） and National 973 Project（No．2011CB706900）．）

2017年1月3日

李俊明，男，教授，博士生導師，清華大學熱能工程系工程熱物理研究所，（010）62771001，E?mail： lijm＠ tsinghua．edu．cn。 研究方向：微細流動與傳熱，凝結換熱，空調制冷系統節能，危險化學品泄露的傳熱傳質過程。