基于岩体爆破累积损伤效应的Hoek-Brown准则修正公式

爆破开挖作业引起的扰动与损伤,特别是频繁爆破产生的累积损伤效应,必然导致岩体完整性降低,岩体力学参数弱化,从而威胁岩体工程稳定性。考虑岩体爆破损伤及其累积效应,指出了Hoek-Brown(赫克-布朗)准则及其改进公式的不足,引入完整性系数K v和损伤因子D,建立了可以表征岩体爆破累积损伤效应、岩体爆破扰动状态及其力学参数弱化程度的mb和s的取值方法。基于声速变化与爆破累积损伤效应之间的联系,建立了以岩体声速降低率η为基准量的岩体爆破累积损伤扩展模型,并将该模型成功应用于提出的Hoek-Brown准则修正公式中。根据现场模拟爆破试验与声波测试数据,对提出的修正公式进行了分析,结果表明,建立的考虑岩体爆破累积损伤效应的Hoek-Brown准则修正公式是合理的。     

采用不同Drucker-Prager屈服准则得到的边坡安全系数的转换

研究了边坡稳定分析的有限元强度折减过程中（平面应变条件下）Drucker-Prager（简称D-P）屈服面的变化特点,从而提出了有限元强度折减法采用不同屈服准则计算得到的边坡安全系数之间的转换关系式。在利用有限元强度折减法计算边坡稳定安全系数时,ANSYS有限元软件采用的岩土材料屈服准则为Mohr-Coulomb（简称M-C）六边形外接圆D-P准则,可先求出外接圆D-P准则条件下的安全系数,然后利用所提出的安全系数转换公式就可直接计算出各D-P准则条件下的安全系数。因此,通过转换就可以在ANSYS程序中实现不同M-C准则,而不需要进行二次开发。采用ANSYS软件通过算例分析,比较了由转换关系式得到的安全系数与实际计算的结果,讨论了转算结果的误差。算例结果表明：通过转换关系式得到的安全系数与计算得到的结果非常接近,具有比较高的计算精度,同时也证明所提出的方法是可行的。     

A fully coupled element‐free Galerkin model for hydro‐mechanical analysis of advancement of fluid‐driven fractures in porous media

Hydraulic fracturing (HF) of underground formations has widely been used in different fields of engineering. Despite the technological advances in techniques of in situ HF, the industry uses semi‐analytical tools to design HF treatment. This is due to the complex interaction among various mechanisms involved in this process, so that for thorough simulations of HF operations a fully coupled numerical model is required. In this study, using element‐free Galerkin (EFG) mesh‐less method, a new formulation for numerical modeling of hydraulic fracture propagation in porous media is developed. This numerical approach, which is based on the simultaneous solution of equilibrium and continuity equations, considers the hydro‐mechanical coupling between the crack and its surrounding porous medium. Therefore, the developed EFG model is capable of simulating fluid leak‐off and fluid lag phenomena. To create the discrete equation system, the Galerkin technique is applied, and the essential boundary conditions are imposed via penalty method. Then, the resultant constrained integral equations are discretized in space using EFG shape functions. For temporal discretization, a fully implicit scheme is employed. The final set of algebraic equations that forms a non‐linear equation system is solved using the direct iterative procedure. Modeling of cracks is performed on the basis of linear elastic fracture mechanics, and for this purpose, the so‐called diffraction method is employed. For verification of the model, a number of problems are solved. According to the obtained results, the developed EFG computer program can successfully be applied for simulating the complex process of hydraulic fracture propagation in porous media. Copyright © 2016 John Wiley & Sons, Ltd.     

地应力对破裂压力和水力裂缝影响的试验研究

为研究水力压裂时地应力对破裂压力和水力裂缝的影响规律,采用水力压裂试验系统,以相似材料作为研究对象,进行了水力压裂物理模拟试验,并借助煤岩断面三维扫描系统、位移计提取了水力裂缝信息。试验结果表明,随着主应力差的增大,破裂压力逐渐降低,破裂时间也逐渐缩短;随着主应力差的增大,破裂面的表面积逐渐增加,破裂面变得逐渐粗糙;围压相近时水力裂缝易出现转向、分叉,压裂液的动力效应越明显;水力裂缝张开度随着主应力差的增大而逐渐减小。研究结果可为水力压裂试验的进一步研究、裂缝网络系统的建立提供参考。     

Fluid overpressure along an Oligocene out-of-sequence thrust in the Shimanto Belt,SW Japan

Out-of-sequence thrusts (OSTs) exposed in ancient accretionary prisms are considered as fossil analogs of present-day megasplay faults in subduction margins and can provide direct information about the conditions of deformation during thrust activity. In modern as well as in ancient accretionary prisms, first-order megasplay faults or OSTs truncate or merge with faults of lesser importance called second-order OSTs. Structural analysis of the Makinokuchi fault, a branch of an Oligocene to lower Miocene second-order OST in the Tertiary Shimanto Belt of central Kyushu, SW Japan, brings information about the conditions of deformation at the time of thrusting. The studied exposure shows that the fault footwall and, to a much lesser extent, the fault hanging-wall, consist of quartz-cemented syntectonic dilatant hydraulic breccias testifying to pore fluid pressures larger than the least principal stress component. The footwall sandstones are crossed by several centimeters thick quartz veins that merge with the footwall breccias. The continuity between the veins and the breccias suggest that the veins acted as conduits which likely collected fluids from the footwall side sandstones upward and toward the fault. Fluid inclusions indicate that the quartz cementing the breccias and that filling the feeder veins crystallized from similar fluids and under similar pressure and temperature conditions (245–285 °C and 5–8 km depth). These similarities suggest that the fluids responsible for syn-tectonic hydraulic brecciation were collected from the footwall through the conduits. The fluid inclusion trapping temperatures are close to the temperatures expected to be reached along the seismogenic zone. Our analysis shows that fluid overpressures can play a key role in the growth and activity of second-order OSTs in accretionary prisms and suggests that fluids collected along second-order OSTs or splay faults may flow upward along first-order OSTs or megasplay faults.     

