地球物理学报  2018, Vol. 61 Issue (7): 2905-2914   PDF    
南海北部琼东南海域活动冷泉特征及形成模式
杨力, 刘斌, 徐梦婕, 刘胜旋, 关永贤, 顾元     
国土资源部海底矿产资源重点实验室, 广州海洋地质调查局, 广州 510760
摘要:近年来,活动冷泉的研究越来越受到关注.本文利用多波束数据、多道地震数据以及底质取样结果研究琼东南海域活动冷泉系统,分析活动冷泉的羽状流特征、海底地貌与底质特征以及流体活动构造特征.多波束水体数据上,观测到多个延伸高度超过750 m的气泡羽状流,海底流体活动非常强烈;多道地震上识别出麻坑、流体运移通道、气烟囱等流体渗漏相关的构造,与其他海域观测到的反射特征不同,羽状流的下方流体运移通道呈强振幅"串珠"反射;重力活塞取样在两个站位上获得浅表层块状天然气水合物.其中一个站位位于活动冷泉附近,天然气水合物赋存于海底以下8 m左右.基于以上三方面的数据,笔者提出了一个用于描述活动冷泉系统的形成模式,游离气通过气烟囱向上运移到达浅层,一部分在天然气水合物稳定带内形成天然气水合物,另一部分穿透天然气水合物稳定带到达海底,形成活动冷泉的羽状流.
关键词: 琼东南海域      活动冷泉系统      羽状流      近海底天然气水合物     
Characteristics of active cold seepages in Qiongdongnan Sea Area of the northern South China Sea
YANG Li, LIU Bin, XU MengJie, LIU ShengXuan, GUAN YongXian, GU Yuan     
Key Laboratory of Marine Mineral Resources, Ministry of Land and Resources, Guangzhou Marine Geological Survey, Guangzhou 510760, China
Abstract: Active cold seepage is a major concern, since the methane released from cold seepage may have an impact on the climate change and is greatly associated with the seafloor ecosystem. Using multi-beam data, sediment cores and multi-channel seismic data, we investigated active cold seep systems in the Qiongdongnan sea area of the northern South China Sea including gas bubble plumes, seafloor features, near-seafloor gas hydrate and the structure of the fluid flow. Multi-beam water column data identify several over 750 m high gas bubble plumes, suggesting that seabed fluid is highly active; Multi-channel seismic data show pockmarks, gas chimneys and fluid migration channel. Unlike the features of fluid migration pathway observed in other oceans, the fluid migration channel manifests itself as a string of "beads" on multi-channel seismic profile; Two massive gas hydrate samples were recovered through piston core, one of which occurs near active cold seepage and lies 8 m below seafloor. Based on these data and analysis, we proposed a model to describe the formation of the cold seep system. Free gases rise to the shallow sediments through gas chimney. Some of them are trapped in gas hydrate stability zone (GHSZ) as gas hydrate, while others enter the ocean water through GHSZ and trigger gas bubble plumes. Our results are significant to understand the migration of seabed fluid and the deposit of gas hydrate near seafloor.
Key words: Qiongdongnan Sea Area    Active cold seepage    Gas plume    Near-seafloor gas hydrate    
0 引言

冷泉是指来自海底沉积界面以下,成分以水、碳氢化合物(甲烷为主)、硫化氢和细粒沉积为主的低温流体(通常只有几摄氏度),以渗漏、喷涌或扩散的方式向海底面运移,并在甲烷渗漏区滋养了大量化能自养生物群落的一种海底环境(Judd and Hovland, 2007).由于其温度相对于海底热液而言属于与周围海水温度相似的低温流体,因此被称为冷泉.

天然气水合物储藏划分为渗漏系统和扩散系统,冷泉正是渗漏系统的典型特征,且较扩散系统而言,渗漏系统里天然气水合物成藏更迅速,成藏储量更大,成藏物化条件更好、更具备开采价值(关进安等, 2009).在墨西哥湾(Kastner et al., 2005Solomon et al., 2005)和Cascadia南北水合物脊(Liu and Flemings, 2006)的研究资料表明,这些地区的天然气渗漏仍在持续,说明在天然气水合物稳定带内(GHSZ)里可能仍有大量的天然气水合物在不断生成.深水冷泉区是天然气水合物产出的理想场所,冷泉系统中发育的天然气水合物具有埋藏浅,品质高的特点.此外,它与全球气候变化及极端环境下的生态系统之间密切相关.因此,冷泉的研究具有非常重要的价值和意义.

