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

    降水變化和人類活動(dòng)對(duì)瀏陽河流域徑流的影響

    2017-11-28 03:34:55,,,
    長江科學(xué)院院報(bào) 2017年11期
    關(guān)鍵詞:人類

    ,, ,

    (1.長沙理工大學(xué) 水利工程學(xué)院,長沙 410114;2.湖南省水文水資源勘測(cè)局,長沙 410007)

    2017,34(11):6-11

    降水變化和人類活動(dòng)對(duì)瀏陽河流域徑流的影響

    于澤興1,胡國華1,陳肖1,顧慶福2

    (1.長沙理工大學(xué) 水利工程學(xué)院,長沙 410114;2.湖南省水文水資源勘測(cè)局,長沙 410007)

    定量評(píng)估降水變化和人類活動(dòng)對(duì)瀏陽河流域徑流的影響,有助于更合理地利用流域水資源。采用線性擬合、滑動(dòng)平均和累積距平等方法分析了1969—2013年瀏陽河流域徑流和降水的演變規(guī)律及其相關(guān)關(guān)系,應(yīng)用累積斜率變化率比較法定量估算了不同時(shí)期降水變化和人類活動(dòng)對(duì)徑流量變化的貢獻(xiàn)率。結(jié)果表明:瀏陽河流域降水量表現(xiàn)為增加趨勢(shì),上升率為3 mm/(10 a);年徑流量表現(xiàn)為下降趨勢(shì),下降率為2.7 (m3·s-1)/(10 a);降水和徑流在1988年和1999年均發(fā)生顯著性突變,經(jīng)歷了“減少—增多—減少”的交替變換。與基準(zhǔn)期1969—1988年相比,1989—1998年和1999—2013年這2個(gè)時(shí)期的降水量變化對(duì)徑流量變化貢獻(xiàn)率依次為56.8%,7.62%,而人類活動(dòng)對(duì)徑流量的貢獻(xiàn)率依次為43.2%,92.38%,可見人類活動(dòng)已成為影響瀏陽河流域徑流變化的主要因素。

    瀏陽河流域;徑流量變化;累積斜率變化率比較法;降水變化;人類活動(dòng)

    1 研究背景

    降水作為氣候變化的主要因素之一,對(duì)水資源的變化起著決定性作用。近幾十年來,全球氣候變暖和大規(guī)模的人類活動(dòng)引起了水資源在時(shí)空上的重新分布,導(dǎo)致河川徑流發(fā)生了明顯變化,對(duì)流域的綜合管理和生態(tài)經(jīng)濟(jì)的可持續(xù)發(fā)展產(chǎn)生了顯著影響,而水資源問題會(huì)引發(fā)一系列的生態(tài)環(huán)境問題[1-4]。對(duì)于不同地區(qū)和流域,受人類生產(chǎn)力和地理?xiàng)l件的限制,影響水資源的主要因素可能會(huì)有所區(qū)別。因此,分析降水變化和徑流的相互關(guān)系以及定量區(qū)分氣候和人類活動(dòng)對(duì)徑流的影響,已成為目前研究的熱點(diǎn)問題之一。

    近年來,研究者對(duì)湘江流域及其支流的降水量、徑流量的演變規(guī)律及其影響因素進(jìn)行了相關(guān)研究[5-8],認(rèn)為湘江流域的徑流量總體呈現(xiàn)減少趨勢(shì),降水量總體呈現(xiàn)增多趨勢(shì),趨勢(shì)的顯著性存在明顯的時(shí)空差異,并且降水量存在明顯的階段性變化。張海斌[9]發(fā)現(xiàn)瀏陽河流域的暴雨中心多在流域上游,小溪河山地降水量多而下游的河口平原地區(qū)降水明顯偏少。張棋等[10]利用L-THIA水文模型分析了瀏陽河流域土地利用變化對(duì)徑流量以及非點(diǎn)源污染負(fù)荷的影響。然而,這些研究很少涉及徑流變化的深層次原因及其主要影響因子貢獻(xiàn)率的定量估算。

    瀏陽河流域經(jīng)濟(jì)發(fā)達(dá),流域內(nèi)長沙縣和瀏陽市位于全國縣域經(jīng)濟(jì)前列。瀏陽河的株樹橋水庫作為長沙市的水源地之一,影響著長沙市區(qū)和瀏陽市數(shù)百萬人口的用水和工農(nóng)業(yè)生產(chǎn)。自20世紀(jì)60年代以來,由于氣候變化和人類活動(dòng)的影響,瀏陽河流域的徑流量也發(fā)生了明顯的階段性變化,并導(dǎo)致了水土流失、濕地退化等嚴(yán)重的生態(tài)環(huán)境問題,流域水資源供需矛盾日益加劇,并且旱澇災(zāi)害頻發(fā),嚴(yán)重影響瀏陽的農(nóng)業(yè)生產(chǎn)和社會(huì)經(jīng)濟(jì)發(fā)展。本文基于線性擬合、滑動(dòng)平均和累積距平等統(tǒng)計(jì)方法,在分析降水和徑流的演變規(guī)律的基礎(chǔ)上,采用累積量斜率變化率比較法[11]定量評(píng)估氣候變化和人類活動(dòng)對(duì)徑流量變化的影響,從而為瀏陽河流域的水資源可持續(xù)利用和管護(hù)提供有意義的參考。

    2 研究區(qū)、數(shù)據(jù)及方法

    2.1 研究區(qū)概況

    瀏陽河是湘江的一級(jí)支流,流經(jīng)瀏陽市、長沙縣及長沙市區(qū),位于112°58′E —114°15′E,27°51′N —28°34′N之間。全長234 km,流域面積4 665 km2,東西走向,形狀狹長,地勢(shì)東北高西南低(圖1)。地貌以山地丘陵為主,屬亞熱帶濕潤季風(fēng)氣候,多年平均氣溫17.6 ℃;多年平均降水量1 601.1 mm,以大圍山一帶為多雨區(qū)。流域出口控制站為郎梨水文站(113°5′E,28°10′N),郎梨站以上集水面積約為3 815 km2,多年平均徑流量為99.7 m3/s。近年來瀏陽市和長沙的經(jīng)濟(jì)飛速發(fā)展,城鎮(zhèn)化進(jìn)程加速和人口增多引起用水量增加,瀏陽河流域下墊面受人類活動(dòng)的影響也越來越強(qiáng)。

    圖1瀏陽河流域位置

    Fig.1LocationofLiuyangRiverwatershed

    2.2 數(shù)據(jù)來源

    本研究利用的原始數(shù)據(jù)為郎梨站1969—2013年實(shí)測(cè)徑流資料和瀏陽河流域內(nèi)12個(gè)雨量站的逐日降水?dāng)?shù)據(jù),這些數(shù)據(jù)來源于湖南省水文與水資源勘測(cè)局。基于ArcGIS 10.2平臺(tái)對(duì)12個(gè)雨量站點(diǎn)逐月降水?dāng)?shù)據(jù),利用泰森多邊形法插值成水文站控制區(qū)域內(nèi)的面降水量,且序列長度與徑流保持一致。

    2.3 研究方法

    采用線性回歸和5 a滑動(dòng)平均法[12]對(duì)距平處理后的降水、徑流時(shí)間序列進(jìn)行了年際、年代際和季節(jié)的變化分析。利用累計(jì)距平[13]和滑動(dòng)t檢驗(yàn)[14]分析標(biāo)準(zhǔn)化后的徑流、降水序列的突變,并結(jié)合分段線性擬合分析二者的相互關(guān)系。

    在定量分析氣候變化和人類活動(dòng)對(duì)徑流的影響時(shí),降水與徑流雙累積曲線法[15]在人類活動(dòng)影響較小、年降水量和徑流量同比例增加或減少時(shí),其雙累積曲線在拐點(diǎn)前后可能仍然保持在同一直線上,從而造成拐點(diǎn)劃分的隨意性,因此本文采用累積量斜率變化率法[11]。其原理是假如徑流量只受到降水變化的影響,則降水和徑流隨年份變化的累積斜率應(yīng)是同倍比變化的,然后根據(jù)各個(gè)影響因素隨時(shí)間累積斜率的變化率占變量累積斜率變化率的比值來推求其影響程度。假設(shè)累積徑流量隨年份的關(guān)系式的斜率在變點(diǎn)前后分別為SRb和SRa(億m3/a), 累積降水量隨年份的關(guān)系式的斜率在變點(diǎn)前后分別為SPb和SPa(mm/a),則累積徑流量斜率變化率RSR為

    RSR=100%×(SRa-SRb)/SRb。

    (1)

    累積降水量斜率變化率為RSP為

    RSP=100%×(SPa-SPb)/SPb。

    (2)

    降水量變化對(duì)徑流量變化的貢獻(xiàn)率CP可以表示為

    CP=100%×RSP/RSR。

    (3)

    氣溫的變化也會(huì)影響流域蒸發(fā)量的改變,從而影響徑流量的變化。如果用CET表示蒸發(fā)量對(duì)流域徑流變化的貢獻(xiàn)率,則人類活動(dòng)對(duì)徑流量變化的貢獻(xiàn)率CH為

    CH=100%-CP-CET。

    (4)

    張婷婷[16]的研究表明,瀏陽位于山地東側(cè)的迎風(fēng)坡,降水豐富,氣候濕潤,風(fēng)速較小導(dǎo)致蒸發(fā)皿蒸發(fā)微弱。從湘江流域蒸發(fā)皿蒸發(fā)量的氣候傾向率空間分布來看,瀏陽站點(diǎn)的蒸發(fā)皿蒸發(fā)量無明顯的變化趨勢(shì),未通過置信度檢驗(yàn)。說明該地區(qū)蒸發(fā)量變化對(duì)徑流量的影響非常有限,因此,假定CET近似為0。

    3 降水量和徑流的變化特征及關(guān)系

    3.1 降水量的年際、年代際和季節(jié)變化

    表1給出了1969—2013年瀏陽河流域年、季節(jié)降水量平均值、標(biāo)準(zhǔn)差和變差系數(shù)。由表1可知,瀏陽河流域1969—2013多年平均降水量為1 601.1 mm,呈上升趨勢(shì),上升率為3 mm/(10 a),標(biāo)準(zhǔn)差259.7 mm,變差系數(shù)0.16。

    表11969—2013年瀏陽河流域年、季節(jié)面降水量平均值、標(biāo)準(zhǔn)差和變差系數(shù)

    Table1Meanvalues,standarddeviationsandvariationcoefficientsofannualandseasonalprecipitationfrom1969to2013inLiuyangRiverwatershed

    統(tǒng)計(jì)時(shí)段平均值/mm標(biāo)準(zhǔn)差/mm變差系數(shù)年度1601.1257.90.16春季594.4126.70.21夏季551.0188.70.34秋季225.999.10.44冬季229.866.00.29

    從降水的季節(jié)分布來看,春、夏季偏多,而秋、冬季偏少。春季和夏季為瀏陽河流域面降水量最多的季節(jié),降水量為594.4,551 mm,分別占年降水量的37.1%,34.4%。而秋、冬季面降水量明顯偏少,分別只有225.9,229.8 mm,分別占年降水量的14.1%,14.4%。從季節(jié)的變差系數(shù)可知,春季最小,最穩(wěn)定;夏、冬季次之;秋季最大,穩(wěn)定性最差。

    由瀏陽河流域1969—2013年降水量距平及5 a滑動(dòng)平均曲線(圖2)可知,瀏陽河流域年降水量有明顯的年際、年代際的變化特征。最大年降水量為1998年的2 241.1 mm,最小降水量為2007年的1 082.8 mm。20世紀(jì)70年代、90年代和21世紀(jì)初年際間的降水量變化較大,變差系數(shù)在0.15~0.23之間;而20世紀(jì)80年代的年際變化較小,變差系數(shù)在0.08左右。20世紀(jì)70年代降水處于波動(dòng)的平穩(wěn)階段,90年代處于降水偏多期,80年代和21世紀(jì)初處于降水偏少期。

    圖2瀏陽河流域1969—2013年降水量距平變化

    Fig.2Variationofannualprecipitationanomaliesfrom1969to2013inLiuyangRiverwatershed

    3.2 徑流的年際、年代際和季節(jié)變化特征

    由表2可知,郎梨水文站1969—2013多年徑流量的平均值為99.7 m3/s,標(biāo)準(zhǔn)差為29.4 m3/s,變差系數(shù)為0.29,呈減少趨勢(shì),下降率為2.7 (m3·s-1)/(10 a)。

    表2 1969—2013年瀏陽河流域年、季節(jié)徑流量平均值、標(biāo)準(zhǔn)差和變差系數(shù)

    從季節(jié)分布來看,春季和夏季的徑流量分別為155.1,146 m3/s,分別占到全年的38.9%,36.6%。秋季和冬季較小,分別為48,49.8 m3/s。從季節(jié)的變差系數(shù)可知,夏季的變差系數(shù)最大,為0.53,波動(dòng)最大,秋季和冬季次之,春季變差系數(shù)最小,為0.35,變化相對(duì)穩(wěn)定。

    由瀏陽河流域1969—2013年徑流量距平及5 a滑動(dòng)平均曲線(圖3)可知,郎梨的年徑流量也有明顯的年際、年代際變化特征。最大年徑流量為1998年的191.5 m3/s,最小年徑流量為2007年的45 m3/s,與面降水量極值年份一致。20世紀(jì)70年代、90年代和21世紀(jì)初年際波動(dòng)較大,變差系數(shù)在0.25~0.33之間;而20世紀(jì)80年代的年際變化較平緩,變差系數(shù)為0.19。徑流量按豐、平、枯水年劃分為5個(gè)等級(jí)[17],頻率分別為(87.5%,100%),(62.5%,87.5%),(37.5%,62.5%),(12.5%,37.5%),(0,12.5%) 的徑流分別對(duì)應(yīng)枯水年、偏枯水年、平水年、偏豐水年和豐水年。20世紀(jì)70年代至90年代未出現(xiàn)枯水年,而21世紀(jì)初有3年出現(xiàn)枯水年,分別為2007,2004,2011年。1969,1973,1998年為豐水年,平水年的年數(shù)為12。偏豐和偏枯的年數(shù)分別為14,13,占總年數(shù)的31.1%,28.9%。從年代際變化來看,20世紀(jì)90年代處于偏大期,80年代和21世紀(jì)初處于偏少期,這與降水的變化基本一致。

    圖3瀏陽河流域1969—2013年徑流量距平變化

    Fig.3Variationofannualrunoffanomaliesfrom1969to2013inLiuyangRiverwatershed

    3.3 降水與徑流的關(guān)系

    對(duì)1969—2013年的瀏陽河流域面降水量和下游郎梨站徑流數(shù)據(jù)作標(biāo)準(zhǔn)化處理后,得出降水和徑流標(biāo)準(zhǔn)化序列后的累積距平曲線(圖4)。徑流與同期降水的演變趨勢(shì)基本一致,徑流處于豐、枯水段時(shí)降水也同樣處于偏多、偏少期。具體來看,1969—1976年,徑流處于偏豐階段;1977—1988年,徑流處于偏枯時(shí)段;1988—1999年,徑流處于總體偏豐階段。這種階段性變化與圖3 的階段性變化基本一致,與降水的年際變化相吻合,呈現(xiàn)一致的下降—上升—下降趨勢(shì)。從圖 4可知,1988—1992年,降水和徑流發(fā)生一次顯著性波動(dòng)后,迅速轉(zhuǎn)為90年代的明顯上升期。徑流與降水的相關(guān)系數(shù)為0.9,通過了0.01的顯著性檢驗(yàn),說明降水的變化是影響徑流變化的一個(gè)重要因素。

    圖4瀏陽河流域1969—2013年標(biāo)準(zhǔn)化的降水與徑流累積距平曲線

    Fig.4Cumulativeanomalycurvesofstandardizedrunoffandprecipitationfrom1969to2013inLiuyangRiverwatershed

    根據(jù)1969—2013年間徑流和降水標(biāo)準(zhǔn)化序列的累積距平的年際變化,在1988年前后徑流量和降水的累積距平分別表現(xiàn)為減少和增大的趨勢(shì),顯然1988年為該時(shí)期的突變年份。從1988—2013年間的徑流量和降水的累積距平的年際變化可以發(fā)現(xiàn),在1999年前后徑流和降水也表現(xiàn)出先增大后減少的變化特征,其突變年份為1999年。利用滑動(dòng)t檢驗(yàn)法對(duì)徑流序列的1988年和1999年作進(jìn)一步檢測(cè),分別取前后長度均為10的2個(gè)子序列進(jìn)行計(jì)算,得到統(tǒng)計(jì)量t分別為-3和2.9,均通過t0.05=2.23的顯著性檢驗(yàn)。

    土地利用覆被變化是影響流域產(chǎn)匯流過程重要因素,從而影響降水徑流之間的關(guān)系。以1988年和1999年為分界劃分為3個(gè)階段,分別對(duì)降水和徑流的標(biāo)準(zhǔn)化數(shù)值進(jìn)行線性擬合(圖5)。在瀏陽河流域,1969—1988年和1989—1998年2個(gè)時(shí)段的降雨與徑流曲線擬合出現(xiàn)明顯的交叉現(xiàn)象,1989—1998年與1969—1988年相比,擬合線斜率增加,這反映流域地表蓄水能力相對(duì)增加,對(duì)水資源涵養(yǎng)能力有所增加,植被覆蓋情況相對(duì)較好。1989—1998年和1999—2013年2個(gè)時(shí)段的降雨與徑流曲線擬合接近平行,后時(shí)段與前時(shí)段相比,等量的降水量所形成的徑流量有所增加,這反映流域的地表蓄水能力相對(duì)減弱了,土地覆被中的不透水面積增加,從而導(dǎo)致水資源涵養(yǎng)能力降低,產(chǎn)流增加,匯流時(shí)間提前。

    圖5各時(shí)期標(biāo)準(zhǔn)化數(shù)值的徑流與降水的關(guān)系

    Fig.5Relationshipofstandardizedrunoffvs.standardizedprecipitationineachperiod

    4 降水和人類活動(dòng)對(duì)徑流影響的定量研究

    由徑流量的2個(gè)突變點(diǎn),即1988年和1999年可以將瀏陽河流域徑流量的變化劃分為3個(gè)時(shí)期:Ra時(shí)期(1969—1988年)、Rb時(shí)期(1989—1998年)和Rc時(shí)期(1999—2013年)。瀏陽河流域在1988年以前,人類活動(dòng)相對(duì)較少,變點(diǎn)1988年以前徑流量的變化主要受到氣候尤其是降水量變化的影響。1988年之后,人類活動(dòng)逐漸加強(qiáng),徑流量的變化既受降水等自然因素的變化影響,還疊加了人類活動(dòng)的影響。如圖6所示,用累積徑流量(圖6(a))和累積降水量(圖6(b))分別與年份進(jìn)行擬合。

    (a) 累積徑流量與年份的關(guān)系

    (b) 累積降水量與年份的關(guān)系

    圖6瀏陽河流域累積徑流量、累積降水量與年份的關(guān)系

    Fig.6Relationshipsofcumulativerunoffandcumulativeprecipitationvs.yearinLiuyangRiverwatershed

    Rb時(shí)期與Ra時(shí)期相比,累積徑流量與年份的線性關(guān)系式的斜率增加了8.4億m3/a,增加率為28.29%(表3)。與同時(shí)期的降水量相比,累積降水量與年份的線性關(guān)系式的斜率增加了244.9 mm/a,增加率為16.07%(表4)。根據(jù)式(3)和式(4)計(jì)算結(jié)果可知,Rb時(shí)期與Ra時(shí)期相比,降水量的變化對(duì)徑流變化的貢獻(xiàn)率為56.8%,在不考慮氣溫變化對(duì)徑流變化的影響時(shí),人類活動(dòng)對(duì)徑流量增加的貢獻(xiàn)率為43.2%。

    表3 瀏陽河流域累積徑流量斜率及其變化率

    Rc時(shí)期與Ra時(shí)期相比,累積徑流量與年份的線性關(guān)系式的斜率減少了2.69億m3/a,減少率為9.06%(表3)。與同時(shí)期的降水量相比,累積降水量與年份的線性關(guān)系式的斜率減少了10.5 mm/a,減少率為0.69%(表4)。根據(jù)式(3)和式(4)計(jì)算結(jié)果可知,Rc時(shí)期與Ra時(shí)期相比,降水量的變化對(duì)徑流變化的貢獻(xiàn)率為7.62%,在不考慮氣溫變化對(duì)徑流變化的影響時(shí),人類活動(dòng)對(duì)徑流量增加的貢獻(xiàn)率為92.38%。

    表4 瀏陽河流域累積降水量斜率及其變化率

    Rc時(shí)期與Rb時(shí)期相比,累積徑流量與年份的線性關(guān)系式的斜率減少了11.09 億m3/a,減少率變?yōu)?9.12%(表3)。與同時(shí)期的降水量相比,累積降水量與年份的線性關(guān)系式的斜率減少了255.4 mm/a,減少率為14.44%(表4)。根據(jù)式(3)和式(4)計(jì)算結(jié)果可知,Rc時(shí)期與Rb時(shí)期相比,降水量的變化對(duì)徑流變化的貢獻(xiàn)率為49.59%,在不考慮氣溫變化對(duì)徑流變化的影響時(shí),人類活動(dòng)對(duì)徑流量增加的貢獻(xiàn)率為50.41%。

    20世紀(jì)80年代以來,水利工程的不斷建設(shè),瀏陽河流域中上游興建大中型水庫2座,中型水庫5座,小型水庫225座。株樹橋引水工程的修建,將提供長沙市每日95萬t的用水量。另外,隨著瀏陽市社會(huì)經(jīng)濟(jì)農(nóng)業(yè)的高速發(fā)展以及人口的增長,工農(nóng)業(yè)和畜牧業(yè)的用水量也隨之增加,會(huì)導(dǎo)致河川徑流減少。1998年以來,受國家政策的影響,瀏陽河流域展開了退耕還林工程,林地的增加導(dǎo)致產(chǎn)流隨之減少??梢娊邓孔兓瘜?duì)徑流量變化的影響程度急劇減少,而劇烈的人類活動(dòng)已成為影響瀏陽河流域徑流變化的主要因素。

    5 結(jié) 論

    采用統(tǒng)計(jì)分析方法分析了瀏陽河流域面降水量與郎梨水文站徑流量的演變規(guī)律及相關(guān)關(guān)系,定量分析了降水變化和人類活動(dòng)對(duì)徑流量變化的貢獻(xiàn)率。主要結(jié)論如下:

    (1)1969—2013年瀏陽河流域降水量表現(xiàn)為增加趨勢(shì),上升率為3 mm/(10 a),而徑流量表現(xiàn)為下降趨勢(shì),下降率為2.7 (m3·s-1)/(10 a)。氣候變化和人類活動(dòng)的共同作用導(dǎo)致瀏陽河流域年徑流量呈下降趨勢(shì),降水量呈增加趨勢(shì)而徑流量呈下降趨勢(shì)從側(cè)面體現(xiàn)了人類活動(dòng)對(duì)徑流的影響大于蒸散發(fā)等其他氣象因素。降水量和徑流量均表現(xiàn)為相似的年際、年代際的變化特征,在20世紀(jì)70年代、90年代和21世紀(jì)初年際變化顯著,變差系數(shù)分別在0.15~0.23和0.25~0.33之間,而年代際之間均表現(xiàn)為“偏多—偏少—偏多”的交替變化。

    (2)從季節(jié)的分布來看,降水量和徑流量的季節(jié)最大值均出現(xiàn)在春季,夏季次之,秋季和冬季相對(duì)較少。春季和夏季的降水量和徑流量均占全年的60%~70%。而秋季的降水量和夏季的徑流波動(dòng)最大,變差系數(shù)分別達(dá)到了0.44和0.53。

    (3)降水量和徑流量序列均在1988和1999年發(fā)生顯著性突變。根據(jù)降水和徑流的標(biāo)準(zhǔn)化數(shù)值分階段線性擬合可知,相對(duì)于基準(zhǔn)期1969—1988年,該流域的水資源涵養(yǎng)能力和植被覆被經(jīng)歷了先增加后下降的2個(gè)階段。

    (4)以1969—1988年為基準(zhǔn)期,在不考慮蒸散發(fā)影響的條件下,定量估算了1989—1998年和1999—2013年降水量的變化和人類活動(dòng)對(duì)徑流量變化的貢獻(xiàn)率。在不同的時(shí)空尺度下,2個(gè)時(shí)期的降水量變化對(duì)徑流量變化貢獻(xiàn)率依次為56.8%和7.62%,而人類活動(dòng)對(duì)徑流量的貢獻(xiàn)率依次為43.2%和92.38%。

    [1] 宋曉猛,張建云,占車生,等.氣候變化和人類活動(dòng)對(duì)水文循環(huán)影響研究進(jìn)展[J].水利學(xué)報(bào),2013, 44(7):779-790.

    [2] 張利平,于松延,段堯彬,等.氣候變化和人類活動(dòng)對(duì)永定河流域徑流變化影響定量研究[J].氣候變化研究進(jìn)展,2013,9(6):391-397.

    [3] 董磊華,熊立華,于坤霞,等.氣候變化與人類活動(dòng)對(duì)水文影響的研究進(jìn)展[J].水科學(xué)進(jìn)展,2012, 23(2): 278-285.

    [4] 褚茜茜,陳 進(jìn),陳廣才.長江源區(qū)1978—2009年徑流極值序列變化趨勢(shì)分析[J].長江科學(xué)院院報(bào),2014,31(2):16-19.

    [5] 肖 義,唐少華,陳 華,等.湘江流域1960—2008年降水氣溫變化趨勢(shì)分析[J].人民長江,2013, 44 (3):10-12.

    [6] 羅 巧,王克林,王勤學(xué),等.基于SWAT模型的湘江流域土地利用變化情景的徑流模擬研究[J].中國生態(tài)農(nóng)業(yè)學(xué)報(bào),2011,19(6):1431-1436.

    [7] 李景保,吳桂生,劉曉清.近數(shù)十年湘江流域河流水文變化規(guī)律[J].熱帶地理,1997,17(3):295-302.

    [8] 杜 鵑,徐 偉,趙智國,等.湘江流域近30年徑流量與水位的長期變化規(guī)律研究[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào),2010,18(3):369-378.

    [9] 張海斌.瀏陽河流域水資源分析與評(píng)價(jià)[J].人民長江,2011,42 (增2):32-35.

    [10] 張 棋,李忠武,曾光明,等.瀏陽河土地利用變化對(duì)非點(diǎn)源污染負(fù)荷的影響[J].環(huán)境工程學(xué)報(bào),2009,3(2):377-380.

    [11] 王隨繼,閆云霞,顏 明,等.皇甫川流域降水和人類活動(dòng)對(duì)徑流量變化的貢獻(xiàn)率分析——累積量斜率變化率比較方法的提出及應(yīng)用[J].地理學(xué)報(bào),2012,67(3):388-397.

    [12] SAID S E, DICKEY D A.Testing for Unit Roots in Autoregressive-moving Average Models of Unknown Order[J]. Biometrika,1984,71(3):599-607.

    [13] WANG S,YAN Y,YAN M,etal.Quantitative Estimation of the Impact of Precipitation and Human Activities on Runoff Change of the Huangfuchuan River Basin[J].Journal of Geographical Sciences,2012, 22(5):906-918.

    [14] 許正福,郭連云.共和盆地貴南地區(qū)近50年降水特征及變化規(guī)律[J].中國農(nóng)學(xué)通報(bào),2016,32(15):144-149.

    [15] 穆興民,張秀勤,高 鵬,等.雙累積曲線方法理論及在水文氣象領(lǐng)域應(yīng)用中應(yīng)注意的問題[J].水文,2010,30(4):47-51.

    [16] 張婷婷.湘江流域蒸發(fā)皿蒸發(fā)量的變化趨勢(shì)及原因分析[D].長沙:湖南師范大學(xué),2013:21-28.

    [17] GB/T 50095—98,水文基本術(shù)語和符號(hào)標(biāo)準(zhǔn)[S]. 北京:中國計(jì)劃出版社,1998.

    (編輯:羅 娟)

    Effects of Precipitation Variation and Human Activities on Runoff in Liuyang River Watershed

    YU Ze-xing1,HU Guo-hua1,CHEN Xiao1,GU Qing-fu2

    (1.School of Hydraulic Engineering, Changsha University of Science amp; Technology, Changsha 410114,China; 2.Hydrology amp; Water Resources Bureau of Hunan Province, Changsha 410007,China)

    Quantitative evaluation of the effects of precipitation variation and human activities on runoff in Liuyang River watershed is conducive to a more rational utilization of water resource. In this paper we analyzed the evolution law of and the relationship between runoff and precipitation in 1969—2013 in Liuyang River basin by means of linear fitting, moving average and accumulative anomaly. We also estimated the contribution rates of precipitation variation and human activities to runoff variation in different periods quantitatively by comparing cumulative slope variance ratio. Results indicated that precipitation displayed an increasing trend in a rate of 3 mm/(10 a), while runoff represented a decreasing trend in a rate of 2.7 (m3·s-1)/(10 a). Precipitation and runoff both witnessed an abrupt change in 1989 and 1999 in a trend of “decrease-increase-decrease”. Compared with the baseline period from 1969 to 1988, the contribution rates of precipitation variation to runoff variation in 1989—1998 and 1999—2013 were 56.8% and 7.62%, while the contribution rates of human activities in 1989—1998 and 1999—2013 were 43.2% and 92.38%, indicating that human activity has become the main factor which affects the runoff variation in Liuyang River watershed.

    Liuyang River watershed; runoff variation; method of comparing cumulative slope variance ratio; precipitation variation; human activities

    2016-07-18;

    2016-08-25

    湖南省重大水利科技項(xiàng)目(2013-243-2);湖南省研究生科研創(chuàng)新項(xiàng)目(CX2015B359)

    于澤興(1991-),男,湖南岳陽人,碩士研究生,從事水文預(yù)報(bào)研究,(電話)14789910820(電子信箱)423552622@qq.com。

    胡國華(1965-),男,湖南新田人,教授,博士,從事水資源演變規(guī)律研究,(電話)13974831708(電子信箱)g.h.hu@163.com。

    10.11988/ckyyb.20160735

    TV21

    A

    1001-5485(2017)11-0006-06

    猜你喜歡
    人類
    顫抖吧,人類
    人類能否一覺到未來?
    人類會(huì)成長起來嗎?
    英語文摘(2021年2期)2021-07-22 07:56:54
    人類沒有進(jìn)化好
    人類過去長啥樣
    出發(fā)!去最北人類定居地
    人類第一殺手
    共建“人類之眼”
    人類,快醒醒
    1100億個(gè)人類的清明
    中文字幕人成人乱码亚洲影| 国产爱豆传媒在线观看| 婷婷精品国产亚洲av| 亚洲va日本ⅴa欧美va伊人久久| 欧美日韩瑟瑟在线播放| 日韩成人在线观看一区二区三区| 嫩草影院新地址| 在线看三级毛片| 九色国产91popny在线| 亚洲经典国产精华液单 | 国产精品亚洲av一区麻豆| 国产探花在线观看一区二区| 日韩欧美 国产精品| 国产蜜桃级精品一区二区三区| 在线观看美女被高潮喷水网站 | 少妇的逼水好多| 永久网站在线| 国产成人福利小说| 欧美xxxx性猛交bbbb| 99久久成人亚洲精品观看| 少妇熟女aⅴ在线视频| av在线观看视频网站免费| 床上黄色一级片| 日日摸夜夜添夜夜添小说| 能在线免费观看的黄片| 观看免费一级毛片| 国产乱人伦免费视频| 一级a爱片免费观看的视频| 噜噜噜噜噜久久久久久91| 能在线免费观看的黄片| 真人做人爱边吃奶动态| 日韩免费av在线播放| 五月玫瑰六月丁香| 久久久色成人| 中文字幕高清在线视频| 日本五十路高清| 精品久久久久久久久av| 午夜激情福利司机影院| 真人做人爱边吃奶动态| 2021天堂中文幕一二区在线观| 久久久久久久午夜电影| 在线免费观看不下载黄p国产 | 亚洲av二区三区四区| 99视频精品全部免费 在线| 亚洲 欧美 日韩 在线 免费| 精品人妻偷拍中文字幕| 亚洲乱码一区二区免费版| 亚洲综合色惰| 国产一区二区在线av高清观看| 乱人视频在线观看| 老鸭窝网址在线观看| 可以在线观看的亚洲视频| 麻豆av噜噜一区二区三区| 国产精品伦人一区二区| 国产在视频线在精品| 国产精品不卡视频一区二区 | 色综合站精品国产| 女人十人毛片免费观看3o分钟| 亚洲成av人片免费观看| 日本黄色片子视频| 黄色丝袜av网址大全| 免费一级毛片在线播放高清视频| 欧美激情国产日韩精品一区| 女同久久另类99精品国产91| 波多野结衣高清无吗| 久久婷婷人人爽人人干人人爱| 亚洲精品粉嫩美女一区| 69人妻影院| 日韩亚洲欧美综合| 久久精品国产99精品国产亚洲性色| 精品人妻一区二区三区麻豆 | 天堂动漫精品| 亚洲欧美日韩卡通动漫| 久9热在线精品视频| 如何舔出高潮| 国产精品一区二区免费欧美| 精品国内亚洲2022精品成人| 亚洲av中文字字幕乱码综合| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 久久久久国产精品人妻aⅴ院| 婷婷精品国产亚洲av| 国内少妇人妻偷人精品xxx网站| 成人三级黄色视频| 88av欧美| 国产av不卡久久| 亚洲七黄色美女视频| 又紧又爽又黄一区二区| 中文资源天堂在线| 啦啦啦观看免费观看视频高清| 免费av观看视频| 欧美日本视频| 三级男女做爰猛烈吃奶摸视频| 国产极品精品免费视频能看的| 亚洲精品亚洲一区二区| 亚洲精品成人久久久久久| 久久精品国产亚洲av香蕉五月| 午夜老司机福利剧场| av在线老鸭窝| 亚洲中文字幕一区二区三区有码在线看| 日本精品一区二区三区蜜桃| 热99re8久久精品国产| 人人妻人人看人人澡| av天堂在线播放| 欧美性感艳星| 亚洲第一电影网av| 国产爱豆传媒在线观看| 日日干狠狠操夜夜爽| 欧美另类亚洲清纯唯美| 亚洲精品一卡2卡三卡4卡5卡| 最新中文字幕久久久久| 在现免费观看毛片| 亚洲无线观看免费| av欧美777| 性色avwww在线观看| 日韩精品中文字幕看吧| 乱码一卡2卡4卡精品| 91九色精品人成在线观看| 国产午夜福利久久久久久| 18禁黄网站禁片免费观看直播| 97人妻精品一区二区三区麻豆| 村上凉子中文字幕在线| 日韩欧美三级三区| 午夜亚洲福利在线播放| 午夜福利视频1000在线观看| 免费看日本二区| 国产 一区 欧美 日韩| 赤兔流量卡办理| 全区人妻精品视频| 精品久久久久久久久亚洲 | 黄色日韩在线| 久久久久久九九精品二区国产| 99久久久亚洲精品蜜臀av| av在线观看视频网站免费| 美女免费视频网站| 男人舔女人下体高潮全视频| 欧美性猛交╳xxx乱大交人| 色播亚洲综合网| 乱码一卡2卡4卡精品| 亚洲成人久久性| 九色成人免费人妻av| 欧美在线一区亚洲| 日韩欧美免费精品| 午夜福利免费观看在线| 在线观看美女被高潮喷水网站 | 成人av在线播放网站| 97超级碰碰碰精品色视频在线观看| 91麻豆av在线| 人人妻人人看人人澡| 小说图片视频综合网站| 此物有八面人人有两片| 日韩欧美三级三区| 精品免费久久久久久久清纯| 无人区码免费观看不卡| 中文字幕人妻熟人妻熟丝袜美| 少妇熟女aⅴ在线视频| 男人舔奶头视频| 黄片小视频在线播放| 亚洲七黄色美女视频| 无人区码免费观看不卡| 人妻丰满熟妇av一区二区三区| 啦啦啦观看免费观看视频高清| 一个人看视频在线观看www免费| 精品国产亚洲在线| 免费av不卡在线播放| 精品一区二区三区人妻视频| 久久亚洲精品不卡| 国产精品女同一区二区软件 | 免费在线观看亚洲国产| 免费人成在线观看视频色| 亚洲成人免费电影在线观看| eeuss影院久久| 极品教师在线免费播放| 熟女人妻精品中文字幕| 十八禁国产超污无遮挡网站| 老熟妇乱子伦视频在线观看| 亚洲美女黄片视频| 欧美一区二区国产精品久久精品| 国产亚洲精品久久久com| 国产av一区在线观看免费| 日日摸夜夜添夜夜添av毛片 | 欧美在线黄色| 欧美一区二区亚洲| 国产中年淑女户外野战色| 99精品在免费线老司机午夜| 久久精品91蜜桃| 他把我摸到了高潮在线观看| 国产精品女同一区二区软件 | 亚洲久久久久久中文字幕| 国产成人freesex在线| 欧美精品人与动牲交sv欧美| 蜜臀久久99精品久久宅男| 又粗又硬又长又爽又黄的视频| 精品久久久精品久久久| 免费黄网站久久成人精品| 少妇人妻一区二区三区视频| 欧美+日韩+精品| 看黄色毛片网站| 国产乱人视频| 国产黄a三级三级三级人| 午夜福利视频1000在线观看| av在线app专区| av天堂中文字幕网| 精品国产乱码久久久久久小说| 嫩草影院新地址| 女人被狂操c到高潮| 亚洲精华国产精华液的使用体验| 麻豆久久精品国产亚洲av| 国产精品国产三级国产av玫瑰| 久久午夜福利片| 熟妇人妻不卡中文字幕| 女人久久www免费人成看片| 神马国产精品三级电影在线观看| 亚洲自偷自拍三级| 日本一二三区视频观看| 亚洲国产精品国产精品| h日本视频在线播放| freevideosex欧美| 综合色av麻豆| 少妇人妻精品综合一区二区| 国产伦精品一区二区三区四那| 国产免费福利视频在线观看| 色吧在线观看| 97热精品久久久久久| 又爽又黄a免费视频| 在线a可以看的网站| 亚洲怡红院男人天堂| 国产精品久久久久久av不卡| 狠狠精品人妻久久久久久综合| 国产毛片a区久久久久| 日本色播在线视频| 日韩成人av中文字幕在线观看| 插阴视频在线观看视频| 亚洲av在线观看美女高潮| 夫妻性生交免费视频一级片| 日韩制服骚丝袜av| 亚洲高清免费不卡视频| 亚洲国产精品国产精品| 人人妻人人澡人人爽人人夜夜| 26uuu在线亚洲综合色| 亚洲欧美日韩卡通动漫| 99热这里只有精品一区| 亚洲成色77777| 婷婷色综合大香蕉| 激情 狠狠 欧美| 国产精品嫩草影院av在线观看| 欧美日韩亚洲高清精品| 18+在线观看网站| 亚洲天堂av无毛| 国产免费视频播放在线视频| 久热久热在线精品观看| 色综合色国产| 国产伦精品一区二区三区视频9| 午夜老司机福利剧场| 国产色爽女视频免费观看| 精品熟女少妇av免费看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 18禁在线播放成人免费| 一级黄片播放器| 特级一级黄色大片| 少妇被粗大猛烈的视频| 嫩草影院入口| 亚洲天堂国产精品一区在线| 一个人看视频在线观看www免费| 久久久午夜欧美精品| 免费大片18禁| 狠狠精品人妻久久久久久综合| av福利片在线观看| 日韩三级伦理在线观看| 男人爽女人下面视频在线观看| 大又大粗又爽又黄少妇毛片口| 永久免费av网站大全| 91久久精品国产一区二区三区| 性色av一级| 伦理电影大哥的女人| 日韩大片免费观看网站| 一区二区三区乱码不卡18| 国产成人a区在线观看| 精品熟女少妇av免费看| 色哟哟·www| 干丝袜人妻中文字幕| 国产伦理片在线播放av一区| 欧美日韩视频精品一区| 亚洲aⅴ乱码一区二区在线播放| 伊人久久国产一区二区| 亚洲第一区二区三区不卡| 2018国产大陆天天弄谢| 成人无遮挡网站| 精品亚洲乱码少妇综合久久| 亚洲经典国产精华液单| 免费av毛片视频| 久久久久久久午夜电影| 国产精品成人在线| 一级毛片久久久久久久久女| 久热久热在线精品观看| 成年版毛片免费区| 九九久久精品国产亚洲av麻豆| 春色校园在线视频观看| 国产欧美日韩精品一区二区| 亚洲av中文av极速乱| 亚洲,一卡二卡三卡| 新久久久久国产一级毛片| 大话2 男鬼变身卡| 嫩草影院精品99| 中文字幕亚洲精品专区| 大陆偷拍与自拍| 亚洲内射少妇av| 性色av一级| 波野结衣二区三区在线| 爱豆传媒免费全集在线观看| 欧美xxxx黑人xx丫x性爽| 国产免费又黄又爽又色| 精品午夜福利在线看| 久久精品国产鲁丝片午夜精品| 欧美日韩视频精品一区| 亚洲精品456在线播放app| 男女啪啪激烈高潮av片| 亚洲欧美中文字幕日韩二区| 91久久精品电影网| 国产 精品1| 大话2 男鬼变身卡| 一区二区三区免费毛片| 国产欧美亚洲国产| 69av精品久久久久久| 在线精品无人区一区二区三 | 国产精品一区www在线观看| 精品国产露脸久久av麻豆| 国产成人a∨麻豆精品| 午夜免费观看性视频| 丰满人妻一区二区三区视频av| 亚洲不卡免费看| 成人黄色视频免费在线看| 久久久久性生活片| freevideosex欧美| 国产精品一二三区在线看| 国内少妇人妻偷人精品xxx网站| 男人和女人高潮做爰伦理| 国产有黄有色有爽视频| 禁无遮挡网站| 亚洲av日韩在线播放| 久久国内精品自在自线图片| 日本-黄色视频高清免费观看| 能在线免费看毛片的网站| 波野结衣二区三区在线| 欧美日韩国产mv在线观看视频 | 日韩制服骚丝袜av| av黄色大香蕉| 免费不卡的大黄色大毛片视频在线观看| 久久精品国产亚洲av天美| 日韩精品有码人妻一区| 久久精品国产鲁丝片午夜精品| 久久久久精品久久久久真实原创| 亚洲图色成人| 日日撸夜夜添| 国产日韩欧美在线精品| 久久久国产一区二区| 国产高清三级在线| 大片电影免费在线观看免费| 亚洲av中文字字幕乱码综合| 国产高清国产精品国产三级 | 成年版毛片免费区| 国国产精品蜜臀av免费| 久久精品久久久久久久性| 伦理电影大哥的女人| 色哟哟·www| 欧美成人午夜免费资源| 交换朋友夫妻互换小说| 男男h啪啪无遮挡| 美女被艹到高潮喷水动态| 黄片wwwwww| 水蜜桃什么品种好| 国产黄片美女视频| 色哟哟·www| 国产91av在线免费观看| 两个人的视频大全免费| 亚洲,一卡二卡三卡| 秋霞伦理黄片| 日韩伦理黄色片| 人妻少妇偷人精品九色| 日韩成人伦理影院| 亚洲一级一片aⅴ在线观看| 超碰97精品在线观看| 寂寞人妻少妇视频99o| 亚洲,一卡二卡三卡| 可以在线观看毛片的网站| 亚洲成人精品中文字幕电影| 欧美xxxx黑人xx丫x性爽| 亚洲精品一二三| 国产精品一及| 国产69精品久久久久777片| 老司机影院毛片| 男女国产视频网站| 丝瓜视频免费看黄片| 国产黄色免费在线视频| 一本久久精品| 精品少妇黑人巨大在线播放| 啦啦啦在线观看免费高清www| 国产亚洲一区二区精品| 制服丝袜香蕉在线| 最近中文字幕高清免费大全6| 国产色婷婷99| 精品国产三级普通话版| 亚洲精品第二区| 国产欧美亚洲国产| 2018国产大陆天天弄谢| 乱系列少妇在线播放| 大话2 男鬼变身卡| 午夜福利视频1000在线观看| 国产精品一区二区三区四区免费观看| 国产黄a三级三级三级人| 久久精品国产亚洲av天美| 只有这里有精品99| 五月天丁香电影| 久久久久久久久久久免费av| 中国国产av一级| 国产精品麻豆人妻色哟哟久久| 97人妻精品一区二区三区麻豆| 久久久久精品性色| 欧美日韩综合久久久久久| 久久综合国产亚洲精品| 大又大粗又爽又黄少妇毛片口| 一区二区av电影网| 国产69精品久久久久777片| av在线蜜桃| 亚洲婷婷狠狠爱综合网| 特级一级黄色大片| 久久精品综合一区二区三区| 在线观看人妻少妇| 成人国产麻豆网| 天堂俺去俺来也www色官网| .国产精品久久| 黄色怎么调成土黄色| 日韩一区二区视频免费看| 色吧在线观看| 成年av动漫网址| 亚洲成人精品中文字幕电影| 国产精品99久久99久久久不卡 | 伦精品一区二区三区| 麻豆国产97在线/欧美| 国产乱人偷精品视频| 精品久久久久久久久亚洲| 国产乱来视频区| 亚洲欧美一区二区三区黑人 | 国产精品久久久久久精品电影小说 | 亚洲伊人久久精品综合| 久久人人爽人人片av| av在线蜜桃| 国产在视频线精品| 九草在线视频观看| 99九九线精品视频在线观看视频| 可以在线观看毛片的网站| 成人国产av品久久久| 亚洲欧美精品专区久久| 下体分泌物呈黄色| 国产熟女欧美一区二区| 免费观看性生交大片5| 国产午夜精品一二区理论片| 国产欧美另类精品又又久久亚洲欧美| 久久久久久久大尺度免费视频| 又大又黄又爽视频免费| 在线观看人妻少妇| 精品一区在线观看国产| 国产一区二区三区综合在线观看 | 国产色爽女视频免费观看| 国产黄a三级三级三级人| 亚洲无线观看免费| 久久人人爽av亚洲精品天堂 | 亚洲国产av新网站| 久久精品熟女亚洲av麻豆精品| 欧美性感艳星| 日本一本二区三区精品| 五月玫瑰六月丁香| 中文资源天堂在线| 欧美xxxx黑人xx丫x性爽| 我的老师免费观看完整版| 99热这里只有是精品在线观看| av播播在线观看一区| 国产久久久一区二区三区| 国产精品国产三级国产av玫瑰| 丝袜喷水一区| 欧美成人一区二区免费高清观看| 嘟嘟电影网在线观看| 日本色播在线视频| 97在线人人人人妻| 亚洲国产欧美在线一区| 麻豆成人av视频| 18禁在线无遮挡免费观看视频| 肉色欧美久久久久久久蜜桃 | 国产精品偷伦视频观看了| 成年av动漫网址| 国产成人精品婷婷| 人人妻人人爽人人添夜夜欢视频 | 久久久成人免费电影| 成人一区二区视频在线观看| 成人漫画全彩无遮挡| 成人亚洲精品一区在线观看 | 国产v大片淫在线免费观看| 色吧在线观看| 国产精品不卡视频一区二区| 97热精品久久久久久| 中文字幕av成人在线电影| 国产精品久久久久久久久免| 爱豆传媒免费全集在线观看| 色5月婷婷丁香| 欧美成人a在线观看| 男的添女的下面高潮视频| 美女视频免费永久观看网站| 亚洲精品色激情综合| av在线亚洲专区| 亚洲高清免费不卡视频| 亚洲精品视频女| 欧美极品一区二区三区四区| 日韩av免费高清视频| 天天一区二区日本电影三级| 在现免费观看毛片| 哪个播放器可以免费观看大片| 尤物成人国产欧美一区二区三区| 国国产精品蜜臀av免费| 国产精品一区二区三区四区免费观看| 亚洲精品影视一区二区三区av| 在线a可以看的网站| 美女cb高潮喷水在线观看| 亚洲成人av在线免费| 亚洲,一卡二卡三卡| 亚洲色图av天堂| 久久99热这里只频精品6学生| 中文字幕久久专区| 国产成人91sexporn| 极品少妇高潮喷水抽搐| 97超视频在线观看视频| 国产成人a区在线观看| 国产精品av视频在线免费观看| 51国产日韩欧美| 国内精品美女久久久久久| 99久久精品国产国产毛片| 国产成人精品福利久久| 国产精品爽爽va在线观看网站| 国产成人午夜福利电影在线观看| 我要看日韩黄色一级片| 国产伦在线观看视频一区| 一级毛片久久久久久久久女| 十八禁网站网址无遮挡 | 91午夜精品亚洲一区二区三区| 国产精品久久久久久精品电影小说 | 亚洲人成网站在线观看播放| 日韩强制内射视频| 欧美人与善性xxx| 91在线精品国自产拍蜜月| 伦理电影大哥的女人| 精品久久久久久久末码| 老司机影院毛片| 99久久精品国产国产毛片| 99视频精品全部免费 在线| 美女cb高潮喷水在线观看| 日本黄大片高清| 午夜爱爱视频在线播放| 蜜臀久久99精品久久宅男| 成人二区视频| 成人黄色视频免费在线看| 欧美成人精品欧美一级黄| 日本wwww免费看| 国产综合懂色| 国产成人免费观看mmmm| 乱码一卡2卡4卡精品| 菩萨蛮人人尽说江南好唐韦庄| 男女那种视频在线观看| 日本猛色少妇xxxxx猛交久久| 久久这里有精品视频免费| 久久精品综合一区二区三区| 大又大粗又爽又黄少妇毛片口| 亚洲精品色激情综合| 午夜精品一区二区三区免费看| 欧美精品国产亚洲| 欧美变态另类bdsm刘玥| 国产高清国产精品国产三级 | 日日摸夜夜添夜夜添av毛片| 久久精品久久精品一区二区三区| 国产白丝娇喘喷水9色精品| 国模一区二区三区四区视频| 亚洲欧美一区二区三区黑人 | 欧美最新免费一区二区三区| 亚洲精品国产av成人精品| 麻豆久久精品国产亚洲av| 久久精品夜色国产| 亚洲精品日韩av片在线观看| 日本三级黄在线观看| 黑人高潮一二区| 少妇的逼水好多| 少妇裸体淫交视频免费看高清| 五月开心婷婷网| 免费观看的影片在线观看| 1000部很黄的大片| 人妻系列 视频| 26uuu在线亚洲综合色| 观看免费一级毛片| 精品久久国产蜜桃| 91精品一卡2卡3卡4卡| 免费观看无遮挡的男女| 你懂的网址亚洲精品在线观看| 在线免费十八禁| 欧美亚洲 丝袜 人妻 在线| 在线亚洲精品国产二区图片欧美 | 国语对白做爰xxxⅹ性视频网站| 白带黄色成豆腐渣| 春色校园在线视频观看| 天天躁夜夜躁狠狠久久av| 国产精品久久久久久久久免| 国产成人91sexporn| 久久这里有精品视频免费| 亚洲精品久久午夜乱码| 九九在线视频观看精品| av天堂中文字幕网| 亚洲精品日本国产第一区| 2021天堂中文幕一二区在线观| 亚洲美女搞黄在线观看|