海南岛及邻区地壳S波Q值成像特征研究

王惠琳　徐晓枫　周龙泉　张晓东沈繁銮　李志雄　占　伟

1)海南省地震局，海口　570203 2)中国地震台网中心，北京　100045 3)中国地震局地震预测研究所，北京　100036 4)中国地震局第一监测中心，天津　300180



海南岛及邻区地壳S波Q值成像特征研究

王惠琳1)徐晓枫1)周龙泉2)张晓东3)沈繁銮1)李志雄1)占伟4)

1)海南省地震局，海口570203 2)中国地震台网中心，北京100045 3)中国地震局地震预测研究所，北京100036 4)中国地震局第一监测中心，天津300180

利用海南岛及邻区2000—2012年精定位后的地震S波资料进行Qs值成像，得到区域Qs值分布和滑动平均Qs值变化，结合地震活动、火山、热流点和GPS基线变化分析区域浅层地壳(深5～15km)S波衰减的时空变化特征，并讨论了区域地壳介质应力条件状态变化的特点。结果表明： 1)衰减成像结果显示研究区Qs值空间上横向不均匀变化显著，总体呈现海南岛中部高、周边低的特征，低Qs值(高衰减)主要分布在王五-文教断裂以北、尖峰-吊罗断裂以南及铺前-清澜断裂以东区域。地震活动多丛集分布在高、低Qs值过渡区域，火山主要分布在王五-文教断裂以北的雷琼坳陷低Qs值区域及其附近，大地热流点分布于低Qs值及其边缘区域。2)滑动平均Qs值的时间变化特征与地震活动性、泸州—琼中GPS基线变化存在较强的相关性。2000—2008年，研究区地震活动性较强，泸州—琼中GPS基线呈缩短趋势，滑动平均Qs值相对较高，区域可能处于 “高”应力介质条件状态；2008—2012年，研究区地震活动性较弱，泸州—琼中GPS基线显示伸长趋势，滑动平均Qs值相对较低，区域可能处于应力调整、松弛的 “低”应力介质条件状态。

Qs成像地震活动性火山热流GPS基线海南岛及邻区

0　引言

地震波在地球介质中传播时产生的能量衰减反映了传播介质的非弹性、非均匀性质(Winkleretal.，1982；陈颙等，2009)。Q值是用来描述地球介质衰减特性的重要参数，其与地壳的结构特征，如构造活动区的断层分布、壳内低速层、火山区的岩浆分布、热流分布等密切相关(苏有锦等，2006；赖晓玲等，2006；马宏生等，2007；周龙泉等，2009)。因地壳应力状态发生改变使得岩石性质产生变化，例如裂隙状态、流体运移、热活动状态等，均可能会导致Q值发生变化(周连庆等，2009；王惠琳，2011)。近年来，地震波衰减层析成像技术被有效用于获取地震波衰减的精细分布特征，如对构造区、地震带、火山裂谷带、陨石坑及水库等地区的衰减结构进行了不同程度的研究，取得了很好的成果(Haydaretal.，1990；洪学海等，2003；Hansenetal.，2004；汪素云等，2007；Zhao，2010；樊计昌等，2010；王惠琳等，2012a，b)。

本文基于2000—2012年海南岛及邻区双差精定位后(徐晓枫等，2014)的ML≥1.0地震S波形资料，利用衰减层析成像方法反演区域地壳Qs值分布和滑动平均Qs值，结合地震活动、火山、热流分布和GPS基线等进行分析，探讨区域地壳S波衰减时空变化特征及大震前后的地壳介质应力状态变化特征。

1　区域断裂构造

海南岛属于新华夏系第2隆起带的南延部分，由于受到南海弧后盆地的多期扩张活动影响，使得雷琼地区分布着多组EW—NEE走向的隆起、坳陷区和右旋扭动断裂(陈运平等，2002)。目前，海南岛及邻区活动断裂大致分布有以岛陆区域为主的近EW向深大断裂、NW向断裂和以海域、近海区域为主的NE向断裂(图1)。

图1　海南岛与邻区台网及断裂分布Fig. 1　Distribution of faults and seismic stations in Hainan Island and adjacent regions.F1 王五-文教断裂；F2 尖峰-吊罗断裂；F3 九所-陵水断裂；F4 铺前-清澜断裂；F5 海口-云龙断裂；F6 长流-仙沟断裂；F7 马袅-铺前断裂

近EW向断裂为贯穿全岛的张性深大断裂，从北向南依次为王五-文教断裂(F1)、尖峰-吊罗断裂(F2)和九所-陵水断裂(F3)。王五-文教断裂(F1)横贯海南岛北部，断裂切割深度达地幔顶部，断裂两端延入海域，是雷琼坳陷和五指山隆起的构造分界，也是新、老地层的分界线(谢振福等，2010)。断裂北部及雷州地区，第四纪火山发育、地震活动强烈，断裂南至尖峰-吊罗断裂(F2)，地震相对平静，晚更新世晚期以来活动不明显(张军龙等，2008)。尖峰-吊罗断裂(F2)新生代以来活动不明显，沿断裂温泉成群出露。九所-陵水断裂(F3)为早-中更新世断裂，在断裂展布区，热泉遍布，海蚀地貌发育(张业明等，2005)。

