海底水道体系沉积构型样式及控制因素:以尼日尔三角洲盆地陆坡区为例*
赵晓明1,2, 刘丽1,2, 谭程鹏1,2, 范廷恩3, 胡光义3, 张迎春3, 张文彪4, 宋来明3
1 天然气地质四川省重点实验室,四川成都 610500
2 西南石油大学地球科学与技术学院,四川成都 610500
3 中海油研究总院有限责任公司,北京 100027
4 中国石化石油勘探开发研究院,北京 100083

第一作者简介 赵晓明,男,1982年生,2011年毕业于中国石油大学(北京),获得博士学位,现为西南石油大学副教授,主要从事油气田开发地质方面的研究工作。E-mail: zhxim98@163.com

摘要

作为当今油气勘探开发的重要目标,海底水道的沉积构型样式复杂多变,其控制因素目前尚无定论,这大大增加了相关油气藏的钻探风险和开发难度。以尼日尔三角洲盆地陆坡区为例,以近海底高频地震信息为原型,揭示了海底水道体系的外部形态样式和内部结构样式,探究了水道体系与复合水道形态特征参数间的相互关系,分析了海底水道体系沉积构型控制因素。研究认为: ( 1)水道体系可分为限制性、半限制性和非限制性 3个大类,限制性水道体系的边界为一明显的大型下切谷界面(或峡谷),半限制性水道体系的边界也发育明显的大型下切谷,但其界面两侧发育大型天然堤沉积,非限制性水道体系则不发育大型下切谷,其中每个大类又可依据其内部复合水道的类型,细分为 2个小类(下切式和包络式)。( 2)水道体系内部复合水道间存在垂向和侧向 2种类型的叠置,前者可细分为孤立式、叠加式和切叠式 3类,后者则包括叠合式和分离式 2类。( 3)地形坡度与水道体系类型和规模具有一定的耦合与关联,在类似油藏规模的局部沉积体内部,陡坡段(坡度大于)主要发育限制性水道体系,缓坡段(坡度 0.5°~1°)主要发育半限制性水道体系,平坦段(坡度小于 0.5°)主要发育非限制性水道体系;随着地形坡度的增大,水道体系的宽度减小、深度增加、宽深比减小,即向“窄深型”形态发展。该项研究不仅对深水沉积学的发展具有较大的理论意义,而且对高效开发海底水道油气藏具有重要的实际意义。

关键词: 海底水道; 沉积构型; 构型样式; 地形坡度; 演化模式
中图分类号:P512.2 文献标志码:A 文章编号:1671-1505(2018)05-0825-16
Styles of submarine-channel architecture and its controlling factors: A case study from the Niger Delta Basin slope
Zhao Xiao-Ming1,2, Liu Li1,2, Tan Cheng-Peng1,2, Fan Ting-En3, Hu Guang-Yi3, Zhang Ying-Chun3, Zhang Wen-Biao4, Song Lai-Ming3
1 Key Laboratory of Natural Gas Geology of Sichuan Province,Chengdu 610500, Sichuan
2 School of Geoscience and Technology,Southwest Petroleum University,Chengdu 610500, Sichuan
3 CNOOC Research Center Co., Ltd.,Beijing 100027
4 Petroleum Exploration & Production Research Institute,SINOPEC,Beijing 100083;

About the first author Zhao Xiao-Ming,born in 1982,is an associate professor of Southwest Petroleum University. He is mainly engaged in research on hydrocarbon development geology. E-mail: zhxim98@163.com.

Abstract

As an important target of the present-day hydrocarbon exploration and development,submarine channels have complex styles of sedimentary architecture and the associated controlling factors are still poorly understood,which significantly raises the risk of drilling of the hydrocarbon reservoirs and the difficulty of development. Taking the Niger Delta slope as an example,the present study,utilizing the high-resolution seismic data from the shallow subsurface,reveals the styles of channel-system external morphology and internal construction,probes the relationships among morphological parameters of the channel system and the channel complex,and analyzes the controlling factors of the submarine-channel sedimentary architecture. It is suggested that: (1)The channel system could be divided into three categories, i.e., confined channel system,semi-confined channel system and unconfined channel system. The boundaries of the confined and semi-confined channel system are obvious erosional surfaces and the latter commonly develops levee deposits;the unconfined channel system,however,does not have the erosional surface. In addition,each kind of the channel system could be further divided into two subgroups(the incised one and the enveloped one). (2)The channel complex within the channel system has the agrradational type as well as the accreted type;the former consists of the isolated,superposed as well as cut-and-fill group and the latter includes the cut-and-fill and detached group. (3)The slope gradient has closed relationship with the channel-system type and scale;in the reaches having steep gradients(>1°)and gentle gradients(0.5°~1°),confined channel system and semi-confined channel systems are respectively dominated,however,unconfined channel system mainly shows up in the flat section of the slope(<0.5°). With the increasing of slope gradients,the widths and heights of the channel system would respectively decrease and increase and thus,resulting the rapid decrease of the ratio of width to height. That is,the channel system tends to be‘narrow and deep’. The current study leads to some new interpretations of deep-water sedimentology and is of great significance for the high-efficient development of associated reservoirs.

Key words: submarine channel; sedimentary architecture; architecture style; slope gradient; evolution pattern

作为深水沉积体系的重要单元(Wynn et al., 2007; Mchargue et al., 2011; 王琪等, 2017), 海底水道是主要的沉积物搬运通道及粗碎屑沉积场所(Khripounoff et al., 2003), 它们是陆缘盆地内重要的油气储集层(Kolla et al., 2001; Mayall et al., 2006; 康洪全等, 2013)。尽管这类储集层具有较高的孔隙度和渗透率, 但是它们的内部结构(连通性、几何形态、岩性)复杂多变, 即便是在很短的侧向距离内, 其储集层厚度和连通性也会有较大的变化(Wood and Mize-Spansky, 2009; 赵晓明等, 2012a, 2012b; 林煜等, 2013; Liu et al., 2013), 从而严重制约了该类油藏的高效开发(Alpak et al., 2013)。因此, 加强海底水道沉积构型研究, 即认知不同级次沉积构型单元的形态、规模、叠置样式及其控制因素, 有助于实现对海底水道砂体的内部结构进行钻前预知, 从而提高井位设计和开发概念设计的质量, 降低钻探风险, 评价储集层非均质性对油藏开发的影响(Mayall et al., 2006; Labourdette and Bez, 2010)。

