張可，吳勝和，馮文杰，鄭定業(yè)，喻宸，劉照瑋

1.中國(guó)石油大學(xué)(北京)地球科學(xué)學(xué)院 油氣資源與探測(cè)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 102249 2.中國(guó)石化石油勘探開(kāi)發(fā)研究院，北京 100083

0 引言

辮狀河是一種十分常見(jiàn)的沉積體系，為河谷較為平直、低彎度、坡降大、洪泛間歇性大、流量變化大、碎屑物粗、以推移質(zhì)為主、多河道分叉合并、不斷遷移改道的河流[1- 7]。根據(jù)沉積物粒度的不同，辮狀河可分為砂質(zhì)辮狀河與礫質(zhì)辮狀河。我國(guó)陸相盆地中廣泛發(fā)育砂質(zhì)辮狀河儲(chǔ)層，如鄂爾多斯盆地[8- 9]、渤海灣盆地[10]以及松遼盆地[11]等。心灘壩作為辮狀河中重要的沉積單元，一直是國(guó)內(nèi)外學(xué)者們研究的重點(diǎn)。

部分學(xué)者通過(guò)露頭、巖芯、測(cè)井以及動(dòng)態(tài)資料進(jìn)行分析、從單井上識(shí)別心灘壩、辮狀水道以及心灘壩內(nèi)部落淤層等[12- 15]，是基于井資料的靜態(tài)分析過(guò)程，并不能建立心灘壩的動(dòng)態(tài)演化過(guò)程及空間分布特征。陳玉琨等人從沉積過(guò)程的角度明確了心灘壩的三種類型及不同類型心灘壩內(nèi)部落淤層的展布樣式，認(rèn)為心灘壩類型受兩側(cè)辮狀水道的影響，但其研究?jī)H限于地質(zhì)分析[16]。近年來(lái)，也有部分學(xué)者逐漸開(kāi)始關(guān)注砂質(zhì)辮狀河以及其內(nèi)部心灘壩沉積演化過(guò)程[17- 18]。Ashworthetal.[19]觀察研究大型砂質(zhì)辮狀河賈木納河心灘壩的沉積演化過(guò)程，認(rèn)為心灘壩發(fā)育過(guò)程分為四個(gè)階段：1)心灘壩的形成及生長(zhǎng)；2)側(cè)向加積作用造成壩的加寬；3)壩的收縮及壩尾的延長(zhǎng)；4)下游新支流的形成，四個(gè)階段均與辮狀河水流強(qiáng)度相關(guān)且具有時(shí)間先后次序。Schuurmanetal.[20- 21]從沉積數(shù)值模擬的角度研究砂質(zhì)辮狀河中辮狀水道的開(kāi)啟與關(guān)閉機(jī)制以及心灘壩演化過(guò)程，但并未明確心灘壩兩側(cè)水道對(duì)心灘壩的影響，未建立心灘壩發(fā)育演化的完整動(dòng)態(tài)過(guò)程。

1 研究方法及研究參數(shù)設(shè)置

砂質(zhì)辮狀河以底負(fù)載搬運(yùn)為主，其他特征基本相似，因此本文以前人模擬數(shù)據(jù)為基礎(chǔ)[20]，設(shè)置沉積數(shù)值模擬的邊界條件等，采用Delft3D軟件模擬心灘壩發(fā)育過(guò)程，分析其沉積特征，并利用現(xiàn)代沉積相互驗(yàn)證，最終建立心灘壩演化模式。

1.1 沉積數(shù)值模擬原理

Delft3D軟件是由荷蘭杜蘭大學(xué)開(kāi)發(fā)的一套三維水動(dòng)力數(shù)值模擬軟件，其以描述水動(dòng)力的納維—斯托克斯方程(N- S方程)為基礎(chǔ)，并結(jié)合物質(zhì)平衡方程，實(shí)現(xiàn)沉積物搬運(yùn)以及沉積地貌演變的三維模擬。其沉積模擬過(guò)程通過(guò)以下方程實(shí)現(xiàn)：

式中，x為X方向坐標(biāo)值(m)；y為Y方向坐標(biāo)值(m)；zw為自由面水位(m)；u為X方向上的流速(m/s)；v為Y方向上的流速(m/s)；h為水深(m)；C為Chezy系數(shù)(m1/2/s)；g為重力加速度(m/s2)；V為水平渦動(dòng)黏度(m2/s)；Fx、Fy分別為X和Y方向上的輻射應(yīng)力梯度(m/s2)。

