Cette thèse se focalise sur le concept du radar MIMO co-localisé. L'acronyme MIMO -- pour Multiple-Input Multiple-Output -- indique l'utilisation de plusieurs émetteurs et de plusieurs récepteurs, tandis que le terme co-localisé signifie que ces éléments sont étroitement espacés. Chaque émetteur envoie une forme d'onde qui lui est propre : un radar MIMO émet donc simultanément un ensemble de signaux.Cette thèse a ainsi pour but d'établir une méthodologie permettant de générer cet ensemble de signaux, tout en respectant certaines contraintes opérationnelles. Cela nous permettra de déterminer les apports éventuels de ce radar. Nous nous sommes intéressés en particulier aux codes de phase, pour des raisons de couplage (qu'on peut traduire ici par la capacité, lors du traitement, à distinguer la position angulaire d'une cible de sa distance).La méthodologie proposée se synthétise simplement en une modélisation sous la forme d'un problème d'optimisation. Contrairement à la littérature et à des précédents résultats théoriques, nous avons décidé d'évaluer l'orthogonalité des signaux émis par le radar en différentes directions, et non l'orthogonalité des signaux élémentaires. Ce problème, plus réaliste, est malheureusement non-convexe et à grande échelle : un benchmark sur différentes méthodes d'optimisation nous a permis de constater l'efficacité des algorithmes basées sur le gradient.Optimiser cette orthogonalité sous-entend l'utilisation de filtres adaptés. Cependant, en pratique, le traitement radar s'effectue à l'aide de filtres désadaptés. Nous suggérons ainsi un problème d'optimisation jointe, permettant de générer de manière simultanée un ensemble de formes d'onde (pour le radar MIMO, entre autres) et les filtres désadaptés associés. Des simulations ont permis de montrer l'efficacité de la méthode. Celle-ci est en particulier préférable aux algorithmes cycliques habituellement utilisés. / This thesis deals with coherent MIMO radars. MIMO stands for Multiple-Input Multiple-Output, meaning that several transmitters and several receivers are used, closely-spaced in a coherent MIMO radar. Each transmitter has its own signal, providing waveform diversity. This thesis aims for defining a way to generate a set of sequences, specific for this radar, while satisfying practical constraints. It may help to determine the potential contribution of a MIMO radar. Only phase codes are concerned here, because they suffer less from the range/angle coupling effect.A simple framework is introduced, based on an optimisation problem.While literature often involves the orthogonality of the elementary signals (because of theoretical aspects), it is suggested to consider the orthogonality of signals from different directions of the surveillance space. Unfortunately, the obtained optimisation problem is non-convex and has a lot of variables. A benchmark on a simpler problem notifies us that gradient-based algorithms are surprisingly efficient.An optimisation of the correlation function corresponds to a processing with matched filters. However, in practice, mismatched filters are usually employed. A joint optimisation problem is suggested accordingly, in order to generate simultaneously a set of sequences (e.g. MIMO radar signals) and their associated mismatched filters. Obtained results are quite promising : as expected, a joint optimisation seems to perform better than a cyclic one, usually employed.

Micro-seismology is attracting more and more attention in the exploration seismology community. The main goal in micro-seismic imaging is to find the source location and the ignition time in order to track the fracture expansion, which will help engineers monitor the reservoirs. Conventional imaging methods work fine in this field but there are many limitations such as manual picking, incorrect migration velocity and low signal to noise ratio (S/N). In traditional surface survey imaging, full waveform inversion (FWI) is widely used. The FWI method updates the velocity model by minimizing the misfit between the observed data and the predicted data. Using FWI to locate and image microseismic events allows for an automatic process (free of picking) that utilizes the full wavefield. Use the FWI technique, and overcomes the difficulties of manual pickings and incorrect velocity model for migration. However, the technique of waveform inversion of micro-seismic events faces its own problems. There is significant nonlinearity due to the unknown source location (space) and function (time). We have developed a source independent FWI of micro-seismic events to simultaneously invert for the source image, source function and velocity model. It is based on convolving reference traces with the observed and modeled data to mitigate the effect of an unknown source ignition time. The adjoint-state method is used to derive the gradient for the source image, source function and velocity updates. To examine the accuracy of the inverted source image and velocity model the extended image for source wavelet in z-axis is extracted. Also the angle gather is calculated to check the applicability of the migration velocity. By inverting for the source image, source wavelet and the velocity model simultaneously, the proposed method produces good estimates of the source location, ignition time and the background velocity in the synthetic experiments with both parts of the Marmousi and the SEG Overthrust model. On the other hand, a new imaging condition of natural Green’s function has been implemented to mitigate the effect of the unknown velocity model. It is based on putting receivers in a horizontal well close to the micro-seismic events so that only a small part of the velocity model is required for the imaging. In order to focus the multi scattering energy to the source location, as well as to suppress the influence of the noise in the data, we introduced a new method to compensate the energy in the receiver wavefield. It is based on reflection waveform inversion (RWI) theory. We simply migrate for the scatters (reflectors) in the medium, and set the image as a secondary source to compensate for the multi scattering energy in the receiver wavefield. By applying the same imaging condition, the energy of those scattering events can be traced to the source location. Thus the source point has higher energy in the source image. A simple two-layer medium test demonstrates the features.

