甘肅東南部水氡中期異常特征研究＊

楊興悅,王 燕,緱亞江,楊 斐,田文通,李旭生,姚玉霞,范志偉

(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)

利用從屬函數(shù)、變差率以及趨勢速率方法提取甘肅東南部水氡觀測資料的地震前兆異常。結(jié)果表明,1984年1月至 2009年12月,7項(xiàng)水氡資料共出現(xiàn) 47次從屬函數(shù)異常、45次變差率異常和 39次趨勢速率異常。其中,33次從屬函數(shù)異常、35次變差率異常和 30次趨勢速率異常很好地對應(yīng)了甘肅及鄰區(qū)震中距在 33～585 km范圍內(nèi)發(fā)生的 5～8級地震。通過分析異常特征及其與地震的關(guān)系,提出中強(qiáng)震和強(qiáng)震時間預(yù)測的指標(biāo)與方法,得出甘肅東南部水氡從屬函數(shù)計(jì)算的經(jīng)驗(yàn)常數(shù),以及變差率和趨勢速率的異常閾值。

水氡;從屬函數(shù);變差率;趨勢速率;前兆異常;甘肅東南部

0 引言

甘肅處于青藏塊體的東北緣,祁連山地震帶和南北地震帶橫穿境內(nèi),是我國大陸強(qiáng)震活動頻繁的場所之一。該地區(qū)歷史上曾發(fā)生過多次強(qiáng)烈地震,現(xiàn)今地震活動仍很頻繁,地震震源淺,頻度與強(qiáng)度高,危害性強(qiáng)。水氡測項(xiàng)是地震預(yù)測的主要依據(jù)之一。甘肅東南部水氡測點(diǎn)多,具有良好的地震水文地球化學(xué)條件,觀測資料連續(xù)可靠(張新基等,2005;楊興悅等,2008)。為了更好地發(fā)揮其重要作用,對其觀測資料進(jìn)行分析總結(jié)是十分必要的。

由于流體資料或多或少地受儀器、室溫等因素的影響,采用一般的“看圖識字”方法,很難將地震前兆異常與觀測系統(tǒng)本身引起的變化區(qū)別開來(楊明波等,2006;范雪芳等,2007,2009)。為此,在本研究中,筆者使用從屬函數(shù)、變差率及趨勢速率 3種數(shù)學(xué)方法對觀測數(shù)據(jù)進(jìn)行處理,以提取明顯的地震異常(楊興悅等,2006;張立等,2006),進(jìn)而系統(tǒng)地剖析前兆特征,總結(jié)異常變化規(guī)律。

1 前兆資料及地震目錄的選取

甘肅省地下流體觀測臺站主要集中在甘肅東南部(圖 1),因此筆者選用甘肅東南部的流體資料作為研究對象。在資料篩選過程中,采用了王吉易等(2003)的“同一化”研究思路,并遵循 4個原則:①觀測數(shù)據(jù)時間長,時間在 10 a以上;②無人為干擾(包括儀器改變、臺站改造、井(泉)裝置改變、無規(guī)律抽水、渠道放水與停水等;)引起的資料臺階狀變化或?qū)ε_階狀變化因素比較清楚,能進(jìn)行校正,以及無其它變化;③具有明顯的年動態(tài)變化;④測點(diǎn)目前仍在正常觀測。

圖1 甘肅東南部地下流體監(jiān)測臺站分布圖Fig.1 Distribution of underground fluid monitoring stations in southeast Gansu Province

經(jīng)過嚴(yán)格篩選,研究區(qū)共有 7個測點(diǎn)的水氡觀測資料滿足以上 4個條件,分別為武山 1號泉、武都殿溝、武山溫泉 22號井、平?jīng)龈郊S、通渭溫泉、西和蘆河和清水溫泉。本文選用的連續(xù)觀測資料最短為 17 a,最長為 26 a。

