王蓓羽, 张健, 艾依飞
中国科学院大学地球与行星科学学院 中国科学院计算地球动力学重点实验室, 北京 100049
2018年11月15日 收稿; 2019年5月20日 收修改稿
基金项目: 中国科学院战略性先导科技专项（XDB42020104）、国家自然科学基金重点项目（91858212）和国家自然科学基金（41906056，U1701245）资助
通信作者: 张健，E-mail:zhangjian@ucas.ac.cn

摘要: 南海盆地位处特提斯、环太平洋和古亚洲三大构造域交汇部，其西南部代表性万安盆地新生代地层发育齐全，但前人对其构造沉降及伸展机制研究比较薄弱，制约了人们对南海扩张和盆地演化的理解。基于构造格局和地层层序等区域地质资料，利用回剥法模拟万安盆地的构造演化史，并通过McKenzie纯剪切拉伸模型反复调整拉张因子拟合构造沉降史，得到万安盆地的伸展因子与伸展速度，最终通过岩石圈的热传导方程求取基底热流。研究结果表明：南海扩张、洋脊跳跃和西南次海盆扩张的共同作用控制万安盆地的沉降、伸展和热异常；万安盆地具有独特的时空演化特征，即沉降作用受断裂带左行和右行活动的影响，整体上具有V型转折特性，而伸展作用以“早东西、后南北”为主要特点；伸展和热异常均受构造活动的影响，呈现出以断裂带为异常轴、向两侧以不同的速率逐渐减弱的特征。
关键词: 万安盆地沉降速率伸展速度基底热流
Analysis of tectonic evolution characteristics of Wan'an Basin in the southwest of the South China Sea
WANG Beiyu, ZHANG Jian, AI Yifei
Key Laboratory of Computational Geodynamics of CAS, College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, Beijing 100049, China


Abstract: The South China Sea Basin is located at the junction of Tethys, Pacific margin, and paleo-Asia tectonic domains. The representative Wan'an Basin in the southwest of the South China Sea Basin has the well-developed Cenozoic strata, but the underappreciated expansion and evolution remains because of a limited understanding of tectonic subsidence and extension mechanism of the basins. Based on the regional geological data, including structural style and stratigraphic sequence, the tectonic evolution history of Wan'an Basin is simulated by utilizing the backstripping method, and the extension factor and extension velocity of Wan'an Basin are obtained after fitting the structural subsidence history through an iterative adjustment of the tensile factor in the McKenzie pure shear tensile model. Furthermore, the heat transfer equation of the lithosphere is used to determine the heat flow of the substrate. The results are given as follows. 1) It is the combined action of the South China Sea expansion, the ridge jump, and the southwest sub-basin expansion that control the sedimentation, extension, and thermal anomalies in Wan'an Basin. 2) Wan'an Basin has unique characteristics of space-time evolution that the left and right rotations of the fracture zone affect the sedimentation which presents a V-shaped transition and that "early east and west, then north and south" is the mainly feature of the extension. 3) Due to the tectonic activity, the extension and thermal anomalies in the basin gradually decline in the directions from the fault to the both sides at different rates.
Keywords: Wan'an Basinsedimentation ratestretch ratebasement heat flow
万安盆地位于南海西南部(图 1(a))，是南沙海区重要的新生代含油气盆地之一。该区地处特提斯、环太平洋和古亚洲3大构造域的交汇区，具有极其独特的大地构造，同时受欧亚、太平洋和印—澳3大板块相互作用的控制，形成了复杂的地质构造，具有丰富的油气资源，是研究和认识南海扩张、资源聚集规律和环境变迁的天然实验室^[1-4]。万安盆地的沉降与伸展均和南海扩张密切相关^[5]，共经历了晚始新世—早渐新世、上新世—第四纪、早中新世3次快速的构造沉降过程^[6]。利用地震剖面对构造沉降量、构造沉积速率与拉伸因子的反演结果，表明万安盆地具有差异性构造沉降的特点^{[5, 7]}。渐新世表现为北快南慢，具有北部坳陷和中部坳陷2个沉降中心；早中新世存在北部坳陷、中部坳陷和东南坳陷3个沉降中心，其中，中部坳陷、北部坳陷沉降中心先后西移，东南坳陷则向东迁移；上新世万安盆地的沉降表现出从西向东和从南向北逐渐增强的特征^{[5, 7]}。由以上可见，前人对万安盆地构造沉降和伸展的研究主要集中在沉降量大小和沉降中心所处位置两个方面^{[5, 8-9]}，但对于相关断裂带性质对万安盆地沉降、伸展、基底热流影响的理解相对薄弱。因此，本文利用回剥法进行相关研究，以期深化对万安盆地构造控制与形成演化的认识。
Fig. 1

	Download: JPG larger image
