摘" 要：結(jié)合寶雞地震臺(tái)的地質(zhì)構(gòu)造和觀測(cè)特點(diǎn)，對(duì)寶雞地震臺(tái)2015—2018年地電阻率觀測(cè)數(shù)據(jù)進(jìn)行系統(tǒng)分析。作者利用長(zhǎng)趨勢(shì)變化分析方法、傅氏滑動(dòng)去年周期方法、歸一化月速率方法、差分方法對(duì)寶雞地震臺(tái)地電阻率數(shù)據(jù)進(jìn)行處理，探討分析九寨溝地震前寶雞地震臺(tái)地電阻率可能存在的異常，并對(duì)寶雞地震臺(tái)地電阻率的觀測(cè)效能及存在的問(wèn)題進(jìn)行分析。

關(guān)鍵詞：寶雞地震臺(tái)；地電阻率；九寨溝Ms7.0級(jí)地震；數(shù)據(jù)處理；異常分析

中圖分類號(hào)：P631.3+22" " " 文獻(xiàn)標(biāo)志碼：A" " " " " 文章編號(hào)：2095-2945（2025）03-0165-04

Abstract： Based on the geological structure and observation characteristics of Baoji Seismic Station， the earth resistivity observation data of Baoji Seismic Station from 2015 to 2018 were systematically analyzed. The author used the long-trend change analysis method， Fourier sliding last year period method， normalized monthly rate method， and difference method to process the resistivity data of the Baoji Seismic Station， discussed and analyzed the possible anomalies of the resistivity of the Baoji Seismic Station before the Jiuzhaigou earthquake， and analyzed the observation efficiency and existing problems of the resistivity of the Baoji Seismic Station.

Keywords： Baoji Seismic Station; earth resistivity; Jiuzhaigou Ms7.0 earthquake; data processing; anomaly analysis

地震在孕育過(guò)程中，震源區(qū)應(yīng)力積累以及巖體內(nèi)的微裂膨脹，均會(huì)引起該區(qū)域地下介質(zhì)電性結(jié)構(gòu)改變。這種改變可以通過(guò)地表的地電阻率觀測(cè)方法觀測(cè)到，地電阻率法即是觀測(cè)地球淺部介質(zhì)的電學(xué)性質(zhì)隨時(shí)間變化來(lái)進(jìn)行預(yù)報(bào)地震的方法[1]。經(jīng)過(guò)多年不斷努力的發(fā)展，地電阻率法用于地震監(jiān)測(cè)、預(yù)報(bào)方面已經(jīng)取得了豐碩的研究成果[2-5]。

筆者根據(jù)寶雞地震臺(tái)（簡(jiǎn)稱“寶雞臺(tái)”）地質(zhì)構(gòu)造及觀測(cè)特點(diǎn)，對(duì)寶雞臺(tái)地電阻率觀測(cè)數(shù)據(jù)進(jìn)行了處理，并對(duì)地電阻率觀測(cè)數(shù)據(jù)的曲線形態(tài)、變化特征、變化幅度等進(jìn)行了分析。為探討分析寶雞臺(tái)地電阻率觀測(cè)的映震效能及提取與地震孕育有關(guān)的地電阻率異常具有非常重要的現(xiàn)實(shí)意義。

1" 寶雞地震臺(tái)基本情況及觀測(cè)簡(jiǎn)述

寶雞地震臺(tái)（磻溪）位于寶雞市陳倉(cāng)區(qū)磻溪鎮(zhèn)新莊村，海拔高度665 m，地殼厚度為45 km。構(gòu)造上處于控制渭河斷陷盆地的秦嶺北緣斷裂以北、渭河斷裂的南側(cè)，渭河?xùn)|西向斷裂與虢鎮(zhèn)-固關(guān)北西向斷裂的交匯部位附近（如圖1所示），地電阻率觀測(cè)布極區(qū)在秦嶺北麓一級(jí)階地上，測(cè)區(qū)內(nèi)地勢(shì)較為平坦，各測(cè)道介質(zhì)電性參數(shù)大體相近，電測(cè)深曲線類型為KQQ型[6]。

寶雞地震臺(tái)深井地電阻率觀測(cè)系統(tǒng)于2013年3月建成并投入觀測(cè)，屬于深井觀測(cè)，布設(shè)3個(gè)測(cè)道，布極方式如圖2所示。電極采用自制鉛塊澆鑄的圓柱筒型電極，外線路全部采用穿管地埋方式鋪設(shè)。

2" 寶雞臺(tái)地電阻率映震效能分析

通過(guò)對(duì)寶雞臺(tái)地電阻率觀測(cè)數(shù)據(jù)的分析可知，寶雞臺(tái)平均連續(xù)率在99%以上，完整率在98%以上，平均觀測(cè)精度在0.03%左右，觀測(cè)質(zhì)量較高。臺(tái)站布極屬于深井觀測(cè)干擾較小，觀測(cè)資料真實(shí)可靠，能夠真實(shí)記錄到地下介質(zhì)電阻率的變化，可能包含著地震前后區(qū)域地下介質(zhì)電性結(jié)構(gòu)的改變，對(duì)提取與地震孕育有關(guān)的地電阻率異常具有非常現(xiàn)實(shí)的意義。

