Search Documents

By P. T. Slatter

The role of viscometers for rheological characterization in the hydraulic transport of solids is discussed, and rheological models and phenomena pertinent to the hydraulic transport of backfill material are reviewed. The two main branches of viscometry (rotational and tube) for the measurement of mining slimes rheology are compared. The effects of solids concentration, particle size and chemical additives and the practical application of rheology in backfill system design are presented.

Jan 1, 1988

By J. Dvorácek

Mine water poses a significant problem in lignite coal mining. The drainage of mine water is the fundamental prerequisite of mining operations. Under the legislation of the Czech Republic, mine water that discharges into surface watercourse is subject to the permission of the state administration body in the water management sector. The permission also stipulates the limits for mine water pollution. Therefore, mine water has to be purified prior to discharge. Although all mining companies operate mine water treatment plants, there is hardly any utilization of the final effluent. The high expenditure involved in the more sophisticated methods of mine water treatment has fuelled the search for the commercial utilization of the effluent.

Jan 1, 2012

In Chapter 1 we mentioned that up to three components are considered to represent the variability of regionalized variables : a deterministic component, a correlated component, and a random component. In Chapter 2 we considered some of the models based on the assumption of randomly distributed values: the random component was assumed to be the only one present. In Chapters 3 to 11, more general but also more complex models were described, based on the assumption that the values were correlated with the possible presence of random noise. In this chapter we introduce a model used to analyse non- stationary data, for which a deterministic component (or drift) must be assumed. This model will be used for estimation purposes, and the corresponding optimal linear estimation procedure is known as universal kriging. The theory of universal kriging is much more complex than the theory of kriging. Fortunately, in the majority of mine valuation problems the drift can be ignored and the simpler form of kriging can be used.

Jan 1, 1978

By A. Cáceres, L. Navarro, N. Altmann, R. Betinol

Copper concentrate is produced in Chile in the high Andes and transported to the coast for shipping through long distance pipelines. The pipeline and fittings installed have to withstand the high abrasion and corrosion environment of the copper concentrate. To protect the inside tube wall against corrosion, high density polyethylene (HDPE) liner is used. While HDPE liners can protect the steel wall against corrosion, it has also been assumed that the low surface roughness of the HDPE liner also inhibited incrustation caused by calcium carbonate and sulphate deposition. Recent events have indicated otherwise in a couple of concentrate pipelines for which BRASS has provided consultant service. This paper delivers a case study which includes (i) formation of incrustation and adhesion to HDPE liner, (ii) the detection of the incrustation through hydraulic analysis, (iii) the mechanism of removal of incrustation, and (iv) the catastrophic effects of their presence in an emergency situation.

Jan 1, 2007

By P. A. Corbeil, C. Swart, F. Miller, V. Falmagne

Automation technologies are increasingly adopted in different mines worldwide. This paper presents the technical aspects of vehicle automation applications at two Noranda operating mines. Teleoperation of 5 LHD’s at Noranda’s Bell Allard mine is used in order to meet aggressive mining requirements under difficult conditions. The equipment is performing under high demand for reliability and performance. In order to further increase performance, automated loading has been commissioned on one of the LHD’s and is being tested. Results to date indicate improved performance with respect to an operator through better filling of the bucket, fewer passes at the muck pile and shorter cycles. Noranda’s Brunswick mine is currently commissioning an automated haulage zone where two optically guided trucks transport ore from a remotely operated chute to two ore passes. Expected operational advantages include increased system capacity and enhanced safety and comfort for the operator. Monitoring to date indicates that the trucks can operate reliably under the conditions imposed by mine production.

Sep 1, 2001

By J. A. L. Napier

Computer programs which rely on the assumption of linear, elastic rockmass behaviour have been used routinely for many years in the South African mining industry to determine the stresses arising in the vicinity of mining excavations. A description is given of the theoretical background and novel solution techniques used in the recently developed MINSIM-D boundary element program. Features of the program include the ability to analyse multiple interacting reef planes at arbitrary orientation to one another, slippage on fault planes and the effects of backfill placement in mined areas. Incremental mining sequences can be modelled, and digitizer and front-end routines have been written to assist in the capture of mining outlines and general problem specifications. A brief indication is given of design problems for which MINSIM-D is appropriate, as well as some future extensions of the system.

Jan 1, 1987

