Journal o f Hydrology, 3 7 ( 1 9 7 8 ) 1 8 5 - - 1 9 4 185

© Elsevier Scientific Publishing C o m p a n y , A m s t e r d a m -- Printed in The Netherlands

[4]

ANALYSE AUTOMATIQUE DES DONNI~ES DE POMPAGE D'ESSAI --

APPLICATION AUX HYPOTHI~SES DE BOULTON ET D'HANTUSH

J. h ; A N I A et M. S U C C H E
Universitd des Sciences et Techniques de Lille, 59650 Villeneuve d'Ascq (France)
(Re~u le 25 avril, 1977; accept~ p o u r publication le 17 aofit, 1977)

ABSTRACT

Mania, J. et Sucche, M., 1978. Analyse a u t o m a t i q u e des donn~es de pompage d'essai --

Application aux hypotheses de B o u l t o n et d'Hantush. ( A u t o m a t i c analysis of p u m p i n g
test data -- A p p l i c a t i o n to the B o u l t o n and Hantush hypotheses.) J. Hydrol., 37:
185--194.

Classical analysis of pumping test curves under unsteady-state flow conditions with
constant discharge rate in a fully penetrating well is based on a manual graphical compar-
ison of drawdowns with theoretical standard curves. This procedure leaves a relative sub-
jectivity to the operator. A m o r e rigorous m e t h o d is proposed: a u t o m a t i c analysis of
d r a w d o w n curves.
The principle of a u t o m a t i c analysis is based on the m i n i m i z a t i o n of differences b e t w e e n
observed d r a w d o w n s and those calculated theoretically according to a nonlinear law. We
use the iterative process in successive degrees of Gauss-Newton applied to the evaluation
of least-squares errors for a nonlinear regression.
The application of this m e t h o d to t w o hydrogeological mecanisms described by N.S.
Boulton and M.S. Hantush permits a quick and direct evaluation of the hydraulic param-
eters of the aquifer by means of partial derivatives with relation to those parameters of
the theoretical function. This latter is d e t e r m i n e d as a generalized integral of a f u n c t i o n
depending on four (first case) or five (second case) parameters. We justify here the
permutations of the derivations and integration.
At the level o f numerical calculation, the difficulty is to evaluate with sufficient
precision the generalized integrals containing Bessel function.
Applications of a u t o m a t i c analysis have been carried out and are discussed.

RESUME
L'analyse classique en r6gime transitoire des eourbes de descente du niveau d ' u n e nappe
d'eau souterraine lors d ' u n p o m p a g e ~ d~bit constant dans un puits parfait repose sur une
interpolation manuelle du trac6 de la courbe des r a b a t t e m e n t s qui laisse une part de sub-
jectivit6 relative ~ l'op~rateur. Une m 6 t h o d e plus rigoureuse est ici propos~e: l'analyse
a u t o m a t i q u e des courbes de descente.
Le principe de l'analyse a u t o m a t i q u e est bas~ sur la minimisation des ~carts entre les
r a b a t t e m e n t s e x p ~ r i m e n t a u x et les r a b a t t e m e n t s calculus t h ~ o r i q u e m e n t selon une loi non
lin~aire. Nous utilisons le processus it~ratif par degr6s successifs de Gauss-Newton appliqu~
l'~valuation des erreurs sur les m o i n d r e s carr~s d ' u n e r6gression non lin~aire.
L'application de cette m ~ t h o d e ~ d e u x schemas hydrog~ologiques d~crits par N.S.
186

Boulton et M.S. Hantush permet d'obtenir rapidement et directement les param~tres

hydrauliques de l'aquif~re en utilisant les d~riv~es partielles par rapport ~ ces param~tres
d'une fonction th4orique. Cette derni~re est d~termin~e comme int4grale g~n~ralis4e d'une
fonction d4pendant de quatre (premier cas) ou de cinq param~tres (second cas). Nous
justifions ici la permutation des d~rivations et de l'int4gration.
Au niveau du calcul num~rique la difficult~ est d'~valuer avec suffisamment de precision
des int~grales g~n~ralis~es contenant une fonction de Bessel.
Des applications de l'analyse automatique ont ~t~ r~alis4es et sont discut~es.

INTRODUCTION

