• <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
    亚洲人成77777在线视频| 国产精品 欧美亚洲| 国产成人av激情在线播放| 欧美国产精品va在线观看不卡| 嫩草影视91久久| 亚洲精品美女久久久久99蜜臀| 国产精品久久久久久人妻精品电影| 久99久视频精品免费| 国产精品永久免费网站| 侵犯人妻中文字幕一二三四区| tube8黄色片| 日韩视频一区二区在线观看| 日本wwww免费看| xxxhd国产人妻xxx| 美女国产高潮福利片在线看| 80岁老熟妇乱子伦牲交| 亚洲精品自拍成人| 最新美女视频免费是黄的| 欧美日韩精品网址| 一区福利在线观看| 性色av乱码一区二区三区2| 精品无人区乱码1区二区| 国产激情欧美一区二区| 乱人伦中国视频| 国产欧美亚洲国产| 在线观看免费日韩欧美大片| 在线播放国产精品三级| 天堂俺去俺来也www色官网| 丰满的人妻完整版| 国产精品偷伦视频观看了| 在线播放国产精品三级| 亚洲精品一卡2卡三卡4卡5卡| 免费少妇av软件| 少妇裸体淫交视频免费看高清 | 两人在一起打扑克的视频| 国产成人免费无遮挡视频| 亚洲精品成人av观看孕妇| 99re在线观看精品视频| 久久亚洲真实| 搡老乐熟女国产| 久久国产精品影院| 99精品欧美一区二区三区四区| 欧美日韩成人在线一区二区| 欧美最黄视频在线播放免费 | 91字幕亚洲| 精品人妻熟女毛片av久久网站| 超色免费av| 中文字幕最新亚洲高清| 国产亚洲精品久久久久久毛片 | 91精品国产国语对白视频| 美女视频免费永久观看网站| videosex国产| 99精品久久久久人妻精品| 久久久久视频综合| 亚洲免费av在线视频| 久久久久久免费高清国产稀缺| 欧美日韩瑟瑟在线播放| 欧美激情久久久久久爽电影 | 麻豆国产av国片精品| 久久久久久亚洲精品国产蜜桃av| √禁漫天堂资源中文www| 成人精品一区二区免费| 咕卡用的链子| 在线观看免费视频日本深夜| 精品一区二区三区视频在线观看免费 | 制服诱惑二区| 久久久精品国产亚洲av高清涩受| 中文字幕最新亚洲高清| 亚洲av成人一区二区三| av免费在线观看网站| 国产精品综合久久久久久久免费 | 十八禁人妻一区二区| 成人手机av| 女警被强在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲成国产人片在线观看| 色综合欧美亚洲国产小说| 丰满的人妻完整版| 天堂√8在线中文| av有码第一页| 国产精品一区二区在线观看99| 69av精品久久久久久| 亚洲欧美一区二区三区黑人| 大香蕉久久网| 久久久久久久精品吃奶| 国产精品成人在线| 久久午夜亚洲精品久久| 女警被强在线播放| 法律面前人人平等表现在哪些方面| √禁漫天堂资源中文www| 黄片播放在线免费| 色播在线永久视频| 777久久人妻少妇嫩草av网站| 亚洲精品在线观看二区| 国产淫语在线视频| 19禁男女啪啪无遮挡网站| 色精品久久人妻99蜜桃| bbb黄色大片| 欧美国产精品一级二级三级| 美国免费a级毛片| 亚洲色图av天堂| 色在线成人网| 人成视频在线观看免费观看| 国产在线精品亚洲第一网站| av中文乱码字幕在线| 国产亚洲精品久久久久5区| 亚洲成a人片在线一区二区| 一本大道久久a久久精品| 国产xxxxx性猛交| a级毛片在线看网站| 黑人欧美特级aaaaaa片| 亚洲美女黄片视频| 国产精华一区二区三区| 亚洲成人手机| 久久久精品免费免费高清| 这个男人来自地球电影免费观看| 国产精品电影一区二区三区 | 一区在线观看完整版| 国产1区2区3区精品| 欧美+亚洲+日韩+国产| 欧美日韩亚洲综合一区二区三区_| 成人国语在线视频| 国产亚洲av高清不卡| 天堂动漫精品| 少妇的丰满在线观看| 国产精品国产高清国产av | 精品少妇久久久久久888优播| 精品电影一区二区在线| 少妇裸体淫交视频免费看高清 | 在线播放国产精品三级| 男人舔女人的私密视频| 黄色视频,在线免费观看| 国产成人一区二区三区免费视频网站| 午夜福利影视在线免费观看| 国产xxxxx性猛交| 亚洲人成电影观看| 好看av亚洲va欧美ⅴa在| 日韩成人在线观看一区二区三区| 亚洲精华国产精华精| 9色porny在线观看| 欧美av亚洲av综合av国产av| 后天国语完整版免费观看| 午夜日韩欧美国产| 亚洲精品久久成人aⅴ小说| 久久久久国内视频| 精品国产国语对白av| 99re在线观看精品视频| 久99久视频精品免费| 中文字幕精品免费在线观看视频| 日韩精品免费视频一区二区三区| 丝袜人妻中文字幕| 咕卡用的链子| 欧美黑人欧美精品刺激| 老司机午夜福利在线观看视频| 麻豆成人av在线观看| 亚洲美女黄片视频| 一区在线观看完整版| 亚洲一区二区三区不卡视频| 一本一本久久a久久精品综合妖精| 国产国语露脸激情在线看| 欧美 亚洲 国产 日韩一| 一区福利在线观看| 久久这里只有精品19| 超碰成人久久| 国产欧美日韩精品亚洲av| 很黄的视频免费| 9热在线视频观看99| 免费在线观看完整版高清| 午夜精品在线福利| 99国产综合亚洲精品| 亚洲精品一二三| 99国产极品粉嫩在线观看| 99国产精品免费福利视频| 黄色视频不卡| 1024视频免费在线观看| 国产成人啪精品午夜网站| 免费在线观看视频国产中文字幕亚洲| 国产精品久久久久久精品古装| 午夜福利欧美成人| 大陆偷拍与自拍| 天天操日日干夜夜撸| 在线av久久热| 99久久99久久久精品蜜桃| 亚洲 国产 在线| 另类亚洲欧美激情| 亚洲aⅴ乱码一区二区在线播放 | 两个人免费观看高清视频| 精品欧美一区二区三区在线| 中出人妻视频一区二区| 亚洲av成人一区二区三| 丝袜在线中文字幕| 欧美午夜高清在线| 男女下面插进去视频免费观看| 两人在一起打扑克的视频| 亚洲国产欧美网| 一区二区三区激情视频| 精品国产一区二区三区四区第35| 十八禁网站免费在线| 制服诱惑二区| 看免费av毛片| 动漫黄色视频在线观看| 1024香蕉在线观看| 18禁裸乳无遮挡免费网站照片 | 50天的宝宝边吃奶边哭怎么回事| 身体一侧抽搐| 亚洲国产看品久久| av欧美777| 亚洲一区中文字幕在线| 精品亚洲成国产av| 另类亚洲欧美激情| 丝袜美腿诱惑在线| 成人黄色视频免费在线看| 欧美成人免费av一区二区三区 | 王馨瑶露胸无遮挡在线观看| 一边摸一边做爽爽视频免费| 亚洲片人在线观看| 亚洲av日韩在线播放| 久久久国产成人精品二区 | 在线十欧美十亚洲十日本专区| 国产野战对白在线观看| 男女之事视频高清在线观看| 国产成人一区二区三区免费视频网站| 99国产精品99久久久久| 不卡av一区二区三区| 欧美在线黄色| 丝袜美腿诱惑在线| 国产精品亚洲av一区麻豆| 国产精品久久久久成人av| 亚洲第一av免费看| 国产极品粉嫩免费观看在线| 日韩中文字幕欧美一区二区| 69av精品久久久久久| 美女福利国产在线| 男女床上黄色一级片免费看| 大片电影免费在线观看免费| 日韩成人在线观看一区二区三区| x7x7x7水蜜桃| 老司机靠b影院| 国产精品久久久久久人妻精品电影| 国产1区2区3区精品| 久久久精品免费免费高清| 亚洲欧美激情综合另类| 午夜视频精品福利| 丝袜人妻中文字幕| 欧美精品亚洲一区二区| 成人亚洲精品一区在线观看| 色尼玛亚洲综合影院| 韩国精品一区二区三区| 麻豆国产av国片精品| 欧美人与性动交α欧美软件| 亚洲av日韩精品久久久久久密| 中文字幕最新亚洲高清| 日韩欧美一区视频在线观看| 女人被狂操c到高潮| 欧美日韩黄片免| 亚洲黑人精品在线| 无限看片的www在线观看| 99在线人妻在线中文字幕 | 亚洲熟妇熟女久久| 精品第一国产精品| 国产欧美日韩一区二区三| 夜夜夜夜夜久久久久| 国产精品久久久久久精品古装| cao死你这个sao货| av不卡在线播放| 成年人免费黄色播放视频| www.熟女人妻精品国产| 国产主播在线观看一区二区| 大码成人一级视频| tube8黄色片| 免费在线观看黄色视频的| 精品福利观看| cao死你这个sao货| 欧美日韩瑟瑟在线播放| 日韩欧美免费精品| 99久久99久久久精品蜜桃| 天堂动漫精品| 久久人人爽av亚洲精品天堂| 久久精品熟女亚洲av麻豆精品| 两人在一起打扑克的视频| 成人av一区二区三区在线看| 制服诱惑二区| 999久久久精品免费观看国产| 久久久久久久精品吃奶| 一本大道久久a久久精品| 亚洲中文av在线| 成人18禁高潮啪啪吃奶动态图| 久久狼人影院| 精品欧美一区二区三区在线| tube8黄色片| 久久久精品国产亚洲av高清涩受| 日韩免费av在线播放| a在线观看视频网站| 美女国产高潮福利片在线看| 亚洲一码二码三码区别大吗| 老司机深夜福利视频在线观看| 久久精品人人爽人人爽视色| 精品高清国产在线一区| 天堂动漫精品| 国产成人av激情在线播放| 国产精品久久久人人做人人爽| 成年女人毛片免费观看观看9 | 亚洲久久久国产精品| 91在线观看av| 精品福利观看| 999精品在线视频| 午夜福利免费观看在线| 欧美黄色淫秽网站| 一边摸一边抽搐一进一出视频| 女人久久www免费人成看片| 国产精品欧美亚洲77777| 亚洲五月色婷婷综合| 国产欧美日韩精品亚洲av| 99热国产这里只有精品6| 国产在视频线精品| 美女午夜性视频免费| 超色免费av| 人妻丰满熟妇av一区二区三区 | 亚洲专区中文字幕在线| 亚洲成a人片在线一区二区| www.精华液| 久热爱精品视频在线9| av视频免费观看在线观看| 操出白浆在线播放| 国产精品国产高清国产av | 日本一区二区免费在线视频| 老熟妇乱子伦视频在线观看| 捣出白浆h1v1| 国产乱人伦免费视频| 亚洲精品美女久久av网站| 国产亚洲欧美98| 亚洲色图av天堂| 亚洲精品美女久久av网站| 亚洲综合色网址| 黄网站色视频无遮挡免费观看| 日韩三级视频一区二区三区| 国产aⅴ精品一区二区三区波| 黄色女人牲交| 久久影院123| 一a级毛片在线观看| 91麻豆av在线| avwww免费| 一级片'在线观看视频| 国产精品久久久av美女十八| 成人三级做爰电影| 久热这里只有精品99| 国产精品久久电影中文字幕 | 亚洲欧美激情综合另类| 国产野战对白在线观看| 老司机午夜福利在线观看视频| 久久久国产欧美日韩av| 成人国语在线视频| 欧美在线黄色| 夫妻午夜视频| 80岁老熟妇乱子伦牲交| 国产日韩一区二区三区精品不卡| 免费在线观看完整版高清| 久久人妻av系列| 黄片小视频在线播放| www.熟女人妻精品国产| 欧美黄色片欧美黄色片| 国产99久久九九免费精品| 91精品三级在线观看| 精品少妇久久久久久888优播| 十八禁网站免费在线| 人妻丰满熟妇av一区二区三区 | 成人影院久久| 高清视频免费观看一区二区| 精品少妇久久久久久888优播| 十八禁网站免费在线| 亚洲欧美日韩高清在线视频| 亚洲久久久国产精品| 高潮久久久久久久久久久不卡| 亚洲熟妇中文字幕五十中出 | 亚洲全国av大片| 黄色片一级片一级黄色片| 精品电影一区二区在线| av片东京热男人的天堂| 国产一区有黄有色的免费视频| 精品久久蜜臀av无| 人妻 亚洲 视频| 国产精品免费大片| 亚洲成人国产一区在线观看| 免费观看精品视频网站| 一个人免费在线观看的高清视频| 日韩免费高清中文字幕av| 中文字幕另类日韩欧美亚洲嫩草| 国产精品国产高清国产av | 亚洲成人免费电影在线观看| 波多野结衣一区麻豆| 少妇粗大呻吟视频| 99久久国产精品久久久| 国产野战对白在线观看| 大片电影免费在线观看免费| 欧美午夜高清在线| 亚洲精品久久午夜乱码| 欧美日韩中文字幕国产精品一区二区三区 | 成人亚洲精品一区在线观看| 亚洲专区字幕在线| 一级片免费观看大全| 亚洲精品成人av观看孕妇| 亚洲五月婷婷丁香| 纯流量卡能插随身wifi吗| 久99久视频精品免费| av线在线观看网站| 精品午夜福利视频在线观看一区| 在线永久观看黄色视频| 精品国产乱子伦一区二区三区| 亚洲精品粉嫩美女一区| 亚洲av成人一区二区三| 在线看a的网站| 亚洲精品一卡2卡三卡4卡5卡| 欧美亚洲 丝袜 人妻 在线| 国产一区二区三区视频了| 在线观看免费午夜福利视频| 精品亚洲成a人片在线观看| 中文字幕色久视频| 丰满的人妻完整版| 欧美最黄视频在线播放免费 | 欧美精品人与动牲交sv欧美| 久久久精品区二区三区| 国产欧美日韩综合在线一区二区| 大型黄色视频在线免费观看| 可以免费在线观看a视频的电影网站| 99re在线观看精品视频| 成人18禁高潮啪啪吃奶动态图| 一二三四社区在线视频社区8| 日日摸夜夜添夜夜添小说| 精品久久久久久久毛片微露脸| 国产精品久久久久成人av| 亚洲黑人精品在线| 欧美日韩黄片免| bbb黄色大片| 亚洲成人免费av在线播放| 少妇的丰满在线观看| 国产主播在线观看一区二区| 50天的宝宝边吃奶边哭怎么回事| 免费高清在线观看日韩| 老熟妇乱子伦视频在线观看| 中文欧美无线码| 大型黄色视频在线免费观看| 欧美黄色淫秽网站| 丰满的人妻完整版| 日日摸夜夜添夜夜添小说| 亚洲九九香蕉| 国产精品免费视频内射| 别揉我奶头~嗯~啊~动态视频| 亚洲专区中文字幕在线| 777米奇影视久久| 国产欧美日韩一区二区三| 一级毛片女人18水好多| 黄色视频,在线免费观看| 亚洲av成人不卡在线观看播放网| 亚洲人成伊人成综合网2020| 国产淫语在线视频| 9色porny在线观看| 久久国产精品大桥未久av| 久久久久国产精品人妻aⅴ院 | 精品国产超薄肉色丝袜足j| 欧美日韩一级在线毛片| 无遮挡黄片免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产蜜桃级精品一区二区三区 | 国内毛片毛片毛片毛片毛片| 国产精品1区2区在线观看. | 午夜两性在线视频| 国产国语露脸激情在线看| 新久久久久国产一级毛片| 三级毛片av免费| 丝袜在线中文字幕| 中文字幕另类日韩欧美亚洲嫩草| 中国美女看黄片| 久久精品亚洲精品国产色婷小说| 一进一出好大好爽视频| 久久香蕉国产精品| 国产日韩一区二区三区精品不卡| 午夜福利欧美成人| 免费在线观看亚洲国产| 每晚都被弄得嗷嗷叫到高潮| 午夜精品在线福利| 亚洲欧美激情在线| 一边摸一边抽搐一进一小说 | 啦啦啦 在线观看视频| 怎么达到女性高潮| 男女床上黄色一级片免费看| 丰满迷人的少妇在线观看| 久久精品国产清高在天天线| 伦理电影免费视频| 日韩欧美一区视频在线观看| 国产亚洲欧美精品永久| 国产精品久久久av美女十八| 18禁裸乳无遮挡动漫免费视频| 亚洲少妇的诱惑av| 女人精品久久久久毛片| 国产在线一区二区三区精| 日韩有码中文字幕| 最新美女视频免费是黄的| 精品一区二区三区视频在线观看免费 | 王馨瑶露胸无遮挡在线观看| 成人亚洲精品一区在线观看| 在线观看免费高清a一片| 欧美日本中文国产一区发布| 人人妻人人添人人爽欧美一区卜| 九色亚洲精品在线播放| 男女下面插进去视频免费观看| 又黄又爽又免费观看的视频| 最近最新中文字幕大全免费视频| 99国产精品一区二区三区| а√天堂www在线а√下载 | 亚洲一码二码三码区别大吗| 757午夜福利合集在线观看| 亚洲av日韩在线播放| 又黄又爽又免费观看的视频| 欧美日韩亚洲综合一区二区三区_| 亚洲国产精品合色在线| 久久久国产成人免费| 搡老岳熟女国产| 黄色a级毛片大全视频| 国产有黄有色有爽视频| 丰满饥渴人妻一区二区三| 99久久精品国产亚洲精品| 国产欧美日韩一区二区三| 亚洲av片天天在线观看| 日韩人妻精品一区2区三区| 成人av一区二区三区在线看| 校园春色视频在线观看| 黄色视频不卡| 国产乱人伦免费视频| bbb黄色大片| 久热这里只有精品99| 国产1区2区3区精品| 久久亚洲真实| 乱人伦中国视频| 亚洲欧洲精品一区二区精品久久久| 深夜精品福利| 欧美人与性动交α欧美软件| 欧美色视频一区免费| 国产成人系列免费观看| 一级毛片女人18水好多| 亚洲aⅴ乱码一区二区在线播放 | 极品教师在线免费播放| 国产精品99久久99久久久不卡| 18禁黄网站禁片午夜丰满| 一区二区日韩欧美中文字幕| 亚洲中文日韩欧美视频| 免费看十八禁软件| av中文乱码字幕在线| 久久人妻av系列| 欧美日韩乱码在线| 老司机靠b影院| 亚洲精品国产一区二区精华液| 大香蕉久久网| 在线观看66精品国产| x7x7x7水蜜桃| 黄片小视频在线播放| 高清欧美精品videossex| 亚洲色图 男人天堂 中文字幕| 欧美在线黄色| 少妇猛男粗大的猛烈进出视频| 一级,二级,三级黄色视频| 人人妻,人人澡人人爽秒播| 最新美女视频免费是黄的| 国产伦人伦偷精品视频| 亚洲国产欧美日韩在线播放| 国产成人啪精品午夜网站| 精品久久久久久久毛片微露脸| 中文亚洲av片在线观看爽 | 亚洲精品一二三| 一本大道久久a久久精品| 久9热在线精品视频| 亚洲片人在线观看| 精品国产国语对白av| 国产精品成人在线| 国产不卡一卡二| 色综合婷婷激情| 男人的好看免费观看在线视频 | 久久久久久久久久久久大奶| 老熟女久久久| 精品久久久久久久久久免费视频 | 一区在线观看完整版| av欧美777| 亚洲欧美一区二区三区久久| 久久性视频一级片| av在线播放免费不卡| 亚洲av片天天在线观看| 久久九九热精品免费| av在线播放免费不卡| 99久久99久久久精品蜜桃| 中文欧美无线码| 一a级毛片在线观看| 亚洲午夜理论影院| 如日韩欧美国产精品一区二区三区| 亚洲国产看品久久| 黑人操中国人逼视频| 两性夫妻黄色片| 在线观看免费视频网站a站| 一级毛片高清免费大全| 大香蕉久久成人网| 国产精品香港三级国产av潘金莲| 精品国产国语对白av| 波多野结衣一区麻豆| 少妇被粗大的猛进出69影院| 成人永久免费在线观看视频| 午夜精品久久久久久毛片777| 80岁老熟妇乱子伦牲交| 一本大道久久a久久精品| 国产极品粉嫩免费观看在线| 亚洲,欧美精品.| 国产精品国产高清国产av | 国产精品国产av在线观看| 欧美乱色亚洲激情| 操出白浆在线播放| 精品乱码久久久久久99久播| 曰老女人黄片| 日本vs欧美在线观看视频|