杨兴悦,王 燕,缑亚江,杨 斐,田文通,李旭生,姚玉霞,范志伟

(1.中国地震局兰州地震研究所,兰州 730000;2.甘肃省地震局武山地震台,甘肃天水 741308;3.甘肃省地震局,兰州 730000;4.河北省地震局,石家庄 050021)

甘肃东南部水氡中期异常特征研究＊

杨兴悦1,2,王 燕3,缑亚江3,杨 斐3,田文通3,李旭生3,姚玉霞3,范志伟4

(1.中国地震局兰州地震研究所,兰州 730000;2.甘肃省地震局武山地震台,甘肃天水 741308;3.甘肃省地震局,兰州 730000;4.河北省地震局,石家庄 050021)

利用从属函数、变差率以及趋势速率方法提取甘肃东南部水氡观测资料的地震前兆异常。结果表明,1984年1月至 2009年12月,7项水氡资料共出现 47次从属函数异常、45次变差率异常和 39次趋势速率异常。其中,33次从属函数异常、35次变差率异常和 30次趋势速率异常很好地对应了甘肃及邻区震中距在 33～585 km范围内发生的 5～8级地震。通过分析异常特征及其与地震的关系,提出中强震和强震时间预测的指标与方法,得出甘肃东南部水氡从属函数计算的经验常数,以及变差率和趋势速率的异常阈值。

水氡;从属函数;变差率;趋势速率;前兆异常;甘肃东南部

0 引言

甘肃处于青藏块体的东北缘,祁连山地震带和南北地震带横穿境内,是我国大陆强震活动频繁的场所之一。该地区历史上曾发生过多次强烈地震,现今地震活动仍很频繁,地震震源浅,频度与强度高,危害性强。水氡测项是地震预测的主要依据之一。甘肃东南部水氡测点多,具有良好的地震水文地球化学条件,观测资料连续可靠(张新基等,2005;杨兴悦等,2008)。为了更好地发挥其重要作用,对其观测资料进行分析总结是十分必要的。

由于流体资料或多或少地受仪器、室温等因素的影响,采用一般的“看图识字”方法,很难将地震前兆异常与观测系统本身引起的变化区别开来(杨明波等,2006;范雪芳等,2007,2009)。为此,在本研究中,笔者使用从属函数、变差率及趋势速率 3种数学方法对观测数据进行处理,以提取明显的地震异常(杨兴悦等,2006;张立等,2006),进而系统地剖析前兆特征,总结异常变化规律。

1 前兆资料及地震目录的选取

甘肃省地下流体观测台站主要集中在甘肃东南部(图 1),因此笔者选用甘肃东南部的流体资料作为研究对象。在资料筛选过程中,采用了王吉易等(2003)的“同一化”研究思路,并遵循 4个原则:①观测数据时间长,时间在 10 a以上;②无人为干扰(包括仪器改变、台站改造、井(泉)装置改变、无规律抽水、渠道放水与停水等;)引起的资料台阶状变化或对台阶状变化因素比较清楚,能进行校正,以及无其它变化;③具有明显的年动态变化;④测点目前仍在正常观测。

图1 甘肃东南部地下流体监测台站分布图Fig.1 Distribution of underground fluid monitoring stations in southeast Gansu Province

经过严格筛选,研究区共有 7个测点的水氡观测资料满足以上 4个条件,分别为武山 1号泉、武都殿沟、武山温泉 22号井、平凉附件厂、通渭温泉、西和芦河和清水温泉。本文选用的连续观测资料最短为 17 a,最长为 26 a。

考虑本区地震活动的周期性和强度起伏的特点,选取距测点 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 2009

2 异常提取方法

2.1 从属函数法

前兆观测资料的异常变化形态虽然是多种多样的,诸如突跳、转折、加速等,但其本质都是观测曲线随时间变化的斜率变化的量。从属函数就是表征一条曲线变化的斜率和跳动程度(贺仲雄,1983;叶秀薇,2004)。通过对资料的分析比较,笔者认为,对于甘肃地区,用 13点滑动平均值作为数据序列计算的效果较好,明显优于直接使用月均值提取异常的方法(杨兴悦等,2006)。中期异常判定方法见表 2。

