Three-dimensional discontinuous deformation analysis (3-D DDA) using a new contact resolution algorithm

This paper presents a new contact calculating algorithm for contacts between two polyhedra with planar boundaries in the three-dimensional discontinuous deformation analysis (3-D DDA). In this algorithm, all six type contacts in 3-D (vertex-to-face, vertex-to-edge, vertex-to-vertex, face-to-face, edge-to-edge, and edge-to-face) are simply transformed into the form of point-to-face contacts. The presented algorithm is a simple and efficient method and it can be easily coded into a computer program. In this paper, formulations of normal contact, shear contact and frictional force submatrices based on the new method are derived and the algorithm has been programmed in VC⁺⁺. Examples are provided to demonstrate the new contact rule between two blocks.     

Hydrochemistry, mineralogy and chemical fractionation of mine and processing wastes associated with porphyry copper mines: A case study from the Sarcheshmeh mine, SE Iran

The Sarcheshmeh is one of the largest Oligo-Miocene porphyry Cu deposits in the world. Comparative hydrochemical, mineralogical and chemical fractionation associated with mining efflorescence salts and processing wastes of this mine are discussed. Hydrochemical results showed that rock waste dumps, reject wastes and old impoundments of tailings are the main sources of acid mine drainage waters (AMD) that contain potentially toxic metals such as Cd, Co, Cu, Mn, Ni and Zn as well as Al. Episodic fluxes of highly contaminated acidic waters were produced in a tailings dam over a short period of time. Secondary soluble minerals provide important controls on the quality of AMD produced, especially in old, dry tailings impoundments. Secondary sulfate minerals such as gypsum, magnesiocopiapite, hydronium jarosite, kornelite and coquimbite were found in rock waste drainages and in old weathered reject wastes. Highly soluble secondary minerals such as gypsum, eriochalcite, and bonattite are also observed in an evaporative layer on old tailings impoundments. Chemical fractionation patterns of potentially toxic elements showed that the geochemical behavior of metals is primarily controlled by the mineralogical composition of waste samples. Elements such as Co, Cr, Cu, Mn, Ni and Zn are readily released into the water soluble fraction from efflorescence salts associated with rock waste drainages, as well as from the evaporative layer of old tailings. Potentially toxic elements, such as As, Mo and Pb, are principally adsorbed or co-precipitated with amorphous and crystalline Fe oxides, but they may also be associated with oxidizing, primary sulfides and residual fractions. Following the development of the dammed tailings pond, the secondary minerals were dissolved, producing acidic waters contaminated by Al (154 mg L⁻¹), Cu (150 mg L⁻¹), Cd (0.31 m gL⁻¹), Co (2.13 mg L⁻¹), Mn (73.7 mg L⁻¹), Ni (1.74 mg L⁻¹), Zn (20.3 mg L⁻¹) and Cl (1690 mg L⁻¹). Therefore, the potential use of recycled water from the Sarcheshmenh dammed tailings pond is diminished by the presence of corrosive ions like Cl⁻ in highly acidic fluids that promote corrosion of pipes and pumps in the water recycling system.     

Probabilistic rainfall thresholds for flood forecasting: evaluating different methodologies for modelling rainfall spatial correlation (or dependence)

Rainfall threshold (RT) method is one of the evolving flood forecasting approaches. When the cumulative rainfall depth for a given initial soil moisture condition intersects the threshold rainfall curve, the peak discharge is expected to be equal or greater than the threshold discharge for flooding at the target site. Besides the total rainfall depth, spatial and temporal distribution of rainfall impacts the flood peak discharge and the time to peak. To revisit a previous study conducted by the authors, in which spatially independent rainfall pattern was assumed, the spatial distribution of rainfall was simulated following a Monte Carlo approach. The structure of the spatial dependence among sub‐watersheds' rainfalls was taken into account under three different scenarios, namely independent, bivariate copula (2copula) and multivariate Gaussian copula (MGC). For each set of generated random dimensionless rainfalls, the probabilistic RT curves were derived for dry moisture condition. Results were evaluated with both historical and simulated events. For the simulated events, threshold curves were assessed by means of categorical statistics, such as hit rate, false rate and critical success index (CSI). Results revealed that the best performance based on the CSI criterion corresponded to 50% curve in 2copula and MGC scenarios as well as 90% curve in the independent scenario. The recognition of 50% curve in 2copula and MGC scenarios is in agreement with our expectations that the mean probable curve should have the best performance. Moreover, the proposed inclusion of spatially dependent rainfall scenario improved the performance of RT curves by about 25% in comparison with the presumed spatially uniform rainfall scenario. Copyright © 2011 John Wiley & Sons, Ltd.     

