童亨茂，陳正樂， 劉瑞珣

①中國石油大學(北京)油氣資源與探測國家重點實驗室，北京 102249；②中國地質(zhì)科學院地質(zhì)力學研究所，北京 100086；③北京大學地球與空間科學學院，北京 100871

廣義剪切活動準則*

童亨茂①?，陳正樂②， 劉瑞珣③

①中國石油大學(北京)油氣資源與探測國家重點實驗室，北京 102249；②中國地質(zhì)科學院地質(zhì)力學研究所，北京 100086；③北京大學地球與空間科學學院，北京 100871

針對經(jīng)典的剪切破裂準則(Coulomb-Mohr準則和以Byerlee律為代表的滑動摩擦律)存在局限，在Coulomb-Mohr準則和Byerlee律的基礎(chǔ)上，從產(chǎn)生剪切破裂的物理本質(zhì)出發(fā)，應用先存構(gòu)造活動性準則和活動趨勢分析理論，提出 “廣義剪切活動準則”。該準則可以用來定量判斷任意介質(zhì)、在任意三軸應力狀態(tài)下、任意方位界面(包括先存薄弱面和非薄弱面)發(fā)生剪切活動的可能性，并把Coulomb-Mohr準則從均勻介質(zhì)擴展到有先存構(gòu)造的非均勻介質(zhì)，從臨界應力狀態(tài)擴展到任意應力狀態(tài)，把Byerlee律為代表的滑動摩擦律從二維應力狀態(tài)擴展到任意三軸應力狀態(tài)，具有廣闊的應用前景。

剪切活動；先存構(gòu)造；Coulomb-Mohr準則；Byerlee律；廣義剪切活動準則

斷層是巖石圈的最基本的構(gòu)造要素之一[1]。地球上幾乎所有的內(nèi)動力地質(zhì)作用和大部分外動力地質(zhì)作用都與斷層有關(guān)，同時礦產(chǎn)資源(如石油和天然氣)和地質(zhì)災害(如地震、火山、山體滑坡等)都直接或間接地與斷層相關(guān)[1-2]。

斷層作用的物理本質(zhì)是巖石的剪切活動。巖石發(fā)生剪切破裂的經(jīng)典準則是Coulomb-Mohr準則[1]。然而，Coulomb-Mohr準則只適用于均勻介質(zhì)，產(chǎn)生的斷層是方位優(yōu)選的[3]，而且只能預測破裂發(fā)生前的瞬間趨勢，不能用來分析破裂發(fā)生后的遞進變形過程。為了克服均勻介質(zhì)的限制，不同學者先后提出了沿先存構(gòu)造面滑動的摩擦滑動定律[4-6]和滑動趨勢分析理論[7]，但這些理論依然不能很好地解決自然界非均勻介質(zhì)斷層作用的力學機理問題。

本文在Coulomb-Mohr準則的基礎(chǔ)上，從產(chǎn)生剪切破裂的物理本質(zhì)(作用在介質(zhì)潛在剪破裂面上的剪應力克服該界面上的內(nèi)聚力和(內(nèi))摩擦力，即可產(chǎn)生剪切破裂)出發(fā)，應用先存構(gòu)造活動性準則(activation criterion of pre-existing fabric)[8]和活動趨勢分析理論(theory of activation tendency analysis)[9]，提出了適用于任意應力狀態(tài)、任意介質(zhì)的脆性剪切活動準則，即“廣義剪切活動準則”，并進一步討論與經(jīng)典剪切破裂準則的關(guān)系及應用前景。

1 Coulomb-Mohr準則

巖石發(fā)生剪切破裂的基本準則是Coulomb-Mohr準則[1,6]。

Coulomb剪切破裂準則是Coulomb于1776年提出的[10]。該準則認為：巖石發(fā)生剪切破壞不僅與作用在截面上的剪應力有關(guān)，而且還與作用在該截面上的正應力有關(guān)，可以用下式表示[1-2]：

