La connaissance de l'environnement marin est nécessaire pour un grand nombre d'applications dans le domaine de l'acoustique sous-marine comme la communication, la localisation et détection sonar et la surveillance des mammifères marins. Il constitue le moyen principal pour éviter les interférences néfastes entre le milieu naturel et les actions industriels et militaires conduites en zones côtières.Notre travail de thèse se place dans un contexte de sonar actif avec des fréquences allant de 1 kHz à 10 kHz pour des distances de propagations allant de 1 km à plusieurs dizaines de kilomètres. Nous nous intéressons particulièrement aux environnements de propagation grands fonds, à l'utilisation des antennes industrielles comme les antennes de flancs, les antennes cylindriques et les antennes linéaires remorquées, et à l'utilisation de signaux large bande afin de travailler avec des résolutions en distance et en vitesse très élevées. Le travail de recherche présenté dans ce mémoire est dédié à la recherche de nouveaux paramètres discriminants pour la classification de cible sous-marine en sonar actif et notamment à l'estimation de l'immersion instantanée.Cette étude présente : (1) les calculs de nouvelles bornes de Cramer-Rao pour la position d'une cible en distance en et en profondeur, (2) l'estimation conjointe de la distance et de l'immersion d'une cible à partir de la mesure des temps d'arrivées et des angles d'élévations sur une antenne surfacique et (3) l'estimation conjointe de la distance, de l'immersion et du gisement d'une cible à partir de la mesure des temps d'arrivées et des pseudo-gisements sur une antenne linéaire remorquée.Les méthodes développées lors de cette étude ont été validées sur des simulations, des données expérimentales à petite échelle et des données réelles en mer. / The knowledge of the marine environment is required for many underwater applications such as communications, sonar localization and detection, and marine mammals monitoring. It enables preventing harmful interference between the natural environment and industrial and military actions in coastal areas.This thesis work concentrates upton the context of active sonar with frequencies from 1 kHz to 10 kHz and long propagation ranges from 1 km to several tens of kilometers. We also concentrates upon deep water environment, the use of industrial arrays such as cylindrical arrays, flank arrays and linear towed arrays, and the use of large time-bandwidth signals in order to obtain high distance and speed resolutions. This research work is dedicated to the research of new features for the underwater target classification in active sonar, and specifically to the instantaneous target-depth estimation.This thesis presents: (1) calculations of new Cramer-Rao bounds for the target-position in range and in depth, (2) the joint estimation of the target-depth and the target-range from the arrival time and elevation angle measures with a surface array, (3) the joint estimation of the target-depth, the target-range and the target-bearing from the arrival time and pseudo-bearing angle measures with a linear towed array.The methods presented in this manuscript have been benchmarked on simulation, on reduced-scale experimental data and real marine data.

Waterborne acoustic signals carry information about the ocean environment. Ocean geoacoustic inversion is the task of estimating environmental parameters from received acoustic signals by matching the measured sound with the predictions of a physics-based model. A lower bound on the uncertainty associated with environmental parameter estimates, the Cramér-Rao bound, can be calculated from the Fisher information, which is dependent on derivatives of a physics-based model. Physics-based preconditioners circumvent the need for variable step sizes when computing numerical derivatives. This work explores the feasibility of using a neural network to perform geoacoustic inversion for environmental parameters and their associated uncertainties from ship noise spectrogram data. To train neural networks, a synthetic dataset is generated and tested for generalizability against 31 measurements taken during the SBCEX2017 study of the New England Mud Patch.

In this thesis, acoustic inversion methods are employed to estimate array element locations and the geoacoustic properties of the seabed using measured acoustic data consisting of noise from a surface ship in the Gulf of Mexico. The array element localization utilizes relative travel-time information obtained by cross-correlating the recorded time series of ship noise received at spatially separated hydrophones. The relative travel-time data are used in an inversion, based on the regularized least-squares method and the acoustic ray tracing equations, to obtain improved estimates of the receiver and source positions and their uncertainties. Optimization and Bayesian matched-field inversion methods are employed to estimate seabed geoacoustic properties and their uncertainties in the vicinity of a bottom-moored vertical line array using the recorded surface ship noise. This study is used to test the feasibility of matched-field methods to detect temporal changes in the geoacoustic properties of the seabed near a known gas hydrate mound in the Gulf of Mexico. Finally, a synthetic study is performed that demonstrates how ignoring environmental range dependence of seabed sound speed and water depth in matched-field inversion can lead to biases in the estimated geoacoustic parameters. The study considers the distributions of optimal parameter estimates obtained from a large number of range-independent inversions of synthetic data generated for random range-dependent environments. Range-independent Bayesian inversions are also performed on selected data sets and the marginal parameter distributions are examined. Both hard- and soft-bottom environments are examined at a number of scales of variability in sound speed and water depth.

