安張輝 詹 艷 陳小斌 姜 峰 高 悅

1)中國(guó)地震局蘭州地震研究所，蘭州 730000 2)中國(guó)地震局地震預(yù)測(cè)研究所，蘭州科技創(chuàng)新基地，蘭州 730000 3)中國(guó)地震局地質(zhì)研究所，地震動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100029

滑動(dòng)自相關(guān)方法在地電阻率觀測(cè)資料分析中的應(yīng)用初探

安張輝1， 2)詹 艷3)陳小斌3)姜 峰3)高 悅1，2)

1)中國(guó)地震局蘭州地震研究所，蘭州 730000 2)中國(guó)地震局地震預(yù)測(cè)研究所，蘭州科技創(chuàng)新基地，蘭州 730000 3)中國(guó)地震局地質(zhì)研究所，地震動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100029

提出了地電阻率滑動(dòng)自相關(guān)數(shù)據(jù)處理方法，利用隨機(jī)時(shí)間序列對(duì)方法進(jìn)行了數(shù)值檢驗(yàn)，驗(yàn)證了方法的有效性。將滑動(dòng)自相關(guān)方法應(yīng)用于3個(gè)地電阻率臺(tái)站實(shí)際觀測(cè)資料，經(jīng)分析研究得到了一些比較有意義的結(jié)果： 1)甘孜臺(tái)和山丹臺(tái)地電阻率滑動(dòng)相關(guān)系數(shù)變化與地震事件有較好的對(duì)應(yīng)關(guān)系，主要表現(xiàn)為，在震前半年到1年時(shí)間范圍內(nèi)，會(huì)出現(xiàn)相關(guān)系數(shù)階段性增大現(xiàn)象，同時(shí)，這2個(gè)臺(tái)站的相關(guān)系數(shù)異常也表現(xiàn)出比較強(qiáng)的各向異性特征； 2)成都臺(tái)地電阻率觀測(cè)雖受干擾較為嚴(yán)重，但相關(guān)系數(shù)異常現(xiàn)象與地震事件仍有比較好的對(duì)應(yīng)關(guān)系，除了說(shuō)明所使用方法的有效性之外，還可能與參加計(jì)算的地震事件震級(jí)比較大，震中距較小，以及發(fā)震時(shí)間比較集中有關(guān)，在該臺(tái)站同樣也存在異常的各向異性特征； 3)通過(guò)對(duì)比不同震級(jí)特征地震的相關(guān)研究結(jié)果得出，當(dāng)選擇的特征地震震級(jí)較小時(shí)，不同地震的震前異常幅度會(huì)有所不同，但當(dāng)震級(jí)較大時(shí)(如MS≥5.0)，其對(duì)研究結(jié)果產(chǎn)生的影響非常有限。此外，還對(duì)地電阻率地震異常的各向異性，以及特征地震的選取問(wèn)題進(jìn)行了簡(jiǎn)單論述。

滑動(dòng)自相關(guān)方法 地電阻率 地震異常 相關(guān)系數(shù) 各向異性

0 引言

將地電阻率觀測(cè)資料應(yīng)用于地震預(yù)測(cè)是中國(guó)的一大特色。地電阻率觀測(cè)始于1966年邢臺(tái)MS7.2地震，其觀測(cè)裝置為對(duì)稱四極觀測(cè)(杜學(xué)彬等，2006，2007； 錢(qián)家棟等，2009)，經(jīng)歷了近50a的觀測(cè)積累和分析研究，得到了與大震和強(qiáng)震突出相關(guān)的多次異?，F(xiàn)象。如： 對(duì)1976年7月28日唐山7.8級(jí)地震進(jìn)行的地電阻率相關(guān)研究(趙玉林等，1978； 錢(qián)復(fù)業(yè)等，1990； 錢(qián)家棟等，1993； Luetal.，1999)，以及對(duì)2008年5月12日汶川MS8.0和2013年4月20日蘆山MS7.0地震地電阻率觀測(cè)開(kāi)展的相關(guān)分析研究(張學(xué)民等，2009； Huang，2011； 錢(qián)家棟等，2013； 朱濤，2013； 杜學(xué)彬等，2015)，都得到了非常好的結(jié)果?？傮w來(lái)看，大地震前的地電阻率變化特征一般表現(xiàn)為： 震前2～3a開(kāi)始出現(xiàn)趨勢(shì)下降異?，F(xiàn)象； 震前地電阻率異常具有階段性，個(gè)別大震前還有短臨和臨震前兆現(xiàn)象； 震前趨勢(shì)異常表現(xiàn)出各向異性(錢(qián)家棟等，1993)等。

