馬 棟，馮向東，張素欣，王曉山，李均良，王 想

（1.河北省地震局，河北 石家莊 050021；2.河北省地震局 易縣地震臺，河北 易縣 074211）

0 引言

2010年7月30日22時21分、22時42分，在河北省易縣相繼發(fā)生了ML3.8、3.3地震，震中分別為：115.35°E、39.42°N 和115.35°E，39.41°N（重新定位結(jié)果）。兩次地震距易縣地震臺僅7km，易縣臺的形變測項有明顯同震響應(yīng)，周邊的赤城臺、懷來臺、天津薊縣臺也有不同程度的同震形變響應(yīng)。連續(xù)兩次3級地震對200km范圍內(nèi)的形變觀測造成同震響應(yīng)，這在中小地震的同震響應(yīng)方面屬于特殊地震事件。本文對同震形變響應(yīng)的性質(zhì)、幅度及其產(chǎn)生原因進行分析，研究形變測項的同震響應(yīng)能力；結(jié)合地震波形資料，計算震源機制解，并通過震源球面投影、震中地震地質(zhì)構(gòu)造，分析有同震響應(yīng)的臺站、無同震響應(yīng)的臺站與斷層錯動的聯(lián)系。

這兩次3級以上地震雖然屬于中小地震，但由于其發(fā)生的地理位置特殊，屬于首都圈重點監(jiān)視地區(qū)，震中周圍易縣伸縮儀、張家口伸縮儀、陽原水平擺傾斜儀等存在變化。因此充分分析這兩次地震的震源機制、同震形變響應(yīng)，對于判斷該區(qū)震情的后續(xù)發(fā)展具有意義。

1 震中周圍地震活動概況

易縣地震發(fā)生在淶水斷裂附近，震中50km范圍歷史上共記錄5級以上地震3次，即1138年8月22日河北易縣5.5級地震、1658年2月3日河北淶水6.0級地震和1923年9月14日河北新城5.5級地震，最近1次ML3.0以上地震是1993年1月6日河北淶源的ML3.1地震。1057－1970年震中100km范圍內(nèi)歷史地震：5.0～5.9級地震10次，6.0～6.9級地震4次，7.0級地震1次；1970－2010年震中100km范圍內(nèi)地震情況：3.0～3.9級地震41次，4.0～4.9級地震5次，5.0～5.9級地震1次（圖1）。2008年1月－2010年7月，易縣臺記錄到易縣地區(qū)0級以上地震共計11次，其中ML0～0.9地震4次，ML1.0～1.9地震4次，ML2.0～2.9地震1次，ML3.0～3.9地震2次。

圖1 易縣周邊震中分布圖（1057－2010）Fig.1 Distribution of earthquake epicenters in Yixian and around region（1057－2010）.

2 同震形變響應(yīng)與地震波形資料對比分析

2.1 同震形變響應(yīng)情況

震中附近電磁和流體的測項較少，只有部分形變測項對這兩次地震有較好的同震響應(yīng)（圖2）。本文收集了河北省、天津市、北京市、山西省的定點形變資料，其中北京市的形變測項對易縣地震無同震響應(yīng)，主要原因是對于中小地震垂直擺傾斜儀的同震形變響應(yīng)較為明顯，而北京市的定點形變測項中垂直擺傾斜儀的測項較少；山西省的形變測項如代縣、太原距震中200～300km，對易縣地震均無同震形變響應(yīng)，其中只有臨汾臺有垂直擺傾斜儀測項（距震中498km），由于距離較遠，也未能記錄到同震形變響應(yīng)；天津市薊縣臺距震中201km，薊縣垂直擺傾斜儀NS向?qū)σ卓h地震有同震響應(yīng)。

河北省的同震形變響應(yīng)測項共計12項，分別為易縣四分量應(yīng)變儀第2、3、4分量、易縣垂直擺傾斜儀EW、NS向、易縣水管傾斜儀EW向、赤城垂直擺傾斜儀EW、NS向、懷來垂直擺傾斜儀EW、NE、NS向、懷來水管傾斜儀NS向。同震響應(yīng)項數(shù)占全省形變測項的16%，主要集中在距震中較近的易縣臺、赤城臺、懷來臺，其他臺站測項無明顯的同震響應(yīng)。同震形變響應(yīng)測項主要集中在垂直擺傾斜儀（7項），其次是四分量應(yīng)變儀（3個分量），另外還有水管傾斜儀（2項）；而伸縮儀、體應(yīng)變儀、水平擺傾斜儀均無明顯的同震響應(yīng)。按照儀器類別分析，有同震響應(yīng)的同類測項占全省該類測項的百分比為：垂直擺傾斜儀39%、四分量應(yīng)變儀75%、水管傾斜儀13%。圖3為易縣垂直擺傾斜儀分鐘值曲線，有較大幅度的同震階變。

