Search Documents

The overall mineral assemblage of the Broken Hill orebody is similar to that of a skarn. However unlike the classical skarns it does not represent contact metamorphism of pre-existing limestone beds. It is a stratiform accumulation of exhalative-sedimentary silica, sulphides, (Ca, Fe, Mn) silicates and carbonates, fluoride and phosphate that has later undergone regional metamorphism. The resulting development of coarse-grained quartz, sphalerite, galena, rhodonite, bustamite, mixed carbonates, fluorite, apatite and many minor phases has produced a sulphide-rich skarn- like mass that like many others of its kind, might best be referred to as a stratiform skarn. Recent investigations of occurrences of this type are showing that the incidence of concentrations of particular metamorphic minerals in and adjacent to such orebodies are closely related to bedding and to former sedimentary facies patterns. From this it appears that the mineral assemblages may be used to deduce not only the nature of the materials laid down in the beds, but also the physico-chemical conditions under which sedimentation took place. Several stratiform deposits, including Broken Hill, are used as examples to show that variations in these conditions may have a profound influence on the subsequent response of the rock to regional metamorphism and hence on the mineral assemblages that finally develop.

Jan 1, 1983

In current mining practices it is desirable to bring about a substantial increase in both rate of face advance and productivity, thus compelling the mine engineer to consider better and more productive ways of mining. Transportation plays a vital role in achieving improved productivity. This paper focuses on horizontal transportation systems. Horizontal transportation is defined from the shaft station area to the end of the stope cross-cut. Traditionally, horizontal transportation is inefficient in comparison to vertical transportation efficiencies and accounts for a high proportion of South African mine accidents. Thus, the role of horizontal transportation has been investigated in detail, as it offers the greatest chance for improvement. The purpose of the paper is to provide a best practice guideline for the design of horizontal underground transport systems. It is intended that the document be used by mine designers and planners as an aide memoir during the design process and benefit may well be gained from retrofit and reorganisation activity. The paper looks at transport systems suitable for deep level gold mining in South Africa and evaluates such items as management, scheduling, and haulage systems. Current transportation systems for personnel, material and rock transportation are evaluated, gaps in technology highlighted, and recommendations made.

Jan 1, 2003

Structure is the most important control on the location of hematite deposits in the Hamersley Province in the Pilbara region of Western Australia. Unravelling the structural history of mineralised areas, however, is complex because the Province underwent three orogenic events, and three regional extension events followed the main episodes of compression. Steeply E to ENE dipping normal faults immediately predated deposition of the Beasley Quartzite of the lower Wyloo Group after the Ophthalmian Orogeny. Steeply NE and SW-S dipping normal faults predated deposition of the McGrath formation of the middle Wyloo Group, following the Panhandle Orogeny. The resulting interpreted rift was subsequently reactivated and tilted during the Capricorn Orogeny. Late NE trending normal faults are confined to the eastern half of the Province. High-grade hematite deposits occur in rocks that range from intensely folded (Mt Whaleback, Giles Mini) to rocks with less, or only little deformation (Paraburdoo, Channar), suggesting that folding had at best limited influence on ore genesis. In contrast, extensional normal faults are always associated with high-grade hematite deposits, and in particular faults that link the underlying dolomites of the Wittenoom Formation with the overlying iron formation, thus providing a pathway for mineralising alkaline fluids.

Jan 1, 2005

A good blasting result is largely achieved by putting the right amount of explosive in the right place and firing it in the right sequence. The best blast design in the world will not break the rock if the planned blast layout is not achieved in the field. One of the critical inputs into a blast design is the diameter of the blasthole. An assumption that the actual blasthole diameter is equal to the nominal drill bit size is generally made. However, a small deviation in the actual hole diameter from that assumed will have a dramatic effect on the hole's volume and, consequently, where the explosive is placed in the blasthole. Any discrepancy in actual blasthole volume from that assumed will be reflected in 'overloading' or 'underloading' of the blasthole with the calculated explosive quantity and, thus, reflect in explosive distribution and stemming lengths. When loading to a 'collar' in an oversized hole extra explosive will be used, resulting in higher than planned powder factors and explosive costs. A blasthole diameter assumed from a nominal drill bit size is an incorrect assumption. This paper will present some actual blasthole diameter measurements illustrating the variation in hole diameters seen in the field and proposes a practical method for achieving a useful blasthole diameter measurement.