图 1 万安盆地地质构造图 Fig. 1 Geological structure of Wan'an Basin 图 1 万安盆地地质构造图 Fig. 1 Geological structure of Wan'an Basin -->

1 地质构造1.1 构造特点万安盆地位于南海西南部，东经107°50′~110°20′，北纬5°40′~9°30′，面积约8.5×10⁴km²，呈纺锤状、近南北向展布，主体水深小于500m，最大水深约1800~2000m，属于新生代走滑拉张盆地，如图 1(b)(图 1中盆地边界、断层位置、构造分区来自前人研究结果^{[5, 8, 10-11]})。沉积盖层为上始新统—第四系，最大厚度达12500m^{[6, 12]}，其构造格局复杂，形成“四坳四隆一阶一斜”的构造格局，自北向南分别为西北断阶带、北部坳陷、北部隆起、中部坳陷、中部隆起、西部隆起、西南斜坡、南部坳陷、东部隆起和东部坳陷10个二级构造单元(图 1(b))^[5]。万安盆地内断裂构造极为发育，以张扭性NE-SW向切割基底的正断层为主，具有多期性、多层次性的特点。根据前人的研究成果，本文在图 1(b)中标出了万安盆地边界内的一级断裂和二级断裂。其中一级断裂与盆地区域构造走向一致，二级断裂控制着盆地二级构造格局的形成与地层的发育。万安断裂位于万安盆地东部边界，呈北北东走向延伸，是一条基底断裂，形成于中生代，并具有多期次的活动和左—右行走滑交替变化的特征^[13-16]，第四纪以来，表现为弱右行走滑的特征。其他各条断裂大多表现为张性正断层性质，少见走滑断层和反转断层^{[5, 8, 10-11]}。
1.2 地层划分根据地震反射与钻井资料，万安盆地自下而上可识别出Tg、T5、T4、T31、T3、T2和T1共7个区域不整合面，分别对应下始新统人骏群(65Ma)、下始新统—渐新统西卫群(37.8Ma)、下中新统万安组(23.03Ma)、中中新统李准组(15.79Ma)、上中新统昆仑组(11.63Ma)、上新统广雅组(5.33Ma)、第四系(2.58Ma)底界面。但不整合面T5和Tg之间的人骏群只在极少处深凹部位可见，大多数地区处于缺失状态^{[5, 8]}，见图 2(a)、2(b)。人骏群形成于初始裂谷期，受礼乐运动的影响，沉积相为河流湖泊沉积相，岩性为粉砂质泥岩和砂砾岩。西卫群主要形成于裂谷期，受西卫运动和南海运动的影响，沉积相为滨海沼泽相、海湾相，岩性为页岩、砂岩、煤层；万安组主要形成于裂谷期，受洋脊跃迁的影响，沉积相为三角洲、滨浅海、碳酸盐岩台地相，岩性为页岩、砂岩、灰岩；李准组形成于走滑改造期，受南沙运动的影响，沉积相为三角洲、滨浅海、碳酸盐岩台地相，岩性为页岩、砂岩、灰岩；昆仑组形成于走滑改造期，受万安运动的影响，沉积相为滨浅海-半深海相，岩性为页岩、砂岩、灰岩；广雅组与第四系形成于裂后加速沉降期，受广雅运动的影响，沉积相为滨浅海-半深海相，岩性为页岩、砂岩^[9]。
Fig. 2

	Download: JPG larger image
图 2 地质剖面图(剖面位置参见图 1) Fig. 2 Geological sections (see Fig. 1 for the location of the sections) 图 2 地质剖面图(剖面位置参见图 1) Fig. 2 Geological sections (see Fig. 1 for the location of the sections) -->