分量應(yīng)變儀的同震響應(yīng)持續(xù)時間為2分鐘，垂直擺傾斜儀的同震響應(yīng)持續(xù)時間為2～10分鐘，水管傾斜儀的同震響應(yīng)持續(xù)時間為2～3分鐘。由于這兩次地震震級較小，能量釋放引起的地震動較弱，因此同震響應(yīng)持續(xù)時間較短。與遠震大震所表現(xiàn)的同震形變波［1－2］不同，易縣地震的同震形變響應(yīng)主要表現(xiàn)為突跳和階躍，其中分量應(yīng)變儀表現(xiàn)為突跳，垂直擺傾斜儀和水管傾斜儀主要表現(xiàn)為階躍（表1）。

表1 易縣ML3.8、3.3地震同震形變響應(yīng)統(tǒng)計表

圖2 易縣地震震中附近形變臺站分布圖Fig.2 Distribution of the deformation observatory stations in the study area.

圖3 2010年7月30－31日易縣垂直擺傾斜儀分鐘值曲線圖Fig.3 Minute value curves of vertical pendulum from July 30to July 31，2010，in Yixian station.

無同震形變響應(yīng)的測項，主要有以下三方面的因素：一是受儀器靈敏度的限制，未能對微小地震動做出響應(yīng)；二是震中距較大；三是部分測項在短時間內(nèi)受到了氣壓等因素的干擾，同震響應(yīng)湮沒在噪聲中，難以對其分析識別。

雖然易縣四分量應(yīng)變儀、垂直擺傾斜儀較其他臺站同類測項的同震響應(yīng)幅度偏大，但表1中的同震響應(yīng)幅度與震中距并無較好的線性統(tǒng)計關(guān)系，甚至出現(xiàn)相反的情況。這是由于近震面波持續(xù)時間短，同震形變波一般難以記錄到準確的最大幅度值［2－5］。遠場大震以彈性形變?yōu)橹?，同震形變響?yīng)表現(xiàn)為形變波，震后恢復(fù)到震前的狀態(tài)和水平；近場中小地震主要以塑性形變?yōu)橹?，引起的同震形變響?yīng)表現(xiàn)為同震階躍，之后形成臺階，未能完全恢復(fù)至震前水平。

2.2 同震形變響應(yīng)與地震波形資料對比分析

兩次地震有兩次對應(yīng)同震響應(yīng)的測項有：易縣四分量應(yīng)變儀的第3分量、赤城垂直擺傾斜儀EW、NS向、懷來垂直擺傾斜儀EW向；其余測項主要是對易縣ML3.8地震有同震響應(yīng)，而對易縣ML3.3地震無同震響應(yīng)。主要因素有兩個：其一是易縣ML3.8地震同震響應(yīng)之后，兩次地震間隔時間短，儀器處于調(diào)整階段，未能對之后的易縣ML3.3地震產(chǎn)生同震響應(yīng)；其二是受儀器靈敏度、震源深度、震級大小等的影響。

同震階變是儀器記錄到的地震期間地球應(yīng)變場的一種準靜態(tài)調(diào)整，地震發(fā)生時，孕震體釋放能量，震源區(qū)應(yīng)變場的調(diào)整是應(yīng)力應(yīng)變的傳遞，震源區(qū)外不能直接記錄到震源區(qū)的應(yīng)變調(diào)整，只能記錄到通過地震動使得各地應(yīng)力場的調(diào)整［6－9］。易縣地震震級小，距震中7km的易縣臺垂直擺傾斜儀有大幅的同震階變，這是記錄到的震源區(qū)應(yīng)變調(diào)整；而赤城臺和懷來臺距震中超過100km、薊縣臺超過200 km，這兩個臺記錄到的同震響應(yīng)不屬于震源區(qū)的應(yīng)力調(diào)整，是地震動使周邊應(yīng)力場產(chǎn)生的調(diào)整，進而記錄到的同震階變。

有同震形變響應(yīng)的4個臺站中，易縣、赤城、懷來、薊縣的震相分析、S波與P波到時差及其同震形變響應(yīng)特征見表2。易縣臺使用的是CTS－1地震計，頻帶50Hz－120s；赤城臺使用的是BBVS－60地震計，頻帶50Hz－60s；懷來臺使用的是FBS－3地震計，頻帶20Hz－20s；上述地震計的采樣率為50～100Hz。易縣垂直擺傾斜儀等形變測項為分鐘值采樣，采樣率為1／60Hz。形變數(shù)據(jù)比地震數(shù)據(jù)的采樣率低幾千倍，中小地震持續(xù)時間短，同震形變響應(yīng)表現(xiàn)為臺階和突跳，未能反映出很好的形變波形，因此不能對同震形變響應(yīng)做震相識別。

