Waste Dump Leaching Copper

Waste Dump Leaching Copper

To help maintain an adequate supply of minerals to meet national economic and strategic needs, the Bureau of Mines developed a mathematical model that describes acidic ferric sulfate leaching of copper from sulfide ores. The leaching model can aid in maximizing minerals and metals recovery from primary sources by giving operators and management a better understanding of the copper leaching process. This will aid in improving dump leaching operations. The model is based on a steady-state approximation of the continuity equation for diffusion of ferric ions in the rock pores and the intrinsic leaching kinetics of the various sulfide mineral particles. Ore characteristics included in the model were ore particle-size distribution, particle shape, grade, mineral content, mineral particle size, porosity, and tortuosity. The continuity equation was used with ferric ion concentration gradients and a fixed value for the diffusion coefficient. Except for effective diffusivity, each of these factors was either measured experimentally or estimated from theoretical considerations. The effective diffusivities were determined by fitting calculated leaching curves with experimental leaching data. Vertical variations in leach solution composition and temperature within a dump are not considered in this model. The model accurately predicted copper extractions for several large-scale tests during nearly 500 days of leaching.

In the domestic mining of million of tons of copper ore, nearly 3X more million tons of waste were moved and discarded. This waste contained copper at a grade too low to permit profitable copper extraction using modern milling techniques. Some of the waste is leached with acidic solutions that percolate through the dump and dissolve a portion of the copper. Dump leaching is an important supplement to open pit-mill production of copper, providing about 1.0 percent of the total copper recovered. Because dump leaching is slow and incomplete, much copper remains in the dumps after leaching has been terminated. Research to improve leaching efficiency is carried out by the Bureau of Mines as part of its mandate to assure adequate domestic supplies of essential metals.

Mechanisms of the dump leaching art are not well known, and a better understanding can improve the operation’s effectiveness. The careful delineation of leaching mechanisms followed by mathematical incorporation in a quantitative model may assist in dump leaching optimization. This report describes the development by the Bureau of a dump leaching model for determining the response of copper sulfide minerals to many of the conditions in large dumps leached with acidic ferric sulfate solution. The model can provide an upper limit for copper extraction from a dump containing copper as sulfide minerals. Hopefully, application of the model will be a step towards achieving more rapid, more complete, and less expensive recovery of copper from leach dumps.

Because of the increasing importance of dump leaching for augmenting copper production, more investigators are attempting to unravel the complex mechanisms involved. Results of these investigations are being presented as mathematical models. For example, a shrinking core model was used by Madsen, Braun, and Roman to describe copper leaching. This model considered a core of unreacted ore that was separated from a reacted crust by an inward-moving reaction zone. The leaching kinetics involved a steady-state diffusion of reactants through the reacted portion of the ore fragment followed by a chemical reaction at the fresh surface of the mineral. When the chemical reactions were rapid compared with the diffusion rate of the reactants, the model, which assumed a constant-thickness reaction zone, satisfactorily described oxygen leaching of chalcopyritic ore at elevated temperatures and pressures, ferric sulfate leaching of chalcocitic ore at ordinary temperatures and pressures, and leaching of oxide copper ores. This model, however, cannot satisfactorily predict copper leaching when chemical reactions are slow relative to the diffusion rate of the oxidant in which the reaction zone thickens with leaching time. Thus, a different model must be used to describe the ferric sulfate leaching of copper at ordinary temperatures and pressures from ores containing refractory minerals such as chalcopyrite.

A more rigorous model using the continuity equation was devised by Bartlett for describing the copper leaching rate from low-grade chalcopyritic ore. This model involves the kinetics of chalcopyrite grains reacting with oxygen and diffusion of dissolved oxygen within the rock pores. The predicted copper extractions and Braun’s experimental results differed. Later, the Bartlett model was modified and experimentally verified by Braithwaite for the autoclave leaching of chalcopyrite ore particles with oxygen.

Both the reaction zone and the continuity models furnished a basis for this study, which extends the theory to incorporate intrinsic leaching characteristics of each of the contained sulfide minerals. The model described in this paper uses the continuity equation and includes kinetic expressions for sulfide minerals including chalcopyrite, chalcocite, covellite, and pyrite. The model was derived after interpreting the results of large-scale, long-term laboratory leaching tests. Other leaching results were used to evaluate the accuracy of predicted copper extractions.

