A Reservoir–Geomechanical Model to Study the Likelihood of Tensile and Shear Failure in the Caprock of Weyburn CCS Project with Regard to Interpretation of Microseismic Data

According to poroelastic theory and also field observations, an increase in reservoir pore pressures can result in a decrease in horizontal stresses in the seal layers. This reduction is in favor of hydrofracture initiation and reactivation of weak planes and has to be studied in caprock integrity analyses. In this paper, a field scale reservoir–geomechanical (GEM-FLAC3D) model is developed for the Phase IB area of the Weyburn (Canada) CCS project that is located in Williston sedimentary basin. A one-way coupling has been conducted between the two codes for the period of CO₂ injection in Phase IB area from 2000 to 2010. Therefore, the reservoir pore pressures at selected timesteps are unidirectionally fed to the FLAC3D. In order to study the likelihood of tensile and shear failure in the seal layer on top as a result of stress transfer due to poroelastic effects, two margin ratios are defined for tensile and shear failure and their variations are studied above the reservoir upon changes in pore pressures within the reservoir. The results show that overall; the likelihood of shear failure has been about 25% greater than that of tensile failure in Weyburn. However, between 2008 and 2009, the pressures were high enough to trigger both tensile and shear mechanisms above the reservoir. A discussion is also presented on relevance of this study for interpretation of microseismic data with regard to failure margin ratios, magnitudes and distribution of events recorded in each year.     

节理化岩体复合式滑移破坏的模型试验研究

针对节理化岩体的自身特点,提出沿岩体层面产生的应力集中以及层面与节理之间的残余完整岩桥破坏,将可能导致岩体产生复合式滑移破坏。并在室内利用岩石三轴刚性伺服机,通过模型试验对岩体的层面和一组与层面正交的节理的情况进行模拟研究,分析了结构面性状和围压对复合式滑移破坏的影响,给出了在贯通性结构面(层面)倾角、非贯通性结构面以及围压等因素共同作用下,节理化岩体发生复合式滑移破坏的条件以及破坏面特征。在与前人的试验资料进行对比分析的基础上得到可靠的结论,期望为在将来的实践工作中正确评估节理化岩体的破坏模式提供指导性意见。      

节理化岩体复合式滑移破坏的模型试验研究

针对节理化岩体的自身特点,提出沿岩体层面产生的应力集中以及层面与节理之间的残余完整岩桥破坏,将可能导致岩体产生复合式滑移破坏。在室内利用岩石三轴刚性伺服机,通过模型试验对岩体的层面和一组与层面正交的节理的情况进行模拟研究,分析了结构面性状和围压对复合式滑移破坏的影响,给出了在贯通性结构面（层面）倾角、非贯通性结构面的以及围压等因素共同作用下,节理化岩体发生复合式滑移破坏的条件以及破坏面特征。最后在与前人的试验资料进行对比分析的基础上得到可靠的结论,期望为在将来的实践工作中正确评估节理化岩体的破坏模式提供指导性意见。     

The Impact of Heterogeneity on the Anisotropic Strength of an Outcrop Shale

Properly accounting for the mechanical anisotropy of shales can be critical for successful drilling of high inclination wells, because shales are known to be weak along bedding planes. To optimize the drilling parameters in such cases, a sufficiently representative, anisotropic rock mechanical model is therefore required. This paper presents such a model developed to better match results from a dedicated, extensive set of uniaxial and triaxial compression tests performed on plugs of Mancos outcrop shale with different orientations relative to the bedding plane. Post failure inspection of the plugs shows that the failure planes are to some extent affected by the orientation of the applied stress relative to the bedding planes, indicating that the bedding planes may represent weak planes which tend to fail before intrinsic failure occurs, whenever the orientation of these planes is suitable. The simple “plane of weakness” model is commonly used to predict strength as function of orientation for such a rock. A comparison of this model to the experimental data shows, however, that the weak planes seem to have an impact on strength even outside the range of orientations where the model predicts such impact. An extension of this model allowing the weak planes to be heterogeneous in terms of patchy weakness was therefore developed. In this model, local shear sliding may occur prior to macroscopic failure, leading to enhanced local stresses and corresponding reduction in strength. The model is found to give better match with strength data at intermediate orientations. The model is also able to partly predict the qualitatively different variation of Young’s modulus with orientation for this data set.     

水下岩质边坡在波浪作用下的稳定性分析

