漾濞MS6.4地震前后的重力变化特征及其孕震含义分析

doi:10.3969/j.issn.0253-4967.2022.06.012

摘要/Abstract

摘要：

2021年5月21日云南省漾濞县发生 M_S6.4 地震, 该地震打破了云南区域自2014年以来近7a的M_S>6.0地震平静。文中采用基于贝叶斯原理的平差方法解算了2015-2021年的相对重力流动观测数据, 获取了研究区内多年的时变重力值, 分析了漾濞 M_S6.4 地震前后的重力变化特征。研究结果显示, 本次地震地质活动的活跃期为2017-2018年, 2018-2020年地质活动相对稳定, 2020-2021年出现地下物质“逃逸”现象, 地震发生于“逃逸”期; 以震中为中心的200km、 100km、 50km、 25km半径范围内, 重力异常值呈四象限分布特征; 震前3a重力正、负变化更加显著, 达80×10^-8m/s², 震前2a重力变化较小, 但依然维持四象限分布特征; 震后研究区内的重力场均发生反向变化, 变化量级达震前3a的累计变化量; 本次地震发生于重力异常四象限中心附近, 且与红河断裂带北段重合, 故推断红河断裂为漾濞地震的孕震构造。文中的研究方法和成果可为时变重力场变化特征研究和震前重力信号分析与解释提供震例参考; 同时, 可为孕震区地质构造运动的孕震含义提供有效解释, 并对下一步可能发生中强震区域的探讨及预判提供参考。

关键词: 时变重力场, 漾濞地震, 贝叶斯平差, 时空域

Abstract:

On May 21, 2021, an earthquake with a magnitude of 6.4 occurred in Yangbi, Yunnan, China, which broke the calm of earthquakes above M_S6.0 in Yunnan Province for nearly seven years since 2015. This paper collects the relative gravity observation data of Yunnan from 2015 to 2021, and adopts the parameters estimation method based on Bayesian principle to obtain the station-value series of multi-year time-varying gravity field. Through years of time-varying gravity data, and according to the guiding idea of “Seeking the source from the field and combining the field and the source”, the series data are divided into different “time-space domains”. The “time domain” takes the epicenter as the center to analyze the one-year difference changes of the measuring points within a radius of 200km since 2015 and the multi-year cumulative change deduced backward from the origin time of earthquakes. The “space domain” takes the epicenter as the center, collects the point value change sequence of the measuring points with the radius of 200km, 100km, 50km and 25km respectively, and uses the method of “layered processing and step-by-step analysis” to analyze the characteristics of gravity change before the earthquake. The research results show that the geological active period of this earthquake is 2017-2018, the geological activity is relatively stable in 2018-2020 and there was the phenomenon of “escape” of underground materials in 2020 to 2021, and the earthquake occurred in the “escape” period. The abnormal gravity near the epicenter is characterized by a four-quadrant distribution, the positive and negative changes of gravity in the three years before the earthquake were more significant, reaching 80×10^-8m/s², the change of gravity in the two years before the earthquake was small, but the four-quadrant distribution characteristics were still maintained. After the earthquake, the gravity field in the study area changed inversely, and the magnitude of change reached the cumulative change of 3 years before the earthquake. This earthquake occurred near the center of the four quadrants of the gravity anomaly and coincided with the northern segment of the Honghe fault zone, so it is inferred that the Honghe Fault is the seismic structure of the Yangbi earthquake. The research methods and results of this paper can provide a reference for earthquake examples for the study of time-varying gravity field variation characteristics and the analysis and interpretation of pre-earthquake gravity signals, as well as for exploring and predicting the areas of future mid-strong earthquakes.
This paper systematically analyzes the characteristics of gravity changes in different “time-space domains” of the 2021-05-21 Yangbi M_S6.4 earthquake, and attempts to analyze its seismogenic meaning. The following conclusions can be drawn:
(1)The gravity anomaly appeared before the Yangbi earthquake and lasted for a long time. The phenomenon of energy accumulation appeared four years before the earthquake.
(2)The variation of gravity point value within the radius of 200km is basically consistent, and the variation of gravity point value within the radius of 100km is highly consistent, indicating that the range of gravity anomaly during the earthquake preparation period reached more than 200km, which is consistent with the research results of previous empirical statistics.
(3)The geological activity associated with this earthquake started at about 4 years before the earthquake, then there was a calm period of up to two years. Half a year before the earthquake, there was an overall decline in the value of gravity points near the epicenter, which is in line with the common sense in seismology that “seismic energy gathers first and then diffuses, and the earthquake occurs in the diffusion period”, but there is a calm period before this earthquake.
(4)The “four quadrant” phenomenon appeared and gradually increased in the variation trend diagram of gravity field in one year, two years and three years before the earthquake.
(5)The “four quadrants” of the variation trend of the gravity field before the earthquake are generally consistent with the focal mechanism solution and the variation of the coseismic displacement field issued by the USGS and Seismology research group, indicating that the phenomenon of “four quadrants” variation in the gravity field will appear before earthquake of strike-slip fault geological structure, and the gravity variation can reveal to a certain extent the earthquake-generating mechanism before earthquake.
(6)The change trend of gravity after the earthquake is completely opposite to that before the earthquake, which is the release and recovery of compression and accumulation of underground materials before the earthquake. Therefore, it is considered that this earthquake can effectively release the earthquake preparation energy in this area, and it is unlikely that earthquakes of M_S≥6.0 will occur again in a short time.

