王飛　陶媛

摘? 要：為研究體應(yīng)變與地電阻率在地震前的相關(guān)變化，以安徽大蜀山地震臺為研究對象，利用臺站體應(yīng)變值和地電阻率值，計算分析兩者數(shù)值的相關(guān)變化系數(shù)及變化曲線，根據(jù)前人大量物理實驗得到的不同巖性下兩者的變化關(guān)系理論模型，確定臺站數(shù)據(jù)的良好對應(yīng)性，為地震短臨預(yù)報提供一種思路。

關(guān)鍵詞：地電阻率? 體應(yīng)變? 相關(guān)變化

中圖分類號：P315 ? ?文獻標識碼：A 文章編號：1672-3791（2019）09（b）-0201-03

“十一五”期間，安徽省地震局為加強短臨跟蹤，沿郯廬斷裂帶中南端安裝了5臺TJ-Ⅱ型鉆孔體應(yīng)變觀測儀，分別安裝在嘉山地震臺、蚌埠地震臺、六安地震臺、合肥地震臺、黃山地震臺，整點值能記錄到完整的固體潮，分鐘值可以清晰地記錄到地震應(yīng)變波、震時的同震應(yīng)變。且在蒙城地震臺、合肥地震臺、黃山地震臺安裝了數(shù)字化地電阻率觀測儀，記錄地震前后電阻率的變化情況，運行狀況良好。

地震發(fā)生前后，兩種手段均有或大或小的數(shù)據(jù)異常，該文根據(jù)在大量物理試驗下得到的巖石體應(yīng)變值及電阻率值隨應(yīng)力變化而變化的理論結(jié)果，對應(yīng)分析大蜀山地震臺的數(shù)據(jù)，確定變化規(guī)律。

1? 原理方法

根據(jù)布雷斯（Brace，W.F，1968）的實驗結(jié)果，當改變軸向加載時，在軸應(yīng)力小于破裂應(yīng)力一半時，巖石的軸向電阻率略有增加;超過這個應(yīng)力值時，巖石的電阻率開始下降，最后在接近破壞應(yīng)力80%時，巖石的電阻率開始大幅度下降。橫向電阻率在較低應(yīng)力下就開始出現(xiàn)明顯下降，與軸向電阻率變化略有不同。

布雷斯在測量巖石電阻率的同時，還測定了相應(yīng)巖石的應(yīng)變。在軸向應(yīng)力達到50%破裂應(yīng)力以前，巖石的體應(yīng)變減小，相當于巖石中裂隙的閉合和空隙的減少;在超過50%破裂應(yīng)力之后，體應(yīng)變反而增大，相當于新生裂隙的產(chǎn)生（擴容）。

通過以上應(yīng)變導致的巖石變形進而引起地電阻率和體應(yīng)變變化的機理，我們可以知道，當受到地下應(yīng)力時，地電阻率和體應(yīng)變值的變化應(yīng)是相反的。

2? 數(shù)據(jù)處理及分析

2.1 數(shù)據(jù)相關(guān)性

選取安徽省大蜀山地震臺2012—2016鉆孔體應(yīng)變及數(shù)字化地電阻率連續(xù)分鐘值預(yù)處理數(shù)據(jù)（2016年后因地鐵運行影響，臺站地電阻率數(shù)據(jù)干擾較大），分析其相關(guān)性，數(shù)據(jù)曲線如圖3所示。

進一步對兩者進行相關(guān)性分析，求取相關(guān)系數(shù)，利用相關(guān)公式：

Cov（X，Y）為X、Y的寫方差，D（X）、D（Y）為X、Y的方差;進而計算出相關(guān)系數(shù)為-0.8，接近-1，則大蜀山地震臺體應(yīng)變值和地電阻率值具有較好的負相關(guān)性。

翻查資料顯示大蜀山地震臺所處地主要巖性為砂巖，由圖1及圖2我們可以發(fā)現(xiàn)，對于砂巖來說，當應(yīng)力低于3×108Pa時，地電阻率值與體應(yīng)變值反向變化;當應(yīng)力在3×108Pa～5×108Pa之間時，地電阻率值與體應(yīng)變值同向變化;大于5×108Pa時，兩者又呈反向變化（Brace，W.F，Electrical Resistivity Changes in Saturated Rocks during Fracture and Frictional Sliding），由此整理出表1。

由此我們可以預(yù)測，因應(yīng)力積累導致巖石破裂這一過程中，兩者的變化關(guān)系應(yīng)當存在這樣的規(guī)律。

2.2 利用地震震例分析

選取2012—2016年合肥周邊發(fā)生的3.0級以上地震作為研究震例，分析在震前震后兩者數(shù)據(jù)間的變化規(guī)律，經(jīng)分析發(fā)現(xiàn)幾乎在每次合肥3.0級地震之前，都存在著地電阻率值與體應(yīng)變值同向變化的異常狀況。這與前面體應(yīng)變與地電阻率變化的規(guī)律，即反向-同向-反向相符合，如圖4所示。

3? 結(jié)語

根據(jù)以上論述可發(fā)現(xiàn)：大蜀山地震臺體應(yīng)變與地電阻率具有較好的負相關(guān)性;兩者數(shù)據(jù)符合特定巖性下反向-同向-反向變化規(guī)律。

但由圖4也可看出，并不是每個地震都有較好的對應(yīng)關(guān)系，且臺站并非都同時具有體應(yīng)變和地電阻率兩種觀測手段，因此該規(guī)律的普遍適用性尚待進一步驗證。該規(guī)律的適用性還受到臺址巖性的影響，對不同的巖性，在應(yīng)力變化時，地電阻率和體應(yīng)變的變化規(guī)律也會不同。

參考文獻

[1] 張?zhí)熘?1.2干巴圍壓下巖樣破裂和摩擦滑動過程中電阻率變化[J].地震學報，1985，7（4）：428-433.

[2] Brace，W.F.Electrical Resistivity Changes in Saturated Rocks during Fracture and Frictional Sliding[J].J.Geophys.Res，1968，73（4）：1433-1445.

[3] 陳鋒，廖發(fā)典，唐獻民.視電阻率法測礦井巖石應(yīng)變實驗的初步總結(jié)[J].地球物理學報，1975（2）：88-99.

[4] 劉川琴.合肥地震臺鉆孔體應(yīng)變干擾分析[J].地震地磁觀測與研究，2014，33（5/6）：208-212.

[5] 李本有.安徽省安慶Ms4.8震前鉆孔體應(yīng)變數(shù)據(jù)變化分析[J].地震地磁觀測與研究，2012（33）：217-220.

[6] 何康，王燚坤.九江M_s5.7地震前地電阻率異常分析[J].地震地磁觀測與研究，2011（3）：62-66

[7] 何康，程鑫.安徽省數(shù)字化地電阻率干擾與短臨異常研究[J].地震地磁觀測與研究，2010（4）.