冷泉广泛发育于活动和被动大陆边缘海底,从热带到两极极区,从浅海陆架到深海海沟均有分布,海底深度从10~3000 m不等(Logan et al., 2010),现已在全球大陆边缘海底发现活动冷泉上千个,如俄勒冈州沿岸、Cascadia水合物脊、日本海、尼罗河深海扇、新西兰大陆边缘、哥斯达黎加边缘海、地中海、黑海、阿留申群岛和中国南海等地区(Judd and Hovland, 2007),其中Cascadia水合物脊(Torres et al., 2004)、俄勒冈州沿岸(Liu and Flemings, 2006)和墨西哥湾附近地区(Duperron et al., 2007Brown et al., 2013)的冷泉活动研究程度最高.探测冷泉的主要手段包括:ROV探测、多波束探测、侧扫声呐探测、浅地层剖面探测、多道地震勘探等.利用这些手段获取冷泉的直接或间接证据,如:羽状流、麻坑、泥火山、丘状体、天然气水合物、碳酸盐岩、生物群落等(Judd and Hovland, 2007).在多道地震剖面上表现为海底火焰、声空白、毯式反射(Hovland and Judd, 1988Taylor, 1992Garcia-Gil et al., 2002)、帘式反射(Taylor, 1992Garcia-Gil et al., 2002Karisiddaiah et al., 1993)、柱式反射(Hovland and Judd, 1988)、浊反射(Garcia et al., 2002Fannin, 1980)、气体通道、泥底辟等.

本文研究区位于南海北部陆坡琼东南海域(图 1a),多道地震上广泛分布的BSR表明该海域存在天然气水合物(王秀娟等, 2010).前人的研究也表明该区域流体活动非常强烈,发育有大量的麻坑、泥火山等冷泉相关的构造(拜阳等, 2014陈江欣等, 2015, 2016).2015年的ROV调查在该海域发现了大面积的冷泉系统,并根据探测器的名字命名为“海马冷泉”.在2016年采集的多波束水体影像数据上,发现四个延伸高度超过750 m的气泡羽状流.本文利用多波束数据,底质取样结果以及多道地震对羽状流区域的冷泉系统进行研究,分析该冷泉系统的特征及其形成模式.

图 1 研究区示意图 (a)研究区位置(红色方框); (b)研究区地震测线; (c)羽状流区域. Fig. 1 Study area schematic (a) Study area (red box); (b) Seismic line in study area; (c) Gas plumes in study area.
1 区域地质背景

琼东南盆地属于新生代被动大陆边缘型盆地(王秀娟等,2010),是南海北部陆缘张裂而形成的一个叠合盆地,经历古近纪的张裂和新近纪以来热沉降过程,沉积了巨厚的地层,最大沉积厚度超过8 km,最大埋深超多9.4 km.盆地新生代沉积最大厚度达12 km,最大生烃凹陷面积超过9000 km2(何家雄等,2015).第四系和第三系上新统海相泥岩具有生物气形成的条件和良好的远景, 第三系煤系地层普遍处于异常高压和高地温条件, 具有形成热成因天然气的地质条件(黄保家,2002).因此, 盆地具有丰富的热成因和生物成因天然气资源(陈多福等,2004).此外,研究区裂隙和断层通道发育,区域存在较多的多边形断层, 即便在缺少连通烃源岩的断裂系统的情况下,流体依然可以运移(吴时国等,2009).研究区(如图 1中红色方框所示)位于琼东南盆地南部,水深约1250~1500 m,海底地形平坦,其地温梯度约为35~42.5 ℃/km(陈多福等,2004).总之,研究区有充足的气源条件,多样的运移通道,适宜的温压条件,适合天然活动冷泉的形成.

2 数据和方法

所用数据包括多波束、底质取样和多道地震,其中多波束数据由海洋六号于2016年采集完成,底质取样由海洋四号于2015年获得,多道地震数据由奋斗四号分别于2004年和2007年采集完成.

多波束利用高频声波,以开角探测的方式,对海底进行调查,主要包括海底地形的刻画,海底底质回波强度的探测和水体的成像等.本次调查使用挪威Kongsberg公司的EM122多波束测深系统,它主要由换能器阵列、发射接收主控单元、Windows XP实时监控工作站组成,其中换能器阵列固定安装于船前部的龙骨处,换能器吃水深度5.5 m.后处理流程包括数据格式转换,加载零潮汐数据,合并数据,建立BASE曲面,导航定位编辑,测线自动滤波,曲面自动滤波,子区编辑,人工数据恢复,曲面插值,数据输出和精度评价.水体数据处理由WCIP软件通过背景噪声消除,水体异常提取,水体影像横向叠加和水体影像纵向叠加完成.

