La liquéfaction des sols qui est déclenchée par des mouvements sismiques forts peut modifier la réponse d’un site. Ceci occasionne des dégâts importants dans les structures comme a été mis en évidence lors des tremblements de terre récents tels que celui de Christchurch, Nouvelle-Zélande et du Tohoku, Japon. L’évaluation du risque sismique des structures nécessite une modélisation robuste du comportement non linéaire de sols et de la prise en compte de l’interaction sol-structure (ISS). En général, le risque sismique est décrit comme la convolution entre l’aléa et la vulnérabilité du système. Cette thèse se pose comme une contribution à l’étude, via une modélisation numérique, de l’apparition de la liquéfaction et à l’utilisation des méthodes pour réduire les dommages induits.A cet effet, la méthode des éléments finis(FEM) dans le domaine temporel est utilisée comme outil numérique. Le modèle principal est composé d’un bâtiment fondé sur un sable liquéfiable. Comme la première étape de l’analyse du risque sismique, la première partie de cette thèse est consacrée à la caractérisation du comportement du sol et à sa modélisation.Une attention particulière est donnée à la sensibilité du modèle à des paramètres numériques. En suite, le modèle est validé pour le cas d’une propagation des ondes 1D avec les mesures issus du benchmark international PRENOLIN sur un site japonais. D’après la comparaison, le modèle arrive à prédire les enregistrements dans un test en aveugle.La deuxième partie, concerne la prise en compte dans la modélisation numérique du couplage de la surpression interstitielle (Δpw)et de la déformation du sol. Les effets favorables ou défavorables de ce type de modélisation ont été évalués sur le mouvement en surface du sol lors de la propagation des ondes et aussi sur le tassement et la performance sismique de deux structures.Cette partie contient des éléments d’un article publié dans Acta Geotechnica (Montoya-Noguera and Lopez-Caballero, 2016). Il a été trouvé que l’applicabilité du modèle dépend à la fois du niveau de liquéfaction et des effets d’ISS.Dans la dernière partie, une méthode est proposée pour modéliser la variabilité spatiale ajoutée au dépôt de sol dû à l’utilisation des techniques pour diminuer le degré de liquéfaction. Cette variabilité ajoutée peut différer considérablement de la variabilité inhérente ou naturelle. Dans cette thèse, elle sera modélisée par un champ aléatoire binaire.Pour évaluer l’efficience du mélange, la performance du système a été étudiée pour différents niveaux d’efficacité, c’est-à-dire,différentes fractions spatiales en allant de non traitées jusqu’à entièrement traitées. Tout d’abord le modèle binaire a été testé sur un cas simple, tel que la capacité portante d’une fondation superficielle sur un sol cohérent.Après, il a été utilisé dans le modèle de la structure sur le sol liquéfiable. Ce dernier cas,en partie, a été publié dans la revue GeoRisk (Montoya-Noguera and Lopez-Caballero,2015). En raison de l’interaction entre les deux types de sols du mélange, une importante variabilité est mise en évidence dans la réponse de la structure. En outre, des théories classiques et avancées d’homogénéisation ont été utilisées pour prédire la relation entre l’efficience moyenne et l’efficacité. En raison du comportement non linéaire du sol, les théories traditionnelles ne parviennent pas à prédire la réponse alors que certaines théories avancées qui comprennent la théorie de la percolation peuvent fournir une bonne estimation. En ce qui concerne l’effet de la variabilité spatiale ajoutée sur la diminution du tassement de la structure, différents séismes ont été testés et la réponse globale semble dépendre de leur rapport de PHV et PHA. / Strong ground motions can trigger soil liquefaction that will alter the propagating signal and induce ground failure. Important damage in structures and lifelines has been evidenced after recent earthquakes such as Christchurch, New Zealand and Tohoku, Japanin 2011. Accurate prediction of the structures’ seismic risk requires a careful modeling of the nonlinear behavior of soil-structure interaction (SSI) systems. In general, seismic risk analysisis described as the convolution between the natural hazard and the vulnerability of the system. This thesis arises as a contribution to the numerical modeling of liquefaction evaluation and mitigation.For this purpose, the finite element method (FEM) in time domain is used as numerical tool. The main numerical model consists of are inforced concrete building with a shallow rigid foundation standing on saturated cohesionless soil. As the initial step on the seismic risk analysis, the first part of the thesis is consecrated to the characterization of the soil behavior and its constitutive modeling. Later on, some results of the model’s validation witha real site for the 1D wave propagation in dry conditions are presented. These are issued from the participation in the international benchmark PRENOLIN and concern the PARI site Sendaiin Japan. Even though very few laboratory and in-situ data were available, the model responses well with the recordings for the blind prediction. The second part, concerns the numerical modeling of coupling excess pore pressure (Δpw) and soil deformation. The effects were evaluated on the ground motion and on the structure’s settlement and performance. This part contains material from an article published in Acta Geotechnica (Montoya-Noguera andLopez-Caballero, 2015). The applicability of the models was found to depend on both the liquefaction level and the SSI effects.In the last part, an innovative method is proposed to model spatial variability added to the deposit due to soil improvement techniques used to strengthen soft soils and mitigate liquefaction. Innovative treatment processes such as bentonite permeations and biogrouting,among others have recently emerged.However, there remains some uncertainties concerning the degree of spatial variability introduced in the design and its effect of the system’s performance.This added variability can differ significantly from the inherent or natural variability thus, in this thesis, it is modeled by coupling FEM with a binary random field. The efficiency in improving the soil behavior related to the effectiveness of the method measured by the amount of soil changed was analyzed. Two cases were studied: the bearing capacity of a shallow foundation under cohesive soil and the liquefaction-induced settlement of a structure under cohesionless loose soil. The latter, in part, contains material published in GeoRisk journal (Montoya-Noguera and Lopez-Caballero, 2015). Due to the interaction between the two soils, an important variability is evidenced in the response. Additionally, traditional and advanced homogenization theories were used to predict the relation between the average efficiency and effectiveness. Because of the nonlinear soil behavior, the traditional theories fail to predict the response while some advanced theories which include the percolation theory may provide a good estimate. Concerning the effect of added spatial variability on soil liquefaction, different input motions were tested and the response of the whole was found to depend on the ratio of PHV and PHA of the input motion.

