陳 睿閆俊崗閆 楷萬(wàn) 娜張予川

1）中國(guó)湖北430074中國(guó)地質(zhì)大學(xué)（武漢）

2）中國(guó)鄭州450016河南省地震局

3）中國(guó)河北056006邯鄲中心地震臺(tái)

利用接收函數(shù)研究太行山前斷裂南部地區(qū)地殼厚度

陳 睿1），2）閆俊崗3）閆 楷2）萬(wàn) 娜2）張予川2）

1）中國(guó)湖北430074中國(guó)地質(zhì)大學(xué)（武漢）

2）中國(guó)鄭州450016河南省地震局

3）中國(guó)河北056006邯鄲中心地震臺(tái)

利用太行山前斷裂南部地區(qū)寬頻帶數(shù)字化地震臺(tái)站的遠(yuǎn)震體波記錄，采用頻率域反褶積方法提取接收函數(shù)，由H—Kappa疊加方法反演得到各臺(tái)站下方地殼厚度和泊松比。結(jié)果顯示： ①該區(qū)域地殼厚度分布為西部較厚，東部較薄，莫霍面呈小角度向西傾斜。西部長(zhǎng)治盆地地殼厚度40.2km，東部平原地區(qū)地殼較薄，厚度31.8km，地殼變化落差約9km；②在太行山與平原交接地帶，特別是沿太行山弧形分布的主峰外緣，即一級(jí)夷平面和二級(jí)夷平面的分界線附近地殼厚度變化較快；③接收函數(shù)反演的地殼厚度與人工地震測(cè)深結(jié)果吻合，太行山隆起區(qū)地殼變厚，湯陰地塹地殼變薄，內(nèi)黃隆起地殼變厚，表明深部莫霍面或上地幔的隆起、凹陷與地面起伏呈鏡像對(duì)稱；④區(qū)域內(nèi)泊松比值變化不大，約0.23—0.27，表明該地區(qū)地下介質(zhì)主要由中酸性巖組成，與較大范圍的花崗巖分布有關(guān)，石英、長(zhǎng)石含量高。

接收函數(shù)；地殼厚度；H—Kappa；疊加

0 引言

太行山山前斷裂帶位于太行山脈與華北平原的過(guò)渡地帶，北起北京懷柔附近，向南經(jīng)房山和河北省的淶水、保定、石家莊、邢臺(tái)、邯鄲及河南省的安陽(yáng)、湯陰至新鄉(xiāng)，總體呈NE—NNE向展布，全長(zhǎng)約620km，是華北及中國(guó)東部地區(qū)的一條重要構(gòu)造帶（徐杰等，2000）。本文主要研究太行山前斷裂帶南部地區(qū)，研究區(qū)域范圍為35°00′—36°54′N(xiāo)，113°00′—115°00′E，主要地質(zhì)構(gòu)造包括邯鄲—安陽(yáng)段和安陽(yáng)—新鄉(xiāng)段。

邯鄲—安陽(yáng)段主要是邯鄲斷裂南段。邯鄲斷裂帶北起邢臺(tái)以北，南至安陽(yáng)南斷裂，長(zhǎng)130km，由梧東、雙西和邯鄲斷裂組成。邯鄲斷裂是其主干斷裂，總體走向NE10°—20°，傾向SEE，傾角一般40°—60°，正斷裂性質(zhì)。此斷裂中、新生代有明顯活動(dòng)，據(jù)石油地質(zhì)勘探資料（石油物探局地質(zhì)研究院，1991），南段與臨漳斷裂一起控制邯鄲中、新生代凹陷的發(fā)育。邯鄲凹陷為雙側(cè)斷陷，堆積厚1 000 m的上侏羅統(tǒng)和白堊系，下第三系厚1 900 m。

