張小濤，韓麗萍，張新東，張雙鳳，暢國平，尹小兵

（1.吉林大學(xué)，吉林 長春 130026；2.河北省地震局邯鄲中心臺(tái)，河北 邯鄲 056001；3.河北省地震局紅山地震臺(tái)，河北 隆堯 055350）

0 引言

在地震孕育過程中，孕震區(qū)介質(zhì)性質(zhì)的變化會(huì)導(dǎo)致波速的變化，研究波速變化可以提供地震孕育的震兆信息。我國自開展波速比研究以來，已經(jīng)取得了很多成果。刁桂苓、周仕勇等討論了岫巖地震序列的平均波速比前兆異常［1］。劉繼祿等分析了張北6.2級(jí)地震前后源區(qū)與場(chǎng)區(qū)波速比的時(shí)空變化特征及其異常［2］。梅世蓉、陸遠(yuǎn)忠等也有利用波速比異常做出預(yù)測(cè)的震例［3－4］。啜永清、宋美琴采用格點(diǎn)嘗試法、和達(dá)法分別計(jì)算了山西大同盆地的平均波速比，探討了該區(qū)平均波速比在中強(qiáng)震前后的變化規(guī)律，得出山西大同盆地地震波速比在主壓應(yīng)力P軸取向穩(wěn)定的條件下同步降低的負(fù)異常與該區(qū)中強(qiáng)地震有一定對(duì)應(yīng)關(guān)系的結(jié)論［5］。趙晉明等從實(shí)驗(yàn)上確證巖石在破裂前，波速比會(huì)發(fā)生明顯變化［6］。此外，波速比也用來分析地殼和地幔的結(jié)構(gòu)。黎明曉、張曉東計(jì)算了首都圈地區(qū)的平均波速比并應(yīng)用多臺(tái)法測(cè)定了華北地區(qū)地殼的平均波速比［7］。馮建剛等利用和達(dá)法計(jì)算了青藏高原東北緣平均波速比［8］。本文開展對(duì)河北省邯鄲—邢臺(tái)地區(qū)波速比和泊松比的空間分布特征研究，能夠?qū)μ崛∫欢ǖ恼鹫仔畔⑻峁┮罁?jù)，對(duì)了解該地區(qū)的孕震環(huán)境和地質(zhì)構(gòu)造介質(zhì)狀況也有參考價(jià)值。

1 方法及原理

假如震中周圍有n個(gè)臺(tái)的資料，則以和達(dá)法為基礎(chǔ)，利用不同地震臺(tái)記錄到的近震直達(dá)Pg、Sg波到時(shí)TPi、TSi，通過最小二乘法計(jì)算平均波速比：

式中ΔTi＝TSi－TPi。

相關(guān)系數(shù)R由下式計(jì)算：

根據(jù)求出的波速比γ，由波速比和泊松比之間的關(guān)系，可得到地殼的泊松比值

γ值的增加或減小反映了巖性或者巖性物理狀態(tài)及化學(xué)成分的變化。

2 觀測(cè)資料情況

本文研究的區(qū)域范圍為 35.0°～38.0°N，113.0°～116.0°E。研究區(qū)位于中、小地震較為活躍的地震帶，該地區(qū)也是邯鄲地震臺(tái)網(wǎng)監(jiān)控最好的地區(qū)。分析了區(qū)域內(nèi)9個(gè)遙測(cè)臺(tái)站所記錄的2001年1月－2008年12月間發(fā)生的535次ML≥1.5地震，震中距都在200km以內(nèi)。要求每個(gè)地震保證至少有4個(gè)臺(tái)站的清晰記錄，去掉波形重疊、難以辨認(rèn)的地震，經(jīng)過認(rèn)真篩選后，用于計(jì)算的ML≥1.5地震413次，震源深度分布在3～25km范圍內(nèi)，平均深度為11.4km。研究區(qū)震中和斷裂帶分布見圖1，各臺(tái)站儀器信息主要參數(shù)見表1。用于分析的地震震源參數(shù)采用雙差精定位結(jié)果［9］。為了研究方便，我們將研究區(qū)按1°×1°分為A1～A9共計(jì)9個(gè)小區(qū)，計(jì)算每個(gè)分區(qū)內(nèi)所有地震的波速比，得到該區(qū)的一個(gè)平均波速比值。

