Search Documents

By D. F. Malan, R. P. Els

The limit equilibrium model, used in displacement discontinuity codes, is a popular method to simulate pillar failure. This paper investigates the use of physical modelling to calibrate this model. For the experiments, an artificial pillar material was prepared and cubes were poured using the standard 100 mm × 100 mm civil engineering concrete moulds. The friction angle between the cubes and the platens of the testing machine was varied by using soap and sandpaper. Different modes of failure were observed depending on the friction angle. Of interest is that significant loadshedding was recorded for some specimens which visually remained mostly intact. This highlights the difficulty of classifying pillars as failed or intact in underground stopes where spalling is observed. The laboratory models enabled a more precise calibration of the limit equilibrium model compared to previous attempts. Guidelines to assist with calibration of the model are given in the paper. The limit equilibrium model appears to be a useful approximation of the pillar failure as it could simulate the stress-strain behaviour of the laboratory models.

May 3, 2023

By W. C. Joughin

A risk evaluation model has been developed, that enables the quantification of expected injuries and economic losses resulting from rockfalls associated with a support system in a given geotechnical environment. This enables the comparison of different support systems based on the expected frequency of injuries and the total cost of the support system, including the expected economic losses. The model takes into consideration the variability in joint orientations and joint strength within the geotechnical environment and the variability and quality of support systems. The quantification of the benefits of other rockfall risk-mitigating strategies, such as the effectiveness of barring or monitoring, are also incorporated in the model. The risk model is described in this paper and the validation of the model and case studies are described in Part II.

Jan 1, 2012

By J. Labuschagne

The primary factor contributing to the profitability of a platinum mine and its associated concentrator is the time taken from the first hoist of material from underground to the first batch of concentrate shipped. To minimize the time it would take to produce concentrate without stopping because of a lack of feed material a project was undertaken to integrate the mining model with a concentrator production schedule () in the most capital-effective manner. A secondary goal was to ensure detailed stockpiling requirements, as the concentrator is in an area where ?visual pollution? is to be avoided. To achieve this a dynamic stockpile model was created utilizing the mine production schedule resolved into a daily production figure. Four different concentrator development models were then proposed. The difference between the mine production and concentrator consumption was integrated over the life of the mine to provide the accumulated stockpile tonnages for both Merensky and UG2 ores. The stockpile levels were never allowed to fall below zero to prevent a situation where the concentrator was operational without feed material,. When a development option was selected, the standard metallurgical performance calculations were used to derive a concentrate production schedule, that was then used as a basis for commercial concentrate sale negotiations.

Jan 1, 2012

By S. Mohanta, B. C. Melkap

"SynopsisDry beneficiation of coal by air dense medium fluidized bed is an emerging trend. It is widely believed that the particle size of the medium has a significant effect on the separation efficiency. This investigation demonstrates that medium particle size has a major effect on separation efficiency. Experimental results show that different size fractions of the same feed respond differently to the same size fraction of medium solids. Furthermore, a particular size fraction of feed coal responds differently with different size fractions of medium solids. The Ep values and overall metallurgical performance parameters, obtained from experimental results, clearly indicate the superior performance of an air dense medium fluidized bed separator when using –150+106 µm magnetite powder as the fluidizing medium. These observations reinforce the importance of sizing the medium particle size for the air dense medium fluidized bed.IntroductionIn an air dense medium fluidized bed (ADMFB), a fluidizing bed of a particular density is created by suspending solid particles in an upward air flow. A feed containing coal particles of different densities is introduced into this fluidizing bed. The particles with relative density less than the bed density float on the top of the bed, while particles denser than the bed density sink through it. Since there is a general correlation between ash content and density of coal particles, it is possible to remove ash-forming impurities from a raw coal to the required level by regulating the density of the fluidizing bed. Recently, the ADMFB process has received significant attention world-wide for its specific advantages of low operating cost and less restrictive for small capacities (Luo and Chen, 2001; Houwelingen and de Jong, 2004; Dwari and Hanumantha Rao, 2007; Sahu et al., 2009; Mohanta et al., 2012). This process is not recognized as a well-established industrial practice because of the availability of efficient wet processing technology. However, with the revival of interest in dry beneficiation of coal in general, and ADMFB in particular, there is a concerted effort to improve the efficiency of the ADMFB process, and successive improvements have made the process a competitive and viable option for beneficiating coal."

