海洋沉积物微量元素地球化学特征对天然气水合物勘探的指示意义

2022-02-18于哲邓义楠陈晨曹珺方允鑫蒋雪筱黄毅

海洋地质与第四纪地质 2022年1期

于哲，邓义楠，陈晨，曹珺，方允鑫，蒋雪筱，黄毅

1. 南方海洋科学与工程广东省实验室（广州），广州 511458

2. 自然资源部海底矿产资源重点实验室，中国地质调查局广州海洋地质调查局，广州 510075

天然气水合物是在高压低温环境下由甲烷、乙烷和二氧化碳等低分子量气体与水分子结合形成的固态物质，是一种新型绿色低碳的高效能源，已被国际社会广泛认为是未来解决能源枯竭和生态恶化等问题的重要战略资源之一。另外由水合物分解所诱发的海底滑坡和温室效应等，使得海洋天然气水合物在能源潜力、海洋工程安全和全球气候变化等方面都受到更多关注[1-4]。因此，天然气水合物及其赋存的沉积物已经成为海洋地质学研究的重点领域[5-9]。

天然气水合物失稳释放常导致大量的甲烷渗漏至孔隙水甚至海水当中，对海洋微量元素分布特征造成重大影响[10]。因此，海洋沉积物微量元素地球化学特征亦常用于天然气水合物勘探的研究[11-12]。例如，早期将沉积物中的“钡锋”作为过去和现在甲烷渗漏的指示[13-14]。近年来，Sato 等[15]、Chen 等[11]和Deng 等[10]发现天然气水合物富集区沉积物中钼（Mo）元素的极大富集，可作为下部可能存在天然气水合物的有效判别指标。Feng 和Chen[16]和Ge 等[17]发现在天然气水合物分解强度较弱的情况下，形成于厌氧条件下的沉积物相对富集部分氧化还原敏感元素，如：Mo、V、Ni、U、Co 和Cd。海马冷泉和墨西哥湾的沉积物还发现As 和Sb 元素的富集，Fe的氢氧化物扮演了重要的赋存介质[18-19]。

然而，前期关于天然气水合物富集区沉积物微量元素地球化学的研究主要集中于浅表层的沉积物[12]，浅层沉积物地球化学特征至今鲜有报道[11]。天然气水合物在沉积物中赋存是多层段的，各部分的微量元素特征对于深入了解甲烷渗漏对海洋微量元素循环的影响及其对指导天然气水合物的勘探具有重要意义。本文对中国南海神狐海域开展天然气水合物钻探，对取得的浅层沉积物样品进行主、微量元素和总有机碳（TOC）地球化学分析，旨在探讨：①受甲烷渗漏影响的钻孔沉积物微量元素地球化学变化特征；②氧化还原状态的变化情况；③沉积物微量元素特征对天然气水合物勘探的指导意义。

1 地质背景

我国在南海北部神狐海域开展过多个航次的天然气水合物调查研究。前期的天然气水合物调查结果表明，南海北部神狐及其邻近海域是重要的水合物富集区，存有丰富的天然气水合物资源，早期研究已经发现地层深部的多处似海底反射层（BSR）[20]，深度在海底以下160～220 m。研究区海底广泛存在断层和褶皱，可充分用于海底甲烷的迁移[20]。南海多数冷泉分布于该海域[21]。

2015 年，广州海洋地质调查局 GMGS-3 航次在南海北部神狐海域进行了23 个站位随钻测井（图1），测井资料表明在每个站位都发现了水合物，局部水合物饱和度高达64%，初步分析认为储层可能为富含有孔虫的黏土质粉砂。在其中的W19、W18、W17 和W11 站位进行了原位测量和取心并获得不可视的高饱和度水合物实物样品（图1b），钻探水深1 249～1 310 m，实际钻探深度为222～315 mbsf（海底以下深度，m），高饱和度水合物产出站位位于钻探区东南角海底海脊下段位置（图1c）。

图1 南海北部神狐海域GMGS-3 航次钻探取样孔位置图及测井曲线图a. 神狐海域研究区位置，b. GMGS-3 钻探区井位分布，c. GMGS-3 钻探区水合物随钻测井曲线[22]。Fig.1 Location and well log of Shenhu area in the South China Seaa. research area in Shenhu area, b. location site in GMGS-3, c. logging curves of GMGS-3[22].

