1.   引言

华南大陆主体由扬子地块和华夏地块构成，NE—NEE走向的钦杭结合带被认为是扬子和华夏两个陆块新元古代的碰撞拼合带（杨明桂等, 2009）。华南地区中、新元古代总体是在全球Rodinia大陆聚合与裂解的构造背景中发展演化的，古生代经历了广泛的陆内造山过程，中生代以来处在全球三大重要板块的汇聚拼合部位，并经历了印支期的碰撞造山和燕山期的陆内造山过程（张国伟等，2013）。由于构造演化历史漫长且十分复杂，扬子地块与华夏地块的接触位置和耦合关系尚有不同的看法（舒良树，2012），块体内部组成有待进一步深化。

泊松比是反映地球物质成分极为重要的参数，测量泊松比是研究地球内部物质成分、构造和物理状态最有效的方法之一（嵇少丞等，2009a）。地壳泊松比的大小取决于岩石中P波和S波的速度比值，即波速比。泊松比和波速比的大小呈正相关，流体的存在能有效地提高岩石的波速比，从而增加泊松比。前人总结了岩石圈中主要造岩矿物和岩石的泊松比-密度的关系（Ji et al., 2002），并以华北地区为例，开展了地壳泊松比异常成因与大地构造意义等有关的研究（嵇少丞等，2009b；杨文采等，2011），前人研究多数使用了地壳天然地震的P、S波走时成像反演方法。

华南天然地震研究成果近年来逐渐丰富，主要是通过远震资料研究华南岩石圈-地幔的速度结构（黄海波等, 2014；He et al., 2015, 2016; 张耀阳等, 2018）。通过近震资料研究华南地壳的速度结构、波速比、泊松比等目前还局限在较小范围，多以单独省份资料独立研究（黄耘等, 2008；李红星等, 2015）。

本文以华南为研究区，系统地搜集、整理了研究区内10多个省份的地壳天然地震观测资料，提取直达Pg、Sg波走时，以及沿Moho面滑行传播的首波Pn、Sn波走时。通过对走时大数据进行拟合分析，得到华南地壳的天然地震平均波速和波速比，并以此值作为地壳三维波速比反演时的初始值，反演得到全地壳的波速比，根据公式转换得到华南地壳的泊松比。在此基础上，分析了华南地壳物质分布、深部岩浆过程及与主要成矿带的关系。

2.   华南天然地震数据分析

2.1   近震数据及筛选

本次研究系统地搜集整理了华南天然地震震相数据，研究数据主要来源于国家地震科学数据共享中心（http://data.earthquake.cn/index.html）、地震专业知识服务平台（http://earthquake.ckcest.cn/index.html）等数据共享平台，以及部分省属台网。

首先，对获取的原始震相报告进行数据结构分析，提出了编写数据处理程序的算法，算法主体由两部分组成：①对包含地震发震时刻、震中参数等起始行进行信息提取；②对该次地震所有记录台站的到时参数、台站参数进行提取，同时计算该台站的Pg、Sg走时和Sg-Pg走时差，如果该台站有Pn、Sn波数据，同步处理Pn、Sn走时数据。在数据处理过程中，采取严格的筛选措施，仅当某一次地震满足4对以上Pg、Sg走时记录的地震震相数据才予以使用。经过大数据分析、提取和筛选，本次研究得到华南研究区内M≥2.0级地震事件11410个，详见图 1。

图  1  华南研究区内地震震中分布及主要断裂

图中断裂及构造块体边界：红色线条据严加永等（2019）；黑色线条据邓起东等（2003）

Figure  1.  Distribution of earthquakes and major faults in South China

Faults and block boundary are adapted: red lines after Yan et al., 2019; black lines after Deng et al., 2003

下载: 全尺寸图片 幻灯片

研究区地貌总体以“盆”“岭”相间为特征，西部的江南造山带和东部武夷山地形相对较高，南北向以南岭为界，向两侧逐渐降低。从地震震中空间分布看，东南沿海福建、广东、广西等省份及邻近海域地震活跃，研究区中部的湖南、贵州等省份地震活动较弱，北部的重庆、湖北、江西、安徽等省份地震活动较活跃。钦杭结合带处于研究区中部，自广西北部湾海域的钦州至浙江杭州湾，钦杭结合带西南段总体呈NE向分布，北东段总体呈NEE向分布。华夏地块位于钦杭结合带的东南侧，扬子地块位于钦杭结合带的北西侧。钦杭结合带的北东段出露绍兴—江山—萍乡和东乡—德兴—歙县两条新元古代早期的蛇绿混杂岩带（舒良树，2012），多数研究者认为它们是碰撞结合带的主要边界断裂；钦杭结合带的西南段主要边界断裂分布、产状等还存在较大争议。

