砂質辮狀河沉積演化機制與沉積構型模式

第一作者簡介，男，1988年生，副教授，碩士生導師，油氣田開發地質學、儲層沉積學，E-mail： fwj1017@yangtzeu.edu.cn

中圖分類號P512.2文獻標志碼A

DOI：10.14027/j.issn.1000-0550.2024.068

CSTR： 32268.14/j.cjxb.62-1038.2024.068

0 引言

砂質辮狀河是一種常見的河流類型，廣泛見于太古代至今的所有地層記錄中1-3]，可形成大規模、寬帶狀、泛連通的優質儲集體。在我國鄂爾多斯盆地、準噶爾盆地、塔里木盆地、松遼盆地、渤海灣盆地及國外多個含油氣盆地內均發現了大規模砂質辮狀河成因的油氣儲層4。大量開發實踐表明：砂質辮狀河儲層雖然規模大且整體泛連通[12]，但其內部構型特征復雜、非均質性強，存在多個層次、多種類型的高滲條帶和夾層[13-15]，導致砂質辮狀河儲層在開發初期往往產量較高，但進入注水開發階段后普遍存在見水快、含水率上升快、開發措施調整效果差等問題，到開發后期，含水率高的同時，大量的剩余油氣仍富集于儲層內部，嚴重制約了砂質辮狀河儲層高效開發[16-21]?；诂F代沉積與露頭觀察解剖、水槽模擬、沉積數值模擬等技術的研究已經表明，辮狀水道的頻繁遷移活動驅動心灘的反復形成、改造、重塑過程，導致辨狀河沉積構型十分復雜，并體現為心灘形態與分布樣式多樣，心灘內部增生體形態各異且復合關系復雜、辮狀水道不同程度侵蝕心灘并形成復雜的空間結構樣式[10.132233]。大量的研究表明，砂質辮狀河儲層內部不同級次的構型界面與構型單元控制的儲層內部單元疊切關系、連通性及儲層質量差異是主導地下流體運動特性與剩余油富集的關鍵因素之—[16-20]。因此，充分考慮砂質辮狀河形成與演化過程的復雜性并建立貼近地下儲層實際的沉積構型模式對于深化辮狀河沉積學理論并指導地下儲層構型解剖與剩余油氣資源挖潛具有重要意義。

長期以來基于現代沉積觀察與探地雷達淺層探測研究及水槽與沉積數值模擬研究主要聚焦于砂質辮狀河表層沉積地貌特征與沉積樣式研究，構建了較為豐富的沉積模式并獲得了大量的構型單元定量規模數據[2426.340]。然而，由于砂質辮狀河水道網絡的高頻活動和強烈的侵蝕、改造、沉積作用，最終保存在地層記錄中的辮狀河沉積構型及其構型單元規模與定量關系明顯有別于尚處于活躍狀態的現代辮狀河淺表部位。差異主要體現為活躍沉積最終保存程度較低（ lt;50% ）心灘保存比例較低、沉積構型要素保存不完整且多以侵蝕改造后的殘余沉積存在[30.41-43]。因而，從現代沉積中觀察到的心灘、辮狀水道難以完整地保存在地下儲層中，基于活躍辮狀河淺表沉積特征或沉積模擬過程中的地貌特征構建的沉積模式及獲取的定量規模信息難以有效地指導地下儲層表征。相比于現代沉積及模擬沉積地貌，野外露頭更接近地下儲層，然而由于露頭往往出露條件有限、系統性差且形成背景復雜，難以充分支撐儲層構型模式與定量規模研究[3，25，448]。

為了建立貼近地下砂質辨狀河儲層的構型模式，并建立構型單元定量規模與關系模型，需從長期演化過程分析的角度出發，明確演化過程中辨狀水道與心灘的遷移演化模式與侵蝕改造及保存規律，據此建立構型模式，并提取已保存在沉積記錄中的沉積構型要素規模以構建可靠的定量規模與相關模型。本文參考現代辮狀河沉積背景條件設定水槽沉積模擬裝置和邊界條件并開展砂質辮狀河水槽模擬實驗研究。實驗過程中，采用三維激光掃描儀按固定時間間隔獲取實驗沉積地貌數據，并據此重構三維數字化沉積構型模型，以分析實驗辮狀河長期連續演化過程及其最終保存的沉積構型特征及定量規模。通過細致的沉積演化過程分析，明確了砂質辨狀河形成過程、水道與心灘協同演化模式、沉積構型模式及構型單元定量規模，為地下砂質辮狀河儲層表征提供可靠的模式指導和定量規模與定量關系約束。

1實驗設計

1.1實驗裝置

為模擬砂質辮狀河的形成與演化過程，設置了一套水槽模擬裝置（圖1），該裝置由砂質底床、循環供水泵、螺桿式定速供砂器、水位控制邊界、儲水槽、三維激光掃描儀及延時攝影相機構成。砂質底床長8.5m ，寬 1.1m ，以天然河砂鋪設并輕微壓實，提供了足夠的模擬空間（圖2），底床高差 5cm 。模擬過程中，循環供水泵從底床末端水槽抽水并泵送到水砂混合漏斗中，與螺桿式定速供砂器提供的砂混合，從水槽底床左側邊界進入砂質底床發生沉積。水位控制邊界限定了砂質底床末端水位高低，以保證模擬過程中不因水位波動影響實驗辮狀河沉積。為了獲取沉積演化過程中的沉積地貌數據，在多個站點利用三維激光掃描儀掃測實驗辮狀河地貌，獲取高精度、高密度三維點云數據。延時攝影相機從多個視角以1s為間隔拍攝實驗辮狀河圖像。

