地球科学进展, 2020, 35(9): 890-901 DOI: 10.11867/j.issn.1001-8166.2020.072

综述与评述

海底周期阶坎研究进展与挑战

王大伟,1,2,3, 孙悦1,2, 司少文1,2, 吴时国1,2,3,4

1.中国科学院深海科学与工程研究所,海南 三亚 572000

2.中国科学院大学,北京 100049

3.南方海洋科学与工程广东省实验室(珠海),广东 珠海 519080

4.青岛海洋科学 与技术国家实验室,山东 青岛 266237

Research Progress and Challenges of Submarine Cyclic Steps

Wang Dawei,1,2,3, Sun Yue1,2, Si Shaowen1,2, Wu Shiguo1,2,3,4

1.Institute of Deep-sea Science and Engineering,Chinese Academy of Sciences,Sanya Hainan 572000,China

2.University of Chinese Academy of Sciences,Beijing 100049,China

3.Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai),Zhuhai Guangdong 519080,China

4.Qingdao National Laboratory for Marine Science and Technology,Qingdao Shandong 266237,China

收稿日期: 2020-08-02   修回日期: 2020-08-25   网络出版日期: 2020-10-27

基金资助: 国家自然科学基金项目“南海珠江口外海底峡谷内底形沉积结构与形成机理”.  41666002
中国科学院前沿科学研究重点计划项目“深水海底峡谷—水道地貌特征及形成机理”.  QYZDB-SSW-SYS025

Received: 2020-08-02   Revised: 2020-08-25   Online: 2020-10-27

作者简介 About authors

王大伟(1976-),男,黑龙江绥化人,研究员,主要从事地震沉积学、深水油气、海洋地质灾害和人工智能识别研究.E-mail:wangdawei@idsse.ac.cn

WangDawei(1976-),male,SuihuaCity,HeilongjiangProvince,Professor.Researchareasincludeseismicsedimentology,deepwateroilandgas,submarinegeohazard,andartificialintelligenceidentification.E-mail:wangdawei@idsse.ac.cn

摘要

超临界流形成的周期阶坎是一种常见的底形,普遍存在于河口三角洲体系、深水峡谷—水道体系和碳酸盐台地等沉积体系中。从沉积环境、底形形态、沉积结构、形成机理和数值模拟等方面介绍了周期阶坎的研究进展,探讨了船测和AUV多波束、船测和AUV浅剖、人工地震、野外露头、数值模拟研究中揭示的周期阶坎分辨率问题,给出了周期阶坎研究的突破方向。利用探地雷达、人工地震、钻井和测井结合的方法,准确确定野外露头的中泓线位置;在水深大于500 m的区域,将AUV和船测多波束、AUV和船测浅地层剖面、人工地震等数据结合,通过载人潜水器获取原位数据,完善周期阶坎的三维精细结构;开发立体化三维数值模拟,建立一套适用于真实周期阶坎形态的动力学模拟方程。

关键词: 周期阶坎 ; 浊流 ; 超临界流 ; 沉积结构 ; 尺度差异

Abstract

Cyclic steps structure derived from the supercritical flows is one of the common bedforms, which is commonly found in sedimentary systems such as delta systems, deep-water canyon-channel systems, and carbonate platforms. In this paper, the research progresses of cyclic steps were introduced from the aspects of sedimentary environment, the features of bedforms, sedimentary structure, formation mechanism, and numerical simulation. The Research scale differences associated with the survey techniques or study methods, such as the ship survey and AUV-based multi-beams, sub-bottom profiles, multichannel seismic, the field outcrop, and numerical simulation were discussed. Finally, the breakthrough directions of the cyclic steps research were given. The combining method of ground penetrating radar, multichannel seismic, drilling and well logging were used to accurately detect the location of field outcrop thalweg. Within the area with water depth greater than 500 meters, the blending of multi-data for the cyclic steps research involved the AUV-based and the ship multi-beams, sub-bottom profiles, and multichannel seismic data. With the in situ samples and observation data obtained by human occupied vehicles, three-dimensional numerical simulation was developed to establish a set of dynamic simulation equations suitable for the real cyclic steps. Therefore, the high resolution three dimensional mode of the deep-water cyclic steps could be obtained more accurately.

Keywords: Cyclic steps ; Turbidity currents ; Supercritical flows ; Sedimentary structure ; Scale difference

PDF (8391KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

王大伟, 孙悦, 司少文, 吴时国. 海底周期阶坎研究进展与挑战. 地球科学进展[J], 2020, 35(9): 890-901 DOI:10.11867/j.issn.1001-8166.2020.072

Wang Dawei, Sun Yue, Si Shaowen, Wu Shiguo. Research Progress and Challenges of Submarine Cyclic Steps. Advances in Earth Science[J], 2020, 35(9): 890-901 DOI:10.11867/j.issn.1001-8166.2020.072

近些年,超临界流逐渐引起人们的关注,在海底上形成一系列的沉积波。在沉积波中,有一些呈向上游迁移的、长波状(波长/波高大于1)、一系列连续的陡峭背流面和平缓迎流面组合,并且经常在背流面与迎流面转换区域形成冲沟或深坑。这种沉积波代表流动状态从超临界流到亚临界流转换的一个过程,称之为“周期阶坎(cyclic steps)1图1)(① 同济大学钟广法教授在南海深部计划2013年度学术研讨会上首次提出“周期阶坎”这一中文译名)。在每个完整的周期阶坎中,以水跃为分割点,完成超临界流到亚临界流的转换,形成背流面到迎流面的底形。在波长超过1 km的周期阶坎中,由于迎流面的超临界流具有高能量、高流速和高浓度的特点,在迎流面容易形成物性较好的砂质油气储层,可以指导精细油气储层预测。在现今海底形成的周期阶坎,由于具有高速的沉积物搬运速率,对海底基础设施构成一定威胁。超临界流对海底光缆、海洋油气开发等海洋工程构成了潜在威胁3~5