L ' i n t e r p r ~ t a t i o n des p o m p a g e s d'essai reste u n des meilleurs m o y e n s d'in-

vestigation en h y d r o g ~ o l o g i e p o u r la r e c h e r c h e des param~tres h y d r a u l i q u e s
d ' u n aquif~re. L ' a n a l y s e m a n u e l l e des valeurs de r a b a t t e m e n t (Forkasiewicz,
1 9 7 0 ; K r u s e m a n et al., 1 9 7 4 ) est sujette ~ la subjectivit~ de l'op~rateur. N o u s
p r ~ s e n t o n s u n e analyse a u t o m a t i q u e bas~e sur la m i n i m i s a t i o n des ~carts
e n t r e les r a b a t t e m e n t s th~oriques et les r a b a t t e m e n t s e x p ~ r i m e n t a u x . Cette
analyse qui a d~j~ ~t~ tent~e avec succ~s dans le cas de l ' h y p o t h ~ s e de Theis
{Mania et S u c c h e , 1 9 7 6 ) est test~e p o u r les h y p o t h e s e s b e a u c o u p plus com-
plexes de B o u l t o n et d ' H a n t u s h .

RAPPEL DES SCHEMAS HYDRODYNAMIQUES DE BOULTON

D e u x t y p e s de schemas h y d r o d y n a m i q u e s r ~ p o n d a n t ~ des crit~res d'hori-

zontalit~, de continuitY, d ' i s o t r o p i e et d ' h o m o g S n ~ i t 8 s o n t envisageables.

(1) L e p r e m i e r s c h d m a c o n c e r n e une n a p p e libre r e p o s a n t sur u n s u b s t r a t u m

h o r i z o n t a l i m p e r m e a b l e ( F i g . l ) . A u c o u r s du p o m p a g e , o n assiste ~: l'appari-

Puit's de production Piezom:ei'r'e log s

- 7 sol

= o

E ,~?

• .
. . . .
i i
I',
,
o
AqulFere , , o o o ~ o o

i o o 11 o 1 ~ 1 -

I " °"°
',l-i ' I °
', o I,
L ' II

- ~- Couche irnper'r'n;able-- / log t

Fig.1. Cas 1 : ~gouttement d'une nappe libre.

Fig.1. Case 1 : delayed yield of unconfined water.

Q = discharge of the well at constant rate; s = drawdown in well; nos - original piezometric
surface (static); nod '= pumping piezometric surface (dynamic); r = distance to observation
well.
 187

tion d'un drainage de la partie sup6rieure de l'aquif~re encore appel~ "l'6gout-