1.2邊界條件

參考現代辮狀河沉積背景條件，設定了模擬邊界條件（表1）。底床成分與模擬中供應的砂質均為天然河砂，其粒度中值為 0.75mm ，兩者成分相似保證模擬過程中不因砂質沉積物成分與粒度對河流沉積演化造成干擾。砂質底床坡降為 0.05m ，坡長 8.5m 坡度約為5.9，經多輪預實驗測試，該坡度適于砂質辮狀河模擬。在模擬過程中，供水量保持在 1.5L/s 實驗實測水流速度大致為 0.1m/s ，經預實驗測試后，將供砂率設定為 3g/s ，以確保模擬過程中上游邊界處不發生砂質堆積或異常下切。為保證激光掃描儀掃測精度，以1h為間隔開展實驗，每個期次結束后排干底床內的明水，開展多站聯合激光掃描，記錄沉積過程中的地貌特征及其演變過程。本實驗共進行48個周期的模擬，在平緩的砂質底床上再現了一條典型砂質辮狀河逐步形成并持續演化的過程。從實驗時長與穩定邊界條件的前提出發，本實驗模擬的辮狀河可類比自然界沖積平原上的單一辮狀河，在地下儲層中可對應單砂層內部的辮狀河單砂體儲層。因而，從油氣藏開發的角度出發，實驗辮狀河構型特征可對應單砂體內部構型。

1.3 實驗數據處理方法

實驗過程中采集的激光點云數據是具有三維坐標、圖像信息的離散點云，無法直接用于沉積學分析和構型解剖。實驗數據的處理分為三個步驟，首先，采用Scene軟件對單期次多站點云數據進行融合處理，形成坐標系統一的點云數據；隨后，提取實驗辮狀河區域內點云，并進行點云數據優選和網格化，形成各期次模擬后的數字化沉積地貌數據與正射投影圖像數據；進一步地，對沉積地貌數據集進行去趨勢處理，排除底床順源高程差異，突出辨狀河心灘與辮狀水道形態、規模信息并用于沉積演化過程分析；最后，參考Straubetal.4提出的數字地層重構方法，開展沉積演化過程回溯，計算經過反復侵蝕改造后最終保存在沉積記錄中的期次邊界，并據此重構模擬辨狀河三維空間結構模型，該模型用于沉積構型模式研究與沉積構型單元定量規模測量分析。

2辮狀河形成與持續演化特征

2.1實驗辮狀河形成過程

在實驗初期（第1＼～16模擬期次），水流攜帶沉積物進人砂質平緩的砂質底床上發生沉積并逐步演化形成一條典型的砂質辮狀河（圖3）。在辮狀河的逐步形成過程中，砂質沉積物主要以底床搬運方式順源推進，首先形成初始的心灘和初步的辮狀水道分流（圖3a、圖4），隨著砂質沉積物向下游遷移，心灘數量增加，辮狀水道流程延長、數量增加、彎曲度變大（圖3b＼～c），這一過程中，已形成的上游心灘也在快速演化。以第9期形成的心灘Bar-A為例，在經歷3個周期的模擬后，該心灘已經被后續的辮狀水道侵蝕切割，由半對稱的大規模復合心灘轉變為多個近對稱的小規模心灘，辮狀水道由單側環繞轉變為多條交織環繞（圖3c＼～d），到第12期模擬結束，該心灘已發生明顯的順源遷移，形成大規模復合心灘。基于辮狀河初始形成過程的觀察表明，底床搬運是砂質辮狀河形成過程中的基本動力機制，這一機制決定了砂質辮狀河整體水淺流急，無法形成長期穩定的水道，其演化總是伴隨著底床沉積物遷移堆積、心灘增生、辮狀水道遷移等復雜的遷移變化過程，且河網與心灘演變速度快、頻次高。

需要注意的是，初始形成的心灘主要由上游辮狀河道集中供給而成，整體呈現朵狀特征（圖3，4），與后期遭受辮狀水道改造后的心灘存在較大差異。典型的初始心灘Bar-B（圖4）在形成過程中先由上游辮狀水道供給形成朵狀雛形并持續順源增大（圖4a1＼～a2，b1＼～b2），隨后，由于心灘頂部堆積沉積物出露水面，心灘增生轉向側緣（圖 4a3～a4，b3～b4， ，最后該心灘發生退積和側積，轉而在心灘上游及側緣發生增生（圖 4a5～a6，b5～b6） 。這類初始心灘在隨后的沉積演化過程中往往被部分破壞，并最后殘存于辮狀河底部。由于心灘形成過程中存在進積、退積及側積，其內部構型較為特殊，有別于辮狀河演化成熟階段形成的心灘。

2.2辮狀水道與心灘連續協同演化機制

不同于辮狀河形成初期沉積物持續順源推進的特征，一旦辮狀河形成并進入持續演化期，辮狀水道與心灘整體樣式大致維持穩定，但辮狀水道與心灘的協同演化仍連續進行，表現為辮狀水道瀕繁遷移演化，心灘也隨之逐步發生增生、遭受改造甚至被徹底侵蝕，從而形成十分復雜的沉積構型，多數心灘和辮狀水道沉積難以被完整保存。因而，明確辮狀水道與心灘協同演化模式，是深人揭示地下儲層沉積構型的必要前提，據分析，辮狀水道與心灘的協同演化機制主要包括三類。

2.2.1 辮狀水道連續側向遷移主導心灘側積增生

