球形地球模型的地震位错理论及其应用＊

孙文科 付广裕 周 新 徐长仪 唐 河董 杰 周江存 杨君妍 王武星 刘 泰

1）中国北京 100049 中国科学院大学地球与行星科学学院

2）中国北京 100083 中国地质大学（北京）地球物理与信息技术学院

3）中国北京 100085 应急管理部国家自然灾害防治研究院

4）中国北京 100029 中国科学院地质与地球物理研究所

5）中国北京 100036 中国测绘科学研究院

6）中国武汉 430071 中国科学院精密测量科学与技术创新研究院

7）中国北京 100036 中国地震局地震预测研究所

引言

地震是地球内部应力场作用下，断层发生摩擦、失稳及破裂的物理过程（Scholz，1998）.从Reid （1910）提出弹性回跳理论、Love （1911）建立经典地球弹性力学到Steketee （1958）将位错理论引入地震学，地球物理学家们逐渐建立了弹性动力学理论（地震波动理论），并将该理论用于解释地球自由震荡和地震波等资料（如Altermanet al，1959；Gilbert，MacDonald，1960；Ben-Menahem，Toksöz，1962；Aki，1964）.与此同时，用于断层破裂完成后静态变形场研究的理论—准静态地震位错理论也逐渐被建立起来，进一步促进了地震循环过程（震间—同震—震后）的研究（如Savage，Burford，1973；Savage，Prescott，1978；Savage，1983）.

早期的准静态地震位错理论假设地球为无限半空间弹性体介质模型，即忽略了地球曲率和介质不均匀性，得到了断层错动引起的地表变形解析式（如Mindlin，Cheng，1950；Chinnery，1961；Maruyama，1964；Weertman，Weertman，1964；Press，1965）.Okada （1985，1992）整理了前人的工作，给出了均匀半无限空间模型中位移场及其空间导数的简洁表达式.在此基础上，Okubo （1991，1992）给出了同震重力和位场变化的解析解.此外，有些研究人员考虑了地球的层状结构但忽略了地球曲率效应，通过汤姆森-哈斯克尔（Thomson-Haskell）传播矩阵方法（Thomson，1950；Haskell，1953）获得了半无限空间的地表变形响应（Ben-Menahem，Singh，1968；Jovanovichet al，1974a，b；Rundle，1980；Roth，1990；Ma，Kusznir，1992；Wanget al，2003，2006）.目前以Okada （1985，1992）、Okubo （1991，1992）和Wang 等（2006）为代表的研究人员所提出的半无限空间介质模型的位错理论已经发展得非常成熟，该理论可被用于计算地震产生的同震和震后变形，包括位移、倾斜、应变、大地水准面和重力变化等，因此广泛用于解释大地测量观测的近场地壳变形，对地震周期变形过程的研究起到了重要作用（如，Rolandoneet al，2006；Amoruso，Crescentini，2009；Daiet al，2014；Jianget al，2014；周硕愚等，2017；Pan，2019）.这些位错理论具有数学上的简洁性和解析性，但是由于半无限空间介质模型与真实地球差异较大，忽略了地球的几何曲率和层状结构效应，其计算结果存在较大的模型误差，不适于研究地震造成的大区域和全球尺度的地震变形、地球动力学变化等科学问题（Sun，Okubo，2002；孙文科，2012b；Donget al，2014）.

