ISSN (0970-2083)

All submissions of the EM system will be redirected to Online Manuscript Submission System. Authors are requested to submit articles directly to Online Manuscript Submission System of respective journal.

PYROMETALLURGICAL PROCESSING OF SULFIDE POLYMETALLIC RAW STUFF IN AUTOGENEOUS FORCED MODE

Georgy Vladimirovich Konovalov*, Tatiana Reonaldovna Kosovtseva and Victor Mikhailovich Sizyakov

Saint Petersburg Mining University, Russian Federation, 199106, St. Petersburg, 21st Line, 2, Russia

Corresponding Author:
Georgy Vladimirovich Konovalov
E-mail: kgv_@mail.ru

Received date: 06 May, 2017; Accepted date: 08 May, 2017

Visit for more related articles at Journal of Industrial Pollution Control

Abstract

An innovative approach to development of pyrometallurgical autogeneous processing of sulfide copper nickel raw stuff is proposed, which is based on radically new air blasting into melt. Spatially oriented air blasting makes it possible to combine areas of heat generation and exchange in operation space of autogeneous assembly to the greatest possible extent. In order to implement the spatially oriented air blasting a new design series of tuyeres is proposed, which provide generation of ultrasound jets outflowing at any preset angle to bath surface. Experiments with cold models and fire assemblies confirmed possibility to increase specific blasting load in two times and higher in comparison with existing methods of air blasting. The following technical effects are assumed after implementation of the scientific research results: decrease in specific power consumptions; increased time of campaign lives of autogeneous assembly; suppression of dust entrainment and reduction of emissions of sulfur dioxide into environment.

Keywords

Air blasting, Melt, Matte, Concentrated product, Sulfides, Copper, Nickel, Tuyere

Introduction

According Pyrometallurgical processes are the major links of production of copper, nickel, and associated metals (Schlesinger, 2011). Autogeneous processes on the basis of chemical energy of sulfide ores and various concentrated products are the preferred approach to development of pyrometallurgy. Potentials of autogeneous processes are not implemented completely, which is confirmed by general features peculiar for all autogeneous assemblies:

1. Air blasting per unit of operation volume is restricted by occurrence of bulk emissions and ranges from 5 to 14 m³/(min·m³) for various alternatives of air blasting (Shalygin, 2003; Shalygin, 2006).

2. Increased content of magnetite and non-ferrous metals in slags of autogeneous smelting evidences incomplete reactions of slag generation upon formation of final products (Schlesinger, 2011; Shalygin, 2006; Sridhar, 1997).

3. Non-stable pattern of gaseous dynamics in furnace and beyond promotes generation of hard melting skulls, preventing recycling of dust gas mixture (Greiver, 1999; Krupnov, 2013);

4. Isolated positions of heat generation and heat exchange lead to forced heat losses with cooling elements and decrease in campaign lives (Shalygin, 2003; Seregin, 2001);

5. Autogeneous assemblies require for additional power sources, which is accompanied by average unrecoverable heat losses into environment in amount of 15-20% (in terms of total input) (Krivandin, 2002; Lisienko, et al., 2005).

In this or that extent these interrelated critical properties are characteristic to Russian autogeneous processes (PV, KIVCET, AAP, and others) (Lisienko, et al., 2005; Reznik, et al., 2003; Vanyukov, 1988; Bulatov, 2014), especially pronounced in decreased quality of processed ore and involvement of low exothermic and hard melting products of metallurgical process (Eroshevich, 2012; Boduen, 2014).

As a consequence of the researches by Shalygin a radically innovative jet swirling apparatus (JSA) was proposed (Shalygin, 2003; Shalygin, 2006; Shalygin, et al., 2006), its main concept is possibility of multiple increase in ultimate amount of air blasting per unit of operation space of apparatus, provided that intensive regular heat and mass exchange in molten bath is maintained created by spatially oriented air blasting.

Experimental Procedure

Based Theoretical studies were based on known physicochemical regularities of processes occurring at the stage of oxidizing pyrometallurgical processing of sulfide polymetallic raw stuff. Mathematical simulation was performed using MathCad and ANSYS software. Experimental studies were performed on cold and fire models in laboratories of Mining University, on large scale pilot facilities at OAO Kola MMC. Chemical composition was determined using X-ray spectral analysis. Composition of gaseous phase was determined by chemical analysis. Technical measurements were performed using regular thermocouples, thermal camera, gauges equipped with disk diaphragms.

