深部不同深度岩石脆延转化力学行为研究

谢和平,高明忠，付成行，鲁义强,6，杨明庆，胡建军，杨本高

(1.深圳大学 广东省深地科学与地热能开发利用重点实验室，广东 深圳 518060; 2.深圳大学 深地科学与绿色能源研究院，广东 深圳 518060; 3.深圳大学 土木与交通工程学院，广东 深圳 518060; 4.深圳大学 深圳市深部工程科学与绿色能源重点实验室，广东 深圳 518060; 5.四川大学 水利水电学院，四川 成都 610065; 6.中铁工程装备集团有限公司，河南 郑州 450000)

向地球深部要空间、要资源是大势所趋，地球深部资源的开发与利用已成为世界各国争先探索的科学制高点，其共性科学基础即是深部岩石力学理论[1-6]。随着能源资源开发由浅部走向深部，深部工程中的“高地应力、高温度、高渗透压”的“三高”地质环境导致岩石的组织结构、基本力学特征和工程响应发生变化，均是深部工程灾害频发且不同于浅部灾害形式的主要原因[7-17]。不同赋存深度岩石力学行为的差异性主要体现在岩石变形、强度特征、破坏特征和脆延性转化等方面[18-21]，其中，岩石的脆延性随深度的演化规律是深部岩石力学研究的重要内容之一。

由于不同深度岩芯难以获取，大多学者仍通过改变围压来近似模拟不同深度的地应力赋存环境，进而分析岩石脆延性状态转化方面的规律，并发现随着围压的增大岩石具有从脆性向延性转化的趋势[22-24]。笔者首创深部岩石原位保真取芯技术，旨在探索研究原位环境下深部岩石的物理力学差异性特征新规律、构建深部原位岩石力学新理论体系[25-26]。VON KARMAN用大理岩进行不同围压条件下的力学实验[27]，发现在围压约69 MPa下，大理岩发生了较大的塑性变形，轴向应变的值达到8%左右时试样还未发生宏观破坏。PATERSON[28]在室温下对Wombegan大理岩做了试验，发现大理岩的围压超过20 MPa后峰值应变显著增加，并且随着围压的增大岩石的力学性质有脆性向延性转变的特性，PATERSON 把这种由低峰值应变下的宏观破裂到高峰值应变下的没有明显破坏迹象的转变称脆延转变。MARTIN和STIMPSON[29]基于大量实验，总结出了一个基于应力和强度关系的脆延转化条件。HANDIN[30]对盐岩在室温下的三轴压缩试验中观测到围压小于20 MPa时会发生脆性向延性转变。MOGI[31]采用大理岩进行实验得出了类似的结果。CLEARY[32]认为岩石的破坏过程在浅部表现为脆性或断裂韧性控制，在深部表现为准静态或延性状态。周宏伟等[33]从微细观角度讨论了岩石内部微结构在脆延转变过程中的作用，认为外部因素为温度和压力，内部因素为内部晶体的微结构，深部岩石脆延转变受到这两类因素共同作用，且它们相互影响。尽管众多学者对深部岩石的脆延转换问题做了大量研究，但大多仍采用通过改变围压来近似模拟不同深度的地应力赋存环境，且多数还集中在定性的讨论阶段[34-35]。

对岩石这种多孔介质来讲，孔隙度对岩石的脆延性转变也有关键控制作用。脆性和延性之间的区别是宏观的，这取决于岩石是否能够承受较大的永久应变而没有宏观的断裂[36]。与岩石的宏观脆性断裂相关的变形机制主要包括晶粒间的微裂纹扩展和摩擦滑动，岩石在延性状态下的非弹性变形机制主要包括晶体可塑性、晶体的转移和分解[36-38]。学者普遍认为，多孔岩石脆延性的转变与初始孔隙度密切相关[36]。而BRACE[39]认为当岩石的孔隙度<5%时，可将其视为少孔岩石，反之则可视为多孔岩石，并将5%这一临界值作为区分少孔岩石和多孔岩石的指标，称为截止孔隙度。基于此，BAUD等[40]研究了Solnhofen石灰岩脆延性的变化，将截止孔隙度临界值降低至3%。BAUD等发现以往关于多孔砂岩脆延性转变的概念过于简单化，岩石的脆延性转变过程可能存在一种中间状态，其破坏模式与压实带的局部发展有关[41-44]。