Orange Free State Branch Minutes of a General Meeting held at the Welkom Mine Canteen on Wednesday, 21st November, 1973, at 4-.00 p.m. Present Mr E. T. Wilson (Chairman) in the Chair, 15 Fellows and Members, 12 Associates, Graduates, and Students, and 11 Visitors (Total 38). Book review Symposium on Environmental Engineering in Coal Mining The 20 papers presented at the symposium cover some 14 subjects related to environmental engineering in coal mines, such as dust, heat, gas emissions, ventilation, noise, and lighting. Competition for students Each year the Institute offers a prize (or prizes should the entries warrant it) of up to Rl00 for the best paper or dissertation on a topic appropriate to the interests of the Institute. A.8. & T.8. Trust The A.S. & T.S. Trust has been formed with the specific object of establishing a worthwhile source of funds which can be used to further the progress of science and technology in South Africa. Nim reports Report No. 1573 A pilot-plant investigation of a flotation model. Report No. 1597 The role of sulphur in the flotation of sulphide minerals. Report No. 1599 The analysis of silicate rocks and metallurgical slags by X -ray-fluorescence spectrometry. Report No. 1607 The analysis of ion-exchange resins by X-ray- fluorescence spectrometry. Report No. 1617 An electronic differential precision thermometer.

Jan 6, 1974

By K. Grossauer, W. Schubert

For tunnelling in poor and heterogeneous ground a reliable on-site short-term prediction is required to allow a safe and economical construction. With the use of absolute displacement monitoring data and advanced software for evaluation of the data in combination with a routine update of the geological model a considerable increase in the reliability of the short-term prediction could be achieved. The paper deals with theoretical considerations, numerical studies and case histories from alpine tunnels in heterogeneous rock. Tunnelbau in schlechtem und heterogenem Gebirge erfordert auf der Baustelle eine verlaessliche Kurzzeitprognose, um ein sicheres und wirtschaftliches Bauwerk gewaehrleisten zu können. Die Verwendung von absoluten Verschiebungsmessdaten und deren Auswertung mit Hilfe von hoch entwickelten Softwareprogrammen in Kombination mit einer laufenden Aktualisierung des Gebirgsmodells ermöglicht eine erhebliche Erhöhung der Verlässlichkeit von Kurzzeitprognosen. Dieser Beitrag beschäftigt sich mit theoretischen Ansätzen, numerischen Modellen und Fallstudien von Tunnelvortrieben in heterogenem Gebirge. La construction de tunnels dans la montagne hétérogène de faible résistance exige un pronostic fiable à court terme sur le chantier, afin de pouvoir garantir un ouvrage sûr et économique. L'utilisation des données de mesure de décalage absolues et de leur évaluation à l'aide de logiciels sophistiqués en combinaison avec une mise à jour courante du modèle de la montagne permet une augmentation considérable de la fiabilité des pronostics à court terme. Cette contribution présente des approches théoriques, des modèles numériques et de cas pratiques de progressions de tunnels dans des montagnes hétérogènes.

Jan 1, 2003

By J. Rowe

The attrition of activated carbon, and the loss of gold associated with it, is of significant economic importance to the operation of a CIP/CIL circuit The focus of this study was on activated carbon recovered from an elution circuit which was deemed too fine for reuse. Results of cyclosizer and laser particle size analysis identified that most of the carbon was contained in the larger size fraction suggesting formation by abrasion. Digestion and AAS analysis of the individual size fractions identified a disproportionate concentration of gold in the finer size fractions which was identified by scanning electron microscopy to be due to the presence of fine metallic gold formed as a result of the acid washing process. Attempts to strip the remaining gold using sodium hydroxide or sodium sulphide based solutions proved unsuccessful due to poor elution efficiencies and re adsorption of gold. Upgrading of the material by froth flotation was also investigated using various conventional flotation reagents which had some success in concentrating the free metallic gold, however, grade and/or recoveries were considered less than desirable. Greater success was ultimately found in the transfer of gold from the fine carbon material onto coarser virgin activated carbon using a caustic cyanide solution. Sodium sulphide, sodium chloride and ammonium chloride based solutions were also trialled but proved less successful owing to their inability to mobilise metallic gold or gold cyanide species. Repetition of the transfer process proved capable of stripping 97% of the gold contained on the carbon fines for abed volume ratio of 6:1.