然而, 当前海底水道构型研究远滞后于河流、冲积扇和三角洲等沉积类型。原因有以下几方面: 一是现代海底水道沉积不利于观察, 与多数陆相沉积不同, 深水现代沉积常位于水深大于300 m的深海环境, 难以用卫星照片分析、实地探槽测量等常规手段开展其沉积学研究; 二是野外露头信息有限, 深水沉积规模较大, 其宽度一般为数千米, 有的可达上万米, 如Amazon扇, 其峡谷最大宽度达15 km(Flood and Piper, 1997), 厚度多介于几十米至数百米之间, 属大尺度单元, 受出露条件的限制, 野外露头资料仅能提供低级次、小尺度构型单元的信息, 难以表征大尺度构型单元。

为此, 作者以尼日尔三角洲盆地陆坡泥岩底辟区为例, 通过对近海底高频地震信息的定量地震地貌学和沉积构型研究, 构建海底水道体系沉积构型样式, 探究其几何形态定量控制参数之间及其与地形坡度之间的耦合和关联, 揭示影响海底水道体系沉积构型发育的关键因素, 并探讨其演化机理。

1 地质背景

研究区位于西非尼日尔三角洲大陆边缘几内亚湾内, 北部距哈科特港200 km, 水深1300~1800 m(图1-A), 有丰富的海底扇沉积, 是研究深海水道的典型区域。

图1 尼日尔三角洲盆地地理位置及构造地层剖面Fig.1 Geographic location of the Niger Delta Basin and structural stratigraphic section

1.1 构造特征

尼日尔三角洲是世界上最大的退积三角洲之一(Doust and Omatsola, 1989), 东为喀麦隆火山岩地区, 西为达荷美盆地, 其内部构造受到沿洋壳的Charcot断裂带等断裂系统控制。尼日尔三角洲盆地为典型的被动大陆边缘盆地, 其演化经历了早白垩世— 晚白垩世裂谷期和古新世以来的漂移期2个阶段, 早始新世以来, 长期海退形成了现今的尼日尔三角洲(张兴, 2000; 邓荣敬等, 2008; 吕明等, 2008)。三角洲前积推进过程中, 受大陆边缘重力作用的影响(Cohen and McClay, 1996), 自北向南依次形成了拉张构造区、底辟构造区和推覆构造区(Damuth, 1994; Morley and Guerin, 1996; Zhang et al., 2016)(图1-B)。

研究区位于尼日尔三角洲区域构造转换带上, 即底辟构造区, 其北侧为拉张构造区, 受铲状断层和滚动背斜影响; 南侧为推覆构造区, 主要受控于叠瓦状逆冲断层(Morley and Guerin, 1996)(图1-B)。而该构造转换带位于陆坡和深海平原地区, 受阿格巴达组逆冲断层错断和阿卡塔组拆离面影响, 欠压实陆坡泥岩缺乏侧向支撑而形成该构造转换带。中新世重力滑动作用开始后, 大量的断层相关褶皱后翼形成背驮式微盆地。重力滑动作用使得该构造带具有复杂的构造样式, 包括断层转折褶皱、逆冲叠瓦系列褶皱等。

1.2 地层特征

尼日尔三角洲盆地发育白垩系到全新世海相碎屑岩地层, 这些地层覆盖在陆壳的海相碎屑岩之上。始新世以来, 沉积物逐渐受到波浪作用的控制, 碎屑沉积物开始在尼日尔三角洲盆地沉积。尼日尔三角洲盆地古新世至现代的地层从下到上分别为阿卡塔组(Akata)、阿格巴达组(Agbada)和贝宁组(Benin)(Short and Stä uble, 1967)(图1-B)。

基底的阿卡塔组是海相地层, 贯穿整个盆地, 岩性主要为灰黑色页岩, 局部含浊流或大陆斜坡深谷沉积, 为尼日尔三角洲盆地最主要的生油层, 厚度为2000~7000 m(Doust and Omatsola, 1989)。

上覆的阿格巴达组是尼日尔三角洲中主要的油气储集单元, 开始沉积于始新世并一直沉积至今, 发育三角洲前缘、三角洲顶积层和河成三角洲沉积。岩性上表现为砂岩与泥岩、海相页岩的互层, 该层总厚度超过3500 m(Doust and Omatsola, 1989; Corredor et al., 2005)。

在海滩和沿岸地区, 贝宁组覆盖在阿格巴达组之上, 贝宁组由始新世到全新世的陆相沉积物组成, 岩性为砂岩夹薄层泥岩(或页岩), 其砂层可能为砂坝、河道和天然堤沉积, 为晚始新世到全新世陆上和浅海地区的沉积物, 厚度达2000 m(邓荣敬等, 2008)。

2 数据与方法
2.1 数据

本次研究采用近海底高频三维地震信息作为原型, 其具有时代较新、沉积遭受后期构造破坏小的特点, 能够真实反映水道原始形态及其充填过程, 利于正确地揭示水道构型发育的控制因素及其作用机理。地震数据处理面元为12.5 m× 12.5 m, 记录长度7.5 s, 采样间隔3 ms。频谱分析表明, 原型模型所在层段优势频宽5~90 Hz, 主频近70 Hz, 分辨率较高, 地层速度若按1900 m/s计算, 可分辨厚度为6 m左右, 这足以满足单一水道层次地貌的定量构型分析。

2.2 方法

由于小级别的构型单元分布受控于大级别构型单元, 因此, 层次划分是进行海底水道沉积构型研究的前提。针对深水体系, 不同学者提出了不同的构型界别(界面)划分方案。Mutti和Normark(1987)将现代浊积体系和古代浊积体系进行了比较, 首次建立了划分深水沉积界面的5级系统; Lamb等(2003)将加利福尼亚San Joaquin盆地Stevens地层富砂型深水沉积体系划分为6个构型级别; 赵晓明等(2010)针对西非尼日尔三角洲盆地陆坡区阿格巴达组深水沉积的特点, 建立了该区深水海底扇的6级划分方案。但无论是哪种划分方案, 海底水道沉积均可由小往大归结为单一水道、复合水道和水道体系3个层次。