Ce travail de thèse propose d’analyser les signaux issus d’une source omnidirectionnelle sphérique réfléchis par un milieu sédimentaire stratifié et enregistré par une antenne d’hydrophones, en vue de caractériser quantitativement les sédiments marins aux moyennes fréquences, i.e. comprises entre 1 et 10 kHz. La recherche développée dans ce manuscrit propose une méthodologie facilitant la recherche des paramètres géoacoustiques du milieu avec la méthode des sources images, ainsi qu’un ensemble de solutions techniques appropriées afin d’améliorer cette méthode d’inversion récemment développée. La méthode des sources images repose sur une modélisation physique de la réflexion des ondes émises par une source sur un milieu stratifié sous l’approximation de Born. Par conséquent, la réflexion de l’onde sur le milieu stratifié peut être représentée par une collection de sources images, symétriques de la source réelle par rapport aux interfaces, dont les positions spatiales sont liées à la vitesse des ondes acoustiques et aux épaisseurs des couches. L’étude se décline en deux volets : traitements des signaux et inversion des paramètres géoacoustiques. La première partie du travail est focalisée sur le développement de la méthode des sources images. La méthode originelle se basait sur la construction de cartes de migration et de semblance de signaux pour déterminer les paramètres d’entrée de l’algorithme d’inversion qui sont des temps de trajet et des angles d’arrivée. Afin d’éviter cette étape, nous détectons les temps d’arrivée avec l’opérateur d’énergie de Teager-Kaiser (TKEO) et nous trouvons les angles par une méthode de triangulation. Le modèle d’inversion a été ensuite intégré en prenant en compte la possibilité de déformation de l’antenne. Cette partie se termine par une nouvelle approche qui combine TKEO et des méthodes temps fréquence afin d’avoir une bonne détection du temps d’arrivée dans le cas de signaux fortement bruités. Sur le plan du modèle et de l’inversion géoacoustique, nous proposons tout d’abord une description précise du modèle direct en introduisant le concept de sources images virtuelles. Cette étape permet de mieux comprendre l’approche développée. Ensuite, nous proposons une extension de la méthode des sources image pour l’inversion de paramètres géoacoustiques supplémentaires : la densité, l’atténuation et la vitesse des ondes de cisaillement. Cette extension est basée sur les résultats de l’inversion originelle (estimation du nombre de strates, de leur épaisseur, et de la vitesse des ondes de compression) ainsi que sur l’utilisation de l’amplitude des signaux réfléchis. Ces améliorations et extensions de la méthode des sources images sont illustrées par leur application sur des signaux synthétiques et des signaux réels issus d’expérimentations en cuve et à la mer. Les résultats obtenus sont très satisfaisants, tant au niveau des performances de calcul que de la qualité des estimations fournies. / This thesis aims at analyzing the signals emitted from a spherical omnidirectional source reflected by a stratified sedimentary environment and recorded by a hydrophone array in order to characterize quantitatively the marine sediments at medium frequencies, i.e. between 1 and 10 kHz. The research developed in this manuscript provides a methodology to facilitate the estimation of medium geoacoustic parameters with the image source method, and some appropriate technical solutions to improve this recently developed inversion method. The image source method is based on a physical modeling of the wave reflection emitted from a source by a stratified medium under the Born approximation. As result, the reflection of the wave on the layered medium can be represented by a set of image sources, symmetrical to the real source with respect to the interfaces, whose spatial positions are related to the sound speeds and the thicknesses of the layers. The study consists of two parts : signal processing and inversion of geoacoustic parameters. The first part of the work is focused on the development of the image source method. The original method was based on migration and semblance maps of the recorded signals to determine the input parameters of the inversion algorithm which are travel times and arrival angles. To avoid this step, we propose to determine the travel times with the Teager-Kaiser energy operator (TKEO) and the arrival angles are estimate with a triangulation approach. The inversion model is then integrated, taking into account the possible deformation of the antenna. This part concludes with a new approach that combines TKEO and time-frequency representations in order to have a good estimation of the travel times in the case of noisy signals. For the modeling and geoacoustic inversion part, we propose first an accurate description of the forward model by introducing the concept of virtual image sources. This idea provides a deeper understanding of the developed approach. Then, we propose an extension of the image sources method to the estimation of supplementary geoacoustic parameters : the density, the absorption coefficient, and the shear wave sound speed. This extension is based on the results of the original inversion (estimation of the number of layers, their thicknesses, and the pressure sound speeds) and on the use of the amplitudes of the reflected signals. These improvements and extents of the image source method are illustrated by their applications on both synthetic and real signals, the latter coming from tank and at-sea measurements. The obtained results are very satisfactory, from a computational point of view as well as for the quality of the provided estimations.