考慮本區(qū)地震活動的周期性和強(qiáng)度起伏的特點(diǎn),選取距測點(diǎn) 33～594 km內(nèi)的甘肅及其鄰區(qū) 5級以上地震作為研究對象。在本文中,所研究的5.0、5.2級地震的最大震中距在 200 km左右;5.8～6.2級 4次地震的震中距在 300 km左右;共和 7.0、6.0級地震和汶川 8.0級地震,個別測點(diǎn)的震中距大于 500 km。依據(jù)本文對資料的分析認(rèn)為,甘肅東南部水氡資料對共和 6級及以上地震是有反應(yīng)的。為方便統(tǒng)計(jì),相鄰兩次地震時間間隔小于 1 a的,只統(tǒng)計(jì)發(fā)生在前的地震。所選地震的參數(shù)見表 1。

表1 1987～2009年甘肅及鄰區(qū)發(fā)生的MS≥5.0地震Tab.1 MS≥5.0 Earthquakes in Gansu Province and the adjacent regions from 1987 to 2009

2 異常提取方法

2.1 從屬函數(shù)法

前兆觀測資料的異常變化形態(tài)雖然是多種多樣的,諸如突跳、轉(zhuǎn)折、加速等,但其本質(zhì)都是觀測曲線隨時間變化的斜率變化的量。從屬函數(shù)就是表征一條曲線變化的斜率和跳動程度(賀仲雄,1983;葉秀薇,2004)。通過對資料的分析比較,筆者認(rèn)為,對于甘肅地區(qū),用 13點(diǎn)滑動平均值作為數(shù)據(jù)序列計(jì)算的效果較好,明顯優(yōu)于直接使用月均值提取異常的方法(楊興悅等,2006)。中期異常判定方法見表 2。

表2 甘肅東南部地區(qū)水氡中期異常判定方法Tab.2 Mid-ter m anomalies decision method of Radon in southeast Gansu Province

表2中,ki為觀測值變化斜率;ri為滑動平均值M(t)與時間 t的相關(guān)系數(shù),它反應(yīng)了觀測序列內(nèi)在質(zhì)量的好壞;α為經(jīng)驗(yàn)常數(shù)。在計(jì)算時,13點(diǎn)滑動平均值、從屬函數(shù)值μi均置于計(jì)算時段的右端點(diǎn)。

2.2 變差率法

在中短期異常判定中,觀測值破年變形態(tài)被認(rèn)為是地下流體參量的重要異?，F(xiàn)象。變差率是表征一條曲線的相對變化幅度。以往分析人員只是通過多年曲線特征對比來分析破年變異常,但不能給出 “破年變”的 “程度”。為此,筆者引入了月變差率方法來定量描述這類異常,也就是用定量值來確定本年相對于上年在同一月的變化程度,如果出現(xiàn)明顯的變化差,則表明破年變程度較大,若變差率較小,則表明本年與上年的年動態(tài)相似①劉耀煒,杜建國 .1990.地下水化學(xué)參量中短期異常識別、標(biāo)志、體系及預(yù)報(bào)方法研究 .。設(shè)定了合適的閾值,就可以判定異常的起始時間和持續(xù)時間。變差率方法適合于分析有年變的觀測資料,但在本文中使用的資料連續(xù)時間較長,有的年份資料年變明顯,有的年份無年變,因此在本研究中筆者先對資料進(jìn)行平滑處理后再作分析。其具體做法如下:

(1)月均值的平滑濾波。取水氡等測項(xiàng)的月均值作為基本的數(shù)據(jù)序列。為抑制水氡的年變化和短期異常變化,突出中短期異常,我們采用 13點(diǎn)滑動平均法,求出月均值的滑動平均值。

(2)變差率的計(jì)算。定義變差率為某個月的水氡值與前一年同月的水氡值之差再除以前一年同月的水氡值。該參量用以作為單點(diǎn)的中期和中短期異常的分析。其計(jì)算式為

計(jì)算結(jié)果表明先平滑后再使用變差率方法提取異常的效果好于直接使用變差率。

2.3 趨勢速率法

