Nadia Lapusta

巖石塊體差異運(yùn)動(dòng)所造成的殼內(nèi)斷層是大部分地震的根源。在一個(gè)地震過(guò)程中，斷層的兩盤(pán)發(fā)生滑動(dòng)，或者發(fā)生斷錯(cuò)，導(dǎo)致地震破裂沿著該斷層傳播。與此同時(shí)，沿著斷層面的摩擦力會(huì)阻止滑動(dòng)。滑動(dòng)過(guò)程中摩擦力的變化是斷層行為隨時(shí)間變化的一個(gè)決定性因素，且從另一方面來(lái)看，摩擦力還影響著地震的震級(jí)大小、烈度大小和地震動(dòng)的分布范圍。人們已經(jīng)知道滑動(dòng)速率在地震過(guò)程中是有變化的，但滑動(dòng)速率是否影響以及如何影響摩擦力和斷層行為，仍是一個(gè)有爭(zhēng)議的話題。在Sone和Shimamoto的實(shí)驗(yàn)中［1］，將地震數(shù)據(jù)推斷出的滑動(dòng)速率應(yīng)用到一條天然斷層中的巖石樣本上，實(shí)驗(yàn)所得的摩擦演化的過(guò)程與典型地震模型中的摩擦行為有所不同。

斷層錯(cuò)動(dòng)是一個(gè)復(fù)雜的過(guò)程，其中有各種因素互相作用。因此，人們對(duì)地震產(chǎn)生、傳播和消減過(guò)程中的很多方面都認(rèn)識(shí)甚少。然而，如果我們要確認(rèn)哪些因素控制了地震成核和終止、地震復(fù)發(fā)，以及受地震影響地震動(dòng)的范圍，對(duì)這些方面的理解又是必要的。目前這方面的地震研究進(jìn)步顯著，人們提出了眾多模型，解決了地震破裂傳播和長(zhǎng)時(shí)間地震滑動(dòng)的問(wèn)題［2－9］，但是現(xiàn)有模型并不能解釋觀測(cè)到的所有現(xiàn)象。因此，這些模型需要利用地震、大地測(cè)量學(xué)、地質(zhì)學(xué)和實(shí)驗(yàn)室數(shù)據(jù)來(lái)不斷進(jìn)行完善。

斷層摩擦，以及摩擦在滑動(dòng)過(guò)程中的演化是地震行為模型中非常重要的部分。依據(jù)早期的理論，摩擦主要取決于滑移量。這是基于以下假定，即摩擦在達(dá)到一個(gè)固定值之前，隨著斷層滑動(dòng)量的增加而呈線性減少特征?；谶@個(gè)假定的模型仍被廣泛應(yīng)用［6］。但是，我們對(duì)斷層摩擦的理解有了進(jìn)一步發(fā)展［2，10］，實(shí)驗(yàn)室試驗(yàn)顯示，摩擦同樣也依賴(lài)于滑動(dòng)速率和滑動(dòng)過(guò)程［2，10］?？焖倩瑒?dòng)，近似于1m／s 的滑動(dòng)速率 即地震中的滑移速度 已被證實(shí)能夠?qū)е履Σ溜@著減弱。但這也可能是由于剪切生熱或是其他因素的作用所致［9，10］。

由Sone和Shimamoto主持的實(shí)驗(yàn)室試驗(yàn)［1］測(cè)試了隨著加載時(shí)間變化的滑動(dòng)速率對(duì)摩擦的作用，將我們對(duì)摩擦的理解往前推進(jìn)了一步。研究人員模擬了兩塊巖石的相對(duì)運(yùn)動(dòng)，這兩塊巖石由當(dāng)中夾著的斷裂帶物質(zhì)區(qū)分為上下盤(pán)，該斷裂帶物質(zhì)采自臺(tái)灣車(chē)籠埔斷裂破碎帶。在實(shí)驗(yàn)中，滑動(dòng)速率設(shè)置為早期隨著時(shí)間增加，后期隨著時(shí)間減少（圖1a），以重現(xiàn)1999年臺(tái)灣集集地震滑動(dòng)速率的變化。這是該試驗(yàn)的獨(dú)特之處，即采用了現(xiàn)實(shí)斷層中的物質(zhì)來(lái)進(jìn)行具有變化的滑動(dòng)速率的地震反演。