Copper Leaching Fundamentals

Copper extraction from waste dumps of low-grade sulfide ores using weakly acidic ferric sulfate solutions proceeds by a series of complicated, not well understood, mechanisms. The mechanics are complicated by the variety of minerals contained in the dump and the changing nature of the dump due to seasonal temperature variations, wet and dry periods, and variations in dump permeability with position and time. Generally, copper leaching requires a supply of both ferric ions to oxidize the sulfides and acid to dissolve copper oxides and to keep dissolved material in solution. Continued copper leaching at a reasonable rate requires ready access to oxygen and certain micro-organisms that oxidize inactive ferrous ion to ferric ion, eliminating the need for continual addition of fresh ferric sulfate. Under favorable leaching conditions, the sulfide minerals react and form soluble cupric sulfate and ferrous sulfate and insoluble sulfur. Leaching ceases when the supply of ferric ions is depleted either through lack of leach liquor circulation or through conditions unfavorable for ferrous oxidation.

Experimental Procedure

Laboratory leach testing was done on a large scale [-6,000 kg (13,227 pounds)] to provide data for samples containing larger fragments than could be leached in bench-scale apparatus. Tests were also conducted using a smaller sample [-200 kg (400 pounds)] for scaleup modeling studies. The large-scale tests were of sufficient size and more readily amenable to scaling up to a commercial venture than bench-scale tests. Leaching results were used to aid in model development and to check the model’s predictions of leaching rates.

Five different copper sulfide samples were leached in the experimental equipment. Ore size distributions, weights, copper analyses, and mineral contents are shown in table 1. The first sample was a monzonite ore containing chalcocite (Cu2S) as the most abundant copper mineral with minor amounts of chalcopyrite (CuFeS2). The second sample was a quartz monzonite ore containing covellite (CuS), bornite (Cu5FeS4), chalcopyrite, and chalcocite. The third sample was also a quartz monzonite ore from the same mine as the second sample and contained the same mineral types. These two samples will be referred to as quartz monzonite I and quartz monzonite II in order to differentiate between them. The fourth sample was a granitic ore containing copper only as chalcopyrite. The fifth sample was a metasediment ore containing chalcocite, native copper (Cu), and trace amounts of cuprite (Cu2O). The weight-percent pyrite for the five samples was 2.7, 6.9, 4.8, 2.9, and 1.2 for the monzonite, quartz monzonite I, quartz monzonite II, granitic, and metasediment samples, respectively. Since the leaching rate of the individual mineral particles is a function of their original size, mineral particle sizes were determined microscopically, and an average size was used for each mineral type in each ore sample. The large ore samples were leached in fiberglass columns 1.37 meters (4.5 feet) in diameter by 3.04 meters (10 feet) high, and the loading capacity was about 8,000 kg (17,637 pounds). In case of the smaller [-200-kg (440-pound)] samples, the ore was leached in 189-liter (50- gallon) stainless steel drums. Solution was evenly distributed over the top of each column by a perforated rotating arm. Figure 1 shows two of the leaching columns.


The samples were leached by downward percolation of dilute ferric sulfate solutions at pH 2.0. Leach solution flow was 20.5 gallons per square foot per day, and solutions were recycled until the copper content was 1 to 2 grams per liter. Then circulation was discontinued, the column was allowed to drain, and the copper was recovered either by cementation on iron or by solvent extraction. After a drying or resting period of about 4 days, leaching was continued with effluent solutions from the copper recovery step. The nearly copper-free pH 2 solutions contained other metallic ions including iron, magnesium, aluminum, calcium, manganese, and zinc. Although natural buffering generally maintained the solution pH at nearly 2, sulfuric acid was added occasionally to prevent ferric iron from precipitating. These conditions typify an ideal portion of the leaching dump where good copper extraction can be achieved. The ionic strength of the leach solutions increased with recycling as other ore components in addition to iron and copper dissolved in the weak acid.


Model Development

