油料倉庫拱頂儲油罐進出油氣相空間壓力模擬計算方法

孟凡芹，奚麗波，于百儉

（空軍勤務學院 航空油料物資系，江蘇 徐州 221000）

油料倉庫拱頂儲油罐進出油氣相空間壓力模擬計算方法

孟凡芹，奚麗波，于百儉

（空軍勤務學院 航空油料物資系，江蘇 徐州 221000）

各類大中型油庫中有為數眾多的露天立式拱頂油罐。該類油罐建造簡單、成本低，維護方便，但在進出油時容易引起油罐的漲罐和吸癟事故。給出拱頂儲油罐進出油簡化計算方法，模擬計算呼吸閥正常工作條件下油罐氣相空間壓力變化情況，為油罐安全運行和油庫科學管理提供參考。

油庫；拱頂儲油罐；氣相空間壓力；模擬計算

1 引言

金屬油罐分為浮頂罐和拱頂罐，其中拱頂金屬油罐由于結構簡單、造價低、維護方便，是各類油庫主要儲存設備，在大宗散裝液體物資流通中起到儲存、中轉作用。拱頂油罐進出油時靠呼吸閥控制罐頂空區（此部分為氣相空間）壓力，進出油流量超限或呼吸閥選擇不當，將發生油罐脹罐或吸癟事故。因此在油罐設計時，要確定呼吸閥規格和數量；收發油作業時，要嚴格控制進出油流量。油罐進出油流量一定時，油罐上部氣相空間壓力取決于呼吸閥、阻火器參數、油罐幾何參數、初始裝油高度等因素，是進出油作業時間的函數。本文利用熱力學方法，模擬計算拱頂油罐大呼吸過程氣相壓力變化，為大中型儲油罐呼吸閥選用及運行維護提供參考。

2 油罐工作條件

立式拱頂油罐如圖1所示。油罐氣相空間壓力、體積和溫度分別為Pg，Vg，Tg。罐頂部為呼吸閥，呼吸閥的參數有通徑d（m），開啟壓力為Ps1（正壓）、Ps2（負壓）（Pa），泄漏量vl等，模擬時取SY/T0511.1-2010標準規定值。呼吸時氣體出罐流量為Qto、吸氣時入罐氣體流量為Qti，單位為m3/s。油罐儲油高度為Ho、油罐立壁高度為H，單位為m。油罐處于進出油狀態，稱為油罐大呼吸狀態。此時設定油罐進油流量為Vin，出罐流量為Vout,單位為m3/s。油罐初始氣相空間為Vs，單位為m3。

圖1 露天立式拱頂金屬油罐

如圖1所示，油罐中油的體積Vo（單位：m3）為：

油罐中氣體的體積為圓柱部分的體積加罐頂部分的體積。其中拱頂管罐頂一般為球形，球形半徑為0.8～1.2D。計算時假設罐頂球形直徑為1.0D[1]，則油罐氣體空間體積Vg（單位：m3）為：

式中D為油罐直徑，這里假設油罐上不和底部直徑相等，為油罐平均直徑，單位為m。

模擬時設定油罐參數包括油罐直徑、高度、總體積、罐頂體積，呼吸閥型號參數。呼吸閥參數和數量決定了呼吸閥泄漏量、通氣量、開啟壓力等。根據以上設定條件，模擬油罐進出油時氣相空間壓力變化。模擬軟件還可以用圖形方式直觀顯示溫度、壓力、通氣量和呼吸閥開啟等動態變化情況。

油罐進出油時將油罐看作是帶有通氣閥的密閉容器。在呼吸閥作用下，油罐進出油過程可以看作是一個壓縮，膨脹和進氣、排氣過程。進出油流量一定，進氣排氣量隨著呼吸閥開啟情況而變化。由于時間短，進出油過程不考慮外部溫度變化的影響。

3 進油過程模擬

3.1 進油時氣相空間壓力計算

在呼吸閥作用下，油罐進油過程可以看作是一個壓縮和排氣過程。進油時液面上升，氣相空間被壓縮，氣相空間壓力上升。當氣相空間壓力小于呼吸閥開啟壓力時，呼吸閥以泄露流量向罐外排氣。隨著進油繼續，氣相空間壓力繼續增大，當大于呼吸閥開啟壓力時，呼吸閥開啟，排氣量增大。氣相空間壓力等于呼吸閥開啟壓力加呼吸閥阻力[1]。

為簡化計算，分兩步考慮。先不考慮呼吸閥通氣，將氣相空間看成一封閉空間，進油時簡化為氣體等溫壓縮過程，此時氣體壓力與時間τ（s）之間的關系可由下式表示：

式中pg0為分布計算時假設的壓力，單位pa。τ為計算周期，單位s。

再考慮閥門通氣量Vv（m3/s），不考慮進油，此時為定溫放氣過程[2]。設空氣密度為ρg=1.293kg/m3，空氣的氣體常數Rg=287.05J/(kg.K)，空氣的摩爾質量為 0.028 965kg/mol。該計算周期壓力變化Δpg（單位：pa）為：

模擬時分兩步計算，過程如下：

首先考慮等溫壓縮過程，得出該計算周期結束時的壓力pg1(n)。

然后再考慮定溫放氣過程壓力變化，得到第n步油罐氣相空間壓力pg(n)為：

3.2 進油時呼吸閥通氣量計算

呼吸閥通氣量Vv由閥門參數確定。根據SY/T0511.1-2010標準，呼吸閥的局部阻力系數只與流速有關，且在大于8m/s時變得恒定；小于6m/s時，隨著速度的減小有增大的趨勢。對于結構一定的閥門，其局部阻力損失系數與通徑無關。根據一般估算方法，正常通氣時，考慮阻火器和呼吸閥的局部阻力系數可選擇為ξ=6.5[3-4]。

