Brevet FR2781944A1 — publié le 20000204

La présente invention concerne une méthode de transmission de données utilisant des techniques de compression sans pertes, en vue d'optimiser l'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 0o 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 centrale 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 FR 2 720 518, 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, l'une à débit de transmission relativement élevé, l'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 et 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 de transmission à bande passante relativement étroite. Par les brevets FR-A-2 692 384, 2 757 641 et la demande de brevet FR97/09.547 du demandeur également, il est connu d'utiliser des boîtiers d'acquisition sismique pourvus notamment de moyens de traitement 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, des pré-traitements des traces sismiques auparavant réalisés à la station centrale après transmission, adaptés également à appliquer aux données sismiques à transmettre des algorithmes de compression, tous pré-traitements visant à diminuer considérablement le volume de données à rapatrier à la station centrale. La tendance actuelle, notamment dans le cadre des méthodes d'exploration sismique dites 3D, est à répartir sur une zone à explorer, à terre, en mer, 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 à transmettre ne cesse de croître. Pour éviter que les problèmes de transmission ne constituent un frein à l'évolution des systèmes sismiques la tendance est de recourir à des procédés de compression de données choisis pour être compatibles avec les exigences propres des géophysiciens. La compression des données sismiques peut apporter à la fois un gain de place appréciable sur les modules de stockage de masse dans les boîtiers locaux d'acquisition et/ou les stations locales de contrôle et de concentration, et aussi un gain considérable en temps de transmission. 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. En géophysique notamment, il est indispensable 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érique effectué sur plusieurs traces. Des pertes éventuelles de précision ne sont tolérables que dans certains cas bien particuliers, si par exemple l'information transmise ne doit servir qu'à contrôler le bon fonctionnement du matériel et à visualiser '"l'allure" des traces échantillonnées. 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, les méthodes de type RLE pour " constant de bits. Une technique connue de compression dite "LPC pour faite à partir de p échantillons précédents, en supposant que le signal est stationnaire. Au lieu de transmettre l'échantillon s(t), on transmet sa prédiction. ?(t) c'est-à-dire les coefficients de prédiction ainsi que les résidus e(t) i.e. 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) = t(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). Le nombre de coefficients servant à calculer,(t) est en général faible et ils occupent peu de place par rapport à s(t). Cependant, avec ce type de codage, il est conventionnel de ranger les échantillons de traces successifs dans des mémoires-tampons (" buffers > ") de taille définie et de calculer des coefficients de prédiction qui sont utilisés pour tous les échantillons d'un même groupe. Il s'avère dans la pratique que ce type de convention est mal adapté aux signaux à grande dynamique comme peuvent l'être les signaux sismiques émis par des sources impulsionnelles (charges de dynamite par exemple), et que des erreurs de codage apparaissent du fait que les coefficients prédictifs choisis ne sont pas bien représentatifs. La méthode selon l'invention permet de réaliser une compression sans pertes de signaux à grande dynamique, tels que des signaux sismiques, en vue de leur transmission. La méthode comporte l'échantillonnage des signaux et leur numérisation, et l'utilisation d'une technique de codage prédictif. Elle est caractérisée en ce qu'elle comporte la détermination d'une plage d'amplitude dépendant d'au moins une valeur seuil (S), qui contient l'essentiel des résidus de prédiction, le codage des résidus de prédiction par des mots binaires ayant un certain nombre de bits suffisant pour le codage de cette plage et le codage des résidus de prédiction dont l'amplitude est en dehors de la plage, par des mots binaires ayant un nombre fixe de bits supérieur à n, et la transmission des seuls résidus de prédiction (les coefficients de prédiction ne sont pas transmis mais recalculés). Suivant un mode de mise en oeuvre, la méthode peut comporter la détermination d'un facteur d'estimation (S') de la valeur-seuil optimale correspondant à une moyenne faite à partir de la position des bits de poids forts respectifs de N résidus de prédiction successivement déterminés. Suivant un autre mode de mise en oeuvre, la méthode peut comporter aussi la détermination de plusieurs plages d'amplitude dépendant de plusieurs valeurs seuils distinctes. La méthode peut aussi comporter par exemple un codage séparé des bits de poids faibles et forts des échantillons de signaux. Suivant un mode préféré de mise en ceuvre, la méthode comporte une détermination des résidus de prédiction pour chaque nouvel échantillon de signal. La méthode selon l'invention est avantageuse dans ses deux étapes prises séparément ou de préférence en combinaison. a) la détermination que l'on effectue du nombre de bits au plus nécessaire au codage de la plus grande partie des erreurs, permet de d'obtenir un taux de compression très important même si l'on tient compte du traitement particulier que l'on est amené à effectuer pour coder les erreurs atypiques restantes (peu nombreuses quand la plage de codage est bien choisie); b) le codage est ici parfaitement adaptatif puisque les coefficients prédictifs sont calculés à chaque échantillon, ce qui garantit une stabilité des résultats; 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 se référant aux dessins annexés o: - la Fig. l illustre la technique connue de prédiction statistique de la valeur d'échantillons d'un signal; - la Fig.2 schématise les opérations de compression d'un signal - la Fig.3 schématise les opérations de décompression symétrique d'un signal; et - la Fig.4 montre dans le cas d'un codage à un seuil, une plage S d'amplitude choisie pour contenir l'essentiel des résidus de prédiction. DESCRIPTION DETAILLEE 1) Codage adaptatif Le procédé selon l'invention comporte l'application aux données à transmettre d'une compression de type sans perte en utilisant la technique prédictive dite LPC o l'amplitude s(t) d'un échantillon de signal à un instant t (Fig.1) est calculée par une prédiction faite à partir d'un certain nombre p échantillons aux temps (t-l), (t-2), (t-3),... (t-p), et la valeur prédite est calculée en accord avec la relation suivante: -1' Au lieu de transmettre l'échantillon, on calcule en ligne les coefficients de prédiction et on transmet 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, cc qui permet de retrouver la valeur s(t) = ?;(t) + e(t) (Fig.2). Au temps t, les valeurs de s aux temps inférieurs ont déjà été calculées, on connaît de ce fait s(t-j) pourj j [1, p]. Les valeurs s(t-j) de l'équation (1) qui ne sont pas connues pour t E [0, p-l] sont remplacées par la valeur 0. L'écart entre le modèle calculé et le signal réel (appelé par la suite erreur de prédiction ou résidu) est codé de façon à réduire la taille moyenne allouée à chaque échantillon (par rapport aux 24 bits initiaux), et à réaliser une compression. Si la prédiction est bonne, les résidus sont faibles et peuvent être codés avec des mots binaires avant un nombre de bits, du moins en moyenne, bien inférieur généralement à celui servant à coder chaque échantillon s(t) (24 bits par exemple). A l'étape de décompression, (Fig.3), le signal est restitué selon la formule suivante: s(t) =,_aj.s(t- j) +e(t); avec t E [0,T] (2) i=l Les résidus compressés e(i) du fichier initial sont décompressés et l'on recalcule les coefficients de prédiction a(k) pour chaque échantillon. De préférence, on réalise une prédiction adaptative: les coefficients a,, a,....,apde prédiction du filtre d'évaluation sont recalculés pour chaque nouvel échantillon, afin de tenir compte de la dynamique souvent très grande des signaux sismiques. Le procédé selon l'invention offre la possibilité de réduire de façon substantielle la taille des mots binaires de codage de ces résidus. On va décrire ci- après une méthode de codage à un seul seuil, permettant de définir une plage d'amplitude contenant l'essentiel des résidus de prédiction. II) Codage des résidus: méthode à un seuil Si l'on suppose que les erreurs peuvent être codées en valeur absolue sur Nma, bits et que Nmax,, correspondant à -2s (1 0... 0). ze, as j fichier comprimé débutera par la séquence suivante: 00000111 0(0000100 10000 00010010 10111 01110 N-,",s.S code en et Sans modélisation simple de la distribution de l'erreur, il est difficile de déterminer a priori la valeur optimale de S, et les performances peuvent être considérablement dégradées par un mauvais choix de cette dernière. Comme on ne peut envisager un balayage de toutes les valeurs possibles pour S, on cherche un bon estimateur de cette valeur. L'estimateur S retenu est défini de la façon suivante: on appelle Nb(e) le nombre de bits qu'occupe un entier positif e (hors signe), ou encore la position de son bit de poids fort. Par convention on considère que Nb(O)= 1. f six =0 Nb(e) =.Llog-(e)j + I.si x > 1, (L J désigne l'entier immédiatement inférieur) On choisit alors: 2 + arrondi( Nb(E )jj, dont le terme final est la moyenne des tailles (en bit) des erreurs de prédiction (en valeur absolue). Cette valeur est partiellement heuristique. Elle a donné le paramètre optimal dans 95 % des cas testés, et sinon une valeur différant d'une unité. On calcule cet estimateur "a la volée", au fur et à mesure de l'arrivée des erreurs. Avantages Cet algorithme de codage présente les avantages suivants: performances en compression: les tests effectués sur les données terrestres (dynamite) produisent, à ordres de prédiction égaux, des résultats meilleurs que ceux des autres méthodes testées, et surtout plus stables, avec des taux de compression situés entre 50 et 70 % (en moyenne aux alentours de 60 %). Pour des données acquises en haute mer ( pas nécessaire de pratiquer à cette étape de coûteuses itérations d'optimisation. Résultats Pour une dizaine de " tirs " donnant lieu chacun à une centaine traces chacun durant deux campagnes d'acquisition sismique, réalisés avec des explosifs (en terrestre) et des canons à air (en marine), chaque trace comprenant chacune 3000 à 4000 échantillons, on a obtenu par exemple les résultats suivants (ordre de prédiction 8): Taux moyen Min./max Ecart-type Terrestre 61 49/71 6,5 Marine 75 65/90 8 On a décrit une méthode de codage à un seul seuil S. Bien entendu, on ne sortirait pas du cadre de l'invention en choisissant un codage à plusieurs valeurs seuils distinctes codables respectivement par plusieurs mots de tailles fixes avec évaluation de la taille du mot suffisante pour coder à chaque instant les résidus de prédiction. Pour améliorer les résultats, il est également possible par exemple a) d'augmenter l'ordre de prédiction (jusqu'à 15 par exemple); b) d'affiner la méthode de codage de l'erreur cn codant séparément les bits de poids faibles (de 1 à 4) et les bits de poids forts du signal. REVENDICATIONS 1) Méthode de compression sans pertes de signaux à grande dynamique, tels que des signaux sismiques, en vue de leur transmission, comportant l'échantillonnage des signaux et leur numérisation, et l'utilisation d'une technique de codage prédictif avec transmission de coefficients de prédiction et de résidus de prédiction, caractérisée en ce qu'elle comporte la détermination d'une plage d'amplitude dépendant d'au moins une valeur seuil (S), qui contient l'essentiel des résidus de prédiction, le codage des résidus de prédiction dans cette plage lé d'amplitude par des mots binaires ayant un certain nombre de bits suffisant pour le codage de cette plage, le codage des résidus de prédiction dont l'amplitude est en dehors de la plage, par des mots binaires ayant un nombre fixe de bits supérieur à n, et la transmission des seuls résidus de prédiction. 2) Méthode selon la revendication l, caractérisée en ce qu'elle comporte la u détermination d'un facteur d'estimation (S) du correspondant à une moyenne faite à partir de la position des bits de poids forts respectifs d'une suite de N résidus de prédiction successivement déterminés. 3) Méthode selon l'une des revendications précédentes, caractérisée en ce qu'elle comporte la détermination de plusieurs plages d'amplitude dépendant de plusieurs valeurs seuils distinctes. 4) Méthode selon l'une des revendications précédentes, caractérisée en ce qu'elle comporte un codage séparé des bits de poids faibles et des bits de poids forts des échantillons de signaux. ) Méthode selon l'une des revendications précédentes, caractérisée en ce qu'elle comporte une détermination des résidus de prédiction pour chaque nouvel échantillon de signal.