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

    WT模型在山區(qū)高速公路風(fēng)預(yù)報訂正中的應(yīng)用

    2022-04-28 09:35:42曲曉黎張杏敏尤琦張中杰吳丹
    氣象科技 2022年2期
    關(guān)鍵詞:崇禮張家口大風(fēng)

    曲曉黎 張杏敏 尤琦,3 張中杰 吳丹

    (1 河北省氣象服務(wù)中心,石家莊 050021; 2 河北省氣象與生態(tài)環(huán)境重點實驗室,石家莊 050021; 3 中國氣象局交通氣象重點開放實驗室,南京 210041;4 河北省氣象科學(xué)研究所,石家莊 050021)

    引言

    近地層風(fēng)的狀況是天氣過程形成的重要因素,也對諸多行業(yè)運營和人們的生產(chǎn)生活有著重要影響。隨著我國氣象觀測站網(wǎng)的精密程度不斷提升以及數(shù)值模擬水平的提高,對風(fēng)的觀測和預(yù)報都有了長足進步,然而復(fù)雜地形風(fēng)的狀況和預(yù)報仍是熱點和難點。宋麗莉等[1-2]多年來基于觀測事實對復(fù)雜山地及沿海地區(qū)近地層風(fēng)的特性進行了詳細(xì)的計算分析,得出了近地層大風(fēng)的演變規(guī)律和特征。賈春暉等[3]利用延慶-張家口一帶的地面氣象觀測數(shù)據(jù)和探空數(shù)據(jù)揭示了這一帶復(fù)雜地形對近地面風(fēng)場的動力和熱力作用及風(fēng)場的精細(xì)化分布特征。馮婧等[4]對京津冀地區(qū)20世紀(jì)50年代以來的日平均風(fēng)速序列進行了非均一性檢驗和均一化訂正。對于風(fēng)的數(shù)值預(yù)報,很多學(xué)者開展了大量的數(shù)值模擬試驗和訂正研究[5-10],結(jié)果表明單純依靠數(shù)值模式還有一定的局限性。近年來,CFD(Computational Fluid Dynamics)流體力學(xué)計算方法也被較多地應(yīng)用于復(fù)雜地形風(fēng)場的模擬研究。李磊等[11]對“中尺度模式/CFD”耦合過程中CFD解域頂高及入流邊界條件設(shè)置2個關(guān)鍵問題進行了詳細(xì)探討。另外,很多學(xué)者分別利用數(shù)值模擬與CFD結(jié)合的方法對我國不同區(qū)域復(fù)雜地形的風(fēng)環(huán)境進行模擬,取得了較好的模擬效果[12-17]。

    針對影響高速公路安全通行的橫風(fēng)或側(cè)風(fēng),已有關(guān)于其時空分布特征、致災(zāi)危險性評價以及行車安全性和安全速度等方面的研究[18-23]。然而,對于復(fù)雜山區(qū)高速公路沿線風(fēng)的預(yù)報少見報道。張承高速公路張家口至崇禮段全長57 km,是2022年北京冬奧會張家口賽區(qū)的重要交通線路,該路段穿過崇禮山區(qū),沿線有多處山谷,大風(fēng)條件下對高速公路行車會產(chǎn)生較大影響,而數(shù)值預(yù)報往往對這些區(qū)域的風(fēng)的預(yù)報尤其是風(fēng)速的預(yù)報能力不足。因此本文試圖應(yīng)用基于計算流體力學(xué)原理的Meteodyn WT模型(以下簡稱WT模型),以及ECMWF數(shù)值模式風(fēng)的預(yù)報產(chǎn)品,模擬輸出張承高速沿線特殊地形的風(fēng)的預(yù)報訂正結(jié)果,以便提升山區(qū)高速公路風(fēng)的精細(xì)化預(yù)報能力。

    1 資料與方法

    1.1 資料來源

    在分析大風(fēng)日數(shù)分布特征和風(fēng)速變化背景時,選取河北省氣象信息中心提供的2001—2020年全省142個國家氣象站逐日最大風(fēng)速資料;驅(qū)動WT模型的風(fēng)的數(shù)值預(yù)報選取ECMWF數(shù)值模式輸出的逐3 h地面10 m風(fēng);在對風(fēng)速進行檢驗時,選取張家口市崇禮區(qū)范圍內(nèi)崇禮國家氣象站、南窩鋪及場地2套區(qū)域氣象觀測站,高家營及西灣子2套交通氣象觀測站的分鐘級風(fēng)的觀測數(shù)據(jù)。

    1.2 WT模型原理

    WT模型是法國美迪公司研發(fā)的基于CFD計算流體力學(xué)的、主要用于野外復(fù)雜地形風(fēng)的評估和模擬的模型。具體原理如下:假設(shè)在給定空間與時間的任意一點處的風(fēng)速是由平均風(fēng)速與隨機分量疊加而成,則瞬時風(fēng)速分量可以表示為:

    (1)

    將式(1)代入大氣連續(xù)方程式(2)進行積分,可得式(3)和式(4),即定常不可壓縮質(zhì)量守恒和平均雷諾納維葉-斯托克斯方程(RANS):

    (2)

    (3)

    Fi=0

    (4)

    采用K-L湍流模型來預(yù)測雷諾應(yīng)力對平均流的影響[24]:

    (5)

    其中,Pk為湍流動能生成項,VT為湍流粘度,u、v、w分別為風(fēng)速在x、y、z方向上的分量。

    1.3 WT模型處理流程

    圖1所示,WT模型可以按照360°等分或指定的風(fēng)向角,根據(jù)指定風(fēng)廓線、地表邊界條件(地面粗糙度、森林冠層模型)、熱穩(wěn)定度等級分別進行定向模擬計算, 由此得到標(biāo)準(zhǔn)入口風(fēng)速、不同風(fēng)向條件下的風(fēng)場分布。再根據(jù)定向計算結(jié)果中測風(fēng)觀測與計算區(qū)域內(nèi)相同位置之間的統(tǒng)計關(guān)系, 推算給定地點或區(qū)域的風(fēng)速值和風(fēng)能參數(shù)。

    圖1 WT模型計算處理流程

    2 河北省大風(fēng)日數(shù)分布及日最大風(fēng)速變化特征

    利用2001—2020年河北省142個國家氣象站逐日最大風(fēng)速資料,繪制過去20年河北省各站6級及以上大風(fēng)日數(shù)分布(圖2),可見張家口西北部及承德西部是6級及以上大風(fēng)日數(shù)最多的區(qū)域,而崇禮及周邊地區(qū)大風(fēng)日數(shù)并不是很多,據(jù)分析,這可能與崇禮氣象觀測站的位置有一定關(guān)系,因此就更有必要對該區(qū)域復(fù)雜地形處的風(fēng)的情況進一步模擬和分析。對利用EOF正交函數(shù)分解的2001—2020年各站逐日最大風(fēng)速距平場進行分析發(fā)現(xiàn),前3個主要模態(tài)方差貢獻較大分別為71.3%、11.7%以及5.4%。第一模態(tài)(圖3a)在空間上呈現(xiàn)全省一致性負(fù)距平的變化特征;第二模態(tài)(圖略)南北反向,省中北部為正距平,東南部為負(fù)距平;第三模態(tài)(圖略)東西反向,西部地區(qū)為負(fù)距平。結(jié)合時間系數(shù)(圖3b)綜合來看,全省絕大部分地區(qū)在過去20年中日最大風(fēng)速呈現(xiàn)明顯減少的趨勢(這與已有研究結(jié)論一致[25-28]),其中張家口、承德北部地區(qū)在第二、第三模態(tài)中略有增大的趨勢,這說明這些區(qū)域的大風(fēng)還是需要格外關(guān)注的。

    圖2 2000—2020年河北省142個國家氣象站6級及以上大風(fēng)日數(shù)分布

    圖3 2001—2020年河北省142個國家氣象站逐日最大風(fēng)速距平場EOF第一空間模態(tài)(a)和第一時間系數(shù)(b)

    3 張承高速公路沿線山區(qū)路段風(fēng)的模擬及檢驗

    3.1 模擬區(qū)域及參數(shù)設(shè)置

    模擬的區(qū)域如圖4所示,位于張家口崇禮區(qū)的南部,處于陰山山脈東段的大馬群山和燕山余脈交接地帶,屬于壩上壩下過渡型山區(qū)的山地地形。該區(qū)域山勢陡峻,山峰海拔多在1500~2000 m之間,最大高差可達(dá)1300 m多。2個解域區(qū)域幾乎均為正方形,邊長約9.3 km,中心位置分別為位于(115.067°E,40.881°N)和(115.219°E,40.948°N),2022年北京冬奧會張家口賽區(qū)的主干要道張承高速公路沿山谷底部邊緣自西南向東北方向穿過該區(qū)域。區(qū)域內(nèi)有5處氣象觀測站,分別為:南窩鋪區(qū)域氣象觀測站(站號:B2963,海拔930 m),距高速公路約1209 m,位于區(qū)域西部較開闊西北—東南走向山谷底部;場地區(qū)域氣象觀測站(站號:B3197,海拔1050 m),距高速公路約191 m,位于區(qū)域東部接近南北走向山谷底部;崇禮國家氣象觀測站(站號:54304,海拔1240 m),距高速公路約159 m;位于張承高速公路沿線K34和K57處的高家營和西灣子交通氣象站。

    圖4 張承高速公路沿線解域區(qū)域地形分布(紅色線條為張承高速公路,藍(lán)框內(nèi)為解域區(qū)域,地形圖來自World Imagery遙感影像)

    在模擬試驗過程中,WT模型采用非結(jié)構(gòu)化網(wǎng)格,解域最小水平分辨率為25 m、水平擴展系數(shù)1.1;最小垂直分辨率為4 m、垂直擴展系數(shù)1.2,解域頂高由WT模型根據(jù)計算區(qū)域內(nèi)地形高度自動確定:

    Ztop=Zmax+max(2800,5(Zmax-Zmin))

    (6)

    其中,Ztop為頂面邊界的海拔高度,Zmax和Zmin分別為解域區(qū)域內(nèi)的地形最高和最低海拔高度。

    選用ECMWF數(shù)值模式每日2次輸出的逐3 h、分辨率為0.125°的地面10 m風(fēng)的預(yù)報產(chǎn)品來驅(qū)動WT模型,選取代表站周邊最近一個ECMWF格點預(yù)報結(jié)果作為該站的環(huán)境風(fēng)輸入WT模型,再運用WT模型進行間隔10°的定向計算,生成各方向風(fēng)加速因數(shù),輸出解域區(qū)域內(nèi)未來24 h預(yù)報時效、逐50 m空間分辨率、3 h時間分辨率的地面10 m風(fēng)產(chǎn)品。進而對張承高速公路沿線山區(qū)路段的風(fēng)的數(shù)值預(yù)報產(chǎn)品進行訂正。

    3.2 天氣過程描述

    模擬的天氣過程包括30個代表日,日極大風(fēng)的最大值出現(xiàn)在崇禮本站2021年4月16日下午,為18.8 m/s(8級)。這些大風(fēng)過程的天氣背景較為相似,500 hPa高度場(圖略)上,烏拉爾山以東地區(qū)均為一脊一槽形勢,高壓脊位于貝加爾湖西側(cè),高空槽處于110°~130°E,且大部分過程有冷渦與之配合,槽脊振幅較大,等高線密集。700 hPa與850 hPa高度場(圖略)上存在鋒區(qū),850 hPa上風(fēng)場與等溫線近乎垂直,有較強的冷平流。地面天氣圖(圖略)上,河北中北部至內(nèi)蒙一帶地面等壓線較為密集,氣壓梯度大,達(dá)到12.5~22.5 hPa/10°。

    3.3 模擬結(jié)果檢驗分析

    將30個代表日中圖4所示的5套氣象觀測站的10 min平均風(fēng)速與WT模型模擬輸出的10 min平均風(fēng)速兩組數(shù)據(jù)進行對比,有效樣本數(shù)為1167個,其相關(guān)系數(shù)達(dá)到0.6225,通過0.001的顯著性水平檢驗。繪制5個觀測站實況和WT模擬結(jié)果的散點圖(圖5)以及崇禮代表站30個大風(fēng)日的實況和模擬結(jié)果對比曲線(圖6)可見,散點相對集中且均勻的分布在趨勢線兩側(cè),模擬的數(shù)值比實際觀測的數(shù)據(jù)略偏高,模擬結(jié)果與實況的誤差在-4.69~4.12 m/s之間。具體分析各站的誤差(表1)發(fā)現(xiàn),崇禮國家站和西灣子交通站的模擬結(jié)果偏高明顯,位于山谷地帶的場地區(qū)域站和高家營交通站的模擬結(jié)果明顯偏小,而地形相對開闊的南窩鋪區(qū)域站模擬結(jié)果相對其他4站誤差范圍明顯要小。整體來看,各站模擬的結(jié)果與實況之間的誤差在±1 m/s之間的均占40%以上(圖7),±2 m/s之間的占比基本能夠達(dá)到80%以上。從模擬的風(fēng)向來看,均是較為一致的偏北風(fēng)或西北風(fēng),與天氣過程相符。

    圖5 2021年2—5月張家口、崇禮區(qū)域5個觀測站10 min平均風(fēng)速實況與WT模擬風(fēng)速

    圖6 2021年2—5月崇禮站30個大風(fēng)日逐3 h間隔10 min平均風(fēng)速實況與WT模擬風(fēng)速

    表1 2021年2—5月張家口、崇禮區(qū)域5個觀測站10 min平均風(fēng)速檢驗

    圖7 2021年2—5月張家口、崇禮區(qū)域各觀測站10 min平均風(fēng)速模擬與實況誤差占比

    由于在開展高速公路氣象服務(wù)過程中更關(guān)注極大風(fēng),因此將這5個氣象觀測站過去10年(其中2套交通氣象站為5年)的歷史觀測資料,通過最小二乘法建立起10 min平均風(fēng)速與該站極大風(fēng)之間的關(guān)系[29],得到5套氣象觀測站的極風(fēng)系數(shù)(表2)。再由WT輸出的10 min平均風(fēng)按照極風(fēng)系數(shù)預(yù)測某一點或路段的極大風(fēng)可能值。

    表2 張家口、崇禮區(qū)域5個觀測站2—5月極風(fēng)系數(shù)

    4 結(jié)論與討論

    (1)分析了2001—2020年河北省142個國家氣象站風(fēng)的特征,發(fā)現(xiàn)張家口西北部及承德西部是6級及以上大風(fēng)日數(shù)最多的區(qū)域。利用EOF正交函數(shù)分解各站逐日最大風(fēng)速距平場并分析發(fā)現(xiàn),在過去20年中日最大風(fēng)速呈現(xiàn)明顯減少的趨勢,但是張家口、承德北部地區(qū)在第二、第三模態(tài)中略有增大的趨勢,這說明這些區(qū)域的大風(fēng)還是需要格外關(guān)注。

    (2)應(yīng)用WT模型對張承高速公路沿線2個區(qū)域的30個代表日的風(fēng)進行模擬,并運用區(qū)域內(nèi)5套氣象站的逐分鐘觀測數(shù)據(jù)進行檢驗。發(fā)現(xiàn)WT模型模擬輸出的10 min平均風(fēng)速與實際觀測值相關(guān)系數(shù)達(dá)到0.6225并通過0.001的顯著性水平檢驗。各站模擬的結(jié)果與實況之間的誤差在±1 m/s之間的均占40%以上,±2 m/s之間的占比基本能夠達(dá)到80%以上。

    (3)本文提供了一種對高速公路沿線山區(qū)復(fù)雜地形風(fēng)的數(shù)值預(yù)報模擬和訂正的可行方法。在高速公路精細(xì)化氣象服務(wù)中,可應(yīng)用本地地形資料以及時空分辨率更高的智能網(wǎng)格預(yù)報產(chǎn)品或其他數(shù)值預(yù)報模式輸出的風(fēng)預(yù)報產(chǎn)品作為驅(qū)動WT模型的數(shù)據(jù)源,進而對特殊地形的風(fēng)速進行訂正。另外,可根據(jù)區(qū)域內(nèi)代表站的極風(fēng)系數(shù)進行極大風(fēng)的預(yù)報。

    (4)復(fù)雜地形的風(fēng)的預(yù)報難度較大,今后需要進一步積累觀測數(shù)據(jù),不斷優(yōu)化WT模型的參數(shù),以提高這些區(qū)域高速公路沿線風(fēng)的預(yù)報和訂正能力。

    猜你喜歡
    崇禮張家口大風(fēng)
    中國農(nóng)業(yè)發(fā)展銀行張家口分行
    崇禮當(dāng)代詞典
    海峽旅游(2019年6期)2019-10-30 02:09:08
    大風(fēng)吹(二)
    幼兒100(2017年31期)2017-11-27 02:37:46
    大風(fēng)吹(一)
    幼兒100(2017年28期)2017-10-27 01:45:49
    人小鬼大狄仁杰
    加大張家口草原生態(tài)建設(shè)支持力度
    行走在張家口的熱土上
    韓國:崇禮門
    可再生能源建筑采暖研討會2016年1月崇禮有約
    太陽能(2015年11期)2015-04-10 12:53:04
    北京至張家口崇禮擬建新高速
    久久精品亚洲熟妇少妇任你| 国产精品自产拍在线观看55亚洲| 国产欧美日韩一区二区三区在线| 亚洲精品一卡2卡三卡4卡5卡| 天堂影院成人在线观看| 老司机午夜十八禁免费视频| 丝袜人妻中文字幕| 琪琪午夜伦伦电影理论片6080| 一本大道久久a久久精品| 电影成人av| 午夜影院日韩av| 99精品在免费线老司机午夜| 亚洲七黄色美女视频| 色哟哟哟哟哟哟| www国产在线视频色| 欧美激情 高清一区二区三区| 一级毛片精品| 最近最新中文字幕大全免费视频| 亚洲成人久久性| 女人被狂操c到高潮| 人人妻人人爽人人添夜夜欢视频| 亚洲一区二区三区不卡视频| 这个男人来自地球电影免费观看| 美女大奶头视频| 精品国产超薄肉色丝袜足j| 搡老熟女国产l中国老女人| av天堂久久9| 久久久国产成人精品二区| 亚洲人成电影免费在线| 亚洲aⅴ乱码一区二区在线播放 | 欧美在线一区亚洲| 啦啦啦 在线观看视频| 男女下面进入的视频免费午夜 | 日本三级黄在线观看| 美女 人体艺术 gogo| 亚洲人成电影观看| 99国产精品一区二区蜜桃av| 美女高潮到喷水免费观看| 1024香蕉在线观看| 97人妻精品一区二区三区麻豆 | 亚洲第一电影网av| 日本撒尿小便嘘嘘汇集6| 如日韩欧美国产精品一区二区三区| 国产亚洲精品久久久久5区| 日本精品一区二区三区蜜桃| 日韩 欧美 亚洲 中文字幕| 午夜视频精品福利| 在线观看免费视频日本深夜| 在线十欧美十亚洲十日本专区| 成人欧美大片| 多毛熟女@视频| 桃色一区二区三区在线观看| 日韩有码中文字幕| av超薄肉色丝袜交足视频| 俄罗斯特黄特色一大片| 国产成年人精品一区二区| 国产又爽黄色视频| 老熟妇仑乱视频hdxx| 一二三四在线观看免费中文在| 久久精品国产综合久久久| 50天的宝宝边吃奶边哭怎么回事| 黄色丝袜av网址大全| 欧美一级毛片孕妇| 国产麻豆69| 好男人在线观看高清免费视频 | 久久久久久久午夜电影| 欧美午夜高清在线| 亚洲第一av免费看| 自线自在国产av| 国产精品亚洲一级av第二区| 欧美日本视频| 欧美亚洲日本最大视频资源| 美女国产高潮福利片在线看| 男男h啪啪无遮挡| 国产精品一区二区免费欧美| 97超级碰碰碰精品色视频在线观看| www.www免费av| 国产一区二区激情短视频| 最好的美女福利视频网| 18美女黄网站色大片免费观看| 亚洲自拍偷在线| 午夜精品国产一区二区电影| 亚洲三区欧美一区| 亚洲五月色婷婷综合| 欧美不卡视频在线免费观看 | 久久精品国产亚洲av香蕉五月| 久久国产乱子伦精品免费另类| 此物有八面人人有两片| 成人国产一区最新在线观看| 好男人电影高清在线观看| 成人av一区二区三区在线看| 老司机靠b影院| 日韩中文字幕欧美一区二区| 午夜两性在线视频| 亚洲人成网站在线播放欧美日韩| 国语自产精品视频在线第100页| 不卡av一区二区三区| 香蕉久久夜色| 免费在线观看视频国产中文字幕亚洲| 看片在线看免费视频| 精品国内亚洲2022精品成人| 真人做人爱边吃奶动态| 热re99久久国产66热| 亚洲 欧美一区二区三区| 国产精品电影一区二区三区| 脱女人内裤的视频| 一级毛片精品| 亚洲国产欧美网| 国产亚洲精品综合一区在线观看 | 男女做爰动态图高潮gif福利片 | 性色av乱码一区二区三区2| 中文字幕另类日韩欧美亚洲嫩草| 中文字幕另类日韩欧美亚洲嫩草| 性色av乱码一区二区三区2| 欧美日韩乱码在线| 亚洲精品中文字幕一二三四区| 亚洲无线在线观看| 变态另类丝袜制服| 免费高清在线观看日韩| 久久久精品欧美日韩精品| 国产精品自产拍在线观看55亚洲| 日本免费a在线| 亚洲av熟女| 久久久久久久久久久久大奶| 1024视频免费在线观看| 亚洲精品久久国产高清桃花| 日韩高清综合在线| 91国产中文字幕| 如日韩欧美国产精品一区二区三区| 欧美人与性动交α欧美精品济南到| 亚洲一区二区三区色噜噜| 久久久久久久久免费视频了| 欧美一级a爱片免费观看看 | 日本欧美视频一区| 久久狼人影院| 亚洲精品久久成人aⅴ小说| 一区二区日韩欧美中文字幕| 久久婷婷人人爽人人干人人爱 | 成人三级做爰电影| 亚洲男人的天堂狠狠| 国产人伦9x9x在线观看| 九色国产91popny在线| 怎么达到女性高潮| 亚洲精品av麻豆狂野| 一级毛片高清免费大全| 91在线观看av| 欧美日韩精品网址| 国产99白浆流出| 亚洲成人久久性| 人人妻人人澡欧美一区二区 | 亚洲人成电影免费在线| 老司机午夜十八禁免费视频| 两人在一起打扑克的视频| 亚洲av电影不卡..在线观看| www.熟女人妻精品国产| 91字幕亚洲| 国产成人欧美| 人人妻人人澡欧美一区二区 | bbb黄色大片| av片东京热男人的天堂| 欧美成人免费av一区二区三区| 校园春色视频在线观看| 国产亚洲欧美精品永久| xxx96com| 精品久久久精品久久久| 色老头精品视频在线观看| 国产精品亚洲一级av第二区| 黄色 视频免费看| 波多野结衣高清无吗| 搡老岳熟女国产| 欧美激情 高清一区二区三区| 国产一区二区在线av高清观看| 国产精品自产拍在线观看55亚洲| 国产成人精品无人区| 久久精品成人免费网站| 国产成人av教育| 成人国语在线视频| 久久久久国产一级毛片高清牌| 熟女少妇亚洲综合色aaa.| 两个人看的免费小视频| 在线观看舔阴道视频| 曰老女人黄片| 日韩中文字幕欧美一区二区| 18禁观看日本| 一级a爱片免费观看的视频| 国产精品,欧美在线| 国产精品日韩av在线免费观看 | 亚洲精品av麻豆狂野| 亚洲av成人一区二区三| 欧美最新免费一区二区三区| 22中文网久久字幕| 亚洲精品日韩av片在线观看| 欧美黑人欧美精品刺激| 琪琪午夜伦伦电影理论片6080| 欧美成人性av电影在线观看| 日韩av在线大香蕉| 中国美女看黄片| 国产精品美女特级片免费视频播放器| 亚洲一区高清亚洲精品| 99久久成人亚洲精品观看| 赤兔流量卡办理| 伦理电影大哥的女人| 老熟妇仑乱视频hdxx| 国产不卡一卡二| 波多野结衣高清作品| 最近中文字幕高清免费大全6 | 校园春色视频在线观看| 欧美潮喷喷水| 少妇熟女aⅴ在线视频| 午夜福利在线观看免费完整高清在 | 欧美+亚洲+日韩+国产| 黄片wwwwww| 欧美日韩中文字幕国产精品一区二区三区| 69人妻影院| 1024手机看黄色片| 99久国产av精品| 男女那种视频在线观看| 小蜜桃在线观看免费完整版高清| 欧美不卡视频在线免费观看| 亚洲成人久久性| www日本黄色视频网| 国产黄a三级三级三级人| 女人十人毛片免费观看3o分钟| 久99久视频精品免费| 精品久久久久久成人av| 国产精品电影一区二区三区| 伦理电影大哥的女人| 大型黄色视频在线免费观看| 男女下面进入的视频免费午夜| 国产国拍精品亚洲av在线观看| 俄罗斯特黄特色一大片| 亚洲最大成人手机在线| 免费看av在线观看网站| 欧美一级a爱片免费观看看| 国产成人一区二区在线| 夜夜看夜夜爽夜夜摸| 美女免费视频网站| 欧美黑人欧美精品刺激| 91久久精品国产一区二区成人| 亚洲va在线va天堂va国产| 国产日本99.免费观看| 观看免费一级毛片| 不卡视频在线观看欧美| av.在线天堂| 午夜福利欧美成人| 亚洲av成人av| 日本a在线网址| 毛片女人毛片| 欧美成人免费av一区二区三区| 中文字幕av在线有码专区| 国产黄片美女视频| 人人妻,人人澡人人爽秒播| 黄色女人牲交| 九九热线精品视视频播放| 欧美xxxx黑人xx丫x性爽| 欧美色欧美亚洲另类二区| 老师上课跳d突然被开到最大视频| 中文资源天堂在线| 91久久精品电影网| 免费av毛片视频| 嫩草影院新地址| 精品欧美国产一区二区三| 欧洲精品卡2卡3卡4卡5卡区| 身体一侧抽搐| 亚洲精品一区av在线观看| 亚洲av成人av| 最后的刺客免费高清国语| 亚洲成a人片在线一区二区| 国产探花在线观看一区二区| 色综合色国产| 日韩人妻高清精品专区| 国产精品av视频在线免费观看| 噜噜噜噜噜久久久久久91| 两个人视频免费观看高清| 精品无人区乱码1区二区| 久久久久精品国产欧美久久久| 老司机深夜福利视频在线观看| 91久久精品电影网| 国产成年人精品一区二区| av视频在线观看入口| 日韩亚洲欧美综合| 欧美色视频一区免费| 级片在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 内射极品少妇av片p| 波野结衣二区三区在线| 身体一侧抽搐| 亚洲乱码一区二区免费版| 国产精品一区二区三区四区久久| 欧美成人a在线观看| 国内精品久久久久久久电影| 精品久久久噜噜| 国产高清不卡午夜福利| 亚洲精品日韩av片在线观看| 嫩草影院精品99| 国产精品一区www在线观看 | ponron亚洲| 亚洲熟妇中文字幕五十中出| 日韩欧美精品免费久久| 欧美国产日韩亚洲一区| 伦理电影大哥的女人| 村上凉子中文字幕在线| 国产国拍精品亚洲av在线观看| 给我免费播放毛片高清在线观看| 国产精品久久久久久久久免| 麻豆一二三区av精品| 不卡一级毛片| 日韩欧美一区二区三区在线观看| 亚洲狠狠婷婷综合久久图片| 久99久视频精品免费| 国产精品电影一区二区三区| 亚洲中文字幕日韩| 国内揄拍国产精品人妻在线| 少妇裸体淫交视频免费看高清| 很黄的视频免费| 国产伦精品一区二区三区四那| 亚洲国产高清在线一区二区三| 日韩精品有码人妻一区| а√天堂www在线а√下载| 成人鲁丝片一二三区免费| 亚洲无线观看免费| av中文乱码字幕在线| 久久久久久久午夜电影| 久99久视频精品免费| av专区在线播放| 一进一出抽搐gif免费好疼| 天堂动漫精品| 男插女下体视频免费在线播放| 亚洲图色成人| 少妇高潮的动态图| eeuss影院久久| 日本a在线网址| 如何舔出高潮| 少妇高潮的动态图| 日韩精品有码人妻一区| 啪啪无遮挡十八禁网站| 欧美极品一区二区三区四区| 黄色配什么色好看| 欧美bdsm另类| 丝袜美腿在线中文| 日韩欧美在线二视频| 一个人看的www免费观看视频| 日韩精品有码人妻一区| 免费人成视频x8x8入口观看| 蜜桃亚洲精品一区二区三区| 搡老熟女国产l中国老女人| 日本五十路高清| 亚洲性夜色夜夜综合| 最新在线观看一区二区三区| 国产精品亚洲美女久久久| 国产一区二区三区av在线 | 精品久久久噜噜| 国产一区二区三区视频了| 亚洲一级一片aⅴ在线观看| 一边摸一边抽搐一进一小说| 久久久久精品国产欧美久久久| 国产老妇女一区| 久久久久久久久久黄片| 亚洲最大成人手机在线| 大又大粗又爽又黄少妇毛片口| 久久久久久九九精品二区国产| 国产成人福利小说| 美女 人体艺术 gogo| 赤兔流量卡办理| 午夜福利在线在线| 极品教师在线免费播放| 色视频www国产| 中国美白少妇内射xxxbb| 波多野结衣高清无吗| 国产精品女同一区二区软件 | 欧美日韩黄片免| 国产精品三级大全| 中文字幕精品亚洲无线码一区| 久久久久久久久中文| 国产三级中文精品| 丰满乱子伦码专区| 性欧美人与动物交配| 看免费成人av毛片| 精品久久久久久久久亚洲 | 在线国产一区二区在线| 人人妻人人澡欧美一区二区| 在线观看美女被高潮喷水网站| 男人狂女人下面高潮的视频| 国内精品一区二区在线观看| eeuss影院久久| 亚洲精品影视一区二区三区av| av黄色大香蕉| 国产老妇女一区| 床上黄色一级片| 51国产日韩欧美| 国产成人一区二区在线| 此物有八面人人有两片| 国产精品一区二区三区四区久久| 最近中文字幕高清免费大全6 | 国产在线男女| 村上凉子中文字幕在线| 亚洲精华国产精华精| 男女那种视频在线观看| 老熟妇仑乱视频hdxx| 九九热线精品视视频播放| 国内精品久久久久久久电影| 特级一级黄色大片| 男女做爰动态图高潮gif福利片| 最后的刺客免费高清国语| 欧美潮喷喷水| 亚洲美女黄片视频| 嫩草影院新地址| 国产在视频线在精品| 身体一侧抽搐| 超碰av人人做人人爽久久| 如何舔出高潮| 亚洲成人免费电影在线观看| 变态另类丝袜制服| 色视频www国产| av福利片在线观看| 国产精品人妻久久久影院| av女优亚洲男人天堂| 久久久久免费精品人妻一区二区| 色噜噜av男人的天堂激情| 久久久久国内视频| a级一级毛片免费在线观看| 久久人人精品亚洲av| 久久久久久久精品吃奶| 真实男女啪啪啪动态图| 午夜福利在线在线| 狠狠狠狠99中文字幕| 成人美女网站在线观看视频| 久久久午夜欧美精品| 69人妻影院| 淫妇啪啪啪对白视频| av在线观看视频网站免费| 狠狠狠狠99中文字幕| 男女那种视频在线观看| 神马国产精品三级电影在线观看| 久久6这里有精品| 99热这里只有是精品50| 午夜影院日韩av| 日韩欧美精品免费久久| 99久久精品一区二区三区| 99国产极品粉嫩在线观看| 欧美黑人巨大hd| 99精品在免费线老司机午夜| 国产精品无大码| 国产成人福利小说| 一边摸一边抽搐一进一小说| 国产精品野战在线观看| 日本与韩国留学比较| 夜夜爽天天搞| 日日夜夜操网爽| 亚洲精品日韩av片在线观看| 午夜福利视频1000在线观看| 国产在线男女| 婷婷六月久久综合丁香| 国产成人aa在线观看| 伊人久久精品亚洲午夜| 国产欧美日韩精品亚洲av| 久久草成人影院| 日本黄大片高清| 国产精品精品国产色婷婷| 乱码一卡2卡4卡精品| 色在线成人网| 嫁个100分男人电影在线观看| 亚洲中文日韩欧美视频| 黄色女人牲交| 三级毛片av免费| av福利片在线观看| 国产精品三级大全| 高清在线国产一区| 欧美一区二区精品小视频在线| 欧美一区二区亚洲| 亚洲第一电影网av| 97碰自拍视频| 99久久无色码亚洲精品果冻| 中文字幕av在线有码专区| 午夜a级毛片| 亚洲人成伊人成综合网2020| 波野结衣二区三区在线| 淫秽高清视频在线观看| 精品99又大又爽又粗少妇毛片 | 夜夜爽天天搞| 天堂av国产一区二区熟女人妻| 九九久久精品国产亚洲av麻豆| 别揉我奶头~嗯~啊~动态视频| 国产白丝娇喘喷水9色精品| 亚洲国产欧美人成| 久久久久精品国产欧美久久久| 午夜福利在线在线| 春色校园在线视频观看| 国产免费男女视频| 亚洲内射少妇av| 亚洲男人的天堂狠狠| 人人妻,人人澡人人爽秒播| 午夜福利在线观看吧| 亚洲无线观看免费| 亚洲自偷自拍三级| 中文字幕免费在线视频6| 国产爱豆传媒在线观看| 亚洲欧美清纯卡通| 熟妇人妻久久中文字幕3abv| 国产国拍精品亚洲av在线观看| 亚洲久久久久久中文字幕| 婷婷精品国产亚洲av在线| 午夜免费激情av| 日本熟妇午夜| a级毛片免费高清观看在线播放| 乱码一卡2卡4卡精品| 国产 一区精品| 别揉我奶头~嗯~啊~动态视频| 三级男女做爰猛烈吃奶摸视频| 美女xxoo啪啪120秒动态图| 最近最新中文字幕大全电影3| 国产午夜精品久久久久久一区二区三区 | 国产精品电影一区二区三区| 国产免费av片在线观看野外av| 精品无人区乱码1区二区| 两个人的视频大全免费| 69人妻影院| 国产亚洲精品综合一区在线观看| 国产精品国产三级国产av玫瑰| 91精品国产九色| 欧美+亚洲+日韩+国产| 久久九九热精品免费| 黄色配什么色好看| or卡值多少钱| 深夜a级毛片| av女优亚洲男人天堂| 在线a可以看的网站| 深爱激情五月婷婷| 精品国内亚洲2022精品成人| www.色视频.com| 99久国产av精品| 日韩精品青青久久久久久| 99久久精品国产国产毛片| 欧美不卡视频在线免费观看| 狂野欧美激情性xxxx在线观看| 亚洲在线观看片| 搡老岳熟女国产| 97超级碰碰碰精品色视频在线观看| 麻豆久久精品国产亚洲av| 国产精品久久久久久av不卡| 一个人看视频在线观看www免费| 欧美一区二区国产精品久久精品| 国产色爽女视频免费观看| 日韩精品有码人妻一区| 极品教师在线视频| 久久欧美精品欧美久久欧美| 男人的好看免费观看在线视频| 国产精品一区二区三区四区免费观看 | 国产精品99久久久久久久久| 精品久久久久久久久av| 啦啦啦观看免费观看视频高清| 精品乱码久久久久久99久播| 精品人妻一区二区三区麻豆 | 国产亚洲精品综合一区在线观看| 欧美精品国产亚洲| 久久精品国产鲁丝片午夜精品 | 亚洲人成伊人成综合网2020| 成人毛片a级毛片在线播放| 国内精品宾馆在线| 全区人妻精品视频| 国模一区二区三区四区视频| 成年版毛片免费区| 欧美zozozo另类| 国产成人一区二区在线| 亚洲内射少妇av| 啦啦啦韩国在线观看视频| 97人妻精品一区二区三区麻豆| 成人鲁丝片一二三区免费| 日日啪夜夜撸| 国产美女午夜福利| 亚洲一区二区三区色噜噜| 久久国内精品自在自线图片| 精品不卡国产一区二区三区| 免费人成视频x8x8入口观看| 伊人久久精品亚洲午夜| 日韩欧美免费精品| 窝窝影院91人妻| 亚洲av成人av| 老师上课跳d突然被开到最大视频| www日本黄色视频网| 蜜桃亚洲精品一区二区三区| 亚洲美女黄片视频| 99热这里只有是精品在线观看| 大型黄色视频在线免费观看| 亚洲国产色片| 日本-黄色视频高清免费观看| videossex国产| 联通29元200g的流量卡| 日本精品一区二区三区蜜桃| 精品一区二区三区av网在线观看| 在线免费观看的www视频| 日韩欧美国产在线观看| 丰满的人妻完整版| 尤物成人国产欧美一区二区三区| 久久精品国产自在天天线| 老司机福利观看| 搡老岳熟女国产| 女的被弄到高潮叫床怎么办 | 亚洲av免费高清在线观看| 亚洲中文字幕一区二区三区有码在线看| 男插女下体视频免费在线播放| 在线观看免费视频日本深夜| 欧美成人一区二区免费高清观看| 国产乱人伦免费视频| 欧美成人a在线观看| 在线观看66精品国产| 久久亚洲真实| 国产精品1区2区在线观看.| 中国美女看黄片| 欧美三级亚洲精品| 久久九九热精品免费| 亚洲中文字幕一区二区三区有码在线看| 免费观看的影片在线观看| 一区二区三区激情视频| 午夜精品在线福利|