Marine Geology

Deep thermal state in the Nansha Trough of South China Sea and its tectonic implications

  • REN Ziqiang , 1, 3 ,
  • SHI Xiaobin , 1, 2 ,
  • WANG Xiaofang 1, 2 ,
  • ZHAO Peng 1, 3 ,
  • SHEN Yongqiang 1, 3
Expand
  • 1. Key Laboratory of Ocean and Marginal Sea Geology, South China Sea Institute of Oceanology, Innovation Academy of South China Sea Ecology and Environmental Engineering, Chinese Academy of Sciences, Guangzhou 510301, China
  • 2. Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou), Guangzhou 511458, China
  • 3. University of Chinese Academy of Science, Beijing 100049, China
SHI Xiaobin. email:

Copy editor: YAO Yantao

Received date: 2020-07-21

  Revised date: 2020-09-01

  Online published: 2020-09-20

Supported by

National Natural Science Foundation of China(41776078)

National Natural Science Foundation of China(42076075)

Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou)(GML2019ZD0104)

Copyright

Copyright reserved © 2021. Office of Acta Agronomica Sinica All articles published represent the opinions of the authors, and do not reflect the official policy of the Chinese Medical Association or the Editorial Board, unless this is clearly specified.

Abstract

The Nansha Trough is a key region where the Proto-South China Sea subducted southward below the Borneo, and the Nansha Block collided with Borneo. Its tectonic evolution contains important evolutional information of previous and present South China Sea region. To further understand the deep thermal state in the Nansha Trough, which has a great effect on its tectonic deformation, with a numerical model of the thermal diffusion equation, we first analyzed the thermal blanketing effect of rapid accumulation of sediment on observed heat flow in the deep-water region of the Nansha Trough, and then calculated the Curie depths in the study area with a magnetic spectrum analysis technique. The results show that due to the thermal blanketing of rapid accumulation of the Brunei Landslide, the observed seafloor heat flow is just 60%~77% of the deep background heat flow in the deep-water region, which is about 77~98 mW·m-2. The Curie depth in the trough region is generally less than 16 km, which is smaller than the Curie depth (18~24 km) in the Nansha Islands and Reefs. Further analyses show that the trough region is of high background heat flow, and the high thermal state in the study area is in accord with its greatly thinned crustal thickness, which is the result of continental rifting of the southern margin of South China, and the evolution of the South China Sea.

Cite this article

REN Ziqiang , SHI Xiaobin , WANG Xiaofang , ZHAO Peng , SHEN Yongqiang . Deep thermal state in the Nansha Trough of South China Sea and its tectonic implications[J]. Journal of Tropical Oceanography, 2021 , 40(4) : 98 -109 . DOI: 10.11978/2020074

南海南部大陆边缘不仅发育有众多的含油气盆地, 而且在南海构造演化中扮演着重要角色。特别是南沙海槽区域, 古南海在此处向婆罗洲之下俯冲消减, 东南向的俯冲拖曳及中新世海沟两侧地块的碰撞对南海陆缘裂陷和南海形成演化具有重要作用(Taylor et al, 1980), 因而该区域的构造变形史记录着南海前世今生构造演化的重要信息。地质体变形和岩石圈流变性质与其热状态密切相关(马辉 等, 2011; 张健 等, 2017), 因而南沙海槽热状态研究不仅有助于评估沉积盆地的油气资源潜力, 而且可为深入研究南沙海槽区域的流变学强度与变形机制提供地热学基础。
大地热流是用以研究区域深部热状态和评价油气资源潜力的重要参数。南海南部大陆边缘热流分布的一个典型特征是以NW向西巴兰线为界, 东侧的地温梯度和大地热流明显低于西侧(图1)(Hall, 2002; Shi et al, 2003; Cullen, 2010; Hutchison, 2010)。西巴兰线以西的曾母盆地热流平均值高达102mW·m-2, 地温梯度介于36~55mW·m-2(Madon, 1999; Shi et al, 2003); 而西巴兰线以东的文莱-沙巴盆地和南沙海槽底部热流平均值仅约60mW·m-2(图1), 其中文莱-沙巴盆地钻井地温梯度一般低于36mW·m-2 (Shi et al, 2003)。南沙海槽为一北东向槽型深水区, 槽内地形平坦, 水深介于2.8~2.9km之间, 海槽东南侧为自婆罗洲向北持续推覆的文莱-沙巴盆地深水冲断带(图1、2)(Ingram et al, 2004; 唐武 等, 2018a, b)。海底表层热流观测容易受到浅层地质过程如沉积物快速堆积作用、流体活动等因素影响, 海槽东南侧冲断带内地形复杂、流体活动频繁, 热流测量值受到多种因素的影响(Zielinski et al, 2007; Morley et al, 2014), 而海槽底部深水区热流测量结果主要受源自东南侧的沉积物快速堆积的影响。因此, 为获取南沙海槽的深部热状态信息, 其观测热流数据还需进行浅层影响因素分析。
图1 南沙海槽位置及热流值

该图基于中国地质调查局广州海洋地质调查局编制出版的南海晕渲地形图制作, 审图号为JS(2015)02-107, 底图无修改。图中洋红色曲线为3000m等深线, 黄色曲线为南沙海槽区域2000m等深线, 黑色方框为图4a位置, 白色阴影区为文莱滑坡[据Gee等(2007)], 黑色直线为图2图3的剖面位置, 红色细线为文莱沙巴盆地深水冲断带[引自唐武等(2018b)], 红色粗虚线为区域断裂及推测的古俯冲位置[据Cullen (2014)]。图中圆点为热流, 热流数据来源于Shi等(2003)徐行等(2018b), 黑色点为礼乐滩礁体区钻井, 空心圆为陈爱华等(2017)测得的热流站位

Fig. 1 Location of the Nansha Trough, showing observed heat flow sites

图2 南沙地块地壳结构

图a地壳结构据丘学林等(2011), 图b地壳结构据Franke等(2008); 图中红色曲线为居里面深度; 剖面位置见图1

Fig. 2 Crustal structure profiles in the Nansha Block

除热流测量外, 可以利用岩石圈磁性层的分布研究南沙海槽区的深部热状态。岩石圈中磁性矿物的温度随深度增加而升高, 当温度超过550~580℃时, 岩石将失去磁化能力(Bhattacharyya et al, 1975)。因此, 可以通过获取磁性层的底界面即居里面深度, 以反映南沙海槽区域的深部热状态信息。
考虑到南沙海槽底部深水区的热流测量结果主要受沉积物快速堆积的影响, 影响因素较为单一, 为了获取南沙海槽区可靠的热状态信息, 本文选取该区观测热流值, 利用一维数值模拟软件SlugSed (Hutnak et al, 2007)分析沉积物快速堆积对海槽底部深水区海底热流的影响, 获取反映海槽深部热状态的热流值; 并通过分析研究区磁异常数据, 获取南沙海槽及其附近区域的居里面深度; 利用表层热流和磁居里面的共同约束, 明确南沙海槽的深部热状态, 进而探讨其构造意义。

