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

    滲流與側(cè)向流共同作用下的潮溝邊壁侵蝕研究進(jìn)展*

    2022-07-28 00:58:10王客予張凱麗
    海洋與湖沼 2022年4期
    關(guān)鍵詞:影響實(shí)驗(yàn)研究

    龔 政 王客予, 2 趙 堃 張凱麗

    滲流與側(cè)向流共同作用下的潮溝邊壁侵蝕研究進(jìn)展*

    龔 政1王客予1, 2①趙 堃1張凱麗1

    (1. 河海大學(xué)江蘇省海岸海洋資源開發(fā)與環(huán)境安全重點(diǎn)實(shí)驗(yàn)室 江蘇南京 210098; 2. 中國(guó)電建集團(tuán)西北勘測(cè)設(shè)計(jì)研究院有限公司 陜西西安 710065)

    潮溝的邊壁侵蝕可以分為水流作用下的邊壁沖刷與土塊重力作用下的邊壁崩塌。前人研究表明, 水流流速、水壓及土體干密度等因素均對(duì)底床沖刷率有顯著影響, 然而側(cè)向邊壁的沖刷率量測(cè)以及相關(guān)影響因素問題亟待解決。前人較多關(guān)注側(cè)向流對(duì)于邊壁崩塌的影響, 而對(duì)滲流作用以及滲流與側(cè)向流在邊壁侵蝕過程中的共同作用研究較少。近期研究表明, 滲流已成為邊壁侵蝕過程中的重要一環(huán)。因此, 滲流與側(cè)向流不同動(dòng)力作用下, 邊壁侵蝕的差異性和相似性, 以及兩種動(dòng)力共同作用下的侵蝕機(jī)制仍待進(jìn)一步探究。

    潮溝; 沖刷率; 邊壁侵蝕; 滲流

    潮灘地處海陸之交, 一般發(fā)育在沿海平原外緣, 廣泛分布在開敞式、港灣型和河口灣型海岸, 比如: 中國(guó)東部沿海(Zhang, 2016; Chen, 2017)、英國(guó)西部及東南海岸(Chen, 2012)、荷蘭西北海岸(Kleinhans, 2009)、美國(guó)東海岸(Mariotti, 2013), 具有寬廣性、尺度大、坡度緩、底質(zhì)顆粒細(xì)(由淤泥質(zhì)黏土、粉砂、粉細(xì)砂等組成)等特征(呂亭豫等, 2016)。在淤泥質(zhì)潮灘的潮間帶, 即平均大潮高潮線至平均大潮低潮線之間的潮灘, 由海洋動(dòng)力, 特別是潮汐作用形成的潮溝系統(tǒng)廣泛發(fā)育并呈現(xiàn)樹枝狀、矩形狀、平形狀或羽狀等平面形態(tài)結(jié)構(gòu)(Ichoku, 1994; 呂亭豫等, 2016)。潮溝為水流及泥沙提供輸入與輸出潮灘的通道, 對(duì)潮灘的地貌動(dòng)力過程具有重要影響, 隨水流運(yùn)動(dòng)的營(yíng)養(yǎng)物質(zhì)對(duì)生態(tài)系統(tǒng)穩(wěn)定有重要意義, 同時(shí)潮溝也是潮灘上海陸相互作用最活躍的微地貌單元(Fagherazzi, 2012; Lanzoni, 2015)。

    經(jīng)典理論將邊壁侵蝕從機(jī)理上分為兩個(gè)階段: (1)水流作用下的邊壁沖刷(水力學(xué)過程); (2) 土塊重力引起的邊壁崩塌(土力學(xué)過程) (Rinaldi, 2007)。研究表明, 潮溝邊壁崩塌是海岸地貌動(dòng)力過程中重要的一環(huán), 對(duì)海岸工程建設(shè)構(gòu)成嚴(yán)重威脅, 例如2015年江蘇沿海條子泥海域潮溝曲流向岸擺動(dòng)距離達(dá)1 800 m, 引發(fā)了一系列應(yīng)急搶險(xiǎn)工程。此外, 大江大河崩岸現(xiàn)象也十分普遍(夏軍強(qiáng)等, 2003; 余文疇等, 2008; 盧金友等, 2017; 假冬冬等, 2021), 由側(cè)向流沖刷形成懸臂及邊壁崩塌現(xiàn)場(chǎng)見圖1。據(jù)不完全統(tǒng)計(jì), 近60年來長(zhǎng)江中下游累計(jì)崩岸長(zhǎng)度超過1 600 km, 占岸線總長(zhǎng)的37.6% (褚明華等, 2016; 張幸農(nóng)等, 2021)。此外, 崩塌土塊也是河道重要的供沙來源, 例如在美國(guó)中西部的黃土區(qū)域, 邊壁供沙占整個(gè)流域供沙的80%以上(Simon, 2000)。

    目前研究邊壁侵蝕的主要困難之一是缺乏土壤性質(zhì)數(shù)據(jù), 需要改進(jìn)相關(guān)性質(zhì)的測(cè)量方法。雖然目前用于穩(wěn)定性分析的一些土壤性質(zhì), 如液體、沉積物的密度容易量測(cè), 但其他性質(zhì)和系數(shù)獲取難度大, 也鮮有現(xiàn)場(chǎng)表征這些性質(zhì)和系數(shù)的方法。即使是主要的侵蝕特性(如沖刷系數(shù)、剪切強(qiáng)度、有效黏聚力和臨界起動(dòng)切應(yīng)力)也難測(cè)量, 不能很好地理解它們與其他土壤特性的相互關(guān)系, 因此邊壁沖刷速率的量測(cè)以及相關(guān)公式的推算成為亟待解決的問題之一。研究邊壁侵蝕的另一個(gè)難點(diǎn)是各種機(jī)制之間的相互作用。因此, 本文重點(diǎn)回顧邊壁沖刷速率量測(cè)及其影響因素和滲流和側(cè)向流共同作用下岸壁侵蝕機(jī)制兩方面的研究進(jìn)展, 揭示潮灘-潮溝系統(tǒng)多因子共同作用下的崩岸機(jī)理, 對(duì)于保障海岸工程安全, 以及潮灘系統(tǒng)的科學(xué)開發(fā)、利用和保護(hù)具有現(xiàn)實(shí)意義, 對(duì)促進(jìn)多學(xué)科交叉具有重要的理論和實(shí)際價(jià)值。

    圖1 側(cè)向流沖刷形成懸臂及邊壁崩塌

    注: a. 側(cè)向流沖刷岸底形成懸臂; b. 邊壁崩塌, 改自文獻(xiàn)(Samadi, 2013)

    1 邊壁沖刷

    邊壁沖刷是指水流在特定條件下對(duì)側(cè)向土體進(jìn)行沖刷的現(xiàn)象, 其主要影響因素為水流條件和土體本身的性質(zhì)。目前對(duì)于側(cè)向邊壁的土體起動(dòng)及沖刷研究較少, 大部分研究集中于底床起動(dòng)和沖刷方面, 其實(shí)驗(yàn)方法及量測(cè)手段可以在側(cè)向邊壁起動(dòng)沖刷的研究中加以借鑒?,F(xiàn)有兩種主流的底床土體沖刷率計(jì)算公式, 其中一種在曹叔尤公式(曹叔尤等, 1986)的基礎(chǔ)上將沖刷率與過量剪應(yīng)力確定成指數(shù)關(guān)系, 即:

    另一種沖刷率公式認(rèn)為水流切應(yīng)力大于土體的起動(dòng)切應(yīng)力時(shí), 土體才會(huì)被水流沖刷(Hanson, 2001):

    =d?(0–c), (2)

    式中,為沖刷厚度, 單位為m/s;d為土體的沖刷系數(shù), 單位為m3/(N?s); ?為沖刷時(shí)間, 單位為s;0為水流的平均切應(yīng)力, 單位為Pa;c為臨界起動(dòng)切應(yīng)力, 單位為Pa, 相應(yīng)的沖刷率為:

    =d(0–c). (3)

    黏性土的抗沖特性主要涉及沖刷率和起動(dòng)條件等, 而土體起動(dòng)除水流作用外, 主要受到土體本身物理及力學(xué)性質(zhì)的影響, 因此水流條件與土體本身性質(zhì)均應(yīng)在黏性土的抗沖研究中加以考慮(宗全利等, 2014)。在邊壁沖刷過程中, 水流入滲會(huì)導(dǎo)致土體本身密度、含水率及液塑限等物理性質(zhì)發(fā)生改變, 進(jìn)一步影響土體的內(nèi)摩擦角及黏聚力等抗剪強(qiáng)度指標(biāo), 土體臨界起動(dòng)切應(yīng)力隨之變化。

    1.1 流速對(duì)邊壁沖刷的影響

    黃建維(1989)在水槽實(shí)驗(yàn)的基礎(chǔ)上, 結(jié)合前人的研究成果, 認(rèn)為均勻容重底床的沖刷率與水流流速近似呈四次方關(guān)系。李華國(guó)等(1995)通過水槽實(shí)驗(yàn)得出了類似的結(jié)論, 認(rèn)為淤泥沖刷率與水流切應(yīng)力近似呈二次方關(guān)系, 與水流流速呈四次方關(guān)系。孫志林等(2011)在錢塘江河口河床獲取2個(gè)鉆孔共28個(gè)巖芯, 進(jìn)行起動(dòng)和沖刷水槽實(shí)驗(yàn), 發(fā)現(xiàn)未充分固結(jié)的淤泥沖刷率與相對(duì)剩余切應(yīng)力的二次方呈正比(與流速的四次方近似相關(guān)), 而固結(jié)已久的淤泥沖刷率與相對(duì)剩余切應(yīng)力呈線性關(guān)系(與流速的二次方近似相關(guān)), 二者的沖刷常數(shù)相差一個(gè)數(shù)量級(jí)。

    1.2 水壓對(duì)邊壁沖刷的影響

    萬兆惠等(1990)通過有壓管道實(shí)驗(yàn)表明, 粗顆粒泥沙在較大水深條件下的起動(dòng)流速幾乎不受水壓的影響; 而黏性細(xì)顆粒泥沙的起動(dòng)流速與水壓呈正比。甘甜(2017)比較了水壓在不同水力條件下(例如管道和明渠)對(duì)黏性泥沙起動(dòng)的影響, 實(shí)驗(yàn)結(jié)果表明: 在有壓管道中, 泥沙淤積固結(jié)初期, 黏結(jié)力對(duì)黏性泥沙的起動(dòng)影響不顯著, 阻礙黏性泥沙起動(dòng)的主要因素為管道壓力, 在后期, 黏性泥沙黏結(jié)力增強(qiáng), 對(duì)其起動(dòng)的影響增大。在明渠大水深條件下, 相對(duì)水深變化小, 無法對(duì)黏性泥沙的起動(dòng)造成明顯影響, 然而在明渠水槽中, 相對(duì)水深成為影響?zhàn)ば阅嗌称饎?dòng)的主要因素。明渠中水深的增大改變了土體微團(tuán)所受薄膜水附加壓力, 起動(dòng)流速和起動(dòng)切應(yīng)力隨水深的增大而增大。

    1.3 土體密度對(duì)邊壁沖刷的影響

    干密度是黏性土顆粒在沉積過程中對(duì)外表現(xiàn)的宏觀結(jié)果, 易于量測(cè)獲得, 能綜合反映出土體在沉降密實(shí)過程中土顆粒間黏結(jié)力的變化。干密度的差異反映出土體固結(jié)度的不同, 土體的固結(jié)度在宏觀層面影響對(duì)外表現(xiàn)的抗沖刷能力, 在微觀層面上對(duì)應(yīng)不同的顆粒間黏結(jié)力。顆粒間黏結(jié)力的大小影響床面泥沙的起動(dòng)現(xiàn)象: 土體干密度較小時(shí), 土顆粒間黏結(jié)力小, 床面泥沙以單顆粒起動(dòng)為主; 孔隙水逐漸排除的過程中土體干密度逐漸增大, 土顆粒的骨架發(fā)生壓縮, 黏結(jié)力的增大導(dǎo)致床面起動(dòng)由單顆粒向微團(tuán)發(fā)生過渡; 當(dāng)土體干密度(或土顆粒間黏結(jié)力)增大到一定程度時(shí), 床面泥沙起動(dòng)則完全表現(xiàn)出微團(tuán)的運(yùn)動(dòng)形式(舒彩文等, 2007)。

    Krone(1999)通過對(duì)前人的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行分析后, 認(rèn)為底床抗沖刷能力隨著床層結(jié)構(gòu)沖刷的進(jìn)展而增加, 底床沖刷率與土體密度呈線性關(guān)系, 當(dāng)顆粒結(jié)構(gòu)幾乎完全沖刷到密度較高的結(jié)構(gòu)時(shí), 該線性關(guān)系的斜率及截距均發(fā)生較大變化。舒彩文等(2007)在前人實(shí)驗(yàn)資料的基礎(chǔ)上研究了土體干容重及粒徑對(duì)臨界起動(dòng)切應(yīng)力的影響, 結(jié)果表明, 黏性泥沙的沖刷率與干容重及顆粒間的黏性正相關(guān), 與粒徑呈負(fù)相關(guān)關(guān)系; 對(duì)于淤泥, 土體沖刷率與干容重則可以量化為六次方的關(guān)系, 基于此發(fā)現(xiàn)建立了相關(guān)的泥沙臨界起動(dòng)切應(yīng)力公式(舒彩文等, 2007)。

    1.4 沖刷系數(shù)與起動(dòng)切應(yīng)力的計(jì)算

    從式(2)和式(3)可以看出, 黏性土的沖刷率不僅受到水流流速的影響, 而且與土體的沖刷系數(shù)和臨界起動(dòng)切應(yīng)力密切相關(guān)。土體本身的特性在很大程度上影響沖刷系數(shù)和臨界起動(dòng)切應(yīng)力, 其中前者是決定土體沖刷率的最主要因素。

    國(guó)內(nèi)外學(xué)者分別選取不同地區(qū)及植被覆蓋條件的土樣, 通過實(shí)驗(yàn)得到了相應(yīng)的沖刷系數(shù)與土體臨界起動(dòng)切應(yīng)力的數(shù)量關(guān)系式, 可歸納為d=·c的形式, 其中、均為常數(shù)。Hanson等(2001)在進(jìn)行了83組水下射流實(shí)驗(yàn)之后, 得到土體沖刷系數(shù)與臨界起動(dòng)切應(yīng)力的定量關(guān)系式為d=2×10–7c–0.5; Wynn等(2004)針對(duì)位于美國(guó)弗吉尼亞州西南部的河岸土體進(jìn)行了142組實(shí)驗(yàn), 得到該地區(qū)受到植被覆蓋的土體沖刷系數(shù)與臨界起動(dòng)切應(yīng)力的定量關(guān)系式為d=3.1× 10–6c–0.37; Karmaker等(2011)在印度雅魯藏布江河岸上進(jìn)行58次現(xiàn)場(chǎng)淹沒射流實(shí)驗(yàn)后, 得到土體沖刷系數(shù)與臨界起動(dòng)切應(yīng)力的定量關(guān)系式為d=3.1× 10–6c–0.185; 宗全利等(2014)對(duì)荊江段河岸黏性土進(jìn)行沖刷實(shí)驗(yàn)后, 通過數(shù)據(jù)擬合得到該河段邊壁黏性土沖刷系數(shù)與臨界起動(dòng)切應(yīng)力關(guān)系式為d=7.677× 10–6c–1.949。

    不難發(fā)現(xiàn), 黏性土體的沖刷系數(shù)與臨界起動(dòng)切應(yīng)力之間確實(shí)存在某種數(shù)量關(guān)系, 但不同的植被覆蓋條件與土體本身性質(zhì)決定了該關(guān)系式中、兩個(gè)常數(shù)的大小。含水率增大導(dǎo)致黏性土體塑性增強(qiáng), 土體沖刷難度增大, 這是Karmaker與宗全利等實(shí)驗(yàn)結(jié)果中沖刷系數(shù)差異大的原因之一。此外, 宗全利等(2014)實(shí)驗(yàn)土體黏粒含量c=24.6%, Hanson、Wynn、Karmaker等實(shí)驗(yàn)土體黏粒含量c均在50%以上, 黏粒含量減小則土體黏性減小, 越容易沖刷, 所以沖刷系數(shù)越大。另外, 土樣來源及實(shí)驗(yàn)條件也會(huì)影響實(shí)驗(yàn)結(jié)果, Hanson、Wynn、Karmaker等實(shí)驗(yàn)均在現(xiàn)場(chǎng)進(jìn)行, 但Hanson等實(shí)驗(yàn)關(guān)注床面土體, Wynn、Karmaker等實(shí)驗(yàn)針對(duì)河岸土體起動(dòng), 且水流對(duì)河岸的沖刷作用受到覆蓋在其上的植被影響, 在一定程度上導(dǎo)致了關(guān)系式中的沖刷系數(shù)計(jì)算結(jié)果偏小; 宗全利等在實(shí)驗(yàn)前將黏性土與水?dāng)嚢杈鶆? 靜置一段時(shí)間使土體在自重作用下固結(jié), 以此得到實(shí)驗(yàn)土樣, 而Hanson、Wynn等實(shí)驗(yàn)結(jié)果均是針對(duì)原狀土體現(xiàn)場(chǎng)測(cè)得, 不同的土樣獲取方式在一定程度上影響計(jì)算結(jié)果。

    通過對(duì)土體沖刷及臨界起動(dòng)切應(yīng)力影響因素的回顧, 發(fā)現(xiàn)大部分研究集中于底床起動(dòng)和沖刷方面, 對(duì)于側(cè)向邊壁的土體起動(dòng)及沖刷研究較少, 且大部分研究對(duì)于流速、水壓、土體密度等因素對(duì)于土體沖刷速率的影響意見不一。鑒于側(cè)向邊壁土顆粒受力與底床土顆粒受力不盡相同(側(cè)壁顆粒由于重量分量的作用沿坡面向下滑動(dòng), 是邊壁坍塌的驅(qū)動(dòng)力之一; 而底床顆粒則受到沿水流方向的拖曳力驅(qū)動(dòng)作用(Aldefae, 2020), 建議后續(xù)對(duì)側(cè)向邊壁沖刷率的影響因素進(jìn)行實(shí)驗(yàn)分析, 提出相關(guān)沖刷率計(jì)算公式。此外, 應(yīng)盡量采用現(xiàn)場(chǎng)取樣及測(cè)試, 以便得到更加貼近實(shí)際的實(shí)驗(yàn)結(jié)果。

    2 滲流對(duì)邊壁侵蝕的影響

    在較長(zhǎng)一段時(shí)間里, 滲流的隱蔽性和復(fù)雜性是研究其影響邊壁侵蝕的難點(diǎn)所在。Hagerty和Fox等列舉了一些原因: 滲流侵蝕是一種非常復(fù)雜的機(jī)理, 在一次暴雨后, 許多不同的條件都會(huì)造成嚴(yán)重的滲流侵蝕, 而在另一次暴雨之后, 滲流侵蝕的規(guī)??梢院雎圆挥?jì)(Hagerty, 1991; Fox, 2010)。此外, 其他邊壁侵蝕的機(jī)制可能會(huì)將滲流造成的侵蝕掩蓋, 例如河岸的侵蝕可能被誤解為完全由側(cè)向流侵蝕造成。在近期的研究中, 與地下水流動(dòng)相關(guān)的孔隙水含量和壓力的變化被認(rèn)為是控制堤岸坍塌發(fā)生的最重要因素之一(Simon, 2000; Darby, 2007; Rinaldi, 2008, 2013)。然而迄今為止發(fā)展起來的數(shù)學(xué)模型雖可用于獲得河床和河岸變形, 但僅適用于有限的邊界條件。此外, 由于邊壁侵蝕現(xiàn)象缺乏局部的詳細(xì)數(shù)據(jù), 對(duì)于坡度和水流等相關(guān)條件的數(shù)值模擬存在一定困難。因此, 需要更多的現(xiàn)場(chǎng)或?qū)嶒?yàn)室研究來標(biāo)定和驗(yàn)證多種條件下與河岸侵蝕相關(guān)的動(dòng)態(tài)特征(Shu, 2019)。

    2.1 滲流影響邊壁侵蝕的實(shí)驗(yàn)研究

    不少學(xué)者通過滲漏儀實(shí)驗(yàn)(Fox, 2006, 2007; Lindow, 2009; Karmaker, 2013)揭示了滲流導(dǎo)致邊壁侵蝕的基本過程: (1) 滲流底切; (2) 在岸坡或岸頂表面形成拉裂; (3) 岸坡坍塌。滲漏儀實(shí)驗(yàn)?zāi)M滲流作用下的邊壁侵蝕過程見圖2。Fox等(2007)通過實(shí)驗(yàn)發(fā)現(xiàn), 滲流侵蝕的作用略大于土壤孔隙水壓力增加在邊壁侵蝕中的作用。拉裂縫是滲流過程中土體基質(zhì)吸力的降低和底切誘導(dǎo)彎矩共同作用的結(jié)果。邊壁坍塌后形成的岸趾對(duì)邊壁起到一定的保護(hù)作用, 阻礙了側(cè)向流的進(jìn)一步侵蝕和滲流的底切過程。只有側(cè)向流或人工將岸趾移除之后, 滲流底切才會(huì)繼續(xù)進(jìn)行, 這項(xiàng)工作表明了今后將河流侵蝕和滲流侵蝕同時(shí)進(jìn)行研究的必要性(Midgley, 2013)。邊壁崩塌土體在河床上的輸移范圍及程度, 關(guān)系到邊壁的二次崩塌以及河床的沖淤變形, 而后者的改變又會(huì)反過來影響邊壁的穩(wěn)定(Fox, 2014)。

    圖2 滲漏儀模擬滲流作用下的邊壁侵蝕過程

    注: 改自文獻(xiàn)(Fox, 2010)

    上述實(shí)驗(yàn)研究大多采用二維滲漏儀進(jìn)行, 忽略了崩岸沿寬度方向的差異, 且二維模型無法模擬非均勻懸臂破壞、黏性岸坡附近和內(nèi)部復(fù)雜的水流特性、泥沙運(yùn)動(dòng)的縱向坡度、孔隙水壓力和崩塌塊體效應(yīng)(Fox, 2007; Patsinghasanee, 2017)。此外, SLOPE/W數(shù)值模型預(yù)測(cè)的邊壁坍塌時(shí)間比在滲漏儀實(shí)驗(yàn)中觀察到的更早, 可能原因是土體受到了滲漏儀兩側(cè)邊壁的壓縮力(Fox, 2007)。Chu-Agor等(2008)在Fox等(2006)二維模型實(shí)驗(yàn)的基礎(chǔ)上進(jìn)行改進(jìn), 使用三維模型進(jìn)行滲漏實(shí)驗(yàn), 在此基礎(chǔ)上建立了泥沙輸移函數(shù), 但仍存在滲流入口集中, 模型比尺較小, 無法真實(shí)還原土體應(yīng)力的變化過程等問題。因此, 需建立滲流影響邊壁侵蝕的三維原型實(shí)驗(yàn)系統(tǒng)進(jìn)行進(jìn)一步研究。

    2.2 滲流作用下的邊壁侵蝕機(jī)制

    滲流可以通過三種不同的機(jī)制導(dǎo)致邊壁侵蝕:

    (1) 增加土壤孔隙水壓力, 降低土體的抗剪強(qiáng)度(Fox, 2014)。室內(nèi)水槽研究表明, 當(dāng)非飽和孔隙水壓力增加到接近飽和時(shí), 土壤最初的臨界剪應(yīng)力或抗侵蝕能力有所提高, 而抗剪強(qiáng)度明顯下降, 地表侵蝕速率增加了一個(gè)數(shù)量級(jí)(Khanal, 2020)。

    (2) 滲流梯度力(Fox, 2014)。當(dāng)滲流壓力大于土的抗剪強(qiáng)度時(shí), 土體會(huì)發(fā)生拉伸或“突出”破壞, 而土體的抗剪強(qiáng)度因含水量增大而降低, 滲流力的增大和土體抗剪強(qiáng)度降低是實(shí)驗(yàn)中觀察到的拉伸或“突出”破壞的原因(Chu-Agor, 2008)??紫端畨毫突|(zhì)吸力均隨地下水位的變化而變化, 前者起到侵蝕作用, 后者起到表觀黏聚力的作用。在洪水退水或落潮期, 靜水壓力迅速下降, 平均地下水位的滯后變化導(dǎo)致孔隙水壓力相對(duì)較高, 不利于岸坡穩(wěn)定。因此, 堤岸坍塌大多發(fā)生在水位下降期(Deng, 2019)。相關(guān)實(shí)驗(yàn)發(fā)現(xiàn), 增大滲透水頭, 黏性土可蝕性隨之增加, 土體內(nèi)濕潤(rùn)鋒推進(jìn)速度加快, 崩岸速度加快; 對(duì)于同一滲透水頭, 隨著岸高的降低, 崩岸泥沙質(zhì)量減少, 但滲流侵蝕對(duì)岸坡泥沙總量的貢獻(xiàn)率增加(Fox, 2006; Al-Madhhachi, 2011)。

    (3) 滲流顆粒的移動(dòng)和底切(Fox, 2014)。滲流梯度力小于土塊的初始抗力時(shí), 顆粒移動(dòng)(即滲透侵蝕)和土體下切導(dǎo)致邊壁崩塌, 這是由滲流和孔隙水壓力累積的合力造成的(Chu-Agor, 2008)。

    2.3 滲流與側(cè)向流共同作用下的邊壁侵蝕機(jī)制

    滲流與其他導(dǎo)致邊壁侵蝕(即河流侵蝕、圍壓和植被)的因子間復(fù)雜相互作用, 使得人們很難充分了解滲流對(duì)邊壁侵蝕的作用(Fox, 2006)。復(fù)雜的原因主要在于難以量化影響崩岸這一現(xiàn)象的土壤、水文和巖土機(jī)制, 以及理解它們之間的相互作用。相關(guān)研究表明, 滲流本身似乎不是造成侵蝕的主要因素。然而, 當(dāng)有其他因素(如降雨)引起侵蝕時(shí), 滲流會(huì)影響侵蝕速率。滲流對(duì)侵蝕的影響最能通過土體可蝕性的變化來確定, 表明滲流與土體可蝕性有直接關(guān)系。這表明土體可蝕性既受土體本身特性的影響, 又受土體中水流條件的影響(Owoputi, 2001)。

    近年來, 一些學(xué)者研究了側(cè)向流對(duì)河岸邊壁侵蝕的影響。Samadi等(2013)采用人工開挖的方式代替?zhèn)认蛄鬟M(jìn)行實(shí)驗(yàn), 并將數(shù)值模型得到的應(yīng)力和形變分布與實(shí)驗(yàn)室觀測(cè)結(jié)果進(jìn)行了比較, 研究發(fā)現(xiàn)邊壁發(fā)生傾倒破壞的可能性大于剪切破壞。Patsinghasanee等(2017)通過對(duì)泰國(guó)U-Tapao河岸的實(shí)驗(yàn)結(jié)果和現(xiàn)場(chǎng)數(shù)據(jù)進(jìn)行數(shù)值模擬, 對(duì)懸臂破壞的機(jī)理進(jìn)行了分析研究, 開發(fā)了可以較準(zhǔn)確地模擬邊壁破壞機(jī)制與懸臂形狀的數(shù)值模型。Thi等(2019)關(guān)注側(cè)向流的漲落速率與邊壁土體性質(zhì)的影響, 研究發(fā)現(xiàn)水位不變時(shí), 土體密度更大、黏聚力更高的邊壁容易出現(xiàn)更大更深的崩塌破壞, 水位變化率低時(shí), 側(cè)向流的入滲成為邊壁侵蝕的主要因素, 而水位變化率高時(shí), 邊壁侵蝕以懸臂的傾倒破壞為主。Aldefae等(2021)發(fā)現(xiàn)邊壁土體侵蝕率隨著側(cè)向流流速的下降而降低。

    Van Eerdt(1985)和Gabet(1998)對(duì)潮汐環(huán)境中的崩岸進(jìn)行了開創(chuàng)性的工作, 用梁破壞分析方法研究了懸臂形鹽沼堤岸和懸崖的穩(wěn)定性。Ginsberg等(1990)證明了崩岸是如何導(dǎo)致侵蝕尖頂形狀的發(fā)展的, 其研究大坡角是潮汐槽岸崩塌的主要原因。Zhao等(2020)主要關(guān)注潮溝邊壁高度與近岸水深的作用, 結(jié)果表明邊壁的后退速率隨邊壁高度與近岸水深之比的增大而先增大后減小。Zhang等(2021)對(duì)潮溝邊壁破壞面的三維結(jié)構(gòu)進(jìn)一步分析后發(fā)現(xiàn)破壞面角與岸坡后退距離呈線性負(fù)相關(guān)。

    前人對(duì)于邊壁侵蝕的研究多傾向于孤立地看待單個(gè)過程, 但邊壁侵蝕是滲流、側(cè)向流等多因素相互作用的綜合結(jié)果(Wolman, 1959; Thorne, 1982; Lawler, 1997)。簡(jiǎn)而言之, 孤立地看待側(cè)向流或滲流等因素在邊壁侵蝕中的作用, 會(huì)導(dǎo)致對(duì)相關(guān)過程的預(yù)測(cè)出現(xiàn)偏差(Darby, 2007)。因此, 多位學(xué)者提出建議, 今后的研究應(yīng)嘗試將滲流侵蝕與側(cè)向流侵蝕相結(jié)合, 對(duì)崩岸機(jī)理的認(rèn)識(shí)應(yīng)涵蓋側(cè)向流平均水位漲落速率和滲流水頭等多方面因素(Fox, 2010; Karmaker, 2013; Thi, 2019)。

    對(duì)滲流與側(cè)向流共同作用下邊壁侵蝕機(jī)制的研究較少, 部分學(xué)者通過對(duì)比研究發(fā)現(xiàn), 在只有側(cè)向流作用于邊壁的情況下, 土體內(nèi)孔隙水壓力從試驗(yàn)開始到結(jié)束一直保持穩(wěn)定(Ning, 2011)。側(cè)向流一方面沖刷岸底促進(jìn)懸臂的形成, 加速邊壁坍塌的過程; 另一方面通過對(duì)邊壁由外向內(nèi)的靜水壓力在一定程度上支撐邊壁, 在此過程中水流的入滲加強(qiáng)了邊壁土體的抗滑移能力(葉威等, 2019), 水流持續(xù)時(shí)間長(zhǎng)短影響沖刷崩塌過程。在側(cè)向流和滲流的共同作用下, 土體內(nèi)孔隙水壓力開始時(shí)增長(zhǎng)較快, 然后增長(zhǎng)速度逐漸減小, 最后達(dá)到穩(wěn)定狀態(tài)。相比于只存在側(cè)向流的情況, 兩種因素共同作用下邊壁崩塌更容易發(fā)生。在相同流速、不同滲透壓力的條件下, 滲流水力坡度越大, 破壞發(fā)生的時(shí)間越早, 破壞的范圍也越大(Ning, 2011)。

    3 河岸與潮溝邊壁侵蝕的差異性

    河岸與潮溝系統(tǒng)擁有相似之處, 二者均具有復(fù)雜的樹枝網(wǎng)狀結(jié)構(gòu)及蜿蜒曲流的特性(Coco, 2013), 在潮溝系統(tǒng)的研究中借鑒河岸侵蝕的方法基本可行, 但需要考慮兩個(gè)系統(tǒng)的差異所帶來的影響。主要的區(qū)別在于二者水流驅(qū)動(dòng)力的不同: 潮溝內(nèi)的水流受潮汐及波浪影響較大(Hibma, 2004; De Swart, 2009; Bendoni, 2014), 潮流以潮周期為時(shí)間尺度, 在潮溝內(nèi)往復(fù)運(yùn)動(dòng), 流速變化劇烈, 其最大流速時(shí)刻為漲潮初期及落潮后期(Zhang, 2016); 河流主要受地形和上游徑流影響(Bendoni, 2014), 以季節(jié)為時(shí)間尺度, 單向流態(tài)穩(wěn)定, 其最大流速與高水位同時(shí)出現(xiàn)(Bayliss-Smith, 1979)。Xin等(2011)的數(shù)值模型研究強(qiáng)調(diào)了潮灘水文研究的復(fù)雜機(jī)制, 發(fā)現(xiàn)在潮周期內(nèi), 滲流存在明顯的不對(duì)稱性, 在由海向陸的方向上潮灘土體內(nèi)孔隙水循環(huán)的特征時(shí)間尺度逐漸減少多個(gè)數(shù)量級(jí)。其研究表明滲流過程在一定的空間尺度下對(duì)于潮溝形態(tài)的演變具有重要作用。往復(fù)流的存在使得潮溝邊壁土壤飽和度始終保持較高水平, 自陸向海的土體淹沒歷時(shí)不同導(dǎo)致土體強(qiáng)度與破壞頻率的空間差異(D’Alpaos, 2006)。Simon等(2000)認(rèn)為邊壁坍塌主要由洪水衰退期的靜水壓力迅速喪失導(dǎo)致。與河流環(huán)境下二元河岸土體各向異性沿邊壁剖面方向不同的是(夏軍強(qiáng)等, 2013; 鄧珊珊等, 2020), 潮汐環(huán)境中潮溝邊壁土體受往復(fù)流影響而導(dǎo)致的各項(xiàng)異性主要表現(xiàn)在垂直于海岸的方向。風(fēng)浪和鹽度等因素亦可通過改變土體的孔隙水壓力分布(Francalanci, 2013)或土壤分散過程(Masoodi, 2019)影響潮溝邊壁的穩(wěn)定性。因此, 漲落潮的頻繁出現(xiàn)對(duì)潮溝邊壁侵蝕產(chǎn)生顯著的影響。

    4 結(jié)論和建議

    目前研究邊壁侵蝕的主要困難之一是缺乏土壤性質(zhì)數(shù)據(jù), 需要改進(jìn)相關(guān)侵蝕特性與土壤性質(zhì)的測(cè)量方法。前人研究表明, 底床泥沙沖刷速率與流速正相關(guān), 而與靜水壓力、土體干密度負(fù)相關(guān), 但就某一種因素對(duì)沖刷率的影響程度意見尚未統(tǒng)一, 且大多數(shù)研究針對(duì)底床沖刷, 較少研究關(guān)注邊壁的側(cè)向沖刷過程, 因此邊壁沖刷侵蝕速率的量測(cè)及其影響因素成為亟待解決的問題之一。邊壁侵蝕的另一個(gè)研究難點(diǎn)是各種機(jī)制之間的相互作用。例如, 前人大多研究側(cè)向流沖刷條件下邊壁的侵蝕坍塌, 而對(duì)滲流作用以及滲流與側(cè)向流在邊壁侵蝕過程中的共同作用研究較少。近期研究表明, 滲流侵蝕是導(dǎo)致邊壁侵蝕破壞的重要因素, 可以通過三種不同的機(jī)制導(dǎo)致邊壁侵蝕: (1) 增加土壤孔隙水壓力, 降低土體的抗剪強(qiáng)度; (2) 滲流梯度力; (3) 滲流顆粒的移動(dòng)和底切。滲流與側(cè)向流不同動(dòng)力作用下, 邊壁侵蝕的差異性和相似性, 以及兩種動(dòng)力共同作用下的侵蝕機(jī)制仍待進(jìn)一步探究。

    邊壁侵蝕的實(shí)驗(yàn)研究較多針對(duì)河流環(huán)境下尺度較大的分層邊壁, 且多使用降比尺的模型實(shí)驗(yàn)進(jìn)行探索, 而對(duì)受潮流影響的尺度相對(duì)小的潮溝邊壁關(guān)注不足。雖然二者的研究方法可以相互借鑒, 但潮溝環(huán)境中往復(fù)流的水力特性及細(xì)顆粒的泥沙性質(zhì)帶來的影響仍不可忽略。針對(duì)潮流、滲流對(duì)泥灘潮溝等海岸環(huán)境的影響, 建議從以下方面進(jìn)行探究: (1) 潮流與滲流不同動(dòng)力下, 潮溝邊壁侵蝕的差異性和相似性。(2) 潮流與滲流共同作用下潮溝邊壁的侵蝕機(jī)制。應(yīng)采用水槽實(shí)驗(yàn)與數(shù)值模擬手段, 考慮水流和土體性質(zhì)的影響, 優(yōu)化邊壁沖刷速率經(jīng)驗(yàn)公式, 基于實(shí)驗(yàn)建立滲流水力坡度與邊壁侵蝕速率的經(jīng)驗(yàn)公式, 分析邊壁侵蝕類型及破壞頻率和距離等指標(biāo), 明確侵蝕過程控制因素及貢獻(xiàn)率。(3) 考慮潮汐環(huán)境中往復(fù)流的水動(dòng)力特性及細(xì)顆粒的泥沙性質(zhì)對(duì)潮溝邊壁侵蝕的影響。(4) 植被根系較高的抗拉強(qiáng)度與吸水的特性對(duì)邊壁穩(wěn)定性影響顯著, 建議可以從植被與地下水的共同作用入手, 將生物因素納入多因子共同作用下的邊壁侵蝕研究中來。除滲流及側(cè)向流外, 植被、土壤性質(zhì)、波浪、底棲生物等因素均對(duì)潮溝邊壁穩(wěn)定性造成不同程度的影響, 對(duì)多因子驅(qū)動(dòng)下潮溝邊壁侵蝕的研究將有助于完善潮灘-潮溝邊壁侵蝕動(dòng)力地貌模型, 從本質(zhì)上闡明潮溝崩岸的力學(xué)機(jī)理與周期性的侵蝕機(jī)制, 對(duì)于保障海岸工程安全, 以及潮灘資源的科學(xué)開發(fā)、利用和保護(hù)具有現(xiàn)實(shí)意義, 對(duì)促進(jìn)多學(xué)科交叉具有重要的理論和實(shí)際價(jià)值。

    萬兆惠, 宋天成, 何青, 1990. 水壓力對(duì)細(xì)顆粒泥沙起動(dòng)流速影響的試驗(yàn)研究[J]. 泥沙研究(4): 62-69.

    鄧珊珊, 夏軍強(qiáng), 宗全利, 等, 2020. 下荊江典型河段蘆葦根系特性及其對(duì)二元結(jié)構(gòu)河岸穩(wěn)定的影響[J]. 泥沙研究, 45(5): 13-19.

    甘甜, 2017. 管道壓力與水深對(duì)黏性泥沙起動(dòng)影響分析[J]. 中國(guó)農(nóng)村水利水電(11): 143-147, 154.

    盧金友, 朱勇輝, 岳紅艷, 等, 2017. 長(zhǎng)江中下游崩岸治理與河道整治技術(shù)[J]. 水利水電快報(bào), 38(11): 6-14.

    葉威, 章光, 胡少華, 等, 2019. 河流沖刷作用下堤岸穩(wěn)定性演化分析[J]. 水利水電快報(bào), 40(1): 42-47.

    呂亭豫, 龔政, 張長(zhǎng)寬, 等, 2016. 粉砂淤泥質(zhì)潮灘潮溝形態(tài)特征及發(fā)育演變過程研究現(xiàn)狀[J]. 河海大學(xué)學(xué)報(bào)(自然科學(xué)版), 44(2): 178-188.

    孫志林, 張翀超, 黃賽花, 等, 2011. 黏性非均勻沙的沖刷[J]. 泥沙研究(3): 44-48.

    李華國(guó), 袁美琦, 張秀芹, 1995. 淤泥臨界起動(dòng)條件及沖刷率試驗(yàn)研究[J]. 水道港口(3): 20-26.

    張幸農(nóng), 假冬冬, 陳長(zhǎng)英, 2021. 長(zhǎng)江中下游崩岸時(shí)空分布特征與規(guī)律[J]. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào), 29(1): 55-63.

    宗全利, 夏軍強(qiáng), 張翼, 等, 2014. 荊江段河岸黏性土體抗沖特性試驗(yàn)[J]. 水科學(xué)進(jìn)展, 25(4): 567-574.

    夏軍強(qiáng), 王光謙, 吳保生, 2003. 黃河下游河床縱向與橫向變形的數(shù)值模擬——Ⅰ二維混合模型的建立[J]. 水科學(xué)進(jìn)展, 14(4): 389-395.

    夏軍強(qiáng), 宗全利, 許全喜, 等, 2013. 下荊江二元結(jié)構(gòu)河岸土體特性及崩岸機(jī)理[J]. 水科學(xué)進(jìn)展, 24(6): 810-820.

    黃建維, 1989. 黏性泥沙在鹽水中沖刷和沉降特性的試驗(yàn)研究[J]. 海洋工程(1): 61-70.

    曹叔尤, 杜國(guó)翰, 1986. 黏性土沖淤的試驗(yàn)研究[J]. 泥沙研究(4): 73-82.

    假冬冬, 楊俊, 陳長(zhǎng)英, 等, 2021. 中國(guó)北方季節(jié)性冰凍河流岸灘崩塌數(shù)值模擬——以松花江為例[J]. 水科學(xué)進(jìn)展, 32(5): 717-726.

    舒彩文, 王軍, 談廣鳴, 2007. 干容重對(duì)黏性淤積物起動(dòng)和沖刷的影響[J]. 武漢大學(xué)學(xué)報(bào)(工學(xué)版), 40(1): 25-28.

    褚明華, 黃先龍, 2016. 2016年長(zhǎng)江中下游崩岸應(yīng)急整治思考[J]. 中國(guó)水利(21): 10-11, 16.

    記得在井岡山革命烈士陵園,我們向長(zhǎng)眠在這里的革命烈士敬獻(xiàn)花圈,聆聽井岡山時(shí)期的艱苦斗爭(zhēng)歷史:“兩年零四個(gè)月的井岡山革命斗爭(zhēng)異常慘烈,犧牲烈士4.8萬人,平均每天就有50多人倒下……”面對(duì)巨大的犧牲,根據(jù)地軍民鐵骨錚錚,“就是剩下我一人,也要舉起紅旗,爬山越嶺干到底!”伍若蘭同志在敵人的嚴(yán)刑拷打面前寧死不屈:“共產(chǎn)黨人從來不怕死……若要我低頭,除非日從西邊出,贛江水倒流。”……一次次淚流滿面,一次次心中激蕩,深深被革命前輩的崇高品格和精神所震撼!這時(shí),真正體會(huì)到黨旗是用烈士的鮮血染紅的,堅(jiān)定的共產(chǎn)主義理想信念是中國(guó)革命紅旗不倒、立于不敗之地最根本的精神支柱。

    ALDEFAE A H, ALKHAFAJI R A, 2021. Experimental and numerical modeling to investigate the riverbank’s stability [J]. SN Applied Sciences, 3(2): 164.

    ALDEFAE A H, AL-KHAFAJI R A, SHAMKHI M S,, 2020. Erosion, sediments transport and riverbank stability: a review [J]. IOP Conference Series: Materials Science and Engineering, 901: 012014.

    AL-MADHHACHI A T, HANSON G J, FOX G A,, 2011. Measuring erodibility of cohesive soils using laboratory jet erosion tests [M] // BEIGHLEY II R E, KILLGORE N W. World Environmental and Water Resources Congress 2011. Palm Springs, California, United States: American Society of Civil Engineers.

    BAYLISS-SMITH T P, HEALEY R, LAILEY R,, 1979. Tidal flows in salt marsh creeks [J]. Estuarine and Coastal Marine Science, 9(3): 235-255.

    BENDONI M, FRANCALANCI S, CAPPIETTI L,, 2014. On salt marshes retreat: experiments and modeling toppling failures induced by wind waves [J]. Journal of Geophysical Research: Earth Surface, 119(3): 603-620.

    CHEN X D, ZHANG C K, PATERSON D M,, 2017. Hindered erosion: the biological mediation of noncohesive sediment behavior [J]. Water Resources Research, 53(6): 4787-4801.

    CHEN Y, THOMPSON C E L, COLLINS M B, 2012. Saltmarsh creek bank stability: biostabilisation and consolidation with depth [J]. Continental Shelf Research, 35: 64-74.

    CHU-AGOR M L, FOX G A, CANCIENNE R M,, 2008. Seepage caused tension failures and erosion undercutting of hillslopes [J]. Journal of Hydrology, 359(3/4): 247-259.

    COCO G, ZHOU Z, VAN MAANEN B,, 2013. Morphodynamics of tidal networks: advances and challenges [J]. Marine Geology, 346: 1-16.

    D'ALPAOS A, LANZONI S, MUDD S M,, 2006. Modeling the influence of hydroperiod and vegetation on the cross-sectional formation of tidal channels [J]. Estuarine, Coastal and Shelf Science, 69(3/4): 311-324.

    DARBY S E, RINALDI M, DAPPORTO S, 2007. Coupled simulations of fluvial erosion and mass wasting for cohesive river banks [J]. Journal of Geophysical Research: Earth Surface, 112(F3): F03022.

    DE SWART H E, ZIMMERMAN J T F, 2009. Morphodynamics of tidal inlet systems [J]. Annual Review of Fluid Mechanics, 41: 203-229.

    DENG S S, XIA J Q, ZHOU M R, 2019. Coupled two-dimensional modeling of bed evolution and bank erosion in the Upper JingJiang Reach of Middle Yangtze River [J]. Geomorphology, 344: 10-24.

    FAGHERAZZI S, KIRWAN M L, MUDD S M,, 2012. Numerical models of salt marsh evolution: ecological, geomorphic, and climatic factors [J]. Reviews of Geophysics, 50(1): RG1002.

    FOX G A, CHU-AGOR M L, WILSON G V, 2007. Erosion of noncohesive sediment by ground water seepage: lysimeter experiments and stability modeling [J]. Soil Science Society of America Journal, 71(6): 1822-1830.

    FOX G A, FELICE R G, 2014. Bank undercutting and tension failure by groundwater seepage: predicting failure mechanisms [J]. Earth Surface Processes and Landforms, 39(6): 758-765.

    FOX G A, WILSON G V, 2010. The role of subsurface flow in hillslope and stream bank erosion: a review [J]. Soil Science Society of America Journal, 74(3): 717-733.

    FOX G A, WILSON G V, PERIKETI R K,, 2006. Sediment transport model for seepage erosion of streambank sediment [J]. Journal of Hydrologic Engineering, 11(6): 603-611.

    FRANCALANCI S, BENDONI M, RINALDI M,, 2013. Ecomorphodynamic evolution of salt marshes: experimental observations of bank retreat processes [J]. Geomorphology, 195: 53-65.

    GABET E J, 1998. Lateral migration and bank erosion in a saltmarsh tidal channel in San Francisco Bay, California [J]. Estuaries, 21(4): 745-753.

    GINSBERG S S, PERILLO G M E, 1990. Channel bank recession in the Bahía Blanca estuary, Argentina [J]. Journal of Coastal Research, 6(4): 999-1010.

    HAGERTY D J, 1991. Piping/sapping erosion. I: basic considerations [J]. Journal of Hydraulic Engineering, 117(8): 991-1008.

    HANSON G J, SIMON A, 2001. Erodibility of cohesive streambeds in the loess area of the midwestern USA [J]. Hydrological Processes, 15(1): 23-38.

    HIBMA A, STIVE M J F, WANG Z B, 2004. Estuarine morphodynamics [J]. Coastal Engineering, 51(8/9): 765-778.

    ICHOKU C, CHOROWICZ J, 1994. A numerical approach to the analysis and classification of channel network patterns [J].Water Resources Research, 30(2): 161-174.

    KARMAKER T, DUTTA S, 2011. Erodibility of fine soil from the composite river bank of Brahmaputra in India [J]. Hydrological Processes, 25(1): 104-111.

    KARMAKER T, DUTTA S, 2013. Modeling seepage erosion and bank retreat in a composite river bank [J]. Journal of Hydrology, 476: 178-187.

    KHANAL A, FOX G A, GUERTAULT L, 2020. Soil moisture impacts linear and nonlinear erodibility parameters from jet erosion tests [J]. Transactions of the ASABE, 63(4): 1123-1131.

    KLEINHANS M G, SCHUURMAN F, BAKX W,, 2009. Meandering channel dynamics in highly cohesive sediment on an intertidal mud flat in the Westerschelde estuary, the Netherlands [J]. Geomorphology, 105(3/4): 261-276.

    KRONE R B, 1999. Effects of bed structure on erosion of cohesive sediments [J]. Journal of Hydraulic Engineering, 125(12): 1297-1301.

    LANZONI S, D'ALPAOS A, 2015. On funneling of tidal channels [J]. Journal of Geophysical Research: Earth Surface, 120(3): 433-452.

    LAWLER D M, THORNE C R, HOOKE J M, 1997. Bank erosion and instability [M] // THORNE C, HEY R, NEWSON M. Applied Fluvial Geomorphology for River Engineering and Management. New York: John Wiley and Sons Ltd: 137-172.

    LINDOW N, FOX G A, EVANS R O, 2009. Seepage erosion in layered stream bank material [J]. Earth Surface Processes and Landforms, 34(12): 1693-1701.

    MARIOTTI G, FAGHERAZZI S, 2013. Critical width of tidal flats triggers marsh collapse in the absence of sea-level rise [J]. Proceedings of the national Academy of Sciences, 110(14): 5353-5356.

    MASOODI A, TABATABAI M R M, NOORZAD A,, 2019. Riverbank stability under the influence of soil dispersion phenomenon [J]. Journal of Hydrologic Engineering, 24(3): 05019001.

    MIDGLEY T L, FOX G A, WILSON G V,, 2013. Seepage-induced streambank erosion and instability: in situ constant-head experiments [J]. Journal of Hydrologic Engineering, 18(10): 1200-1210.

    NING B, WU S Q, TAN Y F,, 2011. Coupling effect of seepage flow and river flow on the bank failure [J]. Journal of Hydrodynamics, 23(6): 834-840.

    OWOPUTI L O, STOLTE W J, 2001. The role of seepage in erodibility [J]. Hydrological Processes, 15(1): 13-22.

    PATSINGHASANEE S, KIMURA I, SHIMIZU Y,, 2017. Cantilever failure investigations for cohesive riverbanks [J]. Proceedings of the Institution of Civil Engineers - Water Management, 170(2): 93-108.

    RINALDI M, DARBY S E, 2007. 9 Modelling river-bank- erosion processes and mass failure mechanisms: progress towards fully coupled simulations [J]. Developments in Earth Surface Processes, 11: 213-239.

    RINALDI M, MENGONI B, LUPPI L,, 2008. Numerical simulation of hydrodynamics and bank erosion in a river bend [J]. Water Resources Research, 44(9): W09428.

    RINALDI M, NARDI L, 2013. Modeling interactions between riverbank hydrology and mass failures [J]. Journal of Hydrologic Engineering, 18(10): 1231-1240.

    SAMADI A, AMIRI-TOKALDANY E, DAVOUDI M H,, 2013. Experimental and numerical investigation of the stability of overhanging riverbanks [J]. Geomorphology, 184: 1-19.

    SHU A P, DUAN G S, RUBINATO M,, 2019. An experimental study on mechanisms for sediment transformation due to riverbank collapse [J]. Water, 11(3): 529.

    SIMON A, CURINI A, DARBY S E,, 2000. Bank and near-bank processes in an incised channel [J]. Geomorphology, 35(3/4): 193-217.

    THI T D, DO MINH D, 2019. Riverbank stability assessment under river water level changes and hydraulic erosion [J]. Water, 11(12): 2598.

    THORNE C R, 1982. Processes and mechanisms of river bank erosion [M] // HEY R D, BATHURST J C, THORNE C R. Gravel Bed Rivers. New York: Wiley: 227-271.

    VAN EERDT M M, 1985. Salt marsh cliff stability in the oosterschelde [J]. Earth Surface Processes and Landforms, 10(2): 95-106.

    WOLMAN M G, 1959. Factors influencing erosion of a cohesive river bank [J]. American Journal of Science, 257(3): 204-216.

    WYNN T M, MOSTAGHIMI S, ALPHIN E F, 2004. The effects of vegetation on stream bank erosion [C] American Society of Agricultural and Biological Engineers. 2004 ASAE Annual Meeting. Paper No.042226. https://doi.org/10. 13031/2013.16423.

    XIN P, YUAN L R, LI L,, 2011. Tidally driven multiscale pore water flow in a creek-marsh system [J]. Water Resources Research, 47(7): W07534.

    ZHANG K L, GONG Z, ZHAO K,, 2021. Experimental and numerical modeling of overhanging riverbank stability [J]. Journal of Geophysical Research: Earth Surface, 126(10): e2021JF006109.

    ZHANG Q, GONG Z, ZHANG C K,, 2016. Velocity and sediment surge: what do we see at times of very shallow water on intertidal mudflats? [J]. Continental Shelf Research, 113: 10-20.

    ZHAO K, GONG Z, ZHANG K L,, 2020. Laboratory experiments of bank collapse: the role of bank height and near-bank water depth [J]. Journal of Geophysical Research: Earth Surface, 125(5): e2019JF005281.

    RESEARCH PROGRESS ON BANK EROSION OF TIDAL CHANNEL UNDER COUPLED EFFECT OF SEEPAGE AND SURFACE FLOW

    GONG Zheng1, WANG Ke-Yu1, 2, ZHAO Kun1, ZHANG Kai-Li1

    (1. Jiangsu Key Laboratory of Coast Ocean Resources Development and Environment Security, Hohai University, Nanjing 210098, China; 2. PowerChina Northwest Engineering Co. Ltd., Xi’an 710065, China)

    The bank erosion of tidal channel can be divided into bank scouring by flow and bank collapse by gravity. Previous studies have shown that flow velocity, water pressure, soil dry density, and other factors would affect the bed scouring rate. However, the measurements of scouring rate of sidewall bank and related influencing factors need remain unsolved. Previous studies paid more attention to the influence of surface flow on bank collapse, but less on seepage and the coupled effect of seepage and surface flow on bank instability. Recent studies show that seepage play an important role in the destabilization of tidal channel bank. In the future, the difference and similarity of bank erosion under different dynamic actions of seepage and surface flow, and the destabilization mechanism of tidal channel bank under the coupled effect of the two forces need to be explored.

    tidal channel; scour rate; bank erosion; seepage

    * 國(guó)家杰出青年科學(xué)基金項(xiàng)目, 51925905號(hào)。龔 政, 博士生導(dǎo)師, 教授, E-mail: gongzheng@hhu.edu.cn

    王客予, 碩士研究生, E-mail: keyu_wang97@163.com

    2021-11-29,

    2022-01-21

    TV149.1

    10.11693/hyhz20211100302

    猜你喜歡
    影響實(shí)驗(yàn)研究
    記一次有趣的實(shí)驗(yàn)
    FMS與YBT相關(guān)性的實(shí)證研究
    是什么影響了滑動(dòng)摩擦力的大小
    遼代千人邑研究述論
    哪些顧慮影響擔(dān)當(dāng)?
    視錯(cuò)覺在平面設(shè)計(jì)中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    做個(gè)怪怪長(zhǎng)實(shí)驗(yàn)
    EMA伺服控制系統(tǒng)研究
    NO與NO2相互轉(zhuǎn)化實(shí)驗(yàn)的改進(jìn)
    實(shí)踐十號(hào)上的19項(xiàng)實(shí)驗(yàn)
    太空探索(2016年5期)2016-07-12 15:17:55
    亚洲av日韩精品久久久久久密| 成人特级黄色片久久久久久久| 欧美最新免费一区二区三区 | 久久久久久久亚洲中文字幕 | 窝窝影院91人妻| 欧美激情久久久久久爽电影| 成年女人看的毛片在线观看| 熟女电影av网| 又紧又爽又黄一区二区| 成年女人看的毛片在线观看| av中文乱码字幕在线| 三级国产精品欧美在线观看| 此物有八面人人有两片| 婷婷精品国产亚洲av在线| 成人特级黄色片久久久久久久| 欧美高清成人免费视频www| 久久久久久久亚洲中文字幕 | 欧美日韩一级在线毛片| 欧美成人一区二区免费高清观看| 女人十人毛片免费观看3o分钟| 亚洲aⅴ乱码一区二区在线播放| 午夜福利在线观看吧| 国产淫片久久久久久久久 | 日韩欧美三级三区| 在线观看一区二区三区| 欧美午夜高清在线| 国产欧美日韩精品亚洲av| 免费看a级黄色片| 国产精品影院久久| 精品一区二区三区人妻视频| 国产色婷婷99| 国产亚洲欧美98| 国产麻豆成人av免费视频| 日韩亚洲欧美综合| 一区二区三区免费毛片| 美女cb高潮喷水在线观看| 久久久久九九精品影院| 嫩草影院入口| 国产中年淑女户外野战色| 亚洲第一电影网av| 亚洲国产精品久久男人天堂| 国产伦一二天堂av在线观看| 欧美最新免费一区二区三区 | 狠狠狠狠99中文字幕| 国产探花极品一区二区| 又黄又粗又硬又大视频| 国产乱人伦免费视频| 91久久精品电影网| 亚洲精品亚洲一区二区| avwww免费| 亚洲av电影在线进入| 丝袜美腿在线中文| 国产淫片久久久久久久久 | 一级作爱视频免费观看| 亚洲avbb在线观看| 亚洲国产日韩欧美精品在线观看 | 在线观看午夜福利视频| 亚洲七黄色美女视频| 久久久久国内视频| 嫩草影院精品99| 性色avwww在线观看| x7x7x7水蜜桃| 成人18禁在线播放| 亚洲午夜理论影院| 国产精品久久久久久久电影 | 成年免费大片在线观看| 日本一本二区三区精品| 国产探花极品一区二区| 一区二区三区高清视频在线| 婷婷精品国产亚洲av| 国产亚洲精品一区二区www| 又紧又爽又黄一区二区| 91字幕亚洲| 国产av在哪里看| 久久精品亚洲精品国产色婷小说| 国产熟女xx| 日韩中文字幕欧美一区二区| 国产精品久久久久久亚洲av鲁大| 亚洲人成网站在线播放欧美日韩| 搞女人的毛片| 国产午夜福利久久久久久| 成人三级黄色视频| 长腿黑丝高跟| 亚洲国产精品成人综合色| 三级毛片av免费| 国产精品久久久人人做人人爽| 一区福利在线观看| 在线免费观看的www视频| 淫秽高清视频在线观看| 国产三级黄色录像| 欧美一区二区亚洲| 久久精品91无色码中文字幕| 别揉我奶头~嗯~啊~动态视频| 色视频www国产| 啦啦啦观看免费观看视频高清| 又黄又粗又硬又大视频| 国产老妇女一区| 99久久久亚洲精品蜜臀av| 国产精品av视频在线免费观看| 日本 av在线| 极品教师在线免费播放| 91在线精品国自产拍蜜月 | 日韩大尺度精品在线看网址| 久久精品国产亚洲av涩爱 | 69人妻影院| 熟女少妇亚洲综合色aaa.| 亚洲av二区三区四区| 国产91精品成人一区二区三区| 黄片大片在线免费观看| 搡老熟女国产l中国老女人| 精华霜和精华液先用哪个| 亚洲精品影视一区二区三区av| 级片在线观看| 中文字幕av成人在线电影| 欧美日韩中文字幕国产精品一区二区三区| 成熟少妇高潮喷水视频| 亚洲专区国产一区二区| 亚洲av日韩精品久久久久久密| 国产97色在线日韩免费| 亚洲成人免费电影在线观看| a级毛片a级免费在线| 国产精品亚洲av一区麻豆| 99久久九九国产精品国产免费| 免费av毛片视频| 成年人黄色毛片网站| 日日夜夜操网爽| 女人高潮潮喷娇喘18禁视频| 最近最新免费中文字幕在线| 午夜精品一区二区三区免费看| 欧美日韩综合久久久久久 | 18禁裸乳无遮挡免费网站照片| xxx96com| 国产男靠女视频免费网站| 51午夜福利影视在线观看| 国产精品亚洲av一区麻豆| 午夜免费男女啪啪视频观看 | h日本视频在线播放| 成人特级黄色片久久久久久久| 男女午夜视频在线观看| 一二三四社区在线视频社区8| 亚洲精品影视一区二区三区av| 欧美黄色淫秽网站| 一级毛片高清免费大全| 熟女电影av网| 草草在线视频免费看| 日本免费a在线| 亚洲人与动物交配视频| 国产成人a区在线观看| 免费无遮挡裸体视频| 国产真实伦视频高清在线观看 | 国产精品99久久99久久久不卡| 亚洲国产精品999在线| 久久亚洲真实| 精品久久久久久久人妻蜜臀av| 91av网一区二区| 日本 欧美在线| 精品国内亚洲2022精品成人| 男女下面进入的视频免费午夜| 给我免费播放毛片高清在线观看| 亚洲国产欧美网| 中文字幕人成人乱码亚洲影| 久99久视频精品免费| 久久精品国产亚洲av涩爱 | 免费看十八禁软件| 免费av观看视频| 精品人妻偷拍中文字幕| 亚洲国产欧美人成| 欧美中文综合在线视频| 岛国在线观看网站| 亚洲片人在线观看| 精品无人区乱码1区二区| 少妇熟女aⅴ在线视频| 美女黄网站色视频| 国产精品久久久久久久久免 | 亚洲av成人不卡在线观看播放网| 一边摸一边抽搐一进一小说| 国产午夜精品论理片| 国产欧美日韩一区二区三| 精品熟女少妇八av免费久了| 亚洲精品一卡2卡三卡4卡5卡| 长腿黑丝高跟| 欧美国产日韩亚洲一区| 男女午夜视频在线观看| 欧美色视频一区免费| 又黄又粗又硬又大视频| 国产一区二区在线av高清观看| 99久久综合精品五月天人人| 精品电影一区二区在线| 国产成人影院久久av| 免费av观看视频| 草草在线视频免费看| 国产精品美女特级片免费视频播放器| 国产乱人伦免费视频| 亚洲精品国产精品久久久不卡| 99久久综合精品五月天人人| 亚洲精品亚洲一区二区| 国产亚洲精品久久久com| 国产伦在线观看视频一区| 国产精品,欧美在线| 99久久精品国产亚洲精品| 黄色日韩在线| 天美传媒精品一区二区| 成人欧美大片| 一进一出好大好爽视频| www国产在线视频色| 欧美三级亚洲精品| 男女之事视频高清在线观看| a级一级毛片免费在线观看| 午夜a级毛片| 女警被强在线播放| 亚洲,欧美精品.| 可以在线观看的亚洲视频| 成年版毛片免费区| 偷拍熟女少妇极品色| 级片在线观看| 亚洲国产精品sss在线观看| 桃红色精品国产亚洲av| 亚洲色图av天堂| 琪琪午夜伦伦电影理论片6080| 国产成人啪精品午夜网站| 精品电影一区二区在线| 一区二区三区免费毛片| 日本熟妇午夜| 亚洲国产欧洲综合997久久,| 精品久久久久久久久久久久久| 婷婷精品国产亚洲av在线| 欧美xxxx黑人xx丫x性爽| 久久精品国产综合久久久| 他把我摸到了高潮在线观看| 欧美成人免费av一区二区三区| 国产高潮美女av| 母亲3免费完整高清在线观看| 美女被艹到高潮喷水动态| 亚洲片人在线观看| 中出人妻视频一区二区| 婷婷精品国产亚洲av在线| 久久精品91无色码中文字幕| 757午夜福利合集在线观看| 麻豆国产97在线/欧美| 欧美国产日韩亚洲一区| 亚洲精品一区av在线观看| 国产亚洲精品久久久com| www日本在线高清视频| 成年女人永久免费观看视频| 国产又黄又爽又无遮挡在线| 男女下面进入的视频免费午夜| 无人区码免费观看不卡| 日本黄大片高清| 91久久精品国产一区二区成人 | 美女高潮喷水抽搐中文字幕| 特大巨黑吊av在线直播| 婷婷精品国产亚洲av在线| 亚洲五月婷婷丁香| 99久久九九国产精品国产免费| 久久久国产精品麻豆| 国产aⅴ精品一区二区三区波| 国产真实伦视频高清在线观看 | 在线观看午夜福利视频| 又黄又爽又免费观看的视频| 国产精品女同一区二区软件 | 亚洲精品乱码久久久v下载方式 | 国产成人av激情在线播放| 欧美在线黄色| 美女 人体艺术 gogo| 国产视频一区二区在线看| 免费观看人在逋| 成人欧美大片| 免费高清视频大片| 中文在线观看免费www的网站| 欧美日本亚洲视频在线播放| 非洲黑人性xxxx精品又粗又长| 中文亚洲av片在线观看爽| 两个人视频免费观看高清| 少妇的逼水好多| 99热只有精品国产| av黄色大香蕉| 日韩欧美一区二区三区在线观看| 亚洲av不卡在线观看| 好男人电影高清在线观看| 高清毛片免费观看视频网站| 成年版毛片免费区| 欧美在线一区亚洲| 在线观看舔阴道视频| 天堂√8在线中文| 日本三级黄在线观看| 久久国产乱子伦精品免费另类| 麻豆国产av国片精品| 女同久久另类99精品国产91| 色精品久久人妻99蜜桃| 免费无遮挡裸体视频| 国内毛片毛片毛片毛片毛片| 在线观看舔阴道视频| 19禁男女啪啪无遮挡网站| 国产aⅴ精品一区二区三区波| 久久精品国产自在天天线| 精品无人区乱码1区二区| 三级国产精品欧美在线观看| 亚洲成av人片免费观看| 亚洲成人免费电影在线观看| 十八禁人妻一区二区| 在线十欧美十亚洲十日本专区| 色综合站精品国产| 免费人成在线观看视频色| 亚洲激情在线av| 精品99又大又爽又粗少妇毛片 | 精品电影一区二区在线| 久99久视频精品免费| 精品一区二区三区视频在线观看免费| 18禁裸乳无遮挡免费网站照片| 欧美激情久久久久久爽电影| 亚洲av熟女| 国产成+人综合+亚洲专区| 亚洲无线在线观看| 在线国产一区二区在线| 久久久久久久午夜电影| 一区二区三区国产精品乱码| 九九热线精品视视频播放| 精品国内亚洲2022精品成人| 99在线人妻在线中文字幕| xxx96com| 亚洲激情在线av| 日韩欧美在线二视频| 亚洲精品成人久久久久久| 欧美色欧美亚洲另类二区| 中文字幕av成人在线电影| 麻豆国产av国片精品| 日韩中文字幕欧美一区二区| 久久性视频一级片| 国产午夜福利久久久久久| 最后的刺客免费高清国语| 国产精品99久久99久久久不卡| 18禁美女被吸乳视频| 亚洲国产日韩欧美精品在线观看 | 婷婷亚洲欧美| 夜夜爽天天搞| 亚洲色图av天堂| 两个人视频免费观看高清| 成人永久免费在线观看视频| 国产黄色小视频在线观看| 好看av亚洲va欧美ⅴa在| 国产精品爽爽va在线观看网站| 可以在线观看毛片的网站| 中文字幕人妻丝袜一区二区| 精品国内亚洲2022精品成人| 国产一区二区亚洲精品在线观看| 成人精品一区二区免费| 欧美性猛交黑人性爽| 亚洲av中文字字幕乱码综合| 亚洲精品粉嫩美女一区| 丁香六月欧美| 九九热线精品视视频播放| 国产激情偷乱视频一区二区| a级毛片a级免费在线| 午夜视频国产福利| 99精品欧美一区二区三区四区| 午夜激情福利司机影院| 好看av亚洲va欧美ⅴa在| 亚洲精品456在线播放app | 亚洲最大成人手机在线| 一本综合久久免费| 老汉色av国产亚洲站长工具| 欧美色欧美亚洲另类二区| 国产毛片a区久久久久| 国产91精品成人一区二区三区| 成人性生交大片免费视频hd| 欧美日韩瑟瑟在线播放| 日本精品一区二区三区蜜桃| 久久精品国产99精品国产亚洲性色| 国产极品精品免费视频能看的| 观看美女的网站| 中文字幕久久专区| 男人舔女人下体高潮全视频| 国产精品美女特级片免费视频播放器| 老熟妇乱子伦视频在线观看| 国产亚洲欧美在线一区二区| 久久伊人香网站| 此物有八面人人有两片| 欧美性猛交╳xxx乱大交人| 欧美国产日韩亚洲一区| 亚洲av成人不卡在线观看播放网| 亚洲精品久久国产高清桃花| 国产伦一二天堂av在线观看| 亚洲av美国av| 欧美极品一区二区三区四区| 一级毛片高清免费大全| 日本黄大片高清| 欧美绝顶高潮抽搐喷水| 少妇人妻一区二区三区视频| 国产亚洲av嫩草精品影院| 成人国产综合亚洲| 露出奶头的视频| 一本精品99久久精品77| 婷婷精品国产亚洲av在线| 午夜福利在线观看吧| 好看av亚洲va欧美ⅴa在| 高清日韩中文字幕在线| 一进一出抽搐动态| 岛国在线观看网站| 一二三四社区在线视频社区8| 一卡2卡三卡四卡精品乱码亚洲| 悠悠久久av| 欧美成人a在线观看| 中文资源天堂在线| 88av欧美| 亚洲久久久久久中文字幕| 精品人妻1区二区| 久久九九热精品免费| 欧美不卡视频在线免费观看| 噜噜噜噜噜久久久久久91| 黑人欧美特级aaaaaa片| eeuss影院久久| 国产成人av教育| 精品国产亚洲在线| 夜夜爽天天搞| 国产探花极品一区二区| 久久久久久国产a免费观看| 最新美女视频免费是黄的| 午夜福利欧美成人| 亚洲人成电影免费在线| 日韩精品中文字幕看吧| 麻豆国产97在线/欧美| 夜夜躁狠狠躁天天躁| 热99在线观看视频| 又爽又黄无遮挡网站| 午夜影院日韩av| 最新在线观看一区二区三区| 亚洲精品亚洲一区二区| 国产真实乱freesex| 不卡一级毛片| 国产亚洲av嫩草精品影院| 欧美精品啪啪一区二区三区| 舔av片在线| 午夜视频国产福利| www.999成人在线观看| 日本在线视频免费播放| 热99在线观看视频| 亚洲国产中文字幕在线视频| 3wmmmm亚洲av在线观看| 色视频www国产| 亚洲精品亚洲一区二区| 在线观看日韩欧美| 色哟哟哟哟哟哟| 男插女下体视频免费在线播放| 国产老妇女一区| 嫩草影院入口| 脱女人内裤的视频| 五月玫瑰六月丁香| 黑人欧美特级aaaaaa片| 欧美黑人欧美精品刺激| 久久这里只有精品中国| 国产精品1区2区在线观看.| 超碰av人人做人人爽久久 | 国产aⅴ精品一区二区三区波| 三级国产精品欧美在线观看| 哪里可以看免费的av片| 成人性生交大片免费视频hd| 中文字幕av在线有码专区| 国产av在哪里看| 搡老妇女老女人老熟妇| 亚洲人成网站在线播| 一级黄色大片毛片| 岛国视频午夜一区免费看| 好男人在线观看高清免费视频| 怎么达到女性高潮| 成人国产一区最新在线观看| 中文字幕熟女人妻在线| 久久久国产成人精品二区| 岛国在线观看网站| 国产蜜桃级精品一区二区三区| 午夜福利视频1000在线观看| 99精品在免费线老司机午夜| 波多野结衣高清无吗| 午夜亚洲福利在线播放| 欧美日韩亚洲国产一区二区在线观看| 老熟妇乱子伦视频在线观看| 91麻豆精品激情在线观看国产| 午夜精品在线福利| 国产伦精品一区二区三区四那| 午夜免费成人在线视频| 亚洲欧美日韩卡通动漫| 中文字幕人妻熟人妻熟丝袜美 | 女人十人毛片免费观看3o分钟| 国产av一区在线观看免费| 亚洲一区二区三区色噜噜| 淫秽高清视频在线观看| 亚洲欧美精品综合久久99| 精品不卡国产一区二区三区| 深夜精品福利| 中文字幕精品亚洲无线码一区| 欧美+日韩+精品| 国产伦精品一区二区三区视频9 | 婷婷丁香在线五月| 在线观看免费视频日本深夜| 中文字幕人妻丝袜一区二区| 老司机午夜十八禁免费视频| 亚洲最大成人中文| 丰满人妻一区二区三区视频av | 制服人妻中文乱码| 亚洲专区中文字幕在线| 欧美日韩综合久久久久久 | 在线免费观看不下载黄p国产 | 欧美日本亚洲视频在线播放| 国产三级黄色录像| 国产亚洲欧美98| 可以在线观看毛片的网站| 人人妻,人人澡人人爽秒播| 欧美zozozo另类| 国产在线精品亚洲第一网站| 99国产极品粉嫩在线观看| 母亲3免费完整高清在线观看| 高潮久久久久久久久久久不卡| 波多野结衣高清无吗| 天天躁日日操中文字幕| av专区在线播放| 一a级毛片在线观看| 国产淫片久久久久久久久 | 99在线视频只有这里精品首页| 国产精品亚洲美女久久久| 香蕉av资源在线| 夜夜夜夜夜久久久久| 在线观看免费视频日本深夜| 一个人免费在线观看的高清视频| 法律面前人人平等表现在哪些方面| 国产高潮美女av| 亚洲欧美日韩卡通动漫| 99热精品在线国产| 欧美bdsm另类| 精品熟女少妇八av免费久了| 又黄又爽又免费观看的视频| 亚洲熟妇熟女久久| 亚洲国产中文字幕在线视频| 亚洲专区中文字幕在线| 夜夜爽天天搞| 国产亚洲精品综合一区在线观看| 男女做爰动态图高潮gif福利片| 日韩有码中文字幕| av天堂中文字幕网| 俄罗斯特黄特色一大片| 国产又黄又爽又无遮挡在线| 无遮挡黄片免费观看| 黄片大片在线免费观看| 特级一级黄色大片| 亚洲熟妇熟女久久| 日本三级黄在线观看| 免费av毛片视频| 啦啦啦观看免费观看视频高清| 国产国拍精品亚洲av在线观看 | 精品国产亚洲在线| 最近最新中文字幕大全免费视频| 国产欧美日韩精品一区二区| 99久久精品热视频| 久久久国产成人精品二区| 三级国产精品欧美在线观看| 色播亚洲综合网| 欧美日本视频| АⅤ资源中文在线天堂| 国产精品 欧美亚洲| 2021天堂中文幕一二区在线观| 精品欧美国产一区二区三| 免费无遮挡裸体视频| 99国产综合亚洲精品| 久久久久久久亚洲中文字幕 | 日韩人妻高清精品专区| 俺也久久电影网| 日韩精品中文字幕看吧| 嫩草影院精品99| 日韩大尺度精品在线看网址| АⅤ资源中文在线天堂| 啦啦啦韩国在线观看视频| 深夜精品福利| 国产精品永久免费网站| 亚洲成人中文字幕在线播放| 在线免费观看不下载黄p国产 | 日本a在线网址| 在线国产一区二区在线| 一二三四社区在线视频社区8| 99久久成人亚洲精品观看| 亚洲成人中文字幕在线播放| 老司机午夜十八禁免费视频| 悠悠久久av| 午夜福利视频1000在线观看| 国产精品98久久久久久宅男小说| 午夜亚洲福利在线播放| 99久久99久久久精品蜜桃| 有码 亚洲区| 午夜精品一区二区三区免费看| 精品国产美女av久久久久小说| 成人亚洲精品av一区二区| 色在线成人网| 亚洲欧美激情综合另类| 国产乱人视频| 欧美日韩一级在线毛片| 国产极品精品免费视频能看的| 日韩高清综合在线| 国产黄片美女视频| 丝袜美腿在线中文| 成人高潮视频无遮挡免费网站| 国产精品女同一区二区软件 | 他把我摸到了高潮在线观看| 亚洲avbb在线观看| 午夜福利欧美成人| 91麻豆av在线| 美女免费视频网站| 又爽又黄无遮挡网站| 村上凉子中文字幕在线| xxxwww97欧美| 波多野结衣高清无吗| 国产成人av教育| 亚洲 国产 在线| 日韩有码中文字幕| 日韩免费av在线播放| 亚洲国产色片| 国内少妇人妻偷人精品xxx网站| 久9热在线精品视频| 亚洲成人久久性|