文中研究对象为水道体系及其内部的复合水道层次。研究过程中, 首先提取近海底浅层高频地震数据中, 能反映浅层地貌特征的地震属性(包括均方根振幅和相干属性等), 然后利用定量地震地貌学和沉积构型分析相结合的方法, 通过井震标定, 明确水道体系的外部形态样式和内部结构样式, 并测量水道体系及其内部保存完整的复合水道的几何形态定量控制参数, 如水道体系和复合水道的宽度、深度和曲率等, 进而结合地貌形态和深泓线坡度, 研究深海水道体系沉积构型发育的关键控制因素及其作用机理。

在地震属性地层切片的约束下, 精细解释垂直物源方向的地震剖面, 并在剖面上确定复合水道的顶面位置和侧向边界, 进而通过测量两者左右交点之间的距离来确定水道体系及其内部复合水道的宽度。在地震剖面测量宽度的同时, 同步测量水道时间深度, 即水道下切最深谷底处与水道顶之间的垂直时间距离, 进而利用目的层1900 m/s的地层速度, 获取水道深度。

3 水道体系构型样式

受侵蚀能力的影响, 不同水道体系和同一水道体系不同位置处, 其沉积构型样式差异较大, 下面从外部形态样式和内部结构样式2方面论述。

3.1 水道体系外部形态样式

依据水道体系在浅层高频地震剖面上的几何外形和内部反射特征, 可将其分为限制性、半限制性和非限制性3大类。

3.1.1 限制性水道体系

限制性水道体系以发育大型下切谷为典型特征, 下切谷两侧不发育或发育不明显天然堤。限制性水道体系边界处地震同相轴明显错断, 推断为大型下切谷界面, 整体呈U型或V型; 水道体系内部以杂乱状或叠瓦状地震反射结构为主, 振幅中— 强, 见弱振幅充填, 同相轴连续性差— 中等(图2), 推断体系内部发育多期砂质和(或)泥质复合水道。

图2 尼日尔三角洲盆地陆坡区限制性下切式(A)和包络式(B)水道体系典型地震响应特征及发育模式Fig.2 Seismic responses and development pattern of the confined aggradational channel complex (A)and laterally-migrated channel complex(B) in the Niger Delta Basin slope

根据大型下切谷(峡谷)内部复合水道的发育样式, 该类沉积体可进一步细分为限制性下切式水道体系(图2-A)和限制性包络式水道体系(图2-B)。前者下切谷内部存在大量块状搬运物、滑塌物和过路沉积, 局部发育下切式复合水道; 后者下切谷内部以发育包络式复合水道为主, 局部发育滑塌物、过路及披覆沉积。

3.1.2 半限制性水道体系

半限制性水道体系也发育下切谷, 与限制性水道体系不同的是, 其下切谷两侧发育大型天然堤。半限制性水道体系边界处地震同相轴明显错断, 为大型下切谷界面, 整体呈U型或V型(图3); 水道体系两侧发育楔状体, 弱振幅弱连续反射, 为大型天然堤沉积, 规模约占整个水道体系的 2/3; 大型下切谷和天然堤组合在一起, 使得该类水道体系几何形状呈“ 海鸥型” ; 水道体系内部以杂乱状和叠瓦状地震反射结构为主, 振幅中— 强, 见弱振幅充填, 同相轴连续性差— 中等, 推断体系内部发育多期砂质和(或)泥质复合水道。

图3 尼日尔三角洲盆地陆坡区半限制性下切式(A)和包络式(B)水道体系典型地震响应特征及发育样式Fig.3 Seismic response and development patterns of the semi-confined the aggradational channel complex (A) and laterally-migrated channel complex (B) in the Niger Delta Basin slope

同样, 根据大型下切谷(峡谷)内部复合水道的发育样式, 该类沉积体可细分为半限制性下切式水道体系(图3-A)和半限制性包络式水道体系(图3-B)。前者在地震剖面上不仅存在大型下切谷, 而且内部有多期小型下切谷, 表明发育下切式复合水道, 同时可见杂乱反射, 推断为滑塌物沉积; 后者复合水道类型主要为包络式, 剖面上见典型的叠瓦状地震响应特征, 其间为弱— 强振幅连续地震反射结构, 推断为披覆、过路等细粒沉积。

3.1.3 非限制性水道体系

非限制性水道体系最大特点在于侵蚀作用比较弱而垂向加积作用比较强, 导致天然堤发育, 使得该类水道体系在地震剖面上往往呈“ 丘型” (图4); 丘状体的两翼主要为弱振幅弱连续地震反射, 局部夹有短波状强振幅反射, 说明水道沉积体两侧天然堤发育, 并见决口扇沉积; 丘状体的中部发育U型或V型水道, 它们之间存在侧向及垂向迁移, 其内部为弱、强振幅充填, 推断为砂质和泥质充填水道并存。

图4 尼日尔三角洲盆地陆坡区非限制性下切式(A)和包络式(B)水道体系典型地震响应特征及发育样式Fig.4 Seismic responses and development patterns of the unconfined channel the aggradational channel complex (A) and laterally-migrated channel complex(B) in the Niger Delta Basin slope

按照复合水道的发育类型, 同样非限制性复合水道细可分为下切式(图4-A)和包络式(图4-B)两种。前者主要由下切式复合水道侧向及垂向迁移形成, 后者主要由单一水道侧向及垂向迁移形成。复合水道两侧发育小型天然堤(水道天然堤), 依据距离水道的远近可将其分为近岸水道天然堤和远岸水道天然堤2种。

3.2 水道体系内部结构样式

由外部形态样式研究可知, 水道体系内部存在多期复合水道, 它们在空间上存在着多种叠置样式(内部结构样式), 这可概括为垂向叠置和侧向叠置2种样式。

3.2.1 垂向叠置样式

垂向叠置样式是指不同期的复合水道在纵向上的接触方式, 反映的是层间砂体的接触关系。复合水道垂向叠置样式可细分为孤立式、叠加式和切叠式3种类型(图5)。

图5 尼日尔三角洲盆地陆坡区复合水道垂向叠置模式的3种类型(据赵晓明等, 2012a)
A— 孤立式; B— 叠加式; C— 切叠式
Fig.5 Three types of vertical aggradation of channel complexes in the Niger Delta Basin slope (Zhao et al., 2012a)

孤立式叠置的2期复合水道在垂向尚未接触, 其间为细粒沉积, 地震剖面上表现为弱振幅充填, 测井曲线存在明显回返, 表现为高伽马、低电阻特征, 可作为流体纵向流动的遮挡层(图5-A)。