以上学者主要以1 000 m以浅的岩石为研究对象，很少考虑深度变化导致的岩石性质的改变，以不同围压模拟深度的变化，导致利用浅部岩石获得的规律不能适用于深部岩石。HERRMANN等[45]研究发现露头页岩和深部页岩的力学性质差异较大，这也验证了采用浅部岩石通过改变围压去模拟深部岩石的力学性质会与深部岩石的真实力学性质存在较大的差异。因此，研究不同深度岩石力学性质的变化，有必要同时考虑不同深度岩石本身的差异性和不同深度环境应力水平的差异性。笔者以松辽盆地1 000～6 400 m深度岩石为研究对象，充分考虑岩石原位应力环境以及初始物性差异，初步探究了不同深度岩石脆延转化力学行为特征。

1 基于松科二井原位岩芯的不同深度岩石三轴力学实验

1.1 松科二井不同深度试样制备

松辽盆地作为中国陆相沉积地层最完整的区域，是研究不同深度岩石物理力学行为差异的最佳靶点。特别是松辽盆地有中国最深(完井深度7 018 m)的科探钻井——松科二井，可以获得用以科学探索的深部岩芯，本研究所取岩芯主要集中在4 800,5 100,5 600,6 400 m深度范围。为了系统的研究松辽盆地不同赋存深度岩石的物理力学规律，建立大深度范围的不同赋存深度规律研究，课题组在大庆油田取芯井(龙17、龙10-08、葡深1、台602、茂61-检89、大424)获得的1 000～3 500 m深部范围内的岩芯，综合松科二井样品开展不同深度岩石的物理力学试验。大庆油田取芯井与松科科探钻井的相对位置如图1所示[46]，两者均位于松辽盆地，且地质条件相似，故在一定程度上可将两者岩芯统一进行研究。

实验采用10个深度的岩芯，赋存深度分别为1 000，1 300，1 600，1 850，2 600，3 500，4 800，5 100，5 600，6 400 m。其中1 000～4 800 m为砂岩，5 100～5 600 m为砾岩，6 400 m为安山岩。将取得的岩芯制作成直径为(25±1) mm，高径比为(2±0.2)的圆柱体试样，试样两端面不平行度应不大于0.01 mm;试样上下端直径偏差不应大于0.1 mm;试件表面应光滑，避免因不规则表面而产生的应力集中现象导致试验数据失真。

1.2 实验设备及试验方案设计

本文考虑深度影响的岩石常规三轴力学实验中，1 000～3 500 m深度实验在四川大学MTS815岩石力学实验机上开展;4 800～6 400 m深度在重庆大学GCTS高温高压岩石力学实验机上开展。实验前采用标准铝块和完整砂岩验证了2台实验机在测试中的误差，确保实验机在测试中产生的误差可以忽略。

为深入探索不同深度岩石脆延转化力学行为规律，共设计3个试验方案:① 同一深度(Same depth)岩石不同应力水平(Different stress)的常规三轴实验，简称SD;② 不同深度(Different depth)岩石同一应力水平(Same stress)的常规三轴实验，简称DS;③ 不同深度(Different depth)岩石真实应力水平(True stress)的常规三轴实验，简称DT。实验参数见表1。

2 不同深度岩石脆延特性评价方法

为了深入研究不同深度岩石的力学特性，充分分析松辽盆地不同深度岩石三轴压缩下脆延性特征差异，采用MOGI[47]和RYBACKI E等[48-49]提出的定义典型应力-应变曲线几个特征力学参数的方法进行表征，如图2所示。

图2中A点为应力应变曲线的屈服点;B点为峰值应力点(试样破坏点);OD代表试样的延展性(延性)，定义为材料在无断裂情况下经受最大永久变形的能力。目前关于延性参数度量方法没有统一，本文给出了一个简单的确定方法，从应力应变曲线的起始点O，做直线(斜率为弹性模量E)，在该直线与应力应变曲线峰值水平线交点处向下做垂线，垂线与应力应变曲线的相交点定义为应力应变曲线的屈服点A。屈服点之前的OA段岩石的变形为弹性变形(包括线弹性和非线弹性)。A点之后开始有塑性变形产生，岩石进入屈服阶段。然后定义岩石应力应变关系的几个参数:E代表岩石的弹性模量，从初始线性部分的曲线斜率中获得;泊松比μ由弹性阶段的环向应变与轴向应变的比值获得;H代表应变硬化模量，由后屈服区曲线的斜率获得。εini,εel,εinel,εmax,εfail分别为剪胀开始时的应变，最大弹性应变(屈服应变)，非弹性应变(延性应变)，峰值应力应变和破坏应变。根据图2中定义的几个参数的计算如式(1)所示:

