<div>Pavement drainage systems are one of the key drivers of pavement function and longevity, and effective drain maintenance can significantly extend a pavement's service life. Maintenance of these drains, however, is often hampered by the challenge of locating the drains. Ground Penetrating Radar (GPR) typically offers a rapid and effective method to detect these underground targets. However, typical detection schema that rely upon the observation of the hyperbolic return from a GPR scan of a buried conduit still tend to miss many of the older drains beneath pavements as they may be partially or fully filled with sediment and/or may be fabricated from clay or other earthen materials, yielding a return signal that is convolved with significant background noise. </div><div><br></div><div>To manage this challenge, this work puts forward an improved background noise and clutter reduction method to enhance the target signals in what amounts to a constructed environment that tends to have more consistent subsurface properties than one might encounter in a general setting. Within this technique, two major algorithms are employed. Algorithm 1 is the core of this method, and plays the role of reducing background noise and clutter. Algorithm 2 is supplementary, and helps eliminate anomalous discontinuous returns generated by the equipment itself, which could otherwise lead to false detection indications in the output of Algorithm 1. Instead of traditional 2-D GPR images, the result of the proposed algorithms is a 1-D plot along the survey line, highlighting a set of “points of interest” that could indicate buried drain locations identified at any given GPR operating frequency. Subsurface exploration using two different operating frequencies, 900 MHz and 400 MHz herein, is then employed to further enhance detection confidence. Points of interest are ultimately coded to define the confidence of the detection. Comparing the final result of proposed algorithms with the original GPR images, the improved algorithm is demonstrated to provide significantly improved detection results, and could potentially be applied to similar problems in other contexts.</div><div><br></div><div>Besides the background reduction methods, a group of simulations performed using GPRMAX2D software are examined to explore the influence of road cross-section designs on sub-pavement drainage conduit GPR signatures, and evaluate the effectiveness of alternate GPR antennae configurations in locating these buried conduits in different ground conditions. Two different models were explored to simulate conduit detection. In addition, different pipe and soil conditions were modeled, such as pipe size, pipe material, soil moisture level, and soil type. Four different quantitative measurements are used to analyze GPR performance based on different key factors. The four measurements are 1) signal to background ratio (SBR) in dB; 2) signal to receiver noise ratio (SNR) in dB; 3) signal energy in Volts; and 4) average signal band power in Watts.</div><div><br></div><div>The water and clay content of subsurface soil can significantly influence the detection results obtained from ground penetrating radar (GPR). Due to the variation of the material properties underground, the center frequency of transmitted GPR signals shifts to a lower range as wave attenuation increases. Examination of wave propagation in the subsurface employing an attenuation filter based on a linear system model shows that received GPR signals will be shifted to lower frequencies than those originally transmitted. The amount of the shift is controlled by a wave attenuation factor, which is determined by the dielectric constant, electric conductivity, and magnetic susceptibility of the transmitted medium. This work introduces a receiver-transmitter-receiver dual-frequency configuration for GPR that employs two operational frequencies for a given test - one higher and one slightly lower - to take advantage of this phenomenon to improve subpavement drain detection results. In this configuration, the original signal is transmitted from the higher frequency transmitter. After traveling through underground materials, the signal is received by two receivers with different frequencies. One of the receivers has the same higher center frequency as the transmitter, and the other receiver has a lower center frequency. This configuration can be expressed as Rx(low-frequency)-Tx(high-frequency)-Rx(high-frequency) and was applied in both laboratory experiments and field tests. Results are analyzed in the frequency domain to evaluate and compare the properties of the signal obtained by both receivers. The laboratory experiment used the configuration of Rx(400MHz)-Tx(900MHz)-Rx(900MHz). The field tests, in addition to the configuration used in the lab tests, employed another configuration of Rx(270MHz)-Tx(400MHz)-Rx(400MHz) to obtain more information about this phenomenon. Both lab and field test results illustrate the frequency-shift phenomenon described by theoretical calculations. Based on the power spectrum for each signal, the lower frequency antenna typically received more energy (higher density values) at its peak frequency than the higher frequency antenna.</div>

No description available.