2 样品与分析

对采集的柱状样沉积物中的常量、微量元素进行了分析，分析方法如下：

海洋沉积物样品经（110±2）℃烘干后细碎至0.074 mm，然后进行各测试项目化学分析。有机碳分析试样在低于80 ℃温度下烘干，磨碎至0.149 mm进行分析。

常量元素分析：SiO2及Al2O3、Fe2O3、CaO、MgO、K2O、Na2O、MnO、P2O5、TiO2等氧化物采用无水四硼酸锂-无水偏硼酸锂-氟化锂为熔剂，将样品熔融成为玻璃片，采用荷兰PANalytical 公司AXIOSX 型X 射线荧光光谱仪进行测试。烧失量（灼失量）采用重量法分析，即试样在1 000 ℃灼烧至恒重，以灼烧减少的量为烧失量。

微量元素分析：试样以氢氟酸、硝酸、高氯酸分解，以稀盐酸溶解盐类，制成酸度10%的试液，采用美国PE 公司Optima 4300DV（ICP-OES）全谱直读等离子体原子发射光谱仪进行测试。部分微量元素（U、Mo、V、Cu、Ni、Cr、Sc）试样经硝酸、氢氟酸、高氯酸分解，蒸发除去氢氟酸，以盐酸加热溶解盐类，转化为硝酸介质，以铑作内标溶液，采用美国Thermo 公司X2 型等离子体质谱仪（ICP-MS）进行测试。

有机碳分析：在浓硫酸介质中，加入一定量的标准重铬酸钾，在加热条件下将样品中有机碳氧化成二氧化碳，剩余的重铬酸钾用硫酸亚铁标准溶液回滴，按照重铬酸钾的消耗量，计算样品中有机碳的含量。沉积物半定量矿物分析参照Deng 等[10]，矿物组成使用X 射线衍射（XRD）方法在Rigaku DXR 3000 上测定，样品通过分段扫描的方式，借助40 mA 电流和40 kV 的加速电压，在5～60° （2θ）区间以0.02 °/s 的速度进行扫描，半定量分析通过程序SIROQUAN T 完成。以上实验均在广州海洋地质调查局完成，分析误差符合标准，分析结果真实可靠。

3 分析结果

本次分析的4 个钻孔站位的结果如下：

W19 样品主量元素以SiO2、Al2O3和CaO 的含量最高，SiO2含量为30.04%～53.65%，Al2O3的含量为5.85%～15.26%，CaO 含量为8.04%～28.64%。在站位顶部出现SiO2和Al2O3的高值，随深度向下含量逐渐下降，而CaO则出现相反的变化趋势。自下而上主量元素含量均呈规律性变化（图2）。

图2 W19 站位主量元素含量垂向变化图Fig.2 Variation of major elements geochemistry at Site W19

上段：52～111 m 层段陆源物质含量高，生源物质含量低。SiO2等主要与粗粒陆源碎屑有关的元素含量高且较稳定，大部分含量均大于45%。Al2O3主要代表了细粒陆源物质沉积，在该段其变化特征与硅质组分类似。生源物质CaO 的含量相对低且变化不大，变化范围为8.04～15.15%；下段：111 m至剖面底部（167.7 m），SiO2和Al2O3含量经历了急剧下降再上升的过程，自生源物质的CaO 含量表现为上升再下降，表明该站位沉积物表层主要为陆源碎屑来源，中段可能具有碳酸盐岩夹层或者富含生物壳体，向下陆源物质开始增多，自生沉积作用减弱。Fe2O3和MnO 在底层出现极高值，K2O 和TiO2变化特征与SiO2和Al2O3类似，同样受到了陆源碎屑物质的干扰。沉积物中矿物含量钙质占绝大多数，在125 m 处含有少量重晶石。未含水合物层段，含有少量硅质组分，黄铁矿含量极少。在钻孔底部水合物赋存区，黄铁矿含量显著升高，最高可达55%。

