日本M9．0级巨震对山东地区地壳活动的影响研究

殷海涛，甘卫军，黄 蓓，肖根如，李 杰，朱成林

1 山东省地震局，济南 250014

2 中国地震局地质研究所 地震动力学国家重点实验室，北京 100029

1 引 言

北京时间2011年3月11日13时46分在日本本州东海岸附近海域发生9．0级巨大地震，并引发了强烈的海啸和严重的核泄漏等次生灾害．中国地震台网测定，该地震震中位于北纬38．1°，东经142．6°，震源深度约20km（http：//www．csi．ac．cn/）．美国地质调查局（USGS）公布的地震构造背景资料显示，此次9．0级地震是由太平洋板块和北美板块边缘的俯冲带附近的逆冲断层运动结果（http：//earthquake．usgs．gov/），该地震是1990 年以来世界上第四大规模的地震．地球物理学家肯·赫德纳特，依据NASA 收集的资料，发现这次强震使日本本州岛向东移动大约2．4m，使地球自转加快1．6μs．所以，此次9．0级巨震必定会给太平洋板块及其边界或邻近板块的活动带来一定的影响．

根据张培震等（2003）活动地块划分结果［1］，山东及其东部沿海区域被划分在华北地块区的次级华北平原地块和鲁东—黄海地块中．日本与山东及其沿海区域同属太平洋板块、欧亚板块及菲律宾板块的交汇区域，地球动力学背景上具有同源性．太平洋板块对欧亚板块的俯冲作用是华北地区地震活动的重要动力来源之一［2］．已有研究结果表明，位于太平洋板块俯冲带的日本地区的强震活动对我国华北地区的地震活动具有重要的影响［3－9］，因此研究相对活跃的日本地震对山东及其东部沿海区域地震活动性的影响是十分必要的．

为了发现此次日本地震对山东所处的鲁东—黄海地块和华北平原地块产生的影响，地震波对本地区块体的实时影响以及该区域在未来一段时间运动趋势将发生何种变化等问题，本文基于山东GPS地壳运动观测网的最新观测成果，对这一系列问题进行分析和探讨，为该区域的地震趋势判定及抗震设防提供了可靠依据．

2 数据收集及处理

沂沭断裂带位于郯（城）庐（江）深断裂的中部，是山东省东部最大最深的断裂带．它控制着山东半岛东部至渤海、黄海的地震活动．历史上曾发生过公元前70年安丘7级和1668年郯城8．5级强烈地震，也是山东地区大地震的主要活动场所［10］．为了监测其地震活动性、研究地震的孕育过程，国家重大科学工程“中国大陆构造环境监测网络”（以下简称陆态网络）项目组和山东省地震局在本地区具有差异运动的活动块体上布设了15个连续观测基准站，组成山东省地壳运动GPS 观测网［11］，站点分布见图1．选用采样间隔为30s和1s的两种数据，分别用于计算日均值形变和地表动态形变过程．数据处理由麻省理工大学研发的GAMIT/GLOBK10．4 版软件，并辅以一系列自主编程的精度统计和形变分析工具来完成．

图1 山东地壳运动GPS观测网站点分布图图中大写英文缩写为站点名称，SDRC－山东荣成，YATI－烟台，SDYT山东烟台，SDQD 山东青岛，SDCY 山东昌邑，RIZH 日照，WUDI无棣，SDZB山东淄博，SDLY山东临沂，CASH 苍山，TAIN 泰安，SDJX 山东嘉祥，LICH 聊城，HEZE菏泽．Fig．1 Distribution of stations of Shandong Crustal Movement GPS Observation Network

2．1 30s采样数据处理

GPS观测信号，从发射、传播到接收，受到了各种地球物理效应和误差源的影响，在GAMIT/GLOBK 软件中，采用双差模式完成相对定位，可以很好地消除接收机钟和卫星钟的钟差影响，同时也可以明显减弱诸如轨道误差、大气折射误差等系统性误差的影响．在数据处理过程中，对误差进行消除或减弱，通常采用模型改正或附加参数估计的方法．这些模型、方法都在经历着一个不断完善、创新的过程．基于目前国际上最新的研究成果［12－16］，我们确定了下列数据处理的模型和方法，以期获得最佳的解算结果．数据处理以单天（24h）为一个处理时段．

