Rozprawy doktorskie na temat „Wake Penetration”

Recent developments in cone penetration testing have resulted in the addition of both pore pressure measurements and seismometers. The seismometers allow shear wave velocity testing to be performed at designated intervals. Both of these additions were researched to improve their application and interpretation. The significant factors effecting the pore pressure generated during cone penetration tests are discussed. The importance of various factors is especially dependent upon permeability, strength, and stiffness. For all sands tested, pore pressures lower than static were recorded behind the tip and higher than static were recorded on the face of the cone. It is believed that the large compressive stresses on the cone face result in positive pore pressures. As the cone tip passes a soil element unloading and continued shearing generate pore pressures lower than static in all sands. The sign of this pore pressure (higher or lower than static) was therefore considered primarily a function of the test equipment. Pore pressure response and the rate of dissipation of excess pore pressures were found useful in distinguishing fine granular soils and explaining soil stratigraphy. In cohesive soils the details of pore pressure measurement were found to be important only in stiff soils. Pore pressures at all measurement locations were found to increase with soil strength in soft to firm clays but may be negative of static in very stiff clays. Pore pressures behind the cone tip were often negative of static in stiff clays. Measurement techniques were refined to improve the accuracy of downhole shear wave velocity measurements. Comparisons of downhole and crosshole measurements were made at three well documented sites validating the technique. At several sites it was found useful to consider the Gmax values determined from shear wave velocity and density to distinguish soil type. Gmax to cone resistance ratios were shown to vary systematically with cone resistance values in sands. A wide range in Gmax to cone resistance was observed in clays. The dependence of both cone penetration resistance and Gmax to increased stress level or overburden stress is discussed.
Applied Science, Faculty of
Civil Engineering, Department of
Graduate

This paper estimates the effects of imports on employment, wages, and the wage share in Austria for the period of 1990-2005 using panel data of manufacturing industry. Imports are disaggregated according to their origin and as final vs. intermediate imports. There is evidence of significant negative effects of imports on employment, wages and the wage share. Particularly workers in high skilled sectors experience negative effects. Offshoring to both Eastern Europe and the developed countries have a negative impact on employment, whereas offshoring to the East has a positive effect on wages, indicating the dominance of scope effects. (author´s abstract)
Series: Department of Economics Working Paper Series

Bibliography: pages 86-88.
Research into stepped frequency continuous wave ground penetrating radar (SFCW GPR) at UCT has been carried out since 1990. A first generation system comprising of Hewlett-Packard test equipment controlled by a PC was assembled. Cavity-backed log spiral antennas were designed and built by the University of Stellenbosch for the specific use of ground penetrating radar. Measurements with the first generation system proved the concept of SFCW GPR and thus a dedicated second generation system was planned. A SFCW GPR system was designed to replace the first generation system. Various designs for transmitter and receiver configurations were investigated and those found most suitable were used for the implementation. The SFCW radar consists of a wideband CW transmitter and a coherent receiver. A 300-1000 MHz transmitter was constructed using varactor-tuned oscillators as frequency sources. A double-sideband, low-IF receiver was constructed for the 300- 1000 MHz signal, to mix it to an IF of 10.7 MHz and I-Q demodulate it. The transmitter was found to operate according to specifications. The receiver was found to operate satisfactorily, but the dynamic range was less than expected. A limiting problem encountered in the first generation GPR was the large direct coupling signal from the transmitter into the receiver. This large signal reduced the effective receiver dynamic range. A method of cancelling this large direct coupling signal was implemented, using a bi-phase modulator to generate the cancelling signal in antiphase to the coupling signal. A 20 dB reduction in coupling was shown. The system was used to measure cable lengths to within the inherent accuracy of the system. A metal plate target was detected by the system feeding two antenna and a concrete floor was detected under 1 m of sand. It was thus shown that the SFCW system could be used as a second generation GPR.

This Master Thesis examines the possibility to apply a magnetometer developed by the Ångstöm space technology center to a small magnetic ground penetrating radar system with dimension in the order of one dm³. The magnetometer is broadband (DC-1GHz) and miniaturized. Loop antennas are used to transmit the signal.    A series of experiments have been performed in order to characterize the system, mainly examining the ability to determine distance to a target, using continuous sine wave signals and pulse trains. Standing wave patterns are formed between antenna and target and can be used for determining distance in the continuous case. When using a pulse train, the echo from the target could not be resolved using the current experiment set up, distance could therefore not be determined.

Knowledge and understanding of wave-induced dynamic soil response are of particular importance for marine geotechnical engineers involved in the design of the foundation for marine infrastructures such as breakwaters, platforms, pipelines and offshore wind turbines. In general, the nonlinearity of deformation occurring in loosely deposited sand foundations with a low relative density, and a low standard penetration test (SPT) value can presumably lead to the build-up of pore water pressures of foundation soil as a result of its plastic volumetric contraction under cyclic shearing. Liquefaction occurs when the accumulated pore pressure overcomes the overburden pressure, making the soil grains become unbound and completely free. The liquefied seabed soil will behave like a kind of heavy fluid without any shear resistance to supported submarine infrastructures, thusleading to catastrophic consequences such as the collapse of platforms and breakwaters, and sinking or floatation of submarine pipelines. In the present study, a theoretical framework is implemented to understand the relevant physical processes associated with wave-seabed-structure interaction, with a focus on the assessment of potential risks for failure of marine structures caused by liquefaction occurring in its loosely deposited sand foundations.
Thesis (PhD Doctorate)
Doctor of Philosophy (PhD)
Griffith School of Engineering
Science, Environment, Engineering and Technology
Full Text

Thesis (Ph.D)--Mechanical Engineering, Georgia Institute of Technology, 2010.
Committee Chair: Ume, Charles; Committee Member: Chen, Ye-Hwa; Committee Member: Michaels, Jennifer; Committee Member: Sadegh, Nader; Committee Member: Vachtsevanos, George. Part of the SMARTech Electronic Thesis and Dissertation Collection.

Thesis (Ph. D.)--Civil and Environmental Engineering, Georgia Institute of Technology, 2006.
Mayne, Paul W., Committee Chair ; Goldsman, David, Committee Member ; Lai, James, Committee Member ; Rix, Glenn J., Committee Member ; Santamarina, J. Carlos, Committee Member.

Includes bibliographical references.
A Ground Penetrating Radar (GPR) sensor is required to provide information that will allow the user to detect, classify and identify the target. This is an extremely tough requirement, especially when one considers the limited amount of information provided by most GPRs to accomplish this task. One way of increasing this information is to capture the complete scattering matrix of the received radar waveform. The objective of this thesis is to develop a signal processing technique to extract polarimetric feature vectors from Stepped Frequency Continuous Wave (SFGWV) GPR data. This was achieved by first developing an algorithm to extract the parameters from single polarization SFCW GPR data and then extending this algorithm to extract target features from fully polarimetric data. A model is required to enable the extraction of target parameters from raw radar data. A single polarization SFCW GPR model is developed based on the radar geometry and linear approximations to the wavenumber in a lossy medium. Assuming high operating frequencies and/or low conductive losses, the model is shown to be equivalent to the exponential model found in signal processing theory. A number of algorithms exist to extract the required target parameters from the measured data in a least squared sense. In this thesis the Matrix Pencil-of-Function Method is used. Numerical simulations are presented to show the performance of this algorithm for increasing model error. Simulations are also provided to compare the standard Inverse Discrete Fourier Transform (IDFT) with the algorithm presented in this thesis. The processing is applied to two sets of measured radar data using the radar developed in the thesis. The technique was able to locate the position of the scatterers for both sets of data, thus demonstrating the success of the algorithm on practical measurements. The single polarization model is extended to a fully polarimetric SFCW GPR model. The model is shown to relate to the multi-dimensional exponential signal processing model, given certain assumptions about the target scattering damping factor. The multi-snapshot Matrix Pencil-of-Function Method is used to extract the scattering matrix parameters from the raw polarimetric stepped frequency data. Those Huynen target parameters that are independent of the properties of the medium, are extracted from the estimated scattering matrices. Simulations are performed to examine the performance of the algorithm for increasing conductive and dielectric losses. The algorithm is also applied to measured data for a number of targets buried a few centimeters below the ground surface, with promising results. Finally, the thesis describes the design and development of a low cost, compact and low power SFCW GPR system. It addresses both the philosophy as well as the technology that was used to develop a 200 - 1600 MHz and a 1 - 2 GHz system. The system is built around a dual synthesizer heterodyne architecture with a single intermediate frequency stage and a novel coherent demodulator system - with a single reference source. Comparison of the radar system with a commercial impulse system, shows that the results are of a similar quality. Further measurements demonstrate the radar performance for different field test cases, including the mapping of the bottom of an outdoor test site down to 1.6 m.

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.

Enhanced Integration of Shear Wave Velocity Profiling in Direct-Push Site Characterization Systems Alexander V. McGillivray 370 Pages Directed by Dr. Paul W. Mayne Shear wave velocity (VS) is a fundamental property of soils directly related to the shear stiffness at small-strains. Therefore, VS should be a routine measurement made during everyday site characterization. There are several lab and field methods for measuring VS, but the seismic piezocone penetration test (SCPTu) and the seismic dilatometer test (SDMT) are the most efficient means for profiling the small-strain stiffness in addition to evaluating large-strain strength, as well as providing evaluations of the geostratigraphy, stress state, and permeability, all within a single sounding. Although the CPT and DMT have been in use for over three decades in the USA, they are only recently becoming commonplace on small-, medium-, and large-size projects as more organizations begin to realize their benefits. Regrettably, the SCPTu and the SDMT are lagging slightly behind their non-seismic counterparts in popularity, in part because the geophysics component of the tests has not been updated during the 25 years since the tests were envisioned. The VS measurement component is inefficient and not cost effective for routine use. The purpose of this research is to remove the barriers to seismic testing during direct-push site characterization with SCPTu and SDMT. A continuous-push seismic system has been developed to improve the integration of VS measurements with SCPTu and SDMT, allowing VS to be measured during penetration without stopping the progress of the probe. A new type of portable automated seismic source, given the name RotoSeis, was created to generate repeated hammer strikes at regularly spaced time intervals. A true-interval biaxial seismic probe and an automated data acquisition system were also developed to capture the shear waves. By not limiting VS measurement to pauses in penetration during rod breaks, it is possible to make overlapping VS interval measurements. This new method, termed frequent-interval, increases the depth resolution of the VS profile to be more compatible with the depth intervals of the near-continuous non-seismic measurements of the SCPTu and the SDMT.

This thesis presents an empirical study comparing the ability of multi-static and bi-static, handheld, ground penetrating radar (GPR) systems, using full wave inversion (FWI), to determine the properties of buried anti-personnel (AP) landmines. A major problem associated with humanitarian demining is the occurrence of many false positives during clearance operations. Therefore, a reduction of the false alarm rate (FAR) and/or increasing the probability of detection (POD) is a key research and technical objective. Sensor fusion has emerged as a technique that promises to significantly enhance landmine detection. This study considers a handheld, combined metal detector (MD) and GPR device, and quantifies the advantages of the use of antenna arrays. During demining operations with such systems, possible targets are detected using the MD and further categorised using the GPR, possibly excluding false positives. A system using FWI imaging techniques to estimate the subsurface parameters is considered in this work. A previous study of multi-static GPR FWI used simplistic, 2D far-field propagation models, despite the targets being 3D and within the near field. This novel study uses full 3D electromagnetic (EM) wave simulation of the antenna arrays and propagation through the air and ground. Full EM simulation allows the sensitivity of radio measurements to landmine characteristics to be determined. The number and configuration of antenna elements are very important and must be optimised, contrary to the 2D sensitivity studies in (Watson, Lionheart 2014, Watson 2016) which conclude that the degree (number of elements) of the multi-static system is not critical. A novel sensitivity analysis for tilted handheld GPR antennas is used to demonstrate the positive impact of tilted antenna orientation on detection performance. A time domain GPR and A-scan data, consistent with a commercial handheld system, the MINEHOUND, is used throughout the simulated experiments which are based on synthetic GPR measurements. Finally, this thesis introduces a novel method of optimising the FWI solution through feature extraction or estimation of the internal air void typically present in pressure activated mines, to distinguish mines from non-mine targets and reduce the incidence of false positives.

Bibliography: leaves 64-68.
A device has been designed that cancels the leakage signal between the transmit and receive antenna in a Stepped Frequency Continuous Wave Ground Penetrating Radar. The front end of the radar operates at high signal levels and, as a result, a large signal is coupled directly from the transmit to the receive antenna. This signal uses a signiï¬ cant part of the dynamic range of the data-capturing device, an analogue-to-digital converter (ADC). The objective of this cancellation is thus to increase the effective instantaneous dynamic range of the radar system. Simulations show that 10-bit amplitude and phase resolution in the digital cancellation circuit would achieve maximum cancellation in the presence of phase noise and other sources of error. This result is conï¬ rmed when the hardware is tested. The device was constructed and operates as intended. Tests show that cancellation exceeding 53dBm is possible through careful calibration. It was concluded that the device could successfully be integrated into the SFCW GPR and that it would achieve an increase in the instantaneous dynamic range.

