喬建永，馬念杰，馬 驥，趙志強，郭曉菲，師皓宇,3

(1.中國礦業(yè)大學(北京) 能源與礦業(yè)學院，北京 100083; 2.北京郵電大學 理學院，北京 100876; 3.華北科技學院 安全工程學院，北京 101601)

地震是自然界中經常發(fā)生的自然現(xiàn)象，是地殼斷裂構造帶演化伴隨的產物，與地質斷裂構造之間在幾何特征和變化過程上具有相似性[1]。已有地震記錄顯示，許多地震區(qū)存在相互交叉的共軛發(fā)震斷層，尤其強震的震源破裂往往呈“X”型共軛破裂，地震震中位于“X”型共軛破裂交匯部位[2-4]。美國加利福尼亞州Superstition Hills地震序列在2組共軛斷層上同時出現(xiàn)地表破裂帶[5];美國蒙大拿州Red Rock Valley體波震級Mb5.3地震震源深度12.4 km，位于Monument Ridge正斷層和Red Rock斷層的共軛交匯處[6];日本Kagoshima里氏震級M6.0地震(1997)發(fā)現(xiàn)共軛斷層的存在，HORIKAWA等(2001)通過反演得出:共軛破裂起始于兩斷層段的交匯處附近并沿兩側傳播[7];YUKUTAKE等(2008)由雙差定位算法反演出主震震源與余震的位置，認為日本Niigata Chuetsu-Oki矩震級Mw6.6地震破裂主要發(fā)生在向東南傾斜的斷層面上，主震的震源周圍西北方向和東南方向傾斜的共軛斷層面幾乎同時破裂[8];CHEN等(2002)總結了臺灣Chi-Chi地震(里氏震級ML7.3)的余震活動和震源機制表明，輕微東向傾斜的車龍埔斷層和陡峭的西向傾斜的深層斷裂帶組成的共軛斷層系統(tǒng)主導了震源區(qū)南部余震的空間格局[9]。

地殼中的斷層活動誘發(fā)了地震，地震發(fā)生又促成了斷層的生成與發(fā)育，通常認為斷層的形成起始于地殼巖體的破裂。以此，學者們應用理論推斷得到的各種震源模型來解析地震發(fā)生時震源區(qū)的力學行為。這些震源模型大體可分為兩類，一類是假設作用在震源處有一組體力系的點源模型，如單力偶、雙力偶和地震矩張量震源模型。雙力偶震源模型普遍應用于定量分析斷層錯動引發(fā)地震產生地震波的特征，求解地震的兩個斷層面的走向、傾向及傾角，以及震源附近主應力方向等參數(shù)[10-12];地震矩張量震源模型可以描述各種形式的震源，如地下爆炸引起的體積膨脹源、地球內部介質體相變產生的快速擴展源或體積突然收縮源，斷層錯動等效的雙力偶震源模型只是地震矩張量描述的各類震源中的一種[13]。點源模型是通過等效力的方法將地震震源假設受到一組力系和力偶矩的作用，缺乏實際的物理意義;另一類震源模型是非點源模型，主要有用震源處某個斷層面的兩側發(fā)生位移或應變的間斷來描述震源的地震位錯模型、有限移動源(Haskell)模型、地震震源的裂紋擴展模型(Starr模型、Knopoff模型)、滑動弱化模型，描述地震震源破裂過程的障礙體(Barrier)與凹凸體(Asperity)模型，還有從摩擦強度的非均勻性入手，基于“彈性回跳假說”提出的彈簧-滑塊動力學模型(Burridge-Knopoff模型)等[10,14-21]。盡管非點源模型對地震同震所引起的位移和應變等變化及其產生后果的分析成果很多，但仍然停留在理論研究的層面，尚未形成統(tǒng)一認識和進行實際應用[22-25]。

筆者從共軛剪切破裂演化過程中地震頻發(fā)的表征現(xiàn)象出發(fā)，以動力系統(tǒng)結構穩(wěn)定性理論中的Leau-Fatou花瓣定理逆定理和蝶形破壞理論[26-34]為理論基礎，構建的“X”型共軛剪切破裂-地震復合模型，用于分析軟弱異性體(如軟弱巖體，破碎固體，融化巖漿等)周圍巖體的應力變化;推導了破裂尺寸計算公式，給出了地震能量的計算方法，并建立起地震、能量、共軛剪切破裂之間符合邏輯的數(shù)學關系;推演出共軛剪切破裂的生成及演化機理、共軛剪切破裂-地震的發(fā)生及活動規(guī)律，明確了“X”型共軛剪切破裂-地震的物理意義，為研究斷層-地震機理及其規(guī)律提供了全新的研究思路和方法，有助于豐富斷層-地震研究的理論體系。