Steady critical discharge rates from vertical and horizontal wells in fresh–saline aquifers with sharp interfaces

In this study, the critical (or maximum) discharge rates before saline water enters a well were determined for vertical and horizontal wells in a freshwater aquifer which is separated from a static saline aquifer by a sharp interface. Flow around the well was solved by integration of a point sink solution along the well axis, and both the critical discharge rate and critical interface rise were determined through a comparison of the heads and vertical gradients at the saline–fresh water interface. The rates were determined for vertical and horizontal wells with various lengths and depths for different aquifer salinities. Results were generalized by drawing dimensionless type curves. The results showed that the dimensionless total critical discharge rates are higher for the longer horizontal wells and longer vertical wells with a certain bottom depth, and they almost linearly decrease with well depth at rates of 0.7–0.9. For the dimensionless well length of 0.2, the dimensionless total discharge rate of a horizontal well is about 0.1 more than that of a vertical well with the same length and well-top depth. Also, the critical discharge rates per unit length of well are inversely proportional to well length and remarkably higher for shallower wells. Additionally, the critical pumping rate is proportional to the salinity difference of the aquifers. These results were confirmed by comparison to existing solutions for vertical wells with dimensionless lengths of 0.2, 0.5 and 0.6, and for critical interface rises in the range of 0.75–0.9.     

Mixed precision parallel preconditioner and linear solver for compositional reservoir simulation

Reservoir simulation role in value creation and strategic management decisions cannot be over emphasized. Simulation of complex challenging reservoirs with millions of grid blocks especially in compositional mode is very time-consuming even with fast modern computers. On the other hand, high price of cluster supercomputers prevents them for being commonly used for fast simulation of such reservoirs. In recent years, the development of many-core processors like cell processors, DSPs, and graphical processing units (GPUs) has provided a very cost-effective hardware platform for fast computational operations. However, programming for such processors is much more difficult than conventional CPUs, and new parallel algorithm design and special parallel implementation methods are needed. Using the computational power of CPUs, GPUs, and/or any other processing unit, Open Computing Language (OpenCL) provides a framework for programming for heterogeneous platforms. In this paper, OpenCL is used to employ the computational power of a GPU to build a preconditioner and solve the linear system arising from compositional formulation of multiphase flow in porous media. The proposed parallel preconditioner is proved to be quite effective, even in heterogeneous porous media. Using data-parallel modules on GPU, the preconditioner/solver runtime reduced at least 1 order of magnitude compared to their serial implementation on CPU.     

The assessment of rock slope instability along the Khosh-Yeylagh Main Road (Iran) using a systems approach

The slope instability is connected to a large diversity of causative and triggering factors, ranging from inherent geological structure to the environmental conditions. Thus, assessment and prediction of slope failure hazard is a difficult and complex multi-parametric problem. In contrast to the analytic approaches, the systems approaches are able to consider infinite number of affecting parameters and assess the interactions of each couple of the parameters in the system. This paper presents a complete application of the rock engineering systems approach in prediction of the instability potential of rock slopes in 15 stations along a 20?km section of the Khosh-Yeylagh Main Road, Iran as the case study of the research. In this research, the main objective has been defining the principal causative and triggering factors responsible for the manifestation of slope instability phenomena, quantify their interactions, obtain their weighted coefficients, and calculate the slope instability index, which refers to the inherent potential instability of each slope of the examined region. The final results have been mapped to highlight the rock slopes susceptible to instability. Finally, as a preliminary validation on the utilization of systems approach in the study region, the stability of investigated rock slopes were analyzed using an empirical method and the results were compared. The comparisons showed a rather good coincidence between the given classes of two methods.     

