第3章 原岩应力及其量测

《第3章 原岩应力及其量测》由会员分享，可在线阅读，更多相关《第3章 原岩应力及其量测（22页珍藏版）》请在装配图网上搜索。

1、第三章原岩应力及其量测3.1地球及其构造的一般概念地球的绝对年龄估计在5055亿年。在4547亿年以前开始形成地壳，就 是说地球诞生在47亿年以前。整个太阳系也是在不到50亿年前由尘埃和大气形 成。我们目前所熟知的地球，具有适于人类生存的大气和丰富的资源，这颗行星 的内部仍在活动。这点已由地震、火山、张开和闭合的大洋及漂移开来的大陆所 证实。地壳400km1000km地幔3000km下地幔外地核5000km6000km内地核地球内部结构示意图根据对深部地带进行地震研究 而得到的现代概念，地球可分为地 壳、上地幔、下地幔、外地核和内地 核。地壳的平均厚度为32km，而且 在大陆上的变化范围是20

2、70km， 在海洋中其变化为515km。地壳是 以莫霍面为分界面，是1909年由南 斯拉夫的莫霍洛维奇契首先发现了 M面。在该面以下，弹性纵波的速度 v突然增长，达到8km/s，而在地壳 中通常是67km/s （最大值为 7.4km/s）。上部地幔物质密度：33 37kN/m3 ；地壳物质密度：27郭里采层（导电 莫霍面（弹性纵波波速从地壳的 率和地震波速迅67km/s迅速增加到8km/s ）速增加的分）km 层）30kN/m3。在地壳范围内，可按地震波特征分为三个主要分层:弹性纵波速度vp =2.05.0km/s，厚度1015kmvp =5.56.0km/s，最大厚度 3040kmv p =

3、6.57.4km/s，其厚度为 10 20km它是两个分层之间弹性波速度变化的地震分界面现在，采矿工作主要是在小于10001800m的深度内进行。在欧洲，有些矿 井的开采深度约达2000m ；在南非及印度，个别金属矿井的开采深度已超过 30003500m。开采石油和天然气的深度达到60007000m。最深的构造钻孔和 勘探钻孔已超过12000m，并开始实现钻孔深度达15000m的计划。上述数字提供了有关地球开发深度的概念及其人类当今已经直接达到和可 能近期达到的深度。显然这些深度属于地壳上部的范围内，其厚度与地球直径相 比微不足道。然而浅部地壳的组成结构及其应力状态是矿山岩石力学和矿压理论 关

4、注的重点问题之一。3.2原岩应力天然状态下地壳中存在地应力，通常在地学中称之为地应力。其主要包括由 岩体重量引起的自重应力和地质构造作用引起的构造应力等。地应力这个概念是 由瑞士地质学者Haim在19051912年间首次提出来的。地应力是在历史地质 作用下发展变化而形成的。它与岩体自重、构造、运动、地下水及温差等有关， 同时又是随时间、空间变化的应力场。但在工程年代，应力场受这种地质作用时 间的影响可以忽略。在采矿工程中，把这种未受采掘扰动影响的岩体原始应力， 又称为原岩应力。在井巷和采场等地下工程结构稳定性分析中，原岩应力是一种初始的应力边 界条件，同时原岩应力是引起地下工程结构变形和破坏的

5、力源。采矿工程中，地下采掘空间对周围岩体内的原岩应力场产生扰动，使得原岩 应力重新分布，并且在井巷和采场的围岩中产生几倍于原岩应力的高值应力（所 谓的二次应力）。围岩随之变形，随着时间的延长，围岩变形继续扩大，甚至引 起围岩破坏或支护物破坏，这就是我们常说的矿山压力显现。由此可见，矿山压 力的来源与原岩应力密切相关，围岩稳定性显然是以原岩应力场为前提条件的。 在计算任何人工开挖的岩体周围的应力分布以前，必须测量或估算开挖前的应力 状态。3.2.1地壳浅部原岩应力实测结果地壳内部的原岩应力场是一个颇为复杂的问题，人们获得原岩应力状态的途 径，主要是通过现场实测来实现。虽然各个国家和地区对原岩应力