1 蝶形花瓣逆定理及其內涵

蝶形花瓣逆定理是孔洞巖體蝶形破壞和Leau-Fatou花瓣逆定理的總稱。孔洞巖體蝶形破壞的概念來自于對圓形巷道圍巖塑性區(qū)分布規(guī)律[28-31]、圓形巷道蝶型沖擊地壓機理[32-33]、巷道蝶型煤與瓦斯突出機理[34]的研究，認為巷道圍巖出現(xiàn)劇烈變形破壞時圍巖的破壞形狀為蝶形。2000年喬建永研究平面復動力系統(tǒng)的結構穩(wěn)定性問題，得到Leau-Fatou花瓣逆定理[26-27]，指出具有花瓣幾何形狀的不穩(wěn)定點的出現(xiàn)，可以很好的描述具有突變特征的非線性復雜動力系統(tǒng)，即結構的突然失穩(wěn)破壞。

1.1 巖體的蝶形破壞

文獻[28-31]指出，處于高偏應力場中含有孔洞或軟弱異性體(圖1)的巖體將出現(xiàn)以軟弱異性體為中心的蝶形破壞區(qū)，如圖2所示，蝶形破壞區(qū)由4個蝶葉組成，分別位于最大主應力和最小主應力的兩側，兩兩對稱分布;兩個蝶葉軸線之間的夾角大體為80°～100°。圖2中，Px，Pz分別為水平主應力與垂直主應力;σ1，σ2，σ3為微小單元體i的最大、中間與最小主應力。

圖1 自然界中真實存在的軟弱異性體Fig.1 Soft anisotropic bodies in nature① —軟弱巖體;② —破碎固體;③ —融化巖漿

圖2 高偏應力場中圓形孔洞周圍的蝶形破壞Fig.2 Butterfly damage around circular holes in high deviated stress field

計算機數(shù)值模擬結果證明，無論軟弱異性體是圓形還是矩形，甚至是各種不規(guī)則的其他形狀，也無論軟弱異性體周圍巖體是均質還是非均質的彈性體，只要存在高偏應力場，都會出現(xiàn)蝶形或者殘缺蝶形的基本破壞形態(tài)，如圖3所示，數(shù)值分析模型為平面應變模型，具體尺寸為2 000 m×2 000 m(截取顯示模型500 m×500 m);軟弱異性體位于地下15 km，上覆巖層平均容重取27 kN/m3;通過采集原位地震波數(shù)據(jù)，實際觀測出露地表的巖層組合剖面，以及實驗室分析噴出火成巖的特性顯示:大陸地殼10～25 km深處以變質巖或火成巖為主，且這一區(qū)域形成由花崗質(強堿性)片麻巖組合巖體逐漸向富含角閃相的鎂鐵質礦物組合過渡[35-36]。圖3(a)，(b)反映不同形狀((a)圓形與(b)矩形)軟弱異性體圍巖塑性破壞特征，且軟弱異性體周圍巖體為均質彈性體，力學參數(shù)與應力狀態(tài)取值分別為彈性模量E=50 GPa;泊松比μ=0.25;黏聚力C=24 MPa，內摩擦角φ=46°;抗拉強度Rm=19 MPa;水平主應力Px=2 025 MPa;圖3(c)，(d)中軟弱異性體位于非均質層狀巖體中，水平主應力Px取值分別為2 050 MPa和2 025 MPa，出現(xiàn)了蝶葉穿層與殘缺蝶形破壞，巖層力學參數(shù)取值分別見表1,2。(注:參照煤的物理力學參數(shù)取值，圖3中軟弱異性體為軟弱夾層，其自身的力學參數(shù)統(tǒng)一取E=1.27 GPa;μ=0.26;C=3 MPa;φ=18°;Rm=0.9 MPa)因此，對于軟弱異性體周圍巖體破壞規(guī)律的討論，可以借用圓形孔洞周圍的蝶形破壞理論結果。

