FR2853085A1 - Une methode de prospection sismique - Google Patents
Une methode de prospection sismique Download PDFInfo
- Publication number
- FR2853085A1 FR2853085A1 FR0304298A FR0304298A FR2853085A1 FR 2853085 A1 FR2853085 A1 FR 2853085A1 FR 0304298 A FR0304298 A FR 0304298A FR 0304298 A FR0304298 A FR 0304298A FR 2853085 A1 FR2853085 A1 FR 2853085A1
- Authority
- FR
- France
- Prior art keywords
- group
- vibrators
- time
- vibrator
- scan
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/003—Seismic data acquisition in general, e.g. survey design
- G01V1/005—Seismic data acquisition in general, e.g. survey design with exploration systems emitting special signals, e.g. frequency swept signals, pulse sequences or slip sweep arrangements
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)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Transducers For Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Une méthode de prospection sismique comprenant les étapes d'activation du ou de chaque vibreur dans un premier groupe de vibreurs au temps T0, et d'activation subséquente du ou de chaque vibreur dans un deuxième groupe de vibreurs au temps T1 qui satisfait T0 < T1 < T0 + S1 + L où S1 est le temps de balayage du premier groupe de vibreurs et L est le temps d'écoute. Au moins un groupe parmi le premier groupe de vibreurs et le deuxième groupe de vibreurs comprend au moins deux vibreurs. Le premier groupe et le deuxième groupe de vibreurs peuvent être le même groupe ou bien ils peuvent être des groupes différents. Cette méthode permet de réduire le temps requis pour achever une étude sismique par comparaison avec les techniques de "tirs simultanés" et de "tirs à balayage glissant" de l'art antérieur.Dans un cas où le premier groupe et le deuxième groupe de vibreurs sont différents, la méthode peut en outre comprendre l'activation du ou de chaque vibreur dans le premier groupe de vibreurs au temps T2, où T1 < T2 < T1 + S2 + L et S2 est le temps de balayage du premier groupe de vibreurs, puis l'activation du ou de chaque vibreur dans le deuxième groupe de vibreurs au temps T3 où T2 < T3 < T2 + S1 + L et où T3 - T2 ≠ T1 - T0. Le temps de retard variable entre un tir du premier groupe de vibreurs et le tir correspondant du deuxième groupe de vibreurs signifie que du bruit harmonique se produira à des temps différents dans les enregistrements de tirs de sorte que le bruit puisse être éliminé en combinant de façon appropriée les enregistrements de tirs.
Description
Une méthode de prospection sismique
La présente invention concerne une méthode de prospection sismique.
En particulier, elle concerne une méthode de prospection sismique dans laquelle deux groupes de sources de vibration ou plus émettent de l'énergie sismique de telle manière que leurs temps de balayage (sweep) se chevauchent les uns les autres.
Le principe de la prospection sismique consiste à faire produire de 10 l'énergie sismique par une source d'énergie sismique, énergie qui se propage vers le bas à travers la Terre. L'énergie sismique se propageant vers le bas est réfléchie par une ou plusieurs structures géologiques à l'intérieur de la Terre qui font office de réflecteurs partiels d'énergie sismique. L'énergie sismique réfléchie est détectée par un ou plusieurs 15 capteurs (appelés généralement "récepteurs"). Il est possible d'obtenir des informations sur la structure géologique de la Terre à partir d'énergie sismique qui subit une réflexion au sein de la Terre et est par la suite acquise au niveau des récepteurs.
En pratique, un agencement de prospection sismique comprend un arrangement de sources d'énergie sismique. Cela est dû au fait qu'il est nécessaire de générer suffisamment d'énergie pour illuminer des structures en profondeur au sein de la Terre, et une source sismique unique ne peut généralement pas réaliser cela. 25 On connaît des sources d'énergie sismique qui émettent de l'énergie sismique à plus d'une fréquence. Parmi des exemples de telles sources sismiques, on trouve des sources de vibration, lesquelles émettent de l'énergie sismique dans une gamme de fréquences allant, par exemple, de 30 5 ou 10 Hz à 100 Hz. Lorsqu'une telle source de vibration est activée, de l'énergie sismique est émise sur une période de temps finie, et la fréquence de l'énergie émise change durant la période sur laquelle de l'énergie sismique est émise. Par exemple, la fréquence de l'énergie émise peut augmenter de façon monotone durant la période sur laquelle de l'énergie sismique est émise. Le processus consistant à faire fonctionner 5 une source de vibration d'énergie sismique afin de provoquer l'émission d'énergie sismique sur la gamme de fréquences du vibreur sera appelé ici "balayage" ("sweeping") du vibreur, et l'étape consistant à initier un balayage de vibreur sera appelée "activation" du vibreur. Chaque émission d'énergie sismique à partir d'un vibreur est connue sous le nom 10 de "tir". La période de temps sur laquelle de l'énergie sismique est émise par la source de vibration sera appelée "temps de balayage", et la "vitesse de balayage" est la vitesse à laquelle change la fréquence au cours du temps de balayage (en pratique, on utilise généralement une vitesse de balayage linéaire).
Une source de vibration sismique destinée à être utilisée sur terre consiste généralement en une embase en contact avec le sol. De l'énergie sismique est transmise dans le sol en appliquant une force vibratoire à la plaque, et cela est réalisé en appliquant une forme d'onde de commande 20 connue sous le nom de "balayage pilote" ("pilot sweep") au mécanisme de commande de vibreur. Le balayage pilote consiste généralement en un signal en fréquence balayée d'amplitude constante, quoiqu'il diminue progressivement à chaque extrémité pour permettre à l'amplitude de la vibration d'être montée et descendue par incréments au début et à la fin 25 du balayage respectivement. En pratique, la forme d'onde appliquée au sol par la plaque n'est pas exactement la même que la forme d'onde pilote; en particulier, en plus d'appliquer une force à la fréquence souhaitée à n'importe quel moment (connue sous le nom de "fréquence fondamentale"), le vibreur applique aussi une force à des multiples entiers 30 de la fréquence fondamentale (connus sous le nom d'harmoniques").
On connaît aussi des sources de vibration d'énergie sismique marines. Là encore, on leur fait faire un balayage afin d'émettre de l'énergie sismique sur une gamme de fréquences.
Lorsqu'une source de vibration sismique est activée pour émettre de l'énergie sismique, l'énergie sismique incidente sur un récepteur est enregistrée pendant une période pré-déterminée depuis le début du temps de balayage de la source. Le laps de temps de la fin du temps de balayage à la fin de la période d'enregistrement est généralement connu sous le 10 nom de "temps d'écoute", et des données sont acquises au niveau d'un récepteur depuis le début du temps de balayage jusqu'à la fin du temps d'écoute. Les données acquises au niveau d'un récepteur en conséquence de l'activation d'une source sont ensuite traitées, par exemple par intercorrélation des données acquises avec le balayage 15 pilote de la source pour produire un enregistrement qui soit de la longueur du temps d'écoute.
La figure 1 est une illustration schématique du processus d'une étude sismique classique qui utilise un arrangement de sources de vibration 20 terrestres. Au temps TO, une source sismique dans l'arrangement de sources est activée pour débuter son balayage. Dans cet exemple, le temps de balayage de vibreur a une durée S, et la fréquence de l'énergie sismique émise par le vibreur augmente de façon monotone depuis une fréquence f0 au temps TO jusqu'à une fréquence f1 (fi > fo) à la conclusion 25 du balayage (au temps TO + S). Le temps de balayage est suivi du temps d'écoute, de sorte que le temps total du processus d'activation et de balayage de la source ainsi que d'écoute au niveau d'un récepteur pour de l'énergie sismique est S + L, o L est la durée du temps d'écoute.
Dans une étude sismique classique, les sources sont activées de telle sorte qu'un récepteur ne recevra de l'énergie sismique que d'une seule source dans une période d'écoute donnée quelconque. Le retard minimum entre le début de deux balayages de vibreur dans une telle étude est donc la somme du temps de balayage S et du temps d'écoute L. Le temps d'écoute L est prévu suffisamment grand pour que toute l'énergie sismique requise au niveau d'un récepteur dans une période d'écoute particulière soit émise durant le temps de balayage précédant immédiatement cette période d'écoute.
Le processus de prospection sismique classique présente l'inconvénient 10 d'être potentiellement lent, du fait qu'il est nécessaire que le temps de retard minimum entre les débuts de deux balayages de vibreur soit la somme du temps de balayage et du temps d'écoute. La technique d'acquisition par "balayage glissant" est une tentative connue de réduction du temps requis pour réaliser une étude sismique. Dans la 15 technique par balayage glissant, le temps de retard minimum entre les débuts de deux balayages de vibreur subséquents est seulement le temps d'écoute, et non la somme du temps de balayage et du temps d'écoute. La longueur d'enregistrement après intercorrélation est la longueur du temps d'écoute.
La technique de balayage glissant est illustrée sur la figure 2. Comme dans la méthode de la figure 1, une source sismique dans l'arrangement de sources est activée afin de débuter son balayage au temps TO, le temps de balayage de vibreur a une durée S, et la période de balayage est 25 suivie d'un temps d'écoute L. Le temps Tl auquel une deuxième source est activée afin de débuter son balayage n'a toutefois pas besoin de satisfaire Tl > TO + S + L, mais doit seulement satisfaire Tl > TO + L. Du fait que le temps de retard minimum entre l'activation de deux sources dans la technique de balayage glissant est seulement le temps d'écoute, 30 et non la somme du temps de balayage et du temps d'écoute, la technique de balayage glissant permet de réduire le temps nécessaire pour achever une étude sismique. Toutefois, elle présente l'inconvénient que des harmoniques de la fréquence fondamentale générées par un vibreur sont présentes sur le sismogramme enregistré par un ou plusieurs vibreurs précédents.
Une technique de prospection connue supplémentaire est la technique de "tirs simultanés". Dans la méthode de tirs simultanés, deux sources sismiques ou plus disposées à des emplacements de tir respectifs sont activées afin de débuter leurs balayages en même temps. L'énergie sismique acquise au niveau d'un récepteur contiendra par conséquent des événements résultant de l'énergie sismique émise par toutes les sources. Afin de permettre aux évènements correspondant à chaque source d'être séparés les uns des autres, chaque vibreur doit effectuer un balayage, à son emplacement de tir, au moins autant de fois qu'il y a de 15 vibreurs dans le groupe, et les données enregistrées sont ensuite manipulées de façon algébrique pour séparer les évènements correspondant à chaque source. Typiquement, chaque vibreur effectuera un balayage pendant le même laps de temps, à la même vitesse de balayage et sur la même gamme de fréquences, mais la relation de phases 20 entre des vibreurs change d'un enregistrement à un autre. Dans le cas d'un groupe de deux vibreurs, par exemple, un schème adéquat serait que les deux vibreurs effectuent un balayage en phase durant le premier enregistrement et effectuent un balayage déphasés de 1800 durant le deuxième enregistrement. La moyenne des deux signaux acquis par un 25 récepteur donne le signal naissant au niveau de ce récepteur suite à l'activation d'un vibreur, et la moitié de la différence des deux signaux acquis par un récepteur donne le signal naissant au niveau de ce récepteur suite à l'activation de l'autre vibreur.
Un premier aspect de la présente invention fournit une méthode de prospection sismique comprenant les étapes consistant: à activer le ou chaque vibreur dans un premier groupe de vibreurs au temps TO; et à activer le ou chaque vibreur dans un deuxième groupe de vibreurs au temps Tl, o TO < Tl < TO + Si + L o Si est le temps de balayage du premier groupe de vibreurs et L est le temps durant lequel des données sont acquises au niveau d'un récepteur; dans laquelle au moins un groupe parmi le premier groupe de vibreurs et le deuxième groupe de vibreurs comprend au moins deux vibreurs.
La présente invention fournit une méthode de prospection sismique qui 10 combine la technique d'acquisition simultanée connue et la technique d'acquisition par balayage glissant connue. Elle fait usage de la technique d'acquisition simultanée en ce qu'au moins l'un des groupes de vibreurs contient deux vibreurs ou plus, mais les différents groupes sont balayés à l'aide d'une technique de balayage glissant. La présente invention fournit 15 une réduction du temps requis pour réaliser une étude sismique par comparaison avec le temps requis par une technique de balayage glissant classique.
L'invention peut être appliquée à un groupe de vibreurs unique, auquel cas le deuxième groupe de vibreurs est le premier groupe de vibreurs, ou bien les premier et deuxième groupes de vibreurs peuvent être distincts.
Dans un mode de réalisation préféré, la méthode comprend l'étape consistant à activer le ou chaque groupe de vibreurs au moins autant de 25 fois qu'il y a de vibreurs dans le groupe respectif et de telle sorte que la contribution de chaque vibreur dans un groupe puisse être déterminée par une opération algébrique sur des enregistrements de récepteur de balayages réalisés par le groupe respectif.
Dans un mode de réalisation préféré, Tl - TO > (n - 1) Sl fo/n (f i - f o) , o n est un entier naturel, f O est la limite de fréquence inférieure du balayage de vibreur et f t est la limite de fréquence supérieure du balayage de vibreur. Cela permet au bruit dans les données acquises résultant de la miè,e harmonique de la fréquence fondamentale d'être estimé pour tout m:S n.
Dans un mode de réalisation préféré, la méthode comprend les étapes consistant: à activer le ou chaque vibreur dans le premier groupe de vibreurs au temps T2, o Tl < T2 < Tl + S2 + L o S2 est le temps de balayage du deuxième groupe de vibreurs; et à activer le ou chaque vibreur dans le deuxième groupe de vibreurs au temps T3 o T2 < T3 < T2 + S + L; le premier groupe de vibreurs est différent du deuxième groupe de vibreurs; et T3 - T2 É Tl - TO.
Ce mode de réalisation permet d'appliquer une technique de réduction 15 supplémentaire de bruit. Le temps de retard variable entre un tir du premier groupe et le tir d'un deuxième groupe signifie que du bruit harmonique se produira à des temps différents dans les enregistrements des deux tirs. Le bruit peut donc être éliminé en faisant de façon appropriée la somme des deux enregistrements de tir, en supposant que 20 chaque enregistrement de tir contienne le même signal.
Un deuxième aspect de la présente invention fournit une méthode de prospection sismique comprenant les étapes consistant: à activer un premier vibreur au temps TO; à activer un deuxième vibreur différent du 25 premier vibreur au temps Tl, o TO < Tl < TO + Si + L o Si est le temps de balayage du premier vibreur et L est le temps durant lequel des données sont acquises au niveau d'un récepteur; à activer le premier vibreur au temps T2, o Tl < T2 < Tl + S2 + L o S2 est le temps de balayage du deuxième vibreur; et à activer le deuxième vibreur au temps 30 T3oT2<T3<T2+S1 +LetoT3-T2 #T1 -TO.
Dans un mode de réalisation préféré T3 - T2 > Tl - TO.
Un troisième aspect de la présente invention fournit un agencement de prospection sismique comprenant: une pluralité de sources de vibration; 5 et un moyen de commande adapté pour activer le ou chaque vibreur dans un premier groupe de vibreurs au temps TO et pour activer le ou chaque vibreur dans un deuxième groupe de vibreurs au temps Tl, o TO < Tl < TO + Si + L o Si est le temps de balayage du premier groupe de vibreurs et L est le temps durant lequel des données sont acquises au niveau d'un 10 récepteur.
Un quatrième aspect de la présente invention fournit un agencement de prospection sismique comprenant: une pluralité de sources de vibration; et un moyen de commande adapté pour: (a) activer un premier vibreur au 15 temps TO; (b) activer un deuxième vibreur différent du premier vibreur au temps Tl, o TO < Ti < TO + Si + L o Si est le temps de balayage du premier vibreur et L est le temps durant lequel des données sont acquises au niveau d'un récepteur; (c) activer le premier vibreur au temps T2, o Tl < T2 < Tl + S2 + L o S2 est le temps de balayage du deuxième vibreur; et (d) activer le deuxième vibreur au temps T3 o T2 < T3 < T2 + Si + L et o T3 - T2 * Tl - TO.
Dans un mode de réalisation préféré, le moyen de commande comprend un processeur de données programmable. 25 Un cinquième aspect de l'invention fournit un support contenant un programme destiné au processeur de données d'un agencement de prospection sismique tel que défini plus haut.
Des modes de réalisation préférés de la présente invention vont maintenant être décrits, à titre d'exemples illustratifs, en référence aux figures qui les accompagnent, sur lesquelles: la figure 1 est un chronogramme schématique d'une étude sismique classique; la figure 2 est un chronogramme schématique d'une étude sismique par balayage glissant classique; la figure 3 est un chronogramme schématique d'une étude sismique selon 10 un premier mode de réalisation de la présente invention; la figure 4 est un organigramme schématique fonctionnel d'un mode de réalisation de la présente invention; la figure 5 est une illustration schématique de données sismiques brutes; la figure 6 illustre le résultat du traitement des données sismiques brutes de la figure 5 selon une première méthode de la présente invention; la figure 7 illustre les résultats du traitement des données sismiques de la figure 5 selon une deuxième méthode de la présente invention; la figure 8 montre le résultat du traitement des données sismiques de la figure 5 selon les première et deuxième méthodes de l'invention; 25 la figure 9(a) est un diagramme schématique d'un agencement de prospection sismique selon la présente invention; et la figure 9(b) est un diagramme schématique d'un moyen de commande 30 de l'agencement de prospection sismique de la figure 9(a).
Le fonctionnement d'un agencement de prospection sismique selon un premier mode de réalisation de la présente invention est illustré schématiquement sur la figure 3, laquelle est un chronogramme de la méthode. Cette méthode permet une réduction du temps pris pour 5 achever une étude sismique par comparaison avec la méthode de balayage glissant classique décrite plus haut.
La méthode suppose que l'étude sismique a un arrangement de sources de vibration sismiques. Afin de mettre en oeuvre la méthode, les vibreurs 10 sont groupés en deux groupes ou plus, un groupe de vibreurs pouvant fonctionner indépendamment du ou de chaque autre groupe de vibreurs.
Le groupement peut être un groupement physique, les vibreurs étant par exemple agencés en un arrangement à deux dimensions, chaque rangée ou colonne constituant un groupe. En alternative, le groupement peut être 15 un groupement fictif, dans lequel des vibreurs dans un arrangement sont divisés de façon fictive en deux groupes ou plus. Il n'est pas nécessaire que chaque groupe ait le même nombre de vibreurs, mais un groupe au moins doit contenir deux vibreurs ou plus.
Au temps TO un premier groupe de vibreurs est activé. C'est à dire que chaque vibreur dans le premier groupe débute son balayage au temps TO.
Tous les vibreurs dans le premier groupe auront le même temps de balayage. (Pour chaque groupe, tous les vibreurs dans le groupe auront le temps de balayage. En principe, différents groupes pourraient avoir différents temps de balayage, bien qu'il soit courant que les groupes aient tous le même temps de balayage les uns les autres.) Si le temps de balayage des vibreurs dans le premier groupe est Si, alors la période de balayage du premier groupe de vibreurs dure du temps TO au temps TO + Si. La période de balayage pour le premier groupe est ensuite suivie par 30 un temps d'écoute, lequel a une durée L et se conclut donc au temps TO + Si + L. il Un deuxième groupe de vibreurs est activé afin de débuter leurs balayages au temps Tl. Le temps Tl est avant la fin de la période d'écoute du premier groupe de vibreurs, c'est à dire que Tl < TO + Si + L. Le groupe de vibreurs activé au temps Tl est appelé deuxième groupe de 5 vibreurs par souci de commodité, mais le deuxième groupe de vibreurs pourrait soit être différent du premier groupe de vibreurs, soit être identique à celui-ci. La période de balayage du deuxième groupe de vibreurs a une durée S2 (qui peut être égale à Si) et donc s'étend de Tl à Tl + S2, et la période d'écoute s'étend alors jusqu'à Tl + S2 + L (les deux 10 groupes auront le même temps d'écoute). Il faut que le temps de retard entre le début du balayage du premier groupe de vibreurs et le début du balayage du deuxième groupe de vibreurs dépasse le temps d'écoute (et si l'on fait effectuer un balayage au même groupe deux fois de suite, alors il faut que le retard soit supérieur au temps de balayage de ce groupe). 15 La figure 3 représente seulement le premier balayage du premier groupe et le premier balayage du deuxième groupe. Toutefois, la méthode de l'invention requiert que l'on fasse effectuer un balayage à chaque groupe de vibreurs au moins autant de fois qu'il y a de vibreurs dans ce groupe, 20 suivant une séquence à partir de laquelle la contribution provenant de chaque vibreur aux données sismiques acquises au niveau d'un récepteur puisse être déterminée, par exemple de façon algébrique.
En général, il y aura M groupes de vibreurs, avec N vibreurs dans chaque 25 groupe. Si chaque groupe procède à K tirs, o K 2 N, il y aura une séquence de MK tirs. Par exemple, s'il y a deux groupes de trois vibreurs (M = 2, N = 3), il faut que chaque groupe procède à 3 tirs au moins (K 2 3).
Si chaque groupe procède à exactement trois tirs, une séquence possible serait: groupe 1; groupe 2; groupe 1; groupe 2; groupe 1; groupe 2. 30 Une certaine différentiation entre les balayages, par exemple un changement de phase, est requise afin de permettre aux enregistrements d'être manipulés de façon algébrique pour isoler la contribution de chaque vibreur. (En principe, si un groupe devait procéder à plus de tirs qu'il n'y a de vibreurs dans le groupe (c'est à dire si K > N), certains des tirs pourraient être identiques, pourvu qu'il y ait au moins N tirs indépendants.) Comme autre exemple, considérons un groupe unique de 4 (M = 1, N = 4).
Dans ce cas, la valeur la plus basse de K est 4, si bien qu'on doit faire effectuer un balayage au groupe au moins 4 fois et la séquence la plus 10 courte possible est: groupe 1; groupe 1; groupe 1; groupe 1.
Comme autre exemple, considérons six vibreurs agencés en trois groupes, de deux vibreurs chacun (M = 3, N = 2). Dans ce cas, la valeur la plus basse de K est 2. Si chaque groupe procède à exactement deux tirs, 15 une séquence possible serait: groupe 1; groupe 2; groupe 3; groupe 1 groupe 2; groupe 3.
Dans cette méthode de l'invention, au moins un des tirs a lieu dans la limite du temps S + L après le début du tir précédent. De préférence, autant de tirs que possible débutent moins de S + L après le tir précédent.
On notera que le tout dernier tir ne contiendra pas de bruit de balayage glissant - par exemple, dans l'exemple simple consistant à tirer avec seulement deux vibreurs, le deuxième vibreur tirerait dans la limite du temps d'écoute du premier vibreur (c'est à dire dans la limite du temps S + 25 L après le début du premier tir). Les deux vibreurs seraient ensuite déplacés jusqu'à leurs prochains points de tir, de sorte que le tir du deuxième vibreur n'aurait pas de bruit de balayage glissant sur lui.
Dans les exemples ci-dessus, chaque groupe de vibreurs contient le 30 même nombre de vibreurs. L'invention n'est pas limitée à cela, et les groupes n'ont pas besoin de tous contenir le même nombre de vibreurs (bien qu'il faille un groupe contenant plus d'un vibreur). Si les groupes ne contiennent pas le même nombre de vibreurs, on doit faire effectuer un balayage au moins une fois de plus qu'il n'y a de vibreurs dans ce groupe.
La figure 3 indique que les vibreurs dans un troisième groupe de vibreurs (qui pourrait être le premier ou le deuxième groupe) sont activés afin de débuter leurs balayages au temps T2 qui satisfait T2 < Tl + S2 + L, o S2 est le temps de balayage du deuxième groupe de vibreurs. Il serait toutefois possible que le troisième groupe soit activé à un temps T2 > Tl 10 + S2 + L puisqu'un tir seulement doit débuter dans la limite de S + L après le début du tir précédent.
Une fois que chaque groupe de vibreurs a été activé un nombre de fois suffisant, les vibreurs peuvent être déplacés à des positions différentes 15 pour permettre à une nouvelle étude d'être menée.
La réponse au niveau de chaque récepteur due à l'énergie sismique générée par chaque vibreur individuel peut être calculée à partir de l'énergie sismique acquise au niveau du (des) récepteur(s) en utilisant 20 n'importe quelle technique appropriée. Une technique appropriée est révélée dans la demande de brevet du Royaume-Uni N0 2 359 363. En principe, la réponse pourrait être calculée immédiatement après que les données aient été acquises, mais les données sont plus souvent stockées, par exemple sur bande ou disque magnétique, en vue d'un 25 traitement ultérieur loin du lieu d'étude.
Un problème lié à l'acquisition par balayage glissant est que les données obtenues au cours de la période d'écoute suivant un balayage contiendront du bruit harmonique qui résulte du balayage subséquent. Il 30 est souhaitable d'éliminer ce bruit harmonique durant le traitement des données. Un mode de réalisation supplémentaire de la présente invention fournit une méthode pour estimer le bruit harmonique dans des données acquises au cours d'une période d'écoute d'après les données acquises.
Le bruit harmonique dans des données acquises au cours d'une période d'écoute est la réponse de la Terre à la sortie harmonique provenant du balayage subséquent. La réponse de la Terre pour les harmoniques du tir subséquent est la même que la réponse de la Terre au balayage fondamental du tir associé à la période d'écoute. La théorie de déconvolution de vibreur est basée sur la connaissance du contenu du 10 balayage fondamental car cela fournit une mesure de la réponse de la Terre. S'il est possible d'estimer aussi la sortie harmonique du vibreur, il est alors possible d'estimer la contribution harmonique au tir précédent en convoluant la réponse de la Terre, telle que déterminée d'après le contenu du balayage fondamental, avec la sortie harmonique du vibreur. 15 Une fois que la contribution harmonique au tir précédent a été trouvée de cette manière, il est possible de la soustraire des données enregistrées.
En principe, il serait possible d'utiliser des mesures sur l'énergie sismique émise au niveau du vibreur pour estimer la sortie harmonique du vibreur. 20 Toutefois, les effets de la non-linéarité du comportement de la Terre au voisinage du vibreur peuvent signifier que le bruit harmonique au niveau du récepteur diffère du bruit harmonique qui serait estimé d'après des mesures sur le vibreur. L'estimation de la sortie harmonique à partir des données elles-mêmes est donc une méthode plus fiable. 25 Un mode de réalisation d'une méthode d'élimination du bruit harmonique de données acquises au cours de la technique d'acquisition par balayage glissant simultanée de la figure 3 va à présent être décrit en référence à la figure 4.
Initialement, à l'étape 1, des sources de vibration sont groupées en deux groupes ou plus. Cette étape de groupement peut être une étape de groupement physique, par exemple durant le déploiement des vibreurs.
En alternative, elle peut consister à définir des groupes fictifs dans un arrangement de vibreurs. Au moins l'un des groupes contient deux vibreurs ou plus.
A l'étape 2, des données de balayage glissant simultanées sont acquises à l'aide d'une technique d'acquisition du type décrit généralement en rapport à la figure 3.
A l'étape 3, la réponse impulsionnelle est calculée pour chaque groupe, à partir des positions des vibreurs dans ce groupe. Cette étape peut être réalisée d'une manière classique à l'aide de l'intercorrélation et d'une matrice d'inversion qui est constante avec la fréquence, ou elle peut en 15 alternative être réalisée de la manière décrite dans le document GB-A2 359 363.
Puis, à l'étape 4, les données acquises sont intercorrélées avec un balayage de fréquence harmonique. Le balayage qui est utilisé dans cette 20 étape d'intercorrélation est essentiellement le même que le balayage fondamental du vibreur, sauf qu'il a une vitesse de balayage qui est un multiple entier de la vitesse de balayage du balayage fondamental. Ainsi, le balayage utilisé dans l'étape d'intercorrélation aurait deux fois la vitesse de balayage du balayage de vibreur fondamental si l'on souhaite 25 retirer la deuxième harmonique, il aurait trois fois la vitesse de balayage du balayage fondamental pour retirer la troisième harmonique, et ainsi de suite. Il n'est pas nécessaire que le balayage utilisé dans l'étape d'intercorrélation s'étende au-dessus de la limite de fréquence supérieure du balayage fondamental - et dans tous les cas, la fréquence supérieure 30 ne devrait pas dépasser la fréquence de Nyquist pour l'échantillonnage si bien qu'un dispositif de diminution progressive est appliqué au balayage harmonique après une fréquence au voisinage de la fréquence supérieure du balayage fondamental. Pour chaque harmonique, lemême balayage harmonique peut être utilisé dans l'étape d'intercorrélation pour tous les récepteurs et pour tous les tirs.
L'effet de l'étape d'intercorrélation est de concentrer dans le temps la réponse de la Terre à une harmonique. La réponse de la Terre à d'autres harmoniques et à la fréquence fondamentale est étalée dans le temps. En particulier, bien que la réponse de la Terre à la fréquence de balayage fondamental reste le principal contributeur aux données acquises, l'effet de l'étape d'intercorrélation est qu'elle arrivera désormais plus tard dans les enregistrements.
A l'étape 5, les données intercorrélées sont transformées par la transformée de Fourier sur une fenêtre de temps qui comporte la première arrivée principale due à l'harmonique souhaitée, mais qui finit avant que la réponse de la Terre à la fréquence de vibreur fondamentale soit observée ou (dans le cas de l'estimation d'harmoniques plus élevées) lorsque la réponse de la Terre à une harmonique d'amplitude plus élevée 20 est observée. La longueur de fenêtre de temps utilisée dans cette étape peut varier d'un récepteur à un autre. Pour chaque fréquence et récepteur, le résultat de la transformée de Fourier est un vecteur Gn de longueur N, o N est le nombre de tirs. Ce vecteur est la partie précoce de la réponse de la Terre à l'harmonique de vibreur choisie. 25 A l'étape 6, la réponse impulsionnelle pour chaque récepteur et position de vibreur est traitée avec une transformée de Fourier de la même longueur que dans l'étape 5. La transformée de Fourier est réalisée sur la même fenêtre de temps, ou une fenêtre de temps légèrement plus courte, 30 que la transformée de Fourier dans l'étape 5. Les résultats de cette étape de transformée de Fourier sont un vecteur Rn pour chaque récepteur et fréquence, Rn ayant une longueur M o M est le nombre de positions de vibreur.
A l'étape 7, la matrice GnR*n est calculée pour chaque fréquence et récepteur. R*n est le nombre complexe conjugué de la transposée de Rn.
Cela donne une matrice N x M. On fait ensuite la moyenne de cette matrice sur des récepteurs pour donner la matrice d'intercorrélation [GnR*,J. Une normalisation basée sur les récepteurs peut être appliquée lors de ce processus de moyenne, afin de tenir compte de l'amplitude de signal 10 variable au niveau de récepteurs différents.
A l'étape 8, la matrice M x M, RnR*n, est calculée pour chaque fréquence et récepteur. On fait ensuite la moyenne de cela sur des récepteurs pour donner la matrice d'autocorrélation [RR*n]. Si une normalisation basée 15 sur des récepteurs était appliquée dans l'étape 7, la même normalisation devrait être appliquée dans l'étape 8.
Dans l'étape 9, on détermine l'estimation de la nième harmonique H, émise par le vibreur. En principe, H1 est donnée par Hn = [GnR*n][RnR*nJ( l).
Toutefois, en pratique, la matrice d'autocorrélation sera mal conditionnée - c'est-à-dire que la plus petite valeur propre sera beaucoup plus petite que la plus grande valeur propre, et prendre un inverse exact entraîne le risque d'une domination par du bruit. Une décomposition en valeurs singulières standard de [RR*nl la décompose sous la forme du produit de 25 trois matrices, U, V et A, o [RnR*n] = U A VI, U et V sont telles que UU* = V* = I o I est la matrice identité et A est diagonale et réelle. U et V sont les matrices de vecteurs propres droite et gauche respectivement, A se compose des valeurs propres. La matrice inverse vraie est donnée par VA-1U*. L'estimation de Hn est donnée par Hr = [GnR*n](VLU*) o L est une 30 matrice diagonale, identique à A(-1), sauf que les plus petits éléments ont été remplacés par des zéros. Une façon de décider combien de valeurs propres doivent être retenues et combien doivent être mises à zéro vient de la comparaison de la dimension moyenne quadratique (rms) des éléments de [RnR*nl. La dimension de la rms des éléments de VLU* multipliée par la dimension de la rms de [RnR*j1 devrait être d'ordre 1.
Souvent, une seule valeur propre est nécessaire. Hn est une matrice N fois M pour chaque fréquence et série de tirs.
A l'étape 10, la matrice Hn est transformée par transformée de Fourier pour retourner au domaine temporel pour obtenir hn. La transformée de 10 Fourier est réalisée pour une fenêtre de temps autour de T = 0 à l'aide d'une diminution progressive continue et ayant une demi largeur d'environ 0,25 seconde ou moins.
A l'étape 11, hn est convolué avec le balayage harmonique utilisé dans 15 l'étape 4. Le résultat de cette convolution est l'estimation de la sortie de vibreur pour la nième harmonique. Cela est une estimation du bruit qui apparaît sur les données acquises durant le tir précédent. Les données acquises pour le tir précédent peuvent être corrigées en ce qui concerne le bruit harmonique, par exemple en soustrayant le bruit harmonique estimé à partir des données brutes acquises au niveau du récepteur. En alternative, si le stade de déconvolution est exécuté sur des données corrélées, le bruit harmonique peut être corrélé avec le balayage pilote approprié, puis soustrait des enregistrements corrélés appropriés.
Les étapes ci-dessus peuvent être répétées pour chaque harmonique que l'on souhaite retirer. De façon typique, les deuxième et troisième harmoniques ont la plus grande amplitude, de sorte que le retrait de ces harmoniques seulement peut être suffisant.
A l'étape 12, les enregistrements de tirs, dont le bruit harmonique a maintenant été retiré, sont séparés pour obtenir leurs composantes de point de tir individuelles. Cette étape peut être réalisée sur des enregistrements de tirs soit non corrélés, soit corrélés. En alternative, le stade de séparation peut être réalisé sur la somme des estimations de bruit seule, et ces estimations de bruit séparées sont ensuite soustraites de composantes de points de tir individuelles calculées auparavant.
La méthode ci-dessus a été décrite en référence à la méthode d'acquisition par sous-balayage simultanée de la figure 3. Elle peut en alternative être appliquée à une méthode d'acquisition par balayage glissant classique du type montré sur la figure 2. La méthode peut être simplifiée lorsqu'elle est appliquée à une telle technique par balayage glissant classique, puisque l'intercorrélation et l'autocorrélation ne sont pas des matrices, mais sont un nombre à chaque fréquence. Ainsi, dans l'étape 9, l'inverse vrai de l'autocorrélation peut être utilisé. 15 L'utilisation de cette méthode requiert un certain temps minimum entre les balayages. Pour retirer la nième harmonique, le temps minimum entre des balayages consécutifs est donné par: (n min n(f, fo)(1 Cela laisse supposer que le balayage fondamental varie de façon linéaire entre la fréquence fo et la fréquence f1 et a un temps de balayage de S. 25 Ainsi, pour retirer la deuxième harmonique, le glissement (slip) entre des balayages consécutifs (par exemple Tl - TO sur la figure 2 ou la figure 3) doit faire légèrement plus de la moitié de la durée de balayage totale; pour retirer la troisième harmonique, il doit faire légèrement plus de deux 30 tiers du temps de balayage, etc. Les figures 5 et 6 illustrent des résultats de la méthode de la figure 4. La figure 5 montre une partie d'un enregistrement de données acquises au niveau d'un récepteur durant l'activation simultanée de trois vibreurs sismiques, chaque vibreur étant activé quatre fois. Bien que les données 5 aient été acquises avec une technique d'acquisition simultanée classique plutôt qu'une technique d'acquisition par balayage glissant simultanée, les enregistrements ont été sommés pour simuler des données acquises par une technique de balayage glissant simultanée de la présente invention. On verra qu'il y a un bruit harmonique de haute fréquence sur 10 une grande partie des données. Les douze traces montrées sur la figure 5 ont été acquises au niveau de douze emplacements de récepteurs distincts.
La figure 6 illustre le résultat de l'application d'une méthode d'estimation 15 et de retrait d'harmoniques du type décrit en référence à la figure 4 aux deuxième et troisième harmoniques dans les données de la figure 5. On peut voir que le bruit harmonique de haute fréquence a été réduit de façon significative.
Une technique alternative destinée à retirer de l'énergie harmonique va être décrite à présent. En contraste avec la méthode de réduction de bruit harmonique décrite ci-dessus, laquelle repose sur la capacité à estimer le bruit harmonique, la technique décrite ci-dessous requiert très peu de connaissances de l'énergie harmonique. La méthode requiert seulement 25 de savoir à quel moment différents ordres d'énergie harmonique arrivent à un récepteur. Cette méthode peut être appliquée à la fois à l'acquisition par balayage glissant simultanée de la figure 3 et à l'acquisition par balayage glissant classique de données sismiques. Une restriction relative à cette technique est qu'il doit être procédé à au moins un balayage de plus que le nombre de points de tir qui sont en cours de séparation. Ainsi, pour une acquisition de données par balayage glissant standard, chaque vibreur doit être activé au moins deux fois dans chaque emplacement de tir. Dans une technique d'acquisition par balayage glissant simultanée dans laquelle chaque groupe comporte deux vibreurs, il doit être procédé à au moins trois balayages dans chaque emplacement.
Considérons initialement le cas d'une technique d'acquisition par balayage glissant standard dans laquelle chaque groupe de vibreurs procède à deux tirs dans chaque emplacement. Après intercorrélation ou déconvolution, les composantes de signaux de l'enregistrement de données pour un tir seront les mêmes que la composante de signaux de l'enregistrement de données pour les autres tirs. Si la composante de bruit harmonique dans les deux enregistrements peut être agencée pour différer, il est alors possible de réduire ou d'éliminer le bruit harmonique en combinant les deux enregistrements de façon appropriée. Une façon 15 de faire cela qui a déjà été proposée consiste à faire varier la phase des tirs. Si la phase de la nième harmonique est n fois plus grande que la phase du balayage fondamental (ce qui est d'habitude à peu près vrai), alors si l'on choisit une différence de phase appropriée et que l'on somme les enregistrements après intercorrélation, une ou plusieurs harmoniques 20 s'annuleront. Si les deux tirs sont déphasés de 900, alors la troisième harmonique s'annulera. Si les deux tirs sont déphasés de 1800, alors les deuxième et quatrième harmoniques s'annuleront. S'il y a trois tirs, chacun étant déphasé de 1200 par rapport aux deux autres tirs, alors les deuxième et troisième harmoniques s'annuleront. 25 La méthode décrite ci-dessous utilise une combinaison d'une technique d'acquisition et d'une technique de traitement, et ne dépend d'aucune relation de phase algébrique.
La méthode décrite ci-dessous est basée sur le principe selon lequel, si le bruit harmonique lors d'un tir apparaît à un temps différent du bruit harmonique lors d'un autre tir, alors un empilement (stacking) approprié peut éliminer ou substantiellement réduire le bruit. La méthode d'empilement qui est utilisée est un empilement à diversité en tempsfréquence. Cette méthode d'empilement fonctionne bien lorsque la composante de signal des deux enregistrements est la même, mais le bruit apparaît dans différents emplacements dans le domaine tempsfréquence.
Pour faire apparaître du bruit harmonique à différents temps dans différents tirs, une méthode consiste à varier le retard entre les tirs. La différence de temps choisie entre un tir dans un emplacement et le tir subséquent (à un autre emplacement) devrait différer du retard de temps entre un autre tir auquel il a été procédé dans le premier lieu et le tir subséquent à celui-ci de telle sorte que, après décomposition en temps15 fréquence, les plus grands pics dans le bruit ne se recouvrent pas substantiellement l'un l'autre. Le retard de temps exact qui est requis dépendra de la vitesse de balayage, mais on a découvert qu'une différence d'entre 1 et 2 secondes est généralement suffisante si le temps de balayage est de moins de 10 secondes. 20 Dans une variation de cette technique, l'effet du retard de temps entre des tirs est imité à l'aide de groupes de vibreurs multiples o le bruit est agencé pour provenir de tirs de groupes de vibreurs différents, largement séparés. La séparation spatiale entre les groupes de vibreurs induira un 25 retard de temps tel que, pour la plupart des récepteurs, le bruit harmonique arrivera à des temps différents et un empilement à diversité en temps-fréquence sera efficace pour retirer le bruit harmonique. Pour certains récepteurs, les pics dans le bruit provenant de différents tirs se recouvriront l'un l'autre, mais ce bruit peut être retiré durant l'empilement 30 à l'aide de la technique révélée dans la demande de brevet du RoyaumeUni N' 2 359 363.
Dans le cas de N vibreurs, il est nécessaire de procéder à au moins (N + 1) tirs de balayages encodés, à partir desquels les n enregistrements de tirs individuels peuvent être extraits. Considérons par exemple deux groupes de deux vibreurs chacun avec un temps de balayage S de 8 secondes et 5 un temps d'écoute L de 5 secondes. Ces vibreurs peuvent exécuter trois tirs pour chaque groupe de la façon la plus rapide comme suit:
Tableau 1
Groupe 1 tir 1 débute à 0 seconde Groupe 1 tir 1 finit à 8 secondes Groupe 1 tir 2 débute à 10 secondes Groupe 1 tir 2 finit à 18 secondes Groupe 1 tir 3 débute à 20 secondes Groupe 1 tir 3 finit à 28 secondes Groupe 2 tir 1 débute à secondes Groupe 2 tir 1 finit à 13 secondes Groupe 2 tir 2 débute à 15 secondes Groupe 2 tir 2 finit à 23 secondes Groupe 2 tir 3 débute à 25 secondes Groupe 2 tir 3 finit à 33 secondes Temps de glissement secondes Temps de glissement secondes Temps de glissement secondes Temps de glissement secondes Temps de glissement secondes Temps de glissement secondes L'addition de retards entre les tirs auxquels le premier groupe a procédé et le tir qui suit chaque tir du premier groupe permet au bruit harmonique d'être placé à différents temps dans les enregistrements de tirs. Par exemple, le schème de tirs ci- dessus peut être modifié en augmentant le 15 retard entre un tir du premier groupe et le tir correspondant du deuxième groupe, comme suit:
Tableau 2
Groupe 1 tir 1 débute à 0 seconde Groupe 2 tir 1 débute à Temps de glissement secondes 5 secondes Groupe 1 tir 1 finit à 8 secondes Groupe 1 tir 2 débute à secondes Groupe 2 tir 1 finit à Temps de glissement 13 secondes 5 secondes Groupe 2 tir 2 débute à Temps de glissement 16 secondes 6 secondes Groupe 1 tir 2 finit à 18 secondes Groupe 1 tir 3 débute à 22 secondes Groupe 2 tir 2 finit à Temps de glissement 24 secondes 6 secondes Groupe 2 tir 3 débute à Temps de glissement 29 secondes 7 secondes Groupe 1 tir 3 finit à secondes Groupe 2 tir 3 finit à Temps de glissement 37 secondes 7 secondes Dans ce schème de tirs modifié du Tableau 2, le retard variable entre un tir du premier groupe et le tir correspondant du deuxième groupe signifie que le bruit harmonique apparaîtra à des temps différents dans les enregistrements de tirs. Il est donc possible d'éliminer le bruit harmonique à l'aide de la "technique d'inversion simultanée de diversité" révélée dans le document GB 2 359 363.
La figure 7 illustre les résultats de l'application de la technique d'inversion simultanée de diversité aux données montrées sur la figure 5. 15 On observera à nouveau qu'une réduction considérable du bruit harmonique a été fournie.
Dans le schème de tirs du Tableau 2, le retard de temps entre un tir du premier groupe de vibreurs et le tir correspondant du deuxième groupe de vibreurs augmente avec le nombre de tirs. L'invention n'est pas limitée à cela, et le retard de temps entre un tir du premier groupe de vibreurs et le tir correspondant du deuxième groupe de vibreurs pourrait en alternative diminuer avec le nombre de tirs.
La technique d'inversion simultanée de diversité peut être combinée avec l'estimation et le retrait d'harmoniques pour les tirs o le retard de temps satisfait l'équation (1) ci-dessus. Par exemple, si le balayage de vibreur va de 10 Hz à 100 Hz, alors la deuxième harmonique peut être estimée et retirée si le temps de glissement est supérieur à 4,5 secondes, la 15 troisième harmonique peut être estimée et retirée si le temps de glissement fait plus de 6 secondes, et la quatrième harmonique peut être retirée si le temps de glissement fait plus de 6,7 secondes. Ainsi, pour le schème d'acquisition du Tableau 2, la deuxième harmonique peut être estimée et retirée du bruit harmonique dans tous les tirs, la troisième 20 harmonique peut être estimée et retirée du bruit dans le deuxième et le troisième tir de chaque groupe, et la quatrième harmonique peut être estimée et retirée du bruit dans le troisième tir de chaque groupe.
La figure 8 illustre les résultats de l'application de la technique d'estimation et de retrait d'harmoniques ainsi que de la technique d'inversion simultanée de diversité. La figure 8 montre les résultats liés à la réalisation de la technique d'inversion simultanée de diversité sur les données montrées sur la figure 6 - qui ont déjà subi une technique d'estimation et de retrait d'harmoniques pour les deuxième et troisième 30 harmoniques. On verra que l'application subséquente de la technique d'inversion simultanée de diversité a eu pour résultat une réduction supplémentaire du bruit.
La manière dont les techniques ci-dessus sont utilisées dépend de facteurs tels que le nombre de groupes, le temps requis pour déplacer un vibreur d'une position de tir souhaitée à une autre et le rapport entre le temps de balayage de vibreur total et le temps d'écoute.
Si le temps d'écoute est égal au temps de balayage complet divisé par le nombre de vibreurs dans un groupe, alors une implémentation efficace 10 consiste à avoir deux groupes. A un temps quelconque, un groupe peut être en cours de repositionnement tandis que l'autre groupe est en cours d'acquisition de données. Le groupe qui est en cours d'acquisition de données tire le même nombre de tirs qu'il y a de vibreurs dans le groupe.
Chaque vibreur effectue des balayages presque en continu, avec seulement une courte pause à la fin de chaque balayage pour réinitialiser l'équipement. Lorsque le premier groupe a terminé de tirer, le deuxième groupe est en position et commence à tirer tandis que le premier groupe se déplace jusqu'au prochain emplacement. Le temps de glissement fait % du temps de balayage, si bien que l'estimation et le retrait d'harmoniques peuvent être appliqués aux données acquises. Toutefois, comme le nombre de tirs est égal au nombre de vibreurs dans chaque groupe, une inversion simultanée de diversité n'est pas possible.
Si, dans un autre exemple, le temps de balayage pour chaque vibreur fait 25 environ deux fois le temps d'écoute, alors une méthode utilisant trois groupes peut être employée à des fins utiles. Dans le cas de neuf vibreurs, il sera possible d'avoir trois groupes de trois vibreurs, les vibreurs dans chaque groupe étant activés quatre fois avec des retards variables entre les tirs. Le modèle suivant lequel les tirs sont effectués 30 serait agencé de façon à ce que les deux premiers tirs de chaque groupe alternent avec les deux derniers tirs du groupe précédent et à ce que les troisième et quatrième tirs de chaque groupe alternent avec les deux premiers tirs du groupe subséquent. Pour les tirs qui satisfont l'équation (1), le bruit harmonique peut être estimé et retiré. Une inversion simultanée de diversité peut ensuite être appliquée aux données acquises pour chaque groupe.
La figure 9(a) est une illustration schématique d'un agencement de prospection sismique selon un mode de réalisation de l'invention.
L'agencement de prospection sismique comprend une pluralité de vibreurs 1, 2, 3, 4 et un moyen de commande 5. Deux groupes de vibreurs A, B sont définis sur la figure 9(a), chaque groupe contenant deux vibreurs dans chaque groupe, mais l'agencement de prospection sismique de l'invention n'est pas limité à ce nombre de groupes particulier ou à ce nombre de vibreurs par groupe. Le moyen de commande 5 est capable d'activer chaque groupe indépendamment de l'autre groupe. Par exemple, le moyen de commande peut être connecté électriquement à chaque groupe de vibreurs de façon à ce qu'il puisse envoyer un signal électrique à un groupe sélectionné pour activer chaque vibreur dans le groupe sélectionné. Le moyen de commande est adapté 20 pour activer les vibreurs selon, par exemple, une méthode de "balayage glissant simultanée" du type décrit en référence à la figure 3 ou une méthode "à retard de temps variable" du type décrit en référence au
Tableau 2.
Quatre vibreurs sont montrés sur la figure 9(a), agencés en deux groupes contenant chacun deux vibreurs, mais un agencement de prospection sismique de l'invention n'est pas limité à ces nombres de vibreurs et de groupes. Un agencement de prospection sismique de l'invention destiné à être utilisé avec la méthode de "balayage glissant simultanée" peut contenir deux groupes ou plus de vibreurs activables indépendamment, un groupe au moins contenant plus d'un vibreur. Un agencement de prospection sismique de l'invention destiné à être utilisé avec la méthode de "retard de temps variable" peut contenir deux groupes ou plus de vibreurs activables indépendamment ou peut, en alternative, contenir deux vibreurs ou plus activables indépendamment.
L'agencement de prospection sismique de la figure 9(a) comprend en outre un arrangement d'un récepteur sismique ou plus (deux récepteurs 6, 7 sont montrés sur la figure 9(a), mais l'agencement de prospection sismique n'est pas limité à deux récepteurs). 10 La figure 9(b) est un organigramme schématique du moyen de commande 5. Le moyen de commande comprend un processeur de données programmable 8 avec une mémoire de programme 9, par exemple sous la forme d'une mémoire morte ROM, stockant un programme destiné à commander le moyen de commande 5 pour activer les vibreurs 1, 2, 3, 4 selon, par exemple, une méthode telle qu'illustrée sur la figure 3 ou dans le Tableau 2, ou telle que définie par l'équation (1) ci-dessus. Le système comprend en outre une mémoire de lecture/écriture non volatile 10 pour stocker, par exemple, toutes données devant être conservées en l'absence de source d'alimentation électrique. Une mémoire "de travail" ou "bloc-notes" pour le processeur de données est fournie par une mémoire vive (RAM) 11. Une interface d'entrée 12 est prévue, par exemple pour recevoir des instructions et des données. Une interface de sortie 13 est prévue, par exemple pour sortir des signaux d'activation jusqu'à un 25 récepteur sélectionné ou à un groupe de récepteurs sélectionné. Un programme définissant la séquence d'activation des récepteurs ou groupes de récepteurs peut être fourni par le biais de l'interface d'entrée 12 ou peut en alternative être fourni par une mémoire de données lisibles par machine 14.
Le programme destiné à faire fonctionner le moyen de commande et à exécuter la méthode décrite précédemment ici est stocké dans la mémoire de programme 9, laquelle peut être mise en oeuvre sous la forme d'une mémoire à semi-conducteurs, par exemple du type ROM, lequel est bien connu. Toutefois, le programme peut être stocké dans n'importe quel 5 autre support de stockage approprié, comme un porteur de données magnétique 9a (tel qu'une "disquette") ou un CD-ROM 9b.
La présente invention est applicable à la fois à des sources de vibration d'énergie sismique terrestres et à des sources de vibration d'énergie 10 sismique marines.
Claims (17)
1. Une méthode de prospection sismique comprenant les étapes consistant: à activer le ou chaque vibreur dans un premier groupe de vibreurs au temps TO; et à activer le ou chaque vibreur dans un deuxième groupe de vibreurs 10 au temps Tl, o TO < T1 < TO + S1 + L o Sl est le temps de balayage du premier groupe de vibreurs et L est le temps d'écoute; dans laquelle au moins un groupe parmi le premier groupe de vibreurs et le deuxième groupe de vibreurs comprend au moins 15 deux vibreurs.
2. Une méthode telle que revendiquée dans la revendication 1 dans laquelle le deuxième groupe de vibreurs est le premier groupe de vibreurs.
3. Une méthode de prospection sismique telle que revendiquée dans la revendication 1 ou la revendication 2 et comprenant l'étape consistant à activer le ou chaque groupe de vibreurs au moins autant de fois qu'il y a de vibreurs dans le groupe respectif et de 25 telle sorte que la contribution de chaque vibreur dans un groupe puisse être déterminée par une opération algébrique sur des enregistrements de récepteur de balayages réalisés par le groupe respectif.
4. Une méthode telle que revendiquée dans la revendication 1, la revendication 2 ou la revendication 3 et comprenant l'étape supplémentaire consistant à activer le ou chaque vibreur dans un troisième groupe de vibreurs au temps T2, o Tl < T2 < Tl + S2 + L o S2 est le temps de balayage du deuxième groupe de vibreurs.
5. Une méthode telle que revendiquée dans n'importe quelle revendication précédente dans laquelle Tl - TO > (n - 1) Si Io! n (fif o), o n est un entier naturel, Jo est la limite de fréquence inférieure du balayage de vibreur et f est la limite de fréquence supérieure du balayage de vibreur.
6. Une méthode telle que revendiquée dans la revendication 5 et comprenant de plus l'estimation à partir des données acquises du bruit dans les données acquises dû à la mIème harmonique, o mSn.
7. Une méthode telle que revendiquée dans n'importe quelle revendication précédente et comprenant de plus les étapes consistant à activer le ou chaque vibreur dans le premier groupe de vibreurs au temps T2, o Tl < T2 < Tl + S2 + L o S2 est le temps de balayage du deuxième groupe de vibreurs; et à activer le ou chaque vibreur dans le deuxième groupe de vibreurs 25 autempsT3oT2<T3<T2+S1 +L; dans laquelle le premier groupe de vibreurs est différent du deuxième groupe de vibreurs; et dans laquelle T3 - T2 É Tl - TO.
8. Une méthode telle que revendiquée dans la revendication 7 dans laquelle le jième groupe de vibreurs contient Nj vibreurs et est activé Ni + 1 fois.
9. Une méthode de prospection sismique comprenant les étapes consistant: à activer un premier vibreur au temps TO; à activer un deuxième vibreur différent du premier vibreur au temps T1, o TO < T1 < TO + Sl + L o Sl est le temps de balayage du 10 premier vibreur et L est le temps d'écoute; à activer le premier vibreur au temps T2, o Tl < T2 < Tl + S2 + L o S2 est le temps de balayage du deuxième vibreur; et à activer le deuxième vibreur au temps T3 o T2 < T3 < T2 + Si + L et o T3-T2ÉTl -TO.
10. Une méthode telle que revendiquée dans la revendication 7, la revendication 8 ou la revendication 9 dans laquelle T3 - T2 > Ti - TO. 20
11. Un agencement de prospection sismique comprenant: une pluralité de vibreurs; et un moyen de commande adapté pour activer le ou chaque vibreur dans un premier groupe de vibreurs au temps TO et pour activer le ou chaque vibreur dans un deuxième groupe de vibreurs au temps Tl, o TO < Tl < TO + Si + L o Si est le temps de balayage du premier groupe de vibreurs et L est le temps d'écoute. 30
12. Un agencement de prospection sismique tel que revendiqué dans la revendication 11 dans lequel le moyen de commande est adapté pour activer le ou chaque groupe de vibreurs au moins autant de fois qu'il y a de vibreurs dans le groupe respectif et de telle sorte que la contribution de chaque vibreur dans un groupe puisse être 5 déterminée par une opération algébrique sur des enregistrements de récepteur de balayages réalisés par le groupe respectif.
13. Un agencement de prospection sismique tel que revendiqué dans la revendication 11 ou la revendication 12 dans lequel le premier groupe de vibreurs est différent du deuxième groupe de vibreurs; et dans lequel le moyen de commande est de plus adapté pour activer le ou chaque vibreur dans le premier groupe de vibreurs au temps T2, o Tl < T2 < Tl + S2 + L o S2 est le temps de balayage du deuxième groupe de vibreurs et pour activer le ou chaque vibreur dans le deuxième groupe de vibreurs au temps T3 o T2 < T3<T2+S1 + LetoT3-T2*T1 -TO.
14. Un agencement de prospection sismique tel que revendiqué dans n'importe lesquelles des revendications 11, 12 et 13 et adapté pour 20 activer les premier et deuxième groupes de vibreurs de telle sorte que Tl - TO > (n - 1) Si fol n (f t - f o), o n est un entier naturel, fo est la limite de fréquence inférieure du balayage de vibreur et f est la limite de fréquence supérieure du balayage de vibreur.
15. Un agencement de prospection sismique comprenant: une pluralité de vibreurs; et un moyen de commande adapté pour: 30 a) activer un premier vibreur au temps TO; b) activer un deuxième vibreur différent du premier vibreur au temps Tl, o TO < T1 < TO + Sl + L o Sl est le temps de balayage du premier vibreur et L est le temps d'écoute; c) activer le premier vibreur au temps T2, o Tl < T2 < Tl + S2 + L o S2 est le temps de balayage du deuxième vibreur; et d) activer le deuxième vibreur au temps T3 o T2 < T3 < T2 + Si + L et o T3 - T2 É Tl - TO.
16. Un agencement de prospection sismique tel que revendiqué dans la revendication 15 et adapté pour activer le deuxième groupe de vibreurs au temps T3 o T3 - T2 > Tl - TO.
17. 20 18.
Un agencement de prospection sismique tel que revendiqué dans n'importe lesquelles des revendications i1 à 16 dans lequel le moyen de commande comprend un processeur de données programmable.
Un support de stockage contenant un programme pour le processeur de données d'un agencement de prospection sismique tel que défini dans la revendication 17.
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
GB0207995A GB2387226C (en) | 2002-04-06 | 2002-04-06 | A method of seismic surveying |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2853085A1 true FR2853085A1 (fr) | 2004-10-01 |
FR2853085B1 FR2853085B1 (fr) | 2010-10-15 |
Family
ID=9934389
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR0304298A Expired - Fee Related FR2853085B1 (fr) | 2002-04-06 | 2003-04-07 | Une methode de prospection sismique |
Country Status (5)
Country | Link |
---|---|
US (5) | US7050356B2 (fr) |
FR (1) | FR2853085B1 (fr) |
GB (1) | GB2387226C (fr) |
NL (2) | NL1023118C2 (fr) |
NO (2) | NO336129B1 (fr) |
Families Citing this family (101)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2387226C (en) * | 2002-04-06 | 2008-05-12 | Westerngeco Ltd | A method of seismic surveying |
WO2004095073A2 (fr) * | 2003-04-01 | 2004-11-04 | Exxonmobil Upstream Research Company | Source vibratoire haute frequence faconnee |
US20060164916A1 (en) * | 2003-08-11 | 2006-07-27 | Krohn Christine E | Method for continuous sweepting and separtion of multiple seismic vibrators |
GB2416033B (en) | 2004-07-10 | 2006-11-01 | Westerngeco Ltd | Seismic vibratory acquisition method and apparatus |
GB2425838B (en) * | 2005-05-03 | 2007-06-27 | Westerngeco Seismic Holdings | Source signature deconvolution method |
US7885143B2 (en) * | 2006-07-05 | 2011-02-08 | Westerngeco L.L.C. | Seismic acquisition system |
US7295490B1 (en) * | 2006-07-20 | 2007-11-13 | Conocophillips Company | System and method of phase encoding for high fidelity vibratory seismic data |
ES2652413T3 (es) | 2006-09-28 | 2018-02-02 | Exxonmobil Upstream Research Company | Inversión iterativa de datos a partir de fuentes geofísicas simultáneas |
DE602007004516D1 (de) * | 2006-11-01 | 2010-03-11 | Eveready Battery Inc | Alkali-batteriezelle mit verminderter gasung und verminderter entfärbung |
RU2450291C2 (ru) * | 2006-12-08 | 2012-05-10 | Конокофиллипс Компани | Способ выбора параметров динамических источников для получения данных сейсмических вибраторов |
US8000168B2 (en) | 2006-12-08 | 2011-08-16 | Conocophillips Company | Dynamic source parameter selection for seismic vibrator data acquisition |
CA2660228A1 (fr) * | 2007-02-08 | 2008-06-12 | Conocophillips Company | Selection dynamique de parametre source pour acquisition de donnees par vibreurs sismiques |
GB2447236B (en) * | 2007-03-09 | 2010-02-24 | Westerngeco Seismic Holdings | Method of estimating harmonic noise within slip-sweep Vibroseis signals |
WO2008123920A1 (fr) * | 2007-04-10 | 2008-10-16 | Exxonmobil Upstream Research Company | Séparation et suppression du bruit de fond de données sismiques de multiples sources vibratoires |
US7859945B2 (en) | 2007-07-06 | 2010-12-28 | Cggveritas Services Inc. | Efficient seismic data acquisition with source separation |
GB2451630B (en) | 2007-08-04 | 2009-12-09 | Westerngeco Seismic Holdings | Composite sweeps of high and low frequency part |
US7864630B2 (en) * | 2007-11-01 | 2011-01-04 | Conocophillips Company | Method and apparatus for minimizing interference between seismic systems |
CA2703588C (fr) * | 2007-12-12 | 2015-12-01 | Exxonmobil Upstream Research Company | Procede et appareil d'evaluation de formations sous-marines |
US20090168600A1 (en) * | 2007-12-26 | 2009-07-02 | Ian Moore | Separating seismic signals produced by interfering seismic sources |
EA017177B1 (ru) | 2008-03-21 | 2012-10-30 | Эксонмобил Апстрим Рисерч Компани | Эффективный способ инверсии геофизических данных |
US7916576B2 (en) * | 2008-07-16 | 2011-03-29 | Westerngeco L.L.C. | Optimizing a seismic survey for source separation |
SG193173A1 (en) * | 2008-08-11 | 2013-09-30 | Exxonmobil Upstream Res Co | Estimation of soil properties using waveforms of seismic surface waves |
US8077547B2 (en) * | 2008-09-26 | 2011-12-13 | Providence technologies, Inc. | Method and apparatus for seismic exploration |
US8681589B2 (en) | 2008-10-03 | 2014-03-25 | Shell Oil Company | Method and system for performing seismic surveys with a low frequency sweep |
US8982665B2 (en) * | 2008-10-08 | 2015-03-17 | Westerngeco L.L.C. | Dithered slip sweep vibroseis acquisition system and technique |
US9213119B2 (en) * | 2008-10-29 | 2015-12-15 | Conocophillips Company | Marine seismic acquisition |
US8467267B2 (en) * | 2008-10-29 | 2013-06-18 | Conocophillips Company | Asynchronous operation of seismic sources in a seismic survey |
US8174927B2 (en) | 2008-12-17 | 2012-05-08 | Westerngeco L.L.C. | Method for optimizing acoustic source array performance |
US9128207B2 (en) * | 2008-12-23 | 2015-09-08 | Westerngeco L.L.C. | Compensating seismic data for source variations |
US9052410B2 (en) * | 2009-02-12 | 2015-06-09 | Conocophillips Company | Multiple seismic signal inversion |
US8395966B2 (en) * | 2009-04-24 | 2013-03-12 | Westerngeco L.L.C. | Separating seismic signals produced by interfering seismic sources |
CA2775150C (fr) * | 2009-12-02 | 2015-05-19 | Conocophillips Company | Extraction d'enregistrements discrets d'enregistrements sismiques continus |
US8079440B2 (en) * | 2010-01-26 | 2011-12-20 | Westerngeco L.L.C. | Determining the polarity of movement of an actuator mass in a seismic vibrator |
US8537638B2 (en) * | 2010-02-10 | 2013-09-17 | Exxonmobil Upstream Research Company | Methods for subsurface parameter estimation in full wavefield inversion and reverse-time migration |
US8400873B2 (en) | 2010-02-22 | 2013-03-19 | Westerngeco L.L.C. | Vibroseis seismic acquisition technique |
US8223587B2 (en) * | 2010-03-29 | 2012-07-17 | Exxonmobil Upstream Research Company | Full wavefield inversion using time varying filters |
US8694299B2 (en) | 2010-05-07 | 2014-04-08 | Exxonmobil Upstream Research Company | Artifact reduction in iterative inversion of geophysical data |
US8756042B2 (en) | 2010-05-19 | 2014-06-17 | Exxonmobile Upstream Research Company | Method and system for checkpointing during simulations |
WO2011144215A2 (fr) | 2010-05-20 | 2011-11-24 | Entreprise Nationale De Geophysique - Enageo- Filiale Du Groupe Sonatrach | Méthode d'atténuation du bruit harmonique en vibrosismique par filtrage temps-variant avec référence |
US8588027B2 (en) | 2010-07-12 | 2013-11-19 | Westerngeco L.L.C. | Vibroseis acquisition technique and system using dynamic source control |
US8767508B2 (en) | 2010-08-18 | 2014-07-01 | Exxonmobil Upstream Research Company | Using seismic P and S arrivals to determine shallow velocity structure |
US8437998B2 (en) | 2010-09-27 | 2013-05-07 | Exxonmobil Upstream Research Company | Hybrid method for full waveform inversion using simultaneous and sequential source method |
SG188191A1 (en) | 2010-09-27 | 2013-04-30 | Exxonmobil Upstream Res Co | Simultaneous source encoding and source separation as a practical solution for full wavefield inversion |
US20120075955A1 (en) * | 2010-09-28 | 2012-03-29 | Timothy Dean | Efficient seismic source operation in connection with a seismic survey |
US8582395B2 (en) * | 2010-11-04 | 2013-11-12 | Westerngeco L.L.C. | Marine vibroseis motion correction |
CA2815054C (fr) | 2010-12-01 | 2017-05-16 | Exxonmobil Upstream Research Company | Inversion de source simultanee pour donnees de flute sismique marines presentant une fonction objective de correlation croisee |
US9134442B2 (en) * | 2010-12-16 | 2015-09-15 | Bp Corporation North America Inc. | Seismic acquisition using narrowband seismic sources |
US9551798B2 (en) * | 2011-01-21 | 2017-01-24 | Westerngeco L.L.C. | Seismic vibrator to produce a continuous signal |
WO2012134621A1 (fr) | 2011-03-30 | 2012-10-04 | Exxonmobil Upstream Research Company | Vitesse de convergence d'une inversion d'un champ d'onde complet utilisant une mise en forme spectrale |
SG193233A1 (en) | 2011-03-31 | 2013-10-30 | Exxonmobil Upstream Res Co | Method of wavelet estimation and multiple prediction in full wavefield inversion |
EP2944980A1 (fr) | 2011-05-13 | 2015-11-18 | Saudi Arabian Oil Company | Conception d'acquisition de source de fréquence couplée balayée en fonction de temps-distance et débruitage de données |
US9588241B2 (en) | 2011-05-13 | 2017-03-07 | Saudi Arabian Oil Company | Frequency-varying filtering of simultaneous source seismic data |
US9158019B2 (en) * | 2011-06-08 | 2015-10-13 | Westerngeco L.L.C. | Enhancing low frequency content in marine simultaneous vibroseis acquisition |
CA2839277C (fr) | 2011-09-02 | 2018-02-27 | Exxonmobil Upstream Research Company | Utilisation d'une projection sur des ensembles convexes pour limiter l'inversion du champ d'ondes complet |
US9176930B2 (en) | 2011-11-29 | 2015-11-03 | Exxonmobil Upstream Research Company | Methods for approximating hessian times vector operation in full wavefield inversion |
US9753163B2 (en) * | 2012-01-12 | 2017-09-05 | Westerngeco L.L.C. | Simultaneous marine vibrators |
US9348041B2 (en) * | 2012-02-15 | 2016-05-24 | Westerngeco L.L.C. | Phase modulation and noise minimization for simultaneous vibroseis acquisition |
US10012745B2 (en) | 2012-03-08 | 2018-07-03 | Exxonmobil Upstream Research Company | Orthogonal source and receiver encoding |
US9341724B2 (en) | 2012-03-22 | 2016-05-17 | Westerngeco L.L.C. | Sweep sequence determination for overlapping sweeps |
JP5965715B2 (ja) * | 2012-04-27 | 2016-08-10 | 古野電気株式会社 | 超音波送受信装置、超音波送受信方法、および超音波送受信プログラム |
US9348050B2 (en) * | 2012-05-23 | 2016-05-24 | Exxonmobil Upstream Research Company | Near-surface noise prediction and removal for data recorded with simultaneous seismic sources |
US20140036623A1 (en) * | 2012-07-31 | 2014-02-06 | Cggveritas Services Sa | Device and method for synchronized marine acquisition with reduced interference noise |
CA2892041C (fr) | 2012-11-28 | 2018-02-27 | Exxonmobil Upstream Research Company | Tomographie fondee sur le facteur de qualite q de donnees sismiques de reflexion |
US9857485B2 (en) | 2013-03-15 | 2018-01-02 | Westerngeco L.L.C. | Methods and systems for marine survey acquisition |
ES2838696T3 (es) | 2013-05-01 | 2021-07-02 | Cgg Services Sa | Aparato y método para la adquisición de datos sísmicos con activación simultánea de vibradores agrupados |
MX346526B (es) | 2013-05-24 | 2017-03-23 | Exxonmobil Upstream Res Co | Inversión multi-parámetro a través de fwi elástica dependiente de compensación. |
US10459117B2 (en) | 2013-06-03 | 2019-10-29 | Exxonmobil Upstream Research Company | Extended subspace method for cross-talk mitigation in multi-parameter inversion |
US10371839B2 (en) | 2013-06-05 | 2019-08-06 | Saudi Arabian Oil Company | Blended land seismic data acquisition employing dispersed source arrays with variable sweep length |
US9702998B2 (en) | 2013-07-08 | 2017-07-11 | Exxonmobil Upstream Research Company | Full-wavefield inversion of primaries and multiples in marine environment |
CN104375165B (zh) * | 2013-08-15 | 2017-02-15 | 中国石油天然气集团公司 | 一种可控震源分区同时扫描激发方法 |
EP3036566B1 (fr) | 2013-08-23 | 2018-04-04 | Exxonmobil Upstream Research Company | Recherche simultanée de sources pendant une acquisition sismique et une inversion sismique |
US10036818B2 (en) | 2013-09-06 | 2018-07-31 | Exxonmobil Upstream Research Company | Accelerating full wavefield inversion with nonstationary point-spread functions |
US9910189B2 (en) | 2014-04-09 | 2018-03-06 | Exxonmobil Upstream Research Company | Method for fast line search in frequency domain FWI |
SG11201608175SA (en) | 2014-05-09 | 2016-11-29 | Exxonmobil Upstream Res Co | Efficient line search methods for multi-parameter full wavefield inversion |
US10185046B2 (en) | 2014-06-09 | 2019-01-22 | Exxonmobil Upstream Research Company | Method for temporal dispersion correction for seismic simulation, RTM and FWI |
CN106662664A (zh) | 2014-06-17 | 2017-05-10 | 埃克森美孚上游研究公司 | 快速粘声波和粘弹性全波场反演 |
US10838092B2 (en) | 2014-07-24 | 2020-11-17 | Exxonmobil Upstream Research Company | Estimating multiple subsurface parameters by cascaded inversion of wavefield components |
US10422899B2 (en) | 2014-07-30 | 2019-09-24 | Exxonmobil Upstream Research Company | Harmonic encoding for FWI |
WO2016044538A1 (fr) * | 2014-09-19 | 2016-03-24 | Conocophillips Company | Extension de bande passante au-delà du signal de balayage de vibreur par l'intermédiaire d'une inversion contrainte de vibreur multiple simultanée |
US10386511B2 (en) | 2014-10-03 | 2019-08-20 | Exxonmobil Upstream Research Company | Seismic survey design using full wavefield inversion |
US9977141B2 (en) | 2014-10-20 | 2018-05-22 | Exxonmobil Upstream Research Company | Velocity tomography using property scans |
EP3234659A1 (fr) | 2014-12-18 | 2017-10-25 | Exxonmobil Upstream Research Company | Planification évolutive de travaux sismiques itératifs parallèles |
US10520618B2 (en) | 2015-02-04 | 2019-12-31 | ExxohnMobil Upstream Research Company | Poynting vector minimal reflection boundary conditions |
SG11201704620WA (en) | 2015-02-13 | 2017-09-28 | Exxonmobil Upstream Res Co | Efficient and stable absorbing boundary condition in finite-difference calculations |
CN107407736B (zh) | 2015-02-17 | 2019-11-12 | 埃克森美孚上游研究公司 | 生成无多次波的数据集的多阶段全波场反演处理 |
WO2016195774A1 (fr) | 2015-06-04 | 2016-12-08 | Exxonmobil Upstream Research Company | Procédé pour générer des images sismiques sans réflexions multiples |
US10838093B2 (en) | 2015-07-02 | 2020-11-17 | Exxonmobil Upstream Research Company | Krylov-space-based quasi-newton preconditioner for full-wavefield inversion |
CN108139499B (zh) | 2015-10-02 | 2020-02-14 | 埃克森美孚上游研究公司 | Q-补偿的全波场反演 |
WO2017065889A1 (fr) | 2015-10-15 | 2017-04-20 | Exxonmobil Upstream Research Company | Empilements d'angles de domaine de modèle d'inversion des formes d'onde complète (fwi) à préservation d'amplitude |
US10564273B2 (en) * | 2015-10-21 | 2020-02-18 | Canon Medical Systems Corporation | Ultrasonic diagnostic apparatus |
US10590758B2 (en) | 2015-11-12 | 2020-03-17 | Schlumberger Technology Corporation | Noise reduction for tubewave measurements |
EP3394642B1 (fr) | 2015-12-22 | 2022-02-16 | Shell Internationale Research Maatschappij B.V. | Procédé et système permettant de générer un relevé sismique |
US10768324B2 (en) | 2016-05-19 | 2020-09-08 | Exxonmobil Upstream Research Company | Method to predict pore pressure and seal integrity using full wavefield inversion |
CN109564296B (zh) | 2016-07-01 | 2021-03-05 | 斯伦贝谢技术有限公司 | 用于检测反射液压信号的井中对象的方法和系统 |
DE102017121346A1 (de) | 2016-09-15 | 2018-03-15 | Osram Opto Semiconductors Gmbh | Messsystem, Verwendung zumindest einer individuell betreibbaren Leuchtdioden-Leuchteinheit als Sendereinheit in einem Messsystem, Verfahren zum Betrieb eines Messsystems und Beleuchtungsquelle mit einem Messsystem |
US10288755B2 (en) | 2017-03-28 | 2019-05-14 | Saudi Arabian Oil Company | Seismic processing workflow for broadband single-sensor single-source land seismic data |
CN107290043B (zh) * | 2017-06-15 | 2023-07-28 | 贵州电网有限责任公司电力科学研究院 | 一种输电线路振动次数在线分布式监测方法 |
US11294083B2 (en) | 2017-08-16 | 2022-04-05 | Pgs Geophysical As | Geophysical survey techniques using different correlation intervals |
US20190317232A1 (en) * | 2018-04-17 | 2019-10-17 | Cgg Services Sas | High-productivity seismic data acquisition using calendar-time-based sweep initiation |
GB2578123B8 (en) * | 2018-10-16 | 2021-08-11 | Darkvision Tech Inc | Overlapped scheduling and sorting for acoustic transducer pulses |
US11703607B2 (en) | 2020-06-15 | 2023-07-18 | Saudi Arabian Oil Company | Determining a seismic quality factor for subsurface formations from a seismic source to a first VSP downhole receiver |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4953657A (en) * | 1987-11-30 | 1990-09-04 | Halliburton Geophysical Services, Inc. | Time delay source coding |
US5410517A (en) * | 1994-05-13 | 1995-04-25 | Exxon Production Research Company | Method for cascading sweeps for a seismic vibrator |
CA2316665A1 (fr) * | 1999-09-10 | 2001-03-10 | Thomas John Fleure | Methode de reduction de l'interference harmonique lors de l'utilisation de techniques d'enregistrement sismique ponctuel de sources chevauchantes |
US6366857B1 (en) * | 1999-06-25 | 2002-04-02 | Trimble Navigation Limited | Noise estimator for seismic exploration |
Family Cites Families (35)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US161379A (en) * | 1875-03-30 | Improvement in machines for making horseshoe-nails | ||
US2348003A (en) * | 1941-01-28 | 1944-05-02 | Gen Electric | Magnetic core |
US2359363A (en) * | 1942-04-11 | 1944-10-03 | American Cyanamid Co | Preparation of sulphonyl guanidines |
GB1423366A (en) * | 1972-07-21 | 1976-02-04 | Seiscom Ltd | Broad line seismic profiling using simultaneously radiating sources |
US4201972A (en) * | 1972-11-07 | 1980-05-06 | United Geophysical Corporation | Seismic prospecting system |
US4004267A (en) * | 1972-11-28 | 1977-01-18 | Geosource Inc. | Discrete frequency seismic exploration using non uniform frequency spectra |
US3815704A (en) * | 1972-12-26 | 1974-06-11 | Texaco Inc | Method for determining optimum seismic pulse |
US3895343A (en) * | 1973-05-07 | 1975-07-15 | Amoco Prod Co | Apparatus for producing adaptive pilot signals |
US4168485A (en) * | 1974-08-12 | 1979-09-18 | Continental Oil Company | Simultaneous use of pseudo-random control signals in vibrational exploration methods |
GB1602117A (en) * | 1977-10-04 | 1981-11-04 | Seismograph Service Englang Lt | Vibration generators particularly for seismic exploration |
US4519053A (en) * | 1981-11-27 | 1985-05-21 | Texas Instruments Incorporated | Force or pressure feedback control for seismic vibrators |
US4823326A (en) * | 1986-07-21 | 1989-04-18 | The Standard Oil Company | Seismic data acquisition technique having superposed signals |
US4750157A (en) * | 1987-05-06 | 1988-06-07 | Standard Oil Production Company | Seismic vibrator earth impedance determination and compensation system |
US4982374A (en) * | 1989-10-23 | 1991-01-01 | Halliburton Geophysical Services, Inc. | Method of source coding and harmonic cancellation for vibrational geophysical survey sources |
US5550786A (en) * | 1995-05-05 | 1996-08-27 | Mobil Oil Corporation | High fidelity vibratory source seismic method |
US5719821A (en) | 1995-09-29 | 1998-02-17 | Atlantic Richfield Company | Method and apparatus for source separation of seismic vibratory signals |
US5721710A (en) * | 1995-09-29 | 1998-02-24 | Atlantic Richfield Company | High fidelity vibratory source seismic method with source separation |
GB2306219B (en) * | 1995-10-12 | 1999-06-23 | Nigel Allister Anstey | 3-d seismic survey using multiple sources simultaneously |
US6181646B1 (en) * | 1997-01-07 | 2001-01-30 | Hyroacoustics, Inc. | Geophysical exploration system using seismic vibrator source which provides a composite sweep |
WO1999026179A1 (fr) * | 1997-11-14 | 1999-05-27 | Western Atlas International, Inc. | Acquisition et traitement de donnees sismiques par utilisation d'une distorsion non lineaire dans un signal de force de sol |
US6035257A (en) * | 1997-12-10 | 2000-03-07 | Pelton Company | Method and apparatus for reducing harmonic distortion |
US6152256A (en) * | 1998-10-29 | 2000-11-28 | Aspect Resources Llc | Method for controlling seismic vibrator ground forces |
GB2348003B (en) | 1999-03-19 | 2001-02-07 | Geco Prakla | Seismic data processing method for data acquired using overlapping vibratory sweeps |
GB9927395D0 (en) * | 1999-05-19 | 2000-01-19 | Schlumberger Holdings | Improved seismic data acquisition method |
GB2359363B (en) | 2000-02-15 | 2002-04-03 | Geco Prakla | Processing simultaneous vibratory seismic data |
CN1254693C (zh) * | 2000-10-17 | 2006-05-03 | 维斯特恩格科地震控股有限公司 | 对源编码和谐波消除采用级联扫描的方法 |
FR2818753B1 (fr) * | 2000-12-21 | 2003-03-21 | Inst Francais Du Petrole | Methode et dispositif de prospection sismique par emission simultanee de signaux sismisques obtenus en codant un signal par des sequences pseudo aleatoires |
US6842701B2 (en) * | 2002-02-25 | 2005-01-11 | Westerngeco L.L.C. | Method of noise removal for cascaded sweep data |
US6603707B1 (en) * | 2002-02-26 | 2003-08-05 | Compagnie Generale De Geophysique | Method of reducing harmonic noise in vibroseis operations |
FR2836723B1 (fr) * | 2002-03-01 | 2004-09-03 | Inst Francais Du Petrole | Methode et dispositif de prospection sismique par emission simultanee de signaux sismiques a base de sequences pseudo aleatoires |
GB2387226C (en) | 2002-04-06 | 2008-05-12 | Westerngeco Ltd | A method of seismic surveying |
US6942059B2 (en) * | 2002-11-13 | 2005-09-13 | Westerngeco, L.L.C. | Composite bandwidth marine vibroseis array |
WO2004095073A2 (fr) * | 2003-04-01 | 2004-11-04 | Exxonmobil Upstream Research Company | Source vibratoire haute frequence faconnee |
GB2416033B (en) * | 2004-07-10 | 2006-11-01 | Westerngeco Ltd | Seismic vibratory acquisition method and apparatus |
US7327633B2 (en) * | 2005-12-12 | 2008-02-05 | Westerneco L.L.C. | Systems and methods for enhancing low-frequency content in vibroseis acquisition |
-
2002
- 2002-04-06 GB GB0207995A patent/GB2387226C/en not_active Expired - Fee Related
-
2003
- 2003-03-31 US US10/403,135 patent/US7050356B2/en not_active Expired - Fee Related
- 2003-04-04 NO NO20031535A patent/NO336129B1/no not_active IP Right Cessation
- 2003-04-07 NL NL1023118A patent/NL1023118C2/nl not_active IP Right Cessation
- 2003-04-07 FR FR0304298A patent/FR2853085B1/fr not_active Expired - Fee Related
-
2005
- 2005-12-12 US US11/299,410 patent/US7672194B2/en not_active Expired - Fee Related
-
2006
- 2006-10-11 US US11/548,521 patent/US7376046B2/en not_active Expired - Lifetime
- 2006-12-29 NL NL1033146A patent/NL1033146C2/nl not_active IP Right Cessation
-
2007
- 2007-11-26 US US11/944,690 patent/US7602670B2/en not_active Expired - Fee Related
-
2009
- 2009-09-01 US US12/551,897 patent/US20100008187A1/en not_active Abandoned
-
2015
- 2015-01-08 NO NO20150043A patent/NO339176B1/no not_active IP Right Cessation
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4953657A (en) * | 1987-11-30 | 1990-09-04 | Halliburton Geophysical Services, Inc. | Time delay source coding |
US5410517A (en) * | 1994-05-13 | 1995-04-25 | Exxon Production Research Company | Method for cascading sweeps for a seismic vibrator |
US6366857B1 (en) * | 1999-06-25 | 2002-04-02 | Trimble Navigation Limited | Noise estimator for seismic exploration |
CA2316665A1 (fr) * | 1999-09-10 | 2001-03-10 | Thomas John Fleure | Methode de reduction de l'interference harmonique lors de l'utilisation de techniques d'enregistrement sismique ponctuel de sources chevauchantes |
Also Published As
Publication number | Publication date |
---|---|
NO20031535L (no) | 2003-10-07 |
GB2387226C (en) | 2008-05-12 |
GB2387226B (en) | 2005-08-10 |
US20060158962A1 (en) | 2006-07-20 |
GB2387226A (en) | 2003-10-08 |
GB0207995D0 (en) | 2002-05-15 |
US20030210609A1 (en) | 2003-11-13 |
NO20031535D0 (no) | 2003-04-04 |
NO336129B1 (no) | 2015-05-18 |
NO20150043L (no) | 2003-10-07 |
US20070091721A1 (en) | 2007-04-26 |
US7376046B2 (en) | 2008-05-20 |
US7672194B2 (en) | 2010-03-02 |
FR2853085B1 (fr) | 2010-10-15 |
NL1023118A1 (nl) | 2003-10-07 |
NL1023118C2 (nl) | 2007-01-12 |
US7050356B2 (en) | 2006-05-23 |
US7602670B2 (en) | 2009-10-13 |
NL1033146A1 (nl) | 2007-01-23 |
US20100008187A1 (en) | 2010-01-14 |
US20080232194A1 (en) | 2008-09-25 |
NL1033146C2 (nl) | 2009-12-04 |
NO339176B1 (no) | 2016-11-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
FR2853085A1 (fr) | Une methode de prospection sismique | |
FR2872922A1 (fr) | Procede et appareil d'acquisition vribratoire sismique | |
CA2520640C (fr) | Source vibratoire haute frequence faconnee | |
FR2719908A1 (fr) | Procédé de génération de données sismiques par la mise en cascade de balayages pour un vibrateur sismique. | |
FR2955396A1 (fr) | Dispositif de traitement de donnees sismiques marines | |
EP1173781A1 (fr) | Methode de surveillance sismique d'une zone souterraine par utilisation simultanee de plusieurs sources vibrosismiques | |
FR2478324A1 (fr) | Procede de prospection et d'exploration sismiques | |
FR2876458A1 (fr) | Perfectionnement aux traitements sismiques pour la suppression des reflexions multiples | |
CA1325840C (fr) | Procede d'acquisition de donnees sismiques et dispositif de mise en oeuvre du procede | |
FR2782389A1 (fr) | Procede et dispositif pour corriger les effets des deplacements du navire dans les mesures de sismique maritime | |
EP1552324B1 (fr) | Méthode de séparation de signaux de microsismicité de signaux sismiques émis par une ou plusieurs sources | |
FR2874270A1 (fr) | Procede d'exploration sismique | |
FR2961316A1 (fr) | Procede de traitement de donnees sismiques marines | |
FR2707115A1 (fr) | Procédé pour produire des formes d'ondelettes sismiques ayant une énergie de lobes latéraux minimale. | |
Wu et al. | Shot repetition: an alternative approach to blending in marine seismic | |
FR2985832A1 (fr) | Dispositif et procede pour le retrait de multiples d'epoques de donnees sismiques | |
FR2595831A1 (fr) | Procede pour retablir l'aspect de fond de donnees sismiques filtrees de facon bidimensionnelle | |
FR2972057A1 (fr) | Dispositif et procede pour un retrait de bruit commande par coherence multidimensionnelle des donnees | |
FR2872919A1 (fr) | Traitement coherent rapide pour codes a spectre de raies periodiques | |
EP4059180B1 (fr) | Procédé de synchronisation d'un système de communication basé sur la retransmission de données | |
EP0204599A1 (fr) | Procédé et dispositif pour la détermination de la signature d'émission à grande distance d'un ensemble d'émission sismique | |
FR2502794A1 (fr) | Procede d'exploration sismique d'un milieu notamment de prospection geophysique par ondes sismiques | |
WO2014191650A1 (fr) | Procedes et dispositifs d'acquisition de donnees sismiques relatives a une zone du sous-sol situee sous la mer | |
EP3459177A1 (fr) | Procédé de traitement d'un signal formé d'une séquence d'impulsions | |
EP2884456A1 (fr) | Procédé de réstitution du mouvement haute fréquence de la ligne de visée d'un imageur |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PLFP | Fee payment |
Year of fee payment: 14 |
|
ST | Notification of lapse |
Effective date: 20171229 |