吳海波 姚運(yùn)生 申學(xué)林 趙凌云

1)中國(guó)地震局地震研究所，地震與大地測(cè)量重點(diǎn)實(shí)驗(yàn)室，武漢 430071

2)湖北省地震局，武漢 430071

0 引言

2014年3月27和3月30日在湖北省秭歸縣分別發(fā)生了MS4.5和MS4.9地震，中國(guó)地震臺(tái)網(wǎng)中心的定位結(jié)果分別為30.92°N，110.80°E和30.91°N，110.82°E;本研究分別稱為秭歸地震1和秭歸地震2。這2次地震均位于長(zhǎng)江三峽水庫(kù)南岸相距不足3km處，是三峽水庫(kù)蓄水后發(fā)生在秭歸地區(qū)的最大地震，也是繼2008年秭歸4.1級(jí)地震后該區(qū)又發(fā)生的2次中等強(qiáng)度地震。2次地震均有較強(qiáng)震感，但未造成人員傷亡和嚴(yán)重的破壞。地震發(fā)生在三峽庫(kù)區(qū)，因而引起了社會(huì)各界的極大關(guān)注。

三峽水庫(kù)蓄水以來(lái)，震區(qū)附近已發(fā)生了5次ML3.5以上地震，這2次地震的孕震背景如何?以及震后構(gòu)造應(yīng)力場(chǎng)如何變化及對(duì)區(qū)域未來(lái)地震的影響如何呢?對(duì)此，需要對(duì)該區(qū)域引發(fā)中等強(qiáng)度地震的活動(dòng)斷裂以及2次地震的震源特征深入討論。目前，對(duì)這2次地震發(fā)震構(gòu)造及震源討論的有關(guān)文獻(xiàn)并不多，湖北省地震局《秭歸地震現(xiàn)場(chǎng)考察報(bào)告》中2次地震烈度等震線呈長(zhǎng)軸SN向分布，并指出與仙女山斷裂活動(dòng)有關(guān)，但余震分布并未完全按仙女山斷裂的走向分布，那么這2次地震的震源是如何分布的，發(fā)震構(gòu)造又是什么呢?本文將從地震學(xué)方面并結(jié)合區(qū)域地質(zhì)構(gòu)造進(jìn)行綜合分析，通過(guò)震源體分布規(guī)律和震源機(jī)制探討2次地震的發(fā)震構(gòu)造及特征，對(duì)于三峽庫(kù)區(qū)水庫(kù)誘發(fā)地震研究和未來(lái)地震災(zāi)害趨勢(shì)評(píng)估具有重要意義。

1 地震構(gòu)造背景

長(zhǎng)江三峽庫(kù)區(qū)位于揚(yáng)子準(zhǔn)地臺(tái)中段北緣的渝東－鄂西臺(tái)褶帶，構(gòu)造上受到南部渝東－八面山弧形構(gòu)造帶與北部大巴山－大洪山弧形構(gòu)造帶相向逆沖推覆，在兩者相碰交接地帶形成奉節(jié)－巴東褶皺帶和秭歸向斜構(gòu)造盆地，因此，該區(qū)域構(gòu)造應(yīng)力場(chǎng)同時(shí)受控于這兩大構(gòu)造單元。

圖1 震中地質(zhì)構(gòu)造簡(jiǎn)圖(據(jù)1︰20萬(wàn)地質(zhì)圖改編)Fig.1 Geotectonic sketch of the epicenter area(adapted from the 1∶200000 geological map).

