Rock Mechanics
E. Khorasani; M. Amini; M. F. Hossaini
Abstract
Bimrocks are complex geomaterials that are defined as mixtures of rocks composed of geotechnically significant blocks within a matrix of finer texture. Bimslopes are made from bimrocks and are usually seen in weathered and shallow environments. Some characteristics of blocks affecting the strength of ...
Read More
Bimrocks are complex geomaterials that are defined as mixtures of rocks composed of geotechnically significant blocks within a matrix of finer texture. Bimslopes are made from bimrocks and are usually seen in weathered and shallow environments. Some characteristics of blocks affecting the strength of bimrocks include VBP (Volumetric Block Proportion), orientation, and arrangement, which have important roles in the stability of bimslopes. Previous studies show that bimrocks usually have a specific block size distribution, and for a bimslope with height of “H”, the size of blocks is changed from 0.05H to 0.75H. In this paper, the influence of large blocks position on bimslope stability was investigated by the physical and numerical models. The blocks that had a dimension larger than 0.5H were considered as “large blocks”. In this work, first, thirty physical models were created and tested using a titling table machine. These models have a specific block size distribution and VBP with ellipsoidal blocks. The main variable of the models is large blocks position, where three categories including lower part of bimslope, upper part of bimslope, and sporadic state are considered. Based on the results of physical trials, thirty numerical models at the laboratory scale were generated using the finite element method. After comparing the physical and numerical models, which showed a good accordance, the numerical models were developed to the natural scale. The theoretical bimslopes investigated in this work showed that the position of large blocks had a significant influence on the stability of bimslopes.
V. Sarfarazi; H. Haeri; M. Fatehi Marji
Abstract
The tensile strength of the anisotropic rock-like material specimens is meastred directly in the laboratory using a new device converting the compressive loading to that of the tensile before the rock breakage. The specially prepared concrete slabs of dimensions 19 cm * 15 cm * 15 cm with a central hole ...
Read More
The tensile strength of the anisotropic rock-like material specimens is meastred directly in the laboratory using a new device converting the compressive loading to that of the tensile before the rock breakage. The specially prepared concrete slabs of dimensions 19 cm * 15 cm * 15 cm with a central hole of 7.5 cm in diameter are tested experimentaly. The specimens are located in the compressive-to-tensile load converting device, and tested under a compressive loading rate of 0.02 MPa/s by the universal testing machine. The cubic slab samples are made in three different configurations to have the directions of 0°, 45°, and -45° with respect to the applied loading direction. In order to compare the direct tensile strength of the concrete samples with that of the indirect measuring tests, some Brazilian tests are also carried out on the concrete disc specimens prepared in the laboratory. By comparing the direct and indirect testing results of the concrete tensile strength, it can be concluded that the direct tensile strength values are somewhat lower than those of the indirect ones. The tensile strength values for the three different configurations of the concrete specimens are nearly the same.
M. Rezaie; S. Moazam
Abstract
Inversion of magnetic data is an important step towards interpretation of the practical data. Smooth inversion is a common technique for the inversion of data. Physical bound constraint can improve the solution to the magnetic inverse problem. However, how to introduce the bound constraint into the inversion ...
Read More
Inversion of magnetic data is an important step towards interpretation of the practical data. Smooth inversion is a common technique for the inversion of data. Physical bound constraint can improve the solution to the magnetic inverse problem. However, how to introduce the bound constraint into the inversion procedure is important. Imposing bound constraint makes the magnetic data inversion a non-linear inverse problem. In this work, a new algorithm is developed for the 3D inversion of magnetic data, which uses an efficient penalization function for imposing the bound constraint and Gauss Newton method to achieve the solution. An adaptive regularization method is used in order to choose the regularization parameter in this inversion approach. The inversion results of synthetic data show that the new method can produce models that adequately match the real location and shape of the synthetic bodies. The test carried out on the field data from Mt. Milligan copper-gold porphyry deposit shows that the new inversion approach can produce the magnetic susceptibility models consistent with the true structures.
S. Mirshrkari; V. Shojaei; H. Khoshdast
Abstract
A coal waste sample loaded with Fe3O4 nanoparticles is employed as an efficient adsorbent to remove Cd from synthetic wastewater. The synthesized nanocomposite is characterized using the Fourier transform-infrared (FT-IR), X-ray diffraction (XRD), and transmission electron microscopy (TEM) techniques. ...
Read More
A coal waste sample loaded with Fe3O4 nanoparticles is employed as an efficient adsorbent to remove Cd from synthetic wastewater. The synthesized nanocomposite is characterized using the Fourier transform-infrared (FT-IR), X-ray diffraction (XRD), and transmission electron microscopy (TEM) techniques. The visual analysis of the microscopic image shows that the mean size of the magnetite nanoparticles is about 10 nm. The effects of the operating variables of the initial solution pH (3-11) and nanocomposite to pollutant ratio (7-233) are evaluated using the response surface methodology on cadmium adsorption. The process is also optimized using the quadratic prediction model based on the central composite design. The statistical analysis reveals that both factors play a significant role in Cd adsorption. The maximum Cd removal of 99.24% is obtained under optimal operating conditions at pH 11 and nanocomposite/cadmium ratio of 90 after 2 h of equilibrium contact time. A study of the adsorption kinetics indicates that the maximum removal could be attained in a short time of about 2 min following a first-order model. The isotherm investigations present that the Cd adsorption on the Fe3O4/coal waste nanocomposite has a linearly descending heat mechanism based on the Temkin isotherm model with the minor applicability parameters than the other isotherm models. The overall removal behaviour is attributed to a two-step mechanism including a rapid adsorption of cadmium ion onto the active sites at the surface of nanocomposite followed by a slow cadmium hydroxide precipitation within the pores over the nanocomposite surface.
E. Farrokh
Abstract
In every tunnel boring machine (TBM) tunneling project, there is an initial low production phase so-called the Learning Phase Period (LPP), in which low utilization is experienced and the operational parameters are adjusted to match the working conditions. LPP can be crucial in scheduling and evaluating ...
Read More
In every tunnel boring machine (TBM) tunneling project, there is an initial low production phase so-called the Learning Phase Period (LPP), in which low utilization is experienced and the operational parameters are adjusted to match the working conditions. LPP can be crucial in scheduling and evaluating the final project time and cost, especially for short tunnels for which it may constitute a major percentage of the total project completion time. The contractors are required to have a better understanding of the initial phase of a project to provide better estimates in the bidding documents. While evaluating and shortening of this phase of low production is important for increasing the productivity and daily advance rate of the machine, there has been limited a direct study and assessment of this period. In this work, we discuss the parameters impacting LPP, and introduce a new methodology for its evaluation. In this regard, an algorithm is introduced for estimation of the approximate extent of LPP based on some TBM tunneling case histories. On the basis of many statistical analyses conducted on the actual data and application of two different shapes of linear and polynomial for the description of LPP, a linear function is proposed for estimation of the learning phase parameters. The major parameters of this function are the learning conditions’ rating and the proportion of LPP to tunnel diameter (X1/D). Analysis of the correlation between these two parameters show a very good coefficient of determination (R2 = 92%). This function can be used for the evaluation of TBM advance rates in LPP and for adjusting the TBM utilization factor in the initial stages of a TBM tunneling project. The learning phase can affect the overall utilization rate and completion time of the tunnels, especially when their lengths are around a couple of kilometers. A true understanding of the LPP characteristics can help the contractors to come up with a more accurate bidding time and cost evaluation, and may also benefit the clients to arrange a better schedule for the final project delivery to the public.
Mohammadhossein Dehghani Firoozabadi; Mohammad Fatehi Marji; Abolfazl Abdollahipour; Alireza Yarahamdi Bafghi; Yousef Mirzaeian
Abstract
The presence of pores and cracks in porous and fractured rocks is mostly accompanied by fluid flow. Poroelasticity can be used for the accurate modeling of many rock structures in the petroleum industry. The approach of the stress to the value of the fracture stress and the effect of pore pressure on ...
Read More
The presence of pores and cracks in porous and fractured rocks is mostly accompanied by fluid flow. Poroelasticity can be used for the accurate modeling of many rock structures in the petroleum industry. The approach of the stress to the value of the fracture stress and the effect of pore pressure on the deformation of rock are among the effects of fluid on the mechanical behavior of the medium. Due to the deformation-diffusion property of porous media, governing equations, strain-displacement, and stress-strain relationships can be changed to each other. In this study, constitutive equations and relationships necessary to investigate the behavior and reaction of rock in a porous environment are stated. Independent and time-dependent differential equations for an impulse and point fluid source are used to obtain the fundamental solutions. Influence functions are obtained by using the shape functions in the formulation of the fundamental solutions and integrating them. To check the validity and correctness of provided formulation, several examples are mentioned. In the first two examples, numerical application and analytical solution are used at different times and in undrained and drained conditions. In times 0 (undrained response of medium) and 4500 seconds (drained response of medium), there is good coordination and agreement between the numerical and analytical results. In the third example, using the numerical application, a crack propagation path in the wellbore wall is obtained, which is naturally in the direction of maximum horizontal stress.
Exploration
Khadijeh Validabadi Bozcheloei; Majid Hashemi Tangestani
Abstract
Evaporites are sediments that chemically precipitate due to the evaporation of an aqueous solution. Most evaporite formations, in addition to evaporite minerals, include detrital rocks such as mudstone, marl, and siltstone. Principal Component Analysis (PCA), Directed Principal Component Analysis (DPCA), ...
Read More
Evaporites are sediments that chemically precipitate due to the evaporation of an aqueous solution. Most evaporite formations, in addition to evaporite minerals, include detrital rocks such as mudstone, marl, and siltstone. Principal Component Analysis (PCA), Directed Principal Component Analysis (DPCA), and Band Ratio methods were applied to Advanced Space-borne Thermal Emission and Reflection Radiometer (ASTER) data for mapping the Gachsaran evaporite formation and distinguishing its lithological units in the Masjed Soleiman oil field, located in southwestern Iran. This oil field was the first recognized oil field in the Middle East. Colour composites of PCs 4, 5, and 2, as RGB images, effectively discriminated this formation from other sedimentary formations. The gypsum spectrum, resampled to the 9 band centres of ASTER, exhibited reflectance in bands 4 and 8 and absorption in bands 6 and 9. As a result, these bands were selected for DPCA application. PC4 effectively highlighted gypsum outcrops as bright pixels, while the band ratio 2/1 accentuated ferric iron, appearing as bright pixels, which correlated with the red marls. The results of this study demonstrate that ASTER image processing is a cost- and time-effective method that can be utilized for mapping evaporite formations and distinguishing their lithological units.
Mineral Processing
M. Maleki Moghaddam; A. R. Hasankhoei; E. Arghavani; A. Haji-Zadeh; M. Yahyaei; S. Banisi
Abstract
Liner design is becoming an increasingly more important tool for the AG/SAG mill performance optimization. The Gol-E-Gohar iron ore concentration plant uses three 9 m × 2.05 m autogenous mills (AG) in parallel in a dry operation. Due to large variations in feed characteristics and inadequate ...
Read More
Liner design is becoming an increasingly more important tool for the AG/SAG mill performance optimization. The Gol-E-Gohar iron ore concentration plant uses three 9 m × 2.05 m autogenous mills (AG) in parallel in a dry operation. Due to large variations in feed characteristics and inadequate blending, the performance of AG mills has been lower than the target value. In order to increase the circuit throughput while maintaining the desired product size, based upon physical and numerical simulations, it was proposed to convert the AG mills to SAG mills. Simulation of the charge trajectory indicated that increasing the liner lifter face angle from 7 to 30° could provide an appropriate charge trajectory in the SAG mode. Installation of the new liners and conversion of AG mill No. 2 to SAG mill, by adding 5% (v/v) balls, resulted in an overall increase of 31% in throughput (from 419 to 548 t/h). Measurement of the wear profiles of shell liners indicated that the wear along the liner length was not uniform. In order to arrive at a uniform wear profile, a new liner design was proposed. Installing the second liner design in AG mill No. 1 and converting it to SAG mill increased the mill throughput by 18% (from 413 to 489 t/h), while the liner life showed a 7% increase. Measurement of the wear profiles of the second liner set indicate that the maximum wear occurs in the centre of the mill. A new liner design was then designed by increasing the width of the lifter top from 12.5 to 15 cm and increasing the lifter height from 16 to 26 cm to enhance the liner life.
Exploitation
S. Mohammadi; M. Ataei; R. Khaloo Kakaie; A. Mirzaghorbanali
Abstract
Immediate roof caving in longwall mining is a complex dynamic process, and it is the core of numerous issues and challenges in this method. Hence, a reliable prediction of the strata behavior and its caving potential is imperative in the planning stage of a longwall project. The span of the main caving ...
Read More
Immediate roof caving in longwall mining is a complex dynamic process, and it is the core of numerous issues and challenges in this method. Hence, a reliable prediction of the strata behavior and its caving potential is imperative in the planning stage of a longwall project. The span of the main caving is the quantitative criterion that represents cavability. In this paper, two approaches are proposed in order to predict the span of the main caving in longwall projects. Cavability index (CI) is introduced based on the hybrid multi-criteria decision-making technique, combining the fuzzy analytical network processes (ANP) and the fuzzy decision-making trial and evaluation laboratory (DEAMTEL). Subsequently, the relationship between the new index and the caving span is determined. In addition, statistical relationships are developed, incorporating the multivariate regression method. The real data for nine panels is used to develop the new models. Accordingly, two models based on CI including the Gaussian and cubic models as well as the linear and non-linear regression models are proposed. The performance of the proposed models is evaluated in various actual cases. The results obtained indicate that the CI-Gaussian model possesses a higher performance in the prediction of the main caving span in actual cases when compared to the other models. These results confirm that it is not possible to consider all the effective parameters in an empirical relationship due to a higher error in the prediction.
Mineral Processing
S. Ghasemi; A. Behnamfard; R. Arjmand
Abstract
The Sangan processing plant consists of four consecutive low-intensity magnetic separation steps with the same magnetic field intensity of 1300 Gauss for upgradation of iron ore. Hence, the iron ore minerals with lower magnetic susceptibility or interlocked with gangue minerals have no opportunity for ...
Read More
The Sangan processing plant consists of four consecutive low-intensity magnetic separation steps with the same magnetic field intensity of 1300 Gauss for upgradation of iron ore. Hence, the iron ore minerals with lower magnetic susceptibility or interlocked with gangue minerals have no opportunity for upgradation, and proceed to the tailing dam. Flotation is a powerful technique for upgradation of these materials, and it is the focus of this research work. A sample of 43.09% Fe and 12.1% FeO was taken from the tailings of second step of magnetic separation. The ore minerals of the sample were determined to be magnetite and hematite. A concentrate of 67% Fe and mass recovery of 50% was produced through the Davis tube test. A reverse flotation route was selected for upgradation of the sample. Fatty acid-based anionic collectors with trade names Alke and Dirol were used in the flotation experiments. The design of experiments was done by resolution IV fractional factorial design with nine factors at two levels per factor. A resolution IV design allows discrimination of all main effects and two-factor interactions. A concentrate of 53.92% Fe at a mass recovery of 60% was obtained at optimum flotation conditions of solid content 20%, pH 12, collector concentration of 1 kg/t, starch as depressant at a concentration of 5 kg/t, Alke/Dirol collector mass ratio of 30/70, conditioning time of 10 min., and concentration of Ca2+ as activator 1 kg/t. In this research work, the concept of natural depression of iron minerals in the reverse flotation was introduced and evaluated.
D. Mohammadi; K. Shahriar; D. Parsapour
Abstract
Tunneling in urban areas has always encountered many uncertainties, which if not considered in both analysis and design of the tunnels, will cause unexpected events during tunnel construction. Obstacles are among the most remarkable uncertainties in tunneling that affect the tunnel construction process. ...
Read More
Tunneling in urban areas has always encountered many uncertainties, which if not considered in both analysis and design of the tunnels, will cause unexpected events during tunnel construction. Obstacles are among the most remarkable uncertainties in tunneling that affect the tunnel construction process. The obstacles in urban tunneling include municipal utilities, surface and sub-surface structures, channels, wells, storages, and unknown cavities. Tehran Metro Line 7 in Iran is no exception to the rule, and has been grappling with the obstacles. In this work, we investigate the effect of the existence of wells and unknown cavities in the zone of influence of excavated tunnels by EPBM. The innovation of this research work is in the EPB tunnel design encountering wells and cavities that are as risky as the adjacent underground structure. In this work, we use a numerical simulation of the 3D finite difference method (FDM) so a series of parametric studies based on the numerical model are examined using the well and unknown cavity geometry and their location relative to the tunnel in alluvium. According to the results obtained, a major disturbance occurs in the near field of the well–tunnel, and the interaction problem happens in front of the tunnel face. The numerical outcome indicates that the most critical state of the ground settlement by EPBM happens when the well and unknown cavity are located in the face of the tunnel. It is also proved that the ground behavior is different for each part of EPBM such as ahead of the face, cutter head, shield, and segmental lining parts.
A. R. Ghanizadeh; A. Yarmahmoudi; H. Abbaslou
Abstract
Due to economical and environmental issues, utilization of mineral wastes, e.g. iron ore mine tailing (IOMT), as road materials can be recommended as a sustainable alternative. In the present study, mechanical properties, as well as resistance to freezing and thawing cycles (F-T) of low plasticity clay ...
Read More
Due to economical and environmental issues, utilization of mineral wastes, e.g. iron ore mine tailing (IOMT), as road materials can be recommended as a sustainable alternative. In the present study, mechanical properties, as well as resistance to freezing and thawing cycles (F-T) of low plasticity clay soil stabilized with different percentages of Portland cement (0, 6, 9, 12 and 15%) and different IOMT content (0, 10, 20, 30 and 40%) has been investigated. To this end, unconfined compressive strength (UCS), initial elastic modulus (E0), and indirect tensile strength (ITS) at different curing times of 7, 14, 18, and 56 days for different admixtures was determined to select optimum mix design for stabilization of clayey subgrade soil. This study shows that by increasing the percentage of cement, strength parameters such as UCS, E0, and ITS increases while increasing IOMT does not show a specific trend to increase strength parameters. Evaluation of strength parameters at different curing time showed that in short-term curing times (7 and 14 days), iron ore mine tailing has a positive effect on the strength parameters, while in long-term curing times (28 and 56 days), iron ore mine tailing has a negative effect on the strength parameters. In total, it was found that 12% of the Portland cement and 10 to 40% of the IOMT passes the UCS and F-T criteria for stabilization of low plasticity clay soils, while clay soil (without IOMT) requires at least 15% of Portland cement for stabilization.
Hassan Vafaie; Seyyed Mohammad Seyyed Alizadeh Ganji
Abstract
The present work is aimed to examine the elimination of cyanide ions from the wastewater derived from the Agh-Darreh gold mine using the Caro’s acid method. The response surface modeling is utilized to evaluate and optimize the influential parameters such as the sulfuric acid/hydrogen peroxide ...
Read More
The present work is aimed to examine the elimination of cyanide ions from the wastewater derived from the Agh-Darreh gold mine using the Caro’s acid method. The response surface modeling is utilized to evaluate and optimize the influential parameters such as the sulfuric acid/hydrogen peroxide ratio, pH, Caro’s acid concentration, and contact time on the elimination process. The results obtained indicate that the increase in the Caro’s acid concentration and contact time has a positive impact on the elimination of the free cyanide ions, while the increment in the weight ratio of sulfuric acid/hydrogen peroxide and pH higher than 9.5 demonstrate a negative impact. Also it is found that the quadratic effect of pH has the highest influence on the removal of cyanide ion, and the linear effect of the ratio of sulfuric acid/hydrogen peroxide has the lowest degree of importance. Additionally, the optimization process is carried out, and about 96.4% of the cyanide ions is eliminated from the wastewater under the optimal conditions including 2 g/L Caro’s acid concentration, 9.3 pH, 8 min contact time, and sulfuric acid to hydrogen peroxide (weight) ratio of 2.
Alireza Dolatshahi; Ali Nouri Qarahasanlou
Abstract
Engineers use various methods to evaluate the performance of concrete structures under dynamic loads, including numerical simulations, laboratory experiments, and field tests. By combining the results of these methods, the engineers can develop a comprehensive understanding of the behavior of concrete ...
Read More
Engineers use various methods to evaluate the performance of concrete structures under dynamic loads, including numerical simulations, laboratory experiments, and field tests. By combining the results of these methods, the engineers can develop a comprehensive understanding of the behavior of concrete structures under dynamic loads and use this information to design more resilient structures capable of withstanding these loads. In this work, four models of the concrete lining of the circular tunnel are simulated to investigate the effect of the pre-cracked in the tunnel's concrete lining under an internal explosion loading. A crack in three different locations at angles of 0, 45, and 90 on the horizontal axis of the tunnel is investigated and analyzed. The coupled Eulerian-Lagrangian method and the constitutive behavior, such as concrete damage plasticity for concrete and Drucker-Prager for soil, allows a more accurate simulation of the internal explosion loading scenario. The selection of Trinitrotoluene and the Jones-Wilkins-Lee equation of state for the explosive provides a realistic representation of the behavior of the explosive material. The modeling results show that in an internal explosion, by examining three different locations of a crack in the concrete, the occurrence of a crack in the crown of the tunnel is more critical than two crack locations. Hence, the existence of a crack with a length of 100 cm and a depth of 15 cm in the crown of the tunnel increases the tensile damage zone by 16.59% compared to the case where there is no crack.
Exploitation
H. Rahimi; O. Asghari; F. Hajizadeh; F. Meysami
Abstract
The purpose of this work is to compare the linear and non-linear kriging methods in the mineral resource estimation of the Qolqoleh gold deposit in Saqqez, NW Iran. Considering the fact that the gold distribution is positively skewed and has a significant difference with a normal curve, a geostatistical ...
Read More
The purpose of this work is to compare the linear and non-linear kriging methods in the mineral resource estimation of the Qolqoleh gold deposit in Saqqez, NW Iran. Considering the fact that the gold distribution is positively skewed and has a significant difference with a normal curve, a geostatistical estimation is complicated in these cases. Linear kriging, as a resource estimation method, can be problematic and gives an unrealistic gold grade. In order to check and correct the errors in the linear methods, the non-linear kriging method has been deployed. One of the applicant's non-linear estimation methods is Indicator Kriging (IK). The IK method converts grade values into binary units of 0 and 1 using multiple thresholds that can be selected by the number-size (N-S) fractal model. The N-S model identifies important and critical thresholds based on the grade distribution. In IK, the Multiple Indicator Kriging (Multiple IK) and Median Indicator Kriging (Median IK) methods could be involved due to the number of indicator thresholds. IK is not sensitive to high values. Here, we make a comparison between Median IK and Multiple IK as well as those with ordinary kriging (OK), which is a linear kriging method. Overall, we conclude that all of these methods are suitable for resource estimation among these methods, although the IK method is better for estimation in different categories of gold grades.
Rock Mechanics
Gh. H Ranjbar; K. Shahriar; K. Ahangari
Abstract
According to the wide application of segmental lining in mechanized tunneling, recognizing the behavior of segmental lining joints is important in tunnels designing. In the structural analysis of the tunnel segmental lining, segmental joints can be considered as elastic joints, and their stiffness characteristics ...
Read More
According to the wide application of segmental lining in mechanized tunneling, recognizing the behavior of segmental lining joints is important in tunnels designing. In the structural analysis of the tunnel segmental lining, segmental joints can be considered as elastic joints, and their stiffness characteristics are affected by the rotational, shear, and axial stiffness. The purpose of this work is to investigate the effect of the rotational, shear, and axial stiffness of segmental lining joints on the internal forces (bending moment and axial force)under the static conditions. For this purpose, a 3D numerical analysis was carried out using the ABAQUS software. The results obtained show that by increasing the rotational stiffness of the segmental joint, the bending moment increases, and for lower values of rotational stiffness, the bending moment variations are higher, while the axial force variations are very slight in comparison with the bending moment. By increasing the axial and shear stiffness of the segmental joint, changes of the bending moment and axial force in segmental lining are negligible.
Mineral Processing
N. Khorasanizadeh; M. Karamoozian; H. Nouri-Bidgoli
Abstract
The bubble diameter effect on the bubble rise velocity profile in a flotation column is studied by the two-phase computational fluid dynamics (CFD) method. The simulations are done in the ANSYS® Fluent® software using a two-phase volume of fluid model. The computational domain is a square cross-section ...
Read More
The bubble diameter effect on the bubble rise velocity profile in a flotation column is studied by the two-phase computational fluid dynamics (CFD) method. The simulations are done in the ANSYS® Fluent® software using a two-phase volume of fluid model. The computational domain is a square cross-section column with a 10 cm width and a 100 cm height, in which air is interred as a single bubble from the lower part of the column by an internal sparger. An experimental test is also performed, the hydrodynamics parameters are recorded, and the simulated results are validated using the values obtained for the bubble rise velocity. The simulation results obtained indicate that CFD can predict the bubble rise velocity profile and its value in the flotation column with less than 5% difference in comparison with the experimental results. Then the simulations are repeated for the other initial bubble diameter in the bubbly flow regime in order to study the bubble diameter effect on the rise velocity profile. The results obtained demonstrate that the larger bubbles reach the maximum velocity faster than the small ones, while the value of maximum velocity decreases by an increase in the bubble diameter. These results can be used to improve the flotation efficiency.
B. Jodeiri Shokri; H. Dehghani; R. Shamsi; F. Doulati Ardejani
Abstract
This work presents a quantitative predicting likely acid mine drainage (AMD) generation process throughout tailing particles resulting from the Sarcheshmeh copper mine in the south of Iran. Indeed, four predictive relationships for the remaining pyrite fraction, remaining chalcopyrite fraction, sulfate ...
Read More
This work presents a quantitative predicting likely acid mine drainage (AMD) generation process throughout tailing particles resulting from the Sarcheshmeh copper mine in the south of Iran. Indeed, four predictive relationships for the remaining pyrite fraction, remaining chalcopyrite fraction, sulfate concentration, and pH have been suggested by applying the gene expression programming (GEP) algorithms. For this, after gathering an appropriate database, some of the most significant parameters such as the tailing particle depths, initial remaining pyrite and chalcopyrite fractions, and concentrations of bicarbonate, nitrite, nitrate, and chloride are considered as the input data. Then 30% of the data is chosen as the training data randomly, while the validation data is included in 70% of the dataset. Subsequently, the relationships are proposed using GEP. The high values of correlation coefficients (0.92, 0.91, 0.86, and 0.89) as well as the low values of RMS errors (0.140, 0.014, 150.301, and 0.543) for the remaining pyrite fraction, remaining chalcopyrite fraction, sulfate concentration, and pH prove that these relationships can be successfully validated. The results obtained also reveal that GEP can be applied as a new-fangled method in order to predict the AMD generation process.
Rahim Mortezaie; Seyed Davoud Mohammadi; Vahab Sarfarazi
Abstract
One of the most important tasks in conducting a laboratory research work is how to make the samples. The purpose of this research work is to create heterogeneous rock-like samples containing non-persistent notches. Regarding that, the molds with dimensions of 250 mm x 200 mm x 50 mm are made. A mixture ...
Read More
One of the most important tasks in conducting a laboratory research work is how to make the samples. The purpose of this research work is to create heterogeneous rock-like samples containing non-persistent notches. Regarding that, the molds with dimensions of 250 mm x 200 mm x 50 mm are made. A mixture of plaster and water with different mixing percentages is used to make the heterogeneous samples. Various techniques are also employed to create non-persistent notches on the samples. One of the methods to create a notch is to insert an aluminum blade into the groove of the mold, and finally, remove it after the plaster slurry has hardened. Due to the displacement of the blade and its tilting during slurring, the notches are out of the vertical position. In addition to the mentioned method, other methods such as water jet, cutting by thread, cutting by diamond wire cutting, cutting by rotary saw, and using hand saw are applied. Finally, using a hand saw to create a notch on the samples is chosen as the best method.
Alireza Dolatshahi; Hamed Molladavoodi
Abstract
The structure's response to the region's prevailing loading conditions guides the engineers in estimating the resilience of the structural materials and their reinforcement. One of the main concerns in designing rock structures is paying attention to the size effect phenomenon. The size effect influences ...
Read More
The structure's response to the region's prevailing loading conditions guides the engineers in estimating the resilience of the structural materials and their reinforcement. One of the main concerns in designing rock structures is paying attention to the size effect phenomenon. The size effect influences the nominal strength, brittleness, load capacity, stress intensity factor, the characteristics of the fracture process zone at the crack tip, and the way and path of crack propagation. Therefore, studying the size effect law will make a guideline for correct decision-making, design, and implementation of efficient support systems. As a comprehensive review, this work investigates specimen size effect on the rock's mechanical and fracture properties. With a comprehensive look at this issue, it explains the essential points that help the engineers design rock structures. During the investigations carried out in this work, it is shown that the specimen size affects the fracture and mechanical properties of the rock. The severity of this phenomenon depends on various factors such as the brittleness index, the shape of the notch or crack length, and the size of the particles that create the rock. In concrete, it depends on the additive boosting materials in the concrete.
Rock Mechanics
vahab sarfarazi; Hadi Haeri; Mohammad Fatehi Marji; Gholamreza Saeedi
Abstract
The mechanical behaviour of transversely isotropic elastic rocks can be numerically simulated by the discrete element method. The successive bedding layers in these rocks may have different mechanical properties. The aim of this research work is to investigate numerically the effect of anisotropy on ...
Read More
The mechanical behaviour of transversely isotropic elastic rocks can be numerically simulated by the discrete element method. The successive bedding layers in these rocks may have different mechanical properties. The aim of this research work is to investigate numerically the effect of anisotropy on the tensile behaviour of transversely isotropic rocks. Therefore, the numerical simulation procedure should be well-calibrated by using the conventional laboratory tests, i.e. tensile (Brazilian), uniaxial, and triaxial compression tests. In this study, two transversely isotropic layers were considered in 72 circular models. These models were prepared with the diameter of 54 mm to investigate the anisotropic effects of the bedding layers on the mechanical behaviour of brittle geo-materials. All these layers were mutually perpendicular in the simulated models, which contained three pairs of thicknesses 5 mm/10 mm, 10 mm/10 mm, and 20 mm/10 mm. Three different diameters for models were chosen, i.e. 5 cm, 10 cm, and 15 cm. These samples were subjected under two different loading rates, i.e. 0.01 mm/min and 10 mm/min. The results gained from these numerically simulated models showed that in the weak layers, the shear cracks with the inclination angles 0° to 90° were developed (considering 15° increment). Also there was no change in the number of shear cracks as the layer thickness was increased. Some tensile cracks were also induced in the intact material of the models. There was no failure in the interface plane toward the layer of higher strength in this research work. The branching was increased by increasing the loading rate. Also the model strength was decreased by increasing the model scale.
Ali Nejati Kalateh; Amin Roshandel kahoo
Abstract
We inverse the surface gravity data to recover subsurface 3D density distribution with two strategy. In the first strategy, we assumed wide density model bound for inverting gravity data and In the second strategy, the inversion procedure have been carried out by limited bound density. Wediscretize the ...
Read More
We inverse the surface gravity data to recover subsurface 3D density distribution with two strategy. In the first strategy, we assumed wide density model bound for inverting gravity data and In the second strategy, the inversion procedure have been carried out by limited bound density. Wediscretize the earth model into rectangular cells of constant andunidentified density. The number of cells is often greater than the number of observation points thus we have an underdetermined inverse problem. The densities are estimated by minimizing a cost function subject to fitting the observed data. The synthetic results show that the recovered model from the first strategy is characterized by broad density distribution around the true model, butthat of the second strategy is closer to true models.We carry out inversion of gravity data taken over chromite deposit located at Hormozgan providence of Iran for estimating of subsurface density distribution. The recovered model obtained from second strategy has appropriate agreement with previous study.
M. J. Babaei; M. A. Molaei; A. Dehghani
Abstract
This study aims to estimate the function of copper consumption using the Johansen approach in time series data, between 1991-2011 in Iran. The literature review of specialized consumption and demand functions shows factors influencing the consumption of copper including copper price variables, aluminum ...
Read More
This study aims to estimate the function of copper consumption using the Johansen approach in time series data, between 1991-2011 in Iran. The literature review of specialized consumption and demand functions shows factors influencing the consumption of copper including copper price variables, aluminum price as a substitute commodity, oil price as a complementary commodity, and industrialization intensity. For this purpose, raw data from the World Bank and International Copper Study Group, were used to extract the initial data needed for the current study, and then research variables were calculated and applied in a seasonal manner. In the next stage, using vector autoregressive, Johansen cointegration test, and vector error correction model, the existence of long-run cointegrated equilibrium relationship was surveyed by vector error correction model. Model estimation results show that there is a positive correlation between between industrialization intensity variable and the price of substitute product (Aluminum) and copper consumption ratio Moreover, there was a significant negative correlation between copper prices and complementary commodity price (Oil), during the review period in Iran. It should be noted that aluminum and copper prices can impact the copper consumption, and an awareness of this can influence making copper sales contracts domestically and abroad.
Exploitation
R. Ghasemi; B. Tokhmechi; G. Borg
Abstract
The known ore deposits and mineralization trends are important key exploration criteria in mineral exploration within a specific region. Fry analysis has conventionally been considered as a suitable method to determine the mineralization trends related to linear structures. Based upon literature sources, ...
Read More
The known ore deposits and mineralization trends are important key exploration criteria in mineral exploration within a specific region. Fry analysis has conventionally been considered as a suitable method to determine the mineralization trends related to linear structures. Based upon literature sources, to date, no investigation has been carried out that includes the Sensitivity Analysis of Feature's Number (SAFN), Sensitivity Analysis of Window Size (SAWS), and Sensitivity Analysis of Spatial Distribution (SASD) of Fry analysis related to mineral locations. In this work, SAFN, SAWS, and SASD are performed by moving several different sub-windows among the main window in order to identify the main trends of mineralization by Fry analysis in the Bavanat region of Iran, which is qualified by its regional and local faults pattern. Based upon our investigation, the effectiveness of the window size and the number of features on Fry analysis are 15-30%. The determined main trends of sub-windows increase, whereas its distribution function of Fry outputs is more similar to the distribution function of Fry outputs of the main window. Moreover, the directions of rose diagrams could be changed due to the edge effects of marginal features around the selected window. However, by selecting an appropriate window, this problem can be solved. Additionally, by an appropriate window selection, the most suitable regional situation is an area that contains the largest number of deposits with a similar metallogenetic origin. Based upon our investigation, the distribution function of the Fry outputs is the main factor that directly controls the identified mineralization pattern of the selected windows.
H. Mahdiyanfar
Abstract
Over the past two decades, the frequency domain (FD) of the geochemical data has been studied by some researchers. Metal zoning is one of the challenging subjects in the mining exploration, where a new scenario has been proposed for solving this problem in FD. Three mineralization areas including the ...
Read More
Over the past two decades, the frequency domain (FD) of the geochemical data has been studied by some researchers. Metal zoning is one of the challenging subjects in the mining exploration, where a new scenario has been proposed for solving this problem in FD. Three mineralization areas including the Dalli (Cu-Au), Zafarghand (Cu-Mo), and Tanurcheh (Au-Cu) mineralization areas are selected for this investigation. After transferring the surface geochemical data to FD, the geochemical signals obtained are filtered using the wavenumber-based filters. The high and moderate frequency signals are removed, and the residual signals are interpreted by the statistical method of principal component analysis (PCA). In order to discriminate the deep metal ore deposits, the principal factors of elemental power spectrum extracted by PCA are depicted in a novel diagram (PC1 vs. PC2). This approach indicates that the geochemical data in the Dalli and Zafarghand deep ore deposits have similar frequency behaviors. The Au, Mo, and Cu elements in these two areas are discriminated from the Au, Mo, and Cu mineralization elements of the Tanurcheh area as a deep non-mineralization zone in this diagram. This new criterion used for distinguishing the buried ore deposits and deep non-mineralization zones is properly confirmed by the exploratory deep drilled boreholes. The geochemical anomaly filtering demonstrates that the strong signatures of deep mineralization are associated with the low frequency geochemical signals at the surface, and the buried mineralization areas with weak surface anomaly can be identified using the geochemical FD data.