張?jiān)鰮Q，馬曉靜，王 時(shí)

河北地質(zhì)大學(xué) a. 研究生學(xué)院、b. 資源學(xué)院，河北 石家莊 050031

拐折斷層地震活動(dòng)性研究進(jìn)展

張?jiān)鰮Qa，馬曉靜b，王 時(shí)a

河北地質(zhì)大學(xué) a. 研究生學(xué)院、b. 資源學(xué)院，河北 石家莊 050031

論文從拐折斷層地震活動(dòng)性特征、斷層黏滑運(yùn)動(dòng)的巖石物理實(shí)驗(yàn)和數(shù)值模擬實(shí)驗(yàn)角度分析總結(jié)了拐折斷層的發(fā)震特征與發(fā)震機(jī)理。對(duì)實(shí)際震例的監(jiān)測與野外觀測結(jié)果表明：拐折斷層的角度和拐點(diǎn)對(duì)地震的發(fā)生、震級(jí)的大小、斷層的破裂方向有重要影響，而室內(nèi)的實(shí)驗(yàn)研究表明斷層帶內(nèi)的物質(zhì)成分、斷層的摩擦系數(shù)、應(yīng)力擾動(dòng)等因素對(duì)斷層黏滑運(yùn)動(dòng)特征有很大影響。作為自然界普遍存在的斷層形式，拐折斷層的黏滑運(yùn)動(dòng)特征及其影響因素的研究有助于我們了解地震的發(fā)生機(jī)理，在未來的研究中應(yīng)進(jìn)一步加強(qiáng)。

拐折斷層；地震活動(dòng)性；黏滑運(yùn)動(dòng)；數(shù)值模擬

一、引言

地震會(huì)給人民的生命財(cái)產(chǎn)造成重大損失，而淺源地震的危害尤為嚴(yán)重。目前，地震預(yù)測模型及其形成、發(fā)展的標(biāo)志性可定量要素尚未找到[1]。而地震發(fā)震機(jī)理的研究是地震預(yù)測的前提，斷層的黏滑運(yùn)動(dòng)作為大震發(fā)生的一種機(jī)制已被越來越多的人接受。然而自然界的斷層面很少以簡單連續(xù)的面狀結(jié)構(gòu)出現(xiàn)，通常都表現(xiàn)為復(fù)雜的斷層帶，大多是不連續(xù)結(jié)構(gòu)，例如拐折斷層、雁列式斷層、階步式斷層、分支斷層等[2]。這些特殊的構(gòu)造形式對(duì)斷層應(yīng)力集中、應(yīng)變分配、斷層滑動(dòng)以及地震孕育發(fā)生及破裂的傳播過程有重要影響。許多震例表明斷層的拐折部位對(duì)地震活動(dòng)起著至關(guān)重要的作用，因此對(duì)拐折斷層黏滑運(yùn)動(dòng)過程的研究具有重要的理論意義和現(xiàn)實(shí)意義。

本文首先總結(jié)了拐折斷層的地震活動(dòng)特征與導(dǎo)致其發(fā)生的原因，對(duì)發(fā)震機(jī)理——斷層的黏滑運(yùn)動(dòng)相關(guān)研究進(jìn)行了分析。闡述了拐折斷層的地震活動(dòng)特性、拐折斷層的地震前兆位置、拐折斷層的破裂特征特征、斷層黏滑活動(dòng)的影響因素。同時(shí)也分析了物理研究方法與數(shù)值模擬方法的優(yōu)勢與局限性。并建議未來研究方向應(yīng)為有更適用地學(xué)研究的物性參數(shù)并涉及多物理場的數(shù)值實(shí)驗(yàn)，以及有更先進(jìn)儀器設(shè)備和更合理的觀測手段的物理實(shí)驗(yàn)室研究，以便為強(qiáng)震發(fā)生的機(jī)制提供更多的理論基礎(chǔ)。

二、拐折斷層的地震活動(dòng)性特征

地震的發(fā)生大多是先存斷裂的活化，實(shí)際震例表明，許多強(qiáng)震大都發(fā)生在拐折斷層上。例如，在San Andreas大斷裂分別于1690年、1630年和1857年發(fā)生的三次地震均位于大拐折的附近；1975年、1976年分別發(fā)生于土耳其的Lice與Caldiran地震和1976年發(fā)生在唐山大地震同樣為拐折斷層的活動(dòng)；而1989年發(fā)生在Loma Prieta的地震和1994年發(fā)生在Lander的地震以及1999年發(fā)生在Izmit的地震也發(fā)生在斷層的拐折附近[3-4]；類似的是2002年發(fā)生在Denali斷層上的地震和1976年發(fā)生的危地馬拉地震亦發(fā)生于拐折斷層活動(dòng)帶上；2008的汶川地震也與拐折斷層相關(guān)。