本次地震序列活動(dòng)位于三峽庫(kù)區(qū)中東段，在黃陵背斜與秭歸坳陷盆地交會(huì)處(圖1)，微構(gòu)造上在仙女山微地塊內(nèi)，由仙女山斷裂和九畹溪斷裂分別控制其東西邊界。這2條斷裂均形成于燕山運(yùn)動(dòng)期，經(jīng)燕山期和喜馬拉雅期強(qiáng)烈構(gòu)造變形，向下切割古生界和白堊系，沿?cái)嗔牙瓘埿纬蓴嘞莶鄣兀瑯?gòu)成了黃陵地塊西南側(cè)分界線的一部分。對(duì)斷層活動(dòng)年代測(cè)試結(jié)果表明，兩者最后1次強(qiáng)烈活動(dòng)時(shí)代均為早、中更新世，最新活動(dòng)年齡為15萬(wàn)年左右。仙女山斷裂走向NNW，傾向SW，自北向南由仙女山斷裂、都鎮(zhèn)灣斷裂和橋溝斷裂組成，本次地震序列活動(dòng)發(fā)生在北段，對(duì)于該段斷裂北端是否延伸穿越長(zhǎng)江一直存在爭(zhēng)議，一般觀點(diǎn)認(rèn)為該段尖滅于長(zhǎng)江南岸荒口一帶，并未穿過(guò)長(zhǎng)江(譚成軒，1991;王瑞江等，1995;畢珉峰等，2012)，但也有學(xué)者通過(guò)斷裂形變特征、力學(xué)性質(zhì)與交接關(guān)系認(rèn)為該斷裂北延過(guò)了長(zhǎng)江(楊淑賢等，1993)。九畹溪斷裂由2條平行近SN走向的斷裂組成，傾向E或W，西支中部在路子口橫穿過(guò)長(zhǎng)江，構(gòu)成仙女山微地塊的東側(cè)邊界。

2 震源機(jī)制解

2.1 理論方法

震源機(jī)制目前有多種求解與反演方法，傳統(tǒng)上采用初動(dòng)符號(hào)法或P，S波振幅比法(Snoke et al.，1984)。初動(dòng)符號(hào)法對(duì)臺(tái)站的分布位置和數(shù)量要求高，并且有些臺(tái)站的初動(dòng)振幅較弱不易清晰分辨，P，S波振幅比方法比初動(dòng)符號(hào)法精度高些，但它對(duì)輻射花樣敏感，易受地下介質(zhì)各向異性的影響，兩者的另一缺陷為求解結(jié)果的多解性問(wèn)題不能很好地解決，因此，基于體波和面波的波形擬合反演矩張量解方法發(fā)展起來(lái)并被廣泛應(yīng)用，Global CMT、USGS、ERI等研究機(jī)構(gòu)均已提出各自的矩張量反演方法并發(fā)布結(jié)果，但大多數(shù)采用遠(yuǎn)震體波或面波數(shù)據(jù)求解全球性大震(MW＞5.5)的震源機(jī)制(Dziewonski et al.，1983;Sipkin，1982，1986;Kawakatsu，1995);對(duì)于使用近臺(tái)甚至地方臺(tái)體波波形資料反演區(qū)域中小地震矩張量解的研究較少，近幾年在國(guó)內(nèi)較流行的CAP方法在對(duì)區(qū)域中小地震矩張量反演時(shí)，要求臺(tái)站的震中距＞100km，震級(jí)＞MS4.0。

對(duì)于區(qū)域中小地震矩張量的反演，希望盡量能反演更小的地震，并盡量使用更近臺(tái)站波形數(shù)據(jù)以提高信噪比，增強(qiáng)反演結(jié)果的準(zhǔn)確性。對(duì)此，德國(guó)地質(zhì)與自然資源研究所(BGR)提出了一種新的地震矩張量求解方法——Kiwi法(Cesca et al.，2010，2011)。該方法基于新提出的“程函源(Eikonal source)”震源模型，對(duì)不同震級(jí)地震使用不同頻率段濾波，適用于全球性大地震和區(qū)域中小地震，其理論波形采用全波擬合理論，分別使用頻率域理論振幅譜和時(shí)間域P或S體波波形反演矩張量，可采用震中距＞20km的近臺(tái)記錄波形，對(duì)震級(jí)＞ML4地震均有效，本文將采用該方法反演2次地震的矩張量解。

通常，理論地震圖可通過(guò)地震矩張量與格林函數(shù)微分的乘積來(lái)計(jì)算，在斷面∑上，矩張量密度mpq、震源產(chǎn)生的位移un可表示為