另一方面，传统大地测量手段（水准、三角网、三边网等）记录了地震引起的地壳变形（Savage，Hastie，1969）.以准静态位错理论为基础，通过反演这些观测数据可以确定发震断层的几何参数和滑动分布，以理解地震的震源过程（Savage，Hastie，1969；陈运泰等，1975，1979）.然而，现代大地测量技术，特别是甚长基线干涉测量（very long baseline interferometry，缩写为VLBI）、全球导航卫星系统（Global Navigation Satellite Systems，缩写为GNSS）、合成孔径雷达干涉（interferometric synthetic aperture radar，缩写为InSAR）、卫星测高、高精度卫星重力测量（包括：challenging mini-satellite payload，缩写为CHAMP；gravity recovery and climate experiment，缩写为GRACE；gravity field and steady-state ocean circulation explorer，缩写 为GOCE；GRACE follow-on，缩写为GFO）等观测技术的出现和发展，使传统大地测量发生了革命性变化，进而催生了动力大地测量学研究，使我们能够以更高的精度和全球尺度研究地球内部结构、动力学过程和地震变形问题，有力地推动了地球科学的发展（许厚泽，王广运，1989；孙文科，1989，2002；胡明城，2000；陈俊勇，2003；Marotta，2003；Chao，2003；李建成等，2006；周硕愚等，2008；Krynski，2012；Bock，Melgar，2016；Fernándezet al，2017；Pepe，Calò，2017；Braitenberg，2018；Giorgiet al，2019；Larson，2019；姚宜斌等，2020；Freymueller，2021；Wiesemeyer，Nothnagel，2021）.现代观测技术克服了传统大地测量方法固有的局限性，具有全球、全天候、实时和连续的观测能力，为科学研究提供了规则、稠密和高质量的测量数据（Wdowinski，Eriksson，2009；Salviet al，2012；Liet al，2015；Zhouet al，2015；Spaans，Hooper，2016；Zhu，2016；姚宜斌等，2020）.特别是，现代大地测量技术可以观测地震周期中不同阶段的变形，弥补了地震学与构造地质学之间的观测空白，逐步成为地壳运动及地震周期变形研究中不可或缺的观测手段（Segall，Davis，1997；Allen，Ziv，2011；Matsuo，Heki，2011；Glennieet al，2013；Meigs，2013；Herringet al，2016；Van Campet al，2017；Biggs，Wright，2020；Wang，Shen，2020）.对于上述大空间尺度的地震周期变形问题，基于半无限空间模型的地震位错理论不再适用，需要一个更适合的理论，如基于球形地球模型的地震位错理论，以开展震源机制、内部结构与断层反演、观测数据解释等科学问题的精细研究.

随着观测技术的不断更新，基于分层球形地球模型的地震位错理论经国内外多个研究团队的研究得以不断发展.例如，考虑地球曲率和径向不均匀性，在球坐标下可以采用传播矩阵（Thomson，1950；Haskell，1953）或龙格-库塔方法（Gilbert，Backus，1966，1968）得到地表同震变形（Wason，Singh，1972；Israel，Ben-Menahem，1974；Dahlen，1974）.但这些工作均对地球模型进行了不同程度的简化，如忽略了地球的固态内核结构、自引力和重力效应及介质的可压缩性等.Sun （1992）以及Sun和Okubo （1993）首次基于球对称、非自转、理想弹性和各向同性（spherically symmetric non-rotating elastic isotropic，缩写为SNREI）地球模型（Dahlen，1968），考虑了介质的自引力和重力、可压缩性等多种因素，解决了同震重力变化和同震位移等问题，建立了球形地球模型位错理论的基本理论框架.同一时期，美国以Pollitz团队为代表，其对自重效应予以近似处理，用简正模方法独立研究了同震和震后变形问题（Pollitz，1992，1996，1997）；欧洲以Sabadini，Vermeersen和Spada团队为代表，他们也使用简正模方法研究了类似的问题，但其早期的工作并未同时考虑介质的可压缩和自重效应（Sabadiniet al，1984；Spadaet al，1992；Piersantiet al，1995；Vermeersenet al，1996；Cambiottiet al，2009）.需要指出的是，在以上研究组中孙文科团队最早针对包含固态内核-液态外核-固体地幔分层结构的地球模型，以严谨的理论同时考虑了介质的可压缩性、自引力和重力效应，提出了球形地球弹性位错理论的基本框架，并经过其学生和合作者得以不断发展和完善，最终构建了一个新的球形地球模型位错理论体系.该团队提出的基于球形地球模型的地震位错理论，开拓了全球地震变形及地震引起的动力学变化研究的新方向，扩展了现代大地测量在地震研究中的应用范围，深化了对地震破裂和震后变形机制的认识，取得了一系列创新性成果，主要包括以下两部分：

1）建立了一套新的地震位错理论，推进了全球地震变形研究.定义了地震位错勒夫（Love）数和格林（Green）函数，构建了解决球形地球模型的地震变形问题的基本理论框架；发展了基于层状弹性地球模型、三维不均匀弹性地球模型和层状黏弹性地球模型的地震位错理论，并给出了不同变形物理量的渐近解；解决了同震和震后重力变化、位移、应变、库仑应力变化及其三维效应等问题.

2）拓展了上述理论的应用范围，开拓了地震引起的动力学变化研究的新方向.发展了基于上述理论的断层滑动、位错勒夫数、黏弹性结构的反演方法；估计了GRACE卫星对地震重力变化的检测潜力；建立了同震垂线偏差变化理论，拓展了GRACE的应用范围；提出了地震引起的地球膨胀和质心偏移理论；定量分析了全球地震产生的地球自转、体积、势能等动力学因子变化；推进了全球地震变形以及相关地球动力学变化研究.

