王同利 崔博聞 王麗紅 顏曉曄 李菊珍

摘要：采用歸一化速率變化方法對(duì)2008年汶川MS8.0地震震中周邊19個(gè)定點(diǎn)臺(tái)站觀測(cè)到的地電阻率的時(shí)空變化過(guò)程進(jìn)行了計(jì)算分析。時(shí)序分析結(jié)果顯示：在汶川地震前，震中周邊多個(gè)地電阻率臺(tái)站準(zhǔn)同步地出現(xiàn)了中短期前兆異常;空間分布過(guò)程顯示，地震孕震過(guò)程中地電阻率在震中區(qū)域有負(fù)異常叢集和由遠(yuǎn)及近的遷移現(xiàn)象，負(fù)異常叢集呈象限分布，長(zhǎng)軸方向與震源機(jī)制解主壓力方向基本吻合，且受孕震斷裂、臺(tái)站位置、觀測(cè)裝置布設(shè)方向等的影響。

關(guān)鍵詞：汶川MS8.0地震;地電阻率;時(shí)序分析;空間分布;歸一化速率變化方法

中圖分類號(hào)：P315.722 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1000-0666（2022）01-0066-09doi：10.20015/j.cnki.ISSN1000-0666.2022.0008

0 引言

2008年5月12日14時(shí)28分，四川省汶川縣發(fā)生MS8.0強(qiáng)烈地震，造成69 000多人死亡，500多萬(wàn)間房屋倒塌，有感面積達(dá)250萬(wàn)km2，主震后又發(fā)生了多次余震。震后眾多的學(xué)者研究了可能與汶川地震孕震有關(guān)的前兆異常變化（杜學(xué)彬，2010;解滔等，2018;Luet al，2016;張學(xué)民等，2009;馮志生，2013;高曙德等，2010;何案華，2012;邱澤華等，2010）。

地電阻率是地下探測(cè)范圍內(nèi)介質(zhì)電阻率的綜合反映，其在地震孕震過(guò)程中表現(xiàn)出的異常變化主要反映構(gòu)造應(yīng)力作用下，測(cè)區(qū)介質(zhì)變形誘發(fā)的微裂隙活動(dòng)引起的電阻率變化，在靠近震中區(qū)的地電阻率異常形態(tài)常以負(fù)異常為主，并伴有年變畸變（杜學(xué)彬等，2007）。震前電阻率異常與震源機(jī)制有關(guān)，其升、降變化取決于臺(tái)址下介質(zhì)的受力狀態(tài)（如膨脹或壓縮等），巖層電阻率對(duì)微小應(yīng)變有放大效應(yīng)，通?？煞糯?個(gè)數(shù)量級(jí)。在地震發(fā)生前的短期階段，由于斷層預(yù)滑，近斷層破裂區(qū)的地下介質(zhì)應(yīng)變加速積累，地電阻率在原有的變化背景上通常會(huì)出現(xiàn)加速變化，而在近震中區(qū)應(yīng)變被部分釋放，地電阻率通常表現(xiàn)為轉(zhuǎn)折的正異常變化（汪志亮等，2002;Du，2011）。對(duì)巖石進(jìn)行應(yīng)力加載實(shí)驗(yàn)也表明，隨著應(yīng)力的增大巖石電阻率呈下降變化，同時(shí)在相同應(yīng)力加載作用下，巖石不同方向的電阻率呈現(xiàn)出各向異性變化（陸陽(yáng)泉等，1990）。