Delft3D軟件中由設(shè)定的邊界條件及底床含泥砂量決定泥砂輸運(yùn)的總量。不同性質(zhì)的泥砂以不同的方式進(jìn)行搬運(yùn)，由不同的輸砂方程描述。

1.2 沉積數(shù)值模擬參數(shù)及邊界條件定義

通過(guò)統(tǒng)計(jì)露頭以及雅魯藏布江等多個(gè)常年流水的現(xiàn)代砂質(zhì)辮狀河資料，辮狀河寬度為300～11 000 m，心灘壩長(zhǎng)度為500～12 000 m，心灘壩寬度為100～4 000 m，辮狀水道寬度為40～800 m，為了能完整再現(xiàn)心灘壩的沉積演化過(guò)程，且同時(shí)保持模擬耗時(shí)在可接受的范圍內(nèi)，本次模擬定義網(wǎng)格大小為50 m×30 m，垂向網(wǎng)格數(shù)為7個(gè)，保證心灘壩及辮狀水道自形成起至少占據(jù)1個(gè)網(wǎng)格。模擬區(qū)大小為80 km×3 km，滿足研究需求。上下兩側(cè)邊界為封閉邊界，左右兩側(cè)為開(kāi)放邊界，左側(cè)邊界為進(jìn)水口，右側(cè)邊界為出水口，坡度為9.30×10-5(圖1)。左側(cè)邊界均勻設(shè)置20個(gè)供水口，每個(gè)供水口之間流量存在差異，保證辮狀河斷面不同位置流速隨機(jī)變化，但總水流量恒定。在模擬網(wǎng)格及邊界定義的基礎(chǔ)上，設(shè)定泥(黏性)與砂(非黏性)兩類沉積物，其中砂質(zhì)沉積物粒度中值為0.2 mm，初始沉積厚度為100 m，初始流量為2 000 m3/s。由于不同的流量、坡度等參數(shù)雖然影響形成心灘壩以及辮狀水道的規(guī)模，但不會(huì)完全改變心灘壩發(fā)育演化過(guò)程[30]，因此參考前人模擬參數(shù)[20]，設(shè)定本次模擬的其他泥沙動(dòng)力學(xué)參數(shù)(表1)。本次模擬時(shí)間為284天，反映2～3年的高流量沉積過(guò)程，能夠再現(xiàn)心灘壩完整的發(fā)育演化過(guò)程。

圖1 沉積數(shù)值模擬邊界及底形設(shè)置Fig.1 Boundary and sedimentary bed form of sedimentary numerical simulation

參數(shù)項(xiàng)模擬設(shè)定值底床坡度9.30×10-5砂質(zhì)沉積物干容重/(kg/m3)1650泥質(zhì)沉積物干容重/(kg/m3)500砂泥比3∶2沉積厚度變化記錄下限/m0.1重力加速度/(m/s2)9.81水體密度/(kg/m3)1000水平渦流黏度/(m2/s)1垂直渦流黏度/(m2/s)0.0001模擬時(shí)間步長(zhǎng)/min0.2模擬時(shí)長(zhǎng)/day284模擬時(shí)間加速倍數(shù)10網(wǎng)格單元大小/m50×30總網(wǎng)格數(shù)/個(gè)1602×162×7

2 心灘壩生長(zhǎng)演化過(guò)程

心灘壩自形成后，不斷受到辮狀水道的沖刷改造，影響生長(zhǎng)及運(yùn)移，因此心灘壩的形成演化及最終形態(tài)樣式主要受心灘壩上游、下游以及兩側(cè)辮狀水道控制，不同的水流條件下，形成的心灘壩也具有不同的內(nèi)部構(gòu)型特征。在辮狀水道的控制下，發(fā)生在心灘壩上的沉積作用主要包括垂向加積、順流加積、側(cè)向加積、漫積以及填積等[31]，而對(duì)心灘壩發(fā)育起控制作用的主要是前三種沉積作用。辮狀水道的流態(tài)差異會(huì)造成不同的沉積過(guò)程：1)當(dāng)兩側(cè)水流為高能對(duì)稱水流時(shí)，形成縱向沙壩，沉積作用主要為順流加積與垂向加積作用；2)當(dāng)兩側(cè)水流為低能對(duì)稱水流時(shí)，形成橫向沙壩，沉積作用主要為順流加積與垂向加積作用；3)當(dāng)兩側(cè)水流為不對(duì)稱水流時(shí)，形成斜列沙壩，沉積作用主要為側(cè)向加積作用。

