基于GPS观测网的山东地区地壳运动特征分析

第３０卷第３期　西北地震学报　Ｖｏ１．３０　Ｎｏ．３　Ｓｅｐｔ．，２００８　２（Ｘ）８年９月　ＮＯＲＴＨＷＥＳＴＥＲＮ　ＳＥＩＳＭ０ＬＯＧＩＣＡＬ　ＪＯＵＲＮＡＬ　基于ＧＰＳ观测网的山东地区地壳运动特征分析　殷海涛　，李　杰　，张　玲　，吴　晨　，董晓娜　（１．山东省地震局，山东济南２５００１４；２．中国地震局地质研究所，北京　１０００２９）　摘要：基于山东地壳运动ＧＰＳ观测网络（ＣＭＯＮＯＳＤ）和周边国际ＩＧＳ站点提供的２００７年全年的　观测数据，利用ＣＡＭＩＴ／ＧＬＯＢＫ软件进行了分析处理，得出各基准站间的基线值和站点速度场，并　对山东省内主要断裂带——沂沭断裂带两侧的地壳运动状况以及山东相对于周边地区的运动状况　进行了分析。结果表明山东地区不同块体的运动状态相对平稳，且随中国大陆东部较均匀地向东　南方向运动。　关键词：山东；地壳运动；ＧＰＳ；沂沭断裂带；基线　中图分类号：Ｐ３１５．７２　文献标识码：Ａ　文章编号：１０００—０８４４（２００８）０３—０２７６—０６　Ａｎａｌｙｓｉｓ　ｏｆ　Ｃｒｕｓｔａｌ　Ｍｏｖｅｍｅｎｔ　Ｆｅａｔｕｒｅｓ　ｉｎ　Ｓｈａｎｄｏｎｇ　Ａｒｅａ　Ｂａｓｅｄ　ｏｎ　ｔｈｅ　Ｄａｔａ　０ｆ　ＧＰＳ　Ｏｂｓｅｒｖａｔｉｏｎ　Ｎｅｔｗｏｒｋ　ＹＩＮ　Ｈａｉ—ｔａｏ　，ＬＩ　Ｊｉｅ　，ＺＨＡＮＧ　Ｌｉｎｇ　，ＷＵ　Ｃｈｅｎ　，ＤＯＮＧ　Ｘｉａｏ—Ｈａ　（１．Ｅａｒｔｈｑｕａｋｅ　Ａｄｍｉｎｉｓｔｒａｔｉｏｎ　ｏｆ　Ｓｈａｎｄｏｎｇ　Ｐｒｏｖｉｎｃｅ．Ｊｉｎａｎ　２５００１４。Ｃｈｉｎａ；　２．Ｉｎｓｔｉｔｕｔｅ　ｏｆ　Ｇｅｏｌｏｇｙ，Ｃｈｉｎａ　Ｅａｒｔｈｑｕａｋｅ　Ａｄｍｉｎｉｓｔｒａｔｉｏｎ，Ｂｅｉｊｉｎｇ　１０００２９，Ｃｈｉｎａ）　Ａｂｓｔｒａｃｔ：Ｂａｓｅｄ　ｏｎ　ｔｈｅ　ＧＰＳ　ｄａｔａ　ｆｒｏｍ　ｔｈｅ　ｃｒｕｓｔａｌ　ｍｏｖｅｍｅｎｔ　ＧＰＳ　ｏｂｓｅｒｖａｔｉｏｎ　ｎｅｔｗｏｒｋ　ｏｆ　Ｓｈａｎｄｏｎｇ（ＣＭ－　０ＮＯＳＤ）ａｎｄ　ＩＧＳ　ｓｔａｔｉｏｎｓ　ｉｎ　２００７，ｔｈｅ　ｍｏｖｅｍｅｎｔ　ｏｆ　ｂｏｔｈ　ｓｉｄｅ　ｏｆ　Ｙｉｓｈｕ　ｆａｕｌｔ　ｚｏｎｅ，ｗｈｉｃｈ　ｉＳ　ｔｈｅ　ｍａｉｎ　ｆａｕｌｔ　ｚｏｎｅ　ｉｎ　Ｓｈａｎｇｄｏｎｇ　ｐｒｏｖｉｎｃｅ，ｉｓ　ａｎａｌｙｚｅｄ．Ｔｈｅ　ｂａｓｅｌｉｎｅｓ　ｖａｌｕｅｓ　ａｍｏｎｇ　ｓｔａｔｉｏｎｓ　ａｎｄ　ｖｅｌｏｃｉｔｙ　ｆｉｅｌｄｓ　ｆｏｒ　ｓｔａｔｉｏｎｓ　ａｒｅ　ｃａｌｃｕｌａｔｅｄ　ｂｙ　ｕｓｉｎｇ　ＧＡＭＩＴ／ＧＬＯＢＫ　ｓｏｆｔｗａｒｅ．Ｔｈｅ　ｒｅｓｕｌｔｓ　ｓｈｏｗ　ｔｈａｔ　ｔｈｅ　ｍｏｖｅｍｅｎｔ　ａｍｏｎｇ　ｄｉｆ－　ｆｅｒｅｎｔ　ｂｌｏｃｋｓ　ｉｎ　Ｓｈａｎｄｏｎｇ　ａｎｄ　ａｄｊａｃｅｎｔ　ａｒｅａ　ｉｓ　ｒｅｌａｔｉｖｅｌｙ　ｓｔａｂｌｅ，ａｎｄ　ｗｈｏｌｅ　ａｒｅａ　ｉｓ　ｍｏｖｉｎｇ　ｔｏ　ｓｏｕｔｈｅａｓｔ　ｓｍｏｏｔｈｌｙ　ａｌｏｎｇ　ｗｉｔｈ　ｔｈｅ　ｅａｓｔｅｒｎ　ｐａｒｔ　ｏｆ　Ｃｈｉｎｅｓｅ　ｍａｉｎｌａｎｄ．　Ｋｅｙ　ｗｏｒｄｓ：Ｓｈａｎｄｏｎｇ　ｉ　Ｃｒｕｓｔａｌ　ｍｏｖｅｍｅｎｔ；ＧＰＳ；Ｙｉｓｈｕ　ｆａｕｌｔ　ｚｏｎｅ；Ｂａｓｅｌｉｎｅ　０引言　大陆岩石圈系统由多个相互作用的块体网络构　成，地震是此系统过程的一种行为　ｊ。大地震的演　化、发震及震后必然与大空间范围的系统过程有关，　义　］。　山东地壳运动ＧＰＳ观测网的建立弥补了山东　地区缺少大范围水平形变测量手段的不足　，为监　控山东地区主要构造块体、主要断裂的活动水平提　供了连续可靠的观测数据，同时为区域连续ＧＰＳ形　变监测系统的建设和数据分析处理积累经验，为区　域地震监测预报研究和防震减灾服务。　本文对山东地壳运动ＧＰＳ观测网络产出的　这已逐步成为中外地学界的共识　。连续ＧＰＳ台　网提供了一种全新的监测工具，实现了“地学望远　镜”（大空间）和“地学显微镜”（精细的时变过程）　的相结合，使我们有可能从系统过程来监测预测大　地震，或为中短期地震预测提供重要信息　。连续　ＧＰＳ网为研究随时间演进的应变积累和应力变化等　“微动态”创造了条件，近几年国际上已有～些初步　研究成果，这对中长期地震危险区判定有重要的意　收稿日期：２００８４３１－２３　２００７年全年的ＧＰＳ观测数据与周边地区的ＧＰＳ观　测数据进行了联合解算，分析了不同基线的变化特　征及其与构造的关系。　基金项目：“十五”山东省防震减灾重点项目（ＳＤ５０１—０４）和２００８年山东省奥运安保项目（２ｏｏ８ＧＧ１Ｏｏ０８Ｏ０４）　作者简介：殷海涛（１９８０一），男，山东聊城人，助理研究员，博士，主要从事ＧＰＳ高精度解算及地壳形变分析研究　维普资讯 http://www.cqvip.com

