Éditeur(s) : HAL CCSD
Résumé : Today, one of the main challenges in the oil industry is the exploitation of new resources in naturally fractured reservoirs often located in structurally complex areas such as plate boundaries, mountain ranges or near salt structures. While numerical models of rock deformation based on continuum mechanics are becoming industry standard in providing efficient means for modelling natural fractures in reservoirs, the knowledge of past perturbed stress fields, which is referred as the type (normal, wrench or reverse), orientation and magnitude of the regional or local tectonic stresses that needs to be applied as boundary condition in numerical simulations, is often unknown or approximated.The essence of this thesis is to improve the recovery of the paleostress field needed to constrain numerical simulations. I first described and discuss through a comprehensive mechanically-based multi-parametric study how far we can rely on the commonly used paleostress inversion methods based on fault slip data and the Wallace and Bott assumptions. I then compared this method with a new generation technique based on geomechanics using field observations. Since slip markers on faults are hardly observed in core or image log, I used observed natural fracture data as main drivers for the inversion of the paleostress using geomechanically-based method. I demonstrated, through various outcrop and subsurface examples, how this can efficiently be done. Finally, I aimed at addressing the problem of polyphase fracture data sets with unknown mechanical type (joints, faults, stylolites ...) and expanded the mechanical stress inversion to the separation of tectonic phases.
Aujourd'hui, l'un des principaux défis dans l'industrie pétrolière est l'exploitation de nouvelles ressources dans les réservoirs naturellement fracturés souvent situés dans des zones structurellement complexes tels que les frontières de plaques, les chaînes de montagnes ou les structures de proches des dômes de sel. Alors que les modèles numériques de déformation des roches basée sur la mécanique des milieux continus fournissant des moyens efficaces pour la modélisation de fractures naturelles dans les réservoirs, sont devenus la norme dans l'industrie, les paléo-contraintes régionales ou locales appliquées comme condition aux limites dans les simulations numériques, définies par le régime (normal, décrochant ou inverse), l'orientation, et la magnitude des contraintes tectoniques, sont souvent inconnues ou approximées.L'essence de cette thèse est d'améliorer la prédiction du champ de paléo-contrainte nécessaire pour contraindre les simulations numériques. J’ai décrit dans un premier temps à travers une étude comparative multi-paramétrique complète basée sur la mécanique, dans quelle mesure nous pouvions compter sur les méthodes d'inversion des paléo-contraintes basées sur les hypothèses de Wallace et Bott et les données de glissement sur les failles. J’ai ensuite comparé ces méthodes avec une technique d’inversion de nouvelle génération basée sur la géomécanique en utilisant les observations de terrain. Puisque les marqueurs de glissement sur les failles sont rarement observés le long du forage, j’ai utilisé l’inversion géomécanique des contraintes pour les données de fractures naturelles observées. J’ai alors démontré, l’efficacité de cette technique à travers divers exemples d'affleurements et de subsurface. Enfin, j’ai proposé une solution afin de résoudre le problème des données de fractures polyphasées à type cinématique inconnu (joints, failles, stylolites ...) et étendu l’inversion de contrainte basée sur la géomécanique à la séparation des phases tectoniques.
https://hal.archives-ouvertes.fr/tel-01324110

tel-01324110
https://hal.archives-ouvertes.fr/tel-01324110
https://hal.archives-ouvertes.fr/tel-01324110/document
https://hal.archives-ouvertes.fr/tel-01324110/file/theselejri.pdf

Éditeur(s) : HAL CCSD
Résumé : Joint inversions are now commonly used in Earth Sciences. They have been developed to better understand Earth structure since they provide more constraints on the inverted parameters. We propose a new process to simultaneously invert several data sets in order to better image 3-D crustal and upper mantle structures. Our inversion uses 3 kind of data which present a good complementarity : (1) P receiver functions to provide Moho depth variations, (2) teleseismic delay times of P waves to retrieve velocity anomalies in the crust and upper mantle, (3) gravity anomalies to image density variations at lithospheric scale. In our scheme, receiver functions are first inverted. The resulting Moho depths are then interpolated and incorporated as a priori information into the joint inversion of teleseismic delay times and gravity anomalies process. In our new approach, we perform a model space search for Moho variations, P-velocity and density structure to find acceptable fit to the three data sets. In order to preferentially sample the good data fit region, we chose the Neighborhood algorithm of Sambridge to optimistically survey the model space. We model the delay times with a 3-D raytracing using velocity nodes evenly spaced and linked with density nodes via a linear relationship. This relationship can evolve through depth, simulating its pressure dependance. We test our scheme on synthetic examples, and apply the method in a continental rifting region where deformation processes are complex and badly known. This joint inversion thus provides more constraints on the inverted parameters, and allows us to better image crust - mantle interactions in the case of the Hangai dome (Central Mongolia). Indeed, in this area, we expect an asthenospheric upwelling to dynamically support the topography, but we image a density - velocity anomaly in the lower crust. The effects on topography and the origin (composition vs. thermal) of this anomaly are still under debate.
Les inversions conjointes sont maintenant communément utilisées en Sciences de la Terre. Elles ont été développées afin d'améliorer notre connaissance et compréhension de la structure interne de la Terre puisqu'elles apportent de plus en plus de contraintes sur les paramètres inversés. Dans cette thèse, nous proposons un nouveau processus d'inversion conjointe qui prend en compte trois paramètres différents et qui mène à l'obtention d'un modèle lithosphérique tridimensionnel. Notre méthode utilise trois types de données gravimétriques et sismologiques qui présentent une bonne complémentarité : (1) Les fonctions récepteur P afin d'obtenir les variations de profondeur du Moho, (2) les délais de temps d'arrivées P des téléséismes pour retrouver les anomalies de vitesse dans la croûte et le manteau supérieur, (3) les anomalies gravimétriques qui donnent accès aux variations de densité à l'échelle lithosphérique. Dans le schéma d'inversion proposé, nous inversons d'abord les fonctions récepteur. Les variations de profondeur du Moho qui en résultent sont alors interpolées puis incorporées comme information a priori dans le processus d'inversion conjointe gravimétrie - tomographie télésismique. Que ce soit pour les fonctions récepteur ou pour la partie conjointe, notre méthode d'inversion se base sur un algorithme de recherche afin de trouver la structure qui minimise l'écart aux trois types de données. Nous avons opté pour l'algorithme de voisinage qui présente l'avantage de concentrer sa recherche dans les régions qui minimisent l'écart aux données sans toutefois abandonner les autres régions afin de limiter les risques de convergence vers un minimum local. Nous utilisons un modèle 3-D constitué de noeuds de vitesse afin de modéliser les délais de temps. Ces noeuds de vitesse sont liés aux noeuds de densité grâce à une relation linéaire entre vitesse et densité. Cette relation peut évoluer avec la profondeur de manière à simuler les effets dus aux variations de pression. Nous avons réalisé plusieurs tests synthétiques afin d'évaluer le comportement de notre nouvelle méthode. Ils ont permis de montrer que les variations du Moho ainsi que les anomalies de vitesse et densité sont bien retrouvées. Ils ont également montré l'importance de la prise en compte de ces trois types de données dans un même schéma d'inversion, notamment en ce qui concerne la faculté de discerner les variations d'interfaces des anomalies volumiques. Nous avons ensuite appliqué la méthode dans une région continentale en phase de rifting où les processus de déformation sont complexes et mal connus. Notre nouvelle méthode a apporté davantage de contraintes sur les paramètres inversés et a permis d'obtenir un modèle correctement résolu de la croûte et du manteau supérieur sous le dôme de Hangai (Mongolie Centrale). Une anomalie asthénosphérique sous le dôme est souvent invoquée pour expliquer sa topographie élevée mais c'est une anomalie (en vitesse et densité) crustale (croûte inférieure) qui s'exprime dans nos modèles. Ses effets sur la topographie et son origine (compositionnelle et/ou thermique) sont encore débattus et nécessiteraient des études plus poussées (paramétriques, VP =VS, ...).
https://tel.archives-ouvertes.fr/tel-00783524