本次實(shí)驗(yàn)?zāi)M了心灘壩形成及演化的完整過(guò)程(圖2)。砂質(zhì)辮狀河沉積物以底負(fù)載搬運(yùn)為主，沉積物隨水流逐漸向下游方向搬運(yùn)，沉積物的搬運(yùn)過(guò)程也是心灘壩的形成及演化過(guò)程。心灘壩在局部高部位堆積形成，接受辮狀水道沖刷、改造，最終復(fù)合，保存下來(lái)的心灘壩都是經(jīng)過(guò)復(fù)合而成的大型心灘壩，而一部分心灘壩隨著演化過(guò)程而消亡。其發(fā)育演化主要包括三個(gè)階段：1)心灘壩的形成、生長(zhǎng)及向下游方向遷移，2)心灘壩的側(cè)向遷移，3)“壩尾沉積”與復(fù)合心灘壩的形成。三個(gè)階段發(fā)育無(wú)時(shí)間先后次序，心灘壩形成初期三個(gè)階段同時(shí)存在，各階段延續(xù)時(shí)間與水深、兩側(cè)辮狀水道的流速以及下游新的心灘壩發(fā)育位置有關(guān)。

2.1 心灘壩的形成、生長(zhǎng)及由于兩側(cè)辮狀水道的沖刷作用向下游方向遷移

心灘壩雛形是受初始底床控制的，由于模擬實(shí)驗(yàn)所設(shè)置的沉積底形高低不平，辮狀水道搬運(yùn)的沉積物此時(shí)在局部高地卸載，不斷堆積，形成心灘壩的雛形，其規(guī)模較小。在水流的持續(xù)作用下，上游沉積物被逐漸搬運(yùn)過(guò)來(lái)，在心灘壩雛形兩側(cè)或頂部堆積，可形成單一心灘壩。當(dāng)壩體達(dá)到一定規(guī)模之后，辮狀水道在心灘壩位置處出現(xiàn)明顯的分叉，心灘壩迎水面兩側(cè)均遭受辮狀水道的沖刷作用，逐漸被侵蝕，而其背水面則不斷接受上游被侵蝕的沉積物而不斷增生，這種現(xiàn)象在宏觀上表現(xiàn)為心灘壩向下游方向遷移。

2014—2017年，采用資料查閱及樣線調(diào)查相結(jié)合的方法，根據(jù)不同海拔高度、坡向、植被類型設(shè)計(jì)調(diào)查路線，設(shè)計(jì)調(diào)查線路60余條，圍繞鴻圖嶂山地共選擇6個(gè)駐點(diǎn)，包括小溪村、龍?zhí)洞?、大峽谷、馬山村、貴人村，對(duì)鴻圖嶂山地進(jìn)行多次的植物調(diào)查、照片拍攝、標(biāo)本采集及生境條件記錄。從最低海拔361 m到最高海拔鴻圖嶂頂峰1277.4 m，對(duì)包括常綠闊葉林、針闊混交林、高山矮林與山頂灌草叢等不同類型的植物群落進(jìn)行調(diào)查記錄。對(duì)采集的標(biāo)本進(jìn)行鑒定，并查閱《中國(guó)景觀植物》[10]、《廣東植物志》[11]等工具書(shū)，結(jié)合野生植物的觀賞特點(diǎn)、生活型及適應(yīng)性作為觀賞植物選取的標(biāo)準(zhǔn)。

沉積模擬實(shí)驗(yàn)135步長(zhǎng)時(shí)，以初始沉積厚度為基準(zhǔn)，心灘壩A最大沉積厚度為5.8 m，沉積寬度為600 m；當(dāng)模擬步長(zhǎng)為145時(shí)，同一位置，最大沉積厚度為4 m，沉積寬度為420 m，10個(gè)時(shí)間步長(zhǎng)內(nèi)，心灘壩A遭受辮狀水道1側(cè)向侵蝕100 m，辮狀水道2側(cè)向侵蝕80 m(圖3a，b)，向下游遷移0.31 km。實(shí)際上，在砂質(zhì)辮狀河中廣泛存在單一心灘壩向下游遷移的現(xiàn)象，在加拿大南薩斯喀徹溫河中十分常見(jiàn)，衛(wèi)星圖像顯示，兩個(gè)心灘壩均向下游方向遷移，一年零兩個(gè)月的時(shí)間內(nèi)遷移距離分別為146.4 m、111.5 m(圖4)。

圖2 沉積模擬心灘壩發(fā)育演化過(guò)程Fig.2 Evolution of bar in sedimentary numerical simulation

圖3 心灘壩A頭部遭受兩側(cè)辮狀水道侵蝕(圖a剖面位置為圖b黃線位置)Fig.3 The erosion of bar due to scouring