式(1)中Gnp，q為格林函數(shù)。在Kiwi方法中，格林函數(shù)已預(yù)先計(jì)算并建立格林函數(shù)庫(kù)，反演時(shí)只需引用該格林函數(shù)庫(kù)即可，本研究中使用Qseis軟件預(yù)先解算了三峽地區(qū)20～600km格林函數(shù)(Wang，1999)。矩張量密度mpq一般有6個(gè)獨(dú)立分量，對(duì)于“程函源”模型，把震源的點(diǎn)源模式和運(yùn)動(dòng)源模式統(tǒng)一為13個(gè)參數(shù)，分別為:中心點(diǎn)時(shí)空位置(x，y，z，t)、標(biāo)量地震矩M0、斷面產(chǎn)狀(走向φ，傾角δ，滑動(dòng)角λ)、上升時(shí)間τ、邊界半徑R、核中心相對(duì)位置(ns，nd)和相對(duì)破裂速率vr/vs，前8個(gè)為點(diǎn)源模式參數(shù)，后5個(gè)為運(yùn)動(dòng)源模式參數(shù)，反演過(guò)程分為3步，只需要前2步就可求得矩張量解，本研究采用的濾波頻帶范圍為0.025～0.5Hz、0.1～0.35Hz。

本方法采用非擬合誤差數(shù)評(píng)價(jià)理論擬合波形與觀測(cè)波形的擬合程度，亦反映了求解結(jié)果的可靠性。如果有N條波形，整體誤差數(shù)M可表示為

圖2 震中及臺(tái)站分布圖Fig.2 The distribution of the epicenters and stations.

2.2 觀測(cè)數(shù)據(jù)及速度模型

本文使用了湖北省測(cè)震臺(tái)網(wǎng)12個(gè)臺(tái)站和重慶測(cè)震臺(tái)網(wǎng)2個(gè)臺(tái)站的近臺(tái)觀測(cè)數(shù)據(jù)(圖2)，前述2次地震的震中距最小的臺(tái)站為興山臺(tái)，震中距38km;最大為萬(wàn)州臺(tái)，223km。所有臺(tái)站均為寬頻帶速度型數(shù)據(jù)，儀器頻帶范圍在50～0.017Hz，滿足矩張量反演的頻帶要求，由于Kiwi方法主要采用振幅譜擬合反演，該方法的一個(gè)特點(diǎn)是對(duì)數(shù)據(jù)的高頻和震相轉(zhuǎn)換細(xì)節(jié)不敏感，因此在反演過(guò)程中直接使用了包含更多高頻成分的速度型波形數(shù)據(jù)而沒(méi)有轉(zhuǎn)換成位移型數(shù)據(jù)。

合適的區(qū)域地殼速度結(jié)構(gòu)模型對(duì)于矩張量反演和地震精定位至關(guān)重要。前人對(duì)三峽及鄰區(qū)地殼速度構(gòu)造研究較多，早在20世紀(jì)80年代末通過(guò)對(duì)以奉節(jié)—觀音檔為主剖面的DSS研究，已初步給出了該區(qū)域地殼結(jié)構(gòu)速度模型(陳學(xué)波等，1994)。三峽水庫(kù)蓄水后利用分布在庫(kù)區(qū)及周邊的大量微小地震，采用近震層析成像技術(shù)，一些學(xué)者對(duì)三峽庫(kù)區(qū)地殼速度結(jié)構(gòu)進(jìn)行了精細(xì)分析(廖武林等，2007;李強(qiáng)等，2009;王小龍等，2013)，也有人采用人工源方式研究三峽及鄰區(qū)近EW向深部剖面結(jié)構(gòu)(Zhang et al.，2009)。不同的研究結(jié)果盡管在一些局部細(xì)節(jié)上存在差別，但總體層次速度結(jié)構(gòu)差異并不明顯，比如，李強(qiáng)等采用8次初始速度結(jié)構(gòu)對(duì)庫(kù)區(qū)周邊的反演結(jié)果顯示，局部速度不均勻性相對(duì)于初始模型＜15%(李強(qiáng)等，2009)，在某些明顯的地質(zhì)體上，速度結(jié)構(gòu)比較一致，比如黃陵背斜高速異常區(qū)和秭歸盆地低速區(qū)，由于反演中臺(tái)站的分布范圍較大，也參考了CRUST2.0模型和人工測(cè)深與重力剖面研究的結(jié)果(Zhang et al.，2009;張毅等，2012)，依據(jù)這些研究該區(qū)域莫霍面自西向東變淺，平均約36km，速度約7.95km/s，對(duì)此，我們給出了6層地殼內(nèi)速度結(jié)構(gòu)模型(表1，圖3)。

2.3 矩張量反演結(jié)果

