地壳形变-地震活动的大地测量研究方法-2014

1、武汉大学武汉大学 地壳形变课程组地壳形变课程组地壳形变地壳形变地壳形变地壳形变 1 1、绪论、绪论2 2、地壳形变测量地壳形变测量 3 3、地球参考系与参考框架、地球参考系与参考框架4 4、板块构造学说与活动地块学说板块构造学说与活动地块学说 5 5、地壳运动监测与数据处理地壳运动监测与数据处理 6 6、地壳应力与应变分析、地壳应力与应变分析 7 7、连续形变、应变观测与数据处理、连续形变、应变观测与数据处理 8 8、地震活动的大地测量研究方法地震活动的大地测量研究方法1,1,中国地震活动的基本特点中国地震活动的基本特点 2,2,地震预报的科学思路地震预报的科学思路3,3,地震地壳形变地震地壳

2、形变 4,GPS4,GPS、InSARInSAR在地震研究中应用在地震研究中应用 5,5,地壳形变监测台网与地壳形变监测台网与地震预报地震预报 地震活动的大地测量研究方法地震活动的大地测量研究方法一、一、 中国地震活动的基本特点中国地震活动的基本特点 我国的地震活动，可用我国的地震活动，可用“多、大、广、浅多、大、广、浅”四个字概括其特点，即地震多四个字概括其特点，即地震多( (频度高频度高) )、强度大、分布广、震源浅、强度大、分布广、震源浅 。表表 10.110.1 中国分区强震频度统计（中国分区强震频度统计（19001900- -20032003） 震级震级 地区地区 7.07.0- -

3、7.47.4 7.57.5- -7.97.9 8.08.0- -8.48.4 8.58.5- -8.98.9 总和总和 大陆东部大陆东部 5 5 1 1 0 0 0 0 6 6 大陆西部大陆西部 2323 1515 5 5 2 2 4545 台湾省台湾省 2323 5 5 2 2 0 0 3030 其它地区其它地区 1 1 1 1 0 0 0 0 2 2 1900-2003 1900-2003年中国大陆强震震中分布图年中国大陆强震震中分布图 地震预测是一个国际性的科学难题，其难点有三：地震预测是一个国际性的科学难题，其难点有三： 一是地震物理过程本身的复杂性；一是地震物理过程本身的复杂性； 二

4、是地球内部的二是地球内部的“不可人性不可人性”即人们无法深入到地球内部在震源区内设即人们无法深入到地球内部在震源区内设置台站、安装观测仪器；置台站、安装观测仪器； 三是大震重复时间的长期性，限制了作为一门观测科学对现象的观测和三是大震重复时间的长期性，限制了作为一门观测科学对现象的观测和对经验认知的进展对经验认知的进展( (陈运泰，陈运泰，1998)1998)。从科学技术上来看要完全实现地震预测预报成功，需要达到两点：从科学技术上来看要完全实现地震预测预报成功，需要达到两点： 第一，对地震机理、地震孕育发生的物理过程的全面的完整的认识；第一，对地震机理、地震孕育发生的物理过程的全面的完整的认识

5、； 第二，要能够有效地监视地震的孕育、发生过程，以便根据孕震过程的第二，要能够有效地监视地震的孕育、发生过程，以便根据孕震过程的发展阶段进行地震预测。发展阶段进行地震预测。 目前的地震预测预报研究还处于探索阶段目前的地震预测预报研究还处于探索阶段 二、地震预报的科学思路二、地震预报的科学思路对一个自然现象的预测，往往有两种科学途径（张国民等，对一个自然现象的预测，往往有两种科学途径（张国民等，20012001）：）： 一一: : 研究并掌握该自然现象的生成机制和受控因素，通过测定有关因子的数研究并掌握该自然现象的生成机制和受控因素，通过测定有关因子的数值，按照该自然现象的成因规律对其作准确的预

