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

    上游來水來沙與洞庭湖入?yún)R頂托對荊江河段沖淤的影響

    2024-01-01 00:00:00尚海鑫胡春宏夏軍強周美蓉
    水科學進展 2024年3期

    DOI:10.14042/j.cnki.32.1309.2024.03.012

    摘要:三峽工程運用后荊江河段處于持續(xù)沖刷狀態(tài),分析上游來水來沙與洞庭湖入?yún)R頂托對荊江河段沖淤的影響,對掌握壩下游河道演變規(guī)律具有重要意義。本研究基于2003—2020年荊江河段實測水沙及沖淤量資料,建立河段累計沖淤量與河段上下游邊界條件之間的函數(shù)關系,并重點分析頂托作用對荊江河段沖淤的影響。結果表明:① 荊江河段累計沖淤量與河段進口水流沖刷強度呈正相關,與進出口水位落差呈負相關,河段沖淤過程是二者共同作用的結果;② 建立進口水流沖刷強度、進出口水位落差與荊江河段累計沖淤量的函數(shù)關系,能較好反映上下游邊界條件對河段沖淤過程的影響(R2gt;0.92);③ 2003—2020年荊江河段累計沖刷量為12.3億m3,在假設無洞庭湖入?yún)R頂托的情況下比實測值增加0.67億m3,即入?yún)R頂托能使荊江河段沖刷量減少約5%;同期上、下荊江累計沖刷量分別為7.3億和5.0億m3,入?yún)R頂托能使沖刷量分別減少約1%和16%。洞庭湖入?yún)R頂托會導致荊江河段的沖刷量減少,對下荊江沖淤影響更顯著。

    關鍵詞:頂托作用;河床沖淤;水流沖刷強度;洞庭湖入?yún)R;荊江河段

    中圖分類號:TV142

    文獻標志碼:A

    文章編號:1001-6791(2024)03-0485-11

    收稿日期:2024-01-18;網(wǎng)絡出版日期:2024-04-29

    網(wǎng)絡出版地址:https:∥link.cnki.net/urlid/32.1309.P.20240429.0909.002

    基金項目:國家自然科學基金資助項目(U2040215;U2240206)

    作者簡介:尚海鑫(1995—),男,甘肅敦煌人,博士研究生,主要從事河流動力學研究。E-mail:shanghx@whu.edu.cn

    通信作者:夏軍強,E-mail:xiajq@whu.edu.cn

    當天然河流的進口水沙條件發(fā)生改變,或由支流側(cè)向入?yún)R頂托導致河流出口邊界條件發(fā)生改變時,河道沖淤過程會隨之變化[1-2]。三峽工程及上游水庫群運行改變了壩下游河道的水沙過程,沙量劇減,引起了長時間、長距離的河床沖刷[3]。而位于荊江河段出口的洞庭湖入?yún)R對荊江河段水位有明顯的頂托作用[4-5]。洞庭湖入?yún)R頂托會影響荊江河段出口的侵蝕基面,也是影響河段沖淤的因素之一。在荊江河段仍處于持續(xù)沖刷狀態(tài)的背景下,厘清壩下游河道沖刷的影響因素,定量分析上游來水來沙與洞庭湖入?yún)R頂托對荊江河段沖淤的影響,對掌握壩下游河道演變規(guī)律具有重要意義。

    三峽工程運用后壩下游河段持續(xù)沖刷受到了廣泛關注。許全喜等[3]深入分析了三峽工程運用以來長江中下游的水沙變化與河道沖淤特征;盧金友等[5]研究了水庫群聯(lián)合作用下三峽壩下游河道的水沙輸移特征及其影響,認為受水庫群運行的影響,進入長江中下游的水沙條件發(fā)生顯著變化,導致壩下游河道發(fā)生長距離的沖刷;許炯心[6]分析了含沙量和懸沙粒徑變化對長江宜昌—漢口河段年沖淤量的影響,建立了河段沖淤量與含沙量、最大流量、三口分流比等影響因素的經(jīng)驗關系??梢?,已有研究主要分析河段沖淤對上游來水來沙量改變的響應[7-8]。針對荊江河段沖淤的研究,較少同時考慮洞庭湖入?yún)R頂托的影響,故難以厘清三峽工程運用后荊江河段沖淤的影響因素。

    目前已有許多針對洞庭湖入?yún)R對長江干流水沙過程影響的研究,包括三峽工程運用后洞庭湖出流過程的變化特征[9-10]、洞庭湖入?yún)R對荊江水位的頂托程度[11-12]、洞庭湖水沙入?yún)R對干流河段造床流量的影響機制等[13]。已有研究成果揭示了洞庭湖入?yún)R頂托作用對荊江河段水文水動力特性的影響,表明洞庭湖入?yún)R頂托影響了荊江河段的水位、流速等,改變了河段出口侵蝕基面,是影響荊江河段沖淤變化的重要因素。但目前仍缺少關于上游來水來沙與洞庭湖入?yún)R頂托的耦合作用對荊江河段沖淤影響的量化研究。荊江河段在持續(xù)沖刷過程中,枝城站水沙條件可視為河段進口邊界條件,洞庭湖入?yún)R頂托則決定了河段出口的侵蝕基面,二者會共同影響河段的沖淤過程。因此,有必要定量分析水沙條件與頂托作用對荊江河段沖淤過程的影響機制。

    本研究通過分析三峽工程運用后荊江河段的水文泥沙及河段沖淤量資料,建立上、下荊江和荊江河段累計沖淤量與河段上游水沙條件、下游出口頂托作用之間的函數(shù)關系,定量計算頂托作用對荊江河段沖淤的影響。

    1" 研究河段概況

    荊江河段進口(枝城)位于三峽大壩下游約102 km處,出口位于城陵磯,全長約347 km,以藕池口為界,分為上、下荊江。上荊江屬于微彎分汊型河道,長約172 km;下荊江屬于典型的彎曲型河道,長約175 km(圖1)。本研究選擇枝城站的實測水沙資料代表上荊江的來水來沙條件;監(jiān)利站位于藕池口下游,至城陵磯之間無分流和支流入?yún)R影響,故選擇監(jiān)利站的實測水沙資料近似代表下荊江的來水來沙條件。新廠水位站位于藕池口上游約7 km,蓮花塘水位站位于城陵磯,選擇枝城、新廠站的實測水位資料近似代表上荊江進、出口水位條件;新廠、蓮花塘站的實測水位資料近似代表下荊江進、出口水位條件。洞庭湖匯集“荊南三口”分流與湘、資、沅、澧“四水”來流,于城陵磯處匯入長江干流[14],七里山站是洞庭湖出口控制站(匯流口上游約5 km),選擇該站的實測流量資料代表洞庭湖出流的水文情勢變化。

    本研究收集了1999—2020年枝城、監(jiān)利和七里山站日均水位、流量、含沙量資料,新廠、蓮花塘水位站的日均水位資料,以及2003—2020年荊江河段(包括上、下荊江)的年沖淤量資料。三峽工程運用后(2003—2020年)枝城站年均徑流量約為4 283億m3,與蓄水前(1990—2002年)相比偏少4%(圖2)。由于水庫上游來沙偏少和三峽水庫攔沙作用,三峽工程運用后進入荊江河段的沙量顯著減小,枝城站年均輸沙量為0.4億t,較蓄水前年均值減少92%。2003—2020年三口合計年均徑流量約479億m3,年均輸沙量約0.084億t,平均分流比和分沙比分別為12%和18%;由洞庭湖入?yún)R長江干流年徑流量的變化范圍為1 955億~4 885億m3,多年平均徑流量為2 550億m3,約為枝城站徑流量的60%。三峽工程運用后,荊江河段基本處于持續(xù)沖刷狀態(tài)。2003—2020年荊江河段平灘河槽累計沖刷量約12.3億m3,其中上、下荊江累計沖刷量分別為7.3億和5.0億m3,分別占荊江河段沖刷量的60%和40%。

    2" 上游水沙與入?yún)R頂托作用與河段沖淤的關系

    河床演變學原理表明河段進口的水沙條件和河段出口侵蝕基面是影響河床演變的主要因素。此處首先分析進口水流沖刷強度和洞庭湖入?yún)R頂托作用對荊江河段沖淤過程的影響機制,然后建立荊江河段累計沖淤量與上述2個影響因素的函數(shù)關系。

    2.1" 上下游邊界條件對荊江段沖淤的影響機理

    當上游建壩后下泄清水時,壩下游河道將會產(chǎn)生沿程沖刷;當河道下游侵蝕基面抬升時,例如洞庭湖入?yún)R對荊江河段水位的頂托作用(圖3,其中Zu為進口水位,Zd為出口水位,ΔZ為水位落差),河道沖刷量會相應減少。下面將具體分析上述2類主要因素對河道沖淤的影響。

    2.1.1" 河段進口水沙條件

    受三峽工程及其上游水庫群運用的影響,進入荊江河段的沙量大幅減小。對于沖積河流,一般采用水流沖刷強度(Ffi)[15]表示河段來水來沙條件:

    Ffi=1Nf∑Nfi=1(Q2fi/Sfi)×10-8(1)

    式中:Qfi和Sfi分別為日均流量(m3/s)和懸移質(zhì)含沙量(kg/m3);Nf為當年天數(shù)。三峽工程運用前,長江中游典型斷面輸沙率(QS)與流量的平方(Q2)成正比關系。在本研究中Q2fi近似代表河段進口的水流挾沙能力,Q2fi/Sfi則代表了某一特定流量下挾沙力與含沙量的比值,F(xiàn)fi越大表示水流沖刷強度越大。

    以枝城站實測水沙資料代表上荊江與荊江河段的水沙條件,監(jiān)利站實測資料代表下荊江水沙條件。通過式(1)計算上、下荊江與荊江的Ffi,點繪河段平灘河槽累計沖淤量(V)與河段Ffi的關系并進行Pearson相關性檢驗,可發(fā)現(xiàn)3個河段的V與Ffi呈顯著的相關關系(p<0.01)(圖4)。

    荊江南岸有松滋口、太平口和藕池口分流分沙,也會影響荊江河段的水沙條件。此處從枝城站日均流量、含沙量中扣除三口的日均流量、含沙量,并由式(1)計算受三口分流影響的河段水流沖刷強度(F*fi)。采用α表示分流分沙對河段水流沖刷強度的影響程度:

    α=F*fi/Ffi(2)

    通過計算三峽工程運用后荊江與上荊江河段的F*fi、α值,可發(fā)現(xiàn)三口分流分沙會導致河段水流沖刷強度呈一定比例減弱。荊江河段的α值為0.69~0.99,上荊江的α值為0.75~0.99(圖5)。故在后續(xù)分析中計算荊江與上荊江河段的水流沖刷強度時,采用F*fi表示受三口分流影響的河段水流沖刷強度。

    此外,累計沖淤量不僅與當前水沙條件相關,還受到前期多年水沙條件的滯后影響[16-17]。此處,前期水沙條件采用前n年平均水流沖刷強度Fnf表示:

    Fnf=1n∑ni=1F*fi(3)

    計算2003—2020年V與Fnf的相關系數(shù),結果表明:河段累計沖刷量與前期水流沖刷強度之間的相關程度隨滑動平均年數(shù)的不同而改變,二者相關系數(shù)在n=5時達到最大。故采用前5 a平均水流沖刷強度(F5f)表征河段來水來沙條件。

    2.1.2" 洞庭湖入?yún)R頂托作用

    荊江河段出口處有洞庭湖入?yún)R,河段內(nèi)沿程水位受到入?yún)R頂托作用的影響。已有研究采用河段出口的水位—流量關系反映頂托程度[10],但河段沖淤的實測資料一般為年沖淤量,本文的研究對象是荊江河段的年沖淤量及累計沖淤量,不易通過關系曲線分析頂托作用對河段沖淤的影響。此處采用洞庭湖入?yún)R(七里山站)與干流(監(jiān)利站)流量的比值定義匯流比(RQ)[18]:

    RQ=QQLS/QJL(4)

    式中:QQLS和QJL分別為七里山站和監(jiān)利站實測日均流量,m3/s。

    干流流量和入?yún)R流量每天都在發(fā)生變化,無法作為反映河段年沖淤量的參數(shù)。采用Mann-Kendall檢驗結果表明,三峽工程運用后年均匯流比的趨勢統(tǒng)計量基本處于95%置信區(qū)間(±1.96),幾乎沒有顯著的趨勢性變化(圖6),故年均匯流比無法準確反映洞庭湖入?yún)R頂托情況。而在特定流量級下,河段上下游水位落差與河段長度之比,即為河段的水面縱比降(量綱一數(shù)),可用來反映侵蝕基面的變化。因研究河段的長度為定值,故水位落差也能間接表示水面縱比降的大小。頂托作用導致侵蝕基面抬高,則水位落差減小,反之則水位落差增大,具有一定的物理意義且便于計算分析。

    入?yún)R頂托作用是干流下泄流量與洞庭湖入?yún)R共同影響的結果[16],而干流流量與入?yún)R流量均隨時間變化,在汛期和枯期會出現(xiàn)明顯的季節(jié)差異。為了準確反映洞庭湖入?yún)R頂托作用,此處以研究河段的年均水位落差(ΔZ)作為代表洞庭湖入?yún)R頂托程度的參數(shù)。ΔZ可表示為

    ΔZ=∑ni=1ΔZi·fi/nQ(5)

    式中:ΔZi和fi分別為各流量級的水位落差和流量級區(qū)間出現(xiàn)的概率;nQ為流量級區(qū)間的個數(shù)。因此式(5)考慮了干流不同流量級出現(xiàn)頻率的影響。

    在計算ΔZ時,首先根據(jù)荊江河段的流量變化范圍將每年枯水至洪水選取6個代表流量級區(qū)間(lt;7 000、7 000~lt;10 000、10 000~lt;20 000、20 000~lt;30 000、30 000~lt;40 000、≥40 000 m3/s);然后計算每個流量級的平均水位落差;最后統(tǒng)計干流(枝城站、監(jiān)利站)各流量級區(qū)間在研究時段內(nèi)每年出現(xiàn)的概率,以此作為權重計算得到ΔZ。

    下荊江河段靠近洞庭湖入?yún)R口,該河段水位明顯受到洞庭湖入?yún)R頂托作用[17-19]。對比水流沖刷強度相近年份的年均水位落差與年沖淤量,可發(fā)現(xiàn)入?yún)R頂托導致河段沖刷量減少的現(xiàn)象。首先由式(1)和式(5)計算2003—2020年下荊江河段的水流沖刷強度與年均水位落差,表1列舉了研究時段內(nèi)4組水流沖刷強度相近的年份,但年均水位落差不同;再對比不同水位落差條件下的年沖刷量,發(fā)現(xiàn)這些年份都有水位落差小則沖刷量少的定性規(guī)律。這表明ΔZ改變對河段年沖淤量的影響規(guī)律與河床演變學中出口侵蝕基面抬高導致沖刷量減少的規(guī)律相符。

    2.2" 河段沖淤量與上下游邊界條件的函數(shù)關系

    從河床演變基本原理的角度,進口水流沖刷強度與出口頂托程度分別影響河段的上下游邊界條件。由于荊江和洞庭湖之間存在復雜的相互影響,在建立2個影響因素與河段沖淤量之間的關系前,還需要證明2個影響因素的獨立性。對荊江及上、下荊江河段的F5f和ΔZ進行Shapiro-Wilk正態(tài)性檢驗(S-W檢驗)和Pearson相關性檢驗(表2)。結果表明,荊江及上、下荊江河段的F5f和ΔZ的S-W檢驗顯著性(p)均大于0.1,說明這些參數(shù)均服從正態(tài)分布;各河段的F5f和ΔZ的相關分析的顯著性均大于0.1,說明2個因素間無顯著相關性?;谏鲜龇治?,認為這2個影響河段沖淤的因素在統(tǒng)計意義上相互獨立,可以建立河段累計沖淤量與這2個因素的函數(shù)關系。

    三峽水庫清水下泄導致荊江河段的水沙條件改變,洞庭湖入?yún)R頂托導致荊江河段的侵蝕基面發(fā)生變化,荊江河段的沖淤過程是二者共同作用的結果,需要同時考慮上、下游邊界條件的影響[20-23]。由于河段累計沖淤量與河段前5 a平均水流沖刷強度、河段上下游平均水位落差呈非線性關系。故采用以下多項式回歸各河段V與F5f、ΔZ的關系:

    V=α1Fβ15f+α2ΔZβ2(6)

    式中:α1、α2為系數(shù);β1、β2為指數(shù)。采用三峽工程運用后(2003—2017年)荊江河段實測沖淤量率定系數(shù)與指數(shù)(表3),采用2018—2020年實測值作為驗證。

    為進一步驗證式(6)計算上、下荊江與荊江河段累計沖淤量的準確性。圖7繪制了2003—2020年上、下荊江和荊江河段累計沖淤量計算值與實測值的對比情況(2018—2020年實測值作為驗證)??梢娛剑?)計算得到的河段累計沖淤過程與實際過程符合較好,表明建立的公式在綜合考慮水沙條件和頂托作用的影響后,能較好地反映三峽工程運行以來荊江河段累計沖淤過程。需要說明的是,式(6)中考慮到河床調(diào)整對水沙條件改變的滯后響應對河床沖刷有較大影響,選取了前5 a平均水流沖刷強度作為水沙條件的代表參數(shù),但實際上不同年份的水沙情勢存在差異,因此會造成個別年份的計算值與實測值出現(xiàn)偏差。

    分析各河段率定得到的函數(shù)關系式,可以看出,式(6)不僅考慮了河段進口的來水來沙條件,還考慮了河段出口洞庭湖入?yún)R頂托作用對河段侵蝕基面的影響,符合河床演變學基本原理。各河段系數(shù)、指數(shù)的符號一致,表明各河段河床沖淤過程對水流沖刷強度和頂托作用的響應一致。各河段α1lt;0、β1>0,表明河床沖淤量(沖刷為負)基本與前期水流沖刷強度成正相關,即水流沖刷強度越強,河床沖刷量越大;α2lt;0、β2>0,表明河床沖刷量與河段水位落差呈正相關。由于水位落差與頂托程度呈負相關關系,故頂托程度越明顯,河床沖刷量越小,可見式(6)反映的情況與實際規(guī)律(表1)相符,說明公式具有一定物理意義。因懸沙含量沿程恢復與三口分流的影響,上、下荊江河段的水流沖刷強度存在一定差距,上荊江水流沖刷強度在21~203之間,而下荊江水流沖刷強度在2.3~29.1之間;同時上荊江和荊江河段率定的α1還包含三口分流分沙的影響,所以出現(xiàn)了α1遠小于α2的情況。

    上述對于荊江河段沖淤過程的定性分析及函數(shù)關系式具有普適性,但定量解釋及在其他河流的應用則需要以相應實測水沙資料與系數(shù)、指數(shù)的率定結果作為基礎。在研究受洞庭湖入?yún)R影響的荊江河段時,應同時考慮河段上游來水來沙條件和頂托作用對河段沖淤過程的綜合影響。

    3" 洞庭湖入?yún)R頂托作用對河段沖淤的影響

    洞庭湖入?yún)R頂托作用會抬高河段出口水位,減緩河段水面縱比降,導致河段上下游水位落差減小。為定量分析頂托作用對荊江河段沖淤的影響,此處先確定假設無洞庭湖入?yún)R條件下的臨界水位落差,再通過式(5)計算假設無洞庭湖入?yún)R條件下荊江河段的沖淤過程,基于此計算入?yún)R頂托的影響。無洞庭湖入?yún)R條件下臨界水位落差的計算方法如下:

    (1) 分析實測資料可發(fā)現(xiàn),相同干流流量條件下的ΔZ與RQ具有顯著的相關性(相關系數(shù)均大于0.80)。通過篩選出每年枯水至洪水6個代表流量級的ΔZ與RQ,并建立2個參數(shù)之間的函數(shù)關系,經(jīng)試算,下式對2個參數(shù)的擬合效果最好:

    ΔZ=aexp(bRQ)(7)

    (2) 由實測資料可知,洞庭湖入?yún)R最小RQ約為0.1。為了確定幾乎無入?yún)R條件下的ΔZ,對RQ取值0.1,以6個代表流量級在每年發(fā)生的頻率為權重,代入式(7)和式(5)即可得到每年無入?yún)R條件下的臨界ΔZ。

    將假設洞庭湖無入?yún)R時各年的臨界ΔZ、已率定的系數(shù)和指數(shù)代入式(6),得到無入?yún)R條件下2003年至后續(xù)年份的河段累計沖淤量。無入?yún)R條件下的累計沖淤量與實際累計沖淤量的差值即為入?yún)R頂托導致的累計沖淤量變化(ΔVi):

    ΔVi=ViW/O-ViACT(8)

    式中:ViW/O為假設無洞庭湖入?yún)R條件下,2003年至后續(xù)年份的河段累計沖淤量,下標i代表年份;ViACT為2003年至后續(xù)年份的河段實際累計沖淤量。

    ΔVi與ViACT的比值,表示洞庭湖入?yún)R頂托對河段沖淤的影響(Pi):

    Pi=ΔVi/ViACT(9)

    式中:下標i代表年份。

    式(6)可計算荊江河段2003年至后續(xù)年份的河段累計沖淤量,由于三峽水庫清水下泄導致荊江河段持續(xù)沖刷,河段累計沖刷量不斷增加,頂托作用對河段累計沖淤量的影響也會隨時間改變。本研究在分析頂托作用對沖淤的影響時,暫不分析ΔVi與Pi的年際變化,而是采用式(8)、式(9)計算2003—2020年頂托作用對河段沖淤量的累計影響(ΔV2020與P2020)。對于荊江河段,在假設無洞庭湖入?yún)R頂托的情況下,2003—2020年的河段累計沖刷量比實測值增加0.67億m3,計算得到入?yún)R頂托能使河段沖刷量減少約5%(圖8(a))。假設無入?yún)R條件下,同期上、下荊江河段比實際累計沖刷量增加約0.06和0.81億m3,即入?yún)R頂托能使兩河段沖刷量分別減少約1%和16%(圖8(b)、圖8(c))。對比入?yún)R頂托對上、下荊江河段沖淤的影響,可發(fā)現(xiàn)頂托作用對下荊江河段沖淤過程的影響更明顯,對上荊江河段沖淤過程的影響較小。

    4" 結" 論

    三峽工程運用后,荊江河段處于持續(xù)沖刷狀態(tài),基于三峽工程運用后(2003—2020年)荊江河段的水文泥沙及河段沖淤量資料,提出了考慮上下邊界影響的河床沖淤過程的定量計算方法,分析了洞庭湖入?yún)R頂托對荊江河段沖淤量的影響。主要結論如下:

    (1) 水沙條件與洞庭湖入?yún)R頂托是影響荊江河段沖淤過程的重要因素。當上游建壩后下泄清水時,壩下游河道將會產(chǎn)生沿程沖刷;當河道下游侵蝕基面抬升時,例如洞庭湖入?yún)R對荊江河段水位的頂托作用,河道沖刷量將相應減少。

    (2) 本研究建立了荊江河段累計沖淤量與進口水流沖刷強度、進出口水位落差的函數(shù)關系式,能較好地反映上下游邊界條件對河段沖淤過程的影響。

    (3) 在假設無洞庭湖入?yún)R頂托的情況下,2003—2020年的河段累計沖刷量比實際情況增加了0.67億m3,即入?yún)R頂托能使荊江河段沖刷量減少約5%。同期無頂托條件下,上、下荊江河段累計沖淤量比實際增加約0.06億和0.81億m3,即入?yún)R頂托能使兩河段沖刷量減少約1%和16%。頂托作用對下荊江沖淤的影響更顯著。

    參考文獻:

    [1]夏軍強,鄧珊珊,周美蓉,等.三峽工程運用對近期荊江段平灘河槽形態(tài)調(diào)整的影響[J].水科學進展,2016,27(3):385-391.(XIA J Q,DENG S S,ZHOU M R,et al.Effects of the Three Gorges Project operation on the recent variation in bankfull channel geometry of the Jingjiang reach[J].Advances in Water Science,2016,27(3):385-391.(in Chinese))

    [2]吳保生.沖積河流平灘流量的滯后響應模型[J].水利學報,2008,39(6):680-687.(WU B S.Delayed response model for bankfull discharge of alluvial rivers[J].Journal of Hydraulic Engineering,2008,39(6):680-687.(in Chinese))

    [3]許全喜,董炳江,袁晶,等.三峽工程運用后長江中下游河道沖刷特征及其影響[J].湖泊科學,2023,35(2):650-661.(XU Q X,DONG B J,YUAN J,et al.Scouring effect of the middle and lower reaches of the Yangtze River and its impact after the impoundment of the Three Gorges Project[J].Journal of Lake Sciences,2023,35(2):650-661.(in Chinese))

    [4]方春明,胡春宏,陳緒堅.三峽水庫運用對荊江三口分流及洞庭湖的影響[J].水利學報,2014,45(1):36-41.(FANG C M,HU C H,CHEN X J.Impacts of Three Georges Reservoir′s operation on outflow of the three outlets of Jingjiang River and Dongting Lake[J].Journal of Hydraulic Engineering,2014,45(1):36-41.(in Chinese))

    [5]盧金友,姚仕明.水庫群聯(lián)合作用下長江中下游江湖關系響應機制[J].水利學報,2018,49(1):36-46.(LU J Y,YAO S M.Response mechanism of the river and lakes in the middle and lower reaches of the Yangtze River under the combined effect of reservoir groups[J].Journal of Hydraulic Engineering,2018,49(1):36-46.(in Chinese))

    [6]許炯心.含沙量和懸沙粒徑變化對長江宜昌—漢口河段年沖淤量的影響[J].水科學進展,2006,17(1):67-73.(XU J X.Influence of variations in suspended sediment concentration and grain size on sediment deposition of Yichang—Hankou reach of Yangtze River[J].Advances in Water Science,2006,17(1):67-73.(in Chinese))

    [7]左利欽,陸永軍,王洪楊,等.荊江河段泥沙沖淤對三峽水庫汛期排沙的響應[J].水科學進展,2023,34(1):115-125.(ZUO L Q,LU Y J,WANG H Y,et al.Responses of erosion and deposition in the Jingjiang reach to sediment delivery of the Three Gorges Reservoir during flood season[J].Advances in Water Science,2023,34(1):115-125.(in Chinese))

    [8]張為,吳美琴,李思璇,等.三峽水庫蓄水后城陵磯至九江段河道沖淤調(diào)整機理[J].水科學進展,2020,31(2):162-171.(ZHANG W,WU M Q,LI S X,et al.Mechanism of adjustment of scouring and silting of Chenglingji—Jiujiang reach in the middle reaches of the Yangtze River after impoundment of the Three Gorges Dam[J].Advances in Water Science,2020,31(2):162-171.(in Chinese))

    [9]周煒興,孫昭華,周坤,等.三峽水庫蓄水前后長江中下游流量頻率分布特征及其對洪水造床作用的影響[J].湖泊科學,2022,34(2):616-629.(ZHOU W X,SUN Z H,ZHOU K,et al.Characteristics of discharge frequency and their effects on flood channel formation in the middle and lower reaches of the Yangtze River before/after the impoundment of the Three Gorges Reservoir[J].Journal of Lake Sciences,2022,34(2):616-629.(in Chinese))

    [10]周建銀,高菲,元媛,等.三峽水庫運行前后長江中下游干流及兩湖的徑流過程變化[J].湖泊科學,2023,35(2):696-708.(ZHOU J Y,GAO F,YUAN Y,et al.Runoff variation in the middle and lower branches of Yangtze River and the two lakes(Dongting and Poyang Lake) before and after Three Gorges Reservoir′s operation[J].Journal of Lake Sciences,2023,35(2):696-708.(in Chinese))

    [11]陳棟,渠庚,郭小虎,等.三峽建庫前后洞庭湖對下荊江的頂托與消落作用研究[J].工程科學與技術,2020,52(2):86-94.(CHEN D,QU G,GUO X H,et al.Study of the supporting and falling impact of Dongting Lake on the lower Jingjiang River before and after construction of Three Gorges Dam[J].Advanced Engineering Sciences,2020,52(2):86-94.(in Chinese))

    [12]柴元方,鄧金運,楊云平,等.長江中游荊江河段同流量—水位演化特征及驅(qū)動成因[J].地理學報,2021,76(1):101-113.(CHAI Y F,DENG J Y,YANG Y P,et al.Evolution characteristics and driving factors of the water level at the same discharge in the Jingjiang reach of Yangtze River[J].Acta Geographica Sinica,2021,76(1):101-113.(in Chinese))

    [13]孫昭華,周煒興,周坤,等.江湖水沙輸移與長江中下游造床流量的關系[J].水利學報,2021,52(5):521-534.(SUN Z H,ZHOU W X,ZHOU K,et al.Relationship between the characteristics of water-sediment transportation in river-lake system and the channel forming discharge of the Middle and Lower Yangtze River[J].Journal of Hydraulic Engineering,2021,52(5):521-534.(in Chinese))

    [14]黃草,陳葉華,李志威,等.洞庭湖區(qū)水系格局及連通性優(yōu)化[J].水科學進展,2019,30(5):661-672.(HUANG C,CHEN Y H,LI Z W,et al.Optimization of water system pattern and connectivity in the Dongting Lake area[J].Advances in Water Science,2019,30(5):661-672.(in Chinese))

    [15]XIA J Q,DENG S S,ZHOU M R,et al.Geomorphic response of the Jingjiang reach to the Three Gorges Project operation[J].Earth Surface Processes and Landforms,2017,42(6):866-876.

    [16]WU B S,ZHENG S,THORNE C R.A general framework for using the rate law to simulate morphological response to disturbance in the fluvial system[J].Progress in Physical Geography:Earth and Environment,2012,36(5):575-597.

    [17]LYU Y W,F(xiàn)AGHERAZZI S,ZHENG S,et al.Enhanced hysteresis of suspended sediment transport in response to upstream damming:an example of the Middle Yangtze River downstream of the Three Gorges Dam[J].Earth Surface Processes and Landforms,2020,45(8):1846-1859.

    [18]楊春瑞,鄧金運,陳立,等.長江中游通江湖泊對干流頂托作用變化規(guī)律[J].水科學進展,2024,35(1):98-111.(YANG C R,DENG J Y,CHEN L,et al.Study on the change of backwater effect of the connected lakes in the middle reaches of the Yangtze River[J].Advances in Water Science,2024,35(1):98-111.(in Chinese))

    [19]孫昭華,周歆玥,范杰瑋,等.考慮回水影響的河道水位—流量關系確定方法[J].水科學進展,2021,32(2):259-270.(SUN Z H,ZHOU X Y,F(xiàn)AN J W,et al.Stage—discharge rating method considering backwater effect in river channel[J].Advances in Water Science,2021,32(2):259-270.(in Chinese))

    [20]ZHOU M R,XIA J Q,SHI X,et al.Hydrologic and geomorphic effects on the reduction of channel discharge capacity in the Middle Yangtze River[J].Science of the Total Environment,2023,866:161279.

    [21]林芬芬,夏軍強,周美蓉,等.下荊江河槽形態(tài)及過流能力調(diào)整對上下游邊界條件的響應[J].水利學報,2019,50(5):641-649.(LIN F F,XIA J Q,ZHOU M R,et al.Recent response of bankfull channel geometry and flood-discharge capacity in the lower Jingjiang reach to upstream and downstream boundary conditions[J].Journal of Hydraulic Engineering,2019,50(5):641-649.(in Chinese))

    [22]程亦菲,夏軍強,周美蓉,等.黃河下游不同河段分組懸沙輸移對河床沖淤的影響[J].水科學進展,2022,33(3):506-517.(CHENG Y F,XIA J Q,ZHOU M R,et al.Effects of grouped suspended sediment transport on channel evolution in the Lower Yellow River[J].Advances in Water Science,2022,33(3):506-517.(in Chinese))

    [23]XIA J Q,ZHOU M R,LIN F F,et al.Variation in reach-scale bankfull discharge of the Jingjiang reach undergoing upstream and downstream boundary controls[J].Journal of Hydrology,2017,547:534-543.

    Effects of the incoming flow-sediment regime and the backwater effect caused by

    the Dongting Lake confluence on channel degradation in the Jingjiang reach

    The study is financially supported by the National Natural Science Foundation of China (No.U2040215;No.U2240206).

    SHANG Haixin1,HU Chunhong1,2,XIA Junqiang1,ZHOU Meirong1

    (1. State Key Laboratory of Water Resources Engineering and Management,Wuhan University,Wuhan 430072,China;

    2. State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin,China Institute of

    Water Resources and Hydropower Research,Beijing 100048,China)

    Abstract:The Jingjiang reach (JR) has experienced continuous and significant channel degradation after the Three Gorges Project operation.It is urgent to quantify the backwater effect caused by the confluence of the Dongting Lake,and the consequent impact of incoming flow-sediment regime on the channel degradation in the JR.Based on the measured hydrological data and data for channel degradation from 2003 to 2020,an empirical function for clarifying the relationship between the cumulative channel evolution amount and boundary conditions for the upstream and downstream of the studied JR reach has been established.In addition,the contribution of the backwater effect to the cumulative channel evolution amount was investigated.The results indicate that:① The cumulative channel evolution amount in JR was positively correlated with the fluvial erosion intensity at the inlet,and it was inversely with the corresponding difference between the average water stages at the inlet and outlet.Therefore,the channel degradation in JR is the result of these two influencing factors.② The empirical function derived could effectively represent the impact of the fluvial erosion intensity and differences between the average water stages of channel inlet and outlet on the cumulative channel evolution amount in JR (R2gt;0.92).③ The cumulative channel evolution amount of the JR section from 2003 to 2020 was 1.23 billion cubic meters.Under hypothetical conditions for the absence of the backwater effect caused by the Dongting Lake confluence,the cumulative channel evolution volume would increase by 67 million cubic meters.This indicated that a potential reduction in the cumulative channel evolution amount of nearly 5% was attributable to the backwater effect.At the same circumstance,the cumulative channel evolution volume in the upstream and downstream JR was 0.73 and 0.5 billion cubic meters,respectively.The backwater effect caused a reduction of the two values by about 1% and 16%,respectively.Therefore,the backwater effect decreased the channel evolution volume in JR and its influence on the channel degradation was more obvious in the downstream than that in the upstream.

    Key words:backwater effect;channel evolution;fluvial erosion intensity;confluence of the Dongting Lake;Jingjiang reach

    av片东京热男人的天堂| 伊人久久大香线蕉亚洲五| 在线播放国产精品三级| 免费搜索国产男女视频| 五月伊人婷婷丁香| 小说图片视频综合网站| 一本精品99久久精品77| e午夜精品久久久久久久| 两性夫妻黄色片| 老鸭窝网址在线观看| 他把我摸到了高潮在线观看| 在线a可以看的网站| 757午夜福利合集在线观看| 精品久久久久久久人妻蜜臀av| 亚洲一区二区三区不卡视频| 色av中文字幕| 欧美黑人巨大hd| 欧美午夜高清在线| 欧美精品亚洲一区二区| 亚洲天堂国产精品一区在线| 国产av又大| 欧美一区二区国产精品久久精品 | 日韩精品青青久久久久久| 在线观看舔阴道视频| 一级黄色大片毛片| 午夜福利视频1000在线观看| 我要搜黄色片| 九九热线精品视视频播放| АⅤ资源中文在线天堂| 嫩草影视91久久| 手机成人av网站| netflix在线观看网站| 桃色一区二区三区在线观看| 最近最新中文字幕大全免费视频| 国产不卡一卡二| 午夜福利18| 欧美+亚洲+日韩+国产| 久久久久久免费高清国产稀缺| 岛国在线观看网站| 国产高清videossex| 欧美成狂野欧美在线观看| 国产一区二区激情短视频| 此物有八面人人有两片| 欧美乱色亚洲激情| 精品午夜福利视频在线观看一区| 久久亚洲真实| 长腿黑丝高跟| 村上凉子中文字幕在线| 国产亚洲精品一区二区www| 国产精品 欧美亚洲| 真人一进一出gif抽搐免费| 少妇裸体淫交视频免费看高清 | 亚洲av美国av| av中文乱码字幕在线| 老司机在亚洲福利影院| 黄色a级毛片大全视频| 日本撒尿小便嘘嘘汇集6| 99精品久久久久人妻精品| 极品教师在线免费播放| 床上黄色一级片| 免费av毛片视频| 啦啦啦免费观看视频1| 巨乳人妻的诱惑在线观看| 精品少妇一区二区三区视频日本电影| 国产精品 欧美亚洲| 国产真实乱freesex| 欧美成人性av电影在线观看| 国产91精品成人一区二区三区| 好男人电影高清在线观看| 不卡av一区二区三区| 久久久久久久久免费视频了| 无人区码免费观看不卡| 精品免费久久久久久久清纯| 亚洲av美国av| 少妇人妻一区二区三区视频| 亚洲中文字幕一区二区三区有码在线看 | 美女大奶头视频| 久久婷婷成人综合色麻豆| 熟女电影av网| 女同久久另类99精品国产91| 99久久精品国产亚洲精品| 狠狠狠狠99中文字幕| 国产免费av片在线观看野外av| 在线国产一区二区在线| 国模一区二区三区四区视频 | 亚洲美女视频黄频| 亚洲熟妇中文字幕五十中出| 色噜噜av男人的天堂激情| 韩国av一区二区三区四区| 欧美日韩国产亚洲二区| 怎么达到女性高潮| 亚洲一区二区三区不卡视频| 成人av一区二区三区在线看| 成人三级黄色视频| 国产精品98久久久久久宅男小说| 大型av网站在线播放| 天堂√8在线中文| 欧美大码av| 午夜精品在线福利| 日本免费a在线| 久久久水蜜桃国产精品网| 日本 av在线| 亚洲av成人不卡在线观看播放网| 少妇人妻一区二区三区视频| 美女免费视频网站| 国产亚洲精品第一综合不卡| 欧美日韩黄片免| 在线观看免费视频日本深夜| 麻豆成人av在线观看| 国产精品免费视频内射| 国产亚洲精品久久久久久毛片| 女同久久另类99精品国产91| 超碰成人久久| 久久久久久久久免费视频了| 精品久久久久久久久久久久久| 国产成人啪精品午夜网站| 成人亚洲精品av一区二区| 精品日产1卡2卡| 日韩av在线大香蕉| 亚洲精品美女久久久久99蜜臀| 99久久精品热视频| 波多野结衣高清无吗| 亚洲人成网站在线播放欧美日韩| 欧美乱码精品一区二区三区| 成人18禁在线播放| 特大巨黑吊av在线直播| 欧美日韩福利视频一区二区| 久9热在线精品视频| 国内精品久久久久久久电影| 可以免费在线观看a视频的电影网站| 校园春色视频在线观看| 亚洲第一欧美日韩一区二区三区| 精品欧美国产一区二区三| 69av精品久久久久久| 91成年电影在线观看| 欧美黑人欧美精品刺激| 亚洲成人中文字幕在线播放| 日本免费一区二区三区高清不卡| 国产精品亚洲一级av第二区| 久久久久国产精品人妻aⅴ院| 久久 成人 亚洲| 最新美女视频免费是黄的| 人妻久久中文字幕网| 免费在线观看黄色视频的| 真人一进一出gif抽搐免费| 国产亚洲av嫩草精品影院| 久久久久九九精品影院| 亚洲精华国产精华精| 国产高清有码在线观看视频 | 午夜亚洲福利在线播放| 久久人妻福利社区极品人妻图片| av视频在线观看入口| 亚洲精华国产精华精| 最近最新中文字幕大全免费视频| 国产亚洲精品久久久久5区| 国产主播在线观看一区二区| 国内揄拍国产精品人妻在线| 久久99热这里只有精品18| www.999成人在线观看| 国产精品av久久久久免费| 亚洲九九香蕉| 中亚洲国语对白在线视频| 精品欧美国产一区二区三| 在线看三级毛片| 国产av一区在线观看免费| 可以在线观看毛片的网站| 一级毛片高清免费大全| 国产av麻豆久久久久久久| 男男h啪啪无遮挡| 99国产精品一区二区三区| 亚洲专区国产一区二区| 老熟妇仑乱视频hdxx| 两个人的视频大全免费| 精品少妇一区二区三区视频日本电影| 757午夜福利合集在线观看| 亚洲专区字幕在线| netflix在线观看网站| 人人妻人人看人人澡| 日日夜夜操网爽| www日本在线高清视频| 身体一侧抽搐| av国产免费在线观看| 最好的美女福利视频网| 男人舔奶头视频| 嫩草影视91久久| 嫩草影视91久久| 国产三级中文精品| 日本免费一区二区三区高清不卡| 我的老师免费观看完整版| 国产一区二区三区视频了| 国产亚洲av高清不卡| 黄色 视频免费看| 2021天堂中文幕一二区在线观| av天堂在线播放| 丰满的人妻完整版| 三级男女做爰猛烈吃奶摸视频| 国产精品乱码一区二三区的特点| 日韩中文字幕欧美一区二区| 全区人妻精品视频| 国产一区二区三区视频了| 一区二区三区高清视频在线| 日本免费一区二区三区高清不卡| 亚洲精品一区av在线观看| 中文字幕人成人乱码亚洲影| av在线播放免费不卡| 久久中文看片网| 黄片小视频在线播放| 长腿黑丝高跟| 欧美日韩乱码在线| 色精品久久人妻99蜜桃| 成人午夜高清在线视频| 国内揄拍国产精品人妻在线| 国产精品一区二区精品视频观看| 午夜福利成人在线免费观看| 免费人成视频x8x8入口观看| 欧美午夜高清在线| 亚洲精品中文字幕在线视频| 免费观看精品视频网站| 黄色片一级片一级黄色片| 国内毛片毛片毛片毛片毛片| 日韩 欧美 亚洲 中文字幕| 在线a可以看的网站| 波多野结衣巨乳人妻| 久久久精品国产亚洲av高清涩受| 成人国语在线视频| 日本在线视频免费播放| 久久 成人 亚洲| 国产乱人伦免费视频| 亚洲欧美日韩无卡精品| 大型黄色视频在线免费观看| 又粗又爽又猛毛片免费看| 国产亚洲精品第一综合不卡| 97人妻精品一区二区三区麻豆| 免费在线观看完整版高清| 国产91精品成人一区二区三区| 国产高清视频在线观看网站| 麻豆国产97在线/欧美 | 亚洲狠狠婷婷综合久久图片| 麻豆av在线久日| 欧美乱妇无乱码| 欧美中文综合在线视频| 制服人妻中文乱码| 超碰成人久久| 男女午夜视频在线观看| 精品电影一区二区在线| 一进一出抽搐gif免费好疼| 无遮挡黄片免费观看| 久久99热这里只有精品18| 波多野结衣高清作品| 女生性感内裤真人,穿戴方法视频| 亚洲欧美一区二区三区黑人| 99riav亚洲国产免费| 国产成人一区二区三区免费视频网站| 婷婷丁香在线五月| 一级毛片女人18水好多| 天天一区二区日本电影三级| 麻豆久久精品国产亚洲av| 草草在线视频免费看| 久久中文字幕一级| 在线观看日韩欧美| 久久草成人影院| 亚洲 欧美一区二区三区| 亚洲一区中文字幕在线| 老司机在亚洲福利影院| 免费看a级黄色片| 亚洲精品美女久久av网站| 国产成人精品久久二区二区91| x7x7x7水蜜桃| 国产av麻豆久久久久久久| 身体一侧抽搐| 国产精品,欧美在线| а√天堂www在线а√下载| 很黄的视频免费| 亚洲av中文字字幕乱码综合| 不卡av一区二区三区| 18禁观看日本| 欧美午夜高清在线| 99久久无色码亚洲精品果冻| 特大巨黑吊av在线直播| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲va日本ⅴa欧美va伊人久久| 久久精品91蜜桃| 天天躁狠狠躁夜夜躁狠狠躁| 国产一区二区三区在线臀色熟女| 免费在线观看亚洲国产| 欧美乱色亚洲激情| 人人妻,人人澡人人爽秒播| 国产蜜桃级精品一区二区三区| 久久久久久久午夜电影| 国产成人精品久久二区二区91| 国产欧美日韩一区二区精品| 色综合欧美亚洲国产小说| 精品国内亚洲2022精品成人| 麻豆成人av在线观看| 村上凉子中文字幕在线| 国产一区二区激情短视频| 好男人在线观看高清免费视频| 最近最新中文字幕大全电影3| 色老头精品视频在线观看| 美女午夜性视频免费| 久久久久久久精品吃奶| 桃色一区二区三区在线观看| 免费看日本二区| 男男h啪啪无遮挡| 精品国产超薄肉色丝袜足j| 午夜a级毛片| 在线观看美女被高潮喷水网站 | 久久中文看片网| 18美女黄网站色大片免费观看| 99热这里只有精品一区 | 母亲3免费完整高清在线观看| 三级毛片av免费| 成人18禁高潮啪啪吃奶动态图| 韩国av一区二区三区四区| 999久久久精品免费观看国产| 2021天堂中文幕一二区在线观| 男女床上黄色一级片免费看| 国产高清有码在线观看视频 | 一级毛片女人18水好多| 亚洲午夜精品一区,二区,三区| 亚洲精品在线观看二区| or卡值多少钱| 真人一进一出gif抽搐免费| 夜夜躁狠狠躁天天躁| 黄色 视频免费看| 村上凉子中文字幕在线| 国产精品,欧美在线| 成人午夜高清在线视频| 精品日产1卡2卡| 成人高潮视频无遮挡免费网站| 婷婷精品国产亚洲av在线| АⅤ资源中文在线天堂| 一级片免费观看大全| 亚洲成a人片在线一区二区| 啦啦啦观看免费观看视频高清| av免费在线观看网站| 国产精品久久久久久亚洲av鲁大| 国产免费av片在线观看野外av| 一夜夜www| xxxwww97欧美| 看片在线看免费视频| 在线观看免费视频日本深夜| 国产一区在线观看成人免费| 嫩草影院精品99| 操出白浆在线播放| 淫妇啪啪啪对白视频| 久久国产精品影院| 国产一区在线观看成人免费| 欧美性猛交╳xxx乱大交人| 精品熟女少妇八av免费久了| 九色国产91popny在线| 亚洲欧美精品综合一区二区三区| 国产成人系列免费观看| 国产精品自产拍在线观看55亚洲| 午夜精品在线福利| 最新美女视频免费是黄的| 成年版毛片免费区| 亚洲成人中文字幕在线播放| 国产精品精品国产色婷婷| 亚洲免费av在线视频| 亚洲av片天天在线观看| 精品人妻1区二区| 99国产综合亚洲精品| 99久久国产精品久久久| 人人妻人人看人人澡| 久久国产乱子伦精品免费另类| 国产黄a三级三级三级人| 欧美乱妇无乱码| 50天的宝宝边吃奶边哭怎么回事| 亚洲中文字幕日韩| 男女下面进入的视频免费午夜| 久久香蕉国产精品| 精品欧美国产一区二区三| 日本一本二区三区精品| 成人特级黄色片久久久久久久| 色综合亚洲欧美另类图片| 国产精品久久久久久人妻精品电影| 免费在线观看成人毛片| 真人做人爱边吃奶动态| 99久久久亚洲精品蜜臀av| 久久精品影院6| 免费在线观看黄色视频的| 精品国内亚洲2022精品成人| 老司机靠b影院| 黄色视频不卡| 国产成人av教育| 日本熟妇午夜| 午夜免费观看网址| 制服诱惑二区| 人人妻人人看人人澡| 搡老岳熟女国产| 免费搜索国产男女视频| 成人亚洲精品av一区二区| 天天一区二区日本电影三级| 国产主播在线观看一区二区| 人成视频在线观看免费观看| 成人亚洲精品av一区二区| av福利片在线观看| 亚洲精品美女久久久久99蜜臀| 无人区码免费观看不卡| 两性午夜刺激爽爽歪歪视频在线观看 | 精品久久久久久,| 午夜福利视频1000在线观看| 国产精品电影一区二区三区| 最好的美女福利视频网| 免费人成视频x8x8入口观看| 黄色丝袜av网址大全| 身体一侧抽搐| 久久人妻福利社区极品人妻图片| 午夜老司机福利片| 中文资源天堂在线| 亚洲真实伦在线观看| 人成视频在线观看免费观看| 亚洲av成人精品一区久久| 大型黄色视频在线免费观看| cao死你这个sao货| 高潮久久久久久久久久久不卡| 成人av一区二区三区在线看| 天堂√8在线中文| 黄色片一级片一级黄色片| 亚洲成av人片在线播放无| 在线视频色国产色| 午夜亚洲福利在线播放| 亚洲精品国产一区二区精华液| 变态另类成人亚洲欧美熟女| 在线永久观看黄色视频| 日韩高清综合在线| 人妻丰满熟妇av一区二区三区| 18禁观看日本| 国内毛片毛片毛片毛片毛片| 成人三级做爰电影| 精品少妇一区二区三区视频日本电影| 亚洲av美国av| 美女 人体艺术 gogo| 欧美性长视频在线观看| 国产av一区二区精品久久| 日本一区二区免费在线视频| 国产精品一区二区三区四区免费观看 | 男女做爰动态图高潮gif福利片| 国产亚洲精品久久久久久毛片| 国产在线观看jvid| 少妇被粗大的猛进出69影院| 精品久久久久久久末码| 欧美日韩国产亚洲二区| 国产精品免费一区二区三区在线| 欧美黄色淫秽网站| 999久久久国产精品视频| 欧美黑人巨大hd| 日韩欧美在线乱码| 久久精品国产清高在天天线| 男插女下体视频免费在线播放| 脱女人内裤的视频| 久久久国产欧美日韩av| 精品久久久久久久末码| 99精品在免费线老司机午夜| 色在线成人网| 日韩大码丰满熟妇| 国产成人精品久久二区二区免费| 国产一区二区在线观看日韩 | 在线看三级毛片| 老司机在亚洲福利影院| 一二三四社区在线视频社区8| 久久精品国产99精品国产亚洲性色| 黑人巨大精品欧美一区二区mp4| 久久久久久久久久黄片| 热99re8久久精品国产| av免费在线观看网站| 欧美一区二区国产精品久久精品 | 麻豆一二三区av精品| 日韩精品免费视频一区二区三区| 欧美黑人巨大hd| 亚洲av成人精品一区久久| 97碰自拍视频| 欧美黑人巨大hd| 亚洲av片天天在线观看| 国产又色又爽无遮挡免费看| 女警被强在线播放| 国产精品一区二区三区四区免费观看 | 人人妻人人看人人澡| 精品久久久久久,| 久久性视频一级片| 又黄又粗又硬又大视频| 在线免费观看的www视频| 亚洲成av人片免费观看| 真人做人爱边吃奶动态| 国产又色又爽无遮挡免费看| 国产真实乱freesex| 超碰成人久久| www.999成人在线观看| 精品熟女少妇八av免费久了| 天天一区二区日本电影三级| 亚洲18禁久久av| 国产单亲对白刺激| 999精品在线视频| 亚洲av熟女| 欧美最黄视频在线播放免费| 午夜福利欧美成人| 国产av又大| 成人国产综合亚洲| 黄片小视频在线播放| 久久久久国内视频| 亚洲av成人一区二区三| 巨乳人妻的诱惑在线观看| 男女之事视频高清在线观看| 黑人欧美特级aaaaaa片| 美女大奶头视频| 操出白浆在线播放| 麻豆成人av在线观看| 欧美日韩亚洲综合一区二区三区_| 精品久久久久久久久久免费视频| 丁香欧美五月| 国内精品久久久久久久电影| 变态另类成人亚洲欧美熟女| 国产成年人精品一区二区| 国产久久久一区二区三区| 日本五十路高清| 免费在线观看黄色视频的| 亚洲一区二区三区色噜噜| 精品久久久久久久毛片微露脸| 成人国语在线视频| 国产精品久久久久久久电影 | 色综合婷婷激情| 欧美黄色片欧美黄色片| 桃色一区二区三区在线观看| 成人特级黄色片久久久久久久| 在线观看美女被高潮喷水网站 | 99精品久久久久人妻精品| 国产精品99久久99久久久不卡| 香蕉丝袜av| 中文字幕高清在线视频| 国产精品一区二区三区四区久久| 日韩欧美在线二视频| 美女黄网站色视频| 亚洲免费av在线视频| 好男人电影高清在线观看| 成年人黄色毛片网站| 麻豆成人av在线观看| cao死你这个sao货| 又黄又粗又硬又大视频| 淫秽高清视频在线观看| 亚洲 欧美一区二区三区| 美女高潮喷水抽搐中文字幕| 国产av一区在线观看免费| 美女午夜性视频免费| 国产日本99.免费观看| 精品少妇一区二区三区视频日本电影| 一个人免费在线观看电影 | 久久久久国产一级毛片高清牌| 91字幕亚洲| 国产欧美日韩一区二区精品| 久久精品91蜜桃| 91国产中文字幕| 久久香蕉精品热| 丝袜人妻中文字幕| 变态另类成人亚洲欧美熟女| 久久亚洲真实| 精品福利观看| av超薄肉色丝袜交足视频| av免费在线观看网站| av超薄肉色丝袜交足视频| 久久久久国产精品人妻aⅴ院| 禁无遮挡网站| 国产精品99久久99久久久不卡| 香蕉丝袜av| 波多野结衣巨乳人妻| 久久久水蜜桃国产精品网| 日本一本二区三区精品| 99久久精品热视频| 美女扒开内裤让男人捅视频| 欧美av亚洲av综合av国产av| 老司机午夜十八禁免费视频| 麻豆久久精品国产亚洲av| 男插女下体视频免费在线播放| 在线观看66精品国产| 欧美在线一区亚洲| 在线观看舔阴道视频| 午夜激情av网站| 69av精品久久久久久| 在线十欧美十亚洲十日本专区| 国内毛片毛片毛片毛片毛片| 99国产精品99久久久久| 国产一区二区三区在线臀色熟女| 国产激情偷乱视频一区二区| 国产一区二区激情短视频| 久久午夜综合久久蜜桃| 免费av毛片视频| 欧美极品一区二区三区四区| 美女午夜性视频免费| 熟女少妇亚洲综合色aaa.| 午夜福利欧美成人| 一卡2卡三卡四卡精品乱码亚洲| 国内揄拍国产精品人妻在线| 欧洲精品卡2卡3卡4卡5卡区| avwww免费| 精品久久久久久成人av| 国产主播在线观看一区二区| 久久伊人香网站| 日韩欧美一区二区三区在线观看| 在线播放国产精品三级| 五月伊人婷婷丁香| 精品欧美一区二区三区在线| 亚洲国产看品久久| 成人三级做爰电影| 两个人免费观看高清视频| 高清在线国产一区| 久久精品国产清高在天天线| 黑人欧美特级aaaaaa片| 国产精品一区二区三区四区久久| 国产成人av激情在线播放| 亚洲精品粉嫩美女一区| 亚洲av成人一区二区三| 制服丝袜大香蕉在线| 曰老女人黄片| 一夜夜www| 久久久久国产精品人妻aⅴ院|