Paleo-environment and paleo-productivity of the hydrate reservoirs in the Shenhu area of northern South China Sea
-
摘要: 天然气水合物作为一种新型的清洁能源,其形成需要稳定的有机质供应。南海北部神狐海域为天然气水合物成藏的有利区域,2017年中国地质调查局在神狐海域水合物试采获得突破性成功。为了进一步了解古环境和古生产率对形成水合物有机质供应的影响,对南海北部神狐海域水合物钻探区SH3站位180~215mbsf(meters below the sea floor)层位,尤其是水合物主要赋存层位190~200mbsf的古环境、古生产率以及陆源碎屑物质的地球化学指标进行分析研究。研究表明水合物主要赋存层位陆源碎屑物质(TDM)输入增加,较高的陆源碎屑物质输入和次氧化的沉积环境共同造就了比较好的有机质的外部保存条件;同时较强的水动力条件,有利于藻类生物的繁殖,因此,生物成因有机质比较丰富, 再加上神狐海域有比较良好的热解气的形成条件,这3个层面共同保证了神狐海域具有比较充足的有机质供应。Abstract: The accumulation of natural gas hydrate requires a stable supply of organic gas. The Shenhu area of the northern part of the South China Sea has been proved a favorable area for natural gas hydrates generation. In the year of 2017, the trial case of hydrate production made a breakthrough there. To further understand the influence of the paleo-environment and paleo- productivity onto organic preservation, this paper is specially devoted to the interval of 180-215 meters below the sea floor at the station of SH3, especially the horizon of 190-200mbsf of the main hydrate storage in the Shenhu area. Studies on the Paleo-environment, redox environment and terrigenous detrital material input suggests that the input of terrigenous detrital material (TDM) increase in the hydrate-bearing horizons and the suboxic depositional environments. The two factors jointly created a favorable external condition for storage and preservation of organic matters. At the same time, the stronger hydrodynamic conditions are conducive to the production of algae, and the biogenic organic gas is relatively abundant. Together with the better conditions for the formation of thermogenic gas, there are enough gas sources for the formation of gas hydrate deposits in the Shenhu area.
-
Keywords:
- paleo-productivity /
- paleo-environments /
- Shenhu area /
- gas hydrate
-
地震台站的背景噪声研究是地震学研究的一个重要课题[1-3]。Peterson选取了全球地震台网中的75个地震台,利用它们记录的背景噪声数据,构建了著名的全球背景噪声高值模型和低值模型[4]。Mcnamara和Buland在传统的噪声功率谱密度分析方法的基础上提出了更全时段地分析地震台站噪声的概率密度函数统计分析方法[5]。自此,通过概率密度函数统计分析方法获取台站数据的功率谱密度的概率分布特征,并与全球背景噪声新的高模型和低模型对比,已经成为评估台站周围环境背景噪声水平的最常用手段[6-9]。
南海东部次海盆一直是海洋地质与地球物理研究的热点地区[10-13]。然而,海域地震台站的缺失导致天然地震学方法在此区域的应用一直存在很大的限制。例如,前人的研究多依赖南海周边的海南岛、吕宋岛和中南半岛上的陆基台站接收到的远场地震数据,故天然地震学手段的应用仅限于层析成像等方法[14-15]。而像接收函数分析、背景噪声分析这样依赖地震台站近场接收数据的研究是无法开展的。海底地震仪的出现大大拓展了直接在深海大洋中进行天然地震观测的能力[16-18]。海底地震仪在海床上记录到的信号包括各种类型的地震信号、背景噪声信号等。因此,海洋背景噪声水平成为影响海底地震仪地震观测质量的一个重要因素[19]。了解海洋背景噪声的分布特征对于评估地震观测质量及后续数据处理如何降低噪声均具有重要指示意义。
国家海洋局第一海洋研究所和同济大学共同实施了南海大规模的被动源海底地震仪台阵长期观测实验。Liu等主要对台站HY02的微震频段的背景噪声进行了概率密度函数统计分析,并将其众数谱曲线与陆地和洋岛台站进行了对比[16]。Xiao等则利用部分数据对微震频带的背景噪声进行了功率谱密度分析,并对其影响因素进行了讨论[20]。可以看出上述研究分析均聚焦于背景噪声的微震频段,缺乏对其他频段、其他台站、特殊地震信号等的讨论和对比。因此,本研究利用上述实验的部分数据对南海的背景噪声在全频段进行了讨论,并与其他台站做了对比;同时对观测过程中出现的地震以及其他信号的概率密度分布进行了归纳总结;最后研究了背景噪声的时间变化规律。
1. 数据与方法
国家海洋局第一海洋研究所和同济大学共同实施了南海第一次大规模的被动源海底地震仪长期观测实验。2012年4月,项目组利用国家自然科学基金船时共享计划,由“东方红2”科学考察船实施了海底地震仪布放航次。该实验围绕黄岩-珍贝海山链呈网状布放了18台海底地震仪,绝大多数台站的布放水深超过了3500 m。2013年4月,项目组利用广东海洋大学的“天龙”科学考察船组织了回收航次,成功地回收了海底地震仪。
本研究利用了该实验中的三台海底地震仪的观测数据(台站编号:HY15、HY16、HY17),三台海底地震仪的型号均为Guralp CMG-40T。布放位置均为残留扩张脊北侧,水深超过3500 m(见图1和表1)。台站的间距约70 km。海底地震仪虽然在海底布放了一年的时间,但受电池容量的影响,它们的有效工作时间均只有两百多天。
表 1 南海东部次海盆台站布放资料Table 1. The deployment data of OBS台站编号 东经 北纬 水深/m 台站型号 HY15 117°32′ 16°30′ 3753 Guralp CMG-40T HY16 118°12′ 16°27′ 3920 Guralp CMG-40T HY17 118°48′ 16°12′ 3870 Guralp CMG-40T 海底地震仪的三分量数据首先去仪器响应转换为速度记录。将连续的速度记录分割为24 h的数据段;再分为长度为1 h的数据段,按照数据段的重合率为50%截取,这样一天将有48段数据;再将1 h的数据段分成长度1000 s的数据段,按照数据段的重合率为80%截取,由此1 h可以得到14段数据;对每小段数据去均值和线性趋势,最后计算这一小段时间序列的速度功率谱密度。重复计算每小段的速度功率谱密度并取平均值,从而得到速度功率谱密度随频率的分布[7]。功率谱密度的计算过程如下。
周期时间序列
$g(t)$ 的傅里叶变换可表示为:$$G(f,n) = \int\limits_0^n {g(t){{\rm{e}}^{ - i2{\rm{\pi }}ft}}{\rm{d}}t} $$ 其中,
$n$ 为时间序列的长度,$f$ 为频率。对于离散频率值${f_k}$ ,可定义为:$${G_k} = \frac{{G\left( {{f_k},n} \right)}}{{\Delta t}}$$ $$ {f_k} = {k / {\left( {N\Delta t} \right)}}\;\;\;\;k = 1,2, \cdots ,N $$ $\Delta t$ 为采样间隔,$N$ 为截取时间段的采样点数。功率谱密度定义为:$${P_k} = \frac{{2\Delta t}}{N}{\left| {{G_k}} \right|^2}$$ 为了对速度功率谱密度进行充分采样,我们按照1/8倍频程间隔对整个频率周期进行重采样处理。功率值在短周期
${T_{\rm{s}}}$ 与长周期${T_{\rm{l}}} = 2*{T_{\rm{s}}}$ 之间进行平均,对应的周期${T_{\rm{c}}}$ 为倍频程内的几何平均值:$${T_{\rm{c}}} = \sqrt {{T_{\rm{s}}}{T_{\rm{l}}}} $$ ${T_{\rm{s}}}$ 按1/8倍频程增加,${T_{\rm{s}}} = {T_{\rm{s}}}*{2^{{1 / 8}}}$ ,用于计算下一个间隔的平均功率。重复这一过程,持续至原始资料时间序列长度最长的有效周期。为了便于后期与NHNM和NLNM进行对比分析,需要将速度功率谱密度转化为加速度功率谱密度,并将加速度功率谱密度的单位转变成分贝:$${\rm{PSD}} = 10\log (4{{\text π}^2}f_k^2{P_k})$$ 最后,以1 dB为步长进行概率统计,统计范围为−200~−50 dB,
${N_f}$ 是以$f$ 为中心频率的所有功率谱密度的个数,${N_{df}}$ 是频率$f$ 处功率谱密度在$d \sim d + {\rm{1}}$ dB范围内的个数,则$\left( {f,d} \right)$ 处的概率密度为$$P\left( {f,d} \right) = \frac{{{N_{df}}}}{{{N_f}}}$$ 由此统计所有的频率,得到功率谱密度在−200~−50 dB上的分布。
2. 背景噪声特征
通过上述计算获得了台站HY15的三分量噪声功率谱密度的概率密度函数分布图,如图2所示。直观来看,在全频带范围内,两水平分量的背景噪声水平是要高于垂直分量的。造成这种差别的主要原因是,本研究所用海底地震仪的水平分量与海底的耦合程度受洋流或地形的影响较垂直分量差[21-23]。Guralp CMG-40T 型海底地震仪记录的频带宽度为50 Hz~60 s,故将其频带宽度划分为3个范围来分析:高频段(0.025~1 s,1~40 Hz)、微震段(1~20 s,0.05~1 Hz)、长周期段(20~50 s,0.02~0.05 Hz)。
在长周期低频段,三分量的噪声功率谱密度都呈现出一个明显高于NHNM的宽缓的峰值,且两个水平分量的水平明显高于垂直分量。这一频带的噪声主要受洋流和极低频海浪的影响。洋流主要通过与海底地震仪相互作用产生噪声,例如,其流过海底地震仪时会形成湍流从而影响水平分量。极低频海浪会作用于海底,使之变形产生柔性噪声从而影响垂直分量,而海底地震仪的倾斜又会使这些噪声泄漏到水平分量中[24]。极低频海浪产生的次重力波仅在波长大于水深h的频率下才是显著的:h=100 m时低于0.12 Hz,h=1000 m时低于0.04 Hz,h=4000 m时低于0.02 Hz[22]。
在微震段,0.07 Hz左右的单频微震峰的主要成因与海浪和近岸浅水区的相互撞击产生的地震能量有关,它的频率与产生它的海浪的频率相同,在远离陆地的公海上较为微弱[5]。台站HY15距离吕宋岛的海岸线较近,故此频率处的背景噪声水平高于NHNM。在0.3 Hz左右高于NHNM的最大峰值被称为双频微震峰,命名规则是因为它的观测频率两倍于产生它的海浪的频率。Longuet-Higgins解释了其成因,沿相反方向传播的等周期的海浪叠加而产生半个标准海浪波周期的重力驻波,静止的重力波在水中引起扰动,这些扰动传播到洋底就是双频微震噪声[25]。双频微震噪声有本地源(短周期,0.2~0.45 Hz)和远源(长周期,0.085~0.2 Hz)[26-27]。本台站在微震段以0.2 Hz为界,小于0.2 Hz的噪声迅速降低意味着远源双频微震的缺失,主要原因是南海作为一个封闭的边缘海与产生远源双频微震的太平洋的联系较弱;大于0.2 Hz的双频微震噪声呈现出大的峰值,这表明本地的海浪是产生微震的主要源头。
在高频段,背景噪声水平处于NHNM和NLNM之间,在一个正常的水平范围内。高频段的背景噪声主要是人类活动引起的噪声,例如公路、铁路、工厂等产生的噪声[7],这些噪声正是深海大洋中极度缺失的。
为了更好地分析南海东部次海盆的背景噪声水平,将台站HY15和与之距离较远的台站HY17进行了对比分析。由图3可知,在全频段上,台站HY17与台站HY15的背景噪声水平相近,表明由海底地震仪得到的背景噪声的覆盖范围是可以代表整个研究区的。但在一些频率处,两个台站仍有微小的差异波动,例如长周期段和微震段。两个台站处于同一个研究区,故台站所处的具体位置成为控制上述差异的主要因素。例如,台站HY17相比台站HY15距离海岸线更近,更容易受海浪拍击海岸的影响;台站HY17位于海山区,故相比于海盆区的台站HY15,其崎岖的海底地形会减弱海底洋流对海底地震仪的影响。
NHNM和NLNM是通过全球陆基地震台站的数据得出的,故再选取中国内陆的陕西安康台进行对比分析。由图4可知,安康台三个分量的单频微震噪声呈现出明显的峰值,双频微震的远源噪声也与NHNM和NLNM的峰值有较好的对应关系。在全频带上,安康台三个分量的背景噪声水平远远低于NHNM。本研究所用的海底地震仪与之相比,背景噪声水平明显偏高,差值最高可达80 dB。另外,海底地震仪的噪声功率谱的变化幅度明显宽于安康台,安康台的变化范围仅在众数谱线的上下10 dB内。这都表明海底地震仪的环境噪声要比陆基地震台站复杂的多,海底地震仪记录数据的质量要比陆基地震台站差的多。作为陆基固定地震台的陕西安康台,其选址、安装、入网是有严格标准的。而海底地震仪布放位置的选择往往只是考虑能否实现预期的科学目的,对于海底的观测环境的选择往往是无能为力的。
3. 典型信号分析
传统的功率谱密度分析方法是选择一段没有地震发生、没有人文噪声产生的平静记录来评估台站的噪声水平[4, 28]。因此,所谓的合理选择往往存在人为判断因素,其结果应是局部时段有代表性的而不是全观测周期的反映台站的噪声水平,例如,非周期性的(地震和人为噪声)和周期性的(季节性的)突发干扰容易被人为剔除掉[6]。而上述概率密度函数分析方法在计算过程中不必排除地震等突发事件,不需要挑选连续平静的噪声记录,而是将记录到的连续时间信号进行统一处理,这些扰动将作为低概率事件存在于最后的统计分布中[7]。
Guralp CMG-40T OBS较宽的频带范围使其记录了不同周期范围的信号。因此,我们筛选了一些典型的时间序列进行了分析讨论。其中,图5展示了属于远震事件的2012年8月14日发生于鄂霍茨克海的7.7级大地震,各个震相(P、S)可以清楚分辨;图6展示了2012年10月8日发生于班达海的一个近震事件。另外,海底地震仪由于自身的原因也会产生一些特殊信号,如图7所示,海底地震仪发生数据丢失现象,波形被填充为零。
由于上述波形示例的采样点数较少,将计算过程的参数进行了调整以获取正确的功率谱密度,用于揭示它们在整个概率密度分布中的位置。如图8所示,远震事件主要影响长周期低频段和微震段,与之相反的是,近震事件主要影响高频段。而数据丢失意味着没有记录任何噪声,其噪声水平在合理的范围之内,且功率谱的变化幅度较窄。
图 8 远震事件、近震事件和数据丢失现象的功率谱密度的概率密度函数分布图其中,上图为远震事件;中图为近震事件;下图代表数据丢失信号;两条蓝色实线分别代表全球背景噪声高值和全球背景噪声低值;黑色实线代表众数值。Figure 8. The PDF of teleseismic, near earthquake and data dropout signal The top panel represents teleseismic, the middle panel represents near earthquake, the bottom panel represents data dropout signal. The two blue solid lines represent NHNM and NLNM. The black solid line represents the mode value.通过上一节的分析可知,海底地震仪的外部观测环境是难以改变的,理想的滤波频带选择成为数据降噪的唯一手段。因此,上述远震、近震事件功率谱的概率分布的确定对于后续滤波处理具有重要指示意义。例如,依据远震事件高于NHNM的频带范围,可以确定用于分离远震事件的带通滤波器的频率为0.02~0.15 Hz。另外,远震、近震和数据丢失现象的出现在整个观测周期内属于小概率发生事件,因此,其对整体的概率密度分布的影响是有限的,这对于数据质量的检查是具有重要意义的。例如,在不顺序检查所有记录波形的情况下,如果海底地震仪观测期内整体的概率密度分布接近数据丢失现象,是有理由怀疑此海底地震仪的数据记录出现了故障。
4. 时间变化特征
将观测期内每一天的噪声功率谱密度按照日期顺序排列,可以获得台站观测期内噪声功率谱密度同时随时间和频率的变化图。图9展示了台站HY15和HY17连续记录了270多天的三分量功率谱密度时间图。可以发现,背景噪声在高频段几乎不随时间发生变化,而在微震段和长周期低频段却都显示出很大的时间变化,例如,儒略日267天存在一个明显的背景噪声增大。这揭示出影响这种变化的因素并不是人文噪声,而是一种可以导致海浪或者洋流发生扰动的自然因素。另外,这种时间变化的持续时间往往是一周左右,而两次变化的间隔则长短不一。对比台站HY15和HY17可知,两个台站的时间变化几乎是同步发生的。从图1可知,两台站的间距非常大,这揭示出造成这种变化的因素的影响范围是非常巨大的。而在海洋中能够造成这种时间变化尺度和空间影响范围的自然因素只有台风。南海是传统的台风区,当台风影响南海或者过境南海的时候带来的巨大风浪增大了背景噪声水平。Xiao等基于南海海盆的海底地震仪台站和周边的陆上地震仪记录的数据以及全球浪高模型,利用功率谱密度分析、相关分析和极化分析,对南海及其周边地区的微地动噪音源进行了定位,对噪音的成分和传播进行了约束[20]。他们的研究结果明确了台风对于微震的影响,而本研究通过对背景噪声时间变化特征的分析讨论也将影响因素聚焦于台风。
图 9 台站HY15和HY17的背景噪声功率谱密度的时间变化图其中上面3幅图为台站HY15,下面3幅图为台站HY17,从左到右分别为垂直分量、水平分量1、水平分量2。Figure 9. The map of time variation of PSD of background nose. The above three maps are results of station HY15 From left to right, they are vertical component, horizontal component 1 and horizontal component 2 respectively. The following three pictures are results of station HY17.5. 结论
本研究基于南海东部次海盆的海底地震仪记录的数据,利用功率谱密度分析方法、概率密度函数统计分析方法对南海东部次海盆的背景噪声进行了研究和讨论。
(1)布放于南海东部次海盆的海底地震仪的背景噪声水平在微震段和长周期低频段高于全球背景噪声高值模型,在全频段上高于陆基地震台的背景噪声水平,这揭示出海底地震仪记录数据质量较陆基地震台站的数据质量差。
(2)受背景噪声的影响,远震事件的优势频带主要为长周期低频段和微震段,近震事件的优势频带则为高频段,数据丢失现象则表现为其功率谱密度位于全球背景噪声高、低值模型之间,这对数据的降噪处理、数据的质量分析提供了重要指示。
(3)南海东部次海盆的背景噪声在微震段和长周期低频段具有明显的时间变化特征,此变化的持续时间、间隔和空间影响范围揭示出台风是造成这一时间变化现象的主要原因。
致谢:感谢在被动源海底地震仪探测实验的布放航次和回收航次给予帮助的单位和个人。
致谢: 本次研究样品由广州海洋地质调查局于2007年在南海北部陆坡神狐海域实施“中国海域天然气水合物钻探”(GMGS1)航次获得,感谢该航次科学家们为研究样品的获得付出的辛勤劳动;感谢中国科学院广州能源研究所分析测试中心、武汉上谱分析科技有限责任公司和中国科学院广州地球化学研究所为本文研究提供的测试。 -
图 1 神狐海域钻探区分布图及SH3站位分布位置地形地质图(改编自文献[15])
Figure 1. Tectonic map of subareas of the Shenhu drilling area. Note the location of SH3
表 1 神狐海域钻探区水合物油气体系(SH2、SH3和SH7站位的地球化学数据据文献[3, 25])
Table 1 Characteristics of the gas-hydrate petroleum system in the Shenhu drilling area highlighting the distinct differences between the three sites
站位 水合物稳定条件 气体组成以及气体来源 气体运移体系 SH2 海底温度:4.84℃
海底压力:12.46MPa
水深:1245mCH4:99.89%;δ13C1:-56.7%,生
物气为主导的混合气气烟囱伴有很多微小的断层;低甲
烷通量SH3 海底温度: 5.53℃
海底压力: 12.61MPa
地热梯度: 49.34℃/kmδ13C1: -61.6%, 生物气为主导的
混合气气烟囱伴有很多微小的断层;低甲
烷通量SH7 海底温度: 6.44℃
海底压力: 11.20MPa
地热梯度: 43.65℃/km生物气为主导的混合气 微小断层;低甲烷通量 表 2 主、微量元素随深度的分布
Table 2 Distribution of major and trace elements with depth
深度/mbsf Al/% Si/% Ti/% P/% TOC/% V/×10-6 Cr/×10-6 Ni/×10-6 Th/×10-6 U/×10-6 Mo/×10-6 Co/×10-6 183.52 6.977 21.285 0.476 0.070 0.181 95.959 69.624 39.928 11.778 2.213 0.242 13.5 190.69 6.584 22.489 0.455 0.080 0.464 85.453 65.464 36.708 10.628 2.246 0.819 12.3 190.86 6.572 21.073 0.436 0.075 0.171 81.188 61.827 71.951 10.444 2.142 0.759 26.2 191.05 6.813 21.706 0.451 0.069 0.175 85.115 64.530 37.070 10.908 2.163 0.278 12.5 193.16 5.490 19.593 0.409 0.095 0.145 70.181 53.829 30.255 9.022 2.304 0.268 10.4 193.49 5.513 20.501 0.403 0.080 0.137 67.611 50.653 34.531 9.225 2.125 0.603 9.9 196.02 6.194 22.299 0.436 0.080 0.166 78.922 61.643 36.370 9.842 2.058 0.668 12.7 201.29 6.943 21.585 0.451 0.073 0.193 90.242 66.371 36.960 10.909 2.134 0.270 13.0 201.33 6.845 20.971 0.470 0.072 0.170 88.183 64.591 40.737 10.650 2.130 0.356 13.9 211.62 7.457 22.312 0.453 0.071 0.197 89.906 68.959 37.797 11.755 2.217 0.300 13.1 211.66 7.405 22.836 0.477 0.071 0.198 91.907 70.742 38.436 11.959 2.200 0.245 13.4 最小值 5.490 19.593 0.403 0.069 0.137 67.611 50.653 30.255 9.022 2.058 0.242 9.880 最大值 7.457 22.836 0.477 0.095 0.464 95.959 70.742 71.951 11.959 2.304 0.819 26.179 平均值 6.618 21.514 0.447 0.076 0.200 84.061 63.476 40.068 10.647 2.176 0.437 13.702 表 3 氧化物指标及各元素比值
Table 3 Oxide indexes and ratios of elements
深度/mbsf Al2O3 /% TiO2 /% SiO2/% Al/Ti P/Al p/Ti SiO2/Al2O3 V/Cr V/Ni U/Th U/Mo Ni/Co 183.52 13.179 0.793 45.676 14.658 0.010 0.147 3.466 1.378 2.403 0.188 9.152 2.947 190.69 12.437 0.758 48.260 14.470 0.012 0.176 3.880 1.305 2.328 0.211 2.742 2.995 190.86 12.414 0.727 45.221 15.073 0.011 0.172 3.643 1.313 1.128 0.205 2.822 2.748 191.05 12.869 0.752 46.579 15.106 0.010 0.153 3.619 1.319 2.296 0.198 7.789 2.976 193.16 10.370 0.682 42.045 13.423 0.017 0.232 4.054 1.304 2.320 0.255 8.596 2.914 193.49 10.414 0.672 43.994 13.680 0.015 0.199 4.225 1.335 1.958 0.230 3.526 3.494 196.02 11.700 0.727 47.852 14.206 0.013 0.183 4.090 1.280 2.170 0.209 3.082 2.870 201.29 13.115 0.752 46.320 15.395 0.011 0.162 3.532 1.360 2.442 0.196 7.907 2.850 201.33 12.930 0.783 45.002 14.564 0.011 0.153 3.481 1.365 2.165 0.200 5.985 2.925 211.62 14.086 0.755 47.880 16.461 0.010 0.157 3.399 1.304 2.379 0.189 7.392 2.890 211.66 13.988 0.795 49.004 15.524 0.010 0.149 3.503 1.299 2.391 0.184 8.980 2.868 min 10.370 0.672 42.045 13.423 0.010 0.147 3.399 1.280 1.128 0.184 2.742 2.748 max 14.086 0.795 49.004 16.461 0.017 0.232 4.225 1.378 2.442 0.255 9.152 3.495 ave 12.500 0.745 46.167 14.778 0.012 0.171 3.717 1.324 2.180 0.206 6.179 2.953 -
[1] Collett T S. Natural gas hydrate of the Prudhoe bay and Kuparuk River area, north slope Alaska [J]. AAPG Bulletin, 1993, 77(5): 793-812.
[2] Collett T S, Dallimore S R. Detailed analysis of gas hydrate roceinduced drilling and production hazardsd [R]. Proceedings of the 4th International Conferenceon Gas Hydrates, 2002, 73: 47-52.
[3] Wu N Y, Zhang H Q, Yang S X, et al. Gas hydrate system of Shenhu area, northern South China Sea: geochemical results [J]. Journal of Geological Research, 2011, 27: 1-10.
[4] Wu J F, Sunda W, Boyle E A, et al. Phosphate depletion in the western North Atlantic Ocean [J]. Science, 2000, 289: 759-762. doi: 10.1126/science.289.5480.759
[5] Miao Q, Thunell R C. Late Pleistocene-Holocene distribution of deep-sea benthic foraminifera in the South China Sea[R]: Paleoceanography implications[J]. Journal of Foraminifera Research, 1996, 26(1): 9-23.
[6] Zhang S L, Chen D F, Huang J Q. Conditions of accumulation of gas hydrate in Baiyun Sag [J]. Natural Gas Industry, 2007, 27: 7-10. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=trqgy200709002
[7] 吴能友, 邬黛黛, 陈雪刚, 等.南海东北部东沙海域沉积物烃类有机地球化学研究及其意义[J].南海地质研究, 2007, 23: 1-14. http://www.cqvip.com/Main/Detail.aspx?id=29623366 WU Nengyou, WU Daidai, CHEN Xuegang, et al. Hydrocarbon organic geochemistry of sediments in the Dongsha area of the northeastern South China Sea and its significance[J]. South China Sea Geological Research, 2007, 23: 1-14. http://www.cqvip.com/Main/Detail.aspx?id=29623366
[8] 陈多福, 姚伯初, 赵振华, 等.珠江口和琼东南盆地天然气水合物形成和稳定分布的地球化学边界条件及其分布区[J].海洋地质与第四纪地质, 2001, 21(4): 73-78. http://hydz.chinajournal.net.cn/WKD/WebPublication/paperDigest.aspx?paperID=b820ad45-6ebd-45ea-b8e7-2f0f5a143172 CHEN Duofu, YAO Bochu, ZHAO Zhenghua, et al. Geochemical boundary conditions and distribution of natural gas hydrate formation in the Pearl River Estuary and Qiongdongnan basin[J]. Marine Geology and Quaternary Geology, 2001, 21(4): 73-78. http://hydz.chinajournal.net.cn/WKD/WebPublication/paperDigest.aspx?paperID=b820ad45-6ebd-45ea-b8e7-2f0f5a143172
[9] 邬黛黛, 吴能友, 张美, 等.东沙海域SMI与甲烷通量的关系及对水合物的指示[J].地球科学-中国地质大学学报, 2013, 38(6): 1309-1320. http://d.old.wanfangdata.com.cn/Periodical/dqkx201306014 WU Daidai, WU Nengyou, ZHANG Mei, et al. Relationship between SMI and methane fluxes in the Dongsha sea and their indications for hydrates [J]. Earth Science-Journal of China University of Geosciences, 2013, 38(6): 1309-1320. http://d.old.wanfangdata.com.cn/Periodical/dqkx201306014
[10] 刘昌玲, 孟庆国, 李承峰, 等.南海北部陆坡天然气水合物及其赋存沉积物特征[J].地学前缘, 2017, 24(4): 41-50. http://d.old.wanfangdata.com.cn/Periodical/dxqy201704005 LIU Changling, MENG Qingguo, LI Chengfeng, et al. Characteristics of gas hydrates and sedimentary sediments in the northern slopes of the South China Sea[J]. Earth Science Frontier, 2017, 24(4): 41-50. http://d.old.wanfangdata.com.cn/Periodical/dxqy201704005
[11] 苏丕波, 梁金强, 付少英, 等.南海北部天然气水合物成藏地质条件及成因模式探讨[J].中国地质, 2017, 44(3): 415-427. http://d.old.wanfangdata.com.cn/Periodical/zgdizhi201703001 SU Pibo, LIANG Jingqiang, FU Shaoying, et al. Geological conditions and genetic models of gas hydrate accumulation in the northern South China Sea[J]. Chinese Geology, 2017. 44(3): 415-427. http://d.old.wanfangdata.com.cn/Periodical/zgdizhi201703001
[12] Archer D, Buffett B, Brovkin V. Ocean methane hydrates as a slow tipping point in the global carbon cycle [J]. Proceedings of the National Academy of Sciences of the United States of America-Physical Sciences, 2009, 106: 20596-20601. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=af0285ca07b3ed347305c74d17785596
[13] Tyson R V. Sedimentation rate, dilution, preservation and total organic carbon: some results of a modelling study [J]. Organic Geochemistry, 2001, 32: 333-339. doi: 10.1016/S0146-6380(00)00161-3
[14] Wan Z X, Xu X, Wang B, et al. Geothermal analysis of boreholes in the Shenhu gas hydrate drilling area, northern South China Sea: Influence of mud diapirs on hydrate occurrence [J]. Journal of Petroleum Science and Engineering, 2017, 158: 424-432. doi: 10.1016/j.petrol.2017.08.053
[15] Ding X J, Liu G D, Lu X J, et al. The impact on organic matter preservation by the degree of oxidation and reduction and sedimentation rates of Erlian Basin [J]. Natural Gas Geoscience, 2014, 25: 810-817. http://d.old.wanfangdata.com.cn/Periodical/trqdqkx201406003
[16] Li W H, Zhang Z H, Li Y C. Some aspects of excellent marine source rock formation: implications on enrichment regularity of organic matter in continental margin basins [J]. Chinese Journal of Geochemistry. 2015, 34: 47-54. doi: 10.1007/s11631-014-0018-2
[17] Tyson R V, Pearson T H. Modern and ancient continental shelf anoxia: An overview [J]. Geological Society of London Special Publication, 1991, 58: 1-24. doi: 10.1144/GSL.SP.1991.058.01.01
[18] 董军社, 陈平富.南海北部陆架晚第三纪δ13C记录与古生产力[J].海洋地质与第四纪地质, 1997, 17(1): 40-44. http://hydz.chinajournal.net.cn/WKD/WebPublication/paperDigest.aspx?paperID=ba13178e-487a-4690-9faa-32257f5c7506 DONG Junshe, CHEN Pingfu. Neogene carbon isotope records from the northern south china sea and their bearing on paleoproductivity. [J]. Marine Geology and Quaternary Geology, 1997, 17(1): 40-44. http://hydz.chinajournal.net.cn/WKD/WebPublication/paperDigest.aspx?paperID=ba13178e-487a-4690-9faa-32257f5c7506
[19] 黄宝琦.南海东北部晚第四纪古生产力变化[J].海洋地质与第四纪地, 2000, 20(2): 65-68. http://d.old.wanfangdata.com.cn/Periodical/hydzydsjdz200002011 HUANG Baoqi. Changes of Late Quaternary productivity in the northeastern South China Sea[J]. Marine Geology and Quaternary Geology, 2000, 20(2): 65-68. http://d.old.wanfangdata.com.cn/Periodical/hydzydsjdz200002011
[20] 李建, 王汝建.南海北部一百万年以来的表层古生产力变化:来自ODP1144站的蛋白石记录[J].地质学报, 2004, 78(2): 228-233. doi: 10.3321/j.issn:0001-5717.2004.02.012 LI Jian, WANG Rujian. Changes in surface productivity in the northern South China Sea over 1 million years: Opal records from ODP 1144 station [J]. Geological Journal, 2004, 78(2): 228-233. doi: 10.3321/j.issn:0001-5717.2004.02.012
[21] 李鹤, 黄宝琦, 王娜.南海北部MD12-3429站位海水古生产力和溶解氧含量特征[J].古生物学报, 2017, 56(2): 238-248. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gswxb201702010 LI He, HUANG Baoqi, WANG Na. Palo- Productivity and dissolved oxygen contents in the seawater of MD12-3429 station in the northern South China Sea[J]. Acta Palaeontologica Sinica, 2017, 56(2): 238-248. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gswxb201702010
[22] 吴聪, 陈炽新, 陈芳, 等.南海北部陆坡SH1B孔MIS6期以来的古生产力和古环境变化的硅藻记录[J].海洋地质与第四纪地质, 2015, 35(5): 95-102. http://hydz.chinajournal.net.cn/WKD/WebPublication/paperDigest.aspx?paperID=09fe8f87-0723-4897-b010-b9c351f3b2c6 WU Cong, CHENG Zhixing, CHENG Fang, et al. The paleo-productivity and paleo-environmental change of diatom records since the MIS6 period of the SH1B hole in the northern slope of the South China Sea [J]. Marine Geology and Quaternary Geology, 2015, 35(5): 95-102. http://hydz.chinajournal.net.cn/WKD/WebPublication/paperDigest.aspx?paperID=09fe8f87-0723-4897-b010-b9c351f3b2c6
[23] Zhu Y H, Su P B. Gas sources of natural gas hydrates in the Shenhu drilling area, South China Sea: geochemical evidence and geological analysis [J]. Journal of the Geological Society of China, 2013, 87(3): 767-776. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dzxb-e201303014
[24] Li S J, Chu F Y, Fang Y X, et al. Associated interpretation of sub-bottom and single-channel seismic profiles from slope of Shenhu Area in the northern South China Sea: Characteristics of gas hydrate sediment [J]. Journal of Tropical Oceanography, 2010, 29(4): 56-62. http://www.jto.ac.cn/en/y2010/v29/i4/56
[25] Su M, Yang R, Wu N Y, et al. Structural characteristics in the Shenhu area, northern continental slope of South China Sea, and their influence on gas hydrate [J]. Acta Geologica Sinica, 2014, 88(3): 318-326. http://en.cnki.com.cn/Article_en/CJFDTOTAL-DZXE201403002.htm
[26] Liang J Q, Wang H B, Su X, et al. Natural gas hydrate formation conditions and the associated controlling factors in the northern slope of the South China Sea [J]. Natural Gas Industry, 2014, 34(7): 128-135. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=trqgy201407022
[27] Zhong Z H, Shi H S, Zhu M, et al. A discussion on the tectonic-stratigraphic framework and its origin mechanism in Pearl River Mouth Basin [J]. China Offshore Oil and Gas, 2014, 26(5): 20-29. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zghsyq-gc201405004
[28] Zhu W L, Huang B J, Mi L J, et al. Geochemistry, origin, and deep-water exploration potential of natural gases in the Pearl River Mouth and Qiongdongnan basins, South China Sea [J]. American Association of Petroleum Geologists Bulletin, 2009, 93: 741-761. doi: 10.1306/02170908099
[29] Lin C C, Lin A S, Liu C S, et al. Canyon-infilling and gas hydrate occurrences in the frontal fold of the offshore accretionary wedge off southern Taiwan [J]. Marine Geophysical Research, 2013, 12: 11-18. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=5ac0b035e8ba97e25d4a8c2a9f5dc2e8
[30] Liu C S, Huang I L, Teng L S. Structural features off southwestern Taiwan [J]. Marine Geology, 1997, 137: 305-319. doi: 10.1016/S0025-3227(96)00093-X
[31] Suess E, Huang Y Y, Wu N Y, et al. Cruise report SO-177, Sino-German cooperative project, South China Sea. in: distribution, formation and effect of methane & gas hydrate on the environment[R]. 2004.
[32] Han X, Suess E, Huang Y, et al. Jiulong methane reef: microbial mediation of seep-carbonates in the South China Sea[J]. Marine Geology, 2008, 248(3-4): 243-256. doi: 10.1016-j.margeo.2007.11.012/
[33] Tyrrell T. The relative influences of nitrogen and phosphorus on oceanic primary production [J]. Nature, 1999, 400: 525-531. doi: 10.1038/22941
[34] Falkowski P G. Evolution of the nitrogen cycle and its influence on the biological sequestration of CO2 in the ocean [J]. Nature, 1997, 387: 272-275. doi: 10.1038/387272a0
[35] Luo Q Y, Zhong N N, Zhu L, et al. Correlation of burial organic carbon and paleoproductivity in the Mesoproterozoic Hongshuizhuang Formation, northern North China [J]. Chinese Science Bulletin, 2013, 58: 1299-1309. doi: 10.1007/s11434-012-5534-z
[36] Li W, Zhang Z. Paleoenvironment and its control of the formation of oligocene marine source rocks in the deep-water area of the northern South China Sea [J]. Energy & Fuels, 2017, 31(10): 10598-10611. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=338a30a1c81365416408fe4406e675a8
[37] Murray R W, Knowlton C, Leinen M, et al. Export production and carbonate dissolution in the central equatorial Pacific Ocean over the past 1 Myr [J]. Paleoceanography, 2000, 15: 570-592. doi: 10.1029/1999PA000457
[38] Algeo T J, Kuwahara K, Sano H, et al. Spatial variation in sediment fluxes, redox conditions, and productivity in the Permian-Triassic Panthalassic Ocean [J]. Paleogeography, Paleoclimatology, Paleoecology, 2011, 308: 65-83. doi: 10.1016/j.palaeo.2010.07.007
[39] Li Y T, Zhang G S, Ellis, et al. Depositional environment and organic matter accumulation of Upper Ordovician-Lower Silurian marine shale in the Upper Yangtze Platform, South China [J]. Paleogeography, Paleoclimatology, Paleoecology, 2017, 466: 252-264. doi: 10.1016/j.palaeo.2016.11.037
[40] Li D, Li R, Zhu Z, et al. Origin of organic matter and paleo-sedimentary environment reconstruction of the Triassic oil shale in Tongchuan City, southern Ordos Basin (China) [J]. Fuel, 2017, 208: 223-235. doi: 10.1016/j.fuel.2017.07.008
[41] Mclennan. relationship between the trace element composition of sedimentary rocks and upper continental crust [J]. Geochemistry, Geophsics, Geosystem, 2011, 73: 181-187.
[42] Li N, Feng D, Chen L, et al. Compositions of foraminifera-rich turbidite sediments from the Shenhu area on the northern slope of the South China Sea: Implication for the presence of deep water bottom currents [J]. Journal of Asian Earth Sciences, 2017, 138: 148-160. doi: 10.1016/j.jseaes.2017.02.010
[43] Murray R W, Leinen M. Scavenged excess aluminum and its relationship to bulk titanium in biogenic sediment from the central equatorial Pacific Ocean [J]. Geochimica et Cosmochimica Acta, 1996, 60: 3869-3878. doi: 10.1016/0016-7037(96)00236-0
[44] 张辉, 卢海龙, 梁金强, 等.南海北部神狐海域沉积物颗粒对于天然气水合物聚集的主要影响[J].科学通报, 2016, 61(3): 388-397. http://www.cqvip.com/QK/94252X/201603/667908763.html ZHANG Hui, LU Jinghai, LIANG Jingqiang, et al. The main influence of sediment particles in the Shenhu area of the northern South China Sea on gas hydrate accumulation [J]. Science Bulletin, 2016, 61(3): 388-397. http://www.cqvip.com/QK/94252X/201603/667908763.html
[45] Qin J Z, Tenger, Fu X D. Study of formation condition on marine excellent source rocks and its evaluation [J]. Petroleum Geology & Experimrnt, 2009, 31(4): 366-372. http://www.sysydz.net/EN/Y2009/V31/I4/366
[46] Bahk J J, Umik, Holland M. Core lithologies and their constraints on gas-hydrate occurrence in the East Sea, offshore Korea: results from the site UBGH1-9 [J]. Marine and Petroleum Geology, 2 011, 28(10): 1943-1952. doi: 10.1016/j.marpetgeo.2010.12.003
[47] 何家雄, 马文宏, 祝有海, 等.南海北部边缘盆地天然气成因类型及成藏时间综合判识与确定[J].中国地质, 2011, 38(1): 145-160. doi: 10.3969/j.issn.1000-3657.2011.01.016 HE Jiangxiong, MA Wenhong, ZHU Youhai, et al. Comprehensive identification and determination of genetic types and formation time of natural gas in the northern marginal basin of the South China Sea [J]. Chinese Geology, 2011, 38(1): 145-160. doi: 10.3969/j.issn.1000-3657.2011.01.016
[48] 苏丕波, 梁金强, 沙志彬, 等.南海北部神狐海域天然气水合物成藏动力学模拟[J].石油学报, 2011, 32(2): 226-232. http://d.old.wanfangdata.com.cn/Periodical/syxb201102006 SU Pibo, LIANG Jingqiang, SHA Zhibing, et al. Dynamics simulation of gas hydrate formation in Shenhu area, northern South China Sea [J]. Acta Petroleum Sinica, 2011, 32(2): 226-232. http://d.old.wanfangdata.com.cn/Periodical/syxb201102006
[49] 苏丕波, 梁金强, 沙志彬, 等.神狐深水海域天然气水合物成藏的起源条件[J].西南石油大学学报:自然科学版, 2014, 36(2): 1-8. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=xnsyxyxb201402001 SU Pibo, LIANG Jingqiang, SHA Zhibing, et al. Gas sources condition of gas hydrate formation in Shenhu deep water sea zone[J]. Journal of Southwest Petroleum University: Natural Science Edition, 2014, 36(2):: 1-8. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=xnsyxyxb201402001
[50] Sadiq M. Thermodynamic solubility relationships of inorganic vanadium in the marine-environment [J]. Marine Chemistry, 1988, 23: 87-96. doi: 10.1016/0304-4203(88)90024-2
[51] Achterberg E P, Boussemart M, Davison W. Speciation and cycling of trace metals in Esthwaite Water: A productive English lake with seasonal deep-water anoxia [J]. Geochemistry & Geophysics, 1997, 61: 5233-5253. https://www.sciencedirect.com/science/article/pii/S0016703797003165
[52] Jones B, Manning D A. Comparison of geochemical indexes used for the interpretation of paleoredox conditions in ancient mudstones[J]. Chemical Geology. 1994, 111(1-4): 111-129. doi: 10.1016/0009-2541(94)90085-X
[53] Bahk J J, Kim D H, Chun J H, et al. Gas hydrate occurrences and their relation to host sediment properties: results from second Ulleung Basin gas hydrate drilling expedition, East Sea [J]. Marine and Petroleum Geology, 2013. 47(11): 21-29. https://www.sciencedirect.com/science/article/abs/pii/S0264817213001293
[54] Meyer E E, Burgreen B N, Lackey H, et al. Evidence for basin restriction during syn-collisional basin formation in the Silurian Arisaig Group, Nova Scotia [J]. Chemical Geology, 2008, 256(1-2): 1-11. doi: 10.1016/j.chemgeo.2008.06.016