FR 3059786 A1 20180608 FR 1757024 A 20170725 1 Procédé d'amélioration des acquisitions sismiques mettant en oeuvre des systèmes actifs ultralégers de détection sismique La présente invention concerne d'une manière générale l'acquisition de 5 données sismiques dans les domaines terrestre et marin. Il est connu de l'art antérieur que l'acquisition de données sismiques permet de générer une image de la structure géophysique du sous-sol. Bien que cette image ne permette pas d'indiquer précisément la localisation des ressources naturelles, telles que du pétrole ou du gaz, elle permet néanmoins, 10 pour une personne du métier, de déterminer la présence ou l'absence de ces ressources. Ainsi, la production en continue de telles images fait partie intégrante du processus d'exploration des sous-sols dans le but de découvrir ces ressources naturelles. Afin de produire ces images, il est connu de l'art antérieur plusieurs méthodes. Pour une exploration en mer, il est possible de 15 mettre en oeuvre un bateau tirant des flûtes sismiques ou streamers en anglais. Ces streamers peuvent être disposés horizontalement, à une profondeur constante par rapport au fond de la mer. Le bateau tire également un ensemble de sources sismiques apte à générer des ondes sismiques. Ces ondes sismiques se propagent vers le fond de la mer et le traverse pour le 20 pénétrer jusqu'à rencontrer une structure réfléchissante pour ainsi les réfléchir. Ces ondes sismiques réfléchies se propagent vers la surface de la mer jusqu'à être détectées par des capteurs sismiques intégrés aux streamers. En se basant sur ces données sismiques, il est alors possible de générer une image du sous-sol. S'il est utilisé un seul streamer, on parle alors d'une image 25 2D. S'il est utilisé plusieurs streamers tirés simultanément par le bateau, on parle d'une image 3D. Une image 4D est obtenue lorsque deux acquisitions sismiques 3D sont effectuées dans la même zone, mais avec un décalage dans le temps. La première acquisition est effectuée à un temps tO et la seconde à un temps tO+t, t pouvant être égal à plusieurs mois ou plusieurs années par 30 exemple. De manière alternative, il est également possible d'enregistrer des données sismiques en utilisant des câbles ou des capteurs posés sur le fond de mer, ces techniques s'appellent ocean bottom cables (OBC), ocean bottom 3059786 2 nodes (OBN) ou ocean bottom seismometers (OBS). Par exemple, dans le cas d'une acquisition OBC, une pluralité de capteurs sont connectés entre eux dans un câble. Plusieurs de ces câbles sont ensuite posés sur le fond de la mer, dans la zone à étudier. Il est également possible d'intégrer un ou plusieurs capteurs dans un véhicule aquatique sous-marin. Ce véhicule, éventuellement autonome, se pose sur le fond de la mer afin d'enregistrer des données sismiques. Il est ensuite récupéré par un bateau pour y transférer ces données sismiques. Dans le cas d'une acquisition sismique terrestre, un système dédié à une acquisition 4D comprend plusieurs capteurs, par exemple des hydrophones ou des accéléromètres, positionnés sur une zone où l'on veut explorer le sous-sol. Ces capteurs sont en contact avec le sol. Plusieurs sources sismiques sont également placées sur le sol, dans la zone à étudier. Un dispositif d'enregistrement est connecté aux capteurs et localisé, par exemple, dans un camion. Chaque source sismique comprend un nombre variable de vibrateurs, typiquement entre un et cinq, et peut comprendre également un contrôleur local. Un contrôleur central peut être présent afin de coordonner les temps de tir des sources sismiques. Un système GPS peut être mis en oeuvre afin de corréler dans le temps le tir des sources et les données acquises par les capteurs. Dans cette configuration, les sources sismiques sont contrôlées pour générer des ondes sismiques, et la pluralité de capteurs enregistre les ondes réfléchies par des réservoirs de pétrole ou de gaz, ou d'autres structures. L'étude sismique peut être répétée à divers moments dans le temps, par exemple tous les mois ou tous les ans, afin d'imager de nouveau le sous-sol de manière à déterminer les changements persistants sur la période dans les réservoirs. Toutes ces techniques peuvent être utilisées pour surveiller un réservoir de production. Pour ces configurations, le but du traitement 4D est de déterminer comment et où les propriétés pétrophysiques changent en évaluant les différences dans les données sismiques acquises à divers moments, en général avant l'exploitation d'un réservoir, il s'agit de l'étude de référence (baseline survey), et après l'exploitation de ce même réservoir, il s'agit de l'étude de surveillance (monitor survey). Actuellement, les solutions 3059786 3 d'étude 4D ont pour objet la mise à jour complète d'une acquisition 3D du sous-sol pour une zone considérée. En d'autres termes, la baseline survey et la monitor survey acquièrent un volume important de données sismiques afin d'obtenir une image 3D complète. Par conséquent, les techniques actuelles 5 prennent beaucoup de temps, de l'ordre de plusieurs semaines à plusieurs mois, sont relativement peu précises d'un point de vue de la répétabilité, car utilisant des systèmes non permanent, et génèrent des coûts importants. Tous ces facteurs représentent une barrière au développement des études 4D. Par conséquent, il existe un besoin pour rendre moins coûteuses les études 4D, 10 mais également les rendre plus rapides et plus précises. Ainsi, l'invention a notamment pour but d'améliorer la sensibilité et par conséquent le niveau de détection des systèmes actifs sismiques ultra légers. Il est décrit sur la figure 1 un exemple de système actif de détection sismique ultraléger. Le terme ultraléger veut dire qu'il est mis en oeuvre un 15 nombre limité de couples source-récepteur par rapport à une acquisition sismique classique pour laquelle plusieurs centaines à plusieurs millions de couples source-récepteur sont utilisés. On peut voir ici un seul couple source 1 et récepteur 2. Il peut également être mis en oeuvre une antenne de sources et une antenne de récepteurs. Une antenne de sources, respectivement de 20 récepteurs, comporte une agglomération de sources, respectivement de récepteurs. Ils sont, par exemple, assemblés par soudage. L'antenne de sources se comporte comme s'il s'agissait d'une source unique et l'antenne de récepteurs comme un récepteur unique. Une antenne de sources permet de focaliser l'émission des ondes et une antenne de récepteurs permet de 25 focaliser la réception de ces ondes. Cette technique est appelée beam steering en anglais. En outre, cela permet de créer un filtre de bruit. Il est également représenté l'aquifère 3, la faille 4, le réservoir de gaz 6, la zone d'intérêt ou spot 7 et le lieu de production 5. La limitation du nombre de couples source-récepteur rend les méthodes spatiales classiques de suppression de bruit peu 30 applicables. Ils peuvent parfois masquer des effets de production faibles, à savoir une modification de moins de 1% des propriétés sismiques dues aux variations des propriétés physiques du spot sur l'échelle de temps considérée. Cela est notamment vrai sur des échelles de temps très courtes, par exemple 3059786 4 environ 3 mois, une journée ou une heure, pendant lesquelles les effets sismiques induits par la production sont faibles. La détection de ces variations permet d'anticiper un effet à venir, comme par exemple un changement de saturation, qui sont normalement visibles sur des périodes plus longues ainsi 5 que leurs effets associés. Les effets, sur une période courte, traduisent l'excitation précoce de la zone d'intérêt sous l'effet d'un changement. Ainsi, il existe une nécessité d'augmenter la sensibilité de détection sur de courtes périodes pour les systèmes actifs sismiques ultralégers. Afin de répondre à ce problème, l'invention a pour objet un procédé de 10 détermination du positionnement optimal de couples source-récepteur aptes à acquérir des données sismiques comprenant les étapes suivantes : une première étape d'identification d'une zone d'intérêt ayant fait l'objet d'une acquisition sismique antérieure afin d'obtenir une image de son sous- 15 sol; une deuxième étape d'obtention des données sismiques acquises lors de l'acquisition sismique antérieure de ladite zone d'intérêt pendant un temps d'intérêt; 20 une troisième étape d'application d'une démigration auxdites données sismiques afin de déterminer les positions de chaque couple source-récepteur ayant contribué à l'image dudit sous-sol de ladite zone d'intérêt pendant ledit temps d'intérêt; 25 une quatrième étape d'obtention des traces brutes pour lesdites positions des couples source-récepteur; une cinquième étape de sélection d'au moins une trace brute optimale parmi lesdites traces brutes ; 3059786 5 - une sixième étape de détermination des positions optimales des couples source-récepteur correspondant à la dite trace brute optimale; De préférence, la démigration de la troisième étape est effectuée par 5 gamme d'octave. De préférence, l'acquisition sismique antérieure est une acquisition 4D. De préférence, la cinquième étape comporte en outre une sélection des traces brutes optimales permettant de détecter les meilleurs effets 4D. De préférence, la cinquième étape comporte en outre l'utilisation de 10 modèles pétro-élastiques afin de sélectionner les traces brutes optimales. De préférence, la cinquième étape comporte en outre une sélection tenant compte des facteurs de surface. De préférence, la cinquième étape comporte en outre une sélection tenant compte d'une zone de référence servant d'étalonnage de variations 15 sismiques. De préférence, la troisième étape comporte en outre l'obtention de la pente et de la direction d'émission et de réception optimales des ondes sismiques pour chaque couple source-récepteur. De préférence, le procédé comprend en outre une septième étape de 20 placement de couples source-récepteur auxdites positions optimales et une huitième étape d'acquisition de nouvelles données sismiques auxdites positions optimales. De préférence, la huitième étape se base sur la pente et la direction d'émission et de réception optimales des ondes sismiques. 25 L'invention sera mieux comprise à la lecture de la description qui suit, faite en référence aux figures annexées, dans lesquelles : - la figure 1 illustre un exemple de système actif de détection sismique ultraléger ; 30 - la figure 2 illustre un mode de réalisation de la présente invention ; - la figure 3 illustre un autre exemple de mise en oeuvre de l'invention ; 3059786 6 - la figure 4 illustre un exemple de mise en oeuvre du procédé de sélection d'au moins une trace optimale objet de la présente invention ; - la figure 5 illustre une onde sismique traversant une zone de référence. 5 - la figure 6 illustre un exemple de traitement d'une trace brute sismique par le procédé mis en oeuvre dans la présente invention. Il est connu de l'état de la technique l'utilisation d'agents chimiques, tels que des polymères, lors de l'extraction de fluides situés dans un sous-sol. 10 Ils permettent l'amélioration de rendement du taux de récupération supérieur à 10%. L'utilisation de ces agents chimiques est encore plus justifiée dans les cas de récupération secondaire de fluide. Ces polymères sont injectés en amont d'une phase de récupération secondaire, c'est-à-dire une nouvelle injection, de manière à créer un moyen de poussée, ayant des effets similaires 15 à celui d'un piston, et qui sera ensuite poussé. Par ailleurs, dans une majorité de cas, une sismique 2D ou 3D a déjà été acquise et traitée. Cela permet d'avoir accès, notamment, à des informations telles qu'un modèle de vitesse, une image sismique migrée, un modèle de multiples, des traces sismiques brutes et prétraitées ainsi que la localisation exacte des sources et récepteurs, 20 notamment utilisées pour l'image existante 2D ou 3D, et un modèle de statique pour les données terrestres. La connaissance et l'interprétation du sous-sol permettent d'obtenir les informations telles que les propriétés élastiques des roches de la cible et l'effet de production sismique attendu, c'est-à-dire la variation des propriétés élastiques sous l'effet de production. La 25 connaissance du terrain permet également de compléter notamment les informations telles qu'une obstruction de surface (route, pipeline, puits, plateforme, usine, habitation, etc.) ou un générateur de bruits sismiques non répétables. L'invention a notamment pour objet l'incorporation d'un agent de 30 contraste sismique dans les injections de fluides, pouvant être liquide ou gazeux afin d'améliorer les acquisitions sismiques mettant en oeuvre des systèmes ultralégers de détection sismique. La plupart des systèmes ultralégers de détection sismique visent à détecter les changements liés à la 3059786 7 production d'un fluide ou d'un gaz. Des agents chimiques, tels que des polymères, sont utilisés pour améliorer le taux de récupération du pétrole à extraire. Un des modes de réalisation de l'invention consiste à injecter, en complément des agents chimiques et/ou du fluide d'injection, un agent de 5 contraste sismique amplifiant ou diminuant de façon sensible la réponse sismique. Cela permet d'augmenter la réponse sismique de l'injection et donc de pouvoir détecter, le plus rapidement possible, l'arrivée du front du fluide injecté dans la zone surveillée. Les agents chimiques étant la plupart du temps à l'avant du front et formant ainsi ce qui peut être appelé un piston chimique. 10 Un tel agent chimique, liquide ou gazeux, est injecté à l'intérieur du piston chimique, à l'avant ou l'arrière de celui-ci. De préférence, l'agent de contraste sismique doit être neutre vis-à-vis du piston chimique, c'est-à-dire qu'il ne doit pas réduire les propriétés de l'agent chimique utilisé. L'agent de contraste sismique peut comprendre des nanoparticules à forte densité ou à faible 15 densité et/ou des nanoparticules à fréquence de résonance compatible avec les acquisitions sismiques. L'agent de contraste sismique peut également comprendre un gaz neutre tel que l'azote. Il peut également comprendre un polymère. Il est montré sur la figure 2 une nappe de pétrole 21 dans laquelle est injectée, par l'intermédiaire d'un injecteur 24, de l'eau 25 et des polymères 20 22. Il est également incorporé au milieu du piston de polymères 22 un agent de contraste sismique 23. L'invention a aussi pour objet un procédé permettant l'amélioration du positionnement des couples source-récepteur et permettant également le traitement sismique d'une trace brute unique, à savoir un enregistrement en 25 sortie de capteur, ou moyenné dans le temps, c'est-à-dire une addition, ou un stack en anglais, dans le temps des traces brutes afin de diminuer le bruit non répétable ou d'augmenter la force du signal, dans le cadre d'un système permanent ou semi-permanent. Dans un mode de réalisation de l'invention, il est identifié la zone soumise à la détection afin de recueillir les informations 30 disponibles pour mettre en oeuvre le procédé. La première étape est l'utilisation des données sismiques existantes dans ou près de la zone d'intérêt, appelée également le spot. Ces données sismiques 2D, 3D, 4D ou de puits sont utilisées comme source d'information pour connaître le trajet des 3059786 8 ondes sismiques dans le sous-sol étudié et en particulier dans la zone à détecter. De préférence, on regarde les points de tir centrés sur la zone à détecter, une première pré-sélection est effectuée en éliminant les offsets, c'est-à-dire les distances entre le point de tir et le récepteur, qui sont trop 5 bruités pour une détection efficace, c'est-à-dire ayant un rapport signal à bruit trop faible. Une démigration, c'est-à-dire un tracé inverse de rayons sismiques, totale ou partielle des données sismiques disponibles est ensuite effectuée. Ce tracé de rayon inverse à partir de la zone d'intérêt permet de retrouver les positions X, Y et Z théoriques de chaque source et récepteur qui 10 ont contribué à l'imagerie du spot d'intérêt, mais également le temps où la zone d'intérêt apparaît sur les traces brutes. Ce temps est appelé temps d'intérêt. Cette étape permet également de déterminer les pentes d'arrivée des ondes sismiques ayant illuminé la zone d'intérêt. La figure 3 est une illustration de la mise en oeuvre de ce procédé. On peut y voir les sources 31 15 et les récepteurs 34, la zone d'intérêt ou spot 32 et les tracés de rayon inverse 33 lors de l'étape de démigration. La démigration peut également être effectuée par gamme d'octave, lorsque, par exemple la géologie est complexe et nécessite une analyse plus fine. Il est ensuite mis en oeuvre une sélection d'au moins une trace 20 optimale. Le procédé s'applique au sein de la zone contribuant à l'imagerie de la zone d'intérêt et comprend une étape de sélection parmi trois paramètres. Le premier paramètre est la contribution effective, c'est-à-dire l'énergie ou amplitude de l'enregistrement brut, des couples source-récepteur qui a permis la construction de l'image de la zone d'intérêt dans les données sismiques 2D, 25 3D ou de puits existantes, mais aussi le rapport signal à bruit global de la trace. En effet, on ne regarde pas seulement le comportement de la trace dans la zone d'intérêt, mais également en dehors de la zone d'intérêt afin d'y détecter des rapports signal à bruit faibles. Une trace peut présenter un rapport signal à bruit élevé dans la zone d'intérêt, mais faible en dehors de 30 celle-ci. De préférence, on applique un second paramètre concernant les informations dynamiques à savoir les données sismiques 4D existantes et/ou les effets de production attendus et la modélisation de leurs impacts sur la ou les traces sismiques. De préférence, on applique un troisième paramètre 3059786 9 concernant la prise en compte des informations de surface, comme par exemple les obstructions qui empêchent un placement de sources ou de récepteurs, les forts générateurs de bruits qu'il faut éviter, la possibilité d'accès à une source de courant ou la couverture internet qu'il faut privilégier. 5 Dans le cas où plusieurs placements optimums ont été identifiés par la prise en compte des deux premiers paramètres, la sélection finale se fera sur les critères de surface. Le premier paramètre est obtenu en recoupant les positions X, Y et Z théoriques des sources et récepteurs lors de l'étape de démigration avec la 10 grille réelle d'acquisition des données sismiques existantes, on récupère les traces brutes c'est-à-dire les mesures effectives du sous-sol depuis ces positions. A partir de ces traces brutes et de l'information de temps pendant lequel la zone d'intérêt apparaît sur les traces brutes, on peut effectuer une sélection qualitative et quantitative sur un sous ensemble de traces 15 potentiellement optimales. Une trace est optimale lorsqu'elle permet la détection d'un changement. L'information de temps est obtenue à l'étape de démigration et est appelée également temps d'intérêt. Les traces fortement bruitées, ou très atténuées, pendant le temps d'intérêt sont écartées. A contrario les traces présentant une forte réflexion, ou réfraction, qu'on appelle 20 les traces propres, sont conservées. Une prise en compte de l'effet des multiples visibles sur les points de tirs et/ou obtenus lors du traitement des données sismiques existantes permet de sélectionner des traces pour lesquelles les forts multiples n'apparaissent pas pendant le temps d'intérêt de chaque trace considérée. La figure 4 montre cinq traces 41, 42, 43, 44 et 45. 25 L'axe 46 en ordonnée représente l'offset, c'est-à-dire la distance source-récepteur. La trace 41 est trop bruitée et les traces 42 et 45 sont trop atténuées pendant le temps d'intérêt. Ces trois traces ne sont donc pas sélectionnées. En revanche les traces 43 et 44 présentent un bon rapport signal sur bruit et une énergie visible pendant le temps d'étude de la zone 30 d'intérêt, elles sont donc considérées comme potentiellement optimales et sont ainsi conservées. Le deuxième paramètre concerne l'information dynamique qui peut être obtenue lorsque les données sismiques 4D sont disponibles. On peut, à partir des images 4D complètes, détecter si un 3059786 10 changement a déjà affecté la zone d'intérêt considérée. Si tel est le cas, l'analyse des traces brutes obtenues et provenant des acquisitions de base et ultérieures permet de déterminer les traces présentant les meilleurs effets 4D pour la zone d'intérêt. Ces traces sont celles qui ont le plus contribuées à 5 l'identification de ce changement et sont ainsi optimales pour y détecter de nouveau les effets 4D. En complément des traces sismiques 4D existantes, ou lorsqu'on ne dispose pas de celles-ci, l'utilisation de modèles pétro-élastiques pour la zone d'intérêt couplés au modèle dynamique de réservoir permet de modéliser une acquisition sismique synthétique estimant les effets sismiques 10 induits par la production. Le troisième paramètre concerne l'information relative aux facteurs de surface obtenue à partir du sous ensemble identifié par le premier paramètre ou la combinaison des deux premiers paramètres, on peut croiser la position des couples source-récepteur sélectionnés avec les informations de surface. Ce 15 croisement d'information permet de vérifier d'une part la faisabilité pratique du positionnement de ces couples source-récepteur, en vérifiant par exemple que les positions choisies ne se trouvent pas dans une zone d'exclusion. Il permet également de prendre en compte les forts générateurs de bruits sismiques pouvant apparaître entre la ou les données sismiques existantes, 20 comme par exemple une route en construction, une unité de compression de gaz ou un couloir maritime. Ces étapes permettent d'analyser la faisabilité de la détection sismique active ultralégère. Le procédé décrit jusqu'à maintenant peut être complété afin 25 d'améliorer le traitement sismique des systèmes actifs ultralégers et afin de détecter les changements dans la zone d'intérêt lorsque les couples source-récepteur ont été placés aux emplacements optimaux, de façon permanente ou semi-permanente. Cela permet de focaliser la détection sur la zone d'intérêt sans construire une image complète du sous-sol. Il s'agit d'une 30 approche trace à trace, où chacun des couples source-récepteur sont observés et traités indépendamment les uns des autres. Pour s'affranchir des changements intervenant au-dessus de la zone d'intérêt, on pose l'hypothèse que sur une échelle de temps courte, de l'ordre de 6 mois maximum, il est 3059786 11 possible d'identifier dans le sous-sol au moins une zone au-dessus de la zone d'intérêt qui n'est pas sujette à un changement. Cette zone est appelée la zone de référence. Elle sert à obtenir une mesure du bruit ambiant associé au système ultraléger de détection. Elle sert également à corriger les variations 5 arrivant au-dessus de cette zone de référence. Afin de corriger les effets 4D situés au-dessus de la zone de référence, où les variations proches de la surface sont les plus importantes, on utilise la zone de référence comme étalonnage en considérant que la réponse sismique de cette zone ne doit pas changer. Il est alors possible de corriger les variations sismiques de cette zone 10 et d'appliquer cette correction de façon à mettre en évidence les variations dans la zone d'intérêt. Il s'agit de variations dans le temps de l'amplitude, du temps de trajet ou des ondelettes. On peut voir sur la figure 5, la source 51, le récepteur 52, l'onde sismique 53, la zone d'intérêt 54 et la zone de référence 55 précédemment décrite. La figure 6 illustre un exemple d'application à une 15 trace brute du procédé objet de la présente invention. En abscisse, il est indiqué le temps en jour et en ordonnée la différence de temps de trajet de l'onde sismique au cours du temps. Le point de départ du trajet étant la source sismique, puis le réservoir et finalement le récepteur. La courbe supérieure correspond à un enregistrement journalier et la courbe inférieure à 20 la dérive par rapport au premier jour d'enregistrement. On constate une variation importante au 25ème jour ce qui signifie qu'il s'est produit un changement géomécanique. Ce changement géomécanique peut être, par exemple, une saturation, une modification de pression ou de température. Par ailleurs, lorsque l'on ne considère qu'une seule trace sismique 25 issue d'une acquisition ultralégère, la notion de cellule, ou bin en anglais, n'a plus de sens. La résolution spatiale de la zone détectée est considérée comme étant la zone de Fresnel, elle-même dépendante de la fréquence. Il est alors possible de filtrer par octave afin d'affiner la sensibilité de la détection. En abaissant la fréquence des données reçues, on augmente ainsi la sensibilité de 30 la détection dans les zones d'intérêt plus lointaines. Pour que cela soit possible, le système d'acquisition ultraléger doit avoir émis et enregistré plusieurs gammes de fréquence. Par exemple, pour une vitesse de 5500 mètres par seconde à la zone d'intérêt, les zones de Fresnel situées à 50, 100, 3059786 12 150, 200 et 275 mètres du centre du spot correspondent respectivement à des fréquences de 27, 14, 9, 7 et 5 Hertz, pour un modèle de vitesse homogène et un modèle d'acquisition zéro offset. Ainsi pour une détection de changement à 5 Hertz et qui n'apparait pas à 14 Hertz, cela permet de déduire que le 5 changement observé dans la zone d'intérêt intervient entre 100 et 275 mètres du centre de la zone d'intérêt. D'autre part, en utilisant les informations obtenues à l'étape de démigration, on peut estimer pour chaque couple source-récepteur la pente d'émission et de réception des ondes sismiques pour illuminer la zone 10 d'intérêt. En émission, il est possible de focaliser l'émission des ondes, ou beam steering en anglais, dans la pente et la direction optimales pour l'illumination de la zone d'intérêt en utilisant plusieurs sources ou une antenne de sources. Ainsi de nombreux facteurs de bruits tels que la diffraction, les échos parasites ou multiples et réflexions parasites sont évités et le seuil de 15 détectabilité est amélioré. En réception, il est possible, en utilisant des récepteurs trois composants, en moyennant une grappe de capteurs proches ou en utilisant une antenne de récepteurs, d'obtenir une trace moyenne ne contenant que les informations sismiques arrivant selon la direction et la pente définies dans l'étape de démigration. Ce filtrage permet de réduire fortement 20 les facteurs de bruit, notamment les multiples, et d'augmenter le seuil de détectabilité. Il s'agit d'augmenter la précision de détection en réduisant la zone de détection au sein de la zone de Fresnel et ainsi mieux focaliser la détection en jouant sur les pentes du beam steering. On peut alors, avec un seul couple source-récepteur, obtenir plusieurs focalisations différentes au sein 25 de cette zone de Fresnel ce qui revient à discriminer dans cette zone plusieurs spots différents. Ces deux filtrages peuvent être mis en oeuvre indépendamment ou ensemble afin d'augmenter encore plus le seuil de détectabilité. Si la démigration a été effectuée par gamme d'octave, il est alors possible d'affiner encore plus ce filtrage en modifiant les pentes d'émission et 30 de réception en fonction des octaves. L'invention permet ainsi de conclure à une détection prenant en compte d'une part les variations dans le temps dans la zone d'intérêt par rapport à une zone de référence et d'autre part les évolutions temporelles de ces variations. 3059786 13 L'invention peut également être mise en oeuvre pour enregistrer une trace complète pouvant être exploitée pour identifier ensuite d'autres spots dits d'opportunité situés au-dessus ou en dessous du spot initial et de jouer pour chacun sur la focalisation dans la zone de Fresnel. 5 REVENDICATIONS1. Procédé de détermination du positionnement optimal de couples source-récepteur aptes à acquérir des données sismiques comprenant les étapes suivantes : une première étape d'identification d'une zone d'intérêt (32) ayant fait l'objet d'une acquisition sismique antérieure afin d'obtenir une image de son sous-sol; une deuxième étape d'obtention des données sismiques acquises lors de l'acquisition sismique antérieure de ladite zone d'intérêt pendant un temps d'intérêt; une troisième étape d'application d'une démigration auxdites données sismiques afin de déterminer les positions de chaque couple source-récepteur (31, 34) ayant contribué à l'image dudit sous-sol de ladite zone d'intérêt pendant ledit temps d'intérêt; une quatrième étape d'obtention des traces brutes (41, 42, 43, 44, 45) pour lesdites positions des couples source-récepteur (31, 34); une cinquième étape de sélection d'au moins une trace brute optimale (43, 44) parmi lesdites traces brutes ; une sixième étape de détermination des positions optimales des couples source-récepteur (31, 34) correspondant à la dite trace brute optimale (43, 44); 2 Procédé selon la revendication 1 caractérisé en ce que la démigration de la troisième étape est effectuée par gamme d'octave. 3059786 15 3. Procédé selon l'une quelconque des revendications précédentes caractérisé en ce que l'acquisition sismique antérieure est une acquisition 4D. 4. Procédé selon la revendication 3 caractérisé en ce que la cinquième 5 étape comporte en outre une sélection des traces brutes optimales (43, 44) permettant de détecter les meilleurs effets 4D. 5. Procédé selon l'une quelconque des revendications précédentes caractérisé en ce que la cinquième étape comporte en outre l'utilisation de modèles pétro-élastiques afin de sélectionner les traces brutes 10 optimales (43, 44). 6. Procédé selon l'une quelconque des revendications précédentes caractérisé en ce que la cinquième étape comporte en outre une sélection tenant compte des facteurs de surface. 7. Procédé selon l'une quelconque des revendications précédentes 15 caractérisé en ce que la cinquième étape comporte en outre une sélection tenant compte d'une zone de référence (55) servant d'étalonnage de variations sismiques. 8. Procédé selon l'une quelconque des revendications précédentes caractérisé en ce que la troisième étape comporte en outre l'obtention 20 de la pente et de la direction d'émission et de réception optimales des ondes sismiques pour chaque couple source-récepteur (31, 34). 9. Procédé selon l'une quelconque des revendications précédentes caractérisé en ce que le procédé comprend en outre une septième étape de placement de couples source-récepteur (31, 34) auxdites positions 25 optimales et une huitième étape d'acquisition de nouvelles données sismiques auxdites positions optimales. 10. Procédé selon la revendication 9 caractérisé en ce que la huitième étape se base sur la pente et la direction d'émission et de réception optimales des ondes sismiques.