6、测和预报。值，按照该自然现象的成因规律对其作准确的预测和预报。 二二: : 根据该自然现象与其他现象之间的关系，应用实践中积累的大量资料，根据该自然现象与其他现象之间的关系，应用实践中积累的大量资料，总结各种现象与预测对象之间的经验性和统计性关系进行预测和预报总结各种现象与预测对象之间的经验性和统计性关系进行预测和预报 。n 地震预报途径之一地震预报途径之一：是关于孕震过程和地震模式的理论和实验研究。孕震过：是关于孕震过程和地震模式的理论和实验研究。孕震过程的研究包括震源物理、地震力学等方面的理论、实验和观测研究，试图通程的研究包括震源物理、地震力学等方面的理论、实验和观测研究，试图通过对震源

7、过程物理力学机制的研究，逐步揭示和掌握地震孕育、发展和发生过对震源过程物理力学机制的研究，逐步揭示和掌握地震孕育、发展和发生的规律，从而达到预报地震的目的。地震模式的研究从一定的理论前提出发，的规律，从而达到预报地震的目的。地震模式的研究从一定的理论前提出发，提出地震发生的模式，从理论上推导各种可能的前兆及不同的关联组合，并提出地震发生的模式，从理论上推导各种可能的前兆及不同的关联组合，并通过实际观测不断检验和修改理论模式。通过实际观测不断检验和修改理论模式。 n 地震预报途径之二地震预报途径之二：是根据在长期实践中积累的大量震例资料，总结出经验：是根据在长期实践中积累的大量震例资料，总结出经

8、验性规律推广应用于预测未来地震。性规律推广应用于预测未来地震。两种地震预报途径两种地震预报途径地震预报方法的分类地震预报方法的分类 地震预报的地震学方法地震预报的地震学方法 地壳形变与地震预报地壳形变与地震预报 地下水微动态与地震预报地下水微动态与地震预报 水文地球化学地震前兆水文地球化学地震前兆 地震的地电前兆地震的地电前兆 地震的地磁前兆地震的地磁前兆 重力预报地震重力预报地震中国的地震预报在上述两条科学途径上探索前进，形成了依据地震异常群中国的地震预报在上述两条科学途径上探索前进，形成了依据地震异常群体特征对孕震过程实行追踪预报的科学思路，即通过大范围、长时间、多体特征对孕震过程实行追踪

9、预报的科学思路，即通过大范围、长时间、多手段前兆的连续观测，监视区域应力场的动态变化，探测其在正常背景上手段前兆的连续观测，监视区域应力场的动态变化，探测其在正常背景上的异常变化，并从场、源和环境相统一的整体观出发，分析异常群体的时的异常变化，并从场、源和环境相统一的整体观出发，分析异常群体的时空强综合特征及其演化过程：应用从大量震例经验和理论、实验研究取得空强综合特征及其演化过程：应用从大量震例经验和理论、实验研究取得的对孕震过程阶段性发展的认识，以及各阶段中异常群体特征的综合判据的对孕震过程阶段性发展的认识，以及各阶段中异常群体特征的综合判据与指标，对孕震过程进行追踪分析，并对地震发生的时

10、间、地点、强度进与指标，对孕震过程进行追踪分析，并对地震发生的时间、地点、强度进行以物理为基础的概率性预报。行以物理为基础的概率性预报。适合板内地震特点，具有中国特色的地震分析预报科学思路适合板内地震特点，具有中国特色的地震分析预报科学思路具有中国特色的地震分析预报科学思路有以下几个方面具有中国特色的地震分析预报科学思路有以下几个方面 1 1 长、中、短、临渐进式预报思路长、中、短、临渐进式预报思路 2 2 源兆与场兆思想源兆与场兆思想 3 3 源的过程追踪与场的动态监视相结合思想源的过程追踪与场的动态监视相结合思想 4 4 “块、带、源、场、兆、触、震块、带、源、场、兆、触、震”协同的思想协