當pg1(n)＞Ps1+1atm時，根據流體力學公式，Vv按下式計算：

式中d為呼吸閥通徑，單位為m。Vl為泄漏量，單位為m3/s。vl0為標準規定的泄漏量（m3/s），即0.75倍開啟壓力時的通氣量，見SY/T0511.1-2010表5。這里加泄漏量是為了計算結果的連續性。

在Ps1和Ps2范圍內，此時通氣量較小且不穩定。設定該階段呼吸閥流量特性曲線如圖3所示，當pg(n)=Ps1+1atm時或pg(n)=Ps2+1atm時，Vv=4vl0。當pg(n)=0.75Ps1+1atm時或pg(n)=0.75Ps2+1atm，Vv=vl0，當pg(n)=1atm時，Vv=0。當罐內氣相空間相對大氣壓力在0，0.75Ps1之間或0，0.75Ps2之間時，按線性插值計算，在0.75Ps1和 Ps1或0.75Ps2和Ps2之間時按拋物線插值計算。

圖2 呼吸閥泄漏量排氣量曲線示意圖

4 出油過程模擬

4.1 出油時氣相空間壓力計算

油罐出油過程與進油過程類似。在呼吸閥作用下，油罐出油過程可以看作是一個膨脹、充氣過程。出油時氣相空間壓力計算與進油時類似，同樣僅考慮出油流量對氣體空間的影響，并分兩步計算。先計算出油時的氣體等溫膨脹過程，再計算閥門吸氣時的等溫充氣過程。出油時計算方法如下：

首先計算單個步長等溫膨脹過程的壓力變化，Pg1(n)計算公式為：

4.2 進油時呼吸閥通氣量計算

出油時呼吸閥通氣量Vv計算與進油時類似。當pg1(n)≤1atm+Ps2時，根據流體力學公式，Vv按下式計算：

對比進氣時的計算方法，當1atm≥pg1(n)≥Ps2+1atm時，可認為沒有穩定的通氣量，此時通氣量較小。當pg(n)=Ps2+1atm時，Vv=4vl0。當pg(n)=0.75Ps1+1atm時，Vv=vl0。當pg(n)=1atm時，Vv=0。當罐內氣相空間相對大氣壓力在0，0.75Ps2之間時，按線性插值計算，在0.75Ps1和Ps1時按拋物線插值計算，如圖2所示。

5 模擬計算實例

假定某2 000m3油罐，安裝有150口徑呼吸閥和阻火器。對其進行進出油模擬計算的結果如圖3和4所示。圖中橫坐標為時間，單位為s，縱坐標為壓力，單位為Pa。圖3為進油流量分別為72m3/h和144m3/h時，罐頂氣相空間壓力隨時間變化曲線。圖4為出油流量分別為72m3/h和144m3/h時，罐頂氣相空間壓力變化曲線。從圖3和圖4可以看出，當流量為144m3/h時，無論出油還是進油，氣相空間壓力將超出一般油罐的設計壓力。進出油流量為72m3/h時，依靠呼吸閥可以滿足進出油排氣和吸氣的要求，氣相空間壓力不會超過油罐頂部設計壓力，這一結果與油罐設計數據一致。上述結果說明，按照標準選擇呼吸閥口徑和數量時，進出油作業通過呼吸閥排氣進氣，油罐壓力不會超過安全范圍。

圖3 油罐進油時氣相空間壓力變化

圖4 油罐出油時氣相空間壓力變化

6 結論

金屬拱頂油罐頂部氣相壓力允許范圍隨著使用年限的增加而逐步減小，收發油作業時油罐進出油流量也應隨著減小。本文所提供的模擬計算方法可以對不同規格油罐和呼吸閥進行模擬計算，得出不同壓力范圍下允許的最大進出油流量。按照最大收發油流量作業，既可保證油罐安全，又能最大程度提高收發油效率。從以上計算實例看，該方法計算簡單，對呼吸閥特性曲線簡化計算合理，計算結果可以作為收發油流量確定的依據。對不同體積、不同規格油罐，其模擬計算步長還可以進一步論證優化，以減小計算誤差提高運算速度。

[1]郭光臣,等.油庫設計與管理[M].東營:石油大學出版社,1991.[2]彥啟森,趙慶珠.建筑熱過程[M].北京:中國建筑工業出版社,1986.

[3]SY/T0511.1-2010,石油儲罐附件第一部分:呼吸閥[S].北京:石油工業出版社,2011.

[4]李征西,等.油品儲運設計手冊[M].北京:石油工業出版社,1997.

Simulation and Calculation of Gas-phase Spatial Pressure of Incoming and Outgoing Oil in Dome-roofed Oil Tanks

Meng Fanqin,Xi Libo,Yu Baijian
(Department of Aviation POL,Air Force Logistics University,Xuzhou 221000,China)

In this paper,we presented the simplified method for the simulation of the incoming and outgoing oil in a dome-roofed oil tank and then calculated the changes in the gas-phase spatial pressure of the oil tank under the normal working condition of the breather valve so as to provide reference for the safe operation and scientific management of the oil tank.

oil depot;dome-roofed oil tank;gas-phase spatial pressure;numerical simulation

F224.0;TE85

A

1005-152X(2017)09-0154-03

10.3969/j.issn.1005-152X.2017.09.035

2017-08-06

孟凡芹（1965-），男，安徽蕭縣人，博士，空軍勤務學院航空油料物資系儲運教研室教授，研究方向：儲運自動化、物流信息化。