图1

图1   周期阶坎示意图[2]

Fig.1   Schematic diagram of cyclic steps[2]


在海底峡谷、水道、深水沉积扇、河口三角洲和碳酸盐台地斜坡等沉积环境中均发育超临界流,控制流经区域的海底底形特征6~9,对于“源—汇”系统研究具有指导意义。通过多波束、浅剖、人工地震和海水测流数据的分析,揭示不同沉积环境中的周期阶坎是由不同类型超临界流形成。浊流形成的周期阶坎主要存在于海底峡谷、水道和河口三角洲10。例如,南海西澎湖峡谷中的周期阶坎,波长为1.2~10.0 m,波长/波高为23~250,长达一百多公里,是目前规模最大的周期阶坎1112。倾泻密度流(cascading density currents)形成的周期阶坎主要存在于碳酸盐台地周缘,例如,墨西哥湾巴哈马台地边缘。倾泻密度流是一种类似于海底瀑布的线性源超临界流,流动过程中发生水跃,形成周期阶坎1314。周期阶坎沉积结构的分析主要通过野外露头描述15~19,例如,Postma等19通过露头和岩心描述了Tabernas盆地周期阶坎的层状构造;Lang等15描述了德国北部古冰川野外露头中的浊积岩体系。超临界底形的发育演化过程主要通过水槽实验进行物理模拟和计算机进行数值模拟。例如,Taki等20和Cartigny等21对超临界流的流动特性以及不同超临界底形之间的形态动力学关系进行了阐述。Kostic22通过使用“平均四方程”建立了周期阶坎的数值演化模型,Vellinga等23利用FLOW-3D模拟了水槽实验的超临界底形发育的动态过程。在理论方面,利用泥沙运动力学理论建立了控制方程24~26,并且,总结了周期阶坎的研究现状与发展趋势1227

然而,随着超临界流和水跃研究的深入,逐渐发现现有的周期阶坎研究手段存在一些不足。本文将就以下相关不足进行讨论:不同观测手段对深水区周期阶坎等超临界底形的观测结果存在一定的偏差;由于野外露头的空间局限性,周期阶坎露头并不能刻画真实的古超临界流路径;短期的数值模拟可能与实际的多次、长期浊流冲刷相驳。在此基础上,对周期阶坎的未来研究方向做出展望。

1 超临界流底形

超临界流(supercritical flow)是指弗洛德数大于1的流体,弗洛德数是流体力学中表征流体惯性力和重力相对大小的一个无量纲参数,记为Fr。它表示惯性力和重力量级的比,即:Fr=U2/gh,式中U为流体运动速度,g为重力加速度,h为流体厚度。超临界流一般产生于坡度大于0.45°的斜坡上,深刻影响着深水沉积体系的形成与演化28~30

超临界流通过水跃转变为亚临界流的过程中形成各种底形31,底形的特征取决于水流条件和泥沙颗粒大小32,可以形成逆行沙丘(antidune)、流槽—冲坑(chute-and-pool)和周期阶坎等底形21。逆行沙丘呈1~2 m长的透镜状沉积,根据稳态和非稳态特征,内部层理分为水平和交错层理,倾向上游,倾角较小,常与平行层理及冲刷构造共存,弗洛德数为0.8~1.833。流槽—冲坑沉积具有陡峭的、向上倾斜的层状纹理,流体的弗洛德数(1.6~2.1)高于逆行沙丘2134,呈现不规则间隔、阶梯式水跃迁移的特征15。周期阶坎很难与逆行沙丘区分开来,两者的形成条件和沉积过程极其相似15。下面列举了周期阶坎与逆行沙丘的异同点35图2)。

图2

图2   逆行沙丘和周期阶坎对比(据参考文献[35]修改)

Fig.2   The comparison between antidunes and cyclic steps (modified after reference [35])


2 周期阶坎的形成环境与底形特征

周期阶坎的形成主要与海底地形和超临界流过程有关。一方面,海底地形影响了周期阶坎的发育,形成具有不同特征的周期阶坎;另一方面,周期阶坎是海底底形的一种具体形式,周期阶坎的特征参数反映了海底地形地貌特征和流体状态。

2.1 沉积环境

在高山、河流、河口三角洲、海底,甚至在火星36上均发现了周期阶坎的存在,不同环境的地形条件、水动力条件和底质条件等因素都控制着周期阶坎的形态特征。海底周期阶坎主要形成于河口三角洲体系、海底峡谷—水道体系和碳酸盐台地中(表1)。

表1   不同沉积环境下的周期阶坎

Table 1  Cycle steps in different sedimentary environments

沉积环境尺度大小形成机理结构特征沉积物组成实例
河口三角洲

波长:30~300 m

