Cette thèse s’intéresse à divers aspects de l’atténuation sismique dans le manteau terrestre et aux implications de celle-ci quant à la structure de ce dernier. L’enjeu est de mieux comprendre les mécanismes d’atténuation ainsi que les mesures que l’on peut effectuer afin d’améliorer les modèles radiaux d’atténuation dont on dispose et in fine l’interprétation des modèles de tomographie. Je me suis concentrée sur deux exemples de mécanismes d’atténuation, appartenant à deux grands types d’atténuation : l’atténuation intrinsèque, liée à l’absorption par le milieu d’une partie de l’énergie sismique dissipée irréversiblement sous forme de chaleur, et l’atténuation extrinsèque, liée à la dispersion de cette énergie par le milieu. Dans le premier cas, j'ai regardé l’effet des transitions de phase sur l’atténuation des ondes sismiques. En appliquant un modèle thermomécanique développé par Ricard et al., 2009, pour prédire l’atténuation des ondes sismiques liée à la transition de phase uniquement et en comparant les valeurs obtenues aux mesures dont on dispose, j'ai pu contraindre la cinétique d’une transition de phase mantellique. Dans le second cas, j'ai testé l’effet de l’anisotropie comme mécanisme d'atténuation apparente, le but étant de trouver une distribution statistique d’orientation d’anisotropie pouvant reproduire la quasi-constance du facteur de qualité Q avec la fréquence, observée en sismologie et lors d’expériences de laboratoire (Knopoff, 1964), et aujourd’hui expliquée par un modèle ad-hoc seulement (Liu, 1976).Enfin, je me intéressée à mesurer cette atténuation sismique sur des enregistrements réels. Après avoir testé la méthode dite de la fréquence instantanée (Ford et al., 2012), je me suis concentrée sur deux régions, l’Amérique centrale et l’Alaska pour l'appliquer. Ces mesures sont ensuite interprétées en termes de modèle radial d’atténuation révélant un manteau inférieur hétérogène atténuant. Je montre aussi qu’une origine compositionnelle est la plus probable pour expliquer ces anomalies d’atténuation. / This thesis is devoted to various aspects of seismic attenuation in the Earth's mantle and the consequences on the mantle structure. The challenge is to better understand the attenuation mechanisms, as well as the measurements that can be done, in order to improve the published radial profiles of attenuation and in fine the interpretation of tomographic models.I focus on two examples of attenuation mechanisms, belonging to two kinds of attenuation: the intrinsic attenuation related to the absorption by the medium of a part of the seismic energy then irreversibly dissipated as heat, and the extrinsic attenuation related to the dispersion of the seismic energy by the medium. In the first case, I investigate the effect of phase transitions upon seismic attenuation. Applying the thermo-mechanical model developped by Ricard et al., 2012, to calculate the attenuation of seismic waves due to the phase transition only and comparing the obtained values to published measurements, I succeed in constraining the kinetics of a mantle phase transition. In the second case, I test the seismic anisotropy as a mechanism of extrinsic attenuation, the aim being to find a statistical distribution of anisotropy orientation and layer thicknesses that can reproduce the observed quasi-frequency independence of Q in seismology and laboratory experiments (Knopoff, 1964), and which is, today, only explained by an ad-hoc model (Liu, 1976).Finally, I was interested in measuring the seismic attenuation on real seismograms. After having tested the method of the instantaneous frequency (Ford et al., 2012), I applied it to seismic records sampling the mantle below Central America and Alaska. These measurements are then inverted for a radial profile of shear attenuation which reveals the existence of an attenuating zone in the lower mantle. I also show that these attenuation anomalies are likely to be of chemical origin.

Le Noyau Terrestre: étude sismologique de quelques structures majeures influençant la dynamique terrestre