表2 甘肃东南部地区水氡中期异常判定方法Tab.2 Mid-ter m anomalies decision method of Radon in southeast Gansu Province

表2中,ki为观测值变化斜率;ri为滑动平均值M(t)与时间 t的相关系数,它反应了观测序列内在质量的好坏;α为经验常数。在计算时,13点滑动平均值、从属函数值μi均置于计算时段的右端点。

2.2 变差率法

在中短期异常判定中,观测值破年变形态被认为是地下流体参量的重要异常现象。变差率是表征一条曲线的相对变化幅度。以往分析人员只是通过多年曲线特征对比来分析破年变异常,但不能给出 “破年变”的 “程度”。为此,笔者引入了月变差率方法来定量描述这类异常,也就是用定量值来确定本年相对于上年在同一月的变化程度,如果出现明显的变化差,则表明破年变程度较大,若变差率较小,则表明本年与上年的年动态相似①刘耀炜,杜建国 .1990.地下水化学参量中短期异常识别、标志、体系及预报方法研究 .。设定了合适的阈值,就可以判定异常的起始时间和持续时间。变差率方法适合于分析有年变的观测资料,但在本文中使用的资料连续时间较长,有的年份资料年变明显,有的年份无年变,因此在本研究中笔者先对资料进行平滑处理后再作分析。其具体做法如下:

(1)月均值的平滑滤波。取水氡等测项的月均值作为基本的数据序列。为抑制水氡的年变化和短期异常变化,突出中短期异常,我们采用 13点滑动平均法,求出月均值的滑动平均值。

(2)变差率的计算。定义变差率为某个月的水氡值与前一年同月的水氡值之差再除以前一年同月的水氡值。该参量用以作为单点的中期和中短期异常的分析。其计算式为

计算结果表明先平滑后再使用变差率方法提取异常的效果好于直接使用变差率。

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 ofmonth mean values of Radon in southeast Gansu Province(a)No.1 spring ofWushan station;(b)Diangou observing spot ofWudu 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

图3 武山 22号井水氡月均值(a)、从属函数(b)、变差率(c)及趋势速率值(d)异常曲线Fig.3 Abnormal curves ofmonth-mean values(a),subordinate function(b),variation rate(c)and trend rate(d)of Radon in No.22 well atWushan station

4 前兆异常特征

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次。

(5)异常变幅:最小为 0.52,最大为 0.83,平均为 0.67。

(6)异常对应率为 70.2%,漏报率为 35.3%,虚报率为 29.8%。

图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 ofWushan station;(b)Diangou observing spot ofWudu 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.2 变差率中期异常特点

笔者用变差率方法分析了 1984年1月至 2009年12月甘肃东南部 7个测项的水氡观测资料,统计结果表明,共出现 45次异常,2009年出现的两次异常视为新异常未统计(图 5),其中 35次异常很好地对应 了甘肃及 邻区距测点 33～594 km的 5～8级地震,其中漏 报 16次,虚报 10次。中 期异常特征见表 4。通 过对甘肃东 南部水氡 资料的研究,笔者总 结出水氡 变差率中期 异常的特 点:

表3 甘肃东南部水氡从属函数中期异常特征统计表Tab.3 Feature of the mid-ter m anomalies of Suberdinate function of Radon in Southeast Gansu Province

(1)持续时间:最短为 2个月,最长为 51个月,平均为 16.9个月。

(2)超前时间Ⅰ:最短为 0个月,最长 35个月,平均为 14.4个月。

(3)超前时间Ⅱ:最短为 -39个月,最长 25个月,平均为 0.5个月。

(4)在异常持续过程中发生地震 18次,异常结束后发生地震 17次。

(5)异常变幅:最小为 1%,最大为 12%,平均为 4.4%。

(6)异常对应率为 77.8%,漏报率为 31.4%,虚报率为 22.2%。

图5 甘肃东南部水氡变差率异常曲线(a)武山 1号泉;(b)武都殿沟;(c)武山 22号井;(d)平凉附件厂;(e)通渭温泉;(f)西和芦河;(g)清水温泉Fig.5 Abno rmal curves of variation rate of Radon in southeast area of Gansu Province(a)No.1 spring ofWushan station;(b)Diangou observing spot ofWudu 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 甘肃东南部水氡趋势速率异常曲线(a)武山 1号泉;(b)武都殿沟;(c)武山 22号井;(d)平凉附件厂;(e)通渭温泉;(f)西和芦河;(g)清水温泉Fig.6 Abno rmal curves of trend rate of Radon in southeast area of Gansu Province(a)No.1 spring ofWushan station;(b)Diangou observing spot ofWudu 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.3 趋势速率中期异常特点

