震源機制解綜述

1、震源機制解綜述1 、引言地震學(xué)是一門以觀測資料為基礎(chǔ)的研究地震的成因及其規(guī)律已成為地震預(yù)報的一種重要手段， 它的發(fā)展奠定了地震預(yù)報的物理基礎(chǔ)。地震震源和地震波傳播介質(zhì)的各種參數(shù)在強震前的變化早就被當(dāng)作地震預(yù)測的地震學(xué)前兆指標(biāo)，隨著地震預(yù)測的深入研究，以及我國 “十五”臺站數(shù)字化改造的完成，我們在進(jìn)一步研究地震時空強分布特征的同時 ,加強對地震波的運動學(xué)和動力學(xué)特征的研究,從中提取震源，我們意識到加強對地震波的運動學(xué)和動力學(xué)的研究，從中提取震源信息，對增強地震預(yù)測的物理基礎(chǔ)，提高地震預(yù)測的水平是十分必要的。地震是地球內(nèi)部物質(zhì)運動的結(jié)果，這種運動反映在地殼上，使得地殼產(chǎn)生破裂，促成了斷層的生成、發(fā)

2、育和活動。地震前后的地形變測量和地震波的觀測研究等結(jié)果確認(rèn)，天然構(gòu)造地震是地下巖層的突然錯動引起的。發(fā)生錯動的巖層可稱為地震斷層。斷層活動誘發(fā)了地震， 地震發(fā)生又促成了斷層的生成與發(fā)育，因此地震與斷層有密切聯(lián)系。地殼中的斷層密如織網(wǎng)。 實際地震斷層的幾何形狀可能很復(fù)雜，但對多數(shù)地震，特別是小地震，作為初級近似，總體上可將地震看成是沿一個平面斷層發(fā)生的突然錯動引起的。2、前人對震源機制解的研究歷程地震震源處地球介質(zhì)的運動方式。通常所說的震源機制是狹義的，即專指研究構(gòu)造地震的機制而言。構(gòu)造地震的機制是震源處介質(zhì)的破裂和錯動。震源機制研究的內(nèi)容包括，確定地震斷層面的方位和巖體的錯動方向，研究震源處巖

3、體的破裂和運動特征，以及這些特征和震源所輻射的地震波之間的關(guān)系。對地震震源的研究開始于20 世紀(jì)初葉。1910 年提出的彈性回跳理論，首次明確表述了地震斷層成因的概念。在地震學(xué)的早期研究中，人們就已注意到P波到達(dá)時地面的初始振動有時是向上的，有時是向下的。20世紀(jì)的1020年代，許多地震學(xué)者在日本和歐洲的部分地區(qū)幾乎同時發(fā)現(xiàn)，同一次地震在不同地點的臺站記錄,所得的 P 波初動方向具有四象限分布。日本的中野廣最早提出了震源的單力偶力系，第一次把斷層的彈性回跳理論和P 波初動的四象限分布聯(lián)系起來。此后，本多弘吉又提出雙力偶力系，事實證明它比單力偶力系更接近實際。美國的拜爾利（PEyerly）發(fā)展了

4、最初的震源機制求解法,1938年第一次利用 P波初動求出完整的地震斷層面解。3、斷層及斷層面參數(shù)3.1、 斷層參數(shù)及分類地震斷層通常用斷層的走向。S傾角8和滑動角入三個參數(shù)來描述（圖 2.1）。按目前國際上常用的描述方法，這些參數(shù)的定義是：走向（f）S:斷層面與水平面交線的方向，但此交線有兩個方向，為唯一確定起見，按以下原則確定其中之一為斷層的走向：人沿走向看去，斷層上盤在右。走向用從正北順時針量至走向方向的角度 。S來表示，0o<（）S<360°。傾角6 :斷層面與水平面的夾角。0o<8 W90°。滑動角入：在斷層面上量度，從走向方向逆時針量至滑動方向的

