南海天然气水合物智能识别方法与应用

田冬梅, 杨胜雄, 刘鑫, 李沅衡, 胡广, 曹荆亚, 周军明, 邓雨恬

田冬梅,杨胜雄,刘鑫,等. 南海天然气水合物智能识别方法与应用[J]. 海洋地质与第四纪地质,2024,44(6): 25-33. DOI: 10.16562/j.cnki.0256-1492.2024092401
引用本文: 田冬梅,杨胜雄,刘鑫,等. 南海天然气水合物智能识别方法与应用[J]. 海洋地质与第四纪地质,2024,44(6): 25-33. DOI: 10.16562/j.cnki.0256-1492.2024092401
TIAN Dongmei,YANG Shengxiong,LIU Xin,et al. Intelligent identification and application of gas hydrate in South China Sea[J]. Marine Geology & Quaternary Geology,2024,44(6):25-33. DOI: 10.16562/j.cnki.0256-1492.2024092401
Citation: TIAN Dongmei,YANG Shengxiong,LIU Xin,et al. Intelligent identification and application of gas hydrate in South China Sea[J]. Marine Geology & Quaternary Geology,2024,44(6):25-33. DOI: 10.16562/j.cnki.0256-1492.2024092401

南海天然气水合物智能识别方法与应用

基金项目: 国家自然科学基金国家重大科研仪器研制项目“海底地震与电磁同步探测系统关键技术及验证样机”(42327901); 国家自然科学基金项目“南海北部高富集天然气水合物储层特征与成藏控制机理研究”(U2244224);广州市基础与应用基础研究项目“基于地震属性海域天然气水合物识别方法研究”(2023A04J0916)
详细信息
    作者简介:

    田冬梅(1995—),女,博士,主要从事海洋地球物理、天然气水合物识别计算研究,E-mail:dongmeitian@gmlab.ac.cn

    通讯作者:

    刘鑫(1983—),男,硕士,高级工程师,主要从事海洋地质环境相关研究,E-mail:lxin05@gmlab.ac.cn

  • 中图分类号: P714