图1 松辽盆地取芯位置[46]Fig.1 Location of cores obtained[46]

表1 不同深度试验样品信息Table 1 Sample information at different depths

图2 三轴压缩实验中岩石的典型应力-应变曲线[47-49]Fig.2 Typical stress-strain curve of rock in triaxial compression test[47-49]

(1)

式中，ΔσL为应力应变曲线直线段两点之间的差值，MPa;ΔεL为直线段两点对应的轴向应变;εd为侧向应变;εL为轴向应变;σel和σmax分别为最大弹性应力和峰值应力;εel和εmax分别最大弹性应变和峰值应力应变。

根据不同深度岩石不同阶段应力应变关系的变化，定义了几个可以表征不同深度岩石脆性特征的参数[48-52]，即

其中，BpreHarden，BiniStrain和BpreStrain均为基于应力-应变曲线特征参数的岩石峰前脆性指标。其中,BpreHarden通过应变硬化模量与弹性模量的比值定义，综合考虑了峰值应力、峰值应力应变、最大弹性应力和最大弹性应变，而BiniStrain和BpreStrain则分别考虑的是岩石发生剪胀时的应变、峰值应力应变和最大弹性应变、峰值应力应变。

不同深度岩石在常规三轴压缩实验中的力学特性的差异不仅存在于弹性阶段，非弹性段和峰后特性的差别也非常显著。为了定量研究岩石在变形过程中的应力-应变关系，考虑不同深度岩石常规三轴压缩试验中的峰后特征，定义出能够表示三轴条件下典型的应力-应变曲线几个特征力学参数的量，如图2所示。其中，峰前特征点的定义与上文相同。这里主要介绍峰后跌落模量M(Descendant modulus)，由应力下降的斜率获得，如式(3)所示。

(3)

式中,εres为残余应变；σres为相应的偏应力值。

从不同深度岩石的应力应变关系中可以看出，在常规三轴压缩实验中，不同深度岩石的脆性特征，特别是峰后脆性特征存在较大的差异。本文在定义峰前脆性特征的基础上，根据图2中的典型应力应变曲线定义能够体现不同深度岩石峰后特性差异的峰后脆性特征指标[49-53]，即

(4)

其中，BresStress，BresStrain和BresER均为基于应力-应变曲线特征参数的岩石峰后脆性指标。其中,BresStress通过应力参数定义，考虑岩石的峰后残余强度和破坏应力;BresStrain通过应变参数定义，考虑岩石的破坏应变和残余应变;BresER通过弹性模量和峰后跌落模量比值定义，其本质是耗散弹性能与可恢复弹性能间的比值，综合考虑了岩石的峰值应力、残余强度、峰值应力应变、残余应变、最大弹性应力和最大弹性应变。BresER脆性评价指标在[-1,0]岩石表现为脆性，BresER为+时岩石表现为完全塑性。岩石峰后脆性变化趋势如图3所示，其展示的是岩石由绝对脆性向延性转换的过程:当跌落模量M<0，且其绝对值大于弹性模量E时，BresER在[-1,0]，岩石表现出很强的脆性特征，其应力应变曲线峰后回弹，如图3中Class Ⅱ所示，岩石峰后无法自稳;若E不变，M值由负无穷大转为正无穷大(应力应变曲线峰后回弹逐渐减弱)后逐渐变小，BresER的值则会逐渐增大，岩石逐渐向延性转换，当岩石表现出延性特征时其应力-应变曲线峰后无明显的强度跌落。

图3 岩石峰后脆性变化趋势示意[52]Fig.3 Schematic diagram of the trend of brittleness after the peak[52]

3 同一深度岩石不同应力水平下脆延差异性特征分析

