FR2858063A1 - Procede d'estimation du taux de couverture d'illumination dans le domaine migre - Google Patents
Procede d'estimation du taux de couverture d'illumination dans le domaine migre Download PDFInfo
- Publication number
- FR2858063A1 FR2858063A1 FR0308862A FR0308862A FR2858063A1 FR 2858063 A1 FR2858063 A1 FR 2858063A1 FR 0308862 A FR0308862 A FR 0308862A FR 0308862 A FR0308862 A FR 0308862A FR 2858063 A1 FR2858063 A1 FR 2858063A1
- Authority
- FR
- France
- Prior art keywords
- seismic
- illumination
- function
- source
- rate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000005286 illumination Methods 0.000 title claims abstract description 81
- 238000013508 migration Methods 0.000 title claims abstract description 56
- 230000005012 migration Effects 0.000 title claims abstract description 55
- 238000000034 method Methods 0.000 title claims abstract description 52
- 239000013598 vector Substances 0.000 claims abstract description 22
- 238000012937 correction Methods 0.000 claims abstract description 8
- 230000000694 effects Effects 0.000 claims abstract description 7
- 238000011161 development Methods 0.000 claims description 13
- 108050005509 3D domains Proteins 0.000 claims description 9
- 230000014509 gene expression Effects 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 5
- 230000018109 developmental process Effects 0.000 description 12
- 238000002310 reflectometry Methods 0.000 description 10
- 238000013459 approach Methods 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 5
- 230000005526 G1 to G0 transition Effects 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 238000013507 mapping Methods 0.000 description 2
- 230000036962 time dependent Effects 0.000 description 2
- 101150106956 Oxtr gene Proteins 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Electroluminescent Light Sources (AREA)
Abstract
L'invention concerne une estimation du taux I(x,p) de couverture d'illumination sismique dans le domaine 3D migré en un point image x, pour un vecteur p de pendage,caractérisé en ce que, pour chaque paire (source s, récepteur r) du relevé sismique, on évalue le taux I(x,p;s,r) de couverture d'illumination relatif à ladite paire (source s, récepteur r) par mise en oeuvre des étapes de :• détermination du temps de trajet de réflexion tr(xr (p) ; s , r) depuis la source s jusqu'au point de réflexion spéculaire xr sur le réflecteur plan passant par le point image x et perpendiculaire au vecteur p de pendage, puis de retour vers le réflecteur r, à partir du temps de trajet de diffraction td(x;s,r) depuis la source s jusqu'audit point image x puis de retour jusqu'au réflecteur r ;• d'incrémentation dudit taux I(x,p;s,r) de couverture d'illumination relatif à ladite paire (source s, récepteur r) en fonction de la différence entre le temps de trajet de diffraction td(x;s,r ) et le temps de trajet de réflexion tr(xr(p);s,r).
Description
<Desc/Clms Page number 1>
Le domaine d'application de l'invention est celui de la prospection sismique.
L'invention concerne l'imagerie sismique, et plus particulièrement l'estimation du taux de couverture d'illumination dans le domaine 3D migré nécessaire pour l'obtention d'une image haute résolution, à amplitude réelle, de la structure du sous-sol.
La prospection sismique consiste d'une manière générale, à émettre dans le sous-sol, à l'aide d'une ou plusieurs sources sismiques, des ondes sismiques, à enregistrer en surface, en fonction du temps, des données sismiques correspondant aux ondes sismiques réfléchies sur les interfaces géologiques du sous-sol (encore appelées réflecteurs) à l'aide de récepteurs (encore appelés géophones ou hydrophones selon que l'on prospecte à terre ou en mer) puis à traiter ces données pour en extraire des informations utiles quant à la géologie du sous-sol.
Le traitement des données (encore appelées traces sismiques) enregistrées par les récepteurs permet en particulier de concevoir des images présentant des informations relatives à la géologie du sous-sol.
L'imagerie sismique à haute résolution, à amplitude réelle, a pour objectif de fournir une image du sous-sol qui est non perturbée notamment par les effets de la propagation du champ d'onde dans le sous-sol et par les effets de l'acquisition (on parle alors d'empreintes d'acquisition).
La représentation en image de la réflectivité acoustique du sous-sol est réalisée par un procédé dit de migration qui permet de restituer les forme correctes des interfacesgéologiques. Les données sismiques sont pour cela converties depuis le domaine dans lequel elles sont acquises vers le domaine migré dans lequel les interfaces géologiques sont représentées à leurs véritables positions.
Dans le domaine migré, l'échelle verticale peut être le temps (on parle alors de migration temps ou time migration selon la terminologie anglo-saxonne) ou la profondeur (on parle alors de migration profondeur ou depth migration selon la terminologie anglo-saxonne).
<Desc/Clms Page number 2>
La migration est généralement réalisée avant sommation (on parle alors de migration avant somme ou pre-stack migration selon la terminologie anglo- saxonne), avec un mappage direct du signal sismique depuis les traces sismiques vers le domaine migré.
L'identification précise de la géologie du sous-sol en un point particulier est directement reliée à la capacité de l'émission sismique à illuminer ce point. Si un point du sous-sol est hautement illuminé, le relevé sismique permet d'évaluer en toute vraisemblance la géologie du sous-sol en ce point. En revanche, si un point du sous-sol est peu illuminé, il est possible que le relevé sismique fournisse une évaluation imprécise, voire incorrecte, de la géologie du sous-sol en ce point.
Afin d'obtenir une représentation précise du sous-sol, les amplitudes réelles de la réflexion de l'onde sismique sur un réflecteur doivent être parfaitement représentées.
Il est donc nécessaire, de connaître avec exactitude si un changement d'amplitude en un point particulier du sous-sol est le résultat d'un changement dans la géologie du sous-sol en ce point ou le résultat d'une illumination particulière en ce point.
Mais l'illumination des réflecteurs est généralement non uniforme. Elle dépend en particulier : - de la configuration de l'acquisition des données sismiques ; - des variations latérales et verticales de la vitesse de propagation des ondes sismiques dans le sous-sol. "
Afin d'obtenir une image haute résolution de la réflectivité du sous-sol, il est nécessaire, lors de la migration des réflexions sismiques, d'effectuer une correction d'amplitude pour prendre en compte la non-uniformité d'illumination des réflecteurs.
Afin d'obtenir une image haute résolution de la réflectivité du sous-sol, il est nécessaire, lors de la migration des réflexions sismiques, d'effectuer une correction d'amplitude pour prendre en compte la non-uniformité d'illumination des réflecteurs.
Comme cela a déjà été exposé précédemment, cette correction (réalisée à l'aide d'un facteur de pondération) est nécessaire car : - pour imager précisément la réflectivité du sous-sol, il faut corriger les distorsions en amplitude liées à la propagation du champ d'onde dans le sous- sol (divergence sphérique, absorption, directivité, etc.) ;
<Desc/Clms Page number 3>
- pour améliorer la résolution de l'image de réflectivité, il faut prendre en compte la variabilité du taux de couverture d'illumination de migration pour chacune des inclinaisons d'un réflecteur qui sont prises en compte lors de la migration.
Le taux de couverture d'illumination correspond, pour un point image donné du domaine migré, au nombre de paires (source, récepteur) dont les contributions pour la constitution de l'amplitude de l'image de réflectivité interfèrent constructivement (condition de phase stationnaire).
Il a été proposé dans le document : Bleistein, N. 1987, On the Imaging of Reflectors in the Earth : Geophysics,
Soc. Of Expl. Geophysics, 931-942 , une méthode permettant d'obtenir, pour une configuration régulière d'acquisition, un facteur de pondération d'amplitude permettant de réaliser une représentation haute résolution, en amplitude réelle de la réflectivité du sous- sol.
Soc. Of Expl. Geophysics, 931-942 , une méthode permettant d'obtenir, pour une configuration régulière d'acquisition, un facteur de pondération d'amplitude permettant de réaliser une représentation haute résolution, en amplitude réelle de la réflectivité du sous- sol.
Cette méthode est une méthode inverse constituée d'une étape de modélisation (modélisation de traces sismiques à partir d'un modèle de réflectivité), suivie d'une étape de migration (depuis les traces sismiques modélisées, vers une réflectivité imagée), le facteur de pondération étant alors choisi de manière à obtenir une concordance parfaite entre la réflectivité imagée et le modèle de réflectivité.
En suivant cette méthode inverse proposée par Bleistein, et comme cela est représenté sur la figure 1, il apparaît que : - pour une configuration -d'acquisition donnée (répartition des paires (source S, récepteur R)), - pour un point image x donné dans le domaine migré, - pour une inclinaison de réflecteur donnée (indiquée par un vecteur p de pendage), le taux de couverture d'illumination I(#,#) est déterminé par le nombre de paires (source, récepteur) tel que leurs contributions de migration interférent constructivement à l'image du réflecteur de pendage p au point image x .
<Desc/Clms Page number 4>
Les expressions développées par Bleistein pour le taux de couverture d'illumination I(#,#) intègrent en une seule fois l'effet de la totalité de la configuration géométrique de l'acquisition.
Par conséquent, la méthode proposée par Bleistein ne peut être utilisée que pour une distribution régulière et paramétrique de l'ensemble des paires (source, récepteur).
Cette condition restrictive d'application en limite par conséquent l'utilisation à un nombre réduit de configurations d'acquisition.
Pour des configurations d'acquisition à géométries irrégulières, il est donc
nécessaire de définir les taux de couverture d'illumination individuels I ( x , p , s , r), pour chacune des paires (source s, récepteur r ) impliquées dans la migration. Lesdits taux individuels contribuent bien entendu au taux de couverture
d'illumination total 1(x, p ) selon I(x, p) ~ I(x, p;s,#) ).
nécessaire de définir les taux de couverture d'illumination individuels I ( x , p , s , r), pour chacune des paires (source s, récepteur r ) impliquées dans la migration. Lesdits taux individuels contribuent bien entendu au taux de couverture
d'illumination total 1(x, p ) selon I(x, p) ~ I(x, p;s,#) ).
#, #
L'étude conduite par Bleistein est toutefois à la base d'un nombre de développements qui ont eu pour but d'estimer le taux de couverture d'illumination dans le domaine migré, pour n'importe quelle configuration géométrique de l'acquisition, et pour n'importe quelle distribution des paires (source, récepteur).
L'étude conduite par Bleistein est toutefois à la base d'un nombre de développements qui ont eu pour but d'estimer le taux de couverture d'illumination dans le domaine migré, pour n'importe quelle configuration géométrique de l'acquisition, et pour n'importe quelle distribution des paires (source, récepteur).
Au cours des années, il a été proposé diverses approximations pour l'onde sismique, qui est l'élément clé pour l'estimation de la condition de phase stationnaire (condition d'interférence constructive). @
Les premiers développements ont consisté à modéliser par lancer de rayons ( ray tracing selon la terminologie anglo-saxonne), pour chacune des paires (source s , récepteur r ), les points de réflexion spéculaires le long d'un horizon (carte 2D) sélectionné dans le domaine migré (3D).
Les premiers développements ont consisté à modéliser par lancer de rayons ( ray tracing selon la terminologie anglo-saxonne), pour chacune des paires (source s , récepteur r ), les points de réflexion spéculaires le long d'un horizon (carte 2D) sélectionné dans le domaine migré (3D).
Cette étape de modélisation de la réflexion permet de déterminer les quantités suivantes : - La localisation des points de réflexion #r(#,#) le long de l'horizon sélectionné ;
- Le temps de trajet tri s , r ) de réflexion spéculaire ;
- Le temps de trajet tri s , r ) de réflexion spéculaire ;
<Desc/Clms Page number 5>
Le temps td de trajet de diffraction le long de l'horizon sélectionné peut alors être estimé par développement en série de Taylor du temps de trajet tr de réflexion
autour du point xr ( s , r ) de réflexion spéculaire (détermination des dérivées première et seconde \\tr(xr) et àxjr(xr)).
autour du point xr ( s , r ) de réflexion spéculaire (détermination des dérivées première et seconde \\tr(xr) et àxjr(xr)).
La figure 2 illustre cette procédure au cours de laquelle le temps de trajet de diffraction td(x) au point image x est estimé à l'aide du temps de trajet de réflexion tr(#r) et plus particulièrement par son développement en série de Taylor autour du point de réflexion spéculaire #r.
La sommation de tous les mappages d'ondelette obtenus pour chacune des paires (source, récepteur) permet ensuite d'estimer le taux de couverture d'illumination dans le domaine migré, le long de l'horizon sélectionné.
En modélisant l'ondelette sismique selon un pic impulsionnel, à bande de fréquence infinie, il a été montré dans le document : Cain, G., Cambois, G., Gehin, M. et Hall, R., 1998, Reducing risk in seismic acquisition and interprétation of complex targets using a Gocad-based 3-D modeling tool, 68th Ann. Internat. Mtg : Soc. of Expl. Geophys., 2072-2075 , qu'il est possible d'obtenir une carte d'impact, seule la contribution spéculaire sur un horizon sélectionné étant alors effectivement prise en compte.
En prenant en compte une approximation réaliste d'une ondelette sismique à bande limitée, il a été montré dans les documents : - Schneider, W., Jr pnd Winbow, G., 1999, Efficient and accurate modeling of 3D seismic illumination, 69th Ann. Internat, Mtg : Soc. Of Expl. Geophys,
Expanded Abstract, 633-636 ; - Laurain, R. et Vinje, V., 2001, Pre-stack migration and illumination maps,
71th Ann. Internat, Mtg : Soc. Of Expl. Geophys, Expanded Abstract, 929-
932 ; qu'il est possible d'améliorer l'estimation du taux de couverture d'illumination le long d'horizons sélectionnés.
Expanded Abstract, 633-636 ; - Laurain, R. et Vinje, V., 2001, Pre-stack migration and illumination maps,
71th Ann. Internat, Mtg : Soc. Of Expl. Geophys, Expanded Abstract, 929-
932 ; qu'il est possible d'améliorer l'estimation du taux de couverture d'illumination le long d'horizons sélectionnés.
<Desc/Clms Page number 6>
Ces approches, basées en particulier sur un balayage d'horizons sélectionnés, se sont montrées relativement fiables et précises. Leur efficacité est cependant limitée, notamment du fait que : la compensation du taux de couverture d'illumination est limitée à un nombre fini d'horizons sélectionnés et l'approche, au contraire de celle développée par
Bleistein, n'est donc pas totalement 3D ; pour chacune des paires (source s, récepteur #) impliquées dans le migration, une étape de modélisation dédiée doit être réalisée après la migration. Or cette étape de modélisation, par lancer de rayons, pour la
détermination des quantités (xr' V)r,J1x,xtr) mentionnées précédemment, le long de chacun des horizons sélectionnés, est particulièrement coûteuse en temps.
Bleistein, n'est donc pas totalement 3D ; pour chacune des paires (source s, récepteur #) impliquées dans le migration, une étape de modélisation dédiée doit être réalisée après la migration. Or cette étape de modélisation, par lancer de rayons, pour la
détermination des quantités (xr' V)r,J1x,xtr) mentionnées précédemment, le long de chacun des horizons sélectionnés, est particulièrement coûteuse en temps.
Afin de surmonter ces limitations, et de se rapprocher d'une approche volumique 3D, une solution, en deux étapes, pour l'estimation puis la compensation du taux de couverture d'illumination a été proposée, en particulier dans les documents : - Albertin, U., Bloor, R., Beasley, C., Chang, W., Jaramillo, H., Mobley, E. and Yingst, D., 1999, Aspects of true amplitude migration, 69th Ann. Internat,
Mtg : Soc. Of Expl. Geophys, Expanded Abstract, 1247-1250 ; - Audebert, F., Froidevaux, P. Huard, I., Nicoletis, L. and Svay-Lucas, J.,
2000, A mutli angle toolbox for restored amplitude images and AVA gathers
70th Ann. Internat, Mtg: Soc. Of Expl. Geophys, Expanded Abstract, 1020-
1023 .
Mtg : Soc. Of Expl. Geophys, Expanded Abstract, 1247-1250 ; - Audebert, F., Froidevaux, P. Huard, I., Nicoletis, L. and Svay-Lucas, J.,
2000, A mutli angle toolbox for restored amplitude images and AVA gathers
70th Ann. Internat, Mtg: Soc. Of Expl. Geophys, Expanded Abstract, 1020-
1023 .
La première étape consiste, dans le domaine migré, en une analyse du taux de couverture d'illumination multi-pendage utilisant les tables pré-établies pour la migration de Kirchhoff des temps de trajet de diffraction td (x ; s , r ).
Lors de la seconde étape, les inverses des taux de couverture d'illumination estimés sont utilisés afin de réaliser la sommation pondérée des données sismiques.
Pour chaque point image x et pour toute une gamme de pendage p de réflecteur, il a été proposé d'incrémenter un compteur d'occurrence d'illumination
<Desc/Clms Page number 7>
I(#,#), chaque fois qu'une contribution de migration atteint le point image x avec la composante p de pendage.
De manière similaire à ce qui avait été proposé par Cain et al (cf. ci-dessus), le critère utilisé pour augmenter le compteur d'occurrence d'illumination prend en compte un modèle impulsionnel, à bande de fréquence infinie, de l'ondelette sismique.
Cette approximation hautes fréquences est telle que, pour une composante de pendage donnée, seules les contributions de migration spéculaires correspondant à ladite composante de pendage sont prises en compte pour l'estimation de l'illumination.
Et cette approximation hautes fréquences présente un certain nombre de limitations : - l'estimation de l'illumination est biaisée puisque la nature à bande limitée du signal sismique n'est pas prise en compte ; - Afin d'éviter la division par des zéros, les taux de couverture d'illumination sont rangés par classes de pendage. La taille des classes est toutefois un élément critique. Si la taille de la classe est trop petite, les risques d'instabilité numérique sont plus grands ; si la taille de la classe est trop grande, la précision quant à l'estimation de la couverture d'illumination est diminuée.
L'invention vise un procédé d'estimation du taux de couverture d'illumination sismique dans le domaine 3D migré non sujet à ces limitations.
Plus précisément/l'invention a pour objectif de proposer un procédé qui prend en compte la nature à bande limitée du signal sismique afin d'estimer des taux de couverture d'illumination plus réalistes, et qui, dans le même temps, résout le problème de la taille des classes.
A cet effet, l'invention propose, selon un premier aspect, un procédé d'estimation du taux I(#,#) de couverture d'illumination sismique dans le domaine 3D migré en un point image x, pour un pendage de vecteur p , caractérisé en ce que, pour chaque paire (source s , récepteur r) du relevé sismique, il évalue le
<Desc/Clms Page number 8>
taux I(x,p;s,#) de couverture d'illumination relatif à ladite paire (source s, récepteur r ) par mise en #uvre des étapes de : #détermination du temps de trajet de réflexion t( xr ( p ) ; s, r ) depuis la source s , jusqu'au point de réflexion spéculaire xr sur le réflecteur plan passant par le point image x et perpendiculaire au vecteur p de pendage, puis de retour vers le réflecteur r , à partir du temps de trajet de diffraction td(#; s, r) depuis la source s jusqu'audit point image x, puis de retour jusqu'au réflecteur r ;
#d'incrémentation dudit taux I(z,p;s,#) de couverture d'illumination relatif à ladite paire (source s , récepteur r ) en fonction de la différence entre le temps de trajet de diffraction td( x ; s, r ) et le temps de trajet de réflexion tr( Xr ( p ) ; s' , r).
Des aspects préférés, mais non limitatifs du procédé selon le premier aspect de l'invention sont les suivants : - le procédé comporte en outre une étape de sommation de chacun des taux I(x,p;s,r) de couverture d'illumination relatif à une paire (source s, récepteur r), de manière à déterminer le taux de couverture d'illumination total
j(xP) ~ (xPS#) , i,f - au cours de l'étape** d'incrémentation, le taux 1(x, p; s, ï) de couverture d'illumination peut être incrémenté en utilisant une fonction d'incrémentation i(td, tr; s,F), selon (x, p;S,r) = I(xpsr) + i(td, tr; s,r), - la fonction d'incrémentation i peut être fonction de l'ondelette sismique s(t), - la fonction d'incrémentation i peut s'exprimer en fonction de la dérivée selon le temps de l'ondelette sismique, - il est envisagé de tenir compte d'une correction a priori w(#;#,#) du taux de couverture d'illumination, correction qui est généralement mise en jeu lors de la migration,
j(xP) ~ (xPS#) , i,f - au cours de l'étape** d'incrémentation, le taux 1(x, p; s, ï) de couverture d'illumination peut être incrémenté en utilisant une fonction d'incrémentation i(td, tr; s,F), selon (x, p;S,r) = I(xpsr) + i(td, tr; s,r), - la fonction d'incrémentation i peut être fonction de l'ondelette sismique s(t), - la fonction d'incrémentation i peut s'exprimer en fonction de la dérivée selon le temps de l'ondelette sismique, - il est envisagé de tenir compte d'une correction a priori w(#;#,#) du taux de couverture d'illumination, correction qui est généralement mise en jeu lors de la migration,
<Desc/Clms Page number 9>
- l'étape de détermination comprend le développement en série de Taylor, du
second ordre, du temps de trajet de diffraction td(x ; s , r ) autour du point image x , - le point de réflexion spéculaire #,(#) est déterminé le long dudit réflecteur de sorte que le temps de trajet de diffraction audit point réflexion spéculaire xr (p) est stationnaire, - l'étape de détermination utilise les cartes des isochrones de migration td(#; s , r) spécifiées pour toute paire (source s, récepteur r) impliquée dans la migration et tout point image x du domaine 3D migré, - l'estimation du taux I(#,#) de couverture d'illumination sismique dans le domaine 3D migré est réalisé lors de la migration par sommation de Kirchhoff des données sismiques enregistrées au cours du relevé sismique 3D.
second ordre, du temps de trajet de diffraction td(x ; s , r ) autour du point image x , - le point de réflexion spéculaire #,(#) est déterminé le long dudit réflecteur de sorte que le temps de trajet de diffraction audit point réflexion spéculaire xr (p) est stationnaire, - l'étape de détermination utilise les cartes des isochrones de migration td(#; s , r) spécifiées pour toute paire (source s, récepteur r) impliquée dans la migration et tout point image x du domaine 3D migré, - l'estimation du taux I(#,#) de couverture d'illumination sismique dans le domaine 3D migré est réalisé lors de la migration par sommation de Kirchhoff des données sismiques enregistrées au cours du relevé sismique 3D.
Selon un autre aspect, l'invention propose un procédé de correction des amplitudes des données sismiques enregistrées au cours d'une prospection sismique 3D afin de compenser l'effet d'une illumination non uniforme des réflecteurs du sous-sol, caractérisé en ce qu'il comporte les étapes consistant à : # estimer le taux de couverture d'illumination I(#,#) par le procédé selon le premier aspect de l'invention,
#utiliser l'inverse I'(x,p ) dudit taux comme facteur de pondération à appliquer à chacune desdites amplitudes des données sismiques.
#utiliser l'inverse I'(x,p ) dudit taux comme facteur de pondération à appliquer à chacune desdites amplitudes des données sismiques.
Selon encore un utre aspect, l'invention propose un procédé pour la sélection, en fonction de la cible d'une prospection sismique 3D, d'une géométrie d'acquisition parmi une pluralité de géométries d'acquisition, caractérisé en ce qu'il comporte les étapes consistant à : # déterminer, pour chacune des géométries d'acquisition considérées, le taux de couverture d'illumination I(#,#) par le procédé selon le premier aspect de l'invention, # sélectionner la géométrie d'acquisition procurant le taux de couverture d'illumination optimum en fonction de la cible.
<Desc/Clms Page number 10>
D'autres aspects, buts et avantages de l'invention apparaîtront mieux à la lecture de la description détaillée suivante d'une forme de réalisation préférée de celle-ci, donnée à titre d'exemple non limitatif et faite en référence aux dessins annexés sur lesquels : - la figure 1, déjà commentée, illustre la contribution constructive ou non, en fonction de s(td - tr), de la trace sismique migrée à l'image du réflecteur penté en x ; - la figure 2, déjà commentée, illustre l'estimation réalisée dans le cadre de l'état de la technique du temps de trajet de diffraction au point image td(x) à partir d'une développement en série de Taylor du temps de trajet de réflexion au point de réflexion spéculaire tr(#r) ; - la figure 3 illustre le procédé selon l'invention selon lequel on estime le point
de réflexion spéculaire x, ainsi que le temps de trajet de réflexion tr( xr) à partir d'un développement en série de Taylor du temps de trajet de diffraction au point image td(x).
de réflexion spéculaire x, ainsi que le temps de trajet de réflexion tr( xr) à partir d'un développement en série de Taylor du temps de trajet de diffraction au point image td(x).
- La figure 4 est un organigramme présentant de manière schématique les principales étapes du procédé selon l'invention.
L'invention propose un procédé d'estimation du taux de couverture d'illumination dans le domaine migré I(x,p;s,r) pour n'importe quelle position migrée x , n'importe quelle composante de pendage p et n'importe quelle distribution des traces sismiques, en utilisant uniquement les quantités les plus basiques mises en jeu dans la migration par sommation de Kirchhoff, à savoir :
#les coordonnées ( s , r ) des paires (source, récepteur), s représentant la position d'une source S et r celle d'un récepteur R ; # les cartes des isochrones de migration td (x ; #, r ) pour chacune des paires (source s , récepteur r ) impliquées dans la migration.
#les coordonnées ( s , r ) des paires (source, récepteur), s représentant la position d'une source S et r celle d'un récepteur R ; # les cartes des isochrones de migration td (x ; #, r ) pour chacune des paires (source s , récepteur r ) impliquées dans la migration.
L'invention s'applique sans distinction à l'imagerie en temps et à l'imagerie en profondeur.
<Desc/Clms Page number 11>
L'invention propose un procédé du type comprenant une étape de modélisation suivie d'une étape de migration.
En ce qui concerne l'étape de modélisation du point de réflexion xr et du temps de trajet de réflexion tr(#r), le procédé objet de l'invention applique la condition de phase stationnaire à un développement en série de Taylor des temps de trajet de diffraction td(# ).
L'ordre dudit développement est un développement du second ordre.
En ce qui concerne l'étape de migration, le procédé objet de l'invention utilise les cartes des isochrones de migration disponibles pour la migration classiquement réalisée par sommation de Kirchhoff.
Pour chacune des paires (source s, récepteur r), le taux de couverture
d'illumination I (z, p ; S , Y ) est mis à jour, en référence à la figure 3, selon l'équation suivante : \(x,p;s,r) = \(x,p;s,r) +i(td,tr; s,r) Equation (1) Avec #x, le point image dans le domaine migré x = (X1, x2)x3) T ; # p , le vecteur de pendage représentatif de la composante de pendage sélectionnée, p = (p1, P2,P3) T ; # s , les coordonnées de la source,s = (si, S2, S3) T ;
#Y, les coordonnées du" récepteur, r = (r-i, r2, r3) T ; # td ( x ; s , r), le temps de trajet de diffraction depuis la source s, jusqu'à la position migrée #, et de retour jusqu'au récepteur r (trajet indiqué par td sur la figure 3 et représenté en trait plein) ;
#tr( Xr (p ); s , r ), le temps de trajet de réflexion depuis la source s, jusqu'au point xr de réflexion spéculaire sur le réflecteur plan passant par le point image x et perpendiculaire au vecteur p de pendage, et de retour jusqu'au récepteur r (trajet indiqué par tr sur la figure 3 et représenté en trait pointillé)
d'illumination I (z, p ; S , Y ) est mis à jour, en référence à la figure 3, selon l'équation suivante : \(x,p;s,r) = \(x,p;s,r) +i(td,tr; s,r) Equation (1) Avec #x, le point image dans le domaine migré x = (X1, x2)x3) T ; # p , le vecteur de pendage représentatif de la composante de pendage sélectionnée, p = (p1, P2,P3) T ; # s , les coordonnées de la source,s = (si, S2, S3) T ;
#Y, les coordonnées du" récepteur, r = (r-i, r2, r3) T ; # td ( x ; s , r), le temps de trajet de diffraction depuis la source s, jusqu'à la position migrée #, et de retour jusqu'au récepteur r (trajet indiqué par td sur la figure 3 et représenté en trait plein) ;
#tr( Xr (p ); s , r ), le temps de trajet de réflexion depuis la source s, jusqu'au point xr de réflexion spéculaire sur le réflecteur plan passant par le point image x et perpendiculaire au vecteur p de pendage, et de retour jusqu'au récepteur r (trajet indiqué par tr sur la figure 3 et représenté en trait pointillé)
<Desc/Clms Page number 12>
# i (td, tr ; #,#) est une fonction d'incrémentation dont l'expression sera détaillée ci-après.
Comme cela a déjà été mentionné précédemment, le procédé selon l'invention utilise avantageusement les cartes des isochrones de migration td (x ; s, r) mises en jeu dans la migration par sommation de Kirchhoff, pour chacune des paires (source s, récepteur r) impliquées dans la migration. Ces cartes sont ainsi des données d'entrée du procédé selon l'invention.
Il utilise également les dérivées première et seconde oxcd(x;s,#), t1x,xtd(X;S,r) des isochrones de migration lorsque celles-ci sont disponibles. Dans le cas contraire, lesdites dérivées ne sont pas des données d'entrée, mais peuvent être calculées dans le cadre du procédé selon l'invention
La figure 4 est un organigramme représentant les différentes étapes du procédé selon un mode de réalisation préféré de l'invention.
La figure 4 est un organigramme représentant les différentes étapes du procédé selon un mode de réalisation préféré de l'invention.
Ledit procédé consiste à réaliser, pour chacune des paires (source s, récepteur r ) impliquées dans la migration, une étape 1 mise en #uvre pour tout point image x du domaine 3D migré.
Ladite étape 1 comprend, pour toute composante de pendage p sélectionnée pour le calcul de l'illumination I(x,p;s,r), deux opérations 2, 3.
La première opération 2 consiste à déterminer le temps de trajet de
réflexion tr( x, ( p ) ; s , r) ou point de réflexion spéculaire x, (p sur le réflecteur plan passant par le point image # et perpendiculaire au vecteur p de pendage.
réflexion tr( x, ( p ) ; s , r) ou point de réflexion spéculaire x, (p sur le réflecteur plan passant par le point image # et perpendiculaire au vecteur p de pendage.
La détermination de l'opération 2 est réalisée à partir du temps de trajet de diffraction td (#; s , #) connu des cartes des isochrones de migration.
La seconde opération 3 consiste à incrémenter le taux de couverture
d'illumination dans le domaine migré 1(x, p;s,#) d'une fonction d'incrémentation i(td, tr ; s,r), ladite fonction d'incrémentation prenant en compte: - la différence ta ( x ; , ;) - tr(;r (p ); s , Y ) }.entre les temps de trajet mentionnés ci-dessus ;
d'illumination dans le domaine migré 1(x, p;s,#) d'une fonction d'incrémentation i(td, tr ; s,r), ladite fonction d'incrémentation prenant en compte: - la différence ta ( x ; , ;) - tr(;r (p ); s , Y ) }.entre les temps de trajet mentionnés ci-dessus ;
<Desc/Clms Page number 13>
l'ondelette sismique s(t), et, éventuellement, le terme d'amplitude de la migration w(#;#,#).
A l'issue des opérations 2 et 3, le taux d'amplitude d'illumination I(#,#;#,#) est connu pour un point image x donné et une composante p de pendage.
En répétant les opérations mentionnées ci-dessus pour chacune des paires (source s , récepteur r), on détermine finalement le taux d'illumination dans le
domaine migré I ( x, p ).
domaine migré I ( x, p ).
On additionne pour cela l'illumination relative à chacune des paires (source s , récepteur r ) selon 7(;c,p) = /(,;,F).
Le point de réflexion #r, ainsi que le temps de trajet de réflexion tr sont déterminés (au cours de l'opération 2 mentionnée ci-dessus) à partir d'une analyse en phase stationnaire, par un développement en série de Taylor, du second ordre, du temps de trajet de diffraction td autour du point image #.
A partir de cette analyse (cf. équation (2) ci-dessous), on obtient les expressions suivantes pour #r et tr (équations (3) et (4)) :
td(x,s,r)-tr(xr,s,r)+(Oxtr)T (x xr}-i-2('x xr)T (x,xtrJlx JCr)-i-.... Equation (2) xr(p)= 5c-M.F'\b Equation (3) t( xr ( P ) ; s , # ) = ta( x ; s , r ) - 1 . b T . F-' . b Equation (4)
Où dans les équations (3) et (4),
M est une matrice (3 x 2) décrite par deux vecteurs s'étendant le long du plan de réflexion. Ces deux vecteurs sont par conséquent perpendiculaires au vecteur de pendage p .
td(x,s,r)-tr(xr,s,r)+(Oxtr)T (x xr}-i-2('x xr)T (x,xtrJlx JCr)-i-.... Equation (2) xr(p)= 5c-M.F'\b Equation (3) t( xr ( P ) ; s , # ) = ta( x ; s , r ) - 1 . b T . F-' . b Equation (4)
Où dans les équations (3) et (4),
M est une matrice (3 x 2) décrite par deux vecteurs s'étendant le long du plan de réflexion. Ces deux vecteurs sont par conséquent perpendiculaires au vecteur de pendage p .
# est un vecteur (2x1) des dérivées du second ordre du temps de trajet de diffraction le long du plan de réflexion,
i =Mr.(V) Equation (5)
i =Mr.(V) Equation (5)
<Desc/Clms Page number 14>
# F est une matrice (2 x 2) des dérivés du premier ordre du temps de trajet de diffraction le long du plan de réflexion,
F = MT .(A,,,td)'M Equation (6)
On décrit ci-après l'obtention des équations (3) et (4).
F = MT .(A,,,td)'M Equation (6)
On décrit ci-après l'obtention des équations (3) et (4).
On précise ici que le vecteur p de pendage représente le pendage apparent d'un réflecteur dans le domaine migré, ce réflecteur passant par le point image #.
Le vecteur p de pendage s'exprime ainsi selon :
p (8x3/8x ôx3 ôx2 ,1) Equation (7a) bzz représente le pendage apparent le long de l'axe Xi et aX31ôx, , le pendage apparent selon l'axe x2.
p (8x3/8x ôx3 ôx2 ,1) Equation (7a) bzz représente le pendage apparent le long de l'axe Xi et aX31ôx, , le pendage apparent selon l'axe x2.
0 et # représentent respectivement l'inclinaison et l'azimut du réflecteur penté passant par le point image x où l'illumination est évaluée.
On a ainsi les expressions suivantes:
cos8= 1+ Pi + P2ti/2 Equation (7b) tgV = P2 /Pi Equation (7c)
A partir du développement en série de Taylor du temps de trajet de diffraction (connu de la carte des isochrones de migration), on a, comme on l'a vu précédemment :
Íct(xr) Íd(x) + \Oxtd )T \'r 'x ) + 2 \ xr 'x )T x,xtd' {'xr x J+... Equation (2) Le point de réflexion , est tout d'abord estimé le long du plan de réflexion pour lequel le temps de trajet de diffraction est stationnaire (condition pour une réflexion spéculaire).
cos8= 1+ Pi + P2ti/2 Equation (7b) tgV = P2 /Pi Equation (7c)
A partir du développement en série de Taylor du temps de trajet de diffraction (connu de la carte des isochrones de migration), on a, comme on l'a vu précédemment :
Íct(xr) Íd(x) + \Oxtd )T \'r 'x ) + 2 \ xr 'x )T x,xtd' {'xr x J+... Equation (2) Le point de réflexion , est tout d'abord estimé le long du plan de réflexion pour lequel le temps de trajet de diffraction est stationnaire (condition pour une réflexion spéculaire).
Afin de ne considérer que les positions appartenant au plan penté, on introduit un nouveau vecteur xp ,tel que :
, 0 x, = x +M.zp = z + 0 1 Xp2 Equation (7d) - ~p 2., 1)
, 0 x, = x +M.zp = z + 0 1 Xp2 Equation (7d) - ~p 2., 1)
<Desc/Clms Page number 15>
On note ici que, comme désiré, les deux vecteurs de la matrice M sont bien perpendiculaires au vecteur p de pendage.
La réécriture de l'équation (2) en utilisant l'équation (7d) fournit l'expression suivante :
td('p)=td(y)+(-v,,t,YM.. p + 2 i..T.MT.A.,^t p d M.JET p Equation (7e) soit : td xp =td x + .xp 2 +... Equation (7f) Dans cette équation (7f) : bT = (xtd )T .M , d'où MT.(,t@)@ et MT 'A,,,,td.m
On en déduit par conséquent la position xrp pour laquelle le temps de trajet de diffraction td(xp) est stationnaire :
Vxptd(xp)1 x", ~ , x,p =-F-1]; => xr =x+M.xrp Equation (7g) En combinant les équations (7g) et (2), on obtient l'expression du temps de trajet de réflexion de l'équation (4) : tr = td (xr) td (x) # bT .F~\b .
td('p)=td(y)+(-v,,t,YM.. p + 2 i..T.MT.A.,^t p d M.JET p Equation (7e) soit : td xp =td x + .xp 2 +... Equation (7f) Dans cette équation (7f) : bT = (xtd )T .M , d'où MT.(,t@)@ et MT 'A,,,,td.m
On en déduit par conséquent la position xrp pour laquelle le temps de trajet de diffraction td(xp) est stationnaire :
Vxptd(xp)1 x", ~ , x,p =-F-1]; => xr =x+M.xrp Equation (7g) En combinant les équations (7g) et (2), on obtient l'expression du temps de trajet de réflexion de l'équation (4) : tr = td (xr) td (x) # bT .F~\b .
Comme cela a été exposé précédemment, la taux de couverture d'illumination est incrémenté, au cours de l'opération 3, en fonction de la différence entre le temps td de trajet de diffraction et temps tr de trajet de réflexion.
Même si le temps de trajet de réflexion tr(3cr(p ) ; s, r) diffère du temps de trajet de diffraction td(x ; #, #) du fait de la nature à bande limitée des données sismiques, la migration de la réflexion ayant lieu à # peut contribuer constructivement à l'image du réflecteur pente en x .
La sélection entre contributions de migration illuminantes (interférences constructives) et contributions de migration non illuminantes (interférences destructives) est réalisée grâce la fonction d'incrémentation i(td, tr ; #,#).
Selon un mode de réalisation de l'invention, la fonction d'incrémentation
s'exprime selon i(td,tr; s,r ) = s(td(x;s,#)-t,.(xr(p);s,#)), où :
s'exprime selon i(td,tr; s,r ) = s(td(x;s,#)-t,.(xr(p);s,#)), où :
<Desc/Clms Page number 16>
s(t) est l'ondelette sismique, qui peut être monochromatique, à bande limitée, ou à bande infinie et qui peut être choisie dépendante du temps et de l'espace,
set) est la dérivée selon le temps de l'ondelette sismique s(t) = d[s(t)]/ dt ,
Une dérivée temporelle est généralement appliquée à chacune des traces sismiques avant la migration de Kirchhoff (en temps ou en profondeur). Ainsi, notamment dans un souci d'uniformité des quantités étudiées, on choisit la dérivée #(t) selon le temps de l'ondelette sismique, au lieu de l'utilisation de l'ondelette sismique s(t) elle-même.
set) est la dérivée selon le temps de l'ondelette sismique s(t) = d[s(t)]/ dt ,
Une dérivée temporelle est généralement appliquée à chacune des traces sismiques avant la migration de Kirchhoff (en temps ou en profondeur). Ainsi, notamment dans un souci d'uniformité des quantités étudiées, on choisit la dérivée #(t) selon le temps de l'ondelette sismique, au lieu de l'utilisation de l'ondelette sismique s(t) elle-même.
Selon un autre mode de réalisation de l'invention, notamment lorsque les traces sismiques ne sont pas dérivées selon le temps, la fonction d'incrémentation prend en compte l'ondelette sismique s(t) au lieu de la dérivée #(t) de celle-ci
selon le temps ; on a alors i(td, tr ; s,r)= s(td (x; s, #) - tr (xr ( p); s, -)) .
selon le temps ; on a alors i(td, tr ; s,r)= s(td (x; s, #) - tr (xr ( p); s, -)) .
Un terme de correction d'amplitude de la migration w(x;s,r) est classiquement calculé lors de la migration, et prend notamment en compte une géométrie régulière d'acquisition.
Ce terme w(#;#,#) peut donc être avantageusement utilisé dans le procédé selon l'invention, et agir comme préconditionnement à la correction réalisée par la détermination du taux de couverture d'illumination précédemment décrite.
La fonction d'incrémentation peut alors être pondérée par ce terme de
préconditionnement, de sorte que le taux I(z,p;s,#) de couverture d'illumination relatif à une paire (source s, récepteur r ) est incrémenté par : (ta, tr; S")'(x;S#).
préconditionnement, de sorte que le taux I(z,p;s,#) de couverture d'illumination relatif à une paire (source s, récepteur r ) est incrémenté par : (ta, tr; S")'(x;S#).
L'évaluation du taux de couverture d'illumination dans le domaine migré selon l'invention est optimum puisque applicable à n'importe quelle géométrie d'acquisition.
Basé sur une approche inverse, le taux de couverture d'illumination I(#,#) estimé par l'invention concorde avec les résultats théoriques établis par Bleistein.
Cela n'est pas le cas par les autres méthodes proposées afin de résoudre un problème similaire.
<Desc/Clms Page number 17>
Comme cela a déjà été exposé précédemment, l'évaluation du taux de couverture d'illumination proposé par l'invention utilise seulement les quantités basiques disponibles pour la migration classiquement réalisée par sommation de Kirchhoff.
De manière avantageuse, l'évaluation du taux de couverture d'illumination selon l'invention peut être réalisée tout en effectuant la migration des traces sismiques.
Cela n'est pas le cas de certaines des autres méthodes proposées (Schneider et al. ; Laurain et Vinje) pour lesquelles ladite évaluation est réalisée post-migration.
Selon un autre aspect, l'invention propose également un procédé de correction des amplitudes des données sismiques enregistrées au cours d'une prospection sismique 3D afin de compenser de l'effet d'une illumination non uniforme des réflecteurs du sous-sol, caractérisé en ce qu'il comporte les étapes consistant à :
estimer le taux de couverture d'illumination I ( x, p ) par le procédé selon le premier aspect de l'invention détaillé ci-dessus, # utiliser l'inverse I-1(#,#) dudit taux comme facteur de pondération à appliquer à chacune desdites amplitudes des données sismiques.
estimer le taux de couverture d'illumination I ( x, p ) par le procédé selon le premier aspect de l'invention détaillé ci-dessus, # utiliser l'inverse I-1(#,#) dudit taux comme facteur de pondération à appliquer à chacune desdites amplitudes des données sismiques.
Selon encore un autre aspect, l'invention propose également un procédé pour la sélection, en fonction de la cible d'une prospection sismique 3D, d'une géométrie d'acquisition parmi une pluralité de géométries d'acquisition, caractérisé en ce qu'il comporte les étapes consistant à : déterminer, pour chacune des géométries d'acquisition considérées, le taux de couverture d'illumination I(#,#) estimé conformément au premier aspect de l'invention détaillé ci-dessus, sélectionner la géométrie d'acquisition procurant le taux de couverture d'illumination optimum en fonction de la cible
Enfin, on notera que le taux de couverture d'illumination I(#,#) estimé conformément au premier aspect de l'invention détaillé ci-dessus, peut aussi être
Enfin, on notera que le taux de couverture d'illumination I(#,#) estimé conformément au premier aspect de l'invention détaillé ci-dessus, peut aussi être
<Desc/Clms Page number 18>
utilisé après la migration, pour éliminer les empreintes d'acquisition, ainsi que pour déconvoluer l'image migrée.
Claims (14)
1. Procédé d'estimation du taux I(x,p ) de couverture d'illumination sismique dans le domaine 3D migré en au moins un point image x , pour au moins un pendage de vecteur p , caractérisé en ce que, pour chaque paire (source s , récepteur r ) du relevé sismique, on évalue le taux I(#,#;#,#) de couverture d'illumination relatif à ladite paire (source s , récepteur r ) par mise en #uvre des étapes de : # détermination du temps de trajet de réflexion tr(#r(#) ; s , r ) depuis la source s jusqu'au point de réflexion spéculaire x, sur le réflecteur plan passant par le point image x et perpendiculaire au vecteur p de pendage, puis de retour vers le réflecteur r ,
à partir du temps de trajet de diffraction td( x ; s , r ) depuis la source s jusqu'audit point image x , puis de retour jusqu'au réflecteur r ; # d'incrémentation dudit taux I(#,#;#,#) de couverture d'illumination relatif à ladite paire (source s, récepteur r ) en fonction de la différence entre le temps de trajet de diffraction td(#; s, r) et le
temps de trajet de réflexion tr( Xr ( p ) ; s , r ).
2. Procédé selon la revendication précédente, caractérisé en ce qu'il comporte en outre une étape de sommation de chacun desdits taux I(#,#;#,#) de couverture d'illumination relatif à une paire (source s, récepteur #), de manière à déterminer le taux de couverture d'illumination total I(x,p) = #I(#,#;#,#).
#. #
<Desc/Clms Page number 20>
3. Procédé selon l'une des revendications précédentes, caractérisé en ce que, au cours de l'étape d'incrémentation, le taux I(#,#;#,#) de couverture d'illumination est incrémenté en utilisant une fonction d'incrémentation i(td,
tr; s,# ), selon I ( x, p )=I ( x, p ) + i(td, tr; s,# ), ladite fonction d'incrémentation prenant en compte la différence entre le temps de trajet de diffraction td(x; s , r ) et le temps de trajet de réflexion tr(#r(#) ; s, r).
4. Procédé selon la revendication précédente, caractérisé en ce que la fonction d'incrémentation i est fonction de rondelette sismique s(t).
5. Procédé selon la revendication précédente, caractérisé en ce que la fonction d'incrémentation i s'exprime en fonction de rondelette sismique s(t) selon :
i(Íd, tr; s,r) = s(td(x;s,f)-tr(xr(p);s,f)).
6. Procédé selon la revendication 4, caractérisé en ce que la fonction d'incrémentation i s'exprime en fonction de la dérivée #(t) selon le temps de l'ondelette sismique s(t) selon :
i(td, tr; s,r) = Sltd('so') trlxrlp)'SY))'
7. Procédé selon l'une des revendications 3 à 6, dans lequel la migration met @ en jeu une correction a priori w(x;s,f) du taux de couverture d'illumination, caractérisé en ce que le taux I(#,#;#,#) de couverture d'illumination relatif à une paire (source s , récepteur r) est incrémenté de i(td, tr ;
s, r ). w(x; s, # ) .
<Desc/Clms Page number 21>
8. Procédé selon l'une des revendications précédentes, caractérisé en ce que l'étape de détermination comprend le développement en série de
Taylor, du second ordre, du temps de trajet de diffraction td(#; s, r) autour du point image x :
td \xrSYI -td ('xSY)+(vxtd lxsr))T 'lxr x)+ 2 (xr 'Xl.x,xtd lx,s,P . xr x
9. Procédé selon la revendication précédente, caractérisé en ce que le point de réflexion spéculaire xr(p) est déterminé le long dudit réflecteur de sorte que le temps de trajet de diffraction audit point réflexion spéculaire #r(#)
soit stationnaire, selon : pT n 0xtd (z;s, r) + x,xtd (z;s, ).(.z, (p) -x))= 0 .
10. Procédé selon les deux revendications précédentes, caractérisé en ce que le point de réflexion spéculaire xr et le temps de trajet de réflexion
tr( x, ( p ) ; s , r ) sont évalués selon les expressions suivantes: xr(p) = x-M.F-'.b; tr(xr(P) S, #)=ta(x s, r) - !. bT . F-1 . b , dans lesquelles, # M, une matrice (3 # 2) décrite par deux vecteurs s'étendant le long du réfelcteur, et donc perpendiculaires au vecteur p de pendage ; # b , un vecteur (2x1) des dérivés du premier ordre du temps de trajet de diffraction le long du plan de réflexion : b = MT .(#xtd) # F est une matrice (2 x 2) des dérivés du second ordre du temps de
trajet de diffraction le long du plan de réflexion : F = M T .( x,xtd ). M.
<Desc/Clms Page number 22>
11. Procédé selon l'une des revendications précédentes, caractérisé en ce que l'étape de détermination utilise les cartes des isochrones de migration td(#; s, r ) spécifiées pour toute paire (source s, récepteur r) impliquée dans la migration en tout point image x du domaine 3D migré.
12. Procédé selon l'une des revendications précédentes, caractérisé en ce que l'estimation du taux I(#,#) de couverture d'illumination sismique dans le domaine 3D migré est réalisé lors de la migration par sommation de
Kirchhoff des données sismiques enregistrées au cours de la prospection sismique 3D.
13. Procédé de correction des amplitudes des données sismiques enregistrées au cours d'une prospection sismique 3D afin de compenser de l'effet d'une illumination non uniforme des réflecteurs du sous-sol, caractérisé en ce qu'il comporte les étapes consistant à : # estimer le taux de couverture d'illumination I(#,#) par le procédé selon l'une des revendications 1 à 12 ; # utiliser l'inverse I-1(#,#) dudit taux comme facteur de pondération à appliquer à chacune desdites amplitudes des données sismiques.
4S
14. Procédé pour la sélection, en fonction de la cible d'une prospection sismique 3D, d'une géométrie d'acquisition parmi une pluralité de géométries d'acquisition, caractérisé en ce qu'il comporte les étapes consistant à : # déterminer, pour chacune des géométries d'acquisition considérées, le taux de couverture d'illumination I(#,#) par le procédé selon l'une des revendications 1 à 12 ; # sélectionner la géométrie d'acquisition procurant le taux de couverture d'illumination optimum en fonction de la cible.
Priority Applications (6)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR0308862A FR2858063B1 (fr) | 2003-07-21 | 2003-07-21 | Procede d'estimation du taux de couverture d'illumination dans le domaine migre |
| US10/565,403 US7406383B2 (en) | 2003-07-21 | 2004-07-20 | Method of estimating the illumination fold in the migrated domain |
| GB0602786A GB2419412B (en) | 2003-07-21 | 2004-07-20 | Method of estimating the illumination fold in the migrated domain |
| PCT/IB2004/002618 WO2005008292A1 (fr) | 2003-07-21 | 2004-07-20 | Procede d'estimation de l'interface d'illumination dans le domaine migre |
| CA2533012A CA2533012C (fr) | 2003-07-21 | 2004-07-20 | Procede d'estimation de l'interface d'illumination dans le domaine migre |
| NO20060818A NO337503B1 (no) | 2003-07-21 | 2006-02-20 | Fremgangsmåte for å estimere seismisk 3D dekningsgrad i det migrerte området |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR0308862A FR2858063B1 (fr) | 2003-07-21 | 2003-07-21 | Procede d'estimation du taux de couverture d'illumination dans le domaine migre |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| FR2858063A1 true FR2858063A1 (fr) | 2005-01-28 |
| FR2858063B1 FR2858063B1 (fr) | 2005-09-16 |
Family
ID=33560960
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR0308862A Expired - Fee Related FR2858063B1 (fr) | 2003-07-21 | 2003-07-21 | Procede d'estimation du taux de couverture d'illumination dans le domaine migre |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US7406383B2 (fr) |
| CA (1) | CA2533012C (fr) |
| FR (1) | FR2858063B1 (fr) |
| GB (1) | GB2419412B (fr) |
| NO (1) | NO337503B1 (fr) |
| WO (1) | WO2005008292A1 (fr) |
Families Citing this family (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007143355A2 (fr) | 2006-05-17 | 2007-12-13 | Nexus Geosciences Inc. | Traitement des données sismiques sur la base de diplets |
| US8194498B2 (en) * | 2008-01-18 | 2012-06-05 | Westerngeco L.L.C. | Using a wave propagator for transversely isotropic media |
| US8699298B1 (en) | 2008-06-26 | 2014-04-15 | Westerngeco L.L.C. | 3D multiple prediction and removal using diplets |
| US9329288B2 (en) | 2009-01-19 | 2016-05-03 | Landmark Graphics Corporation | Data acquisition and prestack migration based on seismic visibility analysis |
| WO2011041782A1 (fr) | 2009-10-02 | 2011-04-07 | Bp Corporation North America Inc. | Détermination d'éclairage basée sur une migration pour évaluation de risque à analyse de valeur moyenne |
| AU2011252076B2 (en) * | 2010-05-12 | 2014-04-10 | Shell Internationale Research Maatschappij B.V. | Seismic P-wave modelling in an inhomogeneous transversely isotropic medium with a tilted symmetry axis |
| IT1404170B1 (it) * | 2011-02-10 | 2013-11-15 | Eni Spa | Metodo di indagine sismica del sottosuolo |
| US9025414B2 (en) | 2011-05-27 | 2015-05-05 | Conocophillips Company | Reciprocal method two-way wave equation targeted data selection for seismic acquisition of complex geologic structures |
| CA2837649A1 (fr) | 2011-09-28 | 2013-04-04 | Conocophillips Company | Procede reciproque de selection de donnees cibles d'equation d'onde a deux voies pour acquisition sismique de structures geologiques complexes |
| WO2013048585A1 (fr) | 2011-09-28 | 2013-04-04 | Conocophillips Company | Sélection ciblée de données par méthode inverse avec équation d'onde bidirectionnelle pour une imagerie améliorée de structures géologiques complexes |
| US9733374B2 (en) | 2012-12-17 | 2017-08-15 | Cgg Services Sas | Repeatability indicator based on shot illumination for seismic acquisition |
| EP2743737A2 (fr) | 2012-12-17 | 2014-06-18 | CGG Services SA | Procédés et systèmes pour le contrôle de la qualité de cartes d'éclairage sismique |
| CN104237939B (zh) * | 2014-09-12 | 2016-10-26 | 中国科学院地质与地球物理研究所 | 一种三维地震观测***多频率偏移噪声获取方法及装置 |
| CN104407383B (zh) * | 2014-12-09 | 2017-01-18 | 中国石油天然气集团公司 | 三维地震照明分析并行实现方法 |
| US10353092B2 (en) * | 2015-12-10 | 2019-07-16 | Pgs Geophysical As | Velocity model update with an inversion gradient |
| CN105467445B (zh) * | 2015-12-15 | 2018-02-02 | 中国石油天然气集团公司 | 一种三维地震观测***的建立方法及装置 |
| CN113866821B (zh) * | 2021-09-26 | 2022-08-02 | 吉林大学 | 一种基于照明方向约束的被动源干涉偏移成像方法和*** |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1995015506A1 (fr) * | 1993-12-01 | 1995-06-08 | Schlumberger Technology Corporation | Procede d'introduction de donnees sismiques et d'equilibrage du volume de calcul de ces donnees a des fins de traitement sur un ordinateur multiprocesseur |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2821677B1 (fr) * | 2001-03-05 | 2004-04-30 | Geophysique Cie Gle | Perfectionnements aux procedes d'inversion tomographique d'evenements pointes sur les donnees sismiques migrees |
| US6665618B1 (en) * | 2002-08-14 | 2003-12-16 | Conocophillips Company | Seismic survey design technique |
-
2003
- 2003-07-21 FR FR0308862A patent/FR2858063B1/fr not_active Expired - Fee Related
-
2004
- 2004-07-20 GB GB0602786A patent/GB2419412B/en not_active Expired - Fee Related
- 2004-07-20 WO PCT/IB2004/002618 patent/WO2005008292A1/fr active Application Filing
- 2004-07-20 US US10/565,403 patent/US7406383B2/en not_active Expired - Lifetime
- 2004-07-20 CA CA2533012A patent/CA2533012C/fr not_active Expired - Fee Related
-
2006
- 2006-02-20 NO NO20060818A patent/NO337503B1/no not_active IP Right Cessation
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1995015506A1 (fr) * | 1993-12-01 | 1995-06-08 | Schlumberger Technology Corporation | Procede d'introduction de donnees sismiques et d'equilibrage du volume de calcul de ces donnees a des fins de traitement sur un ordinateur multiprocesseur |
| US5530646A (en) * | 1993-12-01 | 1996-06-25 | Schlumberger Technology Corporation | Method for inputting seismic data for seismic processing on a multiprocessor computer |
Non-Patent Citations (2)
| Title |
|---|
| BLEISTEIN N. ET AL.: "True-amplitude transformation to zero offset of data from curved reflectors", GEOPHYSICS,SOCIETY OF EXPLORTION GEOPHYSICISTS.TULSA,US, vol. 64, no. 1, 1 January 1999 (1999-01-01) - 28 February 1999 (1999-02-28), pages 112 - 129, XP002276525 * |
| BLEISTEIN: "On the imaging of reflectors in the earth", GEOPHYSICS, SOCIETY OF EXPLORATION GEOPHYSICISTS. TULSA, US, vol. 52, no. 7, 1 July 1987 (1987-07-01), pages 931 - 942, XP002095154, ISSN: 0016-8033 * |
Also Published As
| Publication number | Publication date |
|---|---|
| US20060190179A1 (en) | 2006-08-24 |
| NO337503B1 (no) | 2016-04-25 |
| CA2533012A1 (fr) | 2005-01-27 |
| CA2533012C (fr) | 2014-08-26 |
| US7406383B2 (en) | 2008-07-29 |
| FR2858063B1 (fr) | 2005-09-16 |
| GB0602786D0 (en) | 2006-03-22 |
| GB2419412A (en) | 2006-04-26 |
| WO2005008292A1 (fr) | 2005-01-27 |
| GB2419412B (en) | 2006-11-01 |
| NO20060818L (no) | 2006-04-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP0789849B1 (fr) | Methode pour representer dans un domaine spatio-temporel la trajectoire d'au moins un puits de forage | |
| FR2858063A1 (fr) | Procede d'estimation du taux de couverture d'illumination dans le domaine migre | |
| EP1624321B1 (fr) | Méthode pour construire un modèle d'un milieu hétérogène décrit par plusieurs paramètres à partir de données exprimées dans des échelles de temps différentes | |
| CA2528700C (fr) | Methode pour determiner des informations speculaires apres imagerie sismique avant sommation | |
| CA2670384C (fr) | Methode d'inversion conjointe de donnees sismiques representees sur des echelles de temps differentes | |
| EP3023815B1 (fr) | Procédé de construction d'un modele geologique | |
| EP2253970B1 (fr) | Méthode pour imager une zone cible du sous-sol à partir de données de type walkaway | |
| Wassermann et al. | Toward a single‐station approach for microzonation: Using vertical rotation rate to estimate Love‐wave dispersion curves and direction finding | |
| FR2872584A1 (fr) | Methode pour simuler le depot sedimentaire dans un bassin respectant les epaisseurs des sequences sedimentaires | |
| FR2909775A1 (fr) | Methode de construction d'un modele geologique d'une formation du sous-sol contraint par des donnees sismiques | |
| FR3027688A1 (fr) | ||
| FR2918178A1 (fr) | Methode pour ajuster un modele de vitesse d'ondes sismiques en fonction d'informations relevees aux puits | |
| Wang et al. | Regularized elastic passive equivalent source inversion with full-waveform inversion: Application to a field monitoring microseismic data set | |
| Greer et al. | Matching and merging high-resolution and legacy seismic images | |
| FR2735244A1 (fr) | Methode de traitement pour obtenir des donnees sismiques a deport nul par sommation dans le domaine profondeur | |
| Yablokov et al. | Uncertainty quantification of multimodal surface wave inversion using artificial neural networks | |
| FR2971859A1 (fr) | Procede et dispositif de lissage a preservation du temps de trajet | |
| CA2364938C (fr) | Methode pour obtenir des temps de trajet en reflexion a partir d'une interpretation de donnees sismiques en ondes cylindriques migrees | |
| Liu et al. | Characterizing near-surface structures at the Ostia archaeological site based on instantaneous-phase coherency inversion | |
| Song et al. | A reflection-based efficient wavefield inversion | |
| FR2858064A1 (fr) | Procede de pointe bispectral des parametres de correction d'obliquite anelliptique | |
| Egan et al. | The influence of spatial sampling on resolution | |
| FR2691811A1 (fr) | Procédé de détermination d'un modèle de distribution optimale des vitesses de propagation d'ondes sismiques dans le sous-sol. | |
| Mikhailov et al. | Practical approach to joint imaging of multicomponent data | |
| FR2796726A1 (fr) | Procede de reconstruction tridimensionnelle d'une grandeur physique dans l'environnement d'un forage |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| CA | Change of address | ||
| CD | Change of name or company name | ||
| TP | Transmission of property | ||
| PLFP | Fee payment |
Year of fee payment: 14 |
|
| ST | Notification of lapse |
Effective date: 20180330 |