地電阻率地震異常特征在1993年就已經(jīng)得出(錢(qián)家棟等，1993)，但2～3a的時(shí)間預(yù)測(cè)范圍，所能起到的減災(zāi)實(shí)效比較有限。到目前為止，針對(duì)地電阻率觀測(cè)資料在實(shí)際地震預(yù)測(cè)中進(jìn)行應(yīng)用的方法研究仍顯乏力。杜學(xué)彬等2001年提出的歸一化月速率方法，推進(jìn)了地電阻率觀測(cè)手段在預(yù)測(cè)預(yù)報(bào)中的實(shí)際應(yīng)用，并開(kāi)展了一些較好的分析工作(葉青等，2005)。

以往開(kāi)展的地電阻率觀測(cè)異常分析，很多是關(guān)于特定地震事件進(jìn)行的回溯性研究。雖然也得到了較好的對(duì)應(yīng)結(jié)論，但這些分析研究只是增加了典型震例的數(shù)量，對(duì)地電阻率觀測(cè)在實(shí)際應(yīng)用工作的推動(dòng)有限。本文借鑒地震電磁異?；瑒?dòng)相關(guān)分析方法(Jiangetal.，2016)的研究思想，提出了滑動(dòng)自相關(guān)地震預(yù)測(cè)方法，并初步應(yīng)用于四川省甘孜、 成都、 甘肅省山丹地電阻率臺(tái)的實(shí)際觀測(cè)資料，進(jìn)行一些比較有意義的探索性分析研究。

1 方法簡(jiǎn)介及數(shù)值檢驗(yàn)

1.1 方法簡(jiǎn)介

依據(jù)經(jīng)驗(yàn)性地震預(yù)測(cè)的特點(diǎn)，假設(shè)在觀測(cè)臺(tái)站一定震中距范圍內(nèi)發(fā)生的地震事件，所引起的地電阻率異?，F(xiàn)象具有一定的相似性。根據(jù)這一假設(shè)條件，為突出地電阻率地震異常的中期和短期特點(diǎn)，文中選取某次地震(下文稱 “特征地震”)前一段時(shí)間(180d)的觀測(cè)資料作為特征時(shí)間序列，通過(guò)滑動(dòng)特征時(shí)間序列與原始時(shí)間序列得到相關(guān)系數(shù)曲線。利用相關(guān)系數(shù)曲線與地震事件序列的對(duì)應(yīng)性，分析研究震前相關(guān)系數(shù)的變化特點(diǎn)，從而對(duì)未來(lái)的地震事件進(jìn)行預(yù)測(cè)。相關(guān)系數(shù)的計(jì)算公式：

(1)

計(jì)算分析過(guò)程中采用地電阻率日均值。與地震事件進(jìn)行相關(guān)性計(jì)算分析時(shí)，應(yīng)考慮1d內(nèi)距離臺(tái)站一定范圍內(nèi)所有地震事件的總體響應(yīng)，借鑒Jiang等(2016)的方法，進(jìn)行等效地震計(jì)算。首先，利用式(2)計(jì)算出1d內(nèi)所應(yīng)考慮地震事件的總釋放能量(Gutenbergetal.，1956)，再由式(3)計(jì)算出相應(yīng)的等效震級(jí)(Jiangetal.，2016)。

(2)

(3)

1.2 數(shù)值檢驗(yàn)

為了驗(yàn)證方法的有效性，利用隨機(jī)時(shí)間序列對(duì)該方法進(jìn)行數(shù)值檢驗(yàn)。首先，需要產(chǎn)生1段隨機(jī)時(shí)間序列，選取其中1段時(shí)間序列作為滑動(dòng)時(shí)間序列； 其次，將滑動(dòng)時(shí)間序列，按不同的時(shí)間間隔疊加在原始時(shí)間序列上，從而形成新的時(shí)間序列，模擬地電阻率觀測(cè)中震前異常具有相似性的特點(diǎn)； 最后，通過(guò)滑動(dòng)所選取的特征時(shí)間序列，計(jì)算與新時(shí)間序列之間的滑動(dòng)相關(guān)系數(shù)，通過(guò)比較不同時(shí)間相關(guān)系數(shù)的大小，達(dá)到驗(yàn)證該方法有效性的目的。具體過(guò)程如下：

