FR 2687227 A1 19930813 FR 9201343 A 19920206 I La présente invention concerne un procédé pour déterminer une signature FF(t) en champ lointain d'une pluralité de N sources sismiques Si (izl,, N). L'invention trouve une application particulièrement avanta- geuse dans le domaine de la prospection géophysique par sismique marine. L'exploration géologique des fonds marins utilise générale- ment des navires remorquant des dispositifs d'émission acoustique qui émettent vers les fonds à explorer des ondes acoustiques dans différentes directions Ces ondes acoustiques sont ensuite réfléchies par le sous-sol marin et détectées par des transducteurs disposés le long de câbles, appelés flûtes sismiques, tirés également par le navire de prospection L'analyse des ondes réfléchies par rapport aux ondes émises par les sources renseigne sur la structure, la nature et la composition des fonds étudiés. Les dispositifs d'émission acoustique mis en oeuvre en sismique marine constituent le plus souvent des sous-ensembles de sources immergés, ou canons, formant un réseau régulier de une à trois dimensions Pour émettre une onde acoustique dans une direction donnée, on commande l'émission des sources selon une loi de phase qui régit les retards de phase à appliquer aux sources En particulier, pour émettre dans ces conditions une onde dans la direction verticale, perpendiculaire à la surface de l'eau, la loi de phase nominale est uniforme en ce sens que l'émission des sources se fait en même temps. De façon à pouvoir contrôler que l'émission de l'onde acoustique incidente s'est effectuée sans anomalies, et également pour apporter d'éventuelles corrections aux signaux enregistrés par les hydrophones des flûtes sismiques il est nécessaire, à chaque tir des canons, de disposer d'un signal de référence, représentatif de l'onde acoustique telle qu'elle est reçue à grande distance par les structures géologiques à étudier Ce signal de référence est appelé signature en champ lointain et noté FF(t). De façon pratique et lorsque cela est possible, la signature en champ lointain du système de sources sismiques est mesurée directement en immergeant à une profondeur pouvant atteindre 1000 m un hydrophone situé à l'aplomb du réseau des N sources sismiques. Toutefois, lorsque la prospection géologique concerne des plateaux continentaux situés à faible profondeur, 200 à 300 m par exemple, il n'est plus possible d'acquérir la signature en champ lointain par simple mesure directe On a alors recours à des méthodes indirectes, telles que celle décrite dans la demande de brevet européen N O O 400 769, qui consistent généralement à reconstituer une signature en champ lointain à partir de la mesure d'une signature en champ proche N Fi(t) de chaque source sismique Si Dans ce but, un hydrophone est associé à chaque canon, et le signal détecté par ledit hydrophone constitue la signature en champ proche de l'émetteur correspondant. Pour passer des signatures en champ proche N Fi(t) à la signature en champ lointain FF(t), les procédés connus appliquent des algorithmes de correction à chaque signature en champ proche puis effectuent une sommation des signatures en champ proche ainsi corrigées. Toutefois, ce type de procédé présente l'inconvénient de mettre en oeuvre des algorithmes théoriques reposant sur un certain nombre d'hypothèses qui, souvent mal maîtrisées, jettent un doute sur la précision du résultat obtenu De plus, le traitement des données brutes à savoir les signatures en champ proche, exige des temps de calcul relativement longs et l'emploi d'ordinateurs puissants. Aussi, le problème technique à résoudre par l'objet de la présente invention est de proposer un procédé pour déterminer une signature FF(t) en champ lointain d'une pluralité de N sources sismiques Si (i=I,, N), comprenant une étape de mesure d'une signature N Fi(t) en champ proche de chaque source Si, procédé qui serait de mise en ceuixre beaucoup plus simple que les procédés connus et qui conduirait à dtes résultats plus sûrs. La solution au problème technique posé est, selon l'in\ention, remarquable en ce que ledit procédé comprend les opérations consistant respectivement: dans une étape préalable, à mesurer simultanément une signature \Fi ot) en champ proche initiale de chaque source sismique Si et une signa:,re F Fo(t) en champ lointain initiale de la pluralité des N sources sisn-i 2 te. à déterminer un opérateur de passage G défini par le produ Me convolution N Z N Fio(t) * G F Fo(t) i=l à chaque tir des N sources sismiques Si, à réaliser ladite étape de mesure des signatures N Fi(t) en champ proche, à déterminer ladite signature FF(t) en champ lointain par la relation N FF(t) = J EN Fi(t) * G i=l Ainsi, le procédé selon l'invention repose sur la seule hypothèse, bien vérifiée au demeurant, de la linéarité de la relation entre la signature en champ lointain et les signatures en champ proche Il n'est pas nécessaire en particulier de faire appel à des modèles théoriques dont la validité n'est pas toujours assurée D'autre part, l'opérateur de passage G étant déterminé une fois pour toutes à la suite de l'étape préalable, le calcul de la signature en champ lointain cherchée en fonction des signatures en champ proche peut être effectué dans un temps beaucoup plus court qu'avec les procédures connues. La description qui va suivre en regard du dessin annexé, donné à titre d'exemple non limitatif, fera bien comprendre en quoi consiste l'invention et comment elle peut être réalisée. La figure I est une vue schématique en coupe d'un ensemble de prospection par sismique marine utilisant un procédé conforme à l'invention. L'ensemble de prospection montré sur la figure 1 comprend un navire 100 remorquant au moins une flûte sismique 200 de grande longueur, 3 kilomètres par exemple, le long de laquelle sont disposés des capteurs acoustiques, non représentés, destinés à enregistrer l'onde réfléchie par le fond marin 300 étudié, en réponse à une onde acoustique incidente émise par un réseau 400 comportant une pluralité de N sources sismiques Si a'ec i= 1, N. De manière à déterminer une signature en champ lointain FF(t) de ces sources, un hydrophone Hi est associé à chacune desdites sources afin de mesurer, à chaque tir, une signature en champ proche de la source considérée. Dans une étape préalable, réalisée au cours d'une campagne de prospection o la mesure directe du champ lointain est possible, une signature N Fio(t) en champ proche initiale est enregistrée pour chaque source, tandis qu'un hydrophone Ho, immergé à une grande profondeur, fournit une signature F Fo(t) en champ lointain initiale. On détermine ensuite un opérateur de passage G reliant la N somme N N Fio(t) des signatures en champ proche initiales à la signature i=l F Fo(t) par le produit de convolution N i N Fio(t) * G = F Fo(t) i=l L'opérateur de passage G peut être déterminée de la façon suivante. N La somme N Fo(t) = N Fio(t) est échantillonnée sur i=l m+l 1 valeurs: N Fo(t) = (nfo, nf 1,, nfm) De même, la signature F Fo(t) en champ lointain initiale est échantillonnée sur n+l valeurs: F Fo(t) = (ffo 0, ffl '" ffn) L'opérateur G cherché est lui-même défini par p+ 1 valeurs G = (go' g 1 ' p gp) La méthode de calcul des coefficients de G consiste à rendre le produit N Fo(t) voisin de F Fo(t) au sens des moindres carrés. Le problème revient alors à trouver G tel que N l(nf*g) -ff l 2 j=O soit minimum. En developpant et en dérivant par rapport aux inconnus gk on obtient les p équations suivantes f (A O g O + Alg 1 + + Apgp = CO O )? A Ig Ok + A O g+ + Ap_ 1 gp = C ) (r) Apgo Apg 1 + + Aogp = Cp m avec Ak = Znf nfqk 1 5 -q-Oq q+ m Ck = nfq ffq+k q=o (r) est un système de Toeplitz, se résolvant par l'algorithme de Wiener-Lewison avec second membre. De façon à limiter l'amplification du bruit, on peut introduire un coefficient ú 2 de pondération Dans ce cas, le minimum recherché est celui de: n n 7-lnfg) ff l 2 +2 g 2 j 2 j O J Le système d'équation auquel on aboutit est analogue à (r) si ce n'est que le coefficient A O est remplacé par A O + È 2. L'opérateur de passage G étant ainsi déterminé pour une configuration d'émetteurs donnée, le procédé conforme à l'invention consiste alors à réaliser, lors de la campagne de prospection considérée, pour chaque tir des N sources sismiques Si, la mesure des signatures N Fi(t) en champ proche de chaque source, et à déterminer la signature FF(t) en champ lointain par la relation: N FF(t) =IN Fi(t) * G i::l On dispose de la sorte d'un moyen de contrôle du système de sources En effet, en cas d'anomalie, telle qu'un mauvais fonctionnement d'une source ou une déformation de la géométrie du réseau, la signature en champ lointain aura une allure différente de celle préalablement définie, ce qui provoquerait la non-validation du tir correspondant. La signature en champ lointain est également utilisée pour ramener les signaux détectés par les capteurs situés dans les flûtes à des signaux de phase minimale. REVENDICATION 1 Procédé pour déterminer une signature FF(t) en champ lointain d'une pluralité de N sources sismiques Si (i=l,, N), comprenant une étape de mesure d'une signature N Fi(t) en champ proche de chaque source sismique Si, caractérisé en ce que ledit procédé comprend les opérations consistant respectivement: dans une étape préalable, à mesurer simultanément une signature N Fio(t) en champ proche initiale de chaque source sismique Si et une signature F Fo(t) en champ lointain initiale de la pluralité des N sources sismiques, à déterminer un opérateur de passage G défini par le produit de convolution N >E 7 N Fio(t) * G = F Fo(t) i= 1 à chaque tir des N sources sismiques Si, à réaliser ladite étape de mesure des signatures N Fi(t) en champ proche, à déterminer ladite signature FF(t) en champ lointain par la relation N FF(t) = N Fi(t) * G i= 1