（一）拐折斷層地震特征

Kato[5]的實(shí)驗(yàn)結(jié)果發(fā)現(xiàn)拐折斷層的動(dòng)力學(xué)破裂屬于“雙震型”，而事實(shí)上確實(shí)如此，例如，2006年發(fā)生在日鹽津縣豆沙鎮(zhèn)的兩次5.1級(jí)地震，兩次地震均位于NE向斷裂的同一構(gòu)造區(qū)，這次雙震型地震被認(rèn)為是由于障礙體存在而導(dǎo)致的[6]。云龍[2]的實(shí)驗(yàn)結(jié)果表明單拐折斷層的失穩(wěn)事件往往是時(shí)間相隔很短的雙震事件，而壓性雙拐折斷層的失穩(wěn)事件往往由時(shí)間間隔較長的兩次事件組成。Kato[5]提出拐折斷層的兩次破裂受主壓應(yīng)力夾角的控制，低角度破裂早于高角度，且始于低角度的破裂終止于拐點(diǎn)，而始于高角度的破裂則貫穿整個(gè)斷層，也可以說拐折部位正應(yīng)力場的變化控制著破裂的起點(diǎn)和終點(diǎn)。而初始剪應(yīng)力、正應(yīng)力的分布影響著拐點(diǎn)處的破裂速度[7]。然而，雙震的發(fā)生時(shí)間間隔因?qū)嶒?yàn)中應(yīng)變與觀測手段等的差異并不相同，例如kato等[5]的結(jié)論是10.8ms-3.5s，而郭彥雙等[8]的實(shí)驗(yàn)結(jié)果則是100ms-200ms。此外，馬瑾[9]通過拐折斷層的溫度場實(shí)驗(yàn)研究也同樣觀測到雙震甚至多震而引起的溫度場變化，并提出下斷層先錯(cuò)動(dòng)，上斷層后錯(cuò)動(dòng)，兩者的時(shí)間間隔為0.92s，拐折角度、巖石厚度、巖石類型等條件的不同都會(huì)對(duì)雙震的時(shí)間間隔有影響。

“雙震型”這一特點(diǎn)的原因或許可以歸為拐點(diǎn)的控制作用，而Kato[5]認(rèn)為速度狀態(tài)基礎(chǔ)可能是破裂延遲的原因。郭伶俐[10]研究了三種典型斷層的失穩(wěn)瞬態(tài)過程，認(rèn)為拐折斷層的拐點(diǎn)在斷層的活動(dòng)中起到類似樞紐的作用，這一點(diǎn)與馬瑾[11]的實(shí)驗(yàn)結(jié)論是相同的。正是因?yàn)楣拯c(diǎn)的控制，導(dǎo)致區(qū)域的應(yīng)變積累和積累后的急劇釋放，引起多點(diǎn)錯(cuò)動(dòng)，從而致使拐折斷層多為雙震或者多震。此外，這種斷層錯(cuò)動(dòng)勢必造成斷層幾何發(fā)生變化，甚至導(dǎo)致新斷層的產(chǎn)生及引起區(qū)域應(yīng)力場的重新分布，進(jìn)而使得斷層運(yùn)動(dòng)形式發(fā)生新的變化。Andrews[12]認(rèn)為斷層的拐折部位對(duì)斷層錯(cuò)動(dòng)具有抑制作用，致使拐點(diǎn)兩側(cè)斷層的位移量產(chǎn)生差異，當(dāng)走向上的變化劇烈時(shí)，將會(huì)產(chǎn)生新的斷層，這應(yīng)該是拐折斷層多發(fā)育強(qiáng)震的原因。

（二）拐折斷層的破裂過程