圖1 數(shù)值檢驗(yàn)Fig. 1 The demonstration of numerical examination.縱坐標(biāo)為數(shù)據(jù)大?。?a 隨機(jī)時(shí)間序列； b 隨機(jī)時(shí)間序列中的1段數(shù)據(jù)； c 新合成的時(shí)間序列； d 相關(guān)系數(shù)圖

信號(hào)X(t)(圖1a)是隨機(jī)信號(hào)，X1(t)(圖1b)是從X(t)中隨機(jī)選出的1段信號(hào)，長(zhǎng)度為100個(gè)數(shù)據(jù)點(diǎn)。將X1(t)分別疊加在501～600，1,201～1,300，1,501～1,600和1,901～2,000區(qū)間上，從而形成新的時(shí)間序列Y(t)(圖1c)，通過(guò)滑動(dòng)X1(t)與Y(t)做相關(guān)性計(jì)算，從而得到相關(guān)系數(shù)R曲線(圖1d)，由于開(kāi)始滑動(dòng)的位置為1，所以結(jié)束時(shí)的滑動(dòng)位置為2,000，故相關(guān)系數(shù)只有2,000個(gè)。從圖1d中可以看出，在101位置時(shí)相關(guān)系數(shù)為1，表示為完全相關(guān)，符合實(shí)際情況； 在其他4個(gè)插入滑動(dòng)信號(hào)的地方，其相關(guān)系數(shù)都在0.6以上，呈現(xiàn)出很強(qiáng)的相關(guān)性，從而證實(shí)該方法的正確有效性。

2 實(shí)測(cè)數(shù)據(jù)分析

分析研究過(guò)程中使用地電阻率觀測(cè)日均值，計(jì)算過(guò)程中需對(duì)觀測(cè)資料進(jìn)行一些預(yù)處理。首先，計(jì)算觀測(cè)臺(tái)站每天觀測(cè)值的均方根誤差； 其次，對(duì)>3倍均方根誤差的觀測(cè)數(shù)據(jù)進(jìn)行刪除，然后計(jì)算出相應(yīng)的日均值。此外，在滑動(dòng)計(jì)算時(shí)采用地電阻率日均值的一階差分(前差)，可以降低大階躍和其他干擾帶來(lái)的影響，使計(jì)算結(jié)果更為理想。分析過(guò)程中所使用的地震目錄從中國(guó)地震臺(tái)網(wǎng)中心網(wǎng)站下載，目錄截止日期為2014年11月30日。

2.1 甘孜臺(tái)結(jié)果分析

甘孜地電觀測(cè)臺(tái)站位于四川省甘孜縣城郊區(qū)東南方向斯俄鄉(xiāng)境內(nèi)的布絨郎山下，現(xiàn)有2個(gè)觀測(cè)方向(N30°E和 N60°W)，本次分析研究所使用觀測(cè)資料的時(shí)間段為： 2007年9月26日至2015年1月26日。參加滑移的特征地震事件為2013年1月18日發(fā)生的MS5.5地震(30.95°N，99.4°E)，考慮該方法對(duì)地電阻率短期和臨震異常的實(shí)用性，參加滑動(dòng)的特征長(zhǎng)度取180d，參加等效地震計(jì)算的地震事件選取范圍：MS≥6.5地震的震中距范圍≤350km，MS在5.0～6.5地震的震中距范圍≤250km，MS<5.0地震的震中距范圍≤100km。

圖2 是四川省甘孜地電臺(tái)的滑動(dòng)自相關(guān)結(jié)果圖。從圖中可以看出，地震事件與震前的相關(guān)系數(shù)異常有較好的對(duì)應(yīng)關(guān)系，結(jié)果初步顯示，在震前半年到1年時(shí)間范圍內(nèi)，會(huì)出現(xiàn)相關(guān)系數(shù)增大異?，F(xiàn)象，并且這種異?，F(xiàn)象有一定的階段性，并不是持續(xù)發(fā)生； 從這幾次的對(duì)應(yīng)情況來(lái)看，N60°W測(cè)向比N30°E效果好，這可能與地震震中相對(duì)于臺(tái)站的位置有一定的關(guān)系，同時(shí)，也是地電阻率震前異常各向異性的1個(gè)具體體現(xiàn)。由于該臺(tái)所選時(shí)間范圍內(nèi)，地震具有成簇發(fā)生的特點(diǎn)，導(dǎo)致多次震前異常與震后效應(yīng)可能相互疊加，這個(gè)特點(diǎn)可能給震前異常的判斷帶來(lái)一定的困難。