U、V、Ni、Zn、Co、Cu 和Sc 等微量元素的变化特征与SiO2和Al2O3有着相似的变化趋势（图3），均表现为上段的高值和下段逐渐下降再上升的过程，说明这些元素的来源可能部分为陆源碎屑物质，且富集在细粒和粗粒沉积物中。Mo、Ni、Zn、Co 和Cu 含量分别为0.24～46.2、19.3～135、54.0～231、7.60～36.3 和（10.4～60.4）×10-6，这些元素含量在上部相对较低，但在下部140～141.4 和166.5～167.7 mbsf 层段具极高值，与Fe2O3和MnO 的变化特征一致，意味着部分层段沉积物自生物质的富集作用。

图3 W19 站位微量元素含量垂向变化图Fig.3 Variation of trace elements geochemistry at Site W19

W18 站位特征与W19 站位相似，可分为两段：第一段：50～140 m，该段SiO2和Al2O3含量在表层较高，之后开始逐渐下降至140 m 处达最低值，CaO 浓度则出现明显的上升（图4），表明该段表层主要以陆源碎屑沉积为主，之后碎屑物质减少，与陆源物质相关的TiO2浓度的下降也验证了这一点，与此相反，自生物质的CaO 含量自上而下逐渐升高；第二段：140 m 至站位底部（173.7 m），SiO2和Al2O3含量开始急剧上升，相反的CaO 浓度下降（图4），表明该段与上段特征相反，陆源物质开始增多。矿物特征方面，本站位与W19 较为相似，主要以钙质矿物为主，黄铁矿含量次之，且多集中于下部的水合物赋存区，个别样品（159.5 m）出现重晶石，部分样品有少量石英的分布，含量为1%～2%。

图4 W18 站位主量元素含量垂向变化图Fig.4 Variation of major elements geochemistry at Site W18

该站位的微量元素特征与W19 站位类似，微量元素U、V、Ni、Co、Cu 和Sc 等组分与SiO2和Al2O3变化趋势相同（图5），低值均出现在剖面中段，可能意味着上述元素的来源与碎屑物质关系密切。但是，U、Mo、Ni、Co 和Cu 的浓度范围分别为（1.54～5.67）×10-6、（0.25～33.3）×10-6、（19.4～102）×10-6、（6.3～23.6）×10-6和（16.8～63.5）×10-6。这些元素在站位浅层含量较低，在下半部分逐渐升高，最高值基本位于146 ～159.7 mbsf 层段，均与Fe2O3和MnO 的高值对应，表明了沉积物中微量元素的富集还与沉积环境有关（图5）。

图5 W18 站位微量元素含量垂向变化图Fig.5 Variation of trace elements geochemistry at Site W18

W17 站位样品SiO2和Al2O3含量从剖面自上而下呈逐渐上升的趋势，与TiO2和K2O 相似（图6），表明站位沉积物碎屑物质来源逐渐升高的过程。而CaO、MgO 和P2O5含量逐渐下降，海洋自生物质组分所占比重减少。本站位矿物特征相对复杂，除了以钙质成分为主以外，在52.5～65.95 m 层段有较多的石膏富集，平均含量为9%。石英在125.1～131.5 和198～266.05 m 层段出现显著富集，最高可达85%。在剖面底部的水合物富集区，同样出现了黄铁矿的显著富集。

图6 W17 站位主量元素含量垂向变化图Fig.6 Variation of major elements geochemistry at Site W17

V、Ni、Zn、Co、Cu 和Sc 等微量元素随剖面深度逐渐上升（图7），表示上述元素来源可能多为陆源碎屑物质，而U 与CaO 变化类似，其来源主要为海洋自生物质。Fe2O3在底层出现极高值，Mo、Ni、Co 和Cu 含量变化区间分别为（0.23～119）×10-6、（24.4～44.2）×10-6、（8.80～18.2）×10-6和（12.7～49.0）×10-6，这些元素在剖面浅表层含量较低，随后逐渐升高，在206～241.9 mbsf 层段均出现极度富集。

图7 W17 站位微量元素含量垂向变化图Fig.7 Variation of trace elements geochemistry at Site W17

