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

    烏江流域水沙變化趨勢及驅動力分析

    2021-10-18 00:57:31郭文獻趙瑞超付體江古今用
    人民長江 2021年9期
    關鍵詞:影響

    郭文獻 趙瑞超 付體江 古今用

    摘要:為對烏江流域水文變化趨勢特征及其主要驅動力進行分析,采用Mann-Kendall非參數(shù)檢驗,滑動T檢驗和累積曲線等方法研究了烏江武隆站1965~2017年輸沙量、徑流量、降水量及其變化趨勢與突變特征,采用累積量斜率變化率法定量計算了烏江流域氣候變化和人類活動對輸沙量與徑流量變化的貢獻率,并通過河流影響因子法(River Impact,RI)對烏江水沙情勢改變度進行了分析。結果表明:烏江流域徑流量、降水量、輸沙量均呈下降趨勢,其中輸沙量的下降幅度最為顯著。導致徑流量和輸沙量下降的主要驅動因素分別是氣候變化與人類活動。烏江流域徑流量發(fā)生輕度改變而輸沙量發(fā)生嚴重改變。

    關 鍵 詞:水沙情勢;水文改變度;累積量斜率變化率法;烏江流域;氣候變化

    中圖法分類號:X826

    文獻標志碼:A

    文章編號:1001-4179(2021)09-0071-08

    DOI:10.16232/j.cnki.1001-4179.2021.09.012

    0 引 言

    烏江全長1 050 km,流域總面積為87 920 km2,橫貫貴州省西部、中部和東北部及四川省東部,是長江上游右岸的最大支流,也是三峽水庫最主要的水源補給河流之一[1]。烏江河口距三峽大壩僅489 km[2],也是中國十三大水電基地之一[3]。1982~2013年烏江干流上先后修建了11座大型水庫,這些大型水利工程的修建勢必會在一定程度上引起烏江水沙特性發(fā)生改變[4]。近幾十年來,烏江流域水文情勢改變所產(chǎn)生的負面影響也越來越受到重視。國內(nèi)同時出現(xiàn)了大量的科學文獻研究了人類活動對河流水沙狀況的影響,如吳曉玲等[5]利用水文資料對烏江上游水沙態(tài)勢變異特征及其成因進行系統(tǒng)分析;陳松生等[6]初步評價了烏江流域的水沙特性及其原因;班璇等[7]采用了累積量斜率變化率法對漢江流域氣象水文變化驅動力進行了分析;熊亞蘭等[8]利用Mann-Kendall分析烏江水文特征和變化趨勢,并采用雙累積曲線法分析人類活動和氣候變化對烏江水沙的影響。管曉祥等[9]采用Mann-Kendall、均值差異T等進一步分析了烏江流域水沙演變規(guī)律。但這些研究主要針對河流的水文年際變化特征,缺乏氣候變化和人類活動對河流水沙情勢改變的定量評價。

    本文采用烏江流域武隆站實測徑流、泥沙數(shù)據(jù)和流域內(nèi)多年降水資料,在考慮到氣候變化中降水對于河流水沙情勢的影響最為顯著,而氣溫、蒸發(fā)量等氣候因素對河流水文情勢的影響幾乎在所有的流域中都明顯受到降水差異的中和[10-11],因此本文將降水作為影響烏江水沙情勢的主要氣候因子。在分析近幾十年來烏江流域水沙年際變化特征的基礎上,定量評價了氣候變化和人類活動對烏江水沙情勢的貢獻率及改變程度,本次研究可為烏江流域河流生態(tài)保護研究提供參考和基礎。

    1 數(shù)據(jù)資料和方法

    1.1 研究區(qū)域及數(shù)據(jù)

    研究區(qū)域為烏江流域,流域集水面積87 920 km2,干流全長107 km,天然落差2 124 m,武隆水文站位于烏江下游,為烏江匯入長江的控制水文站,因此研究采用武隆水文站1965~2017年水文數(shù)據(jù)及流域內(nèi)氣象站多年降水數(shù)據(jù)資料分析近幾十年烏江水文情勢變化特征。其中徑流量、輸沙量數(shù)據(jù)來源于長江流域水文年鑒,降水資料來源于中國氣象數(shù)據(jù)網(wǎng)(http:∥data.cma.cn/)。

    1.2 趨勢分析及突變檢驗

    Mann-Kendall趨勢檢驗通過計算時間序列數(shù)據(jù)的標準化變量Z,與某一置信水平α下的臨界變量對比。計算標準變量Z為正表示有上升趨勢,反之表示有下降趨勢;當Z超出臨界值時表明上升或下降趨勢顯著[12]。同時對原時間序列的逆序列進行同樣的統(tǒng)計量計算,使UB=UF,若兩條曲線在95%置信度水平內(nèi)出現(xiàn)交點,表明在該時間點發(fā)生突變,方法具體計算見參考文獻[13-14]。Mann-Kendall非參數(shù)檢驗不受樣本值和分布類型的干擾,但是檢驗過程中可能出現(xiàn)多個突變點,需要對這些突變點進行驗證。而滑動T檢驗法是通過考察兩組樣本平均值的差異是否顯著來檢驗突變。其基本思想是把一氣候序列中兩段子序列均值有無顯著差異視為來自兩個總體均值有無顯著差異的問題來檢驗。如果兩段子序列的均值差異超過了一定的顯著水平,則可認為有突變發(fā)生 [15],即可以計算某一基準點的統(tǒng)計量t并與臨界值對比進行突變性檢驗[16]。本次研究采用Mann-Kendall法、滑動T檢驗法進行突變檢驗。

    1.3 相關性分析

    相關性分析是研究隨機變量間相關關系的統(tǒng)計方法[17]。水文變量之間的相關分析主要通過計算相關系數(shù)揭示不同要素之間的密切程度。將兩組隨機變量xi和yi(i=1,2,…,n)的相關系數(shù)R定義為

    式中:為變量xi的平均值;為變量yi的平均值。當R≥0.8時表示為高度相關;當0.5≤R<0.8時表示為中度相關;當0.3≤R<0.5時表示為低度相關;當R<0.3時表示為相關程度極低。

    1.4 雙累積曲線法

    雙累積曲線是檢驗兩個參數(shù)間關系一致性及變化的方法,通過在直角坐標系中繪制同時期內(nèi)的一個變量連續(xù)累積值與另一個變量的連續(xù)累積值關系曲線,可用于水文氣象要素一致性檢驗、缺值插補及趨勢性變化和強度分析[18-19]。

    1.5 河流影響因子

    河流影響因子法主要考慮水文變量年際變化、極值發(fā)生時間變化、年內(nèi)變化,將這3類影響因子定義為MIF(Magnitude Impact Factor,量級影響因子)、TIF(Timing Impact Factor,時間影響因子)和VIF(Variation Impact Factor,年內(nèi)變化影響因子)并進行量化。MIF作為量級影響因子,其影響等同于TIF和VIF共同作用。河流影響因子RI計算公式為

    VIF可由下式計算:

    其中:

    式中:RRI表示河流水沙情勢變化系數(shù);RRIpre和RRIpost分別表示河流改變前后水沙情勢指數(shù);MRRPk表示河流的月情勢因子,計算過程參考文獻[20-21]。

    時間影響因子TIF主要考慮最大值、最小值和中值出現(xiàn)時間的變化,計算公式如下所示:

    式中:TF表示時間變量指數(shù);DTMax表示輸沙量和徑流量最大值出現(xiàn)時間改變量;DTMin表示輸沙量和徑流量最小值出現(xiàn)時間改變量;DTMedian為輸沙量和徑流量中值出現(xiàn)時間改變量??紤]到中值出現(xiàn)時間選取受水文數(shù)據(jù)統(tǒng)計方法限制,本次研究對式(7)進行精簡,去除原方法中DTMedian作為參數(shù)[22],即修改后時間變量指數(shù)的計算公式為

    公式(2)中,當RI處于0.75~1.00時為低度改變,0.50~0.75為中度改變,0.25~0.50為高度改變,0~0.25為嚴重改變[24]。

    1.6 累積量斜率變化率

    累積量斜率變化率法由王隨繼等[23]提出,可用于定量評估人類活動與氣候變化對河流泥沙改變的貢獻程度。該法假設:如果年輸沙量變化只受降水影響,則降水和輸沙量隨年份的累積曲線斜率為同倍比變化。將變量所有影響因素的綜合定義為1,根據(jù)各種影響因素隨時間累積斜率占變量累積斜率變化率的比值推求各種影響因素對變量的影響程度[24]。假定累積輸沙量隨時間變化在某一年出現(xiàn)拐點,拐點前后變量的斜率分別為YRa和YRb,累積降水量在拐點前后斜率分別為YPa和YPb。則:

    累積水沙量斜率變化率KR(%)為

    需要指出的是,氣候變化主要受到流域內(nèi)降水、蒸散發(fā)、氣溫以及地下水等因素影響。但蒸散發(fā)對徑流量、輸沙量影響較小,同時蒸散發(fā)數(shù)據(jù)獲取較為困難;而氣溫對于河流輸沙量的影響并不明顯[25];并且地下水變化主要體現(xiàn)在年內(nèi)分配,對年際變化影響可忽略不計,故本研究將氣候變化主要歸結為降水量的變化,通過分析降水與徑流量、輸沙量的關系來研究氣候變化對河流水沙情勢的影響。

    2 結果分析

    2.1 徑流-泥沙-降水變化與突變分析

    為揭示烏江流域多年徑流-泥沙-降水變化趨勢,繪制出烏江流域1965~2017年徑流量、輸沙量、降水量變化曲線(見圖1~3),由徑流量變化曲線、輸沙量變化曲線、降水量變化曲線可以看出近50 a烏江流域年徑流量、年輸沙量、年降水量均呈現(xiàn)出減少趨勢,其中徑流量和降水量表現(xiàn)出相似的減少趨勢且均不顯著,年輸沙量表現(xiàn)出顯著的減少趨勢。為定量分析烏江流域徑流量、輸沙量、降水量變化趨勢,采用Mann-Kendall趨勢檢驗計算烏江流域武隆水文站年徑流量、年輸沙量和流域內(nèi)8個氣象站點年均降水量變化趨勢,結果如表1所列。其中年徑流量、年輸沙量、年降水量Mann-Kendall統(tǒng)計量分別為-1.43,-6.95,-1.23,均小于零;年徑流量、年降水量統(tǒng)計量Zc的絕對值均小于顯著水平(取α=0.05)的臨界值1.96,表明烏江流域年徑流量和年降水量均為不顯著的下降趨勢;年輸沙量統(tǒng)計量Zc的絕對值大于顯著水平(取α=0.05)的臨界值1.96,表明烏江流域年輸沙量為顯著下降趨勢。

    由徑流量Mann-Kendall突變檢驗曲線(見圖4)可以看出:UF統(tǒng)計量和UB統(tǒng)計量的交點分別為1971,1986,1994,2003年。表明烏江流域徑流量可能在1971,1986,1994,2003年發(fā)生突變。

    由年輸沙量Mann-Kendall突變檢驗曲線(見圖5)可以看到:UF統(tǒng)計量和UB統(tǒng)計量的交點為1992年,表明年輸沙量在1992年可能發(fā)生突變。

    由降水量Mann-Kendall突變檢驗曲線圖(見圖6)可以看到:UF統(tǒng)計量和UB統(tǒng)計量的交點分別為1986,1994,2001年,說明烏江流域年降水量可能在1986,1994,2001年發(fā)生突變。

    同時采用滑動T檢驗對烏江流域多年水沙、降水數(shù)據(jù)進行突變檢驗,結果顯示烏江流域徑流量分別在1994,2000年發(fā)生突變,輸沙量分別在1975,1991,1995年發(fā)生突變,降水量分別在1990,2000年發(fā)生突變,即統(tǒng)計量t>tα(本文tα=t0.01=2.704)。繪制出烏江流域降水-徑流雙累積曲線、降水-輸沙量雙累積曲線(見圖7~8)。

    如圖7~8所示,烏江流域降水-徑流雙累積曲線未發(fā)生明顯偏轉,說明烏江流域降水徑流關系未發(fā)生改變。烏江流域降水-輸沙量雙累積曲線分別在1975,1992年發(fā)生了明顯的偏轉,說明烏江流域降水量-輸沙量關系在1975,1992年發(fā)生了突變。結合Mann-Kendall檢驗、滑動T檢驗結果和烏江流域降水-徑流雙累積曲線和降水-輸沙量雙累積曲線,最終確定烏江流域徑流量突變年份為1994年,輸沙量突變年份為1975,1992年,降水量突變年份為2000年(見表2)。

    根據(jù)突變年份對烏江年徑流量、年輸沙量變化趨勢進行劃分,將年徑流量劃分為2個時段,Ta基準期(1965~1994年),Tb變異期(1995~2017年);將年輸沙量劃分為3個階段,Ta基準期(1965~1975年),Tb變異Ⅰ期(1976~1992年),Tc變異Ⅱ期(1993~2017年)(見表3)。

    將本節(jié)結論與管曉祥[9]、郭文獻[4]、吳曉玲[5]和邵駿[3]等的部分研究結論進行比較(見表4)。結果發(fā)現(xiàn):由于管曉祥[9]和郭文獻[4]的研究站點相同且時間尺度相似,所以年徑流量與年輸沙量趨勢變化均呈現(xiàn)出下降趨勢,且年輸沙量下降幅度顯著,雖然本文所選取的時間尺度與二人的不同,但是年徑流量與年輸沙量均表現(xiàn)為相似的下降趨勢,且輸沙量下降趨勢明顯,但是徑流量、輸沙量、降水量分別在1994,1975,1992,2000年發(fā)生突變,突變點的不同可能是由于時間尺度選取的不同,導致Mann-Kendall與滑動T檢驗對于突變點的識別產(chǎn)生了變化。而吳曉玲[5]和邵駿[3]的研究站點不同且時間尺度也不相同,但是其所選站點均貫穿烏江上中下游,所以徑流量均表現(xiàn)相似的下降趨勢。這說明時間尺度選取的不同對于流域內(nèi)水文趨勢變化影響不大,但是對異常值的識別影響較為明顯。

    2.2 水沙相關性分析

    為了研究烏江流域徑流量與輸沙量相關性的階段性變化,結合武隆站徑流量與輸沙量變化的階段性劃分,將武隆站劃分為1965~1975年、1976~1992年、1993~2017年,3個階段。繪制武隆站輸沙量徑流量階段相關性曲線(見圖9)。根據(jù)圖9所示結果可以看出:1965~1975年、1976~1992年、1993~2017年期間武隆站年徑流量和年輸沙量在回歸線兩側較為分散,R值分別為0.34,0.01,0.22,呈現(xiàn)出隨時間變化逐漸減少的趨勢,說明武隆站水沙關系隨著時間的增長在逐漸變差,進一步說明了自20世紀90年代后烏江流域上修建運行的一系列水利工程設施對烏江水沙相關性變化的影響是十分顯著的。

    2.3 氣候變化和人類活動對水沙變化貢獻率

    根據(jù)上文徑流量、泥沙量變化階段劃分,采用累積量斜率變化率法分別計算不同時期人類活動和氣候變化對烏江流域徑流量、輸沙量變化的貢獻率,計算結果如表5所列。由貢獻率計算結果可知:在變異期Tb,氣候變化和人類活動對烏江流域徑流量變化的貢獻率分別為66.13%和33.87%,對輸沙量變化的貢獻率分別為11.62%和88.38%;在變異期Tc,氣候變化和人類活動對烏江輸沙量變化的貢獻率分別為6.61%和93.39%。研究結果表明:在變異期Tb,氣候變化是烏江流域徑流量發(fā)生年際變化的主要因素,對徑流量的影響占到了70%左右。在變異期Tb和Tc,人類活動是烏江流域輸沙量減少的主要影響因素,人類通過改變流域下墊面的方式影響流域輸沙是一個漸變的過程,畢竟大規(guī)模的開發(fā)和水土保持措施需要多年才能完成,其對流域產(chǎn)沙、輸沙的影響都是逐漸增加或減少的,而在河道上的水利工程建設卻能夠對輸沙產(chǎn)生立竿見影的影響[26]。在1984~2017年間,烏江上先后修建了烏江渡水電站(1982年)、普定水電站(1995年)、東風水電站(1995年)、引子渡水電站(2003年)、洪家渡水電站(2004年)、索風營水電站(2005年)、彭水水電站(2007年)、構皮灘水電站(2009年)、思林水電站(2009年)、銀盤水電站(2011年)、沙沱水電站(2013年),這些大型水利樞紐的建成運行是導致烏江流域輸沙量下降的主要原因。

    2.4 水沙改變度分析

    采用RI法(河流影響因子評價法)對烏江下游控制總站武隆站徑流量、輸沙量改變程度進行定量研究,根據(jù)表2中烏江水沙的階段性劃分結果來對徑流量、輸沙量階段性改變度進行分析計算,各RI指標計算結果見圖10和表6。結果表明:烏江流域徑流量MIF、VIF、TIF值在變異期Tb分別為0.98,0.44,0.36,說明烏江流域徑流量年內(nèi)最大值、最小值出現(xiàn)的時間和年內(nèi)分布分別發(fā)生了中度和高度改變,表明烏江干流水電站的修建只對烏江流域徑流量年內(nèi)變化影響較大,但是對于徑流量年際變化的影響較小;RI值為0.79,表明烏江流域徑流量在氣候變化的影響下其總體改變度為輕度改變。烏江流域輸沙量MIF、VIF、TIF值在輸沙量變異期Tb分別為0.54,0.45,0.45,表明從1976~1985年烏江年輸沙量年際變化較基準期Ta發(fā)生了高度改變,其年內(nèi)輸沙量最大值、最小值發(fā)生的時間和輸沙量年內(nèi)分布均發(fā)生了高度改變;RI值為0.48,表明在變異期Tb(1976~1992年)烏江流域水電站的修建和水土保持治理措施等人類活動使得烏江年輸沙量總體改變度為高度改變。而在變異期Tc,MIF、VIF、TIF、RI值分別為0.14,0.44,0.40,0.12,表明在1993~2017年輸沙量年際變化較基準期發(fā)生了嚴重改變,輸沙量年內(nèi)變化和輸沙量最大值、最小值出現(xiàn)的時間發(fā)生了高度改變,綜合改變度較基準期為嚴重改變??偠灾?,烏江流域年徑流量在受氣候變化影響方面其綜合改變度為輕度改變,年輸沙量在人類活動的影響下已經(jīng)發(fā)生了嚴重改變。

    3 徑流泥沙演變影響因素

    影響流域內(nèi)徑流量和輸沙量的因素主要包括流域內(nèi)下墊面條件自然災害、氣候降水、人類影響等[27],其中人類活動和降水對河流的產(chǎn)流產(chǎn)沙影響較大。根據(jù)近年的《貴州水資源公報》水資源量統(tǒng)計分析發(fā)現(xiàn):研究區(qū)工農(nóng)業(yè)等取水量未發(fā)生顯著變化[5] 。因此人類活動,例如修建水利工程等對河川徑流進行了一定的調(diào)節(jié),從而改變了河川徑流原有的時空分布。由于烏江上興建的水庫多為日調(diào)節(jié)或不完全年調(diào)節(jié)水庫,對年徑流量的影響相對較小,主要改變的是徑流的年內(nèi)分配過程。而降水變化對徑流的影響較為敏感,能夠直接影響徑流量年際變化[3]。

    人類活動對烏江流域年輸沙量的影響主要表現(xiàn)在水利工程的修建、水土保持治理措施的實施。為了評價水庫蓄水攔沙對烏江年輸沙量的影響,將武隆站作為總控制站分析烏江流域大型水庫建設情況與控制總站武隆站年輸沙量年際變化之間的關系,如表7所列。根據(jù)表7可知年輸沙量的階段性減少與流域內(nèi)水庫累積庫容的增加具有良好的相關性。在1965~1975年間,烏江上修建了一些小型的水電站,但是由于其庫容量較小其蓄水攔沙作用不明顯,所以在20世紀60~80年代初,烏江流域輸沙量變化不明顯,而隨著1982年烏江渡電站(庫容23.0億m3)的修建,1980~1985年庫區(qū)淤積泥沙約1.20億m3(1.39億t),年均淤積量為2 000萬m3(2 310萬t)[28],使得變異期Tb的多年年均輸沙量1.94×107 t較基準期3.61×107 t減少了1.67×107 t,約占基準期年均輸沙量的50%。隨著1991~2017年烏江干流東風水電站(1994年蓄水,總庫容10.25億m3)、普定水電站(1995年建成,總庫容4.20億m3)、引子渡水電站(2003年蓄水,總庫容5.27億m3)、洪家渡水電站(2004年蓄水,總庫容49.47億m3)、烏江干流索風營水電站(2005年蓄水,庫容1.67億m3)、2007年彭水水電站(2007年蓄水,庫容24.65億m3)、思林水電站(2009年蓄水,庫容為12.05億m3)、構皮灘水電站(2009年蓄水,庫容為55.64億m3)等的建成蓄水,使得烏江流域年輸沙量在變異期Tc的多年年均輸沙量0.51×107 t較基準期3.61×107 t減少了3.10×107 t,約占基準期年均輸沙量的90%。所以烏江流域水利樞紐設施的興建與運行對于烏江干流年輸沙量的影響是巨大的。

    4 結 論

    (1)通過對烏江流域水沙、降水量的趨勢分析,發(fā)現(xiàn)武隆站徑流量、降水量其Mann-Kendall標準化變量分別為-1.43,-1.23,均未通過95%置信檢驗,表明徑流量、降水量表現(xiàn)出不明顯的下降趨勢;輸沙量Mann-Kendall標準化變量為-6.95,通過了95%的置信區(qū)間,表明輸沙量下降幅度顯著。

    (2)根據(jù)Mann-Kendall非參數(shù)檢驗法和滑動T檢驗法對武隆站多年水沙序列和流域多年降水序列進行突變檢驗并結合徑流量-降水量雙累積曲、輸沙量-降水量雙累積曲線發(fā)現(xiàn)武隆站年徑流量在1994年發(fā)生突變,輸沙量分別在1975,1992年發(fā)生突變,降水量在2000年發(fā)生突變。并將徑流量變化劃分為1965~1994年、1995~2017年兩個階段,將輸沙量變化劃分為1965~1975年、1976~1992年、1993~2017年3個階段。

    (3)根據(jù)對烏江流域泥沙序列的階段性劃分對烏江流域水沙進行階段相關性分析,結果顯示在輸沙量基準期Ta、變異期Tb、變異期Tc其相關系數(shù)R分別為0.34,0.01,0.22,表明烏江水沙相關性自變異期Tb后均處于極差的狀態(tài)。

    (4)通過累積量斜率變化率法對烏江流域徑流量、輸沙量在不同變異期內(nèi)氣候變化和人類活動的貢獻率進行計算,計算結果顯示在徑流量變異期Tb(1995~2017年)氣候變化和人類活動的貢獻率分別為66.13%,33.87%,說明在此期間氣候變化是年徑流量減少的主要驅動力,輸沙量變異期Tb(1976~1992年)、變異期Tc(1993~2017年)氣候變化和人類活動的貢獻率分別為11.62%和88.38%,6.61%和93.39%,表明烏江流域年輸沙量減少的主要驅動因素是人類活動。

    (5)采用RI法對烏江水沙不同階段的水文改變度進行計算分析,結果顯示年徑流量在變異期Tb(1995~2017年)在年際變化上為輕度改變,年輸沙量在變異期Tb(1976~1992年)、變異期Tc(1993~2017年)分別為高度、嚴重改變。

    (6)通過對烏江流域水庫蓄水攔沙對年輸沙量影響分析可知:年輸沙量在變異期Tb(1976~1992年)較基準期Ta(1965~1975年)減少了1.67×107 t,約占基準期年均輸沙量的50%;在變異期Tc(1993~2017年)較基準期Ta(1965~1975年)減少了3.10×107 t,約占基準期年均輸沙量的90%。說明水庫蓄水攔沙是導致輸沙量顯著減少的重要原因之一。

    參考文獻:

    [1] 張冬青,何騰兵,林昌虎,等.貴州烏江流域喀斯特土壤物理性質的研究[J].貴州科學,2006,24(2):39-43.

    [2] 韋紅波,任紅玉,楊勤科.中國多年平均輸沙模數(shù)的研究[J].泥沙研究,2003(1):39-44.

    [3] 邵駿,范可旭,邴建平,等.烏江干流年徑流變化趨勢及成因分析[J].水文,2012,32(6):86-91.

    [4] 郭文獻,李越,查胡飛,等.近50年來烏江水沙特征變化研究[J].水力發(fā)電,2020,46(2):22-26.

    [5] 吳曉玲,張欣,向小華,等.烏江流域上游水沙特性變化及其水電站建設影響[J].生態(tài)學雜志,2018,37(3):642-650.

    [6] 陳松生,許全喜,陳澤方.烏江流域水沙變化特性及其原因分析[J].泥沙研究,2008(5):43-48.

    [7] 班璇,朱碧瑩,舒鵬,等.漢江流域氣象水文變化趨勢及驅動力分析[J].長江流域資源與環(huán)境,2018,27(12):2817-2829.

    [8] 熊亞蘭,張科利,楊光檄,等.烏江流域水沙特性變化分析[J].生態(tài)環(huán)境,2008,17(5):1942-1947.

    [9] 管曉祥,向小華,李超,等.烏江流域水沙規(guī)律演變及驅動因素定量分析[J].泥沙研究,2019,44(5):36-41.

    [10] 張曉婭.近60年氣候變化和人類活動對長江徑流量影響的研究[D].上海:華東師范大學,2014.

    [11] 李晶瑩,張經(jīng).中國主要河流的輸沙量及其影響因素[J].青島海洋大學學報(自然科學版),2003,33(4):565-573.

    [12] KAHYA E,KALAYC S.Trend analysis of streamflow in Turkey [J].Journal of Hydrology,2004,289:128-144.

    [13] BURN D H,ELNUR,MAH.Detection of hydrologic trends and variability[J].Journal of Hydrology,2002,255:107-122.

    [14] 魏鳳英,曹鴻興.中國、北半球和全球的氣溫突變分析及其趨勢預測研究[J].大氣科學,1995,19(2):140-148.

    [15] 于家瑞,艾萍,袁定波,等.雅礱江流域徑流演變特征及成因[J].水電能源科學,2019,37(10):13-16.

    [16] 丁文峰,張平倉,任洪玉.近50年來嘉陵江流域徑流泥沙演變規(guī)律及驅動因素定量分析[J].長江科學院院報,2008,25(3):23-27.

    [17] 武旭同,李娜,王臘春.近60年來長江干流水沙特征分析[J].泥沙研究,2016(5):40-46.

    [18] 趙玉,穆興民,何毅,等.1950~2011年黃河干流水沙關系變化研究[J].泥沙研究,2014(4):32-38.

    [19] 饒素秋,霍世青,薛建國,等.黃河上中游水沙變化特點分析及未來趨勢展望[J].泥沙研究,2001(2):74-77.

    [20] HAGHIGHI A T,KLVE B .Development of a general river regime index (RRI)for intra-annual flow variation based on the unit river concept and flow variation end-points[J].Journal of Hydrology,2013,503:169-177.

    [21] 劉靜玲,尤曉光,史璇,等.灤河流域大中型閘壩水文生態(tài)效應[J].水資源保護,2016,32(1):23-28.

    [22] HAGHIGHI A T,MARTTILA H,KLVE B.Development of a new index to assess river regime impacts after dam construction[J].Global and Planetary Change,2014,122:186-196.

    [23] 王隨繼,李玲,顏明.氣候和人類活動對黃河中游區(qū)間產(chǎn)流量變化的貢獻率[J].地理研究,2013,32(3):395-402.

    [24] 王彥君,王隨繼,蘇騰.降水和人類活動對松花江徑流量變化的貢獻率[J].自然資源學報,2015,30(2):304-314.

    [25] 田清.近60年來氣候變化和人類活動對黃河、長江、珠江水沙通量影響的研究[D].上海:華東師范大學,2016.

    [26] FRIEDMAN,OSTERKAMP W R,SCOTT M L,et al.Downstream effects of dams on channel geometry and bottomland vegetation:regional patterns in the great plains[J].Wetlands,1998,18(4):619-633.

    [27] 王延貴,史紅玲,劉茜.水庫攔沙對長江水沙態(tài)勢變化的影響[J].水科學進展,2014,25(4):467-476.

    [28] 劉家應,柴方福.烏江渡水庫泥沙淤積觀測研究(二)[J].貴州水力發(fā)電,1994(1):33-43.

    (編輯:江 文)

    猜你喜歡
    影響
    美食網(wǎng)紅如何影響我們吃什么
    英語文摘(2022年4期)2022-06-05 07:45:18
    是什么影響了滑動摩擦力的大小
    哪些顧慮影響擔當?
    當代陜西(2021年2期)2021-03-29 07:41:24
    影響大師
    沒錯,痛經(jīng)有時也會影響懷孕
    媽媽寶寶(2017年3期)2017-02-21 01:22:28
    擴鏈劑聯(lián)用對PETG擴鏈反應與流變性能的影響
    中國塑料(2016年3期)2016-06-15 20:30:00
    基于Simulink的跟蹤干擾對跳頻通信的影響
    如何影響他人
    APRIL siRNA對SW480裸鼠移植瘤的影響
    日韩强制内射视频| 久久久色成人| 人妻夜夜爽99麻豆av| 欧美激情久久久久久爽电影| 亚洲欧美中文字幕日韩二区| 亚洲精品自拍成人| 国产 一区 欧美 日韩| av卡一久久| 看非洲黑人一级黄片| 国产高清不卡午夜福利| 黄片wwwwww| 噜噜噜噜噜久久久久久91| 大香蕉97超碰在线| 又黄又爽又刺激的免费视频.| 国产一区二区三区综合在线观看 | 99热这里只有是精品在线观看| 极品教师在线视频| 亚洲最大成人中文| 亚洲欧美日韩无卡精品| 超碰av人人做人人爽久久| 亚洲欧美中文字幕日韩二区| 啦啦啦韩国在线观看视频| 看免费成人av毛片| 自拍偷自拍亚洲精品老妇| 麻豆国产97在线/欧美| 久久国内精品自在自线图片| 国产成人a区在线观看| 久久久久久久久久人人人人人人| 精品国产三级普通话版| 欧美成人一区二区免费高清观看| 久久久久久久久久久丰满| 国产伦在线观看视频一区| ponron亚洲| 国产精品久久久久久精品电影| 日韩成人av中文字幕在线观看| 51国产日韩欧美| 精品久久久噜噜| 夜夜看夜夜爽夜夜摸| 精品久久久久久久久亚洲| 亚洲av电影在线观看一区二区三区 | 中文字幕av成人在线电影| 国产精品人妻久久久影院| 一区二区三区免费毛片| 亚洲天堂国产精品一区在线| 国产成人精品久久久久久| av免费在线看不卡| 亚洲精品国产av成人精品| 久久精品国产亚洲av天美| 欧美激情久久久久久爽电影| 亚洲av成人精品一二三区| 嫩草影院新地址| 能在线免费观看的黄片| 午夜视频国产福利| 街头女战士在线观看网站| 一级毛片电影观看| 亚洲欧美一区二区三区国产| 免费av观看视频| 亚洲精品视频女| 色哟哟·www| 麻豆成人午夜福利视频| 国产在线男女| 2021天堂中文幕一二区在线观| 久久久精品欧美日韩精品| 国产精品一区www在线观看| 午夜精品一区二区三区免费看| 亚洲国产高清在线一区二区三| 国产成人免费观看mmmm| 亚洲精品国产av成人精品| 国产亚洲av片在线观看秒播厂 | 激情五月婷婷亚洲| 中文乱码字字幕精品一区二区三区 | 亚洲精品,欧美精品| 高清日韩中文字幕在线| 少妇的逼好多水| 日韩成人av中文字幕在线观看| 能在线免费观看的黄片| 久久久a久久爽久久v久久| 十八禁网站网址无遮挡 | 在线观看av片永久免费下载| 日韩制服骚丝袜av| 一区二区三区四区激情视频| av女优亚洲男人天堂| 亚洲国产成人一精品久久久| 成年女人在线观看亚洲视频 | 天堂俺去俺来也www色官网 | 国产在视频线在精品| 亚洲三级黄色毛片| 中文字幕人妻熟人妻熟丝袜美| 国产在线一区二区三区精| 在线观看人妻少妇| 欧美变态另类bdsm刘玥| 熟妇人妻久久中文字幕3abv| 波野结衣二区三区在线| 国产精品一区二区性色av| 中文字幕av成人在线电影| 亚洲熟女精品中文字幕| 草草在线视频免费看| 久久久久久伊人网av| 激情五月婷婷亚洲| 一级毛片我不卡| 最近2019中文字幕mv第一页| 一级爰片在线观看| freevideosex欧美| 日韩av在线大香蕉| 老女人水多毛片| 免费观看在线日韩| 狂野欧美激情性xxxx在线观看| 人妻系列 视频| 久久精品国产亚洲av涩爱| 亚洲成人中文字幕在线播放| 一级毛片我不卡| 国产成人a区在线观看| 一区二区三区高清视频在线| 有码 亚洲区| 少妇人妻精品综合一区二区| 岛国毛片在线播放| 国产 亚洲一区二区三区 | 成人午夜高清在线视频| 久久久a久久爽久久v久久| 两个人的视频大全免费| 综合色丁香网| 中文字幕人妻熟人妻熟丝袜美| 午夜福利在线观看免费完整高清在| 亚洲av成人av| 麻豆av噜噜一区二区三区| 亚洲精华国产精华液的使用体验| 日韩伦理黄色片| 两个人视频免费观看高清| 丝瓜视频免费看黄片| 夜夜看夜夜爽夜夜摸| 2018国产大陆天天弄谢| a级毛色黄片| 成人高潮视频无遮挡免费网站| 搡女人真爽免费视频火全软件| 一区二区三区乱码不卡18| 一级爰片在线观看| 欧美日韩在线观看h| 一个人观看的视频www高清免费观看| 高清在线视频一区二区三区| 男女啪啪激烈高潮av片| 卡戴珊不雅视频在线播放| 久久久久久久国产电影| 国产高潮美女av| 欧美性猛交╳xxx乱大交人| 国产精品蜜桃在线观看| 亚洲第一区二区三区不卡| 色哟哟·www| 天堂俺去俺来也www色官网 | 联通29元200g的流量卡| 日产精品乱码卡一卡2卡三| 久久97久久精品| av线在线观看网站| 麻豆久久精品国产亚洲av| 亚洲最大成人av| 欧美成人精品欧美一级黄| 亚洲av二区三区四区| 国产精品一区二区三区四区久久| 蜜桃亚洲精品一区二区三区| 肉色欧美久久久久久久蜜桃 | av国产久精品久网站免费入址| 亚洲av成人精品一区久久| 亚洲精品自拍成人| 午夜老司机福利剧场| 一级毛片aaaaaa免费看小| 欧美日韩在线观看h| 97热精品久久久久久| 啦啦啦韩国在线观看视频| 黄片无遮挡物在线观看| 亚洲综合精品二区| 91狼人影院| 99久久精品热视频| 亚洲精品第二区| 国产一区二区三区综合在线观看 | 国产不卡一卡二| 亚洲精品国产av成人精品| 精品国内亚洲2022精品成人| 亚洲电影在线观看av| 哪个播放器可以免费观看大片| 一边亲一边摸免费视频| 最近视频中文字幕2019在线8| 一区二区三区乱码不卡18| 色视频www国产| 精品亚洲乱码少妇综合久久| 国产在线男女| 美女国产视频在线观看| 美女高潮的动态| 色网站视频免费| 高清日韩中文字幕在线| 国模一区二区三区四区视频| 99久国产av精品| 日本免费在线观看一区| 69人妻影院| 一级爰片在线观看| 国内揄拍国产精品人妻在线| 99热6这里只有精品| 亚洲成色77777| 亚洲国产欧美人成| av线在线观看网站| 91aial.com中文字幕在线观看| 能在线免费看毛片的网站| 久久久精品免费免费高清| 婷婷色综合大香蕉| 久久人人爽人人片av| 亚洲最大成人中文| 三级男女做爰猛烈吃奶摸视频| 亚洲av日韩在线播放| 欧美bdsm另类| 狂野欧美激情性xxxx在线观看| 久久久久网色| 成人亚洲精品av一区二区| 日韩欧美一区视频在线观看 | 亚洲精品日韩av片在线观看| 亚洲最大成人av| 久久久久网色| 十八禁国产超污无遮挡网站| 女人久久www免费人成看片| 汤姆久久久久久久影院中文字幕 | 99热6这里只有精品| 国产欧美日韩精品一区二区| 男女边摸边吃奶| 国产黄色小视频在线观看| 日韩欧美国产在线观看| 国产 一区精品| 国产亚洲精品久久久com| 伊人久久精品亚洲午夜| 国产亚洲最大av| 尤物成人国产欧美一区二区三区| 美女大奶头视频| 好男人在线观看高清免费视频| 日本熟妇午夜| 听说在线观看完整版免费高清| 国产亚洲av片在线观看秒播厂 | 国产精品一区www在线观看| www.av在线官网国产| 亚洲欧洲日产国产| 欧美成人午夜免费资源| 亚洲国产精品sss在线观看| 久久久色成人| 插逼视频在线观看| 搡老妇女老女人老熟妇| 18+在线观看网站| 亚洲一级一片aⅴ在线观看| 亚洲精品中文字幕在线视频 | 美女xxoo啪啪120秒动态图| 国产综合懂色| 日韩不卡一区二区三区视频在线| 青春草国产在线视频| 高清午夜精品一区二区三区| 国产成年人精品一区二区| 嘟嘟电影网在线观看| 1000部很黄的大片| 插阴视频在线观看视频| 网址你懂的国产日韩在线| 国产亚洲午夜精品一区二区久久 | 国产成人freesex在线| 国产成人精品婷婷| 狠狠精品人妻久久久久久综合| 日韩人妻高清精品专区| av国产久精品久网站免费入址| 日日摸夜夜添夜夜添av毛片| 精品久久久久久久久av| 国产成人a区在线观看| 婷婷色综合大香蕉| 天堂俺去俺来也www色官网 | 在线观看免费高清a一片| 赤兔流量卡办理| 亚洲av男天堂| 内地一区二区视频在线| 五月伊人婷婷丁香| 免费看av在线观看网站| 日韩av在线免费看完整版不卡| 国语对白做爰xxxⅹ性视频网站| 伦精品一区二区三区| 亚洲内射少妇av| 亚洲av在线观看美女高潮| 只有这里有精品99| 国产中年淑女户外野战色| 久久久久久久久久人人人人人人| 日本免费a在线| 国产伦精品一区二区三区四那| 久久99热这里只有精品18| 男人狂女人下面高潮的视频| 天天躁夜夜躁狠狠久久av| 一本久久精品| 成人亚洲欧美一区二区av| 国产成人精品久久久久久| 日日干狠狠操夜夜爽| 国产精品人妻久久久久久| 免费电影在线观看免费观看| 能在线免费看毛片的网站| 国产伦精品一区二区三区四那| 午夜福利在线在线| 一本久久精品| 亚洲欧美成人精品一区二区| 国产免费一级a男人的天堂| 日本猛色少妇xxxxx猛交久久| 别揉我奶头 嗯啊视频| 中文欧美无线码| 乱系列少妇在线播放| 你懂的网址亚洲精品在线观看| 国产高清有码在线观看视频| 亚洲在久久综合| 亚洲欧洲日产国产| 欧美潮喷喷水| av播播在线观看一区| 国模一区二区三区四区视频| 亚洲精品一二三| 2018国产大陆天天弄谢| 亚洲三级黄色毛片| 久久久国产一区二区| 高清在线视频一区二区三区| 婷婷色综合大香蕉| 国产精品蜜桃在线观看| 哪个播放器可以免费观看大片| 少妇熟女欧美另类| 久久久亚洲精品成人影院| 免费观看的影片在线观看| h日本视频在线播放| 日本av手机在线免费观看| 午夜免费观看性视频| 亚洲精品色激情综合| 精品亚洲乱码少妇综合久久| 男人舔女人下体高潮全视频| 晚上一个人看的免费电影| av在线观看视频网站免费| 人妻少妇偷人精品九色| 国产女主播在线喷水免费视频网站 | 成年女人在线观看亚洲视频 | 亚洲最大成人手机在线| 黄色日韩在线| 美女内射精品一级片tv| 搞女人的毛片| 天美传媒精品一区二区| 久久久久免费精品人妻一区二区| 中国国产av一级| 亚洲av中文字字幕乱码综合| 亚洲av福利一区| av在线老鸭窝| 国产男人的电影天堂91| 青春草视频在线免费观看| 亚洲美女视频黄频| 国产在线一区二区三区精| 日本午夜av视频| 免费看av在线观看网站| 国产真实伦视频高清在线观看| 国产精品一区www在线观看| 久久久欧美国产精品| 看非洲黑人一级黄片| 久久精品久久久久久久性| 日韩av在线免费看完整版不卡| 国产淫片久久久久久久久| 亚洲第一区二区三区不卡| 男插女下体视频免费在线播放| av天堂中文字幕网| 能在线免费看毛片的网站| 国产黄片美女视频| 最后的刺客免费高清国语| 成年版毛片免费区| 六月丁香七月| 国产有黄有色有爽视频| 男插女下体视频免费在线播放| 精品国产一区二区三区久久久樱花 | 99久久九九国产精品国产免费| 校园人妻丝袜中文字幕| 日韩,欧美,国产一区二区三区| 美女xxoo啪啪120秒动态图| 久久亚洲国产成人精品v| 淫秽高清视频在线观看| 中文资源天堂在线| 日韩成人av中文字幕在线观看| 麻豆成人av视频| 免费看美女性在线毛片视频| 成人高潮视频无遮挡免费网站| 日本爱情动作片www.在线观看| 一个人看视频在线观看www免费| 亚洲av电影在线观看一区二区三区 | 国产精品99久久久久久久久| 国精品久久久久久国模美| 亚洲欧美日韩东京热| 最近中文字幕高清免费大全6| 伦精品一区二区三区| 少妇裸体淫交视频免费看高清| 麻豆成人av视频| 观看美女的网站| 内射极品少妇av片p| 一个人免费在线观看电影| 高清视频免费观看一区二区 | 国产毛片a区久久久久| 一二三四中文在线观看免费高清| 白带黄色成豆腐渣| 一级毛片我不卡| 人体艺术视频欧美日本| 国产亚洲精品久久久com| 国产一级毛片七仙女欲春2| 日本一二三区视频观看| 欧美三级亚洲精品| 成人性生交大片免费视频hd| 婷婷色麻豆天堂久久| 韩国高清视频一区二区三区| 亚洲熟女精品中文字幕| 亚洲不卡免费看| 亚洲av成人精品一二三区| 又爽又黄a免费视频| 最近最新中文字幕免费大全7| 老女人水多毛片| 精品午夜福利在线看| 日韩电影二区| 亚洲国产欧美人成| av在线亚洲专区| 一夜夜www| 国产精品一区二区三区四区免费观看| 久久午夜福利片| 国产有黄有色有爽视频| 国产精品久久久久久精品电影小说 | 国产精品麻豆人妻色哟哟久久 | 亚洲,欧美,日韩| 久久99热这里只频精品6学生| 少妇熟女aⅴ在线视频| 一个人观看的视频www高清免费观看| 国产精品一区二区在线观看99 | 麻豆成人av视频| 国产黄片视频在线免费观看| 欧美最新免费一区二区三区| 老司机影院成人| 精品人妻偷拍中文字幕| 麻豆精品久久久久久蜜桃| av免费观看日本| 精品久久久久久成人av| av免费在线看不卡| 国产伦精品一区二区三区四那| or卡值多少钱| 国内精品一区二区在线观看| 午夜福利在线观看吧| 亚洲电影在线观看av| 久久久久网色| 久久精品人妻少妇| 91精品国产九色| 国产精品av视频在线免费观看| 精品人妻一区二区三区麻豆| 国产精品一区二区性色av| 国产精品嫩草影院av在线观看| 国产精品一区二区性色av| 在线免费观看不下载黄p国产| av播播在线观看一区| 中文字幕制服av| 免费av毛片视频| 天堂网av新在线| 色综合站精品国产| 中国美白少妇内射xxxbb| 久久6这里有精品| 精品不卡国产一区二区三区| 亚洲精品日本国产第一区| 久久热精品热| 国产单亲对白刺激| 菩萨蛮人人尽说江南好唐韦庄| 99热这里只有是精品50| 91久久精品国产一区二区三区| 又粗又硬又长又爽又黄的视频| 国产成人午夜福利电影在线观看| 日日干狠狠操夜夜爽| 大陆偷拍与自拍| 五月天丁香电影| 日产精品乱码卡一卡2卡三| 欧美成人精品欧美一级黄| 色综合亚洲欧美另类图片| 国产大屁股一区二区在线视频| 亚洲熟妇中文字幕五十中出| 国产精品人妻久久久久久| 91在线精品国自产拍蜜月| 大陆偷拍与自拍| 国产精品久久久久久精品电影小说 | 特大巨黑吊av在线直播| 久久精品久久久久久噜噜老黄| 中文在线观看免费www的网站| 天堂√8在线中文| 国产午夜精品一二区理论片| 一个人免费在线观看电影| 黄色日韩在线| 老司机影院成人| 蜜臀久久99精品久久宅男| 午夜福利成人在线免费观看| 亚洲精品色激情综合| 日韩一区二区视频免费看| 偷拍熟女少妇极品色| 人人妻人人澡人人爽人人夜夜 | 一个人看视频在线观看www免费| 美女内射精品一级片tv| 国产精品一区二区性色av| 非洲黑人性xxxx精品又粗又长| 美女主播在线视频| 国产高清三级在线| 激情 狠狠 欧美| 色综合色国产| 男人和女人高潮做爰伦理| 中文天堂在线官网| 久久草成人影院| 久久久久久久久大av| a级一级毛片免费在线观看| 亚洲av一区综合| 99视频精品全部免费 在线| 国产黄色视频一区二区在线观看| 97人妻精品一区二区三区麻豆| 最近最新中文字幕大全电影3| 丰满乱子伦码专区| 内射极品少妇av片p| 亚洲综合精品二区| 国产黄色免费在线视频| 内地一区二区视频在线| 亚洲内射少妇av| 天天一区二区日本电影三级| 禁无遮挡网站| 黄片无遮挡物在线观看| 日韩欧美国产在线观看| 亚洲最大成人av| 看十八女毛片水多多多| 乱系列少妇在线播放| 美女主播在线视频| 免费观看精品视频网站| 91在线精品国自产拍蜜月| 天堂中文最新版在线下载 | 国产精品久久久久久久久免| 97超视频在线观看视频| 免费看不卡的av| av卡一久久| 国产乱人视频| 免费黄网站久久成人精品| 久久久国产一区二区| 日韩,欧美,国产一区二区三区| 大香蕉97超碰在线| 91av网一区二区| 久久久久久久亚洲中文字幕| 18禁在线播放成人免费| 精品国产露脸久久av麻豆 | 亚洲国产高清在线一区二区三| 麻豆乱淫一区二区| 精品国产三级普通话版| 久久久精品免费免费高清| 91精品伊人久久大香线蕉| 成人毛片a级毛片在线播放| 欧美日韩亚洲高清精品| 精品久久久久久久久av| 色综合色国产| 久久精品夜色国产| 精品久久久久久电影网| 久久久久久久久中文| 天天躁日日操中文字幕| 22中文网久久字幕| 色尼玛亚洲综合影院| 综合色av麻豆| 美女脱内裤让男人舔精品视频| a级毛色黄片| 99久久中文字幕三级久久日本| 在线播放无遮挡| 国产在视频线精品| 91aial.com中文字幕在线观看| 久久久a久久爽久久v久久| 综合色丁香网| av在线蜜桃| 中文字幕免费在线视频6| 国产亚洲午夜精品一区二区久久 | 黄片wwwwww| 午夜免费观看性视频| 一个人观看的视频www高清免费观看| 国产爱豆传媒在线观看| 亚洲美女搞黄在线观看| 精品一区在线观看国产| 精品国产露脸久久av麻豆 | 国产女主播在线喷水免费视频网站 | 纵有疾风起免费观看全集完整版 | 男女啪啪激烈高潮av片| 国产在线一区二区三区精| 舔av片在线| av.在线天堂| 蜜桃亚洲精品一区二区三区| 有码 亚洲区| 国产精品1区2区在线观看.| 51国产日韩欧美| 亚洲av.av天堂| 色5月婷婷丁香| 国产精品嫩草影院av在线观看| 少妇被粗大猛烈的视频| 日韩欧美一区视频在线观看 | 日本黄大片高清| 午夜精品国产一区二区电影 | 一区二区三区四区激情视频| av网站免费在线观看视频 | 亚洲欧美日韩无卡精品| 午夜精品一区二区三区免费看| 久久6这里有精品| 日韩精品青青久久久久久| 国产大屁股一区二区在线视频| 日本av手机在线免费观看| 午夜福利网站1000一区二区三区| 丝瓜视频免费看黄片| 一级二级三级毛片免费看| 国产午夜福利久久久久久| 国产91av在线免费观看| 噜噜噜噜噜久久久久久91| 国产一级毛片七仙女欲春2| 狠狠精品人妻久久久久久综合| 国国产精品蜜臀av免费| 欧美日韩综合久久久久久| 婷婷色av中文字幕| 夫妻午夜视频| 99视频精品全部免费 在线| 99久国产av精品| 久久久精品94久久精品| 99视频精品全部免费 在线| 日韩一区二区三区影片| xxx大片免费视频| 国内精品一区二区在线观看| 亚洲人与动物交配视频| 在线免费观看的www视频| 日本熟妇午夜| 久久精品夜色国产|