Hydroécol. Appl. (1992) Tome 4 Vol 2, pp.

23-31

Sim ulation m athém atique des nappes

d’hydrocarbures et com paraison
avec les observations par télédétection

M athem atical sim ulation of oil spill at sea

and com parison with observations
by remote sensing

K. M a ro ih i
G eohydrodynam ics a n d E nvironm ental Research, (GHER), U niversity o f Liège. B5. S a rt Tilman,
B -4000 Liège, Belgium

E. D e le e rs n ijd e r
Institute o f A stronom y a n d Geophysics, Catholic U niversity o f Louvain, 2, chem in du Cyclotron,
B - 1348 Louvain-la-Neuve, Belgium

A . L o ffe t
W alphot Detection S.A., rue Van O pré 97, 5 100 Namur. Belgium

Résumé. - On présente un modèle mathématique du transport et de l'étalement des

nappes d’hydrocarbures à la surface de la mer.
Le modèle se distingue des modèles antérieurs par sa capacité à prendre en compte
simultanément les effets de gravité, friction, et de tension superficielle, - permettant
ainsi la modélisation des nappes expérimentales de petites dimensions comme des
grandes nappes accidentelles - , et par une nouvelle paramétrisation de la tension su­
perficielle et du frottement interfacial, mieux adaptée aux conditions réelles en mer.
La performance du modèle est évaluée par la modélisation d'expériences in situ où la
position et l’épaisseur de la nappe sont mesurées par télédétection. Un excellent accord
est trouvé entre les prévisions du modèle et les observations.

Abstract. - A mathematical model of oil spill’s transport and spreading is presented.

The model differs from previous ones by its capability of taking simultaneously gravity,
friction and surface tension into account - allowing the simulation of small experimental
spills as well as large accidental spills - and by a new parameterization of surface
24 K. Maroihi et E. Deleersnijder

te n s io n a n d in te rfa c e s tre s s e s , b e tte r a d a p te d to re a l fie ld c o n d itio n s .

T h e m o d e l is te s te d b y th e s im u la tio n o f in s itu e x p e rim e n ts w h e re th e p o s itio n a n d th e
th ic k n e s s o f th e s p ill a re m e a s u re d b y re m o te s e n s in g . A n e x c e lle n t a g re e m e n t is fo u n d
b e tw e e n th e m o d e l’s p re d ic tio n s a n d th e o b s e rv a tio n s .

INTRODUCTION La m odélisation du transport et de

l’étalem ent des nappes d ’hydrocar­
bures est très difficile. Elle requiert
Le développem ent intensif des fo ­ une bonne connaissance des co nd i­
rages « o ff-sho re » et le grand trafic tions m étéorologiques et o céanogra­
des tankers ont augm enté les risques phiques et des caractéristiques
de pollution par hydrocarbures. physico-chim iques des hydrocar­
bures.
L’im portance des dom m ages éco­
logiques et économ iques causés par Le présent m odèle est capable d ’é­
les nappes d ’hydrocarbures, déver­ va lu er le transport général sous
sées a ccidentellem ent ou volo ntaire ­ l’action com binée du vent et des cou­
ment, sont tels q u ’ils ont suscité un rants, et son étalem ent progressif
vaste effort international dans l’am é ­ sous l’action com binée de la gravité,
lioration et le p erfectionnem ent des de la tension superficielle et de la fric­
techniques de protection et de net­ tion. L’évaporation, l’entraînem ent des
toyage et dans l’im plantation des sys­ hydrocarbures dans la colonne d ’eau
tèm es de détection et de pistage. et de l’eau dans la nappe peuvent être
tenus en com pte.
La télédétection peut jo u e r un rôle
m ajeur dans les problèm es de pollu­ On peut syn thé tiser le m odèle en
tion par les hydrocarbures. Ce rôle une équation d’évolution de l’é pais­
consiste en deux p o in ts: seur de la nappe d ’hydrocarbures. Ce
param ètre contient les cara ctéristi­
1) La recherche et la localisation ques essentielles du phénom ène et
des régions polluées. C ette tâ che est présente l’avantage d ’utiliser les don­
principalem ent une tâ che de su rve il­ nées de la télédétection.
lance incluant la détection des o pé ra ­
tions de décharges illégales par les
navires passants, leur identification et
la fourniture de preuves aux trib u ­ LE TRANSPORT
naux.

2) La cartographie et le pistage de À la surface de l’eau, une nappe d ’hy­

