• <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ù)年代際變化特征
    欧美老熟妇乱子伦牲交| www.精华液| 久久国产亚洲av麻豆专区| 下体分泌物呈黄色| 国产野战对白在线观看| 亚洲人成77777在线视频| 欧美人与性动交α欧美精品济南到 | 精品国产一区二区久久| 毛片一级片免费看久久久久| 丝袜在线中文字幕| 国产在线免费精品| 免费在线观看黄色视频的| 国产视频首页在线观看| 一边摸一边做爽爽视频免费| 亚洲第一青青草原| 久久99精品国语久久久| 丰满迷人的少妇在线观看| 日韩欧美精品免费久久| 亚洲精品一区蜜桃| 精品国产一区二区三区四区第35| 满18在线观看网站| 少妇的丰满在线观看| 精品卡一卡二卡四卡免费| 中文字幕亚洲精品专区| 国产免费现黄频在线看| 亚洲av.av天堂| 亚洲成人手机| 老汉色∧v一级毛片| a级片在线免费高清观看视频| 国产女主播在线喷水免费视频网站| 少妇人妻 视频| 免费观看性生交大片5| 国产精品国产三级国产专区5o| 丁香六月天网| 边亲边吃奶的免费视频| 精品一区在线观看国产| 久久久久久免费高清国产稀缺| 亚洲精华国产精华液的使用体验| 中文欧美无线码| 国产精品一区二区在线不卡| 999精品在线视频| 人妻一区二区av| 国产黄色免费在线视频| 国产亚洲最大av| 欧美日韩av久久| 美女午夜性视频免费| 久久人妻熟女aⅴ| 成人国产av品久久久| 中国三级夫妇交换| 国产精品女同一区二区软件| 18禁动态无遮挡网站| 亚洲在久久综合| 久久99蜜桃精品久久| 精品第一国产精品| 久久久久久久国产电影| 人体艺术视频欧美日本| 日韩av免费高清视频| 欧美国产精品va在线观看不卡| 老女人水多毛片| 丝袜美腿诱惑在线| 色吧在线观看| 日韩免费高清中文字幕av| 在线观看人妻少妇| 欧美xxⅹ黑人| 中国国产av一级| 一边摸一边做爽爽视频免费| 欧美人与善性xxx| 国产精品一国产av| 99久久综合免费| 伊人久久国产一区二区| 国产女主播在线喷水免费视频网站| 中国三级夫妇交换| 亚洲五月色婷婷综合| 亚洲人成网站在线观看播放| 国产黄色免费在线视频| 国产1区2区3区精品| 中文字幕亚洲精品专区| 叶爱在线成人免费视频播放| 高清av免费在线| 99久久综合免费| 99久久人妻综合| 黑丝袜美女国产一区| 久久这里只有精品19| 亚洲欧美一区二区三区久久| 男女无遮挡免费网站观看| 高清黄色对白视频在线免费看| 亚洲一级一片aⅴ在线观看| 欧美人与性动交α欧美软件| 在线观看免费视频网站a站| 精品人妻偷拍中文字幕| 久久综合国产亚洲精品| 黑人巨大精品欧美一区二区蜜桃| 免费高清在线观看日韩| 日本猛色少妇xxxxx猛交久久| 最近的中文字幕免费完整| 蜜桃在线观看..| 欧美精品一区二区免费开放| 久久精品国产综合久久久| av在线老鸭窝| 最新的欧美精品一区二区| 久久久久久久久免费视频了| 欧美最新免费一区二区三区| 国产成人精品婷婷| 国产日韩欧美视频二区| www.av在线官网国产| 久久久精品国产亚洲av高清涩受| 99久国产av精品国产电影| 精品国产一区二区三区四区第35| 2022亚洲国产成人精品| 亚洲,一卡二卡三卡| 久久精品国产亚洲av天美| 熟女电影av网| 久久影院123| 另类精品久久| 国产av精品麻豆| xxx大片免费视频| 国产黄色视频一区二区在线观看| 日韩,欧美,国产一区二区三区| 最新中文字幕久久久久| 成年人午夜在线观看视频| 欧美变态另类bdsm刘玥| 搡女人真爽免费视频火全软件| 精品久久久精品久久久| 亚洲视频免费观看视频| 精品久久久久久电影网| 少妇 在线观看| 欧美人与善性xxx| 日韩不卡一区二区三区视频在线| 亚洲四区av| 成人午夜精彩视频在线观看| 欧美日韩一级在线毛片| 女性被躁到高潮视频| 另类亚洲欧美激情| 国产精品女同一区二区软件| 国产在视频线精品| 成人午夜精彩视频在线观看| 少妇人妻 视频| av不卡在线播放| 欧美在线黄色| 如日韩欧美国产精品一区二区三区| 国产av国产精品国产| 女人高潮潮喷娇喘18禁视频| 母亲3免费完整高清在线观看 | 国产精品久久久av美女十八| 国产精品一国产av| 国产片特级美女逼逼视频| 国产精品国产av在线观看| 亚洲美女黄色视频免费看| 在线观看免费高清a一片| 日韩av免费高清视频| 亚洲五月色婷婷综合| 亚洲精品美女久久av网站| 桃花免费在线播放| 日韩中文字幕欧美一区二区 | a级片在线免费高清观看视频| 婷婷色综合大香蕉| 国产精品久久久久成人av| 两个人免费观看高清视频| 久久精品国产综合久久久| 亚洲欧美色中文字幕在线| 1024视频免费在线观看| 热re99久久国产66热| 精品国产一区二区三区四区第35| 欧美人与性动交α欧美精品济南到 | 男人舔女人的私密视频| 久久免费观看电影| 精品卡一卡二卡四卡免费| 亚洲精品久久成人aⅴ小说| 日本av免费视频播放| 久久精品国产综合久久久| 精品一品国产午夜福利视频| 少妇人妻 视频| 亚洲国产精品成人久久小说| 亚洲中文av在线| 成年女人毛片免费观看观看9 | 亚洲精品国产色婷婷电影| 午夜精品国产一区二区电影| 亚洲一区二区三区欧美精品| 午夜老司机福利剧场| 久久这里有精品视频免费| 哪个播放器可以免费观看大片| av电影中文网址| 日韩在线高清观看一区二区三区| 亚洲一级一片aⅴ在线观看| 韩国av在线不卡| 色94色欧美一区二区| 国产激情久久老熟女| 亚洲成人一二三区av| 99国产综合亚洲精品| 男女高潮啪啪啪动态图| 波多野结衣av一区二区av| 超色免费av| 精品久久久精品久久久| 久久鲁丝午夜福利片| 亚洲综合色网址| 丝袜在线中文字幕| 亚洲激情五月婷婷啪啪| 国产精品久久久久久av不卡| 欧美老熟妇乱子伦牲交| 久久久国产一区二区| 丝袜在线中文字幕| 亚洲一区二区三区欧美精品| 亚洲四区av| 91精品三级在线观看| 免费不卡的大黄色大毛片视频在线观看| 久久人人爽人人片av| 色94色欧美一区二区| 丰满乱子伦码专区| 久久久久久久亚洲中文字幕| 一区在线观看完整版| 热99国产精品久久久久久7| 亚洲成av片中文字幕在线观看 | 久久久久久人妻| 国产女主播在线喷水免费视频网站| 久久青草综合色| 国产免费又黄又爽又色| 亚洲美女视频黄频| 一区二区三区激情视频| 精品福利永久在线观看| 婷婷色综合www| 欧美+日韩+精品| 国产一区亚洲一区在线观看| www.自偷自拍.com| 亚洲熟女精品中文字幕| 最近的中文字幕免费完整| 欧美最新免费一区二区三区| av国产精品久久久久影院| 亚洲国产看品久久| 国产成人免费观看mmmm| 三上悠亚av全集在线观看| 纵有疾风起免费观看全集完整版| 欧美日韩国产mv在线观看视频| 亚洲欧美成人精品一区二区| 国产亚洲欧美精品永久| 男女午夜视频在线观看| 国产极品粉嫩免费观看在线| 精品国产超薄肉色丝袜足j| 肉色欧美久久久久久久蜜桃| 国产1区2区3区精品| 少妇的逼水好多| 极品人妻少妇av视频| 久久精品久久久久久久性| 欧美日韩一区二区视频在线观看视频在线| 国产免费视频播放在线视频| 少妇被粗大的猛进出69影院| 国产成人免费观看mmmm| 黄色 视频免费看| 午夜福利网站1000一区二区三区| 美国免费a级毛片| 中文字幕人妻丝袜一区二区 | 九色亚洲精品在线播放| 寂寞人妻少妇视频99o| 亚洲 欧美一区二区三区| 日韩电影二区| 菩萨蛮人人尽说江南好唐韦庄| 国产精品一二三区在线看| 中文字幕人妻丝袜制服| 人人妻人人添人人爽欧美一区卜| 国产成人精品一,二区| 99国产综合亚洲精品| 国产精品欧美亚洲77777| 成年人午夜在线观看视频| 日日撸夜夜添| 日日摸夜夜添夜夜爱| 国产精品国产三级国产专区5o| 考比视频在线观看| 一级毛片黄色毛片免费观看视频| 色吧在线观看| 久久久久久久亚洲中文字幕| 一级片'在线观看视频| 成年人午夜在线观看视频| 欧美在线黄色| 亚洲国产精品一区二区三区在线| 精品一区二区免费观看| 久久精品夜色国产| 美女中出高潮动态图| 视频在线观看一区二区三区| 久久久精品94久久精品| 欧美日韩精品网址| av有码第一页| 中文字幕亚洲精品专区| 午夜精品国产一区二区电影| 成人亚洲欧美一区二区av| 色吧在线观看| 免费观看无遮挡的男女| 国产人伦9x9x在线观看 | 丝袜在线中文字幕| 久久久久久久久久久免费av| 免费不卡的大黄色大毛片视频在线观看| 看免费av毛片| 三级国产精品片| 久久这里有精品视频免费| 午夜福利网站1000一区二区三区| www.精华液| 中文字幕制服av| 99国产综合亚洲精品| 大码成人一级视频| av电影中文网址| av女优亚洲男人天堂| 日韩电影二区| 国产欧美亚洲国产| 欧美日韩视频精品一区| 91精品伊人久久大香线蕉| 精品亚洲成国产av| 交换朋友夫妻互换小说| 黄色怎么调成土黄色| 国产亚洲精品第一综合不卡| 精品人妻偷拍中文字幕| 亚洲色图综合在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲av欧美aⅴ国产| 人人澡人人妻人| 精品亚洲成a人片在线观看| 亚洲第一区二区三区不卡| 亚洲精品国产色婷婷电影| 男人舔女人的私密视频| 交换朋友夫妻互换小说| 色视频在线一区二区三区| 99久久人妻综合| 久久免费观看电影| 久久99热这里只频精品6学生| 精品一品国产午夜福利视频| 男女边吃奶边做爰视频| 亚洲精品久久成人aⅴ小说| 久久久久国产一级毛片高清牌| 国产在线一区二区三区精| 日韩中字成人| 狠狠精品人妻久久久久久综合| 中文欧美无线码| 国产日韩欧美视频二区| 成人手机av| 亚洲精品第二区| 男女下面插进去视频免费观看| 久久国内精品自在自线图片| 黄频高清免费视频| 在线 av 中文字幕| 天天操日日干夜夜撸| 最近手机中文字幕大全| 亚洲欧美精品自产自拍| 女性被躁到高潮视频| 日本-黄色视频高清免费观看| 天天躁日日躁夜夜躁夜夜| 国产免费现黄频在线看| videossex国产| 波野结衣二区三区在线| 国产免费视频播放在线视频| 国产精品一区二区在线观看99| 日日爽夜夜爽网站| 国产成人精品在线电影| 一区二区三区乱码不卡18| 波野结衣二区三区在线| 美国免费a级毛片| 99国产综合亚洲精品| 永久网站在线| 老司机亚洲免费影院| 熟女电影av网| 老司机亚洲免费影院| 天堂中文最新版在线下载| 人妻人人澡人人爽人人| 寂寞人妻少妇视频99o| 欧美人与性动交α欧美软件| 永久免费av网站大全| 男人爽女人下面视频在线观看| 亚洲欧美精品自产自拍| 亚洲一区中文字幕在线| 男女边摸边吃奶| 亚洲欧美一区二区三区国产| 国产极品天堂在线| 另类精品久久| 国产欧美日韩综合在线一区二区| 国产在线一区二区三区精| 菩萨蛮人人尽说江南好唐韦庄| 久久精品国产鲁丝片午夜精品| 日本wwww免费看| 久久人人97超碰香蕉20202| 国产欧美日韩综合在线一区二区| 日日啪夜夜爽| 日韩三级伦理在线观看| 亚洲精品久久成人aⅴ小说| 国产激情久久老熟女| 国产精品香港三级国产av潘金莲 | 日本午夜av视频| 久久久久久人人人人人| 少妇人妻 视频| 国产在线一区二区三区精| 日本vs欧美在线观看视频| 新久久久久国产一级毛片| 18禁裸乳无遮挡动漫免费视频| av在线老鸭窝| 色哟哟·www| 两个人免费观看高清视频| 精品一区在线观看国产| 91精品三级在线观看| 国精品久久久久久国模美| 久久精品人人爽人人爽视色| 一区二区三区乱码不卡18| 中文字幕人妻丝袜制服| 黄色视频在线播放观看不卡| 免费人妻精品一区二区三区视频| 女的被弄到高潮叫床怎么办| 亚洲国产av新网站| 欧美精品高潮呻吟av久久| 又粗又硬又长又爽又黄的视频| 日本欧美视频一区| 精品一区在线观看国产| 老汉色av国产亚洲站长工具| av国产久精品久网站免费入址| 99精国产麻豆久久婷婷| 五月天丁香电影| 老鸭窝网址在线观看| 最近2019中文字幕mv第一页| 成人手机av| 中文字幕亚洲精品专区| 色吧在线观看| 精品人妻在线不人妻| 成人影院久久| 久久鲁丝午夜福利片| 亚洲成人手机| 在线看a的网站| 国产成人aa在线观看| 超色免费av| 欧美成人精品欧美一级黄| 日韩不卡一区二区三区视频在线| 国产精品久久久av美女十八| 亚洲欧洲精品一区二区精品久久久 | 欧美成人午夜免费资源| 最近最新中文字幕免费大全7| 亚洲第一av免费看| 男女边吃奶边做爰视频| 十八禁网站网址无遮挡| 好男人视频免费观看在线| 成年美女黄网站色视频大全免费| 国产黄色视频一区二区在线观看| 亚洲色图 男人天堂 中文字幕| 中文字幕亚洲精品专区| 少妇的丰满在线观看| 午夜免费观看性视频| 丰满迷人的少妇在线观看| 久久久久国产网址| 春色校园在线视频观看| 国产精品 国内视频| 在线观看www视频免费| 汤姆久久久久久久影院中文字幕| 99热全是精品| 下体分泌物呈黄色| 2022亚洲国产成人精品| 狠狠精品人妻久久久久久综合| 一级毛片 在线播放| 国产xxxxx性猛交| 国产视频首页在线观看| 国产一区二区激情短视频 | 欧美精品一区二区大全| 久久久精品区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 国产免费又黄又爽又色| 黑人欧美特级aaaaaa片| 久久久久久人妻| 国产深夜福利视频在线观看| 只有这里有精品99| 天天躁夜夜躁狠狠躁躁| 大码成人一级视频| 欧美黄色片欧美黄色片| videossex国产| 午夜福利网站1000一区二区三区| 精品国产超薄肉色丝袜足j| 人人妻人人澡人人看| 伊人久久国产一区二区| 一区二区av电影网| 国产精品久久久久久精品古装| 国产有黄有色有爽视频| 久久鲁丝午夜福利片| 一级毛片电影观看| 精品人妻一区二区三区麻豆| 青草久久国产| 天美传媒精品一区二区| 日韩欧美精品免费久久| 80岁老熟妇乱子伦牲交| 亚洲综合色网址| 国产日韩欧美亚洲二区| 精品国产一区二区久久| 女人被躁到高潮嗷嗷叫费观| 汤姆久久久久久久影院中文字幕| 亚洲欧美一区二区三区国产| 母亲3免费完整高清在线观看 | 欧美成人午夜精品| 男女免费视频国产| 高清欧美精品videossex| 色94色欧美一区二区| 91精品伊人久久大香线蕉| 国产在线免费精品| 看免费av毛片| 久久鲁丝午夜福利片| 国产一区二区 视频在线| 1024香蕉在线观看| 国产精品欧美亚洲77777| 男女无遮挡免费网站观看| 日韩av免费高清视频| 国产免费又黄又爽又色| 精品一区在线观看国产| 国产成人精品在线电影| 亚洲欧美精品综合一区二区三区 | 亚洲精品成人av观看孕妇| 免费在线观看黄色视频的| 国产色婷婷99| 9色porny在线观看| 黄色配什么色好看| 2021少妇久久久久久久久久久| freevideosex欧美| 十八禁高潮呻吟视频| 免费av中文字幕在线| 亚洲在久久综合| 国产精品久久久久成人av| 欧美日韩亚洲高清精品| av片东京热男人的天堂| 考比视频在线观看| 大片免费播放器 马上看| 日韩人妻精品一区2区三区| 晚上一个人看的免费电影| 亚洲一级一片aⅴ在线观看| 久久亚洲国产成人精品v| 午夜福利视频精品| 日韩一本色道免费dvd| 久久人人97超碰香蕉20202| 国产激情久久老熟女| 岛国毛片在线播放| 人妻少妇偷人精品九色| 99热国产这里只有精品6| 美女中出高潮动态图| 亚洲色图综合在线观看| 欧美精品一区二区大全| 男人舔女人的私密视频| 一级毛片电影观看| 久久国产精品大桥未久av| 侵犯人妻中文字幕一二三四区| 精品国产一区二区久久| 亚洲视频免费观看视频| 一区二区三区乱码不卡18| 只有这里有精品99| 欧美av亚洲av综合av国产av | 999精品在线视频| av又黄又爽大尺度在线免费看| 精品少妇黑人巨大在线播放| 欧美 亚洲 国产 日韩一| 久久99蜜桃精品久久| tube8黄色片| 天美传媒精品一区二区| 国产成人精品无人区| 三上悠亚av全集在线观看| www.精华液| 亚洲经典国产精华液单| 亚洲五月色婷婷综合| 亚洲av.av天堂| 人妻一区二区av| 不卡av一区二区三区| 观看av在线不卡| 亚洲三级黄色毛片| 高清不卡的av网站| 国产精品 国内视频| 久久精品国产鲁丝片午夜精品| 制服人妻中文乱码| 一级毛片黄色毛片免费观看视频| 国产无遮挡羞羞视频在线观看| 中文精品一卡2卡3卡4更新| 九草在线视频观看| 狂野欧美激情性bbbbbb| 美女国产视频在线观看| 免费在线观看完整版高清| 我要看黄色一级片免费的| 久久99蜜桃精品久久| 免费少妇av软件| 国产成人免费无遮挡视频| 午夜日本视频在线| tube8黄色片| av在线播放精品| 日韩一本色道免费dvd| av在线app专区| 观看美女的网站| 成人毛片a级毛片在线播放| 亚洲图色成人| 丰满迷人的少妇在线观看| 国产免费福利视频在线观看| 中文字幕亚洲精品专区| www.熟女人妻精品国产| 丝袜美腿诱惑在线| 香蕉精品网在线| 一级毛片电影观看| 欧美激情高清一区二区三区 | 五月开心婷婷网| 麻豆乱淫一区二区| 亚洲成国产人片在线观看| 欧美日本中文国产一区发布| 欧美日韩av久久| 亚洲av男天堂| 人妻一区二区av| a 毛片基地| 99热网站在线观看| 亚洲在久久综合| 最近2019中文字幕mv第一页| 亚洲一级一片aⅴ在线观看| 黄色配什么色好看| 久久人人爽人人片av| 在线观看三级黄色| 色94色欧美一区二区| 妹子高潮喷水视频| 国产精品免费大片| 99国产综合亚洲精品| 飞空精品影院首页| 伦理电影免费视频| 午夜福利在线观看免费完整高清在| 午夜精品国产一区二区电影| 久久久久人妻精品一区果冻| 满18在线观看网站| 女人高潮潮喷娇喘18禁视频| 91在线精品国自产拍蜜月| 这个男人来自地球电影免费观看 |