采用Kiwi技術(shù)，經(jīng)過(guò)兩個(gè)步驟反演，分別得到2次地震的矩張量解(表2)，第1步中使用全波振幅譜技術(shù)在頻率域反演(圖4)，其中斷面產(chǎn)狀和震源深度均采用格點(diǎn)法搜索最佳解，圖4給出了秭歸地震1矩張量反演結(jié)果，以及最近4個(gè)臺(tái)站觀測(cè)譜和波形與理論譜和波形之間的擬合關(guān)系，可以看出，最近4個(gè)臺(tái)站的理論與觀測(cè)曲線擬合得比較好，全部臺(tái)站的整體擬合差為0.35，震源深度的最佳擬合點(diǎn)在7.7km，第2步使用P波數(shù)據(jù)在時(shí)間域波形反演，整體波形擬合差為0.56，并且大部分臺(tái)站觀測(cè)與理論波形曲線差異很小(圖4)，圖5為秭歸地震2的反演結(jié)果，觀測(cè)與理論曲線也擬合得比較好，說(shuō)明2次地震的矩張量反演結(jié)果是穩(wěn)定可靠的。

表1 地殼速度結(jié)構(gòu)模型Table1 The crustal velocity structure model

表2 2次地震的震源機(jī)制解Table2 Focal mechanisms of two earthquakes

圖3 地殼速度結(jié)構(gòu)模型Fig.3 The crustal velocity structure model.

3 震源與發(fā)震構(gòu)造特征討論

3.1 震源分布與深度特征

截至2014年4月27日，三峽水庫(kù)誘發(fā)地震監(jiān)測(cè)臺(tái)網(wǎng)共記錄到513次余震，采用雙差定位法對(duì)ML0.4以上地震進(jìn)行重新定位，計(jì)算中使用了三峽臺(tái)網(wǎng)15個(gè)子臺(tái)的觀測(cè)數(shù)據(jù)，震中距5～90km，地殼速度結(jié)構(gòu)模型采用前述的6層結(jié)構(gòu)，最后得到了501次高精度定位結(jié)果，其水平誤差總體在30～80m，垂直誤差為50～180m。通過(guò)波形擬合反演和雙差定位法得到2次主震震源深度:地震1為7.7km和7.99km，地震2為9.1km和9.01km，兩者差異很小，均在誤差范圍內(nèi)，表明了2種方法求解結(jié)果是可靠的。

圖4 秭歸地震1矩張量反演結(jié)果Fig.4 The inversion results of the moment tensor of earthquake 1.

圖5 秭歸地震2矩張量反演結(jié)果Fig.5 The inversion results of the moment tensor of earthquake 2.

圖6 秭歸2次地震震中分布圖Fig.6 The epicentral distribution of two Zigui earthquakes.

把本次地震序列精定位結(jié)果分開顯示，27日至29日為地震1的余震序列(圖6a)，30日至4月27日為地震2的余震序列(圖6b)，圖6 c為整個(gè)序列的分布。震中分布圖顯示，地震1的余震比較集中地沿NE向呈條帶狀展布，主震在地震簇中間部分;地震2序列沿NNW和NE向均有分布，但沿NNW向較多，主震偏在震群的東南，整個(gè)序列絕大多數(shù)地震分布在仙女山斷裂(F1)北端，未穿過(guò)長(zhǎng)江。依據(jù)湖北省地震局地震現(xiàn)場(chǎng)工作組的野外調(diào)查結(jié)果，2次地震的宏觀震中均位于秭歸縣屈原鎮(zhèn)西陵峽村，震中烈度均為Ⅴ度，其中地震2存在1個(gè)2km2左右的Ⅵ度異常區(qū)，兩者的等震線也均呈橢圓狀，長(zhǎng)軸均為NNW向，烈度沿長(zhǎng)軸方向衰減較慢，反映了2次地震可能受控于NNW向斷裂。與余震分布比較，地震1的余震分布長(zhǎng)軸方向與烈度等震線長(zhǎng)軸方向不一致;地震2的余震在NE和NNW向均有分布，不過(guò)在NNW向端部的余震與震群簇有點(diǎn)偏離(圖6b，c)，這些事實(shí)似乎表明，2次地震序列既受到NNW向斷裂控制，也與NE向構(gòu)造活動(dòng)有關(guān)。