圖2 甘孜臺(tái)滑動(dòng)自相關(guān)結(jié)果Fig. 2 The results of shifting self-correlation at Ganzi station.虛線表示等效地震的發(fā)生時(shí)間位置； 相關(guān)系數(shù)達(dá)到1.0時(shí)，其所對(duì)應(yīng)的地震事件為該臺(tái)站選定的特征地震

2.2 成都臺(tái)結(jié)果分析

成都地電觀測(cè)臺(tái)站位于四川省都江堰市與郫縣交界的走石山，現(xiàn)有2個(gè)觀測(cè)方向(N58°E和 N49°W)，本次分析研究所使用觀測(cè)資料的時(shí)間段為： 2010年10月29日至2015年1月26日。參加滑移的特征地震事件為2011年9月4日四川發(fā)生的MS4.2地震(31.27°N，103.62°E)，參加滑動(dòng)的特征長(zhǎng)度為180d，參加等效地震計(jì)算的地震事件選取范圍：MS≥6.5地震的震中距范圍≤350km，MS在5.0～6.5地震的震中距范圍≤200km，MS<5.0地震的震中距范圍≤100km。

圖3 是成都地電阻率觀測(cè)的處理結(jié)果圖。臺(tái)站受工業(yè)游散電流、 成青快鐵等干擾較為嚴(yán)重，影響了傳統(tǒng)處理方法的正常使用，但滑動(dòng)自相關(guān)方法仍能得到較理想的結(jié)果。從結(jié)果圖中可以看出，相關(guān)系數(shù)異常與地震事件有比較好的對(duì)應(yīng)關(guān)系，利用滑動(dòng)自相關(guān)方法處理得到的異常現(xiàn)象，出現(xiàn)在震前大概半年時(shí)間內(nèi)，其異常時(shí)段的相關(guān)系數(shù)值與甘孜臺(tái)相比要大許多，這種現(xiàn)象可能是由于，在分析研究的時(shí)間段內(nèi)所選地震事件的震級(jí)比較大，而且相對(duì)比較集中發(fā)生的原因； 此外，還有震中距較小的原因，如蘆山MS7.0地震距離成都臺(tái)約100km。該臺(tái)站N49°W測(cè)向?qū)λ拇ㄌJ山MS7.0地震的異常反映比N58°E好，這種現(xiàn)象與王成虎等(2014)得到的該區(qū)域?yàn)槟鏇_斷層應(yīng)力狀態(tài)，最大水平主壓應(yīng)力方向(N44°～64°W)相一致，說(shuō)明地電阻率異常與主壓應(yīng)力之間可能存在一定的關(guān)系。

2.3 山丹臺(tái)結(jié)果分析

山丹地電觀測(cè)臺(tái)站位于甘肅省山丹縣清泉鎮(zhèn)南灣村，現(xiàn)有3個(gè)觀測(cè)方向(SN、 EW和 N45°W)，本次分析研究所使用觀測(cè)資料的時(shí)間段為： 2007年5月11日至2015年1月26日。參加滑移的特征地震事件為2008年3月30日青海發(fā)生的MS5.2地震(37.97°N，101.92°E)，參加滑動(dòng)的特征長(zhǎng)度為180d，參加等效地震計(jì)算的地震事件選取范圍：MS≥6.5地震的震中距范圍≤300km，MS5.0～6.5地震的震中距范圍≤200km，MS<5.0地震的震中距范圍≤100km。

圖4 為甘肅山丹地電阻率觀測(cè)的處理結(jié)果圖。該臺(tái)站EW測(cè)向映震效能較好，SN和 N45°W測(cè)向次之，相關(guān)系數(shù)異?，F(xiàn)象多發(fā)生在震前半年至1年時(shí)間范圍內(nèi)。其中SN測(cè)向表現(xiàn)出的相關(guān)系數(shù)變化形態(tài)較為特殊，其原因可能是由于該測(cè)向?qū)μ卣鞯卣鸬姆从潮容^弱，參加滑動(dòng)的特征時(shí)間序列包含的異常信息非常少所導(dǎo)致。但這種現(xiàn)象在EW和 N45°W 測(cè)向上卻表現(xiàn)得截然相反，這種現(xiàn)象從側(cè)面反映出，山丹臺(tái)存在明顯的地電阻率各向異性現(xiàn)象。EW和 N45°W 測(cè)向在2015年初出現(xiàn)的這種相關(guān)系數(shù)異常現(xiàn)象，或許與該區(qū)域未來(lái)一定時(shí)期內(nèi)某次MS5.0左右的地震相對(duì)應(yīng)。

