嚴川， 許力生*， 張旭， 李春來， 許康生

1 中國地震局地球物理研究所， 北京 100081 2 甘肅省地震局， 蘭州 730000

?

一種地方與區(qū)域地震震源機制反演技術：廣義極性振幅技術(二)
——對實際震例的應用

嚴川1， 許力生1*， 張旭1， 李春來1， 許康生2

1 中國地震局地球物理研究所， 北京 100081 2 甘肅省地震局， 蘭州 730000

為了檢驗廣義極性振幅技術(GPAT)的實用性，我們利用GPAT反演確定了49次實際地震的震源機制、矩震級和震源深度．為了檢驗GPAT對地方地震、區(qū)域地震以及遠震的實用性，我們選用了震級范圍約為ML0.2～MS7.0之間的地震與震中距范圍約在5～8000 km之間的觀測資料．對反演結果的分析表明，利用GPAT獲取的震源機制結果在合理誤差范圍內是正確的，利用GPAT獲取的矩震級結果是可靠的，而利用GPAT獲取的震源深度總體上似乎比常規(guī)定位深度深約0.6 km．總體而言，GPAT在震源機制、矩震級和震源深度的反演方面表現(xiàn)出良好的實用性．

地方與區(qū)域地震； 震源機制反演； 廣義極性振幅技術； 實際震例

The GPAT employs polarities of the P-wave first-motion and the maximal amplitude with polarities. The phases used here are generalized, which are no longer limited to well-known P-waves, S-waves or surface waves. The maximal amplitudes can be adopted as long as they are clearly recorded in spite of the well-known phases. GPAT uses the correlation coefficient between observed data and synthetic data as the objective function, and sets a weight factor for the first-motion polarities and the maximal amplitudes. In this paper we select forty-nine earthquakes with epicenter distances from 5 km to 8000 km, with magnitudes fromML0.2 toMS7.0, including six far-field events in and surrounding China, twelve regional events in Lushan, Sichuan, eight local events in Yingjiang and twenty-three local events in Qiaojia, Yunnan.

The focal mechanisms of the far-field events determined with GPAT are highly consistent with those from GCMT(Global Centroid Moment Tensor). The focal mechanisms of the regional events determined with GPAT are consistent with those from other studies, and they appear to be similar to each other. For the local events in Yingjiang, the focal mechanism of the main shock is consistent with those from other studies and other events have the mechanisms fairly similar to the main shock. For the local events in Qiaojia, they all have similar strike-slip mechanisms. Moreover, the moment magnitudes from the GPAT are reliable，but the focal depths from the GPAT seem to be about 0.6 km deeper than the routinely determined depths on average.

The GPAT is proved to have a good performance in the inversion for focal mechanism，moment magnitude and focal depth．It is not only suitable for far-field data, but also for regional and local data. It is not only suitable for moderate-size events, but also for small and micro events.Keywords Local and regional earthquakes； Focal mechanism inversion； Generalized polarity and amplitude technique； Real earthquakes

1 引言

對于遠場地震而言，如果地震產(chǎn)生的位移已知，我們可以通過計算描述地震波傳播路徑效應的格林函數(shù)，借助于波形反演技術很容易得到其震源機制(Aki and Richards, 1980; Lay and Wallace, 1995)．這種波形反演技術最初適用于遠場記錄或較大的地震(Dziewonski et al.,1981)，后經(jīng)國內外同行的努力將其應用于區(qū)域地震記錄或中等大小的地震(Helmberger and Engen,1980; Dreger and Helmberger, 1991; Dreger and Helmberger, 1993; Zhao and Helmberger, 1994; Zhu and Helmberger, 1996)．但這種技術對于小震或微震的震源機制的求解仍無能為力．對于小震或微震的震源機制求解，P波初動符號和/或振幅信息是必不可少的資料， 但究竟如何利用這些資料確定震源機制卻一直沒有被廣泛采納的方法(Kisslinger, 1980; Snoke et al., 1984; Reasenberg and Oppenheimer, 1985; Ebel and Bonjer, 1990; Julian and Foulger, 1996; Hardebeck and Shearer, 2002, 2003;許忠淮等，1983; 梁尚鴻等，1984;俞春泉等，2009)．作為對已有研究工作的補充與發(fā)展，我們提出了一種反演技術——廣義極性振幅技術(GPAT)，并已經(jīng)通過一系列數(shù)值實驗檢驗了這種技術的可行性和抗干擾能力(嚴川和許力生，2014)．

2 震例與數(shù)據(jù)篩選