波高:1~8 m

多数由羽状流触发的浊流多为新月形,阶坎较为连续,数量较多砂质、泥质沉积物加拿大斯阔米什三角洲[37]
海底峡谷水道

小型:整体长度小于10 km,波长20~100 m,波高0.5~5.0 m;

中型:整体长度10~50 km,波长200~1 100 m,波高20~50 m;

大型:整体长度大于50 km,波长300~7 200 m,波高5~220 m

海底峡谷—水道内壁垮塌、驻波、内波等触发的浊流以新月形为主,上游尺度较小、对称性较差,下游尺度较大、对称性较好近端粗粒砂岩,远端细粒砂岩、黏土南海东北陆坡西澎湖峡谷[12],非洲Rio Muni盆地[11]
碳酸盐台地

波长:96~117 m

波高:3.5~4.3 m

台地顶部沉积物向深海盆地搬运形成的线性倾泻密度流不仅沿流动方向纵向排列,横向上也成线性排列碳酸盐岩墨西哥湾巴哈马台地[13,14,38]

新窗口打开| 下载CSV


2.2 形态特征

通过高分辨率多波束数据分析,揭示峡谷水道中的周期阶坎外形特征呈现为排成一列的新月形。在被解释为周期阶坎之前,内潮39、蠕动40和超临界流都被认为是新月形底形的成因。近些年,通过形态动力学分析,超临界流成因更加符合新月形底形的形态特征。Covault等41利用自主式水下航行器(Automatic Underwater Vehicle,AUV)多波束和浅地层剖面数据,揭示了周期阶坎的形态特征,呈现长波形、不对称、向上游迁移的新月形(图3)。海底峡谷和水道内周期阶坎都呈现新月形的外形特征,主要原因是超临界流在峡谷、水道横截面上的侵蚀能力差异。峡谷和水道的横截面呈现U形或V形,在中泓线上的超临界流流速最大、侵蚀沉积作用最明显43,超临界流的侵蚀能力由中泓线向峡谷、水道边界逐渐减弱,形成向上游迁移的新月形特征(图3a~c)。

图3

图3   不同环境下的周期阶坎(据参考文献[12,13,41,42]修改)

(a)加拿大斯阔米什三角洲;(b) 美国西海岸圣马特奥水道;(c)中国南海东北陆坡西澎湖峡谷;(d)墨西哥湾巴哈马碳酸盐台地

Fig.3   Cyclic steps in different environments (modified after references[12,13,41,42])

(a) Squamish prodelta of Canada; (b) The San Mateo channel of the west coast of USA; (c) West Penghu Canyon on the northeastern continental slope of the South China Sea; (d) The Great Bahama carbonate platform of the gulf of Mexico


根据超临界流对海底的侵蚀作用差异,周期阶坎分为净侵蚀型和净沉积型周期阶坎。净侵蚀型周期阶坎为一系列向上迁移的连续冲刷槽2244~47,而净沉积型周期阶坎为向上迁移的沉积物波2237434748。净侵蚀型周期阶坎的特点主要表现为闭合的等高线、较好的对称性、不连续的内部结构以及侵蚀基底,净沉积型周期阶坎的特点为不闭合的等高线、对称性较差、层次分明且连续以及不侵蚀基底49。在海底峡谷和水道中,周期阶坎往往从上游的净侵蚀型开始,随着泥沙含量降低,逐渐由净侵蚀型转变为净沉积型周期阶坎,在下游区域形成净沉积型周期阶坎。

2.3 层理结构

超临界流底形的层理结构受泥沙浓度和颗粒大小的影响9。峡谷和水道等限制性区域内超临界流底形具有明显分层的层理结构。在峡谷和水道上游,超临界流的泥沙浓度高、流速快、颗粒大,浊积岩相发育Ta、Tb4、T3a、T3b和Tb2等鲍马层序50图4),形成波长短(0~10 m)、波高较小(1 m)的周期阶坎8。在天然堤漫流区和水道—朵叶体过渡带内,超临界流的泥沙浓度低、流速慢、颗粒小,浊积岩相主要由Ta、T3b、Tb2、Tb1和Tc组成50,形成波长长(103 m)、波高较大(102 m)的周期阶坎44

图4

图4   周期阶坎岩相分布(据参考文献[50]修改)

Fig.4   The distribution of cyclic steps facies (modified after reference[50])


河流入海口快速沉积形成的周期阶坎,往往具有不对称性的泥岩或砂岩内部结构,呈现由超临界流形成的加积型沉积,指示了河流入河口处的超临界流体,在失去河流约束后,先呈扇形的快速展开、然后急剧减速的流体特征15。浊流等超临界流的泥沙含量具有明显的分层结构,上层浓度低,底层浓度高5152。浊流成因的周期阶坎,背流面遭受侵蚀,形成向斜层理;迎流面形成成层性较好的背斜层理2,每一层沉积呈现垂向上从细粒砂岩到中粒砂岩,最后为粗粒砂岩的浊积岩沉积特征18图5)。

图5

图5   周期阶坎粒度分布(据参考文献[18]修改)

Fig.5   The grain size distribution of cyclic steps (modified after reference[18])


3 周期阶坎的形成机理

在泥沙运动力学理论基础之上,进行了水槽实验和数值模拟,建立了周期阶坎的形态动力学模型。