為確定心灘壩向下游遷移位置與速率的變化情況，避免封閉邊界的影響，選擇位于辮狀河中部的ABCD四個(gè)心灘壩(圖5)，分析表明，心灘壩最初遷移速度可達(dá)0.043 km/step，隨著形成時(shí)間增長(zhǎng)，遷移速率雖有小幅度的波動(dòng)，但整體向下遷移速率逐漸減小，最終遷移速率為0，心灘壩穩(wěn)定在一位置處。心灘壩形成并向下遷移至穩(wěn)定共需要25～50個(gè)步長(zhǎng)，對(duì)應(yīng)6—12個(gè)月(圖6)。

心灘壩逐漸趨于穩(wěn)定的原因在于心灘壩兩側(cè)辮狀水道流速的降低。在向下遷移的過(guò)程中，心灘壩不斷生長(zhǎng)，隨著壩頂不斷接近水平面，壩頂處水流速度及輸砂量逐漸降低，而頂面較低洼部位水流開(kāi)始匯聚，不斷沖刷，形成小水道。水道的形成造成心灘壩頭部辮狀水道水流分散，流速減小，從而導(dǎo)致心灘壩向下游遷移速率減小，最終達(dá)到穩(wěn)定狀態(tài)。以心灘壩D為例，圖7a為沉積厚度平面分布圖，7b為相同演化時(shí)間相同位置所對(duì)應(yīng)的水深，從130～150個(gè)時(shí)間步長(zhǎng)，心灘壩不斷增高，壩頂水流速度逐漸降低。在150個(gè)時(shí)間步長(zhǎng)時(shí)，壩頂接近水平面，壩頂水流開(kāi)始匯聚在較低部位，形成小水道，此時(shí)心灘壩頭部的辮狀水道不僅會(huì)分叉到心灘壩兩側(cè)，水流還會(huì)供給壩頂?shù)娜齻€(gè)小水道(圖7a中紅色的箭頭代表水流方向)，到160個(gè)時(shí)間步長(zhǎng)時(shí)，壩頂水道加深加寬，流速變大。在整個(gè)過(guò)程中，由于壩頂新形成的水道的分流作用導(dǎo)致心灘壩頭部?jī)蓚?cè)辮狀水道流速降低，心灘壩位置逐漸穩(wěn)定。值得注意的是，壩頂水道形成之后，不斷加深，最終會(huì)切割心灘壩體，壩頂水道成為兩個(gè)心灘壩之間的辮狀水道。

2.2 心灘壩中下部由于兩側(cè)水道的不對(duì)稱性造成心灘壩的側(cè)向遷移

心灘壩在向下游遷移的過(guò)程中，其中下部由于兩側(cè)水流往往具有不同的流態(tài)，心灘壩靠近主水道一側(cè)受水流沖刷作用強(qiáng)，靠近次水道一側(cè)受水流沖刷作用弱，沉積作用以側(cè)向加積為主，一側(cè)侵蝕，一側(cè)堆積。

圖4 加拿大South Saskatchewan River心灘壩向下游遷移現(xiàn)象紅線代表心灘壩于2012.8.28的位置，藍(lán)線代表心灘壩于2013.10.10的位置，白色箭頭代表水流方向Fig.4 The phenomenon of bar migration in South Saskatchewan River, Canada

圖5 心灘壩向下游遷移Fig.5 The phenomenon of bar migration

圖6 心灘壩遷移距離、速率與演化時(shí)間關(guān)系圖Fig.6 The relationship between evolution and migration distance, migration rate， respectively

以心灘壩A為例，其下游方向兩側(cè)水流速度不同，統(tǒng)計(jì)剖面位置不同時(shí)間步長(zhǎng)的流速，下方辮狀水道1最大流速為1.62 m/s，最小流速為1.44 m/s，平均流速為1.53 m/s，上方辮狀水道2最大流速為1.58 m/s，最小流速為1.37 m/s，平均流速為1.475 m/s。辮狀水道1流速較辮狀水道2流速大(圖8中的箭頭代表流速，流速越大，箭頭越大)，兩側(cè)水流為不對(duì)稱水流(圖8)，辮狀水道1為主水道，辮狀水道2為次水道，辮狀水道1侵蝕心灘壩A，而在辮狀水道2一側(cè)沉積物堆積(圖9)。以初始沉積厚度為基準(zhǔn)，主水道一側(cè)，在沉積模擬140個(gè)步長(zhǎng)至150個(gè)步長(zhǎng)的時(shí)間范圍內(nèi)，心灘壩A側(cè)向侵蝕80 m；次水道一側(cè)，側(cè)向沉積35 m，侵蝕速率大于沉積速率，心灘壩橫向?qū)挾茸冋Ｓ《鹊暮愫右泊嬖谶@種明顯由于不對(duì)稱水流引起的側(cè)向加積作用而導(dǎo)致心灘壩側(cè)向遷移的現(xiàn)象(圖10，圖中紅色輪廓代表先期心灘壩位置)。