6、测量做了大量 工作。但是关于完整应力状态的资料却获得很少，且测量深度也都在3000m之 内，故属地壳浅部。（1 ）原岩应力随深度变化1953年瑞典H.Hast在斯堪的纳维亚半岛首先进行了原岩应力实测工作。此后，欧、美、澳大利亚和我国都先后开展了大规模原岩应力实测工作。E.T.Brown 和E.Hoek（1978）研究了遍及世界不同地区的原岩应力测量，并进行了汇总。在进行资料选择时，对于那些特别反常的地质条件（如近期仍出现构造活动的地 区）的实测结果均略去，只选用了可靠的结果。见下图。上图是铅直应力与深度变化的关系。统计结果表明，铅直应力。与深度的 关系为：。 = 0.027z（MPa） 这是一

7、个重要的铅垂应力估算公式。值得注意的是上式的比例系数与地壳浅部岩石的容重相吻合，通常y= 20 - 30KN/m3。即实测结果说明，铅直应力与上覆岩层的重力相一致。下图是平均水平应力广2 4疽。，）与铅垂应力c 之比K，随埋藏深度Z的变化关系。通过分析发现K值通常取值为:1500Z+ 0.5深度小于500米时，水平应力c人明显大于垂直应力c ；当深度1000米， 水平应力与垂直应力趋于相等，处于静水压力状态。这是因为三个主应力差值很 大时，岩石不可能承受很高应力，否则必然发生破坏，达到新的平衡。3.2.2原岩应力中各应力分量之间的比较(1)平均水平应力.与垂直应力乙的比较。从上面两个图的统计结

8、果看，一般情况下， Z相当于上覆岩层的自重，而 水平应力的波动范围就比较大。且一般大于铅垂应力，其产生原因。一般归结为 地壳的构造运动。据国内外实测资料统计，平均水平应力与 z的比值大部 分在0.81.5之间。见下表统计结果。平均水平应力与铅垂应力的比较表国家名称jhavz叩 zL1.2中国1240282.09澳大利亚022782.95加拿大001002.56美国1841413.29挪威1717665.56瑞典001004.99南非4124352.50前苏联5129204.30其它地区37.537.5251.96(2 )水平应力与间的比较地壳内水平应力中的两个主应力 X与 ,在数值上一般不相等

9、，这一统计结 果反映出了水平应力具有较强的方向性，见下表。水平应力比较实测地区统计数目 /VX1.00.750.750.50.50.250.250(%)斯堪的纳维亚等地5114.067.013.06.0北美22222.046.023.09.0中国3514.345.725.714.3日本3531.048.015.06.0统计值34321.32551.67519.1758.825(3)铅垂应力卢自重应力尸之间的比较岩体上覆岩层的重量是形成岩体初始应力的基本因素之一。一般认为岩体的 铅垂应力大体上相当于上覆岩层的重力P，但并非所有实测结果都如此，从我 国的实测结果表明，铅垂应力 Z与单位面积上的上覆

10、岩层重力P的比例在 0.4319.8之间变化，如果考虑到成果的分散性，以M = 0.81.2作为大体上相Pz等的情况，则仅占8.7%，而 : P 1.2的占69.6%。 这些资料说明，多数的 JPz ；1即铅垂应力多数情况下大于上覆岩体的重量。 这种现象只能解释为某种力场作用的结果。而这种力场不是完全由上覆岩层自重 所引起的。3.2.3自重应力自重应力一由于岩石自重引起的应力称为自重应力。(1) Haim法则(1878年，译为海姆)瑞士地质学家Haim在观察了大型越岭隧道围岩工作状态之后，认为原岩体 铅垂应力为上覆岩体自重。在漫长的地质年代中，由于岩体不能承受较大的差值 应力和与时间有关的变形

11、的影响，使得水平应力与铅垂应力趋于均衡的静水压力 状态。i.e: = = = P由于静水压力下无剪应力，所以任意方向都是主应力方向。0 = = = Px/(2 )金尼克解(苏-A-H-Duhhuk , 1925 )金尼克认为地下岩体为线弹性体，其铅垂应力等于上覆岩体自重：。=p。 在水平方向，岩层内的侧向应力。与。y相等，且水平方向的应变为零： 。=。8 =8 = 0由广义虎克定律：8 = 一* +b I。8 = E y 一 G +b ）L 08 = E y _日 +y 火 0则可解出：y =。= iyi- P令X = -侧向压力系数1 则。有：y =y =Me . = XP一般岩石的泊松比=