tel-00783524
https://tel.archives-ouvertes.fr/tel-00783524
https://tel.archives-ouvertes.fr/tel-00783524/document
https://tel.archives-ouvertes.fr/tel-00783524/file/Basuyau2011.pdf

Éditeur(s) : HAL CCSD
Résumé : Different applications of linear elasticity in structural geology are presented in this thesis through the development of three types of numerical computer codes. The first one uses forward modeling to study displacement and perturbed stress fields around complexly faulted regions. We show that incorporating inequality constraints, such as static Coulomb friction, enables one to explain the angle of initiation of jogs in extensional relays. Adding heterogeneous material properties and optimizations, such as parallelization on multicore architectures and complexity reduction, admits more complex models. The second type deals with inverse modeling, also called parameter estimation. Linear slip inversion on faults with complex geometry, as well as paleo-stress inversion using a geomechanical approach, are developed. The last type of numerical computer code is dedicated to restoration of complexly folded and faulted structures. It is shown that this technique enables one to check balanced cross-sections, and also to retrieve fault chronology. Finally, we show that this code allows one to smooth noisy 3D interpreted faulted and folded horizons using geomechanics.
Différentes applications de l'élasticité linéaire en géologie structurale sont présentées dans cette thèse à travers le développement de trois types de codes numériques. Le premier utilise la modélisation directe pour étudier les déplacements et champs de contraintes autour de zones faillées complexes. On montre que l'ajout de contraintes inégalitaires, telles que la friction de Coulomb, permet d'expliquer l'angle d'initiation des dominos dans les relais extensifs. L'ajout de matériaux hétérogènes et d'optimisations, telles la parallélisation sur processeurs multi-coeurs ainsi que la réduction de complexité des modèles, permettent l'étude de modèles beaucoup plus complexes. Le second type de code numérique utilise la modélisation inverse, aussi appelée estimation de paramètres. L'inversion linéaire de déplacements sur les failles ainsi que la détermination de paléo-contraintes utilisant une approche géomécanique sont développées. Le dernier type de code numérique concerne la restauration de structures complexes plissées et faillées. Il est notamment montré qu'une telle méthode permet de vérifier l'équilibre de coupes géologiques, ainsi que de retrouver la chronologie des failles. Finalement, nous montrons que ce même code permet de lisser des horizons 3D faillés, plissés et bruités en utilisant la géomécanique.
https://tel.archives-ouvertes.fr/tel-00537899

tel-00537899
https://tel.archives-ouvertes.fr/tel-00537899
https://tel.archives-ouvertes.fr/tel-00537899/document
https://tel.archives-ouvertes.fr/tel-00537899/file/thesemaerten.pdf

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : This paper presents an inversion method using the interseismic velocity field to determine effective rigidity of the lithosphere. The method is based on the minimization of a cost function defined as the quadratic measure of the difference between measured and modelled velocity fields on a discrete set of points. The continuous mapping of the rigidity is fulfilled with a limited set of parameters and the forward solution is achieved using a plane stress finite element code. The computation of the cost function gradient in the parameters' space allows one to iteratively find the best parameters set through a suitable optimization algorithm.;We first design a benchmark including an abrupt rigidity variation that cannot be described by a continuous function. For such a case, we show that increasing the number of parameters is a way to accurately describe sharp variations of the rigidity map. Then, we use a dense GPS velocity field over the southwestern United States to estimate the corresponding rigidity variations for different spatial resolutions of the parameters' grid. We analyse the conceptual and practical difficulties associated with our methodology. Finally, rigidity maps obtained by our inversion method in southwestern United States and particularly across the San Andreas Fault System are reviewed and compared to current plate rigidity estimates and geophysical data over this area.
ISSN: 0956-540X

hal-00669897
https://hal.archives-ouvertes.fr/hal-00669897
DOI : 10.1111/j.1365-246X.2011.05192.x