该团队关于地震位错理论和相关应用目前共发表了近100篇科学论文，这些论文分散于众多文献中.为了使广大科研人员，特别是地震科研人员，能够较全面地了解该理论，本文将就该位错理论成果作一个报道性介绍.需要特别声明的是，本文不是对地震位错理论的国内外发展历史进行全方位总结和回顾，而是针对本团队所建立的球形地球模型地震位错理论的研究成果加以概括性介绍，以期阐明本理论的基本面貌及相关应用的最新进展，并对位错理论的发展方向作出展望.

本文首先简要介绍球形位错理论的创建过程，包括弹性位错理论的提出、三维横向非均匀理论的创建、黏弹性位错理论的发展和完善；然后介绍地震位错理论的应用，包括全球地震变形场的解释、断层模型和地下结构的反演以及对空间大地观测技术的应用拓展等；最后对位错理论的发展方向进行简要讨论.

1 球形地球模型地震位错理论的发展

地震位错理论是研究震源机制、反演断层、解释观测数据等的必要理论基础，在地震学、地球动力学、大地测量学中占有核心地位.地震位错理论需要基于一定的地球模型来建立，它相当于地震震源与地震变形之间的系统函数，从理论上给出地球在地震力源作用下的变形响应.由于位错理论的计算精度取决于地球模型，地球模型的准确与否将直接影响理论计算精度及其适用范围.所以，基于合理的地球模型来建立位错理论具有重要科学意义.

为精确地计算同震和震后变形、合理地解释全球大地测量观测数据、更准确地反演断层滑动分布，本团队在研究中采用了球形地球模型，并据此建立了地震位错理论新体系.新位错理论包含地球曲率、层状结构效应、黏弹性效应、三维非均匀效应等，在物理上更严谨、更合理，同时适用范围更宽泛.

而建立球形地球模型的位错理论需要采用与半无限空间位错理论不同的研究思路和数学方法，需要在球坐标系中求解地震变形的偏微分方程组，解决积分变换、级数收敛、积分稳定性等一系列数学问题.根据观测技术以及理论本身的发展，本研究团队先后针对一维分层弹性球形地球模型、三维不均匀弹性球形地球模型、一维分层黏弹性球形地球模型建立了不同的理论分支，形成了较完善的地震位错理论体系.

1.1 弹性球形地球模型的地震位错理论

球形地球模型的地震位错理论首先采用一维分层球形弹性地球模型，即球对称、不旋转、弹性和横向各向同性的地球模型（SNREI模型，Dahlen，1968）.SNREI模型是一大类地球模型，其中1066A模型 （Gilbert，Dziewonski，1975）、PREM模型（Dziewonski，Anderson，1981）和IASP91模型（Kennett，Engdahl，1991）等被广泛使用.基于1066A和PREM模型，孙文科等提出了弹性球形地球模型的地震位错理论.采用半解析半数值方法求解地震变形的微分方程组，定义了地震位错勒夫数，克服了勒夫数计算、求和、截断和加速收敛等数值计算困难，给出了四个独立点震源的格林函数（Sun，1992，2003，2004a，b；Sun，Okubo，1993，1998；Sunet al，2006a，2009），并在此基础上进一步给出了计算任意有限断层在任意位置产生变形的数值积分算法（Sun，Okubo，1998）.利用该理论可以解释实际地震产生的位移、应变、重力位、重力等的同震变化，为地表同震变形解释、断层破裂反演等研究奠定了基础，目前也已被广泛使用.更重要的是，这些研究搭建了球形地球模型位错理论的理论框架，为后续发展三维地球变形和黏弹性地球变形理论奠定了基础.

团队随后对地表破裂源、地震造成的地球动力学变化等问题进行了研究，扩展了地震位错理论的内涵及其适用范围.例如：孙文科和董杰借助互易定理（Okubo，1993）推导了地表破裂源产生的地表变形的地震位错勒夫数和格林函数，解决了地表震源的奇异解问题（Sun，Dong，2013）；徐长仪和周江存等给出了地震造成的地球体积变化、低阶重力场系数变化、地球自转变化和重力势能变化的计算公式（Zhouet al，2013，2014a，2016；Xu，Sun，2014；Xuet al，2014；Xu，Chao，2017；Xu，2021）；孙文科、徐长仪和周江存等给出了特大地震造成的地心移动解（Sun，Dong，2014；Xu，Chao，2015；Zhouet al，2015）.近几年周江存等提出了计算位错勒夫数的近似解析解方法，提高了计算效率，同时解决了内部格林函数的收敛问题（Zhouet al，2019a，b，2020）.

