Journal of Mining and Environement
Volume:11 Issue: 2, Spring 2020

Oxidation of sulfide-containing ores is the main cause of Acid Mine Drainage (AMD), which is an environmental problem associated with both the abandoned and active mines. Iron-bearing sulfide minerals can be oxidized and form mine waters with high sulfate content, low pH, high electrical conductivity, high redox potential, and high concentrations of iron, aluminum, and other heavy metals. In the process of AMD, precipitation of poorly crystalized oxy-hydroxides of iron with a large active surface can occur. On the surface of iron oxy-hydroxide, the precipitated particulate matter, anions, and cations (metals) could be adsorbed. Mine waters can contain a certain amount of precious metals that can also be adsorbed onto an iron particulate matter surface, which is investigated in this research work. In this work, the samples of iron oxy-hydroxide particulate matter at abandoned gold mine waste in Bakovići (Central Bosnia and Herzegovina) are used. Several parameters including pH, water content, particle size distribution, sulfate content, electrical conductivity, redox potential, and amounts of gold, silver, and iron are measured on the selected mine waste samples. The results obtained indicate that significant amounts of gold (average: 6.8 mg/kg) and silver (average: 7.13 mg/kg) are present in the iron precipitate. Adsorption of precious metals onto the iron oxy-hydroxide surface is strongly pH-dependent. At a very low pH value, desorption of precious metals is favorite. Thus, precious metals are only partially adsorbed onto the iron oxy-hydroxide surface.

Sierra Leone is blessed with abundant natural resources but yet prone to environmental degradation due to the mining operations. Most often, the mining communities are faced with social tensions, as a result of the possible trade-off between the expected employment impact and the cost of mining operations to the environment. Over the past decades, the contribution of the mining sector to the development of the country has been overshadowed by the fact that the mining operations have had adverse negative effects, mainly due to the country's weak environmental policies and the failure of the mine monitoring institutions, to supervise the operations of mining company operations. This article seeks to examine both the environmental and the social implication of mining operations on three mining edge communities in Sierra Leone. This paper also utilizes secondary data from the published articles, government’s reports, workshops and conference proceedings, policy documents of non-governmental organizations, newspapers, and the like to generate this writer's view on the topic under review. The thrust of the review will be on the following: Sierra Rutile Limited, Koidu Holdings Limited, and Shandong Iron Ore Mines. The above mining companies have been carefully selected due to the fact that they are located close to dwelling communities, and have been mining in Sierra Leone over a long period of time. The environmental performance index and the mining impact framework are used to clearly show the impact of mining operations on the environment in Sierra Leone. As a result of mining operations, deforestation is skyrocketing, public discomfort and air pollution has worsened, and social unrest has increased as a result of some unacceptable consequences including pollution of water source without recourse to short-term remedy. The literature reviewed by this writer reveals that the mining activities have two faces in Sierra Leone. One is that it serves as a resource curse. An example to this sad reality is the outbreak of civil war, social unrest among others. On the other hand, the mining sector is one of the principal backbones of the economy. It contributes to the livelihood of the country. This paper introduces three-way approaches of mining sector operation remedies that include but not limited to: 1) sound Environmental Impact Assessment (EIA) adoption before mining operations starts; 2) carrying out Strategic Environmental Assessment (SEA); 3) regular engagement with all stakeholders of mining-affected communities. This article recommends that restoration activities by mining companies go along with extraction and adequate compensation.

