邯鄲市周邊斷裂帶的活動特征分析

周俊杰 宋宏利 劉海新 張紅芬 周龍泉

1)河北工程大學，邯鄲 056038

2)河北省資源勘測重點實驗室，邯鄲 056038

3)中國地震臺網中心，北京 100036

邯鄲市周邊斷裂帶的活動特征分析

周俊杰1，2)宋宏利1)劉海新1)張紅芬1)周龍泉3)

1)河北工程大學，邯鄲 056038

2)河北省資源勘測重點實驗室，邯鄲 056038

3)中國地震臺網中心，北京 100036

斷裂活動和地震活動具有互為因果的關系，而不同于一般做法，文中用地震的活動性來分析斷裂帶的活動特征。為了解邯鄲市周邊存在的斷裂帶的活動特征，基于邯鄲市附近地震臺網監測的地震數據，利用地震縱波速度和震源位置聯合反演方法，得到了邯鄲市周邊主要斷裂帶上地震波速度及震中在空間和時間上的分布特征，確定了地下殼質的縱波速度結構特征。綜合分析表明:邯鄲市周邊小震的活動主要由磁縣斷裂和林州-武安斷裂引起，太行山山前斷裂帶南端(邯鄲—磁縣段)幾乎沒有活動;地震活動集中在6～10km的范圍內;活動時間特征是，從開始的2條斷裂帶同時活動逐漸演變為以磁縣斷裂活動為主。

斷裂帶 地震分布 邯鄲市 速度反演

0 引言

地震常發生在活動斷裂上，大量震例表明，地震斷層不僅是產生地震的根源，而且地震時沿斷層廊帶的破壞最為嚴重，人員傷亡也明顯大于斷層兩側的其他地區。近百年來，因地震造成的80%的財產損失、60%的人員傷亡發生在有活斷層的城市中(李樹德，2001;宋和平，2005;王斌，2008)。要減少各方面的損失，一方面要加強建筑物抗震設計并確保建筑質量，另一方面是應迅速查清城市地下隱伏活斷層的分布、性質、活動時代等，為合理地進行城市規劃和社會發展的科學決策提供依據。

據統計，約有95%的6級以上地震分布在新生代以來活動斷裂(帶)以及構造盆地、大型塊體隆起和凹陷邊界活動斷裂帶上，震級越大，這種關系越明顯(李樹德，2001;宋和平，2005)。地震發生在地殼深部，其孕育和發生與震源附近的地殼結構、物質組成和孕震構造等直接相關。依據斷裂區內的地震臺網資料，綜合利用地震波速度成像和構造成像方法，揭示不同深度上的地殼精細結構、深淺構造特征及其與地震分布之間的關系(劉福田，1989;劉建華等，1989;Zhao，2002;郭飚等，2006;胡幸平等，2008)，能更加精確地預測該斷裂構造的活動性。

邯鄲市位于河北省的南部，周邊存在近SN向的太行山斷裂、近EW向的磁縣斷裂和永年斷裂等一系列的斷裂構造，為太行山次級塊體和冀魯次級塊體的結合處，地震活動較為頻繁。研究者認為，1830年6月12日磁縣7.5級強震是近EW向的磁縣斷裂活動導致的。近年來的地震臺網監測數據顯示，在斷裂附近發生了不同級別的中小地震。利用地震波縱波速度和震源位置聯合反演方法，反演出不同地層的結構和震源參數，分析不同位置、不同深度上地震波速度和震源深度的變化，研究從空間和時間上斷裂構造的活動性與地震分布之間的規律，從而得到地震震中分布隨空間和時間的變化規律，預測邯鄲市周邊活動斷裂構造的特征。

1 地震波速度和震源位置聯合反演方法簡介

為了更精準地確定小地震的位置和活動特征，我們采用了地震速度結構和震源位置聯合反演方法(Aki et al.，1976;劉福田，1989;周龍泉等，2005，2006)，利用劉福田提出的正交投影算子，依據地震臺網監測的地震事件，確定研究區的速度結構參數和震源參數。

