主动、被动大陆边缘天然气水合物成藏模式对比
胡高伟1,2, 卜庆涛1,2, 吕万军3, 王家生3, 陈杰1, 李清1,2, 龚建明1,2, 孙建业1,2, 吴能友1,2
1. 自然资源部天然气水合物重点实验室·中国地质调查局青岛海洋地质研究所;
2. 青岛海洋科学与技术试点国家实验室海洋矿产资源评价与探测技术功能实验室 3.中国地质大学(武汉)海洋学院
通信作者:吴能友,1965 年生,研究员,博士研究生导师,本刊第八届编委会委员;主要从事天然气水合物方面的研究工作。地址:(266071)山东省青岛市市南区福州南路62号。ORCID: 0000-0002-1734-6249。E-mail: wuny@ms.giec.ac.cn

作者简介:胡高伟,1982年生,副研究员,博士;主要从事海洋地质学与天然气水合物方面的研究工作。地址:(266071)山东省青岛市市南区福州南路62号。ORCID: 0000-0002-1575-074X。E-mail: hgw-623@163.com

摘要

不同大陆边缘(主动、被动)沉积物中天然气水合物(以下简称水合物)赋存的控制因素与成藏模式有所差异,开展二者之间的对比研究对于指导水合物勘探具有重要的意义。为此,以主动大陆边缘卡斯凯迪亚(Cascadia)和日本南海(Nankai)海槽、被动大陆边缘布莱克海台(Blake Ridge)和尼日尔三角洲盆地(Niger Delta Basin)等典型水合物成藏区为研究对象,借助于综合大洋钻探(IODP)航次资料解剖和数值模拟分析等手段,从应力场的角度探讨了上述两种背景下含甲烷流体的驱动样式,进而对比分析了主动、被动大陆边缘水合物的成藏模式。研究结果表明:①主动大陆边缘以侧向挤压应力为气体垂向运移提供了驱动力和通道,引诱深部游离气和原位生物气沿断裂运移,气体运移通道主要为俯冲—增生产生的断层、断裂和滑塌体;②主动大陆边缘粉砂和砂质粉砂等粗粒浊流沉积孔隙度大、渗透性好,并且增生楔上沉积物厚度大,是水合物成藏较为有利的储集空间;③较之于主动型大陆边缘,被动大陆边缘虽然缺少俯冲带造成的侧向应力,但在其内巨厚沉积层塑性物质及高压流体、陆缘外侧火山活动等的共同作用下,产生垂向加积和拉张应力,形成的扩散型水合物聚集速率主要取决于甲烷的供给速度;④被动大陆边缘有机质含量、产气速率、地温梯度及沉积速率对水合物含量空间分布具有差异性影响,泥火山或底辟构造等为水合物的形成与赋存提供了理想的场所。

关键词: 主动大陆边缘; 被动大陆边缘; 天然气水合物; 成藏模式对比; 卡斯凯迪亚边缘; 日本南海海槽; 布莱克海台; 尼日尔三角洲盆地; 聚集速率
A comparative study on natural gas hydrate accumulation models at active and passive continental margins
HU Gaowei1,2, BU Qingtao1,2, LYU Wanjun3, WANG Jiasheng3, CHEN Jie1, LI Qing1, GONG Jianming1,2, SUN Jianye1,2, WU Nengyou1,2
1. Key Laboratory of Gas Hydrate, Ministry of Natural Resources//Qingdao Institute of Marine Geology, China Geological survey, Qingdao, Shandong 266071, China
2. Laboratory for Marine Mineral Resources, Pilot National Laboratory for Marine Science and Technology-Qingdao, Qingdao, Shandong 266071, China
3. College of Marine Science and Technology, China University of Geosciences-Wuhan, Wuhan, Hubei 430074, China
Abstract

The comparative study on natural gas hydrate accumulation models between active and passive continental margins as well as their controlling factors is of great significance to the guidance of natural gas hydrate exploration. Based on the data and research results of international typical active continental margin hydrate accumulation areas such as the Cascadia margin of the Northeast Pacific, the Nankai trough, etc. and passive continental margin areas like the Blake Ridge, the models of the gas hydrate accumulation system are summarized and numerically simulated, and a preliminary comparison of active and passive continental margin reservoir accumulation models was also carried out. The following results were obtained. (1) The active continental margin provides a driving force and channel for vertical gas migration, which induces deep free gas and in-situ biogas to migrate along the fault. The migration channels are mainly faults, fractures and slumps produced by subduction-accretion. (2) Coarse-grained turbidity sediments such as silt and sandy silt have good porosity and permeability. Moreover, the sediment thickness on the accretionary wedge is large, which provides a good storage space for hydrate accumulation. (3) Numerical simulations of the Blake Ridge, and Niger Delta hydrate accumulation show that the passive continental margin lacks the lateral stress caused by the subduction zone compared with the active continental margin. However, due to the plastic materials in the thick sedimentary layer, high-pressure fluids and volcanic activities outside the continental margin, vertical accretion and tensile stress are generated and the accumulation rate of diffusion-type hydrates mainly depends on the methane supply rate. (4) Organic matter content, gas production rate, geothermal gradient and sedimentation rate at the passive continental margin have different effects on the spatial distribution of hydrate content. Mud volcanoes or diapir structures provide an ideal place for the formation and occurrence of hydrates.

Keyword: Active Continental Margin; Passive Continental Margin; Natural gas hydrate; Accumulation Model; Cascadia margin; Nankai trough; Blake Ridge; Niger Delta Basin; Accumulation rate
1 主动/被动大陆边缘天然气水合物成藏研究概况

天然气水合物(以下简称水合物)广泛分布于世界大陆边缘海底和陆地永久冻土带中[1], 大陆边缘介于海岸线和深海洋盆之间, 包括大陆架、大陆坡和大陆隆(或海沟)等海底地貌单元, 可分为主动大陆边缘(也称太平洋型大陆边缘、汇聚大陆边缘、活动大陆边缘等)和被动大陆边缘(也称大西洋型大陆边缘、稳定大陆边缘、非活动大陆边缘等)[2, 3, 4]。不同大陆边缘沉积物中水合物赋存的控制因素和成藏模式有所差异[5], 开展两者之间的对比研究对于理解水合物富集规律及指导水合物勘探开发具有重要的意义。