L'importance croissante de l'évaluation de la performance des structures soumis au chargement sismique souligne la nécessité d'estimer le risque de liquéfaction. Dans ce scénario extrême de la liquéfaction du sol, des conséquences dévastatrices sont observées, par exemple des tassements excessifs et des instabilités de pentes. Dans le cadre de cette thèse, la réponse dynamique et l'interaction d'un système ouvrage en terre-fondation sont étudiées, afin de déterminer quantitativement le mécanisme de ruine dû à la liquéfaction du sol de la fondation. Par ailleurs, les chargements sismiques peuvent induire dans les ouvrages en terre un mode de rupture générant des bandes de cisaillement. Une étude de sensibilité aux maillages a donc été engagée pour quantifier la dépendance des résultats de l'analyse dynamique. Par conséquent, l'utilisation d'une méthode de régularisation est évaluée au cours des analyses dynamiques. Le logiciel open-source Code_Aster, basé sur la méthode des Eléments Finis et développé par EDF R&D, est utilisé pour les simulations numériques, tandis que le comportement du sol est représenté par le modèle de comportement de l'ECP, développé à CentraleSupélec. En premier lieu, un modèle simplifié de propagation 1D des ondes SH dans une colonne de sol avec comportement hydromécanique couplé non linéaire a été simulé. L'effet des caractéristiques du signal sismique et de la perméabilité du sol sur la liquéfaction est évalué. Le signal sismique d'entrée est un élément important pour l'apparition de la liquéfaction, puisque la durée du choc principal peut conduire à de fortes non linéarités et à un état de liquéfaction étendu. En outre, quand une variation de perméabilité en fonction de l'état de liquéfaction est considérée, des changements significatifs sont observés pendant la phase de dissipation de la surpression interstitielle de l'eau et au comportement du matériau. En revanche, ces changements ne suivent pas une tendance unique. Puis, l'effet d'une méthode de régularisation avec cinématique enrichie, appelée premier gradient de dilatation, sur la propagation des ondes SH est étudié au travers d'une solution analytique. Des problèmes à la réponse dynamique du sol sont observés et discutés quand cette méthode de régularisation est appliquée. Ensuite, un modèle 2D d'un déblai est simulé et sa réponse dynamique est évaluée en conditions sèches, complètement drainées et hydromécanique couplées. Deux critères sont utilisés pour définir le début de la rupture de la structure. Le travail du second ordre est utilisé pour décrire l'instabilité locale à des instants spécifiques du mouvement sismique, tandis que l'estimation d'un facteur de sécurité locale est proposée prenant en compte la résistance résiduelle du sol. En ce qui concerne le mode de ruine, l'effet de la surpression interstitielle de l'eau est de grande importance, puisqu'un déblai stable en conditions sèches et complètement drainées, devient instable lors de l'analyse couplée à cause de la liquéfaction de la fondation. Enfin, un système digue-fondation est simulé et l'influence de la perméabilité du sol, la profondeur de la couche liquéfiable, ainsi que, les caractéristiques du séisme sur la ruine induite par la liquéfaction du sol est évaluée. Pour ce modèle de digue, le niveau de dommages est fortement lié à la fois à l'apparition de la liquéfaction dans la fondation et la dissipation de la surpression d'eau. Une surface d'effondrement circulaire est générée à l'intérieur de la couche du sol liquéfié et se propage vers la crête dans les deux côtés de la digue. Pourtant, lorsque la couche liquéfiée est située en profondeur, la digue n'est pas affectée par la liquéfaction de la fondation pour ce cas particulier de chargement. Ce travail de recherche se concentre sur une étude de cas de référence pour l'évaluation sismique des ouvrages en terre soumis à un séisme et fournit des méthodes et outils de calculs numériques performants accessibles aux ingénieurs. / The increasing importance of performance-based earthquake engineering analysis points out the necessity to assess quantitatively the risk of liquefaction. In this extreme scenario of soil liquefaction, devastating consequences are observed, e.g. excessive settlements, lateral spreading and slope instability. The present PhD thesis discusses the global dynamic response and interaction of an earth structure-foundation system, so as to determine quantitatively the collapse mechanism due to foundation’s soil liquefaction. As shear band generation is a potential earthquake-induced failure mode in such structures, the FE mesh dependency of results of dynamic analyses is thoroughly investigated and an existing regularization method is evaluated. The open-source FE software developed by EDF R&D, called Code_Aster is used for the numerical simulations, while soil behavior is represented by the ECP constitutive model, developed at CentraleSupélec. Starting from a simplified model of 1D SH wave propagation in a soil column with coupled hydromechanical nonlinear behavior, the effect of seismic hazard and soil’s permeability on liquefaction is assessed. Input ground motion is a key component for soil liquefaction apparition, as long duration of mainshock can lead to important nonlinearity and extended soil liquefaction. Moreover, when a variation of permeability as function of liquefaction state is considered, changes in the dissipation phase of excess pore water pressure and material behavior are observed, which do not follow a single trend. The effect of a regularization method with enhanced kinematics approach, called first gradient of dilation model, on 1D SH wave propagation is studied through an analytical solution. Deficiencies of the use of this regularization method are observed and discussed, e.g. spurious waves apparition in the soil’s seismic response. Next, a 2D embankment-type model is simulated and its dynamic response is evaluated in dry, fully drained and coupled hydromechanical conditions. Two criteria are used to define the onset of the structure’s collapse. The second order work is used to describe the local instability at specific instants of the ground motion, while the estimation of a local safety factor is proposed by calculating soil’s residual strength. Concerning the failure mode, the effect of excess pore water pressure is of great importance, as an otherwise stable structure-foundation system in dry and fully drained conditions becomes unstable during coupled analysis. Finally, a levee- foundation system is simulated and the influence of soil’s permeability, depth of the liquefiable layer, as well as, characteristics of input ground motion on the liquefaction-induced failure is evaluated. For the current levee model, its induced damage level (i.e. settlements and deformations) is strongly related to both liquefaction apparition and dissipation of excess pore water pressure on the foundation. A circular collapse surface is generated inside the liquefied region and extends towards the crest in both sides of the levee. Even so, when the liquefied layer is situated in depth, no effect on the levee response is found. This research work can be considered as a reference case study for seismic assessment of embankment-type structures subjected to earthquake and provides a high-performance computational framework accessible to engineers.