Copper sulfide minerals dissolve when oxidized by ferric ions in acidic sulfate solutions. The reaction products include sulfate, elemental sulfur, and both ferrous and cupric ions. Copper dissolution rate primarily is a function of the ferric ion activity at the surface of the unreacted mineral grain. Ferric activity, as a function of position within the ore fragments, depends on the initial ferric ion concentration in the leach solution, ionic strength, temperature, solution chemistry, ferric diffusivity through channels in the porous ore particle, and ferric depletion resulting from reactions with either the copper mineral or with other ore constituents. Ferric depletion is controlled by the reaction kinetic characteristics of the reacting minerals. The hydrogen ion concentration must be maintained to provide favorable conditions for the autotrophic bacteria present, which actively oxidize ferrous iron to ferric maintaining effectively steady-state ferric concentrations in the bulk solution. Thus, many factors are involved in ore leaching, and no simple mathematical model neglecting intrinsic mineral kinetics can adequately describe the process. Attempts to use the mixed kinetic reaction zone model, previously used by Madsen and Braun, failed to describe accurately the ferric sulfate leaching of ores containing a broad array of sulfide minerals.

A model based on the continuity equation for mass transfer was next applied. This model is more complex than the other models, and, because it cannot be solved analytically, numerical evaluation requiring high-speed computer equipment was used. In the model, intrinsic mineral reaction rates were included in a steady-state approximation of the continuity equation to determine the pore diffusion of ferric ions to the mineral surface, thus enabling copper extractions to be calculated. In the calculations, the authors assumed that the ore particles were surrounded by a uniform circulation of leach solution so that bulk solution transport was not a rate-controlling factor. The model did not include temperature effects, calculation of copper extractions from ores containing oxide copper minerals, or permeability of the ore bed. However, by use of the model one can predict an upper limit for copper extraction from a leach dump because the model describes leaching under ideal conditions of adequate air flow and constant temperature.

The Continuity Equation

Leach testing was done at room temperature, and temperature measurements made inside the column showed an essentially constant temperature. With a constant temperature and assuming the leach solution to be uniformly distributed around each ore particle, the modeling problem was reduced to solving the continuity equation for ferric ion distribution within the individual ore fragments.

The continuity equation for ferric ion distribution in an ore particle is


This expression is valid for ideal solutions where the activity coefficients of the diffusing species are equal to unity; that is, ionic activities equal ionic concentrations. However, according to Erdey-Gruz, chemical potential gradients and not concentration gradients should be used when describing diffusional processes. The solutes in a real mixture always diffuse from regions of higher activity towards those of lower activity. The activity effects were included in the effective diffusion coefficient. An empirical activity coefficient was used to relate the effective diffusion coefficient to the diffusion coefficient at infinite dilution according to the equation

Deff = γDeo………………………………………………………………(2)

The ionic strength in solution-filled pores was assumed to be independent of ferric ion reduction to ferrous ion. Therefore, the activity coefficient was assumed to be constant throughout the solution-filled pore volume for a given point in time.

A pore theory proposed by Smith states that the effective diffusion coefficient can be calculated by multiplying the ordinary diffusion coefficient by the porosity and dividing by the tortuosity. This is given by the expression

D eo = εD o/τ……………………………………………………(3)

This equation allows values for the effective diffusivities to be calculated from physical parameters and ordinary diffusion coefficients. Unfortunately, values for the ordinary diffusion coefficients have not been tabulated for the variety of ferric ion complexes present, and tortuosities have not been determined experimentally. Comparison of model predictions with experimental results can be used for determining the effective diffusion coefficient at infinite dilution, Deo.

The operator, ∇², depends on the geometry of the ore particles and is defined by the following equation for the case of spherical ore fragments, as given by


Spherical fragment geometry was used exclusively in this investigation. For spheres, equation 1 becomes


A steady-state approximation of the continuity equation can be applied if the consumption rate of ferric ions, ∑Rk, is much greater than the accumulation of ferric ions, ε. Sohn and Szekely used this approximation when interpreting data involving gaseous reaction in porous solids. The mathematical expression for the approximation is


The assumption of steady state reduces the complex partial differential equation (equation 5) to a simpler ordinary differential expression of the continuity equation


One way to solve the simplified continuity equation (equation 7) using high-speed computer techniques requires finite difference approximations for the derivatives. This was done by the mathematical leaching concept in which the spherical ore particle was divided into concentric shells with a thickness ∆r as described by Braithwaite. There are N numbers of positions located at the shell interfaces. These positions are designated with a subscript j. Position 1 (j=1) corresponds to the center of the fragment, and position N (j=N) corresponds to the outer edge of the fragment.

