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

    渭河流域水沙時(shí)空變化及其對(duì)人類活動(dòng)的響應(yīng)

    2020-10-12 02:45:14徐瑞瑞高鵬穆興民柴雪柯顧朝軍
    人民黃河 2020年3期
    關(guān)鍵詞:人類活動(dòng)渭河流域輸沙量

    徐瑞瑞 高鵬 穆興民 柴雪柯 顧朝軍

    摘 要:基于渭河流域32個(gè)水文站近60 a水沙數(shù)據(jù),采用Mann-Kendall趨勢(shì)檢驗(yàn)、Pettitt突變檢驗(yàn)、雙累計(jì)曲線等方法分析了渭河流域不同區(qū)間水沙變化特征及驅(qū)動(dòng)因素。結(jié)果表明:渭河流域62.5%(50%)區(qū)間徑流量(輸沙量)呈顯著減少的趨勢(shì),34.36%(42.86%)區(qū)間徑流量(輸沙量)無顯著變化趨勢(shì),僅3.12%的區(qū)間徑流量和輸沙量呈顯著增加趨勢(shì);渭河徑流量和輸沙量的突變年份集中在20世紀(jì)80年代末、90年代初;涇河徑流量和輸沙量的突變年份集中在20世紀(jì)90年代;北洛河徑流量和輸沙量在21世紀(jì)初發(fā)生突變;以渭河控制站(頭站和華縣站)水沙之和為基礎(chǔ)量化的人類活動(dòng)和降雨量對(duì)渭河徑流量減少的貢獻(xiàn)率分別為84.86%和15.14%,對(duì)輸沙量減少的貢獻(xiàn)率為92.54%和7.46%;退耕還林(草)工程、淤地壩建設(shè)、水庫攔蓄等是渭河徑流量和輸沙量變化的主要原因。

    關(guān)鍵詞:渭河流域;徑流量;輸沙量;水土保持;人類活動(dòng)

    中圖分類號(hào):P333?? 文獻(xiàn)標(biāo)志碼:A

    doi:10.3969/j.issn.1000-1379.2020.03.003

    Dynamic of Streamflow and Sediment Load and Its Response to Human Activities in the Weihe River Basin

    XU Ruirui1, GAO Peng1,2, MU Xingmin1,2, CHAI Xueke3, GU Chaojun4

    (1.Institute of Soil and Water Conservation, Northwest A & F University, Yangling 712100, China; 2.Institute of Soil

    and Water Conservation, Chinese Academy of Sciences and Ministry of Water Resources, Yangling 712100, China;

    3.Station of Management and Protection of Dongshan Forest Area, Taiyuan 030032, China; 4.Yangtze River Basin

    Soil and Water Conservation Monitoring Center Station, Changjiang Water Resources Commission, Wuhan 430010, China)

    Abstract: This study systematically analyzed the change of streamflow and sediment load and the impact of human activitiesto the Weihe River basin based on the date from 32 hydrological stations with Mann-Kendall test, double mass curve method and Pettitt test. The results show that annual streamflow and sediment load has significant decreasing trends in 62.5% and 50% areas of the Weihe River and only 3.12% areas has significant increasing trends, others without significant trend. Change points of streamflow and sediment load are concentrated between 1980s and 1990s in 25%-28% areas of the Weihe River, while 12%-16% areas of the Jinghe River mainly are appeared in the 1990s and 6.25% areas of Beiluo River are occurred in 2000s. Meanwhile, both streamflow modulus and sediment load modulus exhibit a decreasing tendency. Double mass curve reveals that human activities contribute 84.86% for the reduction of streamflow and 92.54% for the sediment load. Therefore, human activities are the dominant factors for streamflow and sediment reduction in the Weihe River basin.

    Key words: Weihe River basin; runoff; sediment discharge; soil and water conservation; human activities

    自20世紀(jì)50年代以來,受氣候變化和人類活動(dòng)的影響,河川水文情勢(shì)發(fā)生了顯著變化[1-2]。渭河作為黃河的最大支流,被陜西人譽(yù)為“母親河”,其流域水生態(tài)環(huán)境變化影響著整個(gè)區(qū)域農(nóng)業(yè)經(jīng)濟(jì)、洪旱災(zāi)害和水資源利用。近年來,受自然和人為因素的影響,渭河流域水資源分布不均,上中游水土流失嚴(yán)重,下游泥沙大量淤積,水沙問題嚴(yán)重[3],加上洪水災(zāi)害頻發(fā)(如“92·8”“96·7”“03·8”“03·10”等洪水)[4-5],嚴(yán)重影響了當(dāng)?shù)厝嗣竦纳a(chǎn)生活及自然、經(jīng)濟(jì)和社會(huì)的發(fā)展[6]。因此,關(guān)于渭河流域水沙變化的研究已成為當(dāng)前的熱點(diǎn)。

    為揭示渭河流域水沙的演變規(guī)律,我國學(xué)者開展了大量研究,如冉大川等[7]利用渭河流域1970—1996年的徑流泥沙資料,分析了流域水沙時(shí)空演變;胡安焱等[8]利用降雨、水利水保措施等資料分析了渭河流域近50 a來水沙變化規(guī)律,并探討了降水和人類活動(dòng)對(duì)流域水沙的影響;任立良等[9]采用標(biāo)準(zhǔn)化徑流指數(shù)剖析了渭河流域自1991年以來年徑流量顯著減小的影響因素;焦菊英等[10]采用“水文-地貌法”分析了渭河流域侵蝕產(chǎn)沙強(qiáng)度的區(qū)域分異特征。這些研究主要通過個(gè)別站點(diǎn)闡述了渭河流域的水沙變化特征和影響因素(降雨、氣候、水土保持措施、水庫等),主要關(guān)注徑流量和泥沙量,缺乏對(duì)整個(gè)渭河流域水沙特征的分析,難以全面深入理解渭河流域徑流泥沙變化格局及變化的主控因素。本文利用渭河流域32個(gè)水文站的資料,系統(tǒng)研究了渭河流域不同區(qū)間的水沙變化特征,以期全面認(rèn)識(shí)流域水沙變化規(guī)律及其驅(qū)動(dòng)因素,為區(qū)域水資源管理、生產(chǎn)發(fā)展、生態(tài)文明建設(shè)提供科學(xué)參考。

    1 流域概況

    渭河發(fā)源于鳥鼠山,流經(jīng)陜、甘、寧3個(gè)省區(qū),于陜西潼關(guān)附近匯入黃河[11]。流域總面積為134 766 km2,干流全長(zhǎng)為818 km,流域?qū)儆诖箨懶约撅L(fēng)氣候區(qū),降水主要集中在夏季且雨強(qiáng)較大,對(duì)地表侵蝕嚴(yán)重。渭河支流較多,涇河和洛河是北岸支流中最大的兩條,流域面積分別為4.54萬、2.69萬km2。渭河流域年均降水量為400~800 mm,年均蒸發(fā)量為1 000~2 000 mm。流域多年平均徑流量為95億m3,多年平均輸沙量為4.58億t。整個(gè)流域水沙異源,主要產(chǎn)水區(qū)為林家村、張家山和華縣區(qū)間,主要產(chǎn)沙區(qū)為涇河和林家村以上[12]。流域內(nèi)及周邊氣象站、水文站分布見圖1。

    2 數(shù)據(jù)與方法

    2.1 數(shù) 據(jù)

    選用渭河流域內(nèi)32個(gè)具有代表性的水文站點(diǎn),將流域劃分為32個(gè)不同的區(qū)間(見表1)。區(qū)間徑流輸沙量等于下游水文站實(shí)測(cè)水沙量減去上游水文站實(shí)測(cè)水沙量,某個(gè)水文站上游無水文站時(shí),該水文站實(shí)測(cè)值即為該區(qū)間水沙量。數(shù)據(jù)年限為1958—2013年,均來源于《中國河流泥沙公報(bào)》和《黃河流域水文年鑒》等。降雨資料選用渭河流域及周邊24個(gè)氣象站1958—2013年逐日和逐年的降雨資料(見表2),來源于中國國家氣象數(shù)據(jù)共享服務(wù)網(wǎng)(http://data.cma.gov.cn)。

    2.2 分析方法

    2.2.1 趨勢(shì)檢驗(yàn)

    主要采用Mann-Kendall(簡(jiǎn)稱MK)秩次相關(guān)檢驗(yàn)方法[13-14]對(duì)長(zhǎng)時(shí)間的水沙序列進(jìn)行趨勢(shì)性分析[15-16]。通過該法計(jì)算得檢驗(yàn)統(tǒng)計(jì)量Z,若Z值為正則說明序列呈現(xiàn)增大趨勢(shì),若Z為負(fù)則說明序列呈現(xiàn)減小趨勢(shì),0.05和0.01顯著性水平下的Z值分別為±1.96和±2.58。

    2.2.2 突變檢驗(yàn)

    采用Pettitt突變檢驗(yàn)法[17]來診斷水沙的突變年份,該方法是檢驗(yàn)時(shí)間序列均值發(fā)生突變的統(tǒng)計(jì)方法,既能判斷突變點(diǎn)的位置及數(shù)量,也能判斷突變點(diǎn)是否具有統(tǒng)計(jì)意義上的顯著水平,在變點(diǎn)檢測(cè)中有較多的應(yīng)用。通過該方法找到突變點(diǎn)位置,再通過突變點(diǎn)顯著性水平計(jì)算出P值,P值越接近1,突變點(diǎn)的趨勢(shì)越顯著。

    2.2.3 雙累計(jì)曲線法

    雙累計(jì)曲線能直觀、簡(jiǎn)單地分析兩個(gè)變量之間的一致性或長(zhǎng)期趨勢(shì)性變化[18-19],已被廣泛應(yīng)用于診斷流域水沙關(guān)系變異情況。該法通過兩個(gè)累計(jì)變量之間擬合直線斜率的偏離程度來判斷人類活動(dòng)對(duì)流域水沙的影響程度,若斜率偏離明顯,則說明人類活動(dòng)對(duì)流域水沙變化有著顯著的影響,斜率顯著變化的點(diǎn)對(duì)應(yīng)水沙關(guān)系發(fā)生顯著變化的時(shí)間[20]。

    3 結(jié)果與分析

    3.1 渭河流域徑流量、輸沙量變化的基本特征

    3.1.1 渭河流域徑流量、輸沙量年代際變化

    華縣站和頭站為渭河流域和北洛河的把口站,兩個(gè)站徑流量和輸沙量的變化一定程度上能代表整個(gè)流域水沙的變化。分析整個(gè)流域不同年代的水沙特征(見表3)可知:渭河流域徑流量和輸沙量年代際變化趨勢(shì)基本一致,均呈波動(dòng)下降的趨勢(shì),但輸沙量的波動(dòng)幅度明顯大于徑流量的。華縣站多年平均徑流量為64.74億m3,最大年徑流量為187.60億m3(1964年),最小年徑流量為16.83億m3(1997年);輸沙量與徑流量變化趨勢(shì)一致但變幅更大,多年平均輸沙量為2.67億t,最大年輸沙量為10.60億t (1964年),最小值為0.41億t(2012年)。頭站多年平均徑流量為7.09億m3,最大年徑流量20.15億m3 (1964年),最小年徑流量為3.09億m3 (1974年);輸沙量的變化幅度也大于徑流量的,最大年輸沙量為2.63億t(1994年),最小年輸沙量為0.01億t(2009年)。年代際內(nèi),流域的面降水量在20世紀(jì)60年代最大,為592.69 mm,90年代最小,為491.27 mm,年代際變化幅度不大。

    3.1.2 渭河流域徑流量、輸沙量年際變化

    渭河流域水文要素年際變化見圖2。20世紀(jì)60年代以來渭河流域控制站點(diǎn)(華縣、頭和張家山)的徑流量和輸沙量均呈減小趨勢(shì),華縣、頭、張家山的平均年徑流量分別為64.74億、7.09億、15.67億m3,平均年輸沙量為2.67億、0.57億、0.20億t。年徑流量呈顯著減小趨勢(shì),尤其是20世紀(jì)70年代初徑流量減少幅度較大,相比于60年代,華縣減少38.24%、張家山減少19.57%、頭站減小39.92%,這與黃河中游地區(qū)近年來降雨量的減少有直接關(guān)系。其中,涇河張家山站的徑流量減小最為顯著,從21.67億m3(1960—1969年)減小到9.28億m3(2000—2009年),減小57.17%。與徑流量變化情況基本一致,各站年輸沙量主要在20世紀(jì)70年代初呈顯著減小趨勢(shì),從60年代到90年代,各站輸沙量的減幅均在50%以上,其中,北洛河頭站平均年輸沙量的減小最顯著,由60年代的1.00億t減小到2000—2009年的0.14億t,減小86%。輸沙量的減小主要是20世紀(jì)70年代末和80年代初實(shí)施的大規(guī)模水土保持措施導(dǎo)致的,如退耕還林(草)、梯田和淤地壩等,顯著減小了坡面土壤侵蝕量,降低了輸沙能力,進(jìn)一步使得河道的輸沙量減小[21]。

    3.2 不同區(qū)間徑流量和輸沙量的趨勢(shì)特征

    3.2.1 徑流量變化的趨勢(shì)性特征

    圖3為渭河流域32個(gè)區(qū)間徑流量時(shí)序變化的MK統(tǒng)計(jì)分析與檢驗(yàn)結(jié)果,由圖3可知:除渭12和洛6區(qū)間外,渭河流域93.8%的區(qū)間徑流量呈顯著減小趨勢(shì);除渭8、渭14區(qū)間呈非顯著性減小趨勢(shì)外,其余區(qū)間徑流量均呈顯著減小趨勢(shì),尤以渭3區(qū)間減少最為顯著(Z值為-6.71);涇河的涇1、涇2、涇3、涇4、涇10、涇11區(qū)間徑流量減少比較顯著,尤其是涇1區(qū)間變化最為顯著(Z值為-6.84);北洛河的洛1、洛2區(qū)間徑流量減少比較顯著。

    3.2.2 輸沙量變化的趨勢(shì)性特征

    圖4為不同區(qū)間輸沙量變化的MK統(tǒng)計(jì)檢驗(yàn)結(jié)果,由圖4可知:除渭7、渭8、渭12和洛4區(qū)間輸沙量呈增大趨勢(shì)外,渭河流域87.5%的區(qū)間輸沙量均呈顯著性減小趨勢(shì)。除渭河的渭9、渭10區(qū)間呈非顯著性減小趨勢(shì)外,其余區(qū)間的輸沙量均呈顯著性減小趨勢(shì),尤其是渭3區(qū)間減少幅度最為顯著(Z值為-7.32);涇河的涇1、涇2、涇3、涇4區(qū)間和北洛河的洛1、洛2、洛6區(qū)間的輸沙量也呈顯著減小趨勢(shì)。

    3.3 不同區(qū)間徑流量和輸沙量的突變特征

    對(duì)渭河流域不同區(qū)間的年徑流量和輸沙量進(jìn)行對(duì)比分析,結(jié)果見表4。渭河大部分區(qū)間徑流量和輸沙量的突變年份主要集中在20世紀(jì)80年代末和90年代初;涇河部分區(qū)間徑流量和輸沙量出現(xiàn)突變,且突變年份主要集中在20世紀(jì)90年代;北洛河只有部分區(qū)間徑流量和輸沙量發(fā)生突變,且突變主要發(fā)生在20世紀(jì)90年代末和21世紀(jì)初。

    3.4 渭河流域徑流模數(shù)和輸沙模數(shù)的時(shí)空變化

    渭河流域不同年代的徑流模數(shù)變化差異顯著(見圖5,圖中RDM1~RDM4為不同時(shí)段的徑流模數(shù))。20世紀(jì)70年代,徑流模數(shù)>20萬m3/(km2·a)的強(qiáng)來水區(qū)主要分布在渭7、渭9、渭10、渭14;而渭8和洛7區(qū)間徑流模數(shù)<1.0萬m3/(km2·a),為耗水區(qū)。80年代,整個(gè)流域的徑流模數(shù)有所增大,可能與該時(shí)段的降雨變化有關(guān),徑流模數(shù)>20萬m3/(km2·a)的強(qiáng)來水區(qū)未發(fā)生改變,但渭5、渭6區(qū)間的徑流模數(shù)也比較大,為12萬~20萬m3/(km2·a),該時(shí)段未出現(xiàn)徑流模數(shù)<1萬m3/(km2·a)的區(qū)域。90年代,徑流模數(shù)大范圍減小,由圖5可知徑流模數(shù)>20萬m3/(km2·a)的強(qiáng)來水區(qū)面積明顯減小,徑流模數(shù)為1.0萬~2.5萬m3/(km2·a)的產(chǎn)水區(qū)大幅增大。21世紀(jì)初,部分區(qū)域徑流模數(shù)有所增大,其中渭9、渭10、渭11區(qū)域徑流模數(shù)>20萬m3/(km2·a),為主要產(chǎn)水區(qū)。而洛7區(qū)域的徑流模數(shù)<1.0萬m3/(km2·a),為耗水區(qū)。

    渭河流域不同年代的輸沙模數(shù)(SDM)差異顯著(見圖6)。20世紀(jì)70年代,流域輸沙模數(shù)>8 000 t/(km2·a)的主要產(chǎn)沙中心為洛1、洛2和涇2區(qū)間,輸沙模數(shù)呈現(xiàn)從南到北逐漸增大的趨勢(shì)。80年代,整個(gè)流域的輸沙模數(shù)呈減小趨勢(shì)。90年代,渭河流域輸沙模數(shù)有所增大,且產(chǎn)沙區(qū)域變動(dòng)較大,該時(shí)期輸沙模數(shù)>8 000 t/(km2·a)的區(qū)域?yàn)槁?、洛2、洛3、涇5、涇7、涇8區(qū)間。從上游至下游,渭河干流和涇河的輸沙模數(shù)呈先減小后增大的趨勢(shì),而北洛河輸沙模數(shù)呈減小趨勢(shì)。21世紀(jì)初,渭河流域輸沙模數(shù)大幅減小,除了涇5和涇7區(qū)間外,其余區(qū)間輸沙模數(shù)<5 000 t/(km2·a)。不同時(shí)段內(nèi)渭河流域下游的部分區(qū)間輸沙模數(shù)都處于負(fù)值,且20世紀(jì)90年代后處于負(fù)值的區(qū)間增多,輸沙模數(shù)處于負(fù)值表明下游輸沙量小于上游輸沙量,即該區(qū)間處于泥沙淤積狀態(tài)。研究表明,三門峽水庫修建后河道內(nèi)大量泥沙被攔截,渭河下游處于淤積狀態(tài)[2]。90年代淤積區(qū)域達(dá)到最大,可能是降雨減少、河道輸沙能力降低導(dǎo)致區(qū)間輸沙模數(shù)減小。

    綜上可知,渭河流域不同時(shí)段的徑流和輸沙模數(shù)呈現(xiàn)動(dòng)態(tài)分布。1970—2013年,從上游到下游,渭河干流和北洛河徑流模數(shù)呈先增大后減小的趨勢(shì),涇河呈增大趨勢(shì);整個(gè)流域徑流模數(shù)>12萬m3/(km2·a)的強(qiáng)來水區(qū)面積明顯減小。渭河流域的輸沙模數(shù)在不同的年代呈明顯的動(dòng)態(tài)變化,從1970—2013年,整個(gè)流域的輸沙模數(shù)從南到北呈減小趨勢(shì),表明流域侵蝕有所減弱。

    3.5 渭河流域水沙變化對(duì)人類活動(dòng)和降雨量的響應(yīng)

    渭河和北洛河把口站(華縣和頭)的徑流量—降雨量、輸沙量—降雨量的雙累計(jì)曲線見圖7,曲線斜率分別代表單位降雨產(chǎn)生的徑流量和輸沙量。流域下墊面不變時(shí),雙累計(jì)曲線的斜率不變,若雙累計(jì)曲線斜率發(fā)生變化,則表明流域的下墊面狀況發(fā)生了明顯變化,人類活動(dòng)對(duì)流域水沙的影響增大。渭河流域徑流量和輸沙量與降雨量雙累計(jì)曲線分別在1993年和1996年出現(xiàn)了拐點(diǎn),且拐點(diǎn)后的曲線擬合斜率均減小,說明人類活動(dòng)使得流域的水沙量減小。將突變后的降雨資料代入突變前的回歸方程中,即得徑流量和輸沙量的理論值。不同時(shí)段理論值之差為降水的影響量,突變后的理論值與實(shí)際水沙值之差為人類活動(dòng)的影響量。

    3.6 降雨和人類活動(dòng)對(duì)流域水沙的影響

    對(duì)降雨量—徑流量、降雨量—輸沙量的雙累計(jì)曲線進(jìn)行模擬回歸分析,結(jié)果分別見表5和表6。突變年份之后,實(shí)測(cè)累計(jì)徑流量減小39.72%,實(shí)測(cè)累計(jì)輸沙量減小60.90%。降雨量、人類活動(dòng)對(duì)徑流量和輸沙量變化的貢獻(xiàn)率分別為15.14%、84.86%和7.46%、92.54%。由此可以得出:人類活動(dòng)對(duì)徑流量和輸沙量減少的貢獻(xiàn)率遠(yuǎn)大于降水的,說明人類活動(dòng)在流域水沙的減少上占主導(dǎo)作用。

    流域水沙變化與人類活動(dòng)有著密切的關(guān)系。一方面,人口增長(zhǎng)和城鎮(zhèn)化導(dǎo)致流域土地利用變化顯著,進(jìn)而對(duì)流域水沙產(chǎn)生重要影響。20世紀(jì)八九十年代,流域人口數(shù)量增加導(dǎo)致人們對(duì)耕地和建筑用地的需求增加,90年代后期,隨著退耕還林(草)措施的實(shí)施,流域內(nèi)草地、耕地和建筑用地發(fā)生了較大變化,耕地面積大量減少,林地面積大量增加,使得流域內(nèi)的徑流和泥沙在90年代有所減小[22]。20世紀(jì)后期,隨著區(qū)域內(nèi)經(jīng)濟(jì)發(fā)展和城鎮(zhèn)化,耕地和建筑用地面積分別增加453、679 km2[23]。不同程度的土地利用變化都對(duì)流域水沙變化產(chǎn)生了影響,且21世紀(jì)00年代土地利用變化對(duì)徑流的影響程度大于20世紀(jì)90年代[24]。另一方面,改革開放后,渭河流域經(jīng)濟(jì)增長(zhǎng)迅速,工農(nóng)業(yè)生產(chǎn)及生活用水增加。據(jù)統(tǒng)計(jì),渭河流域90年代以前平均年耗水量為27.93億m3,90年代平均年耗水量增加至42.63億m3,耗水量增加了52.6%,增加的耗水量主要來自地下水的過度開采,導(dǎo)致地表徑流也受到影響,流域年入黃水量減少了8.64億m3[25]。截至1996年底,渭河上中游地區(qū)有效灌溉面積達(dá)80 547 hm2(其中引水灌區(qū)占54.9%),中型灌區(qū)33處,其中引水灌區(qū)22處;下游萬畝(15畝=1 hm2)以上灌區(qū)達(dá)99處,總有效灌溉面積590 886 hm2(其中引水灌區(qū)占50%左右),其中河道引水灌區(qū)50處。工農(nóng)業(yè)生產(chǎn)用水的增加進(jìn)一步使得流域內(nèi)的徑流減少,引水灌溉導(dǎo)致的間接引沙使得流域輸沙量減小[26]。

    20世紀(jì)70年代黃土高原開始實(shí)施大范圍的水土保持措施,淤地壩作為一種重要的攔水?dāng)r沙措施,在渭河流域水沙變化中也發(fā)揮了重要作用。截至1996年年底,渭河流域內(nèi)已建成2 336座淤地壩,相對(duì)20世紀(jì)70年代,90年代的流域輸沙量減少866.56萬t。另外,水庫及梯田的修建也使得流域內(nèi)的水沙受到了明顯影響,80年代水利工程平均每年減沙2 836萬t,減少徑流22.7億m3,梯田、造林等減沙1 744萬t,減少徑流1.311億m3[27]。進(jìn)入21世紀(jì)后整個(gè)流域的降雨量增加,但是徑流量卻依然減少,這與當(dāng)時(shí)流域內(nèi)實(shí)施的一系列水土保持措施和政策有關(guān),1999年退耕還林政策的開始實(shí)施以及2005年《渭河流域重點(diǎn)治理規(guī)劃》的正式批復(fù)實(shí)施,均使得流域內(nèi)的水土保持生態(tài)建設(shè)加強(qiáng),水沙量進(jìn)一步減少。另外,由于部分造林等水土保持措施對(duì)流域內(nèi)水沙的控制有一定的滯后效應(yīng),因此部分區(qū)間到了21世紀(jì)水沙量才出現(xiàn)明顯的改變[28]。

    4 結(jié) 論

    (1)渭河流域93.8%的區(qū)間徑流量呈顯著性減小趨勢(shì),87.5%的區(qū)間輸沙量呈顯著性減小趨勢(shì),渭河徑流量和輸沙量的突變主要發(fā)生在20世紀(jì)80年代末和90年代初;涇河區(qū)間徑流量和輸沙量的突變主要發(fā)生在20世紀(jì)90年代;北洛河徑流量和輸沙量的突變主要發(fā)生在21世紀(jì)初。

    (2)渭河流域不同年代徑流和輸沙模數(shù)呈現(xiàn)動(dòng)態(tài)分布,1970—2013年,渭河干流和北洛河徑流模數(shù)呈先增大后減小的趨勢(shì),涇河呈增大趨勢(shì);整個(gè)流域的輸沙模數(shù)從南到北呈逐年減小趨勢(shì)。

    (3)雙累計(jì)曲線法表明人類活動(dòng)對(duì)渭河流域水沙有著顯著影響,且影響量大于降雨對(duì)水沙的影響。人類活動(dòng)對(duì)徑流變化的貢獻(xiàn)率達(dá)84.86%,對(duì)輸沙量變化的貢獻(xiàn)率高達(dá)92.54%;降雨對(duì)徑流量的影響為15.14%,對(duì)輸沙的影響為7.46%。

    參考文獻(xiàn):

    [1] 劉曉瓊,劉彥隨,李同昇,等.基于小波多尺度變換的渭河水沙演變規(guī)律研究[J].地理科學(xué),2015,35(2):211-217.

    [2] GAO P, GEISSEN V, RITSEM C J, et al. Impact of Climate Change and Anthropogenic Activities on Stream Flow and Sediment Discharge in the Wei River Basin, China[J]. Hydrology and Earth System Sciences, 2013, 17: 961-972.

    [3] 李楊俊,段勛年,劉福勤,等.渭河下游河道萎縮特性分析和改善對(duì)策[J].人民黃河,1998,20(7):9-11.

    [4] 陜西救災(zāi)年鑒編委會(huì).陜西救災(zāi)年鑒(1996)[M].西安: 陜西人民出版社,2000:22-23.

    [5] 水利部水土保持司.渭河流域“03·8”特大暴雨洪水及水土保持調(diào)研報(bào)告[R].北京:水利部水土保持司,2003:23-27.

    [6] 王小永.涇、渭河洪水遭遇分析淺議[J].陜西水利,2011(1):115-116.

    [7] 冉大川,劉斌,王宏.水土保持措施對(duì)黃河減水減沙作用的分析[J].中國水土保持,2002(10):38-39,53.

    [8] 胡安焱,劉燕,郭生練,等.渭河流域水沙多年變化及趨勢(shì)分析[J].人民黃河,2007,29(2):39-41.

    [9] 任立良,沈鴻仁,袁飛,等.變化環(huán)境下渭河流域水文干旱演變特征剖析[J].水科學(xué)進(jìn)展,2016,27(4):492-500.

    [10] 焦菊英,馬祥華,王飛,等.渭河流域侵蝕產(chǎn)沙強(qiáng)度的區(qū)域分異特征[J].水土保持研究,2004,11(4):60-63.

    [11] 舒長(zhǎng)先.黃河中游地區(qū)泥沙變化研究進(jìn)展[J].水土保持通報(bào),2000,20(3):10-14.

    [12] 侯欽磊,白紅英,任園園,等.50年來渭河干流徑流變化及其驅(qū)動(dòng)力分析[J].資源科學(xué),2011,33(8):1505-1512.

    [13] YUE S, PILON P, CAVADIAS G. Power of the Mann-Kendall and Spearmans RHO Tests for Detecting Monotonic Trends in Hydrological Series[J]. Journal of Hydrology, 2002, 259(1): 254.

    [14] MANN H B. Nonparametric Tests Against Trend[J]. Econometrica, 1945, 13(3): 245-259.

    [15] 趙廣舉,穆興民,田鵬.近60年黃河中游水沙變化趨勢(shì)及其影響因素分析[J].資源科學(xué),2012,34(6):1070-1078.

    [16] ZHANG S R, LU X X. Hydrological Responses to Precipitation Variation and Diverse Human Activities in a Mountainous Tributary of the Lower Xijiang, China[J]. Catena, 2009, 77(2): 130-142.

    [17] PETTITT A N. A Non 2 Parametric Approach to the Change Point Problem[J]. Applied Statistics, 1979, 28 (2): 126-135.

    [18] SEARCY J K, HARDISON C H. Double-Mass Curves[M]. Washington D C: U.S.Geological Survey Water Supply Paper 1541-B, 1960:37.

    [19] GAO P, LI P, ZHAO B, et al. Use of Double Mass Curves in Hydrologic Benefit Evaluations[J]. Hydrol. Process, 2017, 31 (26): 4639-4646.

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

    [21] 王宏,楊國禮,王瑞芳.渭河流域水利水保措施對(duì)泥沙、徑流影響分析計(jì)算[J].水土保持通報(bào),1994,14(5):48-52.

    [22] 蔣觀滔,高鵬,穆興民,等.退耕還林(草)對(duì)北洛河上游水沙變化的影響[J].水土保持研究,2015,22(6):1-6.

    [23] 柴雪柯,蔣觀滔,高鵬,等.渭河流域土地利用變化特征及驅(qū)動(dòng)力分析[J].人民黃河,2016,38(8):77-81,85.

    [24] 黎云云,暢建霞,王義民,等.渭河流域徑流對(duì)土地利用變化的時(shí)空響應(yīng)[J].農(nóng)業(yè)工程學(xué)報(bào),2016,32(15):232-238.

    [25] 趙雙全,屠新武,楊術(shù)容,等.渭河下游水沙衰減及對(duì)黃河的影響[J].電網(wǎng)與清潔能源,2009,25(5):55-58.

    [26] 趙俊俠,王宏,馬勇,等.1990—1996年渭河流域水沙變化原因初步分析[J].水土保持學(xué)報(bào),2001(增刊2):136-139.

    [27] 張淑蘭,王彥輝,于澎濤,等.涇河流域近50年來的徑流時(shí)空變化與驅(qū)動(dòng)力分析[J].地理科學(xué),2011,31(6):721-727.

    [28] 秦偉,朱清科,劉廣全,等.北洛河上游生態(tài)建設(shè)的水沙調(diào)控效應(yīng)[J].水利學(xué)報(bào),2010,41(11):1325-1332.

    【責(zé)任編輯 翟戌亮】

    收稿日期:2019-02-27

    基金項(xiàng)目:國家科技基礎(chǔ)性工作專項(xiàng)(2014FY210120);國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2016YFC0501707);國家自然科學(xué)基金資助項(xiàng)目(41371277)

    作者簡(jiǎn)介:徐瑞瑞(1994—),女,甘肅臨夏人,碩士研究生,研究方向?yàn)榱饔蛏鷳B(tài)水文和水土保持

    通信作者:高鵬(1976—),男,陜西扶風(fēng)人,研究員,博士,博士生導(dǎo)師,主要從事流域生態(tài)水文和水土保持研究工作

    E-mail:gaopeng@ms.iswc.ac.cn

    猜你喜歡
    人類活動(dòng)渭河流域輸沙量
    遼河干流遼中站河道沖淤變化趨勢(shì)研究
    氣候變化和人類活動(dòng)對(duì)祖厲河輸沙量變化的影響分析
    頻率曲線在荊江三口輸沙量計(jì)算中的應(yīng)用
    氣候變化和人類活動(dòng)對(duì)漢江上游徑流變化影響的定量研究
    氣候變化和人類活動(dòng)對(duì)漢江上游徑流變化影響的定量研究
    基于RS/GIS 渭河流域植被覆蓋時(shí)空變化特征研究
    燕太子回國
    關(guān)于遼河干流河道沖淤量沿程變化規(guī)律探討
    地下水(2015年5期)2015-12-02 02:44:26
    渭河流域香菜夏秋無公害栽培技術(shù)
    人為因素主導(dǎo)下海岸帶生態(tài)系統(tǒng)變遷探討
    綠色科技(2015年6期)2015-08-05 18:01:22
    2018国产大陆天天弄谢| 国产精品久久久久久精品电影小说| 如日韩欧美国产精品一区二区三区| 大陆偷拍与自拍| 国产午夜精品一二区理论片| 天天影视国产精品| 久久人人爽av亚洲精品天堂| 亚洲精品中文字幕在线视频| 无限看片的www在线观看| 狂野欧美激情性xxxx| 日韩欧美一区视频在线观看| 啦啦啦中文免费视频观看日本| 黄色怎么调成土黄色| 午夜福利视频精品| 高潮久久久久久久久久久不卡| 国产高清国产精品国产三级| 搡老乐熟女国产| 国产精品一区二区免费欧美 | 久久久国产一区二区| 国产成人91sexporn| 亚洲欧洲国产日韩| 国产老妇伦熟女老妇高清| 久久精品久久精品一区二区三区| 色综合欧美亚洲国产小说| 91麻豆av在线| 各种免费的搞黄视频| 最新在线观看一区二区三区 | 美女福利国产在线| 国产99久久九九免费精品| 黄色怎么调成土黄色| 亚洲欧美日韩高清在线视频 | 大香蕉久久网| a级片在线免费高清观看视频| 亚洲第一青青草原| 亚洲精品久久成人aⅴ小说| 国产在线观看jvid| 日本91视频免费播放| 精品亚洲乱码少妇综合久久| 男人舔女人的私密视频| 亚洲精品中文字幕在线视频| 亚洲成人免费电影在线观看 | 久久天躁狠狠躁夜夜2o2o | 一本大道久久a久久精品| 91九色精品人成在线观看| 日日夜夜操网爽| 亚洲 欧美一区二区三区| 国产精品一区二区免费欧美 | 欧美久久黑人一区二区| 天堂8中文在线网| 高潮久久久久久久久久久不卡| 精品一区二区三区av网在线观看 | 久久人人97超碰香蕉20202| 女人高潮潮喷娇喘18禁视频| 日本五十路高清| 丝袜喷水一区| 亚洲精品久久成人aⅴ小说| 蜜桃在线观看..| 99国产精品免费福利视频| 性高湖久久久久久久久免费观看| 黑丝袜美女国产一区| 51午夜福利影视在线观看| a级毛片在线看网站| 国产有黄有色有爽视频| 久久精品人人爽人人爽视色| 人人妻人人澡人人爽人人夜夜| 欧美日韩综合久久久久久| 欧美激情高清一区二区三区| 成年美女黄网站色视频大全免费| 肉色欧美久久久久久久蜜桃| cao死你这个sao货| 亚洲av国产av综合av卡| 久久免费观看电影| 制服诱惑二区| 香蕉国产在线看| 视频区图区小说| 一级片'在线观看视频| 亚洲综合色网址| 亚洲欧洲日产国产| 久久国产精品人妻蜜桃| 色视频在线一区二区三区| 啦啦啦在线免费观看视频4| 人妻人人澡人人爽人人| 国产精品熟女久久久久浪| 欧美精品亚洲一区二区| 久久人人爽人人片av| 久久国产精品影院| 少妇粗大呻吟视频| 久久久久久久国产电影| 成人国产一区最新在线观看 | 99国产精品一区二区蜜桃av | 精品一区二区三区四区五区乱码 | 99九九在线精品视频| 老司机影院成人| a级毛片黄视频| 亚洲熟女毛片儿| 蜜桃国产av成人99| 午夜免费成人在线视频| 巨乳人妻的诱惑在线观看| 1024香蕉在线观看| 久久久国产精品麻豆| 国产在线免费精品| 国产三级黄色录像| 欧美激情极品国产一区二区三区| www.av在线官网国产| 久久国产精品大桥未久av| 黑丝袜美女国产一区| 亚洲九九香蕉| 精品福利永久在线观看| 成年人黄色毛片网站| 日本一区二区免费在线视频| 99精品久久久久人妻精品| 国产极品粉嫩免费观看在线| 麻豆国产av国片精品| av片东京热男人的天堂| 国产精品亚洲av一区麻豆| 国产主播在线观看一区二区 | 一级片'在线观看视频| 中国国产av一级| 高潮久久久久久久久久久不卡| 99久久99久久久精品蜜桃| 国产精品久久久久久精品电影小说| 飞空精品影院首页| 国产一区有黄有色的免费视频| 深夜精品福利| 国产精品一区二区在线观看99| 51午夜福利影视在线观看| 国产色视频综合| 久久久久久亚洲精品国产蜜桃av| 久久久久精品国产欧美久久久 | 欧美成狂野欧美在线观看| 成人国产av品久久久| 女人被躁到高潮嗷嗷叫费观| 自线自在国产av| 搡老岳熟女国产| 国产不卡av网站在线观看| 午夜福利,免费看| 精品国产乱码久久久久久男人| 满18在线观看网站| 超色免费av| 日本wwww免费看| 久久久国产一区二区| 人人妻人人添人人爽欧美一区卜| 女人爽到高潮嗷嗷叫在线视频| 亚洲欧美一区二区三区国产| 免费在线观看影片大全网站 | 久久精品国产亚洲av涩爱| 90打野战视频偷拍视频| 久久久亚洲精品成人影院| 欧美日韩成人在线一区二区| 别揉我奶头~嗯~啊~动态视频 | 99热网站在线观看| 日韩中文字幕视频在线看片| 丰满人妻熟妇乱又伦精品不卡| 丰满迷人的少妇在线观看| 天天躁日日躁夜夜躁夜夜| 水蜜桃什么品种好| 麻豆乱淫一区二区| 亚洲精品在线美女| 免费少妇av软件| 1024香蕉在线观看| 国产欧美亚洲国产| 亚洲av电影在线进入| 99国产精品一区二区蜜桃av | 波多野结衣av一区二区av| 无遮挡黄片免费观看| 国产欧美亚洲国产| 亚洲精品日韩在线中文字幕| 国产国语露脸激情在线看| 国产一区亚洲一区在线观看| 国产高清videossex| 一本综合久久免费| 久久久精品区二区三区| 赤兔流量卡办理| 午夜福利免费观看在线| 十八禁人妻一区二区| 国产成人啪精品午夜网站| 99国产精品99久久久久| 婷婷色综合www| 丝袜喷水一区| 色网站视频免费| 最黄视频免费看| av国产精品久久久久影院| 国产亚洲av高清不卡| 少妇精品久久久久久久| 宅男免费午夜| 国产免费福利视频在线观看| 亚洲精品日韩在线中文字幕| 精品第一国产精品| 国产精品熟女久久久久浪| 新久久久久国产一级毛片| 成人亚洲精品一区在线观看| 亚洲美女黄色视频免费看| 大码成人一级视频| 色视频在线一区二区三区| 亚洲成人免费电影在线观看 | 两人在一起打扑克的视频| 脱女人内裤的视频| 一边摸一边做爽爽视频免费| 巨乳人妻的诱惑在线观看| videosex国产| 制服诱惑二区| 欧美精品av麻豆av| 超碰成人久久| 国产在线一区二区三区精| 欧美激情 高清一区二区三区| 国产有黄有色有爽视频| 久久久精品国产亚洲av高清涩受| 久久久精品94久久精品| 一本久久精品| 国产日韩欧美亚洲二区| 美女高潮到喷水免费观看| 亚洲精品久久久久久婷婷小说| 国产亚洲午夜精品一区二区久久| 色播在线永久视频| 国产亚洲av高清不卡| 波多野结衣av一区二区av| 久久女婷五月综合色啪小说| 1024视频免费在线观看| www.精华液| a 毛片基地| 精品一区二区三卡| 国产亚洲精品久久久久5区| 亚洲av成人精品一二三区| 两个人看的免费小视频| 欧美日韩av久久| 丝袜喷水一区| 天天添夜夜摸| 亚洲欧洲国产日韩| 一区二区三区四区激情视频| 激情视频va一区二区三区| 国产91精品成人一区二区三区 | 亚洲成人免费电影在线观看 | 亚洲国产中文字幕在线视频| 女人爽到高潮嗷嗷叫在线视频| 日韩制服丝袜自拍偷拍| 欧美激情极品国产一区二区三区| 后天国语完整版免费观看| 91老司机精品| 亚洲精品在线美女| 欧美人与善性xxx| 最新的欧美精品一区二区| 亚洲成国产人片在线观看| 国产成人a∨麻豆精品| 亚洲中文字幕日韩| 人妻 亚洲 视频| 伊人亚洲综合成人网| 成年女人毛片免费观看观看9 | 精品熟女少妇八av免费久了| 成人18禁高潮啪啪吃奶动态图| 久久人人爽人人片av| 午夜免费男女啪啪视频观看| 男的添女的下面高潮视频| 18禁裸乳无遮挡动漫免费视频| 国产免费一区二区三区四区乱码| 老司机靠b影院| a级片在线免费高清观看视频| 亚洲精品日韩在线中文字幕| 伦理电影免费视频| 国产黄色视频一区二区在线观看| av欧美777| 精品福利永久在线观看| 国产亚洲午夜精品一区二区久久| 午夜激情av网站| 一区二区日韩欧美中文字幕| 电影成人av| 国产日韩一区二区三区精品不卡| 亚洲欧美日韩另类电影网站| 一区二区日韩欧美中文字幕| 天堂中文最新版在线下载| 色综合欧美亚洲国产小说| 777久久人妻少妇嫩草av网站| 亚洲第一av免费看| 黄色一级大片看看| 日韩伦理黄色片| 一区二区三区四区激情视频| 18禁观看日本| 欧美日本中文国产一区发布| 亚洲色图 男人天堂 中文字幕| 在线观看国产h片| 久久久精品区二区三区| www.999成人在线观看| 欧美成狂野欧美在线观看| 中文字幕av电影在线播放| 精品亚洲乱码少妇综合久久| 人体艺术视频欧美日本| 免费不卡黄色视频| 久久精品国产综合久久久| 天天影视国产精品| 亚洲av欧美aⅴ国产| 国产激情久久老熟女| 亚洲精品日本国产第一区| 国产主播在线观看一区二区 | 高清欧美精品videossex| 久久久国产欧美日韩av| 肉色欧美久久久久久久蜜桃| 日韩av免费高清视频| 亚洲精品第二区| 18禁国产床啪视频网站| 亚洲人成电影免费在线| 国产1区2区3区精品| 亚洲av成人不卡在线观看播放网 | 亚洲中文av在线| 性少妇av在线| 99九九在线精品视频| 亚洲av成人不卡在线观看播放网 | 人体艺术视频欧美日本| 精品国产乱码久久久久久小说| 国精品久久久久久国模美| 亚洲国产欧美日韩在线播放| av有码第一页| 亚洲国产av影院在线观看| 咕卡用的链子| 一个人免费看片子| 亚洲黑人精品在线| 日本猛色少妇xxxxx猛交久久| 亚洲精品成人av观看孕妇| 亚洲久久久国产精品| 观看av在线不卡| 99久久综合免费| 亚洲国产看品久久| 亚洲激情五月婷婷啪啪| 青春草亚洲视频在线观看| 亚洲精品在线美女| 90打野战视频偷拍视频| 一本—道久久a久久精品蜜桃钙片| 色婷婷久久久亚洲欧美| 赤兔流量卡办理| 麻豆国产av国片精品| 国产av国产精品国产| 18禁观看日本| 久久中文字幕一级| 午夜激情久久久久久久| 国产视频一区二区在线看| 国产在线视频一区二区| 18禁黄网站禁片午夜丰满| 精品卡一卡二卡四卡免费| 男女无遮挡免费网站观看| 18在线观看网站| 中文字幕高清在线视频| 国产精品麻豆人妻色哟哟久久| 一级,二级,三级黄色视频| 青青草视频在线视频观看| 两个人免费观看高清视频| cao死你这个sao货| 精品久久蜜臀av无| 一级黄片播放器| 亚洲成色77777| 自拍欧美九色日韩亚洲蝌蚪91| a级毛片黄视频| 少妇被粗大的猛进出69影院| 欧美日韩视频精品一区| 久久ye,这里只有精品| 免费看av在线观看网站| www.999成人在线观看| 日本wwww免费看| 9191精品国产免费久久| 国产xxxxx性猛交| 只有这里有精品99| av国产精品久久久久影院| 免费少妇av软件| 汤姆久久久久久久影院中文字幕| 久久久久久久精品精品| 18禁观看日本| 男女免费视频国产| 久久人人爽av亚洲精品天堂| 国产在线一区二区三区精| 久久99一区二区三区| 在线观看国产h片| 黑人巨大精品欧美一区二区蜜桃| 大话2 男鬼变身卡| 搡老乐熟女国产| 国产欧美日韩一区二区三区在线| 天堂中文最新版在线下载| 大片电影免费在线观看免费| 建设人人有责人人尽责人人享有的| 色94色欧美一区二区| 亚洲精品国产一区二区精华液| 亚洲中文字幕日韩| 国产免费现黄频在线看| 爱豆传媒免费全集在线观看| 啦啦啦啦在线视频资源| 波野结衣二区三区在线| 黑人猛操日本美女一级片| 一区二区三区激情视频| 女人久久www免费人成看片| 啦啦啦在线观看免费高清www| 激情视频va一区二区三区| 亚洲av电影在线进入| 国产精品欧美亚洲77777| 亚洲图色成人| 1024视频免费在线观看| 亚洲精品av麻豆狂野| 亚洲欧美色中文字幕在线| 国产精品偷伦视频观看了| 丰满少妇做爰视频| 另类亚洲欧美激情| 中文字幕人妻熟女乱码| 欧美国产精品va在线观看不卡| 一本一本久久a久久精品综合妖精| 日本vs欧美在线观看视频| 日本欧美国产在线视频| 国产伦理片在线播放av一区| 99热网站在线观看| 欧美日韩精品网址| av有码第一页| 国产精品三级大全| 婷婷成人精品国产| 亚洲欧美精品自产自拍| 9色porny在线观看| 精品人妻1区二区| 中国美女看黄片| 日日摸夜夜添夜夜爱| 一个人免费看片子| 免费人妻精品一区二区三区视频| 人妻人人澡人人爽人人| 免费看不卡的av| 性少妇av在线| 国产视频一区二区在线看| 免费高清在线观看日韩| 亚洲av电影在线观看一区二区三区| 欧美成人精品欧美一级黄| 高清黄色对白视频在线免费看| 一级毛片女人18水好多 | 午夜av观看不卡| 中文字幕精品免费在线观看视频| a级片在线免费高清观看视频| e午夜精品久久久久久久| 久久人人爽人人片av| 男女免费视频国产| 久久精品久久精品一区二区三区| 80岁老熟妇乱子伦牲交| 国产麻豆69| 精品少妇一区二区三区视频日本电影| 国产高清不卡午夜福利| 久久影院123| 国产在线观看jvid| 男人操女人黄网站| 亚洲一区中文字幕在线| 黄色一级大片看看| 欧美国产精品一级二级三级| 乱人伦中国视频| 青草久久国产| 国产黄色视频一区二区在线观看| 看十八女毛片水多多多| 成人影院久久| 天天躁夜夜躁狠狠躁躁| 国产精品.久久久| 成年美女黄网站色视频大全免费| 一级,二级,三级黄色视频| 久久精品亚洲av国产电影网| 别揉我奶头~嗯~啊~动态视频 | 少妇被粗大的猛进出69影院| 国产有黄有色有爽视频| 首页视频小说图片口味搜索 | 人成视频在线观看免费观看| 青草久久国产| 2021少妇久久久久久久久久久| 深夜精品福利| 一边亲一边摸免费视频| 在线亚洲精品国产二区图片欧美| 天天躁夜夜躁狠狠久久av| 亚洲欧美清纯卡通| 看免费成人av毛片| 国产精品三级大全| 国产亚洲av片在线观看秒播厂| 国产色视频综合| 超碰成人久久| 精品一区二区三区av网在线观看 | 爱豆传媒免费全集在线观看| 久久精品国产a三级三级三级| 青草久久国产| 精品久久久精品久久久| 捣出白浆h1v1| 欧美大码av| 欧美97在线视频| 国产一区二区在线观看av| 另类精品久久| 伊人亚洲综合成人网| 高清不卡的av网站| 国产不卡av网站在线观看| 啦啦啦在线观看免费高清www| 国产欧美日韩一区二区三区在线| 日韩av不卡免费在线播放| 亚洲欧美清纯卡通| 黄色视频不卡| 精品一区在线观看国产| 国产成人系列免费观看| 日韩伦理黄色片| 国产亚洲av高清不卡| 亚洲av日韩在线播放| 亚洲欧美一区二区三区久久| 国产一区二区三区综合在线观看| 制服人妻中文乱码| 日韩欧美一区视频在线观看| 一本—道久久a久久精品蜜桃钙片| 亚洲精品中文字幕在线视频| 国产精品久久久久久人妻精品电影 | 多毛熟女@视频| 丝袜喷水一区| 操美女的视频在线观看| 18禁观看日本| 欧美日韩视频精品一区| 一级黄色大片毛片| 亚洲av日韩在线播放| 国产成人a∨麻豆精品| 国产精品二区激情视频| 国产成人啪精品午夜网站| 午夜影院在线不卡| 十分钟在线观看高清视频www| 亚洲 国产 在线| 观看av在线不卡| 久久精品久久久久久久性| 久久精品成人免费网站| 美女国产高潮福利片在线看| 中文字幕另类日韩欧美亚洲嫩草| 精品少妇内射三级| 成年人免费黄色播放视频| 黄色怎么调成土黄色| 人人妻人人爽人人添夜夜欢视频| 国产一区亚洲一区在线观看| 极品少妇高潮喷水抽搐| 99香蕉大伊视频| 欧美亚洲日本最大视频资源| 亚洲精品自拍成人| 国产片特级美女逼逼视频| 国产亚洲精品久久久久5区| 9191精品国产免费久久| 国产一区二区三区综合在线观看| 精品第一国产精品| a级片在线免费高清观看视频| 丰满人妻熟妇乱又伦精品不卡| 欧美黑人精品巨大| 啦啦啦在线观看免费高清www| 国产精品久久久av美女十八| 少妇人妻 视频| 久久久久久久国产电影| 亚洲天堂av无毛| 晚上一个人看的免费电影| 99精品久久久久人妻精品| 婷婷成人精品国产| 欧美成狂野欧美在线观看| 欧美日韩黄片免| 亚洲av国产av综合av卡| 黄色视频在线播放观看不卡| 亚洲欧美精品自产自拍| 2021少妇久久久久久久久久久| 国产亚洲精品久久久久5区| 无限看片的www在线观看| 肉色欧美久久久久久久蜜桃| 老司机深夜福利视频在线观看 | 亚洲人成77777在线视频| 中文字幕av电影在线播放| 日韩,欧美,国产一区二区三区| 脱女人内裤的视频| 老汉色av国产亚洲站长工具| 久久人人爽av亚洲精品天堂| 亚洲国产精品999| 亚洲情色 制服丝袜| 深夜精品福利| 国产日韩一区二区三区精品不卡| 国产精品久久久人人做人人爽| av一本久久久久| 久久天堂一区二区三区四区| 激情五月婷婷亚洲| 国产91精品成人一区二区三区 | 老汉色av国产亚洲站长工具| 午夜免费鲁丝| 99国产精品免费福利视频| 校园人妻丝袜中文字幕| 欧美人与善性xxx| 制服人妻中文乱码| 宅男免费午夜| 黑人欧美特级aaaaaa片| 久久精品aⅴ一区二区三区四区| 久久久国产一区二区| h视频一区二区三区| 国产人伦9x9x在线观看| 亚洲av美国av| 黑人欧美特级aaaaaa片| 亚洲成av片中文字幕在线观看| 免费在线观看日本一区| 日本黄色日本黄色录像| 日本av手机在线免费观看| 国产在线视频一区二区| 中文字幕高清在线视频| 国产在线视频一区二区| 亚洲色图综合在线观看| 亚洲中文av在线| 人妻人人澡人人爽人人| 满18在线观看网站| 热99国产精品久久久久久7| 男的添女的下面高潮视频| 9热在线视频观看99| 国产成人精品久久二区二区91| 亚洲一区中文字幕在线| 男女边摸边吃奶| 蜜桃国产av成人99| 久久精品久久久久久噜噜老黄| 婷婷色av中文字幕| 国产成人一区二区三区免费视频网站 | a 毛片基地| 18在线观看网站| 亚洲国产欧美日韩在线播放| 国产成人av教育| 亚洲国产毛片av蜜桃av| 国产亚洲欧美精品永久| 另类亚洲欧美激情| 校园人妻丝袜中文字幕| 日韩熟女老妇一区二区性免费视频| 午夜日韩欧美国产| 亚洲中文av在线| 午夜91福利影院| 亚洲黑人精品在线| 久久精品久久久久久久性| 中文字幕人妻熟女乱码|