Ce travail présente quelques contributions à la modélisation<br />et à la simulation de genèse de bassins sédimentaires.<br />Nous présentons tout d'abord les modèles mathématiques et <br />les schémas numériques mis en oeuvre à l'Institut Français<br />du Pétrole dans le cadre du projet Temis. Cette première partie <br />est illustrée à l'aide de tests numériques portant sur des bassins 1D/2D.<br />Nous étudions ensuite le schéma amont des pétroliers utilisé pour la résolution des équations de Darcy et nous établissons des résultats mathématiques nouveaux<br />dans le cas d'un écoulement de type Dead-Oil.<br />Nous montrons également comment construire un schéma à nombre<br />de Péclet variable en présence de pression capillaire. <br />Là encore, nous effectuons une étude mathématique<br />détaillée et nous montrons la convergence du schéma<br />dans un cas simplifié. Des tests numériques réalisés<br />sur un problème modèle montrent que l'utilisation d'un nombre<br />de Péclet variable améliore la précision des calculs.<br />Enfin nous considérons dans une dernière partie <br />un modèle d'écoulement où les changements de lithologie et <br />les changements de courbes de pression capillaire sont liés.<br />Nous précisons la condition physique que doivent vérifier<br />les solutions en saturation aux interfaces de changement de roche et<br />nous en déduisons une formulation faible originale.<br />L'existence d'une solution à ce problème est obtenue<br />par convergence d'un schéma volumes finis.<br />Des exemples numériques montrent l'influence de la condition<br />d'interface sur le passage ou la retenue des hydrocarbures.

La Mer d'Alboran est un bassin arrière-arc situé à l'extrémité occidentale de la Mer Méditerranée ayant subi une évolution complexe. Basé sur une approche " terre-mer ", pluridisciplinaire, ce travail axé sur deux grandes thématiques de recherche vise à mieux comprendre les modalités d'ouverture et de déformation des bassins sédimentaires au cours du Néogène en étudiant notamment un événement majeur de l'histoire de la Mer Méditerranée, la Crise de salinité messinienne.L'étude tectonique et sédimentaire du bassin Ouest-Alboran, a permis de construire un nouveau schéma d'évolution tectonique et sédimentaire au cours du Miocène. L'initiation de la subsidence s'est effectuée au gré de grandes zones de cisaillement crustales affectant le socle métamorphique. La subsidence est ensuite contrôlée par la traction du panneau plongeant lithosphérique sous-jacent. Dans les Cordillères Bétiques, une seconde génération de bassins, dont celui de Sorbas, s'est développée le long de dômes métamorphiques exhumés au cours d'une phase d'extension. Des reconstitutions paléogéographiques montrent que la formation des bassins de la Mer d'Alboran est en relation avec la subduction sous-jacente.Des études de terrain ont mis en évidence une importante phase d'érosion des bassins périphériques de la Mer d'Alboran en réponse à la Crise de salinité messinienne. Cette érosion subaérienne eut pour conséquence le creusement de canyons fluviatiles à terre qui ont été suivis en mer par l'analyse fine de profils sismiques. Un nouveau scénario est alors proposé afin d'expliquer la réouverture de la connexion avec l'Océan Atlantique après plus de 160.000 ans d'émersion. / The Alboran Sea is a back-arc basin located at the western end of the Mediterranean Sea, and has been affected by complex tectonic settings during its history. Based on a multidisciplinary "land-sea" approach, this work focuses on two major research topics: it aims (1) to better understand the opening and deformation processes of two key sedimentary basins during the Neogene, and (2) to study the effects and imprints of the Messinian Salinity Crisis onto the basins.Tectonic and stratigraphic analysis of the thickest sedimentary accumulation called Western Alboran Basin (WAB) led to build an original tectonic scenario in the Miocene. Inception of extension was triggered by major crustal shear zones favoring the exhumation of the metamorphic basement. Then, the subsidence became most probably controlled by the slab-pull effect of the underlying oceanic slab. In the Betic Cordilleras, a second generation of sedimentary basins, as the Sorbas Basin, developed along metamorphic core complexes exhumed during an extension regime. A set of paleogeographic maps proposes to explain the formation of the Alboran basins in relation with deep geodynamic processes.Field studies evidenced that peripheral basins bounding the Alboran Sea have been affected by a significant subaerial erosional phase in response to the Messinian Salinity Crisis. Stratigraphic analysis of seismic profiles demonstrated the onshore/offshore continuity of fluvial canyon morphologies recognized on land. A new scenario is then proposed to explain the reopening of the worldwide oceanic connection with the Mediterranean realm after more than 160.000 years of emersion.

