Este trabalho trata da simulação da propagação de ondas eletromagnéticas (1-D) em meios geológicos de propriedades físicas conhecidas (espessura, condutividade e permissividade elétrica), do registro das ondas EM refletidas (dados sintéticos) para uma antena receptora, e do uso desse registro em algoritmos de inversão e otimização que procuram estimar as propriedades físicas dos meios geológicos. Adotou-se modelos geológicos estratificados de pequena espessura (< 1,2 m) e altas frequências para a onda eletromagnética (800, 1000 e 1200 MHz). A propagação da onda eletromagnética é modelada por meio da resolução das equações de Maxwell no método FDTD. Os métodos de inversão Quasi-Newton e Otimização Ant Colony modificado são aplicados sobre os dados sintéticos para estimar os parâmetros elétricos para cada camada geológica. Ambos os métodos foram aplicados alternadamente para aumentar a precisão e a convergência ao longo da profundidade. Os métodos de inversão foram capazes de estimar simultaneamente duas propriedades eletromagnéticas do modelo geológico: a permissividade elétrica e a condutividade elétrica. Os métodos de inversão alcançaram bons resultados quando executados simultaneamente sobre os dados sintéticos em 3 diferentes frequências. Exemplos de estimativas dos perfis de condutividade e permissividade elétricas unidimensionais são apresentados, com e sem inserção de ruídos nos dados. Os resultados indicam que a combinação dos métodos de inversão (ACO modificado e Quasi-Newton) pode fornecer bons resultados para as estimativas dos parâmetros físicos de meios geológicos e geotécnicos em meios rasos. Além disso, essa combinação de métodos de inversão abre novas perspectivas para o processamento de dados georradar multi-canais. As investigações também mostram que a simulação dos perfis sintéticos por meio do método FDTD pode propagar erros em espessuras maiores de subsolo, devido ao fato de ser um método explícito. / This work deals with simulation of a 1-D electromagnetic wave propagating into a geological structure of known physical properties (thickness, electric conductivity and permittivity), the record of the reflected EM wave (synthetic data), and the use of these records in inversion and optimization algorithms to estimate back the physical properties of the geological structure. A stratified and thin (< 1.2 m) geological structure was initially constructed in order to be scanned by high frequency EM waves (800, 1000, and 1200 MHz). The EM wave propagation is simulated by Maxwell equations through FDTD method. The Quasi-Newton inversion and Ant Colony Optimization methods were applied into synthetic data to estimate original physical parameters of each geological layer. Both methods were applied in order to increase precision and convergence along depth. These methods were able to simultaneously estimate two physical properties of the geological structure: electrical permittivity and conductivity. The methods showed good results when applied simultaneously upon synthetic data of all three frequencies. Electric permittivity and conductivity profiles are shown with and without noise in the data. The results indicate that combined inversion methods can show good results to estimate physical properties of thin geological and geothecnical structures. The combined inversion methods open new perspectives to multichannel GPR data processing. The investigations also show that simulation of synthetic data through FDTD can propagate errors when dealing with thick geological structures, due to the fact that FDTD is an explicit method.

