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

    長江中游戴家洲河段灘群聯(lián)動(dòng)演變及其對航道工程的響應(yīng)機(jī)制

    2024-01-01 00:00:00楊云平鄭金海王建軍劉萬利朱玲玲
    水科學(xué)進(jìn)展 2024年5期

    摘要: 航道淺灘演變與灘槽形態(tài)調(diào)整具有強(qiáng)關(guān)聯(lián)性,識(shí)別灘槽形態(tài)、河床沖淤與航道工程的響應(yīng)關(guān)系,對評(píng)估航道工程實(shí)施效果具有重要意義。以長江中游戴家洲河段為對象,采用1954—2022 年期間實(shí)測水沙和地形資料,解析灘群演變聯(lián)動(dòng)關(guān)系及對航道工程的響應(yīng)機(jī)制。結(jié)果表明:①自然水沙階段(2008 年以前),巴河邊灘與池湖港邊灘關(guān)聯(lián)性強(qiáng),池湖港邊灘淤漲并減弱了燕磯節(jié)點(diǎn)挑流作用,加速了戴家洲洲頭低灘上延,而巴河邊灘淤漲使得戴家洲洲頭低灘沖刷,但流域來沙量減少增加了洲頭低灘沖刷速率;戴家洲汊道分流關(guān)系不穩(wěn)定且漢口站相同流量條件下的右汊分流比為減少態(tài)勢,來沙量低的大水年份加速戴家洲右緣蝕退并為樂家灣邊灘淤寬提供了空間。②航道工程影響階段(2008 年以來),航道工程實(shí)施使得戴家洲洲頭低灘上延,其右緣沖刷應(yīng)引起重視;戴家洲汊道進(jìn)口段航道工程壓縮了河道橫向空間并加速了巴河邊灘沖失趨勢,池湖港邊灘側(cè)蝕使得分汊段水流分散程度增加,不利于汊道進(jìn)口段淺灘水動(dòng)力的穩(wěn)定或增強(qiáng);戴家洲右緣洲體守護(hù)后的完整性增強(qiáng),樂家灣邊灘由淤漲轉(zhuǎn)為沖刷過程形成的分離體進(jìn)入航槽并礙航,6m 水深航道工程實(shí)施后樂家灣邊灘完整性增強(qiáng),右汊中下段淺區(qū)航道條件改善。

    關(guān)鍵詞: 分流關(guān)系;灘群聯(lián)動(dòng);航道工程;戴家洲河段;長江中游

    中圖分類號(hào): TV148 文獻(xiàn)標(biāo)志碼: A 文章編號(hào): 1001-6791( 2024) 05-0830-13

    長江中下游分汊河段眾多,自然條件下汊道內(nèi)邊心灘沖淤交替,主流隨河床沖淤在多汊間擺動(dòng),不利于航槽邊界穩(wěn)定。流域來沙減少背景下,汊道間河床出現(xiàn)不均衡沖刷,洲灘沖刷后河道展寬[1-2],河槽水動(dòng)力相對減弱后易于泥沙落淤,在放寬段形成礙航淺灘[3-4]。三峽工程運(yùn)行后,長江中下游河段邊灘和江心洲總面積為減少態(tài)勢[5-6],洲灘之間普遍存在聯(lián)動(dòng)性演變特征,如沙市河段[1]、湖廣—羅湖洲河段[7]、口岸直水道[3-4]、福姜沙水道[8-9] 等較為典型。沙市河段洲灘聯(lián)動(dòng)演變中,自然演變階段以臘林洲和三八灘為主導(dǎo),航道工程實(shí)施后以太平口心灘演變?yōu)橹鲗?dǎo)[1],且太平口心灘的萎縮進(jìn)一步引起下游三八灘汊道分流關(guān)系不穩(wěn)定。湖廣—羅湖洲水道灘群聯(lián)動(dòng)演變與水情過程、節(jié)點(diǎn)挑流作用等有關(guān),在航道工程實(shí)施后灘群聯(lián)動(dòng)演變關(guān)系減弱[7]。口岸直水道落成洲和嘶馬彎道的演變及穩(wěn)定性決定著進(jìn)口段灘群聯(lián)動(dòng)關(guān)系,三峽工程運(yùn)行后落成洲沖刷為三益橋邊灘發(fā)育提供了空間條件,在進(jìn)口段形成了交錯(cuò)型的三益橋淺灘[3-4]。福姜沙水道內(nèi)周期性演變的靖江邊灘是灘群聯(lián)動(dòng)演變、汊道分流比關(guān)系調(diào)整的主因,該灘體切割的沙體體積也決定著下游福北水道礙航淤積量大小[8-9]。綜上分析認(rèn)為,沖積型河流灘群演變存在較強(qiáng)的關(guān)聯(lián)性,在航道工程作用下洲灘聯(lián)動(dòng)演變模式出現(xiàn)轉(zhuǎn)變,這一模式調(diào)整對航道淺灘條件的影響需重點(diǎn)關(guān)注。

    戴家洲河段位于長江中游,河段內(nèi)分布有巴河邊灘、池湖港邊灘、戴家洲、樂家灣邊灘等易變洲灘,屬于礙航嚴(yán)重的灘段。同時(shí),戴家洲兩汊分流比接近,作為枯水期主航道直水道進(jìn)口與汊內(nèi)灘槽穩(wěn)定性差,航道條件不穩(wěn)定且存在多處淺區(qū)[10]。在戴家洲河段航道治理初期,通過分流關(guān)系與灘槽格局等分析,提出了航道治理的目標(biāo)河型[11],確定了右汊為枯水期通航主汊道[10],并采用數(shù)學(xué)模型模擬識(shí)別了洲頭分流面位置,支撐了洲頭工程平面布置[12]。戴家洲河段已實(shí)施了4 期航道工程:①一期的戴家洲魚骨壩工程實(shí)施后,穩(wěn)定了戴家洲水道枯水期分流關(guān)系,右汊彎道形態(tài)轉(zhuǎn)優(yōu),進(jìn)口段河寬受控且航道條件改善[13];②二期的戴家洲右緣下段守護(hù)工程實(shí)施后,洲灘邊界條件得到初步穩(wěn)定[14];③三期的戴家洲右汊左岸壩體工程實(shí)施后,右汊中上段航道條件轉(zhuǎn)好且中下段4.5 m 槽貫通[15];④四期的武安段航道工程實(shí)施后,池湖港邊灘和樂家灣邊灘沖刷態(tài)勢得到控制,形成了良好的灘槽形態(tài),直水道航道條件得到改善,基本實(shí)現(xiàn)了6 m 水深貫通[15]。綜上分析認(rèn)為,已有研究關(guān)注了戴家洲河段河勢演變、灘槽形態(tài)、汊道分流關(guān)系、礙航特性等方面,未能系統(tǒng)解析戴家洲河段灘槽形態(tài)演變與航道條件的聯(lián)動(dòng)關(guān)系,為做好已建航道工程效果評(píng)估及航道條件穩(wěn)定性分析工作,應(yīng)深化灘群演變、汊道分流關(guān)系及礙航特性等關(guān)聯(lián)性研究。

    本文以長江中游戴家洲河段為對象,采用1954—2022 年期間的實(shí)測水沙數(shù)據(jù)、河床地形資料,分析灘槽形態(tài)演變與聯(lián)動(dòng)特征,識(shí)別航道工程對洲灘形態(tài)的影響。研究結(jié)果可加深對彎曲分汊河段灘槽演變及礙航特性的認(rèn)識(shí),也可為航道規(guī)劃提供理論參考。

    1 研究區(qū)域與研究方法

    1.1 研究區(qū)域

    戴家洲河段上段的鄂州市至燕磯鎮(zhèn)段為順直放寬河型的巴河水道(長度為14 km),下段燕磯鎮(zhèn)至迴風(fēng)磯段為彎曲分汊河型的戴家洲水道(長度為20 km)(圖1(a))。巴河水道下段左岸為巴河邊灘,右岸為池湖港邊灘;戴家洲水道內(nèi)江心洲為新洲和戴家洲(已合并為戴家洲),右汊內(nèi)右岸側(cè)為樂家灣邊灘。2008 年以來,戴家洲河段已實(shí)施了4 期航道工程, 提升了航道尺度( 圖1( b) ) 。1968 —2009 年航道維護(hù)水深為4.0 m,2008 年以來實(shí)施了戴家洲洲頭守護(hù)工程( 2009—2010 年)、戴家洲右緣下段守護(hù)工程( 2010—2013 年)和右汊右岸側(cè)壩體工程( 2012—2014 年),航道尺度實(shí)現(xiàn)了4.5 m×200 m×1 050 m;2018 年以來實(shí)施了武漢至安慶段航道工程,至2021 年航道維護(hù)水深為5.0 m,2022 —2023 年航道尺度實(shí)現(xiàn)了6.0 m×200 m×1 050 m(大洪水年份為110 m)。

    1.2 研究資料與來源

    收集1954—2022 年漢口站日均流量、輸沙量數(shù)據(jù),用于分析流域來水來沙特征,以及河床沖淤強(qiáng)度、洲灘面積與水沙條件的關(guān)系;收集1957—2022 年戴家洲汊道分流比,目的是比較同流量枯水分流比的階段性特征;收集1958—2021 年共45 套地形數(shù)據(jù),其中2002 年以前為等深線數(shù)據(jù),用于分析洲灘形態(tài)演變與航道條件,2002 年以來為1∶10 000 實(shí)測地形,用于分析河床沖淤分布、洲灘形態(tài)及航道條件,洲灘面積計(jì)算采用等深線數(shù)據(jù),保證了2002 年前后數(shù)據(jù)的一致性和可靠性。以上實(shí)測數(shù)據(jù),主要來自于長江水文年鑒以及長江航道局、長江航道測量中心等機(jī)構(gòu)。

    1.3 灘群聯(lián)動(dòng)關(guān)系的分析方法

    河道邊灘和江心洲分布方式分為并聯(lián)和串聯(lián)2 種,并聯(lián)方式為邊灘處于河道兩岸或多級(jí)分汊段江心洲位置橫向排列[2-3],易形成交錯(cuò)淺灘或正常淺灘;串聯(lián)方式為邊灘或江心洲以上下游形式排列且相對距離較近[16-17],易形成復(fù)式淺灘。戴家洲河段池湖港邊灘和巴河邊灘為并聯(lián)方式的邊灘,形態(tài)演變存在一致性對應(yīng)關(guān)系,并與下游戴家洲洲頭低灘組成串聯(lián)方式的邊灘和江心洲;戴家洲左汊內(nèi)無大型洲灘發(fā)育,戴家洲左緣相對穩(wěn)定,暫不分析洲灘之間聯(lián)動(dòng)關(guān)系;戴家洲右緣與樂家灣邊灘形成了汊道內(nèi)部的并聯(lián)式洲灘,一般表現(xiàn)為戴家洲右緣沖刷蝕退,對應(yīng)樂家灣邊灘發(fā)育的聯(lián)動(dòng)關(guān)系。

    池湖港邊灘面積為SCHGBT,巴河邊灘面積為SBHBT,戴家洲洲頭低灘面積為SDJZDT,樂家灣邊灘面積為SLJWBT, 戴家洲右緣灘體面積為SDJZYYZT, 洲灘之間的聯(lián)動(dòng)關(guān)系為: SBHBT— SCHGBT、SDJZDT— SBHBT、SDJZDT—(SCHGBT/SBHBT)、SDJZYYZT—SLJWBT。

    2 戴家洲河段灘槽演變及洲灘聯(lián)動(dòng)關(guān)系

    2.1 河床沖淤量與代表斷面變化

    選取2001 年、2006 年、2011 年、2016 年和2021 年地形,繪制戴家洲河段河床沖淤分布圖;代表斷面選取1984 年、1989 年、1994 年、1999 年、2004 年、2008 年、2013 年、2018 年和2022 年,分析戴家洲河段的床面形態(tài)變化。

    2.1.1 河床沖淤分布

    宜昌站5 000、10 000 和30 000 m3/s 對應(yīng)水面線以下區(qū)域分別為枯水河槽、基本河槽和平灘河槽,其中,枯水至基本河槽之間的區(qū)域?yàn)榈蜑?,基本河槽至平灘河槽之間的區(qū)域?yàn)楦邽?001—2021 年戴家洲河段枯水河槽、基本河槽和平灘河槽沖刷強(qiáng)度分別為15.87 萬、16.48 萬和15.75 萬m3/(km·a), 沖刷集中在枯水河槽,低灘略有淤積及高灘小幅沖刷(圖2( a) )。2001—2006 年枯水河槽沖刷,低灘和高灘區(qū)域?yàn)橛俜e態(tài)勢;2006—2011 年枯水河槽沖刷,低灘沖刷及高灘淤積;2011—2016 年、2016—2021 年枯水河槽沖刷,低灘和高灘均為沖刷態(tài)勢(圖2(a))。2016—2021 年戴家洲河段進(jìn)口王家灣區(qū)域淤積,池湖港邊灘及戴家洲洲頭淤積, 戴家洲右緣顯著沖刷(圖2( b) ) ; 戴家洲右汊灘地為淤積及深槽沖刷態(tài)勢, 左汊整體為淤積態(tài)勢。2001—2011 年戴家洲河段表現(xiàn)為“沖槽淤灘”,逐漸轉(zhuǎn)變?yōu)?011—2011 年灘槽均沖,且枯水河槽沖刷量占平灘河槽沖刷量的92.4%(圖2(a))。

    2.1.2 代表斷面變化

    1984 —2022 年池湖港邊灘區(qū)域斷面的主槽為刷深態(tài)勢, 最大刷深為22.8 m, 池湖港邊灘經(jīng)歷了1984—2013 年淤高,轉(zhuǎn)變?yōu)?013—2022 年刷低且左緣沖刷,最大蝕退距離為140 m(圖3(a))。龍王磯—燕磯之間斷面的戴家洲洲頭低灘在1984—2006 年為沖刷態(tài)勢,2006—2022 年洲頭低灘頭部已由右側(cè)過渡至河道中部, 伴隨洲頭低灘左擺使得右汊進(jìn)口刷深( 圖3( b) ) 。戴家洲汊道進(jìn)口斷面形態(tài)為“ W” 型,1984—2006 年左汊和右汊的斷面沖淤幅度較小,2006—2022 年左汊沖淤幅度小且右汊逐漸刷深,洲頭低灘為淤高態(tài)勢,且2013—2022 年灘體右緣為側(cè)蝕狀態(tài)(圖3(c)、圖3(d))。

    2.2 汊道分流關(guān)系

    1957—2022 年戴家洲汊道分流關(guān)系變化較大(圖4(a)),枯水期出現(xiàn)了汊道交替現(xiàn)象。2006—2022 年戴家洲右汊分流比為增大態(tài)勢, 由枯水期支汊轉(zhuǎn)變?yōu)橹縻?。以相同漢口站流量— 分流比關(guān)系進(jìn)行分析,1965 —1981 年戴家洲右汊分流比高于1957 —1960 年期間,2000 —2008 年為降低態(tài)勢。2010 —2013 年相同流量下戴家洲右汊分流比高于2000—2008 年期間,2014—2018 年為持續(xù)增大態(tài)勢,至2022 年2 月戴家洲右汊分流比( 69.6%,漢口站流量11 800 m3/s)較2000 —2008 年( 49.1%,依據(jù)圖4( b)曲線換算得到)增加了20.5%( 圖4( b) ) 。依據(jù)汊道分流屬性[18], 戴家洲右汊為洪水傾向性汊道[19], 航道工程實(shí)施后的2014 —2022 年洪水傾向性程度較2000—2008 年減弱(圖4(b))。

    2.3 灘槽形態(tài)演變的聯(lián)動(dòng)關(guān)系

    以0 m 等深線為對象,分析1984—2022 年期間代表年份的洲灘形態(tài)變化,統(tǒng)計(jì)了1968—2022 年洲灘面積,分析洲灘演變的關(guān)聯(lián)性與決定性關(guān)系。

    2.3.1 洲灘演變的形態(tài)特征

    戴家洲河段進(jìn)口段右岸為池湖港邊灘,且該邊灘下段右岸分布有龍王磯和燕磯節(jié)點(diǎn)。左岸側(cè)巴河入?yún)R口存在巴河邊灘,在中小水年份邊灘規(guī)模較小,遇到大水年份巴河入?yún)R水量大,疊加大水趨直作用,在干支交匯口附近泥沙淤積使得邊灘逐漸發(fā)育。經(jīng)歷中洪水年份,巴河口至戴家洲左汊進(jìn)口區(qū)域出現(xiàn)一定規(guī)模的沙洲,如1989 年和1994 年(圖5( a) )。經(jīng)歷了1998 年大水后的1999 年巴河邊灘面積顯著增大,戴家洲洲頭較1994 年后退1 384 m(圖5(b))。三峽工程運(yùn)行后,巴河邊灘持續(xù)沖刷,且這一時(shí)期右岸側(cè)池湖港邊灘淤漲淤寬,當(dāng)池湖港邊灘淤漲后的左緣邊界超過燕磯節(jié)點(diǎn)的前端,加速了巴河邊灘沖刷消失,至2011 年邊灘已沖失(圖5( c) ),在戴家洲洲頭上延工程實(shí)施后束窄了河道橫向空間,至2022 年巴河邊灘未出現(xiàn)恢復(fù)淤積跡象。2013 年以來,池湖港邊灘沖刷,引起的河道展寬不利于汊道分流關(guān)系與航道條件穩(wěn)定。一期工程實(shí)施后,2008—2013 年戴家洲頭部淤積上延,2018—2022 年4 期工程實(shí)施后戴家洲頭部位置基本穩(wěn)定,但右緣仍略有沖刷(圖5(d))。樂家灣邊灘與戴家洲右緣的形態(tài)演變存在對應(yīng)關(guān)系,1994—1999 年戴家洲右緣崩退時(shí),對應(yīng)樂家灣邊灘淤積且以沙舌形式進(jìn)入河槽;戴家洲右緣守護(hù)工程實(shí)施后,樂家灣邊灘為沖刷態(tài)勢,如2013 年和2018 年樂家灣邊灘分離為2 部分;樂家灣邊灘守護(hù)工程實(shí)施后,2018—2022 年該邊灘面積增加115.9% 且完整性增強(qiáng)。

    2.3.2 洲灘面積變化

    1968 —2021 年池湖港邊灘面積先增大后減小,其中,2012 —2016 年以左緣崩退形式使得面積減小,2018 年航道工程實(shí)施后的面積變幅小且逐漸穩(wěn)定( 圖6( a) ) 。1968 —2021 年巴河邊灘面積變幅大, 且1998 —1999 年期間的面積達(dá)到最大,至2011 年邊灘基本沖刷消失,至2022 年未出現(xiàn)復(fù)淤跡象(圖6( a) ) 。1968 —2005 年戴家洲面積減小14.2%, 2005 —2012 年由于洲頭守護(hù)工程實(shí)施后面積增加14.5%, 2012 —2018 年減小3.9%,武漢至安慶段6 m 水深航道工程實(shí)施后面積增加2.5%(圖6(b),圖6(d))。1984—2002 年樂家灣邊灘面積增加160.0%, 2003 —2009 年減小45.3%, 2009 —2013 年增加234.2%, 2013 —2016 年減小52.8%,2016 —2018 年變幅較小,2018 —2021 年由于航道工程的實(shí)施面積增加115.9%(圖6( c) )。1984—2007 年戴家洲中下段右緣面積變化較大, 2007 —2009 年右緣平均蝕退140 m, 右緣守護(hù)工程實(shí)施后的2009—2019 年變幅較小,2019—2021 年樂家灣邊灘守護(hù)工程實(shí)施后,對岸側(cè)的戴家洲右緣小幅后退(平均距離45 m) (圖6( f) )。1984—2021 年由于樂家灣邊灘沖刷、戴家洲右緣崩退等泥沙供給的影響,在河道中部形成泥沙淤積體(圖6(e)),至2016 年淤積體的面積最大,2016—2021 年基本沖刷消失,對應(yīng)戴家洲右汊下段航道條件改善。

    在戴家洲汊道分流比調(diào)整關(guān)系上,2006 —2022 年右汊分流比為增大態(tài)勢,由枯水期支汊轉(zhuǎn)變?yōu)橹縻?。受航道整治一期工程影響?006—2013 年期間戴家洲洲頭低灘面積大幅增加,右汊分流比略有增大,對應(yīng)的2006 —2009 年期間右汊內(nèi)洲灘面積減少, 其后戴家洲右緣守護(hù)及右岸側(cè)壩體工程實(shí)施引起2009 —2013 年期間右汊內(nèi)洲體面積為增大態(tài)勢;2013—2018 年期間戴家洲右汊分流比進(jìn)一步增大,戴家洲洲頭低灘、右汊內(nèi)洲體面積均為減小態(tài)勢, 2018—2022 年期間受武漢至安慶段6 m 水深航道整治工程影響,戴家洲洲頭低灘規(guī)模略有減小及右汊內(nèi)洲體面積為增大態(tài)勢,表明航道整治工程起到了穩(wěn)定洲體形態(tài)的作用。

    2.3.3 洲灘演變的聯(lián)動(dòng)性分析

    上游沙洲水道河勢基本穩(wěn)定,出口段深泓穩(wěn)定使得沙洲—戴家洲河段演變關(guān)聯(lián)性弱[20-21],即戴家洲河段內(nèi)的洲灘演變與上游河勢調(diào)整具有相對獨(dú)立性。戴家洲河段內(nèi)邊灘和江心洲存在聯(lián)動(dòng)關(guān)系,聯(lián)動(dòng)演變關(guān)系轉(zhuǎn)換影響著灘群穩(wěn)定性及航道條件的穩(wěn)定性。

    (1)池湖港邊灘與巴河邊灘的關(guān)系。1968—1991 年期間隨著池湖港邊灘面積增加,巴河邊灘面積為減少態(tài)勢,其后的1992—2021 年兩者關(guān)系不顯著(圖7( a) )。1998—2003 年由于長江干流及巴河大水使得巴河邊灘顯著淤積,濾除1998—2003 年數(shù)據(jù)后的池湖港邊灘與巴河邊灘面積表現(xiàn)出較好的相關(guān)性(R2=0.50)。

    ( 2)戴家洲洲頭低灘與池湖港邊灘、巴河邊灘的關(guān)系。池湖港邊灘和巴河邊灘作為戴家洲洲體上游的對峙邊灘, 以池湖港邊灘與巴河邊灘的面積比作為參數(shù)( α) , 分析該參數(shù)變化對戴家洲洲頭低灘的影響。1984 年以前,戴家洲與新淤洲未合并,新淤洲與池湖港邊灘和巴河邊灘的距離相對較大,洲頭低灘變化主要受水沙過程、寡婦磯節(jié)點(diǎn)等綜合影響。隨著參數(shù)α 增大,1984—2003 年和2004—2010 年戴家洲洲頭低灘面積為增大態(tài)勢( R2=0.60~0.92),即池湖港邊灘淤積削弱了燕磯節(jié)點(diǎn)的挑流功能,巴河邊灘快速萎縮(圖7( b) ) 為戴家洲洲頭淤積上延提供了空間。1984 —2010 年巴河邊灘淤漲引起戴家洲洲頭低灘面積減少(R2=0.65~0.82),2004—2010 年巴河邊灘淤漲相同的面積,戴家洲洲頭低灘面積減幅增大(圖7(c)),主要是流域來沙量的大幅減少加速了洲頭沖刷態(tài)勢。

    (3)樂家灣邊灘與戴家洲右緣的關(guān)系。1984—1999 年期間,隨著樂家灣邊灘面積增加,戴家洲右緣為蝕退態(tài)勢( R2=0.18,圖7( d) ),這一時(shí)期戴家洲河段岸線及洲體邊界線尚未穩(wěn)定,局部崩岸及洲灘坍塌使得樂家灣邊灘與戴家洲右緣面積變化表現(xiàn)為弱相關(guān)。2000—2008 年期間,隨著樂家灣邊灘面積增加,戴家洲右緣為蝕退態(tài)勢( R2=0.72,圖7( d) ),同時(shí)戴家洲右緣蝕退也為樂家灣邊灘發(fā)育提供了橫向空間。2008 年以來,戴家洲河段已實(shí)施了4 期航道工程,受戴家洲右緣守護(hù)工程、右汊分流比增加等綜合影響,先期淤漲的樂家灣邊灘逐漸沖刷并形成淤積體進(jìn)入河槽,對應(yīng)的2009—2019 年戴家洲右緣相對穩(wěn)定。武漢至安慶段航道工程實(shí)施后,樂家灣邊灘形態(tài)得到穩(wěn)定,對應(yīng)的2019—2021 年戴家洲右緣仍出現(xiàn)了小幅沖刷,建議增強(qiáng)戴家洲右緣守護(hù)工程的穩(wěn)定性。

    3 航道工程對戴家洲河段灘槽形態(tài)的影響機(jī)制

    3.1 水沙條件及航道工程對灘群面積的影響

    1968—2022 年漢口站徑流量變化不大,2003—2022 年較1968—2002 年、1984—2002 年輸沙量減幅分別為75.7% 和73.0%(圖8)。1968—2002 年戴家洲洲體面積與漢口站沙量關(guān)系不顯著,2003—2007 年洲灘面積減少與漢口站輸沙量減少表現(xiàn)為同步性;2008 年以來,雖然漢口站輸沙量處于較低水平,但戴家洲洲體面積未出現(xiàn)趨勢性減少,主要受航道工程影響使得2008—2013 年洲體面積為增大態(tài)勢,武漢至安慶段工程實(shí)施后洲體面積逐漸穩(wěn)定。1968 —2005 年戴家洲河段洲體總面積隨著漢口站輸沙量減少表現(xiàn)為同步性減小,2005—2013 年戴家洲河段洲灘面積為增大態(tài)勢,武漢至安慶段工程實(shí)施后洲灘面積基本穩(wěn)定。

    3.2 航道工程對汊道沖淤關(guān)系的影響

    三峽工程運(yùn)行以來,長江中下游分汊河段內(nèi)水流流程短的汊道沖刷發(fā)展[22],右汊河道長度短且河寬大有利于船舶航行,綜合選取了右汊作為枯水期主航道[6]。2006—2021 年戴家洲左汊和右汊枯水河槽沖淤量分別為281.70 萬和 ? 4 126.40 萬m3,平灘河槽沖淤量分別為1 077.85 萬和 ? 4 376.40 萬m3,整體為左汊淤積和右汊沖刷的特征(表1) 。2006 —2021 年戴家洲左汊淤積以低灘為主( 比例為73.9%) , 右汊沖刷且以枯水河槽為主(比例為94.3%)。2006 —2011 年、2011 —2016 年和2016 —2021 年戴家洲右汊枯水河槽沖刷強(qiáng)度分別為82.88 萬、300.20 萬和442.20 萬 m3/a, 平灘河槽沖刷強(qiáng)度分別為111.48 萬、224.40 萬和539.40 萬 m3/a,右汊沖刷強(qiáng)度為增大態(tài)勢。2011—2016 年期間漢口站沙量較2006—2011 年期間減幅為10.4%、枯水期分流比最小值增幅為3.8%(由48.2% 增加至52.0%),而右汊平灘河槽沖刷強(qiáng)度增大101.3%;2016—2021 年期間漢口站沙量較2011—2016 年期間減幅為16.3%、枯水期分流比最小值增幅為8.8%(由52.0% 增加至60.8%),而右汊平灘河槽沖刷強(qiáng)度增大140.4%。2006—2011 年、2011—2016 年和2016—2021 年期間戴家洲右汊航道疏浚維護(hù)量分別為3.0 萬、94.13 萬和7.2 萬 m3,分別占枯水河槽沖刷量的0.72%、6.27% 和0.32%,對戴家洲汊道沖淤強(qiáng)度影響相對較小。綜合分析認(rèn)為,在流域來沙量減少背景下,戴家洲右汊處于沖刷態(tài)勢,航道工程實(shí)施穩(wěn)定了灘槽形態(tài)且顯著增加了右汊分流比, 是引起右汊沖刷強(qiáng)度2016 —2021 年高于2011 —2016 年、2011—2016 年高于2006—2011 年的重要因素。

    3.3 航道工程對水深條件的改善作用

    3.3.1 航道條件變化過程

    戴家洲河段4.0 m 和4.5 m 維護(hù)水深時(shí)期的礙航位置集中在戴家洲右汊進(jìn)口(龍王磯—耀子山)、中段(耀子山—寡婦磯)及出口(寡婦磯—?jiǎng)⒓叶桑﹨^(qū)域,在一期、二期和三期航道工程實(shí)施后,4.5 m 水深航道條件顯著改善[5-7]。2004 年和2008 年測圖顯示,戴家洲河段左汊進(jìn)口6 m 槽斷開,且右汊進(jìn)口、中段及出口段也為斷開態(tài)勢,但右汊斷開距離為減少態(tài)勢;2013 年測圖顯示,戴家洲河段左汊進(jìn)口6 m 槽斷開,且右汊進(jìn)口和中段斷開,下段為貫通態(tài)勢(圖9(a));2018 年水深測圖顯示,戴家洲河段左汊和右汊進(jìn)口段6m 槽斷開,右汊出口段6 m 槽寬度為63 m。2018 年起戴家洲河段實(shí)施了6 m 水深航道工程,2022 年左汊進(jìn)口6 m 槽斷開,右汊6 m 槽通過維護(hù)性疏浚保持了貫通,出口段6 m 槽寬度為484 m(圖9( b) )。在三峽工程運(yùn)行后,戴家洲河段枯水期主航道的右汊航道水深條件轉(zhuǎn)好,2004 —2018 年右汊礙航位置為進(jìn)口、中段及尾部共3 個(gè)區(qū)域,2018 年以來的礙航位置主要為右汊進(jìn)口。

    3.3.2 礙航程度變化

    隨著流域來沙量減少,戴家洲河段河床為沖刷態(tài)勢,為航道水深增加提供了有利基礎(chǔ)。在航道工程實(shí)施前( 1984—2008 年),結(jié)合疏浚維護(hù)實(shí)現(xiàn)了枯水期航道4.0 m 水深暢通(圖10);一期、二期和三期工程實(shí)施后,戴家洲左汊進(jìn)口、右汊進(jìn)口和中下段轉(zhuǎn)為貫通或?qū)挾炔蛔?,水深不足的比例大幅降低(圖10);戴家洲河段四期航道工程的目標(biāo)之一是提升右汊枯水期主航道地位,工程實(shí)施后右汊中下段礙航程度降低,進(jìn)口段通過維護(hù)疏浚措施后航道尺度已達(dá)到6 m×110 m。

    4 結(jié) 論

    本研究以長江中游戴家洲河段為對象,采用了1954—2022 年水沙和地形數(shù)據(jù),解析了灘槽形態(tài)演變過程與聯(lián)動(dòng)關(guān)系,識(shí)別了航道工程對洲灘演變的影響,主要結(jié)論如下:

    ( 1)三峽工程運(yùn)行以來,戴家洲河段枯水河槽和平灘河槽均為沖刷態(tài)勢,其中,低灘區(qū)域表現(xiàn)為沖刷,高灘區(qū)域表現(xiàn)為淤積態(tài)勢,汊道間表現(xiàn)為右汊沖刷和左汊淤積的關(guān)系,其中2011 年以來航道工程實(shí)施后的右汊沖刷強(qiáng)度增加,分流比為增大態(tài)勢。

    (2)在自然水沙階段(1968—2008 年),戴家洲河段進(jìn)口段巴河邊灘與池湖港邊灘具有強(qiáng)聯(lián)動(dòng)關(guān)系,池湖港邊灘持續(xù)淤漲且逐漸降低了燕磯節(jié)點(diǎn)的挑流作用,利于戴家洲洲頭低灘淤積上延,且巴河邊灘為本時(shí)期戴家洲汊道進(jìn)口段灘群演變的主導(dǎo)因素;巴河邊灘的淤漲引起戴家洲洲頭低灘沖刷,但流域來沙量減少加劇了2004—2010 年戴家洲洲頭低灘的沖刷萎縮;右汊內(nèi)戴家洲右緣蝕退為樂家灣邊灘發(fā)育提供空間,枯水期主航道的礙航淺區(qū)分布在汊道進(jìn)口、右汊中段及出口段。

    (3)航道工程影響階段(2008—2022 年),戴家洲河段洲灘穩(wěn)定性增強(qiáng)且洲灘演變關(guān)聯(lián)性減弱,右汊枯水期分流比增加,池湖港邊灘規(guī)模減少增加了進(jìn)口段水流擺動(dòng)空間,4.5 m 和6 m 水深不足區(qū)域集中在右汊進(jìn)口段。

    參考文獻(xiàn):

    [1]趙維陽, 楊云平, 張華慶, 等. 三峽大壩下游近壩段沙質(zhì)河床形態(tài)調(diào)整及洲灘聯(lián)動(dòng)演變關(guān)系[J]. 水科學(xué)進(jìn)展, 2020,31( 6):862-874. (ZHAO W Y,YANG Y P,ZHANG H Q,et al. Adjustment patterns and causes of the morphology of sandy riverbed"downstream of the Three Gorges Dam[J]. Advances in Water Science,2020,31(6):862-874. (in Chinese))

    [2]劉亞, 鄭力, 姚仕明, 等. 分汊河道主支汊交替主導(dǎo)因子的轉(zhuǎn)換模擬[J]. 水科學(xué)進(jìn)展, 2020, 31( 3) : 348-355. (LIU Y,ZHENG L,YAO S M,et al. Simulation of dominant factors transforming on major-minor branches alternation in anabanched rivers[J].Advances in Water Science,2020,31(3):348-355. (in Chinese))

    [3]楊云平,張夏博,鄭金海,等. 長江下游落成洲河段洲灘聯(lián)動(dòng)關(guān)系與航道淺灘礙航機(jī)制[J]. 水科學(xué)進(jìn)展,2023,34(2):250-264. (YANG Y P,ZHANG X B,ZHENG J H,et al. Linkage relationship of beach/central bar and waterway shoal obstruction mechanismin Luochengzhou reach of the lower reaches of Yangtze River[J]. Advances in Water Science,2023,34(2):250-264. (in Chinese))

    [4]楊云平, 鄭金海, 張明進(jìn), 等. 長江下游潮流界變動(dòng)段三益橋邊灘與淺灘演變驅(qū)動(dòng)機(jī)制分析[J]. 水科學(xué)進(jìn)展, 2020,31( 4) : 502-513. ( YANG Y P,ZHENG J H,ZHANG M J,et al. Driving mechanism of Sanyiqiao point bar and shoal evolution influctuation segment of tidal current limit in lower reaches of Yangtze River[J]. Advances in Water Science,2020,31(4):502-513. (inChinese))

    [5]楊云平, 李明, 劉萬利, 等. 長江荊江河段灘槽演變與航道水深資源提升關(guān)系[J]. 水科學(xué)進(jìn)展, 2022, 33( 2) : 240-252.(YANG Y P,LI M,LIU W L,et al. Study on the relationship between beach trough evolution and navigation obstruction characteristics inJingjiang reach of the Yangtze River[J]. Advances in Water Science,2022,33(2):240-252. (in Chinese))

    [6]YANG Y P, ZHENG J H, ZHANG M J, et al. Impacts of human activities on the riverbed morphological in the tidal reaches of theYangtze River[J]. Journal of Hydrology,2024,630:130735.

    [7]劉亞,姚仕明,謝炎,等. 彎曲分汊河型主支汊交替模式及驅(qū)動(dòng)機(jī)制研究[J]. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào),2020,28(1):18-26. ( LIU Y, YAO S M, XIE Y, et al. Mode and domination of major-minor alternation for anabranching bends[J]. Journal of BasicScience and Engineering,2020,28(1):18-26. (in Chinese))

    [8]YANG Y P,ZHENG J H,ZHANG W,et al. Quantitative relationship between channels and bars in a tidal reach of the Lower YangtzeRiver:implications for river management[J]. Journal of Geographical Sciences,2021,31(12):1837-1851.

    [9]WANG J J, YANG Y P, ZHANG M J, et al. Critical threshold of periodic point bar scour and sediment body transport path in tidalreaches:a case study of Fujiangsha reach,Yangtze River[J]. Journal of Geographical Sciences,2023,33(6):1334-1358.

    [10]劉萬利,李旺生,李一兵,等. 長江中游典型微彎分汊河段擇汊問題探討[J]. 水運(yùn)工程,2012(8):141-144. (LIU W L,LI WS, LI Y B, et al. Navigable channel selection of typical braided reach in Middle Yangtze River[J]. Port amp; Waterway Engineering,2012(8):141-144. (in Chinese))

    [11]劉萬利,李旺生,朱玉德,等. 長江中游戴家洲河段航道整治思路探討[J]. 水道港口,2009,30(1):31-36. (LIU W L,LI WS,ZHU Y D,et al. Study on channel regulation of Daijiazhou reach in Middle Yangtze River[J]. Journal of Waterway and Harbor,2009,30(1):31-36. (in Chinese))

    [12]張明進(jìn), 張華慶, 劉萬利, 等. 長江中游戴家洲河段航道整治工程數(shù)值模擬研究[J]. 水道港口, 2008, 29( 3) : 193-198.(ZHANG M J,ZHANG H Q,LIU W L,et al. Numerical study on navigation regulation engineering of Daijiazhou section in the middlereaches of the Yangtze River[J]. Journal of Waterway and Harbor,2008,29(3):193-198. (in Chinese))

    [13]蔡大富. 長江中游戴家洲河段航道整治一期工程效果分析[J]. 水運(yùn)工程, 2011( 2) : 91-96. (CAI D F. Analysis of effects ofDaijiazhou section channel regulation phaseⅠproject in the middle reach of the Yangtze River[J]. Port amp; Waterway Engineering,2011(2):91-96. (in Chinese))

    [14]黃蓓蓓,蔡大富. 長江中游戴家洲河段戴家洲右緣下段守護(hù)工程效果分析[J]. 水運(yùn)工程,2014( 12):110-115. (HUANG B B,CAI D F. Effect analysis of Daijiazhou right border lower section protection engineering of Daijiazhou channel segment in the middle reach ofthe Yangtze River[J]. Port amp; Waterway Engineering,2014(12):110-115. (in Chinese))

    [15]余文鈞,李長鈴. 長江中游戴家洲水道6.0m 航道整治工程效果分析[J]. 水運(yùn)工程,2023(6):144-148,164. (YU W J,LI C L.Effect analysis of 6.0m waterway regulation project of Daijiazhou channel in middle reaches of the Yangtze River[J]. Port amp; WaterwayEngineering,2023(6):144-148,164. (in Chinese))

    [16]孫昭華,馮秋芬,韓劍橋,等. 順直河型與分汊河型交界段洲灘演變及其對航道條件影響:以長江天興洲河段為例[J]. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào),2013,21(4):647-656. (SUN Z H,F(xiàn)ENG Q F,HAN J Q,et al. Fluvial processes of sandbars in the junctionreach of single-threaded channel to anabranching channel and its impact on navigation: a case study of the Tianxingzhou reach of theYangtze River[J]. Journal of Basic Science and Engineering,2013,21(4):647-656. (in Chinese))

    [17]劉林,黃成濤,李明,等. 長江中游典型順直河段交錯(cuò)邊灘復(fù)歸性演變機(jī)理[J]. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào),2014,22( 3):445-456. ( LIU L,HUANG C T,LI M,et al. Periodic evolution mechanism of staggered beach in typical straight reach of the MiddleYangtze River[J]. Journal of Basic Science and Engineering,2014,22(3):445-456. (in Chinese))

    [18]韓劍橋,張為,袁晶,等. 三峽水庫下游分汊河道灘槽調(diào)整及其對水文過程的響應(yīng)[J]. 水科學(xué)進(jìn)展,2018,29(2):186-195.( HAN J Q,ZHANG W,YUAN J,et al. Responses of riverbed morphology to the hydrological regime in anabranching reaches of thedownstream Three Gorges Reservoir[J]. Advances in Water Science,2018,29(2):186-195. (in Chinese))

    [19]劉萬利,李旺生,朱玉德,等. 長江中游戴家洲河段天然水流特性試驗(yàn)研究[J]. 水[19] 道港口,2009,30(2):113-118. (LIU WL,LI W S,ZHU Y D,et al. Study on flow characteristics of braided reach of Daijiazhou reach in Middle Yangtze River[J]. Journal ofWaterway and Harbor,2009,30(2):113-118. (in Chinese))

    [20]蔡大富,李青云. 長江中游戴家洲河段航道整治工程方案研究[J]. 水運(yùn)工程,2008( 12):105-109,113. (CAI D F,LI Q Y.Channel regulation scheme of the Daijiazhou section in the middle reach of the Yangtze River[J]. Port amp; Waterway Engineering,2008(12):105-109,113. (in Chinese))

    [21]蔡大富. 長江中游戴家洲河段河床演變及礙航特性分析[J]. 泥沙研究,2011,36( 2):47-54. (CAI D F. River channel evolutionand charichateristics of navigation-obstruction in Daijiazhou section in Middle Yangtze River[J]. Journal of Sediment Research, 2011,36(2):47-54. (in Chinese))

    [22]方娟娟,付菁,楊云平,等. 長江下游儀征河段汊道分流屬性及航道灘槽演變機(jī)制[J]. 湖泊科學(xué),2022,34(4):1335-1346.(FANG J J,F(xiàn)U J,YANG Y P,et al. Diversion attribute of the branch channels and shoal-channel evolution mechanism in the Yizhengsection of the lower reaches of the Yangtze River[J]. Journal of Lake Sciences,2022,34(4):1335-1346. (in Chinese))

    18禁动态无遮挡网站| 晚上一个人看的免费电影| 精品亚洲乱码少妇综合久久| 久久久久免费精品人妻一区二区| av在线亚洲专区| 两个人视频免费观看高清| av一本久久久久| 日韩大片免费观看网站| 又大又黄又爽视频免费| 99久久中文字幕三级久久日本| 亚洲av免费在线观看| 我要看日韩黄色一级片| 久久久久久国产a免费观看| 淫秽高清视频在线观看| 99九九线精品视频在线观看视频| 激情五月婷婷亚洲| 69av精品久久久久久| 国产亚洲av嫩草精品影院| 亚洲久久久久久中文字幕| 色5月婷婷丁香| 亚洲精品日韩在线中文字幕| 国产精品女同一区二区软件| 人妻系列 视频| 天堂av国产一区二区熟女人妻| 在线天堂最新版资源| 久久久久精品久久久久真实原创| 成人无遮挡网站| 内地一区二区视频在线| 免费大片黄手机在线观看| 国产一区二区三区av在线| 国产视频首页在线观看| 在线免费十八禁| 天堂网av新在线| 你懂的网址亚洲精品在线观看| 精品一区二区免费观看| 久久国内精品自在自线图片| 国产成人a区在线观看| 一级毛片黄色毛片免费观看视频| 国产国拍精品亚洲av在线观看| 久久国内精品自在自线图片| 午夜免费观看性视频| 一区二区三区高清视频在线| 91久久精品国产一区二区三区| 精品人妻熟女av久视频| 亚洲三级黄色毛片| 日韩精品青青久久久久久| 精品久久国产蜜桃| 狂野欧美白嫩少妇大欣赏| 午夜免费观看性视频| 91久久精品电影网| 久久综合国产亚洲精品| 久久久久久久国产电影| 免费看光身美女| 亚洲精品日韩av片在线观看| 18禁在线无遮挡免费观看视频| 精品99又大又爽又粗少妇毛片| 久久久久国产网址| 免费观看a级毛片全部| 欧美一区二区亚洲| 美女高潮的动态| kizo精华| 亚洲精品日韩av片在线观看| 最后的刺客免费高清国语| 又爽又黄a免费视频| 亚洲成人av在线免费| 草草在线视频免费看| 国产成人a∨麻豆精品| 日本三级黄在线观看| 国产高清不卡午夜福利| 亚洲熟妇中文字幕五十中出| 秋霞在线观看毛片| 黄色一级大片看看| 国产精品国产三级国产av玫瑰| 亚洲综合精品二区| 精品久久久久久久人妻蜜臀av| 日韩av不卡免费在线播放| 亚洲精品国产成人久久av| 99热6这里只有精品| 国产午夜精品一二区理论片| 亚洲精品日韩av片在线观看| 亚洲精品国产av成人精品| 国产亚洲av片在线观看秒播厂 | 久久久久久伊人网av| 婷婷色综合www| 国产成人精品一,二区| 日本黄大片高清| 少妇的逼好多水| 乱系列少妇在线播放| 精品久久久久久久人妻蜜臀av| 午夜福利在线在线| 午夜亚洲福利在线播放| 晚上一个人看的免费电影| 成年版毛片免费区| 国产精品一及| 久99久视频精品免费| 九色成人免费人妻av| av在线观看视频网站免费| 国产一区二区三区综合在线观看 | 国产精品无大码| 国产精品国产三级国产av玫瑰| 亚洲最大成人av| 国产高清有码在线观看视频| 成人性生交大片免费视频hd| 免费黄网站久久成人精品| 精品久久久久久久人妻蜜臀av| 国产精品.久久久| 91午夜精品亚洲一区二区三区| 热99在线观看视频| 国产黄色小视频在线观看| 欧美 日韩 精品 国产| 国产成人精品福利久久| 日日撸夜夜添| 男女视频在线观看网站免费| 国产在视频线在精品| 国产有黄有色有爽视频| 噜噜噜噜噜久久久久久91| 青春草国产在线视频| 国产精品99久久久久久久久| 一级毛片 在线播放| 男人和女人高潮做爰伦理| av专区在线播放| 汤姆久久久久久久影院中文字幕 | 美女被艹到高潮喷水动态| 国内少妇人妻偷人精品xxx网站| 国产精品日韩av在线免费观看| 国产精品一二三区在线看| 国产一区二区三区综合在线观看 | 精华霜和精华液先用哪个| 好男人视频免费观看在线| 久久久久久久亚洲中文字幕| 亚洲激情五月婷婷啪啪| 亚洲国产精品成人综合色| 99热这里只有是精品在线观看| 亚洲三级黄色毛片| 22中文网久久字幕| 精品午夜福利在线看| 国产精品嫩草影院av在线观看| 日韩欧美一区视频在线观看 | 高清在线视频一区二区三区| 神马国产精品三级电影在线观看| 69av精品久久久久久| 丝袜喷水一区| 亚洲熟女精品中文字幕| 日本欧美国产在线视频| 日日摸夜夜添夜夜添av毛片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲一区高清亚洲精品| 日本免费在线观看一区| 一级毛片电影观看| 亚洲精品乱码久久久v下载方式| 一级毛片黄色毛片免费观看视频| 青青草视频在线视频观看| 欧美成人午夜免费资源| 精品人妻偷拍中文字幕| 国产精品久久视频播放| 少妇的逼水好多| 亚洲熟女精品中文字幕| 午夜爱爱视频在线播放| 99久久九九国产精品国产免费| 免费电影在线观看免费观看| 日本黄色片子视频| 国产av码专区亚洲av| 日本猛色少妇xxxxx猛交久久| 麻豆av噜噜一区二区三区| 99久久九九国产精品国产免费| 国产成年人精品一区二区| 国产成人aa在线观看| 熟妇人妻久久中文字幕3abv| 亚洲熟妇中文字幕五十中出| 久久草成人影院| 你懂的网址亚洲精品在线观看| a级毛片免费高清观看在线播放| 91精品国产九色| 高清午夜精品一区二区三区| 一个人免费在线观看电影| 亚洲国产精品国产精品| 午夜日本视频在线| videossex国产| 精品一区二区三卡| 国内揄拍国产精品人妻在线| 能在线免费观看的黄片| 夫妻性生交免费视频一级片| 舔av片在线| 欧美高清性xxxxhd video| 九九在线视频观看精品| 听说在线观看完整版免费高清| 女的被弄到高潮叫床怎么办| 亚洲电影在线观看av| 麻豆国产97在线/欧美| 色视频www国产| 黄色一级大片看看| 亚洲最大成人av| 内地一区二区视频在线| 久久99精品国语久久久| 亚洲av免费高清在线观看| 久久亚洲国产成人精品v| 99久久精品热视频| 久久99热6这里只有精品| 岛国毛片在线播放| 国国产精品蜜臀av免费| 精品久久久久久久久久久久久| 国产爱豆传媒在线观看| 最后的刺客免费高清国语| 亚洲人成网站在线观看播放| 亚洲精品成人久久久久久| 99久久精品国产国产毛片| 亚洲无线观看免费| 国产av不卡久久| 久热久热在线精品观看| 天堂中文最新版在线下载 | 欧美成人a在线观看| 日韩一区二区视频免费看| 亚洲精品日韩av片在线观看| 亚洲人与动物交配视频| 亚洲电影在线观看av| 免费大片18禁| 日韩一区二区视频免费看| 少妇丰满av| 久久久久精品性色| 日本熟妇午夜| 综合色av麻豆| 国产亚洲av片在线观看秒播厂 | 久久久久久久午夜电影| 亚洲人与动物交配视频| 小蜜桃在线观看免费完整版高清| 亚洲精品亚洲一区二区| 国产国拍精品亚洲av在线观看| 国产黄色免费在线视频| 国产黄a三级三级三级人| 能在线免费观看的黄片| 久久亚洲国产成人精品v| 免费看a级黄色片| 久久6这里有精品| 男人狂女人下面高潮的视频| 国产淫语在线视频| 超碰97精品在线观看| 久久久久久久久久久免费av| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲,欧美,日韩| 中文欧美无线码| 欧美激情在线99| 久久精品国产鲁丝片午夜精品| h日本视频在线播放| 国产三级在线视频| 中文精品一卡2卡3卡4更新| 一区二区三区乱码不卡18| 中文在线观看免费www的网站| 欧美最新免费一区二区三区| 男女下面进入的视频免费午夜| 少妇猛男粗大的猛烈进出视频 | 熟妇人妻不卡中文字幕| 九草在线视频观看| 亚洲欧洲国产日韩| 欧美成人一区二区免费高清观看| 国产亚洲av片在线观看秒播厂 | 男女边吃奶边做爰视频| 97精品久久久久久久久久精品| 成人午夜精彩视频在线观看| 午夜福利在线在线| 欧美激情在线99| 午夜免费男女啪啪视频观看| 九色成人免费人妻av| 久久精品国产鲁丝片午夜精品| 国产午夜精品久久久久久一区二区三区| 毛片女人毛片| 少妇高潮的动态图| 免费观看性生交大片5| 80岁老熟妇乱子伦牲交| 一级片'在线观看视频| 亚洲精品乱码久久久v下载方式| av在线观看视频网站免费| 五月天丁香电影| 波野结衣二区三区在线| 女人十人毛片免费观看3o分钟| 成人欧美大片| 中文精品一卡2卡3卡4更新| 高清毛片免费看| 欧美日韩在线观看h| 熟女人妻精品中文字幕| 边亲边吃奶的免费视频| 久久久色成人| 亚洲av二区三区四区| 深爱激情五月婷婷| 国产又色又爽无遮挡免| 国产一区二区亚洲精品在线观看| 久久精品国产亚洲av天美| 亚洲自拍偷在线| 简卡轻食公司| 亚洲av日韩在线播放| videossex国产| 亚洲欧美日韩东京热| 久久久精品欧美日韩精品| 亚洲三级黄色毛片| 亚洲av成人av| 大又大粗又爽又黄少妇毛片口| 亚洲伊人久久精品综合| 熟妇人妻不卡中文字幕| 少妇丰满av| 欧美3d第一页| 街头女战士在线观看网站| 国产黄色视频一区二区在线观看| 纵有疾风起免费观看全集完整版 | 人妻制服诱惑在线中文字幕| 嫩草影院新地址| 免费av观看视频| 亚洲av成人av| 精品国内亚洲2022精品成人| 国产亚洲91精品色在线| 联通29元200g的流量卡| 日韩av免费高清视频| 老师上课跳d突然被开到最大视频| 国产不卡一卡二| 啦啦啦啦在线视频资源| 国产真实伦视频高清在线观看| 午夜激情福利司机影院| 国产精品综合久久久久久久免费| av女优亚洲男人天堂| 天堂中文最新版在线下载 | 最近2019中文字幕mv第一页| ponron亚洲| 国产亚洲精品久久久com| 成人亚洲精品一区在线观看 | 日韩不卡一区二区三区视频在线| 免费黄网站久久成人精品| av免费观看日本| 日本黄色片子视频| 国产午夜精品久久久久久一区二区三区| 免费黄色在线免费观看| 男女边摸边吃奶| 麻豆成人午夜福利视频| 成人毛片60女人毛片免费| 2021少妇久久久久久久久久久| 九九在线视频观看精品| 免费观看的影片在线观看| 国产老妇女一区| 毛片一级片免费看久久久久| av福利片在线观看| 性色avwww在线观看| 欧美最新免费一区二区三区| 日韩欧美精品免费久久| 国产有黄有色有爽视频| 国内精品宾馆在线| 亚洲欧美一区二区三区黑人 | 国产精品无大码| 国产视频内射| 亚洲成人精品中文字幕电影| 男人狂女人下面高潮的视频| 中文欧美无线码| 色视频www国产| 免费看a级黄色片| 人妻系列 视频| freevideosex欧美| 欧美丝袜亚洲另类| 街头女战士在线观看网站| 麻豆精品久久久久久蜜桃| 久久久午夜欧美精品| 一级黄片播放器| 男女啪啪激烈高潮av片| 国产乱来视频区| 亚洲精品自拍成人| 午夜激情福利司机影院| 久久99蜜桃精品久久| 日本黄色片子视频| 日韩中字成人| 亚洲不卡免费看| 亚洲av成人精品一二三区| 免费黄网站久久成人精品| 最新中文字幕久久久久| 国产在线一区二区三区精| 亚洲美女视频黄频| 日日啪夜夜爽| 亚洲18禁久久av| 亚洲人成网站高清观看| 国产黄a三级三级三级人| 一级毛片aaaaaa免费看小| 国产国拍精品亚洲av在线观看| 青春草视频在线免费观看| 2018国产大陆天天弄谢| 日韩欧美一区视频在线观看 | 岛国毛片在线播放| 欧美xxxx黑人xx丫x性爽| 欧美日韩综合久久久久久| 一区二区三区免费毛片| 国产熟女欧美一区二区| 国产片特级美女逼逼视频| 欧美精品一区二区大全| 只有这里有精品99| 久久草成人影院| 少妇被粗大猛烈的视频| 国产乱人偷精品视频| 日本wwww免费看| 亚洲图色成人| 五月玫瑰六月丁香| 精品久久久久久久人妻蜜臀av| 国产亚洲最大av| 久久精品熟女亚洲av麻豆精品 | 日本一本二区三区精品| 精品久久久久久久久av| 久久亚洲国产成人精品v| 亚洲av一区综合| 中国美白少妇内射xxxbb| 青青草视频在线视频观看| 国产探花极品一区二区| 少妇猛男粗大的猛烈进出视频 | 日韩大片免费观看网站| 久久精品熟女亚洲av麻豆精品 | 亚洲在线自拍视频| 一级毛片电影观看| 一个人看视频在线观看www免费| 欧美zozozo另类| 三级经典国产精品| 午夜福利视频1000在线观看| 国产女主播在线喷水免费视频网站 | 亚洲一区高清亚洲精品| 亚洲精品乱码久久久久久按摩| 亚洲四区av| 丝袜喷水一区| 成人亚洲精品一区在线观看 | 欧美激情在线99| 小蜜桃在线观看免费完整版高清| 美女黄网站色视频| 国产亚洲一区二区精品| 国产熟女欧美一区二区| 久久精品久久精品一区二区三区| 成人美女网站在线观看视频| 91午夜精品亚洲一区二区三区| 精品久久久久久久末码| 亚洲丝袜综合中文字幕| 国产成人91sexporn| 特大巨黑吊av在线直播| 日日啪夜夜撸| 成年人午夜在线观看视频 | 亚洲精品乱久久久久久| 中文字幕制服av| 亚洲精品,欧美精品| 久99久视频精品免费| 夜夜看夜夜爽夜夜摸| 天天躁日日操中文字幕| 亚洲,欧美,日韩| 日韩强制内射视频| 久久精品熟女亚洲av麻豆精品 | 五月天丁香电影| av福利片在线观看| 亚洲精品久久午夜乱码| 啦啦啦韩国在线观看视频| 菩萨蛮人人尽说江南好唐韦庄| 青春草国产在线视频| 国产乱人视频| 午夜激情欧美在线| 国产乱人视频| 一二三四中文在线观看免费高清| 成人欧美大片| 亚洲内射少妇av| 亚洲最大成人av| 一级毛片aaaaaa免费看小| a级毛片免费高清观看在线播放| 99re6热这里在线精品视频| 丝袜美腿在线中文| av.在线天堂| 国产精品嫩草影院av在线观看| 永久免费av网站大全| 最近手机中文字幕大全| 熟女电影av网| 亚洲欧美清纯卡通| 中文字幕av成人在线电影| 久久99热6这里只有精品| 成人国产麻豆网| 国产高清三级在线| 青春草亚洲视频在线观看| 黄色配什么色好看| 一本久久精品| 久久久午夜欧美精品| 中文字幕制服av| 日韩欧美三级三区| 久久这里有精品视频免费| 亚洲最大成人av| 99久久九九国产精品国产免费| 日本一二三区视频观看| 亚洲精品久久久久久婷婷小说| 成人亚洲精品av一区二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 波多野结衣巨乳人妻| 搡老乐熟女国产| 国产成人精品福利久久| 日韩亚洲欧美综合| 深夜a级毛片| 久久久久久久久久久丰满| 亚洲内射少妇av| 十八禁网站网址无遮挡 | 免费高清在线观看视频在线观看| 丝瓜视频免费看黄片| 亚洲国产欧美在线一区| 熟妇人妻久久中文字幕3abv| 禁无遮挡网站| 蜜臀久久99精品久久宅男| 久久精品国产亚洲av天美| 日韩av免费高清视频| 搞女人的毛片| 日本欧美国产在线视频| 国产真实伦视频高清在线观看| 亚洲天堂国产精品一区在线| 嫩草影院精品99| 女人久久www免费人成看片| 亚洲精品日本国产第一区| 91aial.com中文字幕在线观看| 亚洲天堂国产精品一区在线| 舔av片在线| h日本视频在线播放| 人妻少妇偷人精品九色| 26uuu在线亚洲综合色| 国产色婷婷99| 日韩成人av中文字幕在线观看| 在线a可以看的网站| 少妇人妻精品综合一区二区| 久久99热这里只有精品18| 午夜免费观看性视频| 欧美激情在线99| 欧美日韩综合久久久久久| 在线免费观看不下载黄p国产| 精品国产三级普通话版| 午夜精品一区二区三区免费看| 综合色丁香网| 免费观看av网站的网址| 亚洲成人一二三区av| 免费不卡的大黄色大毛片视频在线观看 | 美女国产视频在线观看| 午夜免费男女啪啪视频观看| 日本爱情动作片www.在线观看| 在线天堂最新版资源| 国产视频内射| 爱豆传媒免费全集在线观看| 伦精品一区二区三区| 高清欧美精品videossex| 亚洲欧美精品专区久久| 亚洲av成人精品一区久久| 视频中文字幕在线观看| 高清午夜精品一区二区三区| 美女高潮的动态| 亚洲av免费在线观看| 特级一级黄色大片| 免费观看的影片在线观看| 日韩欧美精品v在线| 久久久a久久爽久久v久久| 噜噜噜噜噜久久久久久91| 大陆偷拍与自拍| 国产视频内射| 亚洲av免费在线观看| 精品少妇黑人巨大在线播放| 国产老妇伦熟女老妇高清| 欧美bdsm另类| 赤兔流量卡办理| videossex国产| 国产av在哪里看| 国产亚洲av片在线观看秒播厂 | 欧美变态另类bdsm刘玥| 精品久久久久久久末码| av在线观看视频网站免费| 丝袜喷水一区| 啦啦啦中文免费视频观看日本| 三级国产精品欧美在线观看| 亚洲国产成人一精品久久久| 久久久久久久午夜电影| 精品不卡国产一区二区三区| 天堂俺去俺来也www色官网 | 丝瓜视频免费看黄片| 国产精品爽爽va在线观看网站| 白带黄色成豆腐渣| 亚洲国产精品成人久久小说| 午夜福利网站1000一区二区三区| 久久精品夜色国产| 一个人看视频在线观看www免费| 亚洲国产欧美人成| 国产精品不卡视频一区二区| 国产综合精华液| 看十八女毛片水多多多| 午夜福利在线观看吧| 国产av码专区亚洲av| 国产人妻一区二区三区在| 亚洲精品乱码久久久久久按摩| 久久久久性生活片| 91aial.com中文字幕在线观看| 99九九线精品视频在线观看视频| 国产av国产精品国产| 白带黄色成豆腐渣| 一二三四中文在线观看免费高清| 国产精品1区2区在线观看.| 日本一二三区视频观看| 欧美三级亚洲精品| 日韩 亚洲 欧美在线| 亚洲成人精品中文字幕电影| 亚洲精品乱码久久久v下载方式| 国产欧美另类精品又又久久亚洲欧美| 久久久精品免费免费高清| 天美传媒精品一区二区| 色网站视频免费| 水蜜桃什么品种好| 日本色播在线视频| 久久精品国产亚洲av天美| 又粗又硬又长又爽又黄的视频| 91精品一卡2卡3卡4卡| 亚洲国产精品成人综合色| 午夜福利在线观看免费完整高清在| 水蜜桃什么品种好| 欧美日韩在线观看h| 免费黄网站久久成人精品| 国产男人的电影天堂91| 尾随美女入室| 欧美xxxx黑人xx丫x性爽| 久久久久久九九精品二区国产| 国产精品伦人一区二区| 久久久久久伊人网av| 亚洲av中文字字幕乱码综合| 69人妻影院| 亚洲第一区二区三区不卡| 国产黄色视频一区二区在线观看|