Pollution des sols et interaction avec la biosphère

Pollution des sols et interaction avec la biosphère Mazen SAAD Ecole Centrale de Nantes Laboratoire de Mathématiques Jean Leray FRANCE Mathématiques Pour La Planète La Rochelle 17 octobre 2013

Plan 1 Transport d éléments radioactifs 2 Interaction avec la biosphère 3 Stockage du CO2 Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 2 / 31

Transport d éléments radioactifs Stockage des déchets radioactifs GdR MoMaS initié par l ANDRA Simuler la fuite éventuelle des radionucléides Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 3 / 31

Transport d éléments radioactifs Stockage des déchets radioactifs Plan du site de stockage Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 4 / 31

Transport d éléments radioactifs Modélisation de transport des éléments radioactifs Simuler une fuite éventuelle des radionucléides de leur site de stockage Sol : populations Sous sol : radionucleides Site de stockage Une seule phase EAU dans un milieu poreux hétérogène Des composés chimiques miscibles à l eau. Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 5 / 31

Transport d éléments radioactifs Modèle mathématique Déplacement de m-radionuclides de concentration c = (c 1,..., c m) T ( ) ( tg k x, ck + div ck V ) ( ( div D k x, ck, V ) ) c k (t, x) ( ) ( ) + λ k G k x, ck λ l R k,l G l x, cl = fk, k = 1,..., m, (1) l<k div V = 0, et V = 1 ( ) µ(c) K p ρ(c) g avec G k (x, c) = φ(x)c + ( 1 φ(x) ) ρ s(x)f k (x, c) Adsorption : F k (x, c) Transport : c k V Diffusion-Dispersion : D k = d mi d + d s(x, V ) Décroissance radioactive : λ k Filiation : R k,l (structure triangulaire) f k : terme source (2) Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 6 / 31

Transport d éléments radioactifs Stockage des déchets radioactifs Site de Bure Site de stockage : 3240 6 m 2 Plusieurs strates : argile, calcaire 2 composés radioactifs : Iode 129 I, Plutonium 242 Pu Adsorption dans l argile : 0, 001 pour 129 I 2 10 5 pour 242 Pu. Décroissance : T 1 ( 129 I )= 1, 57 10 7 années; T 2 ( 242 Pu)= 3, 76 10 5 années. Ecoulement de droite à gauche Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 7 / 31

Fuite de radioactifs Transport d éléments radioactifs 0.0016 0.0014 0.0012 0.001 0.0008 0.0006 0.0004 0.0002 0 1000 10000 100000 0.0001 9e-05 8e-05 7e-05 6e-05 5e-05 4e-05 3e-05 2e-05 1e-05 0 1000 10000 100000 1e+06 Termes sources pour 129 I (gauche) et 242 Pu (droite) Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 8 / 31

Transport d éléments radioactifs Domaine Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 9 / 31

Transport d éléments radioactifs Maillage Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 10 / 31

Transport d éléments radioactifs Iode 129 au Temps = 10110 années Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 11 / 31

Transport d éléments radioactifs Iode 129 au Temps = 100000 années Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 12 / 31

Transport d éléments radioactifs Iode 129 au Temps = 250000 années Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 13 / 31

Transport d éléments radioactifs Iode 129 au Temps = 10 6 années Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 14 / 31

Transport d éléments radioactifs Plutonium 242 au Temps = 10 6 années Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 15 / 31

Transport d éléments radioactifs Plutonium 242 au Temps = 5 10 6 années Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 16 / 31

Interaction avec la biosphère Interaction avec la biosphère Prévoir l impact sur une population d une exposition à des radionucléides Contamination d une chaîne alimentaire par des radionucléides Couplage entre milieu poreux et dynamique des populations Sol : populations Sous sol : radionucleides Site de stockage Trois espèces : Végétaux, Animaux, Super-prédateur(Humain?) Exposition des Végétaux aux radionucléides par le sol Contamination des animaux sains par prédation des végétaux exposés Contamination des supers prédateurs sains par contact avec les animaux exposés. Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 17 / 31

Interaction avec la biosphère Modèles de Contamination et de prédation Trois espèces : V : Végétaux ; V = V S + V E A : Animaux ; A = A S + A E P : Super-prédateur ; P = P S + P E Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 18 / 31

Interaction avec la biosphère Modèles de Contamination et de prédation Trois espèces : V : Végétaux ; V = V S + V E A : Animaux ; A = A S + A E P : Super-prédateur ; P = P S + P E P S Population Saine ; P E Population Exposée tp S d P S = r P ( 1 P K P )P S β P A E P S tp E d P E = r P ( 1 P K P )P E + β P A E P S β P : taux de contamination des super prédateurs par consommation des animaux infectés Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 18 / 31

Interaction avec la biosphère Modèles de Contamination et de prédation A S Animaux Sains ; A E Animaux Exposés pv S pv E ta S div(d A S ) = m A A S + e A S + (1 β A )e A S δa S P S 1 + qv S 1 + qv E pv S pv E pv E ta E div(d A E ) = m A A E + e A E + e A E + β A e A S δa E P S 1 + qv S 1 + qv E 1 + qv E β A : taux de contamination des animaux par consommation des végétaux Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 19 / 31