迄今为止，上述弹性球形地球模型的地震位错理论已经得到了广泛的应用.利用地震变形格林函数，Sun和Okubo （1998）成功地解释了1964年阿拉斯加MW9.2大地震的同震重力变化，其理论计算结果与观测的重力变化基本吻合，表明新理论可以合理地解释现代大地测量技术观测到的同震变形.此外，为了便于国内科研人员使用该理论，付广裕和孙文科（2012b）撰文详细地介绍了球体位错理论计算程序的总体设计与具体实现过程，并开始在科学社区分享其Fortran程序，得到了学界的认可.更重要的是，该理论因其物理模型的合理性而优于半无限空间位错理论，为地球物理学、大地测量学和地球动力学等应用提供了全新的、可靠的理论保证（孙文科，2012b）.该理论已经成功应用于地震断层破裂反演、地震解释以及火山产生的变形，并且可以用来解释同震变化（Okubo，2020）.例如，Imanishi等（2004）利用我们团队的球形地球位错理论合理解释了超导重力仪记录到的2003年日本十胜冲（Tokachi-Oki）MW8.3地震产生的亚微伽级重力变化，该研究结果发表在国际著名期刊Science上，属于世界上首次报道亚微伽级同震重力变化.

2012年孙文科系统地总结了弹性地震位错理论与部分应用的研究成果，撰写了专著《地震位错理论》（孙文科，2012a），这是目前国内外关于球形地球模型位错理论的唯一专著.

1.2 三维不均匀弹性地球模型的地震位错理论

随着现代大地测量观测精度的日益提高，真实三维地球结构对地震变形的影响逐步被观测到（Huet al，2004；Suito，Freymueller，2009；Freedet al，2017；Diaoet al，2018；Tianet al，2021），然而上述基于球对称地球模型的地震位错理论却无法体现局部的三维效应.为了提升地震变形的理论支撑作用、提高对地震变形观测数据的解析能力，在上述一维球形模型的地震位错理论基础上，付广裕等提出了三维不均匀地球模型的地震位错理论（Fu，Sun，2007，2008，2009；Fuet al，2010；付广裕，孙文科，2012a）.

基于微扰方法的三维地球模型的位错理论，提高了地震变形和重力潮汐的计算精度.具体地说，Fu和Sun （2007）推导了密度的横向增量对重力潮汐影响的计算公式，弥补了Molodenskiy （1980）理论的不足，并计算了半日潮重力潮汐因子的全球分布，为三维地球模型的地震变形研究提供了参考.之后，就三维地球模型的地震变形问题，将实际三维不均匀地球分解成球对称地球模型和对应的横向不均匀增量，分别进行计算得到球对称解和三维响应，之后叠加得到最终变形（Fu，Sun，2007，2008，2009；Fuet al，2010）.球对称解可以直接利用上述球对称地球模型位错理论（Sun，1992）计算得到.三维响应即地球的横向不均匀结构对同震物理场变化的影响可以分为震源的响应和地球横向不均匀结构的响应，可分别由震源函数的扰动及平衡方程式的变分来求解.经过理论推导，给出了六个特殊点源位错引起的位移、重力变化等计算公式（一个垂直走滑位错，两个相互垂直的倾滑位错，三个引张位错）.对这些公式进行线性组合即可计算任意位置、任意类型位错产生的同震变化.依据36阶P波速度模型，Fu和Sun （2008，2009）利用岩石试验经验关系式推导出三维S波速度模型和密度模型，并计算出三种典型类型的点源位错产生的同震重力变化，其结果显示三维响应与位错类型、震源深度均相关，其最大响应约占球对称解的0.5%，且S波速度模型的影响最大.数值结果同时表明：三维响应中震源的响应与地球横向不均匀结构的响应处于同一量级，地球横向不均匀结构对同震位移的贡献达到球对称解的1%—2% （Fuet al，2010；付广裕，孙文科，2012a）.

上述研究成果得到了国内外同行的普遍认可.发表在Science上的Ito和Simons （2011）一文将Fu和Sun （2007）的成果视为潮汐理论研究领域的重要进展加以介绍.国际大地测量联合会2015年度的Bomford奖获得者Yoshiyuki Tanaka博士在多篇论文中引用上述研究（Tanakaet al，2006，2007，2009），将其视为位错理论研究领域的前沿进展加以介绍.在孙文科教授的专著 《地震位错理论》 序言中，许厚泽院士用一句话将上述研究总结如下：“······利用变分方法研究了极为困难的三维地球模型的同震变形问题” .这均说明三维不均匀地球模型的地震位错理论虽然难度很大，但是它在实际应用中具有重要意义和价值.

1.3 黏弹地球模型的地震位错理论