叠加式叠置的2期复合水道在垂向上已接触, 但无明显下切, 地震剖面上复合水道间可见连续稳定同相轴反射, 测井曲线上有回返, 为早期复合水道顶部细粒沉积或末期水道底部滞留沉积, 属于储集层物性差的夹层(图5-B)。

切叠式的2期复合水道在垂向上存在明显下切, 地震剖面上复合水道间纵向上呈明显的叠瓦状结构, 表明先期复合水道的上部被后期水道严重下切侵蚀, 仅仅保留了下部的不完整旋回, 测井曲线上有轻微回返, 为后期复合水道的底部滞留沉积, 该样式叠置关系可形成厚层复合砂体(图5-C)。

3.2.2 侧向叠置样式

侧向叠置样式表现为同层不同位复合水道在横向上的接触方式, 反映的是层内砂体侧向接触关系。复合水道存在叠合式和分离式2种侧向叠置样式(图6)。

图6 尼日尔三角洲坡区复合水道侧向叠置样式的2种类型
A— 复合水道平面分布特征; B— 叠合式侧向叠置样式; C— 分离式侧向叠置样式
Fig.6 Two types of lateral accretion of channel complexes in the Niger Delta Basin slope

叠合式复合水道平面上呈叠合状(图6-B), 同层不同位的2条复合水道剖面上呈侧向切叠关系, 即: 早期水道复合体(Ⅰ )的一侧大部分被后期复合水道(Ⅱ )强烈下切侵蚀, 其顶面可存在一定的高程差异, 但相差不大; 受重力作用和重力流动力学机理的影响, 水道内壁可能存在储集层质量差或非渗透性的滞留沉积和滑塌物, 这会减弱流体在不同复合水道间的侧向流动能力, 造成井间可能表现为弱或中等程度连通。

分离式的2条复合水道剖面上则呈互不接触(完全分离)关系(图6-C), 其间见弱— 中振幅反射, 为水道体系内部的溢岸天然堤、深海披覆泥等细粒沉积物, 造成两复合水道砂体间的互不连通。

4 水道体系定量特征

由水道体系构型样式研究成果可知, 非限制性水道体系不发育大型下切谷(峡谷), 它是由多期复合水道和单一水道空间叠置而成的, 故其定量规模主要取决于单一水道和复合水道的大小。下面重点探讨发育大型下切谷的限制性和半限制性水道体系的定量特征。

4.1 典型水道体系的优选

文中优选的水道体系是一条自北向南流动的浊积体系(图7), 其平面上形态复杂多变; 相干体属性地层切片(图7-A)表明, 水道体系有的部位呈顺直状, 有的则为蛇曲状, 个别地段出现截弯取直现象。

图7 尼日尔三角洲盆地陆坡区样本水道的相干属性和均方根属性图以及相关地震剖面Fig.7 Variance and RMS amplitude extractions and corresponding seismic sections of the target channel system in the Niger Delta Basin slope

水道体系内部包含3期相互切叠的复合水道。如图7-B所示的均方根振幅属性地层切片, 水道体系为中等、强振幅充填, 中等振幅呈条带状分布, 而强振幅呈梭状断续分布, 显然呈现出早期砂质复合水道被后期砂泥混杂质复合水道侵蚀切割的地震响应特征。水道体系两侧弱振幅充填, 推断为天然堤、披覆等泥质细粒沉积, 水道中部弯曲段见扇形中— 强振幅沉积体, 推断为决口扇(图7-B)。

从剖面叠置样式分析, 研究样本发育限制性和半限制性2种水道体系类型。如图7-C所示, 水道体系自北往南均发育明显大型下切谷, 只是顺物源方向不同位置处, 有的水道段下切谷两侧发育楔状体(图7-C, 剖面a-a'和c-c'), 有的下切谷两侧则为弱— 强振幅连续反射(推断为深海披覆细粒沉积)(图7-C, 剖面b-b'和d-d'), 前者为半限制性水道体系, 而后者为限制性水道体系。

利用地震属性平面和剖面互动的方法, 以约600 m的间距(图8-A), 测量了样本水道体系及其内部复合水道Ⅲ 的宽度、深度等数据(图8-B), 共获取了2组共130个数据点。结果表明, 样本水道体系的宽度为428.6~1968.0 m, 深度为80.9~152.0 m, 宽深比为3.96~18.15; 样本水道体系内部的复合水道Ⅲ 宽度为260~1000 m, 深度为50~142 m, 宽深比为3.18~9.31。

图8 尼日尔三角洲盆地陆坡区样本水道体系及其内部水道复合体测量示意图Fig.8 Target channel system and its corresponding measuring methods of internal channel complexes in the Niger Delta Basin slope

4.2 定量特征参数的相互关系

4.2.1 水道体系样本的分段

水道体系地形坡度的测量分析结果表明, 如图9所示顺物源方向, 水道体系样本可划分为A、B、C和D共4个坡度段(测量单元), 其坡度值分别为0.89° 、1.12° 、0.63° 和1.65° 。相比较复合水道坡度测量结果(对应段坡度分别为0.72° 、1.05° 、0.54° 和1.31° ), 尽管两者在测量单元分界点的平面位置相近, 但同一测量单元水道体系的地形坡度明显大于复合水道。其根本原因在于, 水道体系样本形成于早期的峡谷发育阶段, 受泥岩底辟、伴生断层等构造作用的影响, 会形成大小不一、高低起伏的小型隆起和凹陷, 这便使得古地形变化大; 到了水道体系样本发育的中后期, 由于受前期复合水道的填平补齐作用, 古地形的高低起伏程度大大减弱, 从而造成同一测量单元样本复合水道Ⅲ 的坡度相对较小。

图9 尼日尔三角洲盆地陆坡区样本水道体系顺物源方向地形坡度Fig.9 Diagram showing variation of the slope gradient of target channel system along source direction in the Niger Delta Basin slope

4.2.2 宽度与深度的关系

为论证样本水道体系宽度和深度规模的内在关联, 按以上4个坡度段对宽度、深度参数进行了统计, 如图10-A蓝圈中所示, 在坡度相对缓的A段和C段, 研究样本的宽度和深度呈正相关关系; 而在坡度陡的B段和D段, 如图10-A红圈中所示, 两者呈现出负相关关系。这意味着在缓坡区水道体系的宽度越大, 相应位置的下切深度也会增大, 而当地形坡度足够大时, 这一规律会发生逆转, 即: 宽度越大, 水道体系的下切能力越弱, 形成的深度越小。