（1）采用轨道松弛模式，为后期与IGS 全球解的联合解算提供结合点；

（2）使用IGS最终精密星历（延迟约12～18天发布）或快速精密星历（延迟17～41小时发布），并对卫星轨道给予10－9量级的约束；

（3）UT1（经极移改正的世界时）和极移采用USNO 的Bull＿A 值并给予强约束；

（4）IGS站点的EW 向、NS向和U 方向分别以2cm、2cm 和5cm 约束于IGS05，2011年4月7日（DOY＝107，DOY＝Day of Year）后约束于IGS08框架；山东各基准站位置的先验约束设为1m；

（5）地球的重力场模型、固体潮模型和极潮模型均采用IERS2003的最新规范；

（6）利用FES2004全球海潮模型来对海潮引起的测站地壳变形进行修正，同时考虑了海潮造成的地球质心的变化；

（7）电离层折射延迟的消除主要采用LC 线性组合来完成；

（8）大气层干分量的天顶延迟由模型计算获得，湿分量的天顶延迟由参数估计确定，每个站每天估计13个天顶延迟参数（即2h估计一次），映射函数使用的是GMF模型．顾及大气的不均匀性，对每个站的东西向和南北向各附加1 个大气水平梯度参数；

（9）采用绝对的ELEV 模型改正卫星和接收机天线的相位中心变化；

（10）卫星高度截止角设为10°，观测值采用误差模型rms＝a2＋b2/sin（elev）2（mm）定权，模型中的a和b值均由观测值的残差分析拟合获得，elev为高度角；

GAMIT 解算完成后，再利用GLOBK 软件对得到的单天基线结果和SOPAC提供的全球网的单天解算结果（igs1，igs2，…，igs6）进行综合解算．GLOBK 按照Kalman滤波方法，综合全部伪观测值计算得到观测站的坐标和速度值［17］．

2．2 1s采样GPS数据处理

本文使用高频GPS数据处理软件是GAMIT/GLOBK 软 件 中 的TRACK 运 动 学 分 析 模 块［18－19］，该分析模块的核心问题就是单历元整周模糊度的解算．TRACK 模块首先利用伪距观测量得到初始值，然后再采用相对秩方法，利用卡方增量法对载波相位L1和L2模糊度整数值的最优结果与次优结果进行比较，如果最优值的卡方影响小于次优值，则其模糊度可确定为整数，然后使用迭代法直到解算出更多历元的模糊度．最后，利用Kalman滤波对对流层延迟参数和未解决的模糊度参数进行修正，得到最终结果．

在高频GPS数据处理过程中，对于尺度不同的观测网，其处理方法也有所不同：当基线＜10km时，数据解算较为简单，可通过双差法对电离层延迟进行有效消除，L1和L2的模糊度也能用浮点法分别计算出来；而像山东地壳运动GPS观测网这样的中大尺度观测网络，其站点间隔大于100km，解算过程相对复杂，TRACK 模块使用MW－WL 组合模型来消除电离层干扰，从而得到最优的模糊度解．

MW－WL组合，即相位宽巷组合与（MW）伪距窄巷组合（WL）之差，是由Wubbena和Melbourne在1985年［20－21］分别提出来的，其公式如下：

φi为双频相位观测值，Ri为双频伪距观测值，fi为频率（i＝1，2），c为光速．

MW－WL 组合由双频相位观测值和伪码组成，对于GPS观测量，Δf/Σf≈0．124，这表示其距离噪声减小了约1个数量级．从式（1）可看出，MW 组合削弱了几何距离、电离层、对流层及钟差等的影响，但仍受伪距多路径M以及观测噪声e的影响．

选择MW－WL组合的原因主要有两个：一是可免受电离层残差的影响，即不受基线长度的约束．二是可有效提高计算效率，对于实时计算系统，其计算效率也是非常重要的．MW 组合的另一个优点是不必考虑用户的运动状态而进行解算．虽然MW－WL组合可较好地解算整周模糊度，但是对于＞100km的长基线来说，观测数据自身的质量问题将会更直接地影响到结果的精度．

