Rock mass excavatability estimation using artificial neural network

One important decision in design of surface mine is the selection of mine equipment and plant. Demand for mechanical excavation is growing in mining industry because of its high productivity and excavation in large scale with lower costs. Several models have been developed over the years to evaluate the ease of excavation and machine performance against rock mass properties. Due to complexity of excavation process and large number of effective parameters, approaches made for this purpose are essentially empirical. There are many uncertainties in results of these models. An attempt is made in this paper to revise the exisiting models. Neural network models for estimation of rock mass excavatability and production rate of VASM-2D excavating machine at Limestone quarry in Retznei, Austria, is presented. Input parameters of this model are Uniaxial compressive strength, tensile strength and discontinuities spacing of rocks. Output is the specific excavation rate per power consumption (bcm/Kwh) as the productivity indicator. Average of deviation between actual data and results estimated by neural network model was only 15% which is in an acceptable range.     

Application of Principal Component Analysis and Spectral Angle Mapper in the Mapping of Hydrothermal Alteration in the Jebal–Barez Area,Southeastern Iran

The southeastern part of the central Iranian Cenozoic magmatic belt contains many areas with copper mineralization. In an analysis of this region, we used the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) and Enhanced Thematic Mapper plus (ETM⁺) images to map the distribution of hydrothermally altered rocks, based on their mineral assemblages. The spectral measurements, based on the field samples and satellite‐image‐derived spectra, show dominantly Al–OH (sericite and clays) and Fe–O absorption features in the visible–near infrared to shortwave infrared regions of the spectrum. Silica–bearing rocks, as well as silicic alteration, show an intense Si–O vibrational feature in the thermal infrared wavelength region. We analyzed ASTER/ETM⁺ images, performed a directed principal component analysis, and used spectral angle mapper to map areas of hydrothermal alteration and iron oxide/hydroxide minerals. The individual principal component images compiled by directed principal component analysis reveal the distribution of individual alteration minerals such as sericite, kaolinite, chlorite, epidote, and quartz. The best results, in terms of mapping the distribution of alteration, were obtained using the spectral angle mapper method. The altered areas were then sampled and the samples subjected to X–ray diffraction analysis, spectral analysis, and thin sections were observed under a microscope. Field observations reveal that more than 98% of the known copper mineralization occurs within the interpreted alteration areas. The present results indicate the great potential of ASTER and ETM⁺ data for mapping the distribution of alteration and exploring for copper mineralization in areas with a similar climate and geological setting to those of the present study.     

Fracture characterization using converted waves

We present a method for inversion of fracture compliance matrix components from wide‐azimuth noisy synthetic PS reflection data and quantitatively show that reflection amplitude variations with offset and azimuth for converted PS‐waves are more informative than P‐waves for fracture characterization. We consider monoclinic symmetry for fractured reservoir (parameters chosen from Woodford Shale), which can be formed by two or more sets of vertical fractures embedded in a vertically transverse isotropic background. Components of effective fracture compliance matrices for a medium with monoclinic symmetry are related to the characteristics of the fractured medium. Monte Carlo simulation results show that inversion of PS reflection data is more robust than that of PP reflection data to uncertainties in our a priori knowledge (vertically transverse isotropic parameters of unfractured rock) than PP reflection data. We also show that, while inversion of PP reflections is sensitive to contrasts in elastic properties of upper and lower media, inversion of PS reflections is robust with respect to such contrasts.     

Internal inlet for wave generation and absorption treatment

A new method of implementing, in two-dimensional (2-D) Navier–Stokes equations, a numerical internal wave generation in the finite volume formulation is developed. To our knowledge, the originality of this model is on the specification of an internal inlet velocity defined as a source line for the generation of linear and non-linear waves. The use of a single cell to represent the source line and its transformation to an internal boundary condition proved to be an interesting alternative to the common procedure of adding a mass source term to the continuity equation within a multi-cell rectangular region. Given the reduction of the source domain to a one-dimensional region, this simple new type of source introduced less perturbation than the 2-D source type. This model was successfully implemented in the PHOENICS code (Parabolic Hyperbolic Or Elliptic Numerical Integration Code Series). In addition, the volume of fluid (VOF) fraction was used to describe the free surface displacements. A friction force term was added to the momentum transport equation in the vertical direction, in order to enhance wave damping, within relatively limited number of cells representing the sponge layers at the open boundaries. For monochromatic wave, propagating on constant water depth, numerical and analytical results showed good agreements for free surface profiles and vertical distribution of velocity components. For solitary wave simulation, the wave shape and velocity were preserved; while, small discrepancy in the tailing edge of the free surface profiles was observed. The suitability of this new numerical wave generation model for a two source lines extension was investigated and proven to be innovative. The comparisons between numerical, analytical and experimental results showed that the height of the merging waves was correctly reproduced and that the reflected waves do not interact with the source lines.     

The impact of salt diapirs on the quality of carbonate karst waters,Bastak, Iran

The Gavbast karstic aquifer located in southern Iran is in direct contact with an exposed salt diapir. To assess the influence of the diapir on the quality of groundwater in the karstic aquifer, electrical conductivity, total dissolved solids, flow rate, temperature and major ion concentrations were measured at 57 sampling sites, including springs, surface waters and wells. A conceptual model of groundwater flow is proposed for the Gavbast karstic aquifer based on the geological setting, water budget, local base of erosion, and hydrochemistry of the sampling sites. The model suggests two subbasins in the Gavbast Anticline draining into two distinct discharging alluvial sections. Unexpectedly, groundwater discharging from the carbonate Gavbast aquifer is saline or brackish and water is of chloride type. The study indicates that the source of salinity of the Gavbast aquifers is infiltration of surface diapir-derived brine into the aquifer. The contribution of the diapir brine in the Gavbast karst aquifer is calculated about 4 L/s, using chloride mass balance. Construction of salt basins to evaporate brine discharging from the diapir springs is proposed to reduce the salinity of karst water. A row of strategically placed wells in the Gavbast karst aquifer would potentially exploit large volumes of fresh groundwater before it is contaminated by the salt. Such low-cost remediation should allow the agricultural exploitation of 40 km² of currently barren land.     