Jan 1, 1990

By Churchill M

One of the significant changes which has occurred in the flotation sector of the minerals industry, over the last two decades, has been the successful development of the very large flotation machine. The last twenty years has seen the volume of the individual cell, available to the industry, increase by a factor of ten. Today, new mills, and existing mill modernization plans, can consider individual cell volumes as large as 85 m3 with the certainty that the metallurgical performance capability of the large machine is inherently as efficient as the smaller machine. The large flotation cell offers benefits to the user such as reduced floor space, lower total power, process control simplification, and reduced maintenance. However, the flotation machine size selection process, based on overall installation economics, must include such factors as process tonnage, flotation kinetics characteristics, and metallurgical performance objectives. This paper discusses the flotation machine size and circuit arrangement selection based on both metallurgical performance and economic considerations. A computer based flotation kinetic analysis technique, which has proven valuable in quantifying the metallurgical response of flotation circuits having different cell sizes and number, is utilized to illustrate the selection of cell size and arrangement best suited to a client's Requirement.

Jan 1, 1991

To obtain the best economic results in treating the comparatively high-grade ore of the Kalgurli Mine a good roast is essential, and the effort of this paper will be to briefly describe the plant and the conditions prevailing in roasting, which, on an output of about 360 tons daily, subsequently enable an average extraction of slightly over 94% of the gold to be maintained.The percentage of sulphur in the ore varies between 2'5 and 3'5, the lower figure existing when treating oxidized ore. The average gold value is 13 dwt; this occurs chiefly as sylvanite, containing .a lower percentage of silver than is usually present in that mineral. One ore body in the mine contains hessite (Ag2Te) in considerable quantity, associated with gold telluride, and when ore from this stope is being treated the silver contents rise to as high as 20 dwt, but the usual amount is about 9 dwt. per ton, of which only a small proportion is extracted. Other metallic minerals, in addition to pyrite, occurring more or less commonly, are arsenopyrite, magnetite, chalcopyrite, and the tellurides calaverite (AuTez)and coloradoite HgzTe3). A small proportion only of the value occurs as free gold.

Jan 1, 1911

By Dagbert M, Belisle J. M

Several packages for the estimation of porphyry type deposits exist and produce acceptable results. Recently, however, the need for really trying to obtain maximum information from drill holes has prompted the generation of geostatistical solutions to the problem of reserve calculation. It will be shown that an understanding of a few key points is necessary in order for these solutions (variogram analysis and kriging) to he really ptimal. The most important new fact to be ecognized in open pit ore reserve estimation is that block models do not give an accurate rediction of ore tonnage when a cut-off is pplied but often yield an overestimation of p to 25%. When checked against extensive roduction records, geostatistical methods provide a practical solution to this problem. t now seems that the model used here gives he best possible results. Inasmuch as all arge deposits present very similar geostatis- ical characteristics only a few days of work re necessary to produce a complete inventory f these. The techniques used will be demon- trated using examples from five deposits and pecial emphasis will be given to the accurate rediction of tonnages of ore and waste using he new "transfer function" technique of G.

Jan 1, 1977

Railway cuttings require the highest care in design and construction, as the pro- bability of failure concept does not allow any risks to be taken with hundreds of lives being dependent on continued stability. The usual geological influences on stability such as bedding dip and faulting are often discus- sed, but other major problems may come from the presence of sheet joints, cross joints differential weathering and seasonal or daily changes of ground water level. The design of rail cuttings has taken little note of the problems imposed by cons- truction methods, and low cost quarrying methods may seriously weaken a designed shape. The presence of berms and drainage ditches is often difficult to achieve as blasting imposes cumulative stresses and jointing on rocks. Contractors face difficulties in bidding extensive rock cuttings as variations in rock properties are often concealed and overshoot- ing the toe of a batter can easily occur. Modern blasting methods offer some solutions, but the best guide to safe successful cuttings is extensive investigation and the continued involvement of engineering geology at the stages of feasibility, design, construction and maintenance. Two case histories show the influence of folding, faulting and construction methods in the case of slope failure and slope stability.

Jan 1, 1979

By R L. Koenig

