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

    徐圩港區(qū)鄰近海域溢油模擬分析

    2020-07-29 09:41:40林偉波陳曉燕李婧慧
    海洋預(yù)報(bào) 2020年3期
    關(guān)鍵詞:模型

    林偉波,陳曉燕,李婧慧

    (1.江蘇省海涂研究中心,江蘇南京210036;2.江蘇省海洋環(huán)境監(jiān)測(cè)預(yù)報(bào)中心,江蘇南京210036)

    1 引言

    近年來(lái),隨著我國(guó)近海及海上船舶施工和油氣儲(chǔ)運(yùn)設(shè)施增多,油品泄漏事故增加,不僅給企業(yè)帶來(lái)重大經(jīng)濟(jì)損失,而且對(duì)海洋生態(tài)環(huán)境造成巨大威脅。海上溢油事故發(fā)生后,油膜在海域內(nèi)的運(yùn)動(dòng)及變化是一個(gè)極其復(fù)雜的過(guò)程,它是海洋環(huán)境問(wèn)題的難點(diǎn)之一。國(guó)外自20世紀(jì)60年代開(kāi)始溢油輸運(yùn)擴(kuò)散數(shù)值模擬預(yù)報(bào)方法的研究[1]。Reed等[2]對(duì)20世紀(jì)溢油模型的研究和發(fā)展進(jìn)行了回顧,主要經(jīng)歷了油膜擴(kuò)展理論、對(duì)流擴(kuò)散方法和油粒子法3個(gè)階段。“油粒子”模式在溢油模擬發(fā)展過(guò)程中具有劃時(shí)代的意義,也成為了當(dāng)今主流的溢油模式[3]。國(guó)內(nèi)的學(xué)者在海南島南部[4]、福建東吾洋灣[5]、福寧灣[6]、浙江樂(lè)清灣[7]、長(zhǎng)江口[8]、江蘇濱海[9]、膠州灣及渤海灣等[10-12]不同海域運(yùn)用該方法開(kāi)展了大量的溢油數(shù)值模擬預(yù)報(bào)研究工作。

    本文基于丹麥水利研究所(Danish Hydraulic Institute,DHI)開(kāi)發(fā)的MIKE軟件中非結(jié)構(gòu)化網(wǎng)格的水動(dòng)力模型對(duì)徐圩港區(qū)附近海域進(jìn)行潮流數(shù)值模擬,在此基礎(chǔ)上通過(guò)“油粒子”模型對(duì)該區(qū)進(jìn)行溢油模擬,并分析不同情況下溢油油膜的漂移軌跡、掃海面積,為海洋環(huán)境影響評(píng)價(jià)和溢油事故應(yīng)急措施的制定提供科學(xué)依據(jù)。

    2 材料與方法

    2.1 水動(dòng)力模型控制方程

    MIKE21 HD是一個(gè)可用于河口、海灣以及近海潮位及潮流的綜合數(shù)值模型,對(duì)于二維流場(chǎng)模擬已經(jīng)被應(yīng)用于許多工程之中。本文使用該模塊對(duì)徐圩港區(qū)進(jìn)行了水動(dòng)力模擬實(shí)驗(yàn),潮流數(shù)據(jù)驗(yàn)證良好后輸入溢油模塊,作為溢油水動(dòng)力場(chǎng)的驅(qū)動(dòng)條件。水動(dòng)力模型用連續(xù)方程和動(dòng)量守恒方程來(lái)描述流場(chǎng)和潮位的變化。

    連續(xù)方程:

    x向動(dòng)量方程:

    y向動(dòng)量方程:

    式中:x,y為直角坐標(biāo)系坐標(biāo);t為時(shí)間;ζ為相對(duì)某一基面的水位(單位:m);h為相對(duì)某一基面的水深(單位:m);u,v分別為x,y方向上的垂線平均流速分量;Nx為x向水流紊動(dòng)粘性系數(shù)(單位:m2/s);Ny為y向水流紊動(dòng)粘性系數(shù)(單位:m2/s);f為科氏力系數(shù),f=2ωsinφ,φ為緯度;ω為地球自轉(zhuǎn)速度;fb底部摩阻系數(shù),fb=g/C2,其中C為謝才系數(shù),g為重力加速度。

    2.2 溢油控制方程

    (1)隨機(jī)游動(dòng)與擴(kuò)散過(guò)程

    在流場(chǎng)計(jì)算的基礎(chǔ)上,將溢油分成有限個(gè)質(zhì)點(diǎn),每個(gè)質(zhì)點(diǎn)的漂移位置按下式計(jì)算:

    式中:X0、Y0為某質(zhì)點(diǎn)的初始坐標(biāo);U、V為X、Y方向的流速分量;W10為海面10 m高處風(fēng)速;A為風(fēng)向;α為修正系數(shù);r為隨機(jī)擴(kuò)散系數(shù),r=R·E,R為0~1之間的隨機(jī)數(shù),E為擴(kuò)散系數(shù);B為隨機(jī)擴(kuò)散方向,B=2πR。

    (2)延展過(guò)程

    延展過(guò)程決定了表面浮油的面積擴(kuò)展,從而進(jìn)一步影響水面油膜的蒸發(fā)、溶解、擴(kuò)散和光氧化作用。延展是湍流擴(kuò)散以及重力、慣性、黏性和表面張力平衡的聯(lián)合作用結(jié)果。采用修正的Fay理論基礎(chǔ)上的重力-粘力公式計(jì)算油膜擴(kuò)展:

    式中:Aoil為油膜面積,為油膜半徑;Kα為系數(shù)(率定為0.6);t為時(shí)間;油膜體積Voil為:

    式中:hs為油膜初始厚度。

    (3)蒸發(fā)過(guò)程

    蒸發(fā)過(guò)程可導(dǎo)致10%~30%的浮油從水面進(jìn)入大氣,具體百分比取決于油種。簡(jiǎn)單時(shí)間型蒸發(fā)計(jì)算式如下:

    式中:A為油特征常數(shù);B為油的溫度特征常數(shù);T為油溫,單位:℃;t為油齡,單位:min。

    (4)水體攜帶過(guò)程

    水面浮油暴露在風(fēng)和浪中,浮油會(huì)被攜帶或擴(kuò)散進(jìn)入水體。水體攜帶是一種物理過(guò)程,在破碎浪的作用下,油滴從水面遷移到水體中。水體攜帶速率Qd(kg/(m2·s))和油粒子大小之間的關(guān)系見(jiàn)下式:

    式中:C*為與油種和風(fēng)化狀態(tài)相關(guān)的水體攜帶速率經(jīng)驗(yàn)常數(shù);Dd為單位表面積耗散的破碎波能量,單位:J/m2;S為浮油覆蓋的水面面積分?jǐn)?shù);F為受破碎浪侵襲的水面面積分?jǐn)?shù);d為油粒子直徑,單位:m;Δd為油粒子直徑差,單位:m。

    (5)乳化過(guò)程

    水-油乳化物(或稱為乳膠狀物)的形成取決于油的組分和水環(huán)境條件。乳化油可能有80%是以連續(xù)相油存在的微米級(jí)油粒子。一般乳化油的黏度要高于形成乳化油之前的油品黏度。由于水的混入,油/水混合物的體積明顯加大。水混入油相的速率計(jì)算方法見(jiàn)下式:

    式中:U為風(fēng)速,單位:m/s;Kem為乳化率常數(shù);Yw為水分?jǐn)?shù);Ymax為水在油相中的最大比例。

    2.3 水動(dòng)力計(jì)算參數(shù)及網(wǎng)格

    模型計(jì)算域從南至北包含了鹽城、連云港、日照所在海域,南北跨度約130 km,東至30 m等深線處。將其劃分成非結(jié)構(gòu)三角形網(wǎng)格,并在徐圩港區(qū)內(nèi)進(jìn)行了加密處理,在港區(qū)內(nèi)部附近分辨率最大達(dá)到30 m,港池內(nèi)水深在1~5 m之間,整個(gè)計(jì)算區(qū)域包括三角形網(wǎng)格39 544個(gè),網(wǎng)格節(jié)點(diǎn)數(shù)20 554個(gè)(見(jiàn)圖1)。開(kāi)邊界上使用東中國(guó)海模型預(yù)報(bào)得到邊界水位進(jìn)行驅(qū)動(dòng),閉邊界采用干濕邊界處理。模型使用冷啟動(dòng),時(shí)間步長(zhǎng)為60 s,每小時(shí)輸出潮流場(chǎng)數(shù)據(jù),計(jì)算10 d,取后6 d穩(wěn)定流場(chǎng)進(jìn)行驗(yàn)證分析,得到該區(qū)水動(dòng)力場(chǎng)。

    2.4 溢油模型參數(shù)設(shè)置

    圖1 模型計(jì)算網(wǎng)格(單位:m)

    表1 溢油風(fēng)況及潮時(shí)組合

    工程實(shí)施過(guò)程中的施工船舶燃料油泄漏會(huì)造成海域泄漏事故,不同船型、不同事故情況的溢油量都不同,具有較大的隨機(jī)性。此次模擬的溢油油品設(shè)定為燃料油,密度為920 kg/m3,最大含水率為0.85,以10 t燃料油作為泄漏源強(qiáng)。溢油類(lèi)型為固定點(diǎn)源瞬時(shí)溢油,油粒子數(shù)為1 000。油膜漂移過(guò)程中抵達(dá)岸邊視為被攔截,不再計(jì)算溢油量。在本工程實(shí)施時(shí),溢油最大風(fēng)險(xiǎn)可能發(fā)生位置位于四區(qū)導(dǎo)堤龍口南側(cè)區(qū)。根據(jù)連云港海洋站1974—2005年風(fēng)速、風(fēng)頻率資料,夏季常風(fēng)向?yàn)镋,風(fēng)速為4.8 m/s,冬季常風(fēng)向?yàn)镹,風(fēng)速為4.5 m/s,最不利風(fēng)向?yàn)镾E,風(fēng)速為10.7 m/s。本文模擬了夏季常風(fēng)向、冬季常風(fēng)向及最不利風(fēng)向下漲潮和落潮時(shí)刻分別發(fā)生溢油的共計(jì)5種情景,具體溢油風(fēng)況及潮時(shí)組合見(jiàn)表1。利用水動(dòng)力結(jié)果作為溢油模型驅(qū)動(dòng),考慮風(fēng)況及溢油潮時(shí)組合,溢油模型時(shí)間步長(zhǎng)為3 600 s,每小時(shí)輸出溢油結(jié)果,模擬6個(gè)潮周期(72 h)內(nèi)油膜漂移運(yùn)動(dòng)。溢油模型考慮水平擴(kuò)散、蒸發(fā)和乳化過(guò)程,水平擴(kuò)散系數(shù)取0.01 m2/s,蒸發(fā)常數(shù)A和B取2.67和0.06,乳化率常數(shù)取10-6,最大水分?jǐn)?shù)取0.75。溢油蒸發(fā)比例隨時(shí)間變化在10%~26%之間。

    3 結(jié)果與討論

    3.1 模型驗(yàn)證

    潮位驗(yàn)證數(shù)據(jù)潮位采用2015年1月23日13時(shí)—1月29日15時(shí)(北京時(shí),下同)的現(xiàn)場(chǎng)實(shí)測(cè)資料,共設(shè)3個(gè)潮位站,分別為贛榆港區(qū)H1、西連島H2和開(kāi)山島H3站,流速及流向驗(yàn)證數(shù)據(jù)為徐圩港區(qū)V1及灌河口外部V2測(cè)站2015年1月23日13時(shí)—1月24日16時(shí)的大潮實(shí)測(cè)資料,2015年1月28日11時(shí)—1月29日15時(shí)的小潮期實(shí)測(cè)資料,潮流、潮位測(cè)站分布見(jiàn)圖2。圖中流速的正負(fù)分別代表潮流的漲落,落潮流速為正值,漲潮流速為負(fù)值。根據(jù)驗(yàn)證結(jié)果可知,該模型模擬所得潮位及流速流向與實(shí)測(cè)數(shù)據(jù)驗(yàn)證較為吻合,可較好地反映當(dāng)?shù)厮畡?dòng)力形態(tài),計(jì)算所得水動(dòng)力場(chǎng)可用于溢油模型(見(jiàn)圖3—7)。

    圖2 潮位、潮流驗(yàn)證測(cè)站位置

    圖3 贛榆港區(qū)H1潮位站潮位驗(yàn)證過(guò)程線

    圖4 西連島H2潮位站潮位驗(yàn)證過(guò)程線

    圖5 開(kāi)山島H3潮位站潮位驗(yàn)證過(guò)程線

    圖6 V1測(cè)站流速與流向計(jì)算值與實(shí)測(cè)值的比較

    圖7 V2測(cè)站流速與流向計(jì)算值與實(shí)測(cè)值的比較

    3.2 水動(dòng)力模擬結(jié)果

    油類(lèi)在海面漂移,表層流場(chǎng)為其主要的驅(qū)動(dòng)力,圖8為大潮期漲落潮表層流場(chǎng)分布圖,漲潮時(shí),外海潮流基本以NE~SW方向進(jìn)入海州灣;落潮時(shí),潮流則基本以SW~NE向退出海州灣;潮流的流向與等深線或岸線的交角較大,即潮流的沿岸運(yùn)動(dòng)趨勢(shì)較小,而以離岸、向岸的往復(fù)運(yùn)動(dòng)為主。

    3.3 溢油模擬結(jié)果

    表2和圖9—11顯示了當(dāng)船舶分別在漲潮階段及落潮階段發(fā)生泄漏時(shí),油膜粒子經(jīng)72 h后的漂移軌跡及掃海范圍。

    通過(guò)預(yù)測(cè)可以看出,在夏季常風(fēng)條件下,當(dāng)泄漏發(fā)生在漲潮階段,油膜向北側(cè)圍堤內(nèi)漂移,72 h后油膜掃海面積0.6 km2;當(dāng)泄漏發(fā)生在落潮階段,油膜在落潮流和風(fēng)作用下,油膜先向西北漂移出防波堤口門(mén)后,再往西漂移,72 h后油膜掃海面積約15.8 km2。

    表2 溢油風(fēng)險(xiǎn)預(yù)測(cè)分析

    圖8 大潮流場(chǎng)分布(單位:m/s)

    圖9 工程局部大潮流場(chǎng)分布(單位:m/s)

    圖10 夏季常風(fēng)條件下1~72 h掃海面積和漂移軌跡

    在冬季常風(fēng)條件下,當(dāng)泄漏發(fā)生在漲潮階段,油膜向北側(cè)圍堤內(nèi)漂移,72 h后油膜掃海面積2.1 km2;當(dāng)泄漏發(fā)生在落潮階段,油膜在落潮流和風(fēng)作用下,油膜先向西北側(cè)漂移,出口門(mén)后,在N風(fēng)的作用下,往口門(mén)內(nèi)側(cè)移動(dòng),72 h后油膜掃海面積約11.8 km2。

    從夏季和冬季風(fēng)常風(fēng)向的實(shí)驗(yàn)中可以看出漲潮時(shí)溢油僅在圍堤內(nèi)漂移,因此本文對(duì)漲潮時(shí)最不利SE風(fēng)向的情況不再計(jì)算。

    在最不利風(fēng)向條件下,當(dāng)泄漏發(fā)生在落潮階段,油膜在落潮流和風(fēng)作用下,油膜先向西北側(cè)漂移,在3.5 h后抵達(dá)東西防波堤口門(mén),之后在SE風(fēng)和潮流的共同作用下,往西北方向移動(dòng),在11 h后到達(dá)田灣核電站取水工程,此時(shí)油膜距離泄漏位置漂移路程約為25.8 km,72 h后油膜掃海面積約28.2 km2。

    圖11 冬季常風(fēng)條件下1~72 h掃海面積和漂移軌跡

    圖12 落潮期最不利風(fēng)條件下1~72 h掃海面積和漂移軌跡

    4 結(jié)論

    利用MIKE21 HD模塊建立了連云港海域水動(dòng)力模型,較好地模擬了該海域的潮流情況。在此基礎(chǔ)上,利用拉格朗日“油粒子”模型對(duì)徐圩港區(qū)附近溢油情景進(jìn)行了模擬,由數(shù)值模擬結(jié)果可得到以下結(jié)論:

    (1)油類(lèi)在水中運(yùn)動(dòng)受多種因素共同影響,不同潮時(shí)不同風(fēng)況下溢油的漂移軌跡、漂移路程及掃海面積都存在較大差異。

    (2)徐圩港區(qū)附近溢油漂移擴(kuò)散主要受潮流和風(fēng)場(chǎng)影響,72 h內(nèi)油膜最大掃海面積及漂移路程均出現(xiàn)在落潮期最不利風(fēng)時(shí)溢油,分別為28.2 km2及25.8 km。徐圩港區(qū)內(nèi)落潮時(shí)發(fā)生溢油,在夏季風(fēng)E風(fēng)和最不利風(fēng)SE風(fēng)的條件下,油膜會(huì)經(jīng)過(guò)口門(mén)飄向西北側(cè)海域,對(duì)該區(qū)生態(tài)環(huán)境會(huì)造成一定影響。

    (3)不同風(fēng)況下溢油,油膜受風(fēng)及潮流共同驅(qū)使,油膜漂移方向與表面海流和風(fēng)所引起的流速之矢量和的方向大致相同。

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    老汉色∧v一级毛片| 黄片小视频在线播放| 国产私拍福利视频在线观看| 色综合站精品国产| 久久久久九九精品影院| 国产日本99.免费观看| 欧美成人免费av一区二区三区| 蜜桃亚洲精品一区二区三区| 国产激情欧美一区二区| 国产精品1区2区在线观看.| 每晚都被弄得嗷嗷叫到高潮| 一区福利在线观看| 国产伦精品一区二区三区视频9 | 免费看十八禁软件| 成年免费大片在线观看| 99久久成人亚洲精品观看| 日本熟妇午夜| 好男人电影高清在线观看| 亚洲成人中文字幕在线播放| 一二三四社区在线视频社区8| 一边摸一边抽搐一进一小说| 婷婷丁香在线五月| 色综合婷婷激情| 757午夜福利合集在线观看| 村上凉子中文字幕在线| 亚洲,欧美精品.| 噜噜噜噜噜久久久久久91| 一边摸一边抽搐一进一小说| 亚洲av电影不卡..在线观看| 婷婷精品国产亚洲av| 国产三级中文精品| 日本精品一区二区三区蜜桃| 色精品久久人妻99蜜桃| 蜜桃亚洲精品一区二区三区| 免费观看精品视频网站| 国产精品99久久99久久久不卡| 亚洲电影在线观看av| 两个人的视频大全免费| 国产91精品成人一区二区三区| 嫩草影院精品99| 中国美女看黄片| 老鸭窝网址在线观看| tocl精华| 亚洲天堂国产精品一区在线| 老司机在亚洲福利影院| 成人三级黄色视频| 嫩草影院精品99| 久久天躁狠狠躁夜夜2o2o| 亚洲男人的天堂狠狠| 看片在线看免费视频| 网址你懂的国产日韩在线| 日韩免费av在线播放| 丰满的人妻完整版| 69av精品久久久久久| 久久久久九九精品影院| 国产在线精品亚洲第一网站| 免费看光身美女| 一进一出抽搐gif免费好疼| 最新中文字幕久久久久| 女生性感内裤真人,穿戴方法视频| 丝袜美腿在线中文| av片东京热男人的天堂| 婷婷亚洲欧美| 成人三级黄色视频| 少妇裸体淫交视频免费看高清| 青草久久国产| 成熟少妇高潮喷水视频| 国产激情偷乱视频一区二区| 国产亚洲精品久久久久久毛片| 午夜福利在线在线| 19禁男女啪啪无遮挡网站| 亚洲五月天丁香| e午夜精品久久久久久久| 日本三级黄在线观看| 欧美精品啪啪一区二区三区| 色在线成人网| 最近视频中文字幕2019在线8| a在线观看视频网站| 老司机午夜福利在线观看视频| 夜夜夜夜夜久久久久| 人人妻人人看人人澡| 精品一区二区三区视频在线观看免费| 午夜免费观看网址| 日韩高清综合在线| 国产久久久一区二区三区| 白带黄色成豆腐渣| 午夜影院日韩av| 天堂√8在线中文| 欧美另类亚洲清纯唯美| 全区人妻精品视频| 日韩欧美一区二区三区在线观看| 国产精华一区二区三区| 99久久成人亚洲精品观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品一区av在线观看| 十八禁网站免费在线| 日韩亚洲欧美综合| 欧洲精品卡2卡3卡4卡5卡区| 精品福利观看| 搡女人真爽免费视频火全软件 | 一个人免费在线观看的高清视频| 国产淫片久久久久久久久 | 国产高清有码在线观看视频| 一级黄片播放器| 国产伦在线观看视频一区| 日韩欧美 国产精品| 他把我摸到了高潮在线观看| 色精品久久人妻99蜜桃| 91麻豆av在线| 日韩中文字幕欧美一区二区| 宅男免费午夜| 久久久久九九精品影院| 精品人妻一区二区三区麻豆 | 啦啦啦观看免费观看视频高清| 国产精品国产高清国产av| 男人舔奶头视频| 国产国拍精品亚洲av在线观看 | 99热精品在线国产| 国产极品精品免费视频能看的| av欧美777| а√天堂www在线а√下载| 两人在一起打扑克的视频| 99久久精品一区二区三区| 成人特级黄色片久久久久久久| 男人舔女人下体高潮全视频| 免费搜索国产男女视频| 内地一区二区视频在线| 中文字幕高清在线视频| 欧美在线黄色| 欧美日韩乱码在线| 国产精品久久久久久久久免 | 国产精品一区二区免费欧美| 中文字幕高清在线视频| 亚洲av成人不卡在线观看播放网| 一二三四社区在线视频社区8| 黄色片一级片一级黄色片| 亚洲在线观看片| 中文亚洲av片在线观看爽| 国产高清激情床上av| 观看美女的网站| 国产视频内射| 在线观看舔阴道视频| www日本黄色视频网| 国产国拍精品亚洲av在线观看 | 国产单亲对白刺激| 国产精品99久久久久久久久| 亚洲自拍偷在线| 免费搜索国产男女视频| 法律面前人人平等表现在哪些方面| 国产精品亚洲美女久久久| 99国产极品粉嫩在线观看| 国产亚洲欧美在线一区二区| 亚洲人与动物交配视频| 在线观看66精品国产| 最新中文字幕久久久久| 色精品久久人妻99蜜桃| 黑人欧美特级aaaaaa片| 亚洲精品在线美女| 观看美女的网站| 亚洲一区高清亚洲精品| 日韩精品中文字幕看吧| 国产精品香港三级国产av潘金莲| 麻豆成人午夜福利视频| 啦啦啦观看免费观看视频高清| 国产99白浆流出| 亚洲,欧美精品.| 男插女下体视频免费在线播放| 亚洲不卡免费看| 精品国产亚洲在线| 97超视频在线观看视频| 一本一本综合久久| 天堂av国产一区二区熟女人妻| 日本黄大片高清| 国产精品久久久人人做人人爽| 真人一进一出gif抽搐免费| 亚洲欧美日韩高清在线视频| 国产高清视频在线播放一区| 欧美bdsm另类| 亚洲精品久久国产高清桃花| 国产精品一区二区免费欧美| 在线观看av片永久免费下载| 听说在线观看完整版免费高清| 少妇的丰满在线观看| 亚洲国产精品合色在线| 久久亚洲精品不卡| 日本一二三区视频观看| 国产精品野战在线观看| 久久久成人免费电影| 中亚洲国语对白在线视频| 亚洲精品乱码久久久v下载方式 | av在线天堂中文字幕| 亚洲男人的天堂狠狠| 免费无遮挡裸体视频| 精品国产三级普通话版| 国产精华一区二区三区| 国内精品久久久久精免费| 色在线成人网| 久久久久久久精品吃奶| 国产色爽女视频免费观看| 精品国产亚洲在线| 成人特级av手机在线观看| 女人被狂操c到高潮| 757午夜福利合集在线观看| 最新中文字幕久久久久| 久久久色成人| 成人永久免费在线观看视频| 色噜噜av男人的天堂激情| 69av精品久久久久久| 亚洲精品乱码久久久v下载方式 | 精品久久久久久,| www日本黄色视频网| 久久亚洲真实| 欧美日本视频| 亚洲18禁久久av| 国产美女午夜福利| 亚洲最大成人中文| 淫秽高清视频在线观看| 亚洲一区高清亚洲精品| 国产成人a区在线观看| 欧美日韩黄片免| 床上黄色一级片| 久久久国产成人免费| 搞女人的毛片| 久久精品国产自在天天线| 啪啪无遮挡十八禁网站| 观看免费一级毛片| 中文字幕av成人在线电影| 亚洲精华国产精华精| 一本精品99久久精品77| 亚洲专区中文字幕在线| 国内毛片毛片毛片毛片毛片| 99久久精品一区二区三区| 最好的美女福利视频网| 嫩草影院入口| 淫妇啪啪啪对白视频| 无遮挡黄片免费观看| 女人十人毛片免费观看3o分钟| 蜜桃久久精品国产亚洲av| 欧美最新免费一区二区三区 | 午夜免费成人在线视频| 啪啪无遮挡十八禁网站| 亚洲欧美日韩高清在线视频| 国产精品一及| 亚洲人与动物交配视频| 不卡一级毛片| 一本久久中文字幕| 网址你懂的国产日韩在线| 热99在线观看视频| 免费在线观看亚洲国产| 老司机深夜福利视频在线观看| 午夜福利欧美成人| 亚洲国产色片| 男人和女人高潮做爰伦理| av黄色大香蕉| 国产一区二区三区视频了| 欧美一区二区国产精品久久精品| 亚洲久久久久久中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 久久精品夜夜夜夜夜久久蜜豆| 欧美黑人巨大hd| 一本精品99久久精品77| 免费看a级黄色片| 欧美性猛交╳xxx乱大交人| 99热这里只有是精品50| 又粗又爽又猛毛片免费看| 欧美性感艳星| 国产aⅴ精品一区二区三区波| 全区人妻精品视频| 欧美激情在线99| 天堂影院成人在线观看| 国产精品久久久久久亚洲av鲁大| 国产精品 国内视频| 国产久久久一区二区三区| 亚洲美女黄片视频| 99riav亚洲国产免费| 怎么达到女性高潮| av国产免费在线观看| 51国产日韩欧美| 日韩国内少妇激情av| 九九久久精品国产亚洲av麻豆| 99在线视频只有这里精品首页| 免费观看精品视频网站| 丰满人妻一区二区三区视频av | 窝窝影院91人妻| 女人高潮潮喷娇喘18禁视频| 日本一二三区视频观看| 日韩精品中文字幕看吧| 亚洲熟妇熟女久久| 亚洲av不卡在线观看| 亚洲人成网站在线播放欧美日韩| 日本精品一区二区三区蜜桃| 国产精品永久免费网站| 日韩精品青青久久久久久| 亚洲av第一区精品v没综合| 国产 一区 欧美 日韩| 欧美+亚洲+日韩+国产| 久久久久久久精品吃奶| 夜夜躁狠狠躁天天躁| 精品国产三级普通话版| 无人区码免费观看不卡| 亚洲精品一卡2卡三卡4卡5卡| 五月玫瑰六月丁香| 国产色爽女视频免费观看| 国产真实伦视频高清在线观看 | 3wmmmm亚洲av在线观看| 性欧美人与动物交配| 亚洲 国产 在线| 亚洲欧美日韩高清专用| 999久久久精品免费观看国产| 日本 欧美在线| 国产精品乱码一区二三区的特点| 香蕉丝袜av| 午夜日韩欧美国产| 午夜福利免费观看在线| 亚洲国产欧美网| 天堂√8在线中文| 男女那种视频在线观看| 精品国产三级普通话版| 法律面前人人平等表现在哪些方面| 91av网一区二区| 老司机午夜十八禁免费视频| 日本五十路高清| 亚洲人成电影免费在线| 精品电影一区二区在线| 狂野欧美白嫩少妇大欣赏| 亚洲成人久久性| 久久性视频一级片| 午夜福利在线观看免费完整高清在 | 国产av麻豆久久久久久久| 一区福利在线观看| 最近在线观看免费完整版| 国产午夜精品论理片| 宅男免费午夜| 婷婷丁香在线五月| 成人国产一区最新在线观看| 狠狠狠狠99中文字幕| 精品久久久久久久久久免费视频| 欧美bdsm另类| 精品熟女少妇八av免费久了| 两个人的视频大全免费| 国产精品久久视频播放| 亚洲午夜理论影院| 久久久久久国产a免费观看| 欧美中文综合在线视频| 欧美又色又爽又黄视频| 一区二区三区免费毛片| 国产亚洲精品久久久久久毛片| 欧美性猛交╳xxx乱大交人| 国产又黄又爽又无遮挡在线| 嫩草影视91久久| 一进一出抽搐动态| 九九在线视频观看精品| 欧美成人性av电影在线观看| 欧美绝顶高潮抽搐喷水| 久久精品影院6| 色综合欧美亚洲国产小说| 精品一区二区三区人妻视频| 国产精品日韩av在线免费观看| 免费在线观看日本一区| 3wmmmm亚洲av在线观看| 日本黄色片子视频| 精品国产三级普通话版| 一卡2卡三卡四卡精品乱码亚洲| 成人av一区二区三区在线看| 亚洲五月天丁香| 中文资源天堂在线| 亚洲18禁久久av| 久久久久久久亚洲中文字幕 | 性欧美人与动物交配| 久久精品91蜜桃| www.www免费av| 一级黄色大片毛片| 一个人观看的视频www高清免费观看| 一级毛片高清免费大全| 99热只有精品国产| 午夜精品久久久久久毛片777| 国产午夜福利久久久久久| 岛国在线免费视频观看| 在线国产一区二区在线| 成人特级av手机在线观看| 制服人妻中文乱码| 国产真人三级小视频在线观看| 给我免费播放毛片高清在线观看| 最近最新中文字幕大全电影3| 嫩草影视91久久| 亚洲av免费高清在线观看| 激情在线观看视频在线高清| 乱人视频在线观看| 国产成人a区在线观看| 亚洲av成人av| 中文字幕人妻熟人妻熟丝袜美 | 午夜老司机福利剧场| 三级男女做爰猛烈吃奶摸视频| 性色av乱码一区二区三区2| 脱女人内裤的视频| 日韩欧美三级三区| 国产成人a区在线观看| 搡女人真爽免费视频火全软件 | www.熟女人妻精品国产| 亚洲精品一区av在线观看| 一本一本综合久久| 久久香蕉精品热| 国产一区二区亚洲精品在线观看| 亚洲中文日韩欧美视频| 国产成人福利小说| 久久久精品大字幕| 男女下面进入的视频免费午夜| 亚洲国产中文字幕在线视频| 美女免费视频网站| av片东京热男人的天堂| 99国产精品一区二区三区| 成人午夜高清在线视频| 欧美绝顶高潮抽搐喷水| 乱人视频在线观看| 久久精品国产亚洲av涩爱 | ponron亚洲| 精品久久久久久久末码| www日本在线高清视频| 精品乱码久久久久久99久播| 国产成人av激情在线播放| 久久精品国产亚洲av香蕉五月| 两人在一起打扑克的视频| 成人三级黄色视频| www.熟女人妻精品国产| 亚洲,欧美精品.| 日本三级黄在线观看| 美女大奶头视频| 99精品久久久久人妻精品| 在线观看av片永久免费下载| 久9热在线精品视频| 小说图片视频综合网站| 亚洲欧美日韩卡通动漫| 午夜福利高清视频| 成人无遮挡网站| 免费无遮挡裸体视频| 最近视频中文字幕2019在线8| 午夜精品一区二区三区免费看| 久久婷婷人人爽人人干人人爱| 亚洲第一欧美日韩一区二区三区| 国产精品影院久久| 久9热在线精品视频| 亚洲精品在线观看二区| 国产高清有码在线观看视频| 欧美日本视频| 国产一区二区亚洲精品在线观看| 美女被艹到高潮喷水动态| 欧美午夜高清在线| 99久久综合精品五月天人人| 免费av毛片视频| 校园春色视频在线观看| 精品人妻1区二区| 国产精品美女特级片免费视频播放器| 麻豆久久精品国产亚洲av| 99riav亚洲国产免费| 一个人看的www免费观看视频| 听说在线观看完整版免费高清| www.999成人在线观看| 校园春色视频在线观看| 久久人人精品亚洲av| 校园春色视频在线观看| 18禁黄网站禁片免费观看直播| 天美传媒精品一区二区| 亚洲国产欧美网| 色综合站精品国产| 一级毛片高清免费大全| 日日干狠狠操夜夜爽| 欧美乱色亚洲激情| 国产伦人伦偷精品视频| 精品国内亚洲2022精品成人| 国产精品自产拍在线观看55亚洲| 国产精品野战在线观看| 国产三级黄色录像| 国内久久婷婷六月综合欲色啪| 男女床上黄色一级片免费看| 天堂网av新在线| 亚洲av免费在线观看| 99久久精品一区二区三区| 亚洲精品粉嫩美女一区| 99精品在免费线老司机午夜| 日本精品一区二区三区蜜桃| 9191精品国产免费久久| 国产精品自产拍在线观看55亚洲| 天堂av国产一区二区熟女人妻| 少妇裸体淫交视频免费看高清| 国产精品亚洲美女久久久| 国产精品精品国产色婷婷| 久久久久久久久中文| 综合色av麻豆| 欧美高清成人免费视频www| 精品久久久久久久毛片微露脸| a级一级毛片免费在线观看| 男女床上黄色一级片免费看| 国产成人a区在线观看| 国产精品一区二区三区四区久久| 精品人妻一区二区三区麻豆 | 九九在线视频观看精品| 18美女黄网站色大片免费观看| 在线a可以看的网站| 久久久久国内视频| 亚洲国产欧美网| 最近视频中文字幕2019在线8| 国产欧美日韩精品一区二区| 90打野战视频偷拍视频| 真实男女啪啪啪动态图| 天天添夜夜摸| 久久精品亚洲精品国产色婷小说| 99国产综合亚洲精品| 波野结衣二区三区在线 | 欧美另类亚洲清纯唯美| 一a级毛片在线观看| 欧美色欧美亚洲另类二区| 99久久精品热视频| 级片在线观看| 国产精品一及| 看黄色毛片网站| 欧美一区二区国产精品久久精品| 人人妻人人看人人澡| 亚洲性夜色夜夜综合| 国产视频内射| 日韩欧美精品免费久久 | 人人妻人人澡欧美一区二区| 99久久精品国产亚洲精品| 韩国av一区二区三区四区| 欧美日韩乱码在线| 国产精品国产高清国产av| 老鸭窝网址在线观看| 熟妇人妻久久中文字幕3abv| 国产亚洲精品久久久com| 老汉色∧v一级毛片| 国产欧美日韩一区二区精品| 亚洲 国产 在线| 久久久久性生活片| 午夜福利在线在线| 亚洲天堂国产精品一区在线| 日本精品一区二区三区蜜桃| 可以在线观看毛片的网站| av视频在线观看入口| 婷婷精品国产亚洲av| 最后的刺客免费高清国语| 国产老妇女一区| 免费在线观看亚洲国产| 搡老岳熟女国产| 蜜桃久久精品国产亚洲av| 亚洲国产欧美人成| 国产伦在线观看视频一区| 91在线精品国自产拍蜜月 | 欧美黑人欧美精品刺激| 啦啦啦观看免费观看视频高清| 国产一区二区激情短视频| 特级一级黄色大片| 嫩草影院精品99| 又黄又粗又硬又大视频| 久久6这里有精品| 91久久精品电影网| 国产成人av教育| 2021天堂中文幕一二区在线观| 午夜激情欧美在线| 免费av毛片视频| 国产午夜精品论理片| 欧美乱妇无乱码| 搡老岳熟女国产| 男插女下体视频免费在线播放| 精品国产亚洲在线| 老司机午夜福利在线观看视频| 免费高清视频大片| 白带黄色成豆腐渣| x7x7x7水蜜桃| 美女黄网站色视频| 天天一区二区日本电影三级| 亚洲国产色片| 国产伦精品一区二区三区视频9 | 一a级毛片在线观看| 看片在线看免费视频| 国模一区二区三区四区视频| 亚洲精品在线观看二区| 两个人看的免费小视频| 亚洲av一区综合| 身体一侧抽搐| 欧洲精品卡2卡3卡4卡5卡区| 最近视频中文字幕2019在线8| 动漫黄色视频在线观看| 成人亚洲精品av一区二区| 久久久国产精品麻豆| 在线天堂最新版资源| 少妇的丰满在线观看| 麻豆国产97在线/欧美| 在线a可以看的网站| 首页视频小说图片口味搜索| 哪里可以看免费的av片| 免费av不卡在线播放| 欧美黄色淫秽网站| 床上黄色一级片| 在线观看66精品国产| 女人被狂操c到高潮| 麻豆国产97在线/欧美| 欧美区成人在线视频| 99久国产av精品| 亚洲国产欧美人成| 亚洲中文字幕日韩| 91麻豆av在线| 露出奶头的视频| 国产一区二区三区视频了| 日韩欧美在线二视频| 国产精品女同一区二区软件 | 欧美最新免费一区二区三区 | 国产在线精品亚洲第一网站| 久久久久久久精品吃奶| 午夜免费成人在线视频| 人妻丰满熟妇av一区二区三区| 一个人观看的视频www高清免费观看| 亚洲真实伦在线观看| 啦啦啦观看免费观看视频高清| 最近最新免费中文字幕在线| 无遮挡黄片免费观看| a级一级毛片免费在线观看|