根據(jù)2次地震的地表分布特征，按發(fā)震時(shí)間把序列分成6個(gè)時(shí)段:3月27日、28日、29日、30日、31日，4月1日至26日，并分別沿AA'、BB'、CC'和DD'做深度剖面(圖7)，其中AA'、BB'、CC'為同一方向。結(jié)果顯示，地震1序列的深度分布在4.5～9.0km(AA'剖面)，地震2序列的深度主要在5.5～10.0km(BB'剖面)，2次主震均在地震簇的底部，地震2序列總體上比地震1深，地震2序列更靠近仙女山斷裂。對(duì)于整體地震序列，在CC'剖面上，從C到C'地震深度逐漸變淺，呈現(xiàn)1條較明顯的斷面層aa'，傾向SW，傾角50°左右;DD'剖面上，從D到D'地震分布呈三角形狀，形成2條明顯的斷面線，頂點(diǎn)正好位于CC'和DD'線的交會(huì)處，從D至頂點(diǎn)深度逐漸變淺(虛線b);地震2在該斷面線的底部，從頂點(diǎn)至D'深度逐漸變深(虛線c)，地震1在該面底部，b和c的傾角均較小，在60°～80°之間，從CC'和DD'這種分布特征來(lái)看，如果2次地震沿DD'方向錯(cuò)動(dòng)，將與該剖面上余震分布形態(tài)相矛盾，因此從震中分布及震源深度剖面特征分析，我們認(rèn)為2次地震沿CC'向破裂滑動(dòng)更合理。

3.2 矩張量解特征分析

矩張量反演結(jié)果顯示2次地震均以走滑為主兼少量逆沖成分，且2個(gè)斷層節(jié)面的走向比較接近(表2)，均為NE和NNW向，與地震序列精定位分布結(jié)果相符合，但2次地震主壓應(yīng)力P軸方位相差180°，節(jié)面傾向相反，對(duì)于節(jié)面A，地震1傾向NW，傾角達(dá)79°，表現(xiàn)為左旋走滑特征，地震2傾向SE，傾角為68°，表現(xiàn)為右旋走滑特征，與震源深度DD'剖面比較，2次地震節(jié)面A與該剖面上2條斷面線b和c傾向一致，傾角也基本一致，從震中分布、震源深度剖面和矩張量解的分析認(rèn)為，NE向節(jié)面A為2次地震的發(fā)震斷面，而非NNW向的節(jié)面B。

圖7 震源深度剖面圖Fig.7 The profile of focal depth.

另外，反演結(jié)果給出了矩張量解3個(gè)分量的百分比:雙力偶部分(DC)、各向同性分量部分(ISO)和補(bǔ)償線性偶極子部分(CLVD)，通過(guò)比較，地震1的雙力偶分量多于地震2，2次地震各向同性成分均比較強(qiáng)，且地震2的各向同性分量更多，Thomsen(1986)和蔡曉剛等(2011)的研究結(jié)果指出，上地殼巖石和部分材料的ISO分量絕大部分低于10%，2次地震較強(qiáng)的各向同性分量可能與震區(qū)上地殼巖層長(zhǎng)期受庫(kù)水滲透與侵蝕有關(guān)，由于2次地震在同區(qū)域，介質(zhì)的差異相對(duì)較小，兩者DC和CLVD分量的差異可能與斷裂錯(cuò)動(dòng)中的張裂成分有關(guān)(蔡曉剛等，2011)。

3.3 區(qū)域地質(zhì)構(gòu)造特征分析

地質(zhì)上，本次地震序列震中分布在仙女山地塹北端，東西兩側(cè)為九畹溪斷裂與仙女山斷裂切割(圖8)，該地塹經(jīng)燕山和喜馬拉雅兩期運(yùn)動(dòng)，逆沖切割古生界和白堊系并下陷而形成，巖層表現(xiàn)為地塹為下志留統(tǒng)頁(yè)巖(S1)，東西兩側(cè)為寒武統(tǒng)白云質(zhì)灰?guī)r(∈3)，在CC'剖面上，aa'發(fā)震層處在頁(yè)巖層內(nèi)，該層與仙女山斷裂的傾角與傾向均一致，由于仙女山斷裂傾向SW，顯然本次地震序列并未正好發(fā)生在仙女山斷裂的斷層面上，而是在其偏東方向。