為了檢驗GPAT的實用性，我們選擇的地震必須包括具有遠場記錄的大震、具有區(qū)域記錄的中強震以及具有地方記錄的小震或微震．同時需要考慮，所選地震的反演結果盡可能具有可比性．為此，挑選了近年來發(fā)生在我國大陸及周邊的6次地震作為遠震震例(表1, 圖1a)，發(fā)生于2013年4月的蘆山地震序列(蘆山震群)12次地震作為區(qū)域地震震例(表2, 圖1b)，發(fā)生于2014年5月的云南盈江地震序列(盈江震群)的8次地震(表3, 圖1c)以及巧家流動臺陣記錄到的一個地震序列(巧家震群)的23次地震作為地方震的震例(表4, 圖1d)．其中表3和表4中定位結果的不確定性意義已在前文(許力生等，2013a)作過介紹．遠震的數(shù)據(jù)從IRIS(Incorporated Research Institutions for Seismology，美國地震學研究機構聯(lián)合會)數(shù)據(jù)中心獲取，反演結果可以與GCMT結果對比．蘆山震群的數(shù)據(jù)從我國國家臺網(wǎng)備份中心獲取(鄭秀芬等，2009；Zheng et al., 2010)，部分反演結果可以和利用其他反演方法得到的結果對比(呂堅等，2013；趙博等，2013)，也可以在群內地震之間進行對比．盈江震群的數(shù)據(jù)從云南地震臺網(wǎng)數(shù)據(jù)中心獲取，部分反演結果可以和其他機構利用不同反演技術得到的結果對比，同時也可以在群內地震之間進行對比．巧家震群的數(shù)據(jù)從我們布設的巧家流動地震臺陣獲取，反演結果只能在震群內部地震之間進行對比．注意，在這里我們不得不借助于一個并不嚴謹?shù)募僭O，即一個震群內部的主要事件應該具有類似的震源機制．

表1 遠震事件的震源參數(shù)Table 1 Source parameters of the teleseismic events

注: 本表參數(shù)來自IRIS．

表2 蘆山震群的震源參數(shù)Table 2 Source parameters of the Lushan earthquake swarm

注：本表參數(shù)來自國家臺網(wǎng)目錄，5級以上者為MS，5級以下者為ML．

表3 盈江震群的震源參數(shù)Table 3 Source parameters of the Yingjiang earthquake swarm

注：本表參數(shù)來自云南臺網(wǎng)目錄，5級以上者為MS， 5級以下者為ML．

表4 巧家震群的震源參數(shù)Table 4 Source parameters of the Qiaojia earthquake swarm

圖1 震中與臺站分布(a) 6次大震； (b) 蘆山震群； (c) 盈江震群； (d) 巧家震群. 紅色圈為地震震中，藍色三角形為臺站，黑色線為斷層.Fig.1 Distribution of the epicenters and the stations(a) The 6 large events； (b) The Lushan swarm； (c) The Yingjiang swarm； (d) The Qiaojia swarm．The red circles are the epicenters，the blue triangles are the stations，and the black lines are the faults．

圖2 地震記錄的震中距覆蓋范圍(a) 6次大震； (b) 蘆山震群； (c) 盈江震群； (d) 巧家震群.Fig.2 The epicentral distances of the used seismic recordings(a) The 6 large events； (b) The Lushan swarm； (c) The Yingjiang swarm； (d) The Qiaojia swarm．

圖3 GPAT震源機制結果與其他研究結果的比較(a) 6次大震的GPAT結果(紅色)與GCMT結果(黑色)，其中上方的數(shù)字為事件的編號; (b) 蘆山震群的GPAT結果(紅色)與趙博等(2013)的結果(紫色)以及呂堅等(2013)的結果(黑色); (c) 盈江震群的GPAT結果(紅色)與其他機構對19號事件的結果(黑色)以及其他機構對20號事件的結果(紫色); (d) 巧家震群的GPAT結果(紅色)．Fig.3 Comparison of the GPAT mechanisms with the other-source ones(a) The GPAT results (red) and the GCMT ones (black) of the 6 large events，where the numbers are the event codes; (b) The GPAT results (red) and the Zhao et al. (2013)′s ones (purple) as well as Lü et al. (2013)′s ones (black) of the Lushan swarm; (c) The GPAT results (red) and the results of the 19th event from other organizations (black) as well as the results of the 20th event from other organizations (purple) of the Yingjiang swarm; (d) The GPAT results (red) of the Qiaojia swarm．

3 數(shù)據(jù)預處理與格林函數(shù)計算

表1和圖1a所示遠震所使用數(shù)據(jù)的震中距分布在約500 km到8000 km之間(圖2a)，覆蓋了部分區(qū)域距離和遠場距離．在采集反演信息前，我們首先將原始記錄進行采樣，使其采樣率降為2 sps(samples per second)，然后對記錄進行濾波以去除噪聲的影響且保證感興趣的震相特征不發(fā)生明顯變化，所以對較小的事件(事件1、5、6)，采用了0.05～0.2 Hz的濾波窗，而對較大的事件(事件2、3、4)，采用了0.02～0.1 Hz的濾波窗．在進行震源機制反演時，我們采集了垂向記錄的P波初動信息、直達P波最大振幅信息和直達S波最大振幅信息．對于這些遠場地震，均使用IASPEI91全球速度模型，利用反射-折射率方法計算格林函數(shù)(Kennett, 1983; Kennett and Engdahl, 1991; Kennett et al., 1995)．