NW向断裂分布于琼北地区，具有新生性和强烈的活动性，从东向西依次为铺前-清澜断裂(F4)、海口-云龙断裂(F5)、长流-仙沟断裂(F6)(李志雄等，2008)。铺前-清澜断裂(F4)为晚更新世断裂，是琼东北隆起区的主要断裂，控制着东寨港凹陷的生成和发展，构成琼北地区第四纪玄武岩东侧边界。海口-云龙断裂(F5)呈NNW向展布在琼北地区，断裂西侧为起伏丘陵区，东侧为南渡江下游水网沼泽地及滨海平原。长流-仙沟断裂(F6)为北部湾盆地福山凹陷东部边界断裂，控制了西侧福山凹陷厚达数千米的古近系、新近系沉积。

2　方法原理

地震波传播路径造成的衰减中，由于地壳介质的非弹性性质而引起的衰减可用Q值来表示。台站观测到的数字地震波形振幅谱用方程式表示为(Sherbaum，1990)

(1)

式(1)为第j个台站记录到的第i个地震事件的振幅谱，包括地震波传播时所携带的震源谱(Si(f))、仪器响应(Ij(f))、场地响应(Rj(f))、几何扩散(Gij(r))、衰减谱(Bij(f))等综合信息，f为频率。

将震源谱、几何扩散分别采用Brune类型的ω2震源模式(Brune，1970)和地震波在全空间传播时的理论几何衰减函数(Atkinsonetal.，1992)表示，并代入沿整个射线路径的衰减谱。 由于台站建在基岩上，故假定台站场地响应是接近1的常数(周龙泉等，2009)，扣除仪器响应后，式(1)写为

(2)

(3)

2.1二维Q值成像

本文采用二维Q值成像反演，假定地壳介质为S波速度均匀的2维介质，射线路径为直线，则式(3)中的t*扰动主要由沿射线路径的品质因子Q值的扰动而引起，可写为式(4)：

(4)

根据式(4)，假定Q值不变，则t*与震源至台站的走时成正比，因此震源深度对t*数据具有较大的影响。利用t*反演Q值的过程中，必须消除震源深度的影响。在计算理论走时时，震中距用震源至台站的距离来代替。

2.2数据预处理

本文收集到海南岛及邻区2000—2012年共36个数字地震台记录到的地震波形资料，其中包括海南地震台网22个台站、广东地震台网8个台站、广西地震台网6个台站(图1)。台站覆盖了雷琼地区，陆地台站间距平均50km左右。地震总数为1,035个，震级范围为ML1.0～4.5，震源深度0.5～29km，其中分布在5～15km深度的地震占83%，故成像结果主要反映了研究区域浅层地壳(深5～15km)的S波衰减特征。

对于每一条地震波形记录，为得到剪切波频谱的较少畸变的总体形状，采用包含所有可识别S震相的 “S窗”内的信号进行分析，“S窗”是指从S波震相开始到包含其总能量90%的时间段(图2a)。由于地震震中距范围为1～500km，其中96%的地震震中距<300km，故 “S窗”主要包括直达S波及地壳内部间断面、莫霍面的反射波(陈继锋等，2010)。截取后的 “S窗”，为了获得稳定的观测信号傅立叶谱，采用平移窗谱方法(刘杰等，2003；黄玉龙等，2003)。将 “S窗”波形信号分成多个小段，每个小段分别为256个采样点，相邻段间均重叠50%(图3)。挑选至少有3个台观测到且每个台须记录到至少3个地震事件(郭晓等，2010)、信噪比SNR≥2的地震波形用于反演。

图2　噪声窗和 “S窗”(a)与地震在不同台站的观测位移谱和拟合谱(b)Fig. 2　Windows of noise and “S” wave(a) and Observation and fitness spectra of earthquake at different station(b).WET 翁田台；SHP 沙坡台；SAJ 三江台；HSK 火山口台；QXL 七星岭台

图3　平移窗谱技术图(王惠琳，2011)Fig. 3　Translation of window spectrum chart(after Wang hui-lin，2011).

本文对海南岛及邻区地震记录1～15Hz的S波谱，将时间域的波形数据转换到频率域，并将速度振幅谱转换为位移振幅谱。为了降低解的非唯一性和得到全局最优解，采用多台记录谱联合反演方法和遗传算法反演地震在不同台站的位移谱和拟合谱，不同台站位移谱衰减通过消除噪声、几何扩散的影响后，得到不同路径上Qs值的影响。图2a给出了2012年12月28日文昌3.0级地震在翁田台记录到的S波水平分量波形及选取的噪声窗和 “S窗”，图2b为该地震在不同台站的观测位移谱和拟合谱。

3　海南岛及邻区S波Q值成像

3.1检测板分辨率测试

图4　海南岛及邻区台站、事件及射线分布Fig. 4　Distribution of stations，earthquakes and raypaths in Hainan Island and adjacent regions.▲台站；+事件；/射线

图5　检测板分辨率试验结果Fig. 5　Checkerboard resolution tests.矩形为研究区

图6　研究区Qs成像结果Fig. 6　Distribution of Qs in the study area.黑色空心圆： 地震；白色虚线椭圆： 地震丛；黄色实心圆： 热流；红三角： 火山；红五角： 1605年琼山 级地震；F1 王五-文教断裂；F2 尖峰-吊罗断裂；F3 九所-陵水断裂；F4 铺前-清澜断裂；F5 海口-云龙断裂；F6 长流-仙沟断裂；F7马袅-铺前断裂

3.2研究区地壳S波Q值成像结果

图6 为海南岛及邻区地壳S波Q值成像结果，反映了研究区区域地壳空间上的地震波衰减分布特征。

