Discrete element model for crack propagation in brittle materials          下载免费PDF全文

Ba Danh Le  Georg Koval  Cyrille Chazallon 《国际地质力学数值与分析法杂志》2016,40(4):583-595

We propose a discrete element model for brittle rupture. The material consists of a bidimensional set of closed‐packed particles in contact. We explore the isotropic elastic behavior of this regular structure to derive a rupture criterion compatible to continuum mechanics. We introduce a classical criterion of mixed mode crack propagation based on the value of the stress intensity factors, obtained by the analysis of two adjacent contacts near a crack tip. Hence, the toughness becomes a direct parameter of the model, without any calibration procedure. We verify the consistency of the formulation as well as its convergence by comparison with theoretical solutions of tensile cracks, a pre‐cracked beam, and an inclined crack under biaxial stress. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Computational model coupling mode II discrete fracture propagation with continuum damage zone evolution          下载免费PDF全文

Wencheng Jin  Hao Xu  Chloé Arson  Seth Busetti 《国际地质力学数值与分析法杂志》2017,41(2):223-250

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.  相似文献   

The effect of expansion ratio on the critical Reynolds number in single fracture flow with sudden expansion          下载免费PDF全文

Jiazhong Qian  Lei Ma  Hongbin Zhan  Qiankun Luo  Xiao Wang  Mu Wang 《水文研究》2016,30(11):1718-1726

Flow in a single fracture (SF) is an important research subject in groundwater hydrology, hydraulic engineering, radioactive nuclear waste repository and geotechnical engineering. An abruptly changing aperture is a unique type of SF. This study discusses the relation between the values of the critical Reynolds number (Re_c) for the onset of symmetry breaking of flow and the expansion ratio (E) of SF, which is defined as the ratio between the outlet (D) and inlet (d) apertures. This study also investigates the effect of inlet aperture d on Re_c for flow in an SF with abruptly changing apertures (SF‐ACA) using the finite volume method. Earlier numerical and experimental results showed that flow is symmetric in respect to the central plane of the SF‐ACA at small Reynolds number (Re) but becomes asymmetric when Re is sufficiently large. Our simulations show that the value of Re_c decreases with the increasing E, and the relationship between the logarithm of Re_c and E can be described accurately using either a quadratic polynomial function or a logarithmic function. However, the relationship of Re_c and d for a given E value is vague, and Re_c becomes even less sensitive to d when E increases. This study also reveals that the hydraulic gradient (J) and flow velocity (v) follow a super‐linear relationship that can be fitted almost perfectly by the Forchheimer equation. The inertial component (J_i) of J increases monotonically with Re, whereas the viscous component (J_v) of J decreases monotonically with Re. The Re value corresponding to equal inertial and viscous components of J (named as the transitional point Re) decreases when E increases, and such a transitional point Re should be closely related to the critical Reynolds number Re_c, although a rigorous theoretical proof is not yet available. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Large‐scale poroelastic fractured reservoirs modeling using the fast multipole displacement discontinuity method          下载免费PDF全文

A. Verde  A. Ghassemi 《国际地质力学数值与分析法杂志》2016,40(6):865-886

An effective approach to modeling the geomechanical behavior of the network and its permeability variation is to use a poroelastic displacement discontinuity method (DDM). However, the approach becomes rather computationally intensive for an extensive system of cracks, particularly when considering coupled diffusion/deformation processes. This is because of additional unknowns and the need for time‐marching schemes for the numerical integration. The Fast Multipole Method (FMM) is a technique that can accelerate the solution of large fracture problems with linear complexity with the number of unknowns both in memory and CPU time. Previous works combining DDM and FMM for large‐scale problems have accounted only for elastic rocks, neglecting the fluid leak‐off from the fractures into the matrix and its influence on pore pressure and stress field. In this work we develop an efficient geomechanical model for large‐scale natural fracture networks in poroelastic reservoirs with fracture flow in response to injection and production operations. Accuracy and computational performance of the proposed method with those of conventional poroelastic DDM are compared through several case studies involving up to several tens of thousands of boundary elements. The results show the effectiveness of the FMM approach to successfully evaluate field‐scale problems for the design of exploitation strategies in unconventional geothermal and petroleum reservoirs. An example considering faults reveals the impact of reservoir compartmentalization because of sealing faults for both geomechanical and flow variables under elastic and poroelastic rocks. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Displacement velocity effects on rock fracture shear strengths     

《Journal of Structural Geology》2016

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.  相似文献   

Universal hydrofracturing algorithm for shear-thinning fluids: Particle velocity based simulation     

《Computers and Geotechnics》2016