Previous papers at Metplant conferences have covered the options for crushing and grinding circuit selection. High Pressure Grinding Rolls (HPGR) may have considerable advantages over autogenous milling circuits but test work and evaluating different grinding circuits must be done to establish which circuit provides the best option. The Hawsons low grade magnetite project provides a specific example in using thorough test work, careful simulations and detail in design to allow the evaluation of grinding circuit performance. Fortunately tools are available for this evaluation and the test work for the Hawsons iron ore is detailed in this paper. The results of the test work and simulation indicates that although the Hawsons iron ore is competent in terms of autogenous milling given the medium hardness of the ore in the Drop Tests, the HPGR testing indicates that the use of HPGR is likely to use half the energy requirement. This is a characteristic of the Hawsons iron ore and is likely to be valid for other iron ores in the Braemer Iron Formation which stretches from the Hawsons lease in New South Wales to Razor Back and South Dam in South Australia. The Hawsons magnetite iron ore resource has been known from the early 1960s when CRA explored the area. However, the Pilbara region was discovered in this era and this led to development of high grade direct shipping grade hematite ores. Carpentaria Exploration (Capex) rediscovered the magnetite whilst undertaking multi-commodity exploration in the Broken Hill district during 2009. The magnetite grade was found to be low at around 16% Davis Tube mass recovery (DTR) but despite this, the project provides a positive business case as a result of the following key factors:• The mineralisation is particularly soft for crushing and grinding with a Bond Work Index of 6.3.• The mineralised rock fractures easily preferentially along grain boundaries.• The location is close to mining infrastructure and a skilled workforce in Broken Hill, 60km to the North East.• It is close to road and rail infrastructure for access to the resource.• The only contaminant in concentrate is silica.• The resource contains many billions of tonnes of magnetite mineralisation.The low Bond Work index is a key characteristic of the Hawsons iron ore and the use of autogenous milling and HPGR circuits is compared to other iron ores indicating why the plant energy requirement may be just as important as iron ore grades.

Aug 8, 2011

By J J. Eksteen

"In the platinum group metals (PGM) industry, smelters are often perceived as ‘black boxes’ that have to process concentrates of variable mineralogy and originating from multiple mines operating along strike or downdip from each other in the same geographical region, as is the case in the Western Bushveld of South Africa. However, PGM smelter complexes are much more sensitive to concentrate blend variation than most metallurgists or geologists would anticipate.This paper will review the impact of concentrate blends deriving from a number of concentrators processing ores from the Bushveld Igneous Complex (BIC) on smelting and converting behaviour and the constraints they impose upon smelter operation. Significant variation in behaviour is not only noted between the main mineralisation types, namely the Merensky and Upper Group 2 (UG2) reefs, but also within a given reef type. Due to the way PGMs are mineralised and finely disseminated, concentrators have to grind concentrates ultra-fine to liberate the PGMs and provide a sufficiently high-grade concentrate for smelting.The chemistry and mineralogy of the concentrates vary with the type of mine (open pit versus underground), the reef type and varies along strike; this variation impacts numerous key smelter parameters such as the base metals content, the matte fall, the converter aisle capacity, sulfur dioxide gas production, furnace matte and slag temperatures, and the accumulation of high smelting chromite-spinel minerals that has led to significant smelter upset conditions such as catastrophic furnace hearth and sidewall failures.Due to the production intensity of smelting, ie its high rate of processing per unit volume, smelters have been found to be far more sensitive to rapid variation in concentrate properties, highlighting the need for adequate blend control.CITATION:Eksteen, J J, 2016. Impact of Mineralogical Variation on the Smelting of Cu-Ni-PGM Concentrates, in Proceedings The Third AusIMM International Geometallurgy Conference (GeoMet) 2016, pp 321–326 (The Australasian Institute of Mining and Metallurgy: Melbourne)."

Jun 15, 2016