反演中用到的是3個以上地震臺站記錄到震相且有震源參數(經度、緯度和深度)的地震事件，1992—2000年為1072個，2001—2007年為924個，這些地震均屬于Ⅰ類精度。

所有地震波走時資料都取自河北省地震臺網觀測報告，考慮到全部地震射線對研究區的覆蓋以及較大走時殘差對反演的影響，從1996個地震事件共7898條射線中，進一步挑選出參與反演的P波觀測數4825條，參與反演的走時與震中距關系如圖1。

圖1 參與反演的P波走時與震中距關系圖Fig.1 The relationship between the seismic Pwave travel-time and the epicentral distance in the inversion.

綜合多年來發表的該區及其鄰近區域的地殼速度結構研究成果(Huang et al.，2002;張學民等，2006)，選定本區成像的一維參考速度模型如表1。

用網格方法對速度模型參數化，在平面方向上劃分0.25°×0.25°的網格，垂直方向上劃分如表1所示，在模型中速度結構用連續函數表示，網格內任意一點的速度用內插方式計算。

2 反演結果分析

2.1 解的分辨率分析

采用解的分辨矩陣的對角元素來說明反演結果的可靠性(圖2)。

圖2 不同深度分辨矩陣對角元素值等值線分布圖Fig.2 The distinguishability of diagonal element isogram in different depths.

由圖2可見，在0km深度上分辨矩陣對角元素值基本上在0.5左右，只有少數地區達到0.6以上，分辨效果一般;在5km深度上，磁縣斷裂附近分辨矩陣對角元素值基本在0.9以上，分辨效果非常好，在安陽和邯鄲附近也能達到0.7左右，只有邊緣地區分辨較差;在10km深度上，絕大多數地區的分辨矩陣對角元素值≥0.9，邊緣地區值也能達到0.7左右;在15km深度上，分辨矩陣對角元素值分布與10km深度類似，基本上在0.9左右;在20km深度上，研究區的中東部分辨矩陣對角元素值在0.9左右，而在西部對角元素值在0.5左右;在33km深度上，只有中部地區分辨矩陣對角元素值能達到0.7左右，而其他區在0.3左右。

綜合分析認為該地區的速度在5km、10km、15km和20km深度上具有非常好的分辨效果，在0km和33km深度上部分地區的分辨效果比較滿意。

2.2 定位結果分析

圖3是重新定位前(a)后(b)地震震中分布圖。可以看出，重新定位后，地震分布更為集中。

圖3 研究區定位前(a)、后(b)地震震中分布圖Fig.3 The epicenter distribution of earthquakes before(a)and after(b)relocation in the study area.

3 地震震中空間分布特征

3.1 地震震中平面分布特征

研究區小震震中主要分布于太行山山前斷裂帶西側(圖3b)，東側較少。而西側的震中又呈2個帶狀分布:一是沿林州-武安斷裂和太行山山前斷裂之間的太行山隆起區，呈NE向展布，位于剖面線D—D'和F—F'之間，震中集中分布在林州-武安斷裂帶兩側;二是沿近EW向的磁縣斷裂分布，位于磁縣斷裂帶兩側，且在這條帶中，磁縣斷裂被太行山山前斷裂分隔為東西兩段，西段小震分布集中，而東段基本上沒有小震。另外，太行山山前斷裂帶兩側小震分布非常少，其中沿邯鄲-磁縣斷裂幾乎沒有小震分布，反映該段處于平靜狀態;在太行山山前斷裂南段(磁縣斷裂之南)東側有一些小震活動。

3.2 震源深度分布特征