Ces travaux de thèe portent sur le développement d’un outil d’analyse numérique dédié à l’ étude de nouveaux résonateursà quartz à ondes de volume et utilisant les éléments finis. Cette méthode de caractérisation permet la détermination deséléments du schéma électrique équivalent (résistance, inductance et capacité) d’une fréquence de résonance donnée ainsique son facteur de qualité, tout en prenant en compte dans le modèle la sensibilité du cristal de quartz à la températureet aux contraintes induites par le montage. Une étape de validation est d’abord réalisée afin de vérifier nos choix, enterme de modélisation et de calcul, en confrontant les données issues de la simulation aux mesures de résonateurs déjàexistants. Les trois dispositifs analysés (40 MHz, 10 MHz et 100 MHz) montrent une bonne concordance entre théorieet expérience. Pour obtenir de tels résultats, la structure de maintien est prise en compte et modélisée sous forme dezones d’amortissement de Rayleigh lorsque le piégeage de l’énergie n’est pas optimal (présence d’un mode de plaque).Un aspect important des résonateurs est ensuite étudié : le comportement en température. En effet, les contraintes dedilatation thermique ainsi que l’évolution des coefficients élastiques en fonction de la température induisent une dérivefréquentielle. La comparaison entre théorie et expérience nous permet de vérifier l’allure des courbes et de quantifier ledegré de précision du modèle. L’effet d’une contrainte mécanique appliquée sur le pourtour de la lame de quartz est parla suite introduit dans le modèle en utilisant la méthode de perturbation de Sinha-Tiersten. Il est alors possible de définirl’impact des défauts de fabrication sur la fréquence du résonateur. Enfin, la méthode numérique est appliquée à l’étudede structures innovantes dans le cadre du projet FREQUENCE2009. Il s’agit de revisiter le concept du résonateur BVA etd’envisager des procédés de fabrication collective. L’idée consiste ainsi à remplacer le rayon de courbure d’un résonateur,dont la fréquence utile se trouve aux alentours de 9 MHz, par une série de marches, plus compatible avec les procédés dela micro électronique (DRIE : Deep Reactive Ion Etching). Bien que les résultats expérimentaux soient, dans ce cas, loin denos attentes, nous constatons que l’outil d’analyse est parfaitement capable de prédire les caractéristiques de nouvellesstructures. / This work is devoted to the development of a digital analysis tool dedicated to study new bulk acoustic waves quartz resonatorsby using finite elements. This method of characterization allows the calculation of the elements of the equivalentelectrical circuit (resistor, inductance and capacitor) of a given resonant frequency as well as the quality factor, while takinginto account its sensitivity to the temperature and to the stresses induced by the mounting support. Firstly, a validationphase is carried out in order to check our choices, in terms of modeling and computation, by comparing simulation data tothe measures of existing resonators. The three analyzed devices (40MHz, 10 MHz and 100 MHz) show good agreementbetween theory and experiment. To obtain such results, the mounting support is taken into account and modeled thanks toRayleigh damping areas when the trapping of energy is not optimal (presence of a spurious shell vibration mode). Then, animportant aspect of resonators is studied : the temperature behavior of its vibrating modes. Indeed, the thermal expansionstresses as well as the change of stiffness coefficients according to the temperature induce frequency shift. The comparisonbetween theory and experiment allows us to check the shape of curves and to quantify the accuracy of the model.Thereafter, the effect of mechanical stress applied on the edge of the blank of quartz is introduced in the model by usingthe perturbation method developed by Tiersten and Sinha. So, it is possible to define the influence of some manufacturingdefects on the resonant frequency. Finally, the digital method is applied to study innovative structures in the framework ofthe project FREQUENCE2009. The aim is to review the concept of BVA resonator and consider collective manufacturingprocesses. The idea involves replacing the radius of curvature of a resonator, for which the expected frequency is around9 MHz, by several steps, more compatible with microelectronics processes (DRIE : Deep Reactive Ion Etching). Althoughthe results are far from our expectations, we note that the analysis tool is perfectly able to anticipate the characteristics ofnew structures.

La formation du continent sud américain résulte d'une histoire complexe étalée sur plus de 3.5 Ga. Les processus qui ont façonné le continent ont impliqué l'intégralité de la lithosphère. La structure du manteau supérieur a été étudiée à l'échelle du continent en réalisant un modèle tomographique en ondes de surface, tandis qu'une étude du déphasage des ondes de cisaillement a permis d'appréhender la structure du manteau supérieur sous le sud-est du Brésil.<br /> La tomographie sismique anisotrope en ondes de surface donne une image de la structure en vitesse du manteau supérieur sous le continent et les océans environnants. Une bonne corrélation existe entre les grandes structures géologiques et les hétérogénéités de vitesse. L'anisotropie des ondes de Rayleigh suggère une absence de déformation à grande échelle au-delà de 200 km de profondeur, ce qui contraint verticalement la source de l'anisotropie mise en évidence au sud-est du Brésil. L'étude du déphasage des ondes de cisaillement y a révélé une orientation dominante du plan de polarisation de l'onde quasi-S rapide, parallèle aux chaînes péricratoniques. Un déphasage entre les ondes S rapide et lente supérieur à 2s a été mesuré à l'aplomb des décrochements majeurs. Une telle amplitude suggère que : 1) les décrochements traversent la lithosphère, 2) la lithosphère est caractérisée par une très forte anisotropie intrinsèque, 3) l'asthénosphère contribue au déphasage, par la présence de deux couches ou d'une fabrique tectonique cohérente entre lithosphère et asthénosphère, impliquant une absence de découplage depuis le Néoprotérozoïque.<br /> Une modélisation numérique tridimensionnelle de la déformation lithosphérique a permis d'étudier le couplage mécanique entre la croûte et le manteau supérieur, ainsi que la localisation de la déformation en termes de développement de zones de cisaillement.