辮狀水道連續側向遷移現象常見于心灘側緣，主要位于水道彎曲部位，由上游供水與心灘表面水流匯人共同控制形成二次流（secondaryflow），導致在局部河段發生水道單向環流，引發辮狀水道局部段發生凸岸侵蝕和凹岸堆積。在實驗過程中辮狀水道連續側積現象十分常見，并主導心灘的側向增生。以第15＼～19期次模擬過程中局部河段為例（圖5），心灘Bar-C側緣較寬的辮狀水道Ch-A受到心灘分布控制發生連續側積，侵蝕該水道凸岸心灘，并引發凹岸心灘不斷側積，形成側積尾翼（圖5）。在辮狀水道連續側向遷移過程中，其上游段整體維持穩定，但下游快速掃動（圖5a＼～c）并最終取直改道（圖5d），在這一過程中側積尾翼甚至充填了因下游河段快速遷移而廢棄的河道（圖5d）。

辨狀水道連續側積是其一側心灘遭受侵蝕而另一側心灘發生側積的主導因素，這一現象多發生在心灘中游側緣。在連續側向遷移過程中，辮狀水道上游基本維持穩定而下游快速遷移，最終導致下游取直改道，心灘尾部也往往因下游水道遷移廢棄的河道被快速充填而與更下游心灘復合（圖5d）。

不同于曲流河側向遷移過程中河道本身的高度穩定特征，辮狀水道的側積往往導致其下游河段發生快速的遷移和廢棄一充填沉積。受辮狀水道連續側積影響，其凸岸一側心灘多遭受深切而難以完整保存（圖5c＼～d）;而其凹岸一側形成的增生體拼合于原有心灘側緣，且新形成的增生體側積方向明顯有別于原有心灘內部增生體（圖5d）。

2.2.2辮狀水道廢棄充填并與心灘拼接復合

在辮狀水道網絡中，局部河段的變遷易引起下游河段廢棄，形成廢棄充填沉積并與鄰近心灘拼接復合，從而導致心灘迅速增生。以第20＼～21期次模擬過程中局部河段為例（圖6），在第20期模擬結束后心灘Bar-D一側較大規模狀水道Ch-B因其上游水道變遷而流量銳減（圖6a），在第21期次模擬過程中，該辮狀水道被迅速充填并與心灘Bar-D拼合形成更大規模的復合心灘（圖6b）。

大規?；蛑鞲尚娃p狀水道廢棄后，其下游多條水道也可能同步廢棄。如Ch-B廢棄地充填的同時，還存在多條同步廢棄的水道，這類水道主要是心灘中下游中央部位具有匯流特征的小規模水道或水道組合，其主要成因是心灘表面沉積物順水流向下游搬運沉積，這類廢棄水道充填造成心灘順流增生或加積，主要影響心灘中下游中央部位厚度（圖6c），這類廢棄河道充填廣泛存在于心灘中下游部位，圖5所示的連續遷移辮狀水道下游也存在多條這類廢棄河道充填沉積。大規?；蛑鞲尚秃拥缽U棄過程中，也可能形成小規模、短時存續的取直水道（Ch-D）（圖6b，c），這類水道往往形成于大規模或主干型水道廢棄末期，難以持續穩定存在，在后續的沉積過程中易被充填。

2.2.3 辮狀水道匯水沖刷與心灘重構

在辮狀河沉積演化過程中，同一河段多條辮狀水道往往因相向遷移而匯水，由于河床內水流幾乎完全匯集于一條大規模水道內，水動力驟然增強，可大面積侵蝕其下游心灘，并在短時間內局部或整體重塑下游心灘形態（圖7）。

以第25＼～27期模擬過程中局部河段為例，在第25＼～26期次模擬中，心灘Bar-E、Bar-F僅受到辮狀水道遷移或取直沖刷輕微改造，主要的辮狀水道Ch-E和Ch-F維持穩定，心灘Bar-F受辮狀水道Ch-G和Ch-I的廢棄與取直下切而發生形變，但河段整體特征維持大體穩定（圖7a，b）。在第27期次模擬中，心灘Bar-E上游辮狀水道發生匯聚，形成大規模辮狀水道Ch-J并沖刷心灘Bar-E，導致該心灘大部區域被深切沖蝕。被侵蝕的沉積物順水流搬運至下游并導致心灘Bar-G的形成和心灘Bar-F的迎水生長（圖7c）。僅在1個模擬期次內，辮狀水道匯水沖刷導致心灘Bar-E被整體重塑，而心灘Bar-F發生局部重塑（圖7）。

對整個實驗的觀察結果顯示，辮狀水道匯水沖刷作用發生頻次較高、維持時間較短但對辮狀水道網絡與心灘的重塑作用強烈。距離匯水點越近，水流沖刷作用越強，可導致心灘主體被完全沖蝕，沖刷點下游方向則在短時間內迅速堆積被沖刷搬運的沉積物，因而下游河段辮狀水道數量增加、水深變淺（圖7c）。匯水沖刷作用的范圍通常在長1＼～2個復合心灘長度的河段內發生，并不直接影響下游河段沉積特征。因此，辮狀河沉積演化過程中，普遍存在局部河段快速演變，整體水道一心灘樣式維持穩定的特征。

2.3心灘沉積規模變化規律與保存特征

