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

    嘉陵江流域水系形態(tài)特征與侵蝕產(chǎn)沙響應(yīng)

    2016-10-26 02:44:25李進(jìn)林
    水土保持研究 2016年5期
    關(guān)鍵詞:輸沙量水系河段

    李進(jìn)林,韋 杰,2

    (1.重慶師范大學(xué) 地理與旅游學(xué)院,重慶 401331; 2.三峽庫區(qū)地表過程與環(huán)境遙感重慶市重點(diǎn)實(shí)驗(yàn)室,重慶 401331)

    ?

    嘉陵江流域水系形態(tài)特征與侵蝕產(chǎn)沙響應(yīng)

    李進(jìn)林1,韋 杰1,2

    (1.重慶師范大學(xué) 地理與旅游學(xué)院,重慶 401331; 2.三峽庫區(qū)地表過程與環(huán)境遙感重慶市重點(diǎn)實(shí)驗(yàn)室,重慶 401331)

    流域水系形態(tài)指示著流域地貌發(fā)育和水文過程,研究其形態(tài)特征有助于流域水力開發(fā)、生態(tài)保護(hù)、水土流失治理等規(guī)劃的制定?;贒EM數(shù)據(jù),運(yùn)用分形維數(shù)、河段蜿蜒度、灰色關(guān)聯(lián)度和非線性響應(yīng)模型分析了嘉陵江流域水系形態(tài)特征以及侵蝕產(chǎn)沙與水系形態(tài)的關(guān)系。結(jié)果表明:嘉陵江流域內(nèi)各子流域水系形態(tài)差異顯著,其中,渠江子流域水系形態(tài)分形維數(shù)為5.102,水系分叉總數(shù)增加最快,分枝形態(tài)最復(fù)雜,溯源侵蝕能力最強(qiáng)。其余4個(gè)子流域水系形態(tài)分形維數(shù)分別為白龍江子流域3.706,干流上游子流域3.605,干流中下游子流域3.576,涪江子流域3.548。各子流域河段均以低蜿蜒度為主,約占河段總數(shù)的73%~82%;中蜿蜒度河段數(shù)約占總數(shù)的12%~20%。平均輸沙量和泥沙淤積量與水系分形維數(shù)間的關(guān)聯(lián)度分別為0.618,0.641,與河段平均蜿蜒度的關(guān)聯(lián)度分別為0.588,0.591。河段蜿蜒度和分枝結(jié)構(gòu)復(fù)雜度對(duì)輸沙量的相對(duì)貢獻(xiàn)率分別為21%和79%,對(duì)泥沙淤積量的相對(duì)貢獻(xiàn)率分別為18%和82%。

    分形維數(shù); 河段蜿蜒度; 灰色關(guān)聯(lián)度; 相對(duì)貢獻(xiàn)率; 侵蝕產(chǎn)沙

    水系是流域內(nèi)不同等級(jí)河流所構(gòu)成的河道系統(tǒng),是陸地水循環(huán)的主要路徑和活躍的外營力[1]。水系發(fā)育主要受流域地質(zhì)構(gòu)造、地形地貌、巖性、土壤、氣候、植被條件等因子的影響[1-4],反過來,水系又影響著流域水文過程、地表形態(tài)、土壤侵蝕、植被和氣候形成乃至區(qū)域發(fā)展[4-5]。水系形態(tài)特征是決定流域水生態(tài)與水環(huán)境的基礎(chǔ)之一[6],準(zhǔn)確把握水系形態(tài)特征對(duì)理解流域地貌發(fā)育、水文過程和指導(dǎo)制定流域水力開發(fā)、水土流失治理規(guī)劃等具有重要意義[7],但水系具有結(jié)構(gòu)多樣性和空間組合的復(fù)雜性,給科學(xué)精確地量化流域水系形態(tài)帶來了困難。20世紀(jì)70年代美國數(shù)學(xué)家B.Mandelbrot[8]建立了分形理論并將該理論引入水文地理學(xué)研究中,采用分形維數(shù)描述流域水系尺度不變性和形態(tài)的復(fù)雜性[9],為水系形態(tài)特征量化提供了有效工具。隨后該理論進(jìn)一步發(fā)展和完善,Tarboton[10]、Robert[11]、Rosso[12]等推演出了一系列流域水系分形維數(shù)計(jì)算方法。

    常用于表征水系形態(tài)的分形維數(shù)主要有計(jì)盒維數(shù)、霍頓(Horton)比維數(shù)和分枝維數(shù)3種。其中,計(jì)盒維數(shù)能從總體上刻畫水系的復(fù)雜程度[13],該方法因計(jì)算過程簡單和物理意義明確而被廣泛應(yīng)用[14];但計(jì)盒維數(shù)主要刻畫河流的彎曲復(fù)雜特征,無法體現(xiàn)特定流域內(nèi)水系分枝形態(tài)的復(fù)雜性?;纛D比維數(shù)能從一定程度上體現(xiàn)特定流域水系形態(tài)的復(fù)雜程度,但其前提是流域水系要滿足分形樹結(jié)構(gòu)時(shí),其結(jié)果才具有一定的可靠性[15-16],而完全滿足這一前提的自然水系比較少見。Tarboton等[10]應(yīng)用霍頓比法計(jì)算出水系分形維數(shù)高于上限值也證實(shí)了該方法的可靠性不高。分枝維數(shù)廣泛應(yīng)用于交通網(wǎng)絡(luò)、城鎮(zhèn)體系、區(qū)域水系等空間結(jié)構(gòu)研究[4,17],側(cè)重體現(xiàn)研究對(duì)象間的集散狀態(tài)和相對(duì)位置的空間變化關(guān)系。該方法同樣不能很好地刻畫流域不同級(jí)別水系分枝形態(tài)特征,主要原因是水系不存在縱橫交錯(cuò)的結(jié)構(gòu)特征,只是從各發(fā)源地向流域出水口方向逐級(jí)匯流合并,同時(shí),流域水系“樞紐點(diǎn)”的確定也缺乏科學(xué)依據(jù)。

    目前關(guān)于嘉陵江流域的研究主要集中在水土流失及治理、面源污染及控制、水文過程等方面[18-22],而關(guān)于水系形態(tài)特征及其與侵蝕產(chǎn)沙的關(guān)系很少報(bào)道。本文基于DEM數(shù)據(jù)和ArcGIS軟件平臺(tái),利用河段蜿蜒度和改進(jìn)的水系分形維數(shù)探討嘉陵江流域水系形態(tài)特征,并利用灰色關(guān)聯(lián)度和非線性響應(yīng)模型分析水系形態(tài)特征與流域侵蝕產(chǎn)沙量和泥沙淤積量間的關(guān)系,為認(rèn)識(shí)該流域水系形態(tài)和地貌發(fā)育、制定水土流失治理規(guī)劃等提供依據(jù)。

    1 研究區(qū)概況

    嘉陵江流域(29°17′30″—34°28′11″N,102°35′36″—109°01′08″E)東北部以秦巴山脈與漢水為界,東南部以華鎣山與長江相隔,西部及西南部與沱江接壤,西北部以龍門山與岷江相連。流域面積約16萬km2,包括渠江、涪江和嘉陵江干流三大水系,各水系上游均為河床比降大、河谷狹窄的山區(qū),下游出水口在重慶市合川區(qū)附近匯合,構(gòu)成扇形向心水系。流域橫跨三大構(gòu)造單元,地質(zhì)構(gòu)造復(fù)雜,地貌類型多樣,地層分布主要是中生界侏羅系紅色巖層。土壤類型主要有紫色土、黃壤土、沖積土、水稻土、潮土和紅壤。干流長約1 120 km,發(fā)源于陜西省鳳縣秦嶺南麓,向西南流經(jīng)陜、甘、川、渝后匯入長江,是長江上游流域面積最大的一條支流。流域內(nèi)大部分地區(qū)屬于亞熱帶濕潤季風(fēng)氣候,多年平均降雨量和蒸發(fā)量分別為931 mm和513 mm,年內(nèi)降水主要集中在5—10月,占全年降水量的82%,多年平均最高氣溫和最低氣溫分別為19.4℃和4.3℃,多年平均日照時(shí)間為1 450 h[20]。第三次土壤侵蝕遙感調(diào)查結(jié)果顯示,嘉陵江流域水土流失面積51 297.45 km2,占流域總面積的49.24%,其中,輕度、中度、強(qiáng)度、極強(qiáng)度和劇烈水土流失面積分別占流域水土流失面積的32.26%,47.37%,16.36%,3.12%和0.88%。流域多年平均懸移質(zhì)輸沙量和卵石推移質(zhì)輸沙量分別為1.11億t和7萬t,是三峽水庫泥沙的重要來源之一。自1989年實(shí)施“長治”工程以來,流域內(nèi)累計(jì)水土流失治理面積已經(jīng)超過1.65萬km2,水土流失已經(jīng)得到初步控制[22]。

    2 資料與方法

    2.1數(shù)據(jù)來源與處理

    采用的ASTER GDEM V2 (ASTER Global Digital Elevation Model V2 )數(shù)據(jù)來源于中國科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國際科學(xué)數(shù)據(jù)鏡像網(wǎng)站(http:∥www.gscloud.cn),空間分辨率為30 m??紤]到流域不同部位地貌特征的顯著差異,利用ArcGIS軟件的水文分析工具將嘉陵江流域劃分為白龍江子流域、干流上游子流域、涪江子流域、渠江子流域和干流中下游子流域5個(gè)子流域(圖1)。為使各子流域水系結(jié)構(gòu)復(fù)雜程度存在可比性,從DEM數(shù)據(jù)提取各子流域水系時(shí),將匯水量閾值統(tǒng)一設(shè)置為5 000。水系提取完成后,依據(jù)Strahler分級(jí)法將水系河段分級(jí),該分級(jí)法是將流域內(nèi)所有無支流的河段定為第1級(jí)水系,由兩條1級(jí)河段匯流形成的河段定為第2級(jí)水系,以此類推,一直到流域出水口[16];在這種分級(jí)法中,當(dāng)且僅當(dāng)兩條相同級(jí)別的河段匯流成一條時(shí),該河段級(jí)別才會(huì)增加,對(duì)于低級(jí)別河段匯入高級(jí)別河段時(shí),高級(jí)別河段的級(jí)別不會(huì)改變。最終,白龍江流域和渠江流域分為7級(jí),其余3個(gè)流域分為6級(jí);完成分級(jí)后刪除了河網(wǎng)中的偽溝谷和孤河。

    圖1 嘉陵江流域區(qū)位圖及子流域劃分

    多年平均輸沙量和泥沙輸移比是采用馬煉[18]、許全喜[19]和李繼承[20]的部分研究結(jié)果,數(shù)據(jù)來源于研究區(qū)內(nèi)相關(guān)水文站1954—1998年的水沙監(jiān)測(cè)資料。

    2.2水系形態(tài)分形維數(shù)

    采用的水系形態(tài)分形維數(shù)計(jì)算方法是由計(jì)盒維數(shù)改進(jìn)得到,具體的做法是將計(jì)盒維數(shù)中的盒子尺寸改為水系級(jí)別、所對(duì)應(yīng)的非空盒子數(shù)目改為累計(jì)(向高級(jí)別水系累加至最高級(jí)別水系)水系河段數(shù);隨著流域水系級(jí)別的不斷變化,可獲得累計(jì)河段數(shù)目序列,水系級(jí)別序列與累計(jì)河段數(shù)目序列滿足式(1):

    Nr∝r-D

    (1)

    (2)

    式(1)兩邊同時(shí)取對(duì)數(shù)可得式(3):

    lgNr∝-Dlgr

    (3)

    以(lgr,lgNr)為坐標(biāo)系,用(r1,r2,r3,…,)序列和(N1,N2,N3,…,)序列繪制雙對(duì)數(shù)圖,并利用最小二乘法擬合直線(4)。

    lgNr=-Dlgr+A

    (4)

    式(1)—(4)中:r為水系級(jí)別(r =1,2,3,…,n),n為最高級(jí)別水系,Ni為第i級(jí)水系河段數(shù),Nr為r=i級(jí)及其以上級(jí)別水系累計(jì)河段數(shù),A為常數(shù),D為分形維數(shù),其本質(zhì)是由水系河流分枝累計(jì)數(shù)目變化所決定;D值越高,表示水系從流域出水口向各發(fā)源地分叉總數(shù)增加相對(duì)較快,即水系分叉發(fā)育能力和溯源侵蝕能力越強(qiáng)、結(jié)構(gòu)越復(fù)雜;反之則表示水系分叉總數(shù)增加相對(duì)較慢,水系分叉發(fā)育能力和溯源侵蝕能力較弱、結(jié)構(gòu)較簡單。

    2.3河段蜿蜒度

    河段蜿蜒度是指河流中心線長度與河流所處流域中心線長度的比值,計(jì)算方法為式(5):

    S=Lr/Lv

    (5)

    式中:S為河流蜿蜒度;Lr為河流中心線長度,即河段本身的長度;Lv為河流所在流域中心線長度,即河段起止點(diǎn)間的直線距離。

    2.4灰色關(guān)聯(lián)度

    灰色關(guān)聯(lián)分析是根據(jù)兩個(gè)因素的時(shí)間序列發(fā)展趨勢(shì),衡量因素間關(guān)聯(lián)程度的一種方法。由于流域多年平均輸沙量、泥沙淤積量和水系形態(tài)參數(shù)在短時(shí)間尺度內(nèi)不會(huì)發(fā)生變化,本文依據(jù)“時(shí)空對(duì)換”思想,將上述“狀態(tài)數(shù)據(jù)”轉(zhuǎn)化為“過程數(shù)據(jù)”,并以平均輸沙量和淤積量為參考序列,分別計(jì)算二者與比較序列(水系分形維數(shù)和平均蜿蜒度)間的關(guān)聯(lián)系數(shù):

    (6)

    關(guān)聯(lián)系數(shù)是某時(shí)刻比較序列與參考序列的關(guān)聯(lián)程度值,信息過于分散不便于整體比較。因此,將其平均值作為灰色關(guān)聯(lián)度用于衡量因素間的整體關(guān)聯(lián)程度,計(jì)算方法為式(7):

    (7)

    2.5相對(duì)貢獻(xiàn)率分析

    輸沙量和泥沙淤積量對(duì)水系形態(tài)特征的響應(yīng)關(guān)系可以表示為:

    Syi=Wa·Db

    (8)

    Sde=Wc·Dd

    (9)

    式中:Syi表示流域侵蝕產(chǎn)沙量;Sde表示泥沙淤積量;W和D為河段蜿蜒度和水系分形維數(shù),a,b,c和d為響應(yīng)系數(shù),由非線性擬合得到。

    相對(duì)貢獻(xiàn)率可定量描述水系形態(tài)特征對(duì)流域侵蝕產(chǎn)沙和泥沙淤積變化的作用大小,計(jì)算方法為:

    (10)

    式中:ηi為相對(duì)貢獻(xiàn)率(%);αi為非線性響應(yīng)方程中因子指數(shù)絕對(duì)值。

    3 結(jié)果與分析

    3.1嘉陵江各子流域水系形態(tài)分形

    各子流域水系級(jí)別與累計(jì)河段數(shù)目間存在非線性關(guān)系,采用指數(shù)曲線擬合時(shí)(圖2左列),決定系數(shù)R2值均為0.99;從水系級(jí)別與累計(jì)河段數(shù)目的雙對(duì)數(shù)圖(圖2右列)可看出,二者存在較強(qiáng)對(duì)數(shù)線性相關(guān),決定系數(shù)R2值介于0.94到0.97之間;表明水系分枝結(jié)構(gòu)存在分形特征,其他流域的研究結(jié)果也有這樣的特征[16]。

    嘉陵江各子流域水系形態(tài)分形維數(shù)(圖2右列)從大到小依次為渠江子流域5.102,白龍江子流域3.706,干流上游子流域3.605,干流中下游子流域3.576和涪江子流域3.548,表明嘉陵江流域內(nèi)各子流域水系分枝形態(tài)存在差異,渠江子流域水系分叉總數(shù)增加最快,水系分枝形態(tài)最復(fù)雜,溯源侵蝕能力最強(qiáng)。其余4個(gè)子流域間水系分形維數(shù)差異相對(duì)較小,水系分叉發(fā)育能力相當(dāng),水系分枝形態(tài)復(fù)雜程度相似,溯源侵蝕能力接近,原因可能是水系結(jié)構(gòu)發(fā)育過程主要受構(gòu)造運(yùn)動(dòng)、地形地貌、下覆巖性和氣候特征的影響,其中構(gòu)造運(yùn)動(dòng)起控制作用[23]。渠江子流域上游處于大巴山構(gòu)造帶和米蒼山構(gòu)造帶(位于揚(yáng)子地塊北緣、秦嶺構(gòu)造帶和松潘—甘孜褶皺帶間的過渡帶),主要經(jīng)歷了加里東、華力西、印支、燕山和喜馬拉雅運(yùn)動(dòng),地貌條件十分復(fù)雜,中下游位于揚(yáng)子地塊上部的川中盆地丘陵區(qū)和川東平行嶺谷區(qū),雖基底剛硬,但上覆抗侵蝕能力較弱的侏羅系泥巖、頁巖和粉砂巖;再加上區(qū)內(nèi)降雨豐沛,年降雨量為1 014~1 500 mm,使渠江子流域水系結(jié)構(gòu)發(fā)育相對(duì)復(fù)雜。白龍江子流域位于揚(yáng)子地塊、西藏—三江造山帶與秦祁昆造山帶的交匯區(qū),干流上游子流域主要位于西秦嶺構(gòu)造帶,二者同樣經(jīng)歷了多期構(gòu)造演化,但是水系結(jié)構(gòu)復(fù)雜程度與渠江子流域明顯不同;可能是由氣候特征和植被條件的不同所造成,這兩種影響因素作用強(qiáng)度相對(duì)構(gòu)造運(yùn)動(dòng)非常微弱,但長時(shí)間尺度的累計(jì)效應(yīng)卻不可忽視。干流中下游子流域和涪江子流域的氣候特征與植被條件與渠江流域相似,水系復(fù)雜程度與渠江子流域同樣存在差異,主要原因可能是地貌條件和受構(gòu)造運(yùn)動(dòng)的影響不同,前兩者主要位于揚(yáng)子地塊上部的川中盆地丘陵區(qū),水系演化過程主要受揚(yáng)子板塊活動(dòng)的影響[24-26]。

    3.2嘉陵江流域水系河段蜿蜒度

    根據(jù)式(5)計(jì)算出5個(gè)子流域共11 648條河段的蜿蜒度,并參照Rosgen[27]的分類方法對(duì)蜿蜒度進(jìn)行分類,即:S<1.2為低蜿蜒度,1.2≤S≤1.4為中蜿蜒度,S >1.4為高蜿蜒度,結(jié)果見表1。

    子流域河段以低蜿蜒度為主,占河段總數(shù)的比例分別為白龍江子流域78.06%、干流上游子流域82.40%、涪江子流域73.40%、干流中下游子流域73.18%和渠江子流域75.03%,最小蜿蜒度值均接近1。各子流域中蜿蜒度河段數(shù)由多到少依次為渠江子流域544條、涪江子流域481條、干流中下游子流域375條、白龍江子流域290條和干流上游子流域280條,分別占河段總數(shù)的19.35%,19.71%,19.64%,11.85%和13.73%。高蜿蜒度河段數(shù)最多的是白龍江子流域,247條,最少的是干流上游子流域,79條,二者相差168條。進(jìn)一步分析發(fā)現(xiàn),低蜿蜒度河段主要集中于各子流域1級(jí)水系中,原因可能是1級(jí)水系主要分布在山區(qū),河段蜿蜒擺動(dòng)受到限制;除白龍江子流域外,其余各子流域最高級(jí)別水系河段均為高蜿蜒度,其中干流中下游子流域和渠江子流域次高級(jí)別水系河段同樣為高蜿蜒度,但各子流域河段蜿蜒度最高值均未出現(xiàn)在最高和次高級(jí)別水系中,而是出現(xiàn)在最低和次低級(jí)別水系中,蜿蜒度最低值同樣出現(xiàn)在最低和次低級(jí)別水系中,表明特定流域內(nèi)低級(jí)別水系形態(tài)特征相對(duì)復(fù)雜。

    3.3流域侵蝕產(chǎn)沙與水系形態(tài)特征的關(guān)系

    侵蝕產(chǎn)沙和泥沙輸移受流域地質(zhì)地貌、土壤、植被等自然因素和水利工程建設(shè)、水土保持、河道采砂、開發(fā)項(xiàng)目建設(shè)等人為因素的共同影響。其中,自然因素相對(duì)較穩(wěn)定,對(duì)流域不同時(shí)期侵蝕產(chǎn)沙和泥沙輸移量變化的影響較??;人為因素通常具有時(shí)空異質(zhì)性,是影響流域侵蝕產(chǎn)沙和輸沙量變化的重要因素。1989年起,嘉陵江流域被列為長江上游水土保持重點(diǎn)防治區(qū),實(shí)施的“長治”工程等系列水土保持活動(dòng)改變了侵蝕產(chǎn)沙格局。因此,為了從一定程度上剝離人為因素對(duì)侵蝕產(chǎn)沙變化的影響,分析水系形態(tài)特征對(duì)流域侵蝕產(chǎn)沙的貢獻(xiàn)時(shí),采用的輸沙量是各子流域1960—1990年平均值[18-19],并根據(jù)泥沙輸移比[20]推算出各子流域平均泥沙淤積量(表2)。將數(shù)據(jù)標(biāo)準(zhǔn)化處理后,根據(jù)式(6)—(7)分別計(jì)算出平均輸沙量和淤積量與水系形態(tài)參數(shù)間的關(guān)聯(lián)度。

    圖2 嘉陵江子流域累計(jì)河段數(shù)-水系級(jí)別關(guān)系

    表1 嘉陵江各子流域水系河段蜿蜒特征

    表2 嘉陵江子流域1960-1990年平均輸沙量、泥沙淤積量和水系形態(tài)參數(shù)

    平均輸沙量和泥沙淤積量與分形維數(shù)間的關(guān)聯(lián)度分別為0.618,0.641,表明流域水系復(fù)雜程度對(duì)泥沙淤積的影響強(qiáng)于輸沙量,這是因?yàn)榱饔蛩捣种Y(jié)構(gòu)越復(fù)雜,源區(qū)泥沙向流域出水口運(yùn)動(dòng)過程中,沿程淤積的機(jī)會(huì)越大。河段平均蜿蜒度對(duì)泥沙淤積的影響也強(qiáng)于輸沙量,關(guān)聯(lián)度分別為0.588,0.591,河道越蜿蜒曲折,越不利于泥沙遷移,原因可能是河段蜿蜒度由水系高度落差決定[27],河道比降大時(shí),其蜿蜒度低,有利于導(dǎo)水輸沙,當(dāng)河道經(jīng)過地形高度差較小的區(qū)域時(shí),河道開始變得蜿蜒曲折,有利于泥沙淤積。上述結(jié)果還可以認(rèn)為,流域水系形態(tài)特征對(duì)泥沙淤積的影響強(qiáng)于輸沙量。

    研究表明,流域水系形態(tài)與侵蝕產(chǎn)沙間存在顯著的非線性關(guān)系[28],嘉陵江各子流域輸沙量和泥沙淤積量對(duì)水系形態(tài)特征的非線性響應(yīng)方程分別為(11),(12),均通過了擬合檢驗(yàn):

    Syi=W0.164D0.601R2=0.89

    (11)

    Sde=W0.138D0.649R2=0.89

    (12)

    從式(11)可看出,只考慮河段蜿蜒度和水系分枝結(jié)構(gòu)復(fù)雜度對(duì)流域輸沙量的貢獻(xiàn)時(shí),非線性響應(yīng)模型變量系數(shù)分別是0.164,0.601,可見,河段蜿蜒度對(duì)輸沙量的貢獻(xiàn)較小,相對(duì)貢獻(xiàn)率為21%,分枝結(jié)構(gòu)復(fù)雜度的貢獻(xiàn)較大,其相對(duì)貢獻(xiàn)率為79%。泥沙淤積量對(duì)水系形態(tài)特征的響應(yīng)也表現(xiàn)出相似的規(guī)律(式12),其非線性響應(yīng)模型中河段蜿蜒度系數(shù)為0.138,對(duì)泥沙淤積量的相對(duì)貢獻(xiàn)率為18%,水系分枝結(jié)構(gòu)復(fù)雜程度系數(shù)為0.649,對(duì)泥沙淤積量的相對(duì)貢獻(xiàn)率為82%。

    4 討 論

    流域水系形態(tài)特征一直是水文過程和地貌演化研究的熱點(diǎn),研究方法主要有以霍頓定律為代表的經(jīng)驗(yàn)關(guān)系法和以分形理論為基礎(chǔ)的定量計(jì)算法[7-16],定量計(jì)算水系分形維數(shù)的方法多種多樣,不同的計(jì)算方法所得結(jié)果意義也不一致。本文采用一種改進(jìn)的水系分枝維數(shù)描述水系形態(tài),與分叉比不同,該分形維數(shù)由河流分枝數(shù)目變化率所確定,能表征水系分枝發(fā)育能力和流域地貌溯源侵蝕能力強(qiáng)弱以及復(fù)雜的空間變化特征。利用“時(shí)空對(duì)換”思想,將水系形態(tài)分形維數(shù)等“狀態(tài)數(shù)據(jù)”轉(zhuǎn)化為“過程數(shù)據(jù)”以考查流域侵蝕產(chǎn)沙對(duì)水系形態(tài)的響應(yīng)。結(jié)果表明,流域水系形態(tài)特征對(duì)泥沙淤積的影響強(qiáng)于輸沙量。Jin等[28]認(rèn)為,流域水系形態(tài)與侵蝕產(chǎn)沙間存在顯著的非線性關(guān)系,在流域水系發(fā)育之初期,水系分形維數(shù)會(huì)隨侵蝕產(chǎn)沙模數(shù)減小而減小,當(dāng)水系分形維數(shù)跨過臨界最小值后,則隨侵蝕產(chǎn)沙模數(shù)的減小而增大。這主要是在試驗(yàn)條件下,侵蝕產(chǎn)沙和水系發(fā)育都是兩個(gè)連續(xù)的過程,時(shí)間尺度錯(cuò)動(dòng)較小。而自然狀態(tài)下流域水系發(fā)育和一次連續(xù)的侵蝕產(chǎn)沙過程不在同一個(gè)時(shí)間尺度內(nèi)。因此,反過來可以說,流域水系形態(tài)只是侵蝕產(chǎn)沙的一般影響因素,不是決定性或控制性因素,與Zhang等[29]的研究結(jié)果一致。進(jìn)一步分析流域水系形態(tài)特征對(duì)輸沙量和泥沙淤積量的貢獻(xiàn)時(shí)發(fā)現(xiàn),水系分枝結(jié)構(gòu)復(fù)雜度對(duì)輸沙量和泥沙淤積量的貢獻(xiàn)卻大于河段蜿蜒度,原因可能是輸沙量主要受泥沙來源多寡的影響,復(fù)雜水系形態(tài)演化的過程伴隨著強(qiáng)烈的流域侵蝕;而水系河道作為輸沙通道,其沖淤過程對(duì)輸沙量和泥沙淤積量的影響相對(duì)較小,因此,河段蜿蜒度對(duì)輸沙量和泥沙淤積量的貢獻(xiàn)較小。

    需要說明的是,流域侵蝕產(chǎn)沙是一個(gè)十分復(fù)雜的過程,影響因素眾多。本文只研究了流域水系形態(tài)特征與河段蜿蜒度對(duì)流域輸沙量和泥沙淤積量的貢獻(xiàn),沒有考慮流域結(jié)構(gòu)特征等其他因素的影響,也沒有完全剝離人為因素的干擾。在后續(xù)研究中,應(yīng)注重流域水系形態(tài)特征對(duì)輸沙量和泥沙淤積量的影響機(jī)制研究,以及各類輸沙量、泥沙淤積量影響因子的貢獻(xiàn)率定量研究。

    5 結(jié) 論

    (1)嘉陵江流域內(nèi)各子流域水系形態(tài)存在差異,渠江子流域水系分叉總數(shù)增加最快,分枝形態(tài)最復(fù)雜,溯源侵蝕能力最強(qiáng),其余4個(gè)子流域間差異較小。

    (2)各子流域河段均以低蜿蜒度為主,占河段總數(shù)的比例分別為白龍江子流域78.06%、干流上游子流域82.40%、涪江子流域73.40%、干流中下游子流域73.18%和渠江子流域75.03%。

    (3)平均輸沙量和泥沙淤積量與分形維數(shù)間的關(guān)聯(lián)度分別為0.618,0.641,與河段平均蜿蜒度的關(guān)聯(lián)度分別為0.588,0.591。

    (4)河段蜿蜒度和分枝結(jié)構(gòu)復(fù)雜度對(duì)輸沙量的相對(duì)貢獻(xiàn)率分別為21%和79%,對(duì)泥沙淤積量的相對(duì)貢獻(xiàn)率分別為18%和82%。

    [1]Veltri M,Veltri P,Maiolo M.On the fractal description of natural channel networks[J].Journal of Hydrology,1996,187:137-144.

    [2]沈曉華,鄒樂君,陽峰,等.長江河道分形與流域構(gòu)造特征的關(guān)系[J].浙江大學(xué)學(xué)報(bào):理學(xué)版,2001,28(1):107-111.

    [3]Shen X H,Zou L J,Zhang G F,et al.Fractal characteristics of the main channel of Yellow River and its relation to regional tectonic evolution[J].Geomorphology,2011,127:64-70.

    [4]Liu C L,Duan D Z,Zhang H.Relationships Between Fractal Road and Drainage Networks in Wuling Mountainous Area:Another Symmetric Understanding of Human-Environment Relations[J].Journal of Mountain Science,2014,11(4):1060-1069.

    [5]張麗,戴斌祥,王光謙,等.基于Tokunaga網(wǎng)絡(luò)的河網(wǎng)形態(tài)特征量化[J].中國科學(xué)D輯:地球科學(xué),2009,39(10):1413-1420.

    [6]袁雯,楊凱,吳建平.城市化進(jìn)程中平原河網(wǎng)地區(qū)河流結(jié)構(gòu)特征及其分類方法探討[J].地理科學(xué),2007,27(3):401-407.

    [7]王秀春,吳姍,畢曉麗,等.涇河流域水系分維特征及其生態(tài)意義[J].北京師范大學(xué)學(xué)報(bào):自然科學(xué)版,2004,40(3):364-368.

    [8]Mandelbrot B B.How long is the coast of Britain? Statistical self-similarity and fractional dimension[J].Science,1967,155:636-638.

    [9]馮平,馮焱.河流形態(tài)特征的分維計(jì)算方法[J].地理學(xué)報(bào),1997,52(4):324-330.

    [10]Tartobon D G,Bras R L,Rodrguez-iturbe I.The fractal nature of river networks[J].Water Resource Research,1988,24(8):1317-1322.

    [11]Robert A,Roy A G.On the fractal interpretation of the main stream,length drainage area relationship[J].Water Resource Research,1990,26(5):839-842.

    [12]Rosso R,Bacchi B,Barbera P L.Fractal relation of mainstream length to catchment area in river networks[J].Water Resource Research,1991,27(3):381-387.

    [13]孫祝友,杜國云,朱大奎,等.萊州灣東岸河流的分形特征與流域地貌發(fā)育研究[J].地理科學(xué),2010,30(5):755-759.

    [14]蔡凌雁,湯國安,熊禮陽,等.基于DEM的陜北黃土高原典型地貌分形特征研究[J].水土保持通報(bào),2014,34(3):141-144.

    [15]Garcia-Ruiz J M,Otálora F.Fractal trees and Horton's laws[J].Mathematical Geology,1992,24(1):61-71.

    [16]Kim J C,Jung K.Fractal tree analysis of drainage patterns[J].Water Resources Management,2015,29(4):1217-1230.

    [17]柏春廣,蔡先華.南京市交通網(wǎng)絡(luò)的分形特征[J].地理研究,2008,27(6):1419-1426.

    [18]馬煉,張明波,郭海晉,等.嘉陵江流域水保治理前后沿程水沙變化研究[J].水文,2002,22(1):27-31.

    [19]許全喜,陳松生,熊明,等.嘉陵江流域水沙變化特性及原因分析[J].泥沙研究,2008(2):1-8.

    [20]李繼承.嘉陵江流域非點(diǎn)源污染負(fù)荷模擬研究[D].重慶:重慶大學(xué),2007.

    [21]劉臘美,龍?zhí)煊?李崇明.三峽水庫上游流域非點(diǎn)源顆粒態(tài)磷污染負(fù)荷研究[J].長江流域資源與環(huán)境,2009,18(4):320-325.

    [22]韋杰,賀秀斌.人類活動(dòng)對(duì)嘉陵江流域泥沙負(fù)荷的影響[J].長江流域資源與環(huán)境,2010,19(2):196-201.

    [23]張斌,艾南山,黃正文,等.中國嘉陵江河曲的形態(tài)與成因[J].科學(xué)通報(bào),2007,52(22):2671-2682.

    [24]黃汲清,任紀(jì)舜,姜春發(fā),等.中國大地構(gòu)造基本輪廓[J].地質(zhì)學(xué)報(bào),1977(2):117-135.

    [25]潘桂棠,肖慶輝,陸松年,等.中國大地構(gòu)造單元?jiǎng)澐諿J].中國地質(zhì),2009,36(1):1-4.

    [26]賈營營,付碧宏,王巖,等.青藏高原東緣龍門山斷裂帶晚新生代構(gòu)造地貌生長及水系響應(yīng)[J].第四紀(jì)研究,2010,30(4):825-836.

    [27]Rosgen D L.A classification of natural rivers[J].Catena,1994,22(3):169-199.

    [28]Jin D S,Chen H,Guo Q W.Material component to non-linear relation between sediment yield and drainage network development: A flume experimental study[J].Journal of Geographical Sciences,2001,24(3):271-381.

    [29]Zhang H Y,Shi Z H,Fang N F,et al.Linking watershed geomorphic characteristics to sediment yield: Evidence from the Loess Plateau of China[J].Geomorphology,2015,234:19-27.

    Features of River System Morphology and Its Response to Sediment Yield in the Jialing River Basin,China

    LI Jinlin1,WEI Jie1,2

    (1.College of Geography and Tourism,Chongqing Normal University,Chongqing 401331,China; 2.Chongqing Key Laboratory of Surface Process and Environment Remote Sensing in the Three Gorges Reservoir Area,Chongqing 401331,China)

    The features of river system morphology represent the geomorphic development and hydrological processes,which helps make planning for hydropower development,ecological protection and soil and water conservation.However,it is not well understood in the Jialing River Basin,as well as its relation to the sediment yield and deposition.The fractal dimension and the sinuosity of river reaches are used to quantify the river system morphology in the Jialing River Basin.And then,the grey relation degree and nonlinear response model are applied to analyze the relationships between the fractal dimension,the sinuosity of river reaches and the sediment load and deposition,based on the DEM data by GIS software.The results show that the significant differences occur in the river system morphology among sub-basins.Among these sub-basins,the fractal dimension of the Qu River is 5.102,with a fastest bifurcation increment,the most complex branches and the most severe headwater erosion potential.The fractal dimensions for the Bailong River,the upper reach of the main stream,the middle and lower reaches of main stream and the Fu River are 3.706,3.605,3.576 and 3.548,respectively,indicating that no significant differences were found among the rest 4 sub-basins except for the Qu River from a perspective of branches complexity.Most of the reaches in each sub-basin are gentle meandering with a proportion about 73%~82% of the total,and the moderate meandering about 12%~20%.The number of high sinuosity reaches in the upper reaches of main stream,the middle and lower reaches of main stream,the Qu Rivers,the Fu River and the Bailong River is 79,137,158,168 and 247,respectively.The grey relation degrees between the mean sediment loads,the depositions and the fractal dimensions are 0.618 and 0.641,respectively,showing stronger influence of the river system complexity on the sediment deposition than the sediment load.And the grey relation degrees between the mean sediment loads,the depositions and the sinuosity of rive reaches are 0.588 and 0.591,respectively,indicating that a more sensitive response of the sediment deposition to the sinuosity of rive reaches than the sediment load.In addition,the relative contributions of the sinuosity of river reaches and the branches structure complexity to the sediment load is 21% and 79%,while they are 18% and 82% for the sediment deposition.

    fractal dimension; sinuosity of river reach; grey correlation degree; relative contribution; basin erosion and riverine sediment

    2015-12-29

    2016-02-23

    重慶市教委科技項(xiàng)目“三峽庫區(qū)河流泥沙負(fù)荷變化的驅(qū)動(dòng)力判識(shí)研究”(KJ100601)

    李進(jìn)林(1991—),男,重慶開州人,碩士研究生,主要研究方向?yàn)樗帘3峙c生態(tài)建設(shè)。E-mail:cqnulijinlin@sina.com

    韋杰(1979—),男,四川蒼溪人,博士,教授,碩士生導(dǎo)師,主要研究方向?yàn)橥寥狼治g與水土保持。E-mail:wei_jie@mails.ucas.ac.cn

    K903;P931

    A

    1005-3409(2016)05-0007-07

    猜你喜歡
    輸沙量水系河段
    鄱陽湖水系之潦河
    SL流量計(jì)在特殊河段的應(yīng)用——以河源水文站為例
    環(huán)水系旅游方案打造探析——以臨沂市開發(fā)區(qū)水系為例
    石瀘高速公路(紅河段)正式通車
    石油瀝青(2020年1期)2020-05-25 06:54:10
    水系魔法之止水術(shù)
    遼河干流遼中站河道沖淤變化趨勢(shì)研究
    氣候變化和人類活動(dòng)對(duì)祖厲河輸沙量變化的影響分析
    頻率曲線在荊江三口輸沙量計(jì)算中的應(yīng)用
    關(guān)于遼河干流河道沖淤量沿程變化規(guī)律探討
    地下水(2015年5期)2015-12-02 02:44:26
    環(huán)境友好的CLEAN THROUGH水系洗滌劑
    春色校园在线视频观看| 伦理电影大哥的女人| 国产亚洲精品av在线| 你懂的网址亚洲精品在线观看 | 婷婷六月久久综合丁香| 色哟哟·www| 91精品伊人久久大香线蕉| 亚洲人与动物交配视频| 久久久久久久久中文| 精品熟女少妇av免费看| 国产精品一区二区三区四区久久| 丰满少妇做爰视频| 非洲黑人性xxxx精品又粗又长| 丝袜喷水一区| 国产成人精品久久久久久| 午夜a级毛片| 久久久久久大精品| 国产精品熟女久久久久浪| videos熟女内射| 国产精品.久久久| 久久99精品国语久久久| 免费观看人在逋| 久久人人爽人人片av| 成人午夜精彩视频在线观看| 国产高潮美女av| 美女黄网站色视频| 久久久久九九精品影院| 级片在线观看| 国产午夜福利久久久久久| 七月丁香在线播放| 国内精品宾馆在线| 免费一级毛片在线播放高清视频| 欧美精品一区二区大全| 久久久精品94久久精品| 秋霞在线观看毛片| 内射极品少妇av片p| 亚洲欧美精品专区久久| 免费大片18禁| 日本黄色片子视频| 97超碰精品成人国产| 成人欧美大片| 免费观看性生交大片5| 两个人视频免费观看高清| 六月丁香七月| 成人无遮挡网站| 毛片女人毛片| 好男人在线观看高清免费视频| 国产69精品久久久久777片| av黄色大香蕉| 在线播放无遮挡| 色网站视频免费| 纵有疾风起免费观看全集完整版 | 亚洲人成网站在线观看播放| 少妇高潮的动态图| 亚洲中文字幕一区二区三区有码在线看| 人人妻人人澡欧美一区二区| 国产成人免费观看mmmm| 我要看日韩黄色一级片| 嫩草影院新地址| 男人和女人高潮做爰伦理| 啦啦啦啦在线视频资源| kizo精华| 国产不卡一卡二| 日日摸夜夜添夜夜爱| 老司机福利观看| 国产又黄又爽又无遮挡在线| 中文字幕免费在线视频6| 国产精品国产三级国产av玫瑰| 国产成人精品婷婷| 高清av免费在线| 国产黄a三级三级三级人| 听说在线观看完整版免费高清| 女人十人毛片免费观看3o分钟| 搡老妇女老女人老熟妇| 中文字幕免费在线视频6| 综合色丁香网| 中文天堂在线官网| 一级毛片我不卡| 好男人视频免费观看在线| 在线观看av片永久免费下载| 伦理电影大哥的女人| 国产在视频线在精品| 国产真实乱freesex| 日日撸夜夜添| 久久国产乱子免费精品| 波野结衣二区三区在线| 一级av片app| 国产精品久久久久久精品电影| 中文亚洲av片在线观看爽| 国产精品电影一区二区三区| 18禁动态无遮挡网站| 波多野结衣高清无吗| 黑人高潮一二区| 午夜福利成人在线免费观看| 免费电影在线观看免费观看| 免费播放大片免费观看视频在线观看 | 免费播放大片免费观看视频在线观看 | 日本三级黄在线观看| 国产精品一区www在线观看| 久久精品久久久久久久性| 大香蕉97超碰在线| 亚洲熟妇中文字幕五十中出| 国产高清不卡午夜福利| 能在线免费观看的黄片| 久久6这里有精品| 午夜福利在线观看免费完整高清在| 亚洲精品456在线播放app| 国产免费视频播放在线视频 | 国产在线一区二区三区精 | 国内揄拍国产精品人妻在线| av播播在线观看一区| 精品熟女少妇av免费看| 国产探花极品一区二区| 日韩国内少妇激情av| 精品国内亚洲2022精品成人| 亚洲一区高清亚洲精品| 最近的中文字幕免费完整| 久久精品综合一区二区三区| 国产精品乱码一区二三区的特点| 真实男女啪啪啪动态图| 蜜桃亚洲精品一区二区三区| 又爽又黄无遮挡网站| 国产乱来视频区| 久久久久网色| 乱系列少妇在线播放| 欧美成人一区二区免费高清观看| 少妇裸体淫交视频免费看高清| 国国产精品蜜臀av免费| АⅤ资源中文在线天堂| 国产亚洲一区二区精品| 看十八女毛片水多多多| 又粗又硬又长又爽又黄的视频| 国产淫片久久久久久久久| 联通29元200g的流量卡| 久久国产乱子免费精品| 禁无遮挡网站| 亚洲成av人片在线播放无| 日本五十路高清| 欧美性猛交╳xxx乱大交人| 欧美丝袜亚洲另类| 在线播放无遮挡| 只有这里有精品99| 99久久无色码亚洲精品果冻| 国产成人a区在线观看| 两个人的视频大全免费| 少妇裸体淫交视频免费看高清| 中文字幕免费在线视频6| 色网站视频免费| 免费电影在线观看免费观看| 麻豆乱淫一区二区| 白带黄色成豆腐渣| 国产黄a三级三级三级人| 精品无人区乱码1区二区| 亚洲性久久影院| 日本与韩国留学比较| 搡老妇女老女人老熟妇| av在线老鸭窝| 亚洲精品自拍成人| 91久久精品国产一区二区成人| 别揉我奶头 嗯啊视频| 日本-黄色视频高清免费观看| 一边亲一边摸免费视频| 免费看日本二区| 网址你懂的国产日韩在线| 只有这里有精品99| 免费播放大片免费观看视频在线观看 | 97人妻精品一区二区三区麻豆| 九九在线视频观看精品| 欧美+日韩+精品| 精品少妇黑人巨大在线播放 | 淫秽高清视频在线观看| 久久精品久久精品一区二区三区| 在现免费观看毛片| 国产成人aa在线观看| 久久久久久久久中文| 91av网一区二区| 六月丁香七月| 少妇的逼水好多| 亚洲精品亚洲一区二区| 日韩成人av中文字幕在线观看| 中文字幕制服av| 精品国产一区二区三区久久久樱花 | 男人的好看免费观看在线视频| 可以在线观看毛片的网站| 黄片无遮挡物在线观看| 一本久久精品| 内射极品少妇av片p| 亚洲欧洲国产日韩| 水蜜桃什么品种好| 亚洲国产精品合色在线| 床上黄色一级片| 成人鲁丝片一二三区免费| 国产黄片美女视频| 淫秽高清视频在线观看| 内射极品少妇av片p| 国产一区二区在线观看日韩| 在线免费十八禁| av在线亚洲专区| 亚洲美女视频黄频| 午夜激情欧美在线| 在线观看66精品国产| 国产男人的电影天堂91| 国产精品伦人一区二区| 亚洲av日韩在线播放| 国产精品女同一区二区软件| av在线老鸭窝| 日本午夜av视频| 国语对白做爰xxxⅹ性视频网站| 精品国内亚洲2022精品成人| 亚洲av成人精品一二三区| 麻豆成人午夜福利视频| 欧美激情久久久久久爽电影| 男女边吃奶边做爰视频| 国产在线一区二区三区精 | 亚洲精品国产av成人精品| 久久久a久久爽久久v久久| 亚洲激情五月婷婷啪啪| 美女内射精品一级片tv| 久久久久免费精品人妻一区二区| 高清日韩中文字幕在线| 婷婷色av中文字幕| 一区二区三区高清视频在线| 一区二区三区四区激情视频| 久久久欧美国产精品| 国内揄拍国产精品人妻在线| a级毛色黄片| av天堂中文字幕网| 中国国产av一级| 色尼玛亚洲综合影院| 日韩人妻高清精品专区| 一本久久精品| 国产精品.久久久| 国产亚洲午夜精品一区二区久久 | 麻豆精品久久久久久蜜桃| 波野结衣二区三区在线| 一个人免费在线观看电影| 国产大屁股一区二区在线视频| 亚洲精品亚洲一区二区| 国内精品宾馆在线| 看十八女毛片水多多多| 在线观看美女被高潮喷水网站| videos熟女内射| 午夜日本视频在线| 毛片女人毛片| 99久久九九国产精品国产免费| 国产成人福利小说| 国产精品一二三区在线看| 亚洲在线自拍视频| 黄色日韩在线| 久久久成人免费电影| 性插视频无遮挡在线免费观看| 成人一区二区视频在线观看| 深夜a级毛片| 美女xxoo啪啪120秒动态图| 国产精品一区www在线观看| 少妇的逼水好多| 两个人的视频大全免费| 国产在线男女| 国产极品精品免费视频能看的| 亚洲精品色激情综合| 一级毛片我不卡| 精品人妻一区二区三区麻豆| 国产黄片美女视频| 国产成人91sexporn| 美女cb高潮喷水在线观看| 最近最新中文字幕大全电影3| 全区人妻精品视频| 国产一区二区亚洲精品在线观看| 久久国内精品自在自线图片| 老司机影院毛片| 日本免费a在线| 国产日韩欧美在线精品| 午夜a级毛片| 热99re8久久精品国产| 日本黄色视频三级网站网址| 久久韩国三级中文字幕| 国产成人freesex在线| 男女国产视频网站| 免费人成在线观看视频色| 搡女人真爽免费视频火全软件| 狂野欧美白嫩少妇大欣赏| 在线观看66精品国产| 国产视频首页在线观看| 精品国内亚洲2022精品成人| 午夜a级毛片| 美女高潮的动态| 1000部很黄的大片| 99久久精品国产国产毛片| 国产精品.久久久| 超碰av人人做人人爽久久| 秋霞伦理黄片| 中文精品一卡2卡3卡4更新| 久久精品国产鲁丝片午夜精品| 少妇人妻一区二区三区视频| 午夜激情福利司机影院| 边亲边吃奶的免费视频| 国模一区二区三区四区视频| 18+在线观看网站| 国产精品一区二区三区四区久久| 午夜久久久久精精品| 免费播放大片免费观看视频在线观看 | 秋霞在线观看毛片| 久久精品人妻少妇| 国产黄色视频一区二区在线观看 | 精品无人区乱码1区二区| 人妻夜夜爽99麻豆av| 性色avwww在线观看| www.av在线官网国产| 狂野欧美激情性xxxx在线观看| 国产一区有黄有色的免费视频 | 国产精品人妻久久久影院| 国产精品久久电影中文字幕| 看非洲黑人一级黄片| 色视频www国产| ponron亚洲| 99久国产av精品国产电影| 免费av不卡在线播放| 嫩草影院新地址| 亚洲成av人片在线播放无| 国产一区二区亚洲精品在线观看| 成人欧美大片| 国国产精品蜜臀av免费| 波多野结衣巨乳人妻| 少妇的逼好多水| 久久久久久久亚洲中文字幕| 亚洲精品色激情综合| 午夜久久久久精精品| 九九久久精品国产亚洲av麻豆| 欧美成人精品欧美一级黄| 亚洲真实伦在线观看| 人妻少妇偷人精品九色| 三级国产精品欧美在线观看| av女优亚洲男人天堂| 你懂的网址亚洲精品在线观看 | 国产乱人视频| 热99re8久久精品国产| 欧美一区二区精品小视频在线| 午夜精品在线福利| av在线蜜桃| 男人的好看免费观看在线视频| 亚洲成人av在线免费| av视频在线观看入口| 国产成人freesex在线| 青青草视频在线视频观看| 中文资源天堂在线| 成年女人永久免费观看视频| 亚洲人成网站在线观看播放| 成人特级av手机在线观看| 一级毛片aaaaaa免费看小| 国产一区亚洲一区在线观看| 大又大粗又爽又黄少妇毛片口| 五月玫瑰六月丁香| 晚上一个人看的免费电影| 村上凉子中文字幕在线| 成人鲁丝片一二三区免费| 2021天堂中文幕一二区在线观| 免费在线观看成人毛片| 国产极品天堂在线| 国产中年淑女户外野战色| 国产爱豆传媒在线观看| 少妇人妻精品综合一区二区| 欧美最新免费一区二区三区| 少妇猛男粗大的猛烈进出视频 | 91狼人影院| 日韩强制内射视频| 欧美人与善性xxx| 亚洲av电影在线观看一区二区三区 | 麻豆成人午夜福利视频| 纵有疾风起免费观看全集完整版 | 国产精品伦人一区二区| 在线观看一区二区三区| 中文字幕制服av| 男女下面进入的视频免费午夜| 中文字幕免费在线视频6| 成人国产麻豆网| 中文字幕免费在线视频6| 男女下面进入的视频免费午夜| 精品99又大又爽又粗少妇毛片| www日本黄色视频网| 少妇熟女欧美另类| 亚洲精品日韩av片在线观看| 午夜福利成人在线免费观看| 亚洲国产精品专区欧美| 成年av动漫网址| 婷婷六月久久综合丁香| 久久精品国产99精品国产亚洲性色| 国产精品福利在线免费观看| 美女高潮的动态| 国产 一区精品| 国产精品伦人一区二区| 国产伦精品一区二区三区视频9| 国产精品久久电影中文字幕| 免费电影在线观看免费观看| 联通29元200g的流量卡| 成年女人永久免费观看视频| 尾随美女入室| av在线播放精品| 免费无遮挡裸体视频| av卡一久久| 大香蕉久久网| 亚洲精品国产av成人精品| 熟女电影av网| 人妻少妇偷人精品九色| 日本与韩国留学比较| 亚洲精品456在线播放app| 久久精品国产自在天天线| 国产久久久一区二区三区| 久久99热这里只有精品18| 丝袜喷水一区| 有码 亚洲区| 欧美高清成人免费视频www| 偷拍熟女少妇极品色| 久久久久久久久久成人| 亚洲五月天丁香| 欧美区成人在线视频| 国产精品熟女久久久久浪| 国产女主播在线喷水免费视频网站 | 韩国高清视频一区二区三区| av国产免费在线观看| 别揉我奶头 嗯啊视频| 国产精品国产三级国产av玫瑰| 国产色婷婷99| 九九热线精品视视频播放| 精品久久久久久久人妻蜜臀av| 99九九线精品视频在线观看视频| 黄色一级大片看看| 能在线免费看毛片的网站| 寂寞人妻少妇视频99o| 亚洲人与动物交配视频| 亚洲精品日韩av片在线观看| av黄色大香蕉| 国产亚洲午夜精品一区二区久久 | a级一级毛片免费在线观看| 国产精品精品国产色婷婷| 亚洲在久久综合| av卡一久久| 免费观看在线日韩| 三级毛片av免费| 只有这里有精品99| 九九爱精品视频在线观看| 国产精品国产三级国产av玫瑰| 久久精品熟女亚洲av麻豆精品 | 国模一区二区三区四区视频| 日韩欧美 国产精品| 中文字幕精品亚洲无线码一区| 深夜a级毛片| 久久精品国产亚洲av天美| 99久久中文字幕三级久久日本| 免费观看精品视频网站| 国产片特级美女逼逼视频| 老师上课跳d突然被开到最大视频| 国产精品国产三级国产av玫瑰| 国产欧美日韩精品一区二区| 欧美一区二区亚洲| 成人亚洲欧美一区二区av| 99在线视频只有这里精品首页| 国产精品国产三级专区第一集| 国产精品一二三区在线看| 男女下面进入的视频免费午夜| 久久久a久久爽久久v久久| 久久精品综合一区二区三区| 人妻制服诱惑在线中文字幕| 能在线免费看毛片的网站| 一区二区三区高清视频在线| 国产精品伦人一区二区| 免费电影在线观看免费观看| 国产又色又爽无遮挡免| av视频在线观看入口| 全区人妻精品视频| 少妇熟女欧美另类| 婷婷六月久久综合丁香| av专区在线播放| 成人特级av手机在线观看| 日韩,欧美,国产一区二区三区 | 日韩av在线免费看完整版不卡| 日本黄色片子视频| 亚洲怡红院男人天堂| 亚洲av成人精品一区久久| or卡值多少钱| 免费av不卡在线播放| 亚洲国产精品成人久久小说| 91av网一区二区| 国产综合懂色| 床上黄色一级片| 日韩国内少妇激情av| 亚洲va在线va天堂va国产| 中文天堂在线官网| 老司机福利观看| 亚洲自偷自拍三级| 欧美人与善性xxx| 午夜老司机福利剧场| 国内少妇人妻偷人精品xxx网站| 色视频www国产| 美女国产视频在线观看| 人人妻人人看人人澡| 日日摸夜夜添夜夜爱| 国产精品久久久久久久电影| 麻豆乱淫一区二区| 国产精品av视频在线免费观看| 老女人水多毛片| 97在线视频观看| 99久久九九国产精品国产免费| 一本一本综合久久| 2021天堂中文幕一二区在线观| 国产乱人偷精品视频| 精品久久久久久久久av| 成人高潮视频无遮挡免费网站| 国产在视频线精品| 日本与韩国留学比较| 中文字幕人妻熟人妻熟丝袜美| 精品久久国产蜜桃| 欧美不卡视频在线免费观看| 久久久久网色| 黄片wwwwww| 欧美一级a爱片免费观看看| 国产精品精品国产色婷婷| 欧美日韩综合久久久久久| 中文字幕精品亚洲无线码一区| 成人国产麻豆网| 秋霞伦理黄片| 国产精品久久久久久av不卡| 日韩欧美 国产精品| 18禁动态无遮挡网站| 亚洲av成人av| 欧美日韩国产亚洲二区| 日韩高清综合在线| 国产亚洲5aaaaa淫片| 日韩欧美国产在线观看| 日本熟妇午夜| 美女国产视频在线观看| 国产精品蜜桃在线观看| 久久国产乱子免费精品| 狂野欧美白嫩少妇大欣赏| 久久99精品国语久久久| 2021少妇久久久久久久久久久| 亚洲欧美精品综合久久99| 日本-黄色视频高清免费观看| 日韩av在线大香蕉| 少妇的逼水好多| 少妇的逼好多水| 久久久久久久久中文| 国产精品精品国产色婷婷| 精品少妇黑人巨大在线播放 | 插阴视频在线观看视频| 高清午夜精品一区二区三区| 男插女下体视频免费在线播放| 亚洲国产精品合色在线| 午夜福利在线观看吧| av在线蜜桃| 能在线免费观看的黄片| 日产精品乱码卡一卡2卡三| 婷婷色av中文字幕| 日韩中字成人| 亚洲国产精品久久男人天堂| av免费在线看不卡| 丝袜喷水一区| 麻豆国产97在线/欧美| 日韩在线高清观看一区二区三区| 日韩成人av中文字幕在线观看| 国产白丝娇喘喷水9色精品| 男人的好看免费观看在线视频| 麻豆久久精品国产亚洲av| 非洲黑人性xxxx精品又粗又长| 汤姆久久久久久久影院中文字幕 | 春色校园在线视频观看| 三级国产精品欧美在线观看| 韩国av在线不卡| 91aial.com中文字幕在线观看| 午夜精品国产一区二区电影 | 亚洲天堂国产精品一区在线| 精品国产三级普通话版| 国产高清三级在线| 有码 亚洲区| 伊人久久精品亚洲午夜| 18禁裸乳无遮挡免费网站照片| 中文天堂在线官网| 晚上一个人看的免费电影| 爱豆传媒免费全集在线观看| 欧美97在线视频| 美女xxoo啪啪120秒动态图| 久久亚洲国产成人精品v| 午夜亚洲福利在线播放| 亚洲成人av在线免费| 日韩亚洲欧美综合| 精品人妻视频免费看| 国产乱人偷精品视频| 九色成人免费人妻av| 欧美人与善性xxx| 91精品伊人久久大香线蕉| 欧美三级亚洲精品| 建设人人有责人人尽责人人享有的 | 日韩一区二区视频免费看| 欧美日韩在线观看h| 熟女人妻精品中文字幕| av播播在线观看一区| 亚洲久久久久久中文字幕| 亚洲国产欧美在线一区| 18+在线观看网站| 久久韩国三级中文字幕| 夜夜看夜夜爽夜夜摸| 秋霞在线观看毛片| 午夜久久久久精精品| 噜噜噜噜噜久久久久久91| 久久鲁丝午夜福利片| 最后的刺客免费高清国语| 亚洲精品亚洲一区二区| 国产91av在线免费观看| 爱豆传媒免费全集在线观看| 久久久久免费精品人妻一区二区| 天天躁日日操中文字幕| av在线蜜桃| 99久国产av精品| 国产精品,欧美在线| 免费一级毛片在线播放高清视频| 日本av手机在线免费观看|