表2和圖1b所示的區(qū)域地震所使用數(shù)據(jù)的震中距分布在50 km至550 km之間(圖2b)，覆蓋了部分地方距離和區(qū)域距離．在采集反演信息前，首先將原始記錄進行采樣，使其采樣率降為10 sps，然后對記錄同樣進行濾波以去除噪聲的影響且保證感興趣的震相特征不發(fā)生明顯變化，在這里采用了0.5～2 Hz的濾波窗．在進行震源機制反演時，使用了垂向記錄的P波初動信息、直達P波最大振幅信息和直達S波最大振幅信息．格林函數(shù)同樣是根據(jù)IASPEI91全球速度模型，利用反射-折射率方法計算得到(Wang,1999)．

表3和圖1c所示的盈江震群所使用數(shù)據(jù)的震中距分布在約55 km至105 km之間(圖2c)．考慮到地方震群的地震較小，波形記錄的頻率較高，所以資料的采樣率和濾波頻率須根據(jù)具體情況而定．設置濾波頻率的原則是，消除背景噪聲并保證波形主頻不變．根據(jù)這組震群資料的實際情況，我們將這些事件的濾波范圍設置為2～4 Hz；設置采樣率的原則是，采樣頻率為濾波頻率上限的5倍．類似地，計算格林函數(shù)時，仍然使用反射-折射率方法(Wang,1999)，但是速度模型則是王椿鏞等(2002)為川滇建立的水平分層速度模型(表5)．這個模型是我們從已有的多個研究結果中挑選而來的(許力生等，2013b)．

表4和圖1d所示的巧家震群所使用數(shù)據(jù)的震中距分布在約5 km至100 km之間(圖2d)．巧家震群的震級更小，波形記錄的頻率更高．所以我們仍然采用處理盈江震群數(shù)據(jù)的技術流程，不同的地震采用不同的采樣率和頻率窗．同樣，采用王椿鏞等(2002)為川滇建立的水平分層速度模型和反射-折射率方法(Wang,1999)計算格林函數(shù)．

表5 根據(jù)王椿鏞等(2002)建立的速度模型Table 5 Velocity model modified after Wang et al. (2002)

4 GPAT震源機制

4.1 遠震

表6和圖3a展示了利用廣義極性振幅技術反演的6次遠震的震源機制結果(GPAT)和全球矩心矩張量結果(GCMT)．僅從圖3a中的沙灘球幾乎無法辨識GPAT結果與GCMT結果的差異．如果仔細比較對應參數(shù)，兩種結果之間還是存在或多或少的差別．我們注意到，只有一次事件的節(jié)面I的走向相差29°，其余都在15°以內；只有一次事件的節(jié)面II的走向相差29°，其余都在10°以內；所有節(jié)面I的傾角差異都在16°以內；所有節(jié)面II的傾角差異都在10°以內；只有兩次事件的節(jié)面I的滑動角相差17°，其余差別都在7°以內；只有兩次事件的節(jié)面II的滑動角相差24°和22°，其余都在16°以內．雖然GCMT結果是目前國際上的權威結果，但仍不能被認為是“真理”．假如我們把GPAT結果和GCMT結果的中間值作為真理的更好近似，那么GPAT結果的偏差就更?。纱丝磥?，GPAT利用遠場資料和區(qū)域資料反演遠場大震的震源機制是可行的．

4.2 區(qū)域震

表7和圖3b展示了蘆山震群的震源機制反演結果和部分來自其他研究的結果(具體參數(shù)見附表1和2)．通過對比發(fā)現(xiàn)GPAT反演結果與來自其他研究的結果之間的相似性(呂堅等，2013；趙博等，2013)．同時，我們注意到震群內各事件的震源機制也彼此接近．考慮到來自其他研究的結果遠不及GCMT結果那樣權威，我們并沒有仔細計算具體參數(shù)之間的差別．

品牌影響力的路徑系數(shù)為0.73，也說明其對新疆農產(chǎn)品區(qū)域品牌競爭力影響作用顯著，且具有很大提升空間，因此應加強構筑新疆民族特色的品牌形象，注重結合新疆特有的歷史、人文資源，在產(chǎn)品包裝形象、生產(chǎn)環(huán)境形象、經(jīng)營業(yè)績形象、社會貢獻形象、員工層次形象等方面突出新疆特色，增強品牌的文化品位；結合新疆農產(chǎn)品已經(jīng)在消費者心中留下的美好印記，進而充分利用和開發(fā)本地歷史文化資源，緊跟時代特色，為農產(chǎn)品區(qū)域品牌做歷史文化背書。將新疆的區(qū)域品牌培育和產(chǎn)品營銷結合，進行品牌營銷，注重品牌產(chǎn)品在營銷中的形象，進而提升品牌忠誠度。

