李冬圣 蔡玲玲 李小軍 常 亮楊 銳 李雪英 賈 炯

（中國(guó)石家莊 050021河北省地震局）

河北贊皇小震群震源斷層分析

李冬圣 蔡玲玲 李小軍 常 亮楊 銳 李雪英 賈 炯

（中國(guó)石家莊 050021河北省地震局）

基于河北數(shù)字地震臺(tái)網(wǎng)和流動(dòng)臺(tái)站的觀(guān)測(cè)數(shù)據(jù)，利用HYPODD方法對(duì)2013年5月20日河北贊皇小震群進(jìn)行重新定位，定位結(jié)果顯示，地震呈NE向展布，震群震源深度集中在4—8km，采用模擬退火算法和高斯—牛頓算法相結(jié)合的有效算法進(jìn)行斷層參數(shù)反演計(jì)算，計(jì)算結(jié)果與震源機(jī)制解斷層面結(jié)果一致性較好，均為右旋走滑類(lèi)型。

贊皇震群；重新定位；震源機(jī)制；斷層面解

0 引言

根據(jù)河北省遙測(cè)地震臺(tái)網(wǎng)監(jiān)測(cè)，2013年5月20日河北贊皇發(fā)生小震群活動(dòng)，持續(xù)至7月3日，共記錄地震事件220次，其中單臺(tái)記錄地震49次，最大震級(jí)ML3.2，ML2.0以上地震9次。歷史地震資料表明，贊皇縣從未發(fā)生類(lèi)似震群活動(dòng)及中強(qiáng)以上地震。 近年來(lái)，隨著“首都圈示范工程項(xiàng)目”和“‘十五’網(wǎng)絡(luò)化”項(xiàng)目的建設(shè)，河北省地震監(jiān)測(cè)能力不斷提高，記錄到大量小震，且越來(lái)越精確的震相拾取使得小震定位精度越來(lái)越高?；谛≌鹂臻g分布對(duì)地質(zhì)構(gòu)造及地球物理相關(guān)問(wèn)題的研究引起越來(lái)越多地球科學(xué)家關(guān)注（楊智嫻等，2003；朱艾斕等，2005），萬(wàn)永革等（2008）首先提出利用小震分布和區(qū)域應(yīng)力場(chǎng)確定大震斷層面參數(shù)的方法，用于分析唐山地震序列，該方法不僅能在全局搜索最優(yōu)解，而且可以估計(jì)參數(shù)誤差，在求得大震斷層面后還可以根據(jù)區(qū)域構(gòu)造應(yīng)力場(chǎng)參數(shù)估計(jì)斷層面滑動(dòng)角。曾憲偉等（2013）利用該方法研究臨潭—宕昌斷層面參數(shù)，李守勇等（2011）用以確定磁縣1830年7.5級(jí)強(qiáng)震斷層面參數(shù)。本文試圖通過(guò)對(duì)贊皇震群的重新定位、斷層面擬合和震源機(jī)制解的分析，研究2013年贊皇震群的構(gòu)造背景。

1 研究方法

1.1 雙差定位方法

雙差定位法（HYPODD）在相對(duì)定位方法中應(yīng)用較廣，通過(guò)對(duì)一定空間范圍內(nèi)地震事件兩兩組對(duì)，根據(jù)組對(duì)事件之間震相到時(shí)差得到事件對(duì)的相對(duì)空間分布，在一定程度上消除了地殼速度結(jié)構(gòu)橫向不均勻帶來(lái)的定位誤差，進(jìn)而在計(jì)算中消除路徑異常（Waldhauser，Ellsworth，2000；楊智嫻等，2003）。雙差法基本方程為

其中，Δtik和Δtjk分別對(duì)應(yīng)地震i 和j 到臺(tái)站k 的到時(shí)（走時(shí)）殘差，Xi和Xj是兩個(gè)震源位置矢量，Sk是兩個(gè)震源連線(xiàn)中點(diǎn)到臺(tái)站k地震射線(xiàn)的慢度矢量，εik和εjk為地震在臺(tái)站的到時(shí)拾取誤差。