从Qs值分布看，研究区地壳Qs值横向不均匀变化显著，总体呈现海南岛中部高、周边低的特征(图6)。海南岛中部隆起、北部湾盆地中部区域显示为高Qs值(低衰减)特征，而低Qs值(高衰减)区域则主要分布在王五-文教断裂(F1)以北(雷琼断陷区)、尖峰-吊罗断裂(F2)以南和铺前-清澜断裂(F4)以东区域。王五-文教断裂(F1)以北、尖峰-吊罗断裂(F2)以南区域断层汇聚、交错、相互切割，断裂区域岩石破碎程度高，节理、裂隙发育，能量容易被吸收，地震波衰减大，显示为低Qs值异常(樊计昌等，2001)。此外，海南岛速度成像结果显示，在浅部，王五-文教断裂(F1)以北延伸至雷州半岛地区为低波速异常，断裂以南地区为明显的高波速异常(雷建设等，2008；李志雄等，2008；Leietal.，2009)。这与本文中王五-文教断裂以北的雷琼断陷显示低Qs值、断裂以南的海南岛中部隆起显示高Qs值结果相对应，速度和衰减成像结果均一致反映了相同的地壳介质特性。

从雷琼地区火山和大地热流分布看(汪集旸等，1990；胡圣标等，2001；HUetal.，2007；李志雄等，2008)(图6)，雷琼坳陷区内的雷州、琼北2大火山群(40余座火山)和大地热流点(实测热流值范围为61.5～82.0mW/m2，平均热流值为70.0mW/m2，而Davies(2013)计算的全球大陆平均热流值为65mW/m2)基本分布在低Qs值区域或其边缘地带，主要分布在广东湛江、雷州地区和海南海口、临高、万宁地区。刘建华等(2004)、周龙泉等(2009)研究发现华北、云南地区热流点基本与低Q值区域相对应，即热流值与Qs值呈负对应关系。大地热流是地球内部热作用过程在地表最为直接的显示，反映了区域地热状态、地质构造、深部地壳结构的性质。研究区低Qs值可能与地壳内部热作用过程相关，且可能与新生代以来一直延续至第四纪晚期的多期火山岩浆喷发、地幔对地壳底部的底侵、拆沉和地壳仍处于相对温热状态密切相关(陈墨香等，1991；嘉世旭等，2006)。

3.3研究区滑动平均Qs值变化特征

如图7 所示，2000—2008年，研究区各时段平均Qs值逐步升高，在2005—2008时段Qs值达到高值，相应地，地震个数2000—2008年期间处在相对高值阶段。2005—2008年时段后，Qs值开始转折下降，相应地，地震个数也处于下降的低值阶段。裴顺平等(2010)以1a为步长、2a为窗长得到的川滇地区横波平均Q值动态变化与M≥4.5地震个数的关系研究发现，滑动平均Q值与地震活动性存在非常强的正相关关系，即地震个数越多，区域平均Q值越大。本研究结果同样显示了平均Qs值与区域地震个数较强的正相关关系，但本文研究区平均Qs值较地震个数在时间上存在一定的滞后性。川滇地区是中国地震活动较强的地区，地壳运动相对剧烈，其地震频次较高、强度较大；而本文研究区地震频次较低、强度较弱。虽然不同研究区域存在不同的地震活动特点，采用不同的时间窗长和地震震级，但研究结果均表现为平均Qs值与地震个数呈正相关关系。

图7　滑动平均Qs值与ML≥3.0地震的个数对比图Fig. 7　Contrast figure of moving average Qs values with ML≥3.0 earthquake number.

由于每个时段窗长为4a，获得的地震个数为4a的平均结果，难以直观显示地震发生的时间进程与地下介质变化特征的对应情况，故以地震时间序列与Qs值进行对比。如图8 所示，2000—2008年，各时段平均Qs值逐步升高，相应地，ML≥3.0以上地震频发，每年平均发生ML≥3.0地震6次，地震强度逐年升高；2005—2008时段后，平均Qs值达到高值后转折下降，同时，自2008年7月10日北部湾ML3.2地震后，地震数量急剧减少，每年平均仅发生ML≥3.0地震1次，且出现2008年8月—2012年4月共3年8个月的时间里仅发生1次ML3.5地震的异常平静状态。由此可见，平均Qs值升高的时段，地震频率高，震级逐步升高；平均Qs值下降的时段，地震频率低，震级逐步下降。表明研究区滑动平均Qs值的升高—下降趋势恰与地震活动性的增强—减弱趋势相对应。

图8　滑动平均Qs值与ML≥3.0地震时间序列对比图Fig. 8　Contrast figure of moving average Qs values with ML≥3.0 earthquake time series.

对于平均Qs值与地震活动性的正相关关系，裴顺平等(2010)给出了解释，认为由于应力水平的增加，导致岩石中裂隙闭合、裂纹密度减小、流体含量减少，地震波衰减随之减小，表现为Q值较高，而应力积累的外部表现就是地震活动性的增强。因此可认为，2000—2008年研究区平均Qs值相对处于高值，地震活动性较强，表明研究区区域整体应力水平较高；2008—2012年研究区平均Qs值处于持续下降趋势，地震活动性较弱，表明研究区区域整体应力水平较低。

4　滑动平均Qs值与GPS基线关系的讨论

图7、8 所示的研究区平均Qs值升高—下降与地震活动性的增强—减弱的时间转折点大致在2008年前后。众所周知，2008年5月12日，华南地块和青藏高原的接合部——龙门山逆冲推覆断裂上发生了汶川8.0级特大地震，造成的地表破裂长度达(240±5)km，最大垂直位移为(6.5±0.5)m，最大右旋走滑位移4.9m，且在龙门山推覆构造带中段产生了最大约7m的地壳缩短量(徐锡伟等，2010)。研究区2000—2008年的相对高Qs值、强地震活动性和2008—2012年的相对低Qs值、弱地震活动性是否可能受大地震前、后区域应力状态改变的影响？