圖9 心灘壩A的側(cè)向遷移(圖a剖面位置為圖b黃線位置)Fig.9 Lateral migration of bar A

圖10 印度恒河心灘壩由于側(cè)向加積作用造成的側(cè)向遷移現(xiàn)象Fig.10 Lateral accretion in Ganges, India

2.3 壩尾沉積物卸載導(dǎo)致心灘壩復(fù)合

當(dāng)辮狀水道完全越過(guò)心灘壩后，辮狀水道流速降低，在心灘壩頭部及主水道一側(cè)侵蝕的沉積物在心灘壩尾部卸載，發(fā)育“壩尾沉積”。由于主水道水動(dòng)力及輸砂能力較強(qiáng)，主水道一側(cè)“壩尾沉積”更發(fā)育，沉積速度更快，隨著“壩尾沉積”的不斷延長(zhǎng)，上游心灘壩與下游心灘壩之間的辮狀水道被充填，復(fù)合心灘壩形成。為更好的觀察復(fù)合心灘壩的形成過(guò)程，以5個(gè)步長(zhǎng)為單位繪制沉積厚度演化圖(圖11)，心灘壩E發(fā)育“壩尾沉積”，“壩尾沉積”延長(zhǎng)，心灘壩E與F之間的辮狀水道規(guī)?？s小，在150個(gè)時(shí)間步長(zhǎng)時(shí)，辮狀水道被完全充填，復(fù)合心灘壩(E+F)形成。類似的壩尾沉積(圖12a，b)以及與之相關(guān)的復(fù)合心灘壩(圖12c，d)的形成現(xiàn)象廣泛存在于現(xiàn)代沉積中。

值得注意的是，復(fù)合心灘壩的形成不僅取決于壩尾沉積，也與心灘壩的位置有關(guān)。當(dāng)兩個(gè)心灘壩之間距離較遠(yuǎn)，即它們之間的辮狀水道寬度大且穩(wěn)定時(shí)，難以形成穩(wěn)定的復(fù)合心灘壩，即使復(fù)合心灘壩形成之后，其連接處也易在演化過(guò)程中形成新的辮狀水道。

3 心灘壩發(fā)育演化模式

砂質(zhì)辮狀河中沉積物向下游方向搬運(yùn)過(guò)程中，上游底形凸起處逐漸形成單一心灘壩，心灘壩兩側(cè)遭受辮狀水道侵蝕，向下游遷移，遷移速度隨演化時(shí)間的增長(zhǎng)而逐漸減小(圖13a，f)。其向下游遷移的同時(shí)，心灘壩中下部?jī)蓚?cè)由于水道不對(duì)稱性，一側(cè)侵蝕，一側(cè)沉積(圖13a，f)。下游方向，沉積物堆積在心灘壩尾部，出現(xiàn)“壩尾沉積”(圖13c)，隨著“壩尾沉積”的不斷延長(zhǎng)，上游心灘壩與下游心灘壩復(fù)合形成復(fù)合心灘壩，復(fù)合心灘壩之間的辮狀水道被充填(圖13d)。當(dāng)壩頂接近水平面的同時(shí)，心灘壩基本穩(wěn)定，其頂部低洼處水流匯聚，水動(dòng)力增強(qiáng)，侵蝕其頂部，逐漸形成串溝，最后演變?yōu)檗p狀水道，切割心灘壩體(圖13d)。

單一心灘壩存在于心灘壩形成演化早期，上游侵蝕，下游沉積，主水道一側(cè)侵蝕，次水道一側(cè)沉積，因此上游方向心灘壩坡度陡，下游方向緩，主水道一側(cè)陡，次水道一側(cè)緩。隨演化時(shí)間增長(zhǎng)，砂質(zhì)辮狀河中由于心灘壩“壩尾沉積”的逐漸延長(zhǎng)，以穩(wěn)定的復(fù)合心灘壩為主，復(fù)合心灘壩具有與單一心灘壩一致的上游陡、下游緩，主水道一側(cè)陡，次水道一側(cè)緩的特征。當(dāng)單一心灘壩與復(fù)合心灘壩壩頂位置接近水平面時(shí)，壩頂辮狀水道的形成以及后期的消亡將心灘壩內(nèi)部結(jié)構(gòu)復(fù)雜化。