Jan 1, 2015

By J. Loubser

The concept of caring is in most cases confused with material needs and not psychological needs. By giving employees a job, most managers think that employees should be grateful and return the favour by being loyal, productive and working safely. Other managers hold the view that employees can be lucky to have a job. To show and demonstrate that management are caring for the safety of the workforce requires involvement, dedication and commitment. Obtaining the trust and respect of the workforce should be the main thrust of any organization because within trust and respect lie the secret of optimal safety performance. In the last 70 years, all the fatalities and serious injuries that occurred in the mining industry can be traced back to one thing, non-caring. The same safety non-compliances are still observed, reported and planned against as 70 years ago, and strange as it may sound, the same injuries occur the same way. The question could be asked: Does the mining industry need more risk assessments? More codes of conduct? More engineering. Or does the mining industry need more caring?

Jan 1, 2010

By M. Royker

Gold has been demonstrated as a possible catalyst for oxidation reactions. Some evidence for a possible promotion effect of platinum has also been recorded. The influence of platinum as promoter for Au/?-Al2O3 prepared via anionic ion-exchange for the oxidation of glycerol was investigated in a batch reactor at 60°C. It is inferred that the addition of platinum reduces the catalytic activity and the rate of deactivation, resulting in an overall higher final conversion of glycerol with increasing platinum loading. The addition of platinum to the catalyst favours the formation of the desired product, glyceric acid.

Jan 1, 2012

By M. Stachowiak, M. Nowicki, K. Walas, A. Skoczylas, P. Stefaniak

In this paper the authors present a walking robot dedicated to the inspection of a belt conveyor’s transportation system in an underground mine. The developed robot can independently diagnose all of the key rotating elements of the belt conveyor (pulleys, bearings, drums, gears, engines) based on vibration and acoustic signals as well as typical RGB (red, green, blue) and thermal images. The article presents key diagnostic procedures applied to real industrial data obtained during robot testing.The work was developed as part of the ‘THING–subTerranean Haptic INvestiGator’ project. The authors presented the main assumptions of the research methodology developed for the needs of technical diagnostics of conveyor elements typical of inspection work performed by mining staff. Based on industrial data, the use of algorithms for elements such as rollers, drive stations, and the conveyor belt were shown. The main assumption of the research is to develop a robot adapted to the conditions of a deep mine, capable of carrying out an autonomous inspection of mining infrastructure. The conducted tests and the developed fault detection algorithms were aimed at preparing the robot for underground tests in the operating conditions of the mine.

Sep 1, 2021

By O. T. Konina, V. A. Masloboev, A. V. Svetlov, E. A. Krasavtseva, D. V. Makarov, A. A. Goryachev

Dusting phenomena have been examined at the tailings storage facility of Apatit JSC Kirovsk, a mining operation located in the far north of the Russian Federation. The material composition of the tailings was identified. Existing dust suppression methods were examined. Field observations were undertaken and official data was obtained from the Murmansk Administration for Hydrometeorology and Environmental Monitoring, and analysed. A database was built for forecasting dusting according to critical scenarios in order to estimate the peak MPCs of suspended solids in the atmospheric air. Mitigation measures were formulated and recommended to the mining operation. Keywords: Mineral dust particles, apatite-nepheline ore concentration tailings, suspended solids MPC, climate change

Jan 1, 2020

By Arvi Parbo

I would like to use the occasion to reflect on how the world environment for the minerals industry is undergoing substantial change, and how this change may affect the industry. I have given my talk the title 'The Quest For Competitive Edge', because the net effect of the changes taking place is that competition is increasing. To set the scene, let us look back at the situation which has existed in our industry in the past.

Jan 1, 1993

By E. Marè, D. R. Groot, R. Engelbrecht, T. Coetsee

Hydrogen peroxide is often used as an oxidising agent in experimental waste electronic equipment (WEE) leaching work. However, the high cost of peroxide is a negative factor in its consideration in further process development. A few studies have reported on the efficient application of peroxide addition schedules to ensure optimum reagent consumption. Large quantities of unleached substrate material contained in WEE feed material result in a large dead volume of material in the leaching reactor. This factor necessitates fast reaction rates to enable the use of a high feed material throughput reactor to enhance overall project viability via hydrometallurgical processing. In this study the effects of hydrogen peroxide addition schedules on base metal leaching efficiency and solution temperatures were tested at laboratory scale. The results showed that overall peroxide consumption can be optimised to lower levels by matching peroxide additions to the particular feed material’s leaching behaviour. A simplified integrated model for peroxide decomposition is used to explain its effect on solution temperatures and leaching rates. Keywords: Peroxide consumption, waste electronic equipment, copper leaching, integrated model