其中t為臨界剪應力，σn為該截面上的正應力，μ為介質(zhì)的內(nèi)摩擦系數(shù)(常數(shù))。式(1)在t-σn坐標系中為一直線(圖1)，稱為庫侖剪破裂線， μ為該直線的斜率。若μ用庫侖剪破裂線的傾斜角φ表示，則μ= tanφ，φ為巖石的內(nèi)摩擦角。發(fā)生Coulomb剪破裂的臨界應力狀態(tài)被稱為Coulomb應力狀態(tài)。

圖1 Coulomb-Mohr準則的摩爾圓圖解[1-2]

Mohr剪切破裂準則認為，內(nèi)摩擦系數(shù)μ及內(nèi)摩擦角φ并不是一常數(shù)，它與σn相關(guān)，一般隨著σn的增大而減小。這樣，摩爾剪切破裂線是曲線，被稱之為摩爾破裂包絡(luò)線(圖1)[1-2]。

Coulomb準則針對均勻介質(zhì)[3]。Anderson根據(jù)Coulomb剪切破裂準則，指出均勻介質(zhì)的巖石在應力作用下發(fā)生的剪切破裂是方位優(yōu)選的，即剪切破裂面與最大主壓應力軸σ1的夾角(剪裂角)等于±(45°-φ/2)(其中φ為摩擦角)，并與中間主應力軸σ2平行，是一對共軛的破裂面[3]。在此基礎(chǔ)上，他又進一步提出了Anderson斷層模式，成為構(gòu)造地質(zhì)學的世紀經(jīng)典理論，至今還產(chǎn)生重要影響[10]。但該斷層模式基于均勻介質(zhì)，實際應用中具有顯著的局限性[8,11]。

2 滑動摩擦律(Byerlee律)

基于Coulomb-Mohr準則的 Anderson模式被廣泛深入應用后，發(fā)現(xiàn)有很多實際地質(zhì)現(xiàn)象無法用Anderson模式來解釋，主要存在三個方面的問題[8,11]。其中之一，它不符合沿先存構(gòu)造面上滑動的斷層作用。

Anderson斷層模式假定地質(zhì)體是均勻介質(zhì)，預測的斷層是方位優(yōu)選的，只與主應力的方位有關(guān)，與其他(如先存構(gòu)造面)無關(guān)。然而，實際地質(zhì)體是非均質(zhì)的。

地質(zhì)體的非均質(zhì)性存在多種形式，如沉積盆地中沉積巖的層理、基底中先存斷裂構(gòu)造、面理等。砂箱實驗研究表明，在沉積盆地中，先存構(gòu)造(主要是先存斷裂)對非均質(zhì)性的影響一般最為顯著，是造成地質(zhì)體非均質(zhì)的主導因素[8,12]。一般情況下(如露頭尺度)，地殼巖石的非均質(zhì)性是很強的，這是造成Anderson模式存在局限，無法很好地解釋自然界斷層作用現(xiàn)象的根本原因。

為了克服上述問題，先后提出了滑動摩擦律[4-6,13-14]和滑動趨勢分析理論[7]。摩擦律認為，先存斷裂在地應力作用下，當作用在先存斷裂面上的剪應力克服其內(nèi)聚力和摩擦力后，先存斷裂可以重新活動[4-5,14]，其中最著名的是在實驗基礎(chǔ)上提出的Byerlee 律[5]，可以用下式表示：

其中tW為臨界剪應力，CW為先存斷裂面的內(nèi)聚力，σn為先存斷裂面上的正應力，μW為先存斷面上的摩擦系數(shù)。

滑動摩擦律提出沿先存構(gòu)造面滑動的斷層作用克服了Anderson斷層模式中均質(zhì)體的假設(shè)和斜向滑動的問題[4-5,13]，但這些模型都局限于二維應力分析[7]?；瑒于厔莘治隼碚揫7]把沿薄弱面的摩擦滑動分析從二維擴展到三維，在三軸應力已知的情況下(主應力值確定，方向直立或水平，即Anderson應力狀態(tài))，定量分析了先存薄弱面發(fā)生活動的可能性，并用來分析不同先存構(gòu)造地震活動的可能性。但該模型只考慮Anderson應力狀態(tài)(一個主應力軸直立)，而忽略了薄弱面的內(nèi)聚力，也沒有考慮巖石內(nèi)摩擦系數(shù)的變化[9]，因此也具有一定的局限性。