1.2 斷層面擬合方法

采用模擬退火全局搜索和高斯牛頓局部搜索相結(jié)合方法（萬(wàn)永革等，2008），利用贊皇震群重新定位結(jié)果，求解斷層面走向、傾角、位置及誤差估計(jì)，應(yīng)用華北地殼應(yīng)力場(chǎng)參數(shù)，計(jì)算斷層面滑動(dòng)角大小。

1.3 振幅比計(jì)算震源機(jī)制解

梁尚鴻等（1984）提出一種測(cè)定小震震源參數(shù)的方法，利用區(qū)域地震臺(tái)網(wǎng)垂直向位移記錄中Pg波和Sg波最大振幅比值，由計(jì)算層狀介質(zhì)中點(diǎn)源位錯(cuò)的Pg、Sg綜合地震圖測(cè)定小震震源參數(shù)。利用P波和SV波或SH的振幅比，是利用其輻射花樣比值求解震源機(jī)制參數(shù)。振幅比的輻射花樣隨空間方位的變化比單種波的輻射花樣強(qiáng)烈，力偶平面內(nèi)振幅比的輻射花樣呈8瓣分布。由此，只要有正確的直達(dá)波觀(guān)測(cè)振幅比，且觀(guān)測(cè)值歸算到震源球面位置正確，振幅比觀(guān)測(cè)對(duì)震源機(jī)制參數(shù)有較強(qiáng)的約束能力。

使用區(qū)域地震臺(tái)網(wǎng)垂直向地震AP和AS振幅資料計(jì)算震源機(jī)制軟件，配置臺(tái)站位置文件WDSTA、速度結(jié)構(gòu)文件WDHV；讀取事件波形初動(dòng)和振幅值寫(xiě)入WDF文件，進(jìn)行正演和反演計(jì)算，根據(jù)正演和反演輸出的結(jié)果文件，采用GMT繪圖選擇最佳解。

2 資料處理與分析

2.1 河北省遙測(cè)地震臺(tái)網(wǎng)觀(guān)測(cè)資料分析

2013年5月20日河北省遙測(cè)地震臺(tái)網(wǎng)監(jiān)測(cè)到贊皇地區(qū)發(fā)生小震群事件，當(dāng)日記錄地震事件11次，5月24至26日平均每天監(jiān)測(cè)小震達(dá)20次以上，截至7月3日共記錄地震220次，單日記錄事件最多達(dá)24次，最大震級(jí)為6月3日ML3.2（圖1）。參與贊皇地震定位臺(tái)站主要有河北贊皇臺(tái)（HEZAH）、河北元氏臺(tái)（HEYUS）、河北臨城臺(tái)（HELIC），外圍地震臺(tái)山西昔陽(yáng)臺(tái)（SXXIY）。為了進(jìn)一步跟蹤震情，河北省地震監(jiān)測(cè)網(wǎng)絡(luò)中心在贊皇縣以西架設(shè)兩套流動(dòng)臺(tái)（L1301，L1303）。由此，地震監(jiān)測(cè)臺(tái)站均勻分布在贊皇四周，形成較好的地震監(jiān)測(cè)條件，ML0地震至少可被3—4個(gè)臺(tái)站記錄。

以河北省遙測(cè)臺(tái)網(wǎng)產(chǎn)出地震目錄繪圖顯示，贊皇地震呈EW向分布，大部分地震分布在贊皇臺(tái)以西，只有零星幾個(gè)地震分布在贊皇臺(tái)以東位置，見(jiàn)圖2（a）。觀(guān)測(cè)報(bào)告顯示，贊皇臺(tái)以東地震震級(jí)偏小，且發(fā)震時(shí)間在架設(shè)流動(dòng)臺(tái)站前，參與定位臺(tái)站少、張角大，尤其使用hypo 2000和單純形法分別進(jìn)行定位，定位結(jié)果不穩(wěn)定。