Key words: time-varying gravity field, Yangbi earthquake, Bayesian adjustment method, time-space domains

中图分类号:

P315.72⁺6

黄江培, 曹颖, 刘东, 郑秋月, 陈政宇, 王青华. 漾濞M_S6.4地震前后的重力变化特征及其孕震含义分析[J]. 地震地质, 2022, 44(6): 1557-1573.

HUANG Jiang-pei, CAO Ying, LIU Dong, ZHENG Qiu-yue, CHEN Zheng-yu, WANG Qing-hua. ANALYSIS OF GRAVITY VARIATION CHARACTERISTICS AND ITS SEIMOGENIC MEANING BEFORE AND AFTER YANGBI M_S6.4 EARTHQUAKE[J]. SEISMOLOGY AND GEOLOGY, 2022, 44(6): 1557-1573.

图/表 11

图 1 云南地震相对重力观测网

Fig. 1 Seismic relative gravity observation network in Yunnan Province.

表1 重力变化异常范围、异常量级与震级之间的关系统计

Table1 Statistics on the relationship between gravity variation anomaly range, anomaly magnitude and magnitude

震级/M_S	4	5	6	7	8	文献
异常区域	70	140	220	350	660	2019
长直径/km		≥100	≥200	≥400	≥600	2018
异常量级/10^-8m·s^-2	28	50	72	93	115	2019
		≥50	≥80	≥100	≥120	2018

表2 不同半径范围内的测点分布

Table2 Distribution of measuring points in different radius ranges

方位	半径/km
	200	100	50	25
东北	虎跳峡、丽江机场、程海镇、中台、宜就、赵家店A、黄金湾、睦科临、龙潭A、烟地A、丽江、永胜新、白汉场岔口、大厂A、飞崖、仁和镇、大姚井、桥头河、哲理新、荣将	拉渣坡、江昌门、大丽72km、新坪A、剑永2、剑川A、海稍、松桂新、片角新、三庄	农机所、邓川临	大理临
东南	丽白岔口、南华A、沙桥、田房箐新、凤庆、云县新、安召道班、云县区域站、景云桥、邑栖么、景东A、姚安、楚雄A、忙怀A、小定喜、南涧、宝华	南涧区域、清华洞、祥云A、海坝庄、果园新、苴力、弥渡	下关临、感通寺、实验场、下关、飞来寺	平坡
西北	泸水县、维西县、称杆乡、福贡县、通甸	云龙区域	长邑、炼铁	漾濞
西南	腾冲基准台、小桥、腾冲背景场、施甸背景场、扁口、西邑、链子桥、永德、龙陵区域站、新城区域站、梁河	保山基准、老营、昌宁、保山临、瓦窑	永平区域站、哨湾	黄连铺

图 2 1年期重力场差分变化图

Fig. 2 One-year differential change chart of gravity field.

图 3 多年期重力场累积变化图

Fig. 3 Multi-year cumulative change chart of gravity field.

图 4 200km范围内点值的时变序列图

Fig. 4 Time-varying sequence chart of gravity point values within 200km.

图 5 100km范围内点值的时变序列图

Fig. 5 Time-varying sequence chart of gravity point values within 100km.

图 6 50km范围内点值的时变序列图

Fig. 6 Time-varying sequence chart of gravity point values within 50km.

图 7 25km范围内点值的时变序列图

Fig. 7 Time-varying sequence chart of gravity point values within 25km.

图 8 漾濞 MS6.4 地震的震源机制解

Fig. 8 Focal mechanism solution of Yangbi MS6.4 earthquake.

图 9 震前震后(2021-03-2021-06)的变化趋势图

Fig. 9 Variation trend before and after the earthquake(2021-03-2021-06).

参考文献 28

[1]	陈运泰, 顾浩鼎, 卢造勋. 1980. 1975年海城地震与1976年唐山地震前后的重力变化[J]. 地震学报, 2(1): 21-31.
	CHEN Yun-tai, GU Hao-ding, LU Zao-xun. 1980. Variation of the gravity before and after the Haicheng earthquake, 1975 and the Tangshan earthquake, 1976[J]. Acta Seismologica Sinica, 2(1): 21-31. (in Chinese)