地球介质除具有弹性性质之外，其在较长时间尺度上还具有显著的流变特性（Takeuchi，Hasegawa，1965；Anderson，O’Connell，1967；汪汉胜等，1997；Peltier，2021）；地震变形除了呈现明显的同震变形外，还表现出显著的震后变形（Nur，Mavko，1974；Freed，Bürgmann，2004；Wang，2007；Chenet al，2011）.大地震的震后变形已经可以被现代大地测量技术清晰地观测到，例如，GNSS、InSAR和GRACE均可检测到震后信号（Panetet al，2007；Ryderet al，2007；Wanget al，2012；刘泰等，2017；Cambiotti，2020）.为了解释震后变形观测数据，需要建立基于黏弹性地球模型的地震位错理论.

为此，很多学者针对黏弹性地球模型加以研究，提出了相应的计算方法，例如简正模叠加、回路积分和Post-Widder等方法（Pollitz，1997；Wang，1999；Tanakaet al，2006，2007；Meliniet al，2008）.然而，黏弹地球模型的地震位错理论的发展具有一定的挑战性，对于复杂的黏弹性结构地球模型，始终存在数值计算困难.例如：简正模方法无法处理地球的径向连续结构或介质的可压缩性；回路积分方法需要针对特定黏弹性模型，通过数值搜索方法确定积分路径；Post-Widder的实现需要扩展精度函数库.总之，简正模分布复杂、积分核数值震荡、收敛慢、计算量大等问题均未得到完美解决（唐河，孙文科，2021）.为此，唐河和孙文科基于上述弹性位错理论的基本思想，发展和完善了基于黏弹地球模型的地震位错理论（Tang，Sun，2019；Tanget al，2020a，b）.

黏弹性变形理论的核心是在复数域计算复勒夫数，并进行逆拉普拉斯（Laplace）变换以得到时域勒夫数.首先将黏弹性地震变形方程组进行拉普拉斯变换至复数域，转换为拉普拉斯域内的等效弹性问题，从而可以利用上述弹性位错理论的处理方案加以解决，得到复域的地震位错勒夫数和格林函数.该步骤所使用的计算方法已在弹性理论中得到充分验证，可以确保计算的准确性.之后，在将复域解变换回时间域解时，需要处理逆拉普拉斯数值积分，但积分核的震荡特性造成了严重的数值积分困难，极具挑战性.为此，我们对逆拉普拉斯积分核作分式函数近似和级数展开（Valsa，Brančik，1998；Tang，Sun，2019；Tanget al，2020a，b），得到一个交错级数，从而使得时域勒夫数的计算表达为沿虚轴采样的复勒夫数之和.本方法采用的逆拉普拉斯积分方法也已在电气工程等多个领域得到充分验证，具有很高的计算精度.这样，我们构建了一种新的计算黏弹性地震变形的虚轴积分方法，这是黏弹地震变形理论的最新方法.该方法对于任意线性黏弹性地球模型和任意力源类型均适用，且具有一致的计算精度，其计算方法实现简单、收敛快，不需要扩展精度函数库.因此，该方法非常适于黏弹性结构的反演研究.此外，本团队还发展了黏弹性均质球模型的震后变形解析解.例如，Tang和Sun （2018a，b）针对有自重和无自重两种可压缩黏弹均质地球模型，基于对应性原理和解析逆拉普拉斯变换，首次得到了时间域内位错勒夫数和格林函数的解析形式.解析的地震变形格林函数阐明了介质参数与变形物理场之间明确的函数关系，同时也适用于近场地震变形的理论计算.上述理论的计算程序已发布在本团队开发的“地震变形计算平台”（http://124.207.150.85:8085/和http://10.16.24.41:8085/），成为计算震后变形的核心方法，适用于解释大地测量观测数据、理解震后变形机理、反演地球内部黏滞结构等科学问题研究.

2 球形地球模型地震位错理论的应用

上述球形地球模型地震位错理论为解决地震变形的相关问题提供了可靠的理论工具.通过球形地球模型地震位错理论，既可以根据震源参数计算地震变形，也可以根据地表变形反演震源参数和地球介质结构.我们提出的球形地球模型地震位错理论考虑了地球的分层、曲率效应、黏弹性结构、三维非均匀等，可以在全球地球动力学变化、地球内部黏弹性结构、断层反演等方面得到广泛应用，为深化地震周期过程的认识和地球动力学变化研究提供了有力支持.特别是，由于考虑了地球的球形形状，新地震位错理论具有研究地震全球动力学变化的独特潜力，而传统半无限空间理论无法做到.另一方面，现代大地测量观测技术提供的全球性的地震周期变形数据，也需要在球面上进行精细的建模和解释，以便保持模型的一致性.上述球形地球模型地震位错理论正好可以满足现代大地测量学和地震学发展的需求，为地震周期性循环过程的研究提供理论基础.

