南海南部浅表层柱状沉积物孔隙水地球化学特征对甲烷渗漏活动的指示

马晓理，刘丽华，徐行，金光荣，魏雪芹，翟梦月

1. 中国科学院广州能源研究所，中国科学院天然气水合物重点实验室，广州 510640

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

3. 中国地质调查局广州海洋地质调查局，广州 510760

4. 中国地质大学（北京）海洋学院，北京 100083

海底甲烷渗漏活动是指富碳氢化合物（通常以甲烷为主）的流体沿裂隙、断裂等高渗通道及多孔介质以渗流或喷溢的形式从深部向上运移至海底浅部的过程。渗漏的甲烷主要来源于有机碳的原地转化、深部埋藏的天然气水合物分解[1-3]或油气泄漏[4]。甲烷的温室效应强度是二氧化碳的20倍以上，海底中甲烷的大量释放会对全球气候变化产生重大影响。然而研究表明近90%的海底渗漏甲烷被厌氧氧化（Anaerobic Oxidation of Methane, AOM）所消耗而未能释放到大气中[5]。AOM是一种由微生物介导的氧化还原过程，参与的电子受体包括锰氧化物、硝酸盐、铁氧化物、硫酸盐等[6-8]。其中硫酸盐在海底沉积物中广泛存在，通过上覆海水向下扩散或原位储存于海洋沉积物孔隙水中，在沉积物早期成岩过程中发挥重要作用[8]。因此，在甲烷渗漏环境中，硫酸盐作为主要的氧化剂驱动甲烷厌氧氧化作用的发生（Sulfate-Driven Anaerobic Oxidation of Methane, SD-AOM）[6]，简化的化学反应式见式（1）。

SD-AOM发生在一个狭窄的硫酸盐-甲烷过渡带（Sulfate Methane Transition Zone，SMTZ）内[9-14]。反应产物HS−进一步与孔隙水的Fe2+、氧化还原敏感元素（如Mo）等形成自生硫化物矿物；HCO3−提高孔隙水中碱度的同时，与孔隙水中部分碱土金属元素（如Ca2+、Mg2+、Sr2+）结合形成自生碳酸盐岩。这些自生硫化物矿物（如黄铁矿）、自生碳酸盐岩种类以及Mo元素在沉积物中的富集程度可用来推测过去甲烷渗漏的强度和期次[15-19]。相同地，SD-AOM及其相伴的自生矿物形成过程均被记录在沉积物孔隙水中，因此，孔隙水地球化学特征对追踪甲烷渗漏及相关地球化学反应具有一定指示意义[13-14,20-22]。然而，有机质硫酸盐还原作用（Organoclastic Sulfate Reduction，OSR）[23]与SD-AOM同时发生时，可能会使得追踪甲烷渗漏活动变得困难。OSR是沉积物中有机质与硫酸盐发生的氧化还原反应，化学反应式见式（2）：

由于OSR与SD-AOM均是消耗硫酸盐、产生硫化氢还原气体和产生溶解无机碳（DIC）的过程。同时发生时，很难区分二者的相对贡献，进而影响对甲烷渗漏活动的评估。基于此，学者提出了很多识别方法，其中最常见的方法是用DIC产生量与硫酸盐消耗量的比值来辨别孔隙水中OSR与SD-AOM的相对比例[13-14,24]。但若SMTZ深度以下发生了产甲烷作用，则其产生的DIC向上迁移并混合，可能会提高SMTZ及其上部的DIC浓度（如大于30 mM），使得此方法无效[25-29]。因此，又有学者提出稳态反应-传输模型来定量研究OSR和SD-AOM作用过程以及二者消耗硫酸盐的相对占比[21,30]。因为甲烷作用以及甲烷厌氧氧化这两个生物地球化学过程的碳酸氢根产生量与硫酸盐消耗量的比值都是2∶1[29,31]，所以此方法仍可保守地估算OSR和SDAOM对沉积物孔隙水中硫酸盐消耗的相对贡献。此外，OSR和SD-AOM产生的DIC的13C同位素特征差异也可以定性研判二者对硫酸盐还原反应的影响程度[9,24]。

