馬云修 李東澤

(1.中石化長輸油氣管道檢測(cè)有限公司，江蘇 徐州 221000；2.中國石油大學(xué)(華東)，山東 青島 266000)

1 概述

由于油氣回收能力高，浮頂罐在石油石化行業(yè)得到了廣泛應(yīng)用。但是近年來在浮頂油罐推廣的過程中暴露出各種安全隱患，其中油罐起火、爆炸事故頻繁出現(xiàn)[1-4]。油罐起火有三個(gè)必要條件：可燃物、助燃劑和火源。油品蒸氣和空氣混合后，可能形成爆炸性混合氣體，但是只有當(dāng)油品蒸氣在空氣中達(dá)到一定的體積濃度范圍，并遇明火時(shí)，才會(huì)發(fā)生爆炸。

為減少油品蒸發(fā)損耗，浮頂罐常采用兩次密封，這導(dǎo)致一二次密封空間容易存在油蒸氣與空氣的混合物，并很容易處于爆炸區(qū)間內(nèi)，此時(shí)罐內(nèi)若出現(xiàn)靜電、電位差、明火或者溫度過高等，極可能引發(fā)火災(zāi)爆炸事故。

目前國內(nèi)外有許多學(xué)者開展了加強(qiáng)密封的研究。Heisterberg[5]、Kinghorn[6]、李世超[7]等人通過改善機(jī)械結(jié)構(gòu)，增加了浮頂機(jī)械密封效果；李恩田[8-9]指出向浮頂油罐一二次密封間的油氣空間中加入氮?dú)?，可縮小浮頂油罐的一二次密封之間環(huán)形空間油氣爆炸極限范圍，可以有效防止爆炸事故的發(fā)生。

事實(shí)上，液態(tài)油表面油蒸氣濃度高(近似為100%)，大氣中油蒸氣濃度極低(近似為0%)，則必然存在油蒸氣濃度處于爆炸區(qū)間內(nèi)的某個(gè)區(qū)域，若該區(qū)域恰好存在氧氣，則很容易引發(fā)事故。使用氮?dú)獗Ｗo(hù)技術(shù)，將一二次密封體之間的油氣空間中氧氣驅(qū)替，可以降低“爆炸區(qū)間內(nèi)油蒸氣與氧氣共同存在”的可能性，極大程度避免爆炸事故的發(fā)生。考慮到氮?dú)饪赡苡捎诿芊馊Σ粐?yán)等因素發(fā)生逸出泄漏，對(duì)浮頂罐上操作人員的人身安全造成威脅，因此有必要對(duì)浮頂一二次密封空間內(nèi)氮?dú)馓畛?、逸出?guī)律開展研究[10]。

本文對(duì)浮頂一二次密封空間內(nèi)氮?dú)馓畛?、逸出?guī)律開展研究，利用Fluent數(shù)值模擬軟件，計(jì)算收發(fā)油時(shí)浮盤上下移動(dòng)對(duì)氮?dú)庖莩龅挠绊懀治龅獨(dú)獬涮钪脫Q過程中，一二次密封空間垂直方向上不同組分氣體濃度分布規(guī)律。

2 一二次密封空間氮?dú)鈹U(kuò)散規(guī)律數(shù)值仿真

2.1 浮頂移動(dòng)對(duì)充填氮?dú)庖莩龅挠绊?/h3>

當(dāng)一次氮?dú)獬涮钭鳂I(yè)結(jié)束后，氮?dú)庵魅肟谂c浮頂密封空間排氣口關(guān)閉，可以認(rèn)為氮?dú)庠诟№斨芟蛄鲃?dòng)擴(kuò)散較弱，一二次密封環(huán)形空間可以近似認(rèn)為中心對(duì)稱。利用Fluent軟件對(duì)一二次密封環(huán)形空間截面進(jìn)行建模，如圖1(a)所示。其中，下方矩形空間代表罐內(nèi)油氣空間，上方矩形空間代表大氣空間，中間梯形空間代表一二次密封空間。圖1(b)為計(jì)算結(jié)果，展示了不同浮盤移動(dòng)速度下密封空間內(nèi)氮?dú)鈹U(kuò)散規(guī)律?？梢钥闯?，不同浮盤移動(dòng)速度下，密封空間內(nèi)氮?dú)鉂舛茸兓?guī)律略有不同，但基本變化不大?？梢哉J(rèn)為，收發(fā)油產(chǎn)生的浮盤移動(dòng)對(duì)密封空間內(nèi)氮?dú)庖莩鏊俾驶緵]有影響。

圖1 利用Fluent軟件對(duì)一二次密封環(huán)形空間截面進(jìn)行的建模

邊界條件：下方矩形底邊為壓力入口邊界，上方矩形的頂邊為壓力出口邊界，圖中其他邊均為wall邊界，其中左側(cè)wall邊界可以設(shè)置為靜止或一定速度移動(dòng)。初始條件：下方矩形空間充滿油蒸氣，其組分如表1所示，上方矩形空間充滿空氣，中間梯形空間充滿氮?dú)狻?/p>

表1 氣體組分表

2.2 氮?dú)獬涮钸^程中空間垂直方向濃度分布

氮?dú)獬涮钸^程中，因組分密度差、重力協(xié)同作用，無法保證一二次密封空間截面內(nèi)濃度均一。所以，需要監(jiān)控一二次密封空間內(nèi)氮?dú)鉂舛茸畹偷奈恢?。?duì)于現(xiàn)場(chǎng)實(shí)際密封空間，設(shè)有氮?dú)庾⑷肟谂c氣體排出口，顯然排出口附近氮?dú)鉂舛茸畹?，故?duì)排出口附近垂直方向氮?dú)鉂舛冗M(jìn)行監(jiān)控。

以10000m3浮頂罐一二次密封油氣空間為原型，建立如圖2所示的仿真模型，其外環(huán)半徑為14.5m，內(nèi)環(huán)半徑為14.3m，高為1.5m，進(jìn)氣口、出氣口直徑均為0.1m。邊界條件：進(jìn)氣口為速度入口，其進(jìn)氣速度為1m/s；排氣口為壓力出口，其壓力設(shè)置為大氣壓。初始條件：假設(shè)初始時(shí)密封空間充滿油氣，以1m/s的速度向其中通入氮?dú)狻?/p>

圖2 一二次密封環(huán)形空間三維模型

采集注氣開始4h后的出口位置處濃度數(shù)據(jù)，如圖3所示?？梢钥闯?，氮?dú)鉂舛入S空間高度增加而降低，油蒸氣隨著空間高度的增加而增加，這是由于密度差與重力作用造成的。大約4.2h后，空間最低氮?dú)鉂舛冗_(dá)到99%，氮?dú)庵脫Q油氣基本完成。

圖3 不同時(shí)刻排氣口氣體濃度隨垂直高度變化規(guī)律

3 結(jié)語

本文利用Fluent軟件，建立一二次密封環(huán)形空間截面數(shù)值模型，計(jì)算空間內(nèi)氮?dú)鉂舛入S時(shí)間的變化規(guī)律，得到的具體結(jié)論如下：

(1)浮盤移動(dòng)速度對(duì)氮?dú)庖莩鏊俾视绊懖淮螅?/p>

(2)充填氮?dú)膺^程中，一二次密封空間的垂直方向上，氮?dú)鉂舛入S高度增加而降低，油蒸氣反之；

(3)以1m/s流速注入氮?dú)猓蠹s4.2h后，空間最低氮?dú)鉂舛冗_(dá)到99%，氮?dú)庵脫Q油氣基本完成。