2.2   走时数据拟合与射线分布

为了拟合各震相走时数据，首先需要对不同震相的走时数据进行分离，笔者通过编写程序分离出Pg、Sg、Pn和Sn的走时，对各类走时数据计算走时与震中距的关系，分别使用最小二乘法拟合数据点，得到一维拟合直线方程，由拟合斜率即可以求得各震相在地壳或上地幔的平均波速。

首先，笔者对1980—2017年的所有数据进行分析，自2000年以后的数据质量有明显提高，特别是2008年之后由中国地震台网中心统一定位并发布的数据质量很好，走时-震中距关系数据点非常收敛并且离散的飞点数据极少。综合考虑之后，此次研究使用了华南研究区内2000—2017年的M≥2.0级地震走时数据。

此次研究共使用了11410个地震，挑选Pg走时数据140600余条，Sg走时数据164500余条，Pn走时数据64500余条，Sn走时数据27700余条。对不同震相的走时-震中距关系数据点（图 2），分别使用最小二乘法进行直线拟合，通过拟合直线方程的斜率得到华南地壳Pg波平均波速为6.11 km/s, Sg波平均波速为3.60 km/s, Pn波平均波速为8.06 km/s, Sn波平均波速为4.59 km/s。

图  2  华南研究区内走时-震中距关系曲线

Figure  2.  Travel times of different seismic phases observed in South China

下载: 全尺寸图片 幻灯片

图 3是此次研究使用的地震震中、Pg、Sg波射线和观测台站分布。根据数量分布统计，Pg、Sg波射线在东南沿海的福建、广东、广西覆盖最好，其次为研究区西北的重庆、湖北一带，赣东北—皖南—浙西一带以及湖南—贵州交界附近射线覆盖偏少。

图  3  华南Pg、Sg波地震射线及台站分布

Figure  3.  Distribution of the seismic stations and ray paths of Pg and Sg phases of South China

下载: 全尺寸图片 幻灯片

3.   华南地壳平均波速比与泊松比

3.1   华南地壳平均波速比

通过上一节的工作，得到华南地壳的直达Pg波平均波速为6.11 km/s, 直达Sg波平均波速为3.60 km/s, 沿壳幔界面滑行的首波Pn波平均波速为8.06 km/s, Sn波平均波速为4.59 km/s。在此基础上，求得华南地壳平均波速比Vp/Vs为1.697，壳幔界面附近的平均波速比Vp/Vs为1.756。

3.2   华南地壳平均泊松比

由P、S波的平均波速比，按照公式（1）可计算平均泊松比σ：

华南地壳的平均泊松比为0.234，地幔顶部的平均泊松比为0.260。研究表明，虽然泊松比受很多因素的影响，比如流体/熔体的存在、岩石内部的显微结构等，但总体上，泊松比是地壳物质成分的反映。长英质酸性岩组成的地壳其σ≤0.26；铁镁质基性岩组成的地壳其σ值在0.28~0.30；而中性岩组成的地壳或酸性岩和基性岩按照一定的比例混合组成的地壳其σ值在0.26~0.28（Christensen and Mooney, 1995; 刘斌等, 2002; 嵇少丞等, 2009）。华南地壳的平均泊松比σ为0.234，与花岗岩的泊松比接近，总体上与华南出露花岗岩总面积169.7 km²，整个地壳充满长英质酸性物质的事实是吻合的。

4.   华南中、上地壳三维泊松比及地质意义

4.1   泊松比的反演计算

由于Pg、Sg走时数据的质量和数量差异，直接使用地震走时Pg、Sg分别反演P波和S波速度结构，通过简单的Vp和Vs比值来得到波速比会带来很多问题。Thurber等（1993）提出了一个好的计算方法，通过Pg走时和Sg-Pg走时差来获得可靠的波速比。公式如下：

公式（2）中t_ij表示第i个地震事件在第j个台站的Pg或Sg波走时，dt_ij为Sg-Pg的走时差。反演计算流程为：首先通过震源位置和速度结构的联合反演得到三维P波速度结构，以三维P波速度结构和预设的Vp/Vs波速比初始值作为输入参数，通过（2）式计算理论的Sg-Pg走时差dt_ij，求取理论走时差和实际观测数据的走时差dt_ij*的差值△t_ij，以△t_ij扰动Vp/Vs波速比进行迭代反演，多次迭代时的震源位置和三维P波速度结构为固定值，当走时残差△t_ij满足收敛条件时终止迭代计算，得到Vp/Vs波速比反演结果。

