Computational model coupling mode II discrete fracture propagation with continuum damage zone evolution

We propose a numerical method that couples a cohesive zone model (CZM) and a finite element‐based continuum damage mechanics (CDM) model. The CZM represents a mode II macro‐fracture, and CDM finite elements (FE) represent the damage zone of the CZM. The coupled CZM/CDM model can capture the flow of energy that takes place between the bulk material that forms the matrix and the macroscopic fracture surfaces. The CDM model, which does not account for micro‐crack interaction, is calibrated against triaxial compression tests performed on Bakken shale, so as to reproduce the stress/strain curve before the failure peak. Based on a comparison with Kachanov's micro‐mechanical model, we confirm that the critical micro‐crack density value equal to 0.3 reflects the point at which crack interaction cannot be neglected. The CZM is assigned a pure mode II cohesive law that accounts for the dependence of the shear strength and energy release rate on confining pressure. The cohesive shear strength of the CZM is calibrated by calculating the shear stress necessary to reach a CDM damage of 0.3 during a direct shear test. We find that the shear cohesive strength of the CZM depends linearly on the confining pressure. Triaxial compression tests are simulated, in which the shale sample is modeled as an FE CDM continuum that contains a predefined thin cohesive zone representing the idealized shear fracture plane. The shear energy release rate of the CZM is fitted in order to match to the post‐peak stress/strain curves obtained during experimental tests performed on Bakken shale. We find that the energy release rate depends linearly on the shear cohesive strength. We then use the calibrated shale rheology to simulate the propagation of a meter‐scale mode II fracture. Under low confining pressure, the macroscopic crack (CZM) and its damaged zone (CDM) propagate simultaneously (i.e., during the same loading increments). Under high confining pressure, the fracture propagates in slip‐friction, that is, the debonding of the cohesive zone alternates with the propagation of continuum damage. The computational method is applicable to a range of geological injection problems including hydraulic fracturing and fluid storage and should be further enhanced by the addition of mode I and mixed mode (I+II+III) propagation. Copyright © 2016 John Wiley & Sons, Ltd.     

地震工况下坝坡的抗滑安全系数取值标准研究

地震工况下的坝坡抗滑稳定性核算一直是大坝安全性评价的重要环节,采用科学合理的计算方法是评价坝坡安全性的基本保证。本文从概率极限分析的角度探讨了坝坡稳定性可靠度分析的概率意义,提出了将地震发生当做概率事件的坝坡可靠度计算方法。为验证将地震超越概率纳入坝坡稳定分析中的合理性,本文通过两个典型算例和小浪底坝坡典型剖面的抗滑稳定性分析研究了坝坡在地震工况下的单一安全系数与可靠指标的变化规律。研究结果表明在坝坡稳定分析中引入超越概率计算所得结果与工程实际相符,与以往的将地震发生当做确定性事件计算方法相比,该方法显著提高了坝坡的可靠指标值;比较可靠指标与单一安全系数的取值标准显示出在坝坡稳定分析中合理地引入地震概率分析方法能够保证坝坡在地震和正常工况下有相近的可靠指标值     

Fracturing around a preconditioned deep level gold mine stope

The mining faces of the highly stressed tabular stopes in South Africa's deep level gold mines are prone to a type of rock burst termed a face-burst. As a means of ameliorating these face-burst conditions, a destressing technique termed preconditioning has been employed. Part of the studies into the quantification of the effects and mechanisms of preconditioning was a detailed investigation of the fracture pattern around unpreconditioned and preconditioned stopes. Techniques included mapping of the fractures exposed in the mined areas of the stope, measurement of profiles of the stope hangingwall and the use of Ground Penetrating Radar (GPR). These studies indicate that preconditioning does not produce new sets of fractures and stress redistribution ahead of the face occurs by re-activation of specific pre-existing fractures, thereby reducing the potential for face-bursting.     

泥浆护壁灌注桩桩底后压浆技术及其单桩竖向承载力的估算方法

通过国内近年来桩底后压浆的工程实例,介绍了泥浆护壁灌注桩后压浆技术的发展,分析了桩底后压浆技术适应的地质条件和有关技术参数,归纳出单村竖向承载力的估算方法,供初步设计时作为布桩依据。     

桩端加环对PHC桩影响的试验研究

PHC桩在动力打桩过程中,往往会因为锤击数或锤击能量过高而发生桩头或桩身损坏,引起的挤土效应和振动效应也会对周围环境造成不良影响.PHC桩桩端加环可减少沉桩过程中桩侧土体摩阻力的影响,结合某电厂工程综合试桩工作,通过对沉桩总锤击数、最终贯入度、土塞高度、打桩振动测试及高应变检测结果的分析,指出桩端加环对保证PHC桩的完...     