Kase等[7]根據(jù)拐折部位能否抑制斷層的錯(cuò)動(dòng)，將拐折體分為抑制型拐折(Restraining bend)和松弛型拐折(Releasing bend)，并提出拐折對(duì)破裂的傳播有非常重要的作用，拐折附近的位錯(cuò)速度受破裂導(dǎo)致的應(yīng)力變化控制。在San Andreas的Parkf i eld段上的多次地震的震中都在一個(gè)拐折斷層的拐點(diǎn)附近，這些地震的前震分布在拐點(diǎn)的西北段，主震則開始于拐點(diǎn)的附近，而余震位于轉(zhuǎn)折部位的南東段[13]。Custódio[14]通過對(duì)在San Andreas的Parkf i eld段上1966年和2004年的地震的研究指出兩次地震破裂的傳播方向是相反的。馬瑾[11]通過對(duì)5°拐折斷層的物理實(shí)驗(yàn)和數(shù)值模擬結(jié)果表明：大事件并未發(fā)生在拐折位置而是發(fā)生在兩側(cè)的平直斷層位置，并認(rèn)為拐點(diǎn)起到閥門似的作用?？梢?，有些地震的震中位于拐折附近，破裂向拐折的一側(cè)或同時(shí)向兩側(cè)傳播。但也有些地震震中位于拐折兩側(cè)的平直斷層上，這可能與拐折處障礙體的不完全破裂有關(guān)，當(dāng)障礙體全部破裂后大震才會(huì)發(fā)生，這與馬瑾[9]通過實(shí)驗(yàn)研究5°拐折斷層溫度場演化得出的結(jié)論是一致的。

（三）拐折斷層的角度

與強(qiáng)震相關(guān)的拐折斷層多為低角度(5°～30°)，例如，1992年發(fā)生在美國Landers的大地震圣安德列斯斷層的拐角大約為18°，而與唐山大地震相關(guān)的郯廬大斷裂不同段的拐折角度最大也只有25°，而更多與強(qiáng)震相關(guān)的拐折斷層為接近平直的小角度拐折[11]。King等[15]通過對(duì)大量拐折斷層帶的研究，提出很多強(qiáng)震存在于兩個(gè)5°拐折斷層之間。目前大多學(xué)者認(rèn)為在斷層帶上的5°拐折區(qū)域?qū)φ鹪次恢?、斷層破裂傳播過程、前震和余震分布及位錯(cuò)分布有重要影響，所以許多實(shí)驗(yàn)室研究采用了這一角度。而強(qiáng)震多發(fā)育于低角度拐折斷層的原因可能是因?yàn)榈徒嵌葘?duì)破裂的傳播阻礙較小，斷層間更容易發(fā)生較強(qiáng)烈的相互作用。然而，不同拐折角度的對(duì)比研究并沒有學(xué)者分析。

（四）拐折斷層的地震前兆

野外觀測結(jié)果顯示斷層拐折處通常破裂非常復(fù)雜，拐折斷層段上在一定范圍內(nèi)位移量大小與離拐點(diǎn)的距離成反比[12]，所以通常認(rèn)為拐點(diǎn)處的前兆觀測非常重要。馬瑾[16]在研究不同構(gòu)造部位在變形中的作用時(shí)，指出前兆并不一定出現(xiàn)在發(fā)震區(qū)域，會(huì)出現(xiàn)前兆偏離的現(xiàn)象。斷層幾何結(jié)構(gòu)對(duì)地震前兆有重要影響，幾何結(jié)構(gòu)愈復(fù)雜，參與失穩(wěn)過程的構(gòu)造部位愈多，則前兆種類愈多，現(xiàn)象愈復(fù)雜。而其他學(xué)者在對(duì)比研究平直斷層與幾何彎曲斷層時(shí)也得出了相似的結(jié)論[17]。郭彥雙等[8]、云龍[2]則提出觀測的方式方法不同所得出前兆響應(yīng)也存在差異。馬瑾[9]以5°拐折斷層一次粘滑事件中不同變形階段的溫度場演化為例,分析了進(jìn)入亞失穩(wěn)狀態(tài)時(shí)應(yīng)力隨時(shí)間的變化狀態(tài)等問題，提出亞失穩(wěn)到失穩(wěn)階段占整個(gè)時(shí)間比例不到1%。

三、斷層黏滑運(yùn)動(dòng)的實(shí)驗(yàn)室研究