圖8 周坪東南九畹溪河谷兩側(cè)地質(zhì)地貌剖面示意圖(修改自袁登維等，1996)Fig.8 A geomorphic profile across Jiuwanxi Valley at the southeast of Zhouping(modified from YUAN Deng-wei et al.，1996).

在區(qū)域構(gòu)造應(yīng)力場(chǎng)方面，三峽水庫(kù)蓄水前，庫(kù)區(qū)構(gòu)造應(yīng)力場(chǎng)水平最大主壓應(yīng)力方向?yàn)镹EE或近EW向(高士鈞，1992;楊淑賢等，2005;王墩等，2007)，本次序列2次主震的震源機(jī)制解的最大主壓應(yīng)力方向與之基本一致，反映了這2次地震主要受控于區(qū)域構(gòu)造應(yīng)力場(chǎng)，但本次地震序列距三峽水庫(kù)很近，庫(kù)水對(duì)地震活動(dòng)的影響不可忽視，依據(jù)震中區(qū)地殼分層結(jié)構(gòu)，秭歸地震1深度上處于區(qū)域沉積蓋層底部，地震2已處在基底結(jié)晶花崗或變質(zhì)巖層內(nèi)，余震深度主要在5～10km，區(qū)域深部巖層主要為頁(yè)巖和灰?guī)r，這種巖性有利于庫(kù)水向深部巖層滲透與溶蝕，對(duì)局部地下介質(zhì)的物理性質(zhì)將產(chǎn)生影響，在前述矩張量解中，CD、ISO和CLVD分量百分比成分的差異可能是這一結(jié)果的反映，庫(kù)水經(jīng)過(guò)10多年的不斷滲透侵蝕，飽和了巖層中的孔隙，目前已穿透沉積蓋層到達(dá)基底結(jié)晶基巖上部，從而導(dǎo)致地下介質(zhì)ISO分量增加，CD成分減少。

本次地震活動(dòng)在DD'剖面上形成了明顯的NE向斷面，但從現(xiàn)有文獻(xiàn)中并未找到相應(yīng)較明顯的地表地質(zhì)構(gòu)造證據(jù)，在1︰20萬(wàn)地質(zhì)構(gòu)造圖中并未顯示該構(gòu)造線，不過(guò)，在CC'剖面上發(fā)震層的形態(tài)與仙女山斷裂面相一致，在DD'向剖面上2個(gè)斷面的傾角也與仙女山斷裂和九畹溪斷裂的傾角一致，在深度上，2個(gè)斷裂切穿上地殼到達(dá)中地殼，顯然本次地震序列活動(dòng)同時(shí)受到了2個(gè)斷裂的控制與影響，2個(gè)斷裂沿NW向構(gòu)造活動(dòng)，遇到仙女山北端障礙體影響，該部位由于受到庫(kù)水的長(zhǎng)期作用影響局部脆弱，因此在這里產(chǎn)生NE向階列錯(cuò)動(dòng)。本次地震活動(dòng)序列也進(jìn)一步暗示著仙女山斷裂北端可能并未向北穿過(guò)長(zhǎng)江。

4 結(jié)論

依據(jù)本次地震序列精定位后的震中分布、震源深度剖面以及2次主震矩張量解特征，結(jié)合野外現(xiàn)場(chǎng)宏觀調(diào)查結(jié)果和震中區(qū)地質(zhì)構(gòu)造等特點(diǎn)，綜合分析得到以下結(jié)論:

(1)經(jīng)綜合分析推斷仙女山斷裂北端NE向的破裂面為2次主震的發(fā)震斷面，而并非是NNW向。

(2)余震序列的NE和NNW向分布、剖面上發(fā)震層的形態(tài)和深度特征表明，本次地震活動(dòng)受到了仙女山斷裂和九畹溪斷裂活動(dòng)的控制與影響。

(3)2次主震矩張量解中CD成分偏少，ISO成分較多可能是庫(kù)水對(duì)震源區(qū)地下介質(zhì)物性影響的反映。

鑒于三峽水庫(kù)蓄水后震中區(qū)域已多次發(fā)生中小震群活動(dòng)，存在地震活動(dòng)增強(qiáng)的趨勢(shì)，本次地震序列主震發(fā)震斷面可能為新出現(xiàn)的斷面或隱伏構(gòu)造，今后應(yīng)更加關(guān)注該發(fā)震斷面區(qū)的地震活動(dòng)。