表1 層狀巖體的力學參數(shù)Table 1 Physical mechanics parameters of layered rock masses

表2 堅硬-軟弱互層的力學參數(shù)Table 2 Physical mechanics parameters of hard-weak inter-layers

文獻[28-31]以摩爾-庫侖破環(huán)準則為判定條件，推導得出含圓形軟弱異性體圍巖蝶形破壞邊界的隱性方程為

(1)

式中，a為規(guī)則圓形軟弱異性體的半徑;Px(t)，Pz(t)分別為軟弱異性體圍巖所受區(qū)域水平主應力和垂直主應力，是時間t的函數(shù);q為軟弱異性體的反作用力;r，θ為極坐標下的塑性破壞區(qū)邊界半徑和角度，θ∈[0,π]。

1.2 Leau-Fatou花瓣逆定理

花瓣圖像[26-27]是指某點周圍均勻分布若干個花瓣的數(shù)學圖形，如圖4所示?；ò攴譃槲曰ò旰统庑曰ò辏瑐€數(shù)有奇數(shù)和偶數(shù)之分。

圖4 奇數(shù)花瓣和偶數(shù)花瓣(5葉與4葉花瓣)Fig.4 Odd number and even number of petals(5 leaves and 4 leaves)

花瓣逆定理從區(qū)域的形狀判斷吸引軸和排斥軸的分布，同時指出:出現(xiàn)花瓣形狀圖像時將會發(fā)生非線性動力學現(xiàn)象。其核心是應用具有花瓣圖像的數(shù)學函數(shù)準確描述系統(tǒng)對“初始值”的“敏感依賴性”，所謂對“初始條件的敏感依賴性”是非線性動力現(xiàn)象的一個普遍特征，表達的含義是:給定初始值的函數(shù)在復平面空間中的迭代軌跡，即使很微小的改變，如計算中由四舍五入產生的小誤差，經過迭代可能變得放大了，當?shù)侥撤N狀態(tài)時微小的改變就導致整個系統(tǒng)災難性的破壞，而花瓣形狀圖像就是這一迭代軌跡表現(xiàn)出的極端不穩(wěn)定的幾何形態(tài)，也就是說確定性系統(tǒng)內也存在著隨機性的反映，在時間足夠長之后，也可能有不可預測的現(xiàn)象發(fā)生[37]。因此，可以從數(shù)學角度解釋自然界存在的穩(wěn)定狀態(tài)向非穩(wěn)定狀態(tài)的突變現(xiàn)象。

1.3 共軛剪切破裂-地震模型及其地震能計算

1.3.1共軛剪切破裂長度Rmax的計算模型

數(shù)值模擬分析的結果[28,30]表明，當?shù)纹茐膮^(qū)范圍較大時，破壞區(qū)的形狀為“X”形狀，與地質“X”型共軛剪切破裂相似，如圖5所示，模型尺寸2 000 m×2 000 m，軟弱異性體及其圍巖力學參數(shù)與圖3(a),(b)相同，水平主應力Px=2 180 MPa，垂直主應力Pz=405 MPa。所以，如果定義“X”型共軛剪切破裂由中心點向外擴展的長度為破裂長度Rmax，其中，Rmax為θ對應塑性破壞區(qū)半徑r的極大值，則可以借用式(1)來進行計算，改寫成簡式為

Rmax=af(Px(t)，Pz(t)，C，φ，θmax，q)

(2)

1.3.2共軛剪切破裂-地震能量的計算

圖5 軟弱異性體周圍破壞區(qū)的“X”形狀Fig.5 “X” shape of the failure zone around the soft anisotropic body

對比圖2和圖4可知，蝶形破壞區(qū)形狀與4葉花瓣具有相似性。所以根據(jù)Leau-Fatou花瓣逆定理可以推斷，描述共軛剪切破裂長度的式(2)將出現(xiàn)由逐漸擴展到突然擴展的突變現(xiàn)象。由此還可以推斷，基于圖2的蝶形破壞區(qū)力學模型進行能量計算也將出現(xiàn)地震時能量突然釋放的現(xiàn)象。因此，可以在蝶形塑性區(qū)理論模型的基礎上通過能量計算，獲得地震前后整個模型能量的變化情況，即計算在“X”型共軛剪切破裂擴展時釋放的地震能。

