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

    深水浮式結(jié)構(gòu)火災(zāi)爆炸風(fēng)險評價方法及軟件開發(fā)

    2017-01-06 13:42:27王亞瓊余建星盧賀帥
    關(guān)鍵詞:熱輻射后果通量

    王亞瓊,余建星,梁 靜,盧賀帥,李 驍

    (1.天津大學(xué)水利工程仿真與安全國家重點實驗室,天津300072;2.高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心,上海200240)

    深水浮式結(jié)構(gòu)火災(zāi)爆炸風(fēng)險評價方法及軟件開發(fā)

    王亞瓊1,2,余建星1,2,梁 靜1,2,盧賀帥1,2,李 驍1,2

    (1.天津大學(xué)水利工程仿真與安全國家重點實驗室,天津300072;2.高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心,上海200240)

    在美國石油學(xué)會基于風(fēng)險的檢驗技術(shù)(APIRP 581 RISK-BASED INSPECTION TECHNOLOGY)規(guī)范的基礎(chǔ)上,結(jié)合我國南海深海石油開發(fā)實際,編寫了一套功能較全、計算模型更加符合中國南??赡芤l(fā)的火災(zāi)爆炸風(fēng)險評價軟件.軟件首先在內(nèi)核上結(jié)合火災(zāi)爆炸指數(shù)法評價,對待評估對象有整體的危險性認(rèn)識;為簡化API 581規(guī)范中對噴射火和池火危害半徑循環(huán)迭代的計算,對公式進行變形,采用MATLAB畫出曲線,導(dǎo)出危害半徑和熱輻射通量直接的關(guān)系,可以得到任意熱輻射通量下對應(yīng)的目標(biāo)距離.其次在概率和風(fēng)險評價設(shè)計中,實現(xiàn)了事件樹和風(fēng)險矩陣的輸出.最后在輸出結(jié)果上結(jié)合ACCESS數(shù)據(jù)庫,實現(xiàn)了存儲和輸出計算數(shù)據(jù),方便后期的分析使用.通過和PHAST計算算例的結(jié)果比較,軟件可靠、有效,并且在概率計算和風(fēng)險分析上具有優(yōu)勢.

    火災(zāi)爆炸;風(fēng)險評價;API581規(guī)范;軟件開發(fā)

    近年來,隨著經(jīng)濟的高速增長,我國對油氣等燃料的需求日益上升[1].我國南海蘊藏著豐富的石油天然氣和其他礦產(chǎn)資源,南海油氣資源是我國經(jīng)濟與社會可持續(xù)發(fā)展的強大支持和保障因素.近年對南海油氣的開采、運輸日益增多.在可燃性油氣的開采以及儲運的過程中,經(jīng)常會發(fā)生油氣泄漏.可燃性油氣在點火的情況下,會發(fā)生各種火災(zāi)爆炸,例如:噴射火,池火,蒸氣云爆炸等[2].火災(zāi)爆炸會對油氣造成損失,對開采、儲運油氣設(shè)備、財產(chǎn)以及人員造成傷害.海洋工程界通常認(rèn)為300 m以上為深水,1 500 m以上為超深水[3].深水海域,如南海群島,距離大陸較遠,一旦發(fā)生油氣泄漏造成火災(zāi)爆炸,撲救比較困難,因而造成的人員、設(shè)備損失更加嚴(yán)重.因此,對深水浮式結(jié)構(gòu)的火災(zāi)爆炸風(fēng)險評估顯得尤為重要.

    隨著計算機技術(shù)的發(fā)展,運用計算機計算油氣的泄漏以及引發(fā)的事故后果成為可能[4].目前,國內(nèi)外風(fēng)險評價、分析的方法非常之多,各種定量分析的軟件也較多[5-6].挪威船級社(DNV)在安全管理與技術(shù)評價領(lǐng)域內(nèi)具有權(quán)威的地位.最常用的火災(zāi)爆炸風(fēng)險評價軟件是挪威DNV公司開發(fā)的PHAST量化風(fēng)險評估計算軟件.PHAST軟件包實施事故后果、風(fēng)險計算,根據(jù)泄漏介質(zhì)的特性、當(dāng)?shù)貧鉁?、風(fēng)速、點火源位置及人員分布等條件,結(jié)合事故發(fā)生的頻率計算得到各種可能的燃燒性、爆炸性和毒性的后果及風(fēng)險[7].

    但是在現(xiàn)有的風(fēng)險評價軟件中,池火災(zāi)的計算模型一般采用的是點火源模型,精確程度不如API581規(guī)范中介紹的Mudan模型,缺乏對管道或者設(shè)備的整體風(fēng)險程度評價,在事故發(fā)生概率的計算和事件樹的輸出方面不完善,用戶在分析火災(zāi)爆炸計算后果的概率和風(fēng)險等級時不直觀.

    為了編寫一套功能較全、計算模型更加合理的火災(zāi)爆炸風(fēng)險評價軟件,本文以API581規(guī)范和《安全評價方法應(yīng)用指南》為理論算法,針對油氣的火災(zāi)爆炸特性,對規(guī)范的計算公式進行變形,加入了《安全評價方法應(yīng)用指南》中的指數(shù)法評價方法,采用VisualBasic.NET開發(fā)了一套火災(zāi)爆炸風(fēng)險評估軟件,對火災(zāi)爆炸的危險等級進行劃分,對各類事故發(fā)生的可能性、后果以及風(fēng)險等級進行計算評估.

    1 算法改進

    API 581規(guī)范中對噴射火和池火的后果計算部分,采取的是公式中參數(shù)循環(huán)迭代的方法,來計算火災(zāi)的影響區(qū)域.計算復(fù)雜,并且計算量大,不直觀.因此,本文對火災(zāi)的后果區(qū)域的計算公式進行變形,用在編寫的程序中,理論算法更加簡單.

    1.1 噴射火后果區(qū)域算法改進

    噴射火主要是由于對受壓的工藝容器或者管道發(fā)生泄漏后的可燃流體進行立即點火產(chǎn)生的.火災(zāi)不同的熱輻射量對人和物可造成不同程度的損害,如表1所示.

    表1 熱輻射傷害準(zhǔn)則[8]Tab.1 Thermal radiation damage criteria

    API中的規(guī)定熱輻射通量的分界值是人員傷害12.6 kW/m2,設(shè)備損壞是37.8 kW/m2[9].因此,最后計算出來的目標(biāo)處的熱輻射通量要逼近于這兩個臨界值.

    首先取目標(biāo)距離的大概值,計算出大氣輻射系數(shù)和視角系數(shù).根據(jù)噴射火的總輻射能量和計算出來的大氣輻射系數(shù)和視角系數(shù)帶入求解目標(biāo)處的熱輻射通量.將計算值和規(guī)定值 (人員傷害12.6 kW/m2,設(shè)備損壞是37.8 kW/m2)進行比較,如果接近,則輸出目標(biāo)距離,如果相差較大,則改變目標(biāo)距離的值,重新帶入求解,循環(huán)迭代直至目標(biāo)處的熱輻射通量接近規(guī)定值.具體公式參照API581規(guī)范.計算人員和設(shè)備的影響距離步驟如圖1所示.

    從上面的計算過程可以看出,反復(fù)循環(huán)迭代計算量增大,最后計算出來的目標(biāo)處的熱輻射通量是在一定的誤差范圍內(nèi)接近于規(guī)定值,誤差范圍的大小也直接影響到最后計算出來的目標(biāo)距離,導(dǎo)致結(jié)果不準(zhǔn)確.

    圖1 噴射火安全距離計算步驟Fig.1 Jet fire safety distance calculation steps

    為了使計算安全距離更加直接、方便,將API581規(guī)范中的公式進行變形,推導(dǎo)出目標(biāo)處接受的熱通量Ithjnet和安全距離xsn的函數(shù)表達式如下:

    其中:C14為API中給定的系數(shù),可以從附表中查出為1; 為燃燒能量的分?jǐn)?shù),一般取0.35;Wnjet為噴射火泄露速率;HCV為燃燒熱,J/kg;Ithjnet為熱通量,W/m2;Pw為大氣層的水分壓,kPa.

    在公式 (1)中,熱通量的值,API中的規(guī)定分界值是人員傷害12.6 kW/m2,設(shè)備損壞是37.8 kW/m2,即可直接求出人員和設(shè)備的安全距離,變形后的公式計算方法比之前的迭代計算要簡便,大大減小計算量.

    1.2 池火災(zāi)后果區(qū)域算法改進

    當(dāng)易燃液體從設(shè)備或者管道泄露出來,就會形成池火.隨著池火的形成,一些液體會蒸發(fā),并且如果蒸發(fā)的易燃液體遇到高于他的可燃下限(LFL)的點火源時,池火就會形成.對于APIRBI,認(rèn)為池火是由發(fā)生泄漏孔口或者破裂的受壓容器或者管道泄露的可燃液體立即點火形成的.

    池火災(zāi)的計算過程和噴射火大體類似.和噴射火不同的是,池火災(zāi)在計算輻射角系數(shù)時用到的是Mudan模型.視角系數(shù)可以把火焰看做一個垂直的圓柱的傾斜角來計算,運用Mudan提供的公式.為圓柱形的視角系數(shù)Fcyln是由圓柱的垂直視角系數(shù)Fvn和圓柱的水平視角系數(shù)Fhn計算得出[4].其中圓柱的垂直視角系數(shù)和水平視角系數(shù)是火焰傾角、火焰長度、火池半徑和目標(biāo)距離的函數(shù).計算公式非常復(fù)雜,具體可參照API581規(guī)范.由于視角系數(shù)的計算過于復(fù)雜,導(dǎo)致用上述噴射火的循環(huán)迭代計算方法的計算量非常大,在迭代過程中,出現(xiàn)錯誤的幾率也加大.

    為了能夠直接計算出在一定熱輻射通量下的目標(biāo)距離,本文采用MATLAB對上述公式進行變形,畫出距離和熱輻射通量的曲線圖.從曲線圖中可以直觀的找到目標(biāo)距離.

    選用可燃液體的代表性物質(zhì)乙醇進行計算,以火池半徑為10m為例,在MATLAB中將上述基本數(shù)據(jù)以及API581中的池火災(zāi)中的計算垂直輻射角系數(shù)、水平輻射角系數(shù)、輻射角系數(shù)、大氣透射系數(shù)和目標(biāo)處接收到的熱輻射通量公式輸入,即可畫出曲線如圖2.其中橫軸為距離泄漏點的距離,縱軸為對應(yīng)距離的熱輻射通量.

    圖2 安全距離-熱輻射通量曲線圖Fig.2 Safety distance radiation heat flux curve

    在圖2中,y軸為對應(yīng)距離的熱輻射通量,從圖像上選取y軸上任一點的熱輻射通量值,對應(yīng)到曲線上的點,可以讀出點的橫坐標(biāo)即為目標(biāo)距離的值.即可以得到每一個熱輻射通量下對應(yīng)的距離.

    從上面的分析過程可以看出[10-12],在計算后果影響距離時,計算不用迭代求解,根據(jù)圖像的橫縱坐標(biāo)可以直接讀出,借助曲線圖的方法直觀、簡便.

    2 軟件開發(fā)

    軟件的編寫主要以API581規(guī)范和《安全評價方法應(yīng)用

    指南》為依據(jù),具體模塊的算法如下.

    軟件的第1個模塊,火災(zāi)爆炸風(fēng)險程度評價模塊,參照《安全評價方法應(yīng)用指南》中的火災(zāi)爆炸危險指數(shù)(F&EI)評價方法.該方法從設(shè)備或管線內(nèi)流體危險性出發(fā),評價工藝過程中的一般工藝危險與特殊工藝危險,得到考慮補償措施之前的火災(zāi)爆炸危險指數(shù),再引入諸多補償措施的折減系數(shù)(工藝控制、物質(zhì)隔離以及防火措施),得到最終的危險指數(shù)[13].

    軟件的第2、3、4模塊,即事件樹概率分析模塊、火災(zāi)爆炸事故預(yù)測模塊、風(fēng)險評價模塊,參照的是由美國石油學(xué)會(API)于2008年發(fā)布的基于風(fēng)險的檢驗技術(shù)(RISK-BASED INSPECTION TECHNOLOGY),即 API 581規(guī)范2008版.噴射火和池火的后果計算部分,采用上面改進過的公式和算法,避開循環(huán)迭代的求解目標(biāo)距離,減小了計算量,用在程序中更加簡單、直觀.

    在進行深水結(jié)構(gòu)火災(zāi)爆炸事故后果預(yù)測之前,先對待評估對象的火災(zāi)爆炸危險性進行系統(tǒng)預(yù)分析,有的放矢地找出重點評估對象,該過程常用到的方法有火災(zāi)爆炸危險指數(shù)(F&EI)評價方法.然后可以根據(jù)泄漏的判定事件樹,計算事件樹對應(yīng)的后果的發(fā)生概率.根據(jù)對系統(tǒng)的危險程度評價以及不同事故的發(fā)生概率大小,有針對性的對火災(zāi)爆炸事故后果進行計算.最后通過對事故后果發(fā)生的概率和危險區(qū)域、人員傷害范圍進行量化計算和分級,從而評定風(fēng)險等級,在風(fēng)險矩陣中輸出風(fēng)險等級.軟件的計算思路如圖3.

    該系統(tǒng)通過計算火災(zāi)爆炸危險指數(shù)和危險等級、事件樹上對應(yīng)事故的發(fā)生概率、火災(zāi)爆炸范圍計算,可以實現(xiàn)風(fēng)險評估,對用戶的實際工程進行指導(dǎo).編寫的這套軟件主要具有以下特點:

    1)軟件內(nèi)核加入《安全評價方法應(yīng)用指南》中的火災(zāi)爆炸風(fēng)險程度評價,進行設(shè)備或管線火災(zāi)爆炸事故后果預(yù)測之前,根據(jù)用戶輸入的物質(zhì)系數(shù)、工藝危險系數(shù)和補償系數(shù)計算系統(tǒng)的危險指數(shù)(F&EI),對待評估對象進行整體的危險程度評定;2)在事故的發(fā)生概率計算模塊,軟件可以輸出事件樹,并對事件樹發(fā)生事故的對應(yīng)概率計算,用戶使用方便,計算結(jié)果一目了然;3)噴射火和池火的后果影響距離的計算,噴射火計算對API中的公式進行變形,池火計算采用API 581中提供的Mudan模型,對火災(zāi)后果范圍的計算更為直接;4)軟件可以通過對事故發(fā)生概率和后果范圍的計算,實現(xiàn)對待評估對象進行風(fēng)險評價,并且在風(fēng)險矩陣中輸出風(fēng)險等級,對于指導(dǎo)用戶的生產(chǎn)實踐有實際意義;5)軟件融合ACCESS數(shù)據(jù)庫,對用戶輸入的參數(shù)和計算的結(jié)果進行儲存,最后可以輸出Word和PDF報告,方便用戶的后期使用.

    圖3 軟件的計算思路Fig.3 Method of calculation software

    3 算例

    本文選取具體的算例,使用編寫的火災(zāi)爆炸風(fēng)險評價軟件和PHAST軟件分別計算分析.并對二者計算結(jié)果的進行對比分析.

    采用發(fā)生池火災(zāi)的常見物質(zhì)乙醇.環(huán)境溫度為20℃,標(biāo)準(zhǔn)大氣壓下,風(fēng)速為4m/s,火池半徑為10m.乙醇在此條件下的物性參數(shù)如表2所示[14].

    表2 乙醇物性參數(shù)表Tab.2 Physicalparametersofethanol table

    首先對事故的整體危險性進行火災(zāi)爆炸風(fēng)險程度評價和后果概率的計算,結(jié)果如圖4、圖5所示.圖4是對系統(tǒng)的整體危險性進行評價;圖5是輸出事件樹,并且在事件樹對應(yīng)后果后輸出計算的后果概率,用戶可以清楚的看出發(fā)生池火災(zāi)的概率為0.34,是所有后果中最大的,因此需要對池火災(zāi)進行后果計算.

    圖4 火災(zāi)爆炸危險程度評價Fig.4 Fireexplosion hazard degreeevaluation

    下面對池火災(zāi)的后果進行具體的計算,API中的規(guī)定熱輻射通量的分界值是人員傷害12.6 kW/m2,設(shè)備損壞是37.8 kW/m2.池火災(zāi)的相關(guān)后果參數(shù)如表3所示.

    表3 池火災(zāi)后果計算值Tab.3 Pool fire consequence calculation

    最后進行軟件的風(fēng)險評價模塊計算,將前面計算的后果概率和后果半徑分別評級,計算出風(fēng)險等級.界面如圖6,其中將風(fēng)險等級在風(fēng)險矩陣上輸出,用戶使用更加直觀.

    為了和PHAST軟件進行對比,分別使用本文編寫的軟件和PHAST軟件計算池火災(zāi)的后果.池火災(zāi)的相關(guān)后果參數(shù)對比如表4、表5所示.

    表4 池火災(zāi)后果計算值對比Tab.4 Comparison of calculated resultsof pool fire

    表5 安全距離對比Tab.5 Comparison of safety distance

    圖5 概率計算和事件樹的輸出Fig.5 Theoutputof probability calculation and event tree

    圖6 風(fēng)險矩陣計算界面Fig.6 Riskmatrix interface

    從表4、表5中的后果參數(shù)對比可以看出,PHAST計算后果參數(shù)值的比較,二者的計算結(jié)果基本相近本文編寫的火災(zāi)爆炸風(fēng)險評價軟件具有合理性、有效性;此外,本軟件計算出來的傷害半徑和PHAST相比稍大,計算結(jié)果偏保守,對實際工程運用更加安全.

    4 結(jié)語

    1)變形后的噴射火和池火安全距離計算公式簡便,避開了循環(huán)迭代的求解目標(biāo)距離,減小了計算量,可以觀察出熱輻射通量和距離之間的關(guān)系,直接得到任意熱輻射通量對應(yīng)的安全距離值;

    2)軟件在API581規(guī)范的理論算法上,結(jié)合《安全評價方法應(yīng)用指南》的火災(zāi)爆炸指數(shù)法評價,對評價模型多了一層在后果評估之前的危險性評價,對待評估對象有整體的危險性認(rèn)識;軟件在概率和風(fēng)險評價設(shè)計中,使用事件樹和風(fēng)險矩陣的輸出,用戶在使用時對后果認(rèn)識一目了然;軟件可以對用戶輸入的參數(shù)和計算的結(jié)果進行儲存,最后輸出Word和PDF報告,方便用戶的后期分析使用;

    3)通過和PHAST輸出結(jié)果的對比,PHAST偏重于后果值的計算,而本軟件的優(yōu)勢在于能計算后果發(fā)生的概率和風(fēng)險等級,并且在輸出上采取了事件樹和風(fēng)險矩陣,輸出界面清晰,利于用戶的直觀理解和使用,軟件更加人性化,功能更加全面.

    [1]李耀莊,楊志,姜學(xué)鵬.駁船池火災(zāi)的危險性分析 [J].災(zāi)害學(xué),2007,22(4):100-103.

    [2]宋雪飛,陸建輝.海洋平臺池火火災(zāi)模型概述 [J].石油工程建設(shè),2011,37(2):1-4.

    [3]王麗勤,侯金林,龐然,等.深水油氣田開發(fā)工程中的基礎(chǔ)應(yīng)用探討 [J].海洋石油,2011,31(4):88-92.

    [4]Xing Zhixiang,Zhou Yan,DaiChuang.The research and software developmenton simulation assessmentmethod for oilgas storage tank's fire& explosion[J].Computersand Applied Chem istry,2010,27(8):1055-1058.

    [5]陳國華,梁韜,張暉,等.用SAFETI定量評價液氯泄漏事故的風(fēng)險 [J].華南理工大學(xué)學(xué)報(自然科學(xué)版),2006,34(5):103-108.

    [6]楊源,余建星,杜尊峰,等.基于RBFNN的FDPSO火災(zāi)爆炸波及鉆井區(qū)概率分析 [J].天津大學(xué)學(xué)報,2012,45(10):917-923.

    [7]謝忠海.儲備庫的事故后果定量風(fēng)險分析 [J].油氣田地面工程,2010,29(4):75-77.

    [8]蔡賓斌.油池火熱輻射危害模型計算及對比分析研究 [J].中國應(yīng)急救援,2013(2):47-50.

    [9]APIRP 581RISK-BASED INSPECTION TECHNOLOGY[S].

    [10]Zeng Xiaohui,Yu Yang,Zhang Liang,et al.A New Energy-Absorbing Device for Motion Suppression in Deep-Sea Floating Platforms[J].Energies,2015,8(1):111-132.

    [11]Liu Q Q,Shu A P,Singh V P.Closure to Analysis of the vertical profile of concentration in sediment-laden flows[J].Journal of Engineering Mechanics-ASCE,2009,135(7):744-745.

    [12]LiX S,Yang B,Duan LP,etal.Experimentalstudy ongasproduction frommethanehydrate in porousmediaby SAGDmethod[J].ApplEnerg.2013,112(0):1233-40.

    [13]劉鐵民,張興凱,劉功智,等.安全評價方法應(yīng)用指南 [M].北京:化學(xué)工業(yè)出版社,2005.

    [14]胡超,朱國慶,吳維華,等.池火危害模型化計算分析研究 [J].消防科學(xué)與技術(shù),2011,30(7):570-573.

    [責(zé)任編輯 楊 屹]

    The fireand explosion risk assessmentmethod of the deepwater floating structureand the software development

    WANG Yaqiong1,2,YU Jianxing1,2,LIANG Jing1,2,LU Heshuai1,2,LIXiao1,2

    (1.StateKey Laboratory ofHydraulic Engineering Simulationand Safety,Tianjin University,Tianjin300072,China;2.Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration,Shanghai200240,China)

    According to the American Petroleum Institute released of RISK-BASED INSPECTION TECHNOLOGY in 2008,combinedwith theactualoildevelopmentof the South China Sea,asetofwith comprehensive functionsandmore reasonablemodels in linewith the South China Sea fireexplosion risk assessmentsoftwarewasdeveloped.First,the software combined the assessmentof fire and explosion indexmethod to evaluate the risk of object.In order to simplify the iterativecalculation of the jetfireand pool firehazard radiusin API581 specification,formulaswasdeformated and curves was drawn using MATLAB.The function of hazard radius and radiation heat flux was derived.The targetdistance corresponding to any thermal radiation flux can be gotten.Second,in the design of the probability and risk,the software outputevent trees and riskmatrix.Last,combined with ACCESS database in the output,the software for storage and output of data analysis,was convenient for later use.Compared with the PHAST calculation results,this software is reliable and effective,and hasadvantages in calculating the probability and risk analysis.

    fire and explosion;risk assessment;API581 specification;software development

    TE58

    A

    1007-2373(2016)01-0096-06

    10.14081/j.cnki.hgdxb.2016.01.018

    2015-05-25

    國家重點基礎(chǔ)研究發(fā)展計劃(973計劃)(2014CB046804);國家自然科學(xué)基金(51239008,51379145);“十二五”國家科技重大專項子課題(2011ZX05030-006-03)

    王亞瓊(1989-),女(漢族),碩士生.

    梁靜(1988-),女(漢族),博士生,liangjing_0805@126.com.

    猜你喜歡
    熱輻射后果通量
    天津大學(xué)的熱輻射催化乙烷脫氫制乙烯研究獲進展
    冬小麥田N2O通量研究
    “耍帥”的后果
    快樂語文(2021年31期)2022-01-18 05:51:32
    這些行為后果很嚴(yán)重
    熱輻射的危害
    水上消防(2020年5期)2020-12-14 07:16:26
    眾榮的后果8則
    不同水系統(tǒng)阻隔熱輻射研究進展
    上朝遲到了 后果很嚴(yán)重
    學(xué)與玩(2017年3期)2017-02-16 07:05:00
    緩釋型固體二氧化氯的制備及其釋放通量的影響因素
    化工進展(2015年6期)2015-11-13 00:26:29
    春、夏季長江口及鄰近海域溶解甲烷的分布與釋放通量
    精品国内亚洲2022精品成人| 亚洲av电影在线进入| 国产aⅴ精品一区二区三区波| 国产伦人伦偷精品视频| 国产单亲对白刺激| 黑人巨大精品欧美一区二区mp4| 国产精品久久久久久亚洲av鲁大| 99热这里只有是精品50| 麻豆av在线久日| 高潮久久久久久久久久久不卡| 国产黄色小视频在线观看| 麻豆成人午夜福利视频| 丝袜美腿诱惑在线| 国内揄拍国产精品人妻在线| 日本一区二区免费在线视频| 久久这里只有精品中国| 日韩中文字幕欧美一区二区| 1024香蕉在线观看| 亚洲色图av天堂| 日日夜夜操网爽| 看片在线看免费视频| 成人18禁高潮啪啪吃奶动态图| 最近最新免费中文字幕在线| 亚洲一码二码三码区别大吗| 午夜福利在线在线| 久久久久国内视频| 亚洲,欧美精品.| 观看免费一级毛片| videosex国产| 天堂√8在线中文| 久久精品成人免费网站| 久久精品夜夜夜夜夜久久蜜豆 | 99国产综合亚洲精品| 草草在线视频免费看| ponron亚洲| 国产99白浆流出| 全区人妻精品视频| 精品福利观看| 手机成人av网站| 欧美绝顶高潮抽搐喷水| 国产成人精品久久二区二区免费| 日日夜夜操网爽| 亚洲av成人精品一区久久| 在线播放国产精品三级| 色尼玛亚洲综合影院| 国产激情久久老熟女| 成熟少妇高潮喷水视频| 亚洲国产日韩欧美精品在线观看 | 成人国产综合亚洲| 中文资源天堂在线| 窝窝影院91人妻| 亚洲 欧美一区二区三区| 欧美久久黑人一区二区| 欧美国产日韩亚洲一区| 日本三级黄在线观看| 99久久国产精品久久久| 亚洲真实伦在线观看| 久久久久亚洲av毛片大全| 我的老师免费观看完整版| 女人被狂操c到高潮| 国产单亲对白刺激| 法律面前人人平等表现在哪些方面| 欧美国产日韩亚洲一区| 国产亚洲精品av在线| 午夜两性在线视频| 午夜影院日韩av| 亚洲美女视频黄频| 久久精品国产亚洲av香蕉五月| 两个人免费观看高清视频| 久久99热这里只有精品18| 变态另类丝袜制服| av有码第一页| 麻豆久久精品国产亚洲av| 精品久久久久久,| 99久久精品热视频| 视频区欧美日本亚洲| 精品少妇一区二区三区视频日本电影| 国产高清视频在线播放一区| 日韩欧美国产一区二区入口| 丰满人妻一区二区三区视频av | 欧美av亚洲av综合av国产av| 无遮挡黄片免费观看| 天天躁夜夜躁狠狠躁躁| 999久久久国产精品视频| 18禁观看日本| 亚洲真实伦在线观看| 成人国产一区最新在线观看| 欧美av亚洲av综合av国产av| 久久精品aⅴ一区二区三区四区| 亚洲成人中文字幕在线播放| 国产三级在线视频| 波多野结衣高清无吗| 欧美乱色亚洲激情| 在线播放国产精品三级| 欧美 亚洲 国产 日韩一| 亚洲全国av大片| 亚洲精品在线观看二区| 美女大奶头视频| 1024香蕉在线观看| 舔av片在线| 在线观看午夜福利视频| 久久草成人影院| 中文字幕av在线有码专区| 亚洲欧洲精品一区二区精品久久久| 中文字幕av在线有码专区| 无限看片的www在线观看| 妹子高潮喷水视频| 亚洲专区中文字幕在线| 久久久久亚洲av毛片大全| 男女下面进入的视频免费午夜| 在线免费观看的www视频| 国产黄色小视频在线观看| 香蕉国产在线看| 国产69精品久久久久777片 | 欧美日韩一级在线毛片| 又大又爽又粗| 国产精品1区2区在线观看.| 久久久久久免费高清国产稀缺| 久久精品国产综合久久久| 国产成人精品久久二区二区免费| 可以在线观看毛片的网站| 免费观看精品视频网站| 亚洲一区二区三区不卡视频| 国产1区2区3区精品| 欧美日韩精品网址| 成人特级黄色片久久久久久久| 日本黄大片高清| 成人高潮视频无遮挡免费网站| 波多野结衣高清作品| 91麻豆精品激情在线观看国产| 1024香蕉在线观看| 精品人妻1区二区| 亚洲欧美一区二区三区黑人| 在线观看免费日韩欧美大片| 午夜免费激情av| 黄频高清免费视频| 亚洲国产精品合色在线| 91麻豆av在线| 国产又黄又爽又无遮挡在线| 黄色丝袜av网址大全| 黄色毛片三级朝国网站| 国产精品98久久久久久宅男小说| 一区二区三区高清视频在线| 精品久久久久久,| 亚洲国产精品999在线| 1024香蕉在线观看| 99国产精品一区二区蜜桃av| 麻豆国产97在线/欧美 | 舔av片在线| 亚洲欧美日韩高清专用| 亚洲成a人片在线一区二区| 三级毛片av免费| 日韩欧美一区二区三区在线观看| 欧美黑人欧美精品刺激| 国产精品精品国产色婷婷| 欧美乱妇无乱码| 国产v大片淫在线免费观看| 黄色毛片三级朝国网站| 欧美精品亚洲一区二区| 少妇人妻一区二区三区视频| 国产三级黄色录像| 夜夜爽天天搞| 国产成人aa在线观看| 精品日产1卡2卡| 亚洲色图 男人天堂 中文字幕| 久久久久久久久久黄片| 国产91精品成人一区二区三区| 色噜噜av男人的天堂激情| 亚洲av日韩精品久久久久久密| 午夜福利视频1000在线观看| 桃色一区二区三区在线观看| 国产精品精品国产色婷婷| 午夜福利在线观看吧| 国产成人欧美在线观看| 女人爽到高潮嗷嗷叫在线视频| 手机成人av网站| 国产精品影院久久| 19禁男女啪啪无遮挡网站| 免费看十八禁软件| 91麻豆av在线| 波多野结衣巨乳人妻| 亚洲精品中文字幕在线视频| 日韩精品免费视频一区二区三区| 最近视频中文字幕2019在线8| 桃红色精品国产亚洲av| 久久伊人香网站| 亚洲专区字幕在线| 国产三级在线视频| 视频区欧美日本亚洲| 亚洲色图 男人天堂 中文字幕| 久久国产精品人妻蜜桃| 亚洲色图av天堂| 伊人久久大香线蕉亚洲五| 少妇被粗大的猛进出69影院| 黄频高清免费视频| 久久久久久亚洲精品国产蜜桃av| 亚洲av中文字字幕乱码综合| 妹子高潮喷水视频| 国产精品一及| 最近视频中文字幕2019在线8| 久9热在线精品视频| 亚洲第一电影网av| 熟妇人妻久久中文字幕3abv| 国产爱豆传媒在线观看 | 久久精品成人免费网站| 精品久久蜜臀av无| 超碰成人久久| 亚洲美女黄片视频| 欧美日韩中文字幕国产精品一区二区三区| 夜夜爽天天搞| 国产亚洲欧美98| 丁香六月欧美| 999久久久国产精品视频| 久久这里只有精品中国| 久久久水蜜桃国产精品网| 国产私拍福利视频在线观看| 成人午夜高清在线视频| 亚洲欧美激情综合另类| 精品久久久久久久久久久久久| 窝窝影院91人妻| 久久精品成人免费网站| 国产成人av教育| 熟女少妇亚洲综合色aaa.| 天堂av国产一区二区熟女人妻 | 免费在线观看亚洲国产| 最近在线观看免费完整版| 亚洲成a人片在线一区二区| 精品第一国产精品| 午夜激情av网站| 2021天堂中文幕一二区在线观| 亚洲中文av在线| 18禁裸乳无遮挡免费网站照片| 校园春色视频在线观看| 久久久久久久久中文| 国产精品一区二区精品视频观看| 97碰自拍视频| 18禁国产床啪视频网站| 久久久精品欧美日韩精品| 毛片女人毛片| 99热这里只有精品一区 | 午夜精品在线福利| 特大巨黑吊av在线直播| 人妻丰满熟妇av一区二区三区| 在线观看舔阴道视频| 亚洲人成网站在线播放欧美日韩| 女人高潮潮喷娇喘18禁视频| 亚洲熟妇中文字幕五十中出| 国产精品 国内视频| 亚洲成a人片在线一区二区| 国产成人av激情在线播放| 免费无遮挡裸体视频| 婷婷精品国产亚洲av在线| 美女高潮喷水抽搐中文字幕| 亚洲av熟女| 成年免费大片在线观看| 中文字幕久久专区| 国产精品免费视频内射| 一级a爱片免费观看的视频| 琪琪午夜伦伦电影理论片6080| 在线观看免费日韩欧美大片| 亚洲国产精品久久男人天堂| 一本大道久久a久久精品| 国产av在哪里看| 国产av不卡久久| 这个男人来自地球电影免费观看| e午夜精品久久久久久久| 日日夜夜操网爽| 亚洲熟女毛片儿| 精品人妻1区二区| 特级一级黄色大片| 老汉色∧v一级毛片| 亚洲成人精品中文字幕电影| 在线十欧美十亚洲十日本专区| 久久久久久九九精品二区国产 | 男插女下体视频免费在线播放| 91麻豆av在线| 中文字幕av在线有码专区| 一级黄色大片毛片| 午夜福利欧美成人| 免费观看人在逋| 久久久久久亚洲精品国产蜜桃av| 久久久久免费精品人妻一区二区| 国产午夜福利久久久久久| 久久这里只有精品中国| 久久久国产成人精品二区| 国产av麻豆久久久久久久| 三级国产精品欧美在线观看 | 亚洲国产看品久久| 视频区欧美日本亚洲| tocl精华| 啦啦啦观看免费观看视频高清| 三级男女做爰猛烈吃奶摸视频| 精华霜和精华液先用哪个| 听说在线观看完整版免费高清| 欧美成狂野欧美在线观看| 久久伊人香网站| 国产精品一区二区精品视频观看| 国产成人系列免费观看| 婷婷六月久久综合丁香| 99久久无色码亚洲精品果冻| 狂野欧美激情性xxxx| 国产亚洲精品一区二区www| 大型黄色视频在线免费观看| 日本精品一区二区三区蜜桃| 这个男人来自地球电影免费观看| 国产片内射在线| 黄色成人免费大全| 久久精品aⅴ一区二区三区四区| 国产成人精品无人区| 少妇的丰满在线观看| 好看av亚洲va欧美ⅴa在| 此物有八面人人有两片| 欧美性猛交黑人性爽| 久99久视频精品免费| 日日摸夜夜添夜夜添小说| 亚洲精品中文字幕在线视频| 亚洲av美国av| 日韩高清综合在线| 久久伊人香网站| 在线观看66精品国产| 在线永久观看黄色视频| 久久精品91无色码中文字幕| 亚洲第一电影网av| 亚洲中文字幕日韩| 免费在线观看完整版高清| 丝袜人妻中文字幕| 欧美国产日韩亚洲一区| 欧美成人午夜精品| 1024香蕉在线观看| 国产精品 国内视频| 丰满人妻熟妇乱又伦精品不卡| 久久香蕉精品热| 久久精品影院6| 在线看三级毛片| 国产免费男女视频| 亚洲国产看品久久| 一级片免费观看大全| 最近最新中文字幕大全电影3| 亚洲精品一卡2卡三卡4卡5卡| 在线免费观看的www视频| 禁无遮挡网站| 亚洲精品av麻豆狂野| 国产爱豆传媒在线观看 | 久久久久久久久免费视频了| 悠悠久久av| 性色av乱码一区二区三区2| 欧美中文日本在线观看视频| 国产亚洲精品av在线| 亚洲国产看品久久| 99riav亚洲国产免费| 精品高清国产在线一区| 亚洲精品一区av在线观看| 亚洲九九香蕉| 色综合欧美亚洲国产小说| 欧美日韩精品网址| 亚洲男人天堂网一区| 成年人黄色毛片网站| 一级毛片精品| 亚洲真实伦在线观看| 国产成人精品久久二区二区91| 久久伊人香网站| 精品国产超薄肉色丝袜足j| 黄色视频不卡| bbb黄色大片| 一本大道久久a久久精品| 午夜a级毛片| a在线观看视频网站| 人妻丰满熟妇av一区二区三区| 国产精品久久久久久亚洲av鲁大| 制服丝袜大香蕉在线| 午夜福利成人在线免费观看| 777久久人妻少妇嫩草av网站| 在线免费观看的www视频| 国产成人av教育| 国产成+人综合+亚洲专区| 国产免费av片在线观看野外av| 桃红色精品国产亚洲av| 国产1区2区3区精品| 国产精品久久久av美女十八| 一级片免费观看大全| 免费观看精品视频网站| 精品电影一区二区在线| 欧美性猛交╳xxx乱大交人| 国产单亲对白刺激| 久久国产精品影院| 国产黄a三级三级三级人| 最新在线观看一区二区三区| 国产黄片美女视频| 国产精品一及| 国产熟女午夜一区二区三区| 精品久久久久久久人妻蜜臀av| 亚洲欧美精品综合一区二区三区| 午夜免费观看网址| 1024手机看黄色片| 中文资源天堂在线| 国产主播在线观看一区二区| tocl精华| 日韩国内少妇激情av| 精品久久久久久久久久免费视频| 又紧又爽又黄一区二区| 国产高清视频在线观看网站| www.999成人在线观看| 午夜免费观看网址| 国产欧美日韩一区二区精品| 天堂影院成人在线观看| 中文字幕最新亚洲高清| 婷婷精品国产亚洲av在线| 91av网站免费观看| 日本一本二区三区精品| 欧美日本亚洲视频在线播放| 麻豆av在线久日| www.熟女人妻精品国产| 老熟妇仑乱视频hdxx| 女生性感内裤真人,穿戴方法视频| 俺也久久电影网| 黄色视频不卡| 亚洲欧美激情综合另类| 国产aⅴ精品一区二区三区波| 欧美不卡视频在线免费观看 | 激情在线观看视频在线高清| 日韩精品中文字幕看吧| 欧美午夜高清在线| 成年人黄色毛片网站| 国产探花在线观看一区二区| 999久久久国产精品视频| 无遮挡黄片免费观看| 国产一区二区三区在线臀色熟女| 欧美日韩福利视频一区二区| 少妇熟女aⅴ在线视频| 国产免费av片在线观看野外av| 女警被强在线播放| 成人永久免费在线观看视频| 亚洲成人国产一区在线观看| 亚洲欧美日韩高清专用| 国产免费男女视频| 国产真人三级小视频在线观看| 两个人视频免费观看高清| 久久精品夜夜夜夜夜久久蜜豆 | 国产一区二区三区在线臀色熟女| 国产精品久久电影中文字幕| 变态另类成人亚洲欧美熟女| 日本 av在线| 国产成人影院久久av| 久久久久精品国产欧美久久久| 久久精品人妻少妇| 日韩欧美免费精品| 高潮久久久久久久久久久不卡| 国产视频一区二区在线看| 久久香蕉国产精品| 国产69精品久久久久777片 | 亚洲av片天天在线观看| 久久精品国产清高在天天线| 亚洲欧美日韩高清专用| 欧美精品啪啪一区二区三区| 他把我摸到了高潮在线观看| 我的老师免费观看完整版| 亚洲免费av在线视频| 亚洲欧美精品综合一区二区三区| 欧美性猛交╳xxx乱大交人| 成熟少妇高潮喷水视频| 午夜福利高清视频| 老汉色av国产亚洲站长工具| 国产久久久一区二区三区| 999久久久精品免费观看国产| 久久久久久久午夜电影| 首页视频小说图片口味搜索| 亚洲欧美日韩高清在线视频| 日本免费一区二区三区高清不卡| 久久久久久大精品| 国产aⅴ精品一区二区三区波| 欧美一区二区精品小视频在线| 最近最新免费中文字幕在线| 久久热在线av| 久久精品人妻少妇| tocl精华| 精品不卡国产一区二区三区| 国产爱豆传媒在线观看 | 成人av一区二区三区在线看| 啪啪无遮挡十八禁网站| 亚洲人成网站高清观看| 亚洲一区二区三区不卡视频| 欧美激情久久久久久爽电影| 天堂√8在线中文| 成人手机av| 色综合婷婷激情| 天堂√8在线中文| 午夜成年电影在线免费观看| 亚洲成人久久性| 69av精品久久久久久| 亚洲欧美日韩高清在线视频| 搡老岳熟女国产| 久久九九热精品免费| 婷婷精品国产亚洲av| 99在线人妻在线中文字幕| 国产高清videossex| 大型黄色视频在线免费观看| 午夜福利视频1000在线观看| 在线观看舔阴道视频| 麻豆国产97在线/欧美 | 午夜福利高清视频| 宅男免费午夜| 日韩欧美国产一区二区入口| 久久久久九九精品影院| 亚洲在线自拍视频| 亚洲人成电影免费在线| 国产伦一二天堂av在线观看| 国产三级黄色录像| 欧洲精品卡2卡3卡4卡5卡区| 村上凉子中文字幕在线| 91国产中文字幕| 激情在线观看视频在线高清| 久久精品国产清高在天天线| 成人亚洲精品av一区二区| 麻豆久久精品国产亚洲av| 国产成人影院久久av| 亚洲成av人片在线播放无| 欧美黑人欧美精品刺激| 亚洲一区高清亚洲精品| 亚洲人成网站高清观看| 国产又色又爽无遮挡免费看| 免费看日本二区| 欧美av亚洲av综合av国产av| 国产区一区二久久| 亚洲精品一区av在线观看| 黄色片一级片一级黄色片| 亚洲精品中文字幕在线视频| 男插女下体视频免费在线播放| 国语自产精品视频在线第100页| 免费人成视频x8x8入口观看| 韩国av一区二区三区四区| 19禁男女啪啪无遮挡网站| 怎么达到女性高潮| 1024香蕉在线观看| 日韩欧美国产在线观看| 欧美黑人精品巨大| 99热这里只有精品一区 | 久久九九热精品免费| 一边摸一边做爽爽视频免费| 麻豆av在线久日| 成年版毛片免费区| cao死你这个sao货| 欧美 亚洲 国产 日韩一| 免费观看人在逋| 国产精品亚洲美女久久久| 亚洲无线在线观看| 久久 成人 亚洲| 哪里可以看免费的av片| 观看免费一级毛片| 丰满的人妻完整版| 国产久久久一区二区三区| 亚洲片人在线观看| 法律面前人人平等表现在哪些方面| 色哟哟哟哟哟哟| 18禁黄网站禁片免费观看直播| 午夜精品在线福利| 伊人久久大香线蕉亚洲五| 国产成人aa在线观看| 两个人的视频大全免费| 在线观看免费午夜福利视频| 欧美乱色亚洲激情| 岛国在线免费视频观看| 亚洲欧美精品综合一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 日本成人三级电影网站| 国产亚洲精品第一综合不卡| 在线永久观看黄色视频| 亚洲av电影不卡..在线观看| 国产精品综合久久久久久久免费| 1024视频免费在线观看| 在线国产一区二区在线| 久久久精品欧美日韩精品| 看片在线看免费视频| av欧美777| 久99久视频精品免费| 一边摸一边做爽爽视频免费| 国产亚洲精品久久久久5区| 日韩大尺度精品在线看网址| 欧美又色又爽又黄视频| x7x7x7水蜜桃| 亚洲国产精品sss在线观看| 免费在线观看黄色视频的| 色哟哟哟哟哟哟| 在线观看日韩欧美| 小说图片视频综合网站| 亚洲成av人片免费观看| 国产一区二区三区视频了| 久久精品国产亚洲av高清一级| 欧美精品啪啪一区二区三区| 久久久精品欧美日韩精品| 欧美绝顶高潮抽搐喷水| 九色成人免费人妻av| 久久精品国产综合久久久| 国产精品av久久久久免费| 一级毛片女人18水好多| 日韩三级视频一区二区三区| svipshipincom国产片| 99re在线观看精品视频| 操出白浆在线播放| av天堂在线播放| www国产在线视频色| 亚洲成人久久性| 国产成人av教育| 2021天堂中文幕一二区在线观| 好男人在线观看高清免费视频| 天天躁夜夜躁狠狠躁躁| 欧美乱妇无乱码| 精品电影一区二区在线| 两性夫妻黄色片| 亚洲成人久久性| 午夜老司机福利片| 又粗又爽又猛毛片免费看| e午夜精品久久久久久久| 欧美又色又爽又黄视频| 亚洲一区高清亚洲精品| 日韩精品青青久久久久久|