目前，已在中国南海北部的琼东南、西沙、台西南、神狐海域等发现多个甲烷渗漏点[13,19-20,32]；然而对南海南部海域的甲烷渗漏活动仍知之甚少。地质和地球物理研究表明南海南沙海域的北康盆地，南沙海槽油气资源潜力巨大[33-35]。地震记录也表明在南沙海域发育有适宜烃类气体泄漏和迁移的通道[33,35-36]，并识别出了广泛的似海底反射层（BSR）[33]。粗略统计南沙海域北康盆地的石油地质资源量为8.86×108t，天 然 气 地 质 资 源 量 为14 855×108m3[35]。另一方面，在北康盆地发现自生黄铁矿与石膏共生的现象[37]以及较高的磁化率带[38]，可能与深部富烃类流体渗流有关。孔隙水地球化学分析和数值模拟综合研究推测北康盆地海域存在甲烷渗漏的迹象[21]。然而目前对北康盆地是否存在活动的甲烷渗漏尚无明确结论，而孔隙水地球化学特征是认识甲烷渗漏的重要证据。基于此，本文拟对采自南海南部北康盆地的3个重力沉积柱状沉积物孔隙水样品（BH-H75、BH-H13Y和BH-H61）进行甲烷浓度、DIC及其 δ13CDIC、阴离子（SO42−、Cl−、Br−）以及主微量元素（Ca2+、Mg2+、Sr2+、Ba2+）分析，探究SDAOM生物地球化学过程，进一步了解北康盆地研究区甲烷渗漏活动情况。

1 区域地质背景

北康盆地位于南海南部边缘（图1），是南沙中部海域重要的含油气盆地。盆地面积约6×104km2，水深为100～2 000 m。北康盆地位于南沙中部海域大陆坡上，地质构造十分复杂，大量发育北东向、北西向断裂，其中北东向断裂最具规模；盆地从东部到北部相邻的盆地依次为南沙海槽盆地、安渡北盆地、南薇东盆地和南薇西盆地，西南部通过走滑性质的廷贾断裂带与曾母盆地相连；盆地内主要发育中始新统—第四系，地层厚度变化较大，最厚处超过13 km[39]。盆地基底主要为前新生代变质岩和酸性-基性火成岩，其中在东部地区主要分布火成岩[40]。此外，北康盆地中也发育大量烃源岩；分别为中始新统、上始新统—下渐新统湖相泥岩、海陆过渡相泥岩、上渐新统—中中新统的浅海-半深海相泥岩[41-42]。区域烃源岩的大量分布、疏导体系、生物礁以及古地热场结构特征指示北康盆地油气潜力以及天然气水合物资源量巨大[35,41-43]，并在相邻的南沙海槽发现BSRs[44]。构造和气源条件表明北康盆地具备天然气水合物形成、发育的有力条件，为南海水合物成藏一级远景区之一[34,45]。

图 1 南海南部南沙海域北康盆地位置图[46]Fig.1 Location of the Beikang basin in Nansha sea area of the southern South China Sea [46]

2 材料与方法

2.1 孔隙水采集与预处理

2019年8月，广州海洋地质调查局在南海南部（工区为南沙海域）开展了HYSH201905航次的调查，进而评价天然气水合物资源储量。本文研究的柱状沉积物样品利用大型重力活塞取样器获得，样品信息见表1。现场取样后，沉积柱（直径约120 mm）立刻被切成1 m长度，柱子的顶底部使用塑料盖与磁带封闭，并水平放置。以20 cm为间距使用电钻机进行打孔，随后利用0.2 μm过滤膜的Rhizon进行孔隙水采集。在采集BH-H61沉积柱底部孔隙水过程中，闻到强烈的臭鸡蛋气味。采集的所有孔隙水样品均分为两份平行样，第一份平行样取1.8 mL孔隙水样，加入体积为20 μL 的2 M 硝酸酸化处理（pH＜2），用于测试主、微量元素；第二份平行样为剩余的孔隙水样，不作任何处理并全部盛入15 mL样品瓶中，立刻用橡胶塞密封。所有的孔隙水样品均储存于4 ℃的冰箱中，以供进一步研究需要。

表 1 南海南部北康盆地海域采集的3个沉积柱信息Table 1 Information of three sedimentary columns collected from the Beikang Basin in the southern SCS

2.2 孔隙水地球化学分析与扩散通量计算