滑動摩擦律在一定程度上克服了Anderson模式的局限性，但也沒有從根本上解決自然界復雜介質(zhì)條件下斷層作用的力學機制問題。

3 剪切活動的物理本質(zhì)與廣義剪切活動準則

3. 1 剪切活動的物理本質(zhì)

分析Coulomb-Mohr準則和Byerlee律可以看出，不論是均勻介質(zhì)還是非均勻介質(zhì)，產(chǎn)生剪切活動的物理本質(zhì)是共同的，即作用在界面上的剪應力只要能克服該界面上的內(nèi)聚力和(內(nèi))摩擦力，即可產(chǎn)生剪切活動。本文使用“剪切活動”而非“剪切破裂”一詞，是因為沿失去內(nèi)聚力(已經(jīng)破壞)的構(gòu)造面發(fā)生剪切活動不需要再發(fā)生剪切破壞，即“剪切活動”存在兩種情況：①產(chǎn)生剪切活動的潛在界面是有內(nèi)聚力的，則剪切活動通過剪切破壞而產(chǎn)生，需要克服內(nèi)聚力和內(nèi)摩擦力；②產(chǎn)生剪切活動的界面無內(nèi)聚力(如先存斷裂面)，則剪切活動只需克服界面的摩擦力即可產(chǎn)生。

這樣，從剪切活動的物理本質(zhì)出發(fā)，對于任意介質(zhì)中的任意界面，發(fā)生剪切活動的臨界剪應力[tn]可以用下式表示：

式中， [tn]為沿界面的臨界剪應力，CP、μP和σn分別為該界面的內(nèi)聚力、(內(nèi))摩擦系數(shù)和正應力。當界面沒有內(nèi)聚力時，CP取0。

3.2 剪切活動趨勢因子及其數(shù)學表達

通過上述剪切活動物理本質(zhì)的分析，對任意介質(zhì)中的任意界面，定義一參數(shù) fs稱之為剪切活動趨勢因子，用來描述任意界面發(fā)生剪切活動的可能性。fs可表示為

式中，tn為沿界面的剪應力，[tn]為沿該界面產(chǎn)生剪切活動時的臨界剪應力。很顯然，當 fs＝1時，該界面處于臨界剪切活動狀態(tài)；當 fs＞1時，處于已活動狀態(tài)；當 fs＜1時，處于穩(wěn)定狀態(tài)。

結(jié)合式(3)和(4)，得到

在任意三軸應力狀態(tài)下，童亨茂等[8]給出了σn和tn的定量計算公式[8]：

結(jié)合式(5)～(8)，在任意三軸應力狀態(tài)下，得到剪切活動趨勢因子 fs的計算公式：

其中σ1、σ2和σ3分別為最大、中間和最小主壓應力，θ為界面與σ1的夾角，α為界面在σ2-σ3平面上的交線與σ3的夾角。界面的地理方位(如傾向和傾角)與θ和α之間的關(guān)系可以通過坐標轉(zhuǎn)換來建立[9]。

很容易看出，這里的剪切活動趨勢因子fs和童亨茂等[8]提出的先存構(gòu)造活動性準則中的先存構(gòu)造活動性系數(shù) fa具有相同形式的表達式，只是 fa針對的是先存構(gòu)造面，而 fs針對任何界面。用式(9)可以定量計算在任意三軸應力狀態(tài)下、任意方位界面(包括先存薄弱面和非薄弱面)發(fā)生剪切活動的可能性。

3.3 剪切活動趨勢因子在摩爾空間中的表達