3.1 水槽实验

周期阶坎的形成主要是超临界流与海底之间的相互作用。高流态的底形是一种瞬时、不稳定的状态53。沉积物波可以由先成底形产生47,周期阶坎可能由其他底形转换而来。Cartigny等21通过水槽实验揭示逆行沙丘、不稳定沙丘、流槽—冲坑和周期阶坎之间的依次过渡关系。当超临界流速度增加、泥沙浓度增加时,从短波长的逆行沙丘逐渐向更长波长的不稳定逆行沙丘、流槽—冲坑过渡,最后形成长波状周期阶坎。形成周期阶坎的超临界流具有如下特征:在速度结构上,超临界流经过水跃,流速减小,流体厚度增加,形成亚临界流23在湍流动能上,超临界流流速大、动能高,亚临界流流速小、动能低23在泥沙浓度上,迎流面泥沙浓度比背流面泥沙浓度低2343在剪切应力上,超临界流床层剪切应力高,亚临界流床层剪切应力低2354

在一个完整的周期阶坎内部,背流面的超临界流持续加速,导致了床层剪切应力的增加,泥沙夹带速度大于泥沙沉降速度,发生侵蚀作用23(图45)。然而,在迎流面,超临界流的流速降低,转化成亚临界流,剪切应力逐渐降低,侵蚀能力弱23,导致泥沙沉降速度超过泥沙夹带速度,产生沉积55。在背流面与迎流面之间,超临界流会向亚临界流转化,发生湍流动能消散,产生水跃23,较高的剪切应力对海底产生侵蚀23。经多次水跃,逐渐形成了向上迁移的新月形周期阶坎404456

3.2 数值模拟

通过改变数值模拟过程中的不同参量,模拟超临界流和底形的相互关系、周期阶坎发育情况。根据模拟维度,主要有一维的“平均四方程”模型和二维的FLOW-3D软件模拟。

3.2.1 “平均四方程”模型

该模型考虑到泥沙浓度、初始地形的斜率等变量对实际周期阶坎演化的影响。Kostic22通过数值模拟技术,阐明了周期阶坎发育、迁移和结构的控制因素。采用了Kostic等24所提出的“平均四方程”公式,用无量纲参数N描述控制周期阶坎发育、迁移和结构:

N=fFrdRfSCPλpLbreakLvsUΔτrcf

式中:Frd为密度弗洛德数;Rf为颗粒无量纲沉降速率;S为初始底形斜率;C为流体浓度;P为底形抗侵蚀能力的系数;λp为底形孔隙度;Lbreak/L为坡折处长度相对于整个研究区长度的比值;vs/U为颗粒沉降速率与流速比值;Δτ为底形的希尔兹数;r为常数;cf为底形摩擦系数。

Covault等41采用Kostic22的形态动力学模型模拟了圣马特奥水道内的6个底形(图3b),模拟结果表明该系列底形为净沉积周期阶坎(图6)。

图6

图6   沿圣马特奥水道浊流形态动力学模拟结果[41]

(a)底形和浊流自由液面沿着水道中泓线的变化; (b)密度弗洛德数的变化

Fig.6   Results of turbidity current morphodynamic modelling along the San Mateo channel[41]

(a) Bed elevation and turbidity-current interface change along the channel thalweg in response; (b) Variation in the densimetric Froude number


3.2.2 FLOW-3D模拟

Vellinga等23通过使用FLOW-3D建立的计算流体力学模型,模拟了Cartigny等21水槽实验低加积周期阶坎的形态动力学和沉积特征。计算流体力学模型采用了平均雷诺纳韦斯托克流体控制方程,使用单个泥沙运输模型计算了推移质运输和悬移质运输方程。低加积条件下的模拟结果与水槽实验结果相吻合,反映了超临界流底形的发育情况(图7)。

图7

图7   FLOW-3D初始模型图和效果图[23]

Fig.7   FLOW-3D initial model diagram and rendering[23]


4 不同研究手段之间的尺度差异

近些年,随着海底周期阶坎研究的深入,对于其形态和结构的精细化有着更高的要求。目前,对于周期阶坎的研究方法主要有水下观测、野外露头以及数值模拟。不同的观测手段揭示的尺度不同,从而影响对周期阶坎形态结构的认知差异。

4.1 探测设备的差异性