This thesis concerns the study, development and analysis of innovative artificial intelligence (AI)-driven optimization techniques within the System-by-Design (SbD) framework aimed at efficiently addressing the computational complexity inherent in advanced electromagnetic (EM) problems. By leveraging the available a-priori information as well as the proper integration of machine learning (ML) techniques with intelligent exploration strategies, the SbD paradigm enables the effective and reliable solution of the EM problem at hand, with user-defined performance and in a reasonable amount of time. The flexibility of the AI-driven SbD framework is demonstrated in practice with the implementation of two solution strategies to address the fully non-linear inverse scattering problem (ISP) for the detection and imaging of buried objects in ground penetrating radar (GPR)-based applications, and to address the design and optimization of mm-wave automotive radars that comply multiple challenging and contrasting requirements. A comprehensive set of numerical experiments is reported to demonstrate the efficacy and computational efficiency of the SbD-based optimization techniques in solving complex EM problems.

Les premiers mètres à centaines de mètres de la proche surface terrestre sont le siège de processus naturels dont la compréhension requiert une caractérisation fine de la subsurface, via une estimation quantifiée de ses paramètres. Le géoradar est un outil de prospection indirecte à même d'ausculter les milieux naturels et d'en estimer les propriétés électriques (permittivité et conductivité). Basé sur la propagation d'ondes électromagnétiques à des fréquences allant du MHz à quelques GHz, le géoradar est utilisé à des échelles et pour des applications variées concernant la géologie, l'hydrologie ou le génie civil. Dans ce travail de thèse, je propose une méthode d'imagerie quantitative des propriétés électriques sur des sections 2D de la subsurface, à partir de données radar acquises à la surface du sol. La technique mise en oeuvre est l'inversion des formes d'ondes, qui utilise l'intégralité du champ d'ondes enregistré.Dans une première partie, je présente les principes physiques et l'outil de modélisation numérique utilisés pour simuler la propagation des ondes électromagnétiques dans les milieux hétérogènes à deux dimensions. Pour cela, un algorithme de différences finies en domaine fréquentiel développé dans le cadre des ondes visco-acoustiques est adapté au problème électromagnétique 2D grâce à une analogie mathématique.Dans une deuxième partie, le problème d'imagerie est formulé sous la forme d'une optimisation multi-paramètre puis résolu avec l'algorithme de quasi-Newton L-BFGS. Cet algorithme permet d'estimer l'effet de la matrice Hessienne, dont le rôle est crucial pour la reconstruction de paramètres de différents types comme la permittivité et la conductivité. Des tests numériques montrent toutefois que l'algorithme reste sensible aux échelles utilisées pour définir ces paramètres. Dans un exemple synthétique représentatif de la proche surface, il est cependant possible d'obtenir des cartes 2D de permittivité et de conductivité à partir de données de surface, en faisant intervenir des facteurs d'échelle et de régularisation visant à contraindre les paramètres auxquelles l'inversion est la moins sensible. Ces facteurs peuvent être déterminés en analysant la qualité de l'ajustement aux données, sans hypothèse a priori autre que la contrainte de lissage introduite par la régularisation.Dans une dernière partie, la méthode d'imagerie est confrontée à deux jeux de données réelles. Dans un premier temps, l'examen de données expérimentales permet de tester la précision des simulations numériques vis-à-vis de mesures effectuées en environnement contrôlé. La connaissance des cibles à imager permet en outre de valider la méthodologie proposée pour l'imagerie multiparamètre dans des conditions très favorables puisqu'il est possible de calibrer le signal source et de considérer l'espace libre environnant les cibles comme modèle initial pour l'inversion.Dans un deuxième temps, j'envisage le traitement d'un jeu de données radar multi-offsets acquises au sein d'un massif calcaire. L'interprétation de ces données est rendue beaucoup plus difficile par la complexité du milieu géologique environnant, ainsi que par la méconnaissance des caractéristiques précises des antennes utilisées. L'application de la méthode d'inversion des formes d'ondes à ces données requiert donc une étape préliminaire impliquant une analyse de vitesse plus classique, basée sur les arrivées directes et réfléchies, et des simulations numériques dans des modèles hypothétiques à même d'expliquer une partie des données. L'estimation du signal source est effectuée à partir d'arrivées sélectionnées, simultanément avec des valeurs moyennes de conductivité et de hauteur d'antennes de façon à reproduire au mieux les amplitudes observées. Un premier essai d'inversion montre que l'algorithme est capable d'expliquer les données dans la gamme de fréquences considérée et de reconstruire une ébauche des principaux réflecteurs. / The quantitative characterization of the shallow subsurface of the Earth is a critical issue for many environmental and societal challenges. Ground penetrating radar (GPR) is a geophysical method based on the propagation of electromagnetic waves for the prospection of the near subsurface. With central frequencies between 10~MHz and a few GHz, GPR covers a wide range of applications in geology, hydrology and civil engineering. GPR data are sensitive to variations in the electrical properties of the medium which can be related, for instance, to its water content and bring valuable information on hydrological processes. In this work, I develop a quantitative imaging method for the reconstruction of 2D distributions of permittivity and conductivity from GPR data acquired from the ground surface. The method makes use of the full waveform inversion technique (FWI), originating from seismic exploration, which exploits the entire recorded radargrams and has been proved successful in crosshole GPR applications.In a first time, I present the numerical forward modelling used to simulate the propagation of electromagnetic waves in 2D heterogeneous media and generate the synthetic GPR data that are compared to the recorded radargrams in the inversion process. A frequency-domain finite-difference algorithm originally developed in the visco-acoustic approximation is adapted to the electromagnetic problem in 2D via an acoustic-electromagnetic mathematical analogy.In a second time, the inversion scheme is formulated as a fully multiparameter optimization problem which is solved with the quasi-Newton L-BFGS algorithm. In this formulation, the effect of an approximate inverse Hessian is expected to mitigate the trade-off between the impact of permittivity and conductivity on the data. However, numerical tests on a synthetic benchmark of the literature display a large sensitivity of the method with respect to parameter scaling, showing the limits of the L-BFGS approximation. On a realistic subsurface benchmark with surface-to-surface configuration, it has been shown possible to ally parameter scaling and regularization to reconstruct 2D images of permittivity and conductivity without a priori assumptions.Finally, the imaging method is confronted to two real data sets. The consideration of laboratory-controlled data validates the proposed workflow for multiparameter imaging, as well as the accuracy of the numerical forward solutions. The application to on-ground GPR data acquired in a limestone massif is more challenging and necessitates a thorough investigation involving classical processing techniques and forward simulations. Starting permittivity models are derived from the velocity analysis of the direct arrivals and of the reflected events. The estimation of the source signature is performed together with an evaluation of an average conductivity value and of the unknown antenna height. In spite of this procedure, synthetic data do not reproduce the observed amplitudes, suggesting an effect of the radiation pattern of the shielded antennae. In preliminary tests, the inversion succeeds in fitting the data in the considered frequency range and can reconstruct reflectors from a smooth starting model.