5、角度為正，順時針量至滑動方向的角度為負(fù)。滑動方向指斷層上盤相對于下盤的運動方向。-180入<180° o（仰角：力軸與水平面的夾角（小于90 度）方位角：力軸在水平面上的投影線與北方向之間的夾角傾向：節(jié)面的上表面的法線在水平面上的投影線與北方向之間的夾角，順時針量取。）走向jS和傾角8是斷層的幾何參數(shù)，二者規(guī)定了斷層的產(chǎn)狀； 滑動角入是斷層的運動參數(shù)， 由這一參數(shù)的具體數(shù)值，即可描述斷層的各種運動類型（圖2.2） 。0-9090-180左旋逆 右旋逆 凌斷層 沖斷層 左宸連 右旋送斷層 斷層180-270 270-360右旋正左旋走甌層滑斷層090迂沖F層主友走滑-r孱1802

6、70右簧走正斷層 洎L蕓有人用斷層的傾向代替走向，傾向指下盤斷層面向上的法線之水平投影的方向，傾向恒等于走向加90。在地震學(xué)中，通常已較少用傾向描述地震斷層。左旋走涓（0度）圖21 斷層的參數(shù)表示圖2.2滑動角入取不同數(shù)值所描述的斷層類型3.2、 斷層面上的錯動斷層滑動是時間和空間的函數(shù)，斷層面上的錯動主要是平行于斷層面的剪切位錯，描述這些錯動的參數(shù)如下:(1)地震矩：將地震看成斷層面上的突然位錯，則形成地震力矩，定義 0M如下(2)地震能量：TE(3)應(yīng)力降：(T?(4)破裂速度：v3.3、 震源模型圖立3震源坐標(biāo)系芯-y-z和觀測坐標(biāo)系r- 9 -。震源機制解(又稱地震機理)是 指震源區(qū)地

7、震發(fā)生時的力學(xué)過程。鑒于地震機理的研究尚處于探索階段， 目前還屬于推斷性認(rèn)識， 一般采用各 種震源模型進(jìn)行解析，一種是點源模 型，另一種是非點源模型。前者根據(jù) 點源作用力的不同，又進(jìn)一步劃分為 單力偶震源模型和雙力偶震源模型； 非點源模型也劃分為有限移動震源 模型和位錯震源模型兩種。以上震源 模型，在分析求解后，提供兩組力學(xué) 參數(shù)，一組為斷層面走向、傾向和傾 角；另一組為最大主應(yīng)力軸、最小主 應(yīng)力軸和中等主應(yīng)力軸的方位和產(chǎn) 狀。根據(jù)我國境內(nèi)150次地震震源機 制解，P波初動符號資料確定結(jié)果表 明，大多數(shù)主壓應(yīng)力軸(P軸)和主 張應(yīng)力軸都近水平 (T軸)，中等應(yīng)力 軸近于直立。地震學(xué)的震源理論證

8、明，在均勻彈性介質(zhì)中,若在一個小的平面斷層上發(fā)生一個突然的純剪切錯動,則會產(chǎn)生地震波輻射，這樣的剪切錯動震源產(chǎn)生的遠(yuǎn)場地震波與在震源處突然有一個雙力偶的作用產(chǎn)生的地震波相同。即剪切元位錯震源與雙力偶點源在產(chǎn)生遠(yuǎn)場地震波的意義上是等價的。因此，當(dāng)可將震源近似看成點源時，雙力偶點源模型就成為描述發(fā)生了剪切錯動震源的常用模型。力偶由一對大小相等、 方向相反的力偶組成，（圖2.3）是一種合力和合力矩都等于零的集中 力系。這樣的力系作用于剛體時，不會產(chǎn)生任何運動效果，但在彈性體內(nèi)部作用，則會使震源區(qū)介質(zhì)產(chǎn)生突然的變形，從而向外輻射地震波。4、震源機制解的測定利用雙力偶點源模型，根據(jù)地震波觀測（或地震前后

9、的地形變測量資料等）求震源模型參數(shù)的結(jié)果，通常稱為震源機制解答，有人稱作地震的斷層面解。所根據(jù)的觀測資料可以是 P波的初動方向、S波位移的偏振方向、直達(dá) P波和直達(dá)S波 振幅的比值大小，以及 P波和S波的波形資料等。雙力偶點源模型的獨立模型參數(shù)只有3個，例如可以是斷層面（P波兩個節(jié)面中的一個）的走向、傾角和滑動角（圖 2.1）,也可以 是為確定“震源坐標(biāo)架" x-y-z （圖2.3）相對于“地平坐標(biāo)架”（例如可分別選為北、東、下 三個方向）的空間方位所需要的三zyx-個角度值。求解的基本方法是先假定震源模型參數(shù)，計算出在給定地球地震波速度結(jié)構(gòu)時，該震源模型在各觀測臺站所產(chǎn)生的地震波特

10、征，然后與各個臺站的實際觀測地震波資料進(jìn)行對比， 二者擬合最好的模型參數(shù)就作為震源機制的解答。求解過程可以運用反演的數(shù)學(xué)方法來實現(xiàn)，即選定一種使各個臺站的計算結(jié)果與觀測結(jié)果互相擬合的準(zhǔn)則，然后使隨模型參數(shù)變化而變化的準(zhǔn)則函數(shù)（或目標(biāo)函數(shù)）最優(yōu)化，即使其最小或最大，而解出模型參數(shù)來。4.1 P波初動方向法求震源機制解答最簡單的方法是根據(jù)P波初動方向的觀測資料來求解。圖4.1臺站號的觀測量要歸算到震源球面上相應(yīng)的位置3上去地表垂直向地震儀記錄的初至P波的振動方向，有的向上，即壓縮波，記為正號，有的向下，即拉伸波，記為負(fù)號。由于介質(zhì)速度結(jié)構(gòu)的影響，從震源發(fā)出的P波一般不是沿直線到達(dá)每個臺站 Si（i

11、=1,2,）的，如圖3.1所示意表達(dá)的。求震源機制解時，需根據(jù)已知的速 度結(jié)構(gòu)推算出到達(dá)每個臺站Si的P波從震源處是沿什么方位Si'發(fā)出的，即需要將臺站 Si的記錄標(biāo)在震源球面的相應(yīng)位置Si'上去。震源球面是包圍震源的一個球面S,要求球面內(nèi)的射線不再發(fā)生任何彎曲。若將每個臺站Si所記錄到的P波初動方向都標(biāo)在震源球面上的相應(yīng)位置Si'上去后，人們發(fā)現(xiàn)，對于天然構(gòu)造地震，只要記錄足夠多，并且 Si'在球面上 的分布范圍足夠廣，則可以找到過球面中心的兩個互相垂直的平面，將震源球面上的正、負(fù)號分成4個象限，這兩個平面就是上述的雙力偶震源的兩個節(jié)平面。找到兩個節(jié)平面的空間

12、位置后，震源坐標(biāo)架的 x、v、z軸和P、T軸的空間方位也就知道了。上述求解過程可以通過計算機來實現(xiàn)。P注5波ffi 4.2 雙力偶點源沿正x方向移動速度為S波速度的一半時, 力偶平面內(nèi)P波和$波的舞射花樣兩個節(jié)面中有一個是斷層面，但僅卞!據(jù)P波初動方向記錄無法確定哪一個是斷層面，還必須根據(jù)其他資料，例如現(xiàn)場地質(zhì)考察資料、余震的空間分布、地震波輻射輻射花樣的不對 稱性（圖3.2）和地震波輻射的多普勒效應(yīng)等，來從兩個節(jié)面中分析判定實際的斷層面，這 種判定一般只對大地震才能實現(xiàn)。4.2 P波和S波振幅比方法利用P波初動方向記錄反演震源機制解雖簡單易行，但也有缺陷。為能將地震波節(jié)平面的空間位置約束住，

13、 最好要有緊靠節(jié)面位置的初動符號觀測數(shù)據(jù)，然而，由于P波輻射花樣的固有特征，愈靠近節(jié)面，P波愈弱，初動方向愈不易辨認(rèn)。再有，由于地表臺站布局的限制，觀測數(shù)據(jù)點在震源球面上的覆蓋范圍經(jīng)常難以令人滿意，當(dāng)在震源球面上某位置有零散的初動方向數(shù)據(jù)時，它只能告訴你該點應(yīng)處在P波的正象限或負(fù)象限，但不能獲得節(jié)面與此觀測點間的角距離有多大的信息（震源球面上某點觀測到的波的振幅大小或波形資料，則含有節(jié)面離該點角距離有多遠(yuǎn)的信息），因而零散分布的初動方向數(shù)據(jù)對解答起不了多大的約 束作用。此外，經(jīng)常有矛盾的P波初動方向讀數(shù)。上述這些因素常引起解答的不確定性。為解決這些問題，對于較大的遠(yuǎn)地震，目前多利用在不同方位且

14、具有不同震中距的臺站 上獲得的地震波形記錄來反演震源機制解。對地方性的中小地震，有人提出利用從震源向上射出的直達(dá)P波（Pg）和直達(dá)S波（Sg引起的地動位移振幅比求解震源機制的方法5-7。由（2.1）（2.3）式可見，利用直達(dá)P波（）和 是利用它們的輻射花樣的比值來求解震源機制參數(shù)。 單種波的輻射花樣要強烈得多。在力偶平面內(nèi)（圖 此意義來說，只要有正確的直達(dá)波的觀測振幅比，SV波（）或SH（）的振幅ru。u。u比，實際振幅比的輻射花樣隨空間方位的變化比3.2）振幅比的輻射花樣呈8瓣分布。從且觀測值歸算到震源球面上后的位置是正確的話，振幅比觀測對震源機制參數(shù)有較強的約束能力。（但是，振幅比的測量誤

15、差和射線在震源球面上出頭點位置的誤差也容易引起地震節(jié)面解的大誤差。）此外，對于近震，可近似認(rèn)為儀器對直達(dá)P波和S波的頻率響應(yīng)是相同的，求振幅比時可以消去儀器影響。但在實用中，也存在一些困難。 主要困難是直達(dá)S波的識別和結(jié)構(gòu)影響的校正。求震源機制解的理論模型分析的皆是直達(dá)P波和直達(dá)S波的振幅比，對于近震，確定初至S波振幅常常比較困難，特別要使用垂直向的SV波初動振幅時5,6,測量更困難些。觀測量是地動位移，為求震源機制解需要將地動位移的振幅比校正為入射波位移的振幅比。例如，Kisslinger等（1981）的程序5,用的是地方地震的P波和SV波垂直向地動位移的振幅比（USV/UP）乙由于地震較近

16、，用的是上行直達(dá)波引起的地動位移，回避了地殼內(nèi)結(jié)構(gòu) 對波的影響，只考慮了自由表面的作用。若令垂直向地動位移振幅比（USV/UP）Z與入射波垂直向位移振幅比（uSV/uP）Z的比值為G,圖（2.7）給出了G隨入射波入射角的變化曲線。觀測 的振幅比需除以G后才能得到入射波的振幅比。由圖（2.7）可見，當(dāng)入射波的入射角為30多度時，G值變化劇烈，很難從觀測的振幅比求出穩(wěn)定的入射波振幅比。因為大約35。多的角度是SV波入射的臨界角（大于此角度后，入射SV在自由面上反射的P波出射角變?yōu)?0°,質(zhì)點呈橢圓振動，振幅隨深度衰減），反射P波有相移，由入射SM反射SV和反射P波三者合成的地表位移振動變