一些學(xué)者在斷層黏滑運(yùn)動(dòng)的影響因素方面做了很多探討，例如斷層間物質(zhì)對(duì)黏滑運(yùn)動(dòng)的影響，不同滑動(dòng)方式會(huì)在實(shí)驗(yàn)產(chǎn)物中留下不同的顯微結(jié)構(gòu)，呂廣廷[18]等認(rèn)為斷層泥的均勻性與非均勻性是斷層發(fā)生蠕滑和黏滑的直接微觀證據(jù)；斷層泥中水的含量對(duì)斷層黏滑運(yùn)動(dòng)有一定的影響，化學(xué)角度上水的存在增加了其穩(wěn)定性，而力學(xué)角度上水的存在則起到了相反的作用[19]。近年來，隨著相關(guān)研究的不斷進(jìn)行，斷層滑動(dòng)失穩(wěn)的成核、均勻與非均勻斷層黏滑運(yùn)動(dòng)中應(yīng)變以及斷層位錯(cuò)的實(shí)驗(yàn)研究等也為許多學(xué)者所關(guān)注[20]。有的學(xué)者通過實(shí)驗(yàn)研究認(rèn)為斷層黏滑的動(dòng)態(tài)活動(dòng)時(shí)間很短且運(yùn)動(dòng)由高速與低速交錯(cuò)組成[21]，并且加載速率與粘滑事件的周期長短呈負(fù)相關(guān)，而與應(yīng)力降呈正相關(guān)[8]。黃元敏等[22]認(rèn)為剪應(yīng)力擾動(dòng)是影響斷層黏滑運(yùn)動(dòng)的推動(dòng)力，而正應(yīng)力擾動(dòng)則導(dǎo)致斷層接觸狀態(tài)的改變，所以正應(yīng)力擾動(dòng)較剪應(yīng)力對(duì)斷層黏滑運(yùn)動(dòng)的影響大?？姲Ⅺ怺23]以硬石膏和巖鹽等作為斷層泥材料，通過雙軸摩擦實(shí)驗(yàn)研究了摩擦滑動(dòng)的速度依賴性轉(zhuǎn)換，解析了斷層物質(zhì)的微觀結(jié)構(gòu)及相關(guān)地震機(jī)制意義。上述的物理實(shí)驗(yàn)室研究雖不完全是以拐折斷層為模型，但是其影響因素對(duì)拐折斷層黏滑運(yùn)動(dòng)的影響也是類似的。

大多數(shù)物理實(shí)驗(yàn)室研究是通過巖石實(shí)驗(yàn)的方法，所選用的巖石多為上地殼分布最廣的花崗巖類，使用伺服壓機(jī)等對(duì)樣本進(jìn)行實(shí)驗(yàn)并利用聲發(fā)射、應(yīng)變片等物理手段進(jìn)行研究。關(guān)于拐折斷層的研究多為某一單拐折或雙拐折構(gòu)造在一次粘滑事件中變形過程及相應(yīng)的物理場。本文對(duì)比了拐折斷層主要的物理實(shí)驗(yàn)研究，如表1所示。然而，目前的許多研究是某個(gè)粘滑過程中各種物理場的演化過程，缺乏與多個(gè)粘滑過程引起的應(yīng)變場的差異的對(duì)比研究，且缺乏不同拐折角度的對(duì)比研究。此外，實(shí)驗(yàn)研究與野外實(shí)測資料的對(duì)比分析也是很有必要的。

表1 部分典型拐折斷層物理實(shí)驗(yàn)類比表

四、斷層黏滑運(yùn)動(dòng)的數(shù)值實(shí)驗(yàn)研究

邊界元法、有限元法、有限差分法等是主要的數(shù)值模擬方法。邊界元方法廣泛應(yīng)用于各種斷層的黏滑過程模擬研究[24-25]，它可以采用不同摩擦定律，也可以引入斷層的非均質(zhì)性。而有限元通常被認(rèn)為是適用更加廣泛的一種方法[26]。此外，拐折斷層的處理方法主要有兩種，一種是把整個(gè)模型處理成一個(gè)連續(xù)體，把斷層設(shè)置為軟弱地帶；另一種是把整個(gè)幾何模型設(shè)置為2個(gè)接觸的模型，通過接觸分析計(jì)算斷層間的非線性關(guān)系。

