Characteristics of strain release before moderate-strong earthquakes in Yunnan region
-
摘要: 利用云南区域台网地震目录资料，以地震应变能作为地震释放量度，基于应变释放加速模型对前兆地震序列进行非线性模拟. 结果表明: 云南地区1966—2013年发生的115例MS≥5.0地震中，具有加速特征的有79例，占69%；具有减速特征的有28例，占24%；数据太少无法判断加减速特征的有8例，占7%. 主震为孤立型的前兆地震序列，其减速或数据少类型所占比例远大于加速型；主震为走滑型的前兆序列，其加速型序列要比主震为正断型或逆冲型序列多15%左右. 搜索半径与前兆序列持续时间呈反比关系. 对前兆地震序列模拟结果表明，主震的预测震级和时间误差为±0.3和0.5年.Abstract: Referring to the earthquake catalogues produced by the Yunnan Regional Seismic Network，and setting the seismic strain energy as a scale of the seismic energy release，this paper conducts a nonlinear simulation for the precursory earthquake-sequences based on the accelerating strain release model. In the period from 1966 to 2013，115 earthquakes with MS≥5.0 occurred in Yunnan region. Among them，79 events are characterized by accelerating strain release，accounting for 69%. And 28 events are featured as decelerating strain release，accounting for 24%. The strain releases of the rest eight events accounting for 7% fail to show the features of acceleration or deceleration due to the less seismic data. As to the precursory earthquake sequences whose main shocks belong to the single-shock type，the ones whose earthquake strain releases are accelerating，or the ones lacking seismic data，are far more than the ones whose strain releases are accelerating. As to the precursory earthquake sequences whose main shocks belong to the strike-slip type，the ones whose earthquake strain releases are accelerating are 15% more than the ones whose main shocks belong to the normal-fault type，or 15% more than the thrust-fault type ones. The searching radius has an inverse relation with the duration of the precursory earthquake sequence. The simulated results for the precursory earthquake-sequences indicate that the error of the main shock’s magnitude is ±0.30，and the error of the main shock’s occurrence time is ±0.5 year.
-
引言
自1964年电离层异常扰动在美国阿拉斯加大地震震前首次被发现（Leonard，Barnes，1965），地震-电离层耦合效应逐渐得到各国研究人员的广泛关注。对于地震电离层耦合机理的大量研究，异常电场理论能较好地解释震前电离层异常扰动现象。该理论认为，地震孕育过程中由地表岩石破裂引起的一系列复杂的物理化学反应，产生了大气电场异常，进而叠加作用于电离层，引起了电离层电子浓度、电子温度等参量的显著变化（Pulinets et al，2000；Sorokin et al，2001；Kuo et al，2011）。近年来，以全球导航卫星系统（Global Navigation Satellite System，简写为GNSS）连续观测和电离层垂直观测为代表的地基电离层探测技术已被广泛应用于地震前兆监测和研究中（张学民等，2016；吴健等，2016；何宇飞等，2020）。基于该技术，研究人员可获取电离层垂直总电子含量（vertical total electron content，简写为VTEC）、电离层F2层临界频率（foF2）等电离层形态特征的重要指标。其中，电离层VTEC可通过GNSS台站观测数据反演获取，foF2可通过电离层垂测站观测获取。大量研究结果表明，MS5.0以上强震发生前数天到数小时内，上述参量可能会出现不同程度的异常变化（张学民等，2014；王壮凯等，2020），而多参量同步分析可为地震-电离层异常扰动判识提供佐证，对于强震频发地区的震情研判具有重要意义。
川滇地区位于青藏高原东南缘，包括川滇地块、巴颜喀拉地块南部、华南地块西部等二级地块。长期以来，在印度板块的推挤及华南地块的阻挡作用下，川滇地区的构造变形活动强烈，发育了规模、产状和活动速率各异的断裂，致使强震频发。由于地震短临预报的困难性，地震电离层监测预报模式逐渐成为川滇地区地震监测预报的新方向。目前，地震-电离层现象的研究主要处于震例研究和统计研究阶段（闫相相等，2014），区域电离层实时监测系统开发实例较少，电离层多参量异常关联及时空演化分析不够深入，快速获取电离层多参量异常扰动信息及其临震监测预报的应用还有待加强。
鉴于此，本文设计开发川滇地区的电离层多参量异常监测系统，旨在实现各类观测数据的自动化处理、电离层多参量图形及数据结果的自动化产出以及川滇地区电离层多参量异常变化的实时监测，为该区域地震电离层异常扰动判识提供证据。同时，本文拟通过2019年6月17日四川长宁MS6.0地震分析测试该系统图形产出功能的可靠性，以进一步验证川滇地区地震电离层监测预报的应用效能。
1. 系统概述
为了实现川滇地区电离层异常扰动变化的实时监测，本文设计开发川滇地区电离层多参量异常监测系统，快速产出全球和中国区域VTEC异常分布、空间参数、站点VTEC和foF2时序曲线，同时实现网页方式下图形及数据结果的查询、浏览、对比和下载等功能。基于该系统可开展示范性地震电离层监测工作，分析地震发生前后孕震区附近电离层多参量的异常扰动变化，尝试探索地震电离层现象的内在联系，为该地区地震监测预报提供电离层前兆异常的判定依据。
1.1 数据处理
1） 全球和中国区域VTEC异常分布。在中国区域按照均匀分布原则选取陆态网络的70个GNSS基准站。基于该基准站的观测数据，利用球谐函数建立大规模高精度格网电离层模型（林剑等，2009；王泽民等，2016），利用Bernese软件解算处理（Dach et al，2015）获取中国区域电离层VTEC格网数据（经纬度分辨率为1°×1°，时间分辨率为2 h），并结合欧洲轨道确定中心（Center for Orbit Determination in Europe，缩写为CODE）提供的全球电离层格网数据（经纬度分辨率为5°×2.5°，时间分辨率为1 h）自动监测全球和中国区域VTEC异常分布特征变化。
2） 空间参数、站点VTEC和foF2时序曲线。采用川滇地区GNSS基准站观测数据，通过单站点VTEC解算（熊波，2012），获取站点上空的VTEC数据（时间分辨率为15 min）。结合川滇地区垂测站（乐山、道孚、普洱、腾冲）的foF2观测数据（时间分辨率为15 min）、空间环境参数（F10.7，Kp，Dst，AE），自动监测站点VTEC、站点foF2以及其它空间参数的时序曲线特征变化。川滇地区GNSS、垂测台站的分布如图1所示。
![]()
图 1 川滇地区GNSS和垂测台站分布图Figure 1. Distribution of GNSS and vertical ionosonde stations in Sichuan-Yunnan region1.2 异常检验
选用滑动四分位距法提取电离层VTEC和foF2的异常信息，基于前15天相同时刻的观测值求取相应四分位距（IQR）及中位数（M），确定异常检验倍数（N），建立上（下）边界阈值M±N×IQR，观测值超出上（下）边界即视为异常，高于上边界为正异常，低于下边界为负异常。然后，从二维空间分布获取全球和中国区域VTEC异常分布，所用公式如下：
$$ \Delta {\rm{TEC}}\text{\text{＝}}\left\{\begin{array}{l}K{\text{－}}{L}_{1}\qquad K {\text{＞}} {L}_{1}{\text{，}}\\ K{\text{－}}{L}_{2}\qquad K {\text{＜}} {L}_{2}{\text{，}}\\ 0\;\qquad\qquad{L}_{2}{\text{≤}} K{\text{≤}} {L}_{1}{\text{，}}\end{array} \right.$$ (1) 式中，ΔTEC为VTEC异常扰动变化量，L1和L2分别为VTEC背景上、下边界阈值，K为VTEC观测值。站点VTEC和foF2时序曲线的异常变化获取方法同上。
2. 系统设计
2.1 结构设计
为满足川滇地区电离层多参量异常扰动监测的总体需求，结合数据获取、数据解算、图形及数据结果展示等需求，采用CS＋BS架构将该系统分为电离层VTEC异常分布产出子系统（简写为“异常分布子系统”）、GNSS站点VTEC及垂测站foF2曲线产出子系统（简写为“时序曲线子系统”）和电离层多参量异常监测平台（简写为“监测平台”）等三部分。其中：异常分布子系统用于监测全球和中国区域的电离层VTEC异常变化，实时产出全球和中国区域VTEC异常分布；时序曲线子系统用于监测川滇地区GNSS单站点VTEC及垂测站foF2的异常时序变化，实时产出空间参数、站点VTEC和foF2时序曲线；监测平台将各类图形及数据结果以网页方式自动发布，方便用户按照权限查询、浏览、对比和下载。两个子系统布设在不同服务器上，独立运行，方便维护管理。系统总体结构如图2所示。
![]()
图 2 川滇地区电离层多参量异常监测系统的总体结构图Figure 2. Overall structure diagram of the ionospheric multi-parameter anomaly monitoring system in Sichuan-Yunnan region2.2 功能设计
基于系统结构设计、业务需求及日常监测工作的要求，系统功能模块划分如图3所示。可见：异常分布子系统分为参数设置、自动产出、手动产出、数据服务、任务列表等五部分，主要实现全球和中国区域VTEC异常分布图形结果产出，同时提供数据转换、数据导出、数据下载等功能；时序曲线子系统分为参数设置、自动产出、手动产出、数据服务、任务列表等五部分，主要实现站点VTEC和foF2、空间参数时序曲线图形结果的产出，同时提供数据导出、下载等功能。监测平台分为系统管理、图形展示、图形对比、数据下载等四部分，主要实现VTEC异常分布，空间参数、站点VTEC和foF2时序曲线查询、浏览，同时提供电离层多参量数据分析结果的下载功能。
![]()
图 3 川滇地区电离层多参量异常监 测系统功能图Figure 3. The function diagram of the ionospheric multi-parameter anomaly monitoring system in Sichuan-Yunnan region2.3 业务流程
异常分布子系统。该子系统可分别通过自动和手动模式以FTP方式下载各类数据，自动调用Bernese数据处理软件，基于球谐函数建立中国区域电离层格网模型，解算中国区域电离层VTEC格网数据。结合CODE全球电离层格网数据和空间环境参数（F10.7，Kp，Dst，AE），提取全球和中国区域VTEC异常信息，利用ArcGIS平台分别产出全球和中国区域VTEC异常分布以及格网点VTEC时序曲线，自动产出结果延迟1—2天，具体业务流程如图4所示。
时序曲线子系统。该子系统通过自动和手动模式以FTP方式下载各类数据，自动调用GPS-TEC数据处理软件，获取单站点VTEC数据。结合川滇地区垂测站（乐山、道孚、普洱、腾冲） foF2观测数据、空间环境参数（F10.7，Kp，Dst，AE），提取站点VTEC和foF2异常信息，并利用Matlab平台分别产出站点VTEC和foF2、空间参数的时序曲线，自动产出结果延迟1—2天，具体业务流程如图5所示。
监测平台。该平台实时访问异常分布子系统和时序曲线子系统，将各类图形、数据产出结果通过http自动发布，根据用户权限为用户提供查询、浏览、对比、下载等差异化功能。
3. 系统实现和震例分析
3.1 系统实现
该系统性能要求较高，需要实现各类数据的自动化下载、多个第三方软件的自动化运行和监测、各类图形及数据结果的自动化生成与发布。为实现系统运行稳定性、人机交互性及功能扩展性等多方面要求，具体实现内容如下：
1） 异常分布子系统。基于Windows 2008/2012操作系统，采用C/S （Client-Server，即客户端/服务器模式）架构，将面向对象的编程语言C#作为后台开发语言，MySQL作为后台数据库管理系统，ArcGIS作为图形结果产出平台，通过分层技术体系结构实现系统运行的稳定性和交互性。子系统主界面分为功能菜单、用户操作及图形展示等三部分，其中：功能菜单位于主界面左上方，包含系统配置、手动产出、时序曲线、数据转换、数据导出、数据下载、任务列表、站点管理等功能；用户操作位于功能菜单左下方，包含查询条件和查询结果两部分；图形展示位于用户操作界面的右侧，根据查询条件和查询结果查看电离层VTEC异常分布情况，全球VTEC异常分布间隔一小时，中国区域VTEC异常分布间隔两小时（图6a）。
![]()
图 6 系统主界面（a） 异常分布子系统；（b） 时序曲线子系统；（c） 监测平台Figure 6. Main system interfaces（a） Anomaly distribution subsystem；（b） Time series curve subsystem；（c） Monitoring platform2） 时序曲线子系统。基于Windows 2008/2012操作系统，采用C/S架构，将C#作为后台开发语言，MySQL作为后台数据库管理系统，Matlab作为图形结果产出平台，通过分层技术体系结构实现系统运行的稳定性和交互性。子系统主界面分为功能菜单、用户操作及图形展示等三部分，其中：功能菜单位于主界面左上方，包含手动产出、数据导出、数据下载、任务列表、站点设置、系统配置等功能；用户操作位于功能菜单左下方，包含查询条件和查询结果两部分；图形展示位于用户操作界面的右侧，根据查询条件和查询结果分别展示站点VTEC和foF2、空间参数时序曲线（图6b）。
3） 监测平台。基于Windows操作系统，采用B/S （Browser-Server，即浏览器/服务器模式）架构，将C#作为后台开发语言。该平台设置用户权限，通过网页方式登录。功能菜单位于主界面左侧，包含系统管理、VTEC异常分布、空间参数曲线、站点VTEC曲线、站点foF2曲线、图形对比和数据下载等功能。选择功能菜单，打开用户操作界面，根据查询条件和查询结果分别展示VTEC异常分布及空间参数、站点VTEC和foF2的时序曲线。同时，平台提供不同时期的电离层多参量图形结果对比展示功能，方便用户针对电离层前兆异常分析判断（图6c）。
3.2 震例分析
为测试系统图形产出功能的可靠性，验证川滇地区地震电离层监测预报的可行性，本文选择2019年6月17日四川长宁MS6.0地震前后的观测数据，利用该系统手动产出全球和中国区域VTEC异常分布及站点VTEC、foF2和空间参数的时序曲线，综合分析长宁地震震前电离层多参量异常扰动变化特征。鉴于篇幅有限，本文仅展示个别图形结果（图7，8）。
![]()
图 7 异常分布子系统关于长宁地震前2019年6月14日12:00 UT全球（a）和中国区域（b）VTEC异常分布图形结果，图中红色圆点代表震中Figure 7. Global （a） and China regional （b） VTEC anomaly distribution graphic results of the anomaly distribution subsystem at 12:00 UT on June 14，2019 before Changning earthquake，where red dots represent the epicenter从全球和中国区域VTEC异常分布（图7）可以看出，6月14日12:00 UT （北京时间20:00），震中附近区域均出现显著的电离层VTEC正异常，地震电离层耦合的局部性、集中性和显著性特征明显。同时，以个别站点为例，站点VTEC （图8a）和站点foF2 （图8b）的时序曲线显示6月14日四川宁南GNSS站点VTEC出现正异常，云南腾冲垂测站foF2同步出现正异常，该同步异常现象与VTEC异常分布结果基本吻合。结合空间参数时序曲线变化（图8c）可知，6月14日空间环境参数F10.7为68 sfu，Kp最大值为2.3，Dst最小值为−12 nT，AE最大值为338 nT，均未超过参数阈值，可以排除太阳活动、地磁变化对电离层异常扰动的影响。同时，长宁地震发生前后该异常区域并未发生其它强震事件。综上，2019年6月14日震中附近区域电离层多参量异常扰动满足地震电离层异常判识标准，可能是地震前兆异常，与长宁MS6.0地震具有一定的相关性。
![]()
图 8 时序曲线子系统图形结果（2019年6月1—20日）（a） 四川宁南站点VTEC时序曲线；（b） 云南腾冲站点foF2时序曲线；（c） 空间参数时序曲线，图中红色直线为参数阈值Figure 8. Graphic results of time series curve subsystem （June 1 to June 20，2019）（a） VTEC time series curves of Ningnan station in Sichuan Province；（b） foF2 time series curves of Tengchong station in Yunnan Province；（c） Time series curves of spatial parameter，where all red lines are parameter thresholds4. 讨论与结论
本文主要介绍了川滇地区电离层多参量异常监测系统的结构设计、功能设计和业务流程设计，并给出了系统实现结果。该系统实现了全球和中国区域VTEC、站点VTEC和foF2异常变化的实时监测，有助于开展示范性地震-电离层异常监测工作，目前该系统已成功应用于四川省地震局的监测预报日常工作，图形产出结果已参与地震监测预报会商。
该系统虽然较大地提高了川滇地区电离层异常变化监测能力，但对电离层参量的监测手段较为单一，电离层异常对比分析功能较弱。该系统下一步仍需不断补充完善，基于地基和空基电离层探测技术开展天地一体化联合观测，并结合临震电离层多参量同步分析为各参量异常指标提供量化信息，为地震监测任务提供更便捷的技术支撑，为深入研究地震孕育期、临震、主震及余震中的电离层异常扰动机理提供便捷的数据支撑。
CODE提供了全球电离层格网数据，中国地震台网中心和中国地震局地震预测研究所提供了GNSS台站和电离层垂测台站的观测数据，作者在此一并表示衷心的感谢。
-
![]()
图 2 1965年1月—2013年12月云南地区的地震分布
Red circles represent M≥2.5 earthquakes，and blue circles represent M≥5.0 earthquakes
Figure 2. Distribution of earthquakes in Yunnan region from January of 1965 to December of 2013
![]()
图 3 1970年1月—1979年12月滇西和滇西南地区G-R关系
Figure 3. G-R relation in western and southwest Yunnan from January of 1970 to December of 1979
![]()
图 4 1982年2月20日江城MS5.0地震前时间搜索得到的应变释放曲线及其参数实例
(a) 最佳应变释放曲线族; (b) 搜索半径R与起始时间的关系
Figure 4. Strain-release curves obtained by time searching before the Jiangcheng MS5.0 earthquake on February 20，1982
(a) Family of best strain-release curves；(b) Relationship between searching radius and starting time for different range of the index m
![]()
图 5 4次应变释放加速曲线震例
(a) 1982年2月20日江城MS5.0; (b) 1984年4月24日孟连MS6.3; (c) 1998年11月19日宁蒗MS6.2; (d) 2011年3月10日盈江MS5.8
Figure 5. Four examples of accelerating strain release curves
(a) Family of best strain-release curves；(b) Relationship between searching radius and starting time for different range of the index m
![]()
图 6 归一化标准差S等值线影像图 图中标准差S已归一化到1，tf为预测时间，单位为年(十进制)；Mf为预测震级．Ⅰ为直接拟合的标准差；Ⅱ为改善后的标准差 (a) 1984年4月24日孟连MS6.3地震; (b) 1991年4月12日宁蒗MS5.1地震；(c) 1994年12月30日马边MS6.0地震；(d) 2011年3月10日盈江MS5.8地震
Figure 6. Contour image maps of the normalized standard deviation S Standard deviation S has been normalized to 1，tf represents predicted time in years (decimal)，Mf represents predicted magnitude. The left panels (Ⅰ) give the direct fitting standard deviation，and the right panels ( Ⅱ ) give the improved standard deviation. (a) Menglian MS6.3 earthquake on 24 April 1984; (b) Ninglang MS5.1 earthquake on 12 April 1991；(c) Mabian MS6.0 earthquake on 30 December 1994；(d) Yingjiang MS5.8 earthquake on 10 March 2011
![]()
图 7 预测时间和预测震级的误差统计 (a) 预测时间误差Δtf；(b) Δtf 统计百分比; (c) 预测震级误差ΔMf; (d) ΔMf统计百分比
Figure 7. Error statistics of predicted time and magnitude (a) Predicted time error Δtf; (b) The percentage of Δtf; (c) Predicted magnitude error ΔMf; (d) The percentage of ΔMf
![]()
图 8 加减速序列搜索半径R与m关系 (a) 应变释放加速(1982年2月20日江城MS5.0地震); (b) 应变释放减速(1999年11月25日澄江MS5.2地震)；(c) 数据少及干扰(1968年6月13日麻栗坡MS5.2地震)
Figure 8. Examples of relationship between searching radius R and m value in the acceleration or deceleration sequences (a) Accelerating strain release process taking the Jiangcheng MS5.0 earthquake as an example; (b) Decelerating strain release process taking the Chengjiang MS5.2 earthquake as an example; (c) Less seismic data and interference taking the Malipo MS5.2 earthquake as an example
表 1 主震事件参数和模拟结果
Table 1 Parameters of main shock events and simulation results
编号 发震时间 地名 主震 震级 起始 震级 最优 m值 R/km Δtf/a ΔMf 加速 比例 加减速特征 年-月-日 1 1966-01-31 中甸 5.1 3.0 0.18 120 0.23 0.50 0.68 + 2 1966-02-05 东川 6.5 3.5 0.10 300 0.28 0.39 0.25 * 3 1966-09-18 江城 5.6 3.0 0.08 280 0.13 0.67 0.54 + 4 1966-09-19 瑞丽 5.4 3.0 0.23 210 0.18 0.17 0.36 - 5 1966-09-23 宣威 5.0 3.0 0.20 180 0.30 0.61 0.29 - 6 1966-09-28 中甸 6.4 3.4 0.20 270 0.13 0.12 0.43 - 7 1966-10-11 永善 5.2 3.0 0.28 130 0.15 0.24 0.43 - 8 1967-11-04 沧源 5.0 3.0 0.26 250 0.22 0.57 0.25 * 9 1968-03-16 南华 5.1 3.0 0.27 50 0.45 0.22 0.43 - 10 1968-06-13 麻栗 5.2 3.0 0.57 210 0.36 0.03 0.18 * 11 1969-02-09 景洪 5.4 3.0 0.21 220 0.40 0.51 0.32 - 12 1970-01-05 通海 7.8 4.8 0.14 300 0.10 0.28 0.07 * 13 1970-02-07 普洱 6.2 3.2 0.28 130 0.28 0.07 0.46 - 14 1970-07-31 雷波 5.4 3.0 0.21 60 0.35 0.12 0.36 - 15 1971-02-05 保山 5.8 3.0 0.22 60 0.30 0.08 0.25 - 16 1972-01-23 石屏 5.6 3.0 0.10 70 0.37 0.14 0.71 + 17 1973-03-22 景洪 5.5 3.0 0.14 290 0.28 0.63 0.32 - 18 1973-04-22 彝良 5.1 3.0 0.26 70 0.36 0.13 0.54 + 19 1973-06-01 腾冲 5.0 3.0 0.48 50 0.36 0.14 0.86 + 20 1973-08-16 普洱 6.3 3.3 0.23 90 0.17 0.26 0.93 + 21 1974-05-11 大关 7.1 4.1 0.04 230 0.47 -0.33 0.50 + 22 1975-01-12 楚雄 5.4 3.0 0.54 60 0.45 -0.03 0.82 + 23 1975-07-09 开远 5.3 3.0 0.26 80 0.29 0.13 0.96 + 24 1975-09-04 漾濞 5.0 3.0 0.74 30 0.38 -0.17 0.54 + 25 1975-10-28 勐腊 5.7 3.0 0.25 170 0.12 0.18 0.93 + 26 1976-02-16 普洱 5.8 3.0 0.28 90 0.18 0.11 0.96 + 27 1976-05-29 龙陵 7.3 4.3 0.43 280 0.39 -0.33 0.75 + 28 1976-10-09 峨山 5.3 3.0 0.36 40 0.39 -0.01 1.00 + 29 1976-11-07 盐源 6.7 3.7 0.34 260 0.41 -0.22 0.93 + 30 1977-03-17 漾濞 5.4 3.0 0.20 100 0.20 0.27 0.96 + 31 1978-05-19 下关 5.1 3.0 0.27 60 0.25 0.03 0.71 + 32 1978-09-10 普洱 5.8 3.0 0.20 130 0.39 0.33 0.93 + 33 1979-03-15 普洱 6.8 3.8 0.27 160 0.30 0.09 0.79 + 34 1979-08-09 保山 5.0 3.0 0.18 30 0.27 0.05 0.25 - 35 1980-02-02 木里 5.8 3.0 0.35 40 0.39 -0.07 0.86 + 36 1980-06-18 蒙自 5.8 2.8 0.15 130 0.20 0.11 1.00 + 37 1981-05-22 宁蒗 5.0 2.5 0.60 190 0.18 0.36 1.00 + 38 1981-07-07 盈江 5.4 3.0 0.39 50 0.43 -0.06 0.89 + 39 1981-09-19 普洱 6.0 3.0 0.19 120 0.31 0.29 0.93 + 40 1982-02-20 江城 5.0 2.5 0.18 80 0.22 0.23 0.79 + 41 1982-07-03 剑川 5.4 2.5 0.54 40 0.51 -0.10 0.43 - 42 1982-10-08 洱源 5.1 2.5 0.22 180 0.10 0.52 0.82 + 43 1982-10-27 富宁 5.8 2.8 0.10 270 0.36 0.31 0.39 - 44 1982-12-28 景洪 5.3 2.5 0.27 40 0.34 0.12 1.00 + 45 1982-12-28 永平 5.4 2.5 0.35 50 0.31 0.00 0.86 + 46 1983-05-26 景洪 5.1 2.5 0.21 120 0.18 0.57 0.21 * 47 1983-06-04 鲁甸 5.0 2.5 0.30 250 0.12 0.55 0.93 + 48 1984-04-07 景洪 5.1 2.5 0.04 170 0.13 0.40 0.21 * 49 1984-04-24 孟连 6.3 3.3 0.15 80 0.17 0.03 0.93 + 50 1985-04-18 禄劝 6.2 3.2 0.24 110 0.32 -0.01 0.64 + 51 1985-09-02 建水 5.4 2.5 0.36 30 0.26 -0.05 0.36 - 52 1986-03-13 鹤庆 5.4 2.5 0.51 40 0.53 -0.02 0.11 - 53 1986-07-10 盐源 5.0 2.5 0.21 30 0.31 0.10 0.29 - 54 1986-10-07 富民 5.2 2.5 0.22 40 0.25 -0.02 0.68 + 55 1987-05-18 鹤庆 5.0 2.5 0.19 80 0.19 0.51 0.96 + 56 1987-09-26 施甸 5.0 2.5 0.21 270 0.10 0.64 0.82 + 57 1988-01-10 宁蒗 5.5 2.5 0.19 160 0.19 0.49 0.89 + 58 1988-04-15 会东 5.4 2.5 0.50 130 0.29 -0.01 0.75 + 59 1988-08-15 澜沧 5.0 2.5 0.21 280 0.11 0.62 0.54 + 60 1988-11-06 澜沧 7.4 4.4 0.15 130 0.38 0.08 0.46 - 61 1989-02-10 洱源 5.1 2.5 0.25 80 0.24 0.31 0.82 + 62 1989-09-20 嵩明 5.0 2.5 0.20 30 0.09 0.16 0.46 - 63 1991-04-12 宁蒗 5.1 2.5 0.36 40 0.25 0.13 0.93 + 64 1991-07-01 施甸 5.0 2.5 0.27 50 0.24 0.37 0.86 + 65 1992-12-18 永胜 5.4 2.5 0.29 50 0.44 0.12 0.54 + 66 1993-01-27 普洱 6.3 3.3 0.22 150 0.37 0.15 0.75 + 67 1993-02-01 大姚 5.3 2.5 0.19 30 0.41 -0.03 0.75 + 68 1993-07-17 中甸 5.6 2.6 0.29 60 0.26 0.02 0.86 + 69 1993-08-07 犍为 5.4 2.5 0.20 120 0.10 0.22 1.00 + 70 1993-08-14 姚安 5.5 2.5 0.12 150 0.09 0.41 0.96 + 71 1994-09-19 景谷 5.2 2.5 0.24 290 0.10 0.70 0.57 + 72 1994-12-30 马边 6.0 3.0 0.18 50 0.29 0.11 0.96 + 73 1995-02-18 澜沧 5.1 2.5 0.33 190 0.14 0.56 0.79 + 74 1995-05-12 临沧 5.1 2.5 0.20 240 0.05 0.51 0.68 + 75 1995-07-12 孟连 7.3 4.3 0.09 70 0.39 0.27 0.79 + 76 1995-10-24 武定 6.5 3.5 0.33 110 0.42 -0.15 0.82 + 77 1996-02-03 丽江 7.0 4.0 0.24 120 0.47 -0.23 0.82 + 78 1996-06-13 墨江 5.0 2.5 0.43 120 0.03 0.04 0.75 + 79 1997-01-25 景洪 5.1 2.5 0.49 120 0.08 0.15 0.71 + 80 1998-11-19 宁蒗 6.2 3.2 0.20 210 0.04 0.09 0.89 + 81 1998-12-01 宣威 5.1 2.5 0.21 40 0.42 0.00 0.82 + 82 1999-11-25 澄江 5.2 2.5 0.73 50 0.42 0.34 0.43 - 83 2000-01-15 姚安 6.5 3.5 0.22 80 0.24 -0.01 1.00 + 84 2000-01-27 丘北 5.5 2.5 0.30 60 0.21 -0.11 0.57 + 85 2000-08-21 武定 5.1 2.5 0.22 180 0.06 0.08 1.00 + 86 2001-03-12 澜沧 5.0 2.5 0.13 210 0.04 0.33 0.25 - 87 2001-04-12 施甸 5.9 2.9 0.14 300 0.02 0.25 0.75 + 88 2001-05-24 宁蒗 5.8 2.8 0.19 160 0.31 0.21 0.96 + 89 2001-07-10 楚雄 5.3 2.5 0.23 100 0.10 0.19 1.00 + 90 2001-07-15 江川 5.1 2.5 0.21 100 0.07 0.03 0.43 - 91 2001-09-04 景谷 5.0 2.5 0.30 280 0.04 0.21 0.50 + 92 2001-10-27 永胜 6.0 3.0 0.15 140 0.24 0.28 1.00 + 93 2003-07-21 大姚 6.2 3.2 0.21 40 0.40 -0.17 0.64 + 94 2003-08-21 盐源 5.0 2.5 0.28 140 0.06 0.20 0.61 + 95 2003-11-15 鲁甸 5.1 2.5 0.33 140 0.23 0.02 0.57 + 96 2004-08-10 鲁甸 5.6 2.6 0.22 140 0.26 0.35 0.93 + 97 2004-10-19 保山 5.0 2.5 0.20 30 0.30 0.34 0.68 + 98 2004-12-26 双柏 5.0 2.5 0.19 90 0.34 0.27 0.25 * 99 2005-01-26 思茅 5.0 2.5 0.23 40 0.33 0.09 0.43 - 100 2005-08-05 会泽 5.3 2.5 0.21 170 0.26 0.56 0.75 + 101 2005-08-13 文山 5.3 2.5 0.16 80 0.05 0.03 0.64 + 102 2006-01-12 墨江 5.0 2.5 0.20 200 0.06 0.09 0.71 + 103 2006-07-22 盐津 5.1 2.5 0.19 280 0.21 0.50 0.86 + 104 2007-06-03 宁洱 6.4 3.4 0.14 80 0.44 0.15 0.18 - 105 2008-03-21 盈江 5.0 2.5 0.47 220 0.09 0.30 0.46 - 106 2008-08-21 盈江 5.9 2.9 0.49 90 0.43 -0.06 0.64 + 107 2008-08-30 攀枝 6.1 3.1 0.30 60 0.48 -0.24 0.32 - 108 2009-07-09 姚安 6.0 3.0 0.31 240 0.38 0.25 0.79 + 109 2009-11-02 宾川 5.0 2.5 0.32 50 0.33 0.09 0.11 * 110 2010-02-25 元谋 5.1 2.5 0.16 230 0.13 0.75 0.68 + 111 2011-03-10 盈江 5.8 2.8 0.19 40 0.05 0.00 0.86 + 112 2012-06-24 宁蒗 5.7 2.7 0.19 50 0.39 0.08 0.14 - 113 2012-09-07 彝良 5.7 2.7 0.48 50 0.20 -0.27 0.18 - 114 2013-03-03 洱源 5.5 2.5 0.53 100 0.26 -0.22 0.46 - 115 2013-08-31 中甸 5.9 2.9 0.22 50 0.42 -0.16 0.79 + 注： R为搜索半径；Δtf为预测时间误差；ΔMf为预测震级误差；“+”代表序列有加速特征；“-”代表序列有减速特征；“*”代表序列数据少无法判断加减速特征. 
下载: 导出CSV
-
期刊类型引用(7)
1. 郑秋月，黄江培，吴宇琴，陈政宇，刘东，王青华. 2021年漾濞M_S6.4地震前重力段差指标量分析及场源特征反演. 地震研究. 2025(01): 22-31 . 
百度学术
2. 刘代芹，玄松柏，王晓强，陈丽，李瑞，张新林，艾力夏提·玉山，朱治国. 北天山地区近期重力场及地壳密度变化特征研究. 大地测量与地球动力学. 2021(05): 459-465 . 百度学术
3. 陈兆辉，张双喜，刘金钊，张品，朱传东. 地震重力网经典平差与拟稳平差综合分析与研究. 武汉大学学报(信息科学版). 2021(11): 1693-1702 . 百度学术
4. 李玉森，底欣欣，臧石磊，李伟，商占鑫. 辽宁省重力异常指标体系研究. 防灾减灾学报. 2019(01): 13-19 . 百度学术
5. 王同庆，张品，刘金钊，陈兆辉，张双喜，朱传东. 基于重力测点和段差变化的前郭5.8级震群研究与三维反演. 大地测量与地球动力学. 2019(05): 469-475 . 百度学术
6. 陈丽，艾力夏提·玉山，李杰，刘代芹，李桂荣，李瑞，孙小旭，张治广，方伟. 基于流动重力和连续重力的新疆区域重力场时空特征研究. 内陆地震. 2019(03): 236-248 . 百度学术
7. 陈兆辉，王同庆，刘金钊，朱传东，张品，张双喜. 南北地震带南段近期重力场变化与指标量分析. 地震研究. 2018(03): 480-486+490 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 562
- HTML全文浏览量: 224
- PDF下载量: 6
- 被引次数: 7