我國(guó)的地電阻率觀測(cè)是在邢臺(tái)地震后快速發(fā)展起來(lái)的，通常采用對(duì)稱四極裝置，觀測(cè)供電極距在800～2 400 m，同一臺(tái)站在地表布設(shè)2～3個(gè)不同方向的測(cè)道，以整點(diǎn)觀測(cè)的方式進(jìn)行長(zhǎng)期連續(xù)觀測(cè)。在超過(guò)50年的觀測(cè)實(shí)踐中，多次記錄到了發(fā)生在臺(tái)網(wǎng)內(nèi)和臺(tái)網(wǎng)附近的中強(qiáng)地震震中區(qū)的顯著地電阻率中短期異常（錢復(fù)業(yè)等，1982;錢家棟等，1985;桂燮泰等，1989;汪志亮等，2002;杜學(xué)彬等，2015;Lu et al，2016;解滔等，2018;王同利等，2020）。2008年汶川MS8.0地震前，其震中周邊大范圍內(nèi)臺(tái)站記錄的地電阻率出現(xiàn)了較為顯著的異常變化（杜學(xué)彬，2010;Zhang et al，2009;Huang，2011;錢家棟等，2013;朱濤，2013;Lu et al，2016）。為研究汶川MS8.0地震前后地電阻率的時(shí)空動(dòng)態(tài)變化，本文對(duì)其震中周邊19個(gè)定點(diǎn)臺(tái)站的地電阻率數(shù)據(jù)進(jìn)行歸一化速率變化分析，研究地電阻率的時(shí)序變化過(guò)程和空間分布動(dòng)態(tài)過(guò)程，探討地電阻率在強(qiáng)地震發(fā)生前后的異常變化、各向異性以及遷移過(guò)程等。

1 研究區(qū)概況

汶川MS8.0地震發(fā)生在青藏高原東緣的龍門山構(gòu)造帶內(nèi)，地震破裂面長(zhǎng)達(dá)330 km。龍門山斷裂帶呈NE—SW向延伸，自新生代以來(lái)發(fā)生過(guò)多次強(qiáng)烈的地殼變形和破裂，致使該地區(qū)不同走向和規(guī)模的斷裂縱橫交錯(cuò)，其西側(cè)是鮮水河斷裂帶，北側(cè)是西秦嶺北緣斷裂帶，西北側(cè)是東昆侖斷裂帶（圖1）。這些斷裂帶的地震活動(dòng)頻度和強(qiáng)度均較高，其中鮮水河斷裂帶就曾多次發(fā)生7級(jí)以上大地震。

本文研究區(qū)范圍為四川及周邊的甘肅、云南、寧夏和陜西等地區(qū)（25°～40°N，95°～110°E），選取2003—2009年汶川MS8.0地震震中周邊500 km范圍內(nèi)共19個(gè)臺(tái)站（圖1）記錄的地電阻率數(shù)據(jù)。臺(tái)間距為50～200 km，地表觀測(cè)極距為1 000～2 000 m，電極埋深為2 m;井下觀測(cè)極距為80～200 m，電極埋深為100～200 m。

2 研究方法

歸一化速率變化方法（Normalized Variation Rate Method，簡(jiǎn)稱NVRM）主要用于處理原始時(shí)序數(shù)據(jù)曲線，消除年變化曲線等傳統(tǒng)方法不易識(shí)別的“弱”幅度異常（杜學(xué)彬等，2001）。計(jì)算過(guò)程為：首先對(duì)時(shí)序觀測(cè)數(shù)據(jù)中的個(gè)別數(shù)據(jù)進(jìn)行預(yù)處理，然后再按固定步長(zhǎng)ξ個(gè)數(shù)據(jù)組成子序列，對(duì)子序列進(jìn)行線性回歸并求得回歸直線的斜率ki和相關(guān)系數(shù)Ri（第i個(gè)子序列），由ki·Ri得到初始速率值ρ·′i，以ξ為步長(zhǎng)依次向后滑動(dòng)計(jì)算初始速率時(shí)間序列;對(duì)初始速率進(jìn)行均值為“0”“1”的歸一化處理，得到最終的時(shí)間序列（即ρ·值序列）。杜學(xué)彬等（2001）通過(guò)對(duì)196次MS3.2～7.9地震（94%以上為MS≥4.0地震）的計(jì)算表明前兆異常閾值為±2.4，無(wú)量綱。NVRM方法的關(guān)鍵是首先要對(duì)觀測(cè)資料進(jìn)行預(yù)處理，去除突跳數(shù)據(jù)、年變化、非年變周期等，對(duì)存在階段性長(zhǎng)趨勢(shì)變化的資料要根據(jù)實(shí)際異常落實(shí)情況進(jìn)行分時(shí)段處理。該方法已得到廣泛應(yīng)用，眾多學(xué)者獲得了地震前后的地電阻率時(shí)序異常（劉君等，2013;史紅軍等，2014;徐錫泉等，2014;閻睿等，2015;李姜，2019）。

