当涂地震站地方震Sc震相分析及康拉德面特征

陈江琴

1）中国安徽 233500 安徽蒙城地球物理国家野外科学观测研究站

2）中国安徽 243000 马鞍山市应急管理局（地震局）

0 引言

地震波在传播过程中，携带着大量来自震源及传播路径上的地球内部构造信息。研究地震波的性质、特征，是研究地壳结构的有效办法（董奇珍等，1989），而震相分析对深入研究地球内部构造和探查地震活动规律具有重要意义（陈祥开等，2018）。有学者在地震观测中发现地壳中存在康拉德面及其反射波，并用于地震研究，如：董奇珍等（1989）在甘肃成县爆破中发现来自地壳中间界面的反射波PD、SD，且较其他界面反射波阻抗大、反射力强；郭晓然等（2016）在西沙地块地壳中也发现了康拉德面反射波Pc；许健生等（1999）利用兰州地震台记录的永登地区17 个地震的康拉德面反射波Sc 走时特性，测算出该地区地壳中间层深度为（24.9±2.0）km；郭杰等（2013）利用河南濮阳地震台记录，观测到康拉德面反射波Pc，并认为对于震中距约40 km 的地震，地震台站可清晰记录到康拉德面反射波。

康拉德面是指地壳内花岗岩与玄武岩的分界面，其存在与否与研究区地壳结构及其形成和演化的历史密切相关。由于康拉德面是不连续面，在地震观测中来自该面的震相记录并不常见，且易被误认为来自莫霍面，不仅影响地震定位结果，还严重影响对地壳内部结构的认识。

胡德昭等（1989）发现，中国东南部地区地壳内广泛存在康拉德面，其分布深度由内陆向沿海逐渐变浅。戚浩等（2012）在安徽地震台网地方震记录中发现了来自康拉德面的Pb震相，郑晔等（1989）在随县—马鞍山人工源地震探测资料中发现马鞍山—六安剖面上存在Pb 震相，进一步验证了安徽地区地壳中应存在康拉德面。马鞍山测震台网囊括安徽省8 个地震台站（台站代码CHZ、HEF、HSH、MAS、DTT、TOL、JIX、WHU）和江苏省3个地震台站（台站代码GC、LSH、LY），可记录马鞍山地区ML1.0 以上地震。安徽省马鞍山市和芜湖市地处长江南岸，地质构造分区属下扬子断块，位于郯庐断裂带东侧、扬州—铜陵地震带中段。对于该区发生的地震，在满足一定条件下，马鞍山测震台网可记录到康拉德面上的反射波。文中采用当涂地震站对2016 年芜湖市弋江区ML2.1 地震的数据记录，分析Sc 震相特征，并利用马鞍山测震台网记录的6 个具有Sc 震相的典型震例，分析康拉德面特征，为地震研究和地球深部构造研究提供参考。

1 Sc 震相判定

据安徽地震台网测定，2016 年2 月13 日4 时17 分58 秒芜湖市弋江区发生ML2.1 地震，震源深度8 km，当涂地震站（震中距Δ=43.75 km）清晰记录到此次地震波形，且Pg、Sg 震相到时确定，见图1。由图1 可见Pg、Sg 及i 三个明显震相，其中i 震相位于Sg 波段后，见图1 中方框所示，其振幅比Sg 振幅大，而周期相近。

图1 2016 年2 月13 日4 时17 分芜湖市弋江区ML 2.1 地震波形Fig.1 The waveform of the ML 2.1 earthquake in Yijiang District,Wuhu City at 4:17 on February 13,2016

图1 所示地震波似双震波形，由马鞍山测震台网波形记录可知当地未发生双震，排除双震可能。因地方震无面波，仅产生直达波和反射波，由此确定该不明震相应为横波反射波。一般，莫霍面上的反射波在震中距70—80 km 时达到全反射，但此次地震距当涂地震站43.75 km，表明地震波在更小距离上发生了全反射，说明反射波走时缩短，由此推断出在莫霍面上方存在一个使地震波发生全反射的界面，且该界面波阻抗大、上下面速度对比强、反射力强。