4 結(jié)論

本文通過(guò)沉積數(shù)值模擬與大型砂質(zhì)辮狀河現(xiàn)代沉積分析，明確了心灘壩的動(dòng)態(tài)演化過(guò)程，主要認(rèn)識(shí)包括以下四個(gè)方面：

(1) 心灘壩發(fā)育演化包括三個(gè)部分，心灘壩的形成、生長(zhǎng)及遷移，心灘壩的側(cè)向加積，“壩尾沉積”及復(fù)合心灘壩的形成。自心灘壩形成之后，三個(gè)發(fā)育演化過(guò)程無(wú)時(shí)間先后次序，各個(gè)過(guò)程所持續(xù)的時(shí)間取決于心灘壩兩側(cè)辮狀水道水動(dòng)力強(qiáng)度、水深以及下游心灘壩的位置。

圖11 復(fù)合心灘壩的形成Fig.11 The formation of compound bar

圖12 雅魯藏布江“壩尾沉積”及復(fù)合心灘壩形成現(xiàn)象a，b.為心灘壩發(fā)育的“壩尾沉積”，主辮狀水道“壩尾沉積”更發(fā)育；c，d.為由于“壩尾沉積”的發(fā)育形成復(fù)合心灘壩Fig.12 Bar tail and the formation of compound bar in Brahmaputra River

圖13 心灘壩演化模式圖Fig.13 Evolution pattern of bar

(2) 沉積物在局部高地堆積，逐漸生長(zhǎng)成為單一心灘壩，當(dāng)壩體生長(zhǎng)到一定規(guī)模之后，心灘壩頭部遭受兩側(cè)辮狀水道的沖刷作用，逐漸侵蝕，表現(xiàn)為心灘壩向下游方向遷移。遷移速率隨演化時(shí)間增長(zhǎng)而減慢，當(dāng)心灘壩頂部接近水平面時(shí)，向下游遷移速率基本為0，心灘壩穩(wěn)定，同時(shí)，水流在心灘壩頂部低洼處匯聚，沖刷其頂部，形成新的辮狀水道。

(3) 心灘壩向下游遷移的同時(shí)，中下部由于兩側(cè)辮狀水道的不對(duì)稱性，主水道一側(cè)沖刷作用強(qiáng)，心灘壩遭受侵蝕；次水道一側(cè)沖刷作用弱，沉積物在此側(cè)堆積，為側(cè)向加積作用。

(4) 沉積物在心灘壩尾部卸載，發(fā)育“壩尾沉積”，其發(fā)育程度取決于辮狀水道水動(dòng)力強(qiáng)度，主水道一側(cè)“壩尾沉積”發(fā)育程度高，隨著“壩尾沉積”的延長(zhǎng)，心灘壩之間的辮狀水道被充填，復(fù)合心灘壩形成。

)

[1] 裘亦楠. 河流沉積學(xué)中的河型分類[J]. 石油勘探與開(kāi)發(fā)，1985，12(2)：72- 74.[Qiu Yinan. Classification of river sedimentology[J]. Petroleum Exploration and Development, 1985, 12(2): 72- 74.]

[2] 錢(qián)寧. 關(guān)于河流分類及成因問(wèn)題的討論[J]. 地理學(xué)報(bào)，1985，40(1)：1- 10. [Qian Ning. On the classification and causes of formation of different channel patterns[J]. Acta Geographica Sinica, 1985, 40(1): 1- 10.]

[3] 何鯉，舒文震. 長(zhǎng)江上游的心灘：重慶珊瑚壩現(xiàn)代沉積考察[J]. 沉積學(xué)報(bào)，1986，4(1)：118- 125. [He Li, Shu Wenzhen. A channel bar in the upper reaches of the Changjiang river- research on recent sediment of Shanhuba in Chongqing, Sichuan province[J]. Acta Sedimentologica Sinica, 1986, 4(1): 118- 125.]

[4] 李維鋒，肖傳桃，王振奇，等. 湖北荊門(mén)?；蹨仙先B統(tǒng)河流相[J]. 石油與天然氣地質(zhì)，1993，14(4)：340- 345. [Li Weifeng, Xiao Chuantao, Wang Zhenqi, et al. Fluvial facies of upper Triassic at Haihuigou, Jingmen, Hubei[J]. Oil & Gas Geology, 1993, 14(4): 340- 345.]

