Kalkulator powstawania osadów / rozpuszczalności siarczanu wapnia

Metodyka i ograniczenia

Wzory, konwencje dotyczące danych wejściowych i założenia leżące u podstaw tej demonstracji dotyczącej siarczanu wapnia, z wyraźnym rozróżnieniem między interfejsem a termodynamicznym modelem powstawania osadów mineralnych.

1. Zakres, jednostki i dane wejściowe

Uproszczony skład to woda demineralizowana + NaCl + CaSO₄. Wartość wejściowa CaSO₄ to nominalna całkowita ilość dodanej soli, wyrażona w przeliczeniu na bezwodny CaSO₄ na litr roztworu początkowego. Nie jest to dawka gipsu, stężenie wolnych jonów ani zmierzone stężenie substancji rozpuszczonej po osiągnięciu równowagi. Punkt powyżej granicy przedstawia hipotetyczny skład przed wytrąceniem osadu.

Są to zakresy interfejsu, a nie zwalidowane zakresy pracy pod względem chemicznym. Nie wykonuje się obliczeń gęstości ani zmian objętości, które przeliczałyby litry roztworu początkowego na objętość roztworu w stanie równowagi lub kilogramy wody. Obliczenia dla stężonych solanek wymagają jednoznacznego określenia podstawy masowej i sposobu wyrażania stężenia.

2. Dokładna postać funkcji odpowiedzi nasycenia

Dla liczbowych danych wejściowych wyrażonych w powyższych jednostkach wyświetlany próg L (g/L w przeliczeniu na CaSO₄) wynosi:

L = max{0,5;
    (2,05 + 0,009 × T − 0,00015 × T²)
    × [1 + 0,014 × S / (1 + S / 110)]}

Wyświetlany stosunek stężenia do progu (%) = 100 × C / L

Współczynniki i dolną granicę 0,5 g/L przyjęto na potrzeby demonstracji. Kwadratowy człon temperaturowy osiąga maksimum przy 30 °C. Maksimum to jest narzucone przez wzór; nie stanowi przewidywanego przejścia gips–anhydryt ani uniwersalnego, doświadczalnie wyznaczonego maksimum rozpuszczalności.

Czynnik uwzględniający NaCl rośnie monotonicznie, dążąc do skończonej wartości. Nie może odtworzyć odwrócenia rzeczywistego trendu rozpuszczalności przy wysokim zasoleniu. Nie oblicza aktywności jonów, par jonowych wapń–siarczan, aktywności wody ani oddziaływań z innymi rozpuszczonymi solami.

Wartość L jest wyświetlana z dokładnością do dwóch miejsc po przecinku, a stosunek — po zaokrągleniu do najbliższego pełnego procenta. Etykiety stanu wykorzystują niezaokrąglony stosunek: poniżej 85% oznacza „Niedosycenie”, 85–100% — „Zbliżanie się do nasycenia”, a powyżej 100% — „Przesycenie”. Próg 85% przyjęto na potrzeby ostrzeżeń w interfejsie. Etykiety te opisują krzywą demonstracyjną, a nie termodynamiczny stan nasycenia analizowanej wody.

Dla szarej krzywej porównawczej przyjęto S = 0. Wprowadzona wartość CaSO₄ dla próbki pozostaje zachowana; „bez dodatku NaCl” nie oznacza wody pozbawionej jonów.

3. Granica wrzenia i górny limit temperatury

B (°C) = 100 + 28 × ln(P / 1,01325) + 0,025 × S
Tmax (°C) = floor[10 × min(150; B − 0,5)] / 10

Tutaj ln oznacza logarytm naturalny, a floor zaokrągla w dół. Funkcja ta ilustruje wpływ ciśnienia i zasolenia; nie jest korelacją opartą na tablicach parowych ani obliczeniem równowagi para–ciecz dla roztworu. Wpływ CaSO₄ na wrzenie jest pomijany. Margines 0,5 °C i górną granicę 150 °C przyjęto na potrzeby interfejsu.

Zmiana ciśnienia może zawęzić dopuszczalny zakres temperatury, ale ta wersja nie oblicza równowagi mineralnej zależnej od ciśnienia. Nie symuluje się parowania, zmian zasolenia, rozszerzalności cieplnej ani zmian stężenia wywołanych wrzeniem.

4. pH, specjacja i przeciwjony

Ustawienie pH jest w tej wersji jedynie zapisem wartości dla scenariusza, a nie parametrem chemicznym uwzględnianym w obliczeniach. Rzeczywisty model siarczanowy musi uwzględniać podział między jonami siarczanowymi i wodorosiarczanowymi, zwłaszcza w wodzie o odczynie kwaśnym, a także wolne jony wapnia, tworzenie par jonowych i poprawki na aktywność.

W przyszłym modelu dozowania obniżanie pH za pomocą HCl wprowadzałoby chlorki, a podwyższanie go za pomocą NaOH — sód. Dawki tych reagentów i ich przeciwjony należałoby uwzględnić w bilansie roztworu. Kwas siarkowy wprowadzałby również siarczany. Nie oblicza się tutaj tych dawek ani towarzyszących im jonów.

5. Gips, anhydryt i nasycenie termodynamiczne