国内外学者对于主动/被动大陆边缘水合物发育情况已有较多的研究[5, 6, 7, 8, 9, 10, 11, 12], 其对比研究大致可分为点对点对比研究阶段(20世纪60年代— 80年代末, 以单孔钻探点对点对比为主, 偏重于揭示两种不同地质背景下水合物成藏的共同特征[2, 3, 10, 11])、有利成藏构造单元对比研究阶段(20世纪90年代— 21世纪初, 逐渐形成了水合物成藏系统观[13, 14, 15, 16, 17])和成藏系统研究与成藏模式对比阶段(21世纪初至今, 对典型区域的水合物成藏系统进行了深入解剖, 逐步形成成藏模式规律性认识[8, 18, 19, 20])。

水合物成藏系统研究覆盖气体生成、运移、最终聚集成藏的整个过程, 具体集中在气源、输导体系、有利储集空间和高压低温稳定带等四大要素方面。气体来源及类型主要为生物成因气、热解成因气和混合气[21]。最新的研究成果认为:蛇纹岩化生烃作用可提供无机成因气[22]; 气体运移通道主要包括底辟构造、断层、裂隙与高渗透层等, 运移通道的展布、规模及运移效能, 对水合物成藏富集和饱和度大小具有重要的影响[23]; 温度、压力、孔隙水盐度和气体组分决定了水合物稳定带的厚度, 而良好的沉积储集空间是水合物能够富集成藏的重要保证[20]; 依据储集空间的不同主要将水合物藏分为孔隙型和裂隙型[24, 25], 依据甲烷气体运移方式与通量大小可分为低通量下的扩散型(Distributed Low-Flux Gas Hydrate)和高通量下的渗漏型(Focused High-Flux Gas Hydrate)[26]。水合物聚集成藏的整个过程中, 上述要素不可或缺[27, 28]

对世界范围水合物与主动/被动大陆边缘分布范围进行投点分析和统计研究, 结果表明(图1、表1):水合物主要赋存于主动大陆边缘的增生楔和被动大陆边缘的泥火山、泥底辟、麻坑、断层、丘状体及冷泉喷口等地区。总体来说, 构造演化、全球海平面变化、沉积物供给等控制了海底沉积层中物质充填样式[38], 制约着沉积体系内物性空间展布与演化[39, 40], 沉积过程中有机质的输入控制了海底有机质和气源的空间分布[10]; 构造— 沉积体系的演化伴随着空间上温度、压力的变化, 决定了沉积物中有机质— 甲烷— 水合物体系物质相互转化的速率[41], 并控制着物质存在的形式以及沉积物孔隙流体组成与性质的变化。因此可以看出, 构造演化和沉积充填特征既决定了含甲烷流体流动的动力场, 又决定其物理化学场[8, 42], 从而在影响沉积物中水合物形成与分解的动态变化方面具有至关重要的作用[43], 是控制水合物分布及富集的主控因素。然而, 主动、被动大陆边缘因具有显著不同的应力场, 对孔隙流体、物质迁移的控制模式有着很大的差别, 其具体表现形式及对水合物成藏过程和空间分布具有何种影响, 目前尚未取得清晰的认识。

图1 研究区位置及水合物分布示意图

表1 主动/被动大陆边缘主要水合物区域分布表[2, 3, 29, 30, 31, 32, 33, 34, 35, 36, 37]

笔者选取主动大陆边缘卡斯凯迪亚(Cascadia)和日本南海(Nankai)海槽、被动大陆边缘布莱克海台(Blake Ridge)和尼日尔三角洲盆地(Niger Delta Basin)等典型水合物成藏区为研究对象, 通过IODP航次资料解剖和数值模拟分析等手段, 从应力场角度探讨两种背景下含甲烷流体流动样式的控制, 开展主动/被动大陆边缘水合物成藏模式对比研究, 以期为认识我国南海复杂边缘的水合物富集规律及成藏模式提供参考。

2 研究区概况及研究方法
2.1 研究区地质背景与水合物勘探概况

研究区如图1所示。Cascadia主动大陆边缘位于东北太平洋胡安· 德富卡(Juan de Fuca)板块向北美板块俯冲的Cascadia俯冲增生楔上[44], Juan de Fuca板块每年以约45 mm、几乎垂直的方向向北美板块俯冲, 在变形前沿向海方向, Cascadia沉积了前更新统半深海沉积, 上覆快速沉积、总厚度达2500 m的更新统浊积岩。国际大洋钻探ODP204航次和IODP311航次对该区进行了较为深入的研究[45, 46]。ODP204航次在水合物海岭南峰从峰顶到斜坡盆地共钻9个站位, 水深介于700~1 000 m, 代表了通过构造抬升形成的背斜和斜坡盆地不同构造部位[45]。IODP311航次钻探目标是检验该区水合物形成及控制模式的参数[46, 47], 在横切整个大陆边缘、水深900~2 200 m范围内, 共实施钻探了5个站位(U1325、U1326、U1327、U1328、U1329), 代表着水合物在Cascadia边界的不同演化阶段, 从最早出现于西面的增生脊到最后出现在走向东面界限的稍浅水区。

日本Nankai海槽位于菲律宾海西北日本岛弧西南部, 也是典型的主动大陆边缘。Kumano弧前盆地位于日本本州岛的Kii半岛[48, 49], 旁边发育日本Nankai增生楔(中新世— 距今15 Ma), 菲律宾板块以平均每年40.0~58.4± 1.2 mm的速度和300° ~315° 方向向欧亚板块俯冲。Kumano盆地是日本Nankai海槽最大的弧前盆地, 东西延伸100 km, 南北延伸70 km, 水深2 000 m以下, 底部相对平坦并充满海底扇沉积体系, 沉积体系而后被断层切断[50]。由角度不整合可以将盆地分为4个基本单元, 盆地内沉积物由于盆地向海边缘的抬升而下沉, 沉积中心随之变化。在日本Nankai海槽Kumano近海和Kumano盆地向陆斜坡的下方观察到平行于海底的BSR。

