Academic literature on the topic 'Discontinuite sismique'

Abstract From seismic tomography and seismic anisotropy, images of the lithosphere-asthenosphere system beneath France for some remarkable tectonic areas have been computed : a continental rift system (the Upper Rhinegraben), an Hercynian structure reactivated by Neogene volcanism (Massif central), a region of a recent continental collision (Pyrenees) and finally a region of an ancient orogeny (Armorican Massif). These images have a horizontal spatial resolution of the order of 10 km and show not only the geometry of the deep geological structures but will also illustrate the link between surface observations and structures detected at depth. The images demonstrate the passive character of the Rhinegraben mainly because no low-velocity was found below the Moho, show the presence of a thermal anomaly beneath the Massif central interpreted as caused by a mantle plume in the decaying phase of its evolution and prove the lithospheric scale of the North Pyrenean fault and of the South-Armorican shear zone. The anisotropic measurements suggest a lithospheric deformation related to the most recent tectonic event. In the Pyrenees, the Armorican Massif or the Rhinegraben areas, the directions of the fast-polarisation azimuth (the polarisation direction of the fast shear wave) are parallel to the tectonic texture of the last events, but suggest also a reactivation of inherited Hercynian discontinuities. In the Massif central, the splitting parameters distinguish between two lithospheric units regions marked by a distinct fast-polarisation azimuth on each side of the Sillon Houiller fault zone.

SUMMARYWhy is the North American Cordilleran mountain belt high? We expect a thick crust to support high elevations by isostasy but, remarkably, the Cordilleran crust is thin. There is no crustal root. An important recent recognition is that the high elevation is supported by thermal expansion rather than by thickened crust. The elevation of the Cordillera is only one consequence of the Cordillera being uniformly hot and having a thin lithosphere, in common with most current or recent backarcs. Some other consequences of the high temperatures compared to the adjacent cool craton include: (1) The Cordillera and other backarcs are hot, weak mobile belts that can be deformed by available plate-tectonic forces, in contrast to stable cratons that cannot; (2) Most continental seismicity is concentrated in backarcs; (3) In the Cordillera there is widespread sporadic ‘backarc’ volcanism; (4) The high temperatures result in very low strength in the lower crust that allows lower-crust detachment; (5) The lower crust weakness facilitates large-scale crustal oroclines that may be independent of the upper mantle; (6) The lower crust in the Cordillera and other backarcs is in amphibolite- to granulite-facies conditions, ~800–900°C at the Moho; (7) In ancient backarcs globally, regional Barrovian metamorphism is concluded to be the result of high temperatures that predate the orogenic collision and deformation. No "heat of orogeny" is required. Following the termination of subduction, backarcs cool with a time constant of 300–500 m.y.RÉSUMÉPourquoi la chaîne de montagnes de la Cordillère nord-américaine est-elle si haute? On comprend qu’une croûte sur-épaisse puisse expliquer une grande élévation, mais voilà, la croûte de la Cordillère est mince. Il n’existe pas de racine crustale. Or, récemment, une conclusion importante s’est imposée, soit que cette haute élévation s’explique par l’expansion thermique plutôt que par l’existence d’une croûte sur-épaisse. L’élévation de la Cordillère n’est qu’une des conséquences d’une Cordillère uniformément chaude flottant sur une lithosphère mince, caractéristiques communes aux zones d’arrière-arc actuelles ou récentes. Quelques unes des autres conséquences de cette haute température, par opposition aux froids cratons adjacents, comprennent: (1) La Cordillère et d’autres zones d’arrières-arcs sont des zones chaudes et facilement déformables par les forces tectoniques ambiantes, contrairement aux cratons stables; (2) La majorité de l’activité sismique continentale est concentrée dans le zones d’arrière-arc; (3) Dans la Cordillère l’activité volcanique sporadique est généralisé; (4) Ces températures élevées explique la très faible rigidité de la croûte inférieure et les décollements qu’elle subit; (5) La flaccidité de la croûte inférieure facilite la formation d’oroclinaux de grandes magnitudes qui peuvent être indépendants du manteau supérieur; (6) La croûte inférieure de la Cordillère et d’autres zones d’arrière-arc sont dans la zone de faciès amphibolite à granulite, soit 800 à 900oC à la discontinuité Moho; (7) Globalement dans les anciennes zones d’arrière-arc, le métamorphisme régional barrovien s’explique alors comme étant le résultat des hautes températures antérieures à la collision et à la déformation orogénique. Aucune « chaleur orogénique » n’est nécessaire. Après la période de subduction, les zones d’arrière-arc se refroidissent à l’intérieur d’un intervalle de temps de 300 à 500 millions d’années.