Jan 1, 2009

A method, based on the macroscopic population balance equation allowing the dynamic simulation of a multistage adsorption plant given an arbitrary rate expression, is described. The partial differential equations resulting from this method are simplified using the method of characteristics. The simplified equations can then be solved using standard ODE integration techniques. The procedure is illustrated by means of data typical of a CIP application. Requirements for data structures which allow the incorporation of the carbon-in-pulp (CIP) model (and other unit models) into a general metallurgical simulator are described. Based on these considerations, it is concluded that more generality in the description of distributed quantities such as loading and particle size is provided by the use of mass density at fixed class boundaries, rather than discrete masses within fixed classes as has hitherto been used by ore-dressing simulators.

Jan 1, 1987

By D. P. Power

The development of an accounting code for the metallurgical industry has been debated for many years. Following the series of tumultuous economic events in the nineties and earlier this century that have affected business, there was much need for legislation to be passed to minimize potential risk. One of the consequences to affect the mining industry was the enhancement of existing geological codes for reporting. In parallel the sponsors of AMIRA have taken the brave step to develop a generic code for metal accounting that embraces the full gamut of industry. In this paper the code will be reviewed in the context of metallurgical accounting as it applies to coal preparation. This paper provides the background for this development and an insight into the mechanism for the debate for the generation of this code. The key elements of the code as it applies to the coal industry are explained but most importantly the key opportunities are highlighted.

Jan 1, 2010

In the previous chapters we did not concern ourselves with the fact that the geostatistical properties of the regionalized variable x (z) are affected by the size of the support w. Ideally, when analysing an ore body O, we would like to determine the properties of point values in O, i.e. the properties of a regionalized variable with an infinitely small support, say w°. In practice, point values are never known. To measure the value of x at a point z we must take a sample of size w centred at z. For example, if we are interested in the grade at point z, w will be the volume of the sample 'centred in z; if we are interested in the metal accumulation at point z, w will be the area of the cross-section of the sample centred in z. We say that the measured property is a regularized variable, and so the semivariogram (or covariogram) calculated from the values obtained is also regularized. In the majority of the mining problems which can be solved by geostatistical methods, we are interested in the properties of blocks of ore and of samples with a finite non-point support. It will be shown in the following chapters that these problems can be solved by calculating the average value of the point sernivariogram in volumes or areas of finite dimension, or between such volumes or areas. These calculations would be simple if we knew the point semivariogram. This means that the point semivariogram has to be estimated from the observed regularized semivariogram(s).

Jan 1, 1978

By K. -C. Han, S. -Y. Kim

This study was carried out in a bid to survey Ganboke gold mine in Pala, Mayo-Kebbi Province, southwestern Chad. The Ganboke mine is now being prepared for investigation and development, and the research sought to perform rock properties tests of relevant rocks through its geological survey, rock mass evaluation, and in-situ and laboratory tests and analyze the stability by the means of limit equilibrium method, thereby suggest appropriate methods to design the pit slope and individual benches in Ganboke mine. Targenting GS1-GS1' section with the depth of 180 m and GS2-GS2' section with the depth of 150 m in the survey area, we set the height of the bench at 10 m and 15 m, and the angle of bench at 60o and 70o, respectively, and then reviewed the stability for dried condition and saturated condition. The result of examination indicated that it would be appropriate to set the bench height at 10m and the bench angle at 60o in developing the mine slopes concerned. Cette étude est réalisé dans le but de faire des recherches sur Ganboke mine d'or dans Pala, Mayo-Kebbi Provence, Tchad sud-ouest. La Ganboke mine est actuellement en train de étre disposé à la investigation et le développement, et la recherche a cherché à accomplir les tests de la propriétè des roches par l'intermédiaire de la étude géologique, la évaluation de la masse rocheuse, et in-situ et les tests de laboratoire et à faire l'analyse de la stabilité au moyen de la méthode d'équilibre limite, de ce fait à proposer le methodes propres à faire le plan de la pente et les gradins individuels dans la Ganboke mine. Concernant la section de GS1-GS1' avec la profondeur de 180 m et la section de GS2-GS2' avec la profondeur de 150 m dans l'aire d'étude, nous mîmes l'hauteur du gradin à 10 m et 15 m, et l'angle du gradin à 60o and 70o, respectivement, et alors nous avon reconsidéré la stabilité pour la condition sèche et la condition saturé. Le résultat d'examen a indiqué que il pourrait être propre de mettre l'hauteur du gradin à 10m et l'angle du gradin à 60o en développant les pentes de la mine concerné. Diese Forschung wurde um das Gold Bergwerk zu überschauen in Pala, Mayo-Kebbi Provinz, der Südwesten Chad ausgeführen. Das Ganboke Bergwerk ist jetzt zur Erforschung und Entwicklung vorgebereiten worden, und die Forschung suchte zu ausführen die Prüfungen der Felseigenschaft durch geologischen Überblick, Felsmassewertbestimmung, und In-situ und Lavorversuch auf, und analysieren die Stabilität mittelst der Grenzegleichgewichtsmethode, da-durch schlug die angemes-sen Methoden zu entwerfen die Böschungen und die einzelnen Bänke im Ganboke Bergwerk vor. Für den Abschnitt GS1-GS1' mit der Tiefe von 180 m und GS2-GS2' mit der Tiefe von 150 m in der Überblickfläche, wir setzten die Bankhöhe an 10 m und 15 m, und das Bankwinkel an 60o und 70o, beziehungsweise, und dann zurückblickten die Stabilität für die getrockneten Bedingungen und die gesättigten Bedingungen. Das Ergebnis der Untersuchung anzeigte daß es würde angemessen sein zu setzen die Bqnkhöhe an 10m und das Bankwinkel an 60o in der Entwickelungen vom Bergwerk be-teiligt.

