Brevet FR3077645A1 — publié le 20190809

Description Titre de l'invention : SONAR DOPPLER DE POURSUITE À PORTÉE ÉTENDUE Contexte [00011 Domaine d'application [0002] La technologie décrite concerne les systèmes de mesure acoustique sous-marine et, plus particulièrement, un sonar Doppler qui étend la portée de poursuite du fond. [0003] Description de la technologie connexe [0004] Les sonars Doppler, y compris les enregistreurs de vitesse Doppler (DVL) et les profileurs de courant Doppler acoustiques (ADCP), mesurent la vitesse relative entre l'instrument et un groupe de diffuseurs en émettant des impulsions acoustiques le long de plusieurs faisceaux qui pointent dans différentes directions et mesurent le décalage Doppler du signal acoustique qui est rétrodiffusé vers l'instrument dans chaque faisceau. Le groupe de diffuseurs peut être constitué soit de particules en suspension dans la colonne d'eau pour mesurer les courants, soit d'une surface limite, par exemple le fond marin, pour mesurer la vitesse par rapport au sol. Une configuration commune 100 se compose de quatre faisceaux 120a, 120b, 120c et 120d provenant d'un instrument 110 séparés en azimut de 900 et surélevés par rapport à la verticale 140 d'un angle commun désigné l'angle Janus aj 130, tel qu'illustré à la [fig 11. D'autres configurations existent, comme une configuration à trois faisceaux. [0005] La fréquence Doppler radiale mesurée par un faisceau est approchée par [Math.1] [0006] où vr est la vitesse radiale, f, est la fréquence d'émission et c est la vitesse du son. [0007] Les sonars Doppler à bande étroite et à large bande sont deux types courants de sonars Doppler qui se distinguent par les séquences d'émission qu'ils utilisent. Les sonars Doppler à bande étroite émettent une sinusoïde de durée finie. Une représentation mathématique est donnée par [Math.21 s (t) A cos (oct. t Ttx 0 els e [0008] où A est l'amplitude, co, est la fréquence porteuse angulaire, t est le temps et Tt, est la durée d'émission. Les sonars Doppler à large bande, d'autre part, émettent une 2 séquence de codes qui sont répétés Np fois pour former la séquence d'émission complète. Le code se compose généralement d'un signal à compression d'impulsions ou de symboles codés en phase. Un exemple de ces derniers est donné, pour lequel un seul code est composé de L symboles de durée T. La représentation en bande de base complexe est sbb(t) et le signal équivalent passe-bande est [Math.31 (t) Re sbb(t)e eect [0009] La représentation en bande de base de la séquence d'émission complète peut être écrite comme la somme des impulsions codées N. [Math.41 Shb (t) rk% Sc(t ki,Ts) [0010] où s(t) est la représentation en temps continu d'un code à temps discret x,[n]. Un filtre de conformation d'impulsions peut être utilisé pour générer la représentation en temps continu du code. Par exemple, un filtre de conformation d'impulsions rectangulaires génère s(t) à partir de x,[n] selon [Math.51 c[rtirec [0011] Un exemple de code de séquence de longueur maximale de longueur sept est x,[n] = [1, 1, 1, -1, -1, 1, -11, qui est un code couramment utilisé dans les DVL à réseau phasé. Toutefois, il est possible d'utiliser tout type de code de toute longueur et toute autre forme d'onde. [0012] La portée de poursuite de fond maximale réalisable d'un DVL est déterminée par l'équation sonar et le seuil du rapport signal/bruit (RSB) de l'algorithme de traitement. L'équation sonar prédit le rapport entre la puissance moyenne du signal et la puissance moyenne du bruit, à une portée donnée, à partir du niveau de puissance émis, de la géométrie et du rendement du transducteur, de la fréquence acoustique, des conditions environnementales et du niveau de bruit. L'approche évidente pour obtenir une portée maximale longue est d'augmenter la puissance du signal. L'augmentation de la puissance d'émission est coûteuse du point de vue de la consommation d'énergie et il est peu pratique d'obtenir un gain de RSB suffisant pour obtenir une augmentation significative de la portée maximale. De plus, la cavitation et les chocs limitent le RSB maximal que l'on peut obtenir en augmentant la puissance d'émission. Pour une géométrie et un rendement de transducteur donnés, une fréquence acoustique et des 3 conditions environnementales données, les stratégies restantes incluent la réduction de la puissance de bruit et l'amélioration de l'algorithme de traitement pour fonctionner à un RSB inférieur. [0013] La portée des DVL à bande étroite est généralement supérieure à celle des DVL à large bande, en raison de la largeur de bande plus étroite qui réduit la puissance de bruit totale. Cependant, la portée supplémentaire se paye par un écart-type accru et une fidélité moindre de la détection du fond - deux effets secondaires du produit temps/ largeur de bande unitaire. Le développement décrit combine la capacité de portée étendue de la bande étroite avec l'écart-type plus faible de la poursuite du fond à large bande. [0014] Le brevet américain n° RE 35,535 décrit un profileur de courant Doppler acoustique à large bande (ADCP). Le brevet américain n° 5,315,562 et le brevet américain n° 5,422,860 décrivent chacun des aspects d'enregistreurs de vitesse à corrélation. Le brevet américain n° 5,808,967 décrit des réseaux phasés. Le brevet américain n° 7,542,374 décrit pour les réseaux phasés un procédé permettant d'éliminer sensiblement un biais lié à la composante de vitesse orthogonale à la surface. Le brevet américain n° 7,839,720 décrit l'utilisation du codage pour éliminer le biais de couplage des lobes latéraux dans les systèmes à réseaux phasés. Le brevet américain n° 5,923,617 décrit des réseaux blazés. Résumé de CERTAINS ASPECTS INVENTIFS [0015] Dans un aspect, il existe un système de sonar actif sous-marin pour mesurer la vitesse d'un instrument par rapport à une surface limite, le système comprenant un transducteur acoustique configuré pour émettre et recevoir une pluralité de faisceaux acoustiques dans différentes directions ; et un processeur configuré pour : détecter une surface limite à l'intérieur de chaque faisceau émis ; filtrer de manière itérative les signaux acoustiques reçus rétrodiffusés à partir des faisceaux émis avec un filtre adaptatif et la bande passante associée qui est diminuée successivement pour chaque itération ; et mesurer la vitesse d'un instrument par rapport à la surface limite. [0016] Le processeur peut être configuré en outre pour détecter une surface limite à l'intérieur de chaque faisceau en trouvant une réponse de crête dans un signal de coefficient de corrélation. Le système peut en outre comprendre un module émetteur configuré pour émettre une séquence codée en phase comprenant des codes répétitifs de signe alternatif. Le processeur peut être configuré en outre pour former des faisceaux de manière à séparer les signaux acoustiques reçus. Le processeur peut être configuré en outre pour mesurer la vitesse du véhicule par rétrodiffusion de son sur la surface de fond d'une masse d'eau. Le transducteur acoustique peut émettre une ou plusieurs répétitions d'un signal codé en phase ou à compression d'impulsions pour produire un signal à large bande. Le transducteur acoustique peut comprendre au 4 moins un parmi un réseau phasé, un réseau de réseaux phasés, un réseau multicanaux, un réseau blazé, un réseau de réseaux blazés et un ensemble de transducteurs à piston. [0017] Dans un autre aspect, il existe un système de sonar actif sous-marin pour mesurer la vitesse d'un instrument par rapport à une surface limite, le système comprenant un transducteur acoustique configuré pour émettre et recevoir une pluralité de faisceaux acoustiques dans différentes directions ; et un processeur configuré pour : détecter une surface limite à l'intérieur de chaque faisceau émis en trouvant une réponse de crête dans un signal de coefficient de corrélation ; et mesurer la vitesse d'instrument par rapport à la surface limite. [0018] Le processeur peut être configuré en outre pour filtrer de manière itérative les signaux acoustiques reçus rétrodiffusés à partir des faisceaux émis avec un filtre et la largeur de bande associée qui peut être diminuée successivement pour chaque itération. Le système peut en plus comprendre un module émetteur configuré pour émettre une séquence codée en phase comprenant des codes répétitifs de signe alternatif. Le processeur peut être configuré en outre pour former des faisceaux de manière à séparer les signaux acoustiques reçus. Le processeur peut être configuré en outre pour mesurer la vitesse du véhicule par rétrodiffusion de son sur la surface de fond d'une masse d'eau. Le transducteur acoustique peut émettre une ou plusieurs répétitions d'un signal codé en phase ou à compression d'impulsions pour produire un signal à large bande. Le transducteur acoustique peut comprendre au moins un parmi un réseau phasé, un réseau de réseaux phasés, un réseau multicanaux, un réseau blazé, un réseau de réseaux blazés et un ensemble de transducteurs à piston. [0019] Dans un autre aspect, il existe un système de sonar actif sous-marin pour mesurer la vitesse d'un instrument par rapport à une surface limite, le système comprenant un transducteur acoustique configuré pour émettre et recevoir une pluralité de faisceaux acoustiques dans différentes directions ; un module émetteur qui émet une séquence codée en phase comprenant des codes répétés de signe alternatif ; et un processeur configuré pour : détecter une surface limite à l'intérieur de chaque faisceau émis ; et mesurer la vitesse d'instrument par rapport à la surface limite. [0020] Le processeur peut être configuré en outre pour filtrer de manière itérative les signaux acoustiques reçus rétrodiffusés à partir des faisceaux émis avec un filtre et la largeur de bande associée qui peut être diminuée successivement pour chaque itération. Le processeur peut être configuré en outre pour détecter une surface limite à l'intérieur de chaque faisceau en trouvant une réponse de crête dans un signal de coefficient de corrélation. Le processeur peut être configuré en outre pour former des faisceaux de manière à séparer les signaux acoustiques reçus. Le processeur peut être configuré en outre pour mesurer la vitesse du véhicule par rétrodiffusion de son sur la surface de fond d'une masse d'eau. Le transducteur acoustique peut émettre une ou plusieurs ré- pétitions de la séquence codée en phase ou un signal à compression d'impulsions pour produire un signal à large bande. Le transducteur acoustique peut comprendre au moins un parmi un réseau phasé, un réseau de réseaux phasés, un réseau multicanaux, un réseau blazé, un réseau de réseaux blazés et un ensemble de transducteurs à piston. [0021] Dans encore un autre aspect, il existe un procédé pour mesurer la vitesse d'un instrument par rapport à une surface limite à l'aide d'un système de sonar actif sous-marin, le procédé comprenant l'émission et la réception d'une pluralité de faisceaux acoustiques dans différentes directions ; la détection d'une surface limite à l'intérieur de chaque faisceau émis ; le filtrage itératif des signaux acoustiques reçus rétrodiffusés à partir des faisceaux émis et de la bande passante associée qui est diminuée successivement pour chaque itération ; et la mesure de la vitesse de l'instrument par rapport à la surface limite. [0022] L'émission et la réception des faisceaux dans différentes directions peuvent avoir lieu depuis au moins un parmi un réseau phasé, un réseau de réseaux phasés, un réseau multicanaux, un réseau blazé, un réseau de réseaux blazés et un ensemble de transducteurs à piston. Le procédé peut comprendre en outre la détection d'une surface limite à l'intérieur de chaque faisceau en trouvant une réponse de crête dans un signal de coefficient de corrélation. La mesure de la vitesse de l'instrument peut comprendre l'estimation d'un décalage Doppler pour chaque faisceau à partir de segments de signaux correspondant à la surface limite, et le calcul de la vitesse de l'instrument à partir du décalage Doppler. Brève description des dessins [0023] [fig.1] est un schéma d'un exemple de sonar Doppler émettant quatre faisceaux dans une configuration Janus. [0024] [fig.21 est un exemple de tracé du spectre en bande de base de sbb(t) pour des fréquences dans la gamme de [-2, 21 Hz et du code de longueur deux. [0025] [fig.31 est un exemple de schéma du spectre en bande de base de sbb(t) pour des fréquences dans la gamme de [-2, 21 Hz et de la séquence de longueur sept. [0026] [fig.41 est un schéma d'un exemple de sonar Doppler ayant un trajet de signal en amplitude et un trajet de signal en phase. [0027] [fig.51 est un schéma d'un exemple de transducteur à réseau phasé pour un enregistreur de vitesse Doppler (DVL). 6 [0028] [fig.61 est un exemple d'organigramme d'un processus de mesure de la vitesse selon un mode de réalisation du sonar Doppler illustré à [Fig 41. [0029] [fig.71 est un exemple d'organigramme d'un processus de détection de contour selon un mode de réalisation du sonar Doppler illustré à [Fig 41. [0030] [fig.81 est un exemple de tracé du spectre en fonction de la fréquence illustrant deux bandes passantes centrées sur deux pics spectraux de la puissance la plus élevée. [0031] [fig.91 est un exemple de schéma illustrant une probabilité de détection simulée en fonction du rapport signal/bruit, normalisée par rapport au mode de réalisation en bande étroite en comparaison du mode de réalisation en longueur 2. [0032] [fig.10] est un exemple d'organigramme d'un processus d'estimation de la vitesse selon un mode de réalisation du sonar Doppler illustré à [Fig 41. [0033] [fig.11] est un exemple de tracé d'un écart-type de l'erreur d'estimation par rapport à un nombre d'itérations pour un DVL fonctionnant à 150 m d'altitude. [0034] DESCRIPTION DETAILLEE DE CERTAINS MODES DE REALISATION ILLUSTRATIFS [0035] La raison d'être du sonar Doppler de poursuite à portée étendue est de fonctionner à une portée similaire à celle d'un DVL à bande étroite, sans toutefois sacrifier l'écart-type et la performance de détection du fond. La première étape inventive consiste en une séquence d'émission composée de codes répétés de modulation par déplacement binaire de phase. Généralement, les séquences de longueur maximale sont utilisées pour la poursuite du fond, mais un code de longueur deux x,[n] = [1, -11 est utilisé dans certains modes de réalisation. Ce choix de code permet une plus longue durée de symbole et une largeur de bande proportionnellement étroite, pour une vitesse d'ambiguïté radiale donnée U,' définie comme suit [Math.61 4 f,LT, [0036] où fc est la fréquence porteuse. Un spectre de 200 de sbb(t) dans [Math 41 est re- présenté sur [Fig 21 pour le code de longueur deux, Ts = 1 s, et Np = 10 pour les fréquences dans la gamme [-2, 21 Hz. Le spectre se compose de pics étroits 210, 220, répétés à une période de 1 Hz et se poursuit indéfiniment avec une fréquence 7 croissante, mais avec une puissance progressivement décroissante. Les pics spectraux les plus forts sont situés aux fréquences f = ±RI2 = ±1/(2T,), où Rs est le débit de symboles. [0037] Un spectre 300 de sbb(t) dans [Math 41 pour la séquence de longueur maximale conventionnelle de longueur sept contient treize pics 310 dans l'intervalle des deux zéros, situés à f = ±R, comme illustré à [Fig 31. De manière similaire, le spectre de sbb(t ) pour la signalisation en bande étroite, tel que défini dans [Math 21, ne contient qu'un seul pic spectral situé à f = 0 Hz pour un décalage Doppler nul. [0038] Les DVL à réseau phasé utilisent généralement une séquence de longueur maximale de longueur sept. Par conséquent, pour une vitesse d'ambiguïté constante U,' la durée du symbole peut être augmentée d'un facteur 3,5 par rapport à un code de longueur sept, ce qui entraîne une diminution de la largeur de bande du bruit et une augmentation du RSB d'un facteur 3,5. [0039] Le développement décrit ne se limite pas à un type particulier de technologie de transducteur et peut être combiné avec des transducteurs à réseau phasé, des transducteurs à piston, des réseaux blazés ou tout autre type de transducteur. Le développement décrit ne se limite pas à l'utilisation de codes composés de deux codets, mais peut également fonctionner avec d'autres longueurs de code. [0040] La largeur de bande du bruit peut être encore réduite en utilisant un filtre adaptatif à bande étroite dont les bandes passantes sont centrées sur les pics spectraux. Le filtrage peut être effectué de manière adaptative et itérative avec une estimation progressivement affinée du décalage Doppler et une largeur de bande progressivement réduite des bandes passantes pour chaque itération, puisque le décalage Doppler n'est généralement pas connu a priori. Le code de longueur deux est plus avantageux que le code de longueur sept pour les applications à longue portée lorsque le RSB est faible, puisque la largeur de bande totale dans les deux pics spectraux est inférieure à la largeur de bande totale dans les treize pics spectraux correspondant au code de longueur sept. [0041] Le code de longueur deux peut être avantageux par rapport à la signalisation en bande étroite car la signalisation en bande étroite est sensible aux fondus d'amplitude, ce qui entraîne une grande variation du RSB pour chaque impulsion sonar. Pour un milieu diffusant avec un grand nombre de diffuseurs, l'amplitude est distribuée selon la loi de Rayleigh, selon [Math.71 [0042] Par conséquent, l'amplitude des impulsions sonar à bande étroite défavorables est 8 souvent beaucoup plus faible que l'amplitude moyenne. Par exemple, la probabilité d'un fondu de 6 dB d'une impulsion sonar à bande étroite est de 17,8 %. Une faible amplitude réduit la probabilité de détection et augmente également la variance de l'erreur d'estimation - deux propriétés indésirables. La signalisation à large bande, en revanche, contribue à des échantillons d'amplitude non corrélés multiples par impulsion sonar, ce qui entraîne une variation moindre de l'amplitude moyenne quadratique sur les impulsions sonar, réduit le fondu et améliore la détection et l'estimation. [0043] L'homme de l'art devrait à présent réaliser qu'un compromis est nécessaire entre la puissance de bruit totale et la variation de l'amplitude, dont la variable est la longueur du code. Le code de longueur deux permet d'améliorer suffisamment la variation d'amplitude, au prix d'une augmentation de la puissance de bruit d'un facteur deux par rapport à la signalisation en bande étroite. Comparé aux codes plus longs, tels que les séquences de longueur maximale de longueur sept, le code de longueur deux fournit une puissance de bruit très réduite en sacrifiant relativement peu la variation d'amplitude. Le raisonnement ne se limite pas aux codes de longueur deux, mais pourrait être étendu à d'autres codes qui offrent un compromis différent entre la puissance du bruit et la variation d'amplitude moyenne. [0044] Les sonars Doppler 400 contiennent souvent deux trajets de signal distincts : l'amplitude et la phase. Le signal de phase est limité, filtré et numérisé et le signal d'amplitude est associé à une fonction logarithmique comme illustré à [Fig 41. Ce cloisonnement est avantageux pour supprimer la grande plage dynamique qui serait sinon nécessaire dans le circuit du récepteur. Les systèmes conventionnels de détection de fond, comme le procédé décrit par Deines et Maier dans le brevet américain n° 5,122,990, utilisent le signal d'amplitude logarithmique pour détecter le fond. Ces procédés fonctionnent bien avec un RSB élevé à modéré, mais le rendement peut être réduit dans les environnements à faible RSB. [0045] Un schéma fonctionnel d'une partie d'un DVL à réseau phasé typique est illustré à [Fig 51. Un seul exemple de transducteur à réseau phasé 500 est utilisé à la fois comme projecteur et hydrophone et est relié à l'électronique par quatre signaux : co(t), coo(t), ro( t) et 1-90(0 (respectivement 418, 416, 414 et 412 sur la figure 4), où r et c représentent respectivement des lignes et des colonnes. Les signaux de colonne sont illustrés à [Fig 51, où chaque signal de colonne unique est relié à chacune des quatre colonnes d'un côté du transducteur piézo-électrique. Par exemple, une colonne 510 correspond au signal de colonne 515, qui est relié à cinq autres colonnes. De même, la colonne 520 correspond au signal de colonne 525, la colonne 530 correspond au signal de colonne 535 et la colonne 540 correspond au signal de colonne 545. Les paires différentielles co -c180 et C90-C270 sont combinées à des signaux simples par l'intermédiaire de trans- 9 formateurs, et les deux signaux de sortie sont nommés c0 et c90. De même, le côté opposé du transducteur piézoélectrique se compose des signaux de ligne (non représentés) qui sont également combinés à des signaux simples par des transformateurs. La même convention d'appellation est appliquée aux signaux de ligne, où les deux sorties sont nommées ro et r90. De plus, le réseau est découpé suivant les deux dimensions qui sont tournées de 900 pour créer les lignes et les colonnes selon un espacement qui peut être fixe ou variable, mais un choix courant consiste à régler l'espacement à une moitié de la longueur d'onde acoustique correspondant à une vitesse nominale du son. [0046] En référence de nouveau à [Fig 41, un réseau de transducteurs ou transducteur acoustique 410 et quatre canaux sont illustrés, où dans cet exemple, le canal un 426 correspond au signal co(t) 418, le canal deux 424 correspond au signal c 9 40 416, le canal trois 422 correspond au signal r0(t) 414 et le canal quatre 420 correspond au signal r90(t) 412. Dans certains modes de réalisation, le transducteur acoustique peut être au moins un parmi un réseau phasé, un réseau de réseaux phasés, un réseau multicanaux, un réseau blazé, un réseau de réseaux blazés et un ensemble de transducteurs à piston. Les signaux d'émission, st,i(t)-st,4(t), sont créés par le générateur de forme d'onde d'émission 490 et amplifiés par des amplificateurs de puissance 470. Les commutateurs d'émission/de réception 430 réalisent un multiplexage entre les signaux d'émission lors de l'émission et les signaux de réception lors de la réception. En mode réception, les signaux du transducteur 410 sont amplifiés par un préamplificateur 435, conçu avec à l'esprit des caractéristiques à faible bruit. La formation du faisceau 440 s'effectue par addition des signaux 00 et 90° avec un déphasage relatif égal à 90° pour les faisceaux 1 et 3 et -90° pour les faisceaux 2 et 4. La sortie du dispositif de formation de faisceaux est divisée en deux trajets : IQ et amplitude. Le trajet IQ se compose d'un amplificateur limiteur d'amplitude 445, d'un démodulateur Zero-IF 450, d'un filtre passe-bas anti-aliasing 460 et d'une numérisation par un convertisseur analogique-numérique (ADC) 465. Le trajet d'amplitude se compose d'un bloc indicateur de l'intensité du signal reçu (RSSI) 455 qui mesure l'intensité du signal de façon logarithmique et d'une numérisation par l'ADC 465. L'amplificateur limiteur d'amplitude 445 dans le trajet IQ est utilisé pour supprimer la plage dynamique de l'entrée, résultant de la grande variation du RSB. Ce cloisonnement particulier du système 400 n'est qu'un exemple de mode de réalisation. D'autres séparations et arrangements pourraient être créés pour mettre en oeuvre la même fonctionnalité que le mode de réalisation décrit. [0047] Le processeur Doppler 475 de [Fig 41 effectue la tâche de calcul du décalage Doppler et de la vitesse d'instrument associée. La vitesse de l'instrument est fournie à l'interface E/S 480 qui peut fournir une entrée et une sortie à un hôte. L'interface E/S 480 est également reliée à un bloc de configuration 485, qui est également relié au processeur Doppler 475 afin de fournir des informations de configuration. Le bloc de configuration 485 est également en communication avec le générateur de forme d'onde d'émission 490 pour les informations de configuration. Le système 400 comprend également une alimentation 495, qui peut être reliée à une source d'alimentation externe. [0048] [fig.61 est un organigramme d'un processus ou d'une méthode 600 pour mesurer la vitesse selon un mode de réalisation de la technologie décrite. Chacune des étapes peut être réalisée par le processeur Doppler comme illustré à [Fig 41. Dans la description du mode de réalisation, on suppose que la projection est mise en place avec une impulsion émise qui est composée d'une multitude de séquences répétées à un intervalle de temps T L . En outre, on suppose que les faisceaux Janus sont créés soit par une multitude de transducteurs à piston, soit par la formation de faisceaux de signaux de réseau phasé ou par d'autres moyens. [0049] Étape 1. Après un début 605, dans le bloc 610, la méthode 600 filtre les signaux reçus pour chaque faisceau en fonction d'une forme et d'une largeur de bande de filtre désirées. La forme et la largeur de bande du filtre doivent être choisies de manière adaptée au spectre des signaux reçus, y compris tout décalage Doppler. [0050] Étape 2. Dans le bloc 620, la méthode 600 détermine un ou plusieurs contours, tels que le fond marin, la surface ou les calottes glaciaires, pour chaque faisceau et sélectionne les segments de données à utiliser. En utilisant l'intensité du signal, la puissance du signal, le coefficient de corrélation ou la covariance, la méthode 600 détermine les segments appropriés à utiliser pour l'estimation Doppler dans chaque faisceau. La détection d'un contour peut être effectuée en trouvant la valeur maximale, un pic local, une certaine forme, un pic à la sortie d'un filtre adapté à une certaine forme, ou une autre mesure. Selon l'application, un autre type de contour peut être souhaitable pour la détection, ou plusieurs contours peuvent être détectés. [0051] Étape 3. Dans le bloc 630, la méthode 600 estime le décalage Doppler pour chaque faisceau à partir de segments correspondant aux contours détectés. La méthode de covariance est utilisée pour calculer la phase Doppler OD selon [Math. 81 x* [ni Ns FE 11 [Math.91 Lc[n] [0052] où x[m] = xi[m] + jxQ[m] est le signal numérisé complexe reçu pour chaque faisceau, * indique une conjugaison complexe, N est le nombre d'échantillons dont la moyenne est calculée et NspE est le nombre d'échantillons par codet. [0053] Étape 4. Dans le bloc 640, la méthode 600 calcule la vitesse radiale pour chaque faisceau à partir du décalage Doppler, ce qui est accompli par [Math.10] [0054] La vitesse d'un instrument ou d'un référentiel terrestre peut également être calculée en transformant les vitesses radiales en un système de coordonnées instrument-référentiel ou en un système de coordonnées aligné avec la terre. La méthode 600 se termine à une fin 645. [0055] [fig.71 est un organigramme d'un processus ou d'une méthode 700 pour la détection de contour selon un mode de réalisation de la technologie décrite. Chacune des étapes peut être effectuée par le processeur Doppler comme illustré [Fig 41. Dans la description du mode de réalisation, on suppose que la projection est mise en place avec une impulsion émise qui est composée d'une multitude de séquences répétées à un intervalle de temps T L . En outre, on suppose que les faisceaux Janus sont créés soit par une multitude de transducteurs à piston, soit par la formation de faisceaux de signaux de réseau phasé ou par d'autres moyens. Les étapes de l'organigramme de [Fig 71 correspondent à un faisceau, mais sont exécutées pour chaque faisceau indépendamment. [0056] Étape 1. Après un début 705, dans le bloc 710, la méthode 700 filtre les signaux reçus en fonction d'une forme et d'une largeur de bande de filtre désirées. La forme et la largeur de bande du filtre doivent être choisies de manière adaptée au spectre des signaux reçus, y compris tout décalage Doppler. Un avantage du code de longueur deux est que la porteuse est supprimée, ce qui permet de filtrer 800 le signal reçu selon [Fig 81, qui contient deux bandes passantes 810 et 820 centrées sur les deux pics spectraux de puissance maximale 815 et 825, respectivement. Cette technique de filtrage réduit encore la puissance de bruit vue par le processeur Doppler, en réduisant la largeur de bande du bruit. [0057] Étape 2. Dans le bloc 720, la méthode 700 calcule le coefficient de corrélation de 12 l'échantillon sur une fenêtre d'échantillons de longueur W selon [Math.11] r[k] in-FLN-1,& W-1 qi 2 JI [0058] pour k = [1, W, 2W,..., (N-1)W1. Il existe des variations de l'équation pour calculer le coefficient de corrélation de l'échantillon et la version particulière illustrée dans l'équation 11 ne devrait pas limiter la portée du mode de réalisation. [0059] Étape 3. Dans le bloc 730, la méthode 700 calcule la moyenne du coefficient de cor- rélation de l'échantillon sur M échantillons adjacents de r[k] et calcule la valeur absolue. Le calcul de la moyenne peut être mis en oeuvre par un filtre de longueur M, dans lequel les pondérations des prises sont de 11M. L'avantage des étapes 2 et 3 est d'échanger la complexité de calcul requise contre une résolution dans ra[k], la valeur absolue du coefficient de corrélation moyen de l'échantillon. Le choix de W définit la résolution et, pour les applications sonar Doppler, elle peut être réglée grossièrement, puisque le but de la détection de contour est de déterminer un segment d'un signal pour calculer le décalage Doppler au sein d'une plage de mesure, par opposition à une mesure très précise, par rapport au contour. [0060] Étape 4. Dans le bloc 740, la méthode 700 filtre la valeur absolue du coefficient de corrélation moyen de l'échantillon ra[k] avec un filtre accordé sur la forme attendue d'un écho de contour dans l'intervalle de ra[k]. Il est typique de sélectionner W et M de manière à ce que le produit WM corresponde à la durée d'émission Tt, = LNpTs. La ra[k] à forme attendue est alors un triangle et le filtre de choix est un triangle dont la durée est égale à la durée d'émission. [0061] Étape 5. Dans le bloc 750, la méthode 700 recherche la valeur maximale de la sortie du filtre à l'étape 4. L'échantillon correspondant au maximum est affecté à [0062] Étape 6. Dans le bloc 760, la méthode 700 détermine si le coefficient de corrélation de l'échantillon autour de l'indice correspondant à la valeur maximale dépasse un seuil de détection D. Plus précisément, la détection est déclarée si ra[k] dépasse D pour au moins un échantillon k dans la plage kmax-M