运用ANSYS程序对水下的岩质边坡进行数值模拟分析,采用不同的材料强度来模拟岩石和软弱结构面。应用强度折减方法对边坡的强度参数进行折减,直至边坡出现破坏,分析边坡在达到破坏时内部的应力应变情况和位移变形情况。对于波浪力,本文采用简化的波压力方程模拟低频的水波对边坡的波压力影响,并用多个模型比较计算的方法得出最佳的模型尺寸和水下荷载的计算模式。     

含采空区层状岩体微震信号传播规律试验研究

微震信号在层状岩体中传播规律的研究对于实现震源准确定位具有重要的现实意义。在天然岩体中,微震信号在岩层中传播的速度,不仅受其内部因素即岩石自身物理性质的影响,还受外部因素如地质结构面、断层及采空区等条件的影响。通过室内试验的方法对微震信号在不同岩层及采空区处的传播规律进行研究。结果表明：波速随着传播距离的增加而逐渐衰减,且岩层密度越大衰减越慢;穿过的结构面数量越多,微震信号传播速度衰减比例越大;能量的衰减与波速的衰减一致。采空区断面横截面积越大,微震信号通过的用时越长;采空区周围岩石密度、弹性模量等物理性质的数值越大,微震信号通过时衰减越少。边际谱分析结果表明：断层对频率高的微震信号阻碍作用更强。     

The Dynamic Evaluation of Rock Slope Stability Considering the Effects of Microseismic Damage

A state-of-the-art microseismic monitoring system has been implemented at the left bank slope of the Jinping first stage hydropower station since June 2009. The main objectives are to ensure slope safety under continuous excavation at the left slope, and, very recently, the safety of the concrete arch dam. The safety of the excavated slope is investigated through the development of fast and accurate real-time event location techniques aimed at assessing the evolution and migration of the seismic activity, as well as through the development of prediction capabilities for rock slope instability. Myriads of seismic events at the slope have been recorded by the microseismic monitoring system. Regions of damaged rock mass have been identified and delineated on the basis of the tempo-spatial distribution analysis of microseismic activity during the periods of excavation and consolidation grouting. However, how to effectively utilize the abundant microseismic data in order to quantify the stability of the slope remains a challenge. In this paper, a rock mass damage evolutional model based on microseismic data is proposed, combined with a 3D finite element method (FEM) model for feedback analysis of the left bank slope stability. The model elements with microseismic damage are interrogated and the deteriorated mechanical parameters determined accordingly. The relationship between microseismic activities induced by rock mass damage during slope instability, strength degradation, and dynamic instability of the slope are explored, and the slope stability is quantitatively evaluated. The results indicate that a constitutive relation considering microseismic damage is concordant with the simulation results and the influence of rock mass damage can be allowed for its feedback analysis of 3D slope stability. In addition, the safety coefficient of the rock slope considering microseismic damage is reduced by a value of 0.11, in comparison to the virgin rock slope model. Our results demonstrate that microseismic activity induced by construction disturbance only slightly affects the stability of the slope. The proposed feedback analysis technique provides a novel method for dynamically assessing rock slope stability and can be used to assess the slope stability of other similar rock slopes.     

Microseism Induced by Transient Release of In Situ Stress During Deep Rock Mass Excavation by Blasting

During deep rock mass excavation with the method of drill and blast, accompanying the secession of rock fragments and the formation of a new free surface, in situ stress on this boundary is suddenly released within several milliseconds, which is termed the transient release of in situ stress. In this process, enormous strain energy around the excavation face is instantly released in the form of kinetic energy and it inevitably induces microseismic events in surrounding rock masses. Thus, blasting excavation-induced microseismic vibrations in high-stress rock masses are attributed to the combined action of explosion and the transient release of in situ stress. The intensity of stress release-induced microseisms, which depends mainly on the magnitude of the in situ stress and the dimension of the excavation face, is comparable to that of explosion-induced vibrations. With the methods of time–energy density analysis, amplitude spectrum analysis, and finite impulse response (FIR) digital filter, microseismic vibrations induced by the transient release of in situ stress were identified and separated from recorded microseismic signals during a blast of deep rock masses in the Pubugou Hydropower Station. The results show that the low-frequency component in the microseismic records results mainly from the transient release of in situ stress, while the high-frequency component originates primarily from explosion. In addition, a numerical simulation was conducted to demonstrate the occurrence of microseismic events by the transient release of in situ stress, and the results seem to have confirmed fairly well the separated vibrations from microseismic records.     

基于结构面性状对岩体力学性质响应的模拟研究