由于大地震的发生是应力不断积累到释放的过程，空间技术能够观测到地球内部应力场变化在地面的形变响应分布及构造应力缓慢作用过程中地壳中的某些形变和运动特征(张永志等，2000)。程万正等(2012)研究发现大地震造成的地壳大尺度变形呈现为块体间的汇聚挤压，地壳形变呈现为缩短；大地震后则呈现为调整、松弛或拉张状态，地壳形变呈现为伸长。GPS基准站间的基线变化能很好地反映区域地壳相对运动的变化，特别是地块的刚性运动呈现的挤压—拉张状态。同一地块内2个GPS站间的基线变化可反映地块内部沿基线方向的变形(压缩或伸张)动态。为此，本文选取中国地壳运动观测网络基准站中位于华南地块西部的四川泸州(LUZH)站和研究区中部的琼中(QION)站，计算站间的GPS基线变化，并与研究区的滑动平均Qs值进行对比分析。

图9　滑动平均Qs值、泸州—琼中GPS基线对比图(a)与泸州、琼中GPS站位置示意图(b)Fig. 9　Contrast figure of moving average Qs values with Luzhou-Qiongzhong GPS baseline(a)and Luzhou and Qiongzhong GPS stations location map(b).黑点线：平均 Qs值；红五角： 汶川8.0级地震；红三角： GPS站；黄色线与红色线分别为GPS基线的月均值与年均值

如图9 所示，2000年1月—2007年9月琼中—泸州基线年动态特征呈长期线性缩短趋势，沿基线方向地壳运动表现为挤压、应力增强，相应地2000—2008年各时段平均Qs值逐步升高，表明区域处于 “高”应力的介质条件状态；2007年10月后琼中—泸州基线出现转折、反向伸长，沿基线方向地壳运动表现为拉张、应力松弛状态，相应地2008—2012年各时段平均Qs值呈下降趋势，表明区域处于应力调整、松弛的 “低”应力介质条件状态。汶川地震的孕震空间和时间尺度相对较大，且其震源破裂模型为逆冲兼右旋走滑型(王凯英等，2011)，在其能量积累、释放过程中，由于应力在地壳内部的传递作用，对华南地块的地壳变形可能存在加载、卸载作用。而研究区位于华南地块西南端，滇、黔、桂构造带的东南延伸位置，区域整体应力可能随之发生相应变化，从而表现出2000—2008年地震活动性较强、Qs值较高的 “高”应力介质条件状态和2008—2012年地震活动性减弱、Qs值下降的 “低”应力介质条件状态。这可能为进一步理解大地震的孕育发生及其主要控制因素——应力环境的不断变化提供了证据。

5　结论与讨论

通过衰减层析成像方法反演研究区地壳Qs值分布和滑动平均Qs值，结合地震活动、火山、热流分布和泸州—琼中GPS基线变化，分析、讨论区域地壳S波衰减时空变化及区域应力状态特征。

(1)Qs成像结果显示，研究区地壳Qs值横向不均匀变化显著，总体呈现海南岛中部高、周边低的特征，反映了空间上地壳衰减结构存在强烈的横向不均匀特性。高Qs值(低衰减)分布在海南岛中部隆起至北部湾盆地中部区域，低Qs值(高衰减)主要分布在王五-文教断裂(F1)以北(雷琼断陷区)、尖峰-吊罗断裂(F2)以南及铺前-清澜断裂(F4)以东区域。地震多丛集分布在高、低Qs值(低、高衰减)交界或过渡区域，从北向南地震丛集分布在琼西北的北部湾盆地NE向断裂密集区，琼北的NW、EW向断裂交会区，琼南的尖峰-吊罗断裂(F2)和九所-陵水断裂(F3)两端，而中部五指山隆起至北部湾盆地的高Qs值区域地震活动相对较弱。

本文的衰减成像结果与速度成像结果(雷建设等，2008；李志雄等，2008)具有较好的一致性，共同反映了相同的地壳介质结构特性；即在浅部，王五-文教断裂(F1)以北延伸至雷州半岛地区为低波速、低Qs值异常区域，断裂以南地区为明显的高波速、高Qs值异常区域，这可能暗示了深部地壳和上地幔热物质上涌的动力学过程，且可能与“海南地幔柱”(Leietal.，2009，2013)有着密切的相关性。同时，雷琼坳陷区内的雷州、琼北2大火山群(40余座火山)和大地热流点基本分布在低Qs值区域或其边缘地带，同样表明研究区内的低Qs值与地壳内部热作用过程关系密切，可能与新生代以来一直延续至第四纪晚期的多期火山岩浆喷发、地幔对地壳底部的底侵、拆沉和地壳仍处于相对温热状态密切相关(陈墨香等，1991；嘉世旭等，2006)。

(2)滑动平均Qs值与地震活动性、泸州—琼中GPS基线变化存在较强的相关性。2000—2008年，研究区地震活动性较强，平均Qs值相对较高，泸州—琼中GPS基线呈缩短趋势，区域处于 “高”应力介质条件状态；2008—2012年，地震活动性较弱，平均Qs值呈下降趋势，泸州—琼中GPS基线呈伸长趋势，区域处于应力调整、松弛的 “低”应力介质条件状态。可能与大地震的孕育、发生所导致的区域整体应力环境随之改变存在一定的关系，为大地震的主要控制因素——应力环境的不断变化提供了可能的证据。