图10 尼日尔三角洲盆地陆坡区样本水道体系宽度与深度关系统计图Fig.10 Statistical graph of the target channel system in the Niger Delta Basin slope showing the mathematical relations of widths and heights

4.2.3 宽度与宽深比的关系

同理, 文中也统计了样本水道体系宽度和宽深比的关系, 结果表明, 两者呈明显正相关性。如图10-B所示, 4个坡度段的宽度与宽深比均呈良好的正相关性, 且在相同宽度条件下, 坡度较陡的B段和D段宽深比相对较小, 主要分布在4~9之间, 有一半左右取样点宽深比小于6, 而坡度较缓的A段和C段宽深比相对较高, 一般都大于6。这说明同一水道体系在不同地形条件下, 陡坡区往往发育“ 窄深型” 水道, 而缓坡区易发育“ 宽浅型” 水道。

4.2.4 深度与宽深比的关系

受地形坡度变化的影响, 水道体系深度与宽深比关系较为复杂, 陡坡区深度与宽度呈负相关趋势, 而缓坡区两者呈正相关趋势。如图10-C红圈中所示, 对于坡度较陡的B段和D段, 深度增加, 宽深比有减小的趋势, 这说明陡坡区水道体系深度越大, 宽深比越小, 越呈“ 窄深型” 水道; 而在坡度较缓的A段和C段, 如图10-C蓝圈中所示, 呈现出与陡坡区相反的趋势, 即: 随着深度的增加, 宽深比有增加之势, 水道越往“ 宽浅型” 形态发展, 这一规律与下切式复合水道相似。

5 主要控制因素分析

水道体系规模较大, 影响其构型发育的地质因素类型多样(Stow, 1981; Kolla, 2007; 廖计华等, 2016; 陈宇航等, 2017), 诸如海平面的升降(Zhang et al., 2018)、沉积物的供给速度和类型(Gong et al., 2016a, 2016b)、区域构造运动(Clark and Cartwright, 2009)等, 均可造成水道体系的构型样式差异。就同一水道体系局部区域(如油藏范围)而言, 地形坡度是控制其形态特征的主要因素(胡孝林等, 2014; Zhao et al., 2018a, 2018b), 同时, 水道体系的形态也控制了其内部复合水道构型样式的发展演化, 揭示它们之间的相互关系对该类油气藏的高效开发具有重要理论意义。

5.1 坡度对水道体系形态的影响

5.1.1 坡度对宽度、深度的影响

利用前文测量的水道体系宽度、深度数据, 分不同坡度段, 统计了研究样本宽度和深度的平均值和最大值。前者代表了同一坡度条件下水道体系宽度和深度的平均大小, 后者代表了同一坡度条件下水道体系发育的最大延伸宽度和下切深度。如图11-A所示, 水道体系的坡度与宽度呈明显负相关性, 无论是水道体系的平均宽度还是最大宽度, 均随着地形坡度的增加而减小, 这说明同一水道体系局部范围内, 陡坡区的宽度往往小于缓坡区。造成这一现象的原因与复合水道类似, 即: 在缓坡区, 沉积地貌往往呈高低起伏特征, 这容易造成重力流沉积物在流动过程中产生离心力, 从而加大了其侧向迁移能力, 形成宽度较大河型; 而在陡坡区, 尽管也可能存在高低起伏地貌, 但由于地形坡度的增加, 此时重力作用一般远大于离心力作用的影响, 从而使得沉积物局部较强的垂向下切能力, 形成深度大、宽度小的水道体系。

图11 尼日尔三角洲盆地陆坡区样本水道体系宽度、深度、宽深比与坡度拟合图Fig.11 Statistical graph of the target channel system in the Niger Delta Basin slope showing the mathematical relations of slope gradients and widths, heights and ratios of the width to height of the target channel complex

水道体系的坡度与深度呈一定的正相关性。如图11-B所示, 尽管相关性不是很高, 但4个坡度段的深度随着地形坡度的增加呈现出变大的趋势, 这在坡度与平均深度的关系上更为明显。值得注意的是, 研究样本存在异常点, 即: 坡度最大的D段, 其水道体系的最大深度明显小于坡度小的A段和B段, 究其原因在于流动距离对重力流沉积物下切能力的影响; 水道体系样本纵贯研究区南北, 蜿蜒流动距离达40 km有余, 重力流沉积物在流动过程中侵蚀能量不断损失, 当其消耗能量大于由地形坡度增加而补充的能量时, 就会发生水道体系在陡坡区的下切深度小于相邻缓坡区的现象。

5.1.2 坡度对宽深比的影响

同样, 基于不同坡度段的采样点数据, 提取了水道体系宽深比的平均值和最大值, 并分析了坡度对它们的控制作用。结果表明, 随着地形坡度的增加, 水道体系不同段的宽深比平均值和最大值均相应减小(图11-C)。这意味着坡度对水道体系剖面形态具有良好的控制作用, 即: 同一水道体系局部范围内, 陡坡区易形成“ 窄深型” 水道, 而缓坡区则易形成“ 宽浅型” 水道。

5.1.3 坡度对水道体系类型的影响

坡度影响了水道体系的外部形态样式, 研究认为陡坡区容易发育限制性水道体系, 缓坡区易发育半限制性水道体系, 而在平坦区则主要发育非限制性水道体系。如图7-C所示, 水道体系样本在顺物源方向上, 外部形态样式不断发生变化; 在近物源端的A段(坡度0.89° ), 大型下切谷两侧发育天然堤, 为半限制性水道体系; 到了B段(1.12° ), 楔状天然堤消失, 水道体系类型转变为限制性; 再到C段(0.63° ), 水道体系又开始发育天然堤, 类型转化为半限制性; 而在远物源端的D段(1.65° ), 天然堤又消失, 形成限制性水道体系。这说明受地形坡度变化的影响, 水道体系类型会交替变化。

进一步分析认为, 水道体系之所以在缓坡区和陡坡区的外部形态样式存在差异, 与重力流沉积物的流动机理有关。当重力流沿着缓坡流动时, 沉积物在离心力作用下会产生单向流, 此时若能量足够大, 重力流便会在剥离作用下分离出一部分沉积物, 形成天然堤和溢岸沉积(Piper and Normark, 1983); 而当沉积物重力流进入坡度较陡的段时, 在重力作用下, 沉积物的流速不断增加, 侵蚀能力逐渐增强, 水道以下切侵蚀为主, 沉积物易发生过路现象, 不利于天然堤和溢岸沉积的发育(Kneller, 2003), 造成陡坡段多发育限制性水道体系。