11、同的思想 长期预报是依据历史地震的统计，对地质构造活动的地壳形变的观测长期预报是依据历史地震的统计，对地质构造活动的地壳形变的观测分析，以及对近代地震活动图像的分析等多方面研究作出的对某地区分析，以及对近代地震活动图像的分析等多方面研究作出的对某地区今后数年至一二十年地震形势的预报。今后数年至一二十年地震形势的预报。 中期预报是根据地震活动图像、地壳介质的物理性质、地壳形变、地中期预报是根据地震活动图像、地壳介质的物理性质、地壳形变、地下水动态、水化学成分、地电阻率、地球磁场、重力场及地壳应力应下水动态、水化学成分、地电阻率、地球磁场、重力场及地壳应力应变等多方面的监测研究，依据多种趋势性异常

12、所作的一至数年的地震变等多方面的监测研究，依据多种趋势性异常所作的一至数年的地震危险区及地震强度的预报。危险区及地震强度的预报。 短期预报是根据趋势异常加速或转折性变化和短期异常的出现等所作短期预报是根据趋势异常加速或转折性变化和短期异常的出现等所作的数月内的地震预报。的数月内的地震预报。 临震预报则是根据突发性快速变化的异常所作的几天至十几天的地震临震预报则是根据突发性快速变化的异常所作的几天至十几天的地震预报。预报。 阶段性地震预报思想就是使预报过程追踪地震孕育过程的发展，以渐阶段性地震预报思想就是使预报过程追踪地震孕育过程的发展，以渐进的方式向未来地震时空强三要素逐步逼近。进的方式向未来

13、地震时空强三要素逐步逼近。 1 长、中、短、临渐进式预报思路长、中、短、临渐进式预报思路 2 源兆与场兆思想源兆与场兆思想 源即震源，源的研究系指对震源形成及演变过程的研究。源即震源，源的研究系指对震源形成及演变过程的研究。 源兆即为在此过程中震源区及近源地区出现的各种效应。源兆即为在此过程中震源区及近源地区出现的各种效应。 场即区域应力场，地质构造块体在边界力作用下形成区域应力场，由于场即区域应力场，地质构造块体在边界力作用下形成区域应力场，由于块体内部不均匀结构，因而在一些特殊部位形成多个应力集中区，其中块体内部不均匀结构，因而在一些特殊部位形成多个应力集中区，其中有的可能发展成为孕震区，

14、有的则为可能反映应力场变化的敏感点。有的可能发展成为孕震区，有的则为可能反映应力场变化的敏感点。 场兆即为在震源形成及演变过程中，大范围区域应力场在众多敏感点显场兆即为在震源形成及演变过程中，大范围区域应力场在众多敏感点显示的异常现象。示的异常现象。 3 3 源的过程追踪与场的动态监视相结合思想源的过程追踪与场的动态监视相结合思想 源的过程追踪思想基于对孕震过程及其可能产生的效应的研究。源的过程追踪思想基于对孕震过程及其可能产生的效应的研究。孕震过孕震过程可分为弹性变形、非弹性变形和破裂加速阶段程可分为弹性变形、非弹性变形和破裂加速阶段。长期预报阶段主要追。长期预报阶段主要追踪弹性应变和应变能

15、的积累过程。中期预报阶段追踪非弹性变形如微破踪弹性应变和应变能的积累过程。中期预报阶段追踪非弹性变形如微破裂发展裂发展( (微破裂数量和线度增加微破裂数量和线度增加) )或扩容等，以及伴随的效应如流体运移或扩容等，以及伴随的效应如流体运移等导致的中期异常发展过程。在短临阶段主要追踪突发性异常，即由于等导致的中期异常发展过程。在短临阶段主要追踪突发性异常，即由于岩体有效强度降低、破裂扩展加速及贯通、断层加速蠕动和不稳定形变岩体有效强度降低、破裂扩展加速及贯通、断层加速蠕动和不稳定形变区内宏观断裂形成等造成的一系列突发性短临异常。区内宏观断裂形成等造成的一系列突发性短临异常。 场的动态监视思想基于

