Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis     

Ramin Pakzad  Shanyong Wang  Scott W. Sloan 《国际地质力学数值与分析法杂志》2020,44(10):1472-1494

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.  相似文献   

Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law          下载免费PDF全文

Ahmad Hajjar  Luc Scholtès  Constantin Oltéan  Michel Antoine Buès 《水文研究》2018,32(16):2510-2525

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.  相似文献   

鄂尔多斯地块周缘垂直形变场负位错反演与应变积累研究     

张希  刘立炜  郝明  白卓立 《大地测量与地球动力学》2019,39(6):551-556

利用鄂尔多斯地块及其周缘1970～2014年的垂直形变速率场资料，借助负位错反演研究该区域长期应变积累。结果表明，地块东北缘山西断陷带中北段年均能量积累增量、剪应力强度都较高，西南缘六盘山断裂与渭河断裂西段次之；山西断陷带中南段至晋陕交界处年均剪应力强度较高且显示一定程度的能量积累；西秦岭构造区尤其西秦岭北缘断裂西段、晋冀蒙交界区也反映一定程度的能量积累特性。  相似文献   

内蒙古鲁新井田水文地质条件研究          下载免费PDF全文

赵志强  安美艳  高奎锋  王兆亮  王福芝 《山东国土资源》2018,34(6)

矿产资源开发导致了地下水失衡,地下水失衡又给矿产开发造成了极大的安全隐患,近年来,矿产资源开发与地下水环境保护之间的矛盾愈发凸显。通过对内蒙古鲁新井田典型的水文地质条件进行分析研究,分析采矿导致地下水失衡机理,深入剖析矿井开采充水条件及矿山开发对地下水环境的影响,合理提出了促进矿产开发与地下水保护相互协调的对策建议,为实现"采矿保水"协调统一提供了基础地质依据。为类似地区矿产开发过程中遵循自然规律,趋利避害,保障生产安全,保护地下水环境安全,实现资源绿色开发有较好的指导借鉴作用。  相似文献   

南海神狐海域沉积物自生黄铁矿和石膏S同位素特征及对甲烷渗漏强度的示踪   总被引：1，自引：1，他引：0  

张美  陆红锋  管红香  刘丽华  邬黛黛  吴能友 《海洋学报(英文版)》2018,37(7):20-27

The northern slope of the South China Sea is a gas-hydrate-bearing region related to a high deposition rate of organic-rich sediments co-occurring with intense methanogenesis in subseafloor environments.Anaerobic oxidation of methane(AOM) coupled with bacterial sulfate reduction results in the precipitation of solid phase minerals in seepage sediment,including pyrite and gypsum.Abundant aggregates of pyrites and gypsums are observed between the depth of 667 and 850 cm below the seafloor(cmbsf) in the entire core sediment of HS328 from the northern South China Sea.Most pyrites are tubes consisting of framboidal cores and outer crusts.Gypsum aggregates occur as rosettes and spheroids consisting of plates.Some of them grow over pyrite,indicating that gypsum precipitation postdates pyrite formation.The sulfur isotopic values(δ~(34) S) of pyrite vary greatly(from –46.6‰ to –12.3‰ V-CDT) and increase with depth.Thus,the pyrite in the shallow sediments resulted from organoclastic sulfate reduction(OSR) and is influenced by AOM with depth.The relative high abundance and δ~(34) S values of pyrite in sediments at depths from 580 to 810 cmbsf indicate that this interval is the location of a paleo-sulfate methane transition zone(SMTZ).The sulfur isotopic composition of gypsum(from–25‰ to –20.7‰) is much lower than that of the seawater sulfate,indicating the existence of a 34 S-depletion source of sulfur species that most likely are products of the oxidation of pyrites formed in OSR.Pyrite oxidation is controlled by ambient electron acceptors such as MnO_2,iron(Ⅲ) and oxygen driven by the SMTZ location shift to great depths.The δ~(34) S values of gypsum at greater depth are lower than those of the associated pyrite,revealing downward diffusion of 34 S-depleted sulfate from the mixture of oxidation of pyrite derived by OSR and the seawater sulfate.These sulfates also lead to an increase of calcium ions from the dissolution of calcium carbonate mineral,which will be favor to the formation of gypsum.Overall,the mineralogy and sulfur isotopic composition of the pyrite and gypsum suggest variable redox conditions caused by reduced seepage intensities,and the pyrite and gypsum can be a recorder of the intensity evolution of methane seepage.  相似文献   

