REPUBLIQUE TUNISIENNE

Ministère de l’Enseignement Ministère de l’Agriculture

Supérieur de la recherche scientifique et et des Ressources hydrauliques
de la technologie

Institution de la Recherche
et Jendouba
Université de de l’Enseignement
Supérieur Agricoles

Ecole Supérieure des Ingénieurs

de l’Equipement Rural de Medjez El Bab

Elaborée par:

Kheder Sofien
Abdellaoui Aymen

1
EXERCICE I :

I. Présentation du problème  :
Un puits pénètre totalement dans une nappe captive est pompé à raison d’un débit Q.
A l’instant t = 0 le niveau piézométrique est h0, après l’instant t un rabattement est crée à
une distance r du puits. Ce problème est représenté par l’équation de diffusivité suivante :

∂2 h ∂ 2 h S ∂ h
+ =
∂ x2 ∂ y2 T ∂ t

II. Calcul du rabattement par utilisation du programme

II-1 -Données du problème 
Le rayon r = 100 m
La transmissivité T = 300 m2 /j
Le coefficient d’emmagasinement S = 0.002
t0 = 0.01 jours :

II-2-Résultats des calculs

Le calcul du rabattement à 100m du puits de pompage est calculé par l’utilisation du
programme  pour :

 Alpha = 0.5 (Crank-Ncolson ; Méthode implicite)

 Alpha = 1 (Méthode totalement implicite)
 Alpha = 0 (Méthode explicite)

2
 h0-h (m) h0-h (m) h0-h (m) Nbre
T (jours) t (jours) Nbre Itér Nbre Itér
(α =0) (α =0.5) (α =1) Itérations
0.01 0.01 0 1 0.05 3 0.06 3
0.05 0.02 0.46 1 0.41 4 0.38 4
0.13 0.38 -- 1 0.88 5 0.81 6
0.49 0.17 -- 1 1.52 7 1.44 9
1.13 0.38 -- 1 1.92 10 1.79 10
5.82 1.95 -- 1 2.48 9 2.18 5
13.12 4.39 -- 1 2.57 3 2.22 1

Tableau 1 : Valeurs de rabattements calculés par le programme

II-3-Courbes de rabattement en fonction de log (t) pour les différents α

 Courbe de rabattement en fonction de log (t) pour α = 0

rabattement pour alpha=0

0.5
0.45
0.4
0.35
rabattement (m)

0.3
0.25
0.2
0.15
0.1
0.05
0
0.01 0.1
logt

 Courbe de rabattement en fonction de log (t) pour α = 0,5

3
 rabattement pour alpha=0,5
3

2.5

2
Rabattement (m)
1.5

0.5

0
0.01 0.1 1 10 100
logt

 Courbe de rabattement en fonction de log (t) pour α = 1

rabattement pour alpha=1

2.5

2
rabattement (m)

1.5

0.5

0
1 2 4 8
logt

Interprétation des résultats :

4
  La méthode explicite (α =0) donne des résultats instables pour différentes

TΔt
valeurs de t car : SΔx2 <0.25.

TΔt
T (jours) t (jours) Remarque
SΔx2
0.01 0.01 0 .15 Stable
0.05 0.02 0.3 Instable
0.13 0.38 5.7 Instable
0.49 0.17 2.58 Instable
1.13 0.38 5.7 Instable
5.82 1.95 29.25 Instable
13.12 4.39 65.85 Instable

Tableau 2 : Tableau d’interprétation

 Pour réduire l’instabilité, il faut augmenter les valeurs de t .

 Les valeurs déterminées par la méthode implicite et totalement implicite sont trop
proches et donnent presque les mêmes résultats.

III .Méthode de Theis  :

5
Pour un pompage à débit constant Q dans une nappe initialement au repos, dans un
aquifère de caractéristiques hydrauliques, T, S, THEIS a établi la formule suivante donnant
le rabattement à la distance r de forage ou est exercé le pompage, au bout du temps t après
le début de celui-ci :
∞
dx
t 2 ∫ exp(−x )
Q ∫ exp ( −r4 T s ) dθθ Q r2
x
S(r, t) = 4 πT 0 = 4 πT 4 Tt

t=0 ⃗ S(r, 0) = 0
t>0 ⃗ S( r, 0 ) > 0
∞
r2S Q du Q
∫ exp(−u) u
En posant: u = 4 Tt ⇒ S(r, t) = 4 πT U = 4 πT W (u)
W (u) est une fonction exponentielle intégrale
Pour résoudre cet exercice avec la méthode de THEIS, on détermine W (u) à partir de
l’abaque pour différentes valeurs de (1/u) calculés.
Q
Par suite on détermine : s = h0 – h = 4 πT W (u)