赋存深度的变化会导致初始应力水平的改变，而岩石所处的初始应力水平对其力学性质具有显著的影响。岩石材料在较低初始应力水平下往往呈现出脆性破坏特征，而在高应力水平下表现出逐渐向延性过渡的破坏特征，其强度随应力水平变化规律也会表现出非线性[24,28-31]。因此本节从不同深度岩石所处的应力环境入手，研究围压对单一深度岩石力学性质的影响，揭示不同初始应力水平对岩石弹性模量、泊松比、强度、脆延性转变等方面的影响。实验以1 600 m深度的砂岩为研究对象，其不同围压下的实际破坏形态如图4所示，岩石受压主破裂面示意如图5所示，偏应力-应变曲线如图6所示。

图4 1 600 m深度岩石在不同围压下的破坏形态Fig.4 Failure mode of rock at 1 600 m depth under different confining pressure

图5 岩石受压主破裂面示意Fig.5 Schematic diagram of the main fracture surface of rock under pressure

图6 不同围压下1 600 m深度砂岩的偏应力-应变曲线Fig.6 Deviatoric stress-strain curves of sandstone at a depth of 1 600 m under different confining pressures

由图4,5可知，1 600 m深砂岩在三轴应力作用下其破坏形态随着围压的增大经历了低围压下的单剪破坏、较高围压下的双剪破坏以及高围压下的无明显裂纹面的涨鼓破坏3个阶段。低围压下的单剪破坏:围压在13.74～34.14 MPa时，1 600 m深砂岩的破坏形式为典型的单剪破坏。当试样达到屈服应力后，试样的承载能力会较快达到峰值，并开始下降。在围压较小的三轴应力作用下，由于受到侧向应力对裂纹张拉发展的抑制，岩石的破坏过程沿着主破裂面法线方向发生偏转，剪切面法线方向并不总沿着最小主应力方向。试样微观上以沿晶破坏为主，穿晶破坏为辅，晶体没有产生大量的塑性变形从而形成单剪破坏面。较高围压下表现为典型双剪破坏:当围压升高至52.14～73.74 MPa时，岩石破坏形式表现为双剪破坏。产生这种现象的原因是在该围压范围内，试样达到屈服点后其承载能力不会很快达到峰值，会产生一定的塑性变形，之后试样承载能力下降也较为缓慢。由于围压的增大，试样的破坏过程产生大量的穿晶破坏，并且部分晶粒发生错位滑移和塑性变形进而产生双剪破坏。高围压下则在宏观上无明显裂纹面的涨鼓破坏:岩石受高围压作用时，特别是围压大于其单轴抗压强度时，沿最大能量释放率发生方向的弹性应变能小于剪切破坏能量，进而发生膨胀破坏。试样达到屈服点之后进一步产生应变需要更大的应力，产生硬化现象，并在此过程中产生大量的永久变形。试样在微观上以晶粒的塑性变形和错位滑移为主，导致试样产生中部涨鼓的延性破坏特征[53]。

由图6中的偏应力-应变曲线可以直接的观察到:破坏前岩石的峰值应变随着围压的增大而增大;另外，随着围压的增大，岩石的塑性也不断的增大，具有脆性逐渐向延性转化的趋势。当围压在13.74～34.14 MPa时，1 600 m深度砂岩偏应力-应变曲线经历了标准的裂隙压密阶段、线弹性变形阶段、裂纹起裂和稳定扩展阶段、裂隙非稳定发展阶段、残余强度阶段5个阶段。并且随着围压的增大，岩石峰后的塑性越来越强，但是试样整体上有一个峰值强度，破坏后承载迅速下降，表现出一定的脆性特征;当围压升高至52.14～73.74 MPa时，岩石显示出由脆性至延性的转化过渡状态;当围压升高至104.94 MPa时，岩石呈现塑性流动状态，表现出延性特征;当围压升高至143.34 MPa时，试样的偏应力随着围压的增大稳定增长，出现应变硬化的现象。上述结果表明对于1 600 m深度的砂岩，岩石脆延转化并非一个定值，而是存在一个脆延性逐渐转换的应力区间。对于1 600 m深度砂岩，围压在50～70 MPa为其脆延性转化区间，并且在围压为13.74～143.34 MPa范围内，随着围压的增大，岩石的性质出现了从脆性—延性—应变硬化的转变。

