地应力的计算
地应力的计算
《地壳应力随深度的变化规律》
1.水平主应力值随深度的增加而增加，通常比覆盖层静压大几倍，且远大于视岩休为弹性介
质的侧向约束，即按计算水平应力(式中为泊松比，为岩石密度，为重力加速度， Z 为深度)。
地壳中水平应力的另一个特点是其各向异性。也就是说，两个水平主应力(最大水平主应力) 及 (最小水平主应力) 的大小很少是相等的。根据我们的观测结果，在中国大
陆地壳中，最小与最大水平主应力的比值为0.3 一0.7 的约占70%，即一般最大水平主应力是最小水平主应力的1.4 一3.3倍。最大水平主应力与最小水平主应力随深度变化的梯度在不同地区是不同的。
《地壳应力在低渗裂缝砂岩油田开发中的应用》
水平主应力的总和与测点深度的关系式为：:式中
—水平最大主应力，MPa；
—水平最小主应力，MPa；
—测点深度，m；
—地面岩石中水平主应力的总和，MPa；
—应力梯度，MPa
套管抗外挤强度，注水后，注入水窜入泥岩层诱发地应力在井壁产生周向应力。计算公式为：
，最大最小周向应力。
在压裂施工中，当井内压裂液的压力升高到一定数值时，油层即发生破裂，这时油层承受的净压力，称为油层的破裂压力, 表达式:
—油层破裂压力，MPa；
，—油层最小、最大水平主应力，MPa；
—油层岩石抗张强度，MPa；
—油层孔隙压力，MPa；
当停止泵入压裂液，最小主应力将迫使裂缝闭合，当裂缝刚刚张开或恰恰没有闭合时，裂缝中压裂液所承受的净压力称为闭合压力，它近似等于油层的最小主应力。
。
《地形条件对大安山井田地应力的影响》
《断层活动与原地应力状态》——李方全
如果沿断层面的剪应力等于阻碍滑动的摩擦阻力时,在断层面上就会发生摩擦滑动,这就是库伦准则。也可用主应力来改写库伦准则,井引入有效应力概念。对于方位合适的断层面,最大、最小有效主应力之比可表示为摩擦系数μ的函数
若最大、最小有效主应力(式中为孔隙压力)之比小于此值,则断层稳定,不发生滑动.如果
比值等于此值,就会在方位合适的断层上发生滑动。所谓方位合适的断层,就是其断层面的法线方向与最大主应力的夹角为,中间应力位于断层面内,同时与与的关系由下式给
出:
对水平主应力的估计
根据上述关系,可以由岩石密度和摩擦系数及断层活动方式来估计水平应力。在正断层型的应力状态下,由于,得到,其，
其中。在逆断层情况下,由于,
得到。式中下为岩石容重,H为深度。在铅直应力为中间主应力的
情况下,
为最大水平主应力的上限,为最小水平主应力的下限。
断层活动危险性的评估
根据库伦淮则,,如果采用平均应力和最大剪应力,则,
,并考虑到孔隙压力的影响,可得,式中。因此根据原地应力实测资料和断层的凝聚力和摩擦系数,即可利用上述关系讨论一个具体断层的活动方式,并对断层活动的危险性作出估计。
水平应力随深度变化的规律
水平应力-深度曲线为一次函数曲线。其中最大水平主应力为，最小水平应力为。
岩石种类相关系数取值自由项取值岩浆岩最大主应力0.0318 5.8950
岩浆岩最小主应力0.0198 0.2325
沉积岩最大主应力0,。0240 4.9125
沉积岩最小主应力0.0183 1.5673
变质岩最大主应力0.0264 4.0567
变质岩最小主应力0.0194 1.6859
最大剪应力
最大剪应力对岩石的破坏起着重要的作用，一般情况岩石的出现裂缝和裂缝发展方向都由岩石所受最大剪应力决定。由材料力学相关知识可以知道，最大剪应力的确定可由岩块最大
正应力和最小正应力决定：,方向沿最大正应力和最小正应力夹角的角平分线。由于剪应力由正应力公式推导出，所以剪应力的表示形式与正应力表示形式基本相同，
岩性相关系数自由项R
岩浆岩0.0060 2.8313
沉积岩0.0029 1.6726
变质岩0.0035 1.1989
平均水平地应力与垂直地应力比值随深度变化的变化规律
若K为平均水平地应力与垂直地应力比值。平均水平地应力，则。
K的表达式可以简化为
a、b
岩性 a b
岩浆岩113.05 0.9520
沉积岩123.19 0.8042
变质岩96.547 0.7736
Geomech Abstr,1982,19:201～203.
McCutchen将地球假设为理想的球体，按照球队陈问题提出了一种地应力假说，认为水平应
力与垂直应力之比为,式中,,.B为积分常数，取决于计算基面的位置，R为地球半径，取R=6370Km。改假说得出结果为水平应力与垂直应力之比为深度的函数，通过选取计算初始条件，可对Browm和Hoek获得的全球地应力测量的结果给出较满意的解释。
Sheory P R.A theory for in situ stress in isotropic and transversely isotropic rock. Int J Rock Mech Min Sci Geomech Abstr,1994,31(1):23～34
Sheorey建立了一个静弹性热应力模型来估算地应力大小，该模型考虑了地壳曲率和岩层的弹性常数、密度和热膨胀系数的变化，基于此提出了一个估算水平应力与垂直应力之比k
的公式，式中z为埋深，m；为上覆岩层水平方向的平均变形模量，GPa。
（专业文档资料素材和资料部分来自网络，供参考。可复制、编制，期待你的好评与关注）
《岩石力学》地应力及其测量
1. 地壳是静止不动的还是变动的？怎样理解岩体的自然平衡状态？ 答：地壳是变动的。 自然平衡状态是指：岩体中初始应力保持不变的状态。 2. 初始应力、二次应力和应力场的概念。 答：未受影响的应力称为初始应力 工程开挖时，受工程开挖影响而形成的应力称为二次应力 地应力是关于时间和空间的函数，可以用“场”的概念来描述，称之为地应力场。 3. 何谓海姆假说和金尼克假说？ 答：海姆首次提出了地应力的概念，并假定地应力是一种静水应力状态，即地壳中 任意一点的应力在各个方向上均相等，且等于单位面积上覆岩层的重量，即???= ????=???? 金尼克认为地壳中各点的垂直应力等于上覆岩层的重量，而侧向应力（水平应力）是泊松效应的结果，其值应为乘以一个修正系数K。他根据弹性力学理论，认 为这个系数等于?? 1-??，即????=????，???=?? 1-?? ???? 4. 地应力是如何形成的？ 答：地应力的形成主要与地球的各种动力运动过程有关，其中包括：板块边界受压、地幔热对流、地球内应力、地心引力、地球旋转、岩浆侵入和地壳非均匀扩容等。 另外，温度不均、水压梯度、地表剥蚀或其他物理化学变化等也可引起相应的应力 场。 5. 什么是岩体的构造应力？构造应力是怎样产生的？土中有无构造应力？为什么？答：岩体中由于地质构造运动引起的应力称为构造应力。 关于构造应力的形成有两种观点：地质力学观点认为是地球自转速度变比的结果；大地构造学说则认为是出于地球冷却收缩、扩张、脉动、对流等引起的，如板 块边界作用力。 土中没有构造应力，由于土本身是各向同性介质，不存在地质构造。 6. 试述自重应力场与构造应力场的区别和特点。 答：由地心引力引起的应力场称为重力应力场，重力应力场是各种应力场中惟一能 够计算的应力场。地壳中任一点的自重应力等于单位面积的上覆岩层的重量，即????=????。 重力应力为垂直方向应力，它是地壳中所有各点垂直应力的主要组成部分，但 是垂直应力一般并不完全等于自重应力，因为板块移动，岩浆对流和侵入，岩体非 均匀扩容、温度不均和水压梯度均会引起垂直方向应力变化。 构造应力是由地质构造运动形成的。当前的构造应力状态主要由最近一次的构 造运动所控制，但也与历史上的构造运动有关。构造应力主要表现为以水平应力为 主，“在构造应力的作用仅影响地壳上层一定厚度的情况下，水平应力分量的重要性远远超过垂直应力分量。” 7. 岩体原始应力状态与哪些因素有关？ 答：地形地貌；岩体结构；岩石力学性质；地下水。 8. 简述地应力场的分布规律 答：1）地应力场的特性 （1）地应力场是一个以水平应力为主的三向不等压应力场 （2）地应力场是一个具有相对稳定性的非稳定应力场 2）垂直应力的分布规律 在深度为25~~2700m的范围内，????呈线性增长，大致相当于按平均容量??γ等于273kN???-3?计算出来的重力????。 3）水平应力的分布规律
地应力平衡的一个简单例子.
地应力平衡方法 熊志勇陈功奇 第一部分地应力平衡方法简介 地应力平衡有三种方法: (1 *initial conditions,type=stress,input=FileName.csv(或 inp 该方法中的文件 FILENAME.INP 获取方法为 :首先将已知边界条件施加到模型上进行正演计算 , 然后一般是将计算得到的每个单元的应力外插到形心点处并导出6个应力分量 (也可以导出积分点处的应力分量 , 视要求平衡的精确程度而定。其所采用的几何模型可以考虑地表起伏不平的情况以及岩土材料极其不均匀的情况 , 适用范围广。但由于外插的应力有一定误差 , 因此采用弹塑性本构模型时 , 可能会导致某些点的高斯点应力位于屈服面以外 , 当大面积的高斯点上的应力超出屈服面之后 , 应力转移要通过大量的迭代才能完成 , 而且有可能出现解不收敛的情况。在仅考虑自重情况下只能考虑受泊松比的影响带来的侧压力系数效应 , 因此平衡后的效果不一定很理想 , 但无疑其适用性很强。 (2 *initial conditions,type=stress,geostatic 该方法需给出不同材料区域的最高点和最低点的自重应力及其相应坐标。所采用的几何模型一般较规则 , 表面大致水平 , 地应力平衡的好坏一般只受岩体密度的影响 , 无论采用弹性或弹塑性本构模型都能很好的达到平衡 , 可以不必局限于仅受泊松比的影响 , 能够通过考虑水平两个方向的侧压力系数值来施加初始应力场。计算速度快 , 收敛性好。缺点就是不能够很好平衡具有起伏表面的几何模型 , 需知道平整后模型的上覆岩体自重。 (3 *initial conditions,type=stress,geostatic,user
地应力与地应力测量方法简介
地应力与地应力测量方法简介地应力，又称原岩应力，也称岩体初始应力或绝对应力，是在漫长的地质年代里，由于地质构造运动等原因产生的。在一定时间和一定地区内，地壳中的应力状态是各种起源应力的总和。主要由重力应力、构造应力、孔隙压力、热应力和残余应力等耦合而成,重力应力和构造应力是地应力的主要来源。地应力的形成主要与地球的各种动力运动过程有关，其中包括：板块边界受压、地幔热对流、地球内应力、地心引力、地球旋转、岩浆侵入和地壳非均匀扩容等。另外，温度不均、水压梯度、地表剥蚀或其他物理化学变化等也可引起相应的应力场。而重力作用和构造运动是引起地应力的主要原因，其中尤以水平方向的构造运动对地应力的形成影响最大。 地应力测量，就是确定拟开挖岩体及其周围区域的未受扰动的三维应力状态，这种测量通常是通过多个点的量测来完成的。地应力测量是确定工程岩体力学属性、进行围岩稳定性分析、实现岩土工程开挖设计和决策科学化的前提。地应力对矿山开采、地下工程和能源开发等生产实践均起着至关重要的作用，所以地应力研究是当前国际采矿界上的一个前沿性课题，近几十年来，世界上许多国家均开展了地应力的测量及应用研究工作,取得了众多的成果。 随着矿区开采现代化进程的不断提高和开采深度的不断增加，对矿区所处的地质条件和应力环境提出了更进一步的要求。查明矿区深部煤炭资源的开采地质条件和应力环境，为深部矿井的设计、建设和生产提供更加精细可靠的地质资料和数据，以便采取有效技术手段和措施，避免和减少灾害的发生，是实现矿井安全高效生产的重要保障。 地应力是引起采矿工程围岩、支架变形和破坏、产生矿井动力现象的根本作用力，在诸多的影响采矿工程稳定性因素中，地应力是最重要和最根本的因素之一。准确的地应力资料是确定工程岩体力学属性，进行围岩稳定性分析和计算，矿井动力现象区域预测，实现采矿决策和设计科学化的必要前提条件。 采矿规模的不断扩大和开采深度的纵深发展，地应力的影响越加严重，不考虑地应力的影响进行设计和施工往往造成露天边坡的失稳、地下巷道和采场的坍塌破坏、冲击地压等矿井动力现象的发生，致使矿井生产无法进行，并经常引起
地应力及其分布规律
地应力及其分布规律 1 、地应力的基本概念 地应力是存在于地层中的未受工程扰动的天然应力，也称岩体初始应力、绝对应力或原岩应力。广义上也指地球体内的应力。它包括由地热﹑重力﹑地球自转速度变化及其他因素产生的应力。 地应力是各种岩石开挖工程变形和破坏的根本作用力；是确定工程岩体力学属性，进行围岩稳定性分析，实现开挖设计和决策科学化的必要前提条件。此外地应力状态对地震预报、区域地壳稳定性评价、油田油井的稳定性、核废料储存、岩爆、煤和瓦斯突出的研究以及地球动力学的研究等也具有重要意义。 2、地应力的成因 产生地应力的原因是十分复杂的，地应力的形成主要与地球的各种动力运动过程有关，其中包括：板块边界受压、地幔热对流、地球内应力、地心引力、地球旋转、岩浆浸入和地壳非均匀扩容等。另外，温度不均、水压梯度、地表剥蚀或其它物理化学变化等也可引起相应的应力场。其中，构造应力场和自重应力场为现今地应力场的主要组成部分。 当前的地应力状态主要由最近的一次构造运动所控制，但也与历史上的构造运动有关。由于亿万年来，地球经历了无数次大大小小的构造运动，各次构造运动的应力场也经过多次的叠加、牵引和改造，另外，地应力场还受到其他多种因素的影响，造成地应力状态的复杂性和多变性， 地应力成因之一：地幔热对流（图1、图2） 地应力成因之一：板块边界受压（图3）
地应力成因之一：岩浆浸入（图4） 3、地应力的影响因素 地壳深层岩体地应力分布复杂多变，造成这种现象的根本原因在于地应力的多来源性和多因素影响，但主要还是由岩体自重、地质构造运动和剥蚀决定。1)岩体自重的影响 岩体应力的大小等于其上覆岩体自重，研究表明：在地球深部的岩体的地应力分布基本一致。但在初始地应力的研究中人们发现，岩体初始应力场的形成因素众多，剥蚀作用难以合理考虑，在常规的反演分析中，通常只考虑岩体自重和地质构造运动 2)地形地貌和剥蚀作用对地应力的影响 地形地貌对地应力的影响是复杂的，剥蚀作用对地应力也有显著的影响，剥蚀前，岩体内存在一定数量的垂直应力和水平应力，剥蚀后，垂直应力降低较多，但有一部分来不及释放，仍保留一部分应力数量，而水平应力却释放很少，基本上保留为原来的应力数量，这就导致了岩体内部存在着比现有地层厚度所引起的自重应力还要大很多的应力数值。 3)构造运动对地应力的影响 在地壳深层岩体，其地应力分布要复杂很多，此时由于构造运动引起的地应力对地应力的大小起决定性的控制作用。研究表明：岩体的应力状态，一般其铅垂应力分量是由其上覆岩体自重产生的，而水平应力分量则主要由构造应力所控制，其大小比铅垂应力要大得多。 4)岩体的物理力学性质的影响 从能量的角度看，地应力其实是一个能量的积聚和释放的过程。因为岩石中地应力的大小必然受到岩石强度的限制，可以说，在相同的地质构造中。地应力的大小是岩性因素的函数，弹性强度较大的岩体有利于地应力的积累，所以地震和岩爆容易发生在这些部位,而塑性岩体因容易变形而不利于应力的积累。 5)水、温度对地应力的影响 地下水对岩体地应力的大小具有显著的影响，岩体中包含有节理、裂隙等不连通层面，这些裂隙面里又往往含有水，地下水的存在使岩石孔隙中产生孔隙水压力，这些孔隙水压力与岩石骨架的应力共同组成岩体的地应力。温度对地应力的影响主要体现在地温梯度和岩体局部受温度的影响两个方面。由于地温梯度而产生的地温应力，岩体的温度应力场为静压力场，可以与自重应力场进行代数迭加，如果岩体局部寒热不均，就会产生收缩和膨胀，导致岩体内部产生应力。4、地应力的分布规律
地应力测量
地应力测量的国内外研究现状 0 引言 地应力（in-situ stress），又称原岩应力，也称岩体初始应力或绝对应力，是在漫长的地质年代里，由于地质构造运动等原因产生的。在一定时间和一定地区内，地壳中的应力状态是各种起源应力的总和。主要由重力应力、构造应力、孔隙压力、热应力和残余应力等耦合而成,重力应力和构造应力是地应力的主要来源。地应力的形成主要与地球的各种动力运动过程有关，其中包括：板块边界受压、地幔热对流、地球内应力、地心引力、地球旋转、岩浆侵入和地壳非均匀扩容等。另外，温度不均、水压梯度、地表剥蚀或其他物理化学变化等也可引起相应的应力场（雷化南，等译.1976）。而重力作用和构造运动是引起地应力的主要原因，其中尤以水平方向的构造运动对地应力的形成影响最大。因此，岩石中的原地应力是由主动施加的力和积蓄的残余应变两者引起的。 地应力测量（In situ stress measurement），就是确定拟开挖岩体及其周围区域的未受扰动的三维应力状态，这种测量通常是通过多个点的量测来完成的。地应力测量是一项综合性的测试，可以说任何一种单一的方法都不能很好地完成，往往需要几种方法结合起来对比使用，才可以保证结果的可靠性。即使如此，地应力测量中也往往会出现同一测点测量值分散的情况。 地应力测量是确定工程岩体力学属性、进行围岩稳定性分析、实现岩土工程开挖设计和决策科学化的前提。地应力对矿山开采、地下工程和能源开发等生产实践均起着至关重要的作用，所以地应力研究是当前国际采矿界上的一个前沿性课题，近几十年来，世界上许多国家均开展了地应力的测量及应用研究工作,取得了众多的成果。 1 地应力测量在国外发展概况及研究现状 人们最初对地应力概念的认识以及地应力测量技术的发展都源于早期的矿山工程建设，最早的原位地应力测量起始于20世纪30年代。1932年，美国人劳伦斯(Lieurace)在胡佛坝(HooverDam)下面的一个隧道中采用岩体表面应力解除法首次成功地进行了原岩应力的测量。此后，地应力测试技术一直停留在岩体表面应力测量上，发展十分缓慢，在20世纪50年代，哈斯特(Hast)采用应力解
地应力的计算
地应力的计算 《地壳应力随深度的变化规律》 1.水平主应力值随深度的增加而增加，通常比覆盖层静压大几倍，且远大于视岩休为弹性介 质的侧向约束，即按计算水平应力(式中为泊松比，为岩石密度，为重力加速度， Z 为深度)。 地壳中水平应力的另一个特点是其各向异性。也就是说，两个水平主应力(最大水平主应力) 及 (最小水平主应力) 的大小很少是相等的。根据我们的观测结果，在中国大 陆地壳中，最小与最大水平主应力的比值为0.3 一0.7 的约占70%，即一般最大水平主应力是最小水平主应力的1.4 一3.3倍。最大水平主应力与最小水平主应力随深度变化的梯度在不同地区是不同的。 《地壳应力在低渗裂缝砂岩油田开发中的应用》 水平主应力的总和与测点深度的关系式为：:式中 —水平最大主应力，MPa； —水平最小主应力，MPa； —测点深度，m； —地面岩石中水平主应力的总和，MPa； —应力梯度，MPa 套管抗外挤强度，注水后，注入水窜入泥岩层诱发地应力在井壁产生周向应力。计算公式为： ，最大最小周向应力。 在压裂施工中，当井内压裂液的压力升高到一定数值时，油层即发生破裂，这时油层承受的净压力，称为油层的破裂压力, 表达式: —油层破裂压力，MPa； ，—油层最小、最大水平主应力，MPa； —油层岩石抗张强度，MPa；
—油层孔隙压力，MPa； 当停止泵入压裂液，最小主应力将迫使裂缝闭合，当裂缝刚刚张开或恰恰没有闭合时，裂缝中压裂液所承受的净压力称为闭合压力，它近似等于油层的最小主应力。
。 《地形条件对大安山井田地应力的影响》 《断层活动与原地应力状态》——李方全 如果沿断层面的剪应力等于阻碍滑动的摩擦阻力时,在断层面上就会发生摩擦滑动,这就是库伦准则。也可用主应力来改写库伦准则,井引入有效应力概念。对于方位合适的断层面,最大、最小有效主应力之比可表示为摩擦系数μ的函数 若最大、最小有效主应力(式中为孔隙压力)之比小于此值,则断层稳定,不发生滑动.如果 比值等于此值,就会在方位合适的断层上发生滑动。所谓方位合适的断层,就是其断层面的法线方向与最大主应力的夹角为,中间应力位于断层面内,同时与与的关系由下式给 出: 对水平主应力的估计 根据上述关系,可以由岩石密度和摩擦系数及断层活动方式来估计水平应力。在正断层型的应力状态下,由于,得到,其， 其中。在逆断层情况下,由于, 得到。式中下为岩石容重,H为深度。在铅直应力为中间主应力的
地应力计算公式
地应力计算公式 （一）、井中应力场的计算及其应用研究（秦绪英，陈有明，陆黄生 2003年6月） 主应力计算 根据泊松比μ、地层孔隙压力贡献系数V 、孔隙压力0P 及密度测井值b ρ可以计算三个主应力值： ()001H v A VP VP μσσμ??=+-+??-?? ()001h v B VP VP μσσμ??=+-+??-?? H v b dh σρ=?? 相关系数计算： 应用密度声波全波测井资料的纵波、横波时差（p t ?、s t ?）及测井的泥质含量sh V 可以计算泊松比μ、地层孔隙压力贡献系数V 、岩石弹性模量E 及岩石抗拉强度T S 。 ① 泊松比 222 20.52() s p s p t t t t μ?-?= ?-? ② 地层孔隙压力贡献系数 222 22 (34) 12() b s s p m ms mp t t t V t t ρρ??-?=-?-? ③ 岩石弹性模量 22 22234s p b s s p t t E t t t ρ?-?= ? ??-? ④ 岩石抗拉强度 2 2 (34)[(1)]T b s p sh sh S a t t b E V c E V ρ=???-????-+?? 注：,,,m ms mp t t ρρ??分别为密度测井值，地层骨架密度，横波时差和纵波时差值。,,a b c 为地区试验常数。 其它参数 不同地区岩石抗压强度参数是参照岩石抗拉强度数值确定,一般是8～12倍,也可以通过岩心测试获得。岩石内摩擦系数及岩石内聚力是岩石本身固有特性参数,可以通过测试分析获得。地层孔隙压力由地层水密度针对深度积分求取,或者用重复地层测试器RFT 测量。也可以通过地层压裂测试获得,测试时,当井孔压力下降至不再变化时,为储层的孔隙压力。
有限元特辑II 初始地应力平衡
有限元特辑II初始地应力平衡 技术邻作者：闷油瓶 文章所包含相关领域及技术点：应力平衡、abaqus、有限元Motivation 在ABAQUS中，提供了5种定应力平衡方法，分别是 1.(AUTOBALANCE)自动平衡法; 2.*INITIAL CONDITIONS，TYPE=STRESS，GEOSTATIC; 3.*INITIAL CONDITIONS，TYPE=STRESS，FILE=file，INC=i nc; 4.*INITIAL CONDITIONS，TYPE=STRESS，INPUT=XX．DAT; 5.*INITIAL CONDITIONS，TYPE=STRESS，GEOSTATIC，USER。 以上5种方法并不是每一种适用于所有的岩土模型，方法从易到难。方法①的自动平衡法，它省去了自重应力以及生成相应初应力文件和导入的麻烦。在(GEOSTATIC)地应力中选择自动增量步就能使用自动地应力平衡功能，还能指定允许的位移变化容限。不过自动地应力平衡支持有限的几种材料如弹性，塑性等，而其起单元也有一定的要求。 方法②为关键字定义初始地应力法，这种方法需要给出不同材料区域的最高和最低点的自重应力及其相应坐标。所采用的几何模型一般较
规则，表面水平，能够通过考虑水平两个方向的侧压力系数值来施加初始应力场。关键字定义初始地应力法只适合土体表面水平的土体，而初始地应力提取法由于外插的应力有一定的误差，因此对于材料是弹塑性的复杂土体，应力转移要通过大量的迭代才能完成，而其有可能出现不收敛的情况，平衡效果可能不是很理想。 方法③是ODB导入法，这种方法可使用之前算过的ODB文件结果，也就是说提前计算一个初始应力ODB文件，定义初始应力时直接指定ODB文件即可。 方法④初始应力提取法，首先将已知边界条件施加到模型上进行计算，然后是将计算得到的每个单元的应力外插到形心点出，导出S1 1，S22，S33，S12，S13，S23六个应力分量。这种方法是最为通用的方法，可以实用于不同材料，不规则的地形，适用性强。 方法⑤是采用用户子程序SIGINI来定义初始应力场，可以定义其为应力分量为坐标、单元号、积分点号等变量的函数，要达到精确平衡需已知具体边界条件，在实际中应用较少。 N小问 I.为何要施加初始地应力？ II.什么工况下，施加初始地应力？ III.什么时间点施加地应力？ IV.施加过程中及之后需要的注意事项？
地应力平衡方法介绍
6.19 Enhancements to the geostatic procedure Products: Abaqus/Standard Abaqus/CAE Benefits: The geostatic procedure for obtaining the initial equilibrium state has been enhanced so that you no longer have to specify initial stresses that are close to the equilibrium state to obtain a solution corresponding to the original configuration. Description: The geostatic procedure is normally used as the first step of a geotechnical analysis; in such cases gravity loads (and possibly other types of loads) are applied during this step. Ideally, the loads and initial stresses should exactly equilibrate and produce zero deformations. However, in previous releases of Abaqus the geostatic procedure did not enforce this condition. In complex problems it may be difficult to specify initial stresses and loads that equilibrate exactly. Consequently, the displacements corresponding to the equilibrium solution might be large unless a special procedure is used to enforce small displacements. The enhanced geostatic procedure allows you to obtain equilibrium in cases when the initial stress state is unknown or is known only approximately. Abaqus automatically computes the equilibrium corresponding to the initial loads and the initial configuration, allowing only small displacements within user-specified tolerances. The procedure is available with continuum and cohesive elements with pore pressure degrees of freedom and the corresponding stress/displacement elements. The elastic, porous elastic, Cam-clay plasticity, and Mohr-Culomb plasticity material models are supported. Although the list of supported materials includes materials that exhibit inelastic behavior, the procedure is intended to be used in analyses in which the material response is primarily elastic; that is, inelastic deformations are small. The new enhancements are available from the Incrementation tabbed page when you create or edit a geostatic step in Abaqus/CAE. You must select automatic incrementation to access the new controls. The default settings for increment size and maximum displacement change are shown in Figure 6–6. Figure 6–6 The Incrementation options for a geostatic step.
地应力知识
地应力知识 简介 地应力是存在于地层中的未受工程扰动的天然应力，也称岩体初始应力、绝对应力或原岩应力。 随着水利水电、矿山、交通与城建等边坡、洞室及深基坑等事故的明显增加从而使人们对地应力引起较为广泛的注意与重视，所以，地应力研究不但具有重要的实际意义，而且具有重要的理论意义。 一地应力的成因 产生地应力的原因是十分复杂的，也是至今尚不十分清楚的问题。30多年来的实测和理论分析表明，地应力形成主要与地球的各种动力运动过程有关，其中包括： 板块边界受压、地幔热对流、地球内应力、地心引力、地球旋转、岩浆侵入和地壳非均匀扩容等。另外，温度不均、水压梯度、地表剥蚀或其它物理化学等也可引起相应的应力场，其中，构造应力场和重力应力场是现今地应力场的主要组成部分。 1大陆板块边界受压引起的应力场 以中国大陆板块为例，由于受到印度板块和太平洋板块的推挤，推挤速度为每年数厘米，同时受到西伯利亚板块和菲律宾板块的约束。在这样的边界条件下，包括发生变形，产生水平受压应力场。2地幔热对流引起的应力场 由硅镁质组成的地幔因温度很高，具有可塑性，并可以上下对流和蠕动。地幔热对流引起地壳下面的水平切向应力，在亚洲形成由孟加拉湾一直延伸到贝加尔湖的最低重力槽。 3由地心引力引起的应力场（也称为重力场） 重力场，是各种应力场中唯一能够计算的应力场。重力应力为垂直方向应力，是地壳中所有各点垂直应力的主要组成部分，但是垂直应力一般并不完全
等于自重应力，因为板块移动、岩浆对流和侵入、岩体非均匀扩容、温度不均和水压梯度均会引起垂直方向应力变化。 4岩浆侵入引起的应力场 岩浆侵入挤压、冷凝收缩和成岩，均在周围底层中产生相应的应力场，其过程也是相当复杂。熔融状态的岩浆处于静水压力状态，对其周围施加的是各个方向相等均匀压力，但是热的岩浆侵入后逐渐冷凝收缩，并从接触面界面逐渐向内部发展，不同的热膨胀系数及热力学过程会使侵入岩浆自身及其周围岩体应力产生复杂的变化过程。 岩浆侵入引起的应力场是一种局部应力场。 5地温梯度引起的应力场 地层的温度随着深度增加而升高，一般为a=3℃/100m。由于地温梯度引起地层中不同深度不相同的膨胀，从而引起地层中的压应力，其值可达相同深度自重应力的数分之一。6地表剥蚀产生的应力场 地壳上升部分岩体因为风化、侵蚀和雨水冲刷搬运而产生剥蚀作用。剥蚀后，由于岩体内的颗粒结构的变化和应力松弛赶不上这种变化，导致岩体内仍然存在着比由地层厚度引起的自重应力还要大得多的水平应力值。因此，在某些地区，水平应力除与构造应力有关外，还和地表剥蚀有关。 二地应力的研究观点 对地应力的研究已有一百多年的历史了，但总的说来，现在主要有三种观点： 1“静水应力式”分布的观点 它最早是海姆(Heim)于1878年提出的“静水压力”假说。 以后（1905~1912年），又提出相应的应力计算公式。1925年，金尼克也提出了弹性理论计算法及相应的公式。但事实表明，它们只能适用于一定的环境条件下，如，埋深较大的未受到扰动的地层。
第7章 初始地应力场的生成及应用
第7章初始地应力场的生成及应用 在土木工程或采矿工程领域中，初始地应力场的存在和影响不容忽略，它既是影响岩体力学性质的重要控制因素，也是岩体所处环境条件下发生改变时引起变形和破坏的重要力源之一。因此，要想较真实地进行工程模拟仿真，就必须保证初始地应力场的可靠性。初始地应力场生成的主要目的是为了模拟所关注分析阶段之前岩、土体已存在的应力状态。本章即介绍FLAC3D中初始地应力场的生成方法及应用。 本章重点： ?常用的初始地应力场生成方法 ?常见工程初始地应力场的生成 ?路基施工过程的模拟 7.1 初始地应力场生成方法 在FLAC3D中，初始应力场的生成办法较多，但通常用的是以下三种方法，即弹性求解法、改变参数的弹塑性求解法以及分阶段弹塑性求解法。下面将以表7-1所述简单模型为例，介绍这三种生成初始地应力场的方法。 表7-1 模型尺寸、土体密度及变形参数 1×1×2 1×1×2 2000 30 10 0.35 7.1.1 弹性求解法 初始地应力的弹性求解法生成是指将材料的本构模型设置为弹性模型，并将体积模量与剪切模量设置为大值，然后求解生成初始地应力场。例叙述的是采用该法生成上述简单模型的初始地应力场的过程。 例7.1弹性求解生成初始地应力场 new gen zone brick size 1 1 2 model elas prop bulk 3e7 shear 1e7 fix z ran z 0 fix x ran x 0 fix x ran x 1 fix y ran y 0 fix y ran y 1 ;开始一个新的分析 ;生成网格模型 ;设置弹性本构模型 ;设置体积模量和剪切模量 ;固定z=0平面所有节点z向速度;固定x=0平面所有节点x向速度;固定x=1平面所有节点x向速度;固定y=0平面所有节点y向速度;固定y=1平面所有节点y向速度
地应力与地应力测量方法简介
地应力与地应力测量方法简介 地应力，又称原岩应力，也称岩体初始应力或绝对应力，是在漫长的地质年代里，由于地质构造运动等原因产生的。在一定时间和一定地区内，地壳中的应力状态是各种起源应力的总和。主要由重力应力、构造应力、孔隙压力、热应力和残余应力等耦合而成,重力应力和构造应力是地应力的主要来源。地应力的形成主要与地球的各种动力运动过程有关，其中包括：板块边界受压、地幔热对流、地球内应力、地心引力、地球旋转、岩浆侵入和地壳非均匀扩容等。另外，温度不均、水压梯度、地表剥蚀或其他物理化学变化等也可引起相应的应力场。而重力作用和构造运动是引起地应力的主要原因，其中尤以水平方向的构造运动对地应力的形成影响最大。 地应力测量，就是确定拟开挖岩体及其周围区域的未受扰动的三维应力状态，这种测量通常是通过多个点的量测来完成的。地应力测量是确定工程岩体力学属性、进行围岩稳定性分析、实现岩土工程开挖设计和决策科学化的前提。地应力对矿山开采、地下工程和能源开发等生产实践均起着至关重要的作用，所以地应力研究是当前国际采矿界上的一个前沿性课题，近几十年来，世界上许多国家均开展了地应力的测量及应用研究工作,取得了众多的成果。 随着矿区开采现代化进程的不断提高和开采深度的不断增加，对矿区所处的地质条件和应力环境提出了更进一步的要求。查明矿区深部煤炭资源的开采地质条件和应力环境，为深部矿井的设计、建设和生产提供更加精细可靠的地质资料和数据，以便采取有效技术手段和措施，避免和减少灾害的发生，是实现矿井安全高效生产的重要保障。 地应力是引起采矿工程围岩、支架变形和破坏、产生矿井动力现象的根本作用力，在诸多的影响采矿工程稳定性因素中，地应力是最重要和最根本的因素之一。准确的地应力资料是确定工程岩体力学属性，进行围岩稳定性分析和计算，矿井动力现象区域预测，实现采矿决策和设计科学化的必要前提条件。 采矿规模的不断扩大和开采深度的纵深发展，地应力的影响越加严重，不考虑地应力的影响进行设计和施工往往造成露天边坡的失稳、地下巷道和采场的坍塌破坏、冲击地压等矿井动力现象的发生，致使矿井生产无法进行，并经常引起
地应力与地应力测量方法简介
3.1 地应力与地应力测量方法简介 地应力，又称原岩应力，也称岩体初始应力或绝对应力，是在漫长的地质年代里，由于地质构造运动等原因产生的。在一定时间和一定地区，地壳中的应力状态是各种起源应力的总和。主要由重力应力、构造应力、孔隙压力、热应力和残余应力等耦合而成,重力应力和构造应力是地应力的主要来源。地应力的形成主要与地球的各种动力运动过程有关，其中包括：板块边界受压、地幔热对流、地球应力、地心引力、地球旋转、岩浆侵入和地壳非均匀扩容等。另外，温度不均、水压梯度、地表剥蚀或其他物理化学变化等也可引起相应的应力场。而重力作用和构造运动是引起地应力的主要原因，其中尤以水平方向的构造运动对地应力的形成影响最大。 地应力测量，就是确定拟开挖岩体及其周围区域的未受扰动的三维应力状态，这种测量通常是通过多个点的量测来完成的。地应力测量是确定工程岩体力学属性、进行围岩稳定性分析、实现岩土工程开挖设计和决策科学化的前提。地应力对矿山开采、地下工程和能源开发等生产实践均起着至关重要的作用，所以地应力研究是当前国际采矿界上的一个前沿性课题，近几十年来，世界上许多国家均开展了地应力的测量及应用研究工作,取得了众多的成果。 随着矿区开采现代化进程的不断提高和开采深度的不断增加，对矿区所处的地质条件和应力环境提出了更进一步的要求。查明矿区深部煤炭资源的开采地质条件和应力环境，为深部矿井的设计、建设和生产提供更加精细可靠的地质资料和数据，以便采取有效技术手段和措施，避免和减少灾害的发生，是实现矿井安全高效生产的重要保障。 地应力是引起采矿工程围岩、支架变形和破坏、产生矿井动力现象的根本作用力，在诸多的影响采矿工程稳定性因素中，地应力是最重要和最根本的因素之一。准确的地应力资料是确定工程岩体力学属性，进行围岩稳定性分析和计算，矿井动力现象区域预测，实现采矿决策和设计科学化的必要前提条件。 采矿规模的不断扩大和开采深度的纵深发展，地应力的影响越加严重，不考虑地应力的影响进行设计和施工往往造成露天边坡的失稳、地下巷道和采场的坍塌破坏、冲击地压等矿井动力现象的发生，致使矿井生产无法进行，并经常引起
地应力检测(1)
1、地质雷达检测隧道支护情况 包括隧道衬砌厚度是否满足设计要求、钢筋保护层厚度是否满足设计要求、隧道衬砌钢筋布臵是否满足设计要求、隧道衬砌钢架布臵是否满足设计要求、隧道衬砌的密实情况（包括二衬背后脱空及初支背后空洞、不密实）。 评判标准：《公路工程质量检验评定标准》（GTG F80/1-2004）；参考《铁路隧道衬砌质量无损检测规程》（TB10233-2004）。 2、地应力检测 我国地应力测量试验和研究开始于20世纪50年代后期，迄今为止，地应力测量的主要方法虽然很多，但尚未形成统一的分类标准．根据测量数据特点的不同，地应力测量大体分为绝对应力测量和相对应力测量．前者主要是确定地壳应力背景值，即主应力的大小和方向；后者则是观测应力随时间变化的动态变化规律，通常也称为地应力监测．根据测量基本原理的不同，绝对应力测量方法又可分为直接测量法和间接测量法．所谓直接测量法就是利用测量仪器直接测量和记录各种应力量，并由这些应力量和原岩应力的相互关系直接换算得到原岩应力值．间接测量法则是借助某些传感元件或媒介，测量和记录与岩体相关物理量的变化（如密度、泊松比、弹性波速等变化），然后通过相应的公式换算间接得到原岩应力值．目前，较为常用的绝对应力测量方法主要有水压致裂法、声发射法、钻孔崩落法、套芯应力解除法、应变恢复法等．其中，前3种方法属于直接测量方法，后2种方法属于间接测量方法．相对应力测量方法包括压磁法、压容法、体应变法、分量应变法及差应变法等．我们采用水压致裂法 地应力测量存在的问题与展望：随着我国工程建设不断向深部发展，地应力测量及监测正面临着严峻的考验．与发达国家相比，尚存在许多问题与不足．首先，在宏观层面上存在的问题与挑战有：第一，测量和监测深度不足。目前，国际上最大地应力测量深度已达5100m．在德国的KTB深钻及美国的SAFOD计划中，应力测量深度一般达到2000～3000m；日本也建立了数10座深度为1000～3800m的深井观测台站．我国的绝大部分应力测量深度仅数百米，超过1000m的深井观测极为稀少，这严重制约了测量数据在空间上的代表性．第二，缺乏合理系统的地应力监测网络．我国虽然积累了大量的地应力测量数据，但数据分布不均且质量参差不齐，地应力监测台站少、布局不合理，
地应力平衡的一个简单例子
abaqus地应力平衡 先说为什么要施加地应力： 仅仅和密度相关 地应力平衡的正确时间点: 处理地应力平衡的三种方法: 地应力平衡方法 熊志勇陈功奇 第一部分地应力平衡方法简介 地应力平衡有三种方法： （1）*initial conditions,type=stress,input=(或inp) 该方法中的文件获取方法为:首先将已知边界条件施加到模型上进行正演计算,然后一般是将计算得到的每个单元的应力外插到形心点处并导出6个应力分量(也可以导出积分点处的应力分量,视要求平衡的精确程度而定)。其所采用的几何模型可以考虑地表起伏不平的情况以及岩土材料极其不均匀的情况,适用范围广。但由于外插的应力有一定误差,因此采用弹塑性本构模型时,可能会导致某些点的高斯点应力位于屈服面以外,当大面积的高斯点上的应力超出屈服面之后,应力转移要通过大量的迭代才能完成,而且有可能出现解不收敛的情况。在仅考虑自重情况下只能考虑受泊松比的影响带来的侧压力系数效应,因此平衡后的效果不一定很理想,但无疑其适用性很强。 （2）*initial conditions,type=stress,geostatic 该方法需给出不同材料区域的最高点和最低点的自重应力及其相应坐标。所采用的几何模型一般较规则,表面大致水平,地应力平衡的好坏一般只受岩体密度的影响,无论采用弹性或弹塑性本构模型都能很好的达到平衡,可以不必局限于仅受泊松比的影响,
能够通过考虑水平两个方向的侧压力系数值来施加初始应力场。计算速度快,收敛性好。缺点就是不能够很好平衡具有起伏表面的几何模型,需知道平整后模型的上覆岩体自重。 （3）*initial conditions,type=stress,geostatic,user 该方法采用用户子程序SIGINI来定义初始应力场,可以定义其为应力分量为坐标、单元号、积分点号等变量的函数,要达到精确平衡需已知具体边界条件,在实际中应用较少。 第二部分地应力平衡方法实例详解 地应力平衡是岩土工程数值模拟分析的重要的内容，为了让师弟师妹们快点上手，我利用第一种方法做一个较简单的模型，希望对大家有用。 一、模型描述： 3 / 2080m kg = ρ 7 10 1? = E 35 .0 = ν 二、地应力平衡过程 1.启动ABAQUS,单击Create Model Datebase 2.创建部件（Part） 在Part模块，单击创建部件按钮，弹 出如右图的对话框， 按图输入部件名： Part-soil； 采用二维模型选择2D Planar； Type选择可变型（Deformable）； (30,20) (0,20) 12m12m 60kpa 2 m 地基
地应力测量方法共6页
地应力测量方法 1.水压至裂法 水压致裂法地应力测试是通过在钻孔中封隔一小段钻孔，然后向封隔段注入高压流体，从而确定原位地应力的一种方法。水压致裂法的2种方法试验设备相同，都有封隔器、印模器，使用高压泵泵入高压液体使围岩产生新裂隙或使原生裂隙重张。 常规水压致裂法(HF法) HF法是从射井方法移植而来，假定钻孔轴向为1个主应力方向，岩石均质、各向同性、连续、线弹性，采用抗拉破坏准则，在垂直于最小主应力方向出现对称裂缝，其仅能测得垂直于钻孔横截面上的二维应力。在构造作用弱和地形平坦区，垂直孔所测结果可代表2个水平主应力，垂直应力约等于上覆岩体自重，裂缝方位为最大水平主应力方位。 HF法测试周期短，不需要岩石力学参数参与计算，适合工程初勘阶段，不需试验洞，可进行大深度测量，是目前惟一一种可直接进行深部地应力测定的方法。通过对HF法的改进，德国大陆科学深钻计划(KTB)在主孔6 000 m和9 000 m处已成功获得了地应力资料。HF法是一种平面应力测量方法，为获得三维应力，YMizutaI和M KuriyagawaE提出3孔交汇地应力测量，我国长江科学院和地壳所也进行了大量的测试。但研究表明，当钻孔轴向偏离主应力方向，其结果就有疑问，要精确获得三维地应力较困难。为此，文献[7]基于最小主应力破坏准则，对3孔交汇HF法测试理论进行了完善，其有助于提高测量结果的计算精度，但还有待足够的测量数据来验证。
原生裂隙水压致裂法(HTPF法) HTPF法是HF法的发展，其要求在含有原生节理和裂隙的钻孔段进行裂隙重张试验以确定原位应力。HTPF法假定裂隙面是平的，且面上应力一致。对于深孔三维地应力直接测量，HTPF法可进行大尺度的地壳地应力测试，很有发展前途。HTPF法同HF法相比，假设少，不需考虑岩石 破坏准则和孔隙水压力，在单孔中便可获得三维地应力。但用HTPF法测试费时，且裂隙产状和位置的确定误差都可降低计算精度。 2.套钻孔应力解除法 套钻孔应力解除法根据解除方式和传感器的安装部位分为探孔应力解除法、孔底应变解除法和孔壁切割解除法。探孔应力解除法根据传感器的类型可分为孔壁应变法和孔径变形法。 孔壁应变法 孔壁应变法基于岩石各向同性、均质、连续、线弹性的假设，通过孔壁6个以上不同方向的应变值来计算岩体的三维地应力。孔壁应变法又可分为直接粘贴方法和包体方法。CSIR型三轴应变计就是将应变元件直接贴到孔壁中。空心包体是将应变元件贴到薄筒壁中，再用胶将薄筒和孔壁粘结。还有一种实心圆柱式包体技术，由于受包体材料和岩石物理力学性质差异影响大，已基本不用。 孔壁应变法最大的优点是单孔单点可准确测量岩体的三维地应力，缺点是:对岩石的完整性要求高，岩芯解除长度大于40～60 cm，并且在岩芯易饼化时测试很难成功；存在应变元件的粘贴、防潮、全过程测量和定向等问题;受温度变化、岩性差异影响大，测量结果离散性大。
地应力的测量方法
地应力的测量原理 目前地应力测量方法有很多种，根据测量原理可分为三大类： 第一类是以测定岩体中的应变、变形为依据的力学法，如应力恢复法、应力解除法及水压致裂法等； 第二类是以测量岩体中声发射、声波传播规律、电阻率或其他物理量的变化为依据的地球物理方法； 第三类是根据地质构造和井下岩体破坏状况提供的信息确定应力方向。其中，应力解除法与水压致裂法得到比较广泛的应用，其他几种只能作为辅助方法。 1.应力解除法测试原理和技术 1.1应力解除法测试原理 具有初始应力的岩体，用人为的方法卸去其应力，在岩体恢复变形的过程中测试其应变，然后用弹性力学理论计算出地应力的大小，得出其方向、倾角。目前国内外地应力测量普遍采用空心包体应变计测量技术。KX一81型空心包体应变计由A、B、C 3组共12枚应变片嵌埋在1个壁厚约3 mm的空心环氧树脂圆筒中间，圆筒外表面与钻孔壁用专用环氧树脂胶黏结在一起，其是在澳大利亚CSIRO空心包体应变计的基础上研制出来的，是套钻孔应力解除法的一种，只需1个孔就能测量出某点的三维原岩应力，具有使用方便、安装操作简单、成本低、效率高等优点。 1.2完全温度补偿技术 KX一81型空心包体应变计与其他许多应变测量仪器一样，均采用应变计作为敏感元件，并根据惠斯顿电桥的原理13J，将应变的变化转换成电压变化经放大后记录下来。电阻应变计对温度变化是很敏感的，温度发生变化时应变计的电阻值将发生变化，从而产生虚假的附加应变值。因此在现场测试中必须采取温度补偿措施。 惠斯顿电桥原理：平衡时，检流计所在支路电流为零，则有，（1）流过R1和R3的电流相同（记作I1），流过R2和R4的电流相同（记作I2）。（2）B，D两点电位相等，即UB=UD。因而有 I1R1=I2R2；个阻值已知，便可求得第四个电阻。测量时，选择适当的电阻作为R1和R2，用一个可变电阻作为R3，令被测电阻充当R4，调节R3使电桥平衡，而且可利用高灵敏度的检流计来测零，故用电桥测电阻比用欧姆表精确。电桥不平衡时，G的电流IG与R1，R2，R3，R4有关。利用这一关系也可根据IG及三个臂的电阻值 求得第四个臂的阻值，因此不平衡电桥原则上也可测量电阻。在不平衡电桥中，G应从“检流计’改称为“电流计”，其作用而不是检查有无电流而是测量电流的大小。可见，不平衡电桥和平衡电桥的测量原理有原则上的区别。利用电桥还可测量一些非电学量。 1）根据惠斯顿电桥的原理自行设计并制成1个应变一电阻一电压转换装置，在每一桥路中，除工作应变桥臂外，其他3个桥臂均为电阻，其温度系数为1×10.6／℃，这样电阻在温度变化1℃时只产生5 X 10～P变化，从而可以忽略不计。 2)增加1个热敏电阻，在应力解除过程中连续不断地测量测点的温度变化。 3)在每一次应力解除完成后，进行温度、应变标定试验，为计算地应力给出正确的测量数据。 测点的布置 测点布置 测点应布置在裂隙、孔隙少且均匀致密的完整岩体中，且不受开采影响的区域，一般选择在开拓巷道或专门硐室内布置测试钻孔。钻孔要施工到巷道或硐室扰动应力场范围之外，避开巷道和采场的弯、叉拐、顶部等应力增高区，保证应力测点处于原岩应力区，钻孔深度一般