Este trabalho trata da simulação da propagação de ondas eletromagnéticas (1-D) em meios geológicos de propriedades físicas conhecidas (espessura, condutividade e permissividade elétrica), do registro das ondas EM refletidas (dados sintéticos) para uma antena receptora, e do uso desse registro em algoritmos de inversão e otimização que procuram estimar as propriedades físicas dos meios geológicos. Adotou-se modelos geológicos estratificados de pequena espessura (< 1,2 m) e altas frequências para a onda eletromagnética (800, 1000 e 1200 MHz). A propagação da onda eletromagnética é modelada por meio da resolução das equações de Maxwell no método FDTD. Os métodos de inversão Quasi-Newton e Otimização Ant Colony modificado são aplicados sobre os dados sintéticos para estimar os parâmetros elétricos para cada camada geológica. Ambos os métodos foram aplicados alternadamente para aumentar a precisão e a convergência ao longo da profundidade. Os métodos de inversão foram capazes de estimar simultaneamente duas propriedades eletromagnéticas do modelo geológico: a permissividade elétrica e a condutividade elétrica. Os métodos de inversão alcançaram bons resultados quando executados simultaneamente sobre os dados sintéticos em 3 diferentes frequências. Exemplos de estimativas dos perfis de condutividade e permissividade elétricas unidimensionais são apresentados, com e sem inserção de ruídos nos dados. Os resultados indicam que a combinação dos métodos de inversão (ACO modificado e Quasi-Newton) pode fornecer bons resultados para as estimativas dos parâmetros físicos de meios geológicos e geotécnicos em meios rasos. Além disso, essa combinação de métodos de inversão abre novas perspectivas para o processamento de dados georradar multi-canais. As investigações também mostram que a simulação dos perfis sintéticos por meio do método FDTD pode propagar erros em espessuras maiores de subsolo, devido ao fato de ser um método explícito. / This work deals with simulation of a 1-D electromagnetic wave propagating into a geological structure of known physical properties (thickness, electric conductivity and permittivity), the record of the reflected EM wave (synthetic data), and the use of these records in inversion and optimization algorithms to estimate back the physical properties of the geological structure. A stratified and thin (< 1.2 m) geological structure was initially constructed in order to be scanned by high frequency EM waves (800, 1000, and 1200 MHz). The EM wave propagation is simulated by Maxwell equations through FDTD method. The Quasi-Newton inversion and Ant Colony Optimization methods were applied into synthetic data to estimate original physical parameters of each geological layer. Both methods were applied in order to increase precision and convergence along depth. These methods were able to simultaneously estimate two physical properties of the geological structure: electrical permittivity and conductivity. The methods showed good results when applied simultaneously upon synthetic data of all three frequencies. Electric permittivity and conductivity profiles are shown with and without noise in the data. The results indicate that combined inversion methods can show good results to estimate physical properties of thin geological and geothecnical structures. The combined inversion methods open new perspectives to multichannel GPR data processing. The investigations also show that simulation of synthetic data through FDTD can propagate errors when dealing with thick geological structures, due to the fact that FDTD is an explicit method.

Este trabalho trata da simulação da propagação de ondas eletromagnéticas (1-D) em meios geológicos de propriedades físicas conhecidas (espessura, condutividade e permissividade elétrica), do registro das ondas EM refletidas (dados sintéticos) para uma antena receptora, e do uso desse registro em algoritmos de inversão e otimização que procuram estimar as propriedades físicas dos meios geológicos. Adotou-se modelos geológicos estratificados de pequena espessura (< 1,2 m) e altas frequências para a onda eletromagnética (800, 1000 e 1200 MHz). A propagação da onda eletromagnética é modelada por meio da resolução das equações de Maxwell no método FDTD. Os métodos de inversão Quasi-Newton e Otimização Ant Colony modificado são aplicados sobre os dados sintéticos para estimar os parâmetros elétricos para cada camada geológica. Ambos os métodos foram aplicados alternadamente para aumentar a precisão e a convergência ao longo da profundidade. Os métodos de inversão foram capazes de estimar simultaneamente duas propriedades eletromagnéticas do modelo geológico: a permissividade elétrica e a condutividade elétrica. Os métodos de inversão alcançaram bons resultados quando executados simultaneamente sobre os dados sintéticos em 3 diferentes frequências. Exemplos de estimativas dos perfis de condutividade e permissividade elétricas unidimensionais são apresentados, com e sem inserção de ruídos nos dados. Os resultados indicam que a combinação dos métodos de inversão (ACO modificado e Quasi-Newton) pode fornecer bons resultados para as estimativas dos parâmetros físicos de meios geológicos e geotécnicos em meios rasos. Além disso, essa combinação de métodos de inversão abre novas perspectivas para o processamento de dados georradar multi-canais. As investigações também mostram que a simulação dos perfis sintéticos por meio do método FDTD pode propagar erros em espessuras maiores de subsolo, devido ao fato de ser um método explícito. / This work deals with simulation of a 1-D electromagnetic wave propagating into a geological structure of known physical properties (thickness, electric conductivity and permittivity), the record of the reflected EM wave (synthetic data), and the use of these records in inversion and optimization algorithms to estimate back the physical properties of the geological structure. A stratified and thin (< 1.2 m) geological structure was initially constructed in order to be scanned by high frequency EM waves (800, 1000, and 1200 MHz). The EM wave propagation is simulated by Maxwell equations through FDTD method. The Quasi-Newton inversion and Ant Colony Optimization methods were applied into synthetic data to estimate original physical parameters of each geological layer. Both methods were applied in order to increase precision and convergence along depth. These methods were able to simultaneously estimate two physical properties of the geological structure: electrical permittivity and conductivity. The methods showed good results when applied simultaneously upon synthetic data of all three frequencies. Electric permittivity and conductivity profiles are shown with and without noise in the data. The results indicate that combined inversion methods can show good results to estimate physical properties of thin geological and geothecnical structures. The combined inversion methods open new perspectives to multichannel GPR data processing. The investigations also show that simulation of synthetic data through FDTD can propagate errors when dealing with thick geological structures, due to the fact that FDTD is an explicit method.