任何彈性介質在受到外力作用時都會集聚彈性能，已知蝶形破壞區(qū)力學模型，設將該空間有界閉區(qū)域Ω任意分成n個子單元體ΔVi(i=1,2,…,n)，以ΔVi表示第i個子單元體的體積。微小單元體i受到3個主應力的作用，大小分別為(圖6(a)):

σ1i≥σ2i≥σ3i

則該微小單元體的能量Ui為

2μi(σ1iσ2i+σ2iσ3i+σ3iσ1i)]ΔVi

式中，ui為微小單元體i的彈性應變能密度;Ei，μi，σ1i，σ2i，σ3i分別為第i個微小單元體的彈性模量、泊松比，以及最大、中間與最小主應力。

圖6 地震能量的計算模型Fig.6 Seismic energy calculated model

地震發(fā)生前，軟弱異性體及其周圍巖體已經產生一定的塑性破壞，空間閉區(qū)域Ω轉變?yōu)閺椥驭竐和塑性Ωp(應力極限平衡)兩種不同的區(qū)域，則子單元體分成ΔVe,i(i=1,2,…,m)，ΔVp,j(j=m+1,m+2,…,n)(其中Ωe∪Ωp=Ω)，i,j為彈性區(qū)與塑性破壞區(qū)的微小單元體變量，如圖6(b)所示。(假設整個空間區(qū)域范圍內，塑性破壞區(qū)Ωp的微小單元體仍處于彈性狀態(tài)下)則整個空間閉區(qū)域的彈性能為

(3)

地震發(fā)生時，軟弱異性體周圍巖體塑性Ωp區(qū)域出現(xiàn)擴展，擴展時間為Δt，擴展范圍為ΔΩp，整個空間閉區(qū)域的彈性能為

(4)

地震發(fā)生使得整個空間閉區(qū)域Ω的能量比前一時刻減小，這部分消耗的能量一部分用于巖石內結晶晶格錯位，部分生成熱量，還有一部分會引起巖體震動并以地震波的形式傳播出去，即地震能。地震能其表達式為

Wm=β(Uep-U′ep)

(5)

式中，Wm為地震能，即地震發(fā)生時Δt秒釋放的以軟弱異性體圍巖形成的蝶形破壞為震中的震源能(本文假設Δt=1 s)，106J;β為震動能因子，0<β<1，本文中β取值為0.5;Wm為一個以Px,y,z(t)為自變量的復合函數(shù)。

2 “X”型共軛剪切破裂形成機理及演化

2.1 “X”型共軛剪切破裂的形成機理

已建立的“X”型共軛剪切破裂-地震復合模型(圖7)，及其推導得出的“X”型共軛剪切破裂長度的數(shù)學表達式(2)，可用于很好地定量描述共軛剪切破裂形成的力學機理:在高偏應力場作用下，軟弱異性體周圍會形成蝴蝶形狀的高剪切應力區(qū)，該區(qū)域內的巖體與應力狀態(tài)滿足摩爾-庫侖剪切破壞條件而遭到破壞，因而形成以軟弱異性體為中心的“X”型共軛剪切破裂，“X”型共軛剪切破裂長度可通過式(1)，(2)求解。

由于數(shù)學力學分析方法的限制(無法完全按照自然界的實際條件進行解析)，為便于理論計算，假設軟弱異性體為規(guī)則的圓形，位于地下深處h=15 km，半徑a=10 m，縱向長度L=10 km，軟弱異性體周圍2.5 km半徑范圍內的巖體為花崗巖(圖7)。容重取γ=27 kN/m3，軟弱異性體周圍巖體泊松比μ=0.25，彈性模量E=50 GPa，黏聚力C=24 MPa，內摩擦角φ=46°;破壞后塑性巖體泊松比μp=0.35;假設軟弱異性體圍巖處于均勻應力場中，即Px=Pz=γh=405 MPa(等壓狀態(tài));為了分析的方便，設構造應力場變化中垂直主應力Pz(t)=405 MPa，并不隨時間的變化而改變。則“X”型共軛剪切破裂-地震復合模型中只存在水平主應力Px(t)一個與時間有關的自變量，其余參數(shù)都已給出確定值。于是，令Pz(t)=405 MPa，Px(t)=2 583 MPa，a=10 m，C=24 MPa，φ=46°，q=202.5 MPa，可計算得到“X”型共軛剪切破裂的具體長度Rmax=r=303.55 m，如圖7所示。

