0 引言

微地动(Microseisms)是地球背景噪声的主要组成部分(Faulkner et al，2018)，普遍存在于沿海地区宽频带地震仪信号中，信号的位移振幅一般为微米量级。20世纪初有学者观测到微地动现象(Wiechert，1904)，并且发现微地动信号为波浪所激发(Gutenberg，1936；Bernard，1941)，其源主要分布在海洋和湖泊(Koper et al，2009)。Longuet-Higgins(1950)首次提出了微地动产生机制的理论，Hasselmann(1963)对该理论进行了补充，微地动信号以表面波和体波形式进行传播(Lacoss et al，1969；Hasselmann，1963)，信号周期为0.5~20s(McNamara，2004；Gal et al，2019)。根据周期范围和激发机理的不同，微地动通常可以分为两类(Ardhuin et al，2015)：单频微地动(Single Frequency Microseisms，SF)，也称为第一类微地动(Primary Microseisms，PM)，周期为10~20s(Gal et al，2019)，其产生机制为海浪运动产生的压力传导至海底，与海底地形耦合并转换为地震波，与海浪具有相同的周期，因为海浪产生的压力随海水深度呈现指数型递减，因此单频微地动的源区主要分布在浅水区域(Ardhuin，2018)；双频微地动(Double Frequency Microseisms，DF)，又称为第二类微地动(Secondary Microseisms，SM)，其周期范围为0.5~10s(Gal et al，2019)，双频微地动的信号特征明显，产生机理复杂，目前主流观点认为其产生机制为两列传播方向相反、周期相近的海浪相互作用产生驻波，驻波产生二阶压力扰动并传导至海底转换为地震波，其周期为海浪周期的一半(Longuet-Higgins，1950；Gualtieri et al，2013)，水体中的共振效应在双频微地动的激发过程中起了重要的作用(Gualtieri et al，2013)。Kedar等(2008)首次对Longuet-Higgins-Hasselmann理论使用数值模型进行了验证，Ardhuin等(2011)在Kedar等(2008)工作的基础上，对数值模型进行了扩展，提出全球地震背景噪声能量辐射模型(Global Model of Seismic Noise Energy Radiation)，该模型能够模拟计算全球任意位置由海浪所致的双频微地动合成理论功率谱，并加入了海岸线反射波浪的影响。由于双频微地动产生于海浪间的相互作用，因此其既可以产生于远离海岸的深海区，也可以产生于近岸海域(Bromirski et al，2005)。很多学者通过单台站法(Koper et al，2015；Schimmel et al，2011；Bromirski et al，2002；Stutzmann et al，2012)和台阵法(Farra et al，2016；Webb，1998；Euler et al，2014；Liu et al，2016)研究微地动源区位置的空间分布和时间变化。

台风作为一种强对流天气，会引起海面波浪状态发生改变，并激发出微地动信号，国内很多学者在宽频带地震仪和重力仪中观测到了台风激发的微地动信号，且近海附近的信号比内陆地区更为明显(张雁滨等，2010、2013；王梅等，2014、2018；蒋骏等，2012)。国内多位学者分析了微地动信号特征、源区分布及产生机制等，如Xiao等(2018)通过陆地地震仪和海底地震仪分析了中国南海微地动的源区分布及变化；Lin等(2017、2018)分析了台风“Megi”期间不同频段微地动源区的位置分布，根据双频微地动在0.1~0.2Hz频段间的频率弥散特征，分析了海面涌浪的分布和传播过程。

本文使用日本和中国台湾地区的宽频带地震仪数据，通过时频分析和数值模拟方法，分析台风“尼伯特”激发的微地动信号的频谱变化特征及其产生原因，有助于加深我们对台风激发微地动信号的了解。

1 数据与方法 1.1 台风数据

台风“尼伯特”为2016年第1号台风(图 1)，本文使用的台风路径、强度等信息来自日本国立情报研究所^①。2016年7月2日21时左右，“尼伯特”在西北太平洋洋面上生成，7月3日凌晨，美国联合台风预警中心(Joint Typhoon Warning Center，JTWC)将其归为热带低压，之后向西北方向移动，当日9时许，日本气象厅(Japan Meteorological Agency，JMA)将其分类为热带风暴，最大风速18m/s；7月4日18时，“尼伯特”增强为强热带风暴，最大风速达到25m/s；7月5日6时，“尼伯特”进入菲律宾海，最大风速发展至32m/s，JTWC和JMA将“尼伯特”升级为台风级别；至7月6日0时左右，尼伯特发展为超强台风，最大风速47m/s，当日下午，最大风速达到57m/s，台风中心气压900hPa，此后“尼伯特”继续向西北方向移动；7月7日21时左右，台风“尼伯特”在中国台湾台东县登陆，强度开始减弱，7月9日6时左右，“尼伯特”再次在中国福建省石狮市登陆，当日中午减弱为热带低压后消散。