1 地质概况

一般认为, 古近纪期间, 南沙地块位于古南海北侧的被动大陆边缘。随着古南海向婆罗洲之下俯冲消减, 南沙地块与南海现今北部陆缘一起发生裂陷作用, 并在南海海盆扩张过程中逐渐向南漂移(Taylor et al, 1980; Briais et al, 1993; 李家彪 等, 2012)。早中新世末, 南沙地块与婆罗洲发生碰撞作用, 俯冲带南侧的早期俯冲带增生楔等混杂岩呈北西向仰冲于南沙地块的减薄型陆壳之上(图2、3)。中新世以后, 受到婆罗洲南侧印澳板块的北向俯冲影响, 婆罗洲北部经历了多期隆升和剥蚀, 沉积物被运移至沙巴陆缘, 并不断向北推覆在南沙地块之上(Morley, 2007, 2009)。受逆冲推覆体和沉积物的负载作用, 南沙地块南端基底发生快速挠曲下沉(Hall, 2011), 其上发育了文莱-沙巴盆地和南沙海槽盆地(图1)。
图3 南沙海槽深水冲断带地震剖面[修改自Franke等(2008)]

图中黑色线为逆冲断层, 紫色线为冲断带的主滑脱面, 各地层界面年龄据Franke等(2008)。剖面位置见图1

Fig. 3 Seismic profile showing the deep-water thrust zone in the Nansha Trough

南沙海槽位于南沙地块南部, 西南向以西巴兰断裂与曾母盆地相隔, 向东北一直延伸到巴拉巴克断裂, 海槽北侧为南沙岛礁区(包括永暑礁、郑和群礁、礼乐滩和安渡滩等), 南侧为自婆罗洲向北持续推覆的文莱-沙巴盆地深水冲断带(图1、3)(Bell Hamilton, 1979; Hinz et al, 1989; Franke et al, 2008)。冲断带内逆冲断层上盘形成断层传播褶皱, 并突出海底形成背斜脊, 所有逆冲断层向下延伸至同一个滑脱面。冲断带前缘(靠近海槽的一侧)高角度的逆冲断层一直错断至海底(图3)(Clift et al, 2008; Franke et al, 2008), 说明冲断带至今仍在活动。三维地震数据显示, 南沙海槽南端存在可能由背斜顶部崩塌造成的巨型海底滑坡, 称为文莱滑坡(Gee et al, 2007)。该滑坡由婆罗洲陆架前缘延伸至南沙海槽底部(图1、4), 面积超过6300km2(任金锋 等, 2020), 在南沙海槽底部形成约100m厚的沉积(图4b)。根据滑坡上覆沉积层厚度及沉积速率(约6~20cm·ka-1)可估算该滑坡的年龄为2000~7000a(Gee et al, 2007)。同时, 文莱滑坡体之下存在多处透镜状反射体, 说明可能发生过多期滑坡(图4b)。受沉积物不均匀压实的影响, 海槽内局部存在泥火山(Morley et al, 2014)。
图4 文莱滑坡区热流点分布(a)及南沙海槽深水区滑坡体分布(b)[修改自Gee等(2007)]

图a为文莱滑坡区地形, 位置如图1中方框所示; 图中热流点据Zielinski等(2007), 红色点为高热流值点(大于99mW·m-2), 蓝色点为低热流值点(小于40mW·m-2), 绿色点为C182站位, 热流值为27mW·m-2 (Sclater et al, 1976); 黑色长条阴影区为逆冲断层在海底形成的构造脊, 黑色箭头为滑坡滑塌的方向, 黑色曲线为等深线(单位: km), 短虚线为地震识别的滑坡范围; 长虚线为高热流区和低热流区的分界, 其中南侧高热流区平均值为84mW·m-2, 北侧低热流区平均值为59mW·m-2 (Zielinski et al, 2007)。图b为南沙海槽底部穿过文莱滑坡体的剖面[引自Gee等(2007)], 剖面位置如图a所示

Fig. 4 (a) Heat flow stations in the Brunei Landslide and (b) landslide units in the deep-water area of the Nansha Trough

现今南沙地块整体具有减薄的地壳, 其中北部南沙群岛区域岛礁林立、地形复杂, 地壳厚度在10~25km之间(图2), 南部的南沙海槽内水深大于2000m, 地形相对平坦(图1)。穿越南沙海槽的折射地震剖面显示, 地壳厚度自南沙群岛区域向南沙海槽减薄至8~9km(Franke et al, 2008), 南沙海槽基底为强烈减薄的陆壳(Vijayan et al, 2013)。高分辨率地震剖面显示, 海槽内存在高于海底超过1000m的海山, 其上有明显的碳酸盐岩发育特征, 说明海槽区域存在快速沉降(Hutchison, 2010)。

2 热流分布

南沙群岛北侧为南海海盆, 西侧为具有高热流值的曾母盆地(图1)。南沙群岛北部热流平均值约为79mW·m-2(徐行 等, 2018b), 位于永暑礁南侧的ODP1143钻井热流值高达84mW·m-2(Shi et al, 2003)。对南沙海域东部的礼乐滩礁体区的钻井测温数据分析显示, 该区基底热流介于65~75mW·m-2之间(王丽芳 等, 2020; 任自强 等, 2021)。位于南沙海域西部的北康盆地及其邻近海域, 热流值介于43~115mW·m-2之间, 平均值约77mW·m-2(陈爱华 等, 2017)。
根据折射地震解释剖面(图2), 南沙海槽基底为南沙地块东南向延伸的强烈减薄的陆壳。现有热流测量结果显示海槽区域热流值变化很大。Sclater等(1976)曾在该区域获得6个站位的海底热流数据(表1), 其中两个站位(C182和C183)位于海槽两侧, 热流值仅为27mW·m-2和14mW·m-2, 其余4个站位(A176~A179)位于海槽底部, 平均值为61mW·m-2。然而, 其在热流计算过程中所用的沉积物热导率仅为0.70μW·m-3左右(表1), 明显小于南海其他区域的热导率实测值(Xu et al, 2016; 徐行 等, 2018a, b)。如以南海西南次海盆及南沙岛礁区海底沉积物的热导率平均值0.85μW·m-3(徐行 等, 2018b)进行计算, 则海槽底部热流在68~78mW·m-2之间, 平均约74mW·m-2, 与北部南沙群岛区域的热流测量结果相当(表1)。
表1 南沙海槽区域海底热流值

