High-resolution seismic exploration technology with spark source and its application in identification of natural gas hydrates in marine areas
-
摘要:
精准识别天然气水合物对于保障国家能源安全、优化能源结构和环境保护等方面都具有重要的战略意义。地震勘探技术是水合物勘查的重要手段，但海域多道地震勘探多采用气枪震源，震源主频较低，导致地震波的分辨率有限，限制了其在探测浅表层或细微地质结构中的能力。随着水合物勘查精度和难度日益增长，采用气枪震源的地震勘探技术已无法满足精细刻画水合物矿体的需要。因此，提高地震勘探技术的分辨率对于水合物勘查而言至关重要。本文介绍了一种基于电火花震源的小道距高分辨率地震探测技术，通过与气枪震源的对比分析，深入探讨了电火花震源的特点及其在海域水合物勘查中的应用，研究结果表明，电火花震源地震勘探技术在扩散型水合物和浅表层渗漏型水合物识别中都具有良好的应用效果，尤其是对于浅表层渗漏型水合物，电火花震源地震可以更好地识别出海底和浅部地层的振幅异常、流体运移通道和海底微地貌，可有效提高浅表层渗漏型水合物识别的可靠性。
Abstract:Accurate identification of natural gas hydrates is of great strategic importance for ensuring national energy security, optimizing energy structures, and protecting the environment. Seismic exploration technology is a crucial means for hydrate exploration. However, most multi-channel seismic explorations in marine areas use airgun sources with relatively low dominant frequencies, resulting in limited resolution of seismic waves, which limits its ability to detect shallow surface layers or subtle geological structures. With the increasing in demand for accuracy and the difficulty of hydrate exploration, seismic exploration technology using airgun sources can no longer meet the needs of fine characterization of hydrate ore bodies. Therefore, improving the resolution of seismic exploration technology is crucial for hydrate exploration. This paper introduces a high-resolution seismic detection technology based on small-offset electric sparker sources. Through comparative analysis with airgun sources, the characteristics of electric sparker sources and their application in marine hydrate exploration are deeply discussed. Results show that high-resolution seismic exploration technology using electric sparker sources has good application effects in identifying both diffused hydrates and shallow leakage hydrates. Especially for shallow leakage hydrates, using electric sparker source seismology can better identify amplitude anomalies, fluid migration pathways, and subtle seabed morphologies in the seabed and shallow strata, which can improve effectively the reliability of identifying shallow gas hydrate.
-
Keywords:
- high-resolution seismology /
- sparker source /
- gas hydrate /
- seismic attribute
-
菲律宾海板块地处欧亚板块、太平洋板块和印度-澳大利亚板块的交汇处，四周几乎全部被俯冲带所包围，在西太平洋板块构造格局中占据特殊地位。受板块俯冲、海底扩张以及地幔柱活动等多重因素的影响，菲律宾海板块的形成演化过程表现为多阶段性，形成了现今复杂的构造格局，是开展西太平洋板块重建的关键区域菲律宾海蕴含丰富的前沿地球科学问题，是构造地质学、沉积学、古海洋学以及全球变化研究的热点地区。自1969年以来，DSDP-ODP-IODP在菲律宾海及邻区共完成了25个航次的钻探工作。钻探为“弧前增生楔结构和孕震机理”、“菲律宾海板块的形成演化”、“岛弧岩浆作用和大陆边缘陆壳增生机制”以及“西太平洋古气候和古海洋环境演变”等重大科学问题的研究提供了基础。然而，绝大部分钻探航次和站位位于俯冲带发育的主动大陆边缘，集中于日本南海海槽、伊豆-小笠原-马里亚纳沟-弧-盆系统等区域，少量钻孔分布在板块内部的海底高原和海脊地区，深海盆区仅有零星钻孔分布，制约了诸如花东海盆基底属性、俯冲起始过程与岛弧演化等关键科学问题的研究。针对前期大洋钻探工作的局限，预研未来该海域钻探选址与井位建议，不仅可揭示制约菲律宾海构造演化与沉积效应认知的关键地球科学问题，而且对未来我国自有大洋钻探船钻探井位储备具有重要意义。
1. 区域地质背景
菲律宾海板块呈南北向狭长的菱形走向，东部边界为伊豆-小笠原-马里亚纳海沟、雅浦海沟、帕劳海沟和阿玉海槽，西部边界为南海海槽、琉球海沟、马尼拉海沟和菲律宾海沟（图1）。以中央南北走向的九州-帕劳海脊将菲律宾海盆分为东、西两个大型的海盆：西菲律宾海盆和四国-帕里西维拉海盆。总体来说，基于磁异常与钻探数据，西菲律宾主海盆形成时间约为61～33/30 Ma[1]，四国海盆形成时间约为25～15 Ma[2]， 帕里西维拉海盆形成时间约为30～17 Ma[3]，马里亚纳海槽初始扩张时间约为6.5 Ma且尚处扩张时期[4]。
![]()
图 1 西太平洋海底地形和构造单元简图红线代表实测多道地震测线。Figure 1. Sketch map of topography and tectonic units in the Western PacificThe red lines represent new acquired multichannel seismic profiles.菲律宾海范围较大，南北跨越0°～35°N，东西横跨124°～147°E，面积达5.4×106 km2（图1）。海底地形十分复杂，发育有海沟、岛弧、海脊、海山、海盆、裂谷等地形，是一个海岭与海盆并列、具有洋壳基底的大型边缘海，全区水深3 000～6 000 m，平均水深约为4 500 m，总体表现为“西部深、东部浅”的特点（图1）。海盆整体位于菲律宾海板块之上，周边环绕岛弧与深海沟俯冲带，中部九州-帕劳海脊呈反转“S”形南北延伸，长约2 750 km，宽约50～90 km，海脊中段为一系列断续分布的链状海岭，海山与凹地相间分布，两侧呈明显不对称，东侧陡、西侧缓[5]。
西菲律宾海盆和帕里西维拉海盆地区，沉积基底面与海底地形起伏基本平行，其上主要是覆盖了一层第四纪沉积物，平均厚度约500 m，基本为等厚度分布，沉积速率非常低，沉积物来源具有空间差异性，基本以远洋自生沉积[6]、火山成因[4,7]和内陆风成沉积物质为主[8-10]，沉积物直接覆盖在海底玄武岩层之上。
2. 数据来源与研究方法
本文所用多道地震资料取自2020年9月执行的“深海科学钻探井位选址调查”项目菲律宾海调查航次，测量船舶为青岛海洋地质研究所“海洋地质九号”，在菲律宾海九州帕劳海脊和中央海盆断裂带东端共完成二维多道地震和同步多波束测量552.4 km（图1）。
多道地震数据采集系统采用SERCEL公司Sercel Seal 428系统，配以Sercel SSRD固体电缆，数据记录系统采用Seal 428 version1.1，震源系统采用SERCEL公司GII型气枪，空压机系统采用LMF 51S型空压机，拖缆控制系统采用ION公司PCS version 8.30，本次多道地震主要采集参数为接收360道，道间距12.5 m，炮间距37.5 m，震源容量为6 000 立方英寸，记录长度12 s，具体参数见表1。
表 1 地震数据采集系统参数Table 1. Parameters of seismic data acquisition system采集参数 采集参数 接收道数 360道 工作压力 2 000 psi 电缆长度/m 4 500 震源沉放深度/m 7 道间距/m 12.5 电缆沉放深度/m 10（试验后确定） 覆盖次数 60 最小偏移距/m 170 炮间距/m 37.5 低截频 3Hz@6dB/oct （模拟） 采样率/ms 2 高截频 400Hz@370dB/oct （线性） 记录长度/ms 12 030 记录格式 SEG-D 8058 震源容量/cu.in 6 000 记录介质 IBM359盒式磁带 3. 菲律宾海已有钻探科学目标实现与重要科学问题进展
迄今为止，DSDP-ODP-IODP在菲律宾海及邻区共完成了25个航次的钻探工作，在98个站位共完成钻孔218口，绝大部分钻探航次和站位集中于日本南海海槽、伊豆-小笠原-马里亚纳沟-弧-盆系统，少量钻孔分布在板块内部的海底高原和海脊地区（图2）。钻探主要以“弧前增生楔结构和孕震机理”、“菲律宾海板块的形成演化”、“岛弧岩浆作用和大陆边缘陆壳增生机制”以及“西太平洋古气候和古海洋环境演变”等为科学主题。
![]()
图 2 菲律宾海及邻区DSDP-ODP-IODP钻探站位分布图Figure 2. Distribution of DSDP-ODP-IODP drilling holes in the Philippine Sea and adjacent areas3.1 菲律宾海已有钻探航次目标和基本认识
3.1.1 深海钻探计划（DSDP）阶段
为了查明整体上菲律宾海盆的形成时代和盆地演化的历史，科学家们在西太平洋和菲律宾海实施了DSDP- Leg6、31航次的钻探工作，航次成果首次揭示了菲律宾海盆东部的年龄为古近纪或新近纪，证明菲律宾海板块是独立于太平洋板块之外的新生板块[11]；确定了菲律宾海板块形成和演化的年代学框架、古海洋环境变迁和古板块位置，拉开了菲律宾海形成演化过程重建的序幕[12]。为了详细探究不同具体位置的演化信息，科学家开展了DSDP-Leg 58航次，探究的是四国海盆的扩张作用和大东岭海盆起源的工作，四国海盆是由海底扩张作用形成的，而且可能存在单向扩张，大东盆岭区的岩石组成具有很大的非均一性，包括变质岩、火山岩和沉积岩[13]。随后科学家开展了DSDP-Leg 59、60的工作，确定基底年龄和盆地形成的机制，弧后盆地的扩张速率和火山弧的起源与性质[7]。确定了九州-帕劳海脊处的TTR型三联点，推测马里亚纳海沟的俯冲是一个平稳的过程，几乎不涉及来自太平洋板块的物质[14]。
3.1.2 大洋钻探计划（ODP）阶段
1990年前后在伊豆-小笠原-马里亚纳岛弧实施了ODP-Leg125、126航次的钻探工作，目的是为了探究岛弧地区的成因和演化问题。主要的科学成就是：揭示了弧前台地上更年轻的岩浆活动；弧前盆地中心以下火成岩基底为高镁系列安山岩以及低镁系列熔岩；从晚中新世到全新世和始新世—渐新世为火山活动的高峰[15-16]。后期为了探究日本南海海槽与伊豆-马里亚纳岛弧的关系，相继开展了ODP-Leg131[17]、190[18]、195[19]、196[20]航次的工作，对比发现两地区均有流体流动的过程和地震活动的特征。伊豆-马里亚纳岛弧有来自中央火山带的蛇纹岩浆，而日本南海海槽形成了大型逆冲推覆带，表现为粗粒物质从岛弧向海沟运移转变为横向运移，运移系统的转变使得沉积物从日本东部区域沿海沟轴线向下运移，为横断面的沉积物提供了来源。
3.1.3 综合大洋钻探计划（IODP）阶段
该阶段主要研究区域是日本南海海槽，分为多阶段的南海海槽孕震带试验（NanTroSEIZE），第一阶段实施了IODP314航次的工作，运用了最先进的随钻测井（LWD）技术，得到了天然气水合物赋存状态的重要信息；第一阶段后续的IODP315、316航次是在314航次4个站位的基础上增加了新的钻孔工作。主要研究的目标是斜俯冲中增生棱柱体与弧前盆地应变分量的构造研究，根据地层记录重建棱柱体和弧前盆地演化以及阐明巨型帚状断层系统浅部的特征、行为和演化信息[21]。IODP319、322作为第二个阶段，技术更加先进，主要科学目标包括在可能安装观测装置的深度或附近采集岩芯并进行井下试验，推断出四国海盆的粗粒火山沉积物可能是从位于东侧的火山链通过海底等途径运过来的，并成功地采集了构成基底的枕状玄武岩样品。2010年实施了第二阶段的IODP332、333的工作，在原有钻孔的基础上安装了井眼观测台，主要成果有完成了对四国海盆半远洋沉积物上层沉积物理性质转换带的取芯，对沉积物/玄武岩界面下方100 m范围内的玄武岩进行取芯等。南海海槽第三阶段钻探计划主要为IODP326、338和348航次，主要针对深部板块边界及断裂带进行隔水管大口径钻探，获取大型地震发育机理。
3.1.4 大洋发现计划（IODP）阶段
2014年国际大洋发现计划实施了350—352航次的科学钻探工作。IODP350航次是在伊豆弧后钻探的第一航次，以前IODP站位都位于伊豆-小笠原岛弧或弧前地区，导致对伊豆弧岩浆作用的认知并不完整。因此，IODP350航次的主要目的是揭示俯冲工厂“失踪的那一半”的演化历史[22]。第二航次，即IODO351航次将聚焦于九州-帕劳海脊西侧，本区被认为记录了IBM弧早期的基础、起源和演化过程[23]。第三航次，即IODP352航次将通过钻探IBM岛弧外侧的弧前地区，来研究俯冲起始过程[24]。这3个航次重点研究了地球化学不对称的岛弧地壳的起源、弧后岩浆活动、原始地壳与地幔性质、岛弧在古近纪的成分演化、俯冲开始后的地幔熔融演化过程、俯冲开始与初始岛弧地壳形成过程等科学问题。
综上所述，DSDP阶段的航次解决了菲律宾海盆的地质形成年代和演化历史的问题，对四国海盆和大东岭海盆的扩张起源问题有了初步的认识，确定了九州-帕劳海脊火山活动的时间和性质，提出了TTR型三联点：ODP阶段探究了伊豆-小笠原-马里亚纳岛弧的起源和演化问题，并与南海海槽做了对比，发现均有物质的迁移和地震带活动的特征；IODP阶段分为多阶段的南海海槽孕震带试验，探明了天然气水合物和增生棱柱体的信息，推测四国海盆沉积物的来源以及深部取芯；IODP阶段重点研究了地球化学不对称的岛弧地壳的起源、弧后岩浆活动以及俯冲开始与初始岛弧地壳形成过程等科学问题。但同时，仍有一些关键的科学问题制约着对菲律宾海及邻域构造沉积过程的全面认识。
3.2 菲律宾海前沿地球科学问题
3.2.1 花东海盆的年代和构造属性
花东海盆位于台湾岛东侧，是西太平洋地球系统中一个重要的盆地，位于菲律宾海板块与欧亚板块的交汇处及构造运动的转折点，构造位置极为重要[25]，是研究西太平洋板块从白垩纪到新生代构造演化与动力学机制的关键区域[26]（图3）。目前，对花东海盆的年龄和构造属性存在较大的争议，是西太平洋板块重建的关键。
![]()
图 3 琉球俯冲板片结构和花东海盆的构造属性[31]a. 琉球俯冲板片的重建结果，b. 俯冲之前的菲律宾海板块结构，c. 琉球俯冲板片三维结构。Figure 3. The structure of the Ryukyu subduction plate and the tectonic nature of the Huatung basin[31]a. The floated slab mesh connected at the Ryukyu Trench to the Philippine Sea Plate. b. Structure of the Philippine Sea Plate prior to subduction. c. 3D structure of the Ryukyu subducted slab.在新生代时期，菲律宾海板块侵位于欧亚板块与西太平洋板块之间，古地磁数据显示，在白垩纪时期花东海盆和菲律宾海板块位于赤道附近[27-28]，随后相对于欧亚板块沿着巨型剪切带一起向北移动，大位移地（近2 200 km）运行到现今的位置[27]，而花东海盆的洋盆性质最初是基于在海盆中识别出4个E-W向磁异常条带[11]，推断盆地的年龄有：41～36 Ma[1]、15～30 Ma[29]、52～43 Ma[30]等，据此认为它是菲律宾海板块的一部分，是菲律宾海板块楔入体的前锋[31]。后来拖网采得辉长岩样品[32]，单矿物40Ar/39Ar 测年为131～119 Ma，据此推断花东海盆年龄为早白垩世，是板块构造运动遗留给人类的、在菲律宾海盆与欧亚大陆之间唯一的中生代白垩纪洋盆[33]。因此，花东海盆形成时代和构造属性、与菲律宾板海板块的相关作用关系是重建中生代时期太平洋板块与欧亚板块之间的运动学和动力学关系一个非常关键的问题。
3.2.2 俯冲启动和岛弧形成过程
没有俯冲带就没有板块构造，俯冲带成因是理解板块构造的关键，然而长期以来俯冲带的起始俯冲机制一直是板块构造理论中最为薄弱的部分，西太平洋边缘海存在多个初始扩张中心和可能的初始俯冲带[34]，是研究初始俯冲的优良靶区。与被动大陆边缘的初始扩张相比，起始俯冲以及由此产生的洋内弧的研究才刚刚起步[35]。
菲律宾海四周几乎全被俯冲带包围，其东侧伊豆-小笠原-马里亚纳（IBM）俯冲带是由古老的西太平洋板块俯冲到菲律宾海板块与马里亚纳微板块之下而形成的，造就了典型的沟-弧-盆体系[36]。以往该地区的大洋钻探对IBM 50 Ma演化过程中近25 Ma的历史，通过火山灰和火山碎屑研究已取得显著进展。然而，对于IBM形成演化的前半段历史，尤其是岛弧形成的起始过程，却知之甚少。九州-帕劳海脊曾经与现代IBM岛弧连为一体，在约30 Ma的海底扩张形成四国海盆和帕里西维拉盆地之前一起组成了古IBM岛弧，大洋钻探发现计划连续三个航次在IBM展开[22-24]，但九州-帕劳海脊作为俯冲起始和岛弧形成演化的重要起源区[37]，开展的钻探工作明显不足，制约了对上述科学问题的认识。因此，以九州-帕劳海脊区为目标，对俯冲起始机制进行研究，依旧是板块构造学理论体系中的一个重点和难点。
同时，加瓜海脊介于西菲律宾海盆和花东海盆之间，构成了两个结构、性质、形成时代截然不同的海盆的边界。最新的多道地震和OBS探测结果表明[38]，加瓜海脊之下的地壳厚度相对于两侧海盆显著增厚，海脊地壳的增厚是非对称的，可能代表了中新世期间形成的一条夭折俯冲带，俯冲启动之后短时间内终止，俯冲带结构被“定格”在了加瓜海脊之下。因此，加瓜海脊可能也是研究俯冲启动过程的理想场所[34]。
3.2.3 菲律宾海板块的形成演化和板块重建
菲律宾海是东亚三大板块包围而成的最大边缘海。三大主要板块的运动符合基于海底扩张和热点参照系确定的全球板块运动模型[39]。然而菲律宾海板块自早始新世形成以来的运动轨迹，一直是不确定和有争议的。
多年来，不同学者提出了许多截然不同的菲律宾海板块重建方案[39-45]。现今菲律宾板块与太平洋板块均以WNW方向高速移动，需要其与太平洋板块之间存在较强的耦合才能实现。根据菲律宾海板块磁条带确定的海底年龄和扩张模式，为在一定误差范围内部分重建过去的板块几何结构提供了可能[1, 32, 46]，对于渐新世末到现在的菲律宾海扩张过程，已发表的菲律宾海年龄模型之间的差异相对较小。菲律宾海板块过去的运动状态与现今的运动状态截然不同，古地磁倾角数据表明菲律宾海板块晚中生代以来板块漂移中发生了110°的顺时针旋转[43]，但是却很难区分古地磁数据代表的是全板块旋转还是局部块体旋转。Lee和Hilde[40]的年龄模式被广泛用于菲律宾海板块的重建，认为西菲律宾海盆的扩张始于59～56 Ma。相反，基于最新放射性测年结果识别出的冲大东热点遗迹建立的新模型则认为西菲律宾海盆的扩张始于更晚的54～50 Ma[32, 47]。这种年龄与转向差异的不同直接影响了菲律宾海板块形成与演化重建的结果。
同时，西菲律宾海作为规模较大的不活动弧后盆地，其扩张作用终止的动力过程仍不明确。西菲律宾海盆残留的中央扩张中心呈NW-SE向，发育显著的中央裂谷，被近N-S向的古转换断层切割，扩张作用终止的动力学机制被认为是“失去动力”的主动停止过程，但消亡扩张中心在岩石圈结构、扩张末期熔体作用过程和物质组成等方面的差异，是进一步揭示海底扩张作用终止过程和机制的基础。
3.2.4 海洋风成沉积及其古环境记录
海洋风尘沉积是指陆源沉积物通过大气环流传输并在特定的海区沉降而形成的沉积物。每年大约有2 000亿t风尘进入大气层，其中25%落在海洋使得风尘沉积成为许多深海沉积物的重要组成部分[48]，而且亚洲是世界上最大的风尘源区之一，观测结果表明，亚洲风尘对海洋沉积物组成和全球气候变化具有重要的影响[49]。风尘通量和粒度变化是海洋风尘沉积研究中常用的古气候代用指标。ODP885/886孔的风尘粒度记录了11.8 Ma以来西风环流强度的演化历史，风尘粒径在4.6 Ma时增加，反映了晚中新世由于气候变冷所引起的西风加强[50]。风尘通量的显著增加反映了亚洲内陆干旱化程度的加剧[51]，而且干旱程度变化也会在石英通量上有所体现[52-53]，风尘通量的变化表现为明显的冰期-间冰期旋回特征，冰期的高风尘输入可能与亚洲内陆干旱化程度和东亚冬季风强度增强有关[54]。但是，西菲律宾海风尘沉积物的不同源区物质输入通量还无法定量划分，风尘沉积物在海流与西风环流的耦合搬运机制仍不清晰，同时亚洲风尘与大陆沉积尚缺少更长时间尺度和更高分辨率的比对记录，未来在菲律宾海盆钻取关键站位形成南北和东西断面对比，有望揭示新生代亚洲内陆干旱程度历史-西风强弱过程和东亚冬季风演化耦合过程及风尘沉积过程与机制。
4. 科学目标实现与站位选择
针对菲律宾海板块及周边地区的前沿地球科学问题，结合DSDP-ODP-IODP工作程度、实测多道地震数据以及菲律宾海板块调查研究工作进展，并综合考虑通过深海钻探解决科学问题的可行性和难易程度，聚焦3个科学目标，确定了未来可能的钻探井位，主要包括：①俯冲启动和岛弧早期演化；②三联点岩浆作用过程；③西菲律宾海盆年代和古地磁。
4.1 科学目标1—俯冲启动和岛弧早期演化
了解俯冲带的肇始机制，以及洋-洋俯冲产生的洋内火山岛弧的形成和演变过程，是板块构造理论最核心的问题之一。九州-帕劳海脊和加瓜海脊是不同类型俯冲启动的产物，是研究俯冲启动和岛弧早期形成演化的理想场所。在此开展钻探，不仅获取俯冲启动和岛弧演化早期的岩浆活动和沉积记录，甚至钻获夭折俯冲带之下的俯冲洋壳，能够为破解俯冲启动和岛弧形成演化过程提供可能。
九州-帕劳海脊曾经与现代IBM岛弧连为一体，在约30 Ma的海底扩张形成四国海盆和帕里西维拉盆地之前一起组成了古IBM岛弧，是古太平洋板块向菲律宾海板块俯冲的产物，视为古IBM的残留弧。九州-帕劳海脊及周缘原始洋壳基底之上与岛弧侵位同时期的火山碎屑沉积层，记录了俯冲板块对俯冲起始响应过程的证据，是了解板块如何起始俯冲的关键场所。
PSP-1站位的位置由2020年采集的MCS-2多道地震剖面控制，位于九州-帕劳海脊西侧，距离海脊顶部边缘约20 km，站位水深约5 460 m。PSP-1井以钻穿沉积层并揭露九州-帕劳海脊西侧的玄武岩洋底为目标，计划完钻深度6450 m，其中水深约5 460 m，沉积层和基底硬岩钻进990 m。预计将在海底之下790 m，即6250 m深度钻遇玄武岩基底。推测5 460～6 250 m深度的沉积层以深远海沉积物为主，岩性为深海软泥、泥质粉砂岩和粉砂质泥岩，夹少量陆源碎屑和火山碎屑。6 250～6 450 m为九州-帕劳海脊演化早期形成的中-酸性熔岩和火山碎屑岩，该层之下局部可能发育枕状玄武岩和辉长岩岩脉（图4）。
![]()
图 4 菲律宾海深海科学钻井预测图Figure 4. Diagram of sequence and lithology of the proposed drilling holesPSP-2站位的位置由2020年采集的MCS-2多道地震剖面控制，位于九州-帕劳海脊顶部，站位水深约3 380 m。PSP-4井以钻穿沉积层并揭露九州-帕劳海脊顶部的岛弧基底为目标，计划完钻深度4 200 m，其中水深约3 380 m，沉积层和基底硬岩钻进820 m。预计将在海底之下540 m，即3 920 m深度钻遇岛弧基底。推测
3 380～3 920 m深度的沉积层以深远海沉积物为主，岩性为深海软泥、泥质粉砂岩和粉砂质泥岩，夹少量陆源碎屑和火山碎屑。3 920～4 200 m为九州-帕劳海脊演化早期形成的中-酸性熔岩和火山碎屑岩（图4）。
4.2 科学目标2——板块演化与三联点岩浆作用过程
西菲律宾海盆中央断裂带东端与九州-帕劳海脊交汇处存在已经消亡的三联点。中央断裂带是西菲律宾海盆最后一期海底扩张的中心，而九州帕劳海脊代表太平洋板块向西菲律宾海盆之下俯冲形成的洋内弧。由于中央断裂带的形成时代早于俯冲作用且与俯冲带高角度相交，先存洋壳破裂带的存在势必会对俯冲作用引起的熔体作用过程产生显著影响，而俯冲作用也可能导致彼时已经停止活动或即将停止活动的扩张中心再次活化。二者之间的复杂相互作用过程，可能导致不同来源的深部熔体发生大规模的相互作用，并在海底岩石类型和同位素组成上留下显著的痕迹。
位于西菲律宾海盆中央裂谷北侧隆起的PSP-3站位由2020年采集的MCS-1多道地震剖面控制，距离中央裂谷轴部约30 km，站位水深约4 460 m。该井以钻穿沉积层并揭露西菲律宾海盆中央裂谷与九州-帕劳海脊三联点岩石组成为目标，计划完钻深度4 850 m，其中水深约4 460 m，沉积层和基底硬岩钻进390 m。预计将在海底之下220 m，即4 680 m深度钻遇基底(图4)。
位于西菲律宾海盆中央裂谷南侧的PSP-4站位由2020年采集的MCS-1多道地震剖面控制，位于中央裂谷东端南侧隆起之外的一个小型海底火山型凸起之上，距离中央裂谷轴部约80 km，站位水深约4 900 m。该井科学目标与PSP-3井相同，利于裂谷两侧对比，计划完钻深度5 180 m，其中水深约4 900 m，沉积层和基底硬岩钻进280 m。预计将在海底之下180 m，即5 080 m深度钻遇基底，该井层序地层格架与岩性与PSP-3井基本相同（图4）。
4.3 科学目标3—西菲律宾海盆年代和古地磁
古老的西菲律宾海盆年龄和运动轨迹解释存在的显著差异，制约了菲律宾海构造演化与板块重建研究。菲律宾海板块的构造发展史在很大程度上依赖于古地磁研究，然而，迄今为止用于菲律宾海板块古地磁测量的样品几乎全部来自易受局部块体旋转影响的板块边界附近，难以区分古地磁数据代表的是全板块旋转还是局部块体旋转。因此，通过钻探，在板块内部稳定区域尽可能多地采集样品，开展古地磁和放射性测年的研究，是准确开展菲律宾海板块重建的基础。
位于西菲律宾海盆的PSP-5站位由2020年采集的MCS-1和MCS-2多道地震剖面控制，位于MCS-1和MCS-2测线交点附近，距离九州-帕劳海脊约60 km，站位水深约5 550 m。PSP-5井以钻穿沉积层并揭露西菲律宾海盆中央裂谷东端北侧的海盆基底岩石组成为目标，计划完钻深度6 250 m，其中水深约5 550 m，沉积层和基底硬岩钻进700 m。预计将在海底之下550 m，即6 100 m深度钻遇基底。推测5 550～6 100 m深度的沉积层以深远海沉积物为主，岩性为深海软泥、泥质粉砂岩和粉砂质泥岩，夹少量陆源碎屑和火山碎屑。6 100 m之下以西菲律宾海盆海底扩张形成的枕状玄武岩为主，局部可能受九州-帕劳海脊演化早期形成的中-酸性熔岩和火山碎屑岩影响（图4）。
4.4 菲律宾海钻探井位汇总
本文根据目前的调查研究程度与研究进展，选择西菲律宾海盆中央断裂带与九州-帕劳海脊交汇处的三联点关键区域，选取钻探站位5个（表2、图5）。其中， PSP-1、PSP-2站位在九州-帕劳海脊探索可能的俯冲启动过程；PSP-3和PSP-4站位主要针对三联点岩浆作用；PSP-5站位目的在于获得额外的古地磁和年代学样品，为西菲律宾海盆的形成演化提供约束。
表 2 菲律宾海科学钻探建议井位位置Table 2. Proposed well locations of Philippine Sea Scientific Drilling序号 站位编号 站位水深/ m 地球物理测线控制 科学目标 1 PSP-1 5 460 MCS-2多道地震剖面 板块俯冲启动与岛弧演化 2 PSP-2 3 380 MCS-2多道地震剖面 板块俯冲启动与岛弧演化 3 PSP-3 4 460 MCS-1多道地震剖面 板块演化与三联点岩浆作用过程 4 PSP-4 4 900 MCS-1多道地震剖面 板块演化与三联点岩浆作用过程 5 PSP-5 5 550 MCS-1和MCS-2多道地震剖面 海盆年龄与古地磁 ![]()
图 5 菲律宾海深海科学钻探选址关键区域及井位分布图Figure 5. Distribution of the proposed wells in the Philippine SeaPSP-1孔设计终孔深度6 450 m，沉积层厚度约790 m，基底硬岩钻进200 m，采用全取芯钻探；PSP-2孔设计终孔深度4 200 m，沉积层厚度约540 m，基底硬岩钻进280 m，采用全取芯钻探；PSP-3孔设计终孔深度4 850 m，沉积层厚度约220 m，基底硬岩钻进170 m，采用全取芯钻探；PSP-4孔设计终孔深度5 180 m，沉积层厚度约180 m，基底硬岩钻进100 m，采用全取芯钻探；PSP-5孔设计终孔深度6 250 m，沉积层厚度约550 m，基底硬岩钻进150 m，采用全取芯钻探。
5. 结论
（1）系统分析了在菲律宾海及邻区已完成的DSDP-ODP-IODP的25个钻探航次进展与问题，发现绝大部分钻探航次和站位集中于日本南海海槽、伊豆-小笠原-马里亚纳沟-弧-盆系统，少量钻孔分布在板块内部的海底高原和海脊地区，深海盆区仅有零星钻孔分布。钻探目标以弧前增生楔结构和孕震机理、岛弧岩浆作用和大陆边缘陆壳增生机制等为科学主题。
（2）菲律宾海蕴含的花东海盆年龄与构造属性、板块俯冲起始与弧后盆地扩张停止机制、菲律宾海盆年代等科学问题等是制约菲律宾海板块构造演化、板块重建和古气候环境重建的关键。
（3）根据已有研究进展与实测资料，在国际大洋钻探稀少或空白区，在九州-帕劳海脊附近选定区域，聚焦不同科学目标实现甄选了未来深海钻探站位5个PSP1-5。
（4）根据地震资料和前期研究成果，对目标层位的深度和岩性特征进行了初步预测，钻探实施预计将推进以下科学目标的实现：① 进一步揭示俯冲启动机制和岛弧早期演化过程；② 进一步揭示三联点岩浆作用过程；③ 获取更多的西菲律宾海盆年代和古地磁数据，为菲律宾海及邻区的板块重建提供证据。
致谢：在论文评审过程中，审稿专家提出了专业细致的评审意见，使得文稿质量大幅提升；在此表示感谢，一并感谢2020年“海洋地质九号”西太航次调查全体船员及工作人员。
-
![]()
图 1 气枪震源与电火花震源叠加剖面及对应频谱[27]
Figure 1. Stacked profiles and corresponding frequency spectra of airgun and sparker sources
![]()
图 4 气枪震源与电火花震源单炮振幅能量对比
Figure 4. Comparison in amplitude energy between airgun and sparker sources
![]()
图 5 小道距高分辨率地震采集观测系统示意图
Figure 5. The observation system of small-scale high-resolution seismic acquisition
![]()
图 7 气枪震源和电火花震源数值模拟结果对比
Figure 7. Comparison of numerical simulation results between airgun and sparker sources
![]()
图 10 过浅表层水合物地震剖面（下）和沿海底下方振幅分布曲线（上）
Figure 10. Seismic profile (below) and amplitude distribution curve (upper) of shallow gas hydrate.
![]()
图 11 围绕气烟囱环状分布的浅表层渗漏型水合物地震剖面
Figure 11. Seismic profile of shallow gas hydrate distributed annularly around gas chimneys.
-
[1] Rincón-Martínez D, Ruge S M, Arias A S. Seismic analysis of the geological occurrence of gas hydrate in the Colombian Caribbean offshore[J]. Journal of South American Earth Sciences, 2022, 116:103800. doi: 10.1016/j.jsames.2022.103800
[2] Shipley T H, Houston M H, Buffler R T, et al. Seismic evidence for widespread possible gas hydrate horizons on continental slopes and rises[J]. American Association of Petroleum Geologists Bulletin, 1979, 63(12):2204-2213.
[3] Vanneste M, De Batist M, Golmshtok A, et al. Multi-frequency seismic study of gas hydrate-bearing sediments in Lake Baikal, Siberia[J]. Marine Geology, 2001, 172(1-2):1-21. doi: 10.1016/S0025-3227(00)00117-1
[4] Foschi M, Etiope G, Cartwright J A. Seismic evidence of extensive microbial gas migration and trapping in submarine gas hydrates (Rakhine Basin, Bay of Bengal)[J]. Marine and Petroleum Geology, 2023, 149:106100. doi: 10.1016/j.marpetgeo.2023.106100
[5] Yoo D G, Kang N K, Yi B Y, et al. Occurrence and seismic characteristics of gas hydrate in the Ulleung Basin, East Sea[J]. Marine and Petroleum Geology, 2013, 47:236-247. doi: 10.1016/j.marpetgeo.2013.07.001
[6] Pecher I A, Kukowski N, Huebscher C, et al. The link between bottom-simulating reflections and methane flux into the gas hydrate stability zone – new evidence from Lima Basin, Peru Margin[J]. Earth and Planetary Science Letters, 2001, 185(3-4):343-354. doi: 10.1016/S0012-821X(00)00376-9
[7] Zhang W, Liang J Q, Qiu H J, et al. Double bottom simulating reflectors and tentative interpretation with implications for the dynamic accumulation of gas hydrates in the northern slope of the Qiongdongnan Basin, South China Sea[J]. Journal of Asian Earth Sciences, 2022, 229:105151. doi: 10.1016/j.jseaes.2022.105151
[8] 宋海斌, 张岭, 江为为, 等. 海洋天然气水合物的地球物理研究(III): 似海底反射[J]. 地球物理学进展, 2003, 18(2):182-187 doi: 10.3969/j.issn.1004-2903.2003.02.002 SONG Haibin, ZHANG Ling, JIANG Weiwei, et al. Geophysical researches on marine gas Hydrates (III): bottom simulating reflections[J]. Progress in Geophysics, 2003, 18(2):182-187.] doi: 10.3969/j.issn.1004-2903.2003.02.002
[9] Yi B Y, Lee G H, Horozal S, et al. Qualitative assessment of gas hydrate and gas concentrations from the AVO characteristics of the BSR in the Ulleung Basin, East Sea (Japan Sea)[J]. Marine and Petroleum Geology, 2011, 28(10):1953-1966. doi: 10.1016/j.marpetgeo.2010.12.001
[10] Yuan T, Spence G D, Hyndman R D, et al. Seismic velocity studies of a gas hydrate bottom-simulating reflector on the northern Cascadia continental margin: amplitude modeling and full waveform inversion[J]. Journal of Geophysical Research: Solid Earth, 1999, 104(B1):1179-1191. doi: 10.1029/1998JB900020
[11] Crutchley G J, Maslen G, Pecher I A, et al. High-resolution seismic velocity analysis as a tool for exploring gas hydrate systems: an example from New Zealand's southern Hikurangi margin[J]. Interpretation, 2016, 4(1):SA1-SA12. doi: 10.1190/INT-2015-0042.1
[12] Turco F, Crutchley G J, Gorman A R, et al. Seismic velocity and reflectivity analysis of concentrated gas hydrate deposits on the southern Hikurangi Margin (New Zealand)[J]. Marine and Petroleum Geology, 2020, 120:104572. doi: 10.1016/j.marpetgeo.2020.104572
[13] Jeong T, Byun J, Choi H, et al. Estimation of gas hydrate saturation in the Ulleung basin using seismic attributes and a neural network[J]. Journal of Applied Geophysics, 2014, 106:37-49. doi: 10.1016/j.jappgeo.2014.04.006
[14] 沙志彬, 梁金强, 郑涛, 等. 地震属性在天然气水合物预测中的应用[J]. 海洋地质与第四纪地质, 2013, 33(5):185-192 SHA Zhibin, LIANG Jinqiang, ZHENG Tao, et al. The application of seismic attributes to the prediction of gas hydrates[J]. Marine Geology & Quaternary Geology, 2013, 33(5):185-192.]
[15] 王秀娟, 吴时国, 徐宁. 地震属性参数在识别天然气水合物和游离气分布模式中的应用[J]. 海洋与湖沼, 2006, 37(3):271-279 doi: 10.3321/j.issn:0029-814X.2006.03.012 WANG Xiujuan, WU Shiguo, XU Ning. Determining the distribution model of hydrate and free gas occurrence in sediment with seismic attribute parameters[J]. Oceanologia et Limnologia Sinica, 2006, 37(3):271-279.] doi: 10.3321/j.issn:0029-814X.2006.03.012
[16] 邢磊. 海洋小多道地震高精度探测关键技术研究[D]. 中国海洋大学硕士学位论文, 2012 XING Lei. Study of the key technologies of high-precision marine multichannel seismic survey[D]. Master Dissertation of Ocean University of China, 2012.]
[17] 戚宾, 王祥春, 赵庆献. 海洋电火花震源地震勘探研究进展[J]. 物探与化探, 2020, 44(1):107-111 QI Bin, WANG Xiangchun, ZHAO Qingxian. Research on the progress of marine sparker seismic exploration[J]. Geophysical and Geochemical Exploration, 2020, 44(1):107-111.]
[18] Buogo S, Cannelli G B. Implosion of an underwater spark-generated bubble and acoustic energy evaluation using the Rayleigh model[J]. The Journal of the Acoustical Society of America, 2002, 111(6):2594-2600. doi: 10.1121/1.1476919
[19] 卢新培, 潘垣, 张寒虹. 水中脉冲放电的电特性与声辐射特性研究[J]. 物理学报, 2002, 51(7):1549-1553 doi: 10.3321/j.issn:1000-3290.2002.07.024 LU Xinpei, PAN Yuan, ZHANG Hanhong. The electrical and acoustical characteristics of pulsed discharge in water[J]. Acta Physica Sinica, 2002, 51(7):1549-1553.] doi: 10.3321/j.issn:1000-3290.2002.07.024
[20] 秦曾衍, 左公宁, 王永荣, 等. 高压强脉冲放电及其应用[M]. 北京: 北京工业大学出版社, 2000 QIN Zengyan, ZUO Gongning, WANG Yongrong, et al. The Application of High Voltage and Super-Density Pulse Discharge[M]. Beijing: Beijing Industry University Press, 2000.]
[21] 裴彦良, 王揆洋, 李官保, 等. 海洋工程地震勘探震源及其应用研究[J]. 石油仪器, 2007, 21(2):20-23 PEI Yanliang, WANG Kuiyang, LI Guanbao, et al. Application study of marine engineering seismic sources[J]. Petroleum Instruments, 2007, 21(2):20-23.]
[22] 刘怀山, 王文秋, 尹燕欣. 海洋地球物理探测在海岸工程建设中的应用: 以曹妃甸港为例[J]. 海岸工程, 2022, 41(4):313-327 doi: 10.12362/j.issn.1002-3682.20220727001 LIU Huaishan, WANG Wenqiu, YIN Yanxin. Offshore geophysical prospecting on coastal engineering: a case study in Caofeidian[J]. Coastal Engineering, 2022, 41(4):313-327.] doi: 10.12362/j.issn.1002-3682.20220727001
[23] Anitha G, Ramana M V, Ramprasad T, et al. Shallow geological environment of Krishna–Godavari offshore, eastern continental margin of India as inferred from the interpretation of high resolution sparker data[J]. Journal of Earth System Science, 2014, 123(2):329-342. doi: 10.1007/s12040-013-0399-3
[24] Mangipudi V R, Goli A, Desa M, et al. Synthesis of deep multichannel seismic and high resolution sparker data: implications for the geological environment of the Krishna–Godavari offshore, Eastern Continental Margin of India[J]. Marine and Petroleum Geology, 2014, 58:339-355. doi: 10.1016/j.marpetgeo.2014.08.006
[25] 易虎, 詹文欢, 闵伟, 等. 小多道地震震源效果在海域活动断裂探测中的对比研究[J]. 地震地质, 2022, 44(2):333-348 doi: 10.3969/j.issn.0253-4967.2022.02.004 YI Hu, ZHAN Wenhuan, MIN Wei, et al. A comparative study of source effect based on mini-multichannel seismic profile in marine active fault detection[J]. Seismology and Geology, 2022, 44(2):333-348.] doi: 10.3969/j.issn.0253-4967.2022.02.004
[26] 骆迪, 蔡峰, 闫桂京, 等. 浅表层天然气水合物高分辨率地震勘探方法与应用[J]. 海洋地质前沿, 2020, 36(9):101-108 LUO Di, CAI Feng, YAN Guijing, et al. High resolution seismic method for shallow gas hydrates exploration[J]. Marine Geology Frontiers, 2020, 36(9):101-108.]
[27] Luo D, Cai F, Wu Z Q. Numerical simulation for accuracy of velocity analysis in small-scale high-resolution marine multichannel seismic technology[J]. Journal of Ocean University of China, 2017, 16(3):370-382. doi: 10.1007/s11802-017-3145-7
[28] 骆迪, 蔡峰, 吴志强, 等. 海洋短排列高分辨率多道地震高精度成像关键技术[J]. 地球物理学报, 2019, 62(2):730-742 doi: 10.6038/cjg2019M0178 LUO Di, CAI Feng, WU Zhiqiang, et al. The key technologies of marine small scale high resolution multichannel seismic high-precision imaging[J]. Chinese Journal of Geophysics, 2019, 62(2):730-742.] doi: 10.6038/cjg2019M0178
[29] 蔡峰, 吴能友, 闫桂京, 等. 海洋浅表层天然气水合物成藏特征[J]. 海洋地质前沿, 2020, 36(9):73-78 CAI Feng, WU Nengyou, YAN Guijing, et al. Characteristics of shallow gas hydrates accumulation in the sea[J]. Marine Geology Frontiers, 2020, 36(9):73-78.]
[30] Barnes P M, Lamarche G, Bialas J, et al. Tectonic and geological framework for gas hydrates and cold seeps on the Hikurangi subduction margin, New Zealand[J]. Marine Geology, 2010, 272(1-4):26-48. doi: 10.1016/j.margeo.2009.03.012
[31] Ye J L, Wei J G, Liang J Q, et al. Complex gas hydrate system in a gas chimney, South China Sea[J]. Marine and Petroleum Geology, 2019, 104:29-39. doi: 10.1016/j.marpetgeo.2019.03.023
[32] Hovland M, Svensen H. Submarine pingoes: indicators of shallow gas hydrates in a pockmark at Nyegga, Norwegian Sea[J]. Marine Geology, 2006, 228(1-4):15-23. doi: 10.1016/j.margeo.2005.12.005
[33] Ben-Avraham Z, Reshef M, Smith G. Seismic signature of gas hydrate and mud volcanoes of the South African continental margin[C]//Proceedings of the NATO Advanced Research Workshop on Mud Volcanism, Geodynamics and Seismicity. Baku, Azerbaijan: Springer, 2005: 17-27.
[34] Bohrmann G, Ivanov M, Foucher J P, et al. Mud volcanoes and gas hydrates in the Black Sea: new data from Dvurechenskii and Odessa mud volcanoes[J]. Geo-Marine Letters, 2003, 23(3-4):239-249. doi: 10.1007/s00367-003-0157-7
-
期刊类型引用(1)
1. 徐欢，邓帅涛，唐玲，潘婷婷，马志军，高子靖，薛程. 基于卫星追踪技术分析崇明东滩越冬小天鹅迁徙路线及特征. 生态学杂志. 2025(01): 344-352 . 
百度学术
其他类型引用(0)
下载:
计量
- 文章访问数: 27
- HTML全文浏览量: 3
- PDF下载量: 14
- 被引次数: 1
