Numerical simulation on mining-induced water inrushes related to geologic structures using a damage-based hydromechanical model

A large number of statistics indicate that water inrush has a direct relationship with geological structures such as fault and karst collapse columns. Understanding the mechanism of water inrushes controlled by geologic structures is of vital importance for adopting effective measures to prevent their occurrence. The work begins with formulization of a damage-based hydromechanical model based on elastic damage theory. Next, the model is numerically implemented with finite element method by employing a finite element package called COMSOL Multiphysics, and is also validated against some existing experimental observations. Finally, the model is used to simulate the mining-induced groundwater inrushes when the effect of faults and karst collapse columns is considered in the numerical simulation, and some suggestive conclusions for preventing water inrushes and optimizing underground mining operations are drawn.     

Landform‐oriented flow‐routing algorithm for the dual‐structure loess terrain based on digital elevation models

The loess landform in the Loess Plateau of China is with typical dual structure, namely, the upper smooth positive terrain and the lower cliffy negative terrain (P–N terrain for short). Obvious differences in their morphological feature, geomorphological mechanism, and hydrological process could be found in the both areas. Based on the differences, a flow‐routing algorithm that separately addresses the dual‐structure terrain would be necessary to encompass this spatial variation in their hydrological behaviour. This paper proposes a mixed flow‐routing algorithm to address aforementioned problems. First, the loess landform surface is divided into P–N terrains based on digital elevation models. Then, specific catchment area is calculated with the new algorithm to simulate the water flows in both positive and negative terrain areas. The mixed algorithm consists of the multiple flow‐routing algorithm (multiple‐flow direction) for positive areas and the D8 algorithm for negative areas, respectively. The approach is validated in two typical geomorphologic areas with low hills and dense gullies in the northern Shaanxi Loess Plateau. Four indices are used to examine the results, which show that the new algorithm is more suitable for loess terrain in simulating the spatial distribution of water accumulation, as well as in modeling the flow characteristics of the true surface by considering the morphological structures of the terrain. Copyright © 2013 John Wiley & Sons, Ltd.     

Application of the numerical manifold method for stress wave propagation across rock masses

In this paper, the numerical manifold method (NMM) is extended to study wave propagation across rock masses. First, improvements to the system equations, contact treatment, and boundary conditions of the NMM are performed, where new system equations are derived based on the Newmark assumption of the space–time relationship, the edge‐to‐edge contact treatment is further developed for the NMM to handle stress wave propagation across discontinuities, and the viscous non‐reflection boundary condition is derived based on the energy minimisation principle. After the modification, numerical comparisons between the original and improved NMM are presented. The results show that the original system equations result in artificial numerical damping, which can be overcome by the Newmark system equations. Meanwhile, the original contact scheme suffers some calculation problems when modelling stress wave propagation across a discontinuity, which can be solved by the proposed edge‐to‐edge contact scheme. Subsequently, the influence of the mesh size and time step on the improved NMM for stress wave propagation is studied. Finally, 2D wave propagation is modelled, and the model's results are in good agreement with the analytical solution. Copyright © 2013 John Wiley & Sons, Ltd.     

Sensitivity of advective travel time of contaminants to correlated formation porosities

An approach is presented to quantify sensitivity of advective contaminant travel time to porosities of hydrogeologic units (HGUs) along the flowpaths when the porosities of different HGUs are correlated. The approach is an extension of the previous sensitivity analysis technique for independent input porosity cases. It is applicable in situations where a calibrated groundwater flow model exists, but a full contaminant transport model is not available. Three sensitivity indices are introduced based on the decomposition of covariance between the advective contaminant travel time and individual input porosities of HGUs. When the input HGU porosities are correlated, the three sensitivity indices quantify the total, intrinsic and correlated contributions from each individual HGU porosity, which should be considered in order to determine the relative importance to the uncertainty in advective contaminant travel time of the input HGU porosities contributing either independently or correlatedly. Two simple one-dimensional flow examples are presented to illustrate the applications of the approach to scenarios when each HGU has distinct porosity and situations of spatially variable porosity field. The approach is applicable to more complex multi-dimensional cases where advective contaminant travel time can be calculated based on simulated flow results from groundwater flow models.     

一种基于重要性评估的河网分级方法

提出了一种改进的河网分级方法,该方法首先对河段的重要性程度进行评估,并利用树结构实现河网的拓扑关系;之后基于河段重要性指数筛选河网主干,并利用主干抽离算法对河网进行逐层遍历;最后基于遍历结果完成河网地层次划分。试验表明,新方法在算法性能和分级结果上均优于传统方法,因此,新方法能够更好地应用于实际的河网分级工作中。     

SAR影像与光学影像的高斯伽玛型边缘强度特征匹配法

提出了一种基于影像边缘强度图描述的SAR影像与光学影像匹配方法。首先对影像进行粗纠正,消除影像之间的尺度和旋转变化;其次,改进相位一致性特征检测方法,提取对影像相干斑噪声稳健的特征点;然后基于高斯伽玛型双边窗口比值算子计算影像边缘强度图,在此基础上构造不变特征描述符;最后联合几何约束条件,实现SAR影像与光学影像匹配。试验结果证明,与现有方法相比,本文方法能够大幅提高SAR影像与光学影像匹配结果中的正确匹配特征数量以及影像配准精度。     

Regional debris flow susceptibility analysis based on principal component analysis and self-organizing map: a case study in Southwest China

A method was developed to analyze the susceptibilities of 541 regional basins affected by debris flows at the Wudongde Dam site in southwest China. Determining susceptibility requires information on source material quantity and occurrence frequency. However, the large number of debris flows can hinder the individual field investigation in a each small basin. Factors that may trigger debris flows can be identified using remotely sensed interpretation information. Susceptibility analysis can then be conducted based on these factors. In this study, SPOT5 satellite imagery, digital elevation models (DEM), a lithology distribution map, and rainfall monitoring data were used to identify 12 debris flow trigger factors: basin relief ratio, slope gradient in the initiation zone, drainage density, downslope curvature of the main channel, vegetation coverage, main channel aspect, topographic wetness index, Melton’s ruggedness number, lithology, annual rainfall, form factor, and cross-slope curvature of the transportation zone. Principal component analysis was used to obtain the eight principal components of these factors that contribute to susceptibility results. Then, a self-organizing map method was adopted to analyze the principal components, which resulted in a debris flow susceptibility classification. Field validation of 26 debris flow basins was used to evaluate the errors of the susceptibility classification, as well as assess the causes of such errors. The study found that principle component analysis and self-organizing map methodologies are good predictors of basin susceptibility to debris flows.