針對實驗過程的系統觀察表明，心灘沉積演化過程中普遍存在形成、生長、遭受侵蝕、殘存等過程，在這一過程中心灘規模隨之有規律地變化。基于去趨勢沉積地貌圖對實驗辮狀河沉積演化過程中存在的典型心灘進行連續追蹤與規模測量（圖8），取得以下結果：（1）單個心灘的形成過程較為迅速，通常在1＼～6個模擬期次內持續生長至最大規模，隨后在3＼～8個模擬期次內遭受連續侵蝕而規模減?。ū?），最終僅小部分得以最終保存或被完全侵蝕（圖8）；（2）部分心灘形成過程中經歷侵蝕、生長、再侵蝕，或同時遭受侵蝕和生長，其規模存在一定程度的波動，但其最終保存情況不變；（3）心灘長寬及面積變化大體保持一致，但其長寬比存在一定程度的波動，心灘長寬比平均為2.8，與現代辮狀河測量結果相近[24.35.50]，多數心灘在演化過程中維持穩定，但部分心灘在其演化末期受到較大幅度侵蝕后存在長寬比增大現象。針對11個典型心灘演化過程完整追蹤的測量結果顯示，心灘形成后受到辮狀水道的進一步侵蝕改造，僅18.66%～57.41% 的分布區域殘存于沉積體中（表2），平均最終保存率為 36.28% 。

3砂質辨狀河沉積構型特征

在辮狀水道與心灘協同演化模式分析的基礎上，利用序列化沉積地貌數據重構實驗辮狀河三維數字化沉積構型模型并提取典型部位構型剖面，再現了實驗后最終保存在沉積記錄中的砂質辮狀河沉積構型特征（圖9）。需要指出的是，水槽沉積模擬產生的辮狀河規模較小，按固定時間間隔進行三維激光掃描并重構的模擬辮狀河三維構型模型內主要包含了復合砂體、單砂體及單砂體內部增生體（Miall分級為5＼～3級）3個級次的構型單元信息。

3.1 辮狀水道沉積構型特征

根據辮狀水道形態、規模、疊切樣式等特征差異，實驗結束后保存在最終沉積記錄中的辮狀水道（4級構型單元）主要分為三類，包括復雜疊切水道（圖9，類型A）、下切水道（圖9，類型B）孤立水道（圖9，類型C）。

復雜疊切水道發育程度最高，由多期辮狀水道側向遷移沉積與垂向疊切而成的增生體（4級構型單元）構成，其寬度可達單一辮狀水道的2＼～10倍，平均寬 0.28m 深 0.12m 。復雜疊切水道內部一般存在3＼～15期水道疊切過程（圖9b1＼～b4），單期水道整體較厚，在橫切剖面上表現為高強度疊切，而在平面上則表現為多條流向不同、規模各異、河段長度參差的水道在某一河段反復遷移活動（圖7）。復雜疊切水道之間也存在側向疊切或拼接，導致其寬度遠大于剖面上的心灘沉積寬度（圖9b1＼～b4）。在順源方向上，復雜疊切水道延伸疊切強度較低（圖9c1＼～c2），其疊切過程與不同模擬期次內不同辮狀水道深度、遷移活動強度及方向相關。

下切水道多形成于辮狀河模擬的初期階段，由水流匯集深切后充填而成，侵蝕底床并隨后充填（圖9b2，b4，c1，c2，類型B）。下切水道僅局部存在，其深度與復雜疊切水道相當，但寬度明顯小于復雜疊切水道。

孤立水道在切物源剖面上發育程度較低，屬于短時存在的辮狀水道（圖9b2＼～b3，類型C），多由新形成的辮狀水道快速廢棄充填形成。

統計剖面上辮狀水道面積占比可知，以復雜疊切辮狀水道為主的水道沉積占比可達 57.9% ，這表明在最終保存的辨狀河沉積體內，河道占比高于心灘。

3.2 心灘沉積構型特征

剖面構型分析表明，實驗辮狀河內部心灘（4級構型單元）存在5方面特點：其一，心灘多被辮狀水道切割，規模較小，形態破碎（圖9），與現代沉積或實驗地貌中識別出的心灘寬度相比顯著偏??；其二，心灘一般由多個垂向加積體構成（3級構型單元），平均單個心灘內包含3＼～8個期次內形成的多個加積體，加積體頂底界面通常平緩（圖9），僅少量心灘內部增生體存在底平頂凸的特征（圖9b4）；其三，剖面上心灘沉積占比僅為 42.1% ，顯著低于現代辮狀河觀察結果；其四，辮狀河底部心灘往往連續性更好、規模更大，受辮狀水道侵蝕改造程度較低，而中上部心灘沉積比例偏低，且被辮狀水道侵蝕改造程度更高（圖9）；其五，在垂向上，心灘增生體之間往往存在密集發育的增生體界面，顯示該處存在連續多期薄層披覆泥巖沉積，為泥質夾層。

4構型單元及其內部增生體規模與定量關系

以重構的模擬辮狀河三維空間結構模型為基礎，截取切物源剖面并測量最終保存在實驗辮狀河沉積體內部的構型單元及增生體單元的定量規模，建立了構型單元及其內部增生體定量規律。

4.1主要構型單元及其內部增生體規模

以重構的多條構型剖面為基礎，測量相鄰的辨狀水道一心灘組合18個、心灘26個、辮狀水道33個、心灘內部增生體180個、水道內部123個，并按測量剖面統計其平均值（表3），結果表明：（1）辮狀水道沉積占比介于 44.81%～66.20% ，平均為 57.9% ，是辮狀河沉積的主體；（2）心灘沉積占比為 42.1% ，單一心灘平均寬 0.258m 厚 0.013m ，寬厚比為19.8；（3）單一心灘內平均發育5.7＼～12.4期增生體，增生體平均寬0.127m， 厚 0.005m ，寬厚比為25.4；（4）水道沉積平均寬 0.268m 厚 0.019m ，寬厚比為14.1；（5）單一水道內平均包含4.2期增生體，平均寬 0.137m. 厚0.010m ，寬厚比為13.7。

