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

    基于GAMLSS的雅礱江流域極端降水時空特性研究

    2019-04-22 08:32:00
    水力發(fā)電 2019年1期
    關(guān)鍵詞:平穩(wěn)性降雨量極值

    高 潔

    (水電水利規(guī)劃設(shè)計總院,北京100120)

    1 研究概述

    2008年Milly、Betancourt、Falkenmark等在《Science》

    發(fā)表論文[1],對基于平穩(wěn)序列的涉水管理提出質(zhì)疑。文章認為大量人為影響,正在改變降雨、蒸散發(fā)和河川徑流的均值和極值,氣候均值的微小變化可能導致極值的巨大改變,水文氣象領(lǐng)域在研究不確定性的同時,應將非平穩(wěn)性作為另一個重要關(guān)注點。丁一匯[2]通過全國726個地面站1951年至2002年觀測資料研究氣候均值和極值的時間變化趨勢和空間分布特點。研究發(fā)現(xiàn)全國年降水量沒有明顯變化趨勢,20世紀90年代以后極端降水量比例增大,空間分布各有差異。華北地區(qū)極端降水量減小,長江及長江以南地區(qū)極端降水量呈增加趨勢。任國玉等[3]通過1951年至2008年全國極端氣溫、降水指數(shù)研究發(fā)現(xiàn),全國1日和3日最大降水量在南方地區(qū)增加顯著。暴雨和極端強降水強度在西部地區(qū)、東南地區(qū)、長江中下游和華南沿海地區(qū)呈增加趨勢。李洋和楊赤[4]基于GEV(Generalized Extreme Value)分布的平穩(wěn)/非平穩(wěn)極值概率模型,分析全國范圍489個氣象站自建站至2013年極端最高氣溫、最低氣溫和日最大降水量序列,研究氣候極值“常態(tài)”(2年一遇)和“極端態(tài)”(50年一遇)的變化趨勢及空間分布。結(jié)果表明,日最大降水量在全國范圍內(nèi)普遍適用平穩(wěn)模型,如果氣候極值的“常態(tài)”平緩,“極端態(tài)”非平穩(wěn)且變率較大則可能導致劇烈的災害性天氣。根據(jù)IPCC第四次評估報告,我國的極端強降水值有增加趨勢,并呈現(xiàn)不同分區(qū)特點。

    QIN等[5]通過四川、云南、貴州、重慶、廣西共136萬km21960年至2009年116個氣象站降雨和氣溫數(shù)據(jù),分析極端氣溫和降雨指標的空間分布和時間變化。該區(qū)域內(nèi),最大1日降雨量在83站無顯著變化趨勢,26站顯著增加,7站顯著減小;有雨日指標在54站呈顯著減少趨勢,且主要分布在廣西。LIU等[6]通過旋轉(zhuǎn)經(jīng)驗正交函數(shù)法將西南地區(qū)基于降雨量分為5個子區(qū)域。選取9種極端降水指數(shù),通過MK趨勢檢驗,發(fā)現(xiàn)年總降水量呈減少趨勢,但最大1日和5日降雨量均顯著增加,雨型分布更為集中,且趨勢突變點分布在1962年和1996年。在空間分布上,云南和廣西總降雨量減少,但短歷時強降雨增加,加劇了洪水和干旱的風險。佘敦先等[7]通過淮河流域27個氣象站1960年至2009年逐日降水觀測資料,選取年最大降水量序列和超門限峰值序列,分析極端降水的時空變化趨勢,通過GEV、GP和GAMMA函數(shù)擬合降雨極值概率分布。榮艷淑等[8]選取淮河流域25個氣象站1951年至2007年降水資料,分析1日最大降水量序列、年最大過程降水量及最長持續(xù)降水日數(shù)序列的空間分布特征。GAO等[9-10]采用上海、浙江、福建、廣東、海南1960年至2012年86個氣象站日降雨數(shù)據(jù)集,選取日最大降水量和超過日降水量99%分位數(shù)的總降水量序列,通過MK趨勢檢驗,分析序列均值和方差。選取10個大尺度環(huán)流指標表示氣候變化CI,4種溫室氣體排放表示人類活動HI,經(jīng)主成分分析提取主要影響因素,通過GAMLSS模型,建立極端降水量與時間、CI和HI之間的關(guān)系。張東東等[11]研究大渡河流域降水極值的非一致性,采用該流域12個氣象站1960年至2013年最大日降水序列,經(jīng)MK趨勢檢驗,發(fā)現(xiàn)流域內(nèi)年最大日降水序列在上游5個站點呈不顯著減小趨勢,下游7個站點不顯著增加。經(jīng)GAMLSS模型計算,發(fā)現(xiàn)該流域日最大降水主要適用對數(shù)正態(tài)分布,且大部分站點適用平穩(wěn)模型。冶運濤等[12]采用MK趨勢分析法,通過降水發(fā)生率和降水貢獻率兩個指標,研究不同歷時連續(xù)降水(1961年~2005年)的時空演變特征,發(fā)現(xiàn)長江上游流域各歷時降水貢獻率以短歷時降水為主。短歷時(1 d和2 d)和長歷時(6 d和10 d)降水貢獻率分別在1976年和1999年發(fā)生突變。1~2 d降雨量占總降雨量的貢獻率呈不顯著增加趨勢,3~5d降雨貢獻率呈顯著增加趨勢,6~10 d降雨貢獻率呈減小趨勢。研究認為,短歷時降雨貢獻率增加易加劇洪澇災害形勢。

    本研究擬采用MK檢驗和GAMLSS模型,分析流域降水的空間分布特性及時間變化趨勢,探索流域極端降水的非平穩(wěn)性特征。

    2 研究方法

    2.1 建立研究指標

    為全面分析降雨極值的時空變化規(guī)律,本研究基于氣候變化常用指標[3,6]選取9個通用的極端降水指標進行定量分析(見表1)。

    表1 極端降水指標

    2.2 趨勢分析方法

    Mann-Kendall檢驗(以下簡稱“MK檢驗”)是一種非參數(shù)化檢驗方法[13-14],廣泛應用于水文氣象等時間序列的趨勢分析和顯著性檢驗[7,12]。MK檢驗計算樣本序列統(tǒng)計量Z,如果Z>0,趨勢增加,反之Z<0趨勢減小。選定顯著性水平α=0.05,當|Z|>1.96表明變化趨勢顯著。

    (1)計算統(tǒng)計量S

    (1)

    (2)

    式中,xk,xj分別為降雨時間序列中第k年和第j年的值,k>j;n為時間序列樣本總數(shù)。

    (2)計算統(tǒng)計量Z

    (3)

    當n≥10時,Var(S)=[n(n-1)(2n+5)]/18。

    2.3 非平穩(wěn)性分析模型

    針對非平穩(wěn)序列統(tǒng)計參數(shù)的時變特征,GAMLSS模型(Generalized Additive Models for Location, Scale and Shape)[15-16]可以反映降雨極值變量最優(yōu)分布函數(shù)及其均值(mu)、均方差(sigma)隨時間的變化模式,判別時間序列的平穩(wěn)性,并計算非平穩(wěn)序列的頻率設(shè)計值。模型分析通過以下步驟實現(xiàn)。

    (1)初選概率分布模型。根據(jù)水文氣象領(lǐng)域時間序列的分布特征,初選兩參數(shù)的Gamma、Lognormal、Gumbel和Weibull分布函數(shù)作為模型優(yōu)選的基礎(chǔ)。

    表2 概率分布函數(shù)

    (2)初選時間趨勢函數(shù)。針對均值(mu)和均方差(sigma)隨時間的變化模式,分別擬定3種趨勢函數(shù):①不隨時間變化的平穩(wěn)序列(S);②一次線性變化(L);③二次拋物線型變化(P)。

    (3)計算評價指標。根據(jù)赤池信息準則,在不同概率分布模型和各統(tǒng)計參數(shù)隨時間變化趨勢函數(shù)的所有組合中,選擇AIC值最小的模型為最優(yōu)模型。AIC=-2lnML+2k。其中,ML是參數(shù)估計的似然函數(shù)最大值;k為模型參數(shù)個數(shù)。

    (4)進行殘差判別。通過殘差散點圖、worm plot圖等方式,驗證模型的擬合效果。

    3 研究區(qū)域

    雅礱江流域地處青藏高原東南部,是長江上游最大的一級支流,呈南北條帶狀分布,跨越近8個緯度帶,流域面積約13萬km2,天然落差約3 870 m;地形復雜,谷嶺高差大[17],以高山峽谷為主。基于SRTM 90 m數(shù)字高程數(shù)據(jù),海拔高程分布于966~5 888 m,地勢呈現(xiàn)西北高東南低。流域中下游地區(qū)人口密度較大,工農(nóng)業(yè)生產(chǎn)發(fā)展,水資源豐沛,水能資源富集,是全國十四大水電基地之一。根據(jù)水電開發(fā)規(guī)劃,雅礱江中、下游按照“一庫十二級”規(guī)劃,水電站總裝機容量超過2 600萬kW。綜合考慮氣候環(huán)境變化影響,全面建立對流域降水等極端天氣的認知,對流域水電工程發(fā)揮興利除害作用具有重要意義。

    本研究選取流域內(nèi)18個國家基準和基本氣象臺站(見圖1)自建站至2016年12月日值降水數(shù)據(jù)集,計算全流域氣象臺站9種極端降水指標(見表1)的變化趨勢及平穩(wěn)性特征。

    圖1 雅礱江流域DEM、氣象站點分布及MK檢驗趨勢顯著變化站點和指標

    4 結(jié)果分析

    4.1 MK趨勢檢驗

    (1)無顯著變化趨勢的指標。全流域18個氣象臺站的年總降雨量PRCPTOT、3日(Rx3d)和5日連續(xù)最大降雨量(Rx5d)、極強降雨量R99p均無顯著變化趨勢。

    (2)全流域變化趨勢。全流域持續(xù)濕期CWD呈減小趨勢,其中清河站、木里站和冕寧站顯著減小。全流域極強降雨量R99p呈不顯著減小趨勢,但是強降雨量R95p在80%以上站點呈現(xiàn)增加趨勢,其中喜德站顯著增加。表明流域內(nèi)連續(xù)降水在減少,但是較強降水在增多。

    圖2 GAMLSS非平穩(wěn)性檢驗指標及空間分布

    (3)流域分區(qū)變化趨勢。流域上游(清水河至理塘)持續(xù)干期CDD呈減小趨勢,其中甘孜站顯著減小。流域中游和下游(雅江至米易)持續(xù)干期CDD呈增加趨勢,其中木里站和寧蒗站顯著增加。表明流域上游持續(xù)旱情減少,但是中下游反之。

    (4)單站呈現(xiàn)多項顯著趨勢。中下游的木里站1日最大降雨量Rx1d和持續(xù)干期CDD呈顯著增加趨勢,持續(xù)濕期CWD呈顯著減小趨勢;上游的甘孜站中雨日數(shù)R10呈顯著增加趨勢,持續(xù)干期CDD呈顯著減小趨勢;表明木里站的雨型分布更為集中,甘孜站降雨時程分配漸趨平緩。

    (5)總體趨勢。雖然目前氣候變化和人類活動對全流域降水特性無顯著影響;但是,不顯著的變化趨勢揭示,流域上游持續(xù)干期CDD和持續(xù)濕期CWD均減小、年降雨量PRCPTOT增加、5日連續(xù)最大降雨量Rx5d增加、極強降雨量R99p減小但強降雨量R95p增大;流域中下游持續(xù)干期CDD增加、持續(xù)濕期CWD減小、年降雨量PRCPTOT減小、5日最大降雨量Rx5d減小、極強降雨量R99p減小但強降雨量R95p略增大。反映了流域上游連續(xù)干旱和連續(xù)降雨均減少,降水增多且體現(xiàn)為5 d歷時和95%分位數(shù)以上的強降雨;流域中下游連續(xù)干旱增多、連續(xù)降雨減少但強降雨量并未減少,降雨的時程分配更為集中(見圖1)。

    4.2 GAMLSS非平穩(wěn)性分析

    采用GAMLSS模型,計算18個站9種降雨極值指標,最優(yōu)模型以對數(shù)正態(tài)Lognormal和Gamma分布為主。GAMLSS模型成果可反映降雨極值指標分布函數(shù)統(tǒng)計參數(shù)的變化趨勢:mu穩(wěn)態(tài),均值不隨時間變化;mu線性變化,均值隨時間單調(diào)增大或減?。籱u呈拋物線變化,均值隨時間變化存在拐點;sigma穩(wěn)態(tài),離散程度不隨時間變化;sigma線性變化,離散程度具有收斂或發(fā)散的趨勢;sigma拋物線變化,離散程度隨時間局部收斂或發(fā)散。

    (1)流域降雨極值平穩(wěn)性情況。1日(Rx1d)和3日連續(xù)最大降水量Rx3d、極強降水量R99p在流域上游和中游的11至12個站點為平穩(wěn)序列。中雨日數(shù)R10、持續(xù)干期CDD和持續(xù)濕期CWD僅在流域內(nèi)的5個站點呈現(xiàn)平穩(wěn)特征,且分布相對分散,不具有成片或帶狀分布規(guī)律(見圖2)。

    (2)單站平穩(wěn)性特征。西昌站僅2個指標(R10和CWD)為不平穩(wěn)序列,寧蒗和德昌站僅3個指標為不平穩(wěn)序列。但是,鹽源站僅R95p、鹽邊站僅CDD為平穩(wěn)序列,其他所有指標均不平穩(wěn)。觀測系列最平穩(wěn)和最不平穩(wěn)的站點均出現(xiàn)在流域下游地區(qū),該區(qū)域的平穩(wěn)性特征較為復雜。

    (3)站點對比分析。經(jīng)MK檢驗,木里站Rx1d和CDD顯著增加,CWD顯著減小。在GAMLSS非平穩(wěn)性分析中,木里站Rx1d、CDD、CWD、Rx5d、R95p、R99p、PRCPTOT等7個指標均呈現(xiàn)不平穩(wěn)特征,其中Rx1d、CDD、CWD等符合均值線性、均方差穩(wěn)態(tài)的變化特征(圖3a和圖3b)。喜德站R95p指標在MK檢驗中顯著增加,在GAMLSS模型中體現(xiàn)為均值線性增加,均方差拋物線變化,呈現(xiàn)先發(fā)散后收斂的變化規(guī)律(圖3c)。

    (4)非平穩(wěn)性對頻率設(shè)計值的影響。由于降雨極值的非平穩(wěn)性特征,不同時期的頻率設(shè)計值呈現(xiàn)差異。以喜德站R95p為例(圖3d),其100年一遇頻率值在1960年為598.87 mm,1995年上升至809.66 mm。此拐點出現(xiàn)后,2016年逐漸回落至739.87 mm;對于2年一遇常遇頻率值,1960年、1995年和2016年分別為496.43、545.22 mm和583.99 mm,呈現(xiàn)單調(diào)上升趨勢。同時可以發(fā)現(xiàn),喜德站R95p在1960年100年一遇頻率值與2016年2年一遇頻率值均為600 mm左右,量級相當。由此亦可見,曾經(jīng)的極端事件,在近年來頻繁出現(xiàn)。

    圖3 站點及指標非平穩(wěn)性特征

    5 研究結(jié)論

    本研究經(jīng)MK趨勢檢驗發(fā)現(xiàn),雅礱江全流域的總降雨量、3 d和5 d連續(xù)最大降雨量、極強降雨量均無顯著變化趨勢,部分站點呈現(xiàn)降水時程分布漸趨集中的趨勢。經(jīng)GAMLSS非平穩(wěn)性分析發(fā)現(xiàn),降水極值的概率分布函數(shù)以對數(shù)正態(tài)和GAMMA分布為主,1d和3d連續(xù)最大降水量、極強降水量時間序列的平穩(wěn)性較好;由于氣候和環(huán)境變化影響較小,流域上游和中游降水極值的平穩(wěn)性較好。通過GAMLSS非平穩(wěn)性研究也揭示了部分站點近年來2年一遇與40~50年前100年一遇設(shè)計值量級相當,反映了極端現(xiàn)象頻發(fā)的現(xiàn)象。

    猜你喜歡
    平穩(wěn)性降雨量極值
    極值點帶你去“漂移”
    基于非平穩(wěn)性度量的數(shù)字印章信息匹配
    降雨量與面積的關(guān)系
    極值點偏移攔路,三法可取
    基于遞歸量化分析的振動信號非平穩(wěn)性評價
    一類“極值點偏移”問題的解法與反思
    高重合度齒輪傳動的平穩(wěn)性分析及試驗
    洞庭湖區(qū)降雨特性分析
    信貸資源配置與我國經(jīng)濟發(fā)展平穩(wěn)性
    匹配數(shù)為1的極值2-均衡4-部4-圖的結(jié)構(gòu)
    又粗又爽又猛毛片免费看| 欧美日韩亚洲国产一区二区在线观看| x7x7x7水蜜桃| 色尼玛亚洲综合影院| 99久久99久久久精品蜜桃| 在线天堂最新版资源| 久久国产乱子伦精品免费另类| 国产av不卡久久| 18禁裸乳无遮挡免费网站照片| 亚洲成a人片在线一区二区| 国产av麻豆久久久久久久| 免费无遮挡裸体视频| 在线视频色国产色| 99国产综合亚洲精品| 亚洲人成电影免费在线| 日韩大尺度精品在线看网址| 亚洲av成人精品一区久久| 成人国产综合亚洲| 久久久成人免费电影| а√天堂www在线а√下载| x7x7x7水蜜桃| 欧美+日韩+精品| 一区二区三区高清视频在线| 亚洲精品色激情综合| 午夜福利高清视频| 国产精品永久免费网站| 日韩欧美精品免费久久 | 女警被强在线播放| 久久久久久大精品| 日韩欧美三级三区| 国产精品av视频在线免费观看| 精华霜和精华液先用哪个| 最新中文字幕久久久久| 免费电影在线观看免费观看| 亚洲av美国av| 久久久久久久午夜电影| 久久亚洲真实| 亚洲欧美日韩高清专用| xxxwww97欧美| 蜜桃久久精品国产亚洲av| 成人一区二区视频在线观看| 天天一区二区日本电影三级| 三级国产精品欧美在线观看| 欧美在线黄色| 中文字幕人成人乱码亚洲影| 国产精品98久久久久久宅男小说| 十八禁人妻一区二区| 亚洲无线在线观看| 国产一区在线观看成人免费| 色在线成人网| 日韩人妻高清精品专区| 亚洲乱码一区二区免费版| 亚洲片人在线观看| 中文字幕人妻丝袜一区二区| 国产一区二区在线av高清观看| 最近视频中文字幕2019在线8| 亚洲七黄色美女视频| 精品久久久久久久末码| 少妇人妻精品综合一区二区 | 不卡一级毛片| 国产一级毛片七仙女欲春2| 亚洲内射少妇av| 麻豆国产97在线/欧美| 90打野战视频偷拍视频| 久久精品91无色码中文字幕| 99精品久久久久人妻精品| 日韩高清综合在线| 国产成人av教育| 人妻久久中文字幕网| 久久久久精品国产欧美久久久| 看免费av毛片| 在线观看舔阴道视频| 一个人看的www免费观看视频| 欧美日韩黄片免| 婷婷精品国产亚洲av在线| avwww免费| 五月玫瑰六月丁香| 性色av乱码一区二区三区2| 18禁裸乳无遮挡免费网站照片| 亚洲欧美日韩高清专用| 欧美黄色淫秽网站| 757午夜福利合集在线观看| 男女下面进入的视频免费午夜| 欧美日韩乱码在线| 精品人妻一区二区三区麻豆 | 一级a爱片免费观看的视频| 狠狠狠狠99中文字幕| 国产亚洲欧美98| 搞女人的毛片| 国产精品,欧美在线| 日本一二三区视频观看| 脱女人内裤的视频| ponron亚洲| 琪琪午夜伦伦电影理论片6080| 在线观看日韩欧美| 免费看光身美女| 午夜福利欧美成人| 日本一二三区视频观看| 免费在线观看日本一区| 伊人久久大香线蕉亚洲五| 长腿黑丝高跟| 欧美3d第一页| 免费无遮挡裸体视频| 久久精品国产清高在天天线| 无人区码免费观看不卡| 岛国在线观看网站| 国产精品一及| 国产美女午夜福利| 一级毛片女人18水好多| 国产综合懂色| 久久婷婷人人爽人人干人人爱| 亚洲精品一卡2卡三卡4卡5卡| 老熟妇仑乱视频hdxx| 国产一区二区三区视频了| 免费在线观看影片大全网站| h日本视频在线播放| 国产精品爽爽va在线观看网站| www.999成人在线观看| 午夜亚洲福利在线播放| 桃红色精品国产亚洲av| 97超级碰碰碰精品色视频在线观看| www.熟女人妻精品国产| 午夜视频国产福利| 久久香蕉精品热| 亚洲欧美日韩无卡精品| 国产三级黄色录像| 亚洲欧美日韩卡通动漫| 亚洲成av人片免费观看| 成熟少妇高潮喷水视频| 久久久久国产精品人妻aⅴ院| 亚洲电影在线观看av| 变态另类丝袜制服| av片东京热男人的天堂| 亚洲va日本ⅴa欧美va伊人久久| 在线视频色国产色| 国产精品嫩草影院av在线观看 | 国产黄a三级三级三级人| 精品无人区乱码1区二区| 国产精品乱码一区二三区的特点| 亚洲片人在线观看| 一进一出抽搐动态| 看黄色毛片网站| 国产精品影院久久| 亚洲欧美日韩卡通动漫| 精品国内亚洲2022精品成人| 18禁黄网站禁片免费观看直播| 丝袜美腿在线中文| 日韩大尺度精品在线看网址| 男人舔女人下体高潮全视频| 在线观看日韩欧美| 日本熟妇午夜| 天堂√8在线中文| 免费在线观看亚洲国产| 又紧又爽又黄一区二区| 中文字幕人妻丝袜一区二区| 亚洲一区高清亚洲精品| 国产一区二区三区视频了| 久久久精品欧美日韩精品| 色av中文字幕| 日韩免费av在线播放| 天天添夜夜摸| 尤物成人国产欧美一区二区三区| 精品国产亚洲在线| 搡老妇女老女人老熟妇| 久久人妻av系列| 国内少妇人妻偷人精品xxx网站| 久久久久久大精品| 国产高清videossex| 亚洲一区二区三区色噜噜| 淫秽高清视频在线观看| 亚洲人成网站在线播| 欧美三级亚洲精品| 性色av乱码一区二区三区2| 日本三级黄在线观看| 1000部很黄的大片| 欧美在线一区亚洲| 国产色婷婷99| 国产男靠女视频免费网站| 欧美在线一区亚洲| av专区在线播放| 欧美av亚洲av综合av国产av| 亚洲精品久久国产高清桃花| 波多野结衣巨乳人妻| 国产精品久久久久久久电影 | 国产亚洲精品一区二区www| 欧美日韩国产亚洲二区| 日本免费一区二区三区高清不卡| 国产三级在线视频| 国产精品一区二区三区四区久久| 亚洲av第一区精品v没综合| 无人区码免费观看不卡| a级一级毛片免费在线观看| 亚洲 欧美 日韩 在线 免费| av欧美777| 中文字幕熟女人妻在线| 欧美乱妇无乱码| 婷婷丁香在线五月| 免费一级毛片在线播放高清视频| 夜夜看夜夜爽夜夜摸| 色综合亚洲欧美另类图片| 亚洲欧美日韩卡通动漫| 亚洲av不卡在线观看| 午夜免费男女啪啪视频观看 | 九色成人免费人妻av| 欧美3d第一页| 国产精品99久久99久久久不卡| 一进一出抽搐gif免费好疼| 欧美+日韩+精品| 精品99又大又爽又粗少妇毛片 | 激情在线观看视频在线高清| 久久久久久九九精品二区国产| 国产真人三级小视频在线观看| a级毛片a级免费在线| 欧美成狂野欧美在线观看| 欧美中文日本在线观看视频| 中文资源天堂在线| 午夜影院日韩av| а√天堂www在线а√下载| 亚洲精品456在线播放app | 久久婷婷人人爽人人干人人爱| 久久精品综合一区二区三区| 精品一区二区三区视频在线 | 一区福利在线观看| 性色av乱码一区二区三区2| 一本久久中文字幕| www日本在线高清视频| 国产av不卡久久| 亚洲人成网站在线播| 日韩免费av在线播放| 啪啪无遮挡十八禁网站| 91在线观看av| 一二三四社区在线视频社区8| 免费在线观看亚洲国产| 91久久精品电影网| 亚洲av熟女| 日日摸夜夜添夜夜添小说| 精华霜和精华液先用哪个| 精品一区二区三区视频在线 | 99riav亚洲国产免费| 操出白浆在线播放| 麻豆成人午夜福利视频| 国产一区二区在线观看日韩 | 成人三级黄色视频| 琪琪午夜伦伦电影理论片6080| 欧美日韩乱码在线| 美女高潮的动态| 国产在线精品亚洲第一网站| 国产精品99久久99久久久不卡| 无遮挡黄片免费观看| 国产亚洲欧美在线一区二区| 国产亚洲精品久久久com| 真实男女啪啪啪动态图| 国产真实伦视频高清在线观看 | 天天添夜夜摸| 国产色婷婷99| 日韩欧美三级三区| 国产av一区在线观看免费| av中文乱码字幕在线| 一本一本综合久久| 全区人妻精品视频| 中国美女看黄片| 亚洲国产欧美网| 国产欧美日韩一区二区三| 老司机午夜福利在线观看视频| 999久久久精品免费观看国产| 叶爱在线成人免费视频播放| 黑人欧美特级aaaaaa片| 久久久成人免费电影| xxxwww97欧美| 日韩精品青青久久久久久| a在线观看视频网站| 韩国av一区二区三区四区| 亚洲精品美女久久久久99蜜臀| 久久精品国产99精品国产亚洲性色| 日韩中文字幕欧美一区二区| 亚洲av熟女| 国产亚洲精品av在线| 免费高清视频大片| 精品一区二区三区av网在线观看| 在线观看一区二区三区| 久久精品国产99精品国产亚洲性色| 精品国产亚洲在线| 日本 欧美在线| 又黄又粗又硬又大视频| 欧美一区二区亚洲| 欧美高清成人免费视频www| 国产高清视频在线观看网站| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | av国产免费在线观看| 美女 人体艺术 gogo| 日韩中文字幕欧美一区二区| 啪啪无遮挡十八禁网站| 老司机深夜福利视频在线观看| 99精品欧美一区二区三区四区| 蜜桃亚洲精品一区二区三区| 亚洲欧美日韩无卡精品| 性欧美人与动物交配| 日韩欧美 国产精品| 18美女黄网站色大片免费观看| 在线播放无遮挡| 国产国拍精品亚洲av在线观看 | 国产真实伦视频高清在线观看 | 亚洲五月天丁香| 欧美绝顶高潮抽搐喷水| 午夜亚洲福利在线播放| 又黄又粗又硬又大视频| 亚洲av一区综合| a级一级毛片免费在线观看| 亚洲av美国av| 国产av一区在线观看免费| www日本黄色视频网| 老司机午夜十八禁免费视频| 在线观看av片永久免费下载| 国产美女午夜福利| 大型黄色视频在线免费观看| bbb黄色大片| 夜夜看夜夜爽夜夜摸| 欧美成人a在线观看| 国产精品乱码一区二三区的特点| 国产精品爽爽va在线观看网站| 真实男女啪啪啪动态图| 久久久久久久久久黄片| 精品一区二区三区av网在线观看| 国产私拍福利视频在线观看| 高清在线国产一区| 男人舔女人下体高潮全视频| 亚洲男人的天堂狠狠| 亚洲自拍偷在线| 久久久久性生活片| 日本与韩国留学比较| 国产成人av教育| 国产99白浆流出| 俺也久久电影网| 国语自产精品视频在线第100页| 制服丝袜大香蕉在线| 欧美日本视频| 在线观看美女被高潮喷水网站 | 少妇熟女aⅴ在线视频| 嫩草影院精品99| 国产伦精品一区二区三区四那| 91麻豆精品激情在线观看国产| 欧美一级毛片孕妇| 精品久久久久久,| 两人在一起打扑克的视频| 欧美一区二区精品小视频在线| 在线观看午夜福利视频| 一区二区三区激情视频| 日韩高清综合在线| 精品国产亚洲在线| 亚洲最大成人手机在线| 悠悠久久av| 俺也久久电影网| 久久国产精品影院| 国产精品久久视频播放| 久久香蕉精品热| 性色avwww在线观看| 在线a可以看的网站| 午夜福利在线观看吧| 亚洲欧美日韩东京热| 国产精品精品国产色婷婷| 免费av观看视频| 久久久久久大精品| 一本综合久久免费| 黄色视频,在线免费观看| a级一级毛片免费在线观看| 天堂动漫精品| 免费观看精品视频网站| 久久精品综合一区二区三区| 91在线观看av| 18禁美女被吸乳视频| 天天一区二区日本电影三级| 免费人成视频x8x8入口观看| 国产黄片美女视频| 欧美成人一区二区免费高清观看| 又紧又爽又黄一区二区| 成人无遮挡网站| 免费av不卡在线播放| 亚洲黑人精品在线| 国语自产精品视频在线第100页| 国产高清有码在线观看视频| 激情在线观看视频在线高清| 成人无遮挡网站| 女同久久另类99精品国产91| 午夜福利视频1000在线观看| 久久久久久人人人人人| 日本撒尿小便嘘嘘汇集6| 12—13女人毛片做爰片一| 亚洲精品影视一区二区三区av| 国产精品亚洲美女久久久| 极品教师在线免费播放| 老司机深夜福利视频在线观看| 免费大片18禁| www日本在线高清视频| 久久6这里有精品| 日韩精品中文字幕看吧| 日本 欧美在线| 小说图片视频综合网站| 久久久久九九精品影院| 亚洲av二区三区四区| 亚洲五月天丁香| 小说图片视频综合网站| 真人一进一出gif抽搐免费| 一进一出抽搐gif免费好疼| 精品国产超薄肉色丝袜足j| 99riav亚洲国产免费| 动漫黄色视频在线观看| 亚洲成av人片免费观看| 久久精品人妻少妇| 久久久国产成人精品二区| 欧美在线黄色| 精品久久久久久成人av| 精品国产美女av久久久久小说| 蜜桃久久精品国产亚洲av| 免费av毛片视频| 99热这里只有是精品50| 国产一区二区激情短视频| 天堂av国产一区二区熟女人妻| 亚洲av美国av| 精品国内亚洲2022精品成人| 天天添夜夜摸| 两性午夜刺激爽爽歪歪视频在线观看| 国产伦人伦偷精品视频| 国产精品国产高清国产av| 日本 欧美在线| 在线天堂最新版资源| 午夜免费激情av| 精品不卡国产一区二区三区| 91久久精品国产一区二区成人 | 神马国产精品三级电影在线观看| 在线观看美女被高潮喷水网站 | 免费看十八禁软件| 一区二区三区国产精品乱码| 国产精品久久久久久精品电影| 午夜激情福利司机影院| 身体一侧抽搐| 免费人成视频x8x8入口观看| 亚洲一区二区三区不卡视频| 日韩成人在线观看一区二区三区| 性色avwww在线观看| 国产亚洲精品一区二区www| 亚洲18禁久久av| 校园春色视频在线观看| 国产爱豆传媒在线观看| 麻豆国产97在线/欧美| 天堂网av新在线| 午夜福利欧美成人| 非洲黑人性xxxx精品又粗又长| 精品福利观看| 久久久久久久午夜电影| 最近最新中文字幕大全电影3| 五月玫瑰六月丁香| 国产真实伦视频高清在线观看 | 婷婷丁香在线五月| 欧美日韩黄片免| 亚洲七黄色美女视频| 亚洲av中文字字幕乱码综合| 国产91精品成人一区二区三区| 欧美3d第一页| 欧美+亚洲+日韩+国产| 午夜福利视频1000在线观看| 黄色日韩在线| 亚洲欧美日韩卡通动漫| 亚洲一区二区三区不卡视频| 亚洲无线在线观看| 成人欧美大片| 亚洲熟妇中文字幕五十中出| 久久国产精品人妻蜜桃| 男人舔女人下体高潮全视频| 99久国产av精品| 少妇人妻一区二区三区视频| 很黄的视频免费| 伊人久久大香线蕉亚洲五| 久久性视频一级片| 白带黄色成豆腐渣| 欧美成人a在线观看| 欧美日韩精品网址| 听说在线观看完整版免费高清| 国产在线精品亚洲第一网站| 特级一级黄色大片| 国产高清激情床上av| 成人鲁丝片一二三区免费| 一区二区三区高清视频在线| 成人永久免费在线观看视频| 操出白浆在线播放| 老司机在亚洲福利影院| 亚洲av电影在线进入| a级一级毛片免费在线观看| 亚洲18禁久久av| 真实男女啪啪啪动态图| 国产欧美日韩一区二区三| 国产精品久久电影中文字幕| 国产成人系列免费观看| 午夜激情欧美在线| 精品人妻偷拍中文字幕| 欧美日韩一级在线毛片| 精品国产亚洲在线| 一本久久中文字幕| 一区二区三区激情视频| 国产精品 国内视频| 黄色日韩在线| 国产精品久久视频播放| av天堂中文字幕网| 一边摸一边抽搐一进一小说| av天堂中文字幕网| 欧美xxxx黑人xx丫x性爽| 少妇高潮的动态图| 69人妻影院| 无限看片的www在线观看| 天天添夜夜摸| 久久久久久久久中文| 舔av片在线| 成人永久免费在线观看视频| 超碰av人人做人人爽久久 | 波野结衣二区三区在线 | 日本免费一区二区三区高清不卡| 亚洲av第一区精品v没综合| 激情在线观看视频在线高清| 国产精品久久视频播放| 日本黄色视频三级网站网址| 麻豆国产av国片精品| 身体一侧抽搐| 色精品久久人妻99蜜桃| 欧美+日韩+精品| 乱人视频在线观看| 色综合婷婷激情| 女生性感内裤真人,穿戴方法视频| 观看免费一级毛片| 色吧在线观看| 亚洲天堂国产精品一区在线| 色吧在线观看| 国产在线精品亚洲第一网站| 国产精品一区二区三区四区久久| 精品人妻一区二区三区麻豆 | 久久久久亚洲av毛片大全| 午夜激情福利司机影院| av中文乱码字幕在线| 国产单亲对白刺激| 国产一区二区在线观看日韩 | 亚洲av二区三区四区| 中国美女看黄片| 久久国产精品影院| 久久精品91无色码中文字幕| 亚洲人成网站在线播放欧美日韩| www.www免费av| 夜夜看夜夜爽夜夜摸| 久久久色成人| 91字幕亚洲| bbb黄色大片| 国产真实伦视频高清在线观看 | 性欧美人与动物交配| 久久久久久久久久黄片| 国产探花极品一区二区| 国产淫片久久久久久久久 | 亚洲av二区三区四区| 热99re8久久精品国产| av欧美777| 国产精品电影一区二区三区| 男女视频在线观看网站免费| 欧美又色又爽又黄视频| 亚洲天堂国产精品一区在线| 我要搜黄色片| 1000部很黄的大片| 蜜桃久久精品国产亚洲av| 久久久国产精品麻豆| 小说图片视频综合网站| 欧美丝袜亚洲另类 | 亚洲真实伦在线观看| 真实男女啪啪啪动态图| 18禁美女被吸乳视频| 91久久精品电影网| svipshipincom国产片| 久久香蕉国产精品| 波野结衣二区三区在线 | 女同久久另类99精品国产91| 亚洲一区高清亚洲精品| 天堂√8在线中文| 欧美日韩精品网址| 可以在线观看的亚洲视频| 国产免费av片在线观看野外av| 波多野结衣巨乳人妻| 日韩高清综合在线| 成人三级黄色视频| 丁香六月欧美| 天天躁日日操中文字幕| 欧美又色又爽又黄视频| 久久久久久久亚洲中文字幕 | 国产精品亚洲av一区麻豆| 久久精品国产清高在天天线| 老司机深夜福利视频在线观看| 成年人黄色毛片网站| 97人妻精品一区二区三区麻豆| 制服人妻中文乱码| а√天堂www在线а√下载| 日韩免费av在线播放| 午夜久久久久精精品| 精品免费久久久久久久清纯| 欧美一区二区精品小视频在线| 熟女少妇亚洲综合色aaa.| 亚洲精品色激情综合| 国产精品一区二区三区四区免费观看 | 亚洲成人久久爱视频| 一进一出好大好爽视频| 国产探花极品一区二区| 亚洲成人久久爱视频| 久久中文看片网| 又黄又粗又硬又大视频| 熟女人妻精品中文字幕| 高清日韩中文字幕在线| 亚洲欧美日韩东京热| 中文字幕精品亚洲无线码一区| 99久久精品热视频| 内射极品少妇av片p| 亚洲成人久久爱视频| 亚洲,欧美精品.| 亚洲av五月六月丁香网|