Application of the numerical manifold method for stress wave propagation across rock masses

In this paper, the numerical manifold method (NMM) is extended to study wave propagation across rock masses. First, improvements to the system equations, contact treatment, and boundary conditions of the NMM are performed, where new system equations are derived based on the Newmark assumption of the space–time relationship, the edge‐to‐edge contact treatment is further developed for the NMM to handle stress wave propagation across discontinuities, and the viscous non‐reflection boundary condition is derived based on the energy minimisation principle. After the modification, numerical comparisons between the original and improved NMM are presented. The results show that the original system equations result in artificial numerical damping, which can be overcome by the Newmark system equations. Meanwhile, the original contact scheme suffers some calculation problems when modelling stress wave propagation across a discontinuity, which can be solved by the proposed edge‐to‐edge contact scheme. Subsequently, the influence of the mesh size and time step on the improved NMM for stress wave propagation is studied. Finally, 2D wave propagation is modelled, and the model's results are in good agreement with the analytical solution. Copyright © 2013 John Wiley & Sons, Ltd.     

Fault geometry and mechanics of marly carbonate multilayers: An integrated field and laboratory study from the Northern Apennines,Italy

Sealing layers are often represented by sedimentary sequences characterized by alternating strong and weak lithologies. When involved in faulting processes, these mechanically heterogeneous multilayers develop complex fault geometries. Here we investigate fault initiation and evolution within a mechanical multilayer by integrating field observations and rock deformation experiments. Faults initiate with a staircase trajectory that partially reflects the mechanical properties of the involved lithologies, as suggested by our deformation experiments. However, some faults initiating at low angles in calcite-rich layers (θ_i = 5°–20°) and at high angles in clay-rich layers (θ_i = 45°–86°) indicate the important role of structural inheritance at the onset of faulting. With increasing displacement, faults develop well-organized fault cores characterized by a marly, foliated matrix embedding fragments of limestone. The angles of fault reactivation, which concentrate between 30° and 60°, are consistent with the low friction coefficient measured during our experiments on marls (μ_s = 0.39), indicating that clay minerals exert a main control on fault mechanics. Moreover, our integrated analysis suggests that fracturing and faulting are the main mechanisms allowing fluid circulation within the low-permeability multilayer, and that its sealing integrity can be compromised only by the activity of larger faults cutting across its entire thickness.     

Displacement velocity effects on rock fracture shear strengths

Triaxial shear tests are performed to assess the effects of displacement velocity and confining pressure on shear strengths and dilations of tension-induced fractures and smooth saw-cut surfaces prepared in granite, sandstone and marl specimens. A polyaxial load frame is used to apply confining pressures between 1 and 18 MPa with displacement velocities ranging from 1.15 × 10⁻⁵ to 1.15 × 10⁻² mm/s. The results indicate that the shearing resistances of smooth saw-cut surfaces tend to be independent of the displacement velocity and confining pressure. Under each confinement the peak and residual shear strengths and dilation rates of rough fractures increase with displacement velocities. The sheared-off areas increase when the confining pressure increases, and the displacement rate decreases. The velocity-dependent shear strengths tend to act more under high confining pressures for the rough fractures in strong rock (granite) than for the smoother fractures in weaker rocks (sandstone and marl). An empirical criterion that explicitly incorporates the effects of shear velocity is proposed to describe the peak and residual shear strengths. The criterion fits well to the test results for the three tested rocks.     

Modeling discontinuous rock mass based on smoothed finite element method

This paper aims at developing a method for modeling rock mass with preexisting multiple discontinuities within the framework of the smoothed finite element method (SFEM). The discontinuity is simulated by an interface element with zero thickness, the stiffness matrix of which are derived explicitly based on the SFEM. An elastic damage constitutive relation with residual strength is introduced in order to describe the nonlinear mechanical behavior of the discontinuities. The computation codes of the present method were developed. The present method has been verified to be a sound approach for modeling discontinuous rock mass, inheriting the advantages of the SFEM.     

基于岩石图像深度学习的岩性自动识别与分类方法

岩石岩性的识别与分类对于地质分析极为重要,采用机器学习的方法建立识别模型进行自动分类是一条新的途径。基于Inception-v3深度卷积神经网络模型,建立了岩石图像集分析的深度学习迁移模型,运用迁移学习方法实现了岩石岩性的自动识别与分类。采用此方法对所采集的173张花岗岩图像、152张千枚岩图像和246张角砾岩图像进行了学习和识别分类研究,通过训练学习建立岩石图像深度学习迁移模型,并分别采用训练集和测试集中的岩石图像对模型进行了检验分析。对于训练集中的岩石图像,每组岩石分别用3张图像测试,三种岩石的岩性分类均正确,且分类概率值均达到90%以上,显示了模型良好的鲁棒性;对于测试集中的岩石图像,每组岩石分别采用9张图像进行识别分析,三种岩石的岩性分类均正确,并且千枚岩组图像分类概率均高于90%,但是花岗岩组2张图像和角砾岩组的1张图像分类概率值不足70%,概率值较其他岩石图像低,推测其原因是训练集中相同模式的岩石图像较少,导致模型的泛化能力减小。为了提高识别精确度,对准确率较低的岩石图像进行截取,分别取其中的3张图像加入训练集进行再训练,增加与测试图像具有相同模式的训练样本;在新的模型中,对3张图像进行二次检验,测试概率值均达到85%以上,说明在数据足够的状况下模型具有良好的学习能力。与传统的机器学习方法相比,所提出的岩石图像深度学习方法具有以下优点:第一,模型通过搜索图像像素点提取物体特征,不需要手动提取待分类物体特征;第二,对于图像像素大小,成像距离及光照要求低;第三,采用适当的训练集可获得较好的识别分类效果,并具有良好鲁棒性和泛化能力。     