12、0.150.35。.X = = 0.180.541 当=0.5时，X = 1，则金尼克公式与Haim法则一致。3.2.4构造应力构造应力是由于地质构造作用引起的应力。地质构造运动（含地震）归根到 底是一个岩层变形与破坏的力学过程，与之对应的应力场叫构造应力场。在构造应力场研究中，我们只能知道构造运动结果（例如地表或基岩的变形 和破裂情况：地震得震源和震级等），而要寻找的是造成这些结果的力源，这是 一个反序的问题。在构造力场求解中，通常无法知道初始应力状态，不易弄清楚 深部构造的情况和深部地质体的力学性能。只能进行模拟或假想研究。下面是 Vening Meinez构造应力场力学模型。Vening

13、Meinez 模型为了分析地壳上部任何一点应力的作用方式，VeningMeinez采用了一种简 便方法。在地球中，采用球体坐标，从地壳上层取一单元体，以地心为原点，设 所取的单元体的六个面均为主平面。由沿B B方向的力平衡条件：（By 一、-cos + 气 +静.8J sin 80SrABCD80A AB B(x ).+ y +一 -8w sin 8wSw8W)=0A AC C80，轲 1， cos80sin 80808甲 8w R - d 80 = 0r 12S = R 80 R 8w= R2808S = R 8v d = R d 8qSAABB = R .80d = R d80代入平衡方程

14、式一一一 （Qb 一一、一一 一 R280.8w+ b+-80 80 R d .8w+ bw +略去高阶无穷小量：b R2 .8008甲+气 R d 80.8w+b R d .808甲=0b R + b0 d +b d = 0br一 = -d 注：ABCD是地球的水平面。气+b R上式说明，平行于水平面的各个应力分量总和的绝对值与垂直方向应力分量 绝对值之比，等于地球半径与受应力作用岩层的深度d之比。如若受构造应力作 用影响的地壳深度为2km的话，地球半径以6000km计算，则垂直应力分量约占 水平应力分量总和的 1/3000。若受构造应力影响的地壳深度为 10km，则 圣一 =1600，从此

15、可以看出：水平应力分量的重要性远远超过垂直应力分量。 气+b3.2.5影响原岩应力状态的因素(1 )地形和地质条件对自重应力的影响o z z地形的起伏影响山体的自重应力分 布，山体内沿着水平面上自重应力的分布 状况和地表形状完全相似。试验和计算结 果表明，岩层的初始应力方向多数微倾斜 于山顶方向，并且在数值上比按最大覆盖 层厚度(山顶到水平面间距离)计算的自 重应力要小得多。 pz(2)裂隙组及不连续面对构造应力的影响地质构造对自重应力的分布也有影 响，通常在褶曲两翼显示应力增大，而在 褶曲中部应力降低。3.2.6我国地应力分布及量测的基本情况中国地应力量测的试验和研究始于60年代。60年代初

16、在地下矿山的巷道、 桐室表明利用扁千斤顶法测量围岩表明的应力状态。1964年，在陈宗基教授的 带领下，中科院武汉岩土所在湖北大冶铁矿进行了国内首次应力解除法测量，测 量深度为-80m。1966年3月，李四光教授在河北上吴县建立了全国第一个地应 力观测站。60年代后期，中科院武汉岩土所、长沙矿冶研究院、地质所等开展 了地应力测量。70年代中后期，地应力测量在水电部门也得到广泛开展，长江 科学院等都开展了这方面工作。但90%以上的地应力测量则分布在地震研究、 水利水电、采矿行业。通过丰富的地应力测量资料，可对我国大陆浅层地壳应力 的分布规律有了初步认识，并且有以下明显分区特点：(1) 华北地区以太

17、行山为界，东西两个区域有较大差异。太行山以东的华 北平原及周边地区，其最大主压应力的方向为近东西向，而太行山以西 最大主压应力方向侧为近南北向；(2) 秦岭构造带以南的华南地区，最大主压应力的方向为北西西至北西向;(3) 东北地区主压应力方向以北东东向为主；(4) 西部地区测得的最大主压应力以北北东为主，个别为近南北向。(5) 在滇西南北构造带上，小江断裂带附近最大主压应力的方向为近东西 向，从此断裂带向西，包括澜沧江断裂以北。鲜水河断裂以南地区，最 大主压应力的方向逐渐转为北西向或北北西向，图3 1、3 2是我国地应力分布的一般特征。其中，图3-1主要反映的是中国大陆板块受到的其它板块挤压作