[2]	顾功叙, 刘克人, 郑金涵, 等. 1997. 中国京津唐张地区时间上连续的重力变化与地震的孕育和发生[J]. 科学通报, 42(18): 1919-1930.
	GU Gong-xu, LIU Ke-ren, ZHENG Jin-han, et al. 1997. Time continuity of gravity variation in the BJTZ region and earthquake development and occurrence[J]. Chinese Science Bulletin, 42(18): 1919-1930. (in Chinese)
[3]	胡敏章, 郝洪涛, 李辉, 等. 2019. 地震分析预报的重力变化异常指标分析[J]. 中国地震, 35(3): 417-430.
	HU Min-zhang, HAO Hong-tao, LI Hui, et al. 2019. Quantitative analysis of gravity changes for earthquake prediction[J]. Earthquake Research in China, 35(3): 417-430. (in Chinese)
[4]	李传友, 张金玉, 王伟, 等. 2021. 2021年云南漾濞6.4级地震发震构造分析[J]. 地震地质, 43(3): 706-721.
	LI Chuan-you, ZHANG Jin-yu, WANG Wei, et al. 2021. The seismogenic fault of the 2021 Yunnan Yangbi M_S6.4 earthquake[J]. Seismology and Geology, 43(3): 706-721. (in Chinese)
[5]	李辉, 付广裕, 孙少安, 等. 2000. 滇西地区重力场动态变化计算[J]. 地壳形变与地震, 20(1): 60-66.
	LI Hui, FU Guang-yu, SUN Shao-an, et al. 2000. Computation on dynamic gravity changes in the western area of Yunnan Province[J]. Crustal Deformation and Earthquake, 20(1): 60-66. (in Chinese)
[6]	李辉, 申重阳, 孙少安, 等. 2009. 中国大陆近期重力场动态变化图像[J]. 大地测量与地球动力学, 29(3): 1-10.
	LI Hui, SHEN Chong-yang, SUN Shao-an, et al. 2009. Dynamic gravity change in recent years in China continent[J]. Journal of Geodesy and Geodynamics, 29(3): 1-10. (in Chinese)
[7]	卢造勋, 方昌流, 石作亭, 等. 1978. 重力变化与海城地震[J]. 地球物理学报, 21(1): 1-8.
	LU Zao-xun, FANG Chang-liu, SHI Zuo-ting, et al. 1978. Variation of the gravity filed and the Haicheng earthquake[J]. Chinese Journal of Geophysics, 21(1): 1-8. (in Chinese)
[8]	申重阳, 李辉, 付广裕. 2003. 丽江7.0级地震重力前兆模式研究[J]. 地震学报, 25(2): 163-171.
	SHEN Chong-yang, LI Hui, FU Guang-yu. 2003. Study on a gravity precursor mode of Lijiang earthquake with M_S=7.0[J]. Acta Seismologica Sinica, 25(2): 163-171. (in Chinese)
[9]	唐伯雄, 吴国华, 罗剑寒. 1990. 红河断裂带上的重力测量[J]. 地震研究, 13(3): 273-276.
	TANG Bo-xiong, WU Guo-hua, LUO Jian-han. 1990. Gravity measurement along Red River fault zone[J]. Journal of Seismological Research, 13(3): 273-276. (in Chinese)
[10]	王林海, 陈石, 庄建仓, 等. 2020. 精密重力测量中相对重力仪格值系数的贝叶斯估计方法[J]. 测绘学报, 49(12): 1543-1553.
	WANG Lin-hai, CHEN Shi, ZHUANG Jian-cang, et al. 2020. Bayesian estimation of the scale factor of relative gravimeter in precise gravity survey[J]. Acta Geodaetica et Cartographica Sinica, 49(12): 1543-1553. (in Chinese)
[11]	王青华, 郝洪涛, 汪健, 等. 2019. 云南地震流动重力监测网建设与映震能力分析[J]. 大地测量与地球动力学, 39(3): 317-324.
	WANG Qing-hua, HAO Hong-tao, WANG Jian, et al. 2019. Construction and analysis of seismic response capability of gravimetric network in Yunnan[J]. Journal of Geodesy and Geodynamics, 39(3): 317-324. (in Chinese)
[12]	隗寿春, 徐建桥, 周江存. 2016. 重力网的分段线性动态平差[J]. 测绘学报, 45(5): 511-520.
	WEI Shou-chun, XU Jian-qiao, ZHOU Jiang-cun. 2016. Piece-wise linear dynamic adjustment for gravity network[J]. Acta Geodaetica et Cartographica Sinica, 45(5): 511-520. (in Chinese)