The main goal in seismic exploration is to identify locations of hydrocarbons reservoirs and give insights on where to drill new wells. Therefore, estimating an Earth model that represents the right physics of the Earth's subsurface is crucial in identifying these targets. Recent seismic data, with long offsets and wide azimuth features, are more sensitive to anisotropy. Accordingly, multiple anisotropic parameters need to be extracted from the recorded data on the surface to properly describe the model. I study the prospect of applying a scattering integral approach for multi-parameter inversion for a transversely isotropic model with a vertical axis of symmetry. I mainly analyze the sensitivity kernels to understand the sensitivity of seismic data to anisotropy parameters. Then, I use a frequency domain scattering integral approach to invert for the optimal parameterization. The scattering integral approach is based on the explicit computation of the sensitivity kernels. I present a new method to compute the traveltime sensitivity kernels for wave equation tomography using the unwrapped phase. I show that the new kernels are a better alternative to conventional cross-correlation/Rytov kernels. I also derive and analyze the sensitivity kernels for a transversely isotropic model with a vertical axis of symmetry. The kernels structure, for various opening/scattering angles, highlights the trade-off regions between the parameters. For a surface recorded data, I show that the normal move-out velocity vn, ƞ and δ parameterization is suitable for a simultaneous inversion of diving waves and reflections. Moreover, when seismic data is inverted hierarchically, the horizontal velocity vh, ƞ and ϵ is the parameterization with the least trade-off. In the frequency domain, the hierarchical inversion approach is naturally implemented using frequency continuation, which makes vh, ƞ and ϵ parameterization attractive. I formulate the multi-parameter inversion using the scattering integral method. Application to various synthetic and real data examples show accurate inversion results. I show that a good background ƞ model is required to accurately recover vh. For 3-D problems, I promote a hybrid approach, where efficient ray tracing is used to compute the sensitivity kernels. The proposed method highly reduces the computational cost.

The purpose of low probability of intercept (LPI) radar is, on top of the standard requirements on a radar, to remain undetected by hostile electronic warfare (EW) systems. This can be achieved primarily by reducing the amount of radiated power in any given direction at all times and is done by transmitting longer modulated pulses that can then be compressed digitally in order to retain range resolution. There are multiple different methods of performing pulse compression modifying either the phase or frequency of the transmitted waveform. Another method for attaining LPI properties of a radar is to avoid having a large main lobe in the transmit pattern and instead having lower gain patterns. This then results in a need for post-processing of these patterns by summation of weighted combination of these low gain patterns in order to reform the high gain patterns and thus retain angular resolution. In this work a number of pulse compression waveforms are analysed and compared using their ambiguity properties in order to ascertain which ones can be used in a radar system. They are then used in simulation with GO-CFAR detectors using a variety of analysis tools, specifically the short term Fourier transform (STFT), Wigner-Ville distribution (WVD), quadrature mirror filter bank (QMFB) and spectral correlation density (SCD). Their performance against the detector is based on the rate that the waveforms trigger an alarm and the lower the alarm rate the better the performance. The base reference in terms of performance for these evaluations was set as a triangular FMCW waveform. The results show that the polyphase coded waveforms have good radar and LPI properties in comparison to the FMCW. The frequency hopping codes showed good LPI properties with a large number of frequencies in the sequence but suffer from large ACF side lobes and poor Doppler tolerance. The best LPI results were achieved by a phase coded signal with a random order to its phase terms whilst still maintaining a perfect periodic autocorrelation function (PACF). Potential issues remain with high frequency out of band emission that could lead to a mismatch due to receiver bandpass filtering. The low gain patterns investigated were expanded to include two way patterns for a 2D array and array element tapering. The method works and can be further optimised in order to minimise emissions but adds a significant increase to integration times when the array size grows large.