圖2 精定位前后震中分布（a）精定位前震中分布；（b）精定位后震中分布Fig.2 Earthquake distribution before and after the accurate location

2.2 贊皇震群重新定位

使用HYPODD方法（Waldhauser et al，2000）對(duì)贊皇震群進(jìn)行重新定位。截至7月3日，贊皇震群共觀(guān)測(cè)到220次地震，剔除單臺(tái)記錄的49次事件，選取171次地震參與重新定位。震相數(shù)據(jù)包含8 085個(gè)P波到時(shí)、8 635個(gè)S波到時(shí)，重新定位前對(duì)震相數(shù)據(jù)進(jìn)行檢查，修訂定位結(jié)果不穩(wěn)定的地震參數(shù)，最終平均每次地震具有8個(gè)震相數(shù)據(jù)。

重新定位使用的分層速度模型參考邢臺(tái)震源區(qū)地殼上地幔速度結(jié)構(gòu)（祝治平等，1995），模型分6層（表1），波速比設(shè)定為1.72。原始定位采用單純型法，無(wú)水平和垂直誤差，擬合殘差平均值0.15 s。精定位得到147個(gè)震源位置參數(shù)，均方根殘差平均值為0.05 s，震源位置南北向平均誤差0.10km，東西向平均誤差0.08km，深度誤差平均 0.15km。

表1 贊皇地區(qū)地殼速度模型Table 1 Velocity model of the crust in ZAH area

贊皇震群精定位前，地震分布呈EW向，且主要集中在兩個(gè)區(qū)域內(nèi)，精定位后最大震級(jí)ML3.2的地震發(fā)生在震群中間部位，地震分布比較集中，呈NE向展布，震群震源深度平均為5.8km，依照?qǐng)D2（b）對(duì)小震群序列按NW和NE方向進(jìn)行剖面劃分，并將地震震源深度投影到AA′和BB′兩個(gè)剖面，見(jiàn)圖3，深度剖面顯示，震群在震源深度變化上具有一定規(guī)律。圖3（a）為地震條帶沿線(xiàn)AA′的震源擴(kuò)展剖面，可見(jiàn)贊皇小震群總體震源深度較淺，集中分布在4—8km范圍內(nèi)；地震條帶總長(zhǎng)約4km，無(wú)明顯空間幾何特征。圖3（b）為地震條帶沿直線(xiàn)BB′的震源擴(kuò)展剖面，震源深度沿BB′剖面總體傾角較陡，可見(jiàn)震源分布隨距斷裂遠(yuǎn)近存在從淺到深的分布，勾畫(huà)震源深度剖面包絡(luò)線(xiàn)，可以得到傾角約70°的斷面。

圖3 重新定位后震中分布剖面（a）AA′ 剖面；（b）BB′ 剖面Fig.3 Earthquake epicenters after relocation and profle of focal depth

2.3 斷層面求解

采用贊皇震群重新定位數(shù)據(jù)反演斷層參數(shù)，結(jié)果顯示，斷層走向NE，傾向NE［圖4（a），圖4（b）］，傾角較大［圖4（c）］，小震距斷層面在0.5km以?xún)?nèi)精度較高［圖4（d）］，根據(jù)萬(wàn)永革等（2007）提出的將模擬退火算法和高斯—牛頓算法相結(jié)合的有效算法，避免對(duì)初始解的依賴(lài)得到全局最優(yōu)解，求得斷層走向、傾角、到坐標(biāo)原點(diǎn)距離及斷層面4個(gè)頂點(diǎn)位置（表2），并給出相應(yīng)標(biāo)準(zhǔn)差。由贊皇震群反演確定斷層近乎直立，北偏東41°［圖2（b）］。應(yīng)用萬(wàn)永革等（2007）提出的應(yīng)力軸方位較傾角準(zhǔn)確、R值相對(duì)難以確定的準(zhǔn)則，假定P、T軸方位誤差5°，相對(duì)應(yīng)力大小誤差0.1°，依據(jù)大華北構(gòu)造應(yīng)力場(chǎng)參數(shù)（Wan Y G，2010）和地震斷層面幾何參數(shù)及標(biāo)準(zhǔn)差，推算滑動(dòng)角為-155.907 5°，認(rèn)為斷層為右旋走滑斷層，標(biāo)準(zhǔn)差為6.617 4°。