A universal particle velocity based algorithm for simulating hydraulic fractures, previously proposed for Newtonian fluids, is extended to the class of shear-thinning fluids. The scheme is not limited to any particular elasticity operator or crack propagation regime. The computations are based on two dependent variables: the crack opening and the reduced particle velocity. The application of the latter facilitates utilization of the local condition of Stefan type (speed equation) to trace the fracture front. The condition is given in a general explicit form which relates the crack propagation speed (and the crack length) to the solution tip asymptotics. The utilization of a modular structure, and the adaptive character of its basic blocks, result in a flexible numerical scheme. The computational accuracy of the proposed algorithm is validated against a number of analytical benchmark solutions.  相似文献   

渝东南下古生界页岩构造裂缝形成及分布控制因素          下载免费PDF全文

岳锋  程礼军  焦伟伟  王飞 《地质科学》2016,(4):1090-1100

天然裂缝对页岩气储集和渗流有重要影响,但目前对页岩构造裂缝研究不够深入。本文通过露头、岩心裂缝观察和分析,探讨了渝东南下古生界页岩构造裂缝类型、形成机理及分布控制因素。基于裂缝特征及力学成因,将剪切裂缝分为高角度剪切裂缝、倾斜滑脱裂缝和水平滑脱裂缝。倾斜滑脱裂缝是在上覆岩层重力和水平构造应力共同作用下沿应力集中的软弱面发生剪切滑动形成,水平滑脱裂缝是在构造挤压应力作用条件下主要沿页理面方向的剪切或层间滑动形成;页理发育程度及岩层曲率是控制水平滑脱裂缝形成的关键因素,岩石矿物组成、构造作用及岩层厚度是控制其它类型构造裂缝形成和分布的主要因素,页岩岩石力学层的划分需要根据裂缝特征、岩性变化特征、岩石力学参数及沉积界面等综合确定,岩层厚度与层间构造裂缝密度呈负相关关系。  相似文献   

SAR图像海岸线检测的区域距离正则化几何主动轮廓模型   总被引：2，自引：0，他引：2  

姜大伟  范剑超  黄凤荣 《测绘学报》2016,45(9):1096-1103

合成孔径雷达(SAR)卫星遥感图像可以极大地提高全国海岸线覆盖频率,然而受到海洋波浪所引起的随机海水表面粗糙度的影响,海岸目标与海水背景边界易混淆不清,因此本文提出了基于区域距离正则化几何主动轮廓模型(RDRGAC),引入距离正则项,解决重复初始化水平集函数为符号距离函数的问题,提高了算法收敛速度。此外,将区域面积项系数与SAR图像等效视数(ENL)建立非线性拟合关系,实现RDRGAC模型根据不同SAR遥感图像的自适应调整,改善海岸线自动提取精度。通过河北省北戴河和大连市金州湾SAR数据海岸线提取对比试验,验证了所提方法的有效性。  相似文献   

声呐合成孔径干涉测量的研究现状与进展     

蒋廷臣  王秀萍  王晓松  焦明连  董春来 《测绘科学》2016,41(11):1-4,30

针对合成孔径声呐是一种新型高分辨水下成像声呐,尤其是其干涉测量是一种高精度测量海底地形的先进技术,具有很大的应用前景,但是目前尚还没有得到广泛应用的问题,详细分析了声呐合成孔径干涉测量近期的研究现状与研究进展情况。该文在声呐合成孔径干涉测量原理的基础上,比较了与合成孔径雷达干涉的不同,重点阐述了相关理论和设备的研究现状和进展,同时分析了现在声呐合成孔径干涉测量所存在的主要问题,最后介绍了目前声呐合成孔径干涉测量的研究展望。  相似文献   

北大巴山闹阳坪锌萤石矿区断裂特征与找矿预测     

张文高  杨兴科  韩鹏飞  樊培贺  何虎军  Moussounda Kounga Claude 《大地构造与成矿学》2016,(2):323-334

闹阳坪锌萤石矿床位于北大巴山逆冲推覆褶皱带平利隆起东侧,为该区首次发现的受构造裂隙控制的气水–热液型锌萤石矿床。本文从矿区断裂入手,通过对该区矿床地质、断裂特征、矿体特征的研究,总结了该区断裂的演化序列及其对成矿的控制机制与规律。认为NW-SE向断裂组是矿区的主要控矿构造,矿区断裂变形发育演化序列为晚印支–早燕山期形成近EW向断裂F1,早燕山期形成NW-SE向断裂F7和NE-SW向断裂F4,随后的右行走滑作用叠加有张扭性应力,形成有利于成矿热液运移充填成矿的张扭性右行走滑断裂。并明确了成矿期应变椭球体,恢复了成矿期主压应力方向为NNW向(340°~350°)。在此基础上,预测平面上闹阳坪矿区F7与F8断裂之间为下一步找矿重点地段,剖面上K3萤石矿体下伏北东方向可能存在隐伏矿体。  相似文献