Evaluations of stress history and the geostatic state of stress of soils are ascertained on the basis of field geophysical measurements that provide paired complementary types of shear waves. It is well-established that multiple types of shear waves occur in the ground due to their directional and polarization properties. The shear wave velocity (Vs) provides the magnitude of small strain stiffness (G0) which depends on effective stress, void ratio, stress history, and other factors (cementation, age, saturation). Herein, this study examines a hierarchy of shear wave modes with different directions of propagation and particle motion from in-situ geophysical tests (HH, VH, and HV) and laboratory bender element data. A special compiled database from well-documented worldwide sites is assembled where full profiles of stress state, stress history, and several paired modes of Vs profiles have been obtained from crosshole tests (CHT), downhole tests (DHT), and rotary crosshole (RCHT). Reference profiles of the lateral stress coefficient (K0) are available from direct in-situ measurements (self-boring pressuremeter, hydrofracture, and push-in spade cells). Stress history is documented in terms of yield stress ratio (YSR) from consolidation testing and careful engineering geology studies. A methodology is developed that relates both the YSR and K0 to stiffness ratios obtained from directional shear wave velocities. In further efforts, means to extract reliable shear wave profiles from continuous downhole testing via a new GT autosource and seismic piezocone testing are outlined and applied to results from three test sites in Windsor/VA, Norfolk/VA, and Richmond/BC. A driving impetus to this research involves the geologic conditions at the US Dept. of Energy's Savannah River Site (SRS) in South Carolina. Here, the overburden soils in the upper 60 m depths consist of very old Miocene and Eocene sediments, primarily layered deposits of sands, clayey sands, silty sands, and interbedded clays which exhibit an apparent and unusual stress history profile. Special geologic conditions include the dissolutioning of old calcareous sediments (Santee Formation) at depths of 40 to 50 m below grade, similar to karstic limestone deposits. As a consequence, caves, voids, and infilled soft soil zones occur within the soil matrix at these elevations, probably resulting in localized collapse of the overlying soil column. Based on conventional laboratory and in-situ test data conducted during geotechnical investigations at SRS, available interpretative relationships for assessing the soil stress history and geostatic stress states show scattered and inconsistent results. Complications abound in the systematic assessments of these geomaterials due to effects of very old ageing, cementation, desiccation, and diagenesis, as evidenced by unusual in-situ shear wave velocity profiles that decrease in magnitude with depth, as measured by CHT and DHT. Based on the findings of this study, it is recommended that a new set of shear wave velocity measurements be made at SRS to obtain HH waves (and complementary VH waves) needed for an independent assessment of YSR in the upper soil column.

Two new stepped-frequency continuous wave (SFCW) radar sensor prototypes, based on a coherent super-heterodyne scheme, have been developed using Microwave Integrated Circuits (MICs) and Monolithic Millimeter-Wave Integrated Circuits (MMICs) for various surface and subsurface applications, such as profiling the surface and subsurface of pavements, detecting and localizing small buried Anti-Personnel (AP) mines and measuring the liquid level in a tank. These sensors meet the critical requirements for subsurface and surface measurements including small size, light weight, good accuracy, fine resolution and deep penetration. In addition, two novel wideband microstrip quasi-TEM horn antennae that are capable of integration with a seamless connection have also been designed. Finally, a simple signal processing algorithm, aimed to acquire the in-phase (I) and quadrature (Q) components and to compensate for the I/Q errors, was developed using LabView. The first of the two prototype sensors, named as the microwave SFCW radar sensor operating from 0.6-5.6-GHz, is primarily utilized for assessing the subsurface of pavements. The measured thicknesses of the asphalt and base layers of a pavement sample were very much in agreement with the actual data with less than 0.1-inch error. The measured results on the actual roads showed that the sensor accurately detects the 5-inch asphalt layer of the pavement with a minimal error of 0.25 inches. This sensor represents the first SFCW radar sensor operating from 0.6-5.6-GHz. The other sensor, named as the millimeter-wave SFCW radar sensor, operates in the 29.72-35.7-GHz range. Measurements were performed to verify its feasibility as a surface and sub-surface sensor. The measurement results showed that the sensor has a lateral resolution of 1 inch and a good accuracy in the vertical direction with less than  0.04-inch error. The sensor successfully detected and located AP mines of small sizes buried under the surface of sand with less than 0.75 and 0.08 inches of error in the lateral and vertical directions, respectively. In addition, it also verified that the vertical resolution is not greater than 0.75 inches. This sensor is claimed as the first Ka-band millimeter-wave SFCW radar sensor ever developed for surface and subsurface sensing applications.

This project investigates the application of ground penetrating radar (GPR), a remote sensing geophysical survey method, to the archaeological investigation of earthen architecture on the Northwest Coast of North America. The objective of this thesis is to assess the ability of GPR to detect and distinguish between architectural features within an earthen matrix, and to understand the limitations and uncertainties of the method in this and similar contexts. This thesis also assesses the ability of GPR to provide data that are able to contribute to broad anthropological questions of demographic change and socio-political complexity. GPR was used at Welqámex (DiRi-15), a Stó:lō-Coast Salish settlement near Hope, British Columbia, to collect nearly 1,000 m2 of data over a minimum of 11 structures. GPR data were analysed with comparison to surface and subsurface data from Welqámex, including excavation data collected prior to and following GPR survey. The survey identified 157 anomalies that may be useful in guiding future excavations. Direct comparisons of GPR reflection profiles and amplitude slices with excavation stratigraphic profiles and plan views indicate that GPR is moderately successful in detecting sqémél floors, s’iltexwáwtxw floors, and pit features larger than 15 cm in diameter, but is not successful in detecting post and stake mold features larger than 15 cm in diameter, hearth features, and structure boundaries. The anomalies produced from these features, however, are not easily distinguished from one another or from other natural and archaeological features. The results suggest that while GPR is able to identify anomalies that may be useful in guiding archaeological excavation, it is at this time not an ideal method for addressing broader anthropological questions on its own.

The present thesis is devoted to the investigation of novel superconductors by phase-sensitive terahertz transmission and infrared to ultraviolet spectroscopy. In particular, a nominally undoped Pr2CuO4 superconducting thin film, an FeTe0.5Se0.5 thin film, and a LiFeAs single crystal have been investigated. The emphasis is on the low-frequency part of the optical spectrum (i.e., the terahertz and infrared spectrum), as the goal of the study was to shed light on the size and symmetry of the superconducting gaps and also to determine the temperature dependences as well as the absolute values of the penetration depth, which are key input parameters for models applicable for new superconductors. In addition, niobium has been investigated as a reference, so as to see what is expected from conventional superconductors and to clarify the electrodynamics of niobium. A superconducting Nb thin film with Tc of 8.04 K has been investigated by backward wave oscillator-based (BWO-based) and time-domain terahertz (TDT) spectrometers in the frequency range between 4 and 100 cm−1 for temperatures ranging from 2 to 10 K. From these measurements an energy gap of 22.50 cm−1 = 2.79 meV = 4.02kBTc have been determined. The optical conductivity below Tc could nicely be described by calculations according to the Eliashberg theory, with the electron-phonon interaction evaluated from tunneling measurements. Absolute values of the penetration depth have been calculated from phase-sensitive terahertz measurements. The zero-temperature limit of at T = 0 is found to be 115 ± 5 nm. From this value, a London penetration depth of 43 ± 2 nm has been obtained. The overall temperature dependence of the penetration depth follows a behavior typical for conventional s-wave superconductors. A superconducting Pr2CuO4 film with T0 structure and Tc of 27 K has been investigated by use of optical methods in a wide frequency (5 – 55000 cm−1) and temperature (2 – 300 K) range. A Drude-like peak centered at zero frequency is observed in the optical conductivity below 150 K, above which it shifts to finite frequencies. The detailed analysis of the low-frequency conductivity reveals that the Drude peak and a far-infrared (FIR) peak centered at about 300 cm−1 persist at all temperatures. The FIR spectral weight is found to grow at the expense of the Drude spectral weight with increasing temperature. Absolute values of the penetration depth have been obtained from temperature and frequency-dependent measurements. The zero-temperature limit of is estimated to be 1600 ± 100 nm. The overall temperature dependence of follows a behaviour typical for cuprate superconductors. However, a closer look at the penetration depth at T 12 K reveals a flattening in the temperature dependence. A superconducting FeTe0.5Se0.5 thin film with Tc = 19 K has been investigated using a combination of BWO and TDT spectroscopy in the frequency range 4 - 80 cm−1 and between 3 and 150 K. From such measurements, a superconducting energy gap of 30 cm−1, representing a coupling strength = 2.27, is observed. Further, the penetration depth has been derived from the temperature dependence of the imaginary part of complex conductivity with the penetration depth = 530 ± 10 nm at lowest measured temperature. The temperature-dependent normalized superfluid density, just as is the case with most iron-based superconductors, could nicely be described by the so-called two-gap gamma model. Finally, a superconducting LiFeAs single crystal with Tc = 18 K has been investigated by optical spectroscopy in the frequency range 15 - 55000 cm−1 between 5 and 300 K. From these measurements, no clear signature of the superconducting energy-gap opening could be identified in spite of the spectral weight been suppressed in the infrared frequency regime below Tc. This indicates that LiFeAs single crystal is in a clean limit. With the aid of the Ferrell-Glover-Tinkham (FGT) sum rule, an absolute penetration depth of 215 nm has been calculated from the missing area at 5 K.

Since 1976, the General Directorate of Disaster Affairs of Turkey has deployed several strong-motion accelerographs at selected sites. Within the framework of the project entitled Compilation of National Strong Ground Motion Database in Accordance with International Standards, initiated in 2006, site conditions at a total of 153 strong-motion sites were investigated within the uppermost 30 m depth through boreholes including Standard Penetration Testing and surface seismics by means of Multi-channel Analysis of Surface Waves (MASW). In this study, firstly, the assessment of the site characterization was held by making use NEHRP Provisions, EC-8 and Turkish Seismic Design Code. The corrected penetration resistances are calculated and observed how it affects the classification. In addition, the consistency of site classes obtained from either penetration resistance or shear wave velocity criteria is examined. Also the consistency of the boundaries of the site classes in terms of shear wave velocity and penetration resistance data pairs are investigated. Secondly, the liquefaction potential of these sites is examined. Thirdly and finally, the shear wave velocity profiles obtained from MASW technique are contrasted to other seismic tests.

A multi-scale high-resolution geophysical and geological study was conducted in the Sevier Desert, central Utah, found within the Colorado Plateau-Basin and Range Transition Zone. The region is marked by with Quaternary volcanics and faulting as young as 660 yr B.P., with many fault scarps thought to have the potential for 7+ magnitude earthquakes. Three locations within the Sevier Desert which represent three different tectonic expressions of possible faulting at the surface were selected. These include a location found within surface sedimentation, a location with surface sedimentation and sub-surface basalts and a location with basalts, at the surface with very limited sedimentation. A suite of geophysical data were obtained including the use of P-wave, SH-wave, ground-penetrating radar (GPR). Auger holes, microprobe glass analysis, and mapping information were also completed in order to constrain and gain a more complete understanding of the sub-surface structure. These data were used to determine if there are sub-surface expressions of the possible surface scarps and if all the faults within the fault zone have the same structural style. The possible surface fault expressions were found to be connected to sub-surface fault expressions but with differing results within both sediments and basalts. Our data show that a multi-scale approach is needed to obtain a complete view of tectonic activity. The area faulting in the Sevier Desert penetrates at depth involving multiple complex styles that include some faulting that cuts recent lava flows and some that do not. The evidence also indicates that in at least some area faulting was episodic and others may be single events having implications on level of activity and hazard.