Éditeur(s) : Ieee-inst Electrical Electronics Engineers Inc
Résumé : The synthetic aperture radar (SAR) Doppler centroid has been used to estimate the scatter line-of-sight radar velocity. In weak to moderate ocean surface current environment, the SAR Doppler centroid is dominated by the directionality and strength of wave-induced ocean surface displacements. In this paper, we show how this sea state signature can be used to improve surface wind retrieval from SAR. Doppler shifts of C-band radar return signals from the ocean are thoroughly investigated by colocating wind measurements from the ASCAT scatterometer with Doppler centroid anomalies retrieved from Envisat ASAR. An empirical geophysical model function (CDOP) is derived, predicting Doppler shifts at both VV and HH polarization as function of wind speed, radar incidence angle, and wind direction with respect to radar look direction. This function is used into a Bayesian inversion scheme in combination with wind from a priori forecast model and the normalized radar cross section (NRCS). The benefit of Doppler for SAR wind retrieval is shown in complex meteorological situations such as atmospheric fronts or low pressure systems. Using in situ information, validation reveals that this method helps to improve the wind direction retrieval. Uncertainty of the calibration of Doppler shift from Envisat ASAR hampers the inversion scheme in cases where NRCS and model wind are accurate and in close agreement. The method is however very promising with respect of future SAR missions, in particular Sentinel-1, where the Doppler centroid anomaly will be more robustly retrieved.
Ieee Transactions On Geoscience And Remote Sensing (0196-2892) (Ieee-inst Electrical Electronics Engineers Inc), 2012-07 , Vol. 50 , N. 7 , P. 2901-2909

Droits : 2012 IEEE
http://archimer.ifremer.fr/doc/00088/19896/17554.pdf
DOI:10.1109/TGRS.2011.2174998
http://archimer.ifremer.fr/doc/00088/19896/

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Our knowledge and understanding of the 3-D lithospheric structure of the Himalayas and the Tibetan Plateau is still challenging although numerous geophysical studies have been performed in the region. The GOCE satellite mission has the ambitious goal of mapping Earth's gravity field with unprecedented precision (i.e. an accuracy of 1-2 mGal for a spatial resolution of 100 km) to observe the lithosphere and upper-mantle structure. Consequently, it gives new insights in the lithospheric structure beneath the Himalayas and the Tibetan Plateau. Indeed, the GOCE gravity data now allow us to develop a new strategy for joint gravimetry-seismology inversion. Combined with teleseismic data over a large region in a joint inversion scheme, they will lead to lithospheric velocity-density models constrained in two complementary ways. We apply this joint inversion scheme to the Hi-CLIMB (Himalayan-Tibetan Continental Lithosphere during Mountain Building) seismological network which was deployed in South Tibet and the Himalayas for a 3-yr period. The large size of the network, the high quality of the seismological data and the new GOCE gravity data set allow us to image the entire lithosphere of this active area in an innovative way. We image 3-D low velocity and density structures in the middle crust that fit the location of discontinuous low S-velocity zones revealed by receiver functions in previous geophysical studies. In the deeper parts of our velocity model we image a positive anomaly interpreted to be the heterogenous Indian lithosphere vertically descending beneath the centre of the Tibetan Plateau.
ISSN: 0956-540X

hal-00858024
https://hal.archives-ouvertes.fr/hal-00858024
DOI : 10.1093/gji/ggs110

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
In this paper, we describe an outcrop to characterize the effect of fracture spacing and type on larger scale effective elasticity, which is measured for the first time in-situ with a Schmidt hammer. The outcrop is dominated by lime mudstones and belongs to the deformation zone of the St Clement fault, in southern France. Our results suggest that small spacing of faults, open fractures and styolites leads to lesser effective Young's modulus, whereas small sealed fracture spacing leads to greater effective Young's modulus. These relationships are compatible with theoretical models of effective elasticity. Using Amadei and Savage (1993) approach, we define a non-linear model that relates Schmidt hammer rebound to spacing by fracture type. A hemisphere with a radius of 40 to similar to 200 cm is the rheological volume characterized by the Schmidt hammer. Results of model inversion demonstrate that variations of Schmidt hammer rebound over the outcrop can be used to estimate fracture type and stiffness. Stiffness of sealed fractures is 2-3 orders of magnitude greater than the stiffness of faults, stylolites and open fractures. This result is consistent with an increase of the rate of interseismic stress build-up of major faults with sealing of fractures in their damage zone.
ISSN: 0191-8141

hal-00669954
https://hal.archives-ouvertes.fr/hal-00669954
DOI : 10.1016/j.jsg.2011.09.009

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Four decades after their introduction, the validity of fault slip inversion methods based on Wallace (1951) and Bott (1959) hypothesis, which states that the slip on each fault surface has the same direction and sense as the maximum resolved shear stress, is still a subject of debate. According to some authors, this hypothesis is questionable since fault mechanical interactions induce slip reorientations as confirmed by geomechanical models. This leads us to ask as to what extent the Wallace-Bott simplifications are reliable as a basis hypothesis for stress inversion from fault slip data. In this paper, we compare two inversion methods; the first is based on the Wallace-Bott hypothesis, and the second relies on geomechanics and mechanical effects on fault heterogeneous slip distribution. In that context, a multi-parametric stress inversion study covering (i) the friction coefficients (μ), (ii) the full range of Andersonian state of stress and (iii) slip data sampling along the faults is performed. For each tested parameter, the results of the mechanical stress inversion and the Wallace-Bott (WB) based stress inversion for slip are compared in order to understand their respective effects. The predicted discrepancy between the solutions of both stress inversion methods (based on WB and mechanics) will then be used to explain the stress inversions results for the chimney Rock case study. It is shown that a high solution discrepancy is not always correlated with the misfit angle (ω ) and can be found under specific configurations (View the MathML sourceR−, θ, μ, geometry) invalidating the WB solutions. We conclude that in most cases the mechanical stress inversion and the WB based stress inversion are both valid and complementary depending on the fault friction. Some exceptions (i.e. low fault friction, simple fault geometry and pure regimes) that may lead to wrong WB based stress inversion solutions are highlighted.
ISSN: 0040-1951