Tab. 1 Heat flow in the Nansha Trough

站位 纬度 经度 水深/m 地温梯度/(℃·km-1) 热导率1
/(μW·m-3)
热导率2
/(μW·m-3)
热流1
/(mW·m-2)
热流2
/(mW·m-2)
A176 6°34′54″N 113°50′30″E 2350 92 0.707 0.85 65 78
A177 6°46′18″N 114°15′18″E 2863 80 0.700 0.85 56 68
A178 6°58′18″N 114°39′42″E 2882 92 0.696 0.85 64 78
A179 7°7′42″N 115°7′E 2844 83 0.699 0.85 58 71
C182 6°N 114°1′E 2305 30 0.909 × 27 ×
C183 7°27′N 114°32′60″E 1476 × × × 14 ×

注: 原始数据来自Sclater等(1976), 其中C183站位曾因刊印错误而被北移至北纬9.45°(Taylor et al, 1983); 1为原始数据, 2为以西南次海盆及南沙岛礁区海底沉积物热导率平均值改正后的数据; ×表示无有效数据

南沙海槽南端文莱滑坡区域海底探针热流测量结果显示(图4a), 海槽东南侧冲断带内的热流值平均约为84mW·m-2, 而靠近南沙海槽的深水区海底热流平均值明显较低, 仅为59mW·m-2 (Zielinski et al, 2007)。海槽底部滑坡区热流测量值明显低于滑坡区以外热流值, 与其强烈减薄的陆壳特征并不一致, 很可能是受到了浅层地质过程的扰动。

3 研究方法与参数

3.1 沉积物热披覆效应影响分析

南沙海槽及其邻近区域的热流站位多集中于南侧的文莱-沙巴盆地, 海槽底部深水区热流站位很少(图1)。文莱-沙巴盆地早期为逆冲推覆堆积, 近期发育重力滑脱体、滑坡体, 区内地质活动复杂, 存在较多影响热流测量的因素, 而南沙海槽底部热流测量值的影响因素主要为滑坡沉积物的快速堆积作用。因此本文选取南沙海槽底部文莱滑坡范围内的热流站位(图4a), 利用SlugSed程序(Hutnak et al, 2007)对滑坡沉积物快速堆积的热披覆效应进行分析, 以获取能够反映深部热背景的海槽区热流值。
SlugSed方法(Hutnak et al, 2007)采用变换有限差分网格计算沉积作用及流体活动对地温场的影响。考虑到本文主要研究沉积物堆积过程对表层热流影响, 因此忽略流体传热及沉积物生热率的影响, 因而热传递方程可简化为:
$\overline{\rho C}\frac{\partial \text{T}}{\partial t}=\frac{\partial }{\partial z}\lambda \frac{\partial \text{T}}{\partial z}-\overline{\rho C\upsilon }\frac{\partial \text{T}}{\partial z}$
式中: ρCT分别为沉积物的密度、比热和温度, t为时间, z为深度, $\lambda $为沉积物热导率, $\upsilon $为沉积速率。$\overline{\rho C}$与$\overline{\rho C\upsilon }$分别表示沉积物骨架及孔隙水的平均性质:
$\overline{\rho C}=\rho {{C}_{\text{w}}}\phi +\rho {{C}_{s}}(1-\phi )$
$\overline{\rho C\nu }=\rho {{C}_{\text{w}}}{{\nu }_{\text{w}}}\phi +\rho {{C}_{s}}{{\nu }_{s}}(1-\phi )$
式中: 下标w和s分别代表孔隙水和沉积物骨架, $\phi $为沉积物孔隙度。深度上相邻两节点z1z2间的速度v由公式(4)和(5)计算获得:
${{v}_{s}}({{z}_{1}},t)\left[ 1-\phi ({{z}_{1}}) \right]={{v}_{s}}({{z}_{2}},t)\left[ 1-\phi ({{z}_{2}}) \right]$
${{v}_{\text{w}}}({{z}_{1}},t)\phi ({{z}_{1}})={{v}_{\text{w}}}({{z}_{2}},t)\phi ({{z}_{2}})$
沉积物的热导率$\lambda $依据沉积物骨架热导率及孔隙水热导率按下式计算:
$\lambda =\lambda _{\text{w}}^{\phi }\lambda _{s}^{\text{1}-\phi }$
式中: λwλs分别为孔隙水和沉积物骨架的热导率。假设沉积物孔隙度仅随深度减小, 并且忽略孔隙水及泥沙骨架的可压缩性时, 孔隙度的计算方程为:
$\phi ={{\phi }_{0}}{{e}^{-cz}}$
式中: ${{\phi }_{0}}$为海底表层沉积物的孔隙度, c为压实系数。
计算模型分为两层, 上层为不断增厚的沉积层, 下层为沉积基底(图5a)。顶边界为定温度边界, 温度取3℃, 底部为定热流边界Q0。计算过程中以模型顶部(海底表层)热流Qb与底边界热流Q0的比值Qb/Q0计量沉积物热披覆作用对海底热流的影响程度。经计算发现, 相同热物理性质及沉积速率下, 底边界热流Q0的取值对Qb/Q0无影响(图5b)。为方便计算, 本研究的底边界热流Q0取100mW·m-2
图5 沉积物热披覆效应计算模型示意图(a)[修改自Hutnak等(2007)]和模型顶、底热流比值随底边界热流值的变化(b)

图a中加号和向下箭头表示沉积层向下增厚。图b中不同直线表示为图a中模型在不同沉积速率下持续沉积100a, 之后沉积作用结束7000a后顶、底边界热流比值随底边界热流值的变化

Fig. 5 (a) Schematic diagram of sediment thermal blanketing model, and (b) ratio of heat flow at the top and bottom of the model versus the bottom boundary conditions

SlugSed方法(Hutnak et al, 2007)可以依据不同基底的热物理性质条件, 计算沉积物热披覆效应对海底热流的影响。考虑到基底热导率对计算结果影响较小, 一般介于2%~4%(Hutnak et al, 2007), 故本文忽略模型基底热导率变化的影响, 其热导率取值2.0W·m-1·K-1。根据地震剖面显示(图4b), 海槽南端的文莱滑坡在海槽底部形成了约100m厚的沉积层, 由于滑坡体滑塌的时间普遍较短, 本文把滑坡体滑塌过程中沉积物在海槽底部的堆积速率近似为恒定常数。如果该100m厚的沉积物分别在1a、10a和100a时间段内堆积完成, 那么对应的平均堆积速率分别为100m·a-1、10m·a-1及1m·a-1。相对于滑坡沉积物的快速堆积, 滑坡滑塌结束后的正常沉积速率很小, 约为6~20cm·ka-1 (Gee et al, 2007), 经计算发现其对滑塌结束后海底热流恢复的影响小于1 %, 因而计算过程中忽略滑坡滑塌结束后的正常沉积过程。模型假定滑坡体沉积物主要为砂质沉积物(Ingram et al, 2004), 模型参数取值见表2
表2 滑坡体热披覆效应影响分析模型的参数与取值