经过上述GPS解算，其结果中还包含有未模型化的系统误差和多路径效应等重复性误差，为此，我们使用了改进的恒星日滤波［22］和空间滤波［23］来进行处理，使结果更加可靠．

3 高频GPS获取的同震动态形变过程

近年来，随着GPS接收机技术和存储能力的提高，以及数据处理方法的不断改进，高频GPS技术得到了飞速的发展，大大提升了GPS的观测能力，可直接监测到地壳或建筑物的瞬时变化状态，国内外学者做了大量的研究工作［24－31］，并将此项技术逐步应用于震害防御、滑坡监测、火山监测以及气象学等领域．

本文利用1s采样的高频GPS数据获取了此次日本地震的地震波对山东所在地块的影响方式，通过研究其动态形变过程，可为本地区的抗震设防提供可靠依据，其解算结果见图2和图3．

图3 SDLY（a）和YATI（b）站点震时水平运动轨迹Fig．3 Horizontal movement track of SDLY（a）and YATI（b）station during the earthquake

从图2可以看出，山东各GPS基准站震动状态基本相似，地震波动时间和振幅随震中距的增大而逐渐减小．根据USGS 发布的结果：发震时间为：13∶46∶23．GPS记录到的地震波触发烟台（YATI）GPS 基准站地面发生亚厘米级波动的时间为13∶51∶26，距地震发生时刻约为306s．

以YATI站为例，其EW 向地壳运动状态为：首先向东运动，13∶51∶59到达波峰，随后便开始起伏波动，最大振幅约为0．224m，最大振幅持续时间约为160s．NS向地壳运动状态为：首先向南加速运动，触发时间约为13∶54∶32，距地震发生时刻约489s．随后加速向北运动，并接连出现两次北向波峰（随震中距的增大逐渐衰减为单波峰），最大波幅约为0．06m，持续时间约为130s．U 方向运动方式为：首先出现小幅的起伏波动，然后随着主震波的到来，地面呈现抬升并持续约120s，（持续时间各站点不同，随站点震中距的增大而减小），最大振幅约为0．24m，经振荡回落，基本恢复到震前水平．地震波造成的YATI站EW 向地面波动的持续时间约为700s，NS和U 方向约为510s．

在图3中，SDLY 站的NS和EW 向运动幅度相当，其轨迹近似为‘’形状，而YATI站EW 向运动幅度明显大于NS向的，故呈横向的‘λ’形运动．图3展现的是GPS站点受到地震波影响后的真实水平运动轨迹，反应了地表的真实运动方式．

4 同震位移场特征

采用30s采样的GPS数据来计算日均值，将地震发生当天（3月11日，DOY＝70）的结果与震前3天平均值结果进行比较，得到同震位移．GPS 观测到的不同站点的水平同震位移较好地揭示了日本地震对山东地区的影响程度（图4和表1），由于垂向形变较小，且与垂向误差（GPS 测定的垂向形变分辨率为厘米级）相当，故不再显示．

表1 山东地区GPS站点同震位移Table 1 Co－seismic displacement of GPS stations in Shandong

由图4可见，此次日本9．0级巨震对我国的东北和华北地区均产生了较大的影响，其水平同震位移均为朝震中方向运动，这也符合逆冲型地震的运动特征．山东各GPS基准站也发生了不同程度的同震形变，各站点位移量随震中距增加而递减，最大位移为超过10mm，NS和垂直方向变化相对较小，在误差范围内波动．在错动方向上，除处于胶东半岛东端的YATI，SDYT 和SDCY 站点外，其他各站同震位移方向基本一致，为东南方向．而此3个站点的位移方向为近正东方向，主要原因是这3个站点向东位移量明显大于其他各站，由于根据以往对沂沭带两侧站间基线的分析来看，沂沭带N 段呈现挤压状态［32］，故可初步判断此次日本地震对沂沭断裂带的同震效应使N 段产生了一定的拉张作用．

5 时间序列分析

为了更好地分析日本地震对山东所处地块的影响以及以后的发展趋势，首先来分析站点三维坐标时间序列，以处于华东地块的CASH 站和处于鲁东—黄海地块的RIZH 站为例（见图5），受日本地震的影响，使原本平稳向东运动的地块产生了一个加速东移，根据弹性回跳理论，并结合地震前后的实际运动速率，可以计算出震后3～5个月左右将恢复到震前的运动轨迹．