據(jù)式(9)，剪切活動趨勢因子 fs與應力狀態(tài)(3個主應力值)、界面相對主應力的方位(θ和α)及其界面的力學性質(zhì)(CP，μP)有關(guān)，但由于 fs的計算公式(式9)比較復雜，不方便對各影響因素進行定量分析。然而，利用Tong和Yin[9]提出的“摩爾空間”圖解，則可以很直觀地解決這一問題，也就是可以很直觀地分析剪切活動趨勢因子fs與應力狀態(tài)、界面方位和界面力學性質(zhì)等之間的定量關(guān)系。

在任一應力狀態(tài)下，應用摩爾空間可以十分直觀地表達任一確定方位界面的正應力、剪應力以及界面的力學性質(zhì)(圖2)[15]：①按照界面的方位(θ和α)找到該界面在摩爾空間中的極點，該極點的橫坐標和縱坐標分別對應界面上正應力(σn)和剪應力(tn)；②界面的臨界剪切活動線(式3的圖形表達)的起點(CP, 0)決定于界面的內(nèi)聚力 CP,斜率決定于(內(nèi))摩擦系數(shù)μP。這樣，當應力狀態(tài)(決定摩爾空間)、界面的力學性質(zhì)(決定界面的臨界活動線)和界面的方位(決定界面在摩爾空間中極點的位置)確定時，應用摩爾空間就可以十分方便、直觀地表達和預測任一界面的剪切活動性(圖2)。剪切活動趨勢因子 fs由其極點P(σn, tn)與其相對應的臨界活動點(σn,)的相對位置決定，等于tn和的比值(圖2)，即當界面的極點位于：①臨界剪切活動線上時(如圖2中的P3)，界面處于臨界活動狀態(tài)；②臨界剪切活動線之外時(圖2中有橫線的黃色區(qū)域，如P1)，界面已經(jīng)活動；③臨界剪切活動線之內(nèi)(圖2中沒有橫線的黃色區(qū)域，如P2)，界面處于穩(wěn)定狀態(tài)。這種判別方法適用于任何界面(包括薄弱面和非薄弱面)：當界面為非薄弱界面時，臨界剪切活動線為Coulomb-Mohr包絡(luò)線；當界面為先存薄弱面時，臨界剪切活動線為先存構(gòu)造活動線(如圖2中的a和a′)[8]。

圖2 “摩爾空間”及剪切活動趨勢因子fs的表達(據(jù)文獻[15]修改)(σ1-σ3、σ1-σ2和 σ2-σ3三個摩爾圓以及黃色區(qū)域構(gòu)成“摩爾空間”；綠色虛線和粉色虛線為界面方位角θ和φ的等值線；a和a′為先存構(gòu)造活動線；P1、P2和P3分別為3個先存構(gòu)造界面(其力學性質(zhì)CW、μW一致)在摩爾空間中的極點(σ, t)、(σ, t)和(σ, t)，tW、n1n1n2n2n3n3n1tn2W和分別為對應的臨界剪應力，對應的剪切活動趨勢因子fs分別為fS1＞1.0、fS2＜1.0和 fS3=1.0；CW為先存構(gòu)造界面的內(nèi)聚力；A和A′為Coulomb剪切破裂點)

3.4 廣義剪切活動準則及其與Coulomb-Mohr準則、Byerlee律之間的關(guān)系

3.4.1 剪切活動趨勢因子 fs算式在實際地質(zhì)體中的適用性

從前面的分析可以看出，剪切活動趨勢因子fs的算式(式9)理論上適用于任何介質(zhì)，但式(9)隱含了兩個條件：計算區(qū)域內(nèi)的應力場是均勻的；界面是一平面，而且其力學性質(zhì)是一樣的。只有包含平面(板)狀先存構(gòu)造的非均勻介質(zhì)(先存構(gòu)造為平面(板)狀，且不同位置的力學性質(zhì)一致；母質(zhì)為均勻介質(zhì))，才能滿足上述兩個條件。那么，式(9)是否適用于實際的非均勻地質(zhì)體呢？