为研究结构面性状对锦屏大奔流料场边坡岩体性质的影响,在对该边坡砂岩进行室内三轴压缩试验的基础上,开展含4组不同结构面性状岩体的力学模拟测试,最后对该边坡岩体的开挖进行数值模拟。试验结果表明,岩样破坏过程可以分为4个典型的阶段: 弹性变形阶段,屈服阶段,残余强度阶段及塑性流动阶段。随着围压的增大,岩石破坏所需偏应力逐渐增大,且岩石逐渐表现出延性特征。结构面性状模拟结果表明,含软弱充填结构面比不含者对岩体力学性质影响(弱化)大,有倾角结构面比无倾角结构面对岩体性质影响较大,大倾角结构面总体上比小倾角的影响大,软弱充填厚度对岩体内摩擦角影响较大。实例模拟结果表明,开挖卸荷使得边坡后缘岩体产生卸荷回弹变形,产生X向拉伸位移,煌斑岩脉的存在使得塑性区域沿此软弱带发展,开挖第8步时,塑性区贯通,边坡发生失稳。研究结果对陡高边坡岩体开挖的稳定性评价具有一定指导意义。     

A three-dimensional equivalent continuum constitutive model for jointed rock masses containing up to three random joint sets

A three dimensional constitutive model is formulated for deformation analysis of jointed rock masses containing up to three joint sets with arbitrary spatial configurations. A representative elementary volume (REV) that represents the deformational response of the rock mass is defined and the constitutive relationships are developed based on superposition of deformations of the REV components. By representing the constitutive relationships in a tensorial form, the model is able to implement deformation anisotropy of jointed rock masses. The Mohr-Coulomb failure criterion with tension cut-off is used for the intact rock and the joint sets. The model is implemented in FLAC^3D and the deformations and strength values calculated by the model are compared with the results from a 3DEC model and analytical solutions. The model results are in good agreement with those obtained from 3DEC.     

Critical stress difference and orientation of faults in rocks with strength anisotropies: the two-dimensional case

Equations are derived for the critical stress difference on thrust, normal and strike-slip faults with finite cohesive strength, both in homogeneous, isotropic rock, and along pre-existing strength anisotropies with different cohesion and coefficient of friction, subject to the limitation that the plane of anisotropy contains the intermediate axis of stress. The range of orientations for which sliding occurs along pre-existing planes of weakness rather than along a new fault is given as a function of material parameters and critical stress difference for the intact rock. Given the principal stress directions, the equations allow the direct calculation of the stress and orientation conditions for the three faulting regimes.     

大岗山水电站厂房断层控制区域微震监测分析

大岗山水电站地下厂房断层发育,两条陡倾角断层穿过厂房中后部,控制着厂房的围岩稳定性。拱顶处发育的断层由于开挖而出现塌方,需要对塌空区以及地下厂房进行稳定性监测。2010年10月建立南非综合微震系统（ISS）微震监测系统,对断层控制的厂房区域进行微震监测。采用能级分布法、时空分布法、能量指数法、位错法等对断层控制下的塌空区、厂房边墙、断层周边岩体进行了稳定性评价。监测表明,塌空区微震事件少,能级小;边墙岩体位移较小,微震事件数目较少;断层面间滑动位移分布不均,整体位移量不大,断层面岩体活动性不强。因此,断层区域在进行加固后,岩体整体稳定性有所提高,受断层影响较小,故地下厂房稳定性较好。微震监测对水电厂房开挖及运行期的安全稳定性评价有较强的指导意义。     

乌东德水电站地下厂房层状岩体稳定性及变形机制

乌东德水电站右岸地下厂房洞室群密集、开挖规模大、地质条件复杂,开挖期间变形破坏问题突出。在地应力实测资料、工程地质特性、开挖过程分析的基础上,利用离散元程序对右岸地下厂房开挖过程进行数值模拟,再现不同模型的应力场、位移场演化过程,分析陡倾角层面对于主厂房变形的影响作用。同时,构建高精度微震监测系统,对开挖诱发的微震活动实时在线监测;分析微震事件聚集规律、S波及P波能量比值Es /Ep分布特征,识别和圈定潜在风险区域,揭示围岩变形机制。微震监测与数值模拟结果表明：在开挖卸荷作用下,陡倾角层面附近岩石破裂,微震事件聚集,能量释放,引起围岩变形破坏。围岩以结构面控制的应力驱动型张拉破坏为主,并伴随少量重力驱动型变形。该研究结果为乌东德水电站地下厂房的开挖和支护设计方案提供重要依据。     

粉砂质板岩力学特性及各向异性特性