安陽(yáng)—新鄉(xiāng)段為湯西和湯東斷裂及其控制的湯陰凹陷。湯東斷裂、湯西斷裂近平行分布于太行山山前斷裂帶南段，正斷裂性質(zhì)，共同控制湯陰凹陷的發(fā)育。湯西斷裂走向NE30°，傾向SEE，長(zhǎng)約70km。該斷裂南段早第三紀(jì)活動(dòng)強(qiáng)烈，湯陰凹陷的汲縣洼陷靠近斷裂下降盤(pán)一側(cè)，下第三系厚3 000 m。晚第三紀(jì)斷裂有過(guò)強(qiáng)烈活動(dòng)，控制了千余米厚的上第三系沉積。湯東斷裂走向NE30°，傾向NWW，長(zhǎng)90多千米。斷裂下降盤(pán)湯陰洼陷堆積厚2 800 m的下第三系，上第三系厚700余米。

該地區(qū)地震活動(dòng)相對(duì)頻繁，歷史記載的最大破壞性地震是1830年6月12日河北磁縣71/2級(jí)地震；1970年以來(lái)現(xiàn)代地震記錄顯示，磁縣、林州地震活動(dòng)相對(duì)密集，呈團(tuán)簇狀分布。綜合認(rèn)為，該區(qū)域發(fā)震構(gòu)造集中，地震危險(xiǎn)性較高，進(jìn)行地殼深部構(gòu)造特征研究對(duì)于認(rèn)識(shí)該地區(qū)發(fā)震機(jī)制具有重要意義。

1 數(shù)據(jù)處理

1.1 資料選取

研究區(qū)域位于河南、山西、河北交界地區(qū)，分布11個(gè)數(shù)字化測(cè)震臺(tái)站（圖1，表1）。地震臺(tái)站經(jīng)“九五”、“十五”數(shù)字化改造，主要利用8個(gè)寬頻帶地震儀（頻帶寬度 60 s，采樣率100次/s）記錄的地震事件進(jìn)行地殼深度反演。因此，安陽(yáng)臺(tái)（AN）、臨漳臺(tái)（LZT）和延津臺(tái)（YJ）的地震記錄不參與反演計(jì)算。

表1 地震臺(tái)站參數(shù)Table1 The parameter of the stations

選取2009年1月至2012年9月發(fā)生的695個(gè)5.5級(jí)以上遠(yuǎn)震，震中距30°—90°（在遠(yuǎn)場(chǎng)條件下，臺(tái)站下方體波入射角較小，可以看作近似垂直入射地震臺(tái)站下方），篩選65個(gè)均勻分布于研究區(qū)的地震事件（圖2），進(jìn)行接收函數(shù)反演計(jì)算。

圖1 區(qū)域地質(zhì)構(gòu)造及臺(tái)站分布Fig.1 The regional geological structure and the distribution of the stations

圖2 地震事件分布Fig.2 Epicenter distribution of the earthquakes

1.2 接收函數(shù)分析

接收函數(shù)方法最早由 Phinney提出，目的是用遠(yuǎn)震體波研究臺(tái)站下方地殼及上地幔結(jié)構(gòu)。接收函數(shù)不依賴臺(tái)站密度和空間分布，適用于臺(tái)站稀疏區(qū)域，可用遠(yuǎn)震數(shù)據(jù)反演臺(tái)站下方地殼結(jié)構(gòu)，有效解決某些地區(qū)地震較少而無(wú)法用近震震相進(jìn)行反演的問(wèn)題。遠(yuǎn)震體波徑向接收函數(shù)反演法是地震研究的一個(gè)重要方法，在研究地殼上地幔結(jié)構(gòu)中得到廣泛應(yīng)用（Owens J et al，1984；Ammon C J et al，1990； Langston C A，1997； Rai S S et al，2003），中國(guó)學(xué)者也取得一定成果（劉啟元等，1997，2000；吳慶舉等，1998；徐鳴潔等，2005）。

1.3 方法和原理

根據(jù)選取的遠(yuǎn)震事件，在臺(tái)站記錄中截取P波前50 s及后200 s波形數(shù)據(jù)。數(shù)據(jù)一般是SEED格式，由于計(jì)算程序使用SAC格式文件，需進(jìn)行數(shù)據(jù)格式轉(zhuǎn)換，并在頭文件中寫(xiě)入事件信息。

