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

    基于WRF-波浪風(fēng)暴潮耦合模型的深圳市極端風(fēng)暴潮安全分析

    2022-07-01 06:24:44林澤輝
    黑龍江水利科技 2022年5期
    關(guān)鍵詞:風(fēng)暴潮風(fēng)場(chǎng)波浪

    林澤輝

    (惠州市惠陽(yáng)區(qū)永良堤圍管理所,廣東 惠州 516007)

    我國(guó)東南沿海地區(qū)氣象水文條件復(fù)雜,水文方面,海水水位受天體引潮力作用而周期性漲落。氣象方面,西太平洋熱帶地區(qū)氣壓受到擾動(dòng)作用成渦,并不斷發(fā)展移動(dòng),這些氣旋登陸后會(huì)對(duì)人類社會(huì)造成巨大的破壞。另一方面,臺(tái)風(fēng)所造成的大氣壓強(qiáng)場(chǎng)變化也會(huì)導(dǎo)致海面異常起伏,產(chǎn)生風(fēng)暴潮,當(dāng)風(fēng)暴潮恰好疊加上天文大潮,海平面會(huì)巨幅上升,這對(duì)海岸防護(hù)、河口擋潮等水利工程帶來(lái)了巨大挑戰(zhàn)。因此研究臺(tái)風(fēng)過(guò)程前后的波浪以及增水具有重要意義。文章利用…模型對(duì)深圳地區(qū)水域在臺(tái)風(fēng)作用下的水動(dòng)力特征進(jìn)行了數(shù)值模擬研究。

    1 研究方法

    1.1 區(qū)域概況

    深圳市海岸線長(zhǎng)約250km,東起大亞灣,西至珠江口,增水>50cm的風(fēng)暴潮過(guò)程年平均出現(xiàn)2.1次。沿岸分布有機(jī)場(chǎng)、重點(diǎn)港口、大型住宅區(qū)、核電站等大型企業(yè),人口密度極高。

    深圳市年平均風(fēng)速2.6m/s,冬季各月風(fēng)速較大為3.0m/s左右,夏季各月較小約2.0m/s。大風(fēng)日數(shù)(風(fēng)速≥17m/s,即風(fēng)力相當(dāng)于8級(jí)或8級(jí)以上)年平均為7.3d。全年各月均出現(xiàn)大風(fēng),其中以7-9月居多,共占全年大風(fēng)日數(shù)的44%,其次為6月和10月,由于受臺(tái)風(fēng)或熱帶低壓影響,常出現(xiàn)大風(fēng),因此,大風(fēng)出現(xiàn)多的月份,也是臺(tái)風(fēng)季節(jié)。

    1.2 數(shù)據(jù)資料

    影響風(fēng)暴增水量級(jí)的因素有很多,包括臺(tái)風(fēng)氣壓場(chǎng)分布、臺(tái)風(fēng)風(fēng)場(chǎng)的方向與大小以及近岸地形變化因素。在從我國(guó)東南沿海登陸的臺(tái)風(fēng)里,影響到陽(yáng)江至汕尾海域的臺(tái)風(fēng)均可引起深圳海域的風(fēng)暴潮增水。目前在深圳海域,國(guó)家海洋局共設(shè)有兩個(gè)長(zhǎng)期驗(yàn)潮站,分別是赤灣和鹽田站,緊鄰深圳大亞灣岸線的有位于惠州市的惠州海洋站,各站位置見(jiàn)圖1。收集上述站點(diǎn)自建站以來(lái)臺(tái)風(fēng)期間的潮位資料統(tǒng)計(jì)歷史風(fēng)暴潮增水情況。廣東省水文局1964年在深圳市赤灣港內(nèi)建設(shè)有長(zhǎng)期水文觀測(cè)站,1986年國(guó)家海洋局南海分局在相距不到1km處建設(shè)了赤灣海洋站,并于1986年開(kāi)始投入業(yè)務(wù)化觀測(cè)。鹽田海洋站位于大鵬灣鹽田港海域,由國(guó)家海洋局南海分局建設(shè)和管理,2003年投入使用?;葜莺Q笳疚挥诖髞啚碁稠斘鱾?cè)荃灣港區(qū),由國(guó)家海洋局南海分局建設(shè)和管理,2006年建成投入使用,他們的分布位置如圖1所示。

    圖1 資料站位位置示意圖

    1.3 研究方法

    早期人們主要通過(guò)經(jīng)驗(yàn)總結(jié)歸納的方法研究風(fēng)暴潮,觀察臺(tái)風(fēng)天時(shí)海塘外潮位變化研究風(fēng)暴潮的成因。20世紀(jì)中期,隨著計(jì)算機(jī)技術(shù)的發(fā)展,人們開(kāi)始利用數(shù)值模擬研究風(fēng)暴潮。1956年W·Hansen首次利用計(jì)算機(jī)對(duì)歐洲北海的最大風(fēng)暴潮增水進(jìn)行了研究。隨著數(shù)值模擬技術(shù)的進(jìn)步,Jelesnianski建立了研究風(fēng)暴潮的SPLASH模式以及基于二維水動(dòng)力的SLOSH模式,在不斷地摸索研究后,該模式的預(yù)報(bào)精度不斷提高。近年來(lái),不斷發(fā)展的模型計(jì)算實(shí)現(xiàn)了對(duì)復(fù)雜地形影響的考慮,提高了計(jì)算的準(zhǔn)確性[1]。其次,有關(guān)風(fēng)暴增水與天文潮的非線性擬合機(jī)制也變得更加清晰,這也提高了風(fēng)暴潮數(shù)學(xué)模型的準(zhǔn)確性。我國(guó)自20世紀(jì)60年代也開(kāi)展了關(guān)于風(fēng)暴潮的相關(guān)研究。馮士筰所著的《風(fēng)暴潮導(dǎo)論》系統(tǒng)闡述了風(fēng)暴潮的概念,形成機(jī)制,并介紹了早期的風(fēng)暴潮數(shù)值預(yù)報(bào)方法。趙長(zhǎng)進(jìn)等利用ADCIRC潮流與SWAN波浪耦合模式對(duì)長(zhǎng)江口附近海域的風(fēng)暴潮過(guò)程進(jìn)行了模擬。黃本勝等用SELFE模型構(gòu)建了南?!榻涌陔p重嵌套模型,研究了“山竹”臺(tái)風(fēng)過(guò)程中珠江河口的風(fēng)暴增水與波浪增水分布。

    在波浪模擬方面,目前工程界使用較為廣泛的數(shù)學(xué)模型有兩類,一類是波浪流體動(dòng)力學(xué)模型,另一類是動(dòng)譜密度方程。緩坡方程模型與Boussinesq方程模型是主要的波浪流體力學(xué)模型。經(jīng)典的Boussinesq僅僅適用于淺水,并不適用于水深變化大的水域和不規(guī)則波的情況。因此在工程中適用范圍不廣。動(dòng)譜密度方程的基礎(chǔ)是能量平衡理論,在模擬過(guò)程中可提高海岸波浪場(chǎng)的動(dòng)力作用過(guò)程的計(jì)算效率。該方程以動(dòng)譜密度來(lái)描述波浪變形,考慮了對(duì)波-波非線性相互作用、波流相互作用、波浪破碎和底部摩擦引起的能量耗散等過(guò)程,通過(guò)將源項(xiàng)添加到方程中,可以全面考慮動(dòng)力作用的復(fù)雜過(guò)程,使其能完成大尺度區(qū)域的波浪產(chǎn)生、成長(zhǎng)、傳播和變形的計(jì)算過(guò)程。其中對(duì)物理源項(xiàng)的處理和研究是波作用量守恒的關(guān)鍵,波作用量守恒方程發(fā)展到現(xiàn)在已經(jīng)到了第三代,逐漸開(kāi)發(fā)出WAVEWATCH、SWAN、MIKE21SW、TOMAWAC等模型。

    1.4 模型設(shè)置

    1.4.1 模型介紹

    文章采用的是DELFT3D-FLOW與SWAN的耦合模型,DELFT3D是一款由荷蘭DELTRA開(kāi)發(fā)的開(kāi)源水動(dòng)力模擬軟件,在風(fēng)暴潮模擬方面已經(jīng)廣泛應(yīng)用。

    Delft3D-Flow是基于淺水方程和Boussinesq假設(shè),采用有限差分法求解斜壓Navier-Stokes方程的模型,模型控制方程組分別為:

    質(zhì)量守恒方程:

    (1)

    式中:ζ為水位(總水深);d為到參考平面的當(dāng)前水深(凈水深);U、V-x、y為方向上的速度分量;Q為單位面積上質(zhì)量強(qiáng)度。

    動(dòng)量守恒方程:

    (2)

    (3)

    式中:f為科氏力參數(shù);g為重力加速度;νh為動(dòng)態(tài)水平渦流黏度;τsx,τsy為x、y方向上的海平面風(fēng)壓分量;τbx,τby為x、y方向上的底部的剪應(yīng)力分量;ρ0為參考密度;ρ'為不規(guī)則密度。

    SWAN模型的控制方程為:

    (4)

    式中:Cx,Cy,Cσ,Cθ為二維動(dòng)譜密度在x、y、σ、θ為空間上的傳播速度;σ為相對(duì)頻率;θ為波向;S為能量源項(xiàng)。

    1.4.2 模型設(shè)置

    1)模型范圍:

    文章所采用的模型囊括整個(gè)珠江三角洲區(qū)域,見(jiàn)圖2。

    圖2 模型區(qū)域示意圖

    2)地形水深:

    模擬范圍內(nèi)的地形水深采用的是中華人民共和國(guó)海事局2016年繪制的海圖插值數(shù)據(jù)。將地形水深數(shù)據(jù)轉(zhuǎn)化到統(tǒng)一基面(平均海平面)后,用IDW插值方法插值到網(wǎng)格點(diǎn)上,得到模型的水深地形。

    3)開(kāi)邊界設(shè)置:

    開(kāi)邊界采用天文潮水位,主要考慮M2、N2、S2、K2、K1、O1、P1、Q1等八個(gè)天文分潮以及MS4、M4、M6等淺水分潮共11個(gè)分潮的作用。分潮的調(diào)和分析常數(shù)由2016-2019年的潮位實(shí)測(cè)數(shù)據(jù)經(jīng)MATLAB中的T_TIDE工具包調(diào)和分析所得。模型上部采用流量邊界,輸入數(shù)據(jù)為梧州、石角、博羅三個(gè)水文站的2016-2019洪季平均流量。

    4)底摩擦設(shè)置:

    模型根據(jù)水深范圍設(shè)置不同大小的糙率,水深越大糙率越小。

    1.4.3 風(fēng)場(chǎng)輸入

    文章的風(fēng)場(chǎng)采用美國(guó)國(guó)家大氣研究中心(NCAR)等機(jī)構(gòu)共同開(kāi)發(fā)的中尺度大氣預(yù)報(bào)系統(tǒng)。WRF模式適用于臺(tái)風(fēng)模擬,可以通過(guò)嵌套網(wǎng)格運(yùn)算來(lái)達(dá)到高精度模擬的目的。本次采用了WSM6微物理方案、KainFritsch卷積云方案、Monin-Obukhov表面層方案、Noah陸面層方案、MYJ大氣邊界層方案、RRTMG輻射方案模擬臺(tái)風(fēng)。

    2 模型驗(yàn)證

    2.1 風(fēng)場(chǎng)驗(yàn)證

    首先利用實(shí)測(cè)資料對(duì)所建立的風(fēng)暴潮-波浪耦合模型進(jìn)行率定驗(yàn)證。選取1713號(hào)臺(tái)風(fēng)(2017年第13號(hào)臺(tái)風(fēng)—天鴿,其移行路徑如圖3所示)作為模型邊界條件,為了證明WRF模式的準(zhǔn)確性其在風(fēng)暴潮模擬方面的優(yōu)勢(shì),文章也用其他論文常用的經(jīng)驗(yàn)風(fēng)場(chǎng)模型做了驗(yàn)證對(duì)比。經(jīng)驗(yàn)風(fēng)場(chǎng)采用的是宮崎正衛(wèi)模型作為移動(dòng)風(fēng)場(chǎng),HOLLAND模型作為梯度風(fēng)場(chǎng),他們的公式分別為:

    圖3 “天鴿”臺(tái)風(fēng)移行路徑

    (5)

    (5)

    式中:B為HOLLAND模型氣壓參數(shù);f為科氏力參數(shù),與模擬區(qū)域緯度有關(guān);R為最大風(fēng)速半徑,與中心氣壓有關(guān)。

    圖4為WRF模型模擬風(fēng)場(chǎng)值以及經(jīng)驗(yàn)風(fēng)場(chǎng)模擬值與實(shí)測(cè)風(fēng)場(chǎng)的對(duì)比,可以看出,WRF模型模擬的風(fēng)場(chǎng)在峰值以及過(guò)程上都比經(jīng)驗(yàn)風(fēng)場(chǎng)模型更貼近實(shí)際,因此利用WRF模型生成的臺(tái)風(fēng)風(fēng)場(chǎng)作為風(fēng)暴潮-波浪耦合模型的邊界條件更為合理。

    圖4 WRF風(fēng)場(chǎng)及模型風(fēng)場(chǎng)與實(shí)測(cè)風(fēng)場(chǎng)對(duì)比

    2.2 風(fēng)暴潮及波浪率定驗(yàn)證

    以1713號(hào)臺(tái)風(fēng)(2017年第13號(hào)臺(tái)風(fēng))為邊界條件分別驅(qū)動(dòng)珠江口風(fēng)暴潮模型以及珠江口波浪-風(fēng)暴潮耦合模型,得到如圖5所示的赤灣站與鹽田站水位驗(yàn)證曲線,可以看出,考慮波浪后的水位過(guò)程曲線更符合實(shí)際,最大誤差絕對(duì)值為13cm,平均相對(duì)誤差為11.3%,標(biāo)準(zhǔn)差為2.94,符合模擬精度要求。

    圖5(a) 赤灣站不同模型下水位曲線

    圖5(b) 鹽田站不同模型下水位曲線

    3 深圳市沿岸漫堤風(fēng)險(xiǎn)評(píng)估

    3.1 可能最大風(fēng)暴潮模型計(jì)算

    為了尋求可能影響深圳市的最大風(fēng)暴潮對(duì)深圳市沿岸的影響,采用概率論法進(jìn)行計(jì)算可能最大臺(tái)風(fēng)中的近中心最大風(fēng)速Vmax,臺(tái)風(fēng)中心氣壓P0。

    利用中國(guó)氣象局(CMA)1949-2015年熱帶氣旋最佳路徑集(Best-track)數(shù)據(jù),選取出研究區(qū)域400km直徑范圍內(nèi)對(duì)本研究區(qū)域有重要影響的臺(tái)風(fēng),以所選取臺(tái)風(fēng)的最小臺(tái)風(fēng)中心氣壓值作樣本。采用極值I型分布計(jì)算1000年一遇的臺(tái)風(fēng)中心氣壓值作為可能最大臺(tái)風(fēng)的中心氣壓,計(jì)算結(jié)果如圖6(a)所示,即深圳區(qū)域1000a遇的臺(tái)風(fēng)中心氣壓P0=880hPa。

    圖6(a) 影響深圳的臺(tái)風(fēng)中心最低氣壓極值I型分布

    圖6(b) 影響深圳的臺(tái)風(fēng)中心最大風(fēng)速極值I型分布

    利用中國(guó)氣象局(CMA)1949-2015年熱帶氣旋最佳路徑集(Best-track)數(shù)據(jù),取深圳市海灣區(qū)域400km直徑范圍內(nèi)歷年路經(jīng)本海區(qū)的臺(tái)風(fēng)最大風(fēng)速值作樣本,對(duì)于沒(méi)有臺(tái)風(fēng)經(jīng)過(guò)研究區(qū)域的年份,該年Vmax取為臺(tái)風(fēng)Vmax系列中的平均值。采用極值I型分布計(jì)算1000a一遇的最大風(fēng)速值作為可能最大臺(tái)風(fēng)的最大風(fēng)速值,計(jì)算結(jié)果如圖6(b)所示,深圳區(qū)域1000a一遇的臺(tái)風(fēng)中心最大風(fēng)速Vmax=85m/s。另外,臺(tái)風(fēng)路徑選取為東南向西北由大鵬灣登陸的路徑。

    輸出堤前代表點(diǎn)(具體見(jiàn)圖7)的最高潮位與堤頂高程進(jìn)行比對(duì),了解沿岸淹沒(méi)的情況,由于模式計(jì)算結(jié)果是相對(duì)平均海平面的,為便于比對(duì),需要根據(jù)當(dāng)?shù)仄骄F矫媾c85高程的關(guān)系,將計(jì)算結(jié)果轉(zhuǎn)成以85黃海基面為起算面,其中西段提防采用赤灣站高程關(guān)系,即計(jì)算結(jié)果加上0.54m;東段提防采用鹽田站高程關(guān)系進(jìn)行轉(zhuǎn)換,即計(jì)算結(jié)果加上0.66m。表1列出堤前統(tǒng)計(jì)的最高潮位及堤頂高程。從代表點(diǎn)的最高潮位和海堤高程比較來(lái)看,在可能最大熱帶氣旋影響下,深圳沿岸均可能發(fā)生漫堤。

    圖7 堤前代表點(diǎn)位置示意圖

    表1 可能最大風(fēng)暴潮漫堤分析 (單位:m,基面:85高程)

    續(xù)表1 可能最大風(fēng)暴潮漫堤分析 (單位:m,基面:85高程)

    3.2 海堤防御標(biāo)準(zhǔn)

    對(duì)于大部分沿海地區(qū)來(lái)說(shuō),1000a一遇的標(biāo)準(zhǔn)是過(guò)高的。但是深圳作為特大型城市,其風(fēng)暴潮安全具有重要意義,因此考慮1000a一遇情形下的風(fēng)暴潮安全是有必要的,在此情形下得出的淹沒(méi)情況有利于應(yīng)急管理部門采取相應(yīng)的防范措施。

    4 結(jié) 論

    1)利用WRF中長(zhǎng)尺度氣象預(yù)報(bào)模式作為波浪-風(fēng)暴潮耦合模型的風(fēng)場(chǎng)邊界條件輸入,提高了風(fēng)暴潮模型的精確度

    2)在珠江口波浪-風(fēng)暴潮模型驗(yàn)證實(shí)測(cè)資料準(zhǔn)確的基礎(chǔ)上,利用極值I型分布確定1000a一遇臺(tái)風(fēng)參數(shù)對(duì)影響深圳市海岸堤防未來(lái)發(fā)生的最大可能風(fēng)暴潮進(jìn)行計(jì)算分析。文章研究成果可為深圳市風(fēng)暴潮安全防范提供重要技術(shù)支撐,同時(shí)也為海洋水利防災(zāi)應(yīng)急部門提供決策參考。

    猜你喜歡
    風(fēng)暴潮風(fēng)場(chǎng)波浪
    波浪谷和波浪巖
    基于FLUENT的下?lián)舯┝魅S風(fēng)場(chǎng)建模
    2012年“蘇拉”和“達(dá)維”雙臺(tái)風(fēng)影響的近海風(fēng)暴潮過(guò)程
    防范未來(lái)風(fēng)暴潮災(zāi)害的綠色海堤藍(lán)圖
    科學(xué)(2020年4期)2020-11-26 08:27:00
    基于多變量LSTM神經(jīng)網(wǎng)絡(luò)模型的風(fēng)暴潮臨近預(yù)報(bào)
    波浪谷隨想
    去看神奇波浪谷
    “最美風(fēng)場(chǎng)”的贏利法則
    能源(2017年8期)2017-10-18 00:47:39
    側(cè)向風(fēng)場(chǎng)中無(wú)人機(jī)的飛行研究
    廣東省風(fēng)暴潮時(shí)空分布特征及重點(diǎn)城市風(fēng)暴潮風(fēng)險(xiǎn)研究
    婷婷精品国产亚洲av在线| av在线播放免费不卡| 99在线人妻在线中文字幕| 九九热线精品视视频播放| 十八禁人妻一区二区| 老司机午夜福利在线观看视频| 国内久久婷婷六月综合欲色啪| 久久精品亚洲精品国产色婷小说| 婷婷精品国产亚洲av在线| 哪里可以看免费的av片| 久久午夜亚洲精品久久| 亚洲第一电影网av| 一本大道久久a久久精品| 国产午夜福利久久久久久| 欧美av亚洲av综合av国产av| 免费看美女性在线毛片视频| 国产av又大| 午夜久久久久精精品| 超碰成人久久| 神马国产精品三级电影在线观看 | 久久久国产欧美日韩av| 精品一区二区三区视频在线观看免费| 成人av一区二区三区在线看| 天堂√8在线中文| 亚洲性夜色夜夜综合| 最近在线观看免费完整版| 蜜桃久久精品国产亚洲av| 特级一级黄色大片| 一级a爱片免费观看的视频| a级毛片a级免费在线| 亚洲欧美日韩高清专用| 精品国内亚洲2022精品成人| 亚洲一区二区三区色噜噜| 亚洲美女视频黄频| 久久久久国内视频| 日本三级黄在线观看| netflix在线观看网站| 午夜视频精品福利| 国产亚洲欧美在线一区二区| 99久久国产精品久久久| 在线国产一区二区在线| 午夜福利高清视频| 露出奶头的视频| 成人av在线播放网站| 精品第一国产精品| 在线观看66精品国产| 欧美在线黄色| 男女之事视频高清在线观看| 91麻豆av在线| 亚洲色图av天堂| 丁香欧美五月| 国产精品久久久久久人妻精品电影| 久久久久久亚洲精品国产蜜桃av| 中文字幕av在线有码专区| 99久久久亚洲精品蜜臀av| 日韩欧美三级三区| 别揉我奶头~嗯~啊~动态视频| 久久人妻av系列| 精品不卡国产一区二区三区| 88av欧美| 最近在线观看免费完整版| 欧美一级a爱片免费观看看 | 老司机午夜十八禁免费视频| 99国产精品99久久久久| 人妻夜夜爽99麻豆av| 一进一出好大好爽视频| 最近在线观看免费完整版| 一级毛片女人18水好多| 两个人的视频大全免费| 88av欧美| 在线观看日韩欧美| 亚洲欧美日韩高清在线视频| 国产一级毛片七仙女欲春2| 亚洲av美国av| 午夜精品在线福利| 亚洲天堂国产精品一区在线| 国产视频一区二区在线看| 男女做爰动态图高潮gif福利片| 亚洲性夜色夜夜综合| 日日干狠狠操夜夜爽| 最近最新免费中文字幕在线| 久久久精品国产亚洲av高清涩受| 久久伊人香网站| 九色国产91popny在线| 99精品在免费线老司机午夜| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品在线美女| 亚洲乱码一区二区免费版| 午夜福利欧美成人| 亚洲激情在线av| 变态另类丝袜制服| 亚洲国产精品成人综合色| 久久精品91无色码中文字幕| 在线十欧美十亚洲十日本专区| 国产成+人综合+亚洲专区| 欧美乱码精品一区二区三区| 国产精品98久久久久久宅男小说| 三级国产精品欧美在线观看 | 日本五十路高清| 五月伊人婷婷丁香| 国产精品,欧美在线| 日韩高清综合在线| 日日夜夜操网爽| 每晚都被弄得嗷嗷叫到高潮| 中文字幕熟女人妻在线| 老司机福利观看| 欧美不卡视频在线免费观看 | 国内精品久久久久精免费| 黄色 视频免费看| 十八禁网站免费在线| 亚洲色图 男人天堂 中文字幕| 欧美久久黑人一区二区| 国内少妇人妻偷人精品xxx网站 | 国产视频内射| 老司机在亚洲福利影院| 中文字幕人妻丝袜一区二区| 99在线视频只有这里精品首页| 亚洲精品av麻豆狂野| 老熟妇仑乱视频hdxx| 黄色视频不卡| av片东京热男人的天堂| 午夜激情av网站| 在线永久观看黄色视频| 国产亚洲av嫩草精品影院| 国产探花在线观看一区二区| 欧美一区二区精品小视频在线| 国产黄a三级三级三级人| 操出白浆在线播放| 熟妇人妻久久中文字幕3abv| 亚洲国产精品999在线| 亚洲av第一区精品v没综合| 特大巨黑吊av在线直播| 白带黄色成豆腐渣| 国产精品一区二区三区四区久久| 国产麻豆成人av免费视频| 国产成人精品久久二区二区免费| 99国产极品粉嫩在线观看| 极品教师在线免费播放| 国产成人啪精品午夜网站| av福利片在线观看| 18禁国产床啪视频网站| 午夜老司机福利片| 女人爽到高潮嗷嗷叫在线视频| 色在线成人网| 两个人视频免费观看高清| 欧美成人免费av一区二区三区| 亚洲av五月六月丁香网| 99久久精品热视频| 国产成人一区二区三区免费视频网站| 中文字幕高清在线视频| 国产三级在线视频| 嫁个100分男人电影在线观看| 午夜精品一区二区三区免费看| 18禁观看日本| 国产成人精品无人区| 法律面前人人平等表现在哪些方面| 日本成人三级电影网站| 亚洲狠狠婷婷综合久久图片| 精品久久蜜臀av无| 亚洲专区国产一区二区| 一个人免费在线观看的高清视频| 又大又爽又粗| 最新在线观看一区二区三区| 亚洲一区高清亚洲精品| 久久人妻av系列| 亚洲va日本ⅴa欧美va伊人久久| 香蕉久久夜色| 国产精品1区2区在线观看.| 99国产极品粉嫩在线观看| 窝窝影院91人妻| 成人手机av| 欧美黄色淫秽网站| 亚洲专区字幕在线| 此物有八面人人有两片| 国产91精品成人一区二区三区| 国产亚洲精品久久久久5区| 午夜精品在线福利| 亚洲 欧美一区二区三区| 亚洲最大成人中文| 久久久久久九九精品二区国产 | 精品欧美国产一区二区三| 桃色一区二区三区在线观看| 国产aⅴ精品一区二区三区波| 亚洲av片天天在线观看| 不卡一级毛片| 亚洲片人在线观看| 最近视频中文字幕2019在线8| 日韩av在线大香蕉| 后天国语完整版免费观看| 成人亚洲精品av一区二区| 日日爽夜夜爽网站| 一区二区三区国产精品乱码| 国产亚洲精品一区二区www| 三级国产精品欧美在线观看 | av国产免费在线观看| 国产一区二区在线av高清观看| 国产精品久久视频播放| 国产精品美女特级片免费视频播放器 | 国内精品一区二区在线观看| 欧美精品亚洲一区二区| 99国产极品粉嫩在线观看| 88av欧美| 久久精品91无色码中文字幕| 女人被狂操c到高潮| 成人欧美大片| 变态另类成人亚洲欧美熟女| 无遮挡黄片免费观看| 亚洲片人在线观看| 91成年电影在线观看| 国产亚洲av嫩草精品影院| 午夜成年电影在线免费观看| 麻豆一二三区av精品| 精品国产超薄肉色丝袜足j| 黄色毛片三级朝国网站| 中文亚洲av片在线观看爽| 国产精品久久久av美女十八| 一进一出好大好爽视频| av视频在线观看入口| 午夜激情福利司机影院| 亚洲欧美精品综合久久99| 亚洲欧美日韩高清在线视频| 一二三四在线观看免费中文在| 国产三级在线视频| 欧美黄色片欧美黄色片| 国产精品久久久久久人妻精品电影| 国产亚洲欧美98| 精品久久久久久久久久免费视频| 中文字幕av在线有码专区| 日本黄色视频三级网站网址| 久久久国产欧美日韩av| av有码第一页| 十八禁人妻一区二区| 久久久久国产一级毛片高清牌| 国产激情偷乱视频一区二区| 亚洲 欧美 日韩 在线 免费| av国产免费在线观看| 久久久国产精品麻豆| av视频在线观看入口| 97超级碰碰碰精品色视频在线观看| 久久午夜综合久久蜜桃| 熟女电影av网| 国产成年人精品一区二区| 免费人成视频x8x8入口观看| 女警被强在线播放| 精品午夜福利视频在线观看一区| 免费观看精品视频网站| 成人永久免费在线观看视频| 久久久久久九九精品二区国产 | xxxwww97欧美| 亚洲天堂国产精品一区在线| 757午夜福利合集在线观看| 国产爱豆传媒在线观看 | 成人亚洲精品av一区二区| 两个人看的免费小视频| 日韩精品中文字幕看吧| 午夜免费成人在线视频| 亚洲男人天堂网一区| 亚洲国产欧美一区二区综合| 精品久久久久久久末码| xxxwww97欧美| 一级毛片女人18水好多| 日韩精品中文字幕看吧| 天天躁夜夜躁狠狠躁躁| 欧美午夜高清在线| 午夜亚洲福利在线播放| 99在线视频只有这里精品首页| 在线播放国产精品三级| 天天一区二区日本电影三级| 一进一出抽搐gif免费好疼| 欧美黄色片欧美黄色片| 亚洲黑人精品在线| 免费看日本二区| 窝窝影院91人妻| 一区二区三区激情视频| 国产精品久久久av美女十八| 久久香蕉国产精品| 国产亚洲av嫩草精品影院| 国产av一区在线观看免费| 国产免费男女视频| 亚洲全国av大片| 国产精品免费一区二区三区在线| а√天堂www在线а√下载| 国产欧美日韩一区二区精品| 欧美国产日韩亚洲一区| 亚洲最大成人中文| 久久欧美精品欧美久久欧美| 中文字幕人妻丝袜一区二区| a级毛片在线看网站| 女人被狂操c到高潮| 黄色毛片三级朝国网站| 国产亚洲av高清不卡| 美女午夜性视频免费| 亚洲av日韩精品久久久久久密| 午夜精品在线福利| 日韩欧美在线二视频| 在线十欧美十亚洲十日本专区| 亚洲精品久久成人aⅴ小说| or卡值多少钱| 性欧美人与动物交配| 欧洲精品卡2卡3卡4卡5卡区| 69av精品久久久久久| 精品熟女少妇八av免费久了| 18禁黄网站禁片午夜丰满| 搞女人的毛片| 看黄色毛片网站| 亚洲精品中文字幕一二三四区| 亚洲精华国产精华精| 亚洲国产欧美人成| 免费电影在线观看免费观看| avwww免费| 久久久久亚洲av毛片大全| 一进一出抽搐gif免费好疼| 欧美一级毛片孕妇| av有码第一页| 国产一区二区在线观看日韩 | 此物有八面人人有两片| 国产真人三级小视频在线观看| 亚洲精品国产精品久久久不卡| 国产欧美日韩一区二区精品| 久久热在线av| 一二三四在线观看免费中文在| 日韩欧美三级三区| 夜夜夜夜夜久久久久| 精品久久久久久久人妻蜜臀av| 国产精品亚洲一级av第二区| 九色成人免费人妻av| 日韩高清综合在线| 中文在线观看免费www的网站 | 午夜亚洲福利在线播放| 久久草成人影院| 老熟妇乱子伦视频在线观看| avwww免费| 日本一二三区视频观看| 制服人妻中文乱码| 色哟哟哟哟哟哟| 床上黄色一级片| 三级毛片av免费| 看片在线看免费视频| 久久久国产欧美日韩av| 国产成年人精品一区二区| 十八禁人妻一区二区| 啦啦啦韩国在线观看视频| 亚洲全国av大片| 日本a在线网址| 国产精品久久电影中文字幕| 亚洲自偷自拍图片 自拍| 国产三级在线视频| а√天堂www在线а√下载| 欧洲精品卡2卡3卡4卡5卡区| 小说图片视频综合网站| 黄片小视频在线播放| 亚洲色图 男人天堂 中文字幕| 一级毛片精品| 国产精品影院久久| 叶爱在线成人免费视频播放| 久99久视频精品免费| 99久久99久久久精品蜜桃| 在线观看一区二区三区| 午夜精品在线福利| 成人欧美大片| ponron亚洲| 国产一区二区在线av高清观看| 免费看日本二区| 成人永久免费在线观看视频| 国产又色又爽无遮挡免费看| 久久精品成人免费网站| 九色国产91popny在线| 成人国产一区最新在线观看| 亚洲国产精品成人综合色| 九色国产91popny在线| 亚洲精品在线观看二区| 嫩草影院精品99| 欧美黄色淫秽网站| 国产精品99久久99久久久不卡| 久久精品影院6| 久久久国产成人精品二区| 一级黄色大片毛片| 嫩草影院精品99| 国产成人啪精品午夜网站| 真人做人爱边吃奶动态| 亚洲国产中文字幕在线视频| 成人18禁在线播放| 国产精品一区二区免费欧美| 夜夜看夜夜爽夜夜摸| 欧美日韩黄片免| 男男h啪啪无遮挡| 国产一级毛片七仙女欲春2| 天堂动漫精品| 日本一区二区免费在线视频| 99久久综合精品五月天人人| 国产真实乱freesex| 亚洲国产欧美网| 精品人妻1区二区| 亚洲五月天丁香| 欧美最黄视频在线播放免费| 大型黄色视频在线免费观看| 国产成人精品久久二区二区91| 亚洲av成人精品一区久久| 又粗又爽又猛毛片免费看| 国产高清有码在线观看视频 | 亚洲欧美日韩高清在线视频| 亚洲av成人不卡在线观看播放网| 免费观看人在逋| 日日干狠狠操夜夜爽| 免费在线观看影片大全网站| 国产三级黄色录像| 婷婷丁香在线五月| 成人特级黄色片久久久久久久| 90打野战视频偷拍视频| 色av中文字幕| 每晚都被弄得嗷嗷叫到高潮| 久久久久久久精品吃奶| 国产精品电影一区二区三区| 国产三级在线视频| 69av精品久久久久久| 首页视频小说图片口味搜索| 国产视频一区二区在线看| 亚洲无线在线观看| 欧美国产日韩亚洲一区| 伊人久久大香线蕉亚洲五| 午夜亚洲福利在线播放| 久久久久久久久免费视频了| 久久精品人妻少妇| 国产成人精品无人区| 欧美一级a爱片免费观看看 | 国产伦人伦偷精品视频| 欧美一级a爱片免费观看看 | 首页视频小说图片口味搜索| 在线观看www视频免费| 最近视频中文字幕2019在线8| 制服人妻中文乱码| 精品久久久久久,| 久久精品成人免费网站| 一本综合久久免费| 91麻豆精品激情在线观看国产| 久久久久久亚洲精品国产蜜桃av| 一级作爱视频免费观看| 亚洲 欧美一区二区三区| 女人高潮潮喷娇喘18禁视频| 91国产中文字幕| 99国产精品一区二区三区| 中文字幕人妻丝袜一区二区| 老司机福利观看| 日韩欧美国产一区二区入口| 91国产中文字幕| 国产成人啪精品午夜网站| 麻豆一二三区av精品| 99久久精品热视频| 欧美日韩福利视频一区二区| 99国产精品一区二区三区| 丰满人妻一区二区三区视频av | 亚洲,欧美精品.| 日本一二三区视频观看| 国产午夜精品论理片| 欧美黑人精品巨大| 日本五十路高清| 免费看十八禁软件| √禁漫天堂资源中文www| 精品一区二区三区视频在线观看免费| АⅤ资源中文在线天堂| 麻豆成人午夜福利视频| 黄色毛片三级朝国网站| 国产精品影院久久| 欧美日韩乱码在线| 国产成人av教育| 三级男女做爰猛烈吃奶摸视频| 俄罗斯特黄特色一大片| 久久人妻av系列| www日本黄色视频网| 黄色片一级片一级黄色片| 欧美三级亚洲精品| 国产精品 国内视频| 搡老熟女国产l中国老女人| 99久久久亚洲精品蜜臀av| 人成视频在线观看免费观看| 欧美乱码精品一区二区三区| 欧美黑人欧美精品刺激| 日本三级黄在线观看| 国产三级在线视频| 无人区码免费观看不卡| 夜夜看夜夜爽夜夜摸| 90打野战视频偷拍视频| 亚洲第一电影网av| 99re在线观看精品视频| 制服诱惑二区| 熟女电影av网| 亚洲精华国产精华精| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩免费av在线播放| 久久精品91无色码中文字幕| 国产免费男女视频| 一个人免费在线观看的高清视频| 怎么达到女性高潮| 黄片小视频在线播放| 蜜桃久久精品国产亚洲av| 亚洲精品一区av在线观看| 国产av又大| 夜夜看夜夜爽夜夜摸| 麻豆av在线久日| 久久精品国产综合久久久| 亚洲熟妇熟女久久| 国产人伦9x9x在线观看| 丁香欧美五月| 99热6这里只有精品| 中文亚洲av片在线观看爽| 又大又爽又粗| 亚洲色图av天堂| 一级片免费观看大全| 免费无遮挡裸体视频| 久久久久久久精品吃奶| 中文资源天堂在线| 久久中文字幕一级| 久久精品国产综合久久久| 90打野战视频偷拍视频| 丁香六月欧美| 日韩三级视频一区二区三区| 亚洲av成人精品一区久久| 99热只有精品国产| 日本一本二区三区精品| 欧美午夜高清在线| 午夜成年电影在线免费观看| 97人妻精品一区二区三区麻豆| 欧美色欧美亚洲另类二区| 男人舔女人的私密视频| 久久午夜亚洲精品久久| 变态另类成人亚洲欧美熟女| 色播亚洲综合网| 亚洲国产精品合色在线| 在线观看美女被高潮喷水网站 | 中国美女看黄片| 成人欧美大片| 一a级毛片在线观看| 午夜福利视频1000在线观看| 看黄色毛片网站| 在线观看一区二区三区| 1024手机看黄色片| 日韩大尺度精品在线看网址| 欧美+亚洲+日韩+国产| 国产精品av久久久久免费| √禁漫天堂资源中文www| 好男人电影高清在线观看| 免费搜索国产男女视频| 国产激情偷乱视频一区二区| 午夜影院日韩av| 这个男人来自地球电影免费观看| 好看av亚洲va欧美ⅴa在| 色综合亚洲欧美另类图片| 真人一进一出gif抽搐免费| 国产伦在线观看视频一区| 宅男免费午夜| 女人被狂操c到高潮| 精品国产亚洲在线| 国产精品av久久久久免费| 丰满人妻熟妇乱又伦精品不卡| 一区二区三区高清视频在线| 51午夜福利影视在线观看| 午夜精品一区二区三区免费看| 亚洲美女黄片视频| 久久九九热精品免费| 欧美 亚洲 国产 日韩一| 免费电影在线观看免费观看| 国产亚洲精品久久久久久毛片| 91av网站免费观看| 国内精品一区二区在线观看| 久久中文字幕人妻熟女| 特级一级黄色大片| 色播亚洲综合网| 麻豆成人av在线观看| 国模一区二区三区四区视频 | 国产熟女xx| 亚洲九九香蕉| 国产亚洲精品av在线| 国产精品九九99| 亚洲人成网站高清观看| 亚洲精品美女久久久久99蜜臀| 成人一区二区视频在线观看| 亚洲熟妇中文字幕五十中出| 国产精品一区二区免费欧美| 十八禁网站免费在线| 国产av麻豆久久久久久久| 欧美 亚洲 国产 日韩一| 亚洲av美国av| 九色国产91popny在线| 18禁国产床啪视频网站| 欧美另类亚洲清纯唯美| 久久精品国产亚洲av香蕉五月| 国产精品久久电影中文字幕| 亚洲成人中文字幕在线播放| 校园春色视频在线观看| 久久伊人香网站| 97人妻精品一区二区三区麻豆| 夜夜夜夜夜久久久久| 午夜免费成人在线视频| 中文字幕精品亚洲无线码一区| 亚洲电影在线观看av| 国产av又大| 亚洲人成网站高清观看| 亚洲成人久久性| 女警被强在线播放| 欧美又色又爽又黄视频| 国产精品亚洲av一区麻豆| 天天躁狠狠躁夜夜躁狠狠躁| 久久草成人影院| 国产真人三级小视频在线观看| 99久久综合精品五月天人人| 黄色女人牲交| 久久热在线av| 日韩三级视频一区二区三区| 男男h啪啪无遮挡| 久久中文字幕人妻熟女| 在线十欧美十亚洲十日本专区| 狠狠狠狠99中文字幕| 欧美不卡视频在线免费观看 | 50天的宝宝边吃奶边哭怎么回事| www.999成人在线观看|