Brevet FR2805897A1 — publié le 20010907

La présente invention se rapporte, d'une manière générale, à l'analyse de signaux et au traitement de données sismiques pour des relevés géophysiques et, particulièrement, à un procédé de séparation et de synthèse signaux sismiques en composantes harmoniques autres. L'industrie de l'exploration sismique utilise des impulsions acoustiques pour induire des vibrations sonores dans la terre pour dessiner la structure souterraine pour exploration et le développement miniers. Ces impulsions acoustiques peuvent provenir d'une source explosive, implosive ou à balayage (ou fluctuation) de frequence. on utilise, pour produire un enregistrement de -champ sismique, enregistrement des fronts d'onde de réflexion et de refraction acoustiques qui se propagent depuis source jusqu'à un récepteur. Les variations des temps parcours d événements de réflexion et de réfraction dans ces enregistrements de champ indiquent la position surfaces réflexion à l'intérieur de la terre. L'analyse et la corrélation des événements dans un ou plusieurs enregistrements de champ dans un traitement données sismique produisent une image acoustique qui montre la structure souterraine. Les images acoustiques s utilisent pour trouver des gisements miniers de valeur. La source sismique du type à balayage ou fluctuation fréquence peut utiliser un signal pilote relativement long, comme de 2 à 15 secondes, pour garantir que l'on induise une énergie suffisante dans la terre. procédé à source du type à balayage ou fluctuation de fréquence se à la compression de signal pour comprimer signal et garantir une définition verticale suffisante pour trouver position de réflecteurs souterrains. La compression de signal est appelée, généralement, déconvolution, dans de nombreuses techniques bien connues dans l'art traitement données sismiques. La déconvolution de signaux à balayage ou à fluctuation de fréquence comprime le signal source en un signal beaucoup plus court représentatif d'une frontière réfléchissante souterraine. La précision et l'efficacité de toute technique de déconvolution est directement liée au degré de connaissance ou de compréhension l'on a du signal de source. La plupart des opérateurs de déconvolution proviennent d'estimations statistiques forme d'onde réelle de source. En pratique la forme d'onde de source à balayage ou fluctuation fréquence qui est réellement engendrée et induite dans la terre ne ressemble pas souvent au signal de source idéalisé le praticien tente d'induire dans la terre. Des effets non linéaires dans la terre et dans le système de production d'énergie de la source, par exemple un vibreur sismique, peuvent créer des harmoniques puissants et d'autres phénomènes dynamiques qui introduisent des distorsions imprévisibles dans le signal réel engendré induit dans la terre. Ces effets non linéaires ne sont pas prévisibles et dépendent de l'emplacement de la source. L'intensité et les caractéristiques des distorsions non linéaires varient d'une position balayage de source à la suivante. Ces effets non linéaires interfèrent souvent avec la déconvolution, suffisamment pour dégrader sensiblement les images acoustiques. Il serait très souhaitable de trouver un procédé qui décompose de manière précise le signal à balayage de fréquence réellement engendré en ses diverses composantes, en permettant une déconvolution plus efficace et un autre filtrage. Il y aurait bien des avantages à décomposer de manière précise les signaux sismiques en leurs éléments harmoniques composants. La déconvolution, une technique qui est fonction l'analyse précise des signaux, serait beaucoup plus précise. Divers harmoniques séparables, ou diverses composantes de bruit, de signaux sismiques peuvent être renforcés, isolés ou atténués en utilisant des filtres appropriés. exemple, le bruit instrumental, le bruit extrinsèque et d'autres bruits insensibles au filtrage classique peuvent, dans bien des circonstances, être ciblés de manière précise, comme résultat de ce procédé de décomposition. La présente invention propose un procédé, destiné à décomposer et à synthétiser un signal sismique en composantes harmoniques et autres composantes extrinsèques, qui est précis à des instants discrets ou en des points discrets de données d'échantillon. On forme un ensemble de bases de fonctions de composantes de signal en utilisant l'amplitude et la phase instantanées de la source pilote et des composantes associées signal observées dans des données réelles, par exemple, un ensemble de bases comprenant des harmoniques balayage pilote. Ces bases s'utilisent pour former des fonctions complexes de modulation pour chaque composante ou harmonique de signal à analyser et à synthétiser par la présente invention. Le signal sismique qui fait 1 et de l'étude est utilisé pour moduler la fonction de modulation complexe obtenue et la sortie est filtrée passe- . Cette sortie filtrée passe-bas sutilise pour calculer l'amplitude instantanée et la phase différentielle composante de signal en cours d'analyse. À partir de ces calculs, on synthétise les composantes particulières signal qui sont ciblées. Ces données de composante de signal synthétisées sont alors isolées, par soustraction, données cibles d'origine, en laissant les séries temporelles résultantes contenant les composantes restantes de signal. Chaque composante de signal est ensuite synthétisée de la même manière, extraite et ainsi isolée à son tour des données restantes. Le résultat est une décomposition et une synthèse précises du signal sismique d'origine en ses diverses composantes de signal. Les caractéristiques et avantages de l'invention ressortiront d'ailleurs de la description qui va suivre, à titre d'exemple, en se référant aux dessins annexés, dans lesquels : figure 1A représente un système d'enregistrement sismique typique pour acquisition de données ; figure 1B représente un organigramme pour un mode de réalisation préféré de l'invention ; figure 2 représente une forme d'onde typique, avec transition, de source à balayage ou fluctuation de fréquence sans distorsion, associée au balayage pilote ; figure 3 représente une forme d'onde typique à balayage ou fluctuation de fréquence avec distorsion, associee à une énergie réelle induite dans la terre ; figure 4 représente des composantes réelles de bases, typiques, du balayage pilote de la figure 2 figure 5 représente la sortie de la modulation, à l'aide de la fonction de modulation complexe obtenue du balayage déformé de la figure 3, pour le premier harmonique provenant de la figure 4 ; la figure 6 représente les séries temporelles la figure 5 après filtrage passe-bas ; et figure 7 représente la composante de premier harmonique, synthétisée, du signal sismique cible de la figure qui a été calculée en utilisant la sortie la figure La présente invention décrit un nouveau traitement destiné, à partir de données sismiques, à analyser, décomposer, synthétiser et extraire des composantes de signaux sismiques, comme la fondamentale d'un balayage pilote ou ses harmoniques. Cette synthèse se fait en modulant, à l'aide des données qui sont analysées des fonctions de bases obtenues d'un balayage pilote. Les données résultantes sont filtrées passe-bas de façon à pouvoir calculer l'amplitude et la phase différentielle instantanées de la composante particulière de signal à synthétiser. Chaque composante ciblée est alors synthetisée et le résultat est soustrait des données restantes. La présente invention est basée sur le concept que n'importe quelle fonction réelle R(t) peut être représentée par un signal analytique complexe Z (t) constitué de la fonction d'origine en quadrature avec sa transformée Hilbert H [R(t)] I(t) de sorte que Z(t) R(t) + jI(t), où et R(t) et I(t) sont déphasés de n/2 radians à toutes les fréquences, et sont donc orthogonaux. description de la transformée de Hilbert peut se trouver dans Oppenheim. A., et Schafer R., "Discrete-Time Signal Processing'l, Prentice-Hall, Englewood Cliffs, New Jersey, 1975, chapitre 10. ,amplitude A(t), la phase (D(t) et la fréquence 0(t), instantanées, sont décrites dans Taner, M., Koehler, F., Shériff, R., I'Complex seismic trace analysisl' Geophysics, volume 44, numéro 6, juin 1979, pages 1041 a 1063, sont définies comme : A(t) [R2(t) + 12(t)]1/2# (2) (D(t) tan-1 [I(t)/R(t)], et (3) E) (t) (D 1 (t) d [(D (t) 1 /dt (4) - partir des définitions ci-dessus, on peut facilement démontrer que : (t) A (t) cos [(D (t) (5). Pour la présente invention, on modélise un signal (sinusoïdal) à bande étroite avec une amplitude et fréquence variant "lentement". On suppose que l'on définit l'amplitude instantanée ou moyenne comme étant constante plus une certaine perturbation, en tant que fonction du temps, ainsi A(t) ao + a(t) (6) et de f açon similaire pour la phase instantanée ou moyenne (D (t) #O + # (t) (7). En particulier, dans un mode de réalisation de l'invention, la phase représente un signal à fluctuation ou balayage linéaire défini par : C01 = (Ù)T - (00) /T, (9) et 0 - t :9 T. CI est un balayage commençant à wo à t 0, et se terminant à COT à t T. À partir des équations (4) et (8), fréquence instantanée est une fonction linéaire et elle est donnée par : (t) (0() + colt En outre, cette analyse sera limitée au cas où 0 (t) varie lentement avec chaque échantillon, c'est-à-dire : 91 (t) At 1 où At est l'intervalle d'échantillon. À partir de l'équation (10), on voit facilement que la dérivée de 0(t) par rapport au temps est 01(t) (ol. Pour un signal de balayage sismique typique échantillonné à 500 échantillons/sec- en commençant à 10 Hz et en terminant à 90 Hz et en retardant de 8 sec., on a col At 27c(90 - 10)/(8 . 500) 7z/25 radians par échantillon. À partir des équations (5) et (8), nous pouvons maintenant définir un balayage linéaire modulé à amplitude idéale comme : P (t) A (t) cos [(D (t) (11), et avec une certaine distorsion ajoutée modélisée par le kéme harmonique à balayage sans transition : Pk (t) cos [k (V (t) 1 (12), où A (t) est la fonction d'amplitude instantanée de modulation ou ce que l'on appelle une fonction avec transition. Cette fonction avec transition habituelle ment de la forme : (t) #O # 1 + cos [n (l - t /T) 2, 0 :#9 T, (13) (t) t T T, et (14) A (t) () (l cos [ir (l t/T) /2, T - c :## t :9 T, (15) où T est la largeur de la région de transition de début et de fin. Une telle transition est appelée transition cosinus. on pourrait utiliser d'autres transitions qui comprennent une transition linéaire et une transition de Blackman : de telles transitions devraient être connues de ceux qui sont versés dans 1 art et ne seront pas décrites davantage. De plus, il n'est nécessaire pour la mise en pratique de la présente invention que la transition au début du signal pilote soit même que la transition à la fin du signal pilote. Dans un mode de réalisation préféré de l'invention, on analyse un signal S(t) de balayage déformé de vibreur en utilisant, comme référence, signal P(t) d'entrée sans transition ou de balayage pilote. Le balayage déformé S(t) peut être approché par Sn( constitué d'une somme de n harmoniques modulés individuellement en amplitude du balayage pilote sans transition, et pour le'kème harmonique ayant chacun une petite fonction 80t) d'erreur de phase ajoutée à la phase instantanee (Dk(t) d'origine. Ceci peut s'exprimer comme : Ceci peut aussi être Où Sn (t) S (t) (Dk (t) représenté en utilisant la forme quadratique Pour le kème harmonique reste à calculer (Xk(t) : son amplitude instantanée, 8k(t) : sa phase différen tielle instantanée. Pour le cas ordinaire d'analyse du signal pilote P(t) non déformé, on aurait seulement des valeurs pour la fondamentale, c'est-à-dire, n Dans ce cas, on analyserait facilement (t) en ce qui concerne son amplitude A(t) et sa phase (D(t), instantanées. Par conséquent, à partir de llequation (16), on aurait #J (t) (t) et 81 (t) 0 puisque P (t) (t) cos (D (t) Cependant, afin d'analyser un balayage S(t) déformé ayant de nombreux harmoniques, l'analyse directe instantanée d'amplitude et de phase ne fonctionnera pas en raison de l'effet d'interférence à large bande des harmoniques mêmes. présente invention isole chaque harmonique, autre composante, choisi de signal pour produire une amplitude instantanée et une phase différentielle pour cet harmonique, ou autre composante, de signal. On utilise un filtre hétérodyne en quadrature pour décaler chaque harmonique, ou chaque autre composante, de signal vers le bas façon que le kème harmonique instantané d'origine ou la phase de composante de signal (Dk(t) Ceci s'obtient pour harmoniques d'ordre plus élevé en modulant le balayage déformé S(t) à l'aide d'une sinusoïde complexe telle définie par : Hk(t) S(t)#c0s[(Dk(t)] + j sinl(Dk(t)]) (18). L'effet de cette modulation peut se voir substituant Sn(t) à S(t) de l'équation (16) dans léquation (18) en donnant : Celle-ci peut être développée pour donner : regardant maintenant seulement le cas où i cela donne : . -k(t)= #k(t)#(cOsl3k(t)] + cos[2(Dk(t) + Ôk(t)ll + j(sin[2(Dk(t) + 80t)] - sin[8k(t)11/2) (21). Puisque 8k(t) est petit par rapport à (Dk (t) on limite bande du signal complexe Hk(t) en utilisant un filtre passe-bas. Ceci maintient les termes de Fourier près d'une frequence nulle (courant continu) et rejette les termes plus grands de [2(Dk (t) 1 . Cette version f iltrée passe-bas de Hk(t) peut se définir comme : fik(t)= :Pk(t) + jik(t) & k(t)#c0s[8k(t)] j sinlbk(t)])/2, (22). où composantes réelle et imaginaire sont montrées de manière explicite. on a aussi supposé que, puisque l'amplitude instantanée varie légèrement, ix-k(t) (Xk(t) On obtient maintenant les fonctions voulues uk(t) et Ôk(t). Comme pour l'amplitude instantanée et la phase différentielle à partir de l'équation (22) : #k (t) #_# 2 1 Hk (t) 2 [:Ê k 2(t) + Îk 2(t)]1/2 et Ôk(t) z-t tan-' lik(t)/ f#k(t)] (23). En ce qui concerne notre représentation en quadrature de l'équation (17) nous avons Pk (t) 2:#, k (t) , et _Yk (t) 2 1 k (t) (24). Ces approximations sont valables seulement si fonctions de distorsion varient légèrement et si coo, la fréquence de balayage de départ, n'est pas trop petite. Cette dernière restriction est nécessaire pour éviter 1 interférence précoce d'harmoniques dans un signal de balayage ascendant. Dans un mode de réalisation préféré de l'invention, la coupure de filtre passe-bas est entre 1 % et 10 % de la fréquence des Nyquist. Le filtrage passe-bas se fait en utilisant des TFR (transformées de Fourier rapides) dans le domaine fréquenciel. La largeur d'une zone de transition d'environ 25 % de la valeur de coupure du filtre passe-bas est efficace pour réduire les effets de suroscillation du filtrage. L'interférence d'harmoniques est réduite en démarrant l'analyse à k = 1 et en allant jusqu'à k = n. Chaque harmonique peut être synthétisé à partir soit de Sk(t) cck(t) cOs[(Dk(t) + 80t)], soit de Sk(t) Pk(t) cOsl(Dk(t)] + Yk(t) sin [(Dk(t)] (25) et du résultat soustrait du signal de balayage déformé restant. D'autres fonctions de base à côté des harmoniques multiples entiers peuvent également être isolées et soustraites avec ce procédé, par exemple des sous- harmoniques fractionnaires rationnels, des ultra-sous- harmoniques ou d'autres oscillations dynamiques du signal de source de référence pilote. On prend soin de garantir que tous les retards temporels entre le balayage pilote d'entrée P(t) et le balayage de sortie déformé S(t) soient éliminés ou compensés d'une autre manière. Ceci peut se faire en utilisant une corrélation croisée entre balayage et le pilote, et en recherchant alors la valeur de pic. L'emplacement du pic doit être à retard ou sinon les signaux doivent être déphasés de façon appropriée pour garantir cette condition. La précision l'analyse de distorsion peut se mesurer en calculant le résidu de synthèse pour la somme des n premiers harmoniques en utilisant l'équation (16) et elle est définie par : En (t) (t) - Sn (t) (26) en évaluant alors la norme vectorielle relative en IlEn (t) Il / IIS (t) Il Ainsi, en est une mesure de l'erreur de la synthèse du balayage déformé en utilisant les n premiers harmoniques autres composantes de signal provenant de l'analyse. Les données sont acquises au moyen la production, une source, d'un signal sismique approprié, comme le montre la figure 1A. Une source appropriee peut être un vibreur sismique qui engendre des ondes acoustiques sismiques qui se propagent dans la terre. Les vibrations sont enregistrées au niveau d'un récepteur, ou de récepteurs, appropriés comme un géophone ou des géophones, enterrés dans le sol ou en surface, et l'on peut alors les analyser en ce qui concerne les composantes de signal auxquelles on s'intéresse. L'exemple de source de la figure 1A un vibreur 30 monté sur camion pour des relevés de terrain mais d'autres combinaisons sont bien connues dans l'art, exemple un vibreur marin pour l'acquisition de données dans des zones recouvertes d'eau. Les vibrations acoustiques du signal sortant de la source peuvent être enregistrées pour analyse sur la source elle-même comme c'est habituel dans la pratique de l'art, ou très près de la source avec un géophone 32 enterré, dans certains programmes d'acquisition sur le terrain, avec un hydrophone proche de la source dans relevés marins ou n'importe quel autre capteur approprié proche de la zone, ou avec des combinaisons de capteurs. Les vibrations acoustiques se propagent, dans la terre, depuis la source 30 le long de trajets sismiques rayonnants 34 jusqu'à un récepteur ou des récepteurs par exemple un géophone 36. Des données enregistrées un certain récepteur approprié, sur la source elle-même dans zone immédiate de la source ou après avoir parcouru une certaine distance dans la terre, peuvent être analysées a 'aide de la présente invention. La figure 1B présente un organigramme montrant un mode de réalisation préféré de l'invention. Pour la présente invention, un mode de réalisation propose un vibreur sismique propre à engendrer un signal acoustique arbitraire sur une durée de temps et une plage de fréquences prédéterminées. Ce signal acoustique engendré est alors appelé le balayage pilote et il contient des caractéristiques de balayage prédéterminées choisies le praticien. La figure 2 montre un exemple d'un balayage pilote. En pratique, l'énergie acoustique engendrée un vibreur sismique contient de nombreuses autres composantes signal à côté des caractéristiques de balayage pilote predéterminées. La figure 3 montre un exemple d'un balayage reel ou déformé mesuré par un récepteur comme celui montré à figure 1A, contenant ces autres caractéristiques de signal. C'est sur l'analyse et la synthèse ces composantes de signal que porte la présente invention. La première étape consiste à déterminer phase instantanée (D (t) du signal P(t) de balayage pilote d'entrée montré en 62 la figure 1B et, dans le cas habituel, on pose k = 1. La fonction de modulation complexe de la composante de signal à synthétiser s'obtiendra du balayage pilote P(t). La figure 2 montre un balayage pilote typique. Pour la représentation sur la figure, et non pas à titre de limitation, on a représenté seulement une plage limitée de fréquences. Les balayages pilotes des données réelles sont, généralement, entre 5 et 200 Hertz, et peuvent aller plus haut. Dans le cas habituel, la fondamentale (ou premier harmonique) ser a la première composante de signal à être synthétisée, bien que n'importe quelle composante de signal pourrait être synthétisée en premier. L'étape suivante consiste à former la fonction de modulation complexe pour le filtre hétérodyne pour la composante de signal à analyser, comme montré en 64 de la figure 1B, par exemple dans le cas du kI#me harmonique : (COS l(Dk (t) + j sin I(Dk (t) 1 . La figure 4 montre des exemples de bases typiques d'un balayage pilote. la figure 4, on a montré seulement la composante réelle des bases, et ici encore à des fins de représentation de la figure, on a montré seulement des plages limitées de fréquences. Ceux qui pratiquent l'art comprendront que les plages de fréquences de balayage réelles et leurs bases peuvent couvrir plusieurs octaves. Les données dans les données réelles déformées sorties du vibreur, le signal sortant ou de champ proche du champ d'ondes acoustique 52, peuvent être acquises à partir de n'importe lequel des signaux habituels mesurés dans, ou autour du vibreur, ou acquises par n'importe quel capteur, récepteur ou géophone approprié. Les signaux typiques de sortie de vibreur sont des accélérations de masses ou des forces de plaque d'appui et de réaction, des signaux d'autres capteurs de champ proche ou d'une combinaison quelconque. Les signaux de champ proche peuvent se mesurer en utilisant des accéléromètres, des géophones ou d'autres capteurs et combinaisons, comme il convient. Les données ciblées sont alors filtrées comme montré en 66, par exemple le balayage déformé de 52 pour le premier cas lorsque 54 s'applique, par la fonction modulation complexe dérivée, de 64. Dans le cas habituel la première synthèse entraînerait la multiplication de fonction de base sans transition du balayage pilote fondamental par les données de balayage déformées. première série temporelle à la figure 4 (dont la partie imaginaire n'est pas représentée à la figure 4) serait multipliée par la série temporelle de la figure 3. Les séries temporelles complexes résultantes contiennent l'amplitude instantanée uk(t) et la phase instantanée différentielle Ôk(t) auxquelles on s'intéresse. La figure 5 montre un exemple de la sortie de cette opération. Les composantes d'amplitude et de phase de composante d harmonique ou de signal que l'on synthétise sont calculées à partir de cette série temporelle apres application un filtrage passe-bas approprié comme montré en 68. La figure 6 montre une représentation de série temporelle la sortie de filtre passe-bas. Lorsque les composantes de phase et d'amplitude été déterminees comme en 70, on synthétise alors comme en 72 les données ciblées de composantes d'harmonique d'autre signal. Cette synthèse, dans le cas du kème harmonique, peut s'obtenir comme dans l'équation (25) dessus. La figure 7 montre la composante synthétisée de la fondamentale (ou premier harmonique) qui a été extraite balayage déformé d'origine. Après synthèse, on soustrait ces données du balayage déformé en . Dans le premier cas habituel, on soustrait du balayage déformé, la composante synthétisée du balayage pilote associé à la fondamentale. Après les premières synthèse et soustraction, on choisit, 76 avec 78, la prochaine composante de signal à analyser. Dans le cas habituel, on synthétisera, avec ce procédé, le deuxième harmonique du balayage pilote. On peut synthétiser des composantes de signal ultérieures, qui peuvent être des harmoniques d'ordre plus élevé, des sous- harmoniques ou un autre comportement oscillatoire dynamique. Après avoir soustrait, de la série temporelle restante, chaque composante du signal déformé origine, on peut analyser le résultat par le procédé présenté ci-dessus dans les équations (26) et (27) jusqu'à ce 1 on ait atteint un niveau arbitrairement petit ou acceptable. L'invention a été décrite ci-dessus en ilisant un exemple dans lequel la fréquence du signal est fonction linéaire temps. Le procédé de la présente invention peut aussi s'utiliser lorsque la fréquence est fonction continue temps et il est bien entendu qu'une telle variante tombe à l'intérieur de la portée de l'invention. REVENDICATIONS 1. Procédé de prospection sismique utilisant une source vibratoire (30), le procédé étant caractérise en ce qui comprend, en combinaison : (a) l'activation de la source vibratoire en utilisant un signal pilote (P(t)) de référence, ledit signal pilote (P(t de référence ayant une frequence initiale à un instant initial et une fréquence finale à un instant final ; (b) la représentation d'une phase instantanée du signal pilote (P(t de référence en fonction desdites fréquences initiale et finale et desdits instants initial et final ; (c) la mesure du signal sortant ou de champ proche de ladite source vibratoire (30) ; (d) la définition d'une pluralité de signaux de base chacun desdits signaux de base étant une sinusoïde ayant une fréquence linéaire et une phase différentielle, ladite fréquence linéaire étant un multiple d'une frequence instantanée du signal pilote (P(t de référence ; et (e) la représentation dudit signal sortant ou de champ proche, mesuré, en tant que somme des signaux de base pondérée par une fonction associée de pondération. 2 Procédé selon la revendication 1, caractérisé en ce que la fonction associée de pondération s'obtient par un processus qui comprend en outre la définition d'un filtre hétérodyne en quadrature pour au moins l'un de la pluralité de signaux de base. 3. Procédé selon la revendication 2, caractérisé en ce qu'il comprend en outre : (i) l'application d'au moins l'un desdits filtres hétérodynes en quadrature au signal sortant ou champ proche mesuré pour donner un signal sortant modulé - (ii) le filtrage passe-bas du signal sortant modulé et (iii) la détermination de la fonction pondération de la phase différentielle pour au moins l'une de la pluralité de fonctions de base provenant du signal sortant de champ proche, modulé, filtré passe-bas. 4. Procédé selon la revendication caractérisé ce qu'au moins l'un de la pluralité de signaux de base comprend en outre une pluralité ordonnée d'harmoniques du signal pilote (P(t de référence, le procédé comprenant en outre le fait de soustraire de manière itérative, du signal sortant ou de champ proche mesuré, la fonction de base pondérée comprenant la phase différentielle correspondante. 5. Procédé selon la revendication-1 caractérisé en ce que ledit multiple de phase linéaire la fonction de base est choisi entre (i) un multiple entier et (ii) un multiple fractionnaire rationnel. 6. Procédé selon la revendication caractérisé en ce qu'il comprend en outre l'utilisation d'une transition de l'amplitude du signal pilote (P(t)) de référence. 7. Procédé selon la revendication caractérisé ce que la transition est choisie dans le groupe constitué (i) d'une transition cosinus, (ii) d'une transition de Blackman et (iii) d'une transition linéaire. 8. Procédé selon la revendication caractérisé ce que la transition au début du signal pilote (P(t de reférence est différente de la transition à la fin du signal pilote (P(t de référence. 9. Procédé de prospection sismique utilisant une source vibratoire (30), le procédé étant caractérisé en ce 'il comprend, en combinaison : (a) l'activation d'une source vibratoire (30) en utilisant un signal pilote (P(t de référence, ledit signal pilote (P(t de référence ayant une fréquence initiale à un instant initial et une fréquence finale à un instant final ; (b) représentation d'une phase instantanée du signal pilote (P(t)) de référence en fonction desdites fréquences initiale et finale et desdits instants initial et final ; (c) mesure d'un signal sortant ou de champ proche de ladite source vibratoire (30) ; (d) définition d'une pluralité de signaux de base, chacun desdits signaux de base étant une sinusoïde ayant une frequence variant en continu avec le temps et une phase différentielle, ladite fréquence étant un multiple d'une fréquence instantanée du signal pilote (P(t de référence ; et (e) représentation dudit signal sortant ou de champ proche mesuré, en tant que somme des signaux de base, pondéree par une fonction associée de pondération. 10 Procédé selon la revendication 9, caractérisé en ce que la fonction associée de pondération s'obtient par un processus qui comprend en outre la définition d'un filtre hétérodyne en quadrature pour au moins l'un de pluralité de signaux de base. il. Procédé selon la revendication 10, caractérise en ce qu'il comprend en outre : (i) application d'au moins l'un desdits filtres hétérodynes en quadrature au signal sortant ou de champ proche mesuré pour donner un signal sortant modulé ; (ii) le filtrage passe-bas du signal sortant modulé et (iii) la détermination de la fonction de pondération et de la phase différentielle pour au moins l'une de pluralité de fonctions de base provenant du signal sortant modulé, filtré passe-bas. 12. Procédé selon la revendication il, caractérise en ce qu'au moins l'un de la pluralité de signaux de base comprend en outre une pluralité ordonnée d'harmoniques signal pilote (P(t)) de référence, le procédé comprenant en outre le fait de soustraire de manière itérative, du signal sortant mesuré, la fonction de base pondérée comprenant la phase différentielle correspondante. 13. Procédé selon la revendication 9, caractérisé en ce que ledit multiple est choisi entre (i) un multiple entier et ii) un multiple fractionnaire rationnel. 14. Procédé selon la revendication 9, caractérisé en ce il comprend en outre l'utilisation d'une transition de l'amplitude du signal pilote (P(t) de référence. 15. Procédé selon la revendication 14, caracterisé en ce la transition est choisie dans le groupe constitué (i) d'une transition cosinus-, (ii) lune transition de Blackman et (iii) d'une transition linéaire. 16. Procédé selon la revendication 14, caracterisé en ce que la transition au début du signal pilote (P(t de référence est différente de la transition à la fin du signal pilote (P(t de référence.