• <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
    色婷婷久久久亚洲欧美| 亚洲情色 制服丝袜| 丰满少妇做爰视频| 最后的刺客免费高清国语| 啦啦啦中文免费视频观看日本| 久久 成人 亚洲| 久久人人爽人人爽人人片va| 大片免费播放器 马上看| 欧美日韩国产mv在线观看视频| 如何舔出高潮| 亚洲国产日韩一区二区| 亚洲色图 男人天堂 中文字幕 | 最近的中文字幕免费完整| 国产精品麻豆人妻色哟哟久久| 黄色配什么色好看| 少妇被粗大猛烈的视频| 国产成人免费无遮挡视频| 黄色毛片三级朝国网站| 婷婷色av中文字幕| 99国产精品免费福利视频| 交换朋友夫妻互换小说| 午夜久久久在线观看| 亚洲美女黄色视频免费看| 男女国产视频网站| 韩国高清视频一区二区三区| 色5月婷婷丁香| 亚洲精品自拍成人| 最新的欧美精品一区二区| 欧美精品国产亚洲| 男女高潮啪啪啪动态图| 91精品国产国语对白视频| 18禁在线无遮挡免费观看视频| 男男h啪啪无遮挡| 精品国产国语对白av| 我要看黄色一级片免费的| 黄色视频在线播放观看不卡| 欧美丝袜亚洲另类| 久久精品国产a三级三级三级| 国产xxxxx性猛交| 欧美变态另类bdsm刘玥| 国产精品不卡视频一区二区| 国产高清不卡午夜福利| 99热网站在线观看| 国产精品熟女久久久久浪| 老女人水多毛片| 日韩熟女老妇一区二区性免费视频| 亚洲精品aⅴ在线观看| 少妇高潮的动态图| 捣出白浆h1v1| 欧美xxⅹ黑人| 午夜福利网站1000一区二区三区| 美女大奶头黄色视频| 午夜av观看不卡| 一区二区av电影网| 亚洲综合色网址| 搡老乐熟女国产| av卡一久久| 热99国产精品久久久久久7| 久久人人爽人人片av| 久久久久久人人人人人| 免费少妇av软件| 国产免费现黄频在线看| 日韩不卡一区二区三区视频在线| 人妻人人澡人人爽人人| 国产色婷婷99| 日韩av不卡免费在线播放| 国国产精品蜜臀av免费| 另类精品久久| 亚洲久久久国产精品| 国产精品成人在线| 国产1区2区3区精品| 纯流量卡能插随身wifi吗| 九草在线视频观看| 日韩大片免费观看网站| 久久精品夜色国产| 亚洲在久久综合| 亚洲色图 男人天堂 中文字幕 | 国产色爽女视频免费观看| 国产av精品麻豆| 久久久久久伊人网av| 一区二区三区乱码不卡18| 涩涩av久久男人的天堂| 男女啪啪激烈高潮av片| 国精品久久久久久国模美| 久久久久精品久久久久真实原创| 中文字幕精品免费在线观看视频 | 国产精品欧美亚洲77777| 人妻少妇偷人精品九色| 80岁老熟妇乱子伦牲交| 母亲3免费完整高清在线观看 | 有码 亚洲区| 国产精品99久久99久久久不卡 | 国产高清三级在线| 男女下面插进去视频免费观看 | 欧美+日韩+精品| 成人手机av| 久久热在线av| av线在线观看网站| 亚洲欧美精品自产自拍| 亚洲av电影在线观看一区二区三区| 欧美3d第一页| 国产免费现黄频在线看| 99久久中文字幕三级久久日本| 国内精品宾馆在线| 成年动漫av网址| 国产黄色免费在线视频| 国精品久久久久久国模美| 99久久人妻综合| 精品国产露脸久久av麻豆| 日韩大片免费观看网站| 久久久久精品性色| 亚洲国产日韩一区二区| 各种免费的搞黄视频| 久久久久网色| 欧美激情 高清一区二区三区| 亚洲美女黄色视频免费看| 一二三四中文在线观看免费高清| 又黄又爽又刺激的免费视频.| 精品福利永久在线观看| 日日摸夜夜添夜夜爱| 黑人巨大精品欧美一区二区蜜桃 | 麻豆乱淫一区二区| 亚洲成色77777| 成人国产av品久久久| 国产极品粉嫩免费观看在线| 国产亚洲精品久久久com| 男人添女人高潮全过程视频| 日本av免费视频播放| 久久 成人 亚洲| 国产亚洲av片在线观看秒播厂| 永久免费av网站大全| 久久久精品94久久精品| 久久久久久久国产电影| 免费观看在线日韩| 国国产精品蜜臀av免费| 丝瓜视频免费看黄片| 一级毛片黄色毛片免费观看视频| 男女高潮啪啪啪动态图| 亚洲精品一区蜜桃| 亚洲av中文av极速乱| 观看美女的网站| 久久久久久人妻| 中文字幕人妻丝袜制服| 男女国产视频网站| 十八禁高潮呻吟视频| 最近最新中文字幕大全免费视频 | 亚洲成av片中文字幕在线观看 | 国产欧美另类精品又又久久亚洲欧美| 97精品久久久久久久久久精品| 高清视频免费观看一区二区| 最近最新中文字幕免费大全7| 9热在线视频观看99| 亚洲精品一二三| 日韩中字成人| 免费日韩欧美在线观看| 不卡视频在线观看欧美| 爱豆传媒免费全集在线观看| 青春草国产在线视频| 成人免费观看视频高清| 高清av免费在线| 三级国产精品片| 国产一区二区激情短视频 | 亚洲高清免费不卡视频| 另类亚洲欧美激情| 国产一区二区在线观看av| 亚洲av免费高清在线观看| 另类亚洲欧美激情| 日韩不卡一区二区三区视频在线| 欧美xxxx性猛交bbbb| 精品少妇黑人巨大在线播放| 亚洲成国产人片在线观看| 熟女人妻精品中文字幕| 久久女婷五月综合色啪小说| 国产片内射在线| 午夜精品国产一区二区电影| 久久久久久久久久久久大奶| 久久国产亚洲av麻豆专区| 精品亚洲成国产av| 两性夫妻黄色片 | 久久鲁丝午夜福利片| 日本vs欧美在线观看视频| 亚洲国产色片| 久久久国产精品麻豆| 自拍欧美九色日韩亚洲蝌蚪91| 久久99热6这里只有精品| 性高湖久久久久久久久免费观看| 激情五月婷婷亚洲| 在线观看人妻少妇| 两性夫妻黄色片 | 久久鲁丝午夜福利片| 爱豆传媒免费全集在线观看| 中文字幕制服av| 精品熟女少妇av免费看| 青青草视频在线视频观看| 亚洲少妇的诱惑av| 中文欧美无线码| 欧美日韩国产mv在线观看视频| 午夜视频国产福利| 七月丁香在线播放| 久久久久久人人人人人| 2022亚洲国产成人精品| 欧美国产精品va在线观看不卡| 亚洲精品一二三| 日本vs欧美在线观看视频| 狂野欧美激情性bbbbbb| 最近中文字幕2019免费版| 欧美少妇被猛烈插入视频| 国产爽快片一区二区三区| 永久免费av网站大全| 亚洲欧美日韩卡通动漫| 夜夜爽夜夜爽视频| 国产精品一国产av| 国产精品一区www在线观看| 久久久久国产网址| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av综合色区一区| 香蕉丝袜av| 久久热在线av| 精品一品国产午夜福利视频| 久久这里只有精品19| 欧美精品亚洲一区二区| 成人免费观看视频高清| 在线观看美女被高潮喷水网站| 国产日韩欧美亚洲二区| 亚洲天堂av无毛| www日本在线高清视频| 色5月婷婷丁香| 少妇的丰满在线观看| 99久久精品国产国产毛片| 97超碰精品成人国产| 少妇人妻久久综合中文| 999精品在线视频| 亚洲欧美一区二区三区黑人 | 天堂俺去俺来也www色官网| 你懂的网址亚洲精品在线观看| 欧美性感艳星| 国产精品麻豆人妻色哟哟久久| 欧美日韩视频精品一区| 9热在线视频观看99| 久久久久久久久久久免费av| 国产片内射在线| 国产日韩欧美亚洲二区| 最新中文字幕久久久久| 少妇熟女欧美另类| 亚洲一码二码三码区别大吗| 9191精品国产免费久久| 免费在线观看完整版高清| 精品久久蜜臀av无| 纯流量卡能插随身wifi吗| 一区二区三区精品91| 欧美+日韩+精品| 亚洲情色 制服丝袜| 久久国产精品大桥未久av| 性高湖久久久久久久久免费观看| 三上悠亚av全集在线观看| 18禁动态无遮挡网站| 免费看av在线观看网站| 久久久久久久大尺度免费视频| 九色亚洲精品在线播放| 国产精品一区www在线观看| 在线观看www视频免费| a级毛色黄片| 十八禁高潮呻吟视频| 国产色婷婷99| 国产精品不卡视频一区二区| 老司机影院毛片| av网站免费在线观看视频| 只有这里有精品99| 精品一区二区免费观看| 大片电影免费在线观看免费| xxxhd国产人妻xxx| 夜夜骑夜夜射夜夜干| 美女国产高潮福利片在线看| 久久久精品区二区三区| 91成人精品电影| 亚洲人成网站在线观看播放| 在线天堂中文资源库| 边亲边吃奶的免费视频| 超色免费av| 人人妻人人澡人人看| 亚洲av福利一区| 综合色丁香网| 九色亚洲精品在线播放| 国产精品久久久久久精品古装| 久久久国产欧美日韩av| 97超碰精品成人国产| 伦精品一区二区三区| 伊人亚洲综合成人网| 日本wwww免费看| 午夜久久久在线观看| 久久 成人 亚洲| 日韩 亚洲 欧美在线| 老熟女久久久| 宅男免费午夜| 国产av精品麻豆| 汤姆久久久久久久影院中文字幕| 国产国语露脸激情在线看| 伦理电影免费视频| 最近手机中文字幕大全| 亚洲国产成人一精品久久久| 美女中出高潮动态图| 国产精品.久久久| 美女大奶头黄色视频| 免费播放大片免费观看视频在线观看| 美女xxoo啪啪120秒动态图| 在现免费观看毛片| 久久精品国产亚洲av天美| 赤兔流量卡办理| 成人毛片60女人毛片免费| 精品国产乱码久久久久久小说| 亚洲精品国产色婷婷电影| 国内精品宾馆在线| 女人被躁到高潮嗷嗷叫费观| 国产淫语在线视频| 波野结衣二区三区在线| 国产黄频视频在线观看| 精品午夜福利在线看| 高清不卡的av网站| 亚洲人与动物交配视频| 精品亚洲成国产av| 99国产精品免费福利视频| 美女大奶头黄色视频| 国产男人的电影天堂91| 色视频在线一区二区三区| 少妇 在线观看| 日日啪夜夜爽| 亚洲久久久国产精品| 亚洲精品aⅴ在线观看| 国产一区二区在线观看av| 青青草视频在线视频观看| 国产成人aa在线观看| 亚洲三级黄色毛片| 交换朋友夫妻互换小说| 国产精品三级大全| 欧美成人午夜精品| 日本欧美国产在线视频| 精品99又大又爽又粗少妇毛片| 日韩,欧美,国产一区二区三区| 美女国产视频在线观看| 精品第一国产精品| 熟女人妻精品中文字幕| 国产xxxxx性猛交| 欧美日韩一区二区视频在线观看视频在线| 国产麻豆69| 午夜91福利影院| 国产日韩欧美在线精品| 国产永久视频网站| 青春草亚洲视频在线观看| 成人二区视频| 亚洲国产日韩一区二区| 波多野结衣一区麻豆| 日本爱情动作片www.在线观看| 亚洲精品国产av蜜桃| 久久久国产精品麻豆| av片东京热男人的天堂| 欧美人与性动交α欧美软件 | 蜜臀久久99精品久久宅男| 国产精品一区二区在线不卡| 蜜桃国产av成人99| videosex国产| 少妇的丰满在线观看| 精品一区二区三区视频在线| 日本91视频免费播放| 男女下面插进去视频免费观看 | 男女午夜视频在线观看 | 亚洲婷婷狠狠爱综合网| 亚洲国产毛片av蜜桃av| 国产永久视频网站| 欧美性感艳星| 纯流量卡能插随身wifi吗| 亚洲成色77777| 国产国拍精品亚洲av在线观看| 国产精品99久久99久久久不卡 | 免费大片18禁| 在线观看www视频免费| 亚洲,欧美,日韩| 国产国语露脸激情在线看| 美国免费a级毛片| 国产xxxxx性猛交| 一区二区三区精品91| 亚洲,一卡二卡三卡| 内地一区二区视频在线| 在线天堂最新版资源| 大片电影免费在线观看免费| 亚洲精品第二区| 亚洲成av片中文字幕在线观看 | 亚洲精品美女久久久久99蜜臀 | 最黄视频免费看| 91精品国产国语对白视频| 极品少妇高潮喷水抽搐| 全区人妻精品视频| 色哟哟·www| 日韩一区二区三区影片| 亚洲中文av在线| 亚洲一级一片aⅴ在线观看| 在线 av 中文字幕| 亚洲精华国产精华液的使用体验| 亚洲熟女精品中文字幕| 亚洲成人手机| 国产日韩欧美在线精品| 超色免费av| av在线播放精品| 国产色婷婷99| 欧美日韩视频精品一区| 成人无遮挡网站| 日韩电影二区| 成人黄色视频免费在线看| 一二三四中文在线观看免费高清| 亚洲国产精品国产精品| 国产女主播在线喷水免费视频网站| 97在线视频观看| 中文字幕制服av| 日本免费在线观看一区| 有码 亚洲区| 香蕉精品网在线| 狠狠婷婷综合久久久久久88av| 丰满乱子伦码专区| 国产精品蜜桃在线观看| 插逼视频在线观看| 丰满乱子伦码专区| 韩国精品一区二区三区 | 午夜福利影视在线免费观看| 两个人免费观看高清视频| 波多野结衣一区麻豆| 色婷婷久久久亚洲欧美| 人妻一区二区av| 亚洲一级一片aⅴ在线观看| 妹子高潮喷水视频| 免费观看av网站的网址| 久久人人爽av亚洲精品天堂| 91成人精品电影| 大香蕉久久成人网| 一本大道久久a久久精品| 高清视频免费观看一区二区| 桃花免费在线播放| 秋霞在线观看毛片| 国产av精品麻豆| 在线观看美女被高潮喷水网站| 欧美精品高潮呻吟av久久| 欧美人与性动交α欧美精品济南到 | 一二三四在线观看免费中文在 | 在线观看免费日韩欧美大片| 亚洲高清免费不卡视频| 午夜福利乱码中文字幕| 国产一区二区在线观看日韩| a级毛片黄视频| 看非洲黑人一级黄片| 99九九在线精品视频| a级片在线免费高清观看视频| 美女国产视频在线观看| 久久精品久久久久久久性| 最近中文字幕2019免费版| 国产亚洲午夜精品一区二区久久| 男人舔女人的私密视频| 亚洲av国产av综合av卡| 国产在视频线精品| 日韩在线高清观看一区二区三区| 色婷婷av一区二区三区视频| 久久久精品区二区三区| 边亲边吃奶的免费视频| 18禁在线无遮挡免费观看视频| 国产精品国产av在线观看| 色婷婷av一区二区三区视频| 国产69精品久久久久777片| 在线看a的网站| 精品少妇内射三级| 麻豆精品久久久久久蜜桃| 色网站视频免费| 日本爱情动作片www.在线观看| 国产午夜精品一二区理论片| 亚洲一级一片aⅴ在线观看| 国产精品久久久久久久久免| 亚洲,欧美精品.| 国产免费一级a男人的天堂| av卡一久久| 精品少妇黑人巨大在线播放| 欧美精品国产亚洲| 乱人伦中国视频| 中文字幕人妻熟女乱码| 国产成人aa在线观看| 精品熟女少妇av免费看| 91在线精品国自产拍蜜月| 另类亚洲欧美激情| 亚洲精品乱久久久久久| 中文字幕免费在线视频6| 久久99蜜桃精品久久| 最近2019中文字幕mv第一页| 肉色欧美久久久久久久蜜桃| 日日爽夜夜爽网站| 亚洲,欧美精品.| 精品久久国产蜜桃| 精品人妻一区二区三区麻豆| 哪个播放器可以免费观看大片| 午夜激情av网站| 在线天堂最新版资源| 亚洲精品av麻豆狂野| 18+在线观看网站| 国产激情久久老熟女| 精品人妻一区二区三区麻豆| 一区在线观看完整版| 一边摸一边做爽爽视频免费| 免费观看在线日韩| 午夜91福利影院| 在线天堂中文资源库| av天堂久久9| 男女无遮挡免费网站观看| 国产国语露脸激情在线看| 欧美丝袜亚洲另类| 国产精品 国内视频| 老司机影院成人| 在线天堂最新版资源| 国产精品久久久久久av不卡| 亚洲国产精品999| 少妇的丰满在线观看| 亚洲国产精品国产精品| 国产国语露脸激情在线看| 欧美成人精品欧美一级黄| 汤姆久久久久久久影院中文字幕| 久久久久国产精品人妻一区二区| 国产免费现黄频在线看| 亚洲,一卡二卡三卡| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 免费看光身美女| 国国产精品蜜臀av免费| 国产日韩欧美在线精品| 色吧在线观看| 国产亚洲精品第一综合不卡 | 性色av一级| 亚洲精品日韩在线中文字幕| 国产黄色视频一区二区在线观看| 亚洲欧美一区二区三区黑人 | 我要看黄色一级片免费的| 国产av精品麻豆| 99九九在线精品视频| 国产成人精品在线电影| 校园人妻丝袜中文字幕| 国产伦理片在线播放av一区| 91久久精品国产一区二区三区| 国产一区有黄有色的免费视频| 国产成人免费观看mmmm| 老熟女久久久| 国产一区二区在线观看av| 老熟女久久久| 亚洲美女视频黄频| 亚洲av在线观看美女高潮| 制服诱惑二区| 日韩一区二区三区影片| 最新的欧美精品一区二区| 人妻系列 视频| 精品福利永久在线观看| 人妻少妇偷人精品九色| 中文字幕人妻丝袜制服| 美国免费a级毛片| 中文字幕另类日韩欧美亚洲嫩草| 午夜福利影视在线免费观看| 中文字幕人妻丝袜制服| 中文字幕精品免费在线观看视频 | 国精品久久久久久国模美| 五月开心婷婷网| 毛片一级片免费看久久久久| 色94色欧美一区二区| 免费观看无遮挡的男女| 91精品三级在线观看| 久久99热这里只频精品6学生| 在线天堂中文资源库| 久久久久久久精品精品| 国产在线视频一区二区| 一二三四中文在线观看免费高清| 五月开心婷婷网| 国产免费一区二区三区四区乱码| 久久久国产精品麻豆| 丰满少妇做爰视频| 久久久久久久久久人人人人人人| 国产免费一区二区三区四区乱码| 大陆偷拍与自拍| 高清欧美精品videossex| 亚洲欧美一区二区三区黑人 | 卡戴珊不雅视频在线播放| 不卡视频在线观看欧美| 亚洲成国产人片在线观看| 亚洲人与动物交配视频| 18+在线观看网站| 街头女战士在线观看网站| 看免费av毛片| 亚洲精品国产av蜜桃| 久久久久精品人妻al黑| 国产极品天堂在线| 啦啦啦中文免费视频观看日本| 免费大片黄手机在线观看| 成人漫画全彩无遮挡| 久久av网站| 一区二区三区四区激情视频| 国产成人精品在线电影| 午夜老司机福利剧场| 亚洲一码二码三码区别大吗| 欧美精品亚洲一区二区| 午夜福利网站1000一区二区三区| 狂野欧美激情性bbbbbb| 久久99蜜桃精品久久| 狠狠精品人妻久久久久久综合| 精品久久久久久电影网| 久久久久久久亚洲中文字幕| 亚洲av福利一区| 久热这里只有精品99| 青春草亚洲视频在线观看| 亚洲国产日韩一区二区| 男男h啪啪无遮挡| 成人亚洲欧美一区二区av| 欧美人与性动交α欧美软件 | 99久久精品国产国产毛片| 边亲边吃奶的免费视频| 少妇熟女欧美另类| 国产极品粉嫩免费观看在线| 国产在线一区二区三区精| xxx大片免费视频|