事實(shí)上，數(shù)值模擬的方法正在越來越多的應(yīng)用在地震危險(xiǎn)性與預(yù)測的相關(guān)研究中。Sasani[27]通過對(duì)區(qū)域應(yīng)力分布的計(jì)算來評(píng)估區(qū)域地震危險(xiǎn)潛力。馬曉靜[28]通過數(shù)值模擬研究了拉伸速率、摩擦系數(shù)、正斷層傾角對(duì)斷層黏滑運(yùn)動(dòng)的影響，大震發(fā)震間隔隨拉伸速率（＜5mm/a時(shí)）增大而變短，摩擦系數(shù)越?。ǎ?.4時(shí)），發(fā)震間隔越短，傾角大約60°時(shí)，斷層發(fā)震最強(qiáng)烈。姜輝[29]模擬研究了太平洋板塊俯沖帶上的斷層黏滑運(yùn)動(dòng)過程，分析了俯沖帶地震的發(fā)生模式以及俯沖角度的影響，有角度變化的模型發(fā)生的黏滑事件較多。Cruzatienza等[30]利用邊界積分方程法（BIEM）求解動(dòng)力問題以探討拐折斷層的影響，強(qiáng)調(diào)了節(jié)點(diǎn)處剪切應(yīng)力降分析的重要性。焦明若[31]利用RFPA-2D分析了巖石的均質(zhì)程度對(duì)其破裂程度以及微震序列的影響，并提出不同均質(zhì)程度發(fā)育不同的地震序列類型。雖然基于拐折斷層的數(shù)值模擬研究有限，如表2所示，但也是目前很有前景的研究方法，而參數(shù)的設(shè)置應(yīng)有更多適用地球科學(xué)研究的改進(jìn)。薛霆虓[17]利用有限元法模擬了二維彈性介質(zhì)的彎曲斷層，并從前兆和余震、能量釋放、發(fā)震周期等方面與平直斷層進(jìn)行了對(duì)比，認(rèn)為彎曲斷層由于抑制了地震活動(dòng)的頻率而更易導(dǎo)致大震的發(fā)生，且彎曲斷層的前兆較平直斷層明顯。Xing等[32]通過有限元模擬提出拐折斷層大多是“雙震型”。王學(xué)濱[33]通過對(duì)FLAC-3D編程開發(fā)，統(tǒng)計(jì)分析了Z字形復(fù)雜拐折斷層系統(tǒng)中不同斷層帶上與剪切、拉伸相關(guān)的6種量的時(shí)間演變規(guī)律，同樣得出位于標(biāo)本兩側(cè)的拐折斷層和平直斷層交替活動(dòng)，而觸發(fā)下一條平直斷層的活動(dòng)并導(dǎo)致整個(gè)系統(tǒng)的失穩(wěn)。

孫楠[34]利用數(shù)值模擬的方法分析了正應(yīng)力、靜態(tài)應(yīng)力擾動(dòng)及動(dòng)態(tài)應(yīng)力擾動(dòng)等對(duì)斷層黏滑運(yùn)動(dòng)的影響。王學(xué)濱[27]從數(shù)值模型的角度分析了實(shí)驗(yàn)室尺度摩擦弱化、摩擦強(qiáng)化-摩擦弱化兩種模型（Wang等，2012,2013），前者適用于破壞和前兆的模擬，后者適合于斷層黏滑，并指出實(shí)驗(yàn)室尺度模型的模擬有助于充分了解計(jì)算模型的特點(diǎn)及優(yōu)劣。近期，Barbot等[35]國外學(xué)者使用數(shù)值模擬的方法研究了斷層黏滑運(yùn)動(dòng)及主震、同震、震間、震后等現(xiàn)象；黏滑運(yùn)動(dòng)發(fā)生的時(shí)間間隔、長期以及短期現(xiàn)象的數(shù)值模擬也取得一定的研究成果[36-37]，這些成果無疑給地震的預(yù)測與研究提供了新的思路。伴隨著對(duì)地球動(dòng)力學(xué)更深入的認(rèn)知，數(shù)值模擬技術(shù)正在發(fā)揮著越來越重要的作用[38]。因此拐折斷層黏滑運(yùn)動(dòng)過程的研究尤其是數(shù)值研究方面是未來地震機(jī)理研究的一個(gè)重要方向。

五、小結(jié)

拐折斷層地震活動(dòng)性目前的研究已經(jīng)很多，但缺乏完整而系統(tǒng)的研究。作為自然界最普遍存在的一種斷層形式，拐折斷層黏滑運(yùn)動(dòng)過程無論是實(shí)驗(yàn)室研究還是數(shù)值模擬研究都需要進(jìn)一步加強(qiáng)。

（一）方法的優(yōu)勢與局限性

