高级搜索

留言板

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

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

构造掀斜主导的断陷湖盆缓坡层序“源—汇”正演模拟定量研究

谈明轩 朱筱敏 张自力 孙浩南

谈明轩, 朱筱敏, 张自力, 孙浩南. 构造掀斜主导的断陷湖盆缓坡层序“源—汇”正演模拟定量研究[J]. 沉积学报, 2022, 40(6): 1481-1493. doi: 10.14027/j.issn.1000-0550.2022.061
引用本文: 谈明轩, 朱筱敏, 张自力, 孙浩南. 构造掀斜主导的断陷湖盆缓坡层序“源—汇”正演模拟定量研究[J]. 沉积学报, 2022, 40(6): 1481-1493. doi: 10.14027/j.issn.1000-0550.2022.061
TAN MingXuan, ZHU XiaoMin, ZHANG ZiLi, SUN HaoNan. Source-to-sink Quantitative Stratigraphic Forward Modeling on the Tilted Hanging-wall Sequence Architecture of a Tectonically-driven Lacustrine Rift Basin[J]. Acta Sedimentologica Sinica, 2022, 40(6): 1481-1493. doi: 10.14027/j.issn.1000-0550.2022.061
Citation: TAN MingXuan, ZHU XiaoMin, ZHANG ZiLi, SUN HaoNan. Source-to-sink Quantitative Stratigraphic Forward Modeling on the Tilted Hanging-wall Sequence Architecture of a Tectonically-driven Lacustrine Rift Basin[J]. Acta Sedimentologica Sinica, 2022, 40(6): 1481-1493. doi: 10.14027/j.issn.1000-0550.2022.061

构造掀斜主导的断陷湖盆缓坡层序“源—汇”正演模拟定量研究

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

江苏省自然科学基金 BK20200529

中央高校基本业务费 B220202047

详细信息
  • 中图分类号: P618.13

Source-to-sink Quantitative Stratigraphic Forward Modeling on the Tilted Hanging-wall Sequence Architecture of a Tectonically-driven Lacustrine Rift Basin

Funds: 

National Natural Science Foundation of Jiangsu Province BK20200529