孔隙水中的主要阳离子（Na+、K+、Ca2+、Mg2+）采用电感耦合等离子体光谱仪（ICAP-7200系列，Thermo，USA）测试。在测试前，样品稀释10倍，随后取8 mL转移至进样管中进行阳离子测试，测试精度大于5%。微量元素（Ba2+、Sr2+）采用电感耦合等离子体质谱仪（iCAP Q系列，Thermo，USA）进行测试，检出最低值约 0.06 μg/kg, 精度及准确度均优于5%。Cl−和SO42−浓度采用离子色谱仪（ICS-1100型，Thermo，USA）进行测试。平行试验结果显示Cl−和SO42−浓度的RSD（n=5）分别为0.35%和0.28%，均小于1%，说明实验具有较好的重复性。DIC及其δ13CDIC采 用 连 续 流 质 谱 仪（Delta V Advantage，Thermo，USA）进行测试。测试方法如下：首先在12 mL的样品瓶（Labco）中加入适量无水H3PO4，随后采用氦气进行排空处理，去除瓶内空气对样品碳同位素比值测定的影响；其次在经过氦气排空处理后的样品瓶中加入0.2 mL水样，将产生的CO2驱赶并转移至质谱仪中，测定同位素分馏平衡后的CO2同位素比值。δ13C值以VPDB国际标准作为参考标准。DIC浓度及其δ13CDIC的测试精度分别优于2%和0.2‰。甲烷浓度采用Agilent 6850系列气相色谱仪（Thermo，USA）进行测试，测试精度优于2%。所有离子的分析测试工作均在自然资源部第三海洋研究所分析测试中心完成。

沉积物孔隙水中硫酸盐的分子扩散通量采用菲克第一定律计算，公式见式（3）和式（4）。

其中，J代表硫酸盐扩散通量（mmol•m−2•a−1）；φ代表孔隙率，假设为75%[47]；Ds代表沉积物扩散系数（m2•s−1）；C为硫酸盐浓度（mmol•L−1）；x为沉积物深度（m）；Do为海水中硫酸盐的扩散系数，文中Do取5.72×10−10m2•s−1[48]。

3 结果

3.1 阴离子组分、CH4、DIC、δ13CDIC

3.2 部分碱土金属元素浓度

部分碱土金属元素浓度（Ca2+、Mg2+、Sr2+、Ba2+）深度剖面见图3（表2）。研究区3个沉积柱的Ca2+、Mg2+、Sr2+浓度均随深度逐渐降低，在BH-H75和BH-H13Y中仅观察到轻微的降低趋势，对应的Sr/Ca、Mg/Ca比值在深度剖面上轻微升高。而在BH-H61沉积柱中的Ca2+、Mg2+、Sr2+浓度的降低趋势更为显著，其中Ca2+、Mg2+浓度分别降低至沉积柱底部的3.24和39.84 mM，较典型海水值（分别为10.3和53.2 mM）分别贫化了68%和25%。Sr/Ca、Mg/Ca摩尔比值在BH-H61沉积柱底部达到峰值。Ba2+浓度在BH-H75和BH-H13Y沉积柱中均维持较低的水平（0.436～0.813 μM）。在BH-H61沉积柱上部（20～280 cmbsf），Ba2+浓度较低（0.523～0.86 μM）；而在BH-H61沉积柱下部（300 cmbsf以下），Ba2+浓度随 深度逐渐升高，在沉积柱底部达到峰值（9.31 μM）。

图 2 研究区BH-H75、BH-H13Y、BH-H61站位沉积物孔隙水中阴离子组分、CH4、DIC及其δ13CDIC的深度剖面图Fig.2 Depth profile of anion components, methane, DIC and δ13CDICvalues in sediment pore water of BH-H75,BH-H13Y and BH-H61 sites

4 讨论

4.1 OSR与SD-AOM硫酸盐消耗比例