为研究鄂西北地区板岩的力学特性及各向异性特征,开展了志留系粉砂质板岩的单轴压缩和巴西劈裂试验,分析了试样的力学特性和各向异性特征以及在不同荷载作用下的变形破坏模式,揭示了不同破坏模式的力学机制,并通过数值分析研究了不同板理面角度板岩边坡破坏模式及力学机制。研究结果表明,粉砂质板岩中的板理面是影响岩体力学行为的弱面,使得粉砂质板岩表现出明显的各向异性的特征。在单轴压缩条件下,粉砂质板岩在垂直板理面方向比平行板理面方向更易变形,变形量更大;平行板理面方向加载时试样破坏类型为竖向劈裂型张拉破坏,其实质是压杆失稳;垂直板理面加载时试样的破坏类型为穿切板理面的劈裂型剪切破坏;所测得的力学参数各向异性也较明显。在劈裂荷载作用下,粉砂质板岩的破坏模式主要有张拉劈裂破坏和沿板理面剪切破坏两种,所得抗拉强度在平行板理面方向上最大,在垂直板理面方向上最小,两个方向上的抗拉强度均小于抗压强度。由于板理面间的抗拉强度极低,在受到与板理面呈小角度相交的劈裂荷载作用时,容易产生沿板理面的张拉劈裂或拉剪破坏,在实际工程中应尽量避免板理面间的受拉破坏和沿板理面的拉剪破坏。在边坡工程中,板理面倾向、倾角对边坡破坏模式及力学机制有较大影响,边坡防护治理需充分考虑这一影响。上述研究为粉砂质板岩区岩质边坡防护治理以及其他岩石工程设计与施工提供理论依据与技术基础。     

沉积结构面及其对岩体力学性质的影响

基于沉积结构面的成因类型,通过理论分析和实验研究,分析了沉积结构面的变形与破坏特征,研究了沉积结构面力学特征及其力学效应,提出了弱面型层理、层系或层系组界面和岩层面以及不整合面均为沉积构造弱面,这些弱面对沉积岩体力学性质产生重要的影响。分析了沉积岩体沿层面方向和垂直层面方向岩体力学性质的差异性。研究表明,垂直层面方向的纵波速度和弹性模量要比平行层面方向低;平行层面方向的抗拉强度要大于垂直层面方向;而平行层面方向的抗压强度与凝聚力小于垂直层面方向。沉积岩体易于产生垂直于结构面方向的张性破坏和沿结构面方向的剪切破坏,在采动影响下顶板岩体易于沿岩层层面产生离层破坏现象,这进一步揭示了沉积结构面变形与破坏机制。      

层理面特性对砂岩断裂力学行为的影响研究

层理弱面对层状岩石的力学特性影响较显著，为了研究层理面特性对岩石断裂力学特性的影响，开展了具有不同层理方向的砂岩试样三点弯试验，探讨了砂岩断裂韧度及断裂模式的各向异性。之后基于有限元中的黏聚单元建立了数值模型，采用数值模拟方法研究了层理面强度对各层理角度试样断裂力学行为的影响规律。结果表明：层理方向影响下砂岩的断裂韧度及模式存在各向异性；同一层理方向试样的断裂韧度随层理面强度的增大而增大，且试样的层理面与加载方向夹角越小，断裂韧度受层理面强度变化影响越明显；试样的断裂模式不仅与层理面强度有关，还受层理倾角的控制，层理面与加载方向夹角θ = 0o试样断裂模式基本不受层理面强度影响，θ = 30o试样主要沿层理面张拉或剪切破坏，且沿层理面的破裂长度随层理面强度的降低逐渐增大；层理面强度较大时，θ = 45o试样主要沿层理面张拉破坏，θ = 60o～90o试样主要以贯穿层理的张拉破坏为主；层理面强度较小时，θ = 45o～90o试样均以沿层理面的剪切破坏为主，其中θ = 45o试样沿层理剪切长度最大。另外，通过数值模拟结果分析了层理面强度及方向对试样的起裂角及裂纹扩展路径产生的影响。该研究成果可作为层状岩石断裂力学理论的有益补充。     

Nanoseismic monitoring of gravity-induced slope instabilities for the risk management of an aqueduct infrastructure in Central Apennines (Italy)