需要特別說明的是，我們不但利用遠場記錄反演了蘆山主震的震源機制，還分別利用區(qū)域記錄和地方記錄反演了這次地震的震源機制．表6中事件4的GPAT結果為遠場資料的反演結果，表7中的結果7和7*分別為使用區(qū)域資料和地方資料的反演結果．對于節(jié)面I而言，三種結果的走向分別為21°、22°和358°，最大差別24°；三種結果的傾角分別為56°、34°和58°，最大差別為24°；三種結果的滑動角分別為64°、86°和58°，最大差別為28°．對于節(jié)面II而言，三種結果的走向分別為242°、207°和228°，最大差別35°；三種結果的傾角分別為42°，56°和44°，最大差別為14°；三種結果的滑動角分別為123°、93°和130°，最大差別為37°．我們對同一地震使用了不同的資料，盡管反演結果之間存在差別，但三種不同資料的反演結果仍具有一致性．

表6 遠震事件的GPAT與GCMT震源機制Table 6 GPAT and GCMT mechanisms of the teleseismic events

表7 蘆山震群的GPAT震源機制Table 7 GPAT mechanisms of the Lushan earthquake swarm

綜合以上分析，我們認為GAPT利用區(qū)域地震記錄能夠反演得到中強地震的震源機制．

4.3 地方震

4.3.1 盈江震群

表8和圖3c展示了盈江震群的震源機制參數(shù)，其中包括發(fā)生于2014年5月23日和5月30日的兩次主震事件．從反演結果看，盡管這些地震的震源機制不盡相同，但它們皆以陡傾角走滑為主．

對于盈江震群中的余震事件，我們無法為所有事件找到像GCMT那樣的權威結果作比較，但從已有的研究中找到其中部分事件的震源機制(圖1c和附表3、4)．盡管他們使用的方法和資料與我們使用的方法和資料截然不同，但這些事件的震源機制彼此接近．另外，我們已經(jīng)利用事件21至26成功提取了兩次主震的震源時間函數(shù)(許力生等，2014)，該結果表明余震和主震的震源機制的確很相似．

4.3.2 巧家震群

表9和圖3d展示了巧家震群的震源機制結果．清楚地可以看到，盡管這些地震的震源機制彼此存在差異，但它們都以走滑為主．我們無法為這些地震找到可以比照的結果，但這些事件震源機制間的相似性為求解結果的準確提供了保障.對于事件的時空分布而言，這些事件在時間上比較集中，大部分發(fā)生于2013年11月16日至23日之間，但是在空間上它們并不呈線性排列.對于事件的震級大小而言，這個震群的地震都比較小，最大的4.6級，最小的僅為0.2級，這樣小的地震完全可以在局部應力場作用下發(fā)生，因此過高地渴求震源機制的一致性顯然是不合理的．

綜合考慮盈江震群與巧家震群的反演，我們完全有理由認可GPAT利用地方地震記錄反演小微地震震源機制的實用性．

需要說明的是，表7、8和9中反演參數(shù)誤差(不確定性)的意義已在前文作過介紹(嚴川和許力生，2014)．

表8 盈江震群的GPAT震源機制Table 8 GPAT mechanisms of the Yingjiang earthquake swarm

表9 巧家震群的GPAT震源機制Table 9 GPAT mechanisms of the Qiaojia earthquake swarm

5 GPAT矩震級

矩震級的概念誕生于對大地震的研究，是由標量地震矩的大小決定的(Kanamori，1977；Hanks and Kanamori, 1979)．地方震級的概念誕生于對小地震的研究，是由某個特定周期的振幅的大小決定的(陳運泰和劉瑞豐，2004)．由于二者的定義不同， 二者之間存在差別便不難理解．不同震級之間的關系一直以來都是一個重要的研究方向，矩震級與地方震級的關系是其中重要的組成部分(Bakun and Lindh, 1977; Bakun, 1984; Fletcher et al., 1984; Hanks and Boore, 1984; Kim et al., 1989; Vidal and Mungúia, 1991; Ben-Zion and Rice, 1993; Ben-Zion and Zhu, 2002; Clinton et al., 2006; Tan and Helmberger, 2007)．目前的研究表明，兩者之間呈現(xiàn)分段關系，對于較大地震，二者比較一致，地震越小，二者差別越大．

單獨比較某個或某幾個地震的矩震級和地方震級是無意義的，在這里我們把本文中涉及的所有地震事件的矩震級和常規(guī)定位時確定的地方震級(注意：部分較大事件實際上是面波震級)進行對比分析．在圖4中，把同一事件的兩種不同震級的平均值作為“真值”(假定為正確的震級)，并依此真值從小到大的順序對震級進行排列．可以看出，矩震級與地方震級總體上具有相同的趨勢特征，但是，大于約3.5時，矩震級與地方震級一致性更好(注：沒有找到那6次遠震的ML震級)，二者幾乎重合，而小于約3.5時，矩震級與地方震級明顯分裂開來，且震級越小分裂越明顯，在地方震級為0時，二者的差異可達到2級．這種分段特征無疑是對已有研究結果的強力支持．

圖4 矩震級MW與地方震級ML 或者MS的比較縱軸和橫軸均為事件震級，其中藍色圈為來自常規(guī)定位的地方震級；紅色圈為來自震源機制反演的矩震級；空心圈為二者的平均值．Fig.4 Comparison of the moment magnitudes MWwith the local magnitudes ML or MSHorizontal and vertical axes both are for magnitudes of the events．The blue circles are from the routine locating and the red ones are from the mechanism inversion, and the empty ones，the average of the above two．