模擬結(jié)果顯示，在不同的滑動(dòng)速率作用下，摩擦系數(shù)首先增大，隨后大幅減小，最后又有所增加（圖1b）。這個(gè)摩擦演化過(guò)程與許多地震現(xiàn)象的理論和數(shù)值模擬獲得的簡(jiǎn)單摩擦行為（圖1c）非常不同。這些簡(jiǎn)單摩擦模型可能需要根據(jù)該實(shí)驗(yàn)室觀測(cè)結(jié)果進(jìn)行進(jìn)一步的修正。

研究人員認(rèn)為，他們?cè)囼?yàn)中所發(fā)現(xiàn)的摩擦系數(shù)在早期的增加可能阻礙了斷層的運(yùn)動(dòng)。這種摩擦行為可能阻止了小地震變?yōu)榇蟮卣穑部赡軐?duì)現(xiàn)實(shí)中眾多的斷層滑動(dòng)特征和地震震級(jí)的生成有一定的貢獻(xiàn)。此外，斷層摩擦的峰值似乎還依賴(lài)于滑動(dòng)速率的變化，這進(jìn)一步體現(xiàn)了斷層滑動(dòng)的復(fù)雜性。

圖1 摩擦的起伏變化。Sone和Shimamoto為了理解地震中滑動(dòng)速率的影響而設(shè)計(jì)了斷層摩擦實(shí)驗(yàn)［1］（a）；他們發(fā)現(xiàn)摩擦系數(shù)經(jīng)歷了一個(gè)早期增加，之后顯著降低，最后又有所升高的過(guò)程（b）；這種摩擦行為與地震研究中常用的摩擦行為（c）非常不同；地震過(guò)程中滑動(dòng)速率的真實(shí)變化可能要更為復(fù)雜［7－9］（d），且更難以在一個(gè)實(shí)驗(yàn)中實(shí)現(xiàn)

在這些實(shí)驗(yàn)中，摩擦系數(shù)這種早期增大，之后迅速減小，最后又增大的變化，導(dǎo)致了斷層強(qiáng)化。這種摩擦行為與地震期間短時(shí)間內(nèi)滑動(dòng)脈沖沿?cái)鄬觽鞑ツＰ鸵恢拢?－11］。地震反演［11］顯示，這種短期滑動(dòng)脈沖在真實(shí)斷裂上是存在的。對(duì)速率依賴(lài)摩擦模型來(lái)說(shuō)，若要產(chǎn)生短期滑動(dòng)，斷層僅需要相對(duì)于摩擦力峰值而言較低的平均剪切力［9，10］。由于地震僅能在局部應(yīng)力大于靜態(tài)摩擦力的情況下發(fā)生，因此該種行為造成的一個(gè)重要后果就是：地震只能在特殊部位發(fā)生。

實(shí)驗(yàn)結(jié)果有力地支持了斷層摩擦力依賴(lài)滑動(dòng)速率這一觀點(diǎn)。實(shí)驗(yàn)結(jié)果還顯示，在摩擦系數(shù)與滑動(dòng)速率之間并不存在一一對(duì)應(yīng)關(guān)系，這表明在斷層行為中還存在著其他因素的作用。正如研究人員提到的，他們的發(fā)現(xiàn)有助于定性理解所謂的狀態(tài)依賴(lài)摩擦準(zhǔn)則［2，8－10］。在該關(guān)系中，摩擦依賴(lài)于滑動(dòng)狀態(tài)，以及滑動(dòng)狀態(tài)的變化。通過(guò)修改該摩擦準(zhǔn)則的狀態(tài)變化和演化方程，能夠?qū)⒃搶?shí)驗(yàn)結(jié)果融入到數(shù)值模型中。這些數(shù)值模型在某種程度上還是與早期的實(shí)驗(yàn)工作相吻合的。

Sone和Shimamoto同時(shí)也提出了一個(gè)摩擦力、滑動(dòng)量與滑動(dòng)速率之間的數(shù)學(xué)關(guān)系［1］。但是，該數(shù)學(xué)關(guān)系在關(guān)于斷層滑動(dòng)的數(shù)值模擬中需要謹(jǐn)慎使用，因?yàn)樵摂?shù)學(xué)關(guān)系是與特殊的滑動(dòng)速率演化，以及其他實(shí)驗(yàn)中的特殊參數(shù)緊密相連，因此不是普遍適用的。尤其在一般的地震深度（譬如，加利福尼亞3～15km）處，斷層附近的擠壓力比這些實(shí)驗(yàn)中所展示的要大得多。而且在該實(shí)驗(yàn)中滑動(dòng)速率變化是基于低頻次地震資料設(shè)置的。實(shí)際上，滑動(dòng)開(kāi)始時(shí)滑動(dòng)速率的增長(zhǎng)速度要比實(shí)驗(yàn)室數(shù)據(jù)快得多，滑動(dòng)速率值也比實(shí)驗(yàn)室數(shù)據(jù)要高得多（圖1d）。