Blake Ridge东南延伸方向与北美大陆边缘成正交, 是典型的被动大陆边缘。ODP164航次在该区最早组织了水合物航次调查, 共在3个不同区域钻获岩心样品:布莱克海脊、开普菲尔底辟以及布莱克海脊底辟。其中991、992和993站位位于开普菲尔底辟之上或附近; 994、995和997站位位于布莱克海脊的东部边缘之上; 996站位则位于布莱克海脊之上[51]

Niger Delta Basin是早白垩世开始发育的被动大陆边缘盆地, 包括裂谷期和漂移期两个演化阶段, 盆地形成发育与冈瓦纳大陆裂解和南大西洋、赤道大西洋张开有关。始新世以来长期海退形成了现今的Niger Delta Basin。

2.2 研究方法

采用资料分析和数值模拟方法, 分别以Cascadia大陆边缘、日本Nankai海槽和Blake Ridge、Niger Delta Basin为主动/被动大陆边缘实例开展应力场对水合物成藏影响研究对比。主动、被动大陆边缘所具有的构造挤压背景和正常沉积压实背景具有不同的应力场, 对沉积层孔隙度、渗透性的演化起着重要作用, 进而可能影响流体流动样式和含甲烷流体的运移与聚集形式。海底沉积物中甲烷从生成到汇聚在稳定域内部固结为水合物, 流体的流动是输送水合物聚集所需甲烷的关键过程。流体流动的驱动因素和通道受构造应力场的控制。因此, 面向水合物成藏机理与分布规律的数值模拟, 必须将宏观地质演化模拟与微观有机质— 甲烷— 水合物物质能量演化相结合, 才有可能从本质上揭示水合物成藏的实际过程。

宏观地质演化模拟, 基于介质连续性方程、流体流动方程、热流方程、有效应力作用及渗流方程[42, 52, 53, 54, 55, 56, 57], 模拟沉积压实及应力场、温度场、压力场、流体流动。微观物质迁移— 转化模拟, 建立沉积单元体孔隙流体中溶解甲烷浓度的物质守恒模型, 即甲烷浓度随时间上的变化量与以下物理过程的甲烷浓度变化量之间守恒:①空间上的扩散、对流、有机质热解和乙酸发酵; ②CO2还原作用生成的量、被硫酸盐氧化的量; ③与水合物和游离气的转化量。从自由能的角度研究水合物成核的驱动力和物质在空间迁移(扩散)的驱动力, 建立了反映微观孔隙尺度水合物形成演化机理的中热质空间迁移— 转化模型, 从而从机理上反映温度梯度、盐度梯度、浓度梯度影响下物质的迁移趋势和规律, 将这些微观物质— 能量演化过程与宏观的地质演化过程通过温度和压力进行耦合。基于这些模型开发的水合物成藏机理数值模拟分析系统, 宏观过程与微观过程采用分步迭代法(SIP)计算, 偏微分方程组采用有限元求解。有关模型其他相关参数的计算参见本文参考文献[57]。

3 不同大陆边缘水合物成藏模式
3.1 应力作用下主动大陆边缘水合物成藏模式

3.1.1 Cascadia大陆边缘水合物成藏特征

IODP311航次5个钻探站位采集的水合物样品代表了水合物形成的3个时期, 即早期(U1326、U1325站位, 靠近俯冲带一侧, 沉积物年龄不超过1 Ma)、中期(U1327、U1328站位, 沉积物年龄不超过2 Ma)和消亡期(U1329站位, 靠近陆地最浅水区, 沉积物年龄介于2~9 Ma)[58]。从5个站位的C1/C2随深度变化和横跨大陆边缘的变化特征来看, 除活动喷口U1328站位的甲烷浓度明显偏高外, 不同成藏阶段气源无明显变化, 碳同位素值显示气源主要为生物成因气[46]。而不同站位的沉积物粒度特征表明, U1326— U1329沉积环境从远洋单一物源的细粒沉积向近远洋、近陆的复杂物源细— 中粒沉积环境转变, 电阻率显示水合物饱和度存在明显差异(5%~40%), 横向上U1327站位饱和度最大[46], 可见沉积物粒度对饱和度具有一定的控制作用[47]

BSR位置(可视为水合物带底界)随着成藏早期、中期、晚期逐渐变浅(图2), 而甲烷— 硫酸盐界面(SMI, 可视为水合物带顶界)逐渐变深, 表明水合物稳定带厚度呈逐渐减薄、地温梯度呈逐渐变大的趋势。应力作用下热流变化和物质差异分布, 使不同站位BSR稳定程度各异(图2), 其中U1326站位BSR沿脊顶延伸稳定、侧向不稳定, U1325站位BSR在盆地东侧延伸稳定, U1327、U1329站位BSR表现清晰且延伸相对稳定, 而U1328站位BSR具有大量盲带、指示活动喷流状态。较高的水合物饱和度和有限的原位有机碳供给, 指示水合物气源一部分源自原位、一部分源自水合物稳定带之下。增生楔形成过程产生的断裂、褶皱、底辟等构造组合, 为流体流动和水合物成藏提供了良好条件。

图2 IODP311站位的BSR总体表现和横向变化特征图[46]

综合上述特征, 笔者修正了前人基于流体流动的水合物成藏模型(图3), 认为在Cascadia边缘, 沉积有机质通过远洋沉积输入或浊流沉积输入之后进行沉积和埋藏, 甲烷通过原位产生和由下部通道运移而来。富甲烷的孔隙流体和游离甲烷气沿着一系列断裂/断层向上运移, 由于海底附近甲烷溶解度大大降低, 大量气体生成水合物。在Juan de Fuca板块俯冲的侧向应力作用下, 增生楔形成过程中沉积埋葬和横向运输, 随着变形前端距离和深度变化而富集不同浓度的甲烷, 为水合物形成富集提供了气源。水合物稳定带底界(BGHSZ)会形成“ 盖层” , 对下部游离气产生暂时阻碍作用, 但由于构造抬升等作用, 会使底界水合物分解而发生再循环。

图3 Cascadia大陆边缘水合物成藏模型图
(修编自本文参考文献[18])

3.1.2 Kumano盆地水合物成藏特征