This article describes the process of improving risk management practices in a foundry of the ALFET Company (Algeria). The proposed methodology is based on the decision matrix risk assessment technique. This technique allows making a risk assessment for each source of risk (machine, man, environment, and management), which leads to the determination of the overall risk rate during the activity by a new concept. The latter giving a comprehensive vision of occupational health and safety, and compares it with the ALARP principle to determine the acceptability of risk. The main goal of this work is to inculcate a culture on the effects of changing behaviors and attitudes, to disseminate the culture of continuous and sustainable progress within the enterprise, and to ensure that a good atmosphere is maintained in the workplace. It aims to protect and promote the health and safety of workers and the working environment in order to promote a safe and sustainable development company. Our work shows that the working environment is tolerable in terms of health and safety at work. However, to promote a safe and sustainable development in company, an action plan based on the evaluation of the field and feedback through priority actions is recommended for continuous improvement in OSH. Toward the workplace should be continuously monitored to detect risk factors as early as possible before they have negative effects.

Geometallurgy tries to predict the instability the behavior of ores caused by variability in the geological settings, and to optimize the mineral value chain. Understanding the ore variability and subsequently the process response are considered to be the most important functions of an accurate geometallurgical study. In this paper, the geometallurgical index is presented as a new tool to optimize the mining activities. Geometallurgical index is described as any geological feature that makes a footprint on the process performance of ores. In a comprehensive research work at the Sarcheshmeh porphyry copper mine, the geological features that affect the main process responses including the product grade and recovery and plant’s throughput are subjected to investigation. In the current report, the rock hardness variability in terms of semi-autogenous grinding power index (SPI) and its effects on the mill throughput and energy consumption are presented. Ninety samples are collected based on the geological features including lithology, hydrothermal alteration, and geological structures. The samples are mineralogically characterized using XRD, XRF, and electron and optical microscopy. The Starkey laboratory mill, commercialized by Minnovex, is used to perform the SPI comminution test. The SPI results show a wide range of hardness, varying from 12 to 473 minutes. The correlation between the SPI results and the geological features show that lithology is a key geological feature that defines the hardness variability. In addition, the hydrothermal alteration would be an effective parameter in the period that the plant is fed with a single lithology.

The shell liner type, rotation speed, and ball filling percent are the key factors influencing the charge behavior inside the SAG mills, and consequently, their performance. In this work, the milling operation of industrial SAG mills is investigated using the Discrete Element Method (DEM). First, an industrial SAG mill with dimensions of 9.50 m × 4.42 m that has a Smooth-type liner is simulated. Then by changing the liner types, i.e. Wave, Rib, Ship-lap, Lorain, Osborn, and Step liners, six other independent simulations are performed. In order to investigate the impact mechanism and improve the mill performance, two new parameters called ‘head height’ and ‘impact zone length’ are introduced. Then the effects of the mill shell liner type on those parameters at two different mill speeds, i.e. 70% and 80% of its critical speed (CS), are evaluated. Also for validation of the simulation results, a laboratory-scale SAG mill with dimensions of 57.3 cm × 16.0 cm is simulated. The results obtained indicate that the Osborn liner, due to the angularity of its lifters and their proper number and thickness, performs best because it increases both parameters more than the other liners. Thus this liner is recommended as the best and optimal liner in this research work and is suggested for installation inside the industrial SAG mills. Also the Wave liner, due to its specific geometrical shape and its wavy lifters as well as their low number and inadequate thickness, provides the lowest charge ‘head height’. Therefore, it is not recommended to install this liner inside the industrial SAG mills. Meanwhile, comparison of the simulations related to the laboratory-scale SAG mill with the experimental results demonstrates a good agreement that validates the DEM simulations and the software used.

