Exploration
Parnian Javadi Sharif; Alireza Arab Amiri; Behzad Tokhmechi; Fereydoun Sharifi
Abstract
The technique referred to as Complex Resistivity (CR) or Spectral Induced Polarization (SIP) possesses the capability to distinguish between various kinds of minerals or the sources of induced polarization by utilizing the physical characteristics of minerals or polarizable inclusions. The Generalized ...
Read More
The technique referred to as Complex Resistivity (CR) or Spectral Induced Polarization (SIP) possesses the capability to distinguish between various kinds of minerals or the sources of induced polarization by utilizing the physical characteristics of minerals or polarizable inclusions. The Generalized Effective Medium Theory of Induced Polarization (GEMTip) model is utilized to derive physical characteristics from SIP data. Different inversion methods are applied for this task, though they encounter difficulties such as computational costs, non-linearity, and the intricacy of the inverse issue. To tackle this, a new inversion approach based on Deep Learning (DL) via Convolutional Neural Network (CNN) is proposed for predicting the parameters of polarizable particles from SIP data. The CNN is trained on 20000 synthetic datasets produced using the GEMTip forward model. While DL networks address non-linearities, specific modifications are applied to synthetic datasets to evaluate the influence of non-linearity and correlation on the results. In the Kervian region southwest of Saqqez city, gold mineralization is linked to quartz and pyrite minerals, with two types of pyrite recognized - coarse-grained barren and fine-grained auriferous. The existence of sulfide mineral pyrite, along with variations in pyrite sizes, presents an attractive target for SIP exploration in the investigated area. The trained network is also validated on Gravian data and effectively retrieves parameters as evidenced by the data. The proposed methodology simplifies the inversion process by estimating parameters in one step, enabling a direct and efficient procedure.
Exploration
Shaghayegh Esmaeilzadeh; Ali Moradzadeh; omid Asghari; Reza Mohebian
Abstract
Seismic inversion is a critical technique for estimating the spatial distribution of petro-elastic properties in the subsurface, based on the seismic reflection data. This work introduces an iterative geostatistical seismic inversion method, designed to address challenges in complex geological settings ...
Read More
Seismic inversion is a critical technique for estimating the spatial distribution of petro-elastic properties in the subsurface, based on the seismic reflection data. This work introduces an iterative geostatistical seismic inversion method, designed to address challenges in complex geological settings by incorporating self-updating local variogram models. Unlike the conventional approaches that rely on a single global variogram or fixed local variograms, the proposed method dynamically updates the spatial continuity models at each iteration using automatic variogram modeling and clustering of variogram parameters. The optimal number of clusters is determined using three cluster validity indices: Silhouette Index (SI), Davies-Bouldin Index (DB), and Calinski-Harabasz Index (CH). The method’s effectiveness was evaluated using a three-dimensional non-stationary synthetic dataset, demonstrating robust convergence when employing the SI and CH indices, with both achieving a high global correlation coefficient of 0.9 between the predicted and true seismic data. Among these, the CH index provided the best balance between the computational efficiency and inversion accuracy. The results highlight the method’s ability to effectively capture local spatial variability, while maintaining a reasonable computational cost, making it a promising approach for seismic inversion in complex sub-surface environments.
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.
Rock Mechanics
Amirreza Kavandi; Ramin Doostmohammadi
Abstract
So far, limited research has been conducted on the swelling behavior of Marlstone in the presence of cations. In this study, swelling pressure experiments were performed on rock samples obtained from the Marash Dam, located in northwest Iran. The specimens underwent wetting and drying cycles to achieve ...
Read More
So far, limited research has been conducted on the swelling behavior of Marlstone in the presence of cations. In this study, swelling pressure experiments were performed on rock samples obtained from the Marash Dam, located in northwest Iran. The specimens underwent wetting and drying cycles to achieve an equilibrium condition before cation infiltration. Rock specimens were infiltrated with distilled water and with 1, 2, and 3 mol/L solutions of sodium chloride (NaCl) and calcium chloride (CaCl2). The findings suggest that as the concentration of the solutions rises, the swelling pressure of Marlstone diminishes. Furthermore, at the same concentrations, the swelling pressure of samples soaked in CaCl2 solutions was less than that of those treated with NaCl solutions. Additionally, Marlstone saturated with Ca2+ ions exhibited greater resistance to leaching compared to those saturated with Na+ ions. The findings of this research can be applied to control the swelling pressure of weak rocks in proximity to support systems.
Rock Mechanics
Dariush Kaveh Ahangaran; Kaveh Ahangari; Mosleh Eftekhari
Abstract
Blast damage on the stability of the slopes plays an important role in the profitability and safety of mines. Determination of this damage is also revealed in the widely used Hoek-Brown failure criterion. Of course, this damage is used as a moderating factor in this failure criterion, and its accurate ...
Read More
Blast damage on the stability of the slopes plays an important role in the profitability and safety of mines. Determination of this damage is also revealed in the widely used Hoek-Brown failure criterion. Of course, this damage is used as a moderating factor in this failure criterion, and its accurate determination is considered an important challenge in rock engineering. This study aims to investigate the effect of geological structures in blast damage factor using 3D discrete element modeling of two slopes with different directions of geological discontinuities. The dynamic pressure of the explosion is also simulated in three blastholes. To ensure the modeling results, other dynamic properties of the model have been selected based on the proven studies. An analytical analysis was conducted based on the failure zones (blast damage area), and quantitative and qualitative analyses were performed using the recorded PPV values during the blasting simulation. The results show that the geological discontinuities control, damp, and reduce blast damage. The expansion of blast damage is reduced by 75% along with the increase in rock mass strength, and the blast damage can expand up to 33 meters along with the decrease in strength. By reducing the distance of discontinuities, the role of discontinuities in damping becomes greater than other properties of the rock mass and the discontinuities further away from the blasting hole create more damping. The relation between the distance from the Hole and PPV values shows that for more realistic slope stability analysis results, the values of the damage factor in the Hoek-Brown failure criterion should be applied gradually and decreasingly in layers parallel to the slope surface.
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.
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.
Exploration
Bardiya Sadraeifar; Maysam Abedi; Seyed Hossein Hosseini
Abstract
The Shavaz iron deposit, located in the southwest Yazd province in Central Iranian Block, near The Bafq metallogenic belt, is a significant and economically valuable iron oxide-apatite resource. It features hematite and a minor content of magnetite, detectable through potential field geophysical ...
Read More
The Shavaz iron deposit, located in the southwest Yazd province in Central Iranian Block, near The Bafq metallogenic belt, is a significant and economically valuable iron oxide-apatite resource. It features hematite and a minor content of magnetite, detectable through potential field geophysical surveys. This study aimed to target magnetite mineralization within the deposit using constrained susceptibility inversion. We began by simulating a multi-source synthetic model with three identical cubes at different depths to evaluate the sparse norm inversion approach. The method was then applied to the case study after the essential magnetic data corrections. To refine the interpretation of residual magnetic anomalies and gain insights into their source and depth, the analytic signal and upward continuation methods were employed. Inversion results across different cross-sections revealed two distinct, shallow, lens-shaped magnetite mineralizations with an average vertical extent of 60 meters. Notably, one magnetite body lies approximately 30 meters deeper due to the Dehshir-Baft fault influence. Low normalized mis-fit values confirmed the successful minimization of the objective function during inversion. Additionally, the reconstructed susceptibility models align well with the previous geological studies and borehole data, demonstrating the efficiency of the sparse norm inversion algorithm.
Exploration
shirin Jahanmiri; Ali Aalianvari; Malihehe Abbaszadeh
Abstract
Groundwater inflow is a critical subject within the domains of hydrology, hydraulic engineering, hydrogeology, rock engineering, and related disciplines. Tunnels excavated below the groundwater table, in particular, face the inherent risk of groundwater seepage during both the excavation process and ...
Read More
Groundwater inflow is a critical subject within the domains of hydrology, hydraulic engineering, hydrogeology, rock engineering, and related disciplines. Tunnels excavated below the groundwater table, in particular, face the inherent risk of groundwater seepage during both the excavation process and subsequent operational phases. Groundwater inflows, often perceived as rare geological hazards, can induce instability in the surrounding rock formations, leading to severe consequences such as injuries, fatalities, and substantial financial expenditures. The primary objective of this research is to explore the application of machine learning techniques to identify the most accurate method of forecasting tunnel water seepage. The prediction of water loss into the tunnel during the forecasting phase employed a tree equation based on gene expression programming (GEP). These results were compared with those obtained from a hybrid model comprising particle swarm optimization (PSO) and artificial neural networks (ANN). The Whale Optimization Algorithm (WOA) was selected and developed during the optimization phase. Upon contrasting the aforementioned methods, the Whale Optimization Algorithm demonstrated superior performance, precisely forecasting the volume of water lost into the tunnel with a correlation coefficient of 0.99. This underscores the effectiveness of advanced optimization techniques in enhancing the accuracy of groundwater inflow predictions and mitigating potential risks associated with tunneling activities.
Exploration
Hamid Reza Baghzendani; Hamid Aghajani; Gholam Hossein Karami
Abstract
Karsts are important sources of groundwater, and it is crucial to determine their water volume and quality. The Ravansar Karst spring in the Kermanshah province is a significant water resource with a substantial water volume in the area. The source of this spring is the carbonate rock unit from the Cretaceous ...
Read More
Karsts are important sources of groundwater, and it is crucial to determine their water volume and quality. The Ravansar Karst spring in the Kermanshah province is a significant water resource with a substantial water volume in the area. The source of this spring is the carbonate rock unit from the Cretaceous period and is affected by tectonic changes and faulting caused by movements related to the Zagros folding. In this work, geophysical methods of microgravity, electrical resistivity, and induced polarization have been utilized to identify the extent of karst development in the limestone units. The minimum residual gravity values are associated with karstification. The field dataset comprised two electrical profiles with the dipole- dipole and pole-dipole arrays. The resistivity and gravity data were inverted using a 2D algorithm based on the least square’s technique with a smoothing constraint. According to the processing and 3D modelling of gravity data; not only cavity-shaped voids and spacious cavity chambers were identified but also sub-structures and micro-karstification in carbonate rocks were detected. The most significant finding from the field survey is the detection of low gravimetric values, indicating relatively large holes and chambers that were previously unknown and inaccessible from ground level. These findings are consistent with known collapse and sediment infill features, as seen in surface sinkholes, cavities, and karstification systems. Geophysical surveys and field surveys show that the holes and karsts in the area are related to tectonic phenomena and faulting and are conduits for transporting water to the Ravansar spring.
Environment
Shima Rahimi; Mehdi Irannajad
Abstract
In this study, Red Mud (RM) as a byproduct in alumina production process from bauxite was used as an adsorbent for sulfate contaminant adsorption from acid mine drainage (AMD). AMD discharge led to the acidification of water which has detrimental effects on aquatic life and human health. Analytical methods, ...
Read More
In this study, Red Mud (RM) as a byproduct in alumina production process from bauxite was used as an adsorbent for sulfate contaminant adsorption from acid mine drainage (AMD). AMD discharge led to the acidification of water which has detrimental effects on aquatic life and human health. Analytical methods, laboratory studies and molecular simulations were conducted to investigate sulfate adsorption on RM. Thermodynamic calculations were performed after optimizing of existing metal oxide in RM structure with the Material Studio software using the dmol3 and DFT method. The adsorption energy results by Adsorption locator module determined -819.09, -561.7, -268.8, -105.4, and -314.7 kcal/mol for Fe2O3, Al2O3, CaCO3, TiO2 and SiO2, respectively. The most active compounds in RM structure (iron and aluminum oxides) account for 22.5% and 13.3% in the red mud structure, respectively. In addition, seawater washing was employed as RM modification methods, and it could decrease high rates of pH and improve the sorption capacity of raw RM. The effect of this modification was investigated by simulation of solvent in adsorption environment of sulfate on RM and the dielectric constant selection. For water as the primary solvent with a dielectric constant of 78.54, adsorption energy for RM was calculated to be -35.68 kcal/mol and it was increased to -56.69 kcal/mol for the seawater medium with a dielectric constant of 86. Therefore, RM can be considered as a potential sulfate adsorbent because of cost-effectiveness and alkaline pH that can lead to the neutralization of AMD.
Exploitation
saeideh Qaedrahmat; Javad Gholamnejad; Ali dabagh
Abstract
The scheduling of short-term production in open-pit mining requires determining an optimal extraction sequence for blocks to fulfill multiple goals over a short-term monthly, weekly and daily planning horizon. These goals include meeting required limits on ore grade, production tonnage, waste removal, ...
Read More
The scheduling of short-term production in open-pit mining requires determining an optimal extraction sequence for blocks to fulfill multiple goals over a short-term monthly, weekly and daily planning horizon. These goals include meeting required limits on ore grade, production tonnage, waste removal, and slope constraints. One of the key objectives of Short-Term Production Scheduling (STPS) is to ensure a stable and continuous supply of ore to the processing plant, while minimizing operating costs through measures such as reducing unnecessary equipment movements and variation in feed quality. However, one of the major obstacles to the operational feasibility of STPS is the limited working space available for equipment, as well as the excessive equipment movement between benches within each scheduling period. To tackle these challenges, this paper employs an Integer Goal Programming (IGP) with a new constraint that limits active benches per period, enhancing the practicality of production schedules. Unlike previous GP-based STPS models, it improves operational feasibility by ensuring extraction continuity and minimizing equipment movement. The model was tested on a copper deposit using GAMS software. The results show that by applying this new constraint, the average number of active benches per month was reduced from 14 to 10 )36% reduction) and the number of extraction periods per bench from 6 to 4 (33% reduction) without violating the existing constraints such as ore grade, tonnage, or slope. This approach improves equipment efficiency, reduces fuel consumption, reducing equipment relocation costs, promoting operational continuity of extraction and enhances operational feasibility in real conditions.
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.