Jan 1, 2003

By J. E. Quintanilha de Menezes, L. Ribeiro e Sousa, I. Ferreira

The Paradela dam has a spillway that makes it possible to discharge a maximum flow of 400 m3/s, the water discharges being on the river Sela, about 2km upstream the river Cávado. After the first use of the spillway, some landslides have occurred on the right bank of the river Sela as a result of the erosion generated on its base about one hundred meters near some houses of a small village. This paper presents a study performed for the slope stability where finite element results are taken into account. A monitoring program was established that essentially consisted of measuring sub-horizontal displacements in points located in boreholes at different depths and also phreatic levels. Some monitoring results are presented. Le barrage de Paradela a une décharge d’eau avec un débit maximal de 400 m3/s sur le fleuve Sela dans un endroit situé environ 2km sur l'amont du fleuve Cávado. Après la première utilisation de cette décharge, quelques glissements se sont vérifiés sur le versant droit du fleuve Sela résultant de l’érosion générée sur sa base dans une zone à 100 mètres de quelques maisons d’un petit village. Ce travail présente une étude réalisée concernant la stabilité du versant avec la considération des résultats des éléments finis. Un programme d’auscultation a été mis en place pour mesurer les déplacements sub-horizontales sur des points situés en forages à des différentes profondeurs. Quelques résultats de cette monitorization sont présentés. Der Ueberlauf des Padela Dammes ermoeglicht das Ablassen einer max. Wassermenge on 400 m3/s in den Sela Fluss, ca. 2 km flussaufwaerts des Cávado Flusses. Nach Inbetriebnahme des Ueberlaufes traten in ca. hundert Meter Entfernung zur Wohnbebauung eine Dorfes Erdrutschungen am rechten Ufer des Sela Flusses auf, die auf Erosionen am Flussbett zurückzuführen waren. In der vorliegenden Studie wurde die Boeschungsstabilitaet unter Beruecksichtigung von Berechnungsergebnissen der Finiten Element Methode untersucht. Es wurde ein essprogramm ausgearbeiten, das im Wesentlichen auf Messungen der Untergrundverformungen in verschiedenen Bohrlochteufen sowie auf Grundwassermessungen beruht. Im Folgenden sind einige Messergebnisse dargestellt.

Jan 1, 2003

The relationship between the aspect ratio and physical fields in aluminium reduction cells wase studied numerically. The aspect ratio was firstly defined and 7 kinds of 320 kA cells with different aspect ratios were put forward. By using numerical simulation, the relationship between the aspect ratio and the electric-magnetic flow fields, magnetohydrodynamic stabilities, and thermal stability was discussed. It is concluded that the electric-magnetic flow field distributions are greatly affected by the aspect ratio. From the perspective of magnetohydrodynamics stability, the larger the aspect ratio is, the more stable the cell will be. A larger aspect ratio is more beneficial for the thermal stability under the same current density.