Ce travail vise à améliorer la compréhension des signaux sismiques dérivés des fonctions de corrélation inter-récepteur du bruit sismique, ce qui est critique pour une imagerie fiable de la Terre profonde basée sur le bruit. La thèse comprend sept chapitres. Le chapitre 1 introduit les connaissances de base sur le bruit sismique, de la terminologie à ses origines diverses. Le chapitre 2 fournit une vue d'ensemble de la littérature sur l'historique et le développement de la méthode récente de corrélation de bruit, et passe en revue diverses techniques pour le prétraitement des données de bruit sismique et le post-traitement des fonctions de corrélation de bruit. Des méthodes de traitement du bruit basées sur les statistiques et un schéma modifié pour calculer la fonction de corrélation sont développés dans ce chapitre. Le chapitre 3 propose plusieurs techniques basées sur la transformée de Radon pour mesurer les lenteurs des champs d'ondes corrélés et analyser en termes de phases sismiques les signaux dérivés du bruit. Le chapitre 4 montre que les ondes de volume sondant la Terre profonde peuvent être extraites des corrélations de bruit à des distances télésismiques, avec des enregistrements de bruit provenant de deux réseaux sismiques régionaux. Le chapitre 5 applique les techniques proposées au chapitre 3 aux corrélations de bruit entre deux réseaux calculées au chapitre 4, et permet de comprendre l’origine de la phase précoce non-physique observée dans les données. Le chapitre 6 discute des conditions dans lesquelles apparaissent des phases sans correspondance dans la réponse physique de la Terre qui peuvent fausser les analyses des structures profondes basées sur le bruit.. Le dernier chapitre fournit un résumé sur les contributions de cette thèse et une perspective de plusieurs travaux soit en cours soit envisagés pour le futur. / This work aims toward an improved understanding of the seismic signals derived from the inter-receiver correlation functions of seismic noise, which is valuable and critical for a reliable noise-based deep Earth imaging. The thesis consists of seven chapters. Chapter 1 introduces background knowledge on seismic noise, from its classifications to various origins. Chapter 2 provides a literature overview on the history and development of the emerging noise correlation method, and reviews various techniques for the pre-processing of seismic noise data and post-processing of noise correlation functions. Statistics-based noise processing methods and a modified scheme for computing correlation function are developed in this chapter. Chapter 3 proposes several Radon-based techniques to analyze the slownesses of correlated wavefields and to unveil the origin of noise-derived seismic signals. Chapter 4 shows that body waves penetrating into deep Earth can be extracted from noise correlations at teleseismic distances, with noise records from two regional seismic networks. Chapter 5 applies the techniques proposed in chapter 3 to the double-array nose correlations computed in chapter 4, and accordingly reveals the origin of an early spurious phase observed in chapter 4. Chapter 6 discusses several situations that bring ambiguities into the noise-derived seismic signals and can potentially bias the noise-based imaging of subsurface structure. The last chapter provides a summarization over the contributions of this thesis and an outlook of several ongoing and prospected works.