Cette thèse a pour objectif d’établir une vision globale et continue mer-terre-mer d’une barrière littorale paraglaciaire, en combinant une approche morphologique et architecturale. L’originalité de l’isthme de Miquelon-Langlade (NO Atlantique) entre les deux îles éponymes réside en sa diversité et complexité de systèmes (flèches, cordons littoraux, lagune) qui s’emboitent formant une barrière composite évoluant le long de deux façades ouvertes sur l’océan. Ce travail a été rendu possible grâce à l’association de données topographiques (GPS-RTK), de géophysique marine (sismique HR) et terrestre (géo-radar), sédimentaires (carottes, échantillons de surface) ainsi que des datations par OSL. L’étude de cette barrière a permis 1) d’identifier l’ensemble des dépôts d’une séquence paraglaciaire (dépôts de chute, bas niveau, augmentation et haut niveau marin); 2) de pointer l’importance de la géométrie des îles, du substratum rocheux ainsi que de la topographie héritée dans la formation de la barrière; 3) de définir l’architecture interne de beach ridges mixte sable-galet; et 4) de pointer le potentiel de ces systèmes en tant marqueur du paléo-niveau marin tout en proposant la première courbe de tendance de l’évolution du RSL pour l’Archipel de Saint-Pierre-et-Miquelon sur les derniers 3000 ans (+1mm/an). L’association de l’ensemble de ces informations a permis de proposer un schéma de mise en place de la barrière, depuis le retrait de la calotte (13 700 ans) jusqu’à sa forme actuelle, évoluant dans un contexte d’augmentation du RSL et de fort apport sédimentaire, tout en proposant une vision de détail de la formation de la plaine de beach ridges au cours des derniers 3000 ans. / The purpose of this PhD thesis was to define a complete and continuous sea-land-sea view of a paraglacial coastal barrier system, combining morphological and architectural approaches. Particularity of the Miquelon-Langlade isthmus (NO Atlantic) lies in its diversity and complexity of systems (spits, beach ridges, and lagoon) that fit together to form a composite barrier prograding along two open coasts. Investigations of this barrier included topographic (RTK-GPS), ground-penetrating radar (GPR) seismic, sedimentary (cores, augers, samples) and chronologic (OSL) data collections. This study of this paraglacial composite barrier make possible 1) to characterize all deposits of the paraglacial sequence (FSST, LST, TST, HST); 2) to emphasize the primary influence of islands geometry, bedrock and inherited topography on the coastal progradational pattern; 3) to provide a better understanding of the internal architecture of depositional processes and driving mechanisms of mixed sand-and-gravel beach ridges; 4) to examine the potential of these mixed beach-ridge systems to record and preserve paleo-sea-level information, to produce the first sea level trends for the Saint-Pierre-et-Miquelon archipelago over the last 3000 years. The combination of all these data provide a model of development of the barrier since the ice sheet retreat (13 700 yrs) to its present form, in a context of RSL rise and high local glaciogenic sediment supply, and a detailed view of the formation of the beach ridge plain over the last 3000 yrs.