Obliczenie termodynamiczne ocenia każdy minerał osobno, wykorzystując aktywności i stałą równowagi dla tej samej reakcji rozpuszczania:

Gips:   CaSO₄·2H₂O(s) ⇌ Ca²⁺ + SO₄²⁻ + 2H₂O
IAPgypsum = aCa²⁺ × aSO₄²⁻ × aw²

Anhydryt: CaSO₄(s) ⇌ Ca²⁺ + SO₄²⁻
IAPanhydrite = aCa²⁺ × aSO₄²⁻

SIphase = log₁₀[IAPphase / Kphase(T, P)]

Tutaj a oznacza aktywność danej formy chemicznej, a aw — aktywność wody; aktywności czystych minerałów przyjmuje się jako równe jedności. SI = 0 oznacza równowagę z określoną fazą, SI > 0 — przesycenie, a SI < 0 — niedosycenie. Obecny kalkulator nie oblicza tych wielkości ani nie wskazuje fazy mineralnej kontrolującej równowagę.

Termodynamiczny stosunek nasycenia IAP/K na ogół nie jest równy C/L. Całkowite stężenia soli nie mogą zastąpić aktywności wolnych jonów w tych wyrażeniach. Nawet prawidłowo obliczona dodatnia wartość SI wskazuje na tendencję termodynamiczną, a nie na masę osadu, szybkość tworzenia osadów mineralnych czy pewność ich powstania.

Przykład USGS dotyczący gipsu i anhydrytu w PHREEQC przedstawia obliczenia rozpuszczalności dla poszczególnych faz oraz porównanie z danymi doświadczalnymi. Zastosowanie do solanki wymaga odpowiedniego modelu aktywności, na przykład implementacji modelu Pitzera, oraz parametrów sprawdzonych dla jej składu i zakresu temperatury oraz ciśnienia. Demonstracja nie uruchamia PHREEQC.

6. Status walidacji i przykład umożliwiający odtworzenie wyników

Dla tej krzywej demonstracyjnej nie ustalono zbioru danych kalibracyjnych, termodynamicznego punktu odniesienia ani zakresu dokładności dla zastosowań inżynierskich. Sprawdzenie obliczeń arytmetycznych potwierdza zgodność z kodem; nie potwierdza zgodności z rzeczywistością fizyczną.

Przy ustawieniach początkowych (25 °C, NaCl 50 g/L, CaSO₄ 1,5 g/L i ciśnienie bezwzględne 1,013 bar) funkcje dają L = 3,2309766 g/L (wyświetlane jako 3,23), stosunek 46,4256% (wyświetlany jako 46%) oraz B = 101,2431 °C (wyświetlane jako 101,2). Górny limit temperatury wynosi 100,7 °C. Zmiana wyłącznie pH scenariusza nie zmienia tych wartości.

Przyszły model inżynierski należy zweryfikować przez porównanie z niezależnymi obliczeniami równowagi i odpowiednimi pomiarami doświadczalnymi. Należy zapisać wersję oprogramowania/bazy danych, skład roztworu, podstawę wyrażania stężeń, dopuszczone fazy mineralne i tolerancje akceptacji. Należy sprawdzić zarówno gips, jak i anhydryt w całym przewidzianym zakresie zasolenia i temperatury, zamiast polegać na jednym przykładzie z ustawieniami domyślnymi.

7. Znane ograniczenia i właściwe zastosowanie

• Brak obliczeń specjacji w roztworze wodnym, poprawek na aktywność, aktywności wody, rozwiązywania bilansu ładunków ani wyboru faz mineralnych.
• Brak niezależnej analizy wapnia i siarczanów oraz uwzględnienia magnezu lub innych jonów tła poza uproszczonymi danymi wejściowymi dotyczącymi soli.
• Brak bilansu masy dla parowania lub odzysku wody, obliczeń masy wytrąconego osadu, wpływu antyskalantów, zarodkowania ani kinetyki narastania osadów mineralnych.
• Brak zwalidowanych prognoz dotyczących minerałów dla wysokiej temperatury, wysokiego zasolenia lub zależności od ciśnienia.

Demonstracja w tym kalkulatorze szacuje wartość progową dla siarczanu wapnia i ilustruje tendencję do tworzenia osadów mineralnych; nie zastępuje pełnego modelowania geochemicznego ani badań pilotażowych w systemach o wysokim odzysku wody lub wysokiej temperaturze. Obliczenia mogą zawierać błędy lub pominięcia. Przed projektowaniem, wyborem metody uzdatniania lub podejmowaniem decyzji eksploatacyjnych należy niezależnie zweryfikować wyniki.

8. Źródła i dokumentacja modelu

• USGS — Równowagi fazowe i zależności rozpuszczalności: informacje o wpływie zasolenia i faz mineralnych.
• USGS — Dokumentacja PHREEQC w wersji 3: metody obliczania specjacji, równowagi mineralnej i modelowania aktywności.
• USGS — Przykład równowagi gipsu i anhydrytu.

Dokumentacja metodologii v1.0 · opracowana 24 października 2026 r. Dokumentuje istniejące funkcje demonstracyjne; nie wprowadza nowego modelu równowagi ani nie zmienia działania numerycznego kalkulatora.

Powrót do kalkulatora