MIM Holdings Limited (MIM) committed to a long-term future in copper in the Carpentaria Mineral Province in 1996 with the development of the Enterprise and Ernest Henry mines, and an upgrade of MIMsÆ copper smelter in Mount Isa, and copper refinery in Townsville. The upgrade of the copper smelter will enable the smelter to produce 250 000 tpa of copper with a single primary smelting unit, the copper Isasmelt furnace installed in 1992. The copper Isasmelt furnace currently smelts concentrate containing 145 000 tpa copper. To reduce the risk associated with the copper smelter upgrade, and shorten the ramp up time of the upgraded smelter, two high-rate trials were undertaken in 1997 to demonstrate that the copper Isasmelt furnace could, with an additional 525 tpd of tonnage oxygen, treat concentrates containing 250 000 tpa copper. During the trials the concentrate feed rate to the Isasmelt furnace was increased from its average rate of 100 tph of concentrate to 150 tph of concentrate. The increase in the capacity of the furnace was achieved by: supplementing the installed tonnage oxygen capacity of 595 tpd with an additional 300 tpd supply; reducing variability in the furnace feed rate from ¦25 per cent to less than ¦5 per cent; and reducing tramp air. The Isasmelt furnace smelted concentrate at a rate of 140 to 150 tph for four hours 25 minutes during Trial 1, and at a rate of 150 tph of concentrate for five hours and 30 minutes during Trial 2. The oxygen utilisation was 100 per cent during the trials with the process remaining under control. As a result of the development work for the trials the average treatment rate of the furnace under normal conditions increased by ten per cent.

Jan 1, 1998

By A R. Cruden, M S. Egglseder

The formation of banded iron formation (BIF)-hosted high-grade iron ore deposits in the Hamersley Province, Australia has been debated for decades. Although a structural control on their formation is widely accepted, only the roles of large-scale structures (ie faults and folds) formed by brittle and ductile deformation have been considered to date. However, the deformation of multilayered sequences in the Hamersley Province (ie the Mount Bruce Supergroup), with their contrasting rheological behaviours, creates challenges for understanding the complex structural framework at regional to deposit scales, including the significance for mineralisation and orebody location.Our approach to understanding these complex relationships is based on an established principle of structural geology, which is to extrapolate microscale structures to deposit-scale structures and vice versa, supported by state-of-the-art 3D modelling techniques. Complexly deformed hand samples associated with the Mount Tom Price iron ore deposit (Hamersley Province, operated by Rio Tinto Ltd) have been analysed and visualised with 3D explicit (Move™, Midland Valley Ltd) and implicit modelling software packages (Leapfrog™, AranzGeo Ltd) and supported by microstructural investigations (optical and X-ray microscopy) and fieldwork.Our study shows that the deformation of BIFs in the central Hamersley Province is controlled by microscale chemical processes (ie dissolution precipitation creep) due to rheological differences between iron oxide and chert layers, contrasting fluid flow properties and layer thickness variations. These processes account for the removal of vast amounts of silica, its transportation by diffusion in the presence of a fluid and reprecipitation at low-pressure sites. Consequently, this structural framework was used by different upgrading fluids to form iron orebodies.CITATION:Egglseder, M S and Cruden, A R, 2015. Evolution of a structural framework for high-grade iron ore deposits, in Proceedings Iron Ore 2015, pp 255–262 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Jul 13, 2015

By L Keeney

Geometallurgy is not a mechanism designed to replace current geological and metallurgical deposit characterisation but an added component required to understand orebody variability. Current geological characterisation approaches focus on understanding deposit genesis aimed at identifying new mineralised regions and metallurgical characterisation determines parameters required for engineering design of processing plants. Geometallurgical integration has the task of drawing on aspects of each stream to enable spatial mapping of performance indices across a deposit. It is a complimentary activity to current practice and adds a new dimension to orebody knowledge. Although the notion of conducting an integrated geological and metallurgical analysis is not a new concept, limited publications exist in the field proposing integrated methods, with several technical issues existing in current geometallurgical integration approaches. An iterative method drawing on existing geological and metallurgical characterisation and new measurement tools and inputs is proposed. The integrated method comprises of discrete components dynamically linked together through iterative feedback loops with each component having a specific purpose achieved by applying a range of data acquisition technologies and statistical tests. The key outcomes of the method are:overcoming the lack of optimised geometallurgical domains relevant to processing performance,development of a rigorous approach to identifying parameters controlling processability,improved protocols for relevant metallurgical test data and sample selection, andability to generate greater data support to enable rigorous geostatistical modelling of geometallurgical attributes. This paper outlines the integrated mapping and modelling method developed and draws on multiple data acquisition technologies applied within the AMIRA P843 GeMIII Project.

Sep 5, 2011

