vendredi, 24 avril, 2026
Spéciation et propagation
La spéciation, la cinétique et le transport ne sont pas des processus indépendants : un contaminant ne se déplace pas de la même façon selon qu'il est sous forme mobile ou adsorbée, et il ne se dégrade pas à la même vitesse selon les conditions redox locales. Cette page détaille comment IsoFind couple ces trois dimensions dans le moteur de simulation, avec un accent sur les situations où ce couplage change qualitativement le résultat prédit.
Le problème du couplage
Un modèle de transport pur suppose que le contaminant est chimiquement inerte : sa seule évolution est son déplacement. Un modèle réactif sans transport suppose que le contaminant reste en place et ne fait qu'évoluer chimiquement. La réalité combine les deux, et le couplage produit des comportements qu'aucun des deux modèles simples ne capture.
L'exemple le plus démonstratif est le panache de chrome en nappe. Le Cr(VI) est soluble et mobile, il se déplace avec l'écoulement. Le Cr(III) est peu soluble et fortement adsorbé, il reste quasiment en place. Un même stock de chrome peut donc produire deux panaches radicalement différents selon les conditions redox locales. La transformation Cr(VI) vers Cr(III) par réduction microbienne ou abiotique est donc un mécanisme d'atténuation massive du panache, à condition qu'elle soit effective dans les zones traversées.
La capacité à coupler spéciation, cinétique et transport est ce qui rend la Visualisation 3D d'IsoFind différente d'un outil de simulation classique. Les paramètres locaux (pH, Eh, matière organique) modulent les cinétiques maille par maille, pas seulement à l'échelle globale.
Spéciation redox maille par maille
Pour les composés dont le comportement dépend de l'état redox, le moteur calcule à chaque pas de temps la répartition entre espèces dans chaque maille à partir du pH et du Eh locaux. Cette répartition est dérivée de diagrammes Eh-pH tabulés dans la base de référence, ajustés en fonction des concentrations des espèces complexantes présentes.
| Élément | Espèces dominantes | Contrôle principal |
|---|---|---|
| Chrome | Cr(VI) mobile, Cr(III) adsorbé | Eh surtout, pH secondaire |
| Fer | Fe(II) mobile, Fe(III) précipité en oxyhydroxydes | Eh, pH, alcalinité |
| Arsenic | As(V) oxyanion adsorbé, As(III) moins adsorbé et mobile | Eh, présence d'oxyhydroxydes de Fe |
| Uranium | U(VI) uranyle soluble, U(IV) uraninite insoluble | Eh, complexation par carbonates |
| Manganèse | Mn(II) soluble, Mn(IV) oxydes solides | Eh, pH, activité microbienne |
| Sélénium | Se(VI) très mobile, Se(IV) moins mobile, Se(0) insoluble | Eh, présence de Fe(II) |
Couplage cinétique
Pour les composés organiques, la constante cinétique de dégradation k n'est pas uniforme dans le domaine : elle dépend de la voie dominante localement, elle-même fonction des conditions géochimiques. Le moteur calcule donc un k effectif par maille, à partir de plusieurs sources d'information.
| Source d'information | Contribution au calcul local de k |
|---|---|
| Base de référence moléculaire | Liste des voies possibles, k de référence, conditions favorables |
| Conditions géochimiques locales | pH, Eh, O₂ dissous, nitrates, sulfates, fer, méthane interpolés |
| Température locale | Correction d'Arrhenius appliquée au k de référence |
| Matière organique locale | Proxy de la biomasse microbienne potentielle |
Le résultat est un champ de k variable dans l'espace qui peut passer du simple au décuple entre une zone aérobie riche en carbone organique et une zone anaérobie oligotrophe. Cette variabilité est ce qui génère les formes de panache observées sur le terrain : zones de dégradation rapide en périphérie aérobie, zone de persistance centrale plus réductrice.
Override du bridge Nexus
Pour piloter explicitement le choix de la voie ou du mode de fractionnement, la requête de simulation accepte plusieurs champs qui court-circuitent l'auto-sélection du bridge. Ces overrides sont utiles quand le contexte du site contredit ce que le bridge déduirait automatiquement, ou pour tester des hypothèses alternatives.
| Champ de la requête | Rôle | Valeurs typiques |
|---|---|---|
| nexus_process_hint | Force la catégorie de processus à interroger en base de fractionnement | redox, adsorption, biological, dissolution, precipitation, evaporation |
| prefer_fractionation_mode | Force le mode de fractionnement | equilibrium (Rayleigh) ou kinetic |
| oxidation_state_initial et _final | Affine la requête pour les éléments à spéciation | Ex. Cr(VI) → Cr(III) |
| csia_pathway | Voie de dégradation explicite pour une molécule (écrase l'auto-sélection) | Nom exact depuis molecule_degradation_pathways.pathway_name |
| csia_element | Élément isotopique privilégié pour le suivi | C (défaut), N, Cl, H |
Sans ces overrides, le bridge sélectionne automatiquement la voie dominante selon les conditions interpolées sur chaque maille. Les overrides sont donc utiles surtout pour documenter explicitement une hypothèse dans un rapport, ou pour comparer deux scénarios qui diffèrent uniquement par la voie invoquée.
Chaînes parent-métabolite dans la simulation
Lorsque la base de référence déclare une chaîne parent-métabolite pour le composé simulé, le moteur propage simultanément le parent et ses métabolites. Le couplage fonctionne en cascade : la dégradation du parent au taux k_parent produit du métabolite, qui à son tour se dégrade au taux k_métabolite. Pour les solvants chlorés, cette cascade peut comporter quatre étages (PCE vers TCE vers cis-DCE vers VC vers éthylène) simulés en parallèle.
Chaque métabolite a ses propres paramètres de transport : Kd différent, cinétique propre, signature isotopique distincte. La simulation produit donc un empilement de cartes 3D de concentration, une par composé de la chaîne, qui permet de visualiser où chaque composé domine spatialement.
Pour lire un panache en cascade, le mode de visualisation le plus informatif est la superposition semi-transparente des isosurfaces du parent et des principaux métabolites. On voit alors immédiatement si le parent et ses produits coexistent (dégradation active) ou s'ils sont ségrégés dans l'espace (front de dégradation clair entre zone source et zone périphérique).
Propagation de la signature isotopique
La signature isotopique est propagée en parallèle de la concentration. Dans chaque maille où la dégradation se produit, l'enrichissement isotopique est calculé selon l'équation de Rayleigh avec l'epsilon de la voie dominante locale. Le champ final est donc une carte 3D des δ attendus, directement comparable aux mesures CSIA mesurées sur les échantillons du projet.
Cette comparaison est l'un des usages les plus décisifs de la simulation : un modèle qui reproduit quantitativement non seulement les concentrations mais aussi les signatures isotopiques résiduelles est beaucoup plus robuste qu'un modèle qui n'ajuste que les concentrations. Deux modèles radicalement différents peuvent produire les mêmes concentrations, mais ils produiront rarement les mêmes signatures.
Adsorption non linéaire et effets de mémoire
Dans les cas où la source de contamination a été active pendant longtemps puis s'est arrêtée, la matrice solide peut avoir accumulé une quantité significative de contaminant adsorbé. Cette masse se désorbe ensuite progressivement, alimentant le panache sur des durées potentiellement longues après l'arrêt de la source. IsoFind représente ce phénomène avec une isotherme de Langmuir optionnelle, plus réaliste que le modèle Kd linéaire pour les fortes concentrations.
q = (q_max × K × C) / (1 + K × C)
Avec q la masse adsorbée par masse de solide, q_max la capacité maximale de la matrice, K une constante d'affinité et C la concentration en solution. Pour C très petit devant 1/K, l'isotherme de Langmuir se ramène à un Kd linéaire ; pour C grand, elle sature. Cette saturation explique pourquoi les fortes concentrations migrent plus vite que ce qu'un modèle Kd classique prédit.
Précipitation et coprécipitation
Pour certains éléments, le mécanisme dominant n'est pas l'adsorption mais la précipitation sous forme d'une phase solide distincte. Le chrome réduit précipite en oxyhydroxyde Cr(OH)₃, le fer ferreux s'oxyde en oxyhydroxyde de Fe(III), l'uranium réduit précipite en uraninite. Ces précipitations sont traitées comme des puits de masse dans le transport, avec un taux contrôlé par la sursaturation locale.
La coprécipitation, où un élément est entraîné par la précipitation d'une autre phase, est également modélisée pour les cas courants : As coprécipité avec les oxyhydroxydes de Fe, métaux traces coprécipités avec CaCO₃, etc. Ces mécanismes expliquent pourquoi les panaches métalliques s'arrêtent souvent à des interfaces redox marquées sans nécessité d'une dégradation biologique.
Limites du couplage
- Le couplage spéciation-transport suppose l'équilibre thermodynamique local à chaque pas de temps. Pour les réactions lentes, cette hypothèse peut être fausse et le moteur propose alors une cinétique de rééquilibrage avec un taux configurable.
- La représentation par base de référence moléculaire suppose que les voies tabulées couvrent les mécanismes réels du site. Une voie non répertoriée (souche microbienne exotique, condition particulière) n'est pas prise en compte.
- Les hétérogénéités à petite échelle (fissures, micropores, agrégats) qui génèrent des comportements non fickéens ne sont pas résolues par le schéma sur grille régulière.
- Les interactions entre plusieurs contaminants concomitants (compétition d'adsorption, inhibition cinétique mutuelle) ne sont modélisées que dans les cas où elles sont explicitement déclarées.
Le couplage produit des résultats impressionnants, mais l'incertitude sur les paramètres d'entrée (surtout les k de dégradation et les constantes de spéciation) reste souvent l'élément limitant. Présenter les résultats avec une analyse de sensibilité est plus honnête que de présenter une simulation unique comme une vérité.
Pour aller plus loin
- Principe du moteur : équations de base sans couplage.
- Voies de dégradation : conditions favorables à chaque voie.
- Spéciation redox : cadre général pour les équilibres Cr, Fe, As, U.
- Cas : panache Cr(VI) : application concrète du couplage.