High cost is a current problem with modern radar systems. Software-defined radios (SDRs) offer a possible solution for low-cost customizable radar systems. An SDR is a radio communi- cation system where, instead of the traditional radio components implemented in hardware, many of the components are implemented in software on a computer or embedded system. Although SDRs were originally designed for wireless communication systems, the firmware of an SDR can be configured into a radar system. With new companies entering the market, various types of low- cost SDRs have emerged. This thesis explores the use of a LimeSDR-Mini in a short-range radar through open software tools and custom code. The LimeSDR-Mini is successfully shown to detect targets at a short range. However, due to the instability of the LimeSDR-Mini, the consistent detection of a target is not possible. This thesis shows how the LimeSDR is characterized and how timing synchronization and instability issues are mitigated. The LimeSDR-Mini falls short of operating reliable in a radar system and other SDR boards need to be explored as viable options. Test setups using coaxial cables and test setups using antennas in an outdoor environment show the instability of the LimeSDR-Mini. The transmitter and the receiver are asynchronous. The timing difference varies slightly from run to run, which results in issues that are exacerbated in a short-range radar. The bleed-through signal is the signal leakage from the transmitter to the receiver. The bleed-through signal prevents the detection of targets at a short-range. Feed-through nulling is a signal processing technique used to eliminate the bleed-through signal so that short- range targets can be detected. The instability of the LimeSDR-Mini reduces the effectiveness of feed-through nulling techniques.

The theme of this dissertation is to advance safety hazard mitigation by detecting and characterizing hidden targets of concern. Ground-penetrating radar (GPR) is used to detect and characterize hidden targets that pose safety hazards at Earth's surface, within shallow soil, and within rock. The resulting images detect unexploded ordnance (UXO) and detect fractures that pose collapse hazards in a mine. Detecting and characterizing fractures and voids within rock prior to excavation can enable mitigation of mine collapse hazards. GPR data were acquired on the wall of a pillar in an underground mine. Strong radar reflections in the field data correlate with fractures and a cave exposed on the pillar walls. Pillar wall roughness was included in migration, a wavefield imaging algorithm, to quantitatively locate fractures and voids and map their spatial relationships within the rock. Quantifying the radar reflection amplitudes enabled mapping the distance between fracture walls. Detecting and characterizing UXO and landmines from a safe distance can enable de-mining. Migration was used to improve GPR imaging for unmanned aerial vehicle (UAV) data acquisitions. Existing algorithms were adapted for UAV flight irregularities and surface topography, and a new algorithm was developed that does not depend on the unknown soil wavespeed. Errors associated with wavespeed and raypath assumptions were quantified. The algorithms were tested with real and synthetic datasets. The improved and new algorithms are more successful than previous algorithms. To detect linear targets at all orientations, fully polarized GPR data are needed. Polarity combinations were investigated to optimize the detection of surface and subsurface small targets and linear targets. Scattering caused by topographic roughness is the primary shallow subsurface noise. For subsurface targets, detection is optimized by migration plus a polarity combination that captures all scattered energy. Strong reflection and scattering from the air-ground boundary can hide surface targets. Detection is optimized by removing the strong isotropic surface scattering, imaging targets by their anisotropic scattering. / Doctor of Philosophy / The theme of this dissertation is to advance safety hazard mitigation by detecting and characterizing hidden targets of concern. Ground-penetrating radar (GPR) is used to detect and characterize hidden targets that pose safety hazards at Earth's surface, within shallow soil, and within rock. The resulting images detect unexploded ordnance (UXO)/landmines and detect fractures that pose collapse hazards in a mine. Detecting and characterizing fractures and voids within rock prior to mining can enable mitigation of mine collapse hazards. GPR data were acquired on the wall of a pillar in an underground mine. Strong radar reflections in the field data correlate with fractures and a cave exposed on the pillar walls. Pillar wall roughness was included in migration, a wavefield imaging algorithm, to quantitatively locate fractures and voids and map their spatial relationships within the rock. Quantifying the radar reflection amplitudes enabled mapping the distance between fracture walls. Detecting and characterizing UXO, landmines from a safe distance can enable de-mining. Migration was used to improve GPR imaging for an unmanned aerial vehicle (drone) data acquisition. Existing algorithms were adapted for drone flight irregularities and surface topography, and a new algorithm was developed that does not depend on the unknown soil properties. Errors associated with the algorithms' assumptions were quantified. The algorithms were tested with real and computer-generated datasets. The improved and new algorithms are more successful than previous algorithms. To detect all targets regardless of their orientation, GPR data need to be acquired with antenna pointing in multiple directions (different polarities). Polarity combinations were investigated to optimize the detection of surface and subsurface small targets and linear targets. Scattering caused by topographic roughness is the primary shallow subsurface noise. For subsurface targets, detection is optimized by migration plus a polarity combination that captures all scattered energy. Strong radar reflection from the air-ground boundary can hide surface targets. Detection is optimized by removing the strong ground surface from the data, and imaging targets by differences in their radar scattering.