la nappe et la classification par é pa is­ d rocarbures a tendance à dériver
seur des hydrocarbures en vue d ’as­ sous l’action com binée du vent et des
sister les opérations de nettoyage. courants m arins (m arées, tem pêtes,
 Simulation mathématique des nappes d’hydrocarbures 25

circulation résiduelle...). A insi, le vec­ te rve n ir la viscosité m oléculaire. Pour

te u r vite sse de d érive du ce ntre de cette raison, ce m odèle n’est apte à
m asse de la nappe e st don né p ar: d écrire que des expériences en labo­
Vd = V v + Vc (1) ratoire ou en m er trè s calm e.
La param étrisation utilisée dans le
où Vv est le v e cte u r vite sse de dérive
m odèle de Nihoul est identique à celle
dû au vent et V c celui dû aux courants
proposée par W arluzel et Benque
m arins. L’équation du m ouvem ent du (1981). Elle e st la m ieux adaptée aux
centre de m asse s ’é crit a lo rs: co nd itio ns réelles en m er (Nihoul,
1983, 1984).
<2 >

D ans le m odèle de Nihoul (ce m o­

Vv est dévié à d roite de la direction d èle sera d écrit briève m en t par la
du vent d’un angle a (dérive d ’Ekm an) suite), la surface (S) de la nappe est
et sa g ra nd eu r peut être reliée au vent reliée au te m p s (t) par les relations
à 10 m de la su rfa ce par d e s relations suivantes (N ihoul, 1983):
se m i-em piriq ue s (N ihoul, 1983). L’an­ R égim e inertie-gravité
g le de dérive p e u t être relié à l’inten­
S - g '1/ 2 Q 1 /2 t (3)
sité du vent et à la latitude.
Les co uran ts m arins so n t d éte rm i­ R égim e g ravité-friction
nés par les obse rva tion s ou par des S - g' 1/3 0 2/3 k~1/3 p 1/ 3 11/3 (4)
m odèles h ydrodynam iques régionaux.
R ég im e friction-tension superficielle
S ~ y k~1 1 (5)
L’ÉTALEMENT Q e st le vo lum e déversé, k le co effi­
cien t de friction à l’interface entre les
hydrocarbures et l’eau, y la tension
Fay et Nihoul o nt déve lo pp é chacun
superficielle, p la densité des hydro­
un m odèle de l’é talem ent d 'u n e nappe
carbures et
d ’hydrocarbures axi-sym étrique.
L’é talem ent se décom pose en trois
régim es :
1) Inertie-gravité. pw est la densité de l’eau.
2) G ravité-friction. Les tem ps a uxquels se produisent
3) F riction-tension superficielle. les tran sition s d ’un régim e à l’autre
s ’o btiennent en é galisant les relations
D ans chacun des régim es, une re­
(3) et (4) d ’une part e t (4) et (5) d ’a u ­
lation entre la surface de la nappe et
tre part. A in si:
le tem ps, basée s u r une étude dim en-
sionnelle, a été donnée. R égim e ine rtie -gravité -» R ég im e g ra ­
vité-friction
La param étrisation de la friction uti­
lisée par Fay dans son m odèle fa it in­ t c = £21/4 g'- 1/4 k - 1/2 p 1/2 (7)
26 K. Maroihi et E. Deleersnijder