La compréhension des écoulements préférentiels dans la zone non saturée en fond de bassins d’infiltration est une condition préalable à une gestion optimale de ces ouvrages. Le basin d’infiltration étudié repose sur un dépôt fluvioglaciaire hétérogène constitué d’un lithofaciès majoritaire contenant des lentilles de sable et de gravier. Une étude géophysique (radar) a permis de caractériser cette hétérogénéité sédimentologique le long de radar-grammes (sections de quelques mètres de profondeurs et de plusieurs dizaines de mètre de long). Les écoulements ont été modélisés pour toutes ces sections ainsi qu’une section fictive remplie uniquement du lithofaciès majoritaire (témoin) et pour divers évènements pluvieux type incluant une chronique météorologique complète. Les résultats montrent clairement que l'hétérogénéité lithologique a un effet significatif sur l’écoulement avec la mise en place d’écoulements hétérogènes du fait de phénomènes de barrière capillaire au niveau des lentilles de sable et de gravier en non-saturé. Cet effet se répercute sur les chemins de l’eau dans le dépôt et sur les flux en surface (infiltration et ruissellement). Les résultats montrent aussi que cet effet dépend fortement de la section considérée. La variabilité spatiale de l’hétérogénéité lithologique induit donc une variabilité spatiale des écoulements préférentiels. Certains paramètres architecturaux et géométriques relatifs à l’hétérogénéité lithologique ont été proposés pour expliquer l’effet de l’hétérogénéité lithologique sur les écoulements préférentiels, de sorte à proposer une première tentative de modélisation des écoulements à méso-échelle. / Understanding preferential flows in the vadose zone is a prerequisite for a proper management of infiltration basins. The studied infiltration basin lies on a strongly heterogeneous glaciofluvial deposit made of a predominant lithofacies that embeds lenses of sand and gravel. Geophysics (GPR) were used to characterize the lithological heterogeneity along radar-grammes (sections few meters deep and dozens of meters long). Flow was modelled for all these sections and compared to a synthetic section filled only with the predominant lithofacies (control). Numerical results show that, under unsaturated conditions, the lithological heterogeneity has a strong impact on flow with the establishment of preferential flows induced by capillary-barrier effects in the vicinity of lenses of sand and gravel. These effects have also consequences on water fluxes at surface (infiltration and runoff) and depend strongly on the section considered for numerical modelling. It can thus be concluded that the spatial variability of lithological heterogeneity induces a spatial variability of preferential flow. Geometrical and architectural parameters were defined to quantity lithological heterogeneity and link it to preferential flows, with the aim to propose a first step towards modelling flow in the basin at the meso-scale. / O entendimento do fluxo na zona não saturada em formações heterogêneas é um prérequisito para uma gestão adequada de poluentes em bacias de infiltração inseridas nessas formações. Este estudo aborda o efeito da heterogeneidade litológica nos processo de escoamento numa bacia de infiltração com geologia inserida num depósito fluvioglacial. A bacia de infiltração é a Bacia Django Reinhardt, existente em Lyon. Esta bacia foi objeto de vários estudos que demonstraram algumas avaliações do impacto da heterogeneidade do solo. Os objetivos desses trabalhos foram todos baseados apenas no estudo sedimentológico de uma trincheira. Nenhum estudo abordou a potencial variabilidade espacial dos atributos do solo. Neste trabalho, foi modelizado o fluxo na zona vadosa através de três casos: drenagem interna, infiltração de água durante um evento de precipitação pluviométrica constante e solicitação hidrológica a partir de um fluxo variável escoamento superficial. Estes cálculos foram realizados para várias seções, que foram previamente caracterizadas quanto a sua heterogeneidade estrutural usando a técnica não destrutiva do GPR com uma frequência de 200MHz e 100MHz. O comportamento hidráulico dessas seções foi comparado com a de uma seção homogênea. Foi comparada o efeito da resolução da imagem nas estimativas das variáveis de fluxo. Alguns indicadores arquiteturais, geométricos e texturais foram propostos para esclarecer e avaliar a variabilidade espacial do efeito da heterogeneidade do solo sobre o fluxo em condições não saturadas. Os resultados mostram claramente que a heterogeneidade gera impactos nas variáveis de fluxo em condições não saturadas. Os impactos dependem da seção considerada.Os resultados mostram também que: a) a heterogeneidade estrutural do subsolo na meso-escala apresenta um efeito sobre a geração de escoamento preferencial e isso independe do tipo de solicitação aplicada na seção; b) O efeito da heterogeneidade se mostrou independentende das condições de contorno aplicadas na superfície; c) Nos períodos de ocorrência de chuvas intensas quem controla o comportamento das variáveis de fluxo é a intensidade de precipitação. No entanto, nos períodos em que ocorre ausência de precipitação, as variáveis de fluxo são controladas pela heterogeneidade da zona não saturada; d) A utilização de diferentes frequências influencia na estimativa das variáveis de fluxo. Isso é decorrente do impacto na quantificação da heterogeneidade decorrente do efeito da resolução de cada antena; e) Os indicadores não apresentaram relação direta com as variáveis de fluxo estudadas. Isso demonstra que o efeito da heterogeneidade não pode ser explicado por um único fator, mas por um conjunto de fatores.

