Cas : panache Cr(VI)

Ce cas illustre le couplage entre spéciation redox et transport sur un panache de chrome en nappe. La question pratique à laquelle il répond est la suivante : un site industriel a historiquement rejeté du chrome, une partie a migré dans la nappe, comment prédire son extension aval et identifier les zones où l'atténuation naturelle par réduction redox est effective.

Contexte

Ancienne activité de traitement de surface arrêtée depuis quinze ans. Plusieurs forages de contrôle montrent des concentrations en chrome total entre 200 et 1200 µg/L dans la nappe superficielle, avec un ratio Cr(VI) sur chrome total variable d'un point à l'autre. Le maître d'ouvrage souhaite évaluer si l'atténuation naturelle est suffisante pour considérer le site comme stabilisé, ou si une remédiation active est nécessaire.

Les enjeux typiques de ce type de situation se résument à trois questions : où est le front du panache à l'heure actuelle, comment va-t-il évoluer dans les dix prochaines années, et quelle fraction du chrome migré est déjà réduite en Cr(III) immobile.

Mise en place de la scène

L'emprise choisie couvre 800 mètres dans le sens d'écoulement et 400 mètres transversalement, centrée sur la zone source supposée. La profondeur va de la surface à 15 mètres, incluant la nappe superficielle et la couche argileuse qui la sépare de l'aquifère profond. Sept forages d'observation fournissent des points d'échantillonnage, dont trois dans la zone source et quatre en aval.

Élément de scène	Configuration
Emprise horizontale	800 × 400 m, orientée selon l'axe d'écoulement
Profondeur	15 m, référence au niveau du sol
Grille de calcul	80 × 40 × 30 = 96 000 mailles (résolution intermédiaire)
Couche superficielle (0 à 2 m)	Limon sableux, porosité 0,20, K = 5e-5 m/s
Aquifère superficiel (2 à 12 m)	Sable grossier, porosité 0,25, K = 5e-4 m/s
Écran argileux (12 à 15 m)	Argile, porosité 0,04, K = 1e-9 m/s
Points d'échantillonnage	7 forages, tous équipés de piézomètres multi-niveaux

Paramètres du moteur

Le chrome est simulé sous ses deux formes : Cr(VI) mobile avec un Kd faible dans le sable (0,5 L/kg) et Cr(III) pratiquement immobile (Kd 200 L/kg). La transformation Cr(VI) vers Cr(III) est modélisée comme une réaction cinétique dont le taux dépend du Eh local et de la présence de donneurs d'électrons (fer ferreux, matière organique). Les mesures de Eh sur les sept forages donnent une gamme allant de +250 mV à -50 mV, suggérant une hétérogénéité redox significative.

Paramètre	Valeur	Origine
Source	Zone 50 × 30 m, profondeur 2 à 6 m, flux historique estimé	Bilan d'activité de la zone atelier
Gradient hydraulique	0,004 m/m, résolu par piézométrie en 7 points	Mesures des forages
Vitesse effective moyenne dans sable	environ 25 m/an	Darcy avec porosité 0,25
Kd Cr(VI)	0,5 L/kg	Littérature, nappe sableuse peu adsorbante
Kd Cr(III)	200 L/kg	Littérature, précipitation en hydroxyde
Taux de réduction Cr(VI) vers Cr(III)	variable 0,001 à 0,05 /jour	Fonction de Eh local, calé sur ratios Cr(VI)/Cr total
Durée simulée	20 ans (15 ans passés + 5 ans prospectifs)	Couvre la période depuis la fin d'activité

Le point critique du paramétrage est le taux de réduction Cr(VI) vers Cr(III). Il varie de plusieurs ordres de grandeur selon le Eh, la matière organique et la biomasse microbienne. Le caler sur les ratios Cr(VI)/Cr total mesurés aux forages transforme ce paramètre flou en paramètre contraint par les données.

Résultats

Après calcul de 20 ans de simulation, le champ de concentration Cr total présente trois zones distinctes. La zone source accumule encore une masse résiduelle importante, largement sous forme Cr(III) précipité. Une zone intermédiaire de quelques centaines de mètres en aval est dominée par du Cr(VI) en transit avec des signes de réduction partielle. Une zone aval plus lointaine ne montre que des traces de chrome, à des concentrations proches du bruit analytique.

