Calcul Et Analyse Des Attributs Sismiques
Calcul Et Analyse Des Attributs Sismiques
Calcul Et Analyse Des Attributs Sismiques
- 76 -
CAlCUL ET ANALYSE DES ATTRIBUTS
SISMIQUES
INTRODUCTION :
Lanalyse des attributs sismiques complexes de la trace sismique dbute par un traitement
en amplitudes prserves des donnes sismiques , afin de reprsenter les amplitudes relles
dues la rflectivit .
Les donnes sismiques doivent tre corriges des diffrents effets subis par le signal
dexcitation depuis lmission dnergie jusqu' lenregistrement . Bien quon nait pas
encore bien prouv toutes les applications des attributs sismiques , il en existe quelques uns
plus connus, comme lamplitude instantane , ses variations avec loffset, les vitesses
sismiques, la phase et la frquence instantanes, qui ont t utilises dans linterprtation
sismique et la dtection directe des hydrocarbures.
1)- NOTION ET CLASSIFICATION DES ATTRIBUTS
SISMIQUIES :
1)-1)- DEFINITION :
On appelle attribut sismique toute caractristique, qualitative ou quantitative, mesure,
calcule ou dduite partir des donnes sismiques, reprsentant ainsi tous les paramtres de la
trace complexe, les configurations gomtriques des vnements sismiques et leurs variations
spatiales. Ils sont lis aux informations de base des donnes sismiques : temps, amplitude,
frquence et attnuation. Les attributs offrent des reprsentations alternatives aux donnes
sismiques et mettent ainsi en relief les caractristiques ptrophysiques et gologiques des
paramtres tudis. Ces attributs ont t classifis dune faon qui permet dexprimer au
mieux leur utilit en sismique.
1)-2)- OBJECTIFS DE LUTILISATION DES ATTRIBUTS SISMIQUES :
Etude des caractristiques du signal sismique
Etude de la configuration gomtrique des vnements sismiques .
Interprtation structurale ,stratigraphique et lithologique des prospects.
Estimation des proprits des roches :
Attributs avant sommation et combinaison avec lAVO .
Chapitre : 05 Attributs sismiques 3D
- 77 -
Nous divisions les attributs sismiques gnralement en deux classes . Des attributs de la
classe I sont calculs directement partir des traces .( lenveloppe de trace et ses drivs ,
phase instantane et ses drivs , le facteur Q , etc) .
Les attributs de la classe II sont calculs sur des traces avec une amlioration du rapport
S/B.Tous les attributs de la classe I sont calculs dans la classe II .
Les attributs instantans reprsentent les variations instantanes de divers paramtres.
Le calcul de ce type dattributs se fait chantillon par chantillon .Lenveloppe de trace et ses
drivs , la frquence et la phase peuvent tre dtermins partir des traces complexes .
Plusieurs auteurs ont avanc leur propre classification, nous naborderons dans le
prsent travail, que la plus simple base sur les caractristiques des attributs.
1)-3)- LES DIFFERENTS TYPES DATTRIBUTS SISMIQUES :
1)-3)-1)- Attributs physiques :" Physical Attributes "
Ils sont dtermins partir des mesures sismiques lies directement aux caractristiques
et des estimations de paramtres physiques. Lenveloppe de la trace est proportionnelle au
contraste dimpdance acoustique ; les frquences sont lies aux phnomnes de dispersion et
dabsorption ainsi qu lpaisseur des bancs. Les vitesses instantane et moyenne sont
directement lies aux proprits des roches.
Par consquent, ces attributs sont le plus souvent utiliss pour la classification
lithologique et la caractrisation des rservoirs.
1)-3)-2)- Attributs gomtriques : " Geometrical Attributes "
Ces attributs dcrivent les relations spatiale et temporelle de tous les autres attributs.
Initialement, les attributs gomtriques ont t utiliss aider linterprtation stratigraphique.
Cependant, des expriences ont montr quils dfinissent aussi les caractristiques des
vnements sismiques et leurs relations spatiales, quantifient (mesurent) les traits qui peuvent
mener directement lidentification des modles de dpts, et de la lithologie.
1)-3)-3)-Attributs pre Stack :" Pre-stack Attributes "
Le pre-stack sous-entend que les donnes sismiques sont organises en points-miroir
communs, renfermant chacun toutes les traces juxtaposes. Ces donnes contiennent alors des
informations concernant la direction (azimut) et loffset.
Le calcul dans ce cas gnre des quantits colossales de donnes ; de l ces attributs ne sont
pas pratiques pour des tudes prliminaires.
1)-3)-4)-Attributs post Stack :" Post-stack Attributes "
Aprs sommation en points-miroir commun, les informations lies loffset et
lazimut seront perdues. Pour le calcul des attributs dans ce cas, les donnes peuvent tre
sommes en point miroir commun ou migres en temps.
Une meilleure observation de donnes, pour des investigations primaires, est ralise par les
attributs post-stack. Cependant, pour une tude minutieuse les attributs pre-stack peuvent tre
incorpors.
Chapitre : 05 Attributs sismiques 3D
- 78 -
Attributs sismiques
Attributs physiques
Reli la propagation de londe
sismique , lithologie
Attributs Gomtriques
Reli laspect gomtrique des
couches .
Attributs pr-stack
Calculs partir de la trace complexe dont la
partie imaginaire calcule par la
transformation de Hilbert .
Domaine dapplication : donnes avant
sommation
Attributs post-stack
Calculs partir de la trace complexe dont la
partie imaginaire calcule par la
transformation de Hilbert .
Domaine dapplication : donnes aprs
sommation
Attributs instantans
Calcul chantillon par chantillon , indique le
changement lattributs dans le temps et dans lespace
Exp : Pr-stack
Vitesse RMS des reflecteurs ,
Section sismique en ondes P,
Section sismique en S .
Post- stack
Frquence instantane ,
Phase instantane ,
Amplitude de lenveloppe ou amplitude
instantane .
Attributs de lenveloppe
Calcul de toute les enveloppes de la trace , piquer tous les
minimas des enveloppes , la distance entre deux minimas
adjacent reprsente londelette . Lattributs calcul au
maximums reprsente lattribut de lenveloppe .
Ces attributs indiquent les variations latrales des
caractristiques de londelette .
Fig.5.1
Classification des attributs
sismiques
Chapitre : 05 Attributs sismiques 3D
- 79 -
1)-3)-5)-Attributs instantans :" Instantaneous Attributes "
Ces attributs sont calculs partir de la trace complexe (analytique). Leur calcul se fait
chantillon par chantillon, reprsentant ainsi les variations instantanes de plusieurs
paramtres.
1)-3)-6)-Attrbuts dnodelette : " Wavelet Attributes "
Ces attributs sont, par dfinition, les attributs instantans calculs aux points maximaux
de lenveloppe. Leur reprsentation fournira des aperus sur les caractristiques de londelette
rflchie, do le nom dattribut dondelette.
Du moment quun seul attribut est reprsent sur la longueur de londelette, on doit donc
dfinir cette dernire. Londelette sismique simple ou compose (constitue de plusieurs
ondelettes) est suppose occuper la distance en temps entre deux minima denveloppe
adjacents. Aprs avoir calcul lenveloppe de la trace sismique, on cherche toutes les
positions dextremums. Les positions maximales seront utilises pour obtenir les attributs
instantans et assigner ces valeurs la zone temporelle entre les minima de chaque ct de la
position maximale.
Pour plus de prcision, les positions maximales de lenveloppe sont obtenues par
interpolation quadratique, et les attributs sont aussi calculs par interpolation en ces positions
interpoles.
2)- NOTION DE LA TRACE ANALYTIQUE :
La trace sismique peut tre traite comme tant la partie relle dune trace complexe
calcule sous les conditions usuelles.
La trace complexe (analytique) permet lunique sparation de lamplitude de lenveloppe et
linformation concernant la phase, ainsi que le calcul de la frquence instantane.
Soit la trace sismique enregistre T(t) donne par :
) ( cos ) ( ) ( t t E t T = (5.1)
O E(t) lenveloppe de la trace et ) (t la phase.
Considrons un vecteur perpendiculaire laxe du temps et dont la longueur varie
comme une fonction du temps. Ce vecteur subit une rotation, galement comme une fonction
du temps, autour de cet axe (Voir chapitre 4, figure 4.6).
La projection de ce vecteur en rotation, dans le plan rel donnera la trace sismique
conventionnelle T(t). Cette projection peut tre effectue sur le plan imaginaire, ce qui
donnera ce que lon appelle : la trace en quadrature T
Q
(t). T
Q
(t) peut sexprimer ainsi :
) ( sin ) ( ) ( t t E t T
Q
= (5.2)
Les techniques de la transforme de Hilbert (BRACEWELL, 1965) permettent de gnrer la
trace en quadrature partir de celle enregistre .
Chapitre : 05 Attributs sismiques 3D
- 80 -
Ainsi, la trace analytique est la fonction complexe de la trace sismique relle, exprime de la
manire suivante :
) ( ) ( ) ( t jT t T t T
Q A
+ = (5.3)
( ) ) ( sin ) ( cos ) ( ) ( t j t t E t T
A
+ = (5.4)
) (
) ( ) (
t j
A
e t E t T
= (5.5)
Le passage de londe sismique provoque un dplacement des particules de la subsurface
par rapport leur position dquilibre. Ainsi la trace sismique reprsente une mesure de la
vitesse de mouvement (avec des gophones) ou la variation de pression (avec des
hydrophones). De ce fait, elle peut tre considre comme une mesure de lnergie cintique.
tant donn que les particules rsistent au mouvement au moyen dune force lastique,
lnergie devient ainsi, emmagasine sous forme dune nergie potentielle. Comme les
particules se dplacent en rponse lbranlement sismique, il y a donc constamment change
dnergie entre cintique et potentielle. La trace en quadrature est considre comme une
mesure de lnergie potentielle.
3)- CALCUL DE LA TRACE EN QUADRATURE :
La trace en quadrature est obtenue par la transforme de Hilbert de la trace sismique
relle. Elle peut tre, galement, obtenue au moyen de lintgrale de Fourier (Voir chapitre
4, 4.5.2).
T
Q
(t) = TH [ T(t) ]
t
V t T t T
p Q
1 1
) ( ) (
=
= d
t
s
V t s
p
) (
1
) (
TRACE SISMIQUE EN QUADRATURE
(IMAGINAIRE)
TRACE SISMIQUE ENREGISTRE
(RELLE)
TRACE SISMIQUE
COMPLEXE
TEMPS
Fig. 5.2 Reprsentation de la trace analytique.
Chapitre : 05 Attributs sismiques 3D
- 81 -
Les traces relle et en quadrature doivent avoir le mme spectre damplitude, mais leurs
phases sont altre de 2 .
T
Q
(f) = -j sgn(f) . T(f)
Le filtre de Hilbert, dans la domaine temporel, est de la forme
n
1
:
=
0 0
0
) 2 ( sin
2
) (
2
n
n
n
n
n g
g(n) est la rponse du filtre de Hilbert.
Si ce filtre est correct du point de vue thorique, il nest pas cependant, adapt dans les
applications pratiques, cause de sa longueur infinie et des distorsions excessives causes par
les discontinuits dans les grandes valeurs du spectre de frquence.
Aussi, la trace en quadrature obtenue par la transforme de Hilbert na pas le mme
spectre damplitude que la trace relle.
Pour pallier ce problme, on utilise pour le calcul un filtre passe bande convenable, comme
celui de BUTTERWORTH :
un filtre phase nulle pour gnrer la partie relle, et un autre, de caractristique et de bande
passante identiques que celles du premier mais phase de 90, pour gnrer la partie
imaginaire.
Ainsi, parties relle et imaginaire seront des paires de transforme de Hilbert. De cette
manire on assurera un spectre damplitude identique pour les deux parties, et on vitera les
distorsions excessives.
4)- BASES MATHMATIQUES DES ATTRIBUTS SISMIQUES ET
LEURS APPLICATIONS :
Le calcul des attributs est fait partir dune section sismique en amplitudes prserves
somme en point miroir commun et migre en temps.
Nous donnerons, dans la mesure du possible, pour chaque attribut sa formulation
mathmatique, et on indiquera sa relation directe avec les proprits physiques de la
subsurface.
4)-1)-ATTRIBUTS PHYSIQUES (INSTANTANES) :
4)-1)-1)- Amplitude instantane :
Lamplitude instantane est lattribut particulirement important attache directement aux
conditions gologiques , qui reprsente une mesure des amplitudes indpendamment de la
polarit .
Lamplitude instantane dune trace sismique permet indirectement une attnuation forte des
bruits et le renforcement de la rflexion utile ( Voir figufe.5.3 ) .
Chapitre : 05 Attributs sismiques 3D
- 82 -
0 0
Soit la trace analytique T
A
(t) dfinie par :
T
A
(t) = T(t) + j T
Q
(t)
Lamplitude instantane, dite aussi enveloppe de la trace sismique, est dfinie comme le
module de la trace complexe :
2 2
) ( ) ( ) ( t T t T t E
Q
+ = (5.6)
E(t) reprsente lnergie instantane totale, et elle varie de zro lamplitude
maximale de la trace.
Comme indiqu dans lquation (5.6), lenveloppe est indpendante de la phase, car elle peut
atteindre son maximum en un point autre que le pic ou le creux de la trace relle, et elle est
directement lie au contraste dimpdance acoustique. Lenveloppe est un attribut physique
qui peut tre utilis comme un discriminateur efficace pour les caractristiques
suivantes :
Lenveloppe reprsente le contraste dimpdance acoustique, et donc de la rflectivit,
T
Q
(t)
T(t)
E(t)
0
Fig. 5.4
Reprsentation de
lenveloppe de la trace
analytique.
temps
temps
Trace sismique enregistre
Amplitude instantane
Fig.5.3
Reprsentation de lamplitude instantane dune
trace sismique relle
Chapitre : 05 Attributs sismiques 3D
- 83 -
Indicateur des anomalies damplitudes (Bright spots) ,
discordances, changement lithologique,
accumulation de gaz,
changement majeur de lenvironnement de dpot ,
corrlation spatiale des paramtres ptrophysiques ,
Indicateur des fracturations ,
Les variations latrales de facis ,
Leffet des couches minces Tunning effect .
4)-1)-2)- Phase instantane :
La phase instantane est lattribut indpendant de lamplitude instantane , Elle est dfinie
comme largument de la trace complexe :
(
=
) (
) (
arctan ) (
t T
t T
t
Q
(5.7)
La phase instantane est reprsente en degr et dune manire continue, en utilisant une
palette de couleurs pour laquelle au 0 est associ la couleur jaune, +180 la couleur rouge
et -180 (360) la couleur bleue.
Cet attribut est li la phase de propagation du front donde sismique. Puisque ce dernier,
pour la plupart du temps, est dfini comme des lignes de phase constante.
Le principale avantage de la reprsentation de la phase instantane est que celle i peut tre
reprsente aisment dans les zones de faibles amplitude ainsi que les zones de forte
amplitude .
Cet attribut sismique est galement un attribut physique, et peut tre :
un bon indicateur de la continuit latrale et des discontinuits ,
un indicateur des rflexions horizontales qui rsultent du contact gazeau ou gazhuile
(Flat-spot),
utilis pour le calcul de la vitesse de phase,
un indicateur de discordances, mais il nest pas le meilleur,
utilis pour le calcul de la frquence instantane et lacclration,
identifi les limites de squences ,
T
Q
(t)
T(t)
(t)
0
Fig. 5.5
Reprsentation de la phase
instantane de la trace
analytique.
Chapitre : 05 Attributs sismiques 3D
- 84 -
il permet daffiner les lments structuraux et de dtailler les lments stratigraphiques.
4)-1)-3)- Frquence instantane :
La frquence instantane est un autre attribut qui peut tre attach directement aux
cnditions gologiques et qui reprsente par la variations de la phase . Autrement dit, la
frquence instantane dsigne la vitesse de rotation du vecteur de Fresnel
) (t j
e
.
La frquence instantane est une mesure chaque chantillon de la frquence dominante dans
la section .
Lexpression mathmatique donnant la frquence instantane est :
dt
t d
t f
i
) (
2
1
) (
= (5.8)
Soit la trace analytique dfinie par :
) (
) ( ) (
t j
A
e t E t T
=
En prenant le logarithme des deux cts de lquation, on trouve :
) ( ) ( ln ) ( ln t j t E t T
A
+ =
La drive par rapport au temps donne :
dt
t d
j
dt
t E d
dt
t T d
A
) ( ) ( ln ) ( ln
+ = (5.9)
o le terme
dt
t d ) (
reprsente la frquence instantane.
En dveloppant la drive de la trace complexe, on aura :
( )
2 2
) ( ) (
) ( ) ( ) ( ) (
) ( ) (
) ( ) (
) (
) ( ) ( ln
t T t T
t jT t T t jT t T
t jT t T
t jT t T
t T
t T
dt
t T d
Q
Q Q
Q
Q
A
A A
+
|
.
|
\
|
+
=
+
+
=
2 2
) (
2 2
) (
) ( ) (
) ( ) ( ) (
) ( ) (
) ( ) ( ) (
) (
) (
t T t T
T t T t T t T
j
t T t T
T t T t T t T
t T
t T
Q
t Q Q
Q
t Q Q
A
A
+
+
+
(5.10)
Des expressions (5.9) et (5.10) on peut crire :
2 2
) (
) ( ) (
) ( ) ( ) ( ) (
t T t T
T t T t T t T
dt
t d
Q
t Q Q
+
La frquence instantane aura pour expression donc :
(
(
(
=
2 2
) ( ) (
) (
) ( ) (
) (
2
1
) (
t T t T
t T
dt
t dT
dt
t dT
t T
t f
Q
Q
Q
i
(5.11)
Chapitre : 05 Attributs sismiques 3D
- 85 -
La frquence ainsi dfinie est inapproprie ponctuellement mais est considre plutt
que comme une moyenne sur un intervalle.
La frquence instantane est lie la propagation des ondes et au type de dpt ; par
consquent, comme un attrbut peut tre un indicateur :
de la prsence des hydrocarbures ( des anomalies de basses frquences ),
de limite dunits sdimentaires ou structurales comme les failles,
de zones de fractures ( zones basses frquences ) ,
du rapport argile/ grs ,
des rflexions chaotiques .
Il est noter encore que de basses frquences sont souvent observes au niveau des
rflecteurs situs juste au-dessous dun gisement de gaz. Ceci a priori, a deux explications :
un gisement de gaz joue le rle dun filtre passe-bas,
la vitesse de propagation dans un rservoir gaz est faible. La chute de vitesse est
souvent accompagne dune absorption des hautes frquences, ce phnomne est rencontr
souvent dans les couches trs profondes , ce qui peut nous inciter interprter ces couches
comme tant des gisements de gaz ..
La frquence est aussi sensible aux changements des squences des lits (bedding
sequences) ; elle est donc utile pour indiquer o a eu lieu un changement stratigraphique.
4)-1)-4)- Amplitude normalise :
Lamplitude normalise reprsente par dfinition (TANER, 2000) le cosinus de la phase
instantane. La reprsentation de cet attribut contient donc tous les dtails de la section phase
instantane et indique la continuit latrale.
La reprsentation en amplitudes normalises est utilis pour la dtermination de la direction
de bonnes continuits latrale , et ne contient gure dinformation sur lamplitude de la trace
enregistre.
4)-1)-5)- Impdance acoustique relative :
Si on suppose que la trace sismique a t traite dune manire tre dbarrasse des
diffrents bruits multiples, et contenir les ondelettes phase nulle. Partant de ces
hypothse , cette trace devrait donc, reprsenter la srie des coefficients de rflexions
(rflectivit). Elle peut tre exprime par :
Z
Z
V
V
t T
=
=
2
1
2
1
) ( (5.12)
) ln(
2
1
) ( Z t T = (5.13)
Par intgration de la trace phase nulle, on obtient une estimation du log naturel
dimpdance acoustique :
=
t
dt t T Z
0
) ( 2 ) ln( (5.14)
Limpdance naura pas de valeurs absolues et la section somme est une estimation de la
rflectivit offset nul, do le nom dimpdance relative.
Chapitre : 05 Attributs sismiques 3D
- 86 -
Limpdance acoustique relative est un attribut physique, puisquelle reflte les proprits
physiques du contraste, et peut tre utilise pour des procdures de calibrage.
Cet attribut indique :
des surfaces de discordances,
des contrastes dimpdance acoustique apparente.
4)-2)-ATTRBUTS GEOMTRIQUES :
4)-2)-1)- Attribut cohrence cube :
Le traitement en cohrence cube est une mesure mathmatique de la ressemblance
(Similarity) un indicateur de discontinuits qui calcule des mesures visuellement
interprtables de cette dernire do le nom cohrence.
Les mesures de la ressemblance (Similarit) sont un des plus importants outils dans le
traitement sismique. Dans le domaine de l'exploration sismique, le traitement de la cohrence
cube permet lquipe d'exploration et de production de vrifier leurs interprtations. Cest
un processus numrique qui n'exige pas un horizon interprt. Le traitement de la cohrence
cube est un outil pour les donnes sismique 3D.
a) - Calcul de lattribut cohrence cube :
Les divers algorithmes avancs pour le calcul de la cohrence ont rapidement progress
depuis le dveloppement, par Bahorich et Farmer (1995), de lalgorithme original de la cross
corrlation (C1). Bien quil soit efficace du point de vue calcul, l'algorithme original de la
cross corrlation n'a pas donn de bon rsultats avec de donnes bruites.
Il a t remplac par semblance (C2) (Marfurt et al., 1998). Le concept de la cross
corrlation de 2 traces dans la direction inline et crossline a t gnralis a celle de multi-
trace en 3D sur une fentre d'analyse de j traces, travers la construction mathmatique d'une
matrice de covariance j x j. Les valeurs de la cohrence sont obtenues en calculant la
semblance avec plusieurs tests des couples pendage/azimut partir de la matrice de
covariance.
L'algorithme EigenStructure Modifi (Marfurt et al. 1999) pose le problme du pendage
apparent local et fait appel un algorithme hybride pour calculer la Cohrence le long des
plans inclins.
Plus tard en 1999 est venu haute rsolution caractristique (High-Resolution Eigen ) , avec
Attnuation Structurelle (appel quelquefois Eigen Modifi avec des valeurs de la cohrence
affines'). La sortie (output) de cet algorithme est un cube 3D des coefficients de cohrence
amliors (accentus).
b) -Calcul de la fonction de cohrence:
Soit la trace sismique T(t) dfinie par :
) ( ) ( ) ( ) ( t b t s t r t T + = .
Chapitre : 05 Attributs sismiques 3D
- 87 -
Dans cette expression, r(t) et b(t) (respectivement log de rflectivit et bruit) sont des
fonctions alatoires indpendantes, ce qui fait que leurs autocorrlations sont des fonctions
impulsionnelles et leur intercorrlation est nulle.
Ainsi, en appliquant ses proprits, lautocorrlation de T(t) scrit :
) ( ) ( ) (
, , ,
+ =
b b s r s r T T
) ( ) ( ) ( ) (
, , , ,
+ =
b b s s r r T T
) ( ) ( ) ( ) (
, , ,
+ =
b b s s T T
k
) ( ) ( ) (
, , ,
+ =
b b s s T T
k
Soit deux traces T
1
(t) et T
2
(t) voisines ayant le mme contenu informatif mais de bruits
alatoires diffrents. Les autocorrlations de ces deux traces :
) ( ) ( ) (
1 1 1 1
, , ,
+ =
b b s s T T
k
) ( ) ( ) (
2 2 2 2
, , ,
+ =
b b s s T T
k
Leur intercorrlation est donne par :
) ( ) (
, ,
2 1
=
s s T T
k
et comme ) (
,
s s
est une fonction paire :
) ( ) (
1 2 2 1
, ,
=
T T T T
En a dmettant que lnergie moyenne du bruit est constante sur des traces voisines, alors :
) ( ) (
2 2 1 1
, ,
=
b b b b
et ) ( ) (
2 2 1 1
, ,
=
T T T T
et si on prend les chelles telles que k = 1, les expressions seront alors :
) ( ) (
, ,
2 1
=
s s T T
et ) ( ) ( ) (
, , ,
+ =
b b s s T T
La TF dune autocorrlation, dite spectre de puissance de la fonction, reprsente bien son
nergie dans le domaine spectral. On dfinira la fonction de cohrence par :
2 2
2
, ,
,
, ,
, ,
) ( ) (
) (
) ( ) (
) (
) (
2
1
2 2 1 1
1 2 2 1
f B f S
f S
f f
f
f C
B B S S
S S
T T T T
T T T T
+
=
+
=
|
.
|
\
|
=
Cet attribut a pour objectife :
La dtermination des types de vitesse ,
Suivre les vnements rflchis ,
La mesure de rapport signal sur bruit,
Chapitre : 05 Attributs sismiques 3D
- 88 -
De rvle les formes des rflecteurs de la sub-surface ; tes que : biseaux , les discordances
et les limites du chenal , et traits sdimentaires subtils qui sont difficiles interprter sur
du volume sismique tridimensionnel ,
De dvelopper la continuit latrale des vnements sismiques et de classer le type de
discontinuit ,
De classer les terminaisons latrales des rflexions comme des Top-Lap et On-LAP etc,
pour tre utiles dans l'interprtation stratigraphique .
La cohrence prsente un outil puissant pour l'identification des discontinuits structurales
et stratigraphiques dans le plan horizontal , car la conservation de la rsolution
horizontale est assure ,
L'utilit de cet attribut est la comparaison avec les autres attributs ,
Les valeurs de la cohrence annoncent galement des changements dans la rponse
sismique due la lithologie , ou due la prsence dun fluide saturant , la permabilit ,
et aux proprits lastiques .
La cohrence cube peut donner une estimation de lextension latrale des plans des failles
.(Dans le processus du calcul de la cohrence pour un volume des donnes 3D , les failles
deviennent numriquement spares des donnes qui les entourent. En d'autres mots , les
failles engendrent des surfaces de faible cohrence . En employant les logiciels de
visualisation , les plans failles peuvent tre observs en trois dimensions de plusieurs
points de vue , bien qu'aucun de plan de faille n'ait t prenregistr ) Bahorich TEL 1995
.
4)-2)-2)- Zones de discordances :
Lalgorithme de cohrence est utilis pour la dtection de discordances possibles. En
premier lieu sont localiss les pics de lenveloppe, ce qui reprsente en sismique un interface,
donc une surface de discordance ventuelle. Sa grandeur est fonction du degr du contraste
dimpdance. Ensuite, on procde au calcul de contrastes additionnels. Pour se faire, on utilise
les valeurs de plusieurs attributs comme composants des vecteurs attributs. Trois vecteurs
sont confectionns des deux cts de linterface. Les composantes sont alors les trois
attributs : pendage de similitude maximale (dip of maximum similarity), frquence instantane
et similitude et variation de pendage (similarity and dip variance).
Les vecteurs seront normaliss, donc leurs grandeurs seront comprises entre +1 et 1. A la
valeur +1 correspond une grande similitude, indiquant ainsi quil ny a pas de diffrence des
deux cts de linterface du point de vue pendage, frquence et continuit. De faibles valeurs
indiquent un certain degr de diffrence. Le facteur de dissimilitude est form en soustrayant
1 du vecteur normalis, qui est multipli par la grandeur du pic de lenveloppe et sera en
sortie (output) comme indicateur de zones de discordance. Les grandes valeurs de loutput
reprsentent des degrs levs de contraste.
4)-2)-3)- Indicateur dargilosit :
Lindicateur dargile est un attribut gomtrique qui combine dautres attributs
sismiques pour dtecter de possibles zones dargile dans les roches clastiques. Les argiles
peuvent tre identifis par leur configuration gomtrique comme de bancs minces, parallles
avec une continuit latrale leve. Cet attribut utilise les hautes frquences instantanes
comme indicateur de bancs minces, les faibles valeurs de similitude maximale comme
indicateur de bancs parallles et la similitude et ses variations comme indicateur de la
continuit latrale. Les valeurs importantes indiquent une forte possibilit doccurrence
Chapitre : 05 Attributs sismiques 3D
- 89 -
dargile, tandis que les faibles indiquent des possibilits de lithologie autre que largile,
comme le sable ou les carbonates.
4)-2)-4)-Polarit apparente:
Alors que la mesure des attributs sismiques dpend de la qualit des donnes et de la
qualit denregistrement et du traitement, la mesure de la polarit apparente est spcialement
sensible la qualit des donnes.
Lanalyse de la polarit suppose des rflexions singulires, des ondelettes phase nulle et
labsence dambigut due linversion de la phase. Cependant, puisque la plupart des
vnements rflchis sont composs de plusieurs rflexions, la polarit donne une mauvaise
corrlation avec les coefficients de rflexions, do le nom de polarit apparente.
La polarit apparente fait, parfois, une distinction entre deux sortes de point brillant.
Associs aux accumulations de gaz dans les sdiments clastiques, ces points brillants ont
habituellement de plus faible impdance acoustique que les roches encaissantes, et prsentent
une polarit ngative pour les rflexions au top du gisement, cest ce quon appelle bright
spots, et une autre positive pour les rflexions la base (contacte gaz-huile ou gaz-eau), cest
ce quon appelle flat spots.
4)-3)- ATTRIBUTS DE SURFACES :
4)-3)-1)- Attribut pendage - azimut : (Dip-Azimuth)
Conu pour expliquer le fait que quelques failles apparaissent clairement sur la carte de
pendage mais pas sur celle de l'azimut, ou vice versa.
Montrer des aspects de pendage et dazimut en mme temps en combinant leurs valeurs.
Variations dans les couleurs indiquent lazimut du pendage gnralis dans les quatre
quadrants de la carte de lazimut.
Variations dans lintensit indique la magnitude (valeur) du pendage.
a)-Magnitude du pendage : (Dip magnitude)
La magnitude du gradient du temps dun plan qui relie chaque chantillon d'un horizon
interprt avec deux chantillons adjacents dans les directions orthogonales, exprim en
millisecondes par distance.
Magnitude du pendage =((dt /dx )2 +(dt/dy )2 )1/2
O : x : est la direction Inline ,
y : est la direction Xline .
b)- Azimut : (Azimuth)
La direction du pendage dun plan qui relie chaque chantillon d'un horizon interprt
avec deux chantillons adjacents dans les directions orthogonales, exprim en degr.
Chapitre : 05 Attributs sismiques 3D
- 90 -
Azimut =arctan ((dt/dx) / (dt/dy))
5)- COLORIMETRIE ET INTERPRETATION :
Linterprtateur se trouve souvent confront une masse de donnes importante quil
doit assimiler. De plus, chaque attribut reprsenter implique des valeurs associes chaque
chantillon. Assimiler galement une telle quantit de donnes pose un problme majeur.
Lintroduction des couleurs dans la reprsentation, peut aider donner de bonnes
significations aux mesures et leurs interdpendances. En effet, la couleur aide
linterprtateur en ajoutant dautres perspectives la section sismique.
La reprsentation est faite de manire associer une couleur chaque valeur ou une
gamme de valeurs relativement un attribut. Ceci permettra de mettre en vidence mme les
faibles variations difficilement dcelables sur une section sismique conventionnelle.
6) - ATTRIBUTS SISMIQUES ET LANALYSE LITHOLOGIQUE:
Lefficacit des attributs sismiques dpend de la capacit de dterminer, parmi le groupe
dattributs choisis, ceux qui portent le plus dinformations sur la subsurface. Comme dans le
cas de plusieurs problmes mathmatiques, quelques attributs peuvent tre ncessaires,
dautres suffisants. Ceci peut tre dtermin par des expriences utilisant des combinaisons
logiques de divers attributs. Ces derniers sont slectionns par lutilisateur pour tre utiliss
dans le groupement ou la classification des donnes (Clustering or classification).
Il existe plusieurs mthodes de classification ; lune des plus simples est de grouper les
donnes de faon squentielle. Linconvnient de cette procdure, bien quelle peut tre
rapide, est que les clusters (grappes) sont arbitraires et quaucune relation logique nest
maintenue entre eux.
Une autre approche, plus classique encore, est le dveloppement statistique des
donnes, en formant la matrice de covariance et en calculant les composantes principales
(Analyse par Composantes Principales). Ceci donne des valeurs et des vecteurs propres du
plus important au moins important. Le vecteur propre associ la grande valeur propre,
reprsente la combinaison la plus abondante dans lensemble des donnes. Celui suivant
reprsente la combinaison la moins abondante et ainsi de suite. Par la suite, on peut classer les
donnes en calculant la distance entre ces dernires et les vecteurs propres, et retenir ceux qui
ont la distance euclidienne minimale. Par cette mthode on est contraint de rsoudre un
systme matriciel assez grand, compos de toutes les corrlations des donnes.
La caractristique de maintenir les liaisons entre les diffrents clusters, reprsente
laspect important de la mthode de Kohonen : Kohonens Self-Organizing Map(KOHONEN,
T., 1988), o ces clusters sont forms dune faon tels quils sapparentent dans une manire
organise. Cette mthode prsente une certaine analogie avec le cerveau humain quant la
faon dorganiser les informations.
Lalgorithme SOM peut tre utilis dans linterprtation sismique, suivant deux
manires principales : une premire pour gnrer les clusters et une autre pour la
classification.
Chapitre : 05 Attributs sismiques 3D
- 91 -
La mthode Kohonens Self-Organizing Map peut tre utilise aussi dans
linterprtation lithologique, par plusieurs manires diffrentes. La plus vidente est dutiliser
une portion type dune section 2D ou 3D, sur laquelle le rseau apprend (mmorise) les
caractristiques de cette portion : cest les poids de rfrence (refrence weights), qui seront
comme classificateurs pour le reste des donnes.
Si une colonne lithologique (stratigraphique) est disponible, on peut alors utiliser diffrentes
combinaisons logiques des attributs, pour dterminer celle qui ressemble le plus la colonne
stratigraphique.
CONCLUSION :
Lanalyse de la trace sismique comme une partie relle dun signal complexe a permis la
dtermination aise des attributs sismiques complexes, telles que : lamplitude instantane ou
lenveloppe de la trace sismique, la phase instantane, la frquence instantane
La reprsentation en couleur des attributs sismiques aide dans linterprtation des
donnes relatives la stratigraphie. Lenveloppe dcrit la rflectivit, donc elle procure des
informations sur le contraste dimpdance acoustique. La phase instantane accentue la
cohrence latrale. La frquence instantane est utile pour la corrlation, et apparat parfois
comme indicateur daccumulation dhydrocarbures. La polarit apparente est utilise, parfois,
pour lidentification des accumulations de gaz.
Les variations latrales observes sur toutes les reprsentations peuvent tre utiles pour
la localisation des changements stratigraphiques.