A Mathematical Model For The In-Situ Leaching Of Primary Copper Ore ? Introduction

During the last decade, the feasibility of applying the in situ leaching technology to the recovery of copper from deeply buried deposits containing chalcopyrite has received much attention.(1-4) The major incentive for using this technology is that under certain circumstances, copper in deep underground deposits, which cannot be mined economically by conventional means, can be recovered economically. Basically, this technology involves rubblization of the deeply buried copper ore body followed by in-place leaching of the rubblized copper ores. Rubblization of the ore can be achieved by some suitable method such as the use of block caving mining. The rubblized chimney is then filled with water until the original water level is reached. Leaching occurs when the eater saturated with oxygen is introduced into the chimney- and flows through the fractured rock masses and within pores of rocks. As the dissolved oxygen oxidizes the primary sulfides, heat is generated. Because of low conductivity of the surrounding rocks, the generated heat is able to raise the temperature of the chimney. The increase in the temperature of the chimney and the relatively high content of dissolved oxygen in the leach solution under high hydrostatic pressure will result in an enhanced dissolution of the primary sulfide minerals. To ensure a successful application of the in situ leaching technology, a complete understanding of the effects of the ore particle size distribution, ore grade mineralogy, temperature, pressure, lixiviant concentration, pH, etc. on the rate of copper recovery is necessary so that copper can be recovered efficiently and the amount of copper extracted can be predicted. These can be achieved through mathematical modeling of the in situ leaching process and conducting laboratory tests. In this paper we will present a mathematical model that simulates the in situ leaching process as depicted in Figure 1. To model this process properly, not only the 'kinetic and thermodynamic aspects of leaching chemistry must be known, but the effects of heat generation, the physical aspect of fluid flow and other related multicomponent transport phenomena must be considered. Furthermore, since most factors involved in the leaching process change with time, a leaching system must be described using unsteady state equations.