底质取样包括重力柱状取样,重力活塞取样,箱式取样,拖网取样等取样方式,本次底质取样采用重力活塞取样方式,重力活塞约1000 kg,管长9 m,管径73 cm,内装透明管,用于获取深海底表层柱状样品,从而进行现场沉积物描述和化学分析测试.

多道地震调查利用气枪震源和拖缆,通过地震波的激发和接收,经过地震资料处理,对海底沉积层进行成像.多道地震采集参数见表 1.地震资料处理包括:定义观测系统,投影滤波去大值噪声,高通滤波去低频干扰,速度滤波压制线性干扰,球面扩散补偿,道分选,速度分析,正常时差动校正,拉冬变换去常规多次波,动校正畸变远道切除,近道切除压制多次波,叠加,带通滤波,TX域视速度滤波消除大倾角强绕射,FX域偏移.

表 1 多道地震采集参数 Table 1 Multi-channel seismic acquisition parameters
3 结果 3.1 羽状流多波束水体特征

冷泉区多波束海底地形如图 1c所示.研究区水深在1250~1500 m之间,坡度小于2 °.可以看到,该活动冷泉区海底地形平坦.在羽状流的位置上,也并无明显的异常地貌.

研究区存在3处羽状流,分别为A冷泉、B冷泉和C冷泉(图 1c的A、B、C).如图 2所示,A冷泉直径26 m,底深1380 m,顶深645 m,高差735 m;B冷泉直径28 m,底深1380 m,顶深630 m,高差750 m;C冷泉直径50 m,底深1420 m,顶深650 m,高差770 m.3处冷泉的多波束水体影像图均显示,在水深约800 m处,羽状流的形态扭曲最大,存在一个水层界面,界面以上的羽状流比界面以下的羽状流摆幅程度大.

图 2 冷泉A、B、C多波束水体影像 Fig. 2 Multiple beam water images of A, B and C cold seepage

不同时刻采集数据的成像结果如图 3所示.羽状流A、B距离多波束主测线0063距离最近,而多波束测线0060和0066分布在测线0063两侧.三条测线测量时间间隔约3 h.在三次测量时刻,两处冷泉尺度变化不大,高差在750 m左右,直径为26~50 m,但不同时刻3处冷泉的形态发生变化.

图 3 不同时刻采集数据的A、B羽状流成像结果 (a) 0060测线的成像结果; (b) 0063测线的成像结果; (c) 0066测线的成像结果.测线采集的时间间隔约为3个小时,三次测量时刻羽状流形态发生了较大的变化 Fig. 3 The results of A and B plume collected at different times (a) Result of multi-beam data from line 0060; (b) Result of multi-beam data from line 0063; (c) Result of multi-beam data from line 0066.
3.2 海底地质特征

研究区内,在GH1和GH2两个站位上获得块状天然气水合物,其中GH1站位位于羽状流C附近(图 1),GH1站位上天然气水合物位于7.95~8.20 mbsf,约2 cm3,底部为黏土层,GH2站位天然气水合物位于4.95~5.11 mbsf,其尺寸为5 cm×5 cm×2 cm,底部为黏土层.现场化学分析表明,2个站位天然气水合物样品以甲烷绝对优势为特征,甲烷/乙烷值大于10000,初步判断为Ⅰ型天然气水合物,气体来源主要是生物气.

3.3 流体活动构造

多道地震的成像结果如图 4-7所示.从多道地震剖面上,可以看到此区域存在气烟囱、BSR、振幅空白带和杂乱反射、麻坑、裂隙与断层等与天然气水合物相关的地震反射特征.

图 4 测线Line 1(经过羽状流C)的地震成像剖面 Fig. 4 The Profile of Line 1(across C gas plume)
图 5 测线Line 2(经过羽状流B)的地震成像剖面 Fig. 5 The Profile of Line 2(across B gas plume)
图 6 羽状流C附近部分测线Line 1的地震剖面(图 4虚线方框部分) 测线上分布3个麻坑,分别为pm1、pm2和pm3,3个麻坑海底以下有明显的气体通道,其中pm3的气体通道有明显的串珠反射特征,海底以下约300 ms存在BSR Fig. 6 The section profile of Line 1 across C gas plume(dotted box in Fig. 4)
图 7 羽状流B附近测线Line 2的地震剖面(图 5虚线方框部分) 羽状流B位置上存在一个明显的麻坑,麻坑以下存在连续的“串珠”反射,BSR存在于海底以下约200 ms,BSR以下存在高频衰减和速度下拉现象. Fig. 7 The section profile of Line 2 across B gas plume(dotted box in Fig. 5)