畢珉峰，楚全芝，鄧志輝，等.2012.長(zhǎng)江三峽地區(qū)仙女山斷裂北端延伸問(wèn)題探討［J］.地震地質(zhì)，34(2):294—301.doi:10.3969/j.issn.0253－4967.

BI Min-feng，CHU Quan-zhi，DENG Zhi-hui，et al.2012.Tectonic landform and location of the northern end of Xiannüshan Fault at the Three Gorges area［J］.Seismology and Geology，34(2):294—301(in Chinese).

蔡曉剛，姚陳，陳曉非.2011.各向異性ATI介質(zhì)剪切位錯(cuò)源地震矩張量［J］.地球物理學(xué)報(bào)，54(7):1772—1782.

CAI Xiao-gang，YAO Chen，CHEN Xiao-fei.2011.Seismic moment tensor in anisotropic ATI media:Shear faulting［J］.Chinese J Geophys，54(7):1772—1782(in Chinese).

陳學(xué)波，陳步云，張四維，等.1994.長(zhǎng)江三峽工程壩區(qū)及外圍深度構(gòu)造特征研究［M］.北京:地震出版社.1—29.

CHEN Xue-bo，CHEN Bu-yun，ZHANG Si-wei，et al.1994.Study of Deep Tectonic Characteristics in the Yangtze-Three Gorges Dam and Adjacent Area［M］.Seismological Press，Beijing.1—29(in Chinese).

高士鈞.1992.長(zhǎng)江三峽地區(qū)地殼應(yīng)力場(chǎng)與地震［M］.北京:地震出版社.15—32.

GAO Shi-jun.1992.Crustal Stress Field and Earthquake in the Three Gorges Area［M］.Seismological Press，Beijing.15—32(in Chinese).

李強(qiáng)，趙旭，蔡晉安，等.2009.三峽水庫(kù)壩址及鄰區(qū)中上地殼P波速度結(jié)構(gòu)［J］.中國(guó)科學(xué)(D輯)，39(4):427—436.

LI Qiang，ZHAO Xu，CAI Jin-an，et al.2009.P-wave velocity structure of upper and middle crust beneath the Three Gorges reservoir dam and adjacent region［J］.Science in China(Ser D)，39(4):427—436(in Chinese).

廖武林，姚運(yùn)生，丁志峰，等.2007.三峽地區(qū)P波速度層析成像研究［J］.大地測(cè)量與地球動(dòng)力學(xué)，27(3):80—84.

LIAO Wu-lin，YAO Yun-sheng，DING Zhi-feng，et al.2007.Tomographic imagery of P wave velocity structure in Three Gorges region［J］.Journal of Geodesy and Geodynamics，27(3):80—84(in Chinese).

譚成軒.1991.長(zhǎng)江三峽工程庫(kù)首區(qū)仙女山斷裂北延問(wèn)題研究［J］.水文地質(zhì)工程地質(zhì)，18(5):13—17.

TAN Cheng-xuan.1991.Northward extension of Xiannüshan fault in the front region of Yangtze Gorges project［J］.Hydrogeology and Engineering Geology，18(5):13—17(in Chinese).

王墩，姚運(yùn)生，薛軍蓉，等.2007.三峽水庫(kù)重點(diǎn)監(jiān)視區(qū)蓄水前后震源機(jī)制研究［J］.大地測(cè)量與地球動(dòng)力學(xué)，27(5):103—107.

WANG Dun，YAO Yun-sheng，XUE Jun-rong，et al.2007.Focal mechanism of key monitoring area of Three Gorges reservoir before and after impoundment［J］.Journal of Geodesy and Geodynamics，27(5):103—107(in Chinese).

王瑞江，譚成軒.1995.長(zhǎng)江三峽地區(qū)仙女山斷裂帶構(gòu)造活動(dòng)性及其北延問(wèn)題討論［J］.地球科學(xué):中國(guó)地質(zhì)大學(xué)學(xué)報(bào)，20(6):693—696.

WANG Rui-jiang，TAN Cheng-xuan.1995.A discussion of structural activity of Xiannüshan fault belt and its extension towards north in the Three Gorges on the Yangtze River［J］.Earth Science:Journal of China University of Geosciences，20(6):693—696(in Chinese).