Cette thèse étudie les instabilités dynamiques des milieux granulaires modèles saturés à l’aide d’un appareil triaxial classique. Les instabilités englobent la liquéfaction et les effondrements en compression isotrope drainée, les frottements saccadés en compression triaxiale drainée. Ces instabilités apparaissent spontanément à des contraintes effectives de confinement imprévisibles. Elles sont accompagnées de très rapides et très fortes surpressions interstitielles, malgré un drainage approprié; ce que ne présentent pas les milieux granulaires naturels. En compression isotrope drainée (consolidation), des effondrements locaux naissent instantanément. En compression triaxiale drainée, on observe de larges frottements saccadés quasi-périodiques caractérisés par des déformations volumiques et axiales contractantes. De temps en temps, ces effondrements et frottements saccadés locaux peuvent se développer en liquéfaction menant à une destruction complète de la structure granulaire. Les données à haute résolution temporelle issues de ce travail ont permis la découverte d’une nouvelle famille de liquéfaction dynamique et statique. L’étude des émissions acoustiques passives a permis l’identification de signature spectrale caractéristique. Pour les frottements saccadés, la phase de glissement peut être interprétée comme une consolidation dynamique, limitée par l’unique surface en dessous de la ligne critique de rupture dans le plan des contraintes effectives. La séquence temporelle précise des événements exclut que la pression interstitielle soit la cause principale des instabilités. Cependant, le rôle important de la surpression interstitielle est démontré dans des relations quantitatives entre les incréments de contraintes, et de déformations et l’éphémère surpression interstitielle stabilisée développée pendant la phase de glissement. Cela montre finalement la nature quasi-déterministique de ces instabilités dynamiques. Ces relations empiriques sont basées uniquement sur l’amplitude maximale de l’accélération verticale de très courte durée et sont gouvernées indépendamment par la pression de confinement et par l’indice des vides. La similarité de la surpression interstitielle entre différentes instabilités suggère fortement quelques mécanismes similaires de déclenchement, probablement à partir de ré-arrangements de la micro-structure granulaire. / This thesis reports a laboratory study on the dynamic instabilities of model saturated granular material using a triaxial apparatus. The term instability consists of isotropic collapse and liquefaction under isotropic compression and of stick-slip under triaxial compression in drained condition. The instabilities spontaneously occur at unpredictable effective stress with unexpected buildup of excess pore pressure irrespective of fully drained condition, contrasting with the instability-free behaviour of natural granular materials. In isotropic compression, instantaneous local collapse happens and in triaxial compression, very large and quasi-periodic stick-slip occurs with sudden volumetric compaction and axial contraction. Sometimes, these local failures (collapse and stick-slip) can develop into total liquefaction failure, destroying completely the granular structure. High time-resolved data permit the discovery of a new family of dynamic and static liquefaction. Passive acoustic measurements allow the identification of typical spectral signature. For stick-slip phenomenon, the slip phase with constant duration of stress drop can be interpreted as dynamic consolidation at constant deviatoric stress, limited by a unique boundary inside the critical state line in the effective stress plane. The precise temporal sequence of mechanical measurements excludes the generated pore pressure as the main cause of the instabilities. However, the role of pore pressure is emphasised by consistent quantitative relations between the amplitude of incremental stresses, incremental strains and the ephemeral stabilised excess pore pressure developed during the dynamic event, leading to the quasi-deterministic nature of granular instabilities. These empirical relations are based only on the short-lived maximum vertical acceleration and governed separately by the confining pressure and the initial void ratio. The similarity of pore pressure evolution for different kinds of instability strongly suggests some common speculative triggering mechanisms, probably originated from different rearrangements of the granular micro-structure.