泸州—琼中GPS基线由缩短到反向伸长的转折变化并未与汶川地震同步发生，而是超前约8个月，可能与岩石圈板块并不是纯刚性的运动，而是接近于弹性体或黏弹性体的运动状态有关(McCluskyetal.，2000；占伟等，2013)。

(3)本文反演的Qs与频率无关，其初始平均Qs值数值上比与频率相关的反演结果偏高(马宏生，2007；周龙泉，2009；王惠琳，2012)。虽然与频率无关的Q值数值上有所增大，但并不影响Q值的分布，因此依然能够采用其来分析地壳介质的物性分布特征(Eberhart-Phillipsetal.，2002)。

致谢感谢中国地震局地震预测研究所江在森研究员、武艳强研究员，福建省地震局王紫燕工程师的指导与交流；感谢评审专家为本文提出的宝贵意见。

陈继锋，赵翠萍，杨立明. 2010. 甘肃地区S波非弹性衰减Q值研究 [J]. 地震，30: 125—130.

CHEN Ji-feng，ZHAO Cui-ping，YANG Li-ming. 2010. Q values of S wave inelastic attenuation in Gansu region [J]. Earthquake，30: 125—130(in Chinese).

陈墨香，夏斯高，杨淑贞. 1991. 雷州半岛局部地热异常及其形成机制 [J]. 地质科学，4: 369—383.

CHEN Mo-xiang，XIA Si-gao，YANG Shu-zhen. 1991. Local geothermal anomalies and their formation mechanisms on Leizhou Peninsula，South China [J]. Chinese Journal of Geology，4: 369—383(in Chinese).

陈运平，席道瑛，樊星. 2002. 用分形理论研究海南岛的活动断裂 [J]. 地震研究，25: 351—355.

CHEN Yun-ping，XI Dao-ying，FAN Xing. 2002. The study on the active faults in Hainan Island by way of fractal theory [J]. Journal of Seismological Research，25: 351—355(in Chinese).

陈颙，黄庭芳，刘恩儒. 2009. 岩石物理学 [M]. 合肥： 中国科学技术大学出版社. 496.

CHEN Yong，HUANG Ting-fang，LIU En-ru. 2009. Rock Physics [M]. China Science and Technology University Press，Hefei. 496(in Chinese).

程万正，阮祥，江在森. 2012. 苏门答腊巨大地震前后的区域应力状态 [J]. 四川地震，4: 1—9.

CHENG Wan-zheng，RUAN Xiang，JIANG Zai-sen. 2012. Regional stress state before and after Sumatra earthquake [J]. Earthquake Research in Sichuan，4: 1—9(in Chinese).

樊计昌，刘明军，赵成彬，等. 2010. 岫岩陨石坑三维Q值层析成像 [J]. 地球物理学报，53: 2367—2375.

FAN Ji-chang，LIU Ming-jun，ZHAO Cheng-bin，etal. 2010. Three-dimensionalQtomography for Xiuyan meteorite impact crater [J]. Chinese Journal of Geophysics，53: 2367—2375.

郭晓，张元生，何斌. 2010. 兰州市及其附近地区介质衰减特征研究 [J]. 西北地震学报，32: 244—247.

GUO Xiao，ZHANG Yuan-sheng，HE Bin. 2010. Study on the inelastic attenuation of the medium in Lanzhou City and adjacent regions [J]. Northwestern Seismological Journal，32: 244—247(in Chinese).

洪学海，朱介寿，曹家敏，等. 2003. 中国大陆地壳上地幔S波品质因子三维层析成像 [J]. 地球物理学报，46(5): 642— 651.

HONG Xue-hai，ZHU Jie-shou，CAO Jia-min，etal. 2003. Tomography of the 3-D S-wave quality factor of the crust and upper mantle in China [J]. Chinese J Geophys，46(5): 642— 651(in Chinese).

胡圣标，何丽娟，汪集旸. 2001. 中国大陆地区大地热流数据汇编(第三版)[J]. 地球物理学报，44: 611— 626.

HU Sheng-biao，HE Li-juan，WANG Ji-yang. 2001. Compilation of heat flow data in the China continental area(3rd edition)[J]. Chinese J Geophys，44: 611— 626(in Chinese).

黄海波，丘学林，夏少红. 2012. 海南岛地壳厚度与泊松比结构 [J]. 热带海洋学报，31: 65—70.

HUANG Hai-bo，QIU Xue-lin，XIA Shao-hong. 2012. Crustal structure and Poisson’s ratio beneath Hainan Island [J]. Journal of Tropical Oceanography，31: 65—70(in Chinese).

黄玉龙，郑斯华，刘杰，等. 2003. 广东地区地震动衰减和场地响应的研究 [J]. 地球物理学报，46(1): 54— 61.

HUANG Yu-long，ZHENG Si-hua，LIU Jie，etal. 2003. Attenuation of ground motion and site response in Guangdong region [J]. Chinese J Geophys，46(1): 54— 61(in Chinese).

嘉世旭，李志雄，徐朝繁，等. 2006. 雷琼坳陷地壳结构特征 [J]. 地球物理学报，49: 1385—1394.

JIA Shi-xu，LI Zhi-xiong，XU Zhao-fan，etal. 2006. Crustal structure features of the Leiqiong depression in Hainan Province [J]. Chinese J Geophys，49: 1385—1394(in Chinese).