r2S
avec : u = 4 Tt
T (jours) u 1/u W (u) s (m)
0.01 1.666 0.6 0.08 0.042
0.05 0.333 3 0.8 0.424
0.13 0.128 7.8 1.70 0.902
0.49 0.034 29.4 2.9 1.540
1.13 0.014 67.8 3.9 2.070
5.82 0.0028 349.2 5.5 2.917
13.12 0.0012 787.2 6.5 3.448

Tableau 3 : Calcul de rabattements par méthode de Theis

6
 Méthode de Theis
4
3.5
3
2.5
2
Rabattement (m)
1.5
1
0.5
0
0.01 0.1 1 10 100
log t

Courbe de rabattement en fonction de log (t)

IV. Méthode de Jacob  :

Cette méthode n’est applicable que pour 1/u > 100 ; pour notre cas cette condition n’est
vraie que pour deux valeurs de t, pour les autres valeurs de t la condition n’est satisfaite.
Donc l’erreur de calcul de rabattement augmente pour des 1/u <100.

T (jours) u 1/u Remarque s (m)

0.01 1.666 0.6 1/u <100 --
0.05 0.333 3 1/u < 100 --
0.13 0.128 7.8 1/u < 100 --
0.49 0.034 29.4 1/u < 100 --
1.13 0.014 67.8 1/u < 100 --
5.82 0.0028 349.2 1/u > 100 2.801
13.12 0.0012 787.2 1/u > 100 3.230

Tableau 4 : Calcul de rabattements par méthode de Jacob

7
 Méthode de Jacob
3.3
3.2
3.1
3
2.9
Rabattement (m)
2.8
2.7
2.6
2.5
0.01 0.1 1 10 100
log t

V. Comparaison des résultats  :

h0-h (m) h0-h (m) h0-h (m)
T (jours) t (jours) Theis Jacob
(α =0) (α =0.5) (α =1)
0.01 0.01 0 0.05 0.06 0.04 --
0.05 0.02 0.46 0.41 0.38 0.42 --
0.13 0.38 -- 0.88 0.81 0.902 --
0.49 0.17 -- 1.52 1.44 1.540 --
1.13 0.38 -- 1.92 1.79 2.070 --
5.82 1.95 -- 2.48 2.18 2.917 2.801
13.12 4.39 -- 2.57 2.22 3.448 3.23

Tableau 5 : Valeurs du rabattement calculé par les différentes méthodes

Les valeurs trouvé par:

 méthode totalement implicite (α = 1)
 méthode implicite (Crank-Nicolson) (α = 0.5)
 Sont trop proche, et indique presque les mêmes valeurs de rabattements

8
Les valeurs trouvé par:
 méthode de Theis
 méthode de Jacob
⇒ Donnent des valeurs trop différente de celles donner par la méthode implicite et
totalement implicite ; plus le temps (T) augmente plus la différence entre les deux
méthodes augmente.

 La méthode explicite est pratiquement instable

9
EXERCICE II :
1°/ L’équation différentielle qui gouverne l’écoulement dans cette nappe est la
suivante :

2
∂ h S∂h
2
=
∂x T∂t

2°/Développement de la méthode explicite :

t
∂2 h S h t+1 −ht
( )
∂ x2
=
T ∆t
hi+1 +h i−1 −2h ti S hti +1−hti
t t
=
∆ x2 T ∆t
2 h́ ti−2 hti S hti +1−h ti
hi +1−hi−1=2 h́i ↔
t t
=
∆ x2 T ∆t
t t
∆ t ∙T 2 h́i −2 hi
hti +1−hti =
S
∙( ∆ x2 )
t t
T ∆ t 2 h́ i−2 hi
t +1
h =
i
S (∆x
t
2 )
+hti

T ∆ t 2 h́i T ∆t
hti +1=
S ∆x 2 (
+ hti 1−2
S ∆ x2 )
3°/

x° x1 x2 x3 x4 x5 x6 x7 x8 x9 x10
T ∆t
4°/ Pour avoir une solution stable il faut avoir : <0,5
S ∆ x2

5°/h=f(x) :

10
 h=f(x)
18

16

14

12

t=0
10
t=10
h t=100
8
t=400

0
0 20 40 60 80 100 120

6°/

11
H(T400) numérique H(T400) analytique

12
 16.00
16
68.78
15.5
26.44
15
10.55
14.5
5.19
14
5.04
13.5
9.92
13

24.73
12 .5
64.27
12
168.09
11.5
11.00
11

7°/

13
 h=f(t) pour dt=5
18
16
14
12
10
8
6
4
2
0
0 50 100 150 200 250 300 350 400 450

h=f(t) pour dt=8

18
16
14
12
10
h(m) 8
6
4
2
0
0 50 100 150 200 250 300 350 400 450
t(min)

8°/

14
 h=f(x)
18
16
14
12 t=0
10 t=10
h(m) 8 t=100
6 t=400
4
2
0
0 20 40 60 80 100 120
x(m)

15