En etudiant la propagation des ondes de volume, nous pouvons d'une part caracteriser les variations brutales des proprietes physiques des materiaux et d'autre part evaluer l'anisotropie du manteau. Plusieurs traitements ont ete appliques aux fonctions de transfert obtenues a partir des enregistrements du reseau tgrs installe a partir de 1995 dans le sud des alpes francaises et la corse. Afin de completer les donnees, nous avons deploye des stations temporaires dans le cadre du programme geofrance 3d. L'etude conjointe de la forme des composantes longitudinales et transverses de ces fonctions a permis de mettre en evidence une croute anormalement fine, sous l'extreme sud-est de la france, accompagnee d'une geometrie complexe de l'interface croute/manteau. Nous montrons egalement les fortes variations de profondeur des principales discontinuites du manteau superieur entre le sud de l'europe occidentale et l'afrique du nord. Nous observons ainsi un approfondissement significatif de la discontinuite 410 entre les alpes maritimes et le sud algerien tandis que dans un meme temps, l'interface 670 s'eleve d'une trentaine de kilometres. Les donnees sismologiques acquises lors de la campagne broadves nous ont permis d'observer un pendage vers le nord-est du moho sous la region du vesuve (italie). D'autre part, un deuxieme reflecteur plus profond est interprete comme la trace fossile d'une subduction europeenne. Dans le sud des alpes occidentales, l'analyse de l'anisotropie du manteau est menee a partir de l'observation de la birefringence des ondes sks. Les resultats obtenus aux stations sont tres homogenes mais ne correspondent pas au mouvement absolu des plaques predit par le modele nnr-nuvel1. La direction de l'axe rapide n10\ refleterait ainsi la deformation de la lithosphere et non la direction d'ecoulement du manteau a sa base. Le delai moyen obtenu permet d'estimer l'epaisseur de la zone anisotrope a 130 kilometres.

L'objectif principal de cette thèse est la détection d'informations pertinentes dans les images sismiques. Ce type d'image est rencontré en particulier dans l'industrie pétrolière. Nous avons concentré notre travail sur la détection des horizons sismiques et des discontinuités en complément à la fusion des attributs sismiques. Nous avons ajusté l'utilisation des contours actifs aux images sismiques. Une nouvelle approche automatique a été introduite qui est basée sur l'optimisation d'une fonction énergie reliée à l'amplitude de l'image sismique. Nous avons introduit une généralisation de l'entropie cumulative résiduelle que nous avons appelée Entropie Cumulative Résiduelle Généralisée (GCRE) et qui à la différence de la CRE, introduite récemment, n'est pas restreinte aux variables aléatoires à valeur positive. Nous avons montré l'importance de la GCRE pour la détection d'objets saillants dans les images. Nous avons également montré l'utilité de l'attribut d'entropie pour l'expert. La sensibilité des deux informations mutuelles basées sur la GCRE et sur l'entropie de Shannon a aussi été étudiée dans le but d'utiliser ces mesures pour la comparaison d'images et la caractérisation des réservoirs. Nous avons traité la fusion des attributs sismiques pour intégrer le maximum d'information disponible en utilisant la théorie de Dempster-Shafer. Finalement, nous avons montré l'apport de cette nouvelle méthode de fusion sur des données réelles marines<br>The main objective of this thesis is the detection of useful information in seismic images. This type of images is found in particular in the petroleum industry. This work focuses on the detection of seismic horizons and discontinuities in addition to the fusion of seismic attributes. We adjusted the use of active contours to seismic images. A new automatic approach was introduced, based on the optimization of an energy function linked to the amplitude of the seismic image. We also introduced a generalization of the Cumulative Residual Entropy that we called Generalized Cumulative Residual Entropy (GCRE) which, unlike the recently introduced CRE, is not restricted to positive random variables. We showed the importance of the GCRE for detecting salient features in the images. Also, we showed the usefulness of the attribute entropy to the interpreter. The sensitivity of the mutual information measures based on the GCRE and on the Shannon entropy respectively to rigid transformations was studied in order to be able to use these measures in images comparison and reservoir characterization. We addressed the fusion of attributes to integrate as much available information as possible by using the Dempster-Shafer theory. Finally, we showed the strength of this new method of fusion on real marine data

