FR2824917A1 - Systeme pour l'estimation de variations azimutales dans des donnees sismiques - Google Patents
Systeme pour l'estimation de variations azimutales dans des donnees sismiques Download PDFInfo
- Publication number
- FR2824917A1 FR2824917A1 FR0205945A FR0205945A FR2824917A1 FR 2824917 A1 FR2824917 A1 FR 2824917A1 FR 0205945 A FR0205945 A FR 0205945A FR 0205945 A FR0205945 A FR 0205945A FR 2824917 A1 FR2824917 A1 FR 2824917A1
- Authority
- FR
- France
- Prior art keywords
- seismic data
- traces
- trace
- cluster
- time
- 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
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. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- 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. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/53—Statics correction, e.g. weathering layer or transformation to a datum
Abstract
L'invention comprend un système destiné à traiter des données sismiques pour estimer un décalage temporel dû à une anisotropie de vitesse dans le sous-sol terrestre. Un regroupement de traces de données sismiques est formé et des traces de données sismiques sélectionnées incluses dans ledit regroupement à l'intérieur de fenêtres temporelles sélectionnées sont inter-corrélées afin d'estimer le décalage temporel dans les traces de données sismiques incluses dans ledit regroupement dû à une anisotropie de vitesse dans le sous-sol terrestre.
Description
ce qu'il comprend un déboguage de la puce à semi-conducteur.
SYSTEME POUR L'ESTIM:ATION DE VARIATIONS AZIMUTALES
DAN S DES DONNEES SISMIQUES
Cette invention concerne le traitement de données sismiques. Plus précisément, l'invention porte sur un système de traitement des données sismiques pour détecter des
variations azimutales de vitesse.
Des études sismiques sont utilisées de manière routinière dans la recherche de gisements de pétrole et de gaz dans le sous-sol terrestre. Les études sismiques sont réalisées par transiert d'énergie acoustique dans le sol, soit au niveau d'une surface terrestre soit en environnement marin, puis détection de l'énergie acoustique réfléchie et réfractée. Le retard entre le transtert de l'onde d'énergie acoustique au niveau de la position de source et la détection de la même onde au niveau d'une position de récepteur
indique la profondeur des interfaces géologiques réfléchissantes.
Jusqu'à récemment, seules des études sismiques bidimensionnelles ("2D") étaient réalisées, les positions des sources sismiques étant colinéaires avec une ligne de récepteurs. De récents progrès techniques ont permis de collecter et d' analy ser des données d'études sismiques tridimensionnelles ("3D"). Typiquement, dans des études 3D, des réseaux de récepteurs sismiques sont déployés qui reçoivent une énergie acoustique réfléchie ayant été communiquée en diverses positions pouvant étre choisies
spécifiquement pour fournir un riche assortiment d'azimuts pour une couverture multiple.
Une technique fréquemment utilisée dans l'analyse d'études sismiques est l'analyse VAD, qui signifie variations d'amplitude avec le déport, et sont aussi appelées variations d'amplitude avec l'angle d'incidence dans la présente. Conformément à l'approche VAD, des attributs d'une interface souterraine sont déterminés à partir à la fois de l'amplitude de l'énergie sismique réfléchie sous incidence normale, et de la dépendance des réflexions sismiques détectées à l'angle d'incidence-de l'énergie sismique sur une interface souterraine réfléchissante par rapport à la verticale. Dans une analyse VAD conventionnelle, de multiples traces sismiques ayant un point de réflexion commun sont regroupées, dont l'ensemble est couramment appelé regroupement de traces à milieu commun ou en position miroir ou couverture multiple (CMP ou CDP). Du regroupement par milieu commun (ou en couverture multiple), on peut déduire l' amplitude R d'une onde sismique réfléchie par une interface (c'est-à-dire "l'horizon cible") en fonction de l'angle d'incidence par rapport à la normale, conformément à la relation suivante: R() = A+Bsin20 Dans ce cas, le coefficient A est la réponse à déport nul (également appelée ordonnée VAD à l'origine), tandis que le coefficient B est appelé pente ou gradient VAD, puisqu'il
est représentatif du taux de variation de l'amplitude avec le carré de l'angle d'incidence.
Une analyse de la pente VAD et de l'ordonnée VAD à l'origine peut donner des indications de formations intéressantes, du point de vue de la prospection de pétrole et de gaz. Par exemple, des variations des valeurs de A et B par rapport à une droite de tendance théorique de A par rapport à B pour les séquences stratigraphiques attendues peuvent
indiquer l'emplacement de gisements d'hydrocarbures.
Bien que des modèles simples de géologie souterraine supposent une isotropie azimutale dans la propagation d'énergie acoustique, on a observé en fait la présence d'une anisotropie azimutale dans de nombreuses régions d'études, si bien que la vitesse de l'énergie acoustique dépend de l'azimut du trajet source-récepteur. En présence d'une anisotropie azimutale, la correction d'obliquité conventionnelle ne peut pas aligner correctement les traces sismiques du regroupement, ce qui peut entramer une diminution
de qualité de l'analyse VAD.
Une correction d'obliquité des données sismiques, destinée à compenser à la fois des retards dépendant du déport et une anisotropie azimutale due aux terrains de recouvrement, est donc typiquement effectuée dans la production de traces sommées à rapport signal à bruit amélioré utilisées dans une étude sismique 3D. Par exemple, le brevet des Etats-Unis n 5 5 3 2 97 8 décrit un procédé de calcul et d' application de
corrections d'anisotropie azimutale à des signaux d'études sismiques.
La détection d'une direction azimutale préférée au niveau d'une interface réfléchissante peut aussi fournir des informations importantes relatives à des particularités géologiques. Par exemple, une direction azimutale de réflexion prétérée peut indiquer la présence de fissures verticales alignées. Pour des déports modérément grands (angles d'incidence de 25 35 ), l'onde longitudinale se propageant dans le plan parallèle aux fissures verticales alignées a une vitesse plus élevée que l'onde longitudinale se
propageant dans le plan perpendiculaire aux fssures.
De manière traditionnelle, une analyse azimutale de vitesse était effectuée par utilisation de super regroupements sectorisés en azimut et extraction de maxima de similitude à divers azimuts. Cela réduit le problème à une série de solutions 2D, au lieu de résoudre le problème 3D complet. Dans certains cas, il est possible de ne choisir que deux secteurs, perpendiculaire et parallèle aux axes principaux (moyens) de l'anisotropie azimutale. Si plus de deux secteurs sont utilisés, une ellipse est ajustée aux vitesses trouvoes pour donner les valeurs de grande vitesse et petite vitesse et l'azimut de la grande
vitesse. Ces procédures souffrent de plusieurs inconvénients.
L'extraction de similitude, à la main, de données sectorisées en azimut dépend du processeur/interpréteur et est extrêmement coûteuse en temps. La similitude fonctionne bien pour des données qui ne présentent pas de variations d'amplitude avec le déport ("VAD"), toutefois, si les données contiennent des VAD
significatives, en particulier s'il existe une inversion de polarité, la similitude peut échouer.
Dans ce cas, l'extraction automatique des maxima de similitude va êke erronée.
Si le sous-sol présente des variations azimutales de vitesse ("VAV"), alors elles vont apparaître comme une correction statique dépendant du déport vue sur des regroupements par milieu commun triés par déport. Cela va réduire l'efficacité de n'importe quelle solution statique compatible en surface, donc les super regroupements
sectorisés en azimut vont très probablement être contaminés par une correction statique.
Cela va abaisser significativement la qualité de l'analyse de similitude et peut entrâîner la
présence de plusieurs maxima de similitude.
La similitude est basée sur le fait de donner le plus grand résultat de sommation.
Toutefois, pour une analyse VAV on s'intéresse à la vitesse souterraine réelle et pas simplement à la vitesse qui donne la meilleure somme. Par exemple, si une plus grande amplitude apparâit en un azimut particulier à l'intérieur du secteur, alors la vitesse à cet azimut va étre extraite. Si de plus ces grandes amplitudes sont aux déports moyens à petits et sont contaminces par une correction statique résiduelle, alors une vitesse complètement
erronée peut donner la plus grande sirnilitude.
La sectorisation et la sommation partielle des données font qu'il est extrêmement 2 5 difficile d' obtenir des estimations d' erreurs. Non seulement il est difficile d'attribuer une erreur d' exkaction à partir de l'extraction de similitude, mais de plus, de s erreurs dues à la géomékie d' acquisition ne sont pas représentées. Dans n'importe quelle analyse de ce type, il est important de calculer les erreurs associées aux résultats obtenus. Par exemple, une approche par moindres carrés pondérés a été utilisée pour calculer les erreurs dans une technique destince à inverser les variations azimutales d'amplitude pour des données d'onde transversale. On a aussi observé que la fiabilité de l'analyse des variations d'amplitude avec l'azimut était évaluée par recherche d'une absence de reflet de la
géoméhie d'acquisition dans les cartes d'anisokopie.
On notera que la de scription suivante de l' invention ne do it pas être co mpri se comme limitant l'invention aux exemples et formes de réalisation préférées qui sont représentés et décrits. Les personnes compétentes dans l'art auquel appartient cette invention seront capables de concevoir des variantes de cette invention restant la portée
des revendications annexées.
L'invention comprend un système destiné à traiter des données sismiques pour estimer un décalage temporel dû à une anisokopie de vitesse dans le sous-sol terrestre. Un regroupement de traces de données sismiques est formé et des traces de données sismiques sélectionnées incluses dans ledit regroupement à l'intérieur de fenêtres temporelles sélectionnées sont inter-corrélées afn d'estimer le décalage temporel présent dans les traces de données sismiques incluses dans ledit regroupement en résultat d'une anisotropie
de vitesse dans le sous-sol terrestre.
L ' invention sera mieux compri se à l' étude de quelques exemples ne limitant pas l'invention et illustrés par les dessins joints, sur lesquels: la figure 1 représente une partie simplifiée d'une configuration sources-récepteurs sismiques 3D pour une étude sismique 3D; la figure 2A représente des traces de données sismiques représentatives provenant d'un regroupement par milieu commun avant application d'un ajustement d'obliquité; la figure 2B représente les traces de données sismiques de la figure 2A après application d'un ajustement d'obliquité; la figure 3 est un organigramme représentant une forme de réalisation de l' invention, la figure 4 est un organigramme représentant une autre forme de réalisation de l'invention, la fgure 5 est un organigramme représentant une autre forme de réalisation de l' invention; la figure 6 est un organigramme représentant une autre forme de réalisation de l' invention; la figure 7 est un organigramme représentant une autre forme de réalisation de l'invention, la figure 8 est un schéma représentant une relation spatiale de l'invention; et la figure 9 représente un système informatique destiné à mettre en _uvre
l' invention.
La figure 1 représente une partie simplifiée d'une configuration sourcesrécepteurs sismiques 3D pour une étude sismique 3D. La figure 1 représente une partie d'un réscau de récepteurs 12, placé sur la surface terrestre, comprenant trois colonnes de récepteurs, chaque colonne comprenant huit récepteurs. Un réseau de sources 14 comprend un groupe de positions de source. Typiquement, une source sismique est déplacée le long de la surface terrestre et la source est activée, en séquence, en des positions de source spécifiques. L'énergie acoustique transiérée en chaque position de source se propage dans le sol et, après réflexion par des interfaces géologiques souterraines, est détectée par
chaque récepteur du réseau de récepteurs.
La figure 1 représente aussi un exemple d'un milieu commun MP. Le milieu commun représenté est commun à plusieurs trajets source-récepteur dans cette étude. La figure 1 représente MP comme étant un milieu entre une position de source S11 et une position de récepteur R38; une position de source S21 et une position de récepteur R28; une position de source S38 et une position de récepteur R11; et une position de source S36 et une position de récepteur R13. Les personnes ayant des compétences ordinaires dans l' art comprendront que la position MP pourrait être le milieu commun d'un grand nombre de paires source-récepteur, et aussi qu'il existe un grand nombre de milieux communs entre d'autres paires source- récepteur. Chacun de ces trajets source-récepteur peut être d'une longueur (déport) différente et d'une direction (azimut) différente. Bien que des regroupements par milieu commun soient typiquement utilisés dans la mise en _uvre de la présente invention, si les couches souterraines sont inclinées plutôt qu'horizontales, un traitement supplémentaire connu des personnes ayant des compétences ordinaires dans l' art peut être exécuté sur les données sismiques pour réaliser des regroupements par point
de réflexion commun auxquels la présente invention peut être appliquée.
De manière bien connue dans l'art, des corrections d'obliquité sont typiquement apportées aux traces d'un regroupement par milieu commun afin de corriger le retard supplémentaire correspondant à des traces à plus grand déport, de sorte que le temps de propagation de chaque signal sismique est effectivement normalisé en un temps de propagation à déport nul. Dans des situations dans lesquelles le sol ne présente pas d'anisotropie azimutale, les variations d'azimut ne vont introduire aucune variation dans l'ensemble de donnces. Toutefois, en présence d'une anisotropie azimutale dans des signaux d'études sismiques, la correction d'obliquité conventionnelle ne peut pas corriger de manière adéquate les variations des retards pour des traces provenant de paires source
récepteur de directions différentes.
Des variations dans les traces de données sismiques vont apparâître comme étant des variations d'amplitude en fonction des variations d'azimut, alors qu'en réalité les variations dans les traces de données sismiques sont le résultat de décalages temporels qui dépendent de l'azimut. Une variation de vitesse avec l'azimut de quelques pour-cent seulement peut entrâîner une différence temporelle de dix millisecondes ou plus dans la position d'un phénomène dans une trace de données sismiques. En conséquence, des variations de vitesse inférieures à la résolution de l'analyse conventionnelle de vitesse basée sur la similitude peuvent déformer les variations d' amplitude avec l' angle d'incidence (VAD) et les variations d'amplitude avec l'azimut ("VAA"), enhâînant un
l0 calcul incorrect.
En plus de la correction d'obliquité conventionnelle, tvpiquement, pour une acquisition au sol 3D, une déconvolution, une correction statique de réfraction et deux passages de traitements statiques résiduels compatibles en surface sont appliqués aux données. En préliminaire à l'application de l' invention telle que décrite dans la présente, l 5 une réduction du bruit, une mise à l'échelle de trace et n'importe quel autre kaitement sur une seule trace peuvent étre appliqués, toutefois l'application d'un kaitement sur plusieurs
traces doit être évitée.
MISE A PLAT DE PHENOMENES A L'INTERIEUR D'UNE FENETRE
TEMPORELLE MOBILE
Si la vitesse des signaux sismiques dans le sous-sol varie avec la direction azimutale dans laquelle les signaux sismiques se propagent, la correction d'obliquité conventionnelle ne va pas aligner correctement les traces, et selon une forme de réalisation préférée de l' invention, le traitement suivant peut être utilis é pour obtenir un alignement correct des traces. Les figures 2A et 2B représentent des traces provenant d'un regroupement par milieu commun. Un profil de vitesse initiale pour les traces sismiques peut être généré d'une manière conventionnelle et une correction d'obliquité appliquée aux traces. La figure 2A représente les traces avant application de l'ajustement d'obliquité aux traces et la figure 2B représente les mêmes traces après application de l'ajustement d'obliquité. De manière indiquée sur la figure 3, et à propos de la figure 2B, dans une forme de réalisation prétérée de l'invention une première fenêtre temporelle est sélectionnée à une étape 50. Les fenêtres temporelles utilisées dans la mise en oeuvre de l'invention peuvent typiquement être incluses dans l'intervalle de l 00 à 300 millisecondes, toutefois la fenêtre temporelle utilisée peut varier conformément au jugement du processeur. Une première fenêtre temporelle sélectionnée pour les donnces représentées sur la figure 2B peut s'étaler de 500 à 700 millisecondes. A une étape 52, une fenêke spatiale comprenant un certain nombre de traces en relation spatiale est sélectionnée, et les traces dans cette fenêtre spatiale sont sommées à l'intérieur de la fenêtre temporelle sélectionnée afin de créer une trace "pilote". A une étape 54, une seule trace est sélectionnée à partir de la trace pilote, qui sera appelée la trace "d'entrée". A une étape 56, la trace pilote sommée est intercorrélée à la trace d' entrée à l'intérieur de la fenêtre temporelle sélectionnée, et à une étape 58 le décalage temporel entre la trace d'entrée et la
trace pilote qui donne la corrélation maximale est déterminé et enregiské.
La fenêke spatiale est ensuite déplacée d'une trace à travers le regroupement par milieu commun, à l'intérieur de la fenêtre temporelle sélectionnce, et une nouvelle trace pilote est générce pour cette nouvelle fenêke spatiale. La nouvelle kace pilote est ensuite corrélée à la nouvelle kace d'entrce à l'intérieur de la fenêtre temporelle. Cette opération se poursuit pour la fenêtre temporelle sélectionnée, chaque trace de données à l'intérieur du regroupement étant successivement désignée comme "trace d'entrée" et inter-corrélée à une trace pilote qui comprend une pluralité de traces proches, afin d'achever une séquence de corrélation de traces en fenêke temporelle. De cette manière, la kace pilote représente les caractéristiques locales de phase et d'amplitude des données pour chaque trace "d'enkée". En conséquence, il est décidé à une étape 60 si une séquence de corrélation de traces en fenêtre temporelle est achevée ou non. Si la réponse est non, les étapes 52, 54, 56, 58 et 60 sont répétées pour une kace d'entrée successive. Si la réponse est oui, alors à une étape 62 il est décidé si la séquence d'inter-corrélation tout juste achevée est ou n'est pas la première opération d'inter-corrélation effectuée pour la fenétre temporelle. Dans une mise en ceuvre de l' invention, si la réponse est oui, le décalage temporel calculé ayant donné une inter-corrélation maximale enke la trace d'entrée et la trace pilote pour chaque trace est appliqué à chaque trace à une étape 64 et les étapes 52, 54, 56, 58, 60 et 62 sont répétées. Si la réponse à l'étape 62 est non, alors à une étape 66 il est décidé si n'importe quel phénomène de réflexion dans la fenêtre temporelle est sensiblement aligné parmi les traces du regroupement. Si la réponse est non, les étapes 64, 52, 54, 56, 58, 60, 62 et 66 sont répétées. Cette décision à l'étape 66 est normalement basée sur le fait que les décalages temporels supplémentaires calculés pour les traces d'entrce dans la séquence de corrélation de traces en fenêtre temporelle tout juste achevée sont significatifs ou non. Si la réponse à l'étape 66 est oui, alors à une étape 68 les décalages temporels calculés totaux pour chaque trace sont mémorisés et les kaces sont remises sous la forme qu'elles avaient avant le début du kaitement décrit sur la figure 3. Typiquement, deux itérations seulement du traitement décrit aux étapes 52, 54, 56, 58, 60, 62 et 64 sont effectuées, mais d'autres itérations peuvent être effectuées si, au jugement du processeur, d'autres itérations peuvent
améliorer la qualité des donnces.
Il est décidé à une étape 70 si toutes les fenêtres temporelles pertinentes dans le S regroupement de données sismiques ont été sélectionnées ou non, et si la réponse est non
une nouvelle fenêtre temporelle est sélectionnée à l'étape 50 et la procédure d'inter-
corrélation décrite plus haut à propos des étapes 52, 54, 56, 58, 60, 62, 64, 66, 68 et 70 est répétée pour toutes les fenêtres temporelles pertinentes. Les fenêtres temporelles successives sélectionnées peuvent occuper des positions temporelles successives sur les lO traces de données sismiques ou les fenêtres temporelles peuvent se chevaucher, selon la
qualité des données.
Après achèvement de la procédure d'inter-corrélation pour toutes les fenêtres temporelles et toutes les traces à l'intérieur de chaque fenêtre temporelle, à une étape 72 la valeur de décalage temporel ayant donné la corrélation maximale pour chaque trace à l 5 l'intérieur de chaque fenêtre temporelle est appliquce à chaque trace au milieu de chaque fenêtre temporelle, et, à une étape 74, des décalages temporels pour le reste des traces de
données sont interpolés entre ces milieux.
Chaque trace pilote peut comprendre, par exemple, onze traces. La trace située au centre de la fenêtre spatiale sélectionnée, c' e st-à- dire la sixième trace, peut être dé signée comme trace "d'entrée" et intercorrélée à la trace pilote afin d'obtenir un décalage temporel. La fenêtre spatiale est ensuite déplacée d'une trace à travers le regroupement et une nouvelle trace pilote est formée. De nouveau, la trace située au centre de cette fenêtre, c'est-à-dire la trace suivante dans le regroupement, est désignée comme trace "d'entrée" et inter- corrélée à la trace pilote afin d'obtenir un décalage temporel pour cette trace. Au niveau des bords du regroupement, dans cet exemple les première à sixième traces, la fenêtre spatiale constituant la trace pilote peut être raccourcie de sorte que, si la première trace est la trace "d'entrce", les première à sixTème traces sont sommées pour former la trace pilote, et pour la deuxième trace, les première à septième traces sont sommées et ainsi de suite, jusqu'à ce que le nombre total de traces désirées dans la fenêtre spatiale soit atteint (dans cet exemple, onze). Le nombre de traces sélectionnées pour former la trace pilote peut être choisi sur la base de la valeur des variations d'amplitude avec le déport
dans les données et de la valeur de bruit tel qu'une contamination multiple.
D ans l' exécution de la procédure d'inter- corrélation, de s variations d'amplitude et de phase avec le déport sont prises en compte, incluant le cas dans lequel un phénomène change de polarité pour un certain déport. La trace pilote représente les caractéristiques "locales" des donnces. Si un phénomène change de polarité au niveau de grands déports, alors la trace pilote aux petits déports ne doit pas inclure les traces aux grands déports. De manière similaire, la trace pilote aux grands déports n'inclura pas de traces aux petits déports. D'autres attributs des traces, incluant des valeurs absolues, des valeurs moyennes quadratiques (RMS) ou une enveloppe de trace, peuvent étre utilisés pour effectuer l'inter
corrélation en plus de l'amplitude de réflexion brute de la trace.
Une fois que les décalages temporels ont été appliqués aux traces de données l0 sismiques du regroupement, une analyse VAD ainsi qu'une analyse VAA, telles que décrites dans la présente à propos de la figure 4, peuvent étre effectuces sur les traces ajustées.
CALCUL DES VARIATIONS D'AMPLITUDE AVEC L'ANGLE D'INCIDENCE
(VAD) ET DES VARIATIONS D'AMPLITUDE AVEC L'AZIMUT (VAA)
lS Les personnes ayant des compétences ordinaires dans l'art savent qu'il est possible d'approcher les variations d'amplitude avec l'angle d'incidence (également appelées variations d' amplitude avec le déport), ainsi que les variations d' amplitude avec l' azimut, pour une réflexion par une couche isotrope transversale horizontale du sol qui est recouverte de terrain isotrope, par l'équation suivante: R(0,O=+Gisin2H+ G2sin20cos2(3) (Eq. l) dans laquelle: R(0, f) est le coefficient de réflexion en fonction de 0, l'angle d'incidence de l' énergie sismique sur l' interface réfléchis sante souterraine par rapport à la verticale, et de |, l'azimut du récepteur par rapport à une direction prédéLinie d'azimut zéro (par exemple le nord géométrique); I est le contraste d'impédance d'onde longitudinale entre les couches souterraines par lesquelles le signal est réfléchi; G est le gradient VAD isotrope; G2 est le terme azimutal ou anisotrope; et (voir la fgure 8) est l'angle entre la direction prédéfmie d'azimat zéro (telle que
le nord géométrique) et la direction du gradient VAD maximal.
G et G2 sont donnés par: (; 2 V 2(g) ( VT p (Eq. 2) et G2=/v>+2(l) Ay (Eq. 3) dans lesquelles: Ap, AV et AVT sont respectivement la différence de densité, la différence de vitesse d' onde longitudinale, et la différence de vitesse d' onde transversale, p, V, et VT sont respectivement la densité moyenne, la vitesse moyenne d'onde long itudinale et la vites s e moy enne d' onde transvers ale, g =-est la vitesse moyenne d'onde longitudinale divisée par la vitesse moyenne d'onde kansversale, l 0 o(V) est la différence de 6(V) d'un côté à l'autre de la frontière réfléchissante, et Ay est la différence de paramèhe de séparation d'onde transversale y d'un côté à l'autre de la frontière réfléchissante, avec: y=G22-C3232
2C3232
Les personnes ayant des compétences ordinaires dans l'art savent que pour une matière l 5 linéairement élastique, chaque composante de contrainte aij dépend linéairement de toutes les composantes de déformation EkI, dans lesquelles i, j, k et l sont des indices directionnels qui peuvent prendre les valeurs l, 2 ou 3. La dépendance contrainte déformation est donnée par la loi de Hooke: ()ij = Cijklekl dans laquelle Cijk est le tenseur de module d' élasticité et caractérise entièrement l' élasticité du milieu. La relation entre 6(V) et le tenseur de module d'élasticité est donnée par: (Cl133-C3232)2-(C3333-C3232)2
2C3333(C3333-C3232)2
Toutefois, sans connaissance de 0, il n'est pas possible de résoudre l'équation l par une approche par moindres carrés. On peut cependant réécrire l'équation 1 sous la forme: R(0,)=I+{Gi+G2cos2()}sin20 (Eq. 4) qui peut se réécrire sous la forme: R(0,)=I+{GI +(G2*-Gl*)cos2(-} sin20 (Eq. 5) en posant G = G et G2 = G2 - G. En utilisant l'égalité: Gl+ (G2-Gl*)cos2( - O=G2cos2(3)+Glsin2(B) (Eq. 6) il vient: R(0,)=I+{G2cos2(O+Gsin2(-/3)}sin20 (Eq. 7) Les personnes ayant des compétences ordinaires dans l' art savent que: G2 cos2 (f - p) + G sin2 (I - p) = W cos2 + 22 cosi sin + W22 sin2 (Eq. 8) qui est linéaire par rapport aux inconnues Wi, W2 et Wi3, qui peuvent être reliées aux inconnues G, G2 et, de la manière suivante: G2=0,5(Wi+W2+(Wli-W22)2+4Wl22) (Eq. 9) G:=0,5(Wi.+W24(WW22)2+4W22) (Eq. 10)
W-W22+(W-W22)2+4W22 1 (E. 1 1)
B=Arctg 2Wi2 J q Par conséquent, en combinant l'équation 1 avec les équations 4-11, l'équation 1 peut s'écrire sous la forme: R(0,O=I+[Wcos2+ 2W2cossin(+W22sin2sin20 (Eq. 12) avec (=o'5(Wi+Wi2-(w-w22)2+4w22) (Eq. 13) G2=(Wll-W22)2+4WI22 (Eq. 14) 1S et =ArC'gW-w22+>l(wi-w22)2+4w22 (Eq. 15) Les valeurs du coefficient de réflexion R(H, O pour des valeurs spécifques de l'angle d'incidence et de l'azimut source-récepteur peuvent étre obtenues à partir des données sismiques enregistrées pour chaque phénomène de réflexion, par extraction des amplitudes des kaces de données sismiques en fonction du déport et de l'azimut. A propos de la figure 4, à une étape 80, des valeurs du coefficient de réflexion R(H, 4) et de l'azimut
source-récepteur sont obtenues à partir des données sismiques en cours de traitement.
Pour obtenir la valeur de l'angle d' incidence 0, une version lissée de la vitesse d'intervalle est calculée à une étape 82 d'une manière bien connue des personnes ayant des compétences ordinaires dans l'art, et la vitesse moyenne quadratique est calculée à une étape 84 à partir de la version lissée de la vitesse d'intervalle. A une étape 86, les valeurs de l'angle d'incidence peuvent étre déterminées en utilisant l'équation suivante:H=ArcsiniVnt X l (Eq. 16) lVrms X2+To2Vrrs2 J dans laquelle: X est le déport source à récepteur; To est le temps de propagation aller- retour à déport nul; V=s est la vitesse moyenne quadratique; et
V est la vites se d' intervalle à l' instant pertinent.
A une étape 88, un procédé des moindres carrés est utilisé pour calculer un coefficient de réilexion en fonction de l' angle d' azimut et de l'angle d'incidence pour les traces sismiques constituant le regroupement par milieu commun. L'équation 12 est résolue d'une manière directe par moindres carrés, connue des personnes ayant des compétences ordinaires dans l'art, pour les inconnues I, W, W22 et W2. Des valeurs de G (le gradient VAD isotrope), G2 (le terme azimutal ou anisotrope) et peuvent ensuite être calculées à partir des valeurs calculées pour W, W22 et W2. En conséquence, il est démontré que l'équation 1 est linéaire en I, Gi, G2 et la direction p. On remarquera que, comme indiqué plus haut dans les équations 2 et 3, les gradients calculés G et G2 sont reliés aux propriétés physiques des roches-,-, _ P /\(v) et VT P
CALCUL DES VARIATIONS AZIMUTALES DE VITESSE (VAV)
Puisque le traitement décrit sur la figure 3 détermine le décalage temporel associé aux variations azimutales de vitesse dans des traces sismiques, cette information de décalage temporel déterminé peut être utilisce pour calculer rcellement les variations azimutales de vitesse. Des étapes destinées à calculer les variations azimutales de vitesse sont décrites sur la figure 5. A une étape 90, le temps de propagation total T pour chaque trace est calculé par addition des décalages temporels déterminés dans le traitement décrit à propos de la fgure 3 ayant donné une corrélation maximale pour chaque trace et du décalage temporel obtenu en résultat d'une correction d'obliquité conventionnelle à To, le temps de propagation à déport nul. On peut ensuite utiliser l'équation suivante pour obtenir les variations azimutales de vitesse: T T2+ X2 (Eq. 17) Vnmo ()) dans laquelle: T = temps de propagation total 3 0 To = temps de propagation aller-retour à déport nul X= déport Vnno() = vitesse à variations azimutales en fonction de l'azimut l, avec: 12 =Wcos:+ 2W:cossinf+W sin:) (Eq. l 8) Vnm (/) En conséquence, on peut écrire le temps de propagation total sous la forme: T:=To2+[Wcos:+2W2cossinf+ W2:sin:flX: (Eq. l9) A une étape 92, l'équation l 9 peut être résolue en utilisant un procédé des moindres carrés linéaire connu des personnes ayant des compétences ordinaires dans l'art, en utilisant les décalages temporels extraits par le traitement d'inter-corrélation décrit à propos de la figure 3 ayant donné une corrélation maximale. Les équations 9, l0 et 11 peuvent ensuite étre utilisées pour obtenir G, G2 et p. La plus grande vitesse et la plus l0 petite vitesse sont calculées à partir des valeurs calculées de G et G2. La plus grande vitesse est donnée par Vfas =;, la plus petite vitesse est donnce par Esow =:, et l'azimut de la plus petite vitesse est donné par p. Comme les temps de propagation sont ajustés par la résolution par moindres carrés de l'équation l9, l'azimut qui est calculé en utilisant l'équation l9 est l'azimut du plus lS long temps de propagation. Si le temps de propagation est plus long, la vitesse est donc plus petite. Les variations azimutales de vitesse peuvent ensuite étre calculées à une étape 94 à partir de la relation suivante: 12 = 2COS2(- - 12 sin:(- (Eq. 20) Vnmo (()) Vslow Vfast (()) La valeur du décalage temporel du à une anisotropie azimutale peut étre déterminée pour chaque phénomène de réilexion en fonction de l'angle azimutal, et le décalage temporel approprié peut être appliqué à chaque trace pour compenser le décalage temporel azimutal. Une fois que les décalages temporels ont été appliqués aux traces de données sismiques du regroupement, une analyse VAD ainsi qu'une analyse VAA, telles que décrites dans la présente à propos de la figure 4, peuvent étre effectuées sur les traces ajustées.
CALCUL D'ERREURS ASSOCIEES AU CALCUL DES VARIATIONS DE
DECALAGE TEMPOREL, DES VARIATIONS DE VITESSE ET DES VARIATIONS
D'AMPLITUDE AVEC L'AZIMUT
Des erreurs associées au calcul des variations du décalage temporel avec l'azimut, des variations de vitesse avec l'azimut et des variations d'amplitude avec l'azimut peuvent étre estimées en utilisant une approche par moindres carrés. A une étape l 00, l'approche par moindre s carrés à l' estimation des erreurs ass ocices au calcul des variations du décalage temporel avec l'azimut est formulée en notation matricielle, et peut s'écrire: Ax=b (Eq.21) - dans laquelle b est une matrice lxN (c'est-àdire un vecteur) contenant les données (par S exemple des temps de propagation ou des amplitudes), A est une matrice MxN de coefficients (par exemple sin20) et x est une matrice lxM (c'est-à-dire un vecteur) des paramètres à résoudre. Par exemple, pour l'équation 19: 1 X, cos2, X, cosI sin>, X, sin, 1 X22 COS2 j2 X22 COS2 sin j2 X22 si02)2 A = 1 X32 cos2 f3 X32 cosf3 sinf3 X32 sin2 f3 1 Xn COS 4)n Xn cosqJn sint)n Xn sin ó7n To2 x = W,l W12 wl3 et Tl2 T22 b = T3.2 Tn2 dans lesquelles Tn est le temps de propagation observé pour des données ayant un déport source- récepteur Xn et un azimot 4;\n L'équation matricielle est équivalente à N équations simultanées à M inconnues. Dans l'exemple donné, M est égal à quatre. Pour une forrnulation par moindres carrés, N. le nombre de points de données, doit étre supérieur à M, le nombre d'inconnues. I1 existe divers procédés numériques conventionnels, connus
des personnes ayant des compétences ordinaires dans l'art, pour calculer les M inconnues.
A une étape 102, l'écart type des inconnues est calculé par la racine carrce des éléments diagonaux de la matrice E: E = (ATA)-'a2 (Eq. 22) dans laquelle l' expo s ant T désigne la transposce de la matrice, o2 est la variance (c' est-à dire la somme des carrés des différences entre les données et l'ajustement calculé, divisée par n-4). Les éléments diagonaux de la matrice E représentent les erreurs dans chaque inconnue, de sorte que la racine carrée du mième élément diagonal (c'est-à-dire,/<) est l'écart type de la mième inconnue. Le procédé des moindres carrés permet le calcul d'erreurs contenues dans la makice A, qui incluent à la fois la variance (erreur) due à une mauvaise qualité des données ou à du bruit aléatoire ainsi que l'erreur attendue due à la distribution des données. On remarquera que les éléments de la matrice A sont des fonctions des déports et
des azimuts des données, qui sont ensuite combinces et utilisoes pour calculer les erreurs.
De plus, comme la variance due à une mauvaise qualité des données ou à du bruit aléatoire peut être calculée indépendamment, la variance due à la mauvaise qualité des données ou à du bruit aléatoire: E = (ATA)-' (Eq. 23) et l'erreur due à la distribution des données: = a (Eq. 24) peuvent être séparées à une étape 104, et l'une ou l'autre ou les deux peuvent être comparées à une étape 106 au résultat VAV ou VAA pour confirmer ou infirmer le fait que le résultat porte une "empreinte" d'acquisition-un motif dû à la géométrie d'acquisition. Cette comparaison est typiquement effectuée de manière visuelle, bien que les personnes ayant des compétences ordinaires dans l'art puissent savoir effectuer la comparaison d'une manière mathématique. Un indice de précision des résultats obtenus est l' absence de la géométrie d'acquisition (telle que pli, déport maximal et déport minimal)
dans les résultats obtenus.
2 5 Des erreurs associées au calcul des variations d' amplitude avec l'azimut peuvent être calculées d'une manière analogue au calcul des erreurs associées au calcul des variations de décalage temporel avec l'azimut, mais en appliquant le traitement à l'équation 12 plutôt qu'à l'équation 19. Comme les variations de vitesse sont calculées à partir des variations du temps de propagation, les erreurs calculées pour les variations de décalage
temporel avec l'azimut sont applicables aux variations de vitesse avec l'azimut.
CALCUL D'UNE NOWELLE SOLUTION STATIQUE COMPATIBLE EN SURFACE
Des variations azimutales de vitesse (VAV) vont faire appara^tre des corrections statiques qui dépendent du déport et qui dépendent du temps dans des regroupements triés par milieu commun. Des solutions statiques compatibles en surface pour un regroupement par milieu commun utilisent typiquement une trace pilote comprenant une somme incluant tous les déports et azimuts dans le regroupement. S'il existe une VAV, les grands déports ne vont pas être alignés avec les petits déports et ne s'empileront pas de manière cohérente, par conséquent la trace pilote sera représentative des traces à petit déport, qui sont relativement peu affectées par la VAV. Pour obtenir une seule correction statique par trace, une inter-corrélation entre chaque trace et la trace pilote est recherchée sur une large
fenêtre temporelle (par exemple de deux secondes), typiquement autour d'un horizon cible.
Si le retard dû à la VAV varie avec le temps, alors la trace pilote ne sera pas représentative des traces à grands déports, puisque les traces à grands déports seront étirées et comprimoes par rapport aux petits déports et donc aussi par rapport à la trace pilote. Par exemple, un azimut particulier peut être dans la direction lente (entraînant un décalage vers des instants postérieurs) à un premier instant du kaj et, mais changer vers la direction
rapide (entrâînant un décalage vers des instants antérieurs) à un instant postérieur du kajet.
Par conséquent, lorsque cette trace est corrélée à la trace pilote, une mauvaise inter-
corrélation peut être obtenue, entraînant éventuellement une correction statique incorrectement extraite pour cette kace. Même si la correction statique extraite représente le retard "moyen" dû à la VAV, elle appara^itra comme du bruit dans le calcul de correction statique compatible en surface ("CSCS") car ce calcul suppose que les corrections statiques sont compatibles en surface, alors qu'en fait les retards dus à la VAV
ne le sont pas.
Selon une forme de réalisation de la présente invention, telle que décrite sur la figure 7, un processus itératif est utilisé dans lequel, à une étape 110, des décalages temporels sont calculés comme décrit plus haut à propos de la figure 5, qui sont ensuite appliqués aux traces de données sismiques d'une manière équivalente à une correction d'obliquité, suivie par une étape 112 à laquelle des calculs de corrections statiques compatibles en surface connus dans l'art antérieur sont effectués. L'étape 110 de calcul de décalages temporels est ensuite répétée. Ce processus consistant à calculer et appliquer les décalages temporels qui résultent des variations azimutales de vitesse et à calculer les corrections statiques compatibles en surface est répété jusqu'à ce que le processus converge. Cette détermination de la convergence du processus est effectuée à une étape 114. Un critère pouvant étre utilisé pour décider si le processus a convergé ou non est que les décalages temporels calculés à partir du calcul de corrections statiques compatibles en
surface soient globalement inférieurs à deux millisecondes.
Une fois que les décalages temporels dus à une anisotropie azimutale ont été appliqués aux traces de données sismiques du regroupement, une analyse VAD ainsi qu'une analyse VAA, telles que décrites dans la présente à propos de la figure 4, peuvent
être effectuces sur les traces ajustées.
La manière la plus commode de mettre en _uvre le procédé de l'invention décrite dans la présente consiste à écrire un programme informatique pour exécuter les étapes décrites dans la présente sur une station de travail ou autre système informatique numérique conventionnel d'un type normalement utilisé dans l'industrie. Les personnes ayant des compétence s ordinaires dans l' art pourront créer un tel programme en se basant sur les opérations décrites dans la présente. La figure 9 représente un tel système informatique conventionnel comprenant une unité centrale 122, un dispositif d'affichage 124, un dispositif d'entrée 126 et un traceur 128. Le programme d'ordinateur destiné à mettre en _uvre l' invention sera normalement résident sur un support de mémori sation (non représenté) associé à l'unité centrale. Ce programme d'ordinateur peut être transporté sur un CD-ROM ou autre support de mémorisation représenté symboliquement par un
support de mémorisation 130.
Les résultats des calculs selon cette invention peuvent étre présentés par un logiciel de visualisation disponible dans le commerce. Un tel logiciel est bien connu des personnes ayant des compétences ordinaires dans l'art et ne sera pas plus décrit dans la présente. On comprendra que les résultats des procédés de l'invention peuvent être affichés, tracés ou
les deux.
Claims (26)
1. Procédé de traitement de donnces sismiques pour estimer un décalage temporel du à une anisotropie de vitesse dans le sous-sol terreske, caractérisé en ce qu'il comprend les étapes consistant à: former un regroupement de kaces de données sismiques; et inter-corréler des traces de donnces sismiques sélectionnées incluses dans ledit regroupement à l'intérieur de fenétres temporelles sélectionnées afin d'estimer le décalage temporel dans les traces de donnces sismiques incluses dans ledit regroupement dû à une
anisotropie de vites se dans le sous- sol terrestre.
l0
2. Procédé selon la revendication l, caractérisé en ce qu'il comprend en ouke l'étape consistant à ajuster des kaces de données sismiques incluses dans ledit regroupement par la valeur du décalage temporel estimé dans ladite trace de données sismiques ajustée du à
une anisotropie de vitesse.
3. Procédé selon la revendication 2, caractérisé en ce qu'il comprend en ouke l'étape l5 consistant à effectuer une analyse des variations d'amplitude avec l'angle d'incidence sur
lesdites traces de données sismiques ajustées.
4. Procédé selon la revendication 2, caractérisé en ce qu'il comprend en outre l'étape consistant à effectuer une analyse des variations d'amplitude avec l'azimut sur lesdites
kaces de donnces sismiques ajustées.
5. Procédé selon la revendication 2, caractérisé en ce qu'il comprend en outre les étapes consistant à: déterminer l'angle d'incidence pour des kaces de données sismiques incluses dans ledit regroupement; et appliquer un procédé d'analyse par moindres carrés aux valeurs du coefficient de réflexion, de l'angle d'azimut source-récepteur et de l'angle d'incidence pour des traces de données sismiques incluses dans ledit regroupement, afin de calculer les variations d' amplitude avec l' azimut et les variations d'amplitude avec le déport dans des kaces de
données sismiques incluses dans ledit regroupement.
6. Procédé selon la revendication 2, caractérisé en ce qu'il comprend en oube l'étape consistant à appliquer un procédé d'analyse par moindres carrés aux décalages temporels desdites traces de données sismiques afin de calculer les variations de vitesse avec l'azimut
dans des traces de données sismiques incluses dans ledit regroupement.
7. Procédé de traitement de données sismiques, caractérisé en ce qu'il comprend les étapes consistant à: former un regroupement de traces de données sismiques; effectuer un calcul de corrections statiques compatibles en surface sur des kaces de données sismiques incluses dans ledit regroupement; inter-corréler des traces de données sismiques successivement sélectionnées dans ledit regroupement afin d'estimer le décalage temporel dans des traces de données sismiques incluses dans ledit regroupement du à une anisotropie azimutale de vitesse dans le sous- sol terrestre, en appliquant un procédé d'analyse par moindres carrés aux décalages temporels desdites traces de données sismiques afin de calculer les variations de vitesse avec l'azimut dans des traces de donnces sismiques incluses dans ledit regroupement; utiliser les variations calculées de vitesse avec l'azimut pour calculer des décalages temporels dans des traces de données sismiques incluses dans ledit regroupement; appliquer lesdits décalages temporels calculés auxdites traces de données sismiques incluses dans ledit regroupement; et appliquer un procédé d'analyse par moindres carrés aux valeurs du coefficient de réflexion, de l'angle d'azimut source-récepteur et de l'angle d' incidence pour des traces de données sismiques incluses dans ledit regroupement afin de calculer les variations d'amplitude avec l' azimut et les variations d' amplitude avec le déport dans de s traces de
donnces sismiques incluses dans ledit regroupement.
8. Procédé de traitement de données sismiques pour estimer un décalage temporel dû à une anisotropie de vitesse dans le sous-sol terrestre, caractérisé en ce qu'il comprend les étapes consistant à: former un regroupement de traces de données sismiques; former une kace pilote par combinaison d'une pluralité sélectionnée desdites traces de données sismiques à l'intérieur d'une fenétre temporelle sélectionnée; et intercorréler une trace de données sismiques sélectionnée incluse dans ledit regroupement avec ladite trace pilote afin d'estimer le décalage temporel dans ladite trace de données sismiques sélectionnée du à une anisotropie de vitesse dans le sous-sol terrestre.
9. Procédé de traitement de données sismiques pour estimer un décalage temporel du à une anisotropie de vitesse dans le sous-sol terrestre, caractérisé en ce qu'il comprend les étapes consistant à: (a) former un regroupement de traces de données sismiques; (b) former une trace pilote par combinai son d'une pluralité sélectionnce de sdites traces de données sismiques à l'intérieur d'une fenétre temporelle sélectionnée; (c) inter- corréler une trace de données si smiques sélectionnce incluse dans le dit regroupement avec ladite trace pilote afn d'estimer le décalage temporel dans ladite trace de données sismiques sélectionnée dû à une anisotropie de vitesse dans le sous-sol terrestre; et répéter les étapes (b) et (c) jusqu'à ce que toutes les traces de données sismiques
incluses dans ledit regroupement aient été inter-corrélées avec une trace pilote.
10. Procédé de traitement de données sismiques, caractérisé par les étapes consistant à: (a) former un regroupement de traces de données sismiques; (b) former une trace pilote par combinaison d'une pluralité sélectionnée desdites traces de données sismiques à l'intérieur d'une fenêtre temporelle sélectionnée; (c) inter-corréler une trace de données sismiques sélectionnée incluse dans ledit regroup ement ave c ladite trace p il ote afin d' e stimer l e dé calage temp or el dans ladite trac e de données sismiques sélectionnée dû à une anisotropie de vitesse dans le sous-sol terrestre; répéter les étapes (b) et (c) jusqu'à ce que toutes les traces de données sismiques dans ledit regroupement aient été inter-corrélées avec une trace pilote, et ajuster chaque trace de données sismiques sélectionnée par la valeur du décalage temporel estimé dans chaque trace de données sismiques sélectionnée dû à une anisotropie
de vitesse.
11. Procédé selon la revendication 10, caractérisé en ce qu'il comprend en outre l'étape consistant à effectuer une analyse des variations d'amplitude avec l'angle d'incidence sur
lesdites traces de données sismiques aJustées.
12. Procédé selon la revendication 10, caractérisé en ce qu'il comprend en outre l'étape consistant à effectuer une analyse des variations d'amplitude avec l'azimut sur lesdites
traces de donnces sismiques adustées.
13. Procédé de traitement de données sismiques, caractérisé en ce qu'il comprend les étapes consistant à: (a) former un regroupement de traces de données sismiques; (b) former une trace pilote par combinaison d'une pluralité sélectionnée desdites traces de données sismiques à l'intérieur d'une fenêtre temporelle sélectionnée; (c) inter-corréler une trace de données sismiques sélectionnce incluse dans ledit regroupement avec ladite trace pilote afin d'estimer le décalage temporel dans ladite trace de donnces sismiques sélectionnée dû à une anisotropie de vitesse dans le sous-sol terrestre; répéter les étapes (b) et (c) jusqu'à ce que toutes les traces de données sismiques dans ledit regroupement aient été intercorrélées avec une trace pilote; S ajuster chaque kace de données sismiques sélectionnée par la valeur du décalage temporel estimé dans chaque trace de données sismiques sélectionnce du à une anisotropie de vitesse; déterminer l' angle d' incidence pour chaque trac e de do nnée s si smiques sélectionnée,; et appliquer un procédé d'analyse par moindres carrés aux valeurs du coefficient de réflexion, de l' angle d' azimut source-récepteur et de l' angle d' incidence des dites traces de donnces sismiques afin de calculer les variations d'amplitude avec l'azimut et les variations d' amplitude avec le déport dans des traces de données sismiques incluses dans
ledit regroupement.
14. Procédé de traitement de données sismiques, caractérisé en ce qu'il comprend les étapes consistant à: (a) former un regroupement de traces de données sismiques; (b) former une trace pilote par combinaison d'une pluralité sélectionnée desdites traces de données sismiques à l'intérieur d'une fenêtre temporelle sélectionnée; (c) inter-corréler une trace de données sismiques sélectionnée incluse dans ledit regroupement avec ladite trace pilote afin d'estimer le décalage temporel dans ladite trace de données sismiques sélectionnée dû à une anisotropie de vitesse dans le sous-sol terrestre; répéter les étapes (b) et (c) jusqu'à ce que toutes les traces de donnces sismiques dans ledit regroupement aient été intercorrélées avec une trace pilote; et appliquer un procédé d'analyse par moindres carrés aux décalages temporels desdites traces de données sismiques afin de calculer les variations de vitesse avec l'azimut
dans des traces de données sismiques incluses dans ledit regroupement.
15. Procédé selon la revendication 13, caractérisé en ce qu'il comprend en outre l'étape consistant à utiliser une analyse par moindres carrés pour estimer des erreurs associées au
calcul des variations d'amplitude dans lesdites trace, s de données sismiques sélectionnées.
16. Procédé selon la revendication 14, caractérisé en ce qu'il comprend en outre l'étape consistant à utiliser une analyse par moindres carrés pour estimer des erreurs associces au calcul des variations de décalage temporel dans lesdites traces de données sismiques sélectionnées.
17. Procédé selon la revendication 14, caractérisé en ce qu'il comprend en oube l'étape consistant à utiliser une analyse par moindres carrés pour estimer des erreurs associées au calcul des variations de vitesse dans lesdites traces de données sismiques sélectionnées.
18. Procédé de traitement de donnces sismiques, caractérisé en ce qu'il comprend les étapes consistant à: (a) former un regroupement de traces de donnces sismiques; (b) effectuer un calcul de corrections statiques compatibles en surface sur lesdites traces de données sismiques; (c) inter-corréler des traces de données sismiques successivement sélectionnées dans ledit regroupement avec une trace pilote formée par combinaison d'une pluralité sélectionnée de traces de données sismiques dudit regroupement à l'intérieur de fenétres temporelles sélectionnées, afin d'estimer le décalage temporel dans lesdites traces de données sismiques du à une anisotropie azimutale de vitesse dans le sous-sol terrestre, et appliquer un procédé d'analyse par moindres carrés aux décalages temporels desdites traces de données sismiques afin de calculer les variations de vitesse avec l'azimut dans des traces de donnces sismiques incluses dans ledit regroupement; (d) utiliser les variations calculées de vitesse avec l'azimut pour calculer des décalages temporels dans des kaces de données sismiques incluses dans ledit regroupement; (e) appliquer lesdits décalages temporels calculés auxdites traces de données sismiques incluses dans ledit regroupement; répéter les étapes (b)-(e); et appliquer un procédé d'analyse par moindres carrés aux valeurs du coefficient de réflexion, de l'angle d'azimut source-récepteur et de l'angle d' incidence des dites traces de données sismiques afin de calculer les variations d'amplitude avec l'azimut et les variations d' amplitude avec le déport dans des traces de donnces sismique s incluses dans
ledit regroupement.
19. Procédé de traitement de données sismiques, caractérisé en ce qu'il comprend les étapes consistant à: (a) former un regroupement de traces de données sismiques, (b) sélectionner une fenétre temporelle dudit regroupement; (c) sélectionner une pluralité desdites traces, (d) sommer ladite pluralité de traces à l'intérieur de ladite fenêtre temporelle; (e) sélectionner une trace dans ledit regroupement et inter-corréler ladite trace sélectionnée avec ladite somme de pluralité de traces; (f) déterminer un déplacement temporel de ladite kace sélectionnée qui donne une corrélation sensiblement maximale de ladite trace sélectionnée avec ladite somme de pluralité de traces à l'intérieur de ladite fenêhe temporelle sélectionnce; (g) répéter les étapes (c)-(f) pour des traces successives et des sommes de pluralités sélectionnées de traces jusqu'à ce qu'un déplacement temporel ait été déterminé pour chaque trace dudit regroupement qui donne une corrélation sensiblement maximale de chaque trace dudit regroupement avec une somme de pluralité de traces à l'intérieur de ladite fenêtre temporelle; et (h) répéter les étapes (b)(g) pour des fenêtres temporelles successives à l'intérieur
dudit regroupement.
20. Procédé de traitement de données sismiques, caractérisé en ce qu'il comprend les étapes consistant à: (a) former un regroupement de traces de donnces sismiques; (b) sélectionner une fenêtre temporelle dudit regroupement; (c) sélectionner une pluralité desdites kaces; (d) sommer ladite pluralité de kaces à l'intérieur de ladite fenêtre temporelle; (e) sélectionner une trace dudit regroupement et inter-corréler ladite trace sélectionnce avec ladite somme de pluralité de traces; (f) déterminer un déplacement temporel de ladite trace sélectionnée qui donne une corrélation sensiblement maximale de ladite trace sélectionnée avec ladite somme de pluralité de traces à l'intérieur de ladite fenêtre temporelle sélectionnce; (g) répéter les étapes (c)-(f) pour des traces successives et des sommes de pluralités sélectionnces de traces jusqu'à ce qu'un déplacement temporel ait été déterminé pour chaque trace dudit regroupement qui donne une corrélation sensiblement maximale de chaque trace dudit regroupement avec une somme de pluralité de traces à l'intérieur de ladite fenêtre temporelle; (h) répéter les étapes (b)-(g) pour des fenêtres temporelles successives à l'intérieur dudit regroupement; (i) appliquer le déplacement temporel déterminé pour chaque trace à l'intérieur de ladite fenêtre temporelle à la trace située au milieu de chaque fenêhe temporelle; et (j) interpoler des décalages temporels pour la partie restante de chaque trace de
donnces enke les milieux des fenêtres.
21. Procédé selon la revendication 20, caractérisé en ce qu'il comprend en outre les étapes consistant à: S déterminer le coeffcient de réflexion et l'angle d'azimut pour chaque trace du regroupement; calculer la valeur de l'angle d'incidence pour chaque trace du regroupement, et utiliser un procédé par moindres carrés pour calculer les coefficients de réflexion en fonction de l' angle d' azimut et de l' ang le d'incidence pour les trace s si smique s du
1 0 regroupement.
22. Procédé selon la revendication 20, caractérisé en ce qu'il comprend en outre l'étape consistant à utiliser un procédé linéaire par moindres carrés pour calculer le demi petit axe et le demi grand axe, ainsi que l'angle entre une direction prédéfinie d'azimut zéro et la
direction du demi grand axe, d'une ellipse représentant les variations azimutales de vitesse.
23. Procédé selon la revendication 22, caractérisé en ce qu'il comprend en outre les étapes consistant à: formuler une matrice pour une analyse par moindres carrés destinée à estimer les erreurs associées au calcul des variations du décalage temporel avec l'azimut, calculer une estimation d'erreur par la racine carrée des éléments diagonaux de ladite matrice, et comparer l'erreur estimée au calcul des variations azimutales de vitesse afin de confirmer ou d'infirmer le fait que les variations azimutales de vitesse portent une
empreinte d'acquisition.
24. Ordinateur numérique, caractérisé par le fait d'êke programmé pour utiliser des traces de données sismiques obtenues sur une région du soussol terrestre afin d'exécuter un procédé comprenant les étapes consistant à: former un regroupement de traces de donnces sismiques; former une trace pilote par combinaison d'une pluralité sélectionnée desdites kaces de données sismiques à l'intérieur d'une fenêke temporelle sélectionnée, et inter-corréler une trace de données sismiques sélectionnée incluse dans ledit regroupement avec ladite trace pilote afin d'estimer le décalage temporel dans ladite trace de données sismiques sélectionnée dû à une anisokopie de vitesse dans le sous-sol terreske.
25. Dispositif lisible par un ordinateur numérique, caractérisé en ce qu'il contient des instructions définissant le procédé suivant et des inskuctions destinées à l'ordinateur pour exécuter ledit procédé: former un regroupement de traces de donnces sismiques; former une trace pilote par combinaison d'une pluralité sélectionnce desdites traces de données sismiques à l'intérieur d'une fenétre temporelle sélectionnée; et inter-corréler une trace de donnces sismiques sélectionnée incluse dans ledit regroupement avec ladite trace pilote afin d'estimer le décalage temporel dans ladite trace de données sismiques sélectionnée du à une anisotropie de vitesse dans le sous-sol terrestre.
26. Dispositif selon la revendication 25, caractérisé en ce que ledit dispositif est choisi dans le groupe constitué par une bande magnétique, un disque
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US09/855,925 US6681184B2 (en) | 2001-05-15 | 2001-05-15 | System for estimating azimuthal variations in seismic data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| FR2824917A1 true FR2824917A1 (fr) | 2002-11-22 |
| FR2824917B1 FR2824917B1 (fr) | 2011-06-24 |
Family
ID=25322447
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR0205945A Expired - Lifetime FR2824917B1 (fr) | 2001-05-15 | 2002-05-15 | Systeme pour l'estimation de variations azimutales dans des donnees sismiques |
Country Status (4)
| Country | Link |
|---|---|
| US (2) | US6681184B2 (fr) |
| CA (1) | CA2386568C (fr) |
| FR (1) | FR2824917B1 (fr) |
| GB (1) | GB2377495B (fr) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104483701A (zh) * | 2014-12-08 | 2015-04-01 | 中国石油天然气集团公司 | 基于全方位共成像点道集的各向异性表征方法以及系统 |
Families Citing this family (43)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6760667B1 (en) | 2001-09-14 | 2004-07-06 | Emerald Geoscience Research Corporation | Method for generating P-S and S-S seismic data and attributes from P-P seismic data |
| US6817218B1 (en) * | 2001-09-14 | 2004-11-16 | Emerald Geoscience Research Corp. | Non-linear statistically constrained AVO inversion |
| GB2383414B (en) * | 2001-12-22 | 2005-07-13 | Westerngeco Ltd | A method of and apparatus for processing seismic data |
| GB2387454A (en) * | 2002-04-10 | 2003-10-15 | Westerngeco Ltd | Processing Seismic Data |
| US6904368B2 (en) * | 2002-11-12 | 2005-06-07 | Landmark Graphics Corporation | Seismic analysis using post-imaging seismic anisotropy corrections |
| FR2854247B1 (fr) * | 2003-04-28 | 2005-07-29 | Geophysique Cie Gle | Procede de traitement de cubes sismiques correspondant pour une meme zone au sol, a differentes valeurs de deports source/recepteur et/ou d'angles d'incidence |
| FR2858064B1 (fr) * | 2003-07-21 | 2005-09-16 | Geophysique Cie Gle | Procede de pointe bispectral des parametres de correction d'obliquite anelliptique |
| US20050033518A1 (en) * | 2003-08-07 | 2005-02-10 | Jenner Edward Louis | Method for wavelet-based seismic amplitude inversion |
| CA2586603C (fr) * | 2004-11-18 | 2014-04-01 | Exxonmobil Upstream Research Company | Procede destine a quantifier et a attenuer d'une variation d'amplitude avec decalage (avo) azimut induite par inclinaison |
| FR2884620B1 (fr) * | 2005-04-19 | 2007-07-13 | Cie Generale De Geophysique Sa | Procede de traitement de donnees sismiques en vue d'une caracterisation avo ou avoa |
| US7674003B2 (en) | 2006-04-20 | 2010-03-09 | Streamlight, Inc. | Flashlight having plural switches and a controller |
| WO2008043082A2 (fr) | 2006-10-05 | 2008-04-10 | Splunk Inc. | Moteur de recherche sur séries chronologiques |
| US7460437B2 (en) * | 2007-01-03 | 2008-12-02 | Weinman Geoscience | Seismic data processing method and system for migration of seismic signals incorporating azimuthal variations in the velocity |
| US20090052279A1 (en) * | 2007-08-21 | 2009-02-26 | Charles Sicking | Method and system for determining azimuth parameters for seismic data processing |
| US8548782B2 (en) | 2007-08-24 | 2013-10-01 | Exxonmobil Upstream Research Company | Method for modeling deformation in subsurface strata |
| US8768672B2 (en) * | 2007-08-24 | 2014-07-01 | ExxonMobil. Upstream Research Company | Method for predicting time-lapse seismic timeshifts by computer simulation |
| RU2488146C2 (ru) * | 2008-03-27 | 2013-07-20 | Предрилл Стрессиз Интернэшнл Пти. Лтд. | Прогнозирование напряжений на изучаемой площади |
| GB2460013B (en) * | 2008-03-31 | 2010-10-13 | Statoilhydro Asa | A method for reflection time shift matching a first and a second set of seismic reflection data |
| US8174927B2 (en) | 2008-12-17 | 2012-05-08 | Westerngeco L.L.C. | Method for optimizing acoustic source array performance |
| US8947974B2 (en) * | 2009-06-23 | 2015-02-03 | Baker Hughes Incorporated | Seismic measurements while drilling |
| US8902708B2 (en) * | 2009-11-05 | 2014-12-02 | Chevron U.S.A. Inc. | System and method for seismic beam formation that accounts for equipment misalignment |
| CN102455436B (zh) * | 2010-11-02 | 2013-10-30 | 中国石油大学(北京) | 有限方位角条件下压噪叠前纵波各向异性裂缝检测方法 |
| US10310119B2 (en) * | 2011-06-24 | 2019-06-04 | Ion Geophysical Corporation | Method and apparatus for seismic noise reduction |
| US9784863B2 (en) | 2011-09-26 | 2017-10-10 | Exxonmobil Upstream Research Company | S-wave anisotropy estimate by automated image registration |
| CN102879821B (zh) * | 2012-09-26 | 2015-01-21 | 中国石油天然气股份有限公司 | 一种针对地震叠前道集的同相轴精细拉平处理方法 |
| FR3001301B1 (fr) * | 2013-01-24 | 2015-08-07 | Cggveritas Services Sa | Appareil et procede pour determiner la signature de champ lointain pour une source sismique vibratoire marine |
| US10346357B2 (en) | 2013-04-30 | 2019-07-09 | Splunk Inc. | Processing of performance data and structure data from an information technology environment |
| US10225136B2 (en) | 2013-04-30 | 2019-03-05 | Splunk Inc. | Processing of log data and performance data obtained via an application programming interface (API) |
| US10997191B2 (en) | 2013-04-30 | 2021-05-04 | Splunk Inc. | Query-triggered processing of performance data and log data from an information technology environment |
| US10019496B2 (en) | 2013-04-30 | 2018-07-10 | Splunk Inc. | Processing of performance data and log data from an information technology environment by using diverse data stores |
| US10318541B2 (en) | 2013-04-30 | 2019-06-11 | Splunk Inc. | Correlating log data with performance measurements having a specified relationship to a threshold value |
| US10353957B2 (en) | 2013-04-30 | 2019-07-16 | Splunk Inc. | Processing of performance data and raw log data from an information technology environment |
| US10614132B2 (en) | 2013-04-30 | 2020-04-07 | Splunk Inc. | GUI-triggered processing of performance data and log data from an information technology environment |
| US10386516B2 (en) * | 2013-05-15 | 2019-08-20 | Conocophillips Company | Time-lapse 4D scattering for imaging hydraulically induced fractures |
| US20150226868A1 (en) * | 2014-02-10 | 2015-08-13 | Microseismic, Inc. | Method for microseismic event moment magnitude estimation |
| EP3143249B1 (fr) * | 2014-07-17 | 2023-03-08 | Halliburton Energy Services, Inc. | Élimination de bruit pour des données de détection acoustique répartie |
| CA2989619C (fr) * | 2015-06-17 | 2023-05-16 | Conocophillips Company | Estimation de gradient d'azimut sismique |
| US10067255B2 (en) | 2015-09-04 | 2018-09-04 | Saudi Arabian Oil Company | Automatic quality control of seismic travel time |
| US10386519B2 (en) | 2015-12-18 | 2019-08-20 | Saudi Arabian Oil Company | Automated near surface analysis by surface-consistent refraction methods |
| CN107797147B (zh) * | 2017-10-25 | 2019-07-05 | 中国石油化工股份有限公司 | 一种地震初至异常值快速剔除方法 |
| US10935681B2 (en) | 2018-05-31 | 2021-03-02 | Cgg Services Sas | Methods and systems performing joint VVAz and AVAz inversion |
| CN109143352B (zh) * | 2018-07-04 | 2019-10-01 | 中国石油大学(华东) | 一种各向异性介质地震反射特征方程建立方法 |
| CN110135654B (zh) * | 2019-05-24 | 2021-08-27 | 北京百度网讯科技有限公司 | 用于预测强对流天气的方法和装置 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4613960A (en) * | 1983-09-08 | 1986-09-23 | Mobil Oil Corporation | Method for constructing and using an optimal pilot trace in surface consistent statics estimations |
| GB2254146A (en) * | 1991-03-08 | 1992-09-30 | Geco As | Determining residual static shift in seismic data |
| US6061301A (en) * | 1997-06-30 | 2000-05-09 | Atlantic Richfield Company | Filtering of overburden azimuthal anisotropy effects from 3D seismic survey signals |
Family Cites Families (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE1492425A1 (de) * | 1963-01-31 | 1969-12-04 | Pfrimmer & Co J | Verfahren zum Sterilisieren von chirurgischen Naehmitteln |
| US3417370A (en) * | 1966-12-23 | 1968-12-17 | Texas Instruments Inc | Seismic velocity correlation |
| US3611278A (en) * | 1969-07-17 | 1971-10-05 | Mobil Oil Corp | Interval velocity determination |
| US4206509A (en) * | 1978-03-03 | 1980-06-03 | Mobil Oil Corporation | Method of enhancing seismic reflection signals for nonsurface-consistent static time shifts |
| US4570246A (en) * | 1982-02-01 | 1986-02-11 | Chevron Research Company | Method for the interpretation of statistically-related seismic records to yield valuable characteristics, such as gas-bearing potential and lithology of strata |
| US4779237A (en) * | 1984-08-27 | 1988-10-18 | Amoco Corporation | Method of geophysical exploration including processing and displaying seismic data to obtain a measure of subterranean formation rock properties |
| US4755972A (en) * | 1985-09-26 | 1988-07-05 | Amoco Corporation | Method of seismic exploration including processing and displaying shear wave seismic data |
| US4970696A (en) | 1988-07-13 | 1990-11-13 | Atlantic Richfield Company | Method for conducting three-dimensional subsurface and marine seismic surveys |
| US5258960A (en) | 1992-11-27 | 1993-11-02 | Atlantic Richfield Company | Seismic velocity estimation method |
| GB2278920B (en) | 1993-06-07 | 1996-10-30 | Geco As | Method of determining earth elastic parameters in anistropic media |
| US5532978A (en) | 1993-12-22 | 1996-07-02 | Atlantic Richfield Company | Method of correction of seismic data for the effects of azimuthal anisotropy |
| US5508973A (en) | 1994-06-06 | 1996-04-16 | Western Atlas International, Inc. | Method for determining the principal axes of azimuthal anisotropy from seismic P-wave data |
| US5610875A (en) * | 1995-09-29 | 1997-03-11 | Western Atlas International, Inc. | 3-D converted shear wave rotation with layer stripping |
| US5764516A (en) * | 1995-12-29 | 1998-06-09 | Atlantic Richfield Company | Method and system for surface-consistent phase and time lag correction of seismic data |
| US5933789A (en) * | 1996-11-12 | 1999-08-03 | Atlantic Richfield Company | Method and system for applying dispersive normal moveout corrections to seismic survey signals |
| US6041018A (en) * | 1997-11-13 | 2000-03-21 | Colorado School Of Mines | Method for correcting amplitude and phase differences between time-lapse seismic surveys |
| US6263284B1 (en) * | 1999-04-22 | 2001-07-17 | Bp Corporation North America Inc. | Selection of seismic modes through amplitude characteristics |
-
2001
- 2001-05-15 US US09/855,925 patent/US6681184B2/en not_active Expired - Lifetime
-
2002
- 2002-05-14 GB GB0211026A patent/GB2377495B/en not_active Expired - Lifetime
- 2002-05-15 CA CA2386568A patent/CA2386568C/fr not_active Expired - Lifetime
- 2002-05-15 FR FR0205945A patent/FR2824917B1/fr not_active Expired - Lifetime
-
2003
- 2003-10-20 US US10/689,423 patent/US6898147B2/en not_active Expired - Lifetime
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4613960A (en) * | 1983-09-08 | 1986-09-23 | Mobil Oil Corporation | Method for constructing and using an optimal pilot trace in surface consistent statics estimations |
| GB2254146A (en) * | 1991-03-08 | 1992-09-30 | Geco As | Determining residual static shift in seismic data |
| US6061301A (en) * | 1997-06-30 | 2000-05-09 | Atlantic Richfield Company | Filtering of overburden azimuthal anisotropy effects from 3D seismic survey signals |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104483701A (zh) * | 2014-12-08 | 2015-04-01 | 中国石油天然气集团公司 | 基于全方位共成像点道集的各向异性表征方法以及系统 |
| CN104483701B (zh) * | 2014-12-08 | 2017-04-12 | 中国石油天然气集团公司 | 基于全方位共成像点道集的各向异性表征方法以及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| US20030018435A1 (en) | 2003-01-23 |
| US6681184B2 (en) | 2004-01-20 |
| US20040109387A1 (en) | 2004-06-10 |
| GB2377495A (en) | 2003-01-15 |
| CA2386568A1 (fr) | 2002-11-15 |
| CA2386568C (fr) | 2013-10-15 |
| US6898147B2 (en) | 2005-05-24 |
| GB2377495B (en) | 2005-07-13 |
| FR2824917B1 (fr) | 2011-06-24 |
| GB0211026D0 (en) | 2002-06-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| FR2824917A1 (fr) | Systeme pour l'estimation de variations azimutales dans des donnees sismiques | |
| EP0789849B1 (fr) | Methode pour representer dans un domaine spatio-temporel la trajectoire d'au moins un puits de forage | |
| US6501703B2 (en) | Static correction method for exploration seismic data using first arrivals of seismic waves | |
| EP2369369B1 (fr) | Procédé de surveillance d'un site de stockage géologique de gaz par inversion stratigraphique de données sismiques | |
| US10234582B2 (en) | Joint inversion of seismic data | |
| FR2708350A1 (fr) | Méthode d'analyse de signaux sismiques. | |
| US6061301A (en) | Filtering of overburden azimuthal anisotropy effects from 3D seismic survey signals | |
| FR3027688A1 (fr) | ||
| US5128899A (en) | Method for sharpening prestack depth migration images by removing residual moveout from common image point gathers before stacking | |
| AU2003201448A1 (en) | A method of an apparatus for processing seismic data and determining a calibration filter | |
| CA1229906A (fr) | Procede d'obtention et de traitement informatique de donnees sismiques relevees aupres d'un puits d'exploration | |
| Gimbert et al. | A multi‐physics experiment with a temporary dense seismic array on the Argentière glacier, French Alps: The RESOLVE project | |
| US6834235B2 (en) | Method of and apparatus for processing seismic data | |
| CA2456596C (fr) | Methode pour mesurer les similarites locales entre plusieurs cubes de traces sismiques | |
| FR2971859A1 (fr) | Procede et dispositif de lissage a preservation du temps de trajet | |
| Guo et al. | Becoming effective velocity-model builders and depth imagers, Part 2—The basics of velocity-model building, examples and discussions | |
| Lu et al. | Redatuming through a salt canopy and target-oriented salt-flank imaging | |
| AU2002317974A1 (en) | A method of processing geophysical data | |
| CA2590859C (fr) | Procede de correction dans le traitement de traces sismiques | |
| WO2014023910A2 (fr) | Procédé de restauration de données de puits de forage | |
| EP3140677B1 (fr) | Procédé de traitement d'images sismiques | |
| WO2008157029A2 (fr) | Détermination de contraste de densité au niveau d'une interface souterraine | |
| EP0412013B1 (fr) | Procédé de traitement de données de sismique réflexion pour l'obtention de coupes sismiques | |
| Vesnaver et al. | A workflow for processing mono‐channel Chirp and Boomer surveys | |
| WO2001031364A1 (fr) | Procede de prospection sismique mettant en oeuvre un traitement sur les ondes converties |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| TP | Transmission of property | ||
| CD | Change of name or company name |
Owner name: ION GEOPHYSICAL CORPORATION, US Effective date: 20111130 |
|
| PLFP | Fee payment |
Year of fee payment: 15 |
|
| PLFP | Fee payment |
Year of fee payment: 16 |
|
| PLFP | Fee payment |
Year of fee payment: 17 |
|
| PLFP | Fee payment |
Year of fee payment: 18 |
|
| PLFP | Fee payment |
Year of fee payment: 19 |
|
| PLFP | Fee payment |
Year of fee payment: 20 |