留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

广义破裂活动准则

童亨茂 张宏祥 侯泉林 陈正乐 侯贵廷

童亨茂,张宏祥,侯泉林,等,2024. 广义破裂活动准则[J]. 地质力学学报,30(1):3−14 doi: 10.12090/j.issn.1006-6616.2023180
引用本文: 童亨茂,张宏祥,侯泉林,等,2024. 广义破裂活动准则[J]. 地质力学学报,30(1):3−14 doi: 10.12090/j.issn.1006-6616.2023180
TONG H M,ZHANG H X,HOU Q L,et al.,2024. Generalized fracturing activation criteria[J]. Journal of Geomechanics,30(1):3−14 doi: 10.12090/j.issn.1006-6616.2023180
Citation: TONG H M,ZHANG H X,HOU Q L,et al.,2024. Generalized fracturing activation criteria[J]. Journal of Geomechanics,30(1):3−14 doi: 10.12090/j.issn.1006-6616.2023180

广义破裂活动准则

doi: 10.12090/j.issn.1006-6616.2023180
基金项目: 国家自然科学基金项目(41272160,20772086);国家油气重大专项(2011zx05-006-02-01,2011zx5023-004-012)
详细信息
    作者简介:

    童亨茂(1967—),男,博士,教授,主要从事盆地构造分析,构造物理模拟实验等方面的研究。Email:tonghm@cup.edu.cn

  • 中图分类号: P553

Generalized fracturing activation criteria