[13]	吴国华, 罗增雄, 赖群, 等. 1995. 1988年澜沧-耿马地震与滇西实验场的重力变化[J]. 地壳形变与地震, 15(2): 66-73.
	WU Guo-hua, LUO Zeng-xiong, LAI Qun, et al. 1995. Lancang-Gengma earthquake in 1988 and gravitational variation in western Yunnan seismological experiment site(WYSES)[J]. Crustal Deformation and Earthquake, 15(2): 66-73. (in Chinese)
[14]	吴国华, 罗增雄, 赖群. 1997. 丽江7.0级地震前后滇西实验场的重力异常变化特征[J]. 地震研究, 20(1): 103-109.
	WU Guo-hua, LUO Zeng-xiong, LAI Qun. 1997. The variation characteristics of gravity anomaly in the earthquake prediction test site in western Yunnan before and after the M7.0 Lijiang earthquake[J]. Journal of Seismological Research, 20(1): 103-109. (in Chinese)
[15]	邢乐林, 李辉, 夏正超, 等. 2010. CG-5重力仪零漂特性研究[J]. 地震学报, 32(3): 369-373.
	XING Le-lin, LI Hui, XIA Zheng-chao, et al. 2010. Study on zero drift characteristics of CG-5 gravimeter[J]. Acta Seismologica Sinica, 32(3): 369-373. (in Chinese)
[16]	于书媛, 骆佳骥, 杨源源, 等. 2021. InSAR数据约束的2021年5月21日云南漾濞 M_S6.4 地震发震构造研究[J]. 地震工程学报, 43(4): 777-783, 790.
	YU Shu-yuan, LUO Jia-ji, YANG Yuan-yuan, et al. 2021. Seismogenic structure of the Yangbi, Yunnan M_S6.4 earthquake on May 21, 2021 constrained by InSAR data[J]. China Earthquake Engineering Journal, 43(4): 777-783, 790. (in Chinese)
[17]	祝意青, 申重阳, 张国庆, 等. 2018. 我国流动重力监测预报发展之再思考[J]. 大地测量与地球动力学, 38(5): 441-446.
	ZHU Yi-qing, SHEN Chong-yang, ZHANG Guo-qing, et al. 2018. Rethinking the development of earthquake monitoring and prediction in mobile gravity[J]. Journal of Geodesy and Geodynamics, 38(5): 441-446. (in Chinese)
[18]	祝意青, 申重阳, 刘芳, 等. 2020. 重力观测地震预测应用研究[J]. 中国地震, 36(4): 708-717.
	ZHU Yi-qing, SHEN Chong-yang, LIU Fang, et al. 2020. Application of earthquake prediction based on gravity observation[J]. Earthquake Research in China, 36(4): 708-717. (in Chinese)
[19]	Barnes D F. 1966. Gravity changes during the Alaska earthquake[J]. Journal of Geophysical Research, 71(2): 451-456. DOI URL
[20]	Chen S, Liu M, Xing L L, et al. 2016. Gravity increase before the 2015 M_W7.8 Nepal earthquake[J]. Geophysical Research Letters, 43(1): 111-117. DOI URL
[21]	Chen S, Zhuang J C, Li X Y, et al. 2019. Bayesian approach for network adjustment for gravity survey campaign: Methodology and model test[J]. Journal of Geodesy, 93(5): 681-700. DOI
[22]	Fu G Y, Gao S H, Freymueller J T, et al. 2014. Bouguer gravity anomaly and isostasy at western Sichuan Basin revealed by new gravity surveys[J]. Journal of Geophysical Research: Solid Earth, 119(4): 3925-3938. DOI URL
[23]	Fujii Y. 1966. Gravity changes in the shock area of the Niigate earthquake, 16 June 1964[J]. Zisin, 19(3): 200-216. DOI URL
[24]	Kuo J T, Zheng J H, Song S H, et al. 1999. Determination of earthquake epicentroids by inversion of gravity variation data in the BTTZ region, China[J]. Tectonophysics, 312(2): 267-281. DOI URL
[25]	Shen C Y, Li H, Sun S A, et al. 2012. Temporal gravity changes before the 2008 Yutian M_S7.3 earthquake[J]. Geodesy and Geodynamics, 3(1): 19-26. DOI URL
[26]	Zhu Y Q, Liang W F, Zhang S. 2018. Earthquake precursors: Spatial-temporal gravity changes before the great earthquakes in the Sichuan-Yunnan area[J]. Journal of Seismology, 22(1): 217-227. DOI URL
[27]	Zhu Y Q, Liu F, Cao J P, et al. 2012. Gravity changes before and after the 2010 M_S7.1 Yushu earthquake[J]. Geodesy and Geodynamics, 3(4): 1-6.
[28]	Zhu Y Q, Zhan F B, Zhou J C, et al. 2010. Gravity measurements and their variations before the 2008 Wenchuan earthquake[J]. Bulletin of the Seismological Society of America, 100(5B): 2815-2824. DOI URL