本团队在发展理论本身的同时，也就实际的地震变形问题，特别是特大地震事件，进行了针对性的理论扩展和应用研究.例如，针对地震造成的地球质心、体积、自转和势能变化等相关动力学变化问题，我们给出相应的计算方法，并根据实际的断层模型给出其定量结果.此外，还针对GNSS和GRACE这两种非常重要的现代大地测量观测技术，发展了利用远、近场GNSS数据约束断层反演的新方法，提出了利用GRACE震后重力变化反演区域黏弹性结构的新途径.研究结果表明：远场位移可有效地约束断层反演，综合震后位移以及重力变化可以确定黏弹性参数等.总之，地震位错理论应用研究可以大体归纳为以下几点.

2.1 地震引起的地球动力学变化

地震造成的地球内部质量重新分布不仅会引起地表位移和重力等物理场变化，而且会产生全球动力学变化，包括地心移动、体积变化、自转变化和势能变化等（Anderson，1974；Dahlen，1977；Chao，Gross，1987；Soldati，Spada，1999；Niet al，2005）.这些动力学变化对于理解地球内部物质循环、各圈层相互作用及其时空演化具有重要意义.由于本团队的地震位错理论是基于球形地球模型而建立，因此该理论具备研究全球动力学变化相关问题的优势.

基于上述球形地球模型位错理论和相应的基本物理定律，本团队发展了同震地心、地球体积、自转和势能等动力学参数变化的计算方法（Zhouet al，2013，2014a，2015，2016；Sun，Dong，2014；Xu，Sun，2014；Xuet al，2014，2016；Xu，Chao，2015，2017，2019；Xu，2021）.我们定量分析了几次大地震对地球极移和日长变化的影响，以及对地球参考框架中的质心位移和坐标轴指向变化的影响（Zhouet al，2013，2014a，2015，2016；Xu，Chao，2015，2019）.首次将地震作为部分物理激发源，解释现代大地测量观测的地球膨胀现象，并联合全球板块运动模型给出了不同时间尺度下的地球膨胀速率，确定了地震在板块运动中的角色（Xuet al，2016）；修正了计算同震日长变化的传统公式，并指出地震造成的地球转动惯量变化的张量迹不为零、不能被省略，并且提出了相应的数值计算方法（Zhouet al，2014a）；估计了全球大地震产生的地球自转累积变化，发现地震对地球自转轴长期漂移具有显著影响，是不可忽略的，并基于此提出了地球动力学参数变化新模型（Xu，Chao，2019）.

由于传统的地震学方法无法有效地估计地震产生的势能变化，本团队还提出了同震地球旋转动能和重力势能的计算方法.该方法不但适用于剪切型地震，而且适用于拉张型地震，为有效地研究同震重力势能变化与构造对应关系提供了理论基础（周江存等，2017；Xu，Chao，2017；Xu，2021）.据此，我们还分析了不同类型地震所产生的地球重力势能变化特征，结果表明地震产生的重力势能变化在地球内部的各个深度均有分布（周江存等，2017；Xu，Chao，2017）.Xu和Chao （2017）估计了全球历史地震产生的净重力势能以及对各圈层的影响，发现同震重力势能能流密度变化率约为全球平均大地热流的15%；并提出地震释放的重力势能可能作为区域热异常物理成因的假设，对地幔热机和地球长期演化研究具有启示作用.

2.2 断层滑动和介质结构的反演

地震位错理论不仅可以用于解释全球地震变形和地球动力学变化，还可以用于反演断层滑动、约束震源参数和地下介质结构（Yamasaki，Houseman，2012；Zhouet al，2014b）.利用大地测量和地球物理学观测数据对断层滑动分布、地下介质结构进行反演研究已逐渐成为一项常规性工作（Yabuki，Matsu’ura，1992；Reilingeret al，2000；Beresnev，2003；Ide，2007；Sun，Hartzell，2014；Diaoet al，2019）.然而，前人的反演通常是基于均匀或层状半无限空间模型计算的格林函数来进行.由于地球模型的过分简化必然会把地球模型的误差通过格林函数传递到反演结果中，从而影响对震源或地下结构的正确约束（Masterlark，2003；Williams，Wallace，2015；Gómezet al，2017；Marchandonet al，2021），因此，本团队使用接近于真实地球的球形地球模型来计算地震变形格林函数，并以数学上更为严密的和客观的马尔科夫链-蒙特卡洛（Markov chain Monte Carlo）算法来确定正则化因子，提出了一套基于球形地球位错理论的断层滑动模型反演方法（Zhouet al，2014b，2018；周新，2017）.利用上述反演方法和应用程序，我们研究了2004年苏门答腊MW9.3和2011年日本东北MW9.0等地震的同震-震后位移、空间/地面重力变化、断层滑动分布和震后慢滑动等地震变形问题（Zhouet al，2014b，2018；刘泰等，2017，2019）以及青藏高原动力学变化问题（周新等，2018）.