3 相關(guān)討論

圖5 不同特征地震得到的相關(guān)系數(shù)圖Fig. 5 The correlation coefficient map of different characteristic earthquakes.

在分析研究的過(guò)程中，特征地震(Meq)的選取比較關(guān)鍵，不同震級(jí)的地震事件在震前引起的異常變化不會(huì)完全相同。針對(duì)這種可能情況，以甘肅省山丹地電阻率臺(tái)為例，分析討論了參加滑動(dòng)的特征地震不同時(shí)會(huì)出現(xiàn)何種變化情況。圖5 是山丹臺(tái)地電阻率選取不同特征地震時(shí)得到的相關(guān)系數(shù)圖，圖5a是Meq震級(jí)為3.6的結(jié)果，其震中距約118km； 圖5b是Meq震級(jí)為3.8的結(jié)果，其震中距約54km； 圖5c是Meq震級(jí)為5.3的結(jié)果，其震中距約95km。初步研究發(fā)現(xiàn)，參與滑動(dòng)的特征地震震級(jí)與等效地震震級(jí)相差不大時(shí)，其相關(guān)系數(shù)異常較為明顯，如圖5a和圖5b中震前相關(guān)系數(shù)異常幅度均比圖5c的大； 當(dāng)兩者相差較大時(shí)，仍有相關(guān)系數(shù)異?，F(xiàn)象出現(xiàn)，但異常幅度較小(相關(guān)系數(shù)值較小)。針對(duì)不同的特征地震，所得到的相關(guān)系數(shù)圖存在一定的相似性，反映了假設(shè)條件的正確性和所采用方法的有效性。

從文中3個(gè)臺(tái)站觀測(cè)資料的分析能夠看出，不同測(cè)向?qū)Φ仉娮杪实卣甬惓５姆从掣鞑幌嗤?，除了有震?jí)和震中距的因素以外，可能存在著地電阻率各向異性的原因，關(guān)于這個(gè)現(xiàn)象，杜學(xué)彬等(2007)研究了強(qiáng)地震附近視電阻率各向異性變化的原因，通過(guò)建立真、 視電阻率各向異性變化與裂隙率v、 骨架電阻率r、 飽和裂隙電阻率的本構(gòu)關(guān)系，解釋了視電阻率各向異性的原因是，在強(qiáng)地震孕震晚期階段，在震源區(qū)及附近地殼近地表的較深部分，介質(zhì)微裂隙發(fā)育，其走向沿最大主應(yīng)力方位優(yōu)勢(shì)取向，低阻水填充微裂隙，導(dǎo)電通道連通，引起最大主壓應(yīng)力方向的真電阻率變化最為顯著，導(dǎo)致電阻率的各向異性變化。文中各測(cè)向滑動(dòng)自相關(guān)系數(shù)對(duì)地震事件的反映能力的不同，只是地電阻率各向異性的一個(gè)側(cè)面體現(xiàn)。

4 結(jié)論

本文基于距離地電阻率觀測(cè)臺(tái)站一定范圍內(nèi)發(fā)生的地震，在震前會(huì)出現(xiàn)相關(guān)性較高的異?，F(xiàn)象這一前提假設(shè)條件，進(jìn)行隨機(jī)時(shí)間序列的滑動(dòng)自相關(guān)數(shù)值檢驗(yàn)，其結(jié)果說(shuō)明，通過(guò)相關(guān)系數(shù)的變化能夠?qū)μ囟ǖ臅r(shí)間序列進(jìn)行檢出，從而為后續(xù)特征地震序列的檢測(cè)提供可能。隨后，將滑動(dòng)自相關(guān)分析方法，應(yīng)用于3個(gè)地電阻率臺(tái)站實(shí)際觀測(cè)資料的分析研究，得到了一些比較有意義的結(jié)論：