受限于臺(tái)站周邊的中、強(qiáng)地震較少，無(wú)較多震例對(duì)應(yīng)及分析，將寶雞臺(tái)與鄰區(qū)的地電監(jiān)測(cè)資料聯(lián)合應(yīng)用，必能為陜西乃至全國(guó)的地震預(yù)報(bào)提供有力的數(shù)據(jù)支撐。

3" 數(shù)據(jù)處理及異常分析

3.1" 長(zhǎng)趨勢(shì)變化分析方法

地電阻率的原始曲線上識(shí)別曲線形態(tài)、年變化幅度及相位變化等是探索分析地震前兆異常的常用方法之一，原始數(shù)據(jù)曲線雖然包含著各種干擾變化，但同時(shí)也是最為真實(shí)的。寶雞臺(tái)地電阻率月均值曲線如圖3所示，NS向和NE向在2016年5月同時(shí)發(fā)生趨勢(shì)轉(zhuǎn)折，均由上升轉(zhuǎn)為下降趨勢(shì)，變化明顯；EW向在2016年7月由快速下降轉(zhuǎn)為平穩(wěn)變化，至九寨溝Ms7.0級(jí)地震后未發(fā)生明顯改變。寶雞臺(tái)距九寨溝Ms7.0級(jí)地震的震中距為356 km。分析認(rèn)為，寶雞臺(tái)地電阻率原始數(shù)據(jù)曲線變化可能存在一定的九寨溝Ms7.0級(jí)地震的前兆異常指征。

3.2" 傅氏滑動(dòng)去年周期方法

傅氏滑動(dòng)方法是消除地電阻率年變的數(shù)據(jù)處理方法，該方法利用三角級(jí)數(shù)擬合，在年變化周期T上進(jìn)行時(shí)間序列的遞推，以建立新的消除年變化后的序列函數(shù)[7]。在排除較大干擾后，對(duì)寶雞臺(tái)2015—2018年地電阻率月均值運(yùn)用傅氏方法進(jìn)行分析：2016年12月開(kāi)始，3個(gè)測(cè)道同步出現(xiàn)趨勢(shì)轉(zhuǎn)折異常（如圖4所示），NS向和NE向均在2016年12月出現(xiàn)極大值后開(kāi)始破年變由原來(lái)的趨勢(shì)上升轉(zhuǎn)變?yōu)橼厔?shì)下降，下降趨勢(shì)維持8個(gè)月后發(fā)生九寨溝Ms7.0級(jí)地震；EW向由原來(lái)速率較快的趨勢(shì)下降在2016年12月后轉(zhuǎn)變?yōu)槠椒€(wěn)變化。選擇2015年地電阻率均值為正常背景值，計(jì)算得出NS、EW、NE 3個(gè)測(cè)道的變化幅度分別為2.1%、2.2%、1.5%，且均值超過(guò)了正常背景時(shí)段標(biāo)準(zhǔn)偏差的3倍。分析認(rèn)為用傅氏方法分析寶雞臺(tái)地電阻率數(shù)據(jù)在曲線形態(tài)、變化幅度和時(shí)間尺度等方面可能存在一定的異常指征。

3.3" 歸一化月速率方法

歸一化月速率，是一種無(wú)量綱的地電學(xué)處理方法，其以一定的步長(zhǎng)計(jì)算月均值曲線的斜率，并進(jìn)行歸一化處理[8]。筆者排除年變化和非年變周期成分后對(duì)寶雞臺(tái)地電阻率進(jìn)行歸一化月速率處理（如圖5所示），滑動(dòng)步長(zhǎng)取8個(gè)月，異常判據(jù)為超過(guò)±2.4視為異常。寶雞臺(tái)NS向在2016年2月—7月出現(xiàn)速率正異常，速率異常峰值為3.35，13個(gè)月后發(fā)生九寨溝Ms7.0級(jí)地震；寶雞臺(tái)NE向與NS向的速率曲線具有很大相似性，NE向在2016年4月—6月首次出現(xiàn)速率正異常，速率異常峰值為2.7，異常結(jié)束14個(gè)月發(fā)生九寨溝Ms7.0級(jí)地震；NS向和NE向均在2017年10—11月速率曲線發(fā)生了負(fù)異常，異常峰值分別為-2.69、-2.72，分析認(rèn)為第二次速率異常與測(cè)區(qū)的較大降雨導(dǎo)致地電阻率測(cè)值快速下降有關(guān)。EW向未出現(xiàn)明顯速率異常現(xiàn)象。

3.4" 差分方法