Intelligent identification and application of gas hydrate in South China Sea

  • 摘要:

    天然气水合物是一种重要的能源资源,具有能量高、储量大、分布广和埋藏浅等优势。准确地从地层中识别出天然气水合物储层是应用天然气水合物资源的必要前提。本文围绕水合物勘探识别的难点问题,结合海洋-地质-人工智能学科交叉技术,以具有显示度的地球物理属性参数为基础,研究并提出了有效的含水合物地层识别技术方法,在中国南海东沙海域研究区进行了方法的验证,选择了几种较为常用的机器学习算法,例如随机森林、Bagging、AdaBoost、和最近邻(KNN)算法,对水合物变化灵敏度较高的纵波速度和密度属性进行数据分析,通过训练优化不同算法模型参数,对比不同算法模型的识别分类效果。结果表明,这几种算法都能够较好地对地层中是否含有水合物进行区分,其中KNN性能最好,表明基于机器学习手段能够提高天然气水合物的识别精度和准确性。

    Abstract:

    Gas hydrate is an important ideal energy source, with advantages of high energy, large reserves, wide distribution, and shallow burial. Accurate identification of gas hydrate reservoirs and estimation of hydrate saturation are the prerequisite for the application of gas hydrate resources. This study focuses on the difficult issues of hydrate identification, combining the interdisciplinary technologies of oceanology, geology, and artificial intelligence. Effective methods of hydrate-bearing strata identification were proposed based on the geophysical attributes, and verified in the Dongsha area of South China Sea. Machine-learning algorithms were used to analyze whether the sediment contains gas hydrates. Several commonly used machine-learning algorithms were selected, including random forest, Bagging, AdaBoost, and KNN; and data were analyzed based on the P-wave velocity and density attributes that are more sensitive to hydrate existence. The parameters of different algorithms were trained and optimized, and the effects of different algorithms on the identification and classification were compared. All these algorithms could do good on whether there is hydrate in the sediment, of which KNN algorithm was shown the best. Therefore, machine-learning-based methods could improve the identification accuracy of gas hydrate.

  • 海底冷泉活动中广泛发育自生碳酸盐岩[1], 其往往指示强度很高、同时伴随沉积环境和氧化还原条件剧烈变化的流体渗漏活动[1-2]。不同的沉积条件下形成的冷泉碳酸盐岩有着不同的形态学、矿物学和地球化学特征[3-5],可以很好地反映沉积时的冷泉流体类型、活动强度、流体来源及生物地球化学信息[6-11],是研究深海冷泉环境必不可少的研究对象。冷泉碳酸盐岩烟囱以其独特的流体管道和圆柱形的形貌成为冷泉流体渗漏活动最显著的证据之一,使其受到了特别关注,包括成岩环境、流体来源以及成岩模式的研究都较多[12-17]

    中国南海北部广泛发育冷泉区[18-30],其中不同海域的冷泉碳酸盐岩的碳、氧同位素组成变化范围较大, 反映了冷泉流体的来源、渗漏速率及沉积环境均有较大差异[2]。南海东沙东北陆坡发育大量冷泉碳酸盐岩烟囱[31-36],这些烟囱的形成除了与海底天然气水合物分解释放的流体有关外,其烟囱内外层具有明显不同的矿物学和地球化学特征,呈现由外向内的生长模式,表明冷泉中渗漏流体扩散作用的影响,这使得烟囱状冷泉碳酸盐岩成为了记录时间和空间上冷泉流体渗漏组成和强度等特征变化的良好载体[32-34]

    尽管冷泉碳酸盐岩烟囱样品中渗漏流体强度变化对碳、硫、氧化还原敏感的微量元素的影响已经有相关研究[32-35],但目前对这种强渗漏环境下,渗漏流体强度变化对铁(Fe)和磷(P)元素循环及相关生物地球化学过程的影响还不清楚。在海底甲烷渗漏环境中,碳-铁耦合的相关生物地球化学过程可能会驱动磷元素循环过程,即Fe驱动的甲烷厌氧氧化过程有可能导致了含磷矿物蓝铁矿的形成[37-44]。前人对南海块状冷泉碳酸盐岩的磷组分和磷循环进行了分析和研究[42],但块状冷泉碳酸盐岩样品相比烟囱状冷泉碳酸盐岩样品而言,无法提供甲烷渗漏流体强度变化信息,导致无法探讨渗漏流体强度变化对Fe和P元素循环和相关生物地球化学过程的影响。

    因此,南海东沙海域广泛发育的烟囱状冷泉碳酸盐岩样品是研究不同流体渗漏强度对不同磷组分和铁-磷耦合影响的绝佳材料。本文将以南海东沙海域典型的烟囱状冷泉碳酸盐岩为研究对象,进行相应的磷组分和铁组分的分步提取和其他相关基本地球化学参数分析,鉴别不同形态磷和铁的来源和特征,阐明不同形态磷和铁可能参与的生物地球化学循环过程,揭示不同甲烷渗漏强度环境中甲烷渗漏活动对不同磷组分和铁-磷耦合关系影响,从而加深甲烷渗漏活动对海底磷、铁循环影响的认识。

    本文将以由“海洋四号”科考船于2015年在九龙甲烷礁“Site 3”通过海底拖网采集获得的南海东沙海域典型的烟囱状冷泉碳酸盐岩为研究对象,将采集到的冷泉碳酸盐岩样品用切割机沿着烟囱横截面切割,总体上烟囱样品直径约10 cm,外层呈现的颜色要比内层颜色深,内层烟囱孔洞发育。沿横截面从外向内对冷泉碳酸盐岩烟囱用牙钻取10个样品(图1),进行相应的磷组分和铁组分的分步提取和其他相关地球化学分析。

    图  1  采样位置和样品情况
    a: 位于东沙海域的采样站位九龙甲烷礁site 3;b: 冷泉碳酸盐岩烟囱样品横切面和取样点S1到S10。冷泉碳酸盐岩烟囱的横截面直径约10 cm,用牙钻间隔约0.8 cm取一个点,共取10个点。
    Figure  1.  Sampling location and subsamples
    a: The sampling site located in the Dongsha sea area, known as Methane Reef Site 3; b: the subsample of seep carbonates chimneys from S1to S10. The cross-section diameter of the seep carbonates chimneys is about 10 cm, and one point was taken about 0.8 cm apart for 10 points in total.

    冷泉碳酸盐岩烟囱物相分析(XRD)在上海海洋大学完成。用德国 BRUKER D8 ADVANCE X 射线衍射仪进行分析,工作参数为铜靶 Kα射线,功率为 40 kV,30 mA,扫描范围 2θ = 3°~85°,入射狭缝为 1 mm,扫描速度为 4°/min。矿物百分含量以d[104]峰作为校准曲线计算方解石、镁方解石的比例,以d[111]峰高计算文石的相对含量。通过d[104]峰高计算碳酸盐岩矿物中镁的含量。

    冷泉碳酸盐岩烟囱中碳和硫含量在上海海洋大学Elementar型元素分析仪测定。用超精密天平取其中约 25 mg 粉末样品进行上机测试。得到的数据是冷泉碳酸盐岩烟囱中总碳(TC) 和总硫(TS)含量。采用正磷酸法测定碳酸盐岩中碳酸盐相的碳氧同位素组成,具体做法在 75 ℃ 条件下,利用正磷酸与样品反应释放的CO2,通过Thermo Finnigan Delta V Advantage稳定同位素质谱仪对释放的CO2进行分析获得碳氧同位素组成。以 VPDB 国际标准为参考标准,δ13C 和 δ18O分析精度均优于±0.1‰。

    参照经典文献和最新文献[45-47]中所示方法进行冷泉碳酸盐岩中不同铁组分化学连续提取。不同铁组分提取步骤和相应条件如下:①碳酸盐相的铁(FeCarb),约50 mg样品与1 M乙酸钠溶液反应24 h,pH=4.5;②易还原铁氧化物的铁(FeOx1),上一步离心后的残渣与溶于体积分数为25%乙酸的1 M 盐酸羟胺溶液反应 48 h;③可还原性铁氧化物的铁(FeOx2),上一步离心后的残渣与 50 g/L 连二亚硫酸钠溶液反应2 h,pH=4.8(用0.35 M乙酸/0.2 M乙酸钠调节pH值);④磁铁矿的铁(FeMag),上一步离心后的残渣与 0.2 M 草酸铵/0.17 M 草酸混合溶液反应 6 h,pH=3.2。上述操作均在室温下进行,每步提取完成后离心并用0.2 μm滤膜过滤取上清液待测。FeCarb、FeOx1 和 FeOx2含量测试按照Alcott 等[47]所示方法,使用分光光度法进行测试,FeMag的含量用电感耦合等离子体发射光谱仪(ICP-OES)测试。实验室内部标样和平行样用于整个提取实验的质量监控,结果显示不同铁组分分析的相对误差通常小于5%。

    冷泉碳酸盐样品中不同形态磷的提取及含量测定参照化学连续提取法[48](即SEDEX法)进行。约50 mg样品依次与相应提取液反应,主要步骤如下:①可交换态磷(PEx)提取。使用1 M MgCl2提取液(pH=8),提取0.5 h;②铁氧化物结合态磷(PFe)提取。使用0.3 M柠檬酸钠和25 g/L连二亚硫酸钠作为提取液(pH=7.5),反应8 h。过滤后的样品接着用1 M 氯化镁提取液(pH=8)清洗半小时并离心过滤取上清液;③自生磷酸盐的磷(PAuth)提取。使用1 M醋酸钠-醋酸缓冲液提取(醋酸调节pH至4),反应6 h。接着过滤的沉积物同样用1 M 氯化镁提取液(pH=8)清洗半小时并离心过滤取上清液;④陆源碎屑磷(PDetr)提取。使用1 M盐酸溶液提取,反应24 h;⑤有机磷(POrg)提取。为了防止化学连续提取过程中POrg的损失,采用差值法进行POrg的提取。具体步骤包括:(a)冷泉碳酸盐岩样品与1 M盐酸溶液反应24 h;(b)样品经过550 ℃灼烧2 h后,与1 M盐酸溶液反应24 h。上述步骤(b)与(a)磷含量的差值即POrg。上述所有提取实验均在室温下进行,其中,PEx和PFe的提取应在无氧的、充满氮气的手套箱中进行操作,以防止黄铁矿的氧化对PEx和PFe含量的干扰[43-44]。采用分光光度法测试提取液中PEx、PAuth、PDetr、POrg的含量。为了使得提取液中PFe含量能够用分光光度法测试、消除提取液对分光光度法的干扰,按照Huerta-Diaz等[49]前处理方法,对含PFe的提取液使用Mg(OH)2共沉淀法消除干扰,然后用弱碱性溶液清洗干净后用盐酸溶解,最后用分光光度法测试PFe含量。根据实验室内部标样和平行样结果,不同磷组分含量分析相对误差一般小于10%。

    冷泉碳酸盐岩烟囱的矿物成分主要有4种,包括石英、长石、方解石和文石(图2)。结合取样点,可以把样品分成内外两层,其中编号S1—S3和S8—S10为烟囱外层,编号S4—S7为烟囱内层。除内层的S5和S7碳酸盐岩矿物主要为文石外(文石含量在70%以上),其余编号的碳酸盐矿物主要由方解石组成,方解石含量为70.2%~84.1%。烟囱内外层都含有较高的陆源碎屑矿物,主要包括石英、长石,含量为 10.7%~29.8%,总体上来看,大致具有外层陆源碎屑矿物含量普遍高于内层的规律。

    图  2  冷泉碳酸盐岩烟囱的矿物组成
    Figure  2.  Mineral composition of seep carbonate chimneys

    冷泉碳酸盐岩烟囱样品总碳含量为 0.21%~11.70%,变化较大,平均含量约7.62%(图3a)。随外层向内层有微弱的减少趋势,S5点总碳含量最低,而S2出现了相对较高的总碳含量。样品中总硫(TS)含量变化较小,平均含量约0.11%,其中S1出现了极高的总硫含量,S5和S10 出现了相对较高的总硫含量,总体来看硫含量偏高。

    图  3  冷泉碳酸盐岩烟囱样品的碳、硫含量(a)以及碳和氧稳定同位素组成(b)
    Figure  3.  Carbon and sulfur content (a), and carbon and oxygen stable isotope compositions (b) of the seep carbonate chimney samples

    冷泉碳酸盐岩烟囱的碳氧同位素变化范围较大(图3b)。δ13C 为−61.73‰~−55.53‰,平均值约为−58.78‰,δ18O 为5.45‰~5.83‰,平均值约为5.61‰。其中编号S5、S6、S7的碳同位素组成较其他编号明显偏负,而氧同位素组成则呈波动性,编号S2、S6、S10的氧同位素组成明显偏正。

    碳酸盐岩烟囱样品各铁组分含量分析结果见图4a。结果显示,FeCarb含量总体变化不大,内层略低,其中S5处FeCarb含量仅有0.02 mg/g,而烟囱外层FeCarb平均含量为3.97 mg/g。FeMag含量从烟囱内层到外层逐渐升高,其变化范围为0.05~0.90 mg/g,平均为0.58 mg/g,在S5处含量最低。FeOx含量为FeOx1 和 FeOx2含量之和,其含量靠近烟囱内层同样逐渐减少,S5和S7处FeOx含量很低,其中S7处为0.80 mg/g,外层平均5.99 mg/g。由于冷泉碳酸盐岩样品中的硫一般由黄铁矿的硫构成,因此黄铁矿的铁(Fepy)含量可以根据样品中总硫含量进行粗略估算。Fepy含量除S1、S5和S10以外7个点含量大致相同,平均为1.01 mg/g。而S1处测出含量值8.94 mg/g,明显高于其他点位,S5和S10的含量分别为1.94 mg/g 和1.76 mg/g,略高于另外7个点位。整体来看,Fepy含量外层外围高于内层,向内含量逐渐降低,到内层中心又重新增加。

    图  4  冷泉碳酸盐岩烟囱的不同铁和磷组分含量
    Figure  4.  Contents of different iron and phosphorus component of seep carbonates chimneys

    碳酸盐岩烟囱样品各取样点铁组分相对含量分析结果见图5a。结果显示, FeCarb内层占比略高于外层,但是在S5 处出现极低值(0.995%);FeMag占比内外层相当,都较低;FeOx1内外层差异不大,S5未测出FeOx1,占比为0;FeOx2内层低于外层,向内占比逐渐降低,特别是S5和S7均为0;FePy内层占比高于外层,特别是S5处占比为96.5%,明显高于其他点位,S1(37.8%)和S10(34.9%)也较其他点位占比高。整体而言,S1、S5和S7三点FePy占比更高,其他点位FeOx1和FeCarb占比更高。以FePy与烟囱样品高活性铁FeHR(FeHR = FeCarb + FeOx1 + FeOx2 + FeMag + FePy)比值表征样品黄铁矿化程度,结果显示,烟囱外围黄铁矿化度较高,向内有所降低,内层又有所增高,特别是S5,FePy/FeHR值为0.964。

    图  5  冷泉碳酸盐岩烟囱样品各铁组分和磷组分相对百分占比
    a:铁组分相对百分占比,b:磷组分相对百分占比。
    Figure  5.  Relative proportions of iron and phosphorus component in the seep carbonate chimney samples
    a: Relative proportions of the iron component, b: relative proportions of the phosphorus component.

    烟囱样品中各取样点磷组分含量见图4b。POrg平均为0.57 µmol/g,内外层差异不大。PEx含量外层比内层稍高,S5和S7点PEx含量最低,为0.92 µmol/g。PDetr含量内层明显比外层低,外层平均含量为5.15 µmol/g,而S5处含量最低,仅为0.42 µmol/g。PFe和PAuth含量内层与外围含量相当,特别是S6处PAuth明显较其他点位大、出现极大值8.39 µmol/g。从样品的总无机磷(PInorg,PInorg = PEx + PFe + PAuth + PDetr)含量上看,外层含量相当,平均值为13.73 µmol/g。而内层含量则相对较少,特别是S5(8.87 µmol/g)和S7(9.31 µmol/g),含量明显较其他点位小。

    碳酸盐岩烟囱样品各取样点磷组分相对含量分析结果见图5b。结果显示,PEx和PFe内外层相当,占比均较低;PAuth占比内层高于外层,特别是S5和S6,PAuth占比约60%;而PDetr占比则外层高于内层,S5处PDetr占比最小,仅为4.33%。POrg占比均较低,个别点占比为0。整体而言,烟囱样品内外层磷组分占比中PAuth最高,其次为PDetr

    碳酸盐岩烟囱极端负偏的碳同位素组成(δ13CVPDB<−55.53‰,图3)表明甲烷的厌氧氧化作用导致了该碳酸盐的形成。由于样品的δ13C值均低于–50‰,可以推断形成烟囱的流体中甲烷来源为生物成因。根据前人在相同站位冷泉碳酸盐岩的年代学数据和当时的海底温度(3.3~4.1 ℃)和海水氧同位素组成(0.64‰~0.97‰)[32],假定本文冷泉碳酸盐岩样品也在同期形成,通过方解石/文石-水体系的氧同位素分馏方程计算它们与海水平衡时的氧同位素组成变化范围分别为2.5‰~2.7‰和3.2‰~3.4‰。对比发现冷泉碳酸盐岩烟囱表现出异常高的氧同位素组成(5.45‰~5.83‰),这说明其形成时有富18O 流体的加入。天然气水合物分解和黏土矿物脱水反应均有可能导致冷泉流体中的水富含18O。考虑到本研究区附近表层沉积物发育水合物[50],并且表层沉积物孔隙水中也没有黏土矿物脱水反应导致孔隙水浓度异常的报道[51-53],因此推测烟囱的18O富集主要来源于富含18O 水合物的分解。

    前人研究结果显示本研究区冷泉碳酸盐岩烟囱样品的生长是由外向内,即外层先形成,内层后形成[32-34]。冷泉碳酸盐岩烟囱外层具有较高含量的石英、长石等陆源碎屑物质以及铁氧化物以及碎屑磷(图24),说明其在早期形成过程中受到较大的背景沉积的影响,在形成过程中捕获了更多背景沉积物,这些陆源碎屑物质中的TOC可能会影响烟囱的δ13C值。另外,根据之前对本区冷泉碳酸盐岩烟囱的研究,烟囱由外层向内层形成的过程中,由于渗漏通道空间越来越小, 导致流体渗漏强度越来越大[31]。本文研究的碳酸盐岩烟囱内层S3和S5点位以文石为主,而其余样品均以方解石为主(图2),同时一般认为文石相比方解石形成时的环境中硫酸盐浓度更高[16, 54],因此推测本文碳酸盐岩烟囱形成过程中流体渗漏强度内层要高于外层。这种内层渗漏强度高的特征,会使得形成于内层的自生沉积矿物更容易受到甲烷流体的影响,导致本文和类似的烟囱状碳酸盐岩样品中δ13C组成普遍出现两边高、中心低的分布特征[34, 55]图3b)。

    在大陆边缘环境中,铁、磷的耦合循环是一个动态过程[39, 42-43, 56-61]。铁氧化物还原、有机质氧化或甲烷厌氧氧化等过程中铁氧化物伴生的磷会被释放出来。碎屑磷在沉积于海底之前无法被利用,在埋藏后也不会发生相变[40]。烟囱样品磷组分提取结果显示碎屑磷含量内层比外层低(图5b),表明烟囱形成前期陆源碎屑输入较多,受到陆源碎屑输入的影响,可能导致陆源输入的铁氧化物从外层向内层含量降低。随着流体强度变化,冷泉渗漏环境中较强的微生物硫酸盐还原作用,不仅导致更加还原的环境和更多活性铁(主要是铁氧化物)的消耗,影响了烟囱内外层黄铁矿化度,也抑制了Fe的异化还原[33]。这些叠加因素可能最终导致了碳酸盐岩烟囱样品中铁氧化物含量内层低于外层。同时,铁氧化物的还原也会伴随其吸附无机磷的释放[39],也有可能会影响样品中磷的分布。

    由冷泉碳酸盐岩烟囱磷组分占比(图5b)可知,PAuth所占比例最大,其次是PDetr、PFe、PEx和POrg。烟囱样品中磷的分配表明,流体中最主要的一部分磷以自生磷酸盐(PAuth)的形式沉淀。磷从不稳定形态转移到PAuth,被认为是导致磷在全球海洋沉积物中埋藏的主要过程[62]。PAuth的形成是海洋沉积物中普遍存在的现象,但其具体机制尚不清楚[42]。冷泉渗漏流体中的磷来源于硫酸盐-甲烷过渡带下方,并在自生碳酸盐岩矿物中被化学吸附[42]。硫酸盐驱动的甲烷厌氧氧化可通过增加冷泉渗环境中孔隙水的碱度,导致沉积物或海底自生钙碳酸盐岩的沉淀[39,63]。这些形成的碳酸钙可能是自生PAuth沉淀的成核位点[42]。冷泉碳酸盐岩烟囱中PFe/FeOx比值内层明显高于外层(图6),显示冷泉碳酸盐岩烟囱内层更富集PFe,表明内层有额外的磷加入。PFe的含量内层稍高,可能与铁氧化物吸附磷和铁(II)磷矿物有关,而铁(II)磷矿物主要为蓝铁矿类矿物[39],在铁的连续提取过程在第②步被提取,从而被包含在化学提取的PFe中。结合碳酸盐岩烟囱内层FeOx含量低和相应吸附磷含量也低的特征[39],推测烟囱内层可能有蓝铁矿形成。烟囱样品S6中铁结合态磷含量低可能与后期蓝铁矿发生蚀变转变成自生磷灰石有关。综上所述,通过对烟囱状冷泉碳酸盐岩不同的磷和铁组分分析,发现在强渗漏条件下甲烷渗漏强度的改变会影响样品中铁结合态磷和自生磷灰石含量,更强的甲烷渗漏强度更有利于蓝铁矿和自生磷灰石的形成,从而对冷泉环境中的磷和铁循环产生影响。

    图  6  冷泉碳酸盐岩烟囱样品的晶质铁氧化物(FeOx2)和非晶质氧化铁含量(FeOx1)分别与铁结合态磷(PFe)的投点图
    图中红色点为冷泉碳酸盐岩烟囱内层样品,黑色点为冷泉碳酸盐岩烟囱外层样品。
    Figure  6.  The cross-plot of crystalline iron oxides (FeOx2) and amorphous iron oxide content (FeOx1) from cold seep carbonate chimney samples vs iron-bound phosphorus (PFe), respectively
    The red dots are seep carbonate chimney samples. The black dots are samples from the outer layers of seep carbonates.

    通过对南海东沙海域九龙甲烷礁site 3站位的冷泉碳酸盐岩烟囱进行矿物组成、碳氧同位素、碳和硫占比,以及不同铁组分和磷形态分析,确定了碳酸盐岩烟囱从外向内沉积过程中曾发生过强烈的、可能由天然气水合物分解导致的甲烷渗漏活动。但该甲烷渗漏活动并不稳定,且随着流体通路变窄甲烷渗漏活动增强。较强的甲烷渗漏活动会导致烟囱碳同位素值大幅负偏,氧同位素值轻微偏正,并会导致大量沉淀碳酸钙和黄铁矿的生成。在甲烷强渗漏条件下,烟囱内层明显富集铁结合态磷和自生磷灰石,说明更强的甲烷渗漏强度和更接近海底的环境有利于蓝铁矿和自生磷灰石的形成,进而影响冷泉环境中的磷和铁循环,但对于其中可能形成的蓝铁矿机制和微环境等还有待未来进一步研究。

    致谢:“海洋四号”科考船组织者、船员和科考人员在样品采集过程中给予的帮助。感谢审稿人提出的宝贵意见。

  • 图  1   全球天然气水合物分布图(据文献[3]修改)

    Figure  1.   Global gas hydrate distribution (modified from reference [3])

    图  2   研究区位置

    Figure  2.   Location of the study area

    图  3   GMGS2-05、GMGS2-08和GMGS2-16站位纵波速度和密度测井曲线

    天然气水合物区域用蓝色标记。

    Figure  3.   Logging curves of P-wave velocity and density at sites GMGS2-05, GMGS2-08 and GMGS2-16

    The gas hydrate areas are marked in blue.

    图  4   AdaBoost算法对GMGS2-05、GMGS2-08和GMGS2-16纵波速度和密度数据训练识别水合物结果

    Figure  4.   The hydrate identification results using AdaBoost algorithm based on the P-wave velocity and density data at sites GMGS2-05, GMGS2-08, and GMGS2-16

    图  5   随机森林算法对GMGS2-05、GMGS2-08和GMGS2-16站位纵波速度和密度数据训练识别水合物结果

    Figure  5.   The hydrate identification results using by random forest algorithm based on the P-wave velocity and density data at sites GMGS2-05, GMGS2-08 and GMGS2-16

    图  6   Bagging算法对GMGS2-05、GMGS2-08和GMGS2-16站位纵波速度和密度数据训练识别水合物结果

    Figure  6.   The hydrate identification result using Bagging algorithm based on the P-wave velocity and density data at sites GMGS2-05, GMGS2-08 and GMGS2-16

    图  7   KNN算法对GMGS2-05、GMGS2-08和GMGS2-16站位纵波速度和密度数据训练识别水合物结果

    Figure  7.   The hydrate identification result using KNN algorithm based on the P-wave velocity and density data at sites GMGS2-05, GMGS2-08 and GMGS2-16

    图  8   不同算法模型的ROC曲线

    Figure  8.   The ROC curves of different algorithms

    表  1   评价结果指标意义

    Table  1   Significance of evaluation results indicators

    预测为真预测为负
    实际为真TPFN
    实际为负FPTN
    下载: 导出CSV

    表  2   每个模型的超参数调整

    Table  2   Hyperparameter tuning for each algorithm

    监督模型超参数调整范围最优参数
    AdaBoostN estimators5~2515
    Learning rate0.1~0.90.5
    RFmax leaf nodes5~4510
    n estimators0~15030
    Baggingn estimators10~500100
    Max samples50~500150
    KNNN neighbors5~20010
    weightsuniform/distanceuniform
    下载: 导出CSV

    表  3   各算法评价指标计算结果

    Table  3   Calculation results of evaluation indicators of each algorithm

    Supervision algorithm Accuracy Precision Recall F1 AUC
    AdaBoost 0.9171 0.8625 0.4726 0.6106 0.9146
    Random forest 0.9388 0.8857 0.6370 0.7410 0.9278
    Bagging 0.9341 0.8725 0.6096 0.7177 0.9449
    KNN 0.9379 0.8636 0.6507 0.7422 0.9580
    下载: 导出CSV
  • [1]

    Sloan E D Jr, Koh C A. Clathrate Hydrates of Natural Gases[M]. 3rd ed. Boca Raton: CRC Press, 2007.

    [2]

    Chong Z R, Yang S H B, Babu P, et al. Review of natural gas hydrates as an energy resource: Prospects and challenges[J]. Applied Energy, 2016, 162:1633-1652. doi: 10.1016/j.apenergy.2014.12.061

    [3]

    Collett T S, Johnson A H, Knapp C C, et al. Natural gas hydrates-A review[M]//Collett T, Johnson A, Knapp C, et al. Natural Gas Hydrates—Energy Resource Potential and Associated Geologic Hazards. American Association of Petroleum Geologists, 2009, 89: 146-219.

    [4]

    Boswell R, Hancock S, Yamamoto K, et al. Natural gas hydrates: Status of potential as an energy resource[M]//Letcher T M. Future Energy: Improved, Sustainable and Clean Options for our Planet. 3rd ed. New York: Elsevier, 2020: 111-135.

    [5]

    Yang S X, Liang J Q, Lei Y, et al. GMGS4 gas hydrate drilling expedition in the South China Sea[J]. Fire in the Ice, 2017, 17(1):7-11.

    [6] 雷裕红, 宋颖睿, 张立宽, 等. 海洋天然气水合物成藏系统研究进展及发展方向[J]. 石油学报, 2021, 42(6):801-820

    LEI Yuhong, SONG Yingrui, ZHANG Likuan, et al. Research progress and development direction of reservoir-forming system of marine gas hydrates[J]. Acta Petrolei Sinica, 2021, 42(6):801-820.]

    [7]

    Zhang G X, Liang J Q, Lu J A, et al. Geological features, controlling factors and potential prospects of the gas hydrate occurrence in the east part of the Pearl River Mouth Basin, South China Sea[J]. Marine and Petroleum Geology, 2015, 67:356-367. doi: 10.1016/j.marpetgeo.2015.05.021

    [8]

    Liang J Q, Zhang W, Lu J A, et al. Geological occurrence and accumulation mechanism of natural gas hydrates in the eastern Qiongdongnan Basin of the South China Sea: Insights from Site GMGS5-W9-2018[J]. Marine Geology, 2019, 418:106042. doi: 10.1016/j.margeo.2019.106042

    [9]

    Sha Z B, Liang J Q, Zhang G X, et al. A seepage gas hydrate system in northern South China Sea: Seismic and well log interpretations[J]. Marine Geology, 2015, 366:69-78. doi: 10.1016/j.margeo.2015.04.006

    [10]

    Milkov A V. Global estimates of hydrate-bound gas in marine sediments: how much is really out there?[J]. Earth-Science Reviews, 2004, 66(3-4):183-197. doi: 10.1016/j.earscirev.2003.11.002

    [11] 勾丽敏, 张金华, 王嘉玮. 海洋天然气水合物地震识别方法研究进展[J]. 地球物理学进展, 2017, 32(6):2626-2635

    GOU Limin, ZHANG Jinhua, WANG Jiawei. Progress in seismic identification approach of marine gas hydrate[J]. Progress in Geophysics, 2017, 32(6):2626-2635.]

    [12] 宋海斌, 张岭, 江为为, 等. 海洋天然气水合物的地球物理研究(Ⅲ): 似海底反射[J]. 地球物理学进展, 2003, 18(2):182-187

    SONG Haibin, ZHANG Ling, JIANG Weiwei, et al. Geophysical researches on marine gas hydrates (Ⅲ): bottom simulating reflections[J]. Progress in Geophysics, 2003, 18(2):182-187.]

    [13]

    White R S. Gas hydrate layers trapping free gas in the Gulf of Oman[J]. Earth and Planetary Science Letters, 1979, 42(1):114-120. doi: 10.1016/0012-821X(79)90196-1

    [14]

    Stoll R D, Ewing J, Bryan G M. Anomalous wave velocities in sediments containing gas hydrates[J]. Journal of Geophysical Research, 1971, 76(8):2090-2094. doi: 10.1029/JB076i008p02090

    [15]

    Liu X W, He J, Sun Q L. Gas hydrate identification from △Vp/△Vs[C]//Beijing 2009 International Geophysical Conference & Exposition. Beijing: Society of Exploration Geophysicists, 2009: 150.

    [16]

    Tian D M, Liu X W. Identification of gas hydrate based on velocity cross plot analysis[J]. Marine Geophysical Research, 2021, 42(2):11. doi: 10.1007/s11001-021-09431-3

    [17]

    Chen J Z, You J C, Wei J T, et al. Interpreting XGBoost predictions for shear-wave velocity using SHAP: Insights into gas hydrate morphology and saturation[J]. Fuel, 2024, 364:131145. doi: 10.1016/j.fuel.2024.131145

    [18]

    Luo Y L, Zhang G L, Liang G W, et al. Limited-label multiscale deep-learning multihorizon tracking[J]. IEEE Transactions on Geoscience and Remote Sensing, 2024, 62:5918814.

    [19] 王迪, 张益明, 张繁昌, 等. 利用先验信息约束的深度学习方法定量预测致密砂岩“甜点”[J]. 石油地球物理勘探, 2023, 58(1):65-74

    WANG Di, ZHANG Yiming, ZHANG Fanchang, et al. Quantitative prediction of tight sandstone sweet spots based on deep learning method with prior information constraints[J]. Oil Geophysical Prospecting, 2023, 58(1):65-74.]

    [20]

    Bai Y, Tan M J. Dynamic committee machine with fuzzy-c-means clustering for total organic carbon content prediction from wireline logs[J]. Computers & Geosciences, 2021, 146:104626.

    [21] 陈钢花, 张寓侠, 王军, 等. 双向长短时记忆神经网络在滩坝砂储层岩性识别中的应用[J]. 测井技术, 2023, 47(3):319-325

    CHEN Ganghua, ZHANG Yuxia, WANG Jun, et al. Application of BiLSTM in lithology identification of beach-bar sand reservoir[J]. Well Logging Technology, 2023, 47(3):319-325.]

    [22]

    Chen Y J, Dunn K J, Liu X W, et al. New method for estimating gas hydrate saturation in the Shenhu area[J]. Geophysics, 2014, 79(5):IM11-IM22. doi: 10.1190/geo2013-0264.1

    [23]

    Zhu X Y, Liu T, Ma S, et al. Morphology identification of gas hydrate based on a machine learning method and its applications on saturation estimation[J]. Geophysical Journal International, 2023, 234(2):1307-1325. doi: 10.1093/gji/ggad133

    [24] 杨笑, 王志章, 周子勇, 等. 基于参数优化AdaBoost算法的酸性火山岩岩性分类[J]. 石油学报, 2019, 40(4):457-467

    YANG Xiao, WANG Zhizhang, ZHOU Ziyong, et al. Lithology classification of acidic volcanic rocks based on parameter-optimized AdaBoost algorithm[J]. Acta Petrolei Sinica, 2019, 40(4):457-467.]

    [25]

    Breiman L. Random forests[J]. Machine Learning, 2001, 45(1):5-23. doi: 10.1023/A:1010933404324

    [26] 任涛, 林梦楠, 陈宏峰, 等. 基于Bagging集成学习算法的地震事件性质识别分类[J]. 地球物理学报, 2019, 62(1):383-392

    REN Tao, LIN Mengnan, CHEN Hongfeng, et al. Seismic event classification based on Bagging ensemble learning algorithm[J]. Chinese Journal of Geophysics, 2019, 62(1):383-392.]

    [27]

    Bayar B, Stamm M C. Constrained convolutional neural networks: a new approach towards general purpose image manipulation detection[J]. IEEE Transactions on Information Forensics and Security, 2018, 13(11):2691-2706. doi: 10.1109/TIFS.2018.2825953

    [28]

    Nguyen B, Morell C, De Baets B. Large-scale distance metric learning for k-nearest neighbors regression[J]. Neurocomputing, 2016, 214:805-814. doi: 10.1016/j.neucom.2016.07.005

图(8)  /  表(3)
计量
  • 文章访问数:  33
  • HTML全文浏览量:  1
  • PDF下载量:  16
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-09-23
  • 修回日期:  2024-12-22
  • 网络出版日期:  2024-12-25
  • 刊出日期:  2024-12-27

目录

/

返回文章
返回