• <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
    春、夏季長江口及鄰近海域溶解甲烷的分布與釋放通量
    国产精品亚洲美女久久久| 日本精品一区二区三区蜜桃| 国产亚洲欧美在线一区二区| 久久久国产精品麻豆| 黄片小视频在线播放| 国产精品日韩av在线免费观看 | 国产高清videossex| 两个人看的免费小视频| 啦啦啦韩国在线观看视频| 91在线观看av| 久久精品亚洲精品国产色婷小说| 丁香六月欧美| 欧美一级a爱片免费观看看 | 十八禁人妻一区二区| 首页视频小说图片口味搜索| 50天的宝宝边吃奶边哭怎么回事| 欧美性长视频在线观看| 日本黄色视频三级网站网址| 巨乳人妻的诱惑在线观看| 欧美日本中文国产一区发布| 99国产精品一区二区蜜桃av| 无人区码免费观看不卡| 人人妻人人澡人人看| 中文字幕最新亚洲高清| av电影中文网址| 身体一侧抽搐| 一级a爱片免费观看的视频| 久久国产精品男人的天堂亚洲| 国产野战对白在线观看| 丁香欧美五月| 搡老岳熟女国产| 免费在线观看完整版高清| av有码第一页| 日本免费a在线| 久久久久九九精品影院| 欧美色视频一区免费| 国产欧美日韩综合在线一区二区| √禁漫天堂资源中文www| 亚洲中文日韩欧美视频| 欧美不卡视频在线免费观看 | 久久精品国产99精品国产亚洲性色 | 女生性感内裤真人,穿戴方法视频| 亚洲人成77777在线视频| a级毛片在线看网站| 亚洲国产精品成人综合色| 最近最新中文字幕大全免费视频| 午夜福利在线观看吧| 国产熟女午夜一区二区三区| 国产麻豆成人av免费视频| bbb黄色大片| 亚洲成国产人片在线观看| 国产97色在线日韩免费| 这个男人来自地球电影免费观看| 亚洲精品粉嫩美女一区| 亚洲第一青青草原| 亚洲精品国产精品久久久不卡| 免费高清在线观看日韩| 国产精品久久久久久人妻精品电影| 久久久国产成人免费| 校园春色视频在线观看| 高潮久久久久久久久久久不卡| 高潮久久久久久久久久久不卡| 免费一级毛片在线播放高清视频 | 免费高清视频大片| 黄色视频不卡| 国产一区二区三区在线臀色熟女| 日本a在线网址| 日韩大码丰满熟妇| bbb黄色大片| 久久久久久亚洲精品国产蜜桃av| 欧美性长视频在线观看| 午夜免费观看网址| 欧美成狂野欧美在线观看| 怎么达到女性高潮| av中文乱码字幕在线| av在线天堂中文字幕| 国产亚洲欧美98| 好看av亚洲va欧美ⅴa在| 欧美日本亚洲视频在线播放| 日日干狠狠操夜夜爽| 亚洲av电影不卡..在线观看| 夜夜看夜夜爽夜夜摸| 欧美黑人欧美精品刺激| 国内精品久久久久久久电影| 一本综合久久免费| 黄色女人牲交| 91国产中文字幕| 精品国产乱子伦一区二区三区| 丝袜在线中文字幕| 美女高潮到喷水免费观看| bbb黄色大片| 欧美乱码精品一区二区三区| 一区二区日韩欧美中文字幕| 一区在线观看完整版| 丁香六月欧美| 国产真人三级小视频在线观看| 一级片免费观看大全| 国内久久婷婷六月综合欲色啪| 国产亚洲av高清不卡| 国产精品二区激情视频| 日本撒尿小便嘘嘘汇集6| 亚洲一区二区三区不卡视频| 啦啦啦免费观看视频1| 亚洲aⅴ乱码一区二区在线播放 | 色av中文字幕| 少妇 在线观看| 深夜精品福利| 国产亚洲欧美98| 日韩精品中文字幕看吧| 国产精品二区激情视频| 18禁黄网站禁片午夜丰满| 亚洲av片天天在线观看| 色在线成人网| www.www免费av| 怎么达到女性高潮| 美女免费视频网站| 制服诱惑二区| 波多野结衣巨乳人妻| 国产三级在线视频| 丝袜人妻中文字幕| 午夜日韩欧美国产| 久久精品国产亚洲av香蕉五月| 午夜福利视频1000在线观看 | 搡老妇女老女人老熟妇| 国产在线精品亚洲第一网站| 亚洲 国产 在线| 欧美激情高清一区二区三区| 日本三级黄在线观看| 亚洲成人免费电影在线观看| 国产精品爽爽va在线观看网站 | 亚洲欧洲精品一区二区精品久久久| 在线观看免费视频日本深夜| 亚洲狠狠婷婷综合久久图片| 亚洲第一欧美日韩一区二区三区| 手机成人av网站| 宅男免费午夜| 两人在一起打扑克的视频| 757午夜福利合集在线观看| 久热这里只有精品99| 男人的好看免费观看在线视频 | 9色porny在线观看| 人人妻人人澡欧美一区二区 | 色哟哟哟哟哟哟| 久久精品人人爽人人爽视色| 欧美国产日韩亚洲一区| 日韩欧美国产一区二区入口| 日韩大码丰满熟妇| 无限看片的www在线观看| av视频免费观看在线观看| 男人操女人黄网站| 不卡av一区二区三区| 国产欧美日韩一区二区三| 后天国语完整版免费观看| 在线国产一区二区在线| 啦啦啦免费观看视频1| 高清毛片免费观看视频网站| 国产精品亚洲美女久久久| 精品人妻在线不人妻| 成年人黄色毛片网站| 亚洲午夜理论影院| 亚洲,欧美精品.| 黄色片一级片一级黄色片| 日韩国内少妇激情av| 国产91精品成人一区二区三区| 别揉我奶头~嗯~啊~动态视频| 香蕉久久夜色| 亚洲第一青青草原| 欧美一级毛片孕妇| 夜夜看夜夜爽夜夜摸| 淫秽高清视频在线观看| 老司机福利观看| 黄网站色视频无遮挡免费观看| 亚洲av日韩精品久久久久久密| 自拍欧美九色日韩亚洲蝌蚪91| 一二三四在线观看免费中文在| 长腿黑丝高跟| 欧美中文综合在线视频| 91字幕亚洲| av视频在线观看入口| 国产高清videossex| 一进一出抽搐gif免费好疼| 搡老妇女老女人老熟妇| 免费一级毛片在线播放高清视频 | 精品一区二区三区av网在线观看| 亚洲精品一区av在线观看| 国产精品 国内视频| 高潮久久久久久久久久久不卡| 色在线成人网| 两性夫妻黄色片| 久久精品亚洲精品国产色婷小说| 变态另类丝袜制服| 久久精品国产99精品国产亚洲性色 | 中文字幕人成人乱码亚洲影| 18禁黄网站禁片午夜丰满| 日韩精品中文字幕看吧| 国产精品自产拍在线观看55亚洲| 色尼玛亚洲综合影院| 自线自在国产av| 一个人观看的视频www高清免费观看 | 自线自在国产av| 老鸭窝网址在线观看| 啦啦啦观看免费观看视频高清 | 久久人人爽av亚洲精品天堂| 美女午夜性视频免费| 真人一进一出gif抽搐免费| 超碰成人久久| 久久久久久久午夜电影| 久久人人爽av亚洲精品天堂| 亚洲精品久久成人aⅴ小说| 日韩有码中文字幕| 午夜两性在线视频| 亚洲av第一区精品v没综合| 人人澡人人妻人| 亚洲欧美日韩无卡精品| 欧美色视频一区免费| 757午夜福利合集在线观看| 欧美日本中文国产一区发布| 亚洲天堂国产精品一区在线| 国产不卡一卡二| 别揉我奶头~嗯~啊~动态视频| 欧美激情久久久久久爽电影 | 天天一区二区日本电影三级 | 免费在线观看视频国产中文字幕亚洲| 亚洲久久久国产精品| 婷婷丁香在线五月| 国产av在哪里看| 高清毛片免费观看视频网站| 亚洲av五月六月丁香网| 99国产综合亚洲精品| 九色国产91popny在线| 成人三级黄色视频| 免费在线观看亚洲国产| 老司机靠b影院| 大型av网站在线播放| 日日摸夜夜添夜夜添小说| 一区二区日韩欧美中文字幕| 精品久久蜜臀av无| 日本五十路高清| 午夜免费激情av| 超碰成人久久| 久久中文字幕人妻熟女| 91精品国产国语对白视频| 国产高清激情床上av| 久久久精品国产亚洲av高清涩受| 悠悠久久av| a级毛片在线看网站| 亚洲国产日韩欧美精品在线观看 | 久久午夜亚洲精品久久| 一级黄色大片毛片| 搡老妇女老女人老熟妇| 咕卡用的链子| 久久久久久久午夜电影| 久久青草综合色| 国产主播在线观看一区二区| 精品午夜福利视频在线观看一区| 亚洲伊人色综图| 男人的好看免费观看在线视频 | 91九色精品人成在线观看| 99精品欧美一区二区三区四区| 欧美日韩亚洲国产一区二区在线观看| 97碰自拍视频| 69精品国产乱码久久久| 成人特级黄色片久久久久久久| 狠狠狠狠99中文字幕| 久久国产亚洲av麻豆专区| 亚洲激情在线av| 可以免费在线观看a视频的电影网站| 国产欧美日韩一区二区三区在线| 真人一进一出gif抽搐免费| 亚洲色图av天堂| 欧美日韩亚洲综合一区二区三区_| 国产精品亚洲美女久久久| 国产91精品成人一区二区三区| 成年女人毛片免费观看观看9| 日本免费一区二区三区高清不卡 | 久久欧美精品欧美久久欧美| 操出白浆在线播放| 黄色 视频免费看| 亚洲精华国产精华精| 亚洲自偷自拍图片 自拍| 亚洲精品国产精品久久久不卡| 亚洲天堂国产精品一区在线| 久久国产精品男人的天堂亚洲| 国产精品秋霞免费鲁丝片| 国产亚洲精品久久久久久毛片| 男女之事视频高清在线观看| 国产成人精品久久二区二区免费| 国产黄a三级三级三级人| 国产亚洲欧美精品永久| 亚洲欧美日韩无卡精品| 久久久久久亚洲精品国产蜜桃av| 日本五十路高清| 99精品久久久久人妻精品| 国产亚洲av高清不卡| 欧美在线一区亚洲| 少妇熟女aⅴ在线视频| 视频区欧美日本亚洲| 精品久久久久久久久久免费视频| 一进一出好大好爽视频| 高清在线国产一区| 国产精品久久久久久精品电影 | 色老头精品视频在线观看| 99国产综合亚洲精品| 久久国产亚洲av麻豆专区| 人人妻人人爽人人添夜夜欢视频| 久久国产精品人妻蜜桃| 国产私拍福利视频在线观看| 成人18禁在线播放| 国产精品久久久人人做人人爽| 国产成人欧美在线观看| 午夜福利18| 国产精品美女特级片免费视频播放器 | 日本a在线网址| 亚洲午夜理论影院| 午夜福利免费观看在线| 亚洲国产欧美网| 亚洲精品美女久久av网站| 9色porny在线观看| 99精品欧美一区二区三区四区| 国产片内射在线| av欧美777| 中文字幕色久视频| 丁香欧美五月| 精品日产1卡2卡| 精品人妻1区二区| 一进一出抽搐动态| 亚洲av成人不卡在线观看播放网| 18禁黄网站禁片午夜丰满| 悠悠久久av| 成人手机av| 亚洲欧美激情在线| 成人手机av| 久久国产精品影院| 日韩大码丰满熟妇| 亚洲电影在线观看av| 日韩高清综合在线| 黄网站色视频无遮挡免费观看| 1024视频免费在线观看| 变态另类成人亚洲欧美熟女 | 99久久久亚洲精品蜜臀av| 久久九九热精品免费| 中文亚洲av片在线观看爽| 午夜福利在线观看吧| 欧美精品亚洲一区二区| 久久久久国内视频| 激情视频va一区二区三区| 一本久久中文字幕| 精品国产亚洲在线| 久9热在线精品视频| 欧美日韩黄片免| 国产色视频综合| 久久香蕉精品热| 日韩av在线大香蕉| 在线观看66精品国产| 夜夜爽天天搞| 国产精品亚洲一级av第二区| 国产精华一区二区三区| 老汉色∧v一级毛片| 如日韩欧美国产精品一区二区三区| 夜夜看夜夜爽夜夜摸| 午夜免费激情av| 亚洲三区欧美一区| 久久九九热精品免费| 精品一区二区三区视频在线观看免费| 欧美色视频一区免费| 麻豆久久精品国产亚洲av| 好男人在线观看高清免费视频 | 99精品欧美一区二区三区四区| 欧美久久黑人一区二区| 亚洲精品美女久久av网站| 一级毛片女人18水好多| 视频区欧美日本亚洲| 看黄色毛片网站| 天堂√8在线中文| 脱女人内裤的视频| 亚洲av日韩精品久久久久久密| 一区二区三区国产精品乱码| or卡值多少钱| 亚洲专区字幕在线| 亚洲aⅴ乱码一区二区在线播放 | 欧美日韩瑟瑟在线播放| 欧美在线一区亚洲| 一级a爱视频在线免费观看| 国产精品 国内视频| 亚洲专区国产一区二区| 国产色视频综合| 欧美亚洲日本最大视频资源| 在线观看日韩欧美| 日本五十路高清| 十分钟在线观看高清视频www| 涩涩av久久男人的天堂| 女人精品久久久久毛片| 国产又色又爽无遮挡免费看| 成人18禁在线播放| 亚洲久久久国产精品| 脱女人内裤的视频| 黑丝袜美女国产一区| 亚洲国产日韩欧美精品在线观看 | bbb黄色大片| 大型黄色视频在线免费观看| 亚洲免费av在线视频| 91成人精品电影| 操美女的视频在线观看| 中文字幕人成人乱码亚洲影| 正在播放国产对白刺激| 老汉色av国产亚洲站长工具| 真人做人爱边吃奶动态| 99国产综合亚洲精品| 美国免费a级毛片| 少妇裸体淫交视频免费看高清 | 女人被狂操c到高潮| 999久久久国产精品视频| 国产国语露脸激情在线看| 婷婷精品国产亚洲av在线| 久久久久久亚洲精品国产蜜桃av| 男女下面进入的视频免费午夜 | av天堂久久9| 亚洲熟妇中文字幕五十中出| 精品久久久久久久人妻蜜臀av | 啦啦啦韩国在线观看视频| 欧美日韩中文字幕国产精品一区二区三区 | 999久久久国产精品视频| 久久这里只有精品19| 亚洲精品一区av在线观看| 精品久久蜜臀av无| 亚洲激情在线av| 国内精品久久久久精免费| 欧美日韩中文字幕国产精品一区二区三区 | 乱人伦中国视频| 少妇的丰满在线观看| 国产主播在线观看一区二区| 久久热在线av| 90打野战视频偷拍视频| 又黄又爽又免费观看的视频| 高清毛片免费观看视频网站| 大型黄色视频在线免费观看| 亚洲第一欧美日韩一区二区三区| 久久精品国产清高在天天线| 亚洲色图综合在线观看| 丝袜美腿诱惑在线| 黄色成人免费大全| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日本亚洲视频在线播放| 国产精品 欧美亚洲| 欧美黑人欧美精品刺激| 国产亚洲精品综合一区在线观看 | 很黄的视频免费| 亚洲成人久久性| 免费不卡黄色视频| 精品一品国产午夜福利视频| 两个人看的免费小视频| 露出奶头的视频| 在线十欧美十亚洲十日本专区| 亚洲国产欧美一区二区综合| 黄色 视频免费看| 99久久精品国产亚洲精品| 中文字幕人成人乱码亚洲影| 亚洲情色 制服丝袜| 久久精品人人爽人人爽视色| 亚洲一区中文字幕在线| 青草久久国产| 女人高潮潮喷娇喘18禁视频| 亚洲国产精品久久男人天堂| АⅤ资源中文在线天堂| 国产免费男女视频| 亚洲五月天丁香| 国产成人欧美| 亚洲一区高清亚洲精品| 丰满的人妻完整版| 精品久久久久久久人妻蜜臀av | 色播在线永久视频| 欧美成人性av电影在线观看| 在线观看www视频免费| 久99久视频精品免费| 在线视频色国产色| 国产日韩一区二区三区精品不卡| 国产成人精品在线电影| 老司机在亚洲福利影院| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲国产精品sss在线观看| 亚洲精品国产色婷婷电影| 深夜精品福利| 欧美一级a爱片免费观看看 | 国产成人av激情在线播放| 久久天堂一区二区三区四区| 老司机午夜福利在线观看视频| 亚洲国产毛片av蜜桃av| 十分钟在线观看高清视频www| 国产高清视频在线播放一区| 亚洲五月婷婷丁香| 久久精品91蜜桃| 成人手机av| 动漫黄色视频在线观看| 午夜免费激情av| 亚洲熟女毛片儿| 日本 欧美在线| 老司机深夜福利视频在线观看| 亚洲性夜色夜夜综合| 丝袜美足系列| 国产av精品麻豆| 在线国产一区二区在线| 免费高清视频大片| 国产精品一区二区免费欧美| 精品乱码久久久久久99久播| 性欧美人与动物交配| 国产精品久久久久久亚洲av鲁大| 久久久精品国产亚洲av高清涩受| www.自偷自拍.com| 国产麻豆69| 91字幕亚洲| 操出白浆在线播放| 大型黄色视频在线免费观看| 欧美人与性动交α欧美精品济南到| 精品少妇一区二区三区视频日本电影| 人成视频在线观看免费观看| 搡老熟女国产l中国老女人| 亚洲精品粉嫩美女一区| 看免费av毛片| 天堂动漫精品| 亚洲欧美日韩高清在线视频| 男女床上黄色一级片免费看| 亚洲国产中文字幕在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产麻豆69| 亚洲精品中文字幕在线视频| 国产av又大| 国产视频一区二区在线看| 久久精品成人免费网站| 亚洲熟女毛片儿| 禁无遮挡网站| 妹子高潮喷水视频| a在线观看视频网站| 999久久久精品免费观看国产| 亚洲人成电影观看| 久久久久久免费高清国产稀缺| 麻豆国产av国片精品| 亚洲 国产 在线| 操美女的视频在线观看| 亚洲国产毛片av蜜桃av| 精品久久久久久成人av| 日韩欧美一区二区三区在线观看| 久久伊人香网站| 日韩欧美三级三区| 可以在线观看毛片的网站| 久久久国产精品麻豆| 国产亚洲精品一区二区www| 国产伦人伦偷精品视频| 手机成人av网站| 国产成人影院久久av| 色哟哟哟哟哟哟| 国产精品美女特级片免费视频播放器 | 国产av一区二区精品久久| 成人精品一区二区免费| 一个人观看的视频www高清免费观看 | 天堂√8在线中文| 日韩精品中文字幕看吧| 18禁黄网站禁片午夜丰满| 久热这里只有精品99| 国产成人精品久久二区二区免费| 18禁裸乳无遮挡免费网站照片 | 国产欧美日韩一区二区三| 在线观看舔阴道视频| 一边摸一边抽搐一进一小说| 欧美日本亚洲视频在线播放| 我的亚洲天堂| 69av精品久久久久久| 精品电影一区二区在线| 国产高清videossex| 久久中文字幕一级| 999精品在线视频| 国产单亲对白刺激| 欧美日本中文国产一区发布| 国产精品99久久99久久久不卡| 性少妇av在线| 亚洲久久久国产精品| 亚洲午夜精品一区,二区,三区| 国产成人av教育| 91麻豆精品激情在线观看国产| 香蕉久久夜色| 99在线视频只有这里精品首页| 精品一品国产午夜福利视频| 国产亚洲精品综合一区在线观看 | 国产成人精品无人区| 18禁观看日本| 亚洲精品国产色婷婷电影| 少妇裸体淫交视频免费看高清 | 午夜日韩欧美国产| 午夜福利视频1000在线观看 | 高潮久久久久久久久久久不卡| 人成视频在线观看免费观看| 婷婷精品国产亚洲av在线| 91老司机精品| 久久亚洲真实| 91成人精品电影| 欧美最黄视频在线播放免费| 久久精品91无色码中文字幕| 涩涩av久久男人的天堂| 久久青草综合色| 国产精品久久视频播放| 精品久久久久久,| 日日摸夜夜添夜夜添小说| 99riav亚洲国产免费| 国产又爽黄色视频| 色综合婷婷激情| 午夜精品在线福利| 久久欧美精品欧美久久欧美| 久久人人爽av亚洲精品天堂| 一个人观看的视频www高清免费观看 | 大陆偷拍与自拍| 亚洲国产欧美日韩在线播放| √禁漫天堂资源中文www| 日韩欧美免费精品| 日韩免费av在线播放| 国产日韩一区二区三区精品不卡| cao死你这个sao货|