Insights Results and Discussion

The main task of pyrometallurgical processes is consecutive removal of iron sulfide from complex sulfide system. Oxidation of iron sulfides is accompanied with heat release, hence, the possibility of specified autogeneous pyrometallurgical process in an apparatus of certain design is revealed by solution of extended equation of heat balance (Shalygin, 2003):

equation (1)

where X is the heat consumption by blend; equation is the heat content of process product; Qext is the external heat loss; Qend is the heat of endothermic reactions; Qblend is the heat content of initial blend; Qblast is the heat content of air blasting;

Qex is the heat of exothermic reacting considered in material balance.

Remarkable attention in solution of minimization of power consumption and achievement of autogeneous essence of pyrometallurgical facility should be paid to its external heat losses Qext. External heat losses were determined using thermal imaging of main structures of several pyrometallurgical facilities Table 1 involved in technological flowchart of copper production.

No. Item tmax, °C tmin, °C tavg, °C
1. Vanyukov smelting furnace
(matte smelting)
Slag end 159 47 93
2. 153 94 132
3. Tuyere water jackets 52 49 52
4. 51 49 50
5. Brickwork 52 47 49
6. 58 46 56
7. Water jackets 46 42 44
8. 40 35 39
9. Uptake 41 35 40
10. 205 156 172
11. Roof 211 124 174
12. 213 114 180
13. Uptake 202 165 182
14. 84 70 76
15. Matte mixer 211 156 167
16. 190 155 183
17. Horizontal converter
(conversion)
External wall 326 269 309
18. 333 201 296
19. Hood of tuyere area 230 211 220
20. 249 184 231
21. Anodic furnace
(fire refining)
Lining 283 237 257
22. 375 119 290

Table 1. Temperature distribution over external surface of metallurgic assemblies of copper production process flow (matte smelting–conversion–fire refining)

On the basis of performed measurements total amount of heat losses into environment was determined per one ton of processed blend, taking into consideration peculiar features of pyrometallurgical assembly and its dimensions:

equation

where Fsurf is the uncooled surface area of the apparatus with account of ribbing, m²; Fhole is the surface area of open holes, m²; Wcool is the water consumption for cooling of water jackets, kg/t; Wtuyere is the water consumption for tuyere cooling, kg/t; τeq is the balance time, h.

The equilibrium time is a function of specific performance of pyrometallurgical apparatus and external surface area of furnace:

equation (3)

where a is the specific performance, t/m²⋅day, F is the surface area of representative cross section of furnace, m².

Hence, the balance time as a factor determining heat losses into environment at preset daily performance, other conditions being equal, is inversely proportional to specific performance and representative dimensions of the apparatus.

Specific performance of pyrometallurgical apparatuses

Based on the principle of forced supply of gaseous oxidizer, depends undoubtedly on amount of air blasting supplier per unit time. Possibility to increase performance of autogeneous apparatuses is limited by physical impact of blasting jets on molten medium. Combination of dynamic head, forces of thermal expansion of gases in the melt, and buoyancy forces of gas bulk ascending from the melt result in occurrence of massive melt ejection from operation space of the apparatus. Method of air blasting stipulates certain ultimate blasting load.

Upon experiments on models and industrial facilities the following values of ultimate blasting load were determined per 1 m³ of operation space of apparatus (Table 2) (Shalygin, 2003; Shalygin, 2006; Shalygin, 2006). Possible increase in blasting load and, hence, in specific performance will significantly influence on the structure of heat balance of the apparatus, which in particular will expand possibilities of autogeneous oxidizing processing of sulfide polymetallic raw stuff of non-ferrous metallurgy.

Autogeneous unit Air blasting type Blasting load, m³/(min×m³)
Bessemer converter Bottom distributed 4-5
Vertical oxygen converter Upper non-submerged 6-8
Kaldo converter Inclined non-submerged 6-9
Horizontal converter Side-type submerged unilateral 7,5-13
Vanyukov smelting furnace Side-type submerged bilateral 8-10
Jet swirling unit Spatially oriented 50-85

Table 2. Ultimate specific blasting load on various units of autogeneous type reduced to normal conditions (T = 273 К, P = 101315 Pa)