图5 站点坐标时间序列（m 为均方连差）Fig．5 Time series of station coordinate

为了得到山东主要断裂带——沂沭带受日本地震的影响，了解断裂带的活动状态，我们分析了跨断层基线的时间序列变化，结果见图6．

由图6a长期变化趋势可看出，跨沂沭断裂带的两条基线均具有非常明显的年周期性变化，只是在每年的雨季，受气压变化和相对湿度剧烈变化影响而产生了一定的扰动［33］．从长期变化趋势来看，处于N 段的WUDI－YATI基线一直以来均处于平稳的挤压状态，速率约为1．9mm/a，而S段的CASHRIZH 基线年变化非常平稳，没有明显的挤压和拉张．

从图6b短期变化趋势来看，基线CASH－RIZH此次受地震影响出现了一定的变化，但幅度很小，在误差范围内．而基线WUDI－YATI受此次日本地震的影响产生了较明显拉张，其主要变化发生在EW方向上，高程方向在地震前后并未发现明显的差异变化，也说明了日本地震对山东所在地块的影响主要是水平方向的，对垂直向形变影响不大．

6 结论与讨论

山东地壳运动GPS观测网得到的观测资料揭示了3月11日的日本9．0级巨震造成的山东所在地块的同震地表运动过程和同震位移场，以及震后运动趋势等信息，为本地区的抗震设防工作及震情趋势判定提供了可靠的依据，其初步结论如下：

（1）震时山东所处块体运动状态：EW 向最大振幅约为0．224m，持续时间约为160s．NS向地壳运动0．06 m，持续时间约为130s．地震波造成的NS地面波动的持续时间约为510s，EW 向约为700s．垂直方向地面抬升最大幅度约为0．24m，并持续约120s，持续时间各站点不同，随站点震中距的增大而减小．

（2）地震使山东所处地块产生了东偏南方向的错动．幅度最大的超过10 mm，错动量随震中距增加而递减，NS向和垂直方向变化相对较小，在误差范围内波动．目前各基准站均处于向西运动回调状态，经计算预计震后3～5个月左右将恢复到以前的运动趋势．由于YATI，SDYT 和SDRC基准站受日本地震的影响而产生的错动量明显高于其他站点，打破了原有较稳定的相对运动状态，故需密切关注胶东半岛东部及附近海域地区的震情形势．

（3）通过分析跨断裂带的基线来看，此次日本地震对沂沭断裂带起到了拉张的作用，其北段幅度明显大于南段的．

致 谢 感谢麻省理工大学Thomas Herring教授和Bob King教授在数据处理方面给予的支持与帮助，感谢中国地震局地质研究所张培震研究员给出的良好建议和帮助，感谢中国地震局地壳工程中心给予的大力支持，感谢两位审稿人对于本文提出的宝贵意见和建议！

（References）

［1］ 张培震，邓起东，张国民等．中国大陆的强震活动与活动地块．中国科学（D辑），2003，33（增刊）：12－20．

Zhang P Z，Deng Q D，Zhang G M，et al．Strong earthquake activity and active blocks in Chinese mainland．Sci．inChina（SeriesD）（in Chinese），2003，33（S1）：12－20．

［2］ 茂木清夫．巨大地震论．数理科学，1976，158：16－23．Mogi K．Huge earthquake theory．MathematicalSciences（in Chinese），1976，158：16－23．

［3］ 张国民．关于我国大陆地区地震大形势的研究．北京：地震出版社，1989．

Zhang G M．Study on Seismicity Character of Chinese Mainland（in Chinese）．Beijing：The Earthquake Press，1989．

［4］ 李延兴，张静华，何建坤等．由空间大地测量得到的太平洋板块现今构造运动与板内形变应变场．地球物理学报，2007，50（2）：437－447．

Li Y X，Zhang J H，He J K，et al．Current－day tectonic motion and intraplate deformation－strain field obtained from space geodesy in the Pacific Plate．ChineseJ．Geophys．（in Chinese），2007，50（2）：437－447．