Les méthodes de restauration en coupe ou en carte sont des outils performants pour estimer les déplacements dans les domaines déformés. De plus, elles permettent de tester la cohésion géométrique d'une carte ou d'une coupe, laquelle est considérée comme acceptable d'un point de vu cinématique si son état restauré est géologiquement acceptable. Etant donné que les méthodes de restauration en coupe sont limitées aux domaines de déformation plane, nous avons mis au point une méthode numérique de restauration en carte conçue pour les domaines faillés en extension. Nous utilisons des cartes de niveaux stratigraphiques sur lesquelles sont reportés le champ de failles et la composante horizontale du rejet normal sur chacune de ces failles (carte de rejets horizontaux). Nous considérons que le niveau est restauré lorsque le rejet sur les failles a été inversé. Pour ce faire, nous "découpons" la carte de rejets horizontaux en un ensemble de blocs limités par des failles. Les rejets sur les failles sont inversés en ajustant les blocs les uns par rapport aux autres par une série de rotations et de translations rigides visant à minimiser par moindre carré les vides et les recouvrements entre les blocs. A partir de l'état déformé et de l'état restauré de l'horizon stratigraphique, nous établissons (1) un champ des déplacements finis (2) un champ de rotations finies des blocs et (3) un champ de déformation. Nous avons appliqué cette méthode à trois exemples: deux bassins sédimentaires (le bassin de Campos sur la marge atlantique brésilienne et le delta du Niger sur la marge atlantique nigériane) et un domaine d'extension à l'échelle de la croûte (la dépression Afar à la jonction de la Mer Rouge, du Golfe D'Aden et du Rift Est Africain). Nous avons choisi d'appliquer notre méthode dans le bassin de Campos parce que le glissement gravitaire divergent de la couverture sédimentaire a produit un étirement affectant plusieurs directions horizontales c'est à dire une déformation non plane. Le jeu de données dans le delta du Niger nous permet de calculer des champs de déplacements et de déformation incrémentaux en restaurant plusieurs niveaux stratigraphiques au sein d'un même bassin. Dans les deux cas, les données sont issues d'une couverture sismique. La dépression Afar est un exemple qui nous permet d'illustrer notre méthode lorsqu'elle est appliquée à une échelle différente en utilisant des cartes issues de l'interprétation d'image satellites SPOT. Dans le bassin de Campos nous avons restauré deux marqueurs stratigraphiques: un marqueur oligocène et un marqueur crétacé supérieur. Dans les deux cas nous obtenons un champ de déplacement avec des vecteurs divergents vers le large. Cependant le champ de déplacement et le champ de failles sont différents d'un horizon à l'autre. Nous suggérons que ce soit le résultat de la mise en place au Tertiaire d'un delta avec un front de propagation triangulaire. La surcharge sédimentaire associée serait le moteur du l'étalement gravitaire divergent et aurait induit un nouveau champ de failles. Dans le delta du Niger, nous avons appliqué notre méthode à six marqueurs stratigraphiques au sein d'un système de failles de croissance correspondant à deux failles majeures antithétiques et à un "roll-over" faillé. Nous avons calculé les champs de déplacements correspondant aux incréments de temps entre les dépôts des horizons successifs. Nous montrons que les déplacements horizontaux les plus importants ont lieu au cours des dernières étapes du développement du roll-over, au moment de son effondrement gravitaire. Deux modèles ont été proposés jusqu'à maintenant pour expliquer la structuration du Sud-Est de la dépression Afar. Soit l'extension est considérée comme le résultat du recouvrement entre deux rifts' se propageant dans des directions opposées et est accommodée par un mécanisme de rotations de blocs (bookshelt). Soit elle est produite par la rotation anti-horaire de la bordure Est de la dépression. Nous avons déterminé les champs de déplacements, de rotations et de déformations dans ce secteur. Nous montrons qu'il existe un gradient SW-NE dans l'intensité des rotations de blocs, des déformations et des déplacements. De plus il existe un changement des directions de déformation et de délpacement depuis une direction NNE-SSW dans la partie ouest de la dépression vers une direction NE-SW à l'Est. Ces résultats sont compatibles avec une rotation antihoraire de la bordure Est mais pas avec un mécanisme de bookshelf qui devrait être associé à une distribution uniforme des intensités de rotations et de déformation.

Ce mémoire constitue la synthèse des travaux de recherche effectués dans les Andes de la région de Cusco (Pérou). Cette région s'étend entre la bordure NE de la Cordillère Occidentale et la bordure SW de la Cordillère Orientale et constitue la terminaison NW de l'Altiplano. Les modèles interprétatifs issus de ces études sont importants pour la compréhension de l'évolution géodynamique des Andes au Méso-Cénozoïque. Ainsi une large part est faite aux résultats de la confrontation des méthodes structurales, sédimentologiques et magmatiques. On montre que l'évolution paléozoïque et méso-cénozoïque de la zone étudiée serait en partie contrôlée par une limite paléogéographique avec deux types différents de substratum. Cette limite se situe au niveau de la déflexion d' Abancay et du seuil Cusco-Puno. Le Permo-Trias doit être considéré comme une période transitoire entre le cycle hercynien et le cycle andin. Les principaux traits du bassin Permo-Trias sont des structures héritées (y compris la déflexion d'Abancay), qui ont aussi contrôlé toute l'évolution andine. La région de Cusco est caractérisée par plus de 10.000m de couches rouges tertiaires d'origine continentale qui surmontent des sédiments marins et continentaux d'âge Crétacé à Paléocène. Ces Couches Rouges tertiaires se sont déposées dans des bassins synorogéniques liés aux déformations compressives andines. Des mouvements décrochants dextres apparaissent vers 50-44 Ma le long du seuil Cusco-Puno et en bordure NE du bassin mésozoïque sud-péruvien, conduisant à l'ouverture en transtension des bassins de Couches Rouges du Groupe San Jeronimo de Cusco, Sicuani, Ayaviri, Puno et vraisemblablement du Nord de la Bolivie. A cette époque le grand bassin altiplanique aurait commencé à se structurer. Un important changement du régime de contraintes vers 43-40 Ma provoque la fermeture en compression des bassins de Couches Rouges et le début de la sédimentation des conglomérats de la Formation Anta (événement compressif Inca1). Ce changement induit postérieurement une évolution du comportement de la limite Cordillère Occidentale-Bassin Couches Rouges, qui devient le siège de mouvements décrochants senestres, comme le montrent l'évolution des bassins Oligocènes-Miocènes. Ces mouvements senestres sont compatibles avec des modèles qui tentent d'expliquer la formation de l'Orocline bolivien. Les contraintes compressives (crises Quechua) sont des marqueurs de la fermeture des bassins néogènes. Les données géochronologiques semblent indiquer que l' activité magmatique dans la zone étudiée a été plus ou moins continue depuis l'Eocène (48 Ma) et que la répartition du magmatisme montre, par contre, une nette évolution dans l'espace. Les calculs de la déformation indiquent que les valeurs de raccourcissement obtenues à partir des données géologiques de surface, ne suffisent pas à produire l'épaississement crustal mis en évidence par la géophysique. Pour expliquer cet excès de volume on propose un sous-charriage de matériel continental provenant de la subduction ou l'apport de copeaux litbosphériques résultant de l'érosion tectonique. Cependant, si le doublement de la croute est entièrement dû au raccourcissement,la déformation en surface ne serait bien exprimée qu' en quelques endroits. La cinématique de plaques semble bien avoir un rôle sur la tectonique andine, ses crises et périodes ou événements, et la déformation est contrôlée par les grands accidents paléogéographiques.