第３期　殷海涛等：基于ＧＰＳ观测网的山东地区地壳运动特征分析　２７７　１　山东ＧＰＳ观测网概况及数据处理　１．１　山东ＧＰＳ观测网概况　（ＲＩＺＨ）、苍山（ＣＡＳＨ）、无棣（ＷＵＤＩ）和泰安　（ＴＡＩＮ）基准站（该站由ＣＭＯＮＯＣ建造）（图１），均　配置徕卡ＧＲＸ１２００Ｐｒｏ接收机和ＡＴ５０４扼流圈天　线。ＧＰＳ数据中心设在山东省防震减灾中心。于　２００６年１Ｏ月全部建成，并于２００７年１月１日进入　山东地壳运动ＧＰＳ观测网络（ＣＭＯＮＯＳＤ）由６　个基准站和１个数据中心组成。基准站分别为分布　在山东内陆具有差异运动的构造块体内部、位于沂　沭断裂带两侧的烟台（ＹＡＴＩ）、即墨（ＪＩＭＯ）、日照　１１４。　１１　８。　试运行阶段，２００７年７月１日进入正式运行。　１２２。Ｅ　　｝９　／Ｌ　　／、＾无，麓　　番＼　　’‘　／／建．５　　。　Ｎ　３８。　Ｌ　９一～　●数据中心　▲已建基准站　０拟建基准站　图Ｉ　山东地壳运动ＧＰＳ观测网站点分布图　Ｆｉｇ．１　Ｄｉｓｔｒｉｂｕｔｉｏｎ　ｏｆ　ｓｔａｔｉｏｎｓ　ｏｆ　ｃｒｕｓｔａｌ　ｍｏｖｅｍｅｎｔ　ＧＰＳ　ｏｂｓｅｒｖａｔｉｏｎ　ｎｅｔｗｏｒｋ　ｉｎ　Ｓｈａｎｄｏｎｇ　ｐｒｏｖｉｎｃｅ．　１．２数据处理方法　本文所用的数据处理软件是由美国麻省工学院　和Ｓｃｒｉｐｐｓ海洋研究所研制的、目前国际上三大著名　高精度定位软件包之一ＧＡＭＩＴ／ＧＬ０ＢＫｌ０．３２版，并　辅以一系列自主编程的精度统计、形变分析工具。　◆极移和ＵＴ１（经极移改正的世界时）采用ｕＳ．　ＮＯ的Ｂｕｌｌ—Ａ值并给予强约束；　◆ＩＧＳ站东西向、南北向和垂向分别以２　ａｍ、２　ｃｍ和５　ｃｍ约束于ＩＧＳ０５，山东基准站位置的先验约　束为１　ｍ；　◆地球重力场模型、固体潮模型、极潮模型采用　ＩＥＲＳ２００３的最新规范；　数据解算方法如下：　ＧＡＭＩＴ软件采用双差模式完成相对定位，可以　很好地消除接收机钟和卫星钟的钟差影响，同时也　可以明显减弱诸如轨道误差、大气折射误差等系统　性误差的影响。　◆海潮引发的测站地壳形变采用ＦＥＳ２００４全　球海潮模型进行改正，并同时顾及海潮导致的地球　质心的变化；　ＧＰＳ信号从发射、传播到接收，受到各种地球物　理效应和误差源的影响。在数据处理中，误差的消　◆采用Ｌｃ线性组合消除电离层折射的延迟；　◆大气层干分量的天顶延迟由模型计算获得，　除或减弱通常采用模型改正或附加参数估计的方　法。这些模型、方法经历着一个不断完善、出新的过　程。我们基于目前国际上最新的研究成果¨　Ｈ］，确　定了下列数据处理的方法和模型，以期获得最佳的　结果。数据处理以２４小时为一个处理时段。　◆采用轨道松弛模式，为后期与ＩＧＳ全球解的　湿分量的天顶延迟由参数估计确定，每个站每２个　小时估计１个天顶延迟参数，映射函数采用ＧＭＦ　模型。顾及大气的不均匀性，对每个站的东西向和　南北向各附加１个大气水平梯度参数；　◆采用绝对的ＥＬＥＶ模型改正卫星和接收机天　线的相位中心变化；　◆数据采样率：３０秒数据编辑，１２０秒基线解　算；　结合提供结合点；　◆采用ＩＧＳ快速精密星历（ＩＧＳ精密星历延迟２　周），并对卫星轨道给予ｌ０　Ｉ９量级的约束；　◆卫星截止高度角取为ｌ０。，观测值采用误差　维普资讯 http://www.cqvip.com

