• <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中文字幕在线| 日日摸夜夜添夜夜添小说| 在线观看免费高清a一片| 精品高清国产在线一区| 在线观看舔阴道视频| 久热爱精品视频在线9| 男人舔女人的私密视频| 国产精品久久久av美女十八| 日本精品一区二区三区蜜桃| 精品久久久久久久毛片微露脸| 99国产精品一区二区蜜桃av | 狠狠狠狠99中文字幕| 久久亚洲真实| 中文字幕人妻熟女乱码| 精品久久蜜臀av无| 欧美久久黑人一区二区| 亚洲欧美一区二区三区久久| 一级a爱视频在线免费观看| 久久久久久久国产电影| 精品国产一区二区三区四区第35| 国产精品免费一区二区三区在线 | 男女床上黄色一级片免费看| 叶爱在线成人免费视频播放| a在线观看视频网站| aaaaa片日本免费| 黄色丝袜av网址大全| 一区二区三区国产精品乱码| 下体分泌物呈黄色| www.自偷自拍.com| 他把我摸到了高潮在线观看 | 咕卡用的链子| 777久久人妻少妇嫩草av网站| 国产精品久久久久久精品古装| 丝袜喷水一区| 美女高潮喷水抽搐中文字幕| 久久人妻福利社区极品人妻图片| 亚洲五月婷婷丁香| 国产单亲对白刺激| 人人澡人人妻人| 五月天丁香电影| 老熟妇仑乱视频hdxx| 99re在线观看精品视频| 精品少妇黑人巨大在线播放| 黄频高清免费视频| 麻豆国产av国片精品| 国产成人欧美在线观看 | 不卡av一区二区三区| 老司机靠b影院| 99re在线观看精品视频| 日韩三级视频一区二区三区| 欧美在线一区亚洲| 三级毛片av免费| 2018国产大陆天天弄谢| 中文字幕av电影在线播放| 激情在线观看视频在线高清 | 丝袜喷水一区| 亚洲精品粉嫩美女一区| 国产成人精品久久二区二区91| 亚洲天堂av无毛| 欧美成人免费av一区二区三区 | 狂野欧美激情性xxxx| 久久久国产一区二区| xxxhd国产人妻xxx| 亚洲精品一卡2卡三卡4卡5卡| 天天躁夜夜躁狠狠躁躁| 一进一出好大好爽视频| 香蕉丝袜av| 久9热在线精品视频| 黑人猛操日本美女一级片| 丰满少妇做爰视频| 人妻 亚洲 视频| 亚洲中文av在线| 汤姆久久久久久久影院中文字幕| 老司机午夜福利在线观看视频 | 国产日韩欧美亚洲二区| 美女高潮喷水抽搐中文字幕| 日本a在线网址| 精品国产乱码久久久久久小说| √禁漫天堂资源中文www| 丰满少妇做爰视频| 中文字幕制服av| 考比视频在线观看| 极品教师在线免费播放| 午夜精品久久久久久毛片777| 午夜成年电影在线免费观看| 亚洲欧美一区二区三区久久| 婷婷成人精品国产| 成年动漫av网址| 国产免费现黄频在线看| 超色免费av| 精品国产一区二区久久| 淫妇啪啪啪对白视频| 五月天丁香电影| 亚洲少妇的诱惑av| 欧美黑人欧美精品刺激| 国产精品 国内视频| 最新在线观看一区二区三区| 黄色毛片三级朝国网站| 日韩一区二区三区影片| 日韩 欧美 亚洲 中文字幕| 国产精品.久久久| 国产色视频综合| 国产精品欧美亚洲77777| 天天躁日日躁夜夜躁夜夜| 久久久久久久大尺度免费视频| 1024香蕉在线观看| 亚洲国产精品一区二区三区在线| 热99re8久久精品国产| 久久久久久久久免费视频了| 欧美日韩av久久| 亚洲精品在线观看二区| 99在线人妻在线中文字幕 | 99热网站在线观看| 精品国内亚洲2022精品成人 | 五月天丁香电影| 国产精品影院久久| 国产欧美日韩精品亚洲av| av网站免费在线观看视频| 久久久精品国产亚洲av高清涩受| 久久热在线av| 成年人午夜在线观看视频| 日韩制服丝袜自拍偷拍| 亚洲av美国av| 男女下面插进去视频免费观看| 成人国产av品久久久| 国产xxxxx性猛交| 俄罗斯特黄特色一大片| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲欧洲精品一区二区精品久久久| 99国产精品一区二区三区| 国产精品影院久久| 日韩熟女老妇一区二区性免费视频| 久久性视频一级片| 亚洲精品国产色婷婷电影| 午夜激情av网站| 精品一区二区三卡| 精品国产乱码久久久久久男人| 国产成人一区二区三区免费视频网站| 久久99热这里只频精品6学生| 在线观看免费视频日本深夜| 精品视频人人做人人爽| 久久人妻福利社区极品人妻图片| 亚洲人成伊人成综合网2020| 日日爽夜夜爽网站| 老司机靠b影院| 91精品三级在线观看| 高清毛片免费观看视频网站 | 亚洲精品一二三| 国产精品 国内视频| 久久久国产一区二区| 欧美性长视频在线观看| 十八禁高潮呻吟视频| 日韩熟女老妇一区二区性免费视频| 欧美日韩一级在线毛片| 50天的宝宝边吃奶边哭怎么回事| 天堂俺去俺来也www色官网| 国产免费av片在线观看野外av| 电影成人av| 黄色丝袜av网址大全| 国产精品1区2区在线观看. | 男女午夜视频在线观看| 人人妻人人澡人人爽人人夜夜| 国产色视频综合| 黄片播放在线免费| 亚洲精品在线美女| 人妻一区二区av| 免费看十八禁软件| 在线播放国产精品三级| 黄色成人免费大全| avwww免费| 久久久久久久精品吃奶| 亚洲国产中文字幕在线视频| 老司机亚洲免费影院| 美女午夜性视频免费| 天堂中文最新版在线下载| 99热网站在线观看| 老熟女久久久| 丝袜喷水一区| 亚洲av电影在线进入| 日韩有码中文字幕| 精品少妇久久久久久888优播| 欧美亚洲日本最大视频资源| 人人妻人人爽人人添夜夜欢视频| 高清毛片免费观看视频网站 | 欧美激情 高清一区二区三区| 高清毛片免费观看视频网站 | videosex国产| 国产国语露脸激情在线看| 亚洲 欧美一区二区三区| 黄色 视频免费看| 日韩三级视频一区二区三区| 十八禁人妻一区二区| 淫妇啪啪啪对白视频| 亚洲精品国产区一区二| 十八禁网站免费在线| 免费观看a级毛片全部| 国产av国产精品国产| 三上悠亚av全集在线观看| 国产精品九九99| 午夜老司机福利片| 中文字幕高清在线视频| 男女午夜视频在线观看| 久久人妻熟女aⅴ| 大码成人一级视频| 在线 av 中文字幕| h视频一区二区三区| 最近最新免费中文字幕在线| 成人影院久久| 精品一区二区三区视频在线观看免费 | 99在线人妻在线中文字幕 | 国产精品免费大片| 在线观看66精品国产| 午夜福利,免费看| 成人黄色视频免费在线看| 国产99久久九九免费精品| 丝袜在线中文字幕| 女性生殖器流出的白浆| 女人精品久久久久毛片| 中文字幕精品免费在线观看视频| 巨乳人妻的诱惑在线观看| 纵有疾风起免费观看全集完整版| 亚洲欧洲日产国产| 亚洲欧美激情在线| 国产精品二区激情视频| 操出白浆在线播放| 99久久99久久久精品蜜桃| 日韩成人在线观看一区二区三区| 肉色欧美久久久久久久蜜桃| 两人在一起打扑克的视频| 757午夜福利合集在线观看| 一区二区日韩欧美中文字幕| 一边摸一边抽搐一进一小说 | 黑人操中国人逼视频| 国产男靠女视频免费网站| 国产伦人伦偷精品视频| 欧美中文综合在线视频| 国产亚洲一区二区精品| 老汉色av国产亚洲站长工具| 欧美日韩一级在线毛片| 久久精品国产综合久久久| 99re在线观看精品视频| 午夜福利在线观看吧| 我的亚洲天堂| www.自偷自拍.com| 中文字幕另类日韩欧美亚洲嫩草| 亚洲综合色网址| 亚洲五月色婷婷综合| www.999成人在线观看| 无人区码免费观看不卡 | 欧美精品人与动牲交sv欧美| 亚洲精品国产色婷婷电影| 老司机靠b影院| 久久中文字幕人妻熟女| 亚洲精品国产一区二区精华液| 亚洲成人免费电影在线观看| 啦啦啦 在线观看视频| 99热网站在线观看| 一二三四在线观看免费中文在| 99精品欧美一区二区三区四区| 欧美日韩一级在线毛片| 一本综合久久免费| 人妻 亚洲 视频| 成人黄色视频免费在线看| 精品国产国语对白av| 国产精品香港三级国产av潘金莲| 高清毛片免费观看视频网站 | 亚洲七黄色美女视频| 国产aⅴ精品一区二区三区波| 亚洲全国av大片| 亚洲精品中文字幕一二三四区 | 人人妻人人爽人人添夜夜欢视频| 99久久国产精品久久久| www日本在线高清视频| 久久天堂一区二区三区四区| 久久香蕉激情| 久久精品亚洲精品国产色婷小说| 精品人妻在线不人妻| 波多野结衣一区麻豆| 一区二区三区乱码不卡18| 80岁老熟妇乱子伦牲交| 亚洲色图 男人天堂 中文字幕| 色播在线永久视频| 亚洲欧美色中文字幕在线| 老汉色∧v一级毛片| 国产精品一区二区在线观看99| av网站在线播放免费| 亚洲中文字幕日韩| 99在线人妻在线中文字幕 | 欧美精品av麻豆av| 三上悠亚av全集在线观看| 久久中文字幕一级| 黑人猛操日本美女一级片| 久久久精品区二区三区| 午夜福利在线免费观看网站| 2018国产大陆天天弄谢| 黄网站色视频无遮挡免费观看| 性高湖久久久久久久久免费观看| 免费不卡黄色视频| 下体分泌物呈黄色| 国产精品二区激情视频| 午夜老司机福利片| 国产精品熟女久久久久浪| avwww免费| 国产av又大| 午夜91福利影院| 成年动漫av网址| 久久香蕉激情| 久久精品亚洲精品国产色婷小说| 精品人妻在线不人妻| 一进一出抽搐动态| 欧美 亚洲 国产 日韩一| 国产免费av片在线观看野外av| 一区二区av电影网| 亚洲一区中文字幕在线| 国产精品免费一区二区三区在线 | 免费人妻精品一区二区三区视频| 一级毛片精品| 啦啦啦免费观看视频1| 久久精品亚洲av国产电影网| 亚洲av片天天在线观看| 99re在线观看精品视频| 成年人黄色毛片网站| 美国免费a级毛片| 日韩欧美三级三区| aaaaa片日本免费| 日韩一区二区三区影片| 欧美日本中文国产一区发布| 久久天躁狠狠躁夜夜2o2o| 黄色a级毛片大全视频| 不卡一级毛片| 亚洲av美国av| 日韩一卡2卡3卡4卡2021年| 亚洲欧美精品综合一区二区三区| 日韩熟女老妇一区二区性免费视频| 国产亚洲欧美精品永久| 丝袜在线中文字幕| 交换朋友夫妻互换小说| 久久99热这里只频精品6学生| 美女午夜性视频免费| 国产av又大| 免费日韩欧美在线观看| 日本黄色日本黄色录像| 国产成人欧美在线观看 | 精品亚洲成a人片在线观看| 一级,二级,三级黄色视频| 一区二区三区国产精品乱码| 久久精品国产亚洲av香蕉五月 | 多毛熟女@视频| 韩国精品一区二区三区| 老鸭窝网址在线观看| 美女福利国产在线| 久久久精品免费免费高清| 久久久精品国产亚洲av高清涩受| 一二三四在线观看免费中文在| 女人久久www免费人成看片| 黑人猛操日本美女一级片| 韩国精品一区二区三区| 久久精品91无色码中文字幕| 亚洲欧美一区二区三区久久| kizo精华| 最新在线观看一区二区三区| 老汉色∧v一级毛片| 欧美在线一区亚洲| 国产成人精品久久二区二区免费| 另类亚洲欧美激情| 色婷婷av一区二区三区视频| 精品少妇一区二区三区视频日本电影| 国产精品一区二区在线观看99| √禁漫天堂资源中文www| 久久人人97超碰香蕉20202| 一级,二级,三级黄色视频| 精品国产国语对白av| 国产精品二区激情视频| 高清av免费在线| 男人舔女人的私密视频| 99riav亚洲国产免费| 国产精品久久久久成人av| 99在线人妻在线中文字幕 | 亚洲av成人一区二区三| 久久精品国产亚洲av香蕉五月 | 两个人看的免费小视频| 午夜福利视频在线观看免费| 日本av免费视频播放| 在线 av 中文字幕| 国产免费现黄频在线看| 国产精品一区二区在线不卡| 丝袜喷水一区| 中文字幕人妻熟女乱码| 欧美精品亚洲一区二区| 少妇裸体淫交视频免费看高清 | 国产日韩欧美在线精品| 在线观看免费午夜福利视频| 黑人巨大精品欧美一区二区mp4| 欧美日韩亚洲综合一区二区三区_| 久久热在线av| 成年动漫av网址| 国产精品秋霞免费鲁丝片| 精品国产亚洲在线| 满18在线观看网站| 热re99久久国产66热| 天堂动漫精品| 欧美另类亚洲清纯唯美| 人人妻,人人澡人人爽秒播| 亚洲五月婷婷丁香| 精品福利观看| 叶爱在线成人免费视频播放| 丁香六月欧美| 宅男免费午夜| 大片电影免费在线观看免费| 亚洲 国产 在线| 丝袜人妻中文字幕| 五月天丁香电影| 久久国产精品大桥未久av| 国产精品 欧美亚洲| 在线观看人妻少妇| 99热国产这里只有精品6| 亚洲专区国产一区二区| 亚洲av第一区精品v没综合| 中文字幕高清在线视频| 十八禁网站免费在线| 在线播放国产精品三级| 搡老岳熟女国产| 免费不卡黄色视频| 少妇粗大呻吟视频| 日韩大码丰满熟妇| 老熟妇乱子伦视频在线观看| 黑丝袜美女国产一区| 高潮久久久久久久久久久不卡| 一进一出好大好爽视频| 亚洲综合色网址| 国产野战对白在线观看| 国产欧美日韩综合在线一区二区| 亚洲三区欧美一区| 亚洲中文av在线| 精品第一国产精品| 日韩欧美一区视频在线观看| 在线av久久热| www.自偷自拍.com| 免费在线观看黄色视频的| 激情视频va一区二区三区| 母亲3免费完整高清在线观看| 免费在线观看视频国产中文字幕亚洲| 亚洲av电影在线进入| 精品国内亚洲2022精品成人 | 久久精品aⅴ一区二区三区四区| 一本色道久久久久久精品综合| 日韩免费av在线播放| 久久九九热精品免费| 国产精品一区二区精品视频观看| 亚洲七黄色美女视频| 美国免费a级毛片| 日韩熟女老妇一区二区性免费视频| 脱女人内裤的视频| 人妻一区二区av| 两个人看的免费小视频| 他把我摸到了高潮在线观看 | 国产av一区二区精品久久| 婷婷丁香在线五月| 日本五十路高清| 亚洲人成电影免费在线| 精品国产乱码久久久久久男人| 日韩中文字幕视频在线看片| 亚洲av成人不卡在线观看播放网| 午夜精品久久久久久毛片777| 精品国产超薄肉色丝袜足j| 最近最新免费中文字幕在线| 9热在线视频观看99| 久久人妻熟女aⅴ| 又黄又粗又硬又大视频| 午夜福利免费观看在线| 一级毛片女人18水好多| 亚洲熟妇熟女久久| xxxhd国产人妻xxx| 国产欧美日韩精品亚洲av| 精品第一国产精品| 国产一卡二卡三卡精品| 亚洲国产中文字幕在线视频| 日韩 欧美 亚洲 中文字幕| 日韩中文字幕视频在线看片| av天堂在线播放| 极品人妻少妇av视频| 久久精品国产a三级三级三级| 国产欧美日韩综合在线一区二区| 亚洲国产成人一精品久久久| tocl精华| 午夜精品国产一区二区电影| 波多野结衣一区麻豆| 欧美 日韩 精品 国产| 欧美成狂野欧美在线观看| 婷婷成人精品国产| 亚洲色图 男人天堂 中文字幕| 999久久久精品免费观看国产| 大香蕉久久成人网| 蜜桃国产av成人99| 成人黄色视频免费在线看| 999久久久精品免费观看国产| 国产精品成人在线| 国产成人欧美| 久久九九热精品免费| 精品一区二区三区视频在线观看免费 | 青青草视频在线视频观看| 两个人免费观看高清视频| av网站在线播放免费| 国产xxxxx性猛交| 国产激情久久老熟女| 两个人免费观看高清视频| www.熟女人妻精品国产| 1024香蕉在线观看| 欧美久久黑人一区二区| 久久久久精品人妻al黑| 99精品在免费线老司机午夜| av有码第一页| 色综合婷婷激情| 免费观看人在逋| 久久精品亚洲熟妇少妇任你| 国产av精品麻豆| 欧美国产精品va在线观看不卡| 久久99一区二区三区| 久久精品亚洲熟妇少妇任你| 男女下面插进去视频免费观看| 欧美精品亚洲一区二区| 超色免费av| 久久婷婷成人综合色麻豆| 国产欧美日韩综合在线一区二区| 欧美+亚洲+日韩+国产| 首页视频小说图片口味搜索| 亚洲欧美色中文字幕在线| 欧美国产精品一级二级三级| 一进一出好大好爽视频| xxxhd国产人妻xxx| 国产精品.久久久| 亚洲午夜理论影院| 亚洲国产精品一区二区三区在线| 色老头精品视频在线观看| 免费看十八禁软件| 动漫黄色视频在线观看| 99热国产这里只有精品6| 乱人伦中国视频| 少妇的丰满在线观看| 久久精品熟女亚洲av麻豆精品| 多毛熟女@视频| 99精品欧美一区二区三区四区| 国产精品电影一区二区三区 | 视频在线观看一区二区三区| 黄网站色视频无遮挡免费观看| 欧美人与性动交α欧美精品济南到| 三上悠亚av全集在线观看| 国产成人免费观看mmmm| 国产又爽黄色视频| 国产亚洲欧美精品永久| 久久天堂一区二区三区四区| 国产成人av激情在线播放| 日韩熟女老妇一区二区性免费视频| 丰满少妇做爰视频| 亚洲黑人精品在线| 亚洲中文字幕日韩| 中文字幕人妻熟女乱码| videos熟女内射| 最新在线观看一区二区三区| 午夜福利欧美成人| 久久ye,这里只有精品| 亚洲av美国av| 日日摸夜夜添夜夜添小说| 一个人免费在线观看的高清视频| 91九色精品人成在线观看| 中文字幕另类日韩欧美亚洲嫩草| 精品福利观看| 汤姆久久久久久久影院中文字幕| 999久久久精品免费观看国产| 如日韩欧美国产精品一区二区三区| 国产在线一区二区三区精| 久久精品国产99精品国产亚洲性色 | 99国产极品粉嫩在线观看| av网站在线播放免费| a级片在线免费高清观看视频| 日日摸夜夜添夜夜添小说| 香蕉丝袜av| 国产精品亚洲一级av第二区| 欧美人与性动交α欧美精品济南到| 99国产极品粉嫩在线观看| 午夜久久久在线观看| 高清视频免费观看一区二区| 黑人操中国人逼视频| 欧美精品一区二区大全| 国产精品久久久久久精品古装| 国产麻豆69| 国产主播在线观看一区二区| 黄色a级毛片大全视频| 露出奶头的视频| 激情在线观看视频在线高清 | 91麻豆av在线| 亚洲第一青青草原| 精品少妇内射三级| 国产精品久久久久久人妻精品电影 | 99国产精品免费福利视频| netflix在线观看网站| 久久精品国产亚洲av高清一级| 少妇 在线观看| 热99re8久久精品国产| 欧美av亚洲av综合av国产av| 免费av中文字幕在线| 大香蕉久久网| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品麻豆人妻色哟哟久久| 国产人伦9x9x在线观看| 国产三级黄色录像| 女性被躁到高潮视频| 久久性视频一级片| 国产亚洲精品第一综合不卡| 国产精品欧美亚洲77777| 国产一区二区激情短视频| 热99国产精品久久久久久7| 大香蕉久久成人网| 久久久国产欧美日韩av| 久久久精品94久久精品|