［5］ 李延兴，张静华，李智等．太平洋板块俯冲对中国大陆的影响．测绘学报，2006，35（2）：99－105．

Li Y X，Zhang J H，Li Z，et al．The underthrust of Pacific plate to Eurasian plate and its effect on Chinese mainland．ActaGeodaeticaetCartographicaSinica（in Chinese），2006，35（20）：99－105．

［6］ 金双根 朱文耀．太平洋板块运动和形变及其边缘现今相对运动．大地测量学与测量工程，2002，22（2）：57－60．

Jin S G，Zhu W Y．The motion and deformation of pacific plate and contemporary relative motion at its boundaries．JournalofGeodesyandGeodynamics（in Chinese），2002，22（2）：57－60．

［7］ 马宏生，张国民，刘杰等．中国大陆及其邻区强震活动与活动地块关系研究．地学前缘，2003，10（Z1）：74－80．

Ma H S，Zhang G M，Liu J，et al．Correlation between strong earthquake activity and active crustal－block in China mainland and its adjacent regions．EarthSci．Frontiers（in Chinese），2003，8（10）：74－80．

［8］ 郑建常，蒋海坤．鲁东—黄海地块与日本南部地震活动相关性分析及因果关系检验．地震学报，2007，29（4）：358－368．

Zheng J C，Jiang H K．Correlation analysis and causality test between Ludong－Huanghai Block and South Japan．Acta SeismologicaSinica（in Chinese），2007，29（4）：358－368．

［9］ 张立敏，唐晓明．西太平洋板块俯冲运动与中国东北深震带．地球物理学报，1983，26（4）：331－340．

Zhang L M，Tang X M．The underthrusting movement of the western Pacific plate and the deep focus earthquake zone of Northeast China．ActaGeophysicaSinica（ChineseJ．Geophys．）（in Chinese），1983，26（4）：331－340．

［10］ 王志才，贾荣光，孙昭民等．沂沭断裂带安丘—莒县断裂安丘—朱里段几何结构与活动特征．地震地质，2005，27（2）：212－220．

Wang Z C，Jia R G，Sun Z M，et al．Geometry and activity of the Anqiu－Zhuli segment of the Anqiu－Juxian Fault in the Yishu Fault zone．SeismologyandGeology（in Chinese），2005，27（2）：212－220．

［11］ 殷海涛，李杰，张玲等．基于GPS观测网的山东地区地壳运动特征分析．西北地震学报，2008，30（3）：276－281．

Yin H T，Li J，Zhang L，et al．Analysis of crustal movement features in Shandong area based on the data of GPS observation network．NorthwesternSeismologicalJournal（in Chinese），2008，30（3）：276－281．

［12］ Gendt G．［IGSMAIL5438］：IGS switch to absolute antenna model and ITRF2005．

［13］ McCarthy D，Petit G．IERS Conventions（2003）．Frankfurt，Germany：IERS Technical Note 32．2004，．

［14］ Watson C，Tregoning P，Coleman R．Impact of solid Earth tide models on GPS coordinate and tropospheric time series．GeophysicalResearchLetters，2006，33（8），L08306，doi：10．1029/2005GL025538．

［15］ Lyard F，Lefevre F，Letellier T，et al．Modelling the global ocean tides： modern insights from FES2004．Ocean Dynamics，2006，56（5－6）：394－415．

［16］ Boehm J，Niell A，Tregoning P，et al．Global Mapping Function（GMF）：A new empirical mapping function based on numerical weather model data．GeophysicalResearchLetters，2006，33（7）：L07304，doi：10．1029/2005GL025546．

［17］ Herring T A．GLOBK：Global Kalman filter VLBI and GPS analysis program （Ver．10．0）．Massachusetts Institute of Technology，2002．

［18］ http：//geoweb．mit．edu/～tah/track＿example．

［19］ 殷海涛，甘卫军，肖根如等．利用高频GPS技术进行强震地面运动监测的研究进展．地球物理学进展，2009，24（6）：2012－2019．

Yin H T，Gan W J，Xiao G R，et al．Progress on monitoring strong earthquake ground motions using high－rate GPS．ProgressinGeophysics（in Chinese），2009，24（6）：2012－2019．