① http://www.agora.ex.nii.ac.jp/digital-typhoon

1.2 地震仪波形数据

本文使用的宽频带地震仪数据来自中国台湾地震工程研究中心和日本F-NET网站^②，原始数据采样频率100Hz。首先对台风期间地震仪垂向分量数据进行预处理，包括去仪器响应、去线性趋势、重新采样至1Hz、去突跳信号。图 1(b)~(d)分别为UMJ、KGM、B021台站7月2日0时至7月23日0时的波形图，由图可以发现在台风期间，各台站波形均出现了不同强度的幅值增强，包络呈现纺锤状。对于B021台站和KGM台站，在7月5日0时左右出现了明显的波形幅值增强，对于北部离台风距离相对较远的UMJ台站，在7月6日12时左右出现明显波形幅值增强。

② http://www.fnet.bosai.go.jp/top.php

1.3 波形数据时频特征分析

为分析台“风尼伯”特激发微地动的频谱特征，对台风期间地震仪波形数据进行短时傅里叶变换(Short-Time Fourier Transform，STFT)，数据窗口长度2048s，步长1024s，采样频率1Hz，窗函数为汉明窗。图 2(a)、2(c)、2(e)分别为B021、KGM和UMJ台站7月3日12时至7月9日0时短时傅里叶变换所得的功率谱。

2 数据处理结果

分析图 2(a)、2(c)、2(e)发现，在3个台站的功率谱中，均记录到了由于台风所导致的0.1~0.45Hz频率区间的PSD增强，但不同时刻PSD增强的频率段有所不同。对于B021台站(图 2(a))，在7月5日5时前后，0.4Hz频率左右出现了明显的PSD增强现象，然后由高频转至低频，7月7日0时PSD增强频率变化至约0.2Hz，为分析0.2~0.4Hz频段的频率变化特点，我们绘制了7月3日12时至7月9日0时的功率谱等值线图，如图 2(b)所示，由图可见，自7月5日5时至7月7日0时左右，PSD增强频段的低频值逐渐变小，如图 2(b)中黑色虚线所示。

对于KGM台站(图 2(c))，7月5日0时约0.4Hz频率附近出现PSD增强，7月6日19时PSD增强频率变化至约0.2Hz，其频率变化如图 2(d)中黑色虚线所示。

对于位于相对北部的UMJ台站，其功率谱(图 2(e))和功率谱等值线(图 2(f))中未发现类似B021和KGM台站在约7月5日0时至7月7日0时之间的频率变化特征。

我们分析了中国台湾和日本共117个宽频带地震仪垂向分量波形数据的功率谱，发现在约7月5日至7日之间，存在0.2~0.4Hz频段PSD增强及频率由约0.4Hz向0.2Hz变化特点的台站共33个，主要分布在中国台湾、琉球群岛及屋久岛，如图 1(a)中黑色正三角形所示。日本屋久岛以北地区的台站，在约7月5日至7日之间0.2~0.4Hz频段双频微地动信号不明显，频谱中0.2~0.4Hz无明显PSD增强，台站分布如图 1(a)中黑色空心倒三角形所示。

为了分析约7月5日至7日之间PSD增强频率由高频(约0.4Hz)转至低频(约0.2Hz)现象在不同台站出现的时间分布特征，计算每个时刻对应的0.3~0.4Hz频段PSD值之和(图 3(a))，并得到PSD随时间变化曲线(图 3(b))，选择7月4日12时开始(图 3(b)中竖直虚线)从左至右第一个局部最大值点，作为台风初始影响点(图 3(b)中黑色实心圆点)，记录台风初始影响点所在时间，用该时间作为当前台站频率由高频(约0.4Hz)转至低频(约0.2Hz)现象出现的时间。我们选择了干扰信号较少的14个台站(图 1(a)中黑色实心正三角形)，分别提取了每个台站的台风初始影响点，结果如图 4所示。由图 4可以发现，对大部分台站，当台风与台站距离在1500~2000km之间时，地震仪垂向分量信号出现PSD增强频率由高频(约0.4Hz)转至低频(约0.2Hz)的现象。