Hyperspectral remote sensing records reflectance or emittance data in a large sum of contiguous and narrow spectral bands, and thus has many information in detecting and mapping the mineral zones. On the other hand, the geological and geophysical data gives us some other fruitful information about the physical characteristics of soil and minerals that have been recorded from the surface. The Sarcheshmeh mining area located in the NW-trending Uromieh-Dokhtar magmatic belt within Central Iran is mainly of porphyry type, and is associated with extensive hydrothermal alterations. Due to the semi-arid type of climate with abundant rock exposure, this area is suitable for application of remote sensing techniques. In this work, we focus on generating the alteration maps around Cu porphyry copper deposits using the spectral angle mapper algorithm on Hyperion data by applying two filters named reduction to pole and analytical signal on a total magnetic intensity map and generating the Kd map from radiometry data. What is clear is the high importance of applying the adequate pre-processing on Hyperion data because of low signal-to-noise ratio. By comparing the known deposits in the region with the results obtained by applying the mentioned methods, it is revealed that not all the higher K radiometric values are entirely associated with the hydrothermal alteration zones, and in contrast, the potassic alteration map extracted from Hyperion imagery successfully corresponds to the alteration zones around the Sarcheshmeh mining area. Finally, the results particularly obtained from processing the Hyperion data are confirmed by indices of Cu porphyry deposits in the region.

Displacements around a tunnel, occurring as a result of excavation, consist of the elastic and plastic parts. In this paper, we discuss the elastic part of displacements as a result of excavation, called net displacement. In general, the previous analytical solutions presented for determining the displacements around a circular tunnel in an elastic medium do not give the net displacements directly. The well-known Kirsch solution is the most widely used method for determining the induced stresses and net displacements around a circular opening in a biaxially-loaded plate of homogeneous, isotropic, continuous, linearly elastic material. However, the complete solution for obtaining the net displacements has not been presented or highlighted in the available literature. Using the linear elasticity, this paper reviews and presents three different analytical methods for determining the net displacements directly as well as induced stresses around a circular tunnel. The three solution methods are the Lame' method, airy stress function method, and complex variable method. The tunnel is assumed to be situated in an elastic, continuum, and isotropic medium in the plane strain condition. The solutions are presented for both the hydrostatic and non-hydrostatic in situ stresses in the 2D biaxial loading condition along with an internal pressure. Loading and unloading in tunneling occurring as a result of excavation and stress differences between the induced and initial ones are considered to evaluate the net displacements directly. Finally, some examples are given to demonstrate the complete solution and show the difference between the net elastic displacements as a result of excavation and total elastic displacements that are not real.

One of the methods used to investigate the damaged zone in rock structure is the acoustic emission method. This method is based on receiving the elastic waves that are produced by deformation and cracking of the rock mass around the underground excavation. In this research, a study is conducted on the rock samples by a numerical method to investigate the damaged zone caused by the excavation of circular space on it. For this purpose, 33 cube samples of three different material types including sandstone, concrete, and cement-plaster mortar are prepared. A circular hole is drilled in the center of each sample. The hole diameter is 20 or 25 mm. The samples are loaded uniaxially or biaxially with different stress rates. It is tried to study the acoustic events occurring in the samples during the test, and their locations are investigated. Then the experiments are evaluated by a numerical method using the FLAC3D software and some developed codes. The relation between the sample damaged zone where the acoustic events have occurred during the loading period and the numerical elements that reach a degree of tensile and shear yield is studied. The results obtained show that the amount of cumulative acoustic parameters in cement-plaster mortar specimens is more than the others. In fact, the finer grains, the more amounts of energy and counts will be produced. Also, the results show that with increase in the lateral pressure and loading rate, the amount of cumulative energy and counts decreases.

Detection of mineralized zones based on ores and gangues is important for mine planning and excavation operation. The major goal of this research work was to determine the zones based on ores and gangues by a combination of fractal and factor analysis in the Chah Gaz iron ore (Central Iran). The Concentration-Volume (C-V) fractal method was carried out for Fe, P and S, which indicated that the main mineralized zones consisted of the Fe, S, and P values ≥ 57%, ≤ 0.4%, and ≤0.3%, respectively. Factor analysis categorized variables in two groups including factor 1 (F1) and factor 2 (F2) for ore and gangue, respectively. The C-V fractal modeling on the derived factors showed four zones for F1 and F2. Based on the correlation among the results of fractal modeling on the elements and factors, the first and second zones of F1 were proper for exploitation. Furthermore, the last and first zones of F1 and F2 could be assumed as the main waste for mining excavation.