Etude de petrologie experimentale sur les transformations a hautes pressions et hautes temperatures, de silicates et de germanates a structure olicine et pyroxene. Il a ete determine en fonction de la pression et de la temperature, la distribution du magnesium, du fer et du calcium entre mineraux mantelliques. Un modele thermodynamique coherent est construit et a permis de determiner le champ de stabilite de la perovskite silicatee et de proposer une interpretation de la discontinuite sismique des 670 km. Des mecanismes microscopiques de transformations de phase de l'olivine sont proposes. De quelques composes etudies par spectroscopie raman, a ete tire des parametres thermodynamiques harmoniques et anharmoniques de ces phases

Ce travail est consacré à l'étude des structures du manteau terrestre sousles Amériques en utilisant deux méthodes différentes de la sismologie: latomographie sismique et les fonctions récepteurs S. Dans la premièrepartie, nous avons obtenu par la tomographie les modèles de vitesses Vp etVs dans le manteau sous les Amériques à l'échelle régionale et ceux-cisont interprétés en les combinant avec les résultats de la reconstructionde l'histoire de mouvement des plaques. Nous avons réussit à identifierles structures des anomalies de vitesse rapide observées dans nos modèlestomographiques avec les différents fragments de plaques qui ont subductédans le manteau; nous avons pu ainsi suivre le fragmentation de la plaqueFarallon au cours du temps avec nos images tomographiques. Le bon accordentre les volumes et longueurs des plaques subductées estimés avec lesimages tomographiques et ceux calculés avec la reconstruction del'histoire des plaques confirme nos interprétations. Dans la secondepartie, nous avons utilisé la méthode de fonction récepteur S pour étudierles discontinuités du manteau supérieur sous le Nord-Ouest Américain. Lesstacks des fonctions récepteurs mettent en évidence clairement le Moho etla discontinuité à 410 km, cette dernière connaît une dépression d'environ15 km. La discontinuité de Lehmann à 220 km est observée dans la zone sous la plaque océanique suivant l'azimut autour de 236°, alors qu'elle est peuou pas visible suivant les azimuts autour de 136° et 306°. Enfin, nousobservons également une zone de faible vitesse au-dessus de la "410" et lasommation CCP des fonctions récepteurs montre que sa profondeur se situeautour de 380 km<br>We used seismic tomography and S receiver functions to study structures ofthe Earth's manteau beneath the Americas. In the first part, we combined Pand S-wave velocity models provided by the seismic tomography at regionalscale and plate motion history to investigate slabs of subductedlithosphere in the lower mantle beneath the Americas. We have been able toidentify different fast velocity anomalies observed in our tomographicimages with fragments of slabs that have subducted into the mantle; inparticular, we can follow the fragmentation of the Farallon plate withtime. Our interpretations are consistent with the volume and length ofslabs estimated from the tomographic images and those calculated from theplate history reconstruction. In the second part, we have used S receiverfunction technique to investigate seismic discontinuities in the uppermantle beneath the North-Western America. Stacks of receiver functionsshowed clear evidences of the Moho and the "410-km" discontinuity which isprobably depressed by about 15 km. The Lehmann discontinuity at 220 km isobserved in the area under the oceanic plate in the azimuth range around236°, whereas there is no clear evidence of this discontinuity for theazimuth ranges around 136° and 306°. Finally, we detect also a lowvelocity layer atop the "410-km" discontinuity and the stacking at CommunConversion Point show that its depth is around 380 km

