• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    大型LNG船發(fā)電機(jī)室的燃?xì)夤芫€泄漏分析

    2017-10-13 11:55:01岑卓倫劉濤王磊竇旭
    中國(guó)艦船研究 2017年5期
    關(guān)鍵詞:艙室燃?xì)?/a>發(fā)電機(jī)

    岑卓倫,劉濤,王磊,竇旭

    滬東中華造船(集團(tuán))有限公司開(kāi)發(fā)所,上海 200129

    大型LNG船發(fā)電機(jī)室的燃?xì)夤芫€泄漏分析

    岑卓倫,劉濤,王磊,竇旭

    滬東中華造船(集團(tuán))有限公司開(kāi)發(fā)所,上海 200129

    [目的]目前,由雙燃料發(fā)動(dòng)機(jī)組成的電力推進(jìn)系統(tǒng)是大型液化天然氣(LNG)船的主流推進(jìn)方式,必須對(duì)爆炸性可燃?xì)怏w進(jìn)行安全可靠性的定性、定量評(píng)估,以規(guī)避潛在風(fēng)險(xiǎn)。[方法]以某雙燃料電力推進(jìn)大型LNG船發(fā)電機(jī)室為研究對(duì)象,對(duì)其內(nèi)部不同區(qū)域的燃?xì)猓ㄌ烊粴猓┬孤┕r進(jìn)行模擬分析。根據(jù)泄漏發(fā)生的形式、位置和速率等定義危險(xiǎn)泄漏工況,選擇雷諾應(yīng)力模型為湍流模型,采用計(jì)算流體力學(xué)(CFD)軟件Fluent對(duì)發(fā)電機(jī)室燃?xì)夤?yīng)管線的5個(gè)泄漏點(diǎn)進(jìn)行持續(xù)泄漏模擬計(jì)算,并將泄漏擴(kuò)散結(jié)果與艙室通風(fēng)的流場(chǎng)速度分布相結(jié)合,得到不同區(qū)域發(fā)生泄漏后的天然氣擴(kuò)散趨勢(shì)和濃度分布。[結(jié)果]根據(jù)仿真模擬結(jié)果優(yōu)化了可燃?xì)怏w探測(cè)器布置方案,并明確了排氣風(fēng)機(jī)無(wú)需進(jìn)行防爆設(shè)計(jì)。[結(jié)論]研究結(jié)果可為有限空間內(nèi)通風(fēng)條件下的可燃?xì)怏w泄漏事故分析防范提供參考,并且適用于燃燒爆炸破壞的定量評(píng)估,用以指導(dǎo)結(jié)構(gòu)強(qiáng)度設(shè)計(jì)。

    液化天然氣;管道泄漏;計(jì)算流體力學(xué);定量風(fēng)險(xiǎn)評(píng)估

    Abstract:[Objectives]The electric propulsion dual-fuel engine is becoming dominant in newly built Liquefied Natural Gas(LNG)carriers.To avoid the potential risks that accompany the use of flammable and explosive boil-off gas,the performance of precise safety and reliability assessments is indispensable.[Methods]This research concerns the engine rooms of large LNG carriers which are propelled electrically by a dual-fuel engine.Possible fuel gas(natural gas)leak cases in different areas of the engine room are simulated and analyzed.Five representative leak cases defined by leak form,leak location and leak rate are entered into a Computational Fluid Dynamics(CFD)simulation,in which the Reynolds stress model of Fluent software is adopted as the turbulence model.The results of the leaked gas distribution and ventilation velocity field are analyzed in combination to obtain the diffusion tendency and concentration distribution of leaked gas in different areas.[Results]Based on an analysis of the results,an optimized arrangement of flammable gas detectors is provided for the engine room,and the adoption of an explosion-proof exhaust fan is proven to be unnecessary.[Conclusions]These analysis methods can provide a reference for similar gas leakage scenarios occurring in confined ventilated spaces.In addition,the simulation results can be used to quantitatively assess potential fire or explosion damage in order to guide the design of structural reinforcements.

    Key words:Liquefied Natural Gas(LNG);pipe leakage;Computational Fluid Dynamics(CFD);quantitative risk assessment

    0 引 言

    隨著環(huán)保規(guī)范的完善和排放標(biāo)準(zhǔn)的提高,以液化天然氣(Liquefied Natural Gas,LNG)等清潔氣體為動(dòng)力燃料的船舶得以迅速發(fā)展[1-3]。由于天然氣的成本優(yōu)勢(shì)以及發(fā)達(dá)國(guó)家對(duì)其配套基礎(chǔ)設(shè)施的不斷完善,它將成為未來(lái)航運(yùn)的主流燃料[4-7]。但天然氣供氣管線連接件密封不嚴(yán)、管道振動(dòng)、冷熱交替等造成的疲勞損傷,或高含硫燃料的酸性腐蝕,乃至人為操作不當(dāng)?shù)炔豢深A(yù)見(jiàn)因素,均有可能造成天然氣泄漏。天然氣無(wú)色無(wú)味,發(fā)生泄漏后難以察覺(jué),長(zhǎng)時(shí)間累積后若遇火花則存在燃爆風(fēng)險(xiǎn)[8-9]。當(dāng)天然氣摩爾濃度大于4.9%時(shí)即有燃爆風(fēng)險(xiǎn)(按純甲烷計(jì)算),過(guò)高濃度還可能造成人員窒息,因此,需要準(zhǔn)確評(píng)估和有效防范天然氣的使用風(fēng)險(xiǎn)。

    在氣體流動(dòng)擴(kuò)散的場(chǎng)分布定量計(jì)算方面,計(jì)算流體力學(xué)(Computational Fluid Dynamics,CFD)方法能夠模擬由流體的可壓縮性、湍流流動(dòng)、密度差等引起的自然流動(dòng)、混合物熱力學(xué)運(yùn)動(dòng)和分子擴(kuò)散等物理效應(yīng),適于模擬驗(yàn)證復(fù)雜湍流流動(dòng)的氣體擴(kuò)散過(guò)程[10-11],以及易燃、易爆氣體的泄漏擴(kuò)散風(fēng)險(xiǎn)評(píng)估。

    2006年,DNV[12]利用FLACS和KFX軟件對(duì)MAN公司某型主機(jī)高壓供氣系統(tǒng)在機(jī)艙內(nèi)的泄漏燃爆風(fēng)險(xiǎn)進(jìn)行了危險(xiǎn)性評(píng)估。2011年,程浩力等[13]對(duì)城市街道的狹窄有限空間內(nèi)不同風(fēng)力條件下的天然氣管線擴(kuò)散泄漏進(jìn)行了CFD分析。2014年,王麗華[14]利用Fluent軟件分析了風(fēng)速、風(fēng)向、泄漏方向?qū)Q笃脚_(tái)油氣泄漏擴(kuò)散的影響。2016年,F(xiàn)u等[15]將事件樹(shù)分析法與CFD方法相結(jié)合,研究了LNG船泄漏事故的發(fā)生概率并定量評(píng)估了泄漏事故的嚴(yán)重程度。同年,F(xiàn)iates等[16]利用OpenFOAM軟件對(duì)甲烷和二氧化碳?xì)怏w在海洋平臺(tái)和風(fēng)洞內(nèi)的泄漏擴(kuò)散過(guò)程進(jìn)行了模擬分析。IMO[17]針對(duì)某大型礦砂船的LNG燃料預(yù)處理房間,采用事件樹(shù)分析法和故障樹(shù)分析法評(píng)估了潛在風(fēng)險(xiǎn),并借助CFD方法對(duì)燃爆事故后果進(jìn)行了定量評(píng)估。然而,鮮有通過(guò)泄漏事故模擬計(jì)算來(lái)指導(dǎo)預(yù)防機(jī)制(如可燃?xì)怏w探測(cè)器布置)和設(shè)備防爆設(shè)計(jì)等方面的研究。

    本文將針對(duì)某新型1.74×105m3LNG運(yùn)輸船發(fā)電機(jī)室的泄漏事故進(jìn)行模擬計(jì)算,并進(jìn)行可燃?xì)怏w探測(cè)器布置和排氣風(fēng)機(jī)的防爆需求估算。發(fā)電機(jī)室作為核心動(dòng)力單元,存在燃?xì)庑孤⑷急仁鹿曙L(fēng)險(xiǎn),而發(fā)生燃爆的可能性、劇烈程度與燃?xì)庑孤┝?、燃?xì)鉂舛让芮邢嚓P(guān)。因此,將借助商用CFD軟件Fluent計(jì)算分析艙內(nèi)氣體流動(dòng)和燃?xì)庑孤U(kuò)散情況:首先,分析通風(fēng)流場(chǎng)特性,并對(duì)發(fā)電機(jī)室的燃?xì)夤?yīng)管線進(jìn)行分區(qū);然后,在不同區(qū)域內(nèi)選取典型泄漏點(diǎn),計(jì)算分析泄漏后的天然氣氣云(以下簡(jiǎn)稱“氣云”)的擴(kuò)散趨勢(shì)和濃度分布;最后,合理布置可燃?xì)怏w探測(cè)器,最大限度地降低燃爆風(fēng)險(xiǎn)。

    1 建模與計(jì)算

    1.1 計(jì)算方法

    對(duì)于發(fā)電機(jī)室的內(nèi)部空氣流動(dòng),將采用非定常雷諾平均N-S方程(Navier-Stokes equation)求解。對(duì)于泄漏天然氣的擴(kuò)散流動(dòng),將采用組分輸運(yùn)方程求解。在xyz直角坐標(biāo)系下,時(shí)均后的控制方程為[18-19]

    式中:Φ和Φ'分別為流體某一時(shí)均變量和相應(yīng)的脈動(dòng)量;t為時(shí)間;U為時(shí)均速度矢量,其在直角坐標(biāo)系有3個(gè)方向的速度分量u,v,w;u',v',w'為u,v,w相應(yīng)的脈動(dòng)量;為廣義擴(kuò)散系數(shù);SΦ為廣義源項(xiàng)。

    發(fā)電機(jī)室通風(fēng)系統(tǒng)有21個(gè)不同尺寸的進(jìn)風(fēng)口。將入口進(jìn)風(fēng)視為射流,不同入射速度的進(jìn)風(fēng)會(huì)在流場(chǎng)內(nèi)引起大范圍摻混和動(dòng)量交換,同時(shí)發(fā)電機(jī)室內(nèi)多個(gè)不同形狀的大尺度設(shè)備和結(jié)構(gòu)平臺(tái)將在背風(fēng)側(cè)形成大面積的尾渦。由于整個(gè)發(fā)電機(jī)室存在大量的剪切流動(dòng),故在計(jì)算中選取雷諾應(yīng)力模型(Reynolds Stress Model,RSM)作為湍流模型來(lái)封閉N-S動(dòng)量方程中的雷諾應(yīng)力項(xiàng)。與兩方程k-ε湍流模型相比,RSM通過(guò)直接對(duì)雷諾應(yīng)力分量建立6個(gè)輸運(yùn)方程,附加一個(gè)耗散方程進(jìn)行求解[19],而放棄了利用雷諾應(yīng)力正比于時(shí)均速度的應(yīng)變和各向同性的湍流粘性系數(shù)這2個(gè)假設(shè)來(lái)計(jì)算湍流應(yīng)力的方法。

    在預(yù)測(cè)有限空間內(nèi)非穩(wěn)態(tài)湍流流動(dòng)的擴(kuò)散過(guò)程方面,RSM的應(yīng)用非常廣泛,Cehlin等[20]研究了氣體分布器對(duì)房間內(nèi)空氣流動(dòng)的影響,其利用RSM所得的房間內(nèi)空氣速度分布穩(wěn)態(tài)模擬結(jié)果與試驗(yàn)結(jié)果吻合較好,并指出非穩(wěn)態(tài)模擬的結(jié)果將更加準(zhǔn)確。郭棟鵬等[21]應(yīng)用RNG k-ε模型和RSM計(jì)算了建筑物對(duì)污染物擴(kuò)散的影響,其中RSM對(duì)污染物濃度分布和周圍流場(chǎng)變化的預(yù)測(cè)效果更好。因此,本文選擇RSM作為湍流模型,并設(shè)置與文獻(xiàn)[20]和文獻(xiàn)[21]相似的邊界條件,對(duì)通風(fēng)條件下發(fā)電機(jī)室的燃?xì)庑孤U(kuò)散進(jìn)行非穩(wěn)態(tài)模擬計(jì)算。

    1.2 發(fā)電機(jī)室建模

    1.2.1 供氣系統(tǒng)和可燃?xì)怏w探測(cè)系統(tǒng)

    主機(jī)供氣系統(tǒng)的原理如圖1所示。原動(dòng)機(jī)布置于發(fā)電機(jī)室內(nèi),采用MAN DIESEL雙燃料中速機(jī),航行過(guò)程中主要以LNG液貨艙揮發(fā)氣體作為燃料。

    圖1 供氣系統(tǒng)原理圖Fig.1 Fuel-gas supply system diagram

    發(fā)電機(jī)室作為動(dòng)力系統(tǒng)的核心單元,除發(fā)電機(jī)組之外,還布置了輸送泵、風(fēng)機(jī)、行車等電動(dòng)設(shè)備,這些設(shè)備的劇烈振動(dòng)將增加燃?xì)庑孤┖腿急娘L(fēng)險(xiǎn)。一旦發(fā)生泄漏,外部可燃?xì)怏w探測(cè)器及時(shí)有效地反饋并啟動(dòng)切斷機(jī)制是最后一道安全屏障。該船采用固定式可燃?xì)怏w探測(cè)器,當(dāng)任意一個(gè)傳感器達(dá)到20%最低爆炸下限即報(bào)警并啟動(dòng)預(yù)防機(jī)制。因此,可燃?xì)怏w探測(cè)器的合理布放是探測(cè)系統(tǒng)及時(shí)準(zhǔn)確預(yù)警的基礎(chǔ),應(yīng)在安裝前評(píng)估布置方案。按照傳統(tǒng)的布置方案,一般在發(fā)動(dòng)機(jī)上方空間的4個(gè)角落各布置一個(gè)探測(cè)器,如圖2所示。

    圖2 發(fā)電機(jī)室內(nèi)可燃?xì)怏w探測(cè)器的傳統(tǒng)布置Fig.2 Traditional arrangement of flammable gas detectors for engine room

    1.2.2 發(fā)電機(jī)室?guī)缀谓?/h4>

    圖3所示為發(fā)電機(jī)室的三維模型,艙室縱向長(zhǎng)度約22.4 m,橫向?qū)挾燃s20.65 m,高度約10.14 m。如圖3(a)所示,艙內(nèi)布置了1臺(tái)12V雙燃料發(fā)動(dòng)機(jī)、1臺(tái)8L雙燃料發(fā)動(dòng)機(jī)、3條風(fēng)管管線,艙室中部設(shè)有一層結(jié)構(gòu)平臺(tái)。如圖3(b)所示,艙室右后側(cè)(艉向右舷側(cè))區(qū)域由船體外板封閉。

    圖3 發(fā)電機(jī)室的三維模型Fig.3 Three-dimensional model of engine room

    將艙室內(nèi)布置燃?xì)夤?yīng)管線的區(qū)域劃分為3個(gè)部分,如圖4所示。

    圖4 燃?xì)夤?yīng)管線區(qū)域劃分圖Fig.4 Divided partitions of fuel-gas supply pipeline in engine room

    根據(jù)發(fā)電機(jī)室的三維模型(圖3)建立CFD幾何模型,如圖5所示。本文采用統(tǒng)一的全局坐標(biāo):x軸沿船長(zhǎng)方向,船艏為正;y軸沿船寬方向,左舷為正;z軸沿船高方向,上為正。

    圖5 CFD幾何模型Fig.5 Geometry built for CFD analysis

    1.3 網(wǎng)格劃分

    對(duì)發(fā)電機(jī)室上方無(wú)障礙物的開(kāi)闊空間采用結(jié)構(gòu)化六面體網(wǎng)格,其他區(qū)域采用非結(jié)構(gòu)化四面體網(wǎng)格,以提高幾何表面的建模精度。綜合考慮精度要求和仿真時(shí)間,全局網(wǎng)格尺寸設(shè)定為0.3 m。在泄漏源處、進(jìn)出風(fēng)口處等物理參數(shù)梯度變化較大的區(qū)域,以及發(fā)動(dòng)機(jī)上方通風(fēng)匯流區(qū)域進(jìn)行網(wǎng)格加密,而遠(yuǎn)離泄漏源和進(jìn)出風(fēng)口的位置則采用全局設(shè)置以縮減仿真時(shí)間。網(wǎng)格單元的最大偏斜度(Skewness)小于0.85,其中90%以上網(wǎng)格單元的最大偏斜度均小于0.6。網(wǎng)格單元總數(shù)量在6× 105~1.20×106范圍內(nèi)(不同計(jì)算工況的網(wǎng)格數(shù)量不同),整體計(jì)算域的網(wǎng)格劃分示例如圖6所示。

    圖6 網(wǎng)格劃分示例Fig.6 An example of meshing

    1.4 泄漏工況定義

    根據(jù)發(fā)電機(jī)室的總體布置、燃?xì)夤?yīng)管線的位置、通風(fēng)流動(dòng)模擬結(jié)果和泄漏形式,定義5個(gè)泄漏工況,如表1所示。

    表1 泄漏工況定義Table1 Definition of leak cases

    1.4.1 泄漏點(diǎn)選取

    基于燃?xì)夤芫€劃分區(qū)域(圖4)、艙室內(nèi)空氣流動(dòng)趨勢(shì)和易發(fā)生泄漏處這3個(gè)方面,在3個(gè)燃?xì)夤芫€區(qū)域內(nèi)共選擇了5處泄漏點(diǎn),如圖7所示。

    由于設(shè)備振動(dòng)、連接不牢等原因,燃?xì)夤?yīng)管線與發(fā)動(dòng)機(jī)的柔性軟管連接處最容易成為泄漏點(diǎn),因此分別選取8L和12V發(fā)動(dòng)機(jī)的柔性軟管連接處作為1號(hào)和3號(hào)泄漏點(diǎn),如圖7(a)和圖7(b)所示,即Case 1和Case 3模擬計(jì)算工況。依據(jù)燃?xì)夤芫€區(qū)域內(nèi)的通風(fēng)流場(chǎng),分別選取8L和12V發(fā)動(dòng)機(jī)頂部燃?xì)夤苣┒俗鳛?號(hào)和4號(hào)泄漏點(diǎn),如圖7(a)和圖7(b)所示,即Case 2和Case 4模擬計(jì)算工況。并選取12V發(fā)動(dòng)機(jī)頂部燃?xì)夤苤虚g段作為5號(hào)泄漏點(diǎn),如圖7(c)所示,即Case 5模擬計(jì)算工況。

    圖7 泄漏點(diǎn)位置Fig.7 Locations of leak points

    1.4.2 泄漏條件

    管路損壞泄漏分為完全破損和小尺寸破損2種情況。供氣管路完全破損會(huì)導(dǎo)致其內(nèi)部壓力驟降,進(jìn)而觸發(fā)燃?xì)鈮翰顖?bào)警器并關(guān)閉燃?xì)忾y組[22],所以不會(huì)造成燃?xì)獾拇竺娣e泄漏和長(zhǎng)時(shí)間累積。而小尺寸破損一般不足以觸發(fā)燃?xì)鈮翰顖?bào)警器,即會(huì)形成可燃?xì)怏w的持續(xù)泄漏和累積。無(wú)色無(wú)味的天然氣難以及時(shí)察覺(jué),需要通過(guò)可燃?xì)怏w探測(cè)器及時(shí)預(yù)警,故研究對(duì)象為管路小尺寸破損的泄漏條件。

    設(shè)管道內(nèi)徑為D,破口尺寸為d,當(dāng)d/D≤0.2時(shí)可以采用小尺寸泄漏模型來(lái)計(jì)算氣體的泄漏速率。小尺寸泄漏模型將管道容積視為足夠大,假定管內(nèi)壓力不受泄漏影響并保持恒定[23]。將小尺寸泄漏視為絕熱過(guò)程,天然氣視為純甲烷理想氣體,利用伯努利方程和絕熱方程,同時(shí)考慮氣體實(shí)際泄漏時(shí)的摩擦損失,則氣體的泄漏速度v0為[24]

    式中:φ為孔口流速系數(shù),即實(shí)際流速與理論流速之比,一般取為0.97~0.98;p1為供氣管內(nèi)壓強(qiáng),Pa;p0為供氣管外環(huán)境壓強(qiáng),Pa;k為絕熱指數(shù),天然氣一般取為1.3;R=8.314,為摩爾氣體常數(shù),J/(mol·K);T為氣體溫度,K。

    假定泄漏發(fā)生時(shí)管內(nèi)壓強(qiáng)為2.83×105Pa(設(shè)計(jì)運(yùn)行工況)并保持穩(wěn)定,則根據(jù)式(2)得到出口泄漏流速為66.29 m/s。設(shè)供氣管線內(nèi)徑為102 mm,破損開(kāi)口為20 mm×20 mm的方形開(kāi)口。

    1.5 求解及邊界條件設(shè)置

    發(fā)電機(jī)室頂部有一臺(tái)排風(fēng)風(fēng)機(jī),其通風(fēng)系統(tǒng)包含3條風(fēng)管和21個(gè)進(jìn)風(fēng)口。相對(duì)于艙室內(nèi)部空間而言,2臺(tái)雙燃料發(fā)動(dòng)機(jī)的渦輪增壓器空氣入口也是排風(fēng)出口。設(shè)環(huán)境壓力為101 325 Pa,壁面設(shè)置為絕熱。風(fēng)口氣體入流速度均小于0.3倍聲速,將空氣和天然氣視為不可壓縮的理想氣體,艙室內(nèi)流場(chǎng)為非定常流動(dòng)。邊界條件設(shè)置如下:

    1)風(fēng)管送風(fēng)入口。對(duì)于21個(gè)進(jìn)風(fēng)口,采用質(zhì)量流量入口邊界,并依據(jù)實(shí)船設(shè)計(jì)參數(shù)給定質(zhì)量流量(0.92~25.52 kg/s)。入流方向均垂直于入口截面,環(huán)境溫度為300 K。

    2)艙室頂部風(fēng)機(jī)排風(fēng)出口。采用風(fēng)扇邊界,并依據(jù)運(yùn)行工況將質(zhì)量流量設(shè)為30.625 kg/s。

    3)渦輪增壓器排風(fēng)出口。采用壓力出口邊界,出口壓力設(shè)為0 Pa(表壓)。

    選取壓力耦合求解器和隱式控制方程組求解非定常的質(zhì)量、動(dòng)量、能量和組分輸運(yùn)方程,其中壓強(qiáng)—速度關(guān)聯(lián)形式采用PISO格式。為保證精度,采用二階迎風(fēng)格式對(duì)動(dòng)量、湍流動(dòng)能和湍流耗散率等進(jìn)行數(shù)值離散。計(jì)算過(guò)程中,將連續(xù)方程的迭代殘余量收斂標(biāo)準(zhǔn)設(shè)為1×10-4,組分方程設(shè)為5×10-4,其余控制方程為1×10-3。

    2 計(jì)算結(jié)果與分析

    2.1 通風(fēng)計(jì)算結(jié)果

    由于艙內(nèi)燃?xì)庑孤┌l(fā)生于通風(fēng)條件下,因此艙內(nèi)非穩(wěn)態(tài)通風(fēng)流場(chǎng)模擬計(jì)算的時(shí)間步長(zhǎng)取0.02 s,將通風(fēng)場(chǎng)計(jì)算結(jié)果作為泄漏擴(kuò)散計(jì)算的初場(chǎng)。由于天然氣的比重小于1,其發(fā)生泄漏后是向上飄散,同時(shí)可燃?xì)怏w探測(cè)器布置于艙室上方,故選取z=8 m處的水平截面作為分析截面,如圖8(a)所示。圖8(b)~圖8(h)所示為該水平截面內(nèi)不同時(shí)刻的速度矢量分布圖,其中vmax=3.5 m/s。

    由圖8可知,當(dāng)t=1.2~5.2 s時(shí),處于通風(fēng)初始階段,截面內(nèi)流動(dòng)趨勢(shì)尚不明顯,流場(chǎng)內(nèi)主體范圍的流速在0.1 m/s以下;當(dāng)t=5.2~19.9 s時(shí),截面內(nèi)各點(diǎn)速度開(kāi)始增大,主體范圍的流速達(dá)到1 m/s以上;當(dāng)t=19.9~51.3 s時(shí),截面內(nèi)局部速度大小和方向的變化開(kāi)始減??;在t=51.3~91.3 s時(shí),流場(chǎng)內(nèi)主體范圍的流速大小和方向已無(wú)明顯變化,其中流動(dòng)趨勢(shì)較明顯(流速大于1 m/s)區(qū)域的流場(chǎng)已趨于相同,差異主要集中在流動(dòng)較弱的區(qū)域(流速小于0.5 m/s),故將t=91.3 s時(shí)刻的通風(fēng)場(chǎng)計(jì)算結(jié)果作為泄漏擴(kuò)散計(jì)算的初場(chǎng)。

    圖8 z=8 m截面在不同時(shí)刻的速度矢量分布圖Fig.8 Velocity vectors inz=8 m section at different times

    為進(jìn)一步分析截面內(nèi)流場(chǎng),將趨于穩(wěn)定后的流場(chǎng)(t=91.3 s)主體流動(dòng)趨勢(shì)以紅色箭頭標(biāo)出,如圖9所示。當(dāng)t≥51.3s時(shí),截面內(nèi)的主體流動(dòng)趨勢(shì)均與圖9相似。圖9中,箭頭1,2所示區(qū)域即圖4中燃?xì)夤芫€區(qū)域1上方,整體流動(dòng)先沿+x向,接近端部后逐漸轉(zhuǎn)向-y向;箭頭3所示區(qū)域即燃?xì)夤芫€區(qū)域2上方,整體流動(dòng)沿-x向;對(duì)于燃?xì)夤芫€區(qū)域3上方,流場(chǎng)各點(diǎn)處速度方向較為分散,無(wú)明顯的一致流向,僅在右側(cè)區(qū)域呈現(xiàn)出較弱的+y向匯流趨勢(shì),同時(shí)區(qū)域3上方整體流速均較小,低于1 m/s。

    根據(jù)以上分析可知不同區(qū)域發(fā)生燃?xì)庑孤┖蟮臍饬髦黧w流動(dòng)趨勢(shì),但通過(guò)單一截面的速度分布還不足以準(zhǔn)確判斷泄漏后可燃?xì)怏w的擴(kuò)散趨勢(shì)。具體的擴(kuò)散分布除受當(dāng)?shù)仫L(fēng)速和風(fēng)向影響外,還受射流速度、整體空間流動(dòng)、設(shè)備阻礙及密度差等因素的影響。

    圖9 z=8 m截面內(nèi)主體流動(dòng)趨勢(shì)圖(t=91.3 s)Fig.9 Main flow tendency in z=8 m section at t=91.3 s

    2.2 泄漏擴(kuò)散計(jì)算結(jié)果

    2.2.1 Case 1

    1號(hào)泄漏點(diǎn)發(fā)生天然氣泄漏后,不同時(shí)刻的天然氣摩爾濃度分布如圖10所示。天然氣在初始動(dòng)量和密度差作用下向頂部流動(dòng),同時(shí)在分子熱運(yùn)動(dòng)和氣流湍動(dòng)的疊加作用下向四周擴(kuò)散,難以到達(dá)艙室頂部。隨著時(shí)間的延長(zhǎng),擴(kuò)散范圍不斷增大,氣云的高濃度區(qū)域集中在泄漏口上方。

    圖10 1號(hào)泄漏點(diǎn)發(fā)生連續(xù)泄漏后發(fā)電機(jī)室不同時(shí)刻的天然氣摩爾濃度分布圖Fig.10 Distributions of natural gas mole concentration in engine room at different times for the continuous leak of Case 1

    2.2.2 Case 2

    2號(hào)泄漏點(diǎn)發(fā)生天然氣泄漏后,不同時(shí)刻的天然氣摩爾濃度分布如圖11所示。天然氣向上流動(dòng)并向四周擴(kuò)散,與1號(hào)泄漏點(diǎn)相似。但天然氣的主流方向比較明顯,主要流向艙室中部區(qū)域,并在到達(dá)頂部后開(kāi)始向四周大范圍擴(kuò)散。整個(gè)泄漏過(guò)程中,氣云的高濃度區(qū)域均集中在艙室中部。

    圖11 2號(hào)泄漏點(diǎn)發(fā)生連續(xù)泄漏后發(fā)電機(jī)室不同時(shí)刻的天然氣摩爾濃度分布圖Fig.11 Distributions of natural gas mole concentration in engine room at different times for the continuous leak of Case 2

    Case 1和Case 2的泄漏位置均位于燃?xì)夤芫€區(qū)域3(圖4),結(jié)合z=8 m截面(圖9)的通風(fēng)場(chǎng)流動(dòng)趨勢(shì)可知:當(dāng)區(qū)域3左側(cè)(1號(hào)泄漏點(diǎn)附近)發(fā)生泄漏時(shí),天然氣會(huì)在短時(shí)間內(nèi)向四周擴(kuò)散,氣云的高濃度區(qū)域位于泄漏點(diǎn)上方;當(dāng)區(qū)域3右側(cè)(2號(hào)泄漏點(diǎn)附近)發(fā)生泄漏時(shí),天然氣主體會(huì)向艙室中部流動(dòng),氣云的高濃度區(qū)域位于艙室中上部。

    2.2.3 Case 3

    3號(hào)泄漏點(diǎn)發(fā)生天然氣泄漏后,不同時(shí)刻的天然氣摩爾濃度分布如圖12所示。天然氣主體流向艙室右上方(+x,+z方向),并在到達(dá)艙壁后開(kāi)始向-y方向擴(kuò)散,氣云的高濃度區(qū)域集中在泄漏點(diǎn)右上方。

    圖12 3號(hào)泄漏點(diǎn)發(fā)生連續(xù)泄漏后發(fā)電機(jī)室不同時(shí)刻的天然氣摩爾濃度分布圖Fig.12 Distributions of natural gas mole concentration in engine room at different times for the continuous leak of Case 3

    2.2.4 Case 4

    4號(hào)泄漏點(diǎn)發(fā)生天然氣泄漏后,不同時(shí)刻的天然氣摩爾濃度分布如圖13所示。天然氣主體的流動(dòng)趨勢(shì)與3號(hào)泄漏點(diǎn)相似,流向艙室右上方,到達(dá)艙壁后再轉(zhuǎn)向-y方向,氣云的高濃度區(qū)域同樣集中在泄漏點(diǎn)右上方。

    圖13 4號(hào)泄漏點(diǎn)發(fā)生連續(xù)泄漏后發(fā)電機(jī)室不同時(shí)刻的天然氣摩爾濃度分布圖Fig.13 Distributions of natural gas mole concentration in engine room at different times for the continuous leak of Case 4

    Case 3和Case 4的泄漏點(diǎn)位置均位于燃?xì)夤芫€區(qū)域1(圖4),結(jié)合z=8 m截面(圖9)中的通風(fēng)場(chǎng)流動(dòng)趨勢(shì)可知,在區(qū)域1內(nèi)發(fā)生泄漏,天然氣主體將隨通風(fēng)氣流向右側(cè)流動(dòng)并逐漸擴(kuò)散,氣云的高濃度區(qū)域位于泄漏點(diǎn)右上方。

    2.2.5 Case 5

    5號(hào)泄漏點(diǎn)發(fā)生天然氣泄漏后,不同時(shí)刻的天然氣摩爾濃度分布如圖14所示。天然氣主體較為一致的流向左上方(-x,+z方向),與z=8 m截面的通風(fēng)場(chǎng)流動(dòng)趨勢(shì)一致,但排氣管的空間阻礙會(huì)引起分流。Case 5中天然氣摩爾濃度大于0.1%的氣云總體積和高度均明顯小于前4個(gè)泄漏點(diǎn),因?yàn)镃ase 5的泄漏點(diǎn)位于燃?xì)夤芫€區(qū)域2,該區(qū)域內(nèi)的通風(fēng)作用明顯強(qiáng)于其他區(qū)域(圖9),天然氣泄漏后即被吹散并最終由處于下風(fēng)向的頂部排氣扇排出,故泄漏發(fā)生后的8~118 s內(nèi)氣云濃度并無(wú)明顯變化。該區(qū)域內(nèi)發(fā)生泄漏后,天然氣主體都會(huì)隨艙內(nèi)通風(fēng)向左側(cè)(-x)流動(dòng),氣云的高濃度區(qū)域位于泄漏點(diǎn)左上方。

    圖14 5號(hào)泄漏點(diǎn)發(fā)生連續(xù)泄漏后發(fā)電機(jī)室不同時(shí)刻的天然氣摩爾濃度分布圖Fig.14 Distributions of natural gas mole concentration in engine room at different times for the continuous leak of Case 5

    2.3 可燃?xì)怏w探測(cè)器布置優(yōu)化

    為了提高探測(cè)精度,可燃?xì)怏w探測(cè)器應(yīng)按天然氣擴(kuò)散趨勢(shì)布置于最早出現(xiàn)高濃度氣云的區(qū)域內(nèi)。艙室內(nèi)泄漏的天然氣氣云主體的擴(kuò)散方向主要受艙室內(nèi)通風(fēng)流動(dòng)的影響,但天然氣的濃度梯度分布仍依靠濃度差作為推動(dòng)力,因此氣云的高濃度區(qū)域一般位于泄漏點(diǎn)的下風(fēng)側(cè)。將通風(fēng)流場(chǎng)特性與5個(gè)典型泄漏點(diǎn)的天然氣擴(kuò)散趨勢(shì)相結(jié)合,在距離艙室頂部1 m范圍內(nèi)布置可燃?xì)怏w探測(cè)器,如圖15所示。

    1)燃?xì)夤芫€區(qū)域1。天然氣泄漏后的氣云主體流向+x方向且擴(kuò)散過(guò)程相對(duì)緩慢,故在區(qū)域1中部和右側(cè)端部分別布置1號(hào)和2號(hào)探測(cè)器。

    2)燃?xì)夤芫€區(qū)域2。天然氣泄漏后的氣云主體呈明顯的-x方向流動(dòng)趨勢(shì),并由艙室頂部排氣風(fēng)機(jī)排出,故在區(qū)域2的風(fēng)機(jī)附近布置3號(hào)探測(cè)器,而在傳統(tǒng)布置方案中該區(qū)域未布置探測(cè)器。該區(qū)域處于強(qiáng)通風(fēng)的下風(fēng)側(cè)且靠近頂部排氣風(fēng)機(jī),天然氣泄漏后一般處于氣體的飄散流域內(nèi)。

    3)燃?xì)夤芫€區(qū)域3。艙室左側(cè)上方存在相對(duì)獨(dú)立的弱渦旋流動(dòng),此處發(fā)生泄漏后氣云會(huì)富集于泄漏點(diǎn)上方并向四周擴(kuò)散,右側(cè)發(fā)生泄漏后的氣云主體則會(huì)流向艙室中部,故在區(qū)域3左側(cè)布置探測(cè)器4,艙室中部布置探測(cè)器5,而在傳統(tǒng)布置方案中該區(qū)域未布置探測(cè)器。

    圖15 可燃?xì)怏w探測(cè)器布置圖Fig.15 Arrangement of flammable gas detectors

    2.4 排氣風(fēng)機(jī)防爆需求評(píng)估

    隨著設(shè)備防爆等級(jí)的提高,其采購(gòu)和維護(hù)費(fèi)用往往呈倍數(shù)增長(zhǎng),泄漏擴(kuò)散的模擬結(jié)果可以為設(shè)備防爆等級(jí)和危險(xiǎn)區(qū)域劃分提供定性定量的評(píng)估依據(jù),從而消除潛在的隱患或避免不必要的資源浪費(fèi)。

    以本文發(fā)電機(jī)室頂部的排氣風(fēng)機(jī)為例(圖5(a)),由泄漏擴(kuò)散的模擬結(jié)果可知,即使在危險(xiǎn)泄漏工況下,風(fēng)機(jī)附近累積的天然氣濃度也始終遠(yuǎn)小于天然氣爆炸下限(摩爾濃度4.9%,天然氣爆炸上限為摩爾濃度9.5%)。按照優(yōu)化后的可燃?xì)怏w探測(cè)器布置方案,即使存在長(zhǎng)時(shí)間的持續(xù)泄漏累積,最早接近其爆炸下限的區(qū)域也位于探測(cè)器附近,會(huì)及時(shí)觸發(fā)報(bào)警并切斷燃?xì)?。排氣風(fēng)機(jī)附近的天然氣濃度始終低于天然氣爆炸下限,設(shè)備即使發(fā)生故障也不會(huì)形成燃燒或引發(fā)爆炸,因此該發(fā)電機(jī)室無(wú)需選用防爆型排氣風(fēng)機(jī)。

    3 結(jié) 語(yǔ)

    本文建立了發(fā)電機(jī)室三維計(jì)算模型,首先分析了該艙室內(nèi)的非穩(wěn)態(tài)通風(fēng)流場(chǎng)特性,進(jìn)而對(duì)不同區(qū)域燃?xì)夤?yīng)管線的泄漏擴(kuò)散過(guò)程進(jìn)行模擬,得到不同區(qū)域內(nèi)發(fā)生泄漏后的天然氣擴(kuò)散趨勢(shì)和高濃度氣云所在位置,提出了可燃?xì)怏w探測(cè)器的優(yōu)化布置方案,為排氣風(fēng)機(jī)的防爆需求評(píng)估提供了定量依據(jù),為泄漏事故的預(yù)防和改進(jìn)設(shè)計(jì)提供了參考。后期可進(jìn)一步分析通風(fēng)流場(chǎng)特性,用以優(yōu)化風(fēng)管布置和風(fēng)量分布,實(shí)現(xiàn)降耗增效;進(jìn)一步計(jì)算燃爆風(fēng)險(xiǎn)及燃爆載荷,用以支撐結(jié)構(gòu)的局部加強(qiáng)設(shè)計(jì)。

    [1]WEI L,GENG P.A review on natural gas/diesel dual fuel combustion,emissions and performance[J].Fuel Processing Technology,2016,142:264-278.

    [2]KUMAR S,KWON H T,CHOI K H,et al.LNG:an eco-friendly cryogenic fuel for sustainable development[J].Applied Energy,2011,88(12):4264-4273.

    [3]LIN W,ZHANG N,GU A.LNG(liquefied natural gas):a necessary part in China's future energy infrastructure[J].Energy,2010,35(11):4383-4391.

    [4]THOMSON H,CORBETT J J,WINEBRAKE J J.Natural gas as a marine fuel[J].Energy Policy,2015,87:153-167.

    [5]孫培廷.“十二五”期間航運(yùn)業(yè)和造船業(yè)熱點(diǎn)與難點(diǎn)問(wèn)題[J].世界海運(yùn),2013,36(1):1-8.

    [6]張明鋒,鄧凱.簡(jiǎn)述船用LNG發(fā)展現(xiàn)狀[J].能源研究與管理,2013(3):16-20.ZHANG M F,DENG K.Introduction of LNG fuel ship's present situation[J].Energy Research and Management,2013(3):16-20(in Chinese).

    [7]KUMAR S,KWON H T,CHOI K H,et al.Current status and future projections of LNG demand and supplies:a global prospective[J].Energy Policy,2011,39(7):4097-4104.

    [8]趙秀雯,于力,柴建設(shè).天然氣管道安全[M].北京:化學(xué)工業(yè)出版社,2013.

    [9]方東曉.影響天然氣管道安全的危害因素分析及對(duì)策探討[J].上海煤氣,2008(3):38-41.FANG D X.Discussion on the analysis and its countermeasure of the elements on NG pipeline's security harmful[J].Shanghai Gas,2008(3):38-41(in Chinese).

    [10]OHBA R,KOUCHI A,HARA T,et al.Validation of heavy and light gas dispersion models for the safety analysis of LNG tank[J].Journal of Loss Prevention in the Process Industries,2004,17(5):325-337.

    [11]黃琴,蔣軍成.重氣泄漏擴(kuò)散實(shí)驗(yàn)的計(jì)算流體力學(xué)(CFD)模擬驗(yàn)證[J].中國(guó)安全科學(xué)學(xué)報(bào),2008,18(1):50-55.HUANG Q,JIANG J C.Simulation and verification of CFD on the dispersion of heavy gas leakage[J].China Safety Science Journal,2008,18(1):50-55(in Chinese).

    [12]DNV G L.Dual fuel concept analysis of fires and explosions in engine room[R].Norway:DNV G L,2006.

    [13]程浩力,劉德俊,劉倩倩,等.城燃管道街道峽谷泄漏擴(kuò)散CFD數(shù)值模擬[J].石油化工高等學(xué)校學(xué)報(bào),2011,24(4):60-63.CHENG H L,LIU D J,LIU Q Q,et al.CFD modeling of the releases and dispersion of accidental city gas pipeline in urban street canyons[J].Journal of Petrochemical Universities,2011,24(4):60-63(in Chinese).

    [14]王麗華.基于CFD的海洋平臺(tái)油氣泄漏擴(kuò)散仿真研究[D].青島:中國(guó)海洋大學(xué),2014.WANG L H.Study on the diffusion simulation of leaked oilgas on offshore platforms based on CFD[D].Qingdao:Ocean University of China,2014(in Chinese).

    [15]FU S S,YAN X P,ZHANG D,et al.Framework for the quantitative assessment of the risk of leakage from LNG-fueled vessels by anevent tree-CFD[J].Journal of Loss Prevention in the Process Industries,2016,43:42-52.

    [16]FIATES J,SANTOS R R C,NETO F F,et al.An alternative CFD tool for gas dispersion modelling of heavy gas[J].Journal of Loss Prevention in the Process Industries,2016,44:583-593.

    [17]IMO.Quantitative risk assessment of a fuel preparation room having a high-pressure fuel gas supply system on an LNG-fuelled ship[R].Korea:Sub-Committee on Carriage of Cargoes and Containers,2016.

    [18]VERSTEEG H K,MALALASEKERA W.An introduction to computational fluid dynamics:the finite volume method[M].[s.l.]Longman,1996.

    [19]Ansys Inc.Ansys fluent theory guide release 13.0[Z].Ansys Inc,2010:67-78.

    [20]CEHLIN M,MOSHFEGH B.Numerical modeling of a complex diffuser in a room with displacement ventilation[J].Building and Environment,2010,45(10):2240-2252.

    [21]郭棟鵬,姚仁太,喬清黨.建筑物對(duì)污染物擴(kuò)散影響的數(shù)值與風(fēng)洞模擬研究[J].實(shí)驗(yàn)流體力學(xué),2010,24(6):16-21.GUO D P,YAO R T,QIAO Q D.Numerical and wind-tunnel simulation on the effect of building on pollutant dispersion[J].Journal of Experiments in Fluid Mechanics,2010,24(6):16-21(in Chinese).

    [22]MAN Diesel and Turbo.Safety concept dual-fuel engines[R].Germany:MAN Diesel and Turbo,2012:12-20.

    [23]王新.天然氣管道泄漏擴(kuò)散事故危害評(píng)價(jià)[D].哈爾濱:哈爾濱工業(yè)大學(xué),2010.WANG X.Hazard assessment of natural gas pipeline leakage and diffusion[D].Harbin:Harbin Institute of Technology,2010(in Chinese).

    [24]李又綠,姚安林,李永杰.天然氣管道泄漏擴(kuò)散模型研究[J].天然氣工業(yè),2004,24(8):102-104.LI Y L,YAO A L,LI Y J.Study on diffusion model of gas pipeline leaking[J].Natural Gas Industry,2004,24(8):102-104(in Chinese).

    Leakage analysis of fuel gas pipe in large LNG carrier engine room

    CEN Zhuolun,LIU Tao,WANG Lei,DOU Xu
    R&D Dept.,Hudong-Zhonghua Shipbuilding(Group)Co.,Ltd.,Shanghai 200129,China

    U664.14

    A

    10.3969/j.issn.1673-3185.2017.05.014

    http://kns.cnki.net/kcms/detail/42.1755.TJ.20170926.1032.010.html期刊網(wǎng)址:www.ship-research.com

    岑卓倫,劉濤,王磊,等.大型LNG船發(fā)電機(jī)室的燃?xì)夤芫€泄漏分析[J].中國(guó)艦船研究,2017,12(5):109-119.

    CEN Z L,LIU T,WANG L,et al.Leakage analysis of fuel gas pipe in large LNG carrier engine room[J].Chinese Journal of Ship Research,2017,12(5):109-119.

    2017-03-14< class="emphasis_bold">網(wǎng)絡(luò)出版時(shí)間:

    時(shí)間:2017-9-26 10:32

    岑卓倫(通信作者),男,1988年生,碩士,工程師。研究方向:LNG低溫液貨處理系統(tǒng)開(kāi)發(fā)設(shè)計(jì)。E-mail:cenzhuolun@hotmail.com

    劉濤,男,1983年生,博士,高級(jí)工程師。研究方向:輪機(jī)系統(tǒng)開(kāi)發(fā)設(shè)計(jì)。E-mail:smuliutao@163.com

    猜你喜歡
    艙室燃?xì)?/a>發(fā)電機(jī)
    擰緊燃?xì)狻鞍踩y”
    船舶艙室通風(fēng)控制策略研究
    近期實(shí)施的燃?xì)鈬?guó)家標(biāo)準(zhǔn)
    煤氣與熱力(2022年4期)2022-05-23 12:45:00
    探討燃?xì)夤こ痰脑靸r(jià)控制及跟蹤審計(jì)
    大型發(fā)電機(jī)勵(lì)磁用旋轉(zhuǎn)變換器的開(kāi)發(fā)和應(yīng)用
    隨身攜帶的小發(fā)電機(jī)
    軍事文摘(2016年16期)2016-09-13 06:15:49
    柴油發(fā)電機(jī)負(fù)荷計(jì)算
    艦載導(dǎo)彈意外點(diǎn)火時(shí)艙室內(nèi)噴灑抑制的影響研究
    基于PCS-985B的發(fā)電機(jī)定子接地保護(hù)應(yīng)用及整定
    燃?xì)膺h(yuǎn)傳抄表系統(tǒng)分析
    河南科技(2014年7期)2014-02-27 14:11:17
    国产aⅴ精品一区二区三区波| 99久久国产精品久久久| 国产精品一区二区精品视频观看| 国产精品爽爽va在线观看网站 | 99久久99久久久精品蜜桃| 国产一区二区三区在线臀色熟女| 亚洲最大成人中文| 青草久久国产| 深夜精品福利| 黄网站色视频无遮挡免费观看| 视频在线观看一区二区三区| 少妇粗大呻吟视频| 50天的宝宝边吃奶边哭怎么回事| 中国美女看黄片| 中文字幕人成人乱码亚洲影| 久久精品国产综合久久久| 最新在线观看一区二区三区| 亚洲精品中文字幕一二三四区| 亚洲五月色婷婷综合| 亚洲国产精品久久男人天堂| 亚洲aⅴ乱码一区二区在线播放 | 最近最新中文字幕大全电影3 | 国产亚洲精品久久久久久毛片| 欧美大码av| 在线观看日韩欧美| а√天堂www在线а√下载| 亚洲久久久国产精品| 一级a爱片免费观看的视频| 18禁观看日本| 久久人人爽av亚洲精品天堂| 制服丝袜大香蕉在线| 久久欧美精品欧美久久欧美| 中文字幕高清在线视频| 在线观看免费午夜福利视频| 久久精品影院6| 日本黄色视频三级网站网址| 桃色一区二区三区在线观看| 免费女性裸体啪啪无遮挡网站| 女人高潮潮喷娇喘18禁视频| 国产99久久九九免费精品| 久久精品91无色码中文字幕| 国产三级在线视频| 在线av久久热| 嫩草影院精品99| 亚洲少妇的诱惑av| 脱女人内裤的视频| 精品熟女少妇八av免费久了| 亚洲中文av在线| 亚洲欧美精品综合一区二区三区| 国产麻豆69| 欧美日本中文国产一区发布| 亚洲成人免费电影在线观看| 中文字幕色久视频| 精品免费久久久久久久清纯| 两个人看的免费小视频| 美女高潮到喷水免费观看| 老熟妇仑乱视频hdxx| 级片在线观看| 国产亚洲av嫩草精品影院| 宅男免费午夜| 午夜福利高清视频| 韩国av一区二区三区四区| 亚洲精品国产色婷婷电影| 午夜精品久久久久久毛片777| 成人亚洲精品一区在线观看| 国产男靠女视频免费网站| 亚洲 欧美一区二区三区| 日韩一卡2卡3卡4卡2021年| 欧美激情极品国产一区二区三区| 亚洲人成网站在线播放欧美日韩| 免费一级毛片在线播放高清视频 | 欧美一区二区精品小视频在线| 亚洲精品在线美女| 国产精品国产高清国产av| 日韩欧美国产在线观看| 别揉我奶头~嗯~啊~动态视频| 美女午夜性视频免费| 老司机靠b影院| 好看av亚洲va欧美ⅴa在| 日本三级黄在线观看| 在线观看免费午夜福利视频| 嫩草影院精品99| 美女免费视频网站| 少妇 在线观看| 91九色精品人成在线观看| 又大又爽又粗| 亚洲,欧美精品.| 亚洲va日本ⅴa欧美va伊人久久| 一边摸一边抽搐一进一小说| 久热这里只有精品99| 久久精品91蜜桃| 9191精品国产免费久久| 成人国语在线视频| 国产精品二区激情视频| 午夜精品国产一区二区电影| 国产成+人综合+亚洲专区| 正在播放国产对白刺激| 男人的好看免费观看在线视频 | 夜夜躁狠狠躁天天躁| 51午夜福利影视在线观看| 搡老岳熟女国产| 一级黄色大片毛片| 精品一区二区三区av网在线观看| 无人区码免费观看不卡| 久久久久久久久久久久大奶| 中文字幕人妻丝袜一区二区| 精品人妻1区二区| 成人永久免费在线观看视频| 成人18禁在线播放| 日本一区二区免费在线视频| 国产亚洲精品久久久久5区| 亚洲久久久国产精品| 人人妻,人人澡人人爽秒播| 午夜激情av网站| 一级毛片女人18水好多| www.999成人在线观看| 午夜成年电影在线免费观看| 亚洲精品中文字幕在线视频| 久久婷婷成人综合色麻豆| 伦理电影免费视频| 欧美成狂野欧美在线观看| 日本vs欧美在线观看视频| 怎么达到女性高潮| 亚洲一区高清亚洲精品| 美女午夜性视频免费| 好男人电影高清在线观看| 侵犯人妻中文字幕一二三四区| 亚洲一区中文字幕在线| 在线十欧美十亚洲十日本专区| 99久久99久久久精品蜜桃| 久久久久久久久中文| 91九色精品人成在线观看| 99国产精品99久久久久| 欧美一级a爱片免费观看看 | 97碰自拍视频| 亚洲 国产 在线| 男人舔女人的私密视频| 亚洲一区二区三区色噜噜| 国产亚洲av高清不卡| av免费在线观看网站| 国产精品亚洲一级av第二区| 成年版毛片免费区| 午夜激情av网站| 国产乱人伦免费视频| 侵犯人妻中文字幕一二三四区| 国产精品一区二区精品视频观看| 午夜激情av网站| 大型黄色视频在线免费观看| 老司机靠b影院| 久久久久久免费高清国产稀缺| 又黄又粗又硬又大视频| 久久国产精品人妻蜜桃| 久久午夜综合久久蜜桃| 伦理电影免费视频| 亚洲在线自拍视频| 看免费av毛片| xxx96com| 成人18禁高潮啪啪吃奶动态图| 久久影院123| 久久久久久久久免费视频了| av天堂久久9| 国内久久婷婷六月综合欲色啪| 国产人伦9x9x在线观看| 法律面前人人平等表现在哪些方面| 成人亚洲精品av一区二区| 午夜福利成人在线免费观看| www.熟女人妻精品国产| 国产片内射在线| 国产日韩一区二区三区精品不卡| 国产日韩一区二区三区精品不卡| 黄片播放在线免费| 日韩大尺度精品在线看网址 | 亚洲欧美日韩无卡精品| 色在线成人网| 午夜福利一区二区在线看| 久久婷婷人人爽人人干人人爱 | 69精品国产乱码久久久| 国内精品久久久久久久电影| 大码成人一级视频| 欧洲精品卡2卡3卡4卡5卡区| 国产精品一区二区免费欧美| 在线国产一区二区在线| 国产成+人综合+亚洲专区| 一边摸一边抽搐一进一出视频| 成人特级黄色片久久久久久久| 午夜精品国产一区二区电影| 精品久久久久久,| 少妇被粗大的猛进出69影院| 啦啦啦 在线观看视频| 亚洲无线在线观看| 18禁美女被吸乳视频| 18美女黄网站色大片免费观看| 欧美乱妇无乱码| 久久久精品国产亚洲av高清涩受| 欧美成人免费av一区二区三区| 制服人妻中文乱码| 国产精品电影一区二区三区| 香蕉久久夜色| 国产精品乱码一区二三区的特点 | 黄色视频,在线免费观看| 美女高潮到喷水免费观看| 免费观看人在逋| 亚洲aⅴ乱码一区二区在线播放 | 久久精品亚洲精品国产色婷小说| 超碰成人久久| 在线观看免费视频日本深夜| 日本免费一区二区三区高清不卡 | 亚洲国产精品成人综合色| 国产亚洲av嫩草精品影院| 桃色一区二区三区在线观看| 中文字幕最新亚洲高清| 国产1区2区3区精品| 曰老女人黄片| 老熟妇仑乱视频hdxx| 巨乳人妻的诱惑在线观看| 国产精品,欧美在线| 国产精品自产拍在线观看55亚洲| 99久久国产精品久久久| 老司机福利观看| 久久青草综合色| 不卡av一区二区三区| 久久久国产成人精品二区| 最新在线观看一区二区三区| www.熟女人妻精品国产| 麻豆一二三区av精品| 精品福利观看| 亚洲全国av大片| 涩涩av久久男人的天堂| 欧美午夜高清在线| 国产午夜精品久久久久久| 欧美在线一区亚洲| 麻豆国产av国片精品| 精品欧美一区二区三区在线| 啦啦啦韩国在线观看视频| 757午夜福利合集在线观看| 国产99白浆流出| 韩国精品一区二区三区| 国产主播在线观看一区二区| 亚洲精品美女久久久久99蜜臀| 国产av精品麻豆| 久久久久九九精品影院| 国产单亲对白刺激| 精品一区二区三区av网在线观看| 91在线观看av| 国产亚洲av高清不卡| 久久婷婷成人综合色麻豆| 色在线成人网| 国产精品免费一区二区三区在线| 亚洲精品粉嫩美女一区| 高清黄色对白视频在线免费看| 一边摸一边抽搐一进一小说| 精品欧美国产一区二区三| 18禁黄网站禁片午夜丰满| 中出人妻视频一区二区| 中文字幕久久专区| 亚洲一区二区三区不卡视频| 精品国产亚洲在线| 国产成人系列免费观看| 长腿黑丝高跟| 老鸭窝网址在线观看| 女警被强在线播放| www.精华液| 国产亚洲精品综合一区在线观看 | 88av欧美| 正在播放国产对白刺激| √禁漫天堂资源中文www| 日日摸夜夜添夜夜添小说| 精品福利观看| 女人高潮潮喷娇喘18禁视频| 美女高潮喷水抽搐中文字幕| 少妇裸体淫交视频免费看高清 | 午夜福利成人在线免费观看| 一级毛片女人18水好多| 亚洲欧美日韩另类电影网站| www.自偷自拍.com| 在线免费观看的www视频| cao死你这个sao货| 人人妻人人澡人人看| 亚洲第一电影网av| 亚洲国产日韩欧美精品在线观看 | 少妇裸体淫交视频免费看高清 | 日本在线视频免费播放| 国产伦一二天堂av在线观看| 欧美亚洲日本最大视频资源| 一区福利在线观看| 亚洲性夜色夜夜综合| 99在线视频只有这里精品首页| 99久久精品国产亚洲精品| 色尼玛亚洲综合影院| 看免费av毛片| 手机成人av网站| 大型黄色视频在线免费观看| 国产精品久久久久久精品电影 | 午夜福利欧美成人| 欧美激情 高清一区二区三区| 性色av乱码一区二区三区2| tocl精华| 午夜久久久在线观看| 久久国产亚洲av麻豆专区| 亚洲欧美激情在线| 国产一区在线观看成人免费| 亚洲av电影在线进入| av网站免费在线观看视频| 国产97色在线日韩免费| 久久青草综合色| 国产精品免费视频内射| 美女免费视频网站| 天天添夜夜摸| 别揉我奶头~嗯~啊~动态视频| 久久久久久亚洲精品国产蜜桃av| 我的亚洲天堂| 亚洲欧美激情在线| 婷婷精品国产亚洲av在线| 在线天堂中文资源库| 久久天躁狠狠躁夜夜2o2o| 无人区码免费观看不卡| 国产av一区在线观看免费| 岛国在线观看网站| 日韩欧美一区二区三区在线观看| 一夜夜www| 三级毛片av免费| 精品久久久久久久久久免费视频| 色播在线永久视频| 两性午夜刺激爽爽歪歪视频在线观看 | 黑人操中国人逼视频| 精品国产超薄肉色丝袜足j| 国产成人系列免费观看| 国产一区二区三区在线臀色熟女| 精品熟女少妇八av免费久了| 午夜免费鲁丝| 精品国产亚洲在线| 欧美黑人欧美精品刺激| 国产精华一区二区三区| 正在播放国产对白刺激| 亚洲午夜精品一区,二区,三区| 亚洲国产精品999在线| 天天添夜夜摸| 国产精品秋霞免费鲁丝片| 啦啦啦韩国在线观看视频| 亚洲国产欧美日韩在线播放| 亚洲熟妇中文字幕五十中出| 成年人黄色毛片网站| 久久人人爽av亚洲精品天堂| 久久香蕉精品热| 色综合亚洲欧美另类图片| 亚洲国产精品成人综合色| 亚洲精品国产一区二区精华液| 色老头精品视频在线观看| av欧美777| a级毛片在线看网站| 亚洲国产中文字幕在线视频| 香蕉久久夜色| 亚洲自拍偷在线| 大陆偷拍与自拍| aaaaa片日本免费| 中文字幕av电影在线播放| 欧美一区二区精品小视频在线| 亚洲成av片中文字幕在线观看| 大香蕉久久成人网| 精品第一国产精品| 国产一区二区三区视频了| 美女高潮喷水抽搐中文字幕| 男人的好看免费观看在线视频 | 亚洲国产看品久久| 99国产精品一区二区蜜桃av| 美女扒开内裤让男人捅视频| 免费无遮挡裸体视频| 精品一品国产午夜福利视频| 成人av一区二区三区在线看| 欧美激情 高清一区二区三区| 在线观看免费日韩欧美大片| 中文字幕久久专区| 此物有八面人人有两片| 欧美日韩亚洲综合一区二区三区_| 精品欧美一区二区三区在线| 欧美不卡视频在线免费观看 | e午夜精品久久久久久久| 久久精品亚洲精品国产色婷小说| 亚洲欧美日韩无卡精品| 精品免费久久久久久久清纯| www.精华液| 动漫黄色视频在线观看| 亚洲av成人av| 老司机午夜十八禁免费视频| 色综合婷婷激情| 成人三级黄色视频| cao死你这个sao货| 制服人妻中文乱码| 日本黄色视频三级网站网址| 黄色视频不卡| 亚洲成a人片在线一区二区| 欧美 亚洲 国产 日韩一| 中文字幕人成人乱码亚洲影| 久久天躁狠狠躁夜夜2o2o| 欧美大码av| 午夜视频精品福利| 亚洲中文av在线| 成年女人毛片免费观看观看9| 免费无遮挡裸体视频| 99热只有精品国产| 一级作爱视频免费观看| 伦理电影免费视频| 亚洲精品中文字幕在线视频| 天堂动漫精品| 欧美人与性动交α欧美精品济南到| 怎么达到女性高潮| 久久久国产欧美日韩av| 狠狠狠狠99中文字幕| 天堂影院成人在线观看| 久久久久久人人人人人| 亚洲性夜色夜夜综合| 久久草成人影院| 咕卡用的链子| 国产蜜桃级精品一区二区三区| 视频在线观看一区二区三区| 国产人伦9x9x在线观看| 搞女人的毛片| 熟妇人妻久久中文字幕3abv| 大码成人一级视频| 国产色视频综合| 久久久久久国产a免费观看| 午夜福利欧美成人| 国产精品久久久久久人妻精品电影| 在线播放国产精品三级| 国产成人欧美| 国产成年人精品一区二区| bbb黄色大片| 国产成+人综合+亚洲专区| 精品人妻在线不人妻| 99在线视频只有这里精品首页| 亚洲一区二区三区不卡视频| 亚洲视频免费观看视频| 免费在线观看黄色视频的| 久久精品亚洲精品国产色婷小说| 乱人伦中国视频| 免费搜索国产男女视频| 国产片内射在线| 亚洲全国av大片| 亚洲黑人精品在线| 亚洲欧洲精品一区二区精品久久久| 亚洲国产精品999在线| www.999成人在线观看| 操美女的视频在线观看| 免费无遮挡裸体视频| 亚洲精品粉嫩美女一区| 久久久久国产一级毛片高清牌| 国产精品1区2区在线观看.| 亚洲中文av在线| 成人亚洲精品av一区二区| 97人妻天天添夜夜摸| 在线天堂中文资源库| 女人被狂操c到高潮| 久久国产精品影院| 露出奶头的视频| √禁漫天堂资源中文www| 国产欧美日韩一区二区三区在线| 欧美激情高清一区二区三区| 国产成人啪精品午夜网站| 18禁观看日本| 亚洲黑人精品在线| 真人做人爱边吃奶动态| 久热这里只有精品99| 操出白浆在线播放| 国产精品1区2区在线观看.| 19禁男女啪啪无遮挡网站| 免费久久久久久久精品成人欧美视频| 麻豆久久精品国产亚洲av| 欧美激情高清一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 成人三级黄色视频| 久9热在线精品视频| 又黄又爽又免费观看的视频| 国产一区二区激情短视频| 多毛熟女@视频| 国产精品久久久久久亚洲av鲁大| 成人免费观看视频高清| e午夜精品久久久久久久| 亚洲五月色婷婷综合| 久久中文字幕人妻熟女| 精品久久久久久成人av| 99国产精品一区二区三区| 亚洲精品久久成人aⅴ小说| 操出白浆在线播放| 一级毛片精品| 亚洲精品在线观看二区| 黄片大片在线免费观看| 久久精品亚洲熟妇少妇任你| 亚洲中文av在线| 伦理电影免费视频| 女人高潮潮喷娇喘18禁视频| 久久久久国内视频| 搞女人的毛片| 国产精品日韩av在线免费观看 | 日本免费一区二区三区高清不卡 | 欧美绝顶高潮抽搐喷水| 天堂动漫精品| 多毛熟女@视频| 色在线成人网| 久99久视频精品免费| 一级a爱视频在线免费观看| 丝袜美腿诱惑在线| 天堂影院成人在线观看| 国产xxxxx性猛交| 日韩一卡2卡3卡4卡2021年| 精品国产超薄肉色丝袜足j| 色播亚洲综合网| 午夜福利影视在线免费观看| 国产午夜精品久久久久久| 亚洲天堂国产精品一区在线| av天堂在线播放| 欧美乱妇无乱码| 亚洲精品一卡2卡三卡4卡5卡| 天天添夜夜摸| 日日夜夜操网爽| 免费av毛片视频| 久久草成人影院| 亚洲伊人色综图| 亚洲人成电影免费在线| 动漫黄色视频在线观看| 免费人成视频x8x8入口观看| netflix在线观看网站| 久久久水蜜桃国产精品网| 91大片在线观看| 中文字幕人成人乱码亚洲影| 极品人妻少妇av视频| 亚洲熟女毛片儿| 免费搜索国产男女视频| 在线观看舔阴道视频| av有码第一页| 人人澡人人妻人| av视频免费观看在线观看| 精品久久蜜臀av无| 国产精品久久久久久亚洲av鲁大| 欧美色视频一区免费| 亚洲五月色婷婷综合| av视频在线观看入口| 91在线观看av| 亚洲av日韩精品久久久久久密| 9色porny在线观看| av视频免费观看在线观看| 成人亚洲精品av一区二区| 亚洲自拍偷在线| 亚洲国产欧美一区二区综合| 亚洲最大成人中文| 亚洲九九香蕉| 色播在线永久视频| 久久久久久久午夜电影| 啦啦啦观看免费观看视频高清 | 中出人妻视频一区二区| 亚洲伊人色综图| bbb黄色大片| 午夜精品国产一区二区电影| 视频在线观看一区二区三区| 成年版毛片免费区| 国产高清有码在线观看视频 | 麻豆久久精品国产亚洲av| 国产成人影院久久av| 亚洲专区国产一区二区| 亚洲黑人精品在线| 日韩精品中文字幕看吧| 丝袜美腿诱惑在线| 亚洲中文日韩欧美视频| 亚洲熟妇熟女久久| 国内精品久久久久精免费| 国产成年人精品一区二区| 久久国产精品人妻蜜桃| 国产私拍福利视频在线观看| 国产高清videossex| 日韩欧美一区二区三区在线观看| √禁漫天堂资源中文www| 男女下面进入的视频免费午夜 | 男女之事视频高清在线观看| 欧美绝顶高潮抽搐喷水| 69av精品久久久久久| 老司机深夜福利视频在线观看| 亚洲九九香蕉| 香蕉久久夜色| 一夜夜www| 久久 成人 亚洲| 亚洲精品美女久久av网站| 50天的宝宝边吃奶边哭怎么回事| 在线免费观看的www视频| 日本在线视频免费播放| 欧美黑人精品巨大| 成人18禁在线播放| 97碰自拍视频| 国产极品粉嫩免费观看在线| 天天躁夜夜躁狠狠躁躁| 国产欧美日韩综合在线一区二区| 欧美一区二区精品小视频在线| 午夜精品在线福利| 国产一区二区三区视频了| 制服诱惑二区| 欧美亚洲日本最大视频资源| 一区在线观看完整版| 久久久久久人人人人人| 99香蕉大伊视频| 亚洲专区中文字幕在线| 国产成人一区二区三区免费视频网站| 久久欧美精品欧美久久欧美| 日韩一卡2卡3卡4卡2021年| 女人高潮潮喷娇喘18禁视频| 久久久国产成人精品二区| 美女免费视频网站| 91麻豆av在线| av在线播放免费不卡| www.精华液| 久久天堂一区二区三区四区| 九色国产91popny在线| 88av欧美| 国产精品电影一区二区三区| 久久久久久久午夜电影| 香蕉久久夜色| 国产高清有码在线观看视频 | 国产真人三级小视频在线观看|