17、得比較復(fù)雜。梁尚鴻等（1984）6的程序用的觀測量也是垂直向直達(dá)SV和P波地動位移振幅比，他們用計算層狀介質(zhì)理論地震圖的辦法來考慮介質(zhì)結(jié)構(gòu)的影響， 其中也包含自由表面的影響。需要選擇符合實際的層狀結(jié)構(gòu)才能有效消除結(jié)構(gòu)的影響，否則模型結(jié)構(gòu)的誤差會影響機制參數(shù)的正確測定。當(dāng)只用振幅比的大?。ㄕ穹鹊慕^對值）測定震源機制解時5,6,只能求出兩個地震節(jié)面的空間位置，而不能確定可能斷層面的運動特性（旋性），或不能確定P波的壓縮和膨脹的象限，也即不能確定P、T軸的具體方位。為確定后者，還必須至少要知道一個P波初動方向的可靠讀數(shù)。如果程序中考慮了觀測P波和S波的初動方向?qū)φ穹热≌?fù)值的影響了， 則可以求出

18、完整的震源機制參數(shù)來。圖4.3自由面對垂直向地動位移振幅比的影響因子隨直達(dá)波入射角的變化。U代表地表位移，u代表入射波位移。有人發(fā)展只用垂直向 SV和P的振幅比測定震源機制的方法5,6,最初是因為單分向的垂直向記錄最易得到。如果能獲得三分向記錄，使用直達(dá)SH和P波的振幅比求解震源機制口將會減少結(jié)構(gòu)和自由面對結(jié)果的影響，因為結(jié)構(gòu)界面對SH波的影響簡單，而 SH入射到自由面時，不產(chǎn)生其它轉(zhuǎn)換波， 地動位移振幅就簡單地是入射波位移振幅的2倍。如果使用地動速度的數(shù)字記錄，而理論分析是用地動位移的振幅比，最好要將速度記錄積分成地動位移再求地動位移振幅比，因為一般情況下，地動速度振幅比不一定等于地動位移振