NVRM原理如下，首先針對(duì)時(shí)間軸計(jì)算斜率：

ki=∑nj=1Tj∑nj=1yj-n∑nj=1Tjyj（∑nj=1Tj）2-n∑nj=1Tj2（1）

再計(jì)算自相關(guān)系數(shù)：

Ri=∑nj=1Tjyj-1n（∑nj=1Tj∑nj=1yj）∑nj=1Tj2-1n（∑nj=1Tj）2·∑nj=1y2j-1n（∑nj=1yj）2（2）

最后計(jì)算得出歸一化速率變化法序列為：

Si=Ri×Ki/σn-1 （i=n，n+1，…，N）（3）

式中：n為滑動(dòng)步長(zhǎng);N為資料長(zhǎng)度;σn-1為N-n個(gè)Ri×Ki時(shí)間序列的均方差，{y}是等間隔前兆數(shù)據(jù)時(shí)間序列，{T}是相應(yīng)的等間隔時(shí)間序列。如果利用月均法曲線滑動(dòng)，Si是月速率;如果利用五日均法曲線滑動(dòng)，Si是五日變化速率。

3 地電阻率時(shí)序變化

對(duì)2003—2009年汶川地震震中周邊500 km范圍內(nèi)臺(tái)站定點(diǎn)的地電阻率的原始觀測(cè)曲線進(jìn)行分析，并進(jìn)行NVRM分析，結(jié)果顯示：成都、冕寧、江油、天水、甘孜和武都臺(tái)的地電阻率長(zhǎng)程時(shí)序曲線在震前均出現(xiàn)了中短期的下降變化異常（圖2）。圖2中橫虛線為歸一化閾值線，黑色區(qū)域表示超出閾值部分。

成都臺(tái)NE測(cè)道原始數(shù)據(jù)和NVRM計(jì)算數(shù)據(jù)曲線均在2006年中出現(xiàn)年變畸變，原始曲線出現(xiàn)持續(xù)性下降至地震發(fā)生后發(fā)生折返，下降幅度約5.5%（圖2a-1）;NVRM曲線在2007年5月超過(guò)閾值，6月到達(dá)最低值后開始轉(zhuǎn)折回升，汶川地震就發(fā)生在折返上升的過(guò)程中（圖2a-2）;成都臺(tái)NW與NE測(cè)道原始曲線變化基本同步，但變化較為不明顯（圖2b-1），NW測(cè)道NVRM曲線則更清晰地反映了震前地電阻率的變化過(guò)程（圖2b-2）。

冕寧臺(tái)EW測(cè)道原始曲線自2005年12月至2008年2月出現(xiàn)持續(xù)性下降變化，下降幅度約23%（圖2c-1），NVRM曲線顯示下降變化明顯超過(guò)閾值（圖2c-2）。NS和NW 測(cè)道2種曲線都無(wú)顯著變化。

江油臺(tái)NE測(cè)道原始曲線在震前無(wú)明顯變化，震后則出現(xiàn)快速下降-折返上升變化（圖2d-1），NW測(cè)道自2006年8月開始出現(xiàn)下降變化，下降幅度約1.1%（2e-1）;兩測(cè)道的NVRM曲線顯示，在震前的下降變化均沒有超過(guò)閾值，在震后下降超過(guò)閾值后折返上升（圖2d-2、e-2）。