Optimal water and waste load allocation in reservoir–river systems: a case study

In this paper, a new methodology is developed for optimization of water and waste load allocation in reservoir–river systems considering the existing uncertainties in reservoir inflow, waste loads and water demands. A stochastic dynamic programming (SDP) model is used to optimize reservoir operation considering the inflow uncertainty, and another model called PSO-SA is developed and linked with the SDP model for optimizing water and waste load allocation in downstream river. In the PSO-SA model, a particle swarm optimization technique with a dynamic penalty function for handling the constraints is used to optimize water and waste load allocation policies. Also, a simulated annealing technique is utilized for determining the upper and lower bounds of constraints and objective function considering the existing uncertainties. As the proposed water and waste load allocation model has a considerable run-time, some powerful soft computing techniques, namely, Regression tree Induction (named M5P), fuzzy K-nearest neighbor, Bayesian network, support vector regression and an adaptive neuro-fuzzy inference system, are trained and validated using the results of the proposed methodology to develop real-time water and waste load allocation rules. To examine the efficiency and applicability of the methodology, it is applied to the Dez reservoir–river system in the south-western part of Iran.     

Modeling and assessing the effects of land use changes on runoff generation with the CLUE-s and WetSpa models

Theoretical and Applied Climatology - Land use change is an important determinant of hydrological processes and is known to affect hydrological parameters such as runoff volume, flood frequency,...     

Application of integrated geoelectrical methods in Khenadarreh (Arak,Iran) graphite deposit exploration

This paper is devoted to application of SP, resistivity and mise-a-la-masse method in Khenadarreh (Arak, Iran) graphite deposit exploration. Mise-a-la-masse method is one of the main geophysical exploration methods used to investigate conductive bodies with outcrops on the surface or in exploratory pits. By integrating mise-a-la-masse results with other geoelectrical methods, such as self-potential and resistivity surveys, useful results can be acquired. Powerful survey design, data processing and interpretation methods can make geoelectrical surveying results more reliable. Here, we applied integrated geoelectrical methods to explore the Khenadarreh-bala graphite deposit in the Sarband region of the Markazi province, Iran. We present a new and simple method for processing and interpreting mise-a-la-masse data, which corrects for the effect of distance from the current electrode in the conductive body. After correction, the residual potential values corresponded with the SP results and the apparent resistivity pseudosections. The geophysical results were confirmed by drill holes in the study area.     

Geochemical constraints on the genesis of the volcanic rocks in the southeast of Isfahan area, Iran

Late Miocene–Pliocene to Quaternary calc-alkaline lava flows and domes are exposed in southeast of Isfahan in the Urumieh Dokhtar magmatic belt in the Central Iran structural zone. These volcanic rocks have compositions ranging from basaltic andesites, andesites to dacites. Geochemical studies show these rocks are a medium to high K calc-alkaline suite and meta-aluminous. Major element variations are typical for calc-alkaline rocks. The volcanic rocks have SiO₂ contents ranging between 53.8% and 65.3%. Harker diagrams clearly show that the dacitic rocks did not form from the basaltic andesites by normal differentiation processes. They show large ion lithophile elements- and light rare earth elements (LREE)-enriched normalized multielement patterns and negative Nb, Ti, Ta, and P. Condrite-normalized REE patterns display a steep decrease from LREE to light rare earth elements without any Eu anomaly. These characteristics are consistent with ratios obtained from subduction-related volcanic rocks and in collision setting. The melting of a heterogeneous source is possible mechanism for their magma genesis, which was enriched in incompatible elements situated at the upper continental lithospheric mantle or lower crust. The geochemical characteristics of these volcanic rocks suggested that these volcanic rocks evolved by contamination of a parental magma derived from metasomatized upper lithospheric mantle and crustal melts.