從地震波形分析，易縣、赤城、懷來、薊縣的地震波持續(xù)時間為1～3分鐘，而同震形變響應(yīng)測項中有13項持續(xù)時間在2～3分鐘（表1），這與地震計記錄到的地震波形持續(xù)時間一致；有4項持續(xù)時間為6～10分鐘（表1），其曲線形態(tài)為“同震階躍－持續(xù)－部分恢復(fù)”，持續(xù)和恢復(fù)階段主要是受地震波的作用，形變儀器的加載和卸載作用造成的。如不考慮持續(xù)和恢復(fù)階段，同震階躍的時間為1～3分鐘，這與地震計記錄到的地震波持續(xù)時間一致。

表2 易縣ML3.8、3.3地震震相分析、同震形變響應(yīng)特征

3 易縣 ML3.8、3.3地震重新定位結(jié)果

使用hyp2000程序，對易縣 ML3.8、3.3地震進行重新定位，參與定位的21個臺站基本將震中包圍，定位結(jié)果見表3。結(jié)果顯示易縣ML3.8地震速報目錄與重新定位目錄的偏差較大，經(jīng)、緯度定位偏差分別為0.04°、0.03°，速報目錄震中與重新定位震中的偏差主要是由參與定位臺站的數(shù)量、數(shù)據(jù)質(zhì)量以及定位的算法、速度模型的選取等有關(guān)。

表3 2010年7月30日易縣ML3.8、3.3地震重新定位與速報目錄對比

4 震源機制解及其同震形變響應(yīng)臺站震源球面投影分析

4.1 震源機制解

采用點源位錯模型，層狀介質(zhì)速度結(jié)構(gòu)，通過計算理論地震圖，與觀測的地震垂直向記錄的直達P、S波最大振幅擬合，反演震源機制的方法［10］。該方法適用于區(qū)域地震臺網(wǎng)，而且臺少震小的情況。只要具有4個以上臺站完整的垂直向記錄波形，利用振幅比擬合就可以反演得到震源機制［11］。由區(qū)域地震臺網(wǎng)垂直向地震Ap和As振幅資料計算中小地震的震源機制，參考臺站Pg初動選擇最佳結(jié)果。

利用河北省數(shù)字地震臺網(wǎng)記錄的區(qū)域波形資料，選用震中周圍10個臺站。震中位置分別采用河北臺網(wǎng)eqt目錄和重新定位結(jié)果，用垂直向地震Ap和As振幅資料，計算2010年7月30日河北易縣ML3.8、3.3地震的震源機制解參數(shù)，得出不同的解 算結(jié)果（表4、圖4）。

表4 2010年7月30日易縣ML3.8、3.3地震震源機制解參數(shù)

圖4 易縣ML3.8、3.3地震震源機制解Fig.4 Focal mechanism solutions of Yixian ML3.8、ML3.3earthquakes.

重新定位結(jié)果兩次地震位置更為靠近，計算所得的震源機制基本一致。易縣ML3.8地震震源斷層初步判定為NEE向，斷層面參數(shù)為走向69°、傾角38°、滑動角－129°，斷層表現(xiàn)為正斷兼右旋走滑。易縣ML3.3地震震源斷層初步斷定為NEE向，斷層面參數(shù)為走向74°、傾角38°、滑動角－142°，斷層表現(xiàn)為正斷兼右旋走滑。兩次地震震源機制類型基本一致。1970年以來震中100km范圍內(nèi)ML5.0以上地震只有2006年7月4日文安ML5.5地震（據(jù)全國地震小震目錄）。文安地震震源機制為：P軸方位NEE、T軸方位NNW，接近水平。NNE向的節(jié)面屬于直立（傾角87°）；NWW向的節(jié)面傾斜（傾角43°）。震源斷層面的滑動角－133°，為正斷層兼走滑斷層類型［12］。易縣地震與文安地震都為正斷層兼走滑斷層類型，震源機制解近似。張紅艷等將張渤帶分為六個研究區(qū)，其中易縣地震震中位于張渤帶中南段，易縣地震的P軸方位為SWW、T軸方位為NNE，與張渤帶中南段的平均構(gòu)造應(yīng)力場趨于一致［13－14］。

4.2 同震形變響應(yīng)臺站震源球面投影分析