18、用。图3-2是中国东部地区及日本实测的地应力情况。3.3原岩应力测量3.3.1概述原位测量是目前取得工程需要的不同深度原岩应力可靠资料的唯一方法。因 为尽管原岩应力的各种假说和理论对认识地壳的受力规律有一定的参考价值，但 对于工程而言，都或多或少、或大或小存在各种地质构造和影响原岩应力大小和 方向的错纵复杂的因素，因而没有也不可能有任何一种理论可以完全取代实测方 法而能给出工程需要的可靠的资料。美国、西欧各国、澳大利亚、加拿大、南非、 日本等都普遍开展了原岩应力量测。我国于60年代末进行了这项工作，并于70 年代末开展的比较普遍。1964年，在陈宗基教授带领下，中科院武汉岩土所在 湖北大冶铁矿

19、进行了我国首次应力解除法测量原岩应力，测量深度为-80m。1966 年，李四光教授在河北上吴县建立了我国第一个地应力测量观测站。原岩应力测 量在葛洲坝、二滩水电站、长江三峡建设工程中，发挥了重要作用。许多煤矿尤 其是金属矿都开展了大量的原岩应力测量工作。原岩应力测量是指通过某种测试手段测得岩体某一点上的应力数据，即组成 原岩应力场各个应力分量的大小和方向。应力测量结果对于分析地下工程的稳定 性，进行采矿或地下工程设计以及施工管理都是必不可少的依据之一，因而受到 广泛重视。30年代有人用量测洞壁的应变来计算岩体的初始应力状态。自50年代后， 人们着重研究岩体深部未受扰动的应力状态。相继出现了雷曼

20、的门塞式应变仪、 哈斯特压磁应力计以及三向应变计等。1969年在里斯本召开了 “岩体内应力测 定”国际会议。于1976年和1977年分别召开两次与岩体应力测量有关的国际会 议，促进了岩体应力测量技术的发展。岩体应力测试方法很多，按测试的物理量来分，有直接法（观测应力）和间 接法（由测应变、变形及其它物理量转为求应力）；按量测部位分，有深孔法和 表面法；按测量元件分，有机械式、电阻式、电感式或光学的等等。为了便于了 解，先将它们列成下表。岩体应力测试方法测量方法被测物理量所用仪器、设备用途应力恢复法应变或压力(1) 遍千斤顶；(2) 钢弦或电阻式应变仪(3) 频率仪或应变仪岩体表面应力测量应 力

21、 解 除 法孔底平面应变应变(1) 孔底应变计；(2) 电阻应变仪(1) 已知主应力方向求岩体平面 应力大小和方向；(2) 用三个钻孔汇交测取岩体三 向应力分量孔壁应变应变(1) 。36或。46橡皮叉式 三向应变计；(2) 电阻应变仪或应变 采集系统石体二向应力大小和方向孔径变形变形(1) 36-2型钢环式孔径 变形仪；(2) 应变仪用三孔汇交测求岩体的二向应变分量钻孔应力计法应力(1) 玻璃应力计；(2) 简易光弹仪长期监测岩体内应力变化水压致裂法压力(1) 封堵器；(2) 液压泉等已知一个应力方向的条件下，测 求岩体二向应力物理方法声波法声速声衰减、声 发射(1) SYC3型岩石声波参 数

22、仪；(2) 声发射仪(1) 测量岩体二向应力；(2) 长期观测应力变化；(3) 探测岩体声发射源地震法弹性波(1) 微震仪；(2) 测振仪(1) 测地质构造应力；(2) 测岩体动力特性到目前为止，岩体应力测试都是以弹性理论为依据的，即假定岩体是均质连 续、各向同性的线弹性体。这与实际不符合，有待进一步改进。3.3.2应力解除法原理及现场实施步骤(1)原理o/y+Ayz+Az地壳内有一处于原岩应力状态、。，、c乙的单元体，其尺寸分别是x、y、 z，见上图。若将它与母岩分离，相当于解除单元体上的外力，cx、c,、c乙, 其单元体尺寸分别增大到x + Ax，y +颂，z + Az，如果单元体是弹性的