Near field imaging using microwave in medical applications has gain much attention recently as various researches show its high ability and accuracy in illuminating object comparing to the well-known screening tools such as Magnetic Resonance Imaging (MRI), digital mammography, ultrasound etc. This has encourage and motivate scientists continue to exploit the potential of microwave imaging so that a better and more powerful sensing tools can be developed. This thesis documents the development of antenna design for microwave imaging application such as breast cancer detection. The application is similar to the concept of Ground Penetrating Radar (GPR) but operating at higher frequency band. In these systems a short pulse is transmitted from an antenna to the medium and the backscattered response is investigated for diagnose. In order to accommodate such a short pulse, a very wideband antenna with a minimal internal reflection is required. Printed monopole and planar metal plate antenna is implemented to achieve the necessary operating wide bandwidth. The development of new compact printed planar metal plate ultra wide bandwidth antenna is presented. A generalized parametric study is carried out using two well-known software packages to achieve optimum antenna performance. The Prototype antennas are tested and analysed experimentally, in which a reasonable agreement was achieved with the simulations. The antennas present an excellent relative wide bandwidth of 67% with acceptable range of power gain between 3.5 to 7 dBi. A new compact size air-dielectric microstrip patch-antenna designs proposed for breast cancer detection are presented. The antennas consist of a radiating patch mounted on two vertical plates, fed by coaxial cable. The antennas show a wide bandwidth that were verified by the simulations and also confirmed experimentally. The prototype antennas show excellent performance in terms the input impedance and radiation performance over the target range bandwidth from 4 GHz to 8 GHz. A mono-static model with a homogeneous dielectric box having similar properties to human tissue is used to study the interaction of the antenna with tissue. The numerical results in terms the matching required of new optimised antennas were promising. An experimental setup of sensor array for early-stage breast-cancer detection is developed. The arrangement of two elements separated by short distance that confined equivalent medium of breast tissues were modelled and implemented. The operation performances due to several orientations of the antennas locations were performed to determine the sensitivity limits with and without small size equivalent cancer cells model. In addition, a resistively loaded bow tie antenna, intended for applications in breast cancer detection, is adaptively modified through modelling and genetic optimisation is presented. The required wideband operating characteristic is achieved through manipulating the resistive loading of the antenna structure, the number of wires, and their angular separation within the equivalent wire assembly. The results show an acceptable impedance bandwidth of 100.75 %, with a VSWR &lt; 2, over the interval from 3.3 GHz to 10.0 GHz. Feasibility studies were made on the antenna sensitivity for operation in a tissue equivalent dielectric medium. The simulated and measured results are all in close agreement.