地震位错勒夫数和格林函数描述了地球介质对地震断层滑动的物理响应，反映了地下三维介质的结构信息，由大地测量观测数据反演的地震位错勒夫数和格林函数，可以体现研究区域的三维构造特征.本团队还提出了利用地震变形观测数据反演位错勒夫数和格林函数的方法，为约束区域三维结构提供了新思路.Sun等（2006b）提出了利用GRACE观测数据反演地震引力位变化勒夫数的基本思想；杨君妍等进一步发展完善了该反演方法，并以2011年日本东北MW9.0地震为例，利用GRACE和GNSS观测数据，实际反演了体现局部构造特征的地震位错勒夫数和格林函数，证实了该理论方法的可行性（Yanget al，2015；Yang，Sun，2020；杨君妍，孙文科，2020）.该方法可以用来计算地震多发地区的地震位错勒夫数和格林函数，用于准确计算该区域的理论地震变形或反演地球局部构造.

应用地震位错理论还可以反演地下的黏弹性参数.黏滞结构是地球最基本物理参数之一，是理解板块运动、俯冲带应力状态、冰后均衡调整（glacial isostatic adjustment，缩写为GIA）等科学问题的关键因素，然而如何精确地确定黏滞参数仍然是地球物理学特别是地球动力学亟需解决的科学问题（Cathles，1975；Weertman，Weertman，1975；Forteet al，1991；Fang，Hager，1996；Cornelio，Violay，2020；Arguset al，2021）.为此，王武星和张国庆等利用GRACE观测数据以及上述地震位错理论，反演了2004年苏门答腊MW9.3大地震震源区的黏滞参数，揭示了发震断层上盘与下盘黏滞参数的横向差异（Wanget al，2011；张国庆等，2015）.针对2011年日本MW9.0地震产生的震后变形，刘泰等（2017）提出了基于不同时间尺度的震后变形来分离断层余滑和黏滞变形的方法，其结果表明：可以以震后短期观测数据来反演断层余滑分布，以长期观测数据约束黏滞参数；该地震断层余滑两年后基本消失，四年后地幔黏滞性松弛效应超过断层余滑.该成果加深了对震后余滑时空分布和地幔黏滞性结构横向特征的认识.此外，梁明等（2018）还提出基于GNSS观测数据初步估算区域整体黏滞性结构，再利用GRACE观测数据修正深部黏滞系数，然后综合利用这两种观测数据调整黏滞系数，最终确定研究区域黏滞性结构的新方法.另外，陈飞等（2020）的研究结果也显示远场震后变形观测数据对地幔深部黏滞结构比较敏感，因此实际震后变形模拟时需考虑一个随深度变化的黏滞性模型.

2.3 空间观测技术的地震变形应用

重力卫星GRACE是现代大地测量主要技术之一，在地球科学中已经得到了广泛应用（Adam，2002；Tapleyet al，2019；张岚，孙文科，2022）.然而，GRACE能否检测到地震产生的重力变化，以及能否用来研究重力垂线偏差等是科研人员感兴趣但尚未解答的科学问题.另一方面，大地震通常发生在海陆俯冲带上，海域一般缺乏有效的地表变形观测，这使得断层滑动分布的反演变得非常困难.虽然GRACE具有观测包括海域在内的全球地震重力变化的潜力，但是由于缺少相关的理论支持，GRACE观测数据的应用颇为受限.

为此，Sun和Okubo （2004）最早对GRACE的地震变形检测能力进行了估计和预测.通过对地震位错勒夫数的球谐谱分析，他们发现GRACE可以检测出M＞9.0剪切型地震或M＞7.5张裂型地震的同震重力变化，该结论被随后发生的2004年苏门答腊地震所证实（Hanet al，2006）.本项研究成果首次阐明了GRACE对大地震的观测能力和范围，为重力卫星项目的科学规划提供了理论依据.之后，徐长仪等利用该理论和GRACE数据，成功地观测到2013年鄂霍次克海超深源地震MW8.3产生的同震重力变化，为进一步研究超深源地震的震源机制和变形特征提供了大地测量学途径（Xuet al，2017）.

我们发展了计算同震垂线偏差的新方法，为利用GRACE观测数据研究海域大地震变形问题提供了新途径.Sun和Zhou （2012）在上述球形地球模型地震位错理论的基础上，进一步发展了计算同震垂线偏差的新方法，使得利用GRACE观测数据研究海域地震变形问题极为方便.2011年日本东北大地震的研究结果表明，GRACE完全可以检测到大地震产生的垂线偏差变化（Sun，Zhou，2012），并且该变化对断层滑动模型更为敏感，为利用卫星重力观测数据研究地震断层滑动分布提供了新方法.本研究拓展了GRACE卫星数据的应用范围.