Jan 1, 2011

By A. W. S. Schumann

G. H. HENDERSON INTRODUCTION It is a pleasure to have read Mr Schuman's wide ranging and highly informative paper and also to congratulate him. Before proceeding to give views on the use and preservation of coal in South Africa, it is perhaps apposite to look at the South African scene on the world canvas. Mr Schuman's figures in his Table I are highly revealing. The Table shows expected consumption, export and import demand in the year 1980 by the continents of the world. The Republic of South Africa is shown as one of only four exporting countries at that date, although in the interim, Australia has become an exporter and will continue to grow as an exporter of coal. Huge fields have been proved in recent years. There is reason to think that for some years, coal and its economical uses have been under-estimated. The importation of oil carries a constant question mark and the exciting possibilities expressed in the uses of Dranium a decade ago for power generation, appear to have met with some setback. Present indications are that breederreactors will be delayed until about 1990. The capital cost of such reactors, too, is enormous. By 1980 it has been estimated that 20 per cent of the energy required by the leading ten nations of the world will come from nuclear power, but 80 per cent would then remain for fossil fuels and hydro generation. COAL DEMAND AND USAGE It is of interest to record that coal production in all the countries of Europe has declined markedly in the last 25 years. This is due to competition from oil, gas and nuclear power and also for sociological reasons; but there are signs that efforts will be made to resuscitate, at least in some measure, the mining of coking coals for which there is a strong and growing demand. Possibly all hard coal now mined in Germany has coking properties and that country is self-sufficient in this commodity. Germany might perhaps consider selling coking coal to her neighbours and importing additional coal for power-raising purposes. It has been frequently stated in American publications that the coal reserves of the D.S.A. are sufficient for 1 000 years. What that means, I do not know, but it is at least an indication that America has huge reserves of coal. Domestic and export demand in the U.S.A. in 1968 exceeded production for the first time. In 1969, the United States' output amounted to some 556 million tons, an increase of less than 2 per cent on the previous year. The output achieved fell short of demand, domestic and export, by more than 7 000 000 tons. This has caused and is continuing to cause anxiety in the U.S.A. In many countries young men for a number of reasons, largely environmental and social, are not coming forward in sufficient numbers for mining requirements and the same pattern has revealed itself in the Republic. This is, indeed, alarming because mining men in large numbers will be required to meet the estimated tonnage requirements set out in the following Table: It is difficult to visualise the replacement of these coal requirements either by oil or by nuclear fuels. The coal must be mined and men are required for the purpose. Unless a greater section of the population is allowed to participate in the more responsible aspects of mining, there is a serious likelihood that demand will not be matched by output. In addition to local demand, a lively export trade is possible. Exports are dependent on the future capacity of the South African Railways and Harbours and on Government policy. An export figure of 10000000 tons per year during the decade can be visualised and this figure would raise total estimated sales production to a figure of over 100 000 000 tons in 1980. At this point, it is necessary to examine the desirability or otherwise of exporting coal. How much coal

Jan 11, 1970

PRESENT Minutes of the Annual General Meeting, held on Friday the 20th August, 1971, at 7.45 p.m. in the Boulevard Hotel, Witbank. Mr J. Meintjies in the Chair, 15 Fellows, 7 Members, 2 Associate Members and 10 Visitors. Mr Meintjies extended a welcome to all the members and their wives, and the many visitors and their wives. A special welcome to Mr Robinson the outgoing President of the Institute and Prof. Howat the President Elect and Mrs Howat. The minutes were read and adopted. The Chairman read his Report for the past year and also the Financial Statement. Mr E. Jackson proposed the acceptance of the Financial Statement and Mr Winship seconded it. Mr Meintjies thanked Mr Vic Robinson for all the help he had given throughout the past year. He also thanked the outgoing Committee for their support, and the Secretary for arranging the function. Mr Meintjies congratulated Mr Alex Bain on his election as Chairman and asked him to continue with the meeting. Mr Bain thanked Mr Meintjies for the good work which he had done in the past year and then announced the names of the Committee for the ensuing year. A. Bain Chairman G. Carlyon Vice Chairman J. Meintjies Past Chairman I. Botha E. J. Jackson A. W. L. Brereton B. D. McNiece N. J. P. Coetzee A. N. Turner B. T. Hosking W. D. Winship Mr B. D. McNiece proposed a vote of thanks to the retiring Chairman. The meeting closed at 8.10 p.m. Dinner was served and after coffee Mr Vic Robinson gave a most amusing and informative review of the Institute's activities. Mr Bain thanked Mr Robinson for his address and Prof and Mrs Howat and Mr Robinson for coming down to Wit bank for our Annual General Meeting.

