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

    曲流串溝型江心洲形成機(jī)制與演化探討
    ——以現(xiàn)代松花江為例

    2022-11-16 02:24:44張憲國王涵巍張濤段冬平林承焰黃鑫
    沉積學(xué)報 2022年5期
    關(guān)鍵詞:江心洲流河河道

    張憲國,王涵巍,2,張濤,段冬平,林承焰,黃鑫

    1.中國石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院,山東青島 266580

    2.國家知識產(chǎn)權(quán)局專利局專利審查協(xié)作河南中心,鄭州 450046

    3.山東科技大學(xué)地球科學(xué)與工程學(xué)院,山東青島 266590

    4.中海石油(中國)有限公司上海分公司研究院,上海 200035

    0 引言

    江心洲是一種在河道分汊處出露于水面的沉積體,在辮狀河、曲流河、網(wǎng)狀河等不同河型中普遍發(fā)育,不僅是現(xiàn)代河流治理、防災(zāi)減災(zāi)、港口建設(shè)和生態(tài)保護(hù)的重要考慮因素,也是地下古河流沉積油氣儲層表征的重點砂體成因類型。但是,在水利學(xué)和油氣儲層沉積學(xué)等不同研究領(lǐng)域和行業(yè)中對江心洲的成因與演化認(rèn)識存在差異。在河流沉積學(xué)中,河流上游和中游發(fā)育的江心洲常被視為大型辮狀河心灘[1]。在水利學(xué)研究中,江心洲是游蕩河和分汊河中影響河道演變與穩(wěn)定的重要砂體類型[2]。河流中心區(qū)域沙壩沉積、橫向沙壩轉(zhuǎn)化、多壩分割、溯源侵蝕水道或者幾種作用的迭加被認(rèn)為是江心洲形成的主要機(jī)制[1-4]?,F(xiàn)代沉積觀察顯示,無論是辮狀河還是曲流河等不同河型中,江心洲都是一種普遍存在的砂體類型[5],但形成這些江心洲的地貌、物源、水動力條件等不盡相同[6-9]。曲流串溝型江心洲是平原曲流河中常見的一種沉積單元類型,多發(fā)育在高彎度曲流河中,在密蘇里河、科羅拉多河以及長江干流等的曲流河段都存在,但已有研究主要是從河流治理和工程角度對其形態(tài)和汊道流量等的分析[10-14],對江心洲形成的物理過程和控制因素研究缺失,對復(fù)雜水流、泥沙輸送、陸岸作用等問題缺少精細(xì)和量化的認(rèn)識。對于江心洲內(nèi)部結(jié)構(gòu)、巖相組合、沉積模式等尚無系統(tǒng)性的成果報道。本研究在松花江哈爾濱水文站附近選取江心洲發(fā)育特征典型、衛(wèi)星照片等資料相對齊全、沒有人工堤岸干擾的河段開展研究,利用不同時期衛(wèi)星圖片分析和沉積數(shù)值模擬,探究曲流串溝型江心洲的形成機(jī)制與演化,為現(xiàn)代河流治理和地下古河道砂體表征與預(yù)測提供理論依據(jù)。

    1 研究區(qū)概況

    松花江位于中國東北地區(qū)北部,流域面積大,發(fā)育范圍在41°42′~51°38′N,119°52′~132°31′E,松花江全長1 927 km,流域面積55.68×104km2,松花江干流,從吉林省松原市嫩江注入,至匯入黑龍江的河口,俗稱東流松花江。根據(jù)東流松花江的地形及河道特性,可分為上、中、下三段。三汊河至哈爾濱市為上段,也是研究區(qū)所在江段,其中三汊河至下岱吉坡降較緩,為0.022‰,下岱吉至謝家屯江道坡降0.06‰,謝家屯附近至哈爾濱坡降0.052%,屬于平原低緩地貌背景特征,一般洪水即可淹沒凸岸點壩以及平水期出露在河道中的江心洲。點壩的灘面上溝槽發(fā)育,在垂向上從河床向上巖性依次為細(xì)砂—中砂、粉砂或細(xì)砂、黏土或砂質(zhì)黏土(河岸上層0.5~1.5 m)。

    該區(qū)域的松花江哈爾濱水文站近50年的平均年徑流量402.3×108m3。松花江徑流量月度差異大,季風(fēng)性氣候特征明顯,汛期較短,與雨季緊密關(guān)聯(lián),徑流量在8月份最大,2月最小,統(tǒng)計發(fā)現(xiàn)汛期8、9兩月合計流量可占年流量的30%~51%,即兩月平均流速介于2 470~3 930 m3/s,統(tǒng)計1981—1999 年松花江干流平均兩年發(fā)生一次洪水,一次洪水歷時10余天,最大洪水峰值一般在6 000~12 000 m3/s[15-18]。本次研究河段位于松花江上段的大慶肇源縣,河段長約14 km,寬度在200~700 m,下游端距離哈爾濱約120 km,整個研究區(qū)由兩個河曲構(gòu)成,河流彎曲系數(shù)在2.1,屬于高彎度曲流河(圖1)。研究區(qū)曲流河道內(nèi)部發(fā)育大型江心洲,發(fā)育位置在曲流河的河曲位置,露出水面部分的江心洲長軸和短軸長度分別約為1 400 m和600 m。本研究開展基于沉積水動力的沉積過程數(shù)值模擬,認(rèn)識該江心洲的形成過程和水動力機(jī)制,并通過不同時間衛(wèi)星照片揭示的沉積形態(tài)進(jìn)行比較分析。

    圖1 研究區(qū)河段地理位置Fig.1 Geographic location of study river section

    2 研究方法

    本研究采用基于水動力學(xué)的沉積數(shù)值模擬方法,再現(xiàn)江心洲的形成過程,并定量分析形成過程中水流速度等水動力參數(shù)在空間和時間上的變化。

    2.1 沉積數(shù)值模擬方法

    本研究采用的基于泥沙水動力學(xué)的沉積模擬方法是以Navier-Stokes方程為基礎(chǔ),通過物質(zhì)平衡原理和沉積物搬運、沉積及侵蝕方程求解,實現(xiàn)沉積物搬運與沉積過程的三維模擬。該方法被廣泛應(yīng)用于三角洲地貌的控制機(jī)理、三角洲分流河道分叉、辮狀河心灘的形成過程與機(jī)理等研究[19-23]。

    認(rèn)識流場預(yù)測水流流動是分析沉積過程的重要依據(jù)和基礎(chǔ),需要建立在對水流動力、沉積與剝蝕等物理過程的詳細(xì)描述之上,綜合考慮水動力—泥沙輸運—地貌變化這一完整的物理過程,要求上述過程的數(shù)值模型完整地處理水流對地貌的改造以及地貌改變后對水流的動態(tài)反饋過程[24-26]。

    在本次對曲流河江心洲的模擬涉及的主要水動力方程包括泥沙搬運與沉積方程、非黏性物質(zhì)沉積與侵蝕方程和泥沙沉降速度方程:

    (1)泥沙搬運與沉積方程

    泥沙濃度、水流速度、擴(kuò)散系數(shù)以及沉積物顆粒的沉降速度等因素是控制沉積和搬運的重要參數(shù),這些參數(shù)對沉積物的搬運與沉積過程的控制可以通過(式1)進(jìn)行表征。

    根據(jù)模擬對象的沉積特征,本研究選用Delft 3D模擬軟件進(jìn)行曲流河江心洲形成過程的模擬,圖2為一個計算步長中的模擬環(huán)節(jié)和流程。

    圖2 沉積過程模擬計算流程Fig.2 Workflow of sedimentary process simulation

    2.2 沉積模擬條件

    能夠真實反映沉積環(huán)境和沉積水動力特征的模擬參數(shù)是沉積數(shù)值模擬結(jié)果可靠的基礎(chǔ),本研究的沉積模擬參數(shù)來源于研究區(qū)實際地貌和水文資料。

    2.2.1 模型網(wǎng)格與初始地貌

    研究區(qū)邊界數(shù)據(jù)來源于谷歌地球提供的衛(wèi)星數(shù)據(jù),研究河段模型構(gòu)建采用了多次加密和修改網(wǎng)格的方法,在衛(wèi)星圖片上勾勒出河段計算邊界后,先粗略劃分網(wǎng)格,再轉(zhuǎn)為緊密網(wǎng)格,根據(jù)需要對余弦值大于0.02的網(wǎng)格進(jìn)行正交化調(diào)整,改善網(wǎng)格的正交性,調(diào)整網(wǎng)格數(shù)量和網(wǎng)格分布。在網(wǎng)格化基礎(chǔ)上進(jìn)行地貌模型的構(gòu)建,地貌高程數(shù)據(jù)來源于ASF Data Search網(wǎng)站公開的DEM數(shù)據(jù)庫。

    2.2.2 模擬條件與參數(shù)

    讀書的“浪漫”倘真至于如此,則據(jù)我猜想,蓋屬于情侶同讀沖淡平和、溫情脈脈的“性靈小品”,而非獨自頌讀疾痛慘怛、憂憤難譴的《離騷》!

    綜合考慮模擬的精度要求以及穩(wěn)定性和收斂性,模擬時間步長取0.2 min,模擬時間為兩年,輸出時間步長為2 000 min。

    邊界條件設(shè)置是影響模擬的重要基礎(chǔ)參數(shù),本研究流量序列為輸入邊界,水位序列為流出邊界。河流的沉積建造不是隨著時間均勻沉積或侵蝕的,沉積和侵蝕事件主要集中在全年短暫的幾個月[27],對于本研究區(qū)選取的松花江河段來說,河流地貌變化多集中于雨季(8、9月),且主要的沉積和侵蝕也是發(fā)生在這個季節(jié)。因此,為了降低模擬總耗時,在模擬過程中去掉了非汛期月份,直接將高水位期及其間的洪水期在時間上串聯(lián)進(jìn)行模擬??紤]模擬參數(shù)與實際河流的相似性,設(shè)置高水位期與洪水期兩個周期性流量,高水位期取值參照8、9 月平均流量值4 000 m3/s,洪水期(設(shè)置為12天)取12 000 m3/s,兩周期相對長度約為4∶1,這些數(shù)值均在該河段歷史水文記錄數(shù)據(jù)范圍內(nèi),為了使模擬效果更加顯著,取值采用了上述歷史數(shù)據(jù)范圍內(nèi)的偏上限數(shù)值。哈爾濱水文站記錄數(shù)據(jù)表明大洪水期與汛期水位相差約4 m高程,綜合考慮到模型采用DEM(數(shù)字高程模型)基準(zhǔn)參考系存在差異,最終設(shè)置兩個流量周期的下游邊界水位對應(yīng)概化值為132 m和136 m??紤]模型計算效率,簡化了實際床砂級配,突出級配主要特點,將床砂可動層級配設(shè)置為黏土和三個粒徑級別的砂(粒徑分別為64μm、200μm和500μm),單一厚度均設(shè)置為5 m,床砂混合均勻,按照上述設(shè)置,可侵蝕的河床厚度為20 m。

    渦黏系數(shù)與水深的關(guān)系是拋物線形態(tài),但在大水流過程中,動量方程中的紊動輸運項影響小,為簡化計算,將渦黏系數(shù)近似為線性分布。由于本研究著重關(guān)注河床地貌的演化歷程及沉積結(jié)構(gòu)特征,而地貌演變的時間尺度總遠(yuǎn)大于水動力學(xué)研究的時間尺度,本研究使用地貌加速因子來加速地貌過程(表1)。

    表1 沉積過程模擬的主要參數(shù)Table 1 Main parameters of sedimentary process simulation

    3 模擬實驗結(jié)果及分析

    通過模擬恢復(fù)曲流串溝發(fā)育、汊道廢棄等現(xiàn)象,認(rèn)識曲流河演化過程,揭示串溝截彎機(jī)制,結(jié)合現(xiàn)代沉積,認(rèn)識串溝型江心洲地貌演變及水動力特征。

    3.1 串溝截彎形成機(jī)制

    串溝型截彎過程控制著江心洲及周邊汊道演化。第一組模擬結(jié)果顯示(圖3),在平水期,也就是時間步長(time step,以下簡寫作TS)140 之前,由于二次流和橫向斜坡效應(yīng)綜合作用,在彎曲河道中的水流發(fā)生澭水現(xiàn)象,在河道橫截面上形成凹岸高凸岸低的水面斜坡,并將沉積物不斷地從崩塌蝕退的凹岸搬運至凸岸堆積。這一過程的直接結(jié)果就是形成了河道的側(cè)向遷移和點壩加積,形成不斷側(cè)向生長的點壩單元,同時,河道的彎曲度增加而彎曲半徑減小。模擬結(jié)果顯示,在橫截面N65 處凸岸繼續(xù)推進(jìn)12 m,速率約0.1 m/TS,凹岸向后侵蝕26 m;在橫截面N15 處,凸岸推進(jìn)45 m,凹岸蝕退31 m(圖4)??紤]到水流的泥沙含量較低,而水動力強(qiáng),大量沉積物遭受侵蝕后并未就近沉積而是繼續(xù)搬運到河流下游,故而不足以維持與之均衡的凸岸加積速率。在側(cè)向加積過程中,點壩增生軌跡在平面上形成灘脊和流槽,相間分布于點壩頂部,這種現(xiàn)象常見于現(xiàn)代河流沉積。為便于描述,在圖5中對主要的分汊口及汊道以及幾個不同的角度參數(shù)進(jìn)行了編號。

    圖4 N65 截面堤岸側(cè)向遷移(剖面位置見圖3a1)Fig.4 Lateral migration of bank in profile N65(location in Fig.3a1)

    圖5 曲流河分汊口地貌Fig.5 Geomorphology of bifurcation point in meandering river

    進(jìn)入第一個洪水期(140 TS)后,高峰流量引起河流水位上升,水流漫溢到河岸兩側(cè)區(qū)域,其中,在凸岸的一側(cè),漫溢的洪水在其表面沉積細(xì)粒沉積物,形成覆蓋在早期點壩之上的泛濫平原沉積。模擬結(jié)果顯示,當(dāng)高峰流量通過時,水流速增大且水位上升,水流的主流線發(fā)生偏移,在河曲(凸岸點壩)上游一端,深泓線向凸岸一側(cè)偏移,最大應(yīng)力值為9 N/m2。相應(yīng)地,在河曲上游端的水流高速區(qū)也集中分布于凸岸附近,達(dá)到1.7~1.9 m/s。隨著洪水期水位的急劇上升,高速水流對點壩的改造加劇(圖3b1,c1),水流不斷漫溢出河道,而漫溢出的洪水首先進(jìn)入點壩表面相對低洼的弧形流槽中,這些點壩多期側(cè)積體之間自然形成的流槽在漫溢洪水的沖刷下不斷加深和拓寬,最終形成汊道#2(圖3a2)。隨著洪水持續(xù)不均衡侵蝕沖刷,在研究區(qū)兩個河灣的點壩內(nèi)側(cè)分別發(fā)育新的優(yōu)勢流路,形成汊道#3和#4,最高應(yīng)力集中區(qū)由原來的汊道#1轉(zhuǎn)移至汊道#3,水流快速下切侵蝕,將沉積物向下游搬運(圖3b2,c2)。在這一過程中產(chǎn)生的汊道#2寬約80 m,規(guī)模比原河道小,汊道#4 平面上呈蜿蜒狀,在平水期分流能力顯著弱于原河道。

    洪峰過后(200 TS),洪水消退,河流再次進(jìn)入平水期,水動力整體減弱,無論是最大應(yīng)力還是最大流速,都較洪水期顯著減弱,彎曲段水位回落(圖3b2,c2)。在200~540 TS,曲流河進(jìn)入緩慢演化時期,洪水期形成的河道形態(tài)基本穩(wěn)定,處于近沉積平衡的狀態(tài)。550 TS 之后第二洪水期開始,沉積平衡再次打破,上一期洪水形成的汊道被進(jìn)一步改造。#2 和#3 汊道逐漸加深和拓寬,從最大應(yīng)力分布和流速分布特征上看,均取代了舊河道主導(dǎo)地位(圖3a3~c3)。經(jīng)過了上述兩期洪水期與平水期的交替作用,到了本次模擬的末期(790 TS),上游河曲段最初的河床開始淤塞抬高,水流能量進(jìn)一步減弱,有逐漸廢棄的趨勢。下游河曲區(qū)域,#1、#4汊道環(huán)繞的江心洲逐漸向凹岸靠近。

    圖3 曲流河演化中的參數(shù)變化模擬結(jié)果(a1~a4)河床底形高度;(b1~b4)水流速度;(c1~c4)河床剪切力;(d)流量變化Fig.3 Simulated changes in main parameters in meandering river evolution

    利用遙感數(shù)據(jù)提供的不同時間河流形態(tài)圖像對上述模擬結(jié)果進(jìn)行對照。1994—2017年松花江研究區(qū)段上河曲段衛(wèi)星圖像展示了串溝截彎的過程(圖6)。1996年之前,研究河段在凸岸點壩上不同期次點壩側(cè)積單元之間發(fā)育相對低洼的流槽,這些流槽沿著相鄰點壩側(cè)積體的邊界發(fā)育,在該河段的點壩西側(cè)發(fā)育一條貫穿整個河曲的串溝(圖6a),寬度20~100 m;在2004 年12 月的衛(wèi)星圖像上看到,串溝繼續(xù)加寬,流徑更為順直,原主干河道寬度有變小的趨勢,與數(shù)值模擬結(jié)果吻合。2015 年8 月的衛(wèi)星圖像顯示,西側(cè)串溝已經(jīng)發(fā)育為新的河道,完成了截彎取直的過程,并形成長約2.4 km、寬260~500 m 的串溝型江心洲。到2019 年5 月,舊河道出入口兩端出現(xiàn)塞狀壩,向河道內(nèi)部延伸,舊河道寬度進(jìn)一步變窄,大部分不足百米,有逐漸萎縮廢棄的趨勢。大量現(xiàn)代河流的廢棄實例顯示,截彎末端可以在15~20年內(nèi)被阻塞,盡管河道壅塞、冰塞等都可以使河道逐漸淤塞,但廢棄段完全被填滿是一個漫長的過程,往往需要數(shù)十年甚至百年時間[28]。

    圖6 松花江河彎串溝型截彎過程Fig.6 Process of ditch-scouring cutoff in Songhua River

    3.2 汊道的演化

    在曲流河中,汊道形成和演化與江心洲的形成和演化是同一過程,汊道的穩(wěn)定決定了江心洲的發(fā)育,因此,理解分汊河道的演化對認(rèn)識曲流江心洲有重要意義。汊道的影響因素眾多[29-31],如上游曲率變化、新舊河道相對長度與坡度等。Burgeet al.[32]認(rèn)為主次河道若要同時保持活躍,與主干河道并存很多年,維持必要條件是主槽分流比要小于分沙比,或者說支流分沙比要小于分流比。主河道分流量占上游來水的總流量比例越大,意味著主河道對應(yīng)輸砂量占總量的比例也要隨著增大,才能維持弱水動力的新汊道不被泥沙淤填,簡單概括即“大水大砂,小水小砂”。可見,汊道的演化既受制于水動力狀況,也與地貌特征相關(guān)。本文通過模擬結(jié)果的參數(shù)統(tǒng)計,探討汊道維持與淤堵的可能性。

    串溝形成的汊道是切穿點壩而成的,其彎度要明顯小于主流河道,所以與主河道相比,汊道具有長度小、坡度大的特點,故而對沉積物搬運的能力更強(qiáng)。原河道與新汊道的長度比稱為相對梯度優(yōu)勢。從各汊道梯度優(yōu)勢演化可以看出(圖7),汊道#3相對梯度優(yōu)勢大于其他汊道,模擬后期汊道#3 相對梯度優(yōu)勢減弱,這顯示了汊道逐漸曲流化,長度增加。根據(jù)Burgeet al.[32]關(guān)于Renous河的研究,維持長期存在的分汊河道往往具有相近的流徑長度,這樣才能使兩個河道形成一種平衡,如果兩者的長度相差大,會逐漸打破這種平衡,引起低效輸砂河道的漸棄。在本次研究的河段中,上河曲段的汊道#3 相對梯度優(yōu)勢最大,而且隨著河流的演化,這種分流分沙能力越來越強(qiáng),因此汊道#3 廢棄的概率最小。不同汊道的河床泥沙粒徑中值分布存在顯著差異,汊道#3 床砂粒徑中值最大,水動力更強(qiáng),不易廢棄。從河床高程演化來看,原河道的河床一直保持低高程優(yōu)勢(圖8),但這并沒有改變其廢棄的趨勢,可見僅以低高程河床的優(yōu)勢不能維持汊道的活躍。為了進(jìn)一步分析江心洲不同部位的汊道特征,將原主流方向與切灘水流方向的夾角定義為分汊角。李志威等[33-34]通過統(tǒng)計亞馬遜流域、額爾齊斯河等河共67個截彎特征,發(fā)現(xiàn)分汊角變化區(qū)間是24°~103°,平均值54.8°,現(xiàn)代沉積研究河段的衛(wèi)星圖像統(tǒng)計分汊角范圍為20°~85°,研究區(qū)模擬結(jié)果測算的汊道#2 分汊角區(qū)間為25°~70°,上述結(jié)果與模擬結(jié)果匹配度均良好。

    圖7 各汊道相對梯度優(yōu)勢隨時間演化圖Fig.7 Evolution of relative gradient advantage parameter in different channels

    圖8 各汊道高程隨時間演化圖Fig.8 Evolution of elevation in different channels

    除了分汊角外,江心洲頭部與上游水流軸線的對稱關(guān)系也是影響汊道口形態(tài)對江心洲發(fā)育控制作用的因素。將江心洲一側(cè)的頭部切線偏離上游河道主流線的角度定義為分汊偏移角。一個江心洲頭部兩側(cè)各自形成的汊道中,分汊偏移角較大的一側(cè),下游汊道易淤堵,因此,分汊偏移角可以指示江心洲進(jìn)一步遷移或擴(kuò)張后連岸閉合的可能性。分汊偏移角統(tǒng)計顯示(圖9),上游分汊口BP1 江心洲(圖5)的一側(cè)分汊偏移角a11 總是大于另一側(cè)分汊偏移角a12,該分汊不對稱性較強(qiáng),在450 TS后,#1汊道的流速已經(jīng)無法識別主流線趨勢,而偏移角a12已經(jīng)下降到近于零,表明原河道開始逐漸廢棄,分汊的不對稱性與短期汊道演化之間關(guān)系是明確的。對于上游分汊口BP2(圖3),偏移角a22 長期小于a21,第二期洪水(500 TS)后,偏移角差異進(jìn)一步增大,最大應(yīng)力和最大流速集中分布區(qū)域從兩河彎的轉(zhuǎn)折端逐漸轉(zhuǎn)移至#3、#4汊道一帶。結(jié)合最大應(yīng)力分布、流速分布特征可以看出,#3 汊道在第二期洪水后分流分沙能力已經(jīng)最強(qiáng),可推測#2汊道后期廢棄的概率很高。

    圖9 各分汊偏移角隨時間演化圖Fig.9 Evolution of bifurcation offset angle in different channels

    根據(jù)以上參數(shù)與汊道的關(guān)系可得出,汊道位于河曲底端時,更易發(fā)生主次汊道轉(zhuǎn)化;而發(fā)育于河曲頂端附近,不穩(wěn)定因素更多,汊道的演化方向預(yù)測難度加大。另外,本文只涉及串溝型截彎形成的汊道演化分析,而頂沖截彎、頸口截彎、決口改道等方式作用機(jī)制及對于汊道的影響存在差異[35],在此不做探討。

    3.3 曲流串溝型江心洲演化過程

    對于曲流河而言,江心洲的形成造成分流分沙,使水動力復(fù)雜化,一方面,曲流兩岸淤積沖刷影響河道平面形態(tài);另一方面,由于江心洲的壅水作用和不同汊道的阻力差異形成橫向比降及環(huán)流作用,江心洲本身也會發(fā)生遷移、沉積侵蝕的變化,阻塞河道、沖刷大堤、破壞港口等,對于河流治理者和水利研究人員都是一個挑戰(zhàn)[36-38]。

    下河曲段的模擬結(jié)果切片顯示,一期洪水后,洪水切灘截彎,形成江心洲(圖10)。此時江心洲主體為點壩殘余沉積,側(cè)積層傾角4°~5°,傾角與河流側(cè)向遷移規(guī)模和河床地貌有關(guān)。由于江心洲兩側(cè)水動力軸線偏移形成不對稱水流,江心洲東側(cè)侵蝕作用為主,是點壩發(fā)育的部分沉積物來源;在江心洲西側(cè),洪水期快速加積與平水期緩慢側(cè)積交替,洪水期末,有細(xì)粒隔層發(fā)育。洲頭水動力強(qiáng)(圖11a),而上游輸送泥沙量少,平水期沉積較穩(wěn)定,洪水期沖刷快速蝕退,在沙洲形成較大的坡度,洲頭細(xì)粒沉積被侵蝕和攜帶至下游,無法保存;江心洲呈中間厚、兩翼薄的透鏡狀,洲主體向凹岸下游方向遷移。在這一過程中,江心洲西側(cè)汊道不斷萎縮,東側(cè)汊道不斷擴(kuò)張。江心洲頂部漫灘細(xì)粒緩慢加積(圖11c),沉積層面平緩,與下部斜層理呈大角度接觸,在洪水期頂部有流槽形成或消亡,使江心洲內(nèi)部結(jié)構(gòu)復(fù)雜化。受洲體遮擋,洲尾沖淤強(qiáng)度一般遠(yuǎn)小于洲頭,形成洲尾向下游加積(圖11b)。在江心洲有尾翼發(fā)育時,尾翼延伸規(guī)模與兩側(cè)汊道的分沙分流能力相對強(qiáng)弱密切相關(guān),主汊道一側(cè)坡度陡,次汊道一側(cè)坡度緩。統(tǒng)計發(fā)現(xiàn),模擬演化結(jié)束時該江心洲順流軸長約1 480 m,垂直流向上軸長約630 m,與實際松花江位于該處的江心洲規(guī)模分布范圍較為一致,表明模擬結(jié)果具有較高可信度。

    圖10 串溝型江心洲平面演化特征模擬結(jié)果Fig.10 Simulated planar characteristics of evolving ditch-scouring mid-channel bar

    圖11 江心洲不同位置地貌的演化特征(剖面線位置見圖10)Fig.11 Geomorphology evolution of mid-channel bar at different locations (profile location in Fig.10)

    衛(wèi)星圖像顯示,上述模擬中出現(xiàn)的串溝型江心洲向下游遷移和側(cè)向遷移的現(xiàn)象,在松花江流域廣泛存在。圖12 所示江心洲上仍保留著點壩側(cè)向遷移的弧形線與流槽殘留特征,是曲流河點壩沉積的重要標(biāo)志。從2005 年3 月到2010 年5 月,江心洲的洲頭位置向下游遷移了270 m,江心洲受非對稱水流侵蝕,南側(cè)為主水流經(jīng)過,侵蝕作用為主,平均侵蝕速率約14 m/年;北側(cè)為次水道,江心洲主體中部較為穩(wěn)定,近尾部侵蝕速率逐漸增大,最大侵蝕速率約44 m/年,洲尾淤長約420 m。從2010 年5 月到2011 年10 月,江心洲洲頭向下游遷移90 m,北側(cè)江心洲側(cè)向加積為主或者處于沉積—侵蝕的平衡狀態(tài),江心洲向上游區(qū)域側(cè)積速率最大,約90 m/年;南側(cè)汊道對近洲頭區(qū)域沖刷侵蝕,江心洲主體側(cè)向加積或維持平衡,最大側(cè)向加積速率約110 m/年;下游狹隘化河道使河道水動力發(fā)生變化,尾部沉積延伸緩慢近于停滯。上述模擬結(jié)果顯示的江心洲演化與不同時期衛(wèi)星照片揭示的現(xiàn)代沉積特征一致,從而印證了模擬結(jié)果。在這個過程中,大洪水對于該沉積體的演化意義重大,洪水作為集中釋放的高能量對于原地貌的變形(包括新汊道形成、點壩結(jié)構(gòu)破壞)起主要作用,雙向水流作用下汊道的演化愈活躍,對應(yīng)著側(cè)向加積、洲頭侵蝕、洲尾加積等沉積改造越復(fù)雜,曲流江心洲與一般點壩的差異性越大。

    圖12 松花江江心洲演化過程Fig.12 Evolution of mid-channel bar in Songhua River

    4 曲流串溝型江心洲演化模式

    根據(jù)上述模擬結(jié)果和不同時間的衛(wèi)星照片顯示的河流形態(tài)與參數(shù)特征和演化,高彎度曲流串溝截彎過程可以歸納為4個發(fā)育階段(圖13)。

    (1)曲流河發(fā)生側(cè)向遷移和順流遷移,伴隨著側(cè)向遷移,河流的彎曲度增大,每一期大規(guī)模的側(cè)向遷移都對應(yīng)了一起強(qiáng)水動力事件從產(chǎn)生到結(jié)束的完整過程。這些周期性的水動力變化事件在點壩上形成一期側(cè)積體,每一期側(cè)積體在地表形成一個弧形脊灘,不同期次間發(fā)育流槽(圖13a)。

    (2)在洪水期隨著水位上升,溢岸沖刷使得串溝平面上延伸,垂向下蝕,規(guī)模不斷加深拓寬最終形成新汊道,被汊道分開、孤立存在的沉積體——江心洲雛形初現(xiàn)(圖13b)。

    (3)隨著新汊道的形成,江心洲前端兩側(cè)分汊偏移角的差別使得河道形成不對稱分流,造成兩側(cè)汊道的水動力差異,由于不對稱的分流分沙等原因,在原河道入口淤積形成塞狀壩,汊道狹隘化,阻礙水流進(jìn)入,平水期水體近乎停滯(圖13c)。

    (4)洪水退去后,部分泥沙淤落于舊河道,同時由于不對稱分流分沙導(dǎo)致的舊河道水動力減弱,經(jīng)歷漫長演化,細(xì)粒沉積物垂向加積使得舊河道最終廢棄充填,新河道的遷移也開始進(jìn)行(圖13d)。

    圖13 曲流河串溝型截彎模式演化圖Fig.13 Evolution pattern of ditch-scouring cutoffs in meandering river

    串溝型江心洲演化伴隨著汊道的形成、轉(zhuǎn)化與廢棄,是串溝型截彎演化過程的一部分,大型江心洲一般可以較長期穩(wěn)定存在。根據(jù)模擬結(jié)果和現(xiàn)代沉積地貌,在江心洲發(fā)育的高彎度曲流河中,存在六種沉積單元,包括活躍河道、廢棄河道、天然堤、江心洲、越岸沉積、點壩。結(jié)合串溝截彎和江心洲演化過程及不同成因單元的空間配置,建立了曲流串溝型江心洲的發(fā)育模式(圖14)。

    圖14 曲流串溝型江心洲發(fā)育模式Fig.14 Depositional model of ditch-scouring mid-channel bar in meandering river

    雖然曲流串溝型江心洲與辮狀河心灘都是兩側(cè)水道發(fā)育的砂體沉積單元,而且在洪水期上游端都遭受侵蝕,在尾部加積,但是二者存在本質(zhì)區(qū)別。心灘的演化通常包括心灘生長及向下游遷移、側(cè)向遷移、復(fù)合心灘壩的形成及分割等過程,涉及到垂向加積、側(cè)向加積、填積等多種沉積模式[39-40],主體結(jié)構(gòu)以多期次垂向加積為主,在心灘主體內(nèi)部發(fā)育大量的落淤層,心灘壩尾沉積可以形成復(fù)合壩,橫向壩的雙尾翼結(jié)構(gòu)常見[41]。而串溝型江心洲主體結(jié)構(gòu)不是典型的垂向加積,它是以側(cè)積為主的點壩在改造殘留后的結(jié)果,在發(fā)育過程中有向下游加積、兩翼反向側(cè)積等多種沉積模式,其中側(cè)向加積的沉積模式起主導(dǎo)作用,這控制了曲流串溝型江心洲內(nèi)部砂體的非均質(zhì)性格架,是導(dǎo)致其內(nèi)部結(jié)構(gòu)復(fù)雜化的主要原因。長期改造可能使點壩頂部的脊槽特征逐漸弱化,江心洲頂部有洪水退卻時形成的細(xì)粒沉積,在平水期穩(wěn)定保存,尾部沉積受限于河岸,發(fā)育規(guī)模一般較小,下游延伸范圍有限。當(dāng)江心洲進(jìn)入演化末期,兩側(cè)汊道先后廢棄填充,在這一過程中,通常先在上游汊道發(fā)育塞狀壩,隨后緩慢加積廢棄化,最終與泛濫平原成為一體,或與凸岸點壩結(jié)合,進(jìn)入新的改造階段,這加劇了曲流河點壩砂體的沉積非均質(zhì)性。對于地下的曲流河油氣儲層來說,這一發(fā)育過程導(dǎo)致的儲層強(qiáng)非均質(zhì)性使地下油氣水的滲流復(fù)雜化,給油氣開發(fā)帶來諸多問題。

    5 結(jié)論

    (1)在曲流河中,平水期河曲段不斷凹蝕凸增,洪水期發(fā)生漫灘沖刷,串溝發(fā)育擴(kuò)張,截彎取直,形成分汊河道及江心洲地貌。

    (2)汊道的存廢與分流分沙能力的不對稱性相關(guān),江心洲兩側(cè)汊道的存廢具有“競爭性”,可以通過分汊偏移角及相對梯度優(yōu)勢等參數(shù)評估和預(yù)測汊道消亡的可能性。分汊偏移角小、相對梯度優(yōu)勢大的一側(cè)更有分流分沙優(yōu)勢,維系汊道的活躍。河床高程可以影響分流分沙比,低位河床有更大的運輸通道截面積,但不是必要的條件。

    (3)曲流串溝型江心洲演化是基于對點壩沉積的改造,歷經(jīng)向下游加積、兩翼非對稱側(cè)積與侵蝕、垂向加積等多種沉積作用與過程。歷經(jīng)長期演化,受不穩(wěn)定水文條件和較高彎曲度的地貌影響,江心洲演化結(jié)束后,或并入凹岸與泛濫平原成為一體,或并入凸岸進(jìn)一步加劇點壩內(nèi)部結(jié)構(gòu)的復(fù)雜化。

    猜你喜歡
    江心洲流河河道
    七月
    三峽庫區(qū)不同蓄水階段下江心洲的時空演變分析*
    長江干流江心洲過去35年面積變化及其特征
    嘉陵江中游水電梯級開發(fā)對小龍門江心洲的形態(tài)演變與影響機(jī)制研究
    曲流河復(fù)合點壩砂體構(gòu)型表征及流體運移機(jī)理
    從江心洲到江“新”洲——“紅領(lǐng)巾研學(xué)院”課題研究活動
    少先隊活動(2021年3期)2021-06-11 05:43:34
    陽春高流墟鄉(xiāng)韻醉人
    源流(2018年7期)2018-12-03 13:26:38
    河道里的垃圾
    山溪性河道植物選擇與應(yīng)用
    南京市鼓樓區(qū)黑臭河道的治理
    嫩草影院精品99| 午夜福利视频1000在线观看 | 亚洲va日本ⅴa欧美va伊人久久| 最新美女视频免费是黄的| 在线播放国产精品三级| 免费不卡黄色视频| 人人澡人人妻人| 中文字幕高清在线视频| 看免费av毛片| 国产私拍福利视频在线观看| 亚洲 欧美 日韩 在线 免费| 亚洲精品在线美女| 老熟妇仑乱视频hdxx| 国产成+人综合+亚洲专区| 欧美+亚洲+日韩+国产| av免费在线观看网站| 精品国产超薄肉色丝袜足j| 美女国产高潮福利片在线看| 夜夜看夜夜爽夜夜摸| 日韩欧美国产一区二区入口| 久久狼人影院| 大陆偷拍与自拍| 精品国产乱码久久久久久男人| 欧美激情 高清一区二区三区| 久久香蕉国产精品| 大陆偷拍与自拍| 欧美人与性动交α欧美精品济南到| 欧美人与性动交α欧美精品济南到| 无人区码免费观看不卡| 精品第一国产精品| 久久久久国产精品人妻aⅴ院| 神马国产精品三级电影在线观看 | 天天添夜夜摸| 久久性视频一级片| 亚洲精品国产一区二区精华液| 成人18禁高潮啪啪吃奶动态图| 欧美日本视频| 男人舔女人的私密视频| 精品电影一区二区在线| 在线观看免费视频网站a站| 极品教师在线免费播放| 久久九九热精品免费| 亚洲中文字幕日韩| 熟妇人妻久久中文字幕3abv| 黄网站色视频无遮挡免费观看| 在线观看免费视频网站a站| 国产伦人伦偷精品视频| 欧美色欧美亚洲另类二区 | 身体一侧抽搐| 亚洲精品国产区一区二| 亚洲片人在线观看| 88av欧美| 18禁裸乳无遮挡免费网站照片 | 美女高潮喷水抽搐中文字幕| 国产亚洲欧美在线一区二区| 黄色视频,在线免费观看| 午夜免费激情av| 国产精品永久免费网站| 激情在线观看视频在线高清| 人人妻,人人澡人人爽秒播| 日本精品一区二区三区蜜桃| 国产一区二区三区综合在线观看| 亚洲成av片中文字幕在线观看| 18禁美女被吸乳视频| 欧美色欧美亚洲另类二区 | 日韩欧美一区二区三区在线观看| 一个人免费在线观看的高清视频| 最近最新中文字幕大全免费视频| 欧美日韩黄片免| 久久精品国产综合久久久| a在线观看视频网站| 亚洲国产欧美日韩在线播放| 亚洲自拍偷在线| 久久中文看片网| 在线观看免费视频网站a站| 色老头精品视频在线观看| 色老头精品视频在线观看| 免费在线观看影片大全网站| 日韩 欧美 亚洲 中文字幕| 国产三级在线视频| 亚洲国产精品成人综合色| 纯流量卡能插随身wifi吗| 久久久久久免费高清国产稀缺| 最近最新中文字幕大全免费视频| 午夜福利在线观看吧| 国产麻豆成人av免费视频| 日韩欧美在线二视频| 久久人人爽av亚洲精品天堂| 中文字幕高清在线视频| 久热爱精品视频在线9| 欧美国产精品va在线观看不卡| 在线观看www视频免费| 中文字幕人妻丝袜一区二区| 亚洲av成人一区二区三| 美女午夜性视频免费| 成在线人永久免费视频| 国产激情久久老熟女| 757午夜福利合集在线观看| avwww免费| 国产成人精品久久二区二区91| 国产成人欧美在线观看| 成人特级黄色片久久久久久久| 美女高潮喷水抽搐中文字幕| 亚洲成a人片在线一区二区| 99国产精品一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 欧美在线一区亚洲| 在线播放国产精品三级| 一进一出抽搐gif免费好疼| 亚洲精品美女久久av网站| 69精品国产乱码久久久| 精品高清国产在线一区| 亚洲少妇的诱惑av| 久久亚洲真实| 人妻丰满熟妇av一区二区三区| 一进一出抽搐gif免费好疼| 亚洲狠狠婷婷综合久久图片| 可以免费在线观看a视频的电影网站| 在线十欧美十亚洲十日本专区| 日韩 欧美 亚洲 中文字幕| 久久精品成人免费网站| 久久久久久久久免费视频了| 国产三级在线视频| 嫁个100分男人电影在线观看| 亚洲欧美日韩高清在线视频| 成人三级黄色视频| 国产亚洲av嫩草精品影院| 中文字幕久久专区| 波多野结衣av一区二区av| 18禁美女被吸乳视频| 亚洲精品国产色婷婷电影| 亚洲一区二区三区色噜噜| 亚洲av成人av| 国产精品亚洲av一区麻豆| 一边摸一边抽搐一进一出视频| 亚洲激情在线av| 精品欧美国产一区二区三| 看片在线看免费视频| 男女午夜视频在线观看| 亚洲精品久久国产高清桃花| 国产亚洲精品av在线| 久久国产精品男人的天堂亚洲| 色精品久久人妻99蜜桃| 麻豆一二三区av精品| 日本五十路高清| 黄频高清免费视频| 老司机福利观看| 国产亚洲av高清不卡| 国产精品二区激情视频| 成年女人毛片免费观看观看9| 老熟妇乱子伦视频在线观看| 久久香蕉国产精品| 免费久久久久久久精品成人欧美视频| 午夜视频精品福利| 亚洲一码二码三码区别大吗| 欧美久久黑人一区二区| 伊人久久大香线蕉亚洲五| 欧美日韩中文字幕国产精品一区二区三区 | 国产亚洲精品久久久久5区| 国产精品 欧美亚洲| 黄色 视频免费看| 亚洲中文av在线| 成在线人永久免费视频| av网站免费在线观看视频| 免费在线观看视频国产中文字幕亚洲| 老司机深夜福利视频在线观看| 中文亚洲av片在线观看爽| 男女做爰动态图高潮gif福利片 | 一个人观看的视频www高清免费观看 | 极品教师在线免费播放| 久久久精品国产亚洲av高清涩受| 午夜久久久久精精品| 激情在线观看视频在线高清| 亚洲国产精品sss在线观看| 女生性感内裤真人,穿戴方法视频| 母亲3免费完整高清在线观看| 99在线视频只有这里精品首页| 亚洲人成伊人成综合网2020| 级片在线观看| 999精品在线视频| 极品教师在线免费播放| 一级毛片女人18水好多| 免费少妇av软件| 中国美女看黄片| 国产精品野战在线观看| 欧美午夜高清在线| 可以免费在线观看a视频的电影网站| 亚洲欧美激情在线| 亚洲七黄色美女视频| 精品一区二区三区av网在线观看| 亚洲国产精品999在线| 亚洲精华国产精华精| 国产高清videossex| 女人被躁到高潮嗷嗷叫费观| 国产伦人伦偷精品视频| 日韩成人在线观看一区二区三区| 麻豆国产av国片精品| 中亚洲国语对白在线视频| 午夜成年电影在线免费观看| 一进一出抽搐gif免费好疼| 国产精品av久久久久免费| 亚洲男人的天堂狠狠| 亚洲 欧美 日韩 在线 免费| 亚洲黑人精品在线| 91在线观看av| 性少妇av在线| 久久午夜综合久久蜜桃| 91字幕亚洲| 在线天堂中文资源库| 亚洲人成电影免费在线| 国产亚洲av高清不卡| 国产成人精品在线电影| 成人亚洲精品av一区二区| 国产乱人伦免费视频| 十八禁网站免费在线| 中文字幕色久视频| 人人妻,人人澡人人爽秒播| 18美女黄网站色大片免费观看| 成人精品一区二区免费| av视频免费观看在线观看| 久久香蕉激情| 黄色成人免费大全| 亚洲狠狠婷婷综合久久图片| 级片在线观看| 午夜a级毛片| 国产欧美日韩综合在线一区二区| 亚洲专区国产一区二区| 91成年电影在线观看| √禁漫天堂资源中文www| 午夜福利高清视频| 午夜免费激情av| 国产成人免费无遮挡视频| 国产精品香港三级国产av潘金莲| 精品一区二区三区四区五区乱码| 欧美成人午夜精品| а√天堂www在线а√下载| 黄色a级毛片大全视频| 久久久久久免费高清国产稀缺| 日韩视频一区二区在线观看| 中文字幕久久专区| 欧美午夜高清在线| 欧美另类亚洲清纯唯美| 波多野结衣高清无吗| 久久精品成人免费网站| 曰老女人黄片| 国产亚洲精品久久久久5区| 久久青草综合色| 午夜a级毛片| 久久这里只有精品19| 成人特级黄色片久久久久久久| 亚洲欧美一区二区三区黑人| 欧美日本视频| 午夜福利免费观看在线| 99精品在免费线老司机午夜| 国产私拍福利视频在线观看| 色尼玛亚洲综合影院| 91成人精品电影| 成年人黄色毛片网站| 亚洲欧美日韩高清在线视频| 无限看片的www在线观看| 成年人黄色毛片网站| 91老司机精品| 国产精品综合久久久久久久免费 | 一进一出抽搐动态| 午夜激情av网站| 亚洲性夜色夜夜综合| 两个人看的免费小视频| 久久精品国产亚洲av香蕉五月| aaaaa片日本免费| 一区二区三区激情视频| 母亲3免费完整高清在线观看| 91九色精品人成在线观看| 国产私拍福利视频在线观看| 中文字幕人妻熟女乱码| 色综合婷婷激情| 高清黄色对白视频在线免费看| 后天国语完整版免费观看| 亚洲精品一区av在线观看| www.自偷自拍.com| 一夜夜www| 国产精品久久电影中文字幕| 亚洲国产高清在线一区二区三 | 亚洲精品中文字幕在线视频| 叶爱在线成人免费视频播放| 日韩欧美在线二视频| 亚洲电影在线观看av| 国产亚洲欧美在线一区二区| 亚洲成国产人片在线观看| 亚洲专区国产一区二区| 99在线人妻在线中文字幕| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产毛片av蜜桃av| 日本vs欧美在线观看视频| 亚洲精品中文字幕在线视频| 两性夫妻黄色片| 老司机在亚洲福利影院| 亚洲黑人精品在线| 老熟妇乱子伦视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 免费一级毛片在线播放高清视频 | 亚洲一区二区三区不卡视频| 日韩视频一区二区在线观看| 亚洲午夜理论影院| 久久久精品欧美日韩精品| 淫妇啪啪啪对白视频| 婷婷六月久久综合丁香| 国产精品久久电影中文字幕| 成人三级黄色视频| 十分钟在线观看高清视频www| 午夜福利高清视频| 精品午夜福利视频在线观看一区| 亚洲欧美日韩无卡精品| 中文字幕色久视频| 亚洲少妇的诱惑av| 老司机午夜十八禁免费视频| 精品人妻1区二区| 国产极品粉嫩免费观看在线| 亚洲熟女毛片儿| 久久久久久久久久久久大奶| av免费在线观看网站| 欧美成人性av电影在线观看| 欧美精品啪啪一区二区三区| 51午夜福利影视在线观看| 亚洲国产日韩欧美精品在线观看 | 性少妇av在线| 涩涩av久久男人的天堂| 国产97色在线日韩免费| 国产野战对白在线观看| x7x7x7水蜜桃| av中文乱码字幕在线| 亚洲精品久久国产高清桃花| 国产熟女xx| 天堂√8在线中文| 日韩欧美在线二视频| 久久婷婷成人综合色麻豆| 啪啪无遮挡十八禁网站| 国产精品99久久99久久久不卡| 最近最新中文字幕大全免费视频| 成人国产综合亚洲| 老鸭窝网址在线观看| 天天一区二区日本电影三级 | 女同久久另类99精品国产91| 嫁个100分男人电影在线观看| 麻豆成人av在线观看| 国产精品久久久久久人妻精品电影| 香蕉久久夜色| 黄频高清免费视频| 日韩精品中文字幕看吧| 免费人成视频x8x8入口观看| 午夜激情av网站| 欧美在线一区亚洲| 一级a爱视频在线免费观看| cao死你这个sao货| 亚洲欧美精品综合一区二区三区| 日韩大码丰满熟妇| 精品国产亚洲在线| 又大又爽又粗| 亚洲中文字幕日韩| 欧美色视频一区免费| 国产一区二区三区视频了| 香蕉国产在线看| 一区二区三区激情视频| 制服诱惑二区| 妹子高潮喷水视频| 久99久视频精品免费| 午夜视频精品福利| x7x7x7水蜜桃| 久久国产乱子伦精品免费另类| 女人精品久久久久毛片| 日本欧美视频一区| 午夜福利成人在线免费观看| 久热这里只有精品99| 丰满人妻熟妇乱又伦精品不卡| 色综合亚洲欧美另类图片| 级片在线观看| 亚洲av成人av| 久久久久久久久中文| 国产黄a三级三级三级人| 99精品在免费线老司机午夜| 国产成人精品在线电影| 91九色精品人成在线观看| 俄罗斯特黄特色一大片| 不卡av一区二区三区| 亚洲精品中文字幕一二三四区| 757午夜福利合集在线观看| 黄色视频,在线免费观看| 色播在线永久视频| 18禁美女被吸乳视频| 欧美激情久久久久久爽电影 | 亚洲中文av在线| 九色国产91popny在线| 免费在线观看完整版高清| 欧美久久黑人一区二区| АⅤ资源中文在线天堂| 亚洲av片天天在线观看| 午夜亚洲福利在线播放| 国语自产精品视频在线第100页| 亚洲九九香蕉| 欧美 亚洲 国产 日韩一| 香蕉丝袜av| 免费在线观看视频国产中文字幕亚洲| 老鸭窝网址在线观看| 亚洲五月色婷婷综合| 午夜福利视频1000在线观看 | 欧美中文日本在线观看视频| avwww免费| 成人国语在线视频| 美女扒开内裤让男人捅视频| 麻豆一二三区av精品| 丰满人妻熟妇乱又伦精品不卡| 老汉色∧v一级毛片| 岛国视频午夜一区免费看| 国产精品综合久久久久久久免费 | 男女床上黄色一级片免费看| 午夜精品久久久久久毛片777| 91九色精品人成在线观看| 在线免费观看的www视频| 老司机午夜福利在线观看视频| 老熟妇仑乱视频hdxx| 国产精品爽爽va在线观看网站 | 91九色精品人成在线观看| 国产蜜桃级精品一区二区三区| aaaaa片日本免费| 亚洲精品在线美女| 国产精品久久久久久人妻精品电影| 欧美乱码精品一区二区三区| a在线观看视频网站| www日本在线高清视频| 嫩草影院精品99| 国产激情欧美一区二区| 国产欧美日韩一区二区三| 久久久久国产一级毛片高清牌| 神马国产精品三级电影在线观看 | 一二三四在线观看免费中文在| 欧美老熟妇乱子伦牲交| 精品午夜福利视频在线观看一区| 精品福利观看| 麻豆成人av在线观看| 免费女性裸体啪啪无遮挡网站| 免费久久久久久久精品成人欧美视频| 一区二区三区高清视频在线| 国产伦人伦偷精品视频| 国产人伦9x9x在线观看| 亚洲精品中文字幕在线视频| 午夜亚洲福利在线播放| 免费观看精品视频网站| 在线视频色国产色| 99香蕉大伊视频| 欧美日韩精品网址| 91av网站免费观看| 一区二区三区国产精品乱码| 可以在线观看毛片的网站| 此物有八面人人有两片| 亚洲五月婷婷丁香| 一本综合久久免费| 欧美激情极品国产一区二区三区| 国产片内射在线| 午夜a级毛片| 久久草成人影院| 亚洲美女黄片视频| 精品一区二区三区四区五区乱码| 免费在线观看影片大全网站| a在线观看视频网站| 老司机在亚洲福利影院| 国产高清有码在线观看视频 | 国产精品久久久久久精品电影 | 91老司机精品| 精品乱码久久久久久99久播| 欧美黄色淫秽网站| 高潮久久久久久久久久久不卡| 亚洲中文日韩欧美视频| 久久精品影院6| 亚洲精品国产色婷婷电影| 亚洲片人在线观看| 一个人观看的视频www高清免费观看 | 欧美在线黄色| 欧美日韩中文字幕国产精品一区二区三区 | 别揉我奶头~嗯~啊~动态视频| 在线十欧美十亚洲十日本专区| 亚洲精品中文字幕一二三四区| 国产伦一二天堂av在线观看| 91av网站免费观看| 国产激情欧美一区二区| 色综合站精品国产| 精品高清国产在线一区| 国产一区二区激情短视频| 啦啦啦韩国在线观看视频| 亚洲av美国av| 最近最新中文字幕大全电影3 | 亚洲美女黄片视频| 精品无人区乱码1区二区| 12—13女人毛片做爰片一| 久久午夜亚洲精品久久| 亚洲一区二区三区色噜噜| 99riav亚洲国产免费| 变态另类成人亚洲欧美熟女 | 手机成人av网站| 亚洲性夜色夜夜综合| 国产亚洲精品综合一区在线观看 | 一级黄色大片毛片| 电影成人av| 成人精品一区二区免费| 手机成人av网站| 淫秽高清视频在线观看| 丁香六月欧美| 国产私拍福利视频在线观看| 91麻豆av在线| 国产亚洲精品第一综合不卡| 亚洲成a人片在线一区二区| 97碰自拍视频| 妹子高潮喷水视频| 夜夜爽天天搞| netflix在线观看网站| 变态另类成人亚洲欧美熟女 | 99国产精品99久久久久| 亚洲av日韩精品久久久久久密| 热re99久久国产66热| 成在线人永久免费视频| av电影中文网址| 一夜夜www| 成人国产综合亚洲| av超薄肉色丝袜交足视频| 久热爱精品视频在线9| 久久久久九九精品影院| 亚洲电影在线观看av| 成人手机av| 国产高清视频在线播放一区| 18禁黄网站禁片午夜丰满| www.自偷自拍.com| 很黄的视频免费| 91成人精品电影| 亚洲av电影不卡..在线观看| 国产高清激情床上av| 国产av一区在线观看免费| 99热只有精品国产| 中文字幕久久专区| 成人亚洲精品一区在线观看| 在线观看舔阴道视频| 黄频高清免费视频| 色老头精品视频在线观看| 久久亚洲真实| 又紧又爽又黄一区二区| 99riav亚洲国产免费| 色哟哟哟哟哟哟| 精品熟女少妇八av免费久了| 伊人久久大香线蕉亚洲五| 非洲黑人性xxxx精品又粗又长| 午夜影院日韩av| 亚洲精品国产一区二区精华液| 婷婷六月久久综合丁香| 国产精品一区二区免费欧美| 亚洲精品国产精品久久久不卡| 波多野结衣巨乳人妻| 欧美日韩乱码在线| 色老头精品视频在线观看| 一边摸一边抽搐一进一小说| 欧美成人午夜精品| av网站免费在线观看视频| 国产人伦9x9x在线观看| 黄色视频不卡| 香蕉久久夜色| av视频在线观看入口| 黄色成人免费大全| 亚洲最大成人中文| 色精品久久人妻99蜜桃| 精品不卡国产一区二区三区| 久9热在线精品视频| 人人妻,人人澡人人爽秒播| 久久草成人影院| 国产精品久久久久久精品电影 | 亚洲精品av麻豆狂野| 成在线人永久免费视频| 又黄又爽又免费观看的视频| 一级,二级,三级黄色视频| 日本黄色视频三级网站网址| 国产欧美日韩一区二区精品| 久久久久亚洲av毛片大全| 欧美一级毛片孕妇| 丰满的人妻完整版| 国产精品电影一区二区三区| 香蕉丝袜av| 变态另类成人亚洲欧美熟女 | 午夜视频精品福利| 十分钟在线观看高清视频www| av欧美777| 中文字幕久久专区| 丝袜美腿诱惑在线| 国产成+人综合+亚洲专区| 一个人免费在线观看的高清视频| 丝袜美足系列| 高清黄色对白视频在线免费看| 亚洲人成网站在线播放欧美日韩| 真人一进一出gif抽搐免费| 亚洲少妇的诱惑av| 国产成人欧美在线观看| 日本五十路高清| 99久久精品国产亚洲精品| 久久狼人影院| www.熟女人妻精品国产| 一区福利在线观看| 美女扒开内裤让男人捅视频| 欧美黑人欧美精品刺激| 国产熟女xx| 视频在线观看一区二区三区| 丝袜在线中文字幕| 在线av久久热| 宅男免费午夜| 精品高清国产在线一区| 日本a在线网址| 精品国产超薄肉色丝袜足j| 免费在线观看亚洲国产| 日韩精品免费视频一区二区三区| 黄片播放在线免费| 极品教师在线免费播放| 黄色成人免费大全|