海底周期阶坎形态特征和内部结构的观测数据,主要通过船载多波束、浅地层剖面和人工地震以及AUV或ROV多波束、浅地层剖面获得。人工地震、多波束和浅地层剖面可以观测到周期阶坎的波长、波高等关键形态以及内部层理结构,作业水深、搭载平台和探测设备精度会影响底形观测。对于多波束和浅地层剖面而言,其分辨率一方面取决于设备本身,另一方面取决于观测区域的水深。Clarke57提出典型声纳系统的水平分辨率和垂直分辨率分别为水深的3.5%~20.0%和0.2%~0.8%。因此,在深水区域,船测数据分辨率将降低,不能揭示垂向尺度小于1 m的微观底形特征。例如,平均水深4 000~5 000 m的加拿大东部海域,船载多波束测深数据的垂向分辨率为40 m(表2)。因此,对于形态较小的深水周期阶坎(波长小于100 m,波高小于6 m)的观测可能会存在偏差。人工地震数据的分辨率一方面由频带控制(垂直分辨率)和道间距控制(水平分辨率);另一方面又受沉积层的速度结构影响,在同等频带条件下,沉积层的速度越高,垂向分辨率越低。根据前人研究,深水海底周期阶坎主要是由泥岩和砂岩组成的沉积物,地层速度约为1 500 m/s,如果地震的主要频带为30~70 Hz,地震分辨率为5.36~12.50 m。例如,南海西澎湖峡谷内的周期阶坎,垂直分辨率为11~17 m(表2)。随着水深增加,尤其在水深大于500 m的深水区域,海水的能量吸收作用增强,船测数据的频率都会降低,垂向分辨率会减小。在深水区域,船测多波束、浅地层剖面和人工地震数据可能不足以揭示周期阶坎的外部形态和内部层理结构。随着AUV技术的发展,搭载多波束和浅地层剖面仪器潜入海中,始终与海底保持恒定距离(35~50 m),获取海底浅层数据。AUV浅地层剖面的探测精度能达到分米级别,可以揭示几乎与野外露头尺度一致的周期阶坎(表2)。

表2   对比不同探测设备之间的差异性

Table 2  Comparison of the differences between different detection equipments

研究区域水深/m探测设备及相关参数周期阶坎特征(波长L,波高H
作业方式人工地震多波束浅地层剖面
大巴哈马台地西南斜坡[13,14]600~800船测震源:mini GI枪;垂向分辨率:1.5~2.5 m网格尺寸:7.5 m探测深度:海底以下0~50 m,垂向分辨率:1 m

L:96~117 m

H:3.5~4.3 m

美国San Metro水道[41]200~800AUV搭载,距离海底50 m-垂向分辨率:0.3 m,水平分辨率:1.5 m探测深度:海底以下0~30 m,垂向分辨率:0.1 m

L:300~1 100 m

H:20~50 m

非洲Rio Muni盆地[49]1 660~2 080船测主频:50 Hz,垂向分辨率:10 m--

L:506~1 374 m

H:10~41 m

南海西澎湖峡谷[11]2 145~2 950船测主频:30~45 Hz,垂向分辨率:11~17 m垂向分辨率:优于水深0.5%-

L:1 347~4 730 m

H:17~146 m

加拿大东部海域[47]4 000~5 000船测、深拖深拖,电火花震源主频:1.5 kHz,垂向分辨率:1.5 m;船测,震源:GI枪,垂向分辨率:4~10 m船测,网格尺寸:60 m船测,探测深度:海底以下0~5 m,垂向分辨率:1 m

L:600~4 000 m

H:10~50 m

注:“-”表示未采集数据

新窗口打开| 下载CSV


4.2 露头观测的局限性

周期阶坎的野外露头主要是由古浊流形成的浊积岩,根据剖面可以分析沉积物粒度、恢复周期阶坎的波长和波高等参数,从而估算古浊流的流速、流体厚度等流体参数。通过野外露头观测到的周期阶坎的尺度往往很小,甚至能够达到厘米级别,不能显示中泓线上全部周期阶坎的形态演化,具有片面性,容易引起观测结果与实际情况之间的差异。中泓线是超临界流的流动方向,水道或者海底峡谷的中泓线是研究周期阶坎内部结构的理想场所,更能揭示真实的周期阶坎特征。值得注意的是,周期阶坎的波长、波高等参数和流体流动状态是根据中泓线上的底形特征推测确定。然而,露头剖面受到出露部位和风化程度的限制,极有可能不是沿着河道/水道的中泓线。因此,野外观察到的周期阶坎露头剖面极有可能与中泓线上的周期阶坎之间存在误差(图8)。由此得出的周期阶坎特征参数和流体参数,就与形成周期阶坎底形特征的真实流体之间,存在一定偏差。为了确定真实中泓线位置,可以采用探地雷达或者人工地震,对野外露头的内部结构进行三维几何成像。例如,Lang等17利用高分辨率探地雷达对古冰川环境下的周期阶坎沉积结构三维成像,并且沿流动方向对周期阶坎横向和纵向的时空演化进行了解释。此外,还可以通过结合钻井、测井等方式对周期阶坎内部形成的浊积岩进行岩性分析,从而确定中泓线的位置。

图8

图8   周期阶坎野外露头剖面线和中泓线示意图

Fig.8   Line of outcrop profile and thalweg of cyclic steps


4.3 数值模拟的实效性

数值模拟主要是恢复周期阶坎发育演化特征以及流体状态。包括深水浊流在内的超临界流是由滑塌、大波浪、内潮和驻波等多种因素触发下引起58~60,然而,这些触发因素都有很大的不确定性。“平均四方程”模型可以模拟周期阶坎形成过程中的自由液面与底形的相互作用,模拟的浊流流动状态持续,泥沙浓度等参量不变。然而,“平均四方程”模型没有体现流体内部和底形之间的相互作用(表3)。FLOW-3D软件可以有效模拟水槽实验中形成的周期阶坎,模拟精度可以达到厘米级别(表3)。FLOW-3D软件可以清楚地展示流体内部的浓度、流速以及底形的动态变化,但是,在实际海底周期阶坎的模拟应用中还没有取得进展。不仅如此,目前数值模拟方法的维度仅仅停留在二维(xy)空间,在三维(xyz)空间上仍然没有取得突破。

表3   对比不同模拟方法之间的差异性

Table 3  Comparison of the differences between different simulation methods

模拟方法优缺点
数值模拟平均四方程[22]模拟自由液面与底形的作用关系,但是,无法反馈流体内部的动态变化和实际底形的精细特征
FLOW-3D模拟流体和底形的相互作用,揭示泥沙浓度、流速、湍流动能、弗洛德数、剪切力等参数的动态变化,但是,尚未在实际周期阶坎模拟方面取得突破

新窗口打开| 下载CSV


5 结论与建议

周期阶坎作为超临界流的底形表征,对于深水沉积物的再次搬运、海底底形的重构有着重要作用。本文根据前人的研究,总结了海底周期阶坎不仅可以形成于碎屑岩沉积环境,例如,河口三角洲、峡谷—水道等沉积体系;而且,可以形成于碳酸盐台地周缘斜坡,例如,巴哈马台地。两种沉积环境的共同点是具有一定坡度,容易形成浊流和倾泻密度流等超临界流体。典型的周期阶坎具有新月形的形态特征,并且根据沉积性和侵蚀性可以分为净侵蚀和净沉积两大类。在海底峡谷和水道中,随着泥沙浓度和流体速度的逐渐降低,周期阶坎往往从上游的净侵蚀型向下游的净沉积型转变。采用水槽实验和数值模拟揭示周期阶坎的形成原理。着重讨论了不同观测技术手段之间的差异导致周期阶坎研究尺度的差异。

在未来研究方向上,周期阶坎的研究可能会在以下几个方面有所突破:

(1) 在水深大于500 m的深水区域,融合AUV数据、船测数据和人工地震揭示的不同尺度底形特征,更加准确获得深水周期阶坎的形态特征和内部结构。

(2) 使用载人潜水器对深水海底的周期阶坎进行原位观测和测试,获取类似野外露头的样品和观测数据。

(3) 加强对野外露头的研究,结合探地雷达、人工地震、钻井和测井等多种方法,确定中泓线的准确位置,显示周期阶坎的完整特征,完善浊积岩的沉积模式。

(4) 建立一套适用于一定地质时间内、多期次浊流流动状态下形成的周期阶坎模拟方程,并且模拟维度从二维向三维发展。

参考文献

Wang DaweiBai HongxinWu Shiguo.

The research progress of turbidity currents and related deep-water bedforms

[J]. Advances in Earth Science2018331): 52-65.