圖4 贊皇地震斷層重新定位小震分布（a）水平面；（b）斷層面；（c）垂直于斷層面的橫斷面投影；（d）小震距斷層面距離的分布Fig.4 Distribution of precisely located small earthquakes near the ZAH fault

表2 運(yùn)用贊皇震群精定位小震資料求得的地震斷層面走向、傾角和位置Table 2 Fault plane parameters determined precisely local earthquake for segments in ZAH earthquake

2.4 地震序列震源機(jī)制解

在獲得精確的震源位置后，采用精定位所用速度結(jié)構(gòu)和垂直向直達(dá)P、S波最大振幅比及初動(dòng)計(jì)算震源機(jī)制解的方法（梁尚鴻等，1984），對(duì)贊皇震群ML＞ 2.0地震進(jìn)行震源機(jī)制解研究。為了保證單個(gè)震源機(jī)制解的穩(wěn)定性，參與計(jì)算的臺(tái)站數(shù)目要求不少于6個(gè)，且在震中周?chē)鶆蚍植迹鹬芯嗫刂圃?20km之內(nèi)，避免出現(xiàn)Pn和Pg混淆現(xiàn)象；2013年5月20日至7月3日發(fā)生ML＞ 2.0地震9次，使用位于震中四象限范圍內(nèi)數(shù)字記錄垂直分量P波初動(dòng)、P波和S波最大振幅比法，反演震源機(jī)制解。5月26日22點(diǎn)53分ML2.1地震由于雙震重疊無(wú)法識(shí)別P波初動(dòng)，其余8個(gè)震源機(jī)制類(lèi)型基本一致，均為右旋走滑（圖5）。

圖5 贊皇震群震源機(jī)制解Fig.5 Focal mechanism solutions of ZAH earthquake swarm

由圖5可見(jiàn)，參與計(jì)算震源機(jī)制的臺(tái)站空間分布較好，P波初動(dòng)符號(hào)（向上或向下）呈很好的四象限分布，且觀(guān)測(cè)點(diǎn)沿地震射線(xiàn)在震源球下半球的投影點(diǎn)及初動(dòng)符號(hào)對(duì)2個(gè)節(jié)面的空間位置限制嚴(yán)格。因此，認(rèn)為該震源機(jī)制解比較可靠。根據(jù)余震分布［圖2（b）］確定NE向節(jié)面為斷層面，角度約40°。贊皇臺(tái)、臨城臺(tái)、昔陽(yáng)臺(tái)和井陘臺(tái)及后架設(shè)的兩套流動(dòng)臺(tái)站（L1301，L1303），均記錄到震群中ML＞ 1.0地震，且地震的P波初動(dòng)符號(hào)高度一致，贊皇臺(tái)為“-”，臨城臺(tái)為“+”，L1301為“-”，L1301為“+”。

圖6是震源機(jī)制解參數(shù)統(tǒng)計(jì)，Strike為節(jié)面走向，主要集中于NE45°和NW45°；dip是節(jié)面傾角，集中在75°，和震源斷層擬合結(jié)果一致性較好；slip是節(jié)面滑動(dòng)角，接近0°或180°，反映出震源斷層多為走滑類(lèi)型，與斷層擬合結(jié)果的走滑性質(zhì)基本一致。震源機(jī)制應(yīng)力軸P、T軸傾角Pdip和Tdip均較小，接近水平方向，反映了斷層的在水平主壓應(yīng)力作用下的走滑錯(cuò)動(dòng)狀況，Paz是最大壓應(yīng)力主軸的方位，NEE方向優(yōu)勢(shì)分布明顯，Taz是最小壓應(yīng)力主軸的方位，NNW方向優(yōu)勢(shì)分布明顯和華北地區(qū)區(qū)域應(yīng)力場(chǎng)方向基本一致。

