Capturing the complete stress–strain behaviour of jointed rock using a numerical approach

This paper presents the results of a series of numerical experiments using the synthetic rock mass (SRM) approach to quantify the behaviour of jointed rock masses. Field data from a massive sulphide rock mass, at the Brunswick mine, were used to develop a discrete fracture network (DFN). The constructed DFN model was subsequently subjected to random sampling whereby 40 cubic samples, of height to width ratio of two, and of varying widths (0.05 to 10 m) were isolated. The discrete fracture samples were linked to 3D bonded particle models to generate representative SRM models for each sample size. This approach simulated the jointed rock mass as an assembly of fractures embedded into the rock matrix. The SRM samples were submitted to uniaxial loading, and the complete stress–strain behaviour of each specimen was recorded. This approach provided a way to determine the complex constitutive behaviour of large‐scale rock mass samples. This is often difficult or not possible to achieve in the laboratory. The numerical experiments suggested that higher post‐peak modulus values were obtained for smaller samples and lower values for larger sample sizes. Furthermore, the observed deviation of the recorded post‐peak modulus values decreased with sample size. The ratio of residual strength of rock mass samples per uniaxial compressive strength intact increases moderately with sample size. Consequently, for the investigated massive sulphide rock mass, the pre‐peak and post‐peak representative elemental volume size was found to be the same (7 × 7 × 14 m). Copyright © 2015 John Wiley & Sons, Ltd.     

Prediction of rainfall using artificial neural networks for synoptic station of Mashhad: a case study

In this paper, we have utilized ANN (artificial neural network) modeling for the prediction of monthly rainfall in Mashhad synoptic station which is located in Iran. To achieve this black-box model, we have used monthly rainfall data from 1953 to 2003 for this synoptic station. First, the Hurst rescaled range statistical (R/S) analysis is used to evaluate the predictability of the collected data. Then, to extract the rainfall dynamic of this station using ANN modeling, a three-layer feed-forward perceptron network with back propagation algorithm is utilized. Using this ANN structure as a black-box model, we have realized the complex dynamics of rainfall through the past information of the system. The approach employs the gradient decent algorithm to train the network. Trying different parameters, two structures, M₅₃₁ and M₇₄₁, have been selected which give the best estimation performance. The performance statistical analysis of the obtained models shows with the best tuning of the developed monthly prediction model the correlation coefficient (R), root mean square error (RMSE), and mean absolute error (MAE) are 0.93, 0.99, and 6.02 mm, respectively, which confirms the effectiveness of the developed models.     

Landslide susceptibility assessment at the Wuning area,China: a comparison between multi-criteria decision making,bivariate statistical and machine learning methods