鄂尔多斯盆地宜川地区二叠纪孢粉化石组合特征          下载免费PDF全文

高存山 《山东国土资源》2018,34(11)

通过对延长油矿延223井晚古生代孢粉组合研究以及古植被的恢复,探讨了该区晚古生代孢粉植物化石组合的地层学和古气候意义。根据孢粉化石主要属种的时代分布及其在各样品中百分含量的规律性变化,自下而上建立了3个孢粉组合:Gulisporites cochlearius-Laevigatosporites minimus组合、Sinulatisporites-Florinites组合和Lueckeisporites permianus-Platysaccus radialis组合。根据对划分孢粉组合的特征分析,通过与邻区及华北地区晚古生代孢粉组合的对比,探讨了孢粉组合的时代意义。通过孢粉组合,推测了相应时期的古植被类型,并与前人在相同地质时期所建立的大植物化石带进行了对比,分析了各地质时期植被的兴起、发展、繁盛以及衰退绝灭的历史过程,研究区二叠纪植物群属于典型的华夏植物群,反映了一种温暖潮湿的热带雨林气候,二叠纪晚期随着松柏类植物所占比例的增加,气候有变干旱的趋势。  相似文献   

利用短期静态PPP结果反演海潮负荷位移     

魏国光  王琪洁  彭葳  郭良琦 《大地测量与地球动力学》2019,39(6):607-612

采用香港11个GPS测站的观测资料进行1 h、2 h、3 h和4h静态PPP解算，获得4组PPP坐标序列，利用调和分析求取11个测站处8个主要分潮的负荷位移参数（振幅和相位），将其与海潮模型计算的负荷位移参数进行对比，并比较分析PPP反演值与海潮模型值改正海潮负荷信号的效果。结果表明，垂直和水平方向上，不同PPP结果反演8个分潮的负荷位移分别具有约5 mm和7 mm的差异；PPP反演8个分潮垂向负荷位移优于全球海潮模型，但水平方向上的反演效果稍弱。  相似文献   

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock          下载免费PDF全文

N. Zolfaghari  E. Dontsov  A. P. Bunger 《国际地质力学数值与分析法杂志》2018,42(4):587-617

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.  相似文献   

Micromechanical approach to effective viscoelastic properties of micro‐fractured geomaterials          下载免费PDF全文

Cássio B. Aguiar  Samir Maghous 《国际地质力学数值与分析法杂志》2018,42(16):2018-2046

The aim of this paper is to formulate a micromechanics‐based approach to non‐aging viscoelastic behavior of materials with randomly distributed micro‐fractures. Unlike cracks, fractures are discontinuities that are able to transfer stresses and can therefore be regarded from a mechanical viewpoint as interfaces endowed with a specific behavior under normal and shear loading. Making use of the elastic‐viscoelastic correspondence principle together with a Mori‐Tanka homogenization scheme, the effective viscoelastic behavior is assessed from properties of the material constituents and damage parameters related to density and size of fractures. It is notably shown that the homogenized behavior thus formulated can be described in most cases by means of a generalized Maxwell rheological model. For practical implementation in structural analyses, an approximate model for the isotropic homogenized fractured medium is formulated within the class of Burger models. Although the approximation is basically developed for short‐term and long‐term behaviors, numerical applications indicate that the approximate Burger model accurately reproduce the homogenized viscoelastic behavior also in the transient conditions.  相似文献   

Muon density measurements with the KASCADE central detector     

《Astroparticle Physics》2002,16(4):183-386

Frequency distributions of local muon densities in high-energy extensive air showers (EAS) are presented as signature of the primary cosmic ray energy spectrum in the knee region. Together with the gross shower variables like shower core position, angle of incidence, and the shower sizes, the KASCADE experiment is able to measure local muon densities for two different muon energy thresholds. The spectra have been reconstructed for various core distances, as well as for particular subsamples, classified on the basis of the shower size ratio N_μ/N_e. The measured density spectra of the total sample exhibit clear kinks reflecting the knee of the primary energy spectrum. While relatively sharp changes of the slopes are observed in the spectrum of EAS with small values of the shower size ratio, no such feature is detected at EAS of large N_μ/N_e ratio in the energy range of 1–10 PeV. Comparing the spectra for various thresholds and core distances with detailed Monte Carlo simulations the validity of EAS simulations is discussed.  相似文献