R ég im e ine rtie -gravité —> régim e fric- p aisseu r de la nappe peuvent être
tion-tension superficielle é crites (N ihoul, 1983):
ty = Q g '1/ 2 f 3/ 2 kp1/2 (8 ) d (hÛ)
+ V • (hU U ) =
3t
La figure 1 m ontre les phases su c­
cessives de l’é tale m e nt p ou r différents - h V ^ a + g 'h l+ ( x S - x g )p -1 (9)
vo lu m e s déversés.
La param étrisation de la friction 0h
+ V • (hU ) = - q (10)
a doptée conduit à des lois de l’é tale­ at
m ent qui d iffèrent de ce lles données
par Fay (1969), m ais qui so n t en ex­ où h est l’é paisseur de la nappe, U
cellent a ccord a ve c les observations. la vite sse horizontale m oyenne s u r la
nappe, Pa la pression atm osphérique
(supposée constante sur la région de
LE MODELE la nappe, xo (xo) la tension de surface
dans la nappe à l’interface de celle-ci
Le m odèle a été d écrit dans une p ub li­ e t de l’a ir (de l’eau), q est le taux de
cation antérieure (N ihoul, 1983). Les variation de l’é paisseur de la nappe
é quations du m odèle intég ré es sur l’é ­ p ar évaporation, entraînem ent des

3-5

Inertie-Gravité
Friction-Tension
Superficielle

ni Gravité-Friction

3 .5 4 .0 4 .5 S.O S .5 1 ° 9 .0*

F ig . 1. Variation du rayon de la nappe en fonction du tem ps pour différents volum es déversés.

 Simulation mathématique des nappes d’hydrocarbures 27

hydrocarbures dans la colo nn e d ’eau où Ua = x | / k (14) e st la vitesse d ’ad-

et de l’eau de m er dans la nappe (on vection due au v e n t et à l’action des
supposera q > 0 , et se p la çan t dans co uran ts induits localem ent par celui-
le cas le p lus défavorable, on prendra ci, Uw la vitesse du courant m arin ré­
q = 0), p, pw et g ' o nt été défini pré­
gional (m arée, tem pête, circulation
cédem m ent.
résiduelle),
La fo rce norm ale résultant de la
tension supe rficie lle e st négligée m ais | j = (Ywa - 78a —Y w a ) 7k (15)
la force tangentielle e st prise en À = g ' p k ~1 (16)
com pte. Ainsi (N ihoul, 198 3):
x § - Y& = t | - k(U- LT) Le transport dû au vent est para-
m étrisé, d ’a près les observations,
+ (Ywa - Toa - Tow) V o (12) com m e s u it (N ihoul, 1983):
où x i e st la tension due au vent, k le
Il u all = 10-3/2 II VII (17)
co efficie nt de friction (à l’interface dé- —>

fo rm a ble et perm éable entre les où V est la vite sse du v e n t à 10 m

hydrocarbures et l’eau, on a do pte une de la surface. Le ve cte u r Ua est dévié
représentation linéaire de la tension à d roite de V d ’un a ng le :
Xo, en fonction de la vite sse m oyenne
a = 4 0 - 8 II V II172
U, e.g. P edersen, 1972), Uw la vitesse
pour0< Il VII < 2 5 m / s (18)
du courant sous la nappe, ywa la te n ­
sion superficielle à l’interface eau-air, a = 0 pour II VII > 25 m / s (19)
Yoa (Yow) la tension supe rficie lle à l’in­
P o ur les expériences en M er du
terface h yd roca rb ure s-a ir (hydrocar-
N ord, les courants m arins sont déter­
bures-eau), 0 une fo nctio n traduisant
m inés à p a rtir du m odèle h ydrodyna­
les va riatio n s d e s te nsio ns superfi­
m ique du GHER.
cielles en fonction du m élange pro­
g ressif d e s hydrocarbures e t de l’eau L’effet des forces de tensions su ­
de leurs valeurs yoa et y8w au centre perficielles e st relié aux g ra die nts de
de la nappe (hyd ro carb ures purs), à tensions superficielles. C es derniers
leurs valeurs Ywa et O au bord de résultent de l’entraînem ent de l’eau de
celle-ci. m er dans la nappe et produisent un
pro cessu s d’advection des h ydroca r­
Le régim e inertial e st relativem ent
bures qui se superpose à celui asso ­
court et ne concerne que les prem iers
cié aux ve nts et aux courants. En
instants de l’étalem ent. En le négli­
adm ettant que la nappe répond élas-
geant et en co m b in a nt les équations
tiq ue m e nt à la tension superficielle
(9) et (10), on o btie n t com pte tenu de
l’équation (12) : (B a tche lor 1967, C hia-S hun Yih 1968,
N ihoul 1983), la form e appropriée
U + V • ( ( O + Uc + p V 0 )h) p ou r 0 est alo rs (N ihoul 1983, 1984):

= Ä V • (h2 V h) (13)
28 K. Maroihi et E. Deleersnijder

où r e st la distance au centre de la La v a le u r de B est déduite de la

nappe et s son rayon à l’instant t dé­ condition de co nsta nce de volum e, q
fin i par étant nul, l’équation ( 10) du m odèle
se réduit à :
S = te s2
S est la surface de la nappe à l’instant t.
ƒ hdS = £2 Vt (28)
D ans d e s co nd itio ns de co uran t et
de vent relativem ent ca lm es, l’advec-
tion a sso cié e au ve n t et a u x courants où £2 e st le vo lum e initialem ent dé­
e st sp atia le m en t co nsta nte, et dans versé. En com binant les équations
un systèm e d ’a xes a ssocié au centre (29), (28) et (23), on trou ve:
de la nappe, on peut su p p o se r en pre­
3£2 (29)
m ière approxim ation, une sym étrie cy­ B=
2 n
lindrique de cette dernière. L'équation
(13) s'é crit alo rs en coordonnées po­
laires: L’intégration de l’é quation (26)
d on ne :
3h 2 1 r2 1 1 3h
2 P T? h Xrh2
3t + r 3r r 3r dr
S - - a rctg (aS) = 4 n pt (30)
(2 1 ) â
C ette équation a dm et des solutions
/2
a nalytiques de la fo rm e : 2 pk
ou a =
h = B s - 2/( ii) (2 2 ) 3£2 gp
/

où q = r/s (23) S = us2

et où Be st une constante qui sera
déterm inée par la suite. La fo nctio n ƒ L’équation (30) donne le rayon de
doit sa tisfaire à la co n d itio n : la nappe en fo nctio n du tem ps. Son
com p orte m e nt asym ptotique e st en
ƒ (1) = 0 (24)
p arfait accord a vec les p rédictions de
En substituant la fo nctio n (22) dans l’analyse dim ensionnelle (N ihoul
l’équation (21), on o btie n t deux équa­ 1 98 3):
tio ns d ifférentielles p ou r ƒ et s :
d as2 < 1 (ré gim e gravité-friction)
= o (25)
dt dn
s - X1/6 £21/ 3 11/6 (31)
ds (26)
s ~ - ÀB2s-4 = 2 p
dt as2 > 1 (régim e friction-tension su p e r­
ficielle)
Tenant com pte de la co ndition (24),
l'équation (25) d on ne : , 1 /2
(32)
♦ 1 /2
m ) = (1 - n 2)172 (27) s - 2
 Simulation mathématique des nappes d'hydrocarbures 29