Jan 1, 2020

By S. Krishnannair

Singular Spectrum Analysis (SSA) is an effective data-adaptive multimodal tool for monitoring chemical processes. However, monitoring multiscale signals using principal component analysis obtained from the multilevel decomposition of process signals using basic SSA assume the upper control limits of monitoring charts (" ) to have a Gaussian distribution. This assumption deteriorates the performance of basic SSA when applied to the monitoring of nonlinear processes. To address the issue in this paper, the kernel density estimation (KDE) was used to estimate upper control limits of monitoring charts for basic SSA based nonlinear process monitoring techniques. The monitoring performance of the resulting basic SSA-KDE approach was then compared with a basic SSA approach of which upper control limits of the monitoring charts were determined based on a Gaussian distribution. Application on nonlinear dynamic simulated process and base metal flotation plant chemical processes showed that the SSA-KDE approach provided better overall performance than the SSA approach with Gaussian assumption-based upper control limits for its monitoring charts. Keywords: Fault detection, singular value decomposition, singular spectrum analysis, kernel density estimation

Jan 1, 2020

By I. James, R. J. Durrheim, M. S. D. Manzi, O. Valeshin, J. Lepine, L. V. Socco, C. Colombero, A. Malehmir, M. K. Rapetsoa, M. Sihoyiya

The FUTURE project aims to enhance seismic surveying in South African mining by integrating advanced acquisition technologies. In this paper, we detail surface and underground seismic surveys at the South Deep gold mine as part of the project. Surface seismic acquisition used co-located broadband micro-electromechanical sensors, one and three component recorders connected to geophones, distributed acoustic sensing using a straight fibre optic cable, along with a 6 ton broadband seismic vibrator operating with 4–180 Hz linear sweeps. In addition, underground surveys employed a combination of straight, and helically wound fibre optic cables, cabled-based and nodal sensors, along with a 500 kg drophammer mounted on a skid steer bobcat. The surface and underground seismic arrays were time synchronised using a GPS-time transmitter to allow for accurate seismic data recording in a GPS-denied environment. Results show the efficacy of these techniques in acquiring high-resolution seismic data, revealing detailed subsurface geological structures and mineralisation down to 3500 m below ground surface. Surface data revealed clear seismic reflections corresponding to the Black Reef and Ventersdorp Contact Reef. An underground P-wave velocity model exhibits low velocity anomalies that are associated with fracturing, while a P-wave reflection stack section shows significant reflections and a degree of faulting below the tunnel floor that was not observed on the surface reflection seismic data. We highlight the potential of these technologies to enhance safety, efficiency, and sustainability in mining. This study sets a benchmark for future reflection seismic surveys, paving the way for safer, more efficient, deep mineral exploration and extraction practices.

Jan 12, 2025

By S. A. Villalobos, F. A. Villalobos, and L. E. Aguilera

Rockburst events have been a serious problem for many years in many mines worldwide, and in particular at El Teniente mine in Chile. El Teniente is the largest copper mine in the world, located in the Andes Cordillera where high stress levels are present due to intensing mining activity in addition to complex geology. Consequently, the study and management of the rockburst threat are necessary. In this work, the case of the Diablo Regimiento (DR) mine within El Teniente is studied. The energy capacity of dynamic support systems is determined for different tunnel geometries based on two kinetic methodologies, using data from DR. Initially, rockburst potential is determined by means of a stress analysis around different tunnel geometries through the boundary elements method. In the first methodology a yielding zone (YZ) is estimated for each excavation geometry using the finite element method FEM. The second methodology involves the definition and determination of a critical strain energy (SE) around each excavation geometry using a FEM numerical analysis. In both cases, peak particle velocity PPV is estimated by a scaling law, which is subsequently adjusted due to tunnel amplification effects. According to the results, and knowing the working energy capacity applied in DR mine, it was found that the values of energy capacity for the rock dynamic supports were better estimated by the YZ-PPV approach than by the SE approach.

Sep 1, 2022

By N. Snyders, I. Mwandawande, S. M. Bradshaw, M. Karimi

"Computational fluid dynamics (CFD) was applied to predict the average gas holdup and the axial gas holdup variation in a 13.5 m high cylindrical column 0.91 m diameter. The column was operating in batch mode. A Eulerian-Eulerian multiphase approach with appropriate interphase momentum exchange terms was applied to simulate the gas-liquid flow inside the column. Turbulence in the continuous phase was modelled using the k-_ realizable turbulence model. The predicted average gas holdup values were in good agreement with experimental data. The axial gas holdup prediction was generally good for the middle and top parts of the column, but was over-predicted for the bottom part of the column. Bubble velocity profiles were observed in which the axial velocity of the air bubbles decreased with height in the column. This may be related to the upward increase in gas holdup in the column. Simulations were also conducted to compare the gas holdup predicted with the universal, the Schiller-Naumann, and the Morsi-Alexander drag models. The gas holdup predictions for the three drag models were not significantly different. IntroductionColumn flotation is an important concentration technology that is used in the mineral processing and coal beneficiation industries. The growing interest in the use of column flotation in mineral processing has been attributed to the simpler flotation circuits and improved metallurgical performance compared to conventional flotation cells (Finch and Dobby, 1990). Flotation columns have also found other applications outside mineral processing, such as de-inking of recycled paper (Finch and Hardie, 1999).In column flotation, a rising swarm of air bubbles generated by means of air spargers is employed to collect the valuable mineral particles and separate them from the gangue minerals in a countercurrent process. Wash water, which is continuously fed at the top of the column, is used to eliminate entrained particles and stabilize the froth. The column volume can be divided into two sections – the collection zone in which the bubbles collect the floatable mineral particles, and the cleaning zone (or froth zone) where product upgrading is enhanced through the removal of unwanted particles entrained in the water mixed with bubbles from the collection zone."

Jan 1, 2019

By G. Akdogan, C. Dorfling, L. Bumhira

Hydrometallurgical recovery of rare earth elements from phosphors obtained from end-of-life fluorescent lamps involve leaching followed by extraction of dissolved metals from solution. The similar properties of rare earth elements make selective recovery of individual elements challenging. The aim of the paper was to investigate the extraction of primarily Ce and Tb from a hydrochloric acid leach solution also containing other REEs (Y and Eu) as well as impurity metals (Ca and Al). Batch solvent extraction tests were performed with di-(2-ethylhexyl) phosphoric acid (DEHPA) at various extractant concentrations, pH conditions and organic-to-aqueous ratios. Extraction isotherms and separation factors for the range of conditions tested in the project are presented in the paper. Separation of individual REEs was not possible at the conditions investigated, but a mixed REE product can be obtained in a single extraction stage operated at a pH of 0.5 using 1 M extractant at an O/A ratio of 1; these conditions resulted in complete Y and Eu, 95 % Ce and 99 % Tb extraction with Ca and Al coextraction below 7 %. A better quality mixed REE product can be produced by operating at lower pH and extractant concentration conditions, but this would require a larger number of extraction stages. Keywords: Cerium, fluorescent lamps, phosphors, rare earth elements, solvent extraction, terbium

Jan 1, 2020

By I. Y. Motovilov, S. A. Telkov, M. B. Barmenshinova, Z. S. Abisheva

As a rule, removing waste rock at the beginning of the process improves the conditions for the flotation of minerals, and also reduces the operating costs of crushing, grinding and flotation. Lead-zinc ores in the Shalkiya deposit of the Republic of Kazakhstan are classified as difficult for enrichment, but taking into account their significant reserves, they are one of the main potential ores for the production of lead and zinc. Under industrial conditions, for processing ores of the Shalkiya deposit, a selective flotation concentration scheme with the production of lead and zinc concentrates is recommended. However, according to the recommended scheme, it is not possible to obtain high-quality concentrates of lead and zinc, with high metal recovery. The paper presents the results of studies of preliminary gravitational enrichment of lead-zinc ore in the Shalkiya deposit. Based on the fractional analysis for the particle size range of -40+8 mm, Henri- Reingard enrichment curves were constructed, separation density (2730 kg /m3) was established to isolate the light fraction with the lowest possible lead and zinc content, and gravity enrichment indices (0.7-0.71) were determined. It was found out that when using gravitational enrichment of the original large-sized ore in the Shalkiya deposit, tails are released at the beginning of the technological process with a yield of up to 20% with the extraction of silica up to 30%, calcite up to 11% and carbonaceous matter up to 12%. As a result of the studies, it was found that the beneficiation of crushed lump ore could be conducted in heavy suspensions. One of the positive aspects of the use of gravitational enrichment is the reduction of energy costs for grinding the heavy fraction, which was experimentally proved in this article. Keywords: Lead-zinc ores, Shalkiya deposit, fractional analysis, gravitational enrichment, floating fraction, sunken fraction, separation density

Jan 1, 2020

By S. C. Maji, U. S. Chattopadyay, S. Chaudhuri, T. Gouri Charan

India is devoid of coking coal resources, which is a vital raw material for coke making. As a result it has to import about 80% coking coal for steel making. Indian National Steel Policy, 2017 has envisaged crude steel production around 300 million tonnes Mt) by the year 2030-31 and forecasts coking coal demand of 161 Mtpa (million tonnes per annum). One way of improving the indigenous resources is to study the performance of existing washeries and to improve the clean coal recovery by adopting suitable technology. A study was carried out, mainly to focus on the fine coal circuit of an existing coking coal washery, to generate the desired quality of clean coal. The operating coking coal washery being studied is presently feeding coal from nine different sources with varying feed ash contents. The coal fines varied in ash content ranging from 19.5% - 43.6%. For treating the coal fines, the washery has operated a spiral for the past few years. Performance studies were carried out, and it was found that the spiral is not efficient in recovering the cleans. Detailed investigations were carried out on the coal fines using flotation techniques, and it was found that to achieve a low ash clean coal, the froth flotation is a much better beneficiation option than a spiral, The studies showed that the recovery of clean coal through a spiral is 67.7% at ash content of 27.0%, while flotation tests showed a recovery of 68.1% at 16.5% ash content. Keywords: Coking coal washery, near gravity material, spiral washer, froth flotation