根据松辽盆地1 600 m深度砂岩的应力-应变关系，采用前文提出的不同深度岩石脆延特性评价计算方法，得到其在不同初始应力水平下的典型力学参数，如图7所示。在这里需要注意的是，试样在104.94～143.34 MPa围压下表现出很强的延性特征，所以没有残余强度σres和相对应的应变εres，因此在计算峰后跌落模量M时，M=0，BresER=+∞。

由图7可知，1 600 m深度砂岩的弹性模量随着初始应力水平的增大略有提高，但整体上升高幅度不大(图7(a))，围压增大对弹性模量的影响有限。但1 600 m深度砂岩的泊松比随着初始应力水平的增大而呈现增大趋势(图7(b))，表明随着初始应力水平的增大，岩石的塑性逐渐增大。1 600 m深度砂岩的三轴相对强度随着初始应力水平的变化呈现明显的对数增长趋势(图7(c))，也即随着围压的升高，三轴相对强度的增长趋势逐渐减弱并趋近于0。残余强度随着深度增大也逐渐增大，围压超过100 MPa后，1 600 m深度岩石的残余强度与峰值强度基本相等(图7(c))。应变硬化模量随着围压的增大规律性不强(图7(d))。但峰后跌落模量随着深度的增大呈现急剧减小趋势。围压超过100 MPa后峰后跌落模量的值趋近于0，岩石峰后表现出塑性流动特征(图7(e))。

为了研究围压变化对岩石的弹塑性特征的影响，将1 600 m深度砂岩不同围压下的几种脆性指标的变化趋势绘于图8中。

图7 1 600 m深度砂岩在不同围压下的三轴力学参数Fig.7 Triaxial mechanical parameters of sandstone at a depth of 1 600 m under different confining pressures

图8 1 600 m深度岩石在不同围压下脆性评价指标随围压的变化规律Fig.8 Variation of brittleness evaluation index with confining pressure for 1 600 m deep rock under different confining pressures

由图8(a)可知，1 600 m深度砂岩的峰前脆性指标总体上都是随着围压的增大的而减小。由图8(a)可以看出，随着围压的增大，岩石的非弹性段所占的比例逐渐增大;由图8(b),(c)可以看出1 600 m深砂岩峰后脆性指标BresStress和BresStrain随着围压的增大而减小，BresER随着围压的增大而增大，这都表明岩石的峰后塑性随围压的增大逐渐增强，直到峰后表现为完全塑性。对于单一深度的岩石，其矿物组分是一致的，这种情况下围压是影响其脆性特征的主导因素，随着围压的增大，岩石的脆性会逐渐减弱。

4 不同深度岩石同一应力水平下脆延差异性特征分析

岩石的力学性质受诸多因素的影响，如岩石物性、组分、围压、温度、含水率等。其中岩石物性对岩石的物理力学行为会产生重要的影响，是起主导作用的关键因素之一。通常岩石在漫长的地质作用下，其组分和矿物成分复杂。分析哪些矿物组分是影响岩石力学参数差异及脆延性转化的主要矿物是研究不同深度岩石力学性质差异性的重要内容之一。本节着重研究不同深度岩石在仅考虑岩石物性差异时其力学行为的变化。

从上文单一深度岩石不同围压下的脆延性分析可以看出，围压对岩石的脆延性转化有着直接的影响，并且1 600 m深度砂岩在50～70 MPa处于脆延转化状态。本节在研究岩石物性对岩石力学性质影响时，采用该临界范围的围压52.14 MPa(对应于2 600 m深度的初始应力水平)，研究不同深度岩石的岩石物性对力学性质的影响，得到偏应力-应变曲线随赋存深度的变化趋势如图9所示。

图9 不同深度岩石在52.14 MPa围压下的偏应力-应变曲线Fig.9 Stress-strain curves of rocks at different depths under 52.14 MPa confining pressure

