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

    長江上游水庫群聯(lián)合調(diào)度下的河流水文情勢研究

    2019-10-20 15:00:48張康楊明祥梁藉閆寶偉
    人民長江 2019年2期

    張康 楊明祥 梁藉 閆寶偉

    摘要:長江上游水庫群的聯(lián)合調(diào)度運行必然會對下游河道徑流產(chǎn)生影響。選擇北碚、高場、宜昌站作為水庫群串、并、混聯(lián)運行下的控制斷面,通過IHA法選擇出受水庫群聯(lián)合調(diào)度影響的14個水文指標,并采用CRITIC法改進RVA完成串并混聯(lián)運行模式下的水文指標變化趨勢和整體水文變異程度的研究。分析結(jié)果表明:① 隨著水庫群聯(lián)合運行數(shù)量的增加,長江干支流河道水文情勢改變幅度增大;② 水庫群聯(lián)合運行會使汛期流量顯著降低,非汛期流量增加,同時會造成年最大(?。┝髁砍蕼p少(增加)趨勢;③ 水庫串聯(lián)運行模式會使最大流量發(fā)生時間延遲;并聯(lián)運行模式會使汛期流量坦化現(xiàn)象更加顯著;混聯(lián)運行模式對水文情勢的改變程度主要受調(diào)度方式和氣候變化的影響,上游混聯(lián)模式使宜昌站水文情勢為低改變度。

    關(guān)鍵詞:整體水文改變程度; CRITIC法; 水庫串并混聯(lián)運行; 長江上游

    中圖法分類號:P333.6文獻標志碼: ADOI:10.16232/j.cnki.1001-4179.2019.02.020

    長江上游是我國水電開發(fā)程度最高的水系[1],在1980年之后陸續(xù)建成數(shù)十座大型水電站,形成包含串并混聯(lián)的梯級水庫群。然而,隨著長江上游水庫群數(shù)量增加,使天然狀態(tài)下的極值流量、汛期流量等水文情勢發(fā)生明顯變化,出口控制水文站徑流的時空分布也呈現(xiàn)新情勢,這些變化將對中下游防洪和水資源利用帶來明顯影響。

    在水庫運行對河流水文情勢影響的研究上,Hu等[2]將IHA法用于淮河流域徑流變化方面的研究,發(fā)現(xiàn)水利工程對枯季流量影響很大。Gao等[3]采用長江上游降雨資料模擬宜昌徑流變化,并用生態(tài)指標與IHA因子進行了比較。黎云云等[4]發(fā)現(xiàn)黃河干流控制性水庫越多,水文情勢改變度越高。段唯鑫等[5]用IHA(RVA)法分析了長江上游大型水庫群對宜昌站水文情勢的影響。蔡文君[6]等用RVA法評估了三峽水庫運行對長江中下游水文情勢的影響,結(jié)果表明:在不同的豐枯水期,三峽水庫運行對中下游影響的程度是不同的,呈現(xiàn)沿河流流向逐漸減弱的規(guī)律。李林等[7]對長江上游徑流量變化進行影響因子關(guān)系分析,得出了夏季降水減少是造成徑流量減少的主要原因。然而RVA法在評估河道水體改變度上也存在一些不足。Yin等人認為RVA法未考慮水文特征的年變化,忽略了其對水文變化的影響[8]。周毅等[9]采用CRITIC法結(jié)合歐式距離對傳統(tǒng)RVA進行了改進,解決了傳統(tǒng)RVA方法存在忽略水文年變化特征及中低度改變指標的缺陷。

    目前,對長江上游梯級水電開發(fā)影響的研究較多,其中程根偉等[10]和楊麗虎等[11]側(cè)重于梯級水電開發(fā)對水環(huán)境的影響分析。但是對于水庫串并混聯(lián)運行造成河流水文要素影響的研究仍然不足。本文采用長江上游干支流水庫群和關(guān)鍵水文站的實測資料,分別研究水庫群串并混聯(lián)運行對下游水文情勢的影響程度。具體研究方案如下:采用M-K秩相關(guān)檢驗和滑動t檢驗對北碚、高場、寸灘站的徑流系列進行顯著性分析,并結(jié)合水庫建設(shè)情況劃分受影響徑流序列。采用CRITIC法來改進RVA方法,考慮各個指標相關(guān)性來確定指標權(quán)重并計算水庫群串并混聯(lián)運行下的水文情勢改變度。本文研究長江水庫群運行調(diào)度對下游河流情勢的影響,分析水文循環(huán)演變規(guī)律,對長江上游水資源統(tǒng)一管理、下游防洪調(diào)度、生態(tài)調(diào)度、河道治理、水資源開發(fā)利用等方面具有實用性意義和指導(dǎo)作用。

    1研究區(qū)概況及數(shù)據(jù)

    長江上游水資源豐沛,天然落差大,其干支流河段有大量建成或在建的梯級水庫群,主要的水庫及水文站分布如圖1所示。2000年以后,大型水庫建設(shè)主要集中于金沙江中下游、雅礱江、大渡河、烏江流域等[12]。為分析長江上游梯級水庫群串并混聯(lián)運行對下游河段水文情勢的影響,有針對性地選取關(guān)鍵水文站點十分重要。

    大渡河流域的瀑布溝水庫和岷江流域的紫坪鋪水庫均屬大(一)型水庫,且兩水庫之間是并聯(lián)運行關(guān)系。高場站位于大渡河與岷江交匯處的下游,控制面積為13.8萬km?徑流受大渡河與岷江共同影響,而不受長江沿線其他干支流的干擾,利于分析并聯(lián)水庫運行對河流水文情勢影響。北碚站是嘉陵江干流下游控制站,該站控制嘉陵江干流和上游渠江、涪江來水及掌握匯流后的水情變化規(guī)律,控制面積為15.67萬km?2。該站徑流不受長江沿線其他干支流的干擾,利于分析串聯(lián)水庫運行對河流水文情勢影響。宜昌站是長江上游來水來沙的總控制站,控制面積100萬km?占長江流域面積的55%。宜昌站可以監(jiān)測上游水庫混聯(lián)運行時的河道流量,具有反映水庫混聯(lián)對水文情勢影響的能力。

    本文選取高場站、北碚站、宜昌站分別作為水庫群串、并、混聯(lián)模式下的影響斷面,考慮到資料的可靠性和連續(xù)性,選擇這3個站點1950~2014年的逐日徑流數(shù)據(jù)作為基礎(chǔ)數(shù)據(jù)進行研究分析。

    圖1水庫及水文站分布Fig.1Distribution of reservoirs and hydrological stations

    2研究方法

    2.1突變檢驗

    目前,突變統(tǒng)計分析在應(yīng)用中仍然存在一些問題,如突變點缺乏物理解釋,或者使用突變檢驗的方法不當可能得到錯誤的結(jié)論。為避免出現(xiàn)上述問題,本文使用M-K檢驗和滑動t檢驗進行比較分析,同時給定嚴格的顯著性水平進行檢驗,綜合得到徑流序列的突變節(jié)點。

    2.1.1M-K檢驗

    Mann-Kendall(簡稱M-K法)是一種非參數(shù)檢驗方法,也是世界氣象組織(WMO)推薦廣泛使用的一種方法[13]。M-K檢驗不需要樣本遵從一定的分布,也不受少數(shù)異常值的干擾,適用于如氣象、水文等順序變量。

    2.1.2滑動t檢驗

    若兩段子序列的均值差異超過了一定的顯著水平,則認為均值發(fā)生了質(zhì)變,有突變發(fā)生。對于有n個樣本量的時間序列x,人為地設(shè)置某一時刻為基準點,基準點前后兩端子序列x1和x2的樣本為n1和n2,平均值為x1和x2,方差為S?21和S?22,定義統(tǒng)計量:

    t=x1-x2SW1n1+1n2t(n1+n2-2)(1)

    SW=n1S?21+n2S?22n1+n2-2(2)

    給定顯著水平a,若|ti|>t,則認為發(fā)生了突變。采用滑動辦法連續(xù)設(shè)置基準點,基礎(chǔ)點前后兩個子序列的長度一般相同,n1=n2。依次計算統(tǒng)計量,統(tǒng)計序列為ti,其中i=1,2,…,n-(n1+n2)+1。根據(jù)顯著水平,查t,若|ti|

    2.2水文情勢變化分析

    (1) 水文指標改變度。水文指標改變度的定義如下:

    Di=|Ni-NeNe|×100%(3)

    式中,Di為第i個水文指標改變度,Ni為第i個水文指標變異后仍落于RVA目標范圍內(nèi)的年數(shù),Ne為水文指標變異后預(yù)期落于RVA目標范圍內(nèi)的年數(shù),即水文指標變異后的年數(shù)乘以變異前水文指標落入RVA目標范圍的比例。當Di值介于0~30%時為低度改變;介于33%~67%時為中度改變;介于67%~100%時則為高度改變。

    (2) 基于CRITIC法確定水文權(quán)重。CRITIC法是由Diakoulaki.D提出的通過分析各指標間的對比強度和沖突性來綜合衡量指標客觀權(quán)重的方法[14-16]。該方法能考慮IHA指標間的相關(guān)性,進而確定各指標的權(quán)重。其中對比強度用標準差來衡量,標準差大(?。﹦t該指標反映信息量大(?。?。指標間沖突性通過相關(guān)系數(shù)表示,兩指標間負相關(guān)系數(shù)大(小)則沖突性大(?。?,說明這兩個指標在評價中反映的信息有較大(?。┎煌笜怂x權(quán)重也相應(yīng)大(?。?。

    Cj=σjmj=1(1-τij)j=1,2,…,m(4)

    Wj=Cjmj=1(Cj)j=1,2,…,m(5)

    式中,Cj為第j個指標所表示的全部信息量;σj為j個指標的標準差;τij為第i個指標和第j個指標相關(guān)系數(shù);Wj為第j個指標權(quán)重。

    在各指標權(quán)重確定后,結(jié)合RVA法評估河道水文指標整體改變度,公式為

    D=mj=1WmiDi(6)

    式中,Wmi為各指標的綜合權(quán)重;Di為各指標的改變度;D的取值范圍是從0到1,其值越大表示河流水文情勢變化越明顯。

    3結(jié)果及分析

    3.1研究時段劃分

    RVA法是基于IHA指標的用以評價水利工程對天然徑流水文情勢影響的方法。該方法需將高場、北碚、宜昌站的徑流系列劃分成天然狀態(tài)和受影響狀態(tài)。本文利用M-K法和滑動t檢驗分別對3個水文站60 a的徑流數(shù)據(jù)進行突變點檢驗。通過顯著性水平檢驗,得到3個水文站徑流序列突變點結(jié)果(見表1)。

    本文選取長江上游干支流2014以前建設(shè)完成的大型水庫(見表2),進一步分析發(fā)現(xiàn),總庫容大于50億m?3的水庫主要是1998年之后建成的。

    綜合考慮兩種突變檢驗結(jié)果和大型水庫投入運行時間,識別徑流突變的時間分界點。在高場站控制面積內(nèi),大渡河上龔嘴水庫于1972年建成運行,下游銅街子水庫于1994年建成運行,并形成兩庫串聯(lián)運行模式,對高場站水文序列造成顯著影響,岷江上紫坪鋪水庫于2006年投產(chǎn)運行,形成水庫并聯(lián)運行模式。嘉陵江上主要于1998年和2013年分別建成寶珠寺和亭子口兩座大(一)型串聯(lián)的年調(diào)節(jié)水庫,因此北碚站徑流序列以1998年以前為天然徑流,1998~2013年為寶珠寺單庫運行階段,2014~2015為兩庫串聯(lián)運行階段。宜昌站是長江上游水情的總控制站,由于長江上游于1983~1998年期間在干支流大量建設(shè)水庫,因此將此階段作為庫群建設(shè)期,1998年后長江上游水庫建設(shè)基本完成,并投入運行,其中水庫的混聯(lián)運行對水文要素的影響較為顯著,故將1999~2015年作為混聯(lián)運行階段。整體研究時段劃分見表3。

    3.2研究結(jié)果及分析

    3.2.1高場站結(jié)果分析

    高場站是大渡河與岷江的交叉把口控制站,根據(jù)突變檢驗結(jié)果,選擇1942~1972年作為天然徑流序列,以14個水文指標發(fā)生機率的25 %和75 %作為

    RVA上限和下限。計算出高場站受大渡河上龔嘴、銅街子、瀑布溝水庫和岷江上紫坪鋪水庫并聯(lián)水庫運行模式下的指標改變度,計算結(jié)果見表4。

    在1973~2009年期間,高場站月平均流量改變主要受大渡河上水庫建設(shè)運行的影響,并且隨著大渡河梯級水庫數(shù)量增加而增大。其中,龔嘴水庫運行期間,1~12月流量變化度均為低度改變。在龔嘴水庫和銅

    街子水庫聯(lián)合運行后,4,5,7,8月份流量變化度均變?yōu)橹卸雀淖?。由此可以說明大渡河梯級水庫在調(diào)蓄期內(nèi)蓄水造成對下游控制站的前汛期和汛期徑流影響變大。在2010~2015年期間,大渡河上增加瀑布溝水庫,岷江上增加紫坪鋪水庫,構(gòu)成了水庫并聯(lián)運行模式,造成高場站非汛期徑流量呈中高度改變,而汛期徑流量改變度降低為中低度改變。在水庫并聯(lián)運行影響下,河道水文要素整體改變度為38.2 %,屬于中度改變。

    大渡河、岷江流域的水庫一般在8,9月份開始蓄水,9月底基本完成蓄水任務(wù),10月份開始進入消落期,年底進入快速消落,4,5月份消落至死水位附近庫水位再上升,汛前升至汛限水位左右。分析表4中數(shù)據(jù)可知,在龔嘴水庫單庫運行時高場站4,5月份平均徑流改變度出現(xiàn)跳躍式降低,在兩庫運行和四庫并聯(lián)運行時改變度跳躍降低現(xiàn)象分別發(fā)生在6月份和7月份,說明水庫并聯(lián)運行河道徑流的影響逐漸向汛期增大。

    根據(jù)龔嘴水庫、銅街子水庫、瀑布溝水庫、紫坪鋪水庫不同組合模式運行前后高場站徑流量變化情況,得出高場站平均月流量變化趨勢,見圖2。四庫并聯(lián)運行會顯著性削減汛期峰值,使峰型變得更加矮胖,同時在非汛期因為聯(lián)合調(diào)度發(fā)電使河道流量增加。

    高場站上游水庫并聯(lián)運行造成年極端流量相對于天然狀態(tài)流量變化幅度明顯,其中年最大流量降低幅度分別達到了12.1 %,23.2 %,23.8 %,年最小流量分別增加了1.2 %,3.2 %,46.6 %。分析結(jié)果表明:隨水庫數(shù)量增加,最大(最?。┝髁拷档停ㄔ黾樱┓茸兇?。水庫并聯(lián)運行對極值流量的改變會使天然洪水和年內(nèi)豐枯周期逐步消失,危害下游生態(tài)功能,破壞水生生物的生存環(huán)境的穩(wěn)定性。

    3.2.2北碚站結(jié)果分析

    以1950~1998年為影響前的徑流序列,計算出北碚站受寶珠寺和亭子口兩個串聯(lián)水庫運行模式下的指標改變度,計算結(jié)果見表5。

    北碚站月平均流量的改變程度隨著梯級水庫數(shù)量的增加而加強,寶珠寺水庫單獨運行時,各月流量變化度為中、低度改變;寶珠寺和亭子口水庫串聯(lián)運行后12個月中有4個月都呈高度改變,其中以非汛期變化為主。北碚站月平均流量變化趨勢如圖3所示,可知北碚站月平均流量分配極不均勻,且出現(xiàn)兩次流量峰值,分別出現(xiàn)在7月和9月,在汛期7~9月份平均流量為4 753 m?3/s,非汛期平均流量為1 255 m?3/s。寶珠寺水庫單獨運行時,下游汛期來水峰值坦化程度和非汛期徑流量增加程度均較小。在寶珠寺和亭子口兩座水庫串聯(lián)運行后,兩座水庫的總庫容達到66.17 億m?3,使汛期來水峰值有明顯降低,同時汛期峰值發(fā)生時間后移。

    水庫串聯(lián)運行造成年極值流量相對于天然狀態(tài)流量變化幅度明顯,其中年最大流量降低了17.5%,年最小流量增加了41.9%,極值流量的這種變化有益于河道泥沙的沖刷,但會對下游生態(tài)功能的穩(wěn)定性和水生生物的生存環(huán)境造成危害。

    北碚站水文要素隨上游水庫數(shù)量增加改變度也增大,在水庫串聯(lián)運行影響下,河道水文要素整體改變度為40.5%,屬于中度改變。

    3.2.3宜昌站結(jié)果分析

    宜昌站是長江上游重要控制站,根據(jù)突變檢驗結(jié)果,選擇1950~1982年作為天然徑流序列,計算出宜昌站受長江上游干支流水庫混聯(lián)運行模式下的指標改變度,計算結(jié)果見表6。

    長江上游流域大量水庫多在20世紀80年代開始興建運行,因此本文將1983~1998年作為長江流域水庫集中建設(shè)階段。在此期間,宜昌站各月平均流量變化度均為低度改變。然而,在7,9月份平均流量改變度為1.4%、3.1%,在6,8月份平均流量改變度為21.9%,22.7%,出現(xiàn)了徑流改變度跳躍現(xiàn)象。其原因為:在建設(shè)運行期,6月份龔嘴、銅街子、烏江渡、大河口等水庫開始逐漸進入汛期,需要將水庫水位降至汛限水位以下,造成徑流改變度較大;7月份進入主汛期,入庫流量隨著多場洪水的發(fā)生逐漸加大,水庫實施興利調(diào)度,如水庫發(fā)電以滿發(fā)的狀態(tài)進行,7月份的改變度較小;到8月份,往往會形成較大規(guī)模洪水,水庫防洪任務(wù)較重,開始承擔削峰攔洪的任務(wù),對徑流改變度增大;按照《年度長江上游水庫聯(lián)合調(diào)度方案》,長江上游水庫群一般到9月初開始蓄水到正常蓄水位,用以滿足興利需求,水庫蓄到正常蓄水位后,基本對來水直接泄出,導(dǎo)致9月份徑流改變度較小。

    1999~2015年長江上游水庫混聯(lián)運行后,2,3月份出現(xiàn)高程度改變,1,7,10,11月均呈現(xiàn)中度改變,由此可以說明隨著長江上游梯級聯(lián)合調(diào)度水庫運行,非汛期河道徑流量顯著增加,削弱汛期洪峰影響。

    天然狀態(tài)下,高場站月平均流量分配極不均勻,呈現(xiàn)明顯的季節(jié)變化,在汛期7~10月平均流量為23 919 m?3/s,非汛期平均流量為6 615 m?3/s。將水庫建設(shè)運行期和水庫混聯(lián)運行期各月徑流進行統(tǒng)計分析,得出宜昌站平均月流量變化趨勢,見圖4。

    宜昌站上游水庫混聯(lián)運行造成年極端流量相對于天然狀態(tài)流量變化幅度相對較弱,其中年最大流量降低幅度僅為5.4%,年最小流量增加了33.98%,說明水庫運行使河道徑流極值差減小,不利于構(gòu)建河流渠道地貌,對滯洪區(qū)與河道的交換養(yǎng)分有一定的阻礙作用。

    宜昌站水文要素改變度隨上游水庫數(shù)量增加而變大。在上游水庫混聯(lián)運行影響下,河道水文要素整體改變度為27.6%,屬于低度改變。在長江上游水庫群建設(shè)運行期,宜昌站整體水文情勢屬于低度改變。在1983~1998年,長江干支流上主要大型水庫,如三峽水庫、溪洛渡水庫、向家壩水庫、二灘水庫等均處于施工建設(shè)階段,對于河道徑流影響可忽略。嘉陵江、岷江、雅礱江等支流的徑流量均達不到宜昌站的20%,且各支流控制流域面積均為宜昌站的10%~15%,支流水庫的調(diào)度運行對徑流的影響,經(jīng)干支流交匯后對宜昌站的徑流影響程度就會顯著減弱。同時,全球氣候變化也是長江上游汛期徑流減少和枯季徑流增多的影響因素之一[17]。因此在建設(shè)運行期宜昌站徑流為低改變度是長江支流水庫運行和氣候變化共同引起的。

    此外,CRITIC計算結(jié)果表明,高場、北碚、宜昌站6~10月平均流量和最大月流量等指標所占權(quán)重較大,并且隨著水庫建設(shè)增加,其所占權(quán)重和均為70 %以上,表明其所含信息量較大,也說明水庫的建設(shè)運行主要影響汛期和最大的徑流過程。

    3.2.4整體改變度

    黎云云等人認為傳統(tǒng)RVA方法過于重視變化程度較大的指標參數(shù),而忽略了其他大多數(shù)屬于中低度改變的指標,這將在一定程度上給河流生態(tài)的恢復(fù)和管理帶來一定的盲目性和誤導(dǎo)性。

    本文以高場站水文要素14個指標的變化情況為例,采用CRITIC法計算整體指標改變度。龔嘴水庫單庫運行階段的整體改變度為13.7%;在龔嘴和銅街子水庫兩庫運行階段出現(xiàn)了7個指標的中度改變,7個指標所占權(quán)重為62.17%,整體改變度為34.3%;在新建紫坪鋪和瀑布溝水庫后,四庫并聯(lián)運行階段有4個指標呈高度改變,5個指標呈中度改變,其中高度改變指標權(quán)重占8.73%,小于中度改變度指標權(quán)重37.68%,因此中度改變指標對河道水文要素的影響較大。經(jīng)過改變度權(quán)重計算,得到高場站水文要素受上游水庫并聯(lián)運行影響整體改變度為38.2%,屬于中度改變。這是由于新建的紫坪鋪和瀑布溝水庫形成了并聯(lián)運行模式,使水庫群整體的調(diào)節(jié)能力加強,其防洪調(diào)蓄的作用在汛期尤其突出,大流量過程被攔蓄在水庫中,破壞了天然水流過程。從而引發(fā)水文指標的強烈變動。

    本文依據(jù)指標所含信息量大小來進行整體改變度計算,更加全面地衡量中低度改變指標的影響關(guān)系,在評價水文情勢改變度時更加客觀實際。

    4結(jié) 論

    本文耦合CRITIC法和RVA方法,解決了傳統(tǒng)RVA方法存在忽略中低度改變指標的局限性。根據(jù)長江上游水庫群建設(shè)的實際情況,選擇北碚、高場、宜昌站作為水庫串、并、混聯(lián)運行影響的代表站,通過對比分析14個水文指標在單庫、二庫情況下串并混聯(lián)合運行模式下的變化情況,深入剖析水庫群聯(lián)合調(diào)度對下游河道水文情勢的影響程度。主要結(jié)論如下。

    (1) 隨著長江上游干支流水庫建設(shè)數(shù)量增加,水庫群調(diào)蓄能力增強,對下游河道水文情勢的改變作用逐漸增大,水庫群下游各控制斷面月平均流量呈現(xiàn)非汛期流量增大、汛期流量減少的趨勢,全年流量過程更加平緩,說明水庫群的運行對汛期洪水有顯著的調(diào)節(jié)作用,能在一定程度上減輕河道防洪壓力。

    (2) 嘉陵江上游梯級水庫串聯(lián)運行,導(dǎo)致北碚站平均最大流量發(fā)生時間從7月延遲至9月,最大流量峰值顯著降低;影響包括最大、最小月平均流量在內(nèi)的6個水文指標出現(xiàn)高程度改變;綜合各水文指標所占權(quán)重,北碚站水文情勢整體改變度為40.5%,屬于中度改變。高場站上游水庫群從單庫運行模式發(fā)展到四庫并聯(lián)運行模式的過程,導(dǎo)致流量峰值坦化現(xiàn)象愈加明顯;高場站3個階段的水文指標改變程度逐步提高,整體改變度分別為13.7%,34.3%和38.2%,屬于中度改變。在長江上游水庫群建設(shè)運行前后,宜昌站整體水文情勢均屬于低度改變,在非汛期的流量變幅較為相近,汛期流量減幅為10.4%,宜昌站整體改變度為低程度。

    (3) 本文根據(jù)CRITIC法改進RVA方法,發(fā)現(xiàn)6~10月平均流量和最大月流量等指標所占權(quán)重較大,說明汛期流量和最大月流量的改變程度能主要反映河道水文情勢整體改變度。根據(jù)各指標權(quán)重反映控制斷面整體改變度,有效避免了傳統(tǒng)RVA中低度改變指標被忽略的現(xiàn)象,使整體評價結(jié)果更加科學(xué)合理,可為長江上游水利工程的生態(tài)調(diào)度模式提供參考。

    在本文研究分析過程中,發(fā)現(xiàn)典型斷面多年月平均徑流指標改變度有較大幅度的數(shù)據(jù)跳躍現(xiàn)象。針對這一現(xiàn)象,本文根據(jù)典型斷面上游水庫的調(diào)度規(guī)則進行解釋說明,但在后續(xù)研究中需要結(jié)合氣候變化因素和下墊面變化因素開展分析,更進一步探索河道水文情勢演變趨勢。

    參考文獻:

    [1]姚磊,陳盼盼,胡利利,等.長江上游流域水電開發(fā)現(xiàn)狀與存在的問題[J].綿陽師范學(xué)院學(xué)報,2016,35(2):91-97.

    [2]Hu W W,Wang G X,Deng W,et al.The influence of dams on ecohydrological conditions in the Huaihe River basin,China[J].Ecological Engineering,2008,33(3/4):233-241.

    [3]Gao B,Yang D W,Zhao T T G,et al.Changes in the eco-flowmetrics of the upper Yangtze River from 1961 to 2008[J].Journal of Hydrology,2012(448-449):30-38.

    [4]黎云云,暢建霞,涂歡,等.黃河干流控制性梯級水庫聯(lián)合運行對下游水文情勢的影響[J].資源科學(xué),2014,36(6):1183-1190.

    [5]段唯鑫,郭生練,王俊.長江上游大型水庫群對宜昌站水文情勢影響分析[J].長江流域資源與環(huán)境,2016,25(1):120-130.

    [6]蔡文君,殷峻暹,王浩.三峽水庫運行對長江中下游水文情勢的影響[J].人民長江,2012,43(5):22-25.

    [7]李林,王振宇,秦寧生,等.長江上游徑流量變化及其與影響因子關(guān)系分析[J].自然資源學(xué)報,2004(6):694-700.

    [8]Yin X A,Yang Z.F,Petts G E.A new method to assess the flow regime alterations in riverine ecosystems[J].River Research and Applications,2015,31(4):497-504.

    [9]周毅,崔同,高滿,等.考慮不同水文年及IHA指標相關(guān)性的水文特征評估方法[J].水文,2017,37(2):20-25.

    [10]程根偉,麻澤龍,范繼輝.西南江河梯級水電開發(fā)對河流水環(huán)境的影響及對策[J].中國科學(xué)院院刊,2004(6):433-437.

    [11]楊麗虎,陳進,黃薇,等.長江上游梯級水電開發(fā)情況及對生態(tài)環(huán)境的影響初探[J].中國農(nóng)村水利水電,2007(3):79-81.

    [12]金勇,周建軍,黃國鮮.長江上游大型水庫運行對長江上游水文過程的影響[J].水力發(fā)電學(xué)報,2010,29(2):94-101.

    [13]張建云,章四龍,王金星,等.近50年來中國六大流域年際徑流變化趨勢研究[J].水科學(xué)進展,2007(2):230-234.

    [14]Diakoulaki D,Mavrotas G,Papayannakis L.Determining objective weights in multiple criteria problems:the CRITIC method[J].Computers & Operations Research,1995,22(7):763-770.

    [15]袁合才,辛艷輝.基于AHP和CRITIC方法的水資源綜合效益模型[J].安徽農(nóng)業(yè)科學(xué),2011,39(4):2225-2226.

    [16]張玉,魏華波.基于CRITIC的多屬性決策組合賦權(quán)方法[J].統(tǒng)計與決策,2012(16):75-77.

    [17]趙軍凱,李九發(fā),戴志軍,等.長江宜昌站徑流變化過程分析[J].資源科學(xué),2012,34(12):2306-2315.

    引用本文:張康,楊明祥,梁藉,閆寶偉.長江上游水庫群聯(lián)合調(diào)度下的河流水文情勢研究[J].人民長江,2019,50(2):107-114.

    Study on river hydrology regime under joint operation of cascade reservoirs in upper Yangtze River

    ZHANG Kang?, YANG Mingxiang?2, LIANG Ji?, YAN Baowei

    (1.School of Hydropower and Information Engineering, Huazhong University of Science and Technology, Wuhan 430074, China;2.China Institute of Water Resources and Hydropower Research, Beijing 100044, China)

    Abstract: It is inevitable that joint operation of cascade reservoirs in the upper Yangtze River will exert certain effects on the runoffs in the lower river channel. We selected Beibei, Gaochang and Yichang stations as controlling cross-sections, and used IHA method to identify fourteen hydrological indexes that are influenced by the operation of reservoirs. We studied the varying trends of these indexes and the degree of hydrology variations under cascade reservoirs operation modes of series, parallel and series-parallel. The results show that: ① as the number of cascade reservoirs in operation increases, the variations in hydrology of mainstream Yangtze River tend to increase; ② joint operation of cascade reservoirs significantly decreases water flow during flood season; meanwhile, the maximum discharge throughout the year tends to decrease, and vice versa; ③series operation will delay maximum discharge and parallel operation will reinforce the evenness of discharge process during flood season; the variation of hydrological regime caused by series-parallel operation is mostly affected by the operation mode and climate change, series-parallel operation has low impact on the hydrological regime of Yichang station.

    Key words:overall variations of hydrological regime; CRITIC method; series, parallel and series-parallel operation of reservoirs; upper Yangtze River

    十八禁人妻一区二区| 高清毛片免费观看视频网站| 久久草成人影院| 亚洲 欧美一区二区三区| 成人特级黄色片久久久久久久| 美女午夜性视频免费| 色综合站精品国产| 亚洲人与动物交配视频| 丁香六月欧美| 97碰自拍视频| 成在线人永久免费视频| 成人性生交大片免费视频hd| 美女被艹到高潮喷水动态| 精品一区二区三区视频在线观看免费| 99精品欧美一区二区三区四区| 最近最新中文字幕大全免费视频| www国产在线视频色| 91九色精品人成在线观看| 国产乱人视频| 亚洲av电影在线进入| 欧美色视频一区免费| 国产精品 欧美亚洲| 免费看光身美女| 亚洲美女黄片视频| 香蕉丝袜av| 非洲黑人性xxxx精品又粗又长| 精品久久久久久久久久免费视频| 国产精品久久久久久久电影 | 国产亚洲欧美在线一区二区| 精品欧美国产一区二区三| 欧美中文综合在线视频| 欧美不卡视频在线免费观看| 丰满的人妻完整版| www.www免费av| 亚洲精品乱码久久久v下载方式 | 亚洲欧美日韩无卡精品| 俄罗斯特黄特色一大片| 欧美成人免费av一区二区三区| 韩国av一区二区三区四区| 神马国产精品三级电影在线观看| a级毛片a级免费在线| 熟女人妻精品中文字幕| 日韩欧美在线乱码| 国产成人系列免费观看| 欧美成人性av电影在线观看| 国产精品久久电影中文字幕| 变态另类成人亚洲欧美熟女| 久久久久久久久久黄片| 精品一区二区三区视频在线 | 日本黄大片高清| 亚洲,欧美精品.| 性色avwww在线观看| 黄色丝袜av网址大全| av福利片在线观看| 天堂动漫精品| 99久久国产精品久久久| 少妇的逼水好多| 又粗又爽又猛毛片免费看| 琪琪午夜伦伦电影理论片6080| 麻豆国产av国片精品| 午夜福利在线观看免费完整高清在 | 精品日产1卡2卡| 亚洲av片天天在线观看| 2021天堂中文幕一二区在线观| 午夜福利免费观看在线| 久久久久国产一级毛片高清牌| 免费看光身美女| 丁香欧美五月| 国产欧美日韩一区二区三| 久久婷婷人人爽人人干人人爱| 久久伊人香网站| 最新在线观看一区二区三区| 亚洲精品乱码久久久v下载方式 | 99国产综合亚洲精品| 免费搜索国产男女视频| 久久九九热精品免费| 国产成人欧美在线观看| 最近最新中文字幕大全电影3| 搡老熟女国产l中国老女人| 午夜两性在线视频| 久久人人精品亚洲av| 国产精品久久久人人做人人爽| 综合色av麻豆| 美女cb高潮喷水在线观看 | 欧美又色又爽又黄视频| 男女床上黄色一级片免费看| 国产av不卡久久| 色综合亚洲欧美另类图片| 波多野结衣巨乳人妻| 精品国产乱子伦一区二区三区| 观看免费一级毛片| 久久久久久久久免费视频了| 亚洲无线观看免费| 三级国产精品欧美在线观看 | 日本一二三区视频观看| 国产精品亚洲av一区麻豆| 欧洲精品卡2卡3卡4卡5卡区| 美女被艹到高潮喷水动态| 网址你懂的国产日韩在线| 国产乱人视频| 精品乱码久久久久久99久播| 制服丝袜大香蕉在线| 欧美成狂野欧美在线观看| www.精华液| 国产激情偷乱视频一区二区| 亚洲国产精品sss在线观看| cao死你这个sao货| 日韩大尺度精品在线看网址| 色在线成人网| 又大又爽又粗| 色综合婷婷激情| 亚洲美女黄片视频| 亚洲五月天丁香| 中文字幕最新亚洲高清| 亚洲成人久久性| 女生性感内裤真人,穿戴方法视频| 在线免费观看不下载黄p国产 | 国产精品av久久久久免费| 成年女人毛片免费观看观看9| 白带黄色成豆腐渣| 99精品欧美一区二区三区四区| 麻豆成人午夜福利视频| 一个人看视频在线观看www免费 | 色精品久久人妻99蜜桃| 免费在线观看日本一区| 女人高潮潮喷娇喘18禁视频| 国产精品 国内视频| 黑人欧美特级aaaaaa片| 精品电影一区二区在线| cao死你这个sao货| 久久午夜综合久久蜜桃| 亚洲真实伦在线观看| а√天堂www在线а√下载| 国产精品久久电影中文字幕| av福利片在线观看| a级毛片a级免费在线| 精品久久蜜臀av无| 身体一侧抽搐| av在线天堂中文字幕| 校园春色视频在线观看| 黄色 视频免费看| 天天添夜夜摸| 国产 一区 欧美 日韩| 久久久久久九九精品二区国产| 久久中文字幕人妻熟女| 国内久久婷婷六月综合欲色啪| 99视频精品全部免费 在线 | 国产精品国产高清国产av| 午夜a级毛片| 国产精华一区二区三区| 久久精品亚洲精品国产色婷小说| 深夜精品福利| 国产v大片淫在线免费观看| 九九在线视频观看精品| 在线观看美女被高潮喷水网站 | 欧美成人一区二区免费高清观看 | 男女之事视频高清在线观看| 男女那种视频在线观看| 日韩成人在线观看一区二区三区| 国产亚洲精品久久久com| 99久久精品一区二区三区| 国产高清视频在线播放一区| 天天躁日日操中文字幕| 亚洲av成人一区二区三| 久久精品aⅴ一区二区三区四区| 久久性视频一级片| 久久精品国产综合久久久| 一区福利在线观看| 欧美日韩一级在线毛片| 757午夜福利合集在线观看| 宅男免费午夜| 欧美丝袜亚洲另类 | 99riav亚洲国产免费| 窝窝影院91人妻| 丁香欧美五月| 女警被强在线播放| 在线观看免费视频日本深夜| 久久中文字幕一级| 久久久精品大字幕| av天堂中文字幕网| 99热只有精品国产| 国产视频内射| 国产一区二区三区视频了| 欧美性猛交╳xxx乱大交人| 亚洲男人的天堂狠狠| 99视频精品全部免费 在线 | 免费无遮挡裸体视频| 97人妻精品一区二区三区麻豆| 午夜免费激情av| 欧美高清成人免费视频www| 在线观看一区二区三区| 非洲黑人性xxxx精品又粗又长| 激情在线观看视频在线高清| 老汉色av国产亚洲站长工具| 午夜成年电影在线免费观看| 亚洲av成人不卡在线观看播放网| 成人特级av手机在线观看| 免费看十八禁软件| 99在线视频只有这里精品首页| 伦理电影免费视频| 婷婷亚洲欧美| 99久久精品热视频| 亚洲av熟女| av片东京热男人的天堂| 国产aⅴ精品一区二区三区波| 亚洲色图av天堂| 亚洲国产看品久久| 亚洲欧美日韩卡通动漫| 午夜福利成人在线免费观看| 欧美乱色亚洲激情| 日韩欧美精品v在线| 制服丝袜大香蕉在线| 性欧美人与动物交配| 国产伦一二天堂av在线观看| 国产av在哪里看| 欧美一区二区国产精品久久精品| 成人性生交大片免费视频hd| 午夜精品一区二区三区免费看| 久久人人精品亚洲av| 嫁个100分男人电影在线观看| 久久久久久久久中文| 国产蜜桃级精品一区二区三区| 日韩有码中文字幕| 亚洲 欧美 日韩 在线 免费| 亚洲精品中文字幕一二三四区| 国产探花在线观看一区二区| 日本 av在线| 午夜激情福利司机影院| 国产一级毛片七仙女欲春2| 一夜夜www| 成人一区二区视频在线观看| 午夜免费观看网址| 又爽又黄无遮挡网站| 身体一侧抽搐| 亚洲精品久久国产高清桃花| 日韩精品中文字幕看吧| 亚洲九九香蕉| 看片在线看免费视频| 久久99热这里只有精品18| 免费在线观看成人毛片| 欧美三级亚洲精品| 天天躁日日操中文字幕| 91九色精品人成在线观看| 岛国视频午夜一区免费看| 国产精品亚洲av一区麻豆| 国内精品美女久久久久久| 一夜夜www| 国产精品电影一区二区三区| 99热只有精品国产| 国产av一区在线观看免费| 欧美又色又爽又黄视频| 亚洲五月婷婷丁香| 亚洲最大成人中文| 久久久国产精品麻豆| 国产成人啪精品午夜网站| 男人舔奶头视频| 桃红色精品国产亚洲av| 国语自产精品视频在线第100页| 精品久久久久久久毛片微露脸| 国产综合懂色| 热99在线观看视频| 91麻豆av在线| 色在线成人网| 欧美色欧美亚洲另类二区| 成年女人看的毛片在线观看| 国产精品久久久久久人妻精品电影| 少妇的逼水好多| 18禁国产床啪视频网站| 欧美国产日韩亚洲一区| 一级a爱片免费观看的视频| 久99久视频精品免费| 精品福利观看| 香蕉丝袜av| 99热精品在线国产| 美女高潮喷水抽搐中文字幕| www.精华液| 亚洲真实伦在线观看| 香蕉国产在线看| 一个人看的www免费观看视频| 亚洲av电影在线进入| 国产高清三级在线| 草草在线视频免费看| 桃色一区二区三区在线观看| 国产一区二区三区视频了| 嫩草影院入口| 国产精品99久久久久久久久| 一个人免费在线观看的高清视频| 免费一级毛片在线播放高清视频| 久久精品影院6| 一本精品99久久精品77| 亚洲国产精品sss在线观看| 国产精品 国内视频| 又黄又粗又硬又大视频| 久久久精品欧美日韩精品| 国内精品久久久久精免费| 免费av毛片视频| 成人三级做爰电影| 国产一区二区在线观看日韩 | 国产单亲对白刺激| 美女大奶头视频| 在线看三级毛片| 老熟妇仑乱视频hdxx| 国产精华一区二区三区| 亚洲国产高清在线一区二区三| 夜夜爽天天搞| 免费无遮挡裸体视频| 日本免费a在线| 亚洲国产欧美网| 亚洲精品国产精品久久久不卡| xxx96com| 日本精品一区二区三区蜜桃| 亚洲av成人精品一区久久| 国产蜜桃级精品一区二区三区| 特大巨黑吊av在线直播| 午夜福利欧美成人| 精品99又大又爽又粗少妇毛片 | 国产高清三级在线| 亚洲av免费在线观看| 伦理电影免费视频| 亚洲激情在线av| 国产精品 欧美亚洲| 中文字幕人成人乱码亚洲影| 一卡2卡三卡四卡精品乱码亚洲| 亚洲自拍偷在线| 国产视频内射| 亚洲美女黄片视频| 少妇的丰满在线观看| av中文乱码字幕在线| 三级毛片av免费| 国产高清视频在线播放一区| 操出白浆在线播放| 日韩精品青青久久久久久| 在线观看免费午夜福利视频| 身体一侧抽搐| 国产真实乱freesex| 亚洲中文字幕一区二区三区有码在线看 | 成人欧美大片| 视频区欧美日本亚洲| 丁香六月欧美| 亚洲五月天丁香| 欧美日韩国产亚洲二区| 久久中文字幕一级| 日韩国内少妇激情av| 日韩欧美免费精品| 日本熟妇午夜| 亚洲黑人精品在线| 久久国产精品影院| 一本久久中文字幕| 国内精品美女久久久久久| 禁无遮挡网站| 每晚都被弄得嗷嗷叫到高潮| 一进一出抽搐动态| 精品人妻1区二区| 亚洲人成伊人成综合网2020| 国产野战对白在线观看| 老司机午夜福利在线观看视频| 久久久国产欧美日韩av| 97超级碰碰碰精品色视频在线观看| 熟女电影av网| 久久久精品欧美日韩精品| 国产激情久久老熟女| 欧美日韩综合久久久久久 | 免费在线观看影片大全网站| 久久香蕉国产精品| 99热这里只有精品一区 | 精品久久久久久久毛片微露脸| 精品久久久久久久久久久久久| 国产精品久久久久久亚洲av鲁大| 老鸭窝网址在线观看| 亚洲成av人片在线播放无| 免费电影在线观看免费观看| 国产97色在线日韩免费| 成年女人看的毛片在线观看| а√天堂www在线а√下载| 国产探花在线观看一区二区| 午夜a级毛片| 免费在线观看成人毛片| 国产极品精品免费视频能看的| 最近最新免费中文字幕在线| 亚洲国产高清在线一区二区三| 免费观看的影片在线观看| 国产极品精品免费视频能看的| 成年人黄色毛片网站| 999精品在线视频| 成人av在线播放网站| 男人舔女人的私密视频| 最新在线观看一区二区三区| 国产一区二区在线观看日韩 | 国产激情欧美一区二区| 亚洲狠狠婷婷综合久久图片| 亚洲va日本ⅴa欧美va伊人久久| 亚洲七黄色美女视频| 欧美不卡视频在线免费观看| 性色avwww在线观看| 久久精品91蜜桃| 1024手机看黄色片| 国产精品久久久久久亚洲av鲁大| 狂野欧美激情性xxxx| 五月伊人婷婷丁香| 国产久久久一区二区三区| 欧美午夜高清在线| 亚洲 欧美一区二区三区| 色综合站精品国产| 欧美绝顶高潮抽搐喷水| 日本 欧美在线| 久久香蕉精品热| 国产亚洲精品久久久com| 12—13女人毛片做爰片一| 国产成人av教育| 2021天堂中文幕一二区在线观| www国产在线视频色| 一a级毛片在线观看| 欧美国产日韩亚洲一区| 亚洲国产欧洲综合997久久,| 日韩欧美在线乱码| 日日干狠狠操夜夜爽| 国产精品一区二区精品视频观看| 国产亚洲精品久久久com| 一区二区三区高清视频在线| 成年免费大片在线观看| 亚洲美女视频黄频| 久久久精品大字幕| or卡值多少钱| 婷婷亚洲欧美| 给我免费播放毛片高清在线观看| 大型黄色视频在线免费观看| 久久精品国产99精品国产亚洲性色| 午夜福利欧美成人| 真实男女啪啪啪动态图| 欧美日韩中文字幕国产精品一区二区三区| 又黄又粗又硬又大视频| 欧美成人一区二区免费高清观看 | 最新在线观看一区二区三区| 久久亚洲精品不卡| 两个人看的免费小视频| 90打野战视频偷拍视频| 日韩精品中文字幕看吧| 午夜福利视频1000在线观看| 日韩精品青青久久久久久| 欧美中文日本在线观看视频| 在线永久观看黄色视频| 成人欧美大片| 一夜夜www| 变态另类成人亚洲欧美熟女| 国产真实乱freesex| av中文乱码字幕在线| 成人三级做爰电影| 亚洲av免费在线观看| 麻豆一二三区av精品| 日本免费a在线| 亚洲精品一卡2卡三卡4卡5卡| 久久性视频一级片| 99国产极品粉嫩在线观看| 久久久久国产一级毛片高清牌| 国产蜜桃级精品一区二区三区| 国产综合懂色| 午夜福利免费观看在线| 亚洲午夜精品一区,二区,三区| 俺也久久电影网| 午夜福利高清视频| 99久国产av精品| 亚洲激情在线av| 成人欧美大片| 亚洲七黄色美女视频| 国产亚洲av高清不卡| 欧美日韩福利视频一区二区| 搡老岳熟女国产| 午夜久久久久精精品| 女生性感内裤真人,穿戴方法视频| 国产伦一二天堂av在线观看| 18美女黄网站色大片免费观看| 99久久久亚洲精品蜜臀av| 国产亚洲欧美98| 男女做爰动态图高潮gif福利片| 国产成人av激情在线播放| 亚洲狠狠婷婷综合久久图片| 国产精品99久久99久久久不卡| 久久久色成人| 日韩欧美精品v在线| 在线播放国产精品三级| 欧美一区二区国产精品久久精品| 又黄又爽又免费观看的视频| 久久午夜亚洲精品久久| 亚洲自拍偷在线| 高清在线国产一区| 欧美日韩中文字幕国产精品一区二区三区| 欧美日韩瑟瑟在线播放| 亚洲色图 男人天堂 中文字幕| 亚洲真实伦在线观看| 国产精品 欧美亚洲| 日韩三级视频一区二区三区| 在线十欧美十亚洲十日本专区| 桃红色精品国产亚洲av| 人人妻人人澡欧美一区二区| 久久久久精品国产欧美久久久| 久久香蕉国产精品| 婷婷亚洲欧美| 在线观看日韩欧美| 日韩 欧美 亚洲 中文字幕| 亚洲av成人精品一区久久| 国产精品自产拍在线观看55亚洲| 在线观看66精品国产| 一进一出抽搐动态| 国内揄拍国产精品人妻在线| 美女被艹到高潮喷水动态| 国产高清有码在线观看视频| 人妻夜夜爽99麻豆av| www国产在线视频色| 在线免费观看的www视频| 国产亚洲精品av在线| 麻豆av在线久日| 久久精品人妻少妇| 日本三级黄在线观看| 国产亚洲欧美在线一区二区| 亚洲色图av天堂| 久99久视频精品免费| 日韩 欧美 亚洲 中文字幕| 精品国产乱子伦一区二区三区| 免费大片18禁| 欧美三级亚洲精品| 性欧美人与动物交配| 久久久久久久久久黄片| 无限看片的www在线观看| 日本一二三区视频观看| 亚洲自拍偷在线| 亚洲人成网站高清观看| 日本免费a在线| 国产99白浆流出| 久久久国产成人精品二区| 亚洲av免费在线观看| 久久精品亚洲精品国产色婷小说| 黄片大片在线免费观看| 村上凉子中文字幕在线| 久久久久久久久免费视频了| 一夜夜www| 老汉色∧v一级毛片| 欧美一级毛片孕妇| 变态另类丝袜制服| 成年版毛片免费区| 三级毛片av免费| 日韩欧美一区二区三区在线观看| 国产亚洲精品综合一区在线观看| 18禁裸乳无遮挡免费网站照片| 人人妻人人澡欧美一区二区| 男女那种视频在线观看| 一夜夜www| 每晚都被弄得嗷嗷叫到高潮| 美女扒开内裤让男人捅视频| 成人特级av手机在线观看| 免费看美女性在线毛片视频| 日本黄色视频三级网站网址| 狂野欧美激情性xxxx| 免费高清视频大片| 手机成人av网站| 午夜激情欧美在线| 国产91精品成人一区二区三区| 亚洲一区二区三区不卡视频| 精品欧美国产一区二区三| 国产一区二区在线观看日韩 | 欧美激情久久久久久爽电影| 国产成人aa在线观看| 久久99热这里只有精品18| netflix在线观看网站| 久久久久九九精品影院| 亚洲国产精品sss在线观看| 免费在线观看日本一区| 夜夜爽天天搞| 精品人妻1区二区| 亚洲人成网站高清观看| 免费看美女性在线毛片视频| 国产毛片a区久久久久| 又大又爽又粗| 999精品在线视频| 亚洲色图 男人天堂 中文字幕| 国产高清视频在线播放一区| 精品熟女少妇八av免费久了| 国产乱人视频| 性欧美人与动物交配| 午夜a级毛片| 久久午夜综合久久蜜桃| 床上黄色一级片| 好看av亚洲va欧美ⅴa在| 久久久久久久精品吃奶| 国产精品久久久久久久电影 | 欧美3d第一页| 99久久精品热视频| 精品国产乱子伦一区二区三区| 久久精品综合一区二区三区| 国产免费男女视频| 麻豆av在线久日| 999精品在线视频| 国产视频内射| 午夜a级毛片| 午夜久久久久精精品| 变态另类丝袜制服| 国产亚洲精品av在线| 母亲3免费完整高清在线观看| 国产一级毛片七仙女欲春2| 国产真实乱freesex| 女人高潮潮喷娇喘18禁视频| 老司机福利观看| 国产91精品成人一区二区三区| 99热精品在线国产| 免费观看人在逋| 婷婷亚洲欧美| 久久这里只有精品19| 中文字幕人妻丝袜一区二区| 国产黄片美女视频| 成人鲁丝片一二三区免费| 国产一区二区激情短视频| 亚洲成人免费电影在线观看| 亚洲av电影不卡..在线观看| 亚洲狠狠婷婷综合久久图片| av在线蜜桃| 国产高清videossex|