(1)甘孜臺(tái)地電阻率分析結(jié)果顯示，地震事件與震前相關(guān)系數(shù)異常有較好的對(duì)應(yīng)關(guān)系。主要表現(xiàn)為，在震前半年到1年時(shí)間范圍內(nèi)，會(huì)出現(xiàn)相關(guān)系數(shù)增大的異?，F(xiàn)象，并且這種異?，F(xiàn)象有一定的階段性，并不是持續(xù)發(fā)生。從幾次地震的對(duì)應(yīng)情況來(lái)看，該臺(tái) N60°W 測(cè)向比N30°E向效果好，除了與地震震中相對(duì)于臺(tái)站的位置有一定的關(guān)系外，非?？赡芤彩堑仉娮杪收鹎爱惓８飨虍愋缘木唧w體現(xiàn)。

(2)成都地電阻率臺(tái)站雖受工業(yè)游散電流、 成青快鐵等干擾較為嚴(yán)重，但相關(guān)系數(shù)異常仍與地震事件有比較好的對(duì)應(yīng)關(guān)系。但該臺(tái)站異常時(shí)段的相關(guān)系數(shù)值與甘孜臺(tái)相比要大許多，這種現(xiàn)象可能是由于，在分析研究的時(shí)間段內(nèi)所選地震事件的震級(jí)較大，震中距較小，以及發(fā)震時(shí)間比較集中所導(dǎo)致。此外，該臺(tái)站同樣存在異常的各向異性特征。

(3)與甘孜臺(tái)相類(lèi)似，山丹臺(tái)地電阻率相關(guān)系數(shù)的增大異常也發(fā)生在震前半年至1年時(shí)間范圍內(nèi)。此外，該臺(tái)站也存在明顯的地電阻率異常各向異性現(xiàn)象。

(4)參與計(jì)算的特征地震震級(jí)對(duì)所得到的結(jié)果有一定的影響，當(dāng)選擇的特征地震震級(jí)較小時(shí)，對(duì)震前異常的判斷會(huì)有一定的影響，但當(dāng)震級(jí)較大時(shí)(如MS≥5.0)，其對(duì)研究結(jié)果產(chǎn)生的影響非常有限。MS≥5.0的地震正是我們非常關(guān)心的地震事件，更反映了本文所提方法的實(shí)用性。

從上述研究結(jié)果可看出，滑動(dòng)自相關(guān)方法應(yīng)用于地電阻率觀測(cè)資料，能夠得到比較好的結(jié)果，該方法為提高地電阻率觀測(cè)資料在地震預(yù)測(cè)中的應(yīng)用率提供了新的思路。文中所涉及3個(gè)臺(tái)站都存在異常的各向異性現(xiàn)象，研究出現(xiàn)異常的方位與震中的位置關(guān)系，將是進(jìn)一步工作的重點(diǎn)。未來(lái)可以對(duì)更多觀測(cè)臺(tái)站、 更長(zhǎng)觀測(cè)資料進(jìn)行分析研究，還可以針對(duì)不同觀測(cè)臺(tái)站，從統(tǒng)計(jì)學(xué)的角度分析研究屬于各個(gè)臺(tái)站的特征地震異常和特征異常時(shí)間尺度，逐步提高地電阻率地震異常的判定能力。

致謝 本文在中國(guó)地震局地質(zhì)研究所完成，對(duì)中國(guó)地震局人事教育司和甘肅省地震局人事教育處提供的相關(guān)幫助，以及審稿專家的有益建議表示衷心感謝!

杜學(xué)彬，李寧，葉青，等. 2007. 強(qiáng)地震附近視電阻率各向異性變化的原因 [J]. 地球物理學(xué)報(bào)，50(6): 1802—1810.

DU Xue-bin，LI Ning，YE Qing，etal. 2007. A possible reason for the anisotropic changes in apparent resistivity near the focal region of strong earthquake [J]. Chinese Journal of Geophysics，50(6): 1802—1810(in Chinese).

杜學(xué)彬，劉君，崔騰發(fā)，等. 2015. 2次近距離大震前成都臺(tái)視電阻率重現(xiàn)性、 相似性和各向異性變化 [J]. 地球物理學(xué)報(bào)，58(2)： 576—588.

DU Xue-bin，LIU Jun，CUI Teng-fa，etal. 2015. Repeatability，similarity and anisotropy changes in apparent resistivity of station Chengdu before two great earthquakes in near distance [J]. Chinese Journal of Geophysics，58(2)： 576—588(in Chinese).

杜學(xué)彬，阮?lèi)?ài)國(guó)，范世宏，等. 2001. 強(qiáng)震近震中區(qū)地電阻率變化速率的各向異性 [J]. 地震學(xué)報(bào)，23(3)： 289—297.