Cette thèse a pour objectif de développer des méthodes ultrasonores pour l’END de l’adhésion de joints collés. Pour aborder ce problème, les assemblages sont réalisés avec des substrats en aluminium (isotrope, élastique) et un adhésif de type époxy(isotrope, viscoélastique). Selon la géométrie des assemblages, deux méthodes sont proposées pour obtenir une information quantitative sur le niveau d’adhésion. La première est adaptée à des échantillons de type tri-couches avec recouvrement total.Elle consiste à analyser la transmission d’ondes ultrasonores planes de volume à travers l’assemblage immergé dans l'eau.Les conséquences d'un mauvais traitement de surface des substrats sur la mesure des modules de viscoélasticité du joint adhésif sont étudiées. Il a été montré que des interphases dégradées provoquent une anisotropie apparente des modules mesurés pour le joint de colle. Cette anisotropie a été quantifiée à l'aide de deux paramètres β1 et β2 dont les valeurs permettent de révéler la qualité des interphases. Ensuite, les modules élastiques (ou raideurs kL et kT) des interphases ont été estimés, en supposant les propriétés de l'adhésif connues. On montre que leurs valeurs sont maximales lorsque l'adhésion est nominale, et qu’elles diminuent franchement lorsque l'adhésion est dégradée, mais dans des proportions différentes. Des mesures de la contrainte à rupture, réalisées sur des échantillons préparés dans les mêmes conditions, corroborent la chute des modules des interphases. La seconde méthode est adaptée à des échantillons collés avec un recouvrement partiel. Elle es tbasée sur la mesure du coefficient de transmission d'ondes de Lamb se propageant d'un substrat à l'autre, à travers la zone de recouvrement. Une étude de sensibilité numérique (par éléments finis) des coefficients de transmission des modes de Lamb a montré que les propriétés mécaniques des interphases (modélisées par des raideurs surfaciques) peuvent être évaluées si les autres caractéristiques de l’assemblage sont connues. Des mesures expérimentales de ces coefficients de transmission ont ensuite été réalisées avec deux échantillons. L’un d’eux possède des interphases à adhésion nominale et l’autre des interphases dégradées. Une confrontation entre les résultats des mesures obtenus pour les différents modes et les simulations numériques permet de déterminer les valeurs des raideurs d’interfaces pour chaque échantillon. Là encore, il est observé qu’une mauvaise adhésion se traduit par des valeurs faibles des raideurs d’interfaces, qui peuvent être quantifiées, cette fois,grâce aux ondes ultrasonores guidées. / The aim of this thesis is to develop NDT ultrasonic methods for the evaluation of the adhesion at interfaces in bonded joints.To address this problem, the assemblies are made of aluminum substrates (isotropic, elastic) and epoxy adhesive (isotropic,viscoelastic). Depending on the geometry of the studied assemblies, two methods are proposed to obtain quantitativeinformation on the adhesion level. The first method is suitable for three-layered plate-like samples. It consists in analyzingthe transmission of bulk ultrasonic plane waves through the assembly immersed in water. The consequences of a inapropriatesurface treatment of the substrates on the measurement of viscoelastic moduli of the joint are studied. It has been shown thatdegraded interphases cause an apparent anisotropy of the measured elastic moduli of the joint. This anisotropy was quantifiedusing two parameters β1 and β2 whose values can reveal the quality of the interphases. Then, the elastic moduli (or stiffnesseskL and kT) of the interphases were estimated, assuming that the adhesive layer properties are known. It was shown that theirvalues are higher when the adhesion is nominal, and are strongly decreased when the adhesion is degraded. Measurements ofthe mechanical strength, on samples prepared under the same conditions, confirmed the drops in these elastic moduli. Thesecond method is more suitable for lap joint samples. It is based on the Lamb wave transmission coefficient measured fromone substrate to the other, across the overlap zone. A numerical sensitivity analysis (finite element model) of the transmissioncoefficients of Lamb waves showed that the mechanical properties of the interphases (modeled by interfacial stiffnesses) canbe evaluated if the other characteristics of the assembly are well known. Experimental measurements of the transmissioncoefficients were then performed on two samples. One of them has interphases with nominal adhesion and the other hasdegraded interphases. A comparison between the measured results obtained for the different modes and numericalsimulations was used to determine the values of the interfacial stiffnesses for each sample. Again, it was observed that pooradhesion leads to low values of the interfacial stiffnesses of the interphases, that can be quantified using guided ultrasonicwaves.