Problematic Soft Soil Improvement with Both Polypropylene Fiber and Polyvinyl Acetate Resin

Several methods are used to improve mechanical properties of loose soils including rewetting, soil replacement, compaction control, chemical additives, moisture control, thermal methods, and more recently, discrete fibers. All the methods are applied to soft soil to increase load bearing capacity and to improve other properties such as prevention of erosion and dust generation. In the present study, a new method of soil improvement using both discrete polypropylene (PP) fibers and polyvinyl acetate (PVAc) is introduced. The method is applied to improve load bearing capacity of a problematic sandy soil in both dry and saturated states. Based on the results from CBR tests on various specimens, it has been revealed that the combination of PP fiber and PVAc resin with weight percentages of 0.1 and 0.6 %, respectively, had the optimum effect in increasing the CBR value in both saturated and dry soil specimens. It should be mentioned that this method has caused a great increase in the CBR value in the saturated soil.     

The maximum expected earthquake magnitudes in different future time intervals of six seismotectonic zones of Iran and its surroundings

One of the crucial components in seismic hazard analysis is the estimation of the maximum earthquake magnitude and associated uncertainty. In the present study, the uncertainty related to the maximum expected magnitude μ is determined in terms of confidence intervals for an imposed level of confidence. Previous work by Salamat et al. (Pure Appl Geophys 174:763-777, 2017) shows the divergence of the confidence interval of the maximum possible magnitude m_max for high levels of confidence in six seismotectonic zones of Iran. In this work, the maximum expected earthquake magnitude μ is calculated in a predefined finite time interval and imposed level of confidence. For this, we use a conceptual model based on a doubly truncated Gutenberg-Richter law for magnitudes with constant b-value and calculate the posterior distribution of μ for the time interval T_f in future. We assume a stationary Poisson process in time and a Gutenberg-Richter relation for magnitudes. The upper bound of the magnitude confidence interval is calculated for different time intervals of 30, 50, and 100 years and imposed levels of confidence α?=?0.5, 0.1, 0.05, and 0.01. The posterior distribution of waiting times T_f to the next earthquake with a given magnitude equal to 6.5, 7.0, and 7.5 are calculated in each zone. In order to find the influence of declustering, we use the original and declustered version of the catalog. The earthquake catalog of the territory of Iran and surroundings are subdivided into six seismotectonic zones Alborz, Azerbaijan, Central Iran, Zagros, Kopet Dagh, and Makran. We assume the maximum possible magnitude m_max?=?8.5 and calculate the upper bound of the confidence interval of μ in each zone. The results indicate that for short time intervals equal to 30 and 50 years and imposed levels of confidence 1???α?=?0.95 and 0.90, the probability distribution of μ is around μ?=?7.16???8.23 in all seismic zones.     

On integral approach to regional gravity field modelling from satellite gradiometric data

Solution of the gradiometric boundary value problems leads to three integral formulas. If we are satisfied with obtaining a smooth solution for the Earth’s gravity field, we can use the formulas in regional gravity field modelling. In such a case, satellite gradiometric data are integrated on a sphere at satellite level and continued downward to the disturbing potential (geoid) at sea level simultaneously. This paper investigates the gravity field modelling from a full tensor of gravity at satellite level. It studies the truncation bias of the integrals as well as the filtering of noise of data. Numerical studies show that by integrating T _zz with 1 mE noise and in a cap size of 7°, the geoid can be recovered with an error of 12 cm after the filtering process. Similarly, the errors of the recovered geoids from T _xz,yz and T _{xx-yy, 2xy} are 13 and 21 cm, respectively.     

The 27 August 2010 Mw 5.7 Kuh-Zar earthquake (Iran): field investigation and strong-motion evidence

In this study, seismological aspects and field observation of the 2010 Kuh-Zar earthquake has been investigated. The Kuh-Zar earthquake, of magnitude 5.7 (Mw), occurred in northeastern Iran on August 27, 2010. The area is surrounded by branches of the active faults which are coated by the quaternary alluvium. During the past several decades, this area has been struck by a number of earthquakes. This earthquake with a moderate magnitude caused a higher rate of damage contrasted with previous earthquakes of the same magnitude range in Iran. Fortunately, the source of the Kuh-Zar earthquake was in a sparsely populated area, and therefore, it caused a few loss of life with the highest observed intensity of shaking VII (modified Mercalli intensity) in the Kuh-Zar village. The shock killed 4 people, injured 40, and destroyed more than 12 villages. According to the field observation, the mechanism of this shock is defined as a left-lateral strike slip. We also checked out the properties of strong ground motions in this earthquake using the records availed by Iranian strong motion network. At KUZ station, about 7 km east of the epicenter, the recorded PGA and PGV in both horizontal and vertical components were remarkably large for an event of this size, and visual inspection of the velocity time history reveals a pulse-like shape. Unfortunately no other recording stations were located close enough to the fault to capture a directivity pulse. Finally, according to the strong-motion properties and observed information, ShakeMaps of the earthquake have been generated by the native intensity observations and the recorded strong motions.