W11 站位样品SiO2的平均含量为研究区之最，平均值为57.50%，且Al2O3浓度也较高，表明该站位的陆源碎屑物质含量最高。纵向上碎屑物质含量随剖面的下降呈逐渐上升的趋势；与海洋自生物质相关的CaO 和P2O5含量在表层出现高值，之后开始急剧降低，33.5 m 之后趋于稳定。Fe2O3和MnO 同样在站位底部出现极大值（图8），与之对应的Mo、Ni、Zn、Co 和Cu 含量变化范围分别为（0.24～141）×10-6、（21.44～93.9）×10-6、（53.3～285）×10-6、（8.52～28.7）×10-6和（10.1～56.4） ×10-6，这些微量元素均在浅表层含量较低，而在146.5～188.1 mbsf 层段出现极大值（图9）。该站位矿物特征较为不同，50～53 m 深度，以钙质为主。53.5～123.53 m 层段，出现特征的石膏矿物为主，最高可达94%，黄铁矿含量次之。随后至剖面底部，黄铁矿含量骤增。

图8 W11 站位主量元素含量垂向变化图Fig.8 Variation of major elements geochemistry at Site W11

图9 W11 站位微量元素含量垂向变化图Fig.9 Variation of trace elements geochemistry at Site W11

4 讨论

4.1 钡元素变化特征

海底沉积物中钡循环与向上渗漏的甲烷之间存在很好的地球化学耦合关系[13,23]。海洋“Ba 循环”受沉积物孔隙水中SO42-的浓度梯度制约[24]，如果BaSO4的溶解扩散强于其沉积埋藏，就会在SO42-亏损深度附近的沉积物中形成一个范围较窄的、不稳定的Ba 浓度异常，称之为“Ba 峰”[25-26]。在天然气水合物分布区，沉积物中的CH4厌氧氧化普遍存在，它在海底缺氧带内以硫酸盐为氧化剂，通过微生物氧化以甲烷为主的碳氢化合物，并在海底浅表层形成硫酸盐-碳氢化合物转化带（SHT）[25]，SHT 之上亏损碳氢化合物，之下亏损硫酸盐。富含Ba2+的冷泉渗漏流体向上扩散与孔隙水硫酸盐反应，在SHT 之上沉淀重晶石[27]。随着沉积物的埋藏，先前沉积的重晶石向下埋藏至硫酸盐亏损带，随之发生溶解转化为溶解钡并和冷泉流体一起向上扩散[28]。CH4释放通量的变化制约着SO42-亏损带的迁移，进而影响着“Ba 峰”深度的变化。因而，沉积物中的“古Ba 峰”的发育深度和强度往往与历史时期天然气水合物分解CH4释放的通量有关，进而成为海底天然气水合物勘查的新指标[29]。孕育一个“Ba 锋”至少需要1 万年到100 万年的时间。因此，“Ba 锋”代表了在过去至少1 万年的时间里该区海底一直进行着相对稳定的天然气渗漏活动[13]。如果“Ba 锋”在现今亏损硫酸盐的深度（SHT）以上位置出现，可能是“古Ba 峰”的存在，它所在的位置代表过去硫酸盐浓度为零的深度，表明了天然气渗漏活动在时间上的动态变化。

在甲烷渗漏区域，在沉积物中发生甲烷厌氧氧化反应（the anaerobic oxidation of methane, AOM），消耗孔隙水中大量的硫酸盐[16,19]。在孔隙水亏损硫酸盐时，深部孔隙水中Ba 离子会大量富集，富Ba 孔隙水沿裂隙运移至海水-沉积物界面附近，与富硫酸盐孔隙水或海水相互作用而形成重晶石：Ba2++ SO42-=BaSO4[13,26]。在W19、W18、W17 和W11 站位天然气水合物赋存区域，部分层位出现重晶石矿物，是Ba 元素富集的主要载体，意味着研究区天然气水合物可能存在分解的现象，渗漏的甲烷引起了沉积物矿物成分的改变。

调查区各站位的Ba 元素变化特征如图10 所示，4 个钻孔的底部沉积物均出现不同程度的“Ba 峰”，分别位于W19、W18、W17 和W11 柱状样的144.15、159.5、237 和183 m 处，该区域层段与发现的水合物层段较好对应，进一步表明了研究区域沉积物中有大规模甲烷的泄露。孔隙水中SO42-的亏损受CH4释放通量变化的制约，在较长的时间内，如果CH4释放过程是连续的、且释放通量没有发生显著变化，那么就会在SMI 内发育明显的、且宽度较大的“Ba 峰”[25]。如果CH4的释放通量发生过脉冲式的剧烈变化，那么就会在多个深度上发育多个离散式的“Ba 峰”。调查区呈现的“Ba 峰”宽度较窄，部分钻孔（如W19）呈现出较离散型的Ba 富集，表明神狐海域天然气水合物的泄露是阶段性且不连续的。