研究表明，先存構(gòu)造對地質(zhì)體非均質(zhì)性的影響最為顯著[8,12]。Morley等[16]把先存構(gòu)造分為分隔性構(gòu)造和透入性構(gòu)造：前者是指介質(zhì)中存在一個截然的破裂面，相當于先存斷層，內(nèi)聚力很小或為零；后者是由一系列彌漫性構(gòu)造組成，力學性質(zhì)上為一相對薄弱帶，內(nèi)聚力小于周圍介質(zhì)。

按照有限單元的思想，任何類型的實際非均勻地質(zhì)體，只要進行有限單元的劃分，都可以看作是含平面(板)狀先存構(gòu)造的介質(zhì)和均勻介質(zhì)的組合體。例如：砂泥巖互層、層內(nèi)均勻的沉積巖，可以看作是有先存薄弱帶(泥巖)的均勻砂巖；復雜形態(tài)的先存斷裂面，只要把斷面分割成足夠小的有限單元，單元內(nèi)的斷面可以看作是平面；被不規(guī)則巖脈侵入的巖體，不管侵入接觸的邊界多么復雜，只要把邊界附近的巖體分割成足夠小的有限單元，單元內(nèi)就可以看作是有分隔性先存構(gòu)造(侵入接觸面)的均勻巖體。

從上面的分析可以看出，任何實際的非均勻地質(zhì)體，從力學上均可以看作是含平面(板)狀先存構(gòu)造(薄弱帶)的均勻介質(zhì)(如含一條平面狀斷層的均勻地質(zhì)體)及其組合體。因此從理論上說，式(9)適用于任何實際的地質(zhì)體。

3.4.2 式(9)與Coulomb-Mohr準則的關(guān)系

很容易驗證，當Cp=C，μP=μ，即在均勻介質(zhì)的情況下，若fs＝1(發(fā)生臨界剪切活動)，式(9)可以簡化成式(1)，即Coulomb-Mohr準則。這樣，Coulomb-Mohr準則是式(9)在均勻介質(zhì)和臨界剪切活動情況下的端元。

Coulomb-Mohr準則是式(9)端元的情況還可以用摩爾空間來直觀地說明。圖2中的摩爾空間代表臨界Coulomb應力狀態(tài)，大摩爾圓與Coulomb剪破裂線相切，切點為A和A′。在此摩爾空間中，極點位于先存構(gòu)造活動線(a和a′)上及其以外區(qū)域(圖2中有橫線的黃色區(qū)域)的先存構(gòu)造界面(CW、μW均一致)都達到了剪切活動的條件(fs≥1)，即在臨界Coulomb應力狀態(tài)下，存在無數(shù)個潛在的先存薄弱面滿足剪切活動條件。隨著CW的變大，先存構(gòu)造活動線(圖2中的a和a′)向外移動，滿足剪切活動條件的區(qū)間(圖2中有橫線的黃色區(qū)域)逐漸變??；當CW=C，且μP=μ，即有先存構(gòu)造的介質(zhì)轉(zhuǎn)化為均勻介質(zhì)后，該區(qū)間退縮成兩個點(如圖2中的A和A′，代表一對共軛剪切破裂)，轉(zhuǎn)化為Coulomb-Mohr準則(fs=1)，同時先存構(gòu)造活動線(如圖中的a和a′)轉(zhuǎn)化為Coulomb-Mohr剪破裂包絡(luò)線。這樣，式(9)把Coulomb-Mohr準則從均勻介質(zhì)擴展到有先存構(gòu)造的非均勻介質(zhì)，從臨界應力狀態(tài)擴展到任意應力狀態(tài)。同時，式(9)也揭示了按Coulomb-Mohr準則產(chǎn)生的剪破裂為方位優(yōu)選(optimal)[3]的物理本質(zhì)：在均勻介質(zhì)所有方位的界面中，只有方位優(yōu)選的界面才能達到臨界剪切活動的條件fs=1。

3.4.3 式(9)與滑動摩擦律的關(guān)系