圖5 震源深度的比較(a)中M表示兩種深度的均值. 在(a)圖中藍色圈為常規(guī)定位時得到的震源深度，紅色圈為震源機制反演時得到的震源深度，空心圈為二者的平均值．在(b)圖中縱坐標為GPAT深度與平均值之差，紅色虛線為差的平均值．橫坐標均為事件代碼．Fig.5 Comparison of the focal depthsIn subplot (a) the blue circles are from the routine location，the red circles，from the mechanism inversion，and the empty ones，the average of the above two；in subplot (b) the vertical axis shows the differences between the GAPT depths and the average ones，and the red dashed line shows the average value of the differences. The horizontal axes are the event codes．

考慮到其他研究結果，盡管GPAT矩震級與常規(guī)定位時得到的地方震級不總是一致，但我們仍然有理由相信GPAT矩震級是可以信賴的．

6 GPAT震源深度

通過震源機制反演確定的震源深度(G深度)是否正確，用某個事件或者某幾個事件同樣難以得出結論．類似地，我們把所有事件的G深度與通過常規(guī)定位得到的震源深度(L深度)進行整體的比較．如圖5a，我們把同一事件的兩種震源深度的平均值作為“真值”(假定為正確的震源深度)，并依此真值從小到大的順序排列．可以看出，G深度和L深度的相關性非常好，而且對于相當數(shù)量的事件，二者幾乎是重合的．這說明G深度是一個具有實際意義的可用參數(shù)．為了做出定量評價，我們計算了G深度相對于真值的偏差(圖5b)，最小偏差為0.15 km，最大偏差為2.95 km，平均偏差為0.61 km，其中偏差小于2 km者占約67%．而0.61 km的系統(tǒng)差別很可能來自起始破裂點與“矩心”之間的差別．

由此可見，GPAT深度與常規(guī)定位的深度具有相當好的一致性，GPAT深度能夠反映震源的深度．

7 討論與結論

關于GPAT的原理和技術流程，我們已經(jīng)做過系統(tǒng)的介紹(嚴川和許力生，2014)；我們還利用各種數(shù)值實驗闡述了GPAT的可行性和抗干擾能力(嚴川和許力生，2014)．GPAT是為利用區(qū)域和地方地震記錄反演小微地震的震源機制而提出的，但其原理允許利用所有地震記錄反演所有地震的震源機制．在這項研究中，我們不但選擇了利用區(qū)域和地方地震記錄反演中、小、微地震的震源機制的震例，也選擇了利用遠場和區(qū)域地震記錄反演大地震的震源機制的震例．這不單是為了說明，GPAT不僅適合于小微地震，也適合于中強地震，不僅適用于地方和區(qū)域地震記錄，也適用于遠場記錄；更重要的是，對于中強地震，可以找到用其他方法和其他資料得到的震源機制便于和GPAT結果進行比較，而對于小微地震，難以找到強有力的佐證．

利用區(qū)域和遠震記錄反演大震的震源機制時，我們選擇了6個MS5.7～MS7.0之間的地震．所用記錄的震中距從500 km至8000 km，覆蓋了部分區(qū)域到遠場的距離．所用的地震的震源機制有正斷層、逆斷層、走滑斷層以及其他類型的斷層，覆蓋幾乎所有可能的機制．通過GPAT結果和GCMT結果的對比發(fā)現(xiàn)，二者基本一致，這說明GPAT能夠利用區(qū)域記錄和遠場記錄確定大震的震源機制．

利用地方和區(qū)域地震記錄反演中強地震的震源機制時，我們選擇了ML4.5～MS6.9之間的地震．所用記錄的震中距從50 km至550 km，覆蓋了部分區(qū)域距離和部分地方距離．所用的地震為2013年四川蘆山地震序列的主震和較大余震(蘆山震群)，這些地震的震源機制以逆沖為主．經(jīng)驗表明，一個地震序列的較大余震與主震具有類似的震源機制(盡管并非總如此)，而這個地震序列的主震和部分較大余震的震源機制已經(jīng)有一些結果(呂堅等，2013；趙博等，2013)．GPAT反演結果表明，蘆山震群的震源機制都非常接近．無論是震群內部事件之間震源機制的相似性，還是部分事件的震源機制與其他研究結果的相似性均可表明，GPAT反演結果在合理的誤差范圍內是正確的．這說明GPAT能夠利用區(qū)域和地方記錄確定中強地震的震源機制．

利用地方地震記錄反演小微地震的震源機制時，我們選擇了ML0.2～MS6.1之間的地震．所用記錄的震中距從5 km至105 km，覆蓋了幾乎整個地方距離．所用的地震一部分來自2014年云南盈江地震序列的主震和較大余震(盈江震群)，另一部分來自2013年發(fā)生在云南巧家的一個地震序列(巧家震群)．GPAT反演結果表明，盈江震群的事件都具有走滑機制．在震群內部事件之間的震源機制都彼此接近，主震的GPAT結果和其他機構給出的結果也較為一致(圖7)．同樣，巧家震群的事件也具有走滑機制．對于盈江震群和巧家震群，我們無法找到像GCMT那樣令人信服的結果可以比對，但震群內部事件之間的相似性在相當程度上能夠使我們相信，GPAT能夠利用地方地震記錄確定小微地震的震源機制．

