孔祥宇, 焦文玲, 叢德文 , 王建華 , 馮立德
(1.哈爾濱工業(yè)大學建筑學院寒地城鄉(xiāng)人居環(huán)境科學與技術工業(yè)和信息化部重點實驗室,黑龍江哈爾濱150090;2.中國石油天然氣股份有限公司天然氣銷售廣東分公司,廣東廣州510220;3.中國石油天然氣股份有限公司天然氣銷售分公司,北京100101)
國際海事組織(International Maritime Organization,IMO)宣布于2020年1月1日起實施限硫令,要求全球范圍內船用燃油含硫量不得超過0.5%,并制定了到2050年將溫室氣體排放減半的階段目標,這標志著全球海運業(yè)正式進入“低硫、低碳”時代。液化天然氣作為船用可持續(xù)發(fā)展清潔能源,可各方面均滿足國際船舶排放控制要求,是實現(xiàn)“雙碳”目標的重要保障手段。國際航運權威機構預測,到2025年國際LNG動力船舶將達到3 000艘左右,全球船舶訂單中將有60.3%的新造船采用LNG動力[1-2]。
現(xiàn)階段為LNG動力船加注燃料主要方式為槽車加注、岸站加注、船對船加注。其中船對船加注可提供更為高效、及時和規(guī)模化的加注服務。目前,104箱級受注船、大型礦砂船、超大型油輪等大型遠洋船舶大多采用船對船加注。為了減少靠港口靠泊時間,需要在燃料加注的同時進行貨物的裝卸操作,即同步作業(yè)。一旦加注過程中發(fā)生LNG泄漏,產(chǎn)生可燃氣,氣體擴散至同步作業(yè)區(qū)域將會被作業(yè)所產(chǎn)生的靜電、火花等點燃,存在火災爆炸風險[3-4]。
以深圳鹽田國際港區(qū)16號碼頭為作業(yè)碼頭,以新奧普陀號8 500 m3LNG加注船(簡稱加注船)和法國達飛14800TEU系列LNG受注船(簡稱受注船)為研究對象,基于FLACS軟件求解加注船和受注船同步作業(yè)過程的泄漏后果,定量分析軟管破裂和法蘭接頭處泄漏等典型高風險場景,獲取蒸氣云在受注船壓載和滿載工況下的擴散范圍,據(jù)此劃定安全操作范圍,為加注船同步作業(yè)安全區(qū)域劃分提供重要的技術指導和數(shù)據(jù)支撐。
加注船與受注船的主要參數(shù)見表1?;诩幼⒋褪茏⒋闹饕獏?shù),采用FLACS軟件建立加注船、受注船三維模型,分別見圖1、2。圖2中受注船三維模型右側為船頭,左側為船尾,加注船與受注船同向,兩船凈距4 m。兩船間軟管連接三維模型見圖3。圖2中,BAY2—BAY86為受注船集裝箱具體位置,用貝位表示。其中,BAY2—BAY22為艏樓前貝位,BAY26為加注站所在貝位,BAY26—BAY70為艏樓艉樓間貝位,BAY74—BAY86為艉樓后貝位。鹽田港有大型集裝箱深水泊位16個,自大鵬灣口至鹽田港區(qū)進港航道全長約26 km,航道水深17.4 m,航道寬度400 m,基于FLACS建立鹽田港二期三維模型見圖4。
表1 加注船與受注船的主要參數(shù)
圖1 加注船三維模型(軟件截圖)
圖2 受注船三維模型(軟件截圖)
圖3 兩船間軟管連接三維模型(軟件截圖)
圖4 鹽田港二期三維模型(軟件截圖)
采用FLACS軟件對作業(yè)過程的泄漏過程氣相或液相可能擴散區(qū)域進行網(wǎng)格劃分,將泄漏區(qū)域作為核心區(qū)域,依據(jù)各工況泄漏位置分別設置核心區(qū)域、擴展區(qū)域。核心區(qū)域尺寸約為4 m×4 m×4 m,擴展區(qū)域尺寸約為600 m×100 m×80 m。為同時保證計算速度和計算精度,以核心區(qū)域的網(wǎng)格邊長作為變量,進行網(wǎng)格敏感性試算,結果表明核心區(qū)域的網(wǎng)格邊長為0.5 m時,網(wǎng)格符合。由于核心區(qū)域網(wǎng)格邊長均大于泄漏口半徑,故泄漏采用虛擬噴口技術以保證計算精度,核心區(qū)域與擴展區(qū)域采用Smooth功能過渡,比例系數(shù)為1.2,擴展區(qū)域最大網(wǎng)格邊長為5 m。網(wǎng)格劃分見圖5。
圖5 網(wǎng)格劃分(軟件截圖)
采用計算流體力學軟件FLACS,選取氣相軟管泄漏、液相軟管泄漏、加注船加注站法蘭接頭泄漏、受注船加注站法蘭接頭泄漏(液相)4種場景進行泄漏后果定量分析。分析受注船壓載工況(簡稱壓載工況)、受注船滿載工況(簡稱滿載工況)蒸氣云擴散范圍。
各泄漏場景的參數(shù)如下。
① 泄漏時間
氣相軟管、液相軟管及相關接頭規(guī)格均為DN 200 mm,傳輸壓力控制在1 MPa以內,應急切斷(Emergency Shut Down)系統(tǒng)響應時間取30 s,氣體泄漏探測時間取30 s,累計泄漏時間為60 s。
② 擴散邊界
依據(jù)GB/T 20368—2021《液化天然氣(LNG)生產(chǎn)、儲存和裝運》、 中國船級社《油氣定量風險評估指南(2020版)》、中國船級社《液化天然氣燃料加注作業(yè)指南(2021版)》、NFPA59A—2021《液化天然氣(LNG)生產(chǎn)、儲存和裝運標準》(Standard for the Production, Storage, and Handling of Liquefied Natural Gas )等,選取甲烷體積分數(shù)為2.5%作為蒸氣云擴散邊界。
③ 泄漏口直徑
根據(jù)GB/T 26610.4—2022《承壓設備系統(tǒng)基于風險的檢驗實施導則 第4部分:失效可能性定量分析方法》、英國健康安全執(zhí)行局(Health and Safety Executive)數(shù)據(jù)庫、國際油氣生產(chǎn)者協(xié)會(International Association of Oil &Gas Producers)數(shù)據(jù)庫等,軟管泄漏時,泄漏口直徑取軟管直徑的10%,即20 mm,法蘭接頭泄漏時,泄漏口直徑取30 mm。
父母們以為最大限度尊重孩子的意志,就是制造了民主的家庭氛圍,事實上很容易走向極端,制造出專制的氛圍。孩子成為本質上的獨裁者。孩子一哭,全家就得啟動緊張模式,分工合作滿足孩子當下需求;孩子要睡覺,全家就得保持安靜;孩子怕熱,全家就得跟著開空調,冷也得乖乖忍著;孩子愛看動畫片,全家就得低齡化,每天守著喜羊羊和灰太狼傻笑。
④ 泄漏工質組成
泄漏工質按純甲烷考慮。
⑤ 泄漏時環(huán)境情況
環(huán)境溫度15 ℃,環(huán)境壓力101 kPa,泄漏壓力1 MPa。模擬考慮風向為16個方向,根據(jù)風玫瑰圖設定風向。具體風向風速設定見表2。
氣相泄漏模擬時,考慮6個泄漏方向(垂直向上、垂直向下、水平向左、水平向右、垂直船側外板、負向垂直船側外板),同時考慮風向和泄漏方向后,以泄漏開始為0時刻,泄漏時間為60 s,共模擬96個工況,將96個工況進行疊加,得出蒸氣云最大擴散范圍。
液相泄漏模擬時,考慮1個泄漏方向(垂直向下)。同時考慮風向和泄漏方向后,以泄漏開始為0時刻,泄漏時間為60 s,共模擬16個工況,將16個工況進行疊加,得出蒸氣云最大擴散范圍。
表2 風向風速設定
分析泄漏后果時,4種場景均考慮風玫瑰圖中16個風向的風速。
氣相軟管泄漏位置見圖6。氣相軟管泄漏方向考慮相對于泄漏口處垂直向上、垂直向下、水平向左、水平向右、垂直船側外板、負向垂直船側外板6個泄漏方向。氣相軟管泄漏壓載工況、滿載工況模擬結果分別見圖7、8。
圖6 氣相軟管泄漏位置(軟件截圖)
圖7 氣相軟管泄漏壓載工況模擬結果(軟件截圖)
圖8 氣相軟管泄漏滿載工況模擬結果(軟件截圖)
由圖7可知,壓載工況下,由于受注船吃水較淺,兩船高差較大,蒸氣云一部分蔓延至加注船甲板上,少量蔓延至受注船甲板上,最深蔓延至受注船BAY26(加注站所在貝位)處向內6個箱位以及BAY30(加注站后1貝位)處向內3個箱位。蒸氣云最大擴散范圍為船長方向87 m,船寬方向31 m。
液相軟管泄漏位置見圖9。LNG泄漏方向為垂直向下,液相軟管泄漏壓載工況、滿載工況模擬結果分別見圖10、11。
圖9 液相軟管泄漏位置(軟件截圖)
圖10 液相軟管泄漏壓載工況模擬結果(軟件截圖)
圖11 液相軟管泄漏滿載工況模擬結果(軟件截圖)
由圖10可以看出,壓載工況下,由于受注船吃水較淺,兩船高差較大,故泄漏的蒸氣云全部在加注船側積聚,未蔓延至受注船甲板上。蒸氣云最大擴散范圍為船長方向230 m,船寬方向30 m。
由圖11可以看出,滿載工況下,液相軟管發(fā)生破損,泄漏的LNG會與水進行強烈的換熱,產(chǎn)生大量的天然氣,形成蒸氣云。蒸氣云沿著兩船夾道內擴散,并蔓延至加注船甲板上,未能波及受注船甲板上。蒸氣云最大擴散范圍為船長方向227 m,船寬方向30 m。
加注船加注站法蘭接頭泄漏位置見圖12。LNG泄漏方向為垂直向下。加注船加注站法蘭接頭泄漏壓載工況、滿載工況模擬結果分別見圖13、14。
圖12 加注船加注站法蘭接頭泄漏位置(軟件截圖)
圖13 加注船加注站法蘭接頭泄漏壓載工況模擬結果(軟件截圖)
圖14 加注船加注站法蘭接頭泄漏滿載工況模擬結果(軟件截圖)
由圖13可知,壓載工況下,由于受注船吃水較淺,兩船高差較大,蒸氣云全部在加注船側積聚,未蔓延至受注船甲板。蒸氣云最大擴散范圍為船長方向219 m,船寬方向25.5 m。
由圖14可知,滿載工況下,加注船加注站法蘭接頭發(fā)生泄漏,泄漏形成的蒸氣云會沿兩船夾道內擴散,且少量蔓延至加注船甲板上,未能波及受注船甲板上。蒸氣云最大擴散范圍為船長方向330 m,船寬方向25.5 m。
綜合對比壓載工況和滿載工況,加注船側蒸氣云擴散范圍存在顯著差異。這是由于兩船高差改變后,導致局部風場、風剪切(風速在空間中的分布)等發(fā)生了變化,從而影響了蒸氣云擴散范圍。
受注船加注站法蘭接頭泄漏位置見圖15。LNG泄漏方向為垂直向下。受注船加注站法蘭接頭泄漏壓載工況、滿載工況模擬結果分別見圖16、17。
圖15 受注船加注站法蘭接頭泄漏位置(軟件截圖)
圖16 受注船加注站法蘭接頭泄漏壓載工況模擬結果(軟件截圖)
圖17 受注船加注站法蘭接頭泄漏滿載工況模擬結果(軟件截圖)
由圖16可知,壓載工況下受注船加注站法蘭接頭泄漏后,泄漏的LNG會大量蒸發(fā)產(chǎn)生天然氣,形成蒸氣云。蒸氣云蔓延至受注船甲板,同時蔓延至兩船夾道內及加注船甲板上。蒸氣云最大擴散范圍為船長方向185 m,船寬方向76 m。
由圖17可知,滿載工況下受注船加注站法蘭接頭泄漏后,在受注船側,蒸氣云會蔓延至集裝箱甲板上,在加注船側,少量蔓延至加注船甲板上,其余蒸氣云蔓延至兩船夾道內及加注船船頭的水面方向。蒸氣云最大擴散范圍為船長方向223 m,船寬方向90 m。
① 氣相軟管泄漏受注船壓載工況,蒸氣云最大擴散范圍為87 m×31 m;受注船滿載工況,蒸氣云最大擴散范圍為127 m×21 m。
② 液相軟管泄漏受注船壓載工況,蒸氣云最大擴散范圍為230 m×30 m;受注船滿載工況,蒸氣云最大擴散范圍為227 m×30 m。
③ 加注船加注站法蘭接頭泄漏受注船壓載工況,蒸氣云最大擴散范圍為219.0 m×25.5 m;受注船滿載工況,蒸氣云最大擴散范圍為330.0 m×25.5 m。
④ 受注船加注站法蘭接頭泄漏受注船壓載工況,蒸氣云最大擴散范圍為185 m×76 m;受注船滿載工況,蒸氣云最大擴散范圍為223 m×90 m。