hal-01467165
https://hal.archives-ouvertes.fr/hal-01467165
DOI : 10.1016/j.tecto.2016.11.039

Éditeur(s) : HAL CCSD American Geophysical Union - AGU
Résumé : IODP Expedition 335 "Superfast Spreading Rate Crust 4" returned to ODP Hole 1256D with the intent of deepening this reference penetration of intact ocean crust several hundred meters into cumulate gabbros. This was the fourth cruise of the superfast campaign to understand the formation of oceanic crust accreted at fast spreading ridges, by exploiting the inverse relationship between spreading rate and the depth to low velocity zones seismically imaged at active mid-ocean zones, thought to be magma chambers. Site 1256 is located on 15-million-year-old crust formed at the East Pacific Rise during an episode of superfast ocean spreading (>200 mm/yr full rate). Three earlier cruises to Hole 1256D have drilled through the sediments, lavas and dikes and 100 m into a complex dike-gabbro transition zone. The specific objectives of IODP Expedition 335 were to: (1) test models of magmatic accretion at fast spreading ocean ridges; (2) quantify the vigor of hydrothermal cooling of the lower crust; (3) establish the geological meaning of the seismic Layer 2-3 boundary at Site 1256; and (4) estimate the contribution of lower crustal gabbros to marine magnetic anomalies. It was anticipated that even a shortened IODP Expedition could deepen Hole 1256D a significant distance (300 m) into cumulate gabbros. Operations on IODP Expedition 335 proved challenging from the outset with almost three weeks spent re-opening and securing unstable sections of the Hole. When coring commenced, the destruction of a hard-formation C9 rotary coring bit at the bottom of the hole required further remedial operations to remove junk and huge volumes of accumulated drill cuttings. Hole-cleaning operations using junk baskets returned large samples of a contact-metamorphic aureole between the sheeted dikes and a major heat source below. These large (up to 3.5 kg) irregular samples preserve magmatic, hydrothermal and structural relationships hitherto unseen because of the narrow diameter of drill core and previous poor core recovery. Including the ~60 m-thick zone of granoblastic dikes overlying the uppermost gabbro, the dike-gabbro transition zone at Site 1256 is over 170 m thick, of which more than 100 m are recrystallized granoblastic basalts. This zone records a dynamically evolving thermal boundary layer between the principally hydrothermal domain of the upper crust and a deeper zone of intrusive magmatism. The recovered samples document a sequence of evolving geological conditions and the intimate coupling between temporally and spatially intercalated intrusive, hydrothermal, contact-metamorphic, partial melting and retrogressive processes. Despite the operational challenges, we achieved a minor depth advance to 1522 m, but this was insufficient penetration to complete any of the primary objectives. However, Hole 1256D has been thoroughly cleared of junk and drill cuttings that have hampered operations during this and previous Expeditions. At the end of Expedition 335, we briefly resumed coring and stabilized problematic intervals with cement. Hole 1256D is open to its full depth and ready for further deepening in the near future.
InterRidge News

hal-00688821
https://hal.archives-ouvertes.fr/hal-00688821

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
The interpretation of coseismic surface deformation measurements through inversion techniques is of major importance to understand the mechanical behaviour of a seismic fault. Dense geodetic data sets in the vicinity of the ruptured fault provide unique constraints on detailed fault geometry and slip distribution at depth, making them complementary to seismological data. Bam earthquake (Mw 6.6, 2003 December 26) induced surface deformation has been precisely mapped by Envisat ASAR interferometry and by subpixel correlation techniques applied to Spot-5 and ASAR amplitude images. These oblique and horizontal estimations of deformation have been completed with one levelling profile along the main road crossing the rupture from west to east. We process these data (separately and jointly) in a two-step inversion technique, within the elastic half-space theory framework. Our objective is to determine the dislocation model at depth that satisfies simultaneously all the geodetic constraints. Also, we estimate the relative contribution of each geodetic data set to this inversion process.We first use a stochastic direct approach called neighbourhood algorithm in order to estimate the average characteristics of the rupture, and their relative uncertainty. Constraining in this way the geometry of the ruptured fault, we then linearize the inverse problem and compute the slip distribution on the fault using a standard weighted least-square technique, assuming the solution is smooth to some degree. At each step, we discuss the optimal models, their stability as well as the relative influence of each data set on the derived models parameters. Our preferred model reveals a shallow dislocation on a quasi-vertical fault, slightly dipping towards east. The slip vector has a strike-slip component as high as 2 m, while the dip-slip component seems negligible.However, the estimation of the resolution matrices emphasizes the fact that the details of deep fault slip distribution remain out of the scope of this ill-conditioned inverse problem. Yet, our preferred model suggests a main dislocation limited at depth between 1 and 6 km. By contrast, the aftershocks observed in the months following the earthquake are located just beneath the estimated main shock.
ISSN: 0956-540X

hal-00314173
https://hal.archives-ouvertes.fr/hal-00314173
https://hal.archives-ouvertes.fr/hal-00314173/document
https://hal.archives-ouvertes.fr/hal-00314173/file/GJI-2007-169-849-briole.pdf
DOI : 10.1111/j.1365-246X.2007.03358.x

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Energy dissipation is observed on seismic data when a wave propagates through a porous medium, involving different frequency regimes depending on the nature of rock and fluid types. We focus here on the role of partial fluid saturation in unconsolidated porous media, looking in particular at P-wave phase velocity and attenuation. The study consists in running an experiment in a sand-filled tank partially saturated with water. Seismic propagation in the tank is generated in the kHz range by hitting a steel ball on a granite plate. Seismic data are recorded by buried accelerometers and injecting or extracting water controls the partial saturation. Several imbibition/drainage cycles were performed between the water and gas residual saturations. A Continuous Wavelet Transform applied on seismic records allowed us to extract the direct P-wave at each receiver.We observe an hysteresis in phase velocities and inverse quality factors between imbibition and drainage. Phase velocities and inverse quality factors are then jointly inverted to get a final poro-visco-elastic model of the partially saturated sand that satisfactorily reproduces the data. The model formulation consists in generalizing the Biot theory to effective properties of the fluid and medium (permeability and bulk modulus) in order to properly explain the phase velocity variation as a function of the saturation. The strong level of attenuation measured experimentally is further explained by an anelastic effect due to grain to grain sliding, adding to Biot's losses. This study shows that fluid distribution at microscopic scale has strong influence on the attenuation of direct P-waves at macroscopic scale and confirms that seismic prospection may be a powerful tool for the characterization of transport phenomena in porous media.
ISSN: 0956-540X