圖7 “X”型共軛剪切破裂尺寸的計算結果(地下深處15 km)Fig.7 Calculational results of “X” conjugate shear fracture’s radius (It’s 15 kilometers underground)

2.2 “X”型共軛剪切破裂的演化特性

斷裂活動和地震的發(fā)生都與集中于地殼內部構造應力場的作用密切相關[38]。以水平主應力Px(t)為自變量，通過公式(1)計算得到單位時間每增加1 MPa，“X”型共軛剪切破裂長度擴展變化量ΔR，如圖8所示。以Px=Pz=405 MPa的等壓狀態(tài)為中心，曲線右側:隨水平主應力增大，ΔR呈正指數(shù)型增長;曲線左側:隨水平主應力減小，ΔR呈負指數(shù)型增長;定義曲線右側隨水平主應力增大形成水平高偏應力場為擠壓應力區(qū)域，相應的隨水平主應力減小形成的垂直高偏應力場為張拉應力區(qū)域;以擠壓應力區(qū)為例進一步分析，當Px(t)≥2 583 MPa后(該拐點處曲線的斜率k=1)，對應的破裂長度增量ΔR變化曲線近乎與x軸垂直，此時已形成“X”型共軛剪切破裂(圖5)，即數(shù)學上Leau-Fatou花瓣逆定理描述的4葉花瓣圖形，則地殼軟弱異性體與周圍巖體組成的非線性動力系統(tǒng)對于“初始狀態(tài)”具有敏感依賴性。也就是說，在一定地應力和圍巖環(huán)境中，破裂的擴展對處于“初始狀態(tài)”下的條件變化具有敏感依賴性，地應力和地層強度的微小變化將會引發(fā)“X”型共軛剪切破裂的不斷擴展，有時甚至是劇烈擴展。

圖8 水平應力增量為1 MPa時，破裂尺寸增量ΔR的變化關系Fig.8 Curves of the increment ΔR of fractures as the horizontal stress increment is 1 MPa

2.3 “X”型、“Y”型與“V”型共軛剪切破裂特征的成因

“X”型共軛節(jié)理(破裂)是一種普遍存在的節(jié)理系統(tǒng)，規(guī)模大，長度一般由數(shù)千米到十千米，最長30 km，表3所示理論計算得到的“X”型共軛剪切破裂尺寸Rmax與自然界中真實存在的“X”型共軛剪切破裂長度相吻合[39-40]。地質學家通過露頭地層剖面記錄到復雜的共軛破裂體系整體呈“X”型、“Y”型和“V”型的共軛特征[41]，可由式(1)的計算結果很好地定量解釋。如圖7所示，當A,B和C區(qū)域都為軟弱巖層時，軟弱異性體周圍巖體會呈現(xiàn)出“X”型共軛破壞特征;如果B,C為強度高的堅硬巖層，軟弱異性體仍然位于軟弱巖層A中，則共軛剪切破裂的發(fā)展會受到堅硬巖層的阻斷，從而形成“V”型共軛破裂。位于猶他州漢克斯維爾的Entrada和Curtis地層中共軛破裂的現(xiàn)場觀測充分證實了這一點:“V”型共軛破裂交叉點集中在富含淤泥的軟弱地層中并沿著軟弱巖層發(fā)育明顯[42-43]。數(shù)值模擬反演結果也普遍證實了:破裂交叉點附近的幾何特征和軟弱巖性夾層可能對共軛破裂的生長和發(fā)展產生強烈影響[43];當共軛剪切破裂某個發(fā)展方向上巖體強度比軟弱異性體周圍巖體的強度高(圖7中C區(qū)域為堅硬巖層，A,B區(qū)域為軟弱巖層)，則該方向的破裂發(fā)展就會受到限制，從而呈現(xiàn)出“Y”型的共軛特征;如果軟弱異性體周圍只有一個方向上為軟弱巖層(圖7中B為軟弱巖層，其余區(qū)域為堅硬巖層)，則會形成單個的破裂特征。綜上分析可知，共軛剪切破裂的擴展對于地層強度的微小變化具有敏感依賴性，它既可以引發(fā)“X”型共軛破裂的不斷擴展，也可以使得花瓣缺失或殘缺，出現(xiàn)不完整蝶形的形態(tài)，即軟弱異性體周圍巖體強度特征變化，使得共軛破裂表現(xiàn)出“X”,“V”,“Y”型共軛特征，圖9數(shù)值模擬結果充分驗證了理論分析的真實性。其中，數(shù)值模型、應力條件、規(guī)則圓形軟弱異性體尺寸及圍巖力學參數(shù)(選用圖3(d)的參數(shù))如圖3所示。

