岩石介质细观非均匀性对宏观力学行为的影响.docxVIP

岩石介质细观非均匀性对宏观力学行为的影响 1 非均匀材料介质破裂规律 岩石不仅是一个经历了漫长地质作用并具有良好耐候性的地球环境，也是一个非常不均匀的介质，由裂缝和节理等缺陷组成。在外部载荷的作用下，岩屑内部缺陷的扩展和相互作用决定了岩屑的宏观变形和破坏性质。过去，大部分工作集中在均匀介质假设的基础上，但忽视了岩屑自身的不均匀性，以及变形的局部破坏。 地震能否预测的问题,实际上与地壳介质的非均匀性有直接的关系. 结构均匀与非均匀的介质中产生破裂的过程很不一样. 均匀介质中破裂的产生是突然的,而在非均匀介质中,在主破裂之前可以观测到前兆现象. 因此,要了解地震的发生过程或破裂过程,就必须研究非均匀材料介质破裂的发生方式. 前震和前震序列, 是公认的重要地震前兆. 在中国大陆上有前震和前震序列的地震占有一定的比例. 据李大鹏等统计, 1966～1992年MS≥ 5.5的165次地震中有前震序列的为6 %, 其中15次MS≥ 7.0地震中有前震序列的占20 %. 文献指出,根据对世界范围内MS≥7.0地震的统计显示, 前震(主震前40天,100km范围内)发生概率为42.2 %. 在美国圣安德烈斯断层(走滑段)的16次地震中, 有44 %的地震在主震前1天, 5km范围内有前震发生. 因此对前震序列进行深入研究, 对强震的预测具有重要的意义. 地震序列形成机理的研究引起人们的高度重视,并从岩石破裂实验和数学解析方法对其进行了探讨. Mogi对4种不同均质度的岩石样品进行岩石破裂实验研究,从介质较均匀到非常不均匀分别得到了三种地震序列类型:(1) 主震型; (2) 前震-主震-余震型; (3) 震群型. 唐春安等基于岩石的声发射与岩石的损伤具有一致性的假设, 用解析的方法, 运用统计细观损伤力学, 研究了单轴加载条件下岩石加载系统各种参数的变化, 探讨了岩石在不同加载方式下 (柔性加载、刚性加载、伺服控制加载等)的声发射特征. 得到岩石声发射 3种基本模式:(1)主震型;(2)前震-主震-余震型;(3) 震群型. 然而, 由于解析解的方法只适应于区域形状与边界条件都比较简单的情况, 不能考虑更复杂的工程问题或地质条件, 而且很难考虑单元之间的相互影响, 以及岩石破裂实验条件的限制. 因此, 进一步的工作就要借助于数值计算的方法. 本文对5种具有不同均质度系数和5种具有相同均质度系数的岩石试样模型进行了系列单轴压缩载荷作用下的数值试验,研究了介质细观非均匀性以及细观结构随机性对宏观力学行为和地震序列类型的影响. 2 单元破坏强度 本文所用RFPA2D系统, 是一个能模拟岩石介质逐渐破坏过程的数值模拟工具. RFPA2D包括两个方面的功能:应力分析和破坏分析. RFPA2D的应力分析采用有限元法. 破坏分析则是根据一定的破坏准则来检查材料中是否有单元破坏. 对破坏单元则采用刚度特性退化(处理分离)和刚度重建(处理接触)的办法进行处理. 为了模拟试验机加载情况, 采用位移加载方式. 对于每一步给定的位移增量, 首先进行应力计算. 然后根据破坏准则来检查模型中是否有破坏单元. 如果没有, 继续增加一个位移增量, 进行下一步应力计算. 如果有破坏单元, 则根据单元的拉或剪破坏状态进行刚度退化处理. 在本研究中,认为基元破坏前的力学性质为弹脆性,破坏后则具有残余强度特性. 基元破坏表现为刚度和强度的弱化. 然后重新进行当前步的应力计算. 重复上述过程, 直到整个材料产生宏观破坏. 由于单元的破坏为脆性, 因此认为单元破坏释放的弹性能以声发射的形式释放, 据声发射特性来观察岩石破裂过程. 考虑到岩石类脆性材料的抗拉强度远小于抗压强度, 因此本文采用了修正后的库仑准则、包含拉伸截断作为单元破坏的强度判据. 根据这一准则, 单元的破坏可能是被拉坏或剪坏的,表达式为: ???????σ1?(1+sinφ)(1?sinφ)σ3≥σc(σ3≥0)或σ3≤?σt(σ3≤0){σ1-(1+sinφ)(1-sinφ)σ3≥σc(σ3≥0)或σ3≤-σt(σ3≤0) (1) 其中,σc为单轴抗压强度,σt为单轴抗拉强度,φ为摩擦角,σ1、σ3为最大、最小主应力. 3 岩样基元力学性质 为了研究材料的非均匀性对岩石破裂过程的影响以及对微震序列类型的影响,设定岩石的力学性质(单元的弹性模量和强度性质)服从Weibull分布?c(m,α),其中m为形状参数,反映岩石材料力学性质的均质度. 图1反映了均质度m为1.1,1.5,2,3,5的材料基元力学性质的分布形式. 可以看出m越大,表明岩石的性质越均匀.α为反映岩石材料平均性质的参数,本文指模型材料的弹性模量和强度,采用5种具有不同均质度系数岩样研究非均匀性对破裂过程的影响(表1),以及5种相同试样均质度相同情况下岩石细