大地测量联合反演理论及应用
联合反演模型的研究：赵少荣(1991)系统地研究了“基于固体力学的大地测量反演问题”,给出了基于固体力学模式的多类大地测量反演的解算模型；许才军(1994)给出了大地测量联合反演构造应力场的解算模型; 晁定波等(1997)进一步提出了四维整体大地测量有限单元法,把四维整体大地测量模型(包括GPS,水准,GPS水准和重力监测数据)与固体力学基本方程结合起来,对大地测量和地球物理数据进行有效的整体处理,强化了边界的大地测量约束,具有改善刚度矩阵方程数学性质和降低其阶数的优点,从而提高反演解的稳定性和可靠性。P. Segall等(1997)根据永久GPS网的扩建在时空上为地壳形变量提供了足够的数据,建立了一种网络反演滤波模型,它可以综合频繁收集到的各类大地测量网络数据来估计断层滑动的时空分布,反演获得各种参数,包括观测误差、局部移动、瞬时和空间平滑参数。大地测量反演模型的建立涉及到正演模型和观测数据的类型、反演模型选择及反演参数辨识。正演模型是反演模型的基础,只有清楚正演问题建立正演模型,才能给出反演模型。大地测量反演已由单一的观测数据的反演,发展到多种观测数据的大地测量联合反演问题,这种多种观测数据的大地测量联合反演问题不仅仅指大地测量的多种观测数据,而且也包括不同类型的(地震、地质和大地测量等)多种观测数据。
联合变形、重力和地震资料反演研究地球内部介质参数：地球内部介质的弹性参数μ和λ是利用地表的形变观测资料研究地球的应力场及内部的物质运动所必需的参数。观测资料表明,在地球内部,介质的弹性参数μ和λ存在横向和纵向上的不均匀性。（1）利用变形资料可以研究地球介质弹性参数μ和λ；2）利用地震和重力观测资料可以研究参数μ和λ的不均匀性。由于参数μ、λ与介质的密度、横波速度和纵波速度参数有关,而密度、横波速度和纵波速度参数又是地震观测资料的场源,密度参数还是重力数据的场源,因此,进一步探求参数μ和λ结构的问题,可以转化为利用地震和重力观测资料反演密度、横波速度和纵波速度参数结构的问题。联合变形、重力和地震资料研究地球内部介质参数的联合反演问题可以表示为:
BX = UL NX = g SX = d （1）式(1)可以通过下列目标函数取极小值求解:Φ= ( UL-BX)TWL( UL-BX)+(g-NX)T
? Wg(g-NX)+(d- SX)TWd(d-SX)=min。该式中,X为介质参数； WL, Wg, Wd为3类数据的权阵,观测值可以是相关的。
反演模型与参数的辨识：大地测量联合反演模型有多种多样,但无采用何种模型都是在某一准则下使实际观测数据与反演理论数据达到最佳拟合的条件下确定反演模型的参数。在大地测量联合反演中,模型的选择及模型参数的可确定性,可定义为大地测量联合反演的模型辨识。吕爱钟等(1998)在理论上讨论了反演问题的参数可辨识问题;独知行(2001)基于简单的力学模型具体讨论了反演参数的可辨识性。但大地测量联合反演问题可能是一个十分复杂的问题,它的研究对象可以是全球,也可以是地球的一个特定局部,可以具有复杂的材料性质、边界条件及形状,在进行反演研究时,很难简化成有解析表达式的模型。反演模型与参数的可辨识性问题仍需要进一步的研究。
联合反演算法的研究：大地测量反演方法可分为解析法与数值法两大类。由于解析法只适用于线性问题和简单的非线性问题的反演,因此,大地测量反演方法的研究主要在数值算法方面,其中数值算法尤以优化反演方法(即把反演参数的估计通过建立目标函数转化为一个优化问题进行求解)的研究最为普遍,进展也最快。大地测量随机反演方法,如蒙特卡罗法、模拟退火法、遗传算法等,特别是遗传算法在地球物理大地测量的反演解算中得到了广泛的应用。随机反演方法充分利用模型空间参数的优化及先验信息对反演参数增强约束,改善反演解的稳定性。在利用随机方法获取全局最优解的基础上,可利用基于贝叶斯理论的高斯_牛顿法通过迭代计算获取全局最优解,这种联合方法既考虑了随机方法可以在全局解域内搜索全局解的特点,同时又考虑了基于先验信息利用贝叶斯方法可以获取局部最优解的特点,从而提高了反演计算的效率和可靠性(党亚民,1998)。最近IA算法(区间优化算法)(L.G.Casado等,2001;J.Fernandez 等,2001)、数值流形方法———“将连续体的有限单元算法、非连续变形分析方法DDA和解析方法统一起来的一种更高层次最新数值分析方法”也被提出用于大地测量反演解算中(许才军等,2000)。这些将有助于全面探索中国大陆地壳运动及动力学问题,如数值流形方法可以采用连续和非连续覆盖函数的办法把连续和非连续的力学问题的计算统一起来,可更有效地研究中国大陆地壳运动问题。
联合大地测量资料、活断层运动速率和地震矩张量反演地壳运动速度场：
Haines和Holt (1993)提出利用地震矩张量获得的地壳应变率反演地壳运动速度场的方法。目前,该方法已被发展为利用大地测量、活断层运动速率和地震矩张量资料综合反演地壳运动速度场(Haines等,1995;Jackson等,1995;T.B Shen 等,1995;Tinnon等,1995;Kreemer等,2000;Holt等,2000;Beavan等,2001)。该方法的主要思想是利用“双三次样条函数”通过反演方法在一定厚度的板壳上拟合构造应变率分布,从而求出连续的地壳运动速度场。应变率场和速度场的拟合可通过下列目标函数求解:
∑cells∑αβ,λμ(.εfitαβ-.εobsαβ)V- 1αβ,λμ(.εfitλμ-.εobsλμ)+ ∑points∑α,β(vfitα- vobsα)V-1α,β(vfitβ- vobsβ) =min。该式中,Vαβ,λμ和Vα,β分别是平均应变率和大地测量观测速度的方差_协方差阵。α,β,λ,μ分别是经度、纬度的起止值。T.B.Shen等(1998)应用此方法,利用地震矩张量、甚长基线干涉测量和GPS观测值反演研究了美国西部现今地壳运动速度场;Holt等(2000)应用此方法,利用第四纪以来断层滑动速率和GPS观测值反演研究了亚洲地壳运动速度场;许才军等(2001)也利用此方法采用GPS和地震矩张量反演计算了中国大陆地壳运动速度场,给出了中国大陆地壳运动速度场图像。
联合变形、地震、地质资料反演研究地壳应变_应力场、应变能密度场：由GPS、形变观测或地震矩张量可获得地壳应变率场。由弹性体地壳应力与应变的关系可得到地壳应变能密度如下(C.Straub,1996)：DE=12trace(ζ,ε) =μ((ε
ji+λΘ2μδji)εij)应变能密度DE的单位是J/m3。地壳中积累的应变能密度可能是地震释放能量的来源,是影响地震活动的重要因素(C.Straub,1996;张东宁等,1999)。应变能密度是衡量地壳积累能量高低的重要数量指标,应变能密度越大,地壳积累能量越高,今后发生地震的可能性越大。C.Straub(1996)利用GPS 反演计算了Mar-mara Sea地区的地壳应变_应力场、应变能密度场进行地震危险性区域估计;许才军等(2000,2001)利用华北地区GPS监测网多期观测数据所得的速度场在椭球面上进行了应变率计算,得到了其最大剪切应变率、主应变率的图像;并利用1966年以来25次大于5.0级地震的地震矩张量,由Kostrov公式计算了地震应变率,对两种数据所得的应变率进行了比较分析,结合地质构造背景探讨了华北地区现今形变特征及地震危险性地区,在青藏高原及其邻区也作了类
似的工作;许才军等(2001)进一步利用GPS和地震资料联合反演计算了中国大陆地壳应变_应力场、应变能密度场,并结合地震、地质构造背景可以进行中国大陆地震危险性区域分析
联合卫星测高、声纳测深及海洋激光雷达遥感信息反演海洋环境参量：开发研究海洋,是解决人口膨胀、资源匮乏和环境污染等人类目前所面临的一系列问题的主要途径之一。合理开发利用海洋,离不开基础海洋测绘成果———海底地形图的保障,更主要的是要研究海洋环境参量。海洋基础测绘的传统方法是首先布设海底控制网,再利用船载测深仪以一定的路线进行测深,对获取的数据进行改正、平差后推算海底地形,绘制海底地形图。随着科学技术的发展,传统的基础海洋测绘方法被赋予了新的内容,特别是全球定位系统(GPS)、多波束测深系统(MES)、侧扫声纳(SS)、机载激光测深系统(ALB)和海洋地理信息系统(MGIS)等一系列新技术和新方法的产生给海洋科学研究带来了一场技术革命。这些新技术和手段的应用,提高了海洋空间地理信息量的精度,增加了海洋空间地理信息的数量,也提高了获取海洋空间地理信息的速率。卫星测高利用雷达测量卫星和海面之间所获得的距离、量测时的位置及各种海面高改正来提供海面高。利用卫星测高资料可以恢复海洋重力场,在此基础上,利用重力异常和海底地形存在线性关系的均衡假设,便可推估海底地形。（1）W.H.F.Smith等(1994)利用源于卫星测高数据获得的重力资料,采用最优化理论和经验资料反演确定了南部大洋的海深；（2)希腊的D. Arabelos(1997)运用Moho面深度(12′×12′)、JGP95E全球水准数据库(5′×5′)、卫星测高资料和重力异常资料,应用最小二乘配置法反演确定了地中海和北大西洋的海底地形；(3) G.S.Vergos (2001)在IAG大会上介绍了利用卫星测高获得的重力数据通过联合反演方法改善海底地形估计的结果。随着CHAMP、GRACE和GOCE 3颗关于重力场、地磁场、大气探测和海洋环流的小卫星的发射,联合多种卫星测高资料反演研究海底地形将取得重大突破;卫星测高资料也将为反演研究海洋环境参量作出贡献。海洋环境参量主要包括海面温度、潮汐、海面风场、重力场、海浪方向谱、有效波高、海冰厚度、海面高度等,利用海洋航空遥感遥测技术所获得的信息反演研究海洋环境参量是我国海洋“863”资源环境领域海洋监测技术主题“十五”重点研究课题。
联合INSAR、GPS和水准资料反演研究地壳运动：INSAR利用雷达信号的强
度信息和相位信息获取地球表面的三维空间信息。由于雷达信号能够穿透云雾、雨雪,具有全天候的工作特点,对地物有一定的穿透能力,并且可大面积地测定地面点的精确高程及其变化等优越性,因而该技术得到广泛应用。INSAR技术可以获取高精度(精度可达5m)的数字高程模型(DEM),补充遥远、条件恶劣地区的高程数据;差分INSAR技术能够探测到cm级甚至mm级的地表形变。由于雷达干涉的高空间覆盖能力及相对低廉的成本,近年来日益受到重视,在地震、火山及地壳形变监测、DEM提取、土地资源利用调查等领域建立了独特的地位。雷达干涉已用于多例地震,还有冰川运动和冰川地形研究,热带雨林监测、滑坡监测及地表沉降观测等。自身亦存在限制,例如INSAR对大气影响、卫星轨道误差、地表状况和时变去相关性非常敏感,这些误差很容易导致INSAR图像解释错误,而且这些都无法用SAR数据本身消除。
GPS精密定位可以确定地面离散点上的精确位置和高程变化,可以较为精确地确定电离层、对流层参数,是当前研究地壳形变最为精确、方便、实用的手段。INSAR与GPS两种技术,具有很好的互补性,可以相互结合,取长补短。GPS是一种理想的点定位系统,尤其是采用相对定位工作方式时,定位精度已达10-9~10-10,而IN-SAR提供的是整个区域面上连续的信息。GPS获得的是高精度的绝对坐标,而INSAR仅提供相对坐标。由于入射角的关系,INSAR对高程信息特别敏感,而这恰恰是GPS最薄弱的一环。最为重要的是,GPS允许长时间连续观测,而INSAR可被看作瞬时测量。GPS可提供时间分辨率很高的观测数据(采样率为10Hz乃至20Hz),而SAR卫星通常35d左右的重复周期很难提供足够的时间分辨率(ERS_1和ERS_2共同服役期间重复周期较短)。因此,GPS和INSAR互为补充,可以同时保证时间连续性和空间连续性,特别是利用GPS改善INSAR卫星轨道的确定精度以及减少大气效应的影响,进而改善INSAR图像质量是一种行之有效的手段。
1997年Bock和Williams首先提出INSAR与GPS集成的思想。同年,Linlin Ge等人进一步提出了DIDP(doubleinterpolation and double prediction )方法。他们提出利用GPS数据减弱大气影响(对流层延迟、电离层延迟等)和卫星轨道误差,并通过改正后的INSAR图像去加密GPS观测结果。李振洪在2001年IAG 大会上作了题为“Subsidence of Tianjin City Measured by Integration of Interfero-metric SAR and GPS Measurements”的报告,报告提出了GPS与INSAR
两种集成方案,并提出联合INSAR、GPS和水准资料研究天津地面沉降运动的思想。联合INSAR、GPS和水准资料反演研究地壳运动,预测预报地质灾害也将是国际地学界的一个热点。
大地水准面是地球内部物质结构与运动的物理特性的一种几何表征,它与地球深部构造有着密切的关系,因此利用大地水准面起伏和重力异常的精细结构可以反演地球内部构造、地幔对流。此外利用重力、水准、定位、地磁等精确观测数据,可以探测地球内部局部构造细节,地壳上地幔的结构、地表和岩石圈的密度、莫霍面深度和应力场,为地球物理探矿、油气寻找、地热调查等服务。
在过去的几十年中,中外学者已在这方面作了部分工作,取得了初步成果。联合各种观测数据反演研究地球内部构造是地球科学中一项长期的、艰巨的基础任务,它将为探索地球奥秘、发展地球科学理论作出重要贡献。
浅谈地形变观测的发展与地震预报 地震预报是一个系统工程，需要长期的科学探索和科学积累，地形变观测系统已成为地震预报学科领域不可取代的力学型的基础观测系统。本文简述了地形变观测的发展历程以及其在地震预报中的运用。 标签：地震预报地形变测量 我国地域辽阔，地质构造复杂，是世界上地质灾害较严重的国家之一。地震就是对人们生活和生命、财产安全影响较大经常发生的一种灾害。地震的前兆是可以预见的。如果在地震前能获知一些前兆特征，便能提醒人们引起注意，从而减少地震造成的灾害。 地震预报是对未来破坏性地震发生的时间、地点和震级及地震影响的预测。 地震预报目前或今后很长一段时间都将是以观测为主的试验性科学。物理统计分析方法是研究地震孕育动力学过程或进行预测的最具有潜力的途径之一。地形变测量就是物理统计分析方法的其中一种。尽管2008年汶川和2011年东日本大地震预报失败，但这两次和其他大地震前后GPS和其他观测得到的地壳形变表明，大地震是有前兆的，是可以预报的. 地震孕育的能量来源于地壳运动过程中产生的应变能积累，地形变观测是监测地壳运动与变形、认识地震孕育过程与开展地震预测的重要手段。如2008年汶川地震后，通过对其地形变化监测，监测显示，汶川地震引起震中区域监测点的水平位移量达238厘米，沉降量达到70厘米，隆起量达30厘米。通过数据观测，得出结果：地震造成灾区地形发生重大变化地形发生重大变化。根据这个结论，开展地形变分析后可以有效预报大地震后的余震及其震级。 地形变测量是指对一个地区地面的相对变化进行的重复地形变测量或连续观测。一般是在某个区域布设相当数量的形变观测点，在区域外设立基准点，利用常规大地测量仪器和工程测量中监测变形的方法。通过定期观测，并通过适当的数据处理和物理解释准确确定区域内各测点的空间随时间变化的“绝对”位移与方向，以达到预测地震之目的。地形变观测目前已具有多种手段并形成一定规模。地震地形变观测系统目前所涵盖的主要观测手段有：垂直形变测量网、水平形变测量网即GPS区域复测、重力测量网、跨断层形变测量网、重力与固体潮观测台网、地倾斜与固体潮观测台网、地应变与固体潮观测台网和连续GPS监测网也叫GPS基准网等。所采用的手段和仪器主要有：水准测量、三角测量、倾斜仪、伸缩仪、电阻丝应变仪、激光测距仪、测潮仪；近些年来还发展到应用一系列空间对地观测技术，如GPS（全球空间定位系统）等。 地形变观测系统所测定的物理量包括：位移、旋转、速度，加速度、应变（应力）、蠕变、位错、重力、固体潮汐、地下介质物性参量（密度、勒夫数）、电离层与对流层介质物性参量（电子浓度、湿度…）等的空间分布及其随时间变化。
大地测量联合反演理论及应用 联合反演模型的研究：赵少荣(1991)系统地研究了“基于固体力学的大地测量反演问题”,给出了基于固体力学模式的多类大地测量反演的解算模型；许才军(1994)给出了大地测量联合反演构造应力场的解算模型; 晁定波等(1997)进一步提出了四维整体大地测量有限单元法,把四维整体大地测量模型(包括GPS,水准,GPS水准和重力监测数据)与固体力学基本方程结合起来,对大地测量和地球物理数据进行有效的整体处理,强化了边界的大地测量约束,具有改善刚度矩阵方程数学性质和降低其阶数的优点,从而提高反演解的稳定性和可靠性。P. Segall等(1997)根据永久GPS网的扩建在时空上为地壳形变量提供了足够的数据,建立了一种网络反演滤波模型,它可以综合频繁收集到的各类大地测量网络数据来估计断层滑动的时空分布,反演获得各种参数,包括观测误差、局部移动、瞬时和空间平滑参数。大地测量反演模型的建立涉及到正演模型和观测数据的类型、反演模型选择及反演参数辨识。正演模型是反演模型的基础,只有清楚正演问题建立正演模型,才能给出反演模型。大地测量反演已由单一的观测数据的反演,发展到多种观测数据的大地测量联合反演问题,这种多种观测数据的大地测量联合反演问题不仅仅指大地测量的多种观测数据,而且也包括不同类型的(地震、地质和大地测量等)多种观测数据。 联合变形、重力和地震资料反演研究地球内部介质参数：地球内部介质的弹性参数μ和λ是利用地表的形变观测资料研究地球的应力场及内部的物质运动所必需的参数。观测资料表明,在地球内部,介质的弹性参数μ和λ存在横向和纵向上的不均匀性。（1）利用变形资料可以研究地球介质弹性参数μ和λ；2）利用地震和重力观测资料可以研究参数μ和λ的不均匀性。由于参数μ、λ与介质的密度、横波速度和纵波速度参数有关,而密度、横波速度和纵波速度参数又是地震观测资料的场源,密度参数还是重力数据的场源,因此,进一步探求参数μ和λ结构的问题,可以转化为利用地震和重力观测资料反演密度、横波速度和纵波速度参数结构的问题。联合变形、重力和地震资料研究地球内部介质参数的联合反演问题可以表示为: BX = UL NX = g SX = d （1）式(1)可以通过下列目标函数取极小值求解:Φ= ( UL-BX)TWL( UL-BX)+(g-NX)T
?全球问题：①地球动力现象引起的地震海啸火山等自然灾害，给人类生命财产带来损失②全球气候变暖、海平面上升、局部地层沉降和海上溢油公害等是随着工业发展引起的环境问题③由于人口增加和陆地资源枯竭，需开拓生存空间和寻找新矿产资源。 ?地壳运动：指在地球内部构造应力作用下所引起的地壳一些元素的相对运动。它们可以是⊥运动、水平运动或地倾斜运动，综合表现为大面积的地壳形变。广：地质旋回、狭：构造旋回。分类：水平运动指组成地壳的岩层，沿∥地球表面方向的运动，又称造山运动或褶皱运动；⊥运动又称升降运动或造陆运动，表现为岩层部分区域的隆起和相邻区域的下降。长期运动是在地质时间尺度内的运动，由几千年到几百万年，它与板块运动有关；瞬变运动与地震和火山等活动相联系的。范围：全球板块运动和区域及局部地壳运动。地壳运动监测：测定板块运动参数、大陆板块和海洋板块的内部形变、板块边界与大地震有关的区域形变和局部形变、其他地震活动区的区域形变和局部形变，主要通过建立全球测定板块运动监测网和区域、局部地壳运动监测网实现。 ?地壳形变：指在地球内力和外力作用下，地球的地壳表面产生的升降、倾斜、错动等现象及其相应的变化量。成因：①人类活动，地表形变：离散性、短暂性和局部性； ②地球自转和极移，形变：全球规模特性，可理论计算；③日月等天体对固体地球在引力作用下，形变：固体潮,理论上可严格计算；④大地构造运动（地球内部的构造原因），地壳构造形变：连续性、长期性、区域性、复杂性。地壳形变测量：对一个地区的地壳表面的相对变化进行重复或连续的观测称为地壳形变测量。特点：①以动态观测替代大地测量只以静态观测方式来测定地面点变化，并分析研究其物理意思；②主要在构造带、多震区和具有潜在地震危险的重点地区及在大坝等要害部位进行，而大地测量未考虑这些；③测点设置要求稳定可靠，布网边长短、测量精度高、复测周期密。任务：监测地壳形变运动，具体观测元素是地表点位置变化。种类:①全球板块运动监测：主要用来测定板块运动参数，测定大陆板块和海洋板块的内部形变，VLBI、SLR 、GPS。②全国及区域地壳形变测量：精密水准、高精度流动重力、高精度空间大地。③断层形变测量：短水准、短基线、短边GPS网等。④定点形变测量：有效地监测地壳的连续变动，地倾斜、地应变、重力台站。 全球板块运动监测：意义：①寻找矿产资源，板块运动的边界是生成矿源的地点；②防灾减灾，地震一般发生在板块的边界。方法：VLBI测量、SLR测量、GPS测量等。 ?VLBI：可提供整体运动和地壳运动的丰富信息，短时间测量即可得极高精度。特点：目标源和参考源（点源）距离近；两种源时间间隔小。①纯几何方法，不涉及地球重力场；②不受气候限制，有长期稳定性；③提供以河外射电源为参考的最佳准惯性参考系。?SLR：测定台站的位置变化或站间基线长度的变化率来确定板块运动参数。在一列选定的时间间隔内联合求解测站的坐标、卫星轨道和地球定向参数EOP→测站间基线的时间序列→线性拟合→基线长度变化率。不足：①由于采用可见光，受天气影响大，不能全天候观测；②SLR台站建立和维护费用高。精度：取决于在各所选时间间隔内测得的站间基线精度、基线的观测误差；EOP误差以及卫星的定轨误差均有影响。?GNSS WGS84：世界大地坐标系统，地心系，GPS广播星历。WGS84（G1150）：美国对WGS84第3次精化获得的框架。G框架用GPS资料确定，1150开始用的时间。 ?空间大地测量方法：建立现代板块运动模型。观测值：站坐标和站速度（与参考框架有关）→绝对运动参数，站间基线长度变化率（无关）→相对~。原理：板块构造学说认为相邻两板块之间的相对运动实际上是围绕通过地球中心的一个轴的旋转运动，通常用欧拉定理来表述V=ω×r。相对运动参数：，k、l板块上i、 j测站间基线长度的变化 率，kl ω是k对l的相对运动角速度，R i 、R j：i、j的坐标矢量，测两板块间若干ij B →由上式加权的最小二乘平差得kl ω。绝对运动参数：Ω(ω x ，ωy，ωz)角速度矢量，r(x,y,z)位置矢量，λ、φ经纬度。欧拉定理V=ω×r→地心系→地球近似为球体→站心系→角速度ω、纬度Φ、经度Λ。 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? - - - = ? ? ? ? ? ? ? ? ? ? z y x G z y x x y x z y z V V V ω ω ω ，? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? - - - = ? ? ? ? ? ? ? ? ? ? z y x G z y x r r r r r r V V V ω ω ω λ ? λ ? λ ? ? λ ? ? cos cos sin cos cos cos sin sin cos sin ， ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? - - - = ? ? ? ? ? ? z y x r n e r r r r r V V ω ω ω ? λ ? λ λ ? λ cos cos sin sin sin sin cos ，。 区域地壳形变监测：区域空间尺度由几百到1000km的瞬变运动，GPS/INSAR/GRAVITY+VLBI/SLR；局部几百米到几十公里+LEVELING。 ?2000国家GPS大地控制网：2000国家重力、GPS大地网。缺陷：密度<全国天文大地网；所提供的低密度的3维地心坐标框架不能完整实现中国的3维地心系。解决：~GPS 网+全国天文大地网联合平差，将后者纳入3维地心系，并提高其精度和现势性，使我国的大地坐标框架在密度和分布方面实现我国3维地心大地坐标系的跨越式进步。CGCS2000通过~大地控制网的坐标和速度具体实现。参考历元2000.0。 ?GPS用于地壳⊥形变监测：①有可能利用GPS观测直接得到毫米级的大地高数据。重复GPS观测可以求定大地高的变化（或站心系U分量的变化）。②由于椭球体法线与该点夹角很小，大地高和正常高方向基本重合，∴可用大地高的变化代替正常高的变化，即可利用重复GPS观测取代精密水准以监测地面的升降变化。u分量与正常高变化：?h=?u/cos αAB。∵αAB=1°(S AB=110km)→cos αAB ≈1∴S AB ≈100km，可用?u代替?h。 ?InSAR：合成孔径雷达干涉测量一种以合成孔径雷达天线记录到的回波信号为信息源，利用干涉测量获取地球表面的三维地形、地表形变和地物特征变化等信息的测量技术。处理流程：配准、生成干涉图、平地相位消除、去地形影响、干涉图滤波、相位解缠。星载INSAR优点：①全球范围②成像测绘带宽③全天候④轨道稳定，多次成像重复性好⑤飞行平稳，无运动补偿问题⑥下视角小，图像阴影小。局限：①轨道固定，侦察
第一章绪论 1. 地壳运动与变形的基本内涵是什么？ 地壳运动：是在地球内部构造应力的作用下，所引起的地壳一些元素的相对运动，它可以是垂直运动，水平运动或地倾泻运动。综合表现为大面积的地壳变形。 地壳形变：是指在地球内力和外力作用下，地球的地壳表面产生的升降、倾斜、错动等现象及其相应的变化量。 2. 你对地壳运动与变形的具体表现形式有哪些认识？ 其表现形式主要有：断裂，地震，火山，滑坡，泥石流，地面裂缝与沉降，冰川运动。 3. 现代大地测量学在地球动力学研究中可以发挥哪些作用？ 4. 简述地球物理大地测量学的发展。 第二章地壳形变测量 1. 地壳形变测量有哪几类？各有什么特点？ （1）全球板块运动监测 主要用来测定板块运动参数，测定大陆板块和海洋板块的内部形变，其观测手段主要采用VLBI、SLR 和GPS等空间测量技术。 （2）全国及区域地壳形变测量 测定亚板块及构造块体的地壳形变，给出全国大陆动力学的边界条件，以及全国大陆应力场、形变场变化过程的总体和分区特征；区域地壳形变测量主要测定块体边界与大地震有关的区域形变，它可以给出大陆内部地形变的时空演变图象。 （3）断层形变测量 在各活动构造块体边界上进行的近场构造变形测量。能够直接测定块体边界断裂及其不同段落的现今活动方式、相对位移速率以及它们随时间变化的过程，提供震间、震前、同震与震后滑动等构造活动的微动态信息。 目前以短水准、短基线、短边GPS网以及由水管倾斜仪、伸缩仪、蠕变仪、短边激光测距仪或重力仪组成的台阵等为主要手段。 （4）定点形变测量 主要包括地倾斜、地应变和重力(固体潮汐)台站观测。这种方法可以有效地监测地壳的连续变动，可以通过不同时间间隔的采样，在相当宽的频带范围内对地壳动力学现象进行观测。 2. GPS连续跟踪站测量属于定点形变测量吗?为什么? 属于 第三章：地球参考系统与地球参考框架 1. 什么是惯性参考系？ 其时间是均匀流逝的，空间是均匀和各向同性的；在这样的参考系内，描述运动的方程有着最简单的形式。这样的参考系就是惯性系 牛顿运动定律成立的参考系，称为惯性参考系，简称惯性系。 2. 我国的大地测量坐标基准有哪些？ 1954年北京坐标系；1980年西安坐标系（1980年国家大地坐标系）；2000国家大地坐标系。 3. 大地坐标系建立的基本原理是什么？ 4. 全球最优的协议地球参考架CTRF是如何建立的？ 由IERS组织，遵循以下原则①原点位于包括海洋和大气在内的整个地球的质心；②长度国际单位制为米，尺度为广义相对论下的局域地球框架；③定向为最初由国际时间局（BIH）所给出1984.0定向；④定向的时变通过一个关于全球的水平构造运动的非净旋转条件（地壳无整体旋转）来保证。
一一 第４５卷一第６期测一绘一学一报 V o l ．４５,N o ．６ 一２０１６年６月A c t aG e o d a e t i c ae tC a r t o g r a p h i c aS i n i c a J u n e ,２０１６ 引文格式:X U K e k e ．T h e o r y a n d M e t h o d f o rD e t e c t i n g A b n o r m a l C r u s t a l D e f o r m a t i o nU s i n g G N S S [J ]．A c t aG e o d a e t i c a e t C a r t o g r a p h i c a S i n i c a ,２０１６,４５(６):７５６．(徐克科．G N S S 地壳形变异常检测理论与方法[J ]．测绘学报,２０１６,４５(６):７５６．)D O I :１０．１１９４７/j ．A G C S ．２０１６．２０１６０１５２ G N S S 地壳形变异常检测理论与方法 徐克科１, ２１．河南理工大学测绘与国土信息工程学院,河南焦作４５４０００;２．同济大学测绘与地理信息学院,上海２０００９２ T h e o r y a n dM e t h o d f o r D e t e c t i n g A b n o r m a l C r u s t a l D e f o r m a t i o nU s i n g G N S S X UK e k e １, ２１．S c h o o l o f S u r v e y i n g a n dL a n d I n f o r m a t i o nE n g i n e e r i n g o fH e n a nP o l y t e c h n i cU n i v e r s i t y ,J i a o z u o ４５４０００,C h i n a ;２．C o l l e g eo f S u r v e y i n g a n dG e o ＧI n f o r m a t i o no f T o n g j i U n i v e r s i t y ,S h a n g h a i ,２０００９２,C h i n a 一一地震的孕育过程是一个长期缓慢的过程, 地震发生时断层破裂所释放的能量只是一部分,其中有很大一部分能量在地震前后以无震蠕滑的形式释放.因此,检测这些形变异常信息对于地震危险性评估至关重要.随着G N S S 监测网络数据的持续积累,时空分辨率越来越高,就有可能从中挖掘出更有价值的形变信息.本文利用G N S S 时空数据, 从高精度数据处理二全球板块运动模型构建二多尺度速度场与应变场估计二地壳微动态形变异常检测和断层滑移时空反演等方面,研究了G N S S 地壳形变异常检测理论与方法.主要研究内容与结果如下: (１)重建了G N S S 测站速度估计模型:法方程重构模型二基线向量最小二乘模型二卡尔曼滤波模型和坐标时序模型.前３种模型将坐标二速度二年二半年周期项作为参数,在进行网平差的同时一并求解速度值.坐标时序模型引入白噪声和幂律噪声组合,分析了区域周期项振幅与相位空间分布的一致性.基于４种模型,估计了川滇陆态网G N S S 测站速度.结果差异在１mm /a 内,精度在亚毫米级.从而验证了这些速度估计模型的一致性. (２ )提出了基于统计假设检验和稳健估计的全球板块运动模型构建方法.利用I T R F ２００８V E L 求取的全球板块欧拉参数与先前模型具有较好的一致性.建立了相对欧亚板块背景场的中国大陆速度场,分析了中国大陆现今地壳运动特征. (３)构建了G N S S 多尺度速度场和应变场估计模型. 利用负位错模拟数据作了大量试验,结果表明,不同尺度的应变场具有检测不同空间范围的地壳形变的能力.利用２００９ ２０１１年陆态网数据估计并分析了中国大陆３ ８尺度的应变场及形变特征. (４ )联合卡尔曼滤波和主成分时空分析,集时空滤波与形变检测于一体,构建了瞬态无震蠕滑时空检测模型,进一步提高了数据时空信噪比.利用滇西南陆态网 G N S S 数据,检测其时空分布特征与２０１１年缅甸Mw ７．２级地震相对应.得出了滇西南区域的断层活动可能会受到缅甸地震带的影响. (５)提出了基于G N S S 基线面应变和G N S S 网形的 震前形变异常检测方法.得出２０１３年芦山M s ７．０地震和日本近年来的４次地震在震前数月的时间内均有不同程度的异常偏离.尤其G N S S 基线夹角二基线方位角和第一剪应变的异常变化更为突出.推断震前可能产生了强烈的左旋剪切构造应力变化,加速了芦山地震的孕育发生. (６ )以地壳形变检测与断层滑移反演为一体,构建了G N S S 主成分和卡尔曼滤波时空反演模型.通过模拟试验,得出了正确反演断层滑移时空分布所需要的最低信噪比和最优的台站分布密度.以２００５年苏门答腊Mw ８．７地震震后余滑和２００６年墨西哥慢滑移为例,反演了断层蠕滑时空分布. 中图分类号:P ２２８一一一一文献标识码:D 文章编号:１００１Ｇ１５９５(２０１６)０６Ｇ０７５６Ｇ０１基金项目:国家９７３项目(２０１３C B ７３３３０４ );国家自然科学基金(４１４０４０２３ )收稿日期:２０１６Ｇ０４Ｇ０５ 作者简介:徐克科(１９７９ ),男,副教授,２０１５年７月毕业于同济大学,获工学博士学位(指导教师:伍吉仓教授),研究方向为G N S S 数据处理与地壳形变分析. A u t h o r :X U K e k e (１９７９ ),m a l e ,r e c e i v e dh i sd o c t o r a l d e g r e e f r o mT o n g j iU n i v e r s i t y o nJ u l y ２０１５(P h Da d v i s o r :P r o f ．W U J i c a n g ),m a j o r si n G N S S d a t a p r o c e s s i n g a n d c r u s t a l d e f o r m a t i o na n a l y s i s E Ｇm a i l :１２x k k ＠t o n g j i ．e d u ．c n
GPS地壳形变与地震预测 自20世纪80年代以来,VLBI、GPS、SLR等空间大地测量技术的精度得到大幅度提高,其观测结果已被应用于板块运动、冰后期反弹、火山、地震、地球自转和地球系统内部物质再分布等地球动力学过程的研究[1,2]。尤其是GPS观测技术,由于GPS接收机价格不是特别昂贵,观测比较容易实施,在地壳形变监测中得到了广泛应用。如为了监测地壳形变,日本在境内布设了1 000多个连续GPS观测站点,平均每隔30 km就有一个GPS观测站点[3];美国为了监测南加州圣安德烈斯断层带的运动和形变,布设了超过250个站的连续GPS监测网[4]。我国自20世纪80年代开始,在青藏高原、南北地震带、川滇和华北首都圈等地壳运动活跃区和地震监测重点区布设了一系列GPS观测网,用于研究这些区域的地壳形变和地震以及判断板块内部块体构造运动模型[5]。这些GPS区域网中的部分站点已经并入国家连续GPS观测网或定期复测GPS观测网,列入长期观测计划。随着连续和分期GPS网观测的持续,积累的地表运动观测数据越来越多。这是一类观测精度非常高(上百km范围内可以测量1mm/a的水平方向变化)、分布在监测区域地表上的站点坐标时间序列。利用这类时空观测数据提取精确可靠的地壳形变信息,反演地球动力学机制是现代大地形变测量数据处理的科学目标,为此,很多学者开展了一系列的理论和方法研究,也取得了不少应用成果。 从现有的研究成果来看,通过全球连续GPS网观测数据已经成功地解算了全球板块运动模型参数,其结果与基于地质和地球物理调查数据(百万年平均)得到的板块运动模型参数(NUVEL-1A)非常接近，连续或分期的GPS观测网成功地观测到了同震位移和震后位移[6],并且用于反演地震断层面上的位错分布，或结合地震波数据联合反演,增加人们对地震破裂过程的认识，通过GPS观测结果结合地质调查和地球物理勘探结果,研究划分板块内部次一级的构造运动块体,以及确定各块体之间的相对运动等[7]。目前,国内外研究关注的焦点和难点是如何通过GPS观测得到的站点坐标时间序列提取块体边界或断层带的地壳形变信息,进而研究如地震孕震形变等地球动力学过程[8]。
第30卷第9期 岩 土 力 学 V ol.30 No.9 2009年9月 Rock and Soil Mechanics Sept. 2009 收稿日期：2008-03-13 基金项目：中国科学院知识创新工程重要方向项目资助（No.KZCX2-VW-134）；国家自然科学重大研究项目重点基金（No.90814014）资助；973课题（2004CB418405）资助；地震基本科研业务项目（020*******）资助。 第一作者简介：张晁军，男，1965年生，博士，高级工程师，主要从事慢地震及计算地球动力学的研究。E-mail: zhangchaojun@https://www.doczj.com/doc/3818898088.html, 文章编号：1000－7598 (2009) 09－2552－07 昆仑山大地震震后形变反映的地壳岩石流变特性 张晁军1， 2，石耀霖2，马 丽3 （1. 中国地震台网中心 信息网络部，北京 100045；2. 中国科学院研究生院地球科学学院 计算地球动力学实验室，北京 100049； 3. 中国地震局地震预测研究所，北京 100036） 摘 要：2001年昆仑山地震是我国近50年来最大的地震。用3种模型对震后中国地震局跨断层布设的4个GPS 站点记录到显著的震后形变进行了模拟研究。结果表明：单纯的上地壳10 km 为滞弹性的模型不能解释震后形变的幅度；30 km 的弹性上地壳覆盖在 40 km 的柔性下地壳上的松弛模型可以解释变形速率指数衰减的主要特征；而两个模型的结合不但能解释整 体指数衰减的特征，而且还能更好地拟合震后最初几周的较高形变速率。结果表明，在昆仑山断层两侧存在着流变性质的 差异。 关 键 词：震后形变；黏滞系数；流变模型；滞弹性；标准线性固体；Maxwell 体 中图分类号：P 554 文献标识码：A Numerical simulation of crust rheological property reflected by post-seiemic deformations of Kunlun large earthquake ZHANG Chao-jun 1, 2 , SHI Yao-lin 2, MA Li 3 (1. China Earthquake Network Center, Beijing 100045, China; 2. Laboratory of Computational Geodynamics, Graduate School, Chinese Academy of Sciences, Beijing 100049, China; 3. Institute of Earthquake Science, China Earthquake Administration, Beijing 100036, China) Abstract: The Mw7.8 Kunlun earthquake of 14 November, 2001, in the northern Tibetan Plateau of China, is the largest event in the Chinese continent in the latest 50 years. In this paper, layered viscoelastic models are calculated by PSGRN/PSCMP code; and the results are used to fit post-seismic deformation obtained from four GPS stations of CEA across the earthquake fault. The results show that model of a single anelastic layer of 10 km thick upper crust overlying an elastic lower crust cannot explain the amplitude of deformation; the relaxation model of 30 km elastic upper crust overriding 40 km ductile lower crust can explain the main features of exponential attenuation of post-seismic deformation. Combination of the two models, however, can fit the character better with deformation rate being higher in the first few weeks and slower rate thereafter. The viscous layer of lower crust provide the main control of the post-seismic deformation of long term decay of deformation in months, while the upper anelastic layer may contribute to the high rate of initial few weeks after the main earthquake. The result also suggests that rheological differences exist at two sides of the Kunlun Mountains fault. Key words: post-seismic deformation; viscosity rheological model; anelasticity standard line solid; maxwell body 1 前 言 2001年11月14日17时26分，在昆仑山口地区发生了7.8级强烈地震（美国国家地震信息中心NEIC 地震目录）。震后中国地震局利用GPS 开展了相对密集的地壳形变观测。穿越震区沿青藏公路跨断层布设了4个固定GPS 站点：格尔木（JB30）、昆仑山（KLGD ）、不冻泉（BDGD ）和五道梁 （WDGD ），其中不冻泉和五道梁位于东昆仑断裂南侧，距地震地表破裂带距离分别为30 km 和 130 km ；昆仑山和格尔木位于东昆仑断裂北侧，分别距地震地表破裂28 km 和100 km （图1）。从 JB30、KLGD 、BDGD 和WDGD 站点震后地表蠕动位移东向（E ）运动的时间序列[1]可以看出，在断层南侧虽然WDGD 和BDGD 两者与断层距离不同，东向运动幅度与时间变化特征仍很相似，变形速率
基于GPS数据的几种地壳形变分析方法 张永奇 长安大学地质工程和测绘学院，西安（710054） 摘要：主要介绍了几种在地壳形变分析中的方法，介绍了他们的原理，以及在地壳形变监测分析中的应用，同时也说明了几种方法的优缺点。 关键字：DDA；位错模型；刚体旋转加均匀应变模型；经验正交函数分解方法模型；数值流形法 引言 地震的孕育和发生在本质上是地壳内部的应变能逐渐积累并突然释放的结果。伴随着大震孕育或应变能的显著积累, 岩石圈表层必然会表现出某种形式和量级的地壳形变。基于这样的认识, 地壳形变监测一向是地震监测或地震危险性分析的最重要手段之一[1]。 自20 世纪80 年代以来,VLBI、GPS、SL R 等空间大地测量技术的精度得到大幅度提高,其观测结果已被应用于板块运动、冰后期反弹、火山、地震、地球自转和地球系统内部物质再分布等地球动力学过程的研究。尤其是GPS 观测技术,由于GPS 接收机价格不是特别昂贵,观测比较容易实施,在地壳形变监测中得到了广泛应用。我国自20 世纪80 年代开始,在青藏高原、南北地震带、川滇和华北首都圈等地壳运动活跃区和地震监测重点区布设了一系列GPS 观测网,用于研究这些区域的地壳形变和地震以及判断板块内部块体构造运动模型。这些GPS 区域网中的部分站点已经并入国家连续GPS 观测网或定期复测GPS 观测网,列入长期观测计划。随着连续和分期GPS 网观测的持续,积累的地表运动观测数据越来越多。这是一类观测精度非常高(上百km 范围内可以测量1mm/ a 的水平方向变化) 、分布在监测区域地表上的站点坐标时间序列。利用这类时空观测数据提取精确可靠的地壳形变信息,反演地球动力学机制是现代大地形变测量数据处理的科学目标,为此,很多学者开展了一系列的理论和方法研究,也取得了不少应用成果[3]。 进一步的基于GPS观测数据的地壳运动和形变的分析方法，尚处于研究阶段。目前初步的应用方法大体可以分为两类：1.纯数学方法，包括多面函数拟合法，多项式拟合法，双三次样条函数法，协方差经验函数法等；2.以物理力学为基础的方法，如欧拉矢量法，非连续形变分析法（discontinuous deformation analysis ，DDA）等。前者尽管在数学上是自洽的，但物理力学含有不清，且都是以地壳连续性假设为基础的；后者物理力学意义清楚，且考虑了地壳的非连续性，但描述块体形变不够精细和有效。如DDA对块体形变的描述仅用3个参数，过于简单只适合与小尺度块体系统的变形分析。 本文主要介绍现今在地壳形变分析的一些常用方法，如非连续变形分析方法，位错运动模型，刚体旋转加均匀应变模型，经验正交函数分解方法模型，数值流形法。以及对这些方法的优劣进行比较，在此基础上对一些可能应用在地壳形变分析的新方法进行展望，如小波理论，灰色理论，神经网络等。 1.非连续变形分析(DDA)方法 模拟介质不连续缝的历史可追溯到30年前的Goodman、Taylor和Brekke等教授发展的节理单元。对岩土裂缝的数值计算发展很快,并已在岩石工程中得到广泛应用。Cundall介绍的离散元法现在被广泛应用于节理或块状岩石。两者是用虚拟力来调整滑动和阻止块体重叠的一种方法,有时候可达到稳定。1993年,由石根华博士发展的非连续变形分析是全新的数值