Redox-Controlled Multiple-Species Reactive Chemical Transport 2. Verification and Application

A redox-controlled multiple-species, multidimensional reactive chemical transport model, DYNAMIX, has been presented in part 1 of this paper. In part 2 the model is checked against two other independently developed reactive chemical transport codes. Simulated results from DYNAMIX agree reasonably well with those obtained from the other two models. Two approaches for simulating redox reactions, an external approach (based on hypothetical electron activity) and an effective internal approach (based on conservation of electrons), have been examined. The results show that in the external approach the simulated redox front is smeared out by diffusion and dispersion, whereas in the effective internal approach the redox front is sharp and is controlled by the redox reactions. The evolution of solid phases from the two approaches give markedly different mineral assemblages. It is suggested that the external approach is of relevance in industrial processes such as electrometallurgy, in which mineral dissolution is driven by externally supplied electric power. The effective internal approach is applicable to hydrogeochemical systems (groundwater contamination, diagenesis, ore formation, rock weathering, and soil genesis) in which redox potential is dictated by the states of the redox species. The results also suggest that in the presence of precipitation (and consequent retardation of the concentration front), spreading due to hydrodynamic dispersion may be significantly inhibited. To demonstrate the applicability of DYNAMIX for realistic field problems, two field cases are simulated by this model. In the first case, the supergene enrichment of copper at Butte, Montana, the simulated mineral assemblages and their distributions agree reasonably with the mineral assemblages and the ore grade observed in the field. In the second case, involving two-dimensional contaminant transport in a hypothetic aquifer, the simulated results suggest that DYNAMIX is capable of handling a realistic multi-dimensional field problem with several hundred grid blocks and time scales of practical interest.

[1]  J. Ague,et al.  Geochemical modeling of steady state fluid flow and chemical reaction during supergene enrichment of porphyry copper deposits , 1989 .

[2]  M. J. Baedecker,et al.  Modern marine sediments as a natural analog to the chemically stressed environment of a landfill , 1979 .

[3]  Grant Garven,et al.  Theoretical analysis of the role of groundwater flow in the genesis of stratabound ore deposits; 2, Quantitative results , 1984 .

[4]  C. Carnahan Simulation of uranium transport with variable temperature and oxidation potential: The computer program THCC (Thermo-Hydro-Chemical Coupling) , 1986 .

[5]  Larry W. Lake,et al.  Precipitation and dissolution of solids attending flow through porous media , 1984 .

[6]  Robert L. Street,et al.  A Groundwater Mass Transport and Equilibrium Chemistry Model for Multicomponent Systems , 1985 .

[7]  D. Hydes,et al.  Oxidation Fronts in Pelagic Sediments: Diagenetic Formation of Metal-Rich Layers , 1986, Science.

[8]  F. Ransome The Copper Deposits of Ray and Miami, Arizona , 1919 .

[9]  G. Brimhall Early fracture-controlled disseminated mineralization at Butte, Montana , 1977 .

[10]  T. N. Narasimhan,et al.  Groundwater contamination from an inactive uranium mill tailings pile: 2. Application of a dynamic mixing model , 1984 .

[11]  T. J. Wolery,et al.  Calculation of chemical equilibrium between aqueous solution and minerals: the EQ3/6 software package. [In FORTRAN extended 4. 6 for CDC6600 and 7600] , 1979 .

[12]  Y. Kwong A new look at the Afton copper mine in the light of mineral distributions, host rock geochemistry and irreversible mineral-solution interactions , 1981 .

[13]  G. Garven,et al.  Theoretical analysis of the role of groundwater flow in the genesis of stratabound ore deposits; 1, Mathematical and numerical model , 1984 .

[14]  G. Brimhall Lithologic determination of mass transfer mechanisms of multiple-stage porphyry copper mineralization at Butte, Montana; vein formation by hypogene leaching and enrichment of potassium-silicate protore , 1979 .

[15]  C. Miller,et al.  Simulation of Solute Transport in a Chemically Reactive Heterogeneous System' Model Development and Application , 1983 .

[16]  S. Titley Geological characteristics and environment of some porphyry copper occurrences in the southwestern Pacific , 1975 .

[17]  H. Robberecht,et al.  Selenium in environmental waters: Determination, speciation and concentration levels. , 1982, Talanta.

[18]  David L. Parkhurst,et al.  Phreeqe--A Computer Program for Geochemical Calculations , 1980 .

[19]  C. Meyer An early potassic type of wall-rock alteration at butte, Montana , 1965 .

[20]  C. Alpers,et al.  Analysis of supergene ore-forming processes and ground-water solute transport using mass balance principles , 1985 .

[21]  J. D. Jabro,et al.  Gaseous diffusion equations for porous materials , 1982 .

[22]  B. Wood,et al.  Rates of Hydrothermal Reactions , 1983, Science.

[23]  C. Palmer,et al.  Geochemical evolution of groundwater in sequences of sedimentary rocks , 1984 .

[24]  Randall Brown Selenium Threatens Irrigators, Wildlife , 1987 .

[25]  J. Rimstidt,et al.  Rates of reaction of pyrite and marcasite with ferric iron at pH 2 , 1984 .

[26]  T. N. Narasimhan,et al.  Redox‐controlled multiple‐species reactive chemical transport: 1. Model development , 1989 .