2 数据与方法2.1 地质与地震数据为研究断层活动对万安盆地沉积特征的影响，本文选取穿插于断裂的两条经广海局人员解译后的人工地震剖面，通过时深转换获得两条剖面A-B、C-D的地层层位，如图 2(位置见图 1(b))。图 2横坐标为距离剖面NW向起点的距离，纵坐标为通过时深转换得到的地层深度，方块处的数字表示利用各地层的岩性数据计算得到的各地层的岩性密度，具体参数与计算方法在下文研究方法里详细列出，圆点处的数字表示通过地层厚度与双程反射时间的关系求得的各层平均地震波速。从图 1中可看出剖面A-B穿过中部隆起、南部坳陷、东部隆起3个区域，其中间标出的3条断裂为图 1(b)中穿过断层的3条断裂；剖面C-D穿过中部坳陷和中部隆起2个区域，其左边第2条断裂为图 1中最左侧穿过剖面C-D的断裂。从图 2中可看到两条地质剖面均具有6个地层层位，分别为西卫、万安、李准、昆仑、广雅、第四系；沉积层厚度在2000~9000m，断层处厚度异常增大；不同层位岩石密度随深度增加而增加，与压实程度正相关，数值位于1.65~2.54g/cm³之间；不同层位的地震波速具有明显差距，与密度深度正相关，数值位于2.03~3.53km/s之间(图 2)。
2.2 研究方法将地质剖面数据通过逐层剥去，脱压实，利用迭代求解逐层回剥后的深度，得到地质剖面的构造演化图，再除去沉积物负荷沉降、古水深以及海平面变化等对各沉积时期沉积厚度的影响，最终得到构造沉降量^[17-18]。经过去压实、古水深和古海平面变化等校正后的Airy均衡构造沉降量S_t表达式^[19-20]为

${S_t} = {S^*}\frac{{\rho _{\rm{s}}^* - {\rho _{\rm{w}}}}}{{{\rho _{\rm{a}}} - {\rho _{\rm{w}}}}} - {W_{\rm{d}}}\frac{{{\rho _{\rm{w}}}}}{{{\rho _{\rm{a}}} - {\rho _{\rm{w}}}}} + {W_{\rm{d}}} \pm \Delta {S_{\rm{L}}}.$

(1)

式中:S^*是经压实校正后沉积柱的长度(m)；ρ_s^*为沉积物的平均密度(g/cm³)；S_L为古海平面(m)；W_d为古水深(m)；ΔS_L为海平面的变化。由于缺少钻井和3D地震资料，无法利用沉积相分析万安盆地古水深变化，我们采用由现今的水深和南海张裂初始时的水深进行线性插值的方法得到古水深^[21-22]。其中南海张裂初期的水深为经古海平面变化校正后的全球海平面变化曲线^[23]。ρ_a、ρ_w分别为地幔、水的密度，取3.33、1.00g/cm³；ρ_s^*在压实作用下随地层厚度而变化^[24-25]，可根据下式计算得出：

$\rho _{\rm{s}}^* = \frac{{\sum\limits_{j = 1}^i {\left[ {{\phi _j}{\rho _{\rm{w}}} + \left( {1 - {\phi _j}} \right){\rho _{\rm{b}}}} \right]S_j^*} }}{{{S^*}}}.$

(2)

式中:ρ_b为骨架岩性的密度；S_j^*为经压实校正后的地层厚度；?_j为地层j的孔隙度，服从孔隙度?随深度指数分布的特征^[24](式(3))，该区地层为岩性互层，所以可利用混合比例加权的方式求解

${\phi _j} = {\phi _0}{e^{ - cz}}.$

(3)

式中:?₀为岩石的初始孔隙度；z为深度；c为压实系数。本次计算，各层岩性的参数主要参考研究区20口井的钻井地层岩性的分析^[5]，见表 1，沉积物的孔隙度参考前人在北海区域的实测结果^[26-28]，见表 2。
Table 1

表 1 万安盆地地层岩性[5]Table 1 Formation lithology of Wan'an Basin[5] ? % 页岩 砂岩 灰岩 煤层 第四系 84.68 15.32 广雅组 70.65 29.35 昆仑组 54.23 40.21 5.56 李准组 30.67 55.45 13.88 万安组 39.78 56.82 2.95 0.45 西卫组 32.38 66.12 1.5

表 1 万安盆地地层岩性[5]Table 1 Formation lithology of Wan'an Basin[5] ?

Table 2

表 2 万安盆地岩性物理参数[5]Table 2 Lithological physical parameters of Wan'an Basin[5] 泥岩 砂岩 灰岩 煤层 初始孔隙度 0.63 0.52 0.65 0.92 压实系数 0.52 0.22 0.25 0.71 密度/(g/cm3) 2.42 2.67 2.79 1.83

表 2 万安盆地岩性物理参数[5]Table 2 Lithological physical parameters of Wan'an Basin[5]