A presente tese de doutoramento versa sobre o estudo das propriedades geotécnicas de solos residuais não saturados ocorrentes no município de Ubatuba, litoral norte do Estado de São Paulo. O estudo das propriedades geotécnicas foi realizado a partir de uma abordagem macro e micromorfológica, buscando estabelecer a relação existente entre a água presente no solo e sua variação sazonal. Os procedimentos e técnicas utilizadas permitiram diferenciar os horizontes de intemperismo quanto à microestrutura e porosidade, bem como caracterizar seus comportamentos físico e hídrico. Os ensaios geotécnicos, in situ e em laboratório, forneceram informações importantes sobre as propriedades dos materiais, tais como: índices físicos, distribuição granulométrica, limites de Atterberg, curvas de retenção de água e condutividade hidráulica saturada. Os dados registrados pelos equipamentos instalados no campo (pluviômetros, sensores GMS e FDR) permitiram avaliar a relação existente entre os eventos pluviométricos com a distribuição de umidade e sucção matricial ao longo dos perfis de alteração. Além disso, utilizouse metodologia geofísica de forma inédita na geotecnia brasileira para a determinação de perfis de umidade em solos residuais tropicais. Em termos gerais, os resultados obtidos permitiram estabelecer algumas relações entre as propriedades mecânicas e hidráulicas e as peculiaridades micromorfológicas dos horizontes de solos. Os resultados desse trabalho poderão ainda indicar áreas prioritárias para o monitoramento das variáveis climatológicas e geotécnicas, podendo fornecer apreciáveis subsídios para um futuro aperfeiçoamento dos parâmetros técnicos utilizados em Planos Preventivos de Defesa Civil (PPDC). / This doctoral thesis focusing on the study of the geotechnical properties of unsaturated residual soils occurring in the municipal district of Ubatuba, north coast of the State of Sao Paulo. The study of geotechnical properties was conducted from a macro and micromorfology approach, seeking to establish the relationship between the water present in the residual soil and its seasonal variations. The procedures and techniques enabled differentiate the horizons of weathering in relation to the microstructure and porosity, and characterize its physical-hydraulic behaviours. The laboratory and in situ geotechnical tests, provided important information about the properties of materials, such as: physical index, grain-size distribution, Atterberg limits, water retention curve and saturated hydraulic conductivity. The data registered by equipments installed in the field (pluviometers, FDR and GMS sensors) enabled to evaluate the relationship between the rainfall events with the distribution of moisture and suction matrix along the soil profiles. In addition, was used geophysical methodology in order unprecedented in brazilian geotechnical for the determination of moisture profiles in residual tropical soils. Overall, the results allowed establish some relationship between the hydraulic and mechanical properties and the micromorphological peculiarities of the horizons of soils. The results of this study may also indicate priority areas for monitoring the climatological and geotechnical variables, and may provide appreciable subsidies for a future improvement of the technical parameters used in Civil Defense Prevention Plans (PPDC).

The terms &lsquo / ground-probing radar&rsquo / , &lsquo / ground penetrating radar (GPR)&rsquo / , &lsquo / sub-surface radar&rsquo / or &lsquo / surface-penetrating radar (SPR)&rsquo / refer to various techniques for detecting and imaging of subsurface objects. Among those terms GPR is preferred and used more often. In this thesis, the depth and the position of the buried circular cylinder are determined by a GPR system which comprises of an optical fiber sensor (OFS). The system is a combination of OFS, GPR and optical communication link. In order to determine the depth and the position, first of all the electric field distribution at the OFS is obtained by integrating the Green&rsquo / s function over the induced current distribution. Those distributions are observed for different frequency and depth values. The voltages inside the distribution are measured by OFS. By changing the depth of the cylinder and the frequency of the system, various plots showing x axis displacement versus measured voltages are obtained. Those plots are related to interference fringe patterns. The position and the depth of the cylinder are obtained using interference fringe patterns. All of the studies mentioned are performed in MATLAB R2007b program. The noises of the system due to OFS are extracted using OPTIWAVE OPTISYSTEM 7.0 program. By adding those noises to the measured voltage values, the operating frequency of the system is observed.