赖晓玲，孙译. 2006. 伽师近场观测地震资料的三维Q值成像 [J]. 大地测量与地球动力学，26: 28—32.

LAI Xiao-ling，SUN Yi. 2006. 3-DQvalue tomography of near-field observations of Jiashi strong earthquake swarms [J]. Journal of Geodesy and Geodynamics，26: 28—32(in Chinese).

雷建设，赵大鹏，沈繁銮，等. 2008. 海南地幔柱新的成像结果 [J]. 国际地震动态.11: 31.

LEI Jian-she，ZHAO Da-peng，SHEN Fan-luan，etal. 2008. New imaging results of Hainan mantle plume [J]. Recent Developments in World Seismology，11: 31(in Chinese).

李志雄，雷建设，赵大鹏，等. 2008. 海南岛及邻区地壳三维P波速度结构 [J]. 地震学报，30: 441— 448.

LI Zhi-xiong，LEI Jian-she，ZHAO Da-peng，etal. 2008. 3-D P-wave velocity structure of the crust beneath Hainan Island and adjacent regions [J]. Acta Seismologica Sinica，30: 441— 448(in Chinese).

廖超明，武军郦. 2009. 北部湾现今地壳运动特性研究 [J]. 大地测量与地球动力学，28: 22—26.

LIAO Chao-ming，WU Jun-li. 2009. Study on characteristics of present-day crustal movement in Beibu Gulf area [J]. Journal of Geodesy and Geodynamics，28: 22—26(in Chinese).

刘建华，刘福田，阎晓蔚，等. 2004. 华北地区Lg尾波衰减研究： Lg尾波Q0地震成像 [J]. 地球物理学报，47: 1044—1052.

LIU Jian-hua，LIU Fu-tian，YAN Xiao-wen，etal. 2004. A study of Lg coda attenuation beneath North China: Seismic imaging of Lg coda Q0[J]. Chinese J Geophys，47: 1044—1052(in Chinese).

刘杰，郑斯华，黄玉龙. 2003. 利用遗传算法反演非弹性衰减系数、震源参数和场地响应 [J]. 地震学报，25(2): 211—218.

LIU Jie，ZHENG Si-hua，HUANG Yu-long. 2003. The inversion of nonelasticity coefficient，source parameters，site response using genetic algorithms [J]. Acta Seismologica Sinica，25(2): 211—218(in Chinese).

马宏生，汪素云，裴顺平，等. 2007. 川滇及周边地区地壳横波衰减的成像研究 [J]. 地球物理学报，50(2): 465— 471.

MA Hong-sheng，WANG Su-yun，PEI Shun-ping，etal. 2007. Q0tomography of S wave attenuation in Sichuan-Yunnan and adjacent regions [J]. Chinese J Geophys，50(2): 465— 471(in Chinese).

裴顺平，刘杰，马宏生，等. 2010. 川滇地区横波Q值动态变化 [J]. 地球物理学报，53: 1639—1652.

PEI Shun-ping，LIU Jie，MA Hong-sheng，etal. 2010. Dynamic variation of S-waveQvalue beneath Sichuan-Yunnan，China [J]. Chinese J Geophys，53: 1639—1652(in Chinese).

沈繁銮，杨马陵，李志雄，等. 2010. 华南强震动力分析 [J]. 华南地震，30: 1—5.

SHEN Fan-luan，YANG Ma-ling，LI Zhi-xiong，etal. 2010. Dynamic analysis on strong earthquakes in South China [J]. South China Journal of Seismology，30: 1—5(in Chinese).

苏有锦，刘杰，郑斯华，等. 2006. 云南地区S波非弹性衰减Q值研究 [J]. 地震学报，28(2): 206—212.

SU You-jin，LIU Jie，ZHENG Si-hua，etal. 2006.Qvalue of anelastic S-wave attenuation in Yunnan region [J]. Acta Seismologica Sinica，28(2): 206—212(in Chinese).

王惠琳. 2011. 紫坪铺水库区域地壳S波Q值层析成像研究 [D]. 北京: 中国地震局地震预测研究所. 1—70.

WANG Hui-lin. 2011. Crustal S-wave attenuation tomography in Zipingpu reservoir area [D]. Beijing: Institute of Earthquake Science，China Earthquake Administration. 1—70(in Chinese).

王惠琳，张晓东，周龙泉，等. 2012a. 紫坪铺水库区域地壳Qs成像及其与渗透关系研究 [J]. 地球物理学报，55: 526—537.

WANG Hui-lin，ZHANG Xiao-dong，ZHOU Long-quan，etal. 2012a. Study on the relationship between fluid infiltration andQstomography of the crust in Zipingpu reservoir area [J]. Chinese J Geophys，55: 526—537(in Chinese).

王惠琳，张晓东，周龙泉，等. 2012b. 紫坪铺水库区域地壳Qs动态变化及其与水库蓄水关系的研究 [J]. 地震学报，34: 676— 688.

WANG Hui-lin，ZHANG Xiao-dong，ZHOU Long-quan，etal. 2012b. Relationship between the reservoir storage and crustalQsdynamic variation in Zipingpu reservoir area [J]. Acta Seismologica Sinica，34: 676— 688(in Chinese).

汪集旸，黄少鹏. 1990. 中国大陆地区大地热流数据汇编(第二版)[J]. 地震地质，12: 351—366.

WANG Ji-yang，HUANG Shao-peng. 1990. Compilation of heat flow data in the China continental area(2nd Edition)[J]. Seismology and Geology，12: 351—366(in Chinese).