4.2 構型單元及其內部增生體定量關系

根據測量所得定量規模數據，明確了砂質辨狀河沉積體內部心灘、水道及兩者內部增生體規模的定量關系（圖10）。結果表明：（1）心灘與辮狀水道寬度存在良好的正相關關系，線性擬合優度為0.716（圖 10a ），兩者厚度也存在正相關關系，但由于部分水道垂向疊切或水平疊切拼合，其厚度變化幅度較大，線性擬合優度僅為0.341（圖10b）；（2）心灘、水道的寬度與厚度整體呈正相關關系，但由于部分心灘被侵蝕程度高、疊置期次較多、少數水道側向疊切拼接程度異常高，線性擬合優度僅為0.378、0.227（圖10c＼～d）；（3）心灘內部增生體的寬度、厚度也存在一定的正相關關系，部分增生體因遭受較強的側向侵蝕切割，殘存下來的增生體厚度較大、寬度較??；（4）水道內部增生體大致可分為3類，除主體樣本外，存在沉積動力弱、侵蝕改造強度較低的厚層側向增生體和沉積動力較強、侵蝕改造強度高的兩類規模異常的增生體。

（a）心灘與辮狀水道寬度相關性;（b）心灘與辮狀水道沉積厚度相關性；（c）心灘寬厚相關性;（d）辮狀水道沉積寬厚相關性； Π（e） 心灘內部增生體寬厚相關性：（f）辮 狀水道內部增生體寬厚相關性

（o width;（d）correlationofbraidchannelthicknessandwidth;（e）correlationofchannelbaraccretionthicknessandwidth; （f） correlation of braid channel accretion thickness and width

需要指出的是，水槽模擬形成的辮狀河規模較小，其沉積過程也受限于水槽測量固定邊界，因而測量得到的沉積構型單元定量規模與自然界辨狀河可能存在一定程度的差異，有必要尋找出露狀態較好的大規模露頭群開展對比驗證。

5結論

（1）砂質辮狀河的初期演化主要表現為砂質沉積從源頭向下游連續推進，形成初始心灘并逐漸發生水道分叉匯聚，初始心灘一般由上游水道集中供給形成，經過前積、退積及側積逐步形成，多呈朵狀，規模較大，通常僅發育于辮狀河底部，最終保存程度較高。

（2）在辮狀河形成后，辮狀水道與心灘連續協同演化，主要存在三種機制，包括： ① 辮狀水道連續側積主導心灘側積增生； ② 辮狀水道廢棄充填并與心灘拼接復合； ③ 辮狀水道匯水沖刷導致下游心灘重構。辮狀水道與心灘頻繁、快速且持續的協同演化導致心灘快速形成并隨后持續遭受后期水道侵蝕改造，心灘平均最終保存比例僅約 36.28%

（3）模擬實驗結束后，最終保存的砂質辮狀河沉積體內，辮狀水道占比約 57.9% ，心灘占比約 42.1% 。辮狀水道主要分為復雜疊切水道、下切水道及孤立水道，其中復雜疊切水道發育比例最高。心灘多遭受水道侵蝕改造，規模較小，形態破碎，顯著小于從沉積地貌觀察測量獲得的活躍心灘規模。心灘內部一般包含多個垂向加積體，其頂底界面多平緩。整體上，辮狀河沉積體下部心灘保存程度較高，中上部心灘保存程度低、規模小、與辮狀水道疊切關系復雜。

（4）辮狀水道網絡頻繁變遷主導下形成的砂質辨狀河沉積整體可視為一個泛連通砂體（5級），其內部由遭受復雜疊切過程而殘存的辮狀水道與心灘（4級）構成，其中，辮狀水道主要為側向疊切拼接或孤立分布，整體保存程度較高，形態較為完整，心灘則分布于辮狀水道之間，以殘存狀態有限保存。辮狀水道內部通常包含多期、多類增生體（3級），心灘內部通常包含多期垂向加積增生體（3級），垂向增生體之間往往存在披覆泥巖夾層。

（5）系統的規模測量表明，狀水道與心灘寬厚比為14.1和19.8，單一心灘內部發育5.7＼～12.4期增生體，增生體平均寬厚比為25.4，單一水道內部平均包含4.2期增生體，寬厚比為13.7。辮狀水道、心灘及其內部增生體寬厚具有明顯線性正相關關系，但辮狀水道側向拼接復合程度和辨狀水道侵蝕切割強度差異導致部分水道、心灘及增生體規模異常。

參考文獻（References）

[1]張昌民，朱銳，趙康，等．從端點走向連續：河流沉積模式研究 進展述評[J].沉積學報，2017，35（5）：926-944.[Zhang Changmin，ZhuRui，ZhaoKang，etal.From end memberto continuum： Review of fluvial facies model research[J].Acta Sedimentologica Sinica，2017，35（5）：926-944.]

[2]張昌民，張尚鋒，李少華，等．中國河流沉積學研究20年[J].沉 積學報，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.]

[3]張金亮.河流沉積相類型及相模式[J].新疆石油地質，2019，40 （2）：244-252.[Zhang Jinliang. Fluvial facies styles and their sedimentary facies models[J].Xinjiang Petroleum Geology，2019， 40 （2）： 244-252. ]