鲜水河断裂带乾宁段岩石特征与内部结构及其物理-化学性质

断裂带的变形行为和断层滑移机制是目前地震研究关注的热点,断裂带岩石特征、内部结构与物理化学性质是确定断层蠕滑或粘滑行为以及断层滑移机制的基础和关键。本文以鲜水河断裂带乾宁段地表出露的断裂岩为研究对象,通过野外地质调查、室内光学显微镜、扫描电镜、粒度统计、粉末X射线衍射分析（XRD）和薄片X射线荧光光谱分析（XRF）等多种研究方法,对鲜水河断裂带岩石特征、结构构造、物性、矿物成分及化学元素分布开展了详细的分析,并探讨了相关变形行为和滑移机制。分析表明：（1）断裂带核部主要由黑色断层泥、浅黄色及黄色断层角砾岩和灰色碎裂岩、灰色断层角砾岩组成,呈单核对称结构;（2）黑色断层泥厚3~5cm,具有快速滑动结构特征,表现为断层粘滑行为。断层泥可划分出13个滑移带,最窄滑移带厚约40μm,至少代表13期古地震事件;（3）断层泥主要由伊利石、高岭石和石英等矿物组成,其中边部伊蒙混层含量异常高,为最新一次古地震的主滑移带。由于伊利石和伊蒙混层（或蒙脱石）为主要黏土矿物的断层泥渗透率低、孔隙流体压力大,以及发现断层泥楔入脉,表明地震过程中断层滑动存在热增压弱化机制;（4）从断层泥不同滑动带中碎块蚀变程度和矿物分布特征来看,地震主滑动带有向碎裂岩方向迁移的趋势。推测断层在滑动过程中,更趋向于向弱矿物含量高（如伊利石、伊蒙混层）、强矿物含量低（如方解石、高岭石）的围岩一侧迁移。     

金沙江下虎跳滑石板滑坡破坏机理及稳定性分析

1996年10月28日发生在金沙江下虎跳峡口的滑石板滑坡,因堵塞金沙江、摧毁公路而引起高度重视。野外调查发现,滑石板滑坡是典型结构面控制的顺层基岩滑坡,层面和两组近于垂直层面的节理将边坡厚层灰岩切割成不同规模的菱形块体,长期的风化作用、重力、降雨、人工开挖公路等因素作用下,这些菱形块体从顺层稳定状态演化为稳定性较差的悬挂体,受河谷边坡斜向割切影响,表现为梯级叠次悬挂,最外边的悬挂体稳定性最差,最容易发生快速滑动。因此,该滑坡每次发生破坏滑动的体积不会很大,但是具有多次重复发生的特点,值得高度注意。     

喜马拉雅山东南地区地质灾害发育规律初步研究

利用遥感手段,结合MapGis,研究了喜马拉雅山东南地区地质灾害的发育情况,发现本区发育的主要地质灾害有滑坡、崩塌、泥石流、冰湖以及堰塞湖。其中崩塌、滑坡、泥石流斜坡地质灾害是本区最重要的地质灾害类型,占到总灾害数量的95.3%。在此基础上对喜马拉雅山东南地区地质灾害发育规律初步研究,发现本区地质灾害的发育在空间上的分布并非均匀,而是具有丛集性的特点。滑坡灾害主要发育在隆子和朗县。泥石流灾害比较严重的有米林、隆子和洛扎3县,而崩塌则主要集中在隆子县。研究发现,本区滑坡发育与地层、地形坡度以及土地类型关系密切,其中修康群、日当组和念青唐古拉群是本区的易滑地层。涅如组由于面积大,其中发育的滑坡较多,但是滑坡的发育率只略高于本区的平均水平。统计表明,16～30的坡度范围是滑坡最容易发生的。大于45以上的坡段很少发生滑坡。灌木林和天然草地这两种土地类型滑坡发育率最高。对于泥石流,研究表明,涅如组中泥石流发育面积最大,发育率也最高。泥石流发育的最适宜坡度也是16～30这样一个坡度范围。冰川和永久积雪区则最易发生泥石流。崩塌发育与地层类型、坡度的关系较为密切,崩塌主要发育在涅如组中,并且集中在坡度大于60以上的陡坡段中。这些初步成果的取得,是以后进行该区地质灾害空间预测的基础。     

基于虚功原理的悬索桥隧道式锚碇边坡稳定性分析

本文基于虚功原理,结合工程地质勘察资料,在岩层层面和不利结构面组合切割下,由于锚碇工程荷载作用,研究坝陵河大桥西岸隧道式锚碇边坡的整体稳定性。本文的力学分析模型采用多个块体破坏机构的分析方法,假定块体的底面和侧面均达到了极限平衡状态,安全系数可通过功能平衡的虚功原理表达式来求解。功能平衡方程中仅有安全系数一个未知数,条块底部和界面上的法向压力和切向摩擦力并不出现,求解大为简化。