DU Xue-bin，RUAN Ai-guo，F(xiàn)AN Shi-hong，etal. 2001. Anisotropy of the variation rate of apparent resistivity near the epicentral region of strong earthquakes [J]. Acta Seismological Sinica，23(3)： 289—297(in Chinese).

杜學(xué)彬，趙家騮，譚大誠(chéng)，等. 2006. 地震臺(tái)站建設(shè)規(guī)范 地電臺(tái)站： 第1部分 地電阻率臺(tái)站 [M]. 北京： 地震出版社.

DU Xue-bin，ZHAO Jia-liu，TAN Da-cheng，etal. 2006. Specification for the Construction of Seismic Station-Geoelectrical Station，Part 1: Geoelectrical Resistivity Observatory [M]. Seismological Press，Beijing(in Chinese).

錢(qián)復(fù)業(yè)，趙玉林，劉婕，等. 1990. 唐山7.8級(jí)地震地電阻率臨震功率譜異常 [J]. 地震，3： 33—39.

QIAN Fu-ye，ZHAO Yu-lin，LIU Jie，etal. 1990. Power spectrum anomaly of earth resistivity immediately before Tangshan earthquakeM=7.8 [J]. Earthquake，3： 33—39(in Chinese).

錢(qián)家棟. 1993. 與大震孕育過(guò)程有關(guān)的地電阻率變化研究 [J]. 中國(guó)地震，9(4)： 341—350.

QIAN Jia-dong. 1993. A study on the changes in geoelectrical resistivity associated with preparatory process of great earthquake in China [J]. Earthquake Research in China，9(4)： 341—350(in Chinese).

錢(qián)家棟，杜學(xué)彬，蔡晉安，等. 2009. 地震地電觀測(cè)方法 地電阻率： 第1部分 單極距觀測(cè) [M]. 北京： 地震出版社.

QIAN Jia-dong，DU Xue-bin，CAI Jin-an.etal. 2009. The Method of Earthquake-related Geoelectrical Monitoring-Geoelectrical Resistivity Observation-Part 1: The Single Separation Configuration [M]. Seismological Press，Beijing(in Chinese).

錢(qián)家棟，馬欽忠，李劭秾. 2013. 汶川MS8.0地震前成都臺(tái)NE測(cè)線地電阻率異常的進(jìn)一步分析 [J]. 地震學(xué)報(bào)，35(1)： 4—17.

QIAN Jia-dong，MA Qin-zhong，LI Shao-nong. 2013. Further study on the anomalies in apparent resistivity in the NE configuration at Chengdu station associated with WenchuanMS8.0 earthquake [J]. Acta Seismologica Sinica，35(1)： 4—17(in Chinese).

王成虎，宋成科，郭啟良，等. 2014. 利用原地應(yīng)力實(shí)測(cè)資料分析蘆山地震震前淺部地殼應(yīng)力積累 [J]. 地球物理學(xué)報(bào)，57(1): 102—114.

WANG Cheng-hu，SONG Cheng-ke，GUO Qi-liang，etal. 2014. Stress build-up in the shallow crust before the Lushan earthquake based on the in-situ stress measurements [J]. Chinese Journal of Geophysics，57(1): 102—114(in Chinese).

葉青，杜學(xué)彬，陳軍營(yíng)，等. 2005. 2003年大姚和民樂(lè)-山丹地震1年尺度預(yù)測(cè) [J]. 地震研究，28(3)： 226—230.

YE Qing，DU Xue-bin，CHEN Jun-ying，etal. 2005. One year prediction for the Dayao and Minle-Shandan earthquakes in 2003 [J]. Journal of Seismological Research，28(3)： 226—230(in Chinese).

張學(xué)民，李美，關(guān)華平. 2009. 汶川8.0級(jí)地震前的地電阻率異常分析 [J]. 地震，29(1)： 109—115.

ZHANG Xue-min，LI Mei，GUAN Hua-ping. 2009. Anomaly analysis of earth resistivity observations before the Wenchuan earthquake [J]. Earthquake，29(1)： 109—115(in Chinese).

趙玉林，錢(qián)復(fù)業(yè). 1978. 唐山7.8級(jí)強(qiáng)震前震中周?chē)巫冸娮杪实南陆诞惓?[J]. 地球物理學(xué)報(bào)，21(3)： 181—190.