5.2 水道体系形态对复合水道的影响

本次研究通过对处于同一位置处的水道体系与复合水道Ⅲ 的宽度、深度和宽深比取样数据进行拟合(图12), 发现水道体系与复合水道的宽度、深度以及宽深比各自表现为明显的正相关性, 即随着水道体系宽度、深度和宽深比的增大, 复合水道的宽度、深度和宽深比也相应增大。

图12 尼日尔三角洲盆地陆坡区样本水道体系与复合水道Ⅲ 宽度、深度和宽深比拟合图Fig.12 Statistical graph showing the mathematical relations of widths, heights and ratios of the width to height of the target channel system and channel complex Ⅲ in the Niger Delta Basin slope

深入分析认为, 首先, 这与水道体系和复合水道的坡度变化特征类似有关, 相似的坡度对水道体系和复合水道各形态学参数的控制, 使得二者的形态特征具有一定的相似性; 另外, 从深水水道体系的形成过程方面来说, 水道体系对复合水道的各参数也有一定控制作用。限制性和半限制性水道体系以发育大型下切谷为特征, 在水道体系形成早期, 较强的重力流作用侵蚀海底形成大型下切谷, 该下切谷的形成限制了复合水道横向的宽度发展和纵向的深度加大。因此, 下切谷较宽、较深的地方, 复合水道宽度、深度也相对较大; 相对的, 在下切谷较窄、较浅的地方, 复合水道宽度、深度受到限制, 也表现为窄和浅; 当然, 这也保证了水道体系和内部复合水道的宽深比具有一定的正相关性, 但是可能由于早期水道体系峡谷边界的侧向限制作用相对更为强烈, 所以导致复合水道的宽深比大多小于水道体系的宽深比。

5.3 水道体系构型演化模式

沉积物供应、区域盆地构造和相对海平面变化被认为是影响深水沉积体系的重要控制因素(Stow, 1981; Kolla, 2007)。当沉积物到达陆架边缘或超过陆架边缘时, 沉积物输送体系和汇聚盆地结构就成为控制深水沉积体系形态和分布的主要因素(Piper and Normark, 1983; 刘峰等, 2015; Gong et al., 2016a, 2016b; 陈志鹏等, 2017)。地形坡度是水道体系构型样式演化的关键, 它控制了水道体系在局部范围内的类型及样式的差异, 其成因可为泥岩底劈构造, 也可为逆冲断裂构造或者沉积侵蚀地貌。

基于定量地震地貌学研究成果, 建立了研究区坡度控制下的水道体系形态及其构型样式的演化模式(图13), 它也包含了水道体系形态学特征参数之间的关系, 具体内涵体现在以下几方面:

图13 坡度和曲率控制下水道体系沉积构型样式的演变模式Fig.13 Evolutionary model of channel system architecture under the control of sinuositie and slope gradient

1)坡度控制了水道体系的外部形态样式, 陡坡段(坡度大于1° )主要发育限制性水道体系, 缓坡段(坡度0.5° ~1° )主要发育半限制性水道体系, 平坦段(坡度小于0.5° )主要发育非限制性水道体系。

2)坡度控制了水道体系的宽度和深度, 随着地形坡度的增大, 水道体系的宽度呈减小趋势, 而深度则呈增加趋势。

3)坡度控制了水道体系宽深比, 随着地形坡度的增大, 水道体系的宽深比减小, 意味着水道向“ 窄深型” 形态发展。

4)同一坡度条件下, 水道体系的宽度与深度呈正比, 即: 宽度越大, 深度越大。

5)同一坡度条件下, 水道体系的宽度与宽深比呈明显正相关性。

6)同一坡度条件下, 水道体系的深度和宽深比呈2种截然不同的关系:在相对缓坡区, 两者呈正相关趋势, 而相对陡坡区两者呈负相关关系。

6 结论

1)总结了深水水道体系的外部形态样式和内部结构样式。水道体系的外部形态样式可分为限制性、半限制性和非限制性3个大类。限制性水道体系的边界为一明显的大型下切谷界面(或峡谷), 半限制性水道体系的边界也发育明显的大型下切谷, 但其界面两侧发育大型天然堤沉积, 非限制性水道体系则不发育大型下切谷。

2)综合复合水道的外部形态样式, 水道体系的每一个大类可进一步细分为2个小类, 如限制性水道体系可细分为限制性下切式和限制性包络式2类。水道体系内部结构样式, 可概括为垂向叠置和侧向叠置2种类型; 垂向叠置模式可细分为孤立式、叠加式和切叠式3类, 侧向叠置样式则包括叠合式和分离式2类。

3)揭示了限制性水道体系和半限制性水道体系与地形坡度、复合水道之间的关联。在局部范围内, 地形坡度控制了水道体系的形态特征; 首先, 坡度控制了水道体系的外部形态样式, 陡坡段(坡度大于1° )主要发育限制性水道体系, 缓坡段(坡度0.5° ~1° )主要发育半限制性水道体系, 平坦段(坡度小于0.5° )主要发育非限制性水道体系。其次, 坡度控制了水道体系的宽度、深度和宽深比, 随着地形坡度的增大, 水道体系的宽度会减小, 深度增加, 宽深比减小, 表明水道体系向“ 窄深型” 形态发展。同一坡度条件下, 水道体系形态特征参数之间也存在着内在关联; 一般情况下, 同一坡度的水道体系宽度与深度、深度和宽深比呈2种截然不同的关系: 在相对缓坡区, 两者呈正相关趋势, 而相对陡坡区两者呈负相关关系。至于水道体系的宽度与宽深比, 其在缓坡区和陡坡区均呈正比关系, 即宽度越大, 宽深比越大。

作者声明没有竞争性利益冲突.

