Numerical simulation of a stratigraphic model at geological time scales

Description

Niveau: Secondaire, Lycée, Terminale
Numeri al Simulation of a Stratigraphi Model Veronique Gervais and Roland Masson Institut Fran ais du Petrole 1 et 4 av. de Bois Preau 92852 Rueil Malmaison Cedex, Fran e. email : roland.massonifp.fr Abstra t. A multi-lithology diusive stratigraphi model is onsidered whi h simu- lates at large s ales in spa e and time the inll of sedimentary basins governed by the intera tion between te toni s displa ements, eustati s variations, sediments supply and sediments transport laws. The model a ounts for the mass onservation of ea h sediment lithology resulting in a mixed paraboli , hyperboli system of PDEs for the lithology on entrations and the sediments thi kness. It also takes into a ount a limit on the ro k alteration velo ity modelled as a unilaterality onstraint. To obtain a robust, fast and a urate simulation, fully and semi-impli it nite volume dis retization s hemes are derived for whi h the existen e of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality onstraint, and pre onditioning strategies are dened for the solution of the linear system at ea h Newton iteration. They are based on an algebrai approximate de oupling of the sediments thi kness and the on entrations variables as well as on a proper pre onditioning of ea h variable.

• sediments


• exists many


• single lithology


• fun tion


• stratigraphi model


• velo ity


• pre onditioning methods

Sujets

• Première
• Masson
• Gervais
• Sédiment

Informations

• Publié par: profil-zyan-2012
• Nombre de lectures: 21
• Langue: English

Extrait

Sim
ulation
terms
in
Newton
of
These
a
of

tained
Mo
o
del
hand,
V
build


eronique
tro
Gerv
optimize
ais
ks
and

Roland
dening
Masson

Institut
ell.
F
Publishers.
r
w
an
stabilit


du
w
P
whic

as
etr
delling.
ole
v
1
On
et
a
4
determine
av.
ers
de
measuring
Bois
en
Pr
the

as
eau
describing
92852
10:54;
R
and
ueil
test
Malmaison
v
Ce
lateralit
dex,
e
F
m
r


gas),
e.
basins,
email
ma
:
of
r
a

and

the
A
geometry
m
disp
ulti-lithology
es
diusiv

e
of

at
mo

del
of
is
sedimen

time.
whic
are
h
erties
sim
y
u-
h
lates
tology
at
b
large
data,
scales
tologists
in
del
space
Kluwer
and
the
time
in
the
scalabilit
inll
on
of
real
sedimen
Keyw
tary
delling,
basins
metho
go
analysis,
v

erned
linear
b
er,
y
preconditioner,
the
1.
in
Man

h
b
energy
et
are
w
sedimen
een
understanding

is

issue
ts,
managemen
eustatics
natural
v
ell
ariations,
hallenge
sedimen
ph
ts
umerical
supply
o
and
ro
sedimen
erties
ts
geologists
transp
at
ort
t
la
t
ws.
eld
The
one
mo

del
the


ts
es
for
ground
the
are
mass
order

p
ation

of
geometry

la
h
presen
sedimen
the
t
ell-log
lithology
b
resulting
ysical
in
h
a
or
mixed
a
parab
ell,
olic,
in
h
of
yp
erties
erb
ts
olic
the
system
these
of
w
PDEs
data,
for
able
the
geological
lithology
the


trations

and
d
the
GervaisMasson.tex;
sedimen

ts
terms
thic
robustness,
kness.
y
It

also
t
tak
o
es
basin
in

to
ords:

mo
t
nite
a
olume
limit
d,
on
y
the
uni-
ro
y

t,
k
algorithm,
alteration
iterativ
v
solv
elo
blo

k
y

mo
ultigrid
delled
In
as

a
y
unilateralit

y
as

ater,
t.
(oil,
T
ore,
o

obtain
in
a
tary
robust,
the
fast
of
and
h

a
sim
jor
ulation,
to
fully
the
and
t

ground
nite

v
w
olume
as
discretization

sc
in
hemes
of
are
ysical
deriv
n
ed
mo
for
T
whic
determine
h
ground
the

existence
prop
of
and
stable
,
solutions
ha
is
e
pro
their
v
osal
ed.
w
Then,
main
the
yp
set
of
of
data.
nonlinear
the
equations
hand,
is
data
solv
the
ed
in
using
basin
a

Newton
w
algorithm
v
adapted
emitted
to
the
the
surface.
unilateralit
signals
y
pro

in
t,
to
and
the
preconditioning
osition
strategies
the
are

dened
the
for
of
the
ts
solution
y
of
at
the
t
linear
On
system
other
at
w

data
h
obtained
Newton
y
iteration.
ph
They
prop
are

based
as
on
y
an
resistivit

within
appro
giv
ximate
w

whic
of
are
the
terpreted
sedimen
terms
ts
sedimen
thic
prop
kness
of
and
sedimen
the


y
trations
w
v
Using
ariables
eld
as
as
w
ell
ell

as
sedimen-
on
are
a
to
prop
a
er
mo
preconditioning
of
of
basin


h
2005
v
A
ariable.
ademic
These
Printe
algorithms
in
are
Netherlands.
studied
21/04/2005;
and
p.12
V

eronique
del
rate,
dieren
Gerv
the
ais
y
and
(Granjeon
Roland
and
Masson
sedimen
its
t
inll
10:54;
scenario
diusiv
b
terms,
y
assumes
sedimen
the
tation/erosion
basin.
in
the
order

to
limiter



:
the
transp
sedimen
(Riv
tary
to

disc
of
and
the
the

of
e
dep
sedimen

ts

la
ailabilit
y
b
ers,
and
in
w
terms
(Eymard
of
h
paleo-bath
a
ymetry
the
,
1
dep

ositional
p
en
is
vironmen
in
t,
extended

to
osition
onate
of
with
the
formalization
sedi-
rigorously
men
,
ts,
are
grain
the
sizes,
fraction
p
at
orosit
enden
y

,
ts
...
h

tak
h
b
mo
the
dels
dible
are
basin
a
t
starting

p
h
oin
elo
t

of
ation
oil
uses
exploration
a
studies.
form
Ho
of
w
e.g.
ev
ts
er,
v