Reduction of external heat losses due to increased performance can compensate certain decrease in heat of exothermic reactions related with processing of ore with decreased sulfur content and provide processing of such ore in autogeneous mode without additional supply of any fuel.

Heat and mass exchange of jet swirling apparatus

Separated positioning of heat generation and heat exchange areas leads to thermal overloads of furnace lining near local sites of intensive heat release and decrease in campaign lives (Table 3).

Autogeneous apparatus High temperature area Temperature in heat generation area, °C Campaign lives, days
Vertical oxygen converter Blasting crater > 2500 330
Horizontal converter Near-tuyere area 1500-1750 60
Vanyukov smelting furnace 180-330
Suspension smelting Reactive blend 1700-1800 330

Table 3. Performances of autogeneous units

Air blasting procedure into melt via tuyeres located in copper water jackets of Vanyukov smelting furnace cannot be considered as successful, which is evidenced by short operation lifetime with air blasting according to observations in 9 months (Table 4).

No. Failures Time, hours
1 - scheduled preventive operations 37.58
2 - destruction of water jackets 97.42
3 - power failure 30.17
4 - shutdown of air oxygen mixture supply 27.67
5 - shutdown of gas exhausting system 90.42
6 - technological and transport shutdowns 31.33
7 - major overhaul 974.00
Total 1288.59
Calendar period of observations - 6288
Coefficient of operation with blasting without accounting for major overhaul period - 0.95
Coefficient of operation with blasting in calendar period - 0.80

Table 4. Downtime indices of Vanyukov smelting furnace, Balkhash copper smelter (OAO Kazakhmys)

The main reasons of destruction of water jackets were the factors of chemical corrosion with generation of sulfur containing copper compounds decreasing mechanical strength and heat conductivity as well as thermal, caused by variations in composition of exothermic blend and breakage of skull protection.

Regular ordered mass transfer stipulates uniformity of thermal and chemical fields in the melt bulk and provides the most favorable conditions for technological process.

Mathematical interpretation of swirling effect of liquid as a consequence of application of tangential stresses to its surface:

equation(4)

assuming in the first approximation that the tangential stresses on free surface are concentrated in narrow ring with the radius R0 = r0⋅R and width δ << R, then

equation (5)

Substituting equations (4) and (5) in stationary mode we obtain the velocity field:

equation (6)

where r0 = R0 /R, h= H/R are the dimensionless geometrical characteristics of cylinder; R is the radius, m; R0 is the radius of application of tangential stresses, m; H is the height of liquid layer, m; βk is the zero of Bessel function, υ is the relative velocity.

Using the MathCad mathematical model the maximum velocity of liquid swirling was determined near intensive physiochemical interactions (area in the vicinity of bath surface) at the radius r0 of application of tangential stresses to its surface equaling to 0.7 (Fig. 1) which corresponds to the results of physical simulation (Konovalov, 2002).

icontrolpollution-dimensionless-rate

Fig. 1 Dependence of dimensionless rate in z and r coordinates at h = 0.9.

Experimental validation of the JSA concept in melt was performed for various compositions of copper nickel material in the course of fire experiments in the range from 1.0 dm³ to 1.0 m³ and specific blasting load from 26 to 90 m³/min⋅m³ (Table 5).

No. Material Content, % Blasting unit Volume,
Tip diameter,
mm
Specific blasting load, m³/(min×m³)
Cu Ni Fe S SO2
1 Cu-Ni matte 6.0 10.5 48.3 28.6 9-12 6 mini-tuyeres 0.001 1.6 50
2 8.3 13.0 35.0 N/A 5-11.5 6 tuyeres 1 8 26
3 Cu concentrate 69.2 3.7 3.5 20.5 11-12 6 mini-tuyeres 0.001 1.6 65
4 68.7 3.4 4.2 21 10-12 unit (6 tips) 0.001 1.6 90
5 68.7 3.4 4.2 21 N/A unit (6 tips) 0.02 4 60
6 Ni concentrate 2.7 72.8 4.8 17.6 N/A 6 tuyeres 0.03 4 35

Table 5. Fire experiments in the mode of jet swirling apparatus

The pattern of steady process was visually observed upon intensive swirling of the melt (matte, copper and nickel concentrated products from converter matte plant) and without emissions. Blowing process and content of sulfurous anhydride (5-12% upon air blasting) in gases confirmed possibility of their efficient recycling, herewith, cyclone effect was observed above the melt preventing emissions of melt and dust.