23、，则可 恢复到受载荷前的尺寸，则恢复应变分别为：AxAyAze x=T，=亍如果通过测试得到8x , 8x , 8 ,又已知岩体的弹模E和泊松比R，依据虎 克定律可计算出应力解除前的应力。从这一过程中可以看出，应力解除法的使用 时需要满足下列条件才能成立。（1）岩体是均质、连续、完全弹性体；（2）加载与卸载时应力与应变之间的关系完全相同；（3）单元体自重忽略不计。（2）应力解除法的现场实施步骤应力解除法是目前测量岩体原岩应力使用最普遍和比较成熟的办法。孔径变 形法、孔壁应变法和空心包体应变法的应力解除过程的现场实施步骤是相同的。具体如下：（人）自地面或地下井巷围岩暴露面，用 钻孔钻进应力解除孔

24、（俗称大孔） 至原岩应力区。该地点应不受到钻 孔以外其它井巷工程的影响，否则 便成了井巷围岩内次生应力场的测 量，解除孔直径以D表示。（8）磨平大孔孔底后，再同心地钻进小 直径钻孔穿过待测点，该孔为测量 用孔，直径以d表示。大、小孔应 满足关系D=（35）d。为了尽量减 小大孔孔底的端部效应影响，测量 截面与大孔孔底距离不小于（23） D。（C）于测量孔内安装测量应变或位移的 传感元件。（D）用大直径（D）钻具对测量孔作取 芯钻进，至超过测量孔孔底以后即 可取断岩芯，实现应力解除。这一 过程中岩芯在原岩应力作用下，经过钻进测量孔受到一次扰动并很快达到 新的平衡状态。传感元件感受到的则是取芯钻进

25、，直至卡断岩芯，即第二 次扰动时的应力场的效应，最终读数即为该次实测的原始数据。3.3.3孔壁应变法孔壁应变法是在钻孔壁上粘贴三向应变计，通过测量应力解除前后的应变， 来推算岩体应力，利用单一钻孔可获得一点的空间应力分量。CSIR三轴孔壁应变计是该方法的主要测量仪器。也是ISRM建议方法中推 荐的测量仪器。它是南非的E.R.Leeman于1966年开发成功并第一次只用单孔且 一次就可测定三维全应力的技术。CSIR三轴孔壁应变计的主体是三个测量活塞，直径约为1.5m的活塞头是由 橡胶类物质制成的，端部为圆弧状，其孤度和钻孔孤度相一致，以便和钻孔保持 紧密接触，在端部表面粘贴4支1.5cm电阻应变

26、片，组成一个相互间隔45。的圆 周应变花，三个活塞也即三组应变花位于同一圆周上，最初的设计是不等间距分 布，夹角分别为兀2 女泌。后来格雷（W.M.Gray）和托乌斯（N.A.Toews）分 析了应变花分布对测量精度的影响，认为三组应变花等间距分布最好，故后来的 设计改为三个活塞成120。等间距分布，其外壳由前后两部分组成。在前外壳端部 有要一圆槽，上贴一支应变片，后外壳端部有连接14根电阻应变片导线的插头 （见图5 11），使用时首先将一个直径约1.2cm，厚0.8cm的岩石圆片胶结在前 壳端部的应变片上，供温度补偿用，然后将三个活塞头涂上胶结剂，用专门工具 将应变计送入钻孔中测点部位，启动

27、风动压力。将活塞推出，使其端部和钻孔壁 保持紧密接触，知道胶结剂固化为止。最后进行套孔应力解除。在应力解除前后 各测一次应变读数，根据12支应变片的读数变化值来计算应力值。一个单孔应 力测量即可确定测点的三维应力大小和方向。CSIR孔壁应变计的适用孔为36 38mm。R.利曼和G.赫格特在实验室中进行了 CSIR孔壁应变计的测量精度测量，试 样为钢和吕立方体，施加的为单轴压缩应力，得出的测量误差为5 10%，迈尔 万(A.M.Myrvang)介绍了他们于1968年自行研制的和CSIR孔壁应变计相类似 的三轴孔壁应变计，其在挪威得到广泛应用。瑞典律勒欧大学的利金(B.A.Leijon) 等人介绍