Hybrid finite-discrete element modelling of dynamic fracture and resultant fragment casting and muck-piling by rock blast

A state-of-the-art review is conducted to highlight the fracture mechanism in rock blast and advantages and limitations of various methods in modelling it. A hybrid finite-discrete element method (FEM-DEM) is implemented to simulate rock fracture and resultant fragment muck-piling in various blasting scenarios. The modelled crushed, cracked and long radial crack zones are compared with those in literatures to calibrate the hybrid FEM-DEM. Moreover, the hybrid modelling reproduces the rock fragmentation process during blasting. It is concluded that the hybrid FEM-DEM is superior to continuous and discontinuous methods in terms of modelling dynamic fracture of rock under blast-induced impact load.     

Modeling hydraulic fracture propagation using cohesive zone model equipped with frictional contact capability

We present a new pore pressure cohesive element for modeling the propagation of hydraulically induced fracture. The Park-Paulino-Roesler cohesive zone model has been employed to characterize the fracturing behavior. Coulomb’s frictional contact model has been incorporated into the element to model the possible shear reactivation of pre-existing natural fractures. The developed element has been validated through a series of single-element tests and an available analytical solution. Furthermore, intersection behaviors between the hydraulic fracture and the natural fracture under various conditions have been predicted using the present element, which shows good agreement with experimental results.     

低渗透油藏菱形反九点井网不等缝长优化研究

在低渗透油藏开发中广泛使用菱形反九点井网,但针对不同位置生产井及注水井的裂缝缝长配比关系还缺少广泛深入的研究。基于长庆油田某低渗透油藏参数,建立了数值模拟模型,在菱形反九点井网中分别设置不同的边井、角井和注水井裂缝穿透比,研究不同裂缝配置模式对水驱开发效果的影响。流线模拟及数值模拟计算结果表明,在注水开发中,注水井缝长是影响油田产液量及水驱效果的关键因素;边井裂缝较角井裂缝更易见水,对油田综合含水影响较大,在菱形反九点井网中需要进行不等缝长优化设计;边井裂缝缝长对水驱效果的影响分为早、中、晚期(Ⅰ、Ⅱ、Ⅲ)曲线段,随着注水井裂缝穿透比增加,早期(Ⅰ)曲线段增长,中期(Ⅱ)曲线段缩短,而晚期(Ⅲ)曲线段变化较小;随着角井裂缝穿透比增加,采出程度逐渐增大,当角井裂缝穿透比较大时,存在“反转”现象;随着注水井裂缝穿透比增加,角井裂缝优化穿透比减小,而边井裂缝优化穿透比增加。该研究成果对于指导低渗透油藏菱形反九点井网裂缝优化设计,提高波及系数和开发效果具有重要的指导意义。     

含有硬包裹体分布的非均质岩石水力压裂特性研究

基于完全流-固耦合的弹塑性理论给出了水力压裂数值计算的弥散裂缝模型,其中材料的弹性部分采用线弹性本构关系,塑性部分采用摩尔-库仑破坏准则及强化准则。依据当前的有效应力状态修正渗透系数来模拟压裂液在裂缝中的流动。渗透系数的修改使用双曲正切函数,并采用平均有效应力作为水力裂缝的起裂判据。在ABAQUS软件中通过用户自定义程序添加了该模型。根据岩石的切面照片建立了含有硬包裹体分布的非均质岩石的有限元计算模型,模拟了中心点注水条件下的水力压裂传播过程,讨论了在常应力状态下非均质岩石中开裂区域、典型位置的应力路径变化和裂缝传播范围随时间变化的特点。进行了多种条件下含有硬包裹体分布的岩石材料的数值试验,得出了基岩材料的弹性模量、凝聚力和渗透系数以及注水速率对峰值注水压力、平均注水压力和裂缝开度的影响规律。     

岩体破裂尺度与频率特征关系及其工程实证研究

将S变换时频分析技术引入微震信号波形分析领域,通过震源破裂尺度,研究微震信号频率特征与岩体变形破坏之间的联系。以西南地区3个水电工程作为研究对象,将大量微震监测数据作为统计样本以减小外部因素对参数计算的影响。在该基础上统计、分析岩体破裂尺度与微震信号频率之间的关系,结合大渡河猴子岩水电站地下厂房微震与多点位移计监测资料进行工程实证研究。研究结果表明,当岩体破裂尺度增大时微震信号高频成分减少,表现出低频特征;在围岩变形过程中微震信号高频成分呈现出先减少再增加的演化规律,围岩发生明显变形前微震信号出现由高频向低频转移的现象,因此,可将聚集区域内微震信号频率降低作为该区域岩体宏观变形破坏的前兆信息。研究结果可为基于微震监测的岩体工程稳定性分析及其灾害预测预报提供参考。