在Sone和Shimamoto的實(shí)驗(yàn)中，有幾個(gè)可能明顯影響到斷層摩擦的因素沒(méi)有考慮到［1，2，5，7，9，10］。這些因素包括：斷層破碎帶內(nèi)流體的作用；高壓下斷層物質(zhì)的熔融作用；斷層破碎帶中物性參數(shù)的變化。然而，這些因素是否明顯地影響斷層行為，以及它們?nèi)绾位ハ嘧饔?，需要更富有想象力和更精確的試驗(yàn)，通過(guò)重建地震深度上的斷層破碎帶的物性條件來(lái)進(jìn)行確認(rèn)。為了發(fā)現(xiàn)潛在的動(dòng)態(tài)現(xiàn)象和驗(yàn)證地震模型，在實(shí)驗(yàn)?zāi)Ｐ驮O(shè)置上的創(chuàng)造性仍有很大發(fā)展空間［12］，以達(dá)到重建地震破裂傳播過(guò)程的目的。這些實(shí)驗(yàn)?zāi)軌蛑苯幽M某種地震破裂特征，譬如極快的破裂速度，以及極大的滑動(dòng)脈沖［11］。

綜上所述，實(shí)驗(yàn)、理論和實(shí)際觀測(cè)定性地勾畫(huà)了地震期間的斷層行為，但還有更多的內(nèi)容需要研究。越來(lái)越多種類(lèi)的相關(guān)摩擦實(shí)驗(yàn) 譬如Sone和Shimamoto的實(shí)驗(yàn)［1］

有望帶來(lái)高質(zhì)量的觀測(cè)數(shù)據(jù)，高計(jì)算能力和精確的斷層模型，以助于更好地認(rèn)識(shí)地震現(xiàn)象。

譯自：Nature Geoscience．2009，2：676－677原題：The roller coaster of fault friction

（中國(guó)地震局地質(zhì)研究所，中國(guó)地震臺(tái)網(wǎng)中心姚琪 譯；中國(guó)地震局地質(zhì)研究所 汲云濤校）

（譯者電子信箱，姚琪：yqvoxelgeo＠163．com）

［1］Sone H，Shimamoto T．Frictional resistance of faults during accelerating and decelerating earthquake slip．Nature Geosci．，2009，2：705－708

［2］Dieterich J H．Applications of rate－and state－dependent friction to models of fault slip and earthquake occurrence．Treatise on Geophysics，2007，4：107－129

［3］Harris R A，Day S M．Dynamics of fault interaction：parallel strike－slip faults．J．Geophys．Res．，1993，98：4461－4472

［4］Olsen K B，Madariaga R，Archuleta R J．Three－dimensional dynamic simulation of the 1992Landers earthquake．Science，1997，278：834－838

［5］Andrews D J．Mechanics of fault junctions．J．Geophys．Res．，2005，110：B01307

［6］Harris R A，Barall M，Archuleta R，et al．The SCEC／USGS dynamic earthquake rupture code verification exercise．Seis．Res．Lett．，2009，80：119－126

［7］Ampuero J－P，Ben－Zion Y．Cracks，pulses and macroscopic asymmetry of dynamic rupture on a bimaterial interface with velocity－weakening friction．Geophys．J．Int．，2008，2：674－692

［8］Chen T，Lapusta N．Rate and state friction laws can explain scaling of small repeating earthquakes．J．Geophys．Res．，2009，114：B01311

［9］Noda H，Dunham E M，Rice J R．Earthquake ruptures with thermal weakening and the operation of major faults at low overall stress levels．J．Geophys．Res．，2009，114：B07302

［10］Tullis T E．Friction of rock at earthquake slip rates．Treatise on Geophysics，2007，4：131－152

［11］Heaton T H．Evidence for and implications of self－h(huán)ealing pulses of slip in earthquake rupture．Phys．Earth Planet．In．，1990，64：1－20

［12］Rosakis A J，Xia K W，Lykotrafitis G，et al．Dynamic shear rupture in frictional interfaces：speeds，directionality，and modes．Treatise on Geophysics，2007，4：153－192