趨勢變化特征是水化學(xué)參量中期和中短期異常的主要特征之一。由于不同的流體測點(diǎn)所處的構(gòu)造、含水層巖性、井(泉)結(jié)構(gòu)等物理化學(xué)性質(zhì)的差異,趨勢性異常就有上升、下降和波動之分。趨勢上升型中期異常和轉(zhuǎn)折型中短期異常是強(qiáng)震及大震的前兆異常。本文中提取水氡趨勢速率異常時采用以加速為主的上升異常。盡管趨勢速率方法與從屬函數(shù)方法都表示某一時段內(nèi)時間序列斜率大小,都是對線性相關(guān)性綜合特征的描述,但兩種方法得到的結(jié)果不完全相同,在資料分析運(yùn)用中可以作為參考。具體做法是:

(1)13點(diǎn)滑動平均值的計(jì)算。以月均值為數(shù)據(jù)序列做計(jì)算,計(jì)算值放在右端點(diǎn)上。

(2)趨勢速率值(V)的計(jì)算。選擇線性窗 Y=a+bX擬合模型,擬合窗口長度為 13。

(3)異常判定指標(biāo)。若趨勢速率值 Vi≥0.5,則為異常。

3 典型實(shí)例

甘肅東南部部分測點(diǎn)水氡月均值在本文所選震例的震前出現(xiàn)了明顯的異常,有的測點(diǎn)的異常卻不明顯(圖 2),不管明顯與否,都是定性描述,因人而異,帶有經(jīng)驗(yàn)性,且很難在時間尺度上做出預(yù)測。因此,筆者使用從屬函數(shù)、變差率及趨勢速率 3種數(shù)學(xué)方法對觀測數(shù)據(jù)進(jìn)行處理,提取的異常較為明顯,并對異常進(jìn)行定量分析。

武山溫泉 22號井位于武山縣溫泉鄉(xiāng),海拔1 731 m。構(gòu)造上處于西秦嶺北緣深大斷裂帶(西起青海循化,東至陜西寶雞一帶)的溫泉—甘泉深斷裂、聶河斷裂及馬長莊斷裂的交匯部位。含水巖層為印支期似斑狀花崗巖,裂隙發(fā)育。水質(zhì)類型為微礦化度重碳酸鈉型。該井為自流井,井深 154.39 m,水溫 56.3℃,水氡為主要觀測項(xiàng)目之一(楊興悅等,2008)。

圖3為武山 22號井水氡月均值及異常曲線圖。從圖 3可見,在觀測井周邊 5級以上地震發(fā)生前,難以直接從月均值曲線上提取出異常信息,并且對發(fā)震時間也無法作出定量預(yù)測,但利用從屬函數(shù)、變差率以及趨勢速率方法提取的異常信息則較明顯,能定量地給出時間特征參量。以 2006年6月 21日甘肅文縣 MS5.0地震為例,從屬函數(shù)、變差率、趨勢速率異常持續(xù)時間分別為 16、23及16個月,異常結(jié)束至發(fā)震時間分別為 10、5及 8個月。