Ce travail présente les développements récents réalisés sur le pénétromètre PANDA 3. Il s'agit d'un pénétromètre dynamique instrumenté qui permet à partir de la mesure puis du découplage des ondes créées par l’impact sur l’appareil, d’obtenir pour chaque coup une courbe charge-enfoncement σp-sp du sol ausculté. L’exploitation de cette courbe permet de déterminer des paramètres de résistance (résistance de pointe qd), de déformation (module dynamique Ed P3), des caractéristiques d’amortissement Js et de célérité d'ondes (CsP3 et CpP3) des sols auscultés en fonction de la profondeur tout au long du sondage. Cependant, et bien que la méthode soit très intéressante, celle-ci est restée au stade d’un prototype de laboratoire. Il est donc nécessaire de réaliser une étude plus approfondie sur l'essai lui-même et sur l’information contenue dans la courbe σp-sp en vue de fiabiliser la mesure et d’améliorer son exploitation. Dans un premier temps, nous présentons un bref aperçu sur les techniques de reconnaissance géotechnique et plus particulièrement celle des essais de pénétration dynamique. Le principe général du PANDA 3 est également présenté. La deuxième partie est consacrée au développement d'un nouveau prototype de pénétromètre PANDA 3. Ce développement s’appuie sur plusieurs études visant à valider la qualité des informations recueillies, leur bonne reproductibilité et le traitement des signaux d’acquisition. De même, un modèle numérique discret du battage pénétrométrique développé à l’aide du logiciel Particle Flow Code (Itasca) est présenté permettant de valider la technique de mesure. La troisième partie traite d'une étude comparative des résultats obtenus avec le PANDA 3 et d’autres techniques d'auscultation in situ afin de valider les résultats obtenus et l’utilisation de l’appareil dans des conditions réelles. Par ailleurs l’extension de cette technique de mesure au cas des pénétromètres lourds est appliquée dans le but de mesurer l'énergie transmise et d’étalonner le système de battage. Enfin, la dernière partie est consacrée à l'interprétation et l'exploitation des signaux en pointe afin d'affiner le modèle d'interprétation de la courbe charge-enfoncement. L'analyse de l'ensemble des signaux enregistrés au laboratoire a permis d’approcher une méthodologie d'exploitation de la courbe. L'application de la méthode proposée a été réalisée pour différents sols aussi bien au laboratoire que sur le terrain. Les résultats obtenus ont été confrontés avec d'autres types d’essais
This work presents the recent developments made on the penetrometer PANDA®3. The instrumented dynamic penetrometer allowing, from the measurement and the decoupling of waves created by the impact, to obtain the load-penetration curve σp-sp of the soil. The exploitation of this curve allows determining the failure parameter (tip resistance qd), deformation (dynamic modulus EdP3), damping characteristics (Js) and wave speed (CsP3 and CpP3) of the investigated soil according to depth all along the sounding. However, although the proposed method is very interesting, it has remained at the stage of a laboratory prototype. It is therefore necessary to conduct a more thorough study of the test itself and the information provided from the σp-sp curve in order to obtain reliable measurement and improve their exploitation. First of all, we are presenting a brief overview of the geotechnical in-situ testing particularly that of the dynamic penetration tests as well as the general principle of PANDA 3 is presented. The second part is devoted to the development of a new prototype of the PANDA 3 penetrometer. This development is based on several studies aiming at validating the quality of the information, good reproducibility and treatment of the acquisition signals. Similarly, a discrete numerical model of the penetrometer developed using the Particle Flow Code software (Itasca) is presented to validate the measurement technique. The third part deals with a comparative study of the results obtained with the PANDA 3 and other in situ investigation techniques to validate the obtained results and the use of the device in real conditions. Moreover, the extension of this measurement technique in the case of heavy penetrometer is applied in order to measure the transmitted energy and to calibrate the driving system. The last part is devoted to refining the interpretation and exploitation of the load-penetration curve. The analysis of all the signals recorded in the laboratory allowed to approach a methodology of curve exploitation. The application of the proposed method was carried out for different soils both in the laboratory and on field. The results were confronted with other types of tests

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.

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.

Cette thèse présente une étude expérimentale de l'interaction entre la supraconductivité et une onde de densité de charge (ODC). Dans la théorie standard, la température critique d’un matériau supraconducteur est favorisée principalement par deux paramètres : une grande densité d’états au niveau de Fermi (nF), et un fort couplage électron-phonon. Cependant, un fort couplage électron-phonon favorise aussi l’apparition d’une ODC, ce qui réduit nF et rivalise ainsi avec la supraconductivité.Notre démarche a consisté à étudier deux composés où supraconductivité et ODC coexistent, et dans lesquels on peut faire disparaître l’ODC grâce à un paramètre externe : pression ou substitution. Le premier composé, 2H-NbSe2, présente une ODC en dessous de 33 K à pression ambiante. Celle-ci coexiste avec la supraconductivité en dessous de 7 K. Sous pression, l’ODC disparaît au-dessus de 4.6 GPa, sans que la température critique varie notablement. L’ODC disparaît aussi en remplaçant le sélénium par du soufre : 2H-NbS2 est ainsi un supraconducteur sans ODC (Tc = 6 K), et peut donc servir de composé témoin pour une étude comparative. Dans le second composé, Lu5Ir4Si10, une ODC est présente en dessous de 77 K à pression ambiante. Celle-ci disparaît sous pression au-dessus de 2 GPa, tandis que la température critique saute simultanément de 4 à 9 K. Pour étudier ces composés, j’ai utilisé trois techniques expérimentales : la mesure de la dispersion des phonons à basse température (300-2 K) et sous pression (0-16 GPa) par diffusion inélastique des rayons X, la mesure de la dépendance en température de la longueur de pénétration magnétique grâce à un oscillateur à diode tunnel et la mesure des champs critiques via des microsondes Hall.Dans la première partie, je présente la dépendance en température de la dispersion des phonons dans 2H-NbS2. Nous observons la présence d’un phonon mou dont l’énergie reste toujours positive, même extrapolée à température nulle. Ce composé est ainsi à la limite d'une instabilité ODC. De plus, nous montrons qu’il est relativement unique, car seuls les effets anharmoniques empêchent l’amollissement complet des phonons. Je présente ensuite la dépendance en température et en pression de la dispersion des phonons dans 2H-NbSe2. Ces expériences montrent qu’un mode de phonon mou persiste jusqu’à 16 GPa, même quand l'état à température nulle n'est pas l’ODC. La dépendance en température de ce phonon mou est alors similaire à celle de 2H-NbS2. Dans les deux composés, ces phonons mous semblent liés à la présence d'un couplage électron-phonon à la fois fort et anisotrope. Nous suggérons qu’il s’agit d’un élément essentiel pour expliquer leurs propriétés supraconductrices.Dans la seconde partie, je mesure l'anisotropie et la dépendance en température de la longueur de pénétration magnétique dans l’état supraconducteur de 2H-NbS2 et Lu5Ir4Si10. La dépendance en température de la densité superfluide dans 2H-NbS2 confirme la présence d'un gap supraconducteur réduit dont l'amplitude est très proche de celle mesurée dans 2H-NbSe2. Les phonons mous et le gap réduit étant présents dans 2H-NbS2 et 2H-NbSe2, nous prouvons expérimentalement qu'il faut raisonner en termes de renforcement de la supraconductivité par les phonons mous plutôt qu'en termes d’interaction avec l'état fondamental (ODC ou métal). Nous proposons que ce renforcement soit lié à l'anisotropie du couplage électron-phonon.En revanche, cet effet n’est pas général aux composés où supraconductivité et ODC coexistent. Les propriétés supraconductrices de Lu5Ir4Si10 sont en effet bien décrites par le modèle BCS couplage faible. Ceci est peut être lié aux caractéristiques de l’ODC : la présence d’une hystérésis montre que la transition ODC est du premier ordre. D’autre part, les mesures de diffraction X sous pression et à basse température révèlent que cette ODC est multiple : en plus de la périodicité 1/7, nous observons une seconde périodicité de 1/20
This thesis presents an experimental study of the interaction between superconductivity and a charge density wave (CDW). In the standard theory, the critical temperature of a superconductor is principally enhanced by two parameters: a large density of states at the Fermi level (nF) and a strong electron-phonon coupling. However, a strong electron-phonon coupling also favors the appearance of a CDW, which reduces nF and therefore competes with superconductivity.Our strategy was to study two compounds in which superconductivity and CDW coexist, and in which the CDW can be suppressed through an external parameter: pressure or substitution. The first compound is 2H-NbSe2, it presents a CDW below 33 K at ambient pressure. This CDW coexists with superconductivity below 7 K. Under pressure, the CDW disappears above 4.6 GPa, meanwhile the critical temperature slowly changes. The CDW also disappears when replacing selenium by sulfur: 2H-NbS2 is a superconductor without CDW (Tc=6 K), it can therefore serve as a “test compound” for a comparative study. The second compound is Lu5Ir4Si10, it presents a CDW below 77 K at ambient pressure. Under pressure, this CDW disappears above 2 GPa, meanwhile the critical temperature abruptly jumps from 4 to 9 K.For this study, I used three experimental techniques: inelastic x-ray scattering at low temperature (300-2 K) and under pressure (0-16 GPa) to measure the dispersion of phonons, a tunnel diode oscillator to measure the temperature dependence of the magnetic penetration depth, and Hall microprobes to measure the first and second critical fields. In the first part, I present the temperature dependence of the phonon dispersion in 2H-NbS2. We observe a soft phonon that always remains at positive energies, even extrapolated to zero temperature. Thus, this compound is on the verge of CDW instability. It is also relatively unique, since we show anharmonicity is the only effect that prevents the complete softening of the phonons.Then I present the temperature and pressure dependence of the phonon dispersion in 2H-NbSe2. These experiments show that a soft phonon persists up to 16 GPa, even if the ground state is not a CDW. The temperature dependence of this soft phonon is then similar to that of 2H-NbS2. In both compounds, these soft modes seem to be related to the strength and anisotropy of the electron-phonon coupling. We suggest this is a fundamental element to explain their superconducting properties.In the second part, I measure the anisotropy and temperature dependence of the magnetic penetration depth in the superconducting state of 2H-NbS2 and Lu5Ir4Si10. The temperature dependence of the superfluid density in 2H-NbS2 confirms the presence of a reduced superconducting gap. Its amplitude is very similar to the one measured in 2H-NbSe2. The soft modes and the reduced gap being present in both 2H-NbSe2 and 2H-NbS2, we prove experimentally that the enhancement of superconductivity is related to the soft modes rather than to the nature of the ground state (CDW or metal). We suggest this enhancement is due the anisotropy of the electron-phonon coupling.However, this effect is not general to all compounds where superconductivity and CDW coexist. The superconducting properties of Lu5Ir4Si10 are indeed well fitted by the BCS model in the weak coupling limit. This may be related to the characteristics of the CDW: the presence of hysteresis shows that the CDW transition is first order. In addition, under pressure and at low temperature, x-ray diffraction measurements indicate that the CDW is multiple: aside from the periodicity of 1/7, we observe a second periodicity of 1/20

Observations of modern carbonate depositional environments and their accompanying depositional models have been used for decades in the reconstruction and interpretation of ancient carbonate depositional environments and stratigraphic successions. While these Holocene models are necessary for interpreting their more ancient counterparts, they inherently exclude important factors related to the erosion, diagenesis, and ultimate preservation of sediments and sedimentary structures that are ubiquitous in shallow marine carbonate environments. Andros Island, Bahamas is an ideal location to examine the validity of Holocene conceptual models, where geologically young (Late Pleistocene) limestones can be studied immediately adjacent to their well-documented modern equivalents. For this study, two 3D ground-penetrating radar (GPR) datasets (200 MHz and 400 MHz) were collected at a schoolyard in northwest Andros. These surveys reveal the geometries and internal characteristics of a peloidal-oolitic sand wave and tidal channel in unprecedented detail. These two prominent features are underlain by low-energy lagoonal wackestones and packstones, and are bordered laterally to the northwest by wackestones-packstones intermixed with thin sheets of peloidaloolitic grainstone. A deeper radar surface is observed at approximately 6 m depth dipping gently to the west, and is interpreted to be a karstified exposure surface delineating the base of a complete depositional sequence. Interpretation of the 3D radar volumes is enhanced and constrained by data from three cores drilled through the crest and toe of the sand wave, and through the tidal channel. This study is the first of its kind to capture the complex heterogeneity of a carbonate depositional package in three dimensions, where various depositional environments, sedimentary structures, and textures (mudstone to grainstone) have been preserved within a small volume.The results from this study suggest that the degree of vertical and lateral heterogeneity in preserved carbonate successions is often more complex than what can be observed in modern depositional environments, where sediments can generally only be observed in two dimensions, at an instant in time. Data from this study demonstrate the value of using two overlapping GPR datasets at differing resolutions to image the internal characteristics of a complete carbonate depositional package in three dimensions. From these datasets, a depositional model similar to other Holocene and Pleistocene carbonate depositional models is derived.