實(shí)際震例的研究為地震活動(dòng)研究提供真實(shí)直接的證據(jù)，但由于地震的突發(fā)性，觀測手段的局限性，我們很難在野外直接觀測到斷層的黏滑運(yùn)動(dòng)。

實(shí)際測點(diǎn)往往反映的是斷層附近某些點(diǎn)的物理場，難以反映整個(gè)區(qū)域整體場的情況，相對(duì)而言實(shí)驗(yàn)可以通過現(xiàn)有的技術(shù)方法觀測研究模型整體不同場的變化特征，對(duì)于研究拐折斷層的活動(dòng)特征提供很多科學(xué)數(shù)據(jù)。然而，大多數(shù)實(shí)驗(yàn)室研究是有其不足之處的，例如，實(shí)驗(yàn)室研究巖石類型單一，受實(shí)驗(yàn)儀器和觀測及采樣速率的影響等。

表2 部分典型拐折斷層數(shù)值實(shí)驗(yàn)類比表

數(shù)值模擬的方法可以彌補(bǔ)以上不足，數(shù)值模擬斷層效應(yīng)從尺度和效應(yīng)上更接近于真實(shí)斷層的活動(dòng)，并且經(jīng)濟(jì)快捷，各種參數(shù)提取方便，可以方便地開展參數(shù)的敏感性研究，模擬在當(dāng)前技術(shù)條件下無法開展的物理實(shí)驗(yàn)，計(jì)算更加定量化和可視化。然而參數(shù)的選取和模型的設(shè)置存在隨意性，并不能完全符合實(shí)際地質(zhì)情況，所以更適合地學(xué)研究的物性參數(shù)和模型設(shè)置都有待進(jìn)一步改善。

（二）展望

拐折斷層系統(tǒng)破壞以及黏滑運(yùn)動(dòng)過程的綜合研究，無論從力學(xué)方面還是地質(zhì)學(xué)方面都是一個(gè)相當(dāng)復(fù)雜的過程。無論是實(shí)驗(yàn)室研究還是數(shù)值模擬研究都與科技的進(jìn)步密切相關(guān)，更高采樣速率的觀測儀器與更先進(jìn)的軟硬件計(jì)算機(jī)設(shè)備都會(huì)為研究提供更好的平臺(tái)。展望未來的研究應(yīng)是有更適用地學(xué)研究的物性參數(shù)并涉及多物理場的數(shù)值實(shí)驗(yàn)，以及有更先進(jìn)儀器設(shè)備和更合理的觀測手段的物理實(shí)驗(yàn)室研究，以便為強(qiáng)震發(fā)生的機(jī)制提供更多的理論基礎(chǔ)。

[1] 滕吉文. 強(qiáng)烈地震孕育與發(fā)生的地點(diǎn)、時(shí)間及強(qiáng)度預(yù)測的思考與探討[J]. 地球物理學(xué)報(bào), 2010, 53(8)： 1749-1766.

[2] 云龍, 郭彥雙, 馬瑾. 5°拐折斷層在黏滑過程中物理場演化與交替活動(dòng)的實(shí)驗(yàn)研究[J]. 地震地質(zhì), 2011, 33(2)： 356-368.

[3] SCHWARTZ S Y, ORANGE D L, ANDERSON R S. Complex fault interactions in a restraining bend on the San Andreas Fault, southern Santa Cruz Mountains, California[J]. Geophysical Research Letters, 1990, 17(8)： 1207-1210.

[4] AOCHI HIDEO, RAúL MADARIAGA, EIICHI FUKUYAMA. Effect of normal stress during rupture propagation rupture along nonplanar fault[J]. Journal of Geophysical Research, 2002, 107(B2)： ESE5-1-ESE5-10.

[5] KATO N, SATOH T, LEI X, et al. Effect of fault bend on the rupture propagation process of stick-slip[J]. Tectonophysics, 1999, 310(1)： 81-99.

[6] 葉建慶, 陳慧, 劉學(xué)軍, 等. 2006年云南鹽津5.1級(jí)雙震的監(jiān)測與研究[J]. 地震研究, 2008, 31(2)： 134-141.

[7] KASE Y, DAY S M. Spontenous rupture processes on a bending fault[J]. GEOPHYSICAL RESEARCH LETTERS, 2006, 33：L10302.