很容易驗證，滑動摩擦律的代表Byerlee律是式(9)在二維應力狀態(tài)下(界面與中間主應力σ2軸平行)，先存薄弱面處于臨界剪切活動狀態(tài)(fs=1)的端元。在圖2的摩爾空間中，Byerlee律代表先存構(gòu)造活動線a和a′和大摩爾圓相交的4個點。這表明式(9)把以Byerlee律為代表的滑動摩擦律擴展到任意三軸應力狀態(tài)。

前面的分析表明，剪切活動趨勢因子 fs的算式(9)一方面把經(jīng)典的剪切破裂準則Coulomb-Mohr準則和Byerlee律統(tǒng)一起來，并把它們作為其中的端元；另一方面，理論上式(9)適用于任何非均勻介質(zhì)和任意應力狀態(tài)。為此，稱剪切活動趨勢因子 fs的算式(9)為“廣義剪切活動準則”，可以定量判斷任意介質(zhì)、在任意三軸應力狀態(tài)下、任意方位界面(包括先存薄弱面和非薄弱面)發(fā)生剪切活動的可能性。

4 廣義剪切活動準則的應用前景分析

4.1 斷層作用研究中的應用

斷層作用的物理本質(zhì)是地質(zhì)體沿界面發(fā)生剪切活動。Anderson于1905年基于Coloumb準則，建立了巖石圈在Anderson應力狀態(tài)(其中一個主應力直立)和均勻介質(zhì)的假定條件下的斷層作用模式，即Anderson斷層模式[3]。但實際地質(zhì)體存在不同程度的非均質(zhì)性，造成Anderson模式存在局限，無法很好地解釋自然界斷層作用現(xiàn)象[8]。由于廣義剪切破裂準則適用于任何非均勻介質(zhì)和任意應力狀態(tài)，這樣，基于廣義剪切活動準則建立的“廣義斷層模式”(將在另文中介紹)適用于巖石圈中的任意應力狀態(tài)和任何類型的地質(zhì)體，可以合理地解決自然界的各種斷層作用現(xiàn)象，使得斷層作用的理論模型更加符合地質(zhì)實際，具有更好的理論指導作用。

4.2 其他方面的應用

理論上，廣義剪切活動準則可以應用到與剪切活動相關(guān)的各種作用的分析研究中，如邊坡穩(wěn)定性分析和預測，斷層的地震活動性分析，油氣開發(fā)過程中的注水、注氣和熱采等能否導致斷層的重新活動(斷層的重新活動導致鉆井和套管的破壞)的研究分析，隧道(或坑道)開挖時斷裂的穩(wěn)定性分析，以及與大型工程地基相關(guān)的斷層穩(wěn)定性分析等。具體的分析方法不在這里一一詳述。

5 討論

廣義剪切活動準則隱含了兩個條件：①計算區(qū)域內(nèi)的應力場是均勻的；②界面是一平面，而且其力學性質(zhì)是一致的。在要求十分精確的情況下，只有把非均勻介質(zhì)分割成足夠小的微元后，才能達到廣義剪切活動準則的應用條件。