La mer Ionienne est une zone à la géodynamique active en raison de la convergence entre les plaques Nubie et Eurasie. Elle correspond aux derniers stades de vie d’un océan, la Téthys. De ce fait, la tectonique et la sédimentation y sont très réduites ; et les faibles flux sédimentaires permettent ainsi d’enregistrer une multiplicité de processus sédimentaires.Néanmoins, c’est une zone encore peu étudiée d’un point de vue sédimentaire, notamment sur les échelles de temps de l’ordre de la centaine de milliers d’années. Ce travail est basé sur une étude sédimentologique d’archives marines complétée par des données acoustiques (bathymétrie et multi-faisceau) issue de campagnes océanographiques du SHOM. L’analyse détaillée des faciès et des séquences sédimentaire a permis, dans un premier temps, d’établir un calendrier des risques naturels (séismes, tsunamis, volcanisme), leurs sources et leurs processus de dépôt dans le bassin pour les derniers 330 000 ans. Puis, dans un second temps, de retracer l’évolution sur les derniers 400 000 ans de la circulation et de l’oxygénation des eaux de fond dans le bassin ionien et l’influence du détroit de Sicile, et notamment de la plate-forme de Malte, sur les échanges entre les bassins occidental et oriental. L’intégration des données sédimentologiques dans un modèle géoacoustique développé par le SHOM a finalement permis de déterminer l’impact des variations sédimentaires (distribution spatiale, lithologie, stratification) sur la propagation des ondes acoustique pour différentes gammes de fréquences (300 Hz - 3000 Hz) et d’angle d’incidence (0 -90°) et d’établir une cartographie de la réponse sédimentaire du le signal acoustique. / The Ionian Sea is an active geodynamic area because of the convergence between theNubia and the Eurasia plate. It corresponds to the last stage of the Tethys ocean life. Therefore,the tectonics and the sedimentation are much reduced; and the low sedimentary supply enables torecord a multiplicity of sedimentary processes. Nevertheless, this area still poorly studied from asedimentary point of view, in particular on timescales on the order of hundred thousand years.This work is based on a sedimentological study of marine archives supplemented by acoustic data(bathymetry and multibeam imagery) recovered during oceanographic campaign leaded by theSHOM. The detailed sedimentary analysis of facies and sequences allows, at first, to established acalendar of the natural hazard (earthquakes, tsunamis, volcanism), their origins and theirdepositional processes into the basin over the last 330 000 years. Then, secondly, to reconstructthe evolution over the last 400 000 years of the circulation and the oxygenation of bottom waterthrough the Ionian basin and the influence of the strait of Sicily, in particular of the Malta Plateau,on the exchanges between the western and the eastern basins. The integration of thesedimentological data in a geoacoustic modelling developed by the SHOM finally allowed todetermine the impact of the sedimentary variability (special distribution, lithology, stratification)on the acoustic waves propagation for various frequency bands (300 Hz - 3000 Hz) and incidentangle (0 - 90°) and to established a mapping of the sedimentary answer of the acoustic signal.

