随机化方法在岩石力学中应用及研究综述.docVIP

随机化方法在岩石力学中应用及研究综述 　　摘要： 　　 　　岩石是地球上分布最广，总量最大（地球99%以上的材料构成是岩石）的一种自然地质材料[1]，它构成了地球四大圈层之一的岩石圈，因此与岩石有关的各项工程活动将与人类乃至其他生物的生活环境密切相关， 开展岩石力学及工程应用的研究具有重要的理论及现实意义。 　　关键词：岩石随机模糊 　　中图分类号：P58 文献标识码：A 文章编号： 　　自从1963年意大利Vajont大坝上游左岸的库岸边坡滑坡造成毁灭性的灾难以来，岩石的工程力学性质开始被人们所重视，并由此开始了现代岩石力学的工程应用研究。随着全球工业化的快速发展，世界各地都开展了许多大规模的岩石工程，涉及能源和资源开发的各个方面。在过去的半个世纪，人们更多的是将经典的力学方法及确定性理论应用于岩石力学的研究，在这方面已经取得了相当瞩目的科研成就[3]。由于岩石的组构、矿物成分，赋存的力学和地质环境的不同，所以岩石的物理力学性质相差极大，通常都不能作为均一的材料，除非研究对象的几何尺度远远超过岩石的特征尺度。这个特征尺度可以是岩石的微细观的REV或者是岩体的宏观REV尺度。 　　随着岩石力学及工程应用研究的深入，人们越来越发现确定性理论在描述岩石力学及某些工程问题时存在的不足，因此亟需引入随机模糊理论的研究方法，来描述岩石的几何、物理、力学以及工程性质[4]。这就需要从数学模型的角度，对岩石工程遇到的不确定因素进行分类，这些因素可分为3种[2]：随机性、模糊性以及未确知性。这些不确定性因素主要是由于岩石材料的非均一性、空间组成矿物材料的不确定性、以及反应岩石性质的信息、数据的不全面、不完整而导致的不确定性。目前接触较多、研究相对也比较成熟的是力学随机性问题。由此而专门发展出一套有关岩石力学性质的随机力学。另外，从工程角度进行分类，不确定因素还可体现在以下几方面[2]：荷载的不确定性、材料物理力学参数的不确定性、几何尺寸的不确定性、初始条件和边界条件的不确定性和计算模型的不确定性。 　　针对岩石力学与工程的这些不确定性问题，研究人员已经从多个方面开展了模糊岩石力学与随机岩石力学的研究。本文将总结这些方面的研究成果，以期对随机理论和模糊数学在岩石力学与工程中的应用有一个全面而系统的认识，并通过对文献的解读，指出目前研究的热点，存在的不足以及未来研究的展望。 　　2. 岩石及岩体裂隙结构的随机模拟 　　岩石及岩体都是赋存于自然地质条件中的初始损伤体。在岩石内部通常都会由于矿物颗粒的沉积作用以及动力变质作用而形成初始微裂隙，这些微裂隙的存在极大地改变和决定了岩石的微细观结构以及力学性质（包括岩石的初始各向异性和诱发各向异性）。在描述岩石微细裂隙结构方面，X光衍射仪和扫描电镜[5]是这方面的强有力研究手段。朱珍德 等人[6]引入随机数学及模糊数学理论，在Matlab开发平台上，采用Monte-Carlo随机理论描述大理石细观结构体积表征单元，通过计算机对其细观SEM扫描电镜照片进行二值化分析，提取出微裂隙的长度、方位角、宽度、面积和周长等细观信息，测定了大理石颗粒及矿物组构的几何尺度，并得到大理石微细观裂隙的几何尺寸的概率密度及分布函数形式，最后采用分形理论，得到岩石微结构的分形盒维数，并与Swoboda损伤理论相结合，描述了大理岩考虑损伤效应的应力应变关系。 　　岩体中包含大量的节理裂隙，这些节理裂隙形成的网络系统对岩体的力学性质，如：变形、破坏具有控制性的作用，它们是岩体非连续性、非均匀性和各向异性的根源。要明确工程岩体的变形和破坏模式必须对岩体的节理裂隙网络以及岩体的结构有比较清晰的认识，这就需要对节理裂隙进行产状和几何尺寸的现场统计和调查分析。随着电子计算机的发展，人们期望能够将现场统计的节理裂隙地质资料在计算机上再现，然后将建立的几何模型应用于数值模拟和计算分析。数值分析方法确定节理岩体等效参数的主要特点是充分利用现场地质调查、节理裂隙统计、室内及现场岩样和节理结构面试验等岩体已知信息，建立不同尺度的节理岩体裂隙网络模型[7]。这样就可以实现从现场调查的地质资料到计算机的模拟再现和三维地质模型的建立，进而利用连续等效应变和节理损伤张量理论，可以确定宏观岩体的等效力学参数。采用这种数值手段模拟任意尺度下的岩体力学试验，可以在一定程度上弥补不能开展大尺度岩体力学试验的问题。张贵科 等人[8]在前人的研究基础上，对节理岩体裂隙网络的计算机模拟进行了深入的研究，提出采用平面四边形来模拟岩体结构面的方法，这种处理与实际情况更为接近。采用解析代数的方法，[8]中给出了平面四边形的数学描述及其与空间点的关系，并在此基础上，研究了岩体的REV尺度，开发了岩体三维裂隙网络图形的模拟程序。图1是该软件生成的岩体结构图。张贵