ZHAO Yu-lin，QIAN Fu-ye. 1978. Electrical resistivity anomaly observed in and around the epicentral area prior to the Tangshan earthquake of 1976 [J]. Chinese Journal of Geophysics，21(3)： 181—190(in Chinese).

朱濤. 2013. 汶川MS8.0地震前區(qū)域性視電阻率異常初步研究 [J]. 地震學(xué)報(bào)，35(1)： 18—25.

ZHU Tao. 2013. Preliminary study on regional geo-resistivity anomaly before the WenchuanMS8.0 earthquake [J]. Acta Seismologica Sinica，35(1)： 18—25(in Chinese).

Gutenberg B，Richter C F. 1956. Magnitude and energy of earthquakes [J]. Annals of Geophysics，9: 1—15.

Huang Q H. 2011. Retrospective investigation of geophysical data possible associated with theMS8.0 Wenchuan earthquake in Sichuan，China [J]. Journal of Asian Earth Sciences，41(45): 421— 427.

Jiang F，Chen X B，Zhan Y，etal. 2016. Shifting correlation between earthquakes and electromagnetic signals: A case study of the 2013 Minxian-ZhangxianML6.5(MW6.1)earthquake in Gansu，China [J].Pure and Applied Geophysics，173(1): 269—284.

Lu J，Qian F Y，Zhao Y L. 1999. Sensitivity analysis of the Schlumberger monitoring array: Application to changes of resistivity prior to the 1976 earthquake in Tangshan，China [J]. Tectonophy, 307: 397— 405.

SHIFTING SELF-CORRELATION METHOD INITIALY USED IN THE APPARENT RESISTIVITY OBSERVATION DATA

AN Zhang-hui1，2)ZHAN Yan3)CHEN Xiao-bin3)JIANG Feng3)GAO Yue1，2)

1)LanzhouInstituteofSeismology，ChinaEarthquakeAdministration，Lanzhou730000，China2)LanzhouBaseofInstituteofEarthquakePrediction，ChinaEarthquakeAdministration，Lanzhou730000，China3)StateKeyLaboratoryofEarthquakeDynamics，InstituteofGeology，ChinaEarthquakeAdministration，Beijing100029，China

In this paper，we propose a method of seismic prediction using the geo-electric resistivity shifting self-correlation(SSC)，and a numerical test is carried out using random time series analysis to verify the validity of the method. The SSC method is applied to the actual observation data of three geo-electric resistivity stations，and results are obtained as follows: (1)SSC coefficient changes in Ganzi and Shandan stations have good correspondence to earthquake，which is represented mainly by the phased increase of correlation coefficient appearing six months to a year before the earthquake. At the same time，the correlation coefficient anomalies of the two stations also exhibit strong anisotropy. (2)Although Chengdu geo-electric resistivity station had suffered serious disturbance，the correlation coefficient anomaly also has a good correspondence with earthquake. In addition to the validity of the SSC method，it may also be attributed to the magnitude of the earthquake event，the smaller distance of epicenter，and the time of the earthquake. Anisotropy also exists in the anomaly at Chengdu station. (3)By comparing the characteristics of different magnitudes of earthquakes，the results are obtained that，when the magnitude of the selected characteristic earthquake is relatively small，the amplitude of the anomaly before earthquake is different，but when the magnitude is larger，for exampleMS≥5.0，the impact on the results of this study is very limited. In addition，we briefly discussed the anisotropy of seismic geoelectrical resistivity anomalies and the selection of the characteristic earthquake.

shifting self-correlation method，apparent resistivity，earthquake anomaly，correlation coefficient，anisotropy

10.3969/j.issn.0253- 4967.2016.04.017

2015-09-25收稿，2016-05-21改回。

甘肅省科技計(jì)劃項(xiàng)目(145RJZA030)、 中國(guó)地震局系統(tǒng)交流訪問(wèn)學(xué)者學(xué)習(xí)經(jīng)費(fèi)與中國(guó)地震局地震預(yù)測(cè)研究所基本科研業(yè)務(wù)專項(xiàng)(2012IESLZ04)共同資助。

P631.3+22

A

0253-4967(2016)04-1019-11

安張輝，男，1978年生，2004年畢業(yè)于中國(guó)地震局蘭州地震研究所，獲固體地球物理學(xué)碩士學(xué)位，副研究員，主要從事地震電磁學(xué)研究和全國(guó)地電臺(tái)網(wǎng)技術(shù)管理工作，電話： 0931-8277101，E-mail: anzhanghui5@gsdzj.gov.cn。