In this research work, a 3D numerical modeling technique is proposed based on the 3D particle flow code in order to investigate the failure mechanism of rock foundations. Two series of footings with different geometries and areas are considered in this work. The failure mechanism obtained is similar to that of the Terzaghi’s but there is a negligible difference in between. Lastly, one equation is presented to calculate the bearing capacity based on the results achieved from the numerical model and the Mohr-Coulomb theory. The sensitivity analyses are performed on the friction angle, cohesion, and footing width. The results obtained are compared with the corresponding results given by the equations given by Terzaghi and Meyerhof. This comparison demonstrates a good agreement between them. In the friction angle sensitive analysis, the amounts of the bearing capacity diagram are very close to Meyerhof’s, which overlap with each other.

An accurate modeling of sophisticated geological units has a substantial impact on designing a mine extraction plan. Geostatistical simulation approaches, via defining a variogram model or incorporating a training image (TI), can tackle the construction of various geological units when a sparse pattern of drilling is available. The variogram-based techniques (derived from two-point geostatistics) usually suffer from reproducing complex and non-linear geological units as dyke. However, multipoint geostatistics (MPS) resolves this issue by incorporating a training image from a prior geological information. This work deals with the multi-step Single Normal Equation Simulation (SNESIM) algorithm of dyke structures in the Sungun Porphyry-Cu system, NW Iran. In order to perform a multi-step SNESIM algorithm, the multi-criteria decision-making and MPS approaches are used in a combined form. To this end, two TIs are considered, one for simulating dyke structures in the shallow depth, and two for simulating dyke structures in a deeper depth. In the first step, a TI is produced using geological map, which has been mined out during the previous exploration operations. After producing TI, the 35 realizations are simulated for the shallow depth of deposit in the area under study. To select the best realization (as a TI for the next step) of the simulation results, several statistical criteria are used and the results obtained are compared. To this end, a hybrid multi-criteria decision-making is designed on the basis of a group of statistical criteria. In the next step, the dyke structures in the deeper depth are also simulated by the new TI.

The uniaxial compressive strength of weak rocks (UCSWR) is among the essential parameters involved for the design of underground excavations, surface and underground mines, foundations in/on rock masses, and oil wells as an input factor of some analytical and empirical methods such as RMR and RMI. The direct standard approaches are difficult, expensive, and time-consuming, especially with highly fractured, highly porous, weak, and homogeneous rocks. Numerous endeavors have been made to develop indirect approaches of predicting UCSWR. In this research work, a new intelligence method, namely relevance vector regression (RVR), improved by the cuckoo search (CS) and harmony search (HS) algorithms is introduced to forecast UCSWR. The HS and CS algorithms are combined with RVR to determine the optimal values for the RVR controlling factors. The optimized models (RVR-HS and RVR-CS) are employed to the available data given in the open-source literature. In these models, the bulk density, Brazilian tensile strength test, point load index test, and ultrasonic test are used as the inputs, while UCSWR is the output parameter. The performances of the suggested predictive models are tested according to two performance indices, i.e. mean square error and determination coefficient. The results obtained show that RVR optimized by the HS model can be successfully utilized for estimation of UCSWR with R2 = 0.9903 and MSE = 0.0031203.