测线Line1(图 1b黑线)位于羽状流正上方,其成像结果如图 4所示.从成像结果上可以看到,海底地形平坦,未见明显的异常地貌特征.海底以下存在BSR,其相位与海底反射相位相反,呈强反射同相轴,横向连续性较差,指示存在天然气水合物.在羽状流C位置的正下方存在大规模的气烟囱.

测线Line2(图 1b黑线)位于羽状流B的正上方,其成像结果如图 5所示.从成像结果上可以看到,海底地形起伏不大,存在5个小规模麻坑和一个丘状体,但在这些异常的位置上未发现羽状流,羽状流B出现在海底地形平坦的位置上.测线中部海底以下存在BSR现象,其下方存在大规模的气烟囱,这种由气体引起的模糊反射在BSR以上的浅层同样存在,这表明该地区浅层气充分.

为更好地分析羽状流与海底地形以及地层构造之间的关系,我们对与羽状流B、C附近的地震剖面进行放大显示.羽状流C附近的地震剖面(图 4虚线方框)如图 6所示.海底有三个麻坑,分别为pm1、pm2和pm3.pm1和pm2的直径约300 m,其高差约3.5 m,pm1呈凹形口,pm2呈凸形口.pm3的直径约50 m,其高差不明显,pm3呈喷口特征.A冷泉和B冷泉在Line1与Line2的交线处附近,其近海底存在串珠反射.此外,海底浅层多处存在串珠状反射特征,并且微裂隙运移通道发育.

羽状流B附近的地震剖面(图 5虚线方框)如图 7所示.地震剖面显示在羽状流B位置上存在一个明显的麻坑(pm4).麻坑pm4直径约70 m,高差约3 m,呈凹形口,其下部存在具有串珠状反射特征的垂向连续气体运移通道.海底浅层存在多处串珠状反射以及通道,BSR以上的浅层存在相互连通的垂向微裂隙,BSR以下存在由于气体引起的高频衰减和同相轴下拉的反射特征.

4 讨论

在活动冷泉区海底往往存在羽状流和与天然气水合物相关的地貌特征,如:泥火山、水合物脊、丘状体和麻坑等.研究区海底存在羽状流,但与天然气水合物相关的地貌特征并不明显,因此,是否能探测到天然的羽状流显得尤为关键.相比其他常规调查手段,多波束探测技术一方面对羽状流的探测灵敏,另一方面能实现全覆盖调查,在活动冷泉调查中有着独特的优势.本次研究正是始于多波束探测到3处天然羽状流.

4.1 羽状流随时间变化

冷泉的活动性一般会随时间发生变化,主要体现在气体通量和羽状流的形态上.在斯瓦尔巴德岛近两年的连续观测表明冷泉的气体通量存在不同周期的变化,既包括季节性的变化,也包括短周期的日变化(Berndt et al., 2014). Hsu等(2013)在台湾西南岸气体渗漏点进行了连续24 h的声学和海底地震仪(OBS)监测,结果表明羽状流的形态在短时间内(3 h)会发生较大的变化,并受控于潮汐的变化.在本研究区,通过对不同时刻(间隔约3 h)采集的水体数据进行成像,成像结果显示羽状流的规模比较稳定,尽管形态有些变化(如图 3所示),表明观测时间内气体通量稳定.但由于观测时间较短,连续采集的次数较少,本次研究结果并不能排除气体通量长周期变化的可能性.

4.2 海底地貌特征

多波束地形图显示该冷泉区域海底地形平坦,并没有明显的流体渗漏特征,这与其他海域观测到的现象不一致.在其他海域,冷泉系统一般在海底有明显的地貌特征,比如泥火山,麻坑或丘状体等.但在该研究区,多波束数据基本没有揭示这些特征.而在地震数据上,却可以看到多处麻坑.造成这一差异的原因可能有两个:1.多波束资料和地震资料的采集时间相差近10年,而研究区物源丰富,水体流动活跃,古麻坑可能被沉积物覆盖.2.船载多波束的分辨率有限(马小川等, 2017),不能显示海底地形的小尺度变化.为获得高精度的海底地形,下一步应该开展基于ROV或者AUV的多波束测量.