Adjacent precast, prestressed box beam bridges have a history of poor performance and have been observed to exhibit common types of deterioration including longitudinal cracking, concrete spalling, and deterioration of the concrete top flange. The nature of these types of deterioration leads to uncertainty of the extent and effect of deterioration on structural behavior. Due to limitations in previous research and understanding of the strength of deteriorated box beam bridges, conservative assumptions are being made for the assessment and load rating of these bridges. Furthermore, the design of new box beam bridges, which can offer an efficient and economical solution, is often discouraged due to poor past performance. Therefore, the objective of this research is to develop improved recommendations for the inspection, load rating, and design of adjacent box beam bridges. Through a series of bridge inspections, deteriorated box beams were identified and acquired for experimental testing. The extent of corrosion was determined through visual inspection, non-destructive evaluation, and destructive evaluation. Non-destructive tests (NDT) included the use of connectionless electrical pulse response analysis (CEPRA), ground penetrating radar (GPR), and half-cell potentials. The deteriorated capacity was determined through structural testing, and an analysis procedure was developed to estimate deteriorated behavior. A rehabilitation procedure was also developed to restore load transfer of adjacent beams in cases where shear key failures are suspected. Based on the understanding of deterioration developed through study of deteriorated adjacent box beam bridges, improved inspection and load rating procedure are provided along with design recommendations for the next generation of box beam bridges.