[5] 王俊玲，任紀(jì)舜. 嫩江下游現(xiàn)代河流沉積特征[J]. 地質(zhì)論評(píng)，2001，47(2)：193- 199. [Wang Junling, Ren Jishun. Characteristics of modern fluvial deposits in the lower reaches of the Nenjiang River, Northeast China[J]. Geological Review, 2001, 47(2): 193- 199.]

[6] Kleinhans M G, Van Der Berg J H. River channel and bar patterns explained and predicted by an empirical and a physics- based method[J]. Earth Surface Processes and Landforms, 2011, 36(6): 721- 738.

[7] 張昌民，張尚鋒，李少華，等. 中國(guó)河流沉積學(xué)研究20年[J]. 沉積學(xué)報(bào)，2004，22(2)：183- 192. [Zhang Changmin, Zhang Shangfeng, Li Shaohua, et al. Advances in Chinese fluvial sedimentology from 1983 to 2003[J]. Acta Sedimentologica Sinica, 2004, 22(2): 183- 192.]

[8] 蘭朝利，何順利，門(mén)成全. 利用巖心或露頭的交錯(cuò)層組厚度預(yù)測(cè)辮狀河河道帶寬度：以鄂爾多斯盆地蘇里格氣田為例[J]. 油氣地質(zhì)與采收率，2005，12(2)：16- 18. [Lan Chaoli, He Shunli, Men Chengquan. Prediction of braided channel belt width based on cross- stratum sets thickness measurements of cores or outcrops- Taking Sulige gas field, Ordos Basin as an example[J]. Petroleum Geology and Recovery Efficiency, 2005, 12(2): 16- 18.]

[9] 朱筱敏，劉成林. 蘇里格地區(qū)上古生界有效儲(chǔ)層的確定[J]. 天然氣工業(yè)，2006，26(9)：1- 3. [Zhu Xiaomin, Liu Chenglin. Identification of effective upper Paleozoic reservoirs in Sulige area[J]. Natural Gas Industry, 2006, 26(9): 1- 3.]

[10] 侯加根，劉鈺銘，徐芳，等. 黃驊坳陷孔店油田新近系館陶組辮狀河砂體構(gòu)型及含油氣性差異成因[J]. 古地理學(xué)報(bào)，2008，10(5)：459- 464. [Hou Jiagen, Liu Yuming, Xu Fang, et al. Architecture of braided fluvial sandbody and origin for petroliferous difference of the Guantao Formation in Kongdian oilfield of Huanghua depression[J]. Journal of Palaeogeography, 2008, 10(5): 459- 464.]

[11] 趙翰卿，付志國(guó)，劉波. 應(yīng)用精細(xì)地質(zhì)研究準(zhǔn)確鑒別古代河流砂體[J]. 石油勘探與開(kāi)發(fā)，1995，22(2)：68- 70. [Zhao Hanqing, Fu Zhiguo, Liu Bo. Identification of the paleochannel sandbodies based on detailed geological study[J]. Petroleum Exploration and Development, 1995, 22(2): 68- 70.]

[12] Hjellbakk A. Facies and fluvial architecture of a high- energy braided river: the Upper Proterozoic Seglodden Member, Varanger Peninsula, northern Norway[J]. Sedimentary Geology, 1997, 114(1/2/3/4): 131- 161, 145- 149, 153- 161.

[13] Skelly R L, Bristow C S, Ethridge F G. Architecture of channel- belt deposits in an aggrading shallow sandbed braided river: the lower Niobrara River, northeast Nebraska[J]. Sedimentary Geology, 2003, 158(3/4): 249- 270.

[14] 張昌民，尹太舉，趙磊，等. 辮狀河儲(chǔ)層內(nèi)部建筑結(jié)構(gòu)分析[J]. 地質(zhì)科技情報(bào)，2013，32(4)：7- 13. [Zhang Changmin, Yin Taiju, Zhao Lei, et al. Reservoir architectural analysis of braided channel[J]. Geological Science and Technology Information, 2013, 32(4): 7- 13.]

[15] 金振奎，楊有星，尚建林，等. 辮狀河砂體構(gòu)型及定量參數(shù)研究:以阜康、柳林和延安地區(qū)辮狀河露頭為例[J]. 天然氣地球科學(xué)，2014，25(3)：311- 317. [Jin Zhenkui, Yang Youxing, Shang Jianlin, et al. Sandbody architecture and quantitative parameters of single channel sandbodies of braided river: Cases from outcrops of braided river in Fukang, Liulin and Yanan areas[J]. Natural Gas Geoscience, 2014, 25(3): 311- 317.]