湘黔交界中-下寒武统角砾化白云岩的动力学成因

以实际资料为基础,通过对角砾化白云岩基体特征的论述及形成条件的分析,提出了其动力学成因既不是重力滑塌也不是断层破碎而是高孔隙岩层内异常流体高压所导致的水压破裂作用这一新认识;区内中一个寒武统清虚洞组和熬溪组中的角砾化白云岩,岩性独特,层序稳定,以大量方解石贯入胶结为标志,是湘黔边境著名汞矿带的主要含矿层。     

劈裂注浆复合土体平面等效弹性模型理论研究

控制隧道开挖引起的土体沉降变形是劈裂注浆的主要目的之一,隧道在劈裂注浆后复合土体的等效弹性参数取值直接关系着隧道在劈裂注浆后沉降变形的预测精度。首先在对已有劈裂扩散模型研究的基础上,按面积等效原则提出了隧道劈裂注浆后复合土体的二维简化等效单元体模型,并基于均质化理论按变形协调原则推导了二维简化单元体模型的等效弹性参数解析解。然后采用有限元方法分别计算并分析了模型在简化前后的等效弹性参数;同时把二维简化等效单元体模型的有限元计算结果和解析计算结果也做了对比分析。最后基于解析结果分析了土体和浆液结石体各自的弹性参数以及浆液注入率对等效弹性参数的影响。结果表明：（1）按面积等效原则对模型进行简化处理的方法是可行的,可以按照简化模型进行弹性阶段的理论分析;（2）解析结果与有限元结果具有良好的一致性,说明了解析结果的合理性;（3）复合土体的等效弹性模量和等效泊松比主要受浆液注入率和浆液固结体本身模量的影响;浆液固结体的泊松比对等效弹性模量的影响几乎可以忽略。     

汪清油页岩物理力学性质及裂缝起裂压力的研究

本文主要探究了吉林省汪清地区油页岩的物理力学性质及裂缝起裂压力。首先用DNS300型电子万能试验机对油页岩进行了物理力学性质的实验研究,分析并得到了油页岩的各项基本物理力学性质参数的取值范围。其中在垂直层理方向：油页岩的抗压强度范围是21～30 MPa;抗拉强度范围是1.1～3 MPa;弹性模量范围是7.1～9.9 GPa;泊松比范围是0.21～0.27;剪切模量是3.4 GPa。在平行层理方向：抗压强度范围是13～17 MPa;抗拉强度范围是0.3～0.6 MPa;弹性模量范围是1.2～1.6 GPa;泊松比范围是0.31～0.36;剪切模量是0.5 GPa。依据实验结果,选取合适参数,理论计算了汪清油页岩的裂缝起裂压力,并通过真三轴水力压裂装置进行了油页岩水力压裂实验,表明汪清油页岩裂缝起裂压力为5.34 MPa。根据理论计算与实验结果,建议汪清油页岩现场水力压裂实验的起裂压力控制在4～6 MPa。     

基于离散元法的水力压裂数值模拟

不同条件下水压裂隙的发展特性对有效开采页岩气具有重要的指导作用。针对岩体在微观上为颗粒和孔隙的结构系统,提出离散元水力压裂数值模拟方法,离散元能量转化和能量守恒计算方法,建立了相应的三维离散元模型。采用自主研发的三维离散元模拟软件Mat DEM3D,通过控制模型的竖向应变与颗粒直径,来模拟地层中的应力与压裂速率的变化。模拟结果表明:(1)水力压裂产生裂隙的数量和方向受岩石的各向异性,压力状态和变化速率所影响。(2)裂隙在压缩波传播时发展,当水压力高速增加时,诱发的裂隙数量增多,并且有效能量(断裂热)百分比也随之增加,压裂作用也变得更明显。(3)当竖向应变为零时,50%的裂隙呈垂直状态,当竖向应变为-1×10-4时,裂隙趋于沿着最大压力方向发展,竖向裂隙的百分率增大。数值模拟和能量分析为定量地研究岩石水力压裂过程提供了一个新的方法。     

地下采场动态施工损伤演化的数值试验

采用莫尔－库仑准则作为岩体的剪切屈服条件、极限应变作为剪切损伤判据、最大拉应力作为拉伸损伤判据,建立了适合岩体的弹塑性各向异性损伤数值模型。运用该数值模型研究了施工过程不施作支护结构和施作支护结构情况下地下采场覆岩的动态损伤演化过程,并与不施作支护结构时岩体采用弹性模型的数值结果进行了对比。数值试验重现了施工过程地下采场覆岩的动态损伤演化过程,从而为进一步优化支护结构设计、施工以及预防和控制地下采场塌方等工程事故提供了一种新的数值试验手段。