Localization of mining-induced horizontal fractures along rock layer interfaces in overburden: field measurements and prediction

The locations of mining-induced horizontal fractures along rock interfaces in the overburden of Donetsk Coal Basin were identified using an original experimental device. The device traps methane from horizontal fracture zone (100–fold coal seam thickness) over an active longwall mining excavation. Presence or absence of horizontal fractures along rock layer interfaces is correlated with physical characteristics of the overburden, such as thickness, uniaxial compressive strength of overburden rock layers, location of rock layer interfaces and thickness of extracted coal seams. As a result, a combined criterion based on these physical characteristics is proposed to predict the presence of overburden horizontal fracturing in coal mine operations.     

Probabilistic landslide susceptibility mapping in the Lai Chau province of Vietnam: focus on the relationship between tectonic fractures and landslides

This study evaluates the susceptibility of landslides in the Lai Chau province of Vietnam using Geographic Information System (GIS) and remote sensing data to focus on the relationship between tectonic fractures and landslides. Landslide locations were identified from aerial photographs and field surveys. Topographic, geological data and satellite images were collected, processed, and constructed into a spatial database using GIS data and image-processing techniques. A scheme of the tectonic fracturing of crust in the Lai Chau region was established. Lai Chau was identified as a region with many crustal fractures, where the grade of tectonic fracture is closely related to landslide occurrence. The influencing factors of landslide occurrence were: distance from a tectonic fracture, slope, aspect, curvature, soil, and vegetative land cover. Landslide prone areas were analyzed and mapped using the landslide occurrence factors employing the probability–frequency ratio model. The results of the analysis were verified using landslide location data and showed 83.47% prediction accuracy. That emphasized a strong relationship between the susceptibility map and the existing landslide location data. The results of this study can form a basis stable development and land use planning for the region.     

On Tension Failure of 2-D Rock Specimens and Associated Acoustic Emission

Summary To understand the failure mechanism of quasi-brittle materials like rock under tensile stress, observations on the failure process of granite and marble plate specimens under tension are summarized and presented. Micro- and macro-failure properties of rock plates under uniaxial tension have been characterized by using an acoustic emission technique. Acoustic emission signals associated with micro-fractures are captured to locate the sources. An algorithm based on arrival time difference is developed for this purpose. The results reveal clearly the failure processes of rock which include initiation, nucleation and propagation of micro-fractures when the axial stress is close to the peak strength of rock. It is believed that the difference in heterogeneity between granite and marble specimens leads to different fracture shapes and different behaviors of associated acoustic emissions. Numerical simulation of acoustic emissions for two-dimensional tensile test is also carried out. The simulated characteristics are in good agreement with the experimental results.     

Validating the ability of the discontinuous deformation analysis method to model normal P‐wave propagation across rock fractures

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

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

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

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.     

激光测高数据辅助卫星成像几何模型精化处理

资源三号02星搭载了中国首个对地观测试验性激光测高载荷。借鉴目前较成熟的卫星影像区域网平差理论的基础上,结合近年来激光测高数据精度的大幅提升以及资源三号02星激光测高数据的特点,首次提出了激光测高数据辅助卫星立体影像进行成像几何模型精化处理的通用理论。首先,利用传统的区域网平差算法对所处理影像进行高精度连接点匹配处理,并对其进行无约束的自由网平差处理,获得高精度相对精度及不亚于原始成像几何模型的绝对精度;其次,根据激光测高数据3维坐标和精化后参考影像成像几何模型获取激光数据参考影像坐标;而后将参考影像坐标通过几何模型映射获取目标影像上待匹配影像坐标,通过连接点匹配算法,对待匹配目标影像坐标进行精化获取高精度像方同名点;最后,以同名点作为高程控制进行区域网平差计算,对影像成像几何模型进一步处理,获取高精度补偿参数。通过湖北、青海两测区的试验,以激光测高数据辅助卫星影像几何模型精化精度可分别达到1.97 m、3.23 m,结果表明本文提出的方法可有效提高卫星立体数据测图精度。     

Large‐scale poroelastic fractured reservoirs modeling using the fast multipole displacement discontinuity method

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.