由于地質(zhì)體的非均質(zhì)性和地質(zhì)作用的復雜性，到目前為止，還未有任何一個理論能精確分析和預測地質(zhì)作用過程，但這并不影響相關(guān)地質(zhì)理論的應用和價值。適當?shù)恼`差是允許的，對廣義剪切活動準則也是如此。任何理論應用均需對實際地質(zhì)體和地質(zhì)作用過程作適當和必要的簡化，但只要新的理論所對應的模型相對更加符合地質(zhì)實際，得到的認識更為深化，預測相對更為準確，就能體現(xiàn)該理論的意義和價值。例如：①自然界實際的斷層面都是曲面，但并非需要分割成無限小后才能看作是平面，而是綜合考慮相關(guān)數(shù)據(jù)的誤差情況和工作需要后，把斷層面分割成若干段落，每個段落的斷面即可認為是平面，其內(nèi)的應力場可以認為是均勻的，即可達到廣義剪切活動準則的應用條件。②在巖石圈應力狀態(tài)下，由于重力的作用，主應力值隨深度是不斷變化的，應力場必然是非均勻的；但如果只分析斷層的方位，在介質(zhì)是均勻的情況下，應力場也可以被看作是均勻的[3]。如童亨茂等[9,18]在砂箱實驗時，發(fā)現(xiàn)砂箱內(nèi)按均勻應力場預測得到的結(jié)果與模擬實驗結(jié)果具有驚人的一致性。

當然，在條件具備的情況下，單元分割得越精細，用廣義剪切活動準則預測的結(jié)果就相對越準確。如一條斷層面上的應力分析和力學性質(zhì)分析達到很精細的情況下，就可以應用廣義剪切活動準則實現(xiàn)對斷面上不同區(qū)域的剪切活動趨勢作出更為精細的預測結(jié)果。

值得注意的是，剪切活動一旦產(chǎn)生，剪切面附近就會產(chǎn)生應力降，形成局部應力場[17-20]。另外，在剪切活動過程中(斷層作用過程)，隨著地質(zhì)體的變形和變位，局部應力場也會不斷發(fā)生變化[18]。因此，應力和應變是緊密關(guān)聯(lián)的兩個方面，即應力導致應變，應變又引起應力的變化，構(gòu)成應力-應變體系(stress-strain system)。在遞進變形過程中，應力場的演化過程十分復雜，其研究目前還剛起步。廣義剪切破裂準則是聯(lián)系應力和剪切應變(斷層作用)的橋梁，而褶皺變形往往與斷層作用緊密相關(guān)。隨著觀測資料的不斷精細(如高精度三維地震資料)，變形結(jié)果(構(gòu)造幾何特征)的可觀測性和變形過程(構(gòu)造的運動學過程)的可解析性不斷增強，廣義剪切破裂準則理論有望在遞進變形過程中應力場演化這一研究領(lǐng)域發(fā)揮重要作用。

6 結(jié)論

廣義剪切活動準則，從理論上可以定量判斷任意介質(zhì)、在任意三軸應力狀態(tài)下、任意方位界面(包括先存薄弱面和非薄弱面)發(fā)生剪切活動的可能性，并可對不同區(qū)域進行對比分析。廣義剪切活動準則把Coulomb-Mohr準則從均勻介質(zhì)擴展到有先存構(gòu)造的非均勻介質(zhì)，從臨界應力狀態(tài)擴展到任意應力狀態(tài)，以及把Byerlee律為代表的滑動摩擦律從二維應力狀態(tài)擴展到任意三軸應力狀態(tài)。

在實際應用中，綜合考慮相關(guān)數(shù)據(jù)的誤差情況和工作需要，把要研究的界面(如斷層面)分割成若干段落后，認為滿足①計算區(qū)域內(nèi)的應力場是均勻的；②界面是一平面，而且其力學性質(zhì)是統(tǒng)一的這兩個條件后，即可應用廣義剪切活動準則，來定量分析和預測界面不同段落發(fā)生剪切活動的可能性。

(2015年9月23日收稿)

[1] TWISS R J, MOORES E M. Structural geology [M]. New York: W. H. Freeman and Company, 1992: 532.

[2] 徐開禮, 朱志澄. 構(gòu)造地質(zhì)學[M]. 北京: 地質(zhì)出版社, 1998: 270.

[3] ANDERSON E M. The dynamics of faulting [M]. 2nd ed. Edingburgh: Oliver and Boyd, 1951.

[4] BOTT M H P. The mechanics of oblique slip faulting [J]. Geol Mag, 1959, 96: 109-117.

[5] BYERLEE J D. Friction of rocks [J]. Pure and Application Geophysics, 1978, 116: 615-626.

[6] JAEGER J C, COOK N G W. Fundamentals of rock mechanics [M]. 3rd ed. London: Chapman and Hall, 1979: 593.

[7] MORRIS A, DAVID A F, HENDERSON B. Slip-tendency analysis and fault reactivation [J]. Geology, 1996, 24: 275-278.

[8] TONG H, CAI D, WU Y, et al. Activity criterion of pre-existing fabrics in non-homogeneous deformation domain [J]. Science China (Earth Sciences), 2010, 53(8): 1115-1125.

[9] TONG H, YIN A. Reactivation tendency analysis: A theory for predicting the temporal evolution of preexisting weakness under uniform stress state [J]. Tectonophysics, 2011, 503: 195-200.

[10] CéLéRIER B. Seeking Anderson’s faulting in seismicity: a centennial celebration [J]. Reviews of Geophysics, 2008, 46: 1-34.

[11] 童亨茂. “不協(xié)調(diào)伸展”作用下裂陷盆地斷層的形成和演化模式[J]. 地質(zhì)通報, 2011, 29(11): 1606-1613.

[12] 童亨茂, 孟令箭, 蔡東升, 等. 裂陷盆地斷層的形成和演化——目標砂箱模擬實驗與認識[J]. 地質(zhì)學報, 2009, 83(6): 759-774.

[13] WALLACE R E. Geometry of shearing stress and relation to faulting [J]. J Geol, 1951, 59: 118-130.

[14] MCKENZIE D P. The relation between fault plane solutions for earthquakes and the directions of the principal stresses [J]. Bull Seismol Soc Am, 1969, 59: 591-601.

[15] 童亨茂, 王建君, 趙海濤, 等. “摩爾空間”及其在先存構(gòu)造活動性預測中的應用[J]. 中國科學: 地球科學, 2014, 44(9): 1948-1957.

[16] MORLEY C K, HARANYA C, PHOOSONGSEE W S. Activation of rift oblique and rift parallel pre-existing fabrics during extension and their effect on deformation style: examples from the rifts of Thailand [J]. J Structural Geology, 2004, 26: 1803-1829.

[17] GUDMUNDSSON A, SIMMENES T H, LARSEN B, et al. Effects of internal structure and local stresses on fracture propagation, deflection, and arrest in fault zones [J]. J Structural Geology, 2010, 32: 1643-1655.

[18] TONG H, KOYI H, HUANG S, et al. The effect of multiple preexisting weaknesses on formation and evolution of faults in extendedsandbox models [J]. Tectonophysics, 2014, 626: 197-212.

[19] MORLEY C K, SEUSUTTHIY G K. Fault superimposition and linkage resulting from stress changes during rifting: Examples from 3D seismic data, Phitsanulok Basin, Thailand [J]. J Structural Geology, 2007, 29: 646-663.

[20] YIN A. Mechanics of monoclonal systems in the Colorado plateau during the Laramide orogen [J]. J Geophys Res, 1994, 99: 22043-22058.

Generalized shear activation criterion

TONG Hengmao①, CHEN Zhengle②, LIU Ruixun③
①Key National Laboratories of Oil and Gas Resources and Exploration, China University of Petroleum, Beijing 102249, China;②Institute of Geomechanics, Chinese Academy of Geological Sciences, Beijing 100086, China; ③School of Earth and Space Sciences, Peking University, Beijing 100871, China

There are some limitations in classic shear fracture criteria such as Coulomb-Mohr criterion and sliding friction law represented by Byerlee law. On the base of Coulomb-Mohr criterion and Byerlee law, starting from the physical nature of shear fracture, and applying the activation criterion of pre-existing weakness and the theory of activation tendency analysis, a generalized shear activation criterion is proposed. The criterion can be used to quantitatively determine the probability of shear activation for any medium, any stress state, any orientation plane. It extends not only Coulomb-Mohr criterion from homogeneous media to inhomogeneous media with pre-existing weakness, but the sliding friction law represented by Byerlee law from 2D stress state to any triaxial stress state, and shows broad application prospects.

shear activation, pre-existing weakness, Coulomb-Mohr criterion, Byerlee law, generalized shear activation criterion

(編輯：沈美芳)

10.3969/j.issn.0253-9608.2015.06.007

*國家自然科學基金 (41272160、40772086) 和國家油氣重大專項(2011ZX05006-006-02-01、2011ZX5023-004-012)資助

?通信作者，E-mail：tong-hm@163.com