[16] 陳玉琨，吳勝和，王延杰，等. 常年流水型砂質(zhì)辮狀河心灘壩內(nèi)部落淤層展布樣式探討[J]. 沉積與特提斯地質(zhì)，2015，35(1)：96- 101. [Chen Yukun, Wu Shenghe, Wang Yanjie, et al. Distribution patterns of the fall- silt seams in the channel bar of the perennial sandy braided river: An approach[J]. Sedimentary Geology and Tethyan Geology, 2015, 35(1): 96- 101.]

[17] 廖保方，張為民，李列，等. 辮狀河現(xiàn)代沉積研究與相模式:中國(guó)永定河剖析[J]. 沉積學(xué)報(bào)，1998，16(1)：34- 39. [Liao Baofang, Zhang Weimin, Li Lie, et al. Study on modern deposit of a braided stream and facies model- taking the Yongding River as an example[J]. Acta Sedimentologica Sinica, 1998, 16(1): 34- 39.]

[18] Best J, Woodward J, Ashworth P, et al. Bar- top hollows: a new element in the architecture of sandy braided rivers[J]. Sedimentary Geology, 2006, 190(1/2/3/4): 241- 255.

[19] Ashworth P J, Best J L, Roden J E, et al. Morphological evolution and dynamics of a large, sand braid- bar, Jamuna River, Bangladesh[J]. Sedimentology, 2000, 47(3): 533- 555.

[20] Schuurman F, Marra W A, Kleinhans M G. Physics- based modeling of large braided sand- bed rivers: Bar pattern formation, dynamics, and sensitivity[J]. Journal of Geophysical Research: Earth Surface, 2013, 118(4): 2509- 2527.

[21] Schuurman F, Kleinhans M G. Bar dynamics and bifurcation evolution in a modelled braided sand- bed river[J]. Earth Surface Process and Landforms, 2015, 40(10): 1318- 1333.

[22] Lesser G R, Roelvink J A, Van Kester J A T M, et al. Development and validation of a three- dimensional morphological model[J]. Coastal Engineering, 2004, 51(8/9): 883- 915.

[23] Roelvink J A. Coastal morphodynamic evolution techniques[J]. Coastal Engineering, 2006, 53(2/3): 277- 287.

[24] Van Maren D S. Grain size and sediment concentration effects on channel patterns of silt- laden rivers[J]. Sedimentary Geology, 2007, 202(1/2): 297- 316.

[25] Van Der Wegen M, Roelvink J A. Long- term morphodynamic evolution of a tidal embayment using a two- dimensional, process- based model[J]. Journal of Geophysical Research: Oceans, 2008, 113(C3): C03016.

[26] Crosato A, Saleh M S. Numerical study on the effects of floodplain vegetation on river planform style[J]. Earth Surface Processes and Landforms, 2011, 36(6): 711- 720.

[27] Crosato A, Mosselman E, Desta F B, et al. Experimental and numerical evidence for intrinsic nonmigrating bars in alluvial channels[J]. Water Resources Research, 2011, 47(3): 77- 79.

[28] Crosato A, Desta F B, Cornelisse J, et al. Experimental and numerical findings on the long- term evolution of migrating alternate bars in alluvial channels[J]. Water Resources Research, 2012, 48(6): W06524.

[29] 王楊君，尹太舉，鄧智浩，等. 水動(dòng)力數(shù)值模擬的河控三角洲分支河道演化研究[J]. 地質(zhì)科技情報(bào)，2016，35(1)：44- 52. [Wang Yangjun, Yin Taiju, Deng Zhihao, et al. Terminal distributary channels in fluvial- dominated delta systems from numerical simulation of hydrodynamics[J]. Geological Science and Technology Information, 2016, 35(1): 44- 52.]

[30] 錢(qián)寧，張仁，周志德. 河床演變學(xué)[M]. 北京：科學(xué)出版社，1987：160- 220. [Qian Ning, Zhang Ren, Zhou Zhide. Fluvial processes study[M]. Beijing: Science Press, 1987: 160- 220.]

[31] 于興河，馬興祥，穆龍新，等. 辮狀河儲(chǔ)層地質(zhì)模式及層次界面分析[M]. 北京：石油工業(yè)出版社，2004： 56- 58. [Yu Xinghe, Ma Xingxiang, Mu Longxin, et al. Analysis of braide driver reservoir geological model and hierarchy interface[M]. Beijing: Petroleum Industry Press, 2004: 56- 58.]