由图9可以直接的观察到:在同一围压下不同深度岩石的峰值偏应力随着深度的增大呈现增大趋势;另外随着深度的增加岩石的脆性不断增强，延性不断减弱。具体随深度的变化过程为:浅部岩芯(1 000～1 850 m砂岩)在达到峰值荷载之后基本呈现弹塑性流动状态，表现出很强的延性特征;2 600 m砂岩岩芯在保持一段塑性变形之后应力水平才逐渐降低，然后进入稳定变形阶段;深部岩芯(3 500～4 800 m砂岩、5 100～6 400 m砾岩以及6 400 m安山岩)在达到峰值荷载之后应力水平迅速下降，表现出脆性特征。表明在52.14 MPa的相同围压下，不同深度岩石的力学性质受其自身矿物组分的影响，岩石由浅至深表现出延性—脆延转化—脆性的特征，其主要原因为随着深度的增大，中等相和硬相矿物含量的增加导致的。

同一围压(52.14 MPa)下不同深度岩石的三轴力学参数随深度的变化趋势如图10所示。由图10可知，在同一围压(52.14 MPa)下，不同深度砂岩和砾岩的弹性模量均随着深度的增大呈现增大趋势(图10(a));泊松比随着深度的变化离散性较大，但总体上呈现减小趋势(图10(b));不同深度下砂岩的三轴相对强度和残余强度基本上随着深度的增大而增大，而砾岩则相反(图10(c));不同深度岩石的应变硬化模量随着深度的增大呈现增大趋势(图10(d));峰后跌落模量随着深度的增大变化趋势差异较大(图10(e))，在1 000～1 850 m深度范围内的砂岩峰后跌落很小，表明该范围内岩石的峰后处于延性阶段。在2 600 m深砂岩的峰后模量略有增大，表明处于应变软化阶段。3 500～4 800 m深砂岩和5 100～5 600 m深砾岩的峰后模量绝对值急剧增大，表明岩石峰后强度迅速跌落(其中4 800 m深砂岩的跌落模量为负值，表明该深度下岩石脆性很强，岩石峰后曲线回弹，属于图3中的II型曲线)。6 400 m深度下安山岩的峰后模量有所减小。在同一应力水平下，不同深度岩石的脆性特征主要受其矿物组分影响。整体上，随着深度的增大，由于岩石矿物组分的改变，脆性矿物的影响作用逐渐凸显，岩石的脆性特征逐渐增强。

图10 同一围压(52.14 MPa)下在不同深度岩石学参数随深度变化规律Fig.10 Variation of petrological parameters with depth at different depths under the same confining pressure (52.14 MPa)

由图9可以看出，同一围压(52.14 MPa)下不同深度岩石的脆性特征也存在较大的差异，根据前文定义的脆性指标，同一围压(52.14 MPa)下不同深度岩石的脆性特征随深度的变化趋势如图11所示。

图11 同一围压(52.14 MPa)下不同深度岩石的脆性特征随深度的变化趋势Fig.11 Changing trend of brittleness characteristics of rocks at different depths with depth under the same confining pressure(52.14 MPa)

鲁义强等[46,54]把岩石矿物成分分成硬相矿物、中等相矿物和弱相矿物(硬相矿物主要包括石英、长石和黄铁矿;中等相矿物主要包括碳酸盐类矿物，主要是方解石;弱相矿物主要包括云母、黏土矿物等)，并以松辽盆地1 000～6 400 m深度10个深度的岩石为研究对象研究了其力学行为与硬相矿物、中等相矿物和弱相矿物间的关系，研究结果表明:松辽盆地岩石随着深度的增大，基本上呈现弱相矿物减少，中等相矿物增多，硬相矿物略有增大的趋势;不同深度岩石的脆性随着弱相矿物含量的增加而逐渐减小;不同深度岩石的脆性与中等相矿物和硬相矿物含量呈正相关关系。而由图11(a)可知，在同一围压(52.14 MPa)下，1 000～4 800 m砂岩峰前脆性指标整体上都是随深度的增大而增大的，表明随着深度的增大，岩石的脆性增大。总体上1 000～6 400 m内岩石的峰前脆性指标随着深度的增加而增大，这表明在同一围压下，矿物组分是影响岩石脆性的主导因素，硬相和中等相矿物占比越大，岩石的脆性特征越强。不同深度岩石的峰后脆性指标有一定的差异，由图11中BresStrain和BresER变化趋势可知，在同一围压下1 000～4 800 m深度砂岩以及5 100～5 600 m深度砾岩的峰后脆性随着深度的增大呈现减小趋势。