t e m e n t " qui est ~ l'origine d ' u n palier de stabilisation des rabattements. Ce
palier dispara~t d~s que les variations de pression atteignent l'ordre de gran-
deur de la vitesse d'abaissement du niveau pi~zom~trique.
Boulton (1963) a interpr~t~ un tel sch6ma en proposant une relation
analytique entre les rabattements mesur~s en fonction du temps et de quatre
autres param~tres:
s = (Q/4nT) f(S,S',c~,r/B,t)
o/1
t = temps de pompage
S = coefficient d'emmagasinement instantan~ de la couche aquif~re
S'-- coefficient d'emmagasinement retard~ de la couche aquif~re
T = transmissivit6 de l'aquif~re
Ot -- indice de retard
r -- la distance entre le pi~zom~tre d'observation et l'ouvrage
B = (T/aS')1"2 = le facteur d'~gouttement
Q = le d~bit de pompage {constant)

(2) L e d e u x i ~ m e s c h d m a d'application est relatif ~ une nappe semi-captive

limit6e ~ son toit et ~ son mur par des 6pontes imperm~ables (Fig. 2). Un
aquif~re annexe (sup~rieur ou inf~rieur) se draine vers l'aquif~re principal en
cours de pompage et provoque c o m m e dans le premier cas l'apparition d'un
palier de stabilisation. L'6ponte semi-impermeable permet le transfert et
poss~de un coefficient d'emmagasinement non n~gligeable. I1 est possible
d'utiliser la fonction pr~sent~e auparavant o ~ le coefficient B 6gal/t ( T b ' / K ' ) 1/2
s'appelle facteur de drainance avec: b' = ~paisseur de l'~ponte semi-permeable;
et K' = perm~abilit6 verticale de l'6ponte semi-permeable.
P u i r s de p_r'oduci'ion
P i ~ z orn~i- r, e - - ~Q--Q Pi~zombi'r" e log s
'--'-~ ..... /f e / __ sol ~a

0
It' (] n rl e X e "'" o °o
o
-
C
-o L
A~uiF~re ii I I I
I I I ,
principal ',Iiiii I ',i I Iii /"
.2: . . . . . . .
il i I li OU

"- . . . . . . .
srabilisation

--_-_-_-_-Lc_o_c_h p e5 _%0a2 ,_%-:_---- log r

Fig. 2. C a s 2 : d r a i n a n c e e n m i l i e u s e m i - c a p t i f .

Fig. 2. C a s e 2 : l e a k a g e in a l e a k y a r t e s i a n a q u i f e r .
For symbols used refer to caption of Fig.1.
188

La condition suivante est requise t <~ b ' S ' / l O K' pour appliquer la m~thode.
I1 faut signaler dans ce cas l'existence d'une autre m~thode mise au point
par Hantush (1956, 1959) pour t > 2 b ' ( S ' / K ' ) .

EXPOSI~ DE LA Mt~THODE

La fonction s(t) d~pendant du temps et d'un ensemble de param~tres

indique de fa~on th~orique l'~volution du niveau pi~zom&rique dans un puits
implant~ dans une couche aquif~re r~pondant ~ l'un des schemas d~crits pr&
c~demment. Des observations de ce niveau, de valeurs o(ti) , sont faites pour
des temps ti diff~rents et permettent de d~terminer les param~tres de telle
sorte que s(t) rende compte, au mieux, du ph~nom~ne observe.
Ce probl~me a ~t~ abord~ par Boulton (1963) en comparant directement
un ensemble de courbes th~oriques A celles obtenues exp&imentalement. Pour
am~liorer cette premiere estimation des param~tres, on utilisera la m & h o d e
des moindres carr~s: les meilleurs param~tres seront ceux qui rendent la
somme zn=l [s(ti) - o(ti)] 2 minimum.
Programm~ par Sampson (1971), cette m~thode fait intervenir les d&iv~es
de s par rapport aux diff&ents parambtres: partant d'un ensemble de para-
mbtres obtenu empiriquement ou par la m~thode de Boulton, l'ordinateur
d&ermine un nouvel ensemble de parambtres rendant la somme des carr~s
aussi petite que le d~sire l'utilisateur.

NOTATIONS -- EXPRESSIONS DE s(t) ET DE SES DI~,RIVt~ESPARTIELLES

P a r a m h t r e s ~ d k t e r m i n e r : S, S ' , ~, T dans le premier schema; et S, S', a, K',

T dans le second.

N o t a t i o n s interm~diaires:
1 + S'/S; k = at/2; y = [72(1+x2) 2 - 4~7x2]V~

f = [ch(ky) + 77(1 - x 2 ) / y sh(ky)] e x p [ - k v ( 1 + x2)]; v = r/B

et
p = m71/2(~-1) -'/~

To est la fonction de Bessel d'ordre 0.

E x p r e s s i o n de s(t) e t d e ses d~riv~es partielles:
7o
W =J (2/x)(1-f)To(pX) dx; s(t) = (Q/4~rT) W
0

Ainsi par l'interm~diaire de W, s(t) apparMt comme int~grale g~n&alis~e

d'une fonction d~pendant de param~tres l'existence des d&iv~es partielles
n~cessite une justification th~orique. Leurs expressions ~tant diff~rentes dans
les deux schemas envisages.
 189

Dans le premier schema:

[ r )a___W] as Q s' aw
aT 4~T ~ 2 \ T av ' aS 4eT S 2 aV

aS' 4;r s +-2r 1 et

a~ Q t--+ rt-~] -~v

8~T ak

Dans le deuxi~me schema:

+r K' '~ aW as Q s' aw
aT 4~T 2 2 ' aS 4uT S 2 a~

as Q 1 aw as Q aw 3s Q r aW
- - t - - et - -
aS' 4~T S aT? ' aa 8~rT ak aK' 81rT ( K ' T b ' ) v2 av