圖2 甘肅東南部水氡月均值曲線(a)武山 1號泉;(b)武都殿溝;(c)武山 22號井;(d)平?jīng)龈郊S;(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)、從屬函數(shù)(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 從屬函數(shù)中期異常特點(diǎn)

筆者用從屬函數(shù)法分析 1984年1月至 2009年12月甘肅東南部 7個測項(xiàng)的水氡觀測資料,統(tǒng)計(jì)結(jié)果表明,統(tǒng)計(jì)時段內(nèi)共出現(xiàn) 47次異常(圖 4),其中 33次異常很好地對應(yīng)了表 1中所列的甘肅及其鄰區(qū)距測點(diǎn) 33～594 km內(nèi)的 5～8級地震。發(fā)生了地震而沒有出現(xiàn)異常就算漏報(bào),共 18次;出現(xiàn)了異常而沒發(fā)生地震屬虛報(bào),共 14次。中期異常特征列于表 3。通過對甘肅東南部地區(qū)水氡資料的研究,筆者得到水氡從屬函數(shù)中期異常的特點(diǎn):

(1)持續(xù)時間:最短為 2個月,最長為 20個月,平均為 8.6個月。

(2)超前時間Ⅰ:最短為 1個月,最長 33個月,平均為 15.6個月。

(3)超前時間Ⅱ:最短為 -11個月,最長 24個月,平均為 7.0個月。

(4)在異常持續(xù)過程中發(fā)生地震 9次,異常結(jié)束后發(fā)生地震 24次。

(5)異常變幅:最小為 0.52,最大為 0.83,平均為 0.67。

(6)異常對應(yīng)率為 70.2%,漏報(bào)率為 35.3%,虛報(bào)率為 29.8%。

圖4 甘肅東南部水氡從屬函數(shù)異常曲線(a)武山 1號泉;(b)武都殿溝;(c)武山 22號井;(d)平?jīng)龈郊S;(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 變差率中期異常特點(diǎn)

筆者用變差率方法分析了 1984年1月至 2009年12月甘肅東南部 7個測項(xiàng)的水氡觀測資料,統(tǒng)計(jì)結(jié)果表明,共出現(xiàn) 45次異常,2009年出現(xiàn)的兩次異常視為新異常未統(tǒng)計(jì)(圖 5),其中 35次異常很好地對應(yīng) 了甘肅及 鄰區(qū)距測點(diǎn) 33～594 km的 5～8級地震,其中漏 報(bào) 16次,虛報(bào) 10次。中 期異常特征見表 4。通 過對甘肅東 南部水氡 資料的研究,筆者總 結(jié)出水氡 變差率中期 異常的特 點(diǎn):

表3 甘肅東南部水氡從屬函數(shù)中期異常特征統(tǒng)計(jì)表Tab.3 Feature of the mid-ter m anomalies of Suberdinate function of Radon in Southeast Gansu Province

(1)持續(xù)時間:最短為 2個月,最長為 51個月,平均為 16.9個月。

(2)超前時間Ⅰ:最短為 0個月,最長 35個月,平均為 14.4個月。

(3)超前時間Ⅱ:最短為 -39個月,最長 25個月,平均為 0.5個月。

(4)在異常持續(xù)過程中發(fā)生地震 18次,異常結(jié)束后發(fā)生地震 17次。

(5)異常變幅:最小為 1%,最大為 12%,平均為 4.4%。

(6)異常對應(yīng)率為 77.8%,漏報(bào)率為 31.4%,虛報(bào)率為 22.2%。

圖5 甘肅東南部水氡變差率異常曲線(a)武山 1號泉;(b)武都殿溝;(c)武山 22號井;(d)平?jīng)龈郊S;(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)平?jīng)龈郊S;(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 趨勢速率中期異常特點(diǎn)

筆者用趨勢速率方法分析了 1984年1月至2009年12月甘肅東南部 7個測項(xiàng)的水氡觀測資料,發(fā)現(xiàn)該時段內(nèi)共出現(xiàn) 39次異常(2009年出現(xiàn)的一次異常視為新異常未統(tǒng)計(jì))(圖 6),其中 30次異常很好地對應(yīng)了甘肅及鄰區(qū)距測點(diǎn)33～594 km內(nèi)的 5～8級地震,其中漏報(bào) 21次,虛報(bào) 9次。中期異常特征見表 5。通過對這些資料的研究,筆者總結(jié)出甘肅東南部水氡趨勢速率中期異常的特點(diǎn):

表4 甘肅東南部水氡變差率中期異常特征統(tǒng)計(jì)表Tab.4 Features of the mid-term anomalies of variation rate of Radon in southeast Gansu Province

表5 甘肅東南部水氡趨勢速率中期異常特征統(tǒng)計(jì)表Tab.5 Features of the mid-term anomalies of trend rate of Radon in southeast Gansu Province

(1)持續(xù)時間:最短為 1個月,最長為 28個月,平均為 9.7個月。

(2)超前時間Ⅰ:最短為 0個月,最長 33個月,平均為 13.0個月。

(3)超前時間Ⅱ:最短為 -17個月,最長 21個月,平均為 3.3個月。

(4)在異常持續(xù)過程中發(fā)生地震 11次,異常結(jié)束后發(fā)生地震 18次。

(5)異常變幅:最小為 0.5,最大為 1.00,平均為 0.78。

(6)對應(yīng)率為 76.9%,漏報(bào)率為 41.2%,虛報(bào)率為 23.1%。

5 中期時間預(yù)測指標(biāo)

依據(jù)甘肅東南部水氡的多次異常資料,可得到中期異常的超前時間。將異常開始時間至發(fā)震時間的差記為超前時間Ⅰ,將異常結(jié)束時間至發(fā)震時間的差記為超前時間Ⅱ,從表 3、表 4和表 5中可得出各個測點(diǎn)的超前時間Ⅰ和超前時間Ⅱ。由于本文統(tǒng)計(jì)的是一個地區(qū)的資料,觀測點(diǎn)分散,異常出現(xiàn)和結(jié)束的時間差別較大,因此很難得出一個統(tǒng)一的預(yù)測指標(biāo),但可根據(jù)各測點(diǎn)的超前時間Ⅰ的最小值和超前時間Ⅱ的最大值預(yù)測發(fā)震時間段。

當(dāng)從屬函數(shù)等異常達(dá)到預(yù)測指標(biāo)時,可按下式算出發(fā)震時間段:

式中,T1為預(yù)測的發(fā)震開始時間;t1為異常開始時間;ΔT1為超前時間Ⅰ的最小值;T2為預(yù)測的發(fā)震結(jié)束時間;t2為異常結(jié)束時間,ΔT2為超前時間Ⅱ的最大值。

當(dāng)一個測點(diǎn)出現(xiàn)異常后,可用上述方法得出一個發(fā)震時間段,如果多個測點(diǎn)出現(xiàn)異常就可得到多個發(fā)震時間段。將多個發(fā)震時間段進(jìn)行優(yōu)化,通過驗(yàn)證分析,筆者認(rèn)為舍去 1/2以下測點(diǎn)較早的發(fā)震開始時間和較晚的發(fā)震結(jié)束時間,這樣預(yù)測的發(fā)震時間段就大大縮短了。

同樣,變差率、趨勢速率等異常達(dá)到預(yù)測指標(biāo)后可用上述方法優(yōu)化后再對這幾種方法得到的發(fā)震時段進(jìn)一步優(yōu)化。對應(yīng)的發(fā)震時段越靠近預(yù)測發(fā)震時段的中央就越可能是發(fā)生地震的時間。

下面以 2006年6月 21日甘肅文縣 MS5.0地震為例進(jìn)行檢驗(yàn)性計(jì)算,以得出預(yù)測時間段。首先去掉 7個觀測點(diǎn)中的一半(即 3個)測點(diǎn)的最小發(fā)震開始時間和最大發(fā)震結(jié)束時間,得出優(yōu)化后的預(yù)測發(fā)震時間段,然后再將這三種方法的預(yù)測結(jié)果進(jìn)一步優(yōu)化,得到綜合預(yù)測的發(fā)震時間段是 2005年5月至 2007年10月,在這個時間段甘肅東南部各觀控點(diǎn)的監(jiān)控范圍內(nèi)有可能發(fā)生 5級以上地震。

表6 2006年6月 21日文縣 5.0級地震檢驗(yàn)性預(yù)測結(jié)果Tab.6 Predictive results ofMS5.0 earthquake inWenxian County on June 21,2006

6 結(jié)論與討論

(1)用從屬函數(shù)、變差率、趨勢速率方法提取的甘肅東南部水氡資料的異常比較明顯,對應(yīng)地震效果較好。

(2)趨勢異常本身并不是某一個地震的“專有”前兆,而是多個地震的 “共同前兆”。因?yàn)?從成因來看,趨勢異常不是來自震源的 “源兆”,而是 “外因性”的 “場兆”(張國民等,2001;范雪芳等,2010)。這一特點(diǎn)決定了依據(jù)這類異常提出的發(fā)震時間,只是一個 “可能的發(fā)震時段”,即是說,并不是每一次異常都對應(yīng)地震。一組趨勢異常變化并不只對應(yīng)一個地震,而是對應(yīng)一個地震的活躍時段,即趨勢異常為多個強(qiáng)震、中強(qiáng)震的共同前兆。