２７８　西北地震学报　第３Ｏ卷　模型ｌ＇ｍｓ＝ａ　＋ｂ２／ｓｉｎ（ｅｌｅｖ）　（ｍｍ）定权，模型中的　ａ和ｂ值由观测值的残差分析拟合获得。　对ＧＡＭＩＴ得到的单天基线解算结果和ＳＯＰＡＣ　提供的全球单天解算结果（ｉｇｓｌ，ｉｇｓ２，ｉｇｓ３等子网　络），利用ＧＬＯＢＫ软件进行综合解算。ＧＬＯＢＫ按　照Ｋａｌｍａｎ滤波方法，综合全部伪观测值计算观测站　坐标和速度Ｈ孓　。　设对应历元ｔ　的线性化伪观测值向量伪观察　者向量为６１　＝　一Ｚ。　，其中　为初始值，ｆ　为　ＧＡＭＩＴ得到的伪观测值估值，则Ｋａｌｍａｎ滤波的观　测方程为　６１　＝Ａ　８ｘ　＋占　（１）　其中融　为待估参数；Ａ　为系数矩阵；　为伪观测　值２　对应的观测噪声，并且Ｅ｛　｝＝０，Ｅ｛￣＇ｋｏｏ　＝　，Ｅ｛ｇ　＋　｝＝０＾。Ｋａｌｍａｎ滤波的状态方程为　＋１＝Ｓ　８ｘ　＋ｑ　（２）　其中｜ｓ　为状态转移矩阵；ｑ　为状态噪声，其协方差　矩阵为Ｑ　，且Ｅ｛ｑ　ｓ　＝０。则Ｋａｌｍａｎ滤波从￡　到　＋　历元的预测方程和改进方程分别为ｓ氟　＝　．ｓ　互　，ｃ　＋ｌ　Ｉ　＝　｜ｓ　Ｃ　｜ｓ　。　＋　Ｑ　，　互　＋１　：　＋１１　＋　＋１（６ｆ　＋１　一　Ａ　＋１　＋１ｌ　），Ｃ　＋１　＝　Ｃ　＋１　１一　＋　Ｃ川…其中　Ｃ　为参数　，估值的协方差矩　阵。则Ｋａｌｍａｎ滤波增益矩阵为　Ｋ　＋１＝ｃ　＋１ｌ　Ａ　。＋ｌ（Ｐ　＋１＋Ａ　＋ｌｃ　＋１　１Ａ　＋１）Ｉ１　（３）　２结果及分析　２．１　山东内部各基准站相对基线变化分析　利用上述方法对山东地壳运动ＧＰＳ观测网　２００７年产出的观测数据进行了解算。　２．１．１基线变化分析　从约束基准的整周模糊度整数解，可以获取各　基准站间的基线三分量和基线长度。通过分析基线　时间序列的变化，可以了解基线所跨断裂及基准站　所在不同构造块体的相对运动状态¨　。　众所周知，ＧＰＳ定位的内符精度通常较高，以此　衡量ＧＰＳ成果的精度缺乏可靠性。通常用来衡量　连续观测精度的指标是均方连差ｍ：　（４）　其中，　为相邻两天的观测值之差；ｎ为　的个数。　山东地壳运动ＧＰＳ观测网可围绕沂沭断裂带　产出ｌ５条基线。沂沭断裂带是中国东部著名的郯　（城）庐（江）深断裂在山东境内的区段，也是山东最　主要的深断裂带，带上曾发生过郯城８．５级（１６６８—　０７－２５）大震，也是山东地区大地震的主要活动场　所¨　。根据ＧＰＳ观测成果研究它的现代运动状态，　对地震危险区的判定有着重要意义。基线解算结果　如表１所示。通过对２００７年各条基线的分析，我们　可以初步发现如下特征：沂沭断裂带西侧的ＴＡＩＮ—　ＣＡＳＨ、ＴＡＩＮ—ＷＵＤＩ和ＷＵＤＩ—ＣＡＳＨ呈缩短变　化，沂沭断裂带东侧的ＲＩＺＨ—ＹＡＴＩ、ＹＡＴＩ—ＪＩＭＯ　基线也均呈缩短变化。这说明山东南北断块主要受　压性应力作用。而沂沭断裂带西侧的ＷＵＤＩ、ＴＡＩＮ　和ＣＡＳＨ相对于东侧的ＪＩＭＯ、ＲＩＺＨ的基线变化略　呈伸长趋势变化，但量级不大，说明沂沭断裂带受张　性应力作用。整体来看，山东地区地壳运动较平稳，　沂沭断裂带目前表现为张扭性活动背景。　表１基线变化及精度　２．１．２站点速度场分析　在全球范围内选择较均匀分布的ＩＧＳ０５核心站　作为参考站，求解全球无基准解与ＩＧＳ０５间的相似　变化７参数（３个平移、３个旋转、一个尺度）；通过７　参数将全球无基准解转换到ＩＴＲＦ２００５参考框架下。　为了更好地表征山东地区的内部运动特征，将　ＩＴＲＦ２００５框架下的速率扣除ＮＮＲ—ＮＵＶＥＬ１Ａ地质　模型计算的刚性运动速率，得到山东地区相对于欧　亚板块的速度场（图２，表２）。可见，由于各站向南　运动速率随纬度的减小而减小，说明山东地区受由　维普资讯 http://www.cqvip.com 第３期　殷海涛等：基于ＧＰＳ　Ｙ￣　些垦垫塞壅　壁　坌堑　堋　川　北向南的作用力逐渐减弱。山东各ＧＰＳ基准站相　于鲁西块体主要受来自青藏块体北东向挤压，鲁东　对于欧亚板块的总体运动趋势均呈现向东偏南方向　一黄海块体主要受太平洋块体和菲律宾块体向西北　运动，且与中国大陆东部地区地壳运动方向非常一　推挤，两种作用力的合成效果促成了山东地区主应　致（图３）。从整体上看，山东地区的主应变场呈　变场的这一特性　。　ＮＥＥ—ＳＷＷ向挤压，ＮＷＷ—ＳＥＥ向拉张。这是由　¨　４　５　７　５　３　¨”６　ｍ　图２　山东ＧＰＳ基准站速度矢量图　Ｆｉｇ．２　Ｔｈｅ　ｖｅｌｏｃｉｔｙ　ｖｅｃｔｏｒｓ　ｆｏｒ　ＧＰＳ　ｓｔａｔｉｏｎｓ　ｉｎ　Ｓｈａｎｇｄｏｎｇ　表２山东地壳运动ＧＰＳ观测网的站点速度　序列和线性趋势线，它表明这两个站之间一直处于　（ＩＴＲＦ２００５框架下）　压缩状态，在２００５年１１月本州以东海域（Ｎ３８．２。，　站名　仙　ｖｎ　Ｅ１４４．７。）的７．１级地震以后，先是有所反弹然后就　。　恢复到以前的运动水平（约一２２　ｍｍ／ａ），但是２００６　年１１月千岛群岛（Ｎ４６．６。，Ｅ１５３．３。）８．０级地震发　生以后先是有所反弹然后就呈现加速挤压的状态　（约一３５　ｍｍ／ａ），显示太平洋板块对亚欧板块的俯　冲速度加快，其推挤作用的进一步增强，有利于中国　２．２　山东地区相对周边区域的运动情况分析　大陆东部地壳应力的积累　，由于日本南部与鲁　由上一节的结果可以看出，山东地区的主要断　东一黄海地块地震活动具有一定的相关性和同步　性　２引。从动力学观点来看，需要关注目前加速状态　裂带两侧地壳运动状态较平稳，但是要把握该地区　的整体运动状况，还需要结合周边地区ＧＰＳ站点数　及速率出现转折后对中国大陆东部及海域的地震活　据来分析它们之间的相对运动状况。所以，我们选　动情况。　择了周边的６个ＩＧＳ站点：ＢＪＦＳ、ＫＵＮＭ、ＷＵＨＮ、　３结论与探讨　ＵＲＵＭ、ＴＳＫＢ和ＩＲＫＴ，与山东的ＴＡＩＮ站之间的基　通过以上分析可以看出，山东主要断裂带一沂　线进行了分析，结果见表１，站点分布见图３。２００７　沭断裂带同侧站点间基线呈挤压状态，断裂带两侧　年，ＢＪＦＳ、ＷＵＨＮ和ＩＲＫＴ到ＴＡＩＮ的三条基线变化　站间基线略呈拉张状态，但量级均不大，表明沂沭断　量较小，＜３　ｍｍ；ＫＵＮＭ和ＵＲＵＭ至ＴＡＩＮ的基线长　裂带目前表现为张扭性活动背景，断裂带两侧鲁西　度有明显的拉张，６～１０　ｍｍ。从图中可以看出，山　和鲁东一黄海地块整体运动相对平稳。山东地壳运　东地区的形变场与中国大陆东部形变场的整体背景　动速度场显示块体整体向东南方向运动，与中国东　基本保持一致。　部相邻块体的运动方向基本一致，不存在明显的差　图４给出了日本的ＴＳＫＢ至ＴＡＩＮ基线的时间　异运动。但是由于ＴＳＫＢ—ＴＡＩＮ基线自２００６年１１　维普资讯 http://www.cqvip.com