[本文引用: 2]

王大伟白宏新吴时国.

浊流及其相关的深水底形研究进展

[J]. 地球科学进展, 2018331): 52-65.

[本文引用: 2]

Slootman ACartigny M J B.

Cyclic steps: Review and aggradation-based classification

[J]. Earth-Science Reviews, 2019. DOI:10.1016/j.earscirev.2019.102949.

[本文引用: 4]

Clare M ACartigny M J BNorth L Jet al.

Quantification of near-bed dense layers and implications for seafloor structures: New insights into the most hazardous aspects of turbidity currents

[C]//Offshore Technology Conference2015.

[本文引用: 1]

Andrieux OCooper C KWood J.

Turbidity current measurements in the Congo Canyon

[C]//Offshore Technology Conference2013.

Cooper CWood JImran Jet al.

Designing for turbidity currents in the Congo Canyon

[C]//Offshore Technology Conference2016.

[本文引用: 1]

Postma GKleverlaan KCartigny M J B.

Recognition of cyclic steps in sandy and gravelly turbidite sequences, and consequences for the Bouma facies model

[J]. Sedimentology, 2014617): 2 268-2 290.

[本文引用: 1]

Ventra DCartigny M J BBijkerk J Fet al.

Supercritical-flow structures on a Late Carboniferous delta front: Sedimentologic and paleoclimatic significance

[J]. Geology, 2015438): 731-734.

Covault J AKostic SPaull C Ket al.

Cyclic steps and related supercritical bedforms: Building blocks of deep-water depositional systems, western North America

[J]. Marine Geology, 20173934-20.

[本文引用: 1]

Lang JBrandes CWinsemann J.

Erosion and deposition by supercritical density flows during channel avulsion and backfilling: Field examples from coarse-grained deepwater channel-levée complexes (Sandino Forearc Basin, southern Central America)

[J]. Sedimentary Geology, 201734979-102.

[本文引用: 2]

Lamb M PParsons J DMullenbach B Let al.

Evidence for superelevation, channel incision, and formation of cyclic steps by turbidity currents in Eel Canyon, California

[J]. Geological Society of America Bulletin, 20081203/4): 463-475.

[本文引用: 1]

Li SLi WeiAlves Tiago Met al. Large-scale scours formed by supercritical turbidity currents along the full length of a submarine canyon

northeast South China Sea

[C]//AGU Fall Meeting2019424.

[本文引用: 3]

Zhong GuangfaCartigny M J BKuang Zengguiet al.

Cyclic steps along the South Taiwan Shoal and West Penghu submarine canyons on the northeastern continental slope of the South China Sea

[J]. Geological Society of America Bulletin, 20151275/6): 804-824.

[本文引用: 3]

Betzler CLindhorst SEberli G Pet al.

Periplatform drift: The combined result of contour current and off-bank transport along carbonate platforms

[J]. Geology, 20144210): 871-874.

[本文引用: 4]

Wunsch MBetzler CLindhorst Set al.