Therefore, large scale experiments made it possible to obtain positive response to initial question concerning possibility of non-submerged blasting by means of high pressure spatially oriented jets to create ordered mass and heat transfer, meeting the requirements of kinetics of physicochemical interactions in heterogeneous system.

Radial axial tuyeres of new design

Peripheral spatially oriented air blasting is implemented by means of new design of tuyere, its development was based on the following required conditions ((Konovalov, 2015; Shalygin and Konovalov, n. d.):

– Possibility of multiple reduction of dimensions of tuyere fire zone and decrease in water consumption for its cooling;

– Cooling under depressurizing;

– Possibility to decrease the height in case of horizontal positioning of several tuyeres;

– Creation of regular swirling motion of melt.

Such requirements are met by radically new tuyere design (Fig. 2).

icontrolpollution-Design-radial-axial

Fig. 2 Design of radial-axial tuyere: 1 – cap, 2 – tuyere head, 3 – blasting tip, 4 – central tube, 5 – external tube, 6 – blasting cartridge, 7 – seal, 8 – counter nut, 9– intermediate tube, 10 – water cooling cartridge.

Analytical mathematical description of channel of complex shape is difficult, hence, the coefficients of hydraulic resistance were determined using experimental bench for hydrodynamic tests (Shalygin, et al., n. d).

Experimentally determined hydraulic resistance of the tuyere cooling system was significantly lower than 1 atm which proves the validity of previous assumption of possibility to apply explosion-proof cooling under depressurizing.

Heat loads on tuyere were calculated in ANSYS software (ANSYS CFX-Solver Theory Guide ANSYS, Inc. 2009) using SST model (Shear stress transport) (Menter, et al., 2003). As a consequence, temperature distribution of cooling agent was determined. It is established that under extreme conditions of heat transfer water consumption to tuyere with tip diameter of 20 mm should be about 0.018 m3/s.

Skull protection of lining is the most efficient protection against chemical and thermal destruction of operation surface of pyrometallurgical apparatus contacting with molten medium: both of brickwork and of metal surface.

Contrary to the existing variants of jet air blasting (vertical non-submerged, horizontal submerged) the proposed approach to creation of spatially oriented non-submerged jets distributed over periphery of cylindrical body of the apparatus and creating controlled swirling of melt bulk due to tangential stresses is the most favorable also for formation of steady skull.

Such mode was preliminary estimated using specialized laboratory facility (Shalygin, 2003). The estimation was performed for qualitative comparison of heat exchange conditions between melt and wall upon vertical air blasting by single tuyere and by JSA. Both experiments were performed maintaining similar air pressure and flow rate, the same temperature, the same mode of water supply. The pattern of skull generation was obviously different. Upon vertical air blasting skull is heterogeneous in thickness and structure at various areas of the water jacket. Either it is completely absent, or it is rather thick. This can be attributed to already mentioned randomness of behavior of blasting crater, formed by vertical jet, its oscillations, displacements in transversal direction. In addition, significant vertical melt splashes are observed, emitted above the model body, upon single vertical air blasting. When air blasting is performed in the mode of JSA, vertical splashes of the melt are nearly absent, and skull is distributed over total surface of jacket wall as uniform dense layer.

Melt swirling under the impact of spatially oriented jets creates ordered mode of convective thermal action of melt onto cylindrical apparatus wall. Temperature of the swirling melt contacting with the wall due to bath mixing satisfies certain averaged value, which is significantly below the temperature developed in the area of blasting craters. At sufficient heat conductivity of refractory wall arrangement of water cooling jackets in it can provide generation of steady skull on fire surface, which can promote long operation lifetime at high thermal intensity characteristic to autogeneous process.

Cyclone effect of gaseous phase above melt. On the basis of experiments with cold and fire models of JSA the cyclone effect, visually observed in gaseous phase above melt, was mathematically simulated using ANSYS CFX package. Total hydrodynamic characteristic of JSA is illustrated in (Figs. 3–5) by calculations using control volumes.

icontrolpollution-Behavior-liquid

Fig. 3 Behavior of liquid bath.

icontrolpollution-Dynamics-gaseous

Fig. 4 Dynamics of gaseous phase at the height of radial axial tuyeres.

icontrolpollution-radial-axial-tuyeres

Fig. 5 Dynamics of gaseous phase above radial axial tuyeres.

Conclusions

The Physicochemical concept of thermal balance is developed in detail, enabling estimation of each item of thermal balance and evaluation of possibility of reasonable control of these items.

Interrelation between external heat losses and specific performance of metallurgical apparatus is demonstrated. The factor, which restricts increase in specific performance, is the method of air blasting.

In order to perform autogeneous processing of polymetallic raw stuff in forced mode a radically innovative, spatially oriented air blasting was proposed enabling combining of areas of heat generation and exchange in operation space of autogeneous assembly to the greatest possible extent.

Experiments with cold models and fire assemblies confirmed possibility to increase specific blasting load in two times and higher in comparison with existing methods of air blasting.

Flow dynamics pattern in the proposed apparatus of jet rotation under the action of tangential stresses created by kinetic energy of distributed spatially oriented jets is established.

A new design of autogeneous assembly for processing of sulfide polymetallic raw stuff is proposed on the basis of spatially oriented jets creating swirling motion of melt in overall apparatus space.

A radically new design series of tuyeres is developed, which provide generation of ultrasound jets outflowing at any preset angle to bath surface.

The tuyeres, possessing very moderate hydraulic resistance of cooling channels, correspond to conditions required for underpressure cooling system, and water consumption for cooling of set of 6 tuyeres is lower by three times in comparison with vertical single tuyere.

Possibility to create steady skull protection under regular controlled heat and mass transfer in molten bath occurring under the action of ultrasound spatially oriented jets.

References

  1. Altushkin I.A. 2012. Innovacii v metallurgii medi na primere realizacii proekta rekonstrukcii ZAO “Karabashmed'” Chast' 1, 2. (Innovations in metallurgy of copper on the example of the project of reconstruction "Karabashmed" Part 1 & 2). Nonferrous metals. 8 : 25-40.
  2. ANSYS CFX-Solver Theory Guide ANSYS, Inc. 2009, pp. 625-632. Date Views 23.02.2017 http ://www.ansys.com/
  3. Avtogennye processy v cvetnoj metallurgii. (Autogenous processes in non-ferrous metallurgy). (1991). Moscow : Metallurgy, pp. 413.
  4. Bulatov K. V. 2014. Pererabotka polimetallicheskih koncentratov v plavil'nom agregate “Pobeda” OOO “Mednogorskij medno-sernyj kombinat”. [Processing of polymetallic concentrates for smelting unit "Pobeda" LLC "Mednogorsky copper and sulfur plant"]. Non-ferrous metals. 10 : 39-44.
  5. Boduen A. J. 2014. Vozmozhnost' primeneniya gidrometallurgicheskih metodov dlya povysheniya kachestva sul'fidnyh mednyh koncentratov. [The possibility of using hydrometallurgical techniques to improve the quality of copper sulphide concentrates]. Non-Ferrous Metals. 11 : 37-41.
  6. Eroshevich S.Y. 2012. Analiz izmeneniya sostava sul'fidnogo rudnogo syr'ya, pererabatyvaemogo v processe vzveshennoj plavki, i tekhnologicheskie osobennosti raboty v usloviyah snizheniya ego teplotvornosti. [Analysis of changes in the composition of the sulfide ore raw materials, processed by flash smelting and technological features of operation with decrease in its heat value]. Nonferrous metals. 9 : 13-20.
  7. Greiver, T.N. 1999. Osnovy metodov postanovki i resheniya tekhnologicheskih zadach cvetnoj metallurgii. (Basics of staging methods and technological solutions of non-ferrous metallurgy problems). Moscow : The ore and metals, pp. 147.
  8. Konovalov G.V. 2002. The peculiarity of mass transfer under the influence of Francis blowing jets. Mining Institute Notes. Mining Institute (Technical University). 150 : 120-122.
  9. Konovalov G.V. 2015. Studies of the hydroaerodynamic characteristics of a radial-axial lance. Metallurgist. 59(7-8) : 631-636.
  10. Krivandin V.A. 2002. Teplotekhnika metallurgicheskogo proizvodstva. [Heat metallurgical production]. Volume 2. The construction and operation of furnaces. Textbook for high schools. Moscow : MISA, pp. 736.
  11. Krupnov, L.V. 2013. Mekhanizm formirovaniya tugoplavkoj nastyli v pechah vzveshennoj plavki Nadezhdinskogo metallurgicheskogo zavoda. [The mechanism of formation of refractory flooring in the flash smelting furnace at the Nadezhda Metallurgical Plant]. Nonferrous Metals. 2 : 46-51.
  12. Lisienko V.G., Schelokov Y.M. and Ladygechev M.G. 2005. Hrestomatiya ehnergosberezheniya :. Reader of Power Saving : Reference book : In 2 books. (Book 2). Moscow : Teplotekhnik, pp. 768.
  13. Menter F.R., Kuntz M. and Langtry R. 2003. Ten Years of Experience with the SST Turbulence Model / Turbulence, Heat and Mass Transfer 4. Begell House Inc. 625-632.
  14. Reznik I.D., Ermakov, G.P. and Schneerson J.M. 2003. Nickel : in 3 volumes. Moscow : Science and Technology. 608.
  15. Schlesinger, M.E. 2011. Extractive Metallurgy of Copper. Oxford : Elsevier. 455.
  16. Shalygin L.M., Telyakov N.M., Konovalov G.V., Shalygin L.M., Telyakov N.M. and Konovalov G.V., (n. d.). Pat. 2236474 RF, C22B15 / 06 IPC. Sposob konvertirovaniya medno-nikelevyh shtejnov. [The process of converting copper-nickel matte]. The applicant and the patentee St. Petersburg State Mining Institute. No. 2003120316 / 02, appl. 02.07.03 published. 20.09.04.
  17. Shalygin L.M. 2006. Perspektivnoe napravlenie avtogennoj pererabotki sul'fidnyh rud i koncentratov na osnove prostranstvenno orientirovannogo kislorodnogo dut'ya. [Perspective direction autogenous processing of sulphide ores and concentrates, based on space-based oxygen blowing]. Nonferrous metals. 1 : 12-18.
  18. Shalygin L.M. and Konovalov, G.V. (n. d.) Pat. RF 2244020, C21C5 / 48 IPC. Furma dlya produvki rasplavov. [Lance for blowing the melt]; the applicant and the patentee St. Petersburg State Mining Institute. No. 2003136221/02; appl. 15.12.03; publ. 10.01.05.
  19. Shalygin L.M., Konovalov G.V. and Kosovtseva T.R. (n. d.). Pat. RF 2349648, C21C5 / 48 IPC. Furma dlya produvki rasplavov. [Lance for blowing the melt]; the applicant and the patentee St. Petersburg State Mining Institute. No. 2007109567/02; appl. 15.03.07; publ. 20.09.08.
  20. Shalygin, L.M. 2003. Teplogeneraciya i teploperenos v avtogennyh metallurgicheskih apparatah raznogo tipa. [Heat generation and heat transfer in the autogenous smelting machines of various types]. Nonferrous Metals. 10 : 17-24.
  21. Shalygin, L.M. 2006. Usloviya podachi dut'ya v rasplavy i razrabotka sredstv intensifikacii dut'evogo rezhima. [Terms of supply in the blast melts and the development of means of intensification of the blow mode]. Notes of the Mining Institute. SPMI (TU). 169 : 231-237.
  22. Seregin, P.S. 2001. Issledovanie gazovoj dinamiki vnutripechnogo prostranstva, pylevynosa i nastyleobrazovaniya v pechi Vanyukova s ispol'zovaniem metoda fizicheskogo modelirovaniya. [Research of gas dynamics of interfurnace space, and dust and skull formation in Vanyukov furnace using the method of physical modeling] : Abstract. Doctoral thesis. St. Petersburg, pp. 23.
  23. Sridhar, R. (1997). Copper losses and thermodynamic considerations in copper smelting. Metallurgical and Materials. 28(2) : 191-200.
  24. Vanyukov A.V. (1988). Plavka v zhidkoj vanne. [Melting in a liquid bath]. Moscow : Metallurgy. 208.
  25. Tsymbulov L.B. 2014. Pirometallurgicheskie sposoby pererabotki sul'fidnyh mednyh koncentratov Udokanskogo mestorozhdeniya. [Pyrometallurgical processing of Udokan sulphide copper concentrates]. Non-Ferrous Metals. 9(861) : 29-37.

Copyright © 2024 Research and Reviews, All Rights Reserved