[ES] Los esfuerzos que hacen los países, en conjunto con organizaciones mundiales, tales como IWA (por International Water Association), ONU-Agua y OMS (Organización Mundial de la Salud), para mitigar el impacto ambiental en el campo de la hidráulica urbana son considerados de vital importancia. Sin embargo, la escasez de los recursos hídricos en el mundo aumenta diariamente. Esto viene dado por el aumento constante de la demanda en los sectores industrial, agrícola y urbano, provocado por el aumento poblacional y el cambio climático. Los administradores de los sistemas de abastecimiento de agua (WSSs, por sus siglas en inglés, water supply systems) se han visto desafiados a suplir la creciente demanda de los diferentes sectores con la cantidad, calidad y eficiencia necesarios y, a su vez, reducir el desperdicio y el mal uso del recurso. Desde esta perspectiva, las fugas de agua son el mayor problema que enfrentan los administradores de estas empresas de servicios públicos. Las fugas en una red provocan problemas de salud, de escasez, económicos y medioambientales. El uso de técnicas de inspección no destructivas debe favorecer una rápida identificación de problemas, para realizar acciones posteriores de reparación en la red. Este trabajo hace uso del GPR (siglas en inglés de ground penetrating radar) como técnica de inspección no destructiva porque: favorece la exploración del subsuelo sin causar alteraciones al medio, es de fácil aplicación y, además, permite obtener pseudo imágenes del subsuelo. Uno de los objetivos de este documento es identificar y extraer características de una fuga en un WSS mediante imágenes de GPR, con el fin último de recrear las fugas a través de modelos 3D. Se realizaron ensayos de laboratorio bajo condiciones controladas donde se emuló una parcela en la cual se había enterrado una tubería con una pequeño orifico que simula una fuga de agua; tras introducir agua al sistema, se realizaron prospecciones con el GPR. Una vez finalizada la exploración del subsuelo, dado que las imágenes de GPR en bruto obtenidas no son fácilmente interpretables por personal no experto, tales imágenes fueron sometidas a procesamiento de datos que favorezcan su fácil interpretación. Este documento presenta dos metodologías de procesamiento de datos que permiten obtener imágenes a partir de las cuales es posible identificar tanto los componentes del sistema como la fuga y su alcance. Las metodologías de tratamiento de datos aplicadas en este documento son una metodología basada en sistemas multi-agente y el filtro de varianza, metodología basada en parámetros estadísticos de segundo orden. Posteriormente, tras aplicar estas metodologías de procesamiento a las imágenes, se sometieron los resultados a un análisis que facilitase la mejor elección evitando la subjetividad del experto. Bajo este concepto, este documento propone el uso conjunto de técnicas multicriterio. Se utilizó el Proceso de Jerarquía Analítica Difusa (FAHP, por sus siglas en inglés, Fuzzy Analytical Hierarchy Process), que permite ponderar varios criterios de evaluación, con el propósito de mitigar la incertidumbre que caracterizan los juicios de los expertos, en conjunto con el método ELECTRE III para obtener la clasificación final de alternativas, todo esto de la manera más objetiva posible. Los resultados de este documento son satisfactorios, permitiendo obtener amplio conocimiento de las fugas y su interacción con el subsuelo, proporcionando pautas para desarrollar posteriormente metodologías de automatización que permitan localizar, seguir y predecir problemas en los WSSs. / [CA] Els esforços que fan els països en conjunt amb organitzacions mundials, como ara IWA (per International Water Association), ONU-Agua i OMS (per Organització Mundial de la Salut), per a mitigar l'impacte ambiental en el camp de la hidràulica urbana són considerats de vital importància. No obstant això, l'escassetat dels recursos hídrics en el món augmenta diàriament, donat per l'augment constant de la demanda en els sectors industrial, agrícola i urbà, provocat per l'augment poblacional i el canvi climàtic. Els administradors dels sistemes d'abastiment d'aigua (WSSs, per les seus sigles en anglès, water supply systems) s'han vist desafiats a suplir la creixent demanda dels diferents sectors amb la quantitat, qualitat i eficiència necessaris i, al seu torn, reduir el desaprofitament i el mal ús del recurs. Enfocant aquesta perspectiva, les pèrdues d'aigua són el problema més gran fet front pels directors d'aquestes utilitats. Les pèrdues d'aigua en una xarxa provoquen problemes de salut, d'escassetat, econòmics i mediambientals. L'ús de tècniques d'inspecció no destructives que afavoreixen una ràpida identificació per a realitzar accions de reparació posteriors en la xarxa. Aquest treball fa ús del GPR (sigles en anglès per ground penetrating radar) com a tècnica d'inspecció no destructiva perquè afavoreix l'exploració del subsol sense causar alteracions al entorn, és de fàcil aplicació i a més permet obtenir pseudo imatges del subsol. Un dels objectius d'aquest document és identificar i extraure característiques d'una pèrdua en un WSS mitjançant imatges de GPR, amb la fi última de recrear les pèrdues a través de models 3D. Es van realitzar assajos de laboratori sota condicions controlades on es va emular una parcel¿la en la qual s'ha enterrat una canonada amb una xicotet forat que simula una pèrdua d'aigua; després d'introduir aigua al sistema, s'obtenen prospeccions amb el GPR. Una vegada finalitzada l'exploració del subsol, atès que les imatges de GPR en brut obtingudes no són fàcilment interpretables per personal no expert, són sotmeses a processament de dades que afavorisquen la seua fàcil interpretació. Aquest document presenta dues metodologies de processament de dades que permeten obtenir imatges de les quals és possible identificar tant els components del sistema com la pèrdua i el seu abast. Les metodologies de tractament de dades aplicades en aquest document són una metodologia basada en multi-agents (MABS, per les seves sigles en anglès, multi-agent-based systems) i el filtre de variància, metodologia basada en paràmetres estadístics de segon ordre. Posteriorment, després d'aplicar aquestes metodologies de processament a les imatges se sotmeten els resultats a una anàlisi que faciliti la millor elecció evitant la subjectivitat de l'expert. Sota aquest concepte, aquest document proposa l'ús conjunt de tècniques de decisió multi-criteri (MCDM, per les seves sigles en anglès, multi-criteria decision-making). Es va utilitzar el Procés de Jerarquia Analítica Difusa (FAHP, per les seves sigles en anglès, Fuzzy Analytical Hierarchy Process) el qual s'utilitza per a ponderar diversos criteris d'avaluació, amb el propòsit de mitigar la incertesa que caracteritzen els judicis dels experts, en conjunt amb el mètode ELECTRE III, per a obtenir la classificació final d'alternatives, tot això de la manera més objectiva possible. Els resultats d'aquest document són satisfactoris, permetent obtenir ampli coneixement de les pèrdues d'aigua i la seua interacció amb el subsol, donant-nos la pauta per a desenvolupar posteriorment metodologies d'automatització que permeten localitzar, seguir i predir problemes en els WSSs. / [EN] The efforts made by the countries in collaboration with world organizations, such as IWA (for International Water Association), UN-Water and WHO (for World Health Organization), to mitigate the environmental impact in the field of urban hydraulics are considered of vital importance. However, the scarcity of water resources in the world increases daily, given by the constant increase in demand in the industrial, agricultural and urban sectors, caused by the population increase and the climate change. Managers of water supply systems (WSSs) are challenged to supply the growing demand of different sectors with sufficient quantity, quality and efficiency and, in turn, reduce waste and misuse of the resource. Focusing this perspective, water leaks are the biggest problem faced by the managers of these utilities. Leaks in a network cause health, shortage, economic and environmental problems. The use of non-destructive inspection techniques favors rapid identification to carry out subsequent repair actions on the network. This work makes use of the GPR (ground penetrating radar) as a non-destructive inspection technique because: it favors the exploration of the ground without causing alterations to the environment, it is easy to apply, and also allows to obtain pseudo images of the subsoil. This document presents two data processing methodologies that allow obtaining images from which it is possible to identify both the system components and the leak and its scope. The data treatment methodologies applied in this document are a multi-agent-based system (MABS) methodology and the variance filter, a methodology based on second-order statistical parameters. Subsequently, after applying these processing methodologies to the images, the results are subjected to an analysis that eases the best choice, avoiding expert's subjectivity. Under this concept, this document proposes the joint use of two multi-criteria decision-making (MCDM) methods. The Fuzzy Analytical Hierarchy Process (FAHP) is used first to weight various evaluation criteria, in order to mitigate the uncertainty that characterize the experts' judgments, in conjunction with the ELECTRE III method, to obtain the final classification of alternatives in the most objective way. The results of this document are satisfactory, allowing to obtain extensive knowledge of leaks and their interaction with the subsoil, giving a guideline to subsequently develop automation methodologies that allow locating, monitoring and predicting problems in WSSs. / Part of this work has been developed under the support of Fundación Carolina PhD and short-term scholarship program / Ocaña Levario, SJ. (2021). Análisis, caracterización y modelación 3D de fugas de agua en sistemas de abastecimiento de agua mediante imágenes de GPR [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/163677 / TESIS