主动大陆边缘有效应力是沉积负载压力、侧向构造挤压应力、流体压力等变量的函数[58, 59], 笔者以日本Nankai海槽Kumano盆地为例开展了数值模拟研究。日本Nankai海槽增生楔的气体不仅具有热成因特征, 而且Muroto近岸具有热解与生物混合成因特征、向陆斜坡和增生楔弧前盆地具有生物成因特征[59], TOC(总有机碳含量)约为0.4%。根据孔隙水中碘离子(I-)浓度指示, I-浓度在海槽向陆一侧随着深度的增加而迅速增大, 比海水高数个数量级, 而向海一侧仅略高于海水, 甲烷浓度具有类似规律[60]。同时, 向海一侧沉积物孔隙水中碘年龄不足7 Ma, 而向陆一侧样品孔隙水中碘的年龄超过30 Ma。浓度分布和年龄表明, 碘连同甲烷一起从上部板块老的地层中输送到向陆一侧的流体内, 特定的流体通道可能由楔状体、滑脱构造或断层等组成。

利用文献中有机碳含量、甲烷生成条件、地震剖面和地质格架等资料开展数值模拟, 结果表明:在应力作用下, 构造变形比较强烈, 流体活动、甲烷水合物聚集空间上显著受构造断裂裂隙所控制, 水合物饱和度最高可达约30%。模拟结果与Kumano盆地地震剖面和解释结果基本一致, 表明模拟结果具有可靠性。侧向应力导致构造变形及水合物运聚的变化规律初步揭示了Kumano盆地水合物的成藏过程(图4、5)。

图4 Kumano盆地3个演化阶段的有效应力与垂向流速分布图
注:mbsl是meter below sea level的缩写, 表示海平面以下深度, 单位为m; 1 atm=0.101 325 MPa, 下同。

图5 Kumano盆地3个演化阶段的水合物含量分布图

Cascadia和Kumano盆地水合物成藏研究结果表明, 流体的流动是输送水合物聚集所需甲烷的关键过程, 其驱动因素和流体的通道受构造应力场的控制。主动大陆边缘侧向构造应力驱动非常突出, 控制了含甲烷流体流动和演化以及水合物的聚集分布。同时, 构造挤压作用导致热流通量、流体通量、甲烷通量横向上的显著差异, 也是造成水合物空间赋存分布差异的重要因素之一。

3.2 被动大陆边缘水合物成藏模式

3.2.1 Blake Ridge水合物成藏特征

根据994站位、995站位和997站位的地层厚度及岩性参数(表2), 建立该区典型的二维地质模型, 具体数据包括:994站位, 晚第四纪(0~距今0.46 Ma)、中第四纪(距今0.46~0.98 Ma)和早第四纪(距今0.98~1.80 Ma)沉积速率分别为30 m/Ma、55 m/Ma和98 m/Ma, 平均沉积速率68 m/Ma; 晚上新世沉积速率增加到89 m/Ma, 早上新世沉积速率继续增加至140 m/Ma, 晚中新世沉积速率高达303 m/Ma。995站位, 晚第四纪、中第四纪和早第四纪的沉积速率分别为29 m/Ma、60 m/Ma、80 m/Ma, 平均沉积速率60 m/Ma。晚上新世沉积速率为104 m/Ma。地温梯度的选取根据各站位实际的钻探资料, 取值为35 K/km。

表2 Blake Ridge各站位地层厚度及岩性参数表 单位:m

研究结果表明(图6), 沉积物中有机质丰度越高, 水合物含量越高, 水合物在空间的分布范围越广。有机质含量较低时, 孔隙流体中溶解甲烷的浓度不足以形成高饱和度、厚层的水合物, 水合物层仅会在稳定域底界附近形成, 厚度很薄。沉积速率对水合物的含量分布也会产生影响, 在相对较高的沉积速率条件下, 有利于水合物的形成和分布; 而当沉积速率较低时, 沉积物堆积较慢, 压实过程中孔隙流体排出速度较慢, 形成水合物层的厚度很薄。随着沉积速率的加快, 沉积物压实过程加快, 含甲烷的流体被快速排出, 甲烷供给量大, 水合物层显著增厚。地温梯度越低, 水合物在垂向上形成的位置越深, 此外, 较低的地温梯度有利于水合物的形成及空间展布。

图6 泥岩有机质丰度的差异对水合物空间含量的影响图

根据有效应力原理, 在外力作用下沉积地层应力被沉积物骨架和孔隙中流体共同承担, 但是只有通过沉积物颗粒传递的有效应力(总应力与流体压力的差)才会使沉积物产生变形从而影响沉积物孔隙度及渗透率。对于Blake Ridge, 沉积地层应力主要来自其上覆沉积物和海水的总荷载, 在流体没有超压的情况下, 有效应力等于总荷载应力与静水压力的差, 因此模拟剖面上(图7)有效应力随深度均匀增加而等值线大致与海底地形平行, 在这样的有效应力场下, 沉积层随着时间的推移逐渐被压实, 孔隙度则从海底往下呈近似指数形式变小, 渗透系数的对数与孔隙度近似成正比变化。由于剖面左侧(994站位)比右侧(997站位)泥岩含量高, 因而渗透率比右侧要低, 但第四纪以来沉积速率比右侧更大。因此, 单位时间内压实速率比右侧高, 从而排出流体的速率更大, 富水合物层左侧比右侧薄, 与实际调查结果相吻合。

图7 Blake Ridge有效应力、孔隙度、渗透率、超压及水合物含量分布图

3.2.2 Niger Delta Basin水合物成藏特征

Niger Delta Basin水合物成藏数值模拟研究结果表明(图8), 在300~350 km附近一直是沉降中心, 沉积物从右侧向左侧输运形成向海方向进积沉积层序。在这种沉积充填格架下, 盆地中心巨厚的烃源岩不断生烃, 一定条件下产生超压。流体的压力大于静水压力, 造成同深度下有效应力相对其他地段要小, 从而沉积物相对欠压实而孔隙度比其他地段略高。受岩性和孔隙度空间差异的影响, 渗透系数横向和纵向上也呈不均一变化。从烃源岩中排出的气体随流体主要由沉积中心向两侧侧上方向运移, 右侧地形高而不具备水合物稳定条件, 水合物主要富集于水深600 m以深且在稳定域中离盆地生烃中心较近的地段较为富集。

图8 Niger Delta Basin孔隙度、渗透率、超压、流体压力梯度、有效应力及水合物含量分布图