王凯英，马胜利. 2011. 2008年汶川地震的破裂模型及引起的余震活动分析 [J]. 地球物理学报，54: 1511—1516.

WANG Kai-ying，MA Sheng-li. 2011. Analysis of co-seismic rupture models for 2008 Wenchuan earthquake and relations with aftershock [J]. Chinese J Geophys，54: 1511—1516(in Chinese).

汪素云，裴顺平，许忠淮，等. 2007. 利用Ml振幅研究地壳横波Q值I: 不同构造区的衰减特征 [J]. 地球物理学报，50(6): 1740—1747.

WANG Su-yun，PEI Shun-ping，XU Zhong-huai，etal. 2007. Crustal S-waveQestimated fromMLamplitude I: Attenuation in diferent tectonic regions of China [J]. Chinese J Geophys，50(6)： 1740—1747(in Chinese).

谢振福，谢顺胜，蔡水库，等. 2010. 琼北地区重磁场特征与区域地质构造 [J]. 物探与化探，34: 579—582.

XIE Zheng-fu，XIE Shun-sheng，CAI Shui-ku，etal. 2010. Gravity and magnetic field characteristics and regional geological structures of northern Hainan area [J]. Geophysical and Geochemical Exploration，34: 579—582(in Chinese).

徐德如，陈广浩，夏斌，等. 2003. 海南岛几个重大基础地质问题评述 [J]. 地质科技情报，22: 37— 44.

XU De-ru，CHEN Guang-hao，XIA Bin，etal. 2003. Comment on several important basic geological problems in Hainan Island，China [J]. Geological Science and Technology Information，22: 37— 44(in Chinese).

徐锡伟，陈桂华，于贵华，等.2010.5·12汶川地震地表破裂基本参数的再论证及其构造内涵分析 [J]. 地球物理学报，53: 2321—2336.

XU Xi-wei，CHEN Gui-hua，YU Gui-hua，etal. 2010. Reevaluation of surface rupture parameters of the 5·12 Wenchuan earthquake and its tectonic implication for Tibetan uplift [J]. Chinese J Geophys，53: 2321—2336(in Chinese).

徐晓枫，王惠琳，陈小敏. 2014. 海南岛及邻区地震精确定位及断裂构造分析 [J]. 地震研究，37: 216—221.

XU Xiao-feng，WANG Hui-lin，CHEN Xiao-min. 2014. Analysis on accurate location of earthquakes and fault structure in Hainan Island and its adjacent areas [J]. Journal of Seismological Research，37: 216—221(in Chinese).

徐义刚，魏静娴，邱华宁，等. 2012. 用火山岩制约南海的形成演化： 初步认识与研究设想 [J]. 科学通报，57(20): 1863—1878.

XU Yi-gang，WEI Jin-xian，QIU Hua-ning，etal. 2012. Opening and evolution of South China Sea constrained by studies on volcanic rocks: Preliminary results and aresearch design [J]. Chin Sci Bull，57(20): 1863—1878(in Chinese).

占伟，武艳强，丁晓光，等. 2013. 求解模型对欧拉矢量的影响 [J]. 中国地震，29: 91—96.

ZHAN Wei，WU Yan-qiang，DING Xiao-guang，etal. 2013. Effect of solving models on Euler vector [J]. Earthquake Research in China，29: 91—96(in Chinese).

张军龙，田勤俭，李峰，等. 2008. 海南岛北西部新构造特征及其演化研究 [J]. 地震，28: 85—94.

ZHANG Jun-long，TIAN Qin-jian，LI Feng，etal. 2008. Study on neotectonic characteristics and its evolution in northwestern Hainan Island [J]. Earthquake，28: 85—94(in Chinese).

张业明，谢才富，付太安，等. 2005. 海南岛地质构造演化刍论 [J]. 科学技术与工程，5: 1485—1487.

ZHANG Ye-ming，XIE Cai-fu，FU Tai-an，etal. 2005. Tectonic evolution of Hainan Island [J]. Science Technology and Engineering，5: 1485—1487(in Chinese).

张永志，崔笃信，王琪，等. 2000. 利用GPS资料研究区域地壳应力场变化与地震活动关系 [J]. 地震学报，22: 449— 456.

ZHANG Yong-zhi，CUI Du-xin，WANG Qi，etal. 2000. Relationship between regional stress field variation and earthquake activities from gps data [J]. Acta Seismologica Sinica，22: 449— 456(in Chinese).

张振克，孟红明，王万芳，等. 2008. 海南岛1605年历史地震的海岸沉积记录 [J]. 海洋地质与第四纪地质，28: 9—14.

ZHANG Zhen-ke，MENG Hong-ming，WANG Wan-fang，etal. 2008. Preliminary study on the coastal sediments records about the historical earthquake in the year of 1605，Hainan Island，China [J]. Marine Geology & Quaternary Geology，28: 9—14(in Chinese).

周连庆，赵翠萍，陈章立. 2009. 紫坪铺水库地区尾波Qc值研究 [J]. 地震，29(4): 44—51.

ZHOU Lian-qing，ZHAO Cui-ping，CHEN Zhang-li. 2009. Coda waveQc values in the Zipingpu reservoir area，Sichuan Province [J]. Earthquake，29(4): 44—51(in Chinese).

周龙泉，刘杰，苏有锦，等. 2009. 利用S波高频衰减参数对云南地区地壳Q值成像 [J]. 地球物理学报，52(6): 1500—1507.

ZHOU Long-quan，LIU Jie，SU You-jin，etal. 2009. Tomography forQof Yunnan region from high-frequency attenuation of S wave [J]. Chinese J Geophys，52(6)： 1500—1507(in Chinese).