16、中国板内地震具有异常范围较大，异常群体动态场的动态监视思想基于中国板内地震具有异常范围较大，异常群体动态演化过程与上述震源孕震过程同步起伏等基本事实，固而大面积监视场演化过程与上述震源孕震过程同步起伏等基本事实，固而大面积监视场的动态就可以获得震源孕育过程的相关信息或背景性变化。的动态就可以获得震源孕育过程的相关信息或背景性变化。 由于场和源的相互作用，实现地震预报必须将源的过程追踪和场的动态由于场和源的相互作用，实现地震预报必须将源的过程追踪和场的动态监视两者结合起来监视两者结合起来，地震异常在时间上阶段性发展和在空间分布上集中，地震异常在时间上阶段性发展和在空间分布上集中性特点主要反映源的

17、发展演化过程，而前兆现象在空间上和时间上的离性特点主要反映源的发展演化过程，而前兆现象在空间上和时间上的离散性分布则更多地反映孕震过程中场的变化。散性分布则更多地反映孕震过程中场的变化。 4 4 “块、带、源、场、兆、触、震块、带、源、场、兆、触、震”协同的思想协同的思想 “块块”即地震构造块体。即地震构造块体。大陆地壳是由大大小小的不同层次的块体嵌套而成的。大陆地壳是由大大小小的不同层次的块体嵌套而成的。 “带带”为构造块体之间的边界带，亦称构造带为构造块体之间的边界带，亦称构造带。在地球动力因子作用下，地质在地球动力因子作用下，地质构造块体间，出现构造块体间，出现“压、拉、扭、错压、拉、扭

18、、错”多种力学性质的相对运动。边界带是集多种力学性质的相对运动。边界带是集中反映这种运动的剪切带、形变带、应力应变集中带、地球物理和地球化学等中反映这种运动的剪切带、形变带、应力应变集中带、地球物理和地球化学等异常带。异常带。 “源源”即边界带上摩擦强度大的阻挡构造块体运动的地段即边界带上摩擦强度大的阻挡构造块体运动的地段，显然这里将积累应，显然这里将积累应力应变和能量，是可能孕育地震的震源区。力应变和能量，是可能孕育地震的震源区。 “场场”即区域应力场即区域应力场，随着构造活动的持续，应力应变的积累，形成了不断变，随着构造活动的持续，应力应变的积累，形成了不断变化和增强的构造应力场和震源应力

19、场。化和增强的构造应力场和震源应力场。 “兆兆”就是应力场发展过程中形成的反映地震孕育发展过程的异常变化就是应力场发展过程中形成的反映地震孕育发展过程的异常变化。 “触触”是指在孕震晚期震源处于不稳定状态是指在孕震晚期震源处于不稳定状态，外场，外场( (如天体引力、太阳活动、如天体引力、太阳活动、气压场等气压场等) )的某些微小扰动，可能对地震的发生起触发作用。的某些微小扰动，可能对地震的发生起触发作用。 最后，在上述条件统一作用过程中发生地震。最后，在上述条件统一作用过程中发生地震。 三、三、 地震地壳形变地震地壳形变1 1，地震地壳形变概念，地震地壳形变概念： 由于地壳结构的不均匀性，在内

20、力与外力的作用下就会产生不均匀的地壳由于地壳结构的不均匀性，在内力与外力的作用下就会产生不均匀的地壳 形形变，它导致某些地壳特殊部位上的应力变，它导致某些地壳特殊部位上的应力- -应变积累，当这一累积应变达到地壳应变积累，当这一累积应变达到地壳的极限应变值或已有断层上的应力积累达到了断层的抗剪强度时，地壳便突然的极限应变值或已有断层上的应力积累达到了断层的抗剪强度时，地壳便突然破裂，发生地震，破裂，发生地震，这种与地震的孕育和发生过程直接有关的地壳形变过程称为这种与地震的孕育和发生过程直接有关的地壳形变过程称为地震地壳形变。地震地壳形变。震前地壳形变震前地壳形变：孕震过程中的地震地壳形变，也称

