Magazynowanie ciepła i chłodu w warstwach wodonośnych to bezpośrednia interakcja ze środowiskiem gruntowo-wodnym, gdyż to warstwa wodonośna stanowi tu przestrzeń magazynową, a nośnikiem energii jest woda podziemna. Ze względu na potrzebę ochrony zasobów wód podziemnych, każdy system korzystający z nich bezpośrednio powinien być systemem zrównoważonym - woda pobrana z warstwy wodonośnej powinna w całości zostać zatłoczona do tej samej warstwy w ramach poboru zwrotnego. Pobór, zwłaszcza znaczący, wód do celów innych niż wymiana ciepła, jeżeli jest konieczny, najlepiej żeby odbywał się poza obrębem magazynów, gdzie woda występuje w temperaturze naturalnej i równowadze chemicznej z materiałem budującym warstwę wodonośną. Z tego też powodu systemy ATES nie powinny być lokalizowane w strefach dopływu wód podziemnych do istniejących lub planowanych ujęć wód pitnych.
Różnice temperatury mogą powodować zmiany w składzie chemicznym wód, ponieważ temperatura odgrywa bardzo ważną rolę w rozpuszczalności minerałów, kinetyce reakcji chemicznych, utlenianiu materii organicznej, procesach redoks oraz sorpcji-desorpcji anionów i kationów. Wprawdzie zmiany temperatury stosowane w niskotemperaturowych systemach ATES, nie przekraczające na ogół 15°C, wydają się nie mieć prawie żadnego wpływu na zawartość i formę głównych składników chemicznych wód, jednak wrażliwość warunków redoks na niewielkie zmiany temperatury, a zwłaszcza zwiększona mobilność niektórych mikroskładników, np. arsenu, obserwowana w eksperymentach laboratoryjnych opisywanych w literaturze, wskazują na konieczność dalszych badań i ciągłego monitorowania wód w rejonach funkcjonowania systemów ATES. Ze względu na zupełny brak danych dotyczących funkcjonowania tego typu systemów w Polsce, niezbędne jest zaprojektowanie instalacji pilotowych w lokalizacjach o typowych dla kraju warunkach hydrogeologicznych i przeprowadzenie przynajmniej kilkuletnich obserwacji pracy takich systemów. Wydaje się jednak sensowne i potwierdzają to już obecnie dostępne publikacje, że podczas pracy systemu ATES należy w miarę możliwości unikać kontaktu użytkowanych wód z powietrzem atmosferycznych, co zapobiegnie zmianom potencjału utleniająco-redukcyjnego, które zwłaszcza w powiązaniu ze zmianami temperatury, mogą prowadzić do wytrącania się związków słabo rozpuszczalnych i kolmatacji instalacji.
Kolmatacja w systemach ATES znacząco wpływa na ich wydajność i koszty eksploatacji, może prowadzić do całkowitego zablokowania niezbędnego przepływu wód na wymiennikach ciepła lub/i w strefie filtrowej zwłaszcza otworów zatłaczających. Zrozumienie mechanizmów biologicznych, chemicznych i fizycznych tego procesu jest kluczowe dla opracowania skutecznych strategii zapobiegawczych i konserwacyjnych. Dalsze badania nad interakcjami pomiędzy mikroorganizmami, warunkami środowiskowymi oraz procesami biochemicznymi mogą przyczynić się do bardziej precyzyjnego przewidywania i ograniczania kolmatacji. Zastosowanie zaawansowanych technologii monitorowania i modyfikacji systemów ATES (zarówno nisko i wysoko temperaturowych) pozwoli nie tylko zwiększyć ich efektywność, ale także przyczyni się do upowszechnienia tego typu magazynowania energii cieplnej, co z kolei pozwoli na efektywne wykorzystanie zrównoważonych źródeł energii cieplnej w ciepłownictwie sieciowym oraz redukcję zanieczyszczeń pochodzących ze spalania paliw kopalnych i "alternatywnych" w ciepłownictwie indywidualnym.
Czytaj dalej
Modelowanie wpływu zmian temperatury w ATES na chemizm wód podziemnych - analizy archiwalne
Formy występowania poszczególnych składników chemicznych w badanych wodach podziemnych
Poszczególne składniki rozpuszczone w wodach mogą występować w różnych postaciach; jako uwodnione wolne jony, w formie związków kompleksowych, par jonowych, kompleksowych związków organicznych czy cząsteczek obojętnych (np. rozpuszczony tlen). Zależeć to będzie od rozpuszczalności danej substancji w wodzie, od jej zdolności do tworzenia kompleksów i par jonowych, od pH i Eh wody, jej temperatury, siły jonowej oraz składu chemicznego.
Wykorzystując wyniki modelowania hydrogeochemicznego, wykonanego za pomocą programu Wateq4f, poniżej przedstawiono podstawowe formy występowania głównych składników chemicznych w wytypowanych środowiskach wód podziemnych wykonanego dla wybranych punktów obserwowanych w ramach zadań państwowej służby geologicznej (tabele 1-7, lokalizacja punktów - figura 1).
Tabela 1. Formy występowania głównych składników chemicznych w wodach w punkcie 991
Forma składnika | % zawartość | |||
5.0oC | 10.8oC | 15oC | ||
HCO3 | -1 | 86.23 | 87.41 | 88.03 |
H2CO3 aq | 0 | 12.01 | 10.59 | 9.79 |
Fe | 2 | 2.81 | 2.15 | 1.76 |
Fe(OH)2 | 1 | 48.17 | 41.77 | 37.31 |
Fe(OH)3 | 0 | 47.16 | 54.32 | 59.16 |
Mn | 2 | 63.22 | 61.92 | 60.99 |
MnCO3 aq | 0 | 11.83 | 13.74 | 15.09 |
MnHCO3 | 1 | 24.58 | 23.94 | 23.49 |
Zn | 2 | 47.90 | 45.54 | 43.89 |
ZnCO3 | 0 | 22.51 | 25.38 | 27.28 |
Zn(CO3)2 | -2 | 2.49 | 3.34 | 4.02 |
ZnHCO3 | 1 | 26.31 | 24.88 | 23.88 |
Forma składnika | % zawartość | |||
5.0oC | 9.6oC | 15.0oC | ||
Ba | 2 | 92.45 | 92.25 | 91.97 |
BaHCO3 | 1 | 1.83 | 2.10 | 2.47 |
BaSO4 aq | 0 | 5.67 | 5.59 | 5.48 |
Ca | 2 | 94.91 | 94.40 | 93.83 |
CaHCO3 | 1 | 2.96 | 3.38 | 3.84 |
CaSO4 aq | 0 | 1.92 | 1.98 | 2.05 |
HCO3 | -1 | 83.09 | 84.30 | 85.36 |
H2CO3 aq | 0 | 15.20 | 13.77 | 12.46 |
Cu | 2 | 6.56 | 6.07 | 5.57 |
CuCO3 aq | 0 | 55.60 | 58.96 | 62.37 |
CuHCO3 | 1 | 12.72 | 11.72 | 10.69 |
Cu(OH)2 | 0 | 23.95 | 22.11 | 20.23 |
Fe | 2 | 4.52 | 3.73 | 2.96 |
FeHCO3 | 1 | 1.75 | 1.44 | 1.13 |
Fe(OH)2 | 1 | 53.57 | 48.88 | 43.28 |
Fe(OH)3 | 0 | 39.61 | 45.31 | 51.85 |
Mg | 2 | 94.15 | 93.95 | 93.66 |
MgHCO3 | 1 | 4.16 | 4.14 | 4.14 |
MgSO4 aq | 0 | 1.59 | 1.79 | 2.04 |
Mn | 2 | 67.30 | 66.48 | 65.53 |
MnCO3 aq | 0 | 8.43 | 9.55 | 10.86 |
MnHCO3 | 1 | 23.20 | 22.82 | 22.37 |
Forma składnika |
% zawartość | |||
5.0 oC | 10oC | 16.7oC | ||
Ba | 2 | 85.45 | 85.40 | 85.29 |
BaHCO3 | 1 | 1.46 | 1.70 | 2.08 |
BaSO4 aq | 0 | 13.03 | 12.82 | 12.52 |
Ca | 2 | 92.54 | 91.96 | 91.22 |
CaHCO3 | 1 | 2.50 | 2.89 | 3.35 |
CaSO4 aq | 0 | 4.67 | 4.81 | 5.00 |
HCO3 | -1 | 88.35 | 89.07 | 89.69 |
H2CO3 aq | 0 | 9.42 | 8.40 | 7.40 |
Cu | 2 | 3.94 | 3.66 | 3.33 |
CuCO3 aq | 0 | 48.43 | 52.11 | 56.35 |
CuHCO3 | 1 | 6.60 | 6.09 | 5.51 |
Cu(OH)2 | 0 | 39.73 | 36.82 | 33.45 |
Fe(OH)2 | 1 | 43.26 | 37.58 | 30.64 |
Fe(OH)3 | 0 | 53.73 | 59.67 | 66.69 |
Fe(OH)4 | -1 | 0.64 | 0.89 | 1.33 |
Mg | 2 | 92.38 | 91.86 | 91.06 |
MgHCO3 | 1 | 3.54 | 3.51 | 3.49 |
MgSO4 aq | 0 | 3.88 | 4.39 | 5.14 |
Mn | 2 | 65.80 | 64.50 | 62.80 |
MnCO3 aq | 0 | 11.98 | 13.60 | 15.71 |
MnHCO3 | 1 | 19.61 | 19.12 | 18.47 |
MnSO4 aq | 0 | 2.36 | 2.54 | 2.78 |
Zn | 2 | 49.85 | 47.63 | 44.78 |
ZnCO3 | 0 | 22.79 | 25.23 | 28.15 |
Zn(CO3)2 | -2 | 2.75 | 3.54 | 4.71 |
ZnHCO3 | 1 | 20.99 | 19.94 | 18.60 |
ZnSO4 aq | 0 | 3.01 | 2.96 | 2.88 |
Tabela 4. Formy występowania głównych składników chemicznych w wodach w punkcie 2149
Forma składnika | % zawartość | |||
5.0oC | 12oC | 15oC | ||
Ba | 2 | 94.25 | 94.05 | 93.95 |
BaHCO3 | 1 | 0.99 | 1.22 | 1.33 |
BaSO4 aq | 0 | 4.61 | 4.52 | 4.48 |
Ca | 2 | 96.28 | 95.71 | 95.46 |
CaCO3 aq | 0 | 0.58 | 0.73 | 0.81 |
CaHCO3 | 1 | 1.59 | 1.93 | 2.06 |
CaSO4 aq | 0 | 1.55 | 1.63 | 1.66 |
CO3 | -2 | 0.29 | 0.35 | 0.38 |
HCO3 | -1 | 94.77 | 94.90 | 94.90 |
H2CO3 aq | 0 | 3.37 | 2.86 | 2.69 |
CuCO3 aq | 0 | 17.72 | 20.87 | 22.17 |
Cu(OH)2 | 0 | 80.09 | 76.94 | 75.64 |
Mg | 2 | 96.17 | 95.81 | 95.64 |
MgCO3 aq | 0 | 0.31 | 0.42 | 0.47 |
MgHCO3 | 1 | 2.24 | 2.22 | 2.22 |
MgSO4 aq | 0 | 1.28 | 1.54 | 1.66 |
Mn | 2 | 64.43 | 61.29 | 60.01 |
MnCO3 aq | 0 | 22.97 | 26.71 | 28.22 |
MnHCO3 | 1 | 11.79 | 11.13 | 10.85 |
PbCO3 aq | 0 | 93.17 | 93.57 | 93.66 |
Pb(CO3)2 | -2 | 2.03 | 2.50 | 2.71 |
PbHCO3 | 1 | 1.95 | 1.59 | 1.47 |
PbOH | 1 | 1.55 | 1.27 | 1.18 |
H4SiO4aq | 0 | 99.39 | 99.17 | 99.07 |
Sr | 2 | 96.85 | 96.28 | 96.00 |
SrHCO3 | 1 | 1.55 | 1.94 | 2.13 |
SrSO4 aq | 0 | 1.45 | 1.55 | 1.59 |
Zn | 2 | 41.33 | 36.87 | 35.12 |
ZnCO3 | 0 | 37.01 | 40.36 | 41.49 |
Zn(CO3)2 | -2 | 6.85 | 9.19 | 10.21 |
ZnHCO3 | 1 | 10.68 | 9.45 | 8.97 |
Zn(OH)2 | 0 | 2.71 | 2.41 | 2.30 |
Tabela 5. Formy występowania głównych składników chemicznych w wodach w punkcie 2346
Forma składnika | % zawartość | |||
5.0oC | 9.8oC | 15oC | ||
H3BO3 aq | 0 | 91.55 | 90.52 | 89.41 |
H2BO3 | -1 | 8.45 | 9.48 | 10.59 |
Ba | 2 | 89.51 | 89.42 | 89.31 |
BaSO4 aq | 0 | 9.60 | 9.52 | 9.42 |
Ca | 2 | 94.54 | 94.08 | 93.55 |
CaCO3 aq | 0 | 1.07 | 1.25 | 1.49 |
CaHCO3 | 1 | 1.05 | 1.20 | 1.35 |
CaSO4 aq | 0 | 3.33 | 3.45 | 3.59 |
HCO3 | -1 | 96.71 | 96.55 | 96.32 |
H2CO3 aq | 0 | 1.26 | 1.12 | 1.00 |
CuCO3 aq | 0 | 4.96 | 5.66 | 6.43 |
Cu(OH)2 | 0 | 94.55 | 93.82 | 93.02 |
Fe(OH)2 | 1 | 7.62 | 6.00 | 4.63 |
Fe(OH)3 | 0 | 84.12 | 83.82 | 82.77 |
Fe(OH)4 | -1 | 8.16 | 10.12 | 12.56 |
Mg | 2 | 95.08 | 94.57 | 93.96 |
MgHCO3 | 1 | 1.49 | 1.47 | 1.46 |
MgSO4 aq | 0 | 2.77 | 3.15 | 3.59 |
Mn | 2 | 54.92 | 51.99 | 49.11 |
MnCO3 aq | 0 | 36.76 | 40.03 | 43.22 |
MnHCO3 | 1 | 6.78 | 6.37 | 5.97 |
MnSO4 aq | 0 | 1.39 | 1.45 | 1.51 |
H4SiO4aq | 0 | 98.35 | 97.97 | 97.50 |
H3SiO4 | -1 | 1.65 | 2.03 | 2.50 |
Sr | 2 | 95.56 | 95.13 | 94.62 |
SrHCO3 | 1 | 1.03 | 1.20 | 1.41 |
SrSO4 aq | 0 | 3.12 | 3.29 | 3.47 |
Tabela 6. Formy występowania głównych składników chemicznych w wodach w punkcie 1088
Forma składnika | % zawartość | |||
5.0oC | 9.0oC | 15oC | ||
Ba | 2 | 93.96 | 93.76 | 93.42 |
BaHCO3 | 1 | 1.5 | 1.69 | 2.02 |
BaSO4 aq | 0 | 4.31 | 4.26 | 4.19 |
Ca | 2 | 95.25 | 94.79 | 94.09 |
CaCO3 aq | 0 | 0.92 | 1.04 | 1.27 |
CaHCO3 | 1 | 2.4 | 2.69 | 3.1 |
CaSO4 aq | 0 | 1.44 | 1.48 | 1.54 |
HCO3 | -1 | 94.58 | 94.63 | 94.58 |
H2CO3 aq | 0 | 3.19 | 2.89 | 2.55 |
Fe | 2 | 0.38 | 0.3 | 0.2 |
Fe(OH)2 | 1 | 18.83 | 15.92 | 12.34 |
Fe(OH)3 | 0 | 77.68 | 80.02 | 82.5 |
Fe(OH)4 | -1 | 2.92 | 3.6 | 4.84 |
Mg | 2 | 94.95 | 94.74 | 94.37 |
MgCO3 aq | 0 | 0.48 | 0.58 | 0.74 |
MgHCO3 | 1 | 3.38 | 3.36 | 3.35 |
MgSO4 aq | 0 | 1.19 | 1.32 | 1.54 |
Mn | 2 | 53.87 | 51.89 | 49.14 |
MnCO3 aq | 0 | 30.48 | 33.06 | 36.64 |
MnHCO3 | 1 | 15.03 | 14.4 | 13.53 |
MnSO4 aq | 0 | 0.58 | 0.61 | 0.64 |
Zn | 2 | 29.8 | 27.35 | 24.11 |
ZnCO3 | 0 | 42.36 | 43.76 | 45.15 |
Zn(CO3)2 | -2 | 13.01 | 15.17 | 18.39 |
ZnHCO3 | 1 | 11.75 | 10.72 | 9.37 |
Zn(OH)2 | 0 | 2.1 | 1.92 | 1.69 |
W wyniku analizy form występowania substancji w poszczególnych temperaturach, można stwierdzić, że w badanych wodach wapń i magnez występują głównie w formach Ca2+ i Mg2+, natomiast żelazo głównie w formach Fe(OH)3 i Fe(OH)2+, w zróżnicowanych proporcjach. Mangan w ponad 50-60% występuje w formie Mn2+, a w kilkunastu procentach w formie węglanów. W wyniku symulacji wykonanych dla 3 wartości temperatur (naturalna, ok. 5ºC i ok. 15ºC) można uznać, że w większości punktów dominujące formy substancji w analizowanych wodach występują w bardzo zbliżonych proporcjach. Jedynie znaczące zmiany zaobserwowano w analizie z punktu 1967, gdzie w najwyższej temperaturze (16,7ºC) wyraźnie dominuje forma Fe(OH)3 (66,69%), a w temperaturze 5ºC ta dominacja jest niewielka (53,73 %). Figura 1. Lokalizacja punktów opróbowania |
Tabela 7. Formy występowania głównych
|
Analiza stanu nasycenia badanych wód podziemnych względem faz mineralnych
Za pomocą programu WATEQ4F przeprowadzono również analizę stanów nasycenia wód (SI) względem tych faz mineralnych, które mogą wpływają znacząco na skład chemiczny wód podziemnych w badanych środowiskach wód podziemnych. Wyniki modelowania zestawiono w tabelach 8-14.
Tabela 8. Wskaźniki nasycenia wód (SI) względem znaczących faz mineralnych w wodzie z punktu 991
SI | ||||
Minerał | Wzór chemiczny | 5oC | 10.8oC | 15oC |
Aragonite | CaCO3 | -0.066 | 0.028 | 0.094 |
Barite | BaSO4 | -0.056 | -0.194 | -0.285 |
Calcite | CaCO3 | 0.093 | 0.183 | 0.245 |
Chalcedony | SiO2 | 0.395 | 0.319 | 0.266 |
Dolomite | MgCO3 | -1.142 | -0.209 | -0.011 |
Jarosite(s) | KFe3(SO4)2(OH)6 | 1.532 | 0.509 | -0.225 |
Magnesite | MgCO3 | -1.092 | -0.920 | -0.804 |
Quartz | SiO2 | 0.891 | 0.794 | 0.727 |
Rhodochrs | MnCO3 | -0.428 | -0.340 | -0.284 |
Siderite | FeCO3 | -0.409 | -0.411 | -0.423 |
Silica gel | Si | -0.153 | -0.225 | -0.275 |
SiO2 (a) | SiO2 | -0.517 | -0.571 | -0.609 |
Smithsonitte | ZnCO3 | -1.024 | -0.902 | -0.822 |
ZnCO3, 1H2 | -0.534 | -0.482 | -0.451 |
Tabela 9. Wskaźniki nasycenia wód (SI) względem znaczących faz mineralnych w wodzie z punktu 1154
SI | ||||
Minerał | Wzór chemiczny | 5oC | 9.6oC | 15oC |
Aragonite | CaCO3 | -0.134 | -0.058 | 0.026 |
Barite | BaSO4 | 0.378 | 0.269 | 0.150 |
Calcite | CaCO3 | 0.025 | 0.097 | 0 .177 |
Chalcedony | SiO2 | 0.328 | 0.267 | 0.199 |
Dolomite | MgCO3 | -0.895 | -0.658 | -0.402 |
Greenalite | (Fe2+Fe3+)23Si2O5(OH)4 | -0.734 | -0.985 | -1.292 |
Gypsum | CaSO4 x 2H2O | -1.877 | -1.895 | -1.912 |
Jarosite(s | KFe3(SO4)2(OH)6 | 2.631 | 1.845 | 0.930 |
Magnesite | MgCO3 | -1.414 | -1.277 | -1.127 |
Quartz | SiO2 | 0.823 | 0.747 | 0.660 |
Rhodochrsite | MnCO3 | -0.200 | -0.128 | -0.051 |
Siderite | FeCO3 | -0.587 | -0.611 | -0.650 |
Silica gel | Si | -0.220 | -0.277 | -0.342 |
SiO2 (a) | SiO2 | -0.584 | -0.627 | -0.676 |
Tabela 10. Wskaźniki nasycenia wód (SI) względem znaczących faz mineralnych w wodzie z punktu 1967
SI | ||||
Minerał | Wzór chemiczny | 5oC | 10oC | 16.7oC |
Adularia | KALSi3O8 | 0.368 | -0.050 | -0.596 |
Allophane | Al2O3(SiO2)1,3-2H2O | 0.430 | -0.016 | -0.593 |
Aragonite | CaCO3 | 0.095 | 0.177 | 0.280 |
Barite | BaSO4 | 0.533 | 0.415 | 0.271 |
Calcite | CaCO3 | 0.254 | 0.332 | 0.430 |
Chalcedony | SiO2xH2O | 0.180 | 0.114 | 0.030 |
Dolomite | MgCO3 | -0-.324 | -0.070 | 0.240 |
Gibbsite | Al. (OH)3 | 1.261 | 0.995 | 0.646 |
Jarosite | KFe3(SO4)2(OH)6 | -0.008 | -0.414 | -0.968 |
Magnesite | MgCO3 | -1.072 | -0.925 | -0.745 |
Quartz | SiO2 | 0.676 | 0.592 | 0.485 |
Rhodochrs | MnCO3 | -0.120 | -0.045 | 0.044 |
Siderite | FeCO3 | -0.668 | -0.677 | -0.707 |
Silica gel | Si | -0.368 | -0.430 | -0.510 |
SiO2 (a) | SiO2 | -0.732 | -0.779 | -0.839 |
Tabela 11. Wskaźniki nasycenia wód (SI) względem znaczących faz mineralnych w wodzie z punktu 2149
SI | ||||
Minerał | Wzór chemiczny | 5oC | 12oC | 15oC |
Allophane | Al2O3(SiO2)1,3-2H2O | 0.790 | 0.105 | -0.178 |
Aragonite | CaCO3 | 0.111 | 0.225 | 0.271 |
Barite | BaSO4 | -0.775 | -0.937 | -1.000 |
Calcite | CaCO3 | 0.270 | 0.378 | 0.423 |
Chalcedony | SiO2xH2O | -0.007 | -0.099 | -0.137 |
Dolomite | MgCO3 | -0.320 | 0.033 | 0.172 |
Gibbsite | Al.(OH)3 | 0.427 | 0.042 | -0.116 |
Kaolinite | Al4(OH)8 | 2.568 | 1.597 | 1.198 |
Laumontite | Ca AlSiO2 H2O | 0.418 | -0.352 | -0.672 |
Magnesite | MgCO3 | -1.084 | -0.880 | -0.798 |
Montmoril | Ca | 0.913 | -0.216 | -0.683 |
Quartz | SiO2 | 0.489 | 0.372 | 0.324 |
Silica gel | Si | -0.555 | -.0642 | -0.678 |
SiO2 (a) | SiO2 | -0.919 | -0.985 | -1.012 |
Tabela 12. Wskaźniki nasycenia wód (SI) względem znaczących faz mineralnych w wodzie z punktu 2346
SI | ||||
Minerał | Wzór chemiczny | 5oC | 9.8oC | 15oC |
Allophane | Al2O3(SiO2)1,3-2H2O | 0.296 | -0.212 | -0.742 |
Aragonite | CaCO3 | 0.101 | 0.180 | 0.260 |
Barite | BaSO4 | -0.007 | -0.118 | -0.228 |
Calcite | CaCO3 | 0.261 | 0.335 | 0.411 |
Chalcedony | SiO2 x H2O | -0.372 | -0.437 | -0.505 |
Dolomite | MgCO3 | -0.462 | -0.219 | 0.024 |
Gibbsite | Al.(OH)3 | -0.235 | -0.502 | -0.779 |
Kaolinite | Al4(OH)8 | 0.512 | -0.162 | -0.864 |
Quartz | SiO2 | 0.123 | 0.042 | -0.044 |
Rhodochrs | MnCO3 | -0.366 | -0.310 | -0.257 |
Silica gel | Si | -0.921 | -0.981 | -1.045 |
Tabela 13. Wskaźniki nasycenia wód (SI) względem znaczących faz mineralnych w wodzie z punktu 1088
SI | ||||
Minerał | Wzór chemiczny | 5oC | 9oC | 15oC |
Aragonite | CaCO3 | 0.397 | 0.462 | 0.555 |
Barite | BaSO4 | -0.846 | -0.94 | -1.071 |
Calcite | CaCO3 | 0.556 | 0.618 | 0.706 |
Dolomite | MgCO3 | -0.023 | 0.182 | 0.466 |
Magnesite | MgCO3 | -1.074 | -0.954 | -0.788 |
Rhodochrs | MnCO3 | -0.014 | 0.037 | 0.105 |
Smithsonit | ZnCO3 | -0.846 | -0.784 | -0.7 |
ZnCO3, 1H2 | O | -0.357 | -0.342 | -0.329 |
Na podstawie wyników modelowania można zaobserwować, że jednym z minerałów mogących wpływać na skład chemiczny badanych wód jest kalcyt (CaCO3). We wszystkich badanych środowiskach SI kalcytu jest bliskie zeru, czyli woda dąży do stanu nasycenia. Oznacza to, że jest on jednym z minerałów kontrolujących równowagę chemiczną w wodach o temperaturze ok. 10ºC, ale wartość wskaźnika nasycenia SI>0 oznacza, że kalcyt ma niewielką tendencję do wytrącania. Przy spadku temperatury SI kalcytu minimalnie maleje, a następnie przy wzroście temperatury SI minimalnie wzrasta, ale we wszystkich 3 analizowanych zakresach temperatur pozostaje niewielka tendencja do wytrącania. Przy tak małej różnicy temperatur, kalcyt nie będzie wpływał na chemizm wód. Również aragonit kontroluje równowagę w tych wodach, ale w przypadku studni 991, 1154 i 7209 nawet taki niewielki wzrost czy spadek temperatury spowoduje zmianę tendencji do rozpuszczania w przypadku obniżenia temperatury do 5ºC lub do wytrącania w przypadku jej podwyższenia do 15ºC. Analogicznie będzie się zachowywał dolomit (MgCO3) czy rodochrozyt (MnCO3).
Tabela 14. Wskaźniki nasycenia wód (SI) względem znaczących faz mineralnych w wodzie z punktu 7209
SI | ||||
Minerał | Wzór chemiczny | 5oC | 9.2oC | 15oC |
Aragonite | CaCO3 | -0.02 | 0.049 | 0.141 |
Barite | BaSO4 | -0.423 | -0.522 | -0.649 |
Calcite | CaCO3 | 0.139 | 0.205 | 0.292 |
Dolomite | MgCO3 | -0.73 | -0.514 | -0.237 |
Jarosite | KFe3(SO4)2(OH)6 | 0.428 | -0.354 | -1.419 |
Rhodochrs | MnCO3 | -0.072 | -0.009 | 0.069 |
Siderite | FeCO3 | -1.385 | -1.427 | -1.498 |
Wśród analizowanych minerałów siarczanowych występujących w badanych środowiskach, w stanie bliskim równowagi z wodą jest baryt. W badanych wodach ma on niewielką tendencję do rozpuszczania lub wytrącania (1154 i 1967), ale te tendencję nie zmieniają się wraz ze wzrostem czy spadkiem temperatury.
Krzemiany występują w stanie bliskim do równowagi z badanymi wodami we wszystkich punktach, co potwierdza ich niską rozpuszczalność w panujących w badanym rejonie warunkach Eh i pH oraz w analizowanym zakresie temperatur.
Zawartość żelaza jest związana z mającymi wyraźne tendencje do wytrącania tlenkami i wodorotlenkami żelaza, bez zmiany tego trendu w analizowanym zakresie temperatur. Ale warto tu wspomnieć, że każdy znaczący wzrost pH obniżający kwasowość roztworu oraz zmiana warunków redoks w środowisku wód podziemnych na wyraźnie utleniające pobudzać będą wytrącanie się tlenków żelaza m.in. w postaci getytu czy limonitu. Wśród pozostałych minerałów żelaza, jarosyt wykazuje zmianę tendencji ze wzrostem temperatury, z wytrącania do rozpuszczania, ale ta forma żelaza występuje w tych wodach w minimalnej ilości i nie wpłynie na kształtowanie się składu chemicznego wód.
Minerały manganu (np. rodochrozyt czy manganit) wykazują tendencje do rozpuszczania. Zmiany temperatur w analizowanym zakresie nie powodują zmian form występowania manganu w badanych wodach. Ale każdy znaczący wzrost pH obniżający kwasowość roztworu oraz zmiana warunków redoks w środowisku wód podziemnych na wyraźnie utleniające może pobudzać procesy wytrącania związków manganu w postaci zarówno tlenków, węglanów czy siarczanów.
Wyniki modelowania hydrogeochemicznego wody podziemnej z otworu Musuły
W celu wstępnej analizy wpływu zmian temperatury i potencjału redox na procesy hydrogeochemiczne w wodach podziemnych w zasięgu działania systemu ATES wykonano modelowanie hydrogeochemiczne za pomocą programów WATEQ4F oraz PHREEQC dla analizy fizyko-chemicznej wody z otworu Musuły (tab. 1). Głębokość ujętego poziomu wodonośnego wynosi 95,4−120,8 m (piaski i żwiry), a ustalonego zwierciadła wody 1,97 m.
Tabela 1. Wyniki analizy fizyko-chemicznej otworu Musuły
Składnik | Jednostka | Zawartość | |||
PEW | µS/cm | 500 | |||
Substancje rozpuszczone | mg/l | 495 | |||
Zasadowość ogólna [CaCO3] | mg/l | 280 | |||
Twardość ogólna [CaCO3] | mg/l | 292 | |||
Temperatura | oC | 9,7 | |||
pH | 7,2 | ||||
Eh | V | -0,155 | |||
HCO3 | mg/l | 342 | |||
Cl | mg/l | 4,1 | |||
NO3 | mg/l | 0,04 | |||
SO4 | mg/l | 2,17 | |||
NH4+ | mg/l | 0,22 | |||
Ca | mg/l | 89,9 | |||
Mg | mg/l | 16,5 |
Składnik | Jednostka | Zawartość |
Na | mg/l | 6,3 |
K | mg/l | 1,9 |
Fe | mg/l | 3,25 |
Mn | mg/l | 0,211 |
Li | mg/l | 0,0079 |
Cu | mg/l | 0,00025 |
B | mg/l | 0,02 |
Ba | mg/l | 0,152 |
SiO2 | mg/l | 22,5 |
Sr | mg/l | 0,165 |
ChZT | mg/l | 10,7 |
TOC | mg/l | 3,7 |
CO2 | mg/l | 2,94 |
DOX | mg/l | 0,31 |
Formy występowania poszczególnych składników chemicznych w badanych wodach podziemnych przy różnych temperaturach
Poszczególne składniki rozpuszczone w wodach mogą występować w różnych postaciach; jako uwodnione wolne jony, w formie związków kompleksowych, par jonowych, kompleksowych związków organicznych czy cząsteczek obojętnych (np. rozpuszczony tlen). Zależeć to będzie od rozpuszczalności danej substancji w wodzie, od jej zdolności do tworzenia kompleksów i par jonowych, od pH i Eh wody, jej temperatury, siły jonowej oraz składu chemicznego.
Za pomocą programu Wateq4f przeprowadzono modelowanie hydrogeochemiczne analizy wody z otworu Musuły. W pierwszym etapie wykonano modelowanie zakładając różne wielkości temperatury wody, czyli temperaturę pierwotną (naturalną) 9,7ºC oraz temperaturę schłodzonej wody (5ºC), jak również temperatury podgrzanej wody (20, 25 i 30ºC). W tabeli nr 2, przedstawiono wyniki modelowania, czyli podstawowe formy występowania głównych składników chemicznych w badanej wodzie podziemnej w różnych temperaturach.
W wyniku analizy form występowania substancji w wodach w poszczególnych temperaturach, można stwierdzić, że w badanych wodach główne składniki występują w formach typowych dla wód podziemnych na tej głębokości i te formy nie zmieniają się znacząco wraz ze zmianą temperatur w modelowanym zakresie. Natomiast żelazo głownie w formie Fe2+ (około 70%), FeHCO3+ (około 25%) i w niewielkim procencie FeCO3,, i te proporcje zmieniają się nieznacznie, podobnie jak w przypadku manganu, który w ponad 60% występuje w formie Mn2+, a w ponad 30% w formie węglanów.
W wyniku symulacji wykonanych dla 5 wartości temperatur można uznać, że dominujące formy substancji w analizowanej wodzie występują w bardzo zbliżonych proporcjach.
Analiza stanu nasycenia badanych wód podziemnych względem faz mineralnych przy różnych temperaturach
Za pomocą programu WATEQ4F przeprowadzono również analizę stanów nasycenia wód (SI) względem tych faz mineralnych, które mogą wpływać znacząco na skład chemiczny wód podziemnych w badanym środowisku wód podziemnych. Wyniki modelowania zestawiono w tabeli 3.
Na podstawie wyników modelowania można zaobserwować, że jednym z minerałów mogących wpływać na skład chemiczny badanych wód jest kalcyt (CaCO3). We wszystkich założonych temperaturach SI kalcytu jest bliskie zeru, czyli woda dąży do stanu nasycenia. Oznacza to, że jest on jednym z minerałów kontrolujących równowagę chemiczną w wodzie o temperaturze ok. 10 ºC, ale wartość wskaźnika nasycenia SI>0 oznacza, że kalcyt ma niewielką tendencję do wytrącania. Przy spadku temperatury SI kalcytu minimalnie maleje, a następnie przy wzroście temperatury SI minimalnie wzrasta. Przy wzroście temperatury do 20ºC i ponad, pojawia się tendencja do wytrącania. Oznaczać to może, że przy wprowadzeniu do poziomu wodonośnego zimnej wody kalcyt będzie się rozpuszczał i dostawał do wody, a przy wprowadzeniu ciepłej wody (o temperaturze ponad 20ºC) może się wytrącać powodując osadzanie się wapnia (powstawanie tzw. kamienia – rekrystalizacja). Powstawanie kamienia jest problemem w studniach iniekcyjnych, gdzie wytrącane minerały mogą szybko doprowadzić do wzrostu ciśnienia wtrysku i ewentualnie konieczności prac oczyszczających lub wiercenia dodatkowych otworów.
Tabela 2. Formy występowania głównych składników chemicznych w wodzie przy różnych temperaturach
Forma składnika | % zawartość składnika | |||||
Temperatura oC | 5.0 | 9.7 | 20.0 | 25.0 | 30.0 | |
H3BO3 aq | 0 | 99.36 | 99.28 | 99.09 | 99.00 | 98.91 |
Ba | 2 | 7.64 | 97.35 | 96.56 | 96.08 | 95.53 |
Ca | 2 | 96.73 | 96.27 | 95.32 | 94.92 | 94.56 |
CaHCO3 | 1 | 2.93 | 3.36 | 4.19 | 4.52 | 4.8 |
Cl | -1 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 |
HCO3 | -1 | 82.60 | 83.89 | 85.71 | 86.25 | 86.63 |
H2CO3 aq | 0 | 15.57 | 14.09 | 11.84 | 11.11 | 10.56 |
Cu | 1 | 99.63 | 99.63 | 99.64 | 99.64 | 99.64 |
CuCO3 aq | 0 | 0 | 0 | 0 | 0 | 0 |
CuHCO3 | 1 | 0 | 0 | 0 | 0 | 0 |
Cu(OH)2 | 0 | 0 | 0 | 0 | 0 | 0 |
Fe | 2 | 70.71 | 70.49 | 70.00 | 69.75 | 69.49 |
FeCO3 aq | 0 | 2.56 | 2.94 | 3.77 | 4.16 | 4.52 |
FeHCO3 | 1 | 26.60 | 26.41 | 25.96 | 25.74 | 25.51 |
Fe (OH)2 | 1 | 0 | 0 | 0 | 0 | 0 |
Fe(OH)3 | 0 | 0 | 0 | 0 | 0 | 0 |
K | 1 | 99.99 | 99.99 | 99.99 | 99.99 | 99.99 |
Li | 1 | 99.99 | 99.99 | 99.99 | 99.99 | 99.99 |
Mg | 2 | 95.68 | 95.65 | 95.50 | 95.39 | 95.25 |
Mn | 2 | 68.65 | 67.87 | 66.20 | 65.43 | 64.72 |
MnCO3 aq | 0 | 8.24 | 9.37 | 11.8 | 12.91 | 13.95 |
MnHCO3 | 1 | 23.02 | 22.66 | 21.88 | 21.52 | 21.17 |
Na | 1 | 99.74 | 99.74 | 99.74 | 99.74 | 99.73 |
NH4 | 1 | 99.80 | 99.72 | 99.42 | 99.17 | 98.84 |
NO3 | -1 | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 |
H4SiO4aq | 0 | 99.88 | 99.86 | 99.78 | 99.74 | 99.70 |
SO4 | -2 | 82.30 | 81.40 | 79.39 | 78.38 | 77.36 |
Sr | 2 | 96.97 | 96.48 | 95.13 | 94.31 | 93.36 |
SrHCO3 | 1 | 0.13 | 3.31 | 4.6 | 5.38 | 6.29 |
Tabela 3. Wskaźniki nasycenia wód (SI) względem znaczących faz mineralnych w wodzie przy różnych temperaturach
IFaza mineralna | Wzór | SI | ||||
Temperatura oC | 5.00 | 9.7 | 20.0 | 25.0 | 30.0 | |
Eh V | -0.155 | -0.155 | -0.155 | -0.155 | -0.155 | |
Aragonite | CaCO3 | -0.243 | -0.166 | -0.007 | 0.066 | 0.136 |
Barite | BaSO4 | -0.665 | -0.778 | -0.996 | -1.088 | -01.173 |
Calcite | CaCO3 | -0.084 | -0.010 | 0.141 | 0.210 | 0.276 |
Chalcedony | SiO2 | 0.374 | 0.312 | 0.184 | 0.125 | 0.067 |
Christobalite | SiO2 | 0.451 | 0.379 | 0.229 | 0.161 | 0.094 |
Cuprite | Cu2O | -0.458 | -0.540 | -0.710 | -.789 | -0.865 |
Dolomite (c) | MgCO3 | -0.886 | -0.643 | -0.169 | 0.035 | 0.223 |
Fe(OH)2 | Fe(OH)2 | 0.398 | 0.569 | 0.923 | 1.086 | 1.243 |
Goethite | FeO(OH) | 1.712 | 1.882 | 2.237 | 2.400 | 2.557 |
Gypsum | CaSO4 *H2O | -3.113 | -3.131 | -3.161 | -3.171 | -3.178 |
Hematite | Fe2O3 | 3.806 | 4.550 | 6.097 | 6.808 | 7.496 |
Magnesite | MgCO3 | -1.296 | -1.155 | -0.877 | -0.756 | -0.644 |
Quartz | SiO2 | 0.870 | 0.792 | 0.629 | 0.554 | 0.481 |
Rhodochrs(d) | MnCO3 | -1.009 | -.953 | -0.853 | -0.813 | -0.780 |
Rhodochrs(c) | MnCO3 | -0.344 | -0.270 | -0.130 | -0.073 | -0.023 |
Siderite (d) | FeCO3 | 0.245 | 0.304 | 0.412 | 0.455 | 0.491 |
Siderite (c) | FeCO3 | 0.554 | 0.646 | 0.821 | 0.895 | 0.961 |
Silica gel | Silica gel | -0.174 | -0.232 | -0.353 | -0.408 | -0.462 |
Talc | Mg3(OH)2Si4O10 | -4.378 | -3.777 | -2.528 | -1.954 | -1.400 |
Również aragonit i dolomit kontrolują równowagę w tych wodach, i przy osiągnięciu temperatury 25ºC nastąpi zmiana tendencji podobnie jak w przypadku kalcytu. Nawet niewielki wzrost czy spadek temperatury spowoduje zmianę tendencji do rozpuszczania w przypadku obniżenia temperatury (do 5ºC) lub do wytrącania w przypadku jej podwyższenia (do 15ºC). Analogicznie będzie się zachowywał dolomit (MgCO3) czy rodochrozyt (MnCO3).
Wśród analizowanych minerałów siarczanowych występujących w badanych środowiskach, w stanie bliskim równowagi z wodą jest baryt. W badanych wodach ma on niewielką tendencję do rozpuszczania i ta tendencja nie zmienia się wraz ze wzrostem czy spadkiem temperatury.
Krzemiany występują w stanie bliskim do równowagi z badanymi wodami we wszystkich punktach, co potwierdza ich niską rozpuszczalność w panujących w badanym rejonie warunkach Eh i pH oraz w analizowanym zakresie temperatur.
Zawartość żelaza jest związana z mającymi wyraźne tendencje do wytrącania tlenkami, wodorotlenkami i węglanami żelaza, bez zmiany tego trendu w analizowanym zakresie temperatur.
Minerały manganu (np. rodochrozyt) wykazują tendencje do rozpuszczania. Zmiany temperatur w analizowanym zakresie nie powodują zmian tendencji w badanych wodach.
Formy występowania poszczególnych składników chemicznych w badanych wodach podziemnych przy różnych wartościach potencjału redox
Przeprowadzono również modelowanie hydrogeochemiczne tej analizy przy zmieniających się warunkach redukcyjno-oksydacyjnych. Naturalny potencjał redoks pomierzony in situ wynosi w tych wodach −0,155 V, czyli reprezentuje on słabo redukcyjne warunki – do modelowania wybrano 2 analizy: naturalną o temperaturze 9,7ºC oraz o założonej temperaturze 20ºC oraz założono dodatkową wartość potencjału redoks 0,2 V (warunki słabo utleniające, typowe dla płytkich wód). Wyniki zestawiono w tabelach 4 i 5.
W wyniku symulacji wykonanych dla 2 wartości temperatur i 2 wartości potencjału redox można uznać, że dominujące formy substancji w analizowanych wodach występują w bardzo zbliżonych proporcjach. Jedynie znaczące zmiany zaobserwowano w przypadku miedzi, gdzie przy Eh wynoszącym 0,2 V pojawiły się nowe formy z dominującymi węglanami oraz wodorotlenkami. W przypadku form żelaza, przy wyższym Eh pojawiły się niewielkie ilości wodorotlenków żelaza.
Węglany wapnia czy magnezu (kalcyt aragonit, dolomit, magnezyt) nie zmieniają tendencji przy założonej zmianie warunków redox, podobnie siarczany i krzemiany.
Zmiana warunków redoks (Eh) z −0.155 V (warunki słabo redukcyjne) na 0.20 V (warunki słabo utleniające) powoduje zmianę wskaźnika nasycenia tlenków miedzi, ze słabej tendencji do rozpuszczania (prawie równowaga), na silną tendencję do rozpuszczania.
W wodach głębinowych, gdzie nie ma zbyt dużo tlenu, dominuje forma żelaza na drugim stopniu utleniania. Jest to żelazo w postaci rozpuszczalnej i w takiej formie potrafi osiągać bardzo duże stężenia, jednak nie zmienia barwy i nie przyczynia się do powstawania mętności wody. W analizowanym modelu w przypadku wodorotlenków i tlenków żelaza nastąpiła zmiana ze stanu równowagi lub tendencji do rozpuszczania przy naturalnym Eh wynoszącym −0,155 V, na tendencję do wytrącania przy Eh wynoszącym 0,2 V. Przy czy nie ma znaczenia czy to była przyjęta temperatura 9,7ºC czy 20ºC. W wodach, które mają łatwy dostęp do tlenu, żelazo przybiera trzeci stopień utleniania. W takiej postaci bardzo łatwo wytrąca się z wody w formie nierozpuszczalnych cząsteczek o bardzo charakterystycznym rdzawym zabarwieniu i już przy relatywnie niskich stężeniach ma wpływ na pogorszenie barwy i powstanie mętności wody. Żelazo nie jest pożądane ze względów technicznych. Ma ono tendencję do wytrącania się z wody w formie osadu, który z czasem twardnieje i z łatwością gromadzi się w całej instalacji. W rezultacie prowadzi to do: zmniejszenia światła przewodów bądź ich całkowitego zamknięcia, spadków wydajności, zmniejszenia przepływu, fragmentacji osadów i przedostawania się ich w głąb instalacji, powstawaniem korozji pod osadem oraz drogich w naprawie awarii
Minerały manganu (np. rodochrozyt czy manganit) wykazują tendencje do rozpuszczania. Pomimo, że znacząca zmiana warunków redoks w środowisku wód podziemnych na wyraźnie utleniające może pobudzać procesy wytrącania związków manganu, w przypadku przyjętej w tym modelowaniu amplitudy zmiany Eh nie nastąpiła zmiana tendencji.
Tabela 4. Formy występowania głównych składników chemicznych w wodzie w temperaturze naturalnej (9,70C) i 20 0C przy różnych wartościach potencjału redox (Eh)
Forma składnika | % zawartość składnika | ||||
Temperatura oC | 9.7 | 9.7 | 20.0 | 20.0 | |
Eh [V] | -0.155 | 0.2 | -0.155 | 0.2 | |
H3BO3 aq | 0 | 99.28 | 99.28 | 99.09 | 99.09 |
Ba | 2 | 97.35 | 97.35 | 96.56 | 96.56 |
Ca | 2 | 96.27 | 96.26 | 95.32 | 95.32 |
CaHCO3 | 1 | 3.36 | 3.36 | 4.19 | 4.19 |
Cl | -1 | 100.00 | 100.00 | 100.00 | 100.00 |
HCO3 | -1 | 83.89 | 83.89 | 85.71 | 85.73 |
H2CO3 aq | 0 | 14.09 | 14.09 | 11.84 | 11.84 |
Cu | 1 | 99.63 | 0.58 | 99.64 | 0.73 |
CuCO3 aq | 0 | 0 | 58.15 | 0 | 64.1 |
CuHCO3 | 1 | 0 | 11.7 | 0 | 9.89 |
Cu(OH)2 | 0 | 0 | 22.31 | 0 | 18.94 |
Fe | 2 | 70.49 | 69.34 | 70.00 | 64.49 |
FeCO3 aq | 0 | 2,94 | 2.89 | 3.77 | 3.47 |
FeHCO3 | 1 | 26.41 | 25.98 | 25.96 | 23.93 |
Fe (OH)2 | 1 | 0 | 0.85 | 0 | 3.14 |
Fe(OH)3 | 0 | 0 | 0.78 | 0 | 4.66 |
K | 1 | 99.99 | 99.99 | 99.99 | 99.99 |
Li | 1 | 99.99 | 99.99 | 99.99 | 99.99 |
Mg | 2 | 95.65 | 95.65 | 95.50 | 95.50 |
Mn | 2 | 67.87 | 67.87 | 66.20 | 66.19 |
MnCO3 aq | 0 | 9.37 | 9.37 | 11.8 | 11.8 |
MnHCO3 | 1 | 22.66 | 22.66 | 21.88 | 21.89 |
Na | 1 | 99.74 | 99.74 | 99.74 | 99.74 |
NH4 | 1 | 99.72 | 99.72 | 99.42 | 99.42 |
NO3 | -1 | 100.00 | 100.00 | 100.00 | 100.00 |
H4SiO4aq | 0 | 99.86 | 99.86 | 99.78 | 99.78 |
SO4 | -2 | 81.40 | 81.40 | 79.39 | 79.4 |
Sr | 2 | 96.48 | 96.48 | 95.13 | 95.13 |
SrHCO3 | 1 | 3.31 | 3.31 | 4.6 | 4.6 |
Tabela 5. Wskaźniki nasycenia wód (SI) względem znaczących faz mineralnych w wodzie w temperaturze naturalnej (9,70C) i 20 0C przy różnych wartościach potencjału redox (Eh)
Faza mineralna | Wzór | SI | |||
Temperatura oC | 9.7 | 9.7 | 20.0 | 20.0 | |
Eh [V] | -0.155 | 0.2 | -0.155 | 0.2 | |
Aragonite | CaCO3 | -0.166 | -0.166 | -0.007 | -0.007 |
Barite | BaSO4 | -0.778 | -0.778 | -0.996 | -0.996 |
Calcite | CaCO3 | -0.010 | -0.010 | 0.141 | 0.141 |
Chalcedony | SiO2 | 0.312 | 0.312 | 0.184 | 0.184 |
Christobalite | SiO2 | 0.379 | 0.379 | 0.229 | 0.229 |
Cuprite | Cu2O | -0.540 | -5.01 | -.710 | -4.985 |
Dolomite (c) | MgCO3 | -0.643 | -0.643 | -0.169 | -0.169 |
Fe(OH)2 | Fe(OH)2 | 0.569 | 6.89 | .923 | 6.99 |
Ferrihydrite | FeO-OH | -4.18 | 2.31 | -3.65 | 2.41 |
Goethite | FeO(OH) | 1.882 | 8.20 | 2.237 | 8.30 |
Gypsum | CaSO4 *H2O | -3.131 | -3.131 | -3.161 | -3.161 |
Hematite | Fe2O3 | 4.550 | 17.18 | 6.097 | 18.23 |
Magnesite | MgCO3 | -1.155 | -1.155 | -0.877 | -0.877 |
Quartz | SiO2 | 0.792 | 0.792 | 0.629 | 0.629 |
Rhodochrs(d) | MnCO3 | -0.953 | -0.953 | -0.853 | -0.852 |
Rhodochrs(c) | MnCO3 | -0.270 | -0.270 | -0.130 | -0.130 |
Siderite (d) | FeCO3 | 0.304 | 0.297 | 0.412 | 0.377 |
Siderite (c) | FeCO3 | 0.646 | 0.639 | 0.821 | 0.786 |
Silica gel | SiO2 | -0.232 | -0.232 | -0.353 | -0.353 |
Talc | Mg3(OH)2Si4O10 | -3.777 | -3.777 | -2.528 | -2.528 |
Modelowanie mieszania wód o różnych temperaturach
Za pomocą programu PHREEQC przeprowadzono również modelowanie mieszania wód o temperaturze naturalnej (9,7ºC) i temperaturze wyraźnie podwyższonej (25ºC). Wyniki analizy stanów nasycenia wód (SI) względem tych faz mineralnych, które mogą wpływać znacząco na skład chemiczny wód podziemnych w badanym środowisku wód podziemnych zestawiono w tabeli 6.
Wymieszanie wód o temperaturze 9,7ºC i podgrzanych o 15ºC nie spowodowało zmian tendencji w przypadku prawie wszystkich analizowanych minerałów. Jedynie w przypadku aragonitu i dolomitu, w stosunku do wody o naturalnej temperaturze, nastąpiła zmiana z tendencji do wytracania na tendencję do rozpuszczania.
Tabela 6. Wyniki mieszania wód o różnych temperaturach
Faza mineralna | SI | ||
Solution1 | Solution2 | Mixture 0.5/0.5 | |
Temperatura oC | 25 | 9.7 | |
Anhydrite | -3.40 | -3.40 | -3.41 |
Aragonite | -0.09 | 0.14 | 0.02 |
Barite | -0.78 | -1.09 | -0.95 |
Calcite | 0.07 | 0.29 | 0.17 |
Celestite | -4.17 | -4.21 | -4.20 |
Chalcedony | -0.02 | -0.21 | -0.11 |
Chrysotile | -8.57 | -6.59 | -7.63 |
CO2(g) | -1.69 | -1.60 | -1.63 |
Dolomite | -0.49 | 0.19 | -0.16 |
Fe(OH)3(a) | -3.89 | -3.52 | -4.55 |
Goethite | 1.43 | 2.37 | 1.06 |
Gypsum | -3.14 | -3.18 | -3.17 |
H2(g) | -9.16 | -9.16 | -7.51 |
H2O(g) | -1.93 | -1.51 | -1.71 |
Halite | -9.13 | -9.17 | -9.15 |
Hausmannite | -29.99 | -26.06 | -29.70 |
Hematite | 4.78 | 6.76 | 4.10 |
Faza mineralna | SI | ||
Solution1 | Solution2 | Mixture 0.5/0.5 | |
Temperatura oC | 25 | 9.7 | |
Jarosite-K | -24.79 | -22.46 | -26.11 |
Manganite | -12.14 | -12.16 | -0.08 |
Melanterite | -7.06 | -7.30 | -13.00 |
NH3(g) | -9.48 | -8.67 | -9.25 |
O2(g) | -70.22 | -64.80 | -70.74 |
Pyrochroite | -6.58 | -6.60 | -6.61 |
Pyrolusite | -26.18 | -23.62 | -26.54 |
Quartz | 0.46 | 0.22 | 0.34 |
Rhodochrosite | -0.21 | -0.02 | -0.12 |
Sepiolite | -5.36 | -4.95 | -5.20 |
Sepiolite(d) | -7.84 | -7.85 | -7.89 |
Siderite | 0.71 | 0.95 | 0.83 |
SiO2(a) | -0.91 | -1.04 | -0.98 |
Strontianite | -2.13 | -2.00 | -2.08 |
Talc | -5.12 | -3.30 | -4.26 |
Witherite | -3.01 | -2.93 | -2.99 |
Podsumowanie
W wyniku symulacji wykonanych dla 5 wartości temperatur (5, 9, 7, 20, 25 i 30ºC) można uznać, że dominujące formy substancji w analizowanej wodzie występują w bardzo zbliżonych proporcjach, we wszystkich symulowanych wartościach temperatur, natomiast pojawiły się różnice w zachowaniu (tendencji) faz mineralnych w poszczególnych temperaturach:
- w przypadku węglanów istotna okazuje się granica pomiędzy temperaturą 20 a 25ºC, które z tendencji do rozpuszczania przechodzą w tendencję do wytrącania. Może to być źródłem powstawania „kamienia” w instalacji ATES;
- wraz ze wzrostem temperatury zwiększa się tendencja związków żelaza do wytrącania, ale równomiernie, bez wyraźnych "skoków";
- w przypadku krzemionki czy krzemianów sytuacja jest odwrotna, ze wzrostem temperatury zwiększa się tendencja do rozpuszczania, ale równomiernie, bez wyraźnych "skoków".
Zmiana warunków redoks (Eh) z −0,155 V (warunki słabo redukcyjne) na 0,20 V (warunki słabo utleniające) powoduje zmianę wskaźnika nasycenia tlenków miedzi, ze słabej tendencji do rozpuszczania (prawie równowaga), na silną tendencję do rozpuszczania. W przypadku wodorotlenków i tlenków żelaza następuje zmiana ze stanu równowagi lub tendencji do rozpuszczania na tendencję do wytrącania, co może być przyczyną kolmatacji instalacji ATES.
Węglany wapnia czy magnezu (kalcyt aragonit, dolomit, magnezyt) nie zmieniają tendencji przy założonej zmianie warunków redox, podobnie jak siarczany i krzemiany.
-
Mapa potencjału
-
Baza wiedzy
-
Literatura i materiały
-
Masz pytanie?
Skontaktuj się z nami!
dr Monika Konieczyńska
tel. 22 45 92 298
mkon@pgi.gov.pl
Zakład Geologii Środowiskowej
Państwowy Instytut Geologiczny – Państwowy Instytut Badawczy
ul. Rakowiecka 4, 00-975 Warszawa
www.pgi.gov.pl