hal-00779920
https://hal.archives-ouvertes.fr/hal-00779920
https://hal.archives-ouvertes.fr/hal-00779920/document
https://hal.archives-ouvertes.fr/hal-00779920/file/Barriere_et_al_HAL.pdf
DOI : 10.1111/j.1365-246X.2012.05691.x

Éditeur(s) : HAL CCSD Oxford University Press (OUP): Policy P - Oxford Open Option A Oxford University Press (OUP)
Résumé : International audience
The Pyrenean mountain range is a slowly deforming belt with continuous and moderate seismic activity. To quantify its deformation field, we present the velocity field estimated from a GPS survey of the Pyrenees spanning 18 yr. The PotSis and ResPyr networks, including a total of 85 GPS sites, were installed and first measured in 1992 and 1995–1997, respectively, and remeasured in 2008 and 2010. We obtain a deformation field with velocities less than 1 mm yr−1 across the range. The estimated velocities for individual stations do not differ significantly from zero with 95 per cent confidence. Even so, we estimate a maximum extensional horizontal strain rate of 2.0 ± 1.7 nanostrain per year in a N–S direction in the western part of the range. We do not interpret the vertical displacements due to their large uncertainties. In order to compare the horizontal strain rates with the seismic activity, we analyse a set of 194 focal mechanisms using three methods: (i) the ‘r’ factor relating their P and T axes, (ii) the stress tensors obtained by fault slip inversion and (iii) the strain-rate tensors. Stress and strain-rate tensors are estimated for: (i) the whole data set, (ii) the eastern and western parts of the range separately, and (iii) eight zones, which are defined based on the seismicity and the tectonic patterns of the Pyrenees. Each of these analyses reveals a lateral variation of the deformation style from compression and extension in the east to extension and strike-slip in the west of the range. Although the horizontal components of the strain-rate tensors estimated from the seismic data are slightly smaller in magnitude than those computed from the GPS velocity field, they are consistent within the 2σ uncertainties. Furthermore, the orientations of their principal axes agree with the mapped active faults.
ISSN: 0956-540X

hal-01171752
https://hal.archives-ouvertes.fr/hal-01171752
DOI : 10.1093/gji/ggv052

Éditeur(s) : HAL CCSD
Résumé : The mantle deformation processes that control the thinning and break-up of continental lithosphere remain poorly understood. Our knowledge is restricted to either lithospheric scale thermo-mechanical models --that use experimentally derived flow laws--, geophysical imaging and/or rare xenoliths from active continental rifts, such as the East African Rift System. The originality of this work relies on the study of the two largest outcrops of diamond facies subcontinental lithospheric mantle in the world: the Beni Bousera and Ronda peridotite massifs in N Morocco and S Spain, respectively. The structures and petrologic and metamorphic zoning preserved in these massifs --implying a polybaric and polythermal evolution-- provides a unique opportunity to investigate the thermo-mechanical evolution of thick subcontinental lithospheric mantle in extensional settings. In this thesis we studied the deformation mechanisms in both peridotites and pyroxenites to constrain the modes of exhumation of subcontinental lithospheric mantle from garnet-, to spinel-, and finally, to plagioclase lherzolite facies conditions. We combined field mapping of tectonometamorphic domains and structural mapping of ductile structures, microstructural analysis, crystal preferred orientations (CPO) measurements and conventional thermobarometric calculations and thermodynamic modeling (Perple_X) to unravel the pressure and temperature conditions of deformation. We showed that exhumation from garnet- to spinel lherzolite facies conditions was accommodated by fast shearing --in thermal disequilibrium--along a lithospheric scale transtensional shear zone. In this context, the petrological zoning and the large temperature gradient (ca. 100ºC/km) preserved in the Beni Bousera massif represent the mechanical juxtaposition of progressively deeper and hotter lithospheric levels at depths of ca. 60 km in the latest Oligocene (ca. 25 Ma). Final exhumation from spinel- to plagioclase facies lherzolite and emplacement into the crust is best recorded in the Ronda massif where it occurred by inversion and lithospheric scale folding of the highly attenuated continental lithosphere in a back-arc region, probably in relation with southward slab rollback and subsequent collision with the palaeo-Maghrebien passive margin in the early Miocene (21-23 Ma).
Les processus de déformation contrôlant l'amincissement de la lithosphère continentale sont encore mal contraints. Nos connaissances sont principalement basées sur la modélisation thermomécanique d'extension à l'échelle de la lithosphère--utilisant des lois rhéologiques derivées expérimentalement, l'imagerie géophysique et l'analyse de xénolithes provenant de rift continentaux actifs à ce jour, tels que le Rift Est-Africain. L'originalité de ce travail reside dans l'étude des deux plus grands massifs de péridotites sous-continentales ayant enregistrées des conditions primaires du facies à diamant: les massifs de Beni Bousera au nord du Maroc et de Ronda au sud de l'Espagne, respectivement. Les structures et la zonation petrologique et métamorphique --impliquant une évolution polybarique et polythermique-- préservéees dans ces massifs offrent une opportunité unique pour étudier l'évolution thermo-mécanique du manteau sous-continental dans un contexte extensif. Dans ce travail, nous avons étudié les mécanismes de déformation des péridotites et des pyroxénites afin de contraindre les modes d'exhumation du manteau lithosphérique sous-continental, depuis des conditions du facies des lherzolites à grenat, jusqu'au facies à spinelle et enfin à plagioclase. Nous avons combiné la cartographie des faciès tectono-métamorphiques et des structures ductiles de déformation, l'analyse des microstructures, la mesure d'orientations préférentielles de réseau (OPR), et la géothermobarométrie conventionelle couplée à la modélisation thermodynamique (PerpleX) afin de contraindre les conditions de pression et température de la déformation. Nous avons montré que l'exhumation précoce du facies à grenat au facies à spinelle était accomodée par une faille transtensive affectant le manteau lithosphérique. Dans ce contexte, la zonation tectono-métamorphique et le gradient thermique important (ca. 100ºC/km) préservés à Beni Bousera résultent de la juxtaposition mécanique de domaines lithosphériques initialement équilibrés à différentes pressions et températures, fossilisée à une profondeur de ca. 60 km durant l'Oligocène supérieur (ca. 25 Ma). L'exhumation finale du facies de lherzolite à spinelle au facies à plagioclase et l'emplacement final dans la croûte, mieux enregistrés dans Ronda, se sont produits par inversion et plissement de la section lithosphérique fortement amincie dans un contexte arrière-arc, probablement lors du retrait vers le sud de la lithosphère subduite et la collision de l'arc avec les paléo-marges maghrébines au Miocène inférieur (21-23 Ma).
https://tel.archives-ouvertes.fr/tel-00767700