孔隙水中（△DIC+△Ca2++△Mg2+）/△比率图解是识别研究区沉积物孔隙水中OSR和SDAOM最常用的方法[13-14,19,23-24]。OSR和SD-AOM消耗1 mol硫酸盐的同时，分别产生2和1 mol的碳酸氢根（反应式见式（1）和式（2）），即二者碳酸氢根产生量与硫酸盐消耗量的比值分别为2∶1和1∶1。因此，可根据孔隙水中碳酸氢根产生量与硫酸盐消耗量的不同比值来评估OSR和SD-AOM消耗硫酸盐的相对贡献。由于孔隙水中Ca2+、Mg2+离子也会消耗部分碳酸氢根，所以计算总碳酸氢根产生量时，还需考虑Ca2+和Mg2+离子消耗的碳酸氢根总量。首先，根据孔隙水中、DIC、Ca2+和Mg2+浓度与典型海水浓度（：28.9 mM；DIC：2.1 mM；Ca2+：10.3 mM；Mg2+：53.2 mM）的差异可以分别计算出△SO24−、△DIC、△Ca2+和△Mg2+。随后根据（△DIC+△Ca2++△Mg2+）/△SO42−比值来评估OSR和SD-AOM的相对比例。3个沉积柱的（△DIC+△Ca2++△Mg2+）/△SO42−比值见图4。在BH-H75的20～320 cmbsf深度内，数据在2∶1和1.5∶1实线内分布，表明在此深度内OSR占主导；在340～400 cmbsf深度间隔内，数据在1.5∶1和1∶1线内分布，表明靠近沉积柱底部可能受SD-AOM的影响。在BH-H13Y站位的20～195 cmbsf深度内，数据在2∶1和1.5∶1实线内分布，表明在此深度内OSR占主导；而在220～400 cmbsf深度内，数据在1.5∶1和1∶1线内分布（主要围绕1.25∶1虚线分布），表明靠近沉积柱底部SD-AOM占主导地位（约75%）。BH-H61的数据变化趋势呈现一定规律性。在BH-H61的20～340 cmsbf深度内，数据靠近1.25∶1虚线分布，显示SDAOM占主导（约75%）；在360～520 cmsbf深度内，数据在1.25∶1和1.1∶1虚线内分布，显示SD-AOM占主导（平均82.5%）。

图 3 研究区BH-H75、BH-H13Y和BH-H61沉积物孔隙水中部分碱土金属元素浓度（Ca2+、Mg2+、Sr2+、Ba2+）、Sr/Ca与Mg/Ca比值深度剖面图Fig.3 Depth profiles of partial alkaline-earth metal elements （Ca2+, Mg2+, Sr2+, Ba2+） and Sr/Ca and Mg/Ca ratios of sediment pore water in BH-H75, BH-H13Y and BH-H61 sites

图 4 孔隙水中硫酸盐的消耗量与DIC的产生量（Ca2+、Mg2+离子校正后）的相关关系图图中2∶1、1∶1的实线分别代表OSR和SD-AOM的贡献为100%，在这两个实线内，表示这两个过程均发生。其中1.5∶1实线表示OSR和AOM各占50%；1.25∶1和1.1∶1虚线表示SD-AOM的贡献分别占75%和90%，OSR仅占25%和10%。Fig.4 Relationship between the consumption of sulfate and the production of DIC （After correction by Ca2+and Mg2+ions）The solid lines of 2∶1 and 1∶1 in the figure represent the 100%contribution of OSR and SD-AOM respectively. Within these two solid lines, It indicates both processes occur. The 1.5∶1 solid line indicates OSR and AOM account for 50% respectively;The dashed lines of 1.25∶1 and 1.1∶1 indicate the contribution of SD-AOM accounts for 75% and 90%,while OSR accounts for only 25% and 10%.

表 2 BH-H75、BH-H13Y和BH-H61站位沉积物孔隙水中甲烷浓度、阴离子（SO42−、Cl−）、主微量元素（Na+、K+、Ca2+、Mg2+、Sr2+、Ba2+）、DIC和δ13CDIC、Sr/Ca与Mg/Ca比值特征Table 2 Features of methane concentration, anions （SO42−、Cl−）, major trace elements （Na+、K+、Ca2+、Mg2+、Sr2+、Ba2+）, DIC and δ13CDIC, Sr/Ca and Mg/Ca ratios of sediment pore water in BH-H75, BH-H13Y and BH-H61 sites

续表 2