4.3 流体通道内“串珠”反射特征的形成机制

在其他海域观测到的冷泉系统的流体通道多为断层与裂隙,断层在地震剖面上表现为通道两侧地层为强振幅错开特征,裂隙在地震剖面上表现为弱振幅,其两侧地层表现为强振幅特征(Leven and Middleton, 2007Wenau et al., 2015a陈江欣等, 2017).但是在我们的研究区域发现,流体通道在地震剖面呈现一连串的纵向强振幅(Schroot et al., 2005Wenau et al., 2015bVielst ädte et al., 2015),即“串珠”反射特征(图 8a).这与碳酸盐岩缝洞体储层形成的反射结果类似(图 8b).前人对碳酸盐缝洞体储层形成反射研究比较多,物理模拟和数值模拟表明,一定尺度空洞内存在的低速填充物是形成串珠反射的原因(Wood et al., 2008Bünz et al., 2012马灵伟等, 2015).由此我们推断,在该研究区冷泉系统的流体通道内自上而下存在一系列的低速体,可能对应游离气聚集(图 8c).这种离散的游离气聚集可能是受岩性不同控制的.

图 8 “串珠”反射形成机制示意图 (a)冷泉海底以下“串珠”反射特征; (b)碳酸盐岩缝洞体储层的反射结果; (c)可能的游离气聚集模式. Fig. 8 Schematic of "bead" reflectionform mechanism (a) The "beaded" reflection feature below the cold seepage; (b) The reflection of the reservoir of the carbonate rock seam; (c) a possible free gas accumulation pattern.
4.4 游离气穿过天然气水合物稳定带进入海水层

在水深超过500 m的海域,海底浅部沉积层满足天然气水合物稳定条件.在到达海底并进入海水层以前甲烷气体必须穿过天然气水合物稳定(Gorman et al., 2002).这就产生了一个矛盾,游离气在通过稳定带时会以水合物的形式赋存下来而堵住通道,这样后续的气体难以继续运移形成稳定的羽状流(Tréhu et al., 2004).这就需要某种机制阻止水合物的形成,可能的机制主要有四种:(1)天然气水合物在通道壁上覆盖,形成一个管状通道,将孔隙水与游离气分离(Leifer and MacDonald, 2003); (2)热流体局部对流; (3)孔隙水氯含量和盐度的增大; (4)水供应不足(Milkov, 2004; Ruppel et al., 2005).基于地震成像,推测冷泉B和C处的气体通过高渗透通道到达海底,天然气水合物在通道壁上覆盖,形成一个管状通道,将孔隙水与游离气分离.当然,本次研究结果也不能排除其他三种可能.

4.5 活动冷泉系统的形成模式

综合上述讨论,我们提出研究区活动冷泉形成模式(图 9),天然气通过气烟囱由深部向上运移,天然气以扩散和聚集的方式进入GHSZ.在羽状流下方,游离气通过管状通道向上运移和聚集,在海底处形成麻坑,进入海水层形成羽状流.在GHSZ的其他区域,有明显的纵向气体运移通道,推测为BSR以下游离气向上扩散或GHSZ内天然气水合物分解的结果,这种结果将促使活动冷泉更加活跃并在近海底形成天然气水合物.

图 9 琼东南海域研究区活动冷泉形成模式 Fig. 9 The pattern of cold seepages in Qiongdongnan Sea Area of the northern South China Sea
5 结论

本文基于多波束数据、底质取样结果以及多道地震数据对南海北部陆坡琼东南海域的活动冷泉系统进行研究.该活动冷泉系统存在三处形态随时间变化的羽状流,在多波束水体数据上表现为高度超过735 m、直径大于26 m的火焰状声学异常,其所在海底地形平坦,海底浅表层赋存天然气水合物.在多道地震剖面上,与天然气水合物相关的地质反射特征明显,包括气烟囱、BSR、声空白、高频衰减、同相轴下拉等.在羽状流的正下方存在由游离气产生的强振幅反射、由低速体引起的强振幅“串珠”状反射.

与其他地区活动冷泉相比,南海北部琼东南海域的活动冷泉的特征独特,主要表现在三个方面:(1)在观测时间内冷泉气体通量大且稳定; (2)海底地形特征不明显,麻坑纵向尺度为米级; (3)冷泉的通道在地震剖面上表现为明显的“串珠”状反射,解释为冷泉的流体通道内自上而下存在一系列的低速体,对应游离气聚集.

