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

    觀測(cè)間隔對(duì)極值風(fēng)速的影響研究

    2021-08-27 07:57:32肖鈺川
    工程力學(xué) 2021年8期
    關(guān)鍵詞:原始數(shù)據(jù)風(fēng)向極值

    肖鈺川,全 涌,顧 明

    (同濟(jì)大學(xué)土木工程防災(zāi)國家重點(diǎn)實(shí)驗(yàn)室,上海200092)

    在風(fēng)氣候觀測(cè)數(shù)據(jù)格式方面,不同的國家存在很多的差異,其中在風(fēng)速樣本觀測(cè)間隔的選取上尤為突出。歷史上,在觀測(cè)風(fēng)速時(shí)所采用的觀測(cè)設(shè)備大都是手動(dòng)的,故人們?cè)谟涗洈?shù)據(jù)時(shí)總是隔一段時(shí)間去記錄一次,這樣就導(dǎo)致了不同地區(qū)數(shù)據(jù)記錄間隔的差異。各國的規(guī)范中觀測(cè)間隔多種多樣。以中國的氣象站[1]為例:其記錄數(shù)據(jù)的方式有3種:第1種是每隔6 h 記錄一次風(fēng)速和風(fēng)向,記錄的數(shù)據(jù)是4個(gè)整點(diǎn)(0點(diǎn)、6點(diǎn)、12點(diǎn)和18點(diǎn))前最后2 min 的平均風(fēng)速和對(duì)應(yīng)風(fēng)向;第2種是記錄每天24 h 整點(diǎn)前最后2 min 的平均風(fēng)速和對(duì)應(yīng)風(fēng)向;第3種是連續(xù)觀測(cè),每10 min 取一次平均風(fēng)速和平均風(fēng)向。然后分別在一天的若干個(gè)值中取其日最大值作為日極值??梢姡m然同樣是日極值,但其觀測(cè)間隔不同,這必定會(huì)導(dǎo)致不同重現(xiàn)期下極值風(fēng)速預(yù)測(cè)的差異。全涌等[2]研究了平均時(shí)距對(duì)良態(tài)風(fēng)氣候地區(qū)極值風(fēng)速的影響,給出了不同平均時(shí)距極值風(fēng)速的轉(zhuǎn)換關(guān)系式,解決了風(fēng)速格式不統(tǒng)一的部分問題,但對(duì)風(fēng)速格式中觀測(cè)間隔不同的問題未展開研究。Mike等[3]研究了不同觀測(cè)間隔下風(fēng)速極值分布以及母體分布的區(qū)別,提出對(duì)于風(fēng)速的母體分布,二者相對(duì)較符合;對(duì)于風(fēng)速的極值分布,二者差異隨重現(xiàn)期不同而不同,但未提出修正方法。

    第1種觀測(cè)格式常見于建國初期的風(fēng)氣候觀測(cè)記錄文件;第2種觀測(cè)格式常見于我國七八十年代的風(fēng)氣候記錄文件;第3種則常見于近期的風(fēng)氣候記錄文件。毫無疑問,第3種無論從樣本的廣度還是精度上均好于前2種。而前2種方法均有缺陷:只記錄了每小時(shí)或每6 h 最后2 min 的平均值,卻很容易忽略當(dāng)天內(nèi)其他時(shí)段發(fā)生更大的風(fēng)速值。同時(shí),前2種采樣方式由于自身的觀測(cè)間隔亦不相同,間隔小的采樣方式其樣本量相對(duì)更充分,這對(duì)結(jié)果的估計(jì)必然產(chǎn)生影響。

    然而,這種由于觀測(cè)間隔不同而對(duì)風(fēng)速極值的概率統(tǒng)計(jì)特性及估計(jì)產(chǎn)生的影響尚未被研究,原因很可能是缺少一種完全充分且合理的采樣方法。如何利用上述前2種采樣方式得到的數(shù)據(jù)樣本進(jìn)行極值風(fēng)速預(yù)測(cè),使得結(jié)果接近第3種采樣方式所得的結(jié)果是一個(gè)很有意義的課題。

    而美國的氣象臺(tái)最近公布了一個(gè)相對(duì)新的數(shù)據(jù)集,他提供了從2000年?2016年美國11個(gè)地區(qū)、每分鐘的風(fēng)速值和風(fēng)向角,這大大增加了數(shù)據(jù)的完整性,基于此做出的極值估計(jì)結(jié)果在一定程度上就是風(fēng)速極值的真實(shí)值。因此,可以通過對(duì)比真實(shí)值和基于前述3種觀測(cè)間隔下估計(jì)的風(fēng)速極值之間的差異,將其量化并研究,找到修正方法或找出最合適的概率估計(jì)方法。進(jìn)而使我國的工程設(shè)計(jì)更準(zhǔn)確。

    本文基于互聯(lián)網(wǎng)上獲得的美國良態(tài)風(fēng)氣候地區(qū)11個(gè)氣象站的風(fēng)速風(fēng)向長期連續(xù)觀測(cè)數(shù)據(jù)[4?5],先就數(shù)據(jù)序列的觀測(cè)間隔對(duì)極值風(fēng)速預(yù)測(cè)值的差異進(jìn)行了研究,即分析以觀測(cè)間隔為6 h、1 h 和連續(xù)觀測(cè)的2 min 平均風(fēng)速樣本為基礎(chǔ)的極值風(fēng)速預(yù)測(cè)結(jié)果的差異,從原始數(shù)據(jù)入手找出前2種方式向第3種方式轉(zhuǎn)換的修正方法。所得方法的結(jié)果更加精確和符合實(shí)際,有利于工程推廣。然后可根據(jù)平均時(shí)距導(dǎo)致極值風(fēng)速差異的修正方法[2],將連續(xù)2 min 觀測(cè)間隔對(duì)應(yīng)的結(jié)果轉(zhuǎn)換為連續(xù)10 min觀測(cè)間隔對(duì)應(yīng)的結(jié)果,這樣可以有效的解決風(fēng)速風(fēng)向觀測(cè)數(shù)據(jù)格式(觀測(cè)間隔和平均時(shí)距兩個(gè)方面)對(duì)極值風(fēng)速預(yù)測(cè)的影響問題。

    1 風(fēng)速樣本的選取與預(yù)處理

    美國國家海洋和大氣管理局(ASOS, NOAA)近期公布了一個(gè)相對(duì)新的數(shù)據(jù)集[3?4],從中隨機(jī)選取了美國11座城市的從2000年?2016年這16年間每分鐘平均風(fēng)速數(shù)據(jù)。本文將原始數(shù)據(jù)分別按如下方式處理:1)取每6 h 最后2 min 數(shù)據(jù)的平均風(fēng)速和平均風(fēng)向,從一天中的4個(gè)值中選取最大風(fēng)速值和其對(duì)應(yīng)的風(fēng)向做為日極值風(fēng)速和風(fēng)向;2)取每小時(shí)最后2 min 數(shù)據(jù)的平均風(fēng)速和平均風(fēng)向,從一天中的24個(gè)值中選取最大風(fēng)速值和其對(duì)應(yīng)的風(fēng)向做為日極值風(fēng)速和風(fēng)向;3)取連續(xù)2 min數(shù)據(jù)的平均風(fēng)速和平均風(fēng)向,從一天中的24×30=720個(gè)值中選取最大風(fēng)速值和其對(duì)應(yīng)的風(fēng)向作為日極值風(fēng)速和風(fēng)向;將上述3種觀測(cè)間隔下得到的風(fēng)速日極值樣本進(jìn)行分析,并進(jìn)行全時(shí)程檢驗(yàn),圖1 展示了美國中部城市托皮卡市(Topeka,簡稱TOP)2000年1月1 日?2000年1月7日的日極值風(fēng)速樣本。

    由圖1可以看出,連續(xù)2 min 觀測(cè)樣本所得到的日極值一定不小于每小時(shí)最后2 min 觀測(cè)樣本所得到的日極值,同時(shí)每小時(shí)最后2 min 觀測(cè)樣本所得到的日極值一定大于等于每6 h 最后2 min 觀測(cè)樣本所得到的日極值,符合實(shí)際情況。

    圖1 TOP市不同觀測(cè)間隔日極值對(duì)比Fig.1 Comparison of daily extremevaluesof different observing intervalsin TOP

    2 不同觀測(cè)間隔對(duì)應(yīng)的極值擬合與預(yù)測(cè)

    Harris[6?9]給出了一套相對(duì)完整的基于觀測(cè)數(shù)據(jù)的極值估算方法。與過去的所有極值分析方法相比,這套理論更加著眼于極值的特性,并為其量身打造了相應(yīng)的擬合方法,盡可能地削減了極值分析過程中可能造成的誤差?;诖?,本文采用了改進(jìn)后的獨(dú)立風(fēng)暴法[10]:考慮左截?cái)?,并?duì)簡化變量的計(jì)算方法進(jìn)行改進(jìn)[11],同時(shí)采用加權(quán)最小二乘擬合極值Ⅰ型分布,以使得獨(dú)立風(fēng)暴法得到的擬合結(jié)果更符合實(shí)際。將上述3種觀測(cè)間隔下的日極值通過Harris[6? 9]的極值估算方法,可得到不同重現(xiàn)期的極值估算結(jié)果,如圖2所示。

    圖2可以看出,雖然各城市間存在差異,但3種觀測(cè)間隔的預(yù)測(cè)結(jié)果均不相同。連續(xù)2 min 的采樣方式對(duì)應(yīng)的結(jié)果要大于1 h 最后2 min 的采樣方式的結(jié)果,而1 h 最后2 min 的采樣方式的結(jié)果要大于6 h 最后2 min 的采樣方式的結(jié)果,這剛好驗(yàn)證了前文的猜想。

    圖2 各城市不同觀測(cè)間隔極值風(fēng)速預(yù)測(cè)Fig.2 Prediction of extreme wind speed at different observing intervals in different cities

    3 觀測(cè)間隔為1 h 的數(shù)據(jù)樣本的修正方法

    由每分鐘連續(xù)的原始數(shù)據(jù)得到的每2 min 連續(xù)的原始數(shù)據(jù)和每小時(shí)最后2 min 原始數(shù)據(jù)。用這2種數(shù)據(jù)計(jì)算日極值,共分別得到6000個(gè)數(shù)據(jù),將6000個(gè)數(shù)據(jù)做差,得到約6000個(gè)差值。通過分析這些差值的特性,并將這些差值全部加某一定值C,用瑞利分布擬合,可以得到很好的結(jié)果。通過枚舉0 m/s~3 m/s的含有一位小數(shù)的數(shù)字,并進(jìn)行擬合優(yōu)度檢驗(yàn)(下同),可得出當(dāng)C取0.5 m/s時(shí)擬合效果最佳。即:

    日極值差加0.5后的瑞利概率密度分布如圖3。圖3展示了11座城市的連續(xù)2 min 所得日風(fēng)速極值和每小時(shí)最后2 min 的日風(fēng)速極值差值整體加0.5 m/s后的概率分布情況,差值按區(qū)間分成20份并畫出柱狀圖,并用瑞利分布擬合。所得相關(guān)系數(shù)如圖所示,且均大于0.9,這說明瑞利分布擬合效果較好。

    圖3 各城市不同觀測(cè)間隔風(fēng)速差值擬合Fig.3 Fitting of wind speed difference between different observing intervals in different cities

    表1各城市瑞利分布參數(shù)Table1 Rayleigh distribution parametersof each city

    當(dāng)由每小時(shí)最后2 min 日極值風(fēng)速數(shù)據(jù)轉(zhuǎn)換為連續(xù)2 min 日極值風(fēng)速數(shù)據(jù)時(shí),只需將前者的數(shù)據(jù)分別加上生成的隨機(jī)數(shù)再減去0.5 m/s即可。圖4為經(jīng)此修正后各城市的極值風(fēng)速預(yù)測(cè)圖。

    由圖4可以看出,1 h 最后2 min 觀測(cè)數(shù)據(jù)所對(duì)應(yīng)的日極值風(fēng)速與連續(xù)2 min 觀測(cè)數(shù)據(jù)所對(duì)應(yīng)的日極值風(fēng)速之間的差異明顯減小,尤其是當(dāng)簡化變量Y=?ln(?ln(P))<0時(shí),即獨(dú)立風(fēng)暴法預(yù)測(cè)數(shù)據(jù)的下尾部部分,幾乎是重合的,而上尾部因?yàn)閿?shù)據(jù)量較少,且離散程度大,故仍可能存在一定差異。這可以說明通過對(duì)原始數(shù)據(jù)的修正方法是有效的。表2、表3展示了各城市在10年和50年重現(xiàn)期下極值風(fēng)速值的計(jì)算結(jié)果及誤差。

    圖4 各城市修正后不同觀測(cè)間隔極值風(fēng)速預(yù)測(cè)(1 h)Fig.4 Prediction of extreme wind speed in different observing intervals after correction in each city (per hour)

    由表2和表3可以看出,對(duì)于大部分城市而言,修正后的1 h 最后2 min 的預(yù)測(cè)曲線和連續(xù)2 min的預(yù)測(cè)曲線十分接近。對(duì)于個(gè)別城市而言,二者仍會(huì)存在一定誤差,這是由于本文所采用 σ的值為統(tǒng)計(jì)得到的均值所致,但即使是差別最大的城市,在50年重現(xiàn)期下,二者的絕對(duì)誤差在3%~10.2%,且該誤差隨重現(xiàn)期的增加而減小,使修正后的1 h最后2 min 的預(yù)測(cè)結(jié)果更加接近連續(xù)2 min 的預(yù)測(cè)結(jié)果。

    表2 每小時(shí)最后2 min 觀測(cè)間隔在10年重現(xiàn)期下修正前后風(fēng)速極值對(duì)比Table 2 Comparison of extreme wind speeds before and after correction under 10-year return period of 2-minute per hour

    表3 每小時(shí)最后2 min 觀測(cè)間隔50年重現(xiàn)期下修正前后風(fēng)速極值數(shù)據(jù)對(duì)比Table 3 Comparison of extreme wind speeds before and after correction under 50-year return period of 2-minute per hour

    4 觀測(cè)間隔為6 h 的數(shù)據(jù)樣本的修正方法

    與第3節(jié)方法類似,由每分鐘連續(xù)的原始數(shù)據(jù)得到的連續(xù)2 min 的原始數(shù)據(jù)和每6 h 最后2 min原始數(shù)據(jù)。用這2種數(shù)據(jù)計(jì)算日極值,共分別得到6000個(gè)數(shù)據(jù),將6000個(gè)數(shù)據(jù)做差,得到約6000個(gè)差值。通過分析這些差值的特性,并將這些差值全部加某一定值C,用瑞利分布擬合,可以得到很好的結(jié)果。經(jīng)過計(jì)算C取0.2 m/s時(shí)擬合效果最佳。即:

    由圖5可以看出,每6 h 最后2 min 觀測(cè)數(shù)據(jù)所對(duì)應(yīng)的日極值風(fēng)速與連續(xù)2 min 觀測(cè)數(shù)據(jù)所對(duì)應(yīng)的日極值風(fēng)速之間的差異明顯被縮小,其10年重現(xiàn)期和50年重現(xiàn)期下極值風(fēng)速值的預(yù)測(cè)結(jié)果見表4和表5。

    圖5 各城市修正后不同觀測(cè)間隔極值風(fēng)速預(yù)測(cè)(6 h)Fig.5 Prediction of extreme wind speed in different observing intervalsafter correction in each city (per six hours)

    由表4和表5可以看出,對(duì)于大部分城市而言,修正后的6 h 最后2 min 的預(yù)測(cè)曲線和連續(xù)2 min的預(yù)測(cè)曲線十分接近。對(duì)于個(gè)別城市而言,二者仍會(huì)存在一定誤差,這是由于本章所采用的 σ的值為統(tǒng)計(jì)得到的均值所致,但即使是差別最大的城市,在50年重現(xiàn)期下,二者的絕對(duì)誤差在0.2%~10%,且該誤差隨重現(xiàn)期的增加而減小,使修正后的6 h最后2 min 的預(yù)測(cè)結(jié)果更加接近連續(xù)2 min 的預(yù)測(cè)結(jié)果。因此,這種通過原始數(shù)據(jù)入手的修正方法是可靠且可行的。

    表4 每6 h 最后2 min 觀測(cè)間隔在10年重現(xiàn)期下修正前后風(fēng)速極值對(duì)比Table 4 Comparison of extremewind speeds before and after correction under 10-year return period of 2-minute per hour

    表5 每6 h 最后2 min 觀測(cè)間隔在50年重現(xiàn)期下修正前后風(fēng)速極值對(duì)比Table 5 Comparison of extremewind speeds before and after correction under 50-year return period of 2-minute per hour

    5 結(jié)論

    本文討論了風(fēng)速樣本選取間隔對(duì)極值風(fēng)速預(yù)測(cè)的影響。對(duì)美國11座城市氣象站記錄的風(fēng)速每分鐘連續(xù)觀測(cè)數(shù)據(jù)進(jìn)行處理,處理為連續(xù)2 min 平均值序列、觀測(cè)間隔為1 h 和6 h 的2 min 平均值序列3種觀測(cè)間隔的方式。結(jié)合改進(jìn)的獨(dú)立風(fēng)暴法和考慮風(fēng)向相關(guān)性的風(fēng)速概率模型進(jìn)行不同重現(xiàn)期下極值風(fēng)速預(yù)測(cè),分析其結(jié)果差異,并給出了一種從原始數(shù)據(jù)入手的修正方法,簡便且準(zhǔn)確率。得到如下結(jié)論:

    (1)利用不同觀測(cè)間隔得到的樣本進(jìn)行不同重現(xiàn)期極值風(fēng)速預(yù)測(cè)結(jié)果具有10%~20%的誤差率。

    (2)風(fēng)速的連續(xù)觀測(cè)值與間隔1 h 的觀測(cè)值的日極值風(fēng)速差值加上0.5 m/s后滿足參數(shù)為1.56的瑞利分布,即:

    (3)風(fēng)速的連續(xù)觀測(cè)值與間隔6 h 觀測(cè)值的日極值風(fēng)速的差值加上0.2 m/s后滿足參數(shù)為2.48的瑞利分布,即:

    (4)利用基于瑞利分布的修正方法對(duì)有間隔的觀測(cè)數(shù)據(jù)進(jìn)行修正,可以明顯降低基于觀測(cè)數(shù)據(jù)的極值估算誤差,對(duì)于間隔1 h 的數(shù)據(jù),修正后給出的極值估算誤差率均值降低6%;觀測(cè)間隔為6 h的數(shù)據(jù),這個(gè)誤差可降低到3%。

    猜你喜歡
    原始數(shù)據(jù)風(fēng)向極值
    GOLDEN OPPORTUNITY FOR CHINA-INDONESIA COOPERATION
    極值點(diǎn)帶你去“漂移”
    受特定變化趨勢(shì)限制的傳感器數(shù)據(jù)處理方法研究
    極值點(diǎn)偏移攔路,三法可取
    一類“極值點(diǎn)偏移”問題的解法與反思
    全新Mentor DRS360 平臺(tái)借助集中式原始數(shù)據(jù)融合及直接實(shí)時(shí)傳感技術(shù)實(shí)現(xiàn)5 級(jí)自動(dòng)駕駛
    汽車零部件(2017年4期)2017-07-12 17:05:53
    自然與風(fēng)Feeling Nature
    行業(yè)統(tǒng)計(jì)帶來哪些風(fēng)向?
    匹配數(shù)為1的極值2-均衡4-部4-圖的結(jié)構(gòu)
    風(fēng)向
    風(fēng)能(2015年8期)2015-02-27 10:15:11
    国产又黄又爽又无遮挡在线| 国产私拍福利视频在线观看| 欧美一区二区国产精品久久精品| 一夜夜www| 午夜激情福利司机影院| 人人妻人人看人人澡| 999久久久精品免费观看国产| 亚洲片人在线观看| 久久精品综合一区二区三区| bbb黄色大片| 国产 一区 欧美 日韩| 亚洲av五月六月丁香网| 国产主播在线观看一区二区| 欧美成人性av电影在线观看| 国产亚洲欧美98| 国产精品电影一区二区三区| 国产午夜福利久久久久久| 午夜免费观看网址| 国产成人影院久久av| 99久国产av精品| 中文资源天堂在线| 色在线成人网| 欧美zozozo另类| av片东京热男人的天堂| 中文字幕高清在线视频| 亚洲中文字幕一区二区三区有码在线看 | 精品国产美女av久久久久小说| 亚洲性夜色夜夜综合| 在线观看免费午夜福利视频| 99精品在免费线老司机午夜| 午夜福利高清视频| 国产av不卡久久| 国产亚洲av嫩草精品影院| 亚洲第一电影网av| 精品国产超薄肉色丝袜足j| 国产不卡一卡二| 狠狠狠狠99中文字幕| 欧美中文日本在线观看视频| 亚洲国产日韩欧美精品在线观看 | 长腿黑丝高跟| bbb黄色大片| 午夜久久久久精精品| 亚洲av五月六月丁香网| 久久精品国产亚洲av香蕉五月| 国产精品av久久久久免费| 欧美色视频一区免费| 日本一本二区三区精品| 日韩欧美 国产精品| 国产精品一区二区三区四区久久| 久久九九热精品免费| 亚洲国产精品久久男人天堂| 熟女电影av网| 亚洲av免费在线观看| 高清毛片免费观看视频网站| 精品国内亚洲2022精品成人| 亚洲av成人不卡在线观看播放网| 琪琪午夜伦伦电影理论片6080| 香蕉国产在线看| 亚洲第一电影网av| 国产视频一区二区在线看| 97人妻精品一区二区三区麻豆| 亚洲无线在线观看| 女人被狂操c到高潮| 国产成+人综合+亚洲专区| 亚洲自拍偷在线| 日本免费a在线| 欧美黄色淫秽网站| 久久精品亚洲精品国产色婷小说| 91av网一区二区| а√天堂www在线а√下载| 久久久久国产一级毛片高清牌| 麻豆一二三区av精品| 久久久久久人人人人人| 天天躁狠狠躁夜夜躁狠狠躁| 国产亚洲欧美在线一区二区| 日韩欧美免费精品| 欧美黄色淫秽网站| 亚洲自偷自拍图片 自拍| 国产精品影院久久| 欧美激情在线99| 蜜桃久久精品国产亚洲av| 曰老女人黄片| 久久精品影院6| 两性夫妻黄色片| 免费看美女性在线毛片视频| 香蕉丝袜av| 99国产综合亚洲精品| 精华霜和精华液先用哪个| 宅男免费午夜| 香蕉丝袜av| 午夜亚洲福利在线播放| 成人无遮挡网站| 日韩有码中文字幕| 午夜成年电影在线免费观看| 亚洲国产精品久久男人天堂| 小蜜桃在线观看免费完整版高清| 亚洲av片天天在线观看| 757午夜福利合集在线观看| 好男人电影高清在线观看| 哪里可以看免费的av片| 国产69精品久久久久777片 | 老鸭窝网址在线观看| 91麻豆av在线| 婷婷亚洲欧美| 长腿黑丝高跟| 久久这里只有精品19| 日本一本二区三区精品| 人妻久久中文字幕网| 国产日本99.免费观看| 久久久国产成人免费| 亚洲第一欧美日韩一区二区三区| 亚洲人成电影免费在线| 免费在线观看亚洲国产| 亚洲国产精品久久男人天堂| 成人精品一区二区免费| 网址你懂的国产日韩在线| 嫩草影视91久久| 村上凉子中文字幕在线| 在线免费观看的www视频| 亚洲国产欧美人成| 桃红色精品国产亚洲av| 久久久久久久久中文| 精品久久久久久,| 日韩av在线大香蕉| 中文字幕久久专区| 97人妻精品一区二区三区麻豆| 免费大片18禁| www.熟女人妻精品国产| 国语自产精品视频在线第100页| av福利片在线观看| 午夜福利在线在线| 一进一出好大好爽视频| 51午夜福利影视在线观看| 美女cb高潮喷水在线观看 | 亚洲欧美日韩东京热| 黄色 视频免费看| 亚洲精品一区av在线观看| 免费电影在线观看免费观看| 最近视频中文字幕2019在线8| 国产高清有码在线观看视频| 久久久国产成人免费| 欧美另类亚洲清纯唯美| 成人一区二区视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 日韩免费av在线播放| 99在线人妻在线中文字幕| 热99re8久久精品国产| 国产伦在线观看视频一区| 最新中文字幕久久久久 | 搡老熟女国产l中国老女人| 国产又色又爽无遮挡免费看| 午夜免费激情av| 免费一级毛片在线播放高清视频| 男女午夜视频在线观看| 天天躁日日操中文字幕| 中文字幕久久专区| 亚洲片人在线观看| 精品无人区乱码1区二区| 欧美日韩瑟瑟在线播放| 亚洲精品粉嫩美女一区| 色av中文字幕| 男人舔女人的私密视频| 99精品欧美一区二区三区四区| 午夜免费观看网址| av福利片在线观看| 日韩 欧美 亚洲 中文字幕| 岛国视频午夜一区免费看| 亚洲五月天丁香| 国产激情久久老熟女| 免费在线观看日本一区| 国产99白浆流出| 99精品欧美一区二区三区四区| 长腿黑丝高跟| 午夜精品在线福利| 亚洲在线自拍视频| 国产精品国产高清国产av| 精品乱码久久久久久99久播| av视频在线观看入口| 国产日本99.免费观看| 亚洲国产欧美网| 又黄又粗又硬又大视频| 此物有八面人人有两片| 久久99热这里只有精品18| 18禁美女被吸乳视频| 精品一区二区三区视频在线 | 欧美色视频一区免费| 国产精品一区二区免费欧美| 在线观看66精品国产| 久久国产乱子伦精品免费另类| 色综合婷婷激情| 国产熟女xx| 久久欧美精品欧美久久欧美| АⅤ资源中文在线天堂| 99精品久久久久人妻精品| 亚洲国产色片| 成人欧美大片| av在线蜜桃| 男人舔女人下体高潮全视频| 中国美女看黄片| 老鸭窝网址在线观看| 婷婷丁香在线五月| 欧美一级a爱片免费观看看| 日本成人三级电影网站| 日韩三级视频一区二区三区| 免费大片18禁| 老司机午夜十八禁免费视频| 国产 一区 欧美 日韩| 欧美三级亚洲精品| 黄片小视频在线播放| 狠狠狠狠99中文字幕| 国产成人av教育| 国产又黄又爽又无遮挡在线| 国产精品久久久久久亚洲av鲁大| 亚洲国产看品久久| 听说在线观看完整版免费高清| 亚洲成人精品中文字幕电影| 一本久久中文字幕| 国产激情久久老熟女| 亚洲国产欧美人成| 高清毛片免费观看视频网站| 18禁黄网站禁片午夜丰满| 亚洲成a人片在线一区二区| 啦啦啦韩国在线观看视频| 精品电影一区二区在线| 男女之事视频高清在线观看| 午夜福利成人在线免费观看| 亚洲欧美日韩东京热| av欧美777| 18禁黄网站禁片免费观看直播| 欧美性猛交黑人性爽| 日韩有码中文字幕| 精品久久久久久久人妻蜜臀av| 夜夜看夜夜爽夜夜摸| www.999成人在线观看| 啪啪无遮挡十八禁网站| 国产单亲对白刺激| 岛国在线观看网站| 久久久国产欧美日韩av| а√天堂www在线а√下载| 免费看日本二区| 亚洲一区二区三区不卡视频| 日本与韩国留学比较| 嫩草影院入口| 亚洲最大成人中文| 亚洲九九香蕉| 久久久精品大字幕| 午夜福利欧美成人| 非洲黑人性xxxx精品又粗又长| h日本视频在线播放| 天天添夜夜摸| 精品久久久久久久久久免费视频| 在线观看一区二区三区| 国产淫片久久久久久久久 | 美女黄网站色视频| 久久久国产成人免费| 欧美乱色亚洲激情| 国产91精品成人一区二区三区| 男插女下体视频免费在线播放| 亚洲欧洲精品一区二区精品久久久| 欧美乱色亚洲激情| 男女午夜视频在线观看| 精品久久久久久久毛片微露脸| 精品一区二区三区四区五区乱码| 真人做人爱边吃奶动态| 日韩免费av在线播放| 亚洲精华国产精华精| 国产熟女xx| 国产成人福利小说| 亚洲在线观看片| 久久精品aⅴ一区二区三区四区| 一级作爱视频免费观看| 亚洲人成网站在线播放欧美日韩| 一个人免费在线观看电影 | 九色成人免费人妻av| 在线播放国产精品三级| 亚洲片人在线观看| 熟女少妇亚洲综合色aaa.| 成人av一区二区三区在线看| 桃红色精品国产亚洲av| 麻豆成人av在线观看| а√天堂www在线а√下载| 国产成人精品久久二区二区免费| 国产精品一区二区三区四区久久| 成人18禁在线播放| 一a级毛片在线观看| 国产视频内射| 国产精品久久久久久亚洲av鲁大| 最近在线观看免费完整版| 久久国产精品影院| 18禁美女被吸乳视频| 日韩欧美在线乱码| 九九热线精品视视频播放| 国产成人av激情在线播放| 97人妻精品一区二区三区麻豆| 午夜精品久久久久久毛片777| 久久久久国产精品人妻aⅴ院| 麻豆国产97在线/欧美| 伦理电影免费视频| 长腿黑丝高跟| 中文字幕最新亚洲高清| 婷婷六月久久综合丁香| 日本三级黄在线观看| 最好的美女福利视频网| 桃色一区二区三区在线观看| 1000部很黄的大片| 99在线视频只有这里精品首页| 又黄又爽又免费观看的视频| 欧美成人性av电影在线观看| 少妇熟女aⅴ在线视频| 九九热线精品视视频播放| 怎么达到女性高潮| 天堂√8在线中文| 中文在线观看免费www的网站| 搡老岳熟女国产| 久9热在线精品视频| 搞女人的毛片| 精品国产乱码久久久久久男人| 一级毛片女人18水好多| АⅤ资源中文在线天堂| 国产探花在线观看一区二区| 国产欧美日韩精品亚洲av| 色尼玛亚洲综合影院| 中文字幕人成人乱码亚洲影| www国产在线视频色| 免费看a级黄色片| 国产成人aa在线观看| 身体一侧抽搐| 欧美色视频一区免费| 亚洲国产高清在线一区二区三| 欧美激情在线99| 免费在线观看视频国产中文字幕亚洲| a级毛片a级免费在线| 日日夜夜操网爽| 久久精品国产清高在天天线| 久久久久精品国产欧美久久久| 一级a爱片免费观看的视频| 国内精品久久久久久久电影| 国产精品影院久久| 成年免费大片在线观看| 国产精品久久久av美女十八| e午夜精品久久久久久久| 人妻夜夜爽99麻豆av| 桃红色精品国产亚洲av| 日本黄大片高清| 桃红色精品国产亚洲av| 亚洲一区二区三区色噜噜| 国产亚洲av嫩草精品影院| 亚洲熟女毛片儿| cao死你这个sao货| 搞女人的毛片| 欧美zozozo另类| 成人高潮视频无遮挡免费网站| 成在线人永久免费视频| 国产精品野战在线观看| 国产精品日韩av在线免费观看| 最新美女视频免费是黄的| 一进一出好大好爽视频| 午夜久久久久精精品| 免费av毛片视频| 亚洲第一电影网av| 国产精品久久久久久精品电影| 欧美极品一区二区三区四区| 在线观看一区二区三区| 国产成人啪精品午夜网站| 搡老岳熟女国产| 国产亚洲欧美98| 亚洲精品一区av在线观看| 99久久精品国产亚洲精品| 国产高清视频在线播放一区| tocl精华| 麻豆久久精品国产亚洲av| 毛片女人毛片| 国产一级毛片七仙女欲春2| 青草久久国产| 一区二区三区激情视频| 丝袜人妻中文字幕| 久久香蕉精品热| 这个男人来自地球电影免费观看| 国产亚洲av高清不卡| 嫩草影院入口| 亚洲 欧美 日韩 在线 免费| av女优亚洲男人天堂 | 亚洲第一欧美日韩一区二区三区| 成人精品一区二区免费| 丰满的人妻完整版| 黄色片一级片一级黄色片| 国产欧美日韩精品一区二区| 欧美又色又爽又黄视频| 国产亚洲av嫩草精品影院| 亚洲中文日韩欧美视频| 国产亚洲精品久久久久久毛片| www.999成人在线观看| 午夜免费激情av| 亚洲无线观看免费| 精品电影一区二区在线| 九九热线精品视视频播放| 色视频www国产| 婷婷亚洲欧美| 亚洲18禁久久av| 桃红色精品国产亚洲av| 午夜精品一区二区三区免费看| 日韩欧美 国产精品| 黄片大片在线免费观看| 久久草成人影院| 成人精品一区二区免费| 丁香欧美五月| x7x7x7水蜜桃| 麻豆一二三区av精品| 97超视频在线观看视频| 免费看日本二区| 不卡一级毛片| 性色avwww在线观看| 特级一级黄色大片| 久久这里只有精品中国| 色视频www国产| 欧美日韩黄片免| 久久中文字幕一级| 校园春色视频在线观看| 国内少妇人妻偷人精品xxx网站 | 亚洲欧洲精品一区二区精品久久久| 男人舔奶头视频| 亚洲国产欧美一区二区综合| 亚洲国产欧美网| 欧美黑人欧美精品刺激| 成人亚洲精品av一区二区| 999精品在线视频| 亚洲七黄色美女视频| 99精品在免费线老司机午夜| 久久久成人免费电影| 又爽又黄无遮挡网站| 亚洲av成人不卡在线观看播放网| 久久亚洲真实| 真人做人爱边吃奶动态| 91av网一区二区| 99久久无色码亚洲精品果冻| 久久久精品大字幕| 巨乳人妻的诱惑在线观看| 日本一二三区视频观看| 亚洲va日本ⅴa欧美va伊人久久| 99久久精品热视频| 欧美另类亚洲清纯唯美| 久久欧美精品欧美久久欧美| 18禁美女被吸乳视频| 中文字幕高清在线视频| 小蜜桃在线观看免费完整版高清| 免费看a级黄色片| 夜夜爽天天搞| 国内精品一区二区在线观看| 一个人免费在线观看电影 | 夜夜看夜夜爽夜夜摸| 搡老妇女老女人老熟妇| 高潮久久久久久久久久久不卡| 久久国产精品影院| 日韩欧美三级三区| 在线免费观看不下载黄p国产 | 波多野结衣高清作品| 亚洲精品乱码久久久v下载方式 | 不卡一级毛片| 久久久久久大精品| 少妇的丰满在线观看| 国产亚洲精品久久久久久毛片| 一个人免费在线观看电影 | 日本精品一区二区三区蜜桃| 精品日产1卡2卡| 日韩精品青青久久久久久| 日韩有码中文字幕| 大型黄色视频在线免费观看| 亚洲av日韩精品久久久久久密| 岛国在线免费视频观看| 久久久久久久久中文| 波多野结衣高清作品| 天堂网av新在线| 国产精品久久久久久精品电影| 国产成人欧美在线观看| 国产成人系列免费观看| 色视频www国产| 欧美+亚洲+日韩+国产| 一本久久中文字幕| 久久香蕉精品热| 最近最新免费中文字幕在线| 在线观看舔阴道视频| 久久精品影院6| 黑人欧美特级aaaaaa片| 久久这里只有精品中国| 亚洲avbb在线观看| 美女午夜性视频免费| aaaaa片日本免费| 999久久久国产精品视频| 毛片女人毛片| 日本黄大片高清| 亚洲成a人片在线一区二区| 久久精品aⅴ一区二区三区四区| 观看免费一级毛片| 黄色视频,在线免费观看| 婷婷精品国产亚洲av| 精品一区二区三区视频在线 | 97超视频在线观看视频| 欧美日韩国产亚洲二区| 嫁个100分男人电影在线观看| 我的老师免费观看完整版| 最新在线观看一区二区三区| 婷婷六月久久综合丁香| 麻豆一二三区av精品| 欧美zozozo另类| 我的老师免费观看完整版| 国产单亲对白刺激| 岛国视频午夜一区免费看| 国产熟女xx| www.www免费av| 国产精品久久久久久亚洲av鲁大| 校园春色视频在线观看| 国产一区二区在线av高清观看| 黑人操中国人逼视频| 一个人观看的视频www高清免费观看 | 欧美又色又爽又黄视频| 欧美国产日韩亚洲一区| 成人特级av手机在线观看| 亚洲色图av天堂| 成年人黄色毛片网站| 精品电影一区二区在线| 很黄的视频免费| 精品国产美女av久久久久小说| 可以在线观看毛片的网站| 欧美成人性av电影在线观看| 亚洲成人久久性| 桃色一区二区三区在线观看| 亚洲,欧美精品.| 国产av麻豆久久久久久久| 日韩av在线大香蕉| 男人舔女人的私密视频| 成人亚洲精品av一区二区| 99在线人妻在线中文字幕| 欧美激情久久久久久爽电影| 久久性视频一级片| av中文乱码字幕在线| 亚洲天堂国产精品一区在线| a级毛片a级免费在线| 在线a可以看的网站| 老司机午夜十八禁免费视频| 亚洲片人在线观看| 两个人看的免费小视频| 亚洲中文av在线| 成人国产综合亚洲| 国产麻豆成人av免费视频| 亚洲美女视频黄频| 99国产精品99久久久久| 免费人成视频x8x8入口观看| 国产爱豆传媒在线观看| 久久伊人香网站| 色综合站精品国产| 又紧又爽又黄一区二区| 中文字幕精品亚洲无线码一区| 午夜a级毛片| 一级毛片高清免费大全| 久久精品91无色码中文字幕| 精品国产超薄肉色丝袜足j| e午夜精品久久久久久久| www日本黄色视频网| 国产精品野战在线观看| 欧美xxxx黑人xx丫x性爽| 免费看a级黄色片| 亚洲av成人精品一区久久| 亚洲欧美一区二区三区黑人| 动漫黄色视频在线观看| 成人三级做爰电影| 亚洲国产欧洲综合997久久,| 亚洲午夜精品一区,二区,三区| 久久久国产成人免费| 久久亚洲精品不卡| 一级毛片高清免费大全| 又大又爽又粗| 日韩欧美国产一区二区入口| 色尼玛亚洲综合影院| 极品教师在线免费播放| 在线观看免费视频日本深夜| 首页视频小说图片口味搜索| 国产av不卡久久| 黄色日韩在线| 两人在一起打扑克的视频| 一区福利在线观看| 真人做人爱边吃奶动态| a级毛片在线看网站| 国产三级在线视频| 亚洲在线自拍视频| 国产精品久久久久久人妻精品电影| 一进一出抽搐gif免费好疼| 色尼玛亚洲综合影院| 日韩欧美免费精品| 国产精品综合久久久久久久免费| 亚洲av电影不卡..在线观看| 国产成人福利小说| 黄色日韩在线| 成人鲁丝片一二三区免费| 变态另类丝袜制服| 麻豆成人av在线观看| www国产在线视频色| 亚洲中文av在线| 欧美日韩瑟瑟在线播放| 精品熟女少妇八av免费久了| 国产不卡一卡二| 久久精品国产综合久久久| 日韩欧美免费精品| 两个人的视频大全免费| 高清在线国产一区| 性色av乱码一区二区三区2| 动漫黄色视频在线观看| 99热这里只有精品一区 | 亚洲国产中文字幕在线视频| 欧美一区二区国产精品久久精品| 亚洲av五月六月丁香网| 国产视频内射| 国产成人福利小说| 熟女电影av网| 美女扒开内裤让男人捅视频| 床上黄色一级片|