Jan 1, 2020

By J. A. Theron

"This paper provides guidelines on performing mass and energy balance modelling involving coal and coal derivatives. Usually, the inputs to a pyrometallurgical process would be specified in terms of elements and compounds. Reliable thermochemical data is more widely available for species involving uniquely defined, relatively smaller molecules. However, in the case of coal, the molecules are extremely large and not uniquely defined. Consequently, modelling processes involving coal and its derivatives have several potential pitfalls. These are outlined in the present paper. The main findings are as follows:•,Coal proximate analysis should not be regarded as absolute; it could vary with several parameters, including heating rate •,For modelling, the use of ultimate analyses should be considered a preferable option to proximate analyses, where, ‘fixed carbon’ and ‘volatiles’ are not defined in terms of chemical composition •,Significant errors could be made if the larger molecules are neglected during calculation of the calorific value (CV) of coal gas (the gas liberated when coal is heated in the absence of oxygen) •,For elemental analyses, the oxygen content (which is calculated by balance) should be checked to ensure it is within the expected range •,For representation of sulphur in coal, double-counting due to SO3 in the ash analysis should be avoided •,Potentially, oxygen in coal could be represented as O2, H2O, CO, or CO2; however, use of some of these species without considering the experimentally determined gross CV results in significant errors in the energy balance. If coal enthalpy is calculated from elemental analyses without correction, representation of coal oxygen as H2O(l) gives reasonable accuracy (up to 6.9% error for a set of South African coals). Coal volatiles could be represented by a complex mixture of compounds, even using different oxygen-containing species than the abovementioned four, provided the enthalpy is corrected •,It is recommended that an ‘enthalpy correction value’ is incorporated in energy balances involving combustion, devolatilization, or conversion of coal and coal derivatives, e.g. coke, char, or tar. That would imply that proximate analysis, elemental analysis, as well as the gross CV would be required for all solid or liquid coal-derived substances being modelled. No other correction due to the form of carbon being different than graphite should be used, as that would imply doublecounting some effects."