设震源距地面距离为h，台站与震中距离为Δ，地面与莫霍面距离为H，与莫霍面上方的中间界面C 的距离为HC。台站接收的直达波为Pg、Sg，来自莫霍面的反射波为PmS、SmS，来自地壳中间界面的反射波为Pc、Sc，则直达波和反射波传播路径见图2。

图2 直达波和反射波的传播路径示意Fig.2 Schematic illustration of the propagation path of direct and reflected waves

直达波走时方程为

莫霍面上反射波走时方程为

假设该不明震相为莫霍面上的反射波，因震中距Δ＜48.96 km，由安徽地震台站所用华南走时表可知无对应走时，可按照安徽地区地壳一维速度最优模型（谢石文等，2016）（该模型更符合安徽地区实际情况）计算该该震相走时。该模型地壳中间面上层的vPg=6.14 km/s，vSg=3.57 km/s，中间地壳深19 km，莫霍面上层的vPg=6.60 km/s，vSg=3.84 km/s，莫霍面深34 km，将Δ=43.75 km、h=8 km 代入式（1）、（2），计算得到tPg=7.24 s、tSg=12.46 s。由式（3）、（4）计算得到tPmP=11.25 s、tSmS=19.34 s。利用震相实际到时与发震时刻的差值，计算得到各震相实际走时为：tPg=7.04 s、tSg=12.23 s、ti=16.7 s。则莫霍面理论走时与实际走时差为：ΔtPg=0.20 s，ΔtSg=0.23 s，ΔtSmS=2.64 s。由于地下地质结构和地震射线路径的差异，直达波走时差0.20 s 和0.23 s 均在正常范围内，而莫霍面反射波走时比实际走时长2.64 s，差值太大，说明该不明震相的反射界面不应为莫霍面，而是深度比莫霍面更浅的界面。

在华南走时表中，震中距Δ=39.96 km 时即有对应的康拉德面反射波走时。为此，查询华南走时表可知：tPg=7.39 s，tSg=12.53 s，tPc=9.64 s，tSc=16.39 s，与实际走时差分别为：ΔtPg=0.35 s，ΔtSg=0.30 s，ΔtSc=-0.31 s，可见差值0.31 远小于2.64。地下地质结构复杂，参与计算的速度值和深度值均为平均值，并非绝对精确值，且走时表并非绝对适合本地区地壳结构，对于某次地震而言，其地震波的实际走时和理论走时不可能绝对吻合，由此推断该不明震相i 为地壳中间层界面反射波更恰当。

2 反射纵波位置确定

近震震相一般成对出现，但在该地震波形图上未观测到明显纵波反射波，究其原因可能是：①纵波反射波不发育；②纵波反射波可能在Sg 波之后到达，与直达横波震相混在一起难以区分。当震中距达到一定距离时，反射纵波在直达横波之前出现，小于该距离则被淹没在直达横波里，此距离即为干扰距离Δ′，即纵波反射波走时与直达横波走时相等时的震中距。直达横波Sg 和纵波反射波P反射走时曲线见图3，可见：当Δ＜Δ′时，P反射在Sg 之后出现；当Δ＞Δ′时，P反射在Sg 之前出现；当Δ=Δ′时，有tSg=tPmP，则由式（2）、（3）可得

图3 直达横波Sg 和纵波反射波P反射走时曲线Fig.3 Travel time curves of direct shear wave Sg and reflected longitudinal wave P

因Δ=Δ′，则有

由式（5）可知，干扰距离Δ′与反射界面深度H、震源深度h及波速有关。当H=34 km（莫霍面深度）、h=8 km、vSg=3.57 km/s、vP反射=6.60 km/s 时，Δ′=41.77 km；当H=19 km（康拉德面深度）、h=8 km、vSg=3.57 km/s、vP反射=6.14 km/s 时，Δ′=19.05 km。而当涂地震站震中距Δ=43.75 km ＞ 41.77 km ＞ 19.05 km，因此，莫霍面或康拉德面的纵波反射波，均在直达波之前到达，也就是说，不明震相到时在直达横波之前。