地震震源深度統計結果見圖4。研究區小震震源深度分布在0～40km之間，深度分布的峰值在10km左右，10km深度的地震占總數的18%，6～10km深度地震占總數的66%，4～22km深度地震占總數的93%。這些地震的震源深度分布表明，研究區的地震破裂集中在一定的優勢層位，可以將2km和22km分別作為優勢層位的震源深度上下界限。

3.3 主要斷裂帶上地震活動特征

為了分析主要斷裂帶上小地震活動的特征，沿其走向和垂直走向方向做了6條剖面 (圖3b)。由于該區域斷層傾角較大，所以只將垂直剖面兩側10km范圍內的地震投影到剖面上，判定這些地震為斷層活動引起的。A—A'、B—B'和C—C'近乎垂直于太行山山前斷裂和林州-武安斷裂的走向，與磁縣斷裂平行，且B—B'位于磁縣斷裂帶上。D—D'、E—E'和F—F'平行于太行山山前斷裂和林州-武安斷層，與磁縣斷裂幾乎垂直，其中D—D'位于林州-武安斷裂上，F—F'位于太行山山前斷裂帶上。

圖5為 A—A'、B—B'和 C—C'剖面線上速度分布圖以及剖面兩側10km以內小震沿剖面的投影分布情況。分析認為太行山山前斷裂顯示為向東傾的正斷層，林州-武安斷裂為向西傾的正斷層，認為這2個正斷層是由太行山隆起造成的。小震主要分布在太行山隆起區，且小震震源深度較深，大部分集中在6～20km范圍內，而太行山隆起區兩側小震分布少，震源深度較淺，深度集中在10km以內。沿磁縣斷裂(B—B'剖面)，小震主要分布在太行山山前斷裂與林州-武安斷裂之間，而且小震分布密集，這表明磁縣斷裂帶西段依然處于活動狀態，東段處于相對平靜狀態。

圖4 1992年以來地震震源深度分布頻次Fig.4 The frequency distribution of focal depths since 1992.

圖5 A—A'、B—B'和C—C'剖面上地震波速度和震中分布Fig.5 Distribution of seismic wave velocity and the epicenter on A-A'，B-B'and C-C'sections.

圖6 D—D'、E—E'和F—F'剖面地震波速度和震中分布Fig.6 Distribution of seismic wave velocity and the epicenter on D-D'，E-E'and F-F'sections.

圖7 磁縣斷裂帶上不同時間段小震垂直分布Fig.7 Distribution of small earthquakes at different times on the Cixian Fault zone in vertical direction.