天水臺(tái)EW測(cè)道原始曲線于2008年2月底出現(xiàn)階躍性上升短期變化，上升幅度約1.5%，之后快速下降，汶川地震就發(fā)生在快速下降的過(guò)程中，地震后恢復(fù)至震前的觀測(cè)狀態(tài)，NVRM曲線變化基本同步（圖2f）。

自2006年開始，甘孜臺(tái)NE和NW測(cè)道原始曲線在震前2～3 a就開始出現(xiàn)持續(xù)趨勢(shì)下降，且年變形態(tài)清晰，震前累計(jì)下降幅度分別約為4.9%和8.3%，地震發(fā)生后原有變化趨勢(shì)也未改變（圖2g-1、h-1）;NVRM曲線顯示，在2006年底，兩測(cè)道均出現(xiàn)下降變化，超過(guò)閾值后發(fā)生折返上升，于2007年5—6月上升到最高值后下降，同時(shí)NW測(cè)道年變形態(tài)發(fā)生改變，于2008年1月又開始出現(xiàn)下降變化，略折轉(zhuǎn)上升至3月再次下降，在汶川地震發(fā)生后又下降（圖2g-2、h-2）。

武都臺(tái)EW測(cè)道原始曲線年變形態(tài)清晰穩(wěn)定，自2007年2月開始呈現(xiàn)趨勢(shì)上升，未出現(xiàn)年變畸變現(xiàn)象，地震后其變化趨勢(shì)也未改變，截至2009年6月上升幅度約8%，之后折返（圖2i-1）;NVRM曲線顯示，震前上升超過(guò)閾值現(xiàn)象明顯（圖2i-2），地震發(fā)生在折返下降過(guò)程中，NS和NE測(cè)道震前無(wú)明顯異常變化。

4 空間異常動(dòng)態(tài)分析

空間異常動(dòng)態(tài)過(guò)程，是指異常的空間分布圖像隨時(shí)間遷移變化的過(guò)程。電阻率觀測(cè)的空間異常動(dòng)態(tài)是由區(qū)域性應(yīng)力場(chǎng)加強(qiáng)和介質(zhì)不均勻性引起的，即除震源區(qū)是一個(gè)高應(yīng)力集中區(qū)外，周圍介質(zhì)還會(huì)形成多個(gè)集中點(diǎn)，這些應(yīng)力集中點(diǎn)多分布在斷層附近，介質(zhì)不均勻性較為明顯，且對(duì)震源區(qū)的應(yīng)力變化較為靈敏。因此，當(dāng)震源區(qū)應(yīng)力場(chǎng)加強(qiáng)到一定程度時(shí)，周圍介質(zhì)就會(huì)出現(xiàn)多個(gè)局部地段應(yīng)力集中，形成區(qū)域應(yīng)力場(chǎng)異常，這種反映震源區(qū)域的應(yīng)力場(chǎng)變化通常會(huì)導(dǎo)致地電阻率的明顯異常，一般以地電阻率下降為主。強(qiáng)震發(fā)生時(shí)地電阻率下降的異常區(qū)范圍半徑約一二百千米。桂燮泰等（1987）、金安忠（1982）分析結(jié)果顯示在1976年唐山MS7.8地震前，華北地區(qū)相當(dāng)大的范圍內(nèi)地電阻率觀測(cè)值就出現(xiàn)了負(fù)異常隨時(shí)間、空間的動(dòng)態(tài)變化。空間上常常表現(xiàn)為以震中區(qū)為中心向外擴(kuò)展，異常變化幅度從震中區(qū)向外逐漸由大變小，且分布不均勻;同一觀測(cè)站點(diǎn)的異常變化動(dòng)態(tài)過(guò)程為隨著地震發(fā)生時(shí)間臨近，負(fù)異常變化幅度愈來(lái)愈大。王同利等（2020）分析了2017年九寨溝MS7.0地震前地電阻率的空間動(dòng)態(tài)演化過(guò)程，結(jié)果顯示地電阻率負(fù)異常空間分布隨時(shí)間遷移，且和地震震中區(qū)有很好的關(guān)聯(lián)。