The honey bee is one of the most important crop pollinating insects in the world. Researchers have recently identified a disease that has begun to impact the honey bee population. Colony Collapse Disorder results in the death of many bee colonies every year, but the cause for this remains unknown. Investigating the cause, harmonic radars are being considered to track the foraging patterns of honey bees. This research endeavors to find an optimized waveform for use in tracking foraging bees. Harmonic oscillators were developed for a transmit frequency of 1.2 GHz and various waveforms were tested against the oscillators. Ultimately, the waveform was found to be arbitrary. The amount of power that the harmonic oscillator receives is the determining factor. Given this, a general pulsed waveform can be developed that attempts to provide the maximum possible return for a predetermined maximum range of interest.

No description available.

No description available.

The objective of this thesis is to study the dynamic properties of laser diodes and the compensation for the nonlinearities of laser diodes based on the theory of Volterra series. In the first part of this thesis, an analytical expression in Volterra series is discussed to depict the nonlinear distortion of laser diodes up to the third order. The simulation results of this analytical method show that Volterra series model improves the accuracy of the description of the nonlinearity of laser diodes in comparison with small-signal analysis model. In the second part, the p^th-order inverse theory is introduced to compensate the lasers' nonlinear distortion. The compensation scheme is constructed and the simulation of the system is conducted in this thesis. The result shows that the laser nonlinear distortion can be compensated by using this technique. / Thesis / Master of Applied Science (MASc)

The Chesapeake Bay impact structure is one of the largest and most well preserved impact structures on Earth. It has a unique morphology composed of an inner crater penetrating crystalline basement surrounded by a wider crater in the overlying sediments. In 2004, the U.S. Geological Survey conducted a seismic survey with the goals of constraining crater structure and in support of the drilling of a borehole into the deepest part of the crater. Travel-time and waveform inversion were applied to the data to produce a high-resolution velocity model of the crater. Low-fold reflection processing was also applied. Northeast of the crystalline crater, undeformed, eastward-sloping crystalline basement is ~1.5 km deep. The edge of the inner crater is at ~ 15 km radius and slopes gradually down to a depth of 1.5 - 1.8 km. A central peak of 4-5 km radius rises to a depth of ~0.8 km. Basement velocity in the crystalline crater is much lower than undeformed basement, which suggests ~10% fracturing of the crater floor, and up to 20% fracturing of the central uplift. A basement uplift and lateral change of velocity, interpreted as the edge of the transient crater, occurs at a radius of ~ 11 km. Assuming a 22 km diameter transient crater, scaling laws predict a ~30 km diameter crater and central peak diameter of 8-10 km. This indicates that post-impact collapse processes that created the ~ 30 km diameter crystalline crater were unaffected by the much weaker rheology of the overlying sediments. / Master of Science

With the growing interest in software defined radios (SDRs), cognitive radios, the Joint Tactical Radio System (JTRS), and the Software Communication Architecture (SCA) comes the need for a rapid prototyping approach to radio design. In the past, radios have traditionally been designed to have a static implementation with the express goal of implementing a specific type of communication, such as 802.11b, CDMA voice communication, or just a simple FM tuner. However, when designing an SDR, the developer must not only be able to understand the radio engineering aspects of the design, but also be able to interface correctly with the underlying core software framework. This added software complexity, along with the general need for faster, more economical waveform development, illuminates the need for a rapid prototyping SDR development environment. This thesis takes a fresh look at the task of providing radio designers with a functional, straightforward design tool that enables the developer to focus more on the radio design than the tedious task of interacting with CORBA, IDL, and the SCA Core Framework. The design approach used to create such a tool is investigated along with an overview of general SDR concepts and an introduction to MPRG's open source SCA Core Framework, OSSIE. Discussion on the design methodology behind creating an SCA waveform is provided and the final result of this research, OSSIE Waveform Developer (OWD), is introduced and explored in detail. The code generated using OWD is detailed and waveform design approaches are presented with some suggested modifications. Finally, the improvements gained by using OSSIE Waveform Developer instead of the traditional approach of manually developing waveforms are presented. / Master of Science