Pouvoir déterminer la structure et la composition de l’intérieur de la Terre est un enjeu scientifique fondamental, pour la compréhension de l’organisation de la Terre profonde, des mécanismes des séismes et leur localisation en lien avec la prévention du risque sismique, pour la détection et l’exploitation des ressources naturelles telles que l’eau ou les hydrocarbures, ou encore pour toutes les activités de construction et de prévention associées au génie civil. Pour cela, les ondes sismiques sont un outil de choix. L’utilisation d’une approximation haute fréquence pour la modélisation de la propagation des ondes est avantageuse en termes de coût de calcul dès lors que plusieurs centaines, voire milliers, ou plus de longueurs d’ondes doivent être propagées. À la place de l’équation des ondes linéaire, l’approximation haute fréquence fournit trois équations aux dérivées partielles fondamentales. L’équation Eikonal, non linéaire, permet d’obtenir le temps de trajet. Une deuxième équation fournit l’angle d’émergence. L’équation Eikonal et l’équation des angles appartiennent toutes deux à la grande famille des équations de Hamilton-Jacobi. Enfin, l’équation de transport permet de calculer l’amplitude.Le tracé des rais sismiques est une technique lagrangienne qui utilise la méthode des caractéristiques pour obtenir un ensemble d’équations différentielles ordinaires à partir de ces équations aux dérivées partielles. Ces équations peuvent être intégrées facilement, donnant ainsi accès au temps de trajet et à l’amplitude le long des rais. Très largement utilisés dans la communauté géophysique du fait de leur simplicité, les outils de tracé de rais ne sont pas pour autant les plus efficaces et les plus robustes en pratique pour des applications d’imagerie et d’inversion haute résolution. En lieu et place, il peut être utile de résoudre directement les équations aux dérivées partielles par une méthode eulérienne. Durant les trois dernières décennies, une multitude de solveurs ont été développés pour l’équation Eikonal, la plupart utilisant la méthode des différences finies. Ces différents travaux visent à obtenir le meilleur compromis entre précision, coût de calcul, robustesse, facilité d’implémentation et souplesse d’utilisation.Dans cette thèse, je développe une approche différente, se basant principalement sur une méthode Galerkine discontinue. Dans le champ des mathématiques, cette méthode a été largement utilisée pour résoudre les lois de conservation et les équations de Hamilton-Jacobi. Très peu de travaux ont porté sur l'utilisation de cette méthode pour la résolution de l’équation Eikonal statique dans un contexte géophysique, et ce malgré le haut niveau de précision qu'elle apporte. C’est pourquoi, en me basant sur des travaux mathématiques, je propose un nouveau solveur Eikonal adapté au contexte géophysique. Les milieux hétérogènes complexes, anisotropes, et incluant des variations topographiques sont correctement pris en compte, avec une précision sans précédent. En y intégrant de manière robuste une stratégie de balayage rapide, je montre que ce solveur présente une très grande efficacité en deux comme en trois dimensions.J'utilise également ce solveur pour calculer l’angle d’émergence. Je développe par ailleurs un solveur voisin en volumes finis pour la résolution de l’équation de transport, permettant ainsi le calcul de l’amplitude. La variable d’état adjoint pour la tomographie sismique des temps et des pentes vérifiant une équation de transport semblable, je montre qu'on peut également la calculer à l'aide de ce solveur en volumes finis. En conséquence, je propose et analyse un ensemble consistant de solveurs pour la communauté géophysique. Ces outils devraient s’avérer utiles pour une large palette d’applications. Finalement, en guise d’illustration, je les utilise dans des schémas d’imagerie sismique, dans le but de démontrer le bénéfice apporté par une approximation haute fréquence dans ce type de schémas<br>Recovering information on the structure and the composition of the Earth's interior is a fundamental issue for a large range ofapplications, from planetology to seismology, natural resources assessment, and civil engineering. Seismic waves are a very powerful tool for that purpose. Using a high-frequency approximation for the numerical modeling of seismic wave propagation is computationally advantageous when hundreds, thousands, or more of wavelengths have to be propagated. Instead of the linear wave equation, the high-frequency approximation yields three fundamental partial differential equations. The nonlinear Eikonal equation leads to traveltime. A second equation is derived for the take-off angle. Both Eikonal and angle equations belong to the wide Hamilton-Jacobi family of equations. In addition, the transport equation leads to the amplitude.As a Lagrangian approach, seismic ray tracing employs the method of characteristics to derive a set of ordinary differential equations from these partial differential equations. They can be easily integrated, thus yielding traveltime and amplitude along rays. Widely used in the geophysical community for their simplicity, the ray-tracing tools might not be the most efficient and robust ones for practical high-resolution imaging and inversion applications. Instead, it might be desirable to directly solve the partial differential equations in an Eulerian way. In the three last decades, plenty of Eikonal solvers have been designed, mostly based on finite-difference methods. Successive works try to find the best compromise between accuracy, computational efficiency, robustness, ease of implementation, and versatility.In this thesis, I develop a different approach, mainly based on the discontinuous Galerkin method. This method has been intensively used in the mathematical field for solving conservation laws and time-dependent Hamilton-Jacobi equations. Only few investigations have been done regarding its use for solving the static Eikonal equation in a geophysical context, despite the high level of accuracy allowed by this method. Therefore, improving upon mathematical studies, I propose a new Eikonal solver suitable for the geophysical context. Complex heterogeneous anisotropic media with non-flat topographies are correctly handled, with an unprecedented accuracy. Combined with a fast-sweeping strategy in a robust way, I show that this new solver exhibits a high computational efficiency, in two dimensions as well as in three dimensions.I also employ this solver for the computation of the take-off angle. I design an additional finite-volume solver for solving the transport equation, leading to the computation of amplitude. With this solver, I also consider the computation of the adjoint-state variable for seismic tomography, since it satisfies a similar transport equation. Eventually, I propose a whole set of consistent solvers to the geophysical community. These tools should be useful in a wide range of applications. As an illustration, I finally use them in advanced seismic imaging schemes, in order to demonstrate the benefit brought by the high-frequency approximation in this kind of schemes