GPAT不但反演確定震源機制，同時還反演確定矩震級和震源深度．從某一事件的反演結果，我們無從判斷其正誤．但所有事件的矩震級與地方震級的對比以及GPAT震源深度與常規(guī)定位深度的對比分析表明，GPAT反演得到矩震級和震源深度是可以信賴的．

由此可見，GPAT對實際地震是可用的，它能夠在合理的誤差范圍內反演確定地震事件的震源機制、震源深度和矩震級；它不但適合于小微地震，還適合于中強地震；它不但適用于地方和區(qū)域地震記錄，也適用于遠場地震記錄．

致謝 感謝IRIS數(shù)據(jù)中心提供遠震波形數(shù)據(jù)，感謝云南省地震局提供盈江地震序列波形數(shù)據(jù)，也感謝國家地震臺網(wǎng)中心提供蘆山地震序列波形數(shù)據(jù)．

附錄 來自其他研究的部分震源機制結果

附表1 趙博等(2013)確定的蘆山震群部分地震震源機制Appendix Table 1 Focal mechanisms of some of the Lushan swarm from Zhao et al. (2013)

附表2 呂堅等(2013)確定的蘆山震群部分地震震源機制Appendix Table 2 Focal mechanisms of some of the Lushan swarm from Lü et al. (2013)

Aki K, Richards P G. 1980. Quantitative Seismology Theory and Methods, Vol. I. W. H. San Francisco, California: Freeman and Company.

附表3 19號地震的其他震源機制結果Appendix Table 3 Focal mechanisms of the No.19 event

注：Mww:W-phase moment tensor;Mwb: Body-wave moment tensor;Mwc: Centroid moment tensor.

附表4 20號地震的其他震源機制Appendix

Bakun W H, Lindh A G. 1977. Local magnitudes, seismic moments, and coda durations for earthquakes near Oroville, California.Bull.Seism.Soc.Am., 67(3): 615-629.

Bakun W H. 1984. Seismic moments, local magnitudes, and coda-duration magnitudes for earthquakes in central California.Bull.Seism.Soc.Am., 74(2): 439-458.

Ben-Zion Y, Rice J R. 1993. Earthquake failure sequences along a cellular fault zone in a three-dimensional elastic solid containing asperity and nonasperity regions.J.Geophys.Res., 98(B8): 14109-14131.

Ben-Zion Y, Zhu L P. 2002. Potency-magnitude scaling relations for southern California earthquakes with 1. 0

Clinton J F, Hauksson E, Solanki K. 2006. An evaluation of the SCSN moment tensor solutions: Robustness of the Mw magnitude scale, style of faulting, and automation of the method.Bull.Seism.Soc.Am., 96(5): 1689-1705. Dreger D, Helmberger D. 1991. Source parameters of the Sierra Madre earthquake from regional and local body waves.Geophys.Res.Lett., 18(11): 2015-2018.Dreger D S, Helmberger D V. 1993. Determination of source parameters at regional distances with three-component sparse network data.J.Geophys.Res., 98(B5): 8107-8125.

Dziewonski A M, Chou T A, Woodhouse J H. 1981. Determination of earthquake source parameters from waveform data for studies of global and regional seismicity.J.Geophys.Res., 86(B4): 2825-2852.

Ebel J E, Bonjer K P. 1990. Moment tensor inversion of small earthquakes in southwestern Germany for the fault plane solution.Geophys.J.Int., 101(1): 133-146.

Fletcher J, Boatwright J, Haar L, et al. 1984. Source parameters for aftershocks of the Oroville, California, earthquake.Bull.Seism.Soc.Am., 74(4): 1101-1123.

Hanks T C, Kanamori H. 1979. A moment magnitude scale.J.Geophys.Res., 84(B5): 2348-2350.

Hanks T C, Boore D M. 1984. Moment-magnitude relations in theory and practice.J.Geophys.Res., 89(B7): 6229-6235.

Hardebeck J L, Shearer P M. 2002. A new method for determining first-motion focal mechanisms.Bull.Seism.Soc.Am., 92(6): 2264-2276.

Hardebeck J L, Shearer P M. 2003. Using S/P amplitude ratios to constrain the focal mechanisms of small earthquakes.Bull.Seism.Soc.Am., 93(6): 2434-2444.

Helmberger D V, Engen G R. 1980. Modeling the long-period body waves from shallow earthquakes at regional ranges.Bull.Seism.Soc.Am., 70(5): 1699-1714.

Julian B R, Foulger G R. 1996. Earthquake mechanisms from linear-programming inversion of seismic-wave amplitude ratios.Bull.Seism.Soc.Am., 86(4): 972-980.

Kanamori H. 1977. The energy release in great earthquakes.J.Geophys.Res., 82(20): 2981-2987.