Sedimentary dynamics along carbonate slopes (Bahamas archipelago)

[J]. Sedimentology, 2017643): 631-657.

[本文引用: 3]

Lang JWinsemann J.

Lateral and vertical facies relationships of bedforms deposited by aggrading supercritical flows: From cyclic steps to humpback dunes

[J]. Sedimentary Geology, 201329636-54.

[本文引用: 5]

Dietrich PGhienne J FNormandeau Aet al.

Upslope-migrating bedforms in a proglacial sandur delta: Cyclic steps from river-derived underflows?

[J]. Journal of Sedimentary Research, 2016861): 112-122.

Lang JSievers JLoewer Met al.

3D architecture of cyclic-step and antidune deposits in glacigenic subaqueous fan and delta settings: Integrating outcrop and ground-penetrating radar data

[J]. Sedimentary Geology, 201736283-100.

[本文引用: 1]

Gong ChenglinChen LiuqinWest L.

Asymmetrical, inversely graded, upstream-migrating cyclic steps in marine settings: Late Miocene-early Pliocene Fish Creek-Vallecito Basin, southern California

[J]. Sedimentary Geology, 201736035-46.

[本文引用: 3]

Postma GHoyal D CAbreu Vet al.

Morphodynamics of supercritical turbidity currents in the channel-lobe transition zone

[C]//Submarine Mass Movements and Their Consequences201641469-478.

[本文引用: 2]

Taki KParker G.

Transportational cyclic steps created by flow over an erodible bed. Part 1. Experiments

[J]. Journal of Hydraulic Research, 2005435): 488-501.

[本文引用: 1]

Cartigny M J BVentra DPostma Get al.

Morphodynamics and sedimentary structures of bedforms under supercritical-flow conditions: New insights from flume experiments

[J]. Sedimentology, 2014613): 712-748.

[本文引用: 5]

Kostic S.

Modeling of submarine cyclic steps: Controls on their formation, migration, and architecture

[J]. Geosphere, 201172): 294-304.

[本文引用: 6]

Vellinga A JCartigny M J BEggenhuise J Tet al.

Morphodynamics and depositional signature of low-aggradation cyclic steps: New insights from a depth-resolved numerical model

[J]. Sedimentology, 2018652): 540-560.

[本文引用: 12]

Kostic SParker G.

The response of turbidity currents to a canyon-fan transition: Internal hydraulic jumps and depositional signatures

[J]. Journal of Hydraulic Research, 2006445): 631-653.

[本文引用: 2]

Sun TParker G.

Transportational cyclic steps created by flow over an erodible bed. Part 2. Theory and numerical simulation

[J]. Journal of Hydraulic Research, 2005435): 502-514.

Parker GIzumi N.

Purely erosional cyclic and solitary steps created by flow over a cohesive bed

[J]. Journal of Fluid Mechanics, 2000419203-238.

[本文引用: 1]

Xu XiaoyongFuliang Wang Daweiet al.

Cyclic steps and significance to deep water sedimentation

[J]. Marine Origin Petroleum Geology2018234): 1-14.

[本文引用: 1]

许小勇吕福亮王大伟.

周期性阶坎的特征及其对深水沉积研究的意义

[J]. 海相油气地质, 2018234): 1-14.

[本文引用: 1]

Piper D JNormark W R.

Processes that initiate turbidity currents and their influence on turbidites: A marine geology perspective

[J]. Journal of Sedimentary Research, 2009796): 347-362.

[本文引用: 1]

Talling P JPaull C KPiper D J.

How are subaqueous sediment density flows triggered, what is their internal structure and how does it evolve? Direct observations from monitoring of active flows

[J]. Earth-Science Reviews, 2013125244-287.

Talling P JAllin JArmitage D Aet al.

Key future directions for research on turbidity currents and their deposits

[J]. Journal of Sedimentary Research, 2015852): 153-169.

[本文引用: 1]

Allen J. Sedimentary Structures, Their Character and Physical Basis[M]. New YorkElsevier2011.

[本文引用: 1]

Guy H PSimons D BRichardson E V. Summary of alluvial channel data from flume experiments1956-61[C]//U.S. Geological Survey Professional Paper, 19664621): 11-196.

[本文引用: 1]

Carling P AShvidchenko A B.

A consideration of the dune: Antidune transition in fine gravel

[J]. Sedimentology, 2002496): 1 269-1 282.

[本文引用: 1]

Alexander JBridge J SCheel R Jet al.

Bedforms and associated sedimentary structures formed under supercritical water flows over aggrading sand beds

[J]. Sedimentology, 2001481): 133-152.

[本文引用: 1]

Slootman ADe Boer P LCartigny M J Bet al.

Evolution of a carbonate delta generated by gateway‐funnelling of episodic currents

[J]. Sedimentology, 2019664): 1 302-1 340.

[本文引用: 3]

Kostic SSmith I B.

Water on Mars: Do submarine cyclic steps exist on the red planet?

[J]. Progress in Earth and Planetary Science, 201851): 76.

[本文引用: 1]

Clarke J H.

First wide-angle view of channelized turbidity currents links migrating cyclic steps to flow characteristics

[J]. Nature Communications, 20167. DOI:10.1038/ncomms11896.

[本文引用: 2]

Schnyder J S DEberli G PBetzler Cet al.