Jan 1, 2014

By R. W. O. Kersten

"The aim was to investigate stable pillar design procedures currently in use on chrome and platinum mines, subject them to a critical appraisal, and review some of the alternatives proposed in the literature. It is concluded that these pillar design methodologies suffer from drawbacks that can be detrimental to the mining industry and that an improvement is essential.It was decided that the increased availability of analytical models and failure criteria could be exploited to update the empirical to an analytical/empirical method. The proposed methodology adopted FLAC2D simulations using the Hoek-Brown failure criterion to calculate full stress-deformation curves for typical pillars. The mine stiffness concept was introduced to determine the pillar load, integrating the influence of the pillar and strata stiffness and mining geometry.The proposed methodology was calibrated by comparing experimental and predicted deformations and failure depth of pillars in one hard-rock mine. The conclusion was that the methodology is an improvement over the one currently in use. An improved factor of safety is given by the intersection point of curves for the stiffness of the system and the pillar strength. IntroductionThe paper describes part of the findings of a PhD thesis (Kersten, 2016) on an alternative stable pillar design methodology. This publication contains sufficient information to implement the method for further interpretation and testing. A detailed comparison of measured and predicted data is dealt with in the thesis, and the results show that there is promise in the methodology.The design of stable pillars in mining is of fundamental importance, not only in the Bushveld Complex mines but throughout the mining industry. Wherever mining occurs, pillars are formed at some stage and it is essential to predict their behaviour, i.e. whether they will burst, yield, or remain stable. Although of crucial importance, pillar design still suffers from major shortcomings and weaknesses."