No description available.

This thesis analyses the governing equation of elastic wave propagation by the finite difference method , and considered absorbing boundary condition and the material damping to simulate behavior of wave propagation. Otherwise, we combined with the mechanics of the soil pore water pressure raised by shear stress effected repeatedly and the soil property is changed by water pressure effected to simulate physical phenomenon in half-space, and probe into the soil liquefaction process during different force types. Using the developed numerical wave propagation model probe into reducing vibration by dug trench and filler trench, and analyzed data by 1/3 octave band method. This thesis discuss with reducing vibration effect by different trench disposed¡Bdifferent filler material property, complex filler, and extending the force source pile length.

Turkey is one of the most earthquake prone countries in the world. The study area, Erbaa, is located in a seismically active fault zone known as North Anatolian Fault Zone (NAFZ). Erbaa is one of the towns of Tokat located in the Middle Black Sea Region. According to the Earthquake zoning map of Turkey, the study area is in the First Degree Earthquake Zone. The city center of Erbaa (Tokat) was previously settled on the left embankment of Kelkit River. After the disastrous 1942 Niksar-Erbaa earthquake (Mw = 7.2), the settlement was moved southwards. From the period of 1900s, several earthquakes occurred in this region and around Erbaa. The 1942 earthquake is the most destructive earthquake in the center of Erbaa settlement. In this study, the geological and geotechnical properties of the study area were investigated by detailed site investigations. The Erbaa settlement is located on alluvial and Pliocene deposits. The Pliocene clay, silt, sand, and gravel layers exist in the southern part of Erbaa. Alluvium in Erbaa region consists of gravelly, sandy, silty, and clayey layers. The alluvial deposits are composed of stratified materials of heterogeneous grain sizes, derived from various geological units in the vicinity. The main objective of this study is to prepare a seismic microzonation map of the study area for urban planning purposes since it is getting more essential to plan new settlements considering safe development strategies after the disastrous earthquakes. In this respect, seismic hazard analyses were performed to deterministically assess the seismic hazard of the study area. Afterwards, the essential ground motions were predicted regarding near fault effects as the study area is settled on an active fault zone. 1-D equivalent linear site response analyses were carried out to evaluate the site effects in the study area. Amplification values obtained from site response analyses reveal that the soil layers in the study area is quite rigid. Furthermore, liquefaction potential and post liquefaction effects including lateral spreading and vertical settlement were also delineated for the study area. The above-mentioned parameters were taken into account in order to prepare a final seismic microzonation map of the study area. The layers were evaluated on the basis of overlay methodologies including Multi-Criteria Decision Analysis (MCDA). Two different MCDA techniques, Simple Additive Weighting (SAW) and Analytical Hierarchical Process (AHP), were carried out in GIS environment. The seismic microzonation maps prepared by SAW and AHP methods are compared to obtain a final seismic microzonation map. Finally, the map derived from the AHP method is proposed to be the final seismic microzonation map of Erbaa. As an overall conclusion, the northwestern part of the study area where the loose alluvial units exist is found to be vulnerable to earthquake-induced deformations. On the other hand, the Pliocene units in the southern and alluvial units in the northeastern part are quite resistant to earthquake effects. In addition, the proposed final seismic microzonation map should be considered by urban planners and policy makers during urban planning projects in Erbaa.