A monitoring system is operative in the Peschiera Springs slope (Central Apennines, Italy) to mitigate the landslide risk related to the hosted main drainage plant of Rome aqueducts by providing alert warning. Such a strategy allows to avoid out-of-service episodes so reducing extra-costs of water distribution management. The Peschiera Springs slope is involved in a rock mass creep characterized by an average steady strain rate of 1 mm year⁻¹ and responsible for several landforms including sinkholes, subvertical scarps and trenches. Moreover, an average aquifer discharge of 19 m³ s⁻¹ causes an intense limestone dissolution concentrated in correspondence with release bands and discontinuities that dislodge the jointed rock mass. Since 2008, an accelerometric network has been operating within the slope; about 1300 microseismic local events were recorded up to now, distinguished in failures and collapses. A control index, based on frequency of occurrence and cumulative energy of the recorded microseismic events was defined to provide three levels of alert. In 2013, a temporary nanoseismic Seismic Navigation System (SNS) array was installed inside a tunnel of the drainage plant to integrate the pre-existent seismic monitoring system. This array allowed to record 37 microseismic events, which locations are in good agreement with the evolutionary geological model of the ongoing gravitational slope deformation. In 2014, a permanent nanoseismic SNS array was installed in the plant and allowed to record several sequences of underground collapses including more than 500 events. The nanoseismic monitoring system is allowing to: (1) increase the detection level of the monitoring system; (2) locate hypocentres of the events; and (3) detect precursors of the strongest events.

     

Nanoseismic monitoring of gravity-induced slope instabilities for the risk management of an aqueduct infrastructure in Central Apennines (Italy)

A monitoring system is operative in the Peschiera Springs slope (Central Apennines, Italy) to mitigate the landslide risk related to the hosted main drainage plant of Rome aqueducts by providing alert warning. Such a strategy allows to avoid out-of-service episodes so reducing extra-costs of water distribution management. The Peschiera Springs slope is involved in a rock mass creep characterized by an average steady strain rate of 1 mm year^?1 and responsible for several landforms including sinkholes, subvertical scarps and trenches. Moreover, an average aquifer discharge of 19 m³ s^?1 causes an intense limestone dissolution concentrated in correspondence with release bands and discontinuities that dislodge the jointed rock mass. Since 2008, an accelerometric network has been operating within the slope; about 1300 microseismic local events were recorded up to now, distinguished in failures and collapses. A control index, based on frequency of occurrence and cumulative energy of the recorded microseismic events was defined to provide three levels of alert. In 2013, a temporary nanoseismic Seismic Navigation System (SNS) array was installed inside a tunnel of the drainage plant to integrate the pre-existent seismic monitoring system. This array allowed to record 37 microseismic events, which locations are in good agreement with the evolutionary geological model of the ongoing gravitational slope deformation. In 2014, a permanent nanoseismic SNS array was installed in the plant and allowed to record several sequences of underground collapses including more than 500 events. The nanoseismic monitoring system is allowing to: (1) increase the detection level of the monitoring system; (2) locate hypocentres of the events; and (3) detect precursors of the strongest events.     

结构面对岩质边坡地震动影响的数值模拟研究

[摘 要] 本文采用数值模拟的方法,利用UDEC 软件,开展了结构面对岩质边坡地震动影响的数值 模拟研究,研究内容涉及结构面产状和位置变化对岩体边坡地震动放大系数空间展布的影响。根据岩 体边坡中结构面发育的一般规律,本文模型工况主要考虑与坡面相交的贯穿性结构面,模拟的主要物理 对象是层面、贯穿性长大软弱结构面等。数值模拟结果表明:结构面越靠近坡顶,边坡上部的地震动响 应越强烈;顺倾结构面会使坡肩部形成更强烈的地震动响应;边坡越陡,坡顶和坡肩部的地震动响应越 强烈。上述规律对于岩质边坡地震稳定性分析具有指导意义。     

天福庙拱坝坝肩抗滑稳定性研究

本文介绍了天福庙拱坝坝肩岩体稳定性研究,并提出了相应加固措施的建议。即在详尽的工程地质水文地质调查的基础上分析出各种可能的失稳机制,并组合出各种滑移边界条件;通过现场和室内岩体力学试验,测出构成边界条件的各种地质结构面的抗剪强度和变形常数;通过三维有限元弹性分析计算,得到应力应变和破损区的分布,检验了失稳机制和边界条件的合理性;通过刚体极限平衡法,计算出各种荷载和边界条件组合情况下的抗滑稳定安全系数;尔后针对最不稳定边界条件及其控制稳定性的主要因素,提出有效的加固措施的建议。文中还讨论了岩体和各种软弱结构面抗剪强度的取值原则、试验值的使用原则以及非均质结构面和组合结构面的综合抗剪强度问题。