基金项目:中国地震局“三结合”基金项目(2008042)资助.
(1.中国地震局兰州地震研究所,兰州 730000; 2.甘肃省地震局武山地震台,甘肃 天水 741308; 3.甘肃省地震局,兰州 730000; 4.河北省地震局,石家庄050021)
(1.Lanzhou Institute of Seismology,CEA,Lanzhou 730000,Gansu,China)(2.Wushan Seismic Station,Earthquake Administration of Gansu Province,Tianshui 741308,Gansu,China)(3.Earthquake Administration of Gansu Province,Lanzhou 730000,Gansu,China)(4.Earthquake Administration of Hebei Province,Shijiazhuang 050021,Hebei,China)
Radon; subordinate function; variation rate; trend rate; precursory anomaly; Southeast Gansu Province
备注
基金项目:中国地震局“三结合”基金项目(2008042)资助.
利用从属函数、变差率以及趋势速率方法提取甘肃东南部水氡观测资料的地震前兆异常。结果 表明,1984年1月至2009年12月,7项水氡资料共出现47次从属函数异常、45次变差率异常和39次趋势速率异常。其中,33次从属函数异常、35次变差率异常和30次趋势速率异常很好地对应了甘肃及邻区震中距在33~585 km范围内发生的5~8级地震。通过分析异常特征及其与地震的关系,提出中强震和强震时间预测的指标与方法,得出甘肃东南部水氡从属函数计算的经验常数,以及变差率和趋势速率的异常阈值。
Using the methods of subordinate function,variation rate and trend rate,we extract the earthquake precursory anomalies from observational data of Radon in southeast Gansu Province.The result indicates that 7 items of radon data include 47 subordinate-function anomalies,45 variation-rate anomalies and 39 trend-rate anomalies from January 1984 to December 2009.Among them,33 subordinate-function anomalies,35 variation-rate anomalies and 30 trend-rate anomalies corresponded with the earthquakes MS5-8 that ever hit Gansu Province and the adjacent region within 33 to 585 km.By the analysis of the features of these anomaies and their relation to the earthquakes,the method of time prediction for the moderate-strong and strong earthquake were proposed,and the empirical constant,variation rate and threshold value of trend-rate anomaly are obtained from Radon in southeast Gansu Province.
引言
甘肃处于青藏块体的东北缘,祁连山地震带和南北地震带横穿境内,是我国大陆强震活动频繁的场所之一。该地区历史上曾发生过多次强烈地震,现今地震活动仍很频繁,地震震源浅,频度与强度高,危害性强。水氡测项是地震预测的主要依据之一。甘肃东南部水氡测点多,具有良好的地震水文地球化学条件,观测资料连续可靠(张新基等,2005; 杨兴悦等,2008)。为了更好地发挥其重要作用,对其观测资料进行分析总结是十分必要的。
由于流体资料或多或少地受仪器、室温等因素的影响,采用一般的“看图识字”方法,很难将地震前兆异常与观测系统本身引起的变化区别开来(杨明波等,2006; 范雪芳等,2007,2009)。为此,在本研究中,笔者使用从属函数、变差率及趋势速率3种数学方法对观测数据进行处理,以提取明显的地震异常(杨兴悦等,2006; 张立等,2006),进而系统地剖析前兆特征,总结异常变化规律。
1 前兆资料及地震目录的选取
甘肃省地下流体观测台站主要集中在甘肃东南部(图1),因此笔者选用甘肃东南部的流体资料作为研究对象。在资料筛选过程中,采用了王吉易等(2003)的“同一化” 研究思路,并遵循4个原则:① 观测数据时间长,时间在10 a以上; ② 无人为干扰(包括仪器改变、台站改造、井(泉)装置改变、无规律抽水、渠道放水与停水等; )引起的资料台阶状变化或对台阶状变化因素比较清楚,能进行校正,以及无其它变化; ③ 具有明显的年动态变化; ④ 测点目前仍在正常观测。
经过严格筛选,研究区共有7个测点的水氡观测资料满足以上4个条件,分别为武山1号泉、武
图1 甘肃东南部地下流体监测台站分布图
Fig.1 Distribution of underground fluid monitoring stations in southeast Gansu Province考虑本区地震活动的周期性和强度起伏的特点,选取距测点33~594 km内的甘肃及其邻区5级以上地震作为研究对象。在本文中,所研究的5.0、5.2级地震的最大震中距在200 km左右; 5.8~6.2级4次地震的震中距在300 km左右; 共和7.0、6.0级地震和汶川8.0级地震,个别测点的震中距大于500 km。依据本文对资料的分析认为,甘肃东南部水氡资料对共和6级及以上地震是有反应的。为方便统计,相邻两次地震时间间隔小于1 a的,只统计发生在前的地震。所选地震的参数见表1。
表1 1987~2009年甘肃及邻区发生的MS≥5.0地震
Tab.1 MS≥5.0 Earthquakes in Gansu Province and the adjacent regions from 1987 to 20092 异常提取方法
2.1 从属函数法前兆观测资料的异常变化形态虽然是多种多样的,诸如突跳、转折、加速等,但其本质都是观测曲线随时间变化的斜率变化的量。从属函数就是表征一条曲线变化的斜率和跳动程度(贺仲雄,1983; 叶秀薇,2004)。通过对资料的分析比较,笔者认为,对于甘肃地区,用13点滑动平均值作为数据序列计算的效果较好,明显优于直接使用月均值提取异常的方法(杨兴悦等,2006)。中期异常判定方法见表2。
表2中,ki为观测值变化斜率; ri为滑动平均值M(t)与时间t的相关系数,它反应了观测序列内在质量的好坏; α为经验常数。在计算时,13点滑动平均值、从属函数值μi均置于计算时段的右端点。
2.2 变差率法在中短期异常判定中,观测值破年变形态被认为是地下流体参量的重要异常现象。变差率是表征一条曲线的相对变化幅度。以往分析人员只是通过多年曲线特征对比来分析破年变异常,但不能给出“破年变”的“程度”。为此,笔者引入了月变差率方法来定量描述这类异常,也就是用定量值来确定本年相对于上年在同一月的变化程度,如果出现明显的变化差,则表明破年变程度较大,若变差率较小,则表明本年与上年的年动态相似 刘耀炜,杜建国.1990.地下水化学参量中短期异常识别、标志、体系及预报方法研究.。设定了合适的阈值,就可以判定异常的起始时间和持续时间。变差率方法适合于分析有年变的观测资料,但在本文中使用的资料连续时间较长,有的年份资料年变明显,有的年份无年变,因此在本研究中笔者先对资料进行平滑处理后再作分析。其具体做法如下:
(1)月均值的平滑滤波。取水氡等测项的月均值作为基本的数据序列。为抑制水氡的年变化和短期异常变化,突出中短期异常,我们采用13点滑动平均法,求出月均值的滑动平均值。
(2)变差率的计算。定义变差率为某个月的水氡值与前一年同月的水氡值之差再除以前一年同月的水氡值。该参量用以作为单点的中期和中短期异常的分析。其计算式为
R^-i=(C^-i-C^-i-12)/(C^-i-12)·(1)
式中,R^-i为变差率; C^-i和C^-i-12分别为水氡月滑动平均值和一年前同月的值。
在计算时,13点滑动平均值、变差率R^-i值均置于计算时段的右端点。
计算结果表明先平滑后再使用变差率方法提取异常的效果好于直接使用变差率。
2.3 趋势速率法趋势变化特征是水化学参量中期和中短期异常的主要特征之一。由于不同的流体测点所处的构造、含水层岩性、井(泉)结构等物理化学性质的差异,趋势性异常就有上升、下降和波动之分。趋势上升型中期异常和转折型中短期异常是强震及大震的前兆异常。本文中提取水氡趋势速率异常时采用以加速为主的上升异常。尽管趋势速率方法与从属函数方法都表示某一时段内时间序列斜率大小,都是对线性相关性综合特征的描述,但两种方法得到的结果不完全相同,在资料分析运用中可以作为参考。具体做法是:
(1)13点滑动平均值的计算。以月均值为数据序列做计算,计算值放在右端点上。
(2)趋势速率值(V)的计算。选择线性窗Y=a+bX拟合模型,拟合窗口长度为13。
(3)异常判定指标。若趋势速率值Vi≥0.5,则为异常。
3 典型实例
甘肃东南部部分测点水氡月均值在本文所选震例的震前出现了明显的异常,有的测点的异常却不明显(图2),不管明显与否,都是定性描述,因人而异,带有经验性,且很难在时间尺度上做出预测。因此,笔者使用从属函数、变差率及趋势速率3种数学方法对观测数据进行处理,提取的异常较为明显,并对异常进行定量分析。
武山温泉22号井位于武山县温泉乡,海拔1 731 m。构造上处于西秦岭北缘深大断裂带(西起青海循化,东至陕西宝鸡一带)的温泉—甘泉深断裂、聂河断裂及马长庄断裂的交汇部位。含水岩层为印支期似斑状花岗岩,裂隙发育。水质类型为微矿化度重碳酸钠型。该井为自流井,井深154.39 m,水温56.3 ℃,水氡为主要观测项目之一(杨兴悦等,2008)。
图3为武山22号井水氡月均值及异常曲线图。从图3可见,在观测井周边5级以上地震发生前,难以直接从月均值曲线上提取出异常信息,并且对发震时间也无法作出定量预测,但利用从属函数、变差率以及趋势速率方法提取的异常信息则较明显,能定量地给出时间特征参量。以2006年6月21日甘肃文县MS5.0地震为例,从属函数、变差率、趋势速率异常持续时间分别为16、23及16个月,异常结束至发震时间分别为10、5及8个月。
图2 甘肃东南部水氡月均值曲线(a)武山1号泉;(b)武都殿沟;(c)武山22号井; (d)平凉附件厂;(e)通渭温泉;(f)西和芦河; (g)清水温泉
Fig.2 Curves of month mean values of Radon in southeast Gansu Province (a)No.1 spring of Wushan station;(b)Diangou observing spot of Wudu station;(c)observing spot of No.22 well of Wushan station;(d)Fujianchang observing spot of Pingliang station;(e)observing spot of Tongwei spring;(f)Luhe observing spot of Xihe station;(g)Qingshui spring4 前兆异常特征
4.1 从属函数中期异常特点笔者用从属函数法分析1984年1月至2009年12月甘肃东南部7个测项的水氡观测资料,统计结果表明,统计时段内共出现47次异常(图4),其中33次异常很好地对应了表1中所列的甘肃及其邻区距测点33~594 km内的5~8级地震。发生了地震而没有出现异常就算漏报,共18次; 出现了异常而没发生地震属虚报,共14次。中期异常特征列于表3。通过对甘肃东南部地区水氡资料的研究,笔者得到水氡从属函数中期异常的特点:
(1)持续时间:最短为2个月,最长为20个月,平均为8.6个月。
(2)超前时间Ⅰ:最短为1个月,最长33个月,平均为15.6个月。
(3)超前时间Ⅱ:最短为-11个月,最长24个月,平均为7.0个月。
(4)在异常持续过程中发生地震9次,异常结束后发生地震24次。
图4 甘肃东南部水氡从属函数异常曲线(a)武山1号泉;(b)武都殿沟;(c)武山22号井; (d)平凉附件厂;(e)通渭温泉;(f)西和芦河; (g)清水温泉
Fig.4 Abnormal curves of subordinate function of Radon in southeast Gansu Province (a)No.1 spring of Wushan station;(b)Diangou observing spot of Wudu station;(c)observing spot of No.22 well of Wushan station;(d)Fujianchang observing spot of Pingliang station;(e)observing spot of Tongwei spring;(f)Luhe observing spot of Xihe station;(g)Qingshui spring(6)异常对应率为70.2%,漏报率为35.3%,虚报率为29.8%。
4.2 变差率中期异常特点笔者用变差率方法分析了1984年1月至2009年12月甘肃东南部7个测项的水氡观测资料,统计结果表明,共出现45次异常,2009年出现的两次异常视为新异常未统计(图5),其中35次异常
表3 甘肃东南部水氡从属函数中期异常特征统计表
Tab.3 Feature of the mid-term anomalies of Suberdinate function of Radon in Southeast Gansu Province注:超前时间Ⅰ系异常开始时间与发震时间之间的时差; 超前时间Ⅱ系异常结束时间与发震时间之间的时差。“—”表示异常未结束,例如“-1”表示地震后1月异常结束.
很好地对应了甘肃及邻区距测点33~594 km的5~8级地震,其中漏报16次,虚报10次。中期异常特征见表4。通过对甘肃东南部水氡资料的研究,笔者总结出水氡变差率中期异常的特点:
(1)持续时间:最短为2个月,最长为51个月,平均为16.9个月。
(2)超前时间Ⅰ:最短为0个月,最长35个月,平均为14.4个月。
(3)超前时间Ⅱ:最短为-39个月,最长25
图5 甘肃东南部水氡变差率异常曲线(a)武山1号泉;(b)武都殿沟;(c)武山22号井; (d)平凉附件厂;(e)通渭温泉;(f)西和芦河; (g)清水温泉
Fig.5 Abnormal curves of variation rate of Radon in southeast area of Gansu Provinceobserving spot of Xihe station;(g)Qingshui spring
个月,平均为0.5个月。
(4)在异常持续过程中发生地震18次,异常结束后发生地震17次。
(5)异常变幅:最小为1%,最大为12%,平均为4.4%。
(6)异常对应率为77.8%,漏报率为31.4%,虚报率为22.2%。
4.3 趋势速率中期异常特点笔者用趋势速率方法分析了1984年1月至
图6 甘肃东南部水氡趋势速率异常曲线(a)武山1号泉;(b)武都殿沟;(c)武山22号井; (d)平凉附件厂;(e)通渭温泉;(f)西和芦河; (g)清水温泉
Fig.6 Abnormal curves of trend rate of Radon in southeast area of Gansu Province (a)No.1 spring of Wushan station;(b)Diangou observing spot of Wudu station;(c)observing spot of No.22 well of Wushan station;(d)Fujianchang observing spot of Pingliang station;(e)observing spot of Tongwei spring;(f)Luhe observing spot of Xihe station;(g)Qingshui spring表4 甘肃东南部水氡变差率中期异常特征统计表
Tab.4 Features of the mid-term anomalies of variation rate of Radon in southeast Gansu Province注:超前时间Ⅰ系异常开始时间与发震时间之间的时差; 超前时间Ⅱ系异常结束时间与发震时间之间的时差。“—”表示异常未结束,例如“-1”表示地震后1月异常结束.
表5 甘肃东南部水氡趋势速率中期异常特征统计表
Tab.5 Features of the mid-term anomalies of trend rate of Radon in southeast Gansu Province注:超前时间Ⅰ系异常开始时间与发震时间之间的时差; 超前时间Ⅱ系异常结束时间与发震时间之间的时差。“—”表示异常未结束,例如“-1”表示地震后1月异常结束.
2009年12月甘肃东南部7个测项的水氡观测资料,发现该时段内共出现39次异常(2009年出现的一次异常视为新异常未统计)(图6),其中30次异常很好地对应了甘肃及邻区距测点33~594 km内的5~8级地震,其中漏报21次,虚报9次。中期异常特征见表5。通过对这些资料的研究,笔者总结出甘肃东南部水氡趋势速率中期异常的特点:
(1)持续时间:最短为1个月,最长为28个月,平均为9.7个月。
(2)超前时间Ⅰ:最短为0个月,最长33个月,平均为13.0个月。
(3)超前时间Ⅱ:最短为-17个月,最长21个月,平均为3.3个月。
(4)在异常持续过程中发生地震11次,异常结束后发生地震18次。
(5)异常变幅:最小为0.5,最大为1.00,平均为0.78。
(6)对应率为76.9%,漏报率为41.2%,虚报率为23.1%。
5 中期时间预测指标
依据甘肃东南部水氡的多次异常资料,可得到中期异常的超前时间。将异常开始时间至发震时间的差记为超前时间Ⅰ,将异常结束时间至发震时间的差记为超前时间Ⅱ,从表3、表4和表5中可得出各个测点的超前时间Ⅰ和超前时间Ⅱ。由于本文统计的是一个地区的资料,观测点分散,异常出现和结束的时间差别较大,因此很难得出一个统一的预测指标,但可根据各测点的超前时间Ⅰ的最小值和超前时间Ⅱ的最大值预测发震时间段。
当从属函数等异常达到预测指标时,可按下式算出发震时间段:
T1=t1+ΔT1,(1)
T2=t2+ΔT2.(2)
式中,T1为预测的发震开始时间; t1为异常开始时间; ΔT1为超前时间Ⅰ的最小值; T2为预测的发震结束时间; t2为异常结束时间,ΔT2为超前时间Ⅱ的最大值。
当一个测点出现异常后,可用上述方法得出一个发震时间段,如果多个测点出现异常就可得到多个发震时间段。将多个发震时间段进行优化,通过验证分析,笔者认为舍去1/2以下测点较早的发震开始时间和较晚的发震结束时间,这样预测的发震时间段就大大缩短了。
同样,变差率、趋势速率等异常达到预测指标后可用上述方法优化后再对这几种方法得到的发震时段进一步优化。对应的发震时段越靠近预测发震时段的中央就越可能是发生地震的时间。
下面以2006年6月21日甘肃文县MS 5.0地震为例进行检验性计算,以得出预测时间段。首先去掉7个观测点中的一半(即3个)测点的最小发震开始时间和最大发震结束时间,得出优化后的预测发震时间段,然后再将这三种方法的预测结果进一步优化,得到综合预测的发震时间段是2005年5月至2007年10月,在这个时间段甘肃东南部各观控点的监控范围内有可能发生5级以上地震。
表6 2006年6月21日文县5.0级地震检验性预测结果
Tab.6 Predictive results of MS5.0 earthquake in Wenxian County on June 21,20066 结论与讨论
(1)用从属函数、变差率、趋势速率方法提取的甘肃东南部水氡资料的异常比较明显,对应地震效果较好。
(2)趋势异常本身并不是某一个地震的“专有”前兆,而是多个地震的“共同前兆”。因为,从成因来看,趋势异常不是来自震源的“源兆”,而是“外因性”的“场兆”(张国民等,2001; 范雪芳等,2010)。这一特点决定了依据这类异常提出的发震时间,只是一个“可能的发震时段”,即是说,并不是每一次异常都对应地震。一组趋势异常变化并不只对应一个地震,而是对应一个地震的活跃时段,即趋势异常为多个强震、中强震的共同前兆。
(3)异常最大变幅的大小、异常持续时间与震级的大小没有直接关系。
(4)汶川地震震级大,震源浅,震中距近,但震前甘东南水氡测点出现的异常较少,7个测点中只出现了2项从属函数异常、2项趋势速率异常,3项变差率异常。张培震等(2009)研究发现汶川8.0级地震发生在中国大陆内部的高角度逆冲断层上,震前龙门山断裂带的应变速率很低,形变主要发生在龙门山断裂带以西的川西高原,而同震变形主要发生在龙门山断裂带上。这些现象表明汶川地震孕育和发生的机理与其它特大逆冲型地震可能不同,因此汶川地震前前兆异常少可能与其孕育特征和机理有关。
本文在撰写过程中得到了河北省地震局王吉易研究员和甘肃省地震局杨立明研究员的悉心指导。本文中使用的部分数据由甘肃省地震局预报中心张昱高级工程师提供,在此一并表示感谢!
- 范雪芳,刘耀炜,吴桂娥,等.2010.华北地区水位、水氡中期异常特征研究[J].地震研究,33(2):147-158.
- 范雪芳,王吉易,陆明勇.2009.汶川8.0级地震前典型流体中期前兆异常的初步研究[J].地震,29(1):132-140.
- 范雪芳,王吉易,张淑亮,等.2007.山西及其邻区地下流体强震中期和中短期预测时间方法研究[J].西北地震学报,29(2):177-182.
- 贺仲雄.1983.模糊数学及其应用[M].天津:天津科学技术出版社.
- 王吉易,宋贯一,曹志成,等.2003.地下水诱发的浅层前兆异常及其机理与有关的地震预报问题(3)[J].华北地震科学,21(1):1-10.
- 杨明波,王吉易,刘喜兰,等.2006.北京及邻近地区中等以上地震地下流体异常识别与特征研究[J].地震,26(3):53-63.
- 杨兴悦,王燕,闫万生,等.2006.武山22号井水氡中期异常与地震关系的探讨[J].西北地震学报,28(4):378-380.
- 杨兴悦,王益,王燕,等.2008.武山1号泉流量、O2+Ar、N2、He的异常变化分析[J].西北地震学报,30(1):75-78.
- 叶秀薇.2004.粤闽地区地下流体从属函数异常与地震关系的初步研究[J].防灾减灾工程学报,24(2):195-201.
- 张国民,傅征祥,桂燮泰,等.2001.地震预报引论[M].北京:科学出版社.
- 张立,平建军,苏有锦.2006.云南地区综合地震前兆信息量及其短期映震能力分析[J].地震研究,29(4):325-331.
- 张培震,闻学泽,徐锡伟,等.2009.2008年汶川8.0级特大地震孕育和发生的多单元组合模式[J].科学通报,54(7):944-953.
- 张新基,赵仪全,刘耀炜,等.2005.甘肃省地震监测志[M].甘肃:兰州大学出版社.