將直角坐標(biāo)系下原始波形數(shù)據(jù)旋轉(zhuǎn)變換為極坐標(biāo)系，根據(jù)臺(tái)站相對(duì)地震的方位角φ，對(duì)水平記錄的NS和EW向利用公式（1）進(jìn)行轉(zhuǎn)換，將地震觀測(cè)波形變換到徑向R和切向T（胡家福等，2003）。

式中，φ為相對(duì)于震中的方位角；CN和CE分別為臺(tái)站記錄的NS向和EW向地震觀測(cè)原始波形，CR和CT分別為徑向R和切向T上的相應(yīng)地震觀測(cè)波形。

采用時(shí)間域迭代反褶積方法進(jìn)行接收函數(shù)提取，求解徑向和切向分量的接收函數(shù)。以徑向分量為例進(jìn)行說(shuō)明：ER=t，求解方程最小值。DR（t）為臺(tái)站記錄遠(yuǎn)震波形數(shù)據(jù)的徑向分量，DT（t）為切向分量，DV（t）為垂直分量；DV（t）*ER（t）為不斷迭代更新的接收函數(shù)與觀測(cè)垂直分量地震圖的褶積（Li X Q and Nabelek J K，1999；Langston C A and Hammer J K，2001）。

進(jìn)行接收函數(shù)的H—Kappa迭代搜索，對(duì)于一維水平單層地殼模型，當(dāng)給定地殼內(nèi)介質(zhì)中P波和S波的平均速度vP和vS，由公式（2）即可求得地殼厚度H。式中tPs表示vS震相相對(duì)于初至P波的到時(shí)差，p為入射P波參數(shù)。

該方法需要解決莫霍面深度與波速比之間存在的折中問(wèn)題。

莫霍面多次反射轉(zhuǎn)換波攜帶更多約束信息，聯(lián)合利用莫霍面多次反射轉(zhuǎn)換波震相（PpSs，PpSs+PpSs）與初至P波到時(shí)差，可以有效獲得臺(tái)站下方的地殼厚度H和波速比。

因此，Zhu 等（2000）綜合利用莫霍面一次轉(zhuǎn)換震相Ps及多次轉(zhuǎn)換震相（PpSs，PsSs+PpSs）疊加搜索能量最大值方法，確定地殼厚度（H）和縱橫波速比（κ），具體公式為，其中r（t）為徑向接收函數(shù)的振幅，t1、t2、t3分別表示Ps、Pp Ss、PsSs+PpSs震相在預(yù)測(cè)地殼厚度H和 vP/vS條件下對(duì)應(yīng)的到時(shí)。ωt為Ps、PpSs、Ps Ss+PpSs震相的權(quán)重系數(shù)，Σωi=1。利用網(wǎng)格搜索方法可以確定函數(shù)r（H，κ）的最大值及相應(yīng)地殼厚度H和平均地殼介質(zhì)波速比κ，通過(guò)公式（3），可獲得相應(yīng)的平均地殼介質(zhì)泊松比σ。

對(duì)較好的臺(tái)站進(jìn)行疊加，每個(gè)臺(tái)站參與疊加分析的接收函數(shù)不少于30個(gè)，大部分臺(tái)站的接收函數(shù)顯示比較清晰的多次波，地殼厚度疊加過(guò)程中，地殼厚度H取25—45km，κ取1.5—3.0。

2 計(jì)算結(jié)果

對(duì)選取的65個(gè)遠(yuǎn)震事件1 452條地震記錄分析篩選，獲得301個(gè)接收函數(shù)，臺(tái)站方位角徑向接收函數(shù)見(jiàn)圖3，用H—Kappa方法反演，得到臺(tái)站下方地殼厚度和波速比（表2，圖4）。

圖3 地震臺(tái)站下方接收函數(shù)及H—Kappa反演結(jié)果（a）CHZ；（b）HSTFig.3 The receiver function under the station and the H-Kappa inversion results of the stations

3 分析

研究區(qū)域位于華北斷塊區(qū)二級(jí)構(gòu)造單元太行山斷塊。由于中生代燕山運(yùn)動(dòng)和新生代喜山運(yùn)動(dòng)共同作用，以太行山東麓斷裂為界，西部是整體隆起區(qū)，東部為相對(duì)下沉區(qū)，形成華北斷塊區(qū)西高東低的地貌特征（劉堯興等，2000）。