圖1 研究區(qū)震中與臺(tái)站分布圖Fig.1 Distributions of earthquake epicentres and seismic stations in the study area.

3 結(jié)果與分析

研究范圍內(nèi)在A1～A9區(qū)中地震事件全部在6個(gè)以上，其中1966年邢臺(tái)7.2級(jí)大地震震中（37.5°N ，115.1°E）位于 A2區(qū)，1831年磁縣7.5級(jí)大地震震中（36.4°N，114.2°E）位于 A5區(qū)，1981年菏澤5.9級(jí)地震震中（35.3°°N，115.2°°E）位于 A9區(qū)，這三個(gè)區(qū)內(nèi)地震事件較多，最多的小區(qū)A5達(dá)到了149個(gè)，表2給出各個(gè)小區(qū)計(jì)算的事件個(gè)數(shù)和結(jié)果。9個(gè)小區(qū)波速比值在1.696～1.737之間，相關(guān)系數(shù)在0.9883～0.9991之間，研究區(qū)平均波速比值為1.719，略低于地殼的平均波速比1.732。

根據(jù)式（3）計(jì)算得到了研究區(qū)內(nèi)的泊松比值的分布特征。平均泊松比值為0.244，波速比值和泊松比值的空間分布見圖2。從圖2可以看出，研究區(qū)內(nèi)波速比和泊松比在37°N兩側(cè)呈現(xiàn)反向的變化，上側(cè)的A1、A2、A3三個(gè)小區(qū)內(nèi)，自東向西逐漸減??；下側(cè)自東向西逐漸增大，而37°N緯度線把1966年邢臺(tái)震區(qū)和1831年磁縣震區(qū)分開，這種分布應(yīng)該與地殼的速度結(jié)構(gòu)和巖石的物理性質(zhì)有關(guān)。研究區(qū)有很多斷裂帶分布，其中包括新河斷裂、磁縣—大名斷裂、林縣斷裂、聊城—蘭考斷裂等發(fā)生過大地震的斷裂帶，這些復(fù)雜的地質(zhì)環(huán)境通過平均波速比和泊松比都有所反映，平均波速比較為客觀的反映了地下介質(zhì)在空間位置上的差異。

表1 使用的數(shù)字地震臺(tái)站信息表

表2 各小區(qū)內(nèi)平均波速比及相關(guān)系數(shù)

圖2 研究區(qū)內(nèi)波速比和泊松比的空間分布圖Fig.2 Distributions of wave velocity ratio values and Poisson＇s ratio values in the study areas.

為了直觀地揭示邯鄲－邢臺(tái)地區(qū)地殼介質(zhì)變化的空間分布特征，我們利用從9個(gè)小區(qū)所得出的均值結(jié)果經(jīng)插值得到該區(qū)波速比和泊松比的等值線分布如圖3所示。從圖3和圖1中地震空間位置分布可以看出研究區(qū)內(nèi)波速比和泊松比分布比較復(fù)雜，存在四個(gè)高值區(qū)域，1966年邢臺(tái)地震震源區(qū)位于其中的一個(gè)高值區(qū)域，1831年的磁縣地震震源區(qū)位于四個(gè)高值區(qū)域的交匯處。從等值線的變化程度看出，邢臺(tái)震源區(qū)等值線值為1.726，磁縣震區(qū)為1.721，略高于研究區(qū)的均值1.719。由于地震多發(fā)生在地殼18km以上，所以這種變化是上地殼介質(zhì)特性的反映。

圖3 研究區(qū)內(nèi)波速比和泊松比的等值線分布圖Fig.3 The isolines of wave velocity ratio and Poisson＇s ratio in the study areas.

