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

    基于ANUSPLIN模型的不同分辨率DEM場景下江西省降雨量插值對比分析

    2021-05-11 03:59:12李國文凍芳芳程雪蓉曾憲隆
    江西水利科技 2021年2期
    關(guān)鍵詞:雨量降雨量插值

    劉 林,李國文,凍芳芳,劉 貢,程雪蓉,曾憲隆

    (1.贛江上游水文水資源監(jiān)測中心,江西 贛州,341000;2.江西省水文監(jiān)測中心,江西 南昌,330002;3.信江饒河水文水資源監(jiān)測中心,江西 上饒,334000)

    0 引言

    江西省是全國多雨的省區(qū)之一,省內(nèi)暴雨區(qū)屬長江流域五大暴雨區(qū)之首。境內(nèi)地勢南高北低,邊緣群山環(huán)繞,中部丘陵起伏,北部平原坦蕩,四周漸次向鄱陽湖區(qū)傾斜。全省多年平均降雨量為1638mm,降雨時空分布不均,汛期(4~9月)降雨量占全年總降雨量的68%,其中主汛期(4~6月)占整個汛期的67%;由于地形的差異和距離海洋遠(yuǎn)近的不同,省內(nèi)有懷玉山區(qū)、武夷山區(qū)、九嶺山區(qū)、羅霄山區(qū)等4個多雨區(qū)。雨量站觀測是降雨量觀測的主要手段,據(jù)統(tǒng)計,至2015年6月,全省有雨量觀測項目4 545站,平均雨量站網(wǎng)密度36.7km2/站,基本能滿足應(yīng)用要求,但對于中小河流和山洪災(zāi)害防治,需要更高的站網(wǎng)密度,降雨量插值是經(jīng)濟可行的解決方法。空間內(nèi)插就是基于某種數(shù)理方法采用已知的氣象點數(shù)據(jù)來推求未知的氣象點數(shù)據(jù),從而把矢量數(shù)據(jù)轉(zhuǎn)換成柵格數(shù)據(jù),氣象上常用的插值方法有克里金插值法(Kriging)、反距離加權(quán)法(Inverse Distance Weighting,IDW)和樣條函數(shù)法(Spline)等[1]。

    ANUSPLIN模型是澳大利亞科學(xué)家Hutchinson基于薄盤樣條理論編寫的針對氣候數(shù)據(jù)曲面擬合的專用氣象插值軟件,對于時間序列的氣象數(shù)據(jù)尤其適合[2]。劉志紅等[3]2008年詳細(xì)介紹了ANUSPLIN模型及其插值應(yīng)用。錢永蘭等[4]2010年用ANUSPLIN模型對1961~2006年逐日氣溫、降雨數(shù)據(jù)進行插值,并與反向距離權(quán)重法和普通克呂格法的結(jié)果進行對比,結(jié)果表明ANUSPLIN模型的插值誤差最小。王軍邦等[5]2017年用2000~2012年全國日均溫、日降雨數(shù)據(jù)由ANUSPLIN模型插值生成了全國溫度與降雨的1km網(wǎng)格數(shù)據(jù)集。孟清等[6]2019年用1959~2015年秦嶺地區(qū)31個氣象站降雨資料由ANUSPLIN模型生成降雨柵格數(shù)據(jù),研究降雨時空變化。

    地形對降雨的產(chǎn)生與分布具有重要影響,如迎風(fēng)坡具有強迫氣流抬升的作用,暖濕氣流受動力冷卻而成云致雨,從而使降雨量增加[7]。許多學(xué)者研究表明,在ANUSPLIN模型中考慮高程因素影響,有利于提高降雨插值精度。朱求安等[8]2009年采用反距離加權(quán)、普通克里格、ANUSPLIN模型3種方法進行酸雨空間插值模擬,結(jié)果顯示,前兩種方法僅能反映酸雨水平方向上的分布特點,而ANUSPLIN法將地形作為協(xié)變量因子,插值結(jié)果更合理地體現(xiàn)了酸雨時空分布特征。錢永蘭等[4]2010年用ANUSPLIN模型對1961~2006年逐日氣溫、降雨量進行插值,認(rèn)為適當(dāng)增加站點數(shù)量和提高DEM精度可進一步提高其插值精度。李月晨等[9]2014年探討了多種插值方法的適用性和不足之處,認(rèn)為考慮地形等更多的相關(guān)因子以提高模擬精度是未來發(fā)展趨勢。趙冠華等[10]2016年研究了DEM的不確定性對溫度插值精度的影響。陸福志等[11]2019年采用ANUSPLIN模型,以經(jīng)度、緯度、海拔為獨立變量,通過空間插值,建立秦嶺—大巴山氣候與降雨格點數(shù)據(jù)集,結(jié)果與流行的WorldClim 2.0氣候格點數(shù)據(jù)集具有一致性,且比其更精確、分辨率更高、細(xì)節(jié)更突出。

    研究[4,8,9,10,11,12]表明,ANUSPLIN 模型在降雨量插值上具有優(yōu)勢,但選擇不同分辨率數(shù)字高程模型(DEM)數(shù)據(jù)會影響其插值精度,因此,本文選擇ANUSPLIN模型應(yīng)用到江西省降雨量插值中,主要探究不同分辨率DEM場景下ANUSPLIN模型降雨量插值精度。

    1 研究內(nèi)容及數(shù)據(jù)情況

    本文把ANUSPLIN模型應(yīng)用到江西省降雨量插值中,分別選用500mDEM、1kmDEM、2kmDEM三種高程數(shù)據(jù)分辨率場景,根據(jù)江西省南高北低漸次傾斜的地形特點,將贛州市歸為山丘區(qū),吉安市歸為丘陵區(qū),南昌市歸為平原區(qū),采用均方根誤差、確定性系數(shù)兩個評價指標(biāo),對比分析三種分辨率DEM數(shù)據(jù)場景下三種地形代表區(qū)降雨量插值的精度,選擇適合江西省降雨量插值的DEM數(shù)據(jù),為江西省雨情自動測報系統(tǒng)雨量質(zhì)量提升工作奠定基礎(chǔ)。

    目前可從相關(guān)機構(gòu)免費獲取30mDEM高程數(shù)據(jù),綜合考慮目前氣象上精細(xì)化降雨產(chǎn)品網(wǎng)格尺度,本文選取由30mDEM重采樣得到的500mDEM、1kmDEM、2kmDEM三種地形分辨率場景進行分析。

    江西省通過雨量密度研究和面雨量精度比較,選擇了分布均勻、資料系列可靠完整的1 085個基本雨量站統(tǒng)一計算省、市、縣三級行政區(qū)面雨量。1 085雨量站中,贛州、吉安、南昌市雨量站數(shù)如表1所示,從表中可以看出三地市基本站密度差不多,均為155km2/站左右,用以進行降雨表面分析。按照2 000km2/站,并結(jié)合均勻分布的原則,分別在贛州、吉安、南昌市抽取20、13、4站,用以對比模型分析結(jié)果。選擇2014年至2019年三地市同期降雨資料,其中贛州市選擇了30場降雨,吉安市選擇了30場降雨,南昌市選擇了29場降雨,受篇幅所限,在此不列出。本文基于此進行分析計算,并將ANUSPLIN模型分析得到的雨量數(shù)據(jù)與實測雨量數(shù)據(jù)進行誤差對比,確定誤差程度小的DEM場景分辨率,作為適合江西省的數(shù)字高程分辨率。

    表1 抽取站點數(shù)量統(tǒng)計表

    2 研究方法

    2.1 ANUSPLIN模型基本原理

    ANUSPLIN基于普通薄盤和局部薄盤樣條函數(shù)插值理論。局部薄盤光滑樣條法是對薄盤光滑樣條原型的擴展[2],除普通的樣條自變量外允許引入線性協(xié)變量子模型,如溫度與海拔、降雨與海岸線之間的相關(guān)關(guān)系。

    局部薄盤光滑樣條的理論統(tǒng)計模型為:

    式中:Zi是位于空間i點的因變量;xi為d維樣條獨立變量,f為需要估算的關(guān)于 xi的未知光滑函數(shù);yi為p維獨立協(xié)變量;b為yi的p維系數(shù);ei為具有期望值為0且方差為 wiσ2的自變量隨機誤差,其中wi為作為權(quán)重的已知局部相對變異系數(shù),σ2為誤差方差,在所有數(shù)據(jù)點上為常數(shù),但通常未知。

    劉志紅等[3]闡述了SPLINA模塊中18個待選薄盤光滑樣條函數(shù)。本模型計算分析過程中,采用三變量薄盤光滑樣條函數(shù)TVPTPS2,即將經(jīng)度、緯度、高程作為獨立變量,數(shù)量為3,協(xié)變量為0,樣條次數(shù)為默認(rèn)值2,調(diào)整經(jīng)緯度范圍值,其他模型參數(shù)均為默認(rèn)值,進行模型運算。

    2.2 降雨量插值精度評價指標(biāo)

    2.2.1 均方根誤差

    為了驗證不同分辨率的DE M高程對降雨數(shù)據(jù)分析值的影響,通過對雨量站點觀測數(shù)據(jù)進行驗證,通過計算均方根誤差(RMSE)來評價不同DEM高程分辨率下的降雨預(yù)測值的精度,公式為:

    式中:Zx和Zy分別代表降雨量實測值,mm和預(yù)測值,mm;n為計算序列長度。RMSE的值(mm)越小代表降雨預(yù)測值的精度越高。

    2.2.2 確定性系數(shù)

    在此次研究中將觀測值與預(yù)測值進行相關(guān)性分析,確定性系數(shù)R2表示一元多項式回歸方程擬合度的高低。

    關(guān)于確定性系數(shù)R2的評價標(biāo)準(zhǔn)如表2所示。R2的整體取值分布在0~1之間,當(dāng)它的值越大,表示模擬程度越好。當(dāng)它達(dá)到0.9以上時,便認(rèn)為模擬的程度屬于優(yōu)秀程度。

    表2 確定性系數(shù)評價標(biāo)準(zhǔn)

    3 結(jié)果對比與分析

    500mDEM、1kmDEM、2kmDEM三種分辨率場景下,贛州(山丘區(qū))、吉安(丘陵區(qū))、南昌(平原區(qū))三市所抽基本雨量站各同期場次降雨的模型分析雨量值與實測值散點圖如圖1~9所示,其統(tǒng)計特征如表3~5所示。對比圖、表中數(shù)據(jù)可得:

    表3 不同地區(qū)不同DEM分辨率場景下降雨量插值的均方根誤差RMSE對比表 mm

    表4 不同地區(qū)不同DEM分辨率場景下降雨量插值的確定性系數(shù)R2對比表 mm

    表5 不同地區(qū)不同DEM分辨率場景下降雨量插值的計算時間t對比表 min

    圖1 500mDEM贛州市抽站雨量分析值與實測值散點圖

    圖2 1kmDEM贛州市抽站雨量分析值與實測值散點圖

    圖3 2kmDEM贛州市抽站雨量分析值與實測值散點圖

    圖4 500mDEM吉安市抽站雨量分析值與實測值散點圖

    圖5 1kmDEM吉安市抽站雨量分析值與實測值散點圖

    圖6 2kmDEM吉安市抽站雨量分析值與實測值散點圖

    圖7 500mDEM南昌市抽站雨量分析值與實測值散點圖

    圖8 1kmDEM南昌市抽站雨量分析值與實測值散點圖

    圖9 2kmDEM南昌市抽站雨量分析值與實測值散點圖

    ANUSPLIN模型在江西省降雨量插值應(yīng)用中取得了較好效果。三種場景下降雨量插值的均方根誤差在11.3~21.9mm之間;確定性系數(shù)在0.85~0.95之間,山丘區(qū)為良好水平,丘陵區(qū)和平原區(qū)為優(yōu)秀水平。

    在相同DEM分辨率下的不同地區(qū)中:均方根誤差(RMSE)精度方面,平原區(qū)的插值數(shù)據(jù)精度較好,丘陵區(qū)插值數(shù)據(jù)精度其次,山丘區(qū)插值數(shù)據(jù)精度較差。在確定性系數(shù)(R2)方面,平原區(qū)擬合效果較好,丘陵區(qū)擬合效果其次,山丘區(qū)擬合效果稍差。綜上,在相同DEM分辨率的模擬場景下,模擬精度排序為平原區(qū)>丘陵區(qū)>山丘區(qū)。

    在相同地區(qū)的不同DEM分辨率中:均方根誤差(RMSE)精度方面,在所有地區(qū)中,分辨率為500m的DEM作為高程輸入數(shù)據(jù)計算的降雨量插值精度較高,但試驗時間為12至15min,明顯長于1km分辨率和2km分辨率的場景;分辨率為1km的DEM作為高程輸入數(shù)據(jù)計算的降雨量插值精度略低于500m分辨率的插值精度,計算時間較短,為10至13min;分辨率為2km的DEM作為高程輸入數(shù)據(jù)計算的降雨量插值精度較差,計算時間與1km分辨率場景的計算時間相近。確定性系數(shù)(R2)方面,在所有地區(qū)中,500m分辨率的DEM高程應(yīng)用場景和1km分辨率的DEM高程應(yīng)用場景下,兩者的擬合程度相近,500m分辨率的DEM高程應(yīng)用場景的模擬數(shù)值擬合程度稍好于1km分辨率的擬合程度;2km分辨率的DEM高程應(yīng)用場景的模擬數(shù)值擬合程度稍差。綜上,在相同地區(qū)的不同DEM分辨率中,模擬精度排序為500m分辨率≥1km分辨率>2km分辨率。

    綜合上述分析,在計算江西省整個省份的雨量插值的實際應(yīng)用過程中,應(yīng)當(dāng)在保證數(shù)據(jù)質(zhì)量的同時,還需滿足計算時長較短以提高工作效率的前提下,結(jié)合地形影響和實際計算時長影響的綜合因素,選用分辨率為1km的DEM高程值的應(yīng)用場景最為適合。

    4 結(jié) 論

    本文利用ANUSPLIN模型,在500m、1km、2km三種分辨率DEM場景下分別對贛州市(山丘區(qū))、吉安市(丘陵區(qū))、南昌市(平原區(qū))進行了降雨插值分析,對比分析了三種場景不同地形代表區(qū)降雨插值的精度,得出以下結(jié)論:

    (1)ANUSPLIN模型對江西省降雨量插值有較好的適應(yīng)性。三種場景下降雨量插值的均方根誤差在11.3~21.9mm之間;確定性系數(shù)在0.85~0.95之間,均為良好以上水平。

    (2)在相同DEM分辨率的模擬場景下,ANUSPLIN模型降雨量插值精度排序為平原區(qū)>丘陵區(qū)>山丘區(qū)。在相同地區(qū)的不同DEM分辨率模擬場景下,ANUSPLIN模型降雨量插值精度排序為500m分辨率≥1km分辨率>2km分辨率。

    (3)在江西省降雨量插值的實際應(yīng)用過程中,綜合考慮地形影響和模型計算效率等因素,選用分辨率為1km的DEM高程值的應(yīng)用場景較為合適。

    本文將ANUSPLIN模型應(yīng)用到江西省降雨量插值中,確定了合適分辨率的DEM應(yīng)用場景。下一步可以將ANUSPLIN模型應(yīng)用到江西省雨量遙測站點數(shù)據(jù)質(zhì)量提升工作中。

    猜你喜歡
    雨量降雨量插值
    寧夏紅柳溝流域水沙變化及產(chǎn)沙分析
    降雨量與面積的關(guān)系
    基于小波去噪的稱重雨量數(shù)據(jù)分析
    基于Sinc插值與相關(guān)譜的縱橫波速度比掃描方法
    SL—1 型雨量傳感器故障分析排除和維護
    西藏科技(2016年5期)2016-09-26 12:16:40
    一種改進FFT多譜線插值諧波分析方法
    基于四項最低旁瓣Nuttall窗的插值FFT諧波分析
    洞庭湖區(qū)降雨特性分析
    Blackman-Harris窗的插值FFT諧波分析與應(yīng)用
    羅甸縣各鄉(xiāng)鎮(zhèn)實測降雨量分析及應(yīng)用研究
    成人国产一区最新在线观看| 久久久国产成人精品二区 | 好男人电影高清在线观看| 老汉色∧v一级毛片| 一级a爱视频在线免费观看| 精品久久久精品久久久| 9色porny在线观看| 人人妻人人添人人爽欧美一区卜| 男人的好看免费观看在线视频 | 久久久精品区二区三区| 老司机深夜福利视频在线观看| 麻豆国产av国片精品| 精品福利永久在线观看| 亚洲成av片中文字幕在线观看| av在线播放免费不卡| tube8黄色片| 午夜福利,免费看| 亚洲av美国av| a级片在线免费高清观看视频| 国产午夜精品久久久久久| 黄色视频不卡| 亚洲av日韩精品久久久久久密| 亚洲性夜色夜夜综合| 十八禁网站免费在线| 午夜亚洲福利在线播放| 亚洲精华国产精华精| 午夜福利在线观看吧| 国产亚洲精品久久久久5区| 一夜夜www| 香蕉丝袜av| 欧美老熟妇乱子伦牲交| 免费黄频网站在线观看国产| 亚洲精品中文字幕在线视频| 成年人午夜在线观看视频| 无限看片的www在线观看| 高潮久久久久久久久久久不卡| 高清av免费在线| 国产精品欧美亚洲77777| 久久久国产成人精品二区 | 无人区码免费观看不卡| 成人国语在线视频| 精品国产乱子伦一区二区三区| 欧美久久黑人一区二区| 免费一级毛片在线播放高清视频 | 国产欧美亚洲国产| 国产一区二区三区综合在线观看| 亚洲视频免费观看视频| www.999成人在线观看| 欧美黑人精品巨大| 亚洲片人在线观看| 亚洲精品一二三| 国产男女超爽视频在线观看| 亚洲欧美一区二区三区黑人| 视频在线观看一区二区三区| 国产色视频综合| 久久午夜亚洲精品久久| 男女下面插进去视频免费观看| 超色免费av| 久久午夜亚洲精品久久| 国产极品粉嫩免费观看在线| 777米奇影视久久| 高潮久久久久久久久久久不卡| 色播在线永久视频| 在线国产一区二区在线| 变态另类成人亚洲欧美熟女 | 亚洲五月天丁香| 最新的欧美精品一区二区| 麻豆国产av国片精品| 久久久国产成人精品二区 | 啦啦啦在线免费观看视频4| 亚洲 欧美一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 日韩欧美国产一区二区入口| 精品一区二区三卡| 99热只有精品国产| 国产淫语在线视频| 亚洲国产看品久久| 美女午夜性视频免费| 一a级毛片在线观看| 久久天堂一区二区三区四区| 欧美激情 高清一区二区三区| 十八禁网站免费在线| 黄片小视频在线播放| 十八禁网站免费在线| 午夜福利免费观看在线| 国产精品久久久久久人妻精品电影| 黄片大片在线免费观看| 一级作爱视频免费观看| 免费在线观看日本一区| 午夜福利免费观看在线| 极品教师在线免费播放| 久久久久国内视频| 美女福利国产在线| 一级a爱片免费观看的视频| 18禁美女被吸乳视频| 久久久久国产一级毛片高清牌| 超碰成人久久| 日韩欧美一区视频在线观看| 亚洲国产欧美一区二区综合| 国产精品永久免费网站| 丰满饥渴人妻一区二区三| 亚洲久久久国产精品| 男人操女人黄网站| 精品久久久精品久久久| tube8黄色片| 亚洲国产精品一区二区三区在线| 久久午夜综合久久蜜桃| ponron亚洲| 亚洲精品国产精品久久久不卡| 女警被强在线播放| 老司机深夜福利视频在线观看| 巨乳人妻的诱惑在线观看| 男女免费视频国产| 日韩欧美在线二视频 | 欧美成狂野欧美在线观看| 怎么达到女性高潮| 村上凉子中文字幕在线| 狠狠婷婷综合久久久久久88av| 久久草成人影院| 美女国产高潮福利片在线看| 国产欧美日韩精品亚洲av| 免费在线观看黄色视频的| 中文字幕高清在线视频| 日本五十路高清| 亚洲性夜色夜夜综合| 色综合欧美亚洲国产小说| 少妇 在线观看| 国产av一区二区精品久久| 又大又爽又粗| 免费在线观看日本一区| 巨乳人妻的诱惑在线观看| 母亲3免费完整高清在线观看| 久久影院123| 不卡一级毛片| 久久国产精品大桥未久av| 777米奇影视久久| 亚洲专区国产一区二区| 国产三级黄色录像| 国产精品秋霞免费鲁丝片| 午夜福利影视在线免费观看| 午夜免费观看网址| 精品国产一区二区久久| 不卡av一区二区三区| 国产高清激情床上av| 欧美亚洲 丝袜 人妻 在线| 脱女人内裤的视频| 日韩成人在线观看一区二区三区| 欧美乱妇无乱码| 王馨瑶露胸无遮挡在线观看| 国产精品国产av在线观看| 成年人黄色毛片网站| 日日摸夜夜添夜夜添小说| 亚洲美女黄片视频| 大香蕉久久成人网| 亚洲专区字幕在线| 在线免费观看的www视频| av福利片在线| 国产精品99久久99久久久不卡| 高清视频免费观看一区二区| 色婷婷av一区二区三区视频| 国产精品一区二区精品视频观看| 亚洲精品粉嫩美女一区| 777米奇影视久久| 成人特级黄色片久久久久久久| 国产精品九九99| 成人国语在线视频| 老司机靠b影院| 成在线人永久免费视频| 精品久久久久久久久久免费视频 | 精品久久久久久,| 亚洲av熟女| 久久午夜综合久久蜜桃| 色在线成人网| 亚洲一区二区三区欧美精品| 韩国av一区二区三区四区| 国产淫语在线视频| 视频区图区小说| 午夜91福利影院| av欧美777| 欧美乱妇无乱码| 美国免费a级毛片| 午夜免费观看网址| 精品欧美一区二区三区在线| 黄色 视频免费看| 亚洲九九香蕉| 精品国内亚洲2022精品成人 | 在线观看日韩欧美| 精品福利观看| 国产男女超爽视频在线观看| 夜夜躁狠狠躁天天躁| 色在线成人网| 操出白浆在线播放| 免费日韩欧美在线观看| 五月开心婷婷网| 欧美一级毛片孕妇| 男人舔女人的私密视频| 18在线观看网站| 在线观看免费高清a一片| 巨乳人妻的诱惑在线观看| 亚洲熟女精品中文字幕| 欧美日韩精品网址| 久久人妻福利社区极品人妻图片| xxxhd国产人妻xxx| 精品国产美女av久久久久小说| 午夜老司机福利片| 日韩中文字幕欧美一区二区| 男女免费视频国产| 精品国产国语对白av| 精品国产美女av久久久久小说| 久久这里只有精品19| 欧美日韩av久久| 亚洲熟妇熟女久久| 国产高清videossex| 亚洲 国产 在线| 大型av网站在线播放| 人妻丰满熟妇av一区二区三区 | 视频区图区小说| 水蜜桃什么品种好| 午夜福利在线免费观看网站| 黑人欧美特级aaaaaa片| 成年动漫av网址| 国产成人av教育| 国产精品免费视频内射| 亚洲av第一区精品v没综合| 黑人欧美特级aaaaaa片| 精品人妻在线不人妻| 亚洲人成电影免费在线| 麻豆成人av在线观看| 国产精品自产拍在线观看55亚洲 | 国产亚洲欧美在线一区二区| 国产精品一区二区在线不卡| tube8黄色片| 男人操女人黄网站| 国产精品国产高清国产av | 在线十欧美十亚洲十日本专区| 午夜福利欧美成人| 男人的好看免费观看在线视频 | 韩国精品一区二区三区| 久久性视频一级片| 国产精品1区2区在线观看. | 免费在线观看日本一区| 18禁国产床啪视频网站| 夜夜夜夜夜久久久久| 黄网站色视频无遮挡免费观看| 亚洲九九香蕉| 中文字幕高清在线视频| 成年人免费黄色播放视频| 美女扒开内裤让男人捅视频| 国产亚洲精品第一综合不卡| 成人手机av| 国产免费男女视频| 国产午夜精品久久久久久| 国产精品欧美亚洲77777| 新久久久久国产一级毛片| 桃红色精品国产亚洲av| 成年版毛片免费区| 亚洲欧美一区二区三区黑人| 亚洲一区二区三区不卡视频| 捣出白浆h1v1| svipshipincom国产片| 国产精品影院久久| 免费在线观看黄色视频的| 侵犯人妻中文字幕一二三四区| 久久国产精品大桥未久av| 91老司机精品| 一进一出好大好爽视频| 午夜成年电影在线免费观看| 亚洲少妇的诱惑av| 欧美日韩一级在线毛片| 999精品在线视频| 国产精品久久久久久人妻精品电影| 王馨瑶露胸无遮挡在线观看| 亚洲成av片中文字幕在线观看| 美女扒开内裤让男人捅视频| 午夜福利视频在线观看免费| 精品国产一区二区久久| 天天躁夜夜躁狠狠躁躁| 别揉我奶头~嗯~啊~动态视频| 在线观看66精品国产| 久久人妻熟女aⅴ| 成人精品一区二区免费| 日本黄色日本黄色录像| 欧美久久黑人一区二区| 宅男免费午夜| 黄色丝袜av网址大全| 精品久久久久久,| 一进一出抽搐动态| 三上悠亚av全集在线观看| a在线观看视频网站| 色婷婷久久久亚洲欧美| 亚洲精品中文字幕一二三四区| 国产成人免费无遮挡视频| 精品乱码久久久久久99久播| 丰满人妻熟妇乱又伦精品不卡| 国产一区有黄有色的免费视频| 黄色片一级片一级黄色片| 热re99久久精品国产66热6| 男女下面插进去视频免费观看| 午夜91福利影院| 国产精品九九99| 精品第一国产精品| 在线免费观看的www视频| 午夜久久久在线观看| av中文乱码字幕在线| 欧美日韩一级在线毛片| 欧美日韩亚洲高清精品| tocl精华| 亚洲精品在线美女| 高清毛片免费观看视频网站 | 亚洲精品国产区一区二| 桃红色精品国产亚洲av| 少妇 在线观看| 久久久久久久精品吃奶| 曰老女人黄片| 国产成人精品无人区| 欧美日韩国产mv在线观看视频| 老汉色av国产亚洲站长工具| 亚洲在线自拍视频| 精品视频人人做人人爽| 自拍欧美九色日韩亚洲蝌蚪91| 韩国精品一区二区三区| 精品无人区乱码1区二区| 老汉色∧v一级毛片| 大香蕉久久成人网| 亚洲成人免费电影在线观看| 亚洲五月婷婷丁香| 亚洲熟女毛片儿| 中亚洲国语对白在线视频| 热99国产精品久久久久久7| 久久久国产精品麻豆| 人妻丰满熟妇av一区二区三区 | 熟女少妇亚洲综合色aaa.| 亚洲美女黄片视频| 狂野欧美激情性xxxx| 精品亚洲成a人片在线观看| 免费在线观看黄色视频的| 999精品在线视频| 国产精华一区二区三区| 一区二区三区激情视频| 国产精品偷伦视频观看了| 女同久久另类99精品国产91| 视频在线观看一区二区三区| 国产精品亚洲一级av第二区| 最新的欧美精品一区二区| 久热爱精品视频在线9| 天天操日日干夜夜撸| 在线十欧美十亚洲十日本专区| 中文字幕另类日韩欧美亚洲嫩草| 99国产精品一区二区三区| 日韩欧美免费精品| 国产精品乱码一区二三区的特点 | 久久久国产精品麻豆| 啦啦啦在线免费观看视频4| 国产一区二区三区综合在线观看| 日日摸夜夜添夜夜添小说| 天堂√8在线中文| 欧美久久黑人一区二区| 久久久久精品人妻al黑| 99精品在免费线老司机午夜| 黄色视频不卡| 啦啦啦在线免费观看视频4| 日本黄色日本黄色录像| 在线观看免费高清a一片| 美女高潮到喷水免费观看| 人妻丰满熟妇av一区二区三区 | 亚洲精品国产色婷婷电影| 午夜福利一区二区在线看| 国产精品久久电影中文字幕 | 国产精品亚洲av一区麻豆| 99国产极品粉嫩在线观看| bbb黄色大片| 丝袜在线中文字幕| 交换朋友夫妻互换小说| 亚洲成人手机| 精品国产一区二区三区久久久樱花| 成人18禁高潮啪啪吃奶动态图| 涩涩av久久男人的天堂| aaaaa片日本免费| av天堂在线播放| 亚洲成a人片在线一区二区| 精品亚洲成a人片在线观看| 女人精品久久久久毛片| 精品一区二区三卡| 高清在线国产一区| 日本精品一区二区三区蜜桃| 免费在线观看日本一区| 国产精品一区二区在线观看99| 国产成人免费无遮挡视频| 久久精品人人爽人人爽视色| 咕卡用的链子| 久久99一区二区三区| 色婷婷久久久亚洲欧美| 午夜福利在线观看吧| 午夜激情av网站| 成在线人永久免费视频| 婷婷丁香在线五月| 悠悠久久av| 精品福利永久在线观看| 欧美乱码精品一区二区三区| 一本一本久久a久久精品综合妖精| 精品午夜福利视频在线观看一区| 欧美亚洲日本最大视频资源| 亚洲五月天丁香| 又紧又爽又黄一区二区| 香蕉国产在线看| 久久婷婷成人综合色麻豆| 亚洲人成电影免费在线| 男女午夜视频在线观看| 美女扒开内裤让男人捅视频| 亚洲第一青青草原| 脱女人内裤的视频| 国产一区二区三区在线臀色熟女 | 婷婷成人精品国产| 久久精品91无色码中文字幕| 亚洲一区中文字幕在线| 欧美激情高清一区二区三区| 怎么达到女性高潮| 精品国产一区二区三区四区第35| 黄色片一级片一级黄色片| 久久国产精品影院| 色94色欧美一区二区| 亚洲午夜理论影院| 亚洲性夜色夜夜综合| 91精品国产国语对白视频| 国产一卡二卡三卡精品| 国产精品久久久久久人妻精品电影| 日韩视频一区二区在线观看| 欧美大码av| 日本五十路高清| 午夜福利在线观看吧| 成年版毛片免费区| 狠狠婷婷综合久久久久久88av| 久久青草综合色| 国产深夜福利视频在线观看| 亚洲国产欧美一区二区综合| 精品福利观看| 亚洲欧美一区二区三区黑人| 午夜福利影视在线免费观看| 午夜福利,免费看| 久久久久精品国产欧美久久久| 亚洲美女黄片视频| 日本一区二区免费在线视频| 日本黄色视频三级网站网址 | xxx96com| 免费黄频网站在线观看国产| 1024视频免费在线观看| 欧美精品啪啪一区二区三区| 亚洲九九香蕉| 亚洲欧洲精品一区二区精品久久久| 日韩中文字幕欧美一区二区| 亚洲中文av在线| 久久国产亚洲av麻豆专区| 久久香蕉国产精品| 国产人伦9x9x在线观看| 日韩精品免费视频一区二区三区| 免费在线观看亚洲国产| 天堂俺去俺来也www色官网| 一夜夜www| 亚洲欧美日韩高清在线视频| 国产精品免费一区二区三区在线 | 欧美人与性动交α欧美软件| 国产蜜桃级精品一区二区三区 | 久久久久久久国产电影| 日韩人妻精品一区2区三区| 一本大道久久a久久精品| 叶爱在线成人免费视频播放| 久久 成人 亚洲| 日韩熟女老妇一区二区性免费视频| 亚洲欧美日韩高清在线视频| 日韩欧美一区视频在线观看| 久久久国产成人精品二区 | 久久精品aⅴ一区二区三区四区| 十分钟在线观看高清视频www| 男女午夜视频在线观看| 日本撒尿小便嘘嘘汇集6| 国产成人欧美| 人人妻人人澡人人爽人人夜夜| 欧美乱码精品一区二区三区| 国产精品香港三级国产av潘金莲| 色婷婷av一区二区三区视频| 午夜视频精品福利| 69精品国产乱码久久久| 又紧又爽又黄一区二区| 男女下面插进去视频免费观看| 一边摸一边抽搐一进一出视频| 日韩欧美免费精品| 亚洲午夜精品一区,二区,三区| 亚洲aⅴ乱码一区二区在线播放 | 国产单亲对白刺激| 国产亚洲精品第一综合不卡| 日本撒尿小便嘘嘘汇集6| av线在线观看网站| 欧美性长视频在线观看| 中文亚洲av片在线观看爽 | 日韩有码中文字幕| 在线观看www视频免费| 18在线观看网站| 国产精品亚洲一级av第二区| 国产精华一区二区三区| 老汉色∧v一级毛片| 亚洲人成电影免费在线| 十八禁高潮呻吟视频| 女人被狂操c到高潮| 精品久久久久久,| 首页视频小说图片口味搜索| 大型黄色视频在线免费观看| 这个男人来自地球电影免费观看| 一边摸一边做爽爽视频免费| 满18在线观看网站| 亚洲国产精品合色在线| 久久久久国内视频| 国产亚洲av高清不卡| 最新的欧美精品一区二区| 亚洲中文字幕日韩| 国产精品久久久久久人妻精品电影| 岛国在线观看网站| 亚洲黑人精品在线| 久热爱精品视频在线9| 制服诱惑二区| 成人av一区二区三区在线看| 亚洲av欧美aⅴ国产| 97人妻天天添夜夜摸| 国产成人影院久久av| 久久久久久免费高清国产稀缺| 精品熟女少妇八av免费久了| 少妇粗大呻吟视频| 波多野结衣一区麻豆| 激情视频va一区二区三区| 国产黄色免费在线视频| 国产精品亚洲一级av第二区| 嫁个100分男人电影在线观看| 女同久久另类99精品国产91| 一区二区三区精品91| 久久久久久免费高清国产稀缺| 亚洲色图 男人天堂 中文字幕| 免费在线观看影片大全网站| 99热国产这里只有精品6| 黑人巨大精品欧美一区二区蜜桃| 日本vs欧美在线观看视频| 极品少妇高潮喷水抽搐| 免费不卡黄色视频| 亚洲欧美色中文字幕在线| 91麻豆av在线| 欧美黑人精品巨大| 91大片在线观看| 亚洲av美国av| 国产激情欧美一区二区| 人人妻人人澡人人爽人人夜夜| 韩国av一区二区三区四区| 啪啪无遮挡十八禁网站| 精品电影一区二区在线| www.熟女人妻精品国产| 成在线人永久免费视频| 精品一区二区三区视频在线观看免费 | 人人妻人人爽人人添夜夜欢视频| 老司机深夜福利视频在线观看| 男人操女人黄网站| 国产淫语在线视频| 成人精品一区二区免费| 18在线观看网站| 久久久国产欧美日韩av| 黄色a级毛片大全视频| 亚洲av第一区精品v没综合| 国产成人欧美| 日韩欧美在线二视频 | 国产乱人伦免费视频| 国产精品秋霞免费鲁丝片| 亚洲视频免费观看视频| a级毛片黄视频| 一级作爱视频免费观看| 两个人免费观看高清视频| netflix在线观看网站| 最新的欧美精品一区二区| 色在线成人网| 脱女人内裤的视频| 亚洲精华国产精华精| av中文乱码字幕在线| 国产成人av激情在线播放| 久久精品国产亚洲av高清一级| 午夜精品国产一区二区电影| 中出人妻视频一区二区| 十八禁网站免费在线| 成年版毛片免费区| 亚洲欧美激情综合另类| 亚洲成av片中文字幕在线观看| 久久久久国产精品人妻aⅴ院 | 免费少妇av软件| 日日夜夜操网爽| 免费看十八禁软件| 亚洲国产毛片av蜜桃av| 一级毛片精品| 国产精品久久视频播放| 色在线成人网| 丁香六月欧美| 又大又爽又粗| 欧美在线黄色| 色播在线永久视频| 午夜精品在线福利| 人妻 亚洲 视频| 久久精品亚洲精品国产色婷小说| 午夜精品久久久久久毛片777| 午夜福利一区二区在线看| 国产在视频线精品| 精品熟女少妇八av免费久了| 亚洲视频免费观看视频| 热99国产精品久久久久久7| 人人妻人人澡人人看| 亚洲视频免费观看视频| www日本在线高清视频| 99热国产这里只有精品6| 欧美老熟妇乱子伦牲交| a级毛片黄视频| 国产成人免费观看mmmm| 19禁男女啪啪无遮挡网站| 国产亚洲精品久久久久久毛片 |