May 1, 2019

By U. Atici, R. Comakli

"Understanding rock properties is of vital importance in mining and civil engineering. Rock properties are generally determined through laboratory tests, which require a large number of uniformly shaped samples and technicians who use indirect methods to determine rock properties. The texture coefficient (TC), which represents various textural parameters including the size, shape, orientation, and relative proportions of grains and the matrix, is an efficient means for indirectly evaluating the mechanical behaviour of rocks. In this study, 12 different rocks (ten granites, one diorite, and one gabbro) were collected and their engineering properties. including specific gravity, porosity, P- and S-wave velocities, Schmidt rebound hardness, uniaxial compressive strength, and elasticity modulus, were investigated. Relationships between TC and these measured properties were assessed. Statistically significant (at the 95% confidence level) correlations were found between the TC and compressive strength, elasticity modulus, sound velocity (both compression and shear waves for dry and saturated samples), but no significant correlation was found with Schmidt rebound hardness. The results suggest that it is possible to use the relationships between TC and the mechanical properties of rocks in rock classification. IntroductionIn every aspect of rock engineering, knowing rock properties is essential for determining and predicting the behaviour of rocks and rock masses. Mechanical properties of rocks are important parameters in mining, geology, and engineering, and they are affected by textural characteristics such as grain size, grain shape, cement type, degree of cementation, mineral content, and the ratio of the grain area to the matrix area. The effects of textural properties on the physical and mechanical properties of rocks have been investigated for a long time, and previous studies revealed close relationships between single parameters related to rock texture and mechanical properties of rocks (Akesson et al. 2001; Gunsallus and Kulhawy 1984; Merriam, Rieke, and Kim, 1970: Onodera and Kumara 1980; Prikryl 2006; Ulusay, Tureli, and Bider, 1994; Tugrul and Zarif 1998). Rock properties are usually determined using mechanical tests, which require uniformly shaped and relatively large samples. However, it is not always possible to obtain such samples. Furthermore, rocks with similar textural and structural properties may exhibit different mechanical behaviour. Howarth and Rowlands (1987) developed a texture coefficient (TC) that expresses all the main textural parameters of a rock in a single dimensionless quantity based primarily on statistical assessment. The TC makes it possible to understand the effects of a rock’s textural properties on its mechanical properties using thin sections instead of large, uniform samples. Ersoy and Waller (1995) investigated the relationship between strength and TC for limestone, sandstone, and siltstone; they also examined the effect of rock properties and TC on drilling performance. Ozturk, Nasuf, and Bilgin (2004) explored the relationship between rock quality and TC by evaluating the uniaxial compressive strength (UCS) values of limestone, basaltic andesite, and tuff. Gupta and Sharma (2012) evaluated TCs as well as the petrophysical and mechanical properties of quartzite from the lower and higher Himalayas and found a strong relationship between TC and UCS. Similarly, Ozcelik, Bayram, and Yasitli (2013) investigated the relationship between TC and various parameters in marble by applying the statistical analyses developed by Ozturk, Nasuf, and Kahraman (2014)."

Jan 1, 2019