Although silt and clay mixtures were mostly considered to be resistant to cyclic loading due to cohesional components of their shear strength, ground failure case histories compiled from fine grained soil profiles after recent earthquakes (e.g. 1994 Northridge, 1999 Adapazari, 1999 Chi-Chi) revealed that the responses of low plasticity silt and clay mixtures are also critical under cyclic loading. Consequently, understanding the cyclic response of these soils has become a recent challenge in geotechnical earthquake engineering practice. While most of the current attention focuses on the assessment of liquefaction susceptibility of fine-grained soils, it is believed that cyclic strain and strength assessments of silt and clay mixtures need to be also studied as part of complementary critical research components. Inspired by these gaps, a comprehensive laboratory testing program was designed. As part of the laboratory testing program 64 stress-controlled cyclic triaxial tests, 59 static strain-controlled consolidated undrained triaxial tests, 17 oedometer, 196 soil classification tests including sieve analyses, hydrometer, and consistency tests were performed. Additionally 116 cyclic triaxial test results were compiled from available literature. Based on this data probability-based semi-empirical models were developed to assess liquefaction susceptibility and cyclic-induced shear strength loss, cyclically-induced maximum shear, post-cyclic volumetric and residual shear strains of silt and clay mixtures. Performance comparisons of the proposed model alternatives were studied, and it is shown that the proposed models follow an unbiased trend and produce superior predictions of the observed laboratory test response. Superiority of the proposed alternative models was proven by relatively smaller model errors (residuals).