5 不同深度岩石真实应力水平下脆延特征初探

根据测试结果[55]，把松辽盆地不同深度岩石典型应力-应变曲线关系绘在同一个坐标系中，如图12所示，基本上岩石的三轴强度和残余强度都随着深度的增大，围压的增加而增加。从图12可以看出，不同深度不同岩性岩石的力学性质并非随着深度的增加，延性逐渐增加，也不是随着深度的增加走向动态脆性破坏，而是随着深度的增加呈现3个阶段的特征:峰后应变软化(1 000～3 500 m深砂岩)、峰后脆性(4 800 m深砂岩和5 100～5 600 m深砾岩)、峰后塑性流动(6 400 m以深的安山岩)。其主要原因是岩石的力学特性受深度影响是综合作用的结果，其中岩石的组分结构和原位应力状态是2个主要因素。

图12 不同深度岩石真实应力下偏应力-应变曲线[55]Fig.12 Stress-strain curves diagram of rock at different depths under real stress[55]

由图12可以看出，初始应力的变化主要对不同深度岩石常规三轴强度及峰后特性产生影响。从常规三轴压缩试验结果可以看出，随着深度的增加，初始应力水平增大，偏应力-应变曲线的斜率会明显增大。由此可见，只有认真的考虑初始应力水平的影响，才能充分认识不同深度岩石的力学性质差异性规律。不同深度岩石在相应深度地应力条件下的常规三轴力学参数随深度的变化趋势如图13所示。

图13 松辽盆地不同深度岩石常规三轴力学参数的差异性[55]Fig.13 Differences inconventional triaxial mechanical parameters of rocks at different depths in the Songliao Basin[55]

由图13可知，在常规三轴压缩试验中，不同深度岩石的力学参数随着深度的变化基本都呈现非线性趋势。不同深度的砂岩和砾岩的弹性模量E均随着深度的增加而增大(图13(a))，显然围压的变化对岩石的弹性变形产生了一定的影响;泊松比μ随着深度的增大呈现增大趋势(图13(b))，这是由于围压的影响抑制了岩石的横向变形，并且促进岩石向塑性转化，从而导致泊松比的增大;三轴相对强度和残余强度均随着深度的增长而增大(图13(c));应变硬化模量总体上随着深度的增加而增加(图13(d));峰后跌落模量的非线性更强，对于1 000～4 800 m深度砂岩，其峰后跌落模量总体上随深度的增加逐渐增大，5 100～5 600 m砾岩以及6 400 m深安山岩的跌落模量随深度增加逐渐减小(图13(e))。

根据式(2)和(4)中的脆性指标计算方法，深部不同深度岩石常规三轴压缩试验的脆性特征随深度的变化如图14所示。

由图14(a)可知，不同深度岩石的3类峰前脆性指标随深度的变化趋势同步性较高。1 000～4 800 m深度范围内砂岩3类峰前脆性指标总体上随着深度增加呈现先减后增的趋势，3 500 m深处砂岩塑性最强，其原因推测为围压和矿物组分对砂岩脆延特性影响存在博弈现象。1 000～3 500 m围压对砂岩的脆延特性起主导作用，故随着围压的增大该范围内砂岩的脆性逐渐降低。3 500～4 800 m深砂岩矿物组分对砂岩的脆延特性起主导作用，故随着深度的增加，长石、碳酸盐类矿物等硬相或中等相矿物含量占比增加，岩石的脆性特征增强;5 100～5 600 m砾岩的峰前脆性指标随深度的增加而减小，其脆性特征减弱。由图14(b)可知，不同深度岩石的峰后脆性指标随深度的变化趋势各异。1 000～4 800 m砂岩的BresStress波动较大，规律性不明显。5 100～5 600 m砾岩的BresStress随深度的增加而减小，6 400 m安山岩的BresStress最小，表现出峰后塑性流动特征。

6 考虑应力环境和赋存深度影响的岩石脆延特性对比

为更加直观的分析不同深度应力环境和矿物组分对岩石脆延特征的影响，探索深部岩石脆延性转化机理，将同一深度岩石不同应力水平下、不同深度岩石同一应力水平下以及不同深度岩石真实应力水平下岩石的脆性指标进行对比，如图15所示。