反演计算使用Simulps14软件（Thurber，1993；Zhang，2003）。为保证反演质量，要求一次地震至少有4个以上台站记录到，并有准确的Pg走时和Sg-Pg走时差才能入选计算数据。初始波速比通过走时大数据的拟合得到，地壳初始波速比Vp/Vs为1.697，壳幔界面附近的初始波速比Vp/Vs为1.756。一维初始地壳速度模型参考黄启勋（2000）、徐涛等（2014）等人工地震剖面成果。反演使用自中心点向四周1°×1°的网格划分，图 4是通过三维波速比反演计算得到的泊松比分布。由于研究区内地震射线主要分布在中、上地壳，因此反演结果主要揭示中、上地壳泊松比的分布。

图  4  华南中、上地壳5 km（a）和14 km (b)深度泊松比的分布

Figure  4.  Map of the Poisson's ration of 5 km (a) and 14 km (b) in the mid-upper crust of South China

下载: 全尺寸图片 幻灯片

4.2   解的可靠性评价

Simulps14软件对于解的可靠性评价主要采用分辨率对角元素值（RDE）和射线分布偏导权重总数（DWS），解的可靠性与RDE和DWS值呈正相关。图 5分别是5 km和14 km RDE和DWS的分布。

图  5  不同深度的走时成像DWS（a、b）和RDE（c、d）分布

Figure  5.  DWS (a and b) and RED (c and d) distribution in different depths from the travel-time tomography

下载: 全尺寸图片 幻灯片

由图 5可知，台湾附近海域由于射线没有形成良好的覆盖，解的可靠性较低。湖南和贵州两省交界区域，由于地震活动微弱，台站也比较稀疏，射线密度和解的可靠性都较低。东南沿海的福建、广东、广西南部射线覆盖比较好，解的可靠性较高。文中截取的AA’和BB’两条剖面，综合考虑了射线分布、解的可靠性和地质问题等因素。

4.3   华南中上地壳平均泊松比及地质意义

图 4、图 6分别给出了5 km、14 km深度的地壳的泊松比和20 km以浅深度的地壳平均泊松比，由图可见，5 km深度水平切片上，南岭成矿带、武夷山成矿带下方泊松比未见明显的高值区，随着深度增加到14 km的水平切片上，武夷山成矿带、长江中下游成矿带和南岭成矿带的下方都可见到成片分布的高泊松比异常，而且高泊松比区域与东南沿海火山岩带、华夏地块的对应关系更为明显。高泊松比（局部可达0.28）指示地壳成分以镁铁质为主，与区内广泛分布的基性火山岩特征一致，指示深部有可能发生过镁铁质岩浆的底侵作用（Li et al., 2018；Lin et al., 2018）。扬子地块中、上地壳的泊松比以高值为主，高值接近0.29，可能与扬子克拉通的基性物质组成有关，但仍存在局部（西北角的重庆附近）较低的泊松比，反映局部物质的不均匀特征。

图  6  华南中、上地壳平均泊松比异常水平分布与解释

Figure  6.  Map of the average Poisson's ration in the mid-upper crust of South China and its interpretation

下载: 全尺寸图片 幻灯片

江南造山带处在扬子地块和华夏地块之间，是新元古代时期发生在扬子地块东南缘以岛弧地体拼接为主的增生型造山带，目前围绕江南造山带的构造-岩浆演化还存在较大争议（王孝磊等，2017）。扬子地块和华夏地块在新元古代拼合带的位置在研究区内北东段比较清楚，但是在过湖南萍乡之后，地表植被广泛覆盖，缺少关键地质标志，而使得该拼合带的西延界限难以确定（舒良树，2012）。造山过程往往会改变地壳的泊松比，受到构造挤压，处于相同温度、压力条件下的长英质岩石比基性岩石更容易形成褶皱和推覆构造，其结果势必造成地壳泊松比随地壳厚度增加而减少（嵇少丞等, 2009）。本次研究发现在华夏地块与扬子地块结合带中间存在与钦杭结合带近乎平行的低泊松比区，该低泊松比异常区在浙江、安徽、江西、湖南等省份的边界与已知的江南造山带的地质边界（张国伟等，2013）基本一致，因此，笔者认为这一低泊松比异常区反映了江南造山带的大致范围，并据此勾勒出了江南造山带的大致边界位置（图 6）。He et al.（2015）通过远震接收函数方法得到的地壳平均波速比也有类似的结果。