Seismic hazard and microzonation of cities enable to characterize the potential seismic areas that need to be taken into account when designing new structures or retrofitting the existing ones. Study of seismic hazard and preparation of geotechnical microzonation maps will provide an effective solution for city planning and input to earthquake resistant design of structures in an area. Seismic hazard is the study of expected earthquake ground motions at any point on the earth. Microzonation is the process of sub division of region in to number of zones based on the earthquake effects in the local scale. Seismic microzonation is the process of estimating response of soil layers under earthquake excitation and thus the variation of ground motion characteristic on the ground surface. Geotechnical site characterization and assessment of site response during earthquakes is one of the crucial phases of seismic microzonation with respect to ground shaking intensity, attenuation, amplification rating and liquefaction susceptibility. Microzonation mapping of seismic hazards can be expressed in relative or absolute terms, on an urban block-by-block scale, based on local soil conditions (such as soil types) that affect ground shaking levels or vulnerability to soil liquefaction. Such maps would provide general guidelines for integrated planning of cities and in positioning the types of new structures that are most suited to an area, along with information on the relative damage potential of the existing structures in a region. In the present study an attempt has been made to characterize the site and to study the seismic hazard analysis considering the local site effects and to develop microzonation maps for Bangalore. Seismic hazard analysis and microzonation of Bangalore is addressed in this study in three parts: In the first part, estimation of seismic hazard using seismotectonic and geological information. Second part deals about site characterization using geotechnical and shallow geophysical techniques. An area of 220 sq.km, encompassing Bangalore Municipal Corporation has been chosen as the study area in this part of the investigation. There were over 150 lakes, though most of them are dried up due to erosion and encroachments leaving only 64 at present in an area of 220 sq. km and emphasizing the need to study site effects. In the last part, local site effects are assessed by carrying out one-dimensional (1-D) ground response analysis (using the program SHAKE 2000) using both borehole SPT data and shear wave velocity survey data within an area of 220 sq. km. Further, field experiments using microtremor studies have also been carried out (jointly with NGRI) for evaluation of predominant frequency of the soil columns. The same has been assessed using 1-D ground response analysis and compared with microtremor results. Further, Seed and Idriss simplified approach has been adopted to evaluate the liquefaction susceptibility and liquefaction resistance assessment. Microzonation maps have been prepared for Bangalore city covering 220 sq. km area on a scale of 1:20000. Deterministic Seismic Hazard Analysis (DSHA) for Bangalore has been carried out by considering the past earthquakes, assumed subsurface fault rupture lengths and point source synthetic ground motion model. The seismic sources for region have been collected by considering seismotectonic atlas map of India and lineaments identified from satellite remote sensing images. Analysis of lineaments and faults help in understanding the regional seismotectonic activity of the area. Maximum Credible Earthquake (MCE) has been determined by considering the regional seismotectonic activity in about 350 km radius around Bangalore. Earthquake data are collected from United State Geological Survey (USGS), Indian Metrological Department (IMD), New Delhi; Geological Survey of India (GSI) and Amateur Seismic Centre (ASC), National Geophysical Research Institute (NGRI),Hyderabad; Centre for Earth Science Studies (CESS), Akkulam, Kerala; Gauribindanur (GB) Seismic station and other public domain sites. Source magnitude for each source is chosen from the maximum reported past earthquake close to that source and shortest distance from each source to Bangalore is arrived from the newly prepared seismotectonic map of the area. Using these details, and, attenuation relation developed for southern India by Iyengar and Raghukanth (2004), the peak ground acceleration (PGA) has been estimated. A parametric study has been carried out to find fault subsurface rupture length using past earthquake data and Wells and Coppersmith (1994) relation between the subsurface lengths versus earthquake magnitudes. Further seismological model developed by Boore (1983, 2003) SMSIM program has been used to generate synthetic ground motions from vulnerable sources identified in above two methods. From the above three approaches maximum PGA of 0.15g was estimated for Bangalore. This value was obtained for a maximum credible earthquake (MCE) having a moment magnitude of 5.1 from a source of Mandya-Channapatna-Bangalore lineament. Considering this lineament and MCE, a synthetic ground motion has been generated for 850 borehole locations and they are used to prepare PGA map at rock level. The past seismic data has been collected for almost 200 years from different sources such as IMD, BARC (Gauribidanur array), NGRI, CESS, ASC center, USGS, and other public domain data. The seismic data is seen to be homogenous for the last four decades irrespective of the magnitude. Seismic parameters were then evaluated using the data corresponding to the last four decades and also the mixed data (using Kijko’s analysis) for Bangalore region, which are found to be comparable with the earlier reported seismic parameters for south India. The probabilities of distance, magnitude and peak ground acceleration have been evaluated for the six most vulnerable sources using PSHA (Probabilistic Seismic Hazard Analysis). The mean annual rate of exceedance has been calculated for all the six sources at the rock level. The cumulative probability hazard curves have been generated at the bedrock level for peak ground acceleration and spectral acceleration. The spectral acceleration calculation corresponding to a period of 1sec and 5% damping are evaluated. For the design of structures, uniform hazard response spectrum (UHRS) at rock level is developed for the 5% damping corresponding to 10% probability of exceedance in 50 years. The peak ground acceleration (PGA) values corresponding to 10% probability of exceedance in 50 years are comparable to the PGA values obtained in deterministic seismic hazard analysis (DSHA) and higher than Global Seismic Hazard Assessment Program (GSHAP) maps of Bhatia et.al (1997) for the Indian shield area. The 3-D subsurface model with geotechnical data has been generated for site characterization of Bangalore. The base map of Bangalore city (220sq.km) with several layers of information (such as Outer and Administrative boundaries, Contours, Highways, Major roads, Minor roads, Streets, Rail roads, Water bodies, Drains, Landmarks and Borehole locations) has been generated. GIS database for collating and synthesizing geotechnical data available with different sources and 3-dimensional view of soil stratum presenting various geotechnical parameters with depth in appropriate format has been developed. In the context of prediction of reduced level of rock (called as “engineering rock depth” corresponding to about Vs > 700 m/sec) in the subsurface of Bangalore and their spatial variability evaluated using Artificial Neural Network (ANN). Observed SPT ‘N’ values are corrected by applying necessary corrections, which can be used for engineering studies such as site response and liquefaction analysis. Site characterization has also been carried out using measured shear wave velocity with the help of shear wave velocity survey using MASW. MASW (Multichannel Analysis of Surface Wave) is a geophysical method, which generates a shear-wave velocity (Vs) profile (i.e., Vs versus depth) by analyzing Raleigh-type surface waves on a multichannel record. MASW system consisting of 24 channels Geode seismograph with 24 geophones of 4.5 Hz capacity were used in this investigation. The shear wave velocity of Bangalore subsurface soil has been measured and correlation has been developed for shear wave velocity (Vs) with the standard penetration tests (SPT) corrected ‘N’ values. About 58 one-dimensional (1-D) MASW surveys and 20 two-dimensional (2-D) MASW surveys has been carried out with in 220 sq.km Bangalore urban area. Dispersion curves and shear velocity 1-D and 2-D have been evaluated using SurfSeis software. Using 1-dimensional shear wave velocity, the average shear wave velocity of Bangalore soil has been evaluated for depths of 5m, 10m, 15m, 20m, 25m and 30m (Vs30) depths. The sub soil classification has been carried out for local site effect evaluation based on average shear wave velocity of 30m depth (Vs30) of sites using NEHRP (National Earthquake Hazard Research Programme) and IBC (International Building Code) classification. Bangalore falls into site class D type of soil. Mapping clearly indicates that the depth of soil obtained from MASW is closely matching with the soil layers in the bore logs. The measured shear wave velocity at 38 locations close to SPT boreholes, which are used to generate the correlation between the shear wave velocity and corrected ‘N’ values using a power fit. Also, developed relationship between shear wave velocity and corrected ‘N’ values corresponds well with the published relationships of Japan Road Association. Bangalore city, a fast growing urban center, with low to moderate earthquake history and highly altered soil structure (due to large reclamation of land) is been the focus of this work. There were over 150 lakes, though most of them are dried up due to erosion and encroachments leaving only 64 at present in an area of 220 sq km. In the present study, an attempt has been made to assess the site response using geotechnical, geophysical data and field studies. The subsurface profiles of the study area within 220sq.km area was represented by 170 geotechnical bore logs and 58 shear wave velocity profiles obtained by MASW survey. The data from these geotechnical and geophysical technique have been used to study the site response. These soil properties and synthetic ground motions for each borehole locations are further used to study the local site effects by conducting one-dimensional ground response analysis using the program SHAKE2000. The response and amplification spectrum have been evaluated for each layer of borehole location. The natural period of the soil column, peak spectral acceleration and frequency at peak spectral acceleration of each borehole has been evaluated and presented as maps. Predominant frequency obtained from both methods is compared; the correlation between corrected SPT ‘N’ value and low strain shear modulus has been generated. The noise was recorded at 54 different locations in 220sq.km area of Bangalore city using L4-3D short period sensors (CMG3T) equipped with digital data acquisition system. Predominant frequency obtained from ground response studies and microtremor measurement is comparable. To study the liquefaction hazard in Bangalore, the liquefaction hazard assessment has been carried out using standard penetration test (SPT) data and soil properties. Factor of Safety against liquefaction of soil layer has been evaluated based on the simplified procedure of Seed and Idriss (1971) and subsequent revisions of Seed et al (1983, 1985), Youd et al (2001) and Cetin et al (2004). Cyclic Stress Ratio (CSR) resulting from earthquake loading is calculated by considering moment magnitude of 5.1 and amplified peak ground acceleration. Cyclic Resistant Ratio (CRR) is arrived using the corrected SPT ‘N’ values and soil properties. Factor of safety against liquefaction is calculated using stress ratios and accounting necessary magnitude scaling factor for maximum credible earthquake. A simple spread sheet was developed to carryout the calculation for each bore log. The factor of safety against liquefaction is grouped together for the purpose of classification of Bangalore (220 sq. km) area for a liquefaction hazards. Using 2-D base map of Bangalore city, the liquefaction hazard map was prepared using AutoCAD and Arc GIS packages. The results are grouped as four groups for mapping and presented in the form of 2-dimensional maps. Liquefaction possibilities are also assessed conducting laboratory cyclic triaxial test using undisturbed soil samples collected at few locations.