tel-00767700
https://tel.archives-ouvertes.fr/tel-00767700
https://tel.archives-ouvertes.fr/tel-00767700/document
https://tel.archives-ouvertes.fr/tel-00767700/file/THESE_Frets.pdf

Éditeur(s) : HAL CCSD Hindawi Publishing Corporation
Résumé : During the 1982-1984 bradyseismic crises in the Campi Flegrei area (Italy), the University of Wisconsin deployed a network of seismological stations to record local earthquakes. In order to analyse the potential of the recorded data in terms of tomographic imaging, a blind test was recently set up and carried out in the framework of a research project. A model representing a hypothetical 3D structure of the area containing the Campi Flegrei caldera was also set up, and a synthetic dataset of time arrivals was in turn computed. The synthetic dataset consists of several thousand P- and S-time arrivals, computed at about fourteen stations. The tomographic inversion was performed by four independent teams using different methods. The teams had no knowledge of either the input velocity model or the earthquake hypocenters used to create the synthetic dataset. The results obtained by the different groups were compared and analysed in light of the true model. This work provides a thorough analysis of the earthquake tomography potential of the dataset recording the seismic activity at Campi Flegrei in the 1982-1984 period. It shows that all the tested earthquake tomography methods provide reliable low-resolution images of the background velocity field of the Campi Flegrei area, but with some differences. However, none of them succeeds in detecting the hypothetical structure details (i.e. with a size smaller than about 1.5-2 km), such as a magmatic chamber 4 km deep and especially the smaller, isolated bodies, which represent possible magmatic chimneys and intrusions.
ISSN: 1687-885X

hal-00795584
https://hal.archives-ouvertes.fr/hal-00795584
DOI : 10.1155/2012/505286

Éditeur(s) : HAL CCSD American Geophysical Union (AGU)
Résumé : International audience
This paper presents a new balanced structural cross-section of the Jaca thrust-sheet-top basin of the southern Pyrenees combined with paleo-thermometry and apatite fission track (AFT) thermochronology data. The cross-section, based on field data and interpretation of industrial seismic reflection profiles, allows refinement of previous interpretations of the south-directed thrust system, involving the identification of new thrust faults, and of the kinematic relationships between basement and cover thrusts from the middle Eocene to the early Miocene. AFT analysis shows a southward decrease in the level of fission track resetting, from totally reset Paleozoic rocks and lower Eocene turbidites (indicative of heating to Tmax > ~120 °C), to partially reset middle Eocene turbidites and no/very weak resetting in the upper Eocene-lower Oligocene molasse (Tmax < ~60 °C). AFT results indicate a late Oligocene-early Miocene cooling event throughout the Axial Zone and Jaca Basin. Paleo-maximum temperatures determined by vitrinite reflectance measurements and Raman spectroscopy of carbonaceous material reach up to ~240 °C at the base of the turbidite succession. Inverse modelling of AFT and vitrinite reflectance data with the QTQt software for key samples show compatibility between vitrinite-derived Tmax and the AFT reset level for most of the samples. However, they also suggest that the highest temperatures determined in the lowermost turbidites correspond to a thermal anomaly rather than burial heating, possibly due to fluid circulation during thrust activity. From these results, we propose a new sequential restoration of the south Pyrenean thrust system propagation and related basin evolution.
ISSN: 0278-7407

insu-01312746
https://hal-insu.archives-ouvertes.fr/insu-01312746
https://hal-insu.archives-ouvertes.fr/insu-01312746/document
https://hal-insu.archives-ouvertes.fr/insu-01312746/file/Labaume_et_al-2016-Tectonics.pdf
DOI : 10.1002/2016TC004192

Éditeur(s) : HAL CCSD Wiley
Résumé : International audience
In this paper we investigate finite-frequency effects in crustal tomography. We developed an inversion procedure based on an exact numerical computation of the sensitivity kernels. In this approach we compute the 3D travel-time sensitivity kernels by using (1) graph theory and an additional bending to estimate accurately both rays and travel-times between source/receiver and diffraction points and (2) paraxial ray theory to estimate the amplitude along theses rays. We invert both the velocity and the hypocentre parameters, using these so-called banana-doughnut kernels and the LSQR iterative solver. We compare the ray-theoretical and the finite-frequency tomography to image the intermediate structures beneath the Gulf of Corinth (Greece), which has long been recognized as the most active continental rifting zone in the Mediterranean region. Our dataset consists of 451 local events with 9233 P- firstarrival times recorded in the western part of the Gulf (Aigion area) in the framework of the 3F-Corinth European project. Previous tomographic images showed a complex velocity crustal model and a low-dip surface that may accommodate the deformation. Accurate velocity models will help to better constrain the rifting process, which is still a subject of debate. The main results of this study show that finite-frequency tomography improves crustal tomographic images by providing better resolved images of the 3D complicated velocity structure. Because the kernels spread the information over a volume, finite-frequency tomography results in a sharpening of layer boundaries as we observed for the shallower part of the crust (down to 5 km depth) beneath the Gulf of Corinth.
ISSN: 0016-8025