由图1 所示波形可见纵波反射波不明显，将波形进行W.A.仿真，结果见图4，可见NS 向波形在到时04:18:7.16 处可见较大波动，如图4 箭头所指，垂直向波形在该处同样可见较大振幅，计算得到该震相走时为9.16 s，查询华南走时表可知Pc 震相走时为9.64 s，二者相差0.48 s，可见两震相走时相近，且上述已知莫霍面纵波反射波走时为11.25 s，与实际走时相差2.09 s。由此推断此处应为Pc 震相。

图4 仿真后地震波形Fig.4 Seismic waveforms after simulation

3 康拉德面特征分析

3.1 选取数据

由于Sc 波比Sg 波到时晚，其波形多淹没在Sg 波里，导致清晰的Sc 波初动较少见。除芜湖市弋江区地震外，笔者在马鞍山测震台网记录中统计得到5 个Sc 波初动较清晰的震例，地震参数见表1。利用6 个地震的Sc 震相参数，分析康拉德面特征。受篇幅所限，只给出其中3 个地震的波形，见图5。

表1 6 个含Sc 震相的地震参数Table 1 Parameters of 6 earthquakes with the Sc seismic phase

图5 3 个震例波形（a）2020 年12 月25 日17:54:56 南京江宁地震；（b）2018 年12 月26 日16:16:18 滁州定远地震；（c）2017 年11 月1 日23:43:31 常州溧阳地震Fig.5 Waveforms of three earthquakes

3.2 康拉德面特征

6 个地震的Sc 波到时清晰，根据表1所列走时和震中距，设Δ2=y，t2=x，v2=a，-（2H-h）2=b，采用式（4），利用最小二乘法拟合，得到b=-1 272.6、v=3.59 km/s，拟合结果见图6，可得拟合关系如下

图6 6 个震例的Sc 震相t2—Δ2 拟合图Fig.6 The t2—Δ2 fitting diagram of the Sc seismic phase of 6 earthquakes

所得波速v=3.59 km/s 与谢石文等（2016）所得安徽地区地壳中间面上层的横波速度3.57 km/s 基本一致。根据b值和h值，计算得到每个震例对应的康拉德面深度HC：20.35 km、21.85 km、22.85 km、20.85 km、21.85 km、23.85 km（表1），与滕吉文等（1985）测算的马鞍山—常熟间康拉德面深度（HC=22.38 km）距离接近，与谢石文等（2016）计算的安徽地区地壳中间层深度（HC=19 km）亦相差不大。由图6 可见，6 个地震拟合效果较好，数据点均位于拟合直线上或附近，说明6 个Sc 波均为来自地壳同一个层面的反射波。

4 结论

通过对马鞍山测震台网记录的近震波形进行分析，可以得出以下结论。

（1）地方震震中距在达到70—80 km 时，常会出现振幅较大的莫霍面反射波。但是，当地方震震中距在40 km 上下时，也可能出现振幅较大的康拉德面反射波。

（2）就安徽地区而言，一般震源深度在10 km 左右，当震中距约小于41 km 时，来自莫霍面的反射纵波PmP 比直达横波到时晚，其震相容易淹没在Sg 震相中，大于该震中距时，PmP 震相将在Sg 震相前出现。对于存在上地壳的地区，当震中距小于19 km 时，康拉德面反射波Pc 震相比直达横波Sg 震相到时晚，大于该震中距时，则Pc 波到时早。当震中距大于干扰距离时，Sg震相前会出现一个纵波反射波，其振幅比Pg震相大，但比Sg震相小。

（3）通过对6 个地震波形记录中初动较清晰的Sc 波走时进行拟合，得出：地壳中间面上横波速度为3.59 km/s，康拉德面深度HC在20.35—23.85 km 之间。由2016 年2 月13日芜湖市弋江区ML2.1 地震记录，得出HC=21.85 km，证明该不明震相（横波反射波）为康拉德面反射波。

研究分析过程得到审稿专家的悉心指导，并得到安徽省地震局隆爱军工程师的帮助，在此表示衷心的感谢。