A new reliability-based probabilistic model is developed for the estimation of liquefaction-induced lateral ground spreading, taking into consideration the uncertainties within the model functional form and the descriptive variables as well. The new model is also introduced as performance-based probabilistic engineering tool.

L'objectif de cette thèse est d'étudier la valorisation de la liqueur noire non soufrée par deux procédés hydrothermaux : la gazéification en eau supercritique et la liquéfaction hydrothermale. Ceux-ci seront comparés au procédé actuel de valorisation (évaporation puis combustion dans chaudière Tomlinson) selon 3 critères : le rendement énergétique, la récupération du sodium et la production de molécules aromatiques biosourcées.Lors de la gazéification, il a été montré que la formation de gaz est compétitive à celle de char. Une chauffe rapide et des températures élevées vont favoriser le rendement gaz et donc le rendement énergétique. Cependant les rendements énergétiques sont plus faibles que le procédé actuel car la conversion des composés aromatiques provenant de la lignine est faible dans la gamme de température étudiée. Lors d'un procédé en continu, à plus haute température (700°C) avec une chauffe rapide, le rendement énergétique peut être le double au procédé actuel (simulé à l'équilibre thermodynamique). La préhydrolyse du bois et l'utilisation de bois de résineux vont défavoriser la conversion de la liqueur noire en gaz.La liquéfaction quant-à-elle permet la formation composés phénoliques et d'un biocrude dont la combustion permettant de meilleur rendements énergétique que le procédé actuel. En effet, la lignine de la liqueur noire est hydrolysée en fragments réactifs, pouvant être soit dégradés soit se recombiner pour former le biocrude. Cette dernière est favorisée par la présence des carbohydrates. L'utilisation de bois de feuillus et la préhydrolyse vont améliorer le rendement énergétique.La récupération du sodium est satisfaisante pour les deux procédés, validant la faisabilité de la substitution de la chaudière par ces procédés hydrothermaux. / This thesis aims to study sulfur free black liquor valorization through two hydrothermal processes: supercritical water gasification and hydrothermal liquefaction. These processes will be compared to the industrial process (evaporation and Tomlinson boiler) with 3 mains criteria: energetic yield, sodium recovery and phenolic molecules production.In supercritical conditions, gas formation is competitive with char formation. Fast heating and high temperature permit to increase gas yield, thus energetic yield. However, conversion of phenolic compounds from lignin is low below 500°C, leading to a lower energetic yield than reference. In a continuous process, at high temperatures (700°C) and fast heating, energetic yield should be 2 times higher than industrial process (simulation at thermodynamic equilibrium). Wood prehydrolysis and softwood lead to a lower conversion of black liquor.Hydrothermal liquefaction produces a biocrude which can be burnt and phenolic platform compounds. Indeed, lignin is depolymerized into reactive fragments which can be degraded into platform phenolic molecules. Moreover, the recombination of these fragments, leading to biocrude formation, is favored by the carbohydrates derivatives in black liquor. Wood prehydrolysis and hardwood lead to better energetic and phenolic molecules yields.Sodium recovery is satisfactory for both processes. Substitution of Tomlinson recovery by a hydrothermal process is then possible.