Propagation des ondes sismiques dans les milieux biphasés PDF Download
Are you looking for read ebook online? Search for your book and save it on your Kindle device, PC, phones or tablets. Download Propagation des ondes sismiques dans les milieux biphasés PDF full book. Access full book title Propagation des ondes sismiques dans les milieux biphasés by Guy Kuster. Download full books in PDF and EPUB format.
Author: Guy Kuster
Publisher:
ISBN:
Category :
Languages : fr
Pages : 232
Get Book Here
Book Description
Nous présentons des résultats théoriques et expérimentaux sur la réponse macroscopique des milieux biphasés à des ondes élastiques. Dans le modèle théorique de base on suppose que la matrice est solide, que les inclusions sont sphériques et très petites devant les longueurs d'onde et que les interactions entre inclusions sont négligeables. L'examen de ces quatre hypothèses nous conduit à proposer quatre modèles nouveaux: Quand la matrice est un fluide les résultats statiques ne peuvent être appliqués au cas dynamique en raison des effets inertiels. Quand les inclusions sont ellipsoïdales les propriétés effectives dépendent à la fois de la forme et de la concentration des inclusions. En moyenne les interactions entre inclusions sphériques affectent uniquement le module de rigidité effectif. Quand la longueur d'onde n'est plus très grande devant la taille des inclusions, les propriétés effectives dépendent de la fréquence. Pour une suspension de sphères solides dans un liquide les effets dispersifs peuvent être importants. Les études expérimentales portent sur la vitesse et l'atténuation des ondes ultrasoniques dans trois suspensions de sphères solides dans un liquide en fonction de la concentration. Dans tous les cas les vitesses mesurées sont expliquées au mieux par notre modèle théorique où les interactions et les effets dispersifs dus à l'utilisation d'une longueur d'onde finie sont pris en considération. Les atténuations observées se comparent de manière satisfaisante aux valeurs théoriques calculées en tenant compte des pertes par diffraction, de l'absorption due à la viscosité de la matrice et de l'anélasticité intrinsèque des inclusions sphériques
Author: Guy Kuster
Publisher:
ISBN:
Category :
Languages : fr
Pages : 232
Get Book Here
Book Description
Nous présentons des résultats théoriques et expérimentaux sur la réponse macroscopique des milieux biphasés à des ondes élastiques. Dans le modèle théorique de base on suppose que la matrice est solide, que les inclusions sont sphériques et très petites devant les longueurs d'onde et que les interactions entre inclusions sont négligeables. L'examen de ces quatre hypothèses nous conduit à proposer quatre modèles nouveaux: Quand la matrice est un fluide les résultats statiques ne peuvent être appliqués au cas dynamique en raison des effets inertiels. Quand les inclusions sont ellipsoïdales les propriétés effectives dépendent à la fois de la forme et de la concentration des inclusions. En moyenne les interactions entre inclusions sphériques affectent uniquement le module de rigidité effectif. Quand la longueur d'onde n'est plus très grande devant la taille des inclusions, les propriétés effectives dépendent de la fréquence. Pour une suspension de sphères solides dans un liquide les effets dispersifs peuvent être importants. Les études expérimentales portent sur la vitesse et l'atténuation des ondes ultrasoniques dans trois suspensions de sphères solides dans un liquide en fonction de la concentration. Dans tous les cas les vitesses mesurées sont expliquées au mieux par notre modèle théorique où les interactions et les effets dispersifs dus à l'utilisation d'une longueur d'onde finie sont pris en considération. Les atténuations observées se comparent de manière satisfaisante aux valeurs théoriques calculées en tenant compte des pertes par diffraction, de l'absorption due à la viscosité de la matrice et de l'anélasticité intrinsèque des inclusions sphériques
Author: Bastien Dupuy
Publisher:
ISBN:
Category :
Languages : fr
Pages : 0
Get Book Here
Book Description
La propagation des ondes sismiques dans les milieux poreux multiphasiques présente des enjeux nombreux, tant sur le plan environnemental (risques naturels, géotechnique, pollutions de nappes...) que pour les réservoirs (aquifères, hydrocarbures, stockages de CO2...). L'utilisation des ondes sismiques pour étudier ces milieux se justifie par le fait qu'en se propageant, les ondes sont déformées par le milieu qu'elles traversent et contiennent ainsi des informations aux capteurs sur les phases fluides et solides et sur le squelette poreux. Ce travail de thèse s'intéresse aux caractéristiques des ondes sismiques dans les milieux multiphasiques (plusieurs phases fluides et solides), depuis la description physique jusqu'à la caractérisation des paramètres constitutifs par inversion, en passant par la modélisation numérique 2D de la propagation. La première partie du travail a consisté à décrire la physique des milieux multiphasiques (phase par phase et leurs intéractions dynamiques) en utilisant des méthodes d'homogénéisation pour se ramener à un milieu équivalent défini par sept paramètres. Ainsi, dans des milieux simple porosité saturés et dans des milieux plus complexes (double porosité, partiellement saturés ou visco-poroélastiques), je peux calculer la propagation des ondes sismiques sans approximation. En effet, j'utilise une méthode numérique dans le domaine fréquence-espace qui permet de prendre en compte tous les termes qui dépendent de la fréquence sans approximation. La discrétisation spatiale utilise une méthode d'éléments finis discontinus (Galerkin discontinu) qui permet de considérer des milieux hétérogènes.Je montre notamment que les attributs sismiques (vitesses et atténuations) des milieux poreux complexes sont fortement dispersifs et les formes d'ondes complètes, calculées sans approximation, sont fortement dépendantes de la description physique du milieu. La caractérisation des paramètres poroélastiques s'effectue par inversion. Une méthode en deux étapes a été proposée : la première consiste en une inversion ``classique`` (tomographie, inversion des formes d'ondes complètes) des données (sismogrammes) pour obtenir des paramètres macro-échelles (attributs sismiques). La seconde étape permet de reconstruire, à partir des paramètres macro-échelles, les paramètres poroélastiques micro-échelles. Cette étape d'inversion utilise une méthode d'optimisation semi-globale (algorithme de voisinage). Une analyse de sensibilité montre qu'en connaissant a-priori certains paramètres, on peut inverser avec précision les paramètres du squelette poroélastique ou retrouver la nature du fluide saturant, à partir des vitesses de propagation. En revanche, pour retrouver la saturation en fluide, il est préférable de connaître les atténuations. Deux applications réalistes (monitoring de réservoir et hydrogéophysique) mettent en oeuvre ce type d'inversion en deux étapes et démontrent qu'à partir de données estimées par des méthodes classiques d'imagerie, on peut remonter à certains paramètres poroélastiques constitutifs.
Author: Domitille Anache-Ménier
Publisher:
ISBN:
Category :
Languages : fr
Pages : 0
Get Book Here
Book Description
Cette thèse est consacrée à l'étude théorique et expérimentale de la phase des ondes sismiques et des ultrasons se propageant dans les milieux désordonnés.La théorie des distributions et des corrélations des dérivées spatiales et temporelles de la phase est développée dans l'hypothèse d'un champ scalaire analytique gaussien et circulaire. Ces fonctions statistiques permettent de caractériser les diffuseurs dans les deux dispositifs expérimentaux au coeur de cette thèse. D'une part, les fluctuations temporelles de la phase d'ultrasons sont utilisées pour sonder la dynamique d'une suspension de billes millimétriques sur des échelles de temps allant de la milliseconde à la seconde. D'autre part les fluctuations spatiales de la phase donnent une caractérisation de la diffusion multiple des ondes de flexion dans une plaque de Plexiglas ® perforée aléatoirement. Le comportement asymptotique en loi de puissance des distributions des dérivées de la phase démontre les propriétés gaussiennes des codas dans ces deux dispositifs. Enfin, l'étude de la coda de séismes régionaux en Californie a permis de proposer une application à la détermination du libre parcours moyen des ondes sismiques dans la croûte terrestre : il est montré que c'est la seule échelle caractéristique de la fonction de corrélation de la dérivée spatiale de la phase.
Author: Fateh Bouchaala
Publisher:
ISBN:
Category :
Languages : fr
Pages :
Get Book Here
Book Description
La modélisation de la propagation des ondes dans les milieux complexes 3d, inhomogènes, isotropes et viscoélastiques (milieux I.I.V.), est un sujet intéressant et largement d’actualité. cette modélisation aide pour la connaître la géométrie et les propriétés physiques des milieux géologiques. La viscosité du milieu de propagation est modélisée par la prise en compte de vitesses de propagation complexes. Le phénomène de propagation est modélisé par la théorie des rais avec la prise en compte de l’hypothèse de faible viscosité (hypothèse non restrictive dans le domaine de la géophysique). Lles expressions analytiques des coefficients de réflexion transmission entre deux milieux viscoélastiques ont été établies avec une étude numérique de l’effet de la faible viscosité sur ces coefficients a été menée. Le code anray est, a l’origine, un code de simulation de la propagation des ondes dans les milieux inhomogènes 3d purement élastiques. son développement repose sur la théorie des rais. Le code anray a ensuite été modifie (anrayq), afin de l’adapter a la simulation de la propagation des ondes dans les milieux I.I.V. les résultats numériques du code anrayq, ont été confrontés aux résultats numériques issus du code 5kb, code de simulation base sur la méthode de réflectivité. Enfin, des données sismiques 3d, acquises lors de la campagne hydratech (Ifremer 2002), dans une région caractérisée par la présence des hydrates de gaz et du gaz libre, ont été simulées a l’aide du code anrayq. le profil de l’atténuation viscoélastique a été aussi déterminé a partir de ces données sismiques, afin de localiser les hydrates de gaz et le gaz libre dans les couches sédimentaires.
Author: Domitille Anache-Ménier
Publisher:
ISBN:
Category :
Languages : fr
Pages : 172
Get Book Here
Book Description
Cette thèse est consacrée à l'étude théorique et expérimentale de la phase des ondes sismiques et des ultrasons se propageant dans les milieux désordonnés.La théorie des distributions et des corrélations des dérivées spatiales et temporelles de la phase est développée dans l'hypothèse d'un champ scalaire analytique gaussien et circulaire. Ces fonctions statistiques permettent de caractériser les diffuseurs dans les deux dispositifs expérimentaux au coeur de cette thèse. D'une part, les fluctuations temporelles de la phase d'ultrasons sont utilisées pour sonder la dynamique d'une suspension de billes millimétriques sur des échelles de temps allant de la milliseconde à la seconde. D'autre part les fluctuations spatiales de la phase donnent une caractérisation de la diffusion multiple des ondes de flexion dans une plaque de Plexiglas ® perforée aléatoirement. Le comportement asymptotique en loi de puissance des distributions des dérivées de la phase démontre les propriétés gaussiennes des codas dans ces deux dispositifs. Enfin, l'étude de la coda de séismes régionaux en Californie a permis de proposer une application à la détermination du libre parcours moyen des ondes sismiques dans la croûte terrestre : il est montré que c'est la seule échelle caractéristique de la fonction de corrélation de la dérivée spatiale de la phase.
Author: M. Lavergne
Publisher: Editions TECHNIP
ISBN: 9782710805144
Category : Seismic profiling
Languages : fr
Pages : 222
Get Book Here
Book Description
Author: Sylvie Kelner
Publisher:
ISBN:
Category :
Languages : fr
Pages : 0
Get Book Here
Book Description
La diffraction des ondes sismiques par des milieux fissures en deux dimensions (2d) est etudiee par une methode d'equations integrales aux frontieres ou les fonctions de green sont calculees par la methode des nombres d'ondes discrets (dwbiem : discrete wavenumber boundary integral equations method). Cette methode semi-analytique est particulierement bien adaptee aux problemes de la propagation des ondes sismiques dans un milieu homogene contenant des fissures vides ou remplies de fluide. Toutes les conversions d'ondes sont modelisees en appliquant la dwbiem. En premier lieu, nous avons etudie, par simulations numeriques, comment des milieux fissures pouvaient etre caracterises sismiquement. Nous avons ainsi pu observer des phenomenes d'attenuation et d'anisotropie, selon que la longueur d'onde du champ d'ondes incident est du meme ordre de grandeur ou qu'elle est plus grande que la longueur des fissures. Nous avons retrouve un resultat deja connu qui concerne l'attenuation des ondes elastiques lorsqu'elles traversent un milieu fissure : l'attenuation est maximale lorsque la longueur d'onde incidente est proche de la longueur des fissures. Par ailleurs, nous avons modelise la couche de granite fissuree du site de garner valley, en californie, en nous basant sur la theorie des milieux homogenes equivalents. Plusieurs modeles de milieux fissures restituent bien le taux d'anisotropie observe a garner valley. Une etude comparative de l'attenuation d'ondes s enregistrees la-bas et d'ondes s synthetiques permet de conclure que l'anisotropie s'explique par la presence de fissures verticales (et non horizontales) mais ne permet pas de privilegier un modele plus qu'un autre. Enfin, nous avons simule numeriquement la reponse hydro-mecanique d'un massif fracture a un seisme. Les deformations des fissures et les variations de pression dues au champ d'ondes rayonne par une faille en glissement permettent de connaitre les zones d'expulsion de fluide.
Author: Patricia Rodriguez Tome ((Docteur en physique, université de Paris 6, 1983).)
Publisher:
ISBN:
Category :
Languages : fr
Pages : 182
Get Book Here
Book Description
Author: SYLVIE.. KELNER
Publisher:
ISBN:
Category :
Languages : fr
Pages : 212
Get Book Here
Book Description
LA DIFFRACTION DES ONDES SISMIQUES PAR DES MILIEUX FISSURES EN DEUX DIMENSIONS (2D) EST ETUDIEE PAR UNE METHODE D'EQUATIONS INTEGRALES AUX FRONTIERES OU LES FONCTIONS DE GREEN SONT CALCULEES PAR LA METHODE DES NOMBRES D'ONDES DISCRETS (DWBIEM : DISCRETE WAVENUMBER BOUNDARY INTEGRAL EQUATIONS METHOD). CETTE METHODE SEMI-ANALYTIQUE EST PARTICULIEREMENT BIEN ADAPTEE AUX PROBLEMES DE LA PROPAGATION DES ONDES SISMIQUES DANS UN MILIEU HOMOGENE CONTENANT DES FISSURES VIDES OU REMPLIES DE FLUIDE. TOUTES LES CONVERSIONS D'ONDES SONT MODELISEES EN APPLIQUANT LA DWBIEM. EN PREMIER LIEU, NOUS AVONS ETUDIE, PAR SIMULATIONS NUMERIQUES, COMMENT DES MILIEUX FISSURES POUVAIENT ETRE CARACTERISES SISMIQUEMENT. NOUS AVONS AINSI PU OBSERVER DES PHENOMENES D'ATTENUATION ET D'ANISOTROPIE, SELON QUE LA LONGUEUR D'ONDE DU CHAMP D'ONDES INCIDENT EST DU MEME ORDRE DE GRANDEUR OU QU'ELLE EST PLUS GRANDE QUE LA LONGUEUR DES FISSURES. NOUS AVONS RETROUVE UN RESULTAT DEJA CONNU QUI CONCERNE L'ATTENUATION DES ONDES ELASTIQUES LORSQU'ELLES TRAVERSENT UN MILIEU FISSURE : L'ATTENUATION EST MAXIMALE LORSQUE LA LONGUEUR D'ONDE INCIDENTE EST PROCHE DE LA LONGUEUR DES FISSURES. PAR AILLEURS, NOUS AVONS MODELISE LA COUCHE DE GRANITE FISSUREE DU SITE DE GARNER VALLEY, EN CALIFORNIE, EN NOUS BASANT SUR LA THEORIE DES MILIEUX HOMOGENES EQUIVALENTS. PLUSIEURS MODELES DE MILIEUX FISSURES RESTITUENT BIEN LE TAUX D'ANISOTROPIE OBSERVE A GARNER VALLEY. UNE ETUDE COMPARATIVE DE L'ATTENUATION D'ONDES S ENREGISTREES LA-BAS ET D'ONDES S SYNTHETIQUES PERMET DE CONCLURE QUE L'ANISOTROPIE S'EXPLIQUE PAR LA PRESENCE DE FISSURES VERTICALES (ET NON HORIZONTALES) MAIS NE PERMET PAS DE PRIVILEGIER UN MODELE PLUS QU'UN AUTRE. ENFIN, NOUS AVONS SIMULE NUMERIQUEMENT LA REPONSE HYDRO-MECANIQUE D'UN MASSIF FRACTURE A UN SEISME. LES DEFORMATIONS DES FISSURES ET LES VARIATIONS DE PRESSION DUES AU CHAMP D'ONDES RAYONNE PAR UNE FAILLE EN GLISSEMENT PERMETTENT DE CONNAITRE LES ZONES D'EXPULSION DE FLUIDE.
Author: Ronan Madec
Publisher:
ISBN:
Category :
Languages : en
Pages :
Get Book Here
Book Description
Dans l'industrie pétrolière, de nombreuses campagnes sismiques visant à trouver de nouveaux champs pétroliers sont effectuées en mer. Afin de reproduire numériquement les résultats de ces campagnes d'acquisitions sismiques marine, les méthodes de différences finies ou d'éléments finis sont le plus souvent utilisées. Dans cette thèse, nous travaillons avec la méthode des éléments spectraux (Komatitsch et Tromp, 2002), qui est une méthode d'éléments finis explicite à haut degré d'interpolation polynomiale. Pour de l'acquisition sismique marine, la propagation d'ondes s'effectue dans des modèles composés de deux parties : une partie fluide (la couche d'eau épaisse homogène représentant l'océan) et une partie solide (les roches composant le fond des océans). La partie fluide étant moins intéressante que la partie solide où peuvent se trouver les réservoirs d'hydrocarbures, nous voulons une méthode numérique où la propagation d'ondes dans le fluide est la moins coûteuse possible en temps calcul. Un moyen de gagner ce temps est d'utiliser une méthode intégrant des pas de temps locaux (Diaz et Joly, 2005) où la partie fluide est traitée à l'aide d'un pas de temps plus grand que dans la partie solide. Pour cela, il faut s'assurer que l'augmentation du pas de temps sans modifier le maillage ne change pas la condition de stabilité CFL dans le fluide. Cette condition est calculée à partir de la vitesse maximale de propagation des ondes dans le domaine, et des pas de temps et d'espace. Or, les modèles utilisés présentent généralement un fort contraste des vitesses à l'interface fluide-solide en faveur du solide. C'est donc le pas de temps de ce dernier qui gouverne la condition CFL. Donc augmenter le pas de temps seulement dans le fluide ne changera pas la condition CFL et la simulation numérique pourra s'effectuer sans risquer l'explosion de la solution pour cause d'instabilité du schéma numérique. Mon travail a consisté à intégrer un tel procédé dans la méthode des éléments spectraux du code SPECFEM à 2D. Comme on a des pas de temps différents de part et d'autre de l'interface fluide-solide, les données à échanger entre les deux milieux ne sont pas coordonnées car à des temps différents. Pour palier à cette difficulté, on impose la conservation de l'énergie à notre système pour la construction des conditions de raccord. Cette manoeuvre implique néanmoins la résolution d'un système linéaire l'interface mais son coût est négligeable comparé aux calculs économisés dans la partie fluide. Afin de se rapprocher plus encore de la réalité des modèles pétroliers, j'ai travaillé sur l'implémentation de couches absorbantes de type CPML (convolutional perfectly match layer, voir Komatitsch et Martin 2007 ) qui permettent de simuler un milieu infini en absorbant toutes les ondes y pénétrant et donc évitent des réflexions parasites des bords du domaine. Ces conditions ont été implémentées pour des modèles fluide-solide à 2D.