The relationship of the ferric ion consumption rate in terms of the mineral leaching rates is given by


The term ∑kRjk is the rate of consumption of ferric ions at position j for all k minerals in the fragment. The term (dα/dt)jk (the fractional dissolution rate of mineral k at position j) is a function of ferric ion activity.

Intrinsic Mineral Leaching Kinetics

Knowledge of the leaching kinetics of individual ferric-ion-consuming minerals is necessary to describe the rate term, ∑kRjk, equation 8. Laboratory leaching studies have been conducted by a number of investigators with the data mathematically modeled, and the results were used in this study. The pertinent minerals include chalcocite, secondary covellite formed by the reaction of ferric sulfate and chalcocite, natural covellite, chalcopyrite, and pyrite. Bornite and cuprite were excluded because of the small quantity of these minerals in the ores tested. Native copper leaches very slowly with acid ferric sulfate, therefore in the model calculations it was assumed to leach at the same rate as chalcopyrite.


The dissolution of chalcocite (Cu2S) with ferric sulfate occurs in two stages. The first stage is the conversion of chalcocite to secondary covellite (Cu1.2S), according to the following chemical equation:

5Cu2S + 8Fe+3 → 5Cu1.2S + 4Cu+2 + 8Fe+2………………………………….(9)

Marcantonio described the leaching mechanisms as an electrochemical reaction and showed that at low temperatures the first stage was essentially complete before the second stage started. The leaching kinetics were controlled by the particle size, temperature, fraction of the first-stage leaching completed, and ferric ion activity. Deviations of apparent activities of the ferric ions and activities calculated from thermodynamic and theoretical considerations involving the extended Debye-Huckel theory were accounted for by including a chemisorption term, which effectively reduced the apparent activity of the ferric ions. However, if an apparent activity coefficient is used in place of the chemisorption term in Marcantonio’s model, the following expression describes the leaching kinetics of the first stage of chalcocite leaching:


The determination of values for the apparent activity coefficient, y, is described later in this report.

Secondary Covellite

Marcantonio also observed that after 40 percent of the copper had been leached, the second stage of chalcocite leaching began, which can be represented by the reaction

5Cu1.2S + 12Fe+3 → 6Cu+2 + 5S° + 12Fe+2……………………………………..(11)

The kinetics of this stage followed the relationship


This equation includes Marcantonio’s chemisorption effect in the apparent activity coefficient for ferric ions, γ. The leaching rate for secondary covellite was independent of particle size because secondary covellite is very porous and the transport rate of ferric ions through the porous structure is rapid relative to the leaching rate.

Natural Covellite

Natural covellite is less porous than secondary covellite. Dutrizac and MacDonald determined the leaching rate for polished disks of covellite; when converted to spherical geometry, their rate equation is