应用本团队的地震位错理论可以合理地解释全球地震变形观测数据如GNSS等.GNSS也是现代大地测量的主要技术之一，可以检测出接近全球范围内的同震和震后变形（Abidinet al，2009；Vernant，2015；Gautamet al，2020）.Fu和Sun （2006）基于上述地震位错理论，发现2004年苏门答腊-安达曼地震在6 000 km以外可以产生1 mm的同震水平位移，并可以被GPS观测到.通过理论计算和实际震例研究，王武星等发现中国大陆远场GPS同震位移场对环太平洋俯冲型地震震源机制具有约束能力（Wanget al，2014）.Wang等（2014）利用中国地壳运动观测网的GPS资料提取了2011年日本MW9.0地震引起的远场同震位移，并与理论同震位移进行了比较，其结果表明，近场位移可以更好地约束断层的几何形态，而远场GPS位移场可以约束断层滑动模型的地震矩.该研究结果扩展了远场GPS观测数据在地震变形研究中的应用潜力（Wanget al，2014）.此外，周新等（2018）以主喜马拉雅逆冲断层为例，模拟了同震和震间过程产生的大地水准面、重力和重力梯度变化，其结果显示地表的同震重力变化呈南北两极分布，足以被陆地重力测量所检测到.

3 讨论与结论

本文提纲挈领地介绍了本团队在球形地球模型的地震位错理论及其应用方面的科学贡献.该理论体系同时考虑了地球介质的自重效应、可压缩性、黏弹性、径向分层结构等因素，可以准确地模拟地震产生的位移、重力、应变、大地水准面、库仑应力变化、地球极移、势能变化等物理量，可归纳为两个方面：① 基于分层弹性球形地球、三维不均匀弹性地球模型、分层黏弹性地球模型的位错理论；② 球形地球模型位错理论在研究地震引起的地球动力学变化、反演断层滑动和地下介质结构、解释空间大地测量技术等方面的科学应用.

该理论体系建立了地震震源和变形场之间的理论联系，为研究震间变形、同震破裂和震后变形等提供了理论工具.同时，通过该理论的实际应用，解释了多个大地震的同震-震后变形数据，反演了断层滑动模型和地下黏弹性参数，深化了对地震循环过程的认识，凸显了球形位错理论在地震学和大地测量学研究中的理论支撑作用和广泛应用前景.因此，该理论对于推进地震基础理论进步、提升防震减灾科技支撑能力、推动大地震预测预报、服务我国防震减灾工作等具有重要意义.

限于本文讨论的范围，对于分层半无限空间模型的位错理论（如，Wanget al，2006）以及有限元类地震变形数值计算理论和方法（如，林晓光，孙文科，2014；Chenget al，2019）均未进行更多介绍，感兴趣的读者请参考其它文献.就数值类计算方法而言，它们可以针对更复杂的三维模型给出同震和震后变形的计算结果，但计算过程一般需要大型机器，可以预先计算一些三维模型的地震变形格林函数库，为一般用户进行地震变形观测数据解释以及断层和地下结构反演提供基础数据.

位错理论发展以地球几何模型的逐步完善为主线，目前以一维规则分层模型和基于微扰方法的三维模型为主.本文所介绍的三维不均匀弹性地球模型的位错理论也是基于微扰方法建立，尚不能考虑大幅度的横向非均匀性效应.我们也正在思考位错理论中包含地球的椭球效应、地形效应、内部介面起伏，甚至是复杂的横向非均匀性效应.在一维地震变形理论基础上，借鉴和拓展固体潮理论中的相关方法，应当可以直接考虑椭球效应、地形和内部界面的起伏，并求解地震变形微分方程组、建立相关理论和方法.而针对复杂的横向非均匀性效应，需要从数学上寻求其它新的方法，突破微扰方法的限制，才能基于接近于真实三维地球模型建立相关理论.此外，针对三维复杂地球模型，如何有效结合解析计算和数值计算的优势，在数值计算中充分利用部分解析计算，降低三维模型的计算时长，提出一种“混合模式”的计算理论和方法，也是一个值得探讨的问题.总之，笔者认为基于更精密的三维地球模型的位错理论和相关应用是本理论的前进方向，准静态位错理论将继续推进地震周期过程、地球内部介质结构和地球内部动力学过程等重要科学问题的深入研究.

感谢审稿专家提出建设性的修改意见.