3 双频微地动数值模拟

在台风“尼伯特”的影响下，约7月5日至7日之间，双频微地动信号在0.2~0.4Hz频段PSD增强且频率由高频(约0.4Hz)转至低频(约0.2Hz)，单频微地动信号无明显变化，因此我们对台风“尼伯特”所激发的双频微地动信号进行数值模拟，以探究双频微地动信号产生上述频率变化的可能原因。

我们从法国海洋开发研究院(Institut Français de Recherche Pour L'exploitation de La Mer，IFREMER)下载了由频率相近、传播方向相反的海浪之间相互作用产生的二阶压力的功率谱数据(Frequency Spectrum of the Second Order Pressure，P2L)，P2L为第三代海浪模式(WAVEWATCH Ⅲ)通过输入风场等数据模拟所得(Tolman et al，2016)，空间分辨率0.5°，时间分辨率3h。将P2L数据输入全球地震背景噪声能量辐射模型(Global Model of Seismic Noise Energy Radiation)，模拟得到KGM台站所在位置的双频微地动合成理论功率谱，如图 5所示，其中图 5(a)为KGM台站垂向分量波形数据计算所得功率谱，波形数据为位移型数据；图 5(b)为模型中有海岸线反射效应时，模拟所得的合成理论功率谱；图 5(c)为模型中无海岸线反射效应时，模拟所得合成理论功率谱。

4 讨论与结论 4.1 讨论

通过上述分析可知，对于B021和KGM台站，约7月5日至7日之间双频微地动信号明显，0.2~0.4Hz频段PSD增强，且自约7月5日至7日，PSD增强频率由约0.4Hz转至约0.2Hz(图 2(b)、2(d))，出现这种频率变化特点的台站主要分布在中国台湾、琉球群岛和屋久岛地区，屋久岛以北的地震台站0.2~0.4Hz频率变化现象不明显。通过提取台风初始影响点，我们分析了14个地震台站出现频率变化现象的时间，发现当发生这种频率变化时，台风距离台站的距离集中在1500~2000km之间。

双频微地动的产生机制为频率相近、传播方向相反的两列波浪相互作用产生驻波，驻波产生二阶压力作用于海底所致。Ardhuin等(2011)提出产生双频微地动的海况，其中在单个风暴作用下产生双频微地动的海况有2种：①风暴快速移动，不同时刻产生的频率相近、传播方向相反的波浪相互作用激发微地动；②波浪向岸传播并被海岸线反射，反射的离岸波浪与向岸传播波浪发生相互作用，激发微地动。

为了探究约7月5日至7日之间出现的0.2~0.4Hz频段PSD增强及频率由约0.4Hz转至约0.2Hz现象的产生原因，我们使用全球地震背景噪声能量辐射模型，对KGM台站所在位置的双频微地动信号进行了数值模拟(图 5)。对比图 5(a)、5(b)可以发现，有海岸线反射效应的合成功率谱与实际观测数据的功率谱比较相近。在有海岸线反射效应的合成功率谱中(图 5(b))，在约7月5日至7日之间，可以观察到与图 5(a)相似的频率变化特征；在无海岸线反射效应的合成功率谱中(图 5(c))，在约7月5日至7日之间0.2~0.4Hz频段PSD增强及频率由约0.4Hz转至0.2Hz现象不明显，因此所观测到的约7月5日至7日之间0.2~0.4Hz频段频率变化现象，主要为海岸线反射效应所致，且源区可能主要分布在近岸附近。

4.2 结论

在台风“尼伯特”期间，当台风距离中国台湾、琉球群岛及屋久岛地区的宽频带地震台站1500~2000km时，在地震仪垂向分量信号功率谱中，记录到0.2~0.4Hz频段功率谱密度增强且增强频率由约0.4Hz转至0.2Hz的现象，此现象主要为海岸线反射效应所致，即由台风产生的向岸传播的波浪被海岸线反射，反射波浪与频率相近、传播方向相反的向岸传播波浪发生相互作用，激发出0.2~0.4Hz频段的双频微地动信号。