Tab. 2 Parameters and values used in the landslide thermal blanketing model

参数 取值
砂质沉积物骨架
密度/(kg·m-3) 1030* 2650*
比热/(J·kg-1·K-1) 4183* 1088*
热导率/(W·m-1·K-1) 0.67 4.18*
初始孔隙度 × 0.63#
压实因子/m-1 × 0.51×10-3#

注: *引自Hutchison (1985), #引自王丽芳等(2020); ×表示无数据

3.2 居里面深度计算

岩石圈中磁性矿物因温度升高而失去磁性的深度界面被称为居里面, 通常认为磁性地层的底界面为居里温度界面。因此, 如果能获得磁性层的底界, 则可以对研究区深部热状态进行约束。居里面深度可利用区域磁异常的频谱分析获得, 本文采用Salem等(2014)提出的“去分形法”计算南沙海槽及其邻近区域的居里面深度。
依据Salem等(2014)的去分形方法, 当磁性层磁化强度在xy方向上具有分形特征, 而在竖直方向上表现为完全随机分布时, 其磁异常强度功率谱${{\Phi }_{F}}({{k}_{x}},{{k}_{y}})$可以表示为:
${{\Phi }_{F}}({{k}_{x}},{{k}_{y}})=\Phi ({{k}_{x}},{{k}_{y}})\cdot {{k}^{-\alpha }}$
式中: α为该磁性层平面上的分形指数, $\Phi ({{k}_{x}},{{k}_{y}})$为磁化强度完全随机分布时的磁异常强度功率谱, 其计算公式为:
${{\Phi }_{F}}({{k}_{x}},{{k}_{y}})=A\cdot {{\Phi }_{M}}({{k}_{x}},{{k}_{y}}){{({{e}^{-k{{z}_{\text{t}}}}}-{{e}^{-k{{z}_{\text{b}}}}})}^{2}}$
式中: A是与磁化方向和地磁场方向有关的常数, ${{\Phi }_{M}}({{k}_{x}},{{k}_{y}})$为磁化强度M(x, y)的功率谱, ztzb分别为磁性层的顶、底界深度, kxky分别是xy方向上的波数, 并且有:
$k=\sqrt{{{k}_{x}}^{2}+{{k}_{y}}^{2}}$
当磁化强度完全随机分布时, ${{\Phi }_{M}}({{k}_{x}},{{k}_{y}})$为常数(Tanaka et al, 1999), 这时磁异常强度的径向平均功率谱为:
$\Phi (k)={A}'{{e}^{-2k{{z}_{\text{t}}}}}{{\left[ 1-{{e}^{-k({{z}_{\text{b}}}-{{z}_{\text{t}}})}} \right]}^{2}}$
式中: ${A}'$为常数。在中、高波数域, 即当波长小于两倍磁性层厚度时, 公式(11)可近似变为:
$\ln \left[ \Phi (k) \right]=\ln \text{B}-2{{z}_{\text{t}}}k$
式中: B为常数。
假设z0为磁性层的中间点深度, 2d为磁性层的厚度, 即有z0=( z0+ zb)/2。那么在长波长、低波数域有${{e}^{-k({{z}_{\text{t}}}-{{z}_{0}})}}={{e}^{k({{z}_{\text{b}}}-{{z}_{0}})}}\approx k({{z}_{\text{b}}}-{{z}_{0}})$, 因而公式(11)可改写成:
$\begin{align} & \Phi {{(k)}^{\frac{1}{2}}}=C{{e}^{-k{{z}_{0}}}}({{e}^{-k({{z}_{\text{t}}}-{{z}_{0}})}}-{{e}^{-k({{z}_{\text{b}}}-{{z}_{0}})}}) \\ & \ \ \ \ \ \ \ \ \ \approx C{{e}^{-k{{\text{z}}_{0}}}}2k({{z}_{\text{b}}}-{{z}_{0}}) \\ & \ \ \ \ \ \ \ \ \ =C{{e}^{-k{{z}_{0}}}}2kd \\ \end{align}$
对公式(13)两边取对数可得:
$\ln \left\{ \left[ \Phi {{(k)}^{1/2}} \right]/k \right\}=\ln \text{D}-k{{z}_{0}}$
式中: D为常数。因而可以用最小二乘法通过分别拟合$\ln \left[ \Phi {{(k)}^{1/2}} \right]$的中、高波数域和$\ln \left\{ \left[ \Phi {{(k)}^{1/2}} \right]/k \right\}$的低波数域斜率, 以求取磁性层的顶界和中间点的深度, 并通过下式计算磁性层的底界, 即居里面的埋深(Tanaka et al, 1999; Dong et al, 2016):
zb = 2z0-zt
采用去分形法计算时, 通过不断调整分形指数α, 当实测功率谱和理论功率谱一致时, 即可获得磁性层底界深度。具体步骤为: 首先假设一个较小的分形指数α(Salem et al, 2014), 根据式(8)计算实测磁异常的去分形功率谱$\Phi $; 在给定波数域内利用公式(15)计算得到磁性层底界深度zb, 并计算出理论功率谱${{\Phi }_{\text{Mod}}}$; 对比实测功率谱$\Phi $和理论功率谱${{\Phi }_{\text{Mod}}}$, 当$\Phi $与${{\Phi }_{\text{Mod}}}$一致时所得到的zb即为磁性层的底界深度, 如果两者一致性达不到目标值, 则按步长修正分形指数α; 重复前述过程, 直至观测和理论功率谱的一致性达到目标值。
居里面计算所采用的磁异常数据为最新全球磁异常数据EMAG2v3(Meyer et al, 2017), 分辨率为2'(图6a)。为保证计算结果的准确性, 计算窗口避开了区内无磁异常数据的区域。同时为满足窗口尺寸至少大于居里面深度6~10倍的计算要求(Ravat et al, 2007; Salem et al, 2014), 本次计算窗口大小定为256km×256km, 以保证获取足够深的磁性层底界, 且相邻两个窗口之间有50%的重叠, 即窗口在xy方向的移动步长均为128km(图6a)。
图6 南沙及其邻近区域的磁异常分布(a)和计算获得的居里面深度分布(b)

该图基于中国地质调查局广州海洋地质调查局编制出版的南海晕渲地形图制作, 审图号为JS(2015)02-107, 底图无修改。图a中磁异常数据来源于Meyer等(2017); 图b中洋红色曲线为3000m等深线, 蓝色虚线为南沙海槽区2000m等深线, 红色实线为区域断裂, 黑色直线为图2的剖面位置