(3)異常最大變幅的大小、異常持續(xù)時間與震級的大小沒有直接關(guān)系。

(4)汶川地震震級大,震源淺,震中距近,但震前甘東南水氡測點(diǎn)出現(xiàn)的異常較少,7個測點(diǎn)中只出現(xiàn)了 2項(xiàng)從屬函數(shù)異常、2項(xiàng)趨勢速率異常,3項(xiàng)變差率異常。張培震等(2009)研究發(fā)現(xiàn)汶川 8.0級地震發(fā)生在中國大陸內(nèi)部的高角度逆沖斷層上,震前龍門山斷裂帶的應(yīng)變速率很低,形變主要發(fā)生在龍門山斷裂帶以西的川西高原,而同震變形主要發(fā)生在龍門山斷裂帶上。這些現(xiàn)象表明汶川地震孕育和發(fā)生的機(jī)理與其它特大逆沖型地震可能不同,因此汶川地震前前兆異常少可能與其孕育特征和機(jī)理有關(guān)。

本文在撰寫過程中得到了河北省地震局王吉易研究員和甘肅省地震局楊立明研究員的悉心指導(dǎo)。本文中使用的部分?jǐn)?shù)據(jù)由甘肅省地震局預(yù)報(bào)中心張昱高級工程師提供,在此一并表示感謝!

