Brevet FR2845483A1 — publié le 20040409

PROCEDE DE REDUCTION DU BRUIT HARMONIQUE DANS DES SIGNAUX VIBROSISMIQUES DOMAINE GENERAL ET ETAT DE LA TECHNIQUE L'invention concerne un procédé de traitement pour la réduction du bruit harmonique dans un signal vibrosismique enregistré. On rappelle que pour reconstruire des images du sous-sol, les géologues ou géophysiciens utilisent classiquement des sources sismiques 10 destinées à générer des ondes qui se propagent dans le sous-sol et se réfléchissent au moins partiellement aux interfaces des différentes couches de celui-ci (réflecteurs). Les ondes réfléchies sont enregistrées en fonction du temps par des capteurs qui sont par exemple de type géophones en sismique terrestre. Les traces sismiques ainsi enregistrées sont traitées afin 15 d'en déduire une information sur la géologie du sous-sol. On utilise couramment en tant que sources, en particulier en sismique terrestre, des vibrateurs qui génèrent des signaux vibrosismiques. Les signaux vibrosismiques sont par exemple des signaux communément appelés " sweep " par l'homme du métier, dont la fréquence 20 lors d'une séquence d'émission varie continment de façon à balayer une gamme de fréquence donnée. A cet égard, il a été proposé par H.J Rozemond en 1996 une technique d'acquisition sismique autorisant la séparation de signaux vibrosismiques de type " sweep " émis par différentes sources et se recouvrant 25 dans le temps. On pourra notamment avantageusement se référer à " Slip sweep acquisition" - H.J. Rozemond - SEG 1996 International Conference - Denver. Les vibrateurs posent quant à eux des problèmes de distorsion par 30 rapport au signal théorique que l'on souhaiterait générer. Le signal généré par un vibrateur est en effet généralement constitué d'un balayage fondamental et de ses balayages harmoniques, lesquels sont généralement considérés comme des parasites qui peuvent affecter la qualité des données. En outre, le souci d'optimiser la productivité conduit à rechercher des écarts de temps entre vibrations aussi courts que possible. Dans ces 5 conditions, les enregistrements obtenus présentent un rapport signal à bruit dégradé par rapport à celui qui serait obtenu en utilisant les sources sans recouvrement dans le temps. Les données brutes reçues au niveau d'un capteur sont en particulier bruitées par de l'énergie harmonique due au chevauchement des signaux émis par les différentes sources Il a déjà été proposé des méthodes destinées à limiter le bruit de corrélation en vibrosismique. Notamment, le document GB 2 348 003 décrit une méthode qui s'applique à des ensembles de sismogrammes préalablement traités et regroupés en points miroirs (de façon que les réflexions d'un même point de 15 l'espace se trouvent au même point en temps et en profondeur). Ces sismogrammes sont ensuite décomposés en bandes de fréquences sur lesquelles on opère une discrimination statistique du signal et du bruit. Le traitement décrit dans ce document nécessite toutefois de disposer de sismogrammes enregistrés en parallèle sur plusieurs capteurs. 20 D'autres méthodes encore nécessitent de disposer de plusieurs acquisitions successives dans le temps pour un même capteur. C'est le cas en particulier de celle décrite dans le document US 5.410.517. On connaît par ailleurs par l'article: " Amplitude analysis of harmonics on vibrator generated direct 25 waves " - EAGE 64th Conference & Exhibition - Florence - 27-30 May 2002 une méthode permettant de déterminer les rapports d'amplitude entre les différents harmoniques du signal source lorsque l'on dispose des composantes harmoniques du signal enregistré. A cet effet, il est proposé de mettre en ceuvre des transformées en ondelettes sur ces différentes 30 composantes. L'article montre que les rapports des transformées de Fourier de ces transformées en ondelettes sont indépendants notamment de la réponse du sous sol, de sorte que ces rapports correspondent aux rapports des amplitudes des composantes harmoniques de la source. On notera toutefois qu'il n'est aucunement proposé dans cet article une méthode destinée à permettre de réduire le bruit harmonique d'un signal vibrosismique. PRESENTATION DE L'INVENTION Un but de l'invention est de proposer un procédé de réduction du bruit harmonique dans des données vibrosismiques enregistrées. Le procédé proposé par l'invention présente en particulier l'avantage de pouvoir être appliqué individuellement à tout enregistrement 10 et de ne pas nécessiter plusieurs acquisitions successives ou plusieurs acquisitions en parallèle sur différents capteurs. A cet effet, l'invention propose un procédé de réduction du bruit harmonique dans un signal vibrosismique enregistré sur un capteur, ce signal vibrosismique correspondant à une séquence vibratoire donnée 15 émise par au moins une source et propagée dans un sous-sol dans lequel elle est réfléchie, caractérisé par les étapes selon lesquelles: a) on corrèle le signal vibrosismique ou une partie de celui-ci avec un signal correspondant la composante fondamentale de la séquence de 20 vibrations émise par la source b) on sélectionne sur le signal corrélé une partie qui correspond à l'énergie de la composante fondamentale du signal vibrosismique et on applique à la partie sélectionnée un opérateur d'estimation de la composante harmonique dudit signal corrélé c) on en déduit une estimation du bruit harmonique du signal vibrosismique et l'on soustrait cette estimation dudit signal. Comme on l'aura compris, un tel procédé tire partie - d'une part de ce que la corrélation d'un signal vibrosismique par la composante fondamentale de la séquence théoriquement émise 30 par la source présente la propriété de séparer l'énergie correspondant à la partie fondamentale du signal enregistré de celle correspondant à la partie harmonique; l'énergie fondamentale est en effet répartie autour du temps d'origine (temps d'arrivée sur le capteur), tandis que l'énergie correspondant à la partie harmonique du signal est répartie vers les temps négatifs par rapport à la partie fondamentale du signal, - et d'autre part de ce qu'il est possible d'estimer le rapport entre l'énergie fondamentale et l'énergie harmonique indépendamment de la connaissance du sous-sol, différentes méthodes étant à cet égard envisageables. PRESENTATION DES FIGURES D'autres caractéristiques et avantages ressortiront encore de la description qui suit, laquelle est purement illustrative et non limitative et doit être lue en regard des figures annexées parmi lesquelles: - la figure 1 représente schématiquement l'architecture générale d'un 15 dispositif d'acquisition sismique; - les figures 2a à 2c illustrent respectivement une séquence d'émission de signaux " slip sweep " (figure 2a), un signal d'enregistrement issu de la propagation et de la réflexion de ces signaux dans le sous-sol (figure 2b), ainsi qu'un signal qui 20 correspond à la corrélation de ce signal d'enregistrement par la séquence fondamentale de vibrations correspondant aux signaux de la figure 2a; - la figure 3 est un organigramme qui illustre différentes étapes d'un mode de mise en òuvre possible de l'invention; - la figure 4 illustre la décomposition du signal enregistré en fonction du signal source et de la réfléctivité du sous-sol; - les figures 5a à 5c illustrent les opérations de corrélation et de fenêtrage mises en oeuvre dans le cadre d'un prétraitement de détermination d'un opérateur d'estimation de bruit harmonique. 30 DESCRIPTION D'UN OU PLUSIEURS MODES DE REALISATION Généralités - Rappels On a représenté sur la figure 1 un exemple de dispositif d'acquisition de données vibrosismiques, qui comprend un ensemble de sources vibrosismiques Si, S2,... Sm positionnées à la surface du sol, ainsi qu'un ensemble de géophones G1, G2,...Gp également placés à la surface 5 du sol. Les vibrations émises par les sources Si, S2,...Sm (signaux s1, s2,... m) se propagent dans le sous-sol et sont réfléchies aux interfaces des couches de celui-ci. Les géophones G1, G2,.. .Gp convertissent les signaux qu'ils reçoivent en signaux électriques qui sont transmis à une unité de traitement 1 o ils sont numérisés et enregistrés. Les sismogrammes ainsi 10 enregistrés sont généralement traités non pas sur le terrain, mais dans un centre de traitement - référencé en l'occurrence par 2 - prévu à cet effet. On se place dans la suite de la présente description dans le cas o les signaux émis par les sources sont de type "sweep " et se chevauchent (émission de type " slip sweep " selon la terminologie habituellement 15 utilisée par l'homme du métier). On notera toutefois que l'invention ne se veut aucunement limitée à des signaux sources de type " slip sweep " et s'applique de façon générale au traitement de tout signal vibrosismique issu d'un ou plusieurs signaux sources pour lesquels on connaît le ou les signaux sources théoriquement 20 émis. Des signaux sl(t),...s4(t) de type " slip sweep " ont été représentés sur la figure 2a pour les sources SI à S4. Comme on peut le voir sur cette figure, chaque source émet la même séquence vibratoire, les séquences d'émission étant décalées dans 25 le temps d'une source à l'autre. Lors d'une séquence d'émission, la fréquence du signal vibratoire varie de façon continue, par exemple selon une loi linéaire ou exponentielle en fonction du temps. Dans l'exemple de la figure 2a, la fréquence croît dans le temps; des lois de variations avec lesquelles la fréquence décroît 30 dans le temps sont également possibles. On a porté sur la figure 2a la durée T d'une séquence d'émission et le temps ts qui sépare les séquences d'émission de deux sources émettant successivement dans le temps. Le temps ts est inférieur à la durée d'émission T, de sorte que les vibrations générées par ces deux sources se chevauchent. Les signaux reçus au niveau des géophones G1, G2,...Gp sont des signaux composites - différents d'un géophone à un autre - qui sont la 5 somme de différents signaux issus de la réflexion des signaux émis par les différentes sources Si, S2,..Sm.sur les différents réflecteurs du sous-sol. Un tel signal est par exemple le signal e(t) représenté sur la figure 2b. La corrélation de ce signal avec un signal de référence 10 correspondant à la séquence vibratoire d'émission théoriquement émise est un signal qui est théoriquement, en l'absence de bruit harmonique, un signal du type de celui illustré sur la figure 2c et qui présente des pics aux temps d'arrivée, sur le géophone considéré, des séquences émises par les différentes sources après réflexion dans le sous-sol. 15 Exemple de mise en oeuvre On a représenté sur la figure 3 différentes étapes d'un traitement conforme à un exemple de mise en oeuvre possible de l'invention. Selon ce traitement, on met en oeuvre sur un signal bruité reçu sur un géophone (étape 10 de récupération des données d'entrée) une corrélation avec un signal de référence Hi qui correspond à une séquence théorique d'émission de vibrations (étape 20). Cette étape de corrélation permet de mettre en évidence sur le 25 signal reçu les temps qui correspondent aux arrivées des réflexions par les réflecteurs du sous-sol La corrélation mise en coeuvre à cette étape 20 présente en effet la propriété de répartir l'énergie du signal enregistré de la manière suivante: - l'énergie correspondant à la partie fondamentale du signal 30 enregistré est répartie autour du temps d'origine (temps d'arrivée), - l'énergie correspondant à la partie harmonique du signal, dans le cas o les fréquences de celui-ci augmentent avec le temps, est répartie vers les temps négatifs par rapport à la partie fondamentale du signal. Dans l'exposé qui suit, on considérera un signal dont les fréquences augmentent avec le temps, ce qui ne nuit pas à la généralité de la méthode, étant entendu que si les fréquences suivaient une loi descendante, l'énergie serait répartie vers les temps positifs. Cette corrélation permet par conséquent de mettre en évidence les temps d'arrivée. Pour chacune des sources Sm à S2, on met alors le traitement de débruitage des étapes 30 à 60 qui vont maintenant être décrites. On notera ici que ce traitement n'est pas nécessaire pour la source qui correspond à la 10 première arrivée sur le géophone, dans la mesure o le bruit généré par cette source est rejeté en dehors de la plage de temps conservée pour la suite du traitement.Les étapes 30 à 60 sont les suivantes. Pour une source k donnée (avec k entier compris entre 2 et m), on met en oeuvre sur le signal enregistré reçu sur le géophone un fenêtrage 15 incluant le temps d'arrivée mis en évidence pour cette source par l'étape 20 de corrélation (étape 30), de manière à conserver au maximum l'énergie de la composante fondamentale. Puis on convolue le signal corrélé ainsi filtré avec un opérateur OPk d'estimation de la composante harmonique qui est propre à la source k. Cet 20 opérateur opki aura été préalablement déterminé à l'occasion d'un prétraitement 70 qui est décrit plus loin plus en détails. Cette opération de convolution (étape 40) permet de fournir le bruit Nk qui sur le signal corrélé correspond aux composantes harmoniques de la source k. Une fois ce bruit Nk déterminé, il est soustrait du signal corrélé 25 (étape 50). On dispose ainsi d'un signal débruité par rapport aux harmoniques générés par la source k (étape 60). Un traitement de déconvolution permet de déterminer le signal vibrosismique dépourvu du bruit harmonique correspondant à ladite source. Les étapes 30 à 50 sont ainsi mises en oeuvre successivement pour chacune des sources Sm à S2, le signal corrélé utilisé pour mettre en oeuvre ces étapes pour la source Sk étant, lorsque k est différent de m, le signal débruité obtenu à l'étape 60 du traitement mis en òuvre pour la source Sk +1. Lorsque le traitement des étapes 30 à 60 a été mis en oeuvre pour chacune des sources Sm à S2, on dispose alors d'un signal enregistré sur lequel le bruit harmonique a été fortement réduit. Exemple de détermination de l'opérateur d'estimation du bruit harmonique généré par une source donnée On va maintenant décrire un exemple possible pour le pré-traitement de détermination des opérateurs OPk (étape 70) 1) Eléments théoriques Un signal sk émis par une source Sk donnée peut être modélisé, 15 ainsi que l'illustre la partie gauche de la figure 4, comme une somme pondérée de n signaux hi...hn de référence correspondant à la séquence vibratoire fondamental (hi) et à ses harmoniques (h2,..., hn), ces différents signaux de référence étant pondérés par des fonctions de poids c1...c, On peut donc écrire: n Sk= ici *hi [1] I avec: hl =a cos(2zg) hi= a cos(2zLig) o * est l'opération de convolution, o a est une valeur d'amplitude et o g est une fonction connue qui varie en fonction du temps et qui peut 25 s'écrire: di (o f(t) est la fonction de phase). Par ailleurs, on sait qu'en vibrosismique conventionnelle, le signal ek en provenance d'une source unique Sk et enregistré par un géophone donné est considéré comme pouvant être décrit de la manière suivante: es = (Sk * r + n) * imp [2] o sk est le signal émis par une source Sk, r est une fonction représentative des coefficients de réflexion du signal sk sur les couches souterraines, n est une fonction représentant un bruit aléatoire, et imp est une fonction représentant la réponse impulsionnelle de la chaîne d'acquisition utilisée. Si l'on fait l'hypothèse, qui est illustrée sur la figure 4, que le bruit non harmonique est négligeable et que la réponse impulsionnelle de la chaîne d'acquisition n'intervient pas, on peut alors considérer que l'enregistrement est la convolution du signal source sk et des coefficients de réflexion r, c'est-à-dire: ek =Sk *r [3] 2) Différentes étapes du traitement 70 Conformément à la variante Ilustrée sur la figure 3, le prétraitement de détermination des opérateurs OPk des différentes sources Sk peut être 20 mis en oeuvre de la façon suivante. Ce prétraitement nécessite d'utiliser au moins un enregistrement obtenu en faisant émettre la séquence de vibrations sk par la source Sk. Lorsque l'on dispose du sismogramme ainsi enregistré pour une source donnée (étape 701), on corrèle ce signal avec la séquence de 25 vibrations fondamentale hi et ses harmoniques h2,..., hn. On rappelle que la corrélation avec l'un des harmoniques ou fondamental hi correspond à une convolution avec une fonction hdi dont le spectre d'amplitude est égal au spectre d'amplitude de hi et le spectre de phase est égal à l'opposé du spectre de phase de hi et qu'une telle 30 corrélation présente les propriétés suivantes: - l'autocorrélation de hi est de phase nulle et est centrée autour du temps zéro; - la corrélation de hi et de hj rejette l'énergie aux temps positifs lorsque j est supérieur à i et au temps négatif lorsque i est supérieur à j. C'est ce qu'illustrent les figures 5a à 5c qui représentent 5 respectivement la corrélation avec hi, h2 et h3 d'un signal combinaison pondérée de hi, h2 et h3. Comme on le comprend sur ces figures, l'énergie mise en évidence par ces opérations est principalement centrée autour du temps zéro (temps d'arrivée de la séquence de vibrations), les termes correspondant à des 10 corrélations de hi et hj croisées (c.a.d. avec i et j différents) étant quant à eux d'énergie beaucoup plus faible. Il est à noter que seule la condition relative au spectre de phase est nécessaire pour l'obtention du résultat visé. En mettant en oeuvre (étape 703) un fenêtrage centré autour de la 15 zone o le sismogramme est le plus énergétique, c'est-à-dire autour de ce temps zéro (sans être d'ailleurs nécessairement centré sur celui-ci, l'essentiel étant que le pic d'énergie se trouve dans la fenêtre de filtrage), on récupère alors une approximation de la fonction r*ci*hi*hdi. L'opération de fenêtrage est par exemple réalisée en utilisant une 20 fonction d'apodisation ap choisie pour présenter une largeur couvrant parfaitement l'énergie de l'harmonique considéré. On a alors, en posant f =ap (ek *hd): n f, ap (e, * hdi) = ap (r* Ici * hi) * hd [4 fi r * c * h, * hd1 Cette égalité s'écrit dans le domaine de Fourrier: Fi = R -Ci * HiHDI [51 o Fi, R, Ci, Hi et HDi correspondent aux transformées de Fourier de fi, r, ci, hi, hdi. Les opérations de corrélation, puis de fenêtrage permettent par conséquent de déterminer des fonctions CR1 dites de poids réduites: Ci F i} Hl -HD1 CRó - C / F, //l D 6] Cl FI Hi-H HDi qui sont telles que si l'on note cri les transformées de Fourrier inverses desdites fonctions CRi., le signal source vérifie: Sk = cl (hl + cr2 * h2 + cr3 * h3 +...crn * h.) [7] Après détermination de ces coefficients (étape 704), on dispose donc, à une fonction filtre c1 près, d'une modélisation du signal émis par la source réelle Sk. On notera que dans une mise en oeuvre possible, pour déterminer les fonctions CRi, on se limite uniquement aux fréquences pour lesquelles 10 les valeurs de C1 et Ci ne sont pas nulles. Egalement, pour obtenir une approximation plus fine des fonctions CRy, l'étape 704 peut être appliquée à plusieurs traces correspondant par exemple à plusieurs récepteurs d'une même acquisition. Les valeurs des CR1 peuvent être alors être calculés en faisant des calculs de moyenne ou 15 encore à l'aide d'opérateurs statistiques. Une fois les coefficients de poids harmoniques calculés pour une source Sk considérée, on détermine (étape 705) au moyen de ces coefficients un opérateur d'estimation du bruit harmonique généré par la source Sk. Pour une source Sk donnée, le bruit harmonique correspond dans le domaine fréquentiel à la multiplication - de la composante fondamentale ECfod de la transformée de Fourier de l'enregistrement corrélé ECfond = R.C1.H1. HD1 - par l'opérateur OP égal à n n E CRi.Hi X Ci.Hi OP = i=2 HI CI.= 2845483 12 Si l'on s'en tient en effet à l'équation [1], ainsi qu'au modèle de l'équation [3], on a: n EC = R.S.HD1 =R.(XC,.H,).HD1 i=1 et ECfond = R.C1.H1.HD1 n ECharm = R.( Ci.H,).HD1 = OP.ECfond i=2 Dans le domaine temporel (opérateur op), cette opération s'écrit 10 comme la convolution de la transformée de Fourrier inverse de OP par la partie positive de l'enregistrement corrélé (laquelle est assimilable à la partie fondamentale du signal enregistré). Autres variantes de calcul des fonctions de poids harmoniques De nombreuses autres variantes de calcul des fonctions de poids harmoniques sont envisageables. Notamment, si l'on considère que la corrélation d'un signal 5 vibrosismique par la composante fondamentale théorique hl ou un harmonique hi est telle que l'énergie correspondant à la corrélation de cette composante ou harmonique avec des composantes harmoniques de rangs supérieur est négligeable, on peut décrire un signal vibrosismique e d'une part par la composante fondamentale théorique et d'autre part par les 10 différentes composantes harmoniques que l'on souhaite considérer: cl*r*hl*hdl = e*hdl cl*r*hl*hd2 + c2*r*h2*hd2 = e*hd2 cl*r*hl*hd3 + c2*r*h2*hd3 + c3*r*h3*hd3 = e*hd3 15..... cl*r*hl*hdn + c2*r*h2*hdn + c3*r*h3*hdn +.....cn*r*hn*hdn = e*hdn On dispose par conséquent d'un système de n équations permettant de déterminer les fonctions R.Ci/R.C1, c'est-à-dire les fonctions de poids. Il est en outre à noter que si l'on dispose de données relatives à la signature de la source, telles que le signal force (somme pondérée des accélérations mesurées sur la masse et la plaque du vibrateur), le signal source Sk disponible est avantageusement utilisé dans le pré-traitement de l'étape 70 en tant que signal ek.. Dans ce cas il n'y a pas de réflexion et 25 l'équation (3) se réduit à: ek = Sk Comme on l'aura compris, le procédé de réduction de bruit proposé peut être mis en oeuvre soit après l'enregistrement des données bruitées, soit en temps réel, les données enregistrées étant alors des données 30 débruitées. REVENDICATIONS 1. Procédé de réduction du bruit harmonique dans un signal 5 vibrosismique enregistré sur un capteur, ce signal vibrosismique correspondant à une séquence vibratoire donnée émise par au moins une source et propagée dans un sous-sol dans lequel elle est réfléchie, caractérisé par les étapes selon lesquelles: a) on corrèle le signal vibrosismique ou une partie de celui-ci avec 10 un signal correspondant la composante fondamentale de la séquence de vibrations émise par la source b) on sélectionne sur le signal corrélé une partie qui correspond à l'énergie de la composante fondamentale du signal vibrosismique et on applique à la partie sélectionnée un opérateur d'estimation de la 15 composante harmonique dudit signal corrélé c) on en déduit une estimation du bruit harmonique du signal vibrosismique et l'on soustrait cette estimation dudit signal. 2. Procédé selon la revendication 1, caractérisé en ce que la partie du signal corrélé que l'on sélectionne correspond: - à l'énergie autour du temps d'arrivée sur le capteur de la séquence vibratoire propagée et réfléchie - et/ou à l'énergie dans les temps positifs par rapport audit temps 25 d'arrivée. 3. Procédé selon l'une des revendications 1 ou 2, caractérisé en ce que dans le cas o la séquence vibratoire est émise par plusieurs sources séparées dans l'espace et/ou dans le temps, on corrèle le signal 30 vibrosismique avec la composante fondamentale théorique de vibrations de la séquence vibratoire émise et on déduit de cette corrélation les temps d'arrivée sur le capteur des émissions des différentes sources après propagation et réflexion dans le sous-sol. 4. Procédé selon la revendication 3, caractérisé en ce qu'on isole dans le signal vibrosismique une partie du signal autour du dernier temps d'arrivée mis en évidence par la corrélation avec la composante fondamentale théorique, on met en oeuvre à partir de la partie ainsi isolée 5 les étapes a) à c), et on réitère ce traitement en remontant successivement les différents temps d'arrivée, le signal vibrosismique sur lequel on isole une partie de signal autour du temps d'arrivée considéré étant alors à chaque fois le signal débruité issu du traitement correspondant au temps d'arrivée juste précèdent. 5. Procédé selon l'une des revendications précédentes, caractérisé en ce que l'opération correspondant à l'étape b) est mise en oeuvre dans le domaine fréquentiel et correspond à la multiplication dans le domaine fréquentiel de la partie du signal corrélé qui correspond à l'énergie de la 15 composante fondamentale du signal vibrosismique, par un opérateur OP tel que: n CR '.H, oP =i=2 Hi o: HI est la transformée de Fourrier de la composante fondamentale théorique de vibrations de la séquence émise par la source, Hi est la transformée de Fourrier de la composante harmonique d'ordre i, et o CRi désigne une fonction fréquentielle de poids associée à 25 l'harmonique d'ordre i. 6. Procédé selon l'une des revendications 1 à 4, caractérisé en ce que l'opération correspondant à l'étape b) est mise en oeuvre dans le domaine temporel et correspond à la convolution dans le domaine temporel 30 de la partie du signal corrélé qui correspond à l'énergie de la composante fondamentale du signal vibrosismique, par un opérateur qui correspond à la transformée de Fourrier inverse d'un opérateur OP de multiplication dans le domaine fréquentiel tel que: n E CRI.Hi OP HI o: Hi est la transformée de Fourrier de la composante fondamentale théorique de vibrations de la séquence émise par la source, Hi est la transformée de Fourrier de la composante harmonique d'ordre i, et o CRi désigne une fonction fréquentielle de poids associée à l'harmonique d'ordre i. 7. Procédé selon l'une des revendications 5 ou 6, caractérisé en ce que l'opérateur utilisé lors de l'étape b) est déterminé par un prétraitement 15 mis en oeuvre sur au moins un enregistrement obtenu à partir de la séquence émise par la source. 8 Procédé selon l'une des revendications 5 ou 6, caractérisé en ce que l'opérateur utilisé lors de l'étape b) est déterminé par un prétraitement 20 mis en oeuvre sur une estimation du signal émis obtenue à partir de données acquises directement sur la source.. 9. Procédé selon l'une des revendications 7 ou 8, caractérisé en ce que l'opérateur ainsi déterminé dans le domaine temporel ou dans le 25 domaine fréquentiel est tel que: CRi= Fi x H. HD, F. Hi HDi o HDI est la transformée de Fourrier d'une fonction hdM dont le spectre de phase est l'opposé de celui de la composante fondamentale théorique, HDi est la transformée de Fourrier d'une fonction hdi dont le spectre de phase est l'opposé de celui de l'harmonique i, et o Fi est la transformée de Fourrier de f= ap -(e*hdi), * étant l'opération de convolution, " e " étant un signal vibrosismique sur lequel le 5 prétraitement est mis en oeuvre, " ap " étant une fonction de fenétrage centrée sur centrée autour d('une zone énergétique du sismogramme.. 10. Procédé selon la revendication 9, caractérisé en ce que la fonction " ap " est une fonction d'apodisation. 10 11. Procédé selon la revendication 6, caractérisé en ce que pour déterminer l'opérateur utilisé lors de l'étape b): - on corrèle un signal vibrosismique avec la composante fondamentale théorique hl de la séquence vibratoire et avec celle(s) de ses composantes 15 harmoniques hi que l'on souhaite prendre en compte, on multiplie les données corrélées avec une fenêtre excluant l'énergie associée à la corrélation des harmoniques d'ordre supérieur à celui de la composante avec laquelle le signal vibrosismique a été corrélé, - on résout le système d'équations linéaires ainsi obtenus afin d'en déduire 20 les rapports ci/cl o ci et cl sont des fonctions de pondération tels que n s = c1* hi 1=l o s est la séquence vibratoire émise par la source.