Equation 13 is valid for ferric ion concentration less than 5 x 10 -6 mole/cm³, and the term [Fe+3/5 x 10 -6 can be neglected at ferric concentration exceeding 5 x 10 -6 mole/cm³. This equation can be altered to include the activity coefficient, γ, but the value of the rate constant k¹¹¹ would have to be changed to give the same rates of reaction.


Chalcopyrite dissolves in acidic ferric sulfate solutions by the following equation:

CuFeS2 + 4Fe+3 → Cu+2 + 5Fe+2 + 2S°…………………………………………….(14)

Jones and Peters attributed rate control to a surface reaction, whereas others report the reaction rate is limited by transport of iron through a layer of elemental sulfur. Munoz recently reported the opinion that the rate-limiting step is the transport of electrons through the layer of elemental sulfur. The opinion was supported by the independence of the reaction rate with respect to ferric, ferrous, and cupric concentrations plus a high activation energy of 20 kcal/mole (84 kJ/mole). The high activation energy closely agrees with the value calculated for electron transport through elemental sulfur.

The following theoretical rate expression describes electron transport through the reaction layer of elemental sulfur:


Rate data from the initial stage of the reaction, before formation of an appreciable layer of sulfur, showed an inverse first-order dependence on the initial particle diameter, a half-order dependence on ferric concentration, and an activation energy of 8 kcal/mole (34 kJ/mole). The initial reaction was described by equation 17:


Equation 17 was used in this report because the sulfur layer was negligible. Only a small fraction of the chalcopyrite was dissolved, and the presence of sulfur-oxidizing bacteria minimized formation of elemental sulfur. The constant kv in equation 17 was determined by assuming that the leaching rate for chalcopyritic ore equaled the leaching rate of the mineral. This assumption appeared reasonable after 500 days of leaching where the ferric ion activity profiles through even the largest ore particles were computed to be flat. After prolonged leaching, all minerals in the ore particles were exposed to the same ferric ion concentration, and the fractional leaching rate for individual mineral particles equaled the leaching rate of the entire ore sample. The leaching rate of chalcopyrite was assumed to be independent of ferric ion concentrations above 1 x 10 -5 mole/cm³ as shown by Dutrizac.


An expression for the consumption rate of ferric ions by pyrite particles is


The rate of consumption has the units of moles per cubic centimeter per second, as presented in equation 18. This equation is equivalent to the expression derived by Mathews and Robins, except that the units have been changed to be consistent with those used in this study. The combination of terms (WS/V) in equation 18 is the total surface area of pyrite particles per cubic centimeter of ore particle, and all concentrations are expressed as moles per cubic centimeter. This equation was used in the mathematical leaching model, and the pyrite surface area, S, was calculated by microscopically determining both the amount and particle size of the mineral. The hydrogen ion concentration was assumed to be constant throughout the ore particles because of the uniform distribution of acid-generating pyrite.

Comparison of the Leaching Rates of Copper Minerals

Calculated leaching rates are compared in figure 2 for chalcocite, covellite, secondary covellite, and chalcopyrite. The curves for chalcocite and covellite dissolution were determined using equations 10, 12, and 13 for a temperature of 25° C using a leach solution 0.005 molar in ferric sulfate. The leaching rates of covellite and secondary covellite are intermediate between reactive chalcocite and refractory chalcopyrite. Figure 2 also shows that the first-stage leaching of chalcocite is much faster than that of any of the other minerals. However, as the amount of cuprous copper in the solid diminishes, the rate of copper extraction falls to nearly zero, according to


equation 10. In fact, the mathematical expression never allows the fraction of copper reacted to exceed 0.40. An arbitrary value of 0.39 for the fraction reacted was used as a cutoff point where the first-stage leaching of chalcocite ends and the second-stage leaching commences. Natural covellite leaches somewhat slower than chalcocite or secondary covellite. Equation 13 was used to describe the leaching curve for the natural covellite. The slowest leaching mineral was chalcopyrite. Equation 17 was used to calculate the extraction curve for chalcopyrite. Thus, it might be expected that an ore containing an abundance of chalcocite will leach much faster than ore containing chalcopyrite as the main copper mineral.

Apparent Ferric Ion Activity

Ferric ion activity is an important factor in the rate equations. Mathematically predicted copper extractions were consistently higher than experimental extractions when ferric activity coefficients were derived from thermodynamic data and the Debye-Huckel theory. A similar disparity was noted by Marcantonio. The problem was resolved by using an apparent activity coefficient determined empirically by fitting equations 10 and 12 with Marcantonio’s leaching data. The empirical expression for the apparent ferric ion activity coefficient is

ϒ = exp (-1.373 – 239.7 [Fe+3] ½ + 4870 [Fe+3])…………………………………….(19)

Good agreement was achieved between calculated leaching rates, using the empirical activity coefficients, and experimental results. The ferric ion concentration for this comparison ranged from 0.01 to 1.0 molar. Confidence in the validity of equation 19 was reinforced by the similarity between calculated apparent ferric activities and measured activities for trivalent aluminum and indium in sulfate systems. The value for the activity coefficient as calculated in equation 19 was used in the kinetic leaching expressions (equations 10, 12, and 17) and also in the computation of the effective diffusion coefficient (equation 2).

Late Enhancement of the Leaching Rate

Results of copper ore leach testing showed that slightly enhanced leaching rates occurred during the later stages of the tests. This result led to including a factor in the model to describe the phenomenon. However, it must be emphasized that the inclusion of this factor in the model made only small adjustments in the predicted leaching curves. Braun stated that late enhancement results from the generation of cracks and fissures, and the cracking of the leached particles was noted in this study. Braun proposed that a single parameter applied to both the chemically and diffusion controlled processes can be used to describe the late enhancement of the leaching rate. This parameter is a shape factor, ∅v, that varies systematically with the depth of a reaction zone and approaches some limiting degree of surface roughness, The shape factor was related to the fraction of copper leached as


In equation 20 the term λ is the late enhancement constant. As discussed later in this report, the late enhancement term and the effective diffusion coefficient are the only terms evaluated when fitting calculated curves to experimental leaching data. All other parameters were measured in the laboratory prior to leaching an ore sample.

A lower limit of 0.17 was placed on ∅v since smaller values would indicate undue fragmentation of the ore particle. The factor for late enhancement was entered in the leaching model by relating the term ε/τ in equation 3 to the shape factor as follows:

ε/τ = (ε/τ)o/∅v…………………………………………………(21)

Particle-Size Distribution

The ore fragment sizes of the samples leached ranged from about 51 cm (20 inches) to minus 200 mesh. This range of sizes was not narrow enough to be adequately described by a single average particle size. Therefore, the model was extended to include the actual particle-size distribution and the grade of each size fraction.

A normal distribution of broken ore particles can be described by the Schuhmann relationship,

Y = (X/K)Z………………………………………………………………(22)

where Y is the accumulative fraction passing size X (average diameter of an ore particle). A plot of log Y versus log X should result in a straight line with Z as the slope and K as the intercept. Figure 3 is presented as an example of such a plot and shows the particle-size distribution for the minus 33.02-cm (13-inch) plus 1.27-cm (½-inch) monzonite ore.

Using equation 22 and dividing the size distribution into L number of size intervals, a set of L+1 ordered pairs (Xi, Yi) for Xi = K to Xi = Q can be calculated according to the equations


This technique can be used if Q ≠ 0. If Q = 0, then an approximation can be made by setting Q at a positive value close to zero and allowing the


calculations to be carried out. The accumulative size distribution can then be separated into a set of L pairs (rio, Wi) where Wi is the weight fraction of ore having an average radius rio by the following relationships:


Mathematical Leaching Model

All the factors discussed in the preceding text were included in a computer program consisting of a main program and three subroutines. Appendix B contains a description and the program listing for the leaching model. Also included in appendix B are the necessary instructions for feeding data into the computer. The computer language is Fortran V, and many of the computation techniques used were the same as given by Braithwaite.

Factors included in the leaching model were (1) maximum time and time increments, (2) effective diffusivity, (3) stoichiometry, (4) number of shells per particle, (5) intrinsic mineral leaching rate constants, (6) coefficients of polynominals to describe the ferric ion concentration in the bulk solution as a function of time, (7) molar ratio of pyrite to chalcopyrite, (8) late enhancement constant, (9) number of size intervals, (10) constants to calculate the particle size distribution, (11) mineralogy with respect to copper distribution, and (12) grade and density of the ore fragments. All these factors except the effective diffusivity and the late enhancement constant were determined prior to leaching an ore sample. These two factors were determined by correlating calculated copper extractions with experimental results.

Verification of the Leaching Model

The leaching model was evaluated for accuracy by comparing calculated copper extractions with experimentally determined values attained by leaching the monzonite and quartz monzonite I samples using acidic ferric sulfate. Two samples containing different size distributions were leached for each of these ore samples. The size range, weights, and other data concerning these tests are listed in table 1. In the process of determining the best values for the effective diffusivity, Dε0, and the late enhancement term, λ, the model curve was fitted to the experimental data from the minus 15.24-cm (6-inch) plus 1.27-cm (½-inch) samples. The data and modeling results for these columns are included in figures 4 and 5.



The unknown parameters for the granitic ore also were determined by fitting the calculated curve to the experimental data. Since the ore contained only chalcopyrite and since leaching is controlled by the mineral kinetics, the model was fitted by assuming that the effective diffusivity to be similar to that of the monzonite and quartz monzonite I ore (except for differences attributed to porosity differences). Figure 6 contains the calculated extraction along with experimental copper extractions for this ore.

Curve fitting to obtain leaching parameters is useful when the determined values are reasonable and can be successfully used to predict extractions for other experimental conditions. Testing the model’s ability to accurately predict copper extractions for different ore size distributions was done using minus 33-cm (13-inch) plus 1.27-cm (½-inch) monzonite ore and minus 51-cm (20-inch) quartz monzonite I ore. The monzonite ore was screened to remove fine particles, but the quartz monzonite run-of-mine ore was not screened prior to leaching. Both predicted and experimental results are compared in figures 7 and 8. Predicted and experimental results were in good agreement for the monzonite ore; however, the predicted copper extraction for quartz monzonite I ore exceeded the experimental results. The disparity was caused by poor leaching conditions that developed early in leaching the run-of-mine quartz monzonite ore. Ore bed permeability became drastically diminished because salts precipitated in the ore fines, and portions of the ore were denied access to necessary ferric ions. Leach liquor circulation rate decreased from 835 l/m²/day (20.6 gal/ft²/day) to only 33 l/m²/day (0.82 gal/ ft²/day) by the 461st day of leaching. Thus, good leach liquor circulation



enabling uniform lixiviant concentration around each ore particle, a basic requisite for the model, was not achieved using unscreened quartz monzonite I ore.

The model calculations demonstrated that effects due to mineralogy, porosity, ore type, and size distribution are accurately simulated using the proposed model based on the continuity equation. The correlation between the predicted extraction curves with the experimental results was very good.

Leaching tests were made with the metasediment and quartz monzonite II ore samples to determine the feasibility of using data from a small-scale leach test to determine the required parameters for the leaching model. Both small-scale [-200-kg (440-pound)] and large-scale [-6,000-kg (13,227-pound)] leaching tests were made with each ore sample, and the experimental data from the smaller test were used to establish the effective diffusivity and late enhancement constant for the model. The established parameters were then used to predict the copper extraction from the larger sample, and the prediction was compared with the experimental data. Figure 9 shows the experimental data and the calculated copper extractions for the small-scale test with the meta-sediment ore, and figure 10 shows the predicted extractions and experimental extractions for the large-scale test. Since the predictions in figure 10 near the end of the test period were lower than the experimental copper extractions, it was concluded that the late enhancement constant obtained with the small-scale test was too low. This discrepancy occurs because the smaller sized particles are affected less by the late enhancement term (equation 20)


owing to a smaller volume and fewer imperfections. Also, the large number of particles in the large tests make it more representative of the ore sample than the small test. However, in spite of these limitations, the scaleup predictions were always within 7 percent of the observed copper extraction.

Figure 11 shows the experimental data and the calculated copper extractions for the small-scale test with the quartz monzonite II ore, and figure 12 shows the predicted extractions and experimental extraction for the large-scale test. As with the metasediment ore, the prediction of copper extraction from the quartz monzonite II ore was low near the end of the test period. The same reasons for the discrepancy between predicted and observed results are also valid for this case.

The scaleup modeling studies have shown that the model is valid in predicting copper extractions. Application of the model would make laboratory leaching studies less costly if only the small-scale tests are needed for obtaining the leaching parameters required.



Comparison of Leaching Parameters for Different Ores

Long-term ore leaching data were used to calculate ore tortuosities using equation 3 relating tortuosity to ore porosity and to both the ordinary and effective diffusivities. The effective ferric diffusivities were determined by curve-fitting techniques using long-term leaching data, and the ordinary diffusivity was estimated by an equation given by Perry. Ore porosity was measured by water absorption (ASTM C20-70). Values determined for the five copper ores are shown in table 2. That there was little difference in the tortuosities indicated that there was similar pore structure in each of the ores. After determining ore porosity, effective diffusivities may be easily calculated. This information coupled with the mineral content and size distribution of the ore can be used to predict copper extraction curves for different bulk solution ferric ion concentrations without the need for long-term leach testing.


All ores except the metasediment ore required very little adjustment due to the late enhancement of leaching. With the late enhancement constant, λ, equal to 0.01, the calculated leaching curves are only slightly different from curves generated without adjustment. The value of 0.01 for λ appears to be a suitable average value to use with copper porphyry ores. However, with a sedimentary type ore such as the metasediment ore sample, a larger value must be used to obtain accurate predictions of leaching behavior. The metasediment ore required a late enhancement constant of 0.06 and this high value was justified owing to the weathering noted when the leach residue was examined. The metasediment ore broke down much more extensively than did the other ores tested. A microscopic evaluation of the host rock before leaching could identify the tendency for ores to weather during leaching. Typically such a microscopic examination could show the presence of slate-type host rock formations and of sandstone-type rocks that break down easily when exposed to an acidic solution.

Copper Extraction and Ferric Ion Concentration Profiles

In addition to predicting the composite copper extraction for a large ore sample, the model provides computations that describe the copper extraction and ferric ion concentration profile within individual ore fragments. Analyses of these computed profiles are helpful in determining the mechanisms controlling the leaching kinetics. To illustrate both mechanisms, that is, the intrinsic mineral leaching kinetics and the rate of diffusion of lixiviant through the host rock, copper extraction and ferric ion concentration profiles for fragments of monzonite and granitic ores are depicted graphically in figures 13 and 14, respectively.

In figure 13, the copper extraction profiles for a monzonite ore fragment (which originally contained both chalcocite and chalcopyrite) show that copper extraction abruptly decreases from the outside to the center of ore fragments leached for 72 and 245 days. The steep slope of these copper extraction profiles is due to the rapid leaching of chalcocite and secondary covellite in this zone. The reaction zone, 5, moves inwardly through the ore fragment with time, and the thickness of the zone remains essentially constant. The region on the outward side of the zone is devoid of chalcocite and covellite. In this region chalcopyrite is the only remaining copper mineral, and the slight decrease in the extraction profile reflects the slow leaching character of chalcopyrite. The region on the inward side of the zone is unleached and contains both chalcocite and chalcopyrite. The ferric ion concentration profiles in figure 13 show that the ferric ion concentration is near zero from the inward side of zone 5 to the center of the ore fragment. This explains why there was no leaching in this region. However, after 520 days of leaching, the reaction zone for chalcocite and covellite leaching has advanced to the center of the ore fragment and has disappeared. At this time the only remaining copper mineral in the ore fragment is chalcopyrite, which continues to leach at a slow rate. The shape of the copper extraction profiles indicates the mechanism that is controlling the leaching rate. With steeply dipping copper extraction profiles, as with the 72- and 245-day profiles, ferric ions are rapidly consumed and the leaching kinetics are controlled by the diffusion of ferric ions to the leaching zone, δ. However, when steep copper


extraction profiles are not found, such as for the 520-day profile, the leaching kinetics are controlled by the intrinsic leaching rate at the surface of the individual copper minerals. Thus, the mechanism controlling the leaching rate of the 1.54-cm (½-inch) monzonite ore fragment changed between the 245th and 520th day of leaching.

Figure 14 shows the copper extraction and ferric ion concentration profiles for an ore fragment from the granitic ore sample. Since this ore contains copper as chalcopyrite only, all the copper extraction profiles are nearly horizontal. As mentioned previously, this type of profile indicates that the leaching kinetics are controlled by the chemical reactions at the surface of the mineral particles. The decrease of ferric ion concentration from the outside to the center of the ore fragment is due chiefly to the reaction of ferric ions with the pyrite contained in the ore. However, by 245 days the ferric ion concentration is high enough in every position in the fragment to give the maximum leaching rate for the chalcopyrite mineral. Hence, the chalcopyrite minerals at the center of the ore fragment leach as fast as the chalcopyrite minerals on the outside surface of the ore fragment.

Ores Containing a Variety of Copper Sulfide Minerals

The mathematical model used in this study adequately predicts the rate of copper leaching from a large ore sample using acidic-ferric sulfate solution. The accuracy of the prediction was shown to depend on a basic premise—good permeability in the ore bed permitting uniform ferric ion concentration around the ore fragments.

The model is useful in predicting copper recovery rates for large test samples when using leaching parameters obtained from much smaller scale tests. This application of the model would make laboratory leaching studies less costly and more efficient.

The analysis of copper extraction profiles computed by the leaching model is useful in determining whether diffusion of lixiviant through the host rock or chemical leaching at the mineral surfaces is controlling the leaching kinetics.

The ore particle tortuosities were calculated through computer analysis and were similar for each of the four different porphyry copper sulfide ores tested. If the average tortuosity for these four samples is used in the leaching model, then a first approximation of copper extraction can be predicted for laboratory leaching experiments where bulk solution transport is not rate controlling. Also, an approximation to the upper limit of copper extraction from copper sulfide leaching dumps, where bulk solution transport may be rate controlling, can be obtained by using the model predictions.