圖6 震源機(jī)制解參數(shù)統(tǒng)計(jì)Fig.6 Parameters of focal mechanism solutions

3 結(jié)論

（1）對(duì)贊皇震群進(jìn)行重新定位，結(jié)果顯示：地震分布比較集中，且呈NE向展布，主要集中在寬4km、深4—8km的空間范圍內(nèi)；采用模擬退火全局搜索和高斯—牛頓局部搜素相結(jié)合的方法，估算斷層面走向NE41.28°±2.80°、傾角76.97°±1.02°，基于華北地殼應(yīng)力場(chǎng)參數(shù)，反演斷層面滑動(dòng)角為-155.9075°±6.6174°。

（2）計(jì)算震群中8個(gè)ML2.0以上地震的震源機(jī)制解，斷層面優(yōu)勢(shì)走向?yàn)镹E45°和NW45°。傾角約75°，P、T軸接近水平方向，和華北區(qū)域應(yīng)力場(chǎng)基本一致。

（3）小震精確定位結(jié)果模擬斷層面與震源機(jī)制解斷層面結(jié)果一致性較好，表明此次小震群活動(dòng)在區(qū)域應(yīng)力場(chǎng)作用下可能形成一系列性質(zhì)接近的微小斷裂。

小震群活動(dòng)往往發(fā)生在多震地區(qū)及有大斷層穿過(guò)的地震活動(dòng)帶上，贊皇震群區(qū)屬于弱震區(qū)，周邊30km內(nèi)歷史上未發(fā)生過(guò)5級(jí)以上地震（楊銳等，2014），小震活動(dòng)弱，且震源區(qū)附近無(wú)較大斷層，在此背景下，2013年贊皇震群活動(dòng)顯得更加突出，發(fā)生機(jī)理較為特殊。因附近無(wú)明顯斷裂穿過(guò)，而此次震群的最大震級(jí)僅為ML3.2，即使形成斷層，尺度也較小，如何研究和描述需開(kāi)展進(jìn)一步工作。小震斷層擬合與震源斷層面解相互印證的研究思路需要更多震例檢驗(yàn)，在此希望對(duì)相關(guān)地震研究提供新思路和方法。

陳颙.用震源機(jī)制一致性作為描述地震活動(dòng)性的新參數(shù)［J］.地球物理學(xué)報(bào)，1978，21（02）：142-159.

刁桂苓，王俊國(guó)，張宏志.震源機(jī)制一致性的預(yù)測(cè)意義［J］.國(guó)際地震動(dòng)態(tài)，2004，z1：35.

房立華，吳建平，等.2011年云南盈江MS5.8地震及其余震序列重定位［J］.地震學(xué)報(bào)，2011，33（2）：262-267.

李守勇，張雙鳳，閆俊崗.利用小震分布和區(qū)域應(yīng)力場(chǎng)確定磁縣1830年7.5級(jí)強(qiáng)震斷層面參數(shù)［J］.地震地磁觀(guān)測(cè)與研究，2011，32（3）：20-25.

劉伯禮，蔡華昌，趙軍.石家莊斷裂活動(dòng)性及其對(duì)市區(qū)穩(wěn)定性的影響［J］.華北地震科學(xué)，1991，9（4）：13-23.

劉科，王景明.河北平原構(gòu)造特征、地裂縫分布及成因機(jī)制［J］.南水北調(diào)與水利科技，2005，3（6）：38-42.

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