The Sishen iron-ore-mine produces 3 pro= ducts with size gradings of 25 to 8 mm, 11 to 5 mm and 5 to 0,2 mm and has a product capacity of 20 Mtpa. Ore is produced at an average Fe content of 66,2 % Fe for the two coarser pro= ducts and an average of 65,3 % Fe for the fine ore product. Three stage crushing is practised to re= duce the run-of-mine ore to -90 mm after which it is screened into five different size frac= tions of which four are heavy media beneficiated in static baths and cyclones and the fifth frac= tion discarded to slimes dams. A final three stage crushing and screening process is employed to produce the final three Products which are homogenized by means of blen=ding systems nrinr to despatch.

Jan 1, 1987

By R G. Reddy

The complexity in treating environmentally harmful impurities such as As, Sb, Bi from the base metal matte or metal in the smelting stage are responsible for the high cost of the refining process. The impurity capacities (such as arsenic, antimony and bismuth) of slags were calculated a priori using Reddy-Blander (RB) model. The capacity predictions were for a wide range of matte and slag compositions in copper smelting conditions. The calculated impurities capacities and impurity distribution ratios results are in good agreement with the available experimental and industrial slags data. The a priori knowledge of impurities is useful for reduction of energy consumption and enhanced environmental control in the current and future non-ferrous metal processes.

Aug 21, 2024

By K L. Avarmaa, M A. Rhamdhani, A D. P Putera, H T. B M Petrus, G A. Brooks

Multistage slag treatment was carried out to investigate boron removal from metallurgical grade silicon. Boron is known to be one of the most difficult elements to remove from silicon. Depending on boron initial concentration, boron can be removed from silicon into certain level; most previous studies have attained final 50–100 ppmw (parts per million weight) boron in the silicon using slag treatment. The term ‘multistage’ refers to a process where the resulting silicon from the initial slag treatment is further refined multiple times using a new slag. In this study, an initial slag composition of 45CaO-45SiO2-10Al2O3 (wt per cent) was utilised. Using this technique, boron level in the silicon was expected to be gradually lowered because the new slag will always have more capacity to reach the equilibrium distribution coefficient (LB = (ppmwB)slag/(ppmwB)Si). The current study was carried out through combined thermodynamic calculations (using FactSage™ ver 8.2, by Thermfact Ltd, thermochemical package) and high-temperature experimentation. The thermodynamic calculations were carried out for reactions at 1500, 1550, and 1600°C for three different scenarios: 1. Single stage refining with varied slag to Si ratio. 2. Multistage refining with slag to Si ratio of 1:1. 3. Multistage refining with slag to Si ratio of 2:1. The multistage experimental work was carried out at 1550°C under inert gas atmosphere each time for 2 hrs, with slag to Si ratio of 2:1. The simulation results indicate that the boron distribution coefficient stay essentially constant with varying slag to Si ratio and temperature, around 1.8 to 2, and the experimental results measured LB in the same range, between 1.5 and 2.5. By harnessing this feature, the simulation showed that boron can be removed by approximately 78 per cent in each stage. The experimental result confirmed that the boron concentration in the silicon was reduced by a total of 97.8 per cent from its initial value using four stages of slag treatment, from initial 270 ppmw down to 6 ppmw.

Aug 21, 2024

By B F. J Kelly

Multiple-point simulation is a family of methods used to generate conditional facies simulations that present the same spatial dependence as a user-defined training image. The main advantage of using a training image is the easy integration of geological concepts in numerical subsurface models. The strength of multiple-point methods is that it allows making the best use of geological expertise when data are too sparse for inferring a model of spatial continuity. However, despite being appealing to geologists, the multiple-point methods are often not used in practice due to the difficulty of building complex 3D training images that contain realistic features.A new framework for multiple-point simulation is presented that involves small and simple training images. The main idea is to use random transformations to expand the range of structures available in the simple training image. The training image is not regarded any more as a global conceptual geological model, but rather a basic structural element of the subsurface. Complex geological structures are obtained, whose spatial structure can be parameterised by adjusting the statistics of the random transformations, based on field data or geological context. In most cases, such parameterisation is possible by adjusting two numbers only. One advantage is that the training images are so simple that they can be easily built even in 3D. The method allows building 3D models that reproduce shapes corresponding to a desired prior geological concept and are in phase with different types of field observations, such as orientation, facies or geophysical measurements. The flexibility of the approach allows for naturally integrating different types of data. This is demonstrated on an example involving seismic data assimilation. Examples where complex 3D folded structures are built from a simple training image and a few orientation data are shown.

