FR 2664063 A1 19920103 FR 9008267 A 19900629 1 - L'invention a pour objet une méthode pour optimiser le déclenchement d'un ensemble de sources sismiques marines et un dispositif de contrôLe de tir pour sa mise en oeuvre. Les méthodes classiques de prospection sismique en mer sont conduites habituellement en utilisant un ensemble d'émission et un ensemble de réception d'ondes remorqués par un navire le long d'un plan de profil sismique à étudier Les ondes engendrées par l'ensemble d'émission sont réfléchies par les différents réflecteurs des formations immergées et sont reçues par l'ensemble de réception qui consiste généralement en une flûte sismique de grande longueur le long de laquelle sont répartis un grand nombre de capteurs. L'ensemble d'émission est constitué d'une source et plus communément d'une pluralité de sources remorquées en immersion et reliées au navire par des groupes de câbles multi-fonctions ou ombilicaux On utilise le plus souvent des sources impulsionnelles La forme d'onde produite dépend du type de source Si les sources sont du type à explosion comme les canons à air (ou airguns) par exemple, le pic principal est produit le premier Avec des sources du type à implosion, telles que des canons à eau (waterguns), le pic principal est précédé d'un pic précurseur de moindre amplitude. Ces sources sont immergées selon les cas à des profondeurs sensiblement égales ou bien systématiquement à des profondeurs différentes Les instants de déclenchement sont choisis avec précision compte-tenu de la géométrie particulière choisie pour l'ensemble d'émission (array) de manière à obtenir une source puissante et directionnelle Il s'agit, par un choix des instants de déclenchement des différentes sources, d'obtenir une mise en phase de leurs pics 2 - principaux respectifs dans une certaine direction L'opération est en général complexe car il faut tenir compte de multiples paramètres. Selon son type, sa profondeur d'utilisation et son état mécanique après le nombre de "tirs" qu'elle a effectués auparavant, l'instant effectif o se produit le pic principal de la source peut varier dans des proportions notables. Pour obtenir la mise en phase recherchée, on utilise généralement un appareil de commande appelé séquenceur de "tirs" adaptés à tenir compte des différents paramètres caractérisant l'ensemble d'émission utilisé Des capteurs sont placés à proximité des sources pour déterminer les instants de déclenchement et/ou la forme des impulsions produites. Des systèmes o l'on utilise des séquenceurs pour le contrôle de sources impulsionnelles sont décrites par exemple dans les brevets US No 4 599 712, 4 693 336, 4 718 045 4 739 858, 4 757 482 ou les demandes de brevet européen No 31 196 ou 48623. Le grand nombre de paramètres susceptibles d'influer sur le retard de déclenchement de chacune des sources de l'ensemble ou array, rend très complexe la mise en phase recherchée des pics principaux. D'une façon générale les dispositifs assurant le séquencement d'un ensemble d'émission multi-sources, ne tiennent compte que d'un nombre limité de paramètres, tels que le retard structurel au déclenchement, la profondeur d'immersion etc, et ne prennent pas bien en compte les fluctuations de ces paramètres dans les conditions réelles d'utilisation en mer. La méthode selon l'invention permet de mieux tenir compte des différents paramètres influant sur le retard effectif au déclenchement des différentes sources impulsionnelles d'un ensemble d'émission sismique, et d'éviter les causes d'imprécision résultant des systèmes antérieurs. Elle s'applique à un système de prospection sismique réalisant une succession de cycles d'émission-réception le long d'un profil à étudier au moyen d'un ensemble de sources impulsionnelles 3 - remorquées en immersion par un navire et reliées à un dispositif de contrôLe de tir adapté à gérer les séquences de déclenchement des sources et permet d'obtenir une mise en phase des impulsions émises par les différentes sources. La méthode selon l'invention est caractérisée en ce qu'elle comporte en combinaison, l'acquisition pour chacune des sources de l'ensemble d'émission, d'au moins un signal représentatif de l'impulsion produite par ladite source, la validation du signal acquis pour chaque source, par comparaison avec un signal de référence associé à la même source, extrait d'une table de données préalablement acquises portant sur l'ensemble des sources utilisées, de manière à sélectionner ceux qui présentent un degré de similitude suffisant, l'actualisation des données de la table par incorporation des signaux validés, pour tenir compte des modifications intervenues dans l'ensemble d'émission au cours des cycles d'émission-réception précédents, la détermination des instants de déclenchement optimaux des différentes sources en tenant compte des données actualisées, et le déclenchement des différentes sources auxdits instants de déclenchement optimaux. L'opération de validation comporte par exemple la comparaison des spectres de fréquence respectifs de chaque signal acquis et du signal de référence correspondant. L'opération de validation peut comporter également la comparaison par intervalle de temps, des énergies respectives de chaque signal acquis et du signal de référence correspondant. La méthode peut comporter en outre la détermination de l'instant de déclenchement effectif de chaque source, l'actualisation d'une seconde table de données représentatives des temps de réponse effectifs des différentes sources respectivement et l'ajustement des instants de déclenchement en fonction des données de la seconde table après actualisation. 4 - La détermination de l'instant de déclenchement de chaque source est effectuée notamment par une corrélation entre chaque signal validé et le signal de référence correspondant. Cette détermination de l'instant de déclenchement de chaque source peut encore être déterminée par comparaison de l'amplitude de chaque signal validé avec une valeur-seuil fixée. La méthode peut comporter aussi la mémorisation à intervalles de temps définis, correspondant à un certain nombre de cycles d'émission-réception, d'au moins une donnée représentative de l'impulsion émise par chaque source et la comparaison des données successives associées à chaque impulsion pour suivre son évolution au cours du temps. Le dispositif de contrôle de tir selon l'invention permet d'optimiser le déclenchement d'un ensemble d'émission sismique comportant une pluralité de sources impulsionnelles remorquées par un navire, de manière à obtenir une relation de phase optimale entre elles, chaque source étant associée à au moins un capteur de contrôle délivrant un signal représentatif de l'impulsion émise par la source à son déclenchement Il est caractérisé en ce qu'il comporte: un appareil d'acquisition connecté aux différents capteurs pour collecter les signaux qu'ils produisent à chaque tir les numériser et les mémoriser, un appareil de commande de tir pour commander le déclenchement des sources, cet appareil étant pourvu d'au moins une unité de mémorisation pour des données de référence obtenues au cours de séquences de tir antérieures; et un calculateur programmable pourvu de moyens de mémorisation pour au moins une table de données de référence, le calculateur communiquant avec l'appareil d'acquisition et l'appareil de commande de tir et déterminant les instants de déclenchement optimaux des différentes sources en tenant compte des données mémorisées, et pour transférer dans l'unité de mémorisation de l'appareil de commande avant le déclenchement des différentes sources, les valeurs effectives des retards au déclenchement qu'il faut respectivement leur appliquer. L'appareil d'acquisition comporte par exemple au moins un ensemble d'acquisition constitué d'une pluralité de chaînes d'amplification et de filtrage connectées par l'intermédiaire de moyens de multiplexage, des moyens de mémorisation pour les signaux numérisés issus de chaque ensemble d'acquisition, l'appareil d'acquisition et l'appareil de commande de tir communiquant avec le calculateur par un bus d'adresses et de données. Les moyens de mémorisation du calculateur comportent par exemple une unité de mémoire pour une première table de données de référence représentatives de signaux reçus respectivement par les capteurs associés aux différentes sources et validés, et pour une deuxième table de données de référence représentatives des instants de déclenchement associés aux signaux de référence validés. La référence qui est faite avant chaque tir à des données représentatives des conditions réelles des séquences d'émission précédentes, permet de garantir une adaptation constante aux éventuelles modifications des conditions de fonctionnement de l'ensemble d'émission. D'autres caractéristiques et avantages de la méthode et du dispositif selon l'invention apparaîtront mieux à la lecture de la description ci-après d'un mode de réalisation décrit à titre d'exemple non limitatif, en se référant aux dessins annexés o: la Fig 1 montre schématiquement un système d'émission-réception pour la prospection sismique en mer; la fig 2 montre la forme approximative d'une impulsion émise à l'émission par une source sismique du type explosive telle qu'un canon à air; la Fig 3 montre la forme approximative d'une impulsion correspondante émise par une source implosive du type canon à eau; la Fig 4 montre schématiquement Le dispositif de contrôle se Lon l'invention; la Fig 5 montre schématiquement l'appareil d'acquisition de signaux; la Fig 6 montre schématiquement une unité d'acquisition de signaux; la Fig 7 montre l'organigramme de chaque séquence préparatoire; et 6 - la Fig 8 montre l'organigramme d'analyse de tir incluant l'étape de vérification et de calcul menée par le calculateur. Un système d'émission-réception sismique marine, comme celui schématisé à la Fig 1, comporte un ensemble d'émission impulsionnel 1 remorqué en immersion par un navire 2, ainsi qu'un ensemble de réception constitué d'une flûte sismique 3 de grande longueur L'ensemble d'émission ou "array" 1 comporte généralement un nombre P de sources Si (i= 1,2 p) réparties en plusieurs sous- ensembles ("sub-arrays") constitués chacun de plusieurs sources réparties le long d'un faisceau de câbles multi-fonctions ou ombilicaux 4 Les sous-ensembles sont remorqués à la même profondeur avec des déports latéraux différents par rapport à la trajectoire du navire et/ou selon les cas, à des profondeurs différentes Chaque source Si est plus ou moins décalée. Par des lignes de transmission dans l'ombilical 4 qui la relie au navire, chaque source reçoit d'un système de contrôle 5 sur le navire un signal de déclenchement En retour, le système de commande 5 reçoit d'un capteur cinématique intégré à chaque source, un signal de déclenchement TB définissant l'instant exact de déclenchement Un capteur acoustique disposé à proximité immédiate de chaque source (à un mètre par exemple) ou selon les cas de chaque groupe de sources lorsque plusieurs d'entre elles sont regroupées, produit un signal de proximité NF que l'on désigne usuellement par "signature", qui est aussi transmis au système de contrôle 5 Ce même capteur est utilisé également entre les tirs pour mesurer la profondeur d'immersion de chaque source Les signaux TB sont en avance de quelques millisecondes par rapport aux signaux NF correspondants. Un instant initial ODT définit le début de chaque séquence de "tir" (Fig 2, 3) Par rapport à cet instant initial on définit les différents intervalles de temps significatifs suivants: Tsync: est le temps o les pics principaux des différentes sources de l'ensemble d'émission doivent être en phase dans une direction d'émission choisie E (Fig 1); TMA Xi: est le temps o le pic principal de la source Si est détecté 7- par le capteur proche associé; TRE Ti: est le temps de retard avant l'envoi à la source Si d'un ordre de déclenchement ou de mise à feu (O M F); TRE Si: est le temps de réponse de la source Si calcu Lé lors du "tir" précédent de l'ensemble d'émission. Ce temps est déterminé à la première arrivée d'énergie issue de la source IL dépend de la structure mécanique et/ou hydraulique de la source. ATXY Zi est la correction statique (Fig; 1) qui tient compte du décalage géométrique de la position de la source Si par rapport au centre du dispositif d'émission 1, compte-tenu de la direction d'émission E. TPS Pl est la pseudo-période (pour les sources à implosion) de la source Si La pseudo-période de chaque source implosive est calculée d'après la loi de Raleigh bien connue dans le domaine de la sismique, en fonction de la pression hydrostatique à la profondeur d'immersion réelle Zi de la source. Pour modifier l'instant du pic TMA Xi de chaque source Si, on doit jouer sur TRE Ti Les paramètres TRE Si, ATXY Zi et TPS Pl ne sont connus avec précision que par les mesures faites au cours de chaque cycle d'émission-réception. Après son déclenchement commandé avec le temps de retard TRET (instant tl) une source explosive (Fig 2) produit son pic principal (instant t 2) qui est suivi d'un pic secondaire de moindre amplitude (instant t 3). Le pic principal d'une source du type à implosion Si (Fig 3) est précédé à l'instant t 2 d'un pic précurseur L'écart (t 3-t 2) est la pseudo-période TPSP de la source. Le dispositif de contrôLe 5 selon l'invention permet d'adapter très précisément les intervalles de temps TRET et ATXYZ de chaque source aux conditions effectives de fonctionnement de l'ensemble d'émission, pour respecter les déphasages optimaux des différents pics principaux produits. Le dispositif de contrôle 5 comporte (Fig 4) un appareil 8- d'acquisition 6 recevant les signaux TB et/ou NF fournis par Les capteurs Ci associés à toutes les sources S de l'ensemble d'émission. Il comporte aussi un appareil de commande de tir 7 pour l'adaptation des signaux de déclenchement transmis aux électro-vannes des différentes sources Les appareils d'acquisition et de commande 6 et 7 communiquent par un bus d'adresses et de data 8 avec une carte d'interface 9, elle-même reliée à un calculateur programmable 10 Le dialogue entre les opérateurs et le calculateur 10 s'effectue par l'intermédiaire d'une console de commande 11 et d'un écran 12. L'appareil d'acquisition 6 (Fig 5) comporte un certain nombre de cartes d'acquisition A 1, A 2 An chacune pour un certain nombre p de voies (p= 32 par exemple) Chaque voie comporte (Fig 6) un préamplificateur 13 a 13 p en série avec un filtre passe-bande 14 a 14 p Le gain et les caractéristiques de chaque filtre 14 sont adaptées aux types de capteurs C 1 Cp employés Les signaux filtrés sont appliqués à un multiplexeur 15 à p entrées dont la sortie est connectée à un circuit de numérisation 16 Les signaux numérisés par chaque carte d'acquisition A 1 An sont stockés respectivement dans des cartes de mémoire M 1 à Mn Le calculateur peut venir les lire par l'intermédiaire de la carte d'interface 9 L'ensemble des cartes d'acquisition est piloté par l'intermédiaire d'une carte de synchronisation 17 connectée au bus 8 qui gère le multiplexage, la numérisation et le transfert des données numérisées vers la carte de mémoire M 1 à Mn correspondante. Le temps d'acquisition et la fréquence d'échantillonnage sont adaptés selon que les capteurs connectés aux entrées de chaque carte d'acquisition sont des signaux TB ou NF Les signaux TB sont acquis par exemple pendant une durée de 0,3 S à une fréquence d'échantillonnage de 10 K Hz La fréquence est de 5 Khz pour les signaux NF par exemple et leur durée d'échantillonnage est de 0,4 s. L'appareil de commande des canons comporte de préférence deux blocs de mémoire o le calculateur vient avant chaque séquence de déclenchement inscrire des mots numériques représentatifs des retards au déclenchement (TRET) qu'il a calculés et que l'appareil de commande 9 - va devoir appliquer respectivement aux différentes sources pour obtenir la mise en phase souhaitée au temps TSYNC Le calculateur 10 comporte des blocs de mémoire (non représentés) pour une table de référence Tl constituée de mots numériques représentatifs de La forme de L'impulsion émise par chaque source de L'ensemble d'émission 1. Dans ces blocs de mémoire, on inscrit aussi de préférence une deuxième table de référence T 2 représentative des temps de réponse effectifs (TRES) des différentes sources. Chaque cycle de fonctionnement du dispositif de contrôle comporte une étape d'émission de L'ensemble d'émission, une phase d'acquisition en temps réel géré directement par l'appareil d'acquisition et une phase de mesure et de contrôle Durant cette dernière phase le calculateur, à partir des données numérisées mémorisées au cours de la phase précédente, authentifie Les signaux TB ou NF reçus, par comparaison avec des données de référence, mesure leurs temps de réception et rafraîchit les données de tables de référence Tl et T 2 par inclusion des données authentifiées et s'en sert pour commander La séquence de déclenchement du cycle suivant, comme on va le décrire ci-après. Les données acquises à L'issue d'un cycle quelconque N, sont utilisées pour modifier L'instant OMF de chaque source (cf Fig 3), de manière à compenser les décalages éventuels notés dans La mise en phase des différentes sources. Lorsqu'un tir (N) a eu Lieu, Les cartes d'acquisition A 1 An de l'appareil 6 reçoivent et numérisent tous Les signaux TB et/ou NF captés près des différentes sources Si et Les mémorisent dans les mémoires Ml à Mn. Cette phase d'acquisition étant terminée, Le ca Lculateur 10 vient Lire les mémoires Ml à Mn d'une part et dans ses mémoires, la première table représentative de la forme des impulsions émises par les différentes sources Si. Dans une première étape, il fait une reconnaissance de forme par comparaison de chaque signal NF avec Les données correspondantes de La première table afin de vérifier que chaque - signal reçu est bien causé par la source De cette manière, on peut détecter les cas possibles o une source a mal fonctionné ou n'a pas fonctionné du tout bien qu'un signal d'origine parasite ait été reçu et mémorisé La similitude entre chaque signal et le signal de référence correspondant peut être effectuée par plusieurs méthodes Le calculateur 10 peut comparer les spectres de fréquence des signaux, par exemple ou bien encore comparer les valeurs efficaces ou les énergies des signaux par tranche de temps etc. Si le degré de similitude est supérieur à un seuil fixé, le signal est validé et inscrit dans la première table de référence T 1 mémorisée dans l'appareil de commande 7. Dans ce cas, le calculateur 10 détermine l'instant TMA Xi o a eu lieu le pic principal de chaque source Si Cette opération est faite par corrélation entre le signal de référence et le signal TB ou NF de chaque source par exemple ou bien encore par comparaison avec une amplitude seuil à l'intérieur d'une fenêtre de temps bien déterminée L'instant de tir ainsi calculé est inscrit dans la deuxième table T 2 formée dans la mémoire du calculateur 10. Le calculateur 10 calcule alors les modifications à apporter éventuellement aux temps TRET pour le tir suivant de façon que les pics principaux soient bien synchrones avec le temps TSYNC choisi Le temps TRESN obtenu au cours du cycle N pour chaque source Si, est calculé par la relation TRES = TMAX TRET TPSP 1) L'intervalle TRETN+ 1 pour le tir suivant N+ 1 est obtenu par la relation: TRETN+ 1 = TSYNC TRESN ATXYZN+ 1 TPSPN+ 1 2) N+ 1 N N+ 1 N+ 1 C 2) o ATXYZN+ 1 représente la correction géométrique liée au changement éventuel de la configuration de l'ensemble d'émission qui a pu intervenir entre deux tirs successifs et TPSPN+ 1 la nouvelle pseudo-période de la source consécutive à une variation éventuelle de sa profondeur d'immersion entre deux tirs successifs. Le calculateur va modifier la séquence de tir de manière à tenir compte des nouvelles valeurs calculées pour les différents TRET. il - Les coordonnées réelles Xi et Yi de chaque source Si peuvent être mesurées régulièrement en cours de remorquage et intégrées à chaque fois dans le calcul de chaque TXY Zi On peut aussi Les considérer comme suffisamment stables au cours de chaque séance d'enregistrement sismique et ne tenir compte avant chaque tir que des variations de profondeur éventuelles mesurées par le capteur fournissant le signal TB. Les données inscrites dans les deux tables T 1, T 2 peuvent venir remplacer les données précédentes ou bien encore servir à corriger une moyenne faite sur un certain nombre de tirs précédents. Ainsi, les données de référence sont constamment actualisées en tenant compte des conditions réelles de fonctionnement de l'ensemble d'émission. Au début de chaque série de cycles d'émission-réception sismique, en l'absence de signaux de référence, on peut inclure Les signaux captés dans les Listes de signaux de référence sans analyse préalable. Si la similitude des signaux acquis avec le signal de référence correspondant est insuffisante, le signal est rejeté. L'opérateur peut alors commander la visualisation de ce signal et si sa non-validation est due à une défaillance de la source, il peut commander son inhibition pour les séquences de tir suivantes. Le calculateur 10 peut encore comporter des mémoires auxiliaires pour La mémorisation après validation, des différents signaux produits par les sources à intervalles réguliers de m tirs (tous les 100 ou 1000 tirs par exemple) L'opérateur peut ainsi commander au calculateur de faire des comparaisons des signaux d'une même source et l'on peut suivre l'évolution de son fonctionnement au cours du temps. Le dispositif de contrôle peut encore commander des tirs unitaires, o chaque source est déclenchée seule de façon à déterminer séparément les signaux associés TB ou NF et faire des étalonnages. Les deux organigrammes des Fig 7 et 8 résument Les opérations de vérification et de calcul menées par le calculateur 10. 12 - Sur l'organigramme de La Fig 7 relatif aux opérations précédant chaque tir et l'organigramme général de la Fig 8 incluant l'étape d'analyse postérieure à chaque tir, les abréviations utilisées ont les significations suivantes: FLACT: source donnée dans l'ensemble d'émission 0: inactive 1: active FLTBR: capteur TB étalonnage 0: non fait 1: demandé 2: contrôle 3: fait FLXYZ: calcul de position 0: non positionné 1: demandé 2: positionné FLSTA: mode de fonctionnament 0: manuel 1: automatique demandé 2: automatique FLTIR: tir 0: inhibé 1: autorisé 2: inhibé; tir unitaire 3: autorisé; tir unitaire FLTYP: type de canon 1: implosif 2: explosif FLFFS: constitution du fichier catalogue de F F synthétique 1: catalogue 2: participation au calcul FFS TRETD, TRETM, TRETE: retard de tir donné, mesuré, écart F.F synthétique: étape de détermination de la signature du dispositif d'émission à grande distance dite "far field" 13 On ne sortirait pas de l'invention en constituant d'autres tables de façon à mémoriser L'évolution au cours du temps d'autres paramètres tel que par exemple le ATXYZ pour chaque source de l'ensemble d'émission. 14 - REVENDICATIONS 1) Méthode pour optimiser le déclenchement d'un ensemble de sources sismiques marines impulsionnelles (Si) utilisées pour des opérations de prospection sismique marine comportant une succession de cycles d'émission-réception d'impulsions sismiques, l'ensemble de sources étant remorqué en immersion par un navire et relié à un dispositif de contrôLe de tir ( 5) adapté à gérer les temps de déclenchement des sources de façon à obtenir une mise en phase des impulsions émises par les différentes sources, caractérisée en ce qu'elle comporte en combinaison, l'acquisition pour chacune des sources de l'ensemble d'émission, d'au moins un signal (TB, NF) représentatif de l'impulsion qu'elle a produite, la validation du signal acquis pour chaque source, par comparaison avec un signal de référence associé à la même source extrait d'une table de données T 1 préalablement acquises portant sur l'ensemble des sources utilisées, de manière à sélectionner ceux qui présentent un degré de similitude suffisant, et l'actualisation des données de la table par incorporation des signaux validés, pour tenir compte des modifications intervenues dans l'ensemble d'émission au cours des cycles d'émission-réception précédents, la détermination des instants de déclenchement optimaux des différentes sources en tenant compte des données actualisée, et le déclenchement des différentes sources auxdits instants de déclenchement optimaux. 2) Méthode selon la revendication 1, caractérisée en ce que la validation comporte la comparaison des spectres de fréquence respectifs de chaque signal acquis et du signal de référence correspondant. 3) Méthode selon la revendication 1, caractérisé en ce que la validation comporte la comparaison par intervalle de temps, des énergies respectives de chaque signal acquis et du signal de référence - correspondant. 4) Méthode selon l'une des revendications précédentes, caractérisée en ce qu'elle comporte la détermination de l'instant de déclenchement effectif de chaque source (TMAX), l'actualisation d'une seconde table (T 2) de données représentatives des temps de réponse effectifs (TRES) des différentes sources respectivement et l'ajustement des instants de déclenchement en fonction des données de la seconde table après actualisation. ) Méthode selon la revendication 4, caractérisée en ce que la détermination de l'instant de déclenchement de chaque source est effectuée par une corrélation entre chaque signal validé et le signal de référence correspondant. 6) Méthode selon la revendication 4, caractérisé en ce que la détermination de l'instant de déclenchement de chaque source est déterminée par comparaison de l'amplitude de chaque signal validé avec une valeur-seuil déterminée. 7) Méthode selon l'une des revendications précédentes, caractérisée en ce qu'elle comporte en outre la mesure des variations de la géométrie du dispositif d'émission et de sa position, la détermination des instants optimaux de déclenchement tenant compte de lesdites variations de géométrie et de position. 8) Méthode selon l'une des revendications précédentes, caractérisée en ce qu'elle comporte en outre la mémorisation à intervalles de temps définis, correspondant à un certain nombre de cycles d'émission-réception, d'au moins une donnée représentative de l'impulsion émise par chaque source et la comparaison des données successives associées à chaque impulsion pour suivre son évolution au cours du temps. 9) Dispositif de contrôle de tir pour optimiser le déclenchement d'un ensemble d'émission sismique comportant une pluralité de sources impulsionnelles (Si) remorquées par un navire ( 2), de manière à obtenir une relation de phase optimale entre elles, chaque source étant associée à au moins un capteur de contrôle (Ci) délivrant un signal (TB, NF) représentatif de l'impulsion émise par la 16 - source à son déclenchement, le dispositif de contrôle étant caractérisé en ce qu'il comporte: un appareil d'acquisition ( 6) connecté aux différents capteurs pour collecter les signaux qu'ils produisent à chaque tir, les numériser et les mémoriser, un appareil de commande de tir ( 7) pour commander le déclenchement des sources, cet appareil étant pourvu d'au moins une unité de mémorisation des données de référence obtenues au cours de séquences de tir antérieures; et un calculateur programmable ( 10) pourvu de moyens de mémorisation pour au moins une table de données de référence, le calculateur communiquant avec l'appareil d'acquisition ( 6) et l'appareil de commande de tir ( 7) et déterminant les instants de déclenchement des différentes sources en tenant compte des données mémorisées, et pour transférer dans l'unité de mémorisation de l'appareil de commande avant le déclenchement des différentes sources les valeurs effectives des retards au déclenchement (TRET) qu'il faut respectivement leur appliquer. ) Dispositif selon la revendication 9, caractérisé en ce que l'appareil d'acquisition comporte au moins un ensemble d'acquisition (Ai) constitué d'une pluralité de chaînes d'amplification et de filtrage ( 13 i,14 i) connectées par l'intermédiaire de moyens de multiplexage ( 15), des moyens de mémorisation (Mi) pour les signaux numérisés issus de chaque ensemble d'acquisition (Ai), l'appareil d'acquisition ( 6) et l'appareil de commande de tir ( 7) communiquant avec le calculateur ( 10) par un bus d'adresses et de données ( 8). 11) Dispositif selon la revendication 9 ou 10, caractérisé en ce que les moyens de mémorisation du calculateur ( 10) comportent une unité de mémoire pour une première table de données de référence représentatives de signaux reçus respectivement par les capteurs associés aux différentes sources et validés et pour une deuxième table de données de référence représentatives des instants de déclenchement associés aux signaux de référence validés.