FR 3086064 A1 20200320 FR 1858281 A 20180914 DOMAINE TECHNIQUE AUQUEL SE RAPPORTE L'INVENTION La présente invention concerne de manière générale la détermination d'une différence de profondeur, ou la détermination d'un écart latéral par rapport à la verticale entre deux points d'un milieu subaquatique, par mesure d'un temps de propagation et/ou d'une constante de propagation d'une onde acoustique se propageant entre ces deux points, sur la base d'une célérité moyenne dans ce milieu subaquatique. Elle s'applique en particulier à la mesure de la topographie d'un fond marin, et à la localisation d'un objet ou d'un banc de poissons immergé. Elle s'applique également à la détermination de ladite célérité moyenne, à partir des mesures topographiques elles-mêmes. ARRIERE-PLAN TECHNOLOGIQUE De plus en plus de navires sont équipés aujourd'hui de sonars, qu'il s'agisse de drones ou de navires destinés à accueillir un équipage, de surface ou submersibles. En particulier, de plus en plus de drones de type « AUV » ou « USV » (selon des acronymes anglo-saxons de « Autonomous Underwater Vehicule », c'est-à-dire « navire autonome sous-marin », et de « Unmanned Surface Vehicule », c'est-à-dire navire de surface sans équipage) sont aujourd'hui équipés de sonars. De tels sonars permettent notamment de déterminer un profil bathymétrique du fond marin situé sous le navire (un tel profil rassemblant, pour plusieurs points du fond marin, des données de profondeurs et de positions horizontales repérant ces points du fond), ou de localiser un objet ou un banc de poissons situé dans une colonne d'eau s'étendant sous ce navire. Pour déterminer la profondeur d'un point du fond marin situé à l'aplomb du sonar, le sonar émet une onde acoustique vers ce point. Après réflexion sur le fond, cette onde acoustique est reçue par un ou plusieurs récepteurs équipant le sonar. La durée séparant les instants d'émission et de réception de l'onde acoustique, qui est égale au temps nécessaire pour que l'onde acoustique effectue un trajet aller et un trajet retour entre le sonar et le fond marin, permet alors de déterminer la profondeur à laquelle est située ce point du fond marin. Pour cela, la moitié de cette durée est multipliée par une valeur donnée d'une célérité des ondes acoustiques dans le milieu. Mais une telle détermination de profondeur est souvent imprécise car la célérité des ondes acoustiques dans un milieu subaquatique varie généralement 2 en fonction de la profondeur (du fait notamment de variations de température, de densité ou de salinité avec la profondeur). Le « Manuel d'hydrographie » de l'Organisation Hydrographique Internationale (publication C-13, 1 è`e édition, mise à jour de février 2011, 5 disponible sur le site internet https://www.iho.int/iho_pubs/CB/C13_Index.htm#C- 13F), propose alors, pour déterminer la profondeur d'un point du fond situé à l'aplomb du sonar : - préalablement aux mesures de temps d'échos mentionnées ci-dessus, de relever un profil de célérité, c'est-à-dire un ensemble de valeurs de célérité 10 présentées localement, en différents points situés plus ou moins profondément sous la surface de l'eau, par des ondes acoustiques, puis - d'employer une célérité effective, égale à la moyenne harmonique, le long de la colonne d'eau considérée, des valeurs de célérité précédemment relevées (paragraphe 2.2.4 du chapitre 3 de ce document). 15 Puisqu'il s'agit d'une moyenne harmonique, l'inverse de cette célérité effective est égal à la moyenne arithmétique des inverses desdites valeurs locales de célérité du son. Pour un point du fond qui n'est pas situé à l'aplomb du sonar, l'onde acoustique est déviée par réfraction, tout au long de son trajet entre le sonar et le 20 point sondé (du fait de la variation de la célérité des ondes acoustiques avec la profondeur). C'est le temps de propagation de l'onde, le long de ce trajet non-rectiligne, qui est alors mesuré par le sonar. Pour déterminer la profondeur d'un tel point, et son écart latéral par rapport au sonar (écart par rapport à la verticale), le document précité indique (en 25 page 121) qu'il est nécessaire de tracer le trajet suivi point par point par les rayons acoustiques, déviés sous l'effet de la réfraction mentionnée plus haut. Ce document propose plus précisément (en page 161) de modéliser le milieu subaquatique comme un empilement de plusieurs couches, la célérité des ondes acoustiques ou son gradient étant supposés constants dans chaque couche. Le 30 trajet des rayons acoustiques est alors déterminé de proche en proche, couche par couche. Cela permet alors de relier le temps de propagation le long de ce trajet, à la profondeur du point sondé et à son écart latéral par rapport au sonar. Cependant, le trajet suivi par les rayons acoustiques doit être déterminé à nouveau, entièrement, à chaque changement de la direction de l'onde 3 acoustique reçue par le sonar. La mise en oeuvre de cette méthode implique donc de nombreux calculs, et nécessite des ressources informatiques importantes, en particulier lors de l'utilisation de sonars multifaisceaux (qui génèrent simultanément une multitude 5 d'ondes acoustiques d'inclinaisons différentes). Cette méthode nécessite en outre une mesure préalable du profil de célérité des ondes acoustiques sur toute la gamme de profondeurs allant du sonar jusqu'au point sondé. OBJET DE L'INVENTION Dans ce contexte, la présente invention propose un procédé de 10 détermination d'une différence de profondeur, ou d'un écart latéral par rapport à la verticale entre deux points d'un milieu subaquatique, basé sur un modèle monocouche bien choisi du milieu subaquatique, au lieu d'être basé sur un modèle multicouche tel que décrit dans l'art antérieur précité (« Manuel d'hydrographie » de l'Organisation Hydrographique Internationale). Dans ce 15 modèle monocouche, il est supposé que l'onde acoustique se propage en ligne droite dans le milieu, selon une direction de propagation effective (différente de sa direction de réception), et avec une célérité moyenne qui est indépendante de cette direction de propagation effective. L'emploi de ce modèle monocouche simplifie considérablement les calculs nécessaires à une telle détermination, tout 20 en conduisant à des résultats d'une grande précision. L'invention prévoit en particulier un procédé comprenant les étapes suivantes : - émission d'une onde acoustique dans le milieu subaquatique au moyen d'au moins un émetteur, 25 - réception de ladite onde acoustique au moyen d'une antenne réceptrice comprenant plusieurs récepteurs, ces récepteurs délivrant une pluralité respective de signaux de réception lors de la réception de ladite onde acoustique, - détermination d'une constante de propagation de l'onde acoustique reçue, en fonction desdits signaux de réception, ladite constante de propagation 30 étant égale au sinus d'un angle de réception repérant la direction de réception de l'onde acoustique par rapport à la verticale, divisé par une vitesse de propagation locale des ondes acoustiques à la profondeur de ladite antenne réceptrice, - détermination d'un temps de propagation de l'onde acoustique, en fonction d'une durée séparant les instants d'émission et de réception de l'onde 4 acoustique, - détermination de la différence de profondeur ou de l'écart latéral par rapport à la verticale, entre l'antenne réceptrice et l'émetteur, ou entre l'antenne réceptrice et un élément immergé réfléchissant ladite onde acoustique lors de sa 5 propagation de l'émetteur jusqu'à l'antenne réceptrice, en fonction: - d'un produit dudit temps de propagation par ladite célérité moyenne de l'onde acoustique, et en fonction - d'un angle de propagation moyen, défini entre la verticale et la direction de propagation effective de l'onde acoustique, l'angle de propagation 10 moyen étant déterminé en fonction de ladite constante de propagation. Le déposant a constaté que pour une valeur appropriée de ladite célérité moyenne, la différence de profondeur, ou l'écart latéral déterminés comme indiqué ci-dessus sont, de manière surprenante, suffisamment précis au regard des normes hydrographiques internationales généralement recommandées (en 15 particulier au regard des ordres de précision « 1 b » et « spécial » tels que définis par les « Normes de l'O.H.I pour les levés hydrographiques », Monaco, 5è" édition, ordres de précision repris en pages 8 à 10 du « Manuel d'hydrographie » de l'Organisation Hydrographique Internationale mentionné plus haut). Dans ce procédé, l'influence de la réfraction subie par l'onde acoustique 20 est donc prise en compte précisément, sans pour autant qu'il soit nécessaire de déterminer le détail du trajet suivi par l'onde acoustique. Cela permet d'exécuter le procédé de manière avantageusement rapide et/ou avec des ressources de calcul limitées. Il s'avère en outre, comme déjà indiqué, que la valeur appropriée de la 25 célérité moyenne en question (pour laquelle on obtient une détermination précise de la différence de profondeur et/ou de l'écart latéral entre les points considérés) est indépendante de l'inclinaison de l'onde acoustique reçue. Ainsi, à partir du moment où cette valeur appropriée est connue (ou déterminée), plusieurs mesures de différences de profondeur, réalisées pourtant pour des inclinaisons 30 différentes, peuvent être obtenues au moyen d'une même formule numérique, ce qui réduit notablement le temps de calcul nécessaire pour déterminer ces différences de profondeur. On note par ailleurs que la différence de profondeur et/ou l'écart latéral mentionnés ci-dessus peuvent être déterminés en se basant sur une seule valeur 5 de célérité, au lieu d'avoir à se baser sur de nombreuses célérités différentes mesurées préalablement pour toute la gamme de profondeurs allant de l'antenne réceptrice jusqu'à l'émetteur, ou de l'antenne réceptrice jusqu'audit élément immergé. 5 D'autres caractéristiques non limitatives et avantageuses du procédé de détermination conforme à l'invention, prises individuellement ou selon toutes les combinaisons techniquement possibles, sont les suivantes : - il est prévu de déterminer à la fois ladite différence de profondeur et ledit écart latéral par rapport à la verticale ; 10 - ladite différence de profondeur est déterminée, de manière à présenter un écart relatif inférieur à un pour mille par rapport : - au produit dudit temps de propagation par ladite célérité moyenne, - multiplié par le cosinus dudit angle de propagation moyen ; - ledit écart latéral est déterminé, de manière à présenter un écart relatif 15 inférieur à un pour mille par rapport : - au produit dudit temps de propagation par ladite célérité moyenne, - multiplié par le sinus dudit angle de propagation moyen ; - ladite célérité moyenne est représentative d'une moyenne harmonique d'une pluralité de vitesses de propagation locales, présentées par les ondes 20 acoustiques à une pluralité respective de profondeurs allant d'une profondeur de l'émetteur jusqu'à une profondeur de l'antenne réceptrice, ou d'une profondeur dudit élément immergé jusqu'à la profondeur de l'antenne réceptrice ; - l'angle de propagation moyen est déterminé de sorte que son sinus présente un écart relatif inférieur à un pour mille par rapport au produit de ladite 25 constante de propagation par ladite célérité moyenne ; - l'angle de propagation moyen est déterminé en fonction en outre d'une célérité arithmétique moyenne qui est représentative d'une moyenne arithmétique d'une pluralité de vitesses de propagation locales, présentées par les ondes acoustiques à une pluralité respective de profondeurs allant de la profondeur 30 l'émetteur jusqu'à la profondeur de l'antenne réceptrice, ou de la profondeur dudit élément immergé jusqu'à la profondeur de l'antenne réceptrice ; - l'angle de propagation moyen est déterminé de sorte que son sinus présente un écart relatif inférieur à un point mille par rapport à ladite constante de propagation, multipliée par ladite célérité arithmétique moyenne ; 6 - l'émetteur et l'antenne réceptrice équipent un même système sonar, et, au cours du procédé, ce système sonar détermine la différence de profondeurs ou l'écart latéral par rapport à la verticale entre l'antenne réceptrice et un élément immergé réfléchissant ladite onde acoustique lors de sa propagation depuis 5 l'émetteur jusqu'à l'antenne réceptrice ; - ledit émetteur fait partie d'une antenne émettrice directive dudit système sonar, cette antenne émettrice émettant ladite onde acoustique et comprenant au moins un autre émetteur ; - l'émetteur et l'antenne réceptrice équipent respectivement deux 10 systèmes distincts situés à des profondeurs respectives différentes, et, au cours du procédé, le système muni de l'antenne réceptrice détermine la différence de profondeurs ou l'écart latéral par rapport à la verticale entre l'antenne réceptrice et ledit émetteur ; - ladite antenne réceptrice comprend au moins trois récepteurs non 15 alignés, ou, le système comprenant ladite antenne réceptrice comprend au moins une autre antenne réceptrice, une direction de réception de ladite onde acoustique étant repérée de manière tridimensionnelle sur la base des signaux de réception délivrés par l'ensemble des deux antennes réceptrices ; - l'émetteur équipe une balise immergée ; 20 - l'antenne réceptrice équipe un système sonar ; - le procédé comprend en outre les étapes suivantes : - détermination, au moyen d'une sonde positionnée successivement à plusieurs profondeurs différentes, d'une pluralité de vitesses de propagation locales présentées respectivement, à chacune desdites profondeurs, par les 25 ondes acoustiques, puis - détermination de ladite célérité moyenne en fonction de ladite pluralité de vitesses de propagation locales déterminées précédemment ; - le procédé comprend en outre les étapes suivantes : - après l'émission et la réception de ladite onde acoustique, 30 déplacement du système sonar, puis - émission d'une autre onde acoustique dans le milieu subaquatique au moyen dudit émetteur, de sorte que ladite autre onde atteigne ledit élément immergé, - après réflexion sur ledit élément immergé, réception de ladite autre 7 onde acoustique au moyen de l'antenne réceptrice, - détermination d'une autre constante de propagation, pour ladite autre onde acoustique reçue, en fonction des signaux de réception délivrés par les récepteurs de l'antenne réceptrice suite à la réception de cette autre onde 5 acoustique, ladite autre constante de propagation étant égale au sinus d'un autre angle de réception repérant la direction de réception de ladite autre onde acoustique par rapport à la verticale, divisé par une vitesse de propagation locale des ondes acoustiques à la profondeur de ladite antenne réceptrice, - détermination d'un autre temps de propagation, en fonction d'une 10 autre durée séparant les instants d'émission et de réception de ladite autre onde acoustique, et - détermination de ladite célérité moyenne en fonction : - du temps de propagation et de la constante de propagation de ladite onde acoustique, et 15 - du temps de propagation et de la constante de propagation de ladite autre onde acoustique ; - ladite célérité moyenne est déterminée de manière à ce qu'un écart entre : - ladite différence de profondeur, et 20 - une autre différence de profondeur entre l'antenne réceptrice et l'élément immergé soit inférieur à un seuil donné, ladite autre différence de profondeur étant déterminée en fonction : - d'un produit dudit autre temps de propagation par ladite 25 célérité moyenne, et en fonction - d'un autre angle de propagation moyen, défini entre la verticale et une direction de propagation effective de ladite autre onde acoustique, ledit autre angle de propagation moyen étant déterminé en fonction de ladite autre constante de propagation. 30 L'invention concerne aussi un procédé de détermination d'une vitesse de propagation locale des ondes acoustiques dans un milieu subaquatique, comprenant les étapes suivantes : - pour une pluralité d'éléments immergés situés à des profondeurs différentes, détermination d'une pluralité respective de valeurs de ladite célérité 8 moyenne, chacune desdites valeurs étant déterminée conformément au procédé décrit ci-dessus, puis - détermination d'au moins une vitesse de propagation locale des ondes acoustiques, à une profondeur donnée fixe, en fonction de ladite pluralité de 5 valeurs de la célérité moyenne déterminées précédemment. L'invention concerne aussi un procédé de détermination d'une position d'un émetteur immergé dans un milieu subaquatique, comprenant les étapes suivantes : - émission d'une onde acoustique dans le milieu subaquatique par ledit 10 émetteur, dont la profondeur est connue, - réception de ladite onde acoustique au moyen d'une antenne réceptrice équipant un navire, de surface ou submersible, l'antenne réceptrice comprenant plusieurs récepteurs délivrant une pluralité respective de signaux de réception, - détermination d'une constante de propagation de l'onde acoustique 15 reçue, en fonction desdits signaux de réception, ladite constante de propagation étant égale au sinus d'un angle de réception repérant la direction de réception de l'onde acoustique par rapport à la verticale dans le plan vertical contenant l'émetteur et l'antenne réceptrice, divisé par une vitesse de propagation locale des ondes acoustiques à la profondeur de ladite antenne réceptrice, 20 - détermination d'une direction de propagation effective de l'onde acoustique en fonction du produit de ladite constante de propagation et d'une célérité moyenne de propagation entre la profondeur de l'émetteur et la profondeur de l'antenne réceptrice, et - détermination d'un écart latéral par rapport à la verticale entre 25 l'antenne réceptrice et l'émetteur, en fonction du produit de : - la différence entre la profondeur de l'émetteur et la profondeur de l'antenne réceptrice, et de - la tangente d'un angle de propagation moyen défini entre ladite direction de propagation effective et la verticale. 30 L'émetteur peut en particulier équiper une balise immergée, tandis que l'antenne réceptrice équipe un système sonar. L'invention concerne également un système sonar comprenant : - une antenne réceptrice configurée pour réceptionner une onde acoustique qui a été émise précédemment dans un milieu subaquatique par au 9 moins un émetteur, l'antenne réceptrice comprenant plusieurs récepteurs aptes à délivrer une pluralité respective de signaux de réception lors de la réception de ladite onde acoustique, et - une unité de traitement électronique programmée pour : 5 - déterminer une constante de propagation de l'onde acoustique reçue, en fonction desdits signaux de réception, ladite constante de propagation étant égale au sinus d'un angle de réception repérant la direction de réception de l'onde acoustique par rapport à la verticale, divisé par une vitesse de propagation locale des ondes acoustiques à la profondeur de ladite antenne réceptrice, 10 - déterminer un temps de propagation de l'onde acoustique, en fonction d'une durée séparant les instants d'émission et de réception de l'onde acoustique, et pour - déterminer une différence de profondeur ou un écart latéral par rapport à la verticale, entre l'antenne réceptrice et l'émetteur ou entre l'antenne réceptrice 15 et un élément immergé réfléchissant ladite onde acoustique lors de sa propagation depuis l'émetteur jusqu'à l'antenne réceptrice, en fonction : - d'un produit dudit temps de propagation par une célérité moyenne de l'onde acoustique, et en fonction - d'un angle de propagation moyen, défini entre la verticale et une 20 direction de propagation effective de l'onde acoustique, l'angle de propagation moyen étant déterminé en fonction de ladite constante de propagation. Le système sonar peut comprendre aussi ledit émetteur, ou comprendre une antenne émettrice qui comprend ledit émetteur au nombre de un au moins. Les différentes caractéristiques optionnelles présentées plus haut en 25 termes de procédé peuvent aussi s'appliquer au système sonar qui vient d'être décrit. L'invention concerne aussi un système sonar comprenant : - une antenne réceptrice configurée pour réceptionner une onde acoustique qui a été émise précédemment dans un milieu subaquatique par un 30 émetteur dont on connait la profondeur, l'antenne réceptrice comprenant plusieurs récepteurs aptes à délivrer une pluralité respective de signaux de réception lors de la réception de ladite onde acoustique, et - une unité de traitement électronique programmée pour : - déterminer une constante de propagation de l'onde acoustique 10 reçue, en fonction desdits signaux de réception, ladite constante de propagation étant égale au sinus d'un angle de réception repérant la direction de réception de l'onde acoustique par rapport à la verticale, divisé par une vitesse de propagation locale des ondes acoustiques à la profondeur de ladite antenne réceptrice, 5 déterminer une direction de propagation effective de l'onde acoustique reçue, en fonction du produit de ladite constante de propagation et d'une la célérité moyenne de propagation entre la profondeur de l'émetteur et la profondeur de l'antenne réceptrice, - acquérir ou lire dans une mémoire une donnée représentative de la 10 profondeur de l'émetteur, et pour - déterminer l'écart latéral par rapport à la verticale, entre l'antenne réceptrice et l'émetteur, en fonction du produit de : - la différence entre la profondeur de l'émetteur et la profondeur de l'antenne réceptrice, et de 15 - la tangente d'un angle de propagation moyen défini entre ladite direction de propagation effective et la verticale. L'invention concerne également un navire, autonome ou destiné à accueillir un équipage, de surface ou submersible, muni d'un système sonar tel que décrit ci-dessus. 20 DESCRIPTION DETAILLEE D'UN EXEMPLE DE REALISATION La description qui va suivre en regard des dessins annexés, donnés à titre d'exemples non limitatifs, fera bien comprendre en quoi consiste l'invention et comment elle peut être réalisée. Sur les dessins annexés : 25 - la figure 1 représente schématiquement vu de côté, un navire muni d'un système sonar selon l'invention, - la figure 2 représente schématiquement le navire de la figure 1, vu de dessous, - la figure 3 représente schématiquement un faisceau d'ondes 30 acoustiques émis par le système sonar de la figure 1 au cours d'un premier procédé permettant de relever la topographie d'un fond marin, - la figure 4 représente schématiquement, plus en détail, des récepteurs du système sonar de la figure 1, - la figure 5 représente schématiquement le trajet suivi par une onde 11 acoustique reçue par les récepteurs du système sonar de la figure 1, - la figure 6 représente schématiquement les principales étapes du premier procédé mentionné ci-dessus, - les figures 7A et 7B représentent schématiquement l'évolution d'une 5 vitesse de propagation locale des ondes acoustiques dans le milieu subaquatique considéré, et l'évolution d'une célérité moyenne de ces ondes en fonction de la profondeur dans ce milieu, - la figure 8 représente schématiquement une erreur de détermination portant sur une profondeur, et sur une position latérale d'un point du fond, sondé 10 par le système sonar de la figure 1, - la figure 9 représente schématiquement des faisceaux d'ondes acoustiques, émis successivement par le système sonar de la figure 1 au cours d'un deuxième procédé, permettant à la fois de relever la topographie d'un fond marin et de déterminer un profil de célérité moyenne des ondes acoustiques dans 15 le milieu sondé, - la figure 10 représente schématiquement des trajets suivis par des ondes acoustiques émises au cours de ce deuxième procédé, - la figure 11 représente schématiquement un écart entre deux estimations distinctes d'une même profondeur, en fonction de la célérité moyenne 20 employée pour déterminer ces deux estimations, - la figure 12 représente schématiquement des valeurs de célérité moyenne déterminées au moyen de ce deuxième procédé, pour différentes profondeurs, ainsi que des valeurs de référence de cette célérité connues par ailleurs, 25 - la figure 13 représente schématiquement un profil de vitesse de propagation locale des ondes acoustiques, déterminé au cours de ce deuxième procédé, - la figure 14 représente schématiquement le trajet d'une onde acoustique émise par une balise immergée, et reçue par le système sonar de la 30 figure 1, vu de dessus, - la figure 15 représente schématiquement le trajet de l'onde acoustique de la figure 14, vue de côté. Les figures 1 et 2 représentent schématiquement un système sonar 10 équipant un navire 1, en l'occurrence un navire de surface. 12 Un aspect remarquable de ce système sonar 10 est qu'il est configuré de manière à déterminer des positions de points sondés (profondeurs et positions latérales des points sondés) en tenant compte, de manière à la fois simple et précise, de la réfraction des ondes acoustiques causée par les variations de la 5 vitesse de propagation de ces ondes en fonction de la profondeur. La structure de ce système sonar 10 autorise différents modes de fonctionnement. Elle permet notamment de : - déterminer un profil bathymétrique du fond marin, c'est-à-dire de réaliser un relevé de la topographie de ce fond, sur la base d'un profil de célérité 10 moyenne du son dans le milieu subaquatique relevé au préalable (premier procédé, décrit en référence aux figures 3 à 8), - déterminer un tel profil bathymétrique, ainsi qu'un profil de célérité moyenne du son, ce profil de célérité moyenne étant déterminé sur la base des relevés bathymétriques eux-mêmes (deuxième procédé, décrit en référence aux 15 figures 9 à 13), et - déterminer une position d'une balise, ou d'un autre système immergé, par rapport au système sonar 10 (troisième procédé, décrit en référence aux figures 14 et 15). La structure du système sonar 10 sera décrite dans un premier temps, 20 en référence aux figures 1 à 5, et ces différents modes de fonctionnement seront décrits ensuite. Système sonar Le système sonar 10, représenté schématiquement sur la figure 1, comprend : 25 - une tête sonar 11, comportant plusieurs transducteurs 15, 16 ainsi qu'une unité de commande de ces transducteurs (non représentée), - une unité de pilotage 17 de la tête sonar 11, et - une unité de traitement 18 des signaux acquis par la tête sonar. L'unité de commande des transducteurs 15, 16 peut comprendre des 30 convertisseurs numérique-analogique et analogique-numérique, ainsi que des amplificateurs et filtres électroniques, adaptés pour mettre en forme des signaux électriques alimentant les transducteurs ou pour amplifier et filtrer des signaux de réception captés par certains de ces transducteurs, utilisés en réception. Les transducteurs 15, 16 du système sonar 10 sont disposés en croix 13 (disposition dite en « croix de Mills ») : - certains transducteurs 15 sont disposés en ligne les uns à la suite des autres, le long d'une première branche de la croix, et forment ensemble une première antenne 13, tandis que 5 - les autres transducteurs 16 sont disposés les uns à la suite des autres le long d'une deuxième branche de la croix, perpendiculaire à sa première branche, et forment ensemble une deuxième antenne 14 (figure 2). La distance d entre deux transducteurs 15 voisins de la première antenne 13 est constante le long de cette antenne, et il en est de même pour les 10 transducteurs 16 de la deuxième antenne 14. La tête sonar 11 est ici insérée dans un logement formé dans la coque du navire 1, et elle est orientée de sorte que la première antenne 13 (la première branche de la croix de Mills) soit parallèle à un axe longitudinal x du navire 1, la deuxième antenne 14 étant quant à elle parallèle à un axe transversal y du navire. 15 La coque du navire présente une forme allongée selon l'axe longitudinal x en question, et, lorsque le navire se déplace en ligne droite, sa direction de marche coïncide avec cet axe longitudinal x (effets de dérive mis à part). L'axe transversal y est perpendiculaire à l'axe longitudinal x, et parallèle au pont du navire. 20 Ici, chaque transducteur 15, 16 du système sonar 10 est apte à émettre des ondes acoustiques dans le milieu subaquatique E environnant le navire, mais aussi à recevoir des ondes acoustiques provenant de ce milieu (en particulier des ondes acoustiques réfléchies par le fond marin). Chacun de ces transducteurs 15, 16 est donc apte à fonctionner aussi bien en tant qu'émetteur qu'en tant que 25 récepteur. Chacune des première et deuxième antenne 13 et 14 présente une longueur qui, de préférence, est supérieure à 20 centimètres, voire supérieure à 50 centimètres, grâce à quoi cette antenne peut émettre ou recevoir des ondes acoustiques de manière directive, avec une résolution angulaire élevée. 30 Le caractère directif de ces antennes est présenté ci-dessous, dans le cas où la première antenne 13 est employée en tant qu'antenne émettrice et où la deuxième antenne 14 est employée en tant qu'antenne réceptrice. Dans ce cas, la première antenne 13 émet des ondes acoustiques formant un faisceau 30 d'ondes acoustiques, en quelque sorte aplati, qui : 14 - présente une grande ouverture angulaire, de plusieurs dizaines de degrés, dans un plan d'émission perpendiculaire à la première antenne 13, - mais qui présente une ouverture angulaire réduite (par exemple inférieure à 5 degrés) perpendiculairement à ce plan d'émission (C.f. : figure 3). 5 Ce faisceau 30 d'ondes acoustiques se propage ainsi sous la forme d'une nappe peu épaisse, en forme d'éventail, généralement appelée « fauchée ». Cette émission est directive, en ce sens que le faisceau émis présente, parallèlement à l'antenne émettrice (donc ici parallèlement à l'axe longitudinal x), ladite ouverture angulaire réduite. La zone du fond 4 marin atteinte par ce faisceau 10 30 correspond alors à une bande 32 mince et longue, perpendiculaire à l'antenne émettrice. Cette bande 32 est constituée de plusieurs petites portions élémentaires du fond 4 marin, 33, 33', 33"... réparties le long de cette bande, et centrées chacune sur un point P, P', P"... . Chacun de ces portions élémentaire 33, 33', 33" 15 du fond marin réfléchit une onde acoustique, qui, à la manière d'un écho, est renvoyée vers le système sonar 10. Cette onde réfléchie est ensuite réceptionnée par la deuxième antenne 14. Grâce au caractère directif de la deuxième antenne 14, le système sonar 10 est capable de distinguer l'une de l'autre deux ondes acoustiques reçues, 20 provenant respectivement de portions 33, 33', 33" du fond marin distinctes, qui occupent des positions différentes le long de la bande 32 insonifiée mentionnée ci-dessus. Cela permet alors de déterminer la profondeur et la position latérale de chacune de ces portions élémentaires, c'est-à-dire la profondeur et la position latéral du point P, P', P" sur lequel elles sont centrées. 25 Le caractère directif de la deuxième antenne 14 (employée ici en réception) est expliqué ci-dessous, en référence à la figure 4. Chaque transducteur 16 de cette antenne est apte à délivrer un signal de réception s, représentatif des variations instantanées, à la position où se trouve ce transducteur, d'une onde acoustique reçue par l'antenne. Lorsque cette onde n'est 30 pas reçue en incidence normale, ces signaux de réception s sont décalés temporellement les uns par rapport aux autres, d'une quantité qui dépend directement de l'angle de réception 0,, c'est à dire de l'angle d'incidence sous lequel l'onde acoustique est reçue par l'antenne. En l'occurrence, le décalage temporel At entre deux desdits signaux de réception s est donné par la formule F1 15 suivante : At = al / co - d.sin(0,)= d. a (F1). Go Dans cette formule, c, est la vitesse de propagation locale des ondes acoustiques, à la profondeur z, de la deuxième antenne 14 (profondeur z, qui est 5 ici quasiment nulle). Quant à l'angle de réception 0,, il est formé entre une direction perpendiculaire à l'antenne, qui coïncide ici avec l'axe vertical z, et une direction de propagation présentée par l'onde acoustique au niveau de l'antenne 14 (c'est-à-dire là où est située l'antenne 14). Un traitement approprié des signaux de réception s permet alors 10 d'extraire de ces signaux une composante représentative de l'onde, ou des ondes acoustiques, reçues par l'antenne sous un angle de réception 0, donné. Des ondes reçues respectivement sous des angles d'incidence différents (et provenant donc de points P, P', P" du fond marin différents) peuvent ainsi être distinguées les unes des autres par le système sonar 10. Le traitement en question, qui est un 15 traitement classique de formation de voies, peut par exemple consister à appliquer, à chaque signal de réception s, un décalage temporel proportionnel à la position occupée par le transducteur correspondant le long de l'antenne 14, puis à sommer les signaux de réception ainsi décalés. Comme le montre la formule F1 ci-dessus, c'est sous la forme d'une 20 constante de propagation a que l'angle de réception 0, de l'onde reçue intervient, dans le décalage temporel At entre deux des signaux de réception s. Cette constante de propagation a, généralement appelée constante de Snell-Descartes, est égale au sinus de l'angle de réception 0, divisé par la vitesse de propagation locale c, des ondes acoustiques à la profondeur z, de la deuxième 25 antenne 14 : a = sin(00) /c0. Cette quantité est appelée constante de propagation, car, conformément à la loi dites de Snell-Descartes, le rapport sin(0) /c entre : - le sinus de l'angle d'inclinaison e, formé entre la verticale et la direction de propagation de l'onde acoustique au point considéré, et 30 - la vitesse de propagation locale c des ondes acoustiques en ce point, présente la même valeur en chaque point du chemin 31 (figure 5) suivi par cette onde acoustique dans le milieu subaquatique E (ainsi, la quantité sin(0) /c , constante, est égale à sin(00)/c, quel que soit le point considéré le long 16 de la trajectoire 31 suivie par cette onde). Quoi qu'il en soit, l'unité de traitement 18 du système sonar 10 est programmée pour, après chaque émission d'un faisceau d'ondes acoustiques tel que celui décrit plus haut : 5 - acquérir les différents signaux de réception s, captés par les transducteurs 16 de la deuxième antenne 14 en réponse à l'émission de ce faisceau d'ondes acoustiques, - déterminer, en fonction de ces signaux de réception s, qu'une onde acoustique ayant une constante de propagation a donnée a été reçue à un instant 10 de réception, puis - déterminer une durée T entre un instant d'émission du faisceau d'ondes acoustiques, et l'instant de réception de cette onde acoustique (par la deuxième antenne 14). Comme on le verra dans la suite, le couple de données qui comprend 15 cette durée T et la constante de propagation a de l'onde reçue et permet à l'unité de traitement 18 de déterminer la profondeur et la position latérale de l'élément du milieu subaquatique E ayant réfléchit ladite onde acoustique (lors la propagation de cette onde depuis la première antenne 13 jusqu'à la deuxième antenne 14). Les caractéristiques de fonctionnement des première et deuxième 20 antennes 13 et 14 ont été présentées ci-dessus dans un cas où la première antenne 13 fonctionne en émission tandis que la deuxième antenne 14 fonctionne en réception. Les rôles respectifs de ces deux antennes 13, 14 peuvent toutefois être permutés, la première antenne 13 jouant alors le rôle d'antenne réceptrice, et la deuxième antenne 14 jouant le rôle d'antenne émettrice. 25 Pour ce qui est maintenant des unités de pilotage 17 et de traitement 18 du système sonar 10, elles sont réalisées chacune au moyen d'un circuit électronique comprenant au moins un processeur et une mémoire, et, de manière remarquable, elles sont programmées pour exécuter un ou plusieurs des procédés décrits ci-dessous. Le système comprend également une antenne GPS (selon 30 l'acronyme anglo-saxon de « Global Positioning system ») qui fournit une position absolue du système dans le référentiel terrestre, et une centrale d'attitude comprenant au moins un gyromètre, qui délivre des signaux représentatifs des mouvements angulaires du système dans ce référentiel (tangage et roulis notamment). Les mouvements angulaires du bateau sont alors compensés à 17 l'émission et à la réception des ondes, en modifiant électroniquement les directions d'émission et de réception en fonction des signaux délivrés par la centrale d'attitude. Dans la suite de l'exposé on a considéré que ces mouvements étaient parfaitement compensés. 5 Premier procédé 1 détermination d'un profil bathymétrique sur la base d'un profil de célérité moyenne du son. Comme déjà indiqué, au cours de ce premier procédé, la topographie du fond marin situé sous le navire 1 est déterminée, par mesures de temps de propagation d'ondes acoustiques, sur la base d'un profil de célérité moyenne du 10 son dans le milieu subaquatique, cmoy(z). Au cours de ce premier procédé, les transducteurs 15 de la première antenne 13 fonctionnent en émission (cette antenne est employée en tant qu'antenne émettrice). Les transducteurs 16 de la deuxième antenne 14 fonctionnent quant à eux en réception (la deuxième antenne 14 est employée en 15 tant qu'antenne réceptrice). Ce procédé comprend les étapes suivantes (figure 6) : a) émission d'un train d'ondes acoustiques par la première antenne 13, (l'extension temporelle de ce train d'onde étant par exemple de l'ordre d'une milliseconde ou de quelques millisecondes), sous la forme du faisceau d'ondes 20 acoustiques 30 qui a été décrit plus haut, puis b) acquisition des signaux de réceptions s, captés par les transducteurs 16 de la deuxième antenne 14 en réponse à l'émission de l'étape a). Les signaux de réception ainsi acquis permettent à l'unité de traitement de déterminer les profondeurs et positions latérales de plusieurs points P, P', P" 25 du fond 5 marin situés le long d'une ligne de mesure, perpendiculaire à la première antenne 13 (et donc perpendiculaire ici à l'axe longitudinal x du navire 1) et située à l'aplomb du navire 1 (cette ligne de mesure correspond à la ligne moyenne le long de laquelle s'étend la mince bande 32 du fond 5 marin atteinte par les ondes 30 émises à l'étape a)). 30 Les étapes a) et b) sont suivies par une étape c), au cours de laquelle le navire se déplace légèrement, parallèlement à son axe longitudinal x, jusqu'à une nouvelle position. Les étapes a) et b) sont ensuite exécutées à nouveau, ce qui permet de déterminer les profondeurs et positions latérales de points du fond situés le long 18 d'une autre ligne de mesure, légèrement décalée par rapport à la précédente le long de l'axe longitudinal x. Exécuter ainsi l'ensemble des étapes a), b) et c) plusieurs fois successivement permet de relever, ligne par ligne, les positions de différents 5 points d'une surface tridimensionnelle, représentative de la topographie du fond marin considéré. Les signaux de réceptions s acquis lors de chaque exécution de l'étape b) sont traités lors d'étapes de traitement b') et b"). Au cours de l'étape b') l'unité de traitement 18 détermine, en fonction 10 des signaux de réception s acquis à l'étape b), plusieurs couples de données (a, T) associés chacun à l'une des ondes acoustiques reçues par la deuxième antenne 14 à l'étape b). Chacun de ces couples de données comprend la constante de propagation a de cette onde acoustique, ainsi que la durée T séparant l'instant d'émission et l'instant de réception de cette onde. 15 Au cours de l'étape b"), l'unité de traitement 18 détermine, pour chacun de ces couples de données a et T, une profondeur zp et une position latérale yp de la petite portion 33 du fond 5 marin ayant réfléchit l'onde acoustique considérée, lors de sa propagation (aller et retour) depuis la première antenne 13 d'émission jusqu'à la deuxième antenne 14 de réception. Cette profondeur zp et cette position 20 latérale yp repèrent plus précisément le point P du fond sur lequel est centrée cette petite portion 33 de fond marin. La position latérale yp du point P est la coordonnée de ce point le long de l'axe transversal y (axe longitudinal y dont l'origine est située au centre O de la tête sonar 11). Cette position latérale yp est donc égale à l'écart latéral, par 25 rapport à la verticale, entre le point P considéré et la deuxième antenne 14 (antenne réceptrice). D'autre part, c'est en fait la différence de profondeur Az entre le point P et la tête sonar 11 qui est déterminée par l'unité de traitement 18, sur la base de la durée T. Mais la tête sonar 11 est située ici quasiment en surface du milieu 30 subaquatique E, et, pour la simplicité de l'exposé, la différence de profondeur Az est considérée ici comme étant égale à la profondeur zp du point sondé. Cette profondeur zp, parfois appelée immersion, est égale à la distance entre le point P et la surface du milieu subaquatique E, mesurée verticalement. 19 Il est noté par ailleurs que la durée T, entre les instants d'émission et de réception de l'onde acoustique réfléchie par le point P, est égale au temps mis par cette onde acoustique pour, à l'aller, se propager de la première antenne 13 jusqu'au point P, puis, au retour, pour se propager de ce point P jusqu'à la 5 deuxième antenne 14. En l'occurrence, les trajets respectifs suivis par cette onde, respectivement à l'aller et au retour, sont identiques (figure 5). La durée T est ainsi égale au double du temps de propagation t nécessaire pour que cette onde acoustique se propage, depuis le point P considéré jusqu'à la deuxième antenne réceptrice, 14. 10 La méthode mise en oeuvre par l'unité de traitement 18, pour déterminer la profondeur Zp et la position latérale yp du point P sondé est particulièrement remarquable en ce qu'elle permet de tenir compte de la réfraction des ondes acoustiques, causées par les variations de la vitesse de propagation c de ces ondes en fonction de la profondeur z, de manière à la fois simple et précise. 15 Cette méthode est basée sur un modèle de propagation dans lequel il est supposé que l'onde acoustique reçue s'est propagée : - en ligne droite depuis le point P du fond 5 jusqu'à la deuxième antenne 14 (et d'ailleurs aussi en ligne droite depuis la première antenne 13, jusqu'au point P), 20 - et avec une célérité moyenne cmoy(zp) ne dépendant que de la profondeur Zp du point P considéré (cette célérité moyenne étant notamment indépendante de l'inclinaison de la direction de propagation de cette onde par rapport à la verticale). Dans ce modèle, la profondeur Zp et la position latérale yp du point P 25 s'expriment donc conformément aux formules F2 et F3 suivantes : zp = cmoy(zp). t. cos(0',,y) (F2) yp = cmoy(zp). t. sin(O',,y) (F3) où : - le temps de propagation t est égal à la moitié de la durée T séparant les 30 instants d'émission et de réception de l'onde acoustique : t = T/2, et où - emoy est un angle de propagation moyen, défini entre la verticale et une direction de propagation effective de l'onde acoustique (direction de propagation moyenne de l'onde acoustique, entre le point P et la deuxième antenne 14). 20 Autrement formulé, l'angle de propagation moyen Omo), est défini entre l'axe vertical z et la droite 31' qui relie le point P à la deuxième antenne 14. Cette droite 31' correspond au chemin effectif, rectiligne, qui est supposé être suivi par l'onde acoustique dans ce modèle. 5 Moyennant un choix approprié de la célérité moyenne cmoy(zp) et de l'angle de propagation moyen 0',,y, il s'avère que ce modèle permet une détermination très fiable et précise de la profondeur zp et de la position latérale yp du point P ayant réfléchi l'onde reçue. En l'occurrence, l'angle de propagation moyen Omoy est déterminé en 10 fonction de la constante de propagation a de l'onde acoustique reçue (il dépend en effet de l'angle de réception 0, sous lequel cette onde est reçue par la deuxième antenne), et en fonction de ladite célérité moyenne cmoy(zp), ou éventuellement d'une célérité arithmétique moyenne E(zp) de l'onde acoustique sur ce trajet effectif. Cette célérité arithmétique moyenne E(zp) est égale à la 15 moyenne arithmétique d'une pluralité de valeurs de la vitesse de propagation locale c(z), présentées par les ondes acoustiques à une pluralité respective de profondeurs z allant de la profondeur zp du point P jusqu'à la profondeur z=0 de la deuxième antenne 14 (réceptrice) : E(zp) = Zp LzP c(z) dz (F4). 20 Le fait que la célérité moyenne cmoy(zp) soit indépendante de la direction de propagation effective de l'onde acoustique facilite grandement le traitement des données acquises et la détermination de la profondeur et de la position latérale de chaque point du fond sondé par le système sonar. En particulier, la profondeur zp, ainsi que la valeur de la célérité moyenne 25 cmoy(zp) correspondante peuvent être déterminées dans un premier temps, indépendamment de la position latérale yp du point P, ce qui permet de réduire le temps de calcul nécessaire à cette détermination (la position latérale yp étant déterminée ensuite, sur la base de la valeur de la célérité moyenne ainsi déterminée). Cette réduction de temps de calcul est d'autant plus intéressante ici 30 que les profondeurs et positions latérales de plusieurs points P, P', P" du fond sont déterminées, à chaque répétition de l'étape b"). La célérité moyenne cmoy(zp) est égale à la moyenne harmonique d'une pluralité de valeur de la vitesse de propagation locale c(z), présentées par les 21 ondes acoustiques à une pluralité respective de profondeurs z allant de la profondeur zp du point P jusqu'à la profondeur z=0 de la deuxième antenne 14 (réceptrice), soit : 1 1 f ZP 1 cmay(zp) = zp jo c(z)dz Il est noté que la quantité fz0 P c( 1z) dz = fozP dt'' = t'' est égale à la 5 durée t',-, que mettrait une onde acoustique pour se propager, à la verticale depuis la profondeur zp du point sondé, jusqu'à la profondeur z=0 de l'antenne réceptrice. La célérité moyenne cr',,,(zp) est donc égale à la différence de profondeur zp entre le point P sondé et la deuxième antenne 14, divisée par cette durée tvert. 10 Quant à l'angle de propagation moyen Onicy, il est déterminé ici par l'unité de traitement 18 de sorte que son sinus soit égal au produit de la constante de propagation a de l'onde considérée, multipliée par la célérité moyenne cmoy(zp) sin(Bmoy) = a. cmoy(zp) (F6) 15 L'unité de traitement 18 est donc programmée pour déterminer la profondeur zp et la position latérale yp du point P, soit conformément aux formules F2 et F3 ci-dessus (l'angle de de propagation moyen Omoy étant alors déterminé conformément à la formule F6), soit directement au moyen des formules F7 et F8 ci-dessous : 2 20 zp = cmoy(zp). t.,\11 - (a. c''),(zp)) (F7) yp = cmoy(zp). t. (a. cmoy(zp)) (F8) . Dans ce premier procédé, un profil de la célérité moyenne du son, rassemblant plusieurs valeurs présentées par la célérité moyenne cmoy(z) pour une pluralité respective de profondeurs z dans le milieu subaquatique, est pré- 25 enregistré dans la mémoire de l'unité de traitement 18. Sur la base de ce profil de célérité moyenne, préenregistré, et des données trc/2 et a mesurées, l'unité de traitement 18 cherche la profondeur zp qui satisfait la formule F7 (par exemple de manière itérative, comme expliqué plus bas). Une fois cette profondeur déterminée, la position latérale yp du point P est (F5) 22 déterminée, en calculant directement la quantité cmoy(zp).t. (a. cm,,y(zp)) (formule F8). Le profil de la célérité moyenne cmoy(z) est mesuré puis enregistré dans la mémoire de l'unité de traitement 18 lors d'une phase préliminaire de ce 5 procédé. Lors de cette phase préliminaire, une sonde, munie d'un capteur apte à mesurer la vitesse locale de propagation c des ondes acoustiques à la profondeur z de la sonde, est positionnée successivement à plusieurs profondeurs dans le milieu subaquatique E. Un profil c(z) de la vitesse locale de propagation est ainsi relevé. En variante, cette sonde pourrait par exemple être munie de capteurs 10 aptes à déterminer la température et la salinité de l'eau à la profondeur de la sonde, quantités qui permettent ensuite de déterminer la vitesse locale de propagation c(z) à la profondeur considérée. Le profil de la célérité moyenne cmoy(z) est déterminé ensuite, sur la base du profil c(z) de la vitesse locale de propagation relevé précédemment, par 15 intégration numérique, conformément à la formule F5. La figure 7A représente un exemple de profil c(z) de la vitesse locale de propagation (exprimée en mètres par seconde), mesuré pour des profondeurs (exprimées en mètres) comprises entre 0 et 2500 mètres. Comme on peut le constater sur cet exemple, les variations de la vitesse locale de propagation c 20 avec la profondeur z, sur cette gamme de profondeur, peuvent atteindre 20%. Il est donc effectivement souhaitable, comme c'est le cas ici, de tenir compte de ces variations lors de l'analyse des temps de propagation mesurés. La figure 7B montre le profil de la célérité moyenne cmoy(z) (en mètres par seconde) déduit du profil de la vitesse locale de propagation c(z) de la figure 25 7A. Cet exemple illustre bien que les variations de la célérité moyenne cmoy(z) avec la profondeur z (exprimées en mètre) sont plus progressives que les variations de la vitesse locale de propagation c(z) avec la profondeur. En effet, l'opération d'intégration numérique correspondant à la formule F5 a un effet de filtre passe-bas. Elle élimine certaines variations rapides de la vitesse de 30 propagation c(z), ainsi que différents bruits affectant les mesures de cette vitesse. Le fait de déterminer les profondeurs des points sondés sur la base du profil de célérité moyenne cmoy(z) plutôt que sur la base du profil de vitesse locale de propagation du son c(z) permet donc de s'affranchir de ces bruits, et d'améliorer la 23 précision de détermination de ces profondeurs. Un exemple de méthode itérative, permettant de déterminer la profondeur zp du point P sondé, de manière à ce qu'elle satisfasse la formule F7, est décrite maintenant. 5 Cette méthode comprend une étape de détermination d'une estimation zp,; de la profondeur zp, et d'une estimation cmoy.i de la célérité moyenne correspondante cmoy(zp). Cette étape d'estimation est exécutée plusieurs fois successivement, de manière itérative. Lors de l'exécution numéro i de cette étape, une nouvelle estimation 10 cmoy,i de la célérité moyenne est déterminée, sur la base d'un estimation précédente de la profondeur zp, et sur la base du profil de célérité moyenne préenregistré, conformément à la formule F9 suivante : Crnoy,i = Crnoy(ZP,i-1) (F9). Une nouvelle estimation zp,; de la profondeur zp est ensuite déterminée 15 en remplaçant la célérité moyenne cmoy(zp) par son estimation cmym dans la formule F7, soit : \ 2 Zp j = Cmoy,E. t. j 1 - (a. Cmoy.i)2 (F10) Cette méthode itérative débute par exemple à partir de l'estimation initiale cmoy,i = cm,y(z=0) de la célérité moyenne. 20 Les exécutions successives de l'étape d'estimation en question cessent lorsqu'un écart lcm,y,i+i - cmoy,i1 entre deux estimations successives cmoy,i-vi et cmoy,i de la célérité moyenne devient inférieur à un seuil donné. En pratique, ce seuil peut par exemple être compris entre 1 mètre par seconde et 1 millimètre par seconde. En variante, les exécutions successives de cette étape d'estimation 25 pourraient être arrêtées lorsqu'un écart relatif entre deux estimations successives zp,i+i et zp,; de la profondeur zp devient inférieur à une précision requise donnée. La figure 8 représente schématiquement une erreur cz, représentative d'un écart entre la profondeur réelle d'un point sondé P, et la profondeur zp déterminée sur la base du modèle monocouche du milieu subaquatique E 30 présenté ci-dessus (modèle dans lequel ce milieu est représenté par une couche homogène effective, dans laquelle la propagation est rectiligne). Cette erreur cz est représentée en fonction du sinus de l'angle de réception 0,, pour une profondeur du point sondé de 200 mètres. 24 Une erreur sy, représentative d'un écart entre la position latérale réelle d'un point sondé P, et sa position latérale yp déterminée sur la base de ce modèle, est également représentée sur la figure 8. Les deux erreurs a, et ey représentées sur cette figure ont été 5 déterminées par simulation numérique, sur la base du profil de vitesse de propagation c(z) de la figure 7A. Pour cela, les trajets respectifs suivis par des ondes acoustiques, émises avec directions de propagation initiales différentes, ont été déterminés en tenant compte de la réfraction de proche en proche (couche par couche), jusqu'à une profondeur de 200 mètres. Cette simulation a permis de 10 déterminer le temps de propagation t, qui aurait mesuré en pratique pour une onde suivant un tel trajet, puis d'en déduire la profondeur zp et la position latérales yp, qui aurait été déterminées en pratique par l'unité de traitement sur la base de ce temps de propagation t. L'erreur cz est égale à l'écart entre la profondeur zp, issue de cette simulation d'une mesure de profondeur, et la profondeur de 15 référence employée, de 200 mètres. L'erreur ey est égale quant à elle à l'écart entre la position latérale yp, telle qu'elle aurait été déterminée par l'unité de traitement, et la position latérale issue du tracé du trajet de l'onde acoustique considéré, réalisé de proche en proche, couche par couche. Cette simulation a été réalisée pour des valeurs de l'angle de réception Of, comprises entre 0 et 60 20 degrés (et donc pour des valeurs de sin(610) comprises entre 0 et 0,85 environ). Comme on peut le voir sur cette figure, même lorsque les ondes acoustiques reçues sont très inclinées par rapport à la verticale (situation dans laquelle leurs trajectoires sont fortement courbées par réfraction), les erreurs cz et ey restent extrêmement petites : ces erreurs restent inférieures à 5 centimètres, 25 alors que la profondeur sondée est de 200 mètres. Pour des angles de réception inférieurs à 30 degrés, ces erreurs sont même inférieures à 3 millimètres. On constate ainsi, sur cet exemple, que la méthode de détermination mise en oeuvre par l'unité de traitement 18 conduit à des valeurs très fiables et précises de profondeur et de position latérale, tout en simplifiant considérablement 30 la détermination numérique de ces grandeurs par rapport à une détermination qui serait basée sur un tracé point par point du trajet suivi par l'onde. Le fait que ce modèle monocouche conduise à des résultats précis peut s'expliquer en partie sur la base des arguments théoriques présentés ci-dessous. Tout d'abord, il est connu dans le domaine des mesures sismiques que 25 le temps de propagation t d'une onde sismique entre une point source et un récepteur est donné, en bonne approximation, par la formule F11 ci-dessous : t2 = tvert Y2/dms (F11) où : - y désigne l'écart latéral, par rapport à la verticale, entre le point source 5 considéré et le récepteur, - tvert est le temps que mettrait une onde sismique pour se propager, verticalement, depuis la profondeur z du point source considéré jusqu'à celle du récepteur (grandeur qui a déjà été mentionnée plus haut, lors de la présentation de la formule F5), et où 10 - cRMS est une célérité quadratique moyenne, le long du trajet suivi par l'onde sismique : dms =10 ft c2(t) dt t Cette formule (mentionnée par exemple dans le chapitre 5 de l'ouvrage « Geophysical Signal Processing » de E.A. Robinson, T.S. Durrani et L.G. Peardon, Prentice Hall, 1986, ISBN 0133526674), est utilisée habituellement dans 15 le domaine des mesures sismiques pour aligner temporellement différents signaux reçus, afin de pouvoir ensuite sommer ces signaux pour améliorer le rapport signal sur bruit par effet de moyenne. En l'occurrence, le temps de propagation « à la verticale », t',-t, est égal à la profondeur z, divisée par la moyenne harmonique cH(z) de la vitesse de 20 propagation locale c, entre cette profondeur z et la surface : z_ n rz 1 rin, tvert - (z) L - Jo c(zo Par ailleurs, dans un milieu subaquatique, à la différence d'un milieu géologique, la célérité harmonique moyenne cH est très proche de la célérité quadratique moyenne cRMS. Par exemple, dans le cas du milieu subaquatique E 25 considéré ici, dont le profil de vitesse de propagation est tracé sur la figure 7A, la différence entre la célérité harmonique moyenne chl et la célérité quadratique moyenne cRMS reste inférieure à 2 centimètres par seconde, pour toute la plage de profondeurs considérée. Il est noté d'ailleurs que, dans ce milieu subaquatique, la célérité harmonique moyenne cH et la célérité arithmétique moyenne C. sont elles 30 aussi très proches l'une de l'autre, puisqu'elles sont reliées à la célérité quadratique moyenne CRms par la relation suivante : cRMS = cH. C (la célérité quadratique moyenne est égale à la moyenne géométrique de célérité harmonique 26 moyenne, et de la célérité arithmétique moyenne). Pour des mesures réalisées en milieu liquide, on a donc, en bonne approximation : t2 = (z2 y2)/d 5 Or, dans le modèle de propagation monocouche qui a été présenté plus haut, le temps de propagation t est relié à la profondeur zp et la position latérale yp du point sondé par la relation suivante : t2 = (4 +34)/cLy où c - nmy est la célérité moyenne, effective, dans l'unique couche 10 considérée. Pour que ce modèle de propagation soit le plus représentatif possible du lien physique réel entre le temps de propagation t et les coordonnées zp, yp, il doit donc effectivement être paramétré en choisissant, pour la célérité moyenne cmoy(Zp), la célérité harmonique cH mentionnée ci-dessus, soit, comme déjà 15 indiqué : 1 Crrapy(Zp) = CH(Zp) = \ fz-Pd.z) . Z p C(Z) Par ailleurs, en étendant et modifiant certains résultats présentés dans le document « Approximate Methods For Ray Tracing », M.J. Daintith, SaclantCen Conference proceedings N°5, Septembre 1971, l'inventeur a démontré les 20 résultats suivants : un développement des équations de propagation au premier ordre en E(z), où c(z) = c (1 + e(z)), conduit aux formules suivantes : zp = 1 - {i -2 2 3±}) et yp = aE2t (1 + E2 (2A2 - 1)) avec aE = sin OniGy, A = tan Omoy, , E2 (2) = zp (i) fzP e2(z)dz 25 En négligeant alors les termes en A2 e2 et A' e2 dans ces équations, et en utilisant le fait que, au premier ordre, la célérité harmonique moyenne s'écrit cH = 41- e2 on obtient les formules d'approximations suivantes : zp = (cH t),/1 a2E2, yp = (cH ctE et D= J4+yP=cNt. Ces dernières équations montrent que l'angle de propagation moyen à 30 employer est, en très bonne approximation défini par la formule F12 suivante : sin(Omoy) = a. E(zp) (F12). 27 En approximant la célérité arithmétique moyenne E par la célérité harmonique moyenne cH dans la formule F12, on obtient alors la relation suivante : sin(Omoy) = a.cH sur la base de laquelle l'unité de traitement 18 détermine la profondeur et 5 la position latérale des points sondés. La validité de l'approximation sin(Omoy) = a. ci', montrée ci-dessus d'un point de vue théorique, est justifiée en outre, a posteriori, par la précision des résultats qu'elle permet d'obtenir (C.f. : figure 8). Les calculs ci-dessus montrent toutefois que l'approximation de la 10 formule F12, sin(Omoy) = a. E(z p), est encore meilleure. En variante, l'unité de traitement 18 pourrait donc être programmée de manière à déterminer la profondeur zp et la position latérale yp de chaque point P sondé conformément aux formules F13 et F14 suivantes : zp = cmoy(z p). t - (a. C(zp))2 (F13) 15 yp = c,,,,y(zp). t. (a . E(z p)) (F14). Dans le premier procédé qui vient d'être décrit, les positions de différents points d'une surface tridimensionnelle, représentative de la topographie du fond 5 marin, est relevé ligne par ligne, grâce au déplacement longitudinal du navire 1. En variante, le navire pourrait conserver la même position tout au long de 20 ce procédé, la surface tridimensionnelle en question étant alors relevée en faisant varier l'inclinaison du faisceau d'ondes acoustiques 30 émis, entre une émission et la suivante (c'est-à-dire entre deux exécutions successives de l'étape a)), au lieu de systématiquement émettre ce faisceau à la verticale. En variante encore, au lieu de relever les positions de plusieurs points du 25 fond couvrant toute une région de celui-ci, on pourrait prévoir de ne relever la profondeur et/ou la position latérale que d'un point du fond. D'autre part, ce premier procédé peut être employé pour déterminer les profondeurs et/ou positions latérales d'éléments immergés autres que des portions du fond marin, tels qu'un poisson, un banc de poissons, ou une partie 30 d'un banc de poissons (ce procédé n'est donc pas seulement un procédé de détermination d'un profil bathymétrique). Les profondeurs et/ou positions latérales de tels éléments sont déterminées de la même manière que ce qui a été expliqué 28 ci-dessus pour un point du fond marin. Deuxième procédé : détermination d'un profil bathymétrique, et détermination d'un profil de célérité moyenne du son sur la base des relevés bathymétriques eux-mêmes 5 Dans ce deuxième procédé, comme dans le premier procédé, les profondeurs de différents points du fond marin sont déterminées, à partir de mesures de temps de propagation, sur la base du modèle monocouche (modèle à une couche effective) qui a été décrit plus haut. Mais ici, à la différence du premier procédé, le profil de célérité moyenne 10 cmoy(z) dans le milieu subaquatique E, qui intervient lors de l'analyse des mesures de temps de propagation, n'est pas préenregistré dans une mémoire du système sonar 10 (il n'est pas mesuré préalablement aux relevés de profondeurs). Dans ce deuxième procédé, la profondeur d'un point donné du fond marin est déterminée, sans connaissance préalable du profil de célérité moyenne, 15 grâce à une comparaison d'au moins deux mesures distinctes de temps de propagation entre le système sonar 10 et le point sondé, ces deux mesures étant réalisées pour deux positions distinctes du système sonar 10 (et donc pour deux inclinaisons différentes des ondes acoustiques reçues). La comparaison de ces deux temps de propagation permet en outre de 20 déterminer la valeur de la célérité moyenne, cmoy(zm), à la profondeur zm du point M considéré. Les figures 9 et 10 représentent schématiquement la façon dont est réalisée une série de mesures de temps de propagation permettant, par comparaison de ces mesures les unes avec les autres, de déterminer : 25 - les profondeurs respectives de plusieurs points du fond, M, M', M", situés ici le long d'une ligne de mesure parallèle à l'axe longitudinal x du navire 1, et - les valeurs correspondantes de la célérité moyenne cmy dans ce milieu. Au cours de cette série de mesures, pour une première position 01 du 30 navire 1, la première antennel3 émet un faisceau transversal 40 d'ondes acoustiques, tel que le faisceau 30 décrit plus haut en référence à la figure 3. Ce faisceau est émis ici à la verticale, sous le système sonar 10. Le plan moyen défini par ce faisceau, parallèle à l'axe transversal y, s'étend à la verticale sous le système sonar 10. Ce faisceau atteint une mince bande 42 du fond marin, 29 perpendiculaire à l'axe longitudinal x. En réponse à cette émission, la deuxième antenne 14 reçoit les ondes acoustiques qui ont été réfléchies par les éléments de ce bande 42, et les transducteurs 16 de la deuxième antenne 14 acquièrent les signaux de réception correspondant. Sur la base de ces signaux, l'unité de 5 traitement 18 détermine notamment un premier temps de propagation tut mis par l'une de ces ondes pour se propager depuis le système sonar 11 jusqu'à un point M du fond situé à l'aplomb de ce système (en l'occurrence, cette onde acoustique se propage selon un trajet 41 vertical et rectiligne, avec une première constante de propagation ctiel qui est donc nulle). 10 Le navire 1 se déplace ensuite, parallèlement à son axe longitudinal x, jusqu'à d'autres positions intermédiaires. En chacune de ces positions, le système sonar 11 détermine un temps de propagation ... entre le système sonar 11 et un point M', M", ... du fond situé à l'aplomb de ce système, comme expliqué ci-dessus pour la première position 01. 15 Le navire 1 se déplace ensuite jusqu'à une deuxième position 02. A cette deuxième position, la deuxième antenne 14 émet un faisceau longitudinal 50 d'ondes acoustiques qui se propage sous la forme d'une nappe peu épaisse dont le plan moyen est parallèle à l'axe longitudinal x (puisque la deuxième antenne 14 s'étend perpendiculairement à cet axe). Ce faisceau atteint ainsi une mince bande 20 52 du fond marin, qui cette fois est parallèle à l'axe longitudinal x (au lieu de lui être perpendiculaire). En réponse à cette émission, la première antenne 13 reçoit les ondes acoustiques qui ont été réfléchies par les éléments de cette bande 52, et les transducteurs 15 de la première antenne 13 acquièrent les signaux de réception correspondant. Sur la base de ces signaux, l'unité de traitement 18 25 détermine alors - des deuxième temps de propagation tm.2, tM',2, tM-.2 , , mis par les différentes ondes acoustiques reçues, pour se propager depuis le système sonar 11 jusqu'aux points M, M', M" du fond sondés grâce à l'émission de ce faisceau longitudinal 50 (points qui sont réparties le long de la ligne de mesure, parallèle à 30 l'axe longitudinal x, mentionnée ci-dessus), et - des deuxième constantes de propagation OEM.2, M'.2, aM".2 , qui sont les constantes de propagation de ces ondes acoustiques. La profondeur de chaque point M, M', M" situé le long de cette ligne de mesure est ainsi sondé : 30 - une première fois (depuis la deuxième position 01), au moyen d'une onde acoustique se propageant à la verticale, puis - une deuxième fois (depuis la deuxième position 02), au moyen d'une onde acoustique dont la direction de propagation est inclinée par rapport à la 5 verticale (lors de cette deuxième mesure, le trajet 51 suivi par l'onde acoustique qui sonde la position du point M est donc courbé par la réfraction, comme représenté sur la figure 10). La profondeur de ce point est alors déterminée par l'unité de traitement 18 en comparant cette première mesure et cette deuxième mesure l'une avec 10 l'autre. Pour cela, l'unité de traitement 18 détermine tout d'abord la valeur de la célérité moyenne, cmoy(zm), à la profondeur zm du point M considéré, en fonction du premier temps de propagation tM,1, du deuxième temps de propagation tM,2, et de la deuxième constante de propagation am,2. 15 L'unité de traitement 18 détermine ensuite la profondeur zm du point M, sur la base de cette valeur de la célérité moyenne. Plus précisément, l'unité de traitement 18 détermine la valeur de la célérité moyenne cmoy(zm) comme étant la valeur pour laquelle un écart quadratique £2 entre : 20 - une première estimation zm,i de la profondeur zm, déterminée sur la base du premier temps de propagation tM,1, et - une deuxième estimation zm,2 de la profondeur zm, déterminée sur la base du deuxième temps de propagation tM,2 et de la deuxième constante de propagation am.2, 25 est minimum (ou tout au moins inférieur à un seuil, par exemple inférieur à 10 centimètres), chacune de ces deux estimations étant déterminée sur la base de cette même valeur de la célérité moyenne cmoy(zM). Ainsi : Cmoy(%) = argmin (E2) = argmin en, Oz"( tua, cm) - zm,2(tm,2, am,2, cm)12) (F15). La première estimation ZM.1 et la deuxième estimation zm,2 sont 30 déterminées chacune, sur la base du modèle monocouche présenté plus haut, conformément aux formules F2 et F6, ou conformément à la formule F7, soit : 31 ,\I -,2 Zma = Cm. tm,i et zm,2 = Cm. tM,2- 1 - (Ctm,2. Cm ) . Dans le mode de réalisation particulier du deuxième procédé qui est décrit ici, la minimisation de l'écart quadratique £2 (qui revient ici à annuler cet écart) conduit alors à déterminer la célérité moyenne cmoy(zM) au point M 5 conformément à la formule F16 suivante : Cmoy(ZM) = 1 (F16). aM,2 L'unité de traitement 18 peut par exemple déterminer la célérité moyenne cmoy(zM) conformément à la formule F16, ou par une recherche numérique de la valeur de cm qui minimise l'écart quadratique E2 (conformément à la formule F15). 10 Il est noté que, dans ce procédé, parmi les différents couples de données (am.2, tA1.2), (aM',2, tN1',2), (aM",2, , associés respectivement aux différentes ondes reçues par la première antenne 13 depuis la deuxième position 02, il faut identifier celui, (am,2, tm,2), correspondant à l'onde acoustique qui a été réfléchie par le même point M que celui sondé depuis la première position 01 (puisque 15 cette méthode est basée sur la comparaison de deux estimations distinctes de la profondeur d'un même point). Pour cela, l'unité de traitement 18 peut par exemple sélectionner celui de ces couples de données pour lequel un écart latéral xM, déterminé en fonction du couple de données (a, t) considéré (par exemple selon la formule xm = cloy. t. a), est le plus proche de la distance séparant la première 20 position 01 de la deuxième position 02 (distance qui est connue par ailleurs, par exemple grâce à un système de localisation 2, tel qu'un système de type GPS ou un loch équipant le navire 1). La figure 11 représente schématiquement la valeur de l'écart c entre les deux estimations de profondeur zm,i et zm,2 (exprimé en mètre), en fonction de la 25 variable à ajuster cm (en mètres par seconde), dans une situation où le fond marin est plat, situé à une profondeur de 200 mètres. Les valeurs des premier et deuxième temps de propagation et de la deuxième constante de propagation, dont sont issues ces estimations de profondeurs zm,i et zm,2, sont obtenues ici par simulation numérique, comme expliqué plus haut à propos de la figure 8. 30 Comme on peut le voir sur cette figure, dans cette situation, la plus petite valeur de l'écart s est obtenue pour cm = 1510,8 mètres par seconde. Or, dans le 32 milieu aquatique E considéré, la valeur de la célérité moyenne cmoy, à une profondeur de 200 mètres, est de 1510,7 mètres par seconde. Cela montre que cette méthode permet effectivement, par comparaison de plusieurs mesures bathymétriques, de déterminer une valeur de la célérité moyenne extrêmement 5 proche de la valeur réelle de cette célérité moyenne dans le milieu E sondé. La figure 12 montre d'ailleurs plusieurs valeurs de la célérité moyenne cmoy déterminées à partir des relevés bathymétriques eux-mêmes (ces valeurs, exprimées en mètres par seconde, sont reportées sur ce graphique sous la forme d'étoiles), pour différentes profondeurs z (exprimées en mètres) comprises entre 10 10 et 150 mètres, et, pour les mêmes profondeurs, des valeurs de référence de la célérité moyenne, crfloy,ref (représentées sous la forme d'une ligne continue). Les valeurs de la célérité moyenne, issues des relevés bathymétriques eux-mêmes, sont déterminées à partir de temps de propagation obtenus comme précédemment par simulation numérique, dans une situation où la profondeur du 15 fond marin varie progressivement, de 10 mètres jusqu'à 150 mètres, sur une longueur de 3 kilomètres. Cette figure montre à nouveau que les valeurs de célérité moyenne déduites des relevés bathymétriques sont très proches de valeurs de référence de cette célérité (ces valeurs de référence cmoy.ref étant obtenues par intégration d'un profil de référence de la vitesse de propagation 20 locale cref(z)). On peut prévoir en outre, au cours de ce deuxième procédé, que l'unité de traitement 18, détermine un profil de vitesse de propagation locale des ondes acoustiques, c(z), à partir du profil de célérité moyenne cmoy(z) déterminé comme expliqué ci-dessus (c'est-à-dire à partir de plusieurs relevés bathymétriques 25 comparés les uns aux autres). Pour cela, la vitesse de propagation locale c(z), à la profondeur z, peut par exemple être déterminée conformément à la formule d'inversion F17 suivante (déduite de la formule F5) : ci2,,'y(z) c(z) = (F17) cmoy(z) - z. dc,floy(z)/dz La figure 13 représente schématiquement un profil de vitesse de 30 propagation locale c(z) (exprimée en mètres par seconde), déterminé de cette manière, pour différentes profondeurs z (exprimées en mètres) comprises entre 10 et 150 mètres. Pour ces mêmes profondeurs, cette figure représente aussi un 33 profil de référence de la vitesse de propagation locale c'f(z) (ce profil de référence rassemble les valeurs de la vitesse de propagation locale, telles qu'elles auraient été mesurées par une sonde placée successivement à différentes profondeurs, par exemple ; ce profil n'est connu ici que pour quelques profondeurs réparties 5 dans la gamme allant de 10 à 150mètres, et présente donc une allure de ligne brisée, avec une dérivée constante par morceaux). Ces deux profils sont très proches l'un de l'autre (l'écart relatif entre ces deux profils reste inférieur à 0,5 pour mille), ce qui montre que cette méthode de détermination du profil de vitesse de propagation locale c(z), par comparaison de mesures bathymétriques, conduit 10 à des résultats fiables et précis. Il est noté toutefois que dans la formule 17, la dérivée dc',Q),(z)/dz présente au dénominateur peut avoir pour effet d'amplifier de légers bruits ou erreurs numériques affectant éventuellement le profil de la célérité moyenne cmy(z). 15 Pour y remédier, l'unité de traitement 18 peut, par exemple, mettre en oeuvre une procédure de régularisation telle que celles décrites ci-dessous. Par exemple, en supposant une vitesse de propagation locale constante dans l'intervalle [zi, z2], notée c, les mesures de la célérité moyenne cmoy(zi) en z1, et cm,y(z2) en z2 donnent : 1 1 z2 z1 20 C Z2 (Cyno/Z2) Cmoy (Zi) Plus généralement, on peut supposer que le profil de vitesse de propagation locale c(z) varie linéairement sur l'intervalle de profondeurs [zi, considéré : c(z) = c1 + g(z - z1), pour z E [z1, z2]. Les valeurs des deux variables cl et g, qui caractérisent ce profil de 25 vitesse, sont alors déterminées par l'unité de traitement comme étant les valeurs minimisant la fonction coût FC suivante : FC =XKz, g ccni'y(zt) cinozy(z1)) g1lo ( 1 + (cz: z1)) où la somme Z porte sur les différentes profondeurs z, de I intervalle [zi, z2] pour lesquelles la célérité moyenne criwy(zt) a été déterminée. Différentes variantes peuvent être apportées au deuxième procédé qui 30 vient d'être décrit. 34 Tout d'abord, la mesure du premier temps de propagation tAq.1 pourrait être réalisée au moyen d'une onde dont la direction de propagation soit inclinée par rapport à la verticale (au lieu de se propager verticalement, comme c'est le cas dans l'exemple de réalisation décrit ci-dessus). Dans ce dernier cas, la valeur 5 de la célérité moyenne cm''(zm) au point M est déterminée en tenant compte d'une première constante de propagation aml, avec laquelle cette première onde se propage, en plus de prendre en compte les premier et deuxième temps de propagation tut 42 et la deuxième constante de propagation am,2. Dans le cadre de cette variante, la première estimation zmi, de la profondeur zm est par exemple 10 déterminée comme étant égale à zma = cm. tel. 1 - cm )2 (au lieu d'être déterminée en fonction de cm et tmi seulement). Pour ce qui est de la valeur de la célérité moyenne cmoy(zm), elle est déterminée, là encore, de manière à minimiser l'écart quadratique £2 (ou tout au moins de manière à le rendre inférieur à un seuil donné). 15 En variante encore, on pourrait prévoir, dans ce deuxième procédé, de mesurer des temps de propagation tm entre le système sonar et un même point M du fond marin, pour plus de deux positions distinctes du système sonar. Pour chacune de ces positions, Oh le système sonar détermine alors la constante de propagation am de l'onde sonore ayant sondé la profondeur du point M, en plus 20 du temps de propagation tmi précité. La valeur de la célérité moyenne cmoy(zm) au point M est alors déterminée par l'unité de traitement 18 comme étant la valeur pour laquelle un écart global c', représentatif de la dispersion d'un ensemble d'estimations zmi de la profondeur zm, déterminées chacune sur la base de l'un des couples (ami, tm,i) 25 mesuré depuis l'une desdites positions 0i, est minimal (ou inférieur au seuil précité). L'écart global £' peut par exemple être égal à l'écart-type ou à la variance de cet ensemble d'estimations, ou à la somme des écarts quadratiques entre ces estimations prises deux par deux. Chacune de ces estimations est déterminée, 30 comme précédemment, sur la base du modèle monocouche décrit plus haut (et donc, par exemple, conformément à la formule F7). A titre d'exemple, si trois mesures de temps de propagation tm,i, tm,2, tm.3 sont réalisées pour repérer le point M, pour trois positions distinctes du système sonar, cet écart global s' peut être défini comme suit : F e' = kma(tma, am,i,cm)- Zfig.2 am,2, Cm)]2 [Zm,2 (tm,2, am,2,cm)- Zm,3 (tm,3, am3, 42)12 + [Zm,3 Cm) - zma cm/J2 . Troisième procédé : détermination de la position d'une balise immergée 5 Dans ce troisième procédé, un émetteur 4 d'une balise 3 immergée émet une onde acoustique, qui est ensuite reçue par le système sonar 10 (figures 14 et 15). L'unité de traitement 18 du système sonar 10 détermine alors une position de cette balise 3, par rapport au système sonar 10, sur la base : 10 - de la constante de propagation a de cette onde, mesurée lors de sa réception par le système sonar, et, selon le mode de réalisation considéré, sur la base en outre - d'un temps de propagation tB, mis par cette onde pour se propager depuis cet émetteur 4 jusqu'au système sonar 10. 15 Il est noté que l'onde acoustique en question effectue ici un trajet en aller simple, alors que les ondes acoustiques émises lors des premier et deuxième procédés effectuent des trajets aller puis retour (à la manière d'un écho). Le temps de propagation tB est donc égal à la durée séparant l'instant d'émission de cette onde, et l'instant de réception de cette onde (au lieu d'être égal à la moitié d'une 20 telle durée). Quoiqu'il en soit, la position de cette balise 3 est déterminée là encore sur la base du modèle de propagation monocouche qui a été présenté plus haut, en fonction du profil de la célérité moyenne cmo(z). Dans ce troisième procédé, ce profil est préenregistré dans la mémoire de l'unité de traitement 18. 25 D'autre part, dans ce troisième procédé, les première et deuxième antennes 13 et 14 sont employées toutes deux en réception. Le chemin 61 suivi par l'onde acoustique émise, entre la balise 3 et le système sonar 10, est situé dans un plan de propagation PI vertical qui passe par le centre B de l'émetteur 4 et par le centre O de la tête sonar 11. 30 La direction de propagation de l'onde acoustique, au point O de réception de cette onde, est repéré par : - son angle de réception 60, formé, dans le plan de propagation PI, entre cette direction de propagation et la verticale z, et par 36 - un angle de pointé 9, formé entre le plan de propagation PI et le plan (z,x) (plan vertical parallèle à l'axe longitudinal x). La constante de propagation a de l'onde acoustique est, comme précédemment, égale au sinus de l'angle de réception 0,), divisé par la vitesse de 5 propagation locale c, des ondes acoustiques à la profondeur des antennes réceptrices, 13 et 14. L'unité de traitement 18 est programmée pour déterminer notamment la constante de propagation a, sur la base des signaux de réception captés par les transducteurs 15, 16 de la première antenne 13 et de la deuxième antenne 14 lors de la réception de cette onde acoustique. 10 En l'occurrence, le décalage temporel Atx entre deux signaux de réception délivrés par deux transducteurs 15 de la première antenne 13 distants de d, est donné par at, = d,. a. cos(w). De manière comparable, le décalage temporel Aty entre deux signaux de réception délivrés par deux transducteurs 16 de la deuxième antenne 14 distants de d1 est donné par aty = d y . a. sin(w). La 15 constante de propagation a de l'onde acoustique reçue, et les quantités a.cos(p) et a.sin(p) peuvent alors être obtenues sur la base de ces décalages temporels. La constante de propagation est déterminée par exemple en calculant la quantité \I(Atx/d,32 + (Aty/dy)2 . Dans un premier mode de réalisation de ce troisième procédé, l'unité de 20 traitement 18 détermine la position du centre B de l'émetteur 4, par rapport au système sonar, sur la base à la fois de la constante de propagation a, et du temps de propagation tB. Ce premier mode de réalisation est adapté à une situation dans laquelle la balise 3 et le système sonar 10 partagent une même référence temporelle 25 commune, ce qui permet à l'unité de traitement 18 de déterminer le temps de propagation tB (le système sonar et la balise disposant par exemple chacun d'une d'horloge, ces deux horloges étant synchronisée précisément entre elles). Dans ce premier mode de réalisation, l'unité de traitement 18 détermine des coordonnées xB, yB, zB, qui repèrent la balise 3 dans le référentiel x,y,z lié au 30 système sonar 10, comme suit. Dans un premier temps, La profondeur zB, et la célérité moyenne correspondante est déterminée en résolvant l'équation suivante : 37 2 ZB = Cmoy (ZB). tB. - (a. cmoy(zB)) Puis les coordonnées « horizontales » x6, y6, sont déterminées conformément aux formules suivantes : x.8 = (x8). t8. a. cos(9) et YB = cLy(za).tB. a. sin(p). 5 L'unité de traitement 18 peut en outre être programmée pour déterminer une position de la balise 3 dans un référentiel dans lequel la balise est fixe, tel que le référentiel terrestre XT, YT, A. Pour cela, l'unité de traitement 18 acquiert une position du navire 1, repérée dans le référentiel terrestre, au moyen du système de localisation 2 équipant ce navire, le navire 1 étant situé à la même position que 10 lors de la réception de d'onde acoustique émise par la balise 3. La position de la balise 3 est déterminée ensuite en fonction de la position du navire 1 dans le référentiel terrestre, et de la position (xB,yB,zB) de la balise par rapport au système sonar 10. Dans un deuxième mode de réalisation de ce troisième procédé, l'unité 15 de traitement 18 détermine la position latérale (xB,yB) du centre B de l'émetteur 4, par rapport au système sonar 10, sur la base de la constante de propagation a de l'onde acoustique reçue, et d'une profondeur z1B de la balise connue au préalable. Cette profondeur z'B est connue de l'unité de traitement 18 par exemple parce que des données représentatives de cette profondeur ont été transmises de 20 la balise 3 au système sonar 10 par l'intermédiaire de l'onde acoustique émise par la balise. Ce deuxième mode de réalisation permet avantageusement de déterminer la position de la balise 3, même si le système sonar 10 et cette balise 3 ne partagent pas de référence temporelle commune. 25 Les coordonnées x6, yB, qui repèrent la position latérale du centre B de l'émetteur 4, dans le plan horizontal (x,y), sont déterminées chacune en fonction du produit de : - la différence entre la profondeur de l'émetteur 4, z'6, et la profondeur de la tête sonar, ici considérée comme nulle, et de 30 - la tangente de l'angle de propagation moyen en,',,, qui a été défini plus haut (lors de la description du premier procédé). La tangente tan(emoy) de cet angle de propagation moyen ermy est déterminée en fonction de la constante de propagation a de l'onde reçue et de la 38 célérité moyenne cmoy(z'B) à la profondeur z1B de l'émetteur 4 : tan(0,,,y) = 8). a/ \11 - 8). a) 2 Dans ce deuxième mode de réalisation, l'unité de traitement 18 détermine alors les coordonnées xB, yB conformément aux formules suivantes, par exemple : 5 x 8 = 8. cmoy(z'B). a I (cnioy (z' 8). a)2. cos((p) Ye = 8. cmoy(z'B). a / - (c moy (z' 8). a)2. sin(q)). En variante, la tangente de l'angle de propagation moyen emy pourrait toutefois être déterminée en fonction de la célérité arithmétique moyenne E(z' 8), au lieu d'être déterminée en fonction de la célérité (harmonique) moyenne 10 cmoy(z'B). Dans ce cas, la tangente de l'angle de propagation moyen Ormy est déterminée conformément à la formule suivante : tan(0m,,y) = C.(z' 8). a / 11 - (C. (z' 8). a)2 Par ailleurs, comme dans le premier mode de réalisation, l'unité de traitement 18 peut en outre être programmée pour déterminer une position de la 15 balise 3 dans un référentiel dans lequel la balise est fixe, tel que le référentiel terrestre XT, YT, ZT. Il est noté que ce troisième procédé est appliqué ici à la localisation d'une balise immergée, mais qu'il peut aussi s'appliquer à la détermination de la position d'un autre système émergé émettant des ondes acoustiques, tel qu'un 20 sous-marin distinct du navire 1 précité, par exemple, et cela quel soit le mode de réalisation considéré. Différentes variantes peuvent être apportées aux procédés et au système sonar qui ont été décrits ci-dessus. Tout d'abord, le système sonar pourrait équiper un navire submersible, 25 au lieu d'un navire de surface, que ce navire soit autonome ou qu'il accueille un équipage. Dans ce cas, ce n'est pas directement la profondeur d'un point du fond marin, ou la profondeur de l'émetteur de la balise à localiser qui est déterminée, mais une différence Az=zp-z, entre la profondeur Zp de ce point du fond (ou de cet 30 émetteur) et la profondeur z' de l'antenne réceptrice du système sonar 10. 39 Les premier, deuxième et troisième procédés sont alors exécutés par les unités de pilotage 17 et de traitement 18 de manière complètement analogue à ce qui a été décrit plus haut, mais en déterminant la célérité moyenne cmoy(zp) conformément à la formule suivante : zP 1 cmoy(zp) (zp - zo) c(z)dz (F5') 5 et, de manière optionnelle, en déterminant la célérité arithmétique moyenne conformément à la formule suivante : 1 zp (zp z0) - il C(zp) = c(z) dz (F4') D'autre part, le modèle monocouche employé pour déterminer la position du point sondé ou de la balise émettrice pourrait être paramétré de manière légèrement différente de ce qui a été présenté ci-dessus. Ainsi, à titre d'exemple, 10 l'angle de propagation moyen emy pourrait présenter un léger écart relatif par rapport au cas (optimal) présenté ci-dessus, pour lequel sin(Bmoy) = ce. cntoy(zp), (ou, en variante, pour lequel sin(Omoy) = a. E(zp)). Il est toutefois préférable que cet écart relatif reste inférieur à un pour mille (soit 0,001), car cela permet de conserver, pour les profondeurs et écarts latéraux déterminés par l'unité de 15 traitement, des précisions élevées, comparables à ce qui a été présenté plus haut. Par ailleurs, comme déjà indiqué, le système sonar pourrait présenter une structure plus simple : il pourrait par exemple comporter, au lieu de deux antennes linéaires multi-capteurs complètes dont les récepteurs (par exemple au nombre de quatre au moins) sont espacés de X/2 (où X est la longueur d'onde 20 moyenne des onde acoustiques émises), disposées "en croix de Mills", un système de réception lacunaire (dont les récepteurs sont espacés de plus que 212) composé à minima de 3 récepteurs non alignés. Dans ce cas, l'ambiguïté angulaire peut être résolue (c'est-à-dire levée), à titre d'exemple, en émettant des signaux larges bandes (dont le contenu spectral s'étend sur une bande de 25 fréquences large). En variante encore, les unités de commande, de pilotage et de traitement pourraient être réalisées sous la forme d'une même unité électronique, ou, au contraire, REVENDICATIONS1. Procédé de détermination d'une différence de profondeur (zp, Az ; zm,i, zm,2, zm ; ze), ou d'un écart latéral (yp ; xm ; xB, yB) par rapport à la verticale, entre deux points d'un milieu subaquatique (E), le procédé comprenant les étapes suivantes : - émission d'une onde acoustique dans le milieu subaquatique (E) au moyen d'au moins un émetteur (15 ; 15, 16 ; 4), - réception de ladite onde acoustique au moyen d'une antenne réceptrice (14 ; 14, 13) comprenant plusieurs récepteurs (16 ; 16, 15), ces récepteurs délivrant une pluralité respective de signaux de réception (s) lors de la réception de ladite onde acoustique, - détermination d'une constante de propagation (a ; am,1, am,2) de l'onde acoustique reçue, en fonction desdits signaux de réception (s), ladite constante de propagation étant égale au sinus d'un angle de réception (8°) repérant la direction de réception de l'onde acoustique par rapport à la verticale, divisé par une vitesse de propagation locale (co) des ondes acoustiques à la profondeur de ladite antenne réceptrice, - détermination d'un temps de propagation (t ; tM,1, tm,2 ; tB) de l'onde acoustique, en fonction d'une durée séparant les instants d'émission et de réception de l'onde acoustique, et - détermination de la différence de profondeurs (zp, Az ; zm,i, ZM,2. zm ; z8) ou de l'écart latéral (yp ; xm ; xB, yB) par rapport à la verticale, entre l'antenne réceptrice (13, 14) et l'émetteur (4), ou entre l'antenne réceptrice (14 ; 14, 13) et un élément immergé (P ; M) réfléchissant ladite onde acoustique lors de sa propagation depuis l'émetteur (15 ; 15, 16) jusqu'à l'antenne réceptrice (14 ; 14, 13), en fonction : - d'un produit dudit temps de propagation (t ; tM,1, tm,2 ; tB) par une célérité moyenne (cmoy) de l'onde acoustique, et en fonction - d'un angle de propagation moyen (Omoy), défini entre la verticale et une direction de propagation effective de l'onde acoustique, l'angle de propagation moyen (0m0y) étant déterminé en fonction de ladite constante de propagation (a ; am,1, am2). 2. Procédé selon la revendication 1, dans lequel il est prévu de 41 déterminer à la fois ladite différence de profondeur (zp, Az ; ZM,1, zm,2, zm ; zB) et ledit écart latéral (yp ; xm ; xB, yB) par rapport à la verticale. 3. Procédé selon l'une des revendications 1 ou 2, dans lequel ladite différence de profondeur (Zp, Az ; zm,i, ZM,2, zm ; zB) est déterminée, de manière à 5 présenter un écart relatif inférieur à un pour mille par rapport : - au produit dudit temps de propagation (t ; tut tm.2 ; tB) par ladite célérité moyenne (cmoy), - multiplié par le cosinus dudit angle de propagation moyen (6,0y). 4. Procédé selon l'une des revendications 1 à 3, dans lequel ledit écart 10 latéral est déterminé, de manière à présenter un écart relatif inférieur à un pour mille par rapport : - au produit dudit temps de propagation (t ; tut tm.2 , tB) par ladite célérité moyenne (cmoy), - multiplié par le sinus dudit angle de propagation moyen (ernoy). 15 5. Procédé selon l'une des revendications 1 à 4, dans lequel ladite célérité moyenne (cmoy) est représentative d'une moyenne harmonique d'une pluralité de vitesses de propagation locales (c), présentées par les ondes acoustiques à une pluralité respective de profondeurs (z) allant d'une profondeur de l'émetteur (4) jusqu'à une profondeur de l'antenne réceptrice (13, 14), ou d'une 20 profondeur dudit élément immergé (P ; M) jusqu'à la profondeur de l'antenne réceptrice (14 ; 14, 13). 6. Procédé selon l'une des revendications 1 à 5, dans lequel l'angle de propagation moyen (ernoy) est déterminé de sorte que son sinus présente un écart relatif inférieur à un pour mille par rapport au produit de ladite constante de 25 propagation (a ; am,1, am,2) par ladite célérité moyenne (cmoy). 7. Procédé selon l'une des revendications 1 à 5, dans lequel l'angle de propagation moyen (cmoy) est déterminé en fonction en outre d'une célérité arithmétique moyenne (C) qui est représentative d'une moyenne arithmétique d'une pluralité de vitesses de propagation locales (c), présentées par les ondes 30 acoustiques à une pluralité respective de profondeurs (z) allant de la profondeur de l'émetteur (4) jusqu'à la profondeur de l'antenne réceptrice (13, 14), ou de la profondeur dudit élément immergé (P ; M) jusqu'à la profondeur de l'antenne réceptrice (14 ; 14, 13). 8. Procédé selon la revendication 7, dans lequel l'angle de propagation 42 moyen (8,-noy) est déterminé de sorte que son sinus présente un écart relatif inférieur à un point mille par rapport à ladite constante de propagation (a) multipliée par ladite célérité arithmétique moyenne (E). 9. Procédé selon l'une des revendications 1 à 8, dans lequel l'émetteur 5 (15 ; 15, 16) et l'antenne réceptrice (14 ; 14, 13) équipent un même système sonar (10), et au cours duquel ce système sonar (10) détermine la différence de profondeurs (zp, Az ; zm.t, zm,2, ZM) ou l'écart latéral (yp ; xm) par rapport à la verticale entre l'antenne réceptrice et un élément immergé (P ; M) réfléchissant ladite onde acoustique lors de sa propagation depuis l'émetteur (15 ; 15, 16) 10 jusqu'à l'antenne réceptrice (14 ; 14, 13). 10. Procédé selon l'une des revendications 1 à 8, dans lequel l'émetteur (4) et l'antenne réceptrice (13, 14) équipent respectivement deux systèmes distincts (3, 10) situés à des profondeurs respectives différentes (zB, z0), et au cours duquel le système (10) muni de l'antenne réceptrice (13, 14) détermine la 15 différence de profondeurs (zB) ou l'écart latéral (xB, yB) par rapport à la verticale entre l'antenne réceptrice (13, 14) et ledit émetteur (4). 11. Procédé selon l'une des revendications 1 à 10, comprenant en outre les étapes suivantes : - détermination, au moyen d'une sonde positionnée successivement à 20 plusieurs profondeurs (z) différentes, d'une pluralité de vitesses de propagation locales (c) présentées respectivement, à chacune desdites profondeurs (z), par les ondes acoustiques, puis - détermination de ladite célérité moyenne (c,',y) en fonction de ladite pluralité de vitesses de propagation locales (c) déterminées précédemment. 25 12. Procédé selon la revendication 9, comprenant en outre les étapes suivantes : - après l'émission et la réception de ladite onde acoustique, déplacement du système sonar (10), puis - émission d'une autre onde acoustique dans le milieu subaquatique (E) 30 au moyen de dudit émetteur, ou d'un autre émetteur (16) du système sonar, de sorte que ladite autre onde atteigne ledit élément immergé (M), - après réflexion sur ledit élément immergé (M), réception de ladite autre onde acoustique au moyen de ladite antenne réceptrice, ou d'une autre antenne émettrice (13) du système sonar (10), 43 - détermination d'une autre constante de propagation (am,2), pour ladite autre onde acoustique reçue, en fonction des signaux de réception (s) délivrés par les récepteurs (16) de ladite antenne réceptrice, ou de ladite autre antenne réceptrice (16) suite à la réception de cette autre onde acoustique, ladite autre 5 constante de propagation (am,2) étant égale au sinus d'un autre angle de réception repérant la direction de réception de ladite autre onde acoustique par rapport à la verticale, divisé par une vitesse de propagation locale (c0) des ondes acoustiques à la profondeur de ladite antenne réceptrice ou de ladite autre antenne réceptrice, - détermination d'un autre temps de propagation (42), en fonction d'une 10 autre durée séparant les instants d'émission et de réception de ladite autre onde acoustique, et - détermination de ladite célérité moyenne (cmoy) en fonction : - du temps de propagation (tM,1) et de la constante de propagation (a ) de ladite onde acoustique, et 15 - du temps de propagation (42) et de la constante de propagation (a ) de ladite autre onde acoustique. 13. Procédé selon la revendication 12, dans lequel ladite célérité moyenne (cmoy) est déterminée de manière à ce qu'un écart entre - ladite différence de profondeur (zio,i), et 20 - une autre différence de profondeur (zm,2) entre l'antenne réceptrice (14) et l'élément immergé (M) soit inférieur à un seuil donné, ladite autre différence de profondeur (zm2) étant déterminée en fonction : - d'un produit dudit autre temps de propagation (42) par ladite 25 célérité moyenne (cmoy), et en fonction - d'un autre angle de propagation moyen, défini entre la verticale et une direction de propagation effective de ladite autre onde acoustique, ledit autre angle de propagation moyen étant déterminé en fonction de ladite autre constante de propagation (aM,2). 30 14. Procédé de détermination d'une vitesse de propagation locale (c) des ondes acoustiques dans un milieu subaquatique (E), comprenant les étapes suivantes : - pour une pluralité d'éléments immergés (M, M', M") situés à des profondeurs différentes, détermination d'une pluralité respective de valeurs de 44 ladite célérité moyenne (cm,y), chacune desdites valeurs étant déterminée conformément au procédé défini par la revendication 12 ou 13, puis - détermination d'au moins une vitesse de propagation locale (c) des ondes acoustiques, à une profondeur (z) donnée fixe, en fonction de ladite 5 pluralité de valeurs de la célérité moyenne (cm'y) déterminées précédemment. - détermination en chacun desdits éléments immergés, à partir de ladite vitesse de propagation locale déterminée précédemment et de ladite constante de propagation, d'un angle repérant une direction de propagation de ladite onde acoustique par rapport à la verticale. 10 15. Procédé de détermination d'une position (xB, yB) d'un émetteur (4) immergé dans un milieu subaquatique (E), comprenant les étapes suivantes : - émission d'une onde acoustique dans le milieu subaquatique par ledit émetteur (4), dont la profondeur (zB) est connue, - réception de ladite onde acoustique au moyen d'une antenne réceptrice 15 (13, 14) équipant un navire (1), de surface ou submersible, l'antenne réceptrice (13, 14) comprenant plusieurs récepteurs (15, 16) délivrant une pluralité respective de signaux de réception (s), - détermination d'une constante de propagation (a) de l'onde acoustique reçue, en fonction desdits signaux de réception (s), ladite constante de 20 propagation étant égale au sinus d'un angle de réception (On) repérant la direction de réception de l'onde acoustique par rapport à la verticale dans un plan vertical (PI) contenant l'émetteur (4) et l'antenne réceptrice (13, 14), divisé par une vitesse de propagation locale (c0) des ondes acoustiques à la profondeur (z0) de ladite antenne réceptrice, 25 - détermination d'une direction de propagation effective de l'onde acoustique en fonction du produit de ladite constante de propagation (a) et d'une célérité moyenne (cmoy ; e) de propagation entre la profondeur (zB) de l'émetteur et la profondeur (z0) de l'antenne réceptrice, et - détermination d'un écart latéral (xB, yB) par rapport à la verticale entre 30 l'antenne réceptrice (13, 14) et l'émetteur (4), en fonction du produit de : - la différence entre la profondeur (zB) de l'émetteur et la profondeur (z0) de l'antenne réceptrice, et de - la tangente d'un angle de propagation moyen (6'y) défini entre ladite direction de propagation effective et la verticale. 45 16. Système sonar (10) comprenant : - une antenne réceptrice (13, 14) configurée pour réceptionner une onde acoustique qui a été émise précédemment dans un milieu subaquatique par au moins un émetteur (15, 16, 4), l'antenne réceptrice (13, 14) comprenant plusieurs 5 récepteurs (15, 16) aptes à délivrer une pluralité respective de signaux de réception (s) lors de la réception de ladite onde acoustique, et - une unité de traitement (18) électronique programmée pour : - déterminer une constante de propagation (a ; am.i, am,2) de l'onde acoustique reçue, en fonction desdits signaux de réception (s), ladite constante de 10 propagation étant égale au sinus d'un angle de réception (0,) repérant la direction de réception de l'onde acoustique par rapport à la verticale, divisé par une vitesse de propagation locale (c0) des ondes acoustiques à la profondeur de ladite antenne réceptrice, - déterminer un temps de propagation (t ; tut tm,2 ; tB) de l'onde 15 acoustique, en fonction d'une durée séparant les instants d'émission et de réception de l'onde acoustique, et pour - déterminer une différence de profondeur (4,, 4z ; ZM,1, zm.2. zm ; zB) ou un écart latéral (yr, ; xM ; xB, yB) par rapport à la verticale, entre l'antenne réceptrice (13, 14) et l'émetteur (4) ou entre l'antenne réceptrice (14 ; 14, 13) et un élément 20 immergé (P ; M) réfléchissant ladite onde acoustique lors de sa propagation depuis l'émetteur (15 ; 15, 16) jusqu'à l'antenne réceptrice (14 ; 14, 13), en fonction : - d'un produit dudit temps de propagation (t ; tm,i, tm,2 ; tg) par une célérité moyenne de l'onde acoustique, et en fonction 25 - d'un angle de propagation moyen (6,0y), défini entre la verticale et une direction de propagation effective de l'onde acoustique, l'angle de propagation moyen (6,0y) étant déterminé en fonction de ladite constante de propagation (a ; am,2). 17. Système sonar (10) comprenant : 30 - une antenne réceptrice (13, 14) configurée pour réceptionner une onde acoustique qui a été émise précédemment dans un milieu subaquatique par un émetteur (4) dont la profondeur (zB) est connue, l'antenne réceptrice (13, 14) comprenant plusieurs récepteurs (15, 16) aptes à délivrer une pluralité respective de signaux de réception (s) lors de la réception de ladite onde acoustique, et 46 - une unité de traitement (18) électronique programmée pour : - déterminer une constante de propagation (a) de l'onde acoustique reçue, en fonction desdits signaux de réception (s), ladite constante de propagation étant égale au sinus d'un angle de réception (8e) repérant la direction 5 de réception de l'onde acoustique par rapport à la verticale dans un plan vertical (PI) contenant l'émetteur (4) et l'antenne réceptrice (13,14), divisé par une vitesse de propagation locale (c0) des ondes acoustiques à la profondeur (z0) de ladite antenne réceptrice, - déterminer une direction de propagation effective de l'onde 10 acoustique reçue, en fonction du produit de ladite constante de propagation (a) et d'une la célérité moyenne (c''y ; E) de propagation entre la profondeur (zB) de l'émetteur et la profondeur (z0) de l'antenne réceptrice, - acquérir ou lire dans une mémoire une donnée représentative de la profondeur (zB) de l'émetteur (4), et pour 15 - déterminer un écart latéral (xB, yB) par rapport à la verticale, entre l'antenne réceptrice (13, 14) et l'émetteur (4), en fonction du produit de : - la différence entre la profondeur (zB) de l'émetteur et la profondeur (z0) de l'antenne réceptrice, et de - la tangente d'un angle de propagation moyen (9,,y) défini 20 entre ladite direction de propagation effective et la verticale.