Interaction avec la biosphère Modèles de Contamination et de prédation V S Végétaux Sains ; A E Végétaux Exposés ( tv S = r V 1 V ) V S pv S A ρ(c)v S K V 1 + qv S ( tv E = r V 1 V ) V E pv E A+ρ(c)V S K V 1 + qv E ρ(c) : taux de transmission de la concentration des radionulcéides (présents en surface) aux végétaux. Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 20 / 31

Interaction avec la biosphère Interaction radioactifs Végétaux/Animaux Contamination des végétaux par de l iode radioactive 129 I. La trace de la concentration à la surface est moyennée sur 15m de profondeur. Exercice Couplex 1 apres 800000 ans Sol : populations pendant 300 ans Marne y Calcaire 20 km A & V Argilite Calcaire 300 m z Zone d exposition 25 km y 695 m 40 km x Densités d individus exposés d animaux A E : Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 21 / 31

Interaction avec la biosphère Quantités totales en fonction du temps pour l espèce A (gauche) et l espèce V (droite) 45000 300 40000 35000 30000 250 200 25000 20000 15000 150 100 10000 5000 0 0 50 100 150 200 250 300 50 0 0 50 100 150 200 250 300 La ligne solide représente les individus sains, la ligne pointillée la classe des exposés. Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 22 / 31

Interaction avec la biosphère Radioactifs Végétaux/Animaux/Humains Intrusion d un super-prédateur dans une zone isolée et non infectée. 20 km Animaux + Vegetaux 01 00000 11111 00000 11111 0000000 1111111 00000 11111 11100000 11111 01 Zone du super predateur (K=1000) Zone d exposition 40 km Contamination des humains par déplacement des animaux infectés Densités d individus exposés P E du super-prédateur P, après 100,150, 200, 300 ans de contamination Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 23 / 31

Stockage du CO2 Stockage du CO2 CO2 capture and storage. The process consists of capturing waste CO2 and transporting it to a storage site. Various forms have been conceived for storage CO2 into deep geological formations : CO2 is sometimes injected into declining oil fields to increase oil recovery. This option is attractive because the geology of hydrocarbon reservoirs is generally well understood and storage costs may be partly offset by the sale of additional oil that is recovered. Model 1 : Gas 90% CO2 = Oil Gas model. two compressible flow Model 2 : Dissolution of CO2 in water. Two compressible and partially miscible flow Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 24 / 31

Stockage du CO2 Two-phase two-component flow in porous media TWO PHASES TWO COMPONENTS POROUS MEDIA MODEL Discretize a general model for two phase partially miscible flow in porous media : CO2 storage in a liquid gas reservoir. Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 25 / 31

Stockage du CO2 CO2 in liquid and gas form After the CO2 injection, several different trapping mechanisms lead to an entrapment of the CO2. Shortly after the injection, structural trapping through caprocks is the most important factor. Later solubility trapping, where CO2 is dissolved into water, and residual trapping get more important. After several thousand years, there could also occur mineral trapping caused by geochemical reactions To simulate the process of dissolution of CO2, a multiphase flow equation with equilibrium phase exchange is used. The CO2 storage can be modeled with two components (CO2 and water) in two phases (liquid and gas). Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 26 / 31

Mathematical model Stockage du CO2 We consider herein a porous medium saturated with a fluid composed of : two phases : liquid (α = l) and gas (α = g) two components in each phase : CO2 (β = c) and water (β = w) The component CO2 is present in the two phases : In liquid form : dissolved CO2 In gas form : CO2 in the gas phase The component water exists only in liquid form (no vapor water) Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 27 / 31

Stockage du CO2 Two-phase two-component model: Notation s l, s g : liquid and gas phase saturations p l, p g : liquid and gas phase pressures Phases occupy the entire pore space: s l + s g = 1 ρ l, ρ g : liquid and gas phase densities ρ c l, ρc g : density of dissolved CO2 and density of the CO2 in the gas phase Xl c = ρ c l /ρ l : mass fraction of the CO2 in the liquid phase X w l = ρ w l /ρ l : mass fraction of the water in the liquid phase Xl c + Xl w = 1 Xg c = ρc l /ρ l : mass fraction of the CO2 in the gas phase Capillary pressure between liquid pressure and gas pressure : p c(s l ) = p g p l Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 28 / 31

Stockage du CO2 Two-phase two-component model: Assumption Isothermal flow (T is a constant) Perfect gas law for CO2 gas phase ρ c g = Mc RT pg where M c is the molar mass of the CO2 and R is the perfect gas constant Equilibrium between dissolved CO2 and gas CO2 is prescribed by Henry s law, ρ c l = M c H c p g where H c is the Henry s constant. Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 29 / 31

Mathematical model Stockage du CO2 The mass conservation of each component : CO2: Φ t(ρ c l s l + ρ c g s g ) + div(ρ c l V l + ρ c g Vg + Jc l + J c g ) = Qc Water: Φ t(ρ w l s l ) + div(ρ w l V l + J w l ) = Q w where diffusion of the components β = c, w in the phase α = l, g is modeled by J β α = ΦsαDα ρβ α, J c α + Jw α = 0 Darcy s law for the velocity of each phase V α = K krα(sα) µ α ( p α ρ αg) Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 30 / 31

CO2 field Stockage du CO2 CO2 in the dissolved form. Mazen Saad (ECN) Stockage et interaction La Rochelle 17 octobre 2013 31 / 31