O ensaio de simples reconhecimento do solo (SPT) e usualmente utilizado para estimar a resistência do solo, atraves do índice NSPT, que representa a resistencia a penetração dinâmica do amostrador no solo. Esse índice é usado diretamente em correlações empíricas ou semi-empíricas na determinação da capacidade de suporte e recalque das fundações. Entretanto, esse ensaio tem sido alvo de críticas devido a utilização dessas correlações empíricas, geralmente baseadas em observações práticas e sem nenhum fundamento científico. Críticas também estão relacionadas com a dispersão dos resultados, devido a utilização de diferentes tipos de equipamentos e procedimentos executivos. O índice NSPT depende da parcela de energia que e efetivamente transmitida ao amostrador durante a queda do martelo. Assim, análises racionais dos resultados de ensaios SPT dependem essencialmente da estimativa da quantidade dessa energia. Por esse motivo, tem sido desenvolvidas diversas pesquisas relacionados com a energia envolvida no ensaio SPT. Este trabalho apresenta interpretações dos resultados do ensaio SPT a partir das quantidades de energia envolvidas durante a queda do martelo. Essas quantidades de energia são determinadas indiretamente, através dos sinais de força normal e aceleração durante a propagação da onda de tensão ao longo da composição de hastes. Para isso, foi desenvolvido um equipamento que inclui uma instrumentação capaz de registrar esses sinais de força normal e aceleração. Uma série de ensaios SPT foi realizada com instrumentação instalada tanto no topo como na base da composição de hastes. Esse arranjo permitiu estimar as quantidades de energia disponíveis no topo e na base da composição de hastes. Assim, foi possível determinar a eficiência dos equipamentos SPT com base na quantidade de energia que efetivamente atinge o amostrador. A partir dos sinais registrados, também foi possível determinar a resistência dinâmica mobilizada no sistema solo-amostrador. Os resultados experimentais foram comparados com resultados obtidos através de métodos teóricos, baseados no Princípio de Hamilton, para a determinação da resistência estática e dinâmica do solo. Essas análises mostraram que equações teóricas podem ser adequadas para a determinação da resistência do solo, desde que seja considerada a quantidade de energia que efetivamente atinge o amostrador. Desta forma, essas equações teóricas podem ser facilmente incorporadas na prática de projeto de fundações, contribuindo para o desenvolvimento de métodos baseados na quantidade de energia necessária para a penetração do amostrador no solo.
The standard penetration test (SPT) is usually used to estimate the soil strength through the NSPT index which represents the dynamic reaction force of the sampler penetration into the soil. This index is directly used in empirical or semi-empirical correlations for the determination of the bearing capacity and foundation settlements. However, this test has been criticized due to the use of these empirical correlations, often based on practical observation without any scientific basis. Criticism is also related to result dispersion due to the use of different types of equipment and execution procedures. The NSPT index depends on the amount of energy that is effectively transmitted to the sampler during the hammer fall. Thus, rational analyses of the SPT test results essentially depend on the estimation of this amount of energy. For this reason, several studies related to the energy involved in SPT test have been developed. This research presents interpretations of the SPT test results from the amounts of energy involved during the hammer fall. These amounts of energy are indirectly assessed by means of normal force and acceleration signals during the stress wave propagation along the string of rods. For this reason, a device that includes instrumentation capable of registering these signals of normal force and acceleration was developed. A series of tests was performed with the instrumentation installed at the top and the bottom of the string of rods. This arrangement allowed estimating the amount of energy available at the top and bottom of the string of rods. Thus, it was possible to estimate the efficiency of the SPT equipment basing on the amount of energy that actually reaches the sampler. From the recorded signals, it was also possible to assess the dynamic reaction force mobilized in the soil-sample system. The results were compared with other results obtained by theoretical methods, based on the Hamilton´s Principle, for evaluating the static and dynamic reaction force of the soil. These analyzes have shown that theoretical equations may be suitable for the determination of the soil strength, provided that the amount of energy that actually reaches the sampler be considered. Thus, these theoretical equations can be easily incorporated into the practice of foundation designs contributing to the development of methods based on the amount of energy required for the penetration of the sampler into the soil.

La aplicación del ensayo de penetración estándar (SPT) conlleva a un nivel de detalle alto, así mismo es costosa, requiere de mano de obra calificada para su ejecución y es demorada; su utilización es recomendada en suelos arenosos, restrictiva en suelos finos y plásticos y no recomendada en los demás tipos de suelos. Mientras que en el Análisis multicanal de ondas superficiales (MASW) su ejecución no se restringe, permite abarcar grandes áreas en tiempos cortos, no obstante, es un método indirecto, razón por la cual su aplicación debe ser verificada mediante perforaciones, obtención de muestras, y altos criterios de análisis. En los últimos años la implementación y mejoramiento de las técnicas de prospección de suelos ha permitido investigar la relación de dichos ensayos con diferentes parámetros del suelo. En la presente se desarrolló una correlación empírica entre los ensayos geotécnicos SPT mediante el número de golpes de resistencia a la penetración (N60) con el valor de las velocidades de ondas de corte (Vs) de los ensayos geofísicos MASW, se realizaron en 30 puntos del distrito de Juliaca para los suelos arenosos. Se contó con un filtro final de 110 pares ordenados de datos (Vs; N60); mediante el análisis de regresión potencial por el método de mínimos cuadrados donde se obtuvo la ecuación: Vs = 81.079*N600.2646; y un coeficiente de determinación R2:0.76. Esta ecuación fue puesta en análisis comparativo con las ecuaciones de las correlaciones de la literatura, en el cual se encontró una tendencia acorde a las mismas, lo que permitirá utilizarla con suelos afines a la región.
The application of the standard penetration test (SPT) entails a high level of detail, it is also expensive, requires qualified labor for its execution and is time-consuming; Its use is recommended in sandy soils, restrictive in fine and plastic soils and not recommended in other types of soils. While in the multichannel analysis of surface waves (MASW) its execution is not restricted, it allows covering large areas in short times, however, it is an indirect method, which is why its application must be verified by drilling, obtaining samples, and high analytical criteria. In recent years, the implementation and improvement of soil prospecting techniques has made it possible to carry out studies of the relationship of said tests with different soil parameters. In the present, an empirical correlation was developed between the SPT geotechnical tests by means of the number of penetration resistance blows (N60) with the value of the shear wave velocities (Vs) of the MASW geophysical tests, these tests were carried out in 30 points in the city of Juliaca for sandy soils. There was a final filter of 110 ordered pairs of data (Vs; N60); Through the analysis of potential regression by the method of least squares, the equation was obtained: Vs = 81.079*N600.2646; of the determination coefficient R2:0.76, it was determined that the correlation is acceptable. This equation was put into analysis and comparison with the equations of the correlations of the literature, in which it was found in a trend according to them.
Tesis

Um der Mikrowellenerwärmung als Single- oder Hybridanwendung im industriellen Ofenbau ein weiteres Anwendungsspektrum zu eröffnen, wird dem Anlagenkonstrukteur im Rahmen dieser Arbeit eine Aufstellung von Entwicklungsprioritäten und Konstruktionshinweisen übergeben. Für eine energieeffiziente Ausnutzung des technologischen Potenzials wird hierbei nicht - wie in der Vergangenheit - die Feldstärkenverteilung, sondern das Reflexionsverhalten in den Vordergrund gestellt. Durch die Entwicklung und den Aufbau eines Messplatzes zur Bestimmung temperaturspezifischer Dielektrizitätsparameter konnte die Realitätsnähe der Simulation für die anwendungsorientierte Applikatorentwicklung gesteigert werden. Die Anwendbarkeit der Auslegungskriterien wird im Rahmen der Projektierung und des Aufbaus einer Beispielanlage zur Entbinderung von keramischen Grünkörpern und bei der Entwicklung zahlreicher mikrowellenspezifischer Anlagenbauteile verifiziert.

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

This thesis studies the effects of import competition from China and Eastern Europe on the health and fertility decisions of German individuals working in manufacturing. Individuals are matched with separate measures of exposure to competition from China and Eastern Europe, respectively. To isolate exogenous supply shocks from the origin, instrumental variables for competition from each of China and Eastern Europe are constructed. Results in Chapter 1 suggest that higher import competition worsens individual health via job displacement, wage decline, shortened employment duration, increased reliance on welfare and less future orientation, with Chinese import competition affecting individuals twice as much. Health declines as individuals increase their visits to the doctor, exercise less frequently and have a higher probability of developing chronic illness. Also, there is some evidence that individuals do not tend to become disabled but may be slowly pushed into chronic illness. Findings in Chapter 2 show that import competition negatively affects the individual’s probability of having children via reduced earnings, lower satisfaction with personal income and shortened employment duration. The chapter then investigates effects of import exposure by gender. Results show that male and female fertility choices differ upon rising import competition. Higher import exposure lowers female earnings and job autonomy, which in turn generates a lower opportunity cost of work, to the point where having children would become a more rewarding alternative for female workers. By contrast, increased import exposure negatively affects male workers’ fertility through reduced earnings and employment duration.

Demonstration of utter stability over the design envelope and identi fication of safe ight envelope is a prerequisite for operational clearance of any new aircraft design. An important step involved in flight flutter testing is the proper use of flutter prediction techniques to accurately predict flutter. This study focuses on flutter estimation techniques given flight test data. We review various flutter prediction techniques by simulating them on a three-degree-of-freedom aeroelastic system. We show that flutter margin and an auto-regressive model based approach are robust and reasonably accurate in predicting flutter onset. Using these two techniques, flutter margins are computed at flight test points obtained from flight test data of flexible aircraft. Throughout, we have predicted flutter dynamic pressure at constant Mach number. One of the contributions of this work is in predicting flutter dynamic pressure at transonic Mach numbers. A critical issue in flutter prediction is the lack of information on the flutter instability mode and thereby the number of modes to include in a model. In a novel application of tools from statistical signal processing, we determine the optimum model order to construct an auto-regressive model using the flight flutter test time response data. From this auto-regressive model the frequency and damping values are estimated which in turn is used to estimate aeroelastic stability parameters. High resolution property of auto-regressive technique even with short data records is demonstrated in this thesis. This will provide a quick evaluation of spectral estimate and stability parameter using the same auto-regressive model, facilitating a quick envelope expansion. Aeroelastic stability during wake penetration is an essential part of operational clearance of new aircraft. In this thesis, and perhaps one of the first such study, wake flight test response data is used to assess the stability of the aircraft during wake penetration by modeling the wake response data in an auto-regressive framework. We have compared analytical predictions of incremental load factor based on simulations with flight tests on wake interactions. Estimates for safe wake encounter distance that do not exceed structural limit load factors have been determined.

Part I.

We have developed a technique for measuring the depth time history of rigid body penetration into brittle materials (hard rocks and concretes) under a deceleration of ~ 10⁵ g. The technique includes bar-coded projectile, sabot-projectile separation, detection and recording systems. Because the technique can give very dense data on penetration depth time history, penetration velocity can be deduced. Error analysis shows that the technique has a small intrinsic error of ~ 3-4 % in time during penetration, and 0.3 to 0.7 mm in penetration depth. A series of 4140 steel projectile penetration into G-mixture mortar targets have been conducted using the Caltech 40 mm gas/ powder gun in the velocity range of 100 to 500 m/s.

We report, for the first time, the whole depth-time history of rigid body penetration into brittle materials (the G-mixture mortar) under 10⁵ g deceleration. Based on the experimental results, including penetration depth time history, damage of recovered target and projectile materials and theoretical analysis, we find:

1. Target materials are damaged via compacting in the region in front of a projectile and via brittle radial and lateral crack propagation in the region surrounding the penetration path. The results suggest that expected cracks in front of penetrators may be stopped by a comminuted region that is induced by wave propagation. Aggregate erosion on the projectile lateral surface is < 20% of the final penetration depth. This result suggests that the effect of lateral friction on the penetration process can be ignored.

2. Final penetration depth, P_max, is linearly scaled with initial projectile energy per unit cross-section area, e_s , when targets are intact after impact. Based on the experimental data on the mortar targets, the relation is P_max(mm) 1.15e_s (J/mm² ) + 16.39.

3. Estimation of the energy needed to create an unit penetration volume suggests that the average pressure acting on the target material during penetration is ~ 10 to 20 times higher than the unconfined strength of target materials under quasi-static loading, and 3 to 4 times higher than the possible highest pressure due to friction and material strength and its rate dependence. In addition, the experimental data show that the interaction between cracks and the target free surface significantly affects the penetration process.