在这里, 将万安盆地的构造沉降简化为几次盆地拉张的结果，以回剥法获得的构造沉降数据为约束，通过McKenzie纯剪切拉伸模型反复调整拉张因子拟合构造沉降史^[29]，最终得到盆地在该时间范围内的伸展特征。McKenzie纯剪切拉伸模型中S_t包括初始裂陷量(S_i)和裂后热沉降(S_T)^[30]

${S_t} = {S_i} + {S_T}.$

(4)

${S_i} = \frac{{{h_1}\left\{ {\left( {{\rho _{\rm{a}}} - {\rho _{\rm{c}}}} \right)\frac{{{h_{\rm{c}}}}}{{{h_1}}}\left( {1 - \alpha \frac{{{T_1}{h_{\rm{c}}}}}{{2{h_1}}}} \right) - \frac{{\alpha {T_1}{\rho _{\rm{a}}}}}{2}} \right\}\left( {1 - \frac{1}{\beta }} \right)}}{{{\rho _{\rm{a}}}\left( {1 - \alpha {T_1}} \right) - {\rho _{\rm{w}}}}}{\rm{ }}.$

(5)

${S_{T(t)}} = \frac{{4{h_1}{\rho _{\rm{a}}}\alpha {T_1}}}{{{{\rm{ \mathsf{ π} }}^2}\left( {{\rho _a} - {\rho _w}} \right)}}\left[ {\frac{{\beta }}{{\rm{ \mathsf{ π} }}}{\rm{ sin}}\left( {\frac{{\rm{ \mathsf{ π} }}}{\beta }} \right){\rm{ }}} \right]\left[ {1 - {e^{\frac{{t{{\rm{ \mathsf{ π} }}^2}K}}{{h_1^2}}}}} \right].$

(6)

式中:h₁为岩石圈初始厚；h_c为地壳初始厚度，此处在White和Bellingham^[31]提出的初始值的基础上以现今地壳厚度作为约束条件，分别取120、31km；ρ_c为地壳密度，取2.8g/cm³；T₁为软流圈顶界的温度取1330℃；α为岩石圈的热膨胀系数，取3.28×10^-5℃^-1；β为地壳的初始拉伸系数；K为岩石圈的热扩散率，取8.04×10^-7m²·s^-1[32]。盆地基底热流包括两部分，一部分为地壳岩石本身放射性衰变释放的热，一部分为由地球深处向上传导的地幔热^[33]。其中，地壳放射性衰变产生的热流为地壳厚度与地壳生热率的乘积，本文中地壳岩石的生热率是邻区地壳结构对比结果的平均生热率0.54mW/m^{3 [29]}，而来自地幔的热流则根据热力学的基本原理，通过岩石圈的热传导方程^[30]求取：

$Q\left( t \right) = \frac{{k{T_l}}}{h}\left\{ {1 + 2\sum\limits_{n = 1}^\infty {\frac{\beta }{{n{\rm{ \mathsf{ π} }}}}{\rm{sin}}\left( {\frac{{n{\rm{ \mathsf{ π} }}}}{\beta }} \right)} {\rm{ }}{{\rm{e}}^{\left( {\frac{{ - {n^2}t{\pi ^2}K}}{{h_1^2}}} \right)}}} \right\}.$

(7)

式中:k为岩石圈的热导率，取7.5×10^-3cal·cm^-1·s^-1；h为从地基到岩石圈底界的深度，cm。
3 研究结果3.1 地质剖面发育演化根据剖面位置(图 1(b))与地质剖面发育演化图(图 3)，可以得知出地层厚度分布不均匀，隆起区域地层具有大幅减薄的特征，而坳陷区地层厚度较大。地层厚度的变化除受构造位置的影响，也与断层活动密切相关，沿断层位置向两层变薄，下盘地层厚度减小的速度比上盘地层厚度减小的速度大。11Ma后，由于断层活动的减弱，断层处与其他位置的地层厚度差逐渐变小。万安盆地的构造格局受断层活动程度的影响，在相连地层内具有继承性特征，如西卫组与万安组、李准组与昆仑组、广雅组与第四系两两之间都具有相似的地层特征。其中西卫组与万安组受断层活动影响较大，断陷分块作用明显，李准组与昆仑组为断层活动减弱的过度地层，广雅组与第四系则为裂后沉积阶段，该阶段断层活动几乎停止，沉积物丰富，形成了较厚的覆盖地层。对比两条剖面并结合剖面穿过的二级构造，可发现剖面A-B的断裂数量比剖面C-D的断裂数量多，具有更为明显的裂陷分块特征，由此表明南部坳陷区域较中部坳陷区域受断层活动的影响更严重。综合剖面C-D的沉积厚度自NW至ES减小的趋势与剖面C-D比剖面A-B沉积层厚度大的现象，可发现万安盆地的沉积层厚度具有自盆地中心向盆地周缘递减的趋势。
Fig. 3

	Download: JPG larger image