Kennett B L N. 1983. Seismic Wave Propagation in Stratified Media. Cambridge: Cambridge University Press.

Kennett B L N, Engdahl E R. 1991. Traveltimes for global earthquake location and phase identification.Geophys.J.Int., 105(2): 429-465.

Kennett B L N, Engdahl E R, Buland R. 1995. Constraints on seismic velocities in the earth from traveltimes.Geophys.J.Int., 122(1): 108-124.

Kim W Y, Wahlstr?m R, Uski M. 1989. Regional spectral scaling relations of source parameters for earthquakes in the Baltic Shield.Tectonophysics, 166(1-3): 151-161.

Kisslinger C. 1980. Evaluation of S to P amplitude rations for determining focal mechanisms from regional network observations.Bull.Seism.Soc.Am., 70(4): 999-1014.Lay T, Wallace T C. 1995. Modern Global Seismology. New York: Academic Press.

Liang S H, Li Y M, Shu P Y, et al. 1984. On the determining of source parameters of small earthquakes by using amplitude ratios of P and S from regional network observations.ChineseJ.Geophys. (ActaGeophysicaSinica) (in Chinese), 27(3): 249-257.

Lü J, Wang X S, Su J R, et al. 2013. Hypocentral location and source mechanism of theMS7.0 Lushan earthquake sequence.ChineseJ.Geophys. (in Chinese), 56(5): 1753-1763, doi: 10.6038/cjg20130533.

Reasenberg P A, Oppenheimer D H. 1985. FPFIT, FPPLOT and FPPAGE: Fortran computer programs for calculating and displaying earthquake fault plane solutions. U. S. Geological Survey Open-File Report 85-739, 109.

Snoke J A, Munsey J W, Teague A G, et al. 1984. A program for focal mechanism determination by combined use of polarity and SV-P amplitude ratio data.JournalofEarthquakeNotes, 55(3): 15.

Tan Y, Helmberger D V. 2007. A new method for determining small earthquake source parameters using short-period P waves.Bull.Seism.Soc.Am., 97(4): 1176-1195.

Vidal A, Mungúia L. 1991. Local magnitude and source parameters for earthquakes in the Peninsular Ranges of Baja California, Mexico.Bull.Seism.Soc.Am., 81(6): 2254-2267.

Wang C Y, Mooney W D, Wang X L, et al. 2002. Study on 3-D velocity structure of crust and upper mantle in Sichuan-Yunnan region, China.ActaSeismologicaSinica(in Chinese), 24(1): 1-16.

Wang R J. 1999. A simple orthonormalization method for stable and efficient computation of Green′s functions.Bull.Seism.Soc.Am., 89(3): 733-741.

Xu L S, Du H L, Yan C, et al. 2013a. A method for determination of earthquake hypocentroid: Time-reversal imaging technique I-Principle and numerical tests.ChineseJ.Geophys. (in Chinese), 56(4): 1190-1206, doi: 10.6038/cjg20130414.Xu L S, Yan C, Zhang X, et al. 2013b. A method for determination of earthquake hypocentroid: Time-reversal imaging technique—An examination based on people-made earthquake.ChineseJ.Geophys. (in Chinese), 56(12): 4009-4027, doi: 10.6038/cjg20131207.

Xu L S, Yan C, Zhang X, et al. 2014. Rupture histories of the 2014 Yingjiang double earthquakes.ChineseJ.Geophys. (in Chinese), 57(10): 3270-3284, doi: 10.6038/cjg20141015.Xu Z H, Yan M, Zhao Z H. 1983. Evaluation of the direction of tectonic stress in north China from recorded data of a large number of small earthquakes.ActaSeismologicaSinica(in Chinese), 5(3): 268-279.

Yan C, Xu L S. 2014. An inversion technique for the mechanisms of local and regional earthquakes: generalized polarity and amplitude technique (I)—Principle and numerical tests.ChineseJ.Geophys. (in Chinese), 57(8): 2555-2572, doi: 10.6038/cjg20140816.

Yu C Q, Tao K, Cui X F, et al. 2009. P-wave first-motion focal mechanism solutions and their quality evaluation.ChineseJ.Geophys. (in Chinese), 52(5): 1402-1411, doi: 10.3969/j.issn.0001-5733.2009.05.030.

Zhao B, Gao Y, Huang Z B, et al. 2013. Double difference relocation, focal mechanism and stress inversion of LushanMS7.0 earthquake sequence.ChineseJ.Geophys. (in Chinese), 56(10): 3385-3395, doi: 10.6038/cjg20131014.Zhao L S, Helmberger D V. 1994. Source estimation from broadband regional seismograms.Bull.Seism.Soc.Am., 84(1): 91-104.Zheng X F, Ouyang B, Zhang D N, et al. 2009. Technical system construction of Data Backup Centre for China Seismograph Network and the data support to researches on the Wenchuan earthquake.ChineseJ.Geophys. (in Chinese), 52(5): 1412-1417, doi: 10.3969/j.issn.0001-5733.2009.05.031.