Tout mécanisme de formation d'un bassin sédimentaire porte une signature thermique caractéristique. Inversement, la reconstitution de l'histoire thermique d'un bassin sédimentaire permet d'en tirer des conclusions sur les mécanismes et l'histoire de sa formation. Les méthodes géothermométriques nous permettent d'étudier l'histoire thermique d'un bassin. Les enregistrements réalisés nous renseignent, soit sur les conditions de température (ou conditions T/P) instantanées (températures de fermeture, de cristallisation, etc.), soit sur les effets du couple température-temps (énergie) au cours de la diagenèse (transformation de la matière organique, transformation des minéraux argileux). Ce mémoire présente les résultats d'une étude réalisée dans l'avant-pays alpin de Suisse et de Haute-Savoie (France). Le bassin molassique s'étend le long du front externe de l'arc alpin sur plus de 800 km, de Chambéry à Vienne. Depuis le front alpin la sédimentation de la Molasse, essentiellement détritique, pro grade vers le nord-ouest. Cet élargissement du bassin au cours du temps est matérialisé par la position des biseaux appartenant à des groupes lithostratigraphiques successifs. L'épaisseur des dépôts varie de quelques dizaines de mètres dans la partie distale jurassienne à plus de 4000 m dans la partie proximale subalpine. Sur le plan lithostratigraphique, la Molasse peut être divisée en quatre groupes où les faciès prédominants sont alternativement marins et terrestres: UMM (Molasse marine inférieure, Rupélien-Chattien), USM (Molasse d'eau douce inférieure, Rupélien-Burdigalien?), OMM (Molasse marine supérieure, Burdigalien-Langhien?) et OSM (Molasse d'eau douce supérieure, LanghienSerravalien). Du nord-ouest au sud-est, le bassin molassique est divisé en trois ensembles tectoniques: i) la Molasse du Jura, ii) la Molasse du Plateau et iii) la Molasse subalpine. Par le biais d'une meilleure compréhension des processus qui contrôlent le régime paléogéothermique, l'objectif de cette étude est la reconstitution de l' histoire thermique de quelques périmètres du bassin molassique au cours du Tertiaire. La démarche suivie peut être décomposée en trois phases principales: - La phase d'échantillonnage. Afin de rendre les résultats obtenus les plus significatifs possibles, les échantillons ont été prélevés avec une large dispersion géographique. Du sud-ouest au nord-est, cinq régions ont été explorées: Haute-Savoie (Bornes, Plateau de Bornes, Bassin de Bellegarde, Bassin de Rumilly) et bassin genevois, Suisse occidentale (région entre le lac de Neuchâtel et le lac de Genève), région du lac de Thoune, Suisse orientale (région centrée autour du lac de Zürich), Suisse septentrionale (région située dans la prolongation des failles du fossé rhénan vers le sud, dans le Jura tabulaire et le Jura plissé); - La phase d'analyse. Quatre types de méthodes ont été utilisées: étude de la réflectance de la vitrinite, étude diffractométrique des argiles, étude des inclusions fluides et pyrolyse RockEva!. - La cartographie des résultats et la modélisation de l'histoire thermique.

Afin de cibler une meilleure rentabilité des gisements d'hydrocarbures, la simulation numérique présente un intérêt grandissant dans le secteur pétrolier. Dans ce contexte, deux principaux types d'applications sont distingués : l'ingénierie du réservoir, où les modèles géologiques sont définis comme statiques et l'exploration des gisements par la simulation de la formation des hydrocarbures. Cette dernière application nécessite des modèles de bassins dynamiques. Pour ce type de simulation, la modélisation mathématique et numérique a connu une avancée importante durant les dernières années. En revanche, la construction de maillages indispensables pour les simulations reste une tâche lourde dans le cas de géométries complexes et dynamiques. La difficulté se traduit par la particularité du domaine à mailler, qui est dictée par la mécanique des milieux granulaires, c'est-à-dire des milieux quasi-incompressibles et hétérogènes. De plus, les données sismiques sont des surfaces non-implicites modélisant des blocs volumiques. Dans ce cadre, nous nous intéressons à la génération de maillages lagrangiens hexa-dominants des bassins à géométrie complexe. Le maillage souhaité doit contenir une couche de mailles principalement hexaédriques entre deux horizons et respecter les surfaces failles traversant ces horizons. Une méthodologie originale basée sur la construction d’un espace déplié est proposée. Le maillage souhaité est alors traduit par une grille 3D contrainte dans cet espace. Plusieurs techniques d’optimisation de maillages sont aussi proposées / To target more profitable oil and gas deposits, the numerical simulation is of growing interest in the oil sector. In this context, two main types of applications are distinguished: reservoir engineering where geological models are defined as static and exploration for the simulation of hydrocarbons formation. The latter application requires dynamic basins models. For this type of simulation, mathematical and numerical modeling has been an important advance in recent years. However, the construction of meshes needed for the simulations is a difficult task in the case of complex and dynamic geometries. The difficulty is reflected by the characteristic of the domain to mesh, which is defined from the mechanics of granular media which are almost incompressible and heterogeneous environments. In addition, the seismic data represent non-implicit surfaces modeling volume blocks. In this context, we focus on Lagrangian hex-dominant mesh generation of basins with complex geometry. The desired mesh must contain a layer of almost hexahedral meshes between two horizons and conform to fault surfaces through these horizons. A novel methodology based on the construction of an unfolded space is introduced. The desired mesh is then seen as a constrained 3D grid in this novel space. Several mesh optimization techniques are also proposed

Une étude de caractérisation intégrée de pétrographie et géochimie organique a été réalisée dans les bassins de Sabinas et Chihuahua au Nord-Est du Mexique. Ces informations ont permis une modélisation numérique de la formation du gaz en considérant la subsidence thermique des charbons et des shales à gaz. Les objectifs de cette thèse sont les suivants: mise en place d'une méthodologie de caractérisation des roches étudiées ; estimation du potentiel générateur de gaz et sa distribution régionale ; estimation des ressources en gaz méthane. Les analyses effectuées sur l'ensemble des échantillons ont permis de caractériser la roche, les kérogènes, les huiles et les gaz. La méthodologie utilisée permet de croiser les informations pétrographiques et géochimiques pour pouvoir analyser le système pétrolier par modélisation numérique. Les analyses réalisées sont: réflectance% Ro, analyse élémentaire et immédiate, Rock Eval6 ® (Bulk Rock), [delta]13C, [delta]D, (gaz de charbon), microscopie électronique à balayage, analyse d'images, analyse des macéraux et inclusions fluides. Un programme informatique a été construit afin de croiser l'information des analyses des échantillons avec celle des expériences de maturation artificielle en laboratoire. Cette démarche nous à permis d'estimer les ressources en gaz méthane à la fois généré par les charbons et les shales gaz. La méthodologie utilisée pour l'étude des kérogènes et des gaz dans les deux bassins est à notre avis la plus appropriée pour répondre aux objectifs proposés dans ce travail de thèse. Celle-ci permet de caractériser l'évolution thermique des charbons, shales gas (gaz de schiste), et roches mères pétrolières / This comprehensive characterization study was performed using organic petrology and geochemistry conducted in the Sabinas basin and Chihuahua in northern Mexico. This information allowed a numerical modeling of gas formation, considering the thermal subsidence of coal and carbonaceous shales.The objectives of this thesis are: establish a characterization methodology for the studied rocks ; estimate potential gas generator and its regional distribution ; estimate the methane gas resources. For the development of this project, we conducted an intensive campaign representative sampling of coal, carbonaceous shales and coal gas "in situ". For the Sabinas basin were studied 97 samples and 114 samples in the basin of Chihuahua. The analyses carried out that were used on the samples analyzed allowed to characterize the kerogen and gas. The methodology used to cross petrographic and geochemical information to analyze the petroleum system by numerical modeling. Analyses were: Petrographic, reflectance %Ro, elemental analysis and immediate, Rock Eval6 ® (Bulk rock), isotopic analysis, [delta]13C, [delta]D, (coal gas), scanning electron microscopy, image analysis and analysis of macerals fluid inclusions. The analyzes that were used on the samples allowed to characterize the sample, the kerogen and gas. The methodology used to cross petrographic and geochemical information for analyze the oil system by numerical modeling. Analyses were: Petrographic, reflectance% Ro, elemental analysis and immediate, Rock Eval6 ® (Bulk rock), isotopic analysis, d13C, dD, (coal gas), scanning electron microscopy, image analysis and analysis of macerals fluid inclusions. A computer program was constructed to cross the information with the analysis of samples of artificial maturation experiments in the laboratory. This approach allowed estimation of methane gas resources generated by coal and carbonaceous shales. [...] In conclusion, between the two basins studied, the Sabinas basin, generates more gas for industrial exploitation. However, the rocks of the basin of Chihuahua unconventional types (shale gas) show good potential for hydrocarbon generation. The methodology used for the study of kerogen and gas in the two basins is in our opinion the most appropriate to meet the objectives proposed in this thesis. This allows to characterize the thermal evolution of coal, carbonaceous shale (shale gas), coal gas and hydrocarbon source rocks

Le soufre est présent dans les systèmes pétroliers à la fois sous forme organique et minérale. Il est impliqué dans de nombreuses réactions d'oxydoréduction qui affectent la qualité des huiles (par des réactions de sulfuration ou de désulfuration) et du gaz naturel (par la génération de H2S en contexte de réduction thermochimique des sulfates), ainsi que la porosité des roches réservoirs (par la dissolution de l'anhydrite ou la précipitation de soufre élémentaire ou de pyrobitume). Ces réactions sont gouvernées par la température (et dans une moindre mesure la pression), les conditions d'oxydoréduction et la composition chimique globale du système. La thèse comporte trois parties : (1) une étude thermodynamique des réactions chimiques impliquant le soufre dans les bassins sédimentaires ; (2) une étude expérimentale des réactions de sulfuration et de désulfuration des composés organiques ainsi que de la réduction thermochimique des sulfates; (3) une étude pétrographique et d'inclusions fluides sur des échantillons d'un réservoir carbonaté du bassin Pré-Caspien / Sulfur occurs in petroleum systems as both organic compounds and minerals as well as under different oxidation states. It is involved in a number of redox reactions which may impact the quality of crude oils (through sulfurization or desulfurization reactions) and natural gas (through the generation of H2S during thermochemical sulfate reduction), as well as the petrophysical properties of reservoir rocks (through the dissolution of anhydrite and the precipitation of elemental sulfur and pyrobitumen). These reactions are controlled by temperature (and to a lesser extent pressure), the redox conditions, and the overall chemical composition of the system representing the petroleum reservoir. The thesis consists of three parts: (1) a thermodynamic study of chemical reactions involving sulfur which occur in sedimentary basins; (2) an experimental simulation of sulfurization and desulfurization reactions of organic compounds, as well as of thermochemical sulfate reduction; and (3) a petrographic and fluid inclusion study of carbonate rock samples from a sulfur-rich hydrocarbon reservoir of the northern Caspian Sea

