# phreeqc.dat#Halite = 5.0% = 76.74 mol/kgw#Quartz = 50 % = 746.42 mol/kgw#Illite = 20.0 % = 46.73 mol/kgw#Dolomite = 5.0 % = 24.32 mol/kgw#Anhydrite = 15.0 % = 131.77 mol/kgwUSER_PRINT5 PRINT PHASE_FORMULA$("Halite")10 Halite = 5/100*1000 / GFW(PHASE_FORMULA$("Halite"))20 Quartz = 50/100*1000 / GFW(PHASE_FORMULA$("Quartz"))30 Illite = 20/100*1000 / GFW(PHASE_FORMULA$("Illite"))40 Dolomite = 5/100*1000 / GFW(PHASE_FORMULA$("Dolomite"))50 Anhydrite = 15/100*1000 / GFW(PHASE_FORMULA$("Anhydrite"))60 print "Halite: ", Halite, " moles"70 print "Quartz: ", Quartz , " moles"80 print "Illite: ", Illite , " moles"90 print "Dolomite: ", Dolomite , " moles"100 print "Anhydrite: ", Anhydrite , " moles"110 H2O = Halite / 76.74120 print "H2O, kg, from Halite: ", H2O130 H2O = Quartz / 746.42 140 print "H2O, kg, from Quartz: ", H2OSOLUTION
Mineral mol/L water Density, g/cm^3 Halite 43.9 2.16 Quartz 426.9 2.65 Illite 26.7 2.71 Dolomite 13.9 2.86 Anhydrite 56.5 2.95
#Weight percent#Halite = 5.0% = 76.74 mol/kgw#Quartz = 50 % = 746.42 mol/kgw#Illite = 20.0 % = 46.73 mol/kgw#Dolomite = 5.0 % = 24.32 mol/kgw#Anhydrite = 15.0 % = 131.77 mol/kgwUSER_PRINT10 DATA "Halite", "Quartz", "Illite", "Dolomite", "Anhydrite"20 DATA 5, 50, 20, 5, 1530 DIM min$(5), pct(5), vol_per_kg_rock(5), mol_per_kg_rock(5), \ mol_per_L_water(5), min_density(5)40 RESTORE 1050 FOR i = 1 to 5 60 READ min$(i)70 NEXT i80 RESTORE 2090 FOR i = 1 to 5100 READ pct(i)110 NEXT i120 rock_volume = 19 # L130 water_volume = 1 # L200 REM calculate mol/kg rock, L/kg rock, rock density210 FOR i = 1 to 5 220 mol_per_kg_rock(i) = pct(i)/100*1000 / GFW(PHASE_FORMULA$(min$(i)))230 NEXT i240 rock_density = 2.7 # kg/cm^3250 rock_mass = rock_volume * rock_density300 REM calculate mol / L water310 FOR i = 1 to 5 320 mol_per_L_water(i) = mol_per_kg_rock(i) * rock_mass / water_volume330 min_density(i) = GFW(PHASE_FORMULA$(min$(i))) / PHASE_VM(min$(i)) 340 NEXT i400 REM Print410 PRINT "Rock density: ", rock_density420 PRINT "Rock mass: ", rock_mass430 PRINT440 PRINT "Mineral mol/L water Density, g/cm^3"450 FOR i = 1 to 5 460 PRINT PAD(min$(i), 10), STR_F$(mol_per_L_water(i), 10, 1), STR_F$(min_density(i), 15, 2)470 NEXT iSOLUTION