WO2003003052A1 - Perfectionnements aux techniques de pointe de surface - Google Patents

Perfectionnements aux techniques de pointe de surface Download PDF

Info

Publication number
WO2003003052A1
WO2003003052A1 PCT/FR2002/002278 FR0202278W WO03003052A1 WO 2003003052 A1 WO2003003052 A1 WO 2003003052A1 FR 0202278 W FR0202278 W FR 0202278W WO 03003052 A1 WO03003052 A1 WO 03003052A1
Authority
WO
WIPO (PCT)
Prior art keywords
dimensions
point
propagation
dimension
points
Prior art date
Application number
PCT/FR2002/002278
Other languages
English (en)
Inventor
Abdelkrim Talaalout
Dominique Doledec
Original Assignee
Compagnie Generale De Geophysique
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Compagnie Generale De Geophysique filed Critical Compagnie Generale De Geophysique
Priority to US10/482,583 priority Critical patent/US8014951B2/en
Priority to MXPA04000063A priority patent/MXPA04000063A/es
Priority to GB0330107A priority patent/GB2393515B/en
Priority to AU2002333947A priority patent/AU2002333947B2/en
Priority to CA2452206A priority patent/CA2452206C/fr
Publication of WO2003003052A1 publication Critical patent/WO2003003052A1/fr

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. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/301—Analysis for determining seismic cross-sections or geostructures
    • 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/32—Transforming one recording into another or one representation into another
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/67—Wave propagation modeling