（1）區(qū)域內(nèi)地殼厚度分布總體趨勢(shì)是西部較厚，東部較薄，莫霍面總體呈小角度向西傾斜。西部長(zhǎng)治盆地地殼厚度為40.2km，東部平原地區(qū)地殼較薄為31.8km，地殼變化落差約9km。

（2）在太行山與平原交接地帶，特別是沿太行山弧形分布的主峰外緣，即一級(jí)夷平面和二級(jí)夷平面的分界線附近，地殼厚度變化較快（劉堯興等，2000）。接收函數(shù)反演地殼深度結(jié)果顯示，林州盆地底部地殼厚度相對(duì)較?。?3.1km），盆地邊緣地殼逐漸變厚，東部邊緣地殼深度變化較?。?4.3km），向西莫霍面快速變陡；沿林州—長(zhǎng)治一線，地殼變化梯度為111 m/km；焦作—衛(wèi)輝一線地殼變化梯度較大，地表25km距離地殼深度變化約2.6km，地殼變化梯度為104 m/km。

圖4 地震臺(tái)下方地殼厚度和波速比分布Fig.4 The crustal thickness and the wave velocity ratio distribution of the stations

表2 地震臺(tái)下方地殼厚度和波速比結(jié)果Table2 The crustal thickness and the wave velocityratio of the stations

焦作—衛(wèi)輝—林州一線是太行山隆起邊緣，南部有薄壁斷裂，北部是林州斷裂，薄壁斷裂（北段）上新世活動(dòng)劇烈，到早、中更新世時(shí)，斷裂活動(dòng)性仍較強(qiáng)烈，但自北而南活動(dòng)性逐漸減小，晚更新世以后，斷裂垂直活動(dòng)不明顯。林州斷裂是林州盆地的西部邊界斷裂，西盤(pán)（下盤(pán)）由太古界、下古生界地層組成高聳的太行山脈，標(biāo)高可達(dá) 1.2—1.6km，東盤(pán)（上盤(pán)）為下降盤(pán)。林州斷裂有一個(gè)近 2 000 m 寬的破碎帶，電阻率測(cè)深、激發(fā)極化測(cè)深和鉆孔資料表明，寬大的破碎帶除斷層角礫外，均被紅色粘土充填，該斷裂最新活動(dòng)時(shí)代為早、中更新世（劉堯興等，1987）。

（3）接收函數(shù)反演的地殼深度與人工地震測(cè)深結(jié)果吻合。中國(guó)地震局地球物理勘探中心曾在該區(qū)域布設(shè)一條人工地震測(cè)深剖面（劉堯興等，2000），東起山東省巨野縣，經(jīng)河南濮陽(yáng)市、湯陰、林州至山西省平順、屯留縣，全長(zhǎng)436km，方位角294°50′（圖5），由圖5可見(jiàn)，太行山隆起區(qū)地殼變厚，湯陰地塹地殼變薄，內(nèi)黃隆起地殼變厚，表明深部莫霍面或上地幔的隆起、凹陷與地面起伏呈鏡像對(duì)稱（劉堯興等，2000）。

圖5 菏澤—林州—長(zhǎng)治綜合剖面Fig.5 The comprehensive profile of Changzhi - Linzhou - Changzhi

將長(zhǎng)治（CHZ）、林州（LZH）、河順（HST）臺(tái)和?？h臺(tái)（XX）的莫霍面深度投影到人工地震勘測(cè)的地殼深度剖面上，發(fā)現(xiàn)其與人工地震勘測(cè)的莫霍面深度比較吻合（圖5），表明該區(qū)域地殼梯度帶主要遞進(jìn)方向?yàn)橛蓶|向西，同時(shí)證明該方法計(jì)算的莫霍面深度可以對(duì)人工地震測(cè)深結(jié)果進(jìn)行約束。