Copper smelting and refinery factories are the final stages of a pyrometallurgical processing chain, and they cause many environmental challenges around the world. One of the most common environmental problems of these factories is toxic emissions. These toxic gases have harmful effects on the vegetation, animal species, soils, and water resources around the factories. Phytoremediation can play an important role in the reduction of the adverse effects of environmental pollutions arising from copper smelting and refinery factories. In this paper, we first discuss different types of pollutions caused by copper metallurgical factories, and present the main research approaches and studies conducted on these factories. In the second part, we provide a summary and comparison of different remediation technologies used to reduce the environmental pollutions of these factories. Besides, the advantages and disadvantages of each method is also investigated. In the third part, we review the different aspects of the phytoremediation including the effective mechanisms, different types of plants, application environments, and the effective factors. The next part includes the selection of suitable plants for the phytoremediation process applied for copper metallurgical factories and investigation of the native and cultivated hyperaccumulator plants. In addition, different efficiency indices are introduced for evaluating the phytoremediation efficiency and selecting an appropriate hyperaccumulator plant. At the final stage, some appropriate plant species for various types of phytoremediation are introduced. The effects of different environmental stresses and the possibilities of integrating phytoremediation with other remediation technologies as well as the advantages and disadvantages of phytoremediation are eventually investigated.

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.

Identification of the vertical and horizontal distributions for elemental grades is of an important sign in different mineral exploration stages. The main aim of this work is to determine the vertical distribution directional properties of Cu values in the Milloieh Cu porphyry deposit, Kerman (SE Iran) using the power-law frequency fractal model. This work is carried out based on four mineralized boreholes. The Cu grade vertical distribution in mineralized boreholes indicates a positively skewed distribution in the former and multi-modal distribution in the latter types. The power-law frequency analysis in nature represents that the Cu values in four mineralized boreholes are bifractal. The two sections of these graphs outline a break point of about 0.5-1% for Cu values and fractal dimension range of 1.70-4.97 in the mineralized boreholes.

One of the most crucial factors involved in the optimum design and cost estimation of rock sawing process is the rock abrasivity that could result in a significant cost increase. Various methods including direct and indirect tests have been introduced in order to measure rock abrasivity. The Schimazek’s F-abrasiveness factor ( ) is one of the most common indices to assess rock abrasivity.  is the function of three rock parameters including the Brazilian tensile strength ( ), median grain size ( ), and equivalent quartz content ( ). By considering its formulation, it has been revealed that the coefficient of each parameter is equal, which is not correct because each parameter plays a different role in the rock abrasion process. This work aims to modify the original form of  by introducing three correction factors. To calculate these correction factors, an integrated method based on a combination of the statistical analysis and probabilistic simulation is applied to a dataset of 15 different andesite rocks. Based on the results obtained, the values of -0.36, 0.3, and -0.89 are suggested as the correction factors of ,  and , respectively. The performance of the modified Schimazek’s F-abrasiveness factor ( ) is checked not only by the wear rate of diamond wire but also by the cutting rate of the wire sawing process of Andesite rocks. The results obtained indicate that the wear rate and cutting rate of andesite rocks can be reliably predicted using . However, it should be noted that this work is a preliminary one on the limited rock types and further studies are required by incorporating different rock types.

Waste rock dumping is very important in the production planning of open-pit mines. This subject is more crucial when there is a potential of acid-forming (PAF) by waste rocks. In such a type of mines, to protect the environment, the PAF materials should be encapsulated by non-harmful rocks. Therefore, block sequencing of the mined materials should be in such a way that both the environmental and economic considerations are considered. If non-acid forming (NAF) rocks are not mined in a proper time, then a stockpile is required for the NAF materials, which later on would be re-handled for encapsulation of PAF rocks. In the available models, the focus is on either block sequencing or waste dumping strategy. In this work, an attempt has been made to develop an integrated mathematical model for simultaneous optimization of block sequencing and waste rock dumping. The developed model not only maximizes the net present value (NPV) but also decreases the destructive environmental effects of inappropriate waste dumping. The proposed model, which is solved by a CPLEX engine, is applied to two different iron deposits. Also the performance of the proposed model is cross-checked by applying the available (traditional) models in a two-step manner. According to the results obtained, it can be considered that utilizing the developed model, because of extensive re-handling cost reduction, the NPV improvement is significant, especially when the overall stripping ratio is higher (deposit case A).

