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: HE Feng, HU Zhi-jun, WU Shu-ren, 2004. GRAY PREDICTION MODEL OF DEFORMATION OF THE HUANGTUPO LANDSLIDE. Journal of Geomechanics, 10 (1): 51-56.

GRAY PREDICTION MODEL OF DEFORMATION OF THE HUANGTUPO LANDSLIDE

More Information
  • The essence of the gray prediction model is to fit the known data sequence with an exponent form to a curve, and then extend the curve to the future in order to predict the unknown data. This paper establishes a GM(1,1) model for prediction of the horizontal displacement of the Huangtupo landslide.Actual application shows that the model is characterized by a high precision and wide applicability.

     

  • 岩石断裂既是地壳构造变形和演化的重要产物,也是一种重要的构造作用与构造变形方式,它对油气运移和油气藏形成以及寻找金属矿床等很多方面都具有重要意义[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°,一直延伸至不够提供其断裂所要求的极限应力的地下某点消失,呈现出铲式断裂的形态。这与实际情况中对铲式断裂的认识相一致,能够很好地定量解释其形成机制。

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

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

  • Relative Articles

    TONG Hengmao, ZHANG Hongxiang, HOU Quanlin, CHEN Zhengle, HOU Guiting. 2024: Generalized fracturing activation criteria. Journal of Geomechanics, 30(1): 3-14. doi: 10.12090/j.issn.1006-6616.2023180
    ZHANG Di, LI Jiacun, WU Zhonghai, LIU Shaotang, LU Yan. 2021: Using terrestrial LiDAR to accurately measure the microgeomorphologic geometry of active fault: A case study of fault scarp on the Maoyaba fault zone. Journal of Geomechanics, 27(1): 63-72. doi: 10.12090/j.issn.1006-6616.2021.27.01.007
    FENG Chengjun, QI Bangshen, ZHANG Peng, SUN Dongsheng, MENG Jing, NIU Linlin, WANG Miaomiao, TAN Chengxuan, CHEN Qunce. 2018: CRUSTAL STRESS FIELD AND ITS TECTONIC SIGNIFICANCE NEAR THE LONGMENSHAN FAULT BELT, AFTER THE WENCHUAN Ms8.0 EARTHQUAKE. Journal of Geomechanics, 24(4): 439-451. doi: 10.12090/j.issn.1006-6616.2018.24.04.046
    QIN Xianghui, CHEN Qunce, MENG Wen, TAN Chengxuan, ZHANG Chongyuan, FENG Chengjun. 2018: EVALUATING MEASURED IN-SITU STRESS STATE CHANGES ASSOCIATED WITH EARTHQUAKES AND ITS IMPLICATIONS: A CASE STUDY IN THE LONGMENSHAN FAULT ZONE. Journal of Geomechanics, 24(3): 309-320. doi: 10.12090/j.issn.1006-6616.2018.24.03.033
    MA Qiang, DELEQIATI·Jantay, CHEN Yajun, MA Le, YU Jingwei, QI Liqi. 2018: A STUDY ON FAULT STRESS CHARACTERISTICS AND THE LAW OF OIL AND GAS CONTROLLING BY FAULT DEFLECTING: A CASE STUDY OF THE NANPU DEPRESSION. Journal of Geomechanics, 24(5): 587-597. doi: 10.12090/j.issn.1006-6616.2018.24.05.059
    SONG Sui-hong, HOU Jia-gen, LIU Yu-min, CAO Si-fan, HU Chen-bin. 2017: A STUDY ON THE DISTRIBUTION PATTERN OF SUPRASALT FAULT BASED ON NUMERICAL SIMULATION. Journal of Geomechanics, 23(3): 429-435.
    WANG Shang-yan. 2016: FAULTS AND EARTHQUAKES IN GUIZHOU. Journal of Geomechanics, 22(3): 430-438.
    HAN Shu-qin, DENG Jun, YANG Li-qiang, WANG Jian-ping. 2006: APPLICATION OF THE 3D STRUCTURAL STRESS FIELD ANALYSIS IN THE DEEP GOLD EXPLORATION OF THE CENTRAL SEGMENT OF THE ZHAOYUAN-PINGDU FAULT ZONE, JIAODONG AREA. Journal of Geomechanics, 12(3): 338-344,316.
    XIE Xin-sheng. 2003: VARIATION OF CONJUGATE RUPTURES AND LIMITING PRINCIPAL STRESSES WITH THE CRUSTAL DEPTHS IN THE YANYUAN-NINGLANG EARTHQUAKE. Journal of Geomechanics, 9(2): 161-170.
    ZHANG Yue-qiao, P. Vergely, J. L. Mercier. 1999: PLIOCENE-QUATERNARY FAULTING PATTERN AND LEFT-SLIP PROPAGATION TECTONICS IN NORTH CHINA. Journal of Geomechanics, 5(3): 40-46.
    Shen Shumin, Feng Xiangyang, Liu Wenying. 1995: TECTONIC STRESS DRIVING AND MIGRATION OF GAS AND OIL IN NORTHERN TARIM BASIN. Journal of Geomechanics, 1(2): 11-20.
  • Created with Highcharts 5.0.7Amount of accessChart context menuAbstract Views, HTML Views, PDF Downloads StatisticsAbstract ViewsHTML ViewsPDF Downloads2024-052024-062024-072024-082024-092024-102024-112024-122025-012025-022025-032025-040510152025
    Created with Highcharts 5.0.7Chart context menuAccess Class DistributionFULLTEXT: 35.6 %FULLTEXT: 35.6 %META: 62.2 %META: 62.2 %PDF: 2.2 %PDF: 2.2 %FULLTEXTMETAPDF
    Created with Highcharts 5.0.7Chart context menuAccess Area Distribution其他: 5.1 %其他: 5.1 %United States: 0.2 %United States: 0.2 %北京: 1.4 %北京: 1.4 %南京: 0.2 %南京: 0.2 %哥伦布: 0.2 %哥伦布: 0.2 %天津: 0.4 %天津: 0.4 %广州: 0.2 %广州: 0.2 %张家口: 6.3 %张家口: 6.3 %扬州: 0.2 %扬州: 0.2 %杭州: 0.4 %杭州: 0.4 %武汉: 1.0 %武汉: 1.0 %海口: 0.2 %海口: 0.2 %深圳: 0.6 %深圳: 0.6 %漯河: 0.4 %漯河: 0.4 %石家庄: 0.2 %石家庄: 0.2 %芒廷维尤: 34.1 %芒廷维尤: 34.1 %苏州: 0.2 %苏州: 0.2 %西宁: 47.8 %西宁: 47.8 %西安: 0.2 %西安: 0.2 %长沙: 0.4 %长沙: 0.4 %阿什本: 0.2 %阿什本: 0.2 %其他United States北京南京哥伦布天津广州张家口扬州杭州武汉海口深圳漯河石家庄芒廷维尤苏州西宁西安长沙阿什本

Catalog

    Article Metrics

    Article views (183) PDF downloads(7) Cited by()
    Proportional views
    Related

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return