Rock Mechanics
Aref Jaberi; Shokroallah Zare
Abstract
Unlike the mechanical properties of intact rock, which can be obtained on a laboratory scale, estimating the mechanical properties of the jointed rock mass is very difficult due to the presence of different joints and the complexity of the joints. Therefore, to calculate the mechanical parameters of ...
Read More
Unlike the mechanical properties of intact rock, which can be obtained on a laboratory scale, estimating the mechanical properties of the jointed rock mass is very difficult due to the presence of different joints and the complexity of the joints. Therefore, to calculate the mechanical parameters of the jointed rock mass and use the continuous media theory of the jointed rock mass, it is necessary to calculate the Representative Element Volume (REV) of the rock mass. In this study, the Discrete Element Method (DEM) and the mechanical index of strength were used to investigate the effect of persistent and non-persistent joint angles, as well as model size on the REV in x, y, and z directions. The numerical results showed that by changing the joint angles and side length, both the strength and the REV of the rock mass were affected. The maximum representative side length for the persistent joint in the x and z directions occurred at angles of 60° and 75°, respectively. The minimum strength was obtained for joints in the x and z directions at a 45° angle. Finally, the REV for persistent and non-persistent joints is calculated as 10*0.5*8m and 4*0.5*4m, respectively.
A. Abdollahipour; M. Fatehi Marji; A. R. Yarahmadi Bafghi; J. Gholamnejad
Abstract
Hydraulic fracturing (HF), as a stimulation technique in petroleum engineering, has made possible the oil production from reservoirs with very low permeability. The combination of horizontal drilling and multiple HF with various perforation angles has been widely used to stimulate oil reservoirs for ...
Read More
Hydraulic fracturing (HF), as a stimulation technique in petroleum engineering, has made possible the oil production from reservoirs with very low permeability. The combination of horizontal drilling and multiple HF with various perforation angles has been widely used to stimulate oil reservoirs for economical productions. Despite the wide use of HF, there are still ambiguous aspects that require more investigation. Therefore, optimizing the geometry of the initial fractures using numerical methods is of high importance in a successful HF operation. Different geometrical parameters of the initial HF cracks including patterns, spacings, crack lengths, and perforation phase angles were modeled using the higher order displacement discontinuity method (HODDM) in horizontal and vertical oil wells. Several well-known issues in HF such as crack interference and crack arrest were observed in certain patterns of the HF cracks. Also the best possible arrangements of the HF cracks were determined for a better production. The results obtained were verified by the in-situ measurements existing in the literature. In addition, the best perforation phase angle in vertical wells was investigated and determined.
Exploration
Kamran Mostafaei; Mohammad Nabi Kianpour; Mahyar Yousefi; Meisam Saleki
Abstract
Discrimination of geochemical anomalies from background is a challenge in that elemental dispersion patterns are affected by a variety of geological factors, which vary from one to another area. While statistical and fractal methods are commonly employed for anomaly detection, they struggle with selecting ...
Read More
Discrimination of geochemical anomalies from background is a challenge in that elemental dispersion patterns are affected by a variety of geological factors, which vary from one to another area. While statistical and fractal methods are commonly employed for anomaly detection, they struggle with selecting optimal thresholds. This study proposes the Grey Wolf Optimizer (GWO) algorithm as a novel approach for identifying the optimal boundary between anomalies and background. Stream sediment geochemical data from a copper-mineralized area of the Sarduiyeh-Baft sheets in southeast Iran were utilized for analysis. The Geochemical Mineralization Probability Index (GMPI) was first calculated for Cu-Au, Mo-As, Pb-Zn, and porphyry distributions. Subsequently, fractal methods were used to identify anomalous populations within each GMPI. The GWO algorithm was then applied to these distributions to determine the optimal thresholds. Risk analysis, calculated as the ratio of covered copper occurrences to the covered area, revealed superior reliability for the GWO-derived limit compared to those obtained using fractal methods. For porphyry GMPI values, while the fractal reliability indices are 0.127, 0.44, and 0.5, the GWO limit achieved a value of 0.66. Risk analysis for Cu-Au distribution also caused more desired result for GWO limit rather that fractal ones. This demonstrates the enhanced performance and superior reliability of the GWO algorithm for optimizing anomaly detection thresholds in GMPI data.
M. Mokhtarian Asl; J. Sattarvand
Abstract
Production planning of an open-pit mine is a procedure during which the rock blocks are assigned to different production periods in a way that leads to the highest net present value (NPV) subject to some operational and technical constraints. This process becomes much more complicated by incorporation ...
Read More
Production planning of an open-pit mine is a procedure during which the rock blocks are assigned to different production periods in a way that leads to the highest net present value (NPV) subject to some operational and technical constraints. This process becomes much more complicated by incorporation of the uncertainty existing in the input parameters. The commodity price uncertainty is among the most significant factors, whose effects cannot be mitigated through further exploration or investigation. The present work introduced a new approach for integration of the commodity price uncertainty into long-term production planning of open-pit mines. The procedure involves solving the problem by the integer programming method based on a series of economic block models that are realized based on the sampled prices from commodity price distribution function using the median Latin hypercube sampling method. The results obtained showed that the new methodology is able to reduce the risks and the net present value of the new approach at a confidence level 80% more than the conventional methods.
M. Doustmohammadi; A. Jafari; O. Asghari
Abstract
Water inflow is one of the most important challenges in the underground excavations. In addition to inducing working conditions and environmental problems, it decreases the stability and quality of the surrounding rocks. The direct method of measuring rock mass hydraulic conductivity consists of drilling ...
Read More
Water inflow is one of the most important challenges in the underground excavations. In addition to inducing working conditions and environmental problems, it decreases the stability and quality of the surrounding rocks. The direct method of measuring rock mass hydraulic conductivity consists of drilling the boreholes and observing the rate of fluid lost in the boreholes. Applying this method is still problematic due to the depth of underground spaces, and also the groundwater level covering them. Therefore, many researchers have tried to predict the water inflow indirectly. This paper attempts to predict the groundwater conditions in the Beheshtabad tunnel (in Iran) using the fuzzy inference system based on the datasets acquired from the preliminary exploration studies. 250 datasets for the Beheshtabad tunnel were used out of which, 200 datasets were used to develop the model and 50 were used to validate the results obtained. 90% accuracy was obtained through comparing the fuzzy estimation and actual groundwater conditions. The proposed model can be used with much less degree of complexity for prediction of the groundwater conditions as well as decreasing the overall costs of the exploration measurements, and due to these characteristics, it is applicable for most users.
R. Vahedi; B. Tokhmechi; M. Koneshloo
Abstract
We use a multi-resolution analysis based on a wavelet transform to upscale a 3D fractured reservoir. This paper describes a 3D, single-phase, and black-oil geological model (GM) that is used to simulate naturally-fractured reservoirs. The absolute permeability and porosity of GM is upscaled by all the ...
Read More
We use a multi-resolution analysis based on a wavelet transform to upscale a 3D fractured reservoir. This paper describes a 3D, single-phase, and black-oil geological model (GM) that is used to simulate naturally-fractured reservoirs. The absolute permeability and porosity of GM is upscaled by all the possible combinations of Haar, Bior1.3, and Db4 wavelets in three levels of coarsening. The applied upscaling method creates a non-uniform computational grid, which preserves its resolved structure in the near-well zones as well as in the high-permeability sectors but the data are scaled up in the other regions. To demonstrate the accuracy and efficiency of the method, the values for the oil production rate, mean reservoir pressure, water cut, and total amount of water production are studied, and their mean error is estimated for the upscaled models. Finally, the optimized model is selected based on the computation time and accuracy value.
A. Siamaki; H. Bakhshandeh Amnieh
Abstract
A considerable amount of energy is released in the form of shock wave from explosive charge detonation. Shock wave energy is responsible for the creation of crushing and fracture zone around the blast hole. The rest of the shock wave energy is transferred to rock mass as ground vibration. Ground vibration ...
Read More
A considerable amount of energy is released in the form of shock wave from explosive charge detonation. Shock wave energy is responsible for the creation of crushing and fracture zone around the blast hole. The rest of the shock wave energy is transferred to rock mass as ground vibration. Ground vibration is conveyed to the adjacent structures by body and surface waves. Geological structures like faults, fractures, and fillings play important roles in the wave attenuation. Studying the mechanism of ground wave propagation from blasts gives a better understanding about the stress wave transmission and its effect on the near structures. In this research work, the stress wave transmissions from discontinuities and fillings were evaluated using a field measurement and a Universal Distinct Element Code (UDEC). A single-hole blast was conducted in the Kangir dam, and the resulting vibrations were measured in many points before and after the faults. Numerical simulation shows the effects of geo-mechanical properties of fillings on the reflection and refraction rate of the stress wave. There are more energy reflections in the rock boundaries and soil fillings, and more energy is absorbed by soil fillings compared with rock fillings. Furthermore, there is a close correlation between the ground vibration records for the Kangir dam and the numerical results. The maximum relative error between the actual records and the simulated ones was found to be 18.5%, which shows the UDEC ability for the prediction of blast vibrations.
M. Mohammady Oskouei; S. Babakan
Abstract
This work aims to extract the mineralogical constituents of the Lahroud Hyperion scene situated in the NW of Iran. Like the other push-broom sensors, Hyperion images suffer from spectral distortions, namely the smile effect. The corresponding spectral curvature is defined as an across-track wavelength ...
Read More
This work aims to extract the mineralogical constituents of the Lahroud Hyperion scene situated in the NW of Iran. Like the other push-broom sensors, Hyperion images suffer from spectral distortions, namely the smile effect. The corresponding spectral curvature is defined as an across-track wavelength shift from the nominal central wavelength, and alters the pixel spectra. The common “column mean adjusted in MNF space” method was employed in this work to improve the processing accuracy by minimizing the smile effect before carrying out the atmospheric and topographical corrections. The mineral distributions were mapped by applying the standardized hyperspectral processing methodology developed by analytical imaging and geophysics (AIG). The spectral unmixing of the data resulted in the identification of five indicative minerals including natrolite, opal, analcime, kaolinite, and albite; and their spectra were employed for the generation of their distribution maps. Comparison of the results of the data processing with and without smile correction indicated a better classification performance after the smile correction. Quantitative validation of the final mineralogical map was performed using the 100 k geological map and reports of the region. Therefore, the coverage of the extracted minerals were investigated regarding the location of the lithological units in ArcGIS that implies a high coincidence. The mineral distributions in the final map show a high consistency with the geological map of the studied area, and thus it could be utilized successfully to reveal the mineralization trend in the region.
M. Nikkhah; Seyed S. Mousavi; Sh. Zare; O. Khademhosseini
Abstract
The joints between segmental rings can withstand a certain amount of bending moment as well as axial and shear forces. Generally, in the structural analysis of tunnel segmental lining, the joints can be modeled as elastic hinges or rotational springs, and their rigidity should be demonstrated in terms ...
Read More
The joints between segmental rings can withstand a certain amount of bending moment as well as axial and shear forces. Generally, in the structural analysis of tunnel segmental lining, the joints can be modeled as elastic hinges or rotational springs, and their rigidity should be demonstrated in terms of the rigidity of the joints or their rotational stiffness. Therefore, the bending moment acting on the tunnel lining is reduced. Hence, the tunnel designers are free to choose a lining with a lower cost. In this research work, especially considering the joints, the structural analysis of the segmental lining with variation in the flexural stiffness of the joints ( ), soil resistance coefficient ( ), number of segmental lining joints, and joint arrangement of segmental lining were carried out by the Force-Method equations. The imposed bending moment and axial forces were computed based on the Beam-Spring method, which is widely used to analyze the internal forces of segmental lining, and compared them with the results of the Force-Method equations. Then the effects of joint arrangement patterns and joint rotational spring stiffness on the results of the Beam-Spring analysis were evaluated. Finally, the optimum characteristics of the reinforced concrete segmental lining design were evaluated using the interaction diagram of bending moments and axial forces. The results obtained showed that the presented pattern for the segmental lining at the Chamshir tunnel was imposed against the external pressures on the segmental lining with an acceptable safety factor.
H. Molayemat; F. Mohammad Torab
Abstract
Coalbed methane (CBM) plays an important role in coal mining safety and natural gas production. In this work, The CBM potential of B2 seam in Parvadeh IV coal deposit, in central Iran, was evaluated using a combination of local regression and geostatistical methods. As there were 30 sparse methane sampling ...
Read More
Coalbed methane (CBM) plays an important role in coal mining safety and natural gas production. In this work, The CBM potential of B2 seam in Parvadeh IV coal deposit, in central Iran, was evaluated using a combination of local regression and geostatistical methods. As there were 30 sparse methane sampling points in the Parvadeh IV coal deposit, no valid variogram was achieved for the methane content. A multivariate adaptive regression splines (MARS) model was used to reproduce the methane content data based on seam depth, thickness, and ash content. The MARS model results were used in ordinary kriging to estimate the methane content in all mine blocks. A combination of MARS modeling and ordinary kriging in CBM studies is introduced, for the first time, in this paper. The results obtained show that high methane zones are located in the central and south western parts of the deposit. The in situ CBM potential varies from 6.0 to 16.1 m3/t, and it was estimated to be 1.39 billion m3 at the average depth of 267 m in an area of 86.55 km2. Although this volume is remarkable, little is known as how much of this resource is actually producible. Consequently, high methane-bearing zones are highly recommended for further studies as a source of natural gas for extraction and reducing the hazards and explosion risks of underground coal mining.
A. Majdi; M. Yazdani
Abstract
The hydraulic jacking refers to the process of crack growth of the pre-existing joints in the rock mass under grout pressure above the minimum in-situ stress. Thus it is essential to understand the resistance behavior of the joints and maximum grout pressure. This paper describes a novel method for determining ...
Read More
The hydraulic jacking refers to the process of crack growth of the pre-existing joints in the rock mass under grout pressure above the minimum in-situ stress. Thus it is essential to understand the resistance behavior of the joints and maximum grout pressure. This paper describes a novel method for determining the hydraulic jacking occurrence in anisotropic rock mass based on the principle of fracture mechanics. This method is established on three stage developments: developing an equation in order to calculate the equivalent stress intensity factor at the joint tip, determining the fracture toughness by employing the Brazilian disc test with a loading rate of 0.1 MPa/s on the rock cored samples, and assessing the stability of joints using the maximum tangential stress criterion. By comparing the joint stress intensity factor and fracture toughness in the direction of rock anisotropy, the joint stability is evaluated. Then the maximum allowable grout pressure is analytically formulated as a function of fracture toughness in order to avoid the unwanted deformations in the joints (i.e. jacking) during grouting. In order to validate the proposed method, the data obtained from the boreholes used to construct water curtain at the Sanandaj Azad Dam in phyllite rocks are analyzed. Finally, it is concluded that the growth and expansion of the joints due to the instability under grout pressure leads to an increased cement take and the occurrence of hydraulic jacking. In addition, the proposed equation for computing maximum allowable grout pressure provides an acceptable agreement with the existing empirical rules and the results of the field data.
Mohammadhosein Dehghani Firoozabadi; Mohammad Fatehi Marji; Abolfazl Abdollahipour; Alireza Yarahamdi Bafghi; Yousef Mirzaeian
Abstract
In this work, an effective methodology is introduced for simulation of the crack propagation in linear poroelastic media. The presence of pores and saturated cracks that can be accompanied by fluid flow makes the use of poroelastic media inevitable. In this work, involvement of the time parameter in ...
Read More
In this work, an effective methodology is introduced for simulation of the crack propagation in linear poroelastic media. The presence of pores and saturated cracks that can be accompanied by fluid flow makes the use of poroelastic media inevitable. In this work, involvement of the time parameter in crack propagation is of particular importance. The order of doing the work is such that first, derives the fundamental solutions of a poroelastic higher order displacement discontinuity method (PHODDM). Then will be provided a numerical formulation and implementation for PHODDM in a code named linear element poroelastic DDM (LEP-DDM). Analytical solutions use different times to check the correctness and validity of the proposed solution and the newly developed code. The numerical results show a good agreement and coordination with the analytical results in time zero and 5000 seconds . The code is able to pursue crack-propagation in time and space. This topic is introduced and shown in an example.
Exploitation
Pouya Nobahar; Yashar Pourrahimian; Roohollah Shirani Faradonbeh; Fereydoun Mollaei Koshki
Abstract
Mineral reserve evaluation and ore type detection using data from exploratory boreholes are critical in mine design and extraction. However, preparing core samples and conducting chemical and physical tests is a time-consuming and costly procedure, slowing down the modeling process. This paper presents ...
Read More
Mineral reserve evaluation and ore type detection using data from exploratory boreholes are critical in mine design and extraction. However, preparing core samples and conducting chemical and physical tests is a time-consuming and costly procedure, slowing down the modeling process. This paper presents a novel Deep Learning (DL)-based model to recognize the types of kaolinite samples. For this purpose, a dataset containing the images of drilled cores and their types determined from conventional chemical and physical analyses was used. Eight Convolutional Neural Network (CNN) topologies based on individual features were developed, named A, B, C, D, E, F, G, and H. Six of the eight proposed CNN topologies described above had accuracy below 80%, whereas two of them, model A and H, had higher accuracy than other topologies. Due to their similarity in results, both of them analyzed deeply. Model A was more efficient, with 90% accuracy, than model B, with 84% accuracy. Furthermore, the class detection performance of model A was further evaluated using different indices, including precision, recall, and F1-score, which resulted in values of 92%, 92%, and 90%, respectively, which are acceptable accuracies to identify the type of samples when using this approach on six different types of kaolinite.
Rock Mechanics
Hamed Farajollahi; Mohammad Mohammadi; Mohammad Hossein Khosravi
Abstract
A better understanding of rock mass behavior is an essential part of the design and construction of underground structures. Any improvement in the understanding of the behavior of rock mass will facilitate the improvement of the design in terms of the safety of the working environment, long-term safety ...
Read More
A better understanding of rock mass behavior is an essential part of the design and construction of underground structures. Any improvement in the understanding of the behavior of rock mass will facilitate the improvement of the design in terms of the safety of the working environment, long-term safety of the structure, environmental effects, and sound management of public or private resources. Thus, in step one in this paper the experience gained from the application of the GDE (Geo Data Engineering) multiple graph approach for rock mass classification and assessment of its behavior through the course of excavation of the Alborz tunnel is presented. The predicted hazards are compared with the experienced problems and suggestions are given to be considered in future works of tunnel construction. In step two, the GDE multiple graph approach is compared to the rock mass behavior types proposed by Palmstrom & Stille (2007) in terms of the continuity of rock mass. The result of this comparison together with the data obtained from rock mass classification in the Alborz tunnel are used to develop a system that determines the applicability of the rock bolt supporting factor (RSF) in different rock mass behavior classes.
Exploration
mobin saremi; Abbas Maghsoudi; Reza Ghezelbash; mahyar yousefi; Ardeshir Hezarkhani
Abstract
Mineral prospectivity mapping (MPM) is a multi-step and complex process designed to narrow down the target areas for exploratory activities in subsequent stages. To pinpoint promising zones of porphyry copper mineralization in the Varzaghan district, NW Iran, various exploration evidence layers were ...
Read More
Mineral prospectivity mapping (MPM) is a multi-step and complex process designed to narrow down the target areas for exploratory activities in subsequent stages. To pinpoint promising zones of porphyry copper mineralization in the Varzaghan district, NW Iran, various exploration evidence layers were employed in alignment with the conceptual model of these deposits. These layers encompass fault density, proximity to intrusive rocks, multi-element geochemical anomalies, and distances to phyllic and argillic alterations. The geochemical anomaly maps, recognized as the most effective layers, were generated through staged factor analysis (SFA) and the geochemical mineralization probability index (GMPI). Other layers were weighted using a logistic function, and their values were transformed into 0 -1 interval. Ultimately, to integrate the weighted layers, the fuzzy gamma operator and the geometric average method were applied. The normalized density index and prediction-area (P-A) plot were employed to evaluate the MPM models. The findings indicate that the developed models possess considerable validity and can be effectively utilized for planning future exploration endeavors.
Mineral Processing
Seyyed Mohsen Zamzami; Javad Vazifeh Mehrabani
Abstract
In this research, solid phase settling process from the liquid phase were optimized simultaneously on the different responses, using the response surface methodology (RSM). The effect of solid percentage, flocculant dosage, temperature, and pulp pH were evaluated on the responses of solid settling velocity, ...
Read More
In this research, solid phase settling process from the liquid phase were optimized simultaneously on the different responses, using the response surface methodology (RSM). The effect of solid percentage, flocculant dosage, temperature, and pulp pH were evaluated on the responses of solid settling velocity, water turbidity, viscosity and density of settled pulp. The results showed that by increasing the flocculant dosage from 0.5 to 3.5 g/ton, settled pulp viscosity decreases from 49.05 cSt to 17.54 cSt. The higher values of pulp pH as well as low amount of solid percentage resulted in high water turbidity, which shows the lack of contact between flocs and suspended particles. The results indicated that the pulp solid percentage and the flocculants dosage are the most significant parameters on the responses. Optimum test conditions were obtained in industrial mode by using 5 g/t flocculant, solid percentage 23.96%, pH=7.5 temperature of the pulp 21.5°C in which condition, settling rate, pulp viscosity, pulp density and water turbidity were predicted to be 13.23 cm/min, 5.1 cSt, 1.61 g/cm3 and 15.7 NTU respectively. Repetition test in the model predicted optimum condition was carried out and verified the predicted optimized condition.
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.