[8] 郭彥雙, 馬瑾, 云龍. 拐折斷層黏滑過程的實(shí)驗(yàn)研究[J]. 地震地質(zhì), 2011, 33(1)： 26-35.

[9] 馬瑾, S I SHERMAN, 郭彥雙. 地震前亞失穩(wěn)應(yīng)力狀態(tài)的識(shí)別——以5°拐折斷層變形溫度場演化的實(shí)驗(yàn)為例[J]. 中國科學(xué)：地球科學(xué), 2012, 42(5)： 633-645.

[10] 郭玲莉, 劉力強(qiáng), 馬瑾. 黏滑實(shí)驗(yàn)的震級(jí)評(píng)估和應(yīng)力降分析[J].地球物理學(xué)報(bào), 2014, 57(3)： 867-876.

[11] 馬瑾, 馬文濤, 馬勝利, 等. 5°拐折構(gòu)造變形物理場的實(shí)驗(yàn)研究與數(shù)值模擬[J]. 地震地質(zhì), 1995(4)： 318-326.

[12] ANDREWS D J. Fault geometry and earthquake mechanics[J]. Annals of Geophysics, 1994, 37(6)： 79-86.

[13] HILLERS G, BEN-ZION Y, MAI P M. Seismicity on a fault controlled by rate-and state-dependent friction with spatial variations of the critical slip distance[J]. Journal of Geophysical Research Atmospheres, 2006, 111(B1)： 997-999.

[14] CUSTóDIO S, ARCHULETA R J. Parkfield earthquakes：Characteristic or complementary?[J]. Journal of Geophysical Research Atmospheres, 2007, 112(B5)： 622-634.

[15] KING G C P, NABELEK J F. Role of fault bends in the initiation and termination of earthquake rupture[J]. Science, 1985, 228： 984-987.

[16] 馬瑾, 馬勝利, 劉力強(qiáng), 等. 斷層幾何結(jié)構(gòu)與物理場的演化及失穩(wěn)特征[J]. 地震學(xué)報(bào), 1996, 18(2)： 200-205.

[17] 薛霆虓, 傅容珊, 林峰. 幾何彎曲斷層活動(dòng)性的模擬[J]. 地球物理學(xué)報(bào), 2009, 52(10)： 2509-2518.

[18] 呂廣廷, 湯泉, 趙廣堃. 斷層泥在粘滑機(jī)制中的作用[J]. 國際地震動(dòng)態(tài), 1983(2)： 18-21.

[19] 黃建國, 張流. 水對(duì)斷層摩擦滑動(dòng)穩(wěn)定性的影響[J]. 地震地質(zhì), 2002, 24(3)： 387-399.

[20] MA S, LIU L, MA J, et al. Experimental study on nucleation process of stick-slip instability on homogeneous and nonhomogeneous faults[J]. Science China Earth Sciences, 2003, 46(2)： 56-66.

[21] 宋義敏, 馬少鵬, 楊小彬, 等. 斷層黏滑動(dòng)態(tài)變形過程的實(shí)驗(yàn)研究[J]. 地球物理學(xué)報(bào), 2012, 55(1)： 171-179.

[22] 黃元敏, 馬勝利, 繆阿麗,等. 正應(yīng)力擾動(dòng)對(duì)斷層滑動(dòng)失穩(wěn)影響的實(shí)驗(yàn)研究[J]. 地球物理學(xué)報(bào), 2016, 59(3)： 931-940.

[23] 繆阿麗. 幾種模擬斷層泥摩擦滑動(dòng)速度依賴性轉(zhuǎn)換的實(shí)驗(yàn)研究[D]. 北京： 中國地震局地質(zhì)研究所, 2011.

[24] LAPUSTA N, RICEJ R, BEN-ZION Y, et al. Elastodynamic analysis for slow tectonic loading with spontaneous rupture episodes on faults with rate- and state-dependent friction[J]. Journal of Geophysical Research Atmospheres, 2001, 105(B10)： 23765-23789.

[25] BIZZARRI A, COCCO M, ANDREWS D J, et al. Solving the dynamic rupture problem with different numerical approaches and constitutive laws[J]. Geophysical Journal International, 2010, 144(3)： 656-678.

[26] 王學(xué)濱, 馬冰, 呂家慶. 實(shí)驗(yàn)室尺度典型斷層系統(tǒng)破壞、前兆及粘滑過程數(shù)值模擬[J]. 地震地質(zhì), 2014, 36(3)： 845-861.