Morphometric analysis of plunge pools and sediment wave fields along western Great Bahama Bank

[J]. Marine Geology, 201839715-28.

[本文引用: 1]

Smith D PRuiz GKvitek Ret al.

Semiannual patterns of erosion and deposition in upper Monterey Canyon from serial multibeam bathymetry

[J]. GSA Bulletin, 20051179/10): 1 123-1 133.

[本文引用: 1]

Paull C KUssler WCaress D Wet al.

Origins of large crescent-shaped bedforms within the axial channel of Monterey Canyon, offshore California

[J]. Geology, 201066): 755-774.

[本文引用: 2]

Covault J AKostic SPaull C Ket al.

Submarine channel initiation, filling and maintenance from sea-floor geomorphology and morphodynamic modelling of cyclic steps

[J]. Sedimentology, 2014614): 1 031-1 054.

[本文引用: 6]

Clarke J HBrucker SMuggah Jet al. Temporal progression and spatial extent of mass wasting events on the Squamish prodelta slope

in Landslides and engineered slopes: Protecting society through improved understanding

[M]//Landslides and Engineered Slopes: Protecting Society through Improved Understanding. LondonCRC Press2012.

[本文引用: 1]

Cartigny M J BPostma GVan den Berg J Het al.

A comparative study of sediment waves and cyclic steps based on geometries, internal structures and numerical modeling

[J]. Marine Geology, 20112801/4): 40-56.

[本文引用: 3]

Symons W OSumner E JTalling P Jet al.

Large-scale sediment waves and scours on the modern seafloor and their implications for the prevalence of supercritical flows

[J]. Marine Geology, 2016371130-148.

[本文引用: 3]

Maier K LFildani APaull C Ket al.

Deep‐sea channel evolution and stratigraphic architecture from inception to abandonment from high‐resolution autonomous underwater vehicle surveys offshore central California

[J]. Sedimentology, 2013604): 935-960.

Spinewine BSequeiros O EGarcia M Het al.

Experiments on wegde-shaped deep sea sedimentary deposits in minibasins and/or on channels levees emplaced by turbidity currents. Part II. Morphodynamic evolution of the wedge and of the associated bedforms

[J]. Journal of Sedimentary Research, 2009797/8): 608-628.

Normandeau ACampbell D CCartigny M J B.

The influence of turbidity currents and contour currents on the distribution of deep‐water sediment waves offshore eastern Canada

[J]. Sedimentology, 2019. DOI:10.1111/sed.12557.

[本文引用: 4]

Kostic SSequeiros OSpinewine Bet al.

Cyclic steps: A phenomenon of supercritical shallow flow from the high mountains to the bottom of the ocean

[J]. Journal of Hydro-environment Research, 201034): 167-172.

[本文引用: 1]

Li LGong C.

Gradual transition from net erosional to net depositional cyclic steps along the submarine distributary channel thalweg in the rio muni basin: A joint 3‐D seismic and numerical approach

[J]. Journal of Geophysical Research: Earth Surface, 20181239): 2 087-2 106.

[本文引用: 2]

Postma GCartigny M J B.

Supercritical and subcritical turbidity currents and their deposits—A synthesis

[J]. Geology, 20144211): 987-990.

[本文引用: 4]

Fisher R V.

Flow transformations in sediment gravity flows

[J]. Geology, 1983115): 273-274.

[本文引用: 1]

Postma GNemec WKleinspehn K L.

Large floating clasts in turbidites: A mechanism for their emplacement

[J]. Sedimentary Geology, 1988581): 47-61.

[本文引用: 1]

Kostic SCasalbore DChiocci Fet al.

Role of upper-flow-regime bedforms emplaced by sediment gravity flows in the evolution of deltas

[J]. Journal of Marine Science and Engineering, 201971): 5.

[本文引用: 1]

Izumi NYokokawa MParker G.

Incisional cyclic steps of permanent form in mixed bedrock-alluvial rivers

[J]. Journal of Geophysical Research: Earth Surface, 20171221): 30-152.

[本文引用: 1]

Mutti ENormark W R.

Comparing examples of modern and ancient turbidite systems: Problems and concepts

[M]//Marine Clastic Sedimentology. DordrechtSpringer19871-38.

[本文引用: 1]

Paull CNormark W RUssler Wet al.

Large wave-shaped bedforms in the axial channel of Monterey Submarine Canyon: Monterey Bay, California

[C]//AGU Fall Meeting Abstracts. 2007.

[本文引用: 1]

Clarke D.

The effect of fine scale seabed morphology and texture on fidelity of Swath bathymetric sounding data

[J]. Sea Technol, 19983987-90.

[本文引用: 1]

Paull CUssler WGreene H Get al.

Caught in the act: The 20 December 2001 gravity flow event in Monterey Canyon

[J]. Geo-Marine Letters, 2002224): 227-232.

[本文引用: 1]

Inman D LNordstrom C EFlick R E.

Currents in submarine canyons: An air-sea-land interaction

[J]. Annual Review of Fluid Mechanics, 197681): 275-310.

Konsoer KZinger JParker G.

Bankfull hydraulic geometry of submarine channels created by turbidity currents: Relations between bankfull channel characteristics and formative flow discharge

[J]. Journal of Geophysical Research: Earth Surface, 20131181): 216-228.

[本文引用: 1]

/