（4）大陸地殼的巖石成分可以劃分為酸性、中性、基性及富含流體的空隙/破裂或者部分熔融，對(duì)應(yīng)泊松比分別為σ＜ 0.26、0.26≤σ＜0.28、0.28≤σ＜0.30及σ≥0.30。區(qū)域內(nèi)泊松比值變化不大，數(shù)值范圍為0.23—0.27，表明該地區(qū)主要由中酸性巖組成。區(qū)域東北部永年臺(tái)（YON）及多數(shù)臺(tái)站泊松比值為0.23—0.25，與較大范圍的花崗巖分布有關(guān)，主要是石英、長(zhǎng)石含量高。焦作臺(tái)（JZ）、涉縣臺(tái)（SXT）、?？h臺(tái)（XX）泊松比值為0.26—0.27，表明長(zhǎng)英質(zhì)和鐵鎂質(zhì)成分含量相當(dāng)。作為渤海灣盆地西界的太行山山前斷裂帶基本發(fā)育于上地殼拆離滑脫構(gòu)造，類(lèi)似于美國(guó)盆地——山脈省的構(gòu)造特點(diǎn)：拆離斷裂往往不具穿殼性，而是置根于中、下地殼（Allmendinger R W et al，1987），但在斷裂下盤(pán)會(huì)產(chǎn)生一些高角度的開(kāi)斷裂，成為巖漿通道和巖墻侵入場(chǎng)所。

4 結(jié)束語(yǔ)

利用遠(yuǎn)震地震波接收函數(shù)對(duì)太行山前斷裂帶南段的地殼結(jié)構(gòu)進(jìn)行反演，所得地殼厚度和泊松比與該區(qū)域地質(zhì)環(huán)境一致性較好。由于受到臺(tái)站分布位置的局限，獲取的地殼厚度和泊松比數(shù)據(jù)還不能詳細(xì)反映該地區(qū)地殼結(jié)構(gòu)，有待深入探討。地殼厚度和泊松比參數(shù)對(duì)于分析地殼深處的地質(zhì)演化過(guò)程、淺部的地質(zhì)構(gòu)造類(lèi)型以及探究地殼物質(zhì)成分具有重要作用；其自身梯度變化特征與地震構(gòu)造帶內(nèi)地震活動(dòng)的密切關(guān)聯(lián)，為區(qū)域中強(qiáng)地震危險(xiǎn)性分析提供了判定依據(jù)。該研究方法已獲得很多學(xué)者認(rèn)可，豐富且高質(zhì)量的地殼厚度和泊松比研究成果將為地學(xué)研究、地震預(yù)報(bào)研究工作奠定堅(jiān)實(shí)的基礎(chǔ)。

本研究使用朱露培博士的反演程序，在此表示感謝。感謝朱元清教授、劉堯興研究員、王俊斐老師等給予的指導(dǎo)和幫助，感謝審稿專家提出的寶貴意見(jiàn)。

劉啟元，Painer K ，李順成，等.中國(guó)數(shù)字地震臺(tái)網(wǎng)的接收函數(shù)及其線性反演［J］.地球物理學(xué)報(bào)，1997，40（3）：356-368.

劉啟元，陳九輝，李順成，等.新疆伽師強(qiáng)震群區(qū)三維地殼上地幔S波速度結(jié)構(gòu)及其地震成因的探討［J］.地球物理學(xué)報(bào)，2000，43（3）：356-365.

劉堯興.豫北及其鄰區(qū)新構(gòu)造活動(dòng)特征和中、強(qiáng)地震發(fā)生的構(gòu)造條件［D］.北京：國(guó)家地震局地質(zhì)研究所，1987：90-91.

劉堯興，周慶，荊智國(guó)，等.豫北地區(qū)新構(gòu)造活動(dòng)特征及中長(zhǎng)期地震預(yù)報(bào)研究［M］.西安：陜西地圖出版社，2000.

徐 杰，高戰(zhàn)武，宋長(zhǎng)青，等.太行山山前斷裂帶的構(gòu)造特征［J］.地震地質(zhì)，2000，22（2）：111-113.

徐鳴潔，王良書(shū)，劉建華，等.利用接收函數(shù)研究哀牢山紅河斷裂帶地殼上地幔特征［J］.中國(guó)科學(xué)（D 輯），2005，35（8）：729-737.