Fig. 6 (a) Map of total field magnetic anomalies, and (b) Curies depth in Nansha and its adjacent area

4 结果与讨论

4.1 滑坡沉积物堆积对海底热流的影响

图7为3种堆积速率下, 南沙海槽深水区文莱滑坡沉积物对海底热流的影响程度及其随时间的变化曲线。根据该图显示, 在滑坡体发育时, 由于沉积物堆积速率很大, 新堆积的沉积物不能被及时加热, 海底热流在沉积物堆积过程中明显降低; 滑塌结束后, 模型假定没有新的沉积, 滑塌沉积物被缓慢加热, 海底热流逐渐恢复。如图7a所示, 不同堆积速率对海底热流降低过程有较大的影响。当堆积速率较大时, 如100m·a-1和10m·a-1时(相当于滑坡发育持续1a和10a), 由于沉积物快速堆积, 海底热流迅速降低; 而当堆积速率较小时, 如1m·a-1时(相当于滑坡发育持续100a), 海底热流先快速降低然后再缓慢降低, 滑塌结束时海底热流降至最低。由于不同堆积速率下沉积物被加热的时间不同, 因而滑塌结束后海底热流的恢复过程亦存在差别, 然而这种差别仅存在于滑坡形成后的几百年内。如图7a所示, 在滑塌开始后600a内, 海底热流的恢复曲线有明显的差别, 约600a后不同堆积速率对海底热流的影响程度近于相同。
图7 沉积物堆积速率对海底热流的影响

图a为滑坡开始形成后1000a内海底热流恢复程度随时间的变化, 不同曲线分别代表不同的沉积物堆积速率。图b和图c分别为滑坡开始形成后100ka和10ka内海底热流恢复程度随时间的变化, 图b中灰色区域为图c的时间跨度, 图c中灰色区域为估算的文莱滑坡的年龄范围

Fig. 7 Effects of different accumulation rates on heat flow

总体上, 在深部热状态一定并且海底热流仅受文莱滑坡沉积物快速堆积的热披覆作用影响时, 不论滑塌时间持续多长, 南沙海槽深水区海底热流均会在沉积物快速堆积过程中降低, 在约1200a后恢复至深部热流的50%, 而要恢复至90%则至少需要40ka。Gee等(2007)根据文莱滑坡上覆沉积层的厚度估算其年龄为2000~7000a, 由图7b可推测, 现今海槽底部深水区文莱滑坡体上的海底热流仅恢复至深部热流的60%~77%, 说明该滑坡体形成过程中的沉积物快速堆积对现今南沙海槽底部深水区海底热流仍具有较大的影响。

4.2 南沙地块居里面深度变化

图6b为计算得到的南沙及其邻近区域的居里面深度分布图。图中显示, 南沙及其邻近区域的居里面埋深介于12~24km之间, 其中南沙群岛的居里面较深, 南沙海槽以及西巴兰线的居里面相对较浅。
南沙群岛海区的居里面深度在18~24km之间(图2), 其中郑和群礁及礼乐滩等地壳较厚的区域, 其居里面深度最大(大于22km), 位于永暑礁南侧ODP1143井之下的居里面埋深约20km(图6b)。据ODP Leg 184航次报告显示, 1143钻井水深2772m, 沉积物厚度约1250m(Wang et al, 2000), 地震剖面揭示该区域莫霍面埋深约18km, 上、下地壳厚度分别约7km和7.5km(图2)。一维热传导方程计算结果显示, 当顶部温度为3℃, 底部热流为66mW·m-2, 各层岩石的热物性依据表3取值时, 可拟合钻井实测热流值为84mW·m-2, 此时580℃(居里温度)的温度界面埋深约20km(图8a), 与居里面计算结果相当, 说明本次计算得到的居里面深度与该区温度结构具有良好的对应关系。
图8 ODP1143井(a)和南沙海槽区域(b)的一维稳态地温场

一维模型的顶边界为定温度边界, 温度为3℃, 底边界为热流边界。图中红色线为580℃温度线

Fig. 8 One-dimensional steady state temperature field at (a) ODP 1143 drill and (b) in the Nansha Trough

表3 一维地温场计算参数值

Tab. 3 Parameter values used for 1D temperature calculation

参数
沉积层生热率 1.21µW·m-3&
沉积层热导率 1.3W·m-1·K-1†或1.7W·m-1·K-1‡
上地壳生热率 1.8µW·m-3#
上地壳热导率 λ-1=0.33+0.33×10-3T*
下地壳生热率 0.5µW·m-3#
下地壳热导率 λ-1=0.41+0.29×10-3T*
地幔热导率 3.4W·m-1·K-1

注: *Buntebarth(1984), #施小斌等(2020), 徐行等(2017)为ODP1143钻井区域沉积层热导率, 为南沙海槽区沉积层热导率[据施小斌等(2015)]; 因ODP1143钻井区沉积层较薄, 沉积物压实作用较小, 故其沉积层热导率取值也较小。上地壳热导率和下地壳热导率公式中的λ为热导率, T为以℃为单位的地壳温度

由南沙群岛至南沙海槽区域, 居里面抬升至莫霍面以上, 埋深小于16km, 与地壳的减薄趋势相当(图2)。根据地震剖面显示, 南沙海槽内沉积层较厚, 约为5.5km, 上、下地壳厚度分别为4.7km和8.2km(图2)。在一维热传导条件下, 当海底温度取3℃, 底部边界热流取71mW·m-2时, 580℃等值线埋深约为16km(图8b), 此时顶边界热流约为90mW·m-2, 说明南沙海槽稳态热状态下应具有较高的热流值, 但受沉积物快速堆积的热披覆作用影响, 海槽底部滑坡范围内的热流测量值相对较低。
西巴兰线位于南沙海域的西南缘, 可能是一条切穿地壳的NW向走滑断裂带(熊莉娟 等, 2012; 赵斐宇 等, 2017), 沿断裂带居里面深度在16km以浅。该断裂带是古南海向南俯冲的西侧界限, 断裂两侧地壳流变性质不同(Cullen, 2014)。沿断裂带热流较高(图1), 与该线较浅的居里面深度相符。

4.3 南沙海槽深部热状态