Zheng X F, Yao Z X, Liang J H, et al. 2010. The role played and opportunities provided by IGP DMC of China National Seismic Network in Wenchuan earthquake disaster relief and researches.Bull.Seism.Soc.Am., 100(5B): 2866-2872.

Zhu L P, Helmberger D V. 1996. Advancement in source estimation techniques using broadband regional seismograms.Bull.Seism.Soc.Am., 86(5): 1634-1641.

附中文參考文獻

陳運泰, 劉瑞豐. 2004. 地震的震級. 地震地磁觀測與研究, 25(6): 1-12.

梁尚鴻, 李幼銘, 束沛鎰等. 1984. 利用區(qū)域地震臺網(wǎng)P、S振幅比資料測定小震震源參數(shù). 地球物理學報, 27(3): 249-257.

呂堅, 王曉山, 蘇金蓉等. 2013. 蘆山7.0級地震序列的震源位置與震源機制解特征. 地球物理學報, 56(5): 1753-1763, doi: 10.6038/cjg20130533.

王椿鏞, Mooney W D, 王溪莉等. 2002. 川滇地區(qū)地殼上地幔三維速度結構研究. 地震學報, 24(1): 1-16.

許力生, 杜海林, 嚴川等. 2013a. 一種確定震源中心的方法: 逆時成像技術(一)——原理與數(shù)值實驗. 地球物理學報, 56(4): 1190-1206, doi: 10.6038/cjg20130414.

許力生, 嚴川, 張旭等. 2013b. 一種確定震源中心的方法: 逆時成像技術(二)——基于人工地震的檢驗. 地球物理學報, 56(12): 4009-4027, doi: 10.6038/cjg20131207.

許力生, 嚴川, 張旭等. 2014. 2014年盈江雙震的破裂歷史. 地球物理學報, 57(10): 3270-3284, doi: 10.6038/cjg20141015.

許忠淮, 閻明, 趙仲和. 1983. 由多個小地震推斷的華北地區(qū)構造應力場的方向. 地震學報, 5(3): 268-279.

嚴川, 許力生. 2014. 一種地方與區(qū)域地震震源機制反演技術: 廣義極性振幅技術(一)——原理與數(shù)值實驗. 地球物理學報, 57(8): 2555-2572, doi: 10.6038/cjg20140816.

俞春泉, 陶開, 崔效鋒等. 2009. 用格點嘗試法求解P波初動震源機制解及解的質量評價. 地球物理學報, 52(5): 1402-1411, doi: 10.3969/j.issn.0001-5733.2009.05.030.

趙博, 高原, 黃志斌等. 2013. 四川蘆山MS7.0地震余震序列雙差定位、震源機制及應力場反演. 地球物理學報, 56(10): 3385-3395, doi: 10.6038/cjg20131014.

鄭秀芬, 歐陽飚, 張東寧等. 2009. “國家數(shù)字測震臺網(wǎng)數(shù)據(jù)備份中心”技術系統(tǒng)建設及其對汶川大地震研究的數(shù)據(jù)支撐. 地球物理學報, 52(5): 1412-1417, doi: 10.3969/j.issn.0001-5733.2009.05.031.

(本文編輯 何燕)

An inversion technique for mechanisms of local and regional earthquakes: generalized polarity and amplitude technique (II): An application to real seismic events

YAN Chuan1，XU Li-Sheng1*，ZHANG Xu1，LI Chun-Lai1，XU Kang-Sheng2

1InstituteofGeophysics,ChinaEarthquakeAdministration,Beijing100081,China2GansuEarthquakeAdministration,Lanzhou730000,China

Determining focal mechanisms of small or micro-earthquakes is very important but difficult in seismological research, and no method has been commonly accepted and widely used due to complexity of this problem itself till now. In our previous paper we proposed an inversion technique, named generalized polarity and amplitude technique (GPAT), and the feasibility and capability of resisting disturbance of the technique has been tested with a series of numerical experiments. In this paper we present the examination of the practical ability of the GPAT by means of forty-nine real earthquakes.

10.6038/cjg20151014.

Yan C, Xu L S, Zhang X, et al. 2015. An inversion technique for mechanisms of local and regional earthquakes: generalized polarity and amplitude technique (II): An application to real seismic events.ChineseJ.Geophys. (in Chinese),58(10):3601-3614,doi:10.6038/cjg20151014.

中國地震局地球物理研究所基本業(yè)務費(DQJB14B01)和國家自然科學基金課題(41474046)聯(lián)合資助．

嚴川，男，1989年5月生，在讀博士生，主要從事震源機制與構造應力場反演研究．E-mail: yanc@cea-igp.ac.cn

*通訊作者 許力生，男，1965年3月生，研究員，主要從事地震學研究．E-mail: xuls@cea-igp.ac.cn

10.6038/cjg20151014

P315

2015-04-28，2015-06-02收修定稿

嚴川， 許力生， 張旭等. 2015. 一種地方與區(qū)域地震震源機制反演技術：廣義極性振幅技術(二)——對實際震例的應用.地球物理學報，58(10):3601-3614,