南岭成矿带主要分布矽卡岩型和石英脉型钨锡矿床，钦杭成矿带（北东段）主要分布斑岩-矽卡岩型钨、铜矿床，武夷山成矿带主要分布斑岩-矽卡岩型铜矿床以及浅成低温热液型金矿床（毛景文等，2018），长江中下游成矿带主要为矽卡岩型、斑岩型、陆相火山岩型的铜、铁、金多金属为特色的矿床（常印佛等，2017）。由图 6可知，江南造山带南段整体呈较低的泊松比异常（σ < 0.21），南岭成矿带局部具有相对较高的泊松比异常（σ>0.23），结合地质资料分析，相对高泊松比异常可能与南岭成矿带燕山期大规模花岗岩侵位相关，燕山期成矿花岗岩具有幔源信息（周永章等，2017）。武夷成矿带和钦杭成矿带（北东段）地理位置上毗邻，都是重要的斑岩型铜（钼）矿带，带内的斑岩型矿床具有元古宙+中生代双峰年龄，地球化学特征上带有岛弧俯冲的基因（周永章等，2017），泊松比异常背景值约为0.23，未见到显著的高泊松比异常区。长江中下游成矿带泊松比异常背景值也约为0.23，高泊松比异常主要沿西北侧边界郯庐断裂带附近分布。

为便于从剖面上分析，笔者沿钦杭结合带截取了2条泊松比剖面，分别揭示钦杭结合带及两侧不同地质块体之间的相互关系（图 7）。

图  7  AA’和BB’剖面泊松比垂直切片图

Figure  7.  Vertical slices of Poisson's ratio along AA'and BB'profiles

下载: 全尺寸图片 幻灯片

AA’剖面与钦杭结合带北段实施的地球物理廊带剖面之一的宁国—嵊州剖面位置一致，有所延长（中国地质科学院，2018^❶）。AA’剖面揭示了，在以钦杭结合带为中心的西北（A方向）和东南（A’方向）两翼的中、上地壳泊松比异常迥然不同，西北侧对应的扬子地块具有相对高的泊松比异常，而东南侧对应的华夏地块具有相对低的泊松比异常。垂向上，从浅到深，泊松比明显增高，显示基底（10 km）以深基性物质成分明显增加，笔者解释为扬子与华夏地块碰撞拼贴过程中，幔源物质不断注入改造地壳成分的结果。

BB’剖面显示在剖面中心约500 km以西的中、上地壳以低泊松比异常为主，间或夹有局部高泊松比异常，在其以东的中、上地壳主体为相对高的泊松比异常，东南沿海的地壳具有相对高的泊松比异常（0.26~0.28），笔者解释为地壳铁镁质含量增加，与幔源物质底侵过程有关。BB’剖面左侧的低泊松比异常区域（0.17~0.22），反映中上地壳长英质成分为主，推测与江南造山带地壳增厚有关，该异常与扬子地块也有显著不同。从图 6泊松比平面图和AA’剖面图可知扬子地块相对华夏地块具有更高的泊松比，而这一夹在扬子地块与华夏地块之间的江南造山带具有显著的低泊松比异常，笔者认为江南造山带西段，由于造山过程是中上地壳物质变形增厚，总体降低了地壳平均泊松比；而造山带东段出现的相对高的泊松比，解释为燕山期陆内造山阶段大量幔源物质的注入，改变了地壳的基性程度。

5.   结论

（1）华南地壳及上地幔顶部平均波速：Pg波6.11 km/s, Sg波3.60 km/s, Pn波8.06 km/s, Sn波4.59 km/s。华南地壳平均波速比Vp/Vs为1.697，壳幔界面附近的平均波速比Vp/Vs为1.756。

（2）华南地壳平均泊松比为0.234，指示地壳中长英质矿物含量较高，地幔顶部的平均泊松比为0.260。

（3）东南沿海和武夷山具有相对高的地壳泊松比，认为与幔源基性物质的底侵有关。

（4）江南造山带处于扬子地块与华夏地块之间，具有显著的低泊松比异常，与造山过程导致中上地壳增厚有关。

注释

❶中国地质科学院.2018.钦杭结合带及邻区深部地质调查成果报告[R].