19、幅比。4.3 其他方法震源輻射地震波的許多特征都可以用來測定震源機制，例如，利用直達(dá)P波或S波振幅的大小8,或在頻率域分析直達(dá) P波或S波振幅譜的低頻幅值9。利用數(shù)字地震記錄，這 些量可以觀測得比用筆繪記錄更準(zhǔn)確些。對于中小地震，利用近距離臺站時，對觀測振幅或振幅譜一般只需作相對簡單的自由面、幾何擴散和衰減的校正，而不別作復(fù)雜的理論地震圖的計算。4.4 震源機制參數(shù)的表達(dá)震源機制解答通常是通過給出震源參數(shù)和用圖示來表達(dá)。作為一例，表3.1給出2001年11月14日中國昆侖山 MS8.1地震的震源機制解的參數(shù)表達(dá)（哈佛大學(xué)測定）。力軸的方位（Az）從正北順時針量度，傾角（PL）為力軸從水平面向下

20、彳的角度。表中的12個參數(shù)中，只有3個是獨立的，根據(jù)其中的任意3個，例如節(jié)面I或節(jié)面H的3個參數(shù)，或是 P軸的2個參數(shù)加T軸的一個參數(shù)等，可以計算出其他參數(shù)來。5個力軸的空間位置用方位角和仰角兩個參數(shù)來表示：仰角：力軸與水平面的夾角（小于 90度）方位角：力軸在水平面上的投影線與北方向之間的夾角表4.12001年11月14日中國昆侖山M&&1地震.的震源機制解答芍 而 TSL-S_H P軸 T -N軸" s 5 X s B X Az PL Az PL Az PL96° 54°-7D190° 84D -144059° 29

