Brevet FR2757641A1 — publié le 19980626

La présente invention concerne une méthode de transmission de données utilisant des techniques de compression en vue d'optimiser I'utilisation de voies de transmission disponibles. La méthode selon l'invention trouve des applications notamment dans le domaine de la prospection sismique où il est nécessaire de transférer vers une station centrale telle qu'un camion-laboratoire, une masse souvent considérable de données. Des signaux sont captés par un très grand nombre de récepteurs tels que des géophones disposés couplés avec une formation géologique à étudier, en réponse à des ébranlements émis par une source sismique et renvoyés par les discontinuités du sous-sol. Les signaux captés sont collectés par des unités locales d'acquisition réparties parfois sur une distance de plusieurs kilomètres et destinées chacune à collecter les signaux reçus par un ou plusieurs récepteurs, à les numériser, à leur appliquer des traitements préalables plus ou moins complexes et les stocker dans une mémoire locale avant leur transmission en temps réel ou différé à une station de collecte par une voie de transmission telle qu'un câble, une fibre optique, un canal radio etc ETAT DE LA TECHNIQUE Différents systèmes de transmission de données sismiques sont utilisés pour relier des unités locales d'acquisition à une station central soit directement, soit via des stations intermédiaires pourvues de fonctions plus ou moins complexes de concentration ou de contrôle des unités locales. Les liaisons peuvent être assurées au moyen de câbles, de liaisons radio, via un ou plusieurs relais éventuels. ou encore combiner les liaisons par câbles et par liaison radio comme indiqué par exemple dans les brevets FR2720518. FR 2 696 839, 2.608.780, 2.599.533, 2.538.561,2.511.772 ou 2.627.652 du demandeur. Par le brevet FR-A-2 608 780 du demandeur notamment, il est connu d'utiliser des boîtiers d'acquisition sismique dotés de deux voies de transmission, I'une à débit de transmission relativement élevé, I'autre à bande passante qui peut être relativement étroite selon la disponibilité locale des fréquences de transmission. plus facilement disponibles dans le cadre des réglementations d'émission radioélectriques en vigueur. Les données sismiques collectées au cours des cycles successifs sont stockées sur une mémoire de masse dans chacun des boîtiers ct transférées par intermittence à une station centrale de commande et d'enregistrement. Pour permettre à l'opérateur dans la station centrale de vérifier que l'acquisition des données par chacun des boîtiers d'acquisition se déroule normalement, on procède à une transmission des données partielles ce qui s'accommode bien d'une voie dc transmission à bande passante relativement étroite. Par le brevet FR-A-2 615 627 du demandeur également. il est connu d'utiliser des boîtiers d'acquisition sismique pourvus notamment de moyens de traite nient spécialisés dans le traitement des signaux auxquels on fait réaliser de nombreux contrôles des géophones et éléments de la chaîne d'acquisition ainsi que des pré-traitements des traces sismiques auparavant réalisés à la station centrale après transmission. permettant donc de diminuer considérablement le volume dc données à rapatrier. La tendance actuelle, notamment dans le cadre des méthodes d'exploratioii sisiiiique dites 3D. est à répartir sur une zone à explorer, à terre, en mel, ou dans des zones côtières, souvent sur plusieurs kilomètres, des récepteurs sismiques par centaines, voire par milliers. Le volume des données à collecter et à ti.-insnietti.e ne cesse de croître. Pour éviter que les problèmes de transmission ne constituent um frein à l'évolution des systèmes sismiques la tendance est de recourir à des procédés dc coiiipression de données choisis pour être compatibles avec les cxigenccs propres des géophysiciens. La compression des données sismiques peut apporter à la fois un gain dc place appréciable sur les modules de stockage cie masse clans les boîtiers locaux d'acquisition et/ou les stations locales de contrôle et de concentration, et aussi un pain considérable en temps de transmission. De nombreuses méthodes existent pour compresser des données. On peut ranger en deux grandes classes selon qu'elles entraînent ou non une perte d'information, que l'on choisit suivant que l'on peut ou non tolérer une altcration des données après leur décompression. En géophysique notamment. il est indispcnsable que les pertes dues à la compression restent le plus faible possible car les informations les plus pertinentes sont souvent de très faible amplitude et ne peuvent être isolées du bruit de fond que par un traitement numériques effectué sur plusieurs traces. Des pertes éventuelles de précision ne sont tolérables que dans certains cas bien particuliers. si l'information transmise ne doit servir qu'à contrôler le bon fonctionnement du matériel et à visualiser "l'allure" des traces échantillonnées. Les méthodes connues de compression de données peuvent être classées en deux familles, a) les méthodes de compression sans perte d'information et b) les méthodes conduisant à une perte d'information, les données restituées perdant une part plus ou moins grande de leur précision. a) Parmi les méthodes de la première famille, on peut citer les méthodes qui visent à éliminer la redondance des données ou les méthodes dites à dictionnaire où chaque mot est remplacé par son index dans une table de référence, qui sont d'autant plus intéressantes que les fichiers à compresser contiennent beaucoup de redondance. Une méthode également connue dite RLE pour " Run-Length Encoding" est bien adaptée aux fichiers contenant d'importantes séquences de valeurs identiques. On connaît aussi les méthodes de codage de type statistique où l'on cherche à remplacer les données par un code ayant la même signification mais occupant moins de place. La méthode de codage dite de Huffman par exemple consiste à associer à une donnée un code de longueur variable en fonctions de sa fréquence d'occurrence. Les méthodes de type à codage arithmétique visent à représenter un nombre variable de données par un nombre constant de bits. Une technique connue de compression dite LPC pour Linear Predictive Coding est bien adaptée à la compression des ondes acoustiques. Elle consiste essentiellement à remplacer- un échantillon de signal s(t) par une prédiction faite à partir de p échantillons précédents, en supposant que le signal est stationnaire (les propriétés statistiques du signal sont indépendantes du temps) et centré (d'espérance mathématique nulle), avec des coefficients de prédiction indépendants du temps. Au lieu de transmettre l'échantillon s(t), on transmet sa prédiction sAt) c'est-à-dire les coefficients de prédiction et les résidus e(t) c'est-à-dire l'écart entre la valeur réelle et la prédiction qui en est faite au temps t. ce qui permet de retrouver à la décompression la valeur s(t) = sA(t) + e(t) ). Si la prédiction est bonne. les résidus sont faibles et occupent une place moindre que celle des valeurs initiales s(t). Dans la famille des méthodes de compression avec pertes où l'objectif est de modéliser au mieux un signal. on peut citer par exemple la méthode par quantification où l'on représente une information sur n bits par une autre codée seulement sur p (p Suivant d'autres méthodes telles la "transformée en cosinus discrète" ou DCT (Discrete cosine transform), utilisée très couramment utilisée pour la compression d'images, on découpe le signal en blocs de taille fixe, puis on modélise chacun d'eux par un ensemble de coefficients qui ne représentent pas l'amplitude du signal dans le temps mais plutôt son spectre fréquentiel. La compression apparaît en éliminant les coefficients de plus haute fréquence. D'autres méthodes basées sur les transformées en ondelettes, sont également utilisées et on les trouve décrites par exemple par - Bosman. C. et al, Seixrnis dati conipression Lusin wavelet trzzttxJ{Jn7lx, 83 th Annual meeting SEG, - Vetterli. M.et al, Waielets cind filter bans, IEEE Transactions on signal processing, vol 40, No 9, Septembre 1992 - Coifman, R. et al; Best adLlpted ve rocket bases,, Numerical Algorithms Research Group, Yale University, ou - Daubechies, I. Ten Lectures on wavelets, CBMS-NSF, 1993. Un signal quelconque est représenté clans une base de fonctions qui sont les versions translatées (temporel) et dilatées (fréquence) d'une ondelette. L'objectif désiré est de modéliser un signal en cherchant où peuvent être localisées (dans le temps) les irrégularités du signal (par exemple : des variations brusques ou hautes fréquences), représenter l'ondelette choisie pour appliquer la transformée. On peut choisir parmi plusieurs familles d'ondelettes d'un type connu chacune avec des caractéristiques qui lui sont propres. La méthode de transmission selon l'invention comporte l'utilisation de techniques de codage connues en elles-mêmes, choisies de façon à pouvoir intégrer souplement les phases de transmission de données dans le déroulement classique d'une campagne de prospection sismique aussi bien dans les zones purement terrestres, que dans les zones côtières où les données sismiques collectées par des appareils d'acquisition peuvent être le résultat d'ondes renvoyées par les discontinuités du sous-sol, provenant de sources sismiques remorquées en immersion par un navire, et par exemple de canons à eau ou à air. Dans le premier cas, on utilise par exemple des charges explosives. A chaque fois, l'équipe sur le terrain en installe plusieurs, et les intervalles de temps entre les tirs successifs de ces différentes peuvent être relativement courts, réduits sensiblement à la fenêtre de réception utile, quelques dizaines de seconde en pratique. Entre chaque "salve", il s'écoule un temps plus long, nécessaire au positionnement d'un nouvel ensemble de charges par l'équipe de terrain. Dans le deuxième cas, le navire remorqueur évolue en zigzag suivant des lignes parallèles alternativement dans un sens et dans le sens opposé. Une nouvelle série de cycles d'émission-réception-acquisition a lieu Ic long de chacunc des trajectoires successives. à intervalles relativement courts dépendant du temps de oe'armement de la source utilisée. Entre les séries successives dc cycles, il s'écoule un temps relativement long nécessaire au navire pour virer de bord et se repositionner suivant une nouvelle trajectoire. La méthode de transmission selon l'invention est adaptée à tirer au mieux parti de cette distribution irrégulière des temps forts de collecte de données sismiques dans la pratique usuelle des campagnes d'exploration sismique. pour réaliser aussi bien des transmissions de contrôle destinées à vérifier le bon déroulement des opérations d'acquisition par chacun des appareils d'acquisition sur le terrain, que des rapatriements complets des données sismiques à une station centrale de commande, de façon à restituer sans perte toutes les traces sismiques transmises, en limitant au mieux les durées de transmission nécessaires à cet effet. La méthode selon l'invention convient pour la transmission de données entre une ou plusieurs unités locales de collecte de données sismiques telles que des traces ou portions de traces sismiques par exemple. et une station centrale de commande. par au moins une voie de transmission. Elle est caractérisée en ce qu'elle comporte une transmission en deux temps des données sismiques avec, dans un premier temps, une transmission de données compressées permettant un contrôle de la qualité de fonctionnement des unités locales et, dans un deuxième temps, une transmission de données permettant un enregistrement sans perte des données sismiques à la station distante. La méthode comporte par exemple l'enregistrement local des données, une compression des données avec un taux de compression choisi, suivie d'une première étape de transmission de ces données compressées, une détermination dc la différence entre les données enregistrées localement et les données correspondantes compressées, une deuxième étape postérieure de transfert de cette différence, et une reconstitution à la station centrale des données sismiques collectées par chaque unité locale de collecte, par combinaison des données reçues au cours des deux étapes de transmission et dc transfert. On choisit la technique de compression de données et/ou le taux de coinpression à appliquer aux données en fonction par exemple de la durée d'une lenêtre cte transmission disponible et du débit de transmission permis par la voie de transmission utilisée. Suivant un mode de réalisation, les données sismiques collectées sont enregistrées localement après une compression n'impliquant sensiblement aucune perte et la partie significative est transmise avec un taux de compression impliquant une perte. La méthode selon l'invention convient par exemple pour le transfert à la station centrale par au moins une voie de transmission de traces sismiques collectées par des unités d'acquisition locales au cours de sessions d'acquisition de données comportant chacune une pluralité de cycles d'enregistrement de signaux sismiques renvoyées par des discontinuités du sous-sol, en réponse à des signaux sismiques transmis dans le sol par une source sismique, les différents cycles d'une même session étant séparés les uns des autres par des intervalles de temps déterminés. Les sessions sont par exemple des journées de travail. Les intervalles de temps entre les sessions peuvent aussi être par exemple dus à des interruptions techniques nécessaires pour modifier ou déplacer l'équipement de terrain. ou bien réaligner un bateau boute-feu suivant un nouvel axe d'évolution par exemple. La méthode peut être mise en oeuvre par exemple en réalisant une première compression de chaque trace sismique avec un premier taux de compression suffisant pour qu'elle puisse être transmise durant un des dits intervalle dc temps I'elllegistlement dans une unité locale d'acquisition local de chaque trace sismique compressée sensiblement sans perte avec un deuxième taux de compression, la détermination de la différence entre les traces sismiques compressées respectivement avec le premier et le deuxième taux, le transfert différé de cette différence durant un intervalle de temps entre les dites sessions, soit en utilisant une tleuxièine voie de transmission (câble, fibre optique, canal radio etc.), soit en déplaçant éventuellement une mémoire de masse jusqu'auprès des unités de collecte pour venir lire Icur-s modules de mémorisation), et la reconstitution à la station centrale de chaque trace sismique . dans le but par exemple de contrôler le bon fonctionnement de chaque unité de collecte locale. La première compression est effectuée par exemple en appliquant une technique de compression avec perte en utilisant par exemple une technique dc transformation par ondelettes avec un taux impliquant des pertes. Pour obtenir des données compressées sans perte, on peut appliquer par exemple aux données sismiques une technique de prédiction statistique ou de transformation par ondelettes. Le système pour la mise en oeuvre de l'invention est caractérisé en ce que chaque unité locale de collecte comporte des moyens d'enregistrement de données. un ensemble de traitement comprenant des moyens de compression des données avec un taux de compression choisi, des moyens pour déterminer la différence entre les données enregistrées localement et les données correspondantes compressées, des moyens de transmission des données compressées, des moyens de transfert de la dite différence à la station, et un ensemble de traitement dans la stations distante comprenant des moyens pour reconstituer les données sismiques collectées par chaque unité locale de collecte. par combinaison des données transmises et des données transférées, et éventuellement des moyens dans chaque unité locale de collecte pour appliquer aux données sismiques une deuxième compression de type sans perte avec un deuxième taux de compression, Chaquc unité locale de collecte comporte par exemple un ordinateur pourvu d'un processeur de signaux qui est programmé pour effectuer la compression des don nées D'autres caractéristiques et avantages de la méthode selon l'invention, apparaîtront à la lecture de la description ci-après d'un exemple non limitatif de réalisation, en sc refclant aux dessins annexés où - la Fig. 1 montre schématiquement un dispositif d'acquisition et de transmissioii de données sismiques - la Fig.2 montre schématiquement sous forme de blocs fonctionnels. la structure interne d'une unité locale d'acquisition - la Fig.3 illustre la technique connue de prédiction statistique de la valeur d'échantillons d'un signal - la Fig.4 schématise les opérations de compression d'un signal - la Fig.5 schématise les opérations de décompression symétrique d'un signal - la Fig.6 montre différentes ondelettes permettant d'analyser un signal à comprimer; - les Fig.7, 8. 9 montrent respectivement un signal sinusoïdal et ses transformées par ondelette de Haar et ondelette D20 - la Fig. 10 montre différentes fonctions d'ondelette et d'échelle bien adaptées à la compression de signaux sismiques - la Fig. 1 1 montre un histogramme d'erreur quadratique moyenne (RMS) obtcnu par exemple avec une ondelette S20 pour un signal sismique impulsionnel; et - la Fig. 12 montre un histogramme analogue pour un signal sismique vibratoire. Le dispositif sismique schématisé à la Fig. l comporte un ensemble souvent considérable (de plusieurs centaines à plusieurs milliers) de récepteurs sismiques R répartis à intervalles les uns des autres sur une zone à explorer. suivant une disposition convenant pour le type de prospection 2D ou 3D à effectuer, ces récepteurs captant les ondes sismiques renvoyées par des discontinuités soutel-raines, en réponse à la transmission dans le sol d'ondes sismiques produites par une source S, et une station centrale de commandc ct d'enregistrement 1 où tous les signaux sismiques collectés sont finalement centralisé par le rnoyen du système de transmission qui va être décrit. Chacun de ces récepteurs R est constitu le plus souvent d'une bretelle de capteurs élémentaires alignés qui produisent chacun une "trace sismique" Le dispositif comporte un ensemble d'unités locales A de collecte des données sismiques, chacune adaptée à acquérir au moins une trace sismique. L'ensemble des unités de collecte A communique avec une station centrale CS soit directement (par radio ou par câble), soit par l'intermédiaire cie stations intermédiaires LS avec des fonctions plus ou moins complexes. Il peut s'agir de concentrateurs destinés à organise et séquencer les échanges entre les unités locales RTU et la station centrale comme décrit par exemple dans le brevet EP-A-594 477 du demandeur. Outre ces fonctions de concentration, chaque station intermédiaire LS peut être dotée de fonctions de commande et contrôle de tâches diverses exécutées par les unités locales RTU, sous le contrôle de la station centrale comme décrit dans le brevet FR-A-2 720 518. Chacune de ces stations intermédiaires LS contrôle par exemple un nombre p d'unités locales de collecte A. Elles communiquent avec les unités de leurs groupes respectifs par une liaison radio ou éventuellement par une ligne de transmission 1. Les unités intermédiaires communiquent avec la station centrale SC par des canaux hertziens Fi, F2. ..., Fn. La collecte des données sismiques mémorisées dans les différentes unités d'acquisition peut encore s'effectuer par transfert direct sur une mémoire de masse déplacée sur le terrain successivement jusqu' au voisinage de chacune d'elles. Le dispositif comporte une source sismique S. Selon les cas. il s'agit d'une source impulsionnelle telle qu'une charge explosive par exemple produisant une trace sismique, ou d'un vibrateur. Cette source peut être couplée avec les terrains de la zone à explorer et reliée par radio ou câble de commande avec la station centrale ou bien dans le cas d'une exploration de zones côtières, éventuellement oemorquee en immension par un bateau boute-feu, relié par radio avec la station centrale. Chaque appareil d'acquisition Ai (Fig.2) est adapté à collecter par exemple k traces captés par un certain nombre k de récepteurs sismiques R1, R2. Rk. A cet effet, il comporte par exemple (k > 1) chaînes d'acquisition CAI à CAk recevant respectivement les k signaux et comportant chacune un filtre passe-bas F11. F12,. Fk un pré-amplificateur PAl, PA2 PAk, un filtre passe-haut F21. F27,.. F'k et un convertisseur analogique-numérique (ADC) Cl, C2 ..., Ck pour convertir les signaux analogiques amplifiés et filtrés en mots numériques. Toutes les chaînes sont connectées à un microprocesseur 2 traitant les mots numériques de 16 à 32 bits par exemple, programmé pour gérer l'acquisition et les échanges avec la station centrale 1. Au microprocesseur 2 sont associés deux blocs de mémoire Ml et M2 et une mémoire Mp pour les programmes. Le processeur 2 est connecté à une unité 3 d'émission-réception par radio ou par ligne, adaptée à la voie de transmission employée pour la communication avec la station centrale CS ou la station intermédiaire correspondante LS. S'il s'agit d'une voie hertzienne. l'unité 3 comporte un radio-émetteur RE et un radio-récepteur RR qui communique avec une antenne 4. Une unité d'interface 5 décrite dans le brevet FR-A-2.608.780 précité, permet en outre une communication par rayons infra-rouges avec un boîtier d'initialisation 6 à l'aide duquel un opérateur peut éventuellement communiquer au processeur de gestion 2, des instructions d'adressage et de sélection des paramètres de fonctionnement des chaînes d'acquisition. Le système selon l'invention de prétérellce, comporte dans chaque appareil d'acquisition Ai, un processeur 7 spécialisé dans le traitement des signaux. II peut s'agir par exemple d'un processeur à 32 bits à virgule flottante du type DSP 96002 par exemple, qui est associé à un dispositif du type DMA pour accélérer les traiisferts par blocs de données ente les deux processeurs 2 et 7. A ce dernier, est adjoint une mémoire de travail Mp. Chaque appareil d'acquisition comporte aussi une alimentation électrique autonome 8. Le processeur 2 agit en maître. II a pour fonctions de réaliser le décodage des ordres transmis par la station centrale 1. et de gérer - I'acquisition des signaux des récepteurs R1 à Rk par les différentes chaînes d'acquisition, - les transmissions en relation avec l'unité 3; - les mémoires M1 et M2 pour le stockage temporaire des données: - les entrées-sorties; - les interruptions entre programmes; - les échanges avec le processeur de calcul DSP 7 etc. De par sa structure propre, le processeur de calcul DSP 7 est particulièrement adapté à effectuer à grande vitesse des opérations telles que des conversions de format, des multiplications de nombres complexes, des transformations de Fourier du type FFT, des corrélations entre les signaux reçus et les signaux émis, des filtrages numériques, des sommations de tirs successifs avec élimination des bruits perturbateurs de nature non sismique, des combinaisons entre eux des signaux délivrés par des récepteurs sismiques multi-axes tels que des géophones tri-axiaux par exemple, etc. Les pré-traitements accomplis localement avant transmission contribuent à réduire sensiblement le nombre de tâches dévolues à la station centrale I et donc la puissance de calcul installée qui devient considérable quand le nombre de traces sismiques à acquérir en temps réel est de plusieurs centaines voire même dépasse mille. Comme on va le décrire ci-après, le processeur 7 de chaque unité de collecte A est également programmé pour appliquer aux traces collectées, un ou plusieurs algorithmes de compression de données de façon à réduire le volume de données à et à profiter au mieux des intervalles de temps plus ou moins longs intercalés entre les cycles d'émission-réception ou tirs d'une même série de cycles ou entre les séries de tirs successifs. A chacun des cycles successifs d'émission-oe'ception. une trace sismique est numérisée et enregistrée dans une mémoire de chaque unité de collecte Ai. De prétérence, on appliquc à la trace une compression de type sans perte. II peut S';lgil par exemple de la technique prédictive dite LPC précédemment évoquée où l'amplitude s(t) d'un échantillon de signal à un instant t (Fig.3) est calculée par une prédiction faite à partir d'un certain nombre p échantillons aux temps (t-l), (t-2), (t3), ... (t-p). La valeur prédite est calculée en accord avec la relation Au lieu de transmettre l'échantillon, on transmet sa prédiction sA(t) c'est-àdire les coefficients de prédiction et les résidus e(t) c'est-à-dire l'écart que l'on calcule entre la valeur réelle s(t) et la prédiction s^(t) qui en est faite au temps t, ce qui permet de retrouver la valeur s(t) = sA(t) + e(t) ) (Fig.4). Si la prédiction est bonne, les rési La méthode selon l'invention comporte l'application à chaque trace d'un taux de compression relativement important choisi de façon que cette trace comprimée puisse être transmise sur une voie de transmission relativcmcnt étroite et dans une fenêtre de transmission imposée par le déroulenient des sessions d'émission- réception. La durée d'un cycle d'acquisition est de l'ordre de quelqucs secondes (6s par exemple), supérieure le plus souvent à l'intervalle de temps séparait deux cycles successifs. On choisit donc un type de compression et un taux permettant de transmettre par exemple l'ensemble de la trace dans cette fenêtre imposée Ou au moins une partie significative de celle-ci. Il en résulte généralement une perte d'information qui est tolérable ici car les données transmises servent à des fins de contrôle. L'examen de la trace meme légèrement déformée par la compression, suffit à un opérateur dans la station CS (Fig. I ) pour vérifier que la trace a bien été enregistrée. Cette première transmission à des fins de contrôle est utilisée dans la méthode selon l'invention pour reconstituer à la station chacune des traces telle qu'elle a été enregistrée après compression sans perte dans les mémoires M 1, M2 (Fig.2) de chaque unité de collecte A. Plutôt que de transmettre sans perte la masse considérable de données correspondant à toutes les traces mémorisées au cours d'une même série ou d'une session quotidienne, on calcule dans chaque unité de collecte la différence entre d'une part la trace compressée sans perte (ou éventuellement la trace réelle non compressée), et d'autre part la trace compressée qui a été auparavant transmise à des fins de contrôle. Pour reconstituer chaque trace à la station centrale, on transmet dans un deuxième temps seulement cette différence. Cette procédure dc transfert en deux temps réduit globalement et de façon considérable, le ternps qui autrement serait nécessaire pour rapatrier les traces complètes. Comme technique de compression, on peut choisir aussi par exemple celle basée sur la transformée en ondelettes. Comme on l'a vu, l'analyse par ondelettes consiste à décomposer un signal quelconque sur une base de fonctions particulières d'un sous-espace ayant des propriétés bien déterminées. Une ondelette f ou y est une fonction dont les versions translatées (temporel) et dilatées (fréquence) forment une base de ce sous-espace. L'objectif désiré est de modéliser un signal en cherchant où peuvent être localisées (dans le temps) les irrégularités du signal ( par exemple des variations brusques ou hautes fréquences). On appelle f la fonction "d'échelle" qui permet d'avoir une version filtrée passe-bas du signal original, et y la fonction "ondelette" permettant d'avoir les détails (filtre passe-haut) entre deux versions consécutives filtrées par f. Dès lors, la réunion de ces informations nous permet de reconstituer le signal original suivant certaines conditions. On nomme "filtres d'analyse", les filtres permettant la transformée avant par onde lettes et "filtres de synthèse" ceux utilisés pour la tranformée inverse. Différentes familles d'ondelettes connues existent. On peut utiliser par exemple une famille de filtres appelés "filtres quadratiques miroir" (QMF) qui autorisent un filtre d'analyse identique à celui de synthèse, et un filtre représentant la fonction ondelette y miroir du filtre représentant la fonction d'échelle f (les coefficents sont utilisés en sens inverse). Par cette transformée en ondelettes, on obtient une liste de coefficients, et la compression intervient en annulant les coefficients les plus faibles en amplitude, inférieurs à un seuil déterminé. Le nombre de coefficients à conserver intacts est au choix de l'utilisateur en sachant que plus leur nombre est élevé, meilleure est l'approximation (la totalité des coefficients assure une reconstruction parfaite) lors de la transformée inverse. I1 faut noter que l'ordre des coefficients est primordial pour la reconstruction et que donc même les coefficients nuls sont à préserver. Ces coefficients une fois quantifiés en nombres entiers sont bien adaptés à une compression du type RLE déjà citée. On choisit de préférence une ondelette bien adaptée au filtrage de certaines fréquences du signal. Une ondelette bien adaptée aux irrégularités d'un signal comme I'ondelette de Haar (Fig.6) sur un signal n'ayant aucune irrégularité, une sinusolde par exemple (Fig.7), fournit des coefficients de détails élevés (Fig.8). Alors qu'une ondelette de type D20 (Fig.6) augmente en amplitude les coefficients cic basses Fréquences et minimise mieux les autres (Fig.9), ce qui est nettement mieux adapté à la compression. Finalelllent, pour ce type de transformation, on choisit Fig. 10) l'ondelette qui aura la capacité à concentrer un maximum d'énergie du signal dans les basses-fréquetices sera la meilleure d'un point de vue compression. Pour évaluer la qualité cte la reconstruction on détermine (Fig. 11, 12), l'erreur quadratique moyenne) dont le taux RMS qui mesure l'erreur dc reconstruction au sens de l'énergie du signal, est défini commc suit avec X .un échantillon du signal original, et x, un échantillon du signal reconstruit. Une méthode connue facilitant la sélection de la famille d'ondelettes la Inicux adaptée à la nature des signaux sismiques à traiter, est la Méthode de transformation par paquets d'ondelettes avec constitution d'une "librairie" de bases d'ondelettes, attribution d'une note à chaque base à l'aide d'une fonction de coût, telle que le logarithme de l'énergie ou l'entropie, note qui mesure l'importance de l'information contenue dans cette base; puis sélection des bases les mieux notées. perniettai0t la meilleure reconstitution du signal. La compression s'effectue également cn annulaiit les coefficients les plus faibles. On peut utiliser par exemple une fonction de coût d'un type connu définie par: Avec x" x,. l'ensemble des "j" coefficients d'unc base tonnée et # Ic coût de cette basc. Le nombre de coefficients nécessaire à la reconstruction ciu signal détermine le taux de compression appliqué. En l'occurrence, ce taux est imposé par le cadencement fixé des cycles d'missioll-réception sismique qui délimite les fenêtres de transmission des traces compressées transmises à des fins contrôle. On. n'envoie ainsi dans un premier temps que la partie la plus significative des coefficients d'ondelettes, Quant aux résidus, ils sont mémorisés dans les unités de collecte pour obtenir ultérieurement une reconstruction Dans un deuxième temps, on met à profit les intervalles de temps d'interruption inhérents à la technique de prospection utilisée, pour transmettre la différence entre la trace sismique complète et la partie déjà transmise à des fins de contrôle, en l'occurrence les résidus i.e. la série complète des coefficients de la transformation. Voici, à titre d'exemples, les taux d'erreur quadratique moyenne (RMS) obtenus en fonction du nombre de coefficients En prospection sismique avec source impulsionnelle: Nombre de coefficients Taux de compression RMS 100 97% 1,36% 500 91% 14e-3% En prospection sismique avec vibrateur Nombre de coefficients Taux de compression RNIS 100 83% 3,18% 500 57% 28e-2 Pour 100 coefficients d'ondelettes significatifs ( Ts = 97%, TV = 91% ). on obtient les taux RMS moyens suivants par WPT sur l'ensemble des traces D20 C30 S20 Explosif 1,53 % 1,38% % 1,36% Vibratoire 3,33% 3,24% 3,18N On vérifie que les ondelettes à tendance symétriques (C30 et S20), sont les mieux adaptées II s'avère que la sismique vibratoire en pratique est plus difficile à reconstruire que la sismique impulsionnelle. Voici par exemple à titre de comparaison, la même opération avec 500 coefficients (Ts = 8357o, Tv = 57% ): D20 C30 S20 Explosif 16e-3 % 14e-3% 14e-3% Vibratoire 29e-2% 28e-2% 28e-2% Les erreurs deviennent négligeables pour un taux de compression très intéressant. surtout en explosif. Modes de mise en oeuvre: Différentes possibilités peuvent être utilisées pour réaliser cette transmission en deux temps. 1) On peut combiner deux modes de compression différents, l'un pour compresser avec perte les données devant être transmises dans un premier temps à des fins de contrôle de qualité, l'autre pour compresser les données enregistrées dans chaque appareil d'acquisition en vue de leur transmission différée dans un deuxième temps. Dans ce cas, les données servant au contrôle qualité peuvent ne servir qu'à ccla. La transmission effectuée dans le deuxième temps portant sur la totalité des traces compressees. 2) On peut encore utiliser un seul mode de compression des données et leur appliquer deux taux de compression, le premier étant sans perte pour l'eiiregistoement local des données, le deuxième impliquant des pertes. étant appliqué aux données (I'ensemble de chaque trace par exemple) destinées à etre trallslllises durant le premier temps pour les besoins du controle de qualité. Dans ce cas. on stocke dans les mémoires de chaque unité de collecte A. les résidus traduisant les écarts entre les données complètes et celles transìnises avec perte dans le premier temps, lesquels résidus sont transmis dans un deuxième temps à la station centrale pour la reconstitution des signaux. 3) Dans le cas où l'on utilise une transformation par ondelettes permettant de traduire les données par une série de coefficients de valeurs décroissantes on peut choisir le taux de compression appliqué en sélectionnant le nombre limité dc coefficients que l'on transmet pour le contrôle de qualité, compatible avec la fenêtre de temps dont on dispose pour leur transmission entre deux tirs . Ce nombre est fixé en ne sélectionnant que les coefficients au-dessus d'une certaine valcur-seuil. Le choix peut être imposé par l'opérateur ou éventuellement laissé au système de traitement local. D'autres variantes sont possibles sans sortir du cadre de l'invention si elles sont compatibles avec le principe de la transmission en deux temps des données sismiques avec dans un premier temps une transmission de données compressées permettant un contrôle de la qualité de fonctionnement des unités locales et dans un deuxième temps une transmission de données permettant un enlegistr.ement sans perte des données sismiques à la station distante. On a décrit uti mode d'implémentation de la méthode au moyen d'un ensemble de traitement dans chaque appareil d'acquisition comportant un processeur de signal de type DSP. II est bien évident cependant que l'on peut utiliser à cet effet d'une façon plus générale tout dispositif de calcul adapté à effectuer toutes les opérations de compression requises. REVENDICATIONS I ) Méthode pour le transfert de données entre une ou plusieurs unités locales (A) de collecte de données sismiques, et une station distante (CS), par au moins une voie de transmission, caractérisée en ce qu'elle comporte une transmission en deux temps des données sismiques avec. dans un premier temps. une transmission de données compressées permettant un contrôle de la qualité de fonctionnement des unités locales et, dans un deuxième temps, une transmission de données permettant un enregistrement sans perte des données sismiques à la station distante. 2) Méthode selon la revendication 1, caractérisée en ce qu'elle comporte l'enregistrement local des données sans perte, une détermination de la différence entre les données enregistrées sans perte localement et les données correspondantes compressées émises durant le premier temps, une deuxième étape de transfert de cette différence, et une reconstitution à la station distante des données sismiques collectées par chaque unité locale de collecte, par combinaison des données reçues au cours des deux étapes de transmission et de transfert. 3) Méthode selon la revendication 1, caractérisée en ce que les données compressées comporte au moins une partie significative d'une trace sismique acquise localement par chaque unité de collecte. 4) Méthode selon la revendication 1 ou 2, caractérisée en ce que l'on choisit une technique de compression de données et/ou le taux de compression à appliquer aux données en fonction d'un intervalle de temps de transmission disponible et du débit de transmission permis par la voie de transmission utilisée. 5) Méthode selon l'une des revendications précédentes, caractérisée en ce que les données sismiques collectées sont enregistrées localement après une complession sensiblement sans perte et les données transmises au cours du premier temps de transmission, sont compressées avec un taux de compression impliquant une perte. 6) Méthode selon l'une des revendications précédentes, caractérisée en ce que les données sismiques étant collectées par des unités d'acquisition locales (A) au cours dc sessions d'acquisition de données comportant chacune une pluralité de cycles d'nregistrement de traces sismiques correspondant à des sigilaux sisiiiiqLies renvoyés par des discontinuités du sous-sol, en réponse à des signaux sismiques transmis dans le sol par une source sismique (S). les différents cycles d'une réelle session étant séparés les uns des autres par des intervalles de temps déterminés, lt méthode comporte une première compression de chaque trace sismique avec Liii preniier taux de compression suffisant pour que le premier temps de transmission puisse être intercalé durant un des dits intervalles de temps par UllC voie de transmission disponible, l'enregistrement local de chaque trace sismique dans une unité locale d'acquisition après application d'une compression sans perte avec un deuxième taux de compression. la détermination de la différeiice entre les traces sismiques compressées respectivement avec le premier et le deuxième taux, le deuxième temps de transmission différé étant durant Lin intervalle de temps entre les dites sessions, et la reconstitution à la station centrale de chaque trace sismique collectée. 7) Méthode selon l'unc des revendications précédentes, caractérisée en ce que la deuxième étape de transfert différé est effectuée par tille deuxiènie voie dc transltlissioll. X) Méthode selon l'une des revendications précédentes, caractérisée en ce que l'on applique aux données une technique de transformation par ondelettes po obtenir des données compressées. 9) Méthode selon la revendication 8, caractérisée en ce que l'on compresse les données utiles au contrôle de qualité, en sélectionnant des coefficients de transformation en ondelettes suprieurs à une valeur seuil déterminée. 10) Méthode selon l'une des revendications 1 à 9, caractérisée en ce que l'on applique aux données une techniquc de prédiction statistique poLir obtenir des données conipiessées. 11) Système pour le transfert de données entre une ou plusieurs unités locales (A) de collecte de données sismiques. et une station distante (CS), par au moins une voie de transmission, caractérisé en cc que chaque unité locale de collecte comporte des moyens (CA,M) d'acquisition et d'enregistrement de données, un ensemble de traitement (2, 7) comprenant des moyens de compression des donnces avec un taux dc compression choisi. des rnoyens pour déterminer la diffèretice entre les données enregistrées localement et les données correspondantes comprcssecs. des moyens (3) dc transmission des données compressées, des moyens de transfert de la dite différence à la stations, et un ensemble (10) de traitement dans la station distante (CS), comprenant des moyens pour reconstituer les données sisiiiiqties collectées par chaque unité locale dc collecte (A), par combinaison des données transmises et des données transférées. 12) Système selon la revendication précédente, caractérisé en ce qu'il comporte lensemble de traitement (2, 7) dans chaque unité locale de collecte est adapté à appliquer aux données sismiques une deuxième compression de type sans perte avec un deuxième taux cie compression, 13) Système selon l'une deux revendications Il ou 12, caractérisé en ce que chaque unité locale de collecte (A) comporte un ordinateur pourvu d'un processeur de signaux (7) et programmé pour effectuer la compression des traces sismiques.