28、了他们研制的孔壁应变计，其工作原理和CSIR孔壁应变计也很相似， 只是它的电阻应变花的结构作了一些改变，将原CSIR孔壁应变计中的每组应变 花分成纵向排列的两组，每组应变花有两支应变花，成十字形布置(见图5 12)。 此外哈尔耶恩(L.Hallbjorn)和萨尔(P.S.Sarkka)也道到了他们研制的孔壁应变 计。3.4岩体中的弹性变形能地下岩体处于复杂和强烈的自重应力和构造应力场中。地下煤层和岩层在应 力作用下，体积和形状都会发生变化，这种变形是外力做功的结果。岩体受外力 作用而产生弹性变形时，在岩体内部所储存的能量，称之为弹性应变能。不计变 形过程中的能量损耗时，外力做功会全部转变成弹性

29、能。.岩体的弹性变形能=外力所做的功如果弹性岩体的外力解除，则岩体的弹性能就要释放（在岩体是弹性的假设 下），就会产生一系列矿山压力现象，如冲击矿压等。3.4.1单向应力条件dxX由0 b时，做功为d e dV单元体的应变能仍等于外力所做的功设单元体各面上的应力按同一比例从则单元体x面上的外力为bx dydz。沿x方向的伸长量ex dx，若应力七是 由零增加到b，则b做功为b dy dz de dx dW = j b dy dz d e dx = j b单元体应变能二外力所做的功dVe = dW = j b x d e dV单位体积的应变能称为应变能密度dV eV =e dV 由虎克定律:.v

30、e = j3.4.2空间应力状态时零增加到最终值，在线弹性情况，由迭加原理，空间应力状态下的应变能密度为V = 1 b 8 +b 8 +b 8 821 12 23 3代入广义虎克定律:b - *(b +b )整理得:V8 = 2e b2 +b2 +b2 一 2*(bb +b b +b b )(2-19)一般情况下，外力作用下，岩体会同时发生体积和形状变化， 则 岩体内的弹性能=体积改变能+形状改变能。0、体积均匀变化产生的设单元体边长均为a，则dV= a 3应变能密度（V8）=体积改变能密度（VV）+形状改变能密度（Vd）v8= Vv+Vd则变形后单元体的体积为dV = (1 + 8 ) a

31、g(1 + 8 ) a g(1 + 8 ) a1123展开，略去高阶小量dV （1+8 + 8 +8 ）a31123则单元体的体积变化率0为0 = d = (8 +8 +8 ) dV 1230设为体应变 由虎克定律：0 1 E* （。1+。2+。3）体应变0只与三个主应变之和有关，记主应力平均值只引起体积改变，不 mb = (气 +b2 +b3)，若用 b 代替原 b 1、b2、b3 ;贝 b产生形状改变，且体积改变量仍等于原体积改变值，由（2 19） 体积改变能密度：V = 3 2*)。2 1 2* (b +b +b )2V 2E m 6 E1 2 3形状改变能密度：V =匕匕1+L/、,、

32、,=(b _b )2 + (。一b )2 + (。一b )26 E1223313.4.3岩体中的弹性变形能对于自重应力场条件，采深为H时，则原岩应力b 1 =y Hb =b =丫 H23 1 一日则单位体积的体积改变能U = 1 6E （b +b +b =y 2 H 2 f 1+* f6 E 1 - J（1 - 2日）（1 +日=Y2 H26 E （1-同理：形状改变能密度U =(1-2 产(1 丁)y 2 H 2 d3E (1-从上式可见，U,、Ud是采深H平方的增函数。H大，岩体中储存的的弹性 能迅速增加，发生突然破坏的可能性增加。如果岩体中有塑性变形量，且一般认为岩体的形状改变被塑性变形所吸收， 则岩体的变形能是体积改变能。且考虑构造应力的影响b =yH，b +b = 2人yH人 平均水平应力与垂直应力比值(1 _ 2)(1 + 2X)2则 U =Y 2 H 2V6 E在开挖边界，应力集中引起Uy迅速增加。