FR2791438A1 - QUALITY CONTROL METHOD OF SEISMIC DATA MIGRATION USING GENERALIZED RADON TRANSFORM - Google Patents
QUALITY CONTROL METHOD OF SEISMIC DATA MIGRATION USING GENERALIZED RADON TRANSFORM Download PDFInfo
- Publication number
- FR2791438A1 FR2791438A1 FR9903590A FR9903590A FR2791438A1 FR 2791438 A1 FR2791438 A1 FR 2791438A1 FR 9903590 A FR9903590 A FR 9903590A FR 9903590 A FR9903590 A FR 9903590A FR 2791438 A1 FR2791438 A1 FR 2791438A1
- Authority
- FR
- France
- Prior art keywords
- qcimage
- point
- space
- migrated image
- radius
- 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
- 230000005012 migration Effects 0.000 title claims abstract description 32
- 238000013508 migration Methods 0.000 title claims abstract description 32
- 229910052704 radon Inorganic materials 0.000 title claims abstract description 9
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 title claims abstract description 9
- 238000003908 quality control method Methods 0.000 title claims abstract description 8
- 238000000034 method Methods 0.000 claims abstract description 38
- 238000004364 calculation method Methods 0.000 description 5
- 238000001914 filtration Methods 0.000 description 5
- 230000002452 interceptive effect Effects 0.000 description 2
- 239000003086 colorant Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- 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
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/42—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators in one well and receivers elsewhere or vice versa
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2200/00—Details of seismic or acoustic prospecting or detecting in general
- G01V2200/10—Miscellaneous details
- G01V2200/14—Quality control
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/16—Survey configurations
- G01V2210/161—Vertical seismic profiling [VSP]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
- G01V2210/46—Radon transform
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Image Processing (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
La présente invention concerne un procédé de contrôle de qualité d'une migration de données sismiques utilisant la transformée de Radon généralisée, permettant de passer d'un espace de données sismiques à un espace de l'image migrée, dans lequel on calcule la table des paramètres permettant de connaître pour une onde sismique allant d'un point de l'image à une source ou à un capteur la longueur du trajet, sa durée, et les angles d'incidence de l'onde au départ et à l'arrivée du trajet, et dans lequel l'on établit une correspondance entre au moins une zone d'un premier de ces deux espaces et au moins une zone du second espace, en utilisant cette table de paramètres pour remplir une table de correspondance QCimage.The present invention relates to a method of quality control of a seismic data migration using the generalized Radon transform, making it possible to pass from a seismic data space to a space of the migrated image, in which the table of values is calculated. parameters making it possible to know for a seismic wave going from a point of the image to a source or to a sensor the length of the path, its duration, and the angles of incidence of the wave at the start and at the arrival of the path, and in which a correspondence is established between at least one zone of a first of these two spaces and at least one zone of the second space, by using this table of parameters to fill in a QCimage correspondence table.
Description
PROCEDE DE CONTROLE DE QUALITE D'UNEPROCESS FOR QUALITY CONTROL OF A
MIGRATION DE DONNEES SISMIQUES UTILISANT LA TRANSFORMEE SEISMIC DATA MIGRATION USING THE TRANSFORM
DE RADON GENERALISEEOF GENERALIZED RADON
DESCRIPTIONDESCRIPTION
Domaine technique La présente invention concerne un procédé de contrôle de qualité d'une migration de données sismiques utilisant la transformée de Radon généralisée Technical Field The present invention relates to a method for controlling the quality of a migration of seismic data using the generalized Radon transform
(GRT).(GRT).
Etat de la technique antérieure Il existe de nombreuses techniques utilisées par les géophysiciens pour connaître les couches STATE OF THE PRIOR ART There are many techniques used by geophysicists to know the layers
géologiques qui constituent le sous-sol. which constitute the subsoil.
Une de celles-ci, appelée " sismique de puits ", qui est illustrée sur la figure 1, consiste à émettre des signaux acoustiques à la surface du sol (source S) et à enregistrer les ondes sismiques transmises, après réflexion et/ou réfraction sur les limites entre les différentes couches géologiques, en One of these, called "well seismic", which is illustrated in FIG. 1, consists in emitting acoustic signals on the ground surface (source S) and in recording the seismic waves transmitted, after reflection and / or refraction on the boundaries between the different geological layers, in
utilisant des capteurs Ci disposés dans un puits 10. using sensors Ci placed in a well 10.
Après enregistrement des ondes par les capteurs, celles-ci sont alors filtrées pour déterminer quelles ondes ont été réfléchies par des réflecteurs situés sous les capteurs. La figure 2A représente des signaux enregistrés avant filtrage et la figure 2B les After recording the waves by the sensors, they are then filtered to determine which waves have been reflected by reflectors located under the sensors. FIG. 2A represents signals recorded before filtering and FIG. 2B represents them
signaux correspondants après filtrage. corresponding signals after filtering.
Un procédé de migration appliqué aux signaux filtrés permet alors d'obtenir une image migrée, que l'on peut assimiler à un plan de coupe vertical du sous-sol passant par la(ou les) source(s) et les capteurs dans le puits. Comme illustré sur la figure 3 l'image migrée est constituée par une succession A migration process applied to the filtered signals then makes it possible to obtain a migrated image, which can be likened to a vertical section plane of the basement passing through the source (s) and the sensors in the well. . As illustrated in Figure 3 the migrated image is constituted by a succession
SP 16505 DBSP 16505 DB
d'enregistrements verticaux, appelés " traces sismiques migrées ", constitués de signaux de plus ou moins forte vertical records, called "migrated seismic traces", made up of more or less strong signals
amplitude, composés d'arches positives et négatives. amplitude, composed of positive and negative arches.
Les signaux relatifs à un même réflecteur corrèlent d'une trace à l'autre, faisant apparaître à l'oeil des délinéations qui correspondent aux limites entre les The signals relating to the same reflector correlate from one trace to another, revealing to the eye delineations which correspond to the limits between the
couches géologiques.geological layers.
Le document référencé [1] en fin de description The document referenced [1] at the end of the description
décrit un exemple de procédé de migration qui utilise la transformée de Radon généralisée (GRT). Ce procédé de migration GRT construit l'image migrée point à point. Comme illustré sur la figure 4, chaque point de l'image, Il ou I2, est considéré comme un point de changement de vélocité du milieu, capable de réfléchir l'énergie provenant de la source S. L'énergie diffusée par ce point se propage dans toutes les directions pour être finalement détectée par les capteurs. La durée de parcours de l'onde sismique est calculée en divisant la somme de la distance is de la source à ce point et de la distance lc de ce point au capteur, par la vitesse describes an example of a migration process that uses the generalized Radon transform (GRT). This GRT migration process builds the point-to-point migrated image. As illustrated in figure 4, each point of the image, Il or I2, is considered as a point of change of velocity of the medium, capable of reflecting the energy coming from the source S. The energy diffused by this point propagates in all directions to be finally detected by the sensors. The travel time of the seismic wave is calculated by dividing the sum of the distance is from the source at this point and the distance lc from this point to the sensor, by the speed
de propagation supposée de l'onde dans le sous-sol. assumed propagation of the wave in the basement.
Cette vitesse de propagation est appelée vélocité du milieu. Pour chaque point réflecteur particulier dans l'espace de l'image migrée, on peut obtenir une courbe de temps de trajet, dans l'espace des données, CI1, CI2, en calculant les temps d'arrivée des ondes sismiques en chaque capteur, comme illustré sur la This propagation speed is called the velocity of the medium. For each particular reflecting point in the space of the migrated image, it is possible to obtain a journey time curve, in the data space, CI1, CI2, by calculating the arrival times of the seismic waves at each sensor, as shown in the
figure 5.figure 5.
Inversement à chaque point P de l'espace des données correspond l'amplitude du signal mesurée par un capteur à un instant donné, qui est due à plusieurs ondes, réfléchies par différents points réflecteurs, qui atteignent ce capteur en même temps. Ainsi, tout point de l'image migrée, dont la courbe de temps de Conversely, at each point P in the data space corresponds the amplitude of the signal measured by a sensor at a given instant, which is due to several waves, reflected by different reflecting points, which reach this sensor at the same time. Thus, any point in the migrated image, including the time curve of
SP 16505 DBSP 16505 DB
trajet passe par un point de donnée P, contribue à l'amplitude du signal en ce point. L'ensemble ce ces points réflecteurs constitue une courbe isochrone qui, comme illustré sur la figure 4, forme une ellipse, dont la source S et le capteur C sont les focales, lorsque path passes through a data point P, contributes to the amplitude of the signal at this point. The assembly of these reflecting points constitutes an isochronous curve which, as illustrated in FIG. 4, forms an ellipse, of which the source S and the sensor C are the focal points, when
la vélocité est constante dans le sous-sol. the velocity is constant in the basement.
En application de la transformée de Radon généralisée, chaque point de l'image migrée est obtenu en additionnant les contributions des données le long de la courbe de temps de trajet correspondante, ces contributions étant pondérées pour tenir compte des phénomènes physiques relatifs à la propagation des ondes, par exemple des pertes d'amplitude proportionnelles à la longueur du trajet. Pour calculer les courbes de trajet par tracé de rayon, on utilise un modèle de vitesse initial qui intègre les connaissances a priori du sous-sol, tel que celui illustré sur la In application of the generalized Radon transform, each point of the migrated image is obtained by adding the contributions of the data along the corresponding travel time curve, these contributions being weighted to take account of the physical phenomena relating to the propagation of waves, for example amplitude losses proportional to the length of the path. To calculate the path curves by ray tracing, we use an initial speed model which integrates the a priori knowledge of the subsoil, such as that illustrated on the
figure 6.figure 6.
Dans l'algorithme de migration GRT on a ainsi les étapes suivantes: Pour chaque trace T des données, pour chaque point (x,z) de l'image migrée: on calcule le temps de trajet ts de la source S au point (x,z); - on calcule le temps de trajet tc du capteur C au point (x, z); In the GRT migration algorithm we have the following steps: For each trace T of the data, for each point (x, z) of the migrated image: we calculate the travel time ts from the source S to the point (x , z); - the journey time tc of the sensor C at the point (x, z) is calculated;
- on a alors le temps de trajet source- - we then have the source journey time -
capteur: tsc = ts + tc; - on cherche l'échantillon u de la trace T au temps tsc; - on cherche le facteur de pondération w relatif à la source S, au capteur C et au point (x,z); - on calcule la valeur d'une table Image (x,z) sensor: tsc = ts + tc; - we look for the sample u of the trace T at time tsc; - looking for the weighting factor w relative to the source S, to the sensor C and to the point (x, z); - we calculate the value of an Image table (x, z)
en ce point: Image (x,z) = Image (x,z) + u x w. at this point: Image (x, z) = Image (x, z) + u x w.
SP 16505 DBSP 16505 DB
Une variante de cet algorithme GRT consiste, dans une première étape à calculer la table des paramètres temps de trajet t, longueur de trajet 1, orientation du rayon à ses extrémités d et a; et dans une seconde étape, à effectuer le calcul de migration. Ainsi comme illustré sur la figure 7, pour chaque source et pour chaque capteur C, on calcule la durée du temps de trajet t(P,Ixy) pour aller de ladite source ou dudit capteur vers tous les points de l'image migrée Ixy. On calcule par ailleurs la longueur du trajet l(P,Ixy) ainsi que l'orientation du rayon à ses extrémités a(P,Ixy) et d(P, Ixy). Dans la seconde étape la migration utilise la table calculée pendant la A variant of this GRT algorithm consists, in a first step in calculating the table of parameters journey time t, path length 1, orientation of the radius at its ends d and a; and in a second step, to perform the migration calculation. Thus, as illustrated in FIG. 7, for each source and for each sensor C, the duration of the travel time t (P, Ixy) to go from said source or said sensor to all the points of the migrated image Ixy is calculated. The length of the path l (P, Ixy) is also calculated as well as the orientation of the radius at its ends a (P, Ixy) and d (P, Ixy). In the second step the migration uses the table calculated during the
première étape.first stage.
Avec par exemple des signaux sismiques filtrés, tels que représentés sur la figure 8, la migration GRT permet d'obtenir une image migrée telle que représentée With, for example, filtered seismic signals, as represented in FIG. 8, the GRT migration makes it possible to obtain a migrated image as represented
sur la figure 9.in figure 9.
L'opérateur, avant d'utiliser les résultats ainsi représentés en figure 9, doit en vérifier le bien The operator, before using the results thus represented in figure 9, must verify the good
fondé. Un contrôle de qualité est donc nécessaire. based. Quality control is therefore necessary.
Celui-ci est basé traditionnellement sur des critères plus ou moins subjectifs tels que la correspondance des résultats avec les connaissances a priori du sous-sol, ou la continuité des délinéations représentant les réflecteurs. Il est donc très difficile à effectuer, d'autant plus que les représentations des données avant la migration (voir figure 8), et des données dans This is traditionally based on more or less subjective criteria such as the correspondence of the results with a priori knowledge of the subsoil, or the continuity of the delineations representing the reflectors. It is therefore very difficult to perform, especially since the representations of the data before migration (see Figure 8), and of data in
l'image migrée (voir figure 9) sont très différentes. the migrated image (see Figure 9) are very different.
Une méthode décrite dans le document [2], pour la migration de Kirchhoff, permet d'associer à chaque point de l'image migrée un point dans l'espace des données. Cette méthode, toutefois, ne permet pas d'associer à chaque point de l'image migrée l'ensemble des points contributeurs le long de la courbe de temps A method described in document [2], for the Kirchhoff migration, makes it possible to associate with each point of the migrated image a point in the data space. This method, however, does not make it possible to associate with each point of the migrated image all of the contributing points along the time curve
SP 16505 DBSP 16505 DB
de trajet. Elle ne permet pas non plus de trouver l'isochrone associée à un point dans l'espace des données. La présente invention a pour objectif de résoudre ces problèmes en proposant un procédé de contrôle de qualité pour une migration GRT de données sismiques, par une identification interactive des événements entre l'image obtenue par migration et les path. It also does not find the isochronous associated with a point in the data space. The present invention aims to solve these problems by proposing a quality control method for a GRT migration of seismic data, by an interactive identification of the events between the image obtained by migration and the
données sismiques avant la migration. seismic data before migration.
Exposé de l'invention La présente invention propose un procédé de contrôle de qualité d'une migration de données sismiques utilisant la transformée de Radon généralisée, permettant de passer d'un espace de données sismiques à un espace de l'image migrée, dans lequel on calcule une table de paramètres permettant de connaître pour une onde sismique allant d'un point de l'image à une source ou à un capteur, la longueur du trajet, sa durée et les angles d'incidence de l'onde au départ et à l'arrivée du trajet, caractérisé en ce que l'on établit une correspondance entre au moins une zone d'un premier de ces deux espaces et au moins une zone du second espace, en utilisant cette table des paramètres pour remplir une table de correspondance QCimage. Cette table des paramètres est obtenue habituellement par tracé de rayons, mais elle peut aussi être obtenue par tout autre moyen permettant de connaître, pour une onde sismique allant d'un point de l'image à une source ou à un capteur, la longueur du trajet, sa durée et les angles d'incidence de l'onde au SUMMARY OF THE INVENTION The present invention proposes a method for controlling the quality of a migration of seismic data using the generalized Radon transform, making it possible to pass from a space of seismic data to a space of the migrated image, in which a table of parameters is calculated making it possible to know, for a seismic wave going from a point of the image to a source or to a sensor, the length of the path, its duration and the angles of incidence of the wave at the start and at the end of the journey, characterized in that a correspondence is established between at least one zone of a first of these two spaces and at least one zone of the second space, using this table of parameters to fill a table of QCimage correspondence. This table of parameters is usually obtained by ray tracing, but it can also be obtained by any other means making it possible to know, for a seismic wave going from a point of the image to a source or to a sensor, the length of the path, its duration and the angles of incidence of the wave at
départ et à l'arrivée du trajet.departure and arrival of the journey.
Le procédé de l'invention permet, d'une part, de pointer une zone dans l'espace des données sismiques The method of the invention allows, on the one hand, to point an area in the space of seismic data
SP 16505 DBSP 16505 DB
et de la projeter dans l'espace de l'image migrée par une migration locale relative à la zone pointée (" projection avant "). Ce procédé permet, d'autre part, de pointer une zone dans l'espace de l'image migrée et de la projeter dans l'espace des données sismiques, par le calcul du champ des données sismiques contribuant à la zone pointée dans l'image migrée and to project it into the space of the migrated image by a local migration relative to the pointed area ("front projection"). This method makes it possible, on the other hand, to point an area in the space of the migrated image and to project it into the space of the seismic data, by calculating the field of the seismic data contributing to the area pointed in the migrated image
(" projection arrière ").("rear projection").
Dans le premier mode de réalisation (projection avant) le point (x, z) de l'image migrée est mis en correspondance avec la zone correspondante dans l'espace des données si la valeur QCimage (x,z) de la table de correspondance est supérieure à un seuil donné, par exemple zéro. Avantageusement la zone pointée dans l'espace des données est un ensemble d'intervalles Zi, caractérisés chacun par une trace T, In the first embodiment (front projection) the point (x, z) of the migrated image is mapped to the corresponding zone in the data space if the QCimage (x, z) value of the correspondence table is greater than a given threshold, for example zero. Advantageously, the area pointed in the data space is a set of intervals Zi, each characterized by a trace T,
le début td et la fin tf de l'intervalle correspondant. the beginning td and the end tf of the corresponding interval.
Pour chaque intervalle Zi et pour chaque point (x,z) de l'image migrée, on réalise les étapes suivantes: - pour le rayon de la source S, relative à cet intervalle Zi, au point (x,z), on lit dans la table des paramètres les paramètres temps de trajet ts, longueur du trajet ls, orientation du rayon à ses extrémités ds et as; - pour le rayon du capteur C, relatif à cet intervalle Zi, au point (x,z), on lit dans la table des paramètres les paramètres temps de trajet tc, longueur du trajet lc, orientation du rayon à ses extrémités dc et ac; - on calcule le temps de trajet source-capteur tsc = tc + ts; - si td < tsc < tf, on calcule le facteur de pondération w(ls, ds, as, lc, dc, ac); - on calcule la valeur QCimage (x,z) = QCimage(x,z)+w For each interval Zi and for each point (x, z) of the migrated image, the following steps are carried out: - for the radius of the source S, relative to this interval Zi, at point (x, z), we read in the parameter table the parameters travel time ts, length of the journey ls, orientation of the ray at its ends ds and as; - for the radius of sensor C, relative to this interval Zi, at point (x, z), we read in the parameters table the parameters travel time tc, length of the journey lc, orientation of the radius at its ends dc and ac ; - the source-sensor travel time tsc = tc + ts is calculated; - if td <tsc <tf, we calculate the weighting factor w (ls, ds, as, lc, dc, ac); - we calculate the value QCimage (x, z) = QCimage (x, z) + w
SP 16505 DBSP 16505 DB
- on effectue des boucles sur les points et sur les intervalles pour remplir la table QCimage; - on utilise cette table QCimage pour établir - loops are performed on the points and on the intervals to fill the QCimage table; - we use this QCimage table to establish
la correspondance entre zones.correspondence between zones.
Dans le second mode de réalisation (projection arrière) le point correspondant à une trace T au temps de trajet source-capteur tsc des données est mis en correspondance avec la zone correspondante dans l'image migrée si la valeur QCimage (T, tsc) de la table de correspondance est supérieure à un seuil donné, par exemple zéro. Avantageusement la zone pointée dans l'image migrée est un ensemble d'intervalles Zi, caractérisés par la position de la trace xz, le début zd et la fin zf de l'intervalle correspondant. Pour chaque trace T des données d'entrée, pour chaque intervalle Zi, pour chaque point (x,z) de l'intervalle Zi, on a les étapes suivantes: - pour le rayon de la source S, relative à la trace T, au point (x,z), on lit dans la table de paramètres les paramètres temps de trajet ts, longueur de trajet ls, orientation du rayon à ses extrémités ds et as; - pour le rayon du capteur C, relatif à la trace T, au point (x,z), on lit dans la table de paramètres les paramètres temps de trajet tc, longueur de trajet lc, orientation du rayon à ses extrémités dc et ac; - on calcule le temps de trajet source-capteur tsc = tc + ts; - on calcule le facteur de pondération w(ls, ds, as, ls, dc, ac); - on calcule la valeur QCimage(T,tsc) = QCimage(T,tsc)+w In the second embodiment (rear projection) the point corresponding to a trace T at the source-sensor travel time tsc of the data is mapped to the corresponding zone in the migrated image if the value QCimage (T, tsc) of the correspondence table is greater than a given threshold, for example zero. Advantageously, the area pointed in the migrated image is a set of intervals Zi, characterized by the position of the trace xz, the beginning zd and the end zf of the corresponding interval. For each trace T of the input data, for each interval Zi, for each point (x, z) of the interval Zi, we have the following steps: - for the radius of the source S, relative to the trace T, at point (x, z), the parameters travel time ts, path length ls, orientation of the ray at its ends ds and as are read in the parameter table; - for the radius of the sensor C, relative to the trace T, at the point (x, z), the parameters travel time tc, path length lc, orientation of the radius at its ends dc and ac are read in the parameter table ; - the source-sensor travel time tsc = tc + ts is calculated; - we calculate the weighting factor w (ls, ds, as, ls, dc, ac); - we calculate the value QCimage (T, tsc) = QCimage (T, tsc) + w
SP 16505 DBSP 16505 DB
- on effectue des boucles sur les points, sur les intervalles, et sur les traces pour remplir la - loops are carried out on the points, on the intervals, and on the traces to fill the
table QCimage;QCimage table;
- on utilise cette table QCimage pour établir la correspondance entre zones. Avantageusement, dans ledit procédé, on affecte un paramètre, dont la valeur est proportionnelle à la valeur de QCimage, aux zones correspondantes dans les deux espaces. Ce paramètre peut, par exemple, être une couleur semi-transparente, - we use this QCimage table to establish the correspondence between zones. Advantageously, in said method, a parameter is assigned, the value of which is proportional to the value of QCimage, to the corresponding zones in the two spaces. This parameter can, for example, be a semi-transparent color,
et sa valeur la densité de cette couleur. and its value the density of this color.
Le procédé de l'invention est utilisable en de nombreux domaines, notamment en sismique de puits, et également en sismique de surface " pré-stack ", les traces provenant de diverses sources et capteurs n'étant pas ajoutées les unes aux autres dans ce cas, The method of the invention can be used in many fields, in particular in well seismic, and also in "pre-stack" surface seismic, the traces coming from various sources and sensors not being added to each other in this case,
ce qui permet d'appliquer le procédé de migration GRT. which allows to apply the GRT migration process.
Brève description des dessinsBrief description of the drawings
- La figure 1 illustre la sismique de puits; - la figure 2A et 2B illustrent les données avant la migration, avant et après filtrage; - la figure 3 illustre l'image migrée; - les figures 4 et 5 illustrent la construction de l'image migrée point par point, et les courbes de temps de trajet correspondantes; - la figure 6 illustre un exemple de modèle de vitesse; - la figure 7 illustre une variante de l'algorithme traditionnel de migration GRT; - la figure 8 illustre un exemple de données sismiques obtenues après filtrage; - la figure 9 illustre une image migrée obtenue par migration GRT, à partir des données de la figure 8; - Figure 1 illustrates the well seismic; - Figure 2A and 2B illustrate the data before migration, before and after filtering; - Figure 3 illustrates the migrated image; - Figures 4 and 5 illustrate the construction of the migrated image point by point, and the corresponding travel time curves; - Figure 6 illustrates an example of a speed model; - Figure 7 illustrates a variant of the traditional GRT migration algorithm; - Figure 8 illustrates an example of seismic data obtained after filtering; - Figure 9 illustrates a migrated image obtained by GRT migration, from the data in Figure 8;
SP 16505 DBSP 16505 DB
- la figure 10 illustre le flot des données dans le procédé de l'invention; - la figure 11 illustre des zones correspondantes dans l'espace des données et l'espace de l'image selon le procédé de l'invention; - les figures 12 et 14 illustrent des organigrammes de fonctionnement du procédé de l'invention respectivement pour une projection avant et pour une projection arrière; - la figure 13 illustre le fonctionnement du - Figure 10 illustrates the flow of data in the method of the invention; - Figure 11 illustrates corresponding areas in the data space and the image space according to the method of the invention; - Figures 12 and 14 illustrate flowcharts of operation of the method of the invention respectively for a front projection and for a rear projection; - Figure 13 illustrates the operation of the
procédé de l'invention pour une projection avant. method of the invention for front projection.
Exposé détaillé de modes de réalisation La présente invention concerne un procédé de contrôle de qualité d'une migration GRT de données sismiques, par exemple en sismique de puits, qui permet de passer d'un espace de données sismiques à un espace Detailed description of embodiments The present invention relates to a method for controlling the quality of a GRT migration of seismic data, for example in well seismic, which makes it possible to pass from a space of seismic data to a space
de l'image migrée.of the migrated image.
Lors de la migration GRT, la table des paramètres obtenus par tracé de rayon est sauvegardée sur fichier. Ces résultats intermédiaires, pour un rayon allant d'un point de l'image à une source ou un capteur, sont la longueur du trajet, sa durée, et les During GRT migration, the table of parameters obtained by radius tracing is saved to file. These intermediate results, for a radius going from a point of the image to a source or a sensor, are the length of the path, its duration, and the
angles d'incidence au départ et à l'arrivée du trajet. angles of incidence at the start and finish of the journey.
Cette table peut donc être utilisée dans les étapes du procédé de l'invention. Ainsi, comme illustré sur la figure 10, en sortie des capteurs (acquisition) on obtient les données sismiques brutes qui après filtrage sont soumises au procédé de migration GRT, le procédé de l'invention utilisant les données sismiques This table can therefore be used in the steps of the method of the invention. Thus, as illustrated in FIG. 10, at the output of the sensors (acquisition), the raw seismic data are obtained which after filtering are subjected to the GRT migration method, the method of the invention using the seismic data
filtrées, l'image migrée et la table des paramètres. filtered, the migrated image and the parameter table.
Dans le procédé de l'invention on établit une correspondance entre au moins une zone d'un des deux espaces, espace de données sismiques ou espace de In the method of the invention a correspondence is established between at least one zone of one of the two spaces, seismic data space or space of
l'image migrée, et au moins une zone de l'autre espace. the migrated image, and at least one area of the other space.
SP 16505 DBSP 16505 DB
Comme illustré sur la figure 11, on peut ainsi visualiser des zones correspondantes dans l'espace des données et dans l'espace de l'image migrée, lesdites zones étant pointées interactivement, soit dans l'espace des données, soit dans l'espace de l'image migrée. La projection d'une zone à partir de l'espace des données sur l'espace de l'image migrée est appelée " projection avant ". Inversement, la projection d'une zone à partir de l'image migrée sur l'espace des As illustrated in FIG. 11, it is thus possible to visualize corresponding zones in the data space and in the space of the migrated image, said zones being pointed interactively, either in the data space or in space of the migrated image. The projection of an area from the data space onto the space of the migrated image is called "front projection". Conversely, the projection of an area from the migrated image onto the space of
données est appelée " projection arrière ". data is called "rear projection".
Dans le premier mode de fonctionnement (projection avant) les différentes étapes du procédé de contrôle de qualité, selon l'invention, sont illustrées In the first operating mode (front projection) the different stages of the quality control process according to the invention are illustrated
sur l'organigramme de la figure 12.on the flowchart of figure 12.
La zone pointée est tout d'abord décomposée en NZ intervalles Zi caractérisés par une trace T, un début d'intervalle td et une fin d'intervalle tf, comme The pointed area is first of all decomposed into NZ intervals Zi characterized by a trace T, a beginning of interval td and an end of interval tf, as
illustré sur la figure 13.illustrated in figure 13.
Pour chaque intervalle Zi et pour chaque point (x,z) de l'image migrée, on a les étapes suivantes: - lecture des paramètres dans la table des paramètres pour le rayon de la source S, relative à cet intervalle Zi, au point (x,z), à savoir: durée de temps de trajet ts, longueur du trajet is, orientation du rayon à ses extrémités as et ds; - lecture des paramètres dans la table des paramètres pour le rayon du capteur C, relatif à cet intervalle Zi, au point (x,z) à savoir: durée du temps de trajet tc, longueur du trajet lc, orientation du rayon à ses extrémités ac et dc; calcul du temps de trajet source-capteur: ts = tc + ts; - si ce temps tsc est compris à l'intérieur de l'intervalle Zi défini plus haut, c'est-à-dire: For each interval Zi and for each point (x, z) of the migrated image, there are the following steps: - reading of the parameters in the parameter table for the radius of the source S, relative to this interval Zi, at the point (x, z), namely: duration of travel time ts, length of travel is, orientation of the radius at its ends as and ds; - reading of the parameters in the parameter table for the radius of sensor C, relative to this interval Zi, at point (x, z), namely: duration of the travel time tc, length of the journey lc, orientation of the radius at its ends ac and dc; calculation of the source-sensor travel time: ts = tc + ts; - if this time tsc is included within the interval Zi defined above, that is to say:
SP 16505 DBSP 16505 DB
td < tsc < tf, calcul du facteur de pondération w qui est fonction des paramètres is, ds, as, lc, dc et ac, ainsi que la valeur d'une table Qcimage en ce point QCimage(x,z) = QCimage(x,z) + w; - boucles sur les points et sur les intervalles td <tsc <tf, calculation of the weighting factor w which is a function of the parameters is, ds, as, lc, dc and ac, as well as the value of a table Qcimage at this point QCimage (x, z) = QCimage ( x, z) + w; - loops on the points and on the intervals
pour remplir cette table QCimage.to fill this QCimage table.
Cette table QCimage est ensuite utilisée pour pointer l'image migrée. Le point (x,z) de l'image migrée est mis en correspondance avec la zone pointée dans l'espace des données si la valeur QCimage(x,z) est This QCimage table is then used to point to the migrated image. The point (x, z) of the migrated image is mapped to the area pointed in the data space if the value QCimage (x, z) is
supérieure à un seuil donné, par exemple 0. greater than a given threshold, for example 0.
On peut alors affecter aux zones en correspondance un même paramètre dont la valeur est We can then assign to the corresponding zones the same parameter whose value is
proportionnelle à la valeur QCimage (x,z). proportional to the QCimage value (x, z).
Dans le second mode de fonctionnement (projection arrière) les différentes étapes du procédé de contrôle de qualité selon l'invention sont In the second operating mode (rear projection) the different stages of the quality control process according to the invention are
illustrées sur l'organigramme de la figure 14. illustrated in the flowchart in Figure 14.
La zone pointée dans l'image migrée est tout d'abord décomposée en NZ intervalles Zi caractérisés par une trace (x,z), un début d'intervalle zd et une The area pointed in the migrated image is first broken down into NZ intervals Zi characterized by a trace (x, z), a start of interval zd and a
fin d'intervalle zf.end of interval zf.
Pour chaque trace T des données d'entrée, pour chaque intervalle Zi et pour chaque point (x,z) de l'intervalle Zi de l'image migrée on a les étapes suivantes: - lecture dans la table des paramètre des paramètres pour le rayon de la source S, relative à la trace T, au point (x,z), à savoir: durée de temps de trajet ts, longueur du trajet ls, orientation du rayon à ses extrémités as et ds; - lecture dans la table des paramètres des paramètres pour le rayon du capteur C, relatif à la trace T, au point (x,z) à savoir: durée du temps de For each trace T of the input data, for each interval Zi and for each point (x, z) of the interval Zi of the migrated image, there are the following steps: - reading in the parameter table of the parameters for the radius of the source S, relative to the trace T, at the point (x, z), namely: duration of travel time ts, length of the travel ls, orientation of the radius at its ends as and ds; - reading in the parameter parameters table for the radius of the sensor C, relative to the trace T, at the point (x, z) namely: duration of the time
SP 16505 DBSP 16505 DB
trajet tc, longueur du trajet lc, orientation du rayon à ses extrémités ac et dc; - calcul du temps de trajet source-capteur: tsc = tc + ts; - calcul du facteur de pondération w qui est fonction des paramètres ls, ds, as, lc, dc et ac, ainsi que la valeur d'une table QCimage en ce point QCimage(T,tsc) = QCimage(T,tsc) + w; - boucles sur les points, sur les intervalles Zi et sur les traces T pour remplir cette table QCimage. La table QCimage est ensuite utilisée pour pointer les données d'entrée. Le point correspondant à la trace T au temps tsc des données est mis en correspondance avec la zone pointée dans l'image migrée si la valeur QCimage( T,tsc) est supérieure à un seuil path tc, length of path lc, orientation of the ray at its ends ac and dc; - calculation of the source-sensor travel time: tsc = tc + ts; - calculation of the weighting factor w which is a function of the parameters ls, ds, as, lc, dc and ac, as well as the value of a QCimage table at this point QCimage (T, tsc) = QCimage (T, tsc) + w; - loops on the points, on the intervals Zi and on the traces T to fill this QCimage table. The QCimage table is then used to point to the input data. The point corresponding to the trace T at the time tsc of the data is matched with the area pointed in the migrated image if the QCimage value (T, tsc) is greater than a threshold
donné, par exemple 0.given, for example 0.
On peut alors affecter aux zones en correspondances un paramètre dont la valeur est We can then assign to the correspondence zones a parameter whose value is
proportionnelle à la valeur de QCimage(T,tsc). proportional to the value of QCimage (T, tsc).
Dans ses deux modes de réalisation le procédé de l'invention permet ainsi d'identifier les réflecteurs, et vérifier la qualité de la migration et In its two embodiments, the method of the invention thus makes it possible to identify the reflectors, and to verify the quality of the migration and
des données en entrée de la migration. data as input to the migration.
Dans un mode de réalisation avantageux les zones correspondantes dans l'espace des données et dans l'espace de l'image migrée peuvent être représentées par des taches de couleurs semi- transparentes superposées sur l'image migrée et les données sismiques avant migration. Après le calcul des tables par tracé de rayon et la migration, l'utilisateur effectue un pointé interactif de zones, soit sur les données d'entrée (projection avant), soit sur l'image migrée (projection arrière), la zone pointée étant alors colorée par une couleur définie par l'utilisateur. On In an advantageous embodiment, the corresponding zones in the data space and in the space of the migrated image can be represented by spots of semi-transparent colors superimposed on the migrated image and the seismic data before migration. After calculating the tables by radius tracing and migration, the user performs an interactive pointing of areas, either on the input data (front projection) or on the migrated image (rear projection), the pointed area being then colored by a user-defined color. We
SP 16505 DBSP 16505 DB
peut alors, avantageusement, colorer l'image migrée (en projection avant) ou les données (en projection arrière) avec une valeur dont la densité est can then, advantageously, color the migrated image (in front projection) or the data (in rear projection) with a value whose density is
proportionnelle à la valeur de QCimage. proportional to the value of QCimage.
SP 16505 DBSP 16505 DB
REFERENCESREFERENCES
[1] " The Generalized Radon Transform, A Breakthrough In Seismic Migration " de M Oristaglio, G. Beylkin et D. Miller (Seismics, The Technical Review, Volume 35, numéro 3, 1987, pages 20 à 27) [2] " Automatic Association Of Kinematic Information To Prestack Images " de S. Geoltrain et E. Chovet (61st Annual Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 1991, pages 890-892) [1] "The Generalized Radon Transform, A Breakthrough In Seismic Migration" by M Oristaglio, G. Beylkin and D. Miller (Seismics, The Technical Review, Volume 35, number 3, 1987, pages 20 to 27) [2] " Automatic Association Of Kinematic Information To Prestack Images "by S. Geoltrain and E. Chovet (61st Annual Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 1991, pages 890-892)
SP 16505 DBSP 16505 DB
Claims (11)
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR9903590A FR2791438B1 (en) | 1999-03-23 | 1999-03-23 | QUALITY CONTROL METHOD OF SEISMIC DATA MIGRATION USING GENERALIZED RADON TRANSFORM |
| US09/507,878 US6625544B1 (en) | 1999-03-23 | 2000-02-22 | Method of controlling the quality of seismic data migration by using the generalized radon transform description |
| GB0004915A GB2348285B (en) | 1999-03-23 | 2000-03-02 | Quality control process for the migration of seismic data using the generalized radon transform |
| NO20001475A NO20001475L (en) | 1999-03-23 | 2000-03-22 | Procedure for controlling the quality of seismic data migration using the generalized Radon transform description |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR9903590A FR2791438B1 (en) | 1999-03-23 | 1999-03-23 | QUALITY CONTROL METHOD OF SEISMIC DATA MIGRATION USING GENERALIZED RADON TRANSFORM |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| FR2791438A1 true FR2791438A1 (en) | 2000-09-29 |
| FR2791438B1 FR2791438B1 (en) | 2001-08-24 |
Family
ID=9543524
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR9903590A Expired - Fee Related FR2791438B1 (en) | 1999-03-23 | 1999-03-23 | QUALITY CONTROL METHOD OF SEISMIC DATA MIGRATION USING GENERALIZED RADON TRANSFORM |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US6625544B1 (en) |
| FR (1) | FR2791438B1 (en) |
| GB (1) | GB2348285B (en) |
| NO (1) | NO20001475L (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2818387B1 (en) * | 2000-12-18 | 2003-02-14 | Inst Francais Du Petrole | METHOD FOR OBTAINING REFLECTION TRAVEL TIMES FROM AN INTERPRETATION OF SEISMIC DATA IN MIGRATED CYLINDRIC WAVES |
| US7542373B2 (en) * | 2003-10-28 | 2009-06-02 | Baker Hughes Incorporated | Vector 3-component 3-dimensional kirchhoff prestack migration |
| WO2009082545A1 (en) * | 2007-12-21 | 2009-07-02 | Exxonmobil Upstream Research Company | Detection of features in seismic images |
| US9542507B2 (en) | 2012-02-23 | 2017-01-10 | Schlumberger Technology Corporation | Feature detection in seismic volumes |
| CN111913215B (en) * | 2019-05-10 | 2023-04-25 | 中国石油天然气集团有限公司 | Backscatter amplitude-preserving offset imaging method, apparatus and computer storage medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4760563A (en) * | 1986-01-09 | 1988-07-26 | Schlumberger Technology Corporation | Seismic exploration using exactly invertible discrete transformation into tau-p space |
| EP0515189A2 (en) * | 1991-05-22 | 1992-11-25 | Western Atlas International, Inc. | Noise attenuation method in seismic signals |
| FR2739696A1 (en) * | 1995-10-06 | 1997-04-11 | Elf Aquitaine | TREATMENT METHOD OF RECORDED SEISMIC REFLECTION TRACKS FOR VARIABLE DEPORTS |
| US5677893A (en) * | 1994-07-07 | 1997-10-14 | Schlumberger Technology Corporation | Method of processing seismic data |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5583825A (en) * | 1994-09-02 | 1996-12-10 | Exxon Production Research Company | Method for deriving reservoir lithology and fluid content from pre-stack inversion of seismic data |
| US6446007B1 (en) * | 1998-11-05 | 2002-09-03 | Exxonmobil Upstream Research Company | Method for controlled-amplitude prestack time migration of seismic data |
| US6374185B1 (en) * | 2000-02-18 | 2002-04-16 | Rdsp I, L.P. | Method for generating an estimate of lithological characteristics of a region of the earth's subsurface |
-
1999
- 1999-03-23 FR FR9903590A patent/FR2791438B1/en not_active Expired - Fee Related
-
2000
- 2000-02-22 US US09/507,878 patent/US6625544B1/en not_active Expired - Fee Related
- 2000-03-02 GB GB0004915A patent/GB2348285B/en not_active Expired - Fee Related
- 2000-03-22 NO NO20001475A patent/NO20001475L/en not_active Application Discontinuation
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4760563A (en) * | 1986-01-09 | 1988-07-26 | Schlumberger Technology Corporation | Seismic exploration using exactly invertible discrete transformation into tau-p space |
| EP0515189A2 (en) * | 1991-05-22 | 1992-11-25 | Western Atlas International, Inc. | Noise attenuation method in seismic signals |
| US5677893A (en) * | 1994-07-07 | 1997-10-14 | Schlumberger Technology Corporation | Method of processing seismic data |
| FR2739696A1 (en) * | 1995-10-06 | 1997-04-11 | Elf Aquitaine | TREATMENT METHOD OF RECORDED SEISMIC REFLECTION TRACKS FOR VARIABLE DEPORTS |
Non-Patent Citations (2)
| Title |
|---|
| GEOLTRAIN S ET AL: "automatic association of kinematic information to prestack images", ANNUAL INTERNATIONAL SEG MEETING 61TH, EXPANDED ABSTRACTS,US,TULSA, OK, vol. 2, pages 890-892, XP000863152 * |
| MILLER D ET AL: "A new slant on seismic imaging: migration and integral geometry", GEOPHYSICS, JULY 1987, USA, vol. 52, no. 7, pages 943 - 964, XP002124814, ISSN: 0016-8033 * |
Also Published As
| Publication number | Publication date |
|---|---|
| GB0004915D0 (en) | 2000-04-19 |
| US6625544B1 (en) | 2003-09-23 |
| GB2348285B (en) | 2001-07-04 |
| NO20001475L (en) | 2000-09-25 |
| GB2348285A (en) | 2000-09-27 |
| FR2791438B1 (en) | 2001-08-24 |
| NO20001475D0 (en) | 2000-03-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP1396712B1 (en) | Method and system for polarisation optical time domain reflectometry (POTDR) | |
| CA2514112A1 (en) | Method for using data expressed in differing time scales to construct a model of a heterogeneous environment described by several parameters | |
| EP0734536B1 (en) | Method for analyzing and processing reflection seismic data | |
| FR2465213A1 (en) | APPARATUS FOR DIGITAL COLORING OR COLOR MODIFICATION OF AN OBJECT | |
| EP0663590B1 (en) | Modulated spectroscopic ellipsometer | |
| CN102362146B (en) | Reflectance distribution curve modeling method, thickness measurement scheme and thickness measurement reflectometer using same | |
| EP0574987B1 (en) | Arrangement for processing a measurement signal corresponding to the intensity of X-rays reflected by a multilayer structure on a substrate | |
| EP0797780B1 (en) | Method for processing reflection seismic traces recorded for variable offsets | |
| WO2015067864A1 (en) | Method and device for processing seismic signals | |
| FR2791438A1 (en) | QUALITY CONTROL METHOD OF SEISMIC DATA MIGRATION USING GENERALIZED RADON TRANSFORM | |
| CA2464799C (en) | Method of determining a model for seismic wave velocity in a heterogeneous subsurface formation | |
| EP1797441A1 (en) | Method and installation for analyzing an integrated circuit | |
| FR2729766A1 (en) | METHOD FOR EXPLORING A COMPLEX TECTONIC ENVIRONMENT | |
| EP0374043B1 (en) | Process for analysing optical components, optical fibres or networks of optical guides by time-domain reflectometry, and time-domain reflectometer | |
| EP0944847B1 (en) | Method for seismic processing and in particular for three-dimensional seismic exploration using seismic data migration | |
| WO1998019180A1 (en) | Improved method for migration before summation | |
| FR2898971A1 (en) | Object`s e.g. polymer film, opto-geometric characteristic e.g. thickness, measuring method for e.g. packaging industry, involves determining thickness based on refractive index or index based on thickness from Fourier transform of fringes | |
| FR2858064A1 (en) | Determining parameters of speed and anellipticity, for treating seismic traces from collection at common midpoint including anelliptic obliquity correction | |
| FR2801387A1 (en) | Processing seismic data, when prospecting for oil or gas, comprises identifying maxima in energy matrix lying along path in space produced by parabolic Radon transformation, calculating energy matrix along selected path through Radon space | |
| FR2521298A1 (en) | SET OF ULTRASONIC CONTROL HEADS | |
| CA2452206C (en) | Improvements to surface pointing techniques | |
| LU82230A1 (en) | SURFACE QUALITY CONTROL PROCESS | |
| EP4143505B1 (en) | Method for obtaining the profile of a surface moving in relation to the system | |
| WO2000062013A1 (en) | Device for measuring the dimension and controlling of defects in optical fibres during production | |
| EP0571295A1 (en) | Device for evaluating a lubricant, especially of a car engine |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| ST | Notification of lapse |
Effective date: 20101130 |