COMPARAISON m ent de la tache centrale est

AVEC LES OBSERVATIONS com m an dé p a r la gravité, ta nd is que
la couche m ince périphérique e st ra­
pidem ent sujette aux fo rce s de te n ­
En 1981, les va le u rs obse rvé es par sion superficielle.
Kuipers lors d’une expérience de dé­
C eci e st vérifié s u r la fig ure 3 qui
ve rse m e n t au large des cô tes hollan­
m ontre les va riatio n s tem porelles des
daises (3 m3) se disp ose nt fidèlem ent
su rfa ces de la tache centrale (...) et
s u r la d roite ca ra ctéristiqu e de l’é ta­
de la nappe to ta le (xxx) observées par
lem ent dans la phase friction-tension
te m p s calm e (vitesse du v e n t 2 ,4 m/s,
(fig- 2). H ollinger et M ennella).
U ne série d ’e xp é rie nce s de déver­
se m e nts co nd uite s entre août 71 et Lors de l’expérience Archim èdes
août 72 au large d e s cô tes de Virginie (oct. 21-22, 1983) conduite en M er du
o nt toutes m is en évidence une confi­ N ord, les données de la télédétection
guration ty p iq u e : une région centrale n’é ta ie n t pas su ffisantes p ou r c o m p a ­
d ’é pa isse ur de 1 m m ou plus, co nte ­ rer les prévisions du m odèle a u x o b ­
nant plus de 90 % du vo lum e déversé servations. L’a pplication du m odèle a
p ou r une surface de m oins de 10 % , été lim itée à l’évaluation du transport
entourée d ’une nappe plus vaste, d ’é­ général de la nappe (fig. 4). A la fin
p aisseu r p lus p etite de plusieurs o r­ des cinq prédictions rapportées s u r la
dres de grandeurs (H o llin g e r et figure 4, les écarts entre le s positions
M ennella 1973). L’interprétation de calculées et celles obse rvé es se si­
ces obse rva tion s suggère que l’é tale­ tu e n t entre 0,33 km et 2,5 km. C onsi-

TE N S IO N

10- 10 -

z:

10 - 10-J
G R R V IT E -E R IC T IO N

10: 10 '

T f S EC J

Fig. 2. Variation de la superficie de la nappe Fig. 3. Variation de la surface de la tache centrale

dans le cas d ’un déversem ent expérim ental de (...) et de la surface totale de la nappe (xxx) lors
3 m3. Com paraison entre le m odèle et les ob­ d'une expérience de déversement de 3.38 m 3.
servations. Comparaison entre modèle et observations.
30 K. Maroihi et E. Deleersnijder

l’interpolation d e s données de la télé­

