• <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ù)列
    日本五十路高清| 久久人妻熟女aⅴ| 久久精品亚洲精品国产色婷小说| 久久人妻福利社区极品人妻图片| netflix在线观看网站| 老鸭窝网址在线观看| 国产又色又爽无遮挡免费看| 午夜日韩欧美国产| 免费人成视频x8x8入口观看| 又紧又爽又黄一区二区| 丰满人妻熟妇乱又伦精品不卡| av网站免费在线观看视频| 波多野结衣av一区二区av| 好看av亚洲va欧美ⅴa在| 精品无人区乱码1区二区| 悠悠久久av| 国产97色在线日韩免费| 韩国精品一区二区三区| 热99re8久久精品国产| 日本五十路高清| 国产成+人综合+亚洲专区| 色婷婷av一区二区三区视频| 午夜福利,免费看| 777久久人妻少妇嫩草av网站| 视频区图区小说| 天堂√8在线中文| 日韩中文字幕欧美一区二区| x7x7x7水蜜桃| 黄色视频,在线免费观看| 999久久久国产精品视频| 久久久久久久久免费视频了| 69av精品久久久久久| 啦啦啦免费观看视频1| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品乱久久久久久| 国产精品 欧美亚洲| 多毛熟女@视频| 久久午夜亚洲精品久久| 免费在线观看完整版高清| 国产精品1区2区在线观看. | 亚洲一码二码三码区别大吗| 精品一区二区三卡| 亚洲午夜精品一区,二区,三区| 精品欧美一区二区三区在线| 69av精品久久久久久| 欧美黄色片欧美黄色片| 18禁观看日本| 欧美日韩国产mv在线观看视频| 99精品久久久久人妻精品| 99精品欧美一区二区三区四区| 久久狼人影院| 99re6热这里在线精品视频| 99国产综合亚洲精品| 久久性视频一级片| 成年女人毛片免费观看观看9 | a级毛片黄视频| 黑人巨大精品欧美一区二区mp4| 亚洲成人手机| 精品国内亚洲2022精品成人 | 看免费av毛片| 50天的宝宝边吃奶边哭怎么回事| 夜夜夜夜夜久久久久| 一进一出抽搐动态| 亚洲人成电影观看| 欧美精品一区二区免费开放| 99精国产麻豆久久婷婷| 国产欧美日韩一区二区三| 一区福利在线观看| 国产精品久久久av美女十八| 亚洲全国av大片| 麻豆成人av在线观看| 无限看片的www在线观看| 亚洲欧美日韩高清在线视频| 国产亚洲欧美精品永久| 无人区码免费观看不卡| 又大又爽又粗| 满18在线观看网站| 大陆偷拍与自拍| 高清视频免费观看一区二区| 一本一本久久a久久精品综合妖精| 成人av一区二区三区在线看| 高清欧美精品videossex| 啦啦啦视频在线资源免费观看| 国产精品久久电影中文字幕 | 亚洲成国产人片在线观看| 久久香蕉激情| 欧美亚洲 丝袜 人妻 在线| 久久人妻熟女aⅴ| 成人亚洲精品一区在线观看| 宅男免费午夜| 国产高清激情床上av| tube8黄色片| www.自偷自拍.com| 热re99久久精品国产66热6| 亚洲色图av天堂| 无限看片的www在线观看| 真人做人爱边吃奶动态| 精品国产乱子伦一区二区三区| 亚洲综合色网址| 欧美丝袜亚洲另类 | 中文字幕av电影在线播放| 午夜精品久久久久久毛片777| 精品免费久久久久久久清纯 | 久9热在线精品视频| 日本a在线网址| 黑丝袜美女国产一区| 亚洲久久久国产精品| 久久久精品免费免费高清| 国产精品影院久久| 欧美 日韩 精品 国产| 亚洲美女黄片视频| 搡老岳熟女国产| 捣出白浆h1v1| 国产精品免费一区二区三区在线 | 伊人久久大香线蕉亚洲五| 热99久久久久精品小说推荐| 国产成人精品无人区| 黄色片一级片一级黄色片| 日本撒尿小便嘘嘘汇集6| 成人18禁在线播放| svipshipincom国产片| 法律面前人人平等表现在哪些方面| 亚洲色图 男人天堂 中文字幕| 69av精品久久久久久| 亚洲七黄色美女视频| 国产伦人伦偷精品视频| 俄罗斯特黄特色一大片| 午夜福利免费观看在线| 久久久久久久午夜电影 | 日韩制服丝袜自拍偷拍| 悠悠久久av| 在线观看www视频免费| 国产精品亚洲一级av第二区| 国产乱人伦免费视频| 美女国产高潮福利片在线看| 又黄又粗又硬又大视频| 精品国产国语对白av| 成人精品一区二区免费| 99国产精品99久久久久| 我的亚洲天堂| 我的亚洲天堂| 国产99久久九九免费精品| 久久精品国产综合久久久| 男女高潮啪啪啪动态图| 欧美成人免费av一区二区三区 | 久久久精品国产亚洲av高清涩受| 日韩免费av在线播放| 国产乱人伦免费视频| 欧美乱色亚洲激情| 国产日韩一区二区三区精品不卡| 一区二区三区国产精品乱码| 午夜成年电影在线免费观看| 首页视频小说图片口味搜索| 久99久视频精品免费| 青草久久国产| 亚洲一码二码三码区别大吗| 90打野战视频偷拍视频| 99在线人妻在线中文字幕 | 最新的欧美精品一区二区| 大陆偷拍与自拍| 五月开心婷婷网| 天堂动漫精品| 1024香蕉在线观看| 精品一品国产午夜福利视频| 国产精品影院久久| 国产欧美日韩一区二区三区在线| 757午夜福利合集在线观看| 亚洲av片天天在线观看| 日日爽夜夜爽网站| 91大片在线观看| www.精华液| 精品午夜福利视频在线观看一区| 后天国语完整版免费观看| 大型av网站在线播放| 中亚洲国语对白在线视频| av网站在线播放免费| 一本大道久久a久久精品| av天堂在线播放| 一二三四社区在线视频社区8| 在线永久观看黄色视频| 亚洲第一av免费看| 亚洲欧美色中文字幕在线| 亚洲成人免费av在线播放| 美女高潮到喷水免费观看| 欧美精品一区二区免费开放| 69av精品久久久久久| 免费在线观看亚洲国产| 国产精品久久视频播放| 亚洲一区高清亚洲精品| 香蕉久久夜色| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲熟妇中文字幕五十中出 | 俄罗斯特黄特色一大片| 超色免费av| 久久草成人影院| 丰满迷人的少妇在线观看| 国产高清国产精品国产三级| 国产成人免费无遮挡视频| 伊人久久大香线蕉亚洲五| 久久国产精品大桥未久av| 免费在线观看亚洲国产| 国产99久久九九免费精品| 国产精品1区2区在线观看. | av视频免费观看在线观看| 叶爱在线成人免费视频播放| 两性夫妻黄色片| 搡老乐熟女国产| 变态另类成人亚洲欧美熟女 | 老司机午夜福利在线观看视频| 亚洲欧美激情综合另类| 国产国语露脸激情在线看| 午夜精品久久久久久毛片777| 亚洲精华国产精华精| e午夜精品久久久久久久| 极品人妻少妇av视频| 18禁黄网站禁片午夜丰满| 视频区图区小说| 在线av久久热| 国产片内射在线| 亚洲一区中文字幕在线| 男人舔女人的私密视频| 高清视频免费观看一区二区| 高清黄色对白视频在线免费看| 精品国产乱码久久久久久男人| 成在线人永久免费视频| 啪啪无遮挡十八禁网站| 午夜精品在线福利| 一进一出抽搐动态| 狂野欧美激情性xxxx| 久久影院123| av国产精品久久久久影院| 在线观看免费高清a一片| 又黄又爽又免费观看的视频| 丝瓜视频免费看黄片| 99在线人妻在线中文字幕 | 亚洲精品国产一区二区精华液| 国产单亲对白刺激| 亚洲欧美日韩高清在线视频| 日日夜夜操网爽| 国产无遮挡羞羞视频在线观看| 美女国产高潮福利片在线看| 亚洲专区国产一区二区| 丰满饥渴人妻一区二区三| 国产一区二区三区视频了| 精品人妻1区二区| 丝袜美腿诱惑在线| 热99re8久久精品国产| 亚洲国产欧美网| 亚洲成人免费av在线播放| 国产伦人伦偷精品视频| 视频区图区小说| 国产精品久久久久成人av| 在线观看一区二区三区激情| 国产精品.久久久| 成年人黄色毛片网站| 777米奇影视久久| 不卡av一区二区三区| 一区二区日韩欧美中文字幕| 中国美女看黄片| 丰满人妻熟妇乱又伦精品不卡| 侵犯人妻中文字幕一二三四区| 他把我摸到了高潮在线观看| 在线观看免费日韩欧美大片| 黑人欧美特级aaaaaa片| 每晚都被弄得嗷嗷叫到高潮| 老司机影院毛片| 91麻豆av在线| 好看av亚洲va欧美ⅴa在| 丁香欧美五月| 天天操日日干夜夜撸| 成年女人毛片免费观看观看9 | 欧美不卡视频在线免费观看 | 欧美一级毛片孕妇| 精品一区二区三区av网在线观看| 成人永久免费在线观看视频| videos熟女内射| 国产av精品麻豆| 久久精品国产99精品国产亚洲性色 | 丰满迷人的少妇在线观看| 国产精品亚洲一级av第二区| 他把我摸到了高潮在线观看| 午夜福利在线观看吧| 搡老乐熟女国产| 亚洲精品在线美女| 亚洲第一青青草原| 一边摸一边抽搐一进一小说 | 午夜福利乱码中文字幕| 一级作爱视频免费观看| 亚洲欧美日韩另类电影网站| 国产精品.久久久| 欧美丝袜亚洲另类 | 日韩有码中文字幕| 久久国产精品人妻蜜桃| 飞空精品影院首页| 国产亚洲欧美在线一区二区| 别揉我奶头~嗯~啊~动态视频| √禁漫天堂资源中文www| 亚洲一区中文字幕在线| 精品人妻在线不人妻| 午夜福利在线观看吧| 国产高清videossex| 丰满饥渴人妻一区二区三| 99国产精品免费福利视频| 制服诱惑二区| 一边摸一边抽搐一进一小说 | 亚洲三区欧美一区| tube8黄色片| 国产亚洲精品久久久久5区| xxx96com| 亚洲精品在线美女| 国产精品av久久久久免费| 国产成人啪精品午夜网站| 妹子高潮喷水视频| 亚洲精品粉嫩美女一区| 日本一区二区免费在线视频| 成人黄色视频免费在线看| 人人澡人人妻人| 日韩免费av在线播放| 国产免费现黄频在线看| 在线av久久热| 欧美不卡视频在线免费观看 | 午夜91福利影院| 女人被躁到高潮嗷嗷叫费观| 国产蜜桃级精品一区二区三区 | 亚洲av成人av| 女人久久www免费人成看片| 在线观看免费午夜福利视频| 久久热在线av| 一进一出抽搐gif免费好疼 | 久久久水蜜桃国产精品网| 精品国产一区二区久久| 超碰成人久久| 久热这里只有精品99| 青草久久国产| 丝袜美足系列| 热re99久久精品国产66热6| 日韩三级视频一区二区三区| 女警被强在线播放| 亚洲午夜精品一区,二区,三区| 一进一出好大好爽视频| 久久香蕉激情| 又黄又爽又免费观看的视频| 99在线人妻在线中文字幕 | 丁香欧美五月| 一区在线观看完整版| av不卡在线播放| 国产一区二区三区视频了| 日韩制服丝袜自拍偷拍| cao死你这个sao货| 日日夜夜操网爽| 亚洲,欧美精品.| 女人被狂操c到高潮| 一级毛片高清免费大全| 一进一出抽搐gif免费好疼 | 国产一卡二卡三卡精品| 国产麻豆69| 国产精品国产高清国产av | 国产成人av激情在线播放| 亚洲av成人不卡在线观看播放网| 男女下面插进去视频免费观看| 欧美大码av| 国产一区有黄有色的免费视频| 大型av网站在线播放| 欧美大码av| 国产日韩欧美亚洲二区| 丝袜人妻中文字幕| 国精品久久久久久国模美| 宅男免费午夜| 精品久久蜜臀av无| 亚洲全国av大片| 桃红色精品国产亚洲av| 高清欧美精品videossex| 99热只有精品国产| 少妇的丰满在线观看| 老司机在亚洲福利影院| 人妻一区二区av| 亚洲欧美激情综合另类| 欧美日韩瑟瑟在线播放| 亚洲,欧美精品.| 亚洲欧美激情综合另类| 老司机亚洲免费影院| 高清视频免费观看一区二区| 中出人妻视频一区二区| 精品电影一区二区在线| 国产无遮挡羞羞视频在线观看| 乱人伦中国视频| 电影成人av| 亚洲av片天天在线观看| 18禁国产床啪视频网站| 国产成人免费观看mmmm| 少妇 在线观看| 51午夜福利影视在线观看| 真人做人爱边吃奶动态| 亚洲精品久久成人aⅴ小说| 狂野欧美激情性xxxx| 大型黄色视频在线免费观看| 久久九九热精品免费| 欧美不卡视频在线免费观看 | 亚洲欧美日韩高清在线视频| 亚洲第一av免费看| 精品一区二区三区四区五区乱码| 国产精品一区二区精品视频观看| 99国产综合亚洲精品| 亚洲在线自拍视频| 国产淫语在线视频| 国产欧美日韩综合在线一区二区| 亚洲精品一卡2卡三卡4卡5卡| 在线天堂中文资源库| 国产aⅴ精品一区二区三区波| 村上凉子中文字幕在线| 亚洲人成电影免费在线| 建设人人有责人人尽责人人享有的| 午夜福利免费观看在线| 91老司机精品| 久久亚洲精品不卡| 18禁裸乳无遮挡动漫免费视频| 欧美成人午夜精品| 国产精品久久久久久人妻精品电影| 99国产精品99久久久久| 他把我摸到了高潮在线观看| 欧美另类亚洲清纯唯美| 国产日韩欧美亚洲二区| 国产欧美日韩一区二区精品| 欧美日韩瑟瑟在线播放| 国产av精品麻豆| 91精品国产国语对白视频| 国产一卡二卡三卡精品| 国产成人欧美在线观看 | а√天堂www在线а√下载 | 国产精品国产av在线观看| 亚洲人成电影免费在线| 一边摸一边抽搐一进一出视频| 亚洲熟妇熟女久久| a级毛片在线看网站| 一区二区三区精品91| 欧美 亚洲 国产 日韩一| 美女高潮喷水抽搐中文字幕| 久久久久久亚洲精品国产蜜桃av| 日韩免费高清中文字幕av| 深夜精品福利| 大码成人一级视频| 好男人电影高清在线观看| 99riav亚洲国产免费| 亚洲va日本ⅴa欧美va伊人久久| 国产一卡二卡三卡精品| 可以免费在线观看a视频的电影网站| 国产亚洲精品久久久久久毛片 | 欧美中文综合在线视频| 亚洲专区字幕在线| 久久精品国产亚洲av香蕉五月 | 亚洲熟女毛片儿| 少妇裸体淫交视频免费看高清 | 成人免费观看视频高清| 国产激情久久老熟女| 国产精品 欧美亚洲| 亚洲色图av天堂| 成人免费观看视频高清| 国产成人啪精品午夜网站| 美女视频免费永久观看网站| 国产精品 国内视频| 亚洲 欧美一区二区三区| 精品一品国产午夜福利视频| 日本vs欧美在线观看视频| 精品少妇一区二区三区视频日本电影| 女人被躁到高潮嗷嗷叫费观| 视频区欧美日本亚洲| 精品国内亚洲2022精品成人 | а√天堂www在线а√下载 | 欧美色视频一区免费| 亚洲av成人不卡在线观看播放网| 女人被狂操c到高潮| 在线观看免费视频网站a站| 久久人妻av系列| 五月开心婷婷网| 超碰成人久久| 欧美日韩成人在线一区二区| 久久人人97超碰香蕉20202| 女人被狂操c到高潮| 国产av一区二区精品久久| 精品少妇久久久久久888优播| 国产欧美亚洲国产| 午夜福利影视在线免费观看| 亚洲国产欧美一区二区综合| 午夜视频精品福利| 久久久久国产一级毛片高清牌| 亚洲精品久久成人aⅴ小说| 国产成人av激情在线播放| 免费高清在线观看日韩| 欧美精品亚洲一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 人人妻,人人澡人人爽秒播| 免费久久久久久久精品成人欧美视频| 久久中文看片网| 男男h啪啪无遮挡| 狠狠婷婷综合久久久久久88av| 岛国在线观看网站| 欧美日韩视频精品一区| 亚洲美女黄片视频| 久久影院123| 人人妻人人澡人人爽人人夜夜| 高清欧美精品videossex| 一边摸一边抽搐一进一小说 | 亚洲成国产人片在线观看| 在线永久观看黄色视频| 99国产精品99久久久久| av网站免费在线观看视频| 身体一侧抽搐| 老司机午夜福利在线观看视频| 久热这里只有精品99| 精品国内亚洲2022精品成人 | av国产精品久久久久影院| 日韩欧美免费精品| 在线观看免费视频日本深夜| x7x7x7水蜜桃| 国产精品国产av在线观看| 99国产精品99久久久久| 天天操日日干夜夜撸| 飞空精品影院首页| 精品午夜福利视频在线观看一区| 宅男免费午夜| 香蕉久久夜色| 亚洲五月婷婷丁香| 村上凉子中文字幕在线| 两个人看的免费小视频| 免费观看人在逋| 成人特级黄色片久久久久久久| 国产极品粉嫩免费观看在线| 在线观看免费视频网站a站| 看免费av毛片| 男女免费视频国产| 欧美 日韩 精品 国产| 777米奇影视久久| 亚洲精品美女久久久久99蜜臀| 欧美日韩国产mv在线观看视频| 男女下面插进去视频免费观看| 午夜两性在线视频| 最新在线观看一区二区三区| 免费在线观看完整版高清| 天天躁日日躁夜夜躁夜夜| a在线观看视频网站| 亚洲少妇的诱惑av| 久久人妻av系列| 欧美日韩亚洲综合一区二区三区_| 国产乱人伦免费视频| 操出白浆在线播放| 人妻 亚洲 视频| 日日夜夜操网爽| 免费人成视频x8x8入口观看| 亚洲全国av大片| 国产精品 欧美亚洲| 午夜福利免费观看在线| 变态另类成人亚洲欧美熟女 | 久久久久久免费高清国产稀缺| 人人妻人人澡人人看| 国产精品久久久av美女十八| 国产男女内射视频| 丰满的人妻完整版| 高清av免费在线| 黄色a级毛片大全视频| 中文字幕人妻丝袜一区二区| 欧美久久黑人一区二区| 高清视频免费观看一区二区| 看片在线看免费视频| 一夜夜www| 水蜜桃什么品种好| 麻豆国产av国片精品| videos熟女内射| www.精华液| 天堂√8在线中文| 一本大道久久a久久精品| 国产精品久久视频播放| 在线观看www视频免费| 大香蕉久久成人网| 国产一区二区三区视频了| 久9热在线精品视频| 99久久人妻综合| e午夜精品久久久久久久| 在线视频色国产色| 久久久久视频综合| 99久久人妻综合| 欧美人与性动交α欧美精品济南到| 欧美国产精品va在线观看不卡| 国产免费现黄频在线看| 91老司机精品| 欧美午夜高清在线| 精品久久久久久久久久免费视频 | 一级a爱视频在线免费观看| 男人的好看免费观看在线视频 | 精品第一国产精品| 在线天堂中文资源库| 天天躁日日躁夜夜躁夜夜| 高清黄色对白视频在线免费看| 亚洲一区中文字幕在线| 99riav亚洲国产免费| 每晚都被弄得嗷嗷叫到高潮| 精品久久久久久电影网| 在线视频色国产色| 淫妇啪啪啪对白视频| 交换朋友夫妻互换小说| 亚洲伊人色综图| 亚洲欧洲精品一区二区精品久久久| av视频免费观看在线观看| 久久国产精品人妻蜜桃| 男女午夜视频在线观看| 国产av精品麻豆| 免费在线观看视频国产中文字幕亚洲| 视频区欧美日本亚洲| 91在线观看av| 69av精品久久久久久| 激情视频va一区二区三区| 人人妻人人澡人人看| 一级a爱片免费观看的视频| 男女高潮啪啪啪动态图| 日韩三级视频一区二区三区| 美国免费a级毛片|