A presente tese de doutoramento versa sobre o estudo das propriedades geotécnicas de solos residuais não saturados ocorrentes no município de Ubatuba, litoral norte do Estado de São Paulo. O estudo das propriedades geotécnicas foi realizado a partir de uma abordagem macro e micromorfológica, buscando estabelecer a relação existente entre a água presente no solo e sua variação sazonal. Os procedimentos e técnicas utilizadas permitiram diferenciar os horizontes de intemperismo quanto à microestrutura e porosidade, bem como caracterizar seus comportamentos físico e hídrico. Os ensaios geotécnicos, in situ e em laboratório, forneceram informações importantes sobre as propriedades dos materiais, tais como: índices físicos, distribuição granulométrica, limites de Atterberg, curvas de retenção de água e condutividade hidráulica saturada. Os dados registrados pelos equipamentos instalados no campo (pluviômetros, sensores GMS e FDR) permitiram avaliar a relação existente entre os eventos pluviométricos com a distribuição de umidade e sucção matricial ao longo dos perfis de alteração. Além disso, utilizouse metodologia geofísica de forma inédita na geotecnia brasileira para a determinação de perfis de umidade em solos residuais tropicais. Em termos gerais, os resultados obtidos permitiram estabelecer algumas relações entre as propriedades mecânicas e hidráulicas e as peculiaridades micromorfológicas dos horizontes de solos. Os resultados desse trabalho poderão ainda indicar áreas prioritárias para o monitoramento das variáveis climatológicas e geotécnicas, podendo fornecer apreciáveis subsídios para um futuro aperfeiçoamento dos parâmetros técnicos utilizados em Planos Preventivos de Defesa Civil (PPDC). / This doctoral thesis focusing on the study of the geotechnical properties of unsaturated residual soils occurring in the municipal district of Ubatuba, north coast of the State of Sao Paulo. The study of geotechnical properties was conducted from a macro and micromorfology approach, seeking to establish the relationship between the water present in the residual soil and its seasonal variations. The procedures and techniques enabled differentiate the horizons of weathering in relation to the microstructure and porosity, and characterize its physical-hydraulic behaviours. The laboratory and in situ geotechnical tests, provided important information about the properties of materials, such as: physical index, grain-size distribution, Atterberg limits, water retention curve and saturated hydraulic conductivity. The data registered by equipments installed in the field (pluviometers, FDR and GMS sensors) enabled to evaluate the relationship between the rainfall events with the distribution of moisture and suction matrix along the soil profiles. In addition, was used geophysical methodology in order unprecedented in brazilian geotechnical for the determination of moisture profiles in residual tropical soils. Overall, the results allowed establish some relationship between the hydraulic and mechanical properties and the micromorphological peculiarities of the horizons of soils. The results of this study may also indicate priority areas for monitoring the climatological and geotechnical variables, and may provide appreciable subsidies for a future improvement of the technical parameters used in Civil Defense Prevention Plans (PPDC).

UHF bands because more penetration can be achieved at low frequencies. Consequently, large antennas are required, which limits their use for small airborne applications. This thesis explores various GPR antenna designs for a bi-static system that are at least operational from 225 MHz to 255 MHz and suitable for small airborne applications. The 3D electromagnetic simulation software Ansys high-frequency structure simulator (HFSS) was used to simulate various sizes of strip dipole, triangular bowtie, half elliptical bowtie, and elliptical bowtie antennas. Several physical models were constructed to validate the return loss simulation results. Additionally, simulation data is included for a wire dipole and a helical antenna. The helical antenna proved to be too large for small airborne application, so focus was placed on the dipole and bowtie designs. The performance of the dipole and bowtie antenna models are compared by size, weight, return loss (𝑆11), peak gain, and the transmit-to-receive isolation. Out of the fourteen simulated models, twelve meet the bandwidth requirement with an average weight of 0.23 lbs. It is found that the strip dipole exhibited wider bandwidth characteristics than the triangular, elliptical, and half elliptical bowtie models, while maintaining similar weight and size. The smallest strip dipole model is 50 mm x 528 mm x 1 mm, weighs 0.17 lbs, and is operational from 225 MHz to 283 MHz. Two strip dipole test antennas were fabricated and tested. Test results confirm the simulation predictions.