détection.
(5 2 ° 3 1 '0 a " N, 0 3 °4 2 '3 0 " E)
La connexion e n tre les m odèles
m athém atiques e t la télédétection
1: 6 . 0 0 G .M .T peut jo u e r un rôle m ajeur dans l’am é­
2: 1 0 . 0 0 G.M.T
+5
lioration et le renforcem ent des tech­
3: 1 2 . 0 0 G.M.T
\ niques de nettoyage.
4: 1 4 . 0 0 G.M.T
+

5: 1 7 .0 0 G .M .T I. X
i 5,
2x
+ 2 1 km
4: REFERENCES
3x
f /
í '+ + p o s itio n s calculées
3
x p o s itio n s observées
Batchelor G .K., 1967. An introduction to
fluid dynamics. Cambridge Univ.
F ig . 4. Trajectoire de la nappe dans le cas de Press, 615 pp.
l’expérience Archim èdes. Com paraison entre le
m odèle et les observations. Deleersnijder E., 1991. Revisting Nihoul’s
model for oil slicks transport and
spreading on the sea. Submitted to
d éran t que les co uran ts m arins pré­ Ecological Modelling.
dits, et non ce ux localem ent observés, Fay J.A., 1969. In: Oil on the sea, D.P.
o nt été utilisés dans le m odèle, il est Hoult ED, Plenum Publ., New York,
é vident q ue les résultats concordent 53-63.
rem arquablem ent a vec les obse rva ­ Garrett W.D. & Barger W.R., 1970. Envi­
ron. Sec. Technol. 4 : 123-127.
tions.
Gillot R.M., 1983. North Sea campaign on
remote sensing o f oil spill. Oct. 83,
ACPH W DI, 7 pp.
CONCLUSION Hollinger J.P. & M ennella R.A., 1973.
Measurements o f the distribution and
volume of sea surface oil spills using
Le m odèle pré sen té ci-d essu s est ca­ multifrequency microwave radiometry.
pable de te n ir com pte sim ultaném ent Nat. Res. Lab., W ashington, Rep.,
de la gravité, de la friction , de la ten­ 7512, 5 pp.
sion superficielle, des effets du vent Hoult D.P., 1972. Oil spreading on the sea.
et des courants. Ann. Rev. Fluid Mech, 4 : 341-368.
Kuipers H.D., 1981 . A simulation model for
La com paraison entre le s prévi­ oil slicks at sea. North Sea Directorate,
sions du m odèle et les résultats des Department of Public W orks, The Ne­
e xpériences sus m entionnées m ettent therlands, 180 pp.
en é vidence sa capacité à décrire des Lee R.A.S., 1971. A study of the surface
situ atio ns réelles en mer. tension controlled regime o f oil spread.
SM Thesis MIT.
Les résultats du m odèle peuvent
Massin J.M. & Vollmers R.R., 1979. Télé­
être utilisés d ’une m anière dynam ique détection des rejets d'hydrocarbures
dans l’interprétation, l’extrapolation et dans les eaux côtières. Ministère de
 Simulation mathématique des nappes d’hydrocarbures 31

l’Environnement et du Cadre de Vie, Pedersen F.B., 1972. Gradually varying

Paris, 343 pp. tw o-layer stratified flow. Proc. Ann.
Nihoul J.C.J., 1983. A non-linear mathe­ Soc. Civil Eng., HYI, 9 8 : 257-268.
matical m odel fo r the transport and Smith J.E., 1968. Torrey Canion pollution
spreading of oil slicks. In : Modelling and marine life. C am bridge Univ.
the fate and effect of toxic substances Press.
in the environment, Proc. Coll. Inter­ Venkatesh S., Neralka V.R., Jarvic E.C.,
nat. Soc. o f Ecological Modelling, June Jessup R.G. et De Lorenzis B., 1981.
6-10, 1983, Copenhagen. A numerical model for predicting the
Nihoul J.C.J., 1984. Modèle de transport m ovem ent of oil slicks. In : Mechanics
et d ’étalement des nappes d ’hydrocar­ o f oil slicks, ENCP Ed., P aris: 219-
bures et comparaison avec les obser­ 233.
vations. Société Hydrotechnique de Warluzel A. et Benque J.P., 1981. Un
France, XVIIIo Journées de l’Hydrauli­ modèle m athématique de transport et
que, Marseille, 11, 12 et 13 septembre d ’étalement d'une nappe d ’hydrocar­
1984, Rapport n° 10 IV 10.1-IV bures. In: Mechanises o f oil slicks,
10.4 pp. ENPC Ed., P aris: 199-211.