此外，孔隙水中δ13CDIC变化特征也能定性辨别OSR和SD-AOM[9,24]。研究区DIC的可能来源包括：①底部海水向下扩散或原位储存于孔隙水中的DIC；②有机质降解过程产生的DIC；③AOM产生的DIC；④产甲烷作用残余的DIC向上迁移混合。这些不同来源DIC的δ13CDIC值不尽相同，因此，可用δ13CDIC值的变化特征识别产生DIC的生物地球化学过程。本研究中，3个沉积柱的δ13CDIC均随深度逐渐降低。BH-H75的δ13CDIC值在沉积柱底部降

低至−15.69‰，在海水（δ13C约为0）和南海沉积有机质碳同位素值（δ13C约为−20‰）[50]范围内，表明BH-H75可能主要受OSR影响，且底部受SDAOM的影响较小。δ13CDIC在BH-H13Y沉积柱底部降低至−28.4‰，轻于沉积有机质的δ13C值（约−20‰），推测其DIC可能来自于深部SD-AOM、浅部OSR混合。在BH-H61底部的480 cmbsf处，δ13CDIC降低至最小值（−35.07‰），结合（△DIC+△Ca2++△Mg2+）/△SO42−比率图解结果（图4），认为孔隙水中DIC主要来源于SD-AOM。然而在BH-H61的480 cmbsf深度以下，δ13CDIC值呈轻微正偏（−35.07‰～−23.87‰），这种现象可能与底部微生物产甲烷作用有关。由于产甲烷作用的微生物优先消耗反应物二氧化碳或乙酸中轻的12C同位素，导致残余DIC中富集13C同位素[51]，这种富13C同位素的DIC向上迁移混合，使得上覆沉积物孔隙水中δ13CDIC值正偏[28]。BH-H61站位的地温梯度为36.1 K/km（海底温度约2.6 ℃）（表1），考虑到微生物介导的产甲烷过程发生的温度低于50 ℃[5,52]，因此，微生物产甲烷作用在BH-H61底部发生的可能性 较 大。总 之，（△DIC+△Ca2++△Mg2+）/△SO42−比率图解与δ13CDIC值深度剖面特征显示在BHH13Y沉积柱中OSR和SD-AOM共同作用，而在BH-H75沉积柱中OSR占主导，在BH-H61沉积柱中SD-AOM占主导，且其底部可能存在微生物产甲烷作用。

4.2 SMTZ深度与甲烷通量的估算

SMTZ深度主要由上升的甲烷通量控制，一般相对较高的甲烷通量对应较浅的SMTZ深度，反之则相反[53-54]。SMTZ深度计算一般根据最小二乘线性回归的方法进行拟合，再外推硫酸盐降低为0 mM时对应的深度[13-14,21,31]，但前提是沉积柱是AOM主导的生物地球化学过程。因此，计算SMTZ深度时，不考虑研究区OSR占主导的BH-H75沉积柱以及BH-H13Y沉积柱上部（20～195 cmbsf）。在BH-H13Y沉积柱下部（220～400 cmbsf）SD-AOM主导，根据线性回归的方法推断其SMTZ深度约为700 cmbsf（图5）。此外，浓度、DIC浓度最大值和δ13CDIC最小值也可用于判定SMTZ深度[55]。在BH-H61中，浓度的变化与DIC浓度最大值和δ13CDIC最小值具有很好的一致性，从而推测其SMTZ深度为480 cmbsf（图2，图5），这与Feng等在北康盆地BHCL29站位观察的SMTZ的深度（530 cmbsf）接近[21]。但总的来说，北康盆地研究区观察的SMTZ深度（480和700 cmbsf）均浅于邻近海域南沙海槽的SMTZ深度（800～1 100 cmbsf）[56]。由于取样过程的降压和去气，测试的甲烷浓度远低于真实浓度（图2），所以用测试的甲烷浓度计算甲烷通量十分不准确。但是SD-AOM消耗硫酸盐与甲烷的比例为1∶1，可用硫酸盐通量近似代替进入SMTZ中的甲烷通量。值得注意的是，BH-H61和BH-H13Y沉积柱上部的硫酸盐通量远低于下部（图5），可能指示向上迁移的甲烷通量在时间上逐渐增强，从而反映一种非稳态环境下甲烷通量的变化模式[57]。BHH61上部（20～340 cmbsf）硫酸盐向下扩散的通量为30 mmol•m−2•a−1，且SD-AOM消耗硫酸盐比例为75%，推算早期甲烷渗漏通量为22 mmol•m−2•a−1（图5）；BHH61下部（360～520 cmbsf）硫酸盐向下扩散的通量为72.4 mmol•m−2•a−1，SD-AOM消耗硫酸盐比例为82%，推算现今甲烷渗漏通量为60 mmol•m−2•a−1（图5），比早期甲烷渗漏通量增加了38 mmol•m−2•a−1。BHH13Y上部的硫酸盐向下扩散通量为15 mmol•m−2•a−1，但这些硫酸盐主要被有机质降解消耗，表明前期的微弱甲烷渗漏对孔隙水化学组成的影响较少；在该沉积柱下部（220～400 cmbsf深度间隔内）硫酸盐向下扩散的通量为40 mmol•m−2•a−1，且SD-AOM对硫酸盐消耗的贡献为75%（图4），因此，估算现今的甲烷渗漏通量为30 mmol•m−2•a−1（图5）。本文估算的甲烷通量变化范围为22～60 mmol•m−2•a−1，与Feng等[21]在北康盆地研究区通过反应-传输模型模拟的结果接近（34.6～62.7 mmol•m-2•a−1）。南海北部的东沙海域（35～59 mmol•m−2•a−1）[14,32]和南海北部西沙海域麻坑区（14 mmol•m−2•a−1）[13]以及布莱克海台（20 mmol•m−2•a−1）[58]的计算数据也与本文结果具有很好的一致性。但是，本文估算的甲烷通量低于墨西哥湾海域（20～200 mmol•m−2•a−1）[59]和挪威中部G11麻坑区（300～540 mmol•m−2•a−1）[24]，可能指示相对较低的甲烷渗漏环境。综上所述，BHH61和BH-H13Y沉积柱中记录的较浅的SMTZ深度、上升的DIC浓度和强烈负偏的δ13CDIC值指示研究区存在甲烷渗漏活动，并且，向下变陡的硫酸盐浓度梯度表明甲烷渗漏强度随时间增强。

