• <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在线老鸭窝| 国产精品美女特级片免费视频播放器| 99国产综合亚洲精品| 婷婷六月久久综合丁香| 亚洲一区二区三区色噜噜| 18+在线观看网站| 中文字幕免费在线视频6| 久久久国产成人免费| 男女之事视频高清在线观看| 亚洲av.av天堂| 97碰自拍视频| 特级一级黄色大片| 99热这里只有是精品50| 欧美丝袜亚洲另类 | 亚洲av电影在线进入| 精品久久久久久,| 悠悠久久av| 欧美激情久久久久久爽电影| 亚洲三级黄色毛片| 一卡2卡三卡四卡精品乱码亚洲| 国产淫片久久久久久久久 | 亚洲精品456在线播放app | 色综合婷婷激情| 内射极品少妇av片p| 黄色配什么色好看| 在线国产一区二区在线| 看黄色毛片网站| x7x7x7水蜜桃| 欧美3d第一页| 麻豆成人午夜福利视频| 青草久久国产| 麻豆国产97在线/欧美| 深夜a级毛片| 99视频精品全部免费 在线| 啦啦啦观看免费观看视频高清| 99热这里只有是精品在线观看 | 9191精品国产免费久久| 国产精品久久电影中文字幕| 男人和女人高潮做爰伦理| 深爱激情五月婷婷| 丰满的人妻完整版| 国产乱人视频| 欧美不卡视频在线免费观看| 淫妇啪啪啪对白视频| 夜夜夜夜夜久久久久| 欧美成人一区二区免费高清观看| 久久久久久国产a免费观看| 3wmmmm亚洲av在线观看| 国产极品精品免费视频能看的| 韩国av一区二区三区四区| 99久久成人亚洲精品观看| 午夜精品在线福利| 亚洲片人在线观看| 欧美日本视频| 国产v大片淫在线免费观看| 国产亚洲av嫩草精品影院| 美女xxoo啪啪120秒动态图 | 亚洲av不卡在线观看| 国产亚洲av嫩草精品影院| 亚洲性夜色夜夜综合| 久久婷婷人人爽人人干人人爱| 91字幕亚洲| 一级黄色大片毛片| 亚洲午夜理论影院| 亚洲天堂国产精品一区在线| 亚洲成a人片在线一区二区| 午夜福利在线观看吧| 亚洲狠狠婷婷综合久久图片| 亚洲美女视频黄频| 日韩欧美国产一区二区入口| 亚洲专区国产一区二区| 韩国av一区二区三区四区| 国产 一区 欧美 日韩| 国产69精品久久久久777片| 嫁个100分男人电影在线观看| 欧美一区二区亚洲| 99热只有精品国产| 久久人人精品亚洲av| 国产亚洲精品久久久久久毛片| 91狼人影院| 中文字幕久久专区| 日韩欧美在线二视频| 国内精品久久久久久久电影| 桃红色精品国产亚洲av| 天美传媒精品一区二区| 哪里可以看免费的av片| 久久这里只有精品中国| 99国产极品粉嫩在线观看| 日本黄色视频三级网站网址| 免费看光身美女| 国产欧美日韩精品一区二区| 直男gayav资源| 又爽又黄无遮挡网站| 简卡轻食公司| 欧美xxxx黑人xx丫x性爽| 免费观看人在逋| 麻豆国产97在线/欧美| 一进一出好大好爽视频| 久久国产乱子伦精品免费另类| 麻豆久久精品国产亚洲av| 国产乱人伦免费视频| 色综合婷婷激情| 女人十人毛片免费观看3o分钟| 亚洲精品色激情综合| 最新中文字幕久久久久| 精品人妻1区二区| 精品一区二区三区人妻视频| 婷婷精品国产亚洲av| 特级一级黄色大片| 精品久久久久久久人妻蜜臀av| 老女人水多毛片| 精品99又大又爽又粗少妇毛片 | 99精品久久久久人妻精品| 成人无遮挡网站| 狠狠狠狠99中文字幕| 毛片一级片免费看久久久久 | a在线观看视频网站| 别揉我奶头~嗯~啊~动态视频| 免费观看的影片在线观看| 日韩高清综合在线| 欧美成人一区二区免费高清观看| aaaaa片日本免费| 国产人妻一区二区三区在| 欧美激情在线99| 狠狠狠狠99中文字幕| 男女之事视频高清在线观看| 深爱激情五月婷婷| 我要看日韩黄色一级片| 国产一区二区三区在线臀色熟女| 女生性感内裤真人,穿戴方法视频| 日韩欧美三级三区| 国产私拍福利视频在线观看| 国产高清激情床上av| 亚洲美女视频黄频| 草草在线视频免费看| 婷婷六月久久综合丁香| 婷婷六月久久综合丁香| 国产伦精品一区二区三区视频9| 热99在线观看视频| 热99在线观看视频| 午夜a级毛片| 男女下面进入的视频免费午夜| 午夜福利欧美成人| 久久精品久久久久久噜噜老黄 | 麻豆一二三区av精品| 国产在线男女| 欧美高清成人免费视频www| 午夜免费男女啪啪视频观看 | 一区二区三区四区激情视频 | av女优亚洲男人天堂| 国产人妻一区二区三区在| 熟女人妻精品中文字幕| 国产淫片久久久久久久久 | 夜夜躁狠狠躁天天躁| 在线国产一区二区在线| 欧美日本视频| 国产探花极品一区二区| 757午夜福利合集在线观看| 国产精品国产高清国产av| 婷婷精品国产亚洲av| 九九在线视频观看精品| 国产午夜福利久久久久久| 两个人的视频大全免费| 亚洲专区中文字幕在线| 神马国产精品三级电影在线观看| 免费看a级黄色片| 丁香欧美五月| 日本精品一区二区三区蜜桃| 女生性感内裤真人,穿戴方法视频| 国产91精品成人一区二区三区| 久久精品综合一区二区三区| 搡老妇女老女人老熟妇| 久久性视频一级片| 亚洲成人久久爱视频| 亚洲成人久久爱视频| 欧美色欧美亚洲另类二区| 色av中文字幕| 亚洲美女搞黄在线观看 | 91av网一区二区| 国产精品,欧美在线| 九九久久精品国产亚洲av麻豆| 色噜噜av男人的天堂激情| 少妇丰满av| 99久久精品热视频| 在线免费观看的www视频| 国产精品免费一区二区三区在线| 97超视频在线观看视频| 日韩欧美三级三区| www.999成人在线观看| 国产国拍精品亚洲av在线观看| www.999成人在线观看| 国产一区二区在线av高清观看| 午夜福利免费观看在线| 特大巨黑吊av在线直播| 两人在一起打扑克的视频| 欧美中文日本在线观看视频| 国产v大片淫在线免费观看| 国产精品久久久久久久久免 | 日日夜夜操网爽| 免费av观看视频| 亚洲一区二区三区不卡视频| 欧美绝顶高潮抽搐喷水| 日韩精品中文字幕看吧| h日本视频在线播放| 亚洲人成网站在线播放欧美日韩| 国产一区二区三区视频了| 在线a可以看的网站| 夜夜看夜夜爽夜夜摸| 国产一级毛片七仙女欲春2| 国产精品综合久久久久久久免费| 亚洲 国产 在线| 婷婷精品国产亚洲av在线| 亚洲自拍偷在线| 精品久久久久久久久亚洲 | 我的老师免费观看完整版| 人妻夜夜爽99麻豆av| 日本免费一区二区三区高清不卡| 亚洲 国产 在线| 午夜福利18| 欧美bdsm另类| www.色视频.com| 变态另类丝袜制服| 亚洲av一区综合| 性色avwww在线观看| 精品人妻熟女av久视频| 女人被狂操c到高潮| av黄色大香蕉| 黄色日韩在线| 亚洲美女视频黄频| 亚洲七黄色美女视频| 首页视频小说图片口味搜索| 日韩欧美在线乱码| 日韩人妻高清精品专区| 91麻豆精品激情在线观看国产| 亚洲aⅴ乱码一区二区在线播放| 麻豆av噜噜一区二区三区| 国产色婷婷99| 成年女人看的毛片在线观看| 国产精品av视频在线免费观看| 69人妻影院| 又爽又黄a免费视频| 12—13女人毛片做爰片一| 欧美性猛交╳xxx乱大交人| 中文字幕免费在线视频6| 欧美区成人在线视频| 欧美日韩中文字幕国产精品一区二区三区| 成人鲁丝片一二三区免费| 亚洲经典国产精华液单 | 免费在线观看日本一区| 午夜影院日韩av| 欧美成人a在线观看| 国产探花在线观看一区二区| 最好的美女福利视频网| 国产精品久久久久久人妻精品电影| 真实男女啪啪啪动态图| 国产精品影院久久| 国产精品久久久久久人妻精品电影| 真实男女啪啪啪动态图| 午夜福利在线在线| 久久亚洲真实| 精品国内亚洲2022精品成人| 波多野结衣高清无吗| 99久久无色码亚洲精品果冻| 国产 一区 欧美 日韩| 又粗又爽又猛毛片免费看| 69人妻影院| 亚洲人成网站高清观看| 精品人妻1区二区| 老鸭窝网址在线观看| 精品久久久久久久久久免费视频| av在线观看视频网站免费| 亚洲精品影视一区二区三区av| 啦啦啦韩国在线观看视频| 日本与韩国留学比较| 精品久久久久久久末码| 99精品在免费线老司机午夜| 村上凉子中文字幕在线| 久久久久九九精品影院| 国产精品不卡视频一区二区 | 久久人妻av系列| 亚洲精品在线观看二区| 色精品久久人妻99蜜桃| 麻豆一二三区av精品| 午夜a级毛片| 99精品久久久久人妻精品| 亚洲一区二区三区不卡视频| 久久精品国产亚洲av香蕉五月| 午夜精品久久久久久毛片777| 一个人看视频在线观看www免费| 国产精品,欧美在线| xxxwww97欧美| 又紧又爽又黄一区二区| 午夜福利免费观看在线| 亚洲美女搞黄在线观看 | 高潮久久久久久久久久久不卡| 日韩成人在线观看一区二区三区| 十八禁人妻一区二区| 国产熟女xx| 国产精品,欧美在线| 桃色一区二区三区在线观看| 日本成人三级电影网站| or卡值多少钱| 亚洲第一电影网av| 国产精品影院久久| 久久久久国内视频| 国产亚洲精品av在线| 久久久久久久久大av| aaaaa片日本免费| 欧美区成人在线视频| 国产精品爽爽va在线观看网站| 亚洲国产精品久久男人天堂| 好男人电影高清在线观看| 亚洲黑人精品在线| 真人做人爱边吃奶动态| 草草在线视频免费看| 99热6这里只有精品| 久久国产精品影院| 国产精品精品国产色婷婷| 久久国产精品人妻蜜桃| 免费看a级黄色片| 久久九九热精品免费| 91九色精品人成在线观看| 麻豆久久精品国产亚洲av| 久久午夜亚洲精品久久| 99久久精品国产亚洲精品| 国产激情偷乱视频一区二区| 人妻夜夜爽99麻豆av| 男女床上黄色一级片免费看| 国产精华一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 夜夜躁狠狠躁天天躁| 久久香蕉精品热| 日韩欧美一区二区三区在线观看| 99国产精品一区二区蜜桃av| 色视频www国产| 亚洲午夜理论影院| 好男人电影高清在线观看| 国产成人福利小说| 国产精品嫩草影院av在线观看 | 给我免费播放毛片高清在线观看| 久久久久精品国产欧美久久久| 亚洲av二区三区四区| 日韩中字成人| 一区二区三区四区激情视频 | 国产精品久久久久久久久免 | 久久中文看片网| 老司机午夜十八禁免费视频| 成人国产综合亚洲| 黄色女人牲交| 俺也久久电影网| 又爽又黄无遮挡网站| 麻豆一二三区av精品| 亚洲av二区三区四区| 久久人人精品亚洲av| 极品教师在线免费播放| 欧美色欧美亚洲另类二区| 久久国产精品影院| 亚洲片人在线观看| 免费大片18禁| 成年免费大片在线观看| 欧美最黄视频在线播放免费| 国产野战对白在线观看| 欧美日韩中文字幕国产精品一区二区三区| 国内精品久久久久久久电影| 欧美不卡视频在线免费观看| 久久久精品欧美日韩精品| 一卡2卡三卡四卡精品乱码亚洲| 亚洲色图av天堂| 国内揄拍国产精品人妻在线| 毛片一级片免费看久久久久 | 五月伊人婷婷丁香| 久久人人爽人人爽人人片va | 国产精品久久久久久久久免 | 亚洲av.av天堂| 91久久精品国产一区二区成人| 国产成人av教育| 久久久精品大字幕| 成人高潮视频无遮挡免费网站| 老熟妇仑乱视频hdxx| 亚洲欧美日韩无卡精品| 亚洲精华国产精华精| 极品教师在线视频| 最好的美女福利视频网| av女优亚洲男人天堂| 亚洲av二区三区四区| 乱人视频在线观看| 丰满人妻熟妇乱又伦精品不卡| ponron亚洲| 给我免费播放毛片高清在线观看| 免费无遮挡裸体视频| 欧美一区二区国产精品久久精品| 又黄又爽又刺激的免费视频.| 日韩中字成人| 欧美色欧美亚洲另类二区| 国产精品综合久久久久久久免费| 亚洲最大成人av| 午夜精品一区二区三区免费看| 精品国产三级普通话版| 午夜激情福利司机影院| 中文亚洲av片在线观看爽| 久久久久久大精品| 国产久久久一区二区三区| 午夜福利高清视频| 日韩欧美免费精品| 亚洲男人的天堂狠狠| 国产乱人视频| 观看美女的网站| 欧美一级a爱片免费观看看| 全区人妻精品视频| 久久中文看片网| 一本久久中文字幕| 国产在线精品亚洲第一网站| 亚洲中文字幕一区二区三区有码在线看| 亚州av有码| 免费在线观看成人毛片| 悠悠久久av| 十八禁人妻一区二区| 国产精品国产高清国产av| 淫秽高清视频在线观看| 免费在线观看日本一区| 宅男免费午夜| 男人和女人高潮做爰伦理| 啦啦啦观看免费观看视频高清| 久久精品91蜜桃| 国产黄a三级三级三级人| 国产精品美女特级片免费视频播放器| 免费一级毛片在线播放高清视频| 国产精品综合久久久久久久免费| 搡女人真爽免费视频火全软件 | 欧美成人免费av一区二区三区| 看黄色毛片网站| 亚洲,欧美精品.| 一本一本综合久久| 91狼人影院| 偷拍熟女少妇极品色| 午夜日韩欧美国产| 色综合站精品国产| 乱码一卡2卡4卡精品| 国产免费av片在线观看野外av| 人妻久久中文字幕网| 日韩欧美国产一区二区入口| 久久6这里有精品| 丁香六月欧美| 亚洲精品日韩av片在线观看| 男插女下体视频免费在线播放| 日日夜夜操网爽| 波多野结衣高清无吗| 一夜夜www| 国产又黄又爽又无遮挡在线| 色哟哟哟哟哟哟| 欧美日韩中文字幕国产精品一区二区三区| 蜜桃亚洲精品一区二区三区| 国产免费一级a男人的天堂| 在线观看午夜福利视频| 亚洲经典国产精华液单 | 97超视频在线观看视频| 国产欧美日韩一区二区精品| 久久香蕉精品热| 国产一区二区在线av高清观看| 97超级碰碰碰精品色视频在线观看| 亚洲av二区三区四区| 国产av不卡久久| 成年女人毛片免费观看观看9| 久久香蕉精品热| 69av精品久久久久久| 午夜精品一区二区三区免费看| 精品人妻1区二区| 亚洲黑人精品在线| 亚洲在线观看片| 2021天堂中文幕一二区在线观| 男女视频在线观看网站免费| 国产成人a区在线观看| eeuss影院久久| 亚洲av免费高清在线观看| 久久香蕉精品热| 嫁个100分男人电影在线观看| 麻豆成人av在线观看| 在线观看午夜福利视频| av福利片在线观看| 国模一区二区三区四区视频| .国产精品久久| 最后的刺客免费高清国语| 成年版毛片免费区| 免费观看人在逋| 9191精品国产免费久久| 哪里可以看免费的av片| 国产aⅴ精品一区二区三区波| 久久国产乱子免费精品| 国产精品嫩草影院av在线观看 | 国产精品伦人一区二区| av在线蜜桃| 日韩av在线大香蕉| 成人精品一区二区免费| 国产欧美日韩一区二区三| 两人在一起打扑克的视频| 欧美+亚洲+日韩+国产| 99热6这里只有精品| av天堂在线播放| 美女高潮的动态| 久久国产乱子免费精品| 99久国产av精品| 在线观看一区二区三区| 久久伊人香网站| 91麻豆精品激情在线观看国产| 成人鲁丝片一二三区免费| 麻豆久久精品国产亚洲av| 欧美zozozo另类| 中文资源天堂在线| 亚洲熟妇中文字幕五十中出| 蜜桃久久精品国产亚洲av| 精品人妻1区二区| 久久久国产成人免费| 看十八女毛片水多多多| a级毛片a级免费在线| 国产又黄又爽又无遮挡在线| 亚洲国产精品sss在线观看| 国产成人aa在线观看| 一区二区三区激情视频| 女人十人毛片免费观看3o分钟| 中文字幕高清在线视频| 美女高潮喷水抽搐中文字幕| 91麻豆av在线| 欧美bdsm另类| 日本在线视频免费播放| 亚洲av.av天堂| 久久人妻av系列| 国产欧美日韩精品亚洲av| 波多野结衣巨乳人妻| 免费看美女性在线毛片视频| 成人永久免费在线观看视频| 免费观看人在逋| 特级一级黄色大片| 岛国在线免费视频观看| 国产伦在线观看视频一区| 久久久久久大精品| 久久国产精品人妻蜜桃| 精品乱码久久久久久99久播| 色综合欧美亚洲国产小说| 日韩精品青青久久久久久| 亚洲欧美日韩高清专用| 日本成人三级电影网站| 欧美在线一区亚洲| 特级一级黄色大片| a级毛片a级免费在线| 麻豆一二三区av精品| 亚洲av成人av| 亚洲午夜理论影院| 一夜夜www| 在线a可以看的网站| 在现免费观看毛片| 91午夜精品亚洲一区二区三区 | 成年女人毛片免费观看观看9| 深夜a级毛片| 成年免费大片在线观看| 国产精品久久久久久亚洲av鲁大| 一级作爱视频免费观看| 免费在线观看影片大全网站| 国产亚洲av嫩草精品影院| 国产视频一区二区在线看| 97热精品久久久久久| 久久国产精品影院| 国产精品久久视频播放| 精品福利观看| 老鸭窝网址在线观看| 婷婷亚洲欧美| 色播亚洲综合网| 亚洲经典国产精华液单 | 欧美精品啪啪一区二区三区| 天美传媒精品一区二区| 人妻制服诱惑在线中文字幕| 91在线精品国自产拍蜜月| 又黄又爽又免费观看的视频| 亚洲五月天丁香| 悠悠久久av| 男人狂女人下面高潮的视频| 九九久久精品国产亚洲av麻豆| 亚洲av成人精品一区久久| 国产精品爽爽va在线观看网站| 日本熟妇午夜| 欧美在线一区亚洲| 三级男女做爰猛烈吃奶摸视频| 欧美成狂野欧美在线观看| 亚洲一区二区三区色噜噜| 亚洲人成电影免费在线| 一区二区三区高清视频在线| 国产真实乱freesex| ponron亚洲| 免费在线观看成人毛片| 自拍偷自拍亚洲精品老妇| 成人鲁丝片一二三区免费| 精华霜和精华液先用哪个| 蜜桃久久精品国产亚洲av| 丰满的人妻完整版| 国产成人a区在线观看| 亚洲18禁久久av| 99久久精品一区二区三区| 国产精品精品国产色婷婷| 国语自产精品视频在线第100页| 欧美成狂野欧美在线观看| 人人妻,人人澡人人爽秒播| 欧美高清性xxxxhd video| 国产国拍精品亚洲av在线观看| 国产精品99久久久久久久久| 最后的刺客免费高清国语| 91午夜精品亚洲一区二区三区 | 高清在线国产一区| 男女床上黄色一级片免费看| 成人一区二区视频在线观看| 久久久精品欧美日韩精品| 免费看光身美女| 中国美女看黄片| 动漫黄色视频在线观看| 亚洲成人久久爱视频| 亚洲美女黄片视频|