Funds: This research is financially supported by the National Natural Science Foundation of China (Grants No. 41272160 and 20772086) and the National Oil and Gas Major Projects (Grants No. 2011zx05-006-02-01 and 2011zx5023-004-012).
  • 摘要: 岩石破裂及其后续活动是最基本的构造变形方式,然而,经典的破裂准则(Coulomb-Mohr准则、Griffith准则和以Byerlee律为代表的滑动摩擦律)在实际应用中存在不同方面的局限。文章在经典破裂准则的基础上,从产生破裂(扩张破裂和剪切破裂)的物理本质出发,结合广义剪切活动准则和长期的研究实践,通过理论分析,提出了“广义破裂活动准则”。该准则用统一的无量纲因子(破裂活动趋势因子)来表达扩张破裂活动和剪切破裂活动,可以用来定量判断任意介质、在任意三轴应力状态下、任意方位界面(包括先存薄弱面和非薄弱面)发生破裂活动的可能性以及破裂的类型,把Coulomb-Mohr准则、Byerlee律和Griffith准则统一起来,将破裂(fracturing)扩展到破裂活动(fracturing activity)。该准则在与破裂活动相关的资源勘探开发(如页岩气、干热岩等)、灾害(如地震与滑坡等)预测和防治中具有广阔的应用前景。

     

  • 破裂(fracture,包括断层和节理)是岩石圈最基本的构造要素之一(Twiss and Moores,1992),是脆性变形的基本方式。几乎所有的内动力地质作用和大部分外动力地质作用都与之相关,同时资源和能源的形成(特别是资源,如石油和天然气)和地质灾害的发生(如地震、滑坡等)都直接或间接地与之相关(徐开礼和朱志澄,1989Twiss and Moores,1992)。

    从物理本质上,岩石破裂形成存在两种基本方式:剪切破裂和扩张破裂(Twiss and Moores,1992)。

    岩石发生剪切破裂的经典准则是Coulomb-Mohr准则(Twiss and Moores,1992)。然而,Coulomb-Mohr准则只适用于均匀介质,产生的断层是方位优选的(Anderson,1951),而且只能预测破裂发生前的瞬间趋势(不能用来分析破裂发生后的递进变形过程;Tong et al.,2010),实际应用中存在显著局限。在岩石力学实验基础上提出的、以Byerlee律(Byerlee,1978)为代表的滑动摩擦律克服了均质体的假设和斜向滑动的问题(Wallace,1951Bott,1959Byerlee,1978),但这些模型都局限于二维应力状态。

    广义剪切活动准则(Tong et al.,2010Tong and Yin,2011童亨茂等,2015)从产生剪切破裂的物理本质(作用在介质潜在剪破裂面上的剪应力克服该界面上的内聚力和(内)摩擦力,即可产生剪切破裂)出发,应用“先存构造活动性准则”(Activation criterion of pre-existing fabric;Tong et al.,2010)和“活动趋势分析理论”(Theory of activation tendency analysis;Tong and Yin,2011),提出了适用于任意应力状态、任意介质的脆性剪切活动准则,从理论上建立了适用于实际地质体的剪切活动准则。但该准则只适用于判别剪切活动。

    扩张破裂的经典准则是Griffith准则(Griffith,1921)。该准则考虑了微裂纹的应力集中效应,解决了宏观物体的破裂强度大大低于理论值的问题,但没有涉及具有宏观先存构造(包括具有内聚力的薄弱面和没有内聚力的先存破裂构造)的扩张破裂问题。如有天然裂缝发育的页岩,压裂过程中产生的破裂就不能简单地用Griffith准则来判别,因为压裂过程中,不仅产生垂直最小主应力方向的张裂缝,而且先存的天然裂缝也可能产生扩张活动(童亨茂等,2021a)。

    经典的破裂准则(包括Coulomb-Mohr准则和Griffith准则)均没有考虑宏观先存构造,而实际地质体普遍发育宏观先存构造(薄弱面),如断层、节理、层理、软弱层、面理、不整合面、侵入接触边界等。这些先存构造有的有内聚力(如层理、面理等),有的没有内聚力或内聚力很小(如断层和节理)。在有先存构造的情况下,无论是扩张破裂还是剪切破裂均可能不是方位优先的。针对实际的地质体,经典破裂准则存在显著的局限性。

    另外还存在一个没有得到很好解决的问题:当最小有效主应力为零或负值(张应力)的应力状态时,地质体如何产生脆性破坏,即在这种应力状态下,剪切活动能否产生?是否只能产生扩张破裂?这些问题在理论上一直没有得到很好的解决。

    针对上述问题,文章从脆性破裂(剪切破裂和扩张破裂)的物理本质出发,结合广义剪切活动准则,提出了在任意应力状态下,用破裂趋势因子(fracturing tendency factor,fF,无量纲)定量评价有宏观先存构造地质体中任意界面(包括薄弱面)发生破裂及进一步活动(合称破裂活动)的可能性及破裂活动的类型(扩张破裂活动或剪切破裂活动),即“广义破裂活动准则”。该准则不仅克服了经典的扩张破裂准则和剪切破裂准则存在的局限,而且把它们统一在了一起,并把破裂扩展到破裂活动,理论上适用于任何实际的地质体,应用前景广阔。

    在岩石力学实验中,岩石样品在围压不太高的情况下一旦发生宏观破裂就会失去承载力,实验一般到此结束,因此很少再去关心岩石发生破裂以后的事情,包括已有的破裂准则。然而,地质体则不同,发生破裂后,承载力虽然会有一定不同程度的减小(如地震后发震断层附近产生应力降(Bailey and Ben-Zion,2009),但不会完全失去。破裂产生以后,只要载荷(如构造作用力、工程注入流体压力等)继续存在,破裂作用就会继续。其中,剪切破裂发生后,还会沿破裂面产生剪切滑移(断层作用);张破裂发生后,破裂面还会产生扩张(扩张作用);同时,破裂活动过程中,破裂面均会不同程度地产生扩展(破裂面变大)。

    地质体产生破裂及其后续的活动,文中统称为破裂活动(fracturing activity),并划分为剪切破裂活动(简称剪切活动)和扩张破裂活动(简称扩张活动)。

    在有先存薄弱面(带)的情况下(尤其是含有没有黏聚力的先存破裂),破裂活动与均质体存在很大的区别。在载荷作用下,没有黏聚力的先存破裂的破裂活动直接进入破裂的后续活动阶段(扩张或滑移),破裂面的扩展方式也比较独特(Tong et al.,2014)。经典的Coulomb-Mohr准则和Griffith准则均没有考虑宏观先存构造,不适用于有先存薄弱面地质体,也无法解释和预测其破裂作用过程。

    需要特别指出的是,简单地用“破裂”一词难以准确地描述有先存薄弱面地质体的破裂作用过程,为此,文中采用“破裂活动”一词。

    Coulomb剪破裂准则是Coulomb于1776年提出的(Célérier,2008),可以用公式(1)表示(徐开礼和朱志澄,1989Twiss and Moores,1992):

    $$ \tau =C+\mu {\sigma }_{n} $$ (1)

    其中,τ为临界剪应力,σn为该截面上的正应力,C为介质的内聚力(常数),μ为介质的内摩擦系数(常数)。公式(1)在τ-σ坐标系中为一条直线(图1),称为“库伦剪破裂线”,μ为该直线的斜率。其中μ = tanφφ为岩石的内摩擦角,T为岩石的抗张强度,内聚力C=2T。发生Coulomb剪破裂的临界应力状态被称为Coulomb应力状态。文中的应力符号采用固体力学的定义(压应力为正,张应力为负)。

    图  1  Coulomb-Mohr准则与Griffith准则图解
    φ为岩石的内摩擦角;T为岩石的抗张强度;内聚力C=2T
    Figure  1.  Diagram of Coulomb-Mohr criterion and Griffith criterion
    φ is the internal friction angle of the rock, T is the tensile strength of the rock and the cohesion C=2T

    Mohr在岩石力学实验的基础上,发现岩石内摩擦系数μ(或内摩擦角φ)并不是一个常数,它与σn相关,一般随着σn的增大而减小。这样,Mohr对Coulomb准则进行了修正(把剪破裂线从Coulomb准则的直线变为曲线即“Mohr破裂包络线”;图1),即Mohr准则。这两个经典剪破裂准则合称为Coulomb-Mohr准则。

    Coulomb准则是针对均匀介质提出的(Anderson,1951)。Anderson根据Coulomb剪破裂准则,指出均匀介质的岩石在应力作用下发生的剪破裂是方位优选的,即剪破裂面与最大主压应力轴σ1的夹角(剪裂角)等于±(45°−φ/2)(φ为内摩擦角),并与中间主应力轴σ2平行,是一对共轭的破裂面(Anderson,1951),并进一步提出了Anderson断层模式,成为构造地质学的世纪经典理论(Célérier,2008)。但该断层模式是基于均匀介质的,在实际应用中具有显著的局限性(Tong et al.,2010)。

    基于Coulomb-Mohr准则的Anderson断层模式假定地质体是均匀介质,预测断层是方位优选的,只与主应力的方位有关,与其他因素(如先存构造面)无关。

    然而,实际地质体是非均质的。先存构造(主要是先存断裂)对非均质性的影响一般最为显著,是造成地质体非均质的主导因素(童亨茂等,2009Tong et al.,2010),这是造成Anderson模式存在局限、无法很好地解释自然界断层作用现象的根本原因。

    为了解决在先存构造面上断层复活的问题,不同学者先后提出了摩擦滑动律(Wallace,1951Bott,1959McKenzie,1969Byerlee,1978Jaeger and Cook,1979)和滑动趋势分析理论(Morris et al.,1996)。

    摩擦律认为,先存破裂在地应力作用下,当作用在先存破裂面上的剪应力克服其内聚力和摩擦力后,先存破裂可以重新活动(Bott,1959McKenzie,1969Byerlee,1978),其中最著名的是在实验基础上提出的Byerlee律(Byerlee,1978),可以用公式(2)表示。

    $$ {\tau }_{\mathrm{W}}={C}_{\mathrm{W}}+{\mu }_{\mathrm{W}}{\sigma }_{n} $$ (2)

    式中,τW为临界剪应力,CW为先存破裂面的内聚力,σn为先存破裂面上的正应力,μW为先存破裂面上的摩擦系数。

    摩擦滑动律提出沿先存构造面滑动的断层作用克服了Anderson断层模式中均质体的假设和斜向滑动的问题(Wallace,1951Bott,1959Byerlee,1978),但这些模型都局限于二维应力分析(Morris et al.,1996)。滑动趋势分析理论把沿薄弱面的摩擦滑动分析从二维扩展到三维,但该模型也具有一定的局限性:一是只考虑Anderson应力状态(一个主应力轴直立);二是忽略了薄弱面的内聚力;三是没有考虑岩石内摩擦系数的变化(Tong and Yin,2011)。

    滑动摩擦律在一定程度上克服了Anderson模式的局限性,但都没有从根本上解决自然界复杂介质条件下断层作用的力学机制问题。

    分析Coulomb-Mohr准则和Byerlee律可以看出,不论是均匀介质还是非均匀介质,产生剪切活动的物理本质是共同的,即作用在界面上的剪应力只要能克服该界面上的内聚力和(内)摩擦力,即可产生剪切活动。

    文中使用“剪切活动”而非“剪切破裂”一词,是因为沿失去内聚力(已经破坏)的先存构造面发生剪切活动不需要再发生剪切破坏,即“剪切活动”存在两种情况:一是产生剪切活动的潜在界面是有内聚力的,则剪切活动通过剪切破坏而产生,需要克服内聚力和(内)摩擦力;二是产生剪切活动的界面无内聚力,则剪切活动只需克服界面的摩擦力即可产生。

    童亨茂等(2015)提出的广义剪切活动准则从剪切活动的物理本质出发,应用先存构造面活动准则和活动趋势分析理论,提出用剪切活动趋势因子(fS;公式(3))来表达任意三轴应力状态下,任意界面(方位定义见图2)发生剪切活动的可能性(具体公式推导过程见Tong et al.,2010):

    图  2  界面的空间方位定义(以3个主应力为坐标轴)
    σ1σ2σ3分别为最大、中间和最小主应力;θ为界面与σ1的夹角;α为界面在σ2σ3平面上的交线与σ3的夹角
    Figure  2.  The spatial orientation definition of the interface (with 3 principal stresses as axes)
    σ1, σ2 and σ3 are the maximum, intermediate and minimum principal stresses respectively, θ is the angle between the interface and σ1, and α is the angle between σ3 and the interface's intersection line with the σ2σ3 plane.
    $$ {f}_{\mathrm{S}}=\frac{\sqrt{({{\sigma }_{1}}^{2}{\mathrm{sin}}^{2}\theta +{{\sigma }_{2}}^{2}{\mathrm{cos}}^{2}\theta {\mathrm{cos}}^{2}\alpha +{{\sigma }_{3}}^{2}{\mathrm{cos}}^{2}\theta {\mathrm{sin}}^{2}\alpha )-({\sigma }_{1}{\mathrm{sin}}^{2}\theta +{\sigma }_{2}{\mathrm{cos}}^{2}\theta {\mathrm{cos}}^{2}\alpha +{\sigma }_{3}{\mathrm{cos}}^{2}\theta {\mathrm{sin}}^{2}\alpha {)}^{2}}}{{C}_{\mathrm{P}}+{\mu }_{\mathrm{P}}\times ({\sigma }_{1}{\mathrm{sin}}^{2}\theta +{\sigma }_{2}{\mathrm{cos}}^{2}\theta {\mathrm{cos}}^{2}\alpha +{\sigma }_{3}{\mathrm{cos}}^{2}\theta {\mathrm{sin}}^{2}\alpha )} $$ (3)

    式中,σ1σ2σ3分别为最大、中间和最小主应力;CP为界面的内聚力;μP为界面上的摩擦系数;θ为界面与σ1的夹角;α为界面在σ2σ3平面上的交线与σ3的夹角。需要说明的是,界面的地理方位(如倾向和倾角)与θα之间的关系可以通过坐标转换来建立(Tong and Yin,2011)。

    很容易验证,当CP = CμP = μ,即在均匀介质的情况下,若fS= 1.0(即发生临界剪切活动),公式(3)可以简化成公式(1),即Coulomb-Mohr准则。这样,Coulomb-Mohr准则是公式(3)在均匀介质和临界剪切活动情况下的端元。同样很容易验证,滑动摩擦律的代表Byerlee律是公式(3)在二维应力状态下(界面与中间主应力σ2轴平行),界面处于临界剪切活动状态(fS = 1.0)的端元(童亨茂等,2015)。这样,公式(3)就把Coulomb-Mohr准则和以Byerlee律为代表的滑动摩擦律统一在了一起。

    此外,应用有限单元的思想,只要进行有限单元的划分,任何类型的实际非均匀地质体,都可以看作是平面(板)状先存构造(薄弱带)和均匀介质的组合体(如含一条平面状断层的均匀地质体;童亨茂等,2015)。因此,从理论上,公式(3)适用于任何实际的地质体。

    剪切活动趋势因子fS(公式(3))一方面把经典的剪切破裂准则Coulomb-Mohr准则和Byerlee律统一起来,并把它们作为其中的端元;另一方面,理论上适用于任何非均匀介质和任意应力状态。为此,将公式(3)称为“广义剪切活动准则”,可以定量判断任意介质、在任意三轴应力状态下、任意方位界面(包括先存薄弱面和非薄弱面)发生剪切活动的可能性(童亨茂等,2015)。广义剪切活动准则用无量纲的剪切活动趋势因子定量描述地质体中不同先存构造面发生剪切活动的可能性及活动能力大小,从理论上解决了实际地质体断层作用的力学机制问题,为实际地质体断层作用模型的建立提供了理论基础。

    经典的张破裂准则是Griffith准则。Griffith(1921)指出,产生破裂与微裂纹的存在有关,当σ3 = −Tσ1<−3σ3时,微裂纹产生扩展并最终形成与σ1平行的张破裂;当(σ1σ3)2/(σ1+σ3) ≤8T时(假定材料的抗压强度是抗张强度的8倍),则与σ1大约呈30°角的微裂纹最易扩展并产生破裂,破裂包络线方程为:

    $$ {\tau }_{n}^{2}=4T\left({T+\sigma }_{n}\right) $$ (4)

    式中,τn为临界剪应力,T为岩石的抗张强度,σn为该截面上的正应力。

    Griffith准则考虑了微裂纹的应力集中效应,解决了宏观物体的破裂强度大大低于理论值的问题;同时该准则实质上把张破裂和剪切破裂联合起来了(图1显示公式(4)在应力摩尔圆图解上是一条抛物线型的破裂包络线)。但Griffith准则没有涉及具有宏观先存构造(包括没有内聚力的破裂面和具有内聚力的先存薄弱面)的扩张破裂问题,也没有阐述抛物线型的破裂包络线不同段落所代表的物理本质。

    剪切活动和扩张活动的力学机制存在差异,下面分别加以分析。

    不论是均匀介质还是非均匀介质,产生剪切活动的物理本质是共同的,即只要作用在任意界面(先存薄弱面或假想界面)上的剪应力能克服其内聚力和(内)摩擦力,该界面即可产生剪切活动(童亨茂等,2015);同样,产生扩张活动的物理本质也是共同的,即只要作用在任意界面(先存薄弱面或假想界面)上的有效正应力(σnP0)能克服其扩张强度TP(指界面两侧介质间的黏聚力,对于非薄弱面,TP = TP0为界面上的流体压力),该界面即可产生扩张活动,即当σn−P0 = −TP时,界面即可产生扩张活动。在地壳应力状态下,3个主应力一般都是挤压应力,出现张应力的情况可以忽略(Zoback,2007),地质体内产生扩张活动一般必须有流体压力参与才能实现。因此,一般情况下,地质体内界面产生扩张破裂的物理实质是:流体压力克服了界面上的正应力和扩张强度。

    广义剪切活动准则用无量纲的剪切活动趋势因子(fS;公式(3))定量描述地质体中不同构造面(先存薄弱面或假想界面)发生剪切活动的可能性及活动能力大小。按照类似的思想,文中定义扩张活动趋势因子(fT;公式(5))来定量描述地质体中不同界面(先存薄弱面或假想界面,界面的方位也按图2定义)发生扩张活动的可能性及活动能力大小:

    $$ {f}_{\mathrm{T}}=\left({P}_{0}-{\sigma }_{n}\right)∕{T}_{\mathrm{P}} $$ (5)

    按照界面方位在图2中的定义,可以推导出任意三轴应力状态下、任意方位界面fT的表达式(公式6):

    $$ {f}_{\mathrm{T}}=\left({P}_{0}-({\sigma }_{1}{\mathrm{sin}}^{2}\theta +{\sigma }_{2}{\mathrm{cos}}^{2}\theta {\mathrm{cos}}^{2}\alpha +{\sigma }_{3}{\mathrm{cos}}^{2}\theta {\mathrm{sin}}^{2}\alpha )\right)∕{T}_{\mathrm{P}} $$ (6)

    在给定的应力状态下,当fT = 1.0,界面处于临界扩张活动状态;fT> 1.0,界面已发生扩张活动;fT< 1.0,界面处于稳定状态。

    利用摩尔空间(图3)很容易验证,当fT = 1.0、TP = T(即均匀介质)时,只有与σ3垂直的界面(θ= 0°、α =90°)才能满足扩张破裂的条件(其它方位的界面均无法满足),把fT = 1.0、TP = Tθ= 0°、α =90°代入到公式(6),即可得到:σ3′ = −T(其中σ3′ = σ3P0,表示最小有效主应力),这表明:①经典的Griffith张破裂准则是公式(6)为均匀介质(TP = T)处于临界扩张活动条件(fT = 1.0)时的端元;②均匀介质产生的扩张破裂是方位优选的,即扩张破裂面垂直于σ3。文中把σ3′ = −T的应力状态称为Griffith应力状态,与横坐标轴垂直且截距为−T的直线称为Griffith张破裂线(图3中直线g),同时把与横坐标轴垂直且截距为−TP的直线称为薄弱面扩张活动线(图3中直线gP)。

    图  3  Griffith应力状态下的摩尔空间图解
    T为岩石的抗张强度;2T为岩石的内聚力;TP为薄弱面的扩张强度;g为Griffith张破裂线;gP为薄弱面的扩张活动线;θ 为界面与 σ 1的夹角; α 为界面在 σ 2σ 3平面上的交线与 σ 3的夹角
    Figure  3.  Diagram of the Mohr space under Griffith stress state
    T is the tensile strength of the rock, 2T is the cohesion of the rock, TP is the extension strength of the weakness plane, g is the Griffith tensile fracture line, gP is the extension activation line of the weakness plane, θ is the angle between the interface and σ 1, and α is the angle between σ 3 and the interface's intersection line with the σ 2- σ 3 plane.

    而在有薄弱面存在的情况下(TP<T),当达到Griffith应力状态时,位于图3粉色区域(直线gP左侧(含)的摩尔空间区域)不同方位的薄弱面均已达到了扩张破裂的条件(fT ≥ 1.0)。这表明在有薄弱面存在的情况下,产生的扩张破裂不是方位优选的。

    考虑到公式(6)适用于任意有薄弱面的非均匀介质(包括任何实际的地质体)以及任意方位的界面,且涵盖了经典的张破裂准则,将其称之为广义扩张活动准则。

    由于剪切活动趋势因子(fS)和扩张活动趋势因子(fT)均是无量纲的,数值可直接对比。在给定的应力状态下,通过对比任意给定界面fSfT值的大小,即可判定其是否发生破裂活动及破裂活动的类型。当fS>fT ≥ 1.0或fS ≥ 1.0 >fT时,界面发生剪切活动;当fT>fS ≥ 1.0或fT ≥ 1.0>fS时,界面发生扩张活动;当fS = fT ≥ 1.0时,界面发生扩张活动和剪切活动的可能性一致;当fT< 1.0且fS< 1.0时,界面保持稳定。为此,文中把公式(3)和公式(6)合称为广义破裂活动准则,并用公式(7)表达:

    $$ {f}_{\mathrm{F}}=\mathrm{M}\mathrm{A}\mathrm{X}({f}_{\mathrm{S}},{f}_{\mathrm{T}}) $$ (7)

    式中,fF为破裂活动趋势因子。

    考虑到除先存薄弱面外,介质发生破裂均是方位优选的。因此,广义破裂活动准则中的“界面”实际针对的是各种薄弱面和方位优选的界面。

    需要特别注意的是,由于Coulomb应力状态和Griffith应力状态均是介质承受应力的极限,介质的应力状态不可能超越Coulomb应力状态和Griffith应力状态,即其在摩尔空间中的极点(投影点)不可能越过Coulomb剪破裂线(图4中的直线c和c′)及Griffith张破裂线(图4中的直线g)。这是个简单而特别重要的认识。

    图  4  给定应力状态下摩尔空间中不同极点与破裂活动线的关系图解
    T为岩石的抗张强度;2T为岩石的内聚力;C 为介质的内聚力;C P为界面的内聚力TP为薄弱面的扩张强度;g为Griffith张破裂线;gP为薄弱面的扩张活动线;P1—P9为不同薄弱面在摩尔空间中的极点位置
    Figure  4.  Diagram of relationship between different poles and fracturing activity line in Mohr space under given stress state
    T is the tensile strength of the rock, 2T is the cohesion of the rock, C is the cohesion of the medium, CP is the cohesion of the interface, TP is the extension strength of the weakness plane, g is the Griffith tensile fracture line, gP is the extension activation line of the weakness plane, and P1–P9 are the pole positions of different weakness planes in Mohr-space.

    利用摩尔空间(图4)可以非常直观地表达广义破裂活动准则。① 摩尔空间中界面的极点位于纵坐标轴的右侧时(如P1、P2、P3),不会产生扩张活动,只可能产生剪切活动;当极点在其剪切活动线(图4中的直线cP和cP′)外侧时(如P1),该界面发生剪切活动;当极点在其剪切活动线上时(如P2),该界面处于临界剪切活动状态;当极点在其剪切活动线内侧时(如P3),该界面处于稳定状态。② 摩尔空间中界面的极点位于的纵坐标的左侧时(如P4—P9),极点在其剪切活动线内侧、同时在扩张活动线(图4中的直线gP)右侧时(如P4),既不发生扩张活动,也不发生剪切活动;极点在其剪切活动线(含)外侧、同时在扩张活动线(含)左侧时(如P5—P9),活动类型要通过式7的计算(即比较fSfT的相对大小)来确定:对于极点P5,fS = fT=1.0,该界面发生扩张活动和剪切活动的可能性相同;对于极点P6,fS> 1.0 >fT,该界面发生剪切活动;对于极点P7,fS>fT> 1.0,该界面也发生剪切活动;对于极点P8,fT> 1.0 >fS,该界面发生扩张活动;对于极点P9,fT>fS> 1.0,该界面也发生扩张活动。为方便表述,这里简单认为极点P5—P9的内聚力CP和扩张强度TP相等,它们的剪切活动线和扩张活动线重合。

    广义破裂活动准则从破裂的物理本质出发,在给定应力状态下可以定量分析任意界面(先存薄弱面及假想界面)破裂活动的可能性及破裂活动的类型,理论适用于任何实际的地质体,并把经典的破裂准则统一在了一起。考虑到破裂活动与人类生活和生产密切相关,该准则具有广阔的应用前景。下面就此作一些简单的阐述。

    经典破裂准则包括Coulomb-Mohr剪破裂准则和Griffith破裂准则。

    Anderson应用Coulomb准则指出均匀介质发生剪破裂是方位优选的,且为一对共轭剪切破裂面,但没有阐述其原因。

    应用广义破裂活动准则和摩尔空间(童亨茂等,2014),可以很直观地解释均匀介质发生剪破裂是共轭的,而且是方位优选的:均匀介质在载荷作用下,当应力状态达到Coulomb应力状态时,摩尔空间与Coulomb-Mohr剪破裂线相切(上下有两个切点;图1)。此时,只有这两个切点所代表的界面达到了剪切活动的条件(fS = 1.0),而其它方位界面均未达到剪切活动的条件(极点均位于破裂线内侧,fS< 1.0),而且这两个界面同时达到破裂条件、剪切方向相反,两界面与最大主应力的夹角为±(45°−φ/2),与中间主应力轴平行(两切点位于σ1σ3摩尔大圆上),以共轭形式存在。

    Griffith破裂准则实际上分两部分,分别为张破裂准则和剪破裂准则(公式(4))。

    应用广义破裂活动准则和摩尔空间,可以很直观地解释均匀介质发生张破裂是方位优选的,即垂直于最小主应力轴方向:均匀介质在载荷作用下,当应力状态达到Griffith应力状态时,摩尔空间与扩张破裂线相切(只有一个切点,见图3)。此时,只有这一个切点所代表的界面达到了扩张活动的条件(fT = 1.0),而其他方位界面均未达到扩张活动的条件(极点均位于破裂线右侧,fT< 1.0),且此界面与最小主应力轴垂直。

    当破裂点位于纵坐标轴(含)右侧时,Griffith准则和Coulomb-Mohr准则的破裂包络线形态是类似的(图1),但当破裂点位于纵坐标轴左侧时则差异很大。

    图1所示,Griffith破裂包络线是一条以横坐标轴为对称轴的抛物线(Coulomb剪破裂线则是直线),与横坐标轴的交点为(−T,0),与纵坐标轴的交点(即截距)为(0,±2T)。从纵坐标轴向左,Griffith破裂包络线的斜率(即内摩擦系数)不断增加,并在点(−T,0)时达到无穷大,而实际岩石的内摩擦系数一般都小于1.0,因此,Griffith破裂包络线斜率大于1的段落没有实际的物理意义。另外,按照Griffith准则,岩石的内摩擦系数仅与正应力的变化有关,与岩性无关,这显然也存在问题。

    按照Griffith准则,C = 2T,这一结果同时也得到岩石力学实验结果的支持(Jaeger and Cook,1979);岩石的摩擦系数一般为0.5~0.8(Tong et al.,2010)。由于应力状态不能超越Griffith应力状态,在这些参数的取值情况下,应力摩尔圆不可能与纵坐标轴左侧的Coulomb剪破裂线相切,表明均匀介质中正应力为负值(即处于拉张状态)的界面不可能发生剪切破裂,同时也表明位于纵坐轴左侧的Coulomb剪破裂线没有实际的物理意义。

    σ3′ ≥ 0和σ3′ <0两种情况进行讨论。

    σ3′ ≥ 0(摩尔空间位于纵坐标轴右侧)时,所有界面在摩尔空间中极点均位于纵坐标轴(含)右侧,按照广义破裂准则,只可能发生剪切活动,不可能发生扩张活动(fT< 0),破裂活动的类型是确定的,是剪切活动。

    σ3′ <0(最小有效主应力为张应力),且界面在摩尔空间中极点位于纵坐标轴左侧(σn′ <0)时,按照广义破裂准则,界面既可能发生剪切活动(当界面的极点位于其剪切活动线(含)外侧时),也可能发生扩张活动(当界面的极点位于其扩张活动线(含)左侧时),破裂活动的类型取决于其fSfT值的相对大小。当应力状态、薄弱面的方位和力学性质(扩张强度、内聚力和摩擦系数)已知时,即可实现破裂活动类型的预测。

    广义破裂活动准则中界面(各种薄弱面及方位优选的界面)的破裂活动趋势因子虽然是一个静态参数,但可以分析预测各种界面潜在的破裂活动能力及演化趋势。应用童亨茂等(2014)的方法,利用摩尔空间技术方法分别确定给定界面在Griffith应力状态下的扩张活动趋势因子和Coulomb应力状态下的剪切活动趋势因子后,从大到小顺序排列的次序即能反映各个给定界面潜在的破裂活动能力。

    根据前面的分析,地质体产生扩张活动(包括其中的薄弱面和方位优选的界面)是由流体压力的不断增加引起的。在只考虑扩张活动的情况下,如图5所示,随着流体压力的增加,摩尔空间向左移动,fT最大的极点P1首先接触其扩张活动线,该界面发生扩张(图5b);随着流体压力的进一步增加,fT较小的极点P2接触其扩张活动线,该界面发生扩张(图5c),在此过程中,P1通过界面的扩张和扩展来保持其内部的流体压力平衡(即流体压力保持不变)。而当流体压力继续增大使得区域应力状态达到Griffith应力状态时(图5d),地质体产生垂直于最小主应力轴的扩张破裂,之后便产生应力降(其值为岩石的抗张强度T),使区域应力达到新的稳定状态,摩尔空间向右移到与纵坐标轴相切的位置(图5e)。此后,薄弱面的扩张活动停止,若继续有流体补充,垂直于最小主应力轴的扩张破裂会通过扩张和扩展来保持流体压力平衡,直到流体压力补充结束。而在整个过程中,fT最小的极点P3始终无法接触到其扩张活动线,因此该界面一直保持稳定。同样,为方便表述,简单认为极点P1、P2、P3的扩张强度TP相等,它们的扩张活动线重合。

    图  5  含多个薄弱面地质体在流体压力不断增大过程中扩张活动演化的摩尔空间图解
    T—岩石的抗张强度;2T—岩石的内聚力;TP—薄弱面的扩张强度;g—Griffith张破裂线;gP—薄弱面的扩张活动线a—初始应力状态下的摩尔空间图解;b—极点P1与扩张活动线接触应力状态下的摩尔空间图解;c—极点P2与扩张活动线接触应力状态下的摩尔空间图解;d—Griffith应力状态下的摩尔空间图解;e—应力降后区域应力状态下的摩尔空间图解
    Figure  5.  Mohr space diagram of the evolution of extension activities in geological bodies with multiple weak surfaces in the process of increasing fluid pressure
    (a) Mohr space diagram in the initial stress state; (b) Mohr space diagram under the stress state when pole P1 touches the extension activation line; (c) Mohr space diagram under the stress state when pole P2 touches the extension activation line; (d) Mohr space diagram under Griffith stress state; (e) Mohr space diagram of regional stress state after stress drop T is the tensile strength of the rock, 2T is the cohesion of the rock, TP is the extension strength of the weakness plane, g is the Griffith tensile fracture line, and gP is the extension activation line of the weakness plane.

    在只考虑剪切活动的情况下,对于有多个先存薄弱面的地质体剪切活动的递进变形过程,Tong and Yin(2011)已作了详细阐述,文中不再赘述。

    在同时考虑扩张活动和剪切活动的情况下,对于有多个薄弱面的地质体,不同薄弱面在递进变形过程中的破裂活动方式可能是不同的,取决于其在摩尔空间中的极点优先与其扩张活动线还是剪切活动性接触(即其fSfT值中哪个先达到1.0)。而当区域应力状态演化到Griffith应力状态或者Coulomb应力状态时,地质体最终会产生垂直于最小主应力轴的扩张破裂或者Coulomb破裂,具体会出现哪种破裂类型取决于整个摩尔空间优先与Griffith张破裂线还是Coulomb剪破裂线接触(即区域上的fSfT值中哪个先达到1.0)。

    从上面的分析可以看出,广义破裂活动准则可以实现对破裂活动演化过程的动态分析,突破了经典破裂准则静态分析的局限性,应用领域得到进一步拓展。

    理论上,广义破裂活动准则可以应用到与破裂活动相关的各种作用的分析研究中,如:① 压裂过程中压裂液流动途径预测(童亨茂等,2021a);② 压裂过程中形成的裂缝网络分布预测(童亨茂等,2021a);③ 压裂过程中套管变形风险预测(童亨茂等,2021b);④ 压裂过程中压窜风险预测(童亨茂等,2021a);⑤ 固体矿产开采过程中的突水风险预测;⑥ 为地震和滑坡预测提供基础理论模型;⑦ 隧道(或坑道)开挖时的稳定性分析和评价;⑧ 与破裂相关的大型工程地基安全性分析和评价等。具体的分析方法不在这里一一详述。

    与广义剪切活动一样,广义破裂活动准则隐含了两个假设条件:① 界面是一个平面;② 计算区域内的应力场是均匀的。但这两个假设条件并不影响该准则在实际地质体中的应用。在给定的精度要求下,只要把实际地质体分割成足够小的有限单元,就能达到广义破裂活动准则的应用条件。

    如自然界实际的断层面都是曲面,并非需要分割成无限小后才能看作是平面,而是综合考虑相关数据的误差情况和工作需要后,把断层面分割成若干段落后,每个段落的断面即可认为是平面,其内的应力场可以认为是均匀的,即可达到广义破裂活动准则的应用条件。这种有限分割,一方面满足了广义破裂活动准则的应用条件,另一方面可揭示断层不同段落、不同深度层次活动差异性的机理问题。

    再如在岩石圈应力状态下,由于重力的作用,主应力值随深度是不断变化的,应力场必然是非均匀的。但如果只从宏观上分析断层作用,在介质是均匀的情况下,应力场也可以被看作是均匀的(Anderson,1951),Tong等(Tong and Yin,2011Tong et al.,2014)按均匀应力场计算的砂箱内先存断层的活动趋势与模拟实验结果具有惊人的一致性,进一步支持了该认识。

    此外,破裂活动一旦发生,无论是剪切活动(断层作用)还是扩张活动,破裂面附近均会产生应力降,形成局部应力场(Yin,1994Morley et al.,2007Gudmundsson et al.,2010Tong et al.,2014),其中断层作用形成的局部应力场影响范围相对较大。另外,在破裂活动过程中,伴随着破裂面的不断扩展,地质体的结构也会发生变化,尤其是断层作用,还会伴随地质体的变形和变位。因此,在破裂活动过程中,局部应力场和地质体结构时刻在发生变化。在此过程中,应力和应变是紧密关联的两个方面,即应力导致应变,应变又引起应力的变化,构成了复杂的“应力−应变体系”(stress-strain system;童亨茂等,2015)。但是,目前破裂活动及其伴随的构造变形(如断层作用往往伴随褶皱变形,形成断层相关褶皱),以及在递进变形过程中,应力场的演化特征和规律研究还十分初步,是构造地质学的前沿研究课题,而广义破裂活动准则是联系“应力”和“破裂活动”的桥梁。随着观测资料(如高精度三维地震资料、GPS观测数据等)的不断精细,变形结果(构造几何特征)的可观测性和变形过程(构造的运动学过程)的可解析性不断增强,广义剪切破裂准则理论有望在递进变形过程中应力场演化这一研究领域发挥重要作用,也可为解决地震预测等重大难题提供理论支撑。

    文章从产生破裂的物理本质出发,通过界面(薄弱面和假想界面)破裂活动趋势因子的分析,建立了界面发生破裂活动的准则。该准则从理论上可以定量判断任意介质、在任意三轴应力状态下、任意方位界面发生破裂活动的可能性及破裂活动类型(扩张活动或剪切活动),并把经典的破裂准则(Coulomb-Mohr准则和Griffith准则及Byerlee律为代表的滑动摩擦律)统一在了一起,把破裂(fracturing)扩展到破裂活动(fracturing activity),因此该准则被称之为广义破裂活动准则。

    广义破裂活动准则揭示了Coulomb共轭剪破裂和Griffith张破裂方位优选的物理本质,并指出了其中存在的问题,可以定量评价不同界面破裂活动的演化趋势。

    广义破裂活动准则适用于实际的地质体,在与破裂相关的资源勘探开发、灾害(如地震与滑坡)预测和防治中具有广阔的应用前景。

  • 图  1  Coulomb-Mohr准则与Griffith准则图解

    φ为岩石的内摩擦角;T为岩石的抗张强度;内聚力C=2T

    Figure  1.  Diagram of Coulomb-Mohr criterion and Griffith criterion

    φ is the internal friction angle of the rock, T is the tensile strength of the rock and the cohesion C=2T

    图  2  界面的空间方位定义(以3个主应力为坐标轴)

    σ1σ2σ3分别为最大、中间和最小主应力;θ为界面与σ1的夹角;α为界面在σ2σ3平面上的交线与σ3的夹角

    Figure  2.  The spatial orientation definition of the interface (with 3 principal stresses as axes)

    σ1, σ2 and σ3 are the maximum, intermediate and minimum principal stresses respectively, θ is the angle between the interface and σ1, and α is the angle between σ3 and the interface's intersection line with the σ2σ3 plane.

    图  3  Griffith应力状态下的摩尔空间图解

    T为岩石的抗张强度;2T为岩石的内聚力;TP为薄弱面的扩张强度;g为Griffith张破裂线;gP为薄弱面的扩张活动线;θ 为界面与 σ 1的夹角; α 为界面在 σ 2σ 3平面上的交线与 σ 3的夹角

    Figure  3.  Diagram of the Mohr space under Griffith stress state

    T is the tensile strength of the rock, 2T is the cohesion of the rock, TP is the extension strength of the weakness plane, g is the Griffith tensile fracture line, gP is the extension activation line of the weakness plane, θ is the angle between the interface and σ 1, and α is the angle between σ 3 and the interface's intersection line with the σ 2- σ 3 plane.

    图  4  给定应力状态下摩尔空间中不同极点与破裂活动线的关系图解

    T为岩石的抗张强度;2T为岩石的内聚力;C 为介质的内聚力;C P为界面的内聚力TP为薄弱面的扩张强度;g为Griffith张破裂线;gP为薄弱面的扩张活动线;P1—P9为不同薄弱面在摩尔空间中的极点位置

    Figure  4.  Diagram of relationship between different poles and fracturing activity line in Mohr space under given stress state

    T is the tensile strength of the rock, 2T is the cohesion of the rock, C is the cohesion of the medium, CP is the cohesion of the interface, TP is the extension strength of the weakness plane, g is the Griffith tensile fracture line, gP is the extension activation line of the weakness plane, and P1–P9 are the pole positions of different weakness planes in Mohr-space.

    图  5  含多个薄弱面地质体在流体压力不断增大过程中扩张活动演化的摩尔空间图解

    T—岩石的抗张强度;2T—岩石的内聚力;TP—薄弱面的扩张强度;g—Griffith张破裂线;gP—薄弱面的扩张活动线a—初始应力状态下的摩尔空间图解;b—极点P1与扩张活动线接触应力状态下的摩尔空间图解;c—极点P2与扩张活动线接触应力状态下的摩尔空间图解;d—Griffith应力状态下的摩尔空间图解;e—应力降后区域应力状态下的摩尔空间图解

    Figure  5.  Mohr space diagram of the evolution of extension activities in geological bodies with multiple weak surfaces in the process of increasing fluid pressure

    (a) Mohr space diagram in the initial stress state; (b) Mohr space diagram under the stress state when pole P1 touches the extension activation line; (c) Mohr space diagram under the stress state when pole P2 touches the extension activation line; (d) Mohr space diagram under Griffith stress state; (e) Mohr space diagram of regional stress state after stress drop T is the tensile strength of the rock, 2T is the cohesion of the rock, TP is the extension strength of the weakness plane, g is the Griffith tensile fracture line, and gP is the extension activation line of the weakness plane.

  • [1] ANDERSON E M, 1951. The dynamics of faultingand dyke formation with applications to Britain[M]. 2nd ed. Edinburgh: Oliver and Boyd.
    [2] BAILEY I W, BEN-ZION Y, 2009. Statistics of earthquake stress drops on a heterogeneous fault in an elastic half-space[J]. Bulletin of the Seismological Society of America, 99(3): 1786-1800. doi: 10.1785/0120080254
    [3] BOTT M H P, 1959. The mechanics of oblique slip faulting[J]. Geological Magazine, 96(2): 109-117. doi: 10.1017/S0016756800059987
    [4] BYERLEE J, 1978. Friction of rocks[J]. Pure and AppliedGeophysics, 116(4-5): 615-626.
    [5] CÉLÉRIER B, 2008. Seeking Anderson’s faulting in seismicity: a centennial celebration[J]. Reviews of Geophysics, 46(4): RG4001.
    [6] GRIFFITH A A, 1921. The phenomena of rupture and flow in solids[J]. Philosophical Transactions of the Royal SocietyA: Mathematical, Physical and Engineering Sciences, 221(582-593): 163-198.
    [7] GUDMUNDSSON A, SIMMENES T H, LARSEN B, et al. , 2010. Effects of internal structure and local stresses on fracture propagation, deflection, and arrest in fault zones[J]. Journal of Structural Geology, 32(11): 1643-1655. doi: 10.1016/j.jsg.2009.08.013
    [8] JAEGER J C, COOK N G W, 1979. Fundamentals of rock mechanics[M]. 3rd ed. London: Chapman and Hall.
    [9] MCKENZIE D P, 1969. The relation between fault plane solutions for earthquakes and the directions of the principal stresses[J]. Bulletin of the Seismological Society of America, 59(2): 591-601. doi: 10.1785/BSSA0590020591
    [10] MORLEY C K, GABDI S, SEUSUTTHIYAK, 2007. Fault superimposition and linkage resulting from stress changes during rifting: examples from 3D seismic data, Phitsanulok Basin, Thailand[J]. Journal of Structural Geology, 29(4): 646-663. doi: 10.1016/j.jsg.2006.11.005
    [11] MORRIS A, FERRILL D A, HENDERSON D B, 1996. Slip-tendency analysis and fault reactivation[J]. Geology, 24(3): 275-278. doi: 10.1130/0091-7613(1996)024<0275:STAAFR>2.3.CO;2
    [12] TONG H M, MENG L J, CAI D S, et al. , 2009. Fault formation and evolution in rift basins: sandbox modeling and cognition[J]. Acta Geologica Sinica, 83(6): 759-774. (in Chinese with English abstract)
    [13] TONG H M, CAI D S, WU Y P, et al. , 2010. Activity criterion of pre-existing fabrics in non-homogeneous deformation domain[J]. Science China Earth Sciences, 53(8): 1115-1125. doi: 10.1007/s11430-010-3080-6
    [14] TONG H M, YIN A, 2011. Reactivation tendency analysis: a theory for predicting the temporal evolution of preexisting weakness under uniform stress state[J]. Tectonophysics, 503(3-4): 195-200. doi: 10.1016/j.tecto.2011.02.012
    [15] TONG H M, KOYI H, HUANG S, et al. , 2014. The effect of multiple pre-existing weaknesses on formation and evolution of faults in extended sandbox models[J]. Tectonophysics, 626: 197-212. doi: 10.1016/j.tecto.2014.04.046
    [16] TONG H M, WANG J J, ZHAO H T, et al. , 2014. Mohr space and its application to the activation prediction of pre-existing weakness[J]. Science China Earth Sciences, 57(7): 1595-1604. doi: 10.1007/s11430-014-4860-1
    [17] TONG H M, CHEN Z L, LIU R X, 2015. Generalized shear activation criterion[J]. Chinese Journal of Nature, 37(6): 441-447. (in Chinese with English abstract)
    [18] TONG H M, LIU Z P, ZHANG H X, et al. , 2021. Theory and method of temporary macrofracture plugging to prevent casing deformation in shale gas horizontal wells[J]. Natural Gas Industry, 41(5): 92-100. (in Chinese with English abstract)
    [19] TONG H M, ZHANG P, ZHANG H X, et al. , 2021. Geomechanical mechanisms and prevention countermeasures of casing deformation in shale gas horizontal wells[J]. Natural Gas Industry, 41(1): 189-197. (in Chinese with English abstract)
    [20] TWISS R J, MOORES E M, 1992. Structural geology[M]. San Francisco: W. H. Freeman & Co. : 532.
    [21] WALLACE R E, 1951. Geometry of shearing stress and relation to faulting[J]. The Journal of Geology, 59(2): 118-130. doi: 10.1086/625831
    [22] XU K L, ZHUZ C, 1989. Structural geology[M]. 2nd ed. Beijing: Geology Press: 270. (in Chinese)
    [23] YIN A, 1994. Mechanics of monoclinal systems in the Colorado plateau during the Laramide orogeny[J]. Journal of Geophysical Research: Solid Earth, 99(B11): 22043-22058. doi: 10.1029/94JB01408
    [24] ZOBACK M D, 2007. Reservoir geomechanics[M]. New York: Cambridge University Press.
    [25] 童亨茂, 孟令箭, 蔡东升, 等, 2009. 裂陷盆地断层的形成和演化: 目标砂箱模拟实验与认识[J]. 地质学报, 83(6): 759-774. doi: 10.3321/j.issn:0001-5717.2009.06.002
    [26] 童亨茂, 王建君, 赵海涛, 等, 2014. “摩尔空间”及其在先存构造活动性预测中的应用[J]. 中国科学: 地球科学, 44(9): 1948-1957.
    [27] 童亨茂, 陈正乐, 刘瑞珣, 2015. 广义剪切活动准则[J]. 自然杂志, 37(6): 441-447.
    [28] 童亨茂, 刘子平, 张宏祥, 等, 2021a. 暂堵大裂缝防治页岩气水平井套管变形的理论与方法[J]. 天然气工业, 41(5): 92-100.
    [29] 童亨茂, 张平, 张宏祥, 等, 2021b. 页岩气水平井开发套管变形的地质力学机理及其防治对策[J]. 天然气工业, 41(1): 189-197.
    [30] 徐开礼, 朱志澄, 1989. 构造地质学[M]. 2版. 北京: 地质出版社: 270.
  • 期刊类型引用(1)

    1. 童亨茂. 李四光对地质力学(Geomechanics)的开创性贡献——纪念李四光诞辰135周年. 地质力学学报. 2025(01): 1-7 . 本站查看

    其他类型引用(0)

  • 加载中
图(5)
计量
  • 文章访问数:  665
  • HTML全文浏览量:  189
  • PDF下载量:  127
  • 被引次数: 1
出版历程
  • 收稿日期:  2023-11-06
  • 修回日期:  2023-12-29
  • 录用日期:  2024-01-12
  • 预出版日期:  2024-02-02
  • 刊出日期:  2024-02-01

目录

/

返回文章
返回