4.3 自生碳酸盐岩形成与重晶石溶解的孔隙水地球化学证据

自生碳酸盐岩沉淀能有效减少沉积物向上覆海水中释放的碳通量，在全球海洋碳循环中发挥重要作用[18,60]。有机质降解（产甲烷作用除外）和AOM产生的DIC提高孔隙水中碱度的同时，还能促进自生碳酸盐岩形成[61]。研究区孔隙水DIC浓度的升高（图2）对应Ca2+、Mg2+、Sr2+浓度的降低（图3，表2），表明可能有自生碳酸盐岩沉淀生成。研究区孔隙水中Mg2+/Ca2+摩尔比值随深度逐渐升高（图3），其中在BH-H61升高趋势最为显著，在SMTZ深度均达到最大值，可能指示有低Mg2+/Ca2+比值（＜4.86）的自生碳酸盐岩矿物形成，包括高镁方解石和文石（高镁方解石的Mg2+/Ca2+摩尔比为0.28；文石的Mg2+/Ca2+摩尔比为0.00125）[61-62]。沉淀的碳酸盐矿物种类也可用于指示沉积物中甲烷的传输行为[18]：文石倾向形成于高硫酸盐浓度和中—强甲烷渗漏环境，发育于近沉积物-水界面的沉积物中；而高镁方解石则倾向形成于低硫酸盐浓度和低—中等强度的甲烷渗漏环境，沉淀在SMTZ附近及以上沉积物孔隙中，孔隙水中的高硫酸盐浓度会对其形成产生抑制作用[63]。此外，孔隙水中Sr2+/Ca2+比值变化特征可进一步识别自生碳酸盐岩的矿物种类[25]。Sr2+/Ca2+摩尔比随深度逐渐升高（图3），表明在研究区可能有低Sr2+/Ca2+比值（＜0.011）的高镁方解石矿物形成（高镁方解石的Sr2+/Ca2+摩尔比为0.003）[62]。图6展示了孔隙水中Mg/Ca与 Sr/Ca相关关系图，可看出研究区孔隙水数据靠近高镁方解石沉淀线一侧分布，数据变化趋势与高镁方解石沉淀线以及前人数据趋势一致[62]，进一步指示可能存在高镁方解石沉淀。研究区沉积物孔隙水中出现高镁方解石沉淀的证据，而缺乏文石沉淀的证据，粗略指示研究区的甲烷渗漏可能处于低—中等强度，这与本研究甲烷通量的计算结果相符。