Originally developed in the 1970s for the optimal control of systems governed by partial differential equations, the adjoint method has found several successful applications, e.g. in meteorology with large-scale 3D or 4D atmospheric data assimilation schemes, for carbon cycle data assimilation in biogeochemistry and climate research, or in oceanographic modelling with efficient adjoint codes of ocean general circulation models.<p><p>Despite the variety of applications in these research fields, adjoint methods have only very recently drawn attention from the ocean acoustics community. In ocean acoustic tomography and geoacoustic inversion, where the inverse problem is to recover unknown acoustic properties of the water column and the seabed from acoustic transmission data, the solution approaches are typically based on travel time inversion or standard matched-field processing in combination with metaheuristics for global optimization. <p><p>In order to complement the adjoint schemes already in use in meteorology and oceanography with an ocean acoustic component, this thesis is concerned with the development of the adjoint of a full-field acoustic propagation model for shallow water environments. <p><p>In view of the increasing importance of global ocean observing systems such as the European Seas Observatory Network, the Arctic Ocean Observing System and Maritime Rapid Environmental Assessment (MREA) systems for defence and security applications, the adjoint of an ocean acoustic propagation model can become an integral part of a coupled oceanographic and acoustic data assimilation scheme in the future. <p><p>Given the acoustic pressure field measured on a vertical hydrophone array and a modelled replica field that is calculated for a specific parametrization of the environment, the developed adjoint model backpropagates the mismatch (residual) between the measured and predicted field from the receiver array towards the source.<p><p>The backpropagated error field is then converted into an estimate of the exact gradient of the objective function with respect to any of the relevant physical parameters of the environment including the sound speed structure in the water column and densities, compressional/shear sound speeds, and attenuations of the sediment layers and the sub-bottom halfspace. The resulting environmental gradients can be used in combination with gradient descent methods such as conjugate gradient, or Newton-type optimization methods tolocate the error surface minimum via a series of iterations. This is particularly attractive for monitoring slowly varying environments, where the gradient information can be used to track the environmental parameters continuously over time and space.<p><p>In shallow water environments, where an accurate treatment of the acoustic interaction with the bottom is of outmost importance for a correct prediction of the sound field, and field data are often recorded on non-fully populated arrays, there is an inherent need for observation over a broad range of frequencies. For this purpose, the adjoint-based approach is generalized for a joint optimization across multiple frequencies and special attention is devoted to regularization methods that incorporate additional information about the desired solution in order to stabilize the optimization process.<p><p>Starting with an analytical formulation of the multiple-frequency adjoint approach for parabolic-type approximations, the adjoint method is progressively tailored in the course of the thesis towards a realistic wide-angle parabolic equation propagation model and the treatment of fully nonlocal impedance boundary conditions. A semi-automatic adjoint generation via modular graph approach enables the direct inversion of both the geoacoustic parameters embedded in the discrete nonlocal boundary condition and the acoustic properties of the water column. Several case studies based on environmental data obtained in Mediterranean shallow waters are used in the thesis to assess the capabilities of adjoint-based acoustic inversion for different experimental configurations, particularly taking into account sparse array geometries and partial depth coverage of the water column. The numerical implementation of the approach is found to be robust, provided that the initial guesses are not too far from the desired solution, and accurate, and converges in a small number of iterations. During the multi-frequency optimization process, the evolution of the control parameters displays a parameter hierarchy which clearly relates to the relative sensitivity of the acoustic pressure field to the physical parameters. <p><p>The actual validation of the adjoint-generated environmental gradients for acoustic monitoring of a shallow water environment is based on acoustic and oceanographic data from the Yellow Shark '94 and the MREA '07 sea trials, conducted in the Tyrrhenian Sea, south of the island of Elba.<p> <p>Starting from an initial guess of the environmental control parameters, either obtained through acoustic inversion with global search or supported by archival in-situ data, the adjoint method provides an efficient means to adjust local changes with a couple of iterations and monitor the environmental properties over a series of inversions. <p><p>In this thesis the adjoint-based approach is used, e.g. to fine-tune up to eight bottom geoacoustic parameters of a shallow-water environment and to track the time-varying sound speed profile in the water column. <p><p>In the same way the approach can be extended to track the spatial water column and bottom structure using a mobile network of sparse arrays.<p><p>Work is currently being focused on the inclusion of the adjoint approach into hybrid optimization schemes or ensemble predictions, as an essential building block in a combined ocean acoustic data assimilation framework and the subsequent validation of the acoustic monitoring capabilities with long-term experimental data in shallow water environments. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Detailed information about the oceanic environment is essential for many applications in the field of marine geology, marine biology, coastal engineering, and marine operations. Especially, knowledge of the properties of the sediment body is often required. Acoustic remote sensing techniques have become highly attractive for classifying the sea bottom and for mapping the sediment properties, due to their high coverage capabilities and low costs compared to common sampling methods. In the last decades, a number of different acoustic devices and related techniques for analyzing their signals have evolved. Each sensor has its specific application due to limitations in the frequency range and resolution. In practice, often a single acoustic tool is chosen based on the current application, supported by other non-acoustic data where required. However, different acoustic remote sensing techniques can supplement each other, as shown in this thesis. Even more, a combination of complementary approaches can contribute to the proper understanding of sound propagation, which is essential when using sound for environmental classification purposes. This includes the knowledge of the relation between acoustics and sediment properties, the focus of this thesis. Providing a detailed three dimensional picture of the sea bottom sediments that allows for gaining maximum insight into this relation is aimed at.<p><p><p>Chapters 4 and 5 are adapted from published work, with permission:<p>DOI:10.1121/1.3569718 (link: http://asadl.org/jasa/resource/1/jasman/v129/i5/p2878_s1) and<p>DOI:10.1109/JOE.2010.2066711 (link: http://ieeexplore.ieee.org/xpl/articleDetails.jsp?tp=&arnumber=5618582&queryText%3Dsiemes)<p>In reference to IEEE copyrighted material which is used with permission in this thesis, the IEEE does not endorse any of the Université libre de Bruxelles' products or services.<p> / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished