Beaucoup d'études ont montré que l'équation d'advection-dispersion classique ne permet pas de simuler correctement le transport de solutés dans les sols hétérogènes, ni de prendre en considération la spéciation des solutés dans les systèmes géochimiques que constituent les sols. Dans ce travail, un modèle fractionnaire hydrogéochimique a été proposé pour simuler le transport et la spéciation des métaux lourds dans la zone non saturée des sols, que ce soit en régime permanent ou transitoire. Ce modèle a été proposé pour remédier aux limitations du modèle classique d'advection dispersion. En régime permanent, la solution analytique de l'équation fractionnaire d'advection-dispersion a été couplée sous MATLAB au modèle de réactions géochimiques, et ce nouveau modèle a été validé à l'aide de résultats expérimentaux. En régime non permanent, une nouvelle solution numérique de l'équation fractionnaire d'advection-dispersion est proposée, et couplée avec un modèle d'écoulement et un modèle géochimique. Le modèle résultant, programmé sous MATLAB, a été testé en le comparant à des simulations obtenues avec les codes HYDRUS-1D et HP1. Les résultats de validation ont montré que le nouveau modèle fractionnaire reproduit bien le transfert de solutés dans la zone non saturée des sols et qu'il est capable de donner plus de détails sur les espèces chimiques présentes dans le sol, sur leur migration et leur interaction. Le nouveau modèle a été utilisé pour étudier le transfert de zinc dans la région de Kempen (à la frontière entre la Belgique et les Pays Bas). Il s'agit d'un site fortement pollué par les métaux lourds rejetés par les fonderies de zinc existant dans la région. Une étude paramétrique a été conduite pour déterminer la sensibilité du modèle à une variation de ses paramètres hydrologiques ou géochimiques. La conductivité hydraulique du sol (Ks), la teneur en eau à saturation (θs) et la teneur en eau initiale (θini) du sol sont les paramètres les plus influents pour le modèle d'écoulement d'eau. Le modèle fractionnaire de transport de soluté est sensible à la variation de l'ordre fractionnaire de dérivation (α) et à celle du coefficient de dispersion (D). Le pH est le facteur déterminant pour le modèle géochimique, suivi par la concentration en SO42- et en CO32-. L'effet des cations Al3+, Mn2+ et Fe2+ n'est pas significatif.
Identifer | oai:union.ndltd.org:CCSD/oai:tel.archives-ouvertes.fr:tel-00422621 |
Date | 13 July 2007 |
Creators | Dawood, Ihssan |
Publisher | Université des Sciences et Technologie de Lille - Lille I |
Source Sets | CCSD theses-EN-ligne, France |
Language | English |
Detected Language | French |
Type | PhD thesis |
Page generated in 0.002 seconds