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

    1970~2022年金沙江上中游河段徑流演變規(guī)律研究

    2024-10-25 00:00:00汪琳周波郭衛(wèi)朱迪卜慧左建
    人民長江 2024年10期

    摘要:金沙江上中游河段是中國規(guī)劃建設(shè)的重要水電能源基地,探究變化環(huán)境下金沙江上中游河川徑流演變規(guī)律對于促進(jìn)流域水資源高效利用與管理具有重要意義。采用非參數(shù)檢驗、滑動t檢驗法、小波功率法系統(tǒng)研究了金沙江上中游河段的徑流演變規(guī)律,并結(jié)合不均勻系數(shù)、完全調(diào)節(jié)系數(shù)、集中度和基尼系數(shù)等指標(biāo)綜合分析了梯級水庫建設(shè)運(yùn)行對金沙江中游河段徑流的影響。結(jié)果表明:1970~2022年間,金沙江上中游流域的年均流量呈現(xiàn)出不顯著上升趨勢;流域徑流從上游到中游增加幅度越來越大,中游河段的多年平均流量相較上游增加了24.5%,而徑流的突變點(diǎn)越來越滯后,且從上游到中游的徑流變化周期顯著性降低;梯級水庫的投運(yùn)與水庫的蓄洪補(bǔ)枯調(diào)節(jié)作用不僅使得中游河段來水量更加均勻,波動性減小,且河段徑流在豐水期水量下降,而枯水期水量增加,最小月均流量出現(xiàn)時間提前。 研究成果可為金沙江流域水資源優(yōu)化配置提供理論支撐。

    關(guān) 鍵 詞:河川徑流;演變規(guī)律;水庫運(yùn)行;金沙江

    中圖法分類號:TV121.2

    文獻(xiàn)標(biāo)志碼:ADOI:10.16232/j.cnki.1001-4179.2024.10.016

    0 引 言

    隨著氣候變化和人類活動影響的加劇,河川徑流的特征發(fā)生了一定改變,掌握河川徑流演變規(guī)律對流域水資源管理與高效利用具有重要價值[1-3。金沙江上中游河段河道狹窄,水流湍急,水電開發(fā)條件優(yōu)越,該河段承擔(dān)滇中調(diào)水、兩岸灌溉用水和遠(yuǎn)景“南水北調(diào)”西線調(diào)水任務(wù),是中國規(guī)劃建設(shè)的重要水電基地4。因此,研究金沙江上中游河川徑流的演變規(guī)律是金中河段水電開發(fā)的重要基礎(chǔ)[5。

    近年來,許多學(xué)者對各大流域的河川徑流演變規(guī)律展開了大量研究。王國慶等[6針對1956~2018年間中國主要江河的徑流變化特征進(jìn)行了分析,研究發(fā)現(xiàn)除了長江流域的大通站之外,中國十大流域的代表性水文站年徑流量都呈現(xiàn)下降趨勢。舒章康等[7分析了1961~2018年中國主要江河的枯季徑流演變特征,并剖析了不同河流徑流演變的原因,研究發(fā)現(xiàn)長江流域枯季徑流增加十分顯著,而降水增加是其變化的主要原因。鮑振鑫等[8采用Mann-Kendall法探討了黃河上游、中游、下游流域的不同尺度徑流演變規(guī)律,結(jié)果表明除了黃河源區(qū)的徑流變化幅度較小之外,黃河流域的年尺度、季尺度徑流均呈現(xiàn)顯著的下降趨勢??祦嗢o等[9以寸灘站為代表站研究了長江上游的徑流演變特征,結(jié)果表明長江上游的年均徑流呈現(xiàn)不顯著減小的趨勢。

    金沙江流域作為中國重要的水電能源基地,其河川徑流演變規(guī)律也一直是學(xué)者們關(guān)心的熱點(diǎn)問題。馮勝航等[10采用水量平衡法還原了金沙江流域的天然徑流,并探究了金沙江流域天然徑流量的變化特征及其成因,研究發(fā)現(xiàn)氣候變化是金沙江流域徑流變化的主要影響因素。張小峰等[11系統(tǒng)研究了金沙江屏山站以上區(qū)域的年徑流量變化,結(jié)果表明降水是金沙江流域徑流變化的主要驅(qū)動因子,且石鼓站以上區(qū)段受降水影響最為明顯。郭衛(wèi)等12采用RVA法分析了金沙江下游的水文情勢變化情況,研究表明現(xiàn)有梯級水庫的運(yùn)行會顯著改變金沙江下游的來水,尤其是枯季來水。

    然而,目前針對金沙江上中游流域徑流演變規(guī)律的研究大部分聚焦在氣候變化對徑流的影響方面,系統(tǒng)探討梯級水庫運(yùn)行對金沙江徑流的影響研究相對較少。因此,本研究擬從趨勢性、突變性、周期性等多個方面分析金沙江上中游河段的河川徑流演變特征,剖析徑流演變規(guī)律,并結(jié)合多個指標(biāo)從年際、年內(nèi)等多個角度探討金沙江上中游梯級水庫運(yùn)行對流域徑流變化特征的影響,以期為金沙江流域的水資源高效利用和水資源管理提供參考和依據(jù)。

    1 研究區(qū)域概況與資料來源

    金沙江流域地處中國青藏高原、橫斷山區(qū)、云貴高原及四川盆地西部邊緣,位于東經(jīng)90°23′~104°37′,北緯24°28′~35°46′之間,流域集水面積約50萬km2,占長江上游流域控制總面積的近一半。流域地勢西高東低,由西北逐漸向東南傾斜[13-14。金沙江干流一般分為上、中、下3 段,上游為直門達(dá)至石鼓段,中游為石鼓至攀枝花段,下游為攀枝花至宜賓段(圖1)。流域南北縱越,橫跨十余個經(jīng)度及緯度,地形地貌極為復(fù)雜,氣候特征差異極大[15。流域多年平均降水量約在250~1 200 mm之間,多年平均氣溫在-5~20 ℃之間,降水和冰川融雪是流域徑流的主要來源,水文過程對氣候變化較為敏感[16。

    根據(jù)金沙江流域水系分布特征,選擇統(tǒng)計資料系列較長的石鼓水文站以及攀枝花水文站作為金沙江干流上、中游控制站,水文站的徑流資料來自于長江水利委員會水文局,資料系列長度為1970~2022年,時間尺度為月。

    2 研究方法

    Mann-Kendall(M-K)趨勢檢驗法是世界氣象組織(WMO)推薦并已被廣泛使用的一種非參數(shù)檢驗方法,用于衡量數(shù)據(jù)序列的變化趨勢。其優(yōu)點(diǎn)是不要求樣本遵從一定的分布,也不受少數(shù)異常值的干擾,計算簡便,在水文序列分析中得到了廣泛使用[17?;瑒觮檢驗法是一種通過構(gòu)造統(tǒng)計量分析兩組樣本平均值的差異是否顯著來檢驗突變的方法,其基本思想是如果兩段子序列的特征值差異超過了一定的顯著性水平,則可以認(rèn)為有突變發(fā)生[18。小波功率譜分析是頻譜周期性分析方法中的一種,其在小波分析的基礎(chǔ)上增加了紅噪聲標(biāo)準(zhǔn)譜檢驗,在研究序列的周期性問題上應(yīng)用較為廣泛19。本研究擬采用Mann-Kendall 趨勢檢驗法分析金沙江上中游流域兩個控制站年徑流系列的趨勢性,滑動t檢驗法被用來檢驗控制站年徑流系列的突變性,小波功率譜法被用于分析徑流序列的周期性變化特征。

    水庫的運(yùn)行調(diào)度對流域河川徑流演變特征有一定影響,通過分析徑流的特征值變化可進(jìn)一步深層次地探究水庫運(yùn)行對徑流變化的影響。本研究采用不均勻系數(shù)、完全調(diào)節(jié)系數(shù)、集中度[20和基尼系數(shù)分析金沙江上游梯級水庫建立后對上中游重要控制節(jié)點(diǎn)徑流的影響。

    (1)不均勻系數(shù)C計算公式為

    C=1n∑ni=1(R-R—)2R—(1)

    式中:C為不均勻系數(shù);n表示年內(nèi)徑流或水位的天數(shù);R表示年內(nèi)日徑流或日水位值;R—表示年平均值。C值越大,表明徑流、水位在年內(nèi)不同時期相差越大,分配越不均勻。

    (2)完全調(diào)節(jié)系數(shù)C,與C類似,是另一種反映年內(nèi)分配情況的指標(biāo),其計算公式為

    C=∑ni=1ψ(i)(R-R—)∑ni=1R(2)

    其中:

    ψ(i)=0 R<R—1 R≥R—(3)

    (3)集中度C。集中度是分析徑流、水位等水文變量在年內(nèi)集中程度的指標(biāo),其值越大,表明年內(nèi)分配越為集中?;谌諒搅骱腿账粩?shù)據(jù),將日徑流、水位視為向量,分解成x和y兩個方向上的分量:

    R=∑ni=1Rsinθ

    R=∑ni=1Rcosθ

    θ=2πin

    R=R2+R2(4)

    式中:θ表示水文日變量序列的方向。相應(yīng)地,集中度C計算公式為

    C=R/∑ni=1R(5)

    (4)基尼系數(shù)是用于衡量一個國家或地區(qū)財富、收入分配不平等的指標(biāo),近年來被不少學(xué)者應(yīng)用于水文學(xué)領(lǐng)域,分析水文變量在年內(nèi)分配的均勻度?;嵯禂?shù)越大,說明水文變量年內(nèi)分布均勻度越低,年內(nèi)分配越不均勻[21?;嵯禂?shù)可通過繪制洛倫茲曲線進(jìn)行計算。洛倫茲曲線是外凸曲線,如圖2所示。絕對平等線是與橫坐標(biāo)軸呈45°的直線,表明徑流、水位等水文變量分配絕對均勻。基尼系數(shù)為絕對平等線與洛倫茲曲線所圍成面積的2倍,即計算步驟為:① 將控制站點(diǎn)的各年日徑流、日水位數(shù)據(jù)等水文變量按升序排列,計算累積百分比;② 計算水文站的時間累積百分比;③ 以時間累積百分比為自變量x,水文變量累積百分比為因變量y,擬合洛倫茲曲線;④ 根據(jù)洛倫茲曲線計算水文變量的基尼系數(shù)。

    3 結(jié)果與討論

    3.1 年徑流演變規(guī)律

    金沙江上中游流域1970~2022年主要控制站年均流量的演變規(guī)律如圖3所示。在置信水平為 0.05 的情況下,石鼓站的年平均徑流量序列的M-K趨勢檢驗統(tǒng)計量Z=0.55<1.96,攀枝花站統(tǒng)計量Z=0.53<1.96,表明金沙江上中游流域1970~2022年間的年均流量總體呈現(xiàn)不顯著上升趨勢。由圖3可知,石鼓站1970~2022年間多年平均流量為1 334 m3/s,最大年均流量出現(xiàn)在2005年,年均流量值達(dá)1 714 m3/s,最小年均流量出現(xiàn)在1994年,為930 m3/s。石鼓站年均流量線性趨勢線結(jié)果表明,該站年均流量總體呈現(xiàn)增加趨勢,線性傾向率為41.4(m3·s-1)/10 a,5 a滑動平均線表明,該站年均流量自20世紀(jì)70年代呈趨勢性增加并伴隨幅度波動變化,但在21世紀(jì)初呈現(xiàn)下降趨勢。攀枝花站1970~2022年間多年平均流量為1 801 m3/s,最大年均流量出現(xiàn)在1998年,達(dá)2 406 m3/s,最小年均流量出現(xiàn)在1994年,為1 209 m3/s。攀枝花站年均徑流量總體呈現(xiàn)線性增加趨勢,線性傾向率為45.6(m3·s-1)/10 a,其5 a滑動平均線變化表明,該站年均流量自20世紀(jì)70年代呈增加趨勢,在21世紀(jì)初出現(xiàn)轉(zhuǎn)折點(diǎn),而后呈現(xiàn)下降趨勢。

    總體上看,金沙江上游和中游流域的年徑流量在1970~2022年間呈現(xiàn)出不顯著上升趨勢。其中,上游流域自20世紀(jì)80年代到21世紀(jì)初為流域的豐水期;而中游流域自20世紀(jì)90年代到21世紀(jì)初為流域的豐水期,年代際徑流均值較多年來水量偏高;2010年之后,金沙江中游以上流域的年代際徑流均呈現(xiàn)降低的趨勢,但是仍高于流域多年平均來水量,這與流域的氣候變化和人類活動有關(guān)。相關(guān)研究表明,金沙江流域2010年以來降水量偏少[22-24;另一方面,攀枝花水文站上游的金安橋水電站于2011年12月開始蓄水,梨園、阿海、龍開口、魯?shù)乩⒂^音巖水庫于2012年后陸續(xù)蓄水,水電站的調(diào)度運(yùn)行進(jìn)一步影響了金沙江流域的年徑流量變化。空間上,從金沙江上游到中游,徑流的增加幅度越來越大,位于金沙江中游的攀枝花站多年平均流量比位于金沙江上游的石鼓站多年平均流量增加了24.5%,在汛期增加的趨勢更為顯著。除此之外,金沙江流域上中游的年徑流變化趨勢具有一定的同步性,最小年均流量的出現(xiàn)時間具有一致性,但最大年均流量的出現(xiàn)時間并不匹配,且攀枝花站的年均流量都高于石鼓站,其變化的波動性也相對高于石鼓站,這也與流域的空間分布特征相符合。

    金沙江上中游主要控制站1970~2022年間的年均流量滑動t檢驗與小波功率譜圖如圖4~5所示。根據(jù)年均流量的突變性檢驗結(jié)果,石鼓站年徑流系列的突變點(diǎn)分別出現(xiàn)在1973、1979、1984、1986、1988、1997年,其中1986年為最顯著突變點(diǎn),滑動t檢驗值為2.83。而攀枝花站年徑流系列的突變點(diǎn)分別出現(xiàn)在1984、1986、1988、1997年,其中1997年突變最為顯著,滑動t檢驗值為3.01。選擇兩站最為顯著的兩個突變點(diǎn)(1986年與1997年)分別分析突變前后兩階段的徑流變化特征。結(jié)果表明:石鼓站、攀枝花站的年徑流在突變點(diǎn)之后分別相較于突變前增加了10%和12%。因此,金沙江流域的年徑流從上游到中游增加的幅度越來越大,徑流的突變點(diǎn)越來越滯后。

    圖5是金沙江上中游主要控制站實(shí)測年徑流系列小波功率譜圖,由圖可知,石鼓站和攀枝花站年徑流系列的小波功率譜中標(biāo)準(zhǔn)噪聲譜檢驗區(qū)域十分接近,表明金沙江上中游流域的徑流周期性變化規(guī)律較為接近。從時間尺度通過0.05顯著性水平的紅噪聲標(biāo)準(zhǔn)譜檢驗與右圖對應(yīng)的紅色虛線內(nèi)的小波方差峰值來看,金沙江上中游流域的年徑流具有2~5 a的周期性,且從波譜的有效閾值來看,從上游到中游的波譜顯著性區(qū)域變少,周期顯著性降低。

    3.2 月徑流演變規(guī)律

    圖6是石鼓站與攀枝花站1970~2022年間徑流年內(nèi)分配情況。結(jié)果表明:石鼓站年徑流分配不均,徑流主要集中在6~10月,與流域汛期時段相吻合;1970~2022年間,石鼓站汛期徑流量占全年徑流量的74%左右;最大月均流量出現(xiàn)在8月,為2 899 m3/s,最小月均流量出現(xiàn)在2月,為400 m3/s。同樣地,攀枝花站年徑流分配也不均勻,徑流主要集中在6~10月;1970~2022年間,攀枝花站汛期流量占全年流量的74.5%;最大月均流量出現(xiàn)在8月,為4 110 m3/s,最小月均流量出現(xiàn)在3月,為539 m3/s。

    由圖6可知,石鼓站和攀枝花站的徑流年內(nèi)分配特征基本一致,月最大流量出現(xiàn)時間均為8月,但月最小流量出現(xiàn)時間不一致,攀枝花站最小流量出現(xiàn)時間滯后于石鼓站,且攀枝花站的月均流量均高于石鼓站,符合流域的空間分布特征與徑流傳播的滯后性。

    圖7是金沙江流域主要控制站月徑流M-K檢驗結(jié)果。由圖7可知,金沙江中游的月均徑流均呈現(xiàn)不顯著增加的趨勢,其中,石鼓站10月的月徑流增加趨勢較為顯著,而攀枝花站則是11月的月徑流增加趨勢較為顯著。整體來看,金沙江中游豐水期的月徑流增加趨勢弱于枯水期,但月徑流增加趨勢強(qiáng)于上游。在8月和10月,石鼓站月徑流的增加趨勢高于攀枝花站,且7~10月間,兩站的月徑流增加趨勢較為接近,差距相對低于其他月。這就表明,在金沙江上中游流域的豐水期,上游和中游的徑流增加趨勢具有高度一致性,其他時期上游徑流的增加趨勢弱于中游。

    3.3 梯級水庫建設(shè)運(yùn)行對徑流變化的影響

    考慮到2018年之前石鼓站以上無大型水庫運(yùn)行,攀枝花站上游的金安橋水電站于2011年12月開始蓄水,梨園、阿海、龍開口、魯?shù)乩⒂^音巖水庫于2012年后陸續(xù)蓄水,因此,僅以攀枝花站作為代表進(jìn)行分析,并將研究時間序列劃分為建庫前(1971~2011年)和建庫后(2012~2022年)兩個時間段進(jìn)行對比分析。圖8是梯級水庫建成前后金沙江中游攀枝花站的年徑流變化。

    圖8的結(jié)果表明,攀枝花站建庫前多年平均流量為1 792 m3/s,而建庫后多年平均流量為1 829 m3/s,梯級水庫建成之后,金沙江中游河段年均來水量增加,且變化傾向率由8.3(m3·s-1)/10 a變成35.2(m3·s-1)/10 a,來水量呈現(xiàn)增加的趨勢。結(jié)合天然情況(2012~2022年)年徑流變化趨勢來看,其基本與建庫后年徑流變化過程相吻合,且也表現(xiàn)為增加的趨勢。一方面,由于金沙江上游河段梯級水庫建成運(yùn)行時期較短,其對金沙江中游河段的年徑流影響并不明顯。另一方面,相關(guān)研究表明,金沙江中游近10 a來降水較為豐沛,氣候變化導(dǎo)致徑流增加也是不可忽視的因素之一[25。

    圖9展示了攀枝花站建庫后實(shí)測徑流相對天然徑流的年內(nèi)分配變化(2012~2022年)。結(jié)果表明,攀枝花站建庫后月平均流量相較于天然情況減少了1.5%,改變較為微弱。7~10月,攀枝花站建庫后的月均流量相較于天然情況減少了2%,而11月到次年4月,建庫后月均流量相較天然情況增加了3%。此外,天然情況下攀枝花站最大月均流量出現(xiàn)在9月,為4 093 m3/s,建庫后也出現(xiàn)在9月,為3 953 m3/s;天然情況下最小月均流量出現(xiàn)在3月,為553 m3/s,建庫后最小月均流量出現(xiàn)在2月,為570 m3/s。總體而言,梯級水庫的投運(yùn)對攀枝花站的徑流年內(nèi)分配格局有一定影響,具體表現(xiàn)為豐水期水量下降,而枯水期水量增加。另一方面,建庫后攀枝花站的最小月均流量出現(xiàn)時間提前,且月均流量峰值減少,谷值增加。

    攀枝花站建庫前后的參數(shù)對比如表1所列。由表1可知,金沙江中游梯級水庫建成之后,流域中游河段年徑流不均勻系數(shù)、完全調(diào)節(jié)系數(shù)、集中度、基尼系數(shù)都小于建庫前,表明金沙江中游河段的徑流年內(nèi)分配變得更加均勻。原因是梯級水庫陸續(xù)建成投運(yùn)后,水庫的蓄洪補(bǔ)枯調(diào)節(jié)作用使得中游河段來水量更加均勻,波動性減小。

    4 結(jié) 論

    本研究基于1970~2022年間金沙江上中游流域兩個控制站的實(shí)測徑流,采用多種方法從年際、年內(nèi)變化的角度系統(tǒng)分析了金沙江上中游流域的徑流演變規(guī)律,并根據(jù)金沙江上中游梯級水庫運(yùn)行的時間分析了建庫前后金沙江中游流域的徑流變化特征,主要結(jié)論如下:

    (1)金沙江上中游流域的年徑流在1970~2022年間呈現(xiàn)出不顯著上升趨勢,從上游到中游,徑流的增加幅度越來越大,金沙江中游的多年平均流量比上游增加了24.5%,這種增加的趨勢在汛期更為顯著。金沙江流域上中游的年徑流變化趨勢具有一定的同步性,年均最小流量的出現(xiàn)時間具有一致性,但年均最大流量的出現(xiàn)時間具有滯后性。

    (2)金沙江上游流域自20世紀(jì)80年代到21世紀(jì)初為流域的豐水期,而中游流域自20世紀(jì)90年代到21世紀(jì)初為流域的豐水期,年代際徑流均值較多年來水量偏高,2010年之后,金沙江中游以上流域的年代徑流均呈現(xiàn)降低的趨勢,但是仍高于流域多年平均來水量,這與流域的氣候變化和人類活動有關(guān)。

    (3)金沙江流域的年徑流從上游到中游增加的幅度越來越大,徑流的突變點(diǎn)越來越滯后,且從上游到中游的徑流變化周期性越來越明顯。金沙江上中游流域的徑流年內(nèi)分配特征基本一致,最大月均流量出現(xiàn)時間均為8月,但最小月均流量出現(xiàn)時間具有滯后性。在金沙江上中游流域的豐水期,上游和中游的徑流增加趨勢具有高度一致性,其他時期上游徑流的增加趨勢弱于中游。

    (4)梯級水庫的投運(yùn)對金沙江中游的徑流年內(nèi)分配格局有一定影響,豐水期水量下降,而枯水期水量增加。建庫后金沙江中游月最小流量出現(xiàn)時間提前,且月均流量峰值減少,谷值增加。梯級水庫建成之后,水庫蓄洪補(bǔ)枯的調(diào)節(jié)作用使得中游河段來水量更加均勻,波動性減小。

    金沙江流域是中國重要水電能源基地,其源頭下墊面條件復(fù)雜,存在大量冰川和積雪,在氣候變化和人類活動的雙重影響下,河川徑流的演變規(guī)律愈發(fā)復(fù)雜,如何厘定氣候變化和人類活動對金沙江流域徑流變化的影響程度是需要進(jìn)一步研究的問題。

    參考文獻(xiàn):

    [1] 徐宗學(xué),姜瑤.變化環(huán)境下的徑流演變與影響研究:回顧與展望[J].水利水運(yùn)工程學(xué)報,2022(1):9-18.

    [2] 張金萍,肖宏林,張鑫.水庫運(yùn)行對徑流-泥沙關(guān)系的影響分析[J].水電能源科學(xué),2019,37(9):17-20,50.

    [3] 趙偉,毛繼新,關(guān)見朝,等.三峽水庫運(yùn)用后荊南三口分流演變規(guī)律及驅(qū)動因子[J].水利學(xué)報,2023,54(8):1005-1014.

    [4] 《大壩與安全》編輯部.四川“十四五”重點(diǎn)建設(shè)“三江”水電基地[J].大壩與安全,2021(2):5.

    [5] 陳華,朱迪,陳啟會,等.考慮水庫調(diào)蓄的金沙江徑流模擬研究[J].武漢大學(xué)學(xué)報(工學(xué)版),2022,55(1):13-21.

    [6] 王國慶,張建云,管曉祥,等.中國主要江河徑流變化成因定量分析[J].水科學(xué)進(jìn)展,2020,31(3):313-323.

    [7] 舒章康,張建云,金君良,等.1961~2018年中國主要江河枯季徑流演變特征與成因[J].氣候變化研究進(jìn)展,2021,17(3):340-351.

    [8] 鮑振鑫,嚴(yán)小林,王國慶,等.1956~2016年黃河流域河川徑流演變規(guī)律[J].水資源與水工程學(xué)報,2019,30(5):52-57.

    [9] 康亞靜,黃婷婷,劉宇.1956~2018年長江上游寸灘站徑流演變特征[J].水利水電快報,2023,44(11):17-24.

    [10]馮勝航,王黨偉,秦蕾蕾,等.金沙江流域徑流變化特征及成因[J].南水北調(diào)與水利科技,2023,21(2):248-257.

    [11]張小峰,閆昊晨,岳遙,等.近50年金沙江各區(qū)段年徑流量變化及分析[J].長江流域資源與環(huán)境,2018,27(10):2283-2292.

    [12]郭衛(wèi),徐長江,邵駿,等.基于RVA法的金沙江下游水文情勢研究[J].人民長江,2018,49(18):58-63.

    [13]趙冬林,朱仕榮.2010~2021年金沙江干熱河谷植被覆蓋度時空變化及其影響因素[J/OL].生態(tài)學(xué)雜志,1-11[2023-12-25]http:∥kns.cnki.net/kcms/detail/21.1148.Q.20231122.0947.002.html.

    [14]陳啟會,張俊,侯雨坤,等.近60年氣候與土地利用變化對金沙江徑流的影響[J].人民長江,2018,49(23):47-53.

    [15]祁雪姣,楊穎,韓傳楠,等.氣候變化對長江源區(qū)地表徑流的影響[J].南京師大學(xué)報(自然科學(xué)版),2022,45(4):81-90.

    [16]邵駿,熊瑩,卜慧,等.長江源區(qū)基流變化規(guī)律及其氣象影響因素分析[J].人民長江,2022,53(11):61-65,71.

    [17]謝雅潔,劉曙光,周正正.汶川縣降雨時空分布特征及設(shè)計暴雨頻率分析[J].人民長江,2024,55(1):105-112.

    [18]劉琳,徐宗學(xué).西南5省市極端氣候指數(shù)時空分布規(guī)律研究[J].長江流域資源與環(huán)境,2014,23(2):294-301.

    [19]汪琳,龔麗飛,李楊,等.尼亞加拉河1951~2016年水文氣象要素演變特征及其響應(yīng)性[J].華北水利水電大學(xué)學(xué)報(自然科學(xué)版),2021,42(3):33-40.

    [20]譚寓寧,劉懷湘,陸永軍.金沙江下游區(qū)間來沙驅(qū)動因子分析及產(chǎn)沙預(yù)測模型[J].水科學(xué)進(jìn)展,2023,34(2):277-289.

    [21]劉松,陳立華,丁星臣,等.西江流域主要水文站近40年徑流變化分析研究[J].人民長江,2021,52(增2):52-55.

    [22]邢雯慧,王飛龍,邱輝,等.長江上游關(guān)鍵區(qū)主汛期典型澇年氣候成因辨析[J].人民長江,2024,55(1):91-96.

    [23]ZHANG P P,CAI Y,XIE Y L,et al.Effects of a cascade reservoir system on runoff and sediment yields in a River Basin of southwestern China[J].Ecological Engineering:The Journal of Ecotechnology,2022(179):106616.

    [24]ZHANG D,WANG W S,YU S Y,et al.Assessment of the contributions of climate change and human activities to runoff variation:case study in four subregions of the Jinsha River Basin,China[J].Journal of Hydrologic Engineering,2021,26(9):2119.

    [25]劉曉婉.金沙江流域氣象因子特征及徑流規(guī)律初步研究[D].武漢:長江科學(xué)院,2016.

    (編輯:謝玲嫻)

    Study on evolution law of runoff in upper and middle reaches of

    Jinsha River from 1970~2022WANG Lin,ZHOU Bo,GUO Wei,ZHU Di,BU Hui,ZUO Jian

    (Bureau of Hydrology,Changjiang Water Resources Commission,Wuhan 430010,China)

    Abstract: The upper and middle reaches of Jinsha River,a pivotal hydropower energy hub in Chinese strategic planning and construction,holds immense importance in exploring runoff evolution under environmental variations.This understanding is crucial for enhancing the efficient utilization and management of water resources within the basin.In this study,we delved into the runoff evolution in the upper and middle reaches of Jinsha River Basin using non-parametric test,sliding t-test,and wavelet power method,and provided a comprehensive analysis on the influence of cascade reservoir operation on runoff by examining metrics such as the inhomogeneity coefficient,complete adjustment coefficient,concentration,and Gini coefficient.The findings revealed that the annual runoff in the study area showed a non-significant upward trend.Notably,the runoff increased more and more from upstream to midstream.The average annual flow of the middle reaches increased by 24.5% compared to upstream levels,albeit with a delayed abrupt change.Additionally,the runoff periodicity from upstream to midstream became increasingly subdued.Cascade reservoir operations and flood storage regulations not only stabilized water flow in the river's middle reaches,reduced fluctuations,but also modulated runoff,decreasing in wet seasons and increasing in dry seasons,and this also led to an advance in minimum monthly flows.The research results can provide theoretical support and reference for the optimal allocation of water resources in the Jinsha River Basin.

    Key words: river runoff;evolution law;reservoir operation;Jinsha River

    亚洲av五月六月丁香网| 久久婷婷人人爽人人干人人爱| 色综合站精品国产| 中文字幕精品亚洲无线码一区| 色尼玛亚洲综合影院| 99国产精品一区二区三区| 国产精品影院久久| 国产亚洲欧美在线一区二区| 久久性视频一级片| 亚洲精品乱码久久久v下载方式| 欧美+日韩+精品| 亚洲成人久久性| 亚洲精品色激情综合| 亚洲七黄色美女视频| 亚洲欧美激情综合另类| 欧美潮喷喷水| 男人舔奶头视频| 久久精品影院6| 亚洲欧美清纯卡通| 国产又黄又爽又无遮挡在线| 动漫黄色视频在线观看| 天堂√8在线中文| 久久久久九九精品影院| 成人美女网站在线观看视频| 麻豆成人av在线观看| 国产精品av视频在线免费观看| 观看美女的网站| 少妇人妻精品综合一区二区 | 特级一级黄色大片| 五月伊人婷婷丁香| 老鸭窝网址在线观看| 99久国产av精品| 亚洲,欧美,日韩| 婷婷精品国产亚洲av在线| 国产一区二区在线av高清观看| 男人的好看免费观看在线视频| 精品久久久久久久久久久久久| 国产精品一区二区三区四区免费观看 | 国产精品国产高清国产av| 又紧又爽又黄一区二区| 色综合亚洲欧美另类图片| 长腿黑丝高跟| 制服丝袜大香蕉在线| 51午夜福利影视在线观看| 国产欧美日韩一区二区精品| 日韩欧美国产在线观看| 又紧又爽又黄一区二区| 国产精品一及| 嫁个100分男人电影在线观看| 别揉我奶头 嗯啊视频| 99久久九九国产精品国产免费| 亚洲精品成人久久久久久| 我要搜黄色片| 校园春色视频在线观看| 成人高潮视频无遮挡免费网站| 午夜福利18| 亚洲狠狠婷婷综合久久图片| 国产亚洲欧美98| 亚洲午夜理论影院| 97人妻精品一区二区三区麻豆| 极品教师在线免费播放| 欧美3d第一页| 黄色日韩在线| 我的女老师完整版在线观看| 一级作爱视频免费观看| 午夜免费男女啪啪视频观看 | 亚洲av免费高清在线观看| 精品日产1卡2卡| 成人欧美大片| 亚洲精品在线美女| 男人舔女人下体高潮全视频| 色哟哟·www| 青草久久国产| 日本在线视频免费播放| 999久久久精品免费观看国产| 亚洲,欧美精品.| 一本久久中文字幕| 精品一区二区三区av网在线观看| ponron亚洲| 精品久久久久久久久av| 亚洲成a人片在线一区二区| 91在线观看av| 久久国产乱子免费精品| 一个人观看的视频www高清免费观看| 亚洲久久久久久中文字幕| 亚洲av二区三区四区| 男女做爰动态图高潮gif福利片| 一个人观看的视频www高清免费观看| 日韩欧美在线二视频| 波野结衣二区三区在线| 欧美黑人欧美精品刺激| 精品无人区乱码1区二区| 日日摸夜夜添夜夜添av毛片 | 色哟哟·www| 少妇人妻精品综合一区二区 | 久久中文看片网| 国产亚洲av嫩草精品影院| www.熟女人妻精品国产| 久久久久久久久久成人| 欧美三级亚洲精品| 亚洲乱码一区二区免费版| 日韩欧美三级三区| 欧美日韩乱码在线| 日韩欧美精品免费久久 | 国产高潮美女av| 亚洲欧美日韩无卡精品| 久久这里只有精品中国| 精品久久久久久久久av| 啪啪无遮挡十八禁网站| 一进一出抽搐动态| 人妻久久中文字幕网| 欧美最黄视频在线播放免费| 欧美精品啪啪一区二区三区| 99久久久亚洲精品蜜臀av| 精品久久久久久久久亚洲 | 国产精品自产拍在线观看55亚洲| 十八禁人妻一区二区| 成人鲁丝片一二三区免费| 国产精品一区二区免费欧美| 在线观看舔阴道视频| 午夜福利在线观看免费完整高清在 | 在线天堂最新版资源| 午夜福利在线观看免费完整高清在 | 午夜日韩欧美国产| 国产欧美日韩精品一区二区| 国产午夜精品论理片| 熟妇人妻久久中文字幕3abv| 国产乱人伦免费视频| av天堂在线播放| 18禁黄网站禁片免费观看直播| 亚洲,欧美精品.| 综合色av麻豆| 日本 欧美在线| 日本一二三区视频观看| 国产中年淑女户外野战色| 国产成人啪精品午夜网站| 精品无人区乱码1区二区| 国产精品永久免费网站| 亚洲成av人片在线播放无| 成年女人看的毛片在线观看| 深夜a级毛片| 亚洲av电影在线进入| 国产亚洲精品久久久com| 51国产日韩欧美| 久久久久精品国产欧美久久久| 日韩国内少妇激情av| 中文字幕av成人在线电影| 精品久久久久久,| 美女黄网站色视频| 久久久久久久精品吃奶| 中文字幕精品亚洲无线码一区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲成人免费电影在线观看| 美女大奶头视频| 亚洲国产欧洲综合997久久,| 一边摸一边抽搐一进一小说| 麻豆av噜噜一区二区三区| 中文字幕久久专区| 欧美日韩黄片免| 国产免费av片在线观看野外av| 中国美女看黄片| 脱女人内裤的视频| avwww免费| 青草久久国产| 人妻制服诱惑在线中文字幕| 精品一区二区免费观看| 欧美黑人巨大hd| 成人av在线播放网站| 亚洲av不卡在线观看| 中文字幕av在线有码专区| 人妻夜夜爽99麻豆av| 一级作爱视频免费观看| 精品免费久久久久久久清纯| 可以在线观看的亚洲视频| 丁香六月欧美| 国产主播在线观看一区二区| av女优亚洲男人天堂| 国产在线精品亚洲第一网站| 国产午夜精品久久久久久一区二区三区 | 国产亚洲精品综合一区在线观看| 免费电影在线观看免费观看| 老司机午夜福利在线观看视频| 熟女电影av网| 日本a在线网址| 中文在线观看免费www的网站| 精品人妻偷拍中文字幕| 老鸭窝网址在线观看| 成人一区二区视频在线观看| 国产高清激情床上av| 麻豆成人午夜福利视频| 国产毛片a区久久久久| 在线免费观看的www视频| 欧美xxxx性猛交bbbb| 亚洲国产精品合色在线| 久久久久久国产a免费观看| 国产精品久久久久久久久免 | 欧美zozozo另类| 中文字幕人妻熟人妻熟丝袜美| 精品人妻偷拍中文字幕| 一本久久中文字幕| 天天一区二区日本电影三级| 亚洲熟妇熟女久久| 精品日产1卡2卡| 身体一侧抽搐| 制服丝袜大香蕉在线| 观看免费一级毛片| 国产不卡一卡二| 男人和女人高潮做爰伦理| 色综合婷婷激情| 搡女人真爽免费视频火全软件 | 精品乱码久久久久久99久播| 亚洲第一欧美日韩一区二区三区| 国产乱人视频| 十八禁人妻一区二区| 色综合婷婷激情| 精品一区二区免费观看| 国产v大片淫在线免费观看| 国产午夜福利久久久久久| 亚洲av一区综合| 老女人水多毛片| 亚洲乱码一区二区免费版| 男人的好看免费观看在线视频| 91九色精品人成在线观看| 在线免费观看的www视频| 夜夜夜夜夜久久久久| 女同久久另类99精品国产91| 欧美一区二区亚洲| 国产大屁股一区二区在线视频| 亚洲精品一区av在线观看| 亚洲欧美清纯卡通| 白带黄色成豆腐渣| 欧美性猛交黑人性爽| 国产欧美日韩一区二区精品| 欧美不卡视频在线免费观看| 国产精品免费一区二区三区在线| 欧美色视频一区免费| 91九色精品人成在线观看| 国产私拍福利视频在线观看| 校园春色视频在线观看| 亚洲自拍偷在线| 精品欧美国产一区二区三| 欧美激情久久久久久爽电影| 嫩草影院精品99| 国产av一区在线观看免费| av在线观看视频网站免费| 午夜福利视频1000在线观看| 免费无遮挡裸体视频| 丰满乱子伦码专区| 淫妇啪啪啪对白视频| 91av网一区二区| 欧美日韩福利视频一区二区| 成人特级黄色片久久久久久久| 免费看美女性在线毛片视频| 国产激情偷乱视频一区二区| 久久久久久久久久黄片| 校园春色视频在线观看| 亚洲国产精品999在线| 麻豆成人午夜福利视频| 成年女人看的毛片在线观看| 18美女黄网站色大片免费观看| 可以在线观看毛片的网站| 久久亚洲精品不卡| 首页视频小说图片口味搜索| 91午夜精品亚洲一区二区三区 | 69av精品久久久久久| 成人精品一区二区免费| 一二三四社区在线视频社区8| 亚洲男人的天堂狠狠| 国产精品久久久久久亚洲av鲁大| 国产淫片久久久久久久久 | 十八禁网站免费在线| 特大巨黑吊av在线直播| 99国产极品粉嫩在线观看| 午夜福利在线观看吧| 人人妻人人看人人澡| 亚洲综合色惰| 99热这里只有精品一区| 深夜a级毛片| 成人性生交大片免费视频hd| 成年女人看的毛片在线观看| 久久午夜福利片| 日日摸夜夜添夜夜添小说| 他把我摸到了高潮在线观看| 欧美3d第一页| 最好的美女福利视频网| 美女 人体艺术 gogo| 全区人妻精品视频| 日韩av在线大香蕉| 天天一区二区日本电影三级| 极品教师在线免费播放| 午夜免费男女啪啪视频观看 | 国产高潮美女av| 日韩欧美精品免费久久 | 久久性视频一级片| 国产精品1区2区在线观看.| 国内毛片毛片毛片毛片毛片| 韩国av一区二区三区四区| 变态另类丝袜制服| 亚洲av成人精品一区久久| 亚洲国产精品sss在线观看| 亚洲国产色片| 三级国产精品欧美在线观看| 男人的好看免费观看在线视频| 婷婷亚洲欧美| 69av精品久久久久久| 在线观看舔阴道视频| 精品久久久久久久久av| 国产精品电影一区二区三区| 亚洲av美国av| 亚洲第一电影网av| 热99在线观看视频| 亚洲男人的天堂狠狠| 亚洲欧美日韩无卡精品| 性欧美人与动物交配| 中文字幕av成人在线电影| 一本一本综合久久| 999久久久精品免费观看国产| 色精品久久人妻99蜜桃| 免费无遮挡裸体视频| 51午夜福利影视在线观看| 欧美日本视频| 国产精品自产拍在线观看55亚洲| 免费av观看视频| 国产欧美日韩精品亚洲av| 乱码一卡2卡4卡精品| 热99在线观看视频| 别揉我奶头 嗯啊视频| www日本黄色视频网| 国产亚洲欧美98| 久久精品人妻少妇| 熟女人妻精品中文字幕| 国产高清视频在线播放一区| 真人一进一出gif抽搐免费| 亚洲成人免费电影在线观看| 日本一二三区视频观看| 亚洲第一区二区三区不卡| www.www免费av| 少妇高潮的动态图| 一区二区三区四区激情视频 | 深夜精品福利| 99热这里只有是精品50| 五月玫瑰六月丁香| 国产探花在线观看一区二区| 床上黄色一级片| 黄色丝袜av网址大全| 国产不卡一卡二| 99久国产av精品| 午夜影院日韩av| 欧美色欧美亚洲另类二区| 嫩草影视91久久| 桃红色精品国产亚洲av| 久久精品91蜜桃| 国产精品久久视频播放| 欧美区成人在线视频| 男女视频在线观看网站免费| 欧美不卡视频在线免费观看| 最近最新免费中文字幕在线| av福利片在线观看| 久久久久九九精品影院| 成年女人看的毛片在线观看| 久9热在线精品视频| 亚洲国产高清在线一区二区三| 日韩欧美精品免费久久 | 日韩有码中文字幕| 嫩草影视91久久| 热99re8久久精品国产| 久久热精品热| 麻豆成人午夜福利视频| 大型黄色视频在线免费观看| 在线看三级毛片| 久久久久亚洲av毛片大全| 国产成人aa在线观看| 99久久久亚洲精品蜜臀av| 麻豆国产av国片精品| 国产免费一级a男人的天堂| 一本一本综合久久| 午夜福利高清视频| 国产日本99.免费观看| 精品人妻1区二区| 中文字幕精品亚洲无线码一区| 嫩草影视91久久| 午夜激情福利司机影院| 国产野战对白在线观看| 特大巨黑吊av在线直播| 久久亚洲精品不卡| 国产大屁股一区二区在线视频| 国产91精品成人一区二区三区| 久久久久精品国产欧美久久久| 亚洲综合色惰| 欧美一级a爱片免费观看看| 欧美精品国产亚洲| 99精品在免费线老司机午夜| 欧美一区二区国产精品久久精品| 91九色精品人成在线观看| 深爱激情五月婷婷| 欧美乱色亚洲激情| 69av精品久久久久久| 国产白丝娇喘喷水9色精品| 91久久精品国产一区二区成人| 亚洲国产精品合色在线| 欧美成人一区二区免费高清观看| 亚洲成人免费电影在线观看| 久久久国产成人免费| 中文字幕av成人在线电影| 国产中年淑女户外野战色| 免费看a级黄色片| 性欧美人与动物交配| 搡老岳熟女国产| 欧美极品一区二区三区四区| 丰满的人妻完整版| 日韩大尺度精品在线看网址| 国产精品爽爽va在线观看网站| 麻豆一二三区av精品| 国产精品爽爽va在线观看网站| 搡老熟女国产l中国老女人| 午夜亚洲福利在线播放| 婷婷精品国产亚洲av| 午夜福利在线观看吧| 成人特级黄色片久久久久久久| 日韩亚洲欧美综合| 国产三级中文精品| 蜜桃亚洲精品一区二区三区| 男女做爰动态图高潮gif福利片| 国产在视频线在精品| 美女高潮喷水抽搐中文字幕| 99精品在免费线老司机午夜| 极品教师在线免费播放| 激情在线观看视频在线高清| 十八禁人妻一区二区| 国内精品美女久久久久久| 亚洲国产色片| 欧美日本亚洲视频在线播放| 一卡2卡三卡四卡精品乱码亚洲| 最新中文字幕久久久久| 国产高清三级在线| 国产成人福利小说| 欧美乱色亚洲激情| 久久午夜福利片| 永久网站在线| 免费在线观看成人毛片| 亚洲国产精品sss在线观看| 亚洲男人的天堂狠狠| 国产精品一及| 脱女人内裤的视频| 国内毛片毛片毛片毛片毛片| 婷婷精品国产亚洲av| 国产亚洲av嫩草精品影院| 久久伊人香网站| 亚洲精品成人久久久久久| 直男gayav资源| 男女床上黄色一级片免费看| 看黄色毛片网站| 精品免费久久久久久久清纯| 欧美日韩综合久久久久久 | 51国产日韩欧美| 久久国产乱子伦精品免费另类| 亚洲av电影不卡..在线观看| 色吧在线观看| 波多野结衣巨乳人妻| 日本免费a在线| 老熟妇乱子伦视频在线观看| 免费观看的影片在线观看| 日韩人妻高清精品专区| 久久草成人影院| 九九久久精品国产亚洲av麻豆| 日本一本二区三区精品| 噜噜噜噜噜久久久久久91| 欧美日韩瑟瑟在线播放| 嫩草影院入口| 国内精品美女久久久久久| 国产中年淑女户外野战色| 欧美在线黄色| 十八禁国产超污无遮挡网站| 中文字幕熟女人妻在线| 婷婷亚洲欧美| 1024手机看黄色片| 自拍偷自拍亚洲精品老妇| 亚洲久久久久久中文字幕| 亚洲自偷自拍三级| 国产熟女xx| 欧美成人a在线观看| 色哟哟哟哟哟哟| 久久99热6这里只有精品| eeuss影院久久| 久久国产精品人妻蜜桃| 综合色av麻豆| 精品福利观看| 最近在线观看免费完整版| 久久欧美精品欧美久久欧美| 欧美高清性xxxxhd video| 天堂动漫精品| 欧美激情国产日韩精品一区| 97热精品久久久久久| 国内少妇人妻偷人精品xxx网站| 99热6这里只有精品| 丁香欧美五月| 亚洲av成人精品一区久久| 国产精品自产拍在线观看55亚洲| 亚洲精品一卡2卡三卡4卡5卡| 日韩亚洲欧美综合| 少妇人妻一区二区三区视频| 国产色爽女视频免费观看| 淫秽高清视频在线观看| 在线十欧美十亚洲十日本专区| 毛片一级片免费看久久久久 | 亚洲精品一区av在线观看| 别揉我奶头 嗯啊视频| 99热只有精品国产| netflix在线观看网站| 性插视频无遮挡在线免费观看| 悠悠久久av| 久久草成人影院| 欧美黄色片欧美黄色片| netflix在线观看网站| 最近最新免费中文字幕在线| 亚洲 国产 在线| 婷婷丁香在线五月| 亚洲国产欧洲综合997久久,| 一个人看视频在线观看www免费| 男人的好看免费观看在线视频| 午夜福利在线在线| 内地一区二区视频在线| 国产三级黄色录像| 99久久九九国产精品国产免费| 最近在线观看免费完整版| 亚洲国产精品999在线| 真实男女啪啪啪动态图| 美女xxoo啪啪120秒动态图 | 在现免费观看毛片| 少妇人妻一区二区三区视频| 亚洲在线观看片| 俄罗斯特黄特色一大片| 全区人妻精品视频| 国产69精品久久久久777片| 国产私拍福利视频在线观看| 日日摸夜夜添夜夜添av毛片 | 又爽又黄a免费视频| 露出奶头的视频| .国产精品久久| 国产成人av教育| 热99在线观看视频| 午夜免费男女啪啪视频观看 | 深爱激情五月婷婷| 成人特级黄色片久久久久久久| 亚洲成a人片在线一区二区| 国产精品一及| 国产成+人综合+亚洲专区| 午夜福利视频1000在线观看| 日韩成人在线观看一区二区三区| 99热这里只有是精品50| 69人妻影院| 成人av一区二区三区在线看| 精品久久久久久成人av| 欧美一区二区国产精品久久精品| .国产精品久久| 免费在线观看亚洲国产| 国产精品,欧美在线| 国产成+人综合+亚洲专区| 欧美乱妇无乱码| 久久国产乱子免费精品| 精品国内亚洲2022精品成人| av视频在线观看入口| 国产爱豆传媒在线观看| 亚洲不卡免费看| 亚洲av电影在线进入| 国产精品亚洲一级av第二区| 色视频www国产| 男女床上黄色一级片免费看| 人妻制服诱惑在线中文字幕| 亚洲人成电影免费在线| 黄色女人牲交| 中文资源天堂在线| 欧美乱妇无乱码| 3wmmmm亚洲av在线观看| 真人做人爱边吃奶动态| 嫩草影院入口| 亚洲真实伦在线观看| 日本一本二区三区精品| 欧美成人性av电影在线观看| 在线观看av片永久免费下载| 又爽又黄a免费视频| 超碰av人人做人人爽久久| 精品一区二区三区人妻视频| 特大巨黑吊av在线直播| 国产精品日韩av在线免费观看| 三级男女做爰猛烈吃奶摸视频| 夜夜夜夜夜久久久久| 国产精品一区二区三区四区久久| 久久久久国产精品人妻aⅴ院| 国产真实乱freesex| 黄色女人牲交| 国产伦在线观看视频一区| 欧美乱妇无乱码| 男女下面进入的视频免费午夜| 久久久久久久久中文| 午夜福利视频1000在线观看| 欧美在线黄色| 午夜免费男女啪啪视频观看 | 一区二区三区四区激情视频 | 国产成人欧美在线观看| 国产精品一区二区三区四区久久| 伦理电影大哥的女人| 在线看三级毛片| www.999成人在线观看| 国产黄片美女视频| aaaaa片日本免费| 日韩欧美三级三区| 少妇高潮的动态图| 欧美色欧美亚洲另类二区| www.999成人在线观看| 亚洲av二区三区四区| 最新在线观看一区二区三区| 亚洲熟妇熟女久久| 又紧又爽又黄一区二区| 老司机午夜福利在线观看视频|