The emergence of Ultra-wideband (UWB) technology application has yielded tremendous and vital impacts in the field of microwave wireless communications. These applications include military radar imaging, security screening, and tumour detection, especially for early detection of breast cancer. These indicators have stimulated and inspired many researchers to make the best use of this promising technology. UWB technology challenges such as antenna design, the problem of imaging reconstruction techniques, challenges of severe signal attenuation and dispersion in high loss material. Others are lengthy computational time demand and large computer memory requirements are prevalent constraints that need to be tackled especially in a large scale and complex computational electromagnetic analysis. In this regard, it is necessary to find out recently developed optimisation techniques that can provide solutions to these problems. In this thesis, designing, optimisation, development, measurement, and analysis of UWB antennas for near-field microwave imaging applications are considered. This technology emulates the same concept of surface penetrating radar operating in various forms of the UWB spectrum. The initial design of UWB monopole antennas, including T-slots, rectangular slots, and hexagonal slots on a circular radiating patch, was explicitly implemented for medical imaging applications to cover the UWB frequency ranging from 3.1 GHz to 10.6 GHz. Based on this concept, a new bow-tie and Vivaldi UWB antennas were designed for a through-the-wall imaging application. The new antennas were designed to cover a spectrum on a lower frequency ranging from 1 GHz - 4 GHz to ease the high wall losses that will be encountered when using a higher frequency range and to guarantee deeper penetration of the electromagnetic wave. Finally, both simulated and calculated results of the designed, optimised antennas indicate excellent agreement with improved performance in terms of return loss, gain, radiation pattern, and fidelity over the entire UWB frequency. These breakthroughs provided reduced computational time and computer memory requirement for useful, efficient, reliable, and compact sensors for imaging applications, including security and breast cancer detection, thereby saving more lives. / Tertiary Education Trust Fund (TET Fund) Supported by the Nigerian Defence Academy (NDA)

