Minéralogie, porosité et diffusion des solutés dans l'argilite du ...
Minéralogie, porosité et diffusion des solutés dans l'argilite du ... Minéralogie, porosité et diffusion des solutés dans l'argilite du ...
195 CHAPITRE 5 ⎡ D ⎤ ⎛ ⎞ eα D ⎥ × ⎜ eαt C = ⎢ ⎟ r ( t) C0 exp t erfc 2 ⎢ ⎥ ⎜ ⎟ [5.14] ⎣ H f ⎦ ⎝ H f ⎠ Avec Cr(t) la concentration en traceur dans le réservoir source. La forme de cette équation montre que la décroissance en traceur dépend uniquement du couple De×α. 3.2.2 Sensibilité des paramètres De et α L’expression mathématique du rapport M(x,t) / M(0,t) (eq. 5.13) montre une indépendance des paramètres De et α contrairement à l’expression de la décroissance dans le réservoir (eq. 5.14). Cependant cette expression semble indiquer une dépendance forte des paramètres De et α. Pour savoir si l’on pourra identifier séparément De et α, nous avons réalisé une étude de la sensibilité des paramètres De et α (eq. 5.13) a partir d’un critère de moindres carré. Ce critère est défini par la fonction FO: FO ∑ = L ( M ( x, t) − M ( x, t) ) 1 M 1 2 ( x, t) 2 [5.14] avec M1(x,t) / M(0,t) et M2(x,t) / M(0,t) calculées respectivement pour les paramètres De1, α1 et De2 α2. Cette fonction est définie comme la fonction objective (c. à d. la fonction à minimiser) pour la résolution des problèmes inverses (Sun, 1994). La sensibilité des paramètres est définie par la dérivée de FO vis-à-vis des paramètres variables. Pour illustrer cette approche et estimer la sensibilité des paramètres De et α, FO a été calculée pour De1 = 5×10 -10 m²/s et α1 = 100000 et diverses valeurs de De2 et α2. Sa représentation 2D montre qu’une unique solution satisfait FO = 0, cependant une solution acceptable (FO≈0) est obtenue lorsque De1/α1 = De2/α2 ou De1/α1 = Da1 (Fig. 5.30). L’équation 5.13 montre ainsi une dépendance fortes des paramètres De et α, mais il reste possible de séparer De et α. A présent il convient de savoir dans quelle mesure cette séparation sera effective. D’un point de vue général, l’ajustement d’un profil expérimental par une solution analytique consiste à minimiser la fonction FO (Sun, 1994). Une large variété d’algorithmes de minimisation existe dans la littérature selon la complexité du problème. A partir d’une minimisation par un algorithme de type Gauss-Newton (sous Excel®), la solution exacte (De = 5×10 -10 m²/s et α = 100000) est trouvée seulement si chaque valeur du profil présente une
196 CHAPITRE 5 déviation de moins de 0,0001 %. Au-delà de 0,0001 %, il est uniquement possible d’obtenir une solution de type De/α = Da. En considérant logiquement une erreur expérimentale supérieure à 0,0001 % sur chaque point d’analyse, l’ajustement des profils microLIBS et microsonde par l’équation 5.13 permettra uniquement d’obtenir un Da. A partir de l’ajustement analytique de la concentration dans le réservoir, les paramètres De et α pourront cependant être dissociés. α Fig. 5.30 : Evolution de la fonction FO (eq. 5.14) en fonction de De2 et α2 pour De1= 5×10 - 10 et α1=100000). La ligne pointillée en blanc représente les coordonnées De1/α1=De2/α2. 3.2.3 La méthode Time Domain Diffusion La méthode Time Domain Diffusion (cf. chapitre 3) a fait l’objet d’étude afin de pouvoir interpréter les expériences d’in-diffusion analysées à partir de cartographies élémentaires par une approche de modélisation inverse. Néanmoins, à ce jour, celle-ci ne peut être utilisée dans ce but. Interprétant une expérience d’out-diffusion avec une carte de porosité calculée par autoradiographie, Sardini et al. (2007) combinent la méthode TDD avec une procédure d’inversion afin d’identifier les paramètres de diffusion de différents ensembles poreux. La courbe de restitution en traceur est alors utilisée comme donnée à ajuster. Pour pouvoir reproduire cette approche aux cartes de traceurs obtenues par in-diffusion, la méthode TDD doit considérer : D e (×1.10 -11 m²/s) Log (FO) 3
- Page 146 and 147: 145 CHAPITRE 4 (a) (b) Fig. 4.20: L
- Page 148 and 149: 147 CHAPITRE 4 Fréquence fréquenc
- Page 150 and 151: 149 CHAPITRE 4 Afin de mieux appré
- Page 152 and 153: 151 CHAPITRE 4 • Distribution de
- Page 154 and 155: 153 CHAPITRE 4 L’ensemble des par
- Page 156 and 157: 155 CHAPITRE 4 Calcimétrie, CEC
- Page 158 and 159: 157 CHAPITRE 4 terme de rapport C/I
- Page 160 and 161: 159 CHAPITRE 4 accompagné par une
- Page 162 and 163: 161 CHAPITRE 5 CHAPITRE 5 RELATIO
- Page 164 and 165: 163 CHAPITRE 5 a montré la faisabi
- Page 166 and 167: 165 CHAPITRE 5 2 ∂Ci ∂S i ∂ C
- Page 168 and 169: 167 CHAPITRE 5 la section 2.1.1 de
- Page 170 and 171: 169 CHAPITRE 5 1.5 La cartographie
- Page 172 and 173: 171 CHAPITRE 5 Une concentration in
- Page 174 and 175: 173 CHAPITRE 5 Fig. 5.5 : Visualisa
- Page 176 and 177: 175 CHAPITRE 5 carte de C u apparai
- Page 178 and 179: 177 CHAPITRE 5 Fe-Ca, par D. Prêt.
- Page 180 and 181: 179 CHAPITRE 5 Fig. 5.9 : Procédur
- Page 182 and 183: 181 CHAPITRE 5 2.2.2 La microLIBS
- Page 184 and 185: 183 CHAPITRE 5 Ca Si matrice argile
- Page 186 and 187: 185 CHAPITRE 5 Ca (a) (b) Si Fig. 5
- Page 188 and 189: 187 CHAPITRE 5 (a) (b) Fig. 5.19 :
- Page 190 and 191: 189 CHAPITRE 5 (a) (b) Fig. 5.21 :
- Page 192 and 193: 191 CHAPITRE 5 7500 (680) ppm 12400
- Page 194 and 195: 193 CHAPITRE 5 (a) (b) Fig. 5.28 :
- Page 198 and 199: 197 CHAPITRE 5 - l’inversion de l
- Page 200 and 201: 199 CHAPITRE 5 M/M0 [-] Solution Ec
- Page 202 and 203: 201 CHAPITRE 5 3.3.3 Profils de Cu
- Page 204 and 205: 203 CHAPITRE 5 M/M0 [-] M / M0 1 0,
- Page 206 and 207: 205 CHAPITRE 5 l’échantillon dev
- Page 208 and 209: 207 CHAPITRE 5 4.2 Phénoménologie
- Page 210 and 211: 209 CONCLUSIONS limite sa représen
- Page 212 and 213: 211 CONCLUSIONS carbonates (anti co
- Page 214 and 215: 213 CONCLUSIONS à celui des mécan
- Page 216 and 217: 215 ANNEXES ANNEXE 1 ANNEXE 2 ANNEX
- Page 218 and 219: 217 ANNEXES A��EXE 2 CARTES DE
- Page 220 and 221: 219 ANNEXES 2cm (a) (b) Annexe 2.4
- Page 222 and 223: 221 ANNEXES 2cm (a) (b) profondeur
- Page 224 and 225: 223 ANNEXES Annexe 3.2 : Carte de m
- Page 226 and 227: 225 ANNEXES Annexe 3.6: Carte de la
- Page 228 and 229: 227 ANNEXES Annexe 4.2 : (a) Carte
- Page 230 and 231: 229 ANNEXES ∑ k i→ k / . Les ex
- Page 232 and 233: 231 ANNEXES (a) (b) Annexe 5.1 : (a
- Page 234 and 235: 233 ANNEXES La finalité de la mét
- Page 236 and 237: 235 BIBLIOGRAPHIE A BIBLIOGRAPHIE A
- Page 238 and 239: 237 BIBLIOGRAPHIE Castaing, R., 195
- Page 240 and 241: 239 BIBLIOGRAPHIE Gaucher, E., Robe
- Page 242 and 243: 241 BIBLIOGRAPHIE Lefranc, M., 2007
- Page 244 and 245: 243 BIBLIOGRAPHIE O Ohkubo, T., 200
196 CHAPITRE 5<br />
déviation de moins de 0,0001 %. Au-delà de 0,0001 %, il est uniquement possible d’obtenir<br />
une solution de type De/α = Da. En considérant logiquement une erreur expérimentale<br />
supérieure à 0,0001 % sur chaque point d’analyse, l’ajustement <strong>des</strong> profils microLIBS <strong>et</strong><br />
microsonde par l’équation 5.13 perm<strong>et</strong>tra uniquement d’obtenir un Da. A partir de<br />
l’ajustement analytique de la concentration <strong>dans</strong> le réservoir, les paramètres De <strong>et</strong> α pourront<br />
cependant être dissociés.<br />
α<br />
Fig. 5.30 : Evolution de la fonction FO (eq. 5.14) en fonction de De2 <strong>et</strong> α2 pour De1= 5×10 -<br />
10 <strong>et</strong> α1=100000). La ligne pointillée en blanc représente les coordonnées De1/α1=De2/α2.<br />
3.2.3 La méthode Time Domain Diffusion<br />
La méthode Time Domain Diffusion (cf. chapitre 3) a fait l’obj<strong>et</strong> d’étude afin de pouvoir<br />
interpréter les expériences d’in-<strong>diffusion</strong> analysées à partir de cartographies élémentaires par<br />
une approche de modélisation inverse. Néanmoins, à ce jour, celle-ci ne peut être utilisée <strong>dans</strong><br />
ce but.<br />
Interprétant une expérience d’out-<strong>diffusion</strong> avec une carte de <strong>porosité</strong> calculée par<br />
autoradiographie, Sardini <strong>et</strong> al. (2007) combinent la méthode TDD avec une procé<strong>du</strong>re<br />
d’inversion afin d’identifier les paramètres de <strong>diffusion</strong> de différents ensembles poreux. La<br />
courbe de restitution en traceur est alors utilisée comme donnée à ajuster. Pour pouvoir<br />
repro<strong>du</strong>ire c<strong>et</strong>te approche aux cartes de traceurs obtenues par in-<strong>diffusion</strong>, la méthode TDD<br />
doit considérer :<br />
D e (×1.10 -11 m²/s)<br />
Log (FO)<br />
3