Numerical simulation of a base metal deposit related to a fossil geothermal system

Fossil and active geothermal systems that produce ore deposits are sites of complex physicochemical processes and a favorable combination of factors related to the amount of metal-bearing fluid that flows through the system, ore fluid metal concentrations, depositional efficiency, and the duration of ore deposition. Of all these factors, the length of the mineralizing event is one of the least understood aspects of ore genesis.

We used fluid inclusion data, chemical compositions of base metal sulfides, and fluid flow rates to constrain a reactive-transport model of a fossil geothermal system - the Patricia Zn-Pb-Ag deposit in northern Chile. The Patricia deposit consists of quartz and base metal sulfide veins of hydrothermal origin with structural control, hosted in a volcanic succession with intense propylitic alteration. The fluid inclusions are liquid-rich, with homogenization temperatures ranging from 250 to 150 °C and salinities between 22 and 1 wt% NaCl equiv., with an early fluid mixing trend and no evidence of boiling in the system. Sulfide mineralogy indicates intermediate sulfidation conditions.

To identify the most relevant geochemical and transport parameters controlling the formation of this fossil geothermal system >1000 simulations were performed using the reactive-transport code TOUGHREACT. The paragenesis of the deposit is mimicked by a model of successive stages of fluid circulation consistent with the observed mineral assemblage distribution, the fluid inclusion data, and the estimated resources in the deposit.

The entire geothermal activity of the system was modeled considering 10,000 years of fluid-rock interaction, with periods of circulation of metal-barren fluids followed by metal-rich fluids driving the ore formation. In the initial model, base metal solubility with predominant chloride complexing suggests that the most efficient ore-forming mechanism for the Patricia deposit was the result of the interaction of two different fluids, one fluid transporting metals and another fluid transporting reduced sulfur, mixing in a rock volume of high permeability. Mass balance estimations with this model give a period of 3500 to 5000 years for the ore stage duration in which all the ore resources of the Patricia deposit could have been precipitated by fluid mixing.

In a second model, the previous estimates for the duration of the main ore stage were used to simulate the fluid-rock interaction during the ore stage for 3500 years. The results indicated the importance of the permeability of the host rock enhanced by fractures to concentrate the volume of the mineralization and the role of the hydrothermal alteration assemblage in controlling the circulating fluid acidity. A higher efficiency in forming sulfide minerals appears to coincide with pH values ranging from 5.1 to 5.3.

The results of both models are validated by replicating the system evolution, reproducing the same mineral alteration assemblage, the expected base metal resource distribution, and similar amounts of ore resources to those of the Patricia deposit: a total of 52,602 tons of lead and 157,731 tons of zinc. The models indicate that the hydrothermal event might be two times longer than the ore stage.