参考文献
1 陈宇航, 姚根顺, 福亮, 鲁银涛, 陈亮, 唐鹏程, 曹全斌. 2017. 东非鲁伍马盆地渐新统深水水道—朵体沉积特征及控制因素. 石油学报, 38(9): 1047-1058.
[Chen Y H, Yao G S, Fu L, Lu Y T, Chen L, Tang P C, Cao Q B. 2017. Sedimentary characteristics and controlling factors of Oligocene deep-water channel-lobe in Rovuma Basin of East Africa. Acta Petrolei Sinica, 38(9): 1047-1058] [文内引用:1]
2 陈志鹏, 鲍志东, 任战利, 赵艳军, 于春勇. 2017. 尼日尔三角洲盆地某油田深水海底扇沉积特征. 地质科技情报, 36(3): 174-181.
[Chen Z P, Bao Z D, Ren Z L, Zhao Y J, Yu C Y. 2017. Sedimentary characteristics of a submarine fan of an oilfield in Niger Basin. Geological Science and Technology Information, 36(3): 174-181] [文内引用:1]
3 邓荣敬, 邓运华, 于水, 侯读杰. 2008. 尼日尔三角洲盆地油气地质与成藏特征. 石油勘探与开发, 35(6): 755-762.
[Deng R J, Deng Y H, Yu S, Hou D J. 2008. Hydrocarbon geology and reservoir formation characteristics of Niger Delta Basin. Petroleum Exploration and Development, 35(6): 755-762] [文内引用:2]
4 胡孝林, 于水, 刘新颖. 2014. 尼日尔三角洲盆地过路朵叶体特征及发育模式. 东北石油大学学报, 38(5): 31-40.
[Hu X L, Yu S, Liu X Y. 2014. Characteristic and development model of transient lobe in Niger Delta basin. Journal of Northeast Petroleum University, 38(5): 31-40] [文内引用:1]
5 康洪全, 孟金落, 程涛, 贾怀存, 白博, 李明刚. 2013. 巴西坎波斯盆地深水沉积体系特征. 石油勘探与开发, 45(2): 1-11.
[Kang H Q, Meng J L, Cheng T, Jia H C, Bai B, Li M G. 2013. Characteristics of deep water depositional system in Campos basin, Brazil. Petroleum Exploration and Development, 45(2): 1-11] [文内引用:1]
6 廖计华, 徐强, 陈莹, 王颖, 蔡露露, 邹梦君, 曾清波, 焦振华. 2016. 白云—荔湾凹陷珠江组大型深水水道体系沉积特征及成因机制. 地球科学, 41(6): 1041-1054.
[Liao J H, Xu Q, Chen Y, Wang Y, Cai L L, Zou M J, Zeng Q B, Jiao Z H. 2016. Sedimentary characteristics and genesis of the deepwater channel system in Zhujiang Formation of Baiyun-Liwan Sag. Earth Science, 41(6): 1041-1054] [文内引用:1]
7 吕明, 王颖, 陈莹. 2008. 尼日利亚深水区海底扇沉积模式成因探讨及勘探意义. 中国海上油气, 20(4): 275-282.
[ M, Wang Y, Chen Y. 2008. A discussion on origins of submarine fan deposition model and its exploration significance in Nigeria deep-water area. China Offshore Oil and Gas, 20(4): 275-282] [文内引用:1]
8 林煜, 吴胜和, 王星, 路瑶, 万琼华, 张佳佳, 张义楷. 2013. 深水浊积水道体系构型模式研究: 以西非尼日尔三角洲盆地某深水研究区为例. 地质论评, 59(3): 510-520.
[Lin Y, Wu S H, Wang X, Lu Y, Wan Q H, Zhang J J, Zhang Y K. 2013. Research on architecture model of deepwater turbidity channel system: A case study of a deepwater research area in Niger Delta Basin, West Africa. Geological Review, 59(3): 510-520] [文内引用:1]
9 刘峰, 裴健翔, 汪洋, 高华, 潘光超, 李洋森. 2015. 古地貌对海底扇沉积过程的控制及与油气富集的关系: 以莺歌海盆地东方区黄流组一段为例. 中国海上油气, 27(4): 39-46.
[Liu F, Pei J X, Wang Y, Gao H, Pan G C, Li Y S. 2015. Palaeomorphologic control on sedimentary process of submarine fans and hydrocarbon accumulation: A case study of Member 1 of Huangliu Formation in DF area, Yinggehai basin. China Offshore Oil and Gas, 27(4): 39-46] [文内引用:1]
10 王琪, 田兵, 马晓峰, 牟炜卫, 高丽华. 2017. 珠江口盆地白云深水区深水水道沉积体系及成因模式. 天然气地质学, 28(10): 1497-1505.
[Wang Q, Tian B, Ma X F, Mou W W, Gao L H. 2017. Deposition system and formation mechanism of deepwater channel in Baiyun deepwater area, Pearl River Mouth Basin. Natural Gas Geoscience, 28(10): 1497-1505] [文内引用:1]
11 张兴. 2000. 尼日尔三角洲盆地油气地质特征及潜力分析. 石油学报, 21(1): 40-44.
[Zhang X. 2000. Petroleum geological characteristics and potential analysis in the Niger Delta Basin. Acta Petrolei Sinica, 21(1): 40-44] [文内引用:1]
12 赵晓明, 吴胜和, 岳大力, 刘丽. 2010. 深水海底扇储集层构型浅析: 以西非某油田为例. 见: 第五届油气成藏机理与油气资源评价国际学术研讨会论文集. 北京: 石油工业出版社, 722-728.
[Zhao X M, Wu S H, Yue D L, Liu L. 2010. Simple analysis of the deep-water fan reservoir: A case study from the oilfield in West Africa. In: 5th Internal Symposium on Hydrocarbon Accumulation Mechanisms and Petroleum Resource Evolution. Beijing: Petroleum Industry Press, 722-728] [文内引用:1]
13 赵晓明, 吴胜和, 刘丽. 2012 a. 西非陆坡区深水复合水道沉积构型模式. 中国石油大学学报(自然科学版), 36(6): 1-5.
[Zhao X M, Wu S H, Liu L. 2012 a. Sedimentary architecture model of deep-water channel complexes in slope area of West Africa. Journal of China University of Petroleum(Natural Edition), 36(6): 1-5] [文内引用:1]
14 赵晓明, 吴胜和, 刘丽. 2012 b. 尼日尔三角洲盆地AKPO油田新近系深水浊积水道储集层构型表征. 石油学报, 33(6): 1049-1058.
[Zhao X M, Wu S H, Liu L. 2012 b. Characterization of reservoir architectures for Neogene deepwater turbidity channels of Akpo oilfield, Niger Delta Basin. Acta Petrolei Sinica, 33(6): 1049-1058] [文内引用:1]
15 Alpak F O, Barton M D, Naruk S J. 2013. The impact of fine-scale turbidite channel architecture on deep-water reservoir performance. AAPG Bulletin, 97: 251-284. [文内引用:1]
16 Clark I R, Cartwright J A. 2009. Interactions between submarine channel systems and deformation in deepwater fold belts: Examples from the levant basin, eastern mediterranean sea. Marine & Petroleum Geology, 26(8): 1465-1482. [文内引用:1]
17 Cohen H A, Mcclay K. 1996. Sedimentation and shale tectonics of the northwestern Niger delta front. Marine & Petroleum Geology, 13: 313-328. [文内引用:1]
18 Corredor F, Shaw J H, Bilotti F. 2005. Structural styles in the deep-water fold and thrust belts of the Niger Delta. AAPG Bulletin, 89: 753-780. [文内引用:1]
19 Damuth J E. 1994. Neogene gravity tectonics and depositional processes on the deep Niger Delta continental margin. Marine & Petroleum Geology, 11: 320-346. [文内引用:1]
20 Doust H, Omatsola E. 1989. Niger delta. In: Edwards J D, Santogrossi P A(eds). Divergent/Passive Margin Basins. American Association of Petroleum Geologists, Tulsa, 201-238. [文内引用:3]
21 Flood R D, Piper D J W. 1997. Amazon Fan sedimentation: The relationship to equatorial climate change, continental denudation, and sea-level fluctuations. Proceedings, 155. [文内引用:1]
22 Gong C, Steel R J, Wang Y, Lin C, Olariu C. 2016 a. Grain size and transport regime at shelf edge as fundamental controls on delivery of shelf-edge sand s to deepwater. Earth-Science Reviews, 157: 32-60. [文内引用:2]
23 Gong C, Steel R J, Wang Y, Lin C, Olariu C. 2016b. Shelf-margin architecture variability and its role in sediment-budget partitioning into deep-water areas. Earth-Science Reviews, 154: 72-101. [文内引用:2]
24 Khripounoff A, Wagonergriesheim A, Babonneau N, Crassous P, Dennielou B, Savoye B. 2003. Direct observation of intense turbidity current activity in the Zaire submarine valley at 4000 m water depth. Marine Geology, 194: 151-158. [文内引用:1]
25 Kneller B. 2003. The influence of flow parameters on turbidite slope channel architecture. Marine & Petroleum Geology, 20: 901-910. [文内引用:1]
26 Kolla V, Bourges P, Urruty J M, Safa P. 2001. Evolution of Deep-Water Tertiary Sinuous Channels Offshore Angola(West Africa)and Implications for Reservoir Architecture. AAPG Bulletin, 85: 1373-1405. [文内引用:1]
27 Kolla V. 2007. A review of sinuous channel avulsion patterns in some major deep-sea fans and factors controlling them. Marine & Petroleum Geology, 24(6-9): 450-469. [文内引用:2]
28 Labourdette R, Bez M. 2010. Element migration in turbidite systems: Rand om or systematic depositional processes?AAPG Bulletin, 94: 345-368. [文内引用:1]
29 Lamb M A. 2003. Stratigraphic architecture of a sand -rich, deep-sea depositional system: The stevens sand stone, San Joaquin Basin. California. Oklahoma, AAPG Special Publication. [文内引用:1]
30 Liu L, Zhang T, Zhao X, Wu S, Hu J, Wang X, Zhang Y. 2013. Sedimentary architecture models of deepwater turbidite channel systems in the Niger Delta continental slope, West Africa. Petroleum Science, 10(2): 139-148. [文内引用:1]
31 Mayall M, Jones E, Casey M. 2006. Turbidite channel reservoirs: Key elements in facies prediction and effective development. Marine & Petroleum Geology, 23: 821-841. [文内引用:2]
32 Mchargue T, Pyrcz M J, Sullivan M D, Clark J D, Fildani A, Romans B W, Covault J A, Levy M, Posamentier H W, Drinkwater N J. 2011. Architecture of turbidite channel systems on the continental slope: Patterns and predictions. Marine & Petroleum Geology, 28: 728-743. [文内引用:1]
33 Morley C K, Guerin G. 1996. Comparison of gravity-driven deformation styles and behavior associated with mobile shales and salt. Tectonics, 15: 1154-1170. [文内引用:2]
34 Mutti E, Normark W. 1987. Comparing examples of modern and ancient turbidite systems: Problems and concepts. In: Leggett J K, Zuffa G G(eds). Marine Clastic sedimentology: Concepts and Case Studies. Springer. Netherland s, 1-38. [文内引用:1]
35 Piper D J W, Normark W R. 1983. Turbidite depositional patterns and flow characteristics, Navy Submarine Fan, California Borderland . Sedimentology, 30: 681-694. [文内引用:2]
36 Short K C, Stäuble A J. 1967. Outline of geology of Niger delta. AAPG Bulletin, 51: 761-779. [文内引用:1]
37 Stow D A V. 1981. Laurentian fan: Morphology, sediments, processes, and growth pattern. AAPG Bulletin, 65(3): 375-393. [文内引用:2]
38 Wood L J, Mize-Spansky K L. 2009. Quantitative seismic geomorphology of a Quaternary leveed-channel system, offshore eastern Trinidad and Tobago, northeastern South America. AAPG Bulletin, 93: 101-125. [文内引用:1]
39 Wynn R B, Cronin B T, Peakall J. 2007. Sinuous deep-water channels: Genesis, geometry and architecture. Marine & Petroleum Geology, 24: 341-387. [文内引用:1]
40 Zhang J J, Wu S H, Fan T E, Fan H J, Jiang L, Chen C, Wu Q Y, Lin P. 2016. Research on the architecture of submarine-fan lobes in the Niger Delta Basin, offshore West Africa. Journal of Palaeogeography, 5(3): 185-204. [文内引用:1]
41 Zhang J, Wu S, Hu G, Fan T E, Yu B, Lin P, Jiang S. 2018. Sea-level control on the submarine fan architecture in a deepwater sequence of the Niger Delta basin. Marine & Petroleum Geology, 94: 179-197. [文内引用:1]
42 Zhao X, Qi K, Liu L, Gong C, Mccaffrey W D. 2018a. Development of a partially-avulsed submarine channel on the Niger Delta continental slope: Architecture and controlling factors. Marine & Petroleum Geology, 95: 30-49. [文内引用:1]
43 Zhao X, Qi K, Liu L, Xie T, Li M, Hu G. 2018b. Quantitative characterization and controlling factor analysis of the morphology of bukuma-minor channel on southern Niger Delta slope. Interpretation, 6: 57-69. [文内引用:1]