上述计算结果表明, 南沙海槽底部深水区文莱滑坡范围内较低的热流值明显受到了滑坡沉积物快速堆积热披覆作用的影响。如以区内热流观测平均值计算(图4a), 该区深部热流介于77~98mW·m-2之间, 表现为较高的热背景。如图4b所示, 除文莱滑坡以外, 海槽内还发育有多期的古滑坡, 这些更早期的古滑坡活动对现今海底热流可能也有一定的热披覆效应, 说明海槽底部深水区可能具有更高的深部背景热流。依据南沙海槽区较浅的居里面埋深, 当没有沉积物堆积的热披覆效应影响时, 海槽区的海底热流值可能达到90mW·m-2左右(图8b), 与文莱滑坡区沉积物热披覆效应分析得到的深部背景热流值相当。据图8中的计算结果显示, 南沙海槽现今热流大部分来源于深部地幔(分别为66mW·m-2和71mW·m-2), 也说明南沙海槽区深部具有较高的热背景。
南沙海槽区的高背景热流与其新生代构造演化过程相符。现今南海海底扩张之前, 其南侧的古南海于始新世开始向南俯冲于西北婆罗洲之下(Hall et al, 2017), 南沙地块与现今南海北部陆缘一起经历了强烈的拉张减薄作用。在岩石圈拉张减薄过程中, 深部地幔上涌, 南沙地块热流及深部温度持续升高(赵长煜 等, 2014; Tang et al, 2014; 施小斌 等, 2020)。渐新世至早中新世期间, 南沙岛礁区随着古南海的拖曳作用和南海海盆的扩张逐渐向南漂移, 并于早中新世末与婆罗洲碰撞(Hall et al, 2017)。地震剖面显示, 南沙海槽地壳已经被强烈减薄, 现今海槽区及南侧冲断带之下的地壳厚度仅约8~9km (图2)。因而南沙海槽内较高的深部热状态与其地壳的减薄特征相符, 是新生代华南陆缘裂解及南海形成过程中地壳裂陷作用的结果, 海槽底部滑坡范围内现今较低的热流观测值主要受来自东南侧沉积物的热披覆作用影响。
现今南沙海槽是在强烈减薄陆壳基础上, 受沉积物(尤其是东南侧推覆体)负载作用而挠曲下沉形成的前陆盆地(韩冰 等, 2015; 张健 等, 2017), 其下为中中新世之前发育的被动大陆边缘裂陷盆地。挠曲盆地的宽度和挠曲深度不仅与负载大小有关, 而且与其岩石圈强度密切相关。在负载分布一定的情况下, 岩石圈强度越小, 盆地宽度越小, 盆地则越深。由于岩石圈强度与温度密切相关, 因此具有高背景热流的区域往往具有较低的岩石圈强度(马辉 等, 2011; 陈波, 2013; 张健 等, 2017)。依据重力异常和地形数据谱分析获得的南沙海槽区有效弹性厚度仅约15km (Shi et al, 2017), 这亦表明研究区具有较低的岩石圈强度和较高的背景热流。因此, 现今南沙海槽的宽度、基底埋深以及前隆位置与其较低的岩石圈强度和深部高热背景密切相关。

5 结论

本文利用数值模拟方法分析了文莱滑坡沉积物的堆积对南沙海槽深水区海底热流值的影响, 通过磁异常数据分析获取了研究区的居里面深度, 并在此基础上明确了南沙海槽的深部热状态。主要结论如下:
1) 受沉积物快速堆积的影响, 南沙海槽底部文莱滑坡范围内的海底热流在滑坡发育时明显降低, 在滑坡体活动停止后又逐渐升高, 其现今海底热流仅恢复至深部热流的60%~77%。依据该区的热流测量平均值(59mW·m-2)估算, 海槽深部热流约为77~98mW·m-2
2) 南沙及其邻近区域的居里面埋深介于12~ 24km之间, 其中南沙群岛内居里面深度约18~24km, 南沙海槽内居里面深度相对较浅, 在16km以上。南沙地块内居里面深度与其温度结构有良好的对应关系。
3) 依据南沙海槽居里面深度和实测海底热流的分析结果, 可推测南沙海槽深部具有较高的背景热流。该区较高的深部热状态与其超薄的地壳厚度对应, 是新生代以来区域构造演化过程中地壳裂陷减薄作用的结果。
[1]
陈爱华, 徐行, 罗贤虎, 等, 2017. 南海北康盆地热流分布特征及其构造控制因素探讨[J]. 地质学报, 91(8):1720-1728.

CHEN AIHUA, XU XING, LUO XIANHU, et al, 2017. Heat flow characteristics and controlling factors of the Baikang Basin in South China Sea[J]. Acta Geologica Sinica, 91(8):1720-1728 (in Chinese with English abstract).

[2]
陈波, 2013. 中国及邻区岩石圈有效弹性厚度及其动力学意义[D]. 武汉: 中国地质大学(武汉).

CHEN BO, 2013. The effective elastic thickness over China and surroundings and its lithosphere dynamic implication[D]. Wuhan: China University of Geosciences (Wuhan) (in Chinese with English abstract).

[3]
韩冰, 朱本铎, 万玲, 等, 2015. 南沙海槽东南缘深水逆冲推覆构造[J]. 地质论评, 61(5):1061-1067.

HAN BING, ZHU BENDUO, WAN LING, et al, 2015. Deep-water Fold and Thrust Tectonics in Southeastern Nansha Trough[J]. Geological Review, 61(5):1061-1067 (in Chinese with English abstract).

[4]
李家彪, 丁巍伟, 吴自银, 等, 2012. 南海西南海盆的渐进式扩张[J]. 科学通报, 57(20):1896-1905.

LI JIABIAO, DING WEIWEI, WU ZIYIN, et al, 2012. The propagation of seafloor spreading in the southwestern subbasin, South China Sea[J]. Chinese Science Bulletin, 57(24):3182-3191.

[5]
马辉, 许鹤华, 施小斌, 等, 2011. 南沙海槽前陆盆地热流变结构[J]. 地球科学(中国地质大学学报), 36(5):939-948.

MA HUI, XU HEHUA, SHI XIAOBIN, et al, 2011. Lithospheric thermal-rheological structure of Nansha Trough Foreland basin in South China Sea[J]. Earth Science-Journal of China University of Geosciences, 36(5):939-948 (in Chinese with English abstract).

[6]
丘学林, 赵明辉, 敖威, 等, 2011. 南海西南次海盆与南沙地块的OBS探测和地壳结构[J]. 地球物理学报, 54(12):3117-3128.

QIU XUELIN, ZHAO MINGHUI, AO WEI, et al, 2011. OBS survey and crustal structure of the Southwest Sub-basin and Nansha block, South China sea[J]. Chinese Journal of Geophysics, 54(12):3117-3128 (in Chinese with English abstract).

[7]
任金锋, 孙鸣, 韩冰, 2011. 南海南沙海槽大型海底滑坡的发育特征及成因机制[J]. 地球科学, 46(3):1058-1071.

REN JINFENG, SUN MING, HAN BING, 2011. A Giant Submarine Landslide and Its Triggering Mechanisms on the Nansha Trough Margin, South China Sea[J]. Earth Science-Journal of China University of Geosciences, 46(3):1058-1071 (in Chinese with English abstract).