本文對(duì)汶川地震震中周邊19個(gè)臺(tái)站地電阻率觀測(cè)數(shù)據(jù)進(jìn)行NVRM分析，得到了相對(duì)月變化，再進(jìn)行空間等值線動(dòng)態(tài)分析（圖3、4），計(jì)算中考慮大多數(shù)臺(tái)站的布極為NS和EW兩個(gè)測(cè)道，因此空間等值線動(dòng)態(tài)分析為NS和EW兩個(gè)方向，布極方式?jīng)]有這兩個(gè)方向的臺(tái)站，分別用各臺(tái)和這兩個(gè)方向之間夾角最小的測(cè)道在NS和EW向上的投影計(jì)算得到的地電阻率變化量進(jìn)行替代。

圖3顯示了2007年11月至2008年9月汶川地震震中周邊地區(qū)地電阻率NS測(cè)道的空間動(dòng)態(tài)變化過(guò)程。從圖中可以看出，2007年11月震中區(qū)周邊地電阻率變化平穩(wěn);2008年1月在天水、通渭和武都附近出現(xiàn)了小區(qū)域的負(fù)異常變化;2月負(fù)異常變化范圍逐漸擴(kuò)大，向NE-SW向延伸，SW向延伸至甘孜附近，NE向延伸至寶雞附近;3月負(fù)異常范圍進(jìn)一步擴(kuò)大且異常更加明顯，出現(xiàn)顯著的NE-SW向負(fù)異常條帶;4～5月負(fù)異常在震中區(qū)附近積聚，但相比3月有減弱的趨勢(shì)。5月12日汶川MS8.0地震發(fā)生，6月地電阻率負(fù)異常條帶現(xiàn)象減弱趨于消失，在7～8月出現(xiàn)近NS向的負(fù)異常條帶并向南遷移，2008年8月30日四川仁和發(fā)生MS6.1地震。綜上，隨著地震發(fā)生時(shí)間臨近，負(fù)異?，F(xiàn)象首先在震中區(qū)附近形成，然后面積逐漸變大并向外延伸形成條帶，且以震中區(qū)為中心呈象限分布，地震發(fā)生后異常條帶逐漸變?nèi)跸Щ蛘甙l(fā)生遷移。

圖4顯示了2007年11月至2008年9月汶川MS8.0地震震中周邊地區(qū)EW測(cè)道地電阻率的空間動(dòng)態(tài)變化過(guò)程。從圖中可以看出，2007年11月震中周邊地區(qū)地電阻率的空間變化基本平穩(wěn)，沒有負(fù)異?，F(xiàn)象;2008年1—4月在震中區(qū)附近出現(xiàn)由弱變強(qiáng)的NE-SW向負(fù)異常條帶，負(fù)異常區(qū)域首先出現(xiàn)在通渭和平?jīng)龈浇缓笙騈E和SW兩個(gè)方向延伸，NE向延伸至銀川，SW向經(jīng)過(guò)成都、甘孜延伸至冕寧和紅格附近。2008年5月汶川地震發(fā)生后，負(fù)異常條帶減弱消失，6～9月僅有零散的負(fù)異常，大部分地區(qū)表現(xiàn)為正異常或正常變化。

