高级搜索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

不同流量砂质辫状河沉积演化与构型差异解析——基于沉积数值模拟实验研究

郭颖 陈奎 高青松 任广磊 李少华 冯文杰 刘林松 冉辉

郭颖, 陈奎, 高青松, 任广磊, 李少华, 冯文杰, 刘林松, 冉辉. 不同流量砂质辫状河沉积演化与构型差异解析——基于沉积数值模拟实验研究[J]. 沉积学报, 2026, 44(4): 1426-1441. doi: 10.14027/j.issn.1000-0550.2025.011
引用本文: 郭颖, 陈奎, 高青松, 任广磊, 李少华, 冯文杰, 刘林松, 冉辉. 不同流量砂质辫状河沉积演化与构型差异解析——基于沉积数值模拟实验研究[J]. 沉积学报, 2026, 44(4): 1426-1441. doi: 10.14027/j.issn.1000-0550.2025.011
GUO Ying, CHEN Kui, GAO QingSong, REN GuangLei, LI ShaoHua, FENG WenJie, LIU LinSong, RAN Hui. Sedimentary Evolution and Architecture in Sandy Braided Rivers with Different Discharge: Insights from sedimentary numerical simulation[J]. Acta Sedimentologica Sinica, 2026, 44(4): 1426-1441. doi: 10.14027/j.issn.1000-0550.2025.011
Citation: GUO Ying, CHEN Kui, GAO QingSong, REN GuangLei, LI ShaoHua, FENG WenJie, LIU LinSong, RAN Hui. Sedimentary Evolution and Architecture in Sandy Braided Rivers with Different Discharge: Insights from sedimentary numerical simulation[J]. Acta Sedimentologica Sinica, 2026, 44(4): 1426-1441. doi: 10.14027/j.issn.1000-0550.2025.011

不同流量砂质辫状河沉积演化与构型差异解析——基于沉积数值模拟实验研究

doi: 10.14027/j.issn.1000-0550.2025.011
基金项目: 

中国石化科技攻关项目 P23031

详细信息

Sedimentary Evolution and Architecture in Sandy Braided Rivers with Different Discharge: Insights from sedimentary numerical simulation