差分頻次[9]和差分能量[9]方法也是地電學(xué)中經(jīng)常使用的方法，其優(yōu)點(diǎn)是能夠抑制數(shù)據(jù)曲線中長(zhǎng)周期的變化成分而突出擾動(dòng)變化，用于分析地震短臨異常變化的效果較好。差分能量方法可以兼顧差分值大小和差分頻次這2個(gè)因素[9]。

寶雞臺(tái)二階差分頻次和二階差分能量如圖6所示。NS向二階差分頻次在九寨溝Ms7.0級(jí)地震前出現(xiàn)高峰，九寨溝Ms7.0級(jí)地震前后，NS向的二階差分能量值也較高；EW向二階差分頻次和能量在九寨溝Ms7.0級(jí)地震發(fā)生幾天后出現(xiàn)差分頻次和能量值均增高現(xiàn)象；NE向在九寨溝Ms7.0級(jí)地震前二階差分頻次出現(xiàn)最高峰，能量值則無(wú)明顯異?，F(xiàn)象。分析認(rèn)為，九寨溝Ms7.0級(jí)地震前，寶雞臺(tái)地電阻率的二階差分頻次和能量的升高可能反映了應(yīng)力作用下地下介質(zhì)中激發(fā)了低頻電磁場(chǎng)擾動(dòng)，地震后隨著應(yīng)力狀態(tài)的調(diào)整，差分頻次和差分能量值依然較高。

4" 結(jié)束語(yǔ)

1）寶雞臺(tái)地電阻率觀測(cè)數(shù)據(jù)質(zhì)量較好，精度高，能夠較為真實(shí)地記錄到地下介質(zhì)電阻率的變化，可能包含著地震前后區(qū)域地下介質(zhì)電性結(jié)構(gòu)的改變，對(duì)提取與地震孕育有關(guān)的地電阻率異常具有非?，F(xiàn)實(shí)的意義。受限于陜西及周邊中強(qiáng)地震較不活躍，記錄到的地震信息較少，無(wú)較多震例對(duì)應(yīng)及分析，在以后的數(shù)據(jù)應(yīng)用中，應(yīng)注意與鄰區(qū)的地電監(jiān)測(cè)資料聯(lián)合應(yīng)用。

2）筆者利用長(zhǎng)趨勢(shì)變化分析方法、傅氏滑動(dòng)去年周期方法、歸一化月速率方法、差分方法對(duì)寶雞臺(tái)2015—2018年的地電阻率數(shù)據(jù)進(jìn)行了處理和系統(tǒng)分析，九寨溝Ms7.0級(jí)地震前，寶雞臺(tái)地電阻率可能存在一定的前兆異常。

3）存在的問(wèn)題：由于地電外線路多布設(shè)于農(nóng)田中，干擾較大，地電觀測(cè)資料包含著許多干擾信息，對(duì)資料的應(yīng)用可能會(huì)產(chǎn)生較大影響。

參考文獻(xiàn)：

[1] 程鑫，何康，趙桂寶，等.嘉山臺(tái)地電阻率各類變化的判定與分析[J].防災(zāi)科技學(xué)院學(xué)報(bào)，2011，13（4）：30-34.

[2] 葉青，王曉，杜學(xué)彬，等.中國(guó)地震井下地電阻率研究進(jìn)展[J].吉林大學(xué)學(xué)報(bào)（地球科學(xué)版），2022，52（3）：669-683.

[3] 趙斐，葉青，楊曉鵬，等.甘肅平?jīng)龅仉姳葴y(cè)臺(tái)站建設(shè)及地電阻率數(shù)據(jù)分析[J].中國(guó)地震，2022，38（1）：142-152.

[4] 葉青，張翼，黃興輝，等.2015年內(nèi)蒙古阿拉善Ms5.8地震前地電阻率變化特征研究[J].地震，2017，37（4）：123-133.

[5] HUANG Q H.Retrospective investigation of geophysical data possibly associated with the Ms8.0 Wenchuan earthquake in Sichuan，China[J].J Asian Earth Sci，2011，41（4/5）：421-427.

[6] 馬世虎，王秋寧.寶雞地震臺(tái)地電阻率觀測(cè)效能評(píng)價(jià)[J].地震地磁觀測(cè)與研究，2008，29（6）：65-70.

[7] 王志棟，雷正超，寇俊陽(yáng)，等.張掖MS 5.0地震前山丹地震臺(tái)地電阻率異常分析[J].地震地磁觀測(cè)與研究，2020，41（4）：83-89.

[8] DU X B. Two types of changes in apparent resistivity" in" earthquake prediction[J].SciChina Earth Sci， 2011，54（1）：145-156.

[9] 杜學(xué)彬，孫君嵩，陳軍營(yíng).地震預(yù)測(cè)中的地電阻率數(shù)據(jù)處理方法[J].地震學(xué)報(bào)，2017，39（4）：531-548.