EW測(cè)道的地電阻率空間動(dòng)態(tài)變化過(guò)程與NS測(cè)道基本一致，但異常變化稍顯明顯。EW測(cè)道負(fù)異常集中區(qū)域范圍大于NS測(cè)道，負(fù)異常強(qiáng)度也較NS測(cè)道強(qiáng)。地電阻率負(fù)異?？臻g動(dòng)態(tài)變化過(guò)程與地震孕育中的應(yīng)力積累動(dòng)態(tài)變化有關(guān)，且明顯地表現(xiàn)為各向異性。這種各向異性現(xiàn)象與九寨溝MS7.0、唐山MS7.8地震前地電阻率出現(xiàn)的負(fù)異常集中現(xiàn)象較為相似。汶川地震負(fù)異常集中區(qū)遷移方向、異常幅度增強(qiáng)方式與九寨溝MS7.0地震基本一致，這種異常遷移、震級(jí)不同、異常幅度也不同可能與地震斷層錯(cuò)動(dòng)方式和震中區(qū)巖性有關(guān)，兩次地震都發(fā)生在龍門山斷裂帶附近，汶川MS8.0地震的震源機(jī)制解顯示發(fā)震斷層以逆沖和走滑為主，震中區(qū)以花崗巖分布為主（易桂喜等，2012;王振榮，蘭江華，2008）。

5 結(jié)論

汶川MS8.0地震震中周邊地區(qū)地電阻率觀測(cè)是我國(guó)相對(duì)密集的地區(qū)之一。震中周邊共有19個(gè)地電阻率定點(diǎn)觀測(cè)臺(tái)站，多個(gè)臺(tái)站在地震發(fā)生前2～3 a出現(xiàn)了準(zhǔn)同步下降變化。結(jié)合中強(qiáng)地震前地電阻率的孕震機(jī)理分析，獲得了以下主要結(jié)論：

（1）在汶川地震前，成都、冕寧、江油、天水、甘孜、武都臺(tái)地電阻率出現(xiàn)下降異常變化，從孕震機(jī)理上符合地電阻率地震異常變化機(jī)理，且臺(tái)站距離震中區(qū)域較近，這種變化是汶川地震前兆的可能性較大。

（2）地震前，震中周邊地區(qū)一個(gè)相當(dāng)大的范圍內(nèi)，地電阻率觀測(cè)值形成了以震中區(qū)為中心的電阻率負(fù)異常區(qū)域。異常空間分布表現(xiàn)出比較顯著的條帶各向異性，長(zhǎng)軸方向和震源機(jī)制解P軸、發(fā)震斷層破裂方向等基本一致。

（3）空間異常變化分布形狀取決于斷層方向、臺(tái)站和震中區(qū)的位置，并且與震中區(qū)附近的臺(tái)站數(shù)和臺(tái)站分布有關(guān)，異常變化幅度大小和臺(tái)站與震中區(qū)距離有關(guān)。

（4）汶川地震前地電阻率的負(fù)異常區(qū)域在NE—SW方向較為突出，負(fù)異常變化的長(zhǎng)軸方向與震源機(jī)制解主壓力方向吻合，且EW測(cè)道條帶異常較NS測(cè)道明顯。這可能與在孕震晚期，近震中區(qū)主壓應(yīng)力擠壓占優(yōu)勢(shì)有關(guān)。地震前，孕震區(qū)域地下介質(zhì)內(nèi)部裂隙發(fā)育、裂隙走向沿主壓應(yīng)力方向優(yōu)勢(shì)排列，導(dǎo)電介質(zhì)重新分布，導(dǎo)致地電阻率出現(xiàn)以快速下降為主的各向異性。不同觀測(cè)臺(tái)站或者同一觀測(cè)臺(tái)站不同測(cè)道的各向異性現(xiàn)象，是由于在孕震過(guò)程中震中區(qū)域應(yīng)力加強(qiáng)引起發(fā)震斷層附近不均勻介質(zhì)電阻率的變化所導(dǎo)致，與孕震斷層活動(dòng)密切有關(guān)。