More Information
  • 摘要: 目的 辫状河的沉积演化、形态及规模,受流量这一关键因素的影响,但目前尚未系统分析其对辫状水道、心滩及内部结构的影响程度。 方法 采用Delft3D软件通过5组不同流量三维沉积数值模拟实验,探究不同流量下辫状水道、心滩的演化、形态样式及规模差异,再现不同流量辫状河多沉积类型砂坝内部叠置样式和构型特征,总结砂质辫状河沉积特征和演化规律。 结果 (1)流量与辫状河的沉积演化速率成正比。低流量辫状河演化程度低,砂坝以向下迁移为主,不易被河道分割;而高流量辫状河演化程度高,砂坝以侧向迁移为主,心滩合并、被分割的动态平衡过程较普遍。(2)流量是决定辫状水道平面样式和水深的主导因素。流量越大,河道越宽、直、深,分岔少,演化成熟阶段辫状水道网络复杂;而流量较小的辫状河河道,河道水浅、窄、分岔多,演化成熟阶段辫状水道以相对固定的交织河网为主。(3)流量进一步影响辫状水道和心滩的协同演化过程。流量越大,心滩样式更复杂,出现规模较大的复合改造坝和河边滩,其内部叠切关系更复杂,更易发育窜沟。(4)虽然砂质辫状河心滩内部结构具有很强的空间变异性,但仍存在相似点。相同点在于自下而上均发育顺流加积体、侧积体、垂向加积体,不同点在于流量越大,心滩内部侧积体占比增加。(5)流量与增生体规模呈正相关关系,但与长宽比、长厚比呈负相关关系,其中低流量辫状河心滩长宽比中位数为3.5,高流量辫状河心滩长宽比为4.5。 结论 该研究提出了流量与砂质辫状河的心滩及内部增生体的规模成正比的观点,明确了砂质辫状河受流量影响沉积特征和规模差异,为地下复杂储层构型解剖提供依据。
  • 图  1  砂质辫状河模拟网格与初始水深分布

    Figure  1.  Grid and initial water depths for numerical simulation of sandy braided rivers

    图  2  高、低流量砂质辫状河演化过程对比

    (a-c) low-discharge rivers at steps 1 510, 3 250 and 3 890; (d-f) high-discharge rivers at steps 600, 810 and 1 100

    Figure  2.  Comparision of sedimentary evolution in high⁃and low⁃discharge sandy braided rivers

    Fig.2

    图  3  高、低流量砂质辫状河心滩演化模式对比(Guo et al.,2023

    (a-c) low discharge; (d-f) high discharge

    Figure  3.  Comparision of Channel bar evolution model in high⁃and low⁃discharge sandy braided river evolution (Guo et al., 2023)

    Fig.3

    图  4  不同流量下砂质辫状河演化过程中砂坝规模演变趋势

    (a) average length; (b) average width; (c) average aspect ratio; (d) average area

    Figure  4.  Trend in sandbar size evolution during the evolution of sandy braided rivers under different discharge conditions

    Fig.4

    图  5  不同流量条件对照模拟组砂质辫状河剖面

    Figure  5.  Simulated profiles of sandy braided rivers at different discharge conditions

    图  6  不同流量成熟时期辫状水道分析图

    (a) sedimentary water depth; (b) number, curvature and bifurcation coefficient of braided channels

    Figure  6.  Analysis diagrams of braided channels at different discharge mature stages

    Fig.6

    图  7  不同流量下模拟砂质辫状河演化成熟时期沉积地貌图

    Figure  7.  Depositional geomorphology at mature stages of simulated sandy braided rivers under different discharge

    图  8  不同流量成熟时期砂坝参数

    (a) distribution of area proportion for different types of sandbars; (b) distribution of aspect ratios for different types of sandbars

    Figure  8.  Sandbar parameters at mature stages for different discharge

    Fig.8

    图  9  不同流量下模拟砂质辫状河演化成熟阶段沉积微相分布图

    Figure  9.  Distribution of sedimentary microfacies at mature stages for simulated sandy braided rivers in different flows

    图  10  低流量辫状河演化成熟阶段单元坝沉积构型

    (a) transverse section; (b) longitudinal section

    Figure  10.  Sedimentary architecture of unit bar developed with low discharge braided river at mature stages

    Fig.10

    图  11  高流量辫状河演化成熟阶段单元坝沉积构型

    (a) longitudinal section AA΄; (b) longitudinal section CC΄; (c) plan view of the sandbar; (d) water flow mode

    Figure  11.  Sedimentary architecture of unit bar developed with high discharge braided river at mature stages

    Fig.11

    图  12  低流量辫状河演化成熟阶段复合坝沉积构型

    (a) transverse section AA΄; (b) transverse section BB΄; (c) longitudinal section CC΄; (d) plan view of the sandbar;①-⑤ represents different stages of accretions

    Figure  12.  Sedimentary architecture of compound bar developed with low discharge braided river at mature stages

    Fig.12

    图  13  高流量辫状河演化成熟阶段复合坝沉积构型

    (a) transverse section AA΄; (b) transverse section BB΄; (c) transverse section CC΄; (d) longitudinal section DD΄; (e) plan view of the sandbar;①-⑩ represents different stages of accretions

    Figure  13.  Sedimentary architecture of compound bar developed with high discharge braided river at mature stages

    Fig.13

    图  14  高、低流量辫状河单元坝、复合坝沉积构型规模参数对比

    (a) dip angle; (b) length; (c) width; (d) thickness

    Figure  14.  Comparison of sedimentary architecture scale parameters in unit bar and compound bar of high⁃and low⁃discharge braided river

    Fig.14

    图  15  高流量辫状河单元坝与野外露头沉积构型对比

    (a) plan view of sandbar; (b) transverse section CC’; (c1, c2) field outcrop section and dissection diagram

    Figure  15.  Comparison of sedimentary architecture between highdischarge braided river sandbar and rock outcrop

    Fig.15

    图  16  高流量辫状河单元坝与井下实钻案例沉积构型对比

    (a) transverse section; (b) longitudinal source profile; (c) underground wells at different locations

    Figure  16.  Comparison of sedimentary architecture between high⁃discharge braided river sandbar and drill core data

    Fig.16

    表  1  砂质辫状河流量等参数调研(据刘警阳等,2024修改)

    Table  1.   Hydrological parameter survey of sandy braided river (modified from Liu et al., 2024)

    序号河流位置水文站平均洪峰流量m3/s河道宽度/m宽度测量坐标等效流量m2/s
    1亚马孙河巴西Santo Antonio Do Ica68 192.9910 53259°42'48" W,3°2'52" S6.47
    2巴拉那河阿根廷Corrientes21 917.691 72263°13'51" W,23°47'27" S12.73
    3马更些河加拿大Norman Wells18 130.552 811111°32'40" W,59°26'18" N6.45
    4贾木纳河印度Bahadurabad50 531.083 14084°38'41" E,25°19'54" N16.09
    5恒河印度Farakka44 419.305 13680°29'20" E,18°0'28" N8.65
    6Rakaia river新西兰Fighting Hill398.31606168°18'23" E,45°55'10" N0.66
    7育空河美国Ruby, Alas12 966.214 914131°5'29" W,67°3'4" N2.64
    8维斯瓦河波兰Warsaw1 151.131 35020°21'18" E,52°22'20" N0.85
    9奈厄布拉勒河美国Near Verdel, Nebr98.111 787107°34'12" W,46°6'58" N0.05
    10南萨斯喀彻温河加拿大Saskatoon846.79550111°32'52" W,56°50'2" N1.54
    11松花江中国哈尔滨3 397.063 230129°2'26" E,45°46'11" N1.05
    12普拉特河美国Louisville, Nebr534.692 26677°14'46" W,42°40'20" N0.24
    13Halil Rud river伊朗Hossein Abad Jiroft36.4040757°44'56" E,28°24'20" N0.09
    下载: 导出CSV

    表  2  砂质辫状河沉积数值模拟参数

    Table  2.   Sedimentary numerical simulation parameters of sandy braided river

    参数项模拟设定值说明
    沉积物粒度中值/μm300
    泥质沉积物供应量/(kg/m3)10参考现代沉积数据
    砂质沉积物浓度/(kg/m3)2.5
    沉积物组分细砂中砂粗砂巨砂多种组分混合后中值粒径约为300 μm
    砂质沉积物组分粒度中值/μm1503007501 500
    砂质沉积物组分含量/(质量占比,%)25501510
    沉积物浓度/(kg/m3)1.252.50.750.5
    砂质沉积物密度/(g/cm3)2.65
    底床坡度/(°)9.3×10-5冲积平原坡度平缓
    河道宽度/km3.2
    水流量/(m3/s)5 00010 00020 00040 00080 000变量
    等效流量1.563.126.2512.525.0m2/s
    参考河流南萨斯喀彻温河1.54育空河2.64亚马孙河6.47贾木纳河16.09
    演化成熟步长3 8902 9901 5001 100880
    平均水深/m4~6取决于河宽、底床坡度及水
    控制流速/(m/s)0~4.5
    网格单元大小/m50×20
    平面网格单元总量/个1 600×160
    模拟区域尺寸/m80 000×3 200
    地貌演化系数48沉积模拟加速48倍
    底床糙度-Chezy值/(m1/2/s)45经验参数
    水平涡黏系数/(m2/s)*1 000
    垂向涡黏系数/(m2/s)*0.001
    临近干网格侵蚀系数0.25
    水动力计算时间步长/min0.2
    沉积模拟数据保存时间步长/min15
    注:“—”表示无补充说明。
    下载: 导出CSV
  • [1] 陈东阳,王峰,陈洪德,等. 2019. 鄂尔多斯盆地东部府谷天生桥剖面上古生界下石盒子组8段辫状河储层构型表征[J]. 石油与天然气地质,40(2):335-345.

    Chen Dongyang, Wang Feng, Chen Hongde, et al. 2019. Characterization of braided river reservoir architecture of the Upper Paleozoic He 8 member on Fugu Tian-shengqiao outcrop, eastern Ordos Basin[J]. Oil & Gas Geology, 40(2): 335-345.
    [2] 冯文杰,曹荆楚,张昌民,等. 2025. 砂质辫状河沉积演化机制与沉积构型模式:定量化水槽沉积模拟实验研究[J]. 沉积学报,43(3):782-796.

    Feng Wenjie, Cao Jingchu, Zhang Changmin, et al. 2025. Sedimentary evolution mechanisms and architecture models of sandy braided rivers: A study based on quantitative flume experiments[J]. Acta Sedimentologica Sinica, 43(3): 782-796.
    [3] 郭颖,冯文杰,李少华,等. 2024. 砂质辫状河中河边滩和河心滩沉积构型:基于沉积演化过程的精细解析[J]. 地球科学进展,39(4):405-418.

    Guo Ying, Feng Wenjie, Li Shaohua, et al. 2024. Sedimentary architecture between compound middle bars and side bars: Insights from a numerical model of sandy braided rivers[J]. Advances in Earth Science, 39(4): 405-418.
    [4] 何维领,李少华,王濡岳,等. 2024. 砂质辫状河储层构型特征及沉积演化:以山西大同侏罗系云冈组露头为例[J]. 沉积学报,42(5):1699-1710.

    He Weiling, Li Shaohua, Wang Ruyue, et al. 2024. Sandy braided river architecture characteristics and evolution: A case study from outcrops in the Middle Jurassic Yungang Formation, Datong, Shanxi province[J]. Acta Sedimentologica Sinica, 42(5): 1699-1710.
    [5] 贾爱林,王国亭,孟德伟,等. 2018. 大型低渗-致密气田井网加密提高采收率对策:以鄂尔多斯盆地苏里格气田为例[J]. 石油学报,39(7):802-813.

    Jia Ailin, Wang Guoting, Meng Dewei, et al. 2018. Well pattern infilling strategy to enhance oil recovery of giant low-permeability tight gasfield: A case study of Sulige gasfield, Ordos Basin[J]. Acta Petrolei Sinica, 39(7): 802-813.
    [6] 雷涛,任广磊,李晓慧,等. 2023. 砂质辫状河心滩沉积演化规律与沉积构型特征:基于沉积数值模拟的认识[J]. 石油与天然气地质,44(6):1595-1608.

    Lei Tao, Ren Guanglei, Li Xiaohui, et al. 2023. Sedimentary evolution pattern and architectural characteristics of mid-channel bars in sandy braided rivers: Understanding based on sedimentary numerical simulation[J]. Oil & Natural Gas Geology, 44(6): 1595-1608.
    [7] 李伟,岳大力,王武荣,等. 2023. 辫状河沉积构型研究进展:沉积演化与构型特征[J]. 古地理学报,25(5):1032-1048.

    Li Wei, Yue Dali, Wang Wurong, et al. 2023. Depositional models of braided rivers: Characteristics of sedimentary evolution and architecture[J]. Journal of Palaeogeography, 25(5): 1032-1048.
    [8] 廖保方,张为民,李列,等. 1998. 辫状河现代沉积研究与相模式:中国永定河剖析[J]. 沉积学报,16(1):34-39,50.

    Liao Baofang, Zhang Weimin, Li Lie, et al. 1998. Study on modern deposit of a braided stream and factes model: Yaking the Yongding River as an example[J]. Acta Sedimentologica Sinica, 16(1): 34-39, 50.
    [9] 刘警阳,李伟,赵文智,等. 2024. 流量变化控制的辫状河心滩演化过程[J]. 地球科学,49(10):3516-3528.

    Liu Jingyang, Li Wei, Zhao Wenzhi, et al. 2024. Evolution of bars in braided rivers controlled by discharge variability[J]. Earth Science, 49(10): 3516-3528.
    [10] 牛博,高兴军,赵应成,等. 2015. 古辫状河心滩坝内部构型表征与建模:以大庆油田萨中密井网区为例[J]. 石油学报,36(1):89-100.

    Niu Bo, Gao Xingjun, Zhao Yingcheng, et al. 2015. Architecture characterization and modeling of channel bar in paleo-braided river: A case study of dense well pattern area of Sazhong in Daqing oilfield[J]. Acta Petrolei Sinica, 36(1): 89-100.
    [11] 任晓旭. 2020. 砂质辫状河构型界面附近岩性和物性分布模式研究[D]. 北京:中国石油大学(北京).

    Ren Xiaoxu. 2020. Lithology and physical property distribution patterns of sandy braided river near architectural bounding surfaces[D]. Beijing: China University of Petroleum (Beijing).
    [12] 唐洪,龙冠宇,张章,等. 2026. 基于沉积数值模拟的砂质辫状河三角洲沉积特征与演化规律研究[J]. 沉积学报,44(3):841-855.

    Tang Hong, Long Guanyu, Zhang Zhang, et al. 2026. Sedimentary characteristics and evolution laws of a sandy braided river delta based on sedimentary numerical simulation[J]. Acta Sedimentologica Sinica, 44(3): 841-855.
    [13] 于兴河. 2004. 辫状河储层地质模式及层次界面分析[M]. 北京:石油工业出版社.

    Yu Xinhe. 2004. Braided river reservoir geological model and hierarchical interface analysis [M]. Beijing: Petroleum Industry Press.
    [14] 张昌民,尹太举,赵磊,等. 2013. 辫状河储层内部建筑结构分析[J]. 地质科技情报,32(4):7-13.

    Zhang Changmin, Yin Taiju, Zhao Lei, et al. 2013. Reservoir architectural analysis of braided channel[J]. Geological Science and Technology Information, 32(4): 7-13.
    [15] 张可,吴胜和,冯文杰,等. 2018. 砂质辫状河心滩坝的发育演化过程探讨:沉积数值模拟与现代沉积分析启示[J]. 沉积学报,36(1):81-91.

    Zhang Ke, Wu Shenghe, Feng Wenjie, et al. 2018. Discussion on evolution of bar in sandy braided river: Insights from sediment numerical simulation and modern bar[J]. Acta Sedimentologica Sinica, 36(1): 81-91.
    [16] 张宪国,张育衡,张涛,等. 2020. 基于沉积数值模拟的辫状河心滩演化[J]. 中国石油大学学报(自然科学版),44(2):1-9.

    Zhang Xianguo, Zhang Yuheng, Zhang Tao, et al. 2020. Analysis of braided bar evolution based on numerical simulation of deposition process[J]. Journal of China University of Petroleum (Edition of Natural Science), 44(2): 1-9.
    [17] Ashworth P J, Best J L, Roden J E, et al. 2000. Morphological evolution and dynamics of a large, sand braid‐bar, Jamuna River, Bangladesh[J]. Sedimentology, 47(3): 533-555.
    [18] Ashworth P J, Sambrook Smith G H, Best J L, et al. 2011. Evolution and sedimentology of a channel fill in the sandy braided South Saskatchewan River and its comparison to the deposits of an adjacent compound bar[J]. Sedimentology, 58(7): 1860-1883.
    [19] Best J L, Ashworth P J, Bristow C S, et al. 2003. Three-dimensional sedimentary architecture of a large, mid-channel sand braid bar, Jamuna River, Bangladesh[J]. Journal of Sedimentary Research, 73(4): 516-530.
    [20] Bridge J S, Lunt I A. 2006. Depositional models of braided rivers[M]//Sambrook Smith G H, Best J L, Bristow C S, et al.Braided rivers: Process, deposits, ecology and management. Oxford, UK: International Association of Sedimentologists Special Publication, Blackwell: 11-50.
    [21] Cant D J, Walker R G. 1978. Fluvial processes and facies sequences in the sandy braided South Saskatchewan River, Canada[J]. Sedimentology, 25(5): 625-648.
    [22] Ethridge F G. 2011. Interpretation of ancient fluvial channel deposits: Review and recommendations[M]//Davidson S K, Leleu S, North C P. From river to rock record: The preservation of fluvial sediments and their subsequent interpretation. Society of Economic Paleontologists and Mineralogists Special Publication, 97: 9-35.
    [23] Feng W J, Fei-Xiang G, Zhang C M, et al. 2025. Sedimentary architecture of a sandy braided river: Insights from a flume experiment[J]. Petroleum Science, 22(1): 16-28.
    [24] Fielding C R, Alexander J, Allen J P. 2018. The role of discharge variability in the formation and preservation of alluvial sediment bodies[J]. Sedimentary Geology, 365: 1-20.
    [25] Fielding C R, Alexander J, Argiro A. 2025. Towards a new generation of fluvial facies models for the interpretation of ancient deposits, based on inter‐annual peak discharge variance of modern rivers[J]. Sedimentology, 72(1): 275-329.
    [26] Guo Y, Li S H, He F Q, et al. 2023. Sedimentary characteristics and internal architecture of a sandy braided river: Insights from a numerical model[J]. Water, 15(23): 4119.
    [27] Li W, Colombera L, Yue D L, et al. 2023. Controls on the morphology of braided rivers and braid bars: An empirical characterization of numerical models[J]. Sedimentology, 70(1): 259-279.
    [28] Miall A D. 2022. Facies models[M]//Miall A D. Stratigraphy: A modern synthesis. Cham: Springer: 175-230.
    [29] Nicholas A P. 2013. Modelling the continuum of river channel patterns[J]. Earth Surface Processes and Landforms, 38(10): 1187-1196.
    [30] Puig J M, Cabello P, Howell J, et al. 2019. Three-dimensional characterisation of sedimentary heterogeneity and its impact on subsurface flow behaviour through the braided-to-meandering fluvial deposits of the Castissent Formation (Late Ypresian, Tremp-Graus Basin, Spain)[J]. Marine and Petroleum Geology, 103: 661-680.
    [31] Rice S P, Church M, Wooldridge C L, et al. 2009. Morphology and evolution of bars in a wandering gravel-bed river; Lower Fraser river, British Columbia, Canada[J]. Sedimentology, 56(3): 709-736.
    [32] Rust B R. 1977. A classification of alluvial channel systems[M]//Miall A D. Fluvial Sedimentology-Memoir of the Canadian Society of Petroleum Geologists, 5: 187-198.
    [33] Sambrook Smith G H, Ashworth P J, Best J L, et al. 2006. The sedimentology and alluvial architecture of the sandy braided South Saskatchewan River, Canada[J]. Sedimentology, 53(2): 413-434.
    [34] Sambrook Smith G H, Ashworth P J, Best J L, et al. 2009. The sedimentology and alluvial architecture of a large braid bar, Río Paraná, Argentina[J]. Journal of Sedimentary Research, 79(8): 629-642.
    [35] Schuurman F, Kleinhans M G. 2015. Bar dynamics and bifurcation evolution in a modelled braided sand‐bed river[J]. Earth Surface Processes and Landforms, 40(10): 1318-1333.
    [36] Schuurman F, Kleinhans M G, Middelkoop H. 2016. Network response to disturbances in large sand-bed braided rivers[J]. Earth Surface Dynamics, 4(1): 25-45.
    [37] Schuurman F, Marra W A, Kleinhans M G. 2013. Physics‐based modeling of large braided sand‐bed rivers: Bar pattern formation, dynamics, and sensitivity[J]. Journal of Geophysical Research: Earth Surface, 118(4): 2509-2527.
    [38] Straub K M, Ganti V, Paola C, et al. 2012. Prevalence of exponential bed thickness distributions in the stratigraphic record: Experiments and theory[J]. Journal of Geophysical Research: Earth Surface, 117(F2): F02003.
    [39] Tamminga A D, Eaton B C, Hugenholtz C H. 2015. UAS-based remote sensing of fluvial change following an extreme flood event[J]. Earth Surface Processes and Landforms, 40(11): 1464-1476.
    [40] van de Lageweg W I, van Dijk W M, Kleinhans M G. 2013. Morphological and stratigraphical signature of floods in a braided gravel-bed river revealed from flume experiments[J]. Journal of Sedimentary Research, 83(11): 1033-1046.
  • [1] 唐洪, 黄帅, 龙冠宇, 桂万铃, 卫营, 刘博轩, 陈家驹.  关键边界条件对陆相辫状河三角洲沉积过程的控制作用——基于沉积数值模拟的探讨 . 沉积学报, 2026, (): -. doi: 10.14027/j.issn.1000-0550.2026.015
    [2] 唐洪, 龙冠宇, 张章, 曲炳昌, 蒋成思.  基于沉积数值模拟的砂质辫状河三角洲沉积特征与演化规律研究 . 沉积学报, 2026, 44(3): 841-855. doi: 10.14027/j.issn.1000-0550.2024.088
    [3] 马福康, 徐振华, 吴胜和, 陶叶, 陈亚坤, 黎飞凡, 王武荣, 李伟.  季节性湖平面变化下浅水三角洲指状砂坝沉积特征 . 沉积学报, 2026, 44(2): 646-663. doi: 10.14027/j.issn.1000-0550.2024.091
    [4] 赵俊威, 王恒, 张东伟, 孙海航, 陈恭洋, 印森林.  海相砂质滩坝沉积构型及表征研究进展 . 沉积学报, 2025, 43(4): 1275-1292. doi: 10.14027/j.issn.1000-0550.2023.108
    [5] 彭晨阳, 唐明明, 洪瑞峰, 熊思琛.  基于沉积动力学的潮控三角洲沉积数值模拟研究 . 沉积学报, 2025, 43(3): 846-859. doi: 10.14027/j.issn.1000-0550.2024.081
    [6] 冯文杰, 曹荆楚, 张昌民, 钱其豪, 张涛, 雷涛, 周志成.  砂质辫状河沉积演化机制与沉积构型模式——定量化水槽沉积模拟实验研究 . 沉积学报, 2025, 43(3): 782-796. doi: 10.14027/j.issn.1000-0550.2024.068
    [7] 陈彬滔, 马轮, 洪亮, 史忠生, 庞文珠, 薛罗, 王磊, 代寒松, 赵艳军.  断陷湖盆陡坡带扇三角洲—滑塌扇复合扇体的沉积演化及油气地质意义 . 沉积学报, 2024, 42(3): 1058-1072. doi: 10.14027/j.issn.1000-0550.2022.087
    [8] 赵晓明, 刘飞, 葛家旺, 冯潇飞, 张喜, 张文彪, 杨宝泉, 杨莉.  深水水道沉积构型单元分级与结构样式 . 沉积学报, 2023, 41(1): 37-51. doi: 10.14027/j.issn.1000-0550.2022.048
    [9] 王科, 赵俊峰, 薛锐, 闫占冬, 李旋, 李一凡.  鄂尔多斯盆地延安组河流沉积类型及演变 . 沉积学报, 2022, 40(5): 1367-1377. doi: 10.14027/j.issn.1000-0550.2021.026
    [10] 尹力, 冯文杰, 尹艳树, 雷诚, 徐庆岩, 何一鸣.  波浪作用下砂质滩坝的沉积过程与沉积模式 . 沉积学报, 2022, 40(5): 1393-1405. doi: 10.14027/j.issn.1000-0550.2021.147
    [11] 陈彬滔, 于兴河, 王磊, 史忠生, 马轮, 薛罗, 史江龙, 白洁, 赵艳军.  河流相沉积的河型转换特征与控制因素及其油气地质意义 . 沉积学报, 2021, 39(2): 424-433. doi: 10.14027/j.issn.1000-0550.2020.036
    [12] 吕峻岭, 朱一杰, 夏瑞, 郑云柯, 刘晨虎, 冯文杰, 李国艳, 杜晓峰.  干旱型分支河流体系沉积特征与演化过程 . 沉积学报, 2020, 38(5): 994-1005. doi: 10.14027/j.issn.1000-0550.2019.101
    [13] 张可, 吴胜和, 冯文杰, 郑定业, 喻宸, 刘照玮.  砂质辫状河心滩坝的发育演化过程探讨——沉积数值模拟与现代沉积分析启示 . 沉积学报, 2018, 36(1): 81-91. doi: 10.3969/j.issn.1000-0550.2018.010
    [14] 单敬福, 张彬, 赵忠军, 李浮萍, 王辉, 王博.  复合辫状河道期次划分方法与沉积演化过程分析——以鄂尔多斯盆地苏里格气田西区苏X区块为例 . 沉积学报, 2015, 33(4): 773-785. doi: 10.14027/j.cnki.cjxb.2015.04.016
    [15] 赵俊峰, 屈红军, 林晋炎, 刘祥, 杨扬, 林川.  湖泊三角洲沉积露头精细解剖——以鄂尔多斯盆地裴庄剖面为例 . 沉积学报, 2014, 32(6): 1026-1034.
    [16] 塔里木盆地西北缘乌恰地区晚白垩世—渐新世微相分析及沉积演化 . 沉积学报, 2012, 30(3): 501-510.
    [17] 王随继.  网状河流的构型、流量—宽深比关系和能耗率 . 沉积学报, 2003, 21(4): 565-570.
    [18] 王随继.  大同盆地中侏罗世河流沉积体系及古河型演化 . 沉积学报, 2001, 19(4): 501-505.
    [19] 伍涛, 王建国, 王德发.  辫状河砂体储层沉积学研究──以张家口地区露头砂体为例 . 沉积学报, 1998, 16(1): 27-33.
    [20] 柯保嘉, 陈志明, 黄家宽, 陈昌明, 汪寿松, 陈安宁.  山西二叠系河流沉积特征 . 沉积学报, 1991, 9(1): 11-19.
  • 加载中
图(16) / 表 (2)
计量
  • 文章访问数:  423
  • HTML全文浏览量:  61
  • PDF下载量:  383
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-11-05
  • 修回日期:  2025-03-16
  • 录用日期:  2025-05-12
  • 网络出版日期:  2025-05-12
  • 刊出日期:  2026-08-10

目录

    不同流量砂质辫状河沉积演化与构型差异解析——基于沉积数值模拟实验研究

    doi: 10.14027/j.issn.1000-0550.2025.011
      基金项目:

      中国石化科技攻关项目 P23031

      作者简介:

      郭颖,女,1992年出生,博士研究生,储层沉积学,E-mail: guoying.hbsj@sinopec.com

      通讯作者: 李少华,男,教授,E-mail: lish@yangtzeu.edu.cn
    • 中图分类号: P512.2

    摘要: 目的 辫状河的沉积演化、形态及规模,受流量这一关键因素的影响,但目前尚未系统分析其对辫状水道、心滩及内部结构的影响程度。 方法 采用Delft3D软件通过5组不同流量三维沉积数值模拟实验,探究不同流量下辫状水道、心滩的演化、形态样式及规模差异,再现不同流量辫状河多沉积类型砂坝内部叠置样式和构型特征,总结砂质辫状河沉积特征和演化规律。 结果 (1)流量与辫状河的沉积演化速率成正比。低流量辫状河演化程度低,砂坝以向下迁移为主,不易被河道分割;而高流量辫状河演化程度高,砂坝以侧向迁移为主,心滩合并、被分割的动态平衡过程较普遍。(2)流量是决定辫状水道平面样式和水深的主导因素。流量越大,河道越宽、直、深,分岔少,演化成熟阶段辫状水道网络复杂;而流量较小的辫状河河道,河道水浅、窄、分岔多,演化成熟阶段辫状水道以相对固定的交织河网为主。(3)流量进一步影响辫状水道和心滩的协同演化过程。流量越大,心滩样式更复杂,出现规模较大的复合改造坝和河边滩,其内部叠切关系更复杂,更易发育窜沟。(4)虽然砂质辫状河心滩内部结构具有很强的空间变异性,但仍存在相似点。相同点在于自下而上均发育顺流加积体、侧积体、垂向加积体,不同点在于流量越大,心滩内部侧积体占比增加。(5)流量与增生体规模呈正相关关系,但与长宽比、长厚比呈负相关关系,其中低流量辫状河心滩长宽比中位数为3.5,高流量辫状河心滩长宽比为4.5。 结论 该研究提出了流量与砂质辫状河的心滩及内部增生体的规模成正比的观点,明确了砂质辫状河受流量影响沉积特征和规模差异,为地下复杂储层构型解剖提供依据。

    English Abstract

    郭颖, 陈奎, 高青松, 任广磊, 李少华, 冯文杰, 刘林松, 冉辉. 不同流量砂质辫状河沉积演化与构型差异解析——基于沉积数值模拟实验研究[J]. 沉积学报, 2026, 44(4): 1426-1441. doi: 10.14027/j.issn.1000-0550.2025.011
    引用本文: 郭颖, 陈奎, 高青松, 任广磊, 李少华, 冯文杰, 刘林松, 冉辉. 不同流量砂质辫状河沉积演化与构型差异解析——基于沉积数值模拟实验研究[J]. 沉积学报, 2026, 44(4): 1426-1441. doi: 10.14027/j.issn.1000-0550.2025.011
    GUO Ying, CHEN Kui, GAO QingSong, REN GuangLei, LI ShaoHua, FENG WenJie, LIU LinSong, RAN Hui. Sedimentary Evolution and Architecture in Sandy Braided Rivers with Different Discharge: Insights from sedimentary numerical simulation[J]. Acta Sedimentologica Sinica, 2026, 44(4): 1426-1441. doi: 10.14027/j.issn.1000-0550.2025.011
    Citation: GUO Ying, CHEN Kui, GAO QingSong, REN GuangLei, LI ShaoHua, FENG WenJie, LIU LinSong, RAN Hui. Sedimentary Evolution and Architecture in Sandy Braided Rivers with Different Discharge: Insights from sedimentary numerical simulation[J]. Acta Sedimentologica Sinica, 2026, 44(4): 1426-1441. doi: 10.14027/j.issn.1000-0550.2025.011
      • 砂质辫状河为低弯度、河床坡降大、河道水浅流急易分岔合并、频繁迁移改道的常见河流地貌(Rust,1977Ashworth et al.,20002011于兴河,2004),储层以粗粒岩性为主体、多低渗缓冲带的强非均质性、整体泛连通但内部呈结构复杂的迷宫状结构(张昌民等,2013张可等,2018)。辫状水道与心滩是砂质辫状河的主要沉积单元(Cant and Walker,1978Bridge and Lunt,2006),受异源因素(allogenic)的影响(Mail,2022),如植被、气候控制的流量变化等,辫状河形态和保存的沉积序列存在明显差异,其中流量对河流沉积的影响是至关重要的(Schuurman et al.,2013Tamminga et al.,2015)。即使古气候、古沉积背景相似,仍存在形态规模、内部构型样式差异很大的辫状河,Ethridge(2011)认为是流量差异造成上游的供水条件不同,导致同一个地区出现多条不同类型的辫状河共存的现象。如鄂尔多斯盆地下石盒子组辫状河,其坡度、沉积物成分、沉积物粒度等基本沉积背景相近,但在流量等方面,不同时期、发育于不同位置的辫状河存在较大的差异。在既往对下石盒子组辫状河储层的解剖研究中,前人已发现流量不同的辫状河,河道的宽度、砂体的规模存在一定程度的差异(贾爱林等,2018陈东阳等,2019)。

        许多学者意识到流量对辫状河是至关重要的,但仅定性地提出流量影响心滩的沉积演化和规模,没有系统地对辫状水道、心滩及内部结构定量评价。其中Fielding et al.(20182025)认识到洪峰流量与保存下来的沉积物存在极强的相关性,可作为评价河流流量变化的度量,对沉积物的搬运、沉积构造的发育以及沉积特征的保存具有重要意义。刘警阳等(2024)基于现代沉积对比提出流量的差异会控制心滩的沉积演化,进一步影响其形态,认为低流量辫状河的心滩和辫状水道之间交互作用较弱,形态可以维持相对稳定;而高流量辫状水道对心滩形态的改造作用非常强烈,导致心滩发育位置变化明显,辫状水道常发生迁移改道,多为游荡型辫状河。李伟等(2023)基于多组沉积模拟提出流量对心滩的规模有控制作用。与此同时,大多数现代河流测量时间不够长,任一洪水事件形成的沉积结构会被后续的洪水改造使其不具有代表性。Nicholas(2013)开展了物理模拟、探地雷达与取心对比,使在较大的时间和空间尺度上模拟河道形态动力具有可行性。

        因此,本研究采用Delft3D软件通过多组不同流量三维沉积数值模拟实验进行时空迁移分析,探究不同流量下辫状水道、心滩及内部增生体的类型、规模等差异,明确了流量对砂质辫状河的叠置样式、构型特征与内部非均质性的影响,对高效开发地下油气砂体分布预测提供模式指导。

      • 水动力条件是影响河流沉积演化的关键边界条件,在沉积底床、河床形态、沉积物类型与粒度等边界条件相同的情况下,水动力强度的变化可引起辫状河水深、流速等水动力参数的变化,进而对辫状河沉积模拟产生决定性影响。为了明确水动力强度对砂质辫状河沉积的控制作用,通过Delft3D软件设计不同流量的对照方案开展沉积数值模拟,以探索不同流量条件下辫状河沉积特征及其差异。

      • 模拟参数设置主要包含网格参数、粒度选择和水动力参数三个方面,其中网格参数设计是基于研究区同盆地苏里格气田及府谷天生桥露头同层位辫状河单砂体解剖成果(陈东阳等,2019),结合单期增生体的最小几何尺度约束,将模拟网格尺寸设定为50 m×20 m。该设计确保每期沉积单元可覆盖4~5个网格,满足增生体空间展布的精度需求。粒度数据选取整合了地表露头与地下岩心的对比分析,其中府谷天生桥露头18组样品测试显示,岩心粒度中值为300 µm;鄂尔多斯盆地北缘石盒子组砂质辫状河293个地下样品分析表明,峰值粒径为320 µm,与露头数据趋势吻合。据此采用露头数据为基准,按细砂(25%)、中砂(50%)、粗砂(15%)、巨砂(10%)的组分比例调配沉积物,使中值粒径趋近300 µm。水动力参数上,通过筛选全球典型现代砂质辫状河案例(粒度介于0.2~0.4 mm、缓坡区坡度介于0.05°~0.15°、河道宽度介于1~10 km),结合中国水利部砂泥比参数(0.02),最终确定模拟输沙量为2.5 kg/m³,实现沉积动力学过程定量约束。

        模拟区为宽3.2 km,长80 km的矩形,模拟区左侧边界为进水口,右侧边界为出水口(图1),模拟过程中总流量保持恒定。流量是单位时间内流经封闭管道有效截面的流体量,现代河流中根据流量均衡原理,流量增加往往辫状河河道更宽。但在模拟中采用槽型模拟区固定宽度3.2 km,调研现代沉积13条河流的流量数据,计算单位宽度流量为等效流量,数据范围从0.09~16.9 m2/s(表1)。以此参考设计总流量在5 000~80 000 m3/s为变量对照组试验模拟了5条辫状河,其主要沉积背景条件相同,其他参数不变开展五组沉积模拟(表2)。

        图  1  砂质辫状河模拟网格与初始水深分布

        Figure 1.  Grid and initial water depths for numerical simulation of sandy braided rivers

        表 1  砂质辫状河流量等参数调研(据刘警阳等,2024修改)

        Table 1.  Hydrological parameter survey of sandy braided river (modified from Liu et al., 2024)

        序号河流位置水文站平均洪峰流量m3/s河道宽度/m宽度测量坐标等效流量m2/s
        1亚马孙河巴西Santo Antonio Do Ica68 192.9910 53259°42'48" W,3°2'52" S6.47
        2巴拉那河阿根廷Corrientes21 917.691 72263°13'51" W,23°47'27" S12.73
        3马更些河加拿大Norman Wells18 130.552 811111°32'40" W,59°26'18" N6.45
        4贾木纳河印度Bahadurabad50 531.083 14084°38'41" E,25°19'54" N16.09
        5恒河印度Farakka44 419.305 13680°29'20" E,18°0'28" N8.65
        6Rakaia river新西兰Fighting Hill398.31606168°18'23" E,45°55'10" N0.66
        7育空河美国Ruby, Alas12 966.214 914131°5'29" W,67°3'4" N2.64
        8维斯瓦河波兰Warsaw1 151.131 35020°21'18" E,52°22'20" N0.85
        9奈厄布拉勒河美国Near Verdel, Nebr98.111 787107°34'12" W,46°6'58" N0.05
        10南萨斯喀彻温河加拿大Saskatoon846.79550111°32'52" W,56°50'2" N1.54
        11松花江中国哈尔滨3 397.063 230129°2'26" E,45°46'11" N1.05
        12普拉特河美国Louisville, Nebr534.692 26677°14'46" W,42°40'20" N0.24
        13Halil Rud river伊朗Hossein Abad Jiroft36.4040757°44'56" E,28°24'20" N0.09

        表 2  砂质辫状河沉积数值模拟参数

        Table 2.  Sedimentary numerical simulation parameters of sandy braided river

        参数项模拟设定值说明
        沉积物粒度中值/μm300
        泥质沉积物供应量/(kg/m3)10参考现代沉积数据
        砂质沉积物浓度/(kg/m3)2.5
        沉积物组分细砂中砂粗砂巨砂多种组分混合后中值粒径约为300 μm
        砂质沉积物组分粒度中值/μm1503007501 500
        砂质沉积物组分含量/(质量占比,%)25501510
        沉积物浓度/(kg/m3)1.252.50.750.5
        砂质沉积物密度/(g/cm3)2.65
        底床坡度/(°)9.3×10-5冲积平原坡度平缓
        河道宽度/km3.2
        水流量/(m3/s)5 00010 00020 00040 00080 000变量
        等效流量1.563.126.2512.525.0m2/s
        参考河流南萨斯喀彻温河1.54育空河2.64亚马孙河6.47贾木纳河16.09
        演化成熟步长3 8902 9901 5001 100880
        平均水深/m4~6取决于河宽、底床坡度及水
        控制流速/(m/s)0~4.5
        网格单元大小/m50×20
        平面网格单元总量/个1 600×160
        模拟区域尺寸/m80 000×3 200
        地貌演化系数48沉积模拟加速48倍
        底床糙度-Chezy值/(m1/2/s)45经验参数
        水平涡黏系数/(m2/s)*1 000
        垂向涡黏系数/(m2/s)*0.001
        临近干网格侵蚀系数0.25
        水动力计算时间步长/min0.2
        沉积模拟数据保存时间步长/min15
        注:“—”表示无补充说明。
      • 为了分析辫状河的形态动力学和沉积特征,采用MATLAB编程对沉积模拟数据进行处理,具体分为四个步骤:(1)根据设置沉积物厚度等参数计算初始地形;(2)在模拟域内定义一个剖面,并提取所定义剖面的坐标和地形;(3)通过考虑动态侵蚀和沉积过程,将现步长与前一个步长之间的地形差作为不同步长下泥沙增厚的分布范围,计算每个模拟步骤中最终保存的沉积体的地形,作为暂存的沉降量;(4)采用Straub et al.(2012)van de Lageweg et al.(2013)提出的方法计算每个步长最终保留边界高程,绘制沉积构型剖面(Feng et al.,2025)。此外,沿剖面提取沉积物粒度数据,并绘制在沉积构型图内。根据绘制的剖面,识别和测量不同构型要素、规模、位置和倾斜度,包括增生体、心滩和辫状水道。

      • 在沉积数值模拟中,改变流量等参数得到的辫状河模拟结果虽不同,但辫状河的形态动力学与砂坝和河道的协同演化有关,使演化过程仍具有相似性(Schuurman and Kleinhas,2015)。对比五组沉积模拟,认为辫状河演化可分为四个阶段:(1)底床砂在水流的推动作用下形成了不稳定的菱形坝、舌形坝和不连续的水道;(2)顺水流的加积作用形成朵叶状、短陡迎水面和长缓背水面单元坝及相对规则、宽、浅的辫状水道;(3)多个单元坝开始顺流生长或横向“碰撞”形成不对称条形、规模中等的复合坝,此时辫状水道更加连通,形成相对复杂稳定的河网(Rice et al.,2009Schuurman et al.,2016Guo et al.,2023);(4)一种是辫状河道被充填,多个单元坝、复合坝碰撞形成规模最大的复合改造坝,另一种是复合坝被河道切割形成单元坝,与此同时,河道相对稳定,坝顶受到漫溢水流的侵蚀改造使窜沟较发育(雷涛等,2023)。

      • 流量与辫状河演化速率呈正相关,高流量辫状河演化时间在较短时间内就达到了演化成熟阶段,低流量辫状河达到相同的沉积量需要时间增加,即时长和流量成反比。以总流量为40 000 m3/s的辫状河为例(图2f),在1 100模拟步长后辫状河的形态不再发生变化,达到演化成熟阶段,而总流量为5 000 m3/s辫状河心滩增长速率越慢,需要2 000~4 000个模拟步才能达到演化成熟阶段。

        图  2  高、低流量砂质辫状河演化过程对比

        Figure 2.  Comparision of sedimentary evolution in high⁃and low⁃discharge sandy braided rivers

      • 不同流量下辫状河演化过程有明显差异,低流量的辫状河(5 000 m3/s对应南萨斯喀彻温河)演化以前三个阶段为主,演化成熟后不发育复杂改造坝,而高流量辫状河水动力强,在沉积末期易出现河道被充填或改道,会出现规模较大的复杂改造坝或规模较小的单元坝。这是因为在演化第三个阶段,低流量辫状河是以多个单元坝以顺流生长或“向下游碰撞”形成长条形心滩为主体(图3a~c),此时已形成相对稳定且连通的河道,使后续心滩合并、分割的动态平衡阶段不发育(图2a~c)。高流量的辫状河在演化的第三个阶段,多个单元坝以横向“碰撞”为主(图3d~f),单元坝和单元坝碰撞的位置往往存在次级水流,形成不对称、尾部分岔、拉长的椭圆形心滩,河道宽、深、相对连续,心滩头部受水流侵蚀,继续向下迁移,出现复合改造坝或单元坝,在心滩的中部或尾部往往发育多个窜沟(图2d~f)。

        图  3  高、低流量砂质辫状河心滩演化模式对比(Guo et al.,2023

        Figure 3.  Comparision of Channel bar evolution model in high⁃and low⁃discharge sandy braided river evolution (Guo et al., 2023)

      • 在对照模拟组中发现不同流量模拟辫状河演化过程中心滩的规模存在较明显的差异。结果表明流量越大,则心滩坝越长(图4a)。

        图  4  不同流量下砂质辫状河演化过程中砂坝规模演变趋势

        Figure 4.  Trend in sandbar size evolution during the evolution of sandy braided rivers under different discharge conditions

        不同于平均坝长,平均坝宽虽然存在演化速率的差异,但最终平均坝宽接近(图4b)。通过坝长宽比表明,流量越高,则坝长/宽比值增加越快(图4c)。心滩坝平均面积直接指示了心滩规模(图4d),随着模拟中流量的增加,心滩坝面积也快速增加,表明心滩坝规模与流量存在明显正相关关系。总体来说,演化过程中流量与砂坝的长度、面积成正比,但与宽度无正相关关系,说明流量越大,流速激增使砂坝更易受到切割。与此同时,水流的集中和加速又进一步加快了坝之间侧向迁移的速度,从而生成更多的复合坝,在这种“水道切割—砂坝充填”循环协同演化过程中呈现互相调节的双重效果,使得砂坝总体宽度相对固定,长宽比介于3.0~4.5。

      • 地下储层往往是由辫状河演化自然终止的砂质沉积物组成,其中大部分沉积物已经达到成熟状态。不同流量砂质辫状河模拟的成熟阶段河道和心滩沉积特征差异化研究(图56),对地下储层地质学理论的完善具有现实意义。

        图  5  不同流量条件对照模拟组砂质辫状河剖面

        Figure 5.  Simulated profiles of sandy braided rivers at different discharge conditions

        图  6  不同流量成熟时期辫状水道分析图

        Figure 6.  Analysis diagrams of braided channels at different discharge mature stages

      • 辫状河流量不仅影响辫状水道的平面样式、结构特征,更直接控制辫状水道的宽深。剖面河床形态对比发现,低流量辫状河的辫状水道形态通常为U形、V形(图5a),随着流量增加,流速激增,水流集中和加速使更大型的宽阔、平坦的水道逐步发育(图5e);从平面上看,流量越大,辫状水道形态越顺直,辫状水道数量越少、密度越小,但其宽度也越大。低流量辫状河以窄而多的交织河网为主,高流量的辫状河以宽而数量有限的骨干辫状水道网络为主(图27)。

        图  7  不同流量下模拟砂质辫状河演化成熟时期沉积地貌图

        Figure 7.  Depositional geomorphology at mature stages of simulated sandy braided rivers under different discharge

      • 以水流量为5 000 m3/s辫状河为例,切物源方向剖面上平均存在7.5条辫状水道,单个辫状水道宽度20~50 m,弯曲度可达1.3,辫状水道分岔系数在2左右;随着流量增大,辫状水道数量逐步降低至3条,宽度增大为200~600 m,弯曲度减小至1.1左右,分岔系数降为1.2左右。通过统计模拟每个网格水深数据,分析不同流量演化成熟阶段辫状河所有网格水深占比发现,随着流量增加,河道水深分布由单峰偏态逐步转变为双峰弱偏态(图6a),显示流量越大,水道间流速的差异越大,水流多集中于宽阔、对称的辫状水道内。说明更大的水流量代表更强的水动力条件,直接控制辫状河水深,导致辫状水道深度也与水流量呈现正相关关系。

      • 流量的差异不仅主导了辫状水道规模、样式、结构等特征,同时也塑造了不同样式、规模的心滩。在演化成熟阶段存在四种类型的砂坝:单元坝、复合坝、复合改造坝、河心滩(Sambrook Smith et al.,2006,2009;Puig et al.,2019冯文杰等,2025)。其中低流量辫状河以规模较小的复合坝为沉积主体(图8a),可见河边滩和单元坝;随着流量增加,出现大型的复合改造坝(图8b),但此时规模较小的复合坝仍为主体。当流量为20 000 m3/s,辫状河内部坝侧向迁移横向“碰撞”变得更加连片,复合改造坝更宽、规模更大(图8c);当流量为40 000 m3/s时受水流改造,复合改造坝多为长条形(图8d);当流量为80 000 m3/s,流速激增,大部分坝被改造,水道更宽、连通(图8d)。

        图  8  不同流量成熟时期砂坝参数

        Figure 8.  Sandbar parameters at mature stages for different discharge

        研究表明,砂坝所处的位置不同,水流样式不同,其中单元坝、复合坝位于河床中部,两侧被不均衡的辫状水道环绕;而河边滩则位于河道边部,一侧与河岸相接,另一侧受辫状水道控制。值得注意的是,辫状河的河边滩与曲流河的边滩存在差异性。辫状河的河边滩与曲流河的边滩位置相似,但沉积演化过程存在明显差异,使得内部结构也存在不同。辫状河的河边滩沉积演化早期受双侧水流影响仍显示为辫状河砂坝特征,在演化中后期先侧向迁移形成长条形砂坝,砂坝阻碍水流使得双侧水流转变为单侧水流,之后砂坝向下迁移形成河边滩,因此自下而上发育侧积体、垂向加积体或向心充填的河道沉积构成,内部夹层不连续不稳定。而曲流河边滩始终受到单侧水流影响,发育侧积体,内部夹层稳定,相对连续(郭颖等,2024)。

      • 流量对心滩形态样式有至关重要的作用。对照模拟结果显示:在低能辫状河中(图7a),由于辫状水道网络发育程度高、以窄而多的交织河网为整体样式,在其控制下,辫状河演化末期心滩被河道分岔或合并的动态平衡阶段相对不发育,形成的心滩呈简单、叶状的单元坝和长条状的复合坝均匀分布特征。随着河流流量增加,心滩两侧辫状水道规模出现明显的主次、强弱差异并导致心滩不对称性增强,河道两侧水流量的差异使得单元坝更多以“横向碰撞”形成复合坝、规模较大的复合改造坝、河边滩,心滩面积逐渐增大,形态样式变得更复杂。

      • 流量控制着不同类型心滩的规模与定量关系演变。结果表明,流量越大,砂坝经水流改造动力越强,砂坝总面积占比越小(图8a),其中流量越大,复合改造坝占比越多,长宽比越小,复合坝面积占比越少,长宽比减少,河边滩面积占比越少,长宽比增加(图8b)。例如:当流量从5 000 m3/s增加到10 000 m3/s时,复合坝长宽比变小、面积占比减少,复合改造坝占比增加,说明随着流量增加,更多单元坝和复合坝组合形成复合改造坝;当流量从40 000 m3/s增加到80 000 m3/s时,复合改造坝长宽比急速减小,单元坝、复合坝长宽比缓慢减小,说明流量增加时,更多的复合改造坝被河道分割,单元坝向下迁移或与其他砂坝组合形成复合坝或河边滩。在演化成熟阶段,单元坝、复合坝长宽比主要分布在3~3.5,复合改造坝、河边滩长宽比变化较大。因此,辫状河流量高低影响了心滩的形态、规模、结构样式等特征,在实际储层解剖研究中应加以区分考虑。

      • 低流量辫状河心滩内部一般存在1~3条小规模的窜沟,这些窜沟将心滩“分割”为多个顶部分割、底部相连的次级心滩;随着流量增大,辫状河心滩内部窜沟的发育程度逐渐增大,包括窜沟数量、规模、深度均逐步增大。窜沟作为辫状水道侧缘附属的水流分散通道,在辫状河沉积过程中扮演了重要的角色:窜沟冲刷心滩表面,引起心滩顶部沉积物充分分布,形成复杂的窜沟网络(图9)。在心滩内部,窜沟发育程度越高,则其顶部沉积构型越复杂,在油气资源开发中会形成上下部流场差异。

        图  9  不同流量下模拟砂质辫状河演化成熟阶段沉积微相分布图

        Figure 9.  Distribution of sedimentary microfacies at mature stages for simulated sandy braided rivers in different flows

      • 演化成熟后心滩经历了较为复杂的沉积演化过程,从而形成了复杂的内部构型。分别对总流量为5 000 m3/s的低流量辫状河和总流量为40 000 m3/s的高流量辫状河单元坝、复合坝进行解剖,探究高、低流量辫状河下单元坝和复合坝内部多期次增生体的类型、规模、倾角的差异性,获取内部的结构信息,再现了高低流量辫状河内部单元坝和复合坝的叠置样式和构型特征,对指导地下油气精细开发具有重要的实际意义。

      • 通过对比发现,演化成熟的单元坝垂向自下而上均发育顺流加积体、侧积体、垂向加积体,但低流量辫状河内部单元坝更简单,不发育窜沟,顺流加积体为沉积主体,这与前人提出心滩底部发育大型向下游迁移的顺流加积底型的认识一致(廖保方等,1998Best et al.,2003牛博等,2015张宪国等,2020何维领等,2024唐洪等,2026)。而高流量辫状河内单元坝更复杂、尾部多分支,发育1~3个充填或半充填的小窜沟,将其分割成多个双向水流的次级砂坝。

      • 低流量辫状河演化成熟阶段,平面上单元坝规模小,呈长条状、尾部拉长,不发育窜沟。剖面上不同的颜色代表不同时期形成的增生体,单元坝发育五期增生体,自下而上发育顺流加积体、侧积体、垂向加积体,以顺流加积体为主。早期生成的顺流加积体位于单元坝的底部(图10a①~③号),是沉积的主体。侧积体主要位于单元坝两侧,倾角大(平均2.8°),顺物源方向可观察到沉积后期形成的增生体向下游方向迁移的过程(图10b中④~⑤号),侧积体往往发育在单元坝的尾部和两侧;垂向加积体往往位于单元坝切物源方向的顶部,占比较低。

        图  10  低流量辫状河演化成熟阶段单元坝沉积构型

        Figure 10.  Sedimentary architecture of unit bar developed with low discharge braided river at mature stages

        南萨斯喀彻温砂质辫状河单元坝规模小,最小宽约100 m,最大长度300 m,厚1.5~3.0 m,以向下迁移为主。通过现代沉积和探地雷达对比认为单元坝底部以板状或槽状交错层理的顺流加积体为主体,被形象地称为“沙核”(Best et al.,2003),迎水面顶部发育较少的高角度倾斜的沉积物,沉积模拟图10b③号增生体特征符合这一特征。

      • 高流量辫状河演化成熟阶段,平面上单元坝呈叶状、尾部多分支,发育1条窜沟。剖面上单元坝发育九期增生体,自下而上发育顺流加积体、侧积体、垂向加积体(图11)。通过统计顺流加积体倾角可知水流方向是从单元坝的两侧向中心迁移,迁移的过程中倾角越来越陡,介于0.7°~4.6°。侧积体具有下切特点,倾角大(平均为4.5°),往往发育在单元坝的头尾和水浅一侧(图11c中C),方向以上游向下游、水浅一侧(图11c中C)向水深一侧(图11c中C’)为主。窜沟为典型的向心充填,倾角逐渐减小,呈自上游向下游的倾斜方向。垂向加积体往往是沉积后期以较小的倾角覆盖在单元坝的顶部。前人通过现代沉积、探槽、探地雷达等方式对Jamuna河典型砂质辫状河心滩进行解剖,认为在心滩坝底部发育大型平缓顺流加积底型,横向多期增生体加积的沉积模式(Best et al.,2003牛博等,2015),该观点与本文中心滩垂向自下而上发育顺流加积体、侧积体、垂向加积体一致。不同点在于,受窜沟的影响,窜沟内部发育侧积体,使得心滩在水平方向以顺流加积体到侧积体的形式重复出现,与刘警阳等(2024)提出高流量游荡型辫状河受洪水影响心滩易冲裂,发育窜沟的认识一致。

        图  11  高流量辫状河演化成熟阶段单元坝沉积构型

        Figure 11.  Sedimentary architecture of unit bar developed with high discharge braided river at mature stages

      • 对比发现高流量辫状河复合坝比低流量辫状河复合坝叠切关系更复杂、窜沟更发育。

      • 低流量辫状河复合坝呈尾部多分支的不对称叶状,发育3条窜沟(图12)。剖面上单元坝发育五期增生体,自下而上发育顺流加积体、侧积体、垂向加积体,以顺流加积体为主,顺物源方向可观察到沉积早期形成的①~③期顺流加积体向下游方向迁移的过程(图12c中400~700 m处)。④~⑤号增生体为侧积体,河道下切作用不明显,倾角小(平均2.8°)往往发育在单元坝的四周和顶部窜沟内;垂向加积体往往位于单元坝切物源方向的顶部,占比较少。

        图  12  低流量辫状河演化成熟阶段复合坝沉积构型

        Figure 12.  Sedimentary architecture of compound bar developed with low discharge braided river at mature stages

        南萨斯喀彻温砂质辫状河复合坝长约800 m,宽度约400 m。底部2.5~3 m深可见明显侵蚀面的不连续板状或槽状交错层理的顺流加积体(占比63%)(Sambrook et al.,2006),在砂坝顶部0.5 m内低角度反射的增生体占比高达84%,这与沉积模拟中底部发育顺流加积体、顶部发育低角度加积体一致。除此之外,复合坝相之间存在明显的空间变异性,在背水面靠近深泓线部位(类似图12d中C’位置),存在大角度倾斜沉积体,在沉积模拟中也可观察到相同位置出现高角度侧积体。

      • 高流量辫状河复合坝呈尾部多分支的不对称椭圆状,窜沟发育程度高(图15)。单元坝发育十期增生体,在中部切物源剖面上底部存在大规模的顺流加积体,向上逐渐变为薄层侧积体或“向心充填”河道沉积组合(图13b中④号),说明侧积体往往为沉积晚期发育在单元坝的头尾、水浅一侧(图13e中C’)及窜沟内。因此,高流量复合坝自下而上发育顺流加积体、侧积体或河道充填、垂向加积体或窜沟。对比Jamuna河典型砂质辫状河解剖结果,除了具备单元坝的特点以外,在坝头、尾和水浅一侧发育小型垂厚2 m的多组槽状交错反射雷达相特征(Best et al.,2003),向坝中心倾角逐渐减小,这与我们统计复合坝的侧积体中位数3 m规模相近。

        图  13  高流量辫状河演化成熟阶段复合坝沉积构型

        Figure 13.  Sedimentary architecture of compound bar developed with high discharge braided river at mature stages

        图  14  高、低流量辫状河单元坝、复合坝沉积构型规模参数对比

        Figure 14.  Comparison of sedimentary architecture scale parameters in unit bar and compound bar of high⁃and low⁃discharge braided river

        图  15  高流量辫状河单元坝与野外露头沉积构型对比

        Figure 15.  Comparison of sedimentary architecture between highdischarge braided river sandbar and rock outcrop

        通过对比发现,无论是倾角、高度、长、宽等规模参数,高流量辫状河单元坝、复合坝的规模参数均明显大于低流量辫状河的单元坝、复合坝(图14),其中增生体倾角、高、长、宽等规模参数呈明显相关关系,低流量增生体长宽比、宽厚比中位数分别为5和18.7,高流量增生体长宽比、宽厚比中位数分别为3和13.6,说明流量与增生体规模呈正相关,但长宽比和宽厚比与流量呈负相关,这与Li et al.(2023)提出的坝高、坝面积与流量呈正相关一致。除此之外,同一个砂坝内侧积体、顺流加积体、垂向加积体倾角逐渐减小,高度依次降低。从规模上看,同一个砂坝内顺流加积体和垂向加积体的长度和宽度明显高于侧积体。因此,侧积体位于心滩两侧,具有高倾角、较厚和小规模的特点,垂向加积体位于心滩顶部,具有较长的延伸范围、薄厚度、小倾角的特点。值得注意的是,受水流强弱的影响,内部增生体倾角是变化的,尤其是顺流加积、侧积体倾角离散程度高,分布散,这使得地下解剖难度增大。

      • 通过对比野外露头解剖结果,验证模拟结果的可靠性(图15)。大同云冈组石窟段为典型的砂质辫状河沉积,图15c2剖面发育两期心滩,认为单期心滩的厚度介于1.5~7.5 m,宽50~200 m,底部发育大型顺流加积底型,顶部偶见小规模窜沟,心滩内可见1~3 cm间歇性薄落淤层(任晓旭,2020),这与高流量砂质辫状河内单元坝的内部结构相似。值得注意的是,不同位置落淤层的倾角是变化的。例如,2号和4号增生体之间发育落淤层倾角近水平(图15b中C侧),但在1、2、3增生体之间落淤层倾角倾斜(图15b中C’侧)。

      • 在沉积模拟中提取粒度D50剖面数据和心滩不同位置处D50数据(图16),与单井实际案例分析对比,验证模拟结果的可靠性。心滩的储层非均质性与不同类型增生体的叠置样式有关,进一步明确心滩是以复合反韵律为主,头部发育多套反韵律顺流加积体,心滩中部发育窜沟,偶见正韵律侧积体,内部增生体界面处易发育细粒薄夹层的落淤层。

        图  16  高流量辫状河单元坝与井下实钻案例沉积构型对比

        Figure 16.  Comparison of sedimentary architecture between high⁃discharge braided river sandbar and drill core data

        真实的辫状河发育过程受自然界复杂因素的影响,包括沉积物粒度分选性(如物源供给差异)、古地形坡度及流量变化(洪峰或季节性河流)等多因素耦合影响。分析表明流量变化对砂质辫状河的形态、规模及内部结构呈主导作用。为建立明晰的流量和辫状河规模半定量—定量相关关系,本文重点讨论了以流量为单一变量下的辫状河定量解析,关于其他边界条件例如粒度、坡度的影响,将在后续研究中进行系统性探讨。

      • (1) 流量与辫状河的沉积演化速率成正比。低流量辫状河演化程度低,砂坝以向下迁移为主,不易被河道分割或再次合并;而高流量辫状河演化程度高,砂坝以侧向迁移为主,心滩合并、被分割的动态平衡过程较普遍。

        (2) 流量影响辫状水道平面样式、水深,并进一步影响辫状水道和心滩的协同演化过程。流量越大,河道越宽、直、深,分岔少,演化成熟阶段辫状水道网络更复杂;流量越大,心滩样式更复杂,出现规模较大的复合改造坝和河边滩。由于流量增加,心滩经水流改造动力越强,其内部叠切关系更复杂,更易发育窜沟。除此之外,流量越大的辫状河,在演化成熟阶段复合改造坝占比增加,复合坝、河边滩面积占比减少。

        (3) 砂质辫状河心滩内部结构具有很强的空间变异性。相同点在于自下而上均发育顺流加积体、侧积体、垂向加积体,不同点在于流量越大,心滩内部侧积体占比增加。其中顺流加积体(DA)作为沉积主体单元,侧积体位于两侧或窜沟内,具备“短轴、厚层(3~8 m)、高倾角(>4°)”的特点,垂向加积体位于心滩顶部,具备“长轴、薄层、低倾角(<1°)”的特点。

        (4) 流量与增生体规模呈正相关,但与长宽比、宽厚比呈负相关。流量越大,心滩内部增生体的倾角、高度、长、宽等参数越大,但宽度增长有限,其中低流量辫状河心滩长宽比中位数为3.5,高流量辫状河心滩长宽比中位数为4.5。

    参考文献 (40)

    目录

      /

      返回文章
      返回