[27] SASANI M, YAZDANI M. Numerical simulation of earthquake mechanism based on stick-slip behavior of faults[J]. Proc Spie, 2010, 8670(4)： 126.

[28] 馬曉靜. 正斷層發(fā)生地震的動(dòng)力學(xué)過程數(shù)值模擬研究[D]. 北京：中國地震局地質(zhì)研究所, 2012.

[39] 姜輝, 鄧志輝, 高祥林, 等. 日本俯沖帶地震發(fā)生過程的數(shù)值模擬研究[J]. 地球物理學(xué)報(bào), 2013, 56(7)： 2303-2312.

[30] CRUZATIENZA V M, VIRIEUX J, OPERTO S. Dynamic Rupture Simulation of Bent Faults with a New Finite Difference Approach[C]// EGS-AGU-EUG Joint Assembly, 2003.

[31] 焦明若, 唐春安, 張國民, 等. 細(xì)觀非均勻性對(duì)巖石破裂過程和微震序列類型影響的數(shù)值試驗(yàn)研究[J]. 地球物理學(xué)報(bào), 2003, 46(5)： 659-666.

[32] XING H L, MORA P, MAKINOUCHI A. Finite Element Analysis of Fault Bend Influence onStick-Slip Instability along an Intra-Plate Fault[J]. Pure and Applied Geophysics, 2004, 161(9)： 2091-2102.

[33] 王學(xué)濱, 白雪元, 馬冰. 斷層系統(tǒng)中局部斷層活動(dòng)過程的數(shù)值分析方法——以Z字形斷層為例[J]. 地球物理學(xué)進(jìn)展, 2015(6)：2938-2944.

[34] 孫楠, 吳小平, 解朝娣, 等. 基于速率-狀態(tài)依賴摩擦本構(gòu)關(guān)系研究不同因素對(duì)斷層黏滑運(yùn)動(dòng)的影響[J]. 云南大學(xué)學(xué)報(bào)自然科學(xué)版, 2014, 36(3)： 378-383.

[35] BARBOT S, LAPUSTA N, AVOUAC J P. Under the hood of the earthquake machine： toward predictive modeling of the seismic cycle[J]. Science, 2012, 336(6082)： 707-710.

[36] SHIBAZAKI B, OBARA K, MATSUZAWA T, et al. Modeling of slow slip events along the deep subduction zone in the Kii Peninsula and Tokai regions, southwest Japan[J]. Journal of Geophysical Research Atmospheres, 2012, 117(B6)： 6311.

[37] NODA H, LAPUSTA N. Stable creeping fault segments can become destructive as a result of dynamic weakening[J]. Nature, 2013, 493(7433)： 518-21.

[38] 鄧志輝, 宋鍵, 孫君秀, 等. 數(shù)值模擬方法在地震預(yù)測研究中應(yīng)用的初步探討(Ⅰ)[J]. 地震地質(zhì), 2011, 33(3)： 670-683.

（責(zé)任編輯：劉格云）

A Review on Seismic Activity of Bending Fault

ZHANG Zeng-huan, MA Xiao-jing, WANG Shi

Hebei GEO University, Shijiazhuang, Hebei 050031

This paper analyses and summarizes the seism genic mechanism and the seismicity characteristics of bending fault, from these aspects: the bending fault seismicity characteristics, the rock physics experiments on the stick-slip motion of fault, the numerical simulation experiments on the Stick-slip motion of fault. The actual cases studies and observations show that the angle of bending fault and the bending point have important influence on the occurrence of earthquake, magnitude of earthquake and the fault rupture direction. Rock physics experiments have proved that some factors such as the material composition in fault zone, the friction coefficient of fault and the stress disturbance have a great influence on the characteristics of the fault stick-slip motion. Bending fault as a common form in nature, the study of the characteristics of stick-slip motion and its influencing factors would help us understand the mechanism of earthquakes and should be further strengthened in future research.

bending fault; seismic activity; stick-slip motion; numerical simulation

P548

A

1007-6875(2017)03-0013-06

日期：2017-03-08

10.13937/j.cnki.hbdzdxxb.2017.03.003

河北省研究生創(chuàng)新資助項(xiàng)目（1007703）。

張?jiān)鰮Q（1986—），女，河北石家莊人，碩士研究生，研究方向?yàn)閰^(qū)域地質(zhì)構(gòu)造與構(gòu)造動(dòng)力學(xué)。