Hello guys
As a GWB user to study the fluid rock interaction in hydrothermal ore deposit, one of the most important initial work is to obtain the composition of fluids. Generally, we use LA-ICP-MS analysis to get the chemical composition of ore fluids. But the data have to be processed with SpecE8 before used as the input file for React. However, the processing by SpecE8 is seldom detailed shown in the paper, confusing starters like me a lot.
Here I’d like to show a typical example, and wish to get some help.
According to microthermometry and LA-ICPMS analysis of fluid inclusions, we got the following conditions:
T=145 °C, CH4=300 ppm, Cl-=4.81 mol/kg, Ca2+ =0.5 mol/kg,Mg2+ =0.12 mol/kg, K+=0.1 mol/kg, Na+=3.5 mol/kg
By studying phase diagram, we got other key parameters:
pH=4.4, log fO2=-57, H2S= 2*10-3 mol/kg
Since the deposit contain quartz, feldspar, pyrite, chalcopyrite, sphalerite and galena, we use these minerals to buffer the contents of SiO2(aq), Al3+, Fe2+, Cu2+, Zn2+ and Pb2+ in fluids. In addition, I also use dolomite to balance Ca2+ in the fluids.
However, what I got is: Newton-Raphson did not converge after 999 iterations…
So, how can I solve such a problem?
Pb-Zn deposit-initial fluid.sp8