Definitions

  • the present invention relates to improvements in surface pointing techniques used in seismic processing.
  • the geophysicist uses the techniques of "seismic reflection". These consist of emitting acoustic signals on the surface of the ground and recording them after their successive reflections at the limits of the geological layers. These layers constitute the basement which we seek to create a geometric model.
  • 3D acquisition This consists of distributing transmitters and receivers in a horizontal surface plane along a grid of coordinates X, Y.
  • Such a three-dimensional acquisition makes it possible to obtain a volume whose dimensions are two geometric axes X, Y and a time axis or acquisition depth Z. Besides these three spatial dimensions (X, Y and vertical axis Z), of course other quantities are of interest to the physicist, and in particular:
  • the type of seismic waves can simultaneously record compression waves (P waves) and shear waves (S waves); the relationship between the arrival times of these two types of waves is a direct link with the nature of the medium in which the waves propagate.
  • P waves compression waves
  • S waves shear waves
  • Time or depth migration consists in moving the points of the acquisition volume taking into account a speed field.
  • the speed used in this process can be disturbed and constitute an additional dimension of analysis.
  • the Applicant has proposed an approach consisting in regularly distributing the disturbance factors on a speed model, and being interested for a position X, Y, at the value of the sample associated with the same reflector in each of the migrated cubes.
  • Document WO 99/67660 describes a method for monitoring changes in the subsoil. This method consists in determining a displacement field making it possible to pass from a sampled volume of seismic data measured at a time T 0 to a field of seismic data measured at a later time T. The displacement field thus determined is applied to the first seismic volume so as to obtain an image of the subsoil.
  • the information that can be drawn from such a mode of interpretation is necessarily limited.
  • the invention proposes an effective and original contribution to the resolution of this problem.
  • a technique conventionally used for the analysis of 3D seismic volumes is that of the surface pointing, an example of which is given in patent FR 2 657 695.
  • This document relates to a method of surface pointing in a three-dimensional volume according to which a volume representing the values of a parameter sampled in a two-dimensional image using a Peano-Hilbert curve and an automatic pointing program is applied to the 2D image obtained.
  • a difficulty is that the nature and shape of the reflector to be interpreted can change from one volume to another along the dimension that constitutes said quantity.
  • mapping sample The consistency between the different 3D volumes is not necessarily included in the mapping sample, which leads to problems of choice as to the events interpreted by the geophysicist. This can cause discontinuities and aberrations in the pointing of the surface.
  • the invention overcomes these shortcomings.
  • the tool it offers makes it possible to validate the consistency of the analyzes carried out in relation to the geophysical laws on which the quantities that are taken into account depend.
  • a method of processing seismic data is proposed, according to claim 1.
  • a pointing and / or a pointing propagation is implemented on at least one subset of this data which corresponds to a collection of these sampled data.
  • the magnitude or quantities which correspond to the other dimension (s) then being of fixed values, then a projection law is applied to at least one of the points highlighted by this dot and / or this propagation of dots given in order to determine a point which corresponds to it in another subset corresponding to a collection of data sampled according to the same dimensions, for a value different from at least one quantity which corresponds to at least one other dimension.
  • projection law is meant a law intended to connect the points of one of the subsets to those of another, this law being a function of the modifications that are made to the quantity or quantities corresponding to dimensions other than those in which said subsets are represented.
  • a pointing propagation is implemented from at least one point projected into such another sub-assembly.
  • a succession of alternating dots propagations and projections is implemented.
  • an inverse projection law is applied to at least one of the points determined by point propagation in such another subset in order to determine a point which corresponds to it in at least one subset from which we have implemented a projection.
  • a map of the points pointed and / or projected and / or determined by propagation of points is displayed in the different subsets considered.
  • a quantity which is varied to define another dimension is advantageously a source / receiver offset and / or an angle of incidence and / or an acquisition time and / or a type of seismic waves and / or a speed migration and / or a parameter taking into account at least one of these parameters.
  • the number of dimensions taken into account is greater than four.
  • the invention also relates to a method of cartographic determination on seismic data, characterized in that the aforementioned method is implemented to determine a cartography according to at least four dimensions.
  • a neighborhood of points is determined from a projected and / or pointed and / or determined by propagation of dots a point comprising in addition points in a predetermined neighborhood of said point in the subassembly of dimensions two or three to which said point belongs, one or more points projected from said point into other subsets and in that one determines on the set of points of this neighborhood that of these points which optimizes a given criterion.
  • a method of interpreting seismic data in a space of dimension greater than three is proposed, the first three dimensions being advantageously, but not limited to, the spatial dimensions (X, Y and vertical axis).
  • a characteristic of this process is to alternate phases of propagation, i.e. mapping of surfaces in 3D space
  • FIG. 1 is a representation of a three-dimensional volume in which various seismic traces obtained have been pointed out
  • FIG. 2 represents three-dimensional volumes vi, v 2 , and v 3 corresponding to different values of a parameter Xd, - FIG. 3 illustrates an example of surface pointing guide
  • FIGS. 4 and 5 illustrate examples of dotted propagation in three three-dimensional volumes vi, v 2 , v 3 ,
  • recordings are processed to be used in the form of data (for example, but not necessarily the amplitude values of the recorded waves) depending on a number "n" of parameters which can be much greater than 3.
  • xi, x 2 , x n in the remainder of the text. They can correspond to the geometric dimensions X, Y, Z (position and depth or time) and / or to all the other quantities which one wishes to take into account to analyze the records available: source / receiver distance, angle d incidence, time between acquisition, wave type, speed model increment, etc.
  • a first step we choose three of these parameters (xi, x 2 ... x n for example) to define three dimensions of a subspace three-dimensional on which it is possible to carry out a classic 3D cartography, if the other parameters are fixed.
  • the data that we consider are generally collections of 3D seismic data each generated for a given value of one or more other parameters, which we vary to define additional dimensions.
  • FIG. 1 illustrates an example of surface pointing (surface Si) produced on such a volume Vi.
  • one or more so-called projection relationships or “surface pointing guides” are also defined intended to connect the points of one of these volumes to those of another, according to the modifications that are made to the parameters other than those (xi, X2 ... x n ) which correspond to the dimensions in which these volumes are represented, to pass from one to the other of these two volumes.
  • pointing guides can be determined in different ways and in particular in addition to the values of the samples considered, also according to dynamic or kinematic descriptions of the evolution of the volumes in the directions corresponding to the modified parameters.
  • Pointing guides can also be laws that are given empirically.
  • Figures 2 and 3 illustrate an example of what a pointing guide can be.
  • Vi, v 2 and v 3 correspond to three 3D volumes giving collections of samples as a function of the same three parameters xi, x 2 ... x n for three different values of another parameter Xd (the term "value" being here to understand in a broad way, Xd can just as well correspond to a quantifiable parameter, as to a parameter expressing a condition, other than a quantified condition (type of P waves, S waves, etc.)).
  • a guide for pointing the surfaces is for example a curve ( ⁇ in FIG.
  • P (i ⁇ i n ) be a point on the surface S and in a direction other than those forming the predefined 3D volume.
  • a surface pointing guide associated with P in a direction d is an application Yl d p which makes it possible to determine a part of the points belonging to S lying in the vicinity of P in the plane defined by (P,
  • the application ⁇ describes a curve representing the projection of the image described by S in the plane (P, Xd, X n ).
  • the curve is defined by segments connecting the extrema correlating from S from one column to another in the plane (P, X d , X n ).
  • R the map connecting two neighboring points of S by a segment.
  • the inverse relation R "1 is then applied so as to transform the segment into elementary surface elements.
  • the ⁇ application thus makes it possible to pass from one 3D volume to another, to point in a volume then in another of the events in correspondence having regard to the selected propagation criteria. Consistency tests are advantageously provided to validate the points that we retain in the volumes in which we project using the guide ⁇ : likeness and / or cross-correlation calculation for example. It will be understood that the technique which has just been described makes it possible successively to implement a surface pointing in a first volume, to pass into a second volume by a pointing guide, then to resume a surface pointing in this second.
  • FIG. 4 illustrates in this sense, in the case of an analysis taking into account 4 dimensions, combinations of pointings of conventional 2D or 3D surfaces (small arrows on the grids) with propagation guides making it possible to pass from a volume 3D to another.
  • the techniques proposed are for example used to circumvent a discontinuity (referenced by D) preventing the pursuit of a surface dotting in a volume.
  • the surface tip is produced in a volume of dimension Xi, X 2 , X *.
  • a discontinuity D prevents this surface pointing over the whole of this volume.
  • One or more points highlighted by the surface dotted line are then projected into one or more other volumes corresponding to one or more other values for a fourth parameter (in this case X 3 ).
  • a fourth parameter in this case X 3 .
  • Example 1 Some examples of possible surface pointing guides are given below.
  • Example 1 Some examples of possible surface pointing guides are given below.
  • a surface pointing guide in the direction X 3 can be:
  • Xi, X 2 corresponding to the ground position parameter, resulting from an acquisition whose source / receiver distance is regularly sampled with a step ⁇ 3 .
  • X 3 represents the offset distance
  • X 4 the time, sampled with a step ⁇ .
  • the surface pointing guide depends on the average speed in the basement, which is an external parameter to the volume that we are looking for to image. This parameter provides description information in the direction X 3 .
  • the speed disturbance being linear, we can consider as a first approximation that for a position Xi, X 2 , the depth X 4 of the sample associated with a reflector evolves linearly from one seismic volume to another, in the direction X 3 .
  • the surface pointing guide is therefore defined as a linear function (represented in FIG. 7 of the type
  • ⁇ 3 is the sampling step in the direction X3.
  • the point guide is corrected so that it is consistent with the points obtained. For example, each new point in a three-dimensional volume constitutes a new point in the guide ⁇ .
  • the proposed technique makes it possible to take into account existing 2D and 3D propagation methods and to extend them to “volumes” of dimensions 4 or greater. This increases the number of possible extension paths, in particular for surface dots.
  • the projection relationships or pointing guides can depend as much on information describing the dynamics or the kinematics of the seismic environment as on empirical information obtained from samples of the acquisition or from external information.

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)

Abstract

Procédé de traitement de données sismiques, caractérisé en ce qu'en vue de leur analyse en fonction d'au moins quatre dimensions, on met en oeuvre un pointé et/ou une propagation de pointé sur au moins un sous-ensemble de ces données qui correspond à une collection ces données échantillonnées selon trois ou deux dimensions, la ou les grandeurs qui correspondent à la ou aux autres dimensions étant alors de valeurs fixes, puis on applique à au moins un des points mis en évidence par ce pointé et/ou cette propagation de pointés une loi de projection donnée afin de déterminer un point qui lui correspond dans un autre sous-ensemble correspondant à une collection de données échantillonnées selon les mêmes dimensions, pour une valeur différente d'au moins une grandeur qui correspond à au moins une autre dimension.

Description

PERFECTIONNEMENTS AUX TECHNIQUES DE POINTE DE SURFACE
La présente invention concerne des perfectionnements aux techniques de pointé de surfaces utilisées en traitement sismique. Dans le domaine des sciences de la terre, pour imager la structure du sous-sol, le géophysicien utilise les techniques de « sismique réflexion ». Celles-ci consistent à émettre à la surface du sol des signaux acoustiques et à les enregistrer après leurs réflexions successives aux limites des couches géologiques. Ces couches constituent le sous-sol dont on cherche à créer un modèle géométrique.
Pour ce faire, l'un des procédés le plus couramment utilisé, notamment dans la recherche de gisements pétroliers est le procédé « d'acquisition 3D ». Celui-ci consiste à répartir des émetteurs et des récepteurs dans un plan horizontal en surface le long d'une grille de coordonnées X, Y.
Le résultat d'un tel procédé est une image tridimensionnelle du sous-sol organisée sous la forme d'enregistrements échantillonnés en fonction du temps ou de la profondeur, noté Z. Chacun de ces enregistrements, appelé « trace », est un représentatif de la verticale à l'aplomb d'un point Py de coordonnées Xj, Yj.
Une telle acquisition tridimensionnelle permet d'obtenir un volume dont les dimensions sont deux axes géométriques X, Y et un axe des temps ou de profondeur d'acquisition Z. Outres ces trois dimensions spatiales (X, Y et axe vertical Z), bien d'autres grandeurs ont un intérêt pour le physicien, et notamment :
- La distance source/récepteur ou l'angle d'incidence qui a servi à générer le volume. Jusqu'à il y a peu cette distance (déport) était le plus souvent ramené à 0, les volumes étaient obtenus par sommation de contributions à distance constante de façon à augmenter le rapport signal/bruit. La sommation partielle par classes de déport ou d'angle permet désormais au géophysicien de faire des analyses par classe de distance ou d'angle. - Le temps. Plusieurs acquisitions sismiques peuvent être effectuées successivement dans le temps avec un même réseau de récepteurs /émetteurs ; la comparaison des données acquises permet, par exemple, de comprendre l'évolution d'un régime de fluide dans le cas d'une exploitation de réservoir.
- Le type d'ondes sismigues. Les capteurs peuvent enregistrer simultanément les ondes de compression (ondes P) et les ondes de cisaillement (ondes S) ; le rapport entre les temps d'arrivée de ces deux types d'ondes est un lien direct avec la nature du milieu où les ondes se propagent.
- La vitesse. La migration temps ou profondeur consiste à déplacer les points du volume d'acquisition en tenant compte d'un champ de vitesse. La vitesse utilisée dans ce processus peut être perturbée et constituer une dimension d'analyse complémentaire. Par exemple, pour mettre à jour un modèle de vitesse, la demanderesse a proposé une approche consistant à distribuer de manière régulière les facteurs de perturbation sur un modèle de vitesse, et de s'intéresser pour une position X, Y, à la valeur de l'échantillon associé à un même réflecteur dans chacun des cubes migres.
- etc. La prise en compte de ces autres grandeurs est fondamentale dans nombre de domaines des géo-sciences, de l'analyse de vitesse à la caractérisation de réservoir.
Elle est aujourd'hui accessible grâce à l'augmentation de la puissance de calcul informatique et des capacités de stockage ; il reste à disposer des outils d'analyse et d'interprétation adaptés à cette nouvelle orientation.
Aujourd'hui, lorsque l'on prend en compte avec les outils existants une grandeur en plus des trois dimensions spatiales dans lesquelles les volumes sismiques sont généralement représentés, l'interprétation des données se ramène à la somme des interprétations réalisées de façon indépendantes sur chacun des volumes sismiques 3D obtenus pour différentes valeurs de ladite grandeur. Le document WO 99/67660 décrit un procédé permettant de suivre les évolutions du sous-sol. Ce procédé consiste à déterminer un champ de déplacement permettant de passer d'un volume échantillonné de données sismiques mesurées à un temps T0 à un champ de données sismiques mesurées à un temps T ultérieur. Le champ de déplacement ainsi déterminé est appliqué au premier volume sismique de manière à obtenir une image du sous-sol.
Les informations que l'on peut retirer d'un tel mode d'interprétation sont nécessairement limitées. L'invention propose quant à elle une contribution efficace et originale à la résolution de ce problème.
Une technique classiquement utilisée pour l'analyse de volumes sismiques 3D est celle du pointé de surface, dont un exemple est donné dans le brevet FR 2 657 695. Ce document concerne un procédé de pointé de surface dans un volume tridimensionnel selon lequel on transforme un volume représentant les valeurs d'un paramètre échantillonné en une image à deux dimensions en utilisant une courbe de Peano-Hilbert et on applique un programme de pointé automatique à l'image 2D obtenue.
Lorsque l'on prend en compte une grandeur supplémentaire en plus des dimensions spatiales qui constituent les volumes 3D (X, Y et axe vertical) et que l'on fait varier ladite grandeur, une difficulté tient en ce que la nature et la forme du réflecteur à interpréter peut changer d'un volume à l'autre le long de la dimension que constitue ladite grandeur.
La cohérence entre les différents volumes 3D n'est pas forcément inscrite dans l'échantillon cartographie, ce qui conduit à des problèmes de choix quant aux événements interprétés par le géophysicien. Cela peut entraîner des discontinuités et des aberrations dans le pointé de la surface.
L'invention permet de pallier ces insuffisances.
Elle permet des interprétations prenant en compte l'ensemble des paramètres à analyser et tenant compte des lois qui les caractérisent.
L'outil qu'elle propose permet en outre de valider la cohérence des analyses réalisées par rapport aux lois géophysiques dont dépendent les grandeurs que l'on prend en compte. A cet effet, il est proposé un procédé de traitement de données sismiques, selon la revendication 1.
Selon ce procédé, en vue de leur analyse en fonction d'au moins quatre dimensions, on met en œuvre un pointé et/ou une propagation de pointé sur au moins un sous-ensemble de ces données qui correspond à une collection de ces données échantillonnées selon trois ou deux dimensions, la ou les grandeurs qui correspondent à la ou aux autres dimensions étant alors de valeurs fixes, puis on applique à au moins un des points mis en évidence par ce pointé et/ou cette propagation de pointés une loi de projection donnée afin de déterminer un point qui lui correspond dans un autre sous-ensemble correspondant à une collection de données échantillonnées selon les mêmes dimensions, pour une valeur différente d'au moins une grandeur qui correspond à au moins une autre dimension.
Par loi de projection, on entend une loi destinée à relier les points d'un des sous-ensembles à ceux d'un autre, cette loi étant fonction des modifications que l'on apporte à la ou aux grandeurs correspondant aux dimensions autres que celles dans lesquelles lesdits sous-ensembles sont représentés.
Des exemples de loi de projection sont donnés plus loin de façon détaillée.
Avantageusement, on met en œuvre une propagation de pointé à partir d'au moins un point projeté dans un tel autre sous-ensemble.
Notamment, de façon préférée, on met en œuvre une succession d'alternance de propagations de pointés et de projections. De façon avantageuse également, on applique à au moins un des points déterminé par propagation de pointé dans un tel autre sous- ensemble une loi de projection inverse afin de déterminer un point qui lui correspond dans au moins un sous-ensemble à partir duquel on a mis en œuvre une projection. Egalement, on met préférentiellement on met en œuvre sur un point projeté un test de cohérence pour le valider. En particulier, on affiche une cartographie des points pointés et/ou projetés et/ou déterminés par propagation de pointés dans les différents sous-ensembles considérés.
Une grandeur que l'on fait varier pour définir une autre dimension est avantageusement un déport source/récepteur et/ou un angle d'incidence et/ou un temps d'acquisition et/ou un type d'ondes sismiques et/ou une vitesse de migration et/ou un paramètre prenant en compte au moins l'un de ces paramètres.
De façon particulièrement préférée, le nombre de dimensions prises en compte est supérieur à quatre.
L'invention concerne également un procédé de détermination cartographique sur des données sismiques, caractérisé en ce qu'on met en œuvre le procédé précité pour déterminer une cartographie selon au, moins quatre dimensions. Avantageusement notamment, on détermine à partir d'un point projeté et/ou pointé et/ou déterminé par propagation de pointés un voisinage de points comportant en plus des points dans un voisinage prédéterminé dudit point dans le sous-ensemble de dimensions deux ou trois auquel ledit point appartient, un ou des points projetés à partir dudit point dans d'autres sous-ensembles et en ce qu'on détermine sur l'ensemble des points de ce voisinage celui de ces points qui optimise un critère donné.
Autrement formulé, il est proposé un un procédé d'interprétation des données sismiques dans une espace de dimension supérieure à trois, les trois premières dimensions étant avantageusement, mais non limitativement, les dimensions spatiales (X,Y et axe vertical).
Une caractéristique de ce procédé est d'alterner des phases de propagation, c'est à dire de cartographie des surfaces dans l'espace 3D
(X,Y et axe vertical par exemple) et des phases de projection durant lesquelles on cherche à déplacer l'échantillon suivant une. autre dimension, toutes choses égales
Par ailleurs, une fois trouvé dans l'espace à N dimensions une position acceptable pour l'échantillon - c'est à dire remplissant un ensemble de critères de validation eu égard à l'échantillon d'origine --, la propagation dans l'espace des dimensions spatiale peut reprendre, les (N-3) dimensions restantes étant de nouveau figées.
L'ensemble des échantillons validés durant ces alternances de phases de propagation et de projection est cartographie au fur et à mesure que progresse l'algorithme de cartographie de l'hyper-surface dans l'espace à N dimensions.
D'autres caractéristiques et avantages de l'invention ressortiront encore de la description qui suit, laquelle est purement illustrative et non limitative et doit être lue en regard des figures annexées parmi lesquelles :
- la figure 1 est une représentation d'un volume tridimensionnel dans lequel on a pointé différentes traces sismiques obtenues,
- la figure 2 représente des volumes tridimensionnels v-i, v2, et v3 correspondant à différentes valeurs d'un paramètre Xd, - la figure 3 illustre un exemple de guide de pointé de surface,
- les figures 4 et 5 illustrent des exemples de propagation de pointé dans trois volumes tridimensionnels vi, v2, v3,
- les figures 6 à 8 illustrent des exemples de guides de pointé de surface.
On se place dans la suite du texte dans l'hypothèse où un ensemble de traces sismiques a été acquis à l'occasion d'une ou plusieurs campagnes d'enregistrement, par exemple d'enregistrement 3D.
Ces enregistrements sont traités pour être exploités sous forme de données (par exemple, mais non nécessairement les valeurs d'amplitudes des ondes enregistrées) dépendant d'un nombre « n » de paramètres qui peut être bien supérieur à 3.
Ces paramètres sont appelés x-i, x2, xn dans la suite du texte. Ils peuvent correspondre aux dimensions géométriques X, Y, Z (position et profondeur ou temps) et/ou à- toutes les autres grandeurs que l'on souhaite prendre en compte pour analyser les enregistrements dont on dispose : distance source/récepteur, angle d'incidence, temps entre acquisition, type d'onde, incrément de modèle de vitesse, etc.
Dans une première étape on choisit trois de ces paramètres (x-i, x2...xn par exemple) pour définir trois dimensions d'un sous espace tridimensionnel sur lequel il est possible de réaliser une cartographie 3D classique, si les autres paramètres sont fixés.
Les données que l'on considère sont en général des collections de données sismiques 3D générées chacune pour une valeur donnée d'un ou plusieurs autres paramètres, que l'on fait varier pour définir des dimensions supplémentaires.
Pour chaque volume 3D ainsi obtenu, il est possible de réaliser des cartographies, par exemple par pointé de surfaces.
La figure 1 illustre un exemple de pointé de surface (surface Si) réalisé sur un tel volume Vi .
Par ailleurs, on définit également une ou plusieurs relations dites de projection ou de « guides de pointé de surface » destinées à relier les points d'un de ces volumes à ceux d'un autre, en fonction des modifications que l'on apporte aux paramètres autres que ceux (xi, X2...xn) qui correspondent aux dimensions dans lesquelles ces volumes sont représentés, pour passer de l'un à l'autre de ces deux volumes.
Ces guides de pointé peuvent être déterminés de différentes façons et notamment outre en fonction des valeurs des échantillons considérés, en fonction également de descriptions dynamiques ou cinématiques de l'évolution des volumes dans les directions correspondant aux paramètres modifiés.
Les guides de pointé peuvent également être des lois que l'on se donne de façon empirique.
Les figures 2 et 3 illustrent un exemple de ce que peut être un guide de pointé.
V-i, v2 et v3 correspondent à trois volumes 3D donnant des collections d'échantillons en fonction des trois mêmes paramètres x-i, x2...xn pour trois valeurs différentes d'un autre paramètre Xd (le terme « valeur » étant ici à comprendre de façon large, Xd pouvant tout aussi bien correspondre à un paramètre quantifiable, qu'à un paramètre exprimant une condition, autre qu'une condition quantifiée (type d'ondes ondes P, ondes S, etc.)). Si étant une surface pointée autour du point Pi dans le volume Vi, un guide de pointé de surfaces est par exemple une courbe (π sur la figure 3) qui dans le plan (P, Xd, X1 ) permet d'associer à un point Pi d'un pointé de surface dans le volume vi des points P2 et P3 qui dans les volumes V2 et v3 décrivent théoriquement le même événement, pour les mêmes valeurs de X1 , X2.
Considérons une surface S associe au volume de dimension n. Elle est définie par : v (i13.», ι)e[1 .5Nι]x...x[ JN„-ι]
Figure imgf000010_0001
Soit P(iι in) un point de la surface S et d une direction autre que celles formant le prédéfini volume 3D.
Un guide de pointé de surfaces associé à P dans une direction d est une application Yld p qui permet de déterminer une partie des points appartenant à S se trouvant au voisinage de P dans le plan défini par (P,
Si on sait définir l'application π pour un point P dans une direction d, on sait étendre la surface S à partir de P dans cette direction. On augmente ainsi la dimension du voisinage, dans le procédé de propagation, par le nombre de direction, autres que celles formant le prédéfini volume
3D, où une application π peut être définie.
L'application π décrit une courbe représentant la projection de l'image décrite par S dans le plan (P, Xd, Xn). La courbe est définie par des segments reliant les extrema corrélant de S d'une colonne à l'autre dans le plan (P, Xd, Xn). Notons R l'application reliant deux points voisins de S par un segment. On applique alors la relation inverse R"1 de façon à transformer le segment en éléments de surface élémentaire.
L'application π permet ainsi de passer d'un volume 3D à un autre, pour pointer dans un volume puis dans un autre des événements en correspondance eu égard aux critères de propagation sélectionnés. Des tests de cohérence sont avantageusement prévus pour valider les points que l'on retient dans les volumes dans lesquels on se projette à l'aide du guide π : calcul de semblance et/ou de cross-corrélation par exemple. On comprend que la technique qui vient d'être décrite permet successivement de mettre en œuvre un pointé de surface dans un premier volume, de passer dans un deuxième volume par un guide de pointé, puis de reprendre un pointé de surface dans ce deuxième.
En multipliant ce type de combinaisons, on augmente considérablement le nombre de cheminements possibles à partir d'un échantillon pointé, pour parvenir à la cartographie complète de l'événement sur l'ensemble de l'espace.
La figure 4 illustre en ce sens, dans le cas d'une analyse prenant en compte 4 dimensions, des combinaisons de pointés de surfaces 2D ou 3D classiques (petites flèches sur les grilles) avec des guides de propagation permettant de passer d'un volume 3D à un autre.
Dans cet exemple, on a illustré un cas où l'on pend en couple quatre paramètres, l'un de ces paramètres étant celui que l'on fait varier pour passer d'un volume à un autre. II est bien entendu possible de prendre en compte un nombre de paramètres beaucoup plus important et de faire varier plusieurs paramètres.
Il est également possible de changer les paramètres qui définissent les dimensions des volumes 3D dans lesquels les techniques classiques de pointé de surface sont mises en œuvre. Comme l'illustre la figure 5, les techniques proposées sont par exemple utilisées pour contourner une discontinuité (référencée par D) empêchant de poursuivre un pointé de surface dans un volume.
Ainsi, dans l'exemple illustré sur ladite figure 5, le pointe de surface est réalisé dans un volume de dimension Xi, X2, X*. Une discontinuité D empêche de réaliser ce pointé de surface sur l'ensemble de ce volume. On projette alors un ou plusieurs points mis en évidence par le pointé de surface dans un ou plusieurs autres volumes correspondant à une ou plusieurs autres valeurs pour un quatrième paramètre (en l'occurrence X3). Dans chacun de ces autres volumes, il est possible de mettre en œuvre à nouveau des pointés de surfaces classiques, et ainsi de contourner la faille pour à partir d'un point mis en évidence par un de ces pointés de surface poursuivi dans un de ces autres volumes, revenir dans le premier volume par une projection inverse (point P2 et Jj).
On donne ci-après quelques exemples de guides de pointés de surfaces possibles. Exemple 1.
On considère des sismiques tridimensionnelles corrigées au déport, Xi et X2 étant des dimensions correspondant au paramètre de position en surface, X3 correspondant au paramètre de déport, X4 correspondant au temps corrigé. Un guide de pointé de surface dans la direction X3 peut être :
Up(j) = Adj(i4) où Adj est la fonction d'ajustement au signal sur la trace courante d'un point de la surface et où i4 est une valeur échantillonnée selon la direction X . Exemple 2.
On considère le volume sismique représenté selon deux directions
Xi, X2 correspondant au paramètre de position au sol, issu d'une acquisition dont la distance source/récepteur est échantillonnée régulièrement avec un pas Δ3. X3 représente la distance de déport, X4 le temps, échantillonné avec un pas Δ .
On dispose pour le réflecteur à cartographier d'une carte de vitesse V( , i2), vitesse de propagation de l'onde au point (i1 ; i2, i4) du sous-sol. Si on note Adj la fonction d'ajustement du signal sur la trace courante d'un point de la surface, l'utilisation de la loi de Dix permet de définir, pour un point P(iι, Î2, Î3, Î4) de la surface, un guide de propagation dans la direction X3 des offsets par la fonction (représentée à la figure 6)
Figure imgf000012_0001
Ici, le guide de pointé de surfaces dépend de la vitesse moyenne dans le sous-sol, qui est un paramètre externe au volume que l'on cherche à imager. Ce paramètre apporte une information de description dans la direction X3.
Exemple 3. Dans le cadre de la mise à jour d'un modèle de vitesse par la méthode de balayage ou « scan » de vitesse, on fait varier la vitesse initiale d'un facteur de perturbation α. Après les différentes migrations profondeur, on obtient n volumes sismiques échantillonés suivant la direction X3 correspondant au facteur de perturbation i3=α. La perturbation de vitesse étant linéaire, on peut considérer en première approximation que pour une position Xi, X2, la profondeur X4 de l'échantillon associé à un réflecteur évolue de manière linéaire d'un volume sismique à l'autre, dans la direction X3. Le guide de pointé de surface est donc défini comme une fonction linéaire (représentée à la figure 7 du type
Figure imgf000013_0001
où Δ3 est le pas d'échantillonnage selon la direction X3. Au fur et à mesure que l'on projette les pointés d'un volume sismique à l'autre et que l'on valide les pointés obtenus, on corrige le guide de pointé pour qu'il soit cohérents avec les pointés obtenus. Par exemple, chaque nouveau pointé dans un volume tridimensionnel constitue un nouveau point du guide π. Par interpolation ou par régression, on peut modifier le guide de pointé π pour qu'il corresponde aux points obtenus. Ceci peut permettre dans certains cas, de créer de façon incrémentale le guide de pointé ou d'en déduire certaines variables externes au volume V quantifiables analytiquement et apportant une information descriptive des relations entre les points de surface dans une ou plusieurs des directions.
Comme on l'aura compris, la technique proposée permet de prendre en compte des méthodes de propagation 2D et 3D existantes et de les étendre à des « volumes » de dimensions 4 ou supérieure. On augmente ainsi le nombre de chemins d'extension possibles notamment pour les pointés de surface. On notera que les relations de projection ou guides de pointés peuvent dépendre aussi bien d'informations décrivant la dynamique ou la cinématique de l'environnement sismique que d'informations empiriques obtenues à partir des échantillons de l'acquisition ou d'informations externes.

Claims

REVENDICATIONS
1. Procédé de traitement de données sismiques, caractérisé en ce qu'en vue de leur analyse en fonction d'au moins quatre dimensions, on met en œuvre un pointé et/ou une propagation de pointé sur au moins un sous- ensemble de ces données qui correspond à une collection de ces données échantillonnées selon trois ou deux dimensions, la ou les grandeurs qui correspondent à la ou aux autres dimensions étant alors de valeurs fixes, puis on applique à au moins un des points mis en évidence par ce pointé et/ou cette propagation de pointés une loi de projection donnée afin de déterminer un point qui lui correspond dans un autre sous-ensemble correspondant à une collection de données échantillonnées selon les mêmes dimensions, pour une valeur différente d'au moins une grandeur qui correspond à au moins une autre dimension.
2. Procédé selon la revendication 1 , caractérisé en ce qu'on met en œuvre une propagation de pointé à partir d'au moins un point projeté dans un tel autre sous-ensemble.
3. Procédé selon l'une des revendications 1 ou 2, caractérisé en ce qu'on met en œuvre une succession d'alternance de propagations de pointés et de projections.
4. Procédé selon l'une des revendications 2 ou 3, caractérisé en ce qu'on applique à au moins un des points déterminé par propagation de pointé dans un tel autre sous-ensemble une loi de projection inverse afin de déterminer un point qui lui correspond dans au moins un sous-ensemble à partir duquel on a mis en œuvre une projection.
5. Procédé selon l'une des revendications précédentes, caractérisé en ce qu'on met en œuvre sur un point projeté un test de cohérence pour le valider.
6. Procédé selon l'une des revendications précédentes, caractérisé en ce qu'on affiche une cartographie des points pointés et/ou projetés et/ou déterminés par propagation de pointés dans les différents sous-ensembles considérés.
7. Procédé selon l'une des revendications précédentes, caractérisé en ce qu'un paramètre modifié d'un sous-ensemble à un autre est un déport source/récepteur et/ou un angle d'incidence et/ou un temps d'acquisition et/ou un type d'ondes sismiques et/ou une vitesse de migration et/ou un paramètre prenant en compte au moins l'un de ces paramètres.
8. Procédé selon l'une des revendications précédentes, caractérisé en ce que le nombre de dimensions prises en compte est supérieur à quatre.
9. Procédé de détermination cartographique sur des données sismiques, caractérisé en ce qu'on met en œuvre le procédé selon l'une des revendications précédentes pour déterminer une cartographie selon au, moins quatre dimensions.
10. Procédé selon la revendication 9, caractérisé en ce qu'on détermine à partir d'un point projeté et/ou pointé et/ou déterminé par propagation de pointés un voisinage de points comportant en plus des points dans un voisinage prédéterminé dudit point dans le sous-ensemble de dimensions deux ou trois auquel ledit point appartient, un ou des points projetés à partir dudit point dans d'autres sous-ensembles et en ce qu'on détermine sur l'ensemble des points de ce voisinage celui de ces points qui optimise un critère donné.
11. Procédé selon l'une des revendications précédentes, caractérisé en ce que dans le cas de données sismiques tridimensionnelles corrigées au déport, les dimensions prises en compte comprennent deux dimensions correspondant aux paramètres de position en surface (Xi, X2), une dimension ( X3) correspondant au paramètre de déport, une dimension (XΦ) correspondant au temps corrigé et en ce qu'une loi de propagation selon la direction la direction correspondant au paramètre de déport (X3) est :
Up{j) = Adj( où Adj est la fonction d'ajustement au signal sur la trace courante d'un point de la surface et où i est une valeur échantillonnée selon la dimension (X4) correspondant au temps corrigé.
12. Procédé selon l'une des revendications 1 à 10, caractérisé en ce que les dimensions prises en compte comprennent deux dimensions correspondant aux paramètres de position en surface (Xi, X2), une dimension ( X3) correspondant au paramètre de déport, une dimension (X4) correspondant au temps, et en ce que les données étant issues d'une acquisition avec une distance source/récepteur échantillonnée régulièrement avec un pas Δ3, et un temps échantillonné avec un pas Δ4, la loi de projection dans la direction du déport est
Figure imgf000017_0001
ou V(ι'ι, 12) est une vitesse de propagation de l'onde dont on dispose au point (i-ι, i2, i4) du sous-sol, où Adj est la fonction d'ajustement du signal sur la trace courante d'un point de la surface.
13. Procédé selon l'une des revendications 1 à 10, caractérisé en ce que pour mettre à jour un modèle de vitesse par un balayage de vitesse, on fait varier la vitesse initiale d'un facteur de perturbation (α) et en ce que après les différentes migrations profondeur, on obtient n volumes sismiques échantillonnés suivant une dimension (X3) correspondant au facteur de perturbation, en ce que d'autres dimensions prises en compte comprennent deux dimensions correspondant aux paramètres de position en surface (X-i, X2) et une dimension (X ) correspondant à la profondeur de l'échantillon, et en ce qu'une loi de propagation selon la dimension correspondant au facteur de propagation est du type:
Figure imgf000017_0002
où Δ3 est le pas d'échantillonnage selon cette dimension et où Adj est la fonction d'ajustement du signal sur la trace courante d'un point de la surface.
PCT/FR2002/002278 2001-06-29 2002-07-01 Perfectionnements aux techniques de pointe de surface WO2003003052A1 (fr)

Priority Applications (5)

Application Number Priority Date Filing Date Title
US10/482,583 US8014951B2 (en) 2001-06-29 2002-07-01 Surface pointing techniques
MXPA04000063A MXPA04000063A (es) 2001-06-29 2002-07-01 Mejoramientos a tecnicas de punteo superficial.
GB0330107A GB2393515B (en) 2001-06-29 2002-07-01 Improvements to surface pointing techniques
AU2002333947A AU2002333947B2 (en) 2001-06-29 2002-07-01 Improvements to surface picking techniques
CA2452206A CA2452206C (fr) 2001-06-29 2002-07-01 Perfectionnements aux techniques de pointe de surface

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR01/08703 2001-06-29
FR0108703A FR2826734B1 (fr) 2001-06-29 2001-06-29 Perfectionnements aux techniques de pointe de surface

Publications (1)

Publication Number Publication Date
WO2003003052A1 true WO2003003052A1 (fr) 2003-01-09

Family

ID=8864993

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/FR2002/002278 WO2003003052A1 (fr) 2001-06-29 2002-07-01 Perfectionnements aux techniques de pointe de surface

Country Status (8)

Country Link
US (1) US8014951B2 (fr)
CN (1) CN1547673A (fr)
AU (1) AU2002333947B2 (fr)
CA (1) CA2452206C (fr)
FR (1) FR2826734B1 (fr)
GB (1) GB2393515B (fr)
MX (1) MXPA04000063A (fr)
WO (1) WO2003003052A1 (fr)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1452889B1 (fr) * 2003-02-04 2012-03-14 Institut Français du Pétrole Méthode pour mesurer les similarités locales entre plusieurs cubes de traces sismiques

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2429526A (en) * 2005-08-24 2007-02-28 Statoil Asa Identifying continuous features in a seismic data set
US10390892B2 (en) * 2017-08-24 2019-08-27 Synaptive Medical Inc. System and methods for updating patient registration during surface trace acquisition
US10133854B2 (en) 2016-05-12 2018-11-20 International Business Machines Corporation Compositional three-dimensional surface plots

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2657695A1 (fr) * 1990-01-30 1991-08-02 Elf Aquitaine Procede de pointe de surfaces dans un volume 3d.
US5537320A (en) * 1994-10-31 1996-07-16 Landmark Graphics Corporation Method and apparatus for identifying fault curves in seismic data
US5615171A (en) * 1991-07-09 1997-03-25 Landmark Graphics Corporation Method and apparatus for finding horizons in 3D seismic data
WO1999067660A1 (fr) * 1998-06-25 1999-12-29 Schlumberger Holdings Limited Procede de traitement de signaux de donnees sismiques echelonnes dans le temps

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US1240328A (en) * 1914-04-02 1917-09-18 Submarine Signal Co Method and apparatus for locating ore-bodies.
US1730411A (en) * 1922-11-11 1929-10-08 Herbert grove dorset
US2329721A (en) * 1938-09-30 1943-09-21 Cons Eng Corp Geophysical exploration system
US2388703A (en) * 1941-09-20 1945-11-13 United Geophysical Company Inc Geological prospecting system
US3529282A (en) * 1968-01-29 1970-09-15 Robert J S Brown Method for determining in-line and cross dip employing cross steering of seismic data
US4849887A (en) * 1987-08-28 1989-07-18 Amoco Corporation Horizon velocity analysis
US5008861A (en) * 1989-03-06 1991-04-16 Amoco Corporation Geophysical exploration by automatically picking and associating stacked seismic sections with regional coherency peaks of velocity spectra
US5793357A (en) * 1992-11-14 1998-08-11 Ivey; Peter Anthony Device and method for determining movement of a surface
US5771170A (en) * 1994-02-14 1998-06-23 Atlantic Richfield Company System and program for locating seismic events during earth fracture propagation
US5458126A (en) * 1994-02-24 1995-10-17 General Electric Company Cardiac functional analysis system employing gradient image segmentation
US6041018A (en) * 1997-11-13 2000-03-21 Colorado School Of Mines Method for correcting amplitude and phase differences between time-lapse seismic surveys
US6826483B1 (en) * 1999-10-13 2004-11-30 The Trustees Of Columbia University In The City Of New York Petroleum reservoir simulation and characterization system and method
GB2372567B (en) * 2001-02-22 2003-04-09 Schlumberger Holdings Estimating subsurface subsidence and compaction

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2657695A1 (fr) * 1990-01-30 1991-08-02 Elf Aquitaine Procede de pointe de surfaces dans un volume 3d.
US5615171A (en) * 1991-07-09 1997-03-25 Landmark Graphics Corporation Method and apparatus for finding horizons in 3D seismic data
US5537320A (en) * 1994-10-31 1996-07-16 Landmark Graphics Corporation Method and apparatus for identifying fault curves in seismic data
WO1999067660A1 (fr) * 1998-06-25 1999-12-29 Schlumberger Holdings Limited Procede de traitement de signaux de donnees sismiques echelonnes dans le temps

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1452889B1 (fr) * 2003-02-04 2012-03-14 Institut Français du Pétrole Méthode pour mesurer les similarités locales entre plusieurs cubes de traces sismiques

Also Published As

Publication number Publication date
CN1547673A (zh) 2004-11-17
US20040172198A1 (en) 2004-09-02
MXPA04000063A (es) 2004-05-21
FR2826734B1 (fr) 2004-01-16
GB0330107D0 (en) 2004-02-04
CA2452206C (fr) 2014-08-26
GB2393515B (en) 2004-10-13
GB2393515A (en) 2004-03-31
AU2002333947B2 (en) 2007-11-15
CA2452206A1 (fr) 2003-01-09
FR2826734A1 (fr) 2003-01-03
US8014951B2 (en) 2011-09-06

Similar Documents

Publication Publication Date Title
EP0889331B1 (fr) Méthode pour modéliser en 3D l'impédance d'un milieu hétérogène
CA2528700C (fr) Methode pour determiner des informations speculaires apres imagerie sismique avant sommation
US20130028049A1 (en) Method for handling rough sea and irregular recording conditions in multi-sensor towed streamer data
FR2843202A1 (fr) Methode pour former un modele representatif de la distribution d'une grandeur physique dans une zone souterraine, affranchi de l'effet de bruits correles entachant des donnees d'exploration
FR2979711A1 (fr) Procede de traitement de donnees sismiques marines
WO2006096348A1 (fr) Suppression de bruit de donnees sismiques au moyen de transformees de radon
CN103293551A (zh) 一种基于模型约束的阻抗反演方法及系统
EP2075597A2 (fr) Conditionnement spectral pour données sismiques de surface
EP2537049A2 (fr) Estimation de réflexions multiples internes dans des données sismiques
FR2747477A1 (fr) Methode d'inversion pour donnees sismiques
CN107179550A (zh) 一种数据驱动的地震信号零相位反褶积方法
FR2923312A1 (fr) Procede de traitement d'images sismiques du sous-sol
CA2464799C (fr) Methode pour determiner un modele de vitesse d'ondes sismiques dans une formation souterraine heterogene
WO2003003052A1 (fr) Perfectionnements aux techniques de pointe de surface
Bakulin et al. Evaluating strategies for estimation of local kinematic parameters in noisy land data: quality versus performance trade-offs
FR2884619A1 (fr) Procede de traitement de donnees sismiques en vue de la realisation d'un pointe rmo
Li et al. An efficient deep learning method for VSP wavefield separation: A DAS-VSP case
RU2126984C1 (ru) Способ определения глубинно-скоростных параметров среды и построения ее изображения по сейсмическим данным - система prime
CN114114421A (zh) 基于深度学习的导向自学习地震数据去噪方法及装置
FR2888645A1 (fr) Procede d'acquisition et de traitement de donnees magnetometriques par des mises a jour locales et en temps reel
US10338253B2 (en) Method of suppressing spectral artefacts of wavefield decomposition caused by imperfect extrapolation
FR2801387A1 (fr) Perfectionnements aux procedes de traitement sismique mettant en oeuvre une determination des corrections dynamiques
EP0610123B1 (fr) Procédé pour améliorer l'estimation, par analyse de focalisation, des vitesses de propagation des ondes sismiques
FR3024243A1 (fr) Procede de filtrage adaptatif de reflexions sismiques multiples
Mi et al. Prestack depth imaging from topography with a Fourier method

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NO NZ OM PH PL PT RO RU SD SE SG SI SK SL TJ TM TN TR TT TZ UA UG US UZ VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): GH GM KE LS MW MZ SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR IE IT LU MC NL PT SE SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

Kind code of ref document: A1

Designated state(s): GH GM KE LS MW MZ SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR IE IT LU MC NL PT SE SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG US

ENP Entry into the national phase

Ref document number: 0330107

Country of ref document: GB

Kind code of ref document: A

Free format text: PCT FILING DATE = 20020701

121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 10482583

Country of ref document: US

WWE Wipo information: entry into national phase

Ref document number: 2452206

Country of ref document: CA

WWE Wipo information: entry into national phase

Ref document number: PA/A/2004/000063

Country of ref document: MX

WWE Wipo information: entry into national phase

Ref document number: 2002333947

Country of ref document: AU

WWE Wipo information: entry into national phase

Ref document number: 20028166108

Country of ref document: CN

REG Reference to national code

Ref country code: DE

Ref legal event code: 8642

122 Ep: pct application non-entry in european phase
NENP Non-entry into the national phase

Ref country code: JP

WWW Wipo information: withdrawn in national office

Ref document number: JP

ENP Entry into the national phase

Ref document number: 2002333947

Country of ref document: AU

Date of ref document: 20020701

Kind code of ref document: B