21、地壳形变前兆；：孕震过程中的地震地壳形变，也称地壳形变前兆；同震地壳形变同震地壳形变：震时地壳破裂引起的地壳形变；：震时地壳破裂引起的地壳形变；震后地壳形变震后地壳形变：震后继续调整过程的地壳形变：震后继续调整过程的地壳形变。地震地壳形变地震地壳形变2 2，地震地壳形变获取的方式及分类：，地震地壳形变获取的方式及分类： 地震地壳形变可以通过地壳变形观测获得。地震地壳形变可以通过地壳变形观测获得。 地壳变形观测按时间域分类地壳变形观测按时间域分类：连续观测；离散观测：连续观测；离散观测 地壳变形观测按空间域分类地壳变形观测按空间域分类：定点形变：定点形变( (倾斜、固体潮倾斜、固体潮) )观测观

22、测 ； 断层位移观测（局部）断层位移观测（局部） ；区域形变场观测；区域形变场观测 ；全球板块运动监测；全球板块运动监测 地壳变形观测按观测手段分类：地壳变形观测按观测手段分类：有空间技术有空间技术; ;常规大地测量技术常规大地测量技术; ;洞体钻洞体钻 孔应变、倾斜、各类固体潮汐因子等观测孔应变、倾斜、各类固体潮汐因子等观测 。地震地壳形变地震地壳形变3 3，震前地壳形变，震前地壳形变 表现方式：断层系和块体运动；近场地壳应变积累；扩容与地表隆起表现方式：断层系和块体运动；近场地壳应变积累；扩容与地表隆起 3.1 3.1 断层系和块体运动断层系和块体运动 板块内部有地质块体沿断层的相对运动，

23、导致障碍体上的应变积累，从而板块内部有地质块体沿断层的相对运动，导致障碍体上的应变积累，从而形成地震。形成地震。 断层运动及其近场块体的形变受构造运动的驱动断层运动及其近场块体的形变受构造运动的驱动, ,形变主要集中在断层近场形变主要集中在断层近场, ,远场形变比较平缓远场形变比较平缓, ,并且离断层越远形变越平缓。并且离断层越远形变越平缓。 断层近场块体介质的力学特性可以发生很大的变化断层近场块体介质的力学特性可以发生很大的变化, ,这种变化对于块体的变这种变化对于块体的变形与释放能量的方式和速度起着控制作用。形与释放能量的方式和速度起着控制作用。地震地壳形变地震地壳形变3.2 3.2 近场

24、地壳应变积累近场地壳应变积累 通过对岩石做模拟实验可以模拟地震发生的环境，了解地壳的运动模式通过对岩石做模拟实验可以模拟地震发生的环境，了解地壳的运动模式 。 岩石加压模拟实验结果表明，岩石在应力作用下的应变曲线是分段的，即岩石加压模拟实验结果表明，岩石在应力作用下的应变曲线是分段的，即岩石破裂过程中应力、应变关系呈现出岩石破裂过程中应力、应变关系呈现出阶段性特征阶段性特征。 岩石破裂前后应力、应变曲线大致可分为岩石破裂前后应力、应变曲线大致可分为5 5个阶段个阶段 ：硬化阶段硬化阶段 ，线弹性，线弹性阶段阶段 ，非弹性变形或弱化阶段，非弹性变形或弱化阶段 ，失稳阶段，失稳阶段 ，滑动阶段，滑