Atkinson G M，Mereu R F. 1992. The shape of ground motion attenuation curves in southeastern Canada [J]. Bulletin of the Seismological Society of America，82: 2014—2031.

Brune J N. 1970. Tectonic stress and the spectra of seismic shear waves from earthquakes [J]. J Geophys Res，75: 4997—5009.

Cormier V F. 1982. The effect of attenuation on seismic body waves [J]. Bull Seism Soc Am，72: 169—200.

Davies J H. 2013. Global map of solid Earth surface heat flow [J]. Geochemistry Geophysics Geosystems，14(10)： 4068— 4022.

Eberhart-Phillips D，Chadwick M. 2002. Three-dimensional attenuation model of the shallow Hikuangi subduction zone in the Raukumara Peninsula，New Zealand [J]. J Geophys Res，107: 2033. doi: 10.1029/2000JB000046.

Hansen S，Thurber C，Mandernach M，etal. 2004. Seismic velocity and attenuation structure of the east rift zone and south flank of Kilauea volcano Hawaii [J]. Bulletin of the Seismological Society of America，94: 1430—1440.

Haydar J，AI-Shukri，Mitchell B J. 1990. Three dimensional attenuation structure in and around the new Madrid seismic zone [J]. Bull Seism Soc Ain，80: 615— 632.

HU Jiu-chang，BAI Deng-hai，WANG Wei-hua，etal. 2007. Deep electrical anomaly in theM7.5 Qiongzhou earthquake region and its relation with future seismicity [J]. Acta Seismologica Sinica，20: 273—279.

Lei J S，Furen Xie，Qicheng Fan，etal. 2013. Seismic imaging of the deep structure under the Chinese volcanoes: An overview [J]. Phys Earth Planet Interi，224： 104—123.

Lei J S，Zhao D. 2006. Global P-wave tomography: On the effect of various mantle and core phases [J]. Phys Earth Planet Interi，154(1): 44— 69.

Lei J S，Zhao D P，Steinberger B，etal. 2009. New seismic constraints on the upper mantle structure of the Hainan plume [J]. Phys Earth Planet Inter，173: 33—50.

McClusky S，Boassanian S，Barka A，etal. 2000. Global Positioning System constraints on plate kinematics and dynamics in the eastern Mediterranean and Caucasus [J]. J Geophys Res，15: 5695—5719.

Sherbaum F. 1990. Combined inversion for the three-dimensionalQstructure and source parameters using microearthquake spectra [J]. J Geophys Res，95: 12423—12438.

Winkler K W，Nur A. 1982. Seismic attenuation: Effects of pore fluids and frictional sliding [J]. Geophysics，47: 1—15.

Zhao L F，Xie X B，Wang W M，etal. 2010. Seismic Lg-waveQtomography in and around Northeast China [J]. J Geophys Res，115. B08307. doi: 10.1029/2009JB007157.

CRUSTALQsTOMOGRAPHY IN HAINAN ISLAND AND ADJACENT REGIONS

WANG Hui-lin1)XU Xiao-feng1)ZHOU Long-quan2)ZHANG Xiao-dong3)SHEN Fan-luan1)LI Zhi-xiong1)ZHAN Wei4)

1)EarthquakeAdministrationofHainanProvince，Haikou570203，China2)ChinaEarthquakeNetworksCenter，Beijing100045，China3)InstituteofEarthquakeScience，ChinaEarthquakeAdministration，Beijing100036，China4)FirstCrustMonitoringandApplicationCenter，ChinaEarthquakeAdministration，Tianjin300180，China

Using double-difference precise location results of seismic S wave data from 2000 to 2012，distribution and variation of dynamic moving averageQsvalue are imaged in Hainan Island and adjacent regions through seismic attenuation tomography. We explore regional crustal S-wave attenuation characteristics of temporal and spatial variation combined with seismic activity，volcanoes and heat distribution，and GPS baseline changes. The results show that: (1)Attenuation imaging shows that there is a significant lateral heterogeneity inQsdistribution in the study area and theQsvalues are high in the central region and low in the north and south regions of Hainan Island. Low-Qsareas are baseically converged to the north of Wangwu-Wenjiao Fault，to the south of Jianfeng-Diaoluo Fault，and to the east of Puqian-Qinglan Fault. Earthquakes are basically converged to the edge of the transitional regions of high and lowQsvalues. Heat flow sites and volcanoes zones are located in low-Qsarea in the region. (2)There is a strong correlation between dynamic moving averageQsvalue，seismic activity and Luzhou-Qiongzhong GPS baseline. From 2000 to 2008，the averageQsvalue of the study area is relatively high，the seismic activity is strong，and Luzhou-Qiongzhong GPS baseline is decreasing. From 2008 to 2012，the averageQsvalue of the study area shows a downward trend，the seismic activity is weak，and the Luzhou-Qiongzhong GPS baseline displays an elongation trend.

Qstomography，seismic activity，volcanoes，heat flow，GPS baseline，Hainan Island and adjacent regions

10.3969/j.issn.0253- 4967.2016.03.010

2015-08-06收稿，2015-11-21改回。

地震科技星火计划项目(XH12037Y)资助。

P315.3

A

0253-4967(2016)03-0631-15

王惠琳，女，1982年生，2011年毕业于中国地震局地震预测研究所固体地球物理专业，获硕士学位，助理研究员，主要从事地震层析成像研究、地震活动性工作，电话： 18889827682，E-mail: huilin95@163.com。