21、6;318° 20°1980 54°小s走向，5傾角，入一滑動角；Az一方位角，PL一傾角*將震源球面上表示的數(shù)據(jù)和震源機制解結(jié)果用圖示法表達(dá)出來時，需要借助于某種將球面上的點與平面上的點一一對應(yīng)起來的投影法。常用的有烏爾夫網(wǎng)投影和施密特網(wǎng)投影法。二方法的投影平面通常皆是震源球面過球心的水平大圓面（也可以是其它過球心的某個大圓面）（圖3.4a）,投影點為與大圓面相對的極點。一般只用半個球面的投影，震源機制解通常用下半震源球面的投影，這時投影點為震源球面的上極點zo若某射線P穿過下半震源球面的Q點到達(dá)臺站，則其在投影圖上的投影點為Q如果某射線從震源向上發(fā)出，穿過震源

22、球面上的R點到達(dá)臺站 St （圖3.4a）,可將射線經(jīng)震源 O向反方向延伸至震源球面上 R的對 庶點R1,用R1在投影平面上的投影 R1 '來代R的投影。從雙力偶點源地震波輻射的空間 對稱性考慮，這樣做是可以的。5(b) 兩 圖4.4震源球面及其投影網(wǎng)。球面上的Q點的空間方位是由過該點矢徑的方位角Az和離源角ih確定的，Q在投影網(wǎng)上的投影點Q '的位置也是通過這兩個角度來確定的(圖2.8b), Az從網(wǎng)的正北標(biāo)記順時針量(0360° ),與ih對應(yīng)的線段OQ'的長度d根據(jù)投影規(guī)則定。 對烏爾夫投影，如圖2.8c 所示，將大圓的半徑長度看成為 1,d = tan(ih/2)3;對施密特投影，該長度應(yīng)為d=Sin(ih/2)4。 烏爾夫投影是一種等角投影，球面上曲2線的交角投影到平面上后保持不變，球面上的圓投影到平面上后仍是一個圓。該投影關(guān)系簡單，但在表示球面上的圖形時，網(wǎng)心部分常過于集中，邊緣部分則較稀疏