25、动阶段 岩石差应力岩石差应力应变应变(或位移或位移)曲线示意图曲线示意图轴向应力轴向应力体应变曲线体应变曲线地震地壳形变地震地壳形变3.3 3.3 扩容与地表隆起扩容与地表隆起 在应力达到破裂强度的在应力达到破裂强度的2/32/3左右时，近场区开始扩容，在地表可以观测到左右时，近场区开始扩容，在地表可以观测到隆起，但随着应力的继续增加，主破裂带的逐渐形成，应变将大量集中在一个隆起，但随着应力的继续增加，主破裂带的逐渐形成，应变将大量集中在一个有限的最终出现断裂的部位上。其余地区裂缝将闭合，应力有限的最终出现断裂的部位上。其余地区裂缝将闭合，应力- -应变将下降恢复，应变将下降恢复，膨胀也将逐渐

26、恢复。膨胀也将逐渐恢复。 无震滑动无震滑动: : 地壳表面的大部分位移是由无震滑动产生的地壳表面的大部分位移是由无震滑动产生的( (例：伊豆半岛东北例：伊豆半岛东北部所累积的应变能量，大部分是由无震隆升释放的，以地震形式释放的只占一部所累积的应变能量，大部分是由无震隆升释放的，以地震形式释放的只占一小部分小部分 ) )。4, 4, 同震（震时）地壳形变同震（震时）地壳形变地震的同震运动，主要表现为地震断层和地裂缝地震的同震运动，主要表现为地震断层和地裂缝同震形变模型：断层运动模型。断层是地球内部的一个滑动面，在同震形变模型：断层运动模型。断层是地球内部的一个滑动面，在其两侧发生了不连续的岩体运

27、动（或称为错动）。它在地面下的形其两侧发生了不连续的岩体运动（或称为错动）。它在地面下的形状无法观测到，一般通过地面观测资料（例如断层长度、地震时的状无法观测到，一般通过地面观测资料（例如断层长度、地震时的错动等）来建立断层模型。错动等）来建立断层模型。 断层可分为三类（按两盘相对运动的方向）：正断层、逆断层和平推断层。断层可分为三类（按两盘相对运动的方向）：正断层、逆断层和平推断层。上盘相对下降、下盘相对上升的断层称正断层，断层面倾角一般较陡。上盘相对下降、下盘相对上升的断层称正断层，断层面倾角一般较陡。上盘相对上升、下盘相对下降的断层是逆断层，断层面倾角变化较大，从陡上盘相对上升、下盘相对

28、下降的断层是逆断层，断层面倾角变化较大，从陡倾到近水平。倾到近水平。如果断层两侧的岩石不是沿断层面上下移动而是沿水平方向移动，则称平推如果断层两侧的岩石不是沿断层面上下移动而是沿水平方向移动，则称平推断层。断层。如果把这三类断层与形成的构造应力联系起来，通俗地说，正断层由拉张应如果把这三类断层与形成的构造应力联系起来，通俗地说，正断层由拉张应力引起，逆断层是挤压应力的结果（故常造成地壳的缩短），平推断层则与力引起，逆断层是挤压应力的结果（故常造成地壳的缩短），平推断层则与剪切应力有关，其断层面常近直立。剪切应力有关，其断层面常近直立。 4.1 4.1 同震变形的位错模型同震变形的位错模型位错模

29、式几何图位错模式几何图 duuuuFukjkikjinnijkji1在各向同性介质中，表面在各向同性介质中，表面 由于位错由于位错 产生的位移场产生的位移场 如下：如下： ),(321ju),(321xxxui321,UUU分别表示断层上点的走滑位错，倾滑位错，张裂位错分量分别表示断层上点的走滑位错，倾滑位错，张裂位错分量 L L表示断层的半长度，表示断层的半长度，W W为断层宽度，为断层宽度，d d为断层的深度，为断层的深度，为断层的倾角为断层的倾角表示垂直平面的矢量表示垂直平面的矢量kjk克罗内克（克罗内克（Kronecker）符号）符号断层（震源）坐标系(o-xyz)和测量坐标系(O-X

30、YZ)的位置关系 设地面一点在测量坐标系中的位移为 ，相应的在震源坐标系的位移为 , 则有： ),(ZYXUUU),(zyxuuuzyxZYXuuuuuu1000cossin0sincos321,|,UUUbaLWdOYXUUZYXZYX正演问题：用断层的正演问题：用断层的1010个参数计算地表测量坐标系中的位移个参数计算地表测量坐标系中的位移反演问题：地表位移观测值来确定断层的三个位错分量及其位置参数和几何参数反演问题：地表位移观测值来确定断层的三个位错分量及其位置参数和几何参数 4.2 Okada4.2 Okada（19851985）关于矩形位错的公式）关于矩形位错的公式 |sinsin)