［20］ Wubbena G．Software developments for geodetic positioning with GPS using TI－4100code and carrier measurements．//Goad C C．Proceedings of the First International Symposium on Precise Positioning with the Global Positioning Sysetm，Maryland，1985：15－19．

［21］ Melbourne G．The case for ranging in GPS based geodetic systems．//Clyde G．Proceedings of the First International Symposium on Precise Positioning with the Global Positioning Sysetm，Maryland，1985：373－386．

［22］ 殷海涛，甘卫军，肖根如．恒星日滤波的修正以及对高频GPS定位的影响研究．武汉大学学报（信息科学版），2011，36（5）：609－611．

Yin H T，Gan W J，Xiao G R．Modified sidereal filter and it′s effect on high－rate GPS positioning．GeomaticsandInformation ScienceofWuhanUniversity（in Chinese），2011，36（5）：609－611．

［23］ 殷海涛，甘卫军，熊永良等．PCA 空间滤波在高频GPS定位中的应用研究．武汉大学学报（信息科学版），2011，36（7）：825－829．

Yin H T，Gan W J，Xiong Y L，et al．Study on the effect of PCA Spatial Filtering on high－rate GPS positioning．Geomatics andInformationScienceofWuhanUniversity（in Chinese），2011，36（7）：825－829．

［24］ Vigny C，Socquet A，Peyrat S，et al．The 2010Mw8．8 Maule Megathrust Earthquake of Central Chile，Monitored by GPS．Science，2011，332（6036）：1417－1421，DOI：10．1126/science．1204132．

［25］ Larson K M，Bodin P，Gomberg J．Using 1－Hz GPS data to measure deformations caused by the Denali Fault Earthquake．Science，2003，300（5624）：1421－1424．

［26］ Bock Y，Prawirodirdjo L， Melboume T．Detection of arbitrarily large dynamic ground motions with a dense highrate GPS network．GeophyRes．Lett．，2004，31（6），L06604，doi：10．1029/2003GL01950．

［27］ Miyazaki S，Laeson K，Choi K，et al．Modeling the rupture process of the 2003Septembter 25Tokachi－Oki（Hokkaido）earthquake using 1 Hz GPS data．Geophys．Res．Lett．，2004，31（21），L21603，DOI：10．1029/2004GL021457．

［28］ Ji C，Larson M，Tan Y，et al．Slip history of the 2003San Simon earthquake constrained by combining 1－Hz GPS，strong motion，and teleseismic data．Geophys．Res．Lett．，2004，31（17），L17068，doi：10．1029/2004GL020448．

［29］ Wang G，Boore D M，Tang G，et al．Comparisons of Ground motions from collocated and closely spaced one－sample－persecond Global Position System and accelerograph recordings of the 2003M6．5San Simeon，California earthquake in the Parkfield region．BullSeismol．Soc．Am．，2007，97（1B）：76－90．

［30］ Ohta Y，Meiano I，Sagiya T，et al．Large surface wave of the 2004Sumatra－Andaman earthquake captured by the very long baseline kinematic analysis of 1－Hz GPS data．Earth PlanetsSpace，2006，58（2）：153－157．

［31］ Langbein J，Bock Y．High－rate real－time GPS network at Parkfield： Utility for detecting fault slip and seismic displacements．Geophys．Res．Lett．，2004，31（12），L15S20，doi：10．1029/2003GL019408．

［32］ 殷海涛，李杰，董晓娜等．汶川地震前后山东及周边地区地壳运动特征分析．大地测量与地球动力学，2009，29（3）：23－27．

Yin H T，Li J，Dong X N，et al．On characteristics of crustal movement of Shandong and Surrounding areas before and after WenchuanMs8．0Earthquake．JournalofGeodesyand Geodynamics（in Chinese），2009，29（3）：23－27．

［33］ 王敏，沈正康，董大南．非构造形变对GPS连续站位置时间序列的影响和修正．地球物理学报，2005，48（5）：1045－1052．

Wang M，Shen Z K，Dong D N．Effects of non－tectonic crustal deformation on continuous GPS position time series and correction to them．ChineseJ．Geophys．（in Chinese），2005，48（5）：1045－1052．