[8]
任自强, 施小斌, 杨小秋, 等, 2021. 南海礼乐盆地礁体水热循环及其对地温场的影响[J]. 地球物理学报, 64(2):612-627.

REN ZIQAING, SHI XIAOBIN, YANG XIAOQIU, et al, 2021. Hydrothermal circulation and its effect on geothermal field in the reef developing area of the Liyue Basin, South China Sea[J]. Chinese Journal of Geophysics, 64(2):612-627 (in Chinese with English abstract).

[9]
施小斌, 王振峰, 蒋海燕, 等, 2015. 张裂型盆地地热参数的垂向变化与琼东南盆地热流分布特征[J]. 地球物理学报, 58(3):939-952.

SHI XIAOBIN, WANG ZHENFENG, JIANG HAIYAN, et al, 2015. Vertical variations of geothermal parameters in rifted basins and heat flow distribution features of the Qiongdongnan Basin[J]. Chinese Journal of Geophysics, 58(3):939-952 (in Chinese with English abstract).

[10]
施小斌, 王丽芳, 任自强, 等, 2020. 南海礼乐盆地新生代构造热演化特征及其影响因素[J]. 地球物理学报, 63(7):2682-2696.

SHI XIAOBIN, WANG LIFANG, REN ZIQIANG, et al, 2020. Cenozoic Tectono-thermal evolution features and influence factors of the Liyue Basin, South China Sea[J]. Chinese Journal of Geophysics, 63(7):2682-2696 (in Chinese with English abstract).

[11]
唐武, 赵志刚, 张功成, 等, 2018a. 文莱—沙巴盆地深水褶皱冲断带构造变形特征及成因机制[J]. 地球物理学报, 61(10):4281-4295.

TANG WU, ZHAO ZHIGANG, ZHANG GONGCHENG, et al, 2018a. Structural deformation characteristics and genetic mechanism of deepwater fold and thrust belts in the Brunei-Sabah Basin[J]. Chinese Journal of Geophysics, 61(10):4281-4295 (in Chinese with English abstract).

[12]
唐武, 赵志刚, 谢晓军, 等, 2018b. 西北婆罗洲深水褶皱冲断带相关褶皱构造样式[J]. 地球科学, 43(2):491-501.

TANG WU, ZHAO ZHIGANG, XIE XIAOJUN, et al, 2018b. Structural patterns of thrust-related folds of Deepwater fold and thrust belts in NW Borneo[J]. Earth Science, 43(2):491-501 (in Chinese with English abstract).

[13]
王丽芳, 施小斌, 任自强, 等, 2020. 南海南部礼乐盆地礁体发育区的构造热演化特征[J]. 地球物理学报, 63(8):3050-3062.

WANG LIFANG, SHI XIAOBIN, REN ZIQIANG, et al, 2020. Tectono-thermal evolution features of the reef body developing area in the Liyue Basin, southern South China Sea[J]. Chinese Journal of Geophysics, 63(8):3050-3062 (in Chinese with English abstract).

[14]
熊莉娟, 李三忠, 索艳慧, 等, 2012. 南海南部新生代控盆断裂特征及盆地群成因[J]. 海洋地质与第四纪地质, 32(6):113-127.

XIONG LIJUAN, LI SANZHONG, SUO YANHUI, et al, 2012. Cenozoic basin-controlling faults and their bearing on basin groups formation in the southern South China Sea[J]. Marine Geology & Quaternary Geology, 32(6):113-127 (in Chinese with English abstract).

[15]
徐行, 王先庆, 彭登, 等, 2018a. 南海西北次海盆及其邻区的地热流特征与研究[J]. 地球科学, 43(10):3391-3398.

XU XING, WANG XIANQING, PENG DENG, et al, 2018a. Characteristics and research of heat flow in the Northwest Sub-Basin and its adjacent areas of the South China Sea[J]. Earth Science, 43(10):3391-3398. (in Chinese with English abstract).

[16]
徐行, 姚永坚, 彭登, 等, 2018b. 南海西南次海盆的地热流特征与分析[J]. 地球物理学报, 61(7):2915-2925.

XU XING, YAO YONGJIAN, PENG DENG, et al, 2018b. The characteristics and analysis of heat flow in the Southwest sub-basin of South China Sea[J]. Chinese Journal of Geophysics, 61(7):2915-2925 (in Chinese with English abstract).

[17]
张健, 董淼, 吴时国, 等, 2017. 南沙海槽岩石圈热-流变结构与动力学演化分析[J]. 地学前缘, 24(3):27-40.

ZHANG JIAN, DONG MIAO, WU SHIGUO, et al, 2017. Lithosphere thermal-rheological structure and geodynamic evolution model of the Nansha trough basin, South China Sea[J]. Earth Science Frontiers, 24(3):27-40 (in Chinese with English abstract).

[18]
赵长煜, 宋海斌, 杨振武, 等, 2014. 南海南部边缘沉积盆地构造-热演化历史[J]. 地球物理学报, 57(5):1543-1553.

ZHAO CHANGYU, SONG HAIBIN, YANG ZHENWU, et al, 2014. Tectonic and thermal evolution modeling for the marginal basins of the Southern South China Sea[J]. Chinese Journal of Geophysics, 57(5):1543-1553 (in Chinese with English abstract).

[19]
赵斐宇, 姜素华, 李三忠, 等, 2017. 曾母拉分盆地的油气成藏规律: 廷贾-西巴兰线断裂的控制[J]. 海洋地质与第四纪地质, 37(4):209-220.

ZHAO FEIYU, JIANG SUHUA, LI SANZHONG, et al, 2017. Hydrocarbon accumulation in Zengmu pull-apart basin: insights from the control of Tin jar-West Baram line[J]. Marine Geology & Quaternary Geology, 37(4):209-220 (in Chinese with English abstract).

[20]
BELL HAMILTON W, 1979. Tectonics of the Indonesian region[M]. Washington D C: United States Government Printing Office.

[21]
BHATTACHARYYA B K, LEU L -K, 1975. Spectral analysis of gravity and magnetic anomalies due to two-dimensional structures[J]. Geophysics, 40(6):993-1013.

DOI

[22]
BRIAIS A, PATRIAT P, TAPPONNIER P, 1993. Updated interpretation of magnetic anomalies and seafloor spreading stages in the South China Sea: implications for the tertiary tectonics of Southeast Asia[J]. Journal of Geophysical Research: Solid Earth, 98(B4):6299-6328.

[23]
BUNTEBARTH G, 1984. Geothermics: An Introduction[M]. Berlin: Springer.

[24]
CLIFT P, LEE G H, DUC N A, et al, 2008. Seismic reflection evidence for a Dangerous Grounds miniplate: No extrusion origin for the South China Sea[J]. Tectonics, 27(3): TC3008.