31、(2|sincos)(2|sintan)(24121111IRqRRqdUuIRqRRqyUuIqRRRqUuzyx走滑分量 ),(),(),(),(),(WpLxfpLxfWpLxfpLxff各分量分别代表点在矩形四个位置上的位移， sincosdypyx,d为断层中心坐标值（在震源坐标系中）， 为断层倾角， 为断层深度 由位错引起的应变 走滑分量 sincossin2sincos2sincos2sin2sin2432233122312322331121JRRRqAqRqyUyuJAqRqUxuJAqRdUyuJqAUxuyyxx5, 5, 震后的地壳形变震后的地壳形变地壳的震后运动，一般和

32、震前的运动趋势成反向。地壳的震后运动，一般和震前的运动趋势成反向。震后滑动主要发生在主震后几天至几个月内。震后滑动主要发生在主震后几天至几个月内。其特点是滑动速率成对数衰减，位错量可等于或超过同震滑动。其特点是滑动速率成对数衰减，位错量可等于或超过同震滑动。 震后的地壳形变模型震后的地壳形变模型Pollitz最初在1992年给出基于球体分层粘弹性地球模型的时变震后应力释放的解析表达式。该模型主要计算在球体分层弹性-粘弹性介质内，某个弹性层内发生的断层滑动产生的应力变化所引起的震后响应，这种响应对应着球形运动分量和环形运动分量的球谐展开，每种分量代表一种具有特性的时间延迟和空间变形模式的应力释放

33、，它可以灵活计算在任何深度的时间序列震后变形场-位移场和应变场。lTmlmSmlrSrSrS),(),(),()()(iiiiFEtrutru,),(总的位移场 利用上述参数就可以灵活计算任意点在任何深度的时变震后变形场-位移场和应变场 四、四、GPSGPS、InSARInSAR在地震研究中应用在地震研究中应用 Example: Sumatra SequenceFrom Jeff Freymueller2004 Sumatra-AndamanFrom Jeff FreymuellerGPS Timeseries Reveal Coseismic SignalsFrom Jeff Freymue

34、llerTwo clear offsets in SSIA timeseries: January 13: 15+/-3 mm February 13: 39+/-3 mmSSIA- El Salvador (closest to the seismic activity)MANA-Nicaragua (furthest from the seismic activity)Can we use the GPS data to learnsomething about the earthquakes?2004 Sumatra-AndamanFrom Jeff FreymuellerSteps t

35、o Study Collect post-earthquake data Identify displacements in time series How to estimate step? Deformation between pre-earthquake and post-earthquake survey Postseismic deformation What to do with displacements? Elastic dislocation model(s) to relate slip on fault(s) to displacementsParkfield Eart

36、hquakeFrom Jeff FreymuellerSeparating Coseismic and PostseismicFrom Jeff FreymuellerParkfield eqJ. Langbein, USGSLessons from ParkfieldFrom Jeff Freymueller Postseismic deformation appears to begin immediately after earthquake Coseismic and postseismic mechanisms are different Almost all “coseismic”

37、 estimates include postseismic deformationInSAR vs/+ GPSFrom Jeff Freymueller InSAR does not require pre- earthquake field survey Spatially dense measurements But only one component, and the impact of the look angle makes for strange “look” Minimum 30 day repeat intervalThe 2002 Denali EarthquakeFro

38、m Jeff FreymuellerThe EarthquakeFrom Jeff Freymueller November 3, 2002, about 1:30 pm local time Mostly ruptured the Denali fault Also Susitna Glacier fault (thrust) and Totschunda fault (strike-slip) Initial estimated magnitude MW 7.9 Preceded by MW 6.7 on October 23 Preceded by ML 4.5 foreshockTec

39、tonic settingFrom Jeff FreymuellerMainshock and AftershocksFrom Jeff FreymuellerEpicentral RegionFrom Jeff FreymuellerCoseismic GPS SitesFrom Jeff FreymuellerCoseismic Displacements at 232 GPS SitesEstimation of DisplacementsFrom Jeff Freymueller Continuous GPS Average of 4 days before and after ear

40、thquake Campaign GPS with well-determined velocities Linear fit to pre-earthquake data plus displacement For sites with imprecise velocity, we used a model or interpolation of surrounding sites Campaign GPS with no post-earthquake measurements until summer 2003 Postseismic correction applied based o

41、n interpolation of postseismic displacements of surrounding sites.Coseismic Displacements - HorizontalFrom Jeff FreymuellerTwo scales (3:1)Coseismic Displacements - VerticalFrom Jeff FreymuellerModeled Vertical Displacement Field (from GPS data)Richardson Hwy Profile (and the TAPS corridor)From Jeff

42、 FreymuellerCoseismic Slip ModelFrom Jeff FreymuellerCoseismic GPS SitesFrom Jeff FreymuellerComparison With Surface Offset DataFrom Jeff FreymuellerComparison With Surface Offset DataFrom Jeff FreymuellerSAR Range Offset DataFrom Jeff FreymuellerUnfilteredFiltered, baddata masked outModel Compariso

43、nFrom Jeff FreymuellerFrom Jeff Freymueller Decaying curves typical of right-lateral strike-slip motion Displacement in LOS across the fault is 4 m, which roughly translates to 6 m of right-lateral slip Speckle tracking data agrees well with predictions from GPS model, particularly near the fault wh

44、ere model resolution is highestFrom Jeff Freymueller Total moment of 6.85*1020 Nm, corresponds to an Mw7.89 earthquake Average of 4.5 m of dip-slip on the Susitna Glacier fault On Denali fault, no surface slip around epicenter Three patches of high slip (8.5 m - 10 m) between 60 and 140 km Largest s

45、lip patch 40 km west of DTJ with slip reaching 10 m at surface On Totschunda fault, relatively low slip concentrated at shallow depthsCombined slip distribution ModelSurface slip comparisonFrom Jeff Freymueller五五 地壳形变监测台网与地震预报地壳形变监测台网与地震预报 地壳形变测量能精确定量地监测到地震发生前后地壳的一系列运动、变形地壳形变测量能精确定量地监测到地震发生前后地壳的一系列运

46、动、变形（位移、速度、加速度、应变、倾斜、蠕滑）、重力和介质物性（密度、（位移、速度、加速度、应变、倾斜、蠕滑）、重力和介质物性（密度、勒夫数）的空间分布及其随时间变化。地壳形变测量和测震一样，被科学勒夫数）的空间分布及其随时间变化。地壳形变测量和测震一样，被科学界公认为是力学型的最直接的地震前兆，是地震预报的基础。界公认为是力学型的最直接的地震前兆，是地震预报的基础。我国的地壳形变测量，特别是地壳形变监测台网经过近我国的地壳形变测量，特别是地壳形变监测台网经过近4040多年的努力，逐多年的努力，逐步形成了具有特色的点、线、面和空、地、深相结合的布局，积累了大量步形成了具有特色的点、线、面和空、地、深相结合的布局，积累了大量的观测资料，成为地球动力学和地震预测预报必不可缺的基础与支柱。的观测资料，成为地球动力学和地震预测预报必不可缺的基础与支柱。地壳形变监测台网正