La connaissance des structures internes de la Terre, à différentes échelles, présente des enjeux majeurs d'ordres économiques, humains, environnementaux et scientifiques. Diverses méthodes d'imagerie ont été développées en utilisant les informations contenues dans les ondes sismiques. La méthode d'inversion des formes d'ondes construit des images quantitative haute résolution des paramètres physiques du sous-sol, en exploitant le champ d'onde complet, sous la forme d'un problème d'optimisation. Dans ce travail de thèse, je présente l'application de l'inversion des formes d'ondes en domaine fréquentiel, pour imager les paramètres visco-élastiques dans des géometries à deux dimensions à grands offsets. Dans un premier temps les développements méthodologiques et algorithmiques sont présentés. La modélisation de la propagation des ondes P-SV en domaine fréquentiel, le problème direct du processus d'imagerie, est assurée par une méthode d'éléments finis Galerkin discontinus, assurant une grande flexibilité dans le choix des ordres d'interpolation et dans l'utilisation de maillages triangulaires non-structurés. Le problème inverse est résolu sous une forme linéarisée, afin de limiter le nombre de simulations directes, et utilise l'algorithme quasi-Newton L-BFGS permettant de tirer bénéfice de l'estimation "économique" du Hessien. Le processus global d'imagerie est implémenté sous la forme d'un algorithme massivement parallèle destiné aux calculateurs modernes à mémoire distribuée. Dans un deuxième temps, les algorithmes développés sont appliqués à des cas d'étude. Des applications sont menées dans des modèles synthétiques réalistes représentatifs d'environnements terrestres et marins. Ces études montrent les difficultés associées à la reconstruction des paramètres élastiques à partir de données mettant en jeu des phénomènes de propagations complexes (ondes converties, multiples, ondes de surfaces...). Des solutions sont proposées sous forme de processus hiérarchiques multi-échelles, afin de limiter les effets des non-linéarités du problème inverse et ainsi d'améliorer la convergence du processus vers le minimum global. Enfin, la sensibilité de différentes normes et critères de minimisation est analysée, à partir de données bruités issues de modèles synthétiques réalistes, ainsi que sous l'approximation acoustique pour un jeu de données réelles pétrolière. Ces tests montrent certaines limites du formalisme classique basé sur la norme L2 dans l'espace des données, tandis que la norme L1 apparaît comme alternative robuste pour l'inversion de données décimées en domaine fréquentiel.

Cette thèse présente des résultats de recherche sur des résonateurs acoustiques à ondes de volume travaillant à températures cryogéniques (de 3K à 15K dans un réfrigérateur à tube pulsé) ainsi que sur les systèmes construits autour de ces composants. Le premier aspect concerne le fonctionnement de différents composants, dont en particulier le résonateur. Le comportement de ce dernier a fait l'objet d'une étude plus systématiquement sur la gamme de fréquence [1−90 MHz]. Des facteurs de qualité de 417 * 106 et des produit Q* f (facteur de qualité Q à la fréquence f) pouvant atteindre 3, 07*1016 Hz ont été mesurées pour les résonateurs à quartz, valeurs exceptionnelles qui sont des records mondiaux pour cette classe de composants. Il est montré que le facteur de qualité Q ne dépend pas de la fréquence f à ces températures, conformément à la théorie de Landau-Rumer. Les problèmes et avantages de travailler à de telles températures sont évalués. Les limitations qui s'y rattachent sont discutées. D'autres composants passifs ou actifs, tels que des transistors, sont aussi d'étudiés. Le choix des composants appropriés est fait sur la base de comparaisons de leur comportement à 4K. Les résultats sont confirmés par une modélisation d'amplificateurs cryogéniques, réalisée avec succès. Le deuxième aspect mis en lumière dans ce travail est la modélisation et la simulation des composants et des systèmes étudiés. Le rapport présente un modèle rigoureux du bruit de phase des composants à onde acoustique de volume prenant en compte les non linéarités et basé sur la méthode de la moyenne. Ce modèle du résonateur, de type MIMO (multiple-input and multiple-output) est utilisé pour en déduire l'impact des différents paramètres d'un oscillateur sur son bruit de phase et expliquer les résultats expérimentaux. Les modèles établis pour les différents autres composants sont aussi utilisés pour simuler et optimiser des amplificateurs et oscillateurs cryogéniques. Le troisième aspect concerne des sources de fréquence basées sur ces résonateurs à ondes de volume cryogéniques : référence passive et oscillateur cryogéniques. Les systèmes réalisés ont permis de mesurer et caractériser le bruit du résonateur dans des conditions spécifiques. Certaines sources de bruit ont ainsi été identifiées. Les systèmes à boucle d'asservissement testés ont une stabilité relative de fréquence de 4 * 10−13 à 100 secondes et restent à mieux que 1 * 10−12 entre 1 et 2000 s. L'oscillateur réalisé a une stabilité de 1.5 * 10−12 à 200 s et meilleure que 1 * 10−11 pour des temps d'intégration plus grands que 80 ms. Les limitations de ces systèmes sont discutées sur la base des données obtenues.