insu-00354702
https://hal-insu.archives-ouvertes.fr/insu-00354702
DOI : 10.1111/j.1365-2478.2007.00683.x

Éditeur(s) : HAL CCSD
Résumé : When the upper crust below the mountain ranges shows crustal extension with a direction collinear to the convergence direction, it is generally accepted that it is the result of gravitational collapse. However, recent studies challenge this paradigm by showing that erosion induces uplift and extension in the central part of the low convergent mountain ranges. In this study, we investigate the impact of denudation on the seismotectonic regime of mountain ranges.In order to identify a relationship between seismotectonic regime and erosion, the first part of this work presents a compilation of data in the mountain ranges. Based on these observations, a simple kinematic model is proposed to predict the seismotectonic regime of the study ranges.For low convergence rate and moderate mean elevation ranges, this model predicts an extension regime when the denudation rate is 15% higher than the convergence rate.The second part is devoted to the development of thermomechanical 2D finite element models to study the impact of surface processes on the deformation of the Pyrenees. The results show that the isostatic response to erosion reactivates pre-existing structures. The kinematics of an inherited fault plane can be predicted due to the gradient of the horizontal surface velocity profile. Thus, a plane located in the eroded zone is reactivated in normal fault when near the border of the erosion area, this same plane is reactivated in reverse fault. These results suggest that the current deformation in the North Pyrenean Zone could be the result of surface processes.Given the small number of studies quantifying erosion rates in the Pyrenees, the models developed in the second part suffer from high uncertainty. To remedy this, we study the central Pyrenees through a study that combines two types of data : watershed denudation rates from cosmogenic isotopes concentration, and incision rates from sediments buried in the karst. These results are presented in chapter 3. Denudation profiles obtained are consistent with reactivation of a normal fault plane located in the North Zone of the Central Pyrenees.From our study in the Pyrenees, a good correlation appears between the value of the rate of erosion and the glacier occupation. Given that the glaciers where more developed in the Alps, it raises the question of the impact of the late-würmian deglaciation on the present deformation of this range. A numerical model detailing this relationship is presented in the fourth chapter. Theresults show that deglaciation of the western Alps is controlled by the rheological heterogeneity of the crust. Some of our models predict uplift rates consistent with those highlighted by geodesy.
Lorsque l’extension de la croûte sous les parties hautes des chaînes de montagnes est colinéaire à la direction de convergence, il est traditionnellement admis que le moteur est l’effondrement gravitaire. Pourtant, des études récentes remettent en cause ce paradigme en montrant que l’érosion induit un soulèvement et de l’extension dans la partie centrale des chaînes de montagne à faible taux de convergence. L’objectif de notre étude est d’étudier l’impact de la dénudation de la topographie sur le régime sismo-tectonique des chaînes de montagnes.La première partie de ce travail présente une compilation de données dans les chaînes de montagnes afin de dégager des relations entre régime sismo-tectonique et érosion. Sur la base de ces observations, un modèle cinématique simple permettant de prévoir le régime de la chaîne est proposé. Ainsi, pour les chaînes à faible taux de convergence et d’élévation moyenne, ce modèleprédit de l’extension lorsque le taux de dénudation est 15% plus élevé que le taux de convergence.La deuxième partie est consacrée au développement d’un modèle thermo-mécanique 2D en éléments finis pour étudier l’impact des processus de surface sur la déformation des Pyrénées.Les résultats montrent que la réponse isostatique à l’érosion permet de réactiver des structurespré-existantes. La cinématique d’un plan de faille hérité peut être prédite grâce au gradient du profil des vitesses de surface horizontales. Ainsi, un plan situé dans la zone d’érosion est réactivé en faille normale alors qu’en bordure de cette zone une faille est réactivée en régime inverse. Ces résultats suggèrent que la déformation actuelle des Pyrénées pourrait être la conséquence d’unprocessus d’érosion.Compte tenu du faible nombre d’études quantifiant les taux d’érosion dans les Pyrénées, les modèles développés dans la deuxième partie souffrent d’une forte incertitude. Pour y remédier, nous avons cherché à les quantifier dans les Pyrénées Centrales grâce à une étude qui combine deux types de données : taux de dénudation des bassins versants à partir des isotopes cosmogéniques, et vitesses d’incision à partir des sédiments piégés dans les karsts. Ces résultats sontprésentés dans la troisième partie. Les profils de dénudation obtenus sont compatibles avec unrejeu en faille normale d’un plan situé dans la Zone Nord Pyrénéenne des Pyrénées Centrales.Notre étude dans les Pyrénées, montre une bonne corrélation entre la valeur du taux d’érosion et la cartographie de l’occupation glaciaire. Les glaciers ayant été plus développés dans les Alpes, cela pose la question de l’impact de la déglaciation tardi-Wurmienne sur la déformation actuelle de cette chaîne. Un modèle numérique détaillant cette relation est présenté dans le quatrième chapitre. Les résultats montrent que la déglaciation des Alpes occidentales est contrôléepar l’hétérogénéité rhéologique de la croute. Certains de nos modèles prédisent des vitesses de surrection compatibles avec celles mises en évidence par la géodésie.
https://hal.archives-ouvertes.fr/tel-01262333

tel-01262333
https://hal.archives-ouvertes.fr/tel-01262333
https://hal.archives-ouvertes.fr/tel-01262333/document
https://hal.archives-ouvertes.fr/tel-01262333/file/these-Genti-2015-postjury.pdf