地電阻率時(shí)序分析能清晰地反映出汶川MS8.0地震孕震過(guò)程中斷層的應(yīng)力變化過(guò)程?？臻g動(dòng)態(tài)演化過(guò)程分析是多臺(tái)觀測(cè)數(shù)據(jù)變化的集合，因而單臺(tái)數(shù)據(jù)受周邊環(huán)境的干擾影響可以忽略。時(shí)序分析和空間異常動(dòng)態(tài)分析相結(jié)合，能夠更有效地揭示出中強(qiáng)地震的發(fā)生要素，對(duì)汶川地震的發(fā)震機(jī)理、震中區(qū)域的認(rèn)識(shí)等有很好的幫助，但分析結(jié)果中難免存在因臺(tái)站密度低而造成偏差的情況?？傊?，汶川地震發(fā)生前，地電阻率出現(xiàn)長(zhǎng)時(shí)程負(fù)異常變化、空間叢集現(xiàn)象真實(shí)地反映了強(qiáng)震孕育過(guò)程中發(fā)震斷層應(yīng)力及周邊應(yīng)力應(yīng)變的變化，它和地震活動(dòng)性等變化有很好的關(guān)聯(lián)性。分析地電阻率的時(shí)空變化過(guò)程對(duì)認(rèn)識(shí)中強(qiáng)地震孕育過(guò)程及震源區(qū)應(yīng)力場(chǎng)變化具有很大的意義。

參考文獻(xiàn)：

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

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

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

杜學(xué)彬.2010.在地震預(yù)報(bào)中的兩類視電阻率變化[J].中國(guó)科學(xué)：地球科學(xué)，40（10）：1321-1330.

馮志生，李鴻宇，張秀霞，等.2013.地磁諧波振幅比異常與強(qiáng)地震[J].華南地震，33（3）：9-15.

高曙德，湯吉，杜學(xué)彬，等.2010.汶川8.0級(jí)地震前后電磁場(chǎng)的變化特征[J].地球物理學(xué)報(bào)，53（3）：512-525.

桂燮泰，關(guān)華平，戴經(jīng)安.1989.唐山、松潘地震前視電阻率短臨異常圖像重現(xiàn)性[J].西北地震學(xué)報(bào)，11（4）：71-75.

桂燮泰，關(guān)華平，戴經(jīng)安，等.1987.地電視電阻率法測(cè)報(bào)強(qiáng)震指標(biāo)討論[J].華南地震，7（2）：56-63.

何案華，趙剛，劉成龍，等.2012.青海玉樹與德令哈地?zé)嵊^測(cè)井在汶川與玉樹地震前的異常特征[J].地球物理學(xué)報(bào)，55（4）：1261-1268.

解滔，劉杰，盧軍，等.2018.2008年汶川MS8.0地震前定點(diǎn)觀測(cè)電磁異?；厮菪苑治鯷J].地球物理學(xué)報(bào)，61（5）：1922-1937.

金安忠.1982.唐山地震前近震中區(qū)地電阻率的震前突變現(xiàn)象[J].地震學(xué)報(bào)，3（2）：169-173.

李姜，喬子云，張國(guó)苓，等.2019.華北地區(qū)地電阻率歸一化速率異常分析[J].地震地磁觀測(cè)與研究，40（6）：62-71.

劉君，杜學(xué)彬，范瑩瑩，等.2013.甘肅岷縣漳縣MS6.6地震前的地電阻率變化[J].地震工程學(xué)報(bào)，35（4）：819-826.

陸陽(yáng)泉，錢家棟，劉建毅.1990.大型花崗巖標(biāo)本緩慢膨脹破裂過(guò)程中電阻率和聲發(fā)射前兆特征的實(shí)驗(yàn)研究[J].西北地震學(xué)報(bào)，12（2）：35-41.

錢復(fù)業(yè)，趙玉林，余謀明，等.1982.地震前地電阻率的異常變化[J].中國(guó)科學(xué)：化學(xué)，12（9）：831-839.