胡家福，蘇有錦，朱雄關(guān)，等.云南的地殼S波速度與泊松比結(jié)構(gòu)及其意義［J］.中國(guó)科學(xué)（D） 輯，2003，33（8）：714-722.

吳慶舉，曾融生.用寬頻帶遠(yuǎn)震接 收函數(shù)研究青藏高原的地殼結(jié)構(gòu)［J］.地球物理學(xué)報(bào)，1998，4l（5）：669.

Allmendinger R W，Hange T A，Potter C J，Olive J.Tectonic heredity and layered lower crust in the Basin and Range Province，Western United States.In： Coward M P，Dewey J F，Hancock P L eds，Continental Extensional Tectonics ［J］.Blackwell Scientific Publition 1987：223-246.

Ammon C J，Randall G E，Zandt G.On the nonuniqueness of receiver function inversions［J］.Geophys Res，1990，5（B10）：15 303-15 3l8.

Langston C A.Structure under Mount Rainier，Washington inferred from teleseismic body waves［J］.Geophys Res，1997，84（B9）：4 749-4 762.

Langston C A and Hammer J K.The vertical corupouent P-wave receiver function［J］.Bull Seism Soc Am，2001，91：1 805-1 819.

Li X Q and Nabelek J K.Deconvolutiou of teleseismic body waves for enhancing structure beneath a seismometetarray ［J］.Bull Seism Soc Am，1999，89（1）：190-201.

Rai S S，Priestley K，Suryaprakasam K et al.Crustal shear velocity structure of the south Indian shield［J］.Geophys Res，2003，108（B2）：ESE10-l-ESEl0-12.

Thomas J，Owens Z and Taylor S R.Seismic evidence for an ancient rift beneath the Cumberland Plateau，Tennessee ：A detailed an alysis of broad and teleseismic P waveforms［J］.Geophys Res，1984，89（B9）：7 783-7 795.

Study on the crustal thickness in the south of the piedmont fault zone of Taihang mountains by teleseismic receiver function

Chen Rui1），2），Yan Jungang3），Yan Kai2），Wan Na2）and Zhang Yuchuan2）
1） China University of Geosciences（Wuhan） ，Hubei Province 430074，China
2） Earthquake Administration of Henan Province，Zhengzhou 450016，China
3） Handan Seismic Station，Hebei Province 056006，China

In this paper，we calculate receiver functions of body wave under the stations in the south of the piedmont zone of the piedmont fault zone of Taihang mountains from 3-component digital waveform data of teleseismic earthquake events and obtain the thickness and vP/vSradio in the crust of this area with H-Kappa stacking.The results show that:（1） Crustal thickness in the west of this area is thicker than that in the east.in general； The Moho tilts westward with small angle；In the western Changzhi basin，the crustal thickness is 40.2km，the thinner crust in eastern plain area is 31.8km，changes in the earth's crust gap is about 9km.（2） In Taihang mountains and the plain transition zone，especially along the outer edge of the main peak of Taihang Mountains arc distribution，namely near the boundary between first-level planation surface and second-level planation surface，the crust thickness changes quickly.（3） The crustal thickness calculated by teleseismic receiver function is in agreement with the artificial seismic exploration result，the crust of Taihang mountains uplift area becomes thick，that of Tangyin rift valley becomes thin，and that of Neihuang uplift area becomes thick，which shows that the depth of Moho or upper mantle’s height is symmetrical with the height of the ground.（4） Poisson ratio in the district is about 0.23-0.27，which indicates that the underground material in the area is mainly composed of intermediate acidic rocks，especially the widely distributed granite with much quartz，feldspar.

receiver function，crustal structure，H-Kappa，stacking

10.3969/j.issn.1003-3246.2015.05.001

陳睿（1979—），男，工程師，中國(guó)地質(zhì)大學(xué)在讀碩士，主要從事地震監(jiān)測(cè)預(yù)報(bào)研究工作。E-mail:chr8860@126.com

河南省地震局2014年科研基金資助

本文收到日期：2015-04-20