震中350km內(nèi)的主要形變臺站共計24個（包括河北省12個、山西省5個、天津市1個、北京市6個），震中距為0～100km的臺站3個，101～200 km的臺站8個，201～300km的臺站9個，301～350km的臺站4個。計算出震中至形變臺站的方位角、離源角，并根據(jù)重新定位結(jié)果計算出的震源機制解，對上述形變臺站做震源球面投影（圖5）。投影結(jié)果表明有同震響應(yīng)和無同震響應(yīng)的形變臺站主要集中在壓縮區(qū)域，由于易縣地震震級偏小，引起同震響應(yīng)的臺站相對較少，并且有同震響應(yīng)的臺站也主要集中在垂直擺傾斜儀測項，其中易縣垂直擺傾斜儀的同震響應(yīng)為W傾、S傾（圖3），與斷層錯動造成的地傾斜方向基本一致（圖5）。

圖5 形變臺站震源球面投影Fig.5 Seismic source projections of deformation stations.

5 結(jié)論

（1）易縣ML3.8、3.3地震屬于中小地震，但仍有4個臺站的13個測項記錄到較為明顯的同震形變響應(yīng)，這在中小地震的同震形變響應(yīng)分析中是較為突出的震例。由于近震及中小地震同震形變響應(yīng)的主要表現(xiàn)形式為同震階變和突跳，不同于遠震大震表現(xiàn)的同震形變波，不能通過面波延遲、最大面波幅度等進行分析。

（2）地震波形數(shù)據(jù)比定點形變數(shù)據(jù)的采樣頻率高數(shù)千倍，地震計記錄到的易縣地震的地震波持續(xù)時間僅有1～3分鐘，因此無法從同震形變響應(yīng)數(shù)據(jù)中進行震相分析，但同震形變響應(yīng)的持續(xù)時間與地震計記錄到的地震波持續(xù)時間基本一致。

（3）通過河北區(qū)域臺網(wǎng)的21個臺站的數(shù)字波形資料對易縣地震重新定位，河北臺網(wǎng)速報目錄經(jīng)緯度與重新定位經(jīng)緯度最大偏差分別為0.04°、0.03°。

（4）兩次地震震源機制類型基本一致，發(fā)震斷層表現(xiàn)為正斷兼右旋走滑，與張渤帶中南段平均區(qū)域構(gòu)造應(yīng)力場趨于一致。

（5）有同震響應(yīng)和無同震響應(yīng)的形變臺站主要集中在震源球面投影的壓縮區(qū)域，易縣垂直擺傾斜儀的同震響應(yīng)為 W傾、S傾，與斷層錯動造成的地傾斜方向基本一致。

致謝：本文得到天津市地震局、山西省地震局、中國地震臺網(wǎng)中心等單位提供的部分定點形變觀測資料，在此表示衷心感謝！

［1］牛安福，吉平，高福旺，等.印尼強地震引起的同震形變波［J］.地震，2006，26（1）：131－137.

［2］牛安福，張晶，吉平.強地震引起的同震形變響應(yīng)［J］.內(nèi)陸地震，2005，19（1）：1－7

［3］周克昌，吳翼麟.由同震傾斜階躍應(yīng)用反演方法研究震源機制［J］.地殼形變與地震，1990，10（1）：13－20.

［4］劉序儼.姑咱臺地面周日傾斜響應(yīng)特征分析［J］.中國地震，1993，9（1）：68－73.

［5］劉序儼，洪星.潮汐傾斜的幾何物理含義及其剪應(yīng)變屬性［J］.大地測量與地球動力學(xué)，2007，27（4）：82－86.

［6］李蘭生，劉建華.同震階變異常的特征與機制分析［J］.地殼形變與地震.1992，12（3）：81－87.

［7］國家地震局科技監(jiān)測司.地震形變前兆特征的識別與研究［M］.北京：地震出版社，1994.

［8］沈旭章，常千軍，梅秀蘋.蘭州形變臺水管傾斜儀映震能力的分析研究［J］.西北地震學(xué)報，2004，26（4）：368－370.

［9］胡衛(wèi)建，張俊山，周小川，等.隨州—信陽ML4.2地震前定點形變異常特征［J］.西北地震學(xué)報，2001，23（3）：274－278.

［10］梁尚鴻，李幼銘，束沛鎰，等.利用區(qū)域地震臺網(wǎng)P、S振幅比資料測定小震震源參數(shù)［J］.地球物理學(xué)報，1984，27（3）：247－257.

［11］刁桂苓，高國英，聶小紅，等.伽師序列中小地震震源機制分析［J］.國際地震動態(tài)，2007（12）：1－9.

［12］刁桂苓，張四昌，孫佩卿，等.2006年7月4日文安5.1級地震［J］.地震地質(zhì)，2006，28（3）：497－502.

［13］張紅艷，謝富仁，崔效鋒，等.張渤帶陸地段現(xiàn)代構(gòu)造應(yīng)力場的非均勻特征［J］.中國地震，2009，25（3）：314－324.

［14］方穎，張晶.張家口－渤海斷裂帶分段活動性研究［J］.地震，2009，29（3）：136－140.