Le but de ce travail est de comprendre le fonctionnement de trois zones de failles anciennes et le rôle joué par les fluides au sein de ces structures. La première faille est une faille de détachement située dans le Jura. Elle a d'abord fonctionné dans des conditions de déformation ductile puis cassante, lors de l'effondrement de la chaîne Varisque. Peu de fluides ont circulé lors des premières phases de fonctionnement. En revanche, d'importantes circulations de fluides ont eu lieu lors de la phase de déformation cassante. Les fluides chargés en silice se sont retrouvés piégés le long du plan de faille supérieur. Il s'en est suivi des surpressions de fluides qui ont fracturé la roche et ont permis le glissement le long de cette faille faiblement pentée. Les deux autres failles sont des décrochements situés au Chugoku occidental (Japon). Elles jouent l erole de conduit pour les circulations de fluides dans le cas du décrochement de Kake, les fluides facilitent le glissement. Pour celui de Hikimi, des surpressions de fluides fractures la roche et les fractures sont colmatées par le la laumontite. Les trois zones de failles étudiées montrent différents modes d'interactions fluide-roche qui sont fonction de leur géométrie et de leur structure interne. La caractérisation structurale in-situ de ce type d'objets aide à la compréhension des données issues d'observations indirectes ou de forages sur des zones de failles actives<br>The aim of this study is to understand the functioning of three ancient fault zones and the role played by fluids in these structures. The first studied fault is a detachment fault located in the Jura, near the city of Dôle. It first acted in ductile conditions and then in brittle conditions during the collapse of the Variscan chain. A limited quantity of fluid circulated during the first stages of deformation. On the other hand, a significant amount of fluids flowed during the brittle deformation stage. Silica-rich fluids were allowed trapped along the upper fault plane. Then, fluid overpressures fractured the rock and allowed slip along this low-angle fault. The two other studied faults are strike-slip faults located in the western Chugoku region (Japan). They both acted as conduits for fluid fractures were sealed with laumonite. The three fault zones show different fluid-rock interaction characterization of these kinds of structures helps understanding data obtained from indirect observation or from cores drilled across active fault zones

Nous nous intéressons à l'application de la méthode de Galerkine discontinue avec pénalité intérieure (IPDGM) pour la migration terrestre par équations d'ondes. Il s'agit d'utiliser la propagation d'ondes sismiques pour produire des images du sous-sol et notre contribution se situe au niveau de la résolution du problème direct que nous proposons de faire par une méthode d'éléments finis s'appuyant sur des tétraèdres. Cette méthode est particulièrement bien adaptée pour la migration terrestre où il est nécessaire de prendre en compte la topographie du site qui est de plus en plus souvent accidentée. Tout d'abord, nous analysons les performances de la méthode d'éléments finis que nous avons choisie et concluons qu'elle surpasse les méthodes de différences finies classiquement utilisées pour la Reverse Time Migration (RTM) et qu'elle est aussi performante en terme de précision que la méthode d'élément spectraux qui est la méthode d'éléments finis la plus populaire en géophysique. Nous justifions ensuite notre choix en montrant qu'on peut imager des régions complexes en utilisant notamment des données synthétiques générées avec un autre code que le code IPDGFem développé durant la thèse. Pour finir, nous nous intéressons à la technique de pas de temps locaux qui permet d'adapter le pas de temps de la méthode à la taille des mailles utilisées et proposons une nouvelle discrétisation temporelle qui permet de combiner des ordres différents en temps et en espace. Les premiers résultats obtenus en 1D sont très intéressants et encourageants et le 2D est en cours de développement.