表3 不同水平主應力對應“X”型共軛破裂尺寸的計算結果Table 3 Calculational results of different horizontal principal stresses corresponding to “X” conjugate fracture’s radius

注:假設Pz=405 MPa;q=202.5 MPa。

圖9 不同類型共軛剪切破裂形成的數(shù)值模擬結果Fig.9 Numerical simulation results of shear fractures formation of different types

3 地震與“X”型共軛剪切破裂一次性擴展尺寸的關系

在“X”型共軛剪切破裂-地震復合模型及共軛剪切破裂尺寸數(shù)學表達(式(2))的基礎上，從數(shù)學上推演得到的地震釋放能量計算式(5)，可以明確“X”型共軛剪切破裂一次性擴展長度與地震釋放能量大小的關系。如圖10所示，當水平應力增量ΔPx=1 MPa/s時，應用式(2),(5)分別計算得到共軛剪切破裂尺寸一次性擴展量ΔR、地震釋放能量Wm和相應里氏震級ML。分析表明:隨水平主應力Px的變化，地震釋放能量Wm與地震震級ML具有同共軛剪切破裂一次性擴展量ΔR相同的指數(shù)型增長趨勢;同時可以看出，ΔR和Wm都是水平主應力Px的單值函數(shù)，二者之間是一一對應關系。當圖中的水平應力不斷增加或減小，“X”型共軛剪切破裂也在不斷擴展，“X”型共軛剪切破裂尺寸與地震釋放的能量相對于水平主應力Px的變化具有同步性，地震與“X”型共軛剪切破裂擴展之間為相互伴生關系，它們同時出現(xiàn)，也同時消失。

圖10 破裂尺寸一次性擴展量、釋放能量與地震震級Fig.10 One-time expansion of the size of fractures,release energy and seismic magnitude

表4為不同破裂尺寸條件下水平應力擾動為0.1 MPa時，“X”型共軛剪切破裂尺寸增量ΔR、地震能Wm與對應里氏震級計算結果，圍巖力學參數(shù)與圖3數(shù)值模擬選用參數(shù)保持一致?？梢钥闯觯癤”型共軛剪切破裂擴展一次釋放的能量級別可以涵蓋自然發(fā)生的各種地震震級和每次地震震中釋放的能量;相同的應力擾動在不同的應力狀態(tài)下引發(fā)“X”型共軛剪切破裂一次性擴展的尺度不同，從而引發(fā)不同級別的地震(例如水平應力Px=2 591 MPa時，增加0.1 MPa，破裂一次性擴展2.4 m，釋放的能量量級可達1013，相當于里氏6.0級地震)。對于一般地質力學條件下“X”型共軛剪切破裂引發(fā)地震的能量計算，式(5)同樣適用，而且會得出與本文給出的完全相同的結論。

表4 “X”型共軛剪切破裂尺寸增量ΔR、地震能Wm與對應里氏震級計算結果Table 4 Calculational results of increment ΔR of “X” conjugate fracture’s radius,seismic energy Wm and corresponding Richter Magnitude

4 自然地震對微小應力的敏感依賴