笔者用趋势速率方法分析了 1984年1月至2009年12月甘肃东南部 7个测项的水氡观测资料,发现该时段内共出现 39次异常(2009年出现的一次异常视为新异常未统计)(图 6),其中 30次异常很好地对应了甘肃及邻区距测点33～594 km内的 5～8级地震,其中漏报 21次,虚报 9次。中期异常特征见表 5。通过对这些资料的研究,笔者总结出甘肃东南部水氡趋势速率中期异常的特点:

表4 甘肃东南部水氡变差率中期异常特征统计表Tab.4 Features of the mid-term anomalies of variation rate of Radon in southeast Gansu Province

表5 甘肃东南部水氡趋势速率中期异常特征统计表Tab.5 Features of the mid-term anomalies of trend rate of Radon in southeast Gansu Province

(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为超前时间Ⅰ的最小值;T2为预测的发震结束时间;t2为异常结束时间,ΔT2为超前时间Ⅱ的最大值。

当一个测点出现异常后,可用上述方法得出一个发震时间段,如果多个测点出现异常就可得到多个发震时间段。将多个发震时间段进行优化,通过验证分析,笔者认为舍去 1/2以下测点较早的发震开始时间和较晚的发震结束时间,这样预测的发震时间段就大大缩短了。

同样,变差率、趋势速率等异常达到预测指标后可用上述方法优化后再对这几种方法得到的发震时段进一步优化。对应的发震时段越靠近预测发震时段的中央就越可能是发生地震的时间。

下面以 2006年6月 21日甘肃文县 MS5.0地震为例进行检验性计算,以得出预测时间段。首先去掉 7个观测点中的一半(即 3个)测点的最小发震开始时间和最大发震结束时间,得出优化后的预测发震时间段,然后再将这三种方法的预测结果进一步优化,得到综合预测的发震时间段是 2005年5月至 2007年10月,在这个时间段甘肃东南部各观控点的监控范围内有可能发生 5级以上地震。

表6 2006年6月 21日文县 5.0级地震检验性预测结果Tab.6 Predictive results ofMS5.0 earthquake inWenxian County on June 21,2006

6 结论与讨论

(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].甘肃:兰州大学出版社 .

Study on M id-term Anomalies of Radon in Southeast Gansu Province

YANG X ing-yue1,2,WANG Yan3,GOU Ya-jiang3,YANG Fei3,TIAN W en-tong3,L I Xu-sheng3,YAO Yu-xia3,FAN Zhi-w ei4
(1.Lanzhou Institute of Seismology,CEA,Lanzhou730000,Gansu,China)(2.Wushan Seism ic Station,Earthquake Adm inistration of Gansu Province,Tianshui741308,Gansu,China)(3.Earthquake Adm inistration of Gansu Province,Lanzhou730000,Gansu,China)(4.Earthquake Adm inistration of Hebei Province,Shijiazhuang050021,Hebei,China)

Using the m ethods of subordinate function,variation rate and trend rate,w e extract the earthquake precursory anom alies from observational data of Radon in southeast Gansu Province.The result indicates that7item s of radon data include47subordinate-function anom alies,45variation-rate anom alies and39trend-rate anom alies from January1984to D ecem ber2009.Am ong them,33subordinate-function anom alies,35variation-rate anom alies and30trend-rate anom alies corresponded w ith the earthquakesMS5-8that ever hit Gansu Province and the adjacent region w ithin33to585km.B y the analysis of the features of these anom aies and their relation to the earthquakes,the m ethod of tim e prediction for the m oderate-strong and strong earthquake w ere proposed,and the empirical constant,variation rate and threshold value of trend-rate anom aly are obtained from Radon in southeast Gansu Province.

Radon; subordinate function; variation rate; trend rate; precursory anom aly; Southeast Gansu Province

P315.75

A

1000-0666(2011)01-0008-11

2010-03-02.

中国地震局“三结合”基金项目(2008042)资助 .