Filtrage de l’impédance-sismique à l’aide de … · Comparaison des puits aux traces utilisées lors de ... estimée aux points données (sismique et puits). Journées de Géostatistique,

Journées de Géostatistique, Fontainebleau 2011

D. ABURTO

Filtrage de l’impédance-sismiqueà l’aide de l’impédance-log

Plan� Introduction

� Données� Impédance-sismique, impédance-log� Horizontalisation, trace sismique

2

� Analyse exploratoire

� Estimation� Modèle de variogramme� Résultats

� Conclusions et perspectives

Introduction

� ANDRA (Agence Nationale pour la gestion des Déchets Radioactifs).

� Faisabilité du stockage en

3

� Faisabilité du stockage en profondeur dans une Fm géologique.

� Besoin de modèles de propriétés géomécaniques.

Introduction� 230 km à l’Est de

Paris.

� Formation argileuse du Callovo-oxfordien (Jr).

� Structure monoclinale. � Structure monoclinale. Épaisseur : 130-160 m.

� Encadrée par des formations calcaires.

� Propriétés géomécaniques liées au contenu minéralogique.

4

Données

Impédance-log Impédance-sismique

� Issues de mesures diagraphiques le long des puits.

� Issues de sismique réflexion 2D.

5

� Vp x ρ

� Support 1 ms.

� Même fréquence que la sismique : Imp-log-filt.

� Inversion géophysique à partir de l’imp-log.

� Support 1 ms.

� Fréquence ~70 Hz (résolution de ~12 m).

Données

� Impédance-sismique 10 profils.(chaque profil=millier CDP)

6

� Impédance-log-filtrée4 puits.

� Objectif :Trouver un lien entre l’impédance-sismique et l’impédance-log-filtrée.

Approche� Profil considéré comme un plan vertical.

� Horizontalisation des impédances par rapport à la surface S1.

7

� Comparaison des puits aux traces utilisées lors de l’inversion géophysique .

Analyse exploratoire

8

Variogrammes verticaux

� 4 Puits.

� 4 traces sismiques.

Variogrammes individuelsVariogrammes individuels

Comportement similaire

9

Variogrammes verticaux simples et croisé

44.3 10 (11 )I f Zl Cub msγγ ×= +

3 433.5 10 (12.5 ) 22 10 (42 ),

313 10 cos(8.8 , 40 )

Cub ms Cub msIlf Is

Exp ms ms

γ × += × +

+ ×

4 42 10 cos(10.8 ,90 ) 4 10 (7 )Is Z Exp ms ms Cub msγγ = + × + × +

10

s sI RZ= +

lf lfZI R= +

èHypoth se

2 2 lf snon corrélation à entre Z et R et R

ZIs Rsγ γ γ= +

,Ilf Is Zγ γ=ZIlf Rlfγ γγ= +

416 10 (14.5 )

Is Z

Sph ms+ ×

Validation de l’hypothèse

- -lfs fs lI RIY R= = Rs lfY Rγγ γ= +

11

ObservationsVariogrammes verticaux par tranches le long des profils

� Variogrammes verticaux différents d’un profil à l’autre.

� Certains différents aussi dans le même profil.

12

� Géologie (liée à Z)homogénéité latérale.

� Donc, différences associées aux résidus.

Modélisation du variogramme vertical

,Ilf Is Zγ γ=

,3 433.5 10 (12.5 ) 22 10 (42 )

313 10 cos(8.8 , 40 )

Ilf Is Cub ms Cub ms

Exp ms ms

γ × + × +

+ ×

=

13

NNNN

À partir de chaque profil sismique

S SI RZ= +

ZIs Rsγ γ γ= +

NNNN

Modélisation du variogramme horizontalOn ne dispose que du variogramme de l’imp-sism

14

Ajustement des paramètres à partir des structures utilisées pour la direction verticale.

Estimation de la composante commune� À partir de :

� Estimateur :� Sans Biais� Variance minimale

ZIs Rsγ γ γ= +

*Z Isα αλ= ∑

S SI RZ= +

15

� Système de krigeage :

� Variance d’estimation :

( ) ( ) ( )ZIs αγ λ γ× =

[ * - ] [ , ] ( )SVar Z Z Z Z Var Rα αλ γ µ= + +∑

Voisinage

Voisinageuniquement vertical

VoisinageBidimensionnel

16

Données sismiques

Données brutes v/s données filtrées

Résidu sismique

17

Composante commune estimée (Z)

Conclusions et perspectives� L’estimation de la composante commune (Z), correspond à un filtrage

de la sismique.

� On utilise les puits juste pour connaître la structure de Z, au travers du variogramme croisé entre l’imp-log-filtrée et l’imp-sismique.

� La structure du résidu sismique est déterminée profil par profil ou par CDP.

18

CDP.

� L’estimation avec voisinage bi-dimensionnel :� Très longue (18 h)� Basée sur une interprétation de la structure latérale de Z.� Peu de différences avec voisinage uniquement vertical (CDP par CDP)

� Étude du lien entre les paramètres géomécaniques et la composante Zestimée aux points données (sismique et puits).

Journées de Géostatistique, Fontainebleau 2011

D. ABURTO

Filtrage de l’impédance-sismiqueà l’aide de l’impédance-log