JU S TI F IC ATI ON DE L'EXISTENCE DES DI~RIV]~ES PARTIELLES

Une meilleure connaissance des fonctions ~ int6grer est apport6e par les
majorations suivantes:
Pour x > 0; ~ > 1; k e [a,b]c]0,¢¢[:

0 < 1/y < 1~Ix 2 - 11, 1/(1 + x 2) < ~/Y < l l l x 2 - II

0 < e x p [ - k ~ ( l + x 2) + ky] < exp[-2ax2/(1 +x2)], et

0 < exp[-k~(1 +x 2) - ky] <. e x p [ - 2 a x 2]

Pour x ~> 1:

0<l+~(1-x2)/yg2/(l+x2), 2x21(l+x2)<.l-~(1-x2)/y<~ 2
O<.l+2/y-7?(l+x2)/y<. 2/(x2 - 1), O<1-21y+~?(l+x2)/y<2x2/(x2-1)
0 <. ~(x 2 - 1)/y 3 <- l l ( x 2 - 1) 2
2x 2 ( l + x 2)
0< (~ly 2) [~(x2 +1) 2 - (3x~ +1)] + (11y)[~(x~+1)-1] < (x 2 - 1 ) 2

et
_ 2(x 2 + l ) / ( x 2 _1)2 < (~/y2) [~?(x2 +1)5 _ (3x 2 +1)] - ( l / y ) [~?(x2 +1) - 1 ] < 0

Compte tenu des in6galit6s pr~c6dentes et de l'expression asymptotique

des fonctions de Bessel au voisinage de l'infini (Watson, 1958):
T p ( x ) = (2/nx) 1/2 c o s [ x - (2p+1)(~/4)] + O ( x - m )
190

on m o n t r e la convergence u n i f o r m e des int6grales g6n6ralis6es qui intervien-

n e n t dans OW/Ok, OW/O v e t OW/07? par r a p p o r t aux diff6rents param6tres
dans t o u t segment de R ÷ p o u r k et v, de ] 1,~o[ p o u r 77.
Convergence u n i f o r m e acquise par calcul direct p o u r J Tl(pX)dx sachant
1

que TI(X) = -To(X) et en utilisant la t e c h n i q u e suivante darts les autres cas:

on majore la valeur absolue des fonctions ~ int6grer entre 1 et oo par une
f o n c t i o n i n d 6 p e n d a n t e des param6tres, 6quivalente au voisinage de l'infini a
une expression de la f o r m e (K/xC~)lTp(pX)[ avec K r6el positif, a > 1/2 et p = 0
ou 1.
L'expression a s y m p t o t i q u e de Tp(pX) p e r m e t de conclure:
En r e m a r q u a n t enfin que seule l'int6grale de W e s t g6n6ralis6e en 0 et con-
verge, les th6or6mes classiques de p e r m u t a t i o n des d6rivations et des int6grales
g6n6ralis6es p e r m e t t e n t d'6crire les expressions des d6riv6es partielles:

OW/Ok = 277 f x ( ( e x p [ - k ~ ( 2 + x 2) + ky] )[1 + 2/y - 77(1 +x2)/y] +

0

+ ( e x p [ - k r l ( l + x 2 ) - k y ] ) [1 - 2/y + rl(l+x2)/y] }To(px)dx

OW/Ov = 2~1/2(r/-1)-~/2f0 ( f - 1 ) T , ( p x ) d x et 3 WIOrl = 11 +12 + 13 + 14 + Is

avec

11 = v(7/-1)-1[r/07-1)] -lh f (1-f)Tl(px)dx

0

I2 = 2k S x(exp[-k~7(l+x2) +ky] } {(~?/y2) [r/(x 2 +1)2 _ (3x 2 +1)] -

0

- ( l / y ) [7/(x 2 +1) - 1] } To(px)dx

oo

I3 = - 2 7 f0 [x(x2 -1)/YO] ( e x p [ - k r / ( 1 +x 2) + Icy] } To(px)dx

14 = 2k f x { e x p [ - k r l ( l +x 2) - ky] } ((r~/y2)[rl(x 2 +1) 2 - (3x 2 +1)] +

0

+ (1/y)[rl(x 2 +1) - 1] } To(px)dx

Is = 217 f [x(x 2 - 1 ) / y 3] (exp[-kr/(1 +x 2) - k y ] ) To(px)dx

0
 191

Pour l'6valuation num6rique des d6riv6es partielles, le changement de

variable X = px ram6ne chaque expression ~ la forme:

.f g(X)To(X)dX ou ; g(X)T,(X)dX
0 0

avec g fonction de signe constant. On introduit les valeurs approch6es

s, . . . . . Sn des z6ros de To (ou de T1) et on prend comme valeur approch6e de
l'int6grale une somme partielle du rang n de la s6rie altern6e dont le terme
g6n6ral est:
sn fo
fs g(X)To(X)dX ou g(X)T, (X)dX
n --I 8 n --I

APPLICATION ET CONCLUSION

Les expressions successives des d6riv6es partielles de la fonction de Boulton

(1963) par rapport ~ ses param6tres ont 6t6 programm6es et adapt6es au
programme g6n6ral de Sampson (1971). Les calculs ont 6t6 r6alis6s sur
l'ordinateur CII 10070 de l'Universit6 de Lille. Pour 6valuer les param~tres
hydrauliques d'un aquif6re le programme ainsi r6alis6 poss6de cependant ses
propres limites. L'existence 6ventuelle de minimums locaux pour la somme
minimiser et la complexit6 des fonctions augmentent consid6rablement le
temps de calcul et obligent l'utilisateur ~ introduire dans ce programme un
ensemble de param6tres assez proches de la r6alit6. La m6thode statistique
serf alors ~ affiner les valeurs de ces param6tres connues jusqu'~ pr6sent de
faqon fort subjective.
La m6thode a 6t6 test6e pour les deux sch6mas hydrog6ologiques pr6sent6s
pr6c6demment et sur des pompages d'essai effectu6s dans la nappe des alluvions
de la Garonne ~ St. Pardon de Conques (Forkasiewicz, 1970) et dans la nappe
du S6quanien de la r6gion boulonnaise (Mania et Ricour, 1974).
Les essais ont 6t~ effectu6s ~ un d~bit constant et les puits ont 6t6 r~alis6s
de mani6re ~ capter t o u t l'aquif6re test6.
Nous avons report6 sur le Tableau I e t traduit sous forme de graphiques
sur la Fig. 3 les valeurs des rabattements mesur6s sur le terrain (courbe 1).
Nous les comparons successivement ~ ceux issus du calcul automatique
(courbe 2) et ~ ceux d6duits par calcul en reprenant les param6tres obtenus
par la m6thode manuelle (courbe 3).
On remarquera l'am61ioration de l'ajustement. Les valeurs des param6tres
calcul6es par l'ordinateur (Tableau II) diff6rent de celles obtenues par la
m6thode manuelle. Dans le cas pr6sent6 nous sommes partis de valeurs de
param6tre issues de la m6thode manuelle afin d'6courter le temps de calcul.
L'erreur sur les moindres carr6s, apr6s sept it6rations, atteint le chiffre de
9 , 6 . 1 0 .4 , ce qui est acceptable.
192

TABLEAU I

Valeur des rabattements obtenus par calcul ~ partir de param~tres calculus par ordinateur
et ~ partir de param~tres issus de l'analyse manuelle. Cas de l'~gouttement d'une nappe.
Exemple: St. Pardon de Conques (in Forkasiewicz, 1970)

Rabattements Rabattements calculus (m) Temps

mesur~s sur de pompage
le terrain ~ partir des fi partir des (min)
(m) param~tres param~tres
issus de l'ana- obtenus auto-
lyse manuelle matiquement

0,206 0,722 0,218 10

0,211 0,736 0,220 12
0,216 0,756 0,225 16
0,225 0,787 0,240 40
0,250 0,815 0,260 83,33
0,270 0,862 0,285 166,66
0,316 0,950 0.35 466,66

REMERCIEMENTS

Les frais d'informatique ont ~t~ support~s par le Service G~ologique

R~gional du B.R.G.M. de Lezennes; aussi nous remercions Monsieur G. Souliez,
Directeur de ce service, pour son aide pr~cieuse.
Ces travaux ont ~t~ men~s ~ bien grfice ~ l'aide de ~ m e . R e m y et de M.
Hayez du Centre Interuniversitaire de Traitement Informatique que nous
remercions vivement.

log s

3 ~ ~
manuelle
/-

0,5

experiment al

//~6:6=6_~ 2 ~ autornattque

log t (nn n)
0~I
' ' ' s ' ~ . . . . . . . '~'o ~ ' ...... 1'o'

Fig. 3. Cas 3: comparaison des diverses m~thodes et restitution des rabattements.

Fig. 3. Case 3 : comparison of miscellaneous m e t h o d s and restoration o f d r a w d o w n s .

Curve 1 : experimental or observed data curve; curve 2: automatic computed data; c u r v e 3:
manual curve (field data analysed, w i t h t y p e curves on bi-logarithmic paper).
 193

TABLEAU II

Comparaison des param~tres obtenus par la m~thode manuelle et par la m6thode automa-
tique

(A) Cas de l'~gouttement d'une nappe (correspond aux rabattements du Tableau I)

Exemple: St. Pardon de Conques (in Forkasiewicz, 1970):

Param~tres Transmis- Coefficient d'emmagasinement Indice de retard,

cherch6s sivit6, T (s -1)
(m:/s) S S'
apparent retard~

M~thode automatique 3,6.10 -2 9,9o10 -4 4,9.10 -2 4,9-10 s

septi~me iteration
M~thode manuelle 1,9.10 -2 1,5-10 -3 4,5o10 -2 5,5o10 -s

(B) Cas d'une drainance ~ travers u n semi-permeable d'~paisseur b' (ici 3,50 m). Exemple:
St. L~onard (in Mania et Ricour, 1974):

Param~tres Transmis- Coefficient Indice de Perm~abilit~

cherch~s sivit~, T d'emmagasinement retard, a du semi-
(m2/s) (s -1 ) permeable, k'
S S' (m/s)
apparent retard6

M~thode automatique 8,7.10 -3 1-10 -s 1,8.10 -3 6,8.10 -6 2,9.10 -8

septi~me iteration
M~thode manuelle 7,5-10 -3 1,10 -4 1,9.10 -3 6,8o10 -6 4,6o10 -~

BIBLIOGRAPHIE

Boulton, N.S., 1963. Analysis of data from non-equilibrium pumping tests allowing for
delayed yield from storage. Proc. Inst. Civ. Eng. (London), 26(6): 693.
Forkasiewicz, J., 1970. Guide bibliographique m~thodique sur l'interpr~tation des donn~es
de pompages d'essai (~valuation des param~tres aquif~res). Bull. Bur. Rech. G~ol.
Mini~res (Fr.), Sect. 2, 3'(1): 67--72.
Hantush, M.S., 1956. Analysis of data from pumping tests in leaky aquifers. Trans. Am.
Geophys. Union, 37(6): 702--714.
Hantush, M.S., 1959. Non-steady flow to flowing wells in leaky aquifers. J. Geophys. Res.,
64(8).
Kruseman, G.P., De Ridder, N.A. et Meilhac, A., 1974. Interpretation et discussion des
pompages d'essai. ILRI, Wageningen, 213 pp. (2~me ~d.).
Mania, J., 1976. Un exemple d'analyse automatique par ordinateur des rabattements des
nappes d'eau souterraine -- Application ~ l'hypoth~se de Theis. Ann. Soc. G~ol. Nord,
96: 113--119.
Mania, J. et Ricour, J., 1974. Hydrog~ologie et alimentation artificielle de la nappe du
S~quanien dans la r~gion de St. L~onard (P.-de-C.). Ann. Soc. G ~ l . Nord, 94: 133--147.
Mania, J. et Sucche, M., 1976. Mise au point d ' u n programme d'ordinateur pour l'analyse
automatique des courbes de descente -- Application aux hypotheses de Theis et de
Boulton. 4~me R~un. Annu. des Sciences de la Terre, Paris, 13--16 avril 1976, p. 280.
194

Saleem, Z.A., 1970. A computer method for pumping-test analysis. Ground Water, 8(5):
21 --24.
Sampson, P.F., 1971. Program BMD 07R: Nonlinear least squares. Health Sciences
Computing Facility, University of California, Los Angeles, Calif.
Watson, G.N., 1958. Treatise on the Theory of Bessel's Functions. Cambridge University
Press, London, 900 pp.