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

    基于控制室爆炸載荷的可燃?xì)庠频葍r(jià)尺寸閾值

    2017-12-22 05:37:44齊心歌王海清宋賢生陳國(guó)明
    化工學(xué)報(bào) 2017年12期
    關(guān)鍵詞:火氣控制室等價(jià)

    齊心歌,王海清,宋賢生,陳國(guó)明

    (1中國(guó)石油大學(xué)(華東)安全科學(xué)與工程系,山東 青島 266580;2勞氏瑞安咨詢有限公司,上海 200001)

    基于控制室爆炸載荷的可燃?xì)庠频葍r(jià)尺寸閾值

    齊心歌1,王海清1,宋賢生2,陳國(guó)明1

    (1中國(guó)石油大學(xué)(華東)安全科學(xué)與工程系,山東 青島 266580;2勞氏瑞安咨詢有限公司,上海 200001)

    火氣系統(tǒng)是防控氣體泄漏災(zāi)害的重要安全屏障??扇?xì)怏w泄漏后與空氣混合形成可燃?xì)庠?,其爆炸屬于體積爆炸,具有復(fù)雜性和多變性,因此將其折算為等價(jià)氣云并提出閾值尺寸計(jì)算方法,是實(shí)現(xiàn)火氣系統(tǒng)探測(cè)器網(wǎng)絡(luò)量化布設(shè)的關(guān)鍵輸入指標(biāo)。選取控制室為受體,以載荷作為爆炸沖擊波超壓的臨界值,運(yùn)用多能法倒序計(jì)算,逆推得到對(duì)應(yīng)的等價(jià)氣云尺寸作為探測(cè)閾值,進(jìn)而利用等價(jià)氣云方法與高斯擴(kuò)散模型得到火氣系統(tǒng)探測(cè)臨界時(shí)間。通過某 LNG罐區(qū)案例分析,定量確定了該罐區(qū)可承載的氣云尺寸最大值及擴(kuò)散臨界時(shí)間。數(shù)值計(jì)算表明,等價(jià)氣云尺寸閾值不僅可以作為火氣系統(tǒng)探測(cè)設(shè)計(jì)的量化輸入指標(biāo)依據(jù),并可對(duì)探測(cè)時(shí)間設(shè)置以及氣體泄漏及爆炸的防控措施提供理論支持。

    氣云;爆炸;火氣系統(tǒng);爆炸載荷;等價(jià)氣云尺寸;閾值;安全;模擬

    引 言

    可燃?xì)怏w泄漏后,與空氣混合形成可燃?xì)庠?,遇點(diǎn)火源后易發(fā)生爆炸事故,對(duì)人員、設(shè)備、環(huán)境等造成嚴(yán)重危害。通過分析已經(jīng)發(fā)生的氣云爆炸事故,可以得到氣云爆炸的主要傷害形式為沖擊波超壓與熱輻射。由于沖擊波超壓的傷害半徑較大、損傷更嚴(yán)重,因此本文主要研究沖擊波超壓的影響[1-3]。

    氣云爆炸屬于體積爆炸,與凝聚相點(diǎn)源爆炸在燃燒速率、爆炸形式、點(diǎn)火概率、造成后果的嚴(yán)重程度等方面均有很大區(qū)別[4-5]。氣云爆炸受到多種因素影響,主要包括內(nèi)部因素如可燃?xì)怏w密度、反應(yīng)活性、可燃?xì)怏w與空氣混合的均勻程度、形成可燃?xì)庠频捏w積大小等;外部因素如點(diǎn)火時(shí)間、概率與地點(diǎn)、障礙物的大小與形狀、大氣環(huán)境的改變等[6]。利用火氣系統(tǒng)對(duì)氣體泄漏進(jìn)行探測(cè),達(dá)到對(duì)氣云爆炸預(yù)防與控制的目的;針對(duì)氣云擴(kuò)散情形復(fù)雜且多變的特點(diǎn),提出等價(jià)氣云的概念。等價(jià)氣云就是將一定空間范圍內(nèi),形狀不規(guī)則,密度、濃度不均勻的非均相氣云,用理想的均相單位立方體氣云代替,均相單位立方體組成的氣云即為等價(jià)氣云[7-8]。可燃物料質(zhì)量或體積一定的情況下,不同尺寸的氣云爆炸導(dǎo)致的后果嚴(yán)重程度不同,因此等價(jià)氣云的尺寸對(duì)爆炸后果具有重大影響。

    1 火氣系統(tǒng)對(duì)氣云的探測(cè)現(xiàn)狀

    火氣系統(tǒng)是針對(duì)可燃和有毒氣體泄漏及燃燒等具有探測(cè)、報(bào)警并能采取一定控制措施的完整的安全聯(lián)鎖系統(tǒng),探測(cè)系統(tǒng)檢測(cè)到氣體泄漏,邏輯運(yùn)算系統(tǒng)接收到探測(cè)系統(tǒng)信號(hào),計(jì)算等價(jià)氣云體積或濃度,若超過初始設(shè)定閾值,則信號(hào)傳遞到報(bào)警系統(tǒng)發(fā)出報(bào)警,并啟動(dòng)消防聯(lián)動(dòng)控制系統(tǒng)采取控制措施,防止事故發(fā)生造成嚴(yán)重后果,以此實(shí)現(xiàn)對(duì)人員、裝置、設(shè)備、管道、建筑等在危險(xiǎn)環(huán)境下的監(jiān)測(cè)及保護(hù)作用[9]。

    調(diào)查結(jié)果顯示,氣體泄漏是大多數(shù)火災(zāi)爆炸事故的誘因,因此對(duì)可燃?xì)怏w泄漏進(jìn)行探測(cè)并及時(shí)采取防控措施,對(duì)于火災(zāi)爆炸事故的預(yù)防以及保障正常生產(chǎn)活動(dòng)的順利進(jìn)行具有重要意義。目前已有多種標(biāo)準(zhǔn)規(guī)范對(duì)石化企業(yè)作業(yè)區(qū)的可燃?xì)怏w探測(cè)器以及火災(zāi)探測(cè)器做了明確規(guī)定,如GB 50493—2009、GB 50116—2013等[10-11],對(duì)探測(cè)點(diǎn)的布置、針對(duì)不同氣體的探測(cè)濃度范圍等做了規(guī)定,因此對(duì)FGS的設(shè)計(jì)需依據(jù)上述規(guī)范以及探測(cè)器的探測(cè)半徑,對(duì)探測(cè)器位置進(jìn)行定性布設(shè),并設(shè)定不同等級(jí)的報(bào)警閾值。

    火氣系統(tǒng)的有效性取決于探測(cè)器網(wǎng)絡(luò)、火氣硬件連鎖的安全可用性,以及觸發(fā)減緩系統(tǒng)的有效性[12],其中如何實(shí)現(xiàn)探測(cè)器網(wǎng)絡(luò)的定量化部署設(shè)計(jì)是源頭關(guān)鍵?,F(xiàn)有的研究從多個(gè)層面對(duì)探測(cè)器的定量布設(shè)進(jìn)行了優(yōu)化,如ISA(International Society of Automation)[13]以及 Defriend 等[14]和 Seo 等[15]從風(fēng)險(xiǎn)的角度對(duì)探測(cè)器位置進(jìn)行了優(yōu)化,Gimenes等[16]對(duì)場(chǎng)景采用網(wǎng)格數(shù)字化處理方法,章博等[17]針對(duì)具體場(chǎng)景分析了泄漏概率和后果,王海清等[12,18]、Legg等[19-20]利用數(shù)學(xué)規(guī)劃方法優(yōu)化了探測(cè)器布設(shè)網(wǎng)絡(luò)。

    本研究主要針對(duì)可燃?xì)庠频谋ㄊ鹿?,若氣云濃度在爆炸極限范圍內(nèi),遇到燃燒火源后極易發(fā)生爆炸,因此需要及時(shí)準(zhǔn)確地探測(cè)到可燃?xì)庠?。目前?duì)可燃?xì)庠频奶綔y(cè)主要基于氣體濃度,但是在氣體擴(kuò)散過程中,在非理想狀態(tài)下,可燃?xì)怏w與空氣混合不均勻時(shí),不同位置的探頭探測(cè)到的濃度不同,如依據(jù)濃度探測(cè),不能反映出不同泄漏場(chǎng)景氣云的爆炸載荷。

    為增加氣云探測(cè)的準(zhǔn)確性,本文將等價(jià)氣云理論應(yīng)用于火氣系統(tǒng)的探測(cè)輸入。目前已有研究將等價(jià)氣云的概念用于預(yù)測(cè)爆炸后果,但針對(duì)具體場(chǎng)景下受體載荷的氣云閾值尺寸尚未提出確定的算法。因此通過研究受體的爆炸載荷,利用氣云爆炸危害評(píng)估方法分析等價(jià)氣云的爆炸后果,氣云爆炸產(chǎn)生的沖擊波效應(yīng)達(dá)到受體載荷時(shí)對(duì)應(yīng)的氣云尺寸作為等價(jià)氣云尺寸閾值,進(jìn)而通過氣體泄漏擴(kuò)散模型得到擴(kuò)散臨界時(shí)間,并在氣云到達(dá)閾值尺寸之前采取措施,達(dá)到事故防控的目的。

    2 控制室爆炸載荷計(jì)算

    控制室是石化行業(yè)作業(yè)區(qū)的生產(chǎn)和安全的關(guān)鍵機(jī)構(gòu),是安全生產(chǎn)和操作指揮的中樞,對(duì)于維持設(shè)備正常運(yùn)轉(zhuǎn)、生產(chǎn)作業(yè)正常進(jìn)行、人員正常工作有至關(guān)重要的作用?;饸庀到y(tǒng)的核心裝置——控制系統(tǒng)通常設(shè)置在中心控制室內(nèi),因此控制室的安全是維持整個(gè)作業(yè)區(qū)安全及火氣系統(tǒng)的正常運(yùn)行的基礎(chǔ)。氣云爆炸后果在控制室載荷范圍內(nèi)時(shí)[21-22],可保證控制室不遭受破壞,維持火氣系統(tǒng)及作業(yè)區(qū)繼續(xù)運(yùn)行。本文以控制室能夠承受的爆炸載荷作為事故后果沖擊波閾值,通過事故后果計(jì)算方法逆推計(jì)算得到對(duì)應(yīng)的等價(jià)氣云尺寸。

    對(duì)控制室爆炸載荷的計(jì)算選取控制室設(shè)計(jì)防爆載荷的計(jì)算方法,在國(guó)內(nèi)外標(biāo)準(zhǔn)規(guī)范 ASCE 41088、GB 50779—2012[23-24]中均進(jìn)行了詳細(xì)闡述。

    若把控制室等效為封閉的矩形,作用在控制室前墻、后墻、側(cè)墻及屋面的載荷分別計(jì)算如下。

    (1)前墻載荷

    ① 峰值壓力

    式中,pSO為峰值入射超壓,kPa,應(yīng)根據(jù)石化裝置的性質(zhì)或裝置布局等因素進(jìn)行安全評(píng)估來確定,若未進(jìn)行評(píng)估,也可根據(jù)標(biāo)準(zhǔn)規(guī)范確定[12]。

    ② 停滯壓力

    式中,Cd為阻力系數(shù),取決于對(duì)象的表面形狀及朝向,對(duì)于控制室等矩形建筑物,前墻的阻力系數(shù)取+1.0,側(cè)墻及屋面、后墻取-0.4;q0為峰值動(dòng)壓,kPa,由式(3)計(jì)算得到

    式中,patm為環(huán)境標(biāo)準(zhǔn)大氣壓,kPa。

    ③ 前墻正壓等效作用時(shí)間

    式中,Iw為正壓沖量;S為停滯壓力點(diǎn)至建筑物邊緣的最小值,取H或B/2的較小值,其中H為建筑物高度,B為建筑物長(zhǎng)度;tc為反射壓持續(xù)時(shí)間,s;td為正壓作用時(shí)間,與峰值入射超壓一起經(jīng)安全評(píng)估得到或者依據(jù)標(biāo)準(zhǔn)規(guī)范確定,s;te為前墻正壓等效作用時(shí)間,s。

    (2)側(cè)墻及屋面載荷

    作用在側(cè)墻及屋面上的有效沖擊波超壓及其超壓時(shí)間按式(7)、式(8)計(jì)算

    式中,pa為作用在側(cè)墻及屋面上的有效沖擊波超壓,kPa;tr為側(cè)墻及屋面有效沖擊波超壓升壓時(shí)間,s;Cd為側(cè)墻及屋面阻力系數(shù),取-0.4;U為波速,m·s-1,由式(9)計(jì)算得到

    Ce為側(cè)墻及屋面載荷等效系數(shù),按Lw/L值查文獻(xiàn)[24]5.4節(jié)圖5.4.3,其中Lw(m)為沖擊波波長(zhǎng),由式 (10)計(jì)算得到,L(m)為控制室寬度。

    (3)后墻載荷

    作用在后墻上的有效沖擊波超壓及其作用時(shí)間

    式中,pb為作用在后墻上的有效沖擊波超壓,kPa;ta為沖擊波到達(dá)后墻時(shí)間,s;trb為后墻上的有效沖擊波超壓升壓時(shí)間,s;Cd為后墻阻力系數(shù),取-0.4;Ce為后墻載荷等效系數(shù),按Lw/L值查文獻(xiàn)[24]第5.4節(jié)圖5.4.3;D為沖擊波前進(jìn)方向建筑物寬度,后墻計(jì)算時(shí),取建筑物高度H,m。

    本文對(duì)控制室載荷的計(jì)算分析應(yīng)用于火氣系統(tǒng)探測(cè),為避免控制室受到事故損傷,探測(cè)時(shí)間應(yīng)盡量提前,因此在控制室設(shè)計(jì)載荷的基礎(chǔ)上,應(yīng)考慮載荷保守系數(shù)k。綜合考慮氣體擴(kuò)散速率、大氣環(huán)境影響等各方面因素,k取0.9,因此控制室載荷應(yīng)為

    式中,pload為控制室爆炸載荷,kPa;pcr為取前墻載荷ps或側(cè)墻載荷pa、后墻載荷pb的較小值,kPa。

    3 基于多能法的氣云爆炸后果計(jì)算

    針對(duì)在國(guó)內(nèi)外應(yīng)用爆炸破壞效應(yīng)的模型主要有多能法、TNT當(dāng)量法、SHELL法、CAM法等[25],其中多能法[26]被認(rèn)為是目前國(guó)內(nèi)外較先進(jìn)的適用于氣云爆炸后果預(yù)測(cè)的模型,綜合考慮了爆燃過程中的局部加速、障礙物的局部約束以及氣體活性等各種因素,其基本思想認(rèn)為只有在約束的條件下,氣云燃燒爆炸才能產(chǎn)生具有破壞性的沖擊波超壓,并根據(jù)氣云受到約束的程度確定相應(yīng)的爆炸源強(qiáng)度等級(jí),基于量綱1距離與量綱1峰值超壓得到氣云爆炸的峰值超壓。

    在已知控制室載荷求等價(jià)氣云體積的情況下,采用倒序逆推的方法,利用多能法進(jìn)行計(jì)算,步驟如下。

    (1)將控制室能夠承受的爆炸沖擊波載荷作為峰值超壓,根據(jù)式(15)可得到量綱1比擬峰值超壓

    式中,p'為量綱1比擬峰值超壓;pS為峰值超壓,即控制室載荷,kPa;patm為當(dāng)?shù)卮髿鈮毫?,一般可取?01.325 kPa。

    (2)采用Kinsella方法確定爆炸源強(qiáng)度的指數(shù)(1~10),從1到10依次增強(qiáng),爆炸源強(qiáng)度與氣云受到約束的程度有關(guān),1對(duì)應(yīng)無約束的氣云爆炸,10對(duì)應(yīng)氣云爆轟。3個(gè)主要影響因素分別為:點(diǎn)火能量大小、障礙物情況、約束情況。點(diǎn)火能量若小于100 MJ,則判定為弱,反之則為強(qiáng);障礙物擁塞度的判斷需依據(jù)兩個(gè)條件:一是障礙物體積大于整體擁擠區(qū)域的30%,二是障礙物之間的距離小于3 m,若同時(shí)滿足以下兩者則判定為強(qiáng),若不同時(shí)包含則判定為弱,若無障礙物則判定為不存在;約束存在代表氣云被2~3個(gè)固體表面限制,若只被地表1個(gè)固體表面限制則判定約束不存在。具體分級(jí)標(biāo)準(zhǔn)見文獻(xiàn)[26]。

    (3)根據(jù)步驟(1)得到的量綱1比擬峰值超壓與步驟(2)得到的爆炸源強(qiáng)度,依據(jù)超壓峰值與量綱1距離擬合公式或查文獻(xiàn)[26]中超壓峰值隨量綱 1距離變化曲線圖可得到對(duì)應(yīng)的量綱1距離R'。

    (4)計(jì)算得到當(dāng)量燃料——空氣混合物釋放的能量E(kJ)

    式中,R'為量綱1距離;z為控制室與爆炸中心的距離,m。

    (5)計(jì)算氣云體積。當(dāng)量燃料-空氣混合物所釋放的能量估算方法為:可燃物質(zhì)氣云體積與同化學(xué)計(jì)量濃度下烴-空氣混合物的典型燃燒熱值(約為3.5×103kJ·m-3)相乘,因此在控制室載荷下,氣云體積閾值Vc(m3)即為

    式中各參數(shù)意義同前文。

    (6)在給定的空間內(nèi),判斷氣云是否受到約束。

    式中,VF為場(chǎng)景內(nèi)自由體積,m3;VH為場(chǎng)景整體體積,m3;VE為場(chǎng)景內(nèi)設(shè)備體積,m3。若氣云體積Vc<VF,則認(rèn)為不存在不受約束的氣云,在該場(chǎng)景下氣云均對(duì)爆炸沖擊波產(chǎn)生貢獻(xiàn)。

    4 基于高斯擴(kuò)散模型的氣體泄漏時(shí)間計(jì)算

    根據(jù)等價(jià)氣云閾值體積,需分析氣體泄漏時(shí)間,為火氣系統(tǒng)探測(cè)時(shí)間提供參考。選取高斯擴(kuò)散模型,已知等價(jià)氣云尺寸閾值的情況下,通過高斯模型模擬得到泄漏時(shí)間[27-29]。

    高斯模型平均濃度方程為

    基于高斯擴(kuò)散模型,等價(jià)氣云的體積計(jì)算步驟如下。

    (1)網(wǎng)格劃分:建立氣體泄漏場(chǎng)景,將待分析區(qū)域劃分網(wǎng)格。網(wǎng)格密度可自定義,精度越高,計(jì)算越精確,但耗時(shí)較長(zhǎng),因此需要根據(jù)模擬場(chǎng)景的復(fù)雜程度選擇合適的網(wǎng)格密度。本文選擇10×10×10網(wǎng)格。

    (2)單元網(wǎng)格等價(jià)云體積計(jì)算

    ① 單元網(wǎng)格氣體平均體積濃度計(jì)算 利用高斯模型即式(19)計(jì)算得到網(wǎng)格內(nèi)8個(gè)頂點(diǎn)的濃度值,擬合單元網(wǎng)格體積濃度Cv關(guān)于(x,y,z)的四維方程計(jì)算得到氣體體積濃度均值

    ② 模擬場(chǎng)景內(nèi)等價(jià)氣云整體體積計(jì)算

    式中,Vgas為氣體泄漏擴(kuò)散形成的等價(jià)氣云體積,m3;i為網(wǎng)格數(shù)量,取值為1,2,3,…,n,其中n為網(wǎng)格整體數(shù)量;其他符號(hào)同上。

    綜上,在對(duì)可燃?xì)庠频葍r(jià)尺寸的閾值研究中,采用倒序逆推的方法,選取控制室為受體并計(jì)算其爆炸載荷作為爆炸后果,依次逆推得到等價(jià)氣云閾值尺寸以及火氣系統(tǒng)探測(cè)時(shí)間即氣體泄漏時(shí)間,整體分析流程如圖1所示。

    圖1 整體分析流程Fig.1 Overall analysis process

    5 案例計(jì)算分析

    5.1 案例簡(jiǎn)介

    某液化天然氣(LNG)罐區(qū)為梯形,平面尺寸為:上底180 m,下底240 m,寬度為280 m。左側(cè)為兩個(gè)液態(tài)LNG儲(chǔ)罐,容量為1.6×105m3,裝量系數(shù)不超過80%,底面直徑82 m,高度40 m;右側(cè)為控制室,尺寸為48 m×24 m×4 m;管廊高度約為 10 m。罐區(qū)安全距離設(shè)置需根據(jù) HG/T 20508—2014等標(biāo)準(zhǔn)[31],其中儲(chǔ)罐中心與控制室距離為100 m。具體分布情況如圖2所示,以某一LNG儲(chǔ)罐底部管線破裂導(dǎo)致氣體泄漏為例進(jìn)行計(jì)算模擬。假定其他條件一定時(shí),分別研究泄漏量及風(fēng)速等泄漏條件改變對(duì)等價(jià)氣云尺寸及火氣系統(tǒng)探測(cè)臨界時(shí)間的影響。

    圖2 某LNG罐區(qū)儲(chǔ)罐與控制室分布Fig.2 Layout of storage tank and control room in LNG tank industry

    5.2 控制室爆炸載荷計(jì)算

    控制室載荷設(shè)計(jì)采用的峰值入射超壓及對(duì)應(yīng)的正壓作用時(shí)間,需根據(jù)裝置性質(zhì)及平面布局等因素進(jìn)行安全評(píng)估來確定,未進(jìn)行安全評(píng)估時(shí),應(yīng)根據(jù)GB 50779—2012等標(biāo)準(zhǔn)[24]來確定:沖擊波峰值入射超壓可取為21 kPa,正壓作用時(shí)間可取為100 ms。根據(jù)式(1)~式(13),控制室載荷計(jì)算見表1。

    表1 控制室載荷計(jì)算結(jié)果Table 1 Calculation results of control room load

    計(jì)算得到前墻載荷為22.4 kPa,后墻、側(cè)墻及屋面載荷為12.04 kPa,由于側(cè)墻及后墻載荷較低,因此在氣云爆炸后果的計(jì)算中應(yīng)根據(jù)側(cè)墻及后期載荷保守取值,由式(14)得到控制室載荷為10.84 kPa。

    5.3 氣云尺寸計(jì)算

    (1)量綱1比擬峰值超壓 以控制室載荷作為氣云的實(shí)際峰值超壓,根據(jù)式(15)計(jì)算得到量綱 1比擬峰值超壓為0.221。

    (2)爆炸源強(qiáng)度指數(shù)的確定 基于場(chǎng)景的以下因素來選擇:點(diǎn)火能、爆炸源內(nèi)擁塞度、約束,見表2。

    該場(chǎng)景下選取點(diǎn)火能強(qiáng)度小于100 MJ,判斷點(diǎn)火能為弱;罐區(qū)整體作為擁擠區(qū)域,體積為5.88×105m3(高度取管廊高度為10 m),設(shè)備(儲(chǔ)罐與控制室)體積為1.102×105m3,得到擁塞度為18.74%,判定擁塞程度為弱。由于點(diǎn)火能——弱,擁塞度——弱,約束——存在,根據(jù)文獻(xiàn)[26]得到爆炸源強(qiáng)度等級(jí)為3~5,保守取值為5。

    表2 爆炸源強(qiáng)度指數(shù)影響因素Table 2 Influence factors of explosion source intensity index

    爆炸源強(qiáng)度為5時(shí),量綱1超壓與量綱1距離的曲線擬合方程

    (3)氣云體積閾值計(jì)算 根據(jù)量綱 1比擬峰值超壓與爆炸源強(qiáng)度等級(jí),依據(jù)式(22)或文獻(xiàn)[26]得到量綱1距離為1.11。根據(jù)式(16)、式(17)得到爆炸能量為 7.33×107kJ,對(duì)應(yīng)的氣云體積為2.09×104m3。

    5.4 判斷氣云是否受到約束

    氣云體積2.09×104m3小于擁擠區(qū)域內(nèi)自由體積4.78×105m3,不存在不受約束的氣云,該場(chǎng)景下的氣云均對(duì)爆炸沖擊波產(chǎn)生貢獻(xiàn),因此所求氣云尺寸即為可燃?xì)庠票ê罂刂剖逸d荷能夠承受的閾值尺寸。

    5.5 泄漏時(shí)間計(jì)算

    根據(jù)高斯擴(kuò)散模型及等價(jià)氣云體積計(jì)算模型即式(19)~式(21),應(yīng)用MATLAB軟件,進(jìn)行計(jì)算機(jī)模擬,得到如圖3、圖4所示的等價(jià)氣云體積隨時(shí)間變化曲線。圖3中,假定風(fēng)速恒為2 m·s-1,泄漏量分別為3×105~7×105kg時(shí),對(duì)應(yīng)探測(cè)臨界時(shí)間由 28~44 s逐漸增大;圖 4中,假定泄漏量為3×105kg,當(dāng)風(fēng)速分別由2~5 m·s-1變化時(shí),對(duì)應(yīng)探測(cè)臨界時(shí)間由17~44 s逐漸增大??梢?,受體爆炸載荷一定時(shí),針對(duì)不同的泄漏場(chǎng)景,火氣系統(tǒng)的響應(yīng)時(shí)間隨之改變。

    在上述場(chǎng)景下,氣體擴(kuò)散速度較快,火氣系統(tǒng)需要在臨界時(shí)間內(nèi),探測(cè)到氣云并及時(shí)采取相應(yīng)的防控措施,方可將氣云泄漏導(dǎo)致爆炸的事故后果限制在控制室爆炸載荷范圍內(nèi),實(shí)現(xiàn)對(duì)控制室保護(hù)的目的。氣云易受到氣象、地形、障礙物等各種因素的影響而發(fā)生改變,將其折算為等價(jià)氣云并作為火氣系統(tǒng)的設(shè)計(jì)輸入量,具有重要的工程意義。

    圖3 泄漏量改變時(shí)等價(jià)氣云體積-時(shí)間曲線Fig.3 Volume-time curve of equivalent gas cloud when leakage quality changes

    圖4 風(fēng)速改變時(shí)等價(jià)氣云體積-時(shí)間曲線Fig.4 Volume-time curve of equivalent gas cloud when wind speed changes

    6 結(jié) 論

    利用火氣系統(tǒng)對(duì)氣體泄漏及與空氣混合形成可燃?xì)庠频倪^程進(jìn)行探測(cè),可及時(shí)檢測(cè)到可燃?xì)怏w泄漏并采用防控措施。針對(duì)可燃?xì)庠瓢l(fā)生爆炸事故的復(fù)雜性與多變性,將其折算為等價(jià)氣云尺寸,并作為火氣系統(tǒng)的設(shè)計(jì)輸入。等價(jià)氣云概念為火氣系統(tǒng)的探測(cè)提供了理論依據(jù),等價(jià)氣云尺寸閾值的提出為火氣系統(tǒng)的應(yīng)用提供了技術(shù)設(shè)計(jì)指標(biāo),依據(jù)閾值尺寸得到探測(cè)臨界時(shí)間亦可為火氣系統(tǒng)設(shè)計(jì)提供理論支持,具有一定的工程實(shí)踐價(jià)值。通過對(duì)某LNG罐區(qū)控制室爆炸載荷計(jì)算,得到該罐區(qū)的等價(jià)氣云尺寸閾值以及臨界時(shí)間,證明了提出的方法具有合理性與工程應(yīng)用價(jià)值。

    [1]DíAZ A F,GONZáLEZ F E,SáNCHEZ PéREZ J F,et al.Characteristic overpressure-impulse-distance curves for vapour cloud explosions using the TNO multi-energy model[J].Journal of Hazardous Materials,2006,137(2):734-741.

    [2]何思念,常華偉,文科,等.液化天然氣泄漏和水面擴(kuò)散過程模擬[J].化工學(xué)報(bào),2015,66(S2):213-219.HE S N,CHANG H W,WEN K,et al.Mathematical modeling of LNG spills and pool spreading on water[J].CIESC Journal,2015,66(S2):213-219.

    [3]李兆慈,張文花,鄭梅,等.LNG連續(xù)泄漏擴(kuò)散過程模擬[J].化工學(xué)報(bào),2015,66(S2):158-165.LI Z C,ZHANG W H,ZHENG M,et al.LNG continuous leakage and diffusion process simulation[J].CIESC Journal,2015,66(S2):158-165.

    [4]ANGAL R,DEWAN A,SUBRAMANIAN K A.Computational study of turbulent hydrogen dispersion hazards in a closed space[J].IUP Journal of Mechanical Engineering,2012,5(2):28-42.

    [5]MIDDHA P,ICHARD M,ARNTZEN B J.Validation of CFD modelling of LH2,spread and evaporation against large-scale spill experiments[J].International Journal of Hydrogen Energy,2011,36(3):2620-2627.

    [6]BJERKETVEDT D,BAKKE J R,WINGERDEN K V.Gas explosion handbook[J].Journal of Hazardous Materials,1997,52(1):1-150.

    [7]HANSEN O R,GAVELLI F,DAVIS S G,et al.Equivalent cloud methods used for explosion risk and consequence studies[J].Journal of Loss Prevention in the Process Industries,2013,26(3):511-527.

    [8]ZHANG Y Q,OLDENBURG C M,PAN L H.Fast estimation of dense gas dispersion from multiple continuous CO2,surface leakage sources for risk assessment[J].International Journal of Greenhouse Gas Control,2016,49:323-329.

    [9]International Electrotechnical Commission.Functional Safety-Safety Instrumented Systems for the Process Industry Sector:IEC 61511—2003[S].International Electrotechnical Commission,2003.

    [10]中華人民共和國(guó)建設(shè)部.石油化工可燃?xì)怏w和有毒氣體檢測(cè)報(bào)警設(shè)計(jì)規(guī)范:GB 50493—2009[S].北京:中國(guó)計(jì)劃出版社,2009.Ministry of Construction of the People’s Republic of China.Code for design of combustible gas and toxic gas detection and alarm for petrochemical industry:GB 50493—2009[S].Beijing:China Planning Press,2009.

    [11]中華人民共和國(guó)住房和城鄉(xiāng)建設(shè)部.火災(zāi)自動(dòng)報(bào)警系統(tǒng)設(shè)計(jì)規(guī)范:GB 50116—2013[S].北京:中國(guó)計(jì)劃出版社,2013.Ministry of Housing and Urban-Rural Development of the People’s Republic of China.Code for design of automatic fire alarm system:GB 50116—2013[S].Beijing:China Planning Press,2013.

    [12]王海清,王玉秀,宋賢生,等.采用遺傳算法優(yōu)化氣體探測(cè)器布局策略研究[J].中國(guó)安全生產(chǎn)科學(xué)技術(shù),2016,12(9):86-91.WANG H Q,WANG Y X,SONG X S,et al.Study on layout optimization strategy of gas detectors using genetic algorithm[J].China Safety Science and Technology,2016,12(9):86-91.

    [13]Technical Report Guidance on the evaluation of fire,combustible gas and toxic system effectiveness:ISA TR84.00.07—2010[S].

    [14]DEFRIEND S,DEJMEK M,PORTER L,et al.A risk-based approach to flammable gas detector spacing[J].Journal of Hazardous Materials,2008,159(1):142-151.

    [15]SEO J K,DU C K,HA Y C,et al.A methodology for determining efficient gas detector locations on offshore installations[J].Ships &Offshore Structures,2013,8(5):524-535.

    [16]GIMENES G B C,SANTOS P C D,SOUZA V C D,et al.Development of a methodology for toxic and flammable gases sensors positioning in oil platforms[C]// AIChE Spring National Meeting,2008.

    [17]章博,王磊,王志剛.煉油裝置有害氣體泄漏區(qū)域風(fēng)險(xiǎn)等級(jí)劃分[J].中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,39(5):144-149.ZHANG B,WANG L,WNAG Z G.Area risk level classification for hazardous gas release in petroleum refining installations[J].Journal of China University of Petroleum (Edition of Natural Science),2015,39(5):144-149.

    [18]尚勝美,王海清,劉俊芳.FPSO火災(zāi)氣體探測(cè)器的優(yōu)化布置研究[J].石油與天然氣化工,2013,42(3):311-315.SHANG S M,WANG H Q,LIU J F.Study on optimized layout of fire and gas detectors in FPSO[J].Chemical Engineering of Oil and Gas,2013,42(3):311-315.

    [19]LEGG S W,WANG C,BENAVIDES-SERRANO A J,et al.Optimal gas detector placement under uncertainty considering Conditional-Value-at-Risk[J].Journal of Loss Prevention in the Process Industries,2013,26(3):410-417.

    [20]LEGG S W,BENAVIDES-SERRANO A J,SIIROLA J D,et al.A stochastic programming approach for gas detector placement using CFD-based dispersion simulations[J].Computers & Chemical Engineering,2012,47(12):194-201.

    [21]HANSEN O R,KJELLANDER M T,MARTINI R,et al.Estimation of explosion loading on small and medium sized equipment from CFD simulations[J].Journal of Loss Prevention in the Process Industries,2016,41:382-398.

    [22]AHMED I,BENGHERBIA T,ZHVANSKY R,et al.Validation of geometry modelling approaches for offshore gas dispersion simulations[J].Journal of Loss Prevention in the Process Industries,2016,44:594-600.

    [23]ASCE Petrochemical Energy Committee.Design of Blast Resistant Buildings in Petrochemical Facilities[S].American Society of Civil Engineers,2010.

    [24]中華人民共和國(guó)住房和城鄉(xiāng)建設(shè)部.石油化工控制室抗爆設(shè)計(jì)規(guī)范:GB 50779—2012[S].北京:中國(guó)計(jì)劃出版社,2012.Ministry of Housing and Urban-Rural Development of the People’s Republic of China.Code for design of blast resistant control building in petrochemical industry:GB 50779—2012[S].Beijing:China Planning Press,2012.

    [25]宇德明.易燃、易爆、有毒危險(xiǎn)品儲(chǔ)運(yùn)過程定量風(fēng)險(xiǎn)評(píng)價(jià)[M].北京:中國(guó)鐵道出版社,2000.YU D M.Quantitative Risk Assessment of Flammable,Explosive and Dangerous Goods Storage and Transportation Process[M].Beijing:China Railway Press,2000.

    [26]BOSCH C J H,WETERINGS R A P M.Methods for the Calculation of Physical Effects (Yellow Book)[M].The Hague:Committee for the Prevention of Disasters,2005.

    [27]FALLAH-SHORSHANI M,SHEKARRIZFARD M,HATZOPOULOU M.Integrating a street-canyon model with a regional Gaussian dispersion model for improved characterisation of near-road air pollution[J].Atmospheric Environment,2017,153(3):21-31.

    [28]唐建峰,蔡娜,郭清,等.LNG垂直噴射源連續(xù)泄漏擴(kuò)散的模擬[J].化工學(xué)報(bào),2013,64(3):1124-1131.TANG J F,CAI N,GUO Q,et al.Simulation of LNG diffusion:a continuous vertical jet release[J].CIESC Journal,2013,64(3):1124-1131.

    [29]FIATES J,VIANNA S S V.Numerical modelling of gas dispersion using OpenFOAM[J].Process Safety & Environmental Protection,2016,104:277-293.

    [30]國(guó)家安全生產(chǎn)監(jiān)督管理總局.化工企業(yè)定量風(fēng)險(xiǎn)評(píng)價(jià)導(dǎo)則:AQ/T 3046—2013[S].北京:煤炭工業(yè)出版社,2013.State Administration of Work Safety.Guidelines for quantitative risk assessment of chemical enterprises:AQ/T 3046—2013[S].Beijing:China Coal Industry Publishing House,2013.

    [31]中華人民共和國(guó)工業(yè)和信息化部.控制室設(shè)計(jì)規(guī)范:HG/T 20508—2014[S].北京:石油工業(yè)出版社,2014.Ministry of Industry and Information Technology of the People’s Republic of China.Control room design specification:HG/T 20508—2014[S].Beijing:Petroleum Industry Press,2014.

    date:2017-04-24.

    Prof.WANG Haiqing,wanghaiqing@upc.edu.cn

    supported by the National Science and Technology Major Project of the Ministry of Science and Technology of China(D719-ZGSY-555) and the Central Universities Fundamental Research Funds Project (17CX06004).

    Threshold of equivalent gas cloud size based on explosion load of control room

    QI Xinge1,WANG Haiqing1,SONG Xiansheng2,CHEN Guoming1
    (1Department of Safety Science and Engineering,China University of Petroleum,Qingdao266580,Shandong,China;2Lloyd’s Register Consulting-Energy Incorporated,Shanghai200001,China)

    Fire alarm and gas detector system (FGS) is the important safety barrier to prevent gas leakage.Combustible gas leaks and forms combustible gas cloud mixing with air,and its explosion belongs to volume explosion.The gas cloud has complexity and variability subject to a variety of factors,so it is converted into equivalent gas cloud and the threshold size calculation method is proposed.The threshold is key input index to realize the quantitative layout of detector network of FGS.The control room is selected as a receptor for analysis.Based on the critical value of shock wave overpressure which takes load of control room,the corresponding equivalent gas cloud size is calculated as the detecting threshold of combustible gas cloud by using multi-energy method in reverse order.Then the critical detecting time of the FGS is calculated by using equivalent air cloud computing method and Gaussian diffusion model.Through one LNG tank industry case analysis,the maximum cloud size which can be carried by the industry and the critical time of diffusion can be quantitatively determined.The numerical calculation shows that the equivalent gas cloud size threshold can not only be used as the quantitative input indicators for detector design of FGS,but also provide theoretical support for the detection time setting and the prevention and control measures of gas leakage and explosion.

    gas cloud; explosion; fire alarm and gas detector system; explosive load; equivalent gas cloud size;threshold size; safety; simulation

    X 937

    A

    0438—1157(2017)12—4857—08

    10.11949/j.issn.0438-1157.20170458

    2017-04-24收到初稿,2017-08-28收到修改稿。

    聯(lián)系人:王海清。

    齊心歌(1991—),女,博士研究生。

    國(guó)家重大科技專項(xiàng)子課題項(xiàng)目(D719-ZGSY-555);中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金項(xiàng)目(17CX06004)。

    猜你喜歡
    火氣控制室等價(jià)
    進(jìn)入控制室
    消防控制室 “假人”值班被處罰
    消防界(2021年19期)2021-11-14 12:14:40
    消防控制室 “假人”值班被處罰
    吃出“火氣”怎么辦? 6個(gè)食療小妙招巧應(yīng)對(duì)
    健康月歷
    健康之家(2020年7期)2020-11-02 02:18:28
    遠(yuǎn)離負(fù)面情緒 少些暴躁和火氣
    n次自然數(shù)冪和的一個(gè)等價(jià)無窮大
    中文信息(2017年12期)2018-01-27 08:22:58
    控制室在化工廠總平面布置中的安全研究
    張火丁:莫道染塵煙火氣,依然春閨夢(mèng)里人
    金色年華(2016年2期)2016-02-28 01:38:46
    收斂的非線性迭代數(shù)列xn+1=g(xn)的等價(jià)數(shù)列
    十八禁网站免费在线| 麻豆国产av国片精品| 性欧美人与动物交配| 天堂网av新在线| 国产精品98久久久久久宅男小说| 亚洲成人精品中文字幕电影| 变态另类丝袜制服| 男女视频在线观看网站免费| 亚洲男人的天堂狠狠| 中文字幕av在线有码专区| 国产极品精品免费视频能看的| 精品久久久久久,| 精品人妻一区二区三区麻豆 | 九九热线精品视视频播放| 给我免费播放毛片高清在线观看| 国产高清三级在线| 伦理电影大哥的女人| 自拍偷自拍亚洲精品老妇| 精品午夜福利在线看| 嫩草影院入口| 亚洲天堂国产精品一区在线| 精品国产亚洲在线| 亚洲成人久久性| 久久精品夜夜夜夜夜久久蜜豆| 国产私拍福利视频在线观看| 岛国在线免费视频观看| 又爽又黄无遮挡网站| 中文资源天堂在线| 国产三级中文精品| 精品久久国产蜜桃| 99久久精品国产亚洲精品| 日日摸夜夜添夜夜添小说| av在线蜜桃| 美女高潮喷水抽搐中文字幕| 真人一进一出gif抽搐免费| 九色国产91popny在线| 亚洲精品一卡2卡三卡4卡5卡| 午夜福利欧美成人| 岛国在线免费视频观看| 蜜桃亚洲精品一区二区三区| 好看av亚洲va欧美ⅴa在| 日本黄色视频三级网站网址| 熟妇人妻久久中文字幕3abv| 亚洲天堂国产精品一区在线| 国产探花极品一区二区| 亚洲va日本ⅴa欧美va伊人久久| 久久久久免费精品人妻一区二区| 精品福利观看| 国产极品精品免费视频能看的| 91午夜精品亚洲一区二区三区 | 国产高潮美女av| 成人毛片a级毛片在线播放| 午夜影院日韩av| 国产真实伦视频高清在线观看 | 欧美乱色亚洲激情| 国产亚洲精品久久久久久毛片| av在线老鸭窝| 国产又黄又爽又无遮挡在线| 欧美极品一区二区三区四区| 级片在线观看| 国产精品不卡视频一区二区 | a级一级毛片免费在线观看| 精品日产1卡2卡| 国产视频内射| 欧美xxxx性猛交bbbb| 毛片一级片免费看久久久久 | 极品教师在线免费播放| 久久午夜福利片| 级片在线观看| 国产乱人视频| 嫩草影院新地址| 亚洲精品456在线播放app | 3wmmmm亚洲av在线观看| 久久久久国内视频| 久久国产精品人妻蜜桃| 日韩人妻高清精品专区| 亚洲av日韩精品久久久久久密| 99热6这里只有精品| 小蜜桃在线观看免费完整版高清| av中文乱码字幕在线| 欧美性感艳星| 精品久久久久久久久久久久久| 午夜两性在线视频| 美女黄网站色视频| 午夜精品在线福利| 成人一区二区视频在线观看| 最好的美女福利视频网| 观看美女的网站| 亚洲欧美精品综合久久99| 亚洲国产精品合色在线| 久久中文看片网| 美女cb高潮喷水在线观看| 成人美女网站在线观看视频| 此物有八面人人有两片| 黄色视频,在线免费观看| 国产视频内射| 变态另类成人亚洲欧美熟女| 哪里可以看免费的av片| 欧美xxxx性猛交bbbb| 国产精品不卡视频一区二区 | 最近最新中文字幕大全电影3| 亚洲成人免费电影在线观看| 久久久久久久午夜电影| 亚洲av电影不卡..在线观看| 麻豆久久精品国产亚洲av| 亚洲熟妇熟女久久| 少妇人妻精品综合一区二区 | 亚洲av日韩精品久久久久久密| 国产日本99.免费观看| 国产午夜精品论理片| 97人妻精品一区二区三区麻豆| 女人十人毛片免费观看3o分钟| 亚洲国产精品久久男人天堂| 麻豆久久精品国产亚洲av| 欧美一区二区国产精品久久精品| 蜜桃亚洲精品一区二区三区| 成年女人看的毛片在线观看| 成年女人永久免费观看视频| 99热这里只有精品一区| 熟女电影av网| 久久人人爽人人爽人人片va | 身体一侧抽搐| 欧美一区二区国产精品久久精品| 婷婷丁香在线五月| 我的老师免费观看完整版| 成人特级av手机在线观看| 最新在线观看一区二区三区| 在线观看66精品国产| 一级作爱视频免费观看| 91九色精品人成在线观看| 深夜a级毛片| 91久久精品国产一区二区成人| 国产大屁股一区二区在线视频| 午夜福利在线观看免费完整高清在 | 亚洲国产欧美人成| 久久99热6这里只有精品| 一个人看视频在线观看www免费| eeuss影院久久| 亚洲欧美日韩卡通动漫| 精品一区二区三区av网在线观看| av中文乱码字幕在线| 日韩av在线大香蕉| 最近最新中文字幕大全电影3| 国产高清视频在线播放一区| 少妇裸体淫交视频免费看高清| 内地一区二区视频在线| 永久网站在线| av在线天堂中文字幕| 国产成人啪精品午夜网站| 999久久久精品免费观看国产| 精品久久久久久久久av| 欧美极品一区二区三区四区| 国产爱豆传媒在线观看| 日本 欧美在线| 免费在线观看日本一区| 国产野战对白在线观看| 久久久久九九精品影院| 亚洲一区二区三区不卡视频| 美女黄网站色视频| 老司机午夜十八禁免费视频| 国产精品亚洲一级av第二区| 麻豆国产97在线/欧美| 狂野欧美白嫩少妇大欣赏| 韩国av一区二区三区四区| 国产一区二区在线观看日韩| 亚洲欧美日韩高清专用| 国产精品99久久久久久久久| 国产私拍福利视频在线观看| 欧美中文日本在线观看视频| 成人高潮视频无遮挡免费网站| av在线观看视频网站免费| 国产精品久久久久久久久免 | 精品无人区乱码1区二区| 精品无人区乱码1区二区| 国产91精品成人一区二区三区| 噜噜噜噜噜久久久久久91| 精品一区二区三区视频在线| 中文字幕久久专区| 少妇人妻一区二区三区视频| 欧美色视频一区免费| 欧美+日韩+精品| 免费av毛片视频| 色精品久久人妻99蜜桃| 狠狠狠狠99中文字幕| 午夜激情福利司机影院| 啦啦啦观看免费观看视频高清| 免费看a级黄色片| 国产精品久久视频播放| 免费人成在线观看视频色| 亚洲国产精品sss在线观看| 男人的好看免费观看在线视频| 日韩欧美在线乱码| 色综合亚洲欧美另类图片| 午夜激情欧美在线| 欧美不卡视频在线免费观看| 九九热线精品视视频播放| 亚洲男人的天堂狠狠| 成人毛片a级毛片在线播放| 一区二区三区四区激情视频 | 深爱激情五月婷婷| 日韩欧美国产在线观看| 99国产极品粉嫩在线观看| 色精品久久人妻99蜜桃| 亚洲午夜理论影院| 一进一出好大好爽视频| 在线看三级毛片| 久久香蕉精品热| 亚洲av电影在线进入| 国产av一区在线观看免费| 午夜福利在线在线| 桃红色精品国产亚洲av| 日韩精品中文字幕看吧| 变态另类丝袜制服| 99热这里只有精品一区| 波多野结衣高清无吗| 欧美成人a在线观看| 国产伦精品一区二区三区四那| eeuss影院久久| 国产真实乱freesex| 欧美激情久久久久久爽电影| 乱人视频在线观看| 成人欧美大片| 深夜精品福利| 国产精品永久免费网站| 97人妻精品一区二区三区麻豆| 美女cb高潮喷水在线观看| 国产精品99久久久久久久久| 久久久久久大精品| 午夜久久久久精精品| 国产男靠女视频免费网站| 男人的好看免费观看在线视频| 国产成+人综合+亚洲专区| 黄片小视频在线播放| 欧美3d第一页| 老熟妇仑乱视频hdxx| 一个人观看的视频www高清免费观看| 国产精品野战在线观看| 草草在线视频免费看| 狂野欧美白嫩少妇大欣赏| .国产精品久久| 欧美性感艳星| av专区在线播放| 国产主播在线观看一区二区| 在线看三级毛片| 国产极品精品免费视频能看的| 两个人视频免费观看高清| 国产激情偷乱视频一区二区| 国产精品久久久久久精品电影| 精品无人区乱码1区二区| 国产精品久久久久久久电影| 夜夜爽天天搞| 日本熟妇午夜| av在线天堂中文字幕| 亚洲在线观看片| 国产一区二区在线观看日韩| 又爽又黄无遮挡网站| 国内少妇人妻偷人精品xxx网站| 男女视频在线观看网站免费| 人妻制服诱惑在线中文字幕| 舔av片在线| 亚洲欧美日韩高清在线视频| 久久精品国产清高在天天线| 老司机午夜福利在线观看视频| 国产高清视频在线播放一区| 欧美成狂野欧美在线观看| 网址你懂的国产日韩在线| 嫩草影院精品99| 久久久久久久亚洲中文字幕 | 又黄又爽又刺激的免费视频.| 亚洲国产精品sss在线观看| 国产美女午夜福利| 国产精品一区二区性色av| 男女视频在线观看网站免费| 亚洲国产精品合色在线| 精品免费久久久久久久清纯| 一本综合久久免费| 琪琪午夜伦伦电影理论片6080| av黄色大香蕉| www.色视频.com| 首页视频小说图片口味搜索| 人人妻人人澡欧美一区二区| 一进一出抽搐动态| 国产精品久久久久久久电影| 性色avwww在线观看| 少妇熟女aⅴ在线视频| 色综合站精品国产| 真实男女啪啪啪动态图| 国产主播在线观看一区二区| 亚洲人成网站在线播| 国产三级在线视频| bbb黄色大片| 欧美三级亚洲精品| 日本 欧美在线| 国产免费男女视频| 国产精品爽爽va在线观看网站| 欧美一级a爱片免费观看看| 99久久精品热视频| 欧美高清性xxxxhd video| 亚洲美女黄片视频| 免费在线观看成人毛片| 免费在线观看影片大全网站| 真人一进一出gif抽搐免费| 又紧又爽又黄一区二区| 国产高清有码在线观看视频| 亚洲va日本ⅴa欧美va伊人久久| 中文字幕熟女人妻在线| 可以在线观看毛片的网站| 色吧在线观看| 久久人妻av系列| 久久久成人免费电影| 成熟少妇高潮喷水视频| 成人三级黄色视频| 国产亚洲欧美98| 每晚都被弄得嗷嗷叫到高潮| 国产成人影院久久av| 成人国产综合亚洲| 男女之事视频高清在线观看| 亚洲一区二区三区不卡视频| 一级毛片久久久久久久久女| 久久久久久久久久成人| 搡老妇女老女人老熟妇| 搞女人的毛片| 欧美成人性av电影在线观看| 欧美乱色亚洲激情| 51午夜福利影视在线观看| 国产精品人妻久久久久久| 亚洲自偷自拍三级| 校园春色视频在线观看| 村上凉子中文字幕在线| 国产中年淑女户外野战色| 人妻夜夜爽99麻豆av| 国产亚洲av嫩草精品影院| 欧美激情国产日韩精品一区| 美女黄网站色视频| 最新在线观看一区二区三区| 亚洲第一区二区三区不卡| 久久久久久久午夜电影| x7x7x7水蜜桃| 日韩精品中文字幕看吧| 久久精品国产自在天天线| 十八禁人妻一区二区| 成人无遮挡网站| 国产在视频线在精品| 成人毛片a级毛片在线播放| 久久人人爽人人爽人人片va | 亚洲三级黄色毛片| 此物有八面人人有两片| 男女之事视频高清在线观看| 亚洲av不卡在线观看| 成年女人看的毛片在线观看| АⅤ资源中文在线天堂| 热99re8久久精品国产| 亚洲欧美日韩东京热| 亚洲av.av天堂| 一级黄片播放器| 欧美三级亚洲精品| 一夜夜www| 中文在线观看免费www的网站| 亚洲 国产 在线| 一进一出抽搐动态| 亚洲精品一区av在线观看| 精品无人区乱码1区二区| 成人国产一区最新在线观看| 两人在一起打扑克的视频| 亚洲精品在线观看二区| 国产中年淑女户外野战色| 国内少妇人妻偷人精品xxx网站| 国产淫片久久久久久久久 | a在线观看视频网站| 九九在线视频观看精品| 免费观看的影片在线观看| 日本黄色视频三级网站网址| 久9热在线精品视频| 精品一区二区三区人妻视频| 成人性生交大片免费视频hd| 久久欧美精品欧美久久欧美| 又爽又黄无遮挡网站| 国产真实伦视频高清在线观看 | 赤兔流量卡办理| 他把我摸到了高潮在线观看| 免费看日本二区| 国产一区二区三区视频了| 91在线精品国自产拍蜜月| 中文资源天堂在线| 少妇裸体淫交视频免费看高清| 国产精品99久久久久久久久| 欧美在线一区亚洲| 国产精品1区2区在线观看.| 日本 欧美在线| 欧美日韩综合久久久久久 | 国产视频内射| 欧美在线黄色| 色尼玛亚洲综合影院| 男女做爰动态图高潮gif福利片| 校园春色视频在线观看| 丰满的人妻完整版| 最近最新中文字幕大全电影3| 一区福利在线观看| 久99久视频精品免费| 757午夜福利合集在线观看| 69人妻影院| 国产成+人综合+亚洲专区| 综合色av麻豆| 欧美日韩国产亚洲二区| 亚洲欧美日韩卡通动漫| 久久香蕉精品热| 成年人黄色毛片网站| 两个人的视频大全免费| 美女 人体艺术 gogo| 男人舔奶头视频| 欧美丝袜亚洲另类 | 色综合站精品国产| 欧美激情国产日韩精品一区| 欧美黑人巨大hd| 在线a可以看的网站| 性插视频无遮挡在线免费观看| 在线观看免费视频日本深夜| 午夜精品久久久久久毛片777| 久久热精品热| 成年免费大片在线观看| 欧美在线黄色| 观看免费一级毛片| 精品久久久久久久久av| 在线观看舔阴道视频| 在线观看美女被高潮喷水网站 | 欧美性感艳星| 看十八女毛片水多多多| 亚洲成av人片免费观看| 91午夜精品亚洲一区二区三区 | 亚洲欧美日韩东京热| 天天躁日日操中文字幕| 国内精品美女久久久久久| 欧美xxxx性猛交bbbb| 午夜福利18| 国产精品久久久久久久久免 | 老司机午夜福利在线观看视频| 久久精品夜夜夜夜夜久久蜜豆| 亚洲天堂国产精品一区在线| 色哟哟·www| 亚洲国产日韩欧美精品在线观看| 岛国在线免费视频观看| 欧美黑人欧美精品刺激| 日日摸夜夜添夜夜添小说| a级毛片a级免费在线| 国产伦在线观看视频一区| 91久久精品国产一区二区成人| 国产黄片美女视频| 午夜精品在线福利| 1000部很黄的大片| 久久午夜亚洲精品久久| 999久久久精品免费观看国产| 亚洲成人久久爱视频| 欧美最新免费一区二区三区 | 国产欧美日韩一区二区精品| 婷婷精品国产亚洲av在线| 免费看日本二区| 国产精品久久久久久精品电影| 久久久精品大字幕| 国产精品自产拍在线观看55亚洲| 成人高潮视频无遮挡免费网站| 亚洲精品久久国产高清桃花| 在线观看av片永久免费下载| 极品教师在线视频| 久久久成人免费电影| 99热精品在线国产| 性色av乱码一区二区三区2| av天堂中文字幕网| 可以在线观看的亚洲视频| 日本与韩国留学比较| 日韩欧美在线二视频| 国产成年人精品一区二区| 日韩欧美在线乱码| 啪啪无遮挡十八禁网站| 12—13女人毛片做爰片一| 两人在一起打扑克的视频| 日本 av在线| 成人美女网站在线观看视频| 俄罗斯特黄特色一大片| 久久国产精品影院| 一进一出抽搐动态| 少妇丰满av| 国产aⅴ精品一区二区三区波| 国产精品日韩av在线免费观看| 欧美高清成人免费视频www| 岛国在线免费视频观看| 国内精品美女久久久久久| 亚洲成a人片在线一区二区| avwww免费| 国产探花在线观看一区二区| 性色av乱码一区二区三区2| 亚洲人成网站在线播放欧美日韩| bbb黄色大片| 欧美bdsm另类| 啦啦啦观看免费观看视频高清| 国内毛片毛片毛片毛片毛片| 久久香蕉精品热| 非洲黑人性xxxx精品又粗又长| 男女下面进入的视频免费午夜| 国产精品一区二区免费欧美| 国产精品女同一区二区软件 | 亚州av有码| 国产欧美日韩精品一区二区| 精品久久国产蜜桃| 国产精品野战在线观看| 亚洲内射少妇av| 欧美一级a爱片免费观看看| 深夜精品福利| 国产色爽女视频免费观看| 亚洲激情在线av| 精品久久久久久久久久久久久| 精品一区二区免费观看| 国内久久婷婷六月综合欲色啪| 国产伦精品一区二区三区视频9| 日本一二三区视频观看| av天堂中文字幕网| 国模一区二区三区四区视频| 女同久久另类99精品国产91| 十八禁人妻一区二区| 一区二区三区激情视频| 久久九九热精品免费| 欧美三级亚洲精品| 在现免费观看毛片| 成人毛片a级毛片在线播放| 国产精品美女特级片免费视频播放器| 亚洲专区中文字幕在线| 18禁黄网站禁片午夜丰满| 欧美又色又爽又黄视频| 国内精品久久久久久久电影| 村上凉子中文字幕在线| 极品教师在线视频| 热99在线观看视频| 亚洲久久久久久中文字幕| 亚洲综合色惰| 99久久99久久久精品蜜桃| 中国美女看黄片| 国产精品久久久久久亚洲av鲁大| 又爽又黄a免费视频| 久久久精品欧美日韩精品| 欧美区成人在线视频| 在线看三级毛片| 国产私拍福利视频在线观看| 日日摸夜夜添夜夜添av毛片 | 高清在线国产一区| 俄罗斯特黄特色一大片| 午夜福利在线观看免费完整高清在 | 午夜两性在线视频| 女人被狂操c到高潮| 99热只有精品国产| 天天一区二区日本电影三级| 别揉我奶头~嗯~啊~动态视频| 一a级毛片在线观看| 国产欧美日韩精品一区二区| 天天躁日日操中文字幕| 午夜福利在线在线| 九色国产91popny在线| 午夜亚洲福利在线播放| 两个人的视频大全免费| 亚洲久久久久久中文字幕| 丰满乱子伦码专区| 99视频精品全部免费 在线| 国产亚洲欧美98| 亚洲精品色激情综合| 一个人看视频在线观看www免费| 1024手机看黄色片| 国内揄拍国产精品人妻在线| 欧美高清成人免费视频www| 18禁裸乳无遮挡免费网站照片| 自拍偷自拍亚洲精品老妇| 国产aⅴ精品一区二区三区波| 免费观看人在逋| 身体一侧抽搐| 欧美日韩亚洲国产一区二区在线观看| 国产精品日韩av在线免费观看| 国产黄片美女视频| xxxwww97欧美| 国产成年人精品一区二区| 欧美日韩瑟瑟在线播放| 51国产日韩欧美| 我的女老师完整版在线观看| 韩国av一区二区三区四区| 欧美乱妇无乱码| 一a级毛片在线观看| 在线国产一区二区在线| 禁无遮挡网站| 少妇丰满av| 亚洲成av人片在线播放无| 给我免费播放毛片高清在线观看| 午夜福利欧美成人| 欧美高清成人免费视频www| 国产亚洲av嫩草精品影院| 久久精品国产自在天天线| 在线国产一区二区在线| 国产亚洲精品久久久com| 在线观看一区二区三区| 国产精品不卡视频一区二区 | 亚洲乱码一区二区免费版| 啦啦啦观看免费观看视频高清| 一个人观看的视频www高清免费观看| 午夜福利在线在线| 午夜精品一区二区三区免费看| 丝袜美腿在线中文| av天堂在线播放| 人妻丰满熟妇av一区二区三区| 日韩av在线大香蕉| 最好的美女福利视频网| 成人一区二区视频在线观看| 日韩欧美精品免费久久 | 亚洲不卡免费看| 欧美一区二区亚洲| 国产精品综合久久久久久久免费| 亚洲第一电影网av| 精品免费久久久久久久清纯| 欧美+亚洲+日韩+国产| 久久人妻av系列| 色播亚洲综合网|