實(shí)驗(yàn)研究表明［10］，對(duì)于普通的巖石而言，泊松比的變化范圍是0.20～0.35之間，而且它對(duì)巖石的組成特別敏感，地殼平均泊松比的變化可能由許多地球物理和地球化學(xué)因素所致，且?guī)r石的礦物組成對(duì)泊松比的變化有著重要影響，其中花崗巖泊松比為0.24，閃長巖為0.27，輝長巖0.30［11］。同樣，部分熔融對(duì)P波和S波速度比也有很大的影響，地震波速比值隨著熔融體熔融程度的增加而增大［12］，如部分熔融程度為5% 的花崗巖的泊松比約為0.31，而沒有熔融的花崗巖的泊松比約為0.24。本文研究區(qū)域內(nèi)泊松比平均值0.244，反映平均深度11.4 km，反映了上地殼沒有熔融的介質(zhì)性質(zhì)，本文泊松比計(jì)算結(jié)果與研究區(qū)內(nèi)地質(zhì)構(gòu)造以花崗巖為主一致。

4 結(jié)論與討論

根據(jù)多臺(tái)和達(dá)法，利用邯鄲—邢臺(tái)地區(qū)9個(gè)數(shù)字地震臺(tái)的觀測(cè)數(shù)據(jù)對(duì)該區(qū)波速比和泊松比分布特征進(jìn)行研究，獲得了具有一定參考價(jià)值的結(jié)果：

（1）分析計(jì)算了研究區(qū)的波速比和泊松比值的空間分布特征，9個(gè)小區(qū)波速比值在1.696～1.737之間，相關(guān)系數(shù)在0.988 3～0.999 1之間，研究區(qū)平均波速比值為1.719，平均泊松比值為0.244，本文泊松比計(jì)算結(jié)果與研究區(qū)內(nèi)地質(zhì)構(gòu)造以花崗巖為主一致，同時(shí)也說明計(jì)算結(jié)果是可靠的。

（2）復(fù)雜的地質(zhì)構(gòu)造通過平均波速比和泊松比都有所反映，平均波速比和泊松比較為客觀的反映了地下介質(zhì)在空間位置上的差異，研究區(qū)內(nèi)波速比和泊松比變化比較復(fù)雜，這種變化是上地殼介質(zhì)性質(zhì)變化的體現(xiàn)。

［1］刁桂苓，周仕勇，劉杰，等.岫巖地震序列的平均波速比前兆異常［J］.華北地震科學(xué)，2005，23（4）：1－5.

［2］劉繼錄，張從珍，曹肅朝，等.張北6.2級(jí)地震前后波速比在源區(qū)與場(chǎng)區(qū)的變化特征［J］.華北地震科學(xué)，2001，19（2）：14－20.

［3］梅世蓉，馮益德，張國民，等.中國地震預(yù)報(bào)概論［M］.北京：地震出版社，1993：78－80.

［4］陸遠(yuǎn)忠，陳章立.地震預(yù)報(bào)的地震學(xué)方法［M］.北京：地震出版社，1985：175－196.

［5］啜永清，宋美琴.山西大同盆地震源機(jī)制解矛盾符號(hào)比及波速比變化的研究［J］.山西地震，2005，109（2）：1－3.

［6］趙晉明，胡毅力，楊潤海，等.巖石臨破前波速奇異變化的再研究［J］.地震研究，2001，24（2）：136－139.

［7］黎明曉，張曉東.應(yīng)用多臺(tái)法測(cè)定華北地區(qū)地殼的平均波速比［J］.地震，2004，24（1）：163－169.

［8］馮建剛，代煒，董志平，等.青藏高原東北緣平均波速比的測(cè)定及研究［J］.西北地震學(xué)報(bào)，2009，31（1）：61－65.

［9］張小濤、呂堅(jiān)、馬廣慶，等.雙差地震定位法在邯邢地區(qū)地震精確定位中的初步應(yīng)用［J］.地震研究，2008.31（1）：37－41.

［10］Christensen N I，Mooney W D.Seismic velocity structure and composition of the continental crust：aglobal view［J］.J Geophys.Res.，1995，100：9761－9788.

［11］Tarkov A P，Vavakin V V.Poisson＇s ratio behaviour in crystalline rocks：application to the study of the Earth＇s interior［J］.Phys.Earth Planet.Inter.，1982，29：24－29.

［12］Watanabe T.Effects of water an d melt on seismic velocities and their application to characterization of seismic reflectors［J］.Geophys.Res.Lett.，1993，20：2933－2936.