图 5 BH-H13Y、BH-H61的硫酸盐浓度深度剖面图和SMTZ深度与甲烷通量Fig.5 Depth profiles of sulfate concentrations, the depth of SMTZ and the methane flux in BH-H13Y and BH-H61 sites

一般地，海洋沉积物中的钡元素少量以自生重晶石形式存在，大部分主要赋存于铝硅酸盐中[64]。因为重晶石和铝硅酸盐中硫酸钡溶解度低，所以海水及沉积物孔隙水中的Ba2+浓度也较低。然而，SDAOM对孔隙水中硫酸盐的消耗促进重晶石的溶解，从而释放Ba2+离子进入孔隙水。尤其在硫酸盐浓度严重亏损的SMTZ中，大量Ba2+离子被释放进入SMTZ层位，并在SMTZ深度内形成显著的Ba峰[49,54]。在研究区，BH-H75和BH-H13Y沉积柱以及BH-H61沉积柱上部（20～280 cmbsf），Ba2+浓度变化趋势不明显（0.436～0.86 μM），可能仅继承了原位海水的浓度特征。然而在BH-H61沉积柱下部（300 cmbsf以下），Ba2+浓度随深度逐渐升高，在柱子底部达到最大值9.31 μM（图3），表明重晶石可能在SMTZ附近及下部出现了溶解作用。

图 6 研究区孔隙水中Mg/Ca- Sr/Ca图中两条实线表示当孔隙水中形成文石或者髙镁方解石时，孔隙水相对于海水Mg/Ca和Sr/Ca的变化关系，灰色菱形为前人的数据[62]，绿色菱形代表研究区3个沉积柱的数据。Fig.6 The weight ratio of Mg/Ca vs. Sr/Ca in sediment pore water of the study areaTwo solid lines indicate the change relationship of pore water Mg/Ca and Sr/Ca with respect to that of seawater, when aragonite or high Mg-calcite is formed in sedimentary pore water. Grey diamonds in the figure are theprevious data[62]. Green diamonds represent the data of the three sedimentarycolumns of the study area.

5 结论

（1）BH-H75、BH-H13Y、BH-H61沉积柱中均发生了硫酸盐还原反应（OSR或SD-AOM），但每个沉积柱中OSR和SD-AOM的贡献存在差异。（△DIC+△Ca2++△Mg2+）/△SO42−比值和DIC的13C同位素特征表明在BH-H75沉积柱中OSR占主导，在BH-H61沉积柱中SD-AOM占主导，在BH-H13Y站位，浅部发生了OSR，深部发生了SD-AOM。

（2）根据线性拟合硫酸盐浓度及DIC浓度最大值、δ13CDIC最小值特征推测BH-H13Y和BH-H61的SMTZ深度约为700和480 cmbsf。BH-H61和BHH13Y沉积柱中，较浅的SMTZ深度、上升的DIC浓度和强烈负偏的δ13CDIC值指示研究区存在甲烷渗漏活动。BH-H61和BH-H13Y的硫酸盐浓度梯度向下逐渐变陡，可能表明甲烷渗漏强度随时间增强。根据硫酸盐扩散通量估算，BH-H61沉积柱上部记录的早期甲烷渗漏通量为22 mmol•m−2•a−1；下部记录的现今甲烷渗漏通量为60 mmol•m−2•a−1，较早期增加了38 mmol•m−2•a−1，可能显示甲烷渗漏强度随时间逐渐增强的趋势。BH-H13Y沉积柱上部的硫酸盐主要被有机质降解所消耗，下部记录的甲烷渗漏通量为30 mmol•m−2•a−1。

（3）孔隙水中部分碱土金属离子浓度变化特征及Mg/Ca与 Sr/Ca比值指示在3个沉积柱中存在自生碳酸盐岩矿物的沉淀，且以高镁方解石沉淀为主，与根据硫酸盐扩散通量估算的低—中等强度甲烷渗漏通量相符。在BH-H61沉积柱下部，Ba2+浓度变化特征显示重晶石可能在SMTZ附近及下部出现了溶解作用。