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

    渠江流域降水時空分布特征

    2022-04-08 07:40:32孟澤坤
    水土保持學(xué)報 2022年2期
    關(guān)鍵詞:渠江達(dá)川日數(shù)

    孟澤坤,王 彬

    (北京林業(yè)大學(xué)水土保持學(xué)院,重慶縉云山三峽庫區(qū)森林生態(tài)系統(tǒng)國家定位觀測研究站,北京 100083)

    近年來,在全球氣候變化背景下,降水時空分布更加不均勻,極端天氣事件明顯增多,旱澇災(zāi)害頻發(fā),針對降雨時空變異性的研究逐漸成為研究熱點(diǎn)。我國位于歐亞大陸東部,東臨西太平洋,海陸熱力差異明顯,季風(fēng)氣候是主要的氣候特征之一,降水等氣象要素存在明顯的時空變異性。渠江流域是長江上游水系含沙量最大的河流——嘉陵江左岸的最大支流,流域降水量受季風(fēng)的影響表現(xiàn)出明顯的季節(jié)性差異,同時受地形、地貌、土壤、地質(zhì)等條件的影響,該流域氣象水文要素空間特征差異大,水土流失嚴(yán)重,徑流含沙量高。近年來,渠江流域洪水災(zāi)害頻發(fā)。2011年“9·18”暴雨洪災(zāi)中,渠江流域發(fā)生了歷史上最大洪水,經(jīng)濟(jì)損失高達(dá)141億元。2021年7月,強(qiáng)降雨導(dǎo)致渠江流域部分中小河流發(fā)生超警戒超保證洪水,渠江干流三匯站洪峰流量接近歷史峰值,人民生命財產(chǎn)安全受到嚴(yán)重威脅。此外,長江流域中三峽水系、洞庭湖水系、鄱陽湖水系等34條河流形勢也非常嚴(yán)峻。因此,針對渠江流域降水及暴雨事件的時空分布特征的研究是極有必要的。

    現(xiàn)有研究表明,嘉陵江流域降水量整體呈下降趨勢,但支流渠江流域變化趨勢與其相反。嘉陵江流域年降水量存在35年和21年的周期性變化,各季節(jié)降雨量存在很大的差異,夏季和冬季降水量增加,春季降水量呈微弱減少趨勢,秋季降水量顯著減少。渠江流域位于嘉陵江流域東南部,年均降水量較大,夏季降水量增加趨勢最顯著。此外,杜華明等對嘉陵江流域的研究表明,年平均暴雨日數(shù)與年平均降水量存在顯著正相關(guān)關(guān)系,相關(guān)系數(shù)為0.84(<0.01),年平均降水量的多少是影響年平均暴雨出現(xiàn)日數(shù)的重要因素;羅玉等研究發(fā)現(xiàn),2005年開始四川盆地的極端日降水事件明顯增多;符艷紅等研究指出,嘉陵江流域東南部——渠江流域降水日數(shù)最多;袁夢等對渠江流域1970—2012年降水資料分析研究得出,渠江流域汛期降水量增大,暴雨頻率增加,降水趨于極端的結(jié)論。

    現(xiàn)有的研究多是對嘉陵江流域進(jìn)行整體的研究,針對渠江流域的深入研究較少,且渠江流域位于嘉陵江流域東部,受季風(fēng)影響更為明顯,其降水量在時間和空間上的表現(xiàn)與嘉陵江整體趨勢存在差異。因此,本文從年、季節(jié)2個尺度對渠江流域的降水量進(jìn)行時空分布特征研究,以期揭示渠江流域降水時空分異規(guī)律,并對渠江流域汛期暴雨事件的時空分布特征進(jìn)行了研究,為該地區(qū)乃至長江流域季風(fēng)氣候區(qū)防范全球氣候變化影響下的洪澇災(zāi)害提供理論依據(jù)。

    1 流域概況

    渠江是長江支流嘉陵江左岸最大的支流,流經(jīng)四川、陜西、重慶。流域地跨東經(jīng)106°28′—109°00′,北緯30°00′—32°48′,面積3.7萬km。渠江流域地勢北部高南部低,海拔178~2 669 m。土壤類型以棕壤、黃棕壤、紫色土為主。渠江流域?qū)儆趤啛釒駶櫦撅L(fēng)氣候,氣候溫和、雨量充沛,平均降水量變化于785.65~1 737.01 mm。受季風(fēng)性氣候的影響,流域降水量年內(nèi)分布不均,多集中在5—9月,約占全年降水量的63%。

    2 材料與方法

    2.1 數(shù)據(jù)資料

    本研究收集整理了渠江流域內(nèi)部及外緣共13個基本站點(diǎn)的氣象數(shù)據(jù),數(shù)據(jù)來源于中國氣象數(shù)據(jù)網(wǎng)(http://data.cma.cn/)。其中,流域內(nèi)部氣象站點(diǎn)3個,外緣補(bǔ)充站點(diǎn)10個(表1)。達(dá)川站、巴中站、萬源站、高坪站、閬中站、廣元站6個站點(diǎn)數(shù)據(jù)資料為1961—2020年逐日降雨量,萬州站、沙坪壩站、安康站、石泉站、漢中站、略陽站、佛坪站7個站點(diǎn)數(shù)據(jù)資料為1961—2020年逐月降雨量。

    表1 氣象站基本情況

    本文根據(jù)氣象統(tǒng)計(jì)法劃分各季節(jié),春、夏、秋、冬4個季節(jié)分別為3—5月、6—8月、9—11月、12月至翌年2月。文中的汛期均指各年5—9月。本文將一日內(nèi)24 h累計(jì)降水量≥50 mm記為1個暴雨日。

    2.2 分析方法

    2.2.1 降水傾向率 流域降水的傾向率采用一次線性方程式表示:

    =+(=1,2,…,)

    (1)

    式中:為降水量(mm);為時間(a);為常數(shù)項(xiàng);為線性趨勢的斜率。

    2.2.2 Kriging空間插值法 采用較其他傳統(tǒng)方法更為可靠和精確的Kriging空間插值法,依據(jù)協(xié)方差函數(shù)進(jìn)行空間建模和插值,基于區(qū)域樣本點(diǎn)的實(shí)測數(shù)據(jù),根據(jù)待估點(diǎn)鄰域樣本點(diǎn)的空間位置關(guān)系,對待估點(diǎn)進(jìn)行無偏最優(yōu)估計(jì)。本研究采用普通克里金法,選用球面模型作為變異函數(shù)理論模型,對降雨數(shù)據(jù)進(jìn)行插值計(jì)算。

    (2)

    式中:()為點(diǎn)處雨量估計(jì)值(mm);()為站點(diǎn)的實(shí)測雨量(mm);為站點(diǎn)的權(quán)重;為雨量站數(shù)量。

    2.2.3 累積距平法 距平是指某時段的要素值與同一時段該要素平均值的差值,距平值可以直觀地反映出序列的變化情況。當(dāng)累積距平曲線呈上升趨勢時,表示序列距平值增加,序列呈增加趨勢;當(dāng)累積距平曲線呈下降趨勢時,表示序列距平值減小,序列呈下降趨勢。對于時間序列,累積距平可表示為:

    (3)

    2.2.4 Mann-Kendall趨勢檢驗(yàn)法 Mann-Kendall方法是一種非參數(shù)統(tǒng)計(jì)檢驗(yàn)方法,可用于不符合正態(tài)分布特征的時間序列數(shù)據(jù)的長期趨勢分析。該方法計(jì)算過程簡單,檢測范圍寬,干擾度小,應(yīng)用范圍廣,具體計(jì)算方法見參考文獻(xiàn)[19]。

    2.2.4 Wavelet分析 小波(Wavelet)分析是在傅里葉(Fourier)分析基礎(chǔ)上發(fā)展而來的。本文使用Morlet小波函數(shù)對序列的變化趨勢進(jìn)行分析。

    (4)

    式中:為常數(shù);e為虛數(shù);為時間(a)。

    對于時間序列,其連續(xù)小波變換為:

    (5)

    為了判斷序列的主要周期,進(jìn)行小波方差計(jì)算,即將小波系數(shù)的平方值在域積分,其計(jì)算公式為:

    (6)

    小波方差隨的變化過程即為小波方差圖,對應(yīng)峰值處的尺度稱為該序列的主要時間尺度,曲線最高點(diǎn)所對應(yīng)的時間尺度為第1主周期。

    3 結(jié)果與分析

    3.1 年尺度降雨量時空分布特征

    3.1.1 年尺度降雨時間變化特征 渠江流域近60年平均降水量為1 126.17 mm,年均降水量為785.65~1 737.01 mm,年降水量最大值與最小值分別出現(xiàn)在1983年和2001年。流域年降水量整體呈增加趨勢,趨勢方程為=0.69-251.38,年降水量以6.9 mm/10 a的速率增加,60年增長量為41.40 mm,增加趨勢不顯著?;瑒悠骄幚砗?,年降水量增加趨勢減弱,傾向率為3.3 mm/10 a,60年增長量為19.8 mm(圖2)。1980—1985年渠江流域年降水量呈大幅增加趨勢,1992—2003年,年降水量呈大幅減少趨勢,其他年代波動趨勢較緩,這與段文明等對嘉陵江1960—2009年降水量的研究結(jié)論一致。近60年渠江流域年降水量波動趨勢明顯,年降水量距平值呈正負(fù)值交替變化,降水偏多年28年,偏少年32年,整體呈上升趨勢。

    圖1 研究區(qū)域及氣象站點(diǎn)分布位置

    圖2 1961-2020年渠江年降水量變化情況

    對渠江流域年降水量進(jìn)行M—K突變檢驗(yàn),UF、UB均在U水平線內(nèi),年際變化趨勢不明顯(圖3)。1963—1977年、1981—1992年、2010—2020年,UF>0,降水量在這3個階段呈上升趨勢。1978—1980年、1993—2009年,UF<0,降水量呈下降趨勢。在0.05的置信區(qū)間內(nèi),UF、UB 2條曲線交點(diǎn)較多,分別為1963年、1970年、1971年、1976年、1980年、1983年、2003年、2015年和2018年,表明渠江流域年降水存在多個突變點(diǎn),波動劇烈。1984—2004年,UF和UB曲線無交點(diǎn),渠江流域年降水量變化趨勢比較穩(wěn)定,無較大波動,這與杜華明等在嘉陵江流域所得結(jié)論相似。

    圖3 1961-2020渠江流域降水量M-K突變檢驗(yàn)

    1961—2020年渠江流域年降水存在多個振蕩中心,主振蕩周期分別是27年、15年和7年。小波方差在27年左右達(dá)到最大值,27年是其第1主周期。在27年尺度上,年降水量呈現(xiàn)“偏多—偏少—偏多—偏少—偏多—偏少—偏多”波動變化。15年尺度上,年降水量呈現(xiàn)“偏多—偏少—偏多—偏少—偏多—偏少—偏多—偏少—偏多—偏少—偏多—偏少”12次交替(圖4)。

    圖4 降水量Morlet小波變換實(shí)部時頻分布與小波方差

    3.1.2 年尺度降雨空間變化特征 渠江流域多年平均降水量存在明顯的空間變異性,呈現(xiàn)由西部向東北部遞增的趨勢(圖5),最高降水量位于東北部萬源地區(qū),多年平均降水量為1 230~1 250 mm。降水量最小值位于流域西部巴中地區(qū),多年平均降水量為1 130~1 150 mm,低于最高多年平均降水量8.00%~8.13%。渠江流域地處四川盆地東緣,地勢東北高西南低,受地勢抬升和山脈的阻擋,東南季風(fēng)裹挾的大量水汽集中在流域東南部,形成地形雨,致使渠江流域東部、南部降水較為豐沛,西北部降水量較少。

    圖5 60年平均降水量與年降水量傾向率分布

    渠江流域年降水量變化傾向率也存在明顯的空間差異。多年平均降水傾向率呈由西北部向東南部遞增的趨勢,傾向率為正值的區(qū)域居多,即渠江流域整體年降水呈增加趨勢,東部及南部增加趨勢較快,年降水傾向率為10~15 mm/10 a,60年降水量增長60~90 mm。西北部降水傾向率<-5 mm/10 a,60年降水量減少30~60 mm。在全球氣候變暖背景下,極端強(qiáng)降水事件頻發(fā),渠江流域暴雨頻率增加可能是引起流域東部降水量增加的原因之一。

    3.2 季節(jié)尺度降水時空分布特征

    渠江流域處于亞熱帶季風(fēng)氣候區(qū),降水季節(jié)間差異較大。在季節(jié)尺度上量化降水的時空分布特征,有助于細(xì)化分析渠江流域降水時空分異規(guī)律,區(qū)分流域內(nèi)部各季節(jié)降水量分布重點(diǎn)區(qū)域。

    3.2.1 季節(jié)尺度降水量的時間變化特征 1961—2020年渠江流域各季節(jié)多年平均降水量分別為春季257.23 mm,夏季530.46 mm,秋季321.38 mm,冬季38.59 mm,各季節(jié)降水量呈波動變化。春季降水量隨年份增加呈下降趨勢,趨勢不顯著,傾向率為-4.1 mm/10 a,60年降水量減少24.6 mm;夏季降水量呈上升趨勢,傾向率為16.9 mm/10 a,60年增加101.4 mm;秋季降水量呈下降趨勢,傾向率為-7 mm/10 a,60年減少42 mm;冬季降水量呈增加趨勢,傾向率為0.78 mm/10 a,60年減少4.68 mm(圖6)。5年滑動平均處理后,春、夏、冬季降水量變化趨勢均有所減緩(圖7)。其中,夏季降水量上升趨勢顯著,傾向率為14.1 mm/10 a(<0.05),較處理前減少2.8 mm/10 a;秋季降水量下降趨勢增強(qiáng),傾向率為-7.7 mm/10 a(<0.05)。

    圖6 各季節(jié)降水量年際變化

    圖7 5年滑動平均處理后季節(jié)降水量變化趨勢

    渠江流域春季降水量在1961—1979年間呈上升趨勢,1980后主要呈下降趨勢(圖8),UF、UB均在U水平線內(nèi),年際變化趨勢不明顯。1975—2005年間,無突變點(diǎn),突變點(diǎn)集中在1961—1975年和近15年間。夏季降水量除1970年、1972年和1973年外,皆呈上升趨勢,1982—1997年、2010—2015年上升趨勢尤為顯著(<0.05),1961—1970年突變點(diǎn)較多,但近15年僅有2016年發(fā)生了突變。渠江流域秋季降水量在1977—2010年區(qū)間內(nèi)呈下降趨勢,1997—2003年,下降趨勢最為顯著(<0.05),突變發(fā)生情況同春季相似。冬季降水量在上世紀(jì)90年代前主要呈下降趨勢,90年代后主要呈上升趨勢。上世紀(jì)90年代由于西太平洋副熱帶高壓的增強(qiáng),渠江流域的年降水量整體呈減少趨勢,但由于流域處于亞熱帶季風(fēng)氣候區(qū),季節(jié)間差異較大,各季節(jié)降水量趨勢各異。出現(xiàn)春季和秋季降水量減少,夏季和冬季降水增加的現(xiàn)象。

    圖8 1961-2020年各季節(jié)降水量M-K檢驗(yàn)

    3.2.2 不同季節(jié)降水量的空間分布特征 渠江流域各季節(jié)多年平均降水量差別較大,存在明顯的空間差異。春季多年平均降水量呈現(xiàn)由西北部向東南部遞增的趨勢(圖9),流域東南部達(dá)川地區(qū)降水量為305~320 mm,西部巴中地區(qū)春季多年平均降水量為230~260 mm。流域東部受季風(fēng)影響大,夏季、秋季降水量相對較高,兩季多年平均降水量的空間分布情況相似,都呈由流域北部向流域南部遞減的趨勢。夏季、秋季多年平均降雨量最低值均出現(xiàn)在流域南部達(dá)川地區(qū),夏季降水量為515~530 mm,秋季降水量為310~320 mm。流域東北部萬源地區(qū)夏季和秋季降水量最多,夏季降水量在600 mm左右,秋季降水量在350 mm左右。冬季多年平均降水量呈由東北部向南部遞增的趨勢,平均降水量為38.35 mm,南部達(dá)川地區(qū)雨量較多,高于平均水平14.73%~25.16%。春季和冬季,渠江流域受季風(fēng)的影響較小,雨量高值區(qū)由海拔較高的流域東北部向海拔較低的流域南部移動,形成春、冬季與夏、秋季降水量相異的分布格局。

    圖9 1961-2020年各季節(jié)多年平均降水量分布

    渠江流域春季降水傾向率整體呈下降趨勢,流域東部下降趨勢最快(圖10),降水傾向率在-7.5~-5.5 mm/10 a,流域東部60年降水量減少33~45 mm。夏季降水量呈上升趨勢,傾向率的空間分布情況與降水量的分布情況相似,流域北部上升趨勢最顯著,降水傾向率在18~23 mm/10 a,60年降水量增加108~138 mm。秋季降水整體呈下降趨勢,降水傾向率以巴中地區(qū)為負(fù)值中心和達(dá)川地區(qū)為正值中心輻射分布。流域西部巴中地區(qū)降水傾向率-17.5~-14 mm/10 a,60年降水量減少84~105 mm,東南部達(dá)川地區(qū)秋季降水量以0~3.5 mm/10 a的趨勢緩慢增加。冬季降水量整體上升趨勢在0.3~1.5 mm/10 a,增長趨勢緩慢,流域東南部達(dá)川地區(qū)上升趨勢最緩。

    圖10 60年季節(jié)降水傾向率分布

    3.3 汛期暴雨事件時空分布特征

    3.3.1 暴雨天數(shù)年際變化特征 60年間,渠江流域汛期暴雨日數(shù)平均為4.28日,暴雨日數(shù)介于1.05~8.49天。暴雨日數(shù)整體呈增加趨勢,暴雨日數(shù)以0.21 d/10 a的速率增加,增加趨勢不顯著。滑動平均處理后,暴雨日數(shù)的增加趨勢稍有減弱,傾向率為0.17 d/10 a(圖11)。

    圖11 1961-2020年汛期暴雨日數(shù)年際變化與暴雨日數(shù)距平變化

    1961—2020年渠江流域汛期暴雨日數(shù)波動明顯,其中1961年暴雨日數(shù)最少為1.05天,1983年暴雨日數(shù)達(dá)到60年最大值8.49天。1961—2010年,暴雨日數(shù)累積距平值為負(fù)值,2002年后累積距平值上升趨勢,在-2~3浮動。其中,1979—1985年、2003—2015年暴雨日數(shù)明顯偏多,這與袁夢等所得2003—2012年渠江流域暴雨發(fā)生頻率及汛期最大雨強(qiáng)增大的結(jié)論相匹配。

    渠江流域汛期暴雨日數(shù)年際變化呈不顯著的緩慢上升趨勢(圖12)。汛期暴雨日數(shù)增加的趨勢與渠江流域夏季降水量增加相匹配,同時與嘉陵江流域的變化趨勢一致。1965年、1983年和1984年暴雨日數(shù)有增加趨勢,2014年有減少趨勢(<0.05)。2002—2020年,暴雨日數(shù)有持續(xù)上升趨勢。在U的置信區(qū)間內(nèi),UF、UB 2條曲線交點(diǎn)較多,分別為1962年、1970年、1971年、1972年、1977年、1988年和1996年,渠江暴雨日數(shù)存在多個突變點(diǎn),波動劇烈。

    圖12 1961-2020年汛期暴雨日數(shù)M-K檢驗(yàn)

    1961—2020年渠江流域汛期暴雨日數(shù)有明顯的周期性變化并存在多個振蕩中心(圖13),主振蕩周期分別是14年、21年和7年,振蕩周期與年降水量變化相似。小波方差在14年左右達(dá)到最大值,14年是其第1主周期。在14年尺度上,暴雨日數(shù)呈現(xiàn)“偏少—偏多—偏少—偏多—偏少—偏多—偏少—偏多—偏少—偏多—偏少—偏多—偏少—偏多”14次波動變化。21年尺度上,暴雨日數(shù)呈現(xiàn)“偏多—偏少—偏多—偏少—偏多—偏少—偏多—偏少—偏多”9次交替變化。

    圖13 1961-2020年汛期暴雨日數(shù)小波實(shí)部與小波方差

    3.3.2 暴雨日數(shù)空間分布特征 渠江流域多年平均暴雨日數(shù)呈由西南部向東北部遞增的趨勢(圖14),該趨勢與夏季降水量分布趨勢相似。暴雨日數(shù)最大值位在東北部萬源地區(qū),多年平均暴雨日數(shù)在4.7~4.9天。暴雨日數(shù)最小值位于流域南部達(dá)川地區(qū),多年平均暴雨日數(shù)在3.7~3.9天。

    圖14 渠江流域多年平均汛期暴雨日數(shù)分布與流域各氣象站點(diǎn)汛期暴雨日數(shù)變化趨勢

    渠江流域內(nèi)部站點(diǎn)汛期暴雨日數(shù)均呈上升趨勢,外緣站點(diǎn)廣元站呈緩慢下降趨勢,趨勢不顯著。呈上升趨勢的地區(qū)中,流域南部達(dá)川地區(qū)的上升趨勢最快為0.36 d/10 a,高于流域整體上升趨勢0.15 d/10 a,其他地區(qū)上升趨勢由快到慢分別是流域西部巴中地區(qū)、流域東北部萬源地區(qū)、流域外西南部高坪地區(qū)、流域外西部閬中地區(qū)。渠江流域暴雨發(fā)生頻率整體呈上升趨勢,南部、西部低海拔地區(qū)的趨勢較為明顯。雖東北部上升趨勢較緩,但是暴雨的主要集中區(qū)域,洪澇災(zāi)害的防范工作仍需加強(qiáng)。

    4 結(jié) 論

    (1)渠江流域60年年均降水量為1 126.17 mm,年降水量整體呈緩慢上升趨勢,降水傾向率為6.9 mm/10 a,流域年降水量波動劇烈,存在27年、15年的周期性變化。渠江流域多年平均降水量呈由西部向東北部遞增的趨勢,降水量高值區(qū)位于流域東北部萬源地區(qū),低值區(qū)位于西部巴中地區(qū)。年降水量傾向率呈由西北部向東南部遞增的趨勢,東部及南部降水傾向率為10~15 mm/10 a。

    (2)渠江流域降水呈季節(jié)性變化,春季、秋季降水量呈下降趨勢,夏季、冬季降水量呈上升趨勢。5年滑動平均表明,夏季傾向率為16.9 mm/10 a,秋季降水傾向率為-7.7 mm/10 a。4個季節(jié)降水量波動明顯,突變點(diǎn)較多,春、夏、秋季突變點(diǎn)主要集中于1961—1980年和2010—2020年間,冬季突變點(diǎn)集中于1985—2005年間。渠江流域春、夏、秋、冬四季降水傾向率最高值分別位于流域南部、東北部、東南部、西部。

    (3)渠江流域年均汛期暴雨日數(shù)為4.28天,暴雨日數(shù)整體呈增加趨勢,增加速率為0.21 d/10 a。渠江流域汛期暴雨日數(shù)波動明顯,1965年、1983年、1984年有顯著增加趨勢。汛期暴雨日數(shù)有明顯的周期性變化,并存在14年、21年多個振蕩中心。渠江流域多年平均暴雨日數(shù)呈由西南部向東北部遞增的趨勢,流域南部達(dá)川地區(qū)的上升趨勢最快,高于流域整體上升趨勢0.15 d/10 a。

    猜你喜歡
    渠江達(dá)川日數(shù)
    漢江上游漢中區(qū)域不同等級降水日數(shù)的氣候變化特征分析
    綠色科技(2022年16期)2022-09-15 03:04:46
    達(dá)州市達(dá)川區(qū):開展和諧勞動關(guān)系“春風(fēng)行動”
    天津市濱海新區(qū)塘沽地域雷暴日數(shù)變化規(guī)律及特征分析
    天津科技(2020年2期)2020-03-03 05:09:48
    渠江邊
    渠江邊
    脫貧攻堅(jiān)的渠江畫廊
    基于RMP分析的渠江源茶文化旅游開發(fā)研究
    達(dá)川區(qū)退耕還林典型經(jīng)驗(yàn)調(diào)查報告
    對如何提高達(dá)川區(qū)造林質(zhì)量的思考與建議
    海南省雷暴日數(shù)年代際變化特征
    18禁裸乳无遮挡动漫免费视频 | 欧美成人a在线观看| 亚洲美女视频黄频| 欧美日韩综合久久久久久| 国产视频内射| 免费黄色在线免费观看| 在线 av 中文字幕| 黄色视频在线播放观看不卡| 99久国产av精品国产电影| 亚洲欧美日韩卡通动漫| 国产视频首页在线观看| 成人漫画全彩无遮挡| 天美传媒精品一区二区| 国产精品爽爽va在线观看网站| 精品久久久久久久久av| 一级片'在线观看视频| 爱豆传媒免费全集在线观看| 国产精品成人在线| 国产极品天堂在线| 久久99蜜桃精品久久| 亚洲精品,欧美精品| 永久免费av网站大全| 精品人妻偷拍中文字幕| 国产乱人视频| 人人妻人人看人人澡| av黄色大香蕉| 国产美女午夜福利| 国产美女午夜福利| 国产片特级美女逼逼视频| av在线天堂中文字幕| 精品一区在线观看国产| 九色成人免费人妻av| 亚州av有码| 99久久中文字幕三级久久日本| 18禁裸乳无遮挡动漫免费视频 | 建设人人有责人人尽责人人享有的 | 69人妻影院| 一级毛片电影观看| 色哟哟·www| 99久久中文字幕三级久久日本| 免费电影在线观看免费观看| xxx大片免费视频| 亚洲内射少妇av| 国产一区亚洲一区在线观看| 国产男女内射视频| 亚洲精品乱久久久久久| 亚洲精品日韩在线中文字幕| 女人十人毛片免费观看3o分钟| 国产精品无大码| 波野结衣二区三区在线| 在线观看av片永久免费下载| 97超碰精品成人国产| 国产午夜精品一二区理论片| 九九爱精品视频在线观看| 18禁动态无遮挡网站| 直男gayav资源| 国产色婷婷99| 成年女人在线观看亚洲视频 | 午夜福利高清视频| 欧美亚洲 丝袜 人妻 在线| 中文欧美无线码| 美女主播在线视频| 一级毛片 在线播放| 精品久久久久久久久av| 综合色av麻豆| 美女xxoo啪啪120秒动态图| 校园人妻丝袜中文字幕| 69av精品久久久久久| 禁无遮挡网站| 亚洲人与动物交配视频| 国产视频内射| 男人狂女人下面高潮的视频| 久久精品国产鲁丝片午夜精品| 国产爽快片一区二区三区| 九九爱精品视频在线观看| 日韩视频在线欧美| 国产成人午夜福利电影在线观看| 亚洲精华国产精华液的使用体验| 国产成人一区二区在线| 天天一区二区日本电影三级| 色综合色国产| 老司机影院毛片| 交换朋友夫妻互换小说| 只有这里有精品99| 中文欧美无线码| 国产精品av视频在线免费观看| 97超视频在线观看视频| 欧美少妇被猛烈插入视频| 观看美女的网站| 久久久久九九精品影院| 午夜爱爱视频在线播放| 成人二区视频| 国产 一区 欧美 日韩| 成年女人在线观看亚洲视频 | av免费观看日本| 少妇丰满av| 欧美日韩精品成人综合77777| 亚洲国产精品专区欧美| 久久久国产一区二区| 免费大片18禁| 女人十人毛片免费观看3o分钟| 一区二区三区乱码不卡18| 国产片特级美女逼逼视频| 啦啦啦啦在线视频资源| av在线app专区| 免费少妇av软件| 精品久久久久久久人妻蜜臀av| 国产精品一区二区三区四区免费观看| 大陆偷拍与自拍| 亚洲欧美一区二区三区黑人 | 国产真实伦视频高清在线观看| 搡女人真爽免费视频火全软件| 免费不卡的大黄色大毛片视频在线观看| 高清日韩中文字幕在线| 亚州av有码| 亚洲欧洲日产国产| 听说在线观看完整版免费高清| 久久97久久精品| 欧美激情在线99| 毛片女人毛片| 少妇 在线观看| 国产精品久久久久久久久免| 真实男女啪啪啪动态图| 99久久精品一区二区三区| 亚洲国产精品999| 色5月婷婷丁香| 欧美最新免费一区二区三区| 国产熟女欧美一区二区| av播播在线观看一区| 边亲边吃奶的免费视频| 国产精品福利在线免费观看| 菩萨蛮人人尽说江南好唐韦庄| 深爱激情五月婷婷| av卡一久久| 成年免费大片在线观看| 国产大屁股一区二区在线视频| 中文字幕久久专区| 日日摸夜夜添夜夜添av毛片| 亚洲av二区三区四区| 男女边摸边吃奶| 九九久久精品国产亚洲av麻豆| 久久久国产一区二区| 91在线精品国自产拍蜜月| 亚洲欧美清纯卡通| 久久久久久久大尺度免费视频| 国产日韩欧美在线精品| 亚洲丝袜综合中文字幕| 69av精品久久久久久| 久久97久久精品| 婷婷色综合www| 久久精品久久久久久久性| 日本三级黄在线观看| 亚洲人成网站高清观看| 成人一区二区视频在线观看| 大香蕉久久网| 日本熟妇午夜| 22中文网久久字幕| 国产精品秋霞免费鲁丝片| 日韩av免费高清视频| 久久久久久久精品精品| 久久久色成人| 晚上一个人看的免费电影| 国产成人精品一,二区| 国产一区二区三区综合在线观看 | 成年人午夜在线观看视频| 久久久久性生活片| 日本三级黄在线观看| 久热这里只有精品99| 岛国毛片在线播放| 中文字幕av成人在线电影| 久久久久精品久久久久真实原创| 亚洲精品国产成人久久av| 91久久精品国产一区二区成人| 午夜激情福利司机影院| av线在线观看网站| 赤兔流量卡办理| 欧美日韩一区二区视频在线观看视频在线 | 国产老妇女一区| av免费观看日本| 有码 亚洲区| 午夜福利视频精品| 国产精品一二三区在线看| 国产色爽女视频免费观看| 国产黄片美女视频| 男人舔奶头视频| 亚州av有码| 免费观看无遮挡的男女| 乱系列少妇在线播放| 国产亚洲av片在线观看秒播厂| 99re6热这里在线精品视频| 91久久精品电影网| 免费av毛片视频| 免费大片18禁| 熟女人妻精品中文字幕| 特大巨黑吊av在线直播| 秋霞伦理黄片| 肉色欧美久久久久久久蜜桃 | 99热6这里只有精品| 人妻系列 视频| 成人综合一区亚洲| 不卡视频在线观看欧美| 成年女人看的毛片在线观看| 日韩欧美精品v在线| 亚洲精品国产成人久久av| 国产免费一区二区三区四区乱码| 特大巨黑吊av在线直播| 免费看不卡的av| 日本免费在线观看一区| 亚洲精品日韩av片在线观看| 亚洲精品国产成人久久av| av女优亚洲男人天堂| 亚洲美女搞黄在线观看| 久久久久久久久久成人| 成人毛片a级毛片在线播放| 久久久久久久久大av| 天美传媒精品一区二区| 免费在线观看成人毛片| av在线老鸭窝| 免费av不卡在线播放| 亚洲天堂av无毛| 亚洲精品国产色婷婷电影| 男人和女人高潮做爰伦理| 精品人妻一区二区三区麻豆| 十八禁网站网址无遮挡 | 欧美日韩国产mv在线观看视频 | 性色av一级| 赤兔流量卡办理| 亚洲av欧美aⅴ国产| 国产成人免费观看mmmm| 久久97久久精品| 国产精品久久久久久久久免| 免费看光身美女| 亚洲真实伦在线观看| 欧美另类一区| 精品一区二区三区视频在线| 免费黄频网站在线观看国产| .国产精品久久| 自拍偷自拍亚洲精品老妇| 一级毛片电影观看| 国产 精品1| 黄色欧美视频在线观看| 欧美三级亚洲精品| 99久久中文字幕三级久久日本| 国模一区二区三区四区视频| 中文字幕免费在线视频6| 黄色一级大片看看| 日韩中字成人| 91精品伊人久久大香线蕉| 亚洲美女搞黄在线观看| 国产高清不卡午夜福利| 亚洲怡红院男人天堂| 天天躁日日操中文字幕| 如何舔出高潮| 欧美人与善性xxx| 午夜日本视频在线| 国产男女内射视频| 禁无遮挡网站| 99热这里只有是精品50| 美女脱内裤让男人舔精品视频| 七月丁香在线播放| 国产色爽女视频免费观看| 久久久久久久久久久免费av| 精品久久久久久电影网| 日韩欧美精品免费久久| 日本三级黄在线观看| 精品国产一区二区三区久久久樱花 | 欧美高清性xxxxhd video| 伦理电影大哥的女人| 日韩精品有码人妻一区| 大片免费播放器 马上看| av在线app专区| 亚洲欧美一区二区三区黑人 | 中文字幕免费在线视频6| 久久久久久久久久人人人人人人| 亚洲精品乱久久久久久| 亚洲精品日本国产第一区| 看黄色毛片网站| 国产 一区精品| 欧美xxxx黑人xx丫x性爽| 97在线人人人人妻| 国产午夜福利久久久久久| 又粗又硬又长又爽又黄的视频| 久久久久精品性色| kizo精华| 亚洲精品成人久久久久久| 色播亚洲综合网| 老师上课跳d突然被开到最大视频| 国产熟女欧美一区二区| av又黄又爽大尺度在线免费看| av女优亚洲男人天堂| 国产精品国产av在线观看| 亚洲自拍偷在线| 久久97久久精品| 极品少妇高潮喷水抽搐| 国产女主播在线喷水免费视频网站| 亚洲最大成人av| 国产精品蜜桃在线观看| 黄色日韩在线| 亚洲av男天堂| 好男人在线观看高清免费视频| 99久久精品一区二区三区| 亚洲真实伦在线观看| 综合色丁香网| 99久久人妻综合| 免费观看在线日韩| 又粗又硬又长又爽又黄的视频| 午夜激情福利司机影院| 搞女人的毛片| 少妇人妻一区二区三区视频| 午夜福利视频精品| 久久久a久久爽久久v久久| 日韩人妻高清精品专区| 深爱激情五月婷婷| 亚洲久久久久久中文字幕| 熟女电影av网| 麻豆国产97在线/欧美| 性色av一级| 黄色日韩在线| 久久人人爽人人爽人人片va| 亚洲欧美日韩另类电影网站 | 久久精品国产a三级三级三级| 大话2 男鬼变身卡| 欧美精品国产亚洲| 男人爽女人下面视频在线观看| 久久久精品免费免费高清| 精华霜和精华液先用哪个| 插阴视频在线观看视频| 午夜精品一区二区三区免费看| 一区二区三区精品91| 91精品伊人久久大香线蕉| 国产成人精品一,二区| 国产亚洲5aaaaa淫片| 亚洲精品影视一区二区三区av| 国产日韩欧美在线精品| 国产精品偷伦视频观看了| 亚洲国产精品成人综合色| 日本一本二区三区精品| 亚洲激情五月婷婷啪啪| 超碰av人人做人人爽久久| 午夜精品国产一区二区电影 | 美女主播在线视频| 好男人视频免费观看在线| 天堂网av新在线| 国产精品偷伦视频观看了| 亚洲高清免费不卡视频| 国产探花在线观看一区二区| 亚洲精品日韩在线中文字幕| av在线播放精品| 国产成人a区在线观看| 国产乱人视频| 女人十人毛片免费观看3o分钟| 日韩视频在线欧美| 久久99蜜桃精品久久| 最后的刺客免费高清国语| a级一级毛片免费在线观看| 亚洲成人中文字幕在线播放| 久久国内精品自在自线图片| 日韩欧美一区视频在线观看 | 亚洲人成网站在线播| 欧美日韩一区二区视频在线观看视频在线 | 美女主播在线视频| 国产成人a∨麻豆精品| 精品视频人人做人人爽| 国产精品国产三级专区第一集| 在线观看免费高清a一片| 欧美97在线视频| 日日啪夜夜撸| 美女内射精品一级片tv| 久热这里只有精品99| 亚洲精品日韩av片在线观看| 男插女下体视频免费在线播放| 午夜免费观看性视频| 26uuu在线亚洲综合色| 禁无遮挡网站| 欧美老熟妇乱子伦牲交| 少妇人妻久久综合中文| 午夜爱爱视频在线播放| 亚洲国产日韩一区二区| 2018国产大陆天天弄谢| 18禁动态无遮挡网站| 两个人的视频大全免费| 美女xxoo啪啪120秒动态图| 97在线人人人人妻| 在线观看av片永久免费下载| 色播亚洲综合网| 精品视频人人做人人爽| 久久精品国产自在天天线| 男人和女人高潮做爰伦理| 十八禁网站网址无遮挡 | 久久久欧美国产精品| 日韩制服骚丝袜av| 少妇熟女欧美另类| 最近中文字幕2019免费版| 能在线免费看毛片的网站| 国产成人精品一,二区| 亚洲欧美精品专区久久| 免费观看的影片在线观看| 亚洲欧美成人精品一区二区| 菩萨蛮人人尽说江南好唐韦庄| 国产精品99久久99久久久不卡 | 精品熟女少妇av免费看| 又粗又硬又长又爽又黄的视频| 韩国av在线不卡| 国产毛片a区久久久久| 亚洲精品色激情综合| 免费在线观看成人毛片| 国产精品一区二区三区四区免费观看| 免费电影在线观看免费观看| 国产黄色视频一区二区在线观看| 久久鲁丝午夜福利片| 日韩大片免费观看网站| 国产精品一区www在线观看| av在线亚洲专区| 精品国产三级普通话版| 丝袜脚勾引网站| av在线亚洲专区| 高清日韩中文字幕在线| 人妻系列 视频| 九九在线视频观看精品| 久久综合国产亚洲精品| 国产精品麻豆人妻色哟哟久久| 亚洲精品,欧美精品| 日韩免费高清中文字幕av| 国产毛片在线视频| 亚洲精品乱久久久久久| 国产精品久久久久久精品古装| 亚洲电影在线观看av| 黄色怎么调成土黄色| 80岁老熟妇乱子伦牲交| 久久99热这里只有精品18| 欧美丝袜亚洲另类| 亚洲成色77777| 国产精品精品国产色婷婷| 天堂网av新在线| 日韩中字成人| 少妇 在线观看| 全区人妻精品视频| 偷拍熟女少妇极品色| 欧美97在线视频| 老司机影院成人| 国产毛片在线视频| 99久久中文字幕三级久久日本| 免费看av在线观看网站| 制服丝袜香蕉在线| 亚洲欧美日韩卡通动漫| 亚洲av中文字字幕乱码综合| 国产成人免费观看mmmm| 蜜桃亚洲精品一区二区三区| 国产免费一级a男人的天堂| 午夜免费鲁丝| 久热这里只有精品99| 日日撸夜夜添| 蜜桃亚洲精品一区二区三区| 男女边摸边吃奶| 在现免费观看毛片| 波野结衣二区三区在线| 听说在线观看完整版免费高清| 99热这里只有是精品50| 国产一区二区在线观看日韩| 亚洲欧美一区二区三区国产| 成人毛片60女人毛片免费| 欧美激情在线99| 久久久久国产网址| 一级二级三级毛片免费看| 可以在线观看毛片的网站| 少妇被粗大猛烈的视频| 99热6这里只有精品| 成人毛片60女人毛片免费| 欧美成人a在线观看| 一个人看视频在线观看www免费| 精品一区在线观看国产| 老师上课跳d突然被开到最大视频| 成人黄色视频免费在线看| 青春草国产在线视频| 亚洲精品成人av观看孕妇| 一级毛片黄色毛片免费观看视频| 超碰av人人做人人爽久久| 一级毛片aaaaaa免费看小| 爱豆传媒免费全集在线观看| 久久久久久久大尺度免费视频| 汤姆久久久久久久影院中文字幕| 国产亚洲av嫩草精品影院| h日本视频在线播放| 综合色av麻豆| 精品视频人人做人人爽| 少妇丰满av| 久久精品久久久久久久性| 中文字幕av成人在线电影| 好男人在线观看高清免费视频| 波野结衣二区三区在线| 欧美激情在线99| 国产毛片在线视频| 久久精品综合一区二区三区| 国产色爽女视频免费观看| 久久99热这里只有精品18| 精品99又大又爽又粗少妇毛片| 91aial.com中文字幕在线观看| 少妇 在线观看| 91久久精品电影网| 99久国产av精品国产电影| av天堂中文字幕网| 日日摸夜夜添夜夜爱| 尾随美女入室| 听说在线观看完整版免费高清| 国产伦精品一区二区三区四那| 国产av不卡久久| 国产日韩欧美在线精品| 亚洲国产av新网站| 亚洲精品国产av蜜桃| 日本午夜av视频| a级毛色黄片| 日韩一本色道免费dvd| 色婷婷久久久亚洲欧美| 欧美三级亚洲精品| 最后的刺客免费高清国语| 国产黄片美女视频| 亚洲欧美日韩卡通动漫| 久久鲁丝午夜福利片| 91在线精品国自产拍蜜月| 久久鲁丝午夜福利片| 日本爱情动作片www.在线观看| 国产精品不卡视频一区二区| 精品一区在线观看国产| 免费看不卡的av| 丰满乱子伦码专区| 夜夜爽夜夜爽视频| 亚洲国产欧美在线一区| 男女那种视频在线观看| 国产精品久久久久久精品电影| 九九爱精品视频在线观看| 26uuu在线亚洲综合色| 中文字幕亚洲精品专区| 哪个播放器可以免费观看大片| 久久99热6这里只有精品| 国产精品成人在线| 超碰97精品在线观看| 国产真实伦视频高清在线观看| 成人鲁丝片一二三区免费| freevideosex欧美| 99热全是精品| 国产亚洲5aaaaa淫片| 国产乱来视频区| 深夜a级毛片| 晚上一个人看的免费电影| 免费看a级黄色片| 国产视频内射| av在线观看视频网站免费| 麻豆精品久久久久久蜜桃| 一区二区三区精品91| 麻豆精品久久久久久蜜桃| 蜜臀久久99精品久久宅男| 三级国产精品片| 91aial.com中文字幕在线观看| 国产精品久久久久久久久免| 亚洲精品日韩av片在线观看| 免费看日本二区| 日韩av不卡免费在线播放| 日韩 亚洲 欧美在线| 国产永久视频网站| 国产美女午夜福利| 日韩人妻高清精品专区| 亚洲精品成人av观看孕妇| 少妇人妻一区二区三区视频| 亚洲成人一二三区av| 久久99精品国语久久久| 久久久久性生活片| 日本与韩国留学比较| 女人被狂操c到高潮| 成人黄色视频免费在线看| 欧美极品一区二区三区四区| 亚洲欧洲国产日韩| 亚洲精品一区蜜桃| 色吧在线观看| 日本一本二区三区精品| 2018国产大陆天天弄谢| 成人午夜精彩视频在线观看| 亚洲av电影在线观看一区二区三区 | 成年免费大片在线观看| 五月伊人婷婷丁香| 99久久人妻综合| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 看免费成人av毛片| 亚洲国产精品成人综合色| av天堂中文字幕网| 美女cb高潮喷水在线观看| 国内精品美女久久久久久| 国产成人aa在线观看| 免费看a级黄色片| 99热6这里只有精品| 欧美日韩在线观看h| 欧美高清性xxxxhd video| 国产精品久久久久久精品电影| 夜夜爽夜夜爽视频| 久久久精品94久久精品| 中文资源天堂在线| 国产综合精华液| 麻豆成人av视频| 久久精品久久精品一区二区三区| 欧美激情久久久久久爽电影| 最后的刺客免费高清国语| 亚洲精品aⅴ在线观看| 日韩 亚洲 欧美在线| 亚洲欧美精品专区久久| 成年av动漫网址| 最近最新中文字幕免费大全7| 亚洲av成人精品一区久久| 国产一区二区亚洲精品在线观看| 欧美日本视频| 一区二区三区四区激情视频| 在线天堂最新版资源| 水蜜桃什么品种好| 日本与韩国留学比较| 狠狠精品人妻久久久久久综合| 亚洲国产精品专区欧美| 大又大粗又爽又黄少妇毛片口|