图14 不同深度岩石脆性评价指标随深度的变化规律[56]Fig.14 Variation of the brittleness evaluation index of different depth rocks with depth[56]

图15 3种条件下岩石脆性指标对比Fig.15 Comparison of rock brittleness indexes under three conditions

图15中左纵轴围压与右纵轴深度一一对应。SD为同一深度岩石(1 600 m砂岩)不同应力水平下脆性指标的演化曲线，对应围压纵轴;DS为不同深度岩石同一应力水平下(52.14 MPa)脆性指标的演化曲线，对应深度纵轴;DT为不同深度岩石真实应力水平下脆性指标的演化曲线，与围压纵轴和深度纵轴均对应。由图15可以看出，同一深度砂岩在不同应力水平下，随着围压的增大，其脆性整体上逐渐减小，在这种情况下围压对岩石的脆性起主导性作用，会抑制岩石的脆性;不同深度岩石在同一应力水平下，随深度的增加其脆性矿物占比增加，导致岩石脆性逐渐增大，在这种情况下岩石的矿物组分对其脆性起主导作用;不同深度岩石真实应力水平下，其脆性呈现先增大后减小趋势。在这种情况下围压与矿物组分均对岩石脆延特征有影响，且2者之间存在博弈现象。上述结果表明岩石的脆性受到多种因素的综合影响，其中应力环境和矿物组分是2个关键因素。通过对同一深度岩石不同应力水平，不同深度岩石同一应力水平，不同深度岩石真实应力水平下的三轴力学实验的研究，可以得出一个基本的结论:岩石矿物组分和初始应力水平都是影响不同深度岩石力学性质的重要因素。随着深度的增加，长石、碳酸盐类矿物等脆性矿物含量的增加会导致岩石的脆性增长，但是随着深度增加，初始应力水平的增加又会抑制岩石的脆性增长。这2种控制因素此消彼长，存在博弈现象，具体在哪种条件下，哪个因素起主导作用，还需要大量的数据验证。同时，在岩石脆延性转化过程中，能量转换机制的作用也不容忽视[57-59]。目前已经有从很多能量角度表达脆性指数的指标，下一步可以结合不同深度应力条件，考虑从能量存储-耗散转换的角度，提出能科学表达岩石脆延转化的新指标。

7 结 论

(1)深部岩石脆延转化并非瞬时的，而是存在一个脆延性逐渐转换的应力区间。对于松辽盆地1 600 m赋存深度的砂岩，实验表明50～70 MPa围压为其脆延性转化区间。

(2)同一深度砂岩不同应力水平下其脆性特征主要受围压影响。随着围压的增大，砂岩性质出现了从脆性—延性—应变硬化的转变，其峰后塑性逐渐增强，直到峰后呈现完全塑性。

(3)相同围压下岩石的脆延性特征主要受其本身矿物组分的影响。对于松辽盆地岩石样品，其硬相和中等相矿物含量随深度的增加而增大，导致其脆性增大，由浅至深表现出延性—脆延转化—脆性的特征。

(4)不同深度岩石真实应力水平下其峰后特性随着深度的增加呈现出不同的特征:1 000～3 500 m深砂岩表现出峰后应变软化特征、4 800 m深砂岩和5 100～5 600 m深砾岩表现出峰后脆性特征、6 400 m以深的安山岩表现出峰后塑性流动特征。

(5)深部岩石的脆性特征受诸多因素影响，其中岩石的矿物组分和应力环境对岩石的脆性影响存在博弈现象，中等相、硬相矿物含量的增加会导致其脆性增长，而加载应力水平的增加，又会抑制其脆性增长。

需要特别说明的是本文采用的不同深度岩石均从松辽盆地不同深度获取，旨在更加真实的研究不同深度岩石脆延转化力学行为规律。然而由于不同深度岩石并非同一岩性，故研究结果存在一定缺陷，后续还需结合围压和矿物组分对岩石脆延特性影响的单一变量实验，进一步深入探索不同深度下岩石的脆延特性转换机理。同时，由于现有技术与理论的限制，本文并没有考虑深部岩石取芯过程中的主应力转化、加卸载条件变化等力学过程，以及温度、湿度等因素对岩石力学特性的影响。针对这些影响因素后续还需深入探索研究，这也是未来深部原位岩石力学需要重点探索的方向。