Humanitarian clearance of minefields is most often carried out by hand, conventionally using a a metal detector and a probe. Detection is a very slow process, as every piece of detected metal must treated as if it were a landmine and carefully probed and excavated, while many of them are not. The process can be safely sped up by use of Ground-Penetrating Radar (GPR) to image the subsurface, to verify metal detection results and safely ignore any objects which could not possibly be a landmine. In this thesis, we explore the possibility of using Full Wave Inversion (FWI) to improve GPR imaging for landmine detection. Posing the imaging task as FWI means solving the large-scale, non-linear and ill-posed optimisation problem of determining the physical parameters of the subsurface (such as electrical permittivity) which would best reproduce the data. This thesis begins by giving an overview of all the mathematical and implementational aspects of FWI, so as to provide an informative text for both mathematicians (perhaps already familiar with other inverse problems) wanting to contribute to the mine detection problem, as well as a wider engineering audience (perhaps already working on GPR or mine detection) interested in the mathematical study of inverse problems and FWI.We present the first numerical 3D FWI results for GPR, and consider only surface measurements from small-scale arrays as these are suitable for our application. The FWI problem requires an accurate forward model to simulate GPR data, for which we use a hybrid finite-element boundary-integral solver utilising first order curl-conforming N\'d\'{e}lec (edge) elements. We present a novel `line search' type algorithm which prioritises inversion of some target parameters in a region of interest (ROI), with the update outside of the area defined implicitly as a function of the target parameters. This is particularly applicable to the mine detection problem, in which we wish to know more about some detected metallic objects, but are not interested in the surrounding medium. We may need to resolve the surrounding area though, in order to account for the target being obscured and multiple scattering in a highly cluttered subsurface. We focus particularly on spatial sensitivity of the inverse problem, using both a singular value decomposition to analyse the Jacobian matrix, as well as an asymptotic expansion involving polarization tensors describing the perturbation of electric field due to small objects. The latter allows us to extend the current theory of sensitivity in for acoustic FWI, based on the Born approximation, to better understand how polarization plays a role in the 3D electromagnetic inverse problem. Based on this asymptotic approximation, we derive a novel approximation to the diagonals of the Hessian matrix which can be used to pre-condition the GPR FWI problem.