[4］張憲國，林承焰，張濤，等．大港灘海地區地震沉積學研究[J]. 石油勘探與開發，2011，38（1）：40-46.[ZhangXianguo，Lin Chengyan， Zhang Tao， et al. Seismic sedimentologic research in shallow sea areas，Dagang[J].Petroleum Exploration and Development，2011，38（1）：40-46.]

[5］張昌民，尹太舉，趙磊，等．辮狀河儲層內部建筑結構分析[J]. 地質科技情報，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.]

[6］文華國，鄭榮才，高紅燦，等．蘇里格氣田蘇6井區下石盒子組 盒8段沉積相特征[J].沉積學報，2007，25（1）：90-98.[Wen Huaguo， Zheng Rongcai， Gao Hongcan， et al. Sedimentary facies of the8th member of Lower Shihezi Formation in Su6 area， Suligegas field[J].ActaSedimentologica Sinica，2007，5（1）： 90-98.]

[7］譚程鵬，于興河，李勝利，等．辮狀河—曲流河轉換模式探討： 以準噶爾盆地南緣頭屯河組露頭為例[J]．沉積學報，2014，32 （3）：450-458.[Tan Chengpeng，Yu Xinghe，Li Shengli，etal.Discusson on the model of braided river transform to meandering river：As an example of Toutunhe Formation in southern Junggar Basin[J]. Acta Sedimentologica Sinica，2014，32（3）： 450-458.]

[8］賈愛林，唐俊偉，何東博，等．蘇里格氣田強非均質致密砂巖儲 層的地質建模[J]．中國石油勘探，2007，12（1）：12-16，92.[Jia Ailin，Tang Junwei，He Dongbo，et al.Geological modeling for sandstone reservoirs with low permeability and strong heterogeneity in Sulige Gasfield[J].China Petroleum Exploration，2007，12 （1）： 12-16，92. ]

[9]郭智，孫龍德，賈愛林，等．辮狀河相致密砂巖氣藏三維地質建 模[J]．石油勘探與開發，2015，42（1）：76-83.[Guo Zhi， Sun Longde，JiaAilin，etal.3Dgeologicalmodelingfortight sand gas reservoir of braided river facies[J].Petroleum Exploration and Development，2015， 42（1）： 76-83.]

[10］陳東陽，王峰，陳洪德，等．鄂爾多斯盆地東部府谷天生橋剖 面上古生界下石盒子組8段辮狀河儲層構型表征[J].石油與 天然氣地質，2019，40（2）：335-345.[Chen Dongyang，Wang Feng，Chen Hongde，et al.Characterization of braided river reservoir architecture of the Upper Paleozoic He 8 member on Fugu Tianshengqiao outcrop，eastern Ordos Basin[J].Oil amp; Gas Geology，2019，40（2）： 335-345.]

[11］陳仕臻，林承焰，任麗華，等．砂質辮狀河沉積模式的建立：以 委內瑞拉奧里諾科重油帶H區塊為例[J].沉積學報，2015，33 （5）：965-971.[Chen Shizhen，Lin Chengyan，Ren Lihua，et al. Establishmentof thedepositionalmodelof sandybraided river： A case from the H block in Orinoco Heavy Oil Belt，Venezuela [J].Acta Sedimentologica Sinica，2015，33（5）： 965-971.]

[12]葛云龍，逯徑鐵，廖保方，等.辮狀河相儲集層地質模型：\"泛 連通體\"連通體”石油勘探與開發，1998，25（5）：77-79，5-6， 13.[Ge Yunlong，Lu Jingtie，Liao Baofang，etal. A braided riverreservoir geological model： \"Pan communicated sandbody\" [J].Petroleum Exploration and Development，1998，25（5）：77- 79，5-6，13.]

[13]Li S L， Yu X H， Chen B T， et al. Quantitative characterization of architecture elements and their response to base-level change in a sandy braided fluvial system ata mountain front[J].Journal of SedimentaryResearch，2015，85（10）：1258-1274.

[14］劉鈺銘，侯加根，宋保全，等．辮狀河厚砂層內部夾層表征：以 大慶喇嘛甸油田為例[J].石油學報，2011，32（5）：836-841. [Liu Yuming，Hou JiagenSong Baoquan，etal.Characterization of interlayers within braided-river thick sandstones：A case study on the Lamadian oilfield in Daqing[J].Acta Petrolei Sinica， 2011，32（5）： 836-841. ]

[15］楊麗莎，陳彬滔，李順利，等．基于成因類型的砂質辮狀河泥 巖分布模式：以山西大同侏羅系砂質辮狀河露頭為例[J]．天然 氣地球科學，2013，24（1）：93-98.[Yang Lisha，Chen Bintao，Li Shunli，et al. Pattern of genesis-based mudstone distribution for sandy braided river： Acase studyof sandy braided-river outerop， Datong， Shanxi province， China[J]. Natural Gas Geoscience， 2013，24（1）： 93-98.]

[16］孫煥泉，孫國，程會明，等．勝坨油田特高含水期剩余油分布 仿真模型[J].石油勘探與開發，2002，29（3）：66-68.［Sun Huanquan，Sun Guo，Cheng Huiming，et al. The simulation models of remaining oil distribution at super-high water-cut stage of Shengtuo oil field[J].Petroleum Exploration and Development，2002，29（3）： 66-68.]

[17]杜啟振，侯加根，陸基孟．儲層微相及砂體預測模型[J].石油 學報，199，20（2）：45-50.[Du Qizhen，Hou Jiagen，LuJimeng. A predictable geologic model of sedimentary facies and sands [J].Acta Petrolei Sinica，1999，20（2）： 45-50.]

[18] 林承焰，孫廷彬，董春梅，等．基于單砂體的特高含水期剩余 油精細表征[J].石油學報，2013，34（6）：1131-1136.[Lin Chengyan，unTingbin，ongChunmei，etal.Finecharacterization of remaining oil basedon a single sand body in the high water cut Period[J]. Acta Petrolei Sinica，2013，34（6）：1131- 1136.]

[19］胡永樂，王燕靈，楊思玉，等．注水油田高含水后期開發技術 方針的調整[J].石油學報，2004，25（5）：65-69.[Hu Yongle， Wang Yanling，Yang Siyu， et al. Adjustment of technical policy forwater-flooding oilfieldwith highwatercutin the late stage of development[J].Acta Petrolei Sinica，2004，25（5）： 65-69.]

[20］計秉玉，李彥興．喇薩杏油田高含水期提高采收率的主要技 術對策[J].大慶石油地質與開發，2004，23（5）：47-53.[Ji Bingyu，Li Yanxing.Main technical countermeasuresof enhanced oil recovery during high water cut stage in Lasaxing reservoirs[J].Petroleum Geologyamp; Oilfield Development in Daqing，2004，23（5）： 47-53.]

[21］韓大匡．準確預測剩余油相對富集區提高油田注水采收率研 究[Jl.石油學報.2007，28（2）：73-78.「Han Dakuang.Preciselv predicting abundant remaining oil and improving the secondary recovery of mature oilfields[J]．Acta Petrolei Sinica，2007，28 （2）： 73-78. ]

[22]An H P， Chen S C，Chan HC，et al.Dimension and frequency of bar formation ina braided river[J]. International Journal of Sediment Research，2013，28（3）： 358-367.

[23]Egozi R，Ashmore P.Experimental analysis of braided channel pattern response to increased discharge[J]. Journal of Geophysical Research：Earth Surface，2009，114（F2）：F02012.

[24]Hundey E J， Ashmore PE.Length scale of braided river morphology[J].Water Resources Research，2009，45（8）： W08409.

[25]Li W， Yue D L，DuY S，et al.Controls of accommodation to sediment-supply ratio on sedimentary architecture of continental fluvial successions[J]. Petroleum Science，2023，20（4）：1961- 1977.

[26]LiW，ColomberaL，Yue DL，et al. Controls on the morphology of braided rivers and braid bars：An empirical characterization of numerical models[J]. Sedimentology，2023，70（1）： 259-279.

[27]Lisle TE，Ikeda H，Iseya F.Formation of stationaryalternate bars ina steep channel with mixed-size sediment： A flume experiment[J].Earth Surface Processes and Landforms，1991，16（5）： 463-469.

[28]Schuurman F，MarraWA，Kleinhans MG. 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.

[29]Surian N.Fluvial processes in braided rivers[M]//Rowinski P， Radecki-Pawlik A.Rivers-physical，fluvial and environmental processes. Cham： Springer，2015： 403-425.

[30] Van De LagewegWI， Van Dijk WM，Kleinhans MG.Morphologicaland stratigraphical signatureof floodsina braided gravelbed river revealed from flume experiments[J].Journal of Sedimentary Research，2013，83（11）： 1032-1045.

[31]Zhang K，Wu SH，FengWJ，et al.Bar dynamics ina sandy braided river： Insights from sediment numerical simulations[J]. Sedimentary Geology，2020，396：105557.

[32］張可，吳勝和，馮文杰，等．砂質辮狀河心灘壩的發育演化過 程探討：沉積數值模擬與現代沉積分析啟示[J].沉積學報， 2018，36（1）：81-91.[Zhang Ke，Wu Shenghe，Feng Wenjie， et al.Discussion on evolution of bar in sandy braided river：Insights from sediment numerical simulation and modern bar[J]. Acta Sedimentologica Sinica，2018，36（1）： 81-91.]

[33]雷濤，任廣磊，李曉慧，等．砂質辮狀河心灘沉積演化規律與 沉積構型特征：基于沉積數值模擬的認識[J].石油與天然氣地 質，2023，44（6）：1595-1608.[Lei Tao，Ren Guanglei，Li Xiaohui，et al.Sedimentary evolution pattern and architectural characteristics of mid-channel bars in sandy braided rivers：Understanding based on sedimentary numerical simulation[J].Oil amp; Gas Geology，2023，44（6）： 1595-1608.]

[34]Best JL，Ashworth PJ， Bristow C S，et al. Thre-dimensional sedimentaryarchitectureofa large，mid-channel sand braid bar， JamunaRiver，Bangladesh[J].Journal of SedimentaryResearch， 2003，73（4）： 516-530.

[35]Li HC，Li JH，Li Z D. Quantitative scale analysis of the channel bar ina braided river and its internal architecture[J].Applied Sciences，2023，14（1）：257.

[36]Mumpy AJ，Jol H M，Kean WF， et al. Architecture and sedimentology of anactive braid barin the WisconsinRiver based on3-D ground penetrating radar[M]//Baker G S，Jol HM.Stratigraphic analyses using GPR．Boulder：Geological Society of America，2007：111-131.

[37]Okazaki H，Kwak Y， Tamura T. Depositional and erosional architecturesof gravellybraidbarformedbya flood in theAbe River，central Japan， inferred froma three-dimensional groundpenetrating radar analysis[J].Sedimentary Geology，2015，324： 32-46.

[38]Parker NO，Sambrook Smith G H，Ashworth PJ，et al. Quantification of the relation between surface morphodynamics and subsurface sedimentological product in sandy braided rivers[J].Sedimentology，2013，60（3）： 820-839.

[39]Rice S P， Church M， Wooldridge C L，et al. Morphology and evolution of bars ina wandering gravel-bed river;lower Fraser river，British Columbia，Canada[J].Sedimentology，20o9，56（3）： 709-736.

[40]Sambrook Smith G H，Ashworth P J， Best JL，et al. The sedimentology and alluvial architecture of the sandy braided South Saskatchewan River， Canada[J]. Sedimentology，2006，53（2）： 413-434.

[41]Abdel-Fattah Z A.Fluvial architecture of the Upper Cretaceous Nubia Sandstones：Anancient example of sandy braided rivers in central eastern Desert，Egypt[J].Sedimentary Geology，2021， 420：105923.

[42]Paola C，Borgman L. Reconstructing random topography from preservedstratification[J].Sedimentology，1991，38（4）： 553-565.

[43]Sambrook SmithGH，Ashworth PJ，BestJL，et al.The sedimentologyandalluvialarchitectureofalargebraidbar，RioParana， Argentina[J].Journal of Sedimentary Research，20o9，79（8）： 629-642.

[44]Hossain S， Shekhar H，Rahman N. Facies and architectural element analysis of the Upper Bokabil Sandstone in the Bengal Basin[J]. Sedimentary Geology，2023，453：106433.

[45]Leleu S，HartleyAJ，Williams BPJ.Large-scale alluvial architectureand correlation in a Triassic pebbly braided river system， LowerWolfvilleFormation （FundyBasin，Nova Scotia， Canada） [J].Journal of Sedimentary Research，20o9，79（5）：265-286.

[46]PuigJM，CabellP， HowellJ，etal.Thre-dimensional characterisation ofsedimentary heterogeneityand itsimpacton subsurface flow behaviour through the braided-to-meandering fluvial deposits of the Castissent Formation （Late Ypresian，TrempGraus Basin， Spain）[J].Marine and Petroleum Geology， 2019， 103：661-680.

[47]Zang D S，Bao Z D，LiMY，et al.Sandbody architecture analysisofbraided river reservoirsand their significance for remaining oil distribution：A case study based on a new outcrop in the Songliao Basin，Northeast China[J].Energy Explorationamp; Exploitation，2020，38（6）： 2231-2251.

[48]Zhang CM， Yin TJ，Wu S H，et al. Architectural element analysis of nonmarine oil and gas reservoir in China，the research history，progressand future trend：A review[J].Interpretation， 2023，11（1）： SA127-SA154.

[49]Straub K M， Ganti V， Paola C， et al. Prevalence of exponential bed thickness distributions in the stratigraphic record： Experiments and theory[J]. Journal ofGeophysical Research：Earth Surface，2012，117（F2）：F02003.

[50]Kelly S.Scaling and hierarchy in braided rivers and their deposits：Examples and implications for reservoir modelling[M]//SambrookSmithGH，BestJL，BristowCS，etal.Braidedrivers： Process，deposits，ecology and management. Oxford： Blackwell Publication，2006： 75-106.

Abstract：[Objective] During the sedimentation process of sandy braided rivers，thebraided channels undergo frequent，rapid，andcontinuous shifts that leadto erosion and reworking within the preserved braided river deposits， such that deposited braided bars and braided channels are characterized bya fragmented morphology，relatively smallscale，and undefined quantitative relationships.Traditional sedimentary modelsof sandybraided riversedimentationare inadequate to effctively guide thecharacterizationof subsurface reservoirs.Toclarifythe sedimentary evolution mechanisms of sandy braided rivers and establish areliable sedimentaryarchitecture model with quantitative size relationships，this study conducted aflume experiment using constant boundary conditions to reproduce the formationand evolution of asandy braided river.[Methods] Usinga laser scanner，we obtained topographic data at regulartime intervals andaccurately reconstructed athree-dimensional sedimentary architecture model of the simulatedbraided river.Furthermore，the study analyzed sedimentary evolution mechanisms，dissected the sedimentary architecture，and constructed the quantitative size andrelationship models for architecture elements.[Results] （1）In he initial stage of sandy braided river sedimentation，sediments undergo bedload transportationdownstream， forming initial bars.Thedeflectedflowconverges into channels，furtherdeveloping lobe-shaped initialbars，which are then reshaped and organized intoa stable network of braided channelsand bar paterns through the action of the braided channels.（2）Following the formation of the braided river，the braided channelsand bars continually co-evolve，primarilythrough three mechanisms：lateral acretion of bars driven bybraided channels，abandonment and infilling of braided channels overlaying existing bars，and scouring of theconfluence bybraided channels，which results ndownstreambar reorganization.（3）During simulation，bars form within1to6run steps and grow to their maximum size beforeexperiencingcontinuous erosion over 3to8 simulationperiods，ultimately leading to preservation in only 36.28% of the area.（4） Upon completion of the simulation，the internal structure of the braided river deposit is dominated by braided channel deposits，accounting for approximately 57.9% ，which can be classified into complex stacked，incised，and isolated channels.Bars often undergo erosion and reworking from channels，presentingassmaller，fragmented forms.（5）Within thepreserveddeposit，the average width-to-thickness ratioofbraided channelsis14.1，withaninternalaccretionratioof13.7，whereasforbars，theratiois19.8，withaninternalaccretionratio of 25.4.[Conclusions] This study constructed the complex sedimentary architecture formed within sandy braided riverdeposits after persistentand intensive erosional modification bythe braided channel network，establishing aquantitative model of the size and relationships between internal architectural elements and providing a more geologically realistic and quantitative sedimentary architecture model for characterizing subsurface reservoirs.

Key words：sandy braided river;quantitative flume experiment；depositional evolution；sedimentary architecture; quantitative scale and relationship model