王小龍，馬勝利，郭志，等.2013.利用地震背景噪聲成像技術(shù)反演三峽庫(kù)區(qū)及鄰近地區(qū)地殼剪切波速度結(jié)構(gòu)［J］.地球物理學(xué)報(bào)，56(12):4113—4124.

WANG Xiao-long，MA Sheng-li，GUO Zhi，et al.2013.S-wave velocity of the crust in Three Gorges reservoir and theadjacent region inverted from seismic ambient noise tomography［J］.Chinese J Geophys，56(12):4113—4124(in Chinese).

楊淑賢，高士鈞，蔡永建，等.2005.三峽及鄰區(qū)新構(gòu)造期以來(lái)應(yīng)力場(chǎng)分區(qū)研究［J］.大地測(cè)量與地球動(dòng)力學(xué)，25(4):42—45.

YANG Shu-xian，GAO Shi-jun，CAI Yong-jian，et al.2005.Study on stress field zoning in Three Gorges and its adjacent area since neotectonic age［J］.Journal of Geodesy and Geodynamics，25(4):42—45(in Chinese).

楊淑賢，周明禮，徐孝文，等.1993.秭歸盆地東緣斷裂構(gòu)造覓蹤:再論仙女山斷裂帶北延過(guò)長(zhǎng)江問(wèn)題［J］.地殼形變與地震，13(2):48—54.

YANG Shu-xian，ZHOU Ming-li，XU Xiao-wen，et al.1993.Seeking for east boundary fault structure at Zigui Basin:A further discussion on northward extension through Yangtze River of the Xiannüshan fault zone ［J］.Crustal Deformation and Earthquake，13(2):48—54(in Chinese).

袁登維，梅應(yīng)堂，秦興黎，等.1996.長(zhǎng)江三峽工程壩區(qū)及外圍地殼穩(wěn)定性研究［M］.武漢:中國(guó)地質(zhì)大學(xué)出版社.42—65.

YUAN Deng-wei，MEI Ying-tang，QIN Xing-li，et al.1996.Study of Crustal Stability in the Yangtze-Three Gorges Dam and Adjacent Area［M］.China University Geosciences Press，Wuhan.42—65(in Chinese).

張毅，陳超，梁青，等.2012.三峽地區(qū)中上地殼密度結(jié)構(gòu)［J］.地球科學(xué):中國(guó)地質(zhì)大學(xué)學(xué)報(bào)，37(增):213—221.

ZHANG Yi，CHEN Chao，LIANG Qing，et al.2012.Density structure of upper and middle crust in Three Gorges reservoir area［J］.Earth Science:Journal of China University of Geosciences，37(Suppl):213—221(in Chinese).

Cesca S，Heimann S，Dahm T.2011.Rapid directivity detection by azimuthal amplitude spectra inversion［J］.J Seismol，15(1):147—164(submitted).

Cesca S，Heimann S，Stammler K，et al.2010.Automated procedure for point and kinematic source inversion at regional distances［J］.Journal of Geophysical Research，115(B6):B06304.

Dziewonski A M，Woodhouse J H.1983.An experiment in the systematic study of global seismicity:Centroid moment tensor solutions for 201 moderate and large earthquakes of 1981［J］.J Geophys Res，88:3247—3271.

Kawakatsu H.1995.Automated near-real time CMT inversion［J］.Geophys Res Lett，22:2569—2572.

Sipkin S A.1982.Estimation of earthquake source parameters by the inversion of waveform data:Synthetic waveforms［J］.Phys Earth Planet Inter，30:242—259.

Sipkin S A.1986.Estimation of earthquake source parameters by the inversion of waveform data:Global seismicity，1981－1983［J］.Bull Seismol Soc Am，76:1515—1541.

Snoke J A，Munsey J W，Teague A C，et al.1984.A program for focal mechanism determination by combined use of polarity and SV－P amplitude ratio data［J］.Earthquake Notes，55(3):1—15.

Thomsen L.1986.Weak elastic anisotropy［J］.Geophysics，51(10):1954—1966.

Wang R.1999.A simple orthonormalization method for the stable and efficient computation of Green's functions［J］.Bull Seismol Soc Am，89:733—741.

Zhang Z J，Bai Z M，Waiter M，et al.2009.Crustal structure across the Three Gorges area of the Yangtze platform，central China，from seismic refraction/wide-angle reflection data［J］.Tectonophysics，475(3－4):423—437.