苗春蘭，王曉山，等.2008年喀喇沁旗震群精定位與震源機(jī)制分析［J］.地震地磁觀(guān)測(cè)與研究，2013，34（1/2）：10-13.

萬(wàn)永革，沈正康，刁桂苓，等.利用小震分布和區(qū)域應(yīng)力場(chǎng)確定大震斷層面參數(shù)方法及其在唐山地震序列中的應(yīng)用［J］.地球物理學(xué)報(bào)，2008，51（3）：793-804.

王玉婷，呂悅軍，謝卓娟，等.河北平原地震帶的現(xiàn)金活動(dòng)性分析［J］.地震地磁觀(guān)測(cè)與研究，2012，33（2）：6-14.

楊銳，王曉山，王想，等.2013年河北贊皇小震群活動(dòng)特征及其指示意義［J］.華北地震科學(xué)，2014，32（4）：17-20.

楊智嫻，陳運(yùn)泰，鄭月軍，等.雙差地震定位法在我國(guó)中西部地區(qū)地震精確定位中的應(yīng)用［J］.中國(guó)科學(xué)，2003，33（增刊1）：129-134.

于湘?zhèn)?，陳運(yùn)泰，等.華北地區(qū)地震重新定位結(jié)果分析［J］.大地測(cè)量與地球動(dòng)力學(xué)，2010，30（2）：29-33.

朱艾斕，徐錫偉，等.川西地區(qū)小震重新定位及其活動(dòng)構(gòu)造意義［J］.地球物理學(xué)報(bào)，2005，48（3）：629-636.

祝治平，張先康，蓋玉杰，等.邢臺(tái)震源區(qū)及相鄰地區(qū)地殼上地慢速度結(jié)構(gòu)研究［J］.地震學(xué)報(bào)，1995，17（3）：328-334.

曾憲偉，馮建剛，王曉山，等.利用小震分析和區(qū)域應(yīng)力場(chǎng)確定臨潭—宕昌斷裂斷層面參數(shù)［J］.地震工程學(xué)報(bào)，2013，35（4）：743-750.

Waldhauser F，Ellsworth W L.A double-difference earthquake location algorithm：Method and application to the Northern Hayward Fault，California［J］.Bulletin of the Seismological Society of America，2000，90：1 353-1 368.

WanY G.Contemporary tectonic stress feld in China［J］.Earthq Sci，2010，23：377-386.

Analysis of the focal tomography for Zanhuang small earthquake swarm in Hebei Province

Li Dongsheng，Cai Lingling，Li Xiaojun，Chang Liang，Yang Rui，Li Xueying and Jia Jiong
（Earthquake Administration of Hebei Province，Shijiazhuang 050021，China）

Based on the data from Hebei Digital Seismic Network stations and mobile stations，we relocated the earthquakes of Zanhuang small earthquake swarm which occurred on May 20th，2013 through HYPODD methods.The results show that the distribution of those earthquakes is along NE and the source depths are mainly between 4km to 8km.The result from the method combined simulated annealing algorithm and Gauss - Newton algorithm shows that the focal tomography is strike 41°，dip 77°，rake -156°.The focal mechanism for eight earthquakes with magnitude higher than ML2.0 in this earthquake swarm through inversion all show strike NE and NW tendency nodal plane which are consistent with the focal plane with dextral strike-slip fault type.

Zanhuang earthquake swarm，HYPODD relocation，focal mechanism，fault plane solution

10.3969/j.issn.1003-3246.2015.05.004

李冬圣（1982—），女，工程師，主要從事地震監(jiān)測(cè)及分析工作。E-mail：Lds981112@163.com

2014年測(cè)震臺(tái)網(wǎng)青年骨干項(xiàng)目（20140302）及河北省重點(diǎn)地區(qū)殼幔結(jié)構(gòu)及地震監(jiān)測(cè)預(yù)報(bào)關(guān)鍵技術(shù)研究（132754071）聯(lián)合資助

本文收到日期：2015-01-08