２８０　西北地震学报　第３Ｏ卷　月千岛群岛８．０级地震后，缩短速度明显加快（提　国大陆东部地壳应力的积累，因此要需注意其变化　高了１．５倍），其推挤作用的进一步增强，有利于中　趋势以及中国大陆东部及海域的地震活动情况。　图３　中国大陆及周边地区运动速度图　Ｆｉｇ．３　Ｍｏｖｅｍｅｎｔ　ｖｅｌｏｃｉｔｙ　ｏｆ　ｍａｉｎｌａｎｄ　ｏｆ　Ｃｈｉｎａ　ａｎｄ　ａｄｊａｃｅｎｔ　ａｒｅａ．　蠢２．０５２６５８　蝴　图４　ＴＳＫＢ—ＴＡＩＮ（２００４－－２００７）基线变化图　Ｆｉｇ，４　Ｃｈａｎｇｅｓ　ｏｆ　ｂａｓｅｌｉｎｅ　ｆｒｏｍ　ＴＳＫＢ　ｔｏ　ＴＡＩＮ（２００４－－２００７）　由于该系统运行时间不长，且站点和数据资料　可靠的依据。　相对较少，所以对地壳运动的背景研究受到了一定　山东地壳运动ＧＰＳ观测网数据处理过程中得　的限制。随着国家“十一五”陆态网的建设（山东将　到了中国地震局王敏研究员的帮助和指导，在此深　增加３个基准站），数据资料的不断丰富，有待于进　表感谢！　一步深入的研究，为山东今后的防震减灾工作提供　维普资讯 http://www.cqvip.com 第３期　殷海涛等：基于ＧＰＳ观测网的山东地区地壳运动特征分析　∞　［参考文献］　ｆ１］郑剑东．浅谈大陆构造：板块构造的回顾与反思［Ｊ］．地球科　学进展，１９９１，６（３）：５０－５５．［２］　张培震，邓起东，张国民，等．中国大陆的强震活动与活动地块　２００５ＧＬ０２５５３８・　［１３］Ｌｙｍｘｔ　Ｆ，Ｌｅｆｅｖｒｅ　Ｆ，ＬｅｔｅｌｌｉｅｒＴ，ｅｔ　ａｌ・Ｍｏｄｅｌｉｎｇｔｈｅ　ｇｌｏｂａ１　ｏｅｅａｎ　ｔｉｄｅ８：ｉｎｓｉｇｈｔｓ　ｆｒｏｍ　ＦＥＳ２００４［Ｊ］・Ｏｃｅａｎ　Ｄｙｎａｍｉｃｓ，２００６，５：　３９４４１５・　［Ｊ］．中国科学（Ｄ辑），２００３，３３（Ｂ０４）：１２．２０．　［３ｊ　沈正康．关于地震预报研究发展战略的思考［Ｊ］．国际地震动　态，２００５．（５）：６９－７３．　［４］　王琪，张培震，马宗晋．中国大陆现今构造变形ＧＰＳ观测数据　与速度场［Ｊ］．地学前缘，２００２，９（２）：４１５４２９．　［１４］　Ｂｏｅｈｍ　Ｊ，Ｎｉｅｌｌ　Ａ，Ｓｃｈｕｈ　Ｈ・Ｇｌｏｂ￣Ｍａｐｐｉｎｇ　Ｆｕｎｃｔｉｏｎ（ＧＭＦ）：　Ａ　ｎｅｗ　ｅｍｐｉｉｆｅａｌ　ｍａｐｐｉｎｇ　ｒｕｎｅｔｉ。ｎ　ｂａｓｅｄ　ｏｎ　ｎｕｍｅｒｉｃａｌ　ｗｅａｔｈｅｒ　ｍｏｄｅｌ　ｄａｔａ［Ｊ］．Ｇｅｏｐｈｙｓｉｃａｌ　Ｒｅｓｅａｒｃｈ　Ｌｅｔｔｅｒｓ，２００６，３３．　Ｉ／）７３０４，ｄｏｉ：１０・１０２９／２００５ＧＬ０２５５４６・　［１５］Ｈｅｒｒｉｎｇ　Ｔ　Ａ・ＧＬＯＢＫ：Ｇｌｏｂａｌ　Ｋａｌｍａｎ　ｉｆｌｔｅｒ　ＶＬＢＩ　ａｎｄ　ＧＰＳ　ａｎａｌ一　［５］　江在森，张希．华北地区近期地壳水平运动与应力应变场特征　ｙｓｉｓ　ｐｒｏｇｒａｍ（Ｖｅｒ　１０・０）［Ｒ］・Ｍａｓｓａｃｈｕｓｅｔｔｓ，Ｕ・Ｓ・Ａ・：Ｍａｓｓａ。　［Ｊ］．地球物理学报．２０００，４３（５）：６５７－６６５．　［６］　王敏，沈正康，牛之俊，等．现今中国大陆地壳运动与活动块体　模型［】］．中国科学（Ｄ辑），２００３：３３（Ｂ０４）：２１－３２．　［７］　张永志，崔笃信，王琪，等．利用ＧＰＳ资料研究区域地壳应力　场变化与地震活动关系［Ｊ］．地震学报．２０００，２２（５）：４４９．　４５６．　［８］　李建彪，甘卫军，冉勇康，等．青藏高原东部构造块体的运动学　及形变特征分析［Ｊ］．西北地震学报，２００６，２８（２）：９７　１０３．　［９］　李杰，殷海涛，张玲，等．山东地壳运动ＧＰＳ观测网的建设与初　步结果分析［Ｊ］．大地测量与地球动力学，２００７，２７（增刊）：９—　１３．　［１Ｏ］　Ｇｅｎｄｔ　Ｇ．［ＩＧＳＭＡＩＬ５４３８］：ＩＧＳ　ｓｗｉｔｃｈ　ｔｏ　ａｂｓｏｌｕｔｅ　ａｎｔｅｎｎａ　ｍｏｄ．　ｅｌ　ａｎｄ　ＩＴＲ１７２００５．　［１】］　ＭｅＣｅｄ￣ｈｙ　Ｄ，Ｐｅｔｉｔ　Ｇ．ＩＥＲＳ　Ｃｏｎｖｅｎｔｉｏｎｓ（２００３），ＩＥＲＳ　Ｔｅｃｈｎｉｃａｌ　Ｎｏｔｅ　３２．２００４，Ｆｒａｎｋｆｕｒｔ，Ｇｅｒｍａｎｙ．　［１　２］Ｗａｔｓｏｎ　Ｃ，Ｔｒｅｇｏｎｉｎｇ　Ｐ，Ｃｏｌｅｍａｎ　Ｒ．Ｉｍｐａｃｔ　ｏｆ　ｓｏｌｉｄ　Ｅａｒｔｈ　ｔｉｄｅ　ｍｏｄｅｌｓ　ｏｎ　ＧＰＳ　ｃｏｏｒｄｉｎａｔｅ　ａｎｄ　ｔｒｏｐｏｓｐｈｅｉｒｃ　ｔｉｍｅ　ｓｅｉｒｅｓ［Ｊ］．Ｇｅｏ—　ｐｈｙｓｉｃａｌ　Ｒｅｓｅａｒｃｈ　Ｌｅｔｔｅｒｓ，２００６，３３．１／）８３０６，ｄｏｉ：１０．１０２９／　ｅｈｕｓｅｔｔｓ　Ｉｎｓｔｉｔｕｔｅ　ｏｆＴｅｃｈｎｏｌｏｇｙ，２００２－　［１６］　冯金涛。亚太Ｎ－－Ｎ　ＧＰＳ观测数据处理结果及地壳形变特　征分析［Ｊ］・西北地震学报，２００６，２８（３）：２２１－２２４・　［１７］　王双绪，丁平，张希，等．利用ＧＰＳ资料研究印尼８・７级地震　对我国西部地震趋势的影响［Ｊ］．西北地震学报，２００５，２７　（４）：３２２—３２７．　［１８］　刘兴旺，袁道阳，葛伟鹏．兰州黄河阶地高精度ＧＰＳ测量与　构造变形研究［Ｊ］．西北地震学报，２００７，２９（４）：３４１－３５１．　［１９］　王志才，贾荣光，孙昭民，等．沂沭断裂带安丘一莒县断裂安　丘一朱里段几何结构与活动特征［Ｊ］．地震地质，２００５，２７　（２）：２１２—２２０．　［２Ｏ］　张跃刚，帅平，胡新康，等．从ＧＰＳ观测看华北地区的形变场　演化［Ｊ］．大地测量与地球动力学，２００６，２６（１）：３６４１．　［２１］　章纯．中国东部地区地震活动与构造应力场关系的有限元数　值模拟［Ｊ］．西北地震学报，２００７，２９（３）：２３０－２３４．　［２２］　郑建常，蒋海昆．鲁东一黄海地块与日本南部地震活动相关　性分析及因果关系检验［Ｊ］．地震学报，２００７，２９（４）：３５８．　３６８．