L'évolution tectonique de l'avant-pays de la chaîne alpine de Tunisie comporte une étude paléotectonique et néotectonique qui utilise la Méthode de l'Inversion des données sur les systèmes de failles (Carey & Brunier 1974; Carey 1976) afin de définir les champs de contraintes régionaux à différentes époques. L'étude microtectonique des failles contemporaines de la sédimentation (failles synsédimentaires) a permis de mettre en évidence la succession des états de paléo-contraintes au cours du Mésozoïque et du Cénozoïque. Ainsi, ont été reconnu avec une bonne cohérence : -Au Permo-carbonifère, une extension NE-SW (failles normales E-W et NS et les grabens NW-SE) -Au Trias, une extension NE-SW (Horsts, grabens et demi-grabens NW-SE) -Au Jurassique, une extension N-S (Failles normales proches de E-W) -Au Néocomien-Barrémien, extension NE-SW (Horsts et des grabens NW-SE) -Au cours de l'Albo-Aptien, une compression proche de E-W (paléo-anticlinaux N20E) -Au Vraconien-Cénomanien inférieur, une extension NE-SW (failles normales NE-SW) -Au Cénomanien supérieur, une compression locale proche de E-W (Failles inverses subméridiens) -Au Coniacien-Santonien, une compression locale proche de N-S (décrochements conjugués) -Au Campanien, une distension NE-SW (failles normales NW-SE) -Au Paléocène supérieur-Eocène, une compression NW-SE (paléo-plis et failles inverses NE-SW) -A l'Oligocène, une extension proche de NE-SW (failles normales et grabens NW-SE ) -Au Langhien-Serravallien, une extension ENE-WSW (failles normales et grabens N120 à N140 E). Ces épisodes paléo-tectoniques sont associés à des montées diapiriques et/ou volcaniques. Au cours du plissement tortonien (" phase atlasique "), les failles normales anciennes proches de E-W et N-S sont réactivées en décrochements E-W dextres et N-S sénestres. Ce système de décrochements engendrent des plis en échelon et découpe l'Atlas tunisien selon les réseaux de décrochements de Riedel de Types D, R, R'. Un essai de zonation tectonique néogène a permis de dégager les variations spatiales du régime tectonique à différentes époques du Langhien-Serravallien, Tortonien et Messinien-Pliocène. L'étude néotectonique et sismotectonique a permis de mettre en évidence une succession des états de contraintes du Pléistocène à l'Actuel, avec : -Au Villafranchien (Pléistocène inférieur) une phase de plissement (NW-SE à NNW-SSE) donne des plis dans les bassins molassiques, le Sahel et sur les bordures des décrochements majeurs E-W et N-S. -Au Pléistocène supérieur, la compression NW-SE à N-S se poursuit avec des failles inverses dans le Nord et un système de décrochements E-W dextres et N-S sénestres dans l'Atlas et le Sahel. A partir des données de la sismicité et de l'archéosismicité, une zonation sismotectonique a été proposée. Elle permet de distinguer deux (2) grandes provinces de déformations actuelles : une province à déformation compressive à décrochante-compressive, localisée dans la zone de collision du Nord du pays et une province à déformation distensive à décrochante distensive, située dans le domaine de l'intra-plaque de l'avant-pays (Atlas central et méridional).