图10 调查区各站位Ba 元素含量变化特征Fig.10 Variation of Ba concentration at W11, W17, W18, W19 in the research area

4.2 钼元素变化特征

微量元素钼具有多种价态，包括Mo2-到Mo6+，Mo6+是钼的最常见的氧化态形式，在溶液中以MoO42-的形式迁移，在还原条件下，多以Mo4+和Mo5+的形式存在[30-33]，Mo 是一种对环境的氧化-还原状态变化非常敏感的一种元素[10,34]。在海洋的表层，钼主要以钼酸盐（MoO42-）的形式存在于海水中，它不会因为氧化环境中存在大量浮游生物而富集，与碳酸钙等亲和力非常弱。但是，钼对含Mn 的氧化物和氢氧化物非常敏感，在氧化环境中大量的Mo 因被Mn 的氧化物和氢氧化物吸附在表面而富集[32-33]。硫化氢的含量对沉积物富集Mo 起关键作用。比如海水在次氧化环境中富集的Mo 含量低于硫化环境中的Mo 含量[10]。

研究海域中，各站位含天然气水合物的层段Mo 的含量在站位底部出现明显的峰值（图11），与“Ba 峰”重合，其原因可能与水合物分解产生的硫化氢有关[10]。天然气水合物分解后，释放出大量的甲烷气体，通过扩散或沿断裂向海底表层迁移[35]。在海底浅表层沉积物中，甲烷被氧化，硫酸根被还原，反应将按照下式进行[36]：

图11 调查区各站位Mo 元素含量变化特征Fig.11 Variation of Mo concentration at W11, W17, W18, W19 in the research area

Mo 在氧化的环境中，主要吸附于Mn 的氧化物中[12]。在沉积环境变为硫化环境时，出现自由状态的H2S，Mo 的形态由钼酸盐（MoO42-）变为硫代钼酸根（MoS4-xOx2-），并大量富集于Fe 的硫化物相中[18]。W19、W18、W17 和W11 柱状样水合物赋存区域内分别在141.2～167.7、168.5～173.7、190.5～266.6、150～210.36 mbsf 层段出现Mo 元素富集（图11）的同时均显现出黄铁矿的显著增加，意味着研究区底部可能为硫化环境，黄铁矿为Mo 元素的主要赋存载体。

还原性水体，可以还原吸收水体中的Mo 等变价元素[11]。在同一地区的沉积物中，Mo 含量与沉积物中有机碳的含量往往呈正相关关系，但是如果环境中有HS-存在，Mo 就会大量富集而与TOC 含量无关[33]。在有机碳对Mo 含量的图解中，不同环境下形成的沉积物有不同的斜率[18,37]。如图12 所示，研究海域样品的TOC 大多与Mo 的含量存在一定的相关关系，总体上样品的TOC 与Mo 的含量呈一定的相关关系，指示着样品的沉积环境大部分可能处于厌氧环境[33]。部分样品出现Mo 的高度富集且与TOC 相关性较差（图12），表明此时的沉积环境可能出现硫化环境，这些硫化环境的样品均分布于各站位底部发现天然气水合物的层段（图1）。

图12 调查区沉积物总有机碳与Mo 相关关系图Fig.12 The correlation of Mo and TOC concentrations in the sediments

富含天然气水合物的区域，甲烷与硫酸根离子反应产生了HS-，形成硫化氢气体[38]。故该区域虽然可能并非封闭或静水环境，TOC 含量也相对不高，但是水合物的分解或许能产生一定量的HS-或H2S，形成局部的硫化环境从而促进Mo 等金属离子的富集[10]。海底天然气水合物一般较多存在于有机碳含量较高（＞1%）、相对还原的沉积环境中[11]。研究海域的的总有机碳（TOC）含量平均值为0.78%，有机碳含量相对较低。但是，部分站位底部（如W19）富含天然气水合物的的站位TOC 的值可高达6.13%。研究区沉积物TOC 和还原环境的特征均较好的符合海底天然气水合物的形成条件。