Fundamental Research Funds for the Central Universities B220202047

  • 摘要: 断陷湖盆缓坡带的构造掀斜作用对其地层叠置样式和沉积体系展布具有重要的控制作用。通过基于“源—汇”系统理论的地层正演模拟研究方法,在三级层序地层格架下(以2 Myr为总时长)对断陷盆地缓坡带构造掀斜型层序地层发育样式及其相关地貌演化、沉积物通量变化进行了定量表征。以恒定的绝对湖平面为研究前提,共建立了高/低沉积物供给的背景下4种不同的地层正演模型。两种不同缓坡带掀斜初始位移模型的设定反映了掀斜点和沉积滨线的相对位置,以便于不同掀斜层序样式对比分析。研究表明,掀斜支点位于滨线之上的模型常形成湖侵和正常湖退,而掀斜支点位于滨线之下的模型常形成强制湖退。滨线迁移轨迹和层序叠置样式的不一致性贯穿整个层序发育过程,是沉积物通量及滨线迁移共同作用的结果。该模型中沉积物通量的变化直接响应于源区地貌演变,滨线迁移则受控于新增可容空间和沉积物供给的相互关系。这一系列构造掀斜型层序模型及其相关成因解释对断陷盆地层序地层格架及其叠置样式的精细解剖具有重要的指导意义。该模拟研究对断陷湖盆的高频湖平面变化特征也有较为合理的解释(如海洋同位素3期以来的贝加尔湖)。然而本次正演模拟所设置的初始地貌较为单一,在相同构造、气候条件下和不同沉积物供给背景下所形成地层模型无法涵盖所有缓坡带层序地层类型。由此可见,基于“源—汇”系统思想探讨断陷盆地层序地层多样性的相关模拟工作仍然任重道远。
  • 图  1  构造掀斜作用下恒定绝对湖平面的缓坡带滨线迁移模式(据文献[9]修改)

    Figure  1.  Schematic model showing shoreline migration influenced by tilting with a fixed absolute lake level in the hanging⁃wall dipslope (modified from reference [9])

    Fig.1

    图  2  简化的缓坡带初始三维古地形

    Figure  2.  Simplified initial three⁃dimensional paleo⁃topography of the hanging⁃wall dipslope

    Fig.2

    图  3  两种掀斜构造模型的剖面示意及三维立体特征

    Figure  3.  Schematic profiles of two tectonic tilting models and three⁃dimensional characteristics

    Fig.3

    图  4  模型1地层剖面、平面模型及相关沉积物通量变化特征

    Figure  4.  Stratigraphic profile and planar model of model run 1 and variation of sediment flux

    Fig.4

    图  5  模型2地层剖面、平面模型及相关沉积物通量变化特征

    Figure  5.  Stratigraphic profile and planar model of model run 2 and variation of sediment flux

    Fig.5

    图  6  模型3地层剖面、平面模型及相关沉积物通量变化特征

    Figure  6.  Stratigraphic profile and planar model of model run 3 and variation of sediment flux

    Fig.6

    图  7  模型4地层剖面、平面模型及相关沉积物通量变化特征

    Figure  7.  Stratigraphic profile and planar model of model run 4 and variation of sediment flux

    Fig.7

    图  8  不同模型中滨岸带新增可容空间(δA)与沉积物供给(δS)变化特征

    Figure  8.  Variation of newly⁃created accommodation space (δA) and sediment supply (δS) in the coastal area

    Fig.8

    图  9  构造掀斜作用下缓坡带湖平面变化新模式

    Figure  9.  New schematic model of lake level fluctuation of the hanging⁃wall dipslope driven by tectonic tilting

    Fig.9

    图  10  (a,b)贝加尔湖Selenga三角洲地理位置;(c)中—晚更新世以来贝加尔湖(相对及绝对)湖平面变化特征;(d)基于地震解释成果的地层剖面E—E’;(e)贝加尔湖Selenga三角洲发育区NW区块及SE区块相对湖平面变化特征(修改自文献[33,36])

    Figure  10.  (a, b) Location of Selenga delta in Lake Baikal; (c) (relative and absolute) lake level fluctuation of Lake Baikal since the Middle to Upper Pleistocene; (d) stratigraphic profile E⁃E’ based on seismic interpretation; and (e) relative lake levels of NW and SE blocks of the Selenga delta in Lake Baikal (modified from references [33,36])

    Fig.10

    表  1  地层正演模拟关键参数基本设置

    模型构造模型侵蚀系数/(yr-1地表扩散系数/(m2/yr)水下扩散系数/(m2/yr)河流扩散系数/(m2/yr)
    1A1×10-50.300.505.00
    2B1×10-50.300.505.00
    3A1×10-60.300.505.00
    4B1×10-60.300.505.00
    下载: 导出CSV
  • [1] Leeder M R, Gawthorpe R L. Sedimentary models for extensional tilt-block/half-graben basins[J]. Geological Society, London, Special Publications, 1987, 28(1): 139-152.
    [2] Clarke P, Parnell J. Facies analysis of a back-tilted lacustrine basin in a strike-slip zone, Lower Devonian, Scotland[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 1999, 151(1/2/3): 167-190.
    [3] Biggs J, Nissen E, Craig T, et al. Breaking up the hanging wall of a rift-border fault: The 2009 Karonga earthquakes, Malawi[J]. Geophysical Research Letters, 2010, 37(11): L11305.
    [4] Fraser G D, Witkind I J, Nelson W H. A geological interpretation of the epicentral area-the dual-basin concept[R]. U.S. Geological Survey Professional Paper, 1964, 435: 99-106.
    [5] Gawthorpe R L, Leeder M R. Tectono-sedimentary evolution of active extensional basins[J]. Basin Research, 2000, 12(3/4): 195-218.
    [6] Pechlivanidou S, Cowie P A, Duclaux G, et al. Tipping the balance: Shifts in sediment production in an active rift setting[J]. Geology, 2019, 47(3): 259-262.
    [7] 刘震,曾宪斌,张万选. 构造掀斜对单断湖盆湖平面变化的影响[J]. 沉积学报,1997,15(4):64-66.

    Liu Zhen, Zeng Xianbin, Zhang Wanxuan. Effect of tectonic tilt of faulted blocks on lake level change of single faulted continental basins[J]. Acta Sedimentologica Sinica, 1997, 15(4): 64-66.
    [8] 任建业,陆永潮,张青林. 断陷盆地构造坡折带形成机制及其对层序发育样式的控制[J]. 地球科学:中国地质大学学报,2004,29(5):596-602.

    Ren Jianye, Lu Yongchao, Zhang Qinglin. Forming mechanism of structural slope-break and its control on sequence style in faulted basin[J]. Earth Science: Journal of China University of Geosciences, 2004, 29(5): 596-602.
    [9] Wu H, Ji Y L, Wu C L, et al. Stratigraphic response to spatiotemporally varying tectonic forcing in rifted continental basin: Insight from a coupled tectonic-stratigraphic numerical model[J]. Basin Research, 2019, 31(2): 311-336.
    [10] 阳孝法,林畅松,彭莉,等. 幕式构造沉降过程中可容纳空间的变化规律:以掀斜式半地堑盆地为例[J]. 大庆石油地质与开发,2010,29(4):18-22.

    Yang Xiaofa, Lin Changsong, Peng Li, et al. Evolution laws of the accommodation space in the course of episodic tectonic subsidence: A case from tilted half-graben basin[J]. Petroleum Geology and Oilfield Development in Daqing, 2010, 29(4): 18-22.
    [11] Gawthorpe R L, Fraser A J, Collier R E L. Sequence stratigraphy in active extensional basins: Implications for the interpretation of ancient basin-fills[J]. Marine and Petroleum Geology, 1994, 11(6): 642-658.
    [12] 林畅松,夏庆龙,施和生,等. 地貌演化、源—汇过程与盆地分析[J]. 地学前缘,2015,22(1):9-20.

    Lin Changsong, Xia Qinglong, Shi Hesheng, et al. Geomorphological evolution, source to sink system and basin analysis[J]. Earth Science Frontiers, 2015, 22(1): 9-20.
    [13] 谈明轩,朱筱敏,张自力,等. 古“源—汇”系统沉积学问题及基本研究方法简述[J]. 石油与天然气地质,2020,41(5):1107-1118.

    Tan Mingxuan, Zhu Xiaomin, Zhang Zili, et al. Summary of sedimentological issues and fundamental approaches in terms of ancient “source-to-sink” systems[J]. Oil & Gas Geology, 2020, 41(5): 1107-1118.
    [14] 陆威延,朱红涛,徐长贵,等. 源—汇系统级次划分方法及应用[J]. 地球科学,2020,45(4):1327-1336.

    Lu Weiyan, Zhu Hongtao, Xu Changgui, et al. Methods and applications of level subdivision of source-to-sink system[J]. Earth Science, 2020, 45(4): 1327-1336.
    [15] 邵龙义,王学天,李雅楠,等. 深时源—汇系统古地理重建方法评述[J]. 古地理学报,2019,21(1):67-81.

    Shao Longyi, Wang Xuetian, Li Yanan, et al. Review on palaeogeographic reconstruction of deep-time source-to-sink systems[J]. Journal of Palaeogeography, 2019, 21(1): 67-81.
    [16] 操应长,徐琦松,王健. 沉积盆地“源—汇”系统研究进展[J]. 地学前缘,2018,25(4):116-131.

    Cao Yingchang, Xu Qisong, Wang Jian. Progress in “source-to-sink” system research[J]. Earth Science Frontiers, 2018, 25(4): 116-131.
    [17] Ding X S, Salles T, Flament N, et al. Drainage and sedimentary responses to dynamic topography[J]. Geophysical Research Letters, 2019, 46(24): 14385-14394.
    [18] Ding X S, Salles T, Flament N, et al. Quantitative stratigraphic analysis in a source-to-sink numerical framework[J]. Geoscientific Model Development, 2019, 12(6): 2571-2585.
    [19] Salles T, Ding X S, Brocard G. pyBadlands: A framework to simulate sediment transport, landscape dynamics and basin stratigraphic evolution through space and time[J]. PLoS One, 2018, 13(4): e0195557.
    [20] 刘泽,李三忠, Bukhari S W H,等. 动态古地貌再造:Badlands软件在盆地分析中的应用[J]. 古地理学报,2020,22(1):29-38.

    Liu Ze, Li Sanzhong, Bukhari S W H, et al. Reconstruction of dynamic palaeogeomorphy: Application of Badlands software in basin analysis[J]. Journal of Palaeogeography, 2020, 22(1): 29-38.
    [21] Burgess P M. A brief review of developments in stratigraphic forward modelling, 2000–2009[M]//Roberts D G, Bally A M. Regional geology and tectonics: Principles of geologic analysis. Amsterdam: Elsevier, 2012: 378-404.
    [22] Huang X, Griffiths C M, Liu J. Recent development in stratigraphic forward modelling and its application in petroleum exploration[J]. Australian Journal of Earth Sciences, 2015, 62(8): 903-919.
    [23] Chen A, Darbon J, Morel J M. Landscape evolution models: A review of their fundamental equations[J]. Geomorphology, 2014, 219: 68-86.
    [24] Clevis Q, de Boer P, Wachter M. Numerical modelling of drainage basin evolution and three-dimensional alluvial fan stratigraphy[J]. Sedimentary Geology, 2003, 163(1/2): 85-110.
    [25] 操应长,姜在兴,夏斌. 幕式差异沉降运动对断陷湖盆中湖平面和水深变化的影响[J]. 石油实验地质,2003,25(4):323-327.

    Cao Yingchang, Jiang Zaixing, Xia Bin. Effects of episodic differential subsidence on the changes of lake level and water depth in rift lake basins[J]. Petroleum Geology & Experiment, 2003, 25(4): 323-327.
    [26] Stock J D, Montgomery D R. Geologic constraints on bedrock river incision using the stream power law[J]. Journal of Geophysical Research: Solid Earth, 1999, 104(B3): 4983-4993.
    [27] Syvitski J P M, Kettner A J. Scaling sediment flux across landscapes[C]//Proceedings of the sediment dynamics in changing environments. New Zealand: IAHS, 2008: 149-156.
    [28] Emmel B, de Jager G, Zieba K, et al. A 3D, map based approach to reconstruct and calibrate palaeo-bathymetries – Testing the Cretaceous water depth of the Hammerfest Basin, southwestern Barents Sea[J]. Continental Shelf Research, 2015, 97: 21-31.
    [29] 林畅松,李思田,任建业. 断陷湖盆层序地层研究和计算机模拟:以二连盆地乌里雅斯太断陷为例[J]. 地学前缘,1995,2(3/4):124-132.

    Lin Changsong, Li Sitian, Ren Jianye. Sequence architecture and depositional systems of Erlian lacustrine fault basins, North China[J]. Earth Science Frontiers, 1995, 2(3/4): 124-132.
    [30] Bohacs K M, Carroll A R, Neal J E, et al. Lake-Basin type, source potential, and hydrocarbon character: An integrated sequence-stratigraphic-geochemical framework[M]//Gierlowski-Kordesch E H, Kelts K R. Lake basins through space and time. AAPG Studies in Geology 46, Tulsa: AAPG, 2000: 3-34.
    [31] Cross N E, Bosence D W J. Tectono-sedimentary models for rift-basin carbonate systems[M]//Lukasik J, Simo J A. Controls on carbonate platform and reef development. SEPM Special Publication 89, Tulsa: SEPM, 2008: 83-105.
    [32] Scholz C A, Moore T C, Hutchinson D R, et al. Comparative sequence stratigraphy of low-latitude versus high-latitude lacustrine rift basins: Seismic data examples from the East African and Baikal rifts[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 1998, 140(1/2/3/4): 401-420.
    [33] Mats V D, Fujii S, Mashiko K, et al. Changes in lake baikal water levels and runoff direction in the Quaternary Period[M]//Minoura K. Lake baikal. Amsterdam: Elsevier, 2000: 15-34.
    [34] Mats V D, Perepelova T I. A new perspective on evolution of the Baikal Rift[J]. Geoscience Frontiers, 2011, 2(3): 349-365.
    [35] Colman S M. Water-level changes in Lake Baikal, Siberia: Tectonism versus climate[J]. Geology, 1998, 26(6): 531-534.
    [36] Urabe A, Tateishi M, Inouchi Y, et al. Lake-level changes during the past 100, 000 years at Lake Baikal, southern Siberia[J]. Quaternary Research, 2004, 62(2): 214-222.
    [37] Osipov E Y, Khlystov O M. Glaciers and meltwater flux to Lake Baikal during the Last Glacial Maximum[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2010, 294(1/2): 4-15.
    [38] Katsuta N, Ikeda H, Shibata K, et al. Hydrological and climate changes in southeast Siberia over the last 33 kyr[J]. Global and Planetary Change, 2018, 164: 11-26.
  • [1] 朱筱敏, 刘强虎, 谈明轩, 李顺利, 陈贺贺, 聂银兰.  深时源—汇系统综合研究和沙垒田实例分析 . 沉积学报, 2023, 41(6): 1781-1797. doi: 10.14027/j.issn.1000-0550.2023.104
    [2] 李丹, 陈洪德, 刘磊, 王志伟, 赵伟波, 虎建玲.  鄂尔多斯地区南部早二叠世山西组源-汇系统差异与构造-沉积格局 . 沉积学报, 2023, (): -. doi: 10.14027/j.issn.1000-0550.2023.093
    [3] 刘彦锋, 段太忠, 龚伟, 廉培庆, 张文彪, 黄渊.  基于深度学习的地层沉积正演模拟代理模型构建与应用 . 沉积学报, 2023, 41(3): 791-803. doi: 10.14027/j.issn.1000-0550.2021.152
    [4] 梁杰, 刘杰, 牛胜利, 林秋金, 朱君雪, 刘徐敏.  珠江口盆地惠州25转换带文昌组五段早期低位体系域源—渠—汇耦合关系及勘探潜力探讨 . 沉积学报, 2022, 40(6): 1451-1460. doi: 10.14027/j.issn.1000-0550.2022.016
    [5] 王爱军, 叶翔, 徐晓晖, 谢津剑, 陶舒琴, 赖志坤, 吴水兰, 杨雨欣.  亚热带中小型山溪性河流—宽陆架系统“源—汇”过程 . 沉积学报, 2022, 40(6): 1615-1634. doi: 10.14027/j.issn.1000-0550.2022.135
    [6] 刘炳强, 王伟超, 张文龙, 黄曼, 孙玉琦, 邵龙义.  陆相盆地河—湖沉积源—汇系统收支分析 . 沉积学报, 2022, 40(6): 1494-1512. doi: 10.14027/j.issn.1000-0550.2022.053
    [7] 马艺萍, 王荣华, 戴霜, 马晓军.  北祁连北大河沉积物碎屑组成及物源正演分析 . 沉积学报, 2022, 40(6): 1525-1541. doi: 10.14027/j.issn.1000-0550.2022.086
    [8] “源—汇系统”专栏 . 沉积学报, 2022, 40(6): -.
    [9] 袁静, 钟剑辉, 宋明水, 张宇, 向奎, 赵永福, 俞国鼎, 李欣尧.  沾化凹陷孤岛西部斜坡带沙三段重力流沉积特征与源-汇体系 . 沉积学报, 2018, 36(3): 542-556. doi: 10.14027/j.issn.1000-0550.2018.076
    [10] 王圣柱, 吴倩倩, 程世伟, 薛雁, 陈平.  准噶尔盆地北缘哈山构造带油气输导系统与运聚规律 . 沉积学报, 2017, 35(2): 405-412. doi: 10.14027/j.cnki.cjxb.2017.02.017
    [11] 杨伟利, 王毅, 张众琛, 梁晨, 李京昌.  塔里木盆地顺南地区中-下奥陶统正演沉积数值模拟 . 沉积学报, 2017, 35(1): 117-123. doi: 10.14027/j.cnki.cjxb.2017.01.012
    [12] 远源缓坡型薄层细粒浊积岩沉积规律——以松南西斜坡大布苏地区青一段地层为例 . 沉积学报, 2013, 31(2): 291-301.
    [13]  基于滨线轨迹的古水深定量计算新方法①——以古近系沙三中段东营三角洲为例 . 沉积学报, 2012, 30(3): 443-450.
    [14] 梁宏斌.  断陷湖盆缓坡带高位三角洲体系与油气成藏组合特征分析——以二连盆地吉尔嘎朗图凹陷为例 . 沉积学报, 2011, 29(4): 783-792.
    [15] 唐 勇.  准噶尔盆地腹部缓坡型岩性地层油气藏成藏控制因素分析 . 沉积学报, 2009, 27(3): 567-573.
    [16] 辛仁臣.  湖盆无曝露缓坡带层序界面特征及成因——以松辽盆地他拉哈地区上白垩统为例 . 沉积学报, 2008, 26(1): 77-85.
    [17] 王颖, 王英民, 王晓洲, 辛仁臣, 赵志魁, 赵占银.  松辽盆地西部坡折带的成因演化及其对地层分布模式的控制作用 . 沉积学报, 2005, 23(3): 498-506.
    [18] 胡宗全, 朱筱敏.  具有地形坡折带的坳陷湖盆层序地层模拟 . 沉积学报, 2002, 20(2): 217-221.
    [19] 李凤杰, 王多云, 郑希民, 刘自亮, 王峰, 李树同, 王志坤.  陕甘宁盆地华池地区延长组缓坡带三角洲前缘的微相构成 . 沉积学报, 2002, 20(4): 582-587.
    [20] 刘震, 曾宪斌, 张万选.  构造掀斜对单断湖盆湖平面变化的影响 . 沉积学报, 1997, 15(4): 64-66,71.
  • 加载中
图(10) / 表 (1)
计量
  • 文章访问数:  186
  • HTML全文浏览量:  50
  • PDF下载量:  103
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-10-08
  • 修回日期:  2022-05-09
  • 刊出日期:  2022-12-10

目录

    构造掀斜主导的断陷湖盆缓坡层序“源—汇”正演模拟定量研究

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

      江苏省自然科学基金 BK20200529

      中央高校基本业务费 B220202047

    • 中图分类号: P618.13

    摘要: 断陷湖盆缓坡带的构造掀斜作用对其地层叠置样式和沉积体系展布具有重要的控制作用。通过基于“源—汇”系统理论的地层正演模拟研究方法,在三级层序地层格架下(以2 Myr为总时长)对断陷盆地缓坡带构造掀斜型层序地层发育样式及其相关地貌演化、沉积物通量变化进行了定量表征。以恒定的绝对湖平面为研究前提,共建立了高/低沉积物供给的背景下4种不同的地层正演模型。两种不同缓坡带掀斜初始位移模型的设定反映了掀斜点和沉积滨线的相对位置,以便于不同掀斜层序样式对比分析。研究表明,掀斜支点位于滨线之上的模型常形成湖侵和正常湖退,而掀斜支点位于滨线之下的模型常形成强制湖退。滨线迁移轨迹和层序叠置样式的不一致性贯穿整个层序发育过程,是沉积物通量及滨线迁移共同作用的结果。该模型中沉积物通量的变化直接响应于源区地貌演变,滨线迁移则受控于新增可容空间和沉积物供给的相互关系。这一系列构造掀斜型层序模型及其相关成因解释对断陷盆地层序地层格架及其叠置样式的精细解剖具有重要的指导意义。该模拟研究对断陷湖盆的高频湖平面变化特征也有较为合理的解释(如海洋同位素3期以来的贝加尔湖)。然而本次正演模拟所设置的初始地貌较为单一,在相同构造、气候条件下和不同沉积物供给背景下所形成地层模型无法涵盖所有缓坡带层序地层类型。由此可见,基于“源—汇”系统思想探讨断陷盆地层序地层多样性的相关模拟工作仍然任重道远。

    English Abstract

    谈明轩, 朱筱敏, 张自力, 孙浩南. 构造掀斜主导的断陷湖盆缓坡层序“源—汇”正演模拟定量研究[J]. 沉积学报, 2022, 40(6): 1481-1493. doi: 10.14027/j.issn.1000-0550.2022.061
    引用本文: 谈明轩, 朱筱敏, 张自力, 孙浩南. 构造掀斜主导的断陷湖盆缓坡层序“源—汇”正演模拟定量研究[J]. 沉积学报, 2022, 40(6): 1481-1493. doi: 10.14027/j.issn.1000-0550.2022.061
    TAN MingXuan, ZHU XiaoMin, ZHANG ZiLi, SUN HaoNan. Source-to-sink Quantitative Stratigraphic Forward Modeling on the Tilted Hanging-wall Sequence Architecture of a Tectonically-driven Lacustrine Rift Basin[J]. Acta Sedimentologica Sinica, 2022, 40(6): 1481-1493. doi: 10.14027/j.issn.1000-0550.2022.061
    Citation: TAN MingXuan, ZHU XiaoMin, ZHANG ZiLi, SUN HaoNan. Source-to-sink Quantitative Stratigraphic Forward Modeling on the Tilted Hanging-wall Sequence Architecture of a Tectonically-driven Lacustrine Rift Basin[J]. Acta Sedimentologica Sinica, 2022, 40(6): 1481-1493. doi: 10.14027/j.issn.1000-0550.2022.061
      • 构造掀斜作用是断陷盆地构造沉降和隆升的重要机制之一[1]。掀斜作用通常包括正向掀斜(向陆一侧相对抬升、向盆一侧相对沉降)和反向掀斜(向陆一侧相对沉降、向盆一侧相对抬升)两类[12]。其中,正向掀斜作用主要发育于断陷盆地缓坡带上盘,也是本文的主要研究对象。历史地震监测数据证实现代活动断陷盆地存在短期时间尺度的幕式掀斜作用[34]。构造掀斜作用不仅控制了断陷盆地的水系样式和剥蚀作用[56],而且对湖平面变化具有极为重要的影响[1,78]。前人理论研究表明,在构造主控的断陷湖盆中,气候变化对湖平面影响相对微弱,湖平面变化与构造掀斜枢纽点的相对位置密切相关[9]。枢纽点位于湖平面之上造成缓坡带相对湖平面上升;枢纽点位于湖平面之下时造成缓坡带相对湖平面下降(图1)。然而目前大多数研究停留于几何学模型和地层定性分析层面[7,10],缺乏相关实例的验证,其理论可靠性仍有待进一步检验。

        图  1  构造掀斜作用下恒定绝对湖平面的缓坡带滨线迁移模式(据文献[9]修改)

        Figure 1.  Schematic model showing shoreline migration influenced by tilting with a fixed absolute lake level in the hanging⁃wall dipslope (modified from reference [9])

        可容空间和沉积物供给是断陷盆地充填特征及其地层建造样式的重要控制因素[5,11]。由于构造掀斜作用能够直接影响源区剥蚀量和盆地沉降量,因此是造成掀斜型缓坡层序多样性的根本原因之一。然而断陷盆地缓坡带在掀斜作用下遭受广泛的暴露剥蚀,地层发育通常并不完整,为缓坡带层序地层学研究带来多解性。地层数值模拟的优势在于其能够通过改变边界条件建立不同类型层序地层样式,能够正演地层沉积充填过程,是定量层序地层学分析的一种有效手段。地貌演化与“源—汇”系统近年来是地质学热点议题,是构造地貌学、盆地动力学与沉积学交叉研究内容[1213]。本次数值模拟研究主要基于“源—汇”沉积通量守恒原理的三维地层正演模拟方法[1416],分析构造掀斜作用相关的剥蚀与沉积作用及其对断陷湖盆缓坡带湖平面变化、层序地层充填的直接影响,以“源—汇”系统研究手段探讨掀斜型缓坡层序样式的主要驱动机制。

      • 本次模拟所应用的Badlands软件是Python架构下以河流作用和坡面过程模拟侵蚀作用、泥沙运移和沉积过程的地貌演化开源程序[1720]。与其他基于扩散模型、几何模型和水力学模型的地层正演模拟软件相比[2122],该地层正演模拟代码充分考虑了泥沙沉积作用与流域侵蚀作用的内在关系[18]。因此,该代码所模拟的层序地层发育样式和沉积演化特征不仅能够直接响应于构造掀斜作用,而且能将沉积物通量和流域地貌演化直接纳入“源—汇”系统研究。

      • 该正演地层模型主要基于受“分离”限制(“D-L”模型)的水力侵蚀模型定量表征河流下切作用。岩石剥蚀速率E是源区面积(A)、净降雨量(P)、坡度(S)的函数,即:

        E= εAPmSn (1)

        式中:剥蚀系数ε(yr-1)是气候、底床糙度及岩石强度的函数。无量纲参数mn默认值分别为0.5和0.1[1718]。由此可见,沉积物供给主要受控于源区的汇水面积、地形高差及山体陡度等地貌参数。

        此外,该模型通过利用线性或非线性扩散模拟坡面过程[19],即:

                  zt= κ2z (2)

        式中:z为地形高差,扩散系数κ(m2/yr)与源区岩性及降雨量有关[2324]。地表、水下及河流的扩散系数均有所差异。其中,冲积平原区临界沉积坡度为0.001 m/m。当沉积物进入水下,则沿多方向水流进行扩散搬运。

        此外,Badlands软件后处理工具包能够自动提取正演模型中不同时期滨线位置的新增可容空间深度δA(m)和新增沉积地层厚度δS(m)[17],为本次定量模拟地层学研究提供了新的方法。

      • 本研究将断陷盆地缓坡带简化为平缓的单斜地形,其初始地形设置为100 km(长度)×100 km(宽度)、地表网格间距设置为1 km。该单斜初始地形最大相对高差为200 m,最大水深为-400 m,平均坡降为0.006 m/m(图2)。为更直观地揭示地层垂向叠置特征,该正演数值模型的运行时间设置为2 Myr(三级层序时间尺度),单层沉积时间为10 kyr(准层序时间尺度),累计沉积了200层地层单元。

        图  2  简化的缓坡带初始三维古地形

        Figure 2.  Simplified initial three⁃dimensional paleo⁃topography of the hanging⁃wall dipslope

      • 为突出构造掀斜对层序构型的控制作用,本次模拟研究所设置的模型均简化为恒定的绝对湖平面及模拟研究区均一的降雨量(即2 000 mm/yr),以达到控制单一变量的目的。在湖盆中,绝对湖平面是由相对与一个固定基准面(例如地心)所定义的,相对湖平面则是湖盆底面距湖面的高度[25]。绝对湖平面在自然界实例较少,但也并不是完全不可能实现的。闭流湖盆可通过降雨量和蒸发量的平衡所实现,而敞流湖盆则可通过构造稳定的出水口保持恒定的绝对湖平面。在本研究简化的构造模型中,构造掀斜的枢纽位置固定不变,位移量与距离呈线性关系。基于枢纽点和绝对湖平面的相对位置,本研究共建立两个构造模型:构造模型A(最大隆升量和沉降量分别为100 m和-400 m)(图3a);构造模型B(最大隆升量和沉降量分别为300 m和-300 m)(图3b)。在模拟运行过程中,研究区的抬升量和沉降量呈线性的变化特征。

        图  3  两种掀斜构造模型的剖面示意及三维立体特征

        Figure 3.  Schematic profiles of two tectonic tilting models and three⁃dimensional characteristics

        为系统探讨层序地层样式多样性可能,本次模拟研究还额外考虑沉积物供给对掀斜缓坡带地层的影响,因此共建立了4种地层正演模型(表1)。其中,模型1、2和模型3、4分别为高、低沉积物供给背景,主要通过源区岩性差异所造成的高、低风化强度进行调节。本研究主要参照希腊Corinth湾裂谷盆地母岩区多类型基岩剥蚀能力的相关研究[6],将剥蚀系数ε值设置为1×10-6 yr-1代表源区主体以硬岩性组成,岩石扛风化能力较强,同等气候条件下沉积物供给能力较弱;ε值设置为1×10-5 yr-1代表源区为混合岩性组成,岩石扛风化能力较弱,同等气候条件下沉积物供给能力较强[2627]。其他一系列(地表、水下和河流)扩散系数的取值主要参考了前人三角洲数值模拟研究中所选取的数量级,并在各模型设置中保持不变(表1),便于在数值模拟中通过控制相关变量进行对照分析。通过公式(1)进行运算,模型3、4所产生的总沉积物通量仅为模型1、2计算结果的10%左右。由于该正演模拟程序的水下扩散模型较为简单,并没有分别考虑砂泥扩散特征,因此多种沉积相类型主要依据水深进行划分[1718,28]

        表 1  地层正演模拟关键参数基本设置

        模型构造模型侵蚀系数/(yr-1地表扩散系数/(m2/yr)水下扩散系数/(m2/yr)河流扩散系数/(m2/yr)
        1A1×10-50.300.505.00
        2B1×10-50.300.505.00
        3A1×10-60.300.505.00
        4B1×10-60.300.505.00
      • 模型1中,湖岸线向盆地边缘后退了12 km,表明掀斜作用造成缓坡带的相对湖平面上升(图4a)。湖岸线在初始0.75 Myr向盆缘的平均迁移速率为10.7 km/Myr,而后其平均迁移速率降至4.8 km/Myr(图4b)。该模型中源区地貌最终呈夷平特征(图4c),这表明虽然掀斜作用可造成源区不断抬升,但强烈剥蚀的能力同样可导致源区快速夷平。尽管相对湖平面在不断上升,但该模型中的地层发育样式类似于一个完整的被动陆缘I型三级层序构型特征(图4b,c)。虽然湖平面并没有经历经典层序的海平面变化特征,该层序早期和晚期发育的进积型三角洲与其层序模式中的低位前积楔和高位期进积三角洲相似程度较高(图4b)。

        图  4  模型1地层剖面、平面模型及相关沉积物通量变化特征

        Figure 4.  Stratigraphic profile and planar model of model run 1 and variation of sediment flux

        本层序模型中类“低位体系域”(LST)持续时间约为0.3 Myr(图4c),主要发育进积型地层叠置样式。此时期研究区总沉积物通量达到峰值(约为600 km3/Myr)(图4d),三角洲进积型叠置样式与滨线迁移方向相反,整体符合低位楔的发育特征。类“湖侵体系域”(TST)持续时间达0.45 Myr(图4c),然而此时期研究区总沉积物通量骤减为370 km3/Myr(图4d)。类“高位体系域”(HST)则主要表现为加积型地层叠置样式,该时期研究区总沉积物通量降至约200 km3/Myr(图4d)。

      • 模型2中,湖岸线向盆地方向推进了18 km,表明掀斜作用造成缓坡带相对湖平面下降(图5a)。湖岸线向盆地方向的平均迁移速率高达45 km/Myr(发育于初始0.4 Myr),同时伴随着广泛的陆上暴露剥蚀(图5b),随后1.6 Myr湖岸线迁移基本停止。该模型中源区地貌最终呈夷平的准平原化特征(图5c),表明源区基本已剥蚀殆尽。由于研究区发育较为广泛的湖退作用,因此该层序中地层样式主要发育类“下降期体系域”(FSST)和类“低位体系域”(LST)(图5b)。

        图  5  模型2地层剖面、平面模型及相关沉积物通量变化特征

        Figure 5.  Stratigraphic profile and planar model of model run 2 and variation of sediment flux

        本层序模型中类“下降期体系域”(FSST)主要发育于初始0.4 Myr,强制海退现象明显,造成了显著的陆上暴露剥蚀,甚至形成类似于层序界面的特征(图5a)。该时期滨线迁移轨迹为下降型,同时研究区总沉积物通量在短时间内达到峰值(不低于470.4 km3/Myr)。类“低位体系域”(LST)上超于类“下降期体系域”(FSST)之上,事实上已经形成新的层序界面(图5a)。由此表明,掀斜作用是构造主控型断陷湖盆层序界面成因之一[29]。由于相对湖平面持续下降,滨岸平原也造成持续遭受剥蚀(图5b)。该体系域主要表现为降积、加积—进积型的地层叠置样式,同期该模拟研究区的总沉积物通量进一步降低(图5c,d)。

      • 尽管模型3中的相对湖平面不断上升,缓坡带三角洲却整体保持相对稳定(图6a,b)。在该模型中,湖岸线向盆地边缘后退了11 km,其平均迁移速率达5.5 km/Myr(图6b)。该模型中源区地貌以线性和树枝型沟谷发育为特征(图6c),表明源区仍处于强烈剥蚀阶段。尽管使用了相同的构造模型(表1),低沉积物供给的模型3与高沉积物供给下的模型1在地层样式上已存在较大的差异(图4a、图6a)。该地层正演模型并不具有经典层序三分构型特征,垂向上主要表现为退积—加积型地层样式,仅与高位体系域(HST)更为类似。与高沉积供给条件的模型1对比,该模型中研究区的总沉积物通量整体不超过60 km3/Myr(仅为模型1总沉积通量的近十分之一),其变化趋势却全然不同(图6d)。模型3中模拟研究区的沉积物通量变化在初始阶段达到峰值56.7 km3/Myr,而后缓慢降至55.6 km3/Myr,并最终缓慢上升至56.6 km3/Myr(图6d)。

        图  6  模型3地层剖面、平面模型及相关沉积物通量变化特征

        Figure 6.  Stratigraphic profile and planar model of model run 3 and variation of sediment flux

      • 模型4中源区地貌线性和树枝型沟谷发育(图7c),表明源区仍处于强烈剥蚀阶段。该模型中湖岸线最终向盆地方向推进了模型2中同样的距离,表明相对湖平面下降的幅度是一致的(图5b、图7b)。然而其湖岸线迁移特征却迥然相异,整体具有长期、缓慢的特征。模型4在1.6 Myr中平均迁移速率为11.25 km/Myr,并在最后0.4 Myr迁移速率逐渐近乎停滞(图7b)。

        图  7  模型4地层剖面、平面模型及相关沉积物通量变化特征

        Figure 7.  Stratigraphic profile and planar model of model run 4 and variation of sediment flux

        该模型的地层样式整体上与模型2类似,主要包括类“下降期体系域”(FSST)和类“低位体系域”(LST)(图7a,b),但其总沉积物通量变化特征具有较大的差异(图8c)。类“下降期体系域”(FSST)发育时间为初始的1.6 Myr,主要表现为下降型迁移轨迹和显著陆上剥蚀现象(图7b)。该时期研究区沉积物通量先快速降至52 km3/Myr,并缓慢上升至72 km3/Myr(图7d)。类“低位体系域”(LST)仅在最后0.4 Myr发育,主要发育进积型地层叠置样式。同时期模拟研究区的总沉积物通量逐渐上升,并最终达到峰值75 km3/Myr(图7d)。

        图  8  不同模型中滨岸带新增可容空间(δA)与沉积物供给(δS)变化特征

        Figure 8.  Variation of newly⁃created accommodation space (δA) and sediment supply (δS) in the coastal area

      • 虽然掀斜作用对断陷盆地缓坡带湖平面变化的影响在先前相关研究中均有所涉及[89],但是本次基于简化模型的模拟研究表明其变化特征比概念模式更为复杂。即使在相同构造模型边界条件下,仍然可能具有差异性的滨线迁移特征。模式1和模型3中相对湖平面呈上升趋势,δA为正值。正常情况下,δA与δS差值也为正值(图8a,c),表现为整体湖侵的特征(图9a)。然而实际上δA与δS差值却并不总是正值(图8a,c),表明部分时期沉积滨线处新增沉积物沉积厚度大于新增可容空间的深度。在这种情况下(即δA-δS<0),即使相对湖平面不断上升,但由于三角洲沉积速率相对较快,总可容空间也会不断降低[22],也会产生正常湖退的特征(图9b)。

        图  9  构造掀斜作用下缓坡带湖平面变化新模式

        Figure 9.  New schematic model of lake level fluctuation of the hanging⁃wall dipslope driven by tectonic tilting

        模型2、4整体呈相对湖平面下降的趋势,因而造成强制湖退现象(图9c),因此δA为负值。δA与δS差值基本均为负值,且δA与δS差值逐渐趋近于零(图8b,d)。这与两个模型滨线迁移速率降低密切相关,湖岸线开始趋于稳定,δA逐渐趋近于零值。由此可见,在基于构造模型B建立的模型2、4中,其滨线迁移轨迹主要受控于可容空间变化速率。

      • 与经典被动陆缘层序相比,掀斜型缓坡层序样式的复杂性在于其滨线迁移轨迹与地层叠置样式的不一致性。模型1、3的滨线迁移轨迹属于向陆迁移型,但是地层地层叠置样式并不以退积型为主导;模型2、4的滨线迁移轨迹属于下降型,但其地层叠置样式亦可发育加积型特征。由此表明,迁移轨迹分析方法对于其层序构型的分析效果并不理想。这种滨线迁移和地层叠置样式不一致性现象的存在与沉积物供给的主导作用密切相关。前人认为掀斜作用会使缓坡带坡度增大,增强了沉积物供给能力,并最终造成了整体湖侵趋势下三角洲进积现象[30]。断块翘倾旋转作用使得苏伊士湾海相断陷盆地中新世碳酸盐沉积向混积岩转化,证实这种构造作用具有增强陆源碎屑沉积物供给能力的作用[31]

        构造掀斜作用并不一定造成沉积物供给能力的增强。本模型所设置的初始地貌为低地丘陵(最大地形高差仅为200 m),在高沉积物供给条件下,其地形高差和陡坡会迅速减小,源区逐渐呈夷平化地貌特征。该背景下掀斜作用所能增强的沉积物通量是极为有限的。模型1分析表明,在整体湖侵的背景下,源区的汇水面积不断缩小,沉积物供给能力并不一定提高。与前人认识不同的是,沉积物供给速率在初始0.4 Myr就可能达到峰值,有利于三角洲向盆地方向进积,因而在初始阶段就存在这种不一致性(图4a,b)。随着沉积物供给速率的降低,这种层序样式与迁移轨迹的不一致性逐渐消失。在整体湖退趋势下,模型2沉积物速率逐渐减少,使得后期三角洲沉积作用减弱,因而主要表现为加积样式,而非进积样式(图5a,b)。

        在低沉积物供给条件下,源区风化剥蚀速率较低,模型3、4的源区水系在2 Myr之后仍处于壮年期(图6c、图7c)。该背景下构造掀斜作用则会进一步增强研究区的总沉积物通量(图6d、图7d)。模型3的汇水面积不断缩小,但沉积供给能力却略微上升,表明掀斜作用贡献程度相对更大。其地层叠置样式却并不整体表现为进积特征,也证实沉积物供给速率增幅较为有限。在模型4中,在旋转掀斜和强制湖退的共同作用下,沉积物供给能力快速增强,使得三角洲进积速率增大,因而在研究区与沉积滨线的迁移轨迹保持一致性。

        综上所述,在两种不同沉积物供给条件下,沉积物通量大小和变化特征均有很大差异,是造成地层叠置样式和相关“体系域”类型差异的一个关键因素。

      • 本次正演模拟研究表明断陷盆地层序具有复杂性和多样性,并不是被动陆缘层序的简单套用。虽然通常认为构造掀斜作用直接影响可容空间的变化,但实质上是其枢纽点与湖平面相对位置进而决定了可容空间的增减。响应于源区地貌的沉积物通量变化对地层垂向分布特征具有极为重要的控制作用,尤其是在沉积物供给主导情况下形成一系列与湖平面变化相悖的地层样式(如相对湖平面上升背景下发育类“体系域”三分特征,见图4;相对湖平面下降背景下出现地层超覆特征,见图5)。因此针对这种掀斜型缓坡带层序,简单的体系域划分实际上是较为困难的,需要结合地层叠置样式和滨线迁移轨迹共同系统解释。由于本模拟研究设置的初始地形仅为相对单一的低地丘陵,源区总剥蚀产物是相对有限的。因此,本模拟研究仅代表其中一种或几种缓坡带层序可能性的发育模式,对断陷盆地内部掀斜断块形成的局部物源体系(如凹陷内的低凸起型物源体系)更有参考价值。通过多种地貌初始条件,进一步细化构造掀斜模型(甚至考虑反向掀斜、枢纽点迁移以及幕式构造运动等多因素)及气候模型(如季风影响、雨影效应、地形效应等),能够建立更完善、更具多样性的缓坡带层序,以发展和丰富断陷湖盆层序地层学理论。

        由于本研究所设定的绝对湖平面边界条件十分理想,现实中很难找到实例进行佐证[1]。地处高纬度的贝加尔湖在晚更新世以来可能符合类似的条件(图10a,b)[32],但其发育的高频层序(相当于六级层序)却远低于本次正演模拟的研究精度。然而如果将构造和沉积间隔细化到千年尺度,取得的精细模拟结果至少能够反映枢纽点位置与相对湖平面变化的关系。前人研究表明,贝加尔湖自海洋同位素阶段3期(MIS-3期)以来基本保持绝对湖平面长期稳定的状态(图10c)[3334]。由于湖平面下降超过2 m就使其成为闭流湖盆,因此该时期气候驱动的湖平面下降持续时间并不长[35]。其缓坡带Selenga三角洲发育区地震剖面所表现的地层叠置样式及局部相对湖平面变化曲线表明该区域相对湖平面变化频繁(图10d)[36]。基于前人研究成果,这种整体后退型的滨线迁移轨迹表明研究区的相对湖平面自MIS-3期以来整体呈逐渐上升的趋势(图10e)。虽然气候型层序观点(气候主导的相对湖平面变化)能够解释这种滨线迁移轨迹发育特征[36],然而这种高纬度湖泊对于气候变化的敏感度并不高,本质上属于构造主控型湖盆[32]。在绝对湖平面恒定的条件下,贝加尔湖缓坡带Selenga三角洲自MIS-3期以来发育的地震地层单元IIa、Ib及Ia的垂向叠置样式可用掀斜型缓坡层序模式进行成因解释(图10d)。前人地震地层学研究亦表明,Selenga三角洲发育区具有典型的构造掀斜特征,并伴随显著的湖侵作用[32]。基于本次正演模拟的相关研究成果,其构造掀斜枢纽点推测位于湖平面之上。地震地层单元Ib降积样式的发育是海洋同位素阶段2/3期(MIS-2/3)界限处全球气候恶化所致[37]。该时期气候变化造成贝加尔湖的湖平面在短期内可能迅速下降数十米,并形成闭流湖盆[3738]。在构造掀斜作用下,该缓坡带的相对湖平面在新仙女木事件之后表现为持续上升的趋势(图10e)。

        图  10  (a,b)贝加尔湖Selenga三角洲地理位置;(c)中—晚更新世以来贝加尔湖(相对及绝对)湖平面变化特征;(d)基于地震解释成果的地层剖面E—E’;(e)贝加尔湖Selenga三角洲发育区NW区块及SE区块相对湖平面变化特征(修改自文献[33,36])

        Figure 10.  (a, b) Location of Selenga delta in Lake Baikal; (c) (relative and absolute) lake level fluctuation of Lake Baikal since the Middle to Upper Pleistocene; (d) stratigraphic profile E⁃E’ based on seismic interpretation; and (e) relative lake levels of NW and SE blocks of the Selenga delta in Lake Baikal (modified from references [33,36])

      • 断陷盆地缓坡带的构造掀斜作用使得源区隆升遭受风化剥蚀,汇区沉降形成可容空间,直接影响了其地层发育样式。通过解剖本次正演数值模拟所建立的4个定量层序地层模型,主要得出以下结论。

        (1) 掀斜型缓坡层序具有复杂性,其相对湖平面变化与实际地层叠置样式类型并不具有一致性,所发育的与传统体系域类似地层单元并不能与实际沉积滨线的迁移轨迹匹配。

        (2) 掀斜型缓坡的滨线迁移轨迹整体受控于沉积滨线处可容空间和沉积速率变化。枢纽点位于湖平面之上形成湖侵和正常湖退两种类型,受可容空间和沉积物供给两种因素共同影响;枢纽点位于湖平面之下形成强制湖退,仅受可容空间影响。

        (3) 地层叠置样式的发育与沉积物供给动态变化、滨线迁移轨迹密切相关。沉积物供给直接响应于源区地貌演化,因此基于“源—汇”系统的地层正演模拟对于探讨不同地貌条件下掀斜型缓坡层序多样性具有重要的应用价值。

    参考文献 (38)

    目录

      /

      返回文章
      返回