4. Based on the fact that the penetration duration, t_max, increases slowly with e_s and does not depend on projectile radius approximately, the dependence of t_max on projectile length is suggested to be described by t_max(μs) = 2.08e_s (J/mm² + 349.0 x m/(πR²), in which m is the projectile mass in grams and R is the projectile radius in mm. The prediction from this relation is in reasonable agreement with the experimental data for different projectile lengths.

5. Deduced penetration velocity time histories suggest that whole penetration history is divided into three stages: (1) An initial stage in which the projectile velocity change is small due to very small contact area between the projectile and target materials; (2) A steady penetration stage in which projectile velocity continues to decrease smoothly; (3) A penetration stop stage in which projectile deceleration jumps up when velocities are close to a critical value of ~ 35 m/s.

6. Deduced averaged deceleration, a, in the steady penetration stage for projectiles with same dimensions is found to be a(g) = 192.4v + 1.89 x 10⁴, where v is initial projectile velocity in m/s. The average pressure acting on target materials during penetration is estimated to be very comparable to shock wave pressure.

7. A similarity of penetration process is found to be described by a relation between normalized penetration depth, P/P_max, and normalized penetration time, t/t_max, as P/P_max = f(t/t_max, where f is a function of t/t_max. After f(t/t_max is determined using experimental data for projectiles with 150 mm length, the penetration depth time history for projectiles with 100 mm length predicted by this relation is in good agreement with experimental data. This similarity also predicts that average deceleration increases with decreasing projectile length, that is verified by the experimental data.

8. Based on the penetration process analysis and the present data, a first principle model for rigid body penetration is suggested. The model incorporates the models for contact area between projectile and target materials, friction coefficient, penetration stop criterion, and normal stress on the projectile surface. The most important assumptions used in the model are: (1) The penetration process can be treated as a series of impact events, therefore, pressure normal to projectile surface is estimated using the Hugoniot relation of target material; (2) The necessary condition for penetration is that the pressure acting on target materials is not lower than the Hugoniot elastic limit; (3) The friction force on projectile lateral surface can be ignored due to cavitation during penetration. All the parameters involved in the model are determined based on independent experimental data. The penetration depth time histories predicted from the model are in good agreement with the experimental data.

9. Based on planar impact and previous quasi-static experimental data, the strain rate dependence of the mortar compressive strength is described by σ_f/σ⁰_f = exp(0.0905(log(έ/έ_0) ^1.14, in the strain rate range of 10^-7/s to 10³/s (σ⁰_f and έ are reference compressive strength and strain rate, respectively). The non-dispersive Hugoniot elastic wave in the G-mixture has an amplitude of ~ 0.14 GPa and a velocity of ~ 4.3 km/s.

Part II.

Stress wave profiles in vitreous GeO₂ were measured using piezoresistance gauges in the pressure range of 5 to 18 GPa under planar plate and spherical projectile impact. Experimental data show that the response of vitreous GeO₂ to planar shock loading can be divided into three stages: (1) A ramp elastic precursor has peak amplitude of 4 GPa and peak particle velocity of 333 m/s. Wave velocity decreases from initial longitudinal elastic wave velocity of 3.5 km/s to 2.9 km/s at 4 GPa; (2) A ramp wave with amplitude of 2.11 GPa follows the precursor when peak loading pressure is 8.4 GPa. Wave velocity drops to the value below bulk wave velocity in this stage; (3) A shock wave achieving final shock state forms when peak pressure is > 6 GPa. The Hugoniot relation is D = 0.917 + 1.711u (km/s) using present data and the data of Jackson and Ahrens [1979] when shock wave pressure is between 6 and 40 GPa for ρ₀ = 3.655 gj cm³ . Based on the present data, the phase change from 4-fold to 6-fold coordination of Ge⁺⁴ with O^-2 in vitreous GeO₂ occurs in the pressure range of 4 to 15 ± 1 GPa under planar shock loading. Comparison of the shock loading data for fused SiO₂ to that on vitreous GeO₂ demonstrates that transformation to the rutile structure in both media are similar. The Hugoniots of vitreous GeO₂ and fused SiO₂ are found to coincide approximately if pressure in fused SiO₂ is scaled by the ratio of fused SiO₂to vitreous GeO₂ density. This result, as well as the same structure, provides the basis for considering vitreous Ge0₂ as an analogous material to fused SiO₂ under shock loading. Experimental results from the spherical projectile impact demonstrate: (1) The supported elastic shock in fused SiO₂ decays less rapidly than a linear elastic wave when elastic wave stress amplitude is higher than 4 GPa. The supported elastic shock in vitreous GeO₂ decays faster than a linear elastic wave; (2) In vitreous GeO₂ , unsupported shock waves decays with peak pressure in the phase transition range (4-15 GPa) with propagation distance, x, as α 1/x^-3.35 , close to the prediction of Chen et al. [1998]. Based on a simple analysis on spherical wave propagation, we find that the different decay rates of a spherical elastic wave in fused SiO₂ and vitreous GeO₂ is predictable on the base of the compressibility variation with stress under one-dimensional strain condition in the two materials.

碩士
中華大學
土木工程學系
104
As time passes, an increasing number of buildings in Taiwan are characterized by old structures. In fact, there are more than three million old buildings in Taiwan. In addition, the country suffers from frequent earthquakes as it is situated at the boundary between the Philippine Sea Plate and Eurasian Sea Plate. Therefore, leakage is an inevitable problem in building structures. According to statistics, the majority reason for pipeline ruptures is loaded vibration, followed by corrosion of aging pipelines. Ninety percent of ruptured pipelines are PVC pipes. Using non-invasive ground-penetrating radars, this study detected internal leakage inside building structures. These radars excel in precision and convenience, and allow pipeline maintenance or replacement to be carried out with minimal excavation. This study employed the ground-penetrating radar’s 1 GHz probe to detect test leakages in pipes made of different materials buried in standard sand and PVC pipes in doubly reinforced concrete. The algorithm was applied to investigate real-life cases of leakages inside building structures. By using the basic theory of electromagnetic waves of ground-penetrating radars and by extracting and analyzing their digitized signals, observed variations in equivalent permittivity at the dry-wet interface were used to determine the area and source of leakage. These were rendered in 2D and 3D images to investigate the feasibility of such a method to detect leakages under different site conditions.

碩士
中華大學
土木工程學系碩士班
102
This study uses the 1GHz ground penetrating radar(GPR) to inspect the leakage of the pipeline in sand and mortar with the rebar. The research applies the basic theory of electricmagnetic wave and digital signal acquisition technique to explore the equivalent relative dielectric constant changes of the Ottawa sand and cement mortar containing reinforcement under the situation of pipeline leakages. The experiment analyzes the related leakage signal characteristics and changes according to the electromagnetic wave reflections, time difference and equivalent relative dielectric constant. It then uses the difference of the equivalent relative dielectric constant of GPR in dry Ottawa sand and cement mortar material to determine the range of leakage distribution in the material.

碩士
大同工學院
電機工程學系
84
In this thesis, the penetrative imaging of 2D objects embedded in layeredmedia under multifrequency plane waves illumination is being investigated. Weknow that the Green's function of a line source embedded in layered media can beformulated by using plane wave expansion. Here, the multifrequency plane wave backward propagation is being proposed based on the Green's function mensioned above. The FFT is applied on computation process of the Fourier transform,futher, it's employed in the computation of Green's function, penetrating fieldand image reconstruction. Besides, the performance of the penetrative imageunder the interference of uniform distributed random noise and multiple reflected field is also discussed. An experimental measurement is implemented accordingly. Finally, further developments for practical applications in nondestructive detection are forseen by the experimental results together with their simulations.

博士
大同工學院
電機工程學系
84
ABSTRACT -- A penetrative imaging of 2D or 3D objects embedded in planar multilayer media is investigated under a multifrequency plane wave illumination. Both theoretical and experimental results are presented in this dissertation. A 3D multilayer Green''s function is derived in terms of plane wave superposition by using a 2D Fourier transform (FT). Based on this multilayer Green''s function, a penetrative imaging of 3D embedded objects can be expressed as multifrequency plane wave backward propagation. The 3D penetrative image can be reconstructed by a 2D inverse Fourier transform (IFT) on stacks of parallel planes, and the 2D penetrative imaging can be reconstructed by a 1D IFT on stacks of parallel lines. If the embedded objects are in the last layer, the penetrative imaging can be implemented with a 2D or 3D IFT for 2D and 3D cases, respectively. The penetrative imaging of both embedded conductors and dielectrics are demonstrated. Moreover, the performance of penetrative imaging with the random noise and the interference is discussed. An automatic measurement system is implemented for the penetrative imaging. Experimental results under various multilayer media and embedded objects are presented to demonstrate the performance of the imaging algorithm. The numerical and experimental results reveal that the proposed penetrative imaging can be applied to both the 2D and 3D embedded objects whether they are conductors or dielectrics. Besides, the method is insensitive to the high random noise and the high interference up to 100% of the peak of scattered field. Further developments into practical applications in non-destructive testing and underground detection are foreseen.

碩士
中華大學
土木工程學系碩士班
99
This study applied Ground Penetrating Radar of 1GHz against standard sand (Ottawa sand) with moisture that contained the bearing layer in different depths, and that contained pipeline leak’s experiment for the exploration. Based on the basic theory of electromagnetic wave (EM) and wave speed in connection with digital signal mining technique, the experiment processed non-destructive scanning with Ground Penetrating Radar of 1GHz to explore standard sand (Ottawa sand) with moisture that contained the bearing layer in different depths, and that contained pipeline leak. According to analysis of the relevant leak signal’s feature and change in electromagnetic wave (EM) reflected signal, waveform, time difference, frequency spectrum properties, wave speed and relative permittivity constants, etc. of electromagnetic wave (EM), it is concluded that the corresponding relationship between electromagnetic wave (EM) of Ground Penetrating Radar and leak in the standard sand, suggesting to detach dry-wet sand interface by using time difference in electromagnetic wave (EM) speed of Ground Penetrating Radar, and establishing a judgment method for leak range and leak source.

The objective of this Thesis consists in presenting in a concise and effective way the results achieved during researches and studies on satellite systems interference sources, advanced antenna arrays for satellite systems to mitigate the increasing anti-interference needs and on an innovative way to generates inhomogeneous wave in lossless media for contributing to the design of a novel type of antenna for deep penetrating lossy media.

碩士
國立中央大學
土木工程研究所
95
The times in fast variety, people work to always hope fast and then valid, had better can also save expense in addition, deeply the land mine attain a method opposite open to dig in the clear ditch, ground-penetrating radar to explore, obvious since quick and again save expense, in the civil engineering in carry on in the city underground open to dig to usually concern the problem of underground pipeline, if can not control a pipeline to lay position correctly in times before, will cause to construct in the engineering in the days to come last inconvenience and risk, and if adopt traditional clear ditch to open the way for digging to bring about greater influence on transportation easily, and not agree with to match an economic benefit.High Pin ground-penetrating radar deeply the land mine attain to have a high resolution, the signal repetition is high, and tests a fast advantage. So adoption fast, the economy is high of nondestructive evaluation way concert the related data take in to judge, is already the trend of future development, in order to being easily subjected to direct wave rather great influence when the time area measurement''s analysis reflect a signal, often cause shallow layer signal easily with direct wave mutual interference and hard discretion.This research mainly experiments through the school and now experiment, understand the signal characteristic of radar wave, again through behind handle of way, improve deeply land mine to attain of signal. According to the analytical result manifestation of the experimenting of school and scene, behind handle program ability will direct wave signal with keep on the existent background noise clean, having further control to the analysis that need to be measured a thing signal and the judgment.

yes
This work presents the development of a new approach of modelling the source excitation and the penetration of structures by continuous propagating electromagnetic (EM) plane waves. The technique incorporates the solution of time-dependent Maxwell¿s equations and the initial value problem as the structures are illuminated by the plane waves. The propagation of waves from source excitation is simulated by solving a finite-difference Maxwell's equation in the time domain. Subgridding method is used to condense the lattice at the point of interest locally for observing field distribution in high resolution. The computational burden due to huge number of time steps has been eased by employing quasi-static approach. An example of induced EM fields near an underground pipeline runs parallel to a 132 kV overhead power transmission line (OHTL) has been presented which paves the way in the development of new approach of EM fields interaction modelling.
MSCRC

Site characterization is an important task in Geotechnical Engineering. In situ tests based on standard penetration test (SPT), cone penetration test (CPT) and shear wave velocity survey are popular among geotechnical engineers. Site characterization using any of these properties based on finite number of in-situ test data is an imperative task in probabilistic site characterization. These methods have been used to design future soil sampling programs for the site and to specify the soil stratification. It is never possible to know the geotechnical properties at every location beneath an actual site because, in order to do so, one would need to sample and/or test the entire subsurface profile. Therefore, the main objective of site characterization models is to predict the subsurface soil properties with minimum in-situ test data. The prediction of soil property is a difficult task due to the uncertainities. Spatial variability, measurement ‘noise’, measurement and model bias, and statistical error due to limited measurements are the sources of uncertainities. Liquefaction in soil is one of the other major problems in geotechnical earthquake engineering. It is defined as the transformation of a granular material from a solid to a liquefied state as a consequence of increased pore-water pressure and reduced effective stress. The generation of excess pore pressure under undrained loading conditions is a hallmark of all liquefaction phenomena. This phenomena was brought to the attention of engineers more so after Niigata(1964) and Alaska(1964) earthquakes. Liquefaction will cause building settlement or tipping, sand boils, ground cracks, landslides, dam instability, highway embankment failures, or other hazards. Such damages are generally of great concern to public safety and are of economic significance. Site-spefific evaluation of liquefaction susceptibility of sandy and silty soils is a first step in liquefaction hazard assessment. Many methods (intelligent models and simple methods as suggested by Seed and Idriss, 1971) have been suggested to evaluate liquefaction susceptibility based on the large data from the sites where soil has been liquefied / not liquefied. The rapid advance in information processing systems in recent decades directed engineering research towards the development of intelligent models that can model natural phenomena automatically. In intelligent model, a process of training is used to build up a model of the particular system, from which it is hoped to deduce responses of the system for situations that have yet to be observed. Intelligent models learn the input output relationship from the data itself. The quantity and quality of the data govern the performance of intelligent model. The objective of this study is to develop intelligent models [geostatistic, artificial neural network(ANN) and support vector machine(SVM)] to estimate corrected standard penetration test (SPT) value, Nc, in the three dimensional (3D) subsurface of Bangalore. The database consists of 766 boreholes spread over a 220 sq km area, with several SPT N values (uncorrected blow counts) in each of them. There are total 3015 N values in the 3D subsurface of Bangalore. To get the corrected blow counts, Nc, various corrections such as for overburden stress, size of borehole, type of sampler, hammer energy and length of connecting rod have been applied on the raw N values. Using a large database of Nc values in the 3D subsurface of Bangalore, three geostatistical models (simple kriging, ordinary kriging and disjunctive kriging) have been developed. Simple and ordinary kriging produces linear estimator whereas, disjunctive kriging produces nonlinear estimator. The knowledge of the semivariogram of the Nc data is used in the kriging theory to estimate the values at points in the subsurface of Bangalore where field measurements are not available. The capability of disjunctive kriging to be a nonlinear estimator and an estimator of the conditional probability is explored. A cross validation (Q1 and Q2) analysis is also done for the developed simple, ordinary and disjunctive kriging model. The result indicates that the performance of the disjunctive kriging model is better than simple as well as ordinary kriging model. This study also describes two ANN modelling techniques applied to predict Nc data at any point in the 3D subsurface of Bangalore. The first technique uses four layered feed-forward backpropagation (BP) model to approximate the function, Nc=f(x, y, z) where x, y, z are the coordinates of the 3D subsurface of Bangalore. The second technique uses generalized regression neural network (GRNN) that is trained with suitable spread(s) to approximate the function, Nc=f(x, y, z). In this BP model, the transfer function used in first and second hidden layer is tansig and logsig respectively. The logsig transfer function is used in the output layer. The maximum epoch has been set to 30000. A Levenberg-Marquardt algorithm has been used for BP model. The performance of the models obtained using both techniques is assessed in terms of prediction accuracy. BP ANN model outperforms GRNN model and all kriging models. SVM model, which is firmly based on the theory of statistical learning theory, uses regression technique by introducing -insensitive loss function has been also adopted to predict Nc data at any point in 3D subsurface of Bangalore. The SVM implements the structural risk minimization principle (SRMP), which has been shown to be superior to the more traditional empirical risk minimization principle (ERMP) employed by many of the other modelling techniques. The present study also highlights the capability of SVM over the developed geostatistic models (simple kriging, ordinary kriging and disjunctive kriging) and ANN models. Further in this thesis, Liquefaction susceptibility is evaluated from SPT, CPT and Vs data using BP-ANN and SVM. Intelligent models (based on ANN and SVM) are developed for prediction of liquefaction susceptibility using SPT data from the 1999 Chi-Chi earthquake, Taiwan. Two models (MODEL I and MODEL II) are developed. The SPT data from the work of Hwang and Yang (2001) has been used for this purpose. In MODEL I, cyclic stress ratio (CSR) and corrected SPT values (N1)60 have been used for prediction of liquefaction susceptibility. In MODEL II, only peak ground acceleration (PGA) and (N1)60 have been used for prediction of liquefaction susceptibility. Further, the generalization capability of the MODEL II has been examined using different case histories available globally (global SPT data) from the work of Goh (1994). This study also examines the capabilities of ANN and SVM to predict the liquefaction susceptibility of soils from CPT data obtained from the 1999 Chi-Chi earthquake, Taiwan. For determination of liquefaction susceptibility, both ANN and SVM use the classification technique. The CPT data has been taken from the work of Ku et al.(2004). In MODEL I, cone tip resistance (qc) and CSR values have been used for prediction of liquefaction susceptibility (using both ANN and SVM). In MODEL II, only PGA and qc have been used for prediction of liquefaction susceptibility. Further, developed MODEL II has been also applied to different case histories available globally (global CPT data) from the work of Goh (1996). Intelligent models (ANN and SVM) have been also adopted for liquefaction susceptibility prediction based on shear wave velocity (Vs). The Vs data has been collected from the work of Andrus and Stokoe (1997). The same procedures (as in SPT and CPT) have been applied for Vs also. SVM outperforms ANN model for all three models based on SPT, CPT and Vs data. CPT method gives better result than SPT and Vs for both ANN and SVM models. For CPT and SPT, two input parameters {PGA and qc or (N1)60} are sufficient input parameters to determine the liquefaction susceptibility using SVM model. In this study, an attempt has also been made to evaluate geotechnical site characterization by carrying out in situ tests using different in situ techniques such as CPT, SPT and multi channel analysis of surface wave (MASW) techniques. For this purpose a typical site was selected wherein a man made homogeneous embankment and as well natural ground has been met. For this typical site, in situ tests (SPT, CPT and MASW) have been carried out in different ground conditions and the obtained test results are compared. Three CPT continuous test profiles, fifty-four SPT tests and nine MASW test profiles with depth have been carried out for the selected site covering both homogeneous embankment and natural ground. Relationships have been developed between Vs, (N1)60 and qc values for this specific site. From the limited test results, it was found that there is a good correlation between qc and Vs. Liquefaction susceptibility is evaluated using the in situ test data from (N1)60, qc and Vs using ANN and SVM models. It has been shown to compare well with “Idriss and Boulanger, 2004” approach based on SPT test data. SVM model has been also adopted to determine over consolidation ratio (OCR) based on piezocone data. Sensitivity analysis has been performed to investigate the relative importance of each of the input parameters. SVM model outperforms all the available methods for OCR prediction.

Site characterization is an important task in Geotechnical Engineering. In situ tests based on standard penetration test (SPT), cone penetration test (CPT) and shear wave velocity survey are popular among geotechnical engineers. Site characterization using any of these properties based on finite number of in-situ test data is an imperative task in probabilistic site characterization. These methods have been used to design future soil sampling programs for the site and to specify the soil stratification. It is never possible to know the geotechnical properties at every location beneath an actual site because, in order to do so, one would need to sample and/or test the entire subsurface profile. Therefore, the main objective of site characterization models is to predict the subsurface soil properties with minimum in-situ test data. The prediction of soil property is a difficult task due to the uncertainities. Spatial variability, measurement ‘noise’, measurement and model bias, and statistical error due to limited measurements are the sources of uncertainities. Liquefaction in soil is one of the other major problems in geotechnical earthquake engineering. It is defined as the transformation of a granular material from a solid to a liquefied state as a consequence of increased pore-water pressure and reduced effective stress. The generation of excess pore pressure under undrained loading conditions is a hallmark of all liquefaction phenomena. This phenomena was brought to the attention of engineers more so after Niigata(1964) and Alaska(1964) earthquakes. Liquefaction will cause building settlement or tipping, sand boils, ground cracks, landslides, dam instability, highway embankment failures, or other hazards. Such damages are generally of great concern to public safety and are of economic significance. Site-spefific evaluation of liquefaction susceptibility of sandy and silty soils is a first step in liquefaction hazard assessment. Many methods (intelligent models and simple methods as suggested by Seed and Idriss, 1971) have been suggested to evaluate liquefaction susceptibility based on the large data from the sites where soil has been liquefied / not liquefied. The rapid advance in information processing systems in recent decades directed engineering research towards the development of intelligent models that can model natural phenomena automatically. In intelligent model, a process of training is used to build up a model of the particular system, from which it is hoped to deduce responses of the system for situations that have yet to be observed. Intelligent models learn the input output relationship from the data itself. The quantity and quality of the data govern the performance of intelligent model. The objective of this study is to develop intelligent models [geostatistic, artificial neural network(ANN) and support vector machine(SVM)] to estimate corrected standard penetration test (SPT) value, Nc, in the three dimensional (3D) subsurface of Bangalore. The database consists of 766 boreholes spread over a 220 sq km area, with several SPT N values (uncorrected blow counts) in each of them. There are total 3015 N values in the 3D subsurface of Bangalore. To get the corrected blow counts, Nc, various corrections such as for overburden stress, size of borehole, type of sampler, hammer energy and length of connecting rod have been applied on the raw N values. Using a large database of Nc values in the 3D subsurface of Bangalore, three geostatistical models (simple kriging, ordinary kriging and disjunctive kriging) have been developed. Simple and ordinary kriging produces linear estimator whereas, disjunctive kriging produces nonlinear estimator. The knowledge of the semivariogram of the Nc data is used in the kriging theory to estimate the values at points in the subsurface of Bangalore where field measurements are not available. The capability of disjunctive kriging to be a nonlinear estimator and an estimator of the conditional probability is explored. A cross validation (Q1 and Q2) analysis is also done for the developed simple, ordinary and disjunctive kriging model. The result indicates that the performance of the disjunctive kriging model is better than simple as well as ordinary kriging model. This study also describes two ANN modelling techniques applied to predict Nc data at any point in the 3D subsurface of Bangalore. The first technique uses four layered feed-forward backpropagation (BP) model to approximate the function, Nc=f(x, y, z) where x, y, z are the coordinates of the 3D subsurface of Bangalore. The second technique uses generalized regression neural network (GRNN) that is trained with suitable spread(s) to approximate the function, Nc=f(x, y, z). In this BP model, the transfer function used in first and second hidden layer is tansig and logsig respectively. The logsig transfer function is used in the output layer. The maximum epoch has been set to 30000. A Levenberg-Marquardt algorithm has been used for BP model. The performance of the models obtained using both techniques is assessed in terms of prediction accuracy. BP ANN model outperforms GRNN model and all kriging models. SVM model, which is firmly based on the theory of statistical learning theory, uses regression technique by introducing -insensitive loss function has been also adopted to predict Nc data at any point in 3D subsurface of Bangalore. The SVM implements the structural risk minimization principle (SRMP), which has been shown to be superior to the more traditional empirical risk minimization principle (ERMP) employed by many of the other modelling techniques. The present study also highlights the capability of SVM over the developed geostatistic models (simple kriging, ordinary kriging and disjunctive kriging) and ANN models. Further in this thesis, Liquefaction susceptibility is evaluated from SPT, CPT and Vs data using BP-ANN and SVM. Intelligent models (based on ANN and SVM) are developed for prediction of liquefaction susceptibility using SPT data from the 1999 Chi-Chi earthquake, Taiwan. Two models (MODEL I and MODEL II) are developed. The SPT data from the work of Hwang and Yang (2001) has been used for this purpose. In MODEL I, cyclic stress ratio (CSR) and corrected SPT values (N1)60 have been used for prediction of liquefaction susceptibility. In MODEL II, only peak ground acceleration (PGA) and (N1)60 have been used for prediction of liquefaction susceptibility. Further, the generalization capability of the MODEL II has been examined using different case histories available globally (global SPT data) from the work of Goh (1994). This study also examines the capabilities of ANN and SVM to predict the liquefaction susceptibility of soils from CPT data obtained from the 1999 Chi-Chi earthquake, Taiwan. For determination of liquefaction susceptibility, both ANN and SVM use the classification technique. The CPT data has been taken from the work of Ku et al.(2004). In MODEL I, cone tip resistance (qc) and CSR values have been used for prediction of liquefaction susceptibility (using both ANN and SVM). In MODEL II, only PGA and qc have been used for prediction of liquefaction susceptibility. Further, developed MODEL II has been also applied to different case histories available globally (global CPT data) from the work of Goh (1996). Intelligent models (ANN and SVM) have been also adopted for liquefaction susceptibility prediction based on shear wave velocity (Vs). The Vs data has been collected from the work of Andrus and Stokoe (1997). The same procedures (as in SPT and CPT) have been applied for Vs also. SVM outperforms ANN model for all three models based on SPT, CPT and Vs data. CPT method gives better result than SPT and Vs for both ANN and SVM models. For CPT and SPT, two input parameters {PGA and qc or (N1)60} are sufficient input parameters to determine the liquefaction susceptibility using SVM model. In this study, an attempt has also been made to evaluate geotechnical site characterization by carrying out in situ tests using different in situ techniques such as CPT, SPT and multi channel analysis of surface wave (MASW) techniques. For this purpose a typical site was selected wherein a man made homogeneous embankment and as well natural ground has been met. For this typical site, in situ tests (SPT, CPT and MASW) have been carried out in different ground conditions and the obtained test results are compared. Three CPT continuous test profiles, fifty-four SPT tests and nine MASW test profiles with depth have been carried out for the selected site covering both homogeneous embankment and natural ground. Relationships have been developed between Vs, (N1)60 and qc values for this specific site. From the limited test results, it was found that there is a good correlation between qc and Vs. Liquefaction susceptibility is evaluated using the in situ test data from (N1)60, qc and Vs using ANN and SVM models. It has been shown to compare well with “Idriss and Boulanger, 2004” approach based on SPT test data. SVM model has been also adopted to determine over consolidation ratio (OCR) based on piezocone data. Sensitivity analysis has been performed to investigate the relative importance of each of the input parameters. SVM model outperforms all the available methods for OCR prediction.

Um der Mikrowellenerwärmung als Single- oder Hybridanwendung im industriellen Ofenbau ein weiteres Anwendungsspektrum zu eröffnen, wird dem Anlagenkonstrukteur im Rahmen dieser Arbeit eine Aufstellung von Entwicklungsprioritäten und Konstruktionshinweisen übergeben. Für eine energieeffiziente Ausnutzung des technologischen Potenzials wird hierbei nicht - wie in der Vergangenheit - die Feldstärkenverteilung, sondern das Reflexionsverhalten in den Vordergrund gestellt. Durch die Entwicklung und den Aufbau eines Messplatzes zur Bestimmung temperaturspezifischer Dielektrizitätsparameter konnte die Realitätsnähe der Simulation für die anwendungsorientierte Applikatorentwicklung gesteigert werden. Die Anwendbarkeit der Auslegungskriterien wird im Rahmen der Projektierung und des Aufbaus einer Beispielanlage zur Entbinderung von keramischen Grünkörpern und bei der Entwicklung zahlreicher mikrowellenspezifischer Anlagenbauteile verifiziert.

Mankind is facing the problem due to earthquake hazard since prehistoric times. Many of the developed and developing countries are under constant threats from earthquakes hazards. Theories of plate tectonics and engineering seismology have helped to understand earthquakes and also to predicate earthquake hazards on a regional scale. However, the regional scale hazard mapping in terms of seismic zonation has been not fully implemented in many of the developing countries like India. Agglomerations of large population in the Indian cities and poor constructions have raised the risk due to various possible seismic hazards. First and foremost step towards hazard reduction is estimation of the seismic hazards in regional scale. Objective of this study is to estimate the seismic hazard parameters for Lucknow, a part of Indo-Gangetic Basin (IGB) and develop regional scale microzonation map. Lucknow is a highly populated city which is located close to the active seismic belt of Himalaya. This belt came into existence during the Cenozoic era (40-50 million years ago) and is a constant source of seismic threats. Many of the devastating earthquakes which have happened since prehistoric times such as 1255 Nepal, 1555 Srinagar, 1737 Kolkata, 1803 Nepal, 1833 Kathmandu, 1897 Shillong, 1905 Kangra, 1934 Bihar-Nepal, 1950 Assam and 2005 Kashmir. Historic evidences show that many of these earthquakes had caused fatalities even up to 0.1 million. At present, in the light of building up strains and non-occurrence of a great event in between 1905 Kangra earthquake and 1934 Bihar-Nepal earthquake regions the stretch has been highlighted as central seismic gap. This location may have high potential of great earthquakes in the near future. Geodetic studies in these locations indicate a possible slip of 9.5 m which may cause an event of magnitude 8.7 on Richter scale in the central seismic gap. Lucknow, the capital of Uttar Pradesh has a population of 2.8 million as per Census 2011. It lies in ZONE III as per IS1893: 2002 and can be called as moderate seismic region. However, the city falls within 350 km radial distance from Main Boundary Thrust (MBT) and active regional seismic source of the Lucknow-Faizabad fault. Considering the ongoing seismicity of Himalayan region and the Lucknow-Faizabad fault, this city is under high seismic threat. Hence a comprehensive study of understanding the earthquake hazards on a regional scale for the Lucknow is needed. In this work the seismic microzonation of Lucknow has been attempted. The whole thesis is divided into 11 chapters. A detailed discussion on the importance of this study, seismicity of Lucknow, and methodology adopted for detailed seismic hazard assessment and microzonation are presented in first three chapters. Development of region specific Ground Motion Prediction Equation (GMPE) and seismic hazard estimation at bedrock level using highly ranked GMPEs are presented in Chapters 4 and 5 respectively. Subsurface lithology, measurement of dynamic soil properties and correlations are essential to assess region specific site effects and liquefaction potential. Discussion on the experimental studies, subsurface profiling using geotechnical and geophysical tests results and correlation between shear wave velocity (SWV) and standard penetration test (SPT) N values are presented in Chapter 6. Detailed shear wave velocity profiling with seismic site classification and ground response parameters considering multiple ground motion data are discussed in Chapters 7 and 8. Chapters 9 and 10 present the assessment of liquefaction potential and determination of hazard index with microzonation maps respectively. Conclusions derived from each chapter are presented in Chapter 11. A brief summary of the work is presented below: Attenuation relations or GMPEs are important component of any seismic hazard analysis which controls accurate prediction of the hazard values. Even though the Himalayas have experienced great earthquakes since ancient times, suitable GMPEs which are applicable for a wide range of distance and magnitude are limited. Most of the available regional GMPEs were developed considering limited recorded data and/or pure synthetic ground motion data. This chapter presents development of a regional GMPE considering both the recorded as well as synthetic ground motions. In total 14 earthquakes consisting of 10 events with recorded data and 4 historic events with Isoseismal maps are used for the same. Synthetic ground motions based on finite fault model have been generated at unavailable locations for recorded events and complete range distances for historic earthquakes. Model parameters for synthetic ground motion were arrived by detailed parametric study and from literatures. A concept of Apparent Stations (AS) has been used to generate synthetic ground motion in a wide range of distance as well as direction around the epicenter. Synthetic ground motion data is validated by comparing with available recorded data and peak ground acceleration (PGA) from Isoseismal maps. A new GMPE has been developed based on two step stratified regression procedure considering the combined dataset of recorded and synthetic ground motions. The new GMPE is validated by comparing with three recently recorded earthquakes events. GMPE proposed in this study is capable of predicting PGA values close to recorded data and spectral acceleration up to period of 2 seconds. Comparison of new GMPE with the recorded data of recent earthquakes shows a good matching of ground motion as well as response spectra. The new GMPE is applicable for wide range of earthquake magnitudes from 5 to 9 on Mw scale. Reduction of future earthquake hazard is possible if hazard values are predicted precisely. A detailed seismic hazard analysis is carried out in this study considering deterministic and probabilistic approaches. New seismotectonic map has been generated for Lucknow considering a radial distance of 350 km around the city centre, which also covers active Himalayan plate boundaries. Past earthquakes within the seismotectonic region have been collected from United State Geological Survey (USGS), Northern California Earthquake Data Centre (NCEDC), Indian Meteorological Department (IMD), Seismic Atlas of India and its Environs (SEISAT) etc. A total of 1831 events with all the magnitude range were obtained. Collected events were homogenized, declustered and filtered for Mw ≥ 4 events. A total of 496 events were found within the seismic study region. Well delineated seismic sources are compiled from SEISAT. Superimposing the earthquake catalogue on the source map, a seismotectonic map of Lucknow was generated. A total of 47 faults which have experienced earthquake magnitude of 4 and above are found which are used for seismic hazard analysis. Based on the distribution of earthquake events on the seismotectonic map, two regions have been identified. Region I which shows high density of seismic events in the area in and around of Main Boundary Thrust (MBT) and Region II which consists of area surrounding Lucknow with sparse distribution of earthquake events. Data completeness analysis and estimation of seismic parameter “a” and “b” are carried out separately for both the regions. Based on the analysis, available earthquake data is complete for a period of 80 years in both the regions. Using the complete data set, the regional recurrence relations have been developed. It shows a “b” value of 0.86 for region I and 0.9 for Region II which are found comparable with earlier studies. Maximum possible earthquake magnitude in each source has been estimated using observed magnitude and doubly truncated Gutenberg-Richter relation. The study area of Lucknow is divided into 0.015o x 0.015o grid size and PGA at each grid has been estimated by considering all sources and the three GMPEs. A Matlab code was generated for seismic hazard analysis and maximum PGA value at each grid point was determined and mapped. Deterministic seismic hazard analysis (DSHA) shows that maximum expected PGA values at bedrock level varies from 0.05g in the eastern part to 0.13g in the northern region. Response spectrum at city centre is also developed up to a period of 2 seconds. Further, Probabilistic seismic hazard analysis (PSHA) has been carried out and PGA values for 10 % and 2 % probability of exceedence in 50 years have been estimated and mapped. PSHA for 10 % probability shows PGA variation from 0.035g in the eastern parts to 0.07g in the western and northern parts of Lucknow. Similarly PSHA for 2 % probability of exceedence indicates PGA variation from 0.07g in the eastern parts while the northern parts are expecting PGA of 0.13g. Uniform hazard spectra are also developed for 2 % and 10 % probability for a period of up to 2 seconds. The seismic hazard analyses in this study show that the northern and western parts of Lucknow are more vulnerable when compared to other part. Bedrock hazard values completely change due to subsoil properties when it reaches the surface. A detailed geophysical and geotechnical investigation has been carried out for subsoil profiling and seismic site classification. The study area has been divided into grids of 2 km x 2 km and roughly one geophysical test using MASW (Multichannel Analysis Surface Wave) has been carried out in each grid and the shear wave velocity (SWV) profiles of subsoil layers are obtained. A total of 47 MASW tests have been carried out and which are uniformly distributed in Lucknow. In addition, 12 boreholes have also been drilled with necessary sampling and measurement of N-SPT values at 1.5 m interval till a depth of 30 m. Further, 11 more borelog reports are collected from the same agency hired for drilling the boreholes. Necessary laboratory tests are conducted on disturbed and undisturbed soil samples for soil classification and density measurement. Based on the subsoil informations obtained from these boreholes, two cross-sections up to a depth of 30 m have been generated. These cross-sections show the presence of silty sand in the top 10 m at most of the locations followed by clayey sand of low to medium compressibility till a depth of 30 m. In between the sand and clay traces of silt were also been found in many locations. In addition to these boreholes, 20 deeper boreholes (depth ≥150 m) are collected from Jal Nigam (Water Corporation) Lucknow, Government of Uttar Pradesh. Typical cross-section along the alignment of these deeper boreholes has been generated up to 150 m depth. This cross-section shows the presence of fine sand near Gomati while other locations are occupied by surface clayey sand. Also, the medium sand has been found in the western part of the city at a depth of 110 m which continues till 150 m depth. On careful examination of MASW and boreholes with N-SPT, 17 locations are found very close and SWV and N-SPT values are available up to 30 m depth. These SWV and N-SPT values are complied and used to develop correlations between SWV and N-SPT for sandy soil, clayey soil and all soil types. This correlation is the first correlation for IGB soil deposits considered measured data up to 30 m. The new correlation is verified graphically using normal consistency ratio and standard percentage error with respect to measured N-SPT and SWV. Further, SWV and N-SPT profiles are used Another important earthquake induced hazard is liquefaction. Even though many historic earthquakes caused liquefaction in India, very limited attempt has been made to map liquefaction potential in IGB. In this study, a detailed liquefaction analysis has been carried out for Lucknow a part of Ganga Basin to map liquefaction potential. Initially susceptibility of liquefaction for soil deposits has been assessed by comparing the grain size distribution curve obtained from laboratory tests with the range of grain size distribution for potentially liquefiable soils. Most of surface soil deposits in the study area are susceptible to liquefaction. At all the 23 borehole locations, measured N-SPT values are corrected for (a) Overburden Pressure (CN), (b) Hammer energy (CE), (c) Borehole diameter (CB), (d) presence or absence of liner (CS), (e) Rod length (CR) and (f) fines content (Cfines). Surface PGA values at each borehole locations are used to estimate Cyclic Stress Ratio (CSR). Corrected N-SPT values [(N1)60CS] are used to estimate Cyclic Resistance Ratio (CRR) at each layer. CSR and CRR values are used to estimate Factor of Safety (FOS) against liquefaction in each layer. Least factor safety values are indentified from each location and presented liquefaction factor of safety map for average and maximum amplified PGA values. These maps highlight that northern, western and central parts of Lucknow are very critical to critical against liquefaction while southern parts shows moderate to low critical area. The entire alignment of river Gomati falls in very critical to critical regions for liquefaction. Least FOS shows worst scenario and does not account thickness of liquefiable soil layers. Further, these FOS values are used to determine Liquefaction Potential Index (LPI) of each site and developed LPI map. Based on LPI map, the Gomati is found as high to very high liquefaction potential region. Southern and the central parts of Lucknow show low to moderate liquefaction potential while the northern and western Lucknow has moderate to high liquefaction potential. All possible seismic hazards maps for Lucknow have been combined to develop final microzonation map in terms of hazard index values. Hazard index maps are prepared by combining rock PGA map, site classification map in terms of shear wave velocity, amplification factor map, and FOS map and predominant period map by adopting Analytical Hierarchy Process (AHP). All these parameters have been given here in the order starting with maximum weight of 6 for PGA to lower weight of 1 for predominant frequency. Normalized weights of each parameter have been estimated. Depending upon the variation of each hazard parameter values, three to five ranks are assigned and the normalized ranks are calculated. Final hazard index values have been estimated by multiplying normalized ranks of each parameter with the normalized weights. Microzonation map has been generated by mapping hazard index values. Three maps were generated based on DSHA, PSHA for 2% and 10 % probability of exceedence in 50 years. Hazard index maps from DSHA and PSHA for 2 % probability show similar pattern. Higher hazard index were obtained in northern and western parts of Lucknow and lower values in others. The new microzonation maps can help in dividing the Lucknow into three parts as high area i.e. North western part, moderate hazard area i.e. central part and low hazard area which covers southern and eastern parts of Lucknow. This microzonation is different from the current seismic code where all area is lumped in one zone without detailed assessment of different earthquake hazard parameters. Finally this study brings out first region specific GMPE considering recorded and synthetic ground monitions for wide range of magnitudes and distances. Proposed GMPE can also be used in other part of the Himalayan region as it matches well with the highly ranked GMPEs. Detailed rock level PGA map has been generated for Lucknow considering DSHA and PSHA. A detailed geotechnical and geophysical experiments are carried out in Lucknow. These results are used to develop correction between SWV and N-SPT values for soil deposit in IGB and site classification maps for the study area. Amplification and liquefaction potential of Lucknow are estimated by considering multiple ground motions data to account different earthquake ground motion amplitude, duration and frequency, which is unique in the seismic microzonation study.

Mankind is facing the problem due to earthquake hazard since prehistoric times. Many of the developed and developing countries are under constant threats from earthquakes hazards. Theories of plate tectonics and engineering seismology have helped to understand earthquakes and also to predicate earthquake hazards on a regional scale. However, the regional scale hazard mapping in terms of seismic zonation has been not fully implemented in many of the developing countries like India. Agglomerations of large population in the Indian cities and poor constructions have raised the risk due to various possible seismic hazards. First and foremost step towards hazard reduction is estimation of the seismic hazards in regional scale. Objective of this study is to estimate the seismic hazard parameters for Lucknow, a part of Indo-Gangetic Basin (IGB) and develop regional scale microzonation map. Lucknow is a highly populated city which is located close to the active seismic belt of Himalaya. This belt came into existence during the Cenozoic era (40-50 million years ago) and is a constant source of seismic threats. Many of the devastating earthquakes which have happened since prehistoric times such as 1255 Nepal, 1555 Srinagar, 1737 Kolkata, 1803 Nepal, 1833 Kathmandu, 1897 Shillong, 1905 Kangra, 1934 Bihar-Nepal, 1950 Assam and 2005 Kashmir. Historic evidences show that many of these earthquakes had caused fatalities even up to 0.1 million. At present, in the light of building up strains and non-occurrence of a great event in between 1905 Kangra earthquake and 1934 Bihar-Nepal earthquake regions the stretch has been highlighted as central seismic gap. This location may have high potential of great earthquakes in the near future. Geodetic studies in these locations indicate a possible slip of 9.5 m which may cause an event of magnitude 8.7 on Richter scale in the central seismic gap. Lucknow, the capital of Uttar Pradesh has a population of 2.8 million as per Census 2011. It lies in ZONE III as per IS1893: 2002 and can be called as moderate seismic region. However, the city falls within 350 km radial distance from Main Boundary Thrust (MBT) and active regional seismic source of the Lucknow-Faizabad fault. Considering the ongoing seismicity of Himalayan region and the Lucknow-Faizabad fault, this city is under high seismic threat. Hence a comprehensive study of understanding the earthquake hazards on a regional scale for the Lucknow is needed. In this work the seismic microzonation of Lucknow has been attempted. The whole thesis is divided into 11 chapters. A detailed discussion on the importance of this study, seismicity of Lucknow, and methodology adopted for detailed seismic hazard assessment and microzonation are presented in first three chapters. Development of region specific Ground Motion Prediction Equation (GMPE) and seismic hazard estimation at bedrock level using highly ranked GMPEs are presented in Chapters 4 and 5 respectively. Subsurface lithology, measurement of dynamic soil properties and correlations are essential to assess region specific site effects and liquefaction potential. Discussion on the experimental studies, subsurface profiling using geotechnical and geophysical tests results and correlation between shear wave velocity (SWV) and standard penetration test (SPT) N values are presented in Chapter 6. Detailed shear wave velocity profiling with seismic site classification and ground response parameters considering multiple ground motion data are discussed in Chapters 7 and 8. Chapters 9 and 10 present the assessment of liquefaction potential and determination of hazard index with microzonation maps respectively. Conclusions derived from each chapter are presented in Chapter 11. A brief summary of the work is presented below: Attenuation relations or GMPEs are important component of any seismic hazard analysis which controls accurate prediction of the hazard values. Even though the Himalayas have experienced great earthquakes since ancient times, suitable GMPEs which are applicable for a wide range of distance and magnitude are limited. Most of the available regional GMPEs were developed considering limited recorded data and/or pure synthetic ground motion data. This chapter presents development of a regional GMPE considering both the recorded as well as synthetic ground motions. In total 14 earthquakes consisting of 10 events with recorded data and 4 historic events with Isoseismal maps are used for the same. Synthetic ground motions based on finite fault model have been generated at unavailable locations for recorded events and complete range distances for historic earthquakes. Model parameters for synthetic ground motion were arrived by detailed parametric study and from literatures. A concept of Apparent Stations (AS) has been used to generate synthetic ground motion in a wide range of distance as well as direction around the epicenter. Synthetic ground motion data is validated by comparing with available recorded data and peak ground acceleration (PGA) from Isoseismal maps. A new GMPE has been developed based on two step stratified regression procedure considering the combined dataset of recorded and synthetic ground motions. The new GMPE is validated by comparing with three recently recorded earthquakes events. GMPE proposed in this study is capable of predicting PGA values close to recorded data and spectral acceleration up to period of 2 seconds. Comparison of new GMPE with the recorded data of recent earthquakes shows a good matching of ground motion as well as response spectra. The new GMPE is applicable for wide range of earthquake magnitudes from 5 to 9 on Mw scale. Reduction of future earthquake hazard is possible if hazard values are predicted precisely. A detailed seismic hazard analysis is carried out in this study considering deterministic and probabilistic approaches. New seismotectonic map has been generated for Lucknow considering a radial distance of 350 km around the city centre, which also covers active Himalayan plate boundaries. Past earthquakes within the seismotectonic region have been collected from United State Geological Survey (USGS), Northern California Earthquake Data Centre (NCEDC), Indian Meteorological Department (IMD), Seismic Atlas of India and its Environs (SEISAT) etc. A total of 1831 events with all the magnitude range were obtained. Collected events were homogenized, declustered and filtered for Mw ≥ 4 events. A total of 496 events were found within the seismic study region. Well delineated seismic sources are compiled from SEISAT. Superimposing the earthquake catalogue on the source map, a seismotectonic map of Lucknow was generated. A total of 47 faults which have experienced earthquake magnitude of 4 and above are found which are used for seismic hazard analysis. Based on the distribution of earthquake events on the seismotectonic map, two regions have been identified. Region I which shows high density of seismic events in the area in and around of Main Boundary Thrust (MBT) and Region II which consists of area surrounding Lucknow with sparse distribution of earthquake events. Data completeness analysis and estimation of seismic parameter “a” and “b” are carried out separately for both the regions. Based on the analysis, available earthquake data is complete for a period of 80 years in both the regions. Using the complete data set, the regional recurrence relations have been developed. It shows a “b” value of 0.86 for region I and 0.9 for Region II which are found comparable with earlier studies. Maximum possible earthquake magnitude in each source has been estimated using observed magnitude and doubly truncated Gutenberg-Richter relation. The study area of Lucknow is divided into 0.015o x 0.015o grid size and PGA at each grid has been estimated by considering all sources and the three GMPEs. A Matlab code was generated for seismic hazard analysis and maximum PGA value at each grid point was determined and mapped. Deterministic seismic hazard analysis (DSHA) shows that maximum expected PGA values at bedrock level varies from 0.05g in the eastern part to 0.13g in the northern region. Response spectrum at city centre is also developed up to a period of 2 seconds. Further, Probabilistic seismic hazard analysis (PSHA) has been carried out and PGA values for 10 % and 2 % probability of exceedence in 50 years have been estimated and mapped. PSHA for 10 % probability shows PGA variation from 0.035g in the eastern parts to 0.07g in the western and northern parts of Lucknow. Similarly PSHA for 2 % probability of exceedence indicates PGA variation from 0.07g in the eastern parts while the northern parts are expecting PGA of 0.13g. Uniform hazard spectra are also developed for 2 % and 10 % probability for a period of up to 2 seconds. The seismic hazard analyses in this study show that the northern and western parts of Lucknow are more vulnerable when compared to other part. Bedrock hazard values completely change due to subsoil properties when it reaches the surface. A detailed geophysical and geotechnical investigation has been carried out for subsoil profiling and seismic site classification. The study area has been divided into grids of 2 km x 2 km and roughly one geophysical test using MASW (Multichannel Analysis Surface Wave) has been carried out in each grid and the shear wave velocity (SWV) profiles of subsoil layers are obtained. A total of 47 MASW tests have been carried out and which are uniformly distributed in Lucknow. In addition, 12 boreholes have also been drilled with necessary sampling and measurement of N-SPT values at 1.5 m interval till a depth of 30 m. Further, 11 more borelog reports are collected from the same agency hired for drilling the boreholes. Necessary laboratory tests are conducted on disturbed and undisturbed soil samples for soil classification and density measurement. Based on the subsoil informations obtained from these boreholes, two cross-sections up to a depth of 30 m have been generated. These cross-sections show the presence of silty sand in the top 10 m at most of the locations followed by clayey sand of low to medium compressibility till a depth of 30 m. In between the sand and clay traces of silt were also been found in many locations. In addition to these boreholes, 20 deeper boreholes (depth ≥150 m) are collected from Jal Nigam (Water Corporation) Lucknow, Government of Uttar Pradesh. Typical cross-section along the alignment of these deeper boreholes has been generated up to 150 m depth. This cross-section shows the presence of fine sand near Gomati while other locations are occupied by surface clayey sand. Also, the medium sand has been found in the western part of the city at a depth of 110 m which continues till 150 m depth. On careful examination of MASW and boreholes with N-SPT, 17 locations are found very close and SWV and N-SPT values are available up to 30 m depth. These SWV and N-SPT values are complied and used to develop correlations between SWV and N-SPT for sandy soil, clayey soil and all soil types. This correlation is the first correlation for IGB soil deposits considered measured data up to 30 m. The new correlation is verified graphically using normal consistency ratio and standard percentage error with respect to measured N-SPT and SWV. Further, SWV and N-SPT profiles are used Another important earthquake induced hazard is liquefaction. Even though many historic earthquakes caused liquefaction in India, very limited attempt has been made to map liquefaction potential in IGB. In this study, a detailed liquefaction analysis has been carried out for Lucknow a part of Ganga Basin to map liquefaction potential. Initially susceptibility of liquefaction for soil deposits has been assessed by comparing the grain size distribution curve obtained from laboratory tests with the range of grain size distribution for potentially liquefiable soils. Most of surface soil deposits in the study area are susceptible to liquefaction. At all the 23 borehole locations, measured N-SPT values are corrected for (a) Overburden Pressure (CN), (b) Hammer energy (CE), (c) Borehole diameter (CB), (d) presence or absence of liner (CS), (e) Rod length (CR) and (f) fines content (Cfines). Surface PGA values at each borehole locations are used to estimate Cyclic Stress Ratio (CSR). Corrected N-SPT values [(N1)60CS] are used to estimate Cyclic Resistance Ratio (CRR) at each layer. CSR and CRR values are used to estimate Factor of Safety (FOS) against liquefaction in each layer. Least factor safety values are indentified from each location and presented liquefaction factor of safety map for average and maximum amplified PGA values. These maps highlight that northern, western and central parts of Lucknow are very critical to critical against liquefaction while southern parts shows moderate to low critical area. The entire alignment of river Gomati falls in very critical to critical regions for liquefaction. Least FOS shows worst scenario and does not account thickness of liquefiable soil layers. Further, these FOS values are used to determine Liquefaction Potential Index (LPI) of each site and developed LPI map. Based on LPI map, the Gomati is found as high to very high liquefaction potential region. Southern and the central parts of Lucknow show low to moderate liquefaction potential while the northern and western Lucknow has moderate to high liquefaction potential. All possible seismic hazards maps for Lucknow have been combined to develop final microzonation map in terms of hazard index values. Hazard index maps are prepared by combining rock PGA map, site classification map in terms of shear wave velocity, amplification factor map, and FOS map and predominant period map by adopting Analytical Hierarchy Process (AHP). All these parameters have been given here in the order starting with maximum weight of 6 for PGA to lower weight of 1 for predominant frequency. Normalized weights of each parameter have been estimated. Depending upon the variation of each hazard parameter values, three to five ranks are assigned and the normalized ranks are calculated. Final hazard index values have been estimated by multiplying normalized ranks of each parameter with the normalized weights. Microzonation map has been generated by mapping hazard index values. Three maps were generated based on DSHA, PSHA for 2% and 10 % probability of exceedence in 50 years. Hazard index maps from DSHA and PSHA for 2 % probability show similar pattern. Higher hazard index were obtained in northern and western parts of Lucknow and lower values in others. The new microzonation maps can help in dividing the Lucknow into three parts as high area i.e. North western part, moderate hazard area i.e. central part and low hazard area which covers southern and eastern parts of Lucknow. This microzonation is different from the current seismic code where all area is lumped in one zone without detailed assessment of different earthquake hazard parameters. Finally this study brings out first region specific GMPE considering recorded and synthetic ground monitions for wide range of magnitudes and distances. Proposed GMPE can also be used in other part of the Himalayan region as it matches well with the highly ranked GMPEs. Detailed rock level PGA map has been generated for Lucknow considering DSHA and PSHA. A detailed geotechnical and geophysical experiments are carried out in Lucknow. These results are used to develop correction between SWV and N-SPT values for soil deposit in IGB and site classification maps for the study area. Amplification and liquefaction potential of Lucknow are estimated by considering multiple ground motions data to account different earthquake ground motion amplitude, duration and frequency, which is unique in the seismic microzonation study.