从Blake Ridge和Niger Delta Basin等分析结果可以看出, 被动大陆边缘是正常沉积压实应力场下的水合物聚集, 水合物成藏分布主要受沉积速率、地温梯度与有机质含量等因素控制, 沉积压实、浮力、重力等是驱动流体流动的主要因素, 因而高有机质含量、高沉积速率和低地温梯度有利于水合物的成藏富集。

4 不同大陆边缘水合物成藏影响对比

综合文献资料和本文相关结果, 初步获取了主动/被动大陆边缘背景下水合物的成藏特征(表3)。

表3 主动大陆边缘与被动大陆边缘成藏特点对比表

1)两者水合物藏的气体成因类型无明显区别, 均为生物成因气、热成因气或混合成因气, 因蛇纹岩化生烃具有形成无机成因甲烷的潜景, 可能成为主动大陆边缘的气源新成因类型[22]

2)应力场作用下, 流体流动模式的不同可能是造成主动/被动大陆边缘水合物成藏差异的主要因素之一。侧向应力作用对主动大陆边缘增生楔具有以下影响:①构造变形比较强烈, 楔状体和断层等形成流体通道, 使甲烷发生向陆迁移, 并导致热流通量、流体通量、甲烷通量具有显著的横向差异; ②浊流沉积发育, 沉积物粒度相对较高, 更易形成高饱和度的水合物储层。被动大陆边缘受沉积压实、浮力、重力等作用, 快速沉积压实、超压等可导致甲烷流体富集, 增加水合物层厚度; 沉积物中流体通量和水合物的成藏机会更多地取决于沉积地层的有机质含量、埋藏速度和地热梯度等地质及地球物理因素, 垂向埋藏压力和孔隙流体浮力是水合物成藏的主控因素。

3)两者构造环境有所不同。主动大陆边缘背景下, 板块俯冲增生过程中构造活动发育, 从驱动力和运移通道为泥火山的形成奠定了重要基础, 因而主动大陆边缘发育泥火山及气烟囱等特殊地质体, 也为水合物成藏提供了重要的赋存地质条件。被动陆缘背景断层构造活动发育稍差于主动大陆边缘, 流体超压驱动的富烃类流体往往在地层中形成泥底辟等未能有效刺穿地层的地质体, 且由于盐底辟的作用也可产生泥火山, 泥火山及泥底辟的两翼地层及其顶部亦可以成为水合物的有效赋存层位[37, 61]

4)两者水合物成藏所需的温压场控制因素相同, 受海平面变化、海底温度和地温梯度等控制[20]

综合上述认识, 对主动/被动大陆边缘水合物成藏的宏观特征总结如下(图9)。

图9 主动、被动大陆边缘水合物成藏模式对比分析图

1)主动大陆边缘:①以侧向挤压应力为背景, 断层以逆冲为主, 可引诱深部油气藏、水合物稳定带之下游离气和原位生物气等沿断裂运移, 为水合物形成提供主要气源, 其成因类型为微生物成因气, 深部热解气也具有一定贡献; ②气体运移通道主要为俯冲— 增生产生的断层、断裂和滑塌体, 以及泥火山及滑塌构造等, 储集层主要为粉砂和砂质粉砂等粗粒浊流沉积, 具有良好的孔隙度和渗透性; ③增生楔上沉积物厚度大, 断层和褶皱发育, 有利于流体的运移、聚集, 形成水合物堆积体。

2)被动大陆边缘:①气源同样为微生物成因气和热解气; ②以垂向加积为主要应力背景, 与主动型大陆边缘相比, 由于缺少俯冲带造成的侧向应力, 容易导致部分断层和裂隙等疏导体系的缺失; ③但因其内巨厚沉积层塑性物质及高压流体、陆缘外侧火山活动及张裂作用, 可形成大规模的泥火山或底辟构造, 这些构造能使构造侧翼或顶部的沉积层倾斜, 有利于流体排放, 为水合物形成和赋存提供理想场所[37]

5 结论

1)两种地质背景下气源成因无明显区别, 均具生物成因气、热解气和混合气, 水合物形成的温度、压力控制因素也基本相同。

2)主动大陆边缘较为强烈的构造变形形成楔状体、断层等流体通道, 导致热流通量、流体通量、甲烷通量具有显著的横向差异, 对水合物成藏分布造成影响, 而被动大陆边缘缺乏规模较大的断层和裂隙等疏导体系, 有机质含量、地温梯度及沉积速率对水合物含量空间分布具有差异性影响。

3)不同的构造— 沉积体制控制下有效应力对含气流体具有的不同驱动样式, 这可能是造成主动/被动大陆边缘水合物成藏差异的重要因素。主动大陆边缘以侧向挤压应力为背景, 断层以逆冲为主, 可引诱深部油气藏、水合物稳定带之下游离气和原位生物气等沿断裂运移; 被动大陆边缘以垂向加积和拉张应力为应力背景, 断层以阶梯状正断层为主, 可诱导底辟和泥火山等作用为天然气水合物形成和赋存提供理想场所。

4)我国南海水合物成藏区从东至西具有主动— 被动大陆边缘的复杂特性, 希望通过探讨两种不同地质背景下水合物成藏的主控因素和成因机制差异, 以及水合物成藏相关的演化特征, 能够该区水合物成藏研究提供一定的指导作用。

致谢:感谢李彦龙、万义钊和李昂博士在本文撰写过程提出的宝贵建议, 感谢审稿专家和编辑部的建设性意见与建议。

