Rock Mechanics
Masoud Yazdani; Mohammad Fatehi Marji; Mehdi Najafi; Manouchehr Sanei
Abstract
Around 70% of the world's hydrocarbon fields are situated in reservoirs containing low-strength rocks, such as sandstone. During the production of hydrocarbons from sandstone reservoirs, sand-sized particles may become dislodged from the formation and enter the hydrocarbon fluid flow. Sand production ...
Read More
Around 70% of the world's hydrocarbon fields are situated in reservoirs containing low-strength rocks, such as sandstone. During the production of hydrocarbons from sandstone reservoirs, sand-sized particles may become dislodged from the formation and enter the hydrocarbon fluid flow. Sand production is a significant issue in the oil industry due to its potential to cause erosion of pipes and valves. Separating grains from oil is a costly process. Oil and gas companies are motivated to reduce sand production during petroleum extraction. Hydraulic fracturing is one of the parameters that can influence sand production. However, understanding the complex interactions between hydraulic fracturing mechanisms and sand production around wellbores is critical for optimizing reservoir recovery and ensuring the integrity of production wells. This article explores the integrated simulation approach to model hydraulic fracturing processes and assess their effects on sand production. Two-dimensional models were created using the discrete element method in PFC2D software for this research. The fractures' length in the models varies based on the well's radius. The angle between two fractures at 90 and 180 degrees to each other was also modeled. In the first case, the length of the fracture is less than the radius of the well, in the second case, the values are equal and finally, the fracture length is assumed to exceed the well radius. The calibrated and validated results demonstrate the change in sand production rate in comparison to the unbroken state.
Mineral Processing
Ali Nikouei Mahani; Mohammad Karamoozian; Mohammad Jahani Chegeni; Mohammad Mahmoodi Meymand
Abstract
Generally, mineral processing plants generate a large quantity of waste in the form of fine particles. The flotation speed of mineral microbubbles by coarse bubbles is dramatically higher than that of individual particles. The advantage of microbubbles is due to the increase of binding efficiency of ...
Read More
Generally, mineral processing plants generate a large quantity of waste in the form of fine particles. The flotation speed of mineral microbubbles by coarse bubbles is dramatically higher than that of individual particles. The advantage of microbubbles is due to the increase of binding efficiency of conventional bubbles with fine particles coated with microbubbles. Here, the focus is on reducing chemicals consumption and improving recovery. After preparing a representative sample, XRF, XRD, and mineralogical analyses were performed. Then 50 experiments were selected by experimental design using the response surface method (RSM), and in the form of central Composite design (CCD) by (design expert) DX 13 software. The interactions of collector consumption, frother agent, pH, particle size, and solid percentage were investigated, and 25 experiments using typical flotation and without nano-microbubbles and others with nano-microbubbles were conducted. The laboratory standard limit of the collector used in the pilot plant of the Sarcheshmeh Copper copper complex is 40 g/t (25 g/t of C7240 plus 15 g/t of Z11). Here, by consuming 20 g/t of collector in the absence of nanomicrobubbles, a recovery of 79.96% and in the presence of nanomicrobubbles, a recovery of 80.07% was obtained, that is a 50% reduction in collector consumption and a 0.11% increase in recovery was observed. Also the laboratory standard limit of frother used in the pilot plant of Sarcheshemeh Copper Complex is 30 g/t (15 g/t of MIBC plus 15 g/t of A65). Here, by using 10 g/t of frother in the absence of nanomicrobubbles, a recovery of 78.12%, and in the presence of nanomicrobubbles, a recovery of 82.05% was obtained. In other words, a decrease of 66.6% in the consumption of frother and an increase of 1.93% in recovery was observed.
M. Ataei; F. Sereshki
Abstract
Like most limestone mines, which produce the raw materials required for cement companies, the transportation cost of the raw materials used in the Shahrood Cement Company is high. It has been tried to build the crushing and grinding plant close to the mine as much as possible. On the other hand, blasting ...
Read More
Like most limestone mines, which produce the raw materials required for cement companies, the transportation cost of the raw materials used in the Shahrood Cement Company is high. It has been tried to build the crushing and grinding plant close to the mine as much as possible. On the other hand, blasting has harmful effects, and the impacts of blast-induced damages on the sensitive machinery, equipment, and buildings are considerable. In such mines, among the blasting effects, blast-induced vibrations have a great deal of importance. This research work was conducted to analyze the blasting effects, and to propose a valid and reliable formula to predict the blast-induced vibration impacts in such regions, especially for the Shahrood Cement Company. Up to the present time, different indices have been introduced to quantify the blast vibration effects, among which peak particle velocity (PPV) has been widely considered by a majority of researchers. In order to establish a relationship between PPV and the blast site properties, different formulas have been proposed till now, and their frequently-used versions have been employed in the general form of , where W and D are the maximum charge per delay and the distance from the blast site, respectively, and , , and describe the site specifications. In this work, a series of tests and field measurements were carried out, and the required parameters were collected. Then in order to generalize the relationship between different limestone mines, and also to increase the prediction precision, the related data for similar limestone mines was gathered from the literature. In order to find the best equation fitting the real data, a simple regression model with genetic algorithm was used, and the best PPV predictor was achieved. At last, the results obtained for the best predictor model were compared with the real measured data by means of a correlation analysis.
Rock Mechanics
M. B. Eslami Andargoli; K. Shahriar; A. Ramezanzadeh; K. Goshtasbi
Abstract
During the recent decades, the design and construction of underground spaces into rock salt have been particularly regarded for storing petroleum fluids, natural gas, and compressed air energy, and also for disposing nuclear and chemical wastes. The rock salt hosting such spaces will be subjected to ...
Read More
During the recent decades, the design and construction of underground spaces into rock salt have been particularly regarded for storing petroleum fluids, natural gas, and compressed air energy, and also for disposing nuclear and chemical wastes. The rock salt hosting such spaces will be subjected to various types of monotonic/cyclic, short-term/long-term stresses during the construction and/or operation phases. On this basis, it is necessary to investigate the mechanical behavior of the rock salt under the effects of various monotonic short-term/long-term stresses. Out of the most important factors affecting the creep behavior of rock salt are the composition of minerals and size of the crystals comprising the rock salt, humidity, temperature, time, loading scheme, loading rate, strain rate, and loading period. In the present research work, a loading scheme and a loading period were considered. On this basis, in order to achieve a true understanding of the creep behavior of rock salt, it is necessary to determine the creeping coefficients via laboratory tests. Thus, twenty cylindrical (length to diameter ratio > 2) specimens of rock salt were prepared for conducting the creep tests. Two stepwise short-term creep tests (each at three stress levels, namely 4.4, 10.1, and 11.9 MPa, and 7.5, 12, and 17 MPa, respectively) and eighteen long-term creep tests (at six stress levels, namely 5.5, 7.5, 10, 12, 14, and 18 MPa) were conducted. Then, first, the creep coefficients were determined according to the Lubby 2 constitutive model. These coefficients were adjusted using the results of the creep tests. Afterwards, a creep experimental model was presented using linear and nonlinear regression of the creep test data. For validation of the results obtained, both the adjusted Lubby 2 constitutive model and the proposed experimental model were compared with the results obtained for the creep tests. Both models had fairly good agreements with the data for the creep tests at a determination factor of about 93%.
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
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
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
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.
