Calcolatore della solubilità e delle incrostazioni della silice
Esplora la solubilità della silice amorfa quando l’acqua si riscalda, si raffredda o diventa più concentrata.
Le condizioni della tua acqua
Sposta un cursore o inserisci un valore.
5–95 °C · riferimento per acqua liquida in condizioni prossime a quelle atmosferiche.
Inserisci la somma delle specie monomeriche disciolte dell’acido silicico e dei silicati, espressa come SiO₂. Non inserire la silice totale che comprende colloidi o particelle. In acqua con basse concentrazioni di soluti, a temperature prossime a quella ambiente, mg/L equivale approssimativamente a mg/kg; qui non viene applicata alcuna conversione basata sulla densità.
pH 4–10. Ora il pH modifica la curva di saturazione della silice disciolta totale attraverso la speciazione calcolata di H₄SiO₄ / H₃SiO₄⁻ / H₂SiO₄²⁻. Vengono calcolate le correzioni dell’attività e le quantità minime di ioni necessarie al bilanciamento della carica; vedi i limiti del modello riportati di seguito.
Scenario di concentrazione
1–10×. Si assume che sia la silice sia il NaCl siano completamente trattenuti e che la massa d’acqua si riduca, prima della precipitazione. Il pH rimane al valore selezionato. La concentrazione effettiva di NaCl deve rimanere pari o inferiore a 5000 mg/kg. Non si tratta di una simulazione di una membrana RO o di un evaporatore.
0–5000 mg/kg di acqua come NaCl effettivamente presente, prima della fase di concentrazione. Valore predefinito: zero. Il NaCl modifica i coefficienti di attività; non sostituisce i solidi disciolti totali (TDS) né un’analisi completa dell’acqua. Gli altri sali disciolti non sono inclusi.
L’effetto del sale dipende dal pH e dalla temperatura.
Cosa succede quando l’acqua si riscalda o si raffredda?
Silice disciolta totale calcolata alla saturazione della silice amorfa, espressa come SiO₂ al pH selezionato. Il punto scuro rappresenta l’acqua selezionata dopo la fase di concentrazione ideale.
Blu: silice disciolta totale corretta per pH e NaCl · Grigio tratteggiato: riferimento della specie neutra in acqua pura · Punto scuro: acqua selezionata.
25 °C · pH 7.0 · 1× di concentrazione
Al di sotto della saturazione prevista dal modello
Nelle condizioni predefinite, la concentrazione di silice disciolta è inferiore alla concentrazione di saturazione del modello corretta per l’attività.
SOLUBILITÀ TOTALE PREVISTA DAL MODELLO
117.2 mg/kg espressi come SiO₂
CONCENTRAZIONE DELLO SCENARIO
80.0 mg/kg espressi come SiO₂
CONCENTRAZIONE / LIMITE DEL MODELLO
68%
La percentuale è un rapporto tra concentrazioni, non un indice di saturazione basato sull’attività né un limite operativo consentito.
Come il pH modifica la silice disciolta
Un pH più elevato sposta l’equilibrio dall’acido silicico neutro verso gli ioni silicato. Alla saturazione, il riferimento di attività della specie neutra vincola H₄SiO₄; le specie ioniche contribuiscono alla silice disciolta totale. Le percentuali riportate di seguito sono calcolate per l’acqua selezionata, includendo il NaCl di fondo e gli ioni necessari al bilanciamento del pH. I coefficienti di attività dipendono dalla forza ionica, quindi i calcoli per l’acqua selezionata e per il suo stato di saturazione vengono risolti separatamente.
| Specie monomeriche | Frazione calcolata | Silice nell’acqua selezionata · mg/kg espressi come SiO₂ |
|---|---|---|
| H₄SiO₄ · neutra | 99.85% | 79.9 |
| H₃SiO₄⁻ · carica singola | 0.15% | 0.1 |
| H₂SiO₄²⁻ · carica doppia | <0.01% | 0.0 |
pK₁ calcolato ≈ 9.83 a 25 °C. I coefficienti di attività vengono determinati iterativamente.
Il modello include NaCl, silice monomerica disciolta, H⁺, OH⁻ e le quantità minime di Na⁺ o Cl⁻ necessarie al bilanciamento della carica. Gli altri ioni non sono inclusi in questo modello.
Tre indicazioni utili
Il raffreddamento riduce la solubilità di riferimento
Il valore della correlazione in funzione della temperatura aumenta all’aumentare della temperatura. Un flusso che, quando è caldo, si trova al di sotto del riferimento può superarlo durante il raffreddamento.
Un pH più elevato aumenta la solubilità apparente della silice
Il modello calcola una maggiore quantità di ioni silicato a pH più elevati, quindi la silice disciolta totale alla saturazione della silice amorfa aumenta. Un valore di solubilità maggiore non elimina i rischi di incrostazioni da silicati metallici o di altro tipo.
L’effetto del sale dipende dal pH
Il NaCl modifica in modo diverso le attività dell'acido silicico neutro e delle specie di silicato cariche. In questo modello, può ridurre leggermente la solubilità a pH prossimo alla neutralità e aumentarla a pH più elevati. Non esiste una regola universale secondo cui «più sale significa maggiore solubilità».
La silice disciolta, la silice colloidale e la silice totale sono diverse
Questa pagina somma l'acido silicico neutro e le sue due forme monomeriche dissociate e confronta tale totale con una valutazione preliminare dell'equilibrio della silice amorfa. La silice colloidale è costituita da particelle, mentre la silice totale può includere materiale disciolto e particolato, a seconda del metodo analitico. Un risultato di laboratorio indicato come «silice reattiva» dipende dal metodo utilizzato; verificare su quale base è espresso prima di inserirlo.
Il valore della silice espresso come Si deve essere moltiplicato per circa 2.139 per esprimere la stessa quantità come SiO₂. Il riferimento è la silice amorfa, non il quarzo né un minerale di silicato di magnesio.
Dove è utile?
Concentrato da osmosi inversa (RO) e riutilizzo dell'acqua
Utilizzare un valore misurato di silice disciolta per esplorare gli effetti della temperatura e della concentrazione ideale. Una valutazione reale dell'osmosi inversa richiede anche il pH del concentrato, la composizione chimica dei sali, il passaggio della silice e la concentrazione sulla superficie della membrana.
Evaporazione e ZLD
Esplorare l'effetto della concentrazione di una corrente. I liquidi ad alta salinità degli evaporatori richiedono una modellazione geochimica basata sulla composizione completa; questa valutazione preliminare limitata al NaCl non rappresenta le salamoie.
Raffreddamento di correnti geotermiche o di processo
Mostrare perché il raffreddamento può portare la silice disciolta al di sopra di un valore di riferimento della solubilità. La precipitazione effettiva può essere ritardata e deve essere valutata per l'acqua specifica.
I calcoli sopra riportati hanno valore di riferimento preliminare e possono contenere errori od omissioni. Non viene fornita alcuna garanzia di accuratezza o completezza; i risultati devono essere verificati in modo indipendente prima di essere utilizzati nella progettazione, nella scelta del trattamento o nelle decisioni operative. Questo strumento non sostituisce una modellazione geochimica completa o prove pilota nei sistemi ad alto recupero o ad alta temperatura.
Metodologia e limitazioni
Le equazioni, la base di espressione delle concentrazioni e le ipotesi alla base della valutazione preliminare della silice amorfa, comprese la speciazione in funzione del pH, le correzioni di attività per il NaCl e il confronto esistente con PHREEQC.
1. Dati di ingresso, unità e basi dello scenario
| Dato di ingresso | Intervallo dell'interfaccia | Significato in questo modello |
|---|---|---|
| Temperatura | 5–95 °C | Riferimento per acqua liquida in condizioni prossime a quelle atmosferiche; nessuna correzione per la pressione. |
| Silice disciolta | 0–1000 mg/kg di acqua come SiO₂ | Somma delle specie monomeriche disciolte di acido silicico e silicato prima della concentrazione; esclude particelle e colloidi. |
| pH operativo | 4–10 | Imposto indipendentemente a ciascuna temperatura, anziché previsto durante il riscaldamento o la concentrazione. |
| Fattore di concentrazione ideale, F | 1–10× | Moltiplica sia la silice disciolta sia il NaCl di fondo, ipotizzando una ritenzione completa prima della precipitazione. |
| NaCl di fondo | 0–5000 mg/kg di acqua prima della concentrazione | Solo NaCl effettivamente presente; non i TDS né un sostituto degli altri ioni presenti in un'analisi. |
Le concentrazioni sono riferite a chilogrammi di acqua, non a chilogrammi di soluzione. Il codice converte la silice in molalità utilizzando 60084.3 mg/mol e il NaCl utilizzando 58442.77 mg/mol. Non viene applicata alcuna conversione basata sulla densità. Nell'acqua con basse concentrazioni di soluti, a temperatura prossima a quella ambiente, mg/L equivale approssimativamente a mg/kg di acqua; questa approssimazione non deve essere estesa automaticamente all'acqua calda o alle salamoie.
Silice dello scenario = silice inserita × F NaCl dello scenario = NaCl inserito × F
Il pH operativo inserito rimane fisso. Lo scenario non calcola un recupero reale dell'osmosi inversa, un percorso di evaporazione, il passaggio dei sali o la massa precipitata. Gli intervalli dell'interfaccia non costituiscono limiti di accuratezza convalidati in modo indipendente.
Equazioni, modello di attività e soluzione numerica
1. Riferimento di attività delle specie neutre: Fournier–Rowe (1977): log₁₀ C₀ = 4.52 − 731/T, con T in kelvin e C₀ espresso come SiO₂ in mg/kg di acqua. La correlazione per l'acqua pura è utilizzata come riferimento di attività delle specie neutre in condizioni diluite nell'intervallo 5–95 °C. Non viene applicata alcuna correzione per la pressione.
2. Costanti di dissociazione: H₄SiO₄ ⇌ H₃SiO₄⁻ + H⁺ (K₁) e H₄SiO₄ ⇌ H₂SiO₄²⁻ + 2H⁺ (β₂). β₂ è una costante cumulativa, non la costante del secondo stadio di dissociazione. Le espressioni analitiche di PHREEQC utilizzano log K = A₁ + A₂T + A₃/T + A₄ log₁₀T + A₅/T². Coefficienti per K₁: −302.3724, −0.050698, 15669.69, 108.18466, −1119669; per β₂: −294.0184, −0.07265, 11204.49, 108.18466, −1119669; per la dissociazione dell'acqua: −283.971, −0.05069842, 13323, 102.24447, −1119669.
3. Correzioni di attività: per gli ioni silicato, log₁₀ γᵢ = −Azᵢ²[√I/(1 + √I) − 0,3I] (Davies). Per H₄SiO₄ neutro, log₁₀ γ₀ = 0,1I, in accordo con la convenzione per le specie neutre nella versione fissata della banca dati usata per il confronto di riferimento. H⁺ e OH⁻ utilizzano coefficienti di Debye–Hückel estesi con parametri di dimensione ionica pari a 9.0 e 3.5 Å e b-dot nullo. I parametri A e B, dipendenti dalla temperatura, sono interpolati dai valori delle proprietà dell'acqua di PHREEQC a intervalli di 5 °C, nell'intervallo 5–95 °C, alla pressione atmosferica di riferimento. Il NaCl contribuisce con molalità uguali di Na⁺ e Cl⁻; il valore inserito viene diviso per 58442.77 mg/mol. Queste convenzioni sono specifiche della banca dati e non sono coefficienti universali di solubilità dei sali.
4. Equilibrio e chiusura del bilancio di carica: R₁ = 10pH + log K₁; R₂ = 102pH + log β₂. Alla saturazione, C = C₀(1/γ₀ + R₁/γ₁ + R₂/γ₂). Le frazioni sono proporzionali a 1/γ₀, R₁/γ₁ e R₂/γ₂. Per la concentrazione di silice selezionata, le specie vengono ricalcolate alla forza ionica corrispondente. H⁺ e OH⁻ vengono calcolati a partire dal pH imposto e dalla dissociazione dell'acqua. Il modello aggiunge la quantità minima di Na⁺ o Cl⁻ necessaria per l'elettroneutralità, corrispondente a una regolazione idealizzata con NaOH o HCl. La forza ionica I = ½Σmᵢzᵢ² e i coefficienti di attività vengono determinati iterativamente. Si tratta di uno scenario NaCl–silice semplificato con bilancio di carica; non calcola una dose pratica di acido o soda caustica né una composizione dell'acqua misurata.
5. Intervallo di valutazione preliminare: i calcoli relativi sia all'acqua selezionata sia al suo stato di saturazione previsto devono convergere e fornire una forza ionica modellizzata ≤0.10 mol/kg. La concentrazione effettiva di NaCl di fondo dopo la concentrazione deve essere ≤5000 mg/kg di acqua. In caso contrario, il confronto con la saturazione non viene visualizzato; la curva blu è limitata alle condizioni comprese nell'ambito di applicazione del modello. I risultati relativi alle specie non vengono visualizzati quando l'acqua selezionata supera tale ambito. Questi sono limiti operativi per la valutazione preliminare, non limiti di accuratezza convalidati in modo indipendente. Gli altri sali non possono essere rappresentati semplicemente inserendo il loro TDS come NaCl.
6. Semplificazioni: l'attività dell'acqua è impostata a 1, senza un modello osmotico specifico per il sale. L'intercetta arrotondata della relazione di riferimento per la specie neutra differisce leggermente dalla costante di fase della versione fissata di PHREEQC; questa differenza e l'attività dell'acqua non considerata spiegano gran parte dello scostamento residuo nel confronto di riferimento. Il pH operativo viene imposto indipendentemente a ogni temperatura. Il fattore di concentrazione moltiplica sia la concentrazione della silice sia quella di NaCl, ipotizzando una ritenzione completa e una riduzione della massa d'acqua, prima della precipitazione. Non vengono previsti l'alcalinità naturale, gli altri ioni, i complessi, le fasi di silicati metallici, la silice polimerizzata, i colloidi, la cinetica, le prestazioni degli antincrostanti, la concentrazione sulla superficie della membrana e il recupero ammissibile. La solubilità non indica il tempo necessario perché si verifichi il fouling.
Questo calcolatore stima la saturazione della silice amorfa in uno scenario NaCl–silice limitato. La concordanza con PHREEQC a parità di ipotesi non convalida un progetto basato su un'acqua reale. Verificare i risultati rispetto alla composizione chimica completa dell'acqua e, ove opportuno, ai dati sperimentali.
Interpretazione dei risultati e regole di visualizzazione
Rapporto di concentrazione (%) = 100 × silice dello scenario / concentrazione di saturazione modellata
Il rapporto viene arrotondato a una percentuale intera, mentre le etichette di stato utilizzano il valore non arrotondato: al di sotto dell' 85% lo stato è «Al di sotto della saturazione modellizzata»; tra l' 85–100% lo stato è «In prossimità della saturazione modellizzata»; al di sopra del 100% lo stato è «Al di sopra della saturazione modellizzata». La soglia dell' 85% è una scelta di avviso dell'interfaccia, non un margine di sicurezza stabilito su base chimica. Il rapporto non è un indice generale di saturazione minerale basato sull'attività né un limite operativo consentito.
La curva blu rappresenta la saturazione calcolata a ogni temperatura per il pH selezionato e la concentrazione effettiva di NaCl. La speciazione dell'acqua selezionata viene calcolata separatamente alla concentrazione di silice inserita. La curva grigia tratteggiata è il riferimento in funzione della temperatura per la specie neutra in acqua pura, non la curva completa di speciazione in assenza di NaCl. Anche con NaCl di fondo pari a zero, gli ioni di bilanciamento del pH e gli ioni silicato possono contribuire alla forza ionica.
Il procedimento iterativo per la forza ionica consente fino a 160 iterazioni, applica uno smorzamento a mezzo passo e considera raggiunta la convergenza quando la variazione è inferiore a 10⁻¹² + 10⁻⁹ × max(I precedente, I aggiornato). Si arresta anticipatamente se la forza ionica aggiornata non è finita o supera 2 mol/kg. Tale arresto numerico è distinto dal limite più restrittivo di 0.10 mol/kg per la visualizzazione dei risultati.
I punti non supportati dal modello vengono omessi dalla curva blu. Il rapporto e il marcatore selezionato non vengono visualizzati a meno che sia l'acqua selezionata sia il suo stato di saturazione rientrino nell'ambito di applicazione del modello. Le frazioni delle specie non vengono visualizzate se anche solo l'acqua selezionata non rientra in tale ambito. «Fuori dall'ambito di applicazione del modello» non deve essere interpretato come un basso rischio di incrostazione.
Esempio predefinito riproducibile
Il ripristino imposta 25 °C, pH 7.0, silice disciolta 80 mg/kg di acqua come SiO₂, fattore di concentrazione 1× e NaCl di fondo pari a zero. La concentrazione totale di saturazione modellizzata è di circa 117.2 mg/kg come SiO₂; il rapporto di concentrazione visualizzato è del 68%, con lo stato «Al di sotto della saturazione modellizzata». La speciazione dell'acqua selezionata è costituita prevalentemente da H₄SiO₄ neutro. Un valore di NaCl di fondo pari a zero non elimina gli ioni richiesti dal pH e dal bilancio di carica del modello.
Confronto di riferimento e limiti della validazione
Come si confronta questo modello con PHREEQC?
Un confronto di riferimento con acqua sintetica ha confrontato 396 combinazioni di temperatura, pH e NaCl con il motore PHREEQC e la versione fissata del database phreeqc.dat fornita con IPhreeqcPy 2024.0.0. Entrambi i modelli hanno utilizzato silice amorfa, SiO₂(a), con pH operativo fisso. PHREEQC includeva HCl o NaOH per mantenere il pH e il proprio trattamento dell'attività dell'acqua; questa pagina utilizza la quantità minima di ioni per il bilanciamento della carica e assume che l'attività dell'acqua sia 1.
| Condizioni | NaCl · mg/kg di acqua | Questo modello · mg/kg come SiO₂ | PHREEQC · mg/kg come SiO₂ |
|---|---|---|---|
25 °C · pH 7 | 0 | 117.2 | 116.8 |
25 °C · pH 7 | 5000 | 114.9 | 113.9 |
25 °C · pH 9 | 0 | 134.6 | 134.2 |
25 °C · pH 9 | 5000 | 136.6 | 135.4 |
80 °C · pH 9.5 | 5000 | 966.5 | 957.2 |
95 °C · pH 10 | 5000 | 3961.5 · non visualizzato | 3907.0 |
La griglia comprendeva 5–95 °C, pH 4–10 e 0–5000 mg/kg di NaCl, inclusi 18 casi a temperature intermedie. Di questi casi, 391 hanno soddisfatto il criterio di questo modello relativo alla forza ionica allo stato di saturazione; la maggiore differenza relativa rispetto a PHREEQC è stata dello 0.98%. Cinque casi hanno superato il limite e il relativo confronto con la saturazione non viene visualizzato. Queste differenze rilevate nei casi esaminati non costituiscono limiti di accuratezza garantiti. I due calcoli condividono le costanti di equilibrio e le ipotesi del modello di attività, quindi si tratta di un confronto tra implementazioni, non di una validazione sperimentale indipendente.
L'acqua di processo reale può contenere calcio, magnesio, alluminio, ferro e altri costituenti qui assenti. Per la progettazione, utilizzare un'analisi completa dell'acqua, un database e un modello di attività adeguati; né questo confronto né il controllo del NaCl stabiliscono un recupero sicuro o una velocità di fouling.
Riferimenti e documentazione del modello
Riferimenti: Database PHREEQC dell'USGS · Costanti di equilibrio · PHREEQC · Espressioni in funzione della temperatura e modelli di attività · Fournier & Rowe · Solubilità della silice amorfa · DuPont FilmTec · Prevenzione delle incrostazioni di silice
Documentazione della metodologia v1.0 · redatta il 24 ottobre 2026. Documenta il motore di calcolo esistente e il relativo confronto di riferimento riportato; questa scheda non modifica alcuna funzione numerica né alcun coefficiente.