Éditeur(s) : HAL CCSD Geological Society
Résumé : International audience
Deformation mechanisms, long-term kinematics and evolution of fold and thrust belts subjected to erosion are studied through 2D analogue experiments involving large convergence. First-order parameters tested include (1) décollements and/or plastic layers interbedded at different locations within analogue materials and (2) synconvergence surface erosion. Weak layers, depending on their location in the model, favour deformation partitioning characterized by the simultaneous development of underplating domains in the inner part of the wedge (basal accretion) and frontal accretion where the wedge grows forward. Interaction between tectonics and surface processes influences this behaviour. Development of antiformal thrust stacks controlled by underplating shows small- and large-scale cyclicity. Thin plastic layers induce folding processes, which are studied at wedge scale. Recumbent and overturned folds, with large inverted limbs, develop in a shear-induced asymmetric deformation regime via progressive unrolling of synclinal hinges. Surface erosion and underplating at depth induce further rotation (passive tilting) and horizontalization of fold limbs. Model results give insights to discuss the mechanisms responsible for the large-scale structures (i.e. antiformal nappe stacks, klippen and kilometre-scale recumbent fold-nappes) encountered in several mountain belts such as the Montagne Noire (French Massif Central), the Galicia Variscan belt (Spain) and the northern Apennines (Italy).
ISSN: 0016-7649

hal-00937644
https://hal.archives-ouvertes.fr/hal-00937644
DOI : 10.1144/jgs2013-012

Éditeur(s) : HAL CCSD
Résumé : This thesis has been conducted on a fractured crystalline basement located in Ploemeur (South of Brittany, France). The detailed study of this basement requires a good understanding of their internal structures, namely their intrinsic properties, flow mode, density and fractures distribution. This work is focused on the petrophysical characterization of the granitic formations. Physical properties were determined from downhole geophysical measurements and images, and petrophysical measurements with thin sections made for each mini-core in the laboratory. Fresh-altered schist and fresh-altered granite separated by a normal angle fault are the main lithology of the Ploemeur crystalline basement. The results showed that the petrophysical properties were function of both porosity variation and degree of alteration and that the pore-space topology was found to be porositydependent. This thesis introduces a new formula between electrical resistivity and porosity. In addition, the laboratory measurements of Ploemeur revealed that seismic wave velocity was a function not only of porosity, but also of the predominant pore type. In general, the pore aspect ratio inverted from velocity measurements were classified roughly into crack shape. In the Ploemeur crystalline basement, the downhole geophysical measurements and images were used to measure changes in physical and chemical properties in the borehole fluid and the bedrock surrounding the boreholes. The effective integration of different geophysical methods, hydrogeological tests and lithological description was the best way to characterize the fractures in the crystalline basement, and is critical for a suitable estimation of fluid circulation.
Le socle cristallin fracturé, situé à Ploemeur (Bretagne sud, France) est l'objet de ce travail de thèse. L'étude détaillée du socle nécessite une bonne compréhension de la structure interne, à savoir les propriétés intrinsèques, le mode d'écoulement, la densité et la distribution des fractures. Ce travail de thèse est consacré à la caractérisation pétrophysique des formations granitiques. Les propriétés physiques ont été déterminées à partir de l'interprétation des données de mesures de géophysique en forage, des images de paroi du puits et des mesures pétrophysiques réalisées sur carottes en laboratoire. Afin de replacer les mesures dans leur contexte lithologique, une lame mince a été analysée pour chaque mini-carotte. Les formations lithologiques rencontrées dans le socle cristallin de Ploemeur sont des schistes et granites sains et altérés, ces deux lithologies étant séparées par une zone de faille. Les résultats ont montré que les propriétés pétrophysiques des roches étaient fonction à la fois de la variation de la porosité et du degré d'altération des roches et que la topologie de l'espace poreux était fortement dépendante de la porosité. Cette thèse introduit une nouvelle formule entre la résistivité électrique et la porosité. De plus, les mesures en laboratoire sur les échantillons de Ploemeur révèlent que les vitesses sismiques ne sont pas seulement fonction de la quantité de porosité, mais aussi fonction de facteur de forme. Ce dernier a été calculé sur les échantillons du socle de Ploemeur à partir des vitesses acoustiques. Les résultats ont montré que la porosité est de forme fissurale. L'intégration des données issues des images optiques et acoustiques, des mesures de géophysique en forage, et de la description des carottes a permis d'étudier de façon détaillée la géologie et l'hydrodynamique du socle de Ploemeur. L'intégration effective de ces différentes méthodes avec les tests hydrogéologiques ressort comme la meilleure façon de caractériser les fractures et les circulations en eau dans le socle.
https://tel.archives-ouvertes.fr/tel-00444175

Droits : info:eu-repo/semantics/OpenAccess
tel-00444175
https://tel.archives-ouvertes.fr/tel-00444175
https://tel.archives-ouvertes.fr/tel-00444175/document
https://tel.archives-ouvertes.fr/tel-00444175/file/thesebelghoul.pdf

Éditeur(s) : HAL CCSD Copernicus Publications
Résumé : International audience
During Ocean Drilling Program Leg 134, a mixed sequence of volcaniclastic and carbonate sediments was cored together with basaltic breccias and a single basalt sill in Hole 833B. In this study, data from Schlumberger's Geochemical Logging Tool, together with physical property logs, characterise part of the section and establish a complete log-derived lithostratigraphy. The major transition between the sedimentary rocks and a basaltic sill at 830 mbsf (metres below sea floor) is clearly seen as well as the distinction between the silicate-rich sediments and the carbonate-rich part of the section. The transition between the volcanic sandstone interbedded with basaltic and sed-lithic breccia of Unit III, and the volcanic sandstone of Unit IV is clearly marked by the physical property logs. The accuracy of the logging data was evaluated and discussed with an improved log versus core correlation. The cluster analysis of the complete data set shows a correlation with the data provided by the mechanical caliper device, which in turn indicates the quality of the borehole.
ISSN: 1816-8957

hal-00642165
http://hal.upmc.fr/hal-00642165
http://hal.upmc.fr/hal-00642165/document
http://hal.upmc.fr/hal-00642165/file/Rabaute_et_al-SciDrilling-1997.pdf