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

    化工園區(qū)多米諾事故風(fēng)險(xiǎn)評(píng)價(jià)與防控技術(shù)綜述

    2017-04-07 10:27:16賈梅生陳國華胡昆
    化工進(jìn)展 2017年4期
    關(guān)鍵詞:多米諾概率效應(yīng)

    賈梅生,陳國華,胡昆

    (華南理工大學(xué)安全科學(xué)與工程研究所,廣東 廣州 510640)

    化工園區(qū)多米諾事故風(fēng)險(xiǎn)評(píng)價(jià)與防控技術(shù)綜述

    賈梅生,陳國華,胡昆

    (華南理工大學(xué)安全科學(xué)與工程研究所,廣東 廣州 510640)

    定量風(fēng)險(xiǎn)評(píng)價(jià)是化工廠區(qū)或園區(qū)預(yù)防多米諾效應(yīng)的主要方法,各類事故預(yù)防與控制措施的制定及其有效性檢驗(yàn)也依賴于多米諾效應(yīng)定量風(fēng)險(xiǎn)評(píng)價(jià)的結(jié)果。然而,由于概率風(fēng)險(xiǎn)表征、多米諾場景概率結(jié)構(gòu)、火災(zāi)與爆炸致?lián)p概率等存在多種不一致或解釋不明確的方法模型,使風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果差異性很大,無法更進(jìn)一步對(duì)各類多米諾事故預(yù)防與控制技術(shù)進(jìn)行評(píng)估。本文首先界定了多米諾效應(yīng)的基本概念和原理,然后對(duì)當(dāng)前文獻(xiàn)中相關(guān)多米諾效應(yīng)定量風(fēng)險(xiǎn)評(píng)價(jià)方法進(jìn)行深入剖析,明確了多米諾場景概率以及火災(zāi)熱輻射、爆炸沖擊波、爆炸碎片致?lián)p概率的計(jì)算方法。另外,還從閾值與安全距離、安全容量及布局優(yōu)化、安全屏障與防火隔熱層、安全與應(yīng)急管理決策支持、安防與脆弱性分析及災(zāi)害衍生事故等方面對(duì)多米諾事故防控技術(shù)進(jìn)行了評(píng)述。最后,指出以目標(biāo)設(shè)備為核心的各類致?lián)p概率計(jì)算方法,即設(shè)備易損性,是后續(xù)研究的關(guān)注焦點(diǎn)。

    多米諾效應(yīng);風(fēng)險(xiǎn)評(píng)價(jià);致?lián)p概率;設(shè)備易損性;安全;化工園區(qū)

    化工園區(qū)是石油和化學(xué)工業(yè)發(fā)展的必然趨勢,然而集中的園區(qū)管理模式在帶來產(chǎn)業(yè)規(guī)模效應(yīng)的同時(shí),也帶來了較高的風(fēng)險(xiǎn)水平。多米諾效應(yīng)問題的涌現(xiàn)即是由危險(xiǎn)化學(xué)品儲(chǔ)存與工藝裝置規(guī)模的擴(kuò)大,以及土地使用成本的提高而導(dǎo)致的,很多歷史事故案例已經(jīng)證明了多米諾效應(yīng)發(fā)生的可能性及其非常嚴(yán)重的事故后果[1-9]。

    英國健康安全委員會(huì)(HSC)下屬重大危險(xiǎn)源咨詢委員會(huì)在1976年、1979年、1984年公布的重大危險(xiǎn)源控制系列報(bào)告中首次明確了多米諾效應(yīng)概念[10],并在1978年、1981年開展的坎維島風(fēng)險(xiǎn)評(píng)價(jià)報(bào)告中進(jìn)行了詳細(xì)考慮[11]。歷史上,1984年發(fā)生的墨西哥城國家石油公司LPG儲(chǔ)運(yùn)站連鎖爆炸事故[11],極大地促進(jìn)了多米諾效應(yīng)在工業(yè)界與學(xué)術(shù)界的受關(guān)注程度,科技文獻(xiàn)的涌現(xiàn)即始于此。1989年美國化工過程安全中心(CCPS)第一版化工過程定量風(fēng)險(xiǎn)評(píng)價(jià)指南[12],以及1999年荷蘭應(yīng)用科學(xué)研究院(TNO)定量風(fēng)險(xiǎn)評(píng)價(jià)紫皮書[13]已經(jīng)對(duì)多米諾效應(yīng)進(jìn)行了相關(guān)考慮。

    多米諾效應(yīng)的官方條款主要見于歐盟委員會(huì)1996年《塞韋索II指令》第8條及其2012年更新版第9條,主要規(guī)定:企業(yè)需向管理部門告知信息,管理部門有辨識(shí)多米諾效應(yīng)的責(zé)任,相關(guān)企業(yè)需相互合作、互通信息,在編制重大事故預(yù)防策略、安全管理系統(tǒng)、安全報(bào)告、內(nèi)部應(yīng)急預(yù)案時(shí),需考慮多米諾效應(yīng)的影響,并告知公眾[14]。2012年,國務(wù)院安委辦(2012)37號(hào)文件“關(guān)于進(jìn)一步加強(qiáng)化工園區(qū)安全管理的指導(dǎo)意見”中也強(qiáng)調(diào):化工園區(qū)需合理布局,預(yù)防連鎖事故發(fā)生[15]。

    1 多米諾效應(yīng)定義與原理

    眾多多米諾效應(yīng)的概念研究中,RENIERS和COZZANI的定義更具廣泛適用性,其核心是事故擴(kuò)展傳播與后果影響擴(kuò)大,包含3個(gè)基本要素:①初始事故場景及其物理影響,如火災(zāi)熱輻射、爆炸沖擊波、爆炸碎片等;②潛在的二次或一階擴(kuò)展事故場景,源于初始事故的擴(kuò)展傳播,危險(xiǎn)化學(xué)品發(fā)生泄漏,使后果影響擴(kuò)大;③后果影響擴(kuò)大的目標(biāo)設(shè)備或單元[16-19]。

    另外,火災(zāi)熱輻射、爆炸沖擊波、爆炸碎片是三類主要的物理影響因素,即致?lián)p因子。火災(zāi)熱輻射可導(dǎo)致目標(biāo)設(shè)備的殼壁或結(jié)構(gòu)材料發(fā)生高溫強(qiáng)度下降,使壓力儲(chǔ)存的液化烴或工業(yè)氣體的容器內(nèi)壓升高,同時(shí)火焰還有點(diǎn)火引燃的危險(xiǎn)。爆炸沖擊波可使目標(biāo)設(shè)備發(fā)生屈服、倒塌、破裂、分解、整體位移、連接管線斷裂、安全裝置功能失效等,其傳播路徑還可產(chǎn)生二次碎片危險(xiǎn)。爆炸碎片可穿透容器、切斷支撐或管線,高溫碎片也有點(diǎn)火或加熱的危險(xiǎn),碎片來源包括容器爆炸破裂產(chǎn)生的小碎片、爆炸分解的零部件等。一般而言,有毒物質(zhì)泄漏不會(huì)直接導(dǎo)致進(jìn)一步?;沸孤U(kuò)散、火災(zāi)、爆炸等場景,但可使現(xiàn)場操作人員中毒,喪失正常工作能力,間接導(dǎo)致多米諾效應(yīng)。關(guān)于連續(xù)或批處理工藝裝置的基本控制系統(tǒng)與安全功能系統(tǒng)的連鎖效應(yīng)導(dǎo)致的事故,一般不屬于多米諾效應(yīng)的研究范疇。

    2 多米諾事故統(tǒng)計(jì)

    事故統(tǒng)計(jì)分析可以明確多米諾場景的某些特征。2000年,KOURNIOTIS等[1]對(duì)207起化工事故統(tǒng)計(jì)發(fā)現(xiàn),多米諾場景頻率與危險(xiǎn)化學(xué)品類型相關(guān),蒸汽烴類較液體燃料更危險(xiǎn),并擬合了多米諾效應(yīng)的社會(huì)風(fēng)險(xiǎn)曲線。2003年,RONZA等[2]統(tǒng)計(jì)了港口區(qū)的675起化工事故,發(fā)現(xiàn)直接從火災(zāi)到爆炸的擴(kuò)展序列有30個(gè),頻率可達(dá)4.44×10–2。2008年,GóMEZ-MARES等[3]統(tǒng)計(jì)了84起與噴射火相關(guān)的化工事故,指出噴射火有50%的可能性可以觸發(fā)多米諾場景。2010年,DARBRA等[4]對(duì)與多米諾效應(yīng)相關(guān)的225起化工事故,按階段、區(qū)域、涉及危險(xiǎn)化學(xué)品、原因、后果、場所、擴(kuò)展序列等多個(gè)方面進(jìn)行了詳細(xì)統(tǒng)計(jì)分析,并繪制了社會(huì)風(fēng)險(xiǎn)曲線與擴(kuò)展序列概率事件樹圖。類似研究還包括ABDOLHAMIDZADEH等[5]、CHEN等[6],HEMMATIAN等[8]的統(tǒng)計(jì)分析。另外,ZHANG 和ZHENG[7]對(duì)2006—2010年期間國內(nèi)發(fā)生的1632起危險(xiǎn)化學(xué)品事故統(tǒng)計(jì)發(fā)現(xiàn),固定場所多米諾事故占13.2%,移動(dòng)運(yùn)輸場所多米諾事故占9%。2015年,HEMMATIAN等[9]又專門對(duì)127起與沸騰液體擴(kuò)展蒸汽云爆炸(BLEVE)相關(guān)的多米諾事故進(jìn)行統(tǒng)計(jì)分析,繪制了擴(kuò)展序列概率事件樹圖,并指出事故延時(shí)從一分鐘到幾小時(shí)不等,需重點(diǎn)關(guān)注應(yīng)急疏散與消防安全距離。

    3 定量風(fēng)險(xiǎn)評(píng)價(jià)方法

    CCPS化工過程定量風(fēng)險(xiǎn)評(píng)價(jià)指出,處理多米諾效應(yīng)有兩種邏輯[12]:①針對(duì)源設(shè)備或單元,失效頻率與事故場景概率不變,多米諾效應(yīng)使事故后果影響擴(kuò)大,技術(shù)方法是擴(kuò)展事件樹的輸出;②針對(duì)目標(biāo)設(shè)備或單元,多米諾效應(yīng)使失效頻率增加,事故后果影響不變,技術(shù)方法是將多米諾效應(yīng)處理為事故樹的一類外部事件輸入。但無論哪種邏輯,均需首先計(jì)算火災(zāi)熱輻射、爆炸沖擊波、爆炸碎片的致?lián)p概率,或根據(jù)安全距離確定哪些目標(biāo)設(shè)備或單元需納入多米諾效應(yīng)的考慮范圍。另外,由于多米諾效應(yīng)的擴(kuò)展場景與擴(kuò)展模式比較復(fù)雜,多米諾場景概率結(jié)構(gòu)有各種不同的計(jì)算方法,即多米諾風(fēng)險(xiǎn)依多米諾場景概率結(jié)構(gòu)不同而有所差異。表1對(duì)文獻(xiàn)中構(gòu)建的各種多米諾效應(yīng)定量風(fēng)險(xiǎn)評(píng)價(jià)方法進(jìn)行了對(duì)比分析。

    3.1 多米諾場景概率結(jié)構(gòu)

    場景概率結(jié)構(gòu)表達(dá)多米諾效應(yīng)定義及其風(fēng)險(xiǎn)計(jì)算方法,多場景、高階擴(kuò)展、擴(kuò)展協(xié)同效應(yīng)、應(yīng)急措施干預(yù)、動(dòng)態(tài)過程是其技術(shù)難點(diǎn)。如表1所示,

    即使在概率風(fēng)險(xiǎn)評(píng)價(jià)的框架內(nèi),且明確增加失效頻率與擴(kuò)大事故后果影響的邏輯,多數(shù)多米諾場景概率結(jié)構(gòu)依然自成體系,如COZZANI等[16,24-25]的分析方法、ABDOLHAMIDZADEH等[26-27]的蒙特卡洛模擬方法、BERNECHEA等[31]的事件樹方法以及KHAKZAD等[32-34]的貝葉斯網(wǎng)絡(luò)方法。導(dǎo)致混亂的原因除相關(guān)技術(shù)難點(diǎn)外,另一個(gè)重要原因是多米諾風(fēng)險(xiǎn)評(píng)價(jià)對(duì)象不明確或多米諾場景不統(tǒng)一。一般化工過程定量風(fēng)險(xiǎn)評(píng)價(jià)選擇最壞事故場景,但是多米諾效應(yīng)由于擴(kuò)展場景與擴(kuò)展模式比較復(fù)雜,各種組合情況非常多。因此,首先需要確定一種原則或方法以統(tǒng)一多米諾場景的選擇問題。

    3.2 火災(zāi)熱輻射致?lián)p概率

    火災(zāi)熱輻射致?lián)p概率主要包括兩種計(jì)算方法:①比例方法[11,20-23,29-30];②Probit模型,即應(yīng)急及時(shí)概率[16,24-28,32-34,41]。比例方法缺乏嚴(yán)格的理論基礎(chǔ),隨機(jī)性表達(dá)不明確。Probit模型表達(dá)應(yīng)急耗時(shí)小于事故延時(shí)的可能性,其中,事故自然延時(shí)計(jì)算是一種確定的保守?cái)M合關(guān)系式,所以主要反映的是應(yīng)急措施的隨機(jī)性。但是對(duì)于指定的化工園區(qū)或企業(yè)及其應(yīng)急預(yù)案,應(yīng)急耗時(shí)的概率分布是不一樣的,即需要調(diào)整Probit模型系數(shù)。此外,保守?cái)M合關(guān)系式的應(yīng)用也有限制范圍,而且事故延時(shí)本質(zhì)上也是隨機(jī)的,因此,如何分析與表達(dá)事故延時(shí)的隨機(jī)性也需要另作研究。

    3.3 爆炸沖擊波致?lián)p概率

    Probit模型是爆炸沖擊波致?lián)p概率的主要計(jì)算方法,其核心是根據(jù)目標(biāo)設(shè)備破壞失效的嚴(yán)重程度定義致?lián)p概率,概率單位與沖擊波峰值超壓相關(guān),主要反映的是目標(biāo)設(shè)備抵抗沖擊波攻擊的能力,方法的精確性與統(tǒng)計(jì)樣本數(shù)量以及統(tǒng)計(jì)樣本之間的一致性有關(guān)。方法建立在遠(yuǎn)場假設(shè)基礎(chǔ)上,忽略了沖擊波與目標(biāo)設(shè)備之間復(fù)雜的交互過程,不適用于近場范圍的破壞失效。KHAN和ABBASI[21-23]的Probit模型是建立在EISENBERG等[42]模型基礎(chǔ)上的,以有效壓力代替峰值壓力,其他參數(shù)不變,使不同類型的目標(biāo)設(shè)備可以通過曳引系數(shù)不同取值進(jìn)行區(qū)別。COZZANI和SALZANO[43-47]首先將目標(biāo)設(shè)備分為常壓儲(chǔ)罐、壓力容器、塔設(shè)備、輔助設(shè)備四類,然后定義各類設(shè)備破壞失效的嚴(yán)重程度,在75個(gè)文獻(xiàn)樣本的基礎(chǔ)上進(jìn)行統(tǒng)計(jì)分析得到對(duì)應(yīng)Probit函數(shù)。ZHANG和JIANG[48]在COZZANI和SALZANO樣本數(shù)據(jù)的基礎(chǔ)上,重新定義目標(biāo)設(shè)備破壞失效的嚴(yán)重程度,得到新的Probit模型。

    3.4 爆炸碎片致?lián)p概率

    爆炸碎片致?lián)p概率的計(jì)算方法相對(duì)比較復(fù)雜,按研究方法可以分為兩類:①基于碎片拋射軌跡方程的蒙特卡洛模擬;②比例方法。按研究對(duì)象可以分為三類:①碎片產(chǎn)生概率或碎片初始參數(shù)的概率分布,主要包括碎片數(shù)量、大小、形狀、初始速度、拋射角等;②碎片擊中概率;③碎片穿透概率。2001年,HAUPTMANNS[49-50]首次通過碎片拋射軌跡方程的蒙特卡洛模擬得到了碎片的擊中概率,之后涌現(xiàn)出很多以此方法為基礎(chǔ)的研究,主要集中于擊中概率定義、隨機(jī)參數(shù)概率分布確定、參數(shù)影響分析等[51-62]。比例方法主要是GUBINELLI等[63]提出的以概率幾何原理為基礎(chǔ)的擊中概率模型,TUGNOLI等[64]通過一個(gè)事故案例對(duì)該方法的有效性進(jìn)行了分析與驗(yàn)證。關(guān)于碎片初始參數(shù)及其概率分布的確定包括GUBINELLI等[65-66]、MéBARKI等[67]、孫東亮等[68]、TUGNOLI等[69]的研究。關(guān)于碎片穿透概率的研究主要是NGUYEN 等[70]提出的基于臨界殘余壁厚的方法,以及陳剛等[71]提出的基于臨界剩余強(qiáng)度系數(shù)的方法,由于兩種方法同樣是以碎片拋射軌跡方程的蒙特卡洛模擬為基礎(chǔ)的,因此,可以直接獲得碎片擊中概率與穿透概率的乘積。

    表1 多米諾效應(yīng)定量風(fēng)險(xiǎn)評(píng)價(jià)方法對(duì)比分析

    3.5 概率風(fēng)險(xiǎn)表征

    多米諾場景概率結(jié)構(gòu)與火災(zāi)熱輻射、爆炸沖擊波、爆炸碎片的致?lián)p概率最終都要集成于概率風(fēng)險(xiǎn)的框架內(nèi),以概率風(fēng)險(xiǎn)的原理進(jìn)行結(jié)果表征。對(duì)于化工過程定量風(fēng)險(xiǎn)評(píng)價(jià)而言,個(gè)人風(fēng)險(xiǎn)與社會(huì)風(fēng)險(xiǎn)是常用的風(fēng)險(xiǎn)指標(biāo)與風(fēng)險(xiǎn)結(jié)果描述方法,個(gè)人風(fēng)險(xiǎn)是事故場景發(fā)生概率與事故場景物理影響下個(gè)體傷亡概率的數(shù)學(xué)乘積,社會(huì)風(fēng)險(xiǎn)計(jì)算依賴于個(gè)人風(fēng)險(xiǎn)結(jié)果與人員密度分布情況[12-13]。個(gè)體傷亡概率一般基于Probit模型進(jìn)行確定[12-13]。COZZANI等[16]提出的多米諾效應(yīng)定量風(fēng)險(xiǎn)評(píng)價(jià)方法即應(yīng)用個(gè)人風(fēng)險(xiǎn)與社會(huì)風(fēng)險(xiǎn)描述風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果。KOURNIOTIS等[1]、DARBRA等[4]、HEMMATIAN等[8]從事故統(tǒng)計(jì)角度繪制了多米諾效應(yīng)的社會(huì)風(fēng)險(xiǎn)曲線,可以作為風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果的準(zhǔn)則曲線。然而,對(duì)于更多的多米諾效應(yīng)定量風(fēng)險(xiǎn)評(píng)價(jià)方法而言,多米諾場景概率即直接表示多米諾效應(yīng)風(fēng)險(xiǎn)。

    4 事故預(yù)防與控制技術(shù)

    目標(biāo)設(shè)備或單元在致?lián)p因子影響下的破壞失效主要與設(shè)備之間的安全距離、初始事故場景致?lián)p因子強(qiáng)度以及目標(biāo)設(shè)備的防護(hù)措施有關(guān),而關(guān)于降低源設(shè)備失效頻率或初始事故場景物理影響的安全措施,如結(jié)構(gòu)完整性管理、火源控制、保護(hù)層分析(LOPA)等,由于其在一般的過程安全管理中均有考慮,所以不列入詳細(xì)討論范圍,但是初始事故的預(yù)防與安全措施也會(huì)顯著降低多米諾風(fēng)險(xiǎn),而且更為重要。

    4.1 閾值與安全距離、安全容量及布局優(yōu)化

    安全距離與安全容量是互補(bǔ)的、可以實(shí)現(xiàn)本質(zhì)安全的多米諾事故防控技術(shù),兩者均與致?lián)p因子的閾值有關(guān)。閾值主要根據(jù)目標(biāo)設(shè)備或單元類型及其破壞失效模式與嚴(yán)重程度進(jìn)行確定,安全距離根據(jù)閾值與源設(shè)備事故場景特征確定,安全容量設(shè)置主要是限制事故場景物理響應(yīng)范圍。安全距離主要在工廠的布局設(shè)計(jì)階段實(shí)現(xiàn),而安全容量可以實(shí)現(xiàn)動(dòng)態(tài)管理,并且兩者均是布局優(yōu)化與風(fēng)險(xiǎn)動(dòng)態(tài)管理的基本參數(shù)。ALILECHE等[19]對(duì)標(biāo)準(zhǔn)與文獻(xiàn)中關(guān)于多米諾效應(yīng)的閾值與安全距離研究進(jìn)行了詳細(xì)評(píng)述。一般情況,通過實(shí)驗(yàn)或事故統(tǒng)計(jì)可以確定閾值,再根據(jù)事故后果物理影響的經(jīng)驗(yàn)?zāi)P陀?jì)算安全距離,但數(shù)據(jù)不一致仍然是常見的。對(duì)于爆炸沖擊波而言:①缺乏對(duì)目標(biāo)設(shè)備詳細(xì)特征的描述;②破壞失效模式及其嚴(yán)重程度定義比較混亂。對(duì)于火災(zāi)熱輻射而言,還容易忽略:①火災(zāi)場景不同致?lián)p機(jī)理,如接觸點(diǎn)火或遠(yuǎn)距離熱輻射;②火災(zāi)持續(xù)時(shí)間。相對(duì)比較全面細(xì)致的研究是COZZANI等[19,72]的推薦值,還從本質(zhì)安全的角度,在事故后果物理影響經(jīng)驗(yàn)?zāi)P偷幕A(chǔ)上,對(duì)指定事故場景的安全距離與安全容量關(guān)系進(jìn)行了詳細(xì)分析[38,73]。SPOELSTRA等[74]對(duì)荷蘭國家標(biāo)準(zhǔn)PGS18、PGS19中LPG儲(chǔ)罐區(qū)安全距離的確定方法進(jìn)行了詳細(xì)描述。

    理論上,安全距離實(shí)現(xiàn)后,擴(kuò)展場景或多米諾效應(yīng)就不會(huì)發(fā)生了,但安全距離的確定仍有許多不確定或限制因素,只是達(dá)到了目前最好或最優(yōu)水平,因此,布局優(yōu)化有必要將多米諾風(fēng)險(xiǎn)作為一個(gè)目標(biāo)函數(shù),核心問題則是目標(biāo)函數(shù)的構(gòu)建及其求解方法。LEE等[75]、SO等[76]通過柯西最速下降法確定n個(gè)爆炸危險(xiǎn)源的合理位置,以使構(gòu)建的多米諾場景概率最低,前者針對(duì)新建廠區(qū),后者針對(duì)擴(kuò)建廠區(qū)。KHAN和AMYOTTE[77-78]構(gòu)建的I2SI指數(shù)、TUGNOLI等構(gòu)建的DHI指數(shù)[36-37]、COZZANI等[38]構(gòu)建的DCP、DCA、UDI、TDI指數(shù),均可作為布局優(yōu)化的目標(biāo)函數(shù)。JUNG等[79]、LóPEZ-MOLINA等[80]、DE LIRA-FLORES等[81]將布局優(yōu)化定義為混合整數(shù)非線性規(guī)劃問題求解,不同之處是多米諾風(fēng)險(xiǎn)目標(biāo)函數(shù)的構(gòu)建與使用。BERNECHEA和ARNALDOS[82]、KHAKZAD和RENIERS[83]則通過對(duì)不同備選布局方案進(jìn)行多米諾風(fēng)險(xiǎn)分析與評(píng)價(jià),轉(zhuǎn)化為多目標(biāo)決策問題。

    4.2 安全屏障與防火隔熱層

    安全屏障或安全裝置、安全附件與LOPA中的保護(hù)層概念近似,是對(duì)相關(guān)硬件、軟件、管理措施功能抽象的集合描述,從事故預(yù)防與風(fēng)險(xiǎn)管理的角度理解,既包括降低事故場景發(fā)生概率的措施,也包括減小事故后果影響的措施。LOPA指南將保護(hù)層分為8個(gè)層次[84],即工藝設(shè)計(jì)、基本過程控制系統(tǒng)、關(guān)鍵報(bào)警和人員干預(yù)、安全儀表功能、物理保護(hù)(釋放設(shè)施,如安全閥、爆破片等)、釋放后物理保護(hù)(如防火堤、防爆墻等)、廠區(qū)應(yīng)急響應(yīng)、周圍社區(qū)應(yīng)急響應(yīng)。按作用方式與功能,它們又可分為主動(dòng)的與被動(dòng)的措施,以及阻止性的(釋放前)與減緩性的(釋放后)措施。LOPA指南并沒有提供方法以評(píng)估保護(hù)層在降低多米諾風(fēng)險(xiǎn)方面的作用,或指出如何合理布置保護(hù)層以降低多米諾風(fēng)險(xiǎn)。LANDUCCI等[85]基于LOPA方法,得到了預(yù)防火災(zāi)觸發(fā)多米諾場景時(shí),各類安全屏障措施的有效性與要求時(shí)失效概率數(shù)據(jù)。JANSSENS等[86]通過元啟發(fā)式算法,構(gòu)建了一個(gè)合理布置安全屏障以減小火災(zāi)事故場景后果影響,增加擴(kuò)展事故場景延時(shí),降低多米諾風(fēng)險(xiǎn)的決策模型。

    防火隔熱層的概念與方法早于多米諾效應(yīng),主要指阻隔火焰直接接觸影響或降低殼壁吸收熱量速率以增加擴(kuò)展事故場景延時(shí)的涂層或夾套,是一種被動(dòng)的防火措施。BLEVE是一種后果極其嚴(yán)重的火災(zāi)導(dǎo)致的液化烴壓力容器擴(kuò)展事故場景,尤其在運(yùn)輸行業(yè),由于經(jīng)常不可避免地經(jīng)過一些人員密集場所,風(fēng)險(xiǎn)非常高,因此,早期有許多關(guān)于防火隔熱層的設(shè)計(jì)與實(shí)驗(yàn)研究,但是其應(yīng)用也伴隨著腐蝕檢測與消防射流沖擊的問題[11]。LANDUCCI等[87-88]通過對(duì)LPG儲(chǔ)罐的實(shí)驗(yàn)與有限元模擬,分析驗(yàn)證了防火隔熱層在增加BLEVE事故場景延時(shí)方面的有效性,實(shí)驗(yàn)材料包括環(huán)氧樹脂發(fā)泡材料、蛭石、礦物纖維棉、水泥基無機(jī)材料,并提出3個(gè)指標(biāo)以評(píng)估和指導(dǎo)設(shè)計(jì)。PALTRINIERI等[89]進(jìn)一步指出,防火隔熱層的應(yīng)用可以使LPG供應(yīng)鏈的個(gè)人風(fēng)險(xiǎn)降低一個(gè)數(shù)量級(jí),使社會(huì)風(fēng)險(xiǎn)的期望值提高50%。DI PADOVA等[90]、TUGNOLI等[91]基于風(fēng)險(xiǎn)矩陣原理提出一個(gè)防火分區(qū)辨識(shí)方法,以合理應(yīng)用防火隔熱層,降低火災(zāi)觸發(fā)多米諾效應(yīng)風(fēng)險(xiǎn)。此外,由于防火隔熱層也具有一定的結(jié)構(gòu)強(qiáng)度,因此,可以探索研究新的過程設(shè)備設(shè)計(jì)方法,或者開發(fā)設(shè)計(jì)新的保護(hù)層材料或結(jié)構(gòu),以綜合考慮火災(zāi)熱輻射、爆炸沖擊波、爆炸碎片的致?lián)p作用。

    4.3 安全與應(yīng)急管理決策支持

    風(fēng)險(xiǎn)分析與評(píng)價(jià)本質(zhì)上即是一種分級(jí)的安全管理思想與方法,多米諾效應(yīng)的風(fēng)險(xiǎn)分析與評(píng)價(jià)也不例外。應(yīng)急概念則是一種事故控制思想,即假設(shè)事故場景發(fā)生后如何采取措施以減小事故后果影響,多米諾效應(yīng)概念也包含初始事故場景的假設(shè),所以應(yīng)急研究對(duì)于預(yù)防與控制多米諾事故非常重要。KHAN和ABBASI[92]在多米諾效應(yīng)定量風(fēng)險(xiǎn)評(píng)價(jià)方法研究的基礎(chǔ)上,開發(fā)了DOMIFECT軟件,以輔助決策。COZZANI等[93]將多米諾效應(yīng)定量風(fēng)險(xiǎn)評(píng)價(jià)方法集成于Aripar-GIS軟件。RENIERS等[94-95]提出了集成危險(xiǎn)與可操作性分析(HAZOP)、What-If分析、風(fēng)險(xiǎn)矩陣的Hazwim框架,使化工園區(qū)的企業(yè)之間可以相互合作,共同預(yù)防外部多米諾效應(yīng)。之后,RENIERS和DULLAERT[96-97]從化工園區(qū)范圍內(nèi)擴(kuò)展序列優(yōu)先級(jí)的角度開發(fā)了DomPrevPlanning決策支持軟件。另外,RENIERS等[98-99]還應(yīng)用博弈論方法建模分析了化工園區(qū)企業(yè)在預(yù)防多米諾效應(yīng)時(shí)的投資決策過程。ZHANG和CHEN[100]提出了限制化工園區(qū)或企業(yè)危險(xiǎn)源擴(kuò)展網(wǎng)絡(luò)的離散孤島模型。KHAKZAD等[33]應(yīng)用貝葉斯網(wǎng)絡(luò)與沖突分析方法建模了?;钒踩萘康呐渲脝栴}。關(guān)于應(yīng)急管理,ZHOU等[35,101]分別應(yīng)用事件序列圖法與Petri網(wǎng)絡(luò)建模分析了應(yīng)急措施干預(yù)對(duì)火災(zāi)觸發(fā)多米諾場景的影響,以及大規(guī)模并發(fā)火災(zāi)的應(yīng)急響應(yīng)問題。

    4.4 安防與脆弱性分析及災(zāi)害衍生事故

    觸發(fā)化工園區(qū)或企業(yè)多米諾效應(yīng)的原因,除常見的火災(zāi)、爆炸初始事故場景外,蓄意或恐怖襲擊以及自然災(zāi)害也受到廣泛關(guān)注。RENIERS等[102]通過抽象化工園區(qū)范圍內(nèi)所有危險(xiǎn)源及其相互影響關(guān)系,首次提出了分離節(jié)點(diǎn)網(wǎng)絡(luò)的概念,以輔助安防管理。之后,RENIERS和AUDENAERT[103]在分離節(jié)點(diǎn)網(wǎng)理絡(luò)研究的基礎(chǔ)上,提出單元裝置脆弱性分析方法,以確定預(yù)防恐怖襲擊時(shí)單元裝置的防護(hù)優(yōu)先級(jí)。RENIERS等[104]還應(yīng)用物理效應(yīng)衰減與恢復(fù)能力的概念對(duì)化工園區(qū)與企業(yè)的安防管理進(jìn)行了研究,指出分離節(jié)點(diǎn)網(wǎng)絡(luò)的方法可以提高化工園區(qū)遭受攻擊后的快速恢復(fù)能力。LANDUCCI等[105]基于TNT當(dāng)量法研究分析了自制爆炸裝置觸發(fā)多米諾場景時(shí)單元設(shè)施的脆弱性問題。關(guān)于災(zāi)害衍生事故(NaTech),COZZANI等[24]在多米諾效應(yīng)定量風(fēng)險(xiǎn)評(píng)價(jià)方法的基礎(chǔ)上,構(gòu)建了NaTech定量風(fēng)險(xiǎn)評(píng)價(jià)方法。

    5 研究展望

    (1)由于多米諾效應(yīng)概念還有許多爭議,包括范圍界定、場景定義等,使相關(guān)政策、法規(guī)、標(biāo)準(zhǔn)的制定以及定量風(fēng)險(xiǎn)評(píng)價(jià)方法的構(gòu)建存在很多不一致或不相容的解釋與結(jié)果,所以迫切需要在多米諾效應(yīng)概念范圍以及場景選擇方面達(dá)成一種共識(shí),以更好地促進(jìn)工業(yè)實(shí)踐。

    (2)火災(zāi)熱輻射、爆炸沖擊波、爆炸碎片致?lián)p概率的計(jì)算方法與結(jié)果構(gòu)成多米諾效應(yīng)定量風(fēng)險(xiǎn)評(píng)價(jià)框架的基石,當(dāng)前選擇的比例方法或Probit函數(shù),是由于缺乏過程設(shè)備易損性基礎(chǔ)數(shù)據(jù),為簡化定量風(fēng)險(xiǎn)評(píng)價(jià)過程復(fù)雜性與運(yùn)算量而進(jìn)行的折中,因此,需要探索研究新的、高可靠度、高精確度的致?lián)p概率計(jì)算方法,即設(shè)備易損性分析方法,并且在過程設(shè)備設(shè)計(jì)的交付文檔中提供易損性數(shù)據(jù)及其詳細(xì)說明。

    (3)在過程設(shè)備設(shè)計(jì)階段進(jìn)行的致?lián)p概率計(jì)算或易損性分析,可以提高多米諾效應(yīng)定量風(fēng)險(xiǎn)評(píng)價(jià)的有效性,使布局優(yōu)化、安全屏障布置、安全與應(yīng)急管理決策支持更可靠,也可為過程設(shè)備的本質(zhì)安全設(shè)計(jì)及評(píng)估提供理論與方法支持,以更好地從本質(zhì)安全的角度預(yù)防與控制多米諾效應(yīng)。

    (4)可以探索研究新的過程設(shè)備設(shè)計(jì)方法,或開發(fā)新的保護(hù)層材料或結(jié)構(gòu),以綜合考慮火災(zāi)熱輻射、爆炸沖擊波、爆炸碎片的致?lián)p作用。

    [1]KOURNIOTIS S P,KIRANOUDIS C T,MARKATOS N C.Statistical analysis of Domino chemical accidents[J].Journal of Hazardous Materials,2000,71(1/2/3):239-252.

    [2]RONZA A,F(xiàn)éLEZ S,DARBRA R M,et al.Predicting the frequency of accidents in port areas by developing event trees from historical analysis[J].Journal of Loss Prevention in the Process Industries,2003,16(6):551-560.

    [3]GóMEZ-MARES M,ZáRATE L,CASAL J.Jet fires and the Domino effect[J].Fire Safety Journal,2008,43(8):583-588.

    [4]DARBRA R M,PALACIOS A,CASAL J.Domino effect in chemical accidents:main features and accident sequences[J].Journal of Hazardous Materials,2010,183(1/2/3):565-573.

    [5]ABDOLHAMIDZADEH B,ABBASI T,RASHTCHIAN D,et al.Domino effect in process-industry accidents ——an inventory of past events and identification of some patterns[J].Journal of Loss Prevention in the Process Industries,2011,24(5):575-593.

    [6]CHEN Y T,ZHANG M G,GUO P J,et al.2012 International symposium on safety science and technology:investigation and analysis of historical Domino effects statistic[J].Procedia Engineering,2012:45,152-158.

    [7]ZHANG H D,ZHENG X P.Characteristics of hazardous chemical accidents in China:a statistical investigation[J].Journal of Loss Prevention in the Process Industries,2012,25(4):686-693.

    [8]HEMMATIAN B,ABDOLHAMIDZADEH B,DARBRA R M,et al.The significance of Domino effect in chemical accidents[J].Journal of Loss Prevention in the Process Industries,2014,29:30-38.

    [9]HEMMATIAN B,PLANAS E,CASAL J.Fire as a primary event of accident Domino sequences:the case of BLEVE[J].Reliability Engineering & System Safety,2015,139:141-148.

    [10]Advisory Committee on Major Hazards.The control of major hazards - first/ second/ third report[R].London:Healthy and Safety Commission,1976/1979/1984.

    [11]MANNAN S.Lees’loss prevention in the process industries:hazard identification,assessment and control[M].3rd ed.Amsterdam:Elsevier,2005.

    [12]Center for Chemical Process Safety.Guidelines for chemical process quantitative risk analysis[M].2nd ed.New York:American Institute of Chemical Engineers,2000.

    [13]UIJT DE HAAG P A M,ALE B J M.Guideline for quantitative risk assessment ——‘purple book’[M].The Hague:Advisory Council on Dangerous Substances,2005.

    [14]Seveso - III(Directive 2012/18/EU)[EB/OL].European Parliament and Council,2012.[2016-09-05]. http://eur-lex.europa.eu/ legalcontent/EN/TXT/?uri=CELEX:32012L0018.

    [15]國務(wù)院安委會(huì)辦公室.關(guān)于進(jìn)一步加強(qiáng)化工園區(qū)安全管理的指導(dǎo)意見[EB/OL]. 2012. [2016-09-05]. http://www.chinasafety. gov.cn/newpage/Contents/Channel_4976/2012/0809/174959/conten t_174959.htm.The Work Safety Committee Office of The State Council.Guidance on strengthening the safety management of chemical parks[EB/OL].2012.[2016-09-05].http://www.chinasafety.gov.c n/newpage/Contents/Channel_4976/2012/0809/174959/content_1749 59.htm.

    [16]COZZANI V,GUBINELLI G,ANTONIONI G,et al.The assessment of risk caused by Domino effect in quantitative area risk analysis[J].Journal of Hazardous Materials,2005,127(1/2/3):14-30.

    [17]RENIERS G.An external Domino effects investment approach to improve cross-plant safety within chemical clusters[J].Journal of Hazardous Materials,2010,177(1/2/3):167-174.

    [18]RENIERS G,COZZANI V.Domino effects in the process industries——modelling,prevention and managing[M].Amsterdam:Elsevier,2013.

    [19]ALILECHE N,COZZANI V,RENIERS G,et al.Thresholds for Domino effects and safety distances in the process industry:a review of approaches and regulations[J].Reliability Engineering & System Safety,2015,143(s1):74-84.

    [20]BAGSTER D F,PITBLADO R M.The estimation of Domino incident frequencies-an approach[J]. Process Safety and Environmental Protection,1991,69(4):195-199.

    [21]KHAN F I,ABBASI S A.Models for Domino effect analysis in chemical process industries[J].Process Safety Progress,1998,17(2):107-123.

    [22]KHAN F I,ABBASI S A.Studies on the probabilities and likely impacts of chains of accident(Domino effect)in a fertilizer industry[J].Process Safety Progress,2000,19(1):40-56.

    [23]KHAN F I,ABBASI S A.An assessment of the likelihood of occurrence,and the damage potential of Domino effect(chain of accidents)in a typical cluster of industries[J].Journal of Loss Prevention in the Process Industries,2001,14(4):283-306.

    [24]COZZANI V,ANTONIONI G,LANDUCCI G,et al.Quantitative assessment of Domino and NaTech scenarios in complex industrial areas[J].Journal of Loss Prevention in the Process Industries,2014,28(s1):10-22.

    [25]ANTONIONI G,SPADONI G,COZZANI V.Application of Domino effect quantitative risk assessment to an extended industrial area[J].Journal of Loss Prevention in the Process Industries,2009,22(5):614-624.

    [26]ABDOLHAMIDZADEH B,ABBASI T,RASHTCHIAN D,et al.A new method for assessing Domino effect in chemical process industry[J].Journal of Hazardous Materials,2010,182(1/2/3):416-426.

    [27]RAD A,ABDOLHAMIDZADEH B,ABBASI T,et al.FREEDOM II:an improved methodology to assess Domino effect frequency using simulation techniques[J].Process Safety and Environmental Protection,2014,92(6):714-722.

    [28]KADRI F,CHATELET E,CHEN G P.Method for quantitative assessment of the Domino effect in industrial sites[J].Process Safety and Environmental Protection,2013,91(6):452-462.

    [29]張新梅.化工園區(qū)事故多米諾效應(yīng)風(fēng)險(xiǎn)仿真原理及應(yīng)用研究[D].廣州:華南理工大學(xué),2009.ZHANG X M.Study on simulation principle of accident risk caused by Domino effect and its application in chemical industrial parks[D].Guangzhou:South China University of Technology,2009.

    [30]ZHANG X M,CHEN C.Mechanism analysis and risk assessment ofescalation scenario in chemical industry zones[J].Process Safety and Environmental Protection,2013,91(1/2):79-85.

    [31]BERNECHEA E J,ANTONIO VILCHEZ J,ARNALDOS J.A model for estimating the impact of the Domino effect on accident frequencies in quantitative risk assessments of storage facilities[J].Process Safety and Environmental Protection,2013,91(6):423-437.

    [32]KHAKZAD N,KHAN F,AMYOTTE P,et al.Domino effect analysis using Bayesian networks[J].Risk Analysis,2013,33(2):292-306.

    [33]KHAKZAD N,KHAN F,AMYOTTE P,et al.Risk management of Domino effects considering dynamic consequence analysis[J].Risk Analysis,2014,34(6):1128-1138.

    [34]KHAKZAD N.Application of dynamic Bayesian network to risk analysis of Domino effects in chemical infrastructures[J].Reliability Engineering & System Safety,2015,138:263-272.

    [35]ZHOU J F,RENIERS G,KHAKZAD N.Application of event sequence diagram to evaluate emergency response actions during fire-induced Domino effects[J].Reliability Engineering & System Safety,2016,150:202-209.

    [36]TUGNOLI A,KHAN F,AMYOTTE P,et al.Safety assessment in plant layout design using indexing approach:implementing inherent safety perspective:Part 1 guideword applicability and method description[J].Journal of Hazardous Materials,2008,160(1):100-109.

    [37]TUGNOLI A,KHAN F,AMYOTTE P,et al.Safety assessment in plant layout design using indexing approach:implementing inherent safety perspective:Part 2 - Domino hazard index and case study[J].Journal of Hazardous Materials,2008,160(1):110-121.

    [38]COZZANI V,TUGNOLI A,SALZANO E.The development of an inherent safety approach to the prevention of Domino accidents[J].Accident Analysis & Prevention,2009,41(6):1216-1227.

    [39]KHAKZAD N,RENIERS G.Using graph theory to analyze the vulnerability of process plants in the context of cascading effects[J].Reliability Engineering & System Safety,2015,143(s1):63-73.

    [40]NI Z J,WANG Y Z,YIN Z G.Relative risk model for assessing Domino effect in chemical process industry[J].Safety Science,2016,87:156-166.

    [41]LANDUCCI G,GUBINELLI G,ANTONIONI G,et al.The assessment of the damage probability of storage tanks in Domino events triggered by fire[J].Accident Analysis & Prevention,2009,41(6):1206-1215.

    [42]EISENBERG N A,LYNCH C J,BREEDING R J.Vulnerability model:a simulation system for assessing damage resulting from marine spills[R].Springfield:National Technical Information Service,1975.

    [43]COZZANI V,SALZANO E.The quantitative assessment of Domino effects caused by overpressure:Part I.probit models[J].Journal of Hazardous Materials,2004,107(3):67-80.

    [44]COZZANI V,SALZANO E.The quantitative assessment of Domino effect caused by overpressure:Part II.case studies[J].Journal of Hazardous Materials,2004,107(3):81-94.

    [45]COZZANI V,SALZANO E.Threshold values for Domino effects caused by blast wave interaction with process equipment[J].Journal of Loss Prevention in the Process Industries,2004,17(6):437-447.

    [46]SALZANO E,COZZANI V.The analysis of Domino accidents triggered by vapor cloud explosions[J].Reliability Engineering & System Safety,2005,90(2/3):271-284.

    [47]SALZANO E,COZZANI V.A fuzzy set analysis to estimate loss intensity following blast wave interaction with process equipment[J].Journal of Loss Prevention in the Process Industries,2006,19(4):343-352.

    [48]ZHANG M G,JIANG J C.An improved Probit method for assessment of Domino effect to chemical process equipment caused by overpressure[J].Journal of Hazardous Materials,2008,158(2/3):280-286.

    [49]HAUPTMANNS U.A procedure for analyzing the flight of missiles from explosions of cylindrical vessels[J].Journal of Loss Prevention in the Process Industries,2001,14(5):395-402.

    [50]HAUPTMANNS U.A Monte-Carlo based procedure for treating the flight of missiles from tank explosions[J].Probabilistic Engineering Mechanics,2001,16(4):307-312.

    [51]PULA R,KHAN F I,VEITCH B,et al.A model for estimating the probability of missile impact:missiles originating from bursting horizontal cylindrical vessels[J].Process Safety Progress,2007,26(2):129-139.

    [52]張新梅,陳國華.化工罐區(qū)爆炸碎片多米諾效應(yīng)影響概率計(jì)算模型[J].化工學(xué)報(bào),2008,59(11):2946-2953.ZHANG X M,CHEN G H.Impact probability of Domino effect caused by explosion fragments in chemical tank area[J].CIESC Journal,2008,59(11):2946-2953.

    [53]ZHANG X M,CHEN G H.The analysis of Domino effect impact probability triggered by fragments[J].Safety Science,2009,47(7):1026-1032.

    [54]MéBARKI A,NGUYEN Q B,MERCIER F.Structural fragments and explosions in industrial facilities:Part II projectile trajectory and probability of impact[J].Journal of Loss Prevention in the Process Industries,2009,22(4):417-425.

    [55]錢新明,徐亞博,劉振翼.球罐BLEVE碎片拋射的Monte-Carlo分析[J].化工學(xué)報(bào),2009,60(4):1057-1061.QIAN X M,XU Y B,LIU Z Y.Monte-Carlo analysis of fragments projectile from spherical tank BLEVE[J].CIESC Journal,2009,60(4):1057-1061.

    [56]錢新明,徐亞博,劉振翼.球罐BLEVE碎片拋射的危害性研究[J].高壓物理學(xué)報(bào),2009,23(5):389-394.QIAN X M,XU Y B,LIU Z Y.Research on hazard of fragments from spherical tank BLEVE[J].Chinese Journal of High Pressure Physics,2009,23(5):389-394.

    [57]陳剛,朱霽平,武軍,等.爆炸球罐尺寸對(duì)拋射碎片擊中相鄰罐體概率的影響[J].化工學(xué)報(bào),2010,61(6):1599-1604.CHEN G,ZHU W P,WU J,et al.Influence of size of bursting spherical tank on impact probability between fragments and adjacent vessels[J].CIESC Journal,2010,61(6):1599-1604.

    [58]陳剛,朱霽平,武軍,等.化工儲(chǔ)罐間距和體積對(duì)爆炸碎片多米諾效應(yīng)概率的影響[J].火災(zāi)科學(xué),2011,20(1):37-42.CHEN G,ZHU W P,WU J,et al.Impact of gap and volume of storage tanks on Domino effect of explosion fragments[J].Fire Safety Science,2011,20(1): 37-42.

    [59]SUN D L,JIANG J C,ZHANG M G,et al.Parametric approach of the Domino effect for structural fragments[J].Journal of LossPrevention in the Process Industries,2012,25(1):114-126.

    [60]SUN D L,JIANG J C,ZHANG M G,et al.Influence of the source size on Domino effect risk caused by fragments[J].Journal of Loss Prevention in the Process Industries,2015,35:211-223.

    [61]SUN D L,JIANG J C,ZHANG M G,et al.Investigation of multiple Domino scenarios caused by fragments[J].Journal of Loss Prevention in the Process Industries,2016,40:591-602.

    [62]LISI R,CONSOLO G,MASCHIO G,et al.Estimation of the impact probability in Domino effects due to the projection of fragments[J].Process Safety and Environmental Protection,2015,93:99-110.

    [63]GUBINELLI G,ZANELLI S,COZZANI V.A simplified model for the assessment of the impact probability of fragments[J].Journal of Hazardous Materials,2004,116(3):175-187.

    [64]TUGNOLI A,MILAZZO M F,LANDUCCI G,et al.Assessment of the hazard due to fragment projection:a case study[J].Journal of Loss Prevention in the Process Industries,2014,28(SI):36-46.

    [65]GUBINELLI G,COZZANI V.Assessment of missile hazards:evaluation of the fragment number and drag factors[J].Journal of Hazardous Materials,2009,161(1):439-449.

    [66]GUBINELLI G,COZZANI V.Assessment of missile hazards:identification of reference fragmentation patterns[J].Journal of Hazardous Materials,2009,163(2/3):1008-1018.

    [67]MéBARKI A,MERCIER F,NGUYEN Q B,et al.Structural fragments and explosions in industrial facilities:Part I probabilistic description of the source terms[J].Journal of Loss Prevention in the Process Industries,2009,22(4):408-416.

    [68]孫東亮,蔣軍成,張明廣,等.基于最大熵原理的臥罐爆炸碎片數(shù)量概率分布[J].化工學(xué)報(bào),2011,62(s1):219-224.SUN D L,JIANG J C,ZHANG M G,et al.Probabilistic distribution of fragments number from explosion of horizontal tank based on maximum entropy principle[J].CIESC Journal,2011,62(S1):219-224.

    [69]TUGNOLI A,GUBINELLI G,LANDUCCI G,et al.Assessment of fragment projection hazard:probability distributions for the initial direction of fragments[J].Journal of Hazardous Materials,2014,279:418-427.

    [70]NGUYEN Q B,MéBARKI A,SAADA R A,et al.Integrated probabilistic framework for Domino effect and risk analysis[J].Advances in Engineering Software,2009,40(9):892-901.

    [71]陳剛,朱霽平,武軍,等.儲(chǔ)罐爆炸拋射碎片對(duì)目標(biāo)儲(chǔ)罐的破壞概率研究[J].防災(zāi)減災(zāi)工程學(xué)報(bào),2012,32(2):216-222.CHEN G,ZHU W P,WU J,et al.Damage probability model of target storage tank impacted by explosion fragments[J].Journal of Disaster Prevention and Mitigation Engineering,2012,32(2):216-222.

    [72]COZZANI V,GUBINELLI G,SALZANO E.Escalation thresholds in the assessment of Domino accidental events[J].Journal of Hazardous Materials,2006,129(1/2/3):1-21.

    [73]COZZANI V,TUGNOLI A,SALZANO E.Prevention of Domino effect:from active and passive strategies to inherently safer design[J]. Journal of Hazardous Materials,2007,139(2):209-219.

    [74]SPOELSTRA M,MAHESH S,KOOI E,et al.Domino effects at LPG and propane storage sites in the Netherlands[J].Reliability Engineering & System Safety,2015,143(s1):85-90.

    [75]LEE J Y,LEE J W,KO J W,et al.Optimization for allocating the explosive facilities in order to minimize the Domino effect using nonlinear programming[J].Korean Journal of Chemical Engineering,2005,22(5):649-656.

    [76]SO W,KIM Y H,LEE C J,et al.Optimal layout of additional facilities for minimization of Domino effects based on worst-case scenarios[J].Korean Journal of Chemical Engineering,2011,28(3):656-666.

    [77]KHAN F I,AMYOTTE P R.Integrated inherent safety index(I2SI):a tool for inherent safety evaluation[J].Process Safety Progress,2004,23(2):136-148.

    [78]KHAN F I,AMYOTTE P R.I2SI:a comprehensive quantitative tool for inherent safety and cost evaluation[J].Journal of Loss Prevention in the Process Industries,2005,18(4/5/6):310-326.

    [79]JUNG S H,NG D,DIAZ-OVALLE C,et al.New approach to optimizing the facility sitting and layout for fire and explosion scenarios[J].Industrial & Engineering Chemistry Research,2011,50:3928-3937.

    [80]LOPEZ-MOLINA A,VAZQUEZ-ROMAN R,MANNAN M S,et al.An approach for Domino effect reduction based on optimal layouts[J].Journal of Loss Prevention in the Process Industries,2013,26(s1):887-894.

    [81]DE LIRA-FLORES J,VAZQUEZ-ROMAN R,LOPEZ-MOLINA A,et al.A MINLP approach for layout designs based on the Domino hazard index[J].Journal of Loss Prevention in the Process Industries,2014,30:219-227.

    [82]BERNECHEA E J,ARNALDOS J.Optimizing the design of storage facilities through the application of ISD and QRA[J].Process Safety and Environmental Protection,2014,92(6):598-615.

    [83]KHAKZAD N,RENIERS G.Risk-based design of process plants with regard to Domino effects and land use planning[J].Journal of Hazardous Materials,2015,299:289-297.

    [84]Center for Chemical Process Safety.Layer of protection analysis -simplified process risk assessment[M].New York:American Institute of Chemical Engineers,2001.

    [85]LANDUCCI G,ARGENTI F,TUGNOLI A,et al.Quantitative assessment of safety barrier performance in the prevention of Domino scenarios triggered by fire[J].Reliability Engineering & System Safety,2015,143(SI):30-43.

    [86]JANSSENS J,TALARICO L,RENIERS G,et al.A decision model to allocate protective safety barriers and mitigate Domino effects[J].Reliability Engineering & System Safety,2015,143(s1):44-52.

    [87]LANDUCCI G,MOLAG M,REINDERS J,et al.Experimental and analytical investigation of thermal coating effectiveness for 3m3LPG tanks engulfed by fire[J].Journal of Hazardous Materials,2009,161(2/3):1182-1192.

    [88]LANDUCCI G,MOLAG M,COZZANI V.Modeling the performance of coated LPG tanks engulfed in fires[J].Journal of Hazardous Materials,2009,172(1):447-456.

    [89]PALTRINIERI N,LANDUCCI G,MOLAG M,et al.Risk reduction in road and rail LPG transportation by passive fire protection[J].Journal of Hazardous Materials,2009,167(1/2/3):332-344.

    [90]DI PADOVA A,TUGNOLI A,COZZANI V,et al.Identification of fireproofing zones in Oil&Gas facilities by a risk-based procedure[J].Journal of Hazardous Materials,2011,191(1/2/3):83-93.

    [91]TUGNOLI A,COZZANI V,DI PADOVA A,et al.Mitigation of fire damage and escalation by fireproofing:a risk-based strategy[J].Reliability Engineering & System Safety,2012,105(s1):25-35.

    [92]KHAN F I,ABBASI S A.DOMIFFECT(DOMIno eFFECT):user-friendly software for Domino effect analysis[J].Environmental Modelling & Software,1998,13(2):163-177.

    [93]COZZANI V,ANTONIONI G,SPADONI G.Quantitative assessment of Domino scenarios by a GIS-based software tool[J].Journal of Loss Prevention in the Process Industries,2006,19(5):463-477.

    [94]RENIERS G L L,DULLAERT W,ALE B J M,et al.The use of current risk analysis tools evaluated towards preventing external Domino accidents[J].Journal of Loss Prevention in the Process Industries,2005,18(3):119-126.

    [95]RENIERS G L L,DULLAERT W,ALE B J M,et al.Developing an external Domino accident prevention framework:Hazwim[J].Journal of Loss Prevention in the Process Industries,2005,18(3):127-138.

    [96]RENIERS G L L,DULLAERT W.DomPrevPlanning?:user-friendly software for planning Domino effects prevention[J]. Safety Science,2007,45(10):1060-1081.

    [97]RENIERS G L L,DULLAERT W.Knock-on accident prevention in a chemical cluster[J].Expert Systems with Applications,2008,34(1):42-49.

    [98]RENIERS G,DULLAERT W,SOUDAN K.Domino effects within a chemical cluster:a game-theoretical modeling approach by using Nash-equilibrium[J].Journal of Hazardous Materials,2009,167(1/2/3):289-293.

    [99]RENIERS G,CUYPERS S,PAVLOVA Y.A game-theory based Multi-plant Collaboration Model(MCM)for cross-plant prevention in a chemical cluster[J].Journal of Hazardous Materials,2012,209/210:164-176.

    [100]ZHANG X M,CHEN G H.Modeling and algorithm of Domino effect in chemical industrial parks using discrete isolated island method[J].Safety Science,2011,49(3):463-467.

    [101]ZHOU J F,RENIERS G.Petri-net based simulation analysis for emergency response to multiple simultaneous large-scale fires[J].Journal of Loss Prevention in the Process Industries,2016,40:554-562.

    [102]RENIERS G L L,DULLAERT W,AUDENAERT A,et al.Managing Domino effect-related security of industrial areas[J].Journal of Loss Prevention in the Process Industries,2008,21(3):336-343.

    [103]RENIERS G L L,AUDENAERT A.Preparing for major terrorist attacks against chemical clusters:intelligently planning protection measures w.r.t. Domino effects[J],Process Safety and Environmental Protection,2014,92(6):583-589.

    [104]RENIERS G L L,SOERENSEN K,KHAN F,et al.Resilience of chemical industrial areas through attenuation-based security[J]. Reliability Engineering & System Safety,2014,131:94-101.

    [105]LANDUCCI G,RENIERS G,COZZANI V,et al.Vulnerability of industrial facilities to attacks with improvised explosive devices aimed at triggering Domino scenarios[J].Reliability Engineering & System Safety,2015,143(s1):53-62.

    Review of risk assessment and pre-control of Domino effect in Chemical Industry Park

    JIA Meisheng,CHEN Guohua,HU Kun
    (Institute of Safety Science and Engineering,South China University of Technology,Guangzhou 510640,Guangdong,China)

    Quantitative risk assessment of Domino effects(DMQRA)is the key method to be used to prevent Domino accidents in chemical plants or parks. Also,the various pre-control measures and their availability evaluations are based on the result of DMQRA. However,there are many inconsistent or unclear interpreted modes and methods for presentation of probability risk,probability structure of Domino scenarios and damage probabilities of target equipment exposed to fire or blast such that the obtained DMQRA results are in a great variability. Moreover,it is hard to evaluate the specific pre-control measures. In this paper,firstly,the concept and basic principles of Domino effect are clearly defined. Then through carefully analyzing of the various DMQRA methods in literature,the probability structure of Domino scenario and damage probabilities of target equipment attacked by fire thermal radiation,blast wave overpressure and blast fragments are outlined clearly. In addition,damage thresholds,safety distances,safety inventories,layout optimization,safety barriers,thermal insulations,decision supporting for safety and emergency management,security,resilience and natural hazards triggering for domino effect are reviewed in detail. Finally,it is figured out that the advanced methods for calculating of damage probabilities of target equipment exposed to fire,blast and othertriggering vectors,i.e. equipment vulnerability,will be focused upon in the future.

    Domino effect;risk assessment;damage probability;equipment vulnerability;safety;Chemical Industry Park

    TQ086;X937

    A

    1000–6613(2017)04–1534–10

    10.16085/j.issn.1000-6613.2017.04.050

    2016-09-07;修改稿日期:2016-11-04。

    國家自然科學(xué)基金(21576102)及國家重點(diǎn)研發(fā)計(jì)劃(2016YFC0801500)項(xiàng)目。

    賈梅生(1988—),男,博士研究生。聯(lián)系人:陳國華,教授,博士生導(dǎo)師,從事工業(yè)安全與風(fēng)險(xiǎn)評(píng)價(jià)研究。E-mail:mmghchen @scut.edu.cn。

    猜你喜歡
    多米諾概率效應(yīng)
    第6講 “統(tǒng)計(jì)與概率”復(fù)習(xí)精講
    鈾對(duì)大型溞的急性毒性效應(yīng)
    第6講 “統(tǒng)計(jì)與概率”復(fù)習(xí)精講
    概率與統(tǒng)計(jì)(一)
    概率與統(tǒng)計(jì)(二)
    懶馬效應(yīng)
    以反多米諾02號(hào)——木山
    應(yīng)變效應(yīng)及其應(yīng)用
    用創(chuàng)新聯(lián)接未來——多米諾推出更多產(chǎn)品
    塑料包裝(2015年2期)2015-04-09 03:23:06
    三值絕熱多米諾可逆計(jì)數(shù)器設(shè)計(jì)
    久久久久国产网址| 99热网站在线观看| 欧美+日韩+精品| 国产一级毛片在线| 国产精品久久久久久av不卡| 女人精品久久久久毛片| 国产成人精品在线电影| 亚洲国产精品国产精品| 久久韩国三级中文字幕| 精品久久久久久久久亚洲| 亚洲欧美一区二区三区黑人 | 精品国产露脸久久av麻豆| 日本91视频免费播放| av不卡在线播放| 成年人午夜在线观看视频| 亚洲四区av| 在线免费观看不下载黄p国产| 亚洲国产精品国产精品| 日韩精品免费视频一区二区三区 | 99九九在线精品视频| 国产成人精品婷婷| 有码 亚洲区| 欧美最新免费一区二区三区| 纯流量卡能插随身wifi吗| 秋霞在线观看毛片| 国语对白做爰xxxⅹ性视频网站| 一级a做视频免费观看| 极品人妻少妇av视频| 男女免费视频国产| 欧美日韩综合久久久久久| 日韩视频在线欧美| 韩国高清视频一区二区三区| 国产老妇伦熟女老妇高清| 最后的刺客免费高清国语| 亚洲欧美日韩卡通动漫| 国精品久久久久久国模美| 九色成人免费人妻av| 亚洲怡红院男人天堂| 精品一区二区免费观看| 菩萨蛮人人尽说江南好唐韦庄| 交换朋友夫妻互换小说| 国产成人精品无人区| 永久免费av网站大全| 久久av网站| 久久久久久久亚洲中文字幕| 国产黄色视频一区二区在线观看| 日日摸夜夜添夜夜爱| 国产亚洲欧美精品永久| 午夜久久久在线观看| 亚洲av.av天堂| 91午夜精品亚洲一区二区三区| 国产精品一区二区在线观看99| 另类亚洲欧美激情| 赤兔流量卡办理| 欧美日韩亚洲高清精品| 高清不卡的av网站| 黄片无遮挡物在线观看| 在线观看免费日韩欧美大片 | 亚洲精品乱久久久久久| 18在线观看网站| 成人毛片60女人毛片免费| 久久女婷五月综合色啪小说| 免费大片黄手机在线观看| av在线老鸭窝| 国产精品蜜桃在线观看| 丝袜在线中文字幕| 午夜福利影视在线免费观看| 久久精品久久精品一区二区三区| 搡老乐熟女国产| 欧美日韩视频高清一区二区三区二| 久久99热这里只频精品6学生| 少妇猛男粗大的猛烈进出视频| 青春草亚洲视频在线观看| 一级毛片电影观看| 五月玫瑰六月丁香| 日韩av免费高清视频| 欧美精品国产亚洲| 极品人妻少妇av视频| 日韩大片免费观看网站| 日本午夜av视频| 亚洲精品日韩av片在线观看| av网站免费在线观看视频| 岛国毛片在线播放| 尾随美女入室| 欧美变态另类bdsm刘玥| 尾随美女入室| 另类亚洲欧美激情| 免费少妇av软件| 亚洲第一av免费看| 大又大粗又爽又黄少妇毛片口| av在线app专区| 成年av动漫网址| 国产深夜福利视频在线观看| 在线观看人妻少妇| 亚洲国产精品国产精品| 久久久久久久亚洲中文字幕| 成人综合一区亚洲| 天天影视国产精品| 久久精品国产鲁丝片午夜精品| 香蕉精品网在线| 欧美亚洲日本最大视频资源| 国产精品99久久99久久久不卡 | 伦理电影大哥的女人| av视频免费观看在线观看| 大又大粗又爽又黄少妇毛片口| 美女中出高潮动态图| 人妻 亚洲 视频| 少妇的逼水好多| 欧美激情国产日韩精品一区| 欧美+日韩+精品| 国产成人freesex在线| 国产精品99久久久久久久久| 秋霞伦理黄片| 亚洲中文av在线| 久久毛片免费看一区二区三区| 精品一区二区三卡| 国产精品麻豆人妻色哟哟久久| 国产日韩欧美在线精品| 久久婷婷青草| 麻豆成人av视频| 黑丝袜美女国产一区| 伊人久久精品亚洲午夜| 高清视频免费观看一区二区| 高清视频免费观看一区二区| 免费看光身美女| 在线精品无人区一区二区三| 亚洲性久久影院| 尾随美女入室| 大陆偷拍与自拍| 黄色一级大片看看| 亚洲国产精品999| 午夜福利,免费看| 欧美激情极品国产一区二区三区 | 99久久人妻综合| 国产亚洲精品第一综合不卡 | 伊人久久国产一区二区| 女性被躁到高潮视频| 亚洲精品av麻豆狂野| 久久ye,这里只有精品| 国产高清国产精品国产三级| 亚洲欧洲日产国产| 国产欧美另类精品又又久久亚洲欧美| 色吧在线观看| 久久ye,这里只有精品| 国产成人免费无遮挡视频| 日本wwww免费看| 亚洲综合色惰| 久久午夜综合久久蜜桃| 国语对白做爰xxxⅹ性视频网站| av卡一久久| 日本av手机在线免费观看| 人人妻人人添人人爽欧美一区卜| 欧美性感艳星| 久久ye,这里只有精品| 18禁在线播放成人免费| 日本爱情动作片www.在线观看| 久久青草综合色| 中文字幕av电影在线播放| 国产精品久久久久久av不卡| 一二三四中文在线观看免费高清| 中文字幕av电影在线播放| 亚洲国产最新在线播放| 大香蕉97超碰在线| 简卡轻食公司| 日本欧美国产在线视频| 99热全是精品| 亚洲av不卡在线观看| 国产欧美日韩一区二区三区在线 | 午夜福利在线观看免费完整高清在| 日韩一本色道免费dvd| 中文天堂在线官网| 超色免费av| 精品一区二区免费观看| 久久久国产一区二区| 狂野欧美激情性xxxx在线观看| 91在线精品国自产拍蜜月| 最后的刺客免费高清国语| 丝袜脚勾引网站| 国产69精品久久久久777片| 欧美成人午夜免费资源| 大又大粗又爽又黄少妇毛片口| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 插逼视频在线观看| 男的添女的下面高潮视频| 欧美人与性动交α欧美精品济南到 | 人妻夜夜爽99麻豆av| 婷婷色麻豆天堂久久| 2022亚洲国产成人精品| 成人手机av| 一区二区av电影网| 精品午夜福利在线看| 精品一区在线观看国产| 久久久久人妻精品一区果冻| 国产日韩欧美在线精品| 黄色欧美视频在线观看| 成年女人在线观看亚洲视频| 99re6热这里在线精品视频| 欧美最新免费一区二区三区| 午夜福利网站1000一区二区三区| 亚洲精品自拍成人| 一级毛片黄色毛片免费观看视频| 满18在线观看网站| 2022亚洲国产成人精品| 亚洲美女视频黄频| 18禁裸乳无遮挡动漫免费视频| 午夜久久久在线观看| 啦啦啦中文免费视频观看日本| 国内精品宾馆在线| 最近的中文字幕免费完整| 婷婷成人精品国产| 内地一区二区视频在线| 亚洲国产毛片av蜜桃av| 日韩在线高清观看一区二区三区| 国产高清不卡午夜福利| 成人毛片a级毛片在线播放| 免费高清在线观看视频在线观看| 国产极品粉嫩免费观看在线 | 亚洲欧洲精品一区二区精品久久久 | 国产精品不卡视频一区二区| 国产成人精品福利久久| 国产成人精品无人区| 日本av免费视频播放| 国产精品国产三级国产av玫瑰| 国产精品三级大全| 久久精品久久精品一区二区三区| av电影中文网址| 久久久久国产精品人妻一区二区| 成年av动漫网址| 男人添女人高潮全过程视频| 91精品三级在线观看| 成人亚洲欧美一区二区av| 国产成人精品久久久久久| 国产精品成人在线| www.色视频.com| 夫妻午夜视频| 国产精品久久久久久av不卡| 亚洲精品日本国产第一区| 99热这里只有是精品在线观看| 一个人免费看片子| av视频免费观看在线观看| 国产高清不卡午夜福利| 亚洲av国产av综合av卡| 国产国拍精品亚洲av在线观看| 嫩草影院入口| 国产深夜福利视频在线观看| 熟女电影av网| 日韩中文字幕视频在线看片| 大香蕉97超碰在线| 国产精品麻豆人妻色哟哟久久| 美女脱内裤让男人舔精品视频| 精品久久久久久电影网| 卡戴珊不雅视频在线播放| 高清毛片免费看| 国产精品一区二区三区四区免费观看| 大香蕉久久成人网| 日韩三级伦理在线观看| 国精品久久久久久国模美| 成人二区视频| kizo精华| 亚洲av欧美aⅴ国产| 成人国语在线视频| 国产精品女同一区二区软件| 老司机亚洲免费影院| 自线自在国产av| 在线精品无人区一区二区三| 两个人的视频大全免费| av福利片在线| 色视频在线一区二区三区| 亚洲欧美一区二区三区国产| 黑人巨大精品欧美一区二区蜜桃 | 嘟嘟电影网在线观看| 久久这里有精品视频免费| 欧美最新免费一区二区三区| 久久99精品国语久久久| 国产国拍精品亚洲av在线观看| 大香蕉97超碰在线| 免费黄网站久久成人精品| 亚洲国产精品999| 天天操日日干夜夜撸| www.色视频.com| 精品少妇久久久久久888优播| 亚洲精品第二区| 免费黄频网站在线观看国产| 99国产精品免费福利视频| 国产黄色免费在线视频| 国产精品人妻久久久影院| videossex国产| 成年av动漫网址| 欧美人与善性xxx| 亚洲精品aⅴ在线观看| 少妇丰满av| 国产日韩一区二区三区精品不卡 | 亚洲av不卡在线观看| 一区二区日韩欧美中文字幕 | 少妇人妻精品综合一区二区| 高清在线视频一区二区三区| 国产在线免费精品| 亚洲激情五月婷婷啪啪| 日韩中文字幕视频在线看片| 国产成人a∨麻豆精品| 久久99热6这里只有精品| 亚洲精品日韩在线中文字幕| 国产精品偷伦视频观看了| 九色成人免费人妻av| 夜夜骑夜夜射夜夜干| 亚洲不卡免费看| 熟女av电影| 亚洲精品aⅴ在线观看| 国产一区亚洲一区在线观看| 国产免费一区二区三区四区乱码| 久久久久久久大尺度免费视频| 国产探花极品一区二区| 精品人妻偷拍中文字幕| 成人影院久久| 一二三四中文在线观看免费高清| a级毛片黄视频| 成人二区视频| 伦理电影大哥的女人| 国产伦精品一区二区三区视频9| 欧美日韩国产mv在线观看视频| 亚洲成人手机| 欧美xxxx性猛交bbbb| 国产av一区二区精品久久| 国产精品久久久久久久电影| 国产国语露脸激情在线看| 中文字幕av电影在线播放| 久久精品久久久久久久性| 91精品一卡2卡3卡4卡| 超色免费av| 91久久精品国产一区二区三区| 亚洲人与动物交配视频| 久久毛片免费看一区二区三区| 日韩av在线免费看完整版不卡| 欧美另类一区| 精品一区二区三区视频在线| 免费人成在线观看视频色| 免费观看av网站的网址| 日韩大片免费观看网站| 亚洲精品国产av成人精品| 国产在线一区二区三区精| 夫妻性生交免费视频一级片| 国产伦精品一区二区三区视频9| 插阴视频在线观看视频| 新久久久久国产一级毛片| 亚洲经典国产精华液单| 亚洲伊人久久精品综合| 亚洲高清免费不卡视频| 日日啪夜夜爽| 在线天堂最新版资源| 色哟哟·www| 国产 精品1| 少妇的逼水好多| 亚洲一区二区三区欧美精品| 成人漫画全彩无遮挡| 熟妇人妻不卡中文字幕| 三级国产精品欧美在线观看| 亚洲美女黄色视频免费看| 亚洲欧美日韩卡通动漫| 久久久久久久久久久丰满| 色网站视频免费| 日韩大片免费观看网站| 卡戴珊不雅视频在线播放| 亚洲三级黄色毛片| 一区在线观看完整版| 能在线免费看毛片的网站| av专区在线播放| 国产永久视频网站| 在线观看三级黄色| 美女cb高潮喷水在线观看| 两个人免费观看高清视频| 69精品国产乱码久久久| 午夜老司机福利剧场| 欧美一级a爱片免费观看看| 午夜日本视频在线| 日日爽夜夜爽网站| 国产亚洲精品第一综合不卡 | 伦精品一区二区三区| 最新中文字幕久久久久| 99久国产av精品国产电影| 国产午夜精品一二区理论片| 国产精品成人在线| 男男h啪啪无遮挡| 伊人久久国产一区二区| 成人毛片60女人毛片免费| 欧美另类一区| 午夜精品国产一区二区电影| 夜夜看夜夜爽夜夜摸| 国产精品三级大全| 国产精品久久久久久精品古装| 国产成人精品福利久久| 国产亚洲欧美精品永久| 国产视频内射| 国产白丝娇喘喷水9色精品| 女人精品久久久久毛片| 下体分泌物呈黄色| 看免费成人av毛片| 亚洲av福利一区| 欧美3d第一页| 精品熟女少妇av免费看| 黑人欧美特级aaaaaa片| 日日啪夜夜爽| 永久免费av网站大全| 中国国产av一级| 十八禁高潮呻吟视频| 天天操日日干夜夜撸| 黑人欧美特级aaaaaa片| 亚洲av成人精品一区久久| 亚洲精品久久久久久婷婷小说| 国产成人免费无遮挡视频| 免费少妇av软件| 日本av免费视频播放| 嫩草影院入口| 亚洲丝袜综合中文字幕| 日韩中字成人| 亚洲精品日本国产第一区| 91精品一卡2卡3卡4卡| 亚洲国产精品专区欧美| 亚洲欧美一区二区三区国产| 亚洲中文av在线| 夫妻午夜视频| 午夜精品国产一区二区电影| 男的添女的下面高潮视频| 一二三四中文在线观看免费高清| 亚洲三级黄色毛片| av免费观看日本| 亚洲四区av| 国产在线视频一区二区| 亚洲av福利一区| 亚洲av不卡在线观看| 亚洲第一区二区三区不卡| 黄片播放在线免费| 久久久a久久爽久久v久久| 成人手机av| 国产精品久久久久成人av| 国产精品国产av在线观看| 最新中文字幕久久久久| 亚洲欧美中文字幕日韩二区| 一个人看视频在线观看www免费| 国产国拍精品亚洲av在线观看| 18禁动态无遮挡网站| 国产成人精品无人区| 人人妻人人添人人爽欧美一区卜| 十八禁高潮呻吟视频| 亚洲精品亚洲一区二区| 精品少妇内射三级| 考比视频在线观看| 亚洲精品乱久久久久久| 亚洲精品国产av蜜桃| 一级,二级,三级黄色视频| 亚洲高清免费不卡视频| 热99久久久久精品小说推荐| 亚洲欧洲国产日韩| 国产精品久久久久久精品电影小说| 在线精品无人区一区二区三| 日韩av免费高清视频| 香蕉精品网在线| 22中文网久久字幕| 91久久精品国产一区二区成人| 久久国产精品男人的天堂亚洲 | 黑丝袜美女国产一区| 国产精品一区二区在线不卡| av黄色大香蕉| 母亲3免费完整高清在线观看 | 日韩 亚洲 欧美在线| 成人免费观看视频高清| 自拍欧美九色日韩亚洲蝌蚪91| 蜜桃国产av成人99| 两个人的视频大全免费| av卡一久久| 亚洲无线观看免费| 十分钟在线观看高清视频www| 久久精品久久久久久久性| 熟女电影av网| 欧美日韩视频高清一区二区三区二| av在线老鸭窝| 亚洲天堂av无毛| 黄色视频在线播放观看不卡| 日韩熟女老妇一区二区性免费视频| 亚洲av福利一区| 999精品在线视频| 中文天堂在线官网| 亚洲四区av| 色94色欧美一区二区| 9色porny在线观看| 欧美97在线视频| 中国三级夫妇交换| 国产爽快片一区二区三区| 久久热精品热| 夜夜骑夜夜射夜夜干| 国产毛片在线视频| 多毛熟女@视频| 日本91视频免费播放| 日本av免费视频播放| 欧美xxxx性猛交bbbb| 午夜福利网站1000一区二区三区| 一区二区三区乱码不卡18| 色94色欧美一区二区| 免费大片18禁| 人妻 亚洲 视频| 亚洲精品色激情综合| 久久精品国产鲁丝片午夜精品| 国产熟女午夜一区二区三区 | 在线天堂最新版资源| 国产精品久久久久久av不卡| 免费大片黄手机在线观看| 欧美精品国产亚洲| 国产69精品久久久久777片| 国产成人精品无人区| 中文字幕精品免费在线观看视频 | 大香蕉久久网| 欧美日韩视频高清一区二区三区二| 精品一品国产午夜福利视频| 桃花免费在线播放| 纯流量卡能插随身wifi吗| 久久久午夜欧美精品| 曰老女人黄片| 国产成人精品久久久久久| 成人无遮挡网站| 久久鲁丝午夜福利片| 久久ye,这里只有精品| 久久精品熟女亚洲av麻豆精品| 色吧在线观看| 久久精品国产亚洲av天美| 欧美人与性动交α欧美精品济南到 | 亚洲综合精品二区| 亚洲精品美女久久av网站| 久久午夜综合久久蜜桃| 久热这里只有精品99| 日韩三级伦理在线观看| 99久久人妻综合| av福利片在线| 亚洲国产欧美日韩在线播放| 一区二区三区精品91| 最近的中文字幕免费完整| 久久热精品热| 久久久久精品性色| 日本黄大片高清| 人妻人人澡人人爽人人| 大陆偷拍与自拍| 国产在视频线精品| 满18在线观看网站| 人妻少妇偷人精品九色| 国产成人精品无人区| 国产老妇伦熟女老妇高清| 日本黄大片高清| 亚洲国产色片| 亚洲国产精品一区二区三区在线| 国产精品久久久久成人av| 精品国产露脸久久av麻豆| 一级毛片我不卡| 国产成人精品久久久久久| 热99国产精品久久久久久7| av.在线天堂| 美女cb高潮喷水在线观看| 精品国产露脸久久av麻豆| 亚洲综合精品二区| av.在线天堂| 热99国产精品久久久久久7| 国产精品久久久久久久电影| 91精品伊人久久大香线蕉| 黑丝袜美女国产一区| 日韩免费高清中文字幕av| 高清黄色对白视频在线免费看| 黄片播放在线免费| 国产探花极品一区二区| 人人妻人人添人人爽欧美一区卜| 在线精品无人区一区二区三| 亚洲,一卡二卡三卡| 国产精品一区二区在线观看99| 亚洲高清免费不卡视频| 黑丝袜美女国产一区| 亚洲色图综合在线观看| 亚洲中文av在线| 99九九在线精品视频| 肉色欧美久久久久久久蜜桃| 免费日韩欧美在线观看| 中文天堂在线官网| 少妇 在线观看| 如日韩欧美国产精品一区二区三区 | 人妻人人澡人人爽人人| 久久人妻熟女aⅴ| 免费观看的影片在线观看| 成年美女黄网站色视频大全免费 | 黑丝袜美女国产一区| 99九九在线精品视频| 精品国产露脸久久av麻豆| 美女cb高潮喷水在线观看| 寂寞人妻少妇视频99o| 下体分泌物呈黄色| 亚洲精品自拍成人| 久久精品久久久久久久性| 黄色欧美视频在线观看| 中文精品一卡2卡3卡4更新| 日韩欧美一区视频在线观看| 男女高潮啪啪啪动态图| 亚洲无线观看免费| 亚洲色图综合在线观看| videosex国产| 日韩欧美精品免费久久| 全区人妻精品视频| 国产精品熟女久久久久浪| 亚洲精品自拍成人| 久久精品久久久久久久性| av免费观看日本| 久久久久久久久久久丰满| 一级毛片 在线播放| 成人综合一区亚洲| 精品久久国产蜜桃| 午夜福利,免费看| 卡戴珊不雅视频在线播放| 插逼视频在线观看| 久久久久网色| 亚洲av国产av综合av卡| 欧美3d第一页| 日韩精品有码人妻一区|