留言板

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

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

上拱力背景下正断裂剖面形态及极限应力状态研究

宋随宏 陈书平 何明宇

宋随宏, 陈书平, 何明宇, 2012. 上拱力背景下正断裂剖面形态及极限应力状态研究. 地质力学学报, 18 (2): 149-157.
引用本文: 宋随宏, 陈书平, 何明宇, 2012. 上拱力背景下正断裂剖面形态及极限应力状态研究. 地质力学学报, 18 (2): 149-157.
SONG Sui-hong, CHEN Shu-ping, HE Ming-yu, 2012. THE INVESTIGATION OF NORMAL FAULT UNDER UPWELLING FORCE AND ITS CRITICAL STRESS STATE. Journal of Geomechanics, 18 (2): 149-157.
Citation: SONG Sui-hong, CHEN Shu-ping, HE Ming-yu, 2012. THE INVESTIGATION OF NORMAL FAULT UNDER UPWELLING FORCE AND ITS CRITICAL STRESS STATE. Journal of Geomechanics, 18 (2): 149-157.

上拱力背景下正断裂剖面形态及极限应力状态研究

基金项目: 

国家自然科学基金项目 41172124

详细信息
    作者简介:

    宋随宏(1990-), 男, 中国石油大学(北京)地球科学学院地质工程专业本科生。E-mail:air.sunsuihong@yahoo.com

  • 中图分类号: P554