结合以上特征,我们提出了该活动冷泉系统的形成模式,首先天然气通过气烟囱向上运移到达浅层,受温压条件、地质构造和沉积等因素控制,一部分气体在天然气水合物稳定带内形成天然气水合物,另一部分游离气通过高渗漏通道,通道壁上覆盖天然气水合物,形成一个管状通道,管状通道将游离气与孔隙水分离,游离气通过管状通道到达海底,形成羽状流.本文对琼东南活动冷泉的特征及形成模式提出思考,为更准确和深刻的研究冷泉提供研究思路.

致谢

文鹏飞教授级高工、张宝金教授级高工和张如伟高工给予了指导,马金凤高工对本文作图提供了帮助,陈江欣博士对本文提出了修改建议,对科考队员辛劳付出和匿名评审的细心审阅与宝贵意见在此一并感谢!

References
Bai Y, Song H B, Guan Y X, et al. 2014. Structural characteristics and genesis of pockmarks in the northwest of the South China Sea derived from reflective seismic and multibeam data. Chinese J. Geophys. (in Chinese), 57(7): 2208-2222. DOI:10.6038/cjg20140716
Bünz S, Polyanov S, Vadakkepuliyambatta S, et al. 2012. Active gas venting through hydrate-bearing sediments on the Vestnesa Ridge, offshore W-Svalbard. Marine Geology, 332-334: 189-197. DOI:10.1016/j.margeo.2012.09.012
Berndt C, Feseker T, Treude T, et al. 2014. Temporal constraints on hydrate-controlled methane seepage off Svalbard. Science, 343(6168): 284-287. DOI:10.1126/science.1246298
Brown G Jr, Sleeper K, Johnson M W, et al. 2013. Mercury concentrations, speciation, and isotopic composition in sediment from a cold seep in the northern Gulf of Mexico. Marine Pollution Bulletin, 77(1-2): 308-314. DOI:10.1016/j.marpolbul.2013.09.030
Chen D F, Li X X, Xia B. 2004. Distribution of gas hydrate stable zones and resource prediction in the Qiongdongnan Basin of the South China Sea. Chinese J. Geophys. (in Chinese), 47(3): 483-489.
Chen J X, Guan Y X, Song H B, et al. 2015. Distribution characteristics and geological implications of pockmarks and mud volcanoes in the northern and western continental margins of the South China Sea. Chinese J. Geophys. (in Chinese), 58(3): 919-938. DOI:10.6038/cjg20150319
Chen J X, Bai Y, Guan Y X, et al. 2016. Geophysical analysis of abnormal seismic (oceanography) reflection characteristics of oceanic bottom boundary layer. Chinese J. Geophys. (in Chinese), 59(6): 2148-2161. DOI:10.6038/cjg20160620
Chen J X, Song H B, Guan Y X, et al. 2017. A preliminary study of submarine cold seeps applying Seismic Oceanography techniques. Chinese J. Geophys. (in Chinese), 60(2): 604-616. DOI:10.6038/cjg20170215
Duperron S, Sibuet M, MacGregor B J, et al. 2007. Diversity, relative abundance and metabolic potential of bacterial endosymbionts in three Bathymodiolus mussel species from cold seeps in the Gulf of Mexico. Environmental Microbiology, 9(6): 1423-1438. DOI:10.1111/emi.2007.9.issue-6
Fannin N G T. 1980. The use of regional geological surveys in the North Sea and adjacent areas in the recognition of offshore hazards. //Ardus D A. Offshore Site Investigation. London: Graham & Trotman, 5-21.
Garcia-Gil S, Vilas F, Garcia-Garcia A. 2002. Shallow gas features in incised-valley fills (Ría de Vigo, NW Spain):A case study. Continental Shelf Research, 22(16): 2303-2315. DOI:10.1016/S0278-4343(02)00057-2
Gorman A R, Gorman W S, Holbrook M J, et al. 2002. Migration of methane gas through thehydrate stability zone in a low-flux hydrate province. Geology, 30(4): 327-330. DOI:10.1130/0091-7613(2002)030<0327:MOMGTT>2.0.CO;2
Guan J A, Fan S S, Liang D Q, et al. 2009. Methane hydrate growth velocity in seepage system of Qiongdongnan Basin, South China Sea. Chinese J. Geophys. (in Chinese), 52(3): 765-775.
He J X, Su P B, Lu Z Q, et al. 2015. Prediction of gas sources of natural gas hydrate in the Qiongdongnan Basin, northern South China Sea, and its migration, accumulation and reservoir formation pattern. Natural Gas Industry (in Chinese), 35(8): 19-29.
Hovland M, Judd A G. 1988. Seabed Pockmarks and Seepages:Impact on Geology, Biology and the Marine Environment. London: Graham and Trotman: 294.
Hsu S K, Wang S Y, Liao Y C, et al. 2013. Tide-modulated gas emissions and tremors off SW Taiwan. Earth and Planetary Science Letters, 369-370: 98-107. DOI:10.1016/j.epsl.2013.03.013
Huang B J. 2002. Genetic types and migration-accumulation dynamics of natural gases in the Ying-Qiong Basin, the South China Sea (in Chinese). Guangzhou: Guangzhou Institute of Geochemistry of Chinese Academy of Sciences.
Judd A, Hovland M. 2007. Seabed Fluid Flow:The Impact on Geology, Biology and the Marine Environment. Cambridge: Cambridge University Press: 475.
Karisiddaiah S M, Veerayya M, Vora K H, et al. 1993. Gas-charged sediments on the inner continental shelf off western India. Marine Geology, 110(1-2): 143-152. DOI:10.1016/0025-3227(93)90110-H
Kastner M, Bartlett D, MacDonald I, et al. 2005. CH4 fluxes across the seafloor at three distinct gas hydrate fields: impacts on ocean and atmosphere chemistry. //Proceedings of the Fifth International Conference on Gas Hydrates(ICGH5). Trondheim, Norway.
Leifer I, MacDonald I. 2003. Dynamics of the gas flux from shallow gas hydrate deposits:interaction between oily hydrate bubbles and the oceanic environment. Earth Planetary Science Letter, 210(3-4): 411-424. DOI:10.1016/S0012-821X(03)00173-0
Leven J, Middleton M. 2007. Near-surface seismic expression of gas chimneys in the Perth Basin. //ASEG Extended Abstracts. ASEG.
Liu X L, Flemings P B. 2006. Passing gas through the hydrate stability zone at southern Hydrate Ridge, offshore Oregon. Earth and Planetary Science Letters, 241(1-2): 211-226. DOI:10.1016/j.epsl.2005.10.026
Logan G A, Jones A T, Kennard J M, et al. 2010. Australian offshore natural hydrocarbon seepage studies, a review and re-evaluation. Marine and Petroleum Geology, 27(1): 26-45. DOI:10.1016/j.marpetgeo.2009.07.002
Ma L W, Gu H M, Li Z J, et al. 2015. Simulation of carbonate fracture-cavern reservoir reflection characteristics with forward modeling. OGP (in Chinese), 50(2): 290-297.
Ma X C, Luan Z D, Zhang X, et al. 2017. Near-bottom topography measurement using ROV and its application in a deep-sea hydrothermal field in the Manus Basin. Acta Oceanologica Sinica (in Chinese), 39(3): 76-84.
Milkov A V. 2004. Global estimates of hydrate-bound gas in marine sediments:how much is really out there?. Earth Sci. Rev., 66(3-4): 183-197. DOI:10.1016/j.earscirev.2003.11.002
Ruppel C, Dicken G R, Castellini D G, et al. 2005. Heat and salt inhibition of gas hydrateformation in the northern Gulf of Mexico. Geophys. Res. Lett., 32: L04605. DOI:10.1029/2004GL021909
Schroot B M, Klaver G T, Schüttenhelm R T E. 2005. Surface and subsurface expressions of gas seepage to the seabed-examples from the Southern North Sea. Marine and Petroleum Geology, 22(4): 499-515. DOI:10.1016/j.marpetgeo.2004.08.007
Solomon E A, Kastner M, Robertson G, et al. 2005. Insights into the dynamics of in situ gas hydrate formation and dissociation at the Bush Hill gas hydrate field, Gull of Mexico. //Proceedings of the 5th International Conference on Gas Hydrates (ICGH5). Trondheim, Norway.
Taylor D I. 1992. Nearshore shallow gas around the U.K. coast. Continental Shelf Research, 12(10): 1135-1144. DOI:10.1016/0278-4343(92)90074-T
Torres M E, Wallmann K, Tréhu A M, et al. 2004. Gas hydrate growth, methane transport, and chloride enrichment at the southern summit of Hydrate Ridge, Cascadia margin off Oregon. Earth and Planetary Science Letters, 226(1-2): 225-241. DOI:10.1016/j.epsl.2004.07.029
Tréhu A M, Flemings P B, Bang N L, et al. 2004. Feeding methane vents and gas hydrate depositsat south Hydrate Ridge. Geophysical Research Letters, 31: L23310. DOI:10.1029/2004GL021286
Vielstädte L, Karstens J, Haeckel M, et al. 2015. Quantification of methane emissions at abandoned gas wells in the Central North Sea. Marine and Petroleum Geology, 68: 848-860. DOI:10.1016/j.marpetgeo.2015.07.030
Wang X J, Wu S G, Wang D W, et al. 2010. The role of polygonal faults in fluid migration and gas hydrate reservoir forming in Southeast Hainan Basin. OGP (in Chinese), 45(1): 122-128.
Wenau S, Spiess V, Pape T, et al. 2015a. Cold seeps at the saltfront in the Lower Congo Basin Ⅰ:Current methane accumulation and active seepage. Marine and Petroleum Geology, 67: 894-908. DOI:10.1016/j.marpetgeo.2014.07.032
Wenau S, Spiess V, Pape T, et al. 2015b. Cold seeps at the salt front in the Lower Congo Basin Ⅱ:The impact of spatial and temporal evolution of salt-tectonics on hydrocarbon seepage. Marine and Petroleum Geology, 67: 880-893. DOI:10.1016/j.marpetgeo.2014.09.021
Wood W T, Hart P E, Hutchinson D R, et al. 2008. Gas and gas hydrate distribution around seafloor seeps in Mississippi Canyon, Northern Gulf of Mexico, using multi-resolution seismic imagery. Marine and Petroleum Geology, 25(9): 952-959. DOI:10.1016/j.marpetgeo.2008.01.015
Wu S G, Sun Q L, Wu T Y, et al. 2009. Polygonal fault and oil-gas accumulation in deep-water area of Qiongdongnan Basin. Acta Petrolei Sinica (in Chinese), 30(1): 22-26, 32.
拜阳, 宋海斌, 关永贤, 等. 2014. 利用反射地震和多波束资料研究南海西北部麻坑的结构特征与成因. 地球物理学报, 57(7): 2208-2222. DOI:10.6038/cjg20140716
陈多福, 李绪宣, 夏斌. 2004. 南海琼东南盆地天然气水合物稳定域分布特征及资源预测. 地球物理学报, 47(3): 483-489.
陈江欣, 关永贤, 宋海斌, 等. 2015. 麻坑、泥火山在南海北部与西部陆缘的分布特征和地质意义. 地球物理学报, 58(3): 919-938. DOI:10.6038/cjg20150319
陈江欣, 拜阳, 关永贤, 等. 2016. 海底边界层异常(地震海洋学)反射地震特征的地球物理分析. 地球物理学报, 59(6): 2148-2161. DOI:10.6038/cjg20160620
陈江欣, 宋海斌, 关永贤, 等. 2017. 海底冷泉的地震海洋学初探. 地球物理学报, 60(2): 604-616. DOI:10.6038/cjg20170215
关进安, 樊栓狮, 梁德青, 等. 2009. 南海琼东南盆地渗漏系统甲烷水合物生长速度. 地球物理学报, 52(3): 765-775.
何家雄, 苏丕波, 卢振权, 等. 2015. 南海北部琼东南盆地天然气水合物气源及运聚成藏模式预测. 地质勘探, 35(8): 19-29.
黄保家. 2002. 莺琼盆地天然气成因类型及成藏动力学研究. 广州: 中国科学院广州地球化学研究所.
马灵伟, 顾汉明, 李宗杰, 等. 2015. 正演模拟碳酸盐岩缝洞型储层反射特征. 石油地球物理勘探, 50(2): 290-297.
马小川, 栾振东, 张鑫, 等. 2017. 基于ROV的近海底地形测量及其在马努斯盆地热液区的应用. 海洋学报, 39(3): 76-84.
王秀娟, 吴时国, 王大伟, 等. 2010. 琼东南盆地多边形断层在流体运移和天然气水合物成藏中的作用. 石油地球物理勘探, 45(1): 122-128.
吴时国, 孙启良, 吴拓宇, 等. 2009. 琼东南盆地深水区多边形断层的发现及其油气意义. 石油学报, 30(1): 22-26, 32. DOI:10.7623/syxb200901005