On monitoring precursors of major earthquakes with dense network of borehole strainmeters
-
摘要: 提出了地震前兆研究是地震预测研究的基础,其研究的对象是个案. 地震前兆研究受到的最大限制,是对个案的观测严重不足. 龙门山断裂带西南段仍然存在发生强震的危险性,应该不失时机地在那里建立密集的钻孔应变观测网. 四分量钻孔应变仪是我国发明的、 已经可以与地震仪和GPS相提并论的观测仪器. 对地震前兆研究而言,钻孔应变观测在理论和实践两方面都具有优越性. 本文认为: 钻孔应变观测点应该建在应力集中的构造部位; 应该在钻孔选点和仪器安装两方面采取措施以保证观测点建设的成功; 四分量钻孔应变观测最重要的是数据自洽; 钻孔应变观测不必追求很深; 观测点建设不必进行绝对应力测量.Abstract: This paper argues that the earthquake precursor research is the basis for earthquake prediction. Objects of the precursor research are individual cases. The major limitation on the earthquake precursor study comes from insufficient observation cases. Since high potential of major earthquake occurrence still exists along the southwestern Longmenshan fault zone, now it is due time to establish a dense network with four-gauge borehole strainmeters (FGBS) over the area. FGBS is a Chinese invention and has already been compared to seismometers. In the earthquake precursor research, FGBS has advantage both in theory and in practice. We think that an FGBS site should be built in the site with stress concentration. Successful establishment of an observation site depends on good condition of the borehole location and right installation of the instrument. The key point of observation with an FGBS is self-consistency of the data from the four gauges. It is unnecessary to pursue the deep borehole strain observation and to conduct absolute stress measurements in order to build a site of FGBS.
-
Keywords:
- earthquake prediction /
- earthquake precursor /
- borehole strainmeter /
- dense network /
- case study
-
引言
地震矩张量在地震学领域有着重要和广泛的应用,例如用于研究区域应力场、地震成因、震源过程和海啸成因等(赵翠萍等,2008;吴微微等,2014;Hou,Xu,2015).震后快速准确地得到地震矩张量对地震本身的研究、地震灾害的快速评估、震后趋势判定和地震应急救援等方面都具有重要的意义(Pasyanos et al,1996).
目前,对于全球较大地震(MW≥5.0),美国哈佛大学的全球矩心矩张量(global centroid moment tensor, 简写为GCMT)项目组、美国地质调查局(USGS)、东京大学和中国地震局地球物理研究所等机构能够利用全球台网地震记录或大范围的全国台网记录,在震后较短时间内(一般为2小时)给出快速矩张量目录.此外,赵旭等(2014)使用最新构建的虚拟全球地震台,开发了基于W震相技术的全球强震(MW≥6.5)矩心矩张量自动反演系统,该系统能够实现在震后25-40分钟内自动准确测定全球MW≥6.5地震的矩心矩张量.但是,这些研究都利用震中距大于几千千米的远震波形进行分析,同时要求地震具有足够的能量.
20世纪80年代以来,随着测震台网密度的不断提高,高增益、大动态范围的宽频带数字地震仪被广泛应用于地震监测,使用数字化区域宽频带地震波形在区域范围内(震中距 < 1 000 km)获取震源机制的工作得到了很大的发展,并形成了目前常用的两种地震矩张量反演方法:TDMT_INT (time domain seismic moment tensor inversion)方法(Dreger,Helmberger,1993;Romanowicz et al,1993;Benmenahem,1995)和CAP (cut and paste)方法(Zhao,Helmberger,1994;Nábělek,Xia,1995;Randall et al,1995;Zhu,Helmberger,1996).很多研究人员已利用这两种方法,基于区域宽频带地震波形记录,求解了国内部分中等强度地震的矩张量(谢祖军等,2012;陈继峰,2015;康清清等,2015)和强震的中强余震矩张量(吕坚等,2008).
近年来,随着区域台网建设的快速发展和地震矩张量的深入应用,实现区域地震矩张量的自动反演已成为必然的发展趋势.日本国家地球科学与防灾研究所(National Research Institute for Earth Science and Disaster Prevention,简写为NIED)利用区域台网实现了地震矩张量的近实时自动反演(Fukuyama et al,1998;Fukuyama,Dreger,2000;Kubo et al,2002),瑞士地震研究所(Swiss Seismological Service,简写为SED)、安达卢西亚地理学会(Instituto Andaluz de Geofísica,简写为IAG)和西班牙国家地理学会(Instituto Geográfico Nacional,简写为IGN)等机构也已开展了例行的震源矩张量反演(Daniel et al,2003;Bernardi et al,2004;Rueda,Mezcua,2005). 邱毅等(2014)在Dreger和Helmberger (1993)提出的矩张量反演程序的基础上,建立了一套区域中等强度地震的矩张量快速自动反演系统,并对台湾和福建地区的中等强度地震矩张量进行了自动求解.
本文拟结合全国地震自动速报系统,建立基于区域台网宽频带地震记录波形的地震矩张量快速反演系统.地震矩张量反演的核心计算程序来自美国内华达大学地震实验室开发的矩张量反演程序软件包(Ichinose et al, 2003).
1. 方法原理
1.1 矩张量反演
在双力偶点源的假设条件下,矩张量反演方法可表示为(Jost,Herrmann,1989;Kawakatsu,1998;Tajima et al,2002):
(1) 式中, s表示震源,k表示台站,Gisk(t)和dk(t)分别表示震源s至台站k的理论格林函数和台站k的实际观测记录.这里上标k包含台站k的3个分量,mis表示矩张量的第i个分量.若已知地壳速度模型,格林函数可通过频率波数(f-k)法(Wang,Herrmann,1980; Herrmann,Wang,1985; Saikia,1994)计算得到.一旦获得观测数据,利用最小二乘法将观测波形与不同深度的格林函数进行拟合,再根据目标函数VR (方差减小值)来确定最佳解(Dreger,Helmberger,1993;Fukuyama et al,1998).将理论波形记为ok(t), 则
(2) VR值越大,表明波形拟合得越好,求解出来的地震矩张量越可靠.因此我们认为最大VR值对应的深度即为最佳深度.
1.2 矩张量反演程序软件包
内华达地震实验室矩张量反演程序软件包(Ichinose et al, 2003)可基于区域台网三分量地震记录进行全波形地震矩张量反演,反演过程中每一步数据处理都是一个独立的模块,用户可以根据自己的需要对某个模块进行单独控制,最后通过地震矩张量、震源深度及其与发震时刻的关系等判别标准来确定最佳的地震矩张量解.
美国内华达地震实验室已将该矩张量反演程序包用于内华达地区的地震矩张量实时反演中. Ichinose等(2003)曾使用该方法,基于不同的速度结构模型,对美国加利福尼亚东部和内华达州西部MW≥4.0地震进行了矩张量计算,计算结果显示,在方位角分布较为均匀的情况下,最少使用3个台站即可得到较好的矩张量解,且P轴和T轴方向的偏差可控制在10°以内,说明该软件的稳定性有充分保障.
2. 区域台网地震矩张量快速反演系统的技术路线
本研究研发的区域台网地震矩张量快速反演系统包括自动反演平台、手动反演平台和管理平台等3部分.自动反演平台实现了系统触发、地震相关信息与地震波形记录获取、矩张量反演计算和最佳解挑选等过程的自动化;手动反演平台可以根据用户需求,自主进行反演计算;管理平台用于管理各项反演事件及结果.该快速反演系统操作便捷,最终通过网页形式对外提供数据服务.
2.1 自动矩张量计算平台
本文研发的实时快速矩张量自动反演系统的框架如图 1a所示,其处理流程(图 1b)主要分为以下步骤:
1)地震发生后,全国地震速报系统自动触发快速矩张量自动反演系统,并由地震速报系统自动获取地震三要素;
2)自动从区域测震台网实时流数据库中截取地震事件波形记录,并解压缩为后续步骤需要的SAC格式文件,同时准备好反演所需要的台站参数文件;
3)自动进行矩张量反演计算.由于本文使用的内华达大学地震实验室的矩张量反演程序,其求解过程采用人机交互方式,每次反演计算需要人工分多阶段运行多个程序,且各程序运行过程中需要输入的参数和生成的结果也需要人工修改.因此本文对该反演程序进行改进,自动生成各阶段需要的输入文件及自动读取各阶段程序生成的结果,实现自动矩张量解反演.本系统采用Crust2.0速度模型计算格林函数,采用的滤波频段为0.02-0.08 Hz;
4)自动挑选矩张量最佳解.本文程序选用信噪比大于20 dB的地震波形记录,与不同深度理论震动图进行拟合求解矩张量解,最终选取反演方差减小值VR最大的解作为此次地震的矩张量最佳解;
5)输出矩张量最佳解结果,并在网页上进行发布和网页给出的反演结果包括地震三要素、震中位置、台站分布和矩张量解等信息,以及波形拟合结果.
2.2 手动矩张量计算平台
手动计算平台略去了传统矩张量计算中Linux系统操作、脚本修改及使用SAC、rdseed等专业软件处理数据的繁琐过程,用户只需上传波形,输入计算参数,便可进行地震矩张量计算,计算方法与自动计算所采用的方法一致.计算完成后,根据不同深度的波形拟合图可以确定计算结果的准确性,选择是否重新计算.计算结果通过网页形式对外服务.目前手动计算平台只支持标准SEED格式的波形文件.
2.3 管理平台
管理平台主要是对计算服务和自动处理的数据进行管理.用户可在管理平台中设置全国自动地震速报系统的连接参数、自动反演触发条件等,也可校正自动计算结果、手动添加历史地震、清理用户计算文件等.
3. 自动反演结果
本文研发的区域台网快速矩张量反演系统目前已在陕西省地震台网中心在线运行,区域台网覆盖范围包括陕西省区域台网台站及四川、甘肃、新疆、内蒙、青海、西藏、宁夏、云南等省区的部分台站,这些台站为快速反演系统提供了高质量的地震记录.
该系统自2014年8月开始运行,在运行近两年的时间内对西北片区的地震进行监测,共计得到165次地震的地震矩张量解(图 2),其中震级最大的地震为2014年8月3日云南鲁甸MS6.5地震,最小的地震为2014年9月12日四川青川MS3.4地震.这165次地震的地震矩张量解的自动反演结果均在地震发生后15-25分钟内发布.
综合考虑本反演系统使用的数据资料和算法特征,选取求解最优解的VR值对反演结果进行定量化的论证.根据系统算法和源程序开发者Ichinose的经验和建议,选取VR≥50%作为反演结果可靠性的标准.本文对系统自动反演的165次地震的VR值进行统计,其中VR≥50%的反演结果为131个,所占比例为79%;VR < 50%的反演结果为34个,所占比例为21%.对这34个VR < 50%的结果进行对比分析,发现其中有27次地震发生在区域台网覆盖范围的外围或边缘,台站包络性不佳.
将系统自动反演的165次地震中的14次MS≥5.0地震的矩张量解与USGS和中国地震局地球物理研究所的反演结果进行对比,结果如表 1所示, 可以看出各项参数的一致性较高.
表 1 区域台网地震矩张量自动反演系统MS≥5.0地震反演结果对比Table 1. Contrast analysis of MS≥5.0 earthquakes' inversion solutions from regional fast inversion system发震时刻 震中
位置MW 深度
/km节面Ⅰ 节面Ⅱ 资料来源* 年-月-日 时:分 走向/° 倾角/° 滑动角/° 走向/° 倾角/° 滑动角/° 6.0 8 165 90 17 75 73 180 本文 2014-08-03 16:30 云南
鲁甸6.1 14 160 87 5 69 85 177 USGS 6.3 11 165 87 6 74 87 177 地球物理研究所 2014-08-17 06:07 云南
永善5.0 8 318 64 2 227 84 154 本文结果 4.9 8 319 62 0 229 90 152 地球物理研究所 2014-10-01 09:23 四川
越西5.0 8 352 74 20 256 70 163 本文结果 5.0 10 350 79 21 256 69 168 地球物理研究所 5.0 16 299 40 97 110 50 84 本文 2014-10-02 23:56 青海
乌兰5.1 12 299 60 105 90 34 66 USGS 5.1 17 287 47 80 121 44 101 地球物理研究所 5.8 8 234 80 -176 143 86 -10 本文 2014-11-22 16:55 四川
康定5.9 14 231 78 -170 139 80 -12 USGS 6.1 12 233 90 180 143 90 0 地球物理研究所 2014-11-25 23:19 四川
康定5.3 12 151 88 -7 241 83 -178 本文 5.6 12 152 82 -3 242 87 -172 USGS 2015-01-14 13:21 四川
乐山4.6 10 292 14 68 135 77 95 本文 4.6 12 151 47 80 345 44 101 地球物理研究所 2015-02-22 14:43 新疆
沙湾4.8 8 242 18 58 95 74 100 本文 4.7 17 57 45 65 270 50 113 地球物理研究所 2015-04-15 15:39 内蒙古
阿拉善
左旗5.1 4 358 84 -164 267 74 -6 本文 5.4 22 174 85 -169 83 79 -5 USGS 5.0 20 140 43 83 330 48 97 本文 2015-06-25 03:12 新疆
托克逊5.0 33 133 46 98 302 45 82 USGS 5.0 31 154 58 88 338 32 93 地球物理研究所 4.8 12 108 83 28 14 62 172 本文 2015-10-12 18:04 青海
玛多5.2 20 288 81 5 197 85 171 USGS 5.1 20 112 81 6 21 84 171 地球物理研究所 2015-11-23 05:02 青海
祁连5.0 20 10 64 151 114 64 30 本文 5.0 9 11 78 143 110 54 15 地球物理研究所 4.9 20 242 39 27 131 74 125 本文 2016-01-14 05:18 新疆
轮台5.1 24 273 28 85 99 62 93 USGS 5.0 20 272 23 78 105 68 95 地球物理研究所 5.7 12 342 42 109 137 51 74 本文 2016-01-21 01:13 青海
门源5.9 7 342 55 101 144 36 75 USGS 6.0 7 335 53 98 141 38 79 地球物理研究所 * USGS的反演结果引自http://earthquake.usgs.gov/earthquakes/search/, 中国地震局地球物理研究所的反演结果引自http://www.cea-igp.ac.cn/tpxw/index.shtml/. 4. 讨论与结论
本文结合区域台网现状,对美国内华达大学地震实验室的矩张量反演程序进行改进,研发了一套基于区域台网宽频带地震记录波形的快速矩张量反演系统.该系统以矩张量自动反演为主,由地震速报系统触发,自动获取地震三要素和地震波形记录,自动反演计算地震矩张量并挑选最佳解,最后通过网页发布,实现了整个处理过程的自动化.这种自动化有利于快速得到矩张量反演结果,能有效提高震后应急服务产品的产出速度,并为灾后烈度速报、紧急救援和震后趋势的快速判定提供数据支持.此外,该系统还提供了一套手动矩张量计算平台和管理平台,只需通过网页即可完成用户选定事件的矩张量计算和管理,略去了各种程序操作和脚本修改的过程,更加方便、快捷,易于操作.该系统目前在陕西省地震台网中心运行,已实现西北部分地区中等强度地震的自动快速反演.
区域台网地震矩张量快速反演系统使用美国内华达地震实验室的矩张量反演程序包,该软件包有充分的稳定性保障,即使在稀疏台网下亦能够得到较为精确的地震矩张量解.对系统自动计算的165次地震反演结果的VR值进行统计,结果表明,系统可靠性可达79%,而其余21%的“不可靠”结果主要是由台网覆盖不全面造成的.如果只考虑区域台网覆盖范围内的地震,系统可靠率可达95%.再将系统自动计算的14次MS≥5.0地震的矩张量解与国内外权威研究机构给出的结果进行对比分析,各项参数的一致性较高,证实了该方法的可靠性比较高.
地震发生后,经速报系统触发,区域台网地震快速矩张量自动反演系统能在20分钟左右准确计算地震的矩张量解.如果参与反演的台站数据较多,计算格林函数时间变长,反演所用时间也会增加.该系统目前处于上线试运行阶段,在今后的工作中还可对本文的方法进行优化,以进一步减少得到结果所需的时间.
架设于陕西省、四川省和甘肃省台网中心的区域台网地震快速矩张量自动反演系统,经过一年多的试运行,成功反演了3.5≤MS < 6.5地震165次.除去3次系统升级时间段内发生的地震外,该自动反演系统均经速报系统触发,成功反演出地震矩张量,且无地震事件遗漏.因此,可以认为该系统能够稳定地、全面地监控区域地震.
研究中使用了美国内华达州地震实验室的矩张量反演软件包和SAC波形处理软件、中国地震局地球物理研究所和美国地质调查局(USGS)的相关计算结果;中国地震局地球物理研究所国家测震台网数据备份中心(doi:10.11998/SeisDmc/SN;郑秀芬等,2009)和甘肃、内蒙古、宁夏、青海、四川、陕西、新疆、西藏、云南各省区的地震台网为本研究提供地震波形数据;两名审稿专家提出了非常宝贵的意见及建议,作者在此一并表示感谢.
-
-
毕永强, 郑国栋, 武成智, 韩忠民. 2013. 敦化地震台2010年9月钻孔应变观测曲线异常分析[J]. 地震地磁观测与研究, 34(1/2): 140-144. Bi Y Q, Zheng G D, Wu C Z, Han Z M. 2013. Discussion of the abnormal curve record by borehole strainmeter at Dunhua in September 2010[J]. Seismological and Geomagnetic Observation and Research, 34(1/2): 140-144 (in Chinese).
陈颙. 1981. 不同应力途径三轴压缩下岩石的声发射[J]. 地震学报, 3(1): 41-48. Chen Y. 1981. Acoustic emission of rocks under triaxial compression along various stress paths[J]. Acta Seismologica Sinica, 3(1): 41-48 (in Chinese).
陈沅俊, 杨修信. 1990. 双衬套钻孔应变测量的计算[J]. 华北地震科学, 8(4): 80-89. Chen Y J, Yang X X. 1990. Calculations of strain observation in the borehole with two cases[J]. North China Earthquake Sciences, 8(4): 80-89 (in Chinese).
陈运泰. 2008. 地震预测要知难而进[J]. 求是, (15): 58-60. Chen Y T. 2008. Research on earthquake prediction should advance despite difficulties[J]. Qiu Shi, (15): 58-60 (in Chinese).
陈运泰. 2009. 地震预测: 回顾与展望[J]. 中国科学: D辑, 39(12): 1633-1658. Chen Y T. 2009. Earthquake prediction: Retrospect and prospect[J]. Science in China: Series D, 39(12): 1633-1658 (in Chinese).
陈运泰, 杨智娴, 张勇, 刘超. 2013. 从汶川地震到芦山地震[J]. 中国科学: 地球科学, 43(6): 1064-1072. Chen Y T, Yang Z X, Zhang Y, Liu C. 2013. From 2008 Wenchuan earthquake to 2013 Lushan earthquake[J]. Scientia Sinica Terrae, 43(6): 1064-1072 (in Chinese).
池顺良, 池毅, 邓涛, 廖成旺, 唐小林, 池亮. 2009. 从5·12汶川地震前后分量应变仪观测到的应变异常看建设密集应变观测网络的必要性[J]. 国际地震动态, (1): 1-13. Chi S L, Chi Y, Deng T, Liao C W, Tang X L, Chi L. 2009. The necessity of building national strain-observation network from the strain abnormality before Wenchuan earthquake[J]. Recent Developments in World Seismology, (1): 1-13 (in Chinese).
池顺良, 刘琦, 池毅, 邓涛, 廖成旺, 阳光, 张贵萍, 陈洁. 2013. 2013年芦山MS7.0地震的震前及临震应变异常[J]. 地震学报, 35(3): 296-303. Chi S L, Liu Q, Chi Y, Deng T, Liao C W, Yang G, Zhang G P, Chen J. 2013. Borehole strain anomalies before the 20 April 2013 Lushan MS7.0 earthquake[J]. Acta Seismologica Sinica, 35(3): 296-303 (in Chinese).
丁国瑜, 梅世蓉, 马宗晋. 1981. 地震预报方法[C]//丁国瑜, 马宗晋编. 国际地震预报讨论会论文选. 北京: 地震出版社: 174-179. Ding G Y, Mei S R, Ma Z J. 1981. Methods of earthquake prediction[C]//Ding G Y, Ma Z J eds. Selected Papers of the International Symposium on Earthquake Prediction. Beijing: Seismological Press: 174-179 (in Chinese).
刘力强, 马胜利, 马瑾, 雷兴林, 楠濑勤一郎, 西泽修, Jouniaux L. 1999. 岩石构造对声发射统计特征的影响[J]. 地震地质, 21(4): 377-386. Liu L Q, Ma S L, Ma J, Lei X L, Kusunose K, Nishizawa O, Jouniaux L. 1999. Effect of rocks structure on the statistic characteristics of acoustic emission[J]. Seismology and Geology, 21(4): 377-386 (in Chinese).
骆鸣津, 顾梦林, 李安印, 睢建设. 1989. 用引潮力进行钻孔应变-应力的原地标定[J]. 地壳形变与地震, 9(4): 51-62. Luo M J, Gu M L, Li A Y, Sui J S. 1989. In-situ calibration of borehole strain-stress with tidal force[J]. Crustal Deformation and Earthquake, 9(4): 51-62 (in Chinese).
欧阳祖熙. 1977. RDB-1型电容式地应变计[C]//全国地应力专业会议论文选编(下): 337-348. Ouyang Z X. 1977. RDB-1 type capacitance ground strainmeter[C]//Selected Papers of the National Conference on Stress Measurement, Part 2: 337-348 (in Chinese).
欧阳祖熙, 张宗润. 1988. 钻孔应变仪与井壁耦合方法的研究[G]//中国地震局地壳应力研究所编. 地壳构造与地壳应力文集(2). 北京: 地震出版社: 1-10. Ouyang Z X, Zhang Z R. 1988. Studies of method for coupling straingauge with borehole wall[G]//Institute of Crustal Dynamics, China Earthquake Administration eds. Bulletin of the Institute of Crustal Dynamics (2). Beijing: Seismological Press: 1-10 (in Chinese).
潘立宙. 1977. 关于地应力测量的几个公式的推导及讨论[C]//全国地应力专业会议论文选编(上): 1-41. Pan L Z. 1977. On the formulae of ground stress measurement[C]//Selected Papers of the National Conference on Stress Measurement, Part 1: 1-41 (in Chinese).
邱泽华. 2010a. 关于地震前兆的判据问题[J]. 大地测量与地球动力学, 30(增刊Ⅱ): 1-5. Qiu Z H. 2010a. How to identify earthquake precursor[J]. Journal of Geodesy and Geodynamics, 30(Suppl.Ⅱ): 1-5 (in Chinese).
邱泽华. 2010b. 中国分量钻孔地应力-应变观测发展重要事件回顾[J]. 大地测量与地球动力学, 30(5): 42-47. Qiu Z H. 2010b. A review of the continuous borehole observation of stress-strain in China[J]. Journal of Geodesy and Geodynamics, 30(5): 42-47 (in Chinese).
邱泽华, 池顺良. 2013. YRY-4型钻孔应变仪观测的P波剪应变[J]. 地震, 33(4): 64-70. Qiu Z H, Chi S L. 2013. Shear strains of P-wave observed with an YRY-4 borehole strainmeter[J]. Earthquake, 33(4): 64-70 (in Chinese).
邱泽华, 石耀霖. 2004. 观测应变阶在地震应力触发研究中的应用[J]. 地震学报, 26(5): 481-488. Qiu Z H, Shi Y L. 2004. Application of observed strain steps to the study of remote earthquake stress triggering[J]. Acta Seismologica Sinica, 26(5): 481-488 (in Chinese).
邱泽华, 石耀霖, 欧阳祖熙. 2005a. 四分量钻孔应变观测的实地相对标定[J]. 大地测量与地球动力学, 25(1): 118-122. Qiu Z H, Shi Y L, Ouyang Z X. 2005a. Relative in-situ calibration of 4-component borehole strainmeters[J]. Journal of Geodesy and Geodynamics, 25(1): 118-122 (in Chinese).
邱泽华, 石耀霖, 欧阳祖熙. 2005b. 四分量钻孔应变观测的实地绝对标定[J]. 地震, 25(3): 27-34. Qiu Z H, Shi Y L, Ouyang Z X. 2005b. Absolute in-situ calibration of 4-component borehole strainmeters[J]. Earthquake, 25(3): 27-34 (in Chinese).
邱泽华, 马瑾, 池顺良, 刘厚明. 2007. 钻孔差应变仪观测的苏门答腊大地震激发的地球环型自由振荡[J]. 地球物理学报, 50(3): 797-805. Qiu Z H, Ma J, Chi S L, Liu H M. 2007. Earth's free torsional oscillations of the great Sumatra earthquake observed with borehole shear strainmeter[J]. Chinese Journal of Geophysics, 50(3): 797-805 (in Chinese).
邱泽华, 张宝红, 池顺良, 唐磊, 宋茉. 2010. 汶川地震前姑咱台观测的异常应变变化[J]. 中国科学: 地球科学, 40(8): 1031-1039. Qiu Z H, Zhang B H, Chi S L, Tang L, Song M. 2011. Abnormal strain changes observed at Guzan before the Wenchuan earthquake[J]. Science China: Earth Sciences, 54(2): 233-240.
邱泽华, 唐磊, 张宝红, 宋茉. 2012. 用小波-超限率分析提取宁陕台汶川地震体应变异常[J]. 地球物理学报, 55(2): 538-546. Qiu Z H, Tang L, Zhang B H, Song M. 2012. Extracting anomaly of the Wenchuan earthquake from the dilatometer recording at NSH by means of wavelet-overrun rate analysis[J]. Chinese Journal of Geophysics, 55(2): 538-546 (in Chinese).
邱泽华, 张国宏, 张宝红, 唐磊, 宋乃波, 郭燕平, 周龙寿, 赵京轶, 苏文浩. 2013. 用高密度电法勘探预判钻孔应变观测点岩石完整性[J]. 地震地质, 35(4): 805-816. Qiu Z H, Zhang G H, Zhang B H, Tang L, Song N B, Guo Y P, Zhou L S, Zhao J Y, Su W H. 2013. Prospecting underground condition for borehole strain-monitoring location by means of resistivity imaging survey[J]. Seismology and Geology, 35(4): 805-816 (in Chinese).
苏恺之. 1977. 地应力相对测量方法[C]//全国地应力专业会议论文选编(上): 42-61. Su K Z. 1977. Methods of relative measurement of ground stress[C]//Selected Papers of the National Conference on Stress Measurement, Part 1: 42-61 (in Chinese).
苏恺之. 1982. 液位型体积式应变仪的工作原理[J]. 地震科学研究, 4: 57-62. Su K Z. 1982. How does hydraulic dilatometer work?[J]. Research of Earthquake Science, 4: 57-62 (in Chinese).
唐磊, 邱泽华, 阚宝祥. 2007. 中国钻孔体应变台网观测到的地球球型振荡[J]. 大地测量与地球动力学, 27(6): 37-44. Tang L, Qiu Z H, Kan B X. 2007. Earth's spheroidal oscillations observed with China borehole dilatometer network[J]. Journal of Geodesy and Geodynamics, 27(6): 37-44 (in Chinese).
Armstrong B H. 1969. Acoustic emission prior to rockbursts and earthquakes[J]. Bull Seismol Soc Am, 59(3): 1259-1279.
Bakun W H, Aagaard B, Dost B, Ellsworth W L, Hardebeck J L, Harris R A, Ji C, Johnston M J S, Langbein J, Lienkaemper J J, Michael A J, Murray J R, Nadeau R M, Reasenberg P A, Reichle M S, Roeloffs E A, Shakal A, Simpson R W, Waldhauser F. 2005. Implications for prediction and hazard assessment from the 2004 Parkfield earthquake[J]. Nature, 437(7061): 969-974.
Geller R J, Jackson D D, Kagan Y Y, Mulargia F. 1997. Earthquakes cannot be predicted[J]. Science, 275(5306): 1616-1617.
Gladwin M T. 1984. High precision multi-component borehole deformation monitoring[J]. Rev Sci Instru, 55: 2011-2016.
Haimson B C. 2010. The effect of lithology, inhomogeneity, topography, and faults, on in situ stress measurements by hydraulic fracturing, and the importance of correct data interpretation and independent evidence in support of results[C]//International Symposium on In-Situ Rock Stress. Salzburg: International Society for Rock Mechanics: 11-14.
Hudson J A, Feng X T. 2010. Variability of in situ rock stress[C]//International Symposium on In-Situ Rock Stress. Salzburg: International Society for Rock Mechanics: 3-10.
Ishii H. 2001. Development of new multi-component borehole instrument[J]. Report of Tono Research Institute of Earthquake Science, 6: 5-10 (in Japanese).
Johnston M J S, Borcherdt R D, Linde A T, Gladwin M T. 2006. Continuous borehole strain and pore pressure in the near field of the 28 September 2004 M6.0 Parkfield, California, earthquake: Implications for nucleation, fault response, earthquake prediction, and tremor[J]. Bull Seismol Soc Am, 96(4B): S56-S72.
Jordan T H, Chen Y T, Gasparini P, Madariaga R, Main I, Marzocchi W, Papadopoulos G, Sololev G, Yamaoka K, Zschau J. 2011. Operational earthquake forecasting: State of knowledge and guidelines for utilization[J]. Ann Geophys, 54(4): 315-391.
Linde A T, Gladwin M T, Johnston M J S, Gwyther R L, Bilham R G. 1996. A slow earthquake sequence near San Juan Bautista, California in December 1992[J]. Nature, 383: 65-69.
Mogi K. 1982. Earthquake as acoustic emission: 1980 IZU Peninsula earthquake in particular[J]. J Acoust Emission, 1(1): 37-44.
Qiu Z H, Tang L, Zhang B H, Guo Y P. 2013a. In situ calibration of and algorithm for strain monitoring using four-gauge borehole strainmeters (FGBS)[J]. J Geophys Res, 118(4): 1609-1618.
Qiu Z H, Yang G, Tang L, Guo Y P, Zhang B H. 2013b. Abnormal strain changes observed by a borehole strainmeter at Guzan before the MS7.0 Lushan earthquake[J]. Geodesy and Geodynamics, 4(3): 19-29.
Sacks I S, Evertson D W. 1971. A sensitive borehole strain-rate meter[J]. Carnegie Inst Washington Yearbook, 68: 448-455.
Sacks I S, Linde A T, Suyehiro S, Snoke J A. 1978. Slow earthquakes and stress redistribution[J]. Nature, 275: 599-602.
Stephansson O, Zang A. 2010. How to generate the Final Rock Stress Model (FRSM) at a site or an area[C]//International Symposium on In-Situ Rock Stress. Salzburg: International Society for Rock Mechanics: 15-26.
Wyss M. 1991. Evaluation of Proposed Earthquake Precursors[M]. Washington D C: American Geophysical Union: 1-94.
Wyss M. 1997a. Cannot earthquakes be predicted?[J]. Science, 278(5337): 487-490.
Wyss M. 1997b. Second round of evaluations of proposed earthquake precursors[J]. Pure Appl Geophys, 149(1): 3-16.
-
期刊类型引用(1)
1. 何苗, 吴立新, 崔静, 王威, 齐源, 毛文飞, 苗则朗, 陈必焰, 申旭辉. 汶川地震前多圈层短—临遥感异常回顾及其时空关联性. 遥感学报. 2020(06): 681-700 . 百度学术
其他类型引用(1)