錢家棟，陳有發(fā)，金安忠.1985.地電阻率法在地震預(yù)報(bào)中的應(yīng)用[M].北京：地震出版社，10-51.

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

邱澤華，張寶紅，池順良，等.2010.汶川地震前姑咱臺(tái)觀測(cè)的異常應(yīng)變變化[J].中國(guó)科學(xué)：地球科學(xué)，40（8）：1031-1039.

史紅軍，趙衛(wèi)星，張可佳，等.2014.吉林省地電阻率震兆異常分析研究[J].地震學(xué)報(bào)，36（3）：452-463，532.

汪志亮，鄭大林，余素榮.2002.地震地電阻率前兆異?，F(xiàn)象[M].北京：地震出版社，58-65.

王同利，崔博聞，葉青，等.2020.九寨溝MS7.0地震地電阻率變化時(shí)空演化分析[J].地球物理學(xué)報(bào)，63（6）：2345-2356.

王振榮，蘭江華.2008.四川汶川大地震的構(gòu)造分析[J].礦物巖石，28（2）：1-10.

徐錫泉，高昌志，王亮.2014.內(nèi)蒙古寶昌臺(tái)地電阻率長(zhǎng)期觀測(cè)數(shù)據(jù)研究[J].地震工程學(xué)報(bào)，36（2）：405-412.

閆睿，朱石軍，胡樂銀，等.2015.北京市延慶地震臺(tái)地電阻率長(zhǎng)程觀測(cè)數(shù)據(jù)研究[J].地震工程學(xué)報(bào)，37（3）：724-730，771.

易桂喜，龍鋒，張致偉.2012.汶川MS8.0地震余震震源機(jī)制時(shí)空分布特征[J].地球物理學(xué)報(bào)，55（4）：1213-1227.

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

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

Du X B.2011.Two types of changes in apparent resistivity in earthquake prediction[J].Science China（Earth Science），54（1）：145-156.

Huang Q H.2011.Retrospective investigation oF geophysical data possible associated with the MS8.0 Wenchuan earthquake in Sichuan，China[J].J Asian Earth Sci，41（4-5）：421-427.

Lu J，Xie T，Li M， et al.2016.Monitoring shallow resistivity changes prior to the 12 May 2008 M8.0 Wenchuan earthquake on the Longmen Shan tectonic zone，China[J].Tectonophysics，675：244-257.

Zhang X M，Ding J H，Sheng X H， et al.2009.Electromagnetic perturbations before Wenchuan M8.0 earthquake and stereo electromagnetic observation system[J].Chinese Journal of Radio Science，24（1）：1-8.

Temporal and Spatial Evolution of the Apparent Resistivity beforeand after the 2008 Wenchuan MS8.0 Earthquake

WANG Tongli1，CUI Bowen1，WANG Lihong1，YAN Xiaoye2，LI Juzhen1

（1.Beijing Earthquake Agency，Beijing 100080，China）（2.Sichuan Earthquake Agency，Chengdu 610000，Sichuan，China）

Abstract

In this paper，the Normalized Variation Rate method is used to analyze the sequential and spatial variation of the apparent resistivity observed by 19 stations around the Wenchuan MS8.0 earthquake.The sequential analysis shows that some of the earthquake stations had quasi-synchronously recorded the short- and medium-term precursory anomalies of the apparent resistivity before the Wenchuan earthquake.The spatial analysis shows that in the epicentral region，the negative anomalies of the apparent resistivity had been featured with a strip distribution and with a migration from far to near before the earthquake.The cluster of the anomalies displayed a quadrant pattern，whose long axis was consistent with the direction of principal compressive stress of the focal mechanism solution.The seismogenic fault，the station location，and the layout direction of the observational instruments also influenced the distribution direction of the long axis of the precursory anomalies.

Keywords：the Wenchuan MS8.0 earthquake;apparent resistivity;time sequential analysis;spatial distribution;the Normalized Variation Rate method