“X”型共軛剪切破裂-地震復合模型用定量方法解釋了許多微小應力變化觸發(fā)自然界中較大地震的現(xiàn)象。根據(jù)式(5)可以求得不同應力環(huán)境下，水平應力擾動分別為1，10-1，10-2，10-3，10-4，10-5，10-6，10-7和10-8MPa/s時，引發(fā)地震所釋放的能量與對應的震級(里氏震級通過公式lgE=11.8+1.5Ms換算得到[44],Ms為面波震級，見表5)。需要指出:地應力本身是沒有時間概念的，水平應力擾動是有作用時間的，文中水平應力擾動是以秒為單位，反映出每秒釋放的能量。由表5可以看出，對于相同應力狀態(tài)的“X”型共軛剪切破裂來說，擾動應力ΔPx值越大則發(fā)生地震的震級越高。相同于“X”型共軛剪切破裂擴展對應力狀態(tài)和區(qū)域圍巖強度的敏感依賴性，地震對于“X”型共軛剪切破裂所處的應力狀態(tài)和區(qū)域圍巖強度(震前狀態(tài))也同樣具有敏感依賴性，微小、甚至極其微小的應力擾動即可觸發(fā)大級別的自然地震，很好地解釋了自然界中普遍存在的僅僅0.01 MPa甚至更微小的應力變化觸發(fā)地震的現(xiàn)象[45-47]。例如，蘇門答臘-安達曼地震(M9.0，2004)(the Sumatra-Andaman earthquake)震后僅0.01 MPa的局部應力擾動就引發(fā)附近M8.6地震[48];汶川M8.0地震發(fā)生后，與其震中相隔85 km的廬山地震(Ms=7.0)前后應力變化0.02～0.04 MPa(0.2～0.4 bar)[49];甚至距離3 100 km處Mw=7.9Denali地震震中(美國阿拉斯加)產生的0.16～0.22 MPa動態(tài)觸發(fā)應力，亦可引發(fā)有感地震[50];日本北部Iwate火山巖漿膨脹使得區(qū)域應力場發(fā)生變化，0.33 MPa的應力增加引發(fā)了Ms=6.1級地震[51];此外，平均量級103Pa的固體潮汐應力、地球自轉以及人類活動等產生的應力擾動也可對地震的發(fā)生產生重要影響[52-55]。

表5 不同觸發(fā)應力Px，地震能Wm與對應里氏震級計算結果(圍巖力學參數(shù)選取同表3,4)Table 5 Calculational results of different trigger stresses Px,seismic energy Wm and corresponding Richter Magnitude (The surrounding rock mechanics parameters are shown in Table 3,4)

表5 不同觸發(fā)應力Px，地震能Wm與對應里氏震級計算結果(圍巖力學參數(shù)選取同表3,4)Table 5 Calculational results of different trigger stresses Px,seismic energy Wm and corresponding Richter Magnitude (The surrounding rock mechanics parameters are shown in Table 3,4)

觸發(fā)應力ΔPx/(MPa·s-1)地震能/JPx=1 100 MPaPx=2 420 MPaPx=2 598 MPa里氏震級Px=1 100 MPaPx=2 420 MPaPx=2 598 MPa17.87×1093.56×10123.32×10163.45.27.810-17.86×1083.27×10111.62×10152.74.56.910-27.86×1072.08×10101.54×10142.13.76.310-37.86×1062.08×1091.56×10131.43.05.610-47.85×1052.08×1081.19×10120.72.34.810-57.5×1042.08×1071.47×10110.11.74.210-65.00×1032.07×1061.97×108—1.02.310-7—2.05×1051.96×107—0.31.710-8—2.00×1041.95×106——1.0

5 結 論

(1)應用“X”型共軛剪切破裂-地震復合模型的理論計算結果推演了共軛剪切破裂形成的力學機理:在高偏應力場作用下，軟弱異性體周圍巖體會形成蝴蝶形狀的高剪切應力區(qū)，該區(qū)域內的巖體與應力狀態(tài)滿足摩爾—庫侖剪切破壞條件而遭到破壞，因而形成以軟弱異性體為中心的“X”型共軛剪切破裂。

(2)軟弱異性體周圍巖體強度特征變化，使得共軛破裂表現(xiàn)出“X”,“V”,“Y”型共軛特征;共軛剪切破裂尺寸一次性擴展與地震釋放能量具有同步性，即地震與“X”型共軛剪切破裂是相互伴生關系，它們同時出現(xiàn)，并同時消失。

(3)在實際地層中存在某種地應力環(huán)境，使得地殼軟弱異性體與其周圍巖體組成的非線性動力系統(tǒng)對于“初始狀態(tài)”具有敏感依賴性，“初始狀態(tài)”可以是軟弱異性體自身的形狀、尺寸，也可以是所受應力環(huán)境和周圍巖體物理力學性質。也就是說，某些情況下，地應力的微小變化會引發(fā)“X”型共軛剪切破裂的較大擴展和較大地震。