La coupe verticale le long de l'axe d'écoulement révèle une stratification verticale nette : le chrome circule préférentiellement dans les niveaux sableux les plus perméables, et s'accumule en bas de la nappe juste au-dessus de l'écran argileux. La signature isotopique δ⁵³Cr suit cette structure : les valeurs les plus enrichies (signe d'une réduction poussée) sont au fond de la nappe dans la zone intermédiaire, là où le Eh mesuré est le plus bas.

Zone	Concentration Cr total simulée	Fraction Cr(VI)	δ⁵³Cr simulé
Source (0 à 100 m)	500 à 1500 µg/L	20 à 40 %	+0,5 à +1,2 ‰
Intermédiaire (100 à 400 m)	100 à 500 µg/L	40 à 70 %	+1,5 à +3,5 ‰
Front aval (400 à 700 m)	5 à 50 µg/L	60 à 90 %	+0,8 à +2,5 ‰

Interprétation

Le modèle reproduit globalement les mesures aux sept forages : concentrations et ratios Cr(VI)/Cr total corrects dans la zone source et en aval proche, concentrations prédites inférieures au seuil de détection dans la zone aval lointaine en cohérence avec les mesures. La signature isotopique δ⁵³Cr simulée aux trois points où elle a été mesurée est dans l'intervalle de confiance des mesures, ce qui valide le paramétrage de la réduction.

La projection à cinq ans supplémentaires montre une stabilisation globale du panache. La zone source continue de fournir du Cr(VI) par désorption lente, mais la masse mobilisée est largement compensée par la réduction en Cr(III) dans la zone intermédiaire. Le front aval n'avance plus significativement, et les concentrations résiduelles continuent de décroître lentement.

Cette interprétation doit être nuancée par plusieurs incertitudes identifiées. D'abord, le taux de réduction est calé sur une situation actuelle ; si les conditions redox évoluent (assèchement, modification d'usage du site en amont), la dynamique pourrait changer. Ensuite, le stock source reste mal quantifié et une désorption plus lente que prévu pourrait prolonger l'activité du panache au-delà de ce que la simulation prédit. Enfin, la représentation en grille régulière ne capte pas les éventuels chemins préférentiels liés à des hétérogénéités locales.

L'interprétation « atténuation naturelle suffisante » que suggère la simulation reste conditionnelle aux hypothèses de paramétrage. Une remédiation active (barrière ZVI, injection de réducteurs) peut accélérer la réduction d'un facteur 5 à 10, et reste une option à considérer si la stabilité prédite est jugée insuffisamment rapide par rapport aux enjeux du site.

Prolongements

Ce cas peut être étendu de plusieurs manières utiles pour approfondir la compréhension du site ou pour tester des actions.

Variante « sans réduction » : désactiver la réaction Cr(VI) vers Cr(III) pour voir à quoi ressemblerait le panache en l'absence d'atténuation naturelle. La comparaison quantifie l'effet réel de la réduction sur l'extension du panache.
Variante « barrière réactive » : ajouter une zone de forte capacité réductrice perpendiculaire à l'écoulement pour simuler l'effet d'une installation ZVI. La comparaison oriente le dimensionnement (longueur, épaisseur, position).
Variante « climat sec » : diminuer la recharge météorique de 30 %, ce qui ralentit l'écoulement. La comparaison montre si le panache est contrôlé par le transport (accélération bénéfique) ou par la cinétique de réduction (ralentissement neutre ou néfaste).
Variante « source persistante » : allonger la durée d'apport actif à la source pour tester l'hypothèse d'un relargage plus long que prévu.

Ces variantes peuvent être enregistrées comme scénarios distincts dans le module de comparaison d'hypothèses et confrontées directement sur la même scène.

Pour aller plus loin

Spéciation et propagation : détails du couplage utilisé dans ce cas.
Spéciation redox : Cr(VI) contre Cr(III) en géochimie générale.
Multi-hypothèses : comparer les variantes évoquées en prolongement.
Interpréter les résultats : lecture générale des sorties de simulation.