By Lene Hansen

Sulphuric acid plants now and then see unexpected increases in SO2 emissions. Our experience is that it is not always easy for the acid plants to pin point the exact reason for the SO2 increase. Sometimes a catalyst bed is sacrificed even though reliable emperature readings and gas analyses could have verified satisfactory activity, hence the leaking heat exchanger is never discovered and after the shutdown the SO2 emissions are still too high. This paper will look into some basic catalyst management tools: Good record keeping is a necessity for collection of representative catalyst samples, reliable operating data is an inevitable starting point for all analysis and trending of operating data is a must for plant optimisation and troubleshooting. In addition these catalyst management tools will give a continuous updated knowledge of the catalyst status in each bed providing the basis for a more effective replacement strategy. The most typical reasons for increasing SO2 emissions will be presented together with easy guidelines for cross-checking of operating data, temperature optimisation and good record keeping.

Jan 1, 2009

By Q. G. Reynolds

The problem of direct-current plasma arc behaviour, interaction, and dynamics is considered in the context of metallurgical DC arc furnace applications. Particular attention is paid to the transient flow behaviour of arc systems. A mathematical formulation of the physics used to describe the arc system is presented, and includes the spatial and temporal evolution of fluid flow, heat transfer, and electromagnetism. Based on this formulation, a numerical model is developed using a finite difference approach on a regular cartesian grid in both two and three dimensions, with a special focus on robust stability, high resolution modelling, and high performance. A collection of results produced using the numerical model to study pilot plant-scale furnaces is then presented. These address a range of process and design variables and their effect on the numerical model?s results. Where possible, the qualitative behaviour of the model is compared to available experimental data. A number of novel effects and phenomena are seen in the dynamic behaviour of the DC plasma arc model for both single and multiple arc systems, which may lead to improved understanding, control, and manipulation of such systems where they occur in industrial applications. Keywords Pyrometallurgy, direct-current furnace, plasma arc, modelling.

Jan 1, 2010

By J. Bai, C. Zang, W. Yuan, L. Zhuang

"This study evaluates the efficiency of using mechanochemical reduction to assist the extraction of zinc from zinc ferrites by alkaline leaching. The transition of zinc ferrite into a metastable state after mechanochemical reduction with metallic iron contributes to the ready dissolution of zinc from the activated zinc ferrites in alkaline solution. Zinc ferrites were mechanochemically reduced under conditions of Fe:ZnFe2O4 mole ratio of 2:1, using 5 mm diameter stainless steel balls as the activation medium at a mass ratio of balls to raw materials of 25:1. Subsequent leaching in 6 mol/L NaOH solution at 90°C resulted in more than 70% Zn extraction. These results may be used to further develop a hydrometallurgical process for recovering zinc from zinc ferrites in alkaline solution.IntroductionZinc is an important nonferrous metal required for various applications in batteries, solder, dielectric materials, and piezoelectric materials. High-grade, easily concentrated ores are becoming scarcer: to meet world demands for zinc, it is necessary to develop economical and environmentally safe metallurgical technologies that use secondary materials such as leach residues of roasted zinc concentrates and electric arc furnace dusts. The zinc in the majority of these wastes occurs in the form of zinc ferrite (franklinite, ZnFe2O4), which is very stable and insoluble in most acidic and alkaline solutions (Langova and Matysek, 2009; Youcai and Stanforth, 2000).Several pyrometallurgical and hydrometallurgical processes (Jankovi et al., 2014; Langova and Matysek, 2010; Morcali et al., 2012; Shawabkeh, 2010) have been proposed for the treatment of zinc ferrites. Pyrometallurgical methods are costly owing to high energy consumption and their emissions of toxic gases are considered an environmental problem; hydrometallurgical methods are preferred for process economy and environmental reasons (Copur et al., 2004)."

Jan 1, 2016