图 3 地质剖面发育演化图 Fig. 3 Geological section development and evolution 图 3 地质剖面发育演化图 Fig. 3 Geological section development and evolution -->

3.2 沉降与伸展图 4(a)为利用剖面A-B与剖面C-D分别计算的平均总沉降量、构造沉降量、构造沉降速率绘制的折线柱状图。从图中可以看出剖面C-D的最终平均总构造沉降量约为2300m，平均总沉降量约为4500m，剖面A-B的最终平均总构造沉降量约为1850m，平均总沉降量约为3750m，其中构造沉降仅占总沉降量的50%~60%。构造沉积曲率呈现高低起伏的规律，可将其分为3个阶段：37.8—16Ma，16—5Ma，5Ma—现今，分别对应万安断裂的右旋、左旋、弱右旋3个阶段^[13]，也与南海的扩张和停止响应。5Ma—现今整个万安盆地处于向东倾斜的陆架上，物源丰富，此时万安断裂处于弱右旋，整个盆地构造活动弱，构造沉积速率达到3个阶段中的最大：剖面C-D可达210m/Ma，剖面A-B可达130m/Ma。16-5Ma期间南海扩张停止，周缘陆架受其影响沉降速率减少：剖面C-D约20m/Ma，剖面A-B约10m/Ma。在此期间，构造沉降曲线与总沉降曲线形成“凸”的形状，此时构造沉降量占总沉降量的比重偏低。37.8—16Ma期间，在南海扩张、洋脊跳跃和扩张方向的转换的影响下，周围陆缘构造活动增强，断层活动增加，沉降量增多。
Fig. 4

	Download: JPG larger image
图 4 沉降与伸展概况图 Fig. 4 Settlement and extension diagrams 图 4 沉降与伸展概况图 Fig. 4 Settlement and extension diagrams -->

由于第3次快速构造沉降作用为区域构造作用引起，非盆地本身构造作用所引起^[29]，且5.3Ma后，万安盆地伸展因子变化较小^[5]，所以本文将只计算5.3Ma前的伸展因子与速率。从图 4(b)可以看出在23—16Ma期间，盆地伸展速率达到最大，这与此时万安断裂处于右旋阶段，在研究区产生NW-SE向的扭张作用和该期间南海扩张脊向西南跃迁、西南次海盆由NE往SW方向渐进式的扩张有关。在16—11Ma，万安断裂处于左旋阶段，产生NW-SE向的扭压作用和南海海盆扩张停止，盆地伸展速率达到最小。剖面C-D的平均伸展因子为1.8，平均最大伸展速率为4.2mm/a；剖面A-B的平均伸展因子为1.6，平均最大伸展速率为3.9mm/a，对比剖面C-D、A-B的平均伸展因子、平均伸展速率，可发现剖面C-D比剖面A-B具有更加强烈的伸展活动。这与剖面所处的位置有关，剖面C-D与万安断裂的距离较剖面A-B与万安断裂的距离更短，受其活动的影响更大。并且剖面C-D处于南海西南次海盆扩张脊的延伸线上，受到南海东部次海盆和西南次海盆两次不同方向海底扩张产生的伸张作用。根据沉降量、伸展量与万安断层左右旋的对应关系——右旋期间盆地的沉降量、伸展量大于左旋期间的沉降量、伸展量(图 4)，我们提出:万安断裂左旋活动期间，在研究区产生扭压作用、万安盆地以地壳活动为主; 右旋活动期间，在研究区产生扭张作用、以地幔活动为主。
综合分析图 4(a)、4(b)，在37—23Ma期间剖面A-B与剖面C-D的构造沉降量、沉降速率与伸展速率几乎相同，说明该段时间整个万安盆地都处于伸展阶段，并未出现区域性不均匀伸展的特征。在23Ma后盆地的伸展具有明显的不均匀性，剖面C-D较剖面A-B具有较高的沉降、伸展，这说明盆地中部较南部具有更为剧烈的构造活动，万安盆地发育具有“北早南晚”的特征。结合姚永坚等^[5]对各条测线的伸展系数与伸展速度的计算结果等，得出渐新世初期万安盆地EW向的伸展速率大于NE-SW向的伸展速率。此时，EW向的伸展速率达到整个时间历程的最大值。NE-SW向的伸展速率在早中新世时期突然迅速增大达到峰值后降低，但速率整体上大于东西向。综合上述分析，我们认为万安盆地的伸展作用主要受西南次海盆扩张的影响，并以早东西向伸展为主，晚南北向伸展为主。
3.3 伸展速率与基地热流将计算得到的基底热流与沉积层生热相加，得到该区的表面热流，值为43~100mW/m²，该值与该地区真实的热流值所处区间相符，说明反演结果具有可靠性。为研究断层位置对伸展速率和基底热流的影响，将不同时间的伸展速率与基底热流网格化(基底热流取拉伸后的热流计算值)得到图 5。从图 5可以发现A-B的伸展速率为0~9mm/a、基底热流为20~70mW/m²，C-D的伸展速率为0~12mm/a、基底热流为25~100mW/m²，断层附近基底热流和伸展速率形成圈状异常。随着远离断层的位置，伸展速率、基底热流都逐渐变低。对比图 3可以发现伸展速率与基底热流异常区位于断层的下盘，该侧断层较多，构造活动较活跃，说明伸展速率和基底热流与构造活动正相关。图 5(a)、5(b)中伸展速率在23—16Ma较其他时间段伸展速率具有明显的异常，与南海扩张洋中脊的跳跃^[34]相对应。我们认为西北次海盆的南北扩张转向西南次海盆的NW-SE向扩张对万安盆地的伸展具有更大影响。
Fig. 5

	Download: JPG larger image