圖6為 D—D'、E—E'和 F—F'3條垂直斷裂走向的剖面速度分布圖以及剖面兩側10km以內小震分布情況。結果顯示磁縣斷裂向南傾，且磁縣斷裂的南側地殼厚度較北側小，這與兩側的震源深度基本相符，據此推測太行山在磁縣斷裂的南側隆起并不明顯。在林州-武安斷裂(D—D'剖面)上，小震分布非常密集，這表明該斷裂帶目前處于較活躍狀態。沿太行山山前斷裂帶(F—F'剖面)，只有在磁縣斷裂帶附近小震非常活躍，而根據這些小震分布特征以及磁縣斷裂的傾向認為，這些小震主要是由磁縣斷裂活動造成的，受太行山山前斷裂的影響比較弱，表明太行山山前斷裂的南端(邯鄲—磁縣段)目前并沒有活動。

3.4 活動斷裂上地震的時間分布特征

經過小震的平面分布特征和不同方向上的剖面分布特征分析認為，邯鄲市周邊的小震活動主要是由林州-武安斷裂和磁縣斷裂的活動造成的，二者活動時間主要表現為2個階段，2000年前和2000年后。在這2個時間段上，磁縣斷裂帶和林州-武安斷裂帶附近小震垂直分布情況分別表示在圖7和圖8。

圖8 林州-武安斷裂帶上不同時間段小震垂直分布Fig.8 Distribution of small earthquakes at different times on the Linzhou-Wu'an Fault zone in vertical direction a 1992—1999年;b 2000—2007年;c速度等值線分布

圖7表明，在2000年以前磁縣小震主要分布在磁縣斷裂與太行山山前斷裂和林州-武安斷裂的交界附近，而2000年以來小震主要分布在磁縣斷裂與太行山山前斷裂以西，其交界附近基本沒有地震活動。圖8顯示，2000年以前林州-武安斷裂帶上小震都很活躍，而2000年以來小震只集中在與磁縣斷裂帶交界附近。

4 結論

(1)邯鄲市周邊的主要斷裂為向東傾的太行山山前斷裂、向西傾的林州-武安斷裂和向南傾的磁縣斷裂，這是引起小震活動的主要斷裂構造。

(2)邯鄲市周邊斷裂帶上的地震震中主要分布在太行山山前斷裂帶西側的太行山隆起區，平原帶地震較少，而太行山山前斷裂西側的小震震中又集中分布在林州-武安斷裂帶和磁縣斷裂帶的兩側，表明邯鄲市周邊的地震活動主要受這2個斷裂帶的控制。

(3)地震的震源深度分布表明地震破裂集中在一定的優勢層位，深度分布的峰值在10km左右，范圍是6～10km，可以將2km和22km分別作為優勢層位的震源深度上下界限。

(4)從小震活動情況看出，林州-武安斷裂在研究的時間段內都比較活動，而磁縣斷裂2000年前的活動主要集中在與林州-武安斷裂和太行山山前斷裂交界部位，2000年后逐漸演化為以磁縣斷裂帶的西段活動為主。

郭飚，劉啟元，陳九輝，等.2006.中國境內天山地殼上地幔結構的地震層析成像［J］.地球物理學報，49(6):1693—1700.

GUO Biao，LIU Qi-yuan，CHEN Jiu-hui，et al.2006.Seismic tomography of the crustand uppermantle structure underneath the Chinese Tianshan［J］.Chinese JGeophys，49(6):1693—1700(in Chinese).

胡幸平，俞春泉，陶開，等.2008.利用P波初動資料求解汶川地震及其強余震震源機制解［J］.地球物理學報，51(6):1711—1718.

HU Xing-ping，YU Chun-quan，TAO Kai，etal.2008.Focalmechanism solutions ofWenchuan earthquake and its strong aftershocks obtained from initial Pwave polarity analysis［J］.Chinese JGeophys，51(6):1711—1718(in Chinese).

李樹德.2001.試析地震與斷裂的關系［J］.北京大學學報(自然科學版)，37(3):394—399.

LIShu-de.2001.A brief study on the relationships between earthquakes and active faults［J］.Acta Scientiarum Naturalium Universitatis Pekinensis，37(3):394—399(in Chinese).

劉福田，李強，吳華，等.1989.用于速度圖像重建的層析成像法［J］.地球物理學報，32(1):46—61.

LIU Fu-tian，LIQiang，WU Hua，et al.1989.On the tomography inverse method used in velocity image reconstruction［J］.Chinese Journal of Geophysics，32(1):46—61(in Chinese).

劉建華，劉福田，吳華，等.1989.中國南北帶地殼和上地幔的三維速度結構［J］.地球物理學報，32(2):143—151.

LIU Jian-hua，LIU Fu-tian，WU Hua，etal.1989.The three-dimensional velocity structure of the crust and uppermantle beneath the North-South Belt of China［J］.Chinese Journal of Geophysics，32(1):46—61.

宋和平.2005.論新疆深大斷裂特征與地震的關系［J］.內陸地震，19(3):193—201.

SONG He-ping.2005.The relation between characteristics of deep and big faults and earthquakes in Xinjiang［J］.Inland Earthquake，19(3):193—201(in Chinese).

王斌，梁雪萍，周健.2008.江蘇及其周邊地區斷裂活動性與地震關系的分析［J］.高原地震，20(1):38—43.

WANG Bin，LIANG Xue-ping，ZHOU Jian.2008.Analysis of relationship between faultactivity and earthquakes in Jian-gsu Province and its adjacent areas［J］.Plateau Earthquake Research，20(1):38—43(in Chinese).

張學民，刁桂苓，趙英萍，等.2006.華北地區深部地幔剪切波速度結構研究［J］.地球物理學報，49(6):1709—1719.

ZHANG Xue-min，DIAOGui-ling，ZHAOYing-ping，etal.2006.Study onmantle shearwave velocity structures in North China［J］.Chinese JGeophys，49(6):1709—1719(in Chinese).

周龍泉，劉福田，劉勁松，等.2005.利用τ-p波場反演法確定東沙群島的地殼速度模型［J］.地球物理學進展，20(2):503—506.

ZHOU Long-quan，LIU Fu-tian，LIU Jin-song，etal.2005.Determination of the crustal velocitymodelof Dongsha Islands using the inversion ofτ-p wave field［J］.Progress in Geophysics，20(2):503—506(in Chinese).

周龍泉，劉福田，陳曉非.2006.三維介質中速度結構和界面的聯合成像［J］.地球物理學報，49(4):1062—1067.

ZHOU Long-quan，LIU Fu-tian，CHEN Xiao-fei.2006.Simultaneous tomography of 3-D velocity structure and interface［J］.Chinese JGeophysics，49(4):1062—1067(in Chinese).

Aki K，Lee W H K.1976.Determination of three-dimensional velocity anomalies under array using P arrival times from local earthquakes(1):A homogeneous initialmodel［J］.JGeophys Res，81:4381—4399.

Huang J，Zhao D，Zheng S.2002.Lithosperic structure and its relationship to seismic and volcanic activity in southwest China［J］.JGeophys Res，107(B 10):10，1029/2000 JB000137.

Zhao D，Mishra O，Sanda R.2002.Influence of fluids and magma on earthquakes:Seismological evidence ［J］.Phys Earth Planet Interi，102:249—267.

THE ACTIVE CHARACTERISTICSOF FAULTS AROUND HANDAN CITY

ZHOU Jun-jie1，2)SONG Hong-li1)LIU Hai-xin1)ZHANG Hong-fen1)ZHOU Long-quan3)
1)Hebei University of Engineering，Handan 056038，China
2)Key Laboratory of Resource Survey and Research of Hebei Province，Handan 056038，China
3)China Earthquake Networks Center，Beijing 100036，China