La couche D" est une des zones les plus hétérogène de la Terre, abritant des structures mal comprises, probablement à l'origine de la dynamique globale du manteau. Les méthodes actuelles de tomographie sismique ne résolvent pas correctement cette région. La résolution des images de la structure interne de la Terre est fortement liée à la localisation précise des séismes, ainsi qu'à la connaissance du contenu fréquentiel des ondes sismiques. Cette thèse a pour objectif la construction d'une base de données globale de temps de trajet des ondes, mesures utilisées pour l'élaboration des images tomographiques. Notre méthode utilise une détermination plus précise de la profondeur des séismes permettant d'améliorer la mesure des temps de trajet. Pour une meilleure évaluation du contenu fréquentiel, ces mesures sont réalisées dans plusieurs gammes de fréquence. Une nouvelle approche de calcul de corrections visant à éliminer la complexité de la croûte est également présentée. Enfin, une imagerie préliminaire est présentée aux échelles régionale et globale sous l'approximation de la théorie des rais, avant une description plus complète par l'approximation de Born.

Nous avons caractérisé différents cristaux piézoélectriques de Langatate La3Ga5.5Ta0.5O14 (LGT), de différentes couleurs afin d’évaluer l’influence des défauts sur la qualité de dispositifs du domaine Temps-Fréquence (résonateurs acoustiques).Nous avons étudié la composition de ces cristaux dont les variations autour de la composition stoechiométrique sont surtout dues à l’évaporation du Ga2O au cours de la croissance compensée par un excès de Ga2O3 dans le mélange initial. Pour déterminer cette composition, la technique ICP-AES, précédée par une mise en solution par fusion alcaline, est la plus fiable.La présence d’impuretés chimiques, telles que les éléments métalliques ou les terres rares, liée à la pureté des oxydes de base, au frittage de la charge dans un creuset alumine… est étudiée. Nous expliquons ainsi la nature des centres colorés qui évolue selon l’atmosphère de tirage ou lors d’un recuit particulier.C’est ainsi que nous avons déterminé certaines propriétés de LGT permettant l’obtention de dispositifs comme les résonateurs à ondes de volume dont le produit Q.f est supérieur à celui des résonateurs à quartz. Pour cela, nous montrons notamment qu’il est nécessaire que :- la composition du cristal soit la plus proche possible de la stoechiométrie,- la résistivité électrique soit la plus grande possible,- l’étude spectroscopique ne révèle aucune absorption dans le domaine du visible.Dans ces conditions, le facteur de qualité Q d’une résonance à 10 MHz peut être de 1.44 million au point d’inversion de la courbe fréquence-température (1.35 pour le quartz) mais qu’il diminue notablement pour atteindre 0.35 million si nous pratiquons un recuit sous air à 1000°C pendant 48 h. / We have characterized Langatate piezoelectric crystals La3Ga5.5Ta0.5O14 (LGT), differently colored in order to evaluate the defects influence on the quality of Bulk Acoustic Waves resonators for the Time and Frequency domain.We have analyzed the composition of crystals whose variations around stoichiometric composition are mainly due to the evaporation of the Ga2O during growth compensated by adding an excess of Ga2O3 in the initial mixture. Among different techniques, the ICP-AES spectrometry, preceded by a dissolution by alkaline fusion seems to be the most accurate technique to determine the composition.The inevitable presence of chemical impurities, such as metallic elements, rare earth… linked to the purities of the raw materials, to the sintering of the load in an alumina crucible… is also studied. We try to explain in particular the nature of the color centers that evolve according to the growth atmosphere or during a particular annealing.So, we conclude by a list of necessary properties to obtain BAW resonators exhibiting a Q.f product higher than quartz. For this, we establish that it is necessary that:- the crystal composition is as close as possible to the stoichiometric composition,- the electrical resistivity is the highest possible,- the spectroscopic study does not reveal any absorption band in the visible domain.In these conditions, we were able to highlight that the quality factor of a 10 MHz resonance is of 1.44 million at the inversion temperature of the frequency-temperature curve (1.35 for the quartz crystal in the same conditions), but that this one decreases significantly to reach 0.35 million if we perform an annealing under air at 1000 °C during 48 hours.