Validating the ability of the discontinuous deformation analysis method to model normal P‐wave propagation across rock fractures          下载免费PDF全文

Wenjun Gong  Yunsheng Wang  Youjun Ning  Yonghong Luo  Yujie Wang 《国际地质力学数值与分析法杂志》2017,41(10):1267-1282

The effects of fractures on wave propagation problems are increasingly abstracting the attention of scholars and engineers in rock engineering field. This study aims to fully validate the ability of discontinuous deformation analysis (DDA) to model normal P‐wave propagation across rock fractures. The effects of a single fracture and multiple parallel fractures are all tested. The results indicate that DDA can accurately reflect the fracture effects, including the fractures stiffness, the fracture spacing and the fracture number, and the effects of incident wave frequency on one‐dimensional P‐wave propagation problems. Thus, DDA is able to deal well with normal incident P‐wave propagation problems. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

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

Enhancement of backtracking search algorithm for identifying soil parameters     

Yin-Fu Jin  Zhen-Yu Yin 《国际地质力学数值与分析法杂志》2020,44(9):1239-1261

In this paper, an enhanced backtracking search algorithm (so-called MBSA-LS) for parameter identification is proposed with two modifications: (a) modifying the mutation of original backtracking search algorithm (BSA) considering the contribution of current best individual for accelerating convergence speed and (b) novelly incorporating an efficient differential evolution (DE) as local search for improving the quality of population. The proposed MBSA-LS is first validated with better performance than the original BSA and some other typical state-of-the-art optimization algorithms on a benchmark of soil parameter identification in terms of effectiveness, efficiency, and robustness. Then, the efficiency of the MBSA-LS is further illustrated by two representative cases: identifying soil parameters from both laboratory tests and field measurements. All comparisons demonstrate that the proposed MBSA-LS algorithm can give accurate results in a short time. Finally, to conveniently solve the problems of parameter identification, a practical tool ErosOpt for parameter identification is developed by integrating the proposed MBSA-LS and some other efficient algorithms for readers to conduct the parameter identification using optimisation algorithms.  相似文献   

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

X射线衍射技术应用于宝石鉴定-合成及晶体结构研究进展   总被引：2，自引：2，他引：0  

张妮  林春明 《岩矿测试》2016,35(3):217-228

因宝石检测常附带无损要求、宝石评价具科学性和不确定性等特点,X射线衍射(XRD)在宝石学这门矿物学的分支学科中的应用具有特殊性,其应用领域主要集中在三个方面。①宝石鉴定领域:对单晶宝石和多晶宝石进行物相鉴定、宝石产地特征研究以及宝石矿物多型的种类划分等。例如根据XRD图谱,可将品种繁多的微晶石英隐晶质变种玉髓中不同的SiO₂质矿物相精确区分;通过黏土矿物成分及其含量可将鸡血石中"地"的种类进行划分。②宝石矿物的晶体结构研究:对成分复杂的单晶宝石进行晶体化学分析、获取矿物类质同象替代情况、利用结晶度划分宝石品质。例如晶格参数中的c/a比值可将绿柱石的Al八面体类质同象替代和Be四面体类质同象替代进行区分;软玉中通过XRD数据计算出的结晶度与软玉的品质关系密切。③宝石学研究技术的开发:对宝石的优化处理过程进行监测,在宝石合成过程中检验矿物晶形大小、结晶质量及内应力等,以有机宝石的晶体结构为基础的人工养殖技术研究。例如尖晶石在热处理过程中,Mg离子由四面体中的T位迁移至八面体中的M位,导致M-O和T-O键长随温度变化,并反映在尖晶石的色调中;在托帕石合成过程中,XRD数据显示770~800℃时形成含氟托帕石(与天然托帕石结构相同),而在950℃和1000℃出现含氟托帕石离解成刚玉和莫来石。目前,我国在宝石学的XRD应用方面起步较晚,国外在宝石鉴定、优化与合成方面的应用相对成熟。本文认为,该技术与扫描电镜、质子探针等分析技术的联用仍有很大发展空间。  相似文献   

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

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