L'etude detaillee des discontinuites a 400 et 650 km de profondeur devrait permettre d'apporter des informations pour contraindre les modeles dynamiques du manteau terrestre. Pour cela, on analyse les phases qui sont reflechies ou converties aux discontinuites. On utilise une methode de sommation de sismogrammes ("stacking") pour amplifier ces phases. Cette technique est appliquee a des donnees large bande des reseaux nars (europe) et geoscope (global). Trois etudes sont presentees: (1) l'observation d'une phase convertie anormale sous l'ile de kerguelen montre que la discontinuite a 650 km aurait une topographie anormale sous les points chauds. Les effets de la geometrie de la discontinuite sur la phase pds sont examines. (2) les phases pdp permettent d'echantillonner les discontinuites en un point quelconque. Les phases reflechies vers 300 et 660 km de profondeur sous la siberie, ont ete identifiees grqce a la technique de "stacking". (3) l'etude simultanee de phases reflechies et de phases converties permet une etude locale detaillee des discontinuites. Cette methode est testee avec la discontinuite de mohorovocic sous le massif central

Le contexte scientifique de cette thèse est l'imagerie sismique dont le but est de reconstituer la structure du sous-sol de la Terre. Comme le forage a un coût assez élevé, l'industrie pétrolière s'intéresse à des méthodes capables de reconstituer les images de la structure terrestre interne avant de le faire. La technique d'imagerie sismique la plus utilisée est la technique de sismique-réflexion qui est basée sur le modèle de l'équation d'ondes. L'imagerie sismique est un problème inverse qui requiert de résoudre un grand nombre de problèmes directs. Dans ce contexte, nous nous intéressons dans cette thèse à la résolution du problème direct en régime harmonique, soit à la résolution des équations d'Helmholtz. L'objectif principal est de proposer et de développer un nouveau type de solveur élément fini (EF) caractérisé par un opérateur discret de taille réduite (comparée à la taille des solveurs déjà existants) sans pour autant altérer la précision de la solution numérique. Nous considérons les méthodes de Galerkine discontinues (DG). Comme les méthodes DG classiques sont plus coûteuses que les méthodes EF continues si l'on considère un même problème à cause d'un grand nombre de degrés de liberté couplés, résultat des approximations discontinues, nous développons une nouvelle classe de méthode DG réduisant ce problème : la méthode DG hybride (HDG). Pour valider l'efficacité de la méthode HDG proposée, nous comparons les résultats obtenus avec ceux obtenus avec une méthode DG basée sur des flux décentrés en 2D. Comme l'industrie pétrolière s'intéresse au traitement de données réelles, nous développons ensuite la méthode HDG pour les équations élastiques d'Helmholtz 3D<br>The scientific context of this thesis is seismic imaging which aims at recovering the structure of the earth. As the drilling is expensive, the petroleum industry is interested by methods able to reconstruct images of the internal structures of the earth before the drilling. The most used seismic imaging method in petroleum industry is the seismic-reflection technique which uses a wave equation model. Seismic imaging is an inverse problem which requires to solve a large number of forward problems. In this context, we are interested in this thesis in the modeling part, i.e. the resolution of the forward problem, assuming a time-harmonic regime, leading to the so-called Helmholtz equations. The main objective is to propose and develop a new finite element (FE) type solver characterized by a reduced-size discrete operator (as compared to existing such solvers) without hampering the accuracy of the numerical solution. We consider the family of discontinuous Galerkin (DG) methods. However, as classical DG methods are much more expensive than continuous FE methods when considering steady-like problems, because of an increased number of coupled degrees of freedom as a result of the discontinuity of the approximation, we develop a new form of DG method that specifically address this issue: the hybridizable DG (HDG) method. To validate the efficiency of the proposed HDG method, we compare the results that we obtain with those of a classical upwind flux-based DG method in a 2D framework. Then, as petroleum industry is interested in the treatment of real data, we develop the HDG method for the 3D elastic Helmholtz equations

Le dimensionnement d'ouvrages importants, comme les centrales nucléaires sous sollicitations vibratoires ou sismiques, requiert l'étude du comportement des sols de fondation sous l'effet de chargements dynamiques intenses. Dans le but de mieux connaître les phénomènes d'amortissement du sol, un code de calcul capable de simuler la propagation d'une onde sismique de cisaillement dans un milieu dissipatif a été développé. Ce programme, fondé sur un modèle non linéaire hystérétique utilisant les algorithmes de Newmark-Wilson, Newton-Raphson et une discrétisation spatiale à pas variables, s'affranchit des difficultés liées aux discontinuités accéléromètriques. Ces simulations pourraient permettre l'identification des paramètres du sol par comparaison avec des mesures in situ.