图 5 伸展速率与基地热流分布图 Fig. 5 Stretch rate and base heat flow distribution 图 5 伸展速率与基地热流分布图 Fig. 5 Stretch rate and base heat flow distribution -->

4 讨论印度板块楔入欧亚大陆^[35-36]，导致印支地块向SE方向挤出逃逸^[37-38]，沿红河断裂和莺歌海盆地1号断裂的左旋走滑和印支地块的旋转导致万安盆地的发育。32Ma时期古南海的俯冲控制万安盆地的持续伸展，到渐新世末期，古南海的俯冲导致洋脊跃迁，中央海盆扩张，并向西南方向扩展，西南次海盆逐渐发育^[34]，对万安盆地的伸展具有重要的作用。早中新世末期沙巴造山运动发生，南沙地块和婆罗洲地块及卡加延地块碰撞，南海扩张停止，万安盆地全面进入裂后热沉降阶段，该区域断层在左右旋转化过程中对盆地有改造作用。总之，新生代以来，万安盆地主要划分为：初始裂陷期、断层改造期、裂后沉降期，坳陷区还具有“下断上坳”的特征。
在前人研究的基础上，根据所选剖面有关沉降量、伸展、热流的计算结果，本文绘制万安盆地的沉降、伸展程度与盆地内断层位置的关系的模式示意图，见图 6(参考前人有关万安盆地地层剖面的研究成果^{[11, 39]})。图中虚线代表万安盆地内的断层，黑线为万安盆地的地层层序，黑点的密度代表异常值——盆地沉降量、伸展系数、基底热流(由于三者之间正相关，所以三者统称异常值)的大小。从图中可以看出万安盆地沉积层的平均基底埋深约为7~8km^[11]，大多数地区的T5和Tg之间的人骏群处于缺失状态，断裂构造极为发育，以正断层为主，坳陷区域具有“下断上坳”的特征。万安盆地莫霍面的平均深度约为24km^[11]，且在重力均衡的作用下坳陷区莫霍面深度浅，隆起区域莫霍面深。整体上，万安盆地的沉降与伸展活动则在万安断裂右旋期间T5-T31(37—16Ma)伸展活动较剧烈，与南海扩张和洋脊的跳跃相关，我们认为此时地幔活动强于地壳活动。局部上万安盆地沉降、伸展、热流有以断裂带为中心向两侧逐渐减少的趋势，断裂发育密集的区域异常(沉降、伸展、热流)强，反之较弱。
Fig. 6

	Download: JPG larger image
图 6 万安盆地沉降模式示意图 Fig. 6 Schematic of the settlement mode of Wan'an Basin 图 6 万安盆地沉降模式示意图 Fig. 6 Schematic of the settlement mode of Wan'an Basin -->

