É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 Elsevier
Résumé : International audience
Wallace (1951) and Bott (1959) were the first to introduce the idea that the slip on each fault surface has the same direction and sense as the maximum shear stress resolved on that surface. However, this simplified hypothesis is questionable since fault mechanical interactions may induce slip reorientations. Earlier numerical geomechanical models confirmed that the slickenlines (slip vectors) are not necessarily parallel to the maximum resolved shear stress but are consistent with local stress perturbations. This leads us to ask as to what extent the Wallace and Bott simplifications are reliable as a basis hypothesis for stress inversion from fault slip data. Here, a geomechanical multi-parametric study using a 3D boundary element method, covering (i) fault geometries such as intersected faults or corrugated fault surfaces, (ii) the full range of Andersonian state of stress, (iii) fault friction, (iv) fault fluid pressure, (v) half space effect and (vi), rock properties, is performed in order to understand the effect of each parameter on the misfit angle between geomechanical slip vectors and the resolved shear stresses. It is shown that significant misfit angles can be found under specific configurations invalidating the Wallace and Bott assumptions, even though fault friction tends to minimize the misfit. We therefore conclude that in such cases, stress inversions based on fault slip data should be interpreted with care.
ISSN: 0040-1951

hal-01218418
https://hal.archives-ouvertes.fr/hal-01218418
DOI : 10.1016/j.tecto.2015.06.030

Éditeur(s) : HAL CCSD Springer Verlag
Résumé : Some major challenges for geophysicists and structural geologists using three-dimensional boundary element method codes (3D-BEM) are: (1) reducing the amount of memory required to solve large and dense systems and (2) incorporation of inequality constraints such as traction inequality constraints (TIC) and displacement inequality constraints (DIC). The latter serves two purposes. First, for example, inequality constraints can be used to simulate frictional slip (using TIC). Second, these constraints can prevent element interpenetration while allowing opening mode (using DIC). We have developed a method that simultaneously incorporates both types of functionality of the inequality constraints. We show that the use of an appropriate iterative solver not only avoids the allocation of significant memory for solving the system (allowing very large model computation and simplifying parallelization on multi-core processors), but also admits interesting features such as natural incorporation of TICs and DICs. Compared to other techniques of contact management (e.g., Lagrange multipliers, penalty method, or complementarity problem), this new simple methodology, which does not use any incremental trial-and-error procedures, brings more flexibility, while making the system more stable and less subject to round-off errors without any computational overhead. We provide validations and comparisons of the inequality constraints implementation using 2D analytical and numerical solutions.
ISSN: 1420-0597

hal-00561035
https://hal.archives-ouvertes.fr/hal-00561035
DOI : 10.1007/s10596-009-9170-x

É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