THE INVESTIGATION OF NORMAL FAULT UNDER UPWELLING FORCE AND ITS CRITICAL STRESS STATE

  • 摘要: 上拱力构造背景下形成的正断裂及断裂极限应力状态研究在构造解释和油气运移成藏中具有重要的意义。在抛物线型莫尔包络线理论的基础上,通过合理的假设和数学推导以及理论模拟,进一步探讨上拱力构造背景影响下形成的正断裂剖面形态和断裂的极限应力状态,并定量给出函数表达式。研究结果表明,断裂在均质岩体中由浅到深,倾角从竖直开始由大变小,呈铲状断裂,而变化的速率只与地层物性有关,泥岩比砂岩变化快。断裂分布的深度与上拱力的大小成一次线性关系,上拱力越大,断裂分布越浅。砂泥岩互层地层中发育有典型的坡坪式正断裂。研究结果不仅与普遍认识相一致,而且可以进一步在实际应用中预测断裂的存在,定量解释断裂形态,求取上拱力和断裂极限应力状态。

     

  • 岩石断裂既是地壳构造变形和演化的重要产物,也是一种重要的构造作用与构造变形方式,它对油气运移和油气藏形成以及寻找金属矿床等很多方面都具有重要意义[1]

    关于岩石断裂的机理和准则近百年来国内外众多专家和学者都做过很多探索。经典断裂力学理论中的一个根本假设是将岩石断裂轨迹视为直线型平面模型,然而现场实测和实验观测均表明,无论是在晶粒尺度还是在断层尺度上,岩石的断裂面都是非常不规则和粗糙的,难以用一个平直面来近似模拟[2~3],显示出了经典断裂力学理论的局限性。除此之外常用的断裂准则还包括:库仑断裂准则、水平直线型莫尔包络线理论、抛物线型莫尔包络线理论、格里菲斯断裂准则[4]。以卡门实验为例,抛物线型莫尔包络线理论更符合实际情况[5]

    关于地壳中岩石断裂的剖面形态,国内外学者都做过很多研究和探讨,总结出了常见的主要断裂样式:平直型断裂(层)、铲式断裂(层)和坡坪式断裂(层)[6]。这些认识符合自然界实际情况,得到了普遍的认可,能定性地解释和预测一些现象,但对断裂产生机理的定量解释和断裂形态的定量表达却很少[7~8]。本文将在抛物线型莫尔包络线理论基础上,定量地探讨上拱力构造背景下产生正断裂的形态、断裂极限应力及其与上拱力、地层物性之间的关系。这将对于精确预测解释断裂、认识盆地、认识上拱构造以及预测油气藏、分析油气运移等有重大意义。

    当莫尔包络线以抛物线来近似表示(见图 1)时,其方程可写为[5]

    图  1  抛物线型莫尔包络线[5]
    τ, σ—断裂面上的临界剪应力和临界正应力,MPa;σ1, σ3—发生断裂时最大和最小主应力,MPa;σm—平均应力(σm=(σ1+σ3)/2),MPa;σI—岩石在各向等值拉伸条件下的抗张断裂极限,MPa;θ2—岩石剪裂角,(°)
    Figure  1.  The Parabolic Mohr failure envelope
    τ2=K2σI(σIσ)
    (1)

    式中:σI为岩石在各向等值拉伸条件下的抗张断裂极限,MPa;K为岩石粘结强度,MPa;τσ分别为断裂面上的临界剪应力和临界正应力,MPa。本文中规定,正值表示拉力,负值表示压力。

    以主应力来表达,抛物线型莫尔包络线准则可写成[5]

    (σ1σ3)2+2K2σI(σ1+σ3)=K2(4K2σI)
    (2)

    式中:σ1σ3分别为发生断裂时最大和最小主应力,MPa。

    按照这一准则,剪裂角可表达成[5]

    tanθ=2σIK(1K22σIσ1+σ32σI)12θ=arctan2σIK(1K22σIσ1+σ32σI)12]
    (3)

    随着σ1σ3以及KσI的变化,剪裂角可由0°变到45°。从(3)式可以看出,剪裂角不仅与物性参数(KσI)有关,而且与各点极限应力状态有关。即使在同一种岩石中,由于各处极限应力状态不同,也可以形成不同的剪裂角[5]。下面以此理论为基础讨论上拱力构造背景下正断裂剖面形态和极限应力状态。

    研究表明,地壳深部一般应力状态为:

    σ1=σ2=σ3=ρgh
    (4)

    式中:ρ表示岩体的密度,kg/cm3g为重力加速度,m/s2h表示参考点的深度,m;σ1σ2σ3分别表示最大主应力、中间主应力、最小主应力,MPa。

    裂陷盆地中正断裂的形成常与盆地深部岩浆或地幔柱上升等有关[9~12]。在上拱力构造背景下,岩石除受到一般应力的作用,还受到由上拱构造引起的附加水平拉力和上拱力共同作用。这种情况下,假设岩体受到的上拱力恒定,为σy2,附加水平拉力为σx2,埋深为h点处岩体竖直方向一般应力为σy1,水平方向一般应力为σx1,埋深为h点处岩体竖直方向总应力为σy,水平方向总应力为σx,则有:

    σy2=ρgh2
    (5)
    σy1=ρgh
    (6)
    σx1=ρgh
    (7)
    σy=σy1+σy2
    (8)
    σx=σx1+σx2
    (9)

    式中:σy2为上拱力,MPa;h2为上拱力以一般应力形式表示时对应的深度,m,由σy2导出。

    当该点处岩体承受的竖直和水平方向总应力达到极限应力状态而发生正断裂时,根据抛物线型莫尔包络线理论,将(2)式变形整理可得:

    σ1=σ3K2σI±2K1σ3σI
    (10)

    在正断裂中,由于σxσy,所以有:

    σ1=σx,σ3=σy
    (11)

    联立(5)—(11)式,整理可得深度为h点处岩体发生正断裂的极限应力状态:

    σy=ρg(h+h2)
    (12)
    σx=ρg(h+h2)K2σI+2K1ρg(h+h2)σI
    (13)

    分析(13)式中水平方向极限应力并结合水平方向一般应力大小,σx由两部分组成,分别是地壳水平方向一般应力σx1和由上拱构造引起的附加水平拉力σx2,它们都随深度而变化。其中:

    σx2=ρgh2K2σI+2K1ρg(h+h2)σI
    (14)

    当岩体所承受的总应力达到断裂极限应力状态发生断裂(见图 2)时:

    图  2  岩体正断裂剖面数学分析
    l—衡量地面水平长度的量,m;h—参考点深度,m;θ2—岩石剪裂角,(°);σx, σy—深度为h点处岩体水平方向和竖直方向总应力,MPa
    Figure  2.  Mathematical analysis of geometric shape of normal fault of rock mass
    dldh=tanθ2
    (15)

    式中l为衡量地面水平长度的量,m。

    联立(3)、(11)、(12)、(13)、(15)式,可求得断裂剖面的数学表达式为:

    dldh=1+K24σ2I+ρg(h+h2)σIKσI1+ρg(h+h2)σIK2σI1+ρg(h+h2)σIKσI1+ρg(h+h2)σI
    (16)

    式(16)即为在上拱力背景下,岩体受到上拱力σy2和由上拱构造引起的附加水平拉力σx2而达到极限应力状态发生正断裂时,断裂剖面的数学表达式。

    分析(16)式可知,断裂的剖面形态和极限应力状态只与地层物性参数KσIρ有关;而断裂的分布深度则与上拱力大小相关。

    根据对式(16)的分析,下文将研究上拱力大小为0,即在水平拉张背景下,断裂形态和极限应力状态。

    中国石油大学(北京)对兰聊断裂带附近砂岩进行测定及计算得到:K=50 MPa、σI=16 MPa、ρ=2.32 g/cm3。用此砂岩建立各向同性均质砂岩体模型,模拟水平拉张背景下正断裂剖面形态和极限应力状态。根据理论分析,并将具体数值代入(12)、(13)、(14)及(16)式,可得到在水平拉张背景下,发生断裂时断裂的极限应力状态为:

    σy=0.02273h
    (17)
    σx=0.023h+14h+10000156.25
    (18)

    要达到断裂极限应力状态需要提供的附加水平拉力为:

    σx2=14h+10000156.25
    (19)

    断裂剖面函数表达式为:

    dldh=0.0014h3.1×0.0014h+1+3.441.560.0014h3.1×0.0014h+1+1
    (20)

    根据所求得的断裂剖面函数表达式(20)画出断裂剖面(见图 3),可以看出:断裂从h0=6169 m深度点开始出现,随深度增加,断裂由陡变缓,其内摩擦角由90°减小,无限接近0°,即剪裂角由0°开始,无限接近45°方向,最终在地下不够提供断裂极限应力的某点消失,呈铲式,对应于抛物线型莫尔包络线中应力圆与抛物线相切于两点的情况;对于从地表到地下h0=6169 m深度点之间,从函数角度来讲,该表达式无解,显示出了该表达式的缺陷。而从实际地质力学角度来讲,存在竖直形态的断裂,其剪裂角为0°,对应于抛物线型莫尔包络线中应力圆与抛物线相切于其顶点的情况;深度为h0的点,则对应于应力圆与抛物线顶点相切的极限情况。在本文中,将深度为h0这种应力圆刚好与抛物线顶点相切、断裂开始出现弯曲的深度命名为起弯深度,将这个深度点命名为起弯断点(见图 3)。

    图  3  正断裂剖面及对应的抛物线型莫尔包络线模型
    l—衡量地面水平长度的量,m;h0—起弯深度,m
    Figure  3.  The geometric shape of normal fault and parabolic Mohr failure envelope models

    根据对(16)式的分析,断裂的分布深度与上拱力大小有关,上拱力越大,断裂分布越浅。下文将在上拱力(σy2=-ρgh2)为-140.26 MPa(此值可使下文均质砂岩中的断裂起弯深度为0,其计算方法后文将阐述)的前提下模拟不同物性地层产生的不同断裂形态和极限应力状态。

    3.2.1   均质砂岩中正断裂的剖面形态和极限应力状态模拟

    根据对兰聊断裂带附近砂岩的测定及计算结果(K=50 MPa、σI=16 MPa、ρ=2.32 g/cm3)建立各向同性均质砂岩体模型,模拟此情况下正断裂剖面形态和极限应力状态。

    将具体数值代入公式(12)、(13)、(16),可得此情况下极限应力状态和断裂剖面形态表达式:

    σy=0.022736(h2+6169)
    (21)
    σx=0.022736h+100+0.14×(h+6169)296.5
    (22)

    断裂剖面的函数表达式为:

    dldh=0.0014h0.14h+94.1+12.081.560.0014h0.14h+94.1+9.64
    (23)

    由所求得的断裂剖面函数表达式画出断裂剖面(见图 4a),可以看出:在此上拱力的背景下,断裂的起弯深度恰好为0,即断裂恰好从地表开始角度发生变化,随深度增加,断裂剪裂角由0°开始,无限接近45°方向,最终在地下不够提供其断裂极限应力的某点消失。这与实际情况中铲式断裂能够很好地吻合,能够很好定量解释其形成机制。

    图  4  上拱力构造背景下岩石极限应力及正断裂剖面形态
    Figure  4.  Critical stress state of rock and geometric shape of fault under tectonic setting of upwelling force
    3.2.2   均质泥岩中正断裂的剖面形态和极限应力状态模拟

    取兰聊断裂带附近泥岩,测定并计算其参数,建立均质泥岩体模型。应用同样的步骤可以计算出其极限应力状态和断裂剖面形态的函数表达式,并画出断裂剖面(见图 4b),其形态与上述均质砂岩中断裂类似,不同的是该剖面起弯深度为-152 m,表示断裂从地上152 m开始剪裂角发生变化,而实际情况是断裂剖面只为其地下部分。此外,泥岩断裂剖面比砂岩剖面剪裂角变化更快,显示为更平缓。

    3.2.3   均质砂泥岩互层断裂剖面形态和极限应力状态模拟

    在上述均质砂岩模型中插入一薄层均质泥岩,建立砂泥岩互层模型。此时,断裂的极限应力状态和断裂剖面为同等情况下均质砂岩和泥岩模型中应力状态和断裂的组合(见图 4c),可以看出:该砂泥岩互层模型中,砂岩段断裂面相对较陡,泥岩段断裂面较缓,由深至浅组合而形成典型的“陡—缓—陡”的坡坪式正断裂,能够很好地定量解释实际情况中的坡坪式断裂的形成机制和其形态。

    3.2.4   实际地层断裂剖面形态和极限应力状态模拟

    实际情况下,包括KσIρ在内的各项地层参数都发生变化。资料显示,这3个参数都随地层深度增加而增大,并且增大的速度越来越慢。通过近似模拟可得出:在此上拱力背景下,实际地层的断裂剖面形态与均质岩体模型中的形态类似,呈铲式断裂。

    在已知KσIρ三个参数随深度变化的关系时,可定量得出实际地层断裂的极限应力状态和断裂剖面的函数表达式。

    从数学角度分析(16)式可知,断裂分布的深度或起弯深度与上拱力有一次线性关系。

    由于起弯断点对应于应力圆恰好与抛物线型莫尔包络线相切于其顶点的情况(见图 3),且:

    tanθ>0
    (24)

    所以可联立(3)、(11)、(12)、(13)、(24)式,求解得出:

    h>1ρg(σy2+K2σIσI)
    (25)

    进而可求得起弯深度h0:

    h0=1ρg(σy2+K2σIσI)
    (26)

    式(26)即可表示起弯深度随上拱力的变化关系。

    因此,对于在上拱力背景下产生的断裂,断裂形态只与地层物性有关,断裂分布的深度或起弯深度与上拱力呈一次线性关系,上拱力绝对值越大,断裂分布深度越浅,反之则越深。当上拱力为0时,即为水平拉张背景下产生的断裂(见图 3);当上拱力达到σy2=σI-K2/σI时,断裂起弯深度为0,断裂刚好从地表开始发生弯曲(见图 4a图 5

    图  5  起弯深度随上拱力变化关系
    σy2—上拱力,MPa;h—深度,m
    Figure  5.  The relationship of fault's depth and upwelling force

    在上拱力背景下,假设上拱力恒定,均质地层受到一般应力、上拱力和由上拱构造引起的附加水平拉力,当上拱力和附加水平拉力达到一定值时,地层处于极限应力状态发生正断裂,断裂极限应力状态和剖面函数表达式可以定量得出。这将对于断裂解释预测、极限应力状态和上拱力以及附加水平拉力求取、油气藏预测及分析具有重大意义。

    断裂由浅到深可以分为两个部分:从地表到起弯断点之间,断裂呈现竖直形态;起弯断点以下,随深度增加,断裂由陡变缓,剪裂角由0°无限趋近45°,一直延伸至不够提供其断裂所要求的极限应力的地下某点消失,呈现出铲式断裂的形态。这与实际情况中对铲式断裂的认识相一致,能够很好地定量解释其形成机制。

    断裂的形态只与地层物性有关,泥岩中产生的断裂剪裂角比砂岩角度变化快;断裂分布的深度与上拱力呈一次线性关系,上拱力绝对值越大,断裂分布深度越浅。

    当在较硬均质砂岩中夹入较软泥岩层时,可形成典型的“陡—缓—陡”的坡坪式正断裂,断裂在硬岩层中较陡,软岩层中较缓。断裂剖面形态的函数表达式和极限应力状态亦可定量表达出来。能够很好地解释和预测坡坪式断裂的形态及其形成机理。

  • 图  1  抛物线型莫尔包络线[5]

    τ, σ—断裂面上的临界剪应力和临界正应力,MPa;σ1, σ3—发生断裂时最大和最小主应力,MPa;σm—平均应力(σm=(σ1+σ3)/2),MPa;σI—岩石在各向等值拉伸条件下的抗张断裂极限,MPa;θ2—岩石剪裂角,(°)

    Figure  1.  The Parabolic Mohr failure envelope

    图  2  岩体正断裂剖面数学分析

    l—衡量地面水平长度的量,m;h—参考点深度,m;θ2—岩石剪裂角,(°);σx, σy—深度为h点处岩体水平方向和竖直方向总应力,MPa

    Figure  2.  Mathematical analysis of geometric shape of normal fault of rock mass

    图  3  正断裂剖面及对应的抛物线型莫尔包络线模型

    l—衡量地面水平长度的量,m;h0—起弯深度,m

    Figure  3.  The geometric shape of normal fault and parabolic Mohr failure envelope models

    图  4  上拱力构造背景下岩石极限应力及正断裂剖面形态

    Figure  4.  Critical stress state of rock and geometric shape of fault under tectonic setting of upwelling force

    图  5  起弯深度随上拱力变化关系

    σy2—上拱力,MPa;h—深度,m

    Figure  5.  The relationship of fault's depth and upwelling force

  • [1] 谭凯旋, 谢焱石, 陈广浩, 等.岩石断裂作用的复杂性和混沌动力学[J].大地构造与成矿学, 2003, 27(4):372~377. http://www.cnki.com.cn/Article/CJFDTOTAL-DGYK200304007.htm

    TAN Kai-xuan, XIE Yan-shi, CHEN Guang-hao, et al. Complexity and chaotic dynamics of rock faulting[J]. Geotectonica Et Metallogenia, 2003, 27(4):372~377. http://www.cnki.com.cn/Article/CJFDTOTAL-DGYK200304007.htm
    [2] Perie P J. Determination of fracture mechanism by microscopic observation of crack[J]. International Journal of Rock Mechanics and Mining Sciences and Geomechanics Abstracts, 1991, 28(1):83~84. doi: 10.1016/0148-9062(91)93236-Y
    [3] 谢和平.大理岩微观断裂的分形(fractal)模型研究[J].科学通报, 1989, 34(5):365~368. http://www.cnki.com.cn/Article/CJFDTOTAL-KXTB198905013.htm

    XIE He-ping. The investigation of fractal model of micro-fracture of marble[J]. Chinese Science Bulletin, 1989, 34(5):365~368. http://www.cnki.com.cn/Article/CJFDTOTAL-KXTB198905013.htm
    [4] 朱志澄.构造地质学[M].武汉:中国地质大学出版社, 1999:1~226.

    ZHU Zhi-cheng. Structural geology[M]. Wuhan:China University of Geosciences Press, 1999:1~226.
    [5] 陈庆宣, 王维襄, 孙叶, 等.岩石力学与构造应力场分析[M].北京:地质出版社, 1998:1~241.

    CHEN Qing-xuan, WANG Wei-xiang, SUN Ye, et al.The investigation of rock mechanics and tectonic stress field[M].Beijing:Geological Publishing House, 1998. 241.
    [6] 陆克政.含油气盆地分析[M].东营:石油大学出版社, 2001:1~424.

    LU Ke-zheng. The analysis of petroliferous basin[M]. Dongying:China University of Petroleum Press, 2001:1~424.
    [7] Suppe J, Chou G T, Hook S C. Rates of folding and faulting determined from growth strata[C]//McClay K R. Thrust tectonics. London:Chapman & Hall, 1992:105~121. doi: 10.1007/978-94-011-3066-0_9
    [8] 庄培仁, 常志忠.断层构造研究[M].北京:地质出版社, 1996:1~346.

    ZHUANG Pei-ren, CHANG Zhi-zhong. The investigation of fault structures[M]. Beijing:Geological Publishing House, 1996:1~346.
    [9] 李德生.渤海湾含油气盆地的地质和构造特征[J].石油学报, 1980, 1(1):7~20. http://www.cnki.com.cn/Article/CJFDTOTAL-SYXB198001002.htm

    LI De-sheng. Geology and structural characteristics of the Bohai Bay[J]. Acta Petrolei Sinica, 1980, 1(1):7~20. http://www.cnki.com.cn/Article/CJFDTOTAL-SYXB198001002.htm
    [10] 邓晋福, 赵海岭, 莫宣学, 等.中国大陆根-柱构造——大陆动力学的钥匙[M].北京:地质出版社, 1996:1~110.

    DENG Jin-fu, ZHAO Hai-ling, MO Xuan-xue, et al. Continental roots plume tectonics of China:Key to the continental dynamics[M]. Beijing:Geological Publishing House, 1996:1~110.
    [11] 牛树银, 胡华斌, 毛景文, 等.鲁西地区地质构造特征及其形成机制[J].中国地质, 2004, 31(1):34~39. http://www.cnki.com.cn/Article/CJFDTOTAL-DIZI200401003.htm

    NIU Shu-gen, HU Hua-bin, MAO Jing-wen, et al. Structure in western Shandong and its genetic mechanism[J]. Geology in China, 2004, 31(1):34~39. http://www.cnki.com.cn/Article/CJFDTOTAL-DIZI200401003.htm
    [12] 牛树银, 邵济安, 孙爱群, 等.华北东部盆山耦合与内生成矿作用[J].大地构造与成矿学, 2006, 30(3):331~342. http://www.cnki.com.cn/Article/CJFDTOTAL-DGYK200603008.htm

    NIU Shu-gen, SHAO Ji-an, SUN Ai-qun, et al. Coupling relationship between basin-mountain and nedogenic mineralization in the eastern part of north China[J]. Geotectonica et Metallogenia, 2006, 30(3):331~342. http://www.cnki.com.cn/Article/CJFDTOTAL-DGYK200603008.htm
  • 加载中
图(5)
计量
  • 文章访问数:  305
  • HTML全文浏览量:  173
  • PDF下载量:  11
  • 被引次数: 0
出版历程
  • 收稿日期:  2012-01-22
  • 刊出日期:  2012-06-01

目录

/

返回文章
返回