5 结论万安盆地作为南海西南陆缘的重要沉积盆地，具有丰富的油气资源。本文基于构造格局、地层层序等区域地质资料，利用回剥法和McKenzie模型对万安盆地的沉降、伸展特征、基底热流、动力学来源进行研究，主要结论如下：
1) 万安盆地演化阶段主要划分为：初始裂陷期、断层改造期、裂后沉降期，地层厚度分布不均匀，存在自盆地中心向盆地周缘递减的趋势。
2) 整体上万安盆地的沉降、伸展作用主要受西南次海盆扩张和万安断裂走滑运动的影响，具有独特的时空演化特征：沉降作用整体上呈“V型转折”的特性，而伸展作用以“早东西、后南北”为主要特点；局部上万安盆地的沉降、伸展呈现出以断裂带为异常轴，向两侧以不同速率逐渐减弱的特征。
3) 根据沉降量与万安断层左右旋的对应关系——右旋期间沉降、伸展大于左旋期间，认为左旋活动期间产生扭压作用、万安盆地以地壳活动为主，右旋活动期间产生扭张作用、以地幔活动为主。

参考文献

[1]	任纪舜, 王作勋, 陈炳蔚, 等. 新一代中国大地构造图[J]. 中国区域地质, 1997, 16(3): 2-7.
[2]	Kusumastuti A, Rensbergen P V, Warren J K. Seismic sequence analysis and reservoir potential of drowned miocene carbonate platforms in the Madura Strait, East Java, Indonesia[J]. Aapg Bulletin, 2002, 86(2): 213-232.
[3]	姚伯初. 南海西南海盆的海底扩张及其构造意义[J]. 南海地质研究, 1997(9): 20-36.
[4]	Briais A, Patriat P, Tapponnier P. 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, 1993, 98(B4): 6299-6328. Doi:10.1029/92JB02280
[5]	姚永坚, 吕彩丽, 王利杰, 等. 南沙海区万安盆地构造演化与成因机制[J]. 海洋学报, 2018, 40(5): 62-74.
[6]	施秋华, 夏斌, 万志峰, 等. 万安盆地构造沉降与沉积速率及其油气资源前景[J]. 沉积与特提斯地质, 2014, 34(2): 43-46.
[7]	刘振湖. 南海万安盆地沉降作用与油气前景[J]. 中国海上油气地质, 1998, 12(4): 19-25.
[8]	杨楚鹏, 姚永坚, 李学杰, 等. 万安盆地新生代层序地层格架与岩性地层圈闭[J]. 地球科学——中国地质大学学报, 2011, 36(5): 845-852.
[9]	杨木壮, 王明君, 梁金强, 等. 南海万安盆地构造沉降及其油气成藏控制作用[J]. 海洋地质与第四纪地质, 2003, 23(2): 85-88.
[10]	寻博靖.万安盆地重磁场特征研究与构造分析[D].西安: 西安石油大学, 2017.
[11]	杨胜雄, 邱燕, 朱本铎. 南海地质地球物理图系[M]. 天津: 中国航海图书出版社, 2015.
[12]	张厚和, 赫栓柱, 刘鹏, 等. 万安盆地油气地质特征及其资源潜力新认识[J]. 石油实验地质, 2017, 39(5): 625-632.
[13]	安慧婷, 李三忠, 索艳慧, 等. 南海西部新生代控盆断裂及盆地群成因[J]. 海洋地质与第四纪地质, 2012, 32(6): 95-111.
[14]	Liu H, Yan P, Zhang B, et al. Role of the Wan-Na fault system in the western Nansha Islands (southern South China Sea)[J]. Journal of Asian Earth Sciences, 2004, 23(2): 221-233. Doi:10.1016/S1367-9120(03)00121-4
[15]	Morley C K. A tectonic model for the Tertiary evolution of strike-slip faults and rift basins in SE Asia[J]. Tectonophysics, 2002, 347(4): 189-215. Doi:10.1016/S0040-1951(02)00061-6
[16]	刘海龄, 姚永坚, 沈宝云, 等. 南海西缘结合带的贯通性[J]. 地球科学——中国地质大学学报, 2015, 40(4): 615-632.
[17]	Taylor B, Hayes D E. The tectonic evolution of the South China Basin[J]. Tectonic & Geologic Evolution of Southeast Asian Seas & Islands, 1980, 23: 89-104.
[18]	Watts A B, Ryan W B F. Flexure of the lithosphere and continental margin basins[J]. Tectonophysics, 1976, 36(1): 25-44.
[19]	Bond G C, Kominz M A. Construction of tectonic subsidence curves for the early Paleozoic miogeocline, southern Canadian Rocky Mountains:implications for subsidence mechanisms, age of breakup, and crustal thinning[J]. Geological Society of America Bulletin, 1984, 95(2): 155-173. Doi:10.1130/0016-7606(1984)95<155:COTSCF>2.0.CO;2
[20]	Mcintyre A. Basin analysis:principles and applications[J]. Economic Geology, 2006, 113(2): 409-442.
[21]	张云帆, 廖杰, 孙珍, 等. 南沙海域构造沉降特征[J]. 地球科学——中国地质大学学报, 2011, 36(5): 949-955.
[22]	赵中贤, 孙珍, 陈广浩, 等. 南沙海域新生代构造特征和沉降演化[J]. 地球科学——中国地质大学学报, 2011, 36(5): 815-822.
[23]	Haq B U, Hardenbol J, Vail P R. Chronology of fluctuating sea levels since the Triassic[J]. Science, 1987, 235(4793): 1156-1167. Doi:10.1126/science.235.4793.1156
[24]	Athy L F. Density, porosity, and compaction of sedimentary rocks[J]. Aapg Bulletin, 1930, 14(1): 1-24.
[25]	King H M, Rubey W W. Role of fluid pressure in mechanics of overthrust faulting:I. Mechanics of fluid-filled porous solids and its application to overthrust faulting[J]. Geological Society of America Bulletin, 1959, 70(5): 115-116.
[26]	Clift P, Lin J. Preferential mantle lithospheric extension under the South China margin[J]. Marine & Petroleum Geology, 2001, 18(8): 929-945.
[27]	Sclater J G, Christie P A F. Continental stretching:explanation of post-mid-cretaceous subsidence of Central North Sea Basin:Abstract[J]. Aapg Bulletin, 1980, 85(B7): 3711-3739.
[28]	Salisbury M H, Shinobara M, Richter C, et al. Proceedings of the ocean drilling program, initial reports[R].Texas, College Station: Ocean Drilling Program, 2000: 184-185.
[29]	刘振湖. 北康盆地古地热场与油气远景[J]. 海洋地质与第四纪地质, 2004, 24(2): 79-84.
[30]	McKenzie D. Some remarks on the development of sedimentary basins[J]. Earth and Planetary Science Letters, 1978, 40(1): 25-32.
[31]	White N, Bellingham P. A two-dimensional inverse model for extensional sedimentary basins 1.Theory[J]. Journal of Geophysical Research:Atmospheres, 2002, 107(B10). Doi:10.1029/2001JB000173.:2259-2279
[32]	李思田, 解习农, 王华, 等. 沉积盆地分析基础与应用[M]. 北京: 高等教育出版社, 2004: 275-288.
[33]	汪集旸, 等. 地热学及其应用[M]. 北京: 科学出版社, 2015.
[34]	李家彪, 丁巍伟, 高金耀, 等. 南海新生代海底扩张的构造演化模式:来自高分辨率地球物理数据的新认识[J]. 地球物理学报, 2011, 54(12): 3004-3015.
[35]	Patriat P, Achache J. India|-Eurasia collision chronology has implications for crustal shortening and driving mechanism of plates[J]. Nature, 1984, 311(5987): 615-621. Doi:10.1038/311615a0
[36]	Aitchison J C, Ali J R, Davis A M. When and where did India and Asia collide?[J]. Journal of Geophysical Research Solid Earth, 2007, 112(B5): 51-70.
[37]	Charusiri P, Imsamut S, Zhuang Z, et al. Paleomagnetism of the earliest Cretaceous to early late Cretaceous sandstones, Khorat Group, Northeast Thailand:implications for tectonic plate movement of the Indochina block[J]. Gondwana Research, 2006, 9(3): 310-325. Doi:10.1016/j.gr.2005.11.006
[38]	Takemoto K, Halim N, Otofuji Y I, et al. New paleomagnetic constraints on the extrusion of Indochina:late Cretaceous results from the Song Da terrane, northern Vietnam[J]. Earth & Planetary Science Letters, 2005, 229(3/4): 273-285.
[39]	吴峧岐, 高红芳, 孙桂华. 南沙海域万安盆地地质构造与沉积体系特征[J]. 海洋地质与第四纪地质, 2012, 32(5): 1-11.