参考文献
[1] 吴能友, 黄丽, 胡高伟, . 海域天然气水合物开采的地质控制因素和科学挑战[J]. 海洋地质与第四纪地质, 2017, 37(5): 1-11.
WU Nengyou, HUANG Li, HU Gaowei, et al. Geological controlling factors and scientific challenges for offshore gas hydrate exploitation[J]. Marine Geology & Quaternary Geology, 2017, 37(5): 1-11. [本文引用:1]
[2] 张光学, 黄永样, 祝有海, . 活动大陆边缘水合物分布规律及成藏过程[J]. 海洋地质动态, 2001, 17(7): 3-7.
ZHANG Guangxue, HUANG Yongyang, ZHU Youhai, et al. Distribution law and accumulation process of active continental margin hydrate[J]. Marine Geology Development, 2001, 17(7): 3-7. [本文引用:2]
[3] 张光学, 祝有海, 徐华宁. 非活动大陆边缘的天然气水合物及其成藏过程述评[J]. 地质论评, 2003, 49(2): 181-186.
ZHANG Guangxue, ZHU Youhai, XU Huaning. Gas hydrate on the passive continental margin and its pool-formation process[J]. Geological Review, 2003, 49(2): 181-186. [本文引用:2]
[4] 黎明碧, 李家彪, 方银霞. . 大陆边缘地质特征与200海里以外大陆架界限确定[M]. 北京: 海洋出版社, 2015.
LI Mingbi, LI Jiabiao, FANG Yinxia, et al. The geological characteristics of the continental margin and the determination of the continental shelf beyond 200 nautical miles[M]. Beijing: China Ocean Press, 2015. [本文引用:1]
[5] 王宏斌, 杨木壮, 白志琳, . 全球范围内天然气水合物构造背景及其稳定域研究综述[J]. 海洋地质, 2003(2): 12-21.
WANG Hongbin, YANG Muzhuang, BAI Zhilin, et al. A review of the global gas hydrate tectonic background and its stability domain[J]. Marine Geology, 2003(2): 12-21. [本文引用:2]
[6] RAJPUT S, THAKUR N K. Geological controls for gas hydrates and unconventionals[M]. Amsterdam: Elsevier, 2016. [本文引用:1]
[7] WANG Fangtian, ZHAO Bin, LI Gang. Prevention of potential hazards associated with marine gas hydrate exploitation: A review[J]. Energies, 2018, 11(9): 2384. [本文引用:1]
[8] SUESS E. Marine cold seeps and their manifestations: Geological control, biogeochemical criteria and environmental conditions[J]. International Journal of Earth Sciences, 2014(103): 1889-1916. [本文引用:3]
[9] 何家雄, 颜文, 祝有海, . 全球天然气水合物成矿气体成因类型及气源构成与主控因素[J]. 海洋地质与第四纪地质, 2013, 33(2): 121-128.
HE Jiaxiong, YAN Wen, ZHU Youhai, et al. Genetic types of gas hydrate in the world and their main controlling factors[J]. Marine Geology & Quaternary Geology, 2013, 33(2): 121-128. [本文引用:1]
[10] KVENVOLDEN K A. Comparison of marine gas hydrates in sediments of an active and passive continental margin[J]. Marine and Petroleum Geology, 1985, 2(1): 65-71. [本文引用:3]
[11] ESTEBAN M G. Methane hydrates in the chilean continental margin[J]. Electronic Journal of Biotechnology, 2003, 6(2): 80-84. [本文引用:2]
[12] 张金华, 方念乔, 魏伟, . 天然气水合物成藏条件与富集控制因素[J]. 中国石油勘探, 2018, 23(3): 35-46.
ZHANG Jinhua, FANG Nianqiao, WEI Wei, et al. Accumulation conditions and enrichment controlling factors of natural gas hydrate reservoirs[J]. China Petroleum Exploration, 2018, 23(3): 35-46. [本文引用:1]
[13] HYNDMAN R D, SPENCE G D. A seismic study of methane hydrate marine bottom simulating reflectors[J]. Journal of Geophysical Research, 1992, 97(B5): 6683-6698. [本文引用:1]
[14] WHITE R S, LOUDEN K E. The Makran continental margin: Structure of a thickly sedimented convergent plate boundary: Convergent margins: Field investigations of margin structure and stratigraphy[J]. AAPG Special Volumes: Memoir, 1982, A110: 499-518. [本文引用:1]
[15] 吴能友, 张海啟, 杨胜雄, . 南海神狐海域天然气水合物成藏系统初探[J]. 天然气工业, 2007, 27(9): 1-6.
WU Nengyou, ZHANG Haiqi, YANG Shengxiong, et al. Preliminary discussion on natural gas hydrate (NGH) reservoir system of Shenhu area, north slope of South China Sea[J]. Natural Gas Industry, 2007, 27(9): 1-6. [本文引用:1]
[16] 吴能友, 梁金强, 王宏斌, . 海洋天然气水合物成藏系统研究进展[J]. 现代地质, 2008, 22(3): 356-362.
WU Nengyou, LIANG Jinqiang, WANG Hongbin, et al. Marine gas hydrate system: State of the art[J]. Geoscience, 2008, 22(3): 356-362. [本文引用:1]
[17] 卢振权, 吴能友, 陈建文, . 试论天然气水合物成藏系统[J]. 现代地质, 2008, 22(3): 363-375.
LU Zhenquan, WU Nengyou, CHEN Jianwen, et al. Preliminary discussion on gas hydrate geological system[J]. Geoscience, 2008, 22(3): 363-375. [本文引用:1]
[18] RIEDEL M, COLLETT T S, MALONE M. Expedition 311 synthesis: Scientific findings[R]//Proceedings of the Integrated Ocean Drilling Program volume 311. Washington, DC: Integrated Ocean Drilling Program Management International Inc. , 2010. [本文引用:1]
[19] COLLETT T S. The gas hydrate petroleum system[C]//Proceedings of the 8th International Conference on Gas Hydrates, Beijing, China: China Geological Survey, 2014. [本文引用:1]
[20] 卜庆涛, 胡高伟, 业渝光, . 天然气水合物成藏体系研究进展[J]. 新能源进展, 2015, 3(6): 435-443.
BU Qingtao, HU Gaowei, YE Yuguang, et al. Research progress in natural gas hydrate accumulation system[J]. Advances in New and Renewable Energy, 2015, 3(6): 435-443. [本文引用:3]
[21] 魏伟, 张金华, 吝文, . 天然气水合物成藏机理及主控因素[J]. 新疆石油地质, 2010, 31(6): 563-566.
WEI Wei, ZHANG Jinhua, LIN Wen, et al. Mechanism and controlling factors of natural gas hydrate accumulation[J]. Xinjiang Petroleum Geology, 2010, 31(6): 563-566. [本文引用:1]
[22] FRYER P, WHEAT C G, WILLIAMS T, et al. Mariana serpentinite mud volcanism exhumes subducted seamount materials: Implications for the origin of life[J]. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2020, 378(2165): 20180425. [本文引用:2]
[23] 张金华, 苏明, 魏伟, . 含气流体运移与天然气水合物成藏[J]. 地质科技情报, 2017, 36(2): 176-185.
ZHANG Jinhua, SU Ming, WEI Wei, et al. Relationship between gas-bearing fluids migration and accumulation of natural gas hydrate[J]. Geological Science and Technology Information, 2017, 36(2): 176-185. [本文引用:1]
[24] 梁永兴, 曾溅辉, 郭依群, . 神狐钻探区天然气水合物成藏地质条件分析[J]. 现代地质, 2013, 27(2): 425-434.
LIANG Yongxing, ZENG Jianhui, GUO Yiqun, et al. Analysis of natural gas hydrate accumulation conditions of Shenhu prospect[J]. Geoscience, 2013, 27(2): 425-434. [本文引用:1]
[25] 何家雄, 苏丕波, 卢振权, . 南海北部琼东南盆地天然气水合物气源及运聚成藏模式预测[J]. 天然气工业, 2015, 35(8): 19-29.
HE Jiaxiong, SU Pibo, LU Zhenquan, et al. Prediction of gas sources of natural gas hydrate in the Qiongdongnan Basin, northern South China Sea, and its migration, accumulation and reservoir formation pattern[J]. Natural Gas Industry, 2015, 35(8): 19-29. [本文引用:1]
[26] TRÉHU A, RUPPEL C, HOLLAND M, et al. Gas hydrates in marine sediments: Lessons from Scientific Ocean Drilling[J]. Oceanography, 2006, 19(4): 124-142. [本文引用:1]
[27] 苏丕波, 梁金强, 付少英, . 南海北部天然气水合物成藏地质条件及成因模式探讨[J]. 中国地质, 2017, 44(3): 415-427.
SU Pibo, LIANG Jinqiang, FU Shaoying, et al. Geological background and accumulation models of gas hydrate reservoir in northern South China Sea[J]. Geology in China, 2017, 44(3): 415-427. [本文引用:1]
[28] 苏丕波, 何家雄, 梁金强, . 南海北部陆坡深水区天然气水合物成藏系统及其控制因素[J]. 海洋地质前沿, 2017, 33(7): 1-10.
SU Pibo, HE Jiaxiong, LIANG Jinqiang, et al. Natural gas hydrate migration and accumulation system and its controlling factors on northern deep water slope of the South China Sea[J]. Marine Geology Frontiers, 2017, 33(7): 1-10. [本文引用:1]
[29] LEE M W, COLLETT T S. Scale-dependent gas hydrate saturation estimates in sand reservoirs in the Ulleung Basin, East Sea of Korea[J]. Marine and Petroleum Geology, 2013(47): 195-203. [本文引用:1]
[30] NARITA H. Introduction of MH21(Research Consortium for Methane Hydrate Resources in Japan) and current topics in production method & modeling of methane hydrate[C]//Fifth ISOPE Ocean Mining Symposium, 15-19 September 2003, Tsukuba, Japan. [本文引用:1]
[31] LI Qing, WANG Jiasheng, WANG Xiaoqin, et al. Stable Carbon isotopic response of the benthic foraminifera from IODP311 to the marine methane hydrate geo-system[J]. Advances in Earth Science, 2008, 23(11): 1161-1166. [本文引用:1]
[32] LU Haihong, MOUDRAKOVSKI I, JUDITH S, et al. The characteristics of gas hydrates recovered from northern Cascadian margin by IODP expedition 311[C]//AGU 2006 Fall Meeting. San Francisco: American Geophysical Union, 2006. [本文引用:1]
[33] 高莉玲. 智利三联点附近气水合物带底部存在游离气体[J]. 海洋地质动态, 1994(2): 24.
GAO Liling. There is free gas at the Bottom of the gas hydrate zone near the triple point of Chile[J]. Marine Geology Development, 1994(2): 24. [本文引用:1]
[34] MAX D M, 韩乃明. 大洋甲烷水合物是一种“未勘探”的气资源[J]. 海洋地质译丛, 1996(4): 48-56.
MAX M D, HAN Naiming. Oceanic methane hydrate is an "unexplored" gas resource[J]. Translation of Marine Geology, 1996(4): 48-56. [本文引用:1]
[35] RIEDEL M, COLLETT T S, KUMAR P, et al. Seismic imaging of a fractured gas hydrate system in the Krishna-Godavari Basin offshore India[J]. Marine and petroleum geology, 2010, 27(7): 1476-1493. [本文引用:1]
[36] 杨木壮, 潘安定, 沙志彬. 陆缘地区天然气水合物成藏地质模式[J]. 海洋地质与第四纪地质, 2010, 30(6): 85-90.
YANG Muzhuang, PAN Anding, SHA Zhibin. Geological models of gas hydrates deposits along the continental margin[J]. Marine Geology & Quaternary Geology, 2010, 30(6): 85-90. [本文引用:1]
[37] 蔡峰, 闫桂京, 梁杰, . 大陆边缘特殊地质体与水合物形成的关系[J]. 海洋地质前沿, 2011, 27(6): 11-15.
CAI Feng, YAN Guijing, LIANG Jie, et al. The relationship between special geological bodies and hydrate formation at continental margin[J]. Marine Geology Letters, 2011, 27(6): 11-15. [本文引用:2]
[38] DILLON W P, PAULL C K. Marine gas hydrates the geophysical evidence[M]. Boston: Butterworth, 1983. [本文引用:1]
[39] KHLEBNIKOV V N, ANTONOV S V, MISHIN A S, . 多孔介质中天然气水合物生成的主要影响因素[J]. 天然气工业, 2017, 37(5): 38-45.
KHLEBNIKOV V N, ANTONOV S V, MISHIN A S, et al. Major factors influencing the formation of natural gas hydrates in porous media[J]. Natural Gas Industry, 2017, 37(5): 38-45. [本文引用:1]
[40] 张保勇, 周莉红, 刘昌岭, . 不同粒度沉积物介质对气体水合物成核的影响[J]. 天然气工业, 2018, 38(5): 148-155.
ZHANG Baoyong, ZHOU Lihong, LIU Changling, et al. Influence of sediment media with different particle sizes on the nucleation of gas hydrate[J]. Natural Gas Industry, 2018, 38(5): 148-155. [本文引用:1]
[41] 蔡建超, 夏宇轩, 徐赛, . 含水合物沉积物多相渗流特性研究进展[J]. 力学学报, 2020, 52(1): 208-223.
CAI Jianchao, XIA Yuxuan, XU Sai, et al. Advances in multiphase seepage characteristics of natural gas hydrate sediments[J]. Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(1): 208-223. [本文引用:1]
[42] LIU Jinlong, HAECKEL M, RUTQVIST J, et al. The mechanism of methane gas migration through the gas hydrate stability zone: Insights from numerical simulations[J]. Journal of Geophysical Research, 2019, 124(5): 4399-4427. [本文引用:2]
[43] 胡高伟, 业渝光, 张剑, . 沉积物中天然气水合物微观分布模式及其声学响应特征[J]. 天然气工业, 2010, 30(3): 120-124.
HU Gaowei, YE Yuguang, ZHANG Jian, et al. Micro-models of gas hydrate and their impact on the acoustic properties of the host sediments[J]. Natural Gas Industry, 2010, 30(3): 120-124. [本文引用:1]
[44] ROGERS G, DRAGERT H. Episodic tremor and slip on the Cascadia subduction zone: The chatter of silent slip[J]. Science, 2003, 300(5627): 1942-1943. [本文引用:1]
[45] TRÉHU A M, BOHRMANN G, RACK F R, et al. Proceedings of the Ocean Drilling Program, initial reports volume 204[R]. Texas: ODP, 2003. [本文引用:2]
[46] RIEDEL M, COLLETT T S, MALONE M J, et al. IODP Expedition 311 summary. Proceedings of the Integrated Ocean Drilling Program[R]. Washington, DC: Integrated Ocean Drilling Program Management International Inc, 2006. [本文引用:4]
[47] 王家生, 高钰涯, 李清, . 沉积物粒度对水合物形成的制约: 来自IODP 311航次证据[J]. 地球科学进展, 2007, 22(7): 659-665.
WANG Jiasheng, GAO Yuya, LI Qing, et al. Grain size constraint on gas hydrate occurrence: Evidence from sediment size during IODP 311[J]. Advances in Earth Science, 2007, 22(7): 659-665. [本文引用:2]
[48] BABA K, YAMADA Y. BSRs and associated reflections as an indicator of gas hydrate and free gas accumulation: An example of accretionary prism and forearc basin system along the Nankai trough, off Central Japan[J]. Resource Geology, 2004, 54(1): 11-24. [本文引用:1]
[49] CHOI P Y, NAKAE S, KIM H. Fault tectonic analysis of Kii peninsula, Southwest Japan: Preliminary approach to neogene paleostress sequence near the Nankai subduction zone[J]. Island Arc, 2011, 20(4): 455-476. [本文引用:1]
[50] HAMMERSCHMIDT S B, ANDRE HUPERS E E D, KOPF A. Limitation of fluid flow at the Nankai trough megasplay fault zone[J]. Geo-Marine Letters, 2013(33): 405-418. [本文引用:1]
[51] PAULL C K, MATSUMOTO R, WALLACE P J, et al. Ocean Drilling Program, Leg 164, gas hydrate sampling on the Blake Ridge and Carolina Rise[R]. College Station, Texas: Proceedings of the Ocean Drilling Program, 1996. [本文引用:1]
[52] BETHKE C M. A numerical model of compaction-driven groundwater flow and heat transfer and its application to the paleohydrology of intracratonic sedimentary basins[J]. Journal of Geophysical Research, 1985, 90(B8): 6817-6828. [本文引用:1]
[53] RUTQVIST J, WU Y S, TSANG C F, et al. A modeling approach for analysis of coupled multiphase fluid flow, heat transfer, and deformation in fractured porous rock[J]. International Journal of Rock Mechanics and Mining Sciences, 2002, 39(4): 429-442. [本文引用:1]
[54] YIN Zhenyuan, KHURANA M, TAN H K, et al. A review of gas hydrate growth kinetic models[J]. Chemical Engineering Journal, 2018(342): 9-29. [本文引用:1]
[55] YOU K, FLEMINGS P B, MALINVERNO A, et al. Mechanisms of methane hydrate formation in geological systems[J]. Reviews of Geophysics, 2019, 57(4): 1146-1196. [本文引用:1]
[56] DAI S, SANTAMARINA J C, WAITE W F, et al. Hydrate morphology: Physical properties of sand s with patchy hydrate saturation[J]. Journal of Geophysical Research, 2012, 117(B11): B11205. [本文引用:1]
[57] 吕万军. 天然气水合物形成条件与成藏过程——理论、实验与模拟[D]. 广州: 中国科学院广州地球化学研究所, 2004.
LYU Wanjun. Gas Hydrate formation conditions and process of reservoir formation: Theory, experiment and simulation[D]. Guangzhou: Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, 2004. [本文引用:2]
[58] LI Qing, WANG Jiasheng, CHEN Jianwen, et al. Stable carbon isotopes of benthic foraminifers from IODP Expedition 311 as possible indicators of episodic methane seep events in a gas hydrate geosystem[J]. Palaios, 2010, 25(10): 671-681. [本文引用:2]
[59] WASEDA A, UCHIDA T. The geochemical context of gas hydrate in the eastern Nankai trough[J]. Resource Geology, 2004, 54(1): 69-78. [本文引用:2]
[60] TOMARU H, LU Zunli, SNYDER G T, et al. Origin and age of pore waters in an actively venting gas hydrate field near Sado Island , Japan Sea: Interpretation of halogen and 129I distributions[J]. Chemical Geology, 2007, 236(3/4): 350-366. [本文引用:]
[61] 刘杰, 孙美静, 苏明, . 海底泥底辟(泥火山)对天然气水合物成藏的影响[J]. 地质科技情报, 2015, 34(5): 98-104.
LIU Jie, SUN Meijing, SU Ming, et al. Influence of submarine mud diapirs (mud volcanoes) on gas hydrate accumulation[J]. Geological Science and Technology Information, 2015, 34(5): 98-104. [本文引用:1]