范雪芳,劉耀煒,吳桂娥,等 .2010.華北地區(qū)水位、水氡中期異常特征研究[J].地震研究,33(2):147-158.

范雪芳,王吉易,陸明勇 .2009.汶川 8.0級地震前典型流體中期前兆異常的初步研究[J].地震,29(1):132-140.

范雪芳,王吉易,張淑亮,等 .2007.山西及其鄰區(qū)地下流體強(qiáng)震中期和中短期預(yù)測時間方法研究 [J].西北地震學(xué)報(bào),29(2):177-182.

賀仲雄 .1983.模糊數(shù)學(xué)及其應(yīng)用[M].天津:天津科學(xué)技術(shù)出版社.

王吉易,宋貫一,曹志成,等 .2003.地下水誘發(fā)的淺層前兆異常及其機(jī)理與有關(guān)的地震預(yù)報(bào)問題(3)[J].華北地震科學(xué),21(1):1-10.

楊明波,王吉易,劉喜蘭,等 .2006.北京及鄰近地區(qū)中等以上地震地下流體異常識別與特征研究[J].地震,26(3):53-63.

楊興悅,王燕,閆萬生,等 .2006.武山 22號井水氡中期異常與地震關(guān)系的探討[J].西北地震學(xué)報(bào),28(4):378-380.

楊興悅 ,王益 ,王燕 ,等 .2008.武山 1號泉流量、O2+Ar、N2、He的異常變化分析[J].西北地震學(xué)報(bào),30(1):75-78.

葉秀薇 .2004.粵閩地區(qū)地下流體從屬函數(shù)異常與地震關(guān)系的初步研究[J].防災(zāi)減災(zāi)工程學(xué)報(bào),24(2):195-201.

張國民,傅征祥,桂燮泰,等 .2001.地震預(yù)報(bào)引論[M].北京:科學(xué)出版社.

張立,平建軍,蘇有錦 .2006.云南地區(qū)綜合地震前兆信息量及其短期映震能力分析[J].地震研究,29(4):325-331.

張培震,聞學(xué)澤,徐錫偉,等 .2009.2008年汶川 8.0級特大地震孕育和發(fā)生的多單元組合模式[J].科學(xué)通報(bào),54(7):944-953.

張新基,趙儀全,劉耀煒,等 .2005.甘肅省地震監(jiān)測志[M].甘肅:蘭州大學(xué)出版社 .

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.

中國地震局“三結(jié)合”基金項(xiàng)目(2008042)資助 .