Natural Hazards - The aim of this research is to investigate multi-criteria decision making [spatial multi-criteria evaluation (SMCE)], bivariate statistical methods [frequency ratio (FR), index of...     

Pulsatile vortex generators for low-speed maneuvering of small underwater vehicles

Compact zero-mass pulsatile jet actuators are proposed for low-speed maneuvering and station keeping of small underwater vehicles.¹ The flow field of such jets are initially dominated by vortex ring formation. Pinched-off vortices characterize the extremum impulse accumulated by the leading vortex ring in a vortex ring formation process. Relevant parameters in this process are identified in order to design simple and low cost zero-mass pulsatile jet actuators. Thrust optimization of synthetic jets for maximal thrust generation is achieved by enforcing the jet formation number to be around 4. Prototypes of such actuators are built and tested for underwater maneuvering and propulsion. The actuators could be used in two ways: (i) to improve the low-speed maneuvering and station keeping capabilities of traditional propeller driven underwater vehicles, and (ii) as a synthetic jet for flow control and drag reduction at higher cruising speeds. A model for calculating the rotation rate of the underwater vehicle is also proposed and verified.     

3D numerical stability analysis of multi-lateral well junctions

The most important factors in multi-lateral well stability analysis are the magnitude of in situ stresses, the relation between the amount of in situ stresses and orientation of lateral wellbore. In this research, the stability analysis of multi-lateral junction is carried out using FLAC3D numerical code by considering seven varied stress regimes and different lateral wellbore orientations. The Normalized Yielded Zone Area (NYZA, ratio of surrounding yielded cross-sectional area to initial area of well) is determined for different junction mud pressures as well as diverse orientations of lateral wellbore. Then, the junction optimum mud pressure of each lateral wellbore orientation is calculated; hence, the optimum trajectory of lateral wellbore, in which the junction has got the lowest optimum mud pressure, is selected in each stress regime. The stability analysis of multi-lateral wells by means of finite difference method shows that in each stress regime the required mud pressure for the stability of junction is much more than that of the lateral branch and the main wellbore.     

Physical modeling into outflow hydrographs and breach characteristics of homogeneous earthfill dams failure due to overtopping

The main objective of the present study is to introduce new empirical equations for the determination of breach geometrical dimensions and peak outflow discharge(Q_p).Therefore,a historic failure database of 109 embankments was collected and examined.The most important factors that affect the breach evolution,including grading size,hydraulic,and outflow characteristics are also studied.Some of the parameters used for the determination of Q_p and average breach width(Bave)have a significant effect on the erosion process,but they are less reflected in the technical literature.To study the behavior of noncohesive soils during overtopping,15 physical tests were performed at the laboratory,and the effects of interfering parameters were investigated.The experimental output hydrograph was used to simulate the hydrographs resulting from the failure of real dams,and recent artificial intelligence techniques along with linear and nonlinear regression models were employed.The area-time analysis of the laboratory hydrographs shows that the soil particle size and the characteristics of reservoir-basin significantly affect the rate of breach formation and outflow discharge.New relationships are introduced,based on the breach characteristics,by a combination of historical and experimental data,as well as case studies conducted on the hypothetical failure of 10 operational dams.The mathematical model is also used to simulate the process of breach evaluation.Based on statistical indices,comparison of the results,and sensitivity analysis,the developed equations can better express the susceptibility of materials to erosion and their application can minimize downstream vulnerabilities.     

Numerical investigation of the impact of rock mass properties on propagation of ground vibration

Natural Hazards - Ground vibrations induced by the blasting of explosives can cause damage to the slope stability of mines. The important indicators of rock mass such as the geological...     

Studying the stress redistribution around the longwall mining panel using passive seismic velocity tomography and geostatistical estimation

Generally, knowledge of stress redistribution around the longwall panel causes a better understanding of the mechanisms that lead to ground failure, especially to rockbursts. In this paper, passive seismic velocity tomography is used to demonstrate the state of stress around the longwall mining panel. The mining-induced microseismic events were recorded by mounting an array of receivers on the surface, above the active panel. To determine the location of seismic events and execute the process of tomography, double-difference method is employed as a local earthquake tomography. Since passive sources are used, the ray coverage is insufficient to achieve the quality images required. The wave velocity is assumed to be the regionalized variable and it is therefore estimated in a denser network, by using geostatistical estimation method. Subsequently, the three-dimensional images of wave velocity are created and are sliced into the coal seam. These images clearly illustrate the stressed zones that they are appropriately in compliance with the theoretical models. Such compliance is particularly apparent in the front abutment pressure and the side abutment pressure near the tailgate entry. Movements of the stressed zones along the advancing face are also evident. The research conclusion proves that the combined method, based on double-difference tomography and geostatistical estimation, can potentially be used to monitor stress changes around the longwall mining panel continuously. Such observation could lead to substantial improvement in both productivity and safety of mining operations.     

Introducing the recharge–discharge relationship to evaluate the recharge coefficient of karstic aquifers: thematic and indexing approaches

ABSTRACT

Karst aquifers and springs are important with respect to their potential for supplying drinking water in regions suffering from water scarcity in Iran. Accordingly, it is essential to determine the recharge potential of the catchment and the regions with higher obtainability potential. This study provides a road map for the Sheshpeer catchment in southern Iran. A recharge potential (RP) map was produced from which a recharge index (RI) was computed for several selected springs in the catchment. Furthermore, the unit discharge (q) – defined as the average annual discharge for a given catchment area and unit rainfall depth for each spring – was calculated. The plot of q versus RI for the springs showed a linear positive relationship between the two variables (R ² = 0.9). Applying the trend equation of this plot to the whole Sheshpeer karstic catchment reveals that its long-term recharge coefficient is 0.74.     

Comparison of SEBAL,METRIC, and ALARM algorithms for estimating actual evapotranspiration of wheat crop

Theoretical and Applied Climatology - Evapotranspiration is one of the main components of water balance and its accurate estimation is of great importance in planning and optimizing water...