Activity of fault zones can often lead to earthquake.On the contrary，earthquake can also help people to analyze and forecast the characteristics of faults.To know more about the characteristics of active faults around Handan，and serve better for the urban construction and planning，a joint inversion method，which is based on earthquake data from local seismic network in Handan，was employed in this paper to get the distribution characteristic of earthquake parameters in space and time on themain faults around Handan and determine the characteristics of P-wave velocity structure in the crustalmedium.Integrated analysis shows:The small earthquake activities around Handan are caused mainly by Cixian Fault and the Linzhou-Wu'an Fault;There is almost no activity on the southern end of Taihang Mountains piedmont fault zone(Cixian-Wu'an section);Earthquake activities are concentrated in a range of6～10km;Faultactivity gradually evolved from the faulting of the two faults simultaneously at the beginning to that the Cixian Faultwas dominant.

fault zone，earthquake distribution，Handan City，velocity inversion

P315.2

A

0253-4967(2012)01-0100-10

10.3969/j.issn.0253-4967.2012.01.010

2011-04-05收稿，2011-05-16改回。

國家自然科學基金(40872056)“大別山西段高壓-超高壓榴輝巖變質相平衡關系研究”項目資助。

周俊杰，男，1976年生，2000年畢業于中國地質大學(武漢)，2006年在河北工程大學獲得碩士學位，講師，主要從事地球物理探測和構造地質方面的教學與研究，電話:0310-8579542，E-mail:ansanzhou7684@163.com。