Sep 26, 2011

By S Johnston, L Howe, S Stevens

Regulators have an important role to play in fostering and enabling reliable, effective, environmental management. Given this, it is perhaps surprising that there has been relatively little debate on exactly how regulators can enhance their own approaches, in order to promote improved outcomes overall. The authors argue that applying research on High Reliability Organisations, as well as the University of Pennsylvania’s Penn Program Framework for Regulatory Excellence, to regulators has the potential to facilitate positive environmental results across the mine life cycle.

Aug 2, 2023

By D Craw

The Macraes gold deposit is located in a low angle NE dipping (~15¦) shear zone, the Hyde Macraes Shear Zone (HMSZ) cutting greenschist facies metasediments of the Otago Schist. The shear zone is host to large volumes of mineralised schist and relatively sparse quartz veins. Early in the development of this through going shear the schist experienced an episode of hydro fracturing and fluid redistribution. As a result, presently exposed schists within the Intrashear sequence are pervaded by localised swarms of quartz filled pure extension and extensional shear fractures hosting goldscheelite mineralization. The swarms comprise networks or æstockworksÆ of veins that strike NE, are subvertical and form subperpendicular to the shear zone fabric. Vein frequency across strike is typically ~1 m-1, with most veins five to ten centimetres thick and vertically continuous for up to ~20 m. Many of the veins show delicately banded internal textures that are parallel to the vein margins. These bands/lamina indicate incremental growth and are defined by a variation in grain size or alternatively by thin slivers of wallrock schist, though typically uncommon. The majority of veins are overprinted by a weak plastic deformation, although primary structures are largely preserved. This deformation has produced subgrain structures, migration of grain boundaries with sutured margins and undulose extinction of grains. These microstructural observations possibly indicate that temperatures were around 300¦C during stockwork vein formation. In general, the mode of brittle failure that developed the stockwork veins can be directly related to lithological variation. The geometry of stockwork veins indicate that local vertical shortening and extension occurred in the Intrashear sequence of the HMSZ synkinematically with thrusting. Localised extension was mainly parallel to the structural trend; that is, it took place perpendicular to the direction of inferred thrusting. The development of lateral ramps most likely gave rise to extension in parts of the shear zone during continuous shortening to form swarms of stockwork veins.

Jan 1, 2004

By D Adhikary, M Khanal, B Poulsen, H Guo

During longwall mining, the ability to reliably estimate the height of longwall mining-induced fractures is critically important in order to effectively manage mining-induced water inflows and aquifer interference. The extent of mining-induced subsurface deformation, fracture, surface subsidence and aquifer interference is largely controlled by local lithological conditions, mining methods and mine layouts. Therefore for the purpose of quantifying mining impacts it is imperative to incorporate the local geology and mine plan parameters in a mine scale model. Two different numerical models based on 3D continuum (Finite Element Method) and 2D discrete (Discrete element method) domains, respectively, are developed to incorporate the mine site parameters of a study site in Australia. Both the models are calibrated and validated with measured mine data. The continuum numerical modelling method is used to investigate the change in permeability of the overburden strata due to underground extraction considering the full 3D nature of geology and mine layout. The discrete element modelling method is used to simulate the initiation and propagation of fractures, allowing the caving of overburden strata to be represented in a realistic manner, and the permeability enhancement of the strata overlying the goaf to be calculated from first principles. Permeability enhancement in the discrete element model is a direct function of the fractures induced in the deformation/caving process; fracture location, extent, aperture and connectedness. From the numerical modelling under the simulated condition for the mine studied here, it is predicted that the permeability of the caving zone from the goaf to five times the mining height increases by almost six orders of magnitude. In the fracture zone to approximately ten times the mining height, the change in permeability is predicted to be around two to three orders of magnitude. At the height of more than ten times the mining height the magnitude of permeability increased is predicted to be one to two orders of magnitude. CITATION:Khanal, M, Poulsen, B, Adhikary, D and Guo, H, 2018. Modelling of height of connective fractures during longwall mining, in Proceedings The Fourth Australasian Ground Control in Mining Conference (AusRock), pp 254–265 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Nov 30, 2018