矿山岩体的原岩应力及其重新分布

矿山岩体的原岩应力及其重新分布第一节岩体中的原岩应力地壳中没有受到人类工程活动（如矿井中开掘巷道等）影响的岩体称为原岩体，简称原岩。存在于地层中未受工程扰动的天然应力称为原岩应力，也称为岩体初始应力、绝对应力或地应力。天然存在于原岩内而与人为因素无关的应力场称为原岩应力场。一、自重应力设岩体为半无限体，地面为水平面，在距地表深度为H处，任意取一单元体（图2-1），其上作用的应力为σz，σy，σx，形成岩体单元的自重应力状态。图2-1岩体单元体所在位置及其应力状态单元体上所受的垂直应力σz等于单元体上覆岩层的重量，（2-1）式中γ——上覆岩层的平均重力密度（Kn/m3）；H——单元体距离地表的深度（m）。在均匀岩体内，岩体的自重应力状态为(2-2)式中λ为常数，称为侧压系数。在岩体自重应力场内，垂直应力σz和水平应力σx、σy都是主应力。假设岩体为各向同性的弹性体，据广义虎克定律，单元体各方向的应变为（2-3）由于εx=0，εy=0，σx=σy，σz与σx、σy之间的关系为（2-4）（2-5）岩体初始应力状态的静水压力理论认为，在埋藏较深条件下，垂直压应力相当大，岩石呈现明显的塑性。泊松比μ近似等于0.5，侧压系数λ为1.0，此时（2-6）深部的岩体自重应力场达到静水应力状态。二、构造应力构造应力是由于地壳构造运动在岩体中引起的应力，岩体构造应力可以分为现代构造应力和地质构造残余应力。一般情况下地壳运动以水平运动为主，构造应力主要是水平应力；②构造应力分布不均匀，在地质构造变化比较剧烈的地区，最大主应力的大小和方向往往有很大变化。构造应力具有明显方向性，最大水平主应力和最小水平主应力之值一般相差较大。图2-2由地质特征推断构造应力方向（a）～（e）均为平面图图2-3世界各国垂直应力σv与深度H变化规律图三、原岩应力分布的基本规律(1)实测垂直应力基本上等于上覆岩层重量对全世界有关实测垂直应力的统计资料表明，在深度为25～2700m范围内，垂直应力σv呈线性增长，大致相当于按平均容重γ等于27kN/m3计算出来的重力γH，如图2-3所示。在绝大多数测点都发现确有一个主应力接近垂直方向，其偏差不大于200。(2)水平应力普遍大于垂直应力根据国内外实测资料统计，水平应力σh多数大于垂直应力σv，最大水平应力σh,max与垂直应力σv的比值，一般为0.5～5.5，很多情况下比值大于2，最大可达到30，最大水平主应力和最小水平主应力的平均值σh,av与σv相比，一般仍为0.5～5.5，大部分在0.8～1.5之间。这说明在浅层地壳中平均水平地应力也普遍大于垂直应力，垂直应力在多数情况下为最小主应力，在少数情况下为中间主应力。(3)平均水平应力与垂直应力的比值随深度增加而减小平均水平应力σh,av与垂直应力σv的比值λ是表征地区原岩应力场特征的指标，该值随深度增加而减小。但在不同地区，变化的速度不相同。霍克和布朗用回归出下列公式，表示比值的变化范围：(2-7)式中H为深度，单位为m。(4)水平应力具有很强的方向性最大水平主应力和最小水平主应力一般数值相差较大，最小水平主应力和最大水平主应力之比σh,min/σh,max一般为0.2～0.8，多数情况下为0.4～0.8。我国部分矿井原岩应力的实测值见表2-1。表2-1部分矿井原岩应力实测结果矿名深度/m主应力应力值/MPa方位角/0倾角/0τmax/MPa鹤壁六矿447.66σ132.50258.0712.3514.22σ222.16101.5475.56σ34.07349.215.19焦作九里山矿318.70σ18.85357.5974.803.90σ25.70192.0814.53σ31.06101.128.9涟邵牛马司矿556.57σ121.7887.9918.569.96σ217.43282.2370.89σ31.86170.464.38北票冠山矿989.05σ152.9634.7911.2519.86σ230.85148.2463.45σ313.23299.7723.71新汶孙村矿870σ138.13100.124.2σ228.3579.261.5σ31.6114.814.1第二节岩体中的弹性变形能岩体受外力作用而产生弹性变形时，在岩体内部所储存的能量，称为弹性应变能。一、单向应力条件单元体在单向应力状态下（图2-4a），设单元体各边长分别为dx，dy，dz。在x面上作用的力为σxdydy（对单元体而言看作外力），沿x方向的伸长为εxdx（图2-4b）。当应力有一增量时，相应的变形增量为dεxdx，在单元体上力所作的功为（σxdydz）（dεxdx）。应力由零逐渐增加至最终值，应力在单元体上所作的功可用下式表示：（2-8）单元体的应变能dVε在数值上等于力所作的功(2-9)单位体积的应变能称为应变能密度,记作(2-10)对于线弹性体，其应力与应变成正比，(2-11)图2-4单元体在单向应力作用下变形状态图2-5单元体在三向应力作用下体积变化二、空间应力状态（2-12）单元体的畸变能密度vd等于总的应变能密度vε与体积应变能密度vV之差。（2-13）三、岩体中的弹性变形能在自重应力场中，对于深度为H的开采条件，岩体所受到的应力为单位体积岩体体积改变能为（2-14）单位体积岩体畸变能为（2-15）考虑构造应力影响，在原岩应力场中，对于深度为H的开采条件，一般情况下岩体所受到的应力为式中λ—平均水平主应力与垂直应力的比值。单位体积岩体体积改变能为（2-16） 第三节“孔”周围的应力分布一、双向等压应力场内的圆形孔（一）基本假设围岩为均质，各向同性，线弹性，无蠕变或粘性行为；原岩应力为各向等压（静水压力）状态；巷道断面为圆形，在无限长的巷道长度里，围岩的性质一致。并且埋深H大于或等于20倍的巷道半径R0（或其宽、高）。即有（2-17）研究表明，当埋深H≥20R0时，忽略巷道影响范围（3～5倍的R0）内的岩石自重（图2-6），与原问题的误差不超过10%。于是，水平原岩应力可以简化为均布的。这样，原问题就构成荷载与结构都是轴对称的平面应变圆孔问题（图2-7）。图2-6深埋巷道的力学特点图2-7轴对称圆巷的条件（二）计算结果(2-18)(2-19)式中r1——孔的半径。主要结论：①在双向等压应力场中，圆孔周边全处于压缩应力状态。②应力大小与弹性常数E、μ无关。③σt、σr的分布和角度无关，皆为主应力，即切向和径向平面均为主平面。④双向等压应力场中孔周边的切向应力为最大应力，其最大应力集中系K=2，且与孔径的大小无关。当σt=2γH超过孔周边围岩的弹性限时，围岩将进入塑性状态。⑤其它各点的应力大小则与孔径有关。若定义以σt高于1.05σ1或σr低于0.95σ1为巷道影响圈的边界，则σt的影响半径，工程上有时以10％作为影响半径，则σt的影响半径Ri≈3r1。有限元计算常取5r1的范围作为计算域。⑥由公式（2-32）和（2-33）可知，在双向等压应力场中圆孔周围任意点的切向应力σt与径向应力σr之和为常数，且等于2σ1。二、双向不等压应力场内的圆形孔（一）双向不等压应力场内的圆形孔应力解根据弹性理论，双向应力无限板内圆形孔的应力解为：(2-20)（2-21）（二）讨论(1)若取极限情况λ=0，则有（2-22）（2-23）当时，圆孔的顶.底部出现了拉应力区；在圆孔两侧的最大应力集中系数值达到3。(2)只考虑自重情况下原岩应力状态的侧向应力系数在0与1之间，即0≤λ≤1。考虑λ=0；1/7；1/2；1在θ=00；900；1800；2700时的应力分布，圆孔两侧的切向应力集中系数处于2～3之间。(3)λ=1/3可得切向应力为取θ=900、θ=2700，则周边出现=0，即此时圆孔顶与底部不会出现拉应力。图2-8λ=0，1/7，1/2，1时，圆孔周围应力分布图2-9深埋椭圆巷道由上述讨论可见，>1/3,周边不出现拉应力；<1/3时，将出现拉应力；1/3，圆孔顶部与低部不出现拉应力。λ=0时，θ=900处，拉应力最大。所以，λ=0为最不利情况。1为均匀受压的最有利于稳定情况。三、椭圆形孔周边的应力分布在一般原岩应力状态（图2-9）下，深埋椭圆巷道周边切向应力计算公式为（2-24）(1)等应力轴比（理查兹（R.Richards）,1978）由得（2-25）（2-26） 时，，，最佳断面为圆形（圆是椭圆的特例）； 时，，，最佳断面为的竖椭圆。时，，，最佳断面为的横（卧）椭圆。(2)应力（无拉应力）轴比。孔两侧的最大切向应力将随孔的几何尺寸发生变化，其切向应力集中系数：当时,，当时,。显然，孔越扁，则应力集中系数越大。例如=2:1，则同理，分析θ=π/2、θ=3π/2时的情况。可知在时，即形成拉应力。当时，。将不同侧压系数和不同轴长比的椭圆形孔周边切向应力集中系数列于表2-2。表2-2椭圆形孔周边切向应力集中系数a/bθ1/21/1.51/11.5/12/110901.004.001.333.002.002.003.001.334.001.000.250901.750.252.0802.75-0.253.75-0.424.75-0.5000902.00-1.002.33-1.003.00-1.004.00-1.005.00-1.00四、矩形孔和其它形状巷道周边的应力分布表2-3矩形巷道周边切向应力部分计算结果表a:b=5a:b=3.2a:b=1.8a:b=1(正方形)附图λp0P0λp0P0λp0P0λp0P0001.192-0.9401.342-0.981.200-0.8011.472-0.8084503.3520.8213.0003.0005001.158-0.6442.392-0.1932.7632.7470.9803.8606502.6927.0306.201-0.5995.260900-0.6782.420-0.7702.152-0.3342.030-0.8081.472注：表格内的数字分别表示λp0、P0对该点的应力集中影响系数。五、存在多个孔时，孔周围的应力分布断面相同的相邻两个孔的应力分布以双向等压应力场中的圆形孔为例，若相邻两孔的间距>2Ri，则此两孔就不会产生相互影响，巷道周边的应力分布也将和单一孔的情况基本相同。在这种情况下，即使存在多条巷道，它们之间相互也不产生影响。大小不等的相邻两孔的应力分布图2-10所示为不等径相邻两孔的切向应力分布图。从图中可以看出，小孔周边的切向应力集中系数高达4.26，而大孔周边的应力集中系数仅为2.75。这说明大孔对小孔的应力分布影响较大，而小孔对大孔的影响则甚微。这个特点对于研究回采工作面与邻近巷道的相互影响很有参考价值。在同一水平多孔相互影响条件下的应力分布图2-11所示为λ=0条件下，同一水平多孔的相互影响。由图可以看出，孔周边的应力集中系数是随D／B值的增大而增大的(D为孔径，B为孔周边的间距)。另一方面又受同一水平上孔的数目影响。显然，孔的数目愈多，孔周边的应力集中系数也愈大。图2-10不等径相邻两孔的切向应力分布图图2-11多孔对应力集中系数的影响六、回采空间周围应力重新分布西德埃森采矿研究中心对回采空间周围应力分布进行了模拟研究作出应力分布图(图2-12)。图中表明，两个空间相衔接的拐角处应力值最高。图2-13为这种空间周边切向应力分布的实例。开采深度为1000m，原岩垂直应力为σz=25MPa。综上所述，在假设孔周围都处于弹性状态的条件下，应力重新分布有以下一些特点：(1)孔周围形成了切向应力集中，最大切向应力发生在孔周边。对圆形和椭圆形孔，最大切向应力发生在孔两帮中点和顶底的中部。对矩形孔，则最大切向应力发生在四角处。(2)应力集中系数的大小，对单一孔来说，圆形孔仅与侧压系数有关，其值k=2～3。对椭圆形孔，则不仅与有关，还与孔的轴长比有关，一般当a/b=2，λ=0～1时，k=4～5。对多孔来说，k值升高是由于单孔应力分布迭加作用的结果，其值视孔的大小和间距以及原岩应力场的侧压系数值而定。如图2-21所示，在前后两个回采空间的影响条件下，中间巷道所在地点的应力集中系数可达7，有时可能更大。(3)不论何种形状的孔，它周围的应力重新分布(主要是指切向应力分布)从理论上说影响是无限的,但从影响的剧烈程度来看多都有一定的影响半径。通常，可取切向应力值超过原岩垂直应力5％处作为边界线。(4)孔的影响范围与孔的断面大小有关。图2-12两相邻回采空间周围的应力分