Stations are the main components of the subway systems. Despite the progress in the construction and maintenance, stations have always been exposed to the natural and man-made disasters. In such incidents, the station’s evacuation capability has a direct relation with a passenger's life. Various factors affect the station's evacuation capability. Investigation of these factors and evaluation of the station’s evacuation capability have important roles in protecting a passenger's life. For this purpose, the catastrophic events that lead to the evacuation of a station and the factors affecting the evacuation of the station are identified. Due to the difference in the catastrophic event probabilities at each station, the risk of catastrophic events is evaluated. Then the station score is calculated according to the value and weight of the evacuation factors and the wighted influence of the catastrophic events. Accordingly, the proposed model is implemented in the Tehran subway. Based on the results obtained, uncrowded stations, even though served by a small number of passengers, may also have a low evacuation capacity and lead to casualties in an emergency situation. This is due to the lack of emergency management and safety facilities. Also by assessing the risk of catastrophic events at stations and equipping stations on its basis, the degree of safety and evacuation capability can be improved more effectively. The sensitivity analysis of the evacuation factors show that the most effective way to increase the station’s evacuation capability is to improve its status in management factors. Using the proposed model to evaluate the station's evacuation capability is an appropriate method for identifying the stations that have a poor evacuation capability.

The computer-based 3D modeling of ore bodies is one of the most important steps in the resource estimation, grade determination, and production scheduling of open-pit mines. In the modeling phase, the volume of the orebody model is required to be filled by the blocks and sub-blocks. The determination of Block Size (BS) is important due to the dependence of the geostatistical issues and calculations related to mining capabilities on it. There are some factors effective in the determination of an optimal BS including the metal content, estimation error, recovery percentage, mining ability, safety, and dilution. In this work, an optimal BS is determined using a two-stage approach. In the proposed approach, the Fuzzy Delphi Analytic Hierarchy Process (FDAHP) and Fuzzy Multi-Objective Optimization by Ratio Analysis (FMOORA) methods are used. In the first phase, the weight of each criterion is calculated based on the opinions of the experts using the FDAHP method. In the second phase, the FMOORA method is applied in order to determine a suitable BS for the design and operation of mining considering the extracted weights in the previous phase. The block model of the Sungun copper mine is studied as a case study to evaluate the capability of the proposed approach. The results of implementation of this approach are desirable because of converting the opinions of the experts to fuzzy values, weighing the experts according to the experience and technical knowledge, weighting the criteria by FDAHP, and choosing the optimal option with FMOORA. Furthermore, the 12.5×12.5×12.5 m3 block (A5) is chosen as an appropriate BS, which is compatible with the real conditions of the studied mine.

The uniaxial compressive strength (UCS) of intact rocks is one of the key parameters in the course of site characterizations. The isotropy/anisotropy condition of the UCS of intact rocks is dependent on the internal structure of the rocks. The rocks with a random grain structure exhibit an isotropic behavior. However, the rocks with a linear/planar grain structure generally behave transversely-isotropic. In the latter case, the UCS of intact rocks must be determined by a set of laboratory tests on the oriented rock samples. There are some empirical relations available to describe the strength of these rocks. Though characterization of transversely-isotropic rocks is practically a 3D problem, but these relations provide only a 2D description. In this paper, a method is proposed to provide a 3D description of UCS of transversely-isotropic rocks. By means of this formulation, one can determine UCS along with any arbitrary spatial direction. Also, a representative illustration of UCS is proposed in the form of contour-plots on a lower hemisphere Stereonet. The method is applied to an actual case study from the Kanigoizhan dam site located in the Kurdistan Province (Iran). Application of the proposed method to the phyllite rocks of this site show that the direction perpendicular to the dam axis exhibits the most anisotropic behavior. Hence, it is essential to take the strength anisotropy into account during the relevant analysis. The results obtained, together with the statistical variation of UCS, provide a practical approach to select the proper values of UCS according to the scope of the analysis.