[25]
CULLEN A, 2014. Nature and significance of the West Baram and Tinjar Lines, NW Borneo[J]. Marine and Petroleum Geology, 51:197-209.

DOI

[26]
CULLEN A B, 2010. Transverse segmentation of the Baram-Balabac Basin, NW Borneo: refining the model of Borneo's tectonic evolution[J]. Petroleum Geoscience, 16(1):3-29.

DOI

[27]
DONG MIAO, WU SHIGUO, ZHANG JIAN, 2016. Thinned crustal structure and tectonic boundary of the Nansha Block, southern South China Sea[J]. Marine Geophysical Research, 37(4):281-296.

DOI

[28]
FRANKE D, BARCKHAUSEN U, HEYDE I, et al, 2008. Seismic images of a collision zone offshore NW Sabah/Borneo[J]. Marine and Petroleum Geology, 25(7):606-624.

DOI

[29]
GEE M J R, UY H S, WARREN J, et al, 2007. The Brunei slide: A giant submarine landslide on the North West Borneo Margin revealed by 3D seismic data[J]. Marine Geology, 246(1):9-23.

DOI

[30]
HALL R, 2002. Cenozoic geological and plate tectonic evolution of SE Asia and the SW Pacific: computer-based reconstructions, model and animations[J]. Journal of Asian Earth Sciences, 20(4):353-431.

DOI

[31]
HALL R, 2011. Australia-SE Asia collision: plate tectonics and crustal flow[J]. Geological Society, London, Special Publications, 355(1):75-109.

DOI

[32]
HALL R, BREITFELD H T, 2017. Nature and demise of the Proto-South China Sea[J]. Bulletin of the Geological Society of Malaysia, 63:61-67.

DOI

[33]
HINZ K, FRITSCH J, KEMPTER E H K, et al, 1989. Thrust tectonics along the north-western continental margin of Sabah/Borneo[J]. Geologische Rundschau, 78(3):705-730.

DOI

[34]
HUTCHISON C S, 2010. The north-west Borneo trough[J]. Marine Geology, 271(1-2):32-43.

DOI

[35]
HUTCHISON I, 1985. The effects of sedimentation and compaction on oceanic heat flow[J]. Geophysical Journal of the Royal Astronomical Society, 82(3):439-459.

[36]
HUTNAK M, FISHER A T, 2007. Influence of sedimentation, local and regional hydrothermal circulation, and thermal rebound on measurements of seafloor heat flux[J]. Journal of Geophysical Research: Solid Earth, 112(B12):B12101.

[37]
INGRAM G M, CHISHOLM T J, GRANT C J, et al, 2004. Deepwater North West Borneo: hydrocarbon accumulation in an active fold and thrust belt[J]. Marine and Petroleum Geology, 21(7):879-887.

DOI

[38]
MADON M, 1999. The petroleum geology and resources of Malaysia[M]. Kuala Lumpur: Petronas.

[39]
MEYER B, CHULLIAT A, SALTUS R, 2017. Derivation and error analysis of the earth magnetic anomaly grid at 2 arc min Resolution Version 3 (EMAG2v3)[J]. Geochemistry, Geophysics, Geosystems, 18(12):4522-4537.

DOI

[40]
MORLEY C K, 2007. Interaction between critical wedge geometry and sediment supply in a deep-water fold belt[J]. Geology, 35(2):139-142.

DOI

[41]
MORLEY C K, 2009. Growth of folds in a deep-water setting[J]. Geosphere, 5(2):59-89.

DOI

[42]
MORLEY C K, WARREN J, TINGAY M, et al, 2014. Comparison of modern fluid distribution, pressure and flow in sediments associated with anticlines growing in deepwater (Brunei) and continental environments (Iran)[J]. Marine and Petroleum Geology, 51:210-229.

DOI

[43]
RAVAT D, PIGNATELLI A, NICOLOSI I, et al, 2007. A study of spectral methods of estimating the depth to the bottom of magnetic sources from near-surface magnetic anomaly data[J]. Geophysical Journal International, 169(2):421-434.

DOI

[44]
SALEM A, GREEN C, RAVAT D, et al, 2014. Depth to Curie temperature across the central Red Sea from magnetic data using the de-fractal method[J]. Tectonophysics, 624- 625:75-86.

[45]
SCLATER J G, KARIG D, LAWVER L A, et al, 1976. Heat flow, depth, and crustal thickness of the marginal basins of the South Philippine Sea[J]. Journal of Geophysical Research, 81(2):309-318.

DOI

[46]
SHI XIAOBIN, QIU XUELIN, XIA KANYUAN, et al, 2003. Characteristics of surface heat flow in the South China Sea[J]. Journal of Asian Earth Sciences, 22(3):265-277.

DOI

[47]
TANAKA A, OKUBO Y, MATSUBAYASHI O, 1999. Curie point depth based on spectrum analysis of the magnetic anomaly data in East and Southeast Asia[J]. Tectonophysics, 306(3-4):461-470.

DOI

[48]
TANG XIAOYIN, CHEN LIN, HU SHENGBIAO, et al, 2014. Tectono-thermal evolution of the Reed Bank Basin, Southern South China Sea[J]. Journal of Asian Earth Sciences, 96:344-352.

DOI

[49]
TAYLOR B, HAYES D E, 1980. The tectonic evolution of the South China Basin[M]// HAYES D E. The tectonic and geologic evolution of Southeast Asian Seas and Islands, volume 23. Washington, D. C.: American Geophysical Union.

[50]
TAYLOR B, HAYES D E, 1983. Origin and history of the South China Sea basin[M]// Hayes D E. The tectonic and geologic evolution of Southeast Asian Seas and Islands: part 2, volume 27. Washington, D. C.: American Geophysical Union.

[51]
VIJAYAN V R, FOSS C, STAGG H, 2013. Crustal character and thickness over the Dangerous Grounds and beneath the Northwest Borneo Trough[J]. Journal of Asian Earth Sciences, 76:389-398.

DOI

[52]
WANG PINXIAN, PRELL W L, BLUM P, 2000. Proceedings of the ocean drilling program, initial reports (leg 184) south China sea[R]. Texas A&M University, College Station.

[53]
XU XING, WANG ZHIFENG, WANG XIANQING, et al, 2016. The characteristics of heat flow in the Shenhu gas hydrate drilling area, northern South China Sea[J]. Marine Geophysical Research, 37(4):325-